Doublet identification in single-cell sequencing data using scDblFinder

Transient chromatin decompaction at the start of D. melanogaster male embryonic germline development

Embryonic germ cells develop rapidly to establish the foundation for future developmental trajectories, and in this process, they make critical lineage choices including the configuration of their unique identity and a decision on sex. Here, we use single-cell genomics patterns for the entire embryonic germline in Drosophila melanogaster along with the somatic gonadal precursors after embryonic gonad coalescence to investigate molecular mechanisms involved in the setting up and regulation of the germline program. Profiling of the early germline chromatin landscape revealed sex- and stage-specific features. In the male germline immediately after zygotic activation, the chromatin structure underwent a brief remodeling phase during which nucleosome density was lower and deconcentrated from promoter regions. These findings echoed enrichment analysis results of our genomics data in which top candidates were factors with the ability to mediate large-scale chromatin reorganization. Together, they point to the importance of chromatin regulation in the early germline and raise the possibility of a conserved epigenetic reprogramming-like process required for proper initiation of germline development.

Prdm1 positively regulates liver Group 1 ILCs cancer immune surveillance and preserves functional heterogeneity

Group 1 innate lymphoid cells (ILCs) comprise conventional natural killer (cNK) cells and type 1 innate lymphoid cells (ILC1s). The main functions of liver cNK cells and ILC1s not only include directly killing target cells but also regulating local immune microenvironment of the liver through the secretion of cytokines. Uncovering the intricate mechanisms by which transcriptional factors regulate and influence the functions of liver cNK cells and ILC1s, particularly within the context of liver tumors, presents a significant opportunity to amplify the effectiveness of immunotherapies against liver malignancies. Using Ncr1-drived conditional knockout mouse model, our study reveals the regulatory role of Prdm1 in shaping the composition and maturation of cNK cells. Although Prdm1 did not affect the killing function of cNK cells in an in vivo cytotoxicity model, a significant increase in cancer metastasis was observed in Prdm1 knockout mice. Interferon-gamma (IFN-γ), granzyme B, and perforin secretion decreased significantly in Prdm1-deficient cNK cells and liver ILC1s. Single-cell RNA sequencing (scRNA-seq) data also provided evidences that Prdm1 maintains functional subsets of cNK cells and liver ILC1s and facilitates communications between cNK cells, liver ILC1s, and macrophages. The present study unveiled a novel regulatory mechanism of Prdm1 in cNK cells and liver ILC1s, showing promising potential for developing innovative immune therapy strategies against liver cancer.

Multisensory gamma stimulation mitigates the effects of demyelination induced by cuprizone in male mice

Demyelination is a common pathological feature in a wide range of diseases, characterized by the loss of myelin sheath and myelin-supporting oligodendrocytes. These losses lead to impaired axonal function, increased vulnerability of axons to damage, and result in significant brain atrophy and neuro-axonal degeneration. Multiple pathomolecular processes contribute to neuroinflammation, oligodendrocyte cell death, and progressive neuronal dysfunction. In this study, we use the cuprizone mouse model of demyelination to investigate long-term non-invasive gamma entrainment using sensory stimulation as a potential therapeutic intervention for promoting myelination and reducing neuroinflammation in male mice. Here, we show that multisensory gamma stimulation mitigates demyelination, promotes oligodendrogenesis, preserves functional integrity and synaptic plasticity, attenuates oligodendrocyte ferroptosis-induced cell death, and reduces brain inflammation. Thus, the protective effects of multisensory gamma stimulation on myelin and anti-neuroinflammatory properties support its potential as a therapeutic approach for demyelinating disorders.

A single-nucleus transcriptomic atlas of medium spiny neurons in the rat nucleus accumbens

Neural processing of rewarding stimuli involves several distinct regions, including the nucleus accumbens (NAc). The majority of NAc neurons are GABAergic projection neurons known as medium spiny neurons (MSNs). MSNs are broadly defined by dopamine receptor expression, but evidence suggests that a wider array of subtypes exist. To study MSN heterogeneity, we analyzed single-nucleus RNA sequencing data from the largest available rat NAc dataset. Analysis of 48,040 NAc MSN nuclei identified major populations belonging to the striosome and matrix compartments. Integration with mouse and human data indicated consistency across species and disease-relevance scoring using genome-wide association study results revealed potentially differential roles for MSN populations in substance use disorders. Additional high-resolution clustering identified 34 transcriptomically distinct subtypes of MSNs definable by a limited number of marker genes. Together, these data demonstrate the diversity of MSNs in the NAc and provide a basis for more targeted genetic manipulation of specific populations.

Molecular and cellular dynamics of the developing human neocortex at single-cell resolution

The development of the human neocortex is a highly dynamic process and involves complex cellular trajectories controlled by cell-type-specific gene regulation 1 . Here, we collected paired single-nucleus chromatin accessibility and transcriptome data from 38 human neocortical samples encompassing both the prefrontal cortex and primary visual cortex. These samples span five main developmental stages, ranging from the first trimester to adolescence. In parallel, we performed spatial transcriptomic analysis on a subset of the samples to illustrate spatial organization and intercellular communication. This atlas enables us to catalog cell type-, age-, and area-specific gene regulatory networks underlying neural differentiation. Moreover, combining single-cell profiling, progenitor purification, and lineage-tracing experiments, we have untangled the complex lineage relationships among progenitor subtypes during the transition from neurogenesis to gliogenesis in the human neocortex. We identified a tripotential intermediate progenitor subtype, termed Tri-IPC, responsible for the local production of GABAergic neurons, oligodendrocyte precursor cells, and astrocytes. Remarkably, most glioblastoma cells resemble Tri-IPCs at the transcriptomic level, suggesting that cancer cells hijack developmental processes to enhance growth and heterogeneity. Furthermore, by integrating our atlas data with large-scale GWAS data, we created a disease-risk map highlighting enriched ASD risk in second-trimester intratelencephalic projection neurons. Our study sheds light on the gene regulatory landscape and cellular dynamics of the developing human neocortex.

An iPSC-derived small intestine-on-chip with self-organizing epithelial, mesenchymal, and neural cells

Human induced pluripotent stem cell (hiPSC)-derived intestinal organoids are valuable tools for researching developmental biology and personalized therapies, but their closed topology and relative immature state limit applications. Here, we use organ-on-chip technology to develop a hiPSC-derived intestinal barrier with apical and basolateral access in a more physiological in vitro microenvironment. To replicate growth factor gradients along the crypt-villus axis, we locally expose the cells to expansion and differentiation media. In these conditions, intestinal epithelial cells self-organize into villus-like folds with physiological barrier integrity, and myofibroblasts and neurons emerge and form a subepithelial tissue in the bottom channel. The growth factor gradients efficiently balance dividing and mature cell types and induce an intestinal epithelial composition, including absorptive and secretory lineages, resembling the composition of the human small intestine. This well-characterized hiPSC-derived intestine-on-chip system can facilitate personalized studies on physiological processes and therapy development in the human small intestine.

Transcriptomic analysis of the human habenula in schizophrenia

Pathophysiology of many neuropsychiatric disorders, including schizophrenia (SCZD), is linked to habenula (Hb) function. While pharmacotherapies and deep brain stimulation targeting the Hb are emerging as promising therapeutic treatments, little is known about the cell type-specific transcriptomic organization of the human Hb or how it is altered in SCZD. Here we define the molecular neuroanatomy of the human Hb and identify transcriptomic changes in individuals with SCZD compared to neurotypical controls. Utilizing Hb-enriched postmortem human brain tissue, we performed single nucleus RNA-sequencing (snRNA-seq; n=7 neurotypical donors) and identified 17 molecularly defined Hb cell types across 16,437 nuclei, including 3 medial and 7 lateral Hb populations, several of which were conserved between rodents and humans. Single molecule fluorescent in situ hybridization (smFISH; n=3 neurotypical donors) validated snRNA-seq Hb cell types and mapped their spatial locations. Bulk RNA-sequencing and cell type deconvolution in Hb-enriched tissue from 35 individuals with SCZD and 33 neurotypical controls yielded 45 SCZD-associated differentially expressed genes (DEGs, FDR < 0.05), with 32 (71%) unique to Hb-enriched tissue. eQTL analysis identified 717 independent SNP-gene pairs (FDR < 0.05), where either the SNP is a SCZD risk variant (16 pairs) or the gene is a SCZD DEG (7 pairs). eQTL and SCZD risk colocalization analysis identified 16 colocalized genes. These results identify topographically organized cell types with distinct molecular signatures in the human Hb and demonstrate unique genetic changes associated with SCZD, thereby providing novel molecular insights into the role of Hb in neuropsychiatric disorders.

