A computational tool suite to facilitate single-cell lineage tracing analyses.

Tracking the lineage relationships of cell populations is of increasing interest in diverse biological contexts. In this issue of Cell Reports Methods, Holze et al. present a suite of computational tools to facilitate such analyses and encourage their broader application.

Single-cell multiomics profiling reveals heterogeneous transcriptional programs and microenvironment in DSRCTs.

Desmoplastic small round cell tumor (DSRCT) is a rare, aggressive sarcoma driven by the EWSR1::WT1 chimeric transcription factor. Despite this unique oncogenic driver, DSRCT displays a polyphenotypic differentiation of unknown causality. Using single-cell multi-omics on 12 samples from five patients, we find that DSRCT tumor cells cluster into consistent subpopulations with partially overlapping lineage- and metabolism-related transcriptional programs. In vitro modeling shows that high EWSR1::WT1 DNA-binding activity associates with most lineage-related states, in contrast to glycolytic and profibrotic states. Single-cell chromatin accessibility analysis suggests that EWSR1::WT1 binding site variability may drive distinct lineage-related transcriptional programs, supporting some level of cell-intrinsic plasticity. Spatial transcriptomics reveals that glycolytic and profibrotic states specifically localize within hypoxic niches at the periphery of tumor cell islets, suggesting an additional role of tumor cell-extrinsic microenvironmental cues. We finally identify a single-cell transcriptomics-derived epithelial signature associated with improved patient survival, highlighting the clinical relevance of our findings.

Repeated peripheral infusions of anti-EGFRvIII CAR T cells in combination with pembrolizumab show no efficacy in glioblastoma: a phase 1 trial.

We previously showed that chimeric antigen receptor (CAR) T-cell therapy targeting epidermal growth factor receptor variant III (EGFRvIII) produces upregulation of programmed death-ligand 1 (PD-L1) in the tumor microenvironment (TME). Here we conducted a phase 1 trial (NCT03726515) of CAR T-EGFRvIII cells administered concomitantly with the anti-PD1 (aPD1) monoclonal antibody pembrolizumab in patients with newly diagnosed, EGFRvIII+ glioblastoma (GBM) (n = 7). The primary outcome was safety, and no dose-limiting toxicity was observed. Secondary outcomes included median progression-free survival (5.2 months; 90% confidence interval (CI), 2.9-6.0 months) and median overall survival (11.8 months; 90% CI, 9.2-14.2 months). In exploratory analyses, comparison of the TME in tumors harvested before versus after CAR + aPD1 administration demonstrated substantial evolution of the infiltrating myeloid and T cells, with more exhausted, regulatory, and interferon (IFN)-stimulated T cells at relapse. Our study suggests that the combination of CAR T cells and PD-1 inhibition in GBM is safe and biologically active but, given the lack of efficacy, also indicates a need to consider alternative strategies.

Imaging and multi-omics datasets converge to define different neural progenitor origins for ATRT-SHH subgroups.

Atypical teratoid rhabdoid tumors (ATRT) are divided into MYC, TYR and SHH subgroups, suggesting diverse lineages of origin. Here, we investigate the imaging of human ATRT at diagnosis and the precise anatomic origin of brain tumors in the Rosa26-CreERT2::Smarcb1flox/flox model. This cross-species analysis points to an extra-cerebral origin for MYC tumors. Additionally, we clearly distinguish SHH ATRT emerging from the cerebellar anterior lobe (CAL) from those emerging from the basal ganglia (BG) and intra-ventricular (IV) regions. Molecular characteristics point to the midbrain-hindbrain boundary as the origin of CAL SHH ATRT, and to the ganglionic eminence as the origin of BG/IV SHH ATRT. Single-cell RNA sequencing on SHH ATRT supports these hypotheses. Trajectory analyses suggest that SMARCB1 loss induces a de-differentiation process mediated by repressors of the neuronal program such as REST, ID and the NOTCH pathway.

Multi-omics comparison of malignant and normal uveal melanocytes reveals molecular features of uveal melanoma.

