The single-cell transcriptional landscape of mammalian organogenesis

A multi-omic single-cell landscape of cellular diversification in the developing human cerebral cortex

The vast neuronal diversity in the human neocortex is vital for high-order brain functions, necessitating elucidation of the regulatory mechanisms underlying such unparalleled diversity. However, recent studies have yet to comprehensively reveal the diversity of neurons and the molecular logic of neocortical origin in humans at single-cell resolution through profiling transcriptomic or epigenomic landscapes, owing to the application of unimodal data alone to depict exceedingly heterogeneous populations of neurons. In this study, we generated a comprehensive compendium of the developing human neocortex by simultaneously profiling gene expression and open chromatin from the same cell. We computationally reconstructed the differentiation trajectories of excitatory projection neurons of cortical origin and inferred the regulatory logic governing lineage bifurcation decisions for neuronal diversification. We demonstrated that neuronal diversity arises from progenitor cell lineage specificity and postmitotic differentiation at distinct stages. Our data paves the way for understanding the primarily coordinated regulatory logic for neuronal diversification in the neocortex.

Single-cell RNA sequencing analysis of the retina under acute high intraocular pressure

JOURNAL/nrgr/04.03/01300535-202419110-00032/figure1/v/2024-03-08T184507Z/r/image-tiff High intraocular pressure causes retinal ganglion cell injury in primary and secondary glaucoma diseases, yet the molecular landscape characteristics of retinal cells under high intraocular pressure remain unknown. Rat models of acute hypertension ocular pressure were established by injection of cross-linked hyaluronic acid hydrogel (Healaflow®). Single-cell RNA sequencing was then used to describe the cellular composition and molecular profile of the retina following high intraocular pressure. Our results identified a total of 12 cell types, namely retinal pigment epithelial cells, rod-photoreceptor cells, bipolar cells, Müller cells, microglia, cone-photoreceptor cells, retinal ganglion cells, endothelial cells, retinal progenitor cells, oligodendrocytes, pericytes, and fibroblasts. The single-cell RNA sequencing analysis of the retina under acute high intraocular pressure revealed obvious changes in the proportions of various retinal cells, with ganglion cells decreased by 23%. Hematoxylin and eosin staining and TUNEL staining confirmed the damage to retinal ganglion cells under high intraocular pressure. We extracted data from retinal ganglion cells and analyzed the retinal ganglion cell cluster with the most distinct expression. We found upregulation of the B3gat2 gene, which is associated with neuronal migration and adhesion, and downregulation of the Tsc22d gene, which participates in inhibition of inflammation. This study is the first to reveal molecular changes and intercellular interactions in the retina under high intraocular pressure. These data contribute to understanding of the molecular mechanism of retinal injury induced by high intraocular pressure and will benefit the development of novel therapies.

Stromal Rigidity Stress Accelerates Pancreatic Intraepithelial Neoplasia Progression and Chromosomal Instability via Nuclear Protein Tyrosine Kinase 2 Localization

Because the mechanotransduction by stromal stiffness stimulates the rupture and repair of the nuclear envelope in pancreatic progenitor cells, accumulated genomic aberrations are under selection in the tumor microenvironment. Analysis of cell growth, micronuclei, and phosphorylated Ser-139 residue of the histone variant H2AX (γH2AX) foci linked to mechanotransduction pressure in vivo during serial orthotopic passages of mouse KrasLSL-G12D/+;Trp53flox/flox;Pdx1-Cre (KPC) cancer cells in the tumor and in migrating through the size-restricted 3-μm micropores. To search for pancreatic cancer cell-of-origin, analysis of single-cell data sets revealed that the extracellular matrix shaped an alternate route of acinar-ductal transdifferentiation of acinar cells into topoisomerase II α (TOP2A)-overexpressing cancer cells and derived subclusters with copy number amplifications in MYC-PTK2 (protein tyrosine kinase 2) locus and PIK3CA. High-PTK2 expression is associated with 171 differentially methylated CpG loci, 319 differentially expressed genes, and poor overall survival in The Cancer Genome Atlas-Pancreatic Adenocarcinoma cohort. Abolished RGD-integrin signaling by disintegrin KG blocked the PTK2 phosphorylation, increased cancer apoptosis, decreased vav guanine nucleotide exchange factor 1 (VAV1) expression, and prolonged overall survival in the KPC mice. Reduction of α-smooth muscle actin deposition in the CD248 knockout KPC mice remodeled the tissue stroma and down-regulated TOP2A expression in the epithelium. In summary, stromal stiffness induced the onset of cancer cells-of-origin by ectopic TOP2A expression, and the genomic amplification of MYC-PTK2 locus via alternative transdifferentiation of pancreatic progenitor cells is the vulnerability useful for disintegrin KG treatment.

Alcohol-induced epigenetic changes prevent fibrosis resolution after alcohol cessation in miceresolution

BACKGROUND AND AIMS

Alcohol-associated liver disease is a major cause of alcohol-associated mortality. Recently, we identified hepatic demethylases lysine demethylase (KDM)5B and KDM5C as important epigenetic regulators of alcohol response in the liver. In this study, we aimed to investigate the role of KDM5 demethylases in alcohol-associated liver disease resolution.

APPROACH AND RESULTS

We showed that alcohol-induced liver steatosis rapidly resolved after alcohol cessation. In contrast, fibrosis persisted in the liver for up to 8 weeks after the end of alcohol exposure. Defects in fibrosis resolution were in part due to alcohol-induced KDM5B and KDM5C-dependent epigenetic changes in hepatocytes. Using cell-type-specific knockout mice, we found that adeno-associated virus-mediated knockout of KDM5B and KDM5C demethylases in hepatocytes at the time of alcohol withdrawal promoted fibrosis resolution. Single-cell ATAC sequencing analysis showed that during alcohol-associated liver disease resolution epigenetic cell states largely reverted to control conditions. In addition, we found unique epigenetic cell states distinct from both control and alcohol states and identified associated transcriptional regulators, including liver X receptor (LXR) alpha (α). In vitro and in vivo analysis confirmed that knockout of KDM5B and KDM5C demethylases promoted LXRα activity, likely through regulation of oxysterol biosynthesis, and this activity was critical for the fibrosis resolution process. Reduced LXR activity by small molecule inhibitors prevented fibrosis resolution in KDM5-deficient mice.

CONCLUSIONS

In summary, KDM5B and KDM5C demethylases prevent liver fibrosis resolution after alcohol cessation in part through suppression of LXR activity.

Patient-Derived Tumor Organoids Combined with Function-Associated ScRNA-Seq for Dissecting the Local Immune Response of Lung Cancer

In vitro models coupled with multimodal approaches are needed to dissect the dynamic response of local tumor immune microenvironment (TIME) to immunotherapy. Here the patient-derived primary lung cancer organoids (pLCOs) are generated by isolating tumor cell clusters, including the infiltrated immune cells. A function-associated single-cell RNA sequencing (FascRNA-seq) platform allowing both phenotypic evaluation and scRNA-seq at single-organoid level is developed to dissect the TIME of individual pLCOs. The analysis of 171 individual pLCOs derived from seven patients reveals that pLCOs retain the TIME heterogeneity in the parenchyma of parental tumor tissues, providing models with identical genetic background but various TIME. Linking the scRNA-seq data of individual pLCOs with their responses to anti-PD-1 (αPD-1) immune checkpoint blockade (ICB) allows to confirm the central role of CD8+ T cells in anti-tumor immunity, to identify potential tumor-reactive T cells with a set of 10 genes, and to unravel the factors regulating T cell activity, including CD99 gene. In summary, the study constructs a joint phenotypic and transcriptomic FascRNA-seq platform to dissect the dynamic response of local TIME under ICB treatment, providing a promising approach to evaluate novel immunotherapies and to understand the underlying molecular mechanisms.

