Single-cell Analysis Research Articles

Single-cell analysis of anti-BCMA CAR T cell therapy in patients with central nervous system autoimmunity.

Chimeric antigen receptor (CAR) T cell immunotherapy for the treatment of neurological autoimmune diseases is promising, but CAR T cell kinetics and immune alterations after treatment are poorly understood. Here, we performed single-cell multi-omics sequencing of paired cerebrospinal fluid (CSF) and blood samples from patients with neuromyelitis optica spectrum disorder (NMOSD) treated with anti-B cell maturation antigen (BCMA) CAR T cells. Proliferating cytotoxic-like CD8+ CAR T cell clones were identified as the main effectors in autoimmunity. Anti-BCMA CAR T cells with enhanced features of chemotaxis efficiently crossed the blood-CSF barrier, eliminated plasmablasts and plasma cells in the CSF, and suppressed neuroinflammation. The CD44-expressing early memory phenotype in infusion products was potentially associated with CAR T cell persistence in autoimmunity. Moreover, CAR T cells from patients with NMOSD displayed distinctive features of suppressed cytotoxicity compared with those from hematological malignancies. Thus, we provide mechanistic insights into CAR T cell function in patients with neurological autoimmune disease.

Single-cell analysis defines highly specific leukemia-induced neutrophils and links MMP8 expression to recruitment of tumor associated neutrophils during FGFR1 driven leukemogenesis

BackgroundLeukemias driven by activated, chimeric FGFR1 kinases typically progress to AML which have poor prognosis. Mouse models of this syndrome allow detailed analysis of cellular and molecular changes occurring during leukemogenesis. We have used these models to determine the effects of leukemia development on the immune cell composition in the leukemia microenvironment during leukemia development and progression.MethodsSingle cell RNA sequencing (scRNA-Seq) was used to characterize leukemia associated neutrophils and define gene expression changes in these cells during leukemia progression.ResultsscRNA-Seq revealed six distinct subgroups of neutrophils based on their specific differential gene expression. In response to leukemia development, there is a dramatic increase in only two of the neutrophil subgroups. These two subgroups show specific gene expression signatures consistent with neutrophil precursors which give rise to immature polymorphonuclear myeloid-derived suppressor cells (PMN-MDSCs). Analysis of gene expression in these precursor cells identified pathways that were specifically upregulated, the most pronounced of which involved matrix metalloproteinases Mmp8 and Mmp9, during leukemia progression. Pharmacological inhibition of MMPs using Ilomastat preferentially restricted in vitro migration of neutrophils from leukemic mice and led to a significantly improved survival in vivo, accompanied by impaired PMN-MDSC recruitment. As a result, levels of T-cells were proportionally increased. In clinically annotated TCGA databases, MMP8 was shown to act as an independent indicator for poor prognosis and correlated with higher neutrophil infiltration and poor pan-cancer prognosis.ConclusionWe have defined specific leukemia responsive neutrophil subgroups based on their unique gene expression profile, which appear to be the precursors of neutrophils specifically associated with leukemia progression. An important event during development of these neutrophils is upregulation MMP genes which facilitated mobilization of these precursors from the BM in response to cancer progression, suggesting a possible therapeutic approach to suppress the development of immune tolerance.

Cancer-associated fibroblast-derived gene signature discriminates distinct prognoses by integrated single-cell and bulk RNA-seq analyses in breast cancer.

Cancer-associated fibroblasts (CAFs) are one of the most predominant cellular subpopulations in the tumor stroma and play an integral role in cancer occurrence and progression. However, the prognostic role of CAFs in breast cancer remains poorly understood. We identified a number of CAF-related biomarkers in breast cancer by combining single-cell and bulk RNA-seq analyses. Based on univariate Cox regression as well as Least Absolute Shrinkage and Selection Operator (LASSO) regression analysis, a novel CAF-associated prognostic model was developed. Breast cancer patients were grouped according to the median risk score and further analyzed for outcome, clinical characteristic, pathway activity, genomic feature, immune landscape, and drug sensitivity. A total of 341 CAF-related biomarkers were identified from single-cell and bulk RNA-seq analyses. We eventually screened eight candidate prognostic genes, including CERCAM, EMP1, SDC1, PRKG1, XG, TNN, WLS, and PDLIM4, and constructed the novel CAF-related prognostic model. Grouped by the median risk score, high-risk patients showed a significantly worse prognosis and exhibited distinct pathway activities such as uncontrolled cell cycle progression, angiogenesis, and activation of glycolysis. In addition, the combined risk score and tumor mutation burden significantly improved the ability to predict patient prognosis. Importantly, patients in the high-risk group had a higher infiltration of M2 macrophages and a lower infiltration of CD8+ T cells and activated NK cells. Finally, we calculated the IC50 for a range of anticancer drugs and personalized the treatment regimen for each patient. Integrating single-cell and bulk RNA-seq analyses, we identified a list of compositive CAF-associated biomarkers and developed a novel CAF-related prognostic model for breast cancer. This robust CAF-derived gene signature acts as an excellent predictor of patient outcomes and treatment responses in breast cancer.

