A pan-cancer single-cell transcriptional atlas of tumor infiltrating myeloid cells.

Tumor-infiltrating myeloid cells (TIMs) are key regulators in tumor progression, but the similarity and distinction of their fundamental properties across different tumors remain elusive. Here, by performing a pan-cancer analysis of single myeloid cells from 210 patients across 15 human cancer types, we identified distinct features of TIMs across cancer types. Mast cells in nasopharyngeal cancer were found to be associated with better prognosis and exhibited an anti-tumor phenotype with a high ratio of TNF+/VEGFA+ cells. Systematic comparison between cDC1- and cDC2-derived LAMP3+ cDCs revealed their differences in transcription factors and external stimulus. Additionally, pro-angiogenic tumor-associated macrophages (TAMs) were characterized with diverse markers across different cancer types, and the composition of TIMs appeared to be associated with certain features of somatic mutations and gene expressions. Our results provide a systematic view of the highly heterogeneous TIMs and suggest future avenues for rational, targeted immunotherapies.

Single-Cell Analyses Inform Mechanisms of Myeloid-Targeted Therapies in Colon Cancer.

Single-cell RNA sequencing (scRNA-seq) is a powerful tool for defining cellular diversity in tumors, but its application toward dissecting mechanisms underlying immune-modulating therapies is scarce. We performed scRNA-seq analyses on immune and stromal populations from colorectal cancer patients, identifying specific macrophage and conventional dendritic cell (cDC) subsets as key mediators of cellular cross-talk in the tumor microenvironment. Defining comparable myeloid populations in mouse tumors enabled characterization of their response to myeloid-targeted immunotherapy. Treatment with anti-CSF1R preferentially depleted macrophages with an inflammatory signature but spared macrophage populations that in mouse and human expresses pro-angiogenic/tumorigenic genes. Treatment with a CD40 agonist antibody preferentially activated a cDC population and increased Bhlhe40+ Th1-like cells and CD8+ memory T cells. Our comprehensive analysis of key myeloid subsets in human and mouse identifies critical cellular interactions regulating tumor immunity and defines mechanisms underlying myeloid-targeted immunotherapies currently undergoing clinical testing.

Landscape and Dynamics of Single Immune Cells in Hepatocellular Carcinoma.

The immune microenvironment of hepatocellular carcinoma (HCC) is poorly characterized. Combining two single-cell RNA sequencing technologies, we produced transcriptomes of CD45+ immune cells for HCC patients from five immune-relevant sites: tumor, adjacent liver, hepatic lymph node (LN), blood, and ascites. A cluster of LAMP3+ dendritic cells (DCs) appeared to be the mature form of conventional DCs and possessed the potential to migrate from tumors to LNs. LAMP3+ DCs also expressed diverse immune-relevant ligands and exhibited potential to regulate multiple subtypes of lymphocytes. Of the macrophages in tumors that exhibited distinct transcriptional states, tumor-associated macrophages (TAMs) were associated with poor prognosis, and we established the inflammatory role of SLC40A1 and GPNMB in these cells. Further, myeloid and lymphoid cells in ascites were predominantly linked to tumor and blood origins, respectively. The dynamic properties of diverse CD45+ cell types revealed by this study add new dimensions to the immune landscape of HCC.

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.

Single-cell exome sequencing reveals single-nucleotide mutation characteristics of a kidney tumor.

Clear cell renal cell carcinoma (ccRCC) is the most common kidney cancer and has very few mutations that are shared between different patients. To better understand the intratumoral genetics underlying mutations of ccRCC, we carried out single-cell exome sequencing on a ccRCC tumor and its adjacent kidney tissue. Our data indicate that this tumor was unlikely to have resulted from mutations in VHL and PBRM1. Quantitative population genetic analysis indicates that the tumor did not contain any significant clonal subpopulations and also showed that mutations that had different allele frequencies within the population also had different mutation spectrums. Analyses of these data allowed us to delineate a detailed intratumoral genetic landscape at a single-cell level. Our pilot study demonstrates that ccRCC may be more genetically complex than previously thought and provides information that can lead to new ways to investigate individual tumors, with the aim of developing more effective cellular targeted therapies.

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.

Single-cell RNA sequencing reveals intrahepatic and peripheral immune characteristics related to disease phases in HBV-infected patients.

