A novel memory-like Tfh cell subset is precursor to effector Tfh cells in recall immune responses.

T follicular helper (Tfh) cells, essential for germinal center reactions, are not identical, with different phenotypes reported. Whether, when, and how they generate memory cells is still poorly understood. Here, through single-cell RNA-sequencing analysis of CXCR5+Bcl6+ Tfh cells generated under different conditions, we discovered, in addition to PD-1hi effector Tfh cells, a CD62L+PD1low subpopulation. CD62L-expressing Tfh cells developed independently from PD-1+ cells and not in direct contact with B cells. More importantly, CD62L+ Tfh cells expressed memory- and stemness-associated genes, and with better superior long-term survival, they readily generated PD-1hi cells in the recall response. Finally, KLF2 and IL7R, also highly expressed by CD62L+ Tfh cells, were required to regulate their development. Our work thus demonstrates a novel Tfh memory-like cell subpopulation, which may benefit our understanding of immune responses and diseases.

Single-cell RNA sequencing highlights the role of PVR/PVRL2 in the immunosuppressive tumour microenvironment in hepatocellular carcinoma.

Introduction: The conflict between cancer cells and the host immune system shapes the immune tumour microenvironment (TME) in hepatocellular carcinoma (HCC). A deep understanding of the heterogeneity and intercellular communication network in the TME of HCC will provide promising strategies to orchestrate the immune system to target and eradicate cancers. Methods: Here, we performed single-cell RNA sequencing (scRNA-seq) and computational analysis of 35786 unselected single cells from 3 human HCC tumour and 3 matched adjacent samples to elucidate the heterogeneity and intercellular communication network of the TME. The specific lysis of HCC cell lines was examined in vitro using cytotoxicity assays. Granzyme B concentration in supernatants of cytotoxicity assays was measured by ELISA. Results: We found that VCAN+ tumour-associated macrophages (TAMs) might undergo M2-like polarization and differentiate in the tumour region. Regulatory dendritic cells (DCs) exhibited immune regulatory and tolerogenic phenotypes in the TME. Furthermore, we observed intensive potential intercellular crosstalk among C1QC+ TAMs, regulatory DCs, regulator T (Treg) cells, and exhausted CD8+ T cells that fostered an immunosuppressive niche in the HCC TME. Moreover, we identified that the TIGIT-PVR/PVRL2 axis provides a prominent coinhibitory signal in the immunosuppressive TME. In vitro, antibody blockade of PVR or PVRL2 on HCC cell lines or TIGIT blockade on immune cells increased immune cell-mediated lysis of tumour cell. This enhanced immune response is paralleled by the increased secretion of Granzyme B by immune cells. Discussion: Collectively, our study revealed the functional state, clinical significance, and intercellular communication of immunosuppressive cells in HCC at single-cell resolution. Moreover, PVR/PVRL2, interact with TIGIT act as prominent coinhibitory signals and might represent a promising, efficacious immunotherapy strategy in HCC.

Parallel single-cell and bulk transcriptome analyses reveal key features of the gastric tumor microenvironment.

BACKGROUND: The tumor microenvironment (TME) has been shown to strongly influence treatment outcome for cancer patients in various indications and to influence the overall survival. However, the cells forming the TME in gastric cancer have not been extensively characterized. RESULTS: We combine bulk and single-cell RNA sequencing from tumors and matched normal tissue of 24 treatment-naïve GC patients to better understand which cell types and transcriptional programs are associated with malignant transformation of the stomach. Clustering 96,623 cells of non-epithelial origin reveals 81 well-defined TME cell types. We find that activated fibroblasts and endothelial cells are most prominently overrepresented in tumors. Intercellular network reconstruction and survival analysis of an independent cohort imply the importance of these cell types together with immunosuppressive myeloid cell subsets and regulatory T cells in establishing an immunosuppressive microenvironment that correlates with worsened prognosis and lack of response in anti-PD1-treated patients. In contrast, we find a subset of IFNγ activated T cells and HLA-II expressing macrophages that are linked to treatment response and increased overall survival. CONCLUSIONS: Our gastric cancer single-cell TME compendium together with the matched bulk transcriptome data provides a unique resource for the identification of new potential biomarkers for patient stratification. This study helps further to elucidate the mechanism of gastric cancer and provides insights for therapy.

