Search

BACKGROUND: Chronic exposure to microorganisms inside homes can impact respiratory health. Few studies have used advanced sequencing methods to examine adult respiratory outcomes, especially continuous measures. We aimed to identify metagenomic profiles in house dust related to the quantitative traits of pulmonary function and airway inflammation in adults. Microbial communities, 1264 species (389 genera), in vacuumed bedroom dust from 779 homes in a US cohort were characterized by whole metagenome shotgun sequencing. We examined two overall microbial diversity measures: richness (the number of individual microbial species) and Shannon index (reflecting both richness and relative abundance). To identify specific differentially abundant genera, we applied the Lasso estimator with high-dimensional inference methods, a novel framework for analyzing microbiome data in relation to continuous traits after accounting for all taxa examined together. RESULTS: Pulmonary function measures (forced expiratory volume in one second (FEV1), forced vital capacity (FVC), and FEV1/FVC ratio) were not associated with overall dust microbial diversity. However, many individual microbial genera were differentially abundant (p-value < 0.05 controlling for all other microbial taxa examined) in relation to FEV1, FVC, or FEV1/FVC. Similarly, fractional exhaled nitric oxide (FeNO), a marker of airway inflammation, was unrelated to overall microbial diversity but associated with differential abundance for many individual genera. Several genera, including Limosilactobacillus, were associated with a pulmonary function measure and FeNO, while others, including Moraxella to FEV1/FVC and Stenotrophomonas to FeNO, were associated with a single trait. CONCLUSIONS: Using state-of-the-art metagenomic sequencing, we identified specific microorganisms in indoor dust related to pulmonary function and airway inflammation. Some were previously associated with respiratory conditions; others were novel, suggesting specific environmental microbial components contribute to various respiratory outcomes. The methods used are applicable to studying microbiome in relation to other continuous outcomes. Video Abstract.

Studies using 16S rRNA and shotgun metagenomics typically yield different results, usually attributed to PCR amplification biases. We introduce Greengenes2, a reference tree that unifies genomic and 16S rRNA databases in a consistent, integrated resource. By inserting sequences into a whole-genome phylogeny, we show that 16S rRNA and shotgun metagenomic data generated from the same samples agree in principal coordinates space, taxonomy and phenotype effect size when analyzed with the same tree.

Phylogenetic trees provide a framework for organizing evolutionary histories across the tree of life and aid downstream comparative analyses such as metagenomic identification. Methods that rely on single-marker genes such as 16S rRNA have produced trees of limited accuracy with hundreds of thousands of organisms, whereas methods that use genome-wide data are not scalable to large numbers of genomes. We introduce updating trees using divide-and-conquer (uDance), a method that enables updatable genome-wide inference using a divide-and-conquer strategy that refines different parts of the tree independently and can build off of existing trees, with high accuracy and scalability. With uDance, we infer a species tree of roughly 200,000 genomes using 387 marker genes, totaling 42.5 billion amino acid residues.

In 2020, we identified cancer-specific microbial signals in The Cancer Genome Atlas (TCGA) [1]. Multiple peer-reviewed papers independently verified or extended our findings [2-12]. Given this impact, we carefully considered concerns by Gihawi et al. [13] that batch correction and database contamination with host sequences artificially created the appearance of cancer type-specific microbiomes. (1) We tested batch correction by comparing raw and Voom-SNM-corrected data per-batch, finding predictive equivalence and significantly similar features. We found consistent results with a modern microbiome-specific method (ConQuR [14]), and when restricting to taxa found in an independent, highly-decontaminated cohort. (2) Using Conterminator [15], we found low levels of human contamination in our original databases (~1% of genomes). We demonstrated that the increased detection of human reads in Gihawi et al. [13] was due to using a newer human genome reference. (3) We developed Exhaustive, a method twice as sensitive as Conterminator, to clean RefSeq. We comprehensively host-deplete TCGA with many human (pan)genome references. We repeated all analyses with this and the Gihawi et al. [13] pipeline, and found cancer type-specific microbiomes. These extensive re-analyses and updated methods validate our original conclusion that cancer type-specific microbial signatures exist in TCGA, and show they are robust to methodology.

Multiomics has shown promise in noninvasive risk profiling and early detection of various common diseases. In the present study, in a prospective population-based cohort with ~18 years of e-health record follow-up, we investigated the incremental and combined value of genomic and gut metagenomic risk assessment compared with conventional risk factors for predicting incident coronary artery disease (CAD), type 2 diabetes (T2D), Alzheimer disease and prostate cancer. We found that polygenic risk scores (PRSs) improved prediction over conventional risk factors for all diseases. Gut microbiome scores improved predictive capacity over baseline age for CAD, T2D and prostate cancer. Integrated risk models of PRSs, gut microbiome scores and conventional risk factors achieved the highest predictive performance for all diseases studied compared with models based on conventional risk factors alone. The present study demonstrates that integrated PRSs and gut metagenomic risk models improve the predictive value over conventional risk factors for common chronic diseases.

Microbial breakdown of organic matter is one of the most important processes on Earth, yet the controls of decomposition are poorly understood. Here we track 36 terrestrial human cadavers in three locations and show that a phylogenetically distinct, interdomain microbial network assembles during decomposition despite selection effects of location, climate and season. We generated a metagenome-assembled genome library from cadaver-associated soils and integrated it with metabolomics data to identify links between taxonomy and function. This universal network of microbial decomposers is characterized by cross-feeding to metabolize labile decomposition products. The key bacterial and fungal decomposers are rare across non-decomposition environments and appear unique to the breakdown of terrestrial decaying flesh, including humans, swine, mice and cattle, with insects as likely important vectors for dispersal. The observed lockstep of microbial interactions further underlies a robust microbial forensic tool with the potential to aid predictions of the time since death.

