Single-cell mRNA quantification and differential analysis with Census

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.

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.

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.

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.

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.

Single-Nucleus RNA-Seq Reveals Spermatogonial Stem Cell Developmental Pattern in Shaziling Pigs

Normal testicular development ensures the process of spermatogenesis, which is a complex biological process. The sustained high productivity of spermatogenesis throughout life is predominantly attributable to the constant proliferation and differentiation of spermatogonial stem cells (SSCs). The self-renewal and differentiation processes of SSCs are strictly regulated by the SSC niche. Therefore, understanding the developmental pattern of SSCs is crucial for spermatogenesis. The Shaziling pig is a medium-sized indigenous pig breed originating from central China. It is renowned for its superior meat quality and early male sexual maturity. The spermatogenic ability of the boars is of great economic importance to the pig industry. To investigate testicular development, particularly the pattern of SSC development in Shaziling pigs, we used single-cell transcriptomics to identify gene expression patterns in 82,027 individual cells from nine Shaziling pig testes at three key postnatal developmental stages. We generated an unbiased cell developmental atlas of Shaziling pig testicular tissues. We elucidated the complex processes involved in the development of SSCs within their niche in the Shaziling pig. Specifically, we identified potential marker genes and cellular signaling pathways that regulate SSC self-renewal and maintenance. Additionally, we proposed potential novel marker genes for SSCs that could be used for SSC isolation and sorting in Shaziling pigs. Furthermore, by immunofluorescence staining of testicular tissues of different developmental ages using marker proteins (UCHL1 and KIT), the developmental pattern of the spermatogonia of Shaziling pigs was intensively studied. Our research enhances the comprehension of the development of SSCs and provides a valuable reference for breeding Shaziling pigs.

Overactivated MX1 Positive Natural Killer Cells Promote the Progression of Sepsis-Induced Acute Respiratory Distress Syndrome

Background

Natural killer (NK) cells are key regulators of immune defense in sepsis-induced acute respiratory distress syndrome (ARDS), yet the characteristics of NK cell clusters in ARDS remain poorly understood.

Methods

A prospective and observational study enrolled septic patients with ARDS or not was conducted to determine the percentage of NK cells via flow cytometry. The transcriptomes of peripheral blood mononuclear cells (PBMCs) from healthy controls, patients with sepsis only, and patients with sepsis-induced ARDS were profiled. Vitro experiments were performed to confirm the mechanism mediating MX1+NK cell infiltration.

Results

A total of 115 septic patients were analyzed, among whom 63 patients developed ARDS and 52 patients did not. Decreased NK percentages were found in sepsis with ARDS patients (%, 7.46±4.40 vs 11.65±6.88, P=0.0001) compared with sepsis-only patients. A lower percentage of NK cells showed a significant increase in 28-day mortality. Single-cell sequencing analysis revealed distinct characteristics of NK cells in sepsis-induced ARDS, notably the identification of a unique cluster defined as MX1+NK cells. Flow cytometry analysis showed an elevated percentage of MX1+NK cells specifically in individuals with sepsis-induced ARDS, compared with patients with sepsis only. Pseudo-time analysis showed that MX1+NK cells were characterized by upregulation of type I interferon-induced genes and other pro-inflammatory genes. MX1+NK cells can respond to type I interferons and secrete type I interferons themselves. Ligand-receptor interaction analysis also revealed extensive interaction between MX1+NK cells and T/B cells, leading to an uncontrolled inflammatory response in ARDS.

Conclusion

MX1+NK cells can respond to type I interferons and secrete type I interferons themselves, promoting the development of sepsis-induced ARDS. Interfering with the infiltration of MX1+NK cells could be a therapeutic approach for this disease. Due to the limited sample size, a larger sample size was needed for further exploration.

Macrophages communicate with mesangial cells through the CXCL12/DPP4 axis in lupus nephritis pathogenesis

