TARO: tree-aggregated factor regression for microbiome data integration.

Motivation: Although the human microbiome plays a key role in health and disease, the biological mechanisms underlying the interaction between the microbiome and its host are incompletely understood. Integration with other molecular profiling data offers an opportunity to characterize the role of the microbiome and elucidate therapeutic targets. However, this remains challenging to the high dimensionality, compositionality, and rare features found in microbiome profiling data. These challenges necessitate the use of methods that can achieve structured sparsity in learning cross-platform association patterns. Results: We propose Tree-Aggregated factor RegressiOn (TARO) for the integration of microbiome and metabolomic data. We leverage information on the phylogenetic tree structure to flexibly aggregate rare features. We demonstrate through simulation studies that TARO accurately recovers a low-rank coefficient matrix and identifies relevant features. We applied TARO to microbiome and metabolomic profiles gathered from subjects being screened for colorectal cancer to understand how gut microrganisms shape intestinal metabolite abundances. Availability and implementation: The R package TARO implementing the proposed methods is available online at https://github.com/amishra-stats/taro-package .

Neoadjuvant chemotherapy plus nivolumab with or without ipilimumab in operable non-small cell lung cancer: the phase 2 platform NEOSTAR trial.

Neoadjuvant ipilimumab + nivolumab (Ipi+Nivo) and nivolumab + chemotherapy (Nivo+CT) induce greater pathologic response rates than CT alone in patients with operable non-small cell lung cancer (NSCLC). The impact of adding ipilimumab to neoadjuvant Nivo+CT is unknown. Here we report the results and correlates of two arms of the phase 2 platform NEOSTAR trial testing neoadjuvant Nivo+CT and Ipi+Nivo+CT with major pathologic response (MPR) as the primary endpoint. MPR rates were 32.1% (7/22, 80% confidence interval (CI) 18.7-43.1%) in the Nivo+CT arm and 50% (11/22, 80% CI 34.6-61.1%) in the Ipi+Nivo+CT arm; the primary endpoint was met in both arms. In patients without known tumor EGFR/ALK alterations, MPR rates were 41.2% (7/17) and 62.5% (10/16) in the Nivo+CT and Ipi+Nivo+CT groups, respectively. No new safety signals were observed in either arm. Single-cell sequencing and multi-platform immune profiling (exploratory endpoints) underscored immune cell populations and phenotypes, including effector memory CD8+ T, B and myeloid cells and markers of tertiary lymphoid structures, that were preferentially increased in the Ipi+Nivo+CT cohort. Baseline fecal microbiota in patients with MPR were enriched with beneficial taxa, such as Akkermansia, and displayed reduced abundance of pro-inflammatory and pathogenic microbes. Neoadjuvant Ipi+Nivo+CT enhances pathologic responses and warrants further study in operable NSCLC. (ClinicalTrials.gov registration: NCT03158129 .).

Intestinal toxicity to CTLA-4 blockade driven by IL-6 and myeloid infiltration.

Immune checkpoint blockade (ICB) has revolutionized cancer treatment, yet quality of life and continuation of therapy can be constrained by immune-related adverse events (irAEs). Limited understanding of irAE mechanisms hampers development of approaches to mitigate their damage. To address this, we examined whether mice gained sensitivity to anti-CTLA-4 (αCTLA-4)-mediated toxicity upon disruption of gut homeostatic immunity. We found αCTLA-4 drove increased inflammation and colonic tissue damage in mice with genetic predisposition to intestinal inflammation, acute gastrointestinal infection, transplantation with a dysbiotic fecal microbiome, or dextran sodium sulfate administration. We identified an immune signature of αCTLA-4-mediated irAEs, including colonic neutrophil accumulation and systemic interleukin-6 (IL-6) release. IL-6 blockade combined with antibiotic treatment reduced intestinal damage and improved αCTLA-4 therapeutic efficacy in inflammation-prone mice. Intestinal immune signatures were validated in biopsies from patients with ICB colitis. Our work provides new preclinical models of αCTLA-4 intestinal irAEs, mechanistic insights into irAE development, and potential approaches to enhance ICB efficacy while mitigating irAEs.

Obesity Is Associated with Altered Tumor Metabolism in Metastatic Melanoma.