OBJECTIVE: A comprehensive immune landscape for HBV infection is pivotal to achieve HBV cure. DESIGN: We performed single-cell RNA sequencing of 2 43 000 cells from 46 paired liver and blood samples of 23 individuals, including six immune tolerant, 5 immune active (IA), 3 acute recovery (AR), 3 chronic resolved and 6 HBV-free healthy controls (HCs). Flow cytometry and histological assays were applied in a second HBV cohort for validation. RESULTS: Both IA and AR were characterised by high levels of intrahepatic exhausted CD8+ T (Tex) cells. In IA, Tex cells were mainly derived from liver-resident GZMK+ effector memory T cells and self-expansion. By contrast, peripheral CX3CR1+ effector T cells and GZMK+ effector memory T cells were the main source of Tex cells in AR. In IA but not AR, significant cell-cell interactions were observed between Tex cells and regulatory CD4+ T cells, as well as between Tex and FCGR3A+ macrophages. Such interactions were potentially mediated through human leukocyte antigen class I molecules together with their receptors CANX and LILRBs, respectively, contributing to the dysfunction of antiviral immune responses. By contrast, CX3CR1+GNLY+ central memory CD8+ T cells were concurrently expanded in both liver and blood of AR, providing a potential surrogate marker for viral resolution. In clinic, intrahepatic Tex cells were positively correlated with serum alanine aminotransferase levels and histological grading scores. CONCLUSION: Our study dissects the coordinated immune responses for different HBV infection phases and provides a rich resource for fully understanding immunopathogenesis and developing effective therapeutic strategies.

Resolving the intertwining of inflammation and fibrosis in human heart failure at single-cell level.

Inflammation and fibrosis are intertwined mechanisms fundamentally involved in heart failure. Detailed deciphering gene expression perturbations and cell-cell interactions of leukocytes and non-myocytes is required to understand cell-type-specific pathology in the failing human myocardium. To this end, we performed single-cell RNA sequencing and single T cell receptor sequencing of 200,615 cells in both human dilated cardiomyopathy (DCM) and ischemic cardiomyopathy (ICM) hearts. We sampled both lesion and mild-lesion tissues from each heart to sequentially capture cellular and molecular alterations to different extents of cardiac fibrosis. By which, left (lesion) and right ventricle (mild-lesion) for DCM hearts were harvest while infarcted (lesion) and non-infarcted area (mild-lesion) were dissected from ICM hearts. A novel transcription factor AEBP1 was identified as a crucial cardiac fibrosis regulator in ACTA2+ myofibroblasts. Within fibrotic myocardium, an infiltration of a considerable number of leukocytes was witnessed, especially cytotoxic and exhausted CD8+ T cells and pro-inflammatory CD4+ T cells. Furthermore, a subset of tissue-resident macrophage, CXCL8hiCCR2+HLA-DRhi macrophage was particularly identified in severely fibrotic area, which interacted with activated endothelial cell via DARC, that potentially facilitate leukocyte recruitment and infiltration in human heart failure.

Transcriptional mutagenesis mediated by 8-oxoG induces translational errors in mammalian cells.

Reactive oxygen species formed within the mammalian cell can produce 8-oxo-7,8-dihydroguanine (8-oxoG) in mRNA, which can cause base mispairing during gene expression. Here we found that administration of 8-oxoGTP in MTH1-knockdown cells results in increased 8-oxoG content in mRNA. Under this condition, an amber mutation of the reporter luciferase is suppressed. Using second-generation sequencing techniques, we found that U-to-G changes at preassigned sites of the luciferase transcript increased when 8-oxoGTP was supplied. In addition, an increased level of 8-oxoG content in RNA induced the accumulation of aggregable amyloid ß peptides in cells expressing amyloid precursor protein. Our findings indicate that 8-oxoG accumulation in mRNA can alter protein synthesis in mammalian cells. Further work is required to assess the significance of these findings under normal physiological conditions.

Dynamics of peripheral T cell clones during PD-1 blockade in non-small cell lung cancer.

Understanding of the functional states and clonal dynamics of T cells after immune checkpoint blockade (ICB) is valuable for improving these therapeutic strategies. Here we performed Smart-seq2 single-cell RNA sequencing (scRNA-seq) analysis on 3,110 peripheral T cells of non-small cell lung cancer (NSCLC) patients before and after the initiation of programmed cell death protein 1 (PD-1) blockade. We identified individual peripheral T cell clones based on the full-length T cell receptor (TCR) sequences and monitored their dynamics during immunotherapy. We found a higher cytotoxic activity in the tumor-related CD4+ T cell clones than in the CD8+ T cell clones. Based on a large tumor-related CD4+ T cell clone, we observed a dramatically decreased abundance after progression, as well as a reduction in the percentage of PD-1+ T cells. We also detected 25 genes, such as CXCR4, DUSP2 and ZFP36, that were noticeably upregulated or downregulated following progression. In addition, the pseudotime trajectory of CD8+ T cell clones corresponded to the treatment time points, showing a decreased activity in the "cytokine and cytokine receptor interaction" pathway. These analyses provided an insight into the dynamics of peripheral T cell clones during PD-1 blockade in NSCLC.