Cellular basis of enhanced humoral immunity to SARS-CoV-2 upon homologous or heterologous booster vaccination analyzed by single-cell immune profiling.

SARS-CoV-2 vaccine booster dose can induce a robust humoral immune response, however, its cellular mechanisms remain elusive. Here, we investigated the durability of antibody responses and single-cell immune profiles following booster dose immunization, longitudinally over 6 months, in recipients of a homologous BBIBP-CorV/BBIBP-CorV or a heterologous BBIBP-CorV/ZF2001 regimen. The production of neutralizing antibodies was dramatically enhanced by both booster regimens, and the antibodies could last at least six months. The heterologous booster induced a faster and more robust plasmablast response, characterized by activation of plasma cells than the homologous booster. The response was attributed to recall of memory B cells and the de novo activation of B cells. Expanded B cell clones upon booster dose vaccination could persist for months, and their B cell receptors displayed accumulated mutations. The production of antibody was positively correlated with antigen presentation by conventional dendritic cells (cDCs), which provides support for B cell maturation through activation and development of follicular helper T (Tfh) cells. The proper activation of cDC/Tfh/B cells was likely fueled by active energy metabolism, and glutaminolysis might also play a general role in promoting humoral immunity. Our study unveils the cellular mechanisms of booster-induced memory/adaptive humoral immunity and suggests potential strategies to optimize vaccine efficacy and durability in future iterations.

Single-cell immune profiling reveals functional diversity of T cells in tuberculous pleural effusion.

Orchestration of an effective T lymphocyte response at infection sites is critical for protection against Mycobacterium tuberculosis (Mtb) infection. However, the local T cell immunity landscape in human tuberculosis is poorly defined. Tuberculous pleural effusion (TPE), caused by Mtb, is characterized by an influx of leukocytes to the pleural space, providing a platform suitable for delineating complex tissue responses to Mtb infection. Using single-cell transcriptomics and T cell receptor sequencing, we analyzed mononuclear cell populations in paired pleural fluid and peripheral blood of TPE patients. While all major cell clusters were present in both tissues, their relative proportions varied significantly by anatomic location. Lineage tracking analysis revealed subsets of CD8 and CD4 T cell populations with distinct effector functions specifically expanded at pleural sites. Granzyme K-expressing CD8 T cells were preferentially enriched and clonally expanded in pleural fluid from TPE, suggesting that they are involved in the pathogenesis of the disease. The findings collectively reveal the landscape of local T cell immunity in tuberculosis.

STARTRAC analyses of scRNAseq data from tumor models reveal T cell dynamics and therapeutic targets.

Single-cell RNA sequencing is a powerful tool to examine cellular heterogeneity, novel markers and target genes, and therapeutic mechanisms in human cancers and animal models. Here, we analyzed single-cell RNA sequencing data of T cells obtained from multiple mouse tumor models by PCA-based subclustering coupled with TCR tracking using the STARTRAC algorithm. This approach revealed various differentiated T cell subsets and activation states, and a correspondence of T cell subsets between human and mouse tumors. STARTRAC analyses demonstrated peripheral T cell subsets that were developmentally connected with tumor-infiltrating CD8+ cells, CD4+ Th1 cells, and T reg cells. In addition, large amounts of paired TCRα/ß sequences enabled us to identify a specific enrichment of paired public TCR clones in tumor. Finally, we identified CCR8 as a tumor-associated T reg cell marker that could preferentially deplete tumor-associated T reg cells. We showed that CCR8-depleting antibody treatment provided therapeutic benefit in CT26 tumors and synergized with anti-PD-1 treatment in MC38 and B16F10 tumor models.

SciBet as a portable and fast single cell type identifier.

Fast, robust and technology-independent computational methods are needed for supervised cell type annotation of single-cell RNA sequencing data. We present SciBet, a supervised cell type identifier that accurately predicts cell identity for newly sequenced cells with order-of-magnitude speed advantage. We enable web client deployment of SciBet for rapid local computation without uploading local data to the server. Facing the exponential growth in the size of single cell RNA datasets, this user-friendly and cross-platform tool can be widely useful for single cell type identification.

Stress-glucocorticoid-TSC22D3 axis compromises therapy-induced antitumor immunity.

