Robustness of cancer microbiome signals over a broad range of methodological variation

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.

Signatures of HIV and Major Depressive Disorder in the Plasma Microbiome

Inter-individual differences in the gut microbiome are linked to alterations in inflammation and blood-brain barrier permeability, which may increase the risk of depression in people with HIV (PWH). The microbiome profile of blood, which is considered by many to be typically sterile, remains largely unexplored. We aimed to characterize the blood plasma microbiome composition and assess its association with major depressive disorder (MDD) in PWH and people without HIV (PWoH). In this cross-sectional, observational cohort, we used shallow-shotgun metagenomic sequencing to characterize the plasma microbiome of 151 participants (84 PWH and 67 PWoH), all of whom underwent a comprehensive neuropsychiatric assessment. The microbial composition did not differ between PWH and PWoH or between participants with MDD and those without it. Using the songbird model, we computed the log ratio of the highest and lowest 30% of the ranked classes associated with HIV and MDD. We found that HIV infection and lifetime MDD were enriched in a set of differentially abundant inflammatory classes, such as Flavobacteria and Nitrospira. Our results suggest that the circulating plasma microbiome may increase the risk of MDD related to dysbiosis-induced inflammation in PWH. If confirmed, these findings may indicate new biological mechanisms that could be targeted to improve treatment of MDD in PWH.

Oral mucosal breaks trigger anti-citrullinated bacterial and human protein antibody responses in rheumatoid arthritis

Periodontal disease is more common in individuals with rheumatoid arthritis (RA) who have detectable anti-citrullinated protein antibodies (ACPAs), implicating oral mucosal inflammation in RA pathogenesis. Here, we performed paired analysis of human and bacterial transcriptomics in longitudinal blood samples from RA patients. We found that patients with RA and periodontal disease experienced repeated oral bacteremias associated with transcriptional signatures of ISG15+HLADRhi and CD48highS100A2pos monocytes, recently identified in inflamed RA synovia and blood of those with RA flares. The oral bacteria observed transiently in blood were broadly citrullinated in the mouth, and their in situ citrullinated epitopes were targeted by extensively somatically hypermutated ACPAs encoded by RA blood plasmablasts. Together, these results suggest that (i) periodontal disease results in repeated breaches of the oral mucosa that release citrullinated oral bacteria into circulation, which (ii) activate inflammatory monocyte subsets that are observed in inflamed RA synovia and blood of RA patients with flares and (iii) activate ACPA B cells, thereby promoting affinity maturation and epitope spreading to citrullinated human antigens.

Enhancing untargeted metabolomics using metadata-based source annotation

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.

Pan-cancer analyses reveal cancer-type-specific fungal ecologies and bacteriome interactions

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.

Utilizing stability criteria in choosing feature selection methods yields reproducible results in microbiome data

Feature selection is indispensable in microbiome data analysis, but it can be particularly challenging as microbiome data sets are high dimensional, underdetermined, sparse and compositional. Great efforts have recently been made on developing new methods for feature selection that handle the above data characteristics, but almost all methods were evaluated based on performance of model predictions. However, little attention has been paid to address a fundamental question: how appropriate are those evaluation criteria? Most feature selection methods often control the model fit, but the ability to identify meaningful subsets of features cannot be evaluated simply based on the prediction accuracy. If tiny changes to the data would lead to large changes in the chosen feature subset, then many selected features are likely to be a data artifact rather than real biological signal. This crucial need of identifying relevant and reproducible features motivated the reproducibility evaluation criterion such as Stability, which quantifies how robust a method is to perturbations in the data. In our paper, we compare the performance of popular model prediction metrics (MSE or AUC) with proposed reproducibility criterion Stability in evaluating four widely used feature selection methods in both simulations and experimental microbiome applications with continuous or binary outcomes. We conclude that Stability is a preferred feature selection criterion over model prediction metrics because it better quantifies the reproducibility of the feature selection method.

Multiomic Analyses of Nascent Preterm Infant Microbiomes Differentiation Suggest Opportunities for Targeted Intervention