One Sentence Summary

Transcriptomic analysis of the human habenula and identification of molecular changes associated with schizophrenia risk and illness state.

Maternal obesity and offspring neurodevelopment are associated with hypoxic gene expression in term human placenta

One third of women in the United States are affected by obesity during pregnancy. Maternal obesity (MO) is associated with an increased risk of neurodevelopmental and metabolic disorders in the offspring. The placenta, located at the maternal-fetal interface, is a key organ determining fetal development and likely contributes to programming of long-term offspring health. We profiled the term placental transcriptome in humans (pre-pregnancy BMI 35+ [MO condition] or 18.5-25 [lean condition]) using single-nucleus RNA-seq to compare expression profiles in MO versus lean conditions, and to reveal potential mechanisms underlying offspring disease risk. We recovered 62,864 nuclei of high quality from 10 samples each from the maternal-facing and fetal-facing sides of the placenta. On both sides in several cell types, MO was associated with upregulation of hypoxia response genes. On the maternal-facing side only, hypoxia gene expression was associated with offspring neurodevelopmental measures, in Gen3G, an independent pregnancy cohort with bulk placental tissue RNA-seq. We leveraged Gen3G to determine genes that correlated with impaired neurodevelopment and found these genes to be most highly expressed in extravillous trophoblasts (EVTs). EVTs further showed the strongest correlation between neurodevelopment impairment gene scores (NDIGSs) and the hypoxia gene score. We reanalyzed gene expression of cultured EVTs, and found increased NDIGSs associated with exposure to hypoxia. Among EVTs, accounting for the hypoxia gene score attenuated 44% of the association between BMI and NDIGSs. These data suggest that hypoxia in EVTs may be a key process in the neurodevelopmental programming of fetal exposure to MO.

Co-Transplantation-Based Human-Mouse Chimeric Brain Models to Study Human Glial-Glial and Glial-Neuronal Interactions

Human-mouse chimeric brain models, generated by transplanting human induced pluripotent stem cell (hiPSC)-derived neural cells, are valuable for studying the development and function of human neural cells in vivo. Understanding glial-glial and glial-neuronal interactions is essential for unraveling the complexities of brain function and developing treatments for neurological disorders. To explore these interactions between human neural cells within an intact brain environment, we employe a co-transplantation strategy involving the engraftment of hiPSC-derived neural progenitor cells along with primitive macrophage progenitors into the neonatal mouse brain. This approach creates human-mouse chimeric brains containing human microglia, macroglia (astroglia and oligodendroglia), and neurons. Using super-resolution imaging and 3D reconstruction techniques, we examine the dynamics between human neurons and glia, unveiling human microglia engulfing immature human neurons, microglia pruning synapses of human neurons, and significant interactions between human oligodendrocytes and neurons. Single-cell RNA sequencing analysis of the chimeric brain uncovers a close recapitulation of the human glial progenitor cell population, along with a dynamic stage in astroglial development that mirrors the processes found in the human brain. Furthermore, cell-cell communication analysis highlights significant neuronal-glial and glial-glial interactions, especially the interaction between adhesion molecules neurexins and neuroligins. This innovative co-transplantation model opens up new avenues for exploring the complex pathophysiological mechanisms underlying human neurological diseases. It holds particular promise for studying disorders where glial-neuronal interactions and non-cell-autonomous effects play crucial roles.

An atlas of transcribed enhancers across helper T cell diversity for decoding human diseases

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4+ T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

Single cell regulatory architecture of human pancreatic islets suggests sex differences in β cell function and the pathogenesis of type 2 diabetes

Type 2 and type 1 diabetes (T2D, T1D) exhibit sex differences in insulin secretion, the mechanisms of which are unknown. We examined sex differences in human pancreatic islets from 52 donors with and without T2D combining single cell RNA-seq (scRNA-seq), single nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq), hormone secretion, and bioenergetics. In nondiabetic (ND) donors, sex differences in islet cells gene accessibility and expression predominantly involved sex chromosomes. Islets from T2D donors exhibited similar sex differences in sex chromosomes differentially expressed genes (DEGs), but also exhibited sex differences in autosomal genes. Comparing β cells from T2D vs. ND donors, gene enrichment of female β cells showed suppression in mitochondrial respiration, while male β cells exhibited suppressed insulin secretion. Thus, although sex differences in gene accessibility and expression of ND β cells predominantly affect sex chromosomes, the transition to T2D reveals sex differences in autosomes highlighting mitochondrial failure in females.

Revealing the molecular landscape of human placenta: a systematic review and meta-analysis of single-cell RNA sequencing studies

BACKGROUND

With increasing significance of developmental programming effects associated with placental dysfunction, more investigations are devoted to improving the characterization and understanding of placental signatures in health and disease. The placenta is a transitory but dynamic organ adapting to the shifting demands of fetal development and available resources of the maternal supply throughout pregnancy. Trophoblasts (cytotrophoblasts, syncytiotrophoblasts, and extravillous trophoblasts) are placental-specific cell types responsible for the main placental exchanges and adaptations. Transcriptomic studies with single-cell resolution have led to advances in understanding the placenta's role in health and disease. These studies, however, often show discrepancies in characterization of the different placental cell types.

OBJECTIVE AND RATIONALE

We aim to review the knowledge regarding placental structure and function gained from the use of single-cell RNA sequencing (scRNAseq), followed by comparing cell-type-specific genes, highlighting their similarities and differences. Moreover, we intend to identify consensus marker genes for the various trophoblast cell types across studies. Finally, we will discuss the contributions and potential applications of scRNAseq in studying pregnancy-related diseases.

SEARCH METHODS

We conducted a comprehensive systematic literature review to identify different cell types and their functions at the human maternal-fetal interface, focusing on all original scRNAseq studies on placentas published before March 2023 and published reviews (total of 28 studies identified) using PubMed search. Our approach involved curating cell types and subtypes that had previously been defined using scRNAseq and comparing the genes used as markers or identified as potential new markers. Next, we reanalyzed expression matrices from the six available scRNAseq raw datasets with cell annotations (four from first trimester and two at term), using Wilcoxon rank-sum tests to compare gene expression among studies and annotate trophoblast cell markers in both first trimester and term placentas. Furthermore, we integrated scRNAseq raw data available from 18 healthy first trimester and nine term placentas, and performed clustering and differential gene expression analysis. We further compared markers obtained with the analysis of annotated and raw datasets with the literature to obtain a common signature gene list for major placental cell types.

OUTCOMES

Variations in the sampling site, gestational age, fetal sex, and subsequent sequencing and analysis methods were observed between the studies. Although their proportions varied, the three trophoblast types were consistently identified across all scRNAseq studies, unlike other non-trophoblast cell types. Notably, no marker genes were shared by all studies for any of the investigated cell types. Moreover, most of the newly defined markers in one study were not observed in other studies. These discrepancies were confirmed by our analysis on trophoblast cell types, where hundreds of potential marker genes were identified in each study but with little overlap across studies. From 35 461 and 23 378 cells of high quality in the first trimester and term placentas, respectively, we obtained major placental cell types, including perivascular cells that previously had not been identified in the first trimester. Importantly, our meta-analysis provides marker genes for major placental cell types based on our extensive curation.

WIDER IMPLICATIONS

This review and meta-analysis emphasizes the need for establishing a consensus for annotating placental cell types from scRNAseq data. The marker genes identified here can be deployed for defining human placental cell types, thereby facilitating and improving the reproducibility of trophoblast cell annotation.

Molecular maps of synovial cells in inflammatory arthritis using an optimized synovial tissue dissociation protocol

