Tracking cell-type-specific temporal dynamics in human and mouse brains.

Progenitor cells are critical in preserving organismal homeostasis, yet their diversity and dynamics in the aged brain remain underexplored. We introduced TrackerSci, a single-cell genomic method that combines newborn cell labeling and combinatorial indexing to characterize the transcriptome and chromatin landscape of proliferating progenitor cells in vivo. Using TrackerSci, we investigated the dynamics of newborn cells in mouse brains across various ages and in a mouse model of Alzheimer's disease. Our dataset revealed diverse progenitor cell types in the brain and their epigenetic signatures. We further quantified aging-associated shifts in cell-type-specific proliferation and differentiation and deciphered the associated molecular programs. Extending our study to the progenitor cells in the aged human brain, we identified conserved genetic signatures across species and pinpointed region-specific cellular dynamics, such as the reduced oligodendrogenesis in the cerebellum. We anticipate that TrackerSci will be broadly applicable to unveil cell-type-specific temporal dynamics in diverse systems.

Reversible, tunable epigenetic silencing of TCF1 generates flexibility in the T cell memory decision.

The immune system encodes information about the severity of a pathogenic threat in the quantity and type of memory cells it forms. This encoding emerges from lymphocyte decisions to maintain or lose self-renewal and memory potential during a challenge. By tracking CD8+ T cells at the single-cell and clonal lineage level using time-resolved transcriptomics, quantitative live imaging, and an acute infection model, we find that T cells will maintain or lose memory potential early after antigen recognition. However, following pathogen clearance, T cells may regain memory potential if initially lost. Mechanistically, this flexibility is implemented by a stochastic cis-epigenetic switch that tunably and reversibly silences the memory regulator, TCF1, in response to stimulation. Mathematical modeling shows how this flexibility allows memory T cell numbers to scale robustly with pathogen virulence and immune response magnitudes. We propose that flexibility and stochasticity in cellular decisions ensure optimal immune responses against diverse threats.

A single-cell time-lapse of mouse prenatal development from gastrula to birth.

The house mouse (Mus musculus) is an exceptional model system, combining genetic tractability with close evolutionary affinity to humans1,2. Mouse gestation lasts only 3 weeks, during which the genome orchestrates the astonishing transformation of a single-cell zygote into a free-living pup composed of more than 500 million cells. Here, to establish a global framework for exploring mammalian development, we applied optimized single-cell combinatorial indexing3 to profile the transcriptional states of 12.4 million nuclei from 83 embryos, precisely staged at 2- to 6-hour intervals spanning late gastrulation (embryonic day 8) to birth (postnatal day 0). From these data, we annotate hundreds of cell types and explore the ontogenesis of the posterior embryo during somitogenesis and of kidney, mesenchyme, retina and early neurons. We leverage the temporal resolution and sampling depth of these whole-embryo snapshots, together with published data4-8 from earlier timepoints, to construct a rooted tree of cell-type relationships that spans the entirety of prenatal development, from zygote to birth. Throughout this tree, we systematically nominate genes encoding transcription factors and other proteins as candidate drivers of the in vivo differentiation of hundreds of cell types. Remarkably, the most marked temporal shifts in cell states are observed within one hour of birth and presumably underlie the massive physiological adaptations that must accompany the successful transition of a mammalian fetus to life outside the womb.

Single-cell, whole-embryo phenotyping of mammalian developmental disorders.

Mouse models are a critical tool for studying human diseases, particularly developmental disorders1. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse2. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing3 to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions4,5. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be 'decomposable' through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

The single-cell transcriptional landscape of mammalian organogenesis.

Mammalian organogenesis is a remarkable process. Within a short timeframe, the cells of the three germ layers transform into an embryo that includes most of the major internal and external organs. Here we investigate the transcriptional dynamics of mouse organogenesis at single-cell resolution. Using single-cell combinatorial indexing, we profiled the transcriptomes of around 2 million cells derived from 61 embryos staged between 9.5 and 13.5 days of gestation, in a single experiment. The resulting 'mouse organogenesis cell atlas' (MOCA) provides a global view of developmental processes during this critical window. We use Monocle 3 to identify hundreds of cell types and 56 trajectories, many of which are detected only because of the depth of cellular coverage, and collectively define thousands of corresponding marker genes. We explore the dynamics of gene expression within cell types and trajectories over time, including focused analyses of the apical ectodermal ridge, limb mesenchyme and skeletal muscle.

