Scrublet: Computational Identification of Cell Doublets in Single-Cell Transcriptomic Data

Single-Nucleus Transcriptomic Atlas of Human Pericoronary Epicardial Adipose Tissue in Normal and Pathological Conditions

BACKGROUND

Pericoronary epicardial adipose tissue (EAT) is a unique visceral fat depot that surrounds the adventitia of the coronary arteries without any anatomic barrier. Clinical studies have demonstrated the association between EAT volume and increased risks for coronary artery disease (CAD). However, the cellular and molecular mechanisms underlying the association remain elusive.

METHODS

We performed single-nucleus RNA sequencing on pericoronary EAT samples collected from 3 groups of subjects: patients undergoing coronary bypass surgery for severe CAD (n=8), patients with CAD with concomitant type 2 diabetes (n=8), and patients with valvular diseases but without concomitant CAD and type 2 diabetes as the control group (n=8). Comparative analyses were performed among groups, including cellular compositional analysis, cell type-resolved transcriptomic changes, gene coexpression network analysis, and intercellular communication analysis. Immunofluorescence staining was performed to confirm the presence of CAD-associated subclusters.

RESULTS

Unsupervised clustering of 73 386 nuclei identified 15 clusters, encompassing all known cell types in the adipose tissue. Distinct subpopulations were identified within primary cell types, including adipocytes, adipose stem and progenitor cells, and macrophages. CD83high macrophages and FOSBhigh adipocytes were significantly expanded in CAD. In comparison to normal controls, both disease groups exhibited dysregulated pathways and altered secretome in the primary cell types. Nevertheless, minimal differences were noted between the disease groups in terms of cellular composition and transcriptome. In addition, our data highlight a potential interplay between dysregulated circadian clock and altered physiological functions in adipocytes of pericoronary EAT. ANXA1 (annexin A1) and SEMA3B (semaphorin 3B) were identified as important adipokines potentially involved in functional changes of pericoronary EAT and CAD pathogenesis.

CONCLUSIONS

We built a complete single-nucleus transcriptomic atlas of human pericoronary EAT in normal and diseased conditions of CAD. Our study lays the foundation for developing novel therapeutic strategies for treating CAD by targeting and modifying pericoronary EAT functions.

Single-cell landscape of functionally cured chronic hepatitis B patients reveals activation of innate and altered CD4-CTL-driven adaptive immunity

BACKGROUND & AIMS

Hepatitis B surface antigen (HBsAg) loss or functional cure (FC) is considered the optimal therapeutic outcome for patients with chronic hepatitis B (CHB). However, the immune-pathological biomarkers and underlying mechanisms of FC remain unclear. In this study we comprehensively interrogate disease-associated cell states identified within intrahepatic tissue and matched PBMCs (peripheral blood mononuclear cells) from patients with CHB or after FC, at the resolution of single cells, to provide novel insights into putative mechanisms underlying FC.

METHODS

We combined single-cell transcriptomics (single-cell RNA sequencing) with multiparametric flow cytometry-based immune phenotyping, and multiplexed immunofluorescence to elucidate the immunopathological cell states associated with CHB vs. FC.

RESULTS

We found that the intrahepatic environment in CHB and FC displays specific cell identities and molecular signatures that are distinct from those found in matched PBMCs. FC is associated with the emergence of an altered adaptive immune response marked by CD4 cytotoxic T lymphocytes, and an activated innate response represented by liver-resident natural killer cells, specific Kupffer cell subtypes and marginated neutrophils. Surprisingly, we found MHC class II-expressing hepatocytes in patients achieving FC, as well as low but persistent levels of covalently closed circular DNA and pregenomic RNA, which may play an important role in FC.

CONCLUSIONS

Our study provides conceptually novel insights into the immuno-pathological control of HBV cure, and opens exciting new avenues for clinical management, biomarker discovery and therapeutic development. We believe that the discoveries from this study, as it relates to the activation of an innate and altered immune response that may facilitate sustained, low-grade inflammation, may have broader implications in the resolution of chronic viral hepatitis.

IMPACT AND IMPLICATIONS

This study dissects the immuno-pathological cell states associated with functionally cured chronic hepatitis B (defined by the loss of HBV surface antigen or HBsAg). We identified the sustained presence of very low viral load, accessory antigen-presenting hepatocytes, adaptive-memory-like natural killer cells, and the emergence of helper CD4 T cells with cytotoxic or effector-like signatures associated with functional cure, suggesting previously unsuspected alterations in the adaptive immune response, as well as a key role for the innate immune response in achieving or maintaining functional cure. Overall, the insights generated from this study may provide new avenues for the development of alternative therapies as well as patient surveillance for better clinical management of chronic hepatitis B.

Single-cell transcriptional profiling of clear cell renal cell carcinoma reveals a tumor-associated endothelial tip cell phenotype

Clear cell renal cell carcinoma (ccRCC) is the most prevalent form of renal cancer, accounting for over 75% of cases. The asymptomatic nature of the disease contributes to late-stage diagnoses and poor survival. Highly vascularized and immune infiltrated microenvironment are prominent features of ccRCC, yet the interplay between vasculature and immune cells, disease progression and response to therapy remains poorly understood. Using droplet-based single-cell RNA sequencing we profile 50,236 transcriptomes from paired tumor and healthy adjacent kidney tissues. Our analysis reveals significant heterogeneity and inter-patient variability of the tumor microenvironment. Notably, we discover a previously uncharacterized vasculature subpopulation associated with epithelial-mesenchymal transition. The cell-cell communication analysis reveals multiple modes of immunosuppressive interactions within the tumor microenvironment, including clinically relevant interactions between tumor vasculature and stromal cells with immune cells. The upregulation of the genes involved in these interactions is associated with worse survival in the TCGA KIRC cohort. Our findings demonstrate the role of tumor vasculature and stromal cell populations in shaping the ccRCC microenvironment and uncover a subpopulation of cells within the tumor vasculature that is associated with an angiogenic phenotype.

The centrosomal protein FGFR1OP controls myosin function in murine intestinal epithelial cells

Recent advances in human genetics have shed light on the genetic factors contributing to inflammatory diseases, particularly Crohn's disease (CD), a prominent form of inflammatory bowel disease. Certain risk genes associated with CD directly influence cytokine biology and cell-specific communication networks. Current CD therapies primarily rely on anti-inflammatory drugs, which are inconsistently effective and lack strategies for promoting epithelial restoration and mucosal balance. To understand CD's underlying mechanisms, we investigated the link between CD and the FGFR1OP gene, which encodes a centrosome protein. FGFR1OP deletion in mouse intestinal epithelial cells disrupted crypt architecture, resulting in crypt loss, inflammation, and fatality. FGFR1OP insufficiency hindered epithelial resilience during colitis. FGFR1OP was crucial for preserving non-muscle myosin II activity, ensuring the integrity of the actomyosin cytoskeleton and crypt cell adhesion. This role of FGFR1OP suggests that its deficiency in genetically predisposed individuals may reduce epithelial renewal capacity, heightening susceptibility to inflammation and disease.

An epigenetically distinct HSC subset supports thymic reconstitution

Hematopoietic stem cells (HSCs) with multilineage potential are critical for effective T cell reconstitution and restoration of the adaptive immune system after allogeneic Hematopoietic Cell Transplantation (allo-HCT). The Kit lo subset of HSCs is enriched for multipotential precursors, 1, 2 but their T-cell lineage potential has not been well-characterized. We therefore studied the thymic reconstituting and T-cell potential of Kit lo HSCs. Using a preclinical allo-HCT model, we demonstrate that Kit lo HSCs support better thymic recovery, and T-cell reconstitution resulting in improved T cell responses to infection post-HCT. Furthermore, Kit lo HSCs with augmented BM lymphopoiesis mitigate age-associated thymic alterations, thus enhancing T-cell recovery in middle-aged hosts. We find the frequency of the Kit lo subset declines with age, providing one explanation for the reduced frequency of T-competent HSCs and reduced T-lymphopoietic potential in BM precursors of aged mice. 3, 4, 5 Chromatin profiling revealed that Kit lo HSCs exhibit higher activity of lymphoid-specifying transcription factors (TFs), including Zbtb1 . Deletion of Zbtb1 in Kit lo HSCs diminished their T-cell potential, while reinstating Zbtb1 in megakaryocytic-biased Kit hi HSCs rescued T-cell potential, in vitro and in vivo . Finally, we discover an analogous Kit lo HSC subset with enhanced lymphoid potential in human bone marrow. Our results demonstrate that Kit lo HSCs with enhanced lymphoid potential have a distinct underlying epigenetic program.