Lupus nephritis (LN) occurs in 50% of cases of systemic lupus erythematosus (SLE) and is one of the most serious complications that can occur during lupus progression. Mesangial cells (MCs) are intrinsic cells in the kidney that can regulate capillary blood flow, phagocytose apoptotic cells, and secrete vasoactive substances and growth factors. Previous studies have shown that various types of inflammatory cells can activate MCs for hyperproliferation, leading to disruption of the filtration barrier and impairment of renal function in LN. Here, we characterized the heterogeneity of kidney cells of LN mice by single-nucleus RNA sequencing (snRNA-seq) and revealed the interaction between macrophages and MCs through the CXC motif chemokine ligand 12 (CXCL12)/dipeptidyl peptidase 4 (DPP4) axis. In culture, macrophages modulated the proliferation and migration of MCs through this ligand-receptor interaction. In LN mice, treatment with linagliptin, a DPP4 inhibitor, effectively inhibited MC proliferation and reduced urinary protein levels. Together, our findings indicated that targeting the CXCL12/DPP4 axis with linagliptin treatment may serve as a novel strategy for the treatment of LN via the CXCL12/DPP4 axis.

Long noncoding RNAs heat shock RNA omega nucleates TBPH and promotes intestinal stem cell differentiation upon heat shock

In Drosophila, long noncoding RNA Hsrω rapidly assembles membraneless organelle omega speckles under heat shock with unknown biological function. Here, we identified the distribution of omega speckles in multiple tissues of adult Drosophila melanogaster and found that they were selectively distributed in differentiated enterocytes but not in the intestinal stem cells of the midgut. We mimicked the high expression level of Hsrω via overexpression or intense heat shock and demonstrated that the assembly of omega speckles nucleates TBPH for the induction of ISC differentiation. Additionally, we found that heat shock stress promoted cell differentiation, which is conserved in mammalian cells through paraspeckles, resulting in large puncta of TDP-43 (a homolog of TBPH) with less mobility and the differentiation of human induced pluripotent stem cells. Overall, our findings confirm the role of Hsrω and omega speckles in the development of intestinal cells and provide new prospects for the establishment of stem cell differentiation strategies.

Mendelian Randomization and Transcriptomic Analysis Reveal the Protective Role of NKT Cells in Sepsis

Background

Sepsis is a life-threatening clinical syndrome caused by dysregulated host response to infection. The mechanism underlying sepsis-induced immune dysfunction remains poorly understood. Natural killer T (NKT) cells are cytotoxic lymphocytes that bridge the innate and adaptive immune systems, the role of NKT cells in sepsis is not entirely understood, and NKT cell cluster differences in sepsis remain unexplored.

Methods

Mendelian randomization (MR) analyses were first conducted to investigate the causal relationship between side scatter area (SSC-A) on NKT cells and 28-day mortality of septic patients. A prospective and observational study was conducted to validate the relationship between the percentage of NKT cells and 28-day mortality of sepsis. Then, the single-cell RNA sequencing (scRNA-seq) data of peripheral blood mononuclear cells (PBMCs) from healthy controls and septic patients were profiled.

Results

MR analyses first revealed the protective roles of NKT cells in the 28-day mortality of sepsis. Then, 115 septic patients were enrolled. NKT percentage was significantly higher in survivors (n = 84) compared to non-survivors (n = 31) (%, 5.00 ± 3.46 vs 2.18 ± 1.93, P < 0.0001). Patients with lower levels of NKT cells exhibited a significantly increased risk of 28-day mortality. According to scRNA-seq analysis, NKT cell clusters exhibited multiple distinctive characteristics, including a distinguishing cluster defined as FOS+NKT cells, which showed a significant decrease in sepsis. Pseudo-time analysis showed that FOS+NKT cells were characterized by upregulated expression of crucial functional genes such as GZMA and CCL4. CellChat revealed that interactions between FOS+NKT cells and adaptive immune cells including B cells and T cells were decreased in sepsis compared to healthy controls.

Conclusion

Our findings indicate that NKT cells may protect against sepsis, and their percentage can predict 28-day mortality. Additionally, we discovered a unique FOS+NKT subtype crucial in sepsis immune response, offering novel insights into its immunopathogenesis.

IL-4 Licenses B Cell Activation Through Cholesterol Synthesis