The first week after birth is a critical time for the establishment of microbial communities for infants. Preterm infants face unique environmental impacts on their newly acquired microbiomes, including increased incidence of cesarean section delivery and exposure to antibiotics as well as delayed enteral feeding and reduced human interaction during their intensive care unit stay. Using contextualized paired metabolomics and 16S sequencing data, the development of the gut, skin, and oral microbiomes of infants is profiled daily for the first week after birth, and it is found that the skin microbiome appears robust to early life perturbation, while direct exposure of infants to antibiotics, rather than presumed maternal transmission, delays microbiome development and prevents the early differentiation based on body site regardless of delivery mode. Metabolomic analyses identify the development of all gut metabolomes of preterm infants toward full-term infant profiles, but a significant increase of primary bile acid metabolism only in the non-antibiotic treated vaginally birthed late preterm infants. This study provides a framework for future multi-omic, multibody site analyses on these high-risk preterm infant populations and suggests opportunities for monitoring and intervention, with infant antibiotic exposure as the primary driver of delays in microbiome development.

Abstract 3054: Pan-cancer characterization of the tumor mycobiome and its clinical effects

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.

A comparison of six DNA extraction protocols for 16S, ITS and shotgun metagenomic sequencing of microbial communities

Microbial communities contain a broad phylogenetic diversity of organisms; however, the majority of methods center on describing bacteria and archaea. Fungi are important symbionts in many ecosystems and are potentially important members of the human microbiome, beyond those that can cause disease. To expand our analysis of microbial communities to include data from the fungal internal transcribed spacer (ITS) region, five candidate DNA extraction kits were compared against our standardized protocol for describing bacteria and archaea using 16S rRNA gene amplicon- and shotgun metagenomics sequencing. The results are presented considering a diverse panel of host-associated and environmental sample types and comparing the cost, processing time, well-to-well contamination, DNA yield, limit of detection and microbial community composition among protocols. Across all criteria, the MagMAX Microbiome kit was found to perform best. The PowerSoil Pro kit performed comparably but with increased cost per sample and overall processing time. The Zymo MagBead, NucleoMag Food and Norgen Stool kits were included.

Phylogeny-Aware Analysis of Metagenome Community Ecology Based on Matched Reference Genomes while Bypassing Taxonomy

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.

The Host-Microbiome Response to Hyperbaric Oxygen Therapy in Ulcerative Colitis Patients

BACKGROUND & AIMS

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Salivary bacterial signatures in depression-obesity comorbidity are associated with neurotransmitters and neuroactive dipeptides

BACKGROUND

Depression and obesity are highly prevalent, often co-occurring conditions marked by inflammation. Microbiome perturbations are implicated in obesity-inflammation-depression interrelationships, but how the microbiome mechanistically contributes to pathology remains unclear. Metabolomic investigations into microbial neuroactive metabolites may offer mechanistic insights into host-microbe interactions. Using 16S sequencing and untargeted mass spectrometry of saliva, and blood monocyte inflammation regulation assays, we identified key microbes, metabolites and host inflammation in association with depressive symptomatology, obesity, and depressive symptomatology-obesity comorbidity.

RESULTS

Gram-negative bacteria with inflammation potential were enriched relative to Gram-positive bacteria in comorbid obesity-depression, supporting the inflammation-oral microbiome link in obesity-depression interrelationships. Oral microbiome was more highly predictive of depressive symptomatology-obesity co-occurrences than of obesity or depressive symptomatology independently, suggesting specific microbial signatures associated with obesity-depression co-occurrences. Mass spectrometry analysis revealed significant changes in levels of signaling molecules of microbiota, microbial or dietary derived signaling peptides and aromatic amino acids among depressive symptomatology, obesity and comorbid obesity-depression. Furthermore, integration of the microbiome and metabolomics data revealed that key oral microbes, many previously shown to have neuroactive potential, co-occurred with potential neuropeptides and biosynthetic precursors of the neurotransmitters dopamine, epinephrine and serotonin.

CONCLUSIONS

Together, our findings offer novel insights into oral microbial-brain connection and potential neuroactive metabolites involved.

The microbiome and prostate cancer

There is growing evidence that the microbiome is involved in development and treatment of many human diseases, including prostate cancer. There are several potential pathways for microbiome-based mechanisms for the development of prostate cancer: direct impacts of microbes or microbial products in the prostate or the urine, and indirect impacts from microbes or microbial products in the gastrointestinal tract. Unique microbial signatures have been identified within the stool, oral cavity, tissue, urine, and blood of prostate cancer patients, but studies vary in their findings. Recent studies describe potential diagnostic and therapeutic applications of the microbiome, but further clinical investigation is needed. In this review, we explore the existing literature on the discovery of the human microbiome and its relationship to prostate cancer.

Class Ib MHC-Mediated Immune Interactions Play a Critical Role in Maintaining Mucosal Homeostasis in the Mammalian Large Intestine