Integration mapping of cardiac fibroblast single-cell transcriptomes elucidates cellular principles of fibrosis in diverse pathologies

Single-cell technology has allowed researchers to probe tissue complexity and dynamics at unprecedented depth in health and disease. However, the generation of high-dimensionality single-cell atlases and virtual three-dimensional tissues requires integrated reference maps that harmonize disparate experimental designs, analytical pipelines, and taxonomies. Here, we present a comprehensive single-cell transcriptome integration map of cardiac fibrosis, which underpins pathophysiology in most cardiovascular diseases. Our findings reveal similarity between cardiac fibroblast (CF) identities and dynamics in ischemic versus pressure overload models of cardiomyopathy. We also describe timelines for commitment of activated CFs to proliferation and myofibrogenesis, profibrotic and antifibrotic polarization of myofibroblasts and matrifibrocytes, and CF conservation across mouse and human healthy and diseased hearts. These insights have the potential to inform knowledge-based therapies.

Synthetic DNA barcodes identify singlets in scRNA-seq datasets and evaluate doublet algorithms

Single-cell RNA sequencing (scRNA-seq) datasets contain true single cells, or singlets, in addition to cells that coalesce during the protocol, or doublets. Identifying singlets with high fidelity in scRNA-seq is necessary to avoid false negative and false positive discoveries. Although several methodologies have been proposed, they are typically tested on highly heterogeneous datasets and lack a priori knowledge of true singlets. Here, we leveraged datasets with synthetically introduced DNA barcodes for a hitherto unexplored application: to extract ground-truth singlets. We demonstrated the feasibility of our framework, "singletCode," to evaluate existing doublet detection methods across a range of contexts. We also leveraged our ground-truth singlets to train a proof-of-concept machine learning classifier, which outperformed other doublet detection algorithms. Our integrative framework can identify ground-truth singlets and enable robust doublet detection in non-barcoded datasets.

Gliovascular transcriptional perturbations in Alzheimer's disease reveal molecular mechanisms of blood brain barrier dysfunction

To uncover molecular changes underlying blood-brain-barrier dysfunction in Alzheimer's disease, we performed single nucleus RNA sequencing in 24 Alzheimer's disease and control brains and focused on vascular and astrocyte clusters as main cell types of blood-brain-barrier gliovascular-unit. The majority of the vascular transcriptional changes were in pericytes. Of the vascular molecular targets predicted to interact with astrocytic ligands, SMAD3, upregulated in Alzheimer's disease pericytes, has the highest number of ligands including VEGFA, downregulated in Alzheimer's disease astrocytes. We validated these findings with external datasets comprising 4,730 pericyte and 150,664 astrocyte nuclei. Blood SMAD3 levels are associated with Alzheimer's disease-related neuroimaging outcomes. We determined inverse relationships between pericytic SMAD3 and astrocytic VEGFA in human iPSC and zebrafish models. Here, we detect vast transcriptome changes in Alzheimer's disease at the gliovascular-unit, prioritize perturbed pericytic SMAD3-astrocytic VEGFA interactions, and validate these in cross-species models to provide a molecular mechanism of blood-brain-barrier disintegrity in Alzheimer's disease.

Single-cell transcriptional profile of CD34+ hematopoietic progenitor cells from del(5q) myelodysplastic syndromes and impact of lenalidomide

While myelodysplastic syndromes with del(5q) (del(5q) MDS) comprises a well-defined hematological subgroup, the molecular basis underlying its origin remains unknown. Using single cell RNA-seq (scRNA-seq) on CD34+ progenitors from del(5q) MDS patients, we have identified cells harboring the deletion, characterizing the transcriptional impact of this genetic insult on disease pathogenesis and treatment response. Interestingly, both del(5q) and non-del(5q) cells present similar transcriptional lesions, indicating that all cells, and not only those harboring the deletion, may contribute to aberrant hematopoietic differentiation. However, gene regulatory network (GRN) analyses reveal a group of regulons showing aberrant activity that could trigger altered hematopoiesis exclusively in del(5q) cells, pointing to a more prominent role of these cells in disease phenotype. In del(5q) MDS patients achieving hematological response upon lenalidomide treatment, the drug reverts several transcriptional alterations in both del(5q) and non-del(5q) cells, but other lesions remain, which may be responsible for potential future relapses. Moreover, lack of hematological response is associated with the inability of lenalidomide to reverse transcriptional alterations. Collectively, this study reveals transcriptional alterations that could contribute to the pathogenesis and treatment response of del(5q) MDS.

The Advancement and Application of the Single-Cell Transcriptome in Biological and Medical Research

Single-cell RNA sequencing technology (scRNA-seq) has been steadily developing since its inception in 2009. Unlike bulk RNA-seq, scRNA-seq identifies the heterogeneity of tissue cells and reveals gene expression changes in individual cells at the microscopic level. Here, we review the development of scRNA-seq, which has gone through iterations of reverse transcription, in vitro transcription, smart-seq, drop-seq, 10 × Genomics, and spatial single-cell transcriptome technologies. The technology of 10 × Genomics has been widely applied in medicine and biology, producing rich research results. Furthermore, this review presents a summary of the analytical process for single-cell transcriptome data and its integration with other omics analyses, including genomes, epigenomes, proteomes, and metabolomics. The single-cell transcriptome has a wide range of applications in biology and medicine. This review analyzes the applications of scRNA-seq in cancer, stem cell research, developmental biology, microbiology, and other fields. In essence, scRNA-seq provides a means of elucidating gene expression patterns in single cells, thereby offering a valuable tool for scientific research. Nevertheless, the current single-cell transcriptome technology is still imperfect, and this review identifies its shortcomings and anticipates future developments. The objective of this review is to facilitate a deeper comprehension of scRNA-seq technology and its applications in biological and medical research, as well as to identify avenues for its future development in alignment with practical needs.

Single-cell and spatially resolved interactomics of tooth-associated keratinocytes in periodontitis

Periodontitis affects billions of people worldwide. To address relationships of periodontal niche cell types and microbes in periodontitis, we generated an integrated single-cell RNA sequencing (scRNAseq) atlas of human periodontium (34-sample, 105918-cell), including sulcular and junctional keratinocytes (SK/JKs). SK/JKs displayed altered differentiation states and were enriched for effector cytokines in periodontitis. Single-cell metagenomics revealed 37 bacterial species with cell-specific tropism. Fluorescence in situ hybridization detected intracellular 16 S and mRNA signals of multiple species and correlated with SK/JK proinflammatory phenotypes in situ. Cell-cell communication analysis predicted keratinocyte-specific innate and adaptive immune interactions. Highly multiplexed immunofluorescence (33-antibody) revealed peri-epithelial immune foci, with innate cells often spatially constrained around JKs. Spatial phenotyping revealed immunosuppressed JK-microniches and SK-localized tertiary lymphoid structures in periodontitis. Here, we demonstrate impacts on and predicted interactomics of SK and JK cells in health and periodontitis, which requires further investigation to support precision periodontal interventions in states of chronic inflammation.

The ENCODE mouse postnatal developmental time course identifies regulatory programs of cell types and cell states