Nature of epigenetic aging from a single-cell perspective.

Age-related changes in DNA methylation (DNAm) form the basis of the most robust predictors of age-epigenetic clocks-but a clear mechanistic understanding of exactly which aspects of aging are quantified by these clocks is lacking. Here, to clarify the nature of epigenetic aging, we juxtapose the dynamics of tissue and single-cell DNAm in mice. We compare these changes during early development with those observed during adult aging in mice, and corroborate our analyses with a single-cell RNA sequencing analysis within the same multiomics dataset. We show that epigenetic aging involves co-regulated changes as well as a major stochastic component, and this is consistent with transcriptional patterns. We further support the finding of stochastic epigenetic aging by direct tissue and single-cell DNAm analyses and modeling of aging DNAm trajectories with a stochastic process akin to radiocarbon decay. Finally, we describe a single-cell algorithm for the identification of co-regulated and stochastic CpG clusters showing consistent transcriptomic coordination patterns. Together, our analyses increase our understanding of the basis of epigenetic clocks and highlight potential opportunities for targeting aging and evaluating longevity interventions.

Resolvin E1 improves efferocytosis and rescues severe aplastic anemia in mice

Severe aplastic anemia (SAA) is a rare, fatal disease characterized by severe cytopenias and loss of hematopoietic stem cells (HSCs). Immune-mediated destruction and inflammation are known drivers of SAA, however, the underlying mechanisms driving persistent inflammation are unknown. Current treatments for SAA rely on immunosuppressive therapies or HSC transplantation, however, these treatments are not always effective. Using an established mouse model of SAA, we observed a significant increase in apoptotic cells within the bone marrow (BM) and impaired efferocytosis in SAA mice, relative to radiation controls. Single-cell transcriptomic analysis revealed heterogeneity among BM monocytes and unique populations emerged during SAA characterized by increased inflammatory signatures and significantly increased expression of Sirpa and Cd47. CD47, a “don’t eat me” signal, was increased on both live and apoptotic BM cells, concurrent with markedly increased expression of signal regulatory protein alpha (SIRPα) on monocytes. Functionally, SIRPα blockade improved cell clearance and reduced accumulation of CD47-positive apoptotic cells. Lipidomic analysis revealed a reduction in the precursors of specialized pro-resolving lipid mediators (SPMs) and increased prostaglandins in the BM during SAA, indicative of impaired inflammation resolution. Specifically, 18-HEPE, a precursor of E-series resolvins, was significantly reduced in SAA-induced mice relative to radiation controls. Treatment of SAA mice with Resolvin E1 (RvE1) improved efferocytic function, BM cellularity, platelet output, and survival. Our data suggest that impaired efferocytosis and inflammation resolution contributes to SAA progression and demonstrate that SPMs, such as RvE1, offer new and/or complementary treatments for SAA that do not rely on immune suppression.

Single-cell transcriptomic analyses reveal distinct immune cell contributions to epithelial barrier dysfunction in checkpoint inhibitor colitis.