PURPOSE: Overweight/obese (OW/OB) patients with metastatic melanoma unexpectedly have improved outcomes with immune checkpoint inhibitors (ICI) and BRAF-targeted therapies. The mechanism(s) underlying this association remain unclear, thus we assessed the integrated molecular, metabolic, and immune profile of tumors, as well as gut microbiome features, for associations with patient body mass index (BMI). EXPERIMENTAL DESIGN: Associations between BMI [normal (NL < 25) or OW/OB (BMI ≥ 25)] and tumor or microbiome characteristics were examined in specimens from 782 patients with metastatic melanoma across 7 cohorts. DNA associations were evaluated in The Cancer Genome Atlas cohort. RNA sequencing from 4 cohorts (n = 357) was batch corrected and gene set enrichment analysis (GSEA) by BMI category was performed. Metabolic profiling was conducted in a subset of patients (x = 36) by LC/MS, and in flow-sorted melanoma tumor cells (x = 37) and patient-derived melanoma cell lines (x = 17) using the Seahorse XF assay. Gut microbiome features were examined in an independent cohort (n = 371). RESULTS: DNA mutations and copy number variations were not associated with BMI. GSEA demonstrated that tumors from OW/OB patients were metabolically quiescent, with downregulation of oxidative phosphorylation and multiple other metabolic pathways. Direct metabolite analysis and functional metabolic profiling confirmed decreased central carbon metabolism in OW/OB metastatic melanoma tumors and patient-derived cell lines. The overall structure, diversity, and taxonomy of the fecal microbiome did not differ by BMI. CONCLUSIONS: These findings suggest that the host metabolic phenotype influences melanoma metabolism and provide insight into the improved outcomes observed in OW/OB patients with metastatic melanoma treated with ICIs and targeted therapies. See related commentary by Smalley, p. 5.

Diet-driven microbial ecology underpins associations between cancer immunotherapy outcomes and the gut microbiome.

The gut microbiota shapes the response to immune checkpoint inhibitors (ICIs) in cancer, however dietary and geographic influences have not been well-studied in prospective trials. To address this, we prospectively profiled baseline gut (fecal) microbiota signatures and dietary patterns of 103 trial patients from Australia and the Netherlands treated with neoadjuvant ICIs for high risk resectable metastatic melanoma and performed an integrated analysis with data from 115 patients with melanoma treated with ICIs in the United States. We observed geographically distinct microbial signatures of response and immune-related adverse events (irAEs). Overall, response rates were higher in Ruminococcaceae-dominated microbiomes than in Bacteroidaceae-dominated microbiomes. Poor response was associated with lower fiber and omega 3 fatty acid consumption and elevated levels of C-reactive protein in the peripheral circulation at baseline. Together, these data provide insight into the relevance of native gut microbiota signatures, dietary intake and systemic inflammation in shaping the response to and toxicity from ICIs, prompting the need for further studies in this area.

Negative binomial factor regression with application to microbiome data analysis.

The human microbiome provides essential physiological functions and helps maintain host homeostasis via the formation of intricate ecological host-microbiome relationships. While it is well established that the lifestyle of the host, dietary preferences, demographic background, and health status can influence microbial community composition and dynamics, robust generalizable associations between specific host-associated factors and specific microbial taxa have remained largely elusive. Here, we propose factor regression models that allow the estimation of structured parsimonious associations between host-related features and amplicon-derived microbial taxa. To account for the overdispersed nature of the amplicon sequencing count data, we propose negative binomial reduced rank regression (NB-RRR) and negative binomial co-sparse factor regression (NB-FAR). While NB-RRR encodes the underlying dependency among the microbial abundances as outcomes and the host-associated features as predictors through a rank-constrained coefficient matrix, NB-FAR uses a sparse singular value decomposition of the coefficient matrix. The latter approach avoids the notoriously difficult joint parameter estimation by extracting sparse unit-rank components of the coefficient matrix sequentially, effectively delivering interpretable bi-clusters of taxa and host-associated factors. To solve the nonconvex optimization problems associated with these factor regression models, we present a novel iterative block-wise majorization procedure. Extensive simulation studies and an application to the microbial abundance data from the American Gut Project (AGP) demonstrate the efficacy of the proposed procedure. In the AGP data, we identify several factors that strongly link dietary habits and host life style to specific microbial families.