In this study, we optimized the dissociation of synovial tissue biopsies for single-cell omics studies and created a single-cell atlas of human synovium in inflammatory arthritis. The optimized protocol allowed consistent isolation of highly viable cells from tiny fresh synovial biopsies, minimizing the synovial biopsy drop-out rate. The synovium scRNA-seq atlas contained over 100,000 unsorted synovial cells from 25 synovial tissues affected by inflammatory arthritis, including 16 structural, 11 lymphoid, and 15 myeloid cell clusters. This synovial cell map expanded the diversity of synovial cell types/states, detected synovial neutrophils, and broadened synovial endothelial cell classification. We revealed tissue-resident macrophage subsets with proposed matrix-sensing (FOLR2+COLEC12high) and iron-recycling (LYVE1+SLC40A1+) activities and identified fibroblast subsets with proposed functions in cartilage breakdown (SOD2highSAA1+SAA2+SDC4+) and extracellular matrix remodeling (SERPINE1+COL5A3+LOXL2+). Our study offers an efficient synovium dissociation method and a reference scRNA-seq resource, that advances the current understanding of synovial cell heterogeneity in inflammatory arthritis.

Remodeling of anti-tumor immunity with antibodies targeting a p53 mutant

BACKGROUND

p53, the most frequently mutated gene in cancer, lacks effective targeted drugs.

METHODS

We developed monoclonal antibodies (mAbs) that target a p53 hotspot mutation E285K without cross-reactivity with wild-type p53. They were delivered using lipid nanoparticles (LNPs) that encapsulate DNA plasmids. Western blot, BLI, flow cytometry, single-cell sequencing (scRNA-seq), and other methods were employed to assess the function of mAbs in vitro and in vivo.

RESULTS

These LNP-pE285K-mAbs in the IgG1 format exhibited a robust anti-tumor effect, facilitating the infiltration of immune cells, including CD8+ T, B, and NK cells. scRNA-seq revealed that IgG1 reduces immune inhibitory signaling, increases MHC signaling from B cells to CD8+ T cells, and enriches anti-tumor T cell and B cell receptor profiles. The E285K-mAbs were also produced in the dimeric IgA (dIgA) format, whose anti-tumor activity depended on the polymeric immunoglobulin receptor (PIGR), a membrane Ig receptor, whereas that of IgG1 relied on TRIM21, an intracellular IgG receptor.

CONCLUSIONS

Targeting specific mutant epitopes using DNA-encoded and LNP-delivered mAbs represents a potential precision medicine strategy against p53 mutants in TRIM21- or PIGR-positive cancers.

CHIT1 at diagnosis predicts faster disability progression and reflects early microglial activation in multiple sclerosis

Multiple sclerosis (MS) is characterized by heterogeneity in disease course and prediction of long-term outcome remains a major challenge. Here, we investigate five myeloid markers - CHIT1, CHI3L1, sTREM2, GPNMB and CCL18 - in the cerebrospinal fluid (CSF) at diagnostic lumbar puncture in a longitudinal cohort of 192 MS patients. Through mixed-effects and machine learning models, we show that CHIT1 is a robust predictor for faster disability progression. Integrative analysis of 11 CSF and 26 central nervous system (CNS) parenchyma single-cell/nucleus RNA sequencing samples reveals CHIT1 to be predominantly expressed by microglia located in active MS lesions and enriched for lipid metabolism pathways. Furthermore, we find CHIT1 expression to accompany the transition from a homeostatic towards a more activated, MS-associated cell state in microglia. Neuropathological evaluation in post-mortem tissue from 12 MS patients confirms CHIT1 production by lipid-laden phagocytes in actively demyelinating lesions, already in early disease stages. Altogether, we provide a rationale for CHIT1 as an early biomarker for faster disability progression in MS.

Single cell regulatory architecture of human pancreatic islets suggests sex differences in β cell function and the pathogenesis of type 2 diabetes

Biological sex affects the pathogenesis of type 2 and type 1 diabetes (T2D, T1D) including the development of β cell failure observed more often in males. The mechanisms that drive sex differences in β cell failure is unknown. Studying sex differences in islet regulation and function represent a unique avenue to understand the sex-specific heterogeneity in β cell failure in diabetes. Here, we examined sex and race differences in human pancreatic islets from up to 52 donors with and without T2D (including 37 donors from the Human Pancreas Analysis Program [HPAP] dataset) using an orthogonal series of experiments including single cell RNA-seq (scRNA-seq), single nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq), dynamic hormone secretion, and bioenergetics. In cultured islets from nondiabetic (ND) donors, in the absence of the in vivo hormonal environment, sex differences in islet cell type gene accessibility and expression predominantly involved sex chromosomes. Of particular interest were sex differences in the X-linked KDM6A and Y-linked KDM5D chromatin remodelers in female and male islet cells respectively. Islets from T2D donors exhibited similar sex differences in differentially expressed genes (DEGs) from sex chromosomes. However, in contrast to islets from ND donors, islets from T2D donors exhibited major sex differences in DEGs from autosomes. Comparing β cells from T2D and ND donors revealed that females had more DEGs from autosomes compared to male β cells. Gene set enrichment analysis of female β cell DEGs showed a suppression of oxidative phosphorylation and electron transport chain pathways, while male β cell had suppressed insulin secretion pathways. Thus, although sex-specific differences in gene accessibility and expression of cultured ND human islets predominantly affect sex chromosome genes, major differences in autosomal gene expression between sexes appear during the transition to T2D and which highlight mitochondrial failure in female β cells.

KLRG1 Cell Depletion as a Novel Therapeutic Strategy in Patients with Mature T-Cell Lymphoma Subtypes

PURPOSE

Develop a novel therapeutic strategy for patients with subtypes of mature T-cell and NK-cell neoplasms.

EXPERIMENTAL DESIGN

Primary specimens, cell lines, patient-derived xenograft models, commercially available, and proprietary anti-KLRG1 antibodies were used for screening, target, and functional validation.

RESULTS

Here we demonstrate that surface KLRG1 is highly expressed on tumor cells in subsets of patients with extranodal NK/T-cell lymphoma (ENKTCL), T-prolymphocytic leukemia (T-PLL), and gamma/delta T-cell lymphoma (G/D TCL). The majority of the CD8+/CD57+ or CD3-/CD56+ leukemic cells derived from patients with T- and NK-large granular lymphocytic leukemia (T-LGLL and NK-LGLL), respectively, expressed surface KLRG1. The humanized afucosylated anti-KLRG1 monoclonal antibody (mAb208) optimized for mouse in vivo use depleted KLRG1+ TCL cells by mechanisms of ADCC, ADCP, and CDC rather than apoptosis. mAb208 induced ADCC and ADCP of T-LGLL patient-derived CD8+/CD57+ cells ex vivo. mAb208 effected ADCC of subsets of healthy donor-derived KLRG1+ NK, CD4+, CD8+ Tem, and TemRA cells while sparing KLRG1- naïve and CD8+ Tcm cells. Treatment of cell line and TCL patient-derived xenografts with mAb208 or anti-CD47 mAb alone and in combination with the PI3K-δ/γ inhibitor duvelisib extended survival. The depletion of macrophages in vivo antagonized mAb208 efficacy.

CONCLUSIONS

Our findings suggest the potential benefit of a broader treatment strategy combining therapeutic antibodies with PI3Ki for the treatment of patients with mature T-cell and NK-cell neoplasms. See related commentary by Varma and Diefenbach, p. 2300.

Computational identification of surface markers for isolating distinct subpopulations from heterogeneous cancer cell populations

Intratumor heterogeneity reduces treatment efficacy and complicates our understanding of tumor progression. There is a pressing need to understand the functions of heterogeneous tumor cell subpopulations within a tumor, yet biological systems to study these processes in vitro are limited. With the advent of single-cell RNA sequencing (scRNA-seq), it has become clear that some cancer cell line models include distinct subpopulations. Heterogeneous cell lines offer a unique opportunity to study the dynamics and evolution of genetically similar cancer cell subpopulations in controlled experimental settings. Here, we present clusterCleaver, a computational package that uses metrics of statistical distance to identify candidate surface markers maximally unique to transcriptomic subpopulations in scRNA-seq which may be used for FACS isolation. clusterCleaver was experimentally validated using the MDA-MB-231 and MDA-MB-436 breast cancer cell lines. ESAM and BST2/tetherin were experimentally confirmed as surface markers which identify and separate major transcriptomic subpopulations within MDA-MB-231 and MDA-MB-436 cells, respectively. clusterCleaver is a computationally efficient and experimentally validated workflow for identification and enrichment of distinct subpopulations within cell lines which paves the way for studies on the coexistence of cancer cell subpopulations in well-defined in vitro systems.

