MAST: a flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data

Single-Cell Transcriptomic Analysis of Normal and Malignant B Cells

In the past decade, single-cell (sc) transcriptomics has overcome the limitations of bulk analysis by measuring gene expression in individual cells, not just a population average. This can identify diverse cell types and states within a sample with high resolution, even without prior purification. Various technologies exist, each with its own capture, barcoding, and library preparation methods. This chapter focuses on the analysis of normal and malignant mature B cells using the 10× Genomics 5' sc-gene expression in parallel with B cell immune repertoire profiling. By integrating the gene expression data from similar cells, the complete transcriptome for each population can be reconstructed, while the identification of the expressed immunoglobulin genes allows investigating clonotype evolution and the detection of tumor clones that share the same clonally rearranged B cell receptor sequence. Researchers are guided through both the experimental protocols and data analysis with a comprehensive, step-by-step walkthrough of how to use some of the more popular single-cell software tools.

CTLA-4 expressing innate lymphoid cells modulate mucosal homeostasis in a microbiota dependent manner

The maintenance of intestinal homeostasis is a fundamental process critical for organismal integrity. Sitting at the interface of the gut microbiome and mucosal immunity, adaptive and innate lymphoid populations regulate the balance between commensal micro-organisms and pathogens. Checkpoint inhibitors, particularly those targeting the CTLA-4 pathway, disrupt this fine balance and can lead to inflammatory bowel disease and immune checkpoint colitis. Here, we show that CTLA-4 is expressed by innate lymphoid cells and that its expression is regulated by ILC subset-specific cytokine cues in a microbiota-dependent manner. Genetic deletion or antibody blockade of CTLA-4 in multiple in vivo models of colitis demonstrates that this pathway plays a key role in intestinal homeostasis. Lastly, we have found that this observation is conserved in human IBD. We propose that this population of CTLA-4-positive ILC may serve as an important target for the treatment of idiopathic and iatrogenic intestinal inflammation.

ITGA5+ synovial fibroblasts orchestrate proinflammatory niche formation by remodelling the local immune microenvironment in rheumatoid arthritis

OBJECTIVES

To investigate the phenotypic heterogeneity of tissue-resident synovial fibroblasts and their role in inflammatory response in rheumatoid arthritis (RA).

METHODS

We used single-cell and spatial transcriptomics to profile synovial cells and spatial gene expressions of synovial tissues to identify phenotypic changes in patients with osteoarthritis, RA in sustained remission and active state. Immunohistology, multiplex immunofluorescence and flow cytometry were used to identify synovial fibroblasts subsets. Deconvolution methods further validated our findings in two cohorts (PEAC and R4RA) with treatment response. Cell coculture was used to access the potential cell-cell interactions. Adoptive transfer of synovial cells in collagen-induced arthritis (CIA) mice and bulk RNA sequencing of synovial joints further validate the cellular functions.

RESULTS

We identified a novel tissue-remodelling CD45-CD31-PDPN+ITGA5+ synovial fibroblast population with unique transcriptome of POSTN, COL3A1, CCL5 and TGFB1, and enriched in immunoregulatory pathways. This subset was upregulated in active and lympho-myeloid type of RA, associated with an increased risk of multidrug resistance. Transforming growth factor (TGF)-β1 might participate in the differentiation of this subset. Moreover, ITGA5+ synovial fibroblasts might occur in early stage of inflammation and induce the differentiation of CXCL13hiPD-1hi peripheral helper T cells (TPHs) from naïve CD4+ T cells, by secreting TGF-β1. Intra-articular injection of ITGA5+ synovial fibroblasts exacerbates RA development and upregulates TPHs in CIA mice.

CONCLUSIONS

We demonstrate that ITGA5+ synovial fibroblasts might regulate the RA progression by inducing the differentiation of CXCL13hiPD-1hi TPHs and remodelling the proinflammatory microenvironments. Therapeutic modulation of this subpopulation could therefore be a potential treatment strategy for RA.

Integrative genomic analyses identify neuroblastoma risk genes involved in neuronal differentiation

Genome-Wide Association Studies (GWAS) have been decisive in elucidating the genetic predisposition of neuroblastoma (NB). The majority of genetic variants identified in GWAS are found in non-coding regions, suggesting that they can be causative of pathogenic dysregulations of gene expression. Nonetheless, pinpointing the potential causal genes within implicated genetic loci remains a major challenge. In this study, we integrated NB GWAS and expression Quantitative Trait Loci (eQTL) data from adrenal gland to identify candidate genes impacting NB susceptibility. We found that ZMYM1, CBL, GSKIP and WDR81 expression was dysregulated by NB predisposing variants. We further investigated the functional role of the identified genes through computational analysis of RNA sequencing (RNA-seq) data from single-cell and whole-tissue samples of NB, neural crest, and adrenal gland tissues, as well as through in vitro differentiation assays in NB cell cultures. Our results indicate that dysregulation of ZMYM1, CBL, GSKIP, WDR81 may lead to malignant transformation by affecting early and late stages of normal program of neuronal differentiation. Our findings enhance the understanding of how specific genes contribute to NB pathogenesis by highlighting their influence on neuronal differentiation and emphasizing the impact of genetic risk variants on the regulation of genes involved in critical biological processes.

A longitudinal single-cell atlas of anti-tumour necrosis factor treatment in inflammatory bowel disease

Precision medicine in immune-mediated inflammatory diseases (IMIDs) requires a cellular understanding of treatment response. We describe a therapeutic atlas for Crohn's disease (CD) and ulcerative colitis (UC) following adalimumab, an anti-tumour necrosis factor (anti-TNF) treatment. We generated ~1 million single-cell transcriptomes, organised into 109 cell states, from 216 gut biopsies (41 subjects), revealing disease-specific differences. A systems biology-spatial analysis identified granuloma signatures in CD and interferon (IFN)-response signatures localising to T cell aggregates and epithelial damage in CD and UC. Pretreatment differences in epithelial and myeloid compartments were associated with remission outcomes in both diseases. Longitudinal comparisons demonstrated disease progression in nonremission: myeloid and T cell perturbations in CD and increased multi-cellular IFN signalling in UC. IFN signalling was also observed in rheumatoid arthritis (RA) synovium with a lymphoid pathotype. Our therapeutic atlas represents the largest cellular census of perturbation with the most common biologic treatment, anti-TNF, across multiple inflammatory diseases.

Immunotherapy response induces divergent tertiary lymphoid structure morphologies in hepatocellular carcinoma

Tertiary lymphoid structures (TLS) are associated with improved response in solid tumors treated with immune checkpoint blockade, but understanding of the prognostic and predictive value of TLS and the circumstances of their resolution is incomplete. Here we show that in hepatocellular carcinoma treated with neoadjuvant immunotherapy, high intratumoral TLS density at the time of surgery is associated with pathologic response and improved relapse-free survival. In areas of tumor regression, we identify a noncanonical involuted morphology of TLS marked by dispersion of the B cell follicle, persistence of a T cell zone enriched for T cell-mature dendritic cell interactions and increased expression of T cell memory markers. Collectively, these data suggest that TLS can serve as both a prognostic and predictive marker of response to immunotherapy in hepatocellular carcinoma and that late-stage TLS may support T cell memory formation after elimination of a viable tumor.

Exploring the molecular underpinnings of macrosomia in gestational diabetes mellitus: The role of EGFR signaling and placental syncytiotrophoblast

Macrosomia, which is frequently associated with gestational diabetes mellitus (GDM), is linked to maternal glycemic control during gestation. When GDM is complicated by macrosomia (GDMM), the placenta exhibits increased mass, underscoring its role as a critical nexus for maternal-fetal nutrient exchange. Despite this recognized correlation, the underlying mechanisms propelling placental hypertrophy have remained elusive. Our study leveraged single-cell RNA transcriptome sequencing of GDMM placental tissues to pinpoint the specific syncytiotrophoblast (SCT) subsets that regulate placental dimensions. Notably, we observed pronounced upregulation of the epidermal growth factor receptor (EGFR) and its corresponding ligands, with a particular emphasis on the autoregulatory cascade involving the glycoprotein hormone alpha subunit (CGA), EGFR, and the transcription factor GATA binding protein 2 (GATA2), as well as perturbations in hormonal homeostasis within the SCT. Furthermore, our cell interaction analysis revealed an enhanced interplay between myeloid cells and SCT3, augmenting the EGFR signaling pathway. These molecular exchanges underscore the pivotal role of the placental immune microenvironment in the etiology of macrosomia, shedding light on the pathophysiology of GDMM and paving the way for novel therapeutic approaches.

