Medication effects on the gut microbiome in allo-HCT.

Decreased gut microbiome diversity has been associated with negative outcome in allogeneic hematopoietic stem cell transfer (HCT). A study published in this issue of Cell identifies associations between non-antibiotic drug administration, microbiome state transitions, and response to HCT, highlighting the potential impact of such drugs on microbiome and HCT outcome.

Strain level microbial detection and quantification with applications to single cell metagenomics.

Computational identification and quantification of distinct microbes from high throughput sequencing data is crucial for our understanding of human health. Existing methods either use accurate but computationally expensive alignment-based approaches or less accurate but computationally fast alignment-free approaches, which often fail to correctly assign reads to genomes. Here we introduce CAMMiQ, a combinatorial optimization framework to identify and quantify distinct genomes (specified by a database) in a metagenomic dataset. As a key methodological innovation, CAMMiQ uses substrings of variable length and those that appear in two genomes in the database, as opposed to the commonly used fixed-length, unique substrings. These substrings allow to accurately decouple mixtures of highly similar genomes resulting in higher accuracy than the leading alternatives, without requiring additional computational resources, as demonstrated on commonly used benchmarking datasets. Importantly, we show that CAMMiQ can distinguish closely related bacterial strains in simulated metagenomic and real single-cell metatranscriptomic data.

AXL and Error-Prone DNA Replication Confer Drug Resistance and Offer Strategies to Treat EGFR-Mutant Lung Cancer.

Anticancer therapies have been limited by the emergence of mutations and other adaptations. In bacteria, antibiotics activate the SOS response, which mobilizes error-prone factors that allow for continuous replication at the cost of mutagenesis. We investigated whether the treatment of lung cancer with EGFR inhibitors (EGFRi) similarly engages hypermutators. In cycling drug-tolerant persister (DTP) cells and in EGFRi-treated patients presenting residual disease, we observed upregulation of GAS6, whereas ablation of GAS6's receptor, AXL, eradicated resistance. Reciprocally, AXL overexpression enhanced DTP survival and accelerated the emergence of T790M, an EGFR mutation typical to resistant cells. Mechanistically, AXL induces low-fidelity DNA polymerases and activates their organizer, RAD18, by promoting neddylation. Metabolomics uncovered another hypermutator, AXL-driven activation of MYC, and increased purine synthesis that is unbalanced by pyrimidines. Aligning anti-AXL combination treatments with the transition from DTPs to resistant cells cured patient-derived xenografts. Hence, similar to bacteria, tumors tolerate therapy by engaging pharmacologically targetable endogenous mutators. SIGNIFICANCE: EGFR-mutant lung cancers treated with kinase inhibitors often evolve resistance due to secondary mutations. We report that in similarity to the bacterial SOS response stimulated by antibiotics, endogenous mutators are activated in drug-treated cells, and this heralds tolerance. Blocking the process prevented resistance in xenograft models, which offers new treatment strategies. This article is highlighted in the In This Issue feature, p. 2483.

Neoantigen Identification and Response to Adoptive Cell Transfer in Anti-PD-1 Naïve and Experienced Patients with Metastatic Melanoma.