Colorectal cancer (CRC) is driven by genomic alterations in concert with dietary influences, with the gut microbiome implicated as an effector in disease development and progression. While meta-analyses have provided mechanistic insight into patients with CRC, study heterogeneity has limited causal associations. Using multi-omics studies on genetically controlled cohorts of mice, we identify diet as the major driver of microbial and metabolomic differences, with reductions in α diversity and widespread changes in cecal metabolites seen in high-fat diet (HFD)-fed mice. In addition, non-classic amino acid conjugation of the bile acid cholic acid (AA-CA) increased with HFD. We show that AA-CAs impact intestinal stem cell growth and demonstrate that Ileibacterium valens and Ruminococcus gnavus are able to synthesize these AA-CAs. This multi-omics dataset implicates diet-induced shifts in the microbiome and the metabolome in disease progression and has potential utility in future diagnostic and therapeutic developments.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by heterogeneous cognitive, behavioral and communication impairments. Disruption of the gut-brain axis (GBA) has been implicated in ASD although with limited reproducibility across studies. In this study, we developed a Bayesian differential ranking algorithm to identify ASD-associated molecular and taxa profiles across 10 cross-sectional microbiome datasets and 15 other datasets, including dietary patterns, metabolomics, cytokine profiles and human brain gene expression profiles. We found a functional architecture along the GBA that correlates with heterogeneity of ASD phenotypes, and it is characterized by ASD-associated amino acid, carbohydrate and lipid profiles predominantly encoded by microbial species in the genera Prevotella, Bifidobacterium, Desulfovibrio and Bacteroides and correlates with brain gene expression changes, restrictive dietary patterns and pro-inflammatory cytokine profiles. The functional architecture revealed in age-matched and sex-matched cohorts is not present in sibling-matched cohorts. We also show a strong association between temporal changes in microbiome composition and ASD phenotypes. In summary, we propose a framework to leverage multi-omic datasets from well-defined cohorts and investigate how the GBA influences ASD.

BackgroundThe gut-lung axis is generally recognized, but there are few large studies of the gut microbiome and incident respiratory disease in adults.ObjectiveWe sought to investigate the association and predictive capacity of the gut microbiome for incident asthma and chronic obstructive pulmonary disease (COPD).MethodsShallow metagenomic sequencing was performed for stool samples from a prospective, population-based cohort (FINRISK02; N = 7115 adults) with linked national administrative health register-derived classifications for incident asthma and COPD up to 15 years after baseline. Generalized linear models and Cox regressions were used to assess associations of microbial taxa and diversity with disease occurrence. Predictive models were constructed using machine learning with extreme gradient boosting. Models considered taxa abundances individually and in combination with other risk factors, including sex, age, body mass index, and smoking status.ResultsA total of 695 and 392 statistically significant associations were found between baseline taxonomic groups and incident asthma and COPD, respectively. Gradient boosting decision trees of baseline gut microbiome abundance predicted incident asthma and COPD in the validation data sets with mean area under the curves of 0.608 and 0.780, respectively. Cox analysis showed that the baseline gut microbiome achieved higher predictive performance than individual conventional risk factors, with C-indices of 0.623 for asthma and 0.817 for COPD. The integration of the gut microbiome and conventional risk factors further improved prediction capacities.ConclusionsThe gut microbiome is a significant risk factor for incident asthma and incident COPD and is largely independent of conventional risk factors.

BackgroundUnrelieved chemotherapy-induced nausea (CIN) occurs 48% of patients undergoing chemotherapy and is one of the most debilitating symptoms that patients report.ObjectiveThe aims of this study were to identify subgroups of patients with distinct CIN profiles and determine how these subgroups differed on demographic and clinical characteristics; severity, frequency, and distress of CIN; and the co-occurrence of common gastrointestinal symptoms.MethodsPatients (n = 1343) completed demographic questionnaire and Memorial Symptom Assessment Scale 6 times over 2 cycles of chemotherapy. Latent class analysis was used to identify subgroups of patients with distinct CIN profiles. Differences among these subgroups were evaluated using parametric and nonparametric statistics.ResultsFour distinct CIN profiles were identified: none (40.8%), increasing-decreasing (21.5%), decreasing (8.9%), and high (28.8%). Compared with the none class, patients in the high class were younger, had a lower annual household income, had child care responsibilities, had a lower Karnofsky Performance Status score and a higher Self-administered Comorbidity Questionnaire score, and were more likely to have received chemotherapy on a 14-day cycle and a highly emetogenic chemotherapy regimen. In addition, patients in the high class reported high occurrence rates for dry mouth, feeling bloated, diarrhea, lack of appetite, abdominal cramps, difficulty swallowing, mouth sores, weight loss, and change in the way food tastes.ConclusionsThat 60% of the patients reported moderate to high CIN occurrence rates confirms that this unrelieved symptom is a significant clinical problem.Implications for practiceNurses need to evaluate patients' level of adherence to their antiemetic regimen and make appropriate referrals for physical therapy, psychological services, and dietary counseling.

Human untargeted metabolomics studies annotate only ~10% of molecular features. We introduce reference-data-driven analysis to match metabolomics tandem mass spectrometry (MS/MS) data against metadata-annotated source data as a pseudo-MS/MS reference library. Applying this approach to food source data, we show that it increases MS/MS spectral usage 5.1-fold over conventional structural MS/MS library matches and allows empirical assessment of dietary patterns from untargeted data.

