Single-cell mRNA quantification and differential analysis with Census

Prognostic iron-metabolism signature robustly stratifies single-cell characteristics of hepatocellular carcinoma

Cancer immunotherapy has shown to be a promising method in treating hepatocellular carcinoma (HCC), but suboptimal responses in patients are attributed to cellular and molecular heterogeneity. Iron metabolism-related genes (IRGs) are important in maintaining immune system homeostasis and have the potential to help develop new strategies for HCC treatment. Herein, we constructed and validated the iron-metabolism gene prognostic index (IPX) using univariate Cox proportional hazards regression and LASSO Cox regression analysis, successfully categorizing HCC patients into two groups with distinct survival risks. Then, we performed single-sample gene set enrichment analysis, weighted correlation network analysis, gene ontology enrichment analysis, cellular lineage analysis, and SCENIC analysis to reveal the key determinants underlying the ability of this model based on bulk and single-cell transcriptomic data. We identified several driver transcription factors specifically activated in specific malignant cell sub-populations to contribute to the adverse survival outcomes in the IPX-high subgroup. Within the tumor microenvironment (TME), T cells displayed significant diversity in their cellular characteristics and experienced changes in their developmental paths within distinct clusters identified by IPX. Interestingly, the proportion of Treg cells was increased in the high-risk group compared with the low-risk group. These results suggest that iron-metabolism could be involved in reshaping the TME, thereby disrupting the cell cycle of immune cells. This study utilized IRGs to construct a novel and reliable model, which can be used to assess the prognosis of patients with HCC and further clarify the molecular mechanisms of IRGs in HCC at single-cell resolution.

NTCdb: Single-cell transcriptome database of human inflammatory-associated diseases

With both the advancement of technology and the decline in costs, single-cell transcriptomics sequencing has become widespread in the biomedical area in recent years. It can facilitate the pathogenic characteristics at the single-cell level, which will assist clinical researchers in exploring the mechanism of diseases. As a result, single-cell transcriptome data based on clinical samples grew exponentially. However, there is still a lack of a comprehensive database about immunocytes in inflammatory-associated diseases. To address this deficiency, we propose a human inflammatory-associated disease-based single-cell transcriptome database, NTCdb (www.ntcdb.org.cn). NTCdb integrates the open-source data of 1,023,166 cells derived from 11 tissues of 17 inflammatory-associated diseases in a uniform pipeline. It provides a set of analyzing results, including cell communication analysis, enrichment analysis, and Pseudo-Time analysis, to obtain various characteristics of immune cells in inflammatory-associated disease. Taking COVID-19 as a case study, NTCdb displays important information including potentially significant functions of certain cells, genes, and signaling pathways, as well as the commonalities of specific immunocytes between different inflammatory-associated disease.

Single-cell transcriptomics of blood identified IFIT1+ neutrophil subcluster expansion in NTM-PD patients

OBJECTIVE

Non-tuberculous mycobacterial pulmonary disease (NTM-PD) is caused by an imbalance between pathogens and impaired host immune responses. Mycobacterium avium complex (MAC) and Mycobacterium abscessus (MAB) are the two major pathogens that cause NTM-PD. In this study, we sought to dissect the transcriptomes of peripheral blood immune cells at the single-cell resolution in NTM-PD patients and explore potential clinical markers for NTM-PD diagnosis and treatment.

METHODS

Peripheral blood samples were collected from six NTM-PD patients, including three MAB-PD patients, three MAC-PD patients, and two healthy controls. We employed single-cell RNA sequencing (scRNA-seq) to define the transcriptomic landscape at a single-cell resolution. A comprehensive scRNA-seq analysis was performed, and flow cytometry was conducted to validate the results of scRNA-seq.

RESULTS

A total of 27,898 cells were analyzed. Nine T-cells, six mononuclear phagocytes (MPs), and four neutrophil subclusters were defined. During NTM infection, naïve T-cells were reduced, and effector T-cells increased. High cytotoxic activities were shown in T-cells of NTM-PD patients. The proportion of inflammatory and activated MPs subclusters was enriched in NTM-PD patients. Among neutrophil subclusters, an IFIT1+ neutrophil subcluster was expanded in NTM-PD compared to healthy controls. This suggests that IFIT1+ neutrophil subcluster might play an important role in host defense against NTM. Functional enrichment analysis of this subcluster suggested that it is related to interferon response. Cell-cell interaction analysis revealed enhanced CXCL8-CXCR1/2 interactions between the IFIT1+ neutrophil subcluster and NK cells, NKT cells, classical mononuclear phagocytes subcluster 1 (classical Mo1), classical mononuclear phagocytes subcluster 2 (classical Mo2) in NTM-PD patients compared to healthy controls.

CONCLUSIONS

Our data revealed disease-specific immune cell subclusters and provided potential new targets of NTM-PD. Specific expansion of IFIT1+ neutrophil subclusters and the CXCL8-CXCR1/2 axis may be involved in the pathogenesis of NTM-PD. These insights may have implications for the diagnosis and treatment of NTM-PD.

Single-Cell Analyses Reveal the Metabolic Heterogeneity and Plasticity of the Tumor Microenvironment during Head and Neck Squamous Cell Carcinoma Progression

Metabolic reprogramming is a hallmark of cancer. In addition to metabolic alterations in the tumor cells, multiple other metabolically active cell types in the tumor microenvironment (TME) contribute to the emergence of a tumor-specific metabolic milieu. Here, we defined the metabolic landscape of the TME during the progression of head and neck squamous cell carcinoma (HNSCC) by performing single-cell RNA sequencing on 26 human patient specimens, including normal tissue, precancerous lesions, early stage cancer, advanced-stage cancer, lymph node metastases, and recurrent tumors. The analysis revealed substantial heterogeneity at the transcriptional, developmental, metabolic, and functional levels in different cell types. SPP1+ macrophages were identified as a protumor and prometastatic macrophage subtype with high fructose and mannose metabolism, which was further substantiated by integrative analysis and validation experiments. An inhibitor of fructose metabolism reduced the proportion of SPP1+ macrophages, reshaped the immunosuppressive TME, and suppressed tumor growth. In conclusion, this work delineated the metabolic landscape of HNSCC at a single-cell resolution and identified fructose metabolism as a key metabolic feature of a protumor macrophage subpopulation. Significance: Fructose and mannose metabolism is a metabolic feature of a protumor and prometastasis macrophage subtype and can be targeted to reprogram macrophages and the microenvironment of head and neck squamous cell carcinoma.

Compound models and Pearson residuals for single-cell RNA-seq data without UMIs

Recent work employed Pearson residuals from Poisson or negative binomial models to normalize UMI data. To extend this approach to non-UMI data, we model the additional amplification step with a compound distribution: we assume that sequenced RNA molecules follow a negative binomial distribution, and are then replicated following an amplification distribution. We show how this model leads to compound Pearson residuals, which yield meaningful gene selection and embeddings of Smart-seq2 datasets. Further, we suggest that amplification distributions across several sequencing protocols can be described by a broken power law. The resulting compound model captures previously unexplained overdispersion and zero-inflation patterns in non-UMI data.

Single-cell RNA sequencing reveals dynamics of gene expression for 2D elongation and 3D growth in Physcomitrium patens

The transition from two-dimensional (2D) to 3D growth likely facilitated plants to colonize land, but its heterogeneity is not well understood. In this study, we utilized single-cell RNA sequencing to analyze the moss Physcomitrium patens, whose morphogenesis involves a transition from 2D to 3D growth. We profiled over 17,000 single cells covering all major vegetative tissues, including 2D filaments (chloronema and caulonema) and 3D structures (bud and gametophore). Pseudotime analyses revealed larger numbers of candidate genes that determine cell fates for 2D tip elongation or 3D bud differentiation. Using weighted gene co-expression network analysis, we identified a module that connects β-type carbonic anhydrases (βCAs) with auxin. We further validated the cellular expression patterns of βCAs and demonstrated their roles in 3D gametophore development. Overall, our study provides insights into cellular heterogeneity in a moss and identifies molecular signatures that underpin the 2D-to-3D growth transition at single-cell resolution.

The architecture of silk-secreting organs during the final larval stage of silkworms revealed by single-nucleus and spatial transcriptomics