Uveal melanoma (UM) is a rare cancer resulting from the transformation of melanocytes in the uveal tract. Integrative analysis has identified four molecular and clinical subsets of UM. To improve our molecular understanding of UM, we performed extensive multi-omics characterization comparing two aggressive UM patient-derived xenograft models with normal choroidal melanocytes, including DNA optical mapping, specific histone modifications, and DNA topology analysis using Hi-C. Our gene expression and cytogenetic analyses suggest that genomic instability is a hallmark of UM. We also identified a recurrent deletion in the BAP1 promoter resulting in loss of expression and associated with high risk of metastases in UM patients. Hi-C revealed chromatin topology changes associated with the upregulation of PRAME, an independent prognostic biomarker in UM, and a potential therapeutic target. Our findings illustrate how multi-omics approaches can improve our understanding of tumorigenesis and reveal two distinct mechanisms of gene expression dysregulation in UM.

Cell-Free DNA Extracted from CSF for the Molecular Diagnosis of Pediatric Embryonal Brain Tumors.

BACKGROUND: Liquid biopsies are revolutionary tools used to detect tumor-specific genetic alterations in body fluids, including the use of cell-free DNA (cfDNA) for molecular diagnosis in cancer patients. In brain tumors, cerebrospinal fluid (CSF) cfDNA might be more informative than plasma cfDNA. Here, we assess the use of CSF cfDNA in pediatric embryonal brain tumors (EBT) for molecular diagnosis. METHODS: The CSF cfDNA of pediatric patients with medulloblastoma (n = 18), ATRT (n = 3), ETMR (n = 1), CNS NB FOXR2 (n = 2) and pediatric EBT NOS (n = 1) (mean cfDNA concentration 48 ng/mL; range 4-442 ng/mL) and matched tumor genomic DNA were sequenced by WES and/or a targeted sequencing approach to determine single-nucleotide variations (SNVs) and copy number alterations (CNA). A specific capture covering transcription start sites (TSS) of genes of interest was also used for nucleosome footprinting in CSF cfDNA. RESULTS: 15/25 CSF cfDNA samples yielded informative results, with informative CNA and SNVs in 11 and 15 cases, respectively. For cases with paired tumor and CSF cfDNA WES (n = 15), a mean of 83 (range 1-160) shared SNVs were observed, including SNVs in classical medulloblastoma genes such as SMO and KMT2D. Interestingly, tumor-specific SNVs (mean 18; range 1-62) or CSF-specific SNVs (mean 5; range 0-25) were also observed, suggesting clonal heterogeneity. The TSS panel resulted in differential coverage profiles across all 112 studied genes in 7 cases, indicating distinct promoter accessibility. CONCLUSION: CSF cfDNA sequencing yielded informative results in 60% (15/25) of all cases, with informative results in 83% (15/18) of all cases analyzed by WES. These results pave the way for the implementation of these novel approaches for molecular diagnosis and minimal residual disease monitoring.

Evolution of synchronous female bilateral breast cancers and response to treatment.

Synchronous bilateral breast cancer (sBBC) occurs after both breasts have been affected by the same germline genetics and environmental exposures. Little evidence exists regarding immune infiltration and response to treatment in sBBCs. Here we show that the impact of the subtype of breast cancer on levels of tumor infiltrating lymphocytes (TILs, n = 277) and on pathologic complete response (pCR) rates (n = 140) differed according to the concordant or discordant subtype of breast cancer of the contralateral tumor: luminal breast tumors with a discordant contralateral tumor had higher TIL levels and higher pCR rates than those with a concordant contralateral tumor. Tumor sequencing revealed that left and right tumors (n = 20) were independent regarding somatic mutations, copy number alterations and clonal phylogeny, whereas primary tumor and residual disease were closely related both from the somatic mutation and from the transcriptomic point of view. Our study indicates that tumor-intrinsic characteristics may have a role in the association of tumor immunity and pCR and demonstrates that the characteristics of the contralateral tumor are also associated with immune infiltration and response to treatment.