Despite advances in sequencing, lack of standardization makes comparisons across studies challenging and hampers insights into the structure and function of microbial communities across multiple habitats on a planetary scale. Here we present a multi-omics analysis of a diverse set of 880 microbial community samples collected for the Earth Microbiome Project. We include amplicon (16S, 18S, ITS) and shotgun metagenomic sequence data, and untargeted metabolomics data (liquid chromatography-tandem mass spectrometry and gas chromatography mass spectrometry). We used standardized protocols and analytical methods to characterize microbial communities, focusing on relationships and co-occurrences of microbially related metabolites and microbial taxa across environments, thus allowing us to explore diversity at extraordinary scale. In addition to a reference database for metagenomic and metabolomic data, we provide a framework for incorporating additional studies, enabling the expansion of existing knowledge in the form of an evolving community resource. We demonstrate the utility of this database by testing the hypothesis that every microbe and metabolite is everywhere but the environment selects. Our results show that metabolite diversity exhibits turnover and nestedness related to both microbial communities and the environment, whereas the relative abundances of microbially related metabolites vary and co-occur with specific microbial consortia in a habitat-specific manner. We additionally show the power of certain chemistry, in particular terpenoids, in distinguishing Earth's environments (for example, terrestrial plant surfaces and soils, freshwater and marine animal stool), as well as that of certain microbes including Conexibacter woesei (terrestrial soils), Haloquadratum walsbyi (marine deposits) and Pantoea dispersa (terrestrial plant detritus). This Resource provides insight into the taxa and metabolites within microbial communities from diverse habitats across Earth, informing both microbial and chemical ecology, and provides a foundation and methods for multi-omics microbiome studies of hosts and the environment.

Cancer-microbe associations have been explored for centuries, but cancer-associated fungi have rarely been examined. Here, we comprehensively characterize the cancer mycobiome within 17,401 patient tissue, blood, and plasma samples across 35 cancer types in four independent cohorts. We report fungal DNA and cells at low abundances across many major human cancers, with differences in community compositions that differ among cancer types, even when accounting for technical background. Fungal histological staining of tissue microarrays supported intratumoral presence and frequent spatial association with cancer cells and macrophages. Comparing intratumoral fungal communities with matched bacteriomes and immunomes revealed co-occurring bi-domain ecologies, often with permissive, rather than competitive, microenvironments and distinct immune responses. Clinically focused assessments suggested prognostic and diagnostic capacities of the tissue and plasma mycobiomes, even in stage I cancers, and synergistic predictive performance with bacteriomes.