Natural silks are renewable proteins with impressive mechanical properties and biocompatibility that are useful in various fields. However, the cellular and spatial organization of silk-secreting organs remains unclear. Here, we combined single-nucleus and spatially resolved transcriptomics to systematically map the cellular and spatial composition of the silk glands (SGs) of mulberry silkworms late in larval development. This approach allowed us to profile SG cell types and cell state dynamics and identify regulatory networks and cell-cell communication related to efficient silk protein synthesis; key markers were validated via transgenic approaches. Notably, we demonstrated the indispensable role of the ecdysone receptor (ultraspiracle) in regulating endoreplication in SG cells. Our atlas presents the results of spatiotemporal analysis of silk-secreting organ architecture late in larval development; this atlas provides a valuable reference for elucidating the mechanism of efficient silk protein synthesis and developing sustainable products made from natural silk.

Advancing bladder cancer management: development of a prognostic model and personalized therapy

Background

Bladder cancer (BLCA) was recognized as a significant public health challenge due to its high incidence and mortality rates. The influence of molecular subtypes on treatment outcomes was well-acknowledged, necessitating further exploration of their characterization and application. This study was aimed at enhancing the understanding of BLCA by mapping its molecular heterogeneity and developing a robust prognostic model using single-cell and bulk RNA sequencing data. Additionally, immunological characteristics and personalized treatment strategies were investigated through the risk score.

Methods

Single-cell RNA sequencing (scRNA-seq) data from GSE135337 and bulk RNA-seq data from several sources, including GSE13507, GSE31684, GSE32894, GSE69795, and TCGA-BLCA, were utilized. Molecular subtypes, particularly the basal-squamous (Ba/Sq) subtype associated with poor prognosis, were identified. A prognostic model was constructed using LASSO and Cox regression analyses focused on genes linked with the Ba/Sq subtype. this model was validated across internal and external datasets to ensure predictive accuracy. High- and low-risk groups based on the risk score derived from TCGA-BLCA data were analyzed to examine their immune-related molecular profiles and treatment responses.

Results

Six molecular subtypes were identified, with the Ba/Sq subtype being consistently associated with poor prognosis. The prognostic model, based on basal-squamous subtype-related genes (BSSRGs), was shown to have strong predictive performance across diverse clinical settings with AUC values at 1, 3, and 5 years indicating robust predictability in training, testing, and entire datasets. Analysis of the different risk groups revealed distinct immune infiltration and microenvironments. Generally higher tumor mutation burden (TMB) scores and lower tumor immune dysfunction and exclusion (TIDE) scores were exhibited by the low-risk group, suggesting varied potentials for systemic drug response between the groups. Finally, significant differences in potential systemic drug response rates were also observed between risk groups.

Conclusions

The study introduced and validated a new prognostic model for BLCA based on BSSRGs, which was proven effective in prognosis prediction. The potential for personalized therapy, optimized by patient stratification and immune profiling, was highlighted by our risk score, aiming to improve treatment efficacy. This approach was promised to offer significant advancements in managing BLCA, tailoring treatments based on detailed molecular and immunological insights.

Single-nucleus transcriptomics uncovers a geroprotective role of YAP in primate gingival aging

Aging has a profound impact on the gingiva and significantly increases its susceptibility to periodontitis, a worldwide prevalent inflammatory disease. However, a systematic characterization and comprehensive understanding of the regulatory mechanism underlying gingival aging is still lacking. Here, we systematically dissected the phenotypic characteristics of gingiva during aging in primates and constructed the first single-nucleus transcriptomic landscape of gingival aging, by which a panel of cell type-specific signatures were elucidated. Epithelial cells were identified as the most affected cell types by aging in the gingiva. Further analyses pinpointed the crucial role of YAP in epithelial self-renew and homeostasis, which declined during aging in epithelial cells, especially in basal cells. The decline of YAP activity during aging was confirmed in the human gingival tissues, and downregulation of YAP in human primary gingival keratinocytes recapitulated the major phenotypic defects observed in the aged primate gingiva while overexpression of YAP showed rejuvenation effects. Our work provides an in-depth understanding of gingival aging and serves as a rich resource for developing novel strategies to combat aging-associated gingival diseases, with the ultimate goal of advancing periodontal health and promoting healthy aging.

Advancing laser ablation assessment in hyperspectral imaging through machine learning

Hyperspectral imaging (HSI) is gaining increasing relevance in medicine, with an innovative application being the intraoperative assessment of the outcome of laser ablation treatment used for minimally invasive tumor removal. However, the high dimensionality and complexity of HSI data create a need for end-to-end image processing workflows specifically tailored to handle these data. This study addresses this challenge by proposing a multi-stage workflow for the analysis of hyperspectral data and allows investigating the performance of different components and modalities for ablation detection and segmentation. To address dimensionality reduction, we integrated principal component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE) to capture dominant variations and reveal intricate structures, respectively. Additionally, we employed the Faster Region-based Convolutional Neural Network (Faster R-CNN) to accurately localize ablation areas. The two-stage detection process of Faster R-CNN, along with the choice of dimensionality reduction technique and data modality, significantly influenced the performance in detecting ablation areas. The evaluation of the ablation detection on an independent test set demonstrated a mean average precision of approximately 0.74, which validates the generalization ability of the models. In the segmentation component, the Mean Shift algorithm showed high quality segmentation without manual cluster definition. Our results prove that the integration of PCA, t-SNE, and Faster R-CNN enables improved interpretation of hyperspectral data, leading to the development of reliable ablation detection and segmentation systems.

Synaptic plasticity in human thalamocortical assembloids

Synaptic plasticities, such as long-term potentiation (LTP) and depression (LTD), tune synaptic efficacy and are essential for learning and memory. Current studies of synaptic plasticity in humans are limited by a lack of adequate human models. Here, we modeled the thalamocortical system by fusing human induced pluripotent stem cell-derived thalamic and cortical organoids. Single-nucleus RNA sequencing revealed that >80% of cells in thalamic organoids were glutamatergic neurons. When fused to form thalamocortical assembloids, thalamic and cortical organoids formed reciprocal long-range axonal projections and reciprocal synapses detectable by light and electron microscopy, respectively. Using whole-cell patch-clamp electrophysiology and two-photon imaging, we characterized glutamatergic synaptic transmission. Thalamocortical and corticothalamic synapses displayed short-term plasticity analogous to that in animal models. LTP and LTD were reliably induced at both synapses; however, their mechanisms differed from those previously described in rodents. Thus, thalamocortical assembloids provide a model system for exploring synaptic plasticity in human circuits.

Single-cell RNA Sequencing Identifies a Novel Subtype of Microglia with High Cd74 Expression that Facilitates White Matter Inflammation During Chronic Cerebral Hypoperfusion

Vascular dementia (VaD) causes progressive cognitive decline in the elderly population, but there is short of available therapeutic measures. Microglia-mediated neuroinflammation is vigorously involved in the pathogenesis of VaD, but the traditional classification of microglial M1/M2 phenotypes remains restrictive and controversial. This study aims to investigate whether microglia transform into novel subtypes in VaD. Chronic cerebral hypoperfusion (CCH) rat model was constructed to mimic VaD. Microglia were isolated via magnetic-activated cell sorting and analyzed by single-cell RNA sequencing (scRNA-seq) and bioinformatics. The findings inferred from scRNA-seq and bioinformatics were further validated through in vivo experiments. In this study, microglia were divided into eight clusters. The proportion of MG5 cluster was significantly increased in the white matter of the CCH group compared with the Sham group and was named chronic ischemia-associated microglia (CIAM). Immunity- and inflammation-related genes, including RT1-Db1, RT1-Da, RT1-Ba, Cd74, Spp1, C3, and Cd68, were markedly upregulated in CIAM. Enrichment analysis illustrated that CIAM possessed the function of evoking neuroinflammation. Further studies unveiled that Cd74 is associated with the most abundant GO terms involved in inflammation as well as cell proliferation and differentiation. In addition, microglia-specific Cd74 knockdown mediated by adeno-associated virus decreased the abundance of CIAM in the white matter, thereby mitigating inflammatory cytokine levels, alleviating white matter lesions, and improving cognitive impairment for CCH rats. These findings indicate that Cd74 is the core molecule of CIAM to trigger neuroinflammation and induce microglial differentiation to CIAM, suggesting that Cd74 may be a potential therapeutic target for VaD.

Key factor screening in mouse NASH model using single-cell sequencing combined with machine learning

Aims

To identify and analyze genes closely related to the progression of nonalcoholic steatohepatitis (NASH) by employing a combination of single-cell RNA sequencing and machine-learning algorithms.

Main methods

Single-cell RNA sequencing (scRNA-seq) analysis was performed to find the cell population with the most significant differences between the Chow and NASH groups. This approach was used to validate the developmental trajectory of this cell population and investigate changes in cellular communication and important signaling pathways among these cells. Subsequently, high dimensional Weighted Gene Co-expression Network Analysis (hdWGCNA) was used to find the key modules in NASH. Machine learning analyses were performed to further identify core genes. Deep learning techniques were applied to elucidate the correlation between core genes and immune cells. The accuracy of this correlation was further confirmed using deep learning techniques, specifically Convolutional Neural Networks.