Noncanonical splicing junctions between exons and transposable elements represent a source of immunogenic recurrent neo-antigens in patients with lung cancer.

Although most characterized tumor antigens are encoded by canonical transcripts (such as differentiation or tumor-testis antigens) or mutations (both driver and passenger mutations), recent results have shown that noncanonical transcripts including long noncoding RNAs and transposable elements (TEs) can also encode tumor-specific neo-antigens. Here, we investigate the presentation and immunogenicity of tumor antigens derived from noncanonical mRNA splicing events between coding exons and TEs. Comparing human non-small cell lung cancer (NSCLC) and diverse healthy tissues, we identified a subset of splicing junctions that is both tumor specific and shared across patients. We used HLA-I peptidomics to identify peptides encoded by tumor-specific junctions in primary NSCLC samples and lung tumor cell lines. Recurrent junction-encoded peptides were immunogenic in vitro, and CD8+ T cells specific for junction-encoded epitopes were present in tumors and tumor-draining lymph nodes from patients with NSCLC. We conclude that noncanonical splicing junctions between exons and TEs represent a source of recurrent, immunogenic tumor-specific antigens in patients with NSCLC.

Epigenetically controlled tumor antigens derived from splice junctions between exons and transposable elements.

Oncogenesis often implicates epigenetic alterations, including derepression of transposable elements (TEs) and defects in alternative splicing. Here, we explore the possibility that noncanonical splice junctions between exons and TEs represent a source of tumor-specific antigens. We show that mouse normal tissues and tumor cell lines express wide but distinct ranges of mRNA junctions between exons and TEs, some of which are tumor specific. Immunopeptidome analyses in tumor cell lines identified peptides derived from exon-TE splicing junctions associated to MHC-I molecules. Exon-TE junction-derived peptides were immunogenic in tumor-bearing mice. Both prophylactic and therapeutic vaccinations with junction-derived peptides delayed tumor growth in vivo. Inactivation of the TE-silencing histone 3-lysine 9 methyltransferase Setdb1 caused overexpression of new immunogenic junctions in tumor cells. Our results identify exon-TE splicing junctions as epigenetically controlled, immunogenic, and protective tumor antigens in mice, opening possibilities for tumor targeting and vaccination in patients with cancer.

CD8+T cell responsiveness to anti-PD-1 is epigenetically regulated by Suv39h1 in melanomas.

Tumor-infiltrating CD8 + T cells progressively lose functionality and fail to reject tumors. The underlying mechanism and re-programing induced by checkpoint blockers are incompletely understood. We show here that genetic ablation or pharmacological inhibition of histone lysine methyltransferase Suv39h1 delays tumor growth and potentiates tumor rejection by anti-PD-1. In the absence of Suv39h1, anti-PD-1 induces alternative activation pathways allowing survival and differentiation of IFNγ and Granzyme B producing effector cells that express negative checkpoint molecules, but do not reach final exhaustion. Their transcriptional program correlates with that of melanoma patients responding to immune-checkpoint blockade and identifies the emergence of cytolytic-effector tumor-infiltrating lymphocytes as a biomarker of clinical response. Anti-PD-1 favors chromatin opening in loci linked to T-cell activation, memory and pluripotency, but in the absence of Suv39h1, cells acquire accessibility in cytolytic effector loci. Overall, Suv39h1 inhibition enhances anti-tumor immune responses, alone or combined with anti-PD-1, suggesting that Suv39h1 is an "epigenetic checkpoint" for tumor immunity.

Oncogenic chimeric transcription factors drive tumor-specific transcription, processing, and translation of silent genomic regions.