Postnatal genomic regulation significantly influences tissue and organ maturation but is under-studied relative to existing genomic catalogs of adult tissues or prenatal development in mouse. The ENCODE4 consortium generated the first comprehensive single-nucleus resource of postnatal regulatory events across a diverse set of mouse tissues. The collection spans seven postnatal time points, mirroring human development from childhood to adulthood, and encompasses five core tissues. We identified 30 cell types, further subdivided into 69 subtypes and cell states across adrenal gland, left cerebral cortex, hippocampus, heart, and gastrocnemius muscle. Our annotations cover both known and novel cell differentiation dynamics ranging from early hippocampal neurogenesis to a new sex-specific adrenal gland population during puberty. We used an ensemble Latent Dirichlet Allocation strategy with a curated vocabulary of 2,701 regulatory genes to identify regulatory "topics," each of which is a gene vector, linked to cell type differentiation, subtype specialization, and transitions between cell states. We find recurrent regulatory topics in tissue-resident macrophages, neural cell types, endothelial cells across multiple tissues, and cycling cells of the adrenal gland and heart. Cell-type-specific topics are enriched in transcription factors and microRNA host genes, while chromatin regulators dominate mitosis topics. Corresponding chromatin accessibility data reveal dynamic and sex-specific regulatory elements, with enriched motifs matching transcription factors in regulatory topics. Together, these analyses identify both tissue-specific and common regulatory programs in postnatal development across multiple tissues through the lens of the factors regulating transcription.

Early Developmental Origins of Cortical Disorders Modeled in Human Neural Stem Cells

The implications of the early phases of human telencephalic development, involving neural stem cells (NSCs), in the etiology of cortical disorders remain elusive. Here, we explored the expression dynamics of cortical and neuropsychiatric disorder-associated genes in datasets generated from human NSCs across telencephalic fate transitions in vitro and in vivo. We identified risk genes expressed in brain organizers and sequential gene regulatory networks across corticogenesis revealing disease-specific critical phases, when NSCs are more vulnerable to gene dysfunctions, and converging signaling across multiple diseases. Moreover, we simulated the impact of risk transcription factor (TF) depletions on different neural cell types spanning the developing human neocortex and observed a spatiotemporal-dependent effect for each perturbation. Finally, single-cell transcriptomics of newly generated autism-affected patient-derived NSCs in vitro revealed recurrent alterations of TFs orchestrating brain patterning and NSC lineage commitment. This work opens new perspectives to explore human brain dysfunctions at the early phases of development.

One-sentence summary

The temporal analysis of gene regulatory networks in human neural stem cells reveals multiple early critical phases associated with cortical disorders and neuropsychiatric traits.

Glyoxalase 1: Emerging biomarker and therapeutic target in cervical cancer progression

INTRODUCTION

Cervical cancer presents a significant global health challenge, disproportionately impacting underserved populations with limited access to healthcare. Early detection and effective management are vital in addressing this public health concern. This study focuses on Glyoxalase-1 (GLO1), an enzyme crucial for methylglyoxal detoxification, in the context of cervical cancer.

METHODS

We assessed GLO1 expression in cervical cancer patient samples using immunohistochemistry. In vitro experiments using HeLa cells were conducted to evaluate the impact of GLO1 inhibition on cell viability and migration. Single-cell RNA sequencing (scRNA-seq) and gene set variation analysis were utilized to investigate the role of GLO1 in the metabolism of cervical cancer. Additionally, public microarray data were analyzed to determine GLO1 expression across various stages of cervical cancer.

RESULTS

Our analysis included 58 cervical cancer patients, and showed that GLO1 is significantly upregulated in cervical cancer tissues compared to normal cervical tissues, independent of pathological findings and disease stage. In vitro experiments indicated that GLO1 inhibition by S-p-bromobenzylglutathione cyclopentyl diester decreased cell viability and migration in cervical cancer cell lines. Analyses of scRNA-seq data and public gene expression datasets corroborated the overexpression of GLO1 and its involvement in cancer metabolism, particularly glycolysis. An examination of expression data from precancerous lesions revealed a progressive increase in GLO1 expression from normal tissue to invasive cervical cancer.

CONCLUSIONS

This study highlights the critical role of GLO1 in the progression of cervical cancer, presenting it as a potential biomarker and therapeutic target. These findings contribute valuable insights towards personalized treatment approaches and augment the ongoing efforts to combat cervical cancer. Further research is necessary to comprehensively explore GLO1's potential in clinical applications.

Integrated single-cell analysis defines the epigenetic basis of castration-resistant prostate luminal cells