Lymphocyte activation involves a transition from quiescence and associated catabolic metabolism to a metabolic state with noted similarities to cancer cells such as heavy reliance on aerobic glycolysis for energy demands and increased nutrient requirements for biomass accumulation and cell division 1-3 . Following antigen receptor ligation, lymphocytes require spatiotemporally distinct "second signals". These include costimulatory receptor or cytokine signaling, which engage discrete programs that often involve remodeling of organelles and increased nutrient uptake or synthesis to meet changing biochemical demands 4-6 . One such signaling molecule, IL-4, is a highly pleiotropic cytokine that was first identified as a B cell co-mitogen over 30 years ago 7 . However, how IL-4 signaling mechanistically supports B cell proliferation is incompletely understood. Here, using single cell RNA sequencing we find that the cholesterol biosynthetic program is transcriptionally upregulated following IL-4 signaling during the early B cell response to influenza virus infection, and is required for B cell activation in vivo . By limiting lipid availability in vitro , we determine cholesterol to be essential for B cells to expand their endoplasmic reticulum, progress through cell cycle, and proliferate. In sum, we demonstrate that the well-known ability of IL-4 to act as a B cell growth factor is through a previously unknown rewiring of specific lipid anabolic programs, relieving sensitivity of cells to environmental nutrient availability.

Two sequential gene expression programs bridged by cell division support long-distance collective cell migration

The precise assembly of tissues and organs relies on spatiotemporal regulation of gene expression to coordinate the collective behavior of cells. In Drosophila embryos, the midgut musculature is formed through collective migration of caudal visceral mesoderm (CVM) cells, but how gene expression changes as cells migrate is not well understood. Here, we have focused on ten genes expressed in the CVM and the cis-regulatory sequences controlling their expression. Although some genes are continuously expressed, others are expressed only early or late during migration. Late expression relates to cell cycle progression, as driving string/Cdc25 causes earlier division of CVM cells and accelerates the transition to late gene expression. In particular, we found that the cell cycle effector transcription factor E2F1 is a required input for the late gene CG5080. Furthermore, whereas late genes are broadly expressed in all CVM cells, early gene transcripts are polarized to the anterior or posterior ends of the migrating collective. We show this polarization requires transcription factors Snail, Zfh1 and Dorsocross. Collectively, these results identify two sequential gene expression programs bridged by cell division that support long-distance directional migration of CVM cells.

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 promote robust meiotic gene expression in spermatogenic cells that have initiated meiosis. Here, we 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 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 transcriptional regulator STRA8-MEIOSIN, required for the meiotic G1/S phase transition and meiotic gene expression. We conclude that in mice, mRNA decay mediated by MEIOC-YTHDC2-RBM46 enhances the competence of spermatogenic cells to initiate meiosis.

SUMMARY STATEMENT

RNA-binding complex MEIOC-YTHDC2-RBM46 destabilizes its mRNA targets, including transcriptional repressors. This activity facilitates the retinoic acid-dependent activation of Meiosin gene expression and transition into meiosis.

Data normalization for addressing the challenges in the analysis of single-cell transcriptomic datasets

BACKGROUND

Normalization is a critical step in the analysis of single-cell RNA-sequencing (scRNA-seq) datasets. Its main goal is to make gene counts comparable within and between cells. To do so, normalization methods must account for technical and biological variability. Numerous normalization methods have been developed addressing different sources of dispersion and making specific assumptions about the count data.

MAIN BODY

The selection of a normalization method has a direct impact on downstream analysis, for example differential gene expression and cluster identification. Thus, the objective of this review is to guide the reader in making an informed decision on the most appropriate normalization method to use. To this aim, we first give an overview of the different single cell sequencing platforms and methods commonly used including isolation and library preparation protocols. Next, we discuss the inherent sources of variability of scRNA-seq datasets. We describe the categories of normalization methods and include examples of each. We also delineate imputation and batch-effect correction methods. Furthermore, we describe data-driven metrics commonly used to evaluate the performance of normalization methods. We also discuss common scRNA-seq methods and toolkits used for integrated data analysis.

CONCLUSIONS

According to the correction performed, normalization methods can be broadly classified as within and between-sample algorithms. Moreover, with respect to the mathematical model used, normalization methods can further be classified into: global scaling methods, generalized linear models, mixed methods, and machine learning-based methods. Each of these methods depict pros and cons and make different statistical assumptions. However, there is no better performing normalization method. Instead, metrics such as silhouette width, K-nearest neighbor batch-effect test, or Highly Variable Genes are recommended to assess the performance of normalization methods.

Deciphering Müller cell heterogeneity signatures in diabetic retinopathy across species: an integrative single-cell analysis