Understanding the Genetic Mechanisms of Cancer Drug Resistance Using Genomic Approaches.

A major obstacle in precision cancer medicine is the inevitable resistance to targeted therapies. Tremendous effort and progress has been made over the past few years to understand the biochemical and genetic mechanisms underlying drug resistance, with the goal to eventually overcome such daunting challenges. Diverse mechanisms, such as secondary mutations, oncogene bypass, and epigenetic alterations, can all lead to drug resistance, and the number of known involved genes is growing rapidly, thus providing many possibilities to overcome resistance. The finding of these mechanisms and genes invariably requires the application of genomic and functional genomic approaches to tumors or cancer models. In this review, we briefly highlight the major drug-resistance mechanisms known today, and then focus primarily on the technological approaches leading to the advancement of this field.

GE-mini: a mobile APP for large-scale gene expression visualization.

Summary: The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) projects produced large-scale RNA sequencing data, which provides an opportunity for performing integrated expression analysis for all genes across tens of thousands of tumor and normal tissue specimens. Rapid access to and easy visualization of such valuable data could facilitate research in a wide biological area. Here, we present the GE-mini APP for smart phones, a mobile visualization tool for integrated gene expression data based on both TCGA and GTEx. This gene-centric expression viewer provides a convenient method for displaying expression profiles of all available tumor and tissue types, while allowing drilling down to detailed views for specific tissue types. Availability and Implementation: Both the iOS and Android APPs are freely available to all non-commercial users in App Store and Google Play. The QR codes of App store and Google play are also provided for scanning and download. The GE-mini web server is also available at http://gemini.cancer-pku.cn/ . Contacts: tangzefang@pku.edu.cn or huxueda@pku.edu.cn. Supplementary information: Supplementary data are available at Bioinformatics online.

RNA editing is induced by type I interferon in esophageal squamous cell carcinoma.

In recent years, abnormal RNA editing has been shown to play an important role in the development of esophageal squamous cell carcinoma, as such abnormal editing is catalyzed by ADAR (adenosine deaminases acting on RNA). However, the regulatory mechanism of ADAR1 in esophageal squamous cell carcinomas remains largely unknown. In this study, we investigated ADAR1 expression and its association with RNA editing in esophageal squamous cell carcinomas. RNA sequencing applied to esophageal squamous cell carcinoma clinical samples showed that ADAR1 expression was correlated with the expression of STAT1, STAT2, and IRF9. In vitro experiments showed that the abundance of ADAR1 protein was associated with the induced activation of the JAK/STAT pathway by type I interferon. RNA sequencing results showed that treatment with type I interferon caused an increase in the number and degree of RNA editing in esophageal squamous cell carcinoma cell lines. In conclusion, the activation of the JAK/STAT pathway is a regulatory mechanism of ADAR1 expression and causes abnormal RNA editing profile in esophageal squamous cell carcinoma. This mechanism may serve as a new target for esophageal squamous cell carcinoma therapy.

RNA editing of AZIN1 induces the malignant progression of non-small-cell lung cancers.

RNA editing is a widespread post-transcriptional mechanism that confers specific and reproducible nucleotide changes in selected RNA transcripts and plays a critical role in many human cancers. However, little is known about how RNA editing operates in non-small-cell lung cancers. Here, we measured the sequence and expression level of genes of antizyme inhibitor 1 and adenosine deaminase acting on RNA family in 30 non-small-cell lung cancer patient samples and 13 cell lines and revealed RNA editing S367G in antizyme inhibitor 1 is a high-frequent molecular events. We determined overexpression of antizyme inhibitor 1 with RNA editing, implying the oncogenic function of this alteration. We also detected the association of adenosine deaminase acting on RNA overexpression with RNA editing occurred in antizyme inhibitor 1. Furthermore, the RNA editing could cause a cytoplasmic-to-nuclear translocation of antizyme inhibitor 1 protein and conferred the malignant phenotype of non-small-cell lung cancer cells. The in vivo experiment confirmed that this RNA editing confers higher capacity of tumor migration as well. In conclusion, antizyme inhibitor 1 RNA editing and its involvement in tumorigenesis of non-small-cell lung cancer pave a new way for potential clinical management of non-small-cell lung cancer.