PURPOSE: Immune checkpoint blockade (ICB) agents and adoptive cell transfer (ACT) of tumor-infiltrating lymphocytes (TIL) are prominent immunotherapies used for the treatment of advanced melanoma. Both therapies rely on activation of lymphocytes that target shared tumor antigens or neoantigens. Recent analysis of patients with metastatic melanoma who underwent treatment with TIL ACT at the NCI demonstrated decreased responses in patients previously treated with anti-PD-1 agents. We aimed to find a basis for the difference in response rates between anti-PD-1 naïve and experienced patients. PATIENTS AND METHODS: We examined the tumor mutational burden (TMB) of resected tumors and the repertoire of neoantigens targeted by autologous TIL in a cohort of 112 anti-PD-1 naïve and 69 anti-PD-1 experienced patients. RESULTS: Anti-PD-1 naïve patients were found to possess tumors with higher TMBs (352.0 vs. 213.5, P = 0.005) and received TIL reactive with more neoantigens (2 vs. 1, P = 0.003) compared with anti-PD-1 experienced patients. Among patients treated with TIL ACT, TMB and number of neoantigens identified were higher in ACT responders than ACT nonresponders in both anti-PD-1 naïve and experienced patients. Among patients with comparable TMBs and predicted neoantigen loads, treatment products administered to anti-PD-1 naïve patients were more likely to contain T cells reactive against neoantigens than treatment products for anti-PD-1 experienced patients (2.5 vs. 1, P = 0.02). CONCLUSIONS: These results indicate that decreases in TMB and targeted neoantigens partially account for the difference in response to ACT and that additional factors likely influence responses in these patients. See related commentary by Blass and Ott, p. 2980.

Deconvolving Clinically Relevant Cellular Immune Cross-talk from Bulk Gene Expression Using CODEFACS and LIRICS Stratifies Patients with Melanoma to Anti-PD-1 Therapy.

The tumor microenvironment (TME) is a complex mixture of cell types whose interactions affect tumor growth and clinical outcome. To discover such interactions, we developed CODEFACS (COnfident DEconvolution For All Cell Subsets), a tool deconvolving cell type-specific gene expression in each sample from bulk expression, and LIRICS (Ligand-Receptor Interactions between Cell Subsets), a statistical framework prioritizing clinically relevant ligand-receptor interactions between cell types from the deconvolved data. We first demonstrate the superiority of CODEFACS versus the state-of-the-art deconvolution method CIBERSORTx. Second, analyzing The Cancer Genome Atlas, we uncover cell type-specific ligand-receptor interactions uniquely associated with mismatch-repair deficiency across different cancer types, providing additional insights into their enhanced sensitivity to anti-programmed cell death protein 1 (PD-1) therapy compared with other tumors with high neoantigen burden. Finally, we identify a subset of cell type-specific ligand-receptor interactions in the melanoma TME that stratify survival of patients receiving anti-PD-1 therapy better than some recently published bulk transcriptomics-based methods. SIGNIFICANCE: This work presents two new computational methods that can deconvolve a large collection of bulk tumor gene expression profiles into their respective cell type-specific gene expression profiles and identify cell type-specific ligand-receptor interactions predictive of response to immune-checkpoint blockade therapy. This article is highlighted in the In This Issue feature, p. 873.

Characterization of HLH-like manifestations as a CRS variant in patients receiving CD22 CAR T cells.

Chimeric antigen receptor (CAR) T-cell toxicities resembling hemophagocytic lymphohistiocytosis (HLH) occur in a subset of patients with cytokine release syndrome (CRS). As a variant of conventional CRS, a comprehensive characterization of CAR T-cell-associated HLH (carHLH) and investigations into associated risk factors are lacking. In the context of 59 patients infused with CD22 CAR T cells where a substantial proportion developed carHLH, we comprehensively describe the manifestations and timing of carHLH as a CRS variant and explore factors associated with this clinical profile. Among 52 subjects with CRS, 21 (40.4%) developed carHLH. Clinical features of carHLH included hyperferritinemia, hypertriglyceridemia, hypofibrinogenemia, coagulopathy, hepatic transaminitis, hyperbilirubinemia, severe neutropenia, elevated lactate dehydrogenase, and occasionally hemophagocytosis. Development of carHLH was associated with preinfusion natural killer(NK) cell lymphopenia and higher bone marrow T-cell:NK cell ratio, which was further amplified with CAR T-cell expansion. Following CRS, more robust CAR T-cell and CD8 T-cell expansion in concert with pronounced NK cell lymphopenia amplified preinfusion differences in those with carHLH without evidence for defects in NK cell mediated cytotoxicity. CarHLH was further characterized by persistent elevation of HLH-associated inflammatory cytokines, which contrasted with declining levels in those without carHLH. In the setting of CAR T-cell mediated expansion, clinical manifestations and immunophenotypic profiling in those with carHLH overlap with features of secondary HLH, prompting consideration of an alternative framework for identification and management of this toxicity profile to optimize outcomes following CAR T-cell infusion.