Single nucleus and spatial transcriptomic profiling of healthy human hamstring tendon

The molecular and cellular basis of health in human tendons remains poorly understood. Among human tendons, hamstring tendon has markedly low pathology and can provide a prototypic healthy tendon reference. The aim of this study was to determine the transcriptomes and location of all cell types in healthy hamstring tendon. Using single nucleus RNA sequencing, we profiled the transcriptomes of 10 533 nuclei from four healthy donors and identified 12 distinct cell types. We confirmed the presence of two fibroblast cell types, endothelial cells, mural cells, and immune cells, and identified cell types previously unreported in tendons, including different skeletal muscle cell types, satellite cells, adipocytes, and undefined nervous system cells. The location of these cell types within tendon was defined using spatial transcriptomics and imaging, and potential transcriptional networks and cell-cell interactions were analyzed. We demonstrate that fibroblasts have the highest number of potential cell-cell interactions in our dataset, are present throughout the tendon, and play an important role in the production and organization of extracellular matrix, thus confirming their role as key regulators of hamstring tendon homeostasis. Overall, our findings underscore the complexity of the cellular networks that underpin healthy human tendon function and the central role of fibroblasts as key regulators of hamstring tendon tissue homeostasis.

A single-nucleus transcriptomic atlas of medium spiny neurons in the rat nucleus accumbens

Neural processing of rewarding stimuli involves several distinct regions, including the nucleus accumbens (NAc). The majority of NAc neurons are GABAergic projection neurons known as medium spiny neurons (MSNs). MSNs are broadly defined by dopamine receptor expression, but evidence suggests that a wider array of subtypes exist. To study MSN heterogeneity, we analyzed single-nucleus RNA sequencing data from the largest available rat NAc dataset. Analysis of 48,040 NAc MSN nuclei identified major populations belonging to the striosome and matrix compartments. Integration with mouse and human data indicated consistency across species and disease-relevance scoring using genome-wide association study results revealed potentially differential roles for MSN populations in substance use disorders. Additional high-resolution clustering identified 34 transcriptomically distinct subtypes of MSNs definable by a limited number of marker genes. Together, these data demonstrate the diversity of MSNs in the NAc and provide a basis for more targeted genetic manipulation of specific populations.

Accelerating Single-Cell Sequencing Data Analysis with SciDAP: A User-Friendly Approach