Diabetic retinopathy (DR), a leading cause of visual impairment, demands a profound comprehension of its cellular mechanisms to formulate effective therapeutic strategies. Our study presentes a comprehensive single-cell analysis elucidating the intricate landscape of Müller cells within DR, emphasizing their nuanced involvement. Utilizing scRNA-seq data from both Sprague-Dawley rat models and human patients, we delineated distinct Müller cell clusters and their corresponding gene expression profiles. These findings were further validated through differential gene expression analysis utilizing human transcriptomic data. Notably, certain Müller cell clusters displayed upregulation of the Rho gene, implying a phagocytic response to damaged photoreceptors within the DR microenvironment. This phenomenon was consistently observed across species. Additionally, the co-expression patterns of RHO and PDE6G within Müller cell clusters provided compelling evidence supporting their potential role in maintaining retinal integrity during DR. Our results offer novel insights into the cellular dynamics of DR and underscore Müller cells as promising therapeutic targets for preserving vision in retinal disorders induced by diabetes.

Landscape of gut mucosal immune cells showed gap of follicular or memory B cells into plasma cells in immunological non-responders

BACKGROUND

The gut is an important site for human immunodeficiency virus (HIV) infection and immune responses. The role of gut mucosal immune cells in immune restoration in patients infected with HIV undergoing antiretroviral therapy remains unclear.

METHODS

Ileocytes, including 54 475 immune cells, were obtained from colonoscopic biopsies of five HIV-negative controls, nine immunological responders (IRs), and three immunological non-responders (INRs) and were analyzed using single-cell RNA sequencing. Immunohistochemical assays were performed for validation. The 16S rRNA gene was amplified using PCR in faecal samples to analyze faecal microbiota. Flow cytometry was used to analyze CD4+ T-cell counts and the activation of T cells.

RESULTS

This study presents a global transcriptomic profile of the gut mucosal immune cells in patients infected with HIV. Compared with the IRs, the INRs exhibited a lower proportion of gut plasma cells, especially the IGKC+IgA+ plasma cell subpopulation. IGKC+IgA+ plasma cells were negatively associated with enriched f. Prevotellaceae the INRs and negatively correlated with the overactivation of T cells, but they were positively correlated with CD4+ T-cell counts. The INRs exhibited a higher proportion of B cells than the IRs. Follicular and memory B cells were significantly higher in the INRs. Reduced potential was observed in the differentiation of follicular or memory B cells into gut plasma cells in INRs. In addition, the receptor-ligand pairs CD74_MIF and CD74_COPA of memory B/ follicular helper T cells were significantly reduced in the INRs, which may hinder the differentiation of memory and follicular B cells into plasma cells.

CONCLUSIONS

Our study shows that plasma cells are dysregulated in INRs and provides an extensive resource for deciphering the immune pathogenesis of HIV in INRs.

KEY POINTS

An investigation was carried out at the single-cell-level to analyze gut mucosal immune cells alterations in PLWH after ART. B cells were significantly increased and plasma cells were significantly decreased in the INRs compared to the IRs and NCs. There are gaps in the transition from gut follicular or memory B cellsinto plasma cells in INRs.

Deciphering Tumor Cell Evolution in Cutaneous T-Cell Lymphomas: Distinct Differentiation Trajectories in Mycosis Fungoides and Sézary Syndrome

Cutaneous T-cell lymphomas are a heterogeneous group of neoplasms originating in the skin, with mycosis fungoides (MF) and Sézary syndrome (SS) representing the most common variants. The cellular origin of cutaneous lymphomas has remained controversial owing to their immense phenotypic heterogeneity that obfuscates lineage reconstruction on the basis of classical surface biomarkers. To overcome this heterogeneity and reconstruct the differentiation trajectory of malignant cells in MF and SS, TCR sequencing was performed in parallel with targeted transcriptomics at the single-cell resolution among cutaneous samples in MF and SS. Unsupervised lineage reconstruction showed that Sézary cells exist as a population of CD4+ T cells distinct from those in patch, plaque, and tumor MF. Further investigation of malignant cell heterogeneity in SS showed that Sézary cells phenotypically comprised at least 3 subsets on the basis of differential proliferation potentials and expression of exhaustion markers. A T helper 1-polarized cell type, intermediate cell type, and exhausted T helper 2-polarized cell type were identified, with T helper 1- and T helper 2-polarized cells displaying divergent proliferation potentials. Collectively, these findings provide evidence to clarify the relationship between MF and SS and reveal cell subsets in SS that suggest a possible mechanism for therapeutic resistance.