Key findings

By comparing scRNA-seq data between the Chow and NASH groups, we have observed a notable distinction existing in the Kupffer cell population. Signaling interactions between hepatic macrophages and other cells were significantly heightened in the NASH group. Through subsequent analysis of macrophage subtypes and key modules, we identified 150 genes tightly associated with NASH. Finally, we highlighted the 16 most significant core genes using multiple iterations of machine learning. Furthermore, we pointed out the close relationship between core genes and immune cells.

Significances

Using scRNA-seq analysis and machine learning, we can distinguish NASH-related genes from large genetic datasets, providing theoretical support in finding potential targets for the development of novel therapies.

Histone lactylation-related genes correlate with the molecular patterns and functions of cancer-associated fibroblasts and have significant clinical implications in clear cell renal cell carcinoma

Recent research emphasised the indispensable role of histone lactylation in the activation of hepatic stellate cells. The VHL mutation is extremely common in clear cell renal cell carcinoma, which normally causes a metabolic shift in cancer cells and increases lactate production, eventually creating a lactate-enriched tumour microenvironment. Cancer-associated fibroblasts (CAFs) promote tumour progression, which is also vital in clear cell renal cell carcinoma. Therefore, this study investigated histone lactylation in CAFs and its impact on patient survival. Multiomics technology was employed to determine the role of histone lactylation-related genes in the evolution of CAFs which correlated with the function and molecular signatures of CAFs. The results suggested that TIMP1 was the hub gene of histone lactylation-related genes in clear cell renal cell carcinoma.

Destabilization of mRNAs enhances competence to initiate meiosis in mouse spermatogenic cells

The specialized cell cycle of meiosis transforms diploid germ cells into haploid gametes. In mammals, diploid spermatogenic cells acquire the competence to initiate meiosis in response to retinoic acid. Previous mouse studies revealed that MEIOC interacts with RNA-binding proteins YTHDC2 and RBM46 to repress mitotic genes and to promote robust meiotic gene expression in spermatogenic cells that have initiated meiosis. Here, we have used the enhanced resolution of scRNA-seq and bulk RNA-seq of developmentally synchronized spermatogenesis to define how MEIOC molecularly supports early meiosis in spermatogenic cells. We demonstrate that MEIOC mediates transcriptomic changes before meiotic initiation, earlier than previously appreciated. MEIOC, acting with YTHDC2 and RBM46, destabilizes its mRNA targets, including the transcriptional repressors E2f6 and Mga, in mitotic spermatogonia. MEIOC thereby derepresses E2F6- and MGA-repressed genes, including Meiosin and other meiosis-associated genes. This confers on spermatogenic cells the molecular competence to, in response to retinoic acid, fully activate the transcriptional regulator STRA8-MEIOSIN, which is required for the meiotic G1/S phase transition and for meiotic gene expression. We conclude that, in mice, mRNA decay mediated by MEIOC-YTHDC2-RBM46 enhances the competence of spermatogenic cells to initiate meiosis.

Short-chain fatty acids regulate erastin-induced cardiomyocyte ferroptosis and ferroptosis-related genes

Background

Ferroptosis has been proven to contribute to the progression of myocardial ischemia/reperfusion (I/R) injury and can be inhibited or promoted by ATF3. Short-chain fatty acids (SCFAs) have shown benefits in various cardiovascular diseases with anti-inflammatory and antioxidant effects. However, the impact of SCFAs on ferroptosis in ischemic-stimulated cardiomyocytes remains unknown. This study aimed to investigate the effect of SCFAs on cardiomyocyte ferroptosis, the expression of ATF3, and its potential upstream regulators.

Methods and results

The expression of ATF3, ferroptosis pathway geneset (FPG), and geneset of potential regulators for ATF3 (GPRA, predicted by the PROMO database) was explored in the public human myocardial infarction single-cell RNA-seq (sma) dataset. Cardiomyocyte data was extracted from the dataset and re-clustered to explore the FPG, ATF3, and GPRA expression patterns in cardiomyocyte subclusters. A dose-dependent toxic experiment was run to detect the suitable dose for SCFA treatment. The erastin-induced ferroptosis model and hypoxia-reoxygenation (H/R) model (10 h of hypoxia followed by 6 h of reoxygenation) were adopted to assess the effect of SCFAs via the CCK8 assay. Gene expression was examined via RT-PCR and western blot. Ferroptosis markers, including lipid peroxides and Fe2+, were detected using the liperfluo and ferroOrange probes, respectively. In the sma dataset, upregulated ferroptosis pathway genes were mainly found in the infarction-stimulated cardiac cells (border zone and fibrotic zone), particularly the cardiomyocytes and adipocytes. The ATF3 and some of its potential transcription factors (VDR, EGR3, PAX5, and SP1) can be regulated by SCFA. SCFA can attenuate erastin-induced lipid peroxidation in cardiomyocytes. SCFA treatment can also reverse erastin-induced Fe2+ increase but may strengthen the Fe2+ in the H/R model. We also precisely defined a ferroptosis subcluster of cardiomyocytes (CM09) that highly expressed FPG, ATF3, and GPRA.

Conclusion

The ATF3 and the ferroptosis pathway are elevated in cardiomyocytes of injury-related cardiac regions (border zone, ischemic zone, and fibrotic zone). SCFA can attenuate cardiomyocyte ferroptosis and regulate the expression of ATF3. Our study offers novel insights into the potential targets of SCFAs in the cardiovascular system.

SIGLEC15, negatively correlated with PD-L1 in HCC, could induce CD8+ T cell apoptosis to promote immune evasion

Functional roles of SIGLEC15 in hepatocellular carcinoma (HCC) were not clear, which was recently found to be an immune inhibitor with similar structure of inhibitory B7 family members. SIGLEC15 expression in HCC was explored in public databases and further examined by PCR analysis. SIGLEC15 and PD-L1 expression patterns were examined in HCC samples through immunohistochemistry. SIGLEC15 expression was knocked-down or over-expressed in HCC cell lines, and CCK8 tests were used to examine cell proliferative ability in vitro. Influences of SIGLEC15 expression on tumor growth were examined in immune deficient and immunocompetent mice respectively. Co-culture system of HCC cell lines and Jurkat cells, flow cytometry analysis of tumor infiltrated immune cells and further sequencing analyses were performed to investigate how SIGLEC15 could affect T cells in vitro and in vivo. We found SIGLEC15 was increased in HCC tumor tissues and was negatively correlated with PD-L1 in HCC samples. In vitro and in vivo models demonstrated inhibition of SIGLEC15 did not directly influence tumor proliferation. However, SIGLEC15 could promoted HCC immune evasion in immune competent mouse models. Knock-out of Siglec15 could inhibit tumor growth and reinvigorate CD8+ T cell cytotoxicity. Anti-SIGLEC15 treatment could effectively inhibit tumor growth in mouse models with or without mononuclear phagocyte deletion. Bulk and single-cell RNA sequencing data of treated mouse tumors demonstrated SIGLEC15 could interfere CD8+ T cell viability and induce cell apoptosis. In all, SIGLEC15 was negatively correlated with PD-L1 in HCC and mainly promote HCC immune evasion through inhibition of CD8+ T cell viability and cytotoxicity.

Analysis of single-cell CRISPR perturbations indicates that enhancers act multiplicatively and provides limited evidence for epistatic-like interactions

A single gene may have multiple enhancers, but how they work in concert to regulate transcription is poorly understood. To analyze enhancer interactions throughout the genome, we developed a generalized linear modeling framework, GLiMMIRS, for interrogating enhancer effects from single-cell CRISPR experiments. We applied GLiMMIRS to a published dataset and tested for interactions between 46,166 enhancer pairs and corresponding genes, including 264 'high-confidence' enhancer pairs. We found that enhancer effects combine multiplicatively but with limited evidence for further interactions. Only 31 enhancer pairs exhibited significant interactions (FDR < 0.1), of which none came from the high confidence subset and 20 were driven by outlier expression values. Additional analyses of a second CRISPR dataset and in silico enhancer perturbations with Enformer both support a multiplicative model of enhancer effects without interactions. Altogether, our results indicate that enhancer interactions are uncommon or have small effects that are difficult to detect.