NGCN: Drug-target interaction prediction by integrating information and feature learning from heterogeneous network.

Drug-target interaction (DTI) prediction is essential for new drug design and development. Constructing heterogeneous network based on diverse information about drugs, proteins and diseases provides new opportunities for DTI prediction. However, the inherent complexity, high dimensionality and noise of such a network prevent us from taking full advantage of these network characteristics. This article proposes a novel method, NGCN, to predict drug-target interactions from an integrated heterogeneous network, from which to extract relevant biological properties and association information while maintaining the topology information. It focuses on learning the topology representation of drugs and targets to improve the performance of DTI prediction. Unlike traditional methods, it focuses on learning the low-dimensional topology representation of drugs and targets via graph-based convolutional neural network. NGCN achieves substantial performance improvements over other state-of-the-art methods, such as a nearly 1.0% increase in AUPR value. Moreover, we verify the robustness of NGCN through benchmark tests, and the experimental results demonstrate it is an extensible framework capable of combining heterogeneous information for DTI prediction.

Single-cell transcriptomics suggest distinct upstream drivers of IL-17A/F in hidradenitis versus psoriasis.

BACKGROUND: On the basis of the mounting evidence that type 17 T (T17) cells and increased IL-17 play a key role in driving hidradenitis suppurativa (HS) lesion development, biologic agents used previously in psoriasis that block signaling of IL-17A and/or IL-17F isoforms have been repurposed to treat HS. OBJECTIVE: Our research aimed to characterize the transcriptome of HS T17 cells compared to the transcriptome of psoriasis T17 cells, along with their ligand-receptor interactions with neighborhood immune cell subsets. METHODS: Single-cell data of 12,300 cutaneous immune cells from 8 deroofing surgical HS skin samples including dermal tunnels were compared to single-cell data of psoriasis skin (19,525 cells from 11 samples) and control skin (11,920 cells from 10 samples). All single-cell data were generated by the same protocol. RESULTS: HS T17 cells expressed lower levels of IL23R and higher levels of IL1R1 and IL17F compared to psoriasis T17 cells (P < .05). HS Treg cells expressed higher levels of IL1R1 and IL17F compared to psoriasis Treg cells (P < .05). Semimature dendritic cells were the major immune cell subsets expressing IL1B in HS, and IL-1ß ligand-receptor interactions between semimature dendritic cells and T17 cells were increased in HS compared to psoriasis (P < .05). HS dermal tunnel keratinocytes expressed inflammatory cytokines (IL17C, IL1A, IL1B, and IL6) that differed from the HS epidermis keratinocytes (IL36G) (P < .05). IL6, which synergizes with IL1B to maintain cytokine expression in T17 cells, was mainly expressed by fibroblasts in HS, which also expressed IL11+ inflammatory fibroblast genes (IL11, IL24, IL6, and POSTN) involved in the paracrine IL-1/IL-6 loop. CONCLUSION: The IL-1ß-T17 cell cytokine axis is likely a dominant pathway in HS with HS T17 cells activated by IL-1ß signaling, unlike psoriasis T17 cells, which are activated by IL-23 signaling.

A high-throughput single-cell RNA expression profiling method identifies human pericyte markers.

AIMS: We sought to identify and optimise a universally available histological marker for pericytes in the human brain. Such a marker could be a useful tool for researchers. Further, identifying a gene expressed relatively specifically in human pericytes could provide new insights into the biological functions of this fascinating cell type. METHODS: We analysed single-cell RNA expression profiles derived from different human and mouse brain regions using a high-throughput and low-cost single-cell transcriptome sequencing method called EasySci. Through this analysis, we were able to identify specific gene markers for pericytes, some of which had not been previously characterised. We then used commercially (and therefore universally) available antibodies to immunolabel the pericyte-specific gene products in formalin-fixed paraffin-embedded (FFPE) human brains and also performed immunoblots to determine whether appropriately sized proteins were recognised. RESULTS: In the EasySci data sets, highly pericyte-enriched expression was notable for SLC6A12 and SLC19A1. Antibodies against these proteins recognised bands of approximately the correct size in immunoblots of human brain extracts. Following optimisation of the immunohistochemical technique, staining for both antibodies was strongly positive in small blood vessels and was far more effective than a PDGFRB antibody at staining pericyte-like cells in FFPE human brain sections. In an exploratory sample of other human organs (kidney, lung, liver, muscle), immunohistochemistry did not show the same pericyte-like pattern of staining. CONCLUSIONS: The SLC6A12 antibody was well suited for labelling pericytes in human FFPE brain sections, based on the combined results of single-cell RNA-seq analyses, immunoblots and immunohistochemical studies.