Comprehensive analysis of single cell and bulk RNA sequencing reveals the heterogeneity of melanoma tumor microenvironment and predicts the response of immunotherapy

BACKGROUND

Tumor microenvironment (TME) heterogeneity is an important factor affecting the treatment response of immune checkpoint inhibitors (ICI). However, the TME heterogeneity of melanoma is still widely characterized.

METHODS

We downloaded the single-cell sequencing data sets of two melanoma patients from the GEO database, and used the "Scissor" algorithm and the "BayesPrism" algorithm to comprehensively analyze the characteristics of microenvironment cells based on single-cell and bulk RNA-seq data. The prediction model of immunotherapy response was constructed by machine learning and verified in three cohorts of GEO database.

RESULTS

We identified seven cell types. In the Scissor+ subtype cell population, the top three were T cells, B cells and melanoma cells. In the Scissor- subtype, there are more macrophages. By quantifying the characteristics of TME, significant differences in B cells between responders and non-responders were observed. The higher the proportion of B cells, the better the prognosis. At the same time, macrophages in the non-responsive group increased significantly. Finally, nine gene features for predicting ICI response were constructed, and their predictive performance was superior in three external validation groups.

CONCLUSION

Our study revealed the heterogeneity of melanoma TME and found a new predictive biomarker, which provided theoretical support and new insights for precise immunotherapy of melanoma patients.

Subpopulations of fibroblasts derived from human iPS cells

Organ fibrosis causes collagen fiber overgrowth and impairs organ function. Cardiac fibrosis after myocardial infarction impairs cardiac function significantly, pulmonary fibrosis reduces gas exchange efficiency, and liver fibrosis disturbs the natural function of the liver. Its development is associated with the differentiation of fibroblasts into myofibroblasts and increased collagen synthesis. Fibrosis has organ specificity, defined by the heterogeneity of fibroblasts. Although this heterogeneity is established during embryonic development, it has not been defined yet. Fibroblastic differentiation of induced pluripotent stem cells (iPSCs) recapitulates the process by which fibroblasts acquire diversity. Here, we differentiated iPSCs into cardiac, hepatic, and dermal fibroblasts and analyzed their properties using single-cell RNA sequencing. We observed characteristic subpopulations with different ratios in each organ-type fibroblast group, which contained both resting and distinct ACTA2+ myofibroblasts. These findings provide crucial information on the ontogeny-based heterogeneity of fibroblasts, leading to the development of therapeutic strategies to control fibrosis.

Intravesical instillation-based mTOR-STAT3 dual targeting for bladder cancer treatment

BACKGROUND

Recent intravesical administration of adenoviral vectors, either as a single injection or in combination with immune checkpoint inhibitors, exemplified by cretostimogene grenadenorepvec and nadofaragene firadenovec, has demonstrated remarkable efficacy in clinical trials for non-muscle invasive bladder cancer. Despite their ability to induce an enhanced immune reaction within the lesion, the intracellular survival signaling of cancer cells has not been thoroughly addressed.

METHODS

An analysis of the prognostic data revealed a high probability of therapeutic efficacy with simultaneous inhibition of mTOR and STAT3. Considering the challenges of limited pharmaco-accessibility to the bladder due to its pathophysiological structure and the partially undruggable nature of target molecules, we designed a dual siRNA system targeting both mRNAs. Subsequently, this dual siRNA system was encoded into the adenovirus 5/3 (Ad 5/3) to enhance in vivo delivery efficiency.

RESULTS

Gene-targeting efficacy was assessed using cells isolated from xenografted tumors using a single-cell analysis system. Our strategy demonstrated a balanced downregulation of mTOR and STAT3 at the single-cell resolution, both in vitro and in vivo. This approach reduced tumor growth in bladder cancer xenograft and orthotopic animal experiments. In addition, increased infiltration of CD8+ T cells was observed in a humanized mouse model. We provided helpful and safe tissue distribution data for intravesical therapy of siRNAs coding adenoviruses.

CONCLUSIONS

The bi-specific siRNA strategy, encapsulated in an adenovirus, could be a promising tool to augment cancer treatment efficacy and overcome conventional therapy limitations associated with "undruggability." Hence, we propose that dual targeting of mTOR and STAT3 is an advantageous strategy for intravesical therapy using adenoviruses.

GATA2 heterozygosity causes an epigenetic feedback mechanism resulting in myeloid and erythroid dysplasia

The transcription factor GATA2 has a pivotal role in haematopoiesis. Heterozygous germline GATA2 mutations result in a syndrome characterized by immunodeficiency, bone marrow failure and predispositions to myelodysplastic syndrome (MDS) and acute myeloid leukaemia. Clinical symptoms in these patients are diverse and mechanisms driving GATA2-related phenotypes are largely unknown. To explore the impact of GATA2 haploinsufficiency on haematopoiesis, we generated a zebrafish model carrying a heterozygous mutation of gata2b (gata2b+/-), an orthologue of GATA2. Morphological analysis revealed myeloid and erythroid dysplasia in gata2b+/- kidney marrow. Because Gata2b could affect both transcription and chromatin accessibility during lineage differentiation, this was assessed by single-cell (sc) RNA-seq and single-nucleus (sn) ATAC-seq. Sn-ATAC-seq showed that the co-accessibility between the transcription start site (TSS) and a -3.5-4.1 kb putative enhancer was more robust in gata2b+/- zebrafish HSPCs compared to wild type, increasing gata2b expression and resulting in higher genome-wide Gata2b motif use in HSPCs. As a result of increased accessibility of the gata2b locus, gata2b+/- chromatin was also more accessible during lineage differentiation. scRNA-seq data revealed myeloid differentiation defects, that is, impaired cell cycle progression, reduced expression of cebpa and cebpb and increased signatures of ribosome biogenesis. These data also revealed a differentiation delay in erythroid progenitors, aberrant proliferative signatures and down-regulation of Gata1a, a master regulator of erythropoiesis, which worsened with age. These findings suggest that cell-intrinsic compensatory mechanisms, needed to obtain normal levels of Gata2b in heterozygous HSPCs to maintain their integrity, result in aberrant lineage differentiation, thereby representing a critical step in the predisposition to MDS.

Protocol for multimodal profiling of human kidneys with simultaneous high-throughput ATAC and RNA expression with sequencing

Simultaneous high-throughput ATAC and RNA expression with sequencing (SHARE-seq) profiles transcriptomics and chromatin accessibility in the same cells at high throughput. Here, we present a protocol for multimodal profiling of human kidneys with SHARE-seq. We describe steps for processing fixed nuclei for SHARE-seq split-pool barcoding and library preparation. We also detail how to determine the optimal working concentration of Tn5 transposase for transposition and tagmentation. This protocol allows researchers to generate large-scale single-cell multiomics data at low reagent cost. For complete details on the use and execution of this protocol, please refer to Li et al.1.

