Intratumor heterogeneity and T cell exhaustion in primary CNS lymphoma.

BACKGROUND: Primary central nervous system lymphoma (PCNSL) is a rare lymphoma of the central nervous system, usually of diffuse large B cell phenotype. Stereotactic biopsy followed by histopathology is the diagnostic standard. However, limited material is available from CNS biopsies, thus impeding an in-depth characterization of PCNSL. METHODS: We performed flow cytometry, single-cell RNA sequencing, and B cell receptor sequencing of PCNSL cells released from biopsy material, blood, and cerebrospinal fluid (CSF), and spatial transcriptomics of biopsy samples. RESULTS: PCNSL-released cells were predominantly activated CD19+CD20+CD38+CD27+ B cells. In single-cell RNA sequencing, PCNSL cells were transcriptionally heterogeneous, forming multiple malignant B cell clusters. Hyperexpanded B cell clones were shared between biopsy- and CSF- but not blood-derived cells. T cells in the tumor microenvironment upregulated immune checkpoint molecules, thereby recognizing immune evasion signals from PCNSL cells. Spatial transcriptomics revealed heterogeneous spatial organization of malignant B cell clusters, mirroring their transcriptional heterogeneity across patients, and pronounced expression of T cell exhaustion markers, co-localizing with a highly malignant B cell cluster. CONCLUSIONS: Malignant B cells in PCNSL show transcriptional and spatial intratumor heterogeneity. T cell exhaustion is frequent in the PCNSL microenvironment, co-localizes with malignant cells, and highlights the potential of personalized treatments.

Multiplexed single-cell analysis reveals prognostic and nonprognostic T cell types in human colorectal cancer.

Clinical outcomes in colorectal cancer (CRC) correlate with T cell infiltrates, but the specific contributions of heterogenous T cell types remain unclear. To investigate the diverse function of T cells in CRC, we profiled 37,931 T cells from tumors and adjacent normal colon of 16 patients with CRC with respect to transcriptome, TCR sequence, and cell surface markers. Our analysis identified phenotypically and functionally distinguishable effector T cell types. We employed single-cell gene signatures from these T cell subsets to query the TCGA database to assess their prognostic significance. We found 2 distinct cytotoxic T cell types. GZMK+KLRG1+ cytotoxic T cells were enriched in CRC patients with good outcomes. GNLY+CD103+ cytotoxic T cells with a dysfunctional phenotype were not associated with good outcomes, despite coexpression of CD39 and CD103, markers that denote tumor reactivity. We found 2 distinct Treg subtypes associated with opposite outcomes. While total Tregs were associated with good outcomes, CD38+ Tregs were associated with bad outcomes independently of stage and possessed a highly suppressive phenotype, suggesting that they inhibit antitumor immunity in CRC. These findings highlight the potential utility of these subpopulations in predicting outcomes and support the potential for novel therapies directed at CD38+ Tregs or CD8+CD103+ T cells.

Clonal lineage tracing reveals shared origin of conventional and plasmacytoid dendritic cells.

Developmental origins of dendritic cells (DCs) including conventional DCs (cDCs, comprising cDC1 and cDC2 subsets) and plasmacytoid DCs (pDCs) remain unclear. We studied DC development in unmanipulated adult mice using inducible lineage tracing combined with clonal DNA "barcoding" and single-cell transcriptome and phenotype analysis (CITE-seq). Inducible tracing of Cx3cr1+ hematopoietic progenitors in the bone marrow showed that they simultaneously produce all DC subsets including pDCs, cDC1s, and cDC2s. Clonal tracing of hematopoietic stem cells (HSCs) and of Cx3cr1+ progenitors revealed clone sharing between cDC1s and pDCs, but not between the two cDC subsets or between pDCs and B cells. Accordingly, CITE-seq analyses of differentiating HSCs and Cx3cr1+ progenitors identified progressive stages of pDC development including Cx3cr1+ Ly-6D+ pro-pDCs that were distinct from lymphoid progenitors. These results reveal the shared origin of pDCs and cDCs and suggest a revised scheme of DC development whereby pDCs share clonal relationship with cDC1s.