DRGCNCDA: Predicting circRNA-disease interactions based on knowledge graph and disentangled relational graph convolutional network.

Emerging studies have shown that circular RNA (circRNA) plays a significant role in the diagnosis and prognosis of human disease. Some computational methods have been proposed to predict circRNA-disease associations. However, some methods only use circRNA-disease association and ignore the associations of other biological entities. In addition, these methods do not take into account the latent factors of different kinds of circRNAs and diseases. To solve these limitations of existing computational models, we propose a new computational model (DRGCNCDA) based on disentangled relational graph convolutional network. The circRNA-disease multi-relational graphs are constructed by collecting multiple relational data among circRNA, disease, miRNA and lncRNA. Then, the disentangled relational graph convolutional network is employed to obtain the feature vectors of circRNA and disease. Finally, knowledge graph model is applied to predict the affinity scores of circRNA-disease associations based on the embeddings of circRNA and disease. The 5-fold cross validation is utilized to evaluate the performance of the method. The experimental results show that the DRGCNCDA outperforms other existing models. Moreover, the case study demonstrates that the DRGCNCDA is effective to predict the circRNA-disease association and can provide reliable candidates for biological experiments.

Suppression of the HSF1-mediated proteotoxic stress response by the metabolic stress sensor AMPK.

Numerous extrinsic and intrinsic insults trigger the HSF1-mediated proteotoxic stress response (PSR), an ancient transcriptional program that is essential to proteostasis and survival under such conditions. In contrast to its well-recognized mobilization by proteotoxic stress, little is known about how this powerful adaptive mechanism reacts to other stresses. Surprisingly, we discovered that metabolic stress suppresses the PSR. This suppression is largely mediated through the central metabolic sensor AMPK, which physically interacts with and phosphorylates HSF1 at Ser121. Through AMPK activation, metabolic stress represses HSF1, rendering cells vulnerable to proteotoxic stress. Conversely, proteotoxic stress inactivates AMPK and thereby interferes with the metabolic stress response. Importantly, metformin, a metabolic stressor and popular anti-diabetic drug, inactivates HSF1 and provokes proteotoxic stress within tumor cells, thereby impeding tumor growth. Thus, these findings uncover a novel interplay between the metabolic stress sensor AMPK and the proteotoxic stress sensor HSF1 that profoundly impacts stress resistance, proteostasis, and malignant growth.

Similarity-guided graph contrastive learning for lncRNA-disease association prediction.

The increasing research evidence indicates that long non-coding RNAs (lncRNAs) play important roles in regulating biological processes and are closely associated with many human diseases. Computational methods have emerged as indispensable tools for identifying associations between long non-coding RNA (lncRNA) and diseases, primarily due to the time-consuming and costly nature of traditional biological experiments. Given the scarcity of verified lncRNA-disease associations, the intensifying focus on deep learning is playing a crucial role in refining the accuracy of predictive models. Moreover, the contrastive learning method exhibits a clear advantage in situations where data is scarce or annotation costs are high. In this paper, we leverage the advantages of graph neural networks and contrastive learning to innovatively propose a similarity-guided graph contrastive learning (SGGCL) model for predicting lncRNA-disease associations. In the SGGCL model, we employ a novel similarity-guided graph data augmentation method to generate high-quality positive and negative sample pairs, addressing the scarcity of verified data. Additionally, we utilize the RWR algorithm and a graph convolutional neural network for contrastive learning, facilitating the capture of global topology and high-level node embeddings. The experimental results on several datasets demonstrate the superior predictive performance and scalability of our method in lncRNA-disease association prediction compared to state-of-the-art methods.

A Panoramic View of Cell Population Dynamics in Mammalian Aging.

