Single-cell RNA-Seq of follicular lymphoma reveals malignant B-cell types and coexpression of T-cell immune checkpoints.

Follicular lymphoma (FL) is a low-grade B-cell malignancy that transforms into a highly aggressive and lethal disease at a rate of 2% per year. Perfect isolation of the malignant B-cell population from a surgical biopsy is a significant challenge, masking important FL biology, such as immune checkpoint coexpression patterns. To resolve the underlying transcriptional networks of follicular B-cell lymphomas, we analyzed the transcriptomes of 34 188 cells derived from 6 primary FL tumors. For each tumor, we identified normal immune subpopulations and malignant B cells, based on gene expression. We used multicolor flow cytometry analysis of the same tumors to confirm our assignments of cellular lineages and validate our predictions of expressed proteins. Comparison of gene expression between matched malignant and normal B cells from the same patient revealed tumor-specific features. Malignant B cells exhibited restricted immunoglobulin (Ig) light chain expression (either Igκ or Igλ), as well the expected upregulation of the BCL2 gene, but also downregulation of the FCER2, CD52, and major histocompatibility complex class II genes. By analyzing thousands of individual cells per patient tumor, we identified the mosaic of malignant B-cell subclones that coexist within a FL and examined the characteristics of tumor-infiltrating T cells. We identified genes coexpressed with immune checkpoint molecules, such as CEBPA and B2M in regulatory T (Treg) cells, providing a better understanding of the gene networks involved in immune regulation. In summary, parallel measurement of single-cell expression in thousands of tumor cells and tumor-infiltrating lymphocytes can be used to obtain a systems-level view of the tumor microenvironment and identify new avenues for therapeutic development.

[1]  L. Staudt,et al.  Prediction of survival in follicular lymphoma based on molecular features of tumor-infiltrating immune cells. , 2004, The New England journal of medicine.

[2]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[3]  Grace X. Y. Zheng,et al.  Massively parallel digital transcriptional profiling of single cells , 2016, Nature Communications.

[4]  S. Linnarsson,et al.  Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing , 2014, Nature Neuroscience.

[5]  Drew M. Pardoll,et al.  The blockade of immune checkpoints in cancer immunotherapy , 2012, Nature Reviews Cancer.

[6]  Jenna M. Sullivan,et al.  TIM3+FOXP3+ regulatory T cells are tissue-specific promoters of T-cell dysfunction in cancer , 2013, Oncoimmunology.

[7]  Shawn M. Gillespie,et al.  Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma , 2014, Science.

[8]  V. Kuchroo,et al.  TIGIT predominantly regulates the immune response via regulatory T cells. , 2015, The Journal of clinical investigation.

[9]  Juan Carlos Fernández,et al.  Multiobjective evolutionary algorithms to identify highly autocorrelated areas: the case of spatial distribution in financially compromised farms , 2014, Ann. Oper. Res..

[10]  Sean C. Bendall,et al.  Data-Driven Phenotypic Dissection of AML Reveals Progenitor-like Cells that Correlate with Prognosis , 2015, Cell.

[11]  Sean C. Bendall,et al.  An interactive reference framework for modeling a dynamic immune system , 2015, Science.

[12]  Evan Z. Macosko,et al.  Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets , 2015, Cell.

[13]  Alice Giustacchini,et al.  Distinct myeloid progenitor differentiation pathways identified through single cell RNA sequencing , 2016, Nature Immunology.

[14]  P. Linsley,et al.  MAST: a flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data , 2015, Genome Biology.

[15]  Justin Guinney,et al.  GSVA: gene set variation analysis for microarray and RNA-Seq data , 2013, BMC Bioinformatics.

[16]  Ana C Anderson,et al.  Lag-3, Tim-3, and TIGIT: Co-inhibitory Receptors with Specialized Functions in Immune Regulation. , 2016, Immunity.

[17]  Raul Rabadan,et al.  Genetics of follicular lymphoma transformation. , 2014, Cell reports.

[18]  Paul Hoffman,et al.  Integrating single-cell transcriptomic data across different conditions, technologies, and species , 2018, Nature Biotechnology.

[19]  Aleksandra A. Kolodziejczyk,et al.  Classification of low quality cells from single-cell RNA-seq data , 2016, Genome Biology.

[20]  Ash A. Alizadeh,et al.  Hierarchy in somatic mutations arising during genomic evolution and progression of follicular lymphoma. , 2012, Blood.

[21]  Cole Trapnell,et al.  Defining cell types and states with single-cell genomics , 2015, Genome research.

[22]  Kevin R. Moon,et al.  Recovering Gene Interactions from Single-Cell Data Using Data Diffusion , 2018, Cell.

[23]  Tatyana Chernova,et al.  Th1-specific cell surface protein Tim-3 regulates macrophage activation and severity of an autoimmune disease , 2002, Nature.

[24]  Nci Dream Community A community effort to assess and improve drug sensitivity prediction algorithms , 2014 .

[25]  Sean C. Bendall,et al.  viSNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia , 2013, Nature Biotechnology.

[26]  M. Calaminici,et al.  Integrated genomic analysis identifies recurrent mutations and evolution patterns driving the initiation and progression of follicular lymphoma , 2013, Nature Genetics.

[27]  Chen Xu,et al.  Identification of cell types from single-cell transcriptomes using a novel clustering method , 2015, Bioinform..

[28]  Marko Bajec,et al.  Unfolding communities in large complex networks: Combining defensive and offensive label propagation for core extraction , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[29]  J. Friedberg,et al.  Transformed follicular non-Hodgkin lymphoma. , 2015, Blood.

[30]  Randy D Gascoyne,et al.  The architectural pattern of FOXP3-positive T cells in follicular lymphoma is an independent predictor of survival and histologic transformation. , 2010, Blood.

[31]  B. Nadel,et al.  Human germinal center transcriptional programs are de-synchronized in B cell lymphoma , 2018, Nature Immunology.

[32]  Victor X. Jin,et al.  Single-Cell RNA-seq Reveals a Subpopulation of Prostate Cancer Cells with Enhanced Cell-Cycle-Related Transcription and Attenuated Androgen Response. , 2018, Cancer research.

[33]  Ash A. Alizadeh,et al.  Mutations in early follicular lymphoma progenitors are associated with suppressed antigen presentation , 2015, Proceedings of the National Academy of Sciences.

[34]  Geoffrey E. Hinton,et al.  Visualizing Data using t-SNE , 2008 .