Single-cell (sc) RNA, ATAC and Multiome sequencing became powerful tools for uncovering biological and disease mechanisms. Unfortunately, manual analysis of sc data presents multiple challenges due to large data volumes and complexity of configuration parameters. This complexity, as well as not being able to reproduce a computational environment, affects the reproducibility of analysis results. The Scientific Data Analysis Platform (https://SciDAP.com) allows biologists without computational expertise to analyze sequencing-based data using portable and reproducible pipelines written in Common Workflow Language (CWL). Our suite of computational pipelines addresses the most common needs in scRNA-Seq, scATAC-Seq and scMultiome data analysis. When executed on SciDAP, it offers a user-friendly alternative to manual data processing, eliminating the need for coding expertise. In this protocol, we describe the use of SciDAP to analyze scMultiome data. Similar approaches can be used for analysis of scRNA-Seq, scATAC-Seq and scVDJ-Seq datasets.

Decoding the complexity of delayed wound healing following Enterococcus faecalis infection

Wound infections are highly prevalent and can lead to delayed or failed healing, causing significant morbidity and adverse economic impacts. These infections occur in various contexts, including diabetic foot ulcers, burns, and surgical sites. Enterococcus faecalis is often found in persistent non-healing wounds, but its contribution to chronic wounds remains understudied. To address this, we employed single-cell RNA sequencing (scRNA-seq) on infected wounds in comparison to uninfected wounds in a mouse model. Examining over 23,000 cells, we created a comprehensive single-cell atlas that captures the cellular and transcriptomic landscape of these wounds. Our analysis revealed unique transcriptional and metabolic alterations in infected wounds, elucidating the distinct molecular changes associated with bacterial infection compared to the normal wound healing process. We identified dysregulated keratinocyte and fibroblast transcriptomes in response to infection, jointly contributing to an anti-inflammatory environment. Notably, E. faecalis infection prompted a premature, incomplete epithelial-mesenchymal transition in keratinocytes. Additionally, E. faecalis infection modulated M2-like macrophage polarization by inhibiting pro-inflammatory resolution in vitro, in vivo, and in our scRNA-seq atlas. Furthermore, we discovered macrophage crosstalk with neutrophils, which regulates chemokine signaling pathways, while promoting anti-inflammatory interactions with endothelial cells. Overall, our findings offer new insights into the immunosuppressive role of E. faecalis in wound infections.

Premature skewing of T cell receptor clonality and delayed memory expansion in HIV-exposed infants

While preventing vertical HIV transmission has been very successful, HIV-exposed uninfected infants (iHEU) experience an elevated risk to infections compared to HIV-unexposed and uninfected infants (iHUU). Here we present a longitudinal multimodal analysis of infant immune ontogeny that highlights the impact of HIV/ARV exposure. Using mass cytometry, we show alterations in T cell memory differentiation between iHEU and iHUU being significant from week 15 of life. The altered memory T cell differentiation in iHEU was preceded by lower TCR Vβ clonotypic diversity and linked to TCR clonal depletion within the naïve T cell compartment. Compared to iHUU, iHEU had elevated CD56loCD16loPerforin+CD38+CD45RA+FcεRIγ+ NK cells at 1 month postpartum and whose abundance pre-vaccination were predictive of vaccine-induced pertussis and rotavirus antibody responses post 3 months of life. Collectively, HIV/ARV exposure disrupted the trajectory of innate and adaptive immunity from birth which may underlie relative vulnerability to infections in iHEU.

A single dose of cocaine rewires the 3D genome structure of midbrain dopamine neurons

Midbrain dopamine neurons (DNs) respond to a first exposure to addictive drugs and play key roles in chronic drug usage1-3. As the synaptic and transcriptional changes that follow an acute cocaine exposure are mostly resolved within a few days4,5, the molecular changes that encode the long-term cellular memory of the exposure within DNs remain unknown. To investigate whether a single cocaine exposure induces long-term changes in the 3D genome structure of DNs, we applied Genome Architecture Mapping and single nucleus transcriptomic analyses in the mouse midbrain. We found extensive rewiring of 3D genome architecture at 24 hours past exposure which remains or worsens by 14 days, outlasting transcriptional responses. The cocaine-induced chromatin rewiring occurs at all genomic scales and affects genes with major roles in cocaine-induced synaptic changes. A single cocaine exposure triggers extensive long-lasting changes in chromatin condensation in post-synaptic and post-transcriptional regulatory genes, for example the unfolding of Rbfox1 which becomes most prominent 14 days post exposure. Finally, structurally remodeled genes are most expressed in a specific DN sub-type characterized by low expression of the dopamine auto-receptor Drd2, a key feature of highly cocaine-sensitive cells. These results reveal an important role for long-lasting 3D genome remodelling in the cellular memory of a single cocaine exposure, providing new hypotheses for understanding the inception of drug addiction and 3D genome plasticity.

A human neural crest model reveals the developmental impact of neuroblastoma-associated chromosomal aberrations

Early childhood tumours arise from transformed embryonic cells, which often carry large copy number alterations (CNA). However, it remains unclear how CNAs contribute to embryonic tumourigenesis due to a lack of suitable models. Here we employ female human embryonic stem cell (hESC) differentiation and single-cell transcriptome and epigenome analysis to assess the effects of chromosome 17q/1q gains, which are prevalent in the embryonal tumour neuroblastoma (NB). We show that CNAs impair the specification of trunk neural crest (NC) cells and their sympathoadrenal derivatives, the putative cells-of-origin of NB. This effect is exacerbated upon overexpression of MYCN, whose amplification co-occurs with CNAs in NB. Moreover, CNAs potentiate the pro-tumourigenic effects of MYCN and mutant NC cells resemble NB cells in tumours. These changes correlate with a stepwise aberration of developmental transcription factor networks. Together, our results sketch a mechanistic framework for the CNA-driven initiation of embryonal tumours.

Gene networks and systems biology in Alzheimer's disease: Insights from multi-omics approaches

Despite numerous studies in the field of dementia and Alzheimer's disease (AD), a comprehensive understanding of this devastating disease remains elusive. Bulk transcriptomics have provided insights into the underlying genetic factors at a high level. Subsequent technological advancements have focused on single-cell omics, encompassing techniques such as single-cell RNA sequencing and epigenomics, enabling the capture of RNA transcripts and chromatin states at a single cell or nucleus resolution. Furthermore, the emergence of spatial omics has allowed the study of gene responses in the vicinity of amyloid beta plaques or across various brain regions. With the vast amount of data generated, utilizing gene regulatory networks to comprehensively study this disease has become essential. This review delves into some techniques employed in the field of AD, explores the discoveries made using these techniques, and provides insights into the future of the field.

Characterizing cell type specific transcriptional differences between the living and postmortem human brain

Single-nucleus RNA sequencing (snRNA-seq) is often used to define gene expression patterns characteristic of brain cell types as well as to identify cell type specific gene expression signatures of neurological and mental illnesses in postmortem human brains. As methods to obtain brain tissue from living individuals emerge, it is essential to characterize gene expression differences associated with tissue originating from either living or postmortem subjects using snRNA-seq, and to assess whether and how such differences may impact snRNA-seq studies of brain tissue. To address this, human prefrontal cortex single nuclei gene expression was generated and compared between 31 samples from living individuals and 21 postmortem samples. The same cell types were consistently identified in living and postmortem nuclei, though for each cell type, a large proportion of genes were differentially expressed between samples from postmortem and living individuals. Notably, estimation of cell type proportions by cell type deconvolution of pseudo-bulk data was found to be more accurate in samples from living individuals. To allow for future integration of living and postmortem brain gene expression, a model was developed that quantifies from gene expression data the probability a human brain tissue sample was obtained postmortem. These probabilities are established as a means to statistically account for the gene expression differences between samples from living and postmortem individuals. Together, the results presented here provide a deep characterization of both differences between snRNA-seq derived from samples from living and postmortem individuals, as well as qualify and account for their effect on common analyses performed on this type of data.

A cell atlas of thoracic aortic perivascular adipose tissue: a focus on mechanotransducers

Perivascular adipose tissue (PVAT) is increasingly recognized for its function in mechanotransduction. However, major gaps remain in our understanding of the cells present in PVAT, as well as how different cells contribute to mechanotransduction. We hypothesized that snRNA-seq would reveal the expression of mechanotransducers, and test one (PIEZO1) to illustrate the expression and functional agreement between single-nuclei RNA sequencing (snRNA-seq) and physiological measurements. To contrast two brown tissues, subscapular brown adipose tissue (BAT) was also examined. We used snRNA-seq of the thoracic aorta PVAT (taPVAT) and BAT from male Dahl salt-sensitive (Dahl SS) rats to investigate cell-specific expression mechanotransducers. Localization and function of the mechanostransducer PIEZO1 were further examined using immunohistochemistry (IHC) and RNAscope, as well as pharmacological antagonism. Approximately 30,000 nuclei from taPVAT and BAT each were characterized by snRNA-seq, identifying eight major cell types expected and one unexpected (nuclei with oligodendrocyte marker genes). Cell-specific differential gene expression analysis between taPVAT and BAT identified up to 511 genes (adipocytes) with many (≥20%) being unique to individual cell types. Piezo1 was the most highly, widely expressed mechanotransducer. The presence of PIEZO1 in the PVAT but not the adventitia was confirmed by RNAscope and IHC in male and female rats. Importantly, antagonism of PIEZO1 by GsMTX4 impaired the PVAT's ability to hold tension. Collectively, the cell compositions of taPVAT and BAT are highly similar, and PIEZO1 is likely a mechanotransducer in taPVAT.NEW & NOTEWORTHY This study describes the atlas of cells in the thoracic aorta perivascular adipose tissue (taPVAT) of the Dahl-SS rat, an important hypertension model. We show that mechanotransducers are widely expressed in these cells. Moreover, PIEZO1 expression is shown to be restricted to the taPVAT and is functionally implicated in stress relaxation. These data will serve as the foundation for future studies investigating the role of taPVAT in this model of hypertensive disease.

Single-nucleus and spatial transcriptomic profiling of human temporal cortex and white matter reveals novel associations with AD pathology

Alzheimer's disease (AD) is a neurodegenerative disorder with complex pathological manifestations and is the leading cause of cognitive decline and dementia in elderly individuals. A major goal in AD research is to identify new therapeutic pathways by studying the molecular and cellular changes in the disease, either downstream or upstream of the pathological hallmarks. In this study, we present a comprehensive investigation of cellular heterogeneity from the temporal cortex region of 40 individuals, comprising healthy donors and individuals with differing tau and amyloid burden. Using single-nucleus transcriptome analysis of 430,271 nuclei from both gray and white matter of these individuals, we identified cell type-specific subclusters in both neuronal and glial cell types with varying degrees of association with AD pathology. In particular, these associations are present in layer specific glutamatergic (excitatory) neuronal types, along with GABAergic (inhibitory) neurons and glial subtypes. These associations were observed in early as well as late pathological progression. We extended this analysis by performing multiplexed in situ hybridization using the CARTANA platform, capturing 155 genes in 13 individuals with varying levels of tau pathology. By modeling the spatial distribution of these genes and their associations with the pathology, we not only replicated key findings from our snRNA data analysis, but also identified a set of cell type-specific genes that show selective enrichment or depletion near pathological inclusions. Together, our findings allow us to prioritize specific cell types and pathways for targeted interventions at various stages of pathological progression in AD.

Single-nucleus and spatial landscape of the sub-ventricular zone in human glioblastoma

The sub-ventricular zone (SVZ) is the most well-characterized neurogenic area in the mammalian brain. We previously showed that in 65% of patients with glioblastoma (GBM), the SVZ is a reservoir of cancer stem-like cells that contribute to treatment resistance and emergence of recurrence. Here, we built a single-nucleus RNA-sequencing-based microenvironment landscape of the tumor mass (T_Mass) and the SVZ (T_SVZ) of 15 GBM patients and 2 histologically normal SVZ (N_SVZ) samples as controls. We identified a mesenchymal signature in the T_SVZ of GBM patients: tumor cells from the T_SVZ relied on the ZEB1 regulatory network, whereas tumor cells in the T_Mass relied on the TEAD1 regulatory network. Moreover, the T_SVZ microenvironment was predominantly characterized by tumor-supportive microglia, which spatially co-exist and establish heterotypic interactions with tumor cells. Lastly, differential gene expression analyses, predictions of ligand-receptor and incoming/outgoing interactions, and functional assays revealed that the IL-1β/IL-1RAcP and Wnt-5a/Frizzled-3 pathways are therapeutic targets in the T_SVZ microenvironment. Our data provide insights into the biology of the SVZ in GBM patients and identify specific targets of this microenvironment.

Single-cell multiomic dissection of response and resistance to chimeric antigen receptor T cells against BCMA in relapsed multiple myeloma

Markers that predict response and resistance to chimeric antigen receptor (CAR) T cells in relapsed/refractory multiple myeloma are currently missing. We subjected mononuclear cells isolated from peripheral blood and bone marrow before and after the application of approved B cell maturation antigen-directed CAR T cells to single-cell multiomic analyses to identify markers associated with resistance and early relapse. Differences between responders and nonresponders were identified at the time of leukapheresis. Nonresponders showed an immunosuppressive microenvironment characterized by increased numbers of monocytes expressing the immune checkpoint molecule CD39 and suppressed CD8+ T cell and natural killer cell function. Analysis of CAR T cells showed cytotoxic and exhausted phenotypes in hyperexpanded clones compared to low/intermediate expanded clones. We identified potential immunotherapy targets on CAR T cells, like PD1, to improve their functionality and durability. Our work provides evidence that an immunosuppressive microenvironment causes resistance to CAR T cell therapies in multiple myeloma.

Identification of receptors and factors associated with human coronaviruses in the oral cavity using single-cell RNA sequencing

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) ravaged the world, and Coronavirus Disease 2019 (COVID-19) exhibited highly prevalent oral symptoms that had significantly impacted the lives of affected patients. However, the involvement of four human coronavirus (HCoVs), namely SARS-CoV-2, SARS-CoV, MERS-CoV, and HCoV-229E, in oral cavity infections remained poorly understood. We integrated single-cell RNA sequencing (scRNA-seq) data of seven human oral tissues through consistent normalization procedure, including minor salivary gland (MSG), parotid gland (PG), tongue, gingiva, buccal, periodontium and pulp. The Seurat, scDblFinder, Harmony, SingleR, Ucell and scCancer packages were comprehensively used for analysis. We identified specific cell clusters and generated expression profiles of SARS-CoV-2 and coronavirus-associated receptors and factors (SCARFs) in seven oral regions, providing direction for predicting the tropism of four HCoVs for oral tissues, as well as for dental clinical treatment. Based on our analysis, it appears that various SCARFs, including ACE2, ASGR1, KREMEN1, DPP4, ANPEP, CD209, CLEC4G/M, TMPRSS family proteins (including TMPRSS2, TMPRSS4, and TMPRSS11A), and FURIN, are expressed at low levels in the oral cavity. Conversely, BSG, CTSB, and CTSL exhibit enrichment in oral tissues. Our study also demonstrates widespread expression of restriction factors, particularly IFITM1-3 and LY6E, in oral cells. Additionally, some replication, assembly, and trafficking factors appear to exhibit broad oral tissues expression patterns. Overall, the oral cavity could potentially serve as a high-risk site for SARS-CoV-2 infection, while displaying a comparatively lower degree of susceptibility towards other HCoVs (including SARS-CoV, MERS-CoV and HCoV-229E). Specifically, MSG, tongue, and gingiva represent potential sites of vulnerability for four HCoVs infection, with the MSG exhibiting a particularly high susceptibility. However, the expression patterns of SCARFs in other oral sites demonstrate relatively intricate and may only be specifically associated with SARS-CoV-2 infection. Our study sheds light on the mechanisms of HCoVs infection in the oral cavity as well as gains insight into the characteristics and distribution of possible HCoVs target cells in oral tissues, providing potential therapeutic targets for HCoVs infection in the oral cavity.