Traumatic brain injury promotes neurogenesis at the cost of astrogliogenesis in the adult hippocampus of male mice

Traumatic brain injury (TBI) can result in long-lasting changes in hippocampal function. The changes induced by TBI on the hippocampus contribute to cognitive deficits. The adult hippocampus harbors neural stem cells (NSCs) that generate neurons (neurogenesis), and astrocytes (astrogliogenesis). While deregulation of hippocampal NSCs and neurogenesis have been observed after TBI, it is not known how TBI may affect hippocampal astrogliogenesis. Using a controlled cortical impact model of TBI in male mice, single cell RNA sequencing and spatial transcriptomics, we assessed how TBI affected hippocampal NSCs and the neuronal and astroglial lineages derived from them. We observe an increase in NSC-derived neuronal cells and a concomitant decrease in NSC-derived astrocytic cells, together with changes in gene expression and cell dysplasia within the dentate gyrus. Here, we show that TBI modifies NSC fate to promote neurogenesis at the cost of astrogliogenesis and identify specific cell populations as possible targets to counteract TBI-induced cellular changes in the adult hippocampus.

A single-cell atlas of pig gastrulation as a resource for comparative embryology

Cell-fate decisions during mammalian gastrulation are poorly understood outside of rodent embryos. The embryonic disc of pig embryos mirrors humans, making them a useful proxy for studying gastrulation. Here we present a single-cell transcriptomic atlas of pig gastrulation, revealing cell-fate emergence dynamics, as well as conserved and divergent gene programs governing early porcine, primate, and murine development. We highlight heterochronicity in extraembryonic cell-types, despite the broad conservation of cell-type-specific transcriptional programs. We apply these findings in combination with functional investigations, to outline conserved spatial, molecular, and temporal events during definitive endoderm specification. We find early FOXA2 + /TBXT- embryonic disc cells directly form definitive endoderm, contrasting later-emerging FOXA2/TBXT+ node/notochord progenitors. Unlike mesoderm, none of these progenitors undergo epithelial-to-mesenchymal transition. Endoderm/Node fate hinges on balanced WNT and hypoblast-derived NODAL, which is extinguished upon endodermal differentiation. These findings emphasise the interplay between temporal and topological signalling in fate determination during gastrulation.

A Cell Cycle-Aware Network for Data Integration and Label Transferring of Single-Cell RNA-Seq and ATAC-Seq

In recent years, the integration of single-cell multi-omics data has provided a more comprehensive understanding of cell functions and internal regulatory mechanisms from a non-single omics perspective, but it still suffers many challenges, such as omics-variance, sparsity, cell heterogeneity, and confounding factors. As it is known, the cell cycle is regarded as a confounder when analyzing other factors in single-cell RNA-seq data, but it is not clear how it will work on the integrated single-cell multi-omics data. Here, a cell cycle-aware network (CCAN) is developed to remove cell cycle effects from the integrated single-cell multi-omics data while keeping the cell type-specific variations. This is the first computational model to study the cell-cycle effects in the integration of single-cell multi-omics data. Validations on several benchmark datasets show the outstanding performance of CCAN in a variety of downstream analyses and applications, including removing cell cycle effects and batch effects of scRNA-seq datasets from different protocols, integrating paired and unpaired scRNA-seq and scATAC-seq data, accurately transferring cell type labels from scRNA-seq to scATAC-seq data, and characterizing the differentiation process from hematopoietic stem cells to different lineages in the integration of differentiation data.

Single-cell signatures identify microenvironment factors in tumors associated with patient outcomes

The cellular components of tumors and their microenvironment play pivotal roles in tumor progression, patient survival, and the response to cancer treatments. Unveiling a comprehensive cellular profile within bulk tumors via single-cell RNA sequencing (scRNA-seq) data is crucial, as it unveils intrinsic tumor cellular traits that elude identification through conventional cancer subtyping methods. Our contribution, scBeacon, is a tool that derives cell-type signatures by integrating and clustering multiple scRNA-seq datasets to extract signatures for deconvolving unrelated tumor datasets on bulk samples. Through the employment of scBeacon on the The Cancer Genome Atlas (TCGA) cohort, we find cellular and molecular attributes within specific tumor categories, many with patient outcome relevance. We developed a tumor cell-type map to visually depict the relationships among TCGA samples based on the cell-type inferences.

Reconstructing axial progenitor field dynamics in mouse stem cell-derived embryoids

Embryogenesis requires substantial coordination to translate genetic programs to the collective behavior of differentiating cells, but understanding how cellular decisions control tissue morphology remains conceptually and technically challenging. Here, we combine continuous Cas9-based molecular recording with a mouse embryonic stem cell-based model of the embryonic trunk to build single-cell phylogenies that describe the behavior of transient, multipotent neuro-mesodermal progenitors (NMPs) as they commit into neural and somitic cell types. We find that NMPs show subtle transcriptional signatures related to their recent differentiation and contribute to downstream lineages through a surprisingly broad distribution of individual fate outcomes. Although decision-making can be heavily influenced by environmental cues to induce morphological phenotypes, axial progenitors intrinsically mature over developmental time to favor the neural lineage. Using these data, we present an experimental and analytical framework for exploring the non-homeostatic dynamics of transient progenitor populations as they shape complex tissues during critical developmental windows.

Under pressure: integrated endothelial cell response to hydrostatic and shear stresses

Blood flow within the vasculature is a critical determinant of endothelial cell (EC) identity and functionality, yet the intricate interplay of various hemodynamic forces and their collective impact on endothelial and vascular responses are not fully understood. Specifically, the role of hydrostatic pressure in the EC flow response is understudied, despite its known significance in vascular development and disease. To address this gap, we developed in vitro models to investigate how pressure influences EC responses to flow. Our study demonstrates that elevated pressure conditions significantly modify shear-induced flow alignment and increase endothelial cell density. Bulk and single-cell RNA sequencing analyses revealed that, while shear stress remains the primary driver of flow-induced transcriptional changes, pressure modulates shear-induced signaling in a dose-dependent manner. These pressure-responsive transcriptional signatures identified in human ECs were conserved during the onset of circulation in early mouse embryonic vascular development, where pressure was notably associated with transcriptional programs essential to arterial and hemogenic EC fates. Our findings suggest that pressure plays a synergistic role with shear stress on ECs and emphasizes the need for an integrative approach to endothelial cell mechanotransduction, one that encompasses the effects induced by pressure alongside other hemodynamic forces.

Genetic Implication of Prenatal GABAergic and Cholinergic Neuron Development in Susceptibility to Schizophrenia

BACKGROUND

The ganglionic eminences (GE) are fetal-specific structures that give rise to gamma-aminobutyric acid (GABA)- and acetylcholine-releasing neurons of the forebrain. Given the evidence for GABAergic, cholinergic, and neurodevelopmental disturbances in schizophrenia, we tested the potential involvement of GE neuron development in mediating genetic risk for the condition.

STUDY DESIGN

We combined data from a recent large-scale genome-wide association study of schizophrenia with single-cell RNA sequencing data from the human GE to test the enrichment of schizophrenia risk variation in genes with high expression specificity for developing GE cell populations. We additionally performed the single nuclei Assay for Transposase-Accessible Chromatin with Sequencing (snATAC-Seq) to map potential regulatory genomic regions operating in individual cell populations of the human GE, using these to test for enrichment of schizophrenia common genetic variant liability and to functionally annotate non-coding variants-associated with the disorder.