Understanding prostate response to castration and androgen receptor signaling inhibitors (ARSI) is critical to improving long-term prostate cancer (PCa) patient survival. Here, we use a multi-omics approach on 229,794 single cells to create a mouse single-cell reference atlas for interpreting mouse prostate biology and castration response. Our reference atlas refines single-cell annotations and provides a chromatin context, which, when coupled with mouse lineage tracing, demonstrates that castration-resistant luminal cells are distinct from the pre-existent urethra-proximal stem/progenitor cells. Molecular pathway analysis and therapeutic studies further implicate AP1 (JUN/FOS), WNT/β-catenin, FOXQ1, NF-κB, and JAK/STAT pathways as major drivers of castration-resistant luminal populations with relevance to human PCa. Our datasets, which can be explored through an interactive portal (https://visportal.roswellpark.org/data/tang/), can aid in developing combination treatments with ARSI for advanced PCa patients.

Integrated single cell analysis reveals co-evolution of malignant B cells and tumor micro-environment in transformed follicular lymphoma

Histological transformation of follicular lymphoma (FL) to aggressive forms is associated with poor outcome. Phenotypic consequences of this evolution and its impact on the tumor microenvironment (TME) remain unknown. We perform single-cell whole genome sequencing (scWGS) and transcriptome sequencing (scWTS) of 11 paired pre/post-transformation patient samples and scWTS of additional samples from patients without transformation. Our analysis reveals evolutionary dynamics of transformation at single-cell resolution, highlighting a shifting TME landscape, with an emerging immune-cell exhaustion signature, co-evolving with the shifting malignant B phenotype in a regulatory ecosystem. Integration of scWGS and scWTS identifies malignant cell pathways upregulated during clonal tumor evolution. Using multi-color immunofluorescence, we transfer these findings to a TME-based transformation biomarker, subsequently validated in two independent pretreatment cohorts. Taken together, our results provide a comprehensive view of the combined genomic and phenotypic evolution of malignant cells during transformation and shifting crosstalk between malignant cells and the TME.

Single-cell-resolved interspecies comparison shows a shared inflammatory axis and a dominant neutrophil-endothelial program in severe COVID-19

A key issue for research on COVID-19 pathogenesis is the lack of biopsies from patients and of samples at the onset of infection. To overcome these hurdles, hamsters were shown to be useful models for studying this disease. Here, we further leverage the model to molecularly survey the disease progression from time-resolved single-cell RNA sequencing data collected from healthy and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected Syrian and Roborovski hamster lungs. We compare our data to human COVID-19 studies, including bronchoalveolar lavage, nasal swab, and postmortem lung tissue, and identify a shared axis of inflammation dominated by macrophages, neutrophils, and endothelial cells, which we show to be transient in Syrian and terminal in Roborovski hamsters. Our data suggest that, following SARS-CoV-2 infection, commitment to a type 1- or type 3-biased immunity determines moderate versus severe COVID-19 outcomes, respectively.

Demultiplexing of single-cell RNA-sequencing data using interindividual variation in gene expression

Motivation

Pooled designs for single-cell RNA sequencing, where many cells from distinct samples are processed jointly, offer increased throughput and reduced batch variation. This study describes expression-aware demultiplexing (EAD), a computational method that employs differential co-expression patterns between individuals to demultiplex pooled samples without any extra experimental steps.

Results

We use synthetic sample pools and show that the top interindividual differentially co-expressed genes provide a distinct cluster of cells per individual, significantly enriching the regulation of metabolism. Our application of EAD to samples of six isogenic inbred mice demonstrated that controlling genetic and environmental effects can solve interindividual variations related to metabolic pathways. We utilized 30 samples from both sepsis and healthy individuals in six batches to assess the performance of classification approaches. The results indicate that combining genetic and EAD results can enhance the accuracy of assignments (Min. 0.94, Mean 0.98, Max. 1). The results were enhanced by an average of 1.4% when EAD and barcoding techniques were combined (Min. 1.25%, Median 1.33%, Max. 1.74%). Furthermore, we demonstrate that interindividual differential co-expression analysis within the same cell type can be used to identify cells from the same donor in different activation states. By analysing single-nuclei transcriptome profiles from the brain, we demonstrate that our method can be applied to nonimmune cells.

Availability and implementation

EAD workflow is available at https://isarnassiri.github.io/scDIV/ as an R package called scDIV (acronym for single-cell RNA-sequencing data demultiplexing using interindividual variations).

Single-cell analysis reveals a subpopulation of adipose progenitor cells that impairs glucose homeostasis

Adipose progenitor cells (APCs) are heterogeneous stromal cells and help to maintain metabolic homeostasis. However, the influence of obesity on human APC heterogeneity and the role of APC subpopulations on regulating glucose homeostasis remain unknown. Here, we find that APCs in human visceral adipose tissue contain four subsets. The composition and functionality of APCs are altered in patients with type 2 diabetes (T2D). CD9+CD55low APCs are the subset which is significantly increased in T2D patients. Transplantation of these cells from T2D patients into adipose tissue causes glycemic disturbance. Mechanistically, CD9+CD55low APCs promote T2D development through producing bioactive proteins to form a detrimental niche, leading to upregulation of adipocyte lipolysis. Depletion of pathogenic APCs by inducing intracellular diphtheria toxin A expression or using a hunter-killer peptide improves obesity-related glycemic disturbance. Collectively, our data provide deeper insights in human APC functionality and highlights APCs as a potential therapeutic target to combat T2D. All mice utilized in this study are male.

Isthmus progenitor cells contribute to homeostatic cellular turnover and support regeneration following intestinal injury

The currently accepted intestinal epithelial cell organization model proposes that Lgr5+ crypt-base columnar (CBC) cells represent the sole intestinal stem cell (ISC) compartment. However, previous studies have indicated that Lgr5+ cells are dispensable for intestinal regeneration, leading to two major hypotheses: one favoring the presence of a quiescent reserve ISC and the other calling for differentiated cell plasticity. To investigate these possibilities, we studied crypt epithelial cells in an unbiased fashion via high-resolution single-cell profiling. These studies, combined with in vivo lineage tracing, show that Lgr5 is not a specific ISC marker and that stemness potential exists beyond the crypt base and resides in the isthmus region, where undifferentiated cells participate in intestinal homeostasis and regeneration following irradiation (IR) injury. Our results provide an alternative model of intestinal epithelial cell organization, suggesting that stemness potential is not restricted to CBC cells, and neither de-differentiation nor reserve ISC are drivers of intestinal regeneration.

Cell type mapping of inflammatory muscle diseases highlights selective myofiber vulnerability in inclusion body myositis

Inclusion body myositis (IBM) is the most prevalent inflammatory muscle disease in older adults with no effective therapy available. In contrast to other inflammatory myopathies such as subacute, immune-mediated necrotizing myopathy (IMNM), IBM follows a chronic disease course with both inflammatory and degenerative features of pathology. Moreover, causal factors and molecular drivers of IBM progression are largely unknown. Therefore, we paired single-nucleus RNA sequencing with spatial transcriptomics from patient muscle biopsies to map cell-type-specific drivers underlying IBM pathogenesis compared with IMNM muscles and noninflammatory skeletal muscle samples. In IBM muscles, we observed a selective loss of type 2 myonuclei paralleled by increased levels of cytotoxic T and conventional type 1 dendritic cells. IBM myofibers were characterized by either upregulation of cell stress markers featuring GADD45A and NORAD or protein degradation markers including RNF7 associated with p62 aggregates. GADD45A upregulation was preferentially seen in type 2A myofibers associated with severe tissue inflammation. We also noted IBM-specific upregulation of ACHE encoding acetylcholinesterase, which can be regulated by NORAD activity and result in functional denervation of myofibers. Our results provide promising insights into possible mechanisms of myofiber degeneration in IBM and suggest a selective type 2 fiber vulnerability linked to genomic stress and denervation pathways.

MBNL1 Regulates Programmed Postnatal Switching Between Regenerative and Differentiated Cardiac States

BACKGROUND

Discovering determinants of cardiomyocyte maturity is critical for deeply understanding the maintenance of differentiated states and potentially reawakening endogenous regenerative programs in adult mammalian hearts as a therapeutic strategy. Forced dedifferentiation paired with oncogene expression is sufficient to drive cardiac regeneration, but elucidation of endogenous developmental regulators of the switch between regenerative and mature cardiomyocyte cell states is necessary for optimal design of regenerative approaches for heart disease. MBNL1 (muscleblind-like 1) regulates fibroblast, thymocyte, and erythroid differentiation and proliferation. Hence, we examined whether MBNL1 promotes and maintains mature cardiomyocyte states while antagonizing cardiomyocyte proliferation.

METHODS

MBNL1 gain- and loss-of-function mouse models were studied at several developmental time points and in surgical models of heart regeneration. Multi-omics approaches were combined with biochemical, histological, and in vitro assays to determine the mechanisms through which MBNL1 exerts its effects.

RESULTS

MBNL1 is coexpressed with a maturation-association genetic program in the heart and is regulated by the MEIS1/calcineurin signaling axis. Targeted MBNL1 overexpression early in development prematurely transitioned cardiomyocytes to hypertrophic growth, hypoplasia, and dysfunction, whereas loss of MBNL1 function increased cardiomyocyte cell cycle entry and proliferation through altered cell cycle inhibitor transcript stability. Moreover, MBNL1-dependent stabilization of estrogen-related receptor signaling was essential for maintaining cardiomyocyte maturity in adult myocytes. In accordance with these data, modulating MBNL1 dose tuned the temporal window of neonatal cardiac regeneration, where increased MBNL1 expression arrested myocyte proliferation and regeneration and MBNL1 deletion promoted regenerative states with prolonged myocyte proliferation. However, MBNL1 deficiency was insufficient to promote regeneration in the adult heart because of cell cycle checkpoint activation.

CONCLUSIONS

Here, MBNL1 was identified as an essential regulator of cardiomyocyte differentiated states, their developmental switch from hyperplastic to hypertrophic growth, and their regenerative potential through controlling an entire maturation program by stabilizing adult myocyte mRNAs during postnatal development and throughout adulthood. Targeting loss of cardiomyocyte maturity and downregulation of cell cycle inhibitors through MBNL1 deletion was not sufficient to promote adult regeneration.

Cell-type-specific effects of age and sex on human cortical neurons

Altered transcriptional and epigenetic regulation of brain cell types may contribute to cognitive changes with advanced age. Using single-nucleus multi-omic DNA methylation and transcriptome sequencing (snmCT-seq) in frontal cortex from young adult and aged donors, we found widespread age- and sex-related variation in specific neuron types. The proportion of inhibitory SST- and VIP-expressing neurons was reduced in aged donors. Excitatory neurons had more profound age-related changes in their gene expression and DNA methylation than inhibitory cells. Hundreds of genes involved in synaptic activity, including EGR1, were less expressed in aged adults. Genes located in subtelomeric regions increased their expression with age and correlated with reduced telomere length. We further mapped cell-type-specific sex differences in gene expression and X-inactivation escape genes. Multi-omic single-nucleus epigenomes and transcriptomes provide new insight into the effects of age and sex on human neurons.

Memory/Active T-Cell Activation Is Associated with Immunotherapeutic Response in Fumarate Hydratase-Deficient Renal Cell Carcinoma

PURPOSE

Fumarate hydratase-deficient renal cell carcinoma (FH-deficient RCC) is a rare and lethal subtype of kidney cancer. However, the optimal treatments and molecular correlates of benefits for FH-deficient RCC are currently lacking.

EXPERIMENTAL DESIGN

A total of 91 patients with FH-deficient RCC from 15 medical centers between 2009 and 2022 were enrolled in this study. Genomic and bulk RNA-sequencing (RNA-seq) were performed on 88 and 45 untreated FH-deficient RCCs, respectively. Single-cell RNA-seq was performed to identify biomarkers for treatment response. Main outcomes included disease-free survival (DFS) for localized patients, objective response rate (ORR), progression-free survival (PFS), and overall survival (OS) for patients with metastasis.

RESULTS

In the localized setting, we found that a cell-cycle progression signature enabled to predict disease progression. In the metastatic setting, first-line immune checkpoint inhibitor plus tyrosine kinase inhibitor (ICI+TKI) combination therapy showed satisfactory safety and was associated with a higher ORR (43.2% vs. 5.6%), apparently superior PFS (median PFS, 17.3 vs. 9.6 months, P = 0.016) and OS (median OS, not reached vs. 25.7 months, P = 0.005) over TKI monotherapy. Bulk and single-cell RNA-seq data revealed an enrichment of memory and effect T cells in responders to ICI plus TKI combination therapy. Furthermore, we identified a signature of memory and effect T cells that was associated with the effectiveness of ICI plus TKI combination therapy.

CONCLUSIONS

ICI plus TKI combination therapy may represent a promising treatment option for metastatic FH-deficient RCC. A memory/active T-cell-derived signature is associated with the efficacy of ICI+TKI but necessitates further validation.

Glioblastoma-Infiltrating CD8+ T Cells Are Predominantly a Clonally Expanded GZMK+ Effector Population

Recent clinical trials have highlighted the limited efficacy of T cell-based immunotherapy in patients with glioblastoma (GBM). To better understand the characteristics of tumor-infiltrating lymphocytes (TIL) in GBM, we performed cellular indexing of transcriptomes and epitopes by sequencing and single-cell RNA sequencing with paired V(D)J sequencing, respectively, on TILs from two cohorts of patients totaling 15 patients with high-grade glioma, including GBM or astrocytoma, IDH-mutant, grade 4 (G4A). Analysis of the CD8+ TIL landscape reveals an enrichment of clonally expanded GZMK+ effector T cells in the tumor compared with matched blood, which was validated at the protein level. Furthermore, integration with other cancer types highlights the lack of a canonically exhausted CD8+ T-cell population in GBM TIL. These data suggest that GZMK+ effector T cells represent an important T-cell subset within the GBM microenvironment and may harbor potential therapeutic implications.

SIGNIFICANCE

To understand the limited efficacy of immune-checkpoint blockade in GBM, we applied a multiomics approach to understand the TIL landscape. By highlighting the enrichment of GZMK+ effector T cells and the lack of exhausted T cells, we provide a new potential mechanism of resistance to immunotherapy in GBM. This article is featured in Selected Articles from This Issue, p. 897.

Shiba: A unified computational method for robust identification of differential RNA splicing across platforms

Alternative pre-mRNA splicing (AS) is a fundamental regulatory process that generates transcript diversity and cell type variation. We developed Shiba, a robust method integrating transcript assembly, splicing event identification, read counting, and statistical analysis, to efficiently quantify exon splicing levels across various types of RNA-seq datasets. Compared to existing pipelines, Shiba excels in capturing both annotated and unannotated or cryptic differential splicing events with superior accuracy, sensitivity, and reproducibility. Furthermore, Shiba's unique consideration of junction read imbalance and exon-body read coverage reduces false positives, essential for downstream functional analyses. We have further developed scShiba for single-cell/nucleus (sc/sn) RNA-seq data, enabling the exploration of splicing variations in heterogeneous cell populations. Both simulated and real data demonstrate Shiba's robustness across multiple sample sizes, including n=1 datasets and individual cell clusters from scRNA-seq. Application of Shiba on single replicates of RNA-seq identified new AS-NMD targets, and scShiba on snRNA-seq revealed intricate temporal AS regulation in dopaminergic neurons. Both Shiba and scShiba are provided in Docker/Singularity containers and Snakemake pipeline, enhancing accessibility and reproducibility. The comprehensive capabilities of Shiba and scShiba allow systematic and robust quantification of alternative splicing events, laying a solid foundation for mechanistic exploration of functional complexity in RNA splicing.

Molecular pixelation: spatial proteomics of single cells by sequencing

The spatial distribution of cell surface proteins governs vital processes of the immune system such as intercellular communication and mobility. However, fluorescence microscopy has limited scalability in the multiplexing and throughput needed to drive spatial proteomics discoveries at subcellular level. We present Molecular Pixelation (MPX), an optics-free, DNA sequence-based method for spatial proteomics of single cells using antibody-oligonucleotide conjugates (AOCs) and DNA-based, nanometer-sized molecular pixels. The relative locations of AOCs are inferred by sequentially associating them into local neighborhoods using the sequence-unique DNA pixels, forming >1,000 spatially connected zones per cell in 3D. For each single cell, DNA-sequencing reads are computationally arranged into spatial proteomics networks for 76 proteins. By studying immune cell dynamics using spatial statistics on graph representations of the data, we identify known and new patterns of spatial organization of proteins on chemokine-stimulated T cells, highlighting the potential of MPX in defining cell states by the spatial arrangement of proteins.

Multiplex profiling of developmental cis-regulatory elements with quantitative single-cell expression reporters

The inability to scalably and precisely measure the activity of developmental cis-regulatory elements (CREs) in multicellular systems is a bottleneck in genomics. Here we develop a dual RNA cassette that decouples the detection and quantification tasks inherent to multiplex single-cell reporter assays. The resulting measurement of reporter expression is accurate over multiple orders of magnitude, with a precision approaching the limit set by Poisson counting noise. Together with RNA barcode stabilization via circularization, these scalable single-cell quantitative expression reporters provide high-contrast readouts, analogous to classic in situ assays but entirely from sequencing. Screening >200 regions of accessible chromatin in a multicellular in vitro model of early mammalian development, we identify 13 (8 previously uncharacterized) autonomous and cell-type-specific developmental CREs. We further demonstrate that chimeric CRE pairs generate cognate two-cell-type activity profiles and assess gain- and loss-of-function multicellular expression phenotypes from CRE variants with perturbed transcription factor binding sites. Single-cell quantitative expression reporters can be applied in developmental and multicellular systems to quantitatively characterize native, perturbed and synthetic CREs at scale, with high sensitivity and at single-cell resolution.

An alternative cell cycle coordinates multiciliated cell differentiation

The canonical mitotic cell cycle coordinates DNA replication, centriole duplication and cytokinesis to generate two cells from one1. Some cells, such as mammalian trophoblast giant cells, use cell cycle variants like the endocycle to bypass mitosis2. Differentiating multiciliated cells, found in the mammalian airway, brain ventricles and reproductive tract, are post-mitotic but generate hundreds of centrioles, each of which matures into a basal body and nucleates a motile cilium3,4. Several cell cycle regulators have previously been implicated in specific steps of multiciliated cell differentiation5,6. Here we show that differentiating multiciliated cells integrate cell cycle regulators into a new alternative cell cycle, which we refer to as the multiciliation cycle. The multiciliation cycle redeploys many canonical cell cycle regulators, including cyclin-dependent kinases (CDKs) and their cognate cyclins. For example, cyclin D1, CDK4 and CDK6, which are regulators of mitotic G1-to-S progression, are required to initiate multiciliated cell differentiation. The multiciliation cycle amplifies some aspects of the canonical cell cycle, such as centriole synthesis, and blocks others, such as DNA replication. E2F7, a transcriptional regulator of canonical S-to-G2 progression, is expressed at high levels during the multiciliation cycle. In the multiciliation cycle, E2F7 directly dampens the expression of genes encoding DNA replication machinery and terminates the S phase-like gene expression program. Loss of E2F7 causes aberrant acquisition of DNA synthesis in multiciliated cells and dysregulation of multiciliation cycle progression, which disrupts centriole maturation and ciliogenesis. We conclude that multiciliated cells use an alternative cell cycle that orchestrates differentiation instead of controlling proliferation.

Joint single-cell profiling resolves 5mC and 5hmC and reveals their distinct gene regulatory effects

Oxidative modification of 5-methylcytosine (5mC) by ten-eleven translocation (TET) DNA dioxygenases generates 5-hydroxymethylcytosine (5hmC), the most abundant form of oxidized 5mC. Existing single-cell bisulfite sequencing methods cannot resolve 5mC and 5hmC, leaving the cell-type-specific regulatory mechanisms of TET and 5hmC largely unknown. Here, we present joint single-nucleus (hydroxy)methylcytosine sequencing (Joint-snhmC-seq), a scalable and quantitative approach that simultaneously profiles 5hmC and true 5mC in single cells by harnessing differential deaminase activity of APOBEC3A toward 5mC and chemically protected 5hmC. Joint-snhmC-seq profiling of single nuclei from mouse brains reveals an unprecedented level of epigenetic heterogeneity of both 5hmC and true 5mC at single-cell resolution. We show that cell-type-specific profiles of 5hmC or true 5mC improve multimodal single-cell data integration, enable accurate identification of neuronal subtypes and uncover context-specific regulatory effects on cell-type-specific genes by TET enzymes.

Multimodal decoding of human liver regeneration

The liver has a unique ability to regenerate1,2; however, in the setting of acute liver failure (ALF), this regenerative capacity is often overwhelmed, leaving emergency liver transplantation as the only curative option3-5. Here, to advance understanding of human liver regeneration, we use paired single-nucleus RNA sequencing combined with spatial profiling of healthy and ALF explant human livers to generate a single-cell, pan-lineage atlas of human liver regeneration. We uncover a novel ANXA2+ migratory hepatocyte subpopulation, which emerges during human liver regeneration, and a corollary subpopulation in a mouse model of acetaminophen (APAP)-induced liver regeneration. Interrogation of necrotic wound closure and hepatocyte proliferation across multiple timepoints following APAP-induced liver injury in mice demonstrates that wound closure precedes hepatocyte proliferation. Four-dimensional intravital imaging of APAP-induced mouse liver injury identifies motile hepatocytes at the edge of the necrotic area, enabling collective migration of the hepatocyte sheet to effect wound closure. Depletion of hepatocyte ANXA2 reduces hepatocyte growth factor-induced human and mouse hepatocyte migration in vitro, and abrogates necrotic wound closure following APAP-induced mouse liver injury. Together, our work dissects unanticipated aspects of liver regeneration, demonstrating an uncoupling of wound closure and hepatocyte proliferation and uncovering a novel migratory hepatocyte subpopulation that mediates wound closure following liver injury. Therapies designed to promote rapid reconstitution of normal hepatic microarchitecture and reparation of the gut-liver barrier may advance new areas of therapeutic discovery in regenerative medicine.

Differential Impact In Vivo of Pf4-ΔCre-Mediated and Gp1ba-ΔCre-Mediated Depletion of Cyclooxygenase-1 in Platelets in Mice

BACKGROUND

Low-dose aspirin is widely used for the secondary prevention of cardiovascular disease. The beneficial effects of low-dose aspirin are attributable to its inhibition of platelet Cox (cyclooxygenase)-1-derived thromboxane A2. Until recently, the use of the Pf4 (platelet factor 4) Cre has been the only genetic approach to generating megakaryocyte/platelet ablation of Cox-1 in mice. However, Pf4-ΔCre displays ectopic expression outside the megakaryocyte/platelet lineage, especially during inflammation. The use of the Gp1ba (glycoprotein 1bα) Cre promises a more specific, targeted approach.

METHODS

To evaluate the role of Cox-1 in platelets, we crossed Pf4-ΔCre or Gp1ba-ΔCre mice with Cox-1flox/flox mice to generate platelet Cox-1-/- mice on normolipidemic and hyperlipidemic (Ldlr-/-; low-density lipoprotein receptor) backgrounds.

RESULTS

Ex vivo platelet aggregation induced by arachidonic acid or adenosine diphosphate in platelet-rich plasma was inhibited to a similar extent in Pf4-ΔCre Cox-1-/-/Ldlr-/- and Gp1ba-ΔCre Cox-1-/-/Ldlr-/- mice. In a mouse model of tail injury, Pf4-ΔCre-mediated and Gp1ba-ΔCre-mediated deletions of Cox-1 were similarly efficient in suppressing platelet prostanoid biosynthesis. Experimental thrombogenesis and attendant blood loss were similar in both models. However, the impact on atherogenesis was divergent, being accelerated in the Pf4-ΔCre mice while restrained in the Gp1ba-ΔCres. In the former, accelerated atherogenesis was associated with greater suppression of PGI2 biosynthesis, a reduction in the lipopolysaccharide-evoked capacity to produce PGE2 (prostaglandin E) and PGD2 (prostanglandin D), activation of the inflammasome, elevated plasma levels of IL-1β (interleukin), reduced plasma levels of HDL-C (high-density lipoprotein receptor-cholesterol), and a reduction in the capacity for reverse cholesterol transport. By contrast, in the latter, plasma HDL-C and α-tocopherol were elevated, and MIP-1α (macrophage inflammatory protein-1α) and MCP-1 (monocyte chemoattractant protein 1) were reduced.

CONCLUSIONS

Both approaches to Cox-1 deletion similarly restrain thrombogenesis, but a differential impact on Cox-1-dependent prostanoid formation by the vasculature may contribute to an inflammatory phenotype and accelerated atherogenesis in Pf4-ΔCre mice.

Systematic benchmarking of single-cell ATAC-sequencing protocols

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Unveiling the hidden role of disulfidptosis in kidney renal clear cell carcinoma: a prognostic signature for personalized treatment

The role of disulfidptosis in kidney renal clear cell carcinoma (KIRC) remains unknown. This study investigated disulfidptosis-related biomarkers for KIRC prognosis prediction and individualized treatment. KIRC patients were clustered by disulfidptosis profiles. Differential expression analysis, survival models, and machine learning were used to construct the disulfidptosis-related prognostic signature (DRPS). Characterizations of the tumor immune microenvironment, genetic drivers, drug sensitivity, and immunotherapy response were explored according to the DRPS risk stratification. Markers included in the signature were validated using single-cell, spatial transcriptomics, quantitative RT-qPCR, and immunohistochemistry. In the discovery cohort, we unveiled two clusters of KIRC patients that differed significantly in disulfidptosis regulator expressions and overall survival (OS). After multiple feature selection steps, a DRPS prognostic model with four features (CHAC1, COL7A1, FOXM1, SHOX2) was constructed and validated. Combined with clinical factors, the model demonstrated robust performance in the discovery and external validation cohorts (5-year AUC = 0.793 and 0.846, respectively). KIRC patients with high-risk scores are characterized by inferior OS, less tumor purity, and increased infiltrations of fibroblasts, M1 macrophages, and B cells. High-risk patients also have higher frequencies of BAP1 and AHNAK2 mutation. Besides, the correlation between the DRPS score and the chemotherapy-response signature indicated the potential effect of Gefitinib for high-risk patients. Among the signature genes, FOXM1 is highly expressed in cycling tumor cells and exhibits spatial aggregation, while others are expressed sparsely within tumor samples. The DRPS model enables improved clinical management and personalized KIRC therapy. The identified biomarkers and immune characteristics offer new mechanistic insight into disulfidptosis in KIRC.

The chromatin landscape of pathogenic transcriptional cell states in rheumatoid arthritis

Synovial tissue inflammation is a hallmark of rheumatoid arthritis (RA). Recent work has identified prominent pathogenic cell states in inflamed RA synovial tissue, such as T peripheral helper cells; however, the epigenetic regulation of these states has yet to be defined. Here, we examine genome-wide open chromatin at single-cell resolution in 30 synovial tissue samples, including 12 samples with transcriptional data in multimodal experiments. We identify 24 chromatin classes and predict their associated transcription factors, including a CD8 + GZMK+ class associated with EOMES and a lining fibroblast class associated with AP-1. By integrating with an RA tissue transcriptional atlas, we propose that these chromatin classes represent 'superstates' corresponding to multiple transcriptional cell states. Finally, we demonstrate the utility of this RA tissue chromatin atlas through the associations between disease phenotypes and chromatin class abundance, as well as the nomination of classes mediating the effects of putatively causal RA genetic variants.

Aging represses lung tumorigenesis and alters tumor suppression

Most cancers are diagnosed in persons over the age of sixty, but little is known about how age impacts tumorigenesis. While aging is accompanied by mutation accumulation - widely understood to contribute to cancer risk - it is also associated with numerous other cellular and molecular changes likely to impact tumorigenesis. Moreover, cancer incidence decreases in the oldest part of the population, suggesting that very old age may reduce carcinogenesis. Here we show that aging represses tumor initiation and growth in genetically engineered mouse models of human lung cancer. Moreover, aging dampens the impact of inactivating many, but not all, tumor suppressor genes with the impact of inactivating PTEN, a negative regulator of the PI3K/AKT pathway, weakened to a disproportionate extent. Single-cell transcriptomic analysis revealed that neoplastic cells from tumors in old mice retain many age-related transcriptomic changes, showing that age has an enduring impact that persists through oncogenic transformation. Furthermore, the consequences of PTEN inactivation were strikingly age-dependent, with PTEN deficiency reducing signatures of aging in cancer cells and the tumor microenvironment. Our findings suggest that the relationship between age and lung cancer incidence may reflect an integration of the competing effects of driver mutation accumulation and tumor suppressive effects of aging.

Microglia-astrocyte crosstalk in the amyloid plaque niche of an Alzheimer's disease mouse model, as revealed by spatial transcriptomics

The amyloid plaque niche is a pivotal hallmark of Alzheimer's disease (AD). Here, we employ two high-resolution spatial transcriptomics (ST) platforms, CosMx and Spatial Enhanced Resolution Omics-sequencing (Stereo-seq), to characterize the transcriptomic alterations, cellular compositions, and signaling perturbations in the amyloid plaque niche in an AD mouse model. We discover heterogeneity in the cellular composition of plaque niches, marked by an increase in microglial accumulation. We profile the transcriptomic alterations of glial cells in the vicinity of plaques and conclude that the microglial response to plaques is consistent across different brain regions, while the astrocytic response is more heterogeneous. Meanwhile, as the microglial density of plaque niches increases, astrocytes acquire a more neurotoxic phenotype and play a key role in inducing GABAergic signaling and decreasing glutamatergic signaling in hippocampal neurons. We thus show that the accumulation of microglia around hippocampal plaques disrupts astrocytic signaling, in turn inducing an imbalance in neuronal synaptic signaling.

Advancing skeletal health and disease research with single-cell RNA sequencing

Orthopedic conditions have emerged as global health concerns, impacting approximately 1.7 billion individuals worldwide. However, the limited understanding of the underlying pathological processes at the cellular and molecular level has hindered the development of comprehensive treatment options for these disorders. The advent of single-cell RNA sequencing (scRNA-seq) technology has revolutionized biomedical research by enabling detailed examination of cellular and molecular diversity. Nevertheless, investigating mechanisms at the single-cell level in highly mineralized skeletal tissue poses technical challenges. In this comprehensive review, we present a streamlined approach to obtaining high-quality single cells from skeletal tissue and provide an overview of existing scRNA-seq technologies employed in skeletal studies along with practical bioinformatic analysis pipelines. By utilizing these methodologies, crucial insights into the developmental dynamics, maintenance of homeostasis, and pathological processes involved in spine, joint, bone, muscle, and tendon disorders have been uncovered. Specifically focusing on the joint diseases of degenerative disc disease, osteoarthritis, and rheumatoid arthritis using scRNA-seq has provided novel insights and a more nuanced comprehension. These findings have paved the way for discovering novel therapeutic targets that offer potential benefits to patients suffering from diverse skeletal disorders.

Rapid response of fly populations to gene dosage across development and generations

Although the effects of genetic and environmental perturbations on multicellular organisms are rarely restricted to single phenotypic layers, our current understanding of how developmental programs react to these challenges remains limited. Here, we have examined the phenotypic consequences of disturbing the bicoid regulatory network in early Drosophila embryos. We generated flies with two extra copies of bicoid, which causes a posterior shift of the network's regulatory outputs and a decrease in fitness. We subjected these flies to EMS mutagenesis, followed by experimental evolution. After only 8-15 generations, experimental populations have normalized patterns of gene expression and increased survival. Using a phenomics approach, we find that populations were normalized through rapid increases in embryo size driven by maternal changes in metabolism and ovariole development. We extend our results to additional populations of flies, demonstrating predictability. Together, our results necessitate a broader view of regulatory network evolution at the systems level.

NG2 glia protect against prion neurotoxicity by inhibiting microglia-to-neuron prostaglandin E2 signaling

Oligodendrocyte-lineage cells, including NG2 glia, undergo prominent changes in various neurodegenerative disorders. Here, we identify a neuroprotective role for NG2 glia against prion toxicity. NG2 glia were activated after prion infection in cerebellar organotypic cultured slices (COCS) and in brains of prion-inoculated mice. In both model systems, depletion of NG2 glia exacerbated prion-induced neurodegeneration and accelerated prion pathology. Loss of NG2 glia enhanced the biosynthesis of prostaglandin E2 (PGE2) by microglia, which augmented prion neurotoxicity through binding to the EP4 receptor. Pharmacological or genetic inhibition of PGE2 biosynthesis attenuated prion-induced neurodegeneration in COCS and mice, reduced the enhanced neurodegeneration in NG2-glia-depleted COCS after prion infection, and dampened the acceleration of prion disease in NG2-glia-depleted mice. These data unveil a non-cell-autonomous interaction between NG2 glia and microglia in prion disease and suggest that PGE2 signaling may represent an actionable target against prion diseases.

CCR5-overexpressing mesenchymal stem cells protect against experimental autoimmune uveitis: insights from single-cell transcriptome analysis

Autoimmune uveitis is a leading cause of severe vision loss, and animal models provide unique opportunities for studying its pathogenesis and therapeutic strategies. Here we employ scRNA-seq, RNA-seq and various molecular and cellular approaches to characterize mouse models of classical experimental autoimmune uveitis (EAU), revealing that EAU causes broad retinal neuron degeneration and marker downregulation, and that Müller glia may act as antigen-presenting cells. Moreover, EAU immune response is primarily driven by Th1 cells, and results in dramatic upregulation of CC chemokines, especially CCL5, in the EAU retina. Accordingly, overexpression of CCR5, a CCL5 receptor, in mesenchymal stem cells (MSCs) enhances their homing capacity and improves their immunomodulatory outcomes in preventing EAU, by reducing infiltrating T cells and activated microglia and suppressing Nlrp3 inflammasome activation. Taken together, our data not only provide valuable insights into the molecular characteristics of EAU but also open an avenue for innovative MSC-based therapy.

C5aR1 inhibition reprograms tumor associated macrophages and reverses PARP inhibitor resistance in breast cancer

Although Poly (ADP-ribose) polymerase (PARP) inhibitors (PARPi) have been approved in multiple diseases, including BRCA1/2 mutant breast cancer, responses are usually transient requiring the deployment of combination therapies for optimal efficacy. Here we thus explore mechanisms underlying sensitivity and resistance to PARPi using two intrinsically PARPi sensitive (T22) and resistant (T127) syngeneic murine breast cancer models in female mice. We demonstrate that tumor associated macrophages (TAM) potentially contribute to the differential sensitivity to PARPi. By single-cell RNA-sequencing, we identify a TAM_C3 cluster, expressing genes implicated in anti-inflammatory activity, that is enriched in PARPi resistant T127 tumors and markedly decreased by PARPi in T22 tumors. Rps19/C5aR1 signaling is selectively elevated in TAM_C3. C5aR1 inhibition or transferring C5aR1hi cells increases and decreases PARPi sensitivity, respectively. High C5aR1 levels in human breast cancers are associated with poor responses to immune checkpoint blockade. Thus, targeting C5aR1 may selectively deplete pro-tumoral macrophages and engender sensitivity to PARPi and potentially other therapies.

Single-cell multi-cohort dissection of the schizophrenia transcriptome

The complexity and heterogeneity of schizophrenia have hindered mechanistic elucidation and the development of more effective therapies. Here, we performed single-cell dissection of schizophrenia-associated transcriptomic changes in the human prefrontal cortex across 140 individuals in two independent cohorts. Excitatory neurons were the most affected cell group, with transcriptional changes converging on neurodevelopment and synapse-related molecular pathways. Transcriptional alterations included known genetic risk factors, suggesting convergence of rare and common genomic variants on neuronal population-specific alterations in schizophrenia. Based on the magnitude of schizophrenia-associated transcriptional change, we identified two populations of individuals with schizophrenia marked by expression of specific excitatory and inhibitory neuronal cell states. This single-cell atlas links transcriptomic changes to etiological genetic risk factors, contextualizing established knowledge within the human cortical cytoarchitecture and facilitating mechanistic understanding of schizophrenia pathophysiology and heterogeneity.

Podocyte Pathogenic Bone Morphogenetic Protein-2 Pathway and Immune Cell Behaviors in Primary Membranous Nephropathy

Primary membranous nephropathy (PMN) is one of the leading causes of end-stage renal disease, and the most frequent cause of massive proteinuria in nondiabetic adults, resulting in fatal complications. However, the underlying pathomechanisms of PMN remain largely unclear. Here, single-cell RNA sequencing is employed to analyze kidney biopsies from eleven PMN patients and seven healthy subjects. Profiling 44 060 cells from patients allowed us to characterize the cellular composition and cell-type-specific gene expression in the PMN kidney. The complement-induced BMP2/pSMAD1/COL4 pathway is identified as the pathogenic pathway in podocytes, bridging two key events, i.e., complement system activation and glomerular basement membrane thickening in PMN. Augmented infiltration and activation of myeloid leukocytes and B lymphocytes are found, profiling delicate crosstalk of immune cells in PMN kidneys. Overall, these results provide valuable insights into the roles of podocytes and immune cells in PMN, and comprehensive resources toward the complete understanding of PMN pathophysiology.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Single-cell and spatial transcriptomics analysis of non-small cell lung cancer

Lung cancer is the second most frequently diagnosed cancer and the leading cause of cancer-related mortality worldwide. Tumour ecosystems feature diverse immune cell types. Myeloid cells, in particular, are prevalent and have a well-established role in promoting the disease. In our study, we profile approximately 900,000 cells from 25 treatment-naive patients with adenocarcinoma and squamous-cell carcinoma by single-cell and spatial transcriptomics. We note an inverse relationship between anti-inflammatory macrophages and NK cells/T cells, and with reduced NK cell cytotoxicity within the tumour. While we observe a similar cell type composition in both adenocarcinoma and squamous-cell carcinoma, we detect significant differences in the co-expression of various immune checkpoint inhibitors. Moreover, we reveal evidence of a transcriptional "reprogramming" of macrophages in tumours, shifting them towards cholesterol export and adopting a foetal-like transcriptional signature which promotes iron efflux. Our multi-omic resource offers a high-resolution molecular map of tumour-associated macrophages, enhancing our understanding of their role within the tumour microenvironment.

Systematic evaluation of single-cell RNA-seq analyses performance based on long-read sequencing platforms

INTRODUCTION

The rapid development of next-generation sequencing (NGS)-based single-cell RNA sequencing (scRNA-seq) allows for detecting and quantifying gene expression in a high-throughput manner, providing a powerful tool for comprehensively understanding cellular function in various biological processes. However, the NGS-based scRNA-seq only quantifies gene expression and cannot reveal the exact transcript structures (isoforms) of each gene due to the limited read length. On the other hand, the long read length of third-generation sequencing (TGS) technologies, including Oxford Nanopore Technologies (ONT) and Pacific Biosciences (PacBio), enable direct reading of intact cDNA molecules.

OBJECTIVES

Both ONT and PacBio have been used in conjunction with scRNA-seq, but their performance in single-cell analyses has not been systematically evaluated.

METHODS

To address this, we generated ONT and PacBio data from the same single-cell cDNA libraries containing different amount of cells.

RESULTS

Using NGS as a control, we assessed the performance of each platform in cell type identification. Additionally, the reliability in identifying novel isoforms and allele-specific gene/isoform expression by both platforms was verified, providing a systematic evaluation to design the sequencing strategies in single-cell transcriptome studies.

CONCLUSION

Beyond gene expression analysis, which the NGS-based scRNA-seq only affords, TGS-based scRNA-seq achieved gene splicing analyses, identifying novel isoforms. Attribute to higher sequencing quality of PacBio, it outperforms ONT in accuracy of novel transcripts identification and allele-specific gene/isoform expression.

The Macrophage Landscape Across the Lifespan of a Human Cardiac Allograft

BACKGROUND

Much of our knowledge of organ rejection after transplantation is derived from rodent models.

METHODS

We used single-nucleus RNA sequencing to investigate the inflammatory myocardial microenvironment in human pediatric cardiac allografts at different stages after transplantation. We distinguished donor- from recipient-derived cells using naturally occurring genetic variants embedded in single-nucleus RNA sequencing data.

RESULTS

Donor-derived tissue resident macrophages, which accompany the allograft into the recipient, are lost over time after transplantation. In contrast, monocyte-derived macrophages from the recipient populate the heart within days after transplantation and form 2 macrophage populations: recipient MP1 and recipient MP2. Recipient MP2s have cell signatures similar to donor-derived resident macrophages; however, they lack signatures of pro-reparative phagocytic activity typical of donor-derived resident macrophages and instead express profibrotic genes. In contrast, recipient MP1s express genes consistent with hallmarks of cellular rejection. Our data suggest that recipient MP1s activate a subset of natural killer cells, turning them into a cytotoxic cell population through feed-forward signaling between recipient MP1s and natural killer cells.

CONCLUSIONS

Our findings reveal an imbalance of donor-derived and recipient-derived macrophages in the pediatric cardiac allograft that contributes to allograft failure.

Transcriptomic and spatial dissection of human ex vivo right atrial tissue reveals proinflammatory microvascular changes in ischemic heart disease

Cardiovascular disease plays a central role in the electrical and structural remodeling of the right atrium, predisposing to arrhythmias, heart failure, and sudden death. Here, we dissect with single-nuclei RNA sequencing (snRNA-seq) and spatial transcriptomics the gene expression changes in the human ex vivo right atrial tissue and pericardial fluid in ischemic heart disease, myocardial infarction, and ischemic and non-ischemic heart failure using asymptomatic patients with valvular disease who undergo preventive surgery as the control group. We reveal substantial differences in disease-associated gene expression in all cell types, collectively suggesting inflammatory microvascular dysfunction and changes in the right atrial tissue composition as the valvular and vascular diseases progress into heart failure. The data collectively suggest that investigation of human cardiovascular disease should expand to all functionally important parts of the heart, which may help us to identify mechanisms promoting more severe types of the disease.

Integrated analysis of single-cell and bulk RNA sequencing data reveals a myeloid cell-related regulon predicting neoadjuvant immunotherapy response across cancers

BACKGROUND

Immunotherapy has brought about a paradigm shift in the treatment of cancer. However, the majority of patients exhibit resistance or become refractory to immunotherapy, and the underlying mechanisms remain to be explored.

METHODS

Sing-cell RNA sequencing (scRNA‑seq) datasets derived from 1 pretreatment and 1 posttreatment achieving pathological complete response (pCR) patient with lung adenocarcinoma (LUAD) who received neoadjuvant immunotherapy were collected, and pySCENIC was used to find the gene regulatory network (GRN) between cell types and immune checkpoint inhibitor (ICI) response. A regulon predicting ICI response was identified and validated using large‑scale pan-cancer data, including a colorectal cancer scRNA‑seq dataset, a breast cancer scRNA‑seq dataset, The Cancer Genome Atlas (TCGA) pan-cancer cohort, and 5 ICI transcriptomic cohorts. Symphony reference mapping was performed to construct the myeloid cell map.

RESULTS

Thirteen major cluster cell types were identified by comparing pretreatment and posttreatment patients, and the fraction of myeloid cells was higher in the posttreatment group (19.0% vs. 11.8%). A PPARG regulon (containing 23 target genes) was associated with ICI response, and its function was validated by a colorectal cancer scRNA‑seq dataset, a breast cancer scRNA‑seq dataset, TCGA pan-cancer cohort, and 5 ICI transcriptomic cohorts. Additionally, a myeloid cell map was developed, and cluster I, II, and III myeloid cells with high expression of PPARG were identified. Moreover, we constructed a website called PPARG ( https://pparg.online/PPARG/ or http://43.134.20.130:3838/PPARG/ ), which provides a powerful discovery tool and resource value for researchers.

CONCLUSIONS

The PPARG regulon is a predictor of ICI response. The myeloid cell map enables the identification of PPARG subclusters in public scRNA-seq datasets and provides a powerful discovery tool and resource value.