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.

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.

In vivo CRISPR screens identify the E3 ligase Cop1 as a modulator of macrophage infiltration and cancer immunotherapy target.

Despite remarkable clinical efficacy of immune checkpoint blockade (ICB) in cancer treatment, ICB benefits for triple-negative breast cancer (TNBC) remain limited. Through pooled in vivo CRISPR knockout (KO) screens in syngeneic TNBC mouse models, we found that deletion of the E3 ubiquitin ligase Cop1 in cancer cells decreases secretion of macrophage-associated chemokines, reduces tumor macrophage infiltration, enhances anti-tumor immunity, and strengthens ICB response. Transcriptomics, epigenomics, and proteomics analyses revealed that Cop1 functions through proteasomal degradation of the C/ebpδ protein. The Cop1 substrate Trib2 functions as a scaffold linking Cop1 and C/ebpδ, which leads to polyubiquitination of C/ebpδ. In addition, deletion of the E3 ubiquitin ligase Cop1 in cancer cells stabilizes C/ebpδ to suppress expression of macrophage chemoattractant genes. Our integrated approach implicates Cop1 as a target for improving cancer immunotherapy efficacy in TNBC by regulating chemokine secretion and macrophage infiltration in the tumor microenvironment.

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.

Dendritic cell type 3 arises from Ly6C+ monocyte-dendritic cell progenitors.

Conventional dendritic cells (cDCs) are professional antigen-presenting cells that control the adaptive immune response. Their subsets and developmental origins have been intensively investigated but are still not fully understood as their phenotypes, especially in the DC2 lineage and the recently described human DC3s, overlap with monocytes. Here, using LEGENDScreen to profile DC vs. monocyte lineages, we found sustained expression of FLT3 and CD45RB through the whole DC lineage, allowing DCs and their precursors to be distinguished from monocytes. Using fate mapping models, single-cell RNA sequencing and adoptive transfer, we identified a lineage of murine CD16/32+CD172a+ DC3, distinct from DC2, arising from Ly6C+ monocyte-DC progenitors (MDPs) through Lyz2+Ly6C+CD11c- pro-DC3s, whereas DC2s develop from common DC progenitors (CDPs) through CD7+Ly6C+CD11c+ pre-DC2s. Corresponding DC subsets, developmental stages, and lineages exist in humans. These findings reveal DC3 as a DC lineage phenotypically related to but developmentally different from monocytes and DC2s.

Single-Cell Sequencing of Peripheral Mononuclear Cells Reveals Distinct Immune Response Landscapes of COVID-19 and Influenza Patients.

The coronavirus disease 2019 (COVID-19) pandemic poses a current world-wide public health threat. However, little is known about its hallmarks compared to other infectious diseases. Here, we report the single-cell transcriptional landscape of longitudinally collected peripheral blood mononuclear cells (PBMCs) in both COVID-19- and influenza A virus (IAV)-infected patients. We observed increase of plasma cells in both COVID-19 and IAV patients and XIAP associated factor 1 (XAF1)-, tumor necrosis factor (TNF)-, and FAS-induced T cell apoptosis in COVID-19 patients. Further analyses revealed distinct signaling pathways activated in COVID-19 (STAT1 and IRF3) versus IAV (STAT3 and NFκB) patients and substantial differences in the expression of key factors. These factors include relatively increase of interleukin (IL)6R and IL6ST expression in COVID-19 patients but similarly increased IL-6 concentrations compared to IAV patients, supporting the clinical observations of increased proinflammatory cytokines in COVID-19 patients. Thus, we provide the landscape of PBMCs and unveil distinct immune response pathways in COVID-19 and IAV patients.

Peptide sequencing based on host-guest interaction-assisted nanopore sensing.

Direct protein sequencing technologies with improved sensitivity and throughput are still needed. Here, we propose an alternative method for peptide sequencing based on enzymatic cleavage and host-guest interaction-assisted nanopore sensing. We serendipitously discovered that the identity of any proteinogenic amino acid in a particular position of a phenylalanine-containing peptide could be determined via current blockage during translocation of the peptide through α-hemolysin nanopores in the presence of cucurbit[7]uril. Building upon this, we further present a proof-of-concept demonstration of peptide sequencing by sequentially cleaving off amino acids from C terminus of a peptide with carboxypeptidases, and then determining their identities and sequence with a peptide probe in nanopore. With future optimization, our results point to a different way of nanopore-based protein sequencing.