STUDY RESULTS

Schizophrenia common variant liability was enriched in genes with high expression specificity for developing neuron populations that are predicted to form dopamine D1 and D2 receptor-expressing GABAergic medium spiny neurons of the striatum, cortical somatostatin-positive GABAergic interneurons, calretinin-positive GABAergic neurons, and cholinergic neurons. Consistent with these findings, schizophrenia genetic risk was concentrated in predicted regulatory genomic sequence mapped in developing neuronal populations of the GE.

CONCLUSIONS

Our study implicates prenatal development of specific populations of GABAergic and cholinergic neurons in later susceptibility to schizophrenia, and provides a map of predicted regulatory genomic elements operating in cells of the GE.

Decoding temporal heterogeneity in NSCLC through machine learning and prognostic model construction

BACKGROUND

Non-small cell lung cancer (NSCLC) is a prevalent and heterogeneous disease with significant genomic variations between the early and advanced stages. The identification of key genes and pathways driving NSCLC tumor progression is critical for improving the diagnosis and treatment outcomes of this disease.

METHODS

In this study, we conducted single-cell transcriptome analysis on 93,406 cells from 22 NSCLC patients to characterize malignant NSCLC cancer cells. Utilizing cNMF, we classified these cells into distinct modules, thus identifying the diverse molecular profiles within NSCLC. Through pseudotime analysis, we delineated temporal gene expression changes during NSCLC evolution, thus demonstrating genes associated with disease progression. Using the XGBoost model, we assessed the significance of these genes in the pseudotime trajectory. Our findings were validated by using transcriptome sequencing data from The Cancer Genome Atlas (TCGA), supplemented via LASSO regression to refine the selection of characteristic genes. Subsequently, we established a risk score model based on these genes, thus providing a potential tool for cancer risk assessment and personalized treatment strategies.

RESULTS

We used cNMF to classify malignant NSCLC cells into three functional modules, including the metabolic reprogramming module, cell cycle module, and cell stemness module, which can be used for the functional classification of malignant tumor cells in NSCLC. These findings also indicate that metabolism, the cell cycle, and tumor stemness play important driving roles in the malignant evolution of NSCLC. We integrated cNMF and XGBoost to select marker genes that are indicative of both early and advanced NSCLC stages. The expression of genes such as CHCHD2, GAPDH, and CD24 was strongly correlated with the malignant evolution of NSCLC at the single-cell data level. These genes have been validated via histological data. The risk score model that we established (represented by eight genes) was ultimately validated with GEO data.

CONCLUSION

In summary, our study contributes to the identification of temporal heterogeneous biomarkers in NSCLC, thus offering insights into disease progression mechanisms and potential therapeutic targets. The developed workflow demonstrates promise for future applications in clinical practice.

CellRank 2: unified fate mapping in multiview single-cell data

Single-cell RNA sequencing allows us to model cellular state dynamics and fate decisions using expression similarity or RNA velocity to reconstruct state-change trajectories; however, trajectory inference does not incorporate valuable time point information or utilize additional modalities, whereas methods that address these different data views cannot be combined or do not scale. Here we present CellRank 2, a versatile and scalable framework to study cellular fate using multiview single-cell data of up to millions of cells in a unified fashion. CellRank 2 consistently recovers terminal states and fate probabilities across data modalities in human hematopoiesis and endodermal development. Our framework also allows combining transitions within and across experimental time points, a feature we use to recover genes promoting medullary thymic epithelial cell formation during pharyngeal endoderm development. Moreover, we enable estimating cell-specific transcription and degradation rates from metabolic-labeling data, which we apply to an intestinal organoid system to delineate differentiation trajectories and pinpoint regulatory strategies.

pTDP-43 levels correlate with cell type specific molecular alterations in the prefrontal cortex of C9orf72 ALS/FTD patients

Repeat expansions in the C9orf72 gene are the most common genetic cause of amyotrophic lateral sclerosis and familial frontotemporal dementia (ALS/FTD). To identify molecular defects that take place in the dorsolateral frontal cortex of patients with C9orf72 ALS/FTD, we compared healthy controls with C9orf72 ALS/FTD donor samples staged based on the levels of cortical phosphorylated TAR DNA binding protein (pTDP-43), a neuropathological hallmark of disease progression. We identified distinct molecular changes in different cell types that take place during FTD development. Loss of neurosurveillance microglia and activation of the complement cascade take place early, when pTDP-43 aggregates are absent or very low, and become more pronounced in late stages, suggesting an initial involvement of microglia in disease progression. Reduction of layer 2-3 cortical projection neurons with high expression of CUX2/LAMP5 also occurs early, and the reduction becomes more pronounced as pTDP-43 accumulates. Several unique features were observed only in samples with high levels of pTDP-43, including global alteration of chromatin accessibility in oligodendrocytes, microglia, and astrocytes; higher ratios of premature oligodendrocytes; increased levels of the noncoding RNA NEAT1 in astrocytes and neurons, and higher amount of phosphorylated ribosomal protein S6. Our findings reveal previously unknown progressive functional changes in major cell types found in the frontal cortex of C9orf72 ALS/FTD patients that shed light on the mechanisms underlying the pathology of this disease.

SIMS: A deep-learning label transfer tool for single-cell RNA sequencing analysis

Cell atlases serve as vital references for automating cell labeling in new samples, yet existing classification algorithms struggle with accuracy. Here we introduce SIMS (scalable, interpretable machine learning for single cell), a low-code data-efficient pipeline for single-cell RNA classification. We benchmark SIMS against datasets from different tissues and species. We demonstrate SIMS's efficacy in classifying cells in the brain, achieving high accuracy even with small training sets (<3,500 cells) and across different samples. SIMS accurately predicts neuronal subtypes in the developing brain, shedding light on genetic changes during neuronal differentiation and postmitotic fate refinement. Finally, we apply SIMS to single-cell RNA datasets of cortical organoids to predict cell identities and uncover genetic variations between cell lines. SIMS identifies cell-line differences and misannotated cell lineages in human cortical organoids derived from different pluripotent stem cell lines. Altogether, we show that SIMS is a versatile and robust tool for cell-type classification from single-cell datasets.

Disease activity drives divergent epigenetic and transcriptomic reprogramming of monocyte subpopulations in systemic lupus erythematosus

OBJECTIVES

Systemic lupus erythematosus (SLE) is characterised by systemic inflammation involving various immune cell types. Monocytes, pivotal in promoting and regulating inflammation in SLE, differentiate from classic monocytes into intermediate and non-classic monocytes, assuming diverse roles and changing their proportions in inflammation. In this study, we investigated the epigenetic and transcriptomic profiles of these and novel monocyte subsets in SLE in relation to activity and progression.

METHODS

We obtained the DNA methylomes and transcriptomes of classic, intermediate, non-classic monocytes in patients with SLE (at first and follow-up visits) and healthy donors. We integrated these data with single-cell transcriptomics of SLE and healthy donors and interrogated their relationships with activity and progression.

RESULTS

