The proneural transcription factor atonal basic helix-loop-helix transcription factor 7 (ATOH7) is expressed in early progenitors in the developing neuroretina. In vertebrates, this is crucial for the development of retinal ganglion cells (RGCs), as mutant animals show an almost complete absence of RGCs, underdeveloped optic nerves, and aberrations in retinal vessel development. Human mutations are rare and result in autosomal recessive optic nerve hypoplasia (ONH) or severe vascular changes, diagnosed as autosomal recessive persistent hyperplasia of the primary vitreous (PHPVAR). To better understand the role of ATOH7 in neuroretinal development, we created ATOH7 knockout and eGFP-expressing ATOH7 reporter human induced pluripotent stem cells (hiPSCs), which were differentiated into early-stage retinal organoids. Target loci regulated by ATOH7 were identified by Cleavage Under Targets and Release Using Nuclease with sequencing (CUT&RUN-seq) and differential expression by RNA sequencing (RNA-seq) of wildtype and mutant organoid-derived reporter cells. Additionally, single-cell RNA sequencing (scRNA-seq) was performed on whole organoids to identify cell type-specific genes. Mutant organoids displayed substantial deficiency in axon sprouting, reduction in RGCs, and an increase in other cell types. We identified 469 differentially expressed target genes, with an overrepresentation of genes belonging to axon development/guidance and Notch signaling. Taken together, we consolidate the function of human ATOH7 in guiding progenitor competence by inducing RGC-specific genes while inhibiting other cell fates. Furthermore, we highlight candidate genes responsible for ATOH7-associated optic nerve and retinovascular anomalies, which sheds light to potential future therapy targets for related disorders.

Fibrous dysplasia (FD) is a rare bone disorder characterized by the replacement of normal bone with benign fibro-osseous tissue. Developments in our understanding of the pathophysiology and treatment options are impeded by the lack of suitable research models. In this study, we developed an in vitro organotypic model capable of recapitulating key intrinsic and phenotypic properties of FD. Initially, transcriptomic profiling of individual cells isolated from patient lesional tissues unveiled intralesional molecular and cellular heterogeneity. Leveraging these insights, we established patient-derived organoids (PDOs) using primary cells obtained from patient FD lesions. Evaluation of PDOs demonstrated preservation of fibrosis-associated constituent cell types and transcriptional signatures observed in FD lesions. Additionally, PDOs retained distinct constellations of genomic and metabolic alterations characteristic of FD. Histological evaluation further corroborated the fidelity of PDOs in recapitulating important phenotypic features of FD that underscore their pathophysiological relevance. Our findings represent meaningful progress in the field, as they open up the possibility for in vitro modeling of rare bone lesions in a three-dimensional context and may signify the first step towards creating a personalized platform for research and therapeutic studies.

Diabetic nephropathy (DN) is a severe complication of diabetes, characterized by changes in kidney structure and function. The natural product rosmarinic acid (RA) has demonstrated therapeutic effects, including anti-inflammation and anti-oxidative-stress, in renal damage or dysfunction. In this study, we characterized the heterogeneity of the cellular response in kidneys to DN-induced injury and RA treatment at single cell levels. Our results demonstrated that RA significantly alleviated renal tubular epithelial injury, particularly in the proximal tubular S1 segment and on glomerular epithelial cells known as podocytes, while attenuating the inflammatory response of macrophages, oxidative stress, and cytotoxicity of natural killer cells. These findings provide a comprehensive understanding of the mechanisms by which RA alleviates kidney damage, oxidative stress, and inflammation, offering valuable guidance for the clinical application of RA in the treatment of DN.