Demuxafy: improvement in droplet assignment by integrating multiple single-cell demultiplexing and doublet detection methods

Recent innovations in single-cell RNA-sequencing (scRNA-seq) provide the technology to investigate biological questions at cellular resolution. Pooling cells from multiple individuals has become a common strategy, and droplets can subsequently be assigned to a specific individual by leveraging their inherent genetic differences. An implicit challenge with scRNA-seq is the occurrence of doublets-droplets containing two or more cells. We develop Demuxafy, a framework to enhance donor assignment and doublet removal through the consensus intersection of multiple demultiplexing and doublet detecting methods. Demuxafy significantly improves droplet assignment by separating singlets from doublets and classifying the correct individual.

scRNA-seq characterizing the heterogeneity of fibroblasts in breast cancer reveals a novel subtype SFRP4+ CAF that inhibits migration and predicts prognosis

Introduction

Cancer-associated fibroblasts (CAFs) are a diverse group of cells that significantly impact the tumor microenvironment and therapeutic responses in breast cancer (BC). Despite their importance, the comprehensive profile of CAFs in BC remains to be fully elucidated.

Methods

To address this gap, we utilized single-cell RNA sequencing (scRNA-seq) to delineate the CAF landscape within 14 BC normal-tumor paired samples. We further corroborated our findings by analyzing several public datasets, thereby validating the newly identified CAF subtype. Additionally, we conducted coculture experiments with BC cells to assess the functional implications of this CAF subtype.

Results

Our scRNA-seq analysis unveiled eight distinct CAF subtypes across five tumor and six adjacent normal tissue samples. Notably, we discovered a novel subtype, designated as SFRP4+ CAFs, which was predominantly observed in normal tissues. The presence of SFRP4+ CAFs was substantiated by two independent scRNA-seq datasets and a spatial transcriptomics dataset. Functionally, SFRP4+ CAFs were found to impede BC cell migration and the epithelial-mesenchymal transition (EMT) process by secreting SFRP4, thereby modulating the WNT signaling pathway. Furthermore, we established that elevated expression levels of SFRP4+ CAF markers correlate with improved survival outcomes in BC patients, yet paradoxically, they predict a diminished response to neoadjuvant chemotherapy in cases of triple-negative breast cancer.

Conclusion

This investigation sheds light on the heterogeneity of CAFs in BC and introduces a novel SFRP4+ CAF subtype that hinders BC cell migration. This discovery holds promise as a potential biomarker for refined prognostic assessment and therapeutic intervention in BC.

Respiratory viral infection promotes the awakening and outgrowth of dormant metastatic breast cancer cells in lungs

Breast cancer is the second most common cancer globally. Most deaths from breast cancer are due to metastatic disease which often follows long periods of clinical dormancy1. Understanding the mechanisms that disrupt the quiescence of dormant disseminated cancer cells (DCC) is crucial for addressing metastatic progression. Infection with respiratory viruses (e.g. influenza or SARS-CoV-2) is common and triggers an inflammatory response locally and systemically2,3. Here we show that influenza virus infection leads to loss of the pro-dormancy mesenchymal phenotype in breast DCC in the lung, causing DCC proliferation within days of infection, and a greater than 100-fold expansion of carcinoma cells into metastatic lesions within two weeks. Such DCC phenotypic change and expansion is interleukin-6 (IL-6)-dependent. We further show that CD4 T cells are required for the maintenance of pulmonary metastatic burden post-influenza virus infection, in part through attenuation of CD8 cell responses in the lungs. Single-cell RNA-seq analyses reveal DCC-dependent impairment of T-cell activation in the lungs of infected mice. SARS-CoV-2 infected mice also showed increased breast DCC expansion in lungs post-infection. Expanding our findings to human observational data, we observed that cancer survivors contracting a SARS-CoV-2 infection have substantially increased risks of lung metastatic progression and cancer-related death compared to cancer survivors who did not. These discoveries underscore the significant impact of respiratory viral infections on the resurgence of metastatic cancer, offering novel insights into the interconnection between infectious diseases and cancer metastasis.

Dissecting the spatiotemporal diversity of adult neural stem cells

Adult stem cells are important for tissue turnover and regeneration. However, in most adult systems it remains elusive how stem cells assume different functional states and support spatially patterned tissue architecture. Here, we dissected the diversity of neural stem cells in the adult zebrafish brain, an organ that is characterized by pronounced zonation and high regenerative capacity. We combined single-cell transcriptomics of dissected brain regions with massively parallel lineage tracing and in vivo RNA metabolic labeling to analyze the regulation of neural stem cells in space and time. We detected a large diversity of neural stem cells, with some subtypes being restricted to a single brain region, while others were found globally across the brain. Global stem cell states are linked to neurogenic differentiation, with different states being involved in proliferative and non-proliferative differentiation. Our work reveals principles of adult stem cell organization and establishes a resource for the functional manipulation of neural stem cell subtypes.

De novo antibody discovery in human blood from full-length single B cell transcriptomics and matching haplotyped-resolved germline assemblies