Zebrafish mylipb attenuates antiviral innate immunity through two synergistic mechanisms targeting transcription factor irf3.

IFN regulatory factor 3 (IRF3) is the transcription factor crucial for the production of type I IFN in viral defence and inflammatory responses. The activity of IRF3 is strictly modulated by post-translational modifications (PTMs) to effectively protect the host from infection while avoiding excessive immunopathology. Here, we report that zebrafish myosin-regulated light chain interacting protein b (mylipb) inhibits virus-induced type I IFN production via two synergistic mechanisms: induction of autophagic degradation of irf3 and reduction of irf3 phosphorylation. In vivo, mylipb-null zebrafish exhibit reduced lethality and viral mRNA levels compared to controls. At the cellular level, overexpression of mylipb significantly reduces cellular antiviral capacity, and promotes viral proliferation. Mechanistically, mylipb associates with irf3 and targets Lys 352 to increase K6-linked polyubiquitination, dependent on its E3 ubiquitin ligase activity, leading to autophagic degradation of irf3. Meanwhile, mylipb acts as a decoy substrate for the phosphokinase tbk1 to attenuate irf3 phosphorylation and cellular antiviral responses independent of its enzymatic activity. These findings support a critical role for zebrafish mylipb in the limitation of antiviral innate immunity through two synergistic mechanisms targeting irf3.

Correction: TDP-43 and other hnRNPs regulate cryptic exon inclusion of a key ALS/FTD risk gene, UNC13A.

[This corrects the article DOI: 10.1371/journal.pbio.3002028.].

TDP-43 and other hnRNPs regulate cryptic exon inclusion of a key ALS/FTD risk gene, UNC13A.

A major function of TAR DNA-binding protein-43 (TDP-43) is to repress the inclusion of cryptic exons during RNA splicing. One of these cryptic exons is in UNC13A, a genetic risk factor for amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). The accumulation of cryptic UNC13A in disease is heightened by the presence of a risk haplotype located within the cryptic exon itself. Here, we revealed that TDP-43 extreme N-terminus is important to repress UNC13A cryptic exon inclusion. Further, we found hnRNP L, hnRNP A1, and hnRNP A2B1 bind UNC13A RNA and repress cryptic exon inclusion, independently of TDP-43. Finally, higher levels of hnRNP L protein associate with lower burden of UNC13A cryptic RNA in ALS/FTD brains. Our findings suggest that while TDP-43 is the main repressor of UNC13A cryptic exon inclusion, other hnRNPs contribute to its regulation and may potentially function as disease modifiers.

A comprehensive assessment of cell type-specific differential expression methods in bulk data.

Accounting for cell type compositions has been very successful at analyzing high-throughput data from heterogeneous tissues. Differential gene expression analysis at cell type level is becoming increasingly popular, yielding biomarker discovery in a finer granularity within a particular cell type. Although several computational methods have been developed to identify cell type-specific differentially expressed genes (csDEG) from RNA-seq data, a systematic evaluation is yet to be performed. Here, we thoroughly benchmark six recently published methods: CellDMC, CARseq, TOAST, LRCDE, CeDAR and TCA, together with two classical methods, csSAM and DESeq2, for a comprehensive comparison. We aim to systematically evaluate the performance of popular csDEG detection methods and provide guidance to researchers. In simulation studies, we benchmark available methods under various scenarios of baseline expression levels, sample sizes, cell type compositions, expression level alterations, technical noises and biological dispersions. Real data analyses of three large datasets on inflammatory bowel disease, lung cancer and autism provide evaluation in both the gene level and the pathway level. We find that csDEG calling is strongly affected by effect size, baseline expression level and cell type compositions. Results imply that csDEG discovery is a challenging task itself, with room to improvements on handling low signal-to-noise ratio and low expression genes.

A novel statistical method for decontaminating T-cell receptor sequencing data.

The T-cell receptor (TCR) repertoire is highly diverse among the population and plays an essential role in initiating multiple immune processes. TCR sequencing (TCR-seq) has been developed to profile the T cell repertoire. Similar to other high-throughput experiments, contamination can happen during several steps of TCR-seq, including sample collection, preparation and sequencing. Such contamination creates artifacts in the data, leading to inaccurate or even biased results. Most existing methods assume 'clean' TCR-seq data as the starting point with no ability to handle data contamination. Here, we develop a novel statistical model to systematically detect and remove contamination in TCR-seq data. We summarize the observed contamination into two sources, pairwise and cross-cohort. For both sources, we provide visualizations and summary statistics to help users assess the severity of the contamination. Incorporating prior information from 14 existing TCR-seq datasets with minimum contamination, we develop a straightforward Bayesian model to statistically identify contaminated samples. We further provide strategies for removing the impacted sequences to allow for downstream analysis, thus avoiding any need to repeat experiments. Our proposed model shows robustness in contamination detection compared with a few off-the-shelf detection methods in simulation studies. We illustrate the use of our proposed method on two TCR-seq datasets generated locally.