Method of moments framework for differential expression analysis of single-cell RNA sequencing data

Differential expression analysis of single-cell RNA sequencing (scRNA-seq) data is central for characterizing how experimental factors affect the distribution of gene expression. However, distinguishing between biological and technical sources of cell-cell variability and assessing the statistical significance of quantitative comparisons between cell groups remain challenging. We introduce Memento, a tool for robust and efficient differential analysis of mean expression, variability, and gene correlation from scRNA-seq data, scalable to millions of cells and thousands of samples. We applied Memento to 70,000 tracheal epithelial cells to identify interferon-responsive genes, 160,000 CRISPR-Cas9 perturbed T cells to reconstruct gene-regulatory networks, 1.2 million peripheral blood mononuclear cells (PBMCs) to map cell-type-specific quantitative trait loci (QTLs), and the 50-million-cell CELLxGENE Discover corpus to compare arbitrary cell groups. In all cases, Memento identified more significant and reproducible differences in mean expression compared with existing methods. It also identified differences in variability and gene correlation that suggest distinct transcriptional regulation mechanisms imparted by perturbations.

Single cell analysis of Idh mutant growth plates identifies cell populations responsible for longitudinal bone growth and enchondroma formation

Enchondromas are a common tumor in bone that can occur as multiple lesions in enchondromatosis, which is associated with deformity of the affected bone. These lesions harbor somatic mutations in IDH and driving expression of a mutant Idh1 in Col2 expressing cells in mice causes an enchondromatosis phenotype. Here we compared growth plates from E18.5 mice expressing a mutant Idh1 with control littermates using single cell RNA sequencing. Data from Col2 expressing cells were analysed using UMAP and RNA pseudo-time analyses. A unique cluster of cells was identified in the mutant growth plates that expressed genes known to be upregulated in enchondromas. There was also a cluster of cells that was underrepresented in the mutant growth plates that expressed genes known to be important in longitudinal bone growth. Immunofluorescence showed that the genes from the unique cluster identified in the mutant growth plates were expressed in multiple growth plate anatomic zones, and pseudo-time analysis also suggested these cells could arise from multiple growth plate chondrocyte subpopulations. This data supports the notion that a subpopulation of chondrocytes become enchondromas at the expense of contributing to longitudinal growth.

Progressive plasticity during colorectal cancer metastasis

As cancers progress, they become increasingly aggressive-metastatic tumours are less responsive to first-line therapies than primary tumours, they acquire resistance to successive therapies and eventually cause death1,2. Mutations are largely conserved between primary and metastatic tumours from the same patients, suggesting that non-genetic phenotypic plasticity has a major role in cancer progression and therapy resistance3-5. However, we lack an understanding of metastatic cell states and the mechanisms by which they transition. Here, in a cohort of biospecimen trios from same-patient normal colon, primary and metastatic colorectal cancer, we show that, although primary tumours largely adopt LGR5+ intestinal stem-like states, metastases display progressive plasticity. Cancer cells lose intestinal cell identities and reprogram into a highly conserved fetal progenitor state before undergoing non-canonical differentiation into divergent squamous and neuroendocrine-like states, a process that is exacerbated in metastasis and by chemotherapy and is associated with poor patient survival. Using matched patient-derived organoids, we demonstrate that metastatic cells exhibit greater cell-autonomous multilineage differentiation potential in response to microenvironment cues compared with their intestinal lineage-restricted primary tumour counterparts. We identify PROX1 as a repressor of non-intestinal lineage in the fetal progenitor state, and show that downregulation of PROX1 licenses non-canonical reprogramming.

Single-cell transcriptomics of rectal organoids from individuals with perianal fistulizing Crohn's disease reveals patient-specific signatures

Perianal fistulizing Crohn's disease (CD) is a severe gastrointestinal disorder causing extensive mucosal damage with limited treatment options. Severe manifestations of the disease appear at higher rates in non-Europeans but the genetic and cellular mechanisms driving the disease phenotypes remain poorly understood. Herein, we tested whether pathologic determinants in the epithelial stem cell compartment could be detected at the transcript level in rectal organoids derived from a diverse patient population. Rectal organoid and mucosal cells from endoscopic biopsies of each patient having perianal fistulizing CD or no disease controls were prepared for and sequenced at the single cell level. After cell type annotations based on expressed marker genes, samples were analyzed by principal components, for differential transcript expression, cell type proportions, and pathway enrichment. After QC, we produced 77,044 rectal organoid cells (n = 13 patients; 8 CD, 5 controls) with high quality sequences that identified 10 distinct epithelial subtypes, that we compared to 141,367 mucosal epithelial cells (n = 29 patients; 18 CD, 11 controls). Consistent with mucosal epithelial cells, rectal organoids prominently displayed disease signatures represented by the stem and transit amplifying regions of the rectal crypt, including alterations in transcriptional signatures of metabolic, epigenetic, and proliferating pathways. Organoids also retained their gender- and ancestral-specific gene expression signatures. However, they lacked many of the inflammatory signatures observed in epithelial cells from diseased mucosa. Perianal CD patient derived rectal organoids reflect gene expression signatures related to disease, gender, and ancestry, suggesting they harbor inherent properties amenable to further patient-specific, disease-related experimentation.

Heterogeneous enhancer states orchestrate β cell responses to metabolic stress

Obesity-induced β cell dysfunction contributes to the onset of type 2 diabetes. Nevertheless, elucidating epigenetic mechanisms underlying islet dysfunction at single cell level remains challenging. Here we profile single-nuclei RNA along with enhancer marks H3K4me1 or H3K27ac in islets from lean or obese mice. Our study identifies distinct gene signatures and enhancer states correlating with β cell dysfunction trajectory. Intriguingly, while many metabolic stress-induced genes exhibit concordant changes in both H3K4me1 and H3K27ac at their enhancers, expression changes of specific subsets are solely attributable to either H3K4me1 or H3K27ac dynamics. Remarkably, a subset of H3K4me1+H3K27ac- primed enhancers prevalent in lean β cells and occupied by FoxA2 are largely absent after metabolic stress. Lastly, cell-cell communication analysis identified the nerve growth factor (NGF) as protective paracrine signaling for β cells through repressing ER stress. In summary, our findings define the heterogeneous enhancer responses to metabolic challenges in individual β cells.

Theoretical framework for the difference of two negative binomial distributions and its application in comparative analysis of sequencing data

High-throughput sequencing (HTS) technologies have been instrumental in investigating biological questions at the bulk and single-cell levels. Comparative analysis of two HTS data sets often relies on testing the statistical significance for the difference of two negative binomial distributions (DOTNB). Although negative binomial distributions are well studied, the theoretical results for DOTNB remain largely unexplored. Here, we derive basic analytical results for DOTNB and examine its asymptotic properties. As a state-of-the-art application of DOTNB, we introduce DEGage, a computational method for detecting differentially expressed genes (DEGs) in scRNA-seq data. DEGage calculates the mean of the sample-wise differences of gene expression levels as the test statistic and determines significant differential expression by computing the P-value with DOTNB. Extensive validation using simulated and real scRNA-seq data sets demonstrates that DEGage outperforms five popular DEG analysis tools: DEGseq2, DEsingle, edgeR, Monocle3, and scDD. DEGage is robust against high dropout levels and exhibits superior sensitivity when applied to balanced and imbalanced data sets, even with small sample sizes. We utilize DEGage to analyze prostate cancer scRNA-seq data sets and identify marker genes for 17 cell types. Furthermore, we apply DEGage to scRNA-seq data sets of mouse neurons with and without fear memory and reveal eight potential memory-related genes overlooked in previous analyses. The theoretical results and supporting software for DOTNB can be widely applied to comparative analyses of dispersed count data in HTS and broad research questions.