Immune checkpoint inhibitor (ICI) therapy has revolutionized oncology, but treatments are limited by immune-related adverse events, including checkpoint inhibitor colitis (irColitis). Little is understood about the pathogenic mechanisms driving irColitis, which does not readily occur in model organisms, such as mice. To define molecular drivers of irColitis, we used single-cell multi-omics to profile approximately 300,000 cells from the colon mucosa and blood of 13 patients with cancer who developed irColitis (nine on anti-PD-1 or anti-CTLA-4 monotherapy and four on dual ICI therapy; most patients had skin or lung cancer), eight controls on ICI therapy and eight healthy controls. Patients with irColitis showed expanded mucosal Tregs, ITGAEHi CD8 tissue-resident memory T cells expressing CXCL13 and Th17 gene programs and recirculating ITGB2Hi CD8 T cells. Cytotoxic GNLYHi CD4 T cells, recirculating ITGB2Hi CD8 T cells and endothelial cells expressing hypoxia gene programs were further expanded in colitis associated with anti-PD-1/CTLA-4 therapy compared to anti-PD-1 therapy. Luminal epithelial cells in patients with irColitis expressed PCSK9, PD-L1 and interferon-induced signatures associated with apoptosis, increased cell turnover and malabsorption. Together, these data suggest roles for circulating T cells and epithelial-immune crosstalk critical to PD-1/CTLA-4-dependent tolerance and barrier function and identify potential therapeutic targets for irColitis.

Single-cell BCR and transcriptome analysis reveals peripheral immune signatures in patients with thyroid-associated ophthalmopathy.

Thyroid-associated ophthalmopathy (TAO) is the most prevalent orbital disease in adults caused by an autoimmune disorder, which can lead to disfigurement and vision impairment. Developing effective treatments for this condition presents challenges due to our limited understanding of its underlying immune aberrations. In this study, we profiled the immune components in the peripheral blood of patients with TAO as well as healthy individuals, utilizing single-cell RNA sequencing and B-cell receptor repertoires (BCR) analysis. We observed a significant reduction in the proportions of regulatory B cells (Bregs) and type 2 conventional dendritic cells (DCs) in patients with TAO during the active phase. Conversely, there was a significant increase in the proportion of type 1 DCs. Further analysis of cell differentiation trajectory revealed potential impairment in the transition of B cells towards Breg phenotype during the active phase of TAO. Besides, the activation process of TAO appeared to involve inflammation and immune dysfunction, as indicated by the dynamic changes in the activities of key regulators. The abnormalities in the peripheral immune system, such as the reduced capacity of Bregs to suppress inflammation, were primarily driven by the enhanced interaction among Breg, DCs, and monocytes (i.e., CD22-PTPRC and BTLA-TNFRSF14). Collectively, our findings offer a comprehensive insight into the molecular regulation and cellular reconfiguration during the active phase of TAO at the single-cell level, in order to explore the pathogenesis of TAO and provide new ideas for the future treatment of TAO.

Loss of PBX1 function in Leydig cells causes testicular dysgenesis and male sterility

Leydig cells are essential components of testicular interstitial tissue and serve as a primary source of androgen in males. A functional deficiency in Leydig cells often causes severe reproductive disorders; however, the transcriptional programs underlying the fate decisions and steroidogenesis of these cells have not been fully defined. In this study, we report that the homeodomain transcription factor PBX1 is a master regulator of Leydig cell differentiation and testosterone production in mice. PBX1 was highly expressed in Leydig cells and peritubular myoid cells in the adult testis. Conditional deletion of Pbx1 in Leydig cells caused spermatogenic defects and complete sterility. Histological examinations revealed that Pbx1 deletion impaired testicular structure and led to disorganization of the seminiferous tubules. Single-cell RNA-seq analysis revealed that loss of Pbx1 function affected the fate decisions of progenitor Leydig cells and altered the transcription of genes associated with testosterone synthesis in the adult testis. Pbx1 directly regulates the transcription of genes that play important roles in steroidogenesis (Prlr, Nr2f2 and Nedd4). Further analysis demonstrated that deletion of Pbx1 leads to a significant decrease in testosterone levels, accompanied by increases in pregnenolone, androstenedione and luteinizing hormone. Collectively, our data revealed that PBX1 is indispensable for maintaining Leydig cell function. These findings provide insights into testicular dysgenesis and the regulation of hormone secretion in Leydig cells.Graphical abstract

Knocking down of Xkr8 enhances chemotherapy efficacy through modulating tumor immune microenvironment

