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 its potential roles in association with clinical outcomes, such as survival, in patients with metastatic disease treated with immune checkpoint inhibitors (ICI). Baseline tumor samples were collected from 71 patients with metastatic melanoma before treatment with ICIs. Bulk RNA sequencing (RNA-seq) was conducted on the formalin-fixed, paraffin-embedded and fresh frozen 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 age of the 71 patients with metastatic melanoma ranged 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 nonresponsive tumors. Responders showed significant enrichment of bacteriophages in the phylum Uroviricota, and nonresponders showed enrichment of several bacteria, including Campylobacter jejuni. 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. SIGNIFICANCE: We analyzed the tumor microbiome and interactions with genes and pathways in metastatic melanoma treated with immunotherapy and identified several microbes associated with immunotherapy response and immune-related gene expression signatures. Machine learning models that combined microbe abundances and gene expression outperformed models using either dataset alone in predicting immunotherapy responses.

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.

Nutrient scavenging-fueled growth in pancreatic cancer depends on caveolae-mediated endocytosis under nutrient-deprived conditions.

Pancreatic ductal adenocarcinoma (PDAC) is characterized by its nutrient-scavenging ability, crucial for tumor progression. Here, we investigated the roles of caveolae-mediated endocytosis (CME) in PDAC progression. Analysis of patient data across diverse datasets revealed a strong association of high caveolin-1 (Cav-1) expression with higher histologic grade, the most aggressive PDAC molecular subtypes, and worse clinical outcomes. Cav-1 loss markedly promoted longer overall and tumor-free survival in a genetically engineered mouse model. Cav-1-deficient tumor cell lines exhibited significantly reduced proliferation, particularly under low nutrient conditions. Supplementing cells with albumin rescued the growth of Cav-1-proficient PDAC cells, but not in Cav-1-deficient PDAC cells under low glutamine conditions. In addition, Cav-1 depletion led to significant metabolic defects, including decreased glycolytic and mitochondrial metabolism, and downstream protein translation signaling pathways. These findings highlight the crucial role of Cav-1 and CME in fueling pancreatic tumorigenesis, sustaining tumor growth, and promoting survival through nutrient scavenging.

Hypoxic adaptation of mitochondria and its impact on tumor cell function.

Mitochondria are the major sink for oxygen in the cell, consuming it during ATP production. Therefore, when environmental oxygen levels drop in the tumor, significant adaptation is required. Mitochondrial activity is also a major producer of biosynthetic precursors and a regulator of cellular oxidative and reductive balance. Because of the complex biochemistry, mitochondrial adaptation to hypoxia occurs through multiple mechanisms and has significant impact on other cellular processes such as macromolecule synthesis and gene regulation. In tumor hypoxia, mitochondria shift their location in the cell and accelerate the fission and quality control pathways. Hypoxic mitochondria also undergo significant changes to fundamental metabolic pathways of carbon metabolism and electron transport. These metabolic changes further impact the nuclear epigenome because mitochondrial metabolites are used as enzymatic substrates for modifying chromatin. This coordinated response delivers physiological flexibility and increased tumor cell robustness during the environmental stress of low oxygen.

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.