Physiological and pathogenic T cell autoreactivity converge in type 1 diabetes

Autoimmune diseases result from autoantigen-mediated activation of adaptive immunity; intriguingly, autoantigen-specific T cells are also present in healthy donors. An assessment of dynamic changes of this autoreactive repertoire in both health and disease is thus warranted. Here we investigate the physiological versus pathogenic autoreactive processes in the context of Type 1 diabetes (T1D) and one of its landmark autoantigens, glutamic acid decarboxylase 65 (GAD65). Using single cell gene expression profiling and tandem T cell receptor (TCR) sequencing, we find that GAD65-specific true naïve cells are present in both health and disease, with GAD65-specific effector and memory responses showing similar ratios in healthy donors and patients. Deeper assessment of phenotype and TCR repertoire uncover differential features in GAD65-specific TCRs, including lower clonal sizes of healthy donor-derived clonotypes in patients. We thus propose a model whereby physiological autoimmunity against GAD65 is needed during early life, and that alterations of these physiological autoimmune processes in predisposed individuals trigger overt Type 1 diabetes.

Spatial pattern and differential expression analysis with spatial transcriptomic data

The emergence of spatial transcriptomic technologies has opened new avenues for investigating gene activities while preserving the spatial context of tissues. Utilizing data generated by such technologies, the identification of spatially variable (SV) genes is an essential step in exploring tissue landscapes and biological processes. Particularly in typical experimental designs, such as case-control or longitudinal studies, identifying SV genes between groups is crucial for discovering significant biomarkers or developing targeted therapies for diseases. However, current methods available for analyzing spatial transcriptomic data are still in their infancy, and none of the existing methods are capable of identifying SV genes between groups. To overcome this challenge, we developed SPADE for spatial pattern and differential expression analysis to identify SV genes in spatial transcriptomic data. SPADE is based on a machine learning model of Gaussian process regression with a gene-specific Gaussian kernel, enabling the detection of SV genes both within and between groups. Through benchmarking against existing methods in extensive simulations and real data analyses, we demonstrated the preferred performance of SPADE in detecting SV genes within and between groups. The SPADE source code and documentation are publicly available at https://github.com/thecailab/SPADE.

Age, sex, and cell type-resolved hypothalamic gene expression across the pubertal transition in mice

BACKGROUND

The hypothalamus plays a central role in regulating puberty. However, our knowledge of the postnatal gene regulatory networks that control the pubertal transition in males and females is incomplete. Here, we investigate the age-, sex- and cell-type-specific gene regulation in the hypothalamus across the pubertal transition.

METHODS

We used RNA-seq to profile hypothalamic gene expression in male and female mice at five time points spanning the onset of puberty (postnatal days (PD) 12, 22, 27, 32, and 37). By combining this data with hypothalamic single nuclei RNA-seq data from pre- and postpubertal mice, we assigned gene expression changes to their most likely cell types of origin. In our colony, pubertal onset occurs earlier in male mice, allowing us to focus on genes whose expression is dynamic across ages and offset between sexes, and to explore the bases of sex effects.

RESULTS

Our age-by-sex pattern of expression enriched for biological pathways involved hormone production, neuronal activation, and glial maturation. Additionally, we inferred a robust expansion of oligodendrocytes precursor cells into mature oligodendrocytes spanning the prepubertal (PD12) to peri-pubertal (PD27) timepoints. Using spatial transcriptomic data from postpubertal mice, we observed the lateral hypothalamic area and zona incerta were the most oligodendrocyte-rich regions and that these cells expressed genes known to be involved in pubertal regulation.

CONCLUSION

Together, by incorporating multiple biological timepoints and using sex as a variable, we identified gene and cell-type changes that may participate in orchestrating the pubertal transition and provided a resource for future studies of postnatal hypothalamic gene regulation.

FORMATION OF MALIGNANT, METASTATIC SMALL CELL LUNG CANCERS THROUGH OVERPRODUCTION OF cMYC PROTEIN IN TP53 AND RB1 DEPLETED PULMONARY NEUROENDOCRINE CELLS DERIVED FROM HUMAN EMBRYONIC STEM CELLS

We recently described our initial efforts to develop a model for small cell lung cancer (SCLC) derived from human embryonic stem cells (hESCs) that were differentiated to form pulmonary neuroendocrine cells (PNECs), a putative cell of origin for neuroendocrine-positive SCLC. Although reduced expression of the tumor suppressor genes TP53 and RB1 allowed the induced PNECs to form subcutaneous growths in immune-deficient mice, the tumors did not display the aggressive characteristics of SCLC seen in human patients. Here we report that the additional, doxycycline-regulated expression of a transgene encoding wild-type or mutant cMYC protein promotes rapid growth, invasion, and metastasis of these hESC-derived cells after injection into the renal capsule. Similar to others, we find that the addition of cMYC encourages the formation of the SCLC-N subtype, marked by high levels of NEUROD1 RNA. Using paired primary and metastatic samples for RNA sequencing, we observe that the subtype of SCLC does not change upon metastatic spread and that production of NEUROD1 is maintained. We also describe histological features of these malignant, SCLC-like tumors derived from hESCs and discuss potential uses of this model in efforts to control and better understand this recalcitrant neoplasm.

Myeloid-derived β-hexosaminidase is essential for neuronal health and lysosome function: implications for Sandhoff disease

Lysosomal storage disorders (LSDs) are a large disease class involving lysosomal dysfunction, often resulting in neurodegeneration. Sandhoff disease (SD) is an LSD caused by a deficiency in the β subunit of the β-hexosaminidase enzyme ( Hexb ). Although Hexb expression in the brain is specific to microglia, SD primarily affects neurons. To understand how a microglial gene is involved in maintaining neuronal homeostasis, we demonstrated that β-hexosaminidase is secreted by microglia and integrated into the neuronal lysosomal compartment. To assess therapeutic relevance, we treated SD mice with bone marrow transplant and colony stimulating factor 1 receptor inhibition, which broadly replaced Hexb -/- microglia with Hexb -sufficient cells. This intervention reversed apoptotic gene signatures, improved behavior, restored enzymatic activity and Hexb expression, ameliorated substrate accumulation, and normalized neuronal lysosomal phenotypes. These results underscore the critical role of myeloid-derived β- hexosaminidase in neuronal lysosomal function and establish microglial replacement as a potential LSD therapy.

Post-stroke hippocampal neurogenesis is impaired by microvascular dysfunction and PI3K signaling in cerebral amyloid angiopathy

Ischemic stroke and cerebral amyloid angiopathy (CAA) pose significant challenges in an aging population, particularly in post-stroke recovery. Using the 5xFAD mouse model, we explore the relationship between CAA, ischemic stroke, and tissue recovery. We hypothesize that amyloid-beta accumulation worsens stroke outcomes by inducing blood-brain barrier (BBB) dysfunction, leading to impaired neurogenesis. Our findings show that CAA exacerbates stroke outcomes, with mice exhibiting constricted BBB microvessels, reduced cerebral blood flow, and impaired tissue recovery. Transcriptional analysis shows that endothelial cells and neural progenitor cells (NPCs) in the hippocampus exhibit differential gene expression in response to CAA and stroke, specifically targeting the phosphatidylinositol 3-kinase (PI3K) pathway. In vitro experiments with human NPCs validate these findings, showing that disruption of the CXCL12-PIK3C2A-CREB3L2 axis impairs neurogenesis. Notably, PI3K pathway activation restores neurogenesis, highlighting a potential therapeutic approach. These results suggest that CAA combined with stroke induces microvascular dysfunction and aberrant neurogenesis through this specific pathway.

PI3K/AKT signaling controls ICM maturation and proper epiblast and primitive endoderm specification in mice

