Obesity-associated microbiomes instigate visceral adipose tissue inflammation by recruitment of distinct neutrophils.

Neutrophils are increasingly implicated in chronic inflammation and metabolic disorders. Here, we show that visceral adipose tissue (VAT) from individuals with obesity contains more neutrophils than in those without obesity and is associated with a distinct bacterial community. Exploring the mechanism, we gavaged microbiome-depleted mice with stool from patients with and without obesity during high-fat or normal diet administration. Only mice receiving high-fat diet and stool from subjects with obesity show enrichment of VAT neutrophils, suggesting donor microbiome and recipient diet determine VAT neutrophilia. A rise in pro-inflammatory CD4+ Th1 cells and a drop in immunoregulatory T cells in VAT only follows if there is a transient spike in neutrophils. Human VAT neutrophils exhibit a distinct gene expression pattern that is found in different human tissues, including tumors. VAT neutrophils and bacteria may be a novel therapeutic target for treating inflammatory-driven complications of obesity, including insulin resistance and colon cancer.

The Tumor Microbiome Reacts to Hypoxia and Can Influence Response to Radiation Treatment in Colorectal Cancer.

Tumor hypoxia has been shown to predict poor patient outcomes in several cancer types, partially because it reduces radiation's ability to kill cells. We hypothesized that some of the clinical effects of hypoxia could also be due to its impact on the tumor microbiome. Therefore, we examined the RNA sequencing data from the Oncology Research Information Exchange Network database of patients with colorectal cancer treated with radiotherapy. We identified microbial RNAs for each tumor and related them to the hypoxic gene expression scores calculated from host mRNA. Our analysis showed that the hypoxia expression score predicted poor patient outcomes and identified tumors enriched with certain microbes such as Fusobacterium nucleatum. The presence of other microbes, such as Fusobacterium canifelinum, predicted poor patient outcomes, suggesting a potential interaction between hypoxia, the microbiome, and radiation response. To experimentally investigate this concept, we implanted CT26 colorectal cancer cells into immune-competent BALB/c and immune-deficient athymic nude mice. After growth, in which tumors passively acquired microbes from the gastrointestinal tract, we harvested tumors, extracted nucleic acids, and sequenced host and microbial RNAs. We stratified tumors based on their hypoxia score and performed a metatranscriptomic analysis of microbial gene expression. In addition to hypoxia-tropic and -phobic microbial populations, analysis of microbial gene expression at the strain level showed expression differences based on the hypoxia score. Thus, hypoxia gene expression scores seem to associate with different microbial populations and elicit an adaptive transcriptional response in intratumoral microbes, potentially influencing clinical outcomes. SIGNIFICANCE: Tumor hypoxia reduces radiotherapy efficacy. In this study, we explored whether some of the clinical effects of hypoxia could be due to interaction with the tumor microbiome. Hypoxic gene expression scores associated with certain microbes and elicited an adaptive transcriptional response in others that could contribute to poor clinical outcomes.

Comparative analysis of the tumor microbiome, molecular profiles, and immune cell abundances by HPV status in mucosal head and neck cancers and their impact on survival.

Head and Neck Squamous Cell Carcinoma (HNSCC) comprises a diverse group of tumors with variable treatment response and prognosis. The tumor microenvironment (TME), which includes microbiome and immune cells, can impact outcomes. Here, we sought to relate the presence of specific microbes, gene expression, and tumor immune infiltration using tumor transcriptomics from The Cancer Genome Atlas (TCGA) and associate these with overall survival (OS). RNA sequencing (RNAseq) from HNSCC tumors in TCGA was processed through the exogenous sequences in tumors and immune cells (exotic) pipeline to identify and quantify low-abundance microbes. The detection of the Papillomaviridae family of viruses assessed HPV status. All statistical analyses were performed using R. A total of 499 RNAseq samples from TCGA were analyzed. HPV was detected in 111 samples (22%), most commonly Alphapapillomavirus 9 (90.1%). The presence of Alphapapillomavirus 9 was associated with improved OS [HR = 0.60 (95%CI: 0.40-0.89, p = .01)]. Among other microbes, Yersinia pseudotuberculosis was associated with the worst survival (HR = 3.88; p = .008), while Pseudomonas viridiflava had the best survival (HR = 0.05; p = .036). Microbial species found more abundant in HPV- tumors included several gram-negative anaerobes. HPV- tumors had a significantly higher abundance of M0 (p < .001) and M2 macrophages (p = .035), while HPV+ tumors had more T regulatory cells (p < .001) and CD8+ T-cells (p < .001). We identified microbes in HNSCC tumor samples significantly associated with survival. A greater abundance of certain anaerobic microbes was seen in HPV tumors and pro-tumorigenic macrophages. These findings suggest that TME can be used to predict patient outcomes and may help identify mechanisms of resistance to systemic therapies.

A Bioinformatics Tool for Identifying Intratumoral Microbes from the ORIEN Dataset.