Scramblase Xk-related protein 8 (Xkr8) regulates the externalization of phosphatidylserine (PS) during apoptosis and holds a pivotal role in fostering tumor immunosuppression. Targeting Xkr8 in conjunction with chemotherapy demonstrated a novel avenue for amplifying antitumor immune response and overcoming chemo-immune resistance. Here we further evaluated this strategy by using a clinically relevant orthotopic model and elucidated the mechanism through in-depth single-cell RNA sequencing (scRNA-seq). We found that Xkr8 knockdown exhibited the potential to lead to immunogenic cell death (ICD) by impeding the normal clearance of apoptotic cells. Co-delivery of Xkr8 small interference RNA (siRNA) and a prodrug conjugate of 5-fluorouracil (5-Fu) and oxoplatin (FuOXP) showed remarkable therapeutic efficacy in an orthotopic pancreatic tumor model with increased infiltration of proliferative NK cells and activated macrophages in the tumor microenvironment (TME). Single-cell trajectory analysis further unveiled that tumor infiltrating CD8+ T cells are differentiated favorably to cytotoxic over exhausted phenotype after combination treatment. Our study sheds new light on the impact of Xkr8 knockdown on TME and solidifies the rationale of combining Xkr8 knockdown with chemotherapy to treat various types of cancers.

A pathological joint–liver axis mediated by matrikine-activated CD4+ T cells

The knee joint has long been considered a closed system. The pathological effects of joint diseases on distant organs have not been investigated. Herein, our clinical data showed that post-traumatic joint damage, combined with joint bleeding (hemarthrosis), exhibits a worse liver function compared with healthy control. With mouse model, hemarthrosis induces both cartilage degeneration and remote liver damage. Next, we found that hemarthrosis induces the upregulation in ratio and differentiation towards Th17 cells of CD4+ T cells in peripheral blood and spleen. Deletion of CD4+ T cells reverses hemarthrosis-induced liver damage. Degeneration of cartilage matrix induced by hemarthrosis upregulates serological type II collagen (COL II), which activates CD4+ T cells. Systemic application of a COL II antibody blocks the activation. Furthermore, bulk RNAseq and single-cell qPCR analysis revealed that the cartilage Akt pathway is inhibited by blood treatment. Intra-articular application of Akt activator blocks the cartilage degeneration and thus protects against the liver impairment in mouse and pig models. Taken together, our study revealed a pathological joint–liver axis mediated by matrikine-activated CD4+ T cells, which refreshes the organ-crosstalk axis and provides a new treatment target for hemarthrosis-related disease.Intra-articular bleeding induces cartilage degradation through down-reulation of cartilage Akt pathway. During this process, the soluble COL II released from the damaged cartilage can activate peripheral CD4+ T cells, differention into Th17 cells and secretion of IL-17, which consequently induces liver impairment. Intra-articular application of sc79 (inhibitor of Akt pathway) can prevent the cartilage damage as well as its peripheral influences.

Unraveling the immunogenic cell death pathways in gastric adenocarcinoma: A multi-omics study.

Gastric cancer (GC) is a prevalent malignant tumor of the gastrointestinal (GI) system. However, the lack of reliable biomarkers has made its diagnosis, prognosis, and treatment challenging. Immunogenic cell death (ICD) is a type of programmed cell death that is strongly related to the immune system. However, its function in GC requires further investigation. We used multi-omics and multi-angle approaches to comprehensively explore the prognostic features of ICD in patients with stomach adenocarcinoma (STAD). At the single-cell level, we screened genes associated with ICD at the transcriptome level, selected prognostic genes related to ICD using weighted gene co-expression network analysis (WGCNA) and machine learning, and constructed a prognostic model. In addition, we constructed nomograms that incorporated pertinent clinical features and provided effective tools for prognostic prediction in clinical settings. We also investigated the sensitivity of the risk subgroups to both immunotherapy and drugs. Finally, in addition to quantitative real-time polymerase chain reaction, immunofluorescence was used to validate the expression of ICD-linked genes. Based on single-cell and transcriptome WGCNA analyses, we identified 34 ICD-related genes, of which 11 were related to prognosis. We established a prognostic model using the least absolute shrinkage and selection operator (LASSO) algorithm and identified dissimilarities in overall survival (OS) and progression-free survival (PFS) in risk subgroups. The nomograms associated with the ICD-related signature (ICDRS) demonstrated a good predictive value for clinical applications. Moreover, we detected changes in the tumor microenvironment (TME), including biological functions, mutation landscapes, and immune cell infiltration, between the high- and low-risk groups. We constructed an ICD-related prognostic model that incorporated features related to cell death. This model can serve as a useful tool for predicting the prognosis of GC, targeted prevention, and personalized medicine.

Native N-glycome profiling of single cells and ng-level blood isolates using label-free capillary electrophoresis-mass spectrometry