The inner cell mass (ICM) of early mouse embryos is specified into epiblast (Epi) and primitive endoderm (PrE) lineages during blastocyst formation. The antagonistic transcription factors (TFs) NANOG and GATA-binding protein 6 (GATA6) in combination with fibroblast growth factor (FGF)/extracellular-signal-regulated kinase (ERK) signaling are central actors in ICM fate choice. However, what initiates the specification of ICM progenitors into Epi or PrE and whether other factors are involved in this process has not been fully understood yet. Here, we show that phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT) is constitutively active during preimplantation development. Using pharmacological inhibition, we demonstrate that PI3K/AKT enables the formation of a functional ICM capable of giving rise to both the Epi and the PrE: it maintains the expression of the TF NANOG, which specifies the Epi, and confers responsiveness to FGF4, which is essential for PrE specification. Our work thus identifies PI3K/AKT signaling as an upstream regulator controlling the molecular events required for both Epi and PrE specification.

Mapping the glial transcriptome in Huntington's disease using snRNAseq: selective disruption of glial signatures across brain regions

Huntington's disease (HD) is an autosomal dominant neurodegenerative disease with a fatal outcome. There is accumulating evidence of a prominent role of glia in the pathology of HD, and we investigated this by conducting single nuclear RNA sequencing (snRNAseq) of human post mortem brain in four differentially affected regions; caudate nucleus, frontal cortex, hippocampus and cerebellum. Across 127,205 nuclei from donors with HD and age/sex matched controls, we found heterogeneity of glia which is altered in HD. We describe prominent changes in the abundance of certain subtypes of astrocytes, microglia, oligodendrocyte precursor cells and oligodendrocytes between HD and control samples, and these differences are widespread across brain regions. Furthermore, we highlight possible mechanisms that characterise the glial contribution to HD pathology including depletion of myelinating oligodendrocytes, an oligodendrocyte-specific upregulation of the calmodulin-dependent 3',5'-cyclic nucleotide phosphodiesterase 1 A (PDE1A) and an upregulation of molecular chaperones as a cross-glial signature and a potential adaptive response to the accumulation of mutant huntingtin (mHTT). Our results support the hypothesis that glia have an important role in the pathology of HD, and show that all types of glia are affected in the disease.

Single-cell RNA sequencing reveals the process of CA19-9 production and dynamics of the immune microenvironment between CA19-9 (+) and CA19-9 (-) PDAC

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is one of the main types of malignant tumor of the digestive system, and patient prognosis is affected by difficulties in early diagnosis, poor treatment response, and a high postoperative recurrence rate. Carbohydrate antigen 19-9 (CA19-9) has been widely used as a biomarker for the diagnosis and postoperative follow-up of PDAC patients. Nevertheless, the production mechanism and potential role of CA19-9 in PDAC progression have not yet been elucidated.

METHODS

We performed single-cell RNA sequencing on six samples pathologically diagnosed as PDAC (three CA19-9-positive and three CA19-9-negative PDAC samples) and two paracarcinoma samples. We also downloaded and integrated PDAC samples (each from three CA19-9-positive and CA19-9-negative patients) from an online database. The dynamics of the proportion and potential function of each cell type were verified through immunofluorescence. Moreover, we built an in vitro coculture cellular model to confirm the potential function of CA19-9.

RESULTS

Three subtypes of cancer cells with a high ability to produce CA19-9 were identified by the markers TOP2A , AQP5 , and MUC5AC . CA19-9 production bypass was discovered on antigen-presenting cancer-associated fibroblasts (apCAFs). Importantly, the proportion of immature ficolin-1 positive (FCN1+) macrophages was high in the CA19-9-negative group, and the proportion of mature M2-like macrophages was high in the CA19-9-positive group. High proportions of these two macrophage subtypes were associated with an unfavourable clinical prognosis. Further experiments indicated that CA19-9 could facilitate the transformation of M0 macrophages into M2 macrophages in the tumor microenvironment.

CONCLUSIONS

Our study described CA19-9 production at single-cell resolution and the dynamics of the immune atlas in CA19-9-positive and CA19-9-negative PDAC. CA19-9 could promote M2 polarization of macrophage in the pancreatic tumor microenvironment.

Aging is associated with a modality-specific decline in taste

Deficits in chemosensory processing are associated with healthy aging, as well as numerous neurodegenerative disorders, including Alzheimer's disease (AD). The fruit fly, Drosophila melanogaster, is a powerful model for studying chemosensation, aging, and aging-related pathologies, yet the effects of aging and neurodegeneration on taste function remain largely unexplored. Aging impaired response to sugars, but not medium-chain fatty acids that are sensed by a shared population of neurons. Selective expression of the human amyloid beta (Aβ) peptide phenocopied the effects of aging. Functional imaging of gustatory axon terminals revealed reduced response to sugar, but not fatty acids. Axonal innervation of the fly taste center was largely intact in aged flies; however, axonal innervation was reduced upon expression of Aβ. A comparison of transcript expression within the sugar-sensing taste neurons revealed age-related changes in 66 genes. Together, these findings suggest that different mechanisms underly taste deficits in aged and AD model flies.

Single-cell transcriptomics of pediatric Burkitt lymphoma reveals intra-tumor heterogeneity and markers of therapy resistance

Burkitt lymphoma (BL) is the most frequent B-cell lymphoma in pediatric patients. While most patients are cured, a fraction of them are resistant to therapy. To investigate BL heterogeneity and the features distinguishing therapy responders (R) from non-responders (NR), we analyzed by single-cell (sc)-transcriptomics diagnostic EBV-negative BL specimens. Analysis of the non-tumor component revealed a predominance of immune cells and a small representation of fibroblasts, enriched in NR. Tumors displayed patient-specific features, as well as shared subpopulations that expressed transcripts related to cell cycle, signaling pathways and cell-of-origin signatures. Several transcripts were differentially expressed in R versus NR. The top candidate, Tropomyosin 2 (TPM2), a member of the tropomyosin actin filament binding protein family, was confirmed to be significantly higher in NR both at the transcript and protein level. Stratification of patients based on TPM2 expression at diagnosis significantly correlated with prognosis, independently of TP53 mutations. These results indicate that BL displays transcriptional heterogeneity and identify candidate biomarkers of therapy resistance.

Epigenetic and 3D genome reprogramming during the aging of human hippocampus

Age-related cognitive decline is associated with altered physiology of the hippocampus. While changes in gene expression have been observed in aging brain, the regulatory mechanisms underlying these changes remain underexplored. We generated single-nucleus gene expression, chromatin accessibility, DNA methylation, and 3D genome data from 40 human hippocampal tissues spanning adult lifespan. We observed a striking loss of astrocytes, OPC, and endothelial cells during aging, including astrocytes that play a role in regulating synapses. Microglia undergo a dramatic switch from a homeostatic state to a primed inflammatory state through DNA methylome and 3D genome reprogramming. Aged cells experience erosion of their 3D genome architecture. Our study identifies age-associated changes in cell types/states and gene regulatory features that provide insight into cognitive decline during human aging.

Tissue niche influences immune and metabolic profiles to Staphylococcus aureus biofilm infection

Infection is a devastating post-surgical complication, often requiring additional procedures and prolonged antibiotic therapy. This is especially relevant for craniotomy and prosthetic joint infections (PJI), both of which are characterized by biofilm formation on the bone or implant surface, respectively, with S. aureus representing a primary cause. The local tissue microenvironment likely has profound effects on immune attributes that can influence treatment efficacy, which becomes critical to consider when developing therapeutics for biofilm infections. However, the extent to which distinct tissue niches influence immune function during biofilm development remains relatively unknown. To address this, we compare the metabolomic, transcriptomic, and functional attributes of leukocytes in mouse models of S. aureus craniotomy and PJI complemented with patient samples from both infection modalities, which reveals profound tissue niche-dependent differences in nucleic acid, amino acid, and lipid metabolism with links to immune modulation. These signatures are both spatially and temporally distinct, differing not only between infection sites but evolving over time within a single model. Collectively, this demonstrates that biofilms elicit unique immune and metabolic responses that are heavily influenced by the local tissue microenvironment, which will likely have important implications when designing therapeutic approaches targeting these infections.