Tracing tumorigenesis in a solid tumor model at single-cell resolution.

Characterizing the complex composition of solid tumors is fundamental for understanding tumor initiation, progression and metastasis. While patient-derived samples provide valuable insight, they are heterogeneous on multiple molecular levels, and often originate from advanced tumor stages. Here, we use single-cell transcriptome and epitope profiling together with pathway and lineage analyses to study tumorigenesis from a developmental perspective in a mouse model of salivary gland squamous cell carcinoma. We provide a comprehensive cell atlas and characterize tumor-specific cells. We find that these cells are connected along a reproducible developmental trajectory: initiated in basal cells exhibiting an epithelial-to-mesenchymal transition signature, tumorigenesis proceeds through Wnt-differential cancer stem cell-like subpopulations before differentiating into luminal-like cells. Our work provides unbiased insights into tumor-specific cellular identities in a whole tissue environment, and emphasizes the power of using defined genetic model systems.

Multiplexed detection of proteins, transcriptomes, clonotypes and CRISPR perturbations in single cells.

Multimodal single-cell assays provide high-resolution snapshots of complex cell populations, but are mostly limited to transcriptome plus an additional modality. Here, we describe expanded CRISPR-compatible cellular indexing of transcriptomes and epitopes by sequencing (ECCITE-seq) for the high-throughput characterization of at least five modalities of information from each single cell. We demonstrate application of ECCITE-seq to multimodal CRISPR screens with robust direct single-guide RNA capture and to clonotype-aware multimodal phenotyping of cancer samples.

Single-cell stabilization method identifies gonadotrope transcriptional dynamics and pituitary cell type heterogeneity.

Immediate-early response genes (IEGs) are rapidly and transiently induced following an extracellular signal. Elucidating the IEG response patterns in single cells (SCs) requires assaying large numbers of timed samples at high accuracy while minimizing handling effects. To achieve this, we developed and validated RNA stabilization Buffer for Examination of Single-cell Transcriptomes (RNA-Best), a versatile single-step cell and tissue preservation protocol that stabilizes RNA in intact SCs without perturbing transcription patterns. We characterize for the first time SC heterogeneity in IEG responses to pulsatile gonadotropin-releasing hormone (GnRH) stimuli in pituitary gonadotrope cells. Our study identifies a gene-specific hierarchical pattern of all-or-none transcript induction elicited by increasing concentrations of GnRH. This quantal pattern of gene activation raises the possibility that IEG activation, when accurately resolved at the SC level, may be mediated by gene bits that behave as pure binary switches.

Essential roles of Raf/extracellular signal-regulated kinase/mitogen-activated protein kinase pathway, YY1, and Ca2+ influx in growth arrest of human vascular smooth muscle cells by bilirubin.

The biological effects of bilirubin, still poorly understood, are concentration-dependent ranging from cell protection to toxicity. Here we present data that at high nontoxic physiological concentrations, bilirubin inhibits growth of proliferating human coronary artery smooth muscle cells by three events. It impairs the activation of Raf/ERK/MAPK pathway and the cellular Raf and cyclin D1 content that results in retinoblastoma protein hypophosphorylation on amino acids S608 and S780. These events impede the release of YY1 to the nuclei and its availability to regulate the expression of genes and to support cellular proliferation. Moreover, altered calcium influx and calpain II protease activation leads to proteolytical degradation of transcription factor YY1. We conclude that in the serum-stimulated human vascular smooth muscle primary cell cultures, bilirubin favors growth arrest, and we propose that this activity is regulated by its interaction with the Raf/ERK/MAPK pathway, effect on cyclin D1 and Raf content, altered retinoblastoma protein profile of hypophosphorylation, calcium influx, and YY1 proteolysis. We propose that these activities together culminate in diminished 5 S and 45 S ribosomal RNA synthesis and cell growth arrest. The observations provide important mechanistic insight into the molecular mechanisms underlying the transition of human vascular smooth muscle cells from proliferative to contractile phenotype and the role of bilirubin in this transition.