The development of reliable single-cell dispensers and substantial sensitivity improvement in mass spectrometry made proteomic profiling of individual cells achievable. Yet, there are no established methods for single-cell glycome analysis due to the inability to amplify glycans and sample losses associated with sample processing and glycan labeling. In this work, we present an integrated platform coupling online in-capillary sample processing with high-sensitivity label-free capillary electrophoresis-mass spectrometry for N-glycan profiling of single mammalian cells. Direct and unbiased quantitative characterization of single-cell surface N-glycomes are demonstrated for HeLa and U87 cells, with the detection of up to 100 N-glycans per single cell. Interestingly, N-glycome alterations are unequivocally detected at the single-cell level in HeLa and U87 cells stimulated with lipopolysaccharide. The developed workflow is also applied to the profiling of ng-level amounts (5–500 ng) of blood-derived protein, extracellular vesicle, and total plasma isolates, resulting in over 170, 220, and 370 quantitated N-glycans, respectively.

Cross-species single-cell analysis uncovers the immunopathological mechanisms associated with IgA nephropathy progression.

IgA nephropathy (IgAN) represents the main cause of renal failure, while the precise pathogenetic mechanisms have not been fully determined. Herein, we conducted a cross-species single-cell survey on human IgAN and mouse and rat IgAN models to explore the pathogenic programs. Cross-species single-cell RNA sequencing (scRNA-Seq) revealed that the IgAN mesangial cells (MCs) expressed high levels of inflammatory signatures CXCL12, CCL2, CSF1, and IL-34 and specifically interacted with IgAN macrophages via the CXCL12/CXCR4, CSF1/IL-34/CSF1 receptor, and integrin subunit alpha X/integrin subunit alpha M/complement C3 (C3) axes. IgAN macrophages expressed high levels of CXCR4, PDGFB, triggering receptor expressed on myeloid cells 2, TNF, and C3, and the trajectory analysis suggested that these cells derived from the differentiation of infiltrating blood monocytes. Additionally, protein profiling of 21 progression and 28 nonprogression IgAN samples revealed that proteins CXCL12, C3, mannose receptor C-type 1, and CD163 were negatively correlated with estimated glomerular filtration rate (eGFR) value and poor prognosis (30% eGFR as composite end point). Last, a functional experiment revealed that specific blockade of the Cxcl12/Cxcr4 pathway substantially attenuated the glomerulus and tubule inflammatory injury, fibrosis, and renal function decline in the mouse IgAN model. This study provides insights into IgAN progression and may aid in the refinement of IgAN diagnosis and the optimization of treatment strategies.

IL-2RG as a possible immunotherapeutic target in CRC predicting poor prognosis and regulated by miR-7-5p and miR-26b-5p

Despite advances in treatment strategies, colorectal cancer (CRC) continues to cause significant morbidity and mortality, with mounting evidence a close link between immune system dysfunctions issued. Interleukin-2 receptor gamma (IL-2RG) plays a pivotal role as a common subunit receptor in the IL-2 family cytokines and activates the JAK-STAT pathway. This study delves into the role of Interleukin-2 receptor gamma (IL-2RG) within the tumor microenvironment and investigates potential microRNAs (miRNAs) that directly inhibit IL-2RG, aiming to discern their impact on CRC clinical outcomes. Bioinformatics analysis revealed a significant upregulation of IL-2RG mRNA in TCGA-COAD samples and showed strong correlations with the infiltration of various lymphocytes. Single-cell analysis corroborated these findings, highlighting IL-2RG expression in critical immune cell subsets. To explore miRNA involvement in IL-2RG dysregulation, mRNA was isolated from the tumor tissues and lymphocytes of 258 CRC patients and 30 healthy controls, and IL-2RG was cloned into the pcDNA3.1/CT-GFP-TOPO vector. Human embryonic kidney cell lines (HEK-293T) were transfected with this construct. Our research involved a comprehensive analysis of miRPathDB, miRWalk, and Targetscan databases to identify the miRNAs associated with the 3′ UTR of human IL-2RG. The human microRNA (miRNA) molecules, hsa-miR-7-5p and hsa-miR-26b-5p, have been identified as potent suppressors of IL-2RG expression in CRC patients. Specifically, the downregulation of hsa-miR-7-5p and hsa-miR-26b-5p has been shown to result in the upregulation of IL-2RG mRNA expression in these patients. Prognostic evaluation of IL-2RG, hsa-miR-7-5p, and hsa-miR-26b-5p, using TCGA-COAD data and patient samples, established that higher IL-2RG expression and lower expression of both miRNAs were associated with poorer outcomes. Additionally, this study identified several long non-coding RNAs (LncRNAs), such as ZFAS1, SOX21-AS1, SNHG11, SNHG16, SNHG1, DLX6-AS1, GAS5, SNHG6, and MALAT1, which may act as competing endogenous RNA molecules for IL2RG by sequestering shared hsa-miR-7-5p and hsa-miR-26b-5p. In summary, this investigation underscores the potential utility of IL-2RG, hsa-miR-7-5p, and hsa-miR-26b-5p as serum and tissue biomarkers for predicting CRC patient prognosis while also offering promise as targets for immunotherapy in CRC management.Graphical