Revealing the pathogenesis of keloids based on the status: Active vs inactive

Recently, the pathomechanisms of keloids have been extensively researched using transcriptomic analysis, but most studies did not consider the activity of keloids. We aimed to profile the transcriptomics of keloids according to their clinical activity and location within the keloid lesion, compared with normal and mature scars. Tissue samples were collected (keloid based on its activity (active and inactive), mature scar from keloid patients and normal scar (NS) from non-keloid patients). To reduce possible bias, all keloids assessed in this study had no treatment history and their location was limited to the upper chest or back. Multiomics assessment was performed by using single-cell RNA sequencing and multiplex immunofluorescence. Increased mesenchymal fibroblasts (FBs) was the main feature in keloid patients. Noticeably, the proportion of pro-inflammatory FBs was significantly increased in active keloids compared to inactive ones. To explore the nature of proinflammatory FBs, trajectory analysis was conducted and CCN family associated with mechanical stretch exhibited higher expression in active keloids. For vascular endothelial cells (VECs), the proportion of tip and immature cells increased in keloids compared to NS, especially at the periphery of active keloids. Also, keloid VECs highly expressed genes with characteristics of mesenchymal activation compared to NS, especially those from the active keloid center. Multiomics analysis demonstrated the distinct expression profile of active keloids. Clinically, these findings may provide the future appropriate directions for development of treatment modalities of keloids. Prevention of keloids could be possible by the suppression of mesenchymal activation between FBs and VECs and modulation of proinflammatory FBs may be the key to the control of active keloids.

Age-specific nasal epithelial responses to SARS-CoV-2 infection

Children infected with SARS-CoV-2 rarely progress to respiratory failure. However, the risk of mortality in infected people over 85 years of age remains high. Here we investigate differences in the cellular landscape and function of paediatric (<12 years), adult (30-50 years) and older adult (>70 years) ex vivo cultured nasal epithelial cells in response to infection with SARS-CoV-2. We show that cell tropism of SARS-CoV-2, and expression of ACE2 and TMPRSS2 in nasal epithelial cell subtypes, differ between age groups. While ciliated cells are viral replication centres across all age groups, a distinct goblet inflammatory subtype emerges in infected paediatric cultures and shows high expression of interferon-stimulated genes and incomplete viral replication. In contrast, older adult cultures infected with SARS-CoV-2 show a proportional increase in basaloid-like cells, which facilitate viral spread and are associated with altered epithelial repair pathways. We confirm age-specific induction of these cell types by integrating data from in vivo COVID-19 studies and validate that our in vitro model recapitulates early epithelial responses to SARS-CoV-2 infection.

Single-cell transcriptomics: background, technologies, applications, and challenges

Single-cell sequencing was developed as a high-throughput tool to elucidate unusual and transient cell states that are barely visible in the bulk. This technology reveals the evolutionary status of cells and differences between populations, helps to identify unique cell subtypes and states, reveals regulatory relationships between genes, targets and molecular mechanisms in disease processes, tumor heterogeneity, the state of the immune environment, etc. However, the high cost and technical limitations of single-cell sequencing initially prevented its widespread application, but with advances in research, numerous new single-cell sequencing techniques have been discovered, lowering the cost barrier. Many single-cell sequencing platforms and bioinformatics methods have recently become commercially available, allowing researchers to make fascinating observations. They are now increasingly being used in various industries. Several protocols have been discovered in this context and each technique has unique characteristics, capabilities and challenges. This review presents the latest advancements in single-cell transcriptomics technologies. This includes single-cell transcriptomics approaches, workflows and statistical approaches to data processing, as well as the potential advances, applications, opportunities and challenges of single-cell transcriptomics technology. You will also get an overview of the entry points for spatial transcriptomics and multi-omics.

Identifying tumor type and cell type-specific gene expression alterations in pediatric central nervous system tumors