Combination of platelet count and mean platelet volume (COP-MPV) predicts postoperative prognosis in both resectable early and advanced stage esophageal squamous cell cancer patients.

The aim of this study is to search the most powerful prognostic factor from routine blood test for esophageal squamous cell cancer (ESCC) patients. Multiple laboratory tests were evaluated including those reflecting red blood cell parameters (hemoglobin (Hb), mean corpuscular volume (MCV), mean corpuscular hemoglobin concentration (MCHC), and red blood cell distribution width (RDW)), platelet morphological parameters (mean platelet volume (MPV) and platelet count (PLT)), blood coagulation status (D-dimer), and tumor biomarker (CA19-9). Known inflammatory indices (NLR and PLR) were also calculated. A total of 468 patients who were diagnosed with ESCC between December 2005 and December 2008 were retrospectively analyzed in this study. By utilizing univariate and multivariate Cox proportional hazard analyses, we found that PLT and MPV were significantly associated with overall survival (OS) and disease-free survival (DFS) of ESCC patients, with optimal cutoff values of 212 and 10.6, respectively. Moreover, the combination of the preoperative PLT and MPV (COP-MPV) was calculated as follows: patients with both PLT (≥212 × 10(9) L(-1)) and MPV (≥10.6 fL) elevation were assigned a score of 2, and patients with one or neither were assigned a score of 1 and 0. The COP-MPV was an independent prognostic factor for OS (hazard ratio (HR) 0.378, 95 % confidence interval (CI) 0.241 to 0.593, P < 0.001, 0/2) and DFS (HR 0.341, 95 % CI 0.218 to 0.534, P < 0.001, 0/2) in multivariate analyses. In subgroup analyses for early (stages I and II) and locally (stage III) advanced stage patients, COP-MPV was found significantly associated with OS and DFS in each group (P = 0.025 and P = 0.018 for OS and P = 0.029 and P = 0.002 for DFS). In conclusion, we considered that COP-MPV is a promising predictor for postoperative survival in ESCC patients.

Multilayered molecular profiling supported the monoclonal origin of metastatic renal cell carcinoma.

Primary renal cell carcinomas (pRCCs) have a high degree of intratumoral heterogeneity and are composed of multiple distinct subclones. However, it remains largely unknown that whether metastatic renal cell carcinomas (mRCCs) also have startling intratumoral heterogeneity or whether development of mRCCs is due to early dissemination or late diagnosis. To decipher the evolution of mRCC, we analyzed the multilayered molecular profiles of pRCC, local invasion of the vena cava (IVC), and distant metastasis to the brain (MB) from the same patient using whole-genome sequencing, whole-exome sequencing, DNA methylome profiling, and transcriptome sequencing. We found that mRCC had a lower degree of heterogeneity than pRCC and was likely to result from recent clonal expansion of a rare, advantageous subclone. Consequently, some key pathways that are targeted by clinically available drugs showed distinct expression patterns between pRCC and mRCC. From the genetic distances between different tumor subclones, we estimated that the progeny subclone giving rise to distant metastasis took over half a decade to acquire the full potential of metastasis since the birth of the subclone that evolved into IVC. Our evidence supported that mRCC was monoclonal and distant metastasis occurred late during renal cancer progression. Thus, there was a broad window for early detection of circulating tumor cells and future targeted treatments for patients with mRCCs should rely on the molecular profiles of metastases.

Deciphering Immune Landscape Remodeling Unravels the Underlying Mechanism for Synchronized Muscle and Bone Aging.