Single-cell RNA sequencing reveals special basal cells and fibroblasts in idiopathic pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is the most predominant type of idiopathic interstitial pneumonia and has an increasing incidence, poor prognosis, and unclear pathogenesis. In order to investigate the molecular mechanisms underlying IPF further, we performed single-cell RNA sequencing analysis on three healthy controls and five IPF lung tissue samples. The results revealed a significant shift in epithelial cells (ECs) phenotypes in IPF, which may be attributed to the differentiation of alveolar type 2 cells to basal cells. In addition, several previously unrecognized basal cell subtypes were preliminarily identified, including extracellular matrix basal cells, which were increased in the IPF group. We identified a special population of fibroblasts that highly expressed extracellular matrix-related genes, POSTN, CTHRC1, COL3A1, COL5A2, and COL12A1. We propose that the close interaction between ECs and fibroblasts through ligand-receptor pairs may have a critical function in IPF development. Collectively, these outcomes provide innovative perspectives on the complexity and diversity of basal cells and fibroblasts in IPF and contribute to the understanding of possible mechanisms in pathological lung fibrosis.

Single-cell RNA sequencing uncovers cellular heterogeneity and the progression of heterotopic ossification of the elbow

Heterotopic ossification of the elbow (HOE) is a complicated pathologic process characterized by extra bone formation in the elbow. Bone formation is a complex developmental process involving the differentiation of mesenchymal stem cells into osteoblasts. The aim of this study was to explore the cellular origin and progression of HOE by single-cell RNA sequencing. We identified 13 clusters of cells in HOE and further analyzed the subclusters for 4 of the main cell types. Six subclusters of osteoblasts, nine subclusters of chondrocytes, six subclusters of fibroblasts, and five subclusters of mononuclear phagocytes (MPs) were identified and analyzed. The new findings on osterix (OSX) and SOX9 expression in osteoblast subclusters and chondrocyte subclusters indicate that HOE is mediated through endochondral ossification. Further identification of the corresponding signature gene sets of distinct subclusters indicated that subclusters of osteoblasts_3, osteoblasts_4, osteoblasts_5, and osteoblasts_6 are relatively more mature during the osteoblastic progression of HOE. The trajectory analysis of the osteoblasts demonstrated that some genes were gradually downregulated, such as CRYAB, CCL3, SFRP4, WIF1, and IGFBP3, while other critical genes were upregulated, such as VCAN, IGFBP4, FSTL1, POSTN, MDK, THBS2, and ALPL, suggesting that these factors may participate in HOE progression. Cell-cell communication networks revealed extensive molecular interactions among the 13 HOE clusters. Ligand-receptor pairs for IL6, COL24A1, COL22A1, VWF, FZD6, FGF2, and NOTCH1 were identified, suggesting that multiple signaling pathways may be involved in HOE. In conclusion, this study provided the cellular atlas for HOE. We have established a greater extent of the heterogeneity of HOE cells than previously known through transcriptomic analysis at the single-cell level. We have observed gradual patterns of signature gene expression during the differentiation and maturation progression of osteoblasts from stem cells in HOE with higher resolution. The cell heterogeneity of HOE deserves further investigation to pave the way for identification of potential targets for HOE early diagnosis and therapeutic treatment.

Prostate cancer cancer-associated fibroblasts with stable markers post-androgen deprivation therapy associated with tumor progression and castration resistant prostate cancer

The specificity and clinical relevance of cancer-associated fibroblasts (CAFs) in prostate cancer (PCa), as well as the effect of androgen deprivation therapy (ADT) on CAFs, remain to be fully elucidated. Using cell lineage diversity and weighted gene co-expression network analysis (WGCNA), we pinpointed a unique CAF signature exclusive to PCa. The specificity of this CAF signature was validated through single-cell RNA sequencing (scRNA-seq), cell line RNA sequencing, and immunohistochemistry. This signature associates CAFs with tumor progression, elevated Gleason scores, and the emergence of castration resistant prostate cancer (CRPC). Using scRNA-seq on collected samples, we demonstrated that the CAF-specific signature is not altered by ADT, maintaining its peak signal output. Identifying a PCa-specific CAF signature and observing signaling changes in CAFs after ADT lay essential groundwork for further PCa studies.

Stepwise developmental mimicry generates proximal-biased kidney organoids

The kidney maintains body fluid homeostasis by reabsorbing essential compounds and excreting waste. Proximal tubule cells, crucial for renal reabsorption of a range of sugars, ions, and amino acids, are highly susceptible to damage, leading to pathologies necessitating dialysis and kidney transplants. While human pluripotent stem cell-derived kidney organoids are used for modeling renal development, disease, and injury, the formation of proximal nephron cells in these 3D structures is incomplete. Here, we describe how to drive the development of proximal tubule precursors in kidney organoids by following a blueprint of in vivo human nephrogenesis. Transient manipulation of the PI3K signaling pathway activates Notch signaling in the early nephron and drives nephrons toward a proximal precursor state. These "proximal-biased" (PB) organoid nephrons proceed to generate proximal nephron precursor cells. Single-cell transcriptional analyses across the organoid nephron differentiation, comparing control and PB types, confirm the requirement of transient Notch signaling for proximal development. Indicative of functional maturity, PB organoids demonstrate dextran and albumin uptake, akin to in vivo proximal tubules. Moreover, PB organoids are highly sensitive to nephrotoxic agents, display an injury response, and drive expression of HAVCR1 / KIM1 , an early proximal-specific marker of kidney injury. Injured PB organoids show evidence of collapsed tubules, DNA damage, and upregulate the injury-response marker SOX9 . The PB organoid model therefore has functional relevance and potential for modeling mechanisms underpinning nephron injury. These advances improve the use of iPSC-derived kidney organoids as tools to understand developmental nephrology, model disease, test novel therapeutics, and for understanding human renal physiology.

Mitotic catastrophe heterogeneity: implications for prognosis and immunotherapy in hepatocellular carcinoma

Background and aims

The mitotic catastrophe (MC) pathway plays an important role in hepatocellular carcinoma (HCC) progression and tumor microenvironment (TME) regulation. However, the mechanisms linking MC heterogeneity to immune evasion and treatment response remain unclear.

Methods

Based on 94 previously published highly correlated genes for MC, HCC patients' data from the Cancer Genome Atlas (TCGA) and changes in immune signatures and prognostic stratification were studied. Time and spatial-specific differences for MCGs were assessed by single-cell RNA sequencing and spatial transcriptome (ST) analysis. Multiple external databases (GEO, ICGC) were employed to construct an MC-related riskscore model.

Results

Identification of two MC-related subtypes in HCC patients from TCGA, with clear differences in immune signatures and prognostic risk stratification. Spatial mapping further associates low MC tumor regions with significant immune escape-related signaling. Nomogram combining MC riskscore and traditional indicators was validated great effect for early prediction of HCC patient outcomes.

Conclusion

MC heterogeneity enables immune escape and therapy resistance in HCC. The MC gene signature serves as a reliable prognostic indicator for liver cancer. By revealing clear immune and spatial heterogeneity of HCC, our integrated approach provides contextual therapeutic strategies for optimal clinical decision-making.

Altered tumor microenvironment heterogeneity of penile cancer during progression from non-lymphatic to lymphatic metastasis

BACKGROUND

Lymphatic metastasis is the major challenge in the treatment of penile cancer. The prognosis of individuals with lymphatic metastasis is extremely poor. Therefore, early identification of disease progression and lymphatic metastasis is an urgent task for researchers in penile cancer worldwide.

METHODS

In this study, using single-cell RNA sequencing, an immune landscape was established for the cancer ecosystem based on 46,861 cells from six patients with penile cancer (four with lymphatic metastasis [stage IV] and two without lymphatic metastasis [stage I]). Using bulk RNA sequencing, the discrepancy between the cancers and their respective metastatic lymph nodes was depicted based on seven patients with penile cancer.

RESULTS

The interaction between epithelial cells, fibroblasts, and endothelial cells, and the functional cooperation among invasion, epithelial-mesenchymal transition, and angiogenesis were found to be important landscapes in the penile cancer ecosystem, playing important roles in progression of cancer and lymph node metastasis.

CONCLUSIONS

This study is the first to investigate the altered tumor microenvironment heterogeneity of penile cancer as it evolves from non-lymphatic to lymphatic metastasis and provides insights into the mechanisms underlying malignant progression, the premetastatic niche, and lymphatic metastasis in penile cancer.

Machine learning for the identification of neoantigen-reactive CD8 + T cells in gastrointestinal cancer using single-cell sequencing

BACKGROUND

It appears that tumour-infiltrating neoantigen-reactive CD8 + T (Neo T) cells are the primary driver of immune responses to gastrointestinal cancer in patients. However, the conventional method is very time-consuming and complex for identifying Neo T cells and their corresponding T cell receptors (TCRs).