Single-cell analyses reveal the metabolic heterogeneity and plasticity of the tumor microenvironment during head and neck squamous cell carcinoma progression.

Metabolic reprogramming is a hallmark of cancer. In addition to metabolic alterations in the tumor cells, multiple other metabolically active cell types in the tumor microenvironment (TME) contribute to the emergence of a tumor-specific metabolic milieu. Here, we defined the metabolic landscape of the TME during progression of head and neck squamous cell carcinoma (HNSCC) by performing single-cell RNA sequencing (scRNA-seq) on 26 human patient specimens, including normal tissue, pre-cancerous lesions, early-stage cancer, advanced-stage cancer, lymph node metastases, and recurrent tumors. The analysis revealed substantial heterogeneity at the transcriptional, developmental, metabolic, and functional levels in different cell types. SPP1+ macrophages were identified as a pro-tumor and pro-metastatic macrophage subtype with high fructose and mannose metabolism, which was further substantiated by integrative analysis and validation experiments. An inhibitor of fructose metabolism reduced the proportion of SPP1+ macrophages, reshaped the immunosuppressive TME, and suppressed tumor growth. In conclusion, this work delineated the metabolic landscape of HNSCC at a single-cell resolution and identified fructose metabolism as a key metabolic feature of a pro-tumor macrophage subpopulation.

Salivary Gland Tissue Recombination Can Modify Cell Fate.

Although mesenchyme is essential for inducing the epithelium of ectodermal organs, its precise role in organ-specific epithelial fate determination remains poorly understood. To elucidate the roles of tissue interactions in cellular differentiation, we performed single-cell RNA sequencing and imaging analyses on recombined tissues, where mesenchyme and epithelium were switched ex vivo between two types of embryonic mouse salivary glands: the parotid gland (a serous gland) and the submandibular gland (a predominantly mucous gland). We found partial induction of molecules that define gland-specific acinar and myoepithelial cells in recombined salivary epithelium. The parotid epithelium recombined with submandibular mesenchyme began to express mucous acinar genes not intrinsic to the parotid gland. While myoepithelial cells do not normally line parotid acini, newly induced myoepithelial cells densely populated recombined parotid acini. However, mucous acinar and myoepithelial markers continued to be expressed in submandibular epithelial cells recombined with parotid mesenchyme. Consequently, some epithelial cells appeared to be plastic, such that their fate could still be modified in response to mesenchymal signaling, whereas other epithelial cells appeared to be already committed to a specific fate. We also discovered evidence for bidirectional induction: transcriptional changes were observed not only in the epithelium but also in the mesenchyme after heterotypic tissue recombination. For example, parotid epithelium induced the expression of muscle-related genes in submandibular fibroblasts that began to mimic parotid fibroblast gene expression. These studies provide the first comprehensive unbiased molecular characterization of tissue recombination approaches exploring the regulation of cell fate.

TauFisher predicts circadian time from a single sample of bulk and single-cell pseudobulk transcriptomic data

As the circadian clock regulates fundamental biological processes, disrupted clocks are often observed in patients and diseased tissues. Determining the circadian time of the patient or the tissue of focus is essential in circadian medicine and research. Here we present tauFisher, a computational pipeline that accurately predicts circadian time from a single transcriptomic sample by finding correlations between rhythmic genes within the sample. We demonstrate tauFisher’s performance in adding timestamps to both bulk and single-cell transcriptomic samples collected from multiple tissue types and experimental settings. Application of tauFisher at a cell-type level in a single-cell RNAseq dataset collected from mouse dermal skin implies that greater circadian phase heterogeneity may explain the dampened rhythm of collective core clock gene expression in dermal immune cells compared to dermal fibroblasts. Given its robustness and generalizability across assay platforms, experimental setups, and tissue types, as well as its potential application in single-cell RNAseq data analysis, tauFisher is a promising tool that facilitates circadian medicine and research.