Defective ribosome assembly impairs leukemia progression in a murine model of acute myeloid leukemia

Despite an advanced understanding of disease mechanisms, the current therapeutic regimen fails to cure most patients with acute myeloid leukemia (AML). In the present study, we address the role of ribosome assembly in leukemia cell function. We apply patient datasets and murine models to demonstrate that immature leukemia cells in mixed-lineage leukemia-rearranged AML are characterized by relatively high ribosome biogenesis and protein synthesis rates. Using a model with inducible regulation of ribosomal subunit joining, we show that defective ribosome assembly extends survival in mice with AML. Single-cell RNA sequencing and proteomic analyses reveal that leukemia cell adaptation to defective ribosome assembly is associated with an increase in ribosome biogenesis and deregulation of the transcription factor landscape. Finally, we demonstrate that defective ribosome assembly shows antileukemia efficacy in p53-deficient AML. Our study unveils the critical requirement of a high protein synthesis rate for leukemia progression and highlights ribosome assembly as a therapeutic target in AML.

Accurate identification of locally aneuploid cells by incorporating cytogenetic information in single cell data analysis

Single-cell RNA sequencing is a powerful tool to investigate the cellular makeup of tumor samples. However, due to the sparse data and the complex tumor microenvironment, it can be challenging to identify neoplastic cells that play important roles in tumor growth and disease progression. This is especially relevant for blood cancers, where neoplastic cells may be highly similar to normal cells. To address this challenge, we have developed partCNV and partCNVH, two methods for rapid and accurate detection of aneuploid cells with local copy number deletion or amplification. PartCNV uses an expectation-maximization (EM) algorithm with mixtures of Poisson distributions and incorporates cytogenetic information to guide the classification. PartCNVH further improves partCNV by integrating a hidden Markov model for feature selection. We have thoroughly evaluated the performance of partCNV and partCNVH through simulation studies and real data analysis using three scRNA-seq datasets from blood cancer patients. Our results show that partCNV and partCNVH have favorable accuracy and provide more interpretable results compared to existing methods. In the real data analysis, we have identified multiple biological processes involved in the oncogenesis of myelodysplastic syndromes and acute myeloid leukemia.

Characterization of Cancer Stem Cells in Laryngeal Squamous Cell Carcinoma by Single-cell RNA Sequencing

Cancer stem cells (CSCs) constitute a pivotal element within the tumor microenvironment (TME), driving the initiation and progression of cancer. However, the identification of CSCs and their underlying molecular mechanisms in laryngeal squamous cell carcinoma (LSCC) remains a formidable challenge. Here, we employed single-cell RNA sequencing of matched primary tumor tissues, paracancerous tissues, and local lymph nodes from three LSCC patients to comprehensively characterize the CSCs in LSCC. Two distinct clusters of stem cells originating from epithelial populations were delineated and verified as CSCs and normal stem cells (NSCs), respectively. CSCs were abundant in the paracancerous tissues compared to those in the tumor tissues. CSCs showed high expression of stem cell marker genes such as PROM1, ALDH1A1, and SOX4, and increased the activity of tumor-related hypoxia, Wnt/β-catenin, and Notch signaling pathways. We then explored the intricate crosstalk between CSCs and the TME cells and identified targets within the TME that related with CSCs. We also found eight marker genes of CSCs that were correlated significantly with the prognosis of LSCC patients. Furthermore, bioinformatics analyses showed that drugs such as erlotinib, OSI-027, and ibrutinib selectively targeted the CSC-specifically expressed genes. In conclusion, our results represent the first comprehensive characterization of CSC properties in LSCC at the single-cell level.

Distinct tumor architectures and microenvironments for the initiation of breast cancer metastasis in the brain

Brain metastasis, a serious complication of cancer, hinges on the initial survival, microenvironment adaptation, and outgrowth of disseminated cancer cells. To understand the early stages of brain colonization, we investigated two prevalent sources of cerebral relapse, triple-negative (TNBC) and HER2+ (HER2BC) breast cancers. Using mouse models and human tissue samples, we found that these tumor types colonize the brain, with a preference for distinctive tumor architectures, stromal interfaces, and autocrine programs. TNBC models tend to form perivascular sheaths with diffusive contact with astrocytes and microglia. In contrast, HER2BC models tend to form compact spheroids driven by autonomous tenascin C production, segregating stromal cells to the periphery. Single-cell transcriptomics of the tumor microenvironment revealed that these architectures evoke differential Alzheimer's disease-associated microglia (DAM) responses and engagement of the GAS6 receptor AXL. The spatial features of the two modes of brain colonization have relevance for leveraging the stroma to treat brain metastasis.

IMMUNE AND MOLECULAR CORRELATES OF RESPONSE TO IMMUNOTHERAPY REVEALED BY BRAIN-METASTATIC MELANOMA MODELS

Despite the promising results of immune checkpoint blockade (ICB) therapy, outcomes for patients with brain metastasis (BrM) remain poor. Identifying resistance mechanisms has been hindered by limited access to patient samples and relevant preclinical models. Here, we developed two mouse melanoma BrM models that recapitulate the disparate responses to ICB seen in patients. We demonstrate that these models capture the cellular and molecular complexity of human disease and reveal key factors shaping the tumor microenvironment and influencing ICB response. BR1-responsive tumor cells express inflammatory programs that polarize microglia into reactive states, eliciting robust T cell recruitment. In contrast, BR3-resistant melanoma cells are enriched in neurological programs and exploit tolerance mechanisms to maintain microglia homeostasis and limit T cell infiltration. In humans, BR1 and BR3 expression signatures correlate positively or negatively with T cell infiltration and BrM patient outcomes, respectively. Our study provides clinically relevant models and uncovers mechanistic insights into BrM ICB responses, offering potential biomarkers and therapeutic targets to improve therapy efficacy.

Meningeal lymphatic function promotes oligodendrocyte survival and brain myelination

The precise neurophysiological changes prompted by meningeal lymphatic dysfunction remain unclear. Here, we showed that inducing meningeal lymphatic vessel ablation in adult mice led to gene expression changes in glial cells, followed by reductions in mature oligodendrocyte numbers and specific lipid species in the brain. These phenomena were accompanied by altered meningeal adaptive immunity and brain myeloid cell activation. During brain remyelination, meningeal lymphatic dysfunction provoked a state of immunosuppression that contributed to delayed spontaneous oligodendrocyte replenishment and axonal loss. The deficiencies in mature oligodendrocytes and neuroinflammation due to impaired meningeal lymphatic function were solely recapitulated in immunocompetent mice. Patients diagnosed with multiple sclerosis presented reduced vascular endothelial growth factor C in the cerebrospinal fluid, particularly shortly after clinical relapses, possibly indicative of poor meningeal lymphatic function. These data demonstrate that meningeal lymphatics regulate oligodendrocyte function and brain myelination, which might have implications for human demyelinating diseases.

CD4+ T cells re-wire granuloma cellularity and regulatory networks to promote immunomodulation following Mtb reinfection

Immunological priming-in the context of either prior infection or vaccination-elicits protective responses against subsequent Mycobacterium tuberculosis (Mtb) infection. However, the changes that occur in the lung cellular milieu post-primary Mtb infection and their contributions to protection upon reinfection remain poorly understood. Using clinical and microbiological endpoints in a non-human primate reinfection model, we demonstrated that prior Mtb infection elicited a long-lasting protective response against subsequent Mtb exposure and was CD4+ T cell dependent. By analyzing data from primary infection, reinfection, and reinfection-CD4+ T cell-depleted granulomas, we found that the presence of CD4+ T cells during reinfection resulted in a less inflammatory lung milieu characterized by reprogrammed CD8+ T cells, reduced neutrophilia, and blunted type 1 immune signaling among myeloid cells. These results open avenues for developing vaccines and therapeutics that not only target lymphocytes but also modulate innate immune cells to limit tuberculosis (TB) disease.