In addition to shared DNA methylation and transcriptomic alterations associated with a strong interferon signature, we identified monocyte subset-specific alterations, especially in DNA methylation, which reflect an impact of SLE on monocyte differentiation. SLE classic monocytes exhibited a proinflammatory profile and were primed for macrophage differentiation. SLE non-classic monocytes displayed a T cell differentiation-related phenotype, with Th17-regulating features. Changes in monocyte proportions, DNA methylation and expression occurred in relation to disease activity and involved the STAT pathway. Integration of bulk with single-cell RNA sequencing datasets revealed disease activity-dependent expansion of SLE-specific monocyte subsets, further supported the interferon signature for classic monocytes, and associated intermediate and non-classic populations with exacerbated complement activation.

CONCLUSIONS

Disease activity in SLE drives a subversion of the epigenome and transcriptome programme in monocyte differentiation, impacting the function of different subsets and allowing to generate predictive methods for activity and progression.

Spatiotemporal cell landscape of human embryonic tooth development

Understanding the cellular composition and trajectory of human tooth development is valuable for dentistry and stem cell engineering research. Previous single-cell studies have focused on mature human teeth and developing mouse teeth, but the cell landscape of human embryonic dental development is still unknown. In this study, tooth germ tissues were collected from aborted foetus (17-24 weeks) for single-cell RNA sequence and spatial transcriptome analysis. The cells were classified into seven subclusters of epithelium, and seven clusters of mesenchyme, as well as other cell types such as Schwann cell precursor and pericyte. For epithelium, the stratum intermedium branch and the ameloblast branch diverged from the same set of outer enamel-inner enamel-ALCAM+ epithelial cell lineage, but their spatial distribution of two branches was not clearly distinct. This trajectory received spatially adjacent regulation signals from mesenchyme and pericyte, including JAG1 and APP. The differentiation of pulp cell and pre-odontoblast showed four waves of temporally distinct gene expression, which involved regulation networks of LHX9, DLX5 and SP7, and these genes were regulated by upstream ligands such as the BMP family. This provides a reference landscape for the research on early human tooth development, covering different spatial structures and developmental periods.

An inducible genetic tool to track and manipulate specific microglial states reveals their plasticity and roles in remyelination

Recent single-cell RNA sequencing studies have revealed distinct microglial states in development and disease. These include proliferative-region-associated microglia (PAMs) in developing white matter and disease-associated microglia (DAMs) prevalent in various neurodegenerative conditions. PAMs and DAMs share a similar core gene signature. However, the extent of the dynamism and plasticity of these microglial states, as well as their functional significance, remains elusive, partly due to the lack of specific tools. Here, we generated an inducible Cre driver line, Clec7a-CreERT2, that targets PAMs and DAMs in the brain parenchyma. Utilizing this tool, we profiled labeled cells during development and in several disease models, uncovering convergence and context-dependent differences in PAM and DAM gene expression. Through long-term tracking, we demonstrated microglial state plasticity. Lastly, we specifically depleted DAMs in demyelination, revealing their roles in disease recovery. Together, we provide a versatile genetic tool to characterize microglial states in CNS development and disease.

DOT: a flexible multi-objective optimization framework for transferring features across single-cell and spatial omics

Single-cell transcriptomics and spatially-resolved imaging/sequencing technologies have revolutionized biomedical research. However, they suffer from lack of spatial information and a trade-off of resolution and gene coverage, respectively. We propose DOT, a multi-objective optimization framework for transferring cellular features across these data modalities, thus integrating their complementary information. DOT uses genes beyond those common to the data modalities, exploits the local spatial context, transfers spatial features beyond cell-type information, and infers absolute/relative abundance of cell populations at tissue locations. Thus, DOT bridges single-cell transcriptomics data with both high- and low-resolution spatially-resolved data. Moreover, DOT combines practical aspects related to cell composition, heterogeneity, technical effects, and integration of prior knowledge. Our fast implementation based on the Frank-Wolfe algorithm achieves state-of-the-art or improved performance in localizing cell features in high- and low-resolution spatial data and estimating the expression of unmeasured genes in low-coverage spatial data.

Progenitor-like cells contributing to cellular heterogeneity in the nucleus pulposus are lost in intervertebral disc degeneration

The nucleus pulposus (NP) in the intervertebral disc (IVD) arises from embryonic notochord. Loss of notochordal-like cells in humans correlates with onset of IVD degeneration, suggesting that they are critical for healthy NP homeostasis and function. Comparative transcriptomic analyses identified expression of progenitor-associated genes (GREM1, KRT18, and TAGLN) in the young mouse and non-degenerated human NP, with TAGLN expression reducing with aging. Lineage tracing using Tagln-CreERt2 mice identified peripherally located proliferative NP (PeriNP) cells in developing and postnatal NP that provide a continuous supply of cells to the entire NP. PeriNP cells were diminished in aged mice and absent in puncture-induced degenerated discs. Single-cell transcriptomes of postnatal Tagln-CreERt2 IVD cells indicate enrichment for TGF-β signaling in Tagln descendant NP sub-populations. Notochord-specific removal of TGF-β/BMP mediator Smad4 results in loss of Tagln+ cells and abnormal NP morphologies. We propose Tagln+ PeriNP cells are potential progenitors crucial for NP homeostasis.

Metric multidimensional scaling for large single-cell datasets using neural networks

Metric multidimensional scaling is one of the classical methods for embedding data into low-dimensional Euclidean space. It creates the low-dimensional embedding by approximately preserving the pairwise distances between the input points. However, current state-of-the-art approaches only scale to a few thousand data points. For larger data sets such as those occurring in single-cell RNA sequencing experiments, the running time becomes prohibitively large and thus alternative methods such as PCA are widely used instead. Here, we propose a simple neural network-based approach for solving the metric multidimensional scaling problem that is orders of magnitude faster than previous state-of-the-art approaches, and hence scales to data sets with up to a few million cells. At the same time, it provides a non-linear mapping between high- and low-dimensional space that can place previously unseen cells in the same embedding.

From whole bodies to single cells: A guide to transcriptomic approaches for ecology and evolutionary biology

RNA sequencing (RNAseq) methodology has experienced a burst of technological developments in the last decade, which has opened up opportunities for studying the mechanisms of adaptation to environmental factors at both the organismal and cellular level. Selecting the most suitable experimental approach for specific research questions and model systems can, however, be a challenge and researchers in ecology and evolution are commonly faced with the choice of whether to study gene expression variation in whole bodies, specific tissues, and/or single cells. A wide range of sometimes polarised opinions exists over which approach is best. Here, we highlight the advantages and disadvantages of each of these approaches to provide a guide to help researchers make informed decisions and maximise the power of their study. Using illustrative examples of various ecological and evolutionary research questions, we guide the readers through the different RNAseq approaches and help them identify the most suitable design for their own projects.

Decoding changes in tumor-infiltrating leukocytes through dynamic experimental models and single-cell technologies

The ability to characterize immune cells and explore the molecular interactions that govern their functions has never been greater, fueled in recent years by the revolutionary advance of single-cell analysis platforms. However, precisely how immune cells respond to different stimuli and where differentiation processes and effector functions operate remain incompletely understood. Inferring cellular fate within single-cell transcriptomic analyses is now omnipresent, despite the assumptions typically required in such analyses. Recently developed experimental models support dynamic analyses of the immune response, providing insights into the temporal changes that occur within cells and the tissues in which such transitions occur. Here we will review these approaches and discuss how these can be combined with single-cell technologies to develop a deeper understanding of the immune responses that should support the development of better therapeutic options for patients.

CAbiNet: joint clustering and visualization of cells and genes for single-cell transcriptomics