METHODS

By mapping neoantigen-reactive T cells from the single-cell transcriptomes of thousands of tumour-infiltrating lymphocytes, we developed a 26-gene machine learning model for the identification of neoantigen-reactive T cells.

RESULTS

In both training and validation sets, the model performed admirably. We discovered that the majority of Neo T cells exhibited notable differences in the biological processes of amide-related signal pathways. The analysis of potential cell-to-cell interactions, in conjunction with spatial transcriptomic and multiplex immunohistochemistry data, has revealed that Neo T cells possess potent signalling molecules, including LTA, which can potentially engage with tumour cells within the tumour microenvironment, thereby exerting anti-tumour effects. By sequencing CD8 + T cells in tumour samples of patients undergoing neoadjuvant immunotherapy, we determined that the fraction of Neo T cells was significantly and positively linked with the clinical benefit and overall survival rate of patients.

CONCLUSION

This method expedites the identification of neoantigen-reactive TCRs and the engineering of neoantigen-reactive T cells for therapy.

Progressive meristem and single-cell transcriptomes reveal the regulatory mechanisms underlying maize inflorescence development and sex differentiation

Maize develops separate ear and tassel inflorescences with initially similar morphology but ultimately different architecture and sexuality. The detailed regulatory mechanisms underlying these changes still remain largely unclear. In this study, through analyzing the time-course meristem transcriptomes and floret single-cell transcriptomes of ear and tassel, we revealed the regulatory dynamics and pathways underlying inflorescence development and sex differentiation. We identified 16 diverse gene clusters with differential spatiotemporal expression patterns and revealed biased regulation of redox, programmed cell death, and hormone signals during meristem differentiation between ear and tassel. Notably, based on their dynamic expression patterns, we revealed the roles of two RNA-binding proteins in regulating inflorescence meristem activity and axillary meristem formation. Moreover, using the transcriptional profiles of 53 910 single cells, we uncovered the cellular heterogeneity between ear and tassel florets. We found that multiple signals associated with either enhanced cell death or reduced growth are responsible for tassel pistil suppression, while part of the gibberellic acid signal may act non-cell-autonomously to regulate ear stamen arrest during sex differentiation. We further showed that the pistil-protection gene SILKLESS 1 (SK1) functions antagonistically to the known pistil-suppression genes through regulating common molecular pathways, and constructed a regulatory network for pistil-fate determination. Collectively, our study provides a deep understanding of the regulatory mechanisms underlying inflorescence development and sex differentiation in maize, laying the foundation for identifying new regulators and pathways for maize hybrid breeding and improvement.

Unraveling dynamic immunological landscapes in intracerebral hemorrhage: insights from single-cell and spatial transcriptomic profiling

Intracerebral hemorrhage (ICH) poses a formidable challenge in stroke management, with limited therapeutic options, particularly in the realm of immune-targeted interventions. Clinical trials targeting immune responses post-ICH have encountered setbacks, potentially attributable to the substantial cellular heterogeneity and intricate intercellular networks within the brain. Here, we present a pioneering investigation utilizing single-cell RNA sequencing and spatial transcriptome profiling at hyperacute (1 h), acute (24 h), and subacute (7 days) intervals post-ICH, aimed at unraveling the dynamic immunological landscape and spatial distributions within the cerebral tissue. Our comprehensive analysis revealed distinct cell differentiation patterns among myeloid and lymphocyte populations, along with delineated spatial distributions across various brain regions. Notably, we identified a subset of lymphocytes characterized by the expression of Spp1 and Lyz2, termed macrophage-associated lymphocytes, which exhibited close interactions with myeloid cells. Specifically, we observed prominent interactions between Lgmn+Macro-T cells and microglia through the spp1-cd44 pathway during the acute phase post-ICH in the choroid plexus. These findings represent a significant advancement in our understanding of immune cell dynamics at single-cell resolution across distinct post-ICH time points, thereby laying the groundwork for exploring critical temporal windows and informing the development of targeted therapeutic strategies.

Single-Cell RNA Sequencing Reveals Dysregulated POSTN+WNT5A+ Fibroblast Subclusters in Prurigo Nodularis

Prurigo nodularis (PN) is an intensely pruritic, inflammatory skin disease with a poorly understood pathogenesis. We performed single-cell transcriptomic profiling of 28,695 lesional and nonlesional PN cells. Lesional PN has increased dysregulated fibroblasts (FBs) and myofibroblasts. FBs in lesional PN were shifted toward a cancer-associated FB-like phenotype, with POSTN+WNT5A+ cancer-associated FBs increased in PN and similarly so in squamous cell carcinoma. A multicenter cohort study revealed an increased risk of squamous cell carcinoma and cancer-associated FB-associated malignancies (breast and colorectal) in patients with PN. Systemic fibroproliferative diseases (renal sclerosis and idiopathic pulmonary fibrosis) were upregulated in patients with PN. Ligand-receptor analyses demonstrated an FB neuronal axis with FB-derived WNT5A and periostin interactions with neuronal receptors melanoma cell adhesion molecule and ITGAV. These findings identify a pathogenic and targetable POSTN+WNT5A+ FB subpopulation that may predispose cancer-associated FB-associated malignancies in patients with PN.

PICALM Regulating the Generation of Amyloid β-Peptide to Promote Anthracycline-Induced Cardiotoxicity

Anthracyclines are chemotherapeutic drugs used to treat solid and hematologic malignancies. However, life-threatening cardiotoxicity, with cardiac dilation and heart failure, is a drawback. A combination of in vivo for single cell/nucleus RNA sequencing and in vitro approaches is used to elucidate the underlying mechanism. Genetic depletion and pharmacological blocking peptides on phosphatidylinositol binding clathrin assembly (PICALM) are used to evaluate the role of PICALM in doxorubicin-induced cardiotoxicity in vivo. Human heart tissue samples are used for verification. Patients with end-stage heart failure and chemotherapy-induced cardiotoxicity have thinner cell membranes compared to healthy controls do. Using the doxorubicin-induced cardiotoxicity mice model, it is possible to replicate the corresponding phenotype in patients. Cellular changes in doxorubicin-induced cardiotoxicity in mice, especially in cardiomyocytes, are identified using single cell/nucleus RNA sequencing. Picalm expression is upregulated only in cardiomyocytes with doxorubicin-induced cardiotoxicity. Amyloid β-peptide production is also increased after doxorubicin treatment, which leads to a greater increase in the membrane permeability of cardiomyocytes. Genetic depletion and pharmacological blocking peptides on Picalm reduce the generation of amyloid β-peptide. This alleviates the doxorubicin-induced cardiotoxicity in vitro and in vivo. In human heart tissue samples of patients with chemotherapy-induced cardiotoxicity, PICALM, and amyloid β-peptide are elevated as well.

Regulation of endocrine cell alternative splicing revealed by single-cell RNA sequencing in type 2 diabetes pathogenesis

The prevalent RNA alternative splicing (AS) contributes to molecular diversity, which has been demonstrated in cellular function regulation and disease pathogenesis. However, the contribution of AS in pancreatic islets during diabetes progression remains unclear. Here, we reanalyze the full-length single-cell RNA sequencing data from the deposited database to investigate AS regulation across human pancreatic endocrine cell types in non-diabetic (ND) and type 2 diabetic (T2D) individuals. Our analysis demonstrates the significant association between transcriptomic AS profiles and cell-type-specificity, which could be applied to distinguish the clustering of major endocrine cell types. Moreover, AS profiles are enabled to clearly define the mature subset of β-cells in healthy controls, which is completely lost in T2D. Further analysis reveals that RNA-binding proteins (RBPs), heterogeneous nuclear ribonucleoproteins (hnRNPs) and FXR1 family proteins are predicted to induce the functional impairment of β-cells through regulating AS profiles. Finally, trajectory analysis of endocrine cells suggests the β-cell identity shift through dedifferentiation and transdifferentiation of β-cells during the progression of T2D. Together, our study provides a mechanism for regulating β-cell functions and suggests the significant contribution of AS program during diabetes pathogenesis.

Transcriptomic diversity of innate lymphoid cells in human lymph nodes compared to BM and spleen