Central nervous system (CNS) tumors are the leading cause of pediatric cancer death, and these patients have an increased risk for developing secondary neoplasms. Due to the low prevalence of pediatric CNS tumors, major advances in targeted therapies have been lagging compared to other adult tumors. We collect single nuclei RNA-seq data from 84,700 nuclei of 35 pediatric CNS tumors and three non-tumoral pediatric brain tissues and characterize tumor heterogeneity and transcriptomic alterations. We distinguish cell subpopulations associated with specific tumor types including radial glial cells in ependymomas and oligodendrocyte precursor cells in astrocytomas. In tumors, we observe pathways important in neural stem cell-like populations, a cell type previously associated with therapy resistance. Lastly, we identify transcriptomic alterations among pediatric CNS tumor types compared to non-tumor tissues, while accounting for cell type effects on gene expression. Our results suggest potential tumor type and cell type-specific targets for pediatric CNS tumor treatment. Here we address current gaps in understanding single nuclei gene expression profiles of previously under-investigated tumor types and enhance current knowledge of gene expression profiles of single cells of various pediatric CNS tumors.

Atypical heat shock transcription factor HSF5 is critical for male meiotic prophase under non-stress conditions

Meiotic prophase progression is differently regulated in males and females. In males, pachytene transition during meiotic prophase is accompanied by robust alteration in gene expression. However, how gene expression is regulated differently to ensure meiotic prophase completion in males remains elusive. Herein, we identify HSF5 as a male germ cell-specific heat shock transcription factor (HSF) for meiotic prophase progression. Genetic analyzes and single-cell RNA-sequencing demonstrate that HSF5 is essential for progression beyond the pachytene stage under non-stress conditions rather than heat stress. Chromatin binding analysis in vivo and DNA-binding assays in vitro suggest that HSF5 binds to promoters in a subset of genes associated with chromatin organization. HSF5 recognizes a DNA motif different from typical heat shock elements recognized by other canonical HSFs. This study suggests that HSF5 is an atypical HSF that is required for the gene expression program for pachytene transition during meiotic prophase in males.

Myelin regulatory factor ( Myrf ) is a critical early regulator of retinal pigment epithelial development

Myelin regulatory factor (Myrf) is a critical transcription factor in early retinal and retinal pigment epithelial development, and human variants in MYRF are a cause for nanophthalmos. Single cell RNA sequencing (scRNAseq) was performed on Myrf conditional knockout mice ( Rx>Cre Myrf fl/fl ) at 3 developmental timepoints. Myrf was expressed specifically in the RPE, and expression was abrogated in Rx>Cre Myrf fl/fl eyes. scRNAseq analysis revealed a loss of RPE cells at all timepoints resulting from cell death. GO-term analysis in the RPE revealed downregulation of melanogenesis and anatomic structure morphogenesis pathways, which were supported by electron microscopy and histologic analysis. Novel structural target genes including Ermn and Upk3b , along with macular degeneration and inherited retinal disease genes were identified as downregulated, and a strong upregulation of TGFß/BMP signaling and effectors was observed. Regulon analysis placed Myrf downstream of Pax6 and Mitf and upstream of Sox10 in RPE differentiation. Together, these results suggest a strong role for Myrf in the RPE maturation by regulating melanogenesis, cell survival, and cell structure, in part acting through suppression of TGFß signaling and activation of Sox10 .

SUMMARY STATEMENT

Myrf regulates RPE development, melanogenesis, and is important for cell structure and survival, in part through regulation of Ermn , Upk3b and Sox10, and BMP/TGFb signaling.

The Role of JAK-STAT-SOCS1 Axis in Tumorigenesis, Malignant Progression and Lymphatic Metastasis of Penile Cancer

Background: To uncover the potential significance of JAK-STAT-SOCS1 axis in penile cancer, our study was the pioneer in exploring the altered expression processes of JAK-STAT-SOCS1 axis in tumorigenesis, malignant progression and lymphatic metastasis of penile cancer. Methods: In current study, the comprehensive analysis of JAK-STAT-SOCS1 axis in penile cancer was analyzed via multiple analysis approaches based on GSE196978 data, single-cell data (6 cancer samples) and bulk RNA data (7 cancer samples and 7 metastasis lymph nodes). Results: Our study observed an altered molecular expression of JAK-STAT-SOCS1 axis during three different stages of penile cancer, from tumorigenesis to malignant progression to lymphatic metastasis. STAT4 was an important dominant molecule in penile cancer, which mediated the immunosuppressive tumor microenvironment by driving the apoptosis of cytotoxic T cell and was also a valuable biomarker of immune checkpoint inhibitor treatment response. Conclusions: Our findings revealed that the complexity of JAK-STAT-SOCS1 axis and the predominant role of STAT4 in penile cancer, which can mediate tumorigenesis, malignant progression, and lymphatic metastasis. This insight provided valuable information for developing precise treatment strategies for patients with penile cancer.