Patient-derived organoids have proven to be a highly relevant model for evaluating of disease mechanisms and drug efficacies, as they closely recapitulate in vivo physiology. Colorectal cancer organoids, specifically, exhibit a diverse range of morphologies, which have been analyzed with image-based profiling. However, the relationship between morphological subtypes and functional parameters of the organoids remains underexplored. Here, we identified two distinct morphological subtypes (\"cystic\" and \"solid\") across 31360 bright field images using image-based profiling, which correlated differently with viability and apoptosis level of colorectal cancer organoids. Leveraging object detection neural networks, we were able to categorize single organoids achieving higher viability scores as \"cystic\" than \"solid\" subtype. Furthermore, a deep generative model was proposed to predict apoptosis intensity based on a apoptosis-featured dataset encompassing over 17000 bright field and matched fluorescent images. Notably, a significant correlation of 0.91 between the predicted value and ground truth was achived, underscoring the feasibility of this generative model as a potential means for assessing organoid functional parameters. The underlying cellular heterogeneity of the organoids, i.e., conserved colonic cell types and rare immune components, was also verified with scRNA sequencing, implying a compromised tumor microenvironment. Additionally, the \"cystic\" subtype was identified as a relapse phenotype featuring intestinal stem cell signatures, suggesting that this visually discernible relapse phenotype shows potential as a novel biomarker for colorectal cancer diagnosis and prognosis. In summary, our findings demonstrate that the morphological heterogeneity of colorectal cancer organoids explicitly recapitulate the association of phenotypic features and exogenous perturbations through the image-based profiling, providing new insights into disease mechanisms.

UNASSIGNED: Cuproptosis is a copper-dependent cell death that is connected to the development and immune response of multiple diseases. However, the function of cuproptosis in the immune characteristics of sepsis remains unclear.
UNASSIGNED: We obtained two sepsis datasets (GSE9960 and GSE134347) from the GEO database and classified the raw data with R packages. Cuproptosis-related genes were manually curated, and differentially expressed cuproptosis-related genes (DECuGs) were identified. Afterwards, we applied enrichment analysis and identified key DECuGs by performing machine learning techniques. Then, the immune cell infiltrations and correlation between DECuGs and immunocyte features were created by the CIBERSORT algorithm. Subsequently, unsupervised hierarchical clustering analysis was performed based on key DECuGs. We then constructed a ceRNA network based on key DECuGs by using multi-step computational strategies and predicted potential drugs in the DrugBank database. Finally, the role of these key genes in immune cells was validated at the single-cell RNA level between septic patients and healthy controls.
UNASSIGNED: Overall, 16 DECuGs were obtained, and most of them had lower expression levels in sepsis samples. Afterwards, we obtained six key DECuGs by performing machine learning. Then, the LIPT1-T-cell CD4 memory resting was the most positively correlated DECuG-immunocyte pair. Subsequently, two different subclusters were identified by six DECuGs. Bioinformatics analysis revealed that there were different immune characteristics between the two subclusters. Moreover, we identified the key lncRNA OIP5-AS1 within the ceRNA network and obtained 4 drugs that may represent novel drugs for sepsis. Finally, these key DECuGs were statistically significantly dysregulated in another validation set and showed a major distribution in monocytes, T cells, B cells, NK cells and platelets at the single-cell RNA level.
UNASSIGNED: These findings suggest that cuproptosis might promote the progression of sepsis by affecting the immune system and metabolic dysfunction, which provides a new direction for understanding potential pathogenic processes and therapeutic targets in sepsis.

Epidermal cells are the main avenue for signal and material exchange between plants and the environment. Leaf epidermal cells primarily include pavement cells, guard cells, and trichome cells. The development and distribution of different epidermal cells are tightly regulated by a complex transcriptional regulatory network mediated by phytohormones, including jasmonic acid, and transcription factors. How the fate of leaf epidermal cells is determined, however, is still largely unknown due to the diversity of cell types and the complexity of their regulation. Here, we characterized the transcriptional profiles of epidermal cells in 3-day-old true leaves of Arabidopsis thaliana using single-cell RNA sequencing. We identified two genes encoding BASIC LEUCINE-ZIPPER (bZIP) transcription factors, namely bZIP25 and bZIP53, which are highly expressed in pavement cells and early-stage meristemoid cells. Densities of pavement cells and trichome cells were found to increase and decrease, respectively, in bzip25 and bzip53 mutants, compared with wild-type plants. This trend was more pronounced in the presence of jasmonic acid, suggesting that these transcription factors regulate the development of trichome cells and pavement cells in response to jasmonic acid.

The incidence of Crohn\'s disease (CD) continues to increase worldwide. The contribution of CD4+ cell populations remains to be elucidated. We aimed to provide an in-depth transcriptional assessment of CD4+ T cells driving chronic inflammation in CD.
We performed single-cell RNA-sequencing in CD4+ T cells isolated from ileal biopsies of patients with CD compared with healthy individuals. Cells underwent clustering analysis, followed by analysis of gene signaling networks. We overlapped our differentially expressed genes with publicly available microarray data sets and performed functional in vitro studies, including an in vitro suppression assay and organoid systems, to model gene expression changes observed in CD regulatory T (Treg) cells and to test predicted therapeutics.
We identified 5 distinct FOXP3+ regulatory Treg subpopulations. Tregs isolated from healthy controls represent the origin of pseudotemporal development into inflammation-associated subtypes. These proinflammatory Tregs displayed a unique responsiveness to tumor necrosis factor-α signaling with impaired suppressive activity in vitro and an elevated cytokine response in an organoid coculture system. As predicted in silico, the histone deacetylase inhibitor vorinostat normalized gene expression patterns, rescuing the suppressive function of FOXP3+ cells in vitro.
We identified a novel, proinflammatory FOXP3+ T cell subpopulation in patients with CD and developed a pipeline to specifically target these cells using the US Food and Drug Administration-approved drug vorinostat.

Evading immune surveillance is a hallmark for the development of multiple cancer types. Whether immune evasion contributes to the pathogenesis of high-grade prostate cancer (HGPCa) remains an area of active inquiry.
Through single-cell RNA sequencing and multicolor flow cytometry of freshly isolated prostatectomy specimens and matched peripheral blood, we aimed to characterize the tumor immune microenvironment (TME) of localized prostate cancer (PCa), including HGPCa and low-grade prostate cancer (LGPCa).
HGPCa are highly infiltrated by exhausted CD8+ T cells, myeloid cells, and regulatory T cells (TRegs). These HGPCa-infiltrating CD8+ T cells expressed high levels of exhaustion markers including TIM3, TOX, TCF7, PD-1, CTLA4, TIGIT, and CXCL13. By contrast, a high ratio of activated CD8+  effector T cells relative to TRegs and myeloid cells infiltrate the TME of LGPCa. HGPCa CD8+  tumor-infiltrating lymphocytes (TILs) expressed more androgen receptor and prostate-specific membran antigen yet less prostate-specific antigen than the LGPCa CD8+  TILs. The PCa TME was infiltrated by macrophages but these did not clearly cluster by M1 and M2 markers.
Our study reveals a suppressive TME with high levels of CD8+ T cell exhaustion in localized PCa, a finding enriched in HGPCa relative to LGPCa. These studies suggest a possible link between the clinical-pathologic risk of PCa and the associated TME. Our results have implications for our understanding of the immunologic mechanisms of PCa pathogenesis and the implementation of immunotherapy for localized PCa.

Multiple sclerosis (MS) is an autoimmune disease of the central nervous system still lacking a cure. Treatment typically focuses on slowing the progression and managing MS symptoms. Single-cell transcriptomics allows the investigation of the immune system-the key player in MS onset and development-in great detail increasing our understanding of MS mechanisms and stimulating the discovery of the targets for potential therapies. Still, de novo drug development takes decades; however, this can be reduced by drug repositioning. A promising approach is to select potential drugs based on activated or inhibited genes and pathways. In this study, we explored the public single-cell RNA data from an experiment with six patients on single-cell RNA peripheral blood mononuclear cells (PBMC) and cerebrospinal fluid cells (CSF) of patients with MS and idiopathic intracranial hypertension. We demonstrate that AIM2 inflammasome, SMAD2/3 signaling, and complement activation pathways are activated in MS in different CSF and PBMC immune cells. Using genes from top-activated pathways, we detected several promising small molecules to reverse MS immune cells\' transcriptomic signatures, including AG14361, FGIN-1-27, CA-074, ARP 101, Flunisolide, and JAK3 Inhibitor VI. Among these molecules, we also detected an FDA-approved MS drug Mitoxantrone, supporting the reliability of our approach.

The parabrachial nucleus (PBN) is a major hub that receives sensory information from both internal and external environments. Specific populations of PBN neurons are involved in behaviors including food and water intake, nociceptive responses, breathing regulation, as well as learning and responding appropriately to threatening stimuli. However, it is unclear how many PBN neuron populations exist and how different behaviors may be encoded by unique signaling molecules or receptors. Here we provide a repository of data on the molecular identity, spatial location, and projection patterns of dozens of PBN neuron subclusters. Using single-cell RNA sequencing, we identified 21 subclusters of neurons in the PBN and neighboring regions. Multiplexed in situ hybridization showed many of these subclusters are enriched within specific PBN subregions with scattered cells in several other regions. We also provide detailed visualization of the axonal projections from 21 Cre-driver lines of mice. These results are all publicly available for download and provide a foundation for further interrogation of PBN functions and connections.