Many cancers are characterized by gene fusions encoding oncogenic chimeric transcription factors (TFs) such as EWS::FLI1 in Ewing sarcoma (EwS). Here, we find that EWS::FLI1 induces the robust expression of a specific set of novel spliced and polyadenylated transcripts within otherwise transcriptionally silent regions of the genome. These neogenes (NGs) are virtually undetectable in large collections of normal tissues or non-EwS tumors and can be silenced by CRISPR interference at regulatory EWS::FLI1-bound microsatellites. Ribosome profiling and proteomics further show that some NGs are translated into highly EwS-specific peptides. More generally, we show that hundreds of NGs can be detected in diverse cancers characterized by chimeric TFs. Altogether, this study identifies the transcription, processing, and translation of novel, specific, highly expressed multi-exonic transcripts from otherwise silent regions of the genome as a new activity of aberrant TFs in cancer.

Identification of Tissue of Origin and Guided Therapeutic Applications in Cancers of Unknown Primary Using Deep Learning and RNA Sequencing (TransCUPtomics).

Cancers of unknown primary (CUP) are metastatic cancers for which the primary tumor is not found despite thorough diagnostic investigations. Multiple molecular assays have been proposed to identify the tissue of origin (TOO) and inform clinical care; however, none has been able to combine accuracy, interpretability, and easy access for routine use. We developed a classifier tool based on the training of a variational autoencoder to predict tissue of origin based on RNA-sequencing data. We used as training data 20,918 samples corresponding to 94 different categories, including 39 cancer types and 55 normal tissues. The TransCUPtomics classifier was applied to a retrospective cohort of 37 CUP patients and 11 prospective patients. TransCUPtomics exhibited an overall accuracy of 96% on reference data for TOO prediction. The TOO could be identified in 38 (79%) of 48 CUP patients. Eight of 11 prospective CUP patients (73%) could receive first-line therapy guided by TransCUPtomics prediction, with responses observed in most patients. The variational autoencoder added further utility by enabling prediction interpretability, and diagnostic predictions could be matched to detection of gene fusions and expressed variants. TransCUPtomics confidently predicted TOO for CUP and enabled tailored treatments leading to significant clinical responses. The interpretability of our approach is a powerful addition to improve the management of CUP patients.

Splicing Patterns in SF3B1-Mutated Uveal Melanoma Generate Shared Immunogenic Tumor-Specific Neoepitopes.

Disruption of splicing patterns due to mutations of genes coding splicing factors in tumors represents a potential source of tumor neoantigens, which would be both public (shared between patients) and tumor-specific (not expressed in normal tissues). In this study, we show that mutations of the splicing factor SF3B1 in uveal melanoma generate such immunogenic neoantigens. Memory CD8+ T cells specific for these neoantigens are preferentially found in 20% of patients with uveal melanoma bearing SF3B1-mutated tumors. Single-cell analyses of neoepitope-specific T cells from the blood identified large clonal T-cell expansions, with distinct effector transcription patterns. Some of these expanded T-cell receptors are also present in the corresponding tumors. CD8+ T-cell clones specific for the neoepitopes specifically recognize and kill SF3B1-mutated tumor cells, supporting the use of this new family of neoantigens as therapeutic targets. SIGNIFICANCE: Mutations of the splicing factor SF3B1 in uveal melanoma generate shared neoantigens that are uniquely expressed by tumor cells, leading to recognition and killing by specific CD8 T cells. Mutations in splicing factors can be sources of new therapeutic strategies applicable to diverse tumors.This article is highlighted in the In This Issue feature, p. 1861.

Contribution of resident and circulating precursors to tumor-infiltrating CD8+ T cell populations in lung cancer.

Tumor-infiltrating lymphocytes (TILs), in general, and especially CD8+ TILs, represent a favorable prognostic factor in non-small cell lung cancer (NSCLC). The tissue origin, regenerative capacities, and differentiation pathways of TIL subpopulations remain poorly understood. Using a combination of single-cell RNA and T cell receptor (TCR) sequencing, we investigate the functional organization of TIL populations in primary NSCLC. We identify two CD8+ TIL subpopulations expressing memory-like gene modules: one is also present in blood (circulating precursors) and the other one in juxtatumor tissue (tissue-resident precursors). In tumors, these two precursor populations converge through a unique transitional state into terminally differentiated cells, often referred to as dysfunctional or exhausted. Differentiation is associated with TCR expansion, and transition from precursor to late-differentiated states correlates with intratumor T cell cycling. These results provide a coherent working model for TIL origin, ontogeny, and functional organization in primary NSCLC.