Regional analysis to delineate intrasample heterogeneity with RegionalST.

MOTIVATION: Spatial transcriptomics has greatly contributed to our understanding of spatial and intra-sample heterogeneity, which could be crucial for deciphering the molecular basis of human diseases. Intra-tumor heterogeneity, e.g. may be associated with cancer treatment responses. However, the lack of computational tools for exploiting cross-regional information and the limited spatial resolution of current technologies present major obstacles to elucidating tissue heterogeneity. RESULTS: To address these challenges, we introduce RegionalST, an efficient computational method that enables users to quantify cell type mixture and interactions, identify sub-regions of interest, and perform cross-region cell type-specific differential analysis for the first time. Our simulations and real data applications demonstrate that RegionalST is an efficient tool for visualizing and analyzing diverse spatial transcriptomics data, thereby enabling accurate and flexible exploration of tissue heterogeneity. Overall, RegionalST provides a one-stop destination for researchers seeking to delve deeper into the intricacies of spatial transcriptomics data. AVAILABILITY AND IMPLEMENTATION: The implementation of our method is available as an open-source R/Bioconductor package with a user-friendly manual available at https://bioconductor.org/packages/release/bioc/html/RegionalST.html.

Extracellular RNA Sensing Mediates Inflammation and Organ Injury in a Murine Model of Polytrauma.

Severe traumatic injury leads to marked systemic inflammation and multiorgan injury. Endogenous drivers such as extracellular nucleic acid may play a role in mediating innate immune response and the downstream pathogenesis. Here, we explored the role of plasma extracellular RNA (exRNA) and its sensing mechanism in inflammation and organ injury in a murine model of polytrauma. We found that severe polytrauma-bone fracture, muscle crush injury, and bowel ischemia-induced a marked increase in plasma exRNA, systemic inflammation, and multiorgan injury in mice. Plasma RNA profiling with RNA sequencing in mice and humans revealed a dominant presence of miRNAs and marked differential expression of numerous miRNAs after severe trauma. Plasma exRNA isolated from trauma mice induced a dose-dependent cytokine production in macrophages, which was almost abolished in TLR7-deficient cells but unchanged in TLR3-deficient cells. Moreover, RNase or specific miRNA inhibitors against the selected proinflammatory miRNAs (i.e., miR-7a-5p, miR-142, let-7j, miR-802, and miR-146a-5p) abolished or attenuated trauma plasma exRNA-induced cytokine production, respectively. Bioinformatic analyses of a group of miRNAs based on cytokine readouts revealed that high uridine abundance (>40%) is a reliable predictor in miRNA mimic-induced cytokine and complement production. Finally, compared with the wild-type, TLR7-knockout mice had attenuated plasma cytokine storm and reduced lung and hepatic injury after polytrauma. These data suggest that endogenous plasma exRNA of severely injured mice and ex-miRNAs with high uridine abundance prove to be highly proinflammatory. TLR7 sensing of plasma exRNA and ex-miRNAs activates innate immune responses and plays a role in inflammation and organ injury after trauma.

Altered hydroxymethylome in the substantia nigra of Parkinson's disease.

Parkinson's disease (PD) is the second most common neurodegenerative disorder, and aging and genetic and environmental exposure can contribute to its pathogenesis. DNA methylation has been suggested to play a pivotal role in neurodevelopment and neurodegenerative diseases. 5-hydroxymethylcytosine (5hmC) is generated through 5-methylcytosine (5mC) oxidization by ten-eleven translocation proteins and is particularly enriched in the brain. Although 5hmC has been linked to multiple neurological disorders, little is known about 5hmC alterations in the substantia nigra of patients with PD. To determine the specific alterations in DNA methylation and hydroxymethylation in PD brain samples, we examined the genome-wide profiles of 5mC and 5hmC in the substantia nigra of patients with PD and Alzheimer's disease (ad). We identified 4119 differentially hydroxymethylated regions (DhMRs) and no differentially methylated regions (DMRs) in the postmortem brains of patients with PD compared with those of controls. These DhMRs were PD-specific when compared with the results of AD. Gene ontology analysis revealed that several signaling pathways, such as neurogenesis and neuronal differentiation, were significantly enriched in PD DhMRs. KEGG enrichment analysis revealed substantial alterations in multiple signaling pathways, including phospholipase D (PLD), cAMP and Rap1. In addition, using a PD Drosophila model, we found that one of the 5hmC-modulated genes, PLD1, modulated α-synuclein toxicity. Our analysis suggested that 5hmC may act as an independent epigenetic marker and contribute to the pathogenesis of PD.