Lymphocytes within the intestinal epithelial layer (IEL) in mammals have unique composition compared with their counterparts in the lamina propria. Little is known about the role of some of the key colonic IEL subsets, such as TCRαβ⁺CD8⁺ T cells, in inflammation. We have recently described liver-enriched innate-like TCRαβ⁺CD8αα regulatory T cells, partly controlled by the non-classical MHC molecule, Qa-1^b, that upon adoptive transfer protect from T cell-induced colitis. In this study, we found that TCRαβ⁺CD8αα T cells are reduced among the colonic IEL during inflammation, and that their activation with an agonistic peptide leads to significant Qa-1^b-dependent protection in an acute model of colitis. Cellular expression of Qa-1^b during inflammation and corresponding dependency in peptide-mediated protection suggest that Batf3-dependent CD103⁺CD11b^- type 1 conventional dendritic cells control the protective function of TCRαβ⁺CD8αα T cells in the colonic epithelium. In the colitis model, expression of the potential barrier-protective gene, Muc2, is enhanced upon administration of a Qa-1b agonistic peptide. Notably, in steady state, the mucin metabolizing Akkermansia muciniphila was found in significantly lower abundance amid a dramatic change in overall microbiome and metabolome, increased IL-6 in explant culture, and enhanced sensitivity to dextran sulfate sodium in Qa-1b deficiency. Finally, in patients with inflammatory bowel disease, we found upregulation of HLA-E, a Qa-1^b analog with inflammation and biologic non-response, in silico, suggesting the importance of this regulatory mechanism across species.

Efficient computation of Faith's phylogenetic diversity with applications in characterizing microbiomes

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.

Gastrointestinal Surgery for Inflammatory Bowel Disease Persistently Lowers Microbiome and Metabolome Diversity

BACKGROUND

Many studies have investigated the role of the microbiome in inflammatory bowel disease (IBD), but few have focused on surgery specifically or its consequences on the metabolome that may differ by surgery type and require longitudinal sampling. Our objective was to characterize and contrast microbiome and metabolome changes after different surgeries for IBD, including ileocolonic resection and colectomy.

METHODS