Network-based approaches elucidate differences within APOBEC and clock-like signatures in breast cancer.

BACKGROUND: Studies of cancer mutations have typically focused on identifying cancer driving mutations that confer growth advantage to cancer cells. However, cancer genomes accumulate a large number of passenger somatic mutations resulting from various endogenous and exogenous causes, including normal DNA damage and repair processes or cancer-related aberrations of DNA maintenance machinery as well as mutations triggered by carcinogenic exposures. Different mutagenic processes often produce characteristic mutational patterns called mutational signatures. Identifying mutagenic processes underlying mutational signatures shaping a cancer genome is an important step towards understanding tumorigenesis. METHODS: To investigate the genetic aberrations associated with mutational signatures, we took a network-based approach considering mutational signatures as cancer phenotypes. Specifically, our analysis aims to answer the following two complementary questions: (i) what are functional pathways whose gene expression activities correlate with the strengths of mutational signatures, and (ii) are there pathways whose genetic alterations might have led to specific mutational signatures? To identify mutated pathways, we adopted a recently developed optimization method based on integer linear programming. RESULTS: Analyzing a breast cancer dataset, we identified pathways associated with mutational signatures on both expression and mutation levels. Our analysis captured important differences in the etiology of the APOBEC-related signatures and the two clock-like signatures. In particular, it revealed that clustered and dispersed APOBEC mutations may be caused by different mutagenic processes. In addition, our analysis elucidated differences between two age-related signatures-one of the signatures is correlated with the expression of cell cycle genes while the other has no such correlation but shows patterns consistent with the exposure to environmental/external processes. CONCLUSIONS: This work investigated, for the first time, a network-level association of mutational signatures and dysregulated pathways. The identified pathways and subnetworks provide novel insights into mutagenic processes that the cancer genomes might have undergone and important clues for developing personalized drug therapies.

A Sticky Multinomial Mixture Model of Strand-Coordinated Mutational Processes in Cancer.

The characterization of mutational processes in terms of their signatures of activity relies mostly on the assumption that mutations in a given cancer genome are independent of one another. Recently, it was discovered that certain segments of mutations, termed processive groups, occur on the same DNA strand and are generated by a single process or signature. Here we provide a first probabilistic model of mutational signatures that accounts for their observed stickiness and strand coordination. The model conditions on the observed strand for each mutation and allows the same signature to generate a run of mutations. It can both use known signatures or learn new ones. We show that this model provides a more accurate description of the properties of mutagenic processes than independent-mutation achieving substantially higher likelihood on held-out data. We apply this model to characterize the processivity of mutagenic processes across multiple types of cancer.

Modeling clinical and molecular covariates of mutational process activity in cancer.

MOTIVATION: Somatic mutations result from processes related to DNA replication or environmental/lifestyle exposures. Knowing the activity of mutational processes in a tumor can inform personalized therapies, early detection, and understanding of tumorigenesis. Computational methods have revealed 30 validated signatures of mutational processes active in human cancers, where each signature is a pattern of single base substitutions. However, half of these signatures have no known etiology, and some similar signatures have distinct etiologies, making patterns of mutation signature activity hard to interpret. Existing mutation signature detection methods do not consider tumor-level clinical/demographic (e.g. smoking history) or molecular features (e.g. inactivations to DNA damage repair genes). RESULTS: To begin to address these challenges, we present the Tumor Covariate Signature Model (TCSM), the first method to directly model the effect of observed tumor-level covariates on mutation signatures. To this end, our model uses methods from Bayesian topic modeling to change the prior distribution on signature exposure conditioned on a tumor's observed covariates. We also introduce methods for imputing covariates in held-out data and for evaluating the statistical significance of signature-covariate associations. On simulated and real data, we find that TCSM outperforms both non-negative matrix factorization and topic modeling-based approaches, particularly in recovering the ground truth exposure to similar signatures. We then use TCSM to discover five mutation signatures in breast cancer and predict homologous recombination repair deficiency in held-out tumors. We also discover four signatures in a combined melanoma and lung cancer cohort-using cancer type as a covariate-and provide statistical evidence to support earlier claims that three lung cancers from The Cancer Genome Atlas are misdiagnosed metastatic melanomas. AVAILABILITY AND IMPLEMENTATION: TCSM is implemented in Python 3 and available at https://github.com/lrgr/tcsm, along with a data workflow for reproducing the experiments in the paper. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genome-wide prediction of synthetic rescue mediators of resistance to targeted and immunotherapy.