Immunoglobulin (IGH, IGK, IGL) loci in the human genome are highly polymorphic regions that encode the building blocks of the light and heavy chain IG proteins that dimerize to form antibodies. The processes of V(D)J recombination and somatic hypermutation in B cells are responsible for creating an enormous reservoir of highly specific antibodies capable of binding a vast array of possible antigens. However, the antibody repertoire is fundamentally limited by the set of variable (V), diversity (D), and joining (J) alleles present in the germline IG loci. To better understand how the germline IG haplotypes contribute to the expressed antibody repertoire, we combined genome sequencing of the germline IG loci with single-cell transcriptome sequencing of B cells from the same donor. Sequencing and assembly of the germline IG loci captured the IGH locus in a single fully-phased contig where the maternal and paternal contributions to the germline V, D, and J repertoire can be fully resolved. The B cells were collected following a measles, mumps, and rubella (MMR) vaccination, resulting in a population of cells that were activated in response to this specific immune challenge. Single-cell, full-length transcriptome sequencing of these B cells resulted in whole transcriptome characterization of each cell, as well as highly-accurate consensus sequences for the somatically rearranged and hypermutated light and heavy chain IG transcripts. A subset of antibodies synthesized based on their consensus heavy and light chain transcript sequences demonstrated binding to measles antigens and neutralization of measles live virus.

Repetitive head impacts induce neuronal loss and neuroinflammation in young athletes

Repetitive head impacts (RHI) sustained from contact sports are the largest risk factor for chronic traumatic encephalopathy (CTE). Currently, CTE can only be diagnosed after death and the multicellular cascade of events that trigger initial hyperphosphorylated tau (p-tau) deposition remain unclear. Further, the symptoms endorsed by young individuals with early disease are not fully explained by the extent of p-tau deposition, severely hampering development of therapeutic interventions. Here, we show that RHI exposure associates with a multicellular response in young individuals (<51 years old) prior to the onset of CTE p-tau pathology that correlates with number of years of RHI exposure. Leveraging single nucleus RNA sequencing of tissue from 8 control, 9 RHI-exposed, and 11 low stage CTE individuals, we identify SPP1+ inflammatory microglia, angiogenic and inflamed endothelial cell profiles, reactive astrocytes, and altered synaptic gene expression in excitatory and inhibitory neurons in all individuals with exposure to RHI. Surprisingly, we also observe a significant loss of cortical sulcus layer 2/3 neurons in contact sport athletes compared to controls independent of p-tau pathology. These results provide robust evidence that multiple years of RHI exposure is sufficient to induce lasting cellular alterations that may underlie p-tau deposition and help explain the early clinical symptoms observed in young former contact sport athletes. Furthermore, these data identify specific cellular responses to repetitive head impacts that may direct future identification of diagnostic and therapeutic strategies for CTE.

CDK12 Loss Promotes Prostate Cancer Development While Exposing Vulnerabilities to Paralog-Based Synthetic Lethality

Biallelic loss of cyclin-dependent kinase 12 (CDK12) defines a unique molecular subtype of metastatic castration-resistant prostate cancer (mCRPC). It remains unclear, however, whether CDK12 loss per se is sufficient to drive prostate cancer development-either alone, or in the context of other genetic alterations-and whether CDK12-mutant tumors exhibit sensitivity to specific pharmacotherapies. Here, we demonstrate that tissue-specific Cdk12 ablation is sufficient to induce preneoplastic lesions and robust T cell infiltration in the mouse prostate. Allograft-based CRISPR screening demonstrated that Cdk12 loss is positively associated with Trp53 inactivation but negatively associated with Pten inactivation-akin to what is observed in human mCRPC. Consistent with this, ablation of Cdk12 in prostate organoids with concurrent Trp53 loss promotes their proliferation and ability to form tumors in mice, while Cdk12 knockout in the Pten-null prostate cancer mouse model abrogates tumor growth. Bigenic Cdk12 and Trp53 loss allografts represent a new syngeneic model for the study of androgen receptor (AR)-positive, luminal prostate cancer. Notably, Cdk12/Trp53 loss prostate tumors are sensitive to immune checkpoint blockade. Cdk12-null organoids (either with or without Trp53 co-ablation) and patient-derived xenografts from tumors with CDK12 inactivation are highly sensitive to inhibition or degradation of its paralog kinase, CDK13. Together, these data identify CDK12 as a bona fide tumor suppressor gene with impact on tumor progression and lends support to paralog-based synthetic lethality as a promising strategy for treating CDK12-mutant mCRPC.

Whole organism snRNA-seq reveals systemic peripheral changes in Alzheimer's Disease fly models

Peripheral tissues become disrupted in Alzheimer's Disease (AD). However, a comprehensive understanding of how the expression of AD-associated toxic proteins, Aβ42 and Tau, in neurons impacts the periphery is lacking. Using Drosophila, a prime model organism for studying aging and neurodegeneration, we generated the Alzheimer's Disease Fly Cell Atlas (AD-FCA): whole-organism single-nucleus transcriptomes of 219 cell types from adult flies neuronally expressing human Aβ42 or Tau. In-depth analyses and functional data reveal impacts on peripheral sensory neurons by Aβ42 and on various non-neuronal peripheral tissues by Tau, including the gut, fat body, and reproductive system. This novel AD atlas provides valuable insights into potential biomarkers and the intricate interplay between the nervous system and peripheral tissues in response to AD-associated proteins.

Dendritic cell-targeted therapy expands CD8 T cell responses to bona-fide neoantigens in lung tumors

Cross-presentation by type 1 DCs (cDC1) is critical to induce and sustain antitumoral CD8 T cell responses to model antigens, in various tumor settings. However, the impact of cross-presenting cDC1 and the potential of DC-based therapies in tumors carrying varied levels of bona-fide neoantigens (neoAgs) remain unclear. Here we develop a hypermutated model of non-small cell lung cancer in female mice, encoding genuine MHC-I neoepitopes to study neoAgs-specific CD8 T cell responses in spontaneous settings and upon Flt3L + αCD40 (DC-therapy). We find that cDC1 are required to generate broad CD8 responses against a range of diverse neoAgs. DC-therapy promotes immunogenicity of weaker neoAgs and strongly inhibits the growth of high tumor-mutational burden (TMB) tumors. In contrast, low TMB tumors respond poorly to DC-therapy, generating mild CD8 T cell responses that are not sufficient to block progression. scRNA transcriptional analysis, immune profiling and functional assays unveil the changes induced by DC-therapy in lung tissues, which comprise accumulation of cDC1 with increased immunostimulatory properties and less exhausted effector CD8 T cells. We conclude that boosting cDC1 activity is critical to broaden the diversity of anti-tumoral CD8 T cell responses and to leverage neoAgs content for therapeutic advantage.

Single-cell division tracing and transcriptomics reveal cell types and differentiation paths in the regenerating lung

Understanding the molecular and cellular processes involved in lung epithelial regeneration may fuel the development of therapeutic approaches for lung diseases. We combine mouse models allowing diphtheria toxin-mediated damage of specific epithelial cell types and parallel GFP-labeling of functionally dividing cells with single-cell transcriptomics to characterize the regeneration of the distal lung. We uncover cell types, including Krt13+ basal and Krt15+ club cells, detect an intermediate cell state between basal and goblet cells, reveal goblet cells as actively dividing progenitor cells, and provide evidence that adventitial fibroblasts act as supporting cells in epithelial regeneration. We also show that diphtheria toxin-expressing cells can persist in the lung, express specific inflammatory factors, and transcriptionally resemble a previously undescribed population in the lungs of COVID-19 patients. Our study provides a comprehensive single-cell atlas of the distal lung that characterizes early transcriptional and cellular responses to concise epithelial injury, encompassing proliferation, differentiation, and cell-to-cell interactions.

ScLinear predicts protein abundance at single-cell resolution

Single-cell multi-omics have transformed biomedical research and present exciting machine learning opportunities. We present scLinear, a linear regression-based approach that predicts single-cell protein abundance based on RNA expression. ScLinear is vastly more efficient than state-of-the-art methodologies, without compromising its accuracy. ScLinear is interpretable and accurately generalizes in unseen single-cell and spatial transcriptomics data. Importantly, we offer a critical view in using complex algorithms ignoring simpler, faster, and more efficient approaches.

Autophagy promotes efficient T cell responses to restrict high-dose Mycobacterium tuberculosis infection in mice