Precursor central memory versus effector cell fate and naïve CD4+ T cell heterogeneity

Upon antigenic stimulation, naïve CD4+ T cells can give rise to phenotypically distinct effector T helper cells and long-lived memory T cells. We computationally reconstructed the in vivo trajectory of CD4+ T cell differentiation during a type I inflammatory immune response and identified two distinct differentiation paths for effector and precursor central memory T cells arising directly from naïve CD4+ T cells. Unexpectedly, our studies revealed heterogeneity among naïve CD4+ T cells, which are typically considered homogeneous save for their diverse T cell receptor usage. Specifically, a previously unappreciated population of naïve CD4+ T cells sensing environmental type I IFN exhibited distinct activation thresholds, suggesting that naïve CD4+ T cell differentiation potential may be influenced by environmental cues. This population was expanded in human viral infection and type I IFN response-lined autoimmunity. Understanding the relevance of naïve T cell heterogeneity to beneficial and maladaptive T cell responses may have therapeutic implications for adoptive T cell therapies in cancer immunotherapy and vaccination.

Multi-omic human pancreatic islet endoplasmic reticulum and cytokine stress response mapping provides type 2 diabetes genetic insights

Endoplasmic reticulum (ER) and inflammatory stress responses contribute to islet dysfunction in type 2 diabetes (T2D). Comprehensive genomic understanding of these human islet stress responses and whether T2D-associated genetic variants modulate them is lacking. Here, comparative transcriptome and epigenome analyses of human islets exposed ex vivo to these stressors revealed 30% of expressed genes and 14% of islet cis-regulatory elements (CREs) as stress responsive, modulated largely in an ER- or cytokine-specific fashion. T2D variants overlapped 86 stress-responsive CREs, including 21 induced by ER stress. We linked the rs6917676-T T2D risk allele to increased islet ER-stress-responsive CRE accessibility and allele-specific β cell nuclear factor binding. MAP3K5, the ER-stress-responsive putative rs6917676 T2D effector gene, promoted stress-induced β cell apoptosis. Supporting its pro-diabetogenic role, MAP3K5 expression correlated inversely with human islet β cell abundance and was elevated in T2D β cells. This study provides genome-wide insights into human islet stress responses and context-specific T2D variant effects.

Synergic activity of FGFR2 and MEK inhibitors in the treatment of FGFR2-amplified cancers of unknown primary