Comparative analysis of atherosclerotic cardiovascular disease burden between ages 20-54 and over 55 years: insights from the Global Burden of Disease Study 2019.

BACKGROUND: To systematically analyze differences in atherosclerotic cardiovascular disease (ASCVD) burden between young and older adults. METHODS: We estimated the prevalence, mortality, and disability-adjusted life years (DALYs) of ASCVD, including ischemic heart disease (IHD), ischemic stroke (IS), and peripheral artery disease (PAD), in individuals aged 20-54 and > 55 years from 1990-2019, utilizing data from the 2019 Global Burden of Disease Study. The annual percentage changes (EAPCs) for age-specific prevalence, mortality, or DALY rates were calculated to quantify the temporal trends of ASCVD burden. We also analyzed population attribution fractions (PAF) of premature ASCVD mortality and DALYs for different risk factors and compared the burden of extremely premature, premature, and non-premature ASCVD cases based on clinical classifications. RESULTS: From 1990-2019, the global prevalence rates of IHD, IS, and PAD in the 20-54 years age group increased by 20.55% (from 694.74 to 837.49 per 100,000 population), 11.50% (from 439.48 to 490.03 per 100,000 population), and 7.38% (from 384.24 to 412.59 per 100,000 population), respectively. Conversely, the ASCVD prevalence in > 55years age group decreased. Adverse outcome burdens, including mortality and DALYs, varied among ASCVD subtypes. The decrease in the mortality/DALY burden of IHD and IS was lower in the 20-54 years group than in the > 55 years group. For PAD, DALYs among those aged 20-54 increased but decreased among those aged > 55 years. When grouped according to socio-demographic index (SDI) values, lower SDI regions exhibited a higher proportion of young ASCVD burden. The prevalence of young IHD, IS, and PAD in low SDI regions reached 20.70%, 40.05%, and 19.31% in 2019, respectively, compared with 12.14%, 16.32%, and 9.54%, respectively, in high SDI regions. Metabolic risks were the primary contributors to the ASCVD burden in both age groups. Increased susceptibility to ambient particulate matter pollution and inadequate control of high body-mass index and high fasting plasma glucose in young individuals may partially explain the differing temporal trends between young and older individuals. CONCLUSIONS: The ASCVD burden in young individuals may become a growing global health concern, especially in areas with lower socioeconomic development levels that require more effective primary prevention strategies.

Zebrafish maoc1 Attenuates Spring Viremia of Carp Virus Propagation by Promoting Autophagy-Lysosome-Dependent Degradation of Viral Phosphoprotein.

Spring viremia of carp virus (SVCV) is the causative agent of spring viremia of carp (SVC), an important infectious disease that causes high mortality in aquaculture cyprinids. How the host defends against SVCV infection and the underlying mechanisms are still elusive. In this study, we identify that a novel gene named maoc1 is induced by SVCV infection. maoc1-deficient zebrafish are more susceptible to SVCV infection, with higher virus replication and antiviral gene induction. Further assays indicate that maoc1 interacts with the P protein of SVCV to trigger P protein degradation through the autophagy-lysosomal pathway, leading to the restriction of SVCV propagation. These findings reveal a unique zebrafish defense machinery in response to SVCV infection. IMPORTANCE SVCV P protein plays an essential role in the virus replication and viral immune evasion process. Here, we identify maoc1 as a novel SVCV-inducible gene and demonstrate its antiviral capacity through attenuating SVCV replication, by directly binding to P protein and mediating its degradation via the autophagy-lysosomal pathway. Therefore, this study not only reveals an essential role of maoc1 in fighting against SVCV infection but also demonstrates an unusual host defense mechanism in response to invading viruses.

A comprehensive comparison of supervised and unsupervised methods for cell type identification in single-cell RNA-seq.