To elucidate the aging-associated cellular population dynamics throughout the body, here we present PanSci, a single-cell transcriptome atlas profiling over 20 million cells from 623 mouse tissue samples, encompassing a range of organs across different life stages, sexes, and genotypes. This comprehensive dataset allowed us to identify more than 3,000 unique cellular states and catalog over 200 distinct aging-associated cell populations experiencing significant depletion or expansion. Our panoramic analysis uncovered temporally structured, organ- and lineage-specific shifts of cellular dynamics during lifespan progression. Moreover, we investigated aging-associated alterations in immune cell populations, revealing both widespread shifts and organ-specific changes. We further explored the regulatory roles of the immune system on aging and pinpointed specific age-related cell population expansions that are lymphocyte-dependent. The breadth and depth of our 'cell-omics' methodology not only enhance our comprehension of cellular aging but also lay the groundwork for exploring the complex regulatory networks among varied cell types in the context of aging and aging-associated diseases.

PerturbSci-Kinetics: Dissecting key regulators of transcriptome kinetics through scalable single-cell RNA profiling of pooled CRISPR screens.

Here we described PerturbSci-Kinetics, a novel combinatorial indexing method for capturing three-layer single-cell readout (i.e., whole transcriptomes, nascent transcriptomes, sgRNA identities) across hundreds of genetic perturbations. Through PerturbSci-Kinetics profiling of pooled CRISPR screens targeting a variety of biological processes, we were able to decipher the complexity of RNA regulations at multiple levels (e.g., synthesis, processing, degradation), and revealed key regulators involved in miRNA and mitochondrial RNA processing pathways. Our technique opens the possibility of systematically decoding the genome-wide regulatory network underlying RNA temporal dynamics at scale and cost-effectively.

Dissecting key regulators of transcriptome kinetics through scalable single-cell RNA profiling of pooled CRISPR screens.

We present a combinatorial indexing method, PerturbSci-Kinetics, for capturing whole transcriptomes, nascent transcriptomes and single guide RNA (sgRNA) identities across hundreds of genetic perturbations at the single-cell level. Profiling a pooled CRISPR screen targeting various biological processes, we show the gene expression regulation during RNA synthesis, processing and degradation, miRNA biogenesis and mitochondrial mRNA processing, systematically decoding the genome-wide regulatory network that underlies RNA temporal dynamics at scale.

Hepatocellular carcinoma subtypes based on metabolic pathways reveals potential therapeutic targets.

Introduction: Hepatocellular carcinoma (HCC) is an aggressive malignancy with steadily increasing incidence rates worldwide and poor therapeutic outcomes. Studies show that metabolic reprogramming plays a key role in tumor genesis and progression. In this study, we analyzed the metabolic heterogeneity of epithelial cells in the HCC and screened for potential biomarkers. Methods: The hepatic single-cell RNA sequencing (scRNA-seq) datasets of HCC patients and healthy controls were obtained from the Gene Expression Omnibus (GEO) database. Based on data intergration and measurement of differences among groups, the metabolic epithelial cell subpopulations were identified. The single-cell metabolic pathway was analyzed and the myeloid subpopulations were identified. Cell-cell interaction analysis and single-cell proliferation analysis were performed. The gene expression profiles of HCC patients were obtained from the GSE14520 dataset of GEO and TCGA-LIHC cohort of the UCSC Xena website. Immune analysis was performed. The differentially expressed genes (DEGs) were identified and functionally annotated. Tumor tissues from HCC patients were probed with anti-ALDOA, anti-CD68, anti-CD163, anti-CD4 and anti-FOXP3 antibodies. Results We analyzed the scRNA-seq data from 48 HCC patients and 14 healthy controls. The epithelial cells were significantly enriched in HCC patients compared to the controls (p = 0.011). The epithelial cells from HCC patients were classified into two metabolism-related subpopulations (MRSs) - pertaining to amino acid metabolism (MRS1) and glycolysis (MRS2). Depending on the abundance of these metabolic subpopulations, the HCC patients were also classified into the MRS1 and MRS2 subtype distinct prognoses and immune infiltration. The MRS2 group had significantly worse clinical outcomes and more inflamed tumor microenvironment (TME), as well as a stronger crosstalk between MRS2 cells and immune subpopulations that resulted in an immunosuppressive TME. We also detected high expression levels of ALDOA in the MRS2 cells and HCC tissues. In the clinical cohort, HCC patients with higher ALDOA expression showed greater enrichment of immunosuppressive cells including M2 macrophages and T regulatory cells. Discussion: The glycolytic subtype of HCC cells with high ALDOA expression is associated with an immunosuppressive TME and predicts worse clinical outcomes, providing new insights into the metabolism and prognosis of HCC.