Evidence supports significant interactions among microbes, immune cells, and tumor cells in at least 10%-20% of human cancers, emphasizing the importance of further investigating these complex relationships. However, the implications and significance of tumor-related microbes remain largely unknown. Studies have demonstrated the critical roles of host microbes in cancer prevention and treatment responses. Understanding interactions between host microbes and cancer can drive cancer diagnosis and microbial therapeutics (bugs as drugs). Computational identification of cancer-specific microbes and their associations is still challenging due to the high dimensionality and high sparsity of intratumoral microbiome data, which requires large datasets containing sufficient event observations to identify relationships, and the interactions within microbial communities, the heterogeneity in microbial composition, and other confounding effects that can lead to spurious associations. To solve these issues, we present a bioinformatics tool, microbial graph attention (MEGA), to identify the microbes most strongly associated with 12 cancer types. We demonstrate its utility on a dataset from a consortium of nine cancer centers in the Oncology Research Information Exchange Network. This package has three unique features: species-sample relations are represented in a heterogeneous graph and learned by a graph attention network; it incorporates metabolic and phylogenetic information to reflect intricate relationships within microbial communities; and it provides multiple functionalities for association interpretations and visualizations. We analyzed 2,704 tumor RNA sequencing samples and MEGA interpreted the tissue-resident microbial signatures of each of 12 cancer types. MEGA can effectively identify cancer-associated microbial signatures and refine their interactions with tumors. SIGNIFICANCE: Studying the tumor microbiome in high-throughput sequencing data is challenging because of the extremely sparse data matrices, heterogeneity, and high likelihood of contamination. We present a new deep learning tool, MEGA, to refine the organisms that interact with tumors.

Exogenous Sequences in Tumors and Immune Cells (Exotic): A Tool for Estimating the Microbe Abundances in Tumor RNA-seq Data.

The microbiome affects cancer, from carcinogenesis to response to treatments. New evidence suggests that microbes are also present in many tumors, though the scope of how they affect tumor biology and clinical outcomes is in its early stages. A broad survey of tumor microbiome samples across several independent datasets is needed to identify robust correlations for follow-up testing. We created a tool called {exotic} for "exogenous sequences in tumors and immune cells" to carefully identify the tumor microbiome within RNA sequencing (RNA-seq) datasets. We applied it to samples collected through the Oncology Research Information Exchange Network (ORIEN) and The Cancer Genome Atlas. We showed how the processing removes contaminants and batch effects to yield microbe abundances consistent with non-high-throughput sequencing-based approaches and DNA-amplicon-based measurements of a subset of the same tumors. We sought to establish clinical relevance by correlating the microbe abundances with various clinical and tumor measurements, such as age and tumor hypoxia. This process leveraged the two datasets and raised up only the concordant (significant and in the same direction) associations. We observed associations with survival and clinical variables that are cancer specific and relatively few associations with immune composition. Finally, we explored potential mechanisms by which microbes and tumors may interact using a network-based approach. Alistipes, a common gut commensal, showed the highest network degree centrality and was associated with genes related to metabolism and inflammation. The {exotic} tool can support the discovery of microbes in tumors in a way that leverages the many existing and growing RNA-seq datasets. SIGNIFICANCE: The intrinsic tumor microbiome holds great potential for its ability to predict various aspects of cancer biology and as a target for rational manipulation. Here, we describe a tool to quantify microbes from within tumor RNA-seq and apply it to two independent datasets. We show new associations with clinical variables that justify biomarker uses and more experimentation into the mechanisms by which tumor microbiomes affect cancer outcomes.

Mumps epidemic dynamics in the United States before vaccination (1923-1932).

Mumps is a vaccine-preventable, reemerging, and highly transmissible infectious disease. Widespread vaccination dramatically reduced cases; however, case counts have been increasing over the past 20 years. To provide a quantitative overview of historical mumps dynamics that can act as baseline information to help identify causes of mumps reemergence, we analyzed timeseries of cases reported from 1923 to 1932 in the United States. During that time, 239,230 mumps cases were reported in 70 cities. Larger cities reported annual epidemics and smaller cities reported intermittent, sporadic outbreaks. The critical community size above which transmission continuously occurred was likely between 365,583 and 781,188 individuals but could range as high as 3,376,438 individuals. Mumps cases increased as city size increased, suggesting density-dependent transmission. Using a density-dependent SEIR model, we calculated a mean effective reproductive number (Re) of 1.2. Re varied by city and over time, with periodic high values that could characterize short periods of very high transmission known as superspreading events. Case counts most often peaked in March, with higher-than-average transmission from December through April and showed a correlation with weekly births. While certain city pairs in Midwestern states had synchronous outbreaks, most outbreaks were less synchronous and not driven by distance between cities. This work demonstrates the importance of long-term infectious disease surveillance data and will inform future studies on mumps reemergence and control.

The tumor microbiome as a predictor of outcomes in patients with metastatic melanoma treated with immune checkpoint inhibitors.