Genetic alterations of SUGP1 mimic mutant-SF3B1 splice pattern in lung adenocarcinoma and other cancers.

Genes involved in 3'-splice site recognition during mRNA splicing constitute an emerging class of oncogenes. SF3B1 is the most frequently mutated splicing factor in cancer, and SF3B1 mutants corrupt branchpoint recognition leading to usage of cryptic 3'-splice sites and subsequent aberrant junctions. For a comprehensive determination of alterations leading to this splicing pattern, we performed a pan-TCGA screening for SF3B1-specific aberrant acceptor usage. While the most of aberrant 3'-splice patterns were explained by SF3B1 mutations, we also detected nine SF3B1 wild-type tumors (including five lung adenocarcinomas). Genomic profile analysis of these tumors identified somatic mutations combined with loss-of-heterozygosity in the splicing factor SUGP1 in five of these cases. Modeling of SUGP1 loss and mutations in cell lines showed that both alterations induced mutant-SF3B1-like aberrant splicing. Our study provides definitive evidence that genetic alterations of SUGP1 genocopy SF3B1 mutations in lung adenocarcinoma and other cancers.

Author Correction: Recurrent PTPRT/JAK2 mutations in lung adenocarcinoma among African Americans.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Recurrent PTPRT/JAK2 mutations in lung adenocarcinoma among African Americans.

Reducing or eliminating persistent disparities in lung cancer incidence and survival has been challenging because our current understanding of lung cancer biology is derived primarily from populations of European descent. Here we show results from a targeted sequencing panel using NCI-MD Case Control Study patient samples and reveal a significantly higher prevalence of PTPRT and JAK2 mutations in lung adenocarcinomas among African Americans compared with European Americans. This increase in mutation frequency was validated with independent WES data from the NCI-MD Case Control Study and TCGA. We find that patients carrying these mutations have a concomitant increase in IL-6/STAT3 signaling and miR-21 expression. Together, these findings suggest the identification of these potentially actionable mutations could have clinical significance for targeted therapy and the enrollment of minority populations in clinical trials.

1q21.1 deletion and a rare functional polymorphism in siblings with thrombocytopenia-absent radius-like phenotypes.

Thrombocytopenia-absent radii (TAR) syndrome, characterized by neonatal thrombocytopenia and bilateral radial aplasia with thumbs present, is typically caused by the inheritance of a 1q21.1 deletion and a single-nucelotide polymorphism in RBM8A on the nondeleted allele. We evaluated two siblings with TAR-like dysmorphology but lacking thrombocytopenia in infancy. Family NCI-107 participated in an IRB-approved cohort study and underwent comprehensive clinical and genomic evaluations, including aCGH, whole-exome, whole-genome, and targeted sequencing. Gene expression assays and electromobility shift assays (EMSAs) were performed to evaluate the variant of interest. The previously identified TAR-associated 1q21.1 deletion was present in the affected siblings and one healthy parent. Multiple sequencing approaches did not identify previously described TAR-associated SNPs or mutations in relevant genes. We discovered rs61746197 A > G heterozygosity in the parent without the deletion and apparent hemizygosity in both siblings. rs61746197 A > G overlaps a RelA-p65 binding motif, and EMSAs indicate the A allele has higher transcription factor binding efficiency than the G allele. Stimulation of K562 cells to induce megakaryocyte differentiation abrogated the shift of both reference and alternative probes. The 1q21.1 TAR-associated deletion in combination with the G variant of rs61746197 on the nondeleted allele is associated with a TAR-like phenotype. rs61746197 G could be a functional enhancer/repressor element, but more studies are required to identify the specific factor(s) responsible. Overall, our findings suggest a role of rs61746197 A > G and human disease in the setting of a 1q21.1 deletion on the other chromosome.