A fundamental analysis task for single-cell transcriptomics data is clustering with subsequent visualization of cell clusters. The genes responsible for the clustering are only inferred in a subsequent step. Clustering cells and genes together would be the remit of biclustering algorithms, which are often bogged down by the size of single-cell data. Here we present 'Correspondence Analysis based Biclustering on Networks' (CAbiNet) for joint clustering and visualization of single-cell RNA-sequencing data. CAbiNet performs efficient co-clustering of cells and their respective marker genes and jointly visualizes the biclusters in a non-linear embedding for easy and interactive visual exploration of the data.

Gene regulatory network topology governs resistance and treatment escape in glioma stem-like cells

Poor prognosis and drug resistance in glioblastoma (GBM) can result from cellular heterogeneity and treatment-induced shifts in phenotypic states of tumor cells, including dedifferentiation into glioma stem-like cells (GSCs). This rare tumorigenic cell subpopulation resists temozolomide, undergoes proneural-to-mesenchymal transition (PMT) to evade therapy, and drives recurrence. Through inference of transcriptional regulatory networks (TRNs) of patient-derived GSCs (PD-GSCs) at single-cell resolution, we demonstrate how the topology of transcription factor interaction networks drives distinct trajectories of cell-state transitions in PD-GSCs resistant or susceptible to cytotoxic drug treatment. By experimentally testing predictions based on TRN simulations, we show that drug treatment drives surviving PD-GSCs along a trajectory of intermediate states, exposing vulnerability to potentiated killing by siRNA or a second drug targeting treatment-induced transcriptional programs governing nongenetic cell plasticity. Our findings demonstrate an approach to uncover TRN topology and use it to rationally predict combinatorial treatments that disrupt acquired resistance in GBM.

Modulation of FGF pathway signaling and vascular differentiation using designed oligomeric assemblies

Many growth factors and cytokines signal by binding to the extracellular domains of their receptors and driving association and transphosphorylation of the receptor intracellular tyrosine kinase domains, initiating downstream signaling cascades. To enable systematic exploration of how receptor valency and geometry affect signaling outcomes, we designed cyclic homo-oligomers with up to 8 subunits using repeat protein building blocks that can be modularly extended. By incorporating a de novo-designed fibroblast growth factor receptor (FGFR)-binding module into these scaffolds, we generated a series of synthetic signaling ligands that exhibit potent valency- and geometry-dependent Ca2+ release and mitogen-activated protein kinase (MAPK) pathway activation. The high specificity of the designed agonists reveals distinct roles for two FGFR splice variants in driving arterial endothelium and perivascular cell fates during early vascular development. Our designed modular assemblies should be broadly useful for unraveling the complexities of signaling in key developmental transitions and for developing future therapeutic applications.

Dissecting and Visualizing the Functional Diversity of Cardiac Macrophages

Cardiac macrophages represent a functionally diverse population of cells involved in cardiac homeostasis, repair, and remodeling. With recent advancements in single-cell technologies, it is possible to elucidate specific macrophage subsets based on transcriptional signatures and cell surface protein expression to gain a deep understanding of macrophage diversity in the heart. The use of fate-mapping technologies and parabiosis studies have provided insight into the ontogeny and dynamics of macrophages identifying subsets derived from embryonic and adult definitive hematopoietic progenitors that include tissue-resident and bone marrow monocyte-derived macrophages, respectively. Within the heart, these subsets have distinct tissue niches and functional roles in the setting of homeostasis and disease, with cardiac resident macrophages representing a protective cell population while bone marrow monocyte-derived cardiac macrophages have a context-dependent effect, triggering both proinflammatory tissue injury, but also promoting reparative functions. With the increased understanding of the clinical relevance of cardiac macrophage subsets, there has been an increasing need to detect and measure cardiac macrophage compositions in living animals and patients. New molecular tracers compatible with positron emission tomography/computerized tomography and positron emission tomography/ magnetic resonance imaging have enabled investigators to noninvasively and serially visualize cardiac macrophage subsets within the heart to define associations with disease and measure treatment responses. Today, advancements within this thriving field are poised to fuel an era of clinical translation.

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.

Heterochromatin protein ERH represses alternative cell fates during early mammalian differentiation

During development, H3K9me3 heterochromatin is dynamically rearranged, silencing repeat elements and protein coding genes to restrict cell identity. Enhancer of Rudimentary Homolog (ERH) is an evolutionarily conserved protein originally characterized in fission yeast and recently shown to be required for H3K9me3 maintenance in human fibroblasts, but its function during development remains unknown. Here, we show that ERH is required for proper segregation of the inner cell mass and trophectoderm cell lineages during mouse development by repressing totipotent and alternative lineage programs. During human embryonic stem cell (hESC) differentiation into germ layer lineages, ERH is crucial for silencing naïve and pluripotency genes, transposable elements, and alternative lineage genes. Strikingly, ERH depletion in somatic cells reverts the H3K9me3 landscape to an hESC state and enables naïve and pluripotency gene and transposable element activation during iPSC reprogramming. Our findings reveal a role for ERH in initiation and maintenance of developmentally established gene repression.

Spatiotemporal transcriptomic landscape of rice embryonic cells during seed germination

Characterizing cellular features during seed germination is crucial for understanding the complex biological functions of different embryonic cells in regulating seed vigor and seedling establishment. We performed spatially enhanced resolution omics sequencing (Stereo-seq) and single-cell RNA sequencing (scRNA-seq) to capture spatially resolved single-cell transcriptomes of germinating rice embryos. An automated cell-segmentation model, employing deep learning, was developed to accommodate the analysis requirements. The spatial transcriptomes of 6, 24, 36, and 48 h after imbibition unveiled both known and previously unreported embryo cell types, including two unreported scutellum cell types, corroborated by in situ hybridization and functional exploration of marker genes. Temporal transcriptomic profiling delineated gene expression dynamics in distinct embryonic cell types during seed germination, highlighting key genes involved in nutrient metabolism, biosynthesis, and signaling of phytohormones, reprogrammed in a cell-type-specific manner. Our study provides a detailed spatiotemporal transcriptome of rice embryo and presents a previously undescribed methodology for exploring the roles of different embryonic cells in seed germination.

Single-cell transcriptomic analyses of mouse idh1 mutant growth plate chondrocytes reveal distinct cell populations responsible for longitudinal growth and enchondroma formation

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

A spatial portrait of the human sebaceous gland transcriptional program

Sebaceous glands (SG) and their oily secretion (sebum) are indispensable for maintaining skin structure and function, and their deregulation causes skin disorders including but not limited to acne. Recent studies also indicate that sebum may have important immunomodulatory activities and may influence whole-body energy metabolism. However, the progressive transcriptional changes of sebocytes that lead to sebum production have never been characterized in detail. Here, we exploited the high cellular resolution provided by sebaceous hyperplasia and integrated spatial transcriptomics, pseudotime analysis, RNA velocity, and functional enrichment to map the landscape of sebaceous differentiation by employing spatial transcriptomics. Our results were validated by comparison with published SG transcriptome data and further corroborated by assessing the protein expression pattern of a subset of the transcripts in the public repository Human Protein Atlas. Departing from four sebocyte differentiation stages generated by unsupervised clustering, we demonstrate consecutive modulation of cellular functions associable with specific gene sets, from cell proliferation and oxidative phosphorylation via lipid synthesis to cell death. Both validation methods confirmed the biological significance of our results. Our report is complemented by a freely available and browsable online tool. Our data provide the first high-resolution spatial portrait of the SG transcriptional landscape and deliver starting points for experimentally assessing novel candidate molecules for regulating SG homeostasis in health and disease.