Psychological distress has long been suspected to influence cancer incidence and mortality. It remains largely unknown whether and how stress affects the efficacy of anticancer therapies. We observed that social defeat caused anxiety-like behaviors in mice and dampened therapeutic responses against carcinogen-induced neoplasias and transplantable tumors. Stress elevated plasma corticosterone and upregulated the expression of glucocorticoid-inducible factor Tsc22d3, which blocked type I interferon (IFN) responses in dendritic cell (DC) and IFN-γ+ T cell activation. Similarly, close correlations were discovered among plasma cortisol levels, TSC22D3 expression in circulating leukocytes and negative mood in patients with cancer. In murine models, exogenous glucocorticoid injection, or enforced expression of Tsc22d3 in DC was sufficient to abolish therapeutic control of tumors. Administration of a glucocorticoid receptor antagonist or DC-specific Tsc22d3 deletion reversed the negative impact of stress or glucocorticoid supplementation on therapeutic outcomes. Altogether, these results indicate that stress-induced glucocorticoid surge and Tsc22d3 upregulation can subvert therapy-induced anticancer immunosurveillance.

Cancer IgG, a potential prognostic marker, promotes colorectal cancer progression.

OBJECTIVE: Currently, no satisfactory targets for colorectal cancer or markers for immunotherapy and diagnosis and prognosis are available. Immunoglobulin G (IgG) is widely expressed in many cancers, and it promotes cancer progression. This study explored the role of cancer-derived IgG (CIgG) in colorectal cancer. METHODS: First, using a monoclonal antibody to CIgG, we examined the expression levels of CIgG in colorectal cancer cell lines by western blot and immunofluorescence analyses and in tissue specimens by immunohistochemistry. Second, the variable region gene was amplified by nested polymerase chain reaction (PCR), and PCR products were sequenced and analyzed. Third, we investigated the effect of CIgG on colorectal cancer cells by cell proliferation, wound healing, migration and invasion assays, and colony formation assay. Fourth, we performed in vivo tumorigenicity experiments to explore the effect of CIgG on tumorigenicity. Finally, we used RNA-seq analysis and co-immunoprecipitation experiments to further clarify possible mechanisms of CIgG. RESULTS: We found that CIgG is widely expressed in colorectal cancer cells, and the overexpression of CIgG indicates significantly poor colorectal cancer prognosis. Furthermore, CIgG knockdown significantly inhibits the proliferation, migration and invasion ability of cells, and tumor growth in vivo. RNA-seq analysis indicated that CIgG knockdown results primarily in changes in expression of apical junction and epithelial-mesenchymal transition-related genes. CIgG may be involved in colorectal cancer invasion and metastasis through interacting with E-cadherin. CONCLUSIONS: CIgG is a potential human oncogene in colorectal cancer and that it has potential for application as a novel target in targeted therapy and a marker for prognostic evaluation.

GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis.

Introduced in 2017, the GEPIA (Gene Expression Profiling Interactive Analysis) web server has been a valuable and highly cited resource for gene expression analysis based on tumor and normal samples from the TCGA and the GTEx databases. Here, we present GEPIA2, an updated and enhanced version to provide insights with higher resolution and more functionalities. Featuring 198 619 isoforms and 84 cancer subtypes, GEPIA2 has extended gene expression quantification from the gene level to the transcript level, and supports analysis of a specific cancer subtype, and comparison between subtypes. In addition, GEPIA2 has adopted new analysis techniques of gene signature quantification inspired by single-cell sequencing studies, and provides customized analysis where users can upload their own RNA-seq data and compare them with TCGA and GTEx samples. We also offer an API for batch process and easy retrieval of the analysis results. The updated web server is publicly accessible at http://gepia2.cancer-pku.cn/.

CDK5RAP3, a UFL1 substrate adaptor, is crucial for liver development.