Innate lymphoid cells (ILCs) are largely tissue-resident, mostly described within the mucosal tissues. However, their presence and functions in the human draining lymph nodes (LNs) are unknown. Our study unravels the tissue-specific transcriptional profiles of 47,287 CD127+ ILCs within the human abdominal and thoracic LNs. LNs contain a higher frequency of CD127+ ILCs than in BM or spleen. We define independent stages of ILC development, including EILP and pILC in the BM. These progenitors exist in LNs in addition to naïve ILCs (nILCs) that can differentiate into mature ILCs. We define three ILC1 and four ILC3 sub-clusters in the LNs. ILC1 and ILC3 subsets have clusters with high heat shock protein-encoding genes. We identify previously unrecognized regulons, including the BACH2 family for ILC1 and the ATF family for ILC3. Our study is the comprehensive characterization of ILCs in LNs, providing an in-depth understanding of ILC-mediated immunity in humans.

Single-cell landscape of immunological responses in elderly patients with sepsis

Sepsis is a dysregulated host response to severe infections, and immune dysfunction plays a crucial role in its pathogenesis. Elderly patients, a special population influenced by immunosenescence, are more susceptible to sepsis and have a worse prognosis. However, the immunopathogenic mechanisms underlying sepsis in elderly patients remain unclear. Here, we performed single-cell RNA sequencing of peripheral blood samples from young and old subjects and patients with sepsis. By exploring the transcriptional profiles of immune cells, we analyzed immune cell compositions, phenotype shifts, expression heterogeneities, and intercellular communication. In elderly patients with sepsis, innate immune cells (e.g., monocytes and DCs) exhibit decreased antigen presentation, presenting an overactive inflammatory and senescent phenotype. However, the immunophenotype of T cells shifted to characterize effector, memory, and exhaustion. Moreover, we identified strong interferon-γ responses of T cells in both aging and sepsis groups and a deranged inflammaging status in elderly sepsis patients. Tregs in elderly patients with sepsis showed increased abundance and enhanced immunosuppressive effects. In addition, metabolism-associated pathways were upregulated in T cells in elderly patients with sepsis, and the lysine metabolism pathway was enriched in Tregs. Cell-cell interaction analysis showed that the expression profile of ligand-receptor pairs was probably associated with aggravated immune dysfunction in elderly patients with sepsis. A novel HLA-KIR interaction was observed between Tregs and CD8 + T cells. These findings illustrate the immunological hallmarks of sepsis in elderly patients, and highlight that immunosuppressive and metabolic regulatory pathways may undergo important alterations in elderly patients with sepsis.

Capturing totipotency in human cells through spliceosomal repression

The cleavage of zygotes generates totipotent blastomeres. In human 8-cell blastomeres, zygotic genome activation (ZGA) occurs to initiate the ontogenesis program. However, capturing and maintaining totipotency in human cells pose significant challenges. Here, we realize culturing human totipotent blastomere-like cells (hTBLCs). We find that splicing inhibition can transiently reprogram human pluripotent stem cells into ZGA-like cells (ZLCs), which subsequently transition into stable hTBLCs after long-term passaging. Distinct from reported 8-cell-like cells (8CLCs), both ZLCs and hTBLCs widely silence pluripotent genes. Interestingly, ZLCs activate a particular group of ZGA-specific genes, and hTBLCs are enriched with pre-ZGA-specific genes. During spontaneous differentiation, hTBLCs re-enter the intermediate ZLC stage and further generate epiblast (EPI)-, primitive endoderm (PrE)-, and trophectoderm (TE)-like lineages, effectively recapitulating human pre-implantation development. Possessing both embryonic and extraembryonic developmental potency, hTBLCs can autonomously generate blastocyst-like structures in vitro without external cell signaling. In summary, our study provides key criteria and insights into human cell totipotency.

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.

Tracing Immunological Interaction in Trimethylamine N-Oxide Hydrogel-Derived Zwitterionic Microenvironment During Promoted Diabetic Wound Regeneration

The diabetic wound healing is challenging due to the sabotaged delicate balance of immune regulation via an undetermined pathophysiological mechanism, so it is crucial to decipher multicellular signatures underlying diabetic wound healing and seek therapeutic strategies. Here, this work develops a strategy using novel trimethylamine N-oxide (TMAO)-derived zwitterionic hydrogel to promote diabetic wound healing, and explore the multi-cellular ecosystem around zwitterionic hydrogel, mapping out an overview of different cells in the zwitterionic microenvironment by single-cell RNA sequencing. The diverse cellular heterogeneity is revealed, highlighting the critical role of macrophage and neutrophils in managing diabetic wound healing. It is found that polyzwitterionic hydrogel can upregulate Ccl3+ macrophages and downregulate S100a9+ neutrophils and facilitate their interactions compared with polyanionic and polycationic hydrogels, validating the underlying effect of zwitterionic microenvironment on the activation of adaptive immune system. Moreover, zwitterionic hydrogel inhibits the formation of neutrophil extracellular traps (NETs) and promotes angiogenesis, thus improving diabetic wound healing. These findings expand the horizons of the sophisticated orchestration of immune systems in zwitterion-directed diabetic wound repair and uncover new strategies of novel immunoregulatory biomaterials.

Single-Cell Analysis Uncovers Striking Cellular Heterogeneity of Lung-Infiltrating Regulatory T Cells during Eosinophilic versus Neutrophilic Allergic Airway Inflammation

Allergic airway inflammation results from uncontrolled immune responses to environmental Ags. Although it is well established that allergic immune responses exhibit a high degree of diversity, driven by primary effector cell types such as eosinophils, neutrophils, or CD4 T cells with distinct effector signatures, the mechanisms responsible for such pathogenesis remain elusive. Foxp3+ regulatory T cells (Tregs) are essential immune regulators during chronic inflammation, including allergic airway inflammation. Emerging evidence suggests that Tregs infiltrating inflamed tissues exhibit distinct phenotypes dependent on the specific tissue sites and can display heterogeneity and tissue residency. Whether diverse allergic airway inflammatory responses influence infiltrating Treg heterogeneity or Treg lung residency has not been explored. We employed an unbiased single-cell RNA sequencing approach to investigate lung-infiltrating Tregs in models of eosinophilic and neutrophilic airway inflammation. We found that lung-infiltrating Tregs are highly heterogeneous, and that Tregs displaying lung-resident phenotypes are significantly different depending on the types of inflammation. Treg expression of ST2, a receptor for alarmin IL-33, was predominantly associated with eosinophilic inflammation and tissue residency. Nevertheless, Treg-specific ST2 deficiency did not affect the development of eosinophilic allergic inflammation or the generation of lung-resident Tregs. These results uncover a stark heterogeneity among Tregs infiltrating the lungs during allergic airway inflammation. The results indicate that varying types of inflammation may give rise to phenotypically distinct lung-resident Tregs, underscoring a (to our knowledge) novel mechanism by which inflammatory cues may shape the composition of infiltrating Tregs, allowing them to regulate inflammatory responses through tissue-adapted mechanisms.

Single-cell analysis identifies conserved features of immune dysfunction in simulated microgravity and spaceflight

Microgravity is associated with immunological dysfunction, though the mechanisms are poorly understood. Here, using single-cell analysis of human peripheral blood mononuclear cells (PBMCs) exposed to short term (25 hours) simulated microgravity, we characterize altered genes and pathways at basal and stimulated states with a Toll-like Receptor-7/8 agonist. We validate single-cell analysis by RNA sequencing and super-resolution microscopy, and against data from the Inspiration-4 (I4) mission, JAXA (Cell-Free Epigenome) mission, Twins study, and spleens from mice on the International Space Station. Overall, microgravity alters specific pathways for optimal immunity, including the cytoskeleton, interferon signaling, pyroptosis, temperature-shock, innate inflammation (e.g., Coronavirus pathogenesis pathway and IL-6 signaling), nuclear receptors, and sirtuin signaling. Microgravity directs monocyte inflammatory parameters, and impairs T cell and NK cell functionality. Using machine learning, we identify numerous compounds linking microgravity to immune cell transcription, and demonstrate that the flavonol, quercetin, can reverse most abnormal pathways. These results define immune cell alterations in microgravity, and provide opportunities for countermeasures to maintain normal immunity in space.

CD63 + tumor-associated macrophages drive the progression of hepatocellular carcinoma through the induction of epithelial-mesenchymal transition and lipid reprogramming

BACKGROUND

Tumor-associated macrophages (TAMs) constitute a substantial part of human hepatocellular carcinoma (HCC). The present study was devised to explore TAM diversity and their roles in HCC progression.

METHODS

Through the integration of multiple 10 × single-cell transcriptomic data derived from HCC samples and the use of consensus nonnegative matrix factorization (an unsupervised clustering algorithm), TAM molecular subtypes and expression programs were evaluated in detail. The roles played by these TAM subtypes in HCC were further probed through pseudotime, enrichment, and intercellular communication analyses. Lastly, vitro experiments were performed to validate the relationship between CD63, which is an inflammatory TAM expression program marker, and tumor cell lines.