Molecular and transcriptional basis of bidirectional CD4+ T cell exhaustion in oropharyngeal squamous cell carcinoma

Tumor-infiltrating CD4+ T cells orchestrate the adaptive immune response through remarkable plasticity, and the expression patterns of exhaustion-related inhibitory receptors in these cells differ significantly from those of CD8+ T cells. Thus, a better understanding of the molecular basis of CD4+ T cell exhaustion and their responses to immune checkpoint blockade (ICB) is required. Here, we integrated multiomics approaches to define the phenotypic and molecular profiles of exhausted CD4+ T cells in oropharyngeal squamous cell carcinoma (OPSCC). Two distinct immune-promoting (Module 1) and immunosuppressive (Module 2) functional modules in tumor-infiltrating CD4+ T cells were identified, and both the immune-promoting function of Module 1 cells and immunosuppressive function of Module 2 cells were positively associated with their corresponding exhaustion states. Furthermore, the application of ICBs targeting effector CD4+ T cells in Module 1 (αPD-1) and Treg cells in Module 2 (αCTLA-4) in mouse models could help reinvigorate the effector function of Module 1-exhausted CD4+ T cells and reduce the immunosuppressive function of Module 2-exhausted CD4+ T cells, ultimately promoting OPSCC tumor regression. Taken together, our study provides a crucial cellular basis for the selection of optimal ICB in treating OPSCC.

Transcriptional control of the Cryptosporidium life cycle

The parasite Cryptosporidium is a leading agent of diarrhoeal disease in young children, and a cause and consequence of chronic malnutrition1,2. There are no vaccines and only limited treatment options3. The parasite infects enterocytes, in which it engages in asexual and sexual replication4, both of which are essential to continued infection and transmission. However, their molecular mechanisms remain largely unclear5. Here we use single-cell RNA sequencing to reveal the gene expression programme of the entire Cryptosporidium parvum life cycle in culture and in infected animals. Diverging from the prevailing model6, we find support for only three intracellular stages: asexual type-I meronts, male gamonts and female gametes. We reveal a highly organized program for the assembly of components at each stage. Dissecting the underlying regulatory network, we identify the transcription factor Myb-M as the earliest determinant of male fate, in an organism that lacks genetic sex determination. Conditional expression of this factor overrides the developmental program and induces widespread maleness, while conditional deletion ablates male development. Both have a profound impact on the infection. A large set of stage-specific genes now provides the opportunity to understand, engineer and disrupt parasite sex and life cycle progression to advance the development of vaccines and treatments.

Chronic Opioid Treatment Arrests Neurodevelopment and Alters Synaptic Activity in Human Midbrain Organoids

Understanding the impact of long-term opioid exposure on the embryonic brain is critical due to the surging number of pregnant mothers with opioid dependency. However, this has been limited by human brain inaccessibility and cross-species differences in animal models. Here, a human midbrain model is established that uses hiPSC-derived midbrain organoids to assess cell-type-specific responses to acute and chronic fentanyl treatment and fentanyl withdrawal. Single-cell mRNA sequencing of 25,510 cells from organoids in different treatment groups reveals that chronic fentanyl treatment arrests neuronal subtype specification during early midbrain development and alters synaptic activity and neuron projection. In contrast, acute fentanyl treatment increases dopamine release but does not significantly alter gene expression related to cell lineage development. These results provide the first examination of the effects of opioid exposure on human midbrain development at the single-cell level.

The effect of sevoflurane exposure on cell-type-specific changes in the prefrontal cortex in young mice

Sevoflurane, the predominant pediatric anesthetic, has been linked to neurotoxicity in young mice, although the underlying mechanisms remain unclear. This study focuses on investigating the impact of neonatal sevoflurane exposure on cell-type-specific alterations in the prefrontal cortex (PFC) of young mice. Neonatal mice were subjected to either control treatment (60% oxygen balanced with nitrogen) or sevoflurane anesthesia (3% sevoflurane in 60% oxygen balanced with nitrogen) for 2 hours on postnatal days (PNDs) 6, 8, and 10. Behavioral tests and single-nucleus RNA sequencing (snRNA-seq) of the PFC were conducted from PNDs 31 to 37. Mechanistic exploration included clustering analysis, identification of differentially expressed genes (DEGs), enrichment analyses, single-cell trajectory analysis, and genome-wide association studies (GWAS). Sevoflurane anesthesia resulted in sociability and cognition impairments in mice. Novel specific marker genes identified 8 distinct cell types in the PFC. Most DEGs between the control and sevoflurane groups were unique to specific cell types. Re-defining 15 glutamatergic neuron subclusters based on layer identity revealed their altered expression profiles. Notably, sevoflurane disrupted the trajectory from oligodendrocyte precursor cells (OPCs) to oligodendrocytes (OLs). Validation of disease-relevant candidate genes across the main cell types demonstrated their association with social dysfunction and working memory impairment. Behavioral results and snRNA-seq collectively elucidated the cellular atlas in the PFC of young male mice, providing a foundation for further mechanistic studies on developmental neurotoxicity induced by anesthesia.

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.

Spatiotemporal transcriptomic atlas of rhizome formation in Oryza longistaminata

Rhizomes are modified stems that grow underground and produce new individuals genetically identical to the mother plant. Recently, a breakthrough has been made in efforts to convert annual grains into perennial ones by utilizing wild rhizomatous species as donors, yet the developmental biology of this organ is rarely studied. Oryza longistaminata, a wild rice species featuring strong rhizomes, provides a valuable model for exploration of rhizome development. Here, we first assembled a double-haplotype genome of O. longistaminata, which displays a 48-fold improvement in contiguity compared to the previously published assembly. Furthermore, spatiotemporal transcriptomics was performed to obtain the expression profiles of different tissues in O. longistaminata rhizomes and tillers. Two spatially reciprocal cell clusters, the vascular bundle 2 cluster and the parenchyma 2 cluster, were determined to be the primary distinctions between the rhizomes and tillers. We also captured meristem initiation cells in the sunken area of parenchyma located at the base of internodes, which is the starting point for rhizome initiation. Trajectory analysis further indicated that the rhizome is regenerated through de novo generation. Collectively, these analyses revealed a spatiotemporal transcriptional transition underlying the rhizome initiation, providing a valuable resource for future perennial crop breeding.

Targeting pro-inflammatory T cells as a novel therapeutic approach to potentially resolve atherosclerosis in humans