Most patients with advanced cancer eventually acquire resistance to targeted therapies, spurring extensive efforts to identify molecular events mediating therapy resistance. Many of these events involve synthetic rescue (SR) interactions, where the reduction in cancer cell viability caused by targeted gene inactivation is rescued by an adaptive alteration of another gene (the rescuer). Here, we perform a genome-wide in silico prediction of SR rescuer genes by analyzing tumor transcriptomics and survival data of 10,000 TCGA cancer patients. Predicted SR interactions are validated in new experimental screens. We show that SR interactions can successfully predict cancer patients' response and emerging resistance. Inhibiting predicted rescuer genes sensitizes resistant cancer cells to therapies synergistically, providing initial leads for developing combinatorial approaches to overcome resistance proactively. Finally, we show that the SR analysis of melanoma patients successfully identifies known mediators of resistance to immunotherapy and predicts novel rescuers.

Harnessing synthetic lethality to predict the response to cancer treatment.

While synthetic lethality (SL) holds promise in developing effective cancer therapies, SL candidates found via experimental screens often have limited translational value. Here we present a data-driven approach, ISLE (identification of clinically relevant synthetic lethality), that mines TCGA cohort to identify the most likely clinically relevant SL interactions (cSLi) from a given candidate set of lab-screened SLi. We first validate ISLE via a benchmark of large-scale drug response screens and by predicting drug efficacy in mouse xenograft models. We then experimentally test a select set of predicted cSLi via new screening experiments, validating their predicted context-specific sensitivity in hypoxic vs normoxic conditions and demonstrating cSLi's utility in predicting synergistic drug combinations. We show that cSLi can successfully predict patients' drug treatment response and provide patient stratification signatures. ISLE thus complements existing actionable mutation-based methods for precision cancer therapy, offering an opportunity to expand its scope to the whole genome.

Putative functional genes in idiopathic dilated cardiomyopathy.

Idiopathic dilated cardiomyopathy (DCM) is a complex disorder with a genetic and an environmental component involving multiple genes, many of which are yet to be discovered. We integrate genetic, epigenetic, transcriptomic, phenotypic, and evolutionary features into a method - Hridaya, to infer putative functional genes underlying DCM in a genome-wide fashion, using 213 human heart genomes and transcriptomes. Many genes identified by Hridaya are experimentally shown to cause cardiac complications. We validate the top predicted genes, via five different genome-wide analyses: First, the predicted genes are associated with cardiovascular functions. Second, their knockdowns in mice induce cardiac abnormalities. Third, their inhibition by drugs cause cardiac side effects in human. Fourth, they tend to have differential exon usage between DCM and normal samples. Fifth, analyzing 213 individual genotypes, we show that regulatory polymorphisms of the predicted genes are associated with elevated risk of cardiomyopathy. The stratification of DCM patients based on cardiac expression of the functional genes reveals two subgroups differing in key cardiac phenotypes. Integrating predicted functional genes with cardiomyocyte drug treatment experiments reveals novel potential drug targets. We provide a list of investigational drugs that target the newly identified functional genes that may lead to cardiac side effects.