RESULTS

We found that the inflammatory expression program in TAMs had a more obvious interaction with HCC cells, and CD63, as a marker gene of the inflammatory expression program, was associated with poor prognosis of HCC patients. Both bulk RNA-seq and vitro experiments confirmed that higher TAM CD63 expression was associated with the growth of HCC cells as well as their epithelial-mesenchymal transition, metastasis, invasion, and the reprogramming of lipid metabolism.

CONCLUSIONS

These analyses revealed that the TAM inflammatory expression program in HCC is closely associated with malignant tumor cells, with the hub gene CD63 thus representing an ideal target for therapeutic intervention in this cancer type.

A pan-cancer single-cell transcriptional analysis of antigen-presenting cancer-associated fibroblasts in the tumor microenvironment

Background

Cancer-associated fibroblasts (CAFs) are the primary stromal cells found in tumor microenvironment, and display high plasticity and heterogeneity. By using single-cell RNA-seq technology, researchers have identified various subpopulations of CAFs, particularly highlighting a recently identified subpopulation termed antigen-presenting CAFs (apCAFs), which are largely unknown.

Methods

We collected datasets from public databases for 9 different solid tumor types to analyze the role of apCAFs in the tumor microenvironment.

Results

Our data revealed that apCAFs, likely originating mainly from normal fibroblast, are commonly found in different solid tumor types and generally are associated with anti-tumor effects. apCAFs may be associated with the activation of CD4+ effector T cells and potentially promote the survival of CD4+ effector T cells through the expression of C1Q molecules. Moreover, apCAFs exhibited highly enrichment of transcription factors RUNX3 and IKZF1, along with increased glycolytic metabolism.

Conclusions

Taken together, these findings offer novel insights into a deeper understanding of apCAFs and the potential therapeutic implications for apCAFs targeted immunotherapy in cancer.

Remodeling ceramide homeostasis promotes functional maturation of human pluripotent stem cell-derived β cells

Human pluripotent stem cell-derived β cells (hPSC-β cells) show the potential to restore euglycemia. However, the immature functionality of hPSC-β cells has limited their efficacy in application. Here, by deciphering the continuous maturation process of hPSC-β cells post transplantation via single-cell RNA sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq), we show that functional maturation of hPSC-β cells is an orderly multistep process during which cells sequentially undergo metabolic adaption, removal of negative regulators of cell function, and establishment of a more specialized transcriptome and epigenome. Importantly, remodeling lipid metabolism, especially downregulating the metabolic activity of ceramides, the central hub of sphingolipid metabolism, is critical for β cell maturation. Limiting intracellular accumulation of ceramides in hPSC-β cells remarkably enhanced their function, as indicated by improvements in insulin processing and glucose-stimulated insulin secretion. In summary, our findings provide insights into the maturation of human pancreatic β cells and highlight the importance of ceramide homeostasis in function acquisition.

Single-cell analysis reveals hypoxia-induced immunosuppressive microenvironment in intrahepatic cholangiocarcinoma

The role of hypoxia in the tumor microenvironment of intrahepatic cholangiocarcinoma (iCCA) remains unclear. Here, we generated a comprehensive atlas of the entire tumor microenvironment and delineated the multifaceted cell-cell interactions to decipher hypoxia-induced pro-tumor immune suppression. We discovered hypoxia is significantly associated with iCCA progression via the activation of HIF1A expression. Moreover, hypoxia-dependent PPARγ-mediated fatty acid oxidation in APOE+ TAMs promoted M2 macrophage polarization by activating the HIF1A-PPARG-CD36 axis. These polarized APOE+ TAMs recruited Treg cell infiltration via the CCL3-CCR5 pair to form an immunosuppressive microenvironment. APOE+ TAMs tended to co-localize spatially with Treg cells in the malignant tissue based on spatial transcriptome data and immunofluorescence analysis results. We identified tumor-reactive CXCL13+ CD8-PreTex with specific high expression of ENTPD1 and ITGAE, which acted as precursors of CD8-Tex and had higher cytotoxicity, lower exhaustion, and more vigorous proliferation. Consequently, CXCL13+ CD8-PreTex functioned as a positive regulator of antitumor immunity by expressing the pro-inflammatory cytokines IFNG and TNF, associated with a better survival outcome. Our study reveals the mechanisms involved in hypoxia-induced immunosuppression and suggests that targeting precursor-exhausted CXCL13+CD8+ T cells might provide a pratical immunotherapeutic approach.

Multi-omics Analysis Revealed that the CCN Family Regulates Cell Crosstalk, Extracellular Matrix, and Immune Escape, Leading to a Poor Prognosis of Glioma

The CCN family is a group of matricellular proteins associated with the extracellular matrix. This study aims to explore the role of the CCN family in glioma development and its implications in the tumor microenvironment. Through analysis of bulk RNA-seq cohorts, correlations between CCN family expression and glioma subtypes, patient survival, and bioactive pathway enrichment were investigated. Additionally, single-cell datasets were employed to identify novel cell subgroups, followed by analyses of cell communication and transcription factors. Spatial transcriptomic analysis was utilized to validate the CCN family's involvement in glioma. Results indicate overexpression of CYR61,CTGF, and WISP1 in glioma, associated with unfavorable subtypes and reduced survival. Enrichment analyses revealed associations with oncogenic pathways, while CTGF and WISP1 expression correlated with increased infiltration of regulatory T cells and M2 macrophages. Single-cell analysis identified MES-like cells as the highest CCN expression. Moreover, intercellular signal transduction analysis demonstrated active pathways, including SPP1-CD44, in cell subgroups with elevated CYR61 and CTGF expression. Spatial transcriptomic analysis confirmed co-localization of CYR61,CTGF and SPP1-CD44 with high oncogenic pathway activity. These findings suggest that CCN family members may serve as potential prognostic biomarkers and therapeutic targets for glioma.

Neuronal cell populations in circumoral nerve ring of sea cucumber Apostichopus japonicus: Ultrastructure and transcriptional profile

The echinoderm nervous system has been studied as a model for understanding the evolution of the chordate nervous system. Neuronal cells are essential groups that release a 'cocktail' of messenger molecules providing a spectrum of biological actions in the nervous system. Among echinoderms, most evidence on neuronal cell types has been obtained from starfish and sea urchin. In sea cucumbers, most research has focused on the location of neuronal cells, whereas their transcriptional features have rarely been investigated. Here, we observed the ultrastructure of neuronal cells in the sea cucumber, Apostichopus japonicus. The transcriptional profile of neuronal cells from the circumoral nerve ring (CNR) was investigated using single-cell RNA sequencing (scRNA-seq), and a total of six neuronal cell types were identified. 26 neuropeptide precursor genes (NPPs) and 28 G-protein-coupled receptors (GPCR) were expressed in the six neuronal cell types, comprising five NPP/NP-GPCR pairs. Unsupervised pseudotime analysis of neuronal cells showed their different differentiation status. We also located the neuronal cells in the CNR by immunofluorescence (IF) and identified the potential hub genes of key cell populations. This broad resource serves as a valuable support in the development of cell-specific markers for accurate cell-type identification in sea cucumbers. It also contributes to facilitating comparison across species, providing a deeper understanding of the evolutionary processes of neuronal cells.

Cell communication pathway prognostic model identified detrimental neurodevelopmental pathways in neuroblastoma

Neurodevelopmental cell communication plays a crucial role in neuroblastoma prognosis. However, determining the impact of these communication pathways on prognosis is challenging due to limited sample sizes and patchy clinical survival information of single cell RNA-seq data. To address this, we have developed the cell communication pathway prognostic model (CCPPM) in this study. CCPPM involves the identification of communication pathways through single-cell RNA-seq data, screening of prognosis-significant pathways using bulk RNA-seq data, conducting functional and attribute analysis of these pathways, and analyzing the post-effects of communication within these pathways. By employing the CCPPM, we have identified ten communication pathways significantly influencing neuroblastoma, all related to axongenesis and neural projection development, especially the BMP7-(BMPR1B-ACVR2B) communication pathway was found to promote tumor cell migration by activating the transcription factor SMAD1 and regulating UNK and MYCBP2. Notably, BMP7 expression was higher in neuroblastoma samples with distant metastases. In summary, CCPPM offers a novel approach to studying the influence of cell communication pathways on disease prognosis and identified detrimental communication pathways related to neurodevelopment.

Unravelling infiltrating T-cell heterogeneity in kidney renal clear cell carcinoma: Integrative single-cell and spatial transcriptomic profiling