Association of preoperative aspartate aminotransferase to platelet ratio index with outcomes and tumour microenvironment among colorectal cancer with liver metastases

This study aims to investigate applicable robust biomarkers that can improve prognostic predictions for colorectal liver metastasis (CRLM) patients receiving simultaneous resection. A total of 1323 CRLM patients from multiple centres were included. The preoperative aspartate aminotransferase to platelet ratio index (APRI) level from blood of patients were obtained. Patients were stratified into a high APRI group and a low APRI group, and comparisons were conducted by analyzing progression-free survival (PFS), overall survival (OS) and postoperative early recurrence. Tumour samples of CRLM were collected to perform single-cell RNA sequencing and multiplex immunohistochemistry/immunofluorescence (mIHC/IF) to investigate the association of APRI levels and the tumour microenvironment of CRLM. Compared with APRI <0.33, PFS disadvantage (IPTW-adjusted HR = 1.240, P = 0.015) and OS disadvantage (IPTW- adjusted HR = 1.507, P = 0.002) of APRI ≥0.33 were preserved in the IPTW-adjusted Cox hazards regression analyses. An APRI ≥0.25 was associated with a significantly increased risk of postoperative early recurrence after adjustment (IPTW-adjusted OR = 1.486, P = 0.001). The external validation showed consistent results with the training cohort. In the high APRI group, the single-cell RNA sequencing results revealed a heightened malignancy of epithelial cells, the enrichment of inflammatory-like cancer-associated fibroblasts and SPP1+ macrophages associated with activation of malignant cells and fibrotic microenvironment, and a more suppressed-function T cells; mIHC/IF showed that PD1+ CD4+ T cells, FOXP3+ CD4+ T cells, PD1+ CD8+ T cells, FOXP3+ CD8+ T cells, SPP1+ macrophages and iCAFs were significantly increased in the intratumoral region and peritumoral region. This study contributed valuable evidence regarding preoperative APRI for predicting prognoses among CRLM patients receiving simultaneous resection and provided underlying clues supporting the association between APRI and clinical outcomes by single-cell sequencing bioinformatics analysis and mIHC/IF.

Cancer SLC6A6-mediated taurine uptake transactivates immune checkpoint genes and induces exhaustion in CD8+ T cells

Taurine is used to bolster immunity, but its effects on antitumor immunity are unclear. Here, we report that cancer-related taurine consumption causes T cell exhaustion and tumor progression. The taurine transporter SLC6A6 is correlated with aggressiveness and poor outcomes in multiple cancers. SLC6A6-mediated taurine uptake promotes the malignant behaviors of tumor cells but also increases the survival and effector function of CD8+ T cells. Tumor cells outcompete CD8+ T cells for taurine by overexpressing SLC6A6, which induces T cell death and malfunction, thereby fueling tumor progression. Mechanistically, taurine deficiency in CD8+ T cells increases ER stress, promoting ATF4 transcription in a PERK-JAK1-STAT3 signaling-dependent manner. Increased ATF4 transactivates multiple immune checkpoint genes and induces T cell exhaustion. In gastric cancer, we identify a chemotherapy-induced SP1-SLC6A6 regulatory axis. Our findings suggest that tumoral-SLC6A6-mediated taurine deficiency promotes immune evasion and that taurine supplementation reinvigorates exhausted CD8+ T cells and increases the efficacy of cancer therapies.

Single-cell transcriptomics reveal distinct immune-infiltrating phenotypes and macrophage-tumor interaction axes among different lineages of pituitary neuroendocrine tumors

BACKGROUND