Abstract While the study of the tumor microbiome and its effects on cancer biology has expanded considerably over the last few years, most of this research focused on bacteria and viruses, leaving behind the fungal kingdom. Recently, a few studies have demonstrated that specific fungi may promote tumor progression, stressing the importance of comprehensively studying the tumor mycobiome and its effects. To address this, we have characterized the mycobiome in 1,183 human tumors and their adjacent tissues, originating from eight major solid tumor types. Staining and imaging demonstrated the presence of fungi in both cancer and immune cells, with tumor-type specific distribution patterns. Quantitative PCR of the fungal 5.8s rDNA revealed the presence of fungal DNA in all tumor types. To characterize the tumor mycobiome and address potential contamination during tissue handling and processing, we subjected all samples, as well as 295 negative controls of different types, to sequencing of the ITS2 region that is situated between fungal rRNA genes. We found cancer-type specific mycobial signatures with relatively high similarity between tumors and their adjacent tissues. While the fungal mycobiome had a lower species richness as compared to the bacterial microbiome of the same tumors, fungi showed significant co-occurrences with specific bacteria, suggesting the existence of ecological niches within the tumors. We also found significant correlations with clinical parameters such as patient’s age, tumor stage, progression-free survival, overall survival, and response to immune checkpoint blockade therapy. Characterization of the tumor mycobiome may add a biologically relevant, previously overlooked, component to be considered in the study of cancer, including its effects on tumor initiation, progression, diagnosis, and response to therapy. Citation Format: Lian Narunsky Haziza, Gregory D. Sepich-Poore, Ilana Livyatan, Omer Asraf, Cameron Martino, Deborah Nejman, Nancy Gavert, Jason E. Stajich, Guy Amit, Antonio González, Stephen Wandro, Gili Perry, Ruthie Ariel, Arnon Meltser, Justin P. Shaffer, Qiyun Zhu, Nora Balint-Lahat, Iris Barshack, Maya Dadian, Einav N. Gal-Yam, Sandip P. Pate, Amir Bashan, Austin D. Swafford, Yitzhak Pilpel, Rob Knight, Ravid Straussman. Pan-cancer characterization of the tumor mycobiome and its clinical effects [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3054.

BackgroundUnrelieved chemotherapy-induced nausea (CIN) is a significant problem for patients with breast cancer (BC).ObjectiveIn a sample of patients with BC who were assessed before their second or third cycle of chemotherapy, study purposes were to evaluate for the occurrence, severity, frequency, and distress associated with CIN; evaluate for differences in demographic and clinical characteristics and gastrointestinal (GI) symptom occurrence rates between patients who did and did not report CIN; and determine which demographic, clinical, and symptom characteristics were associated with the occurrence of CIN.MethodsPatients completed demographic and clinical questionnaires and the Memorial Symptom Assessment Scale for nausea and common GI symptom assessments. Univariate analyses evaluated for differences in demographic and clinical characteristics and GI symptom occurrence between patients who did and did not report CIN. Multiple logistic regression analysis evaluated for characteristics associated with CIN.ResultsOf the 532 patients with BC, 47.2% reported CIN occurrence. Characteristics associated with CIN group membership were poorer functional status, receipt of chemotherapy on a 14-day cycle, and higher occurrence rates of 5 GI symptoms (ie, dry mouth, vomiting, constipation, change in the way food tastes, and lack of appetite; all P < .001).ConclusionsUnrelieved CIN is a common symptom in patients with BC. This study is the first to demonstrate that 5 co-occurring GI symptoms were associated with CIN occurrence.Implications for practiceThis study identified new risk factors for CIN occurrence in patients with BC. Clinicians may be able to initiate additional interventions to alleviate CIN.

Increasing data volumes on high-throughput sequencing instruments such as the NovaSeq 6000 leads to long computational bottlenecks for common metagenomics data preprocessing tasks such as adaptor and primer trimming and host removal. Here, we test whether faster recently developed computational tools (Fastp and Minimap2) can replace widely used choices (Atropos and Bowtie2), obtaining dramatic accelerations with additional sensitivity and minimal loss of specificity for these tasks. Furthermore, the taxonomic tables resulting from downstream processing provide biologically comparable results. However, we demonstrate that for taxonomic assignment, Bowtie2's specificity is still required. We suggest that periodic reevaluation of pipeline components, together with improvements to standardized APIs to chain them together, will greatly enhance the efficiency of common bioinformatics tasks while also facilitating incorporation of further optimized steps running on GPUs, FPGAs, or other architectures. We also note that a detailed exploration of available algorithms and pipeline components is an important step that should be taken before optimization of less efficient algorithms on advanced or nonstandard hardware. IMPORTANCE In shotgun metagenomics studies that seek to relate changes in microbial DNA across samples, processing the data on a computer often takes longer than obtaining the data from the sequencing instrument. Recently developed software packages that perform individual steps in the pipeline of data processing in principle offer speed advantages, but in practice they may contain pitfalls that prevent their use, for example, they may make approximations that introduce unacceptable errors in the data. Here, we show that differences in choices of these components can speed up overall data processing by 5-fold or more on the same hardware while maintaining a high degree of correctness, greatly reducing the time taken to interpret results. This is an important step for using the data in clinical settings, where the time taken to obtain the results may be critical for guiding treatment.

We introduce the operational genomic unit (OGU) method, a metagenome analysis strategy that directly exploits sequence alignment hits to individual reference genomes as the minimum unit for assessing the diversity of microbial communities and their relevance to environmental factors. This approach is independent of taxonomic classification, granting the possibility of maximal resolution of community composition, and organizes features into an accurate hierarchy using a phylogenomic tree. The outputs are suitable for contemporary analytical protocols for community ecology, differential abundance, and supervised learning while supporting phylogenetic methods, such as UniFrac and phylofactorization, that are seldom applied to shotgun metagenomics despite being prevalent in 16S rRNA gene amplicon studies. As demonstrated in two real-world case studies, the OGU method produces biologically meaningful patterns from microbiome data sets. Such patterns further remain detectable at very low metagenomic sequencing depths. Compared with taxonomic unit-based analyses implemented in currently adopted metagenomics tools, and the analysis of 16S rRNA gene amplicon sequence variants, this method shows superiority in informing biologically relevant insights, including stronger correlation with body environment and host sex on the Human Microbiome Project data set and more accurate prediction of human age by the gut microbiomes of Finnish individuals included in the FINRISK 2002 cohort. We provide Woltka, a bioinformatics tool to implement this method, with full integration with the QIIME 2 package and the Qiita web platform, to facilitate adoption of the OGU method in future metagenomics studies. IMPORTANCE Shotgun metagenomics is a powerful, yet computationally challenging, technique compared to 16S rRNA gene amplicon sequencing for decoding the composition and structure of microbial communities. Current analyses of metagenomic data are primarily based on taxonomic classification, which is limited in feature resolution. To solve these challenges, we introduce operational genomic units (OGUs), which are the individual reference genomes derived from sequence alignment results, without further assigning them taxonomy. The OGU method advances current read-based metagenomics in two dimensions: (i) providing maximal resolution of community composition and (ii) permitting use of phylogeny-aware tools. Our analysis of real-world data sets shows that it is advantageous over currently adopted metagenomic analysis methods and the finest-grained 16S rRNA analysis methods in predicting biological traits. We thus propose the adoption of OGUs as an effective practice in metagenomic studies.

Integration of sensory and molecular inputs from the environment shapes animal behaviour. A major site of exposure to environmental molecules is the gastrointestinal tract, in which dietary components are chemically transformed by the microbiota1 and gut-derived metabolites are disseminated to all organs, including the brain2. In mice, the gut microbiota impacts behaviour3, modulates neurotransmitter production in the gut and brain4,5, and influences brain development and myelination patterns6,7. The mechanisms that mediate the gut-brain interactions remain poorly defined, although they broadly involve humoral or neuronal connections. We previously reported that the levels of the microbial metabolite 4-ethylphenyl sulfate (4EPS) were increased in a mouse model of atypical neurodevelopment8. Here we identified biosynthetic genes from the gut microbiome that mediate the conversion of dietary tyrosine to 4-ethylphenol (4EP), and bioengineered gut bacteria to selectively produce 4EPS in mice. 4EPS entered the brain and was associated with changes in region-specific activity and functional connectivity. Gene expression signatures revealed altered oligodendrocyte function in the brain, and 4EPS impaired oligodendrocyte maturation in mice and decreased oligodendrocyte-neuron interactions in ex vivo brain cultures. Mice colonized with 4EP-producing bacteria exhibited reduced myelination of neuronal axons. Altered myelination dynamics in the brain have been associated with behavioural outcomes7,9-14. Accordingly, we observed that mice exposed to 4EPS displayed anxiety-like behaviours, and pharmacological treatments that promote oligodendrocyte differentiation prevented the behavioural effects of 4EPS. These findings reveal that a gut-derived molecule influences complex behaviours in mice through effects on oligodendrocyte function and myelin patterning in the brain.

Human genetic variation affects the gut microbiota through a complex combination of environmental and host factors. Here we characterize genetic variations associated with microbial abundances in a single large-scale population-based cohort of 5,959 genotyped individuals with matched gut microbial metagenomes, and dietary and health records (prevalent and follow-up). We identified 567 independent SNP-taxon associations. Variants at the LCT locus associated with Bifidobacterium and other taxa, but they differed according to dairy intake. Furthermore, levels of Faecalicatena lactaris associated with ABO, and suggested preferential utilization of secreted blood antigens as energy source in the gut. Enterococcus faecalis levels associated with variants in the MED13L locus, which has been linked to colorectal cancer. Mendelian randomization analysis indicated a potential causal effect of Morganella on major depressive disorder, consistent with observational incident disease analysis. Overall, we identify and characterize the intricate nature of host-microbiota interactions and their association with disease.

Ulcerative colitis (UC) is driven by disruptions in host-microbiota homoeostasis, but current treatments exclusively target host inflammatory pathways. To understand how host-microbiota interactions become disrupted in UC, we collected and analysed six faecal- or serum-based omic datasets (metaproteomic, metabolomic, metagenomic, metapeptidomic and amplicon sequencing profiles of faecal samples and proteomic profiles of serum samples) from 40 UC patients at a single inflammatory bowel disease centre, as well as various clinical, endoscopic and histologic measures of disease activity. A validation cohort of 210 samples (73 UC, 117 Crohn's disease, 20 healthy controls) was collected and analysed separately and independently. Data integration across both cohorts showed that a subset of the clinically active UC patients had an overabundance of proteases that originated from the bacterium Bacteroides vulgatus. To test whether B. vulgatus proteases contribute to UC disease activity, we first profiled B. vulgatus proteases found in patients and bacterial cultures. Use of a broad-spectrum protease inhibitor improved B. vulgatus-induced barrier dysfunction in vitro, and prevented colitis in B. vulgatus monocolonized, IL10-deficient mice. Furthermore, transplantation of faeces from UC patients with a high abundance of B. vulgatus proteases into germfree mice induced colitis dependent on protease activity. These results, stemming from a multi-omics approach, improve understanding of functional microbiota alterations that drive UC and provide a resource for identifying other pathways that could be inhibited as a strategy to treat this disease.

Background & aimsHyperbaric oxygen therapy (HBOT) is a promising treatment for moderate-to-severe ulcerative colitis. However, our current understanding of the host and microbial response to HBOT remains unclear. This study examined the molecular mechanisms underpinning HBOT using a multi-omic strategy.MethodsPre- and post-intervention mucosal biopsies, tissue, and fecal samples were collected from HBOT phase 2 clinical trials. Biopsies and fecal samples were subjected to shotgun metaproteomics, metabolomics, 16s rRNA sequencing, and metagenomics. Tissue was subjected to bulk RNA sequencing and digital spatial profiling (DSP) for single-cell RNA and protein analysis, and immunohistochemistry was performed. Fecal samples were also used for colonization experiments in IL10-/- germ-free UC mouse models.ResultsProteomics identified negative associations between HBOT response and neutrophil azurophilic granule abundance. DSP identified an HBOT-specific reduction of neutrophil STAT3, which was confirmed by immunohistochemistry. HBOT decreased microbial diversity with a proportional increase in Firmicutes and a secondary bile acid lithocholic acid. A major source of the reduction in diversity was the loss of mucus-adherent taxa, resulting in increased MUC2 levels post-HBOT. Targeted database searching revealed strain-level associations between Akkermansia muciniphila and HBOT response status. Colonization of IL10-/- with stool obtained from HBOT responders resulted in lower colitis activity compared with non-responders, with no differences in STAT3 expression, suggesting complementary but independent host and microbial responses.ConclusionsHBOT reduces host neutrophil STAT3 and azurophilic granule activity in UC patients and changes in microbial composition and metabolism in ways that improve colitis activity. Intestinal microbiota, especially strain level variations in A muciniphila, may contribute to HBOT non-response.

BackgroundCrohn's disease (CD), an inflammatory bowel disease (IBD) subtype, results from pathologic interactions between host cells and its resident gut microbes. CD manifests in both isolated disease locations (ileum or colon) or a combination of locations (ileocolonic). To date, a comprehensive understanding of how isolated CD subtypes influence molecular profiles remains outstanding. To address this, we sought to define CD location signatures by leveraging a large cross-sectional feature set captured from the stool of over 200 IBD patients and healthy controls using metaproteomics, shotgun metagenomics, 16S rRNA sequencing, metabolomic profiling, and host genetics paired with clinical endoscopic assessments.ResultsNeither metagenomic nor host genetics alone distinguished CD location subtypes. In contrast, ileal and colonic CD were distinguished using mass spectrometry-based methods (metabolomics or metaproteomics) or a combined multi-omic feature set. This multi-omic feature set revealed colonic CD was strongly associated with neutrophil-related proteins. Additionally, colonic CD displayed a disease-severity-related association with Bacteroides vulgatus. Colonic CD and ulcerative colitis profiles harbored strikingly similar feature enrichments compared to ileal CD, including neutrophil-related protein enrichments. Compared to colonic CD, ileal CD profiles displayed increased primary and secondary bile acid levels and concomitant shifts in taxa with noted sensitivities such as Faecalibacterium prausnitzii or affinities for bile acid-rich environments, including Gammaproteobacteria and Blautia sp. Having shown robust molecular and microbial distinctions tied to CD locations, we leveraged these profiles to generate location-specific disease severity biomarkers that surpass the performance of Calprotectin.ConclusionsWhen compared using multi-omics features, colonic- and ileal-isolated CD subtypes display striking differences that suggest separate location-specific pathologies. Colonic CD's strong similarity to ulcerative colitis, including neutrophil and Bacteroides vulgatus involvement, is also evidence of a shared pathology for colonic-isolated IBD subtypes, while ileal CD maintains a unique, bile acid-driven profile. More broadly, this study demonstrates the power of multi-omics approaches for IBD biomarker discovery and elucidating the underlying biology. Video Abstract.

BackgroundImproving probiotic engraftment in the human gut requires a thorough understanding of the in vivo adaptive strategies of probiotics in diverse contexts. However, for most probiotic strains, these in vivo genetic processes are still poorly characterized. Here, we investigated the effects of gut selection pressures from human, mice, and zebrafish on the genetic stability of a candidate probiotic Lactiplantibacillus plantarum HNU082 (Lp082) as well as its ecological and evolutionary impacts on the indigenous gut microbiota using shotgun metagenomic sequencing in combination with isolate resequencing methods.ResultsWe combined both metagenomics and isolate whole genome sequencing approaches to systematically study the gut-adaptive evolution of probiotic L. plantarum and the ecological and evolutionary changes of resident gut microbiomes in response to probiotic ingestion in multiple host species. Independent of host model, Lp082 colonized and adapted to the gut by acquiring highly consistent single-nucleotide mutations, which primarily modulated carbohydrate utilization and acid tolerance. We cultivated the probiotic mutants and validated that these gut-adapted mutations were genetically stable for at least 3 months and improved their fitness in vitro. In turn, resident gut microbial strains, especially competing strains with Lp082 (e.g., Bacteroides spp. and Bifidobacterium spp.), actively responded to Lp082 engraftment by accumulating 10-70 times more evolutionary changes than usual. Human gut microbiota exhibited a higher ecological and genetic stability than that of mice.ConclusionsCollectively, our results suggest a highly convergent adaptation strategy of Lp082 across three different host environments. In contrast, the evolutionary changes within the resident gut microbes in response to Lp082 were more divergent and host-specific; however, these changes were not associated with any adverse outcomes. This work lays a theoretical foundation for leveraging animal models for ex vivo engineering of probiotics to improve engraftment outcomes in humans. Video abstract.

The number of publicly available microbiome samples is continually growing. As data set size increases, bottlenecks arise in standard analytical pipelines. Faith's phylogenetic diversity (Faith's PD) is a highly utilized phylogenetic alpha diversity metric that has thus far failed to effectively scale to trees with millions of vertices. Stacked Faith's phylogenetic diversity (SFPhD) enables calculation of this widely adopted diversity metric at a much larger scale by implementing a computationally efficient algorithm. The algorithm reduces the amount of computational resources required, resulting in more accessible software with a reduced carbon footprint, as compared to previous approaches. The new algorithm produces identical results to the previous method. We further demonstrate that the phylogenetic aspect of Faith's PD provides increased power in detecting diversity differences between younger and older populations in the FINRISK study's metagenomic data.

Graves' Disease is the most common organ-specific autoimmune disease and has been linked in small pilot studies to taxonomic markers within the gut microbiome. Important limitations of this work include small sample sizes and low-resolution taxonomic markers. Accordingly, we studied 162 gut microbiomes of mild and severe Graves' disease (GD) patients and healthy controls. Taxonomic and functional analyses based on metagenome-assembled genomes (MAGs) and MAG-annotated genes, together with predicted metabolic functions and metabolite profiles, revealed a well-defined network of MAGs, genes and clinical indexes separating healthy from GD subjects. A supervised classification model identified a combination of biomarkers including microbial species, MAGs, genes and SNPs, with predictive power superior to models from any single biomarker type (AUC = 0.98). Global, cross-disease multi-cohort analysis of gut microbiomes revealed high specificity of these GD biomarkers, notably discriminating against Parkinson's Disease, and suggesting that non-invasive stool-based diagnostics will be useful for these diseases.

Carbon fixation by chemoautotrophic microbes such as homoacetogens had a major impact on the transition from the inorganic to the organic world. Recent reports have shown the presence of genes for key enzymes associated with the Wood-Ljungdahl pathway (WLP) in the phylum Actinobacteria, which adds to the diversity of potential autotrophs. Here, we compiled 42 actinobacterial metagenome-assembled genomes (MAGs) from new and existing metagenomic datasets and propose three novel classes, Ca. Aquicultoria, Ca. Geothermincolia and Ca. Humimicrobiia. Most members of these classes contain genes coding for acetogenesis through the WLP, as well as a variety of hydrogenases (NiFe groups 1a and 3b-3d; FeFe group C; NiFe group 4-related hydrogenases). We show that the three classes acquired the hydrogenases independently, yet the carbon monoxide dehydrogenase/acetyl-CoA synthase complex (CODH/ACS) was apparently present in their last common ancestor and was inherited vertically. Furthermore, the Actinobacteria likely donated genes for CODH/ACS to multiple lineages within Nitrospirae, Deltaproteobacteria (Desulfobacterota), and Thermodesulfobacteria through multiple horizontal gene transfer events. Finally, we show the apparent growth of Ca. Geothermincolia and H2-dependent acetate production in hot spring enrichment cultures with or without the methanogenesis inhibitor 2-bromoethanesulfonate, which is consistent with the proposed homoacetogenic metabolism.

Accurate microbial identification and abundance estimation are crucial for metagenomics analysis. Various methods for classification of metagenomic data and estimation of taxonomic profiles, broadly referred to as metagenomic profilers, have been developed. Nevertheless, benchmarking of metagenomic profilers remains challenging because some tools are designed to report relative sequence abundance while others report relative taxonomic abundance. Here we show how misleading conclusions can be drawn by neglecting this distinction between relative abundance types when benchmarking metagenomic profilers. Moreover, we show compelling evidence that interchanging sequence abundance and taxonomic abundance will influence both per-sample summary statistics and cross-sample comparisons. We suggest that the microbiome research community pay attention to potentially misleading biological conclusions arising from this issue when benchmarking metagenomic profilers, by carefully considering the type of abundance data that were analyzed and interpreted and clearly stating the strategy used for metagenomic profiling.

Standard workflows for analyzing microbiomes often include the creation and curation of phylogenetic trees. Here we present EMPress, an interactive web tool for visualizing trees in the context of microbiome, metabolome, and other community data scalable to trees with well over 500,000 nodes. EMPress provides novel functionality-including ordination integration and animations-alongside many standard tree visualization features and thus simplifies exploratory analyses of many forms of 'omic data.IMPORTANCE Phylogenetic trees are integral data structures for the analysis of microbial communities. Recent work has also shown the utility of trees constructed from certain metabolomic data sets, further highlighting their importance in microbiome research. The ever-growing scale of modern microbiome surveys has led to numerous challenges in visualizing these data. In this paper we used five diverse data sets to showcase the versatility and scalability of EMPress, an interactive web visualization tool. EMPress addresses the growing need for exploratory analysis tools that can accommodate large, complex multi-omic data sets.

Understanding the proactive roles of receivers in peer feedback processes is crucial because proactive recipience carries great potential in enhancing the effectiveness of feedback and supporting self-regulated (SRL) and co-regulated learning (CoRL). However, receiver's proactivity has been insufficiently explored and the field lacks a clear understanding of how peer feedback receivers could aid academic self-regulation and co-regulation. This study unpacks the black box through examining different receiver roles in peer feedback dialogue and receiver-triggered SRL and CoRL behaviours in an undergraduate writing course for first-year English majors in China. Data were collected through audio-taped peer feedback dialogue, stimulated recall interviews and journals. Findings revealed a variety of increasingly active receiver roles: respondent, verifier, explicator, negotiator, seeker and generator. Assuming these roles, receivers not only regulated their own learning by self-monitoring works, evaluating the quality of received comments and co-producing feedback but also improved feedback providers' writing and evaluative skills. The study challenges the stereotypical image of passive receivers and argues that receiver proactivity could turn peer feedback into a mutually beneficial learning activity for receivers and providers. Implications for developing receiver proactivity are discussed.

Untargeted mass spectrometry is employed to detect small molecules in complex biospecimens, generating data that are difficult to interpret. We developed Qemistree, a data exploration strategy based on the hierarchical organization of molecular fingerprints predicted from fragmentation spectra. Qemistree allows mass spectrometry data to be represented in the context of sample metadata and chemical ontologies. By expressing molecular relationships as a tree, we can apply ecological tools that are designed to analyze and visualize the relatedness of DNA sequences to metabolomics data. Here we demonstrate the use of tree-guided data exploration tools to compare metabolomics samples across different experimental conditions such as chromatographic shifts. Additionally, we leverage a tree representation to visualize chemical diversity in a heterogeneous collection of samples. The Qemistree software pipeline is freely available to the microbiome and metabolomics communities in the form of a QIIME2 plugin, and a global natural products social molecular networking workflow.

Microbial natural products are important for the understanding of microbial interactions, chemical defense and communication, and have also served as an inspirational source for numerous pharmaceutical drugs. Tropical marine cyanobacteria have been highlighted as a great source of new natural products, however, few reports have appeared wherein a multi-omics approach has been used to study their natural products potential (i.e., reports are often focused on an individual natural product and its biosynthesis). This study focuses on describing the natural product genetic potential as well as the expressed natural product molecules in benthic tropical cyanobacteria. We collected from several sites around the world and sequenced the genomes of 24 tropical filamentous marine cyanobacteria. The informatics program antiSMASH was used to annotate the major classes of gene clusters. BiG-SCAPE phylum-wide analysis revealed the most promising strains for natural product discovery among these cyanobacteria. LCMS/MS-based metabolomics highlighted the most abundant molecules and molecular classes among 10 of these marine cyanobacterial samples. We observed that despite many genes encoding for peptidic natural products, peptides were not as abundant as lipids and lipopeptides in the chemical extracts. Our results highlight a number of highly interesting biosynthetic gene clusters for genome mining among these cyanobacterial samples.

BackgroundThe discovery of membrane-enclosed, metabolically functional organelles in Bacteria has transformed our understanding of the subcellular complexity of prokaryotic cells. Biomineralization of magnetic nanoparticles within magnetosomes by magnetotactic bacteria (MTB) is a fascinating example of prokaryotic organelles. Magnetosomes, as nano-sized magnetic sensors in MTB, facilitate cell navigation along the local geomagnetic field, a behaviour referred to as magnetotaxis or microbial magnetoreception. Recent discovery of novel MTB outside the traditionally recognized taxonomic lineages suggests that MTB diversity across the domain Bacteria are considerably underestimated, which limits understanding of the taxonomic distribution and evolutionary origin of magnetosome organelle biogenesis.ResultsHere, we perform the most comprehensive metagenomic analysis available of MTB communities and reconstruct metagenome-assembled MTB genomes from diverse ecosystems. Discovery of MTB in acidic peatland soils suggests widespread MTB occurrence in waterlogged soils in addition to subaqueous sediments and water bodies. A total of 168 MTB draft genomes have been reconstructed, which represent nearly a 3-fold increase over the number currently available and more than double the known MTB species at the genome level. Phylogenomic analysis reveals that these genomes belong to 13 Bacterial phyla, six of which were previously not known to include MTB. These findings indicate a much wider taxonomic distribution of magnetosome organelle biogenesis across the domain Bacteria than previously thought. Comparative genome analysis reveals a vast diversity of magnetosome gene clusters involved in magnetosomal biogenesis in terms of gene content and synteny residing in distinct taxonomic lineages. Phylogenetic analyses of core magnetosome proteins in this largest available and taxonomically diverse dataset support an unexpectedly early evolutionary origin of magnetosome biomineralization, likely ancestral to the origin of the domain Bacteria.ConclusionsThese findings expand the taxonomic and phylogenetic diversity of MTB across the domain Bacteria and shed new light on the origin and evolution of microbial magnetoreception. Potential biogenesis of the magnetosome organelle in the close descendants of the last bacterial common ancestor has important implications for our understanding of the evolutionary history of bacterial cellular complexity and emphasizes the biological significance of the magnetosome organelle. Video Abstract.

SummaryThe software pipeline SHOGUN profiles known taxonomic and gene abundances of short-read shotgun metagenomics sequencing data. The pipeline is scalable, modular and flexible. Data analysis and transformation steps can be run individually or together in an automated workflow. Users can easily create new reference databases and can select one of three DNA alignment tools, ranging from ultra-fast low-RAM k-mer-based database search to fully exhaustive gapped DNA alignment, to best fit their analysis needs and computational resources. The pipeline includes an implementation of a published method for taxonomy assignment disambiguation with empirical Bayesian redistribution. The software is installable via the conda resource management framework, has plugins for the QIIME2 and QIITA packages and produces both taxonomy and gene abundance profile tables with a single command, thus promoting convenient and reproducible metagenomics research.Availability and implementationhttps://github.com/knights-lab/SHOGUN.

QIIME 2 is a completely re-engineered microbiome bioinformatics platform based on the popular QIIME platform, which it has replaced. QIIME 2 facilitates comprehensive and fully reproducible microbiome data science, improving accessibility to diverse users by adding multiple user interfaces. QIIME 2 can be combined with Qiita, an open-source web-based platform, to re-use available data for meta-analysis. The following basic protocol describes how to install QIIME 2 on a single computer and analyze microbiome sequence data, from processing of raw DNA sequence reads through generating publishable interactive figures. These interactive figures allow readers of a study to interact with data with the same ease as its authors, advancing microbiome science transparency and reproducibility. We also show how plug-ins developed by the community to add analysis capabilities can be installed and used with QIIME 2, enhancing various aspects of microbiome analyses-e.g., improving taxonomic classification accuracy. Finally, we illustrate how users can perform meta-analyses combining different datasets using readily available public data through Qiita. In this tutorial, we analyze a subset of the Early Childhood Antibiotics and the Microbiome (ECAM) study, which tracked the microbiome composition and development of 43 infants in the United States from birth to 2 years of age, identifying microbiome associations with antibiotic exposure, delivery mode, and diet. For more information about QIIME 2, see https://qiime2.org. To troubleshoot or ask questions about QIIME 2 and microbiome analysis, join the active community at https://forum.qiime2.org. © 2020 The Authors. Basic Protocol: Using QIIME 2 with microbiome data Support Protocol: Further microbiome analyses.

Microbial genomes are available at an ever-increasing pace, as cultivation and sequencing become cheaper and obtaining metagenome-assembled genomes (MAGs) becomes more effective. Phylogenetic placement methods to contextualize hundreds of thousands of genomes must thus be efficiently scalable and sensitive from closely related strains to divergent phyla. We present PhyloPhlAn 3.0, an accurate, rapid, and easy-to-use method for large-scale microbial genome characterization and phylogenetic analysis at multiple levels of resolution. PhyloPhlAn 3.0 can assign genomes from isolate sequencing or MAGs to species-level genome bins built from >230,000 publically available sequences. For individual clades of interest, it reconstructs strain-level phylogenies from among the closest species using clade-specific maximally informative markers. At the other extreme of resolution, it scales to large phylogenies comprising >17,000 microbial species. Examples including Staphylococcus aureus isolates, gut metagenomes, and meta-analyses demonstrate the ability of PhyloPhlAn 3.0 to support genomic and metagenomic analyses.

Systematic characterization of the cancer microbiome provides the opportunity to develop techniques that exploit non-human, microorganism-derived molecules in the diagnosis of a major human disease. Following recent demonstrations that some types of cancer show substantial microbial contributions1-10, we re-examined whole-genome and whole-transcriptome sequencing studies in The Cancer Genome Atlas11 (TCGA) of 33 types of cancer from treatment-naive patients (a total of 18,116 samples) for microbial reads, and found unique microbial signatures in tissue and blood within and between most major types of cancer. These TCGA blood signatures remained predictive when applied to patients with stage Ia-IIc cancer and cancers lacking any genomic alterations currently measured on two commercial-grade cell-free tumour DNA platforms, despite the use of very stringent decontamination analyses that discarded up to 92.3% of total sequence data. In addition, we could discriminate among samples from healthy, cancer-free individuals (n = 69) and those from patients with multiple types of cancer (prostate, lung, and melanoma; 100 samples in total) solely using plasma-derived, cell-free microbial nucleic acids. This potential microbiome-based oncology diagnostic tool warrants further exploration.

Urbanization represents a profound shift in human behaviour, and has considerable cultural and health-associated consequences1,2. Here, we investigate chemical and microbial characteristics of houses and their human occupants across an urbanization gradient in the Amazon rainforest, from a remote Peruvian Amerindian village to the Brazilian city of Manaus. Urbanization was found to be associated with reduced microbial outdoor exposure, increased contact with housing materials, antimicrobials and cleaning products, and increased exposure to chemical diversity. The degree of urbanization correlated with changes in the composition of house bacterial and microeukaryotic communities, increased house and skin fungal diversity, and an increase in the relative abundance of human skin-associated fungi and bacteria in houses. Overall, our results indicate that urbanization has large-scale effects on chemical and microbial exposures and on the human microbiota.