The UC San Diego IBD Biobank was used to prospectively collect 332 stool samples from 129 subjects (50 ulcerative colitis; 79 Crohn's disease). Of these, 21 with Crohn's disease had ileocolonic resections, and 17 had colectomies. We used shotgun metagenomics and untargeted liquid chromatography followed by tandem mass spectrometry metabolomics to characterize the microbiomes and metabolomes of these patients up to 24 months after the initial sampling.

RESULTS

The species diversity and metabolite diversity both differed significantly among groups (species diversity: Mann-Whitney U test P value = 7.8e-17; metabolomics, P-value = 0.0043). Escherichia coli in particular expanded dramatically in relative abundance in subjects undergoing surgery. The species profile was better able to classify subjects according to surgery status than the metabolite profile (average precision 0.80 vs 0.68).

CONCLUSIONS

Intestinal surgeries seem to reduce the diversity of the gut microbiome and metabolome in IBD patients, and these changes may persist. Surgery also further destabilizes the microbiome (but not the metabolome) over time, even relative to the previously established instability in the microbiome of IBD patients. These long-term effects and their consequences for health outcomes need to be studied in prospective longitudinal trials linked to microbiome-involved phenotypes.

A comparison of DNA/RNA extraction protocols for high-throughput sequencing of microbial communities

One goal of microbial ecology researchers is to capture the maximum amount of information from all organisms in a sample. The recent COVID-19 pandemic, caused by the RNA virus SARS-CoV-2, has highlighted a gap in traditional DNA-based protocols, including the high-throughput methods the authors previously established as field standards. To enable simultaneous SARS-CoV-2 and microbial community profiling, the authors compared the relative performance of two total nucleic acid extraction protocols with the authors' previously benchmarked protocol. The authors included a diverse panel of environmental and host-associated sample types, including body sites commonly swabbed for COVID-19 testing. Here the authors present results comparing the cost, processing time, DNA and RNA yield, microbial community composition, limit of detection and well-to-well contamination between these protocols.

Context-aware dimensionality reduction deconvolutes gut microbial community dynamics

The translational power of human microbiome studies is limited by high interindividual variation. We describe a dimensionality reduction tool, compositional tensor factorization (CTF), that incorporates information from the same host across multiple samples to reveal patterns driving differences in microbial composition across phenotypes. CTF identifies robust patterns in sparse compositional datasets, allowing for the detection of microbial changes associated with specific phenotypes that are reproducible across datasets.

Feasibility of using alternative swabs and storage solutions for paired SARS-CoV-2 detection and microbiome analysis in the hospital environment

BACKGROUND

Determining the role of fomites in the transmission of SARS-CoV-2 is essential in the hospital setting and will likely be important outside of medical facilities as governments around the world make plans to ease COVID-19 public health restrictions and attempt to safely reopen economies. Expanding COVID-19 testing to include environmental surfaces would ideally be performed with inexpensive swabs that could be transported safely without concern of being a source of new infections. However, CDC-approved clinical-grade sampling supplies and techniques using a synthetic swab are expensive, potentially expose laboratory workers to viable virus and prohibit analysis of the microbiome due to the presence of antibiotics in viral transport media (VTM). To this end, we performed a series of experiments comparing the diagnostic yield using five consumer-grade swabs (including plastic and wood shafts and various head materials including cotton, synthetic, and foam) and one clinical-grade swab for inhibition to RNA. For three of these swabs, we evaluated performance to detect SARS-CoV-2 in twenty intensive care unit (ICU) hospital rooms of patients including COVID-19+ patients. All swabs were placed in 95% ethanol and further evaluated in terms of RNase activity. SARS-CoV-2 was measured both directly from the swab and from the swab eluent.

RESULTS

Compared to samples collected in VTM, 95% ethanol demonstrated significant inhibition properties against RNases. When extracting directly from the swab head as opposed to the eluent, RNA recovery was approximately 2-4× higher from all six swab types tested as compared to the clinical standard of testing the eluent from a CDC-approved synthetic (SYN) swab. The limit of detection (LoD) of SARS-CoV-2 from floor samples collected using the consumer-grade plastic (CGp) or research-grade plastic The Microsetta Initiative (TMI) swabs was similar or better than the SYN swab, further suggesting that swab type does not impact RNA recovery as measured by the abundance of SARS-CoV-2. The LoD for TMI was between 0 and 362.5 viral particles, while SYN and CGp were both between 725 and 1450 particles. Lastly microbiome analyses (16S rRNA gene sequencing) of paired samples (nasal and floor from same patient room) collected using different swab types in triplicate indicated that microbial communities were not impacted by swab type, but instead driven by the patient and sample type.

CONCLUSIONS

Compared to using a clinical-grade synthetic swab, detection of SARS-CoV-2 from environmental samples collected from ICU rooms of patients with COVID was similar using consumer-grade swabs, stored in 95% ethanol. The yield was best from the swab head rather than the eluent and the low level of RNase activity and lack of antibiotics in these samples makes it possible to perform concomitant microbiome analyses. Video abstract.

Feasibility of using alternative swabs and storage solutions for paired SARS-CoV-2 detection and microbiome analysis in the hospital environment

Background Determining the role of fomites in the transmission of SARS-CoV-2 is essential in the hospital setting and will likely be important outside of medical facilities as governments around the world make plans to ease COVID-19 public health restrictions and attempt to safely reopen economies. Expanding COVID-19 testing to include environmental surfaces would ideally be performed with inexpensive swabs that could be transported safely without concern of being a source of new infections. However, CDC-approved clinical-grade sampling supplies and techniques using a synthetic swab are expensive, potentially expose laboratory workers to viable virus and prohibit analysis of the microbiome due to the presence of antibiotics in viral transport media (VTM). To this end, we performed a series of experiments comparing the diagnostic yield using five consumer-grade swabs (including plastic and wood shafts and various head materials including cotton, synthetic, and foam) and one clinical grade swab for inhibition to RNA. For three of these swabs, we evaluated performance to detect SARS-CoV-2 in twenty intensive care unit (ICU) hospital rooms of patients including COVID-19+ patients. All swabs were placed in 95% ethanol and further evaluated in terms of RNase activity. SARS-CoV-2 was measured both directly from the swab and from the swab eluent. Results Compared to samples collected in VTM, 95% ethanol demonstrated significant inhibition properties against RNases. When extracting directly from the swab head as opposed to the eluent, RNA recovery was approximately 2-4x higher from all six swab types tested as compared to the clinical standard of testing the eluent from a CDC-approved synthetic (SYN) swab. The limit of detection (LoD) of SARSSARS-CoV-2 from floor samples collected using the consumer-grade plastic (CGp) or research-grade plastic The Microsetta Initiative (TMI) swabs was similar or better than the SYN swab, further suggesting that swab type does not impact RNA recovery as measured by the abundance of SARSSARS-CoV-2. The LoD for TMI was between 0-362.5 viral particles while SYN and CGp were both between 725-1450 particles. Lastly microbiome analyses (16S rRNA gene sequencing) of paired samples (nasal and floor from same patient-room) collected using different swab types in triplicate indicated that microbial communities were not impacted by swab type, but instead driven by the patient and sample type. Conclusions Compared to using a clinical-grade synthetic swab, detection of SARS-CoV-2 from environmental samples collected from ICU rooms of patients with COVID was similar using consumer grade swabs, stored in 95% ethanol. The yield was best from the swab head rather than the eluent and the low level of RNase activity and lack of antibiotics in these samples makes it possible to perform concomitant microbiome analyses.

Feasibility of using alternative swabs and storage solutions for paired SARS-CoV-2 detection and microbiome analysis in the hospital environment

Background

Determining the role of fomites in the transmission of SARS-CoV-2 is essential in the hospital setting and will likely be important outside of medical facilities as governments around the world make plans to ease COVID-19 public health restrictions and attempt to safely reopen economies. Expanding COVID-19 testing to include environmental surfaces would ideally be performed with inexpensive swabs that could be transported safely without concern of being a source of new infections. However, CDC-approved clinical-grade sampling supplies and techniques using a synthetic swab are expensive, potentially expose laboratory workers to viable virus and prohibit analysis of the microbiome due to the presence of antibiotics in viral transport media (VTM). To this end, we performed a series of experiments comparing the diagnostic yield using five consumer-grade swabs (including plastic and wood shafts and various head materials including cotton, synthetic, and foam) and one clinical grade swab for inhibition to RNA. For three of these swabs, we evaluated performance to detect SARS-CoV-2 in twenty intensive care unit (ICU) hospital rooms of patients including COVID-19+ patients. All swabs were placed in 95% ethanol and further evaluated in terms of RNase activity. SARS-CoV-2 was measured both directly from the swab and from the swab eluent.

Results

Compared to samples collected in VTM, 95% ethanol demonstrated significant inhibition properties against RNases. When extracting directly from the swab head as opposed to the eluent, RNA recovery was approximately 2-4x higher from all six swab types tested as compared to the clinical standard of testing the eluent from a CDC-approved synthetic (SYN) swab. The limit of detection (LoD) of SARSSARS-CoV-2 from floor samples collected using the consumer-grade plastic (CGp) or research-grade plastic The Microsetta Initiative (TMI) swabs was similar or better than the SYN swab, further suggesting that swab type does not impact RNA recovery as measured by the abundance of SARSSARS-CoV-2. The LoD for TMI was between 0-362.5 viral particles while SYN and CGp were both between 725-1450 particles. Lastly microbiome analyses (16S rRNA gene sequencing) of paired samples (nasal and floor from same patient-room) collected using different swab types in triplicate indicated that microbial communities were not impacted by swab type, but instead driven by the patient and sample type.

Conclusions

Compared to using a clinical-grade synthetic swab, detection of SARS-CoV-2 from environmental samples collected from ICU rooms of patients with COVID was similar using consumer grade swabs, stored in 95% ethanol. The yield was best from the swab head rather than the eluent and the low level of RNase activity and lack of antibiotics in these samples makes it possible to perform concomitant microbiome analyses.

A comparison of DNA/RNA extraction protocols for high-throughput sequencing of microbial communities

One goal among microbial ecology researchers is to capture the maximum amount of information from all organisms in a sample. The recent COVID-19 pandemic, caused by the RNA virus SARS-CoV-2, has highlighted a gap in traditional DNA-based protocols, including the high-throughput methods we previously established as field standards. To enable simultaneous SARS-CoV-2 and microbial community profiling, we compare the relative performance of two total nucleic acid extraction protocols and our previously benchmarked protocol. We included a diverse panel of environmental and host-associated sample types, including body sites commonly swabbed for COVID-19 testing. Here we present results comparing the cost, processing time, DNA and RNA yield, microbial community composition, limit of detection, and well-to-well contamination, between these protocols.

Accession numbers

Raw sequence data were deposited at the European Nucleotide Archive (accession#: ERP124610) and raw and processed data are available at Qiita (Study ID: 12201). All processing and analysis code is available on GitHub ( github.com/justinshaffer/Extraction_test_MagMAX ).

Methods summary

To allow for downstream applications involving RNA-based organisms such as SARS-CoV-2, we compared the two extraction protocols designed to extract DNA and RNA against our previously established protocol for extracting only DNA for microbial community analyses. Across 10 diverse sample types, one of the two protocols was equivalent or better than our established DNA-based protocol. Our conclusion is based on per-sample comparisons of DNA and RNA yield, the number of quality sequences generated, microbial community alpha- and beta-diversity and taxonomic composition, the limit of detection, and extent of well-to-well contamination.

Bacterial modification of the host glycosaminoglycan heparan sulfate modulates SARS-CoV-2 infectivity

The human microbiota has a close relationship with human disease and it remodels components of the glycocalyx including heparan sulfate (HS). Studies of the severe acute respiratory syndrome coronavirus (SARS-CoV-2) spike protein receptor binding domain suggest that infection requires binding to HS and angiotensin converting enzyme 2 (ACE2) in a codependent manner. Here, we show that commensal host bacterial communities can modify HS and thereby modulate SARS-CoV-2 spike protein binding and that these communities change with host age and sex. Common human-associated commensal bacteria whose genomes encode HS-modifying enzymes were identified. The prevalence of these bacteria and the expression of key microbial glycosidases in bronchoalveolar lavage fluid (BALF) was lower in adult COVID-19 patients than in healthy controls. The presence of HS-modifying bacteria decreased with age in two large survey datasets, FINRISK 2002 and American Gut, revealing one possible mechanism for the observed increase in COVID-19 susceptibility with age. In vitro , bacterial glycosidases from unpurified culture media supernatants fully blocked SARS-CoV-2 spike binding to human H1299 protein lung adenocarcinoma cells. HS-modifying bacteria in human microbial communities may regulate viral adhesion, and loss of these commensals could predispose individuals to infection. Understanding the impact of shifts in microbial community composition and bacterial lyases on SARS-CoV-2 infection may lead to new therapeutics and diagnosis of susceptibility.

Feasibility of using alternative swabs and storage solutions for paired SARS-CoV-2 detection and microbiome analysis in the hospital environment

Background Determining the role of fomites in the transmission of SARS-CoV-2 is essential in the hospital setting and will likely be important outside of medical facilities as governments around the world make plans to ease COVID-19 public health restrictions and attempt to safely reopen economies. Expanding COVID-19 testing to include environmental surfaces would ideally be performed with inexpensive swabs that could be transported safely without concern of being a source of new infections. However, CDC-approved clinical-grade sampling supplies and techniques using a synthetic swab are expensive, potentially expose laboratory workers to viable virus and prohibit analysis of the microbiome due to the presence of antibiotics in viral transport media (VTM). To this end, we performed a series of experiments comparing the diagnostic yield using five consumer-grade swabs (including plastic and wood shafts and various head materials including cotton, synthetic, and foam) and one clinical grade swab for inhibition to RNA. For three of these swabs, we evaluated performance to detect SARS-CoV-2 in twenty intensive care unit (ICU) hospital rooms of patients with 16 COVID-19+. All swabs were placed in 95% ethanol and further evaluated in terms of RNase activity. SARS-CoV-2 was measured both directly from the swab and from the swab eluent. Results Compared to samples collected in VTM, 95% ethanol demonstrated significant inhibition properties against RNases. When extracting directly from the swab head as opposed to the eluent, RNA recovery was approximately 2-4x higher from all six swab types tested as compared to the clinical standard of testing the eluent from a CDC-approved synthetic swab. The limit of detection (LoD) of SARs-CoV-2 from floor samples collected using the CGp or TMI swabs was similar or better than the CDC standard, further suggesting that swab type does not impact RNA recovery as measured by SARs-CoV-2. The LoD for TMI was between 0-362.5 viral particles while SYN and CGp were both between 725-1450 particles. Lastly microbiome analyses (16S rRNA) of paired samples (e.g., environment to host) collected using different swab types in triplicate indicated that microbial communities were not impacted by swab type but instead driven by the patient and sample type (floor or nasal). Conclusions Compared to using a clinical-grade synthetic swab, detection of SARS-CoV-2 from environmental samples collected from ICU rooms of patients with COVID was similar using consumer grade swabs, stored in 95% ethanol. The yield was best from the swab head rather than the eluent and the low level of RNase activity in these samples makes it possible to perform concomitant microbiome analysis. Keywords: COVID-19, SARS-CoV-2, RT-qPCR, swab, global health.

Microbiome analyses of blood and tissues suggest cancer diagnostic approach

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.

Human Skin, Oral, and Gut Microbiomes Predict Chronological Age

Human gut microbiomes are known to change with age, yet the relative value of human microbiomes across the body as predictors of age, and prediction robustness across populations is unknown. In this study, we tested the ability of the oral, gut, and skin (hand and forehead) microbiomes to predict age in adults using random forest regression on data combined from multiple publicly available studies, evaluating the models in each cohort individually. Intriguingly, the skin microbiome provides the best prediction of age (mean ± standard deviation, 3.8 ± 0.45 years, versus 4.5 ± 0.14 years for the oral microbiome and 11.5 ± 0.12 years for the gut microbiome). This also agrees with forensic studies showing that the skin microbiome predicts postmortem interval better than microbiomes from other body sites. Age prediction models constructed from the hand microbiome generalized to the forehead and vice versa, across cohorts, and results from the gut microbiome generalized across multiple cohorts (United States, United Kingdom, and China). Interestingly, taxa enriched in young individuals (18 to 30 years) tend to be more abundant and more prevalent than taxa enriched in elderly individuals (>60 yrs), suggesting a model in which physiological aging occurs concomitantly with the loss of key taxa over a lifetime, enabling potential microbiome-targeted therapeutic strategies to prevent aging.IMPORTANCE Considerable evidence suggests that the gut microbiome changes with age or even accelerates aging in adults. Whether the age-related changes in the gut microbiome are more or less prominent than those for other body sites and whether predictions can be made about a person's age from a microbiome sample remain unknown. We therefore combined several large studies from different countries to determine which body site's microbiome could most accurately predict age. We found that the skin was the best, on average yielding predictions within 4 years of chronological age. This study sets the stage for future research on the role of the microbiome in accelerating or decelerating the aging process and in the susceptibility for age-related diseases.

Challenges in the construction of knowledge bases for human microbiome-disease associations

The last few years have seen tremendous growth in human microbiome research, with a particular focus on the links to both mental and physical health and disease. Medical and experimental settings provide initial sources of information about these links, but individual studies produce disconnected pieces of knowledge bounded in context by the perspective of expert researchers reading full-text publications. Building a knowledge base (KB) consolidating these disconnected pieces is an essential first step to democratize and accelerate the process of accessing the collective discoveries of human disease connections to the human microbiome. In this article, we survey the existing tools and development efforts that have been produced to capture portions of the information needed to construct a KB of all known human microbiome-disease associations and highlight the need for additional innovations in natural language processing (NLP), text mining, taxonomic representations, and field-wide vocabulary standardization in human microbiome research. Addressing these challenges will enable the construction of KBs that help identify new insights amenable to experimental validation and potentially clinical decision support.

Author Correction: Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Age- and Sex-Dependent Patterns of Gut Microbial Diversity in Human Adults

Gut microbial diversity changes throughout the human life span and is known to be associated with host sex. We investigated the association of age, sex, and gut bacterial alpha diversity in three large cohorts of adults from four geographical regions: subjects from the United States and United Kingdom in the American Gut Project (AGP) citizen-science initiative and two independent cohorts of Colombians and Chinese. In three of the four cohorts, we observed a strong positive association between age and alpha diversity in young adults that plateaued after age 40 years. We also found sex-dependent differences that were more pronounced in younger adults than in middle-aged adults, with women having higher alpha diversity than men. In contrast to the other three cohorts, no association of alpha diversity with age or sex was observed in the Chinese cohort. The association of alpha diversity with age and sex remained after adjusting for cardiometabolic parameters in the Colombian cohort and antibiotic usage in the AGP cohort. We further attempted to predict the microbiota age in individuals using a machine-learning approach for the men and women in each cohort. Consistent with our alpha-diversity-based findings, U.S. and U.K. women had a significantly higher predicted microbiota age than men, with a reduced difference being seen above age 40 years. This difference was not observed in the Colombian cohort and was observed only in middle-aged Chinese adults. Together, our results provide new insights into the influence of age and sex on the biodiversity of the human gut microbiota during adulthood while highlighting similarities and differences across diverse cohorts. IMPORTANCE Microorganisms in the human gut play a role in health and disease, and in adults higher gut biodiversity has been linked to better health. Since gut microorganisms may be pivotal in the development of microbial therapies, understanding the factors that shape gut biodiversity is of utmost interest. We performed large-scale analyses of the relationship of age and sex to gut bacterial diversity in adult cohorts from four geographic regions: the United States, the United Kingdom, Colombia, and China. In the U.S., U.K., and Colombian cohorts, bacterial biodiversity correlated positively with age in young adults but plateaued at about age 40 years, with no positive association being found in middle-aged adults. Young, but not middle-aged, adult women had higher gut bacterial diversity than men, a pattern confirmed via supervised machine learning. Interestingly, in the Chinese cohort, minimal associations were observed between gut biodiversity and age or sex. Our results highlight the patterns of adult gut biodiversity and provide a framework for future research.

Differential Activation of Hepatic Invariant NKT Cell Subsets Plays a Key Role in Progression of Nonalcoholic Steatohepatitis

Innate immune mechanisms play an important role in inflammatory chronic liver diseases. In this study, we investigated the role of type I or invariant NKT (iNKT) cell subsets in the progression of nonalcoholic steatohepatitis (NASH). We used α-galactosylceramide/CD1d tetramers and clonotypic mAb together with intracytoplasmic cytokine staining to analyze iNKT cells in choline-deficient l-amino acid-defined (CDAA)-induced murine NASH model and in human PBMCs, respectively. Cytokine secretion of hepatic iNKT cells in CDAA-fed C57BL/6 mice altered from predominantly IL-17⁺ to IFN-γ⁺ and IL-4⁺ during NASH progression along with the downmodulation of TCR and NK1.1 expression. Importantly, steatosis, steatohepatitis, and fibrosis were dependent upon the presence of iNKT cells. Hepatic stellate cell activation and infiltration of neutrophils, Kupffer cells, and CD8⁺ T cells as well as expression of key proinflammatory and fibrogenic genes were significantly blunted in Jα18^-/- mice and in C57BL/6 mice treated with an iNKT-inhibitory RAR-γ agonist. Gut microbial diversity was significantly impacted in Jα18^-/- and in CDAA diet-fed mice. An increased frequency of CXCR3⁺IFN-γ⁺T-bet⁺ and IL-17A⁺ iNKT cells was found in PBMC from NASH patients in comparison with nonalcoholic fatty liver patients or healthy controls. Consistent with their in vivo activation, iNKT cells from NASH patients remained hyporesponsive to ex-vivo stimulation with α-galactosylceramide. Accumulation of plasmacytoid dendritic cells in both mice and NASH patients suggest their role in activation of iNKT cells. In summary, our findings indicate that the differential activation of iNKT cells play a key role in mediating diet-induced hepatic steatosis and fibrosis in mice and its potential involvement in NASH progression in humans.

Qiita: rapid, web-enabled microbiome meta-analysis

Multi-omic insights into microbiome function and composition typically advance one study at a time. However, to understand relationships across studies, they must be aggregated into meta-analyses. This makes it possible to generate new hypotheses by finding features that are reproducible across biospecimens and data layers. Qiita dramatically accelerates such integration tasks in a web-based microbiome comparison platform, which we demonstrate with Human Microbiome Project and iHMP data.

Detection and correction of artefacts in estimation of rare copy number variants and analysis of rare deletions in type 1 diabetes

Copy number variants (CNVs) have been proposed as a possible source of ‘missing heritability’ in complex human diseases. Two studies of type 1 diabetes (T1D) found null associations with common copy number polymorphisms, but CNVs of low frequency and high penetrance could still play a role. We used the Log-R-ratio intensity data from a dense single nucleotide polymorphism (SNP) array, ImmunoChip, to detect rare CNV deletions (rDELs) and duplications (rDUPs) in 6808 T1D cases, 9954 controls and 2206 families with T1D-affected offspring. Initial analyses detected CNV associations. However, these were shown to be false-positive findings, failing replication with polymerase chain reaction. We developed a pipeline of quality control (QC) tests that were calibrated using systematic testing of sensitivity and specificity. The case–control odds ratios (OR) of CNV burden on T1D risk resulting from this QC pipeline converged on unity, suggesting no global frequency difference in rDELs or rDUPs. There was evidence that deletions could impact T1D risk for a small minority of cases, with enrichment for rDELs longer than 400 kb (OR = 1.57, P = 0.005). There were also 18 de novo rDELs detected in affected offspring but none for unaffected siblings (P = 0.03). No specific CNV regions showed robust evidence for association with T1D, although frequencies were lower than expected (most less than 0.1%), substantially reducing statistical power, which was examined in detail. We present an R-package, plumbCNV, which provides an automated approach for QC and detection of rare CNVs that can facilitate equivalent analyses of large-scale SNP array datasets.