Insights Gained from Single-Cell Analysis of Immune Cells in the Tumor Microenvironment.

Understanding tumor immune microenvironments is critical for identifying immune modifiers of cancer progression and developing cancer immunotherapies. Recent applications of single-cell RNA sequencing (scRNA-seq) in dissecting tumor microenvironments have brought important insights into the biology of tumor-infiltrating immune cells, including their heterogeneity, dynamics, and potential roles in both disease progression and response to immune checkpoint inhibitors and other immunotherapies. This review focuses on the advances in knowledge of tumor immune microenvironments acquired from scRNA-seq studies across multiple types of human tumors, with a particular emphasis on the study of phenotypic plasticity and lineage dynamics of immune cells in the tumor environment. We also discuss several imminent questions emerging from scRNA-seq observations and their potential solutions on the horizon.

Pan-cancer single-cell dissection reveals phenotypically distinct B cell subtypes.

Characterizing the compositional and phenotypic characteristics of tumor-infiltrating B cells (TIBs) is important for advancing our understanding of their role in cancer development. Here, we establish a comprehensive resource of human B cells by integrating single-cell RNA sequencing data of B cells from 649 patients across 19 major cancer types. We demonstrate substantial heterogeneity in their total abundance and subtype composition and observe immunoglobulin G (IgG)-skewness of antibody-secreting cell isotypes. Moreover, we identify stress-response memory B cells and tumor-associated atypical B cells (TAABs), two tumor-enriched subpopulations with prognostic potential, shared in a pan-cancer manner. In particular, TAABs, characterized by a high clonal expansion level and proliferative capacity as well as by close interactions with activated CD4 T cells in tumors, are predictive of immunotherapy response. Our integrative resource depicts distinct clinically relevant TIB subsets, laying a foundation for further exploration of functional commonality and diversity of B cells in cancer.

Accelerating the understanding of cancer biology through the lens of genomics.

A core mission of cancer genomics is to comprehensively chart molecular underpinnings of cancer-driving events and to provide personalized therapeutic strategies. Primarily focused on cancer cells, cancer genomics studies have successfully uncovered many drivers for major cancer types. Since the emergence of cancer immune evasion as a critical cancer hallmark, the paradigm has been elevated to the holistic tumor ecosystem, with distinct cellular components and their functional states elucidated. We highlight the milestones of cancer genomics, depict the evolving path of the field, and discuss future directions in completing the understanding of the tumor ecosystem and in advancing therapeutic strategies.

A pan-cancer single-cell panorama of human natural killer cells.

Natural killer (NK) cells play indispensable roles in innate immune responses against tumor progression. To depict their phenotypic and functional diversities in the tumor microenvironment, we perform integrative single-cell RNA sequencing analyses on NK cells from 716 patients with cancer, covering 24 cancer types. We observed heterogeneity in NK cell composition in a tumor-type-specific manner. Notably, we have identified a group of tumor-associated NK cells that are enriched in tumors, show impaired anti-tumor functions, and are associated with unfavorable prognosis and resistance to immunotherapy. Specific myeloid cell subpopulations, in particular LAMP3+ dendritic cells, appear to mediate the regulation of NK cell anti-tumor immunity. Our study provides insights into NK-cell-based cancer immunity and highlights potential clinical utilities of NK cell subsets as therapeutic targets.

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.

Decoding the genetic basis of anti-tumor immunity.

Substantial variations of tumor immune properties exist among cancer patients, but the contributing factors underlying those variations are poorly understood. In this issue of Immunity, Sayaman et al. uncover associations between germline genetic variants and tumor immune properties, revealing candidate causal genes.

SARS-CoV-2 exacerbates proinflammatory responses in myeloid cells through C-type lectin receptors and Tweety family member 2.

Despite mounting evidence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) engagement with immune cells, most express little, if any, of the canonical receptor of SARS-CoV-2, angiotensin-converting enzyme 2 (ACE2). Here, using a myeloid cell receptor-focused ectopic expression screen, we identified several C-type lectins (DC-SIGN, L-SIGN, LSECtin, ASGR1, and CLEC10A) and Tweety family member 2 (TTYH2) as glycan-dependent binding partners of the SARS-CoV-2 spike. Except for TTYH2, these molecules primarily interacted with spike via regions outside of the receptor-binding domain. Single-cell RNA sequencing analysis of pulmonary cells from individuals with coronavirus disease 2019 (COVID-19) indicated predominant expression of these molecules on myeloid cells. Although these receptors do not support active replication of SARS-CoV-2, their engagement with the virus induced robust proinflammatory responses in myeloid cells that correlated with COVID-19 severity. We also generated a bispecific anti-spike nanobody that not only blocked ACE2-mediated infection but also the myeloid receptor-mediated proinflammatory responses. Our findings suggest that SARS-CoV-2-myeloid receptor interactions promote immune hyperactivation, which represents potential targets for COVID-19 therapy.