Pituitary neuroendocrine tumors (PitNETs) are common gland neoplasms demonstrating distinctive transcription factors. Although the role of immune cells in PitNETs has been widely recognized, the precise immunological environment and its control over tumor cells are poorly understood.

METHODS

The heterogeneity, spatial distribution, and clinical significance of macrophages in PitNETs were analyzed using single-cell RNA sequencing (scRNA-seq), bulk RNA-seq, spatial transcriptomics, immunohistochemistry, and multiplexed quantitative immunofluorescence (QIF). Cell viability, cell apoptosis assays, and in vivo subcutaneous xenograft experiments have confirmed that INHBA-ACVR1B influences the process of tumor cell apoptosis.

RESULTS

The present study evaluated scRNA-seq data from 23 PitNET samples categorized into 3 primary lineages. The objective was to explore the diversity of tumors and the composition of immune cells across these lineages. Analyzed data from scRNA-seq and 365 bulk RNA sequencing samples conducted in-house revealed the presence of three unique subtypes of tumor immune microenvironment (TIME) in PitNETs. These subtypes were characterized by varying levels of immune infiltration, ranging from low to intermediate to high. In addition, the NR5A1 lineage is primarily associated with the subtype characterized by limited infiltration of immune cells. Tumor-associated macrophages (TAMs) expressing CX3CR1+, C1Q+, and GPNMB+ showed enhanced contact with tumor cells expressing NR5A1 + , TBX19+, and POU1F1+, respectively. This emphasizes the distinct interaction axes between TAMs and tumor cells based on their lineage. Moreover, the connection between CX3CR1+ macrophages and tumor cells via INHBA-ACVR1B regulates tumor cell apoptosis.

CONCLUSIONS

In summary, the different subtypes of TIME and the interaction between TAM and tumor cells offer valuable insights into the control of TIME that affects the development of PitNET. These findings can be utilized as prospective targets for therapeutic interventions.

DPSCs regulate epithelial-T cell interactions in oral submucous fibrosis

BACKGROUND

Oral submucous fibrosis (OSF) is a precancerous lesion characterized by fibrous tissue deposition, the incidence of which correlates positively with the frequency of betel nut chewing. Prolonged betel nut chewing can damage the integrity of the oral mucosal epithelium, leading to chronic inflammation and local immunological derangement. However, currently, the underlying cellular events driving fibrogenesis and dysfunction are incompletely understood, such that OSF has few treatment options with limited therapeutic effectiveness. Dental pulp stem cells (DPSCs) have been recognized for their anti-inflammatory and anti-fibrosis capabilities, making them promising candidates to treat a range of immune, inflammatory, and fibrotic diseases. However, the application of DPSCs in OSF is inconclusive. Therefore, this study aimed to explore the pathogenic mechanism of OSF and, based on this, to explore new treatment options.

METHODS

A human cell atlas of oral mucosal tissues was compiled using single-cell RNA sequencing to delve into the underlying mechanisms. Epithelial cells were reclustered to observe the heterogeneity of OSF epithelial cells and their communication with immune cells. The results were validated in vitro, in clinicopathological sections, and in animal models. In vivo, the therapeutic effect and mechanism of DPSCs were characterized by histological staining, immunohistochemical staining, scanning electron microscopy, and atomic force microscopy.

RESULTS

A unique epithelial cell population, Epi1.2, with proinflammatory and profibrotic functions, was predominantly found in OSF. Epi1.2 cells also induced the fibrotic process in fibroblasts by interacting with T cells through receptor-ligand crosstalk between macrophage migration inhibitory factor (MIF)-CD74 and C-X-C motif chemokine receptor 4 (CXCR4). Furthermore, we developed OSF animal models and simulated the clinical local injection process in the rat buccal mucosa using DPSCs to assess their therapeutic impact and mechanism. In the OSF rat model, DPSCs demonstrated superior therapeutic effects compared with the positive control (glucocorticoids), including reducing collagen deposition and promoting blood vessel regeneration. DPSCs mediated immune homeostasis primarily by regulating the numbers of KRT19 + MIF + epithelial cells and via epithelial-stromal crosstalk.

CONCLUSIONS

Given the current ambiguity surrounding the cause of OSF and the limited treatment options available, our study reveals that epithelial cells and their crosstalk with T cells play an important role in the mechanism of OSF and suggests the therapeutic promise of DPSCs.