Integrative analysis of single-cell and bulk RNA sequencing unveils a machine learning-based pan-cancer major histocompatibility complex-related signature for predicting immunotherapy efficacy

Major histocompatibility complex (MHC) could serve as a potential biomarker for tumor immunotherapy, however, it is not yet known whether MHC could distinguish potential beneficiaries. Single‐cell RNA sequencing datasets derived from patients with immunotherapy were collected to elucidate the association between MHC and immunotherapy response. A novel MHCsig was developed and validated using large‐scale pan‐cancer data, including The Cancer Genome Atlas and immunotherapy cohorts. The therapeutic value of MHCsig was further explored using 17 CRISPR/Cas9 datasets. MHC-related genes were associated with drug resistance and MHCsig was significantly and positively associated with immunotherapy response and total mutational burden. Remarkably, MHCsig significantly enriched 6% top‐ranked genes, which were potential therapeutic targets. Moreover, we generated Hub-MHCsig, which was associated with survival and disease-special survival of pan-cancer, especially low-grade glioma. This result was also confirmed in cell lines and in our own clinical cohort. Later low-grade glioma-related Hub-MHCsig was established and the regulatory network was constructed. We provided conclusive clinical evidence regarding the association between MHCsig and immunotherapy response. We developed MHCsig, which could effectively predict the benefits of immunotherapy for multiple tumors. Further exploration of MHCsig revealed some potential therapeutic targets and regulatory networks.

Boosting humoral and cellular immunity with enhanced STING activation by hierarchical mesoporous metal-organic framework adjuvants

Vaccination is essential for preventing and controlling infectious diseases, along with reducing mortality. Developing safe and versatile adjuvants to enhance humoral and cellular immune responses to vaccines remains a key challenge in vaccine development. Here, we designed hierarchical mesoporous MOF-801 (HM801) using a Cocamidopropyl betaine (CAPB) and a Pluronics F127 in an aqueous phase system. Meanwhile, we synthesized a novel SARS-CoV-2 nanovaccine (R@M@HM801) with a high loading capacity for both the STING agonist (MSA-2) and the Delta receptor binding domain (Delta-RBD) antigen. R@M@HM801 enhanced MSA-2 and RBD utilization and effectively co-delivered MSA-2 and RBD antigens to antigen-presenting cells in the draining lymph nodes, thereby promoting the activation of both T and B cells. Lymphocyte single-cell analysis showed that R@M@HM801 stimulated robust CD11b+CD4+ T cells, CXCR5+CD4+ T follicular helper (Tfh), and durable CD4+CD44+CD62L−, CD8+CD44+CD62L− effector memory T cell (TEM) immune responses, and promoted the proliferative activation of CD26+ B cells in vivo. Meanwhile, R@M@HM801 induced stronger specific antibodies and neutralization of pseudovirus against Delta compared to the RBD + MAS-2 and RBD + MAS-2 + Alum vaccines. Our study demonstrated the efficacy of a hierarchical mesoporous HM801 and its potential immune activation mechanism in enhancing adaptive immune responses against viruses and other diseases.

E4BP4 in macrophages induces an anti-inflammatory phenotype that ameliorates the severity of colitis

Macrophages are versatile cells of the innate immune system that work by altering their pro- or anti-inflammatory features. Their dysregulation leads to inflammatory disorders such as inflammatory bowel disease. We show that macrophage-specific upregulation of the clock output gene and transcription factor E4BP4 reduces the severity of colitis in mice. RNA-sequencing and single-cell analyses of macrophages revealed that increased expression of E4BP4 leads to an overall increase in expression of anti-inflammatory genes including Il4ra with a concomitant reduction in pro-inflammatory gene expression. In contrast, knockout of E4BP4 in macrophages leads to increased proinflammatory gene expression and decreased expression of anti-inflammatory genes. ChIP-seq and ATAC-seq analyses further identified Il4ra as a target of E4BP4, which drives anti-inflammatory polarization in macrophages. Together, these results reveal a critical role for E4BP4 in regulating macrophage inflammatory phenotypes and resolving inflammatory bowel diseases.