Liver tumour immune microenvironment subtypes and neutrophil heterogeneity.

The heterogeneity of the tumour immune microenvironment (TIME), organized by various immune and stromal cells, is a major contributing factor of tumour metastasis, relapse and drug resistance1-3, but how different TIME subtypes are connected to the clinical relevance in liver cancer remains unclear. Here we performed single-cell RNA-sequencing (scRNA-seq) analysis of 189 samples collected from 124 patients and 8 mice with liver cancer. With more than 1 million cells analysed, we stratified patients into five TIME subtypes, including immune activation, immune suppression mediated by myeloid or stromal cells, immune exclusion and immune residence phenotypes. Different TIME subtypes were spatially organized and associated with chemokine networks and genomic features. Notably, tumour-associated neutrophil (TAN) populations enriched in the myeloid-cell-enriched subtype were associated with an unfavourable prognosis. Through in vitro induction of TANs and ex vivo analyses of patient TANs, we showed that CCL4+ TANs can recruit macrophages and that PD-L1+ TANs can suppress T cell cytotoxicity. Furthermore, scRNA-seq analysis of mouse neutrophil subsets revealed that they are largely conserved with those of humans. In vivo neutrophil depletion in mouse models attenuated tumour progression, confirming the pro-tumour phenotypes of TANs. With this detailed cellular heterogeneity landscape of liver cancer, our study illustrates diverse TIME subtypes, highlights immunosuppressive functions of TANs and sheds light on potential immunotherapies targeting TANs.

Electrifying HCOOH synthesis from CO2 building blocks over Cu-Bi nanorod arrays.

Precise electrochemical synthesis of commodity chemicals and fuels from CO2 building blocks provides a promising route to close the anthropogenic carbon cycle, in which renewable but intermittent electricity could be stored within the greenhouse gas molecules. Here, we report state-of-the-art CO2-to-HCOOH valorization performance over a multiscale optimized Cu-Bi cathodic architecture, delivering a formate Faradaic efficiency exceeding 95% within an aqueous electrolyzer, a C-basis HCOOH purity above 99.8% within a solid-state electrolyzer operated at 100 mA cm-2 for 200 h and an energy efficiency of 39.2%, as well as a tunable aqueous HCOOH concentration ranging from 2.7 to 92.1 wt%. Via a combined two-dimensional reaction phase diagram and finite element analysis, we highlight the role of local geometries of Cu and Bi in branching the adsorption strength for key intermediates like *COOH and *OCHO for CO2 reduction, while the crystal orbital Hamiltonian population analysis rationalizes the vital contribution from moderate binding strength of η2(O,O)-OCHO on Cu-doped Bi surface in promoting HCOOH electrosynthesis. The findings of this study not only shed light on the tuning knobs for precise CO2 valorization, but also provide a different research paradigm for advancing the activity and selectivity optimization in a broad range of electrosynthetic systems.

Genomic basis for RNA alterations in cancer.

Transcript alterations often result from somatic changes in cancer genomes1. Various forms of RNA alterations have been described in cancer, including overexpression2, altered splicing3 and gene fusions4; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)5. Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis, of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed 'bridged' fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

The OsMYC2-JA feedback loop regulates diurnal flower-opening time via cell wall loosening in rice.

Diurnal flower-opening time (DFOT), the time of spikelet opening during the day, is an important trait for hybrid rice (Oryza sativa L.) seed production. Hybrids between indica and japonica rice varieties have strong heterosis, but the parental lines usually have different, nonoverlapping DFOTs. This reduces the success of hybrid seed production in crosses between indica and japonica subspecies, thus hindering the utilization of indica and japonica inter-subspecies heterosis. However, little is known about the molecular mechanisms regulating DFOT in rice. Here, we obtained japonica rice lines with a DFOT 1.5 h earlier than the wild type by overexpressing OsMYC2, a gene encoding a key transcription factor in the jasmonate (JA) signaling pathway. OsMYC2 is activated by JA signaling and directly regulates the transcription of genes related to JA biosynthesis and cell wall metabolism. Overexpressing OsMYC2 led to significantly increased JA contents and decreased cellulose and hemicellulose contents in lodicule cells, as well as the softening of lodicule cell walls. This may facilitate the swelling of lodicules, resulting in early diurnal flower-opening. These results suggest that the OsMYC2-JA feedback loop regulates DFOT in rice via cell wall remodeling. These findings shed light on the understanding of regulatory mechanism of the DFOT of plants, which should promote the development of indica and japonica varieties suitable for hybrid rice breeding.