Patients with cancer of unknown primary (CUP) carry the double burden of an aggressive disease and reduced access to therapies. Experimental models are pivotal for CUP biology investigation and drug testing. We derived two CUP cell lines (CUP#55 and #96) and corresponding patient-derived xenografts (PDXs), from ascites tumor cells. CUP cell lines and PDXs underwent histological, immune-phenotypical, molecular, and genomic characterization confirming the features of the original tumor. The tissue-of-origin prediction was obtained from the tumor microRNA expression profile and confirmed by single-cell transcriptomics. Genomic testing and fluorescence in situ hybridization analysis identified FGFR2 gene amplification in both models, in the form of homogeneously staining region (HSR) in CUP#55 and double minutes in CUP#96. FGFR2 was recognized as the main oncogenic driver and therapeutic target. FGFR2-targeting drug BGJ398 (infigratinib) in combination with the MEK inhibitor trametinib proved to be synergic and exceptionally active, both in vitro and in vivo. The effects of the combined treatment by single-cell gene expression analysis revealed a remarkable plasticity of tumor cells and the greater sensitivity of cells with epithelial phenotype. This study brings personalized therapy closer to CUP patients and provides the rationale for FGFR2 and MEK targeting in metastatic tumors with FGFR2 pathway activation.

Unveiling the dynamics and molecular landscape of a rare chronic lymphocytic leukemia subpopulation driving refractoriness: insights from single-cell RNA sequencing

Early identification of resistant cancer cells is currently a major challenge, as their expansion leads to refractoriness. To capture the dynamics of these cells, we made a comprehensive analysis of disease progression and treatment response in a chronic lymphocytic leukemia (CLL) patient using a combination of single-cell and bulk genomic methods. At diagnosis, the patient presented with unfavorable genetic markers, including notch receptor 1 (NOTCH1) mutation and loss(11q). The initial and subsequent treatment lines did not lead to a durable response and the patient developed refractory disease. Refractory CLL cells featured substantial dysregulation in B-cell phenotypic markers such as human leukocyte antigen (HLA) genes, immunoglobulin (IG) genes, CD19 molecule (CD19), membrane spanning 4-domains A1 (MS4A1; previously known as CD20), CD79a molecule (CD79A) and paired box 5 (PAX5), indicating B-cell de-differentiation and disease transformation. We described the clonal evolution and characterized in detail two cell populations that emerged during the refractory disease phase, differing in the presence of high genomic complexity. In addition, we successfully tracked the cells with high genomic complexity back to the time before treatment, where they formed a rare subpopulation. We have confirmed that single-cell RNA sequencing enables the characterization of refractory cells and the monitoring of their development over time.

Non-endothelial expression of endomucin in the mouse and human choroid

Endomucin (EMCN) is a 261 amino acid transmembrane glycoprotein that is highly expressed by venous and capillary endothelial cells where it plays a role in VEGF-mediated angiogenesis and regulation of immune cell recruitment. However, it is better known as a histological marker, where it has become widespread due to the commercial availability of high-quality antibodies that work under a wide range of conditions and in many tissues. The specificity of EMCN staining has been well-validated in retinal vessels, but while it has been used extensively as a marker in other tissues of the eye, including the choroid, the pattern of expression has not been described in detail. Here, in addition to endothelial expression in the choriocapillaris and deeper vascular layers, we characterize a population of EMCN-positive perivascular cells in the mouse choroid that did not co-localize with cells expressing other endothelial markers such as PECAM1 or PODXL. To confirm that these cells represented a new population of EMCN-expressing stromal cells, we then performed single cell RNA sequencing in choroids from adult wild-type mice. Analysis of this new dataset confirmed that, in addition to endothelial cells, Emcn mRNA expression was present in choroidal pericytes and a subset of fibroblasts, but not vascular smooth muscle cells. Besides Emcn, no known endothelial gene expression was detected in these cell populations, confirming that they did not represent endothelial-stromal doublets, a common technical artifact in single cell RNA seq datasets. Instead, choroidal Emcn-expressing fibroblasts exhibited high levels of chemokine and interferon signaling genes, while Emcn-negative fibroblasts were enriched in genes encoding extracellular matrix proteins. Emcn expressing fibroblasts were also detected in published datasets from mouse brain and human choroid, suggesting that stromal Emcn expression was not unique to the choroid and was evolutionarily conserved. Together, these findings highlight unique fibroblast and pericyte populations in the choroid and provide new context for the role of EMCN in the eye.

Molecular programs guiding arealization of descending cortical pathways

Layer 5 extratelencephalic (ET) neurons are present across neocortical areas and send axons to multiple subcortical targets1-6. Two cardinal subtypes exist7,8: (1) Slco2a1-expressing neurons (ETdist), which predominate in the motor cortex and project distally to the pons, medulla and spinal cord; and (2) Nprs1- or Hpgd-expressing neurons (ETprox), which predominate in the visual cortex and project more proximally to the pons and thalamus. An understanding of how area-specific ETdist and ETprox emerge during development is important because they are critical for fine motor skills and are susceptible to spinal cord injury and amyotrophic lateral sclerosis9-12. Here, using cross-areal mapping of axonal projections in the mouse neocortex, we identify the subtype-specific developmental dynamics of ET neurons. Whereas subsets of ETprox emerge by pruning of ETdist axons, others emerge de novo. We outline corresponding subtype-specific developmental transcriptional programs using single-nucleus sequencing. Leveraging these findings, we use postnatal in vivo knockdown of subtype-specific transcription factors to reprogram ET neuron connectivity towards more proximal targets. Together, these results show the functional transcriptional programs driving ET neuron diversity and uncover cell subtype-specific gene regulatory networks that can be manipulated to direct target specificity in motor corticofugal pathways.

Single-cell transcriptomic analysis of flowering regulation and vernalization in Chinese cabbage shoot apex

In Chinese cabbage development the interplay between shoot apex activity and vernalization is pivotal for flowering timing. The intricate relationship between various cell types in the shoot apex meristem and their roles in regulating flowering gene expression in Chinese cabbage is not yet fully understood. A thorough analysis of single-cell types in the Chinese cabbage shoot apex and their influence on flowering genes and vernalization is essential for deeper insight. Our study first established a single-cell transcriptomic atlas of Chinese cabbage after 25 days of non-vernalization. Analyzing 19 602 single cells, we differentiated them into 15 distinct cell clusters using established marker genes. We found that key genes in shoot apex development and flowering were primarily present in shoot meristematic cells (SMCs), companion cells (CCs), and mesophyll cells (MCs). MADS-box protein FLOWERING LOCUS C 2 (BrFLC2), a gene suppressing flowering, was observed in CCs, mirroring patterns found in Arabidopsis. By mapping developmental trajectories of SMCs, CCs, and MCs, we elucidated the evolutionary pathways of crucial genes in shoot apex development and flowering. The creation of a single-cell transcriptional atlas of the Chinese cabbage shoot apex under vernalization revealed distinct alterations in the expression of known flowering genes, such as VERNALIZATION INSENSITIVE 3 (VIN3), VERNALIZATION 1 (VRN1), VERNALIZATION 2 (VRN2), BrFLC, and FLOWERING LOCUS T (FT), which varied by cell type. Our study underscores the transformative impact of single-cell RNA sequencing (scRNA-seq) for unraveling the complex differentiation and vernalization processes in the Chinese cabbage shoot apex. These insights are pivotal for enhancing breeding strategies and cultivation management of this vital vegetable.

A Prime-Boost Vaccination Approach Induces Lung Resident Memory CD8+ T Cells Derived from Central Memory T Cells That Prevent Tumor Lung Metastasis

Memory T cells play a key role in immune protection against cancer. Vaccine-induced tissue-resident memory T (TRM) cells in the lung have been shown to protect against lung metastasis. Identifying the source of lung TRM cells can help to improve strategies, preventing tumor metastasis. Here, we found that a prime-boost vaccination approach using intramuscular DNA vaccine priming, followed by intranasal live-attenuated influenza-vectored vaccine (LAIV) boosting induced higher frequencies of lung CD8+ TRM cells compared with other vaccination regimens. Vaccine-induced lung CD8+ TRM cells, but not circulating memory T cells, conferred significant protection against metastatic melanoma and mesothelioma. Central memory T (TCM) cells induced by the DNA vaccination were major precursors of lung TRM cells established after the intranasal LAIV boost. Single-cell RNA sequencing analysis indicated that transcriptional reprogramming of TCM cells for differentiation into TRM cells in the lungs started as early as day 2 post the LAIV boost. Intranasal LAIV altered the mucosal microenvironment to recruit TCM cells via CXCR3-dependent chemotaxis and induced CD8+ TRM-associated transcriptional programs. These results identified TCM cells as the source of vaccine-induced CD8+ TRM cells that protect against lung metastasis. Significance: Prime-boost vaccination shapes the mucosal microenvironment and reprograms central memory T cells to generate lung resident memory T cells that protect against lung metastasis, providing insights for the optimization of vaccine strategies.

The stress response regulator HSF1 modulates natural killer cell anti-tumour immunity

Diverse cellular insults converge on activation of the heat shock factor 1 (HSF1), which regulates the proteotoxic stress response to maintain protein homoeostasis. HSF1 regulates numerous gene programmes beyond the proteotoxic stress response in a cell-type- and context-specific manner to promote malignancy. However, the role(s) of HSF1 in immune populations of the tumour microenvironment remain elusive. Here, we leverage an in vivo model of HSF1 activation and single-cell transcriptomic tumour profiling to show that augmented HSF1 activity in natural killer (NK) cells impairs cytotoxicity, cytokine production and subsequent anti-tumour immunity. Mechanistically, HSF1 directly binds and regulates the expression of key mediators of NK cell effector function. This work demonstrates that HSF1 regulates the immune response under the stress conditions of the tumour microenvironment. These findings have important implications for enhancing the efficacy of adoptive NK cell therapies and for designing combinatorial strategies including modulators of NK cell-mediated tumour killing.

Building, benchmarking, and exploring perturbative maps of transcriptional and morphological data

The continued scaling of genetic perturbation technologies combined with high-dimensional assays such as cellular microscopy and RNA-sequencing has enabled genome-scale reverse-genetics experiments that go beyond single-endpoint measurements of growth or lethality. Datasets emerging from these experiments can be combined to construct perturbative "maps of biology", in which readouts from various manipulations (e.g., CRISPR-Cas9 knockout, CRISPRi knockdown, compound treatment) are placed in unified, relatable embedding spaces allowing for the generation of genome-scale sets of pairwise comparisons. These maps of biology capture known biological relationships and uncover new associations which can be used for downstream discovery tasks. Construction of these maps involves many technical choices in both experimental and computational protocols, motivating the design of benchmark procedures to evaluate map quality in a systematic, unbiased manner. Here, we (1) establish a standardized terminology for the steps involved in perturbative map building, (2) introduce key classes of benchmarks to assess the quality of such maps, (3) construct 18 maps from four genome-scale datasets employing different cell types, perturbation technologies, and data readout modalities, (4) generate benchmark metrics for the constructed maps and investigate the reasons for performance variations, and (5) demonstrate utility of these maps to discover new biology by suggesting roles for two largely uncharacterized genes.

ScRNAbox: empowering single-cell RNA sequencing on high performance computing systems

BACKGROUND

Single-cell RNA sequencing (scRNAseq) offers powerful insights, but the surge in sample sizes demands more computational power than local workstations can provide. Consequently, high-performance computing (HPC) systems have become imperative. Existing web apps designed to analyze scRNAseq data lack scalability and integration capabilities, while analysis packages demand coding expertise, hindering accessibility.

RESULTS

In response, we introduce scRNAbox, an innovative scRNAseq analysis pipeline meticulously crafted for HPC systems. This end-to-end solution, executed via the SLURM workload manager, efficiently processes raw data from standard and Hashtag samples. It incorporates quality control filtering, sample integration, clustering, cluster annotation tools, and facilitates cell type-specific differential gene expression analysis between two groups. We demonstrate the application of scRNAbox by analyzing two publicly available datasets.

CONCLUSION

ScRNAbox is a comprehensive end-to-end pipeline designed to streamline the processing and analysis of scRNAseq data. By responding to the pressing demand for a user-friendly, HPC solution, scRNAbox bridges the gap between the growing computational demands of scRNAseq analysis and the coding expertise required to meet them.

Tumor editing suppresses innate and adaptive antitumor immunity and is reversed by inhibiting DNA methylation

Cancer cells edit gene expression to evade immunosurveillance. However, genome-wide studies of gene editing during early tumorigenesis are lacking. Here we used single-cell RNA sequencing in a breast cancer genetically engineered mouse model (GEMM) to identify edited genes without bias. Late tumors repressed antitumor immunity genes, reducing infiltrating immune cells and tumor-immune cell communications. Innate immune genes, especially interferon-stimulated genes, dominated the list of downregulated tumor genes, while genes that regulate cell-intrinsic malignancy were mostly unedited. Naive and activated CD8+ T cells in early tumors were replaced with exhausted or precursor-exhausted cells in late tumors. Repression of immune genes was reversed by inhibiting DNA methylation using low-dose decitabine, which suppressed tumor growth and restored immune control, increasing the number, functionality and memory of tumor-infiltrating lymphocytes and reducing the number of myeloid suppressor cells. Decitabine induced important interferon, pyroptosis and necroptosis genes, inflammatory cell death and immune control in GEMM and implanted breast and melanoma tumors.

Single-cell transcriptomic profiling of human pancreatic islets reveals genes responsive to glucose exposure over 24 h

AIMS/HYPOTHESIS

Disruption of pancreatic islet function and glucose homeostasis can lead to the development of sustained hyperglycaemia, beta cell glucotoxicity and subsequently type 2 diabetes. In this study, we explored the effects of in vitro hyperglycaemic conditions on human pancreatic islet gene expression across 24 h in six pancreatic cell types: alpha; beta; gamma; delta; ductal; and acinar. We hypothesised that genes associated with hyperglycaemic conditions may be relevant to the onset and progression of diabetes.

METHODS

We exposed human pancreatic islets from two donors to low (2.8 mmol/l) and high (15.0 mmol/l) glucose concentrations over 24 h in vitro. To assess the transcriptome, we performed single-cell RNA-seq (scRNA-seq) at seven time points. We modelled time as both a discrete and continuous variable to determine momentary and longitudinal changes in transcription associated with islet time in culture or glucose exposure. Additionally, we integrated genomic features and genetic summary statistics to nominate candidate effector genes. For three of these genes, we functionally characterised the effect on insulin production and secretion using CRISPR interference to knock down gene expression in EndoC-βH1 cells, followed by a glucose-stimulated insulin secretion assay.

RESULTS

In the discrete time models, we identified 1344 genes associated with time and 668 genes associated with glucose exposure across all cell types and time points. In the continuous time models, we identified 1311 genes associated with time, 345 genes associated with glucose exposure and 418 genes associated with interaction effects between time and glucose across all cell types. By integrating these expression profiles with summary statistics from genetic association studies, we identified 2449 candidate effector genes for type 2 diabetes, HbA1c, random blood glucose and fasting blood glucose. Of these candidate effector genes, we showed that three (ERO1B, HNRNPA2B1 and RHOBTB3) exhibited an effect on glucose-stimulated insulin production and secretion in EndoC-βH1 cells.

CONCLUSIONS/INTERPRETATION

The findings of our study provide an in-depth characterisation of the 24 h transcriptomic response of human pancreatic islets to glucose exposure at a single-cell resolution. By integrating differentially expressed genes with genetic signals for type 2 diabetes and glucose-related traits, we provide insights into the molecular mechanisms underlying glucose homeostasis. Finally, we provide functional evidence to support the role of three candidate effector genes in insulin secretion and production.

DATA AVAILABILITY

The scRNA-seq data from the 24 h glucose exposure experiment performed in this study are available in the database of Genotypes and Phenotypes (dbGap; https://www.ncbi.nlm.nih.gov/gap/ ) with accession no. phs001188.v3.p1. Study metadata and summary statistics for the differential expression, gene set enrichment and candidate effector gene prediction analyses are available in the Zenodo data repository ( https://zenodo.org/ ) under accession number 11123248. The code used in this study is publicly available at https://github.com/CollinsLabBioComp/publication-islet_glucose_timecourse .

Novel microglial transcriptional signatures promote social and cognitive deficits following repeated social defeat

Chronic stress is associated with anxiety and cognitive impairment. Repeated social defeat (RSD) in mice induces anxiety-like behavior driven by microglia and the recruitment of inflammatory monocytes to the brain. Nonetheless, it is unclear how microglia communicate with other cells to modulate the physiological and behavioral responses to stress. Using single-cell (sc)RNAseq, we identify novel, to the best of our knowledge, stress-associated microglia in the hippocampus defined by RNA profiles of cytokine/chemokine signaling, cellular stress, and phagocytosis. Microglia depletion with a CSF1R antagonist (PLX5622) attenuates the stress-associated profile of leukocytes, endothelia, and astrocytes. Furthermore, RSD-induced social withdrawal and cognitive impairment are microglia-dependent, but social avoidance is microglia-independent. Furthermore, single-nuclei (sn)RNAseq shows robust responses to RSD in hippocampal neurons that are both microglia-dependent and independent. Notably, stress-induced CREB, oxytocin, and glutamatergic signaling in neurons are microglia-dependent. Collectively, these stress-associated microglia influence transcriptional profiles in the hippocampus related to social and cognitive deficits.

Hypoxia induces ROS-resistant memory upon reoxygenation in vivo promoting metastasis in part via MUC1-C

Hypoxia occurs in 90% of solid tumors and is associated with metastasis and mortality. Breast cancer cells that experience intratumoral hypoxia are 5x more likely to develop lung metastasis in animal models. Using spatial transcriptomics, we determine that hypoxic cells localized in more oxygenated tumor regions (termed 'post-hypoxic') retain expression of hypoxia-inducible and NF-kB-regulated genes, even in the oxygen-rich bloodstream. This cellular response is reproduced in vitro under chronic hypoxic conditions followed by reoxygenation. A subset of genes remains increased in reoxygenated cells. MUC1/MUC1-C is upregulated by both HIF-1α and NF-kB-p65 during chronic hypoxia. Abrogating MUC1 decreases the expression of superoxide dismutase enzymes, causing reactive oxygen species (ROS) production and cell death. A hypoxia-dependent genetic deletion of MUC1, or MUC1-C inhibition by GO-203, increases ROS levels in circulating tumor cells (CTCs), reducing the extent of metastasis. High MUC1 expression in tumor biopsies is associated with recurrence, and MUC1+ CTCs have lower ROS levels than MUC1- CTCs in patient-derived xenograft models. This study demonstrates that therapeutically targeting MUC1-C reduces hypoxia-driven metastasis.

Loss of synovial tissue macrophage homeostasis precedes rheumatoid arthritis clinical onset

This study performed an in-depth investigation into the myeloid cellular landscape in the synovium of patients with rheumatoid arthritis (RA), "individuals at risk" of RA, and healthy controls (HC). Flow cytometric analysis demonstrated the presence of a CD40-expressing CD206+CD163+ macrophage population dominating the inflamed RA synovium, associated with disease activity and treatment response. In-depth RNA sequencing and metabolic analysis demonstrated that this macrophage population is transcriptionally distinct, displaying unique inflammatory and tissue-resident gene signatures, has a stable bioenergetic profile, and regulates stromal cell responses. Single-cell RNA sequencing profiling of 67,908 RA and HC synovial tissue cells identified nine transcriptionally distinct macrophage clusters. IL-1B+CCL20+ and SPP1+MT2A+ are the principal macrophage clusters in RA synovium, displaying heightened CD40 gene expression, capable of shaping stromal cell responses, and are importantly enriched before disease onset. Combined, these findings identify the presence of an early pathogenic myeloid signature that shapes the RA joint microenvironment and represents a unique opportunity for early diagnosis and therapeutic intervention.

Profiling Long Noncoding RNA in Psoriatic Skin Using Single-Cell RNA Sequencing

The expressions of long noncoding RNAs (lncRNAs) and their roles in epidermal differentiation have been previously defined using bulk RNA sequencing. Despite their tissue-specific expression profiles, most lncRNAs are not well-annotated at the single-cell level. In this study, we evaluated the use of single-cell RNA sequencing to profile and characterize lncRNAs using data from 6 patients with psoriasis with paired uninvolved and lesional psoriatic skin. Despite their overall lower expression, we were able to detect >7000 skin-expressing lncRNAs and their cellular sources. Differential gene expression analysis revealed 137 differentially expressed lncRNAs in lesional psoriasis skin and identified 169 cell-type-specific lncRNAs. Keratinocytes had the highest number of differentially expressed lncRNA in psoriatic skin, which we validated using spatial transcriptomic data. We further showed that expression of the keratinocyte-specific lncRNA, AC020916.1, upregulated in lesional skin, is significantly correlated with expressions of genes participating in cell proliferation/epidermal differentiation, including SPRR2E and transcription factor ZFP36, particularly in the psoriatic skin. Our study highlights the potential for using single-cell RNA sequencing to profile skin-expressing lncRNA transcripts and to infer their cellular origins, providing a crucial approach that can be applied to the study of other inflammatory skin conditions.