Protein modification by ubiquitin and ubiquitin-like proteins (UBLs) regulates numerous biological functions. The UFM1 system, a novel UBL conjugation system, is implicated in mouse development and hematopoiesis. However, its broad biological functions and working mechanisms remain largely elusive. CDK5RAP3, a possible ufmylation substrate, is essential for epiboly and gastrulation in zebrafish. Herein, we report a crucial role of CDK5RAP3 in liver development and hepatic functions. Cdk5rap3 knockout mice displayed prenatal lethality with severe liver hypoplasia, as characterized by delayed proliferation and compromised differentiation. Hepatocyte-specific Cdk5rap3 knockout mice suffered post-weaning lethality, owing to serious hypoglycemia and impaired lipid metabolism. Depletion of CDK5RAP3 triggered endoplasmic reticulum stress and activated unfolded protein responses in hepatocytes. We detected the in vivo interaction of CDK5RAP3 with UFL1, the defined E3 ligase in ufmylation. Notably, loss of CDK5RAP3 altered the ufmylation profile in liver cells, suggesting that CDK5RAP3 serves as a novel substrate adaptor for this UBL modification. Collectively, our study identifies CDK5RAP3 as an important regulator of ufmylation and suggests the involvement of ufmylation in mammalian development.

Distinct epigenetic features of tumor-reactive CD8+ T cells in colorectal cancer patients revealed by genome-wide DNA methylation analysis.

BACKGROUND: Tumor-reactive CD8+ tumor-infiltrating lymphocytes (TILs) represent a subtype of T cells that can recognize and destroy tumor specifically. Understanding the regulatory mechanism of tumor-reactive CD8+ T cells has important therapeutic implications. Yet the DNA methylation status of this T cell subtype has not been elucidated. RESULTS: In this study, we segregate tumor-reactive and bystander CD8+ TILs, as well as naïve and effector memory CD8+ T cell subtypes as controls from colorectal cancer patients, to compare their transcriptome and methylome characteristics. Transcriptome profiling confirms previous conclusions that tumor-reactive TILs have an exhausted tissue-resident memory signature. Whole-genome methylation profiling identifies a distinct methylome pattern of tumor-reactive CD8+ T cells, with tumor-reactive markers CD39 and CD103 being specifically demethylated. In addition, dynamic changes are observed during the transition of naïve T cells into tumor-reactive CD8+ T cells. Transcription factor binding motif enrichment analysis identifies several immune-related transcription factors, including three exhaustion-related genes (NR4A1, BATF, and EGR2) and VDR, which potentially play an important regulatory role in tumor-reactive CD8+ T cells. CONCLUSION: Our study supports the involvement of DNA methylation in shaping tumor-reactive and bystander CD8+ TILs, and provides a valuable resource for the development of novel DNA methylation markers and future therapeutics.

Understanding tumor ecosystems by single-cell sequencing: promises and limitations.

Cellular heterogeneity within and across tumors has been a major obstacle in understanding and treating cancer, and the complex heterogeneity is masked if bulk tumor tissues are used for analysis. The advent of rapidly developing single-cell sequencing technologies, which include methods related to single-cell genome, epigenome, transcriptome, and multi-omics sequencing, have been applied to cancer research and led to exciting new findings in the fields of cancer evolution, metastasis, resistance to therapy, and tumor microenvironment. In this review, we discuss recent advances and limitations of these new technologies and their potential applications in cancer studies.

Lineage tracking reveals dynamic relationships of T cells in colorectal cancer.

T cells are key elements of cancer immunotherapy1 but certain fundamental properties, such as the development and migration of T cells within tumours, remain unknown. The enormous T cell receptor (TCR) repertoire, which is required for the recognition of foreign and self-antigens2, could serve as lineage tags to track these T cells in tumours3. Here we obtained transcriptomes of 11,138 single T cells from 12 patients with colorectal cancer, and developed single T cell analysis by RNA sequencing and TCR tracking (STARTRAC) indices to quantitatively analyse the dynamic relationships among 20 identified T cell subsets with distinct functions and clonalities. Although both CD8+ effector and 'exhausted' T cells exhibited high clonal expansion, they were independently connected with tumour-resident CD8+ effector memory cells, implicating a TCR-based fate decision. Of the CD4+ T cells, most tumour-infiltrating T regulatory (Treg) cells showed clonal exclusivity, whereas certain Treg cell clones were developmentally linked to several T helper (TH) cell clones. Notably, we identified two IFNG+ TH1-like cell clusters in tumours that were associated with distinct IFNγ-regulating transcription factors -the GZMK+ effector memory T cells, which were associated with EOMES and RUNX3, and CXCL13+BHLHE40+ TH1-like cell clusters, which were associated with BHLHE40. Only CXCL13+BHLHE40+ TH1-like cells were preferentially enriched in patients with microsatellite-instable tumours, and this might explain their favourable responses to immune-checkpoint blockade. Furthermore, IGFLR1 was highly expressed in both CXCL13+BHLHE40+ TH1-like cells and CD8+ exhausted T cells and possessed co-stimulatory functions. Our integrated STARTRAC analyses provide a powerful approach to dissect the T cell properties in colorectal cancer comprehensively, and could provide insights into the dynamic relationships of T cells in other cancers.