Although autophagy sequesters Mycobacterium tuberculosis (Mtb) in in vitro cultured macrophages, loss of autophagy in macrophages in vivo does not result in susceptibility to a standard low-dose Mtb infection until late during infection, leaving open questions regarding the protective role of autophagy during Mtb infection. Here we report that loss of autophagy in lung macrophages and dendritic cells results in acute susceptibility of mice to high-dose Mtb infection, a model mimicking active tuberculosis. Rather than observing a role for autophagy in controlling Mtb replication in macrophages, we find that autophagy suppresses macrophage responses to Mtb that otherwise result in accumulation of myeloid-derived suppressor cells and subsequent defects in T cell responses. Our finding that the pathogen-plus-susceptibility gene interaction is dependent on dose has important implications both for understanding how Mtb infections in humans lead to a spectrum of outcomes and for the potential use of autophagy modulators in clinical medicine.

The single-cell transcriptome program of nodule development cellular lineages in Medicago truncatula

Legumes establish a symbiotic relationship with nitrogen-fixing rhizobia by developing nodules. Nodules are modified lateral roots that undergo changes in their cellular development in response to bacteria, but the transcriptional reprogramming that occurs in these root cells remains largely uncharacterized. Here, we describe the cell-type-specific transcriptome response of Medicago truncatula roots to rhizobia during early nodule development in the wild-type genotype Jemalong A17, complemented with a hypernodulating mutant (sunn-4) to expand the cell population responding to infection and subsequent biological inferences. The analysis identifies epidermal root hair and stele sub-cell types associated with a symbiotic response to infection and regulation of nodule proliferation. Trajectory inference shows cortex-derived cell lineages differentiating to form the nodule primordia and, posteriorly, its meristem, while modulating the regulation of phytohormone-related genes. Gene regulatory analysis of the cell transcriptomes identifies new regulators of nodulation, including STYLISH 4, for which the function is validated.

ISL1 and POU4F1 Directly Interact to Regulate the Differentiation and Survival of Inner Ear Sensory Neurons

The inner ear sensory neurons play a pivotal role in auditory processing and balance control. Though significant progresses have been made, the underlying mechanisms controlling the differentiation and survival of the inner ear sensory neurons remain largely unknown. During development, ISL1 and POU4F transcription factors are co-expressed and are required for terminal differentiation, pathfinding, axon outgrowth and the survival of neurons in the central and peripheral nervous systems. However, little is understood about their functional relationship and regulatory mechanism in neural development. Here, we have knocked out Isl1 or Pou4f1 or both in mice of both sexes. In the absence of Isl1, the differentiation of cochleovestibular ganglion (CVG) neurons is disturbed and with that Isl1-deficient CVG neurons display defects in migration and axon pathfinding. Compound deletion of Isl1 and Pou4f1 causes a delay in CVG differentiation and results in a more severe CVG defect with a loss of nearly all of spiral ganglion neurons (SGNs). Moreover, ISL1 and POU4F1 interact directly in developing CVG neurons and act cooperatively as well as independently in regulating the expression of unique sets of CVG-specific genes crucial for CVG development and survival by binding to the cis-regulatory elements including the promoters of Fgf10, Pou4f2, and Epha5 and enhancers of Eya1 and Ntng2 These findings demonstrate that Isl1 and Pou4f1 are indispensable for CVG development and maintenance by acting epistatically to regulate genes essential for CVG development.

Massive-scale single-cell chromatin accessibility sequencing using combinatorial fluidic indexing

Single-cell ATAC-seq has emerged as a powerful approach for revealing candidate cis-regulatory elements genome-wide at cell-type resolution. However, current single-cell methods suffer from limited throughput and high costs. Here, we present a novel technique called single-cell combinatorial fluidic indexing ATAC-sequencing ("scifi-ATAC-seq"), which combines a barcoded Tn5 pre-indexing step with droplet-based single-cell ATAC-seq using a widely commercialized microfluidics platform (10X Genomics). With scifi-ATAC-seq, up to 200,000 nuclei across multiple samples in a single emulsion reaction can be indexed, representing a ~20-fold increase in throughput compared to the standard 10X Genomics workflow.

Multiomic single-cell sequencing defines tissue-specific responses in Stevens-Johnson Syndrome and Toxic epidermal necrolysis

Stevens-Johnson syndrome and toxic epidermal necrolysis (SJS/TEN) is a rare but life-threatening cutaneous drug reaction mediated by human leukocyte antigen (HLA) class I-restricted CD8+ T-cells. To obtain an unbiased assessment of SJS/TEN cellular immunopathogenesis, we performed single-cell (sc) transcriptome, surface proteome, and TCR sequencing on unaffected skin, affected skin, and blister fluid from 17 SJS/TEN patients. From 119,784 total cells, we identified 16 scRNA-defined subsets, confirmed by subset-defining surface protein expression. Keratinocytes upregulated HLA and IFN-response genes in the affected skin. Cytotoxic CD8+ T-cell subpopulations of expanded and unexpanded TCRαβ clonotypes were shared in affected skin and blister fluid but absent or unexpanded in SJS/TEN unaffected skin. SJS/TEN blister fluid is a rich reservoir of oligoclonal CD8+ T-cells with an effector phenotype driving SJS/TEN pathogenesis. This multiomic database will act as the basis to define antigen-reactivity, HLA restriction, and signatures of drug-antigen-reactive T-cell clonotypes at a tissue level.

An atlas of cells in the human tonsil

Palatine tonsils are secondary lymphoid organs (SLOs) representing the first line of immunological defense against inhaled or ingested pathogens. We generated an atlas of the human tonsil composed of >556,000 cells profiled across five different data modalities, including single-cell transcriptome, epigenome, proteome, and immune repertoire sequencing, as well as spatial transcriptomics. This census identified 121 cell types and states, defined developmental trajectories, and enabled an understanding of the functional units of the tonsil. Exemplarily, we stratified myeloid slan-like subtypes, established a BCL6 enhancer as locally active in follicle-associated T and B cells, and identified SIX5 as putative transcriptional regulator of plasma cell maturation. Analyses of a validation cohort confirmed the presence, annotation, and markers of tonsillar cell types and provided evidence of age-related compositional shifts. We demonstrate the value of this resource by annotating cells from B cell-derived mantle cell lymphomas, linking transcriptional heterogeneity to normal B cell differentiation states of the human tonsil.

Cell atlas of the Atlantic salmon spleen reveals immune cell heterogeneity and cell-specific responses to bacterial infection

The spleen is a conserved secondary lymphoid organ that emerged in parallel to adaptive immunity in early jawed vertebrates. Recent studies have applied single cell transcriptomics to reveal the cellular composition of spleen in several species, cataloguing diverse immune cell types and subpopulations. In this study, 51,119 spleen nuclei transcriptomes were comprehensively investigated in the commercially important teleost Atlantic salmon (Salmo salar L.), contrasting control animals with those challenged with the bacterial pathogen Aeromonas salmonicida. We identified clusters of nuclei representing the expected major cell types, namely T cells, B cells, natural killer-like cells, granulocytes, mononuclear phagocytes, endothelial cells, mesenchymal cells, erythrocytes and thrombocytes. We discovered heterogeneity within several immune lineages, providing evidence for resident macrophages and melanomacrophages, infiltrating monocytes, several candidate dendritic cell subpopulations, and B cells at distinct stages of differentiation, including plasma cells and an igt + subset. We provide evidence for twelve candidate T cell subsets, including cd4+ T helper and regulatory T cells, one cd8+ subset, three γδT subsets, and populations double negative for cd4 and cd8. The number of genes showing differential expression during the early stages of Aeromonas infection was highly variable across immune cell types, with the largest changes observed in macrophages and infiltrating monocytes, followed by resting mature B cells. Our analysis provides evidence for a local inflammatory response to infection alongside B cell maturation in the spleen, and upregulation of ccr9 genes in igt + B cells, T helper and cd8+ cells, and monocytes, consistent with the recruitment of immune cell populations to the gut to deal with Aeromonas infection. Overall, this study provides a new cell-resolved perspective of the immune actions of Atlantic salmon spleen, highlighting extensive heterogeneity hidden to bulk transcriptomics. We further provide a large catalogue of cell-specific marker genes that can be leveraged to further explore the function and structural organization of the salmonid immune system.