Optimized single-nucleus transcriptional profiling by combinatorial indexing.

Single-cell combinatorial indexing RNA sequencing (sci-RNA-seq) is a powerful method for recovering gene expression data from an exponentially scalable number of individual cells or nuclei. However, sci-RNA-seq is a complex protocol that has historically exhibited variable performance on different tissues, as well as lower sensitivity than alternative methods. Here, we report a simplified, optimized version of the sci-RNA-seq protocol with three rounds of split-pool indexing that is faster, more robust and more sensitive and has a higher yield than the original protocol, with reagent costs on the order of 1 cent per cell or less. The total hands-on time from nuclei isolation to final library preparation takes 2-3 d, depending on the number of samples sharing the experiment. The improvements also allow RNA profiling from tissues rich in RNases like older mouse embryos or adult tissues that were problematic for the original method. We showcase the optimized protocol via whole-organism analysis of an E16.5 mouse embryo, profiling ~380,000 nuclei in a single experiment. Finally, we introduce a 'Tiny-Sci' protocol for experiments in which input material is very limited.

Multi-omics segregate different transcriptomic impacts of anti-IL-17A blockade on type 17 T-cells and regulatory immune cells in psoriasis skin.

Durable psoriasis improvement has been reported in a subset of psoriasis patients after treatment withdrawal of biologics blocking IL-23/Type 17 T-cell (T17) autoimmune axis. However, it is not well understood if systemic blockade of the IL-23/T17 axis promotes immune tolerance in psoriasis skin. The purpose of the study was to find translational evidence that systemic IL-17A blockade promotes regulatory transcriptome modification in human psoriasis skin immune cell subsets. We analyzed human psoriasis lesional skin 6 mm punch biopsy tissues before and after systemic IL-17A blockade using the muti-genomics approach integrating immune cell-enriched scRNA-seq (n = 18), microarray (n = 61), and immunohistochemistry (n = 61) with repository normal control skin immune cell-enriched scRNA-seq (n = 10) and microarray (n = 8) data. For the T17 axis transcriptome, systemic IL-17A blockade depleted 100% of IL17A + T-cells and 95% of IL17F + T-cells in psoriasis skin. The expression of IL23A in DC subsets was also downregulated by IL-17A blockade. The expression of IL-17-driven inflammatory mediators (IL36G, S100A8, DEFB4A, and DEFB4B) in suprabasal keratinocytes was correlated with psoriasis severity and was downregulated by IL-17A blockade. For the regulatory DC transcriptome, the proportion of regulatory semimature DCs expressing regulatory DC markers of BDCA-3 (THBD) and DCIR (CLEC4A) was increased in posttreatment psoriasis lesional skin compared to pretreatment psoriasis lesional skin. In addition, IL-17A blockade induced higher expression of CD1C and CD14, which are markers of CD1c+ CD14+ dendritic cell (DC) subset that suppresses antigen-specific T-cell responses, in posttreatment regulatory semimature DCs compared to pretreatment regulatory semimature DCs. In conclusion, systemic IL-17A inhibition not only blocks the entire IL-23/T17 cell axis but also promotes regulatory gene expression in regulatory DCs in human psoriasis skin.

A global view of aging and Alzheimer's pathogenesis-associated cell population dynamics and molecular signatures in human and mouse brains.

Conventional methods fall short in unraveling the dynamics of rare cell types related to aging and diseases. Here we introduce EasySci, an advanced single-cell combinatorial indexing strategy for exploring age-dependent cellular dynamics in the mammalian brain. Profiling approximately 1.5 million single-cell transcriptomes and 400,000 chromatin accessibility profiles across diverse mouse brains, we identified over 300 cell subtypes, uncovering their molecular characteristics and spatial locations. This comprehensive view elucidates rare cell types expanded or depleted upon aging. We also investigated cell-type-specific responses to genetic alterations linked to Alzheimer's disease, identifying associated rare cell types. Additionally, by profiling 118,240 human brain single-cell transcriptomes, we discerned cell- and region-specific transcriptomic changes tied to Alzheimer's pathogenesis. In conclusion, this research offers a valuable resource for probing cell-type-specific dynamics in both normal and pathological aging.