Publisher Correction: Global characterization of T cells in non-small-cell lung cancer by single-cell sequencing.

In the version of this article originally published, the P statistic described in Fig. 3d was incorrect. It was described as "P < 22 × 10-16". It should have been "P < 2.2 × 10-16". Also, the "CD8+ Treg" label in Fig. 4f was incorrect. It should have been "CD4+ Treg". The errors have been corrected in the HTML and PDF versions of this article.

Global characterization of T cells in non-small-cell lung cancer by single-cell sequencing.

Cancer immunotherapies have shown sustained clinical responses in treating non-small-cell lung cancer1-3, but efficacy varies and depends in part on the amount and properties of tumor infiltrating lymphocytes4-6. To depict the baseline landscape of the composition, lineage and functional states of tumor infiltrating lymphocytes, here we performed deep single-cell RNA sequencing for 12,346 T cells from 14 treatment-naïve non-small-cell lung cancer patients. Combined expression and T cell antigen receptor based lineage tracking revealed a significant proportion of inter-tissue effector T cells with a highly migratory nature. As well as tumor-infiltrating CD8+ T cells undergoing exhaustion, we observed two clusters of cells exhibiting states preceding exhaustion, and a high ratio of "pre-exhausted" to exhausted T cells was associated with better prognosis of lung adenocarcinoma. Additionally, we observed further heterogeneity within the tumor regulatory T cells (Tregs), characterized by the bimodal distribution of TNFRSF9, an activation marker for antigen-specific Tregs. The gene signature of those activated tumor Tregs, which included IL1R2, correlated with poor prognosis in lung adenocarcinoma. Our study provides a new approach for patient stratification and will help further understand the functional states and dynamics of T cells in lung cancer.

Landscape of Infiltrating T Cells in Liver Cancer Revealed by Single-Cell Sequencing.

Systematic interrogation of tumor-infiltrating lymphocytes is key to the development of immunotherapies and the prediction of their clinical responses in cancers. Here, we perform deep single-cell RNA sequencing on 5,063 single T cells isolated from peripheral blood, tumor, and adjacent normal tissues from six hepatocellular carcinoma patients. The transcriptional profiles of these individual cells, coupled with assembled T cell receptor (TCR) sequences, enable us to identify 11 T cell subsets based on their molecular and functional properties and delineate their developmental trajectory. Specific subsets such as exhausted CD8+ T cells and Tregs are preferentially enriched and potentially clonally expanded in hepatocellular carcinoma (HCC), and we identified signature genes for each subset. One of the genes, layilin, is upregulated on activated CD8+ T cells and Tregs and represses the CD8+ T cell functions in vitro. This compendium of transcriptome data provides valuable insights and a rich resource for understanding the immune landscape in cancers.

GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses.

Tremendous amount of RNA sequencing data have been produced by large consortium projects such as TCGA and GTEx, creating new opportunities for data mining and deeper understanding of gene functions. While certain existing web servers are valuable and widely used, many expression analysis functions needed by experimental biologists are still not adequately addressed by these tools. We introduce GEPIA (Gene Expression Profiling Interactive Analysis), a web-based tool to deliver fast and customizable functionalities based on TCGA and GTEx data. GEPIA provides key interactive and customizable functions including differential expression analysis, profiling plotting, correlation analysis, patient survival analysis, similar gene detection and dimensionality reduction analysis. The comprehensive expression analyses with simple clicking through GEPIA greatly facilitate data mining in wide research areas, scientific discussion and the therapeutic discovery process. GEPIA fills in the gap between cancer genomics big data and the delivery of integrated information to end users, thus helping unleash the value of the current data resources. GEPIA is available at http://gepia.cancer-pku.cn/.