Evidence from numerous studies has revealed the synchronous progression of aging in bone and muscle; however, little is known about the underlying mechanisms. To this end, human muscles and bones are harvested and the aging-associated transcriptional dynamics of two tissues in parallel using single-cell RNA sequencing are surveyed. A subset of lipid-associated macrophages (triggering receptor expressed on myeloid cells 2, TREM2+ Macs) is identified in both aged muscle and bone. Genes responsible for muscle dystrophy and bone loss, such as secreted phosphoprotein 1 (SPP1), are also highly expressed in TREM2+ Macs, suggesting its conserved role in aging-related features. A common transition toward pro-inflammatory phenotypes in aged CD4+ T cells across tissues is also observed, activated by the nuclear factor kappa B subunit 1 (NFKB1). CD4+ T cells in aged muscle experience Th1-like differentiation, whereas, in bone, a skewing toward Th17 cells is observed. Furthermore, these results highlight that degenerated myocytes produce BAG6-containing exosomes that can communicate with Th17 cells in the bone through its receptor natural cytotoxicity triggering receptor 3 (NCR3). This communication upregulates CD6 expression in Th17 cells, which then interact with TREM2+ Macs through CD6-ALCAM signaling, ultimately stimulating the transcription of SPP1 in TREM2+ Macs. The negative correlation between serum exosomal BCL2-associated athanogene 6 (BAG6) levels and bone mineral density further supports its role in mediating muscle and bone synchronization with aging.

Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome.

Understanding the dynamics of eukaryotic transcriptome is essential for studying the complexity of transcriptional regulation and its impact on phenotype. However, comprehensive studies of transcriptomes at single base resolution are rare, even for modern organisms, and lacking for rice. Here, we present the first transcriptome atlas for eight organs of cultivated rice. Using high-throughput paired-end RNA-seq, we unambiguously detected transcripts expressing at an extremely low level, as well as a substantial number of novel transcripts, exons, and untranslated regions. An analysis of alternative splicing in the rice transcriptome revealed that alternative cis-splicing occurred in approximately 33% of all rice genes. This is far more than previously reported. In addition, we also identified 234 putative chimeric transcripts that seem to be produced by trans-splicing, indicating that transcript fusion events are more common than expected. In-depth analysis revealed a multitude of fusion transcripts that might be by-products of alternative splicing. Validation and chimeric transcript structural analysis provided evidence that some of these transcripts are likely to be functional in the cell. Taken together, our data provide extensive evidence that transcriptional regulation in rice is vastly more complex than previously believed.

The DNA methylome of human peripheral blood mononuclear cells.

DNA methylation plays an important role in biological processes in human health and disease. Recent technological advances allow unbiased whole-genome DNA methylation (methylome) analysis to be carried out on human cells. Using whole-genome bisulfite sequencing at 24.7-fold coverage (12.3-fold per strand), we report a comprehensive (92.62%) methylome and analysis of the unique sequences in human peripheral blood mononuclear cells (PBMC) from the same Asian individual whose genome was deciphered in the YH project. PBMC constitute an important source for clinical blood tests world-wide. We found that 68.4% of CpG sites and <0.2% of non-CpG sites were methylated, demonstrating that non-CpG cytosine methylation is minor in human PBMC. Analysis of the PBMC methylome revealed a rich epigenomic landscape for 20 distinct genomic features, including regulatory, protein-coding, non-coding, RNA-coding, and repeat sequences. Integration of our methylome data with the YH genome sequence enabled a first comprehensive assessment of allele-specific methylation (ASM) between the two haploid methylomes of any individual and allowed the identification of 599 haploid differentially methylated regions (hDMRs) covering 287 genes. Of these, 76 genes had hDMRs within 2 kb of their transcriptional start sites of which >80% displayed allele-specific expression (ASE). These data demonstrate that ASM is a recurrent phenomenon and is highly correlated with ASE in human PBMCs. Together with recently reported similar studies, our study provides a comprehensive resource for future epigenomic research and confirms new sequencing technology as a paradigm for large-scale epigenomics studies.

Single-cell meta-analyses reveal responses of tumor-reactive CXCL13+ T cells to immune-checkpoint blockade.

Immune-checkpoint blockade (ICB) therapies represent a paradigm shift in the treatment of human cancers; however, it remains incompletely understood how tumor-reactive T cells respond to ICB across tumor types. Here, we demonstrate that measuring CXCL13 expression could effectively identify both precursor and terminally differentiated tumor-reactive CD8+ T cells within tumors. Applying this approach, we performed meta-analyses of published single-cell data for CXCL13+CD8+ T cells in 225 samples from 102 patients treated with ICB across five cancer types. We found that CXCL13+CD8+ T cells were correlated with favorable responses to ICB, and the treatment further increased such cells in responsive tumors. In addition, CXCL13+ tumor-reactive subsets exhibited variable responses to ICB in distinct contexts, likely due to different degrees of exhaustion-related immunosuppression. Our integrated analyses provide insights into mechanisms underlying ICB and suggest that bolstering precursor tumor-reactive CD8+ T cells might provide an effective therapeutic approach to improve cancer treatment.