Kidney renal clear cell carcinoma (KIRC) pathogenesis intricately involves immune system dynamics, particularly the role of T cells within the tumour microenvironment. Through a multifaceted approach encompassing single-cell RNA sequencing, spatial transcriptome analysis and bulk transcriptome profiling, we systematically explored the contribution of infiltrating T cells to KIRC heterogeneity. Employing high-density weighted gene co-expression network analysis (hdWGCNA), module scoring and machine learning, we identified a distinct signature of infiltrating T cell-associated genes (ITSGs). Spatial transcriptomic data were analysed using robust cell type decomposition (RCTD) to uncover spatial interactions. Further analyses included enrichment assessments, immune infiltration evaluations and drug susceptibility predictions. Experimental validation involved PCR experiments, CCK-8 assays, plate cloning assays, wound-healing assays and Transwell assays. Six subpopulations of infiltrating and proliferating T cells were identified in KIRC, with notable dynamics observed in mid- to late-stage disease progression. Spatial analysis revealed significant correlations between T cells and epithelial cells across varying distances within the tumour microenvironment. The ITSG-based prognostic model demonstrated robust predictive capabilities, implicating these genes in immune modulation and metabolic pathways and offering prognostic insights into drug sensitivity for 12 KIRC treatment agents. Experimental validation underscored the functional relevance of PPIB in KIRC cell proliferation, invasion and migration. Our study comprehensively characterizes infiltrating T-cell heterogeneity in KIRC using single-cell RNA sequencing and spatial transcriptome data. The stable prognostic model based on ITSGs unveils infiltrating T cells' prognostic potential, shedding light on the immune microenvironment and offering avenues for personalized treatment and immunotherapy.

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.

Single-cell transcriptome atlas reveals somatic cell embryogenic differentiation features during regeneration

Understanding somatic cell totipotency remains a challenge facing scientific inquiry today. Plants display remarkable cell totipotency expression, illustrated by single-cell differentiation during somatic embryogenesis (SE) for plant regeneration. Determining cell identity and exploring gene regulation in such complex heterogeneous somatic cell differentiation have been major challenges. Here, we performed high-throughput single-cell sequencing assays to define the precise cellular landscape and revealed the modulation mode of marker genes during embryogenic differentiation in cotton (Gossypium hirsutum L.) as the crop for biotechnology application. We demonstrated that nonembryogenic calli (NEC) and primary embryogenic calli (PEC) tissues were composed of heterogeneous cells that could be partitioned into four broad populations with six distinct cell clusters. Enriched cell clusters and cell states were identified in NEC and PEC samples, respectively. Moreover, a broad repertoire of new cluster-specific genes and associated expression modules were identified. The energy metabolism, signal transduction, environmental adaptation, membrane transport pathways, and a series of transcription factors were preferentially enriched in cell embryogenic totipotency expression. Notably, the SE-ASSOCIATED LIPID TRANSFER PROTEIN (SELTP) gene dose-dependently marked cell types with distinct embryogenic states and exhibited a parabolic curve pattern along the somatic cell embryogenic differentiation trajectory, suggesting that SELTP could serve as a favorable quantitative cellular marker for detecting embryogenic expression at the single-cell level. In addition, RNA velocity and Scissor analysis confirmed the pseudo-temporal model and validated the accuracy of the scRNA-seq data, respectively. This work provides valuable marker-genes resources and defines precise cellular taxonomy and trajectory atlases for somatic cell embryogenic differentiation in plant regeneration.

A roadmap to human hippocampal neurogenesis in adulthood, aging and AD

In the rodent, hippocampal neurogenesis plays critical roles in learning and memory1,2, is tightly regulated by inhibitory neurons3-7 and contributes to memory dysfunction in Alzheimer's disease (AD) mouse models8-10. In contrast, the mechanisms regulating neurogenesis in the adult human hippocampus, the dynamic shifts in the transcriptomic and epigenomic profiles in aging and AD and putative niche interactions within the cellular environment, remain largely unknown. Using single nuclei multi-omics of postmortem human hippocampi we map the molecular mechanisms of hippocampal neurogenesis across aging, cognitive decline, and AD neuropathology. Transcriptomic and epigenetic profiling of neural stem cells (NSCs), neuroblasts and immature neurons suggests that the earliest shift in the characteristics of neurogenesis takes place in NSCs in aging. Cognitive impairment was associated with changes in neuroblast profile. In AD, there was a widespread cessation of the transcription machinery in immature neurons, with robust downregulation of genes regulating ribosomal and mitochondrial function. Further, there was substantial loss of parvalbumin+ inhibitory neurons in the hippocampus in aging. The number of the rest of inhibitory neurons were reduced as a function of age and diagnosis. Notably, a similar system-level effect was observed between immature and inhibitory neurons in the transition from aging to AD, manifested by common molecular pathways that were ultimately lost in AD. The numbers of neuroblasts, immature and GABAergic neurons inversely correlated with extent of neuropathology. Using CellChat and NeuronChat, we inferred the ligands and receptors by which neurogenic cells communicate with their cellular environment. Loss of synaptic adhesion molecules and neurotransmitters, either sent or received by neurogenic cells, was observed in AD. Together, this study delineates the molecular mechanisms and dynamics of human neurogenesis, functional association with inhibitory neurons and a mechanism of hippocampal hyperexcitability in AD.

Single-cell RNA sequencing reveals the MIF/ACKR3 receptor-ligand interaction between neutrophils and nucleus pulposus cells in intervertebral disc degeneration

OBJECTIVES

To unravel the heterogeneity and function of microenvironmental neutrophils during intervertebral disc degeneration (IDD).

METHODS

Single-cell RNA sequencing (scRNA-seq) was utilized to dissect the cellular landscape of neutrophils in intervertebral disc (IVD) tissues and their crosstalk with nucleus pulposus cells (NPCs). The expression levels of macrophage migration inhibitory factor (MIF) and ACKR3 in IVD tissues were detected. The MIF/ACKR3 axis was identified and its effects on IDD were investigated in vitro and in vivo.

RESULTS

We sequenced here 71520 single cells from 5 control and 9 degenerated IVD samples using scRNA-seq. We identified a unique cluster of neutrophils abundant in degenerated IVD tissues that highly expressed MIF and was functionally enriched in extracellular matrix organization (ECMO). Cell-to-cell communication analyses showed that this ECMO-neutrophil subpopulation was closely interacted with an effector NPCs subtype, which displayed high expression of ACKR3. Further analyses revealed that MIF was positively correlated with ACKR3 and functioned via directly binding to ACKR3 on effector NPCs. MIF inhibition attenuated degenerative changes of NPCs and extracellular matrix, which could be partially reversed by ACKR3 overexpression. Clinically, a significant correlation of high MIF/ACKR3 expression with advanced IDD grade was observed. Furthermore, we also found a positive association between MIF+ ECMO-neutrophil counts and ACKR3+ effector NPCs density as well as higher expression of the MIF/ACKR3 signaling in areas where these two cell types were neighbors.

CONCLUSIONS

These data suggest that ECMO-neutrophil promotes IDD progression by their communication with NPCs via the MIF/ACKR3 axis, which may shed light on therapeutic strategies.

Ly6Chi Monocytes Are Metabolically Reprogrammed in the Blood during Inflammatory Stimulation and Require Intact OxPhos for Chemotaxis and Monocyte to Macrophage Differentiation

Acute inflammation is a rapid and dynamic process involving the recruitment and activation of multiple cell types in a coordinated and precise manner. Here, we investigate the origin and transcriptional reprogramming of monocytes using a model of acute inflammation, zymosan-induced peritonitis. Monocyte trafficking and adoptive transfer experiments confirmed that monocytes undergo rapid phenotypic change as they exit the blood and give rise to monocyte-derived macrophages that persist during the resolution of inflammation. Single-cell transcriptomics revealed significant heterogeneity within the surface marker-defined CD11b+Ly6G-Ly6Chi monocyte populations within the blood and at the site of inflammation. We show that two major transcriptional reprogramming events occur during the initial six hours of Ly6Chi monocyte mobilisation, one in the blood priming monocytes for migration and a second at the site of inflammation. Pathway analysis revealed an important role for oxidative phosphorylation (OxPhos) during both these reprogramming events. Experimentally, we demonstrate that OxPhos via the intact mitochondrial electron transport chain is essential for murine and human monocyte chemotaxis. Moreover, OxPhos is needed for monocyte-to-macrophage differentiation and macrophage M(IL-4) polarisation. These new findings from transcriptional profiling open up the possibility that shifting monocyte metabolic capacity towards OxPhos could facilitate enhanced macrophage M2-like polarisation to aid inflammation resolution and tissue repair.