Atherosclerosis (AS), a leading cause of cardio-cerebrovascular disease worldwide, is driven by the accumulation of lipid contents and chronic inflammation. Traditional strategies primarily focus on lipid reduction to control AS progression, leaving residual inflammatory risks for major adverse cardiovascular events (MACEs). While anti-inflammatory therapies targeting innate immunity have reduced MACEs, many patients continue to face significant risks. Another key component in AS progression is adaptive immunity, but its potential role in preventing AS remains unclear. To investigate this, we conducted a retrospective cohort study on tumor patients with AS plaques. We found that anti-programmed cell death protein 1 (PD-1) monoclonal antibody (mAb) significantly reduces AS plaque size. With multi-omics single-cell analyses, we comprehensively characterized AS plaque-specific PD-1+ T cells, which are activated and pro-inflammatory. We demonstrated that anti-PD-1 mAb, when captured by myeloid-expressed Fc gamma receptors (FcγRs), interacts with PD-1 expressed on T cells. This interaction turns the anti-PD-1 mAb into a substitute PD-1 ligand, suppressing T-cell functions in the PD-1 ligands-deficient context of AS plaques. Further, we conducted a prospective cohort study on tumor patients treated with anti-PD-1 mAb with or without Fc-binding capability. Our analysis shows that anti-PD-1 mAb with Fc-binding capability effectively reduces AS plaque size, while anti-PD-1 mAb without Fc-binding capability does not. Our work suggests that T cell-targeting immunotherapy can be an effective strategy to resolve AS in humans.

Pro-CRISPR PcrIIC1-associated Cas9 system for enhanced bacterial immunity

The CRISPR system is an adaptive immune system found in prokaryotes that defends host cells against the invasion of foreign DNA1. As part of the ongoing struggle between phages and the bacterial immune system, the CRISPR system has evolved into various types, each with distinct functionalities2. Type II Cas9 is the most extensively studied of these systems and has diverse subtypes. It remains uncertain whether members of this family can evolve additional mechanisms to counter viral invasions3,4. Here we identify 2,062 complete Cas9 loci, predict the structures of their associated proteins and reveal three structural growth trajectories for type II-C Cas9. We found that novel associated genes (NAGs) tended to be present within the loci of larger II-C Cas9s. Further investigation revealed that CbCas9 from Chryseobacterium species contains a novel β-REC2 domain, and forms a heterotetrameric complex with an NAG-encoded CRISPR-Cas-system-promoting (pro-CRISPR) protein of II-C Cas9 (PcrIIC1). The CbCas9-PcrIIC1 complex exhibits enhanced DNA binding and cleavage activity, broader compatibility for protospacer adjacent motif sequences, increased tolerance for mismatches and improved anti-phage immunity, compared with stand-alone CbCas9. Overall, our work sheds light on the diversity and 'growth evolutionary' trajectories of II-C Cas9 proteins at the structural level, and identifies many NAGs-such as PcrIIC1, which serves as a pro-CRISPR factor to enhance CRISPR-mediated immunity.

Pan-Cancer Single-Cell and Spatial-Resolved Profiling Reveals the Immunosuppressive Role of APOE+ Macrophages in Immune Checkpoint Inhibitor Therapy

The heterogeneity of macrophages influences the response to immune checkpoint inhibitor (ICI) therapy. However, few studies explore the impact of APOE+ macrophages on ICI therapy using single-cell RNA sequencing (scRNA-seq) and machine learning methods. The scRNA-seq and bulk RNA-seq data are Integrated to construct an M.Sig model for predicting ICI response based on the distinct molecular signatures of macrophage and machine learning algorithms. Comprehensive single-cell analysis as well as in vivo and in vitro experiments are applied to explore the potential mechanisms of the APOE+ macrophage in affecting ICI response. The M.Sig model shows clear advantages in predicting the efficacy and prognosis of ICI therapy in pan-cancer patients. The proportion of APOE+ macrophages is higher in ICI non-responders of triple-negative breast cancer compared with responders, and the interaction and longer distance between APOE+ macrophages and CD8+ exhausted T (Tex) cells affecting ICI response is confirmed by multiplex immunohistochemistry. In a 4T1 tumor-bearing mice model, the APOE inhibitor combined with ICI treatment shows the best efficacy. The M.Sig model using real-world immunotherapy data accurately predicts the ICI response of pan-cancer, which may be associated with the interaction between APOE+ macrophages and CD8+ Tex cells.

Heterogeneity generating capacity in tumorigenesis and cancer therapeutics

Cells of multicellular organisms generate heterogeneity in a controlled and transient fashion during embryogenesis, which can be reactivated in pathologies such as cancer. Although genomic heterogeneity is an important part of tumorigenesis, continuous generation of phenotypic heterogeneity is central for the adaptation of cancer cells to the challenges of tumorigenesis and response to therapy. Here I discuss the capacity of generating heterogeneity, hereafter called cell hetness, in cancer cells both as the activation of hetness oncogenes and inactivation of hetness tumor suppressor genes, which increase the generation of heterogeneity, ultimately producing an increase in adaptability and cell fitness. Transcriptomic high hetness states in therapy-tolerant cell states denote its importance in cancer resistance to therapy. The definition of the concept of hetness will allow the understanding of its origins, its control during embryogenesis, its loss of control in tumorigenesis and cancer therapeutics and its active targeting.

Coordinated chemokine expression defines macrophage subsets across tissues

Lung-resident macrophages, which include alveolar macrophages and interstitial macrophages (IMs), exhibit a high degree of diversity, generally attributed to different activation states, and often complicated by the influx of monocytes into the pool of tissue-resident macrophages. To gain a deeper insight into the functional diversity of IMs, here we perform comprehensive transcriptional profiling of resident IMs and reveal ten distinct chemokine-expressing IM subsets at steady state and during inflammation. Similar IM subsets that exhibited coordinated chemokine signatures and differentially expressed genes were observed across various tissues and species, indicating conserved specialized functional roles. Other macrophage types shared specific IM chemokine profiles, while also presenting their own unique chemokine signatures. Depletion of CD206hi IMs in Pf4creR26EYFP+DTR and Pf4creR26EYFPCx3cr1DTR mice led to diminished inflammatory cell recruitment, reduced tertiary lymphoid structure formation and fewer germinal center B cells in models of allergen- and infection-driven inflammation. These observations highlight the specialized roles of IMs, defined by their coordinated chemokine production, in regulating immune cell influx and organizing tertiary lymphoid tissue architecture.

Transient induction of actin cytoskeletal remodeling associated with dedifferentiation, proliferation, and redifferentiation stimulates cardiac regeneration

The formation of new and functional cardiomyocytes requires a 3-step process: dedifferentiation, proliferation, and redifferentiation, but the critical genes required for efficient dedifferentiation, proliferation, and redifferentiation remain unknown. In our study, a circular trajectory using single-nucleus RNA sequencing of the pericentriolar material 1 positive (PCM1+) cardiomyocyte nuclei from hearts 1 and 3 days after surgery-induced myocardial infarction (MI) on postnatal Day 1 was reconstructed and demonstrated that actin remodeling contributed to the dedifferentiation, proliferation, and redifferentiation of cardiomyocytes after injury. We identified four top actin-remodeling regulators, namely Tmsb4x, Tmsb10, Dmd, and Ctnna3, which we collectively referred to as 2D2P. Transiently expressed changes of 2D2P, using a polycistronic non-integrating lentivirus driven by Tnnt2 (cardiac-specific troponin T) promoters (Tnnt2-2D2P-NIL), efficiently induced transiently proliferative activation and actin remodeling in postnatal Day 7 cardiomyocytes and adult hearts. Furthermore, the intramyocardial delivery of Tnnt2-2D2P-NIL resulted in a sustained improvement in cardiac function without ventricular dilatation, thickened septum, or fatal arrhythmia for at least 4 months. In conclusion, this study highlights the importance of actin remodeling in cardiac regeneration and provides a foundation for new gene-cocktail-therapy approaches to improve cardiac repair and treat heart failure using a novel transient and cardiomyocyte-specific viral construct.