High-quality lateral monolayer-multilayer graphene junction formed by selective laser thinning for self-powered photodetection.

The formation of an asymmetric junction is key to graphene-based photodetectors of high-sensitive photodetectability, because such a junction can not only facilitate the diffusion or drift of photogenerated carriers but also realize a self-powered operation. Here, a monolayer-multilayer graphene junction photodetector is accomplished by selectively thinning part of a multilayer graphene to a high-quality monolayer. Benefiting from the large photoabsorption cross section of multilayer graphene and strong asymmetry caused by the significant differences in optoelectronic properties between monolayer and multilayer graphene, the monolayer-multilayer graphene junction shows a 7-fold increase in short-circuit photocurrent as compared with that at the monolayer graphene-metal contact in scanning photocurrent images. The asymmetric configuration also enables the photodetector to work at zero bias with minimized dark current noise and stand-by power consumption. Under global illumination with visible light, a photoswitching ratio of 3.4 × 103, a responsivity of 8.8 mA W-1, a specific detectivity of 1.3 × 108Jones and a response time of 11 ns can be obtained, suggesting a promising photoresponse. Moreover, it is worth mentioning that such a performance enhancement is achieved without compromising the broadband spectral response of graphene photodetector and it is hence applicable for long wavelength spectral range including infrared and terahertz.

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.

Effects of carbonization temperature and time on the characteristics of carbonized sludge.

To investigate the influence of carbonization process parameters on the characteristics of municipal sludge carbonization products, this study selected carbonization temperatures of 300-700 °C and carbonization times of 0.5-1.5 h to carbonize municipal sludge. The results showed that with an increase in temperature and carbonization time, the sludge was carbonized more completely, and the structure and performance characteristics of the sludge changed significantly. Organic matter was continuously cracked, the amorphous nature of the material was reduced, its morphology was transformed into an increasing number of regular crystalline structures, and the content of carbon continued to decrease, from the initial 52.85 to 38.77%, while the content of inorganic species consisting continued to increase. The conductivity was reduced by 87.8%, and the degree of conversion of salt ions into their residual and insoluble states was significant. Natural water absorption in the sludge decreased from 8.13 to 1.29%, and hydrophobicity increased. The dry-basis higher calorific value decreased from 8,703 to 3,574 kJ/kg. Heavy metals were concentrated by a factor of 2-3, but the content of the available state was very low. The results of this study provide important technological support for the selection of suitable carbonization process conditions and for resource utilization.

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.

Single-Cell DNA Sequencing Reveals Punctuated and Gradual Clonal Evolution in Hepatocellular Carcinoma.

BACKGROUND & AIMS: Copy number alterations (CNAs), elicited by genome instability, are a major source of intratumor heterogeneity. How CNAs evolve in hepatocellular carcinoma (HCC) remains unknown. METHODS: We performed single-cell DNA sequencing (scDNA-seq) on 1275 cells isolated from 10 patients with HCC, ploidy-resolved scDNA-seq on 356 cells from 1 additional patient, and single-cell RNA sequencing on 27,344 cells from 3 additional patients. Three statistical fitting models were compared to investigate the CNA accumulation pattern. RESULTS: Cells in the tumor were categorized into the following 3 subpopulations: euploid, pseudoeuploid, and aneuploid. Our scDNA-seq analysis revealed that CNA accumulation followed a dual-phase copy number evolution model, that is, a punctuated phase followed by a gradual phase. Patients who exhibited prolonged gradual phase showed higher intratumor heterogeneity and worse disease-free survival. Integrating bulk RNA sequencing of 17 patients with HCC, published datasets of 1196 liver tumors, and immunohistochemical staining of 202 HCC tumors, we found that high expression of CAD, a gene involved in pyrimidine synthesis, was correlated with rapid tumorigenesis and reduced survival. The dual-phase copy number evolution model was validated by our single-cell RNA sequencing data and published scDNA-seq datasets of other cancer types. Furthermore, ploidy-resolved scDNA-seq revealed the common clonal origin of diploid- and polyploid-aneuploid cells, suggesting that polyploid tumor cells were generated by whole genome doubling of diploid tumor cells. CONCLUSIONS: Our work revealed a novel dual-phase copy number evolution model, showed HCC with longer gradual phase was more severe, identified CAD as a promising biomarker for early recurrence of HCC, and supported the diploid origin of polyploid HCC.