Cancer treatment by immune checkpoint blockade (ICB) can bring long-lasting clinical benefits, but only a fraction of patients respond to treatment. To predict ICB response, we developed TIDE, a computational method to model two primary mechanisms of tumor immune evasion: the induction of T cell dysfunction in tumors with high infiltration of cytotoxic T lymphocytes (CTL) and the prevention of T cell infiltration in tumors with low CTL level. We identified signatures of T cell dysfunction from large tumor cohorts by testing how the expression of each gene in tumors interacts with the CTL infiltration level to influence patient survival. We also modeled factors that exclude T cell infiltration into tumors using expression signatures from immunosuppressive cells. Using this framework and pre-treatment RNA-Seq or NanoString tumor expression profiles, TIDE predicted the outcome of melanoma patients treated with first-line anti-PD1 or anti-CTLA4 more accurately than other biomarkers such as PD-L1 level and mutation load. TIDE also revealed new candidate ICB resistance regulators, such as SERPINB9, demonstrating utility for immunotherapy research.

Summary Cherubism is an autosomal dominant syndrome characterized by inflammatory destructive bony lesions resulting in symmetrical deformities of the facial bones. Cherubism is caused by mutations in Sh3bp2, the gene that encodes the adaptor protein 3BP2. Most identified mutations in 3BP2 lie within the peptide sequence RSPPDG. A mouse model of cherubism develops hyperactive bone remodelling osteoclasts and systemic inflammation characterized by expansion of the myelomonocytic lineage. The mechanism by which cherubism mutations alter 3BP2 function has remained obscure. Here we show that Tankyrase, a member of the poly(ADP-ribose)polymerase (PARP) family, regulates 3BP2 stability through ADP-ribosylation and subsequent ubiquitylation by the E3-ubiquitin ligase RNF146 in osteoclasts. Cherubism mutations uncouple 3BP2 from Tankyrase-mediated protein destruction, which results in its stabilization and subsequent hyperactivation of the SRC, SYK and VAV signalling pathways.

Tumor-infiltrating B cells are an important component in the microenvironment with unclear anti-tumor impacts. We enhanced our previous computational algorithm TRUST to extract the B cell immunoglobulin (Ig) hypervariable regions from bulk tumor RNA-seq data. TRUST assembled over 30 million complementarity-determining region 3 (CDR3s) of the B cell heavy chain (IgH) from The Cancer Genome Atlas (TCGA). Widespread B cell clonal expansions and Ig subclass switch events were observed in diverse human cancers. Prevalent somatic copy number alterations in MICA and MICB genes related to antibody-dependent cell mediated cytotoxicity (ADCC) were identified in tumors with elevated B cell activity. IgG3-1 subclass switch interacts with the B cell receptor affinity maturation and defects in the ADCC pathway. Comprehensive pan-cancer analyses of tumor-infiltrating B cell receptor repertoires identified novel tumor immune evasion mechanisms through genetic alterations. The IgH sequences identified here are potentially useful resources for future development of immunotherapies.