The cell type identification is among the most important tasks in single-cell RNA-sequencing (scRNA-seq) analysis. Many in silico methods have been developed and can be roughly categorized as either supervised or unsupervised. In this study, we investigated the performances of 8 supervised and 10 unsupervised cell type identification methods using 14 public scRNA-seq datasets of different tissues, sequencing protocols and species. We investigated the impacts of a number of factors, including total amount of cells, number of cell types, sequencing depth, batch effects, reference bias, cell population imbalance, unknown/novel cell type, and computational efficiency and scalability. Instead of merely comparing individual methods, we focused on factors' impacts on the general category of supervised and unsupervised methods. We found that in most scenarios, the supervised methods outperformed the unsupervised methods, except for the identification of unknown cell types. This is particularly true when the supervised methods use a reference dataset with high informational sufficiency, low complexity and high similarity to the query dataset. However, such outperformance could be undermined by some undesired dataset properties investigated in this study, which lead to uninformative and biased reference datasets. In these scenarios, unsupervised methods could be comparable to supervised methods. Our study not only explained the cell typing methods' behaviors under different experimental settings but also provided a general guideline for the choice of method according to the scientific goal and dataset properties. Finally, our evaluation workflow is implemented as a modularized R pipeline that allows future evaluation of new methods. Availability: All the source codes are available at https://github.com/xsun28/scRNAIdent.

Redox proteomics of PANC-1 cells reveals the significance of HIF-1 signaling protein oxidation in pancreatic ductal adenocarcinoma pathogenesis.

BACKGROUND: Protein cysteine oxidation is substantially involved in various biological and pathogenic processes, but its implications in pancreatic cancer development remains poorly understood. METHODS AND RESULTS: In this study, we performed a global characterization of protein oxidation targets in PDAC cells through iodoTMT-based quantitative proteomics, which identified over 4300 oxidized cysteine sites in more than 2100 proteins in HPDE6c7 and PANC-1 cells. Among them, 1715 cysteine residues were shown to be differentially oxidized between HPDE6c7 and PANC-1 cells. Also, charged amino acids including aspartate, glutamate and lysine were significantly overrepresented in flanking sequences of oxidized cysteines. Differentially oxidized proteins in PANC-1 cells were enriched in multiple cancer-related biological processes and signaling pathways. Specifically, the HIF-1 signaling proteins exhibited significant oxidation alterations in PANC-1 cells, and the reduced PHD2 oxidation in human PDAC tissues was correlated with lower survival time in pancreatic cancer patients. CONCLUSION: These investigations provided new insights into protein oxidation-regulated signaling and biological processes during PDAC pathogenesis, which might be further explored for pancreatic cancer diagnosis and treatment.

The Phenotypic and Genotypic Spectrum of CSF1R-Related Disorder in China.

BACKGROUND: Colony-stimulating factor 1 receptor (CSF1R)-related disorder (CRD) is a rare autosomal dominant disease. The clinical and genetic characteristics of Chinese patients have not been elucidated. OBJECTIVE: The objective of the study is to clarify the core features and influence factors of CRD patients in China. METHODS: Clinical and genetic-related data of CRD patients in China were collected. Mini-Mental State Examination (MMSE), Montreal Cognitive Assessment (MoCA), and Sundal MRI Severity Score were evaluated. Whole exome sequencing was used to analyze the CSF1R mutation status. Patients were compared between different sexes, mutation types, or mutation locations. RESULTS: A total of 103 patients were included, with a male-to-female ratio of 1:1.51. The average age of onset was (40.75 ± 8.58). Cognitive impairment (85.1%, 86/101) and parkinsonism (76.2%, 77/101) were the main clinical symptoms. The most common imaging feature was bilateral asymmetric white matter changes (100.0%). A total of 66 CSF1R gene mutants (22 novel mutations) were found, and 15 of 92 probands carried c.2381 T > C/p.I794T (16.30%). The MMSE and MoCA scores (17.0 [9.0], 11.90 ± 7.16) of female patients were significantly lower than those of male patients (23.0 [10.0], 16.36 ± 7.89), and the white matter severity score (20.19 ± 8.47) of female patients was significantly higher than that of male patients (16.00 ± 7.62). There is no statistical difference in age of onset between male and female patients. CONCLUSIONS: The core manifestations of Chinese CRD patients are progressive cognitive decline, parkinsonism, and bilateral asymmetric white matter changes. Compared to men, women have more severe cognitive impairment and imaging changes. c.2381 T > C/p.I794T is a hotspot mutation in Chinese patients. © 2024 International Parkinson and Movement Disorder Society.