Emerging evidence supports the important role of the tumor microbiome in oncogenesis, cancer immune phenotype, cancer progression, and treatment outcomes in many malignancies. In this study, we investigated the metastatic melanoma tumor microbiome and potential roles in association with clinical outcomes, such as survival, in patients with metastatic disease treated with immune checkpoint inhibitors (ICIs). Baseline tumor samples were collected from 71 patients with metastatic melanoma before treatment with ICIs. Bulk RNA-seq was conducted on the formalin-fixed paraffin-embedded (FFPE) tumor samples. Durable clinical benefit (primary clinical endpoint) following ICIs was defined as overall survival ≥24 months and no change to the primary drug regimen (responders). We processed RNA-seq reads to carefully identify exogenous sequences using the {exotic} tool. The 71 patients with metastatic melanoma ranged in age from 24 to 83 years, 59% were male, and 55% survived >24 months following the initiation of ICI treatment. Exogenous taxa were identified in the tumor RNA-seq, including bacteria, fungi, and viruses. We found differences in gene expression and microbe abundances in immunotherapy responsive versus non-responsive tumors. Responders showed significant enrichment of several microbes including Fusobacterium nucleatum, and non-responders showed enrichment of fungi, as well as several bacteria. These microbes correlated with immune-related gene expression signatures. Finally, we found that models for predicting prolonged survival with immunotherapy using both microbe abundances and gene expression outperformed models using either dataset alone. Our findings warrant further investigation and potentially support therapeutic strategies to modify the tumor microbiome in order to improve treatment outcomes with ICIs.

A bioinformatics tool for identifying intratumoral microbes from the ORIEN dataset.

Evidence supports significant interactions among microbes, immune cells, and tumor cells in at least 10-20% of human cancers, emphasizing the importance of further investigating these complex relationships. However, the implications and significance of tumor-related microbes remain largely unknown. Studies have demonstrated the critical roles of host microbes in cancer prevention and treatment responses. Understanding interactions between host microbes and cancer can drive cancer diagnosis and microbial therapeutics (bugs as drugs). Computational identification of cancer-specific microbes and their associations is still challenging due to the high dimensionality and high sparsity of intratumoral microbiome data, which requires large datasets containing sufficient event observations to identify relationships, and the interactions within microbial communities, the heterogeneity in microbial composition, and other confounding effects that can lead to spurious associations. To solve these issues, we present a bioinformatics tool, MEGA, to identify the microbes most strongly associated with 12 cancer types. We demonstrate its utility on a dataset from a consortium of 9 cancer centers in the Oncology Research Information Exchange Network (ORIEN). This package has 3 unique features: species-sample relations are represented in a heterogeneous graph and learned by a graph attention network; it incorporates metabolic and phylogenetic information to reflect intricate relationships within microbial communities; and it provides multiple functionalities for association interpretations and visualizations. We analyzed 2704 tumor RNA-seq samples and MEGA interpreted the tissue-resident microbial signatures of each of 12 cancer types. MEGA can effectively identify cancer-associated microbial signatures and refine their interactions with tumors.

Association between Tumor Microbiome and Hypoxia across Anatomic Subsites of Head and Neck Cancers.

Purpose/Objective(s): Microbiome has been shown to affect tumorigenesis by promoting inflammation. However, the association between the upper aerodigestive microbiome and head and neck squamous cell carcinoma (HNSCC) is not well established. Hypoxia is a modifiable factor associated with poor radiation response. Our study analyzed the HNSCC tumor samples from The Cancer Genome Atlas (TCGA) to investigate the relationship between different HNSCC tumor subsites, hypoxia, and local tumor microbiome composition. Results: A total of 357 patients were included [Oral cavity (OC) = 226, Oropharynx (OPx) = 53, and Larynx/Hypopharynx (LHPx) = 78], of which 12.8%, 71.7%, and 10.3%, respectively, were HPV positive. The mean (SD) hypoxia scores were 30.18 (11.10), 24.31 (14.13), and 29.53 (12.61) in OC, OPx, and LHPx tumors, respectively, with higher values indicating greater hypoxia. The hypoxia score was significantly higher for OC tumors compared to OPx (p = 0.044) and LHPx (p = 0.002). There was no significant correlation between hypoxia and HPV status. Pseudomonas sp. in OC, Actinomyces sp. and Sulfurimonas sp. in OPx, and Filifactor, Pseudomonas and Actinomyces sp. in LHPx had the strongest association with the hypoxia score. Materials/Methods: Tumor RNAseq samples from TCGA were processed, and the R package "tmesig" was used to calculate gene expression signature, including the Buffa hypoxia (BH) score, a validated hypoxia signature using 52 hypoxia-regulated genes. Microbe relative abundances were modeled with primary tumor location and a high vs. low tertile BH score applying a gamma-distributed generalized linear regression using the "stats" package in R, with adjusted p-value < 0.05 considered significant. Conclusions: In our study, oral cavity tumors were found to be more hypoxic compared to other head and neck subsites, which could potentially contribute to their radiation resistance. For each subsite, distinct microbial populations were over-represented in hypoxic tumors in a subsite-specific manner. Further studies focusing on an association between microbiome, hypoxia, and patient outcomes are warranted.