Search

Methamphetamine is a highly addictive stimulant known to cause neurotoxicity, cognitive deficits, and immune dysregulation in the brain. Despite significant research, the molecular mechanisms driving methamphetamine-induced neurotoxicity and glial cell dysfunction remain poorly understood. This study investigates how methamphetamine disrupts glial cell function and contributes to neurodevelopmental and neurodegenerative processes. Using single-cell RNA sequencing (scRNA-seq), we analyzed the transcriptomes of 4000 glial cell-associated genes from the cortical regions of mice chronically administered methamphetamine. Methamphetamine exposure altered the key pathways in astrocytes, including the circadian rhythm and cAMP signaling; in microglia, affecting autophagy, ubiquitin-mediated proteolysis, and mitophagy; and in oligodendrocytes, disrupting lysosomal function, cytoskeletal regulation, and protein processing. Notably, several transcription factors, such as Zbtb16, Hif3a, Foxo1, and Klf9, were significantly dysregulated in the glial cells. These findings reveal profound methamphetamine-induced changes in the glial transcriptomes, particularly in the cortical regions, highlighting potential molecular pathways and transcription factors as targets for therapeutic intervention. This study provides novel insights into the glial-mediated mechanisms of methamphetamine toxicity, contributing to our understanding of its effects on the central nervous system and laying the groundwork for future strategies to mitigate its neurotoxic consequences. [ABSTRACT FROM AUTHOR]

(1) Background: Hepatoblastoma and medulloblastoma are two types of pediatric tumors with embryonic origins. Both tumor types can exhibit genetic alterations that affect the β-catenin and Wnt pathways; (2) Materials and Methods: This study used bioinformatics and integrative analysis of multi-omics data at both the tumor and single-cell levels to investigate two distinct pediatric tumors: medulloblastoma and hepatoblastoma; (3) Results: The cross-transcriptome analysis revealed a commonly regulated expression signature between hepatoblastoma and medulloblastoma tumors. Among the commonly upregulated genes, the transcription factor LEF1 was significantly expressed in both tumor types. In medulloblastoma, LEF1 upregulation is associated with the WNT-subtype. The analysis of LEF1 genome binding occupancy in H1 embryonic stem cells identified 141 LEF1 proximal targets activated in WNT medulloblastoma, 13 of which are involved in Wnt pathway regulation: RNF43, LEF1, NKD1, AXIN2, DKK4, DKK1, LGR6, FGFR2, NXN, TCF7L1, STK3, YAP1, and NFATC4. The ROC curve analysis of the combined expression of these 13 WNT-related LEF1 targets yielded an area under the curve (AUC) of 1.00, indicating 100% specificity and sensitivity for predicting the WNT subtype in the PBTA medulloblastoma cohort. An expression score based on these 13 WNT-LEF1 targets accurately predicted the WNT subtype in two independent medulloblastoma transcriptome cohorts. At the single-cell level, the WNT-LEF1 expression score was exclusively positive in WNT-medulloblastoma tumor cells. This WNT-LEF1-dependent signature was also confirmed as activated in the hepatoblastoma tumor transcriptome. At the single-cell level, the WNT-LEF1 expression score was higher in tumor cells from both human hepatoblastoma samples and a hepatoblastoma patient-derived xenotransplant model; (4) Discussion: This study uncovered a shared transcriptional activation of a LEF1-dependent embryonic program, which orchestrates the regulation of the Wnt signaling pathway in tumor cells from both hepatoblastoma and medulloblastoma. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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.

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. [ABSTRACT FROM AUTHOR]

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. Single-cell RNA sequencing (scRNA-seq) analyses have provided a novel insight into cell-specific gene expression changes in diabetic nephropathy (DP) mice with rosmarinic acid (RA) treatment. Here, this study reveals RA could reduce oxidative stress and inflammation in the PT segment and podocytes, and mitigate inflammation in M1-like macrophages and natural killer cells. This study further deepens the understanding of novel mechanisms and potential therapeutic targets of RA on DP. [Display omitted] [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Abstract Molecular profiles of neurons influence neural development and function but bridging the gap between genes, circuits, and behavior has been very difficult. Here we used single cell RNAseq to generate a complete gene expression atlas of the Drosophila larval central nervous system composed of 131,077 single cells across three developmental stages (1 h, 24 h and 48 h after hatching). We identify 67 distinct cell clusters based on the patterns of gene expression. These include 31 functional mature larval neuron clusters, 1 ring gland cluster, 8 glial clusters, 6 neural precursor clusters, and 13 developing immature adult neuron clusters. Some clusters are present across all stages of larval development, while others are stage specific (such as developing adult neurons). We identify genes that are differentially expressed in each cluster, as well as genes that are differentially expressed at distinct stages of larval life. These differentially expressed genes provide promising candidates for regulating the function of specific neuronal and glial types in the larval nervous system, or the specification and differentiation of adult neurons. The cell transcriptome Atlas of the Drosophila larval nervous system is a valuable resource for developmental biology and systems neuroscience and provides a basis for elucidating how genes regulate neural development and function.

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. [ABSTRACT FROM AUTHOR]

Cancer-associated fibroblasts (CAFs) represent a central cell population of the tumor microenvironment (TME). Recently, single-cell RNA-sequencing (scRNA-seq) analyses of primary tumors of different cancer entities yielded different classifications of CAF subsets underscoring the heterogeneity of CAFs within the TME. Here, we analyzed the transcriptional signatures of approximately 8400 CAFs and normal fibroblasts by scRNA-seq and compared genetic profiles of CAFs from murine melanoma primary tumors to CAFs from corresponding melanoma lung metastases. This revealed distinct subsets for primary tumor and metastasis-specific CAF populations, respectively. Combined with the spatial characterization of metastasis CAFs at the RNA and protein level, scRNA analyses indicate tumor-dependent crosstalk between neutrophils and CAFs, mediated via SAA3 and IL1b-related signaling pathways, which can be recapitulated in vitro. Analyzing tissue sections of human patient samples, this interaction was found to be present in human melanoma metastasis. Taken together, our data highlight unique characteristics of metastasis CAFs with potential therapeutic impact for melanoma metastasis., (© 2025 The Author(s). Pigment Cell & Melanoma Research published by John Wiley & Sons Ltd.)

Two of the most common neurodegenerative disorders – Alzheimer’s and Parkinson’s diseases – are characterized by synaptic dysfunction and degeneration that culminate in neuronal loss due to abnormal protein accumulation. The intracellular aggregation of hyper-phosphorylated tau and the extracellular aggregation of amyloid beta plaques form the basis of Alzheimer’s disease pathology. The major hallmark of Parkinson’s disease is the loss of dopaminergic neurons in the substantia nigra pars compacta, following the formation of Lewy bodies, which consists primarily of alpha-synuclein aggregates. However, the discrete mechanisms that contribute to neurodegeneration in these disorders are still poorly understood. Both neuronal loss and impaired adult neurogenesis have been reported in animal models of these disorders. Yet these findings remain subject to frequent debate due to a lack of conclusive evidence in post mortem brain tissue from human patients. While some publications provide significant findings related to axonal regeneration in Alzheimer’s and Parkinson’s diseases, they also highlight the limitations and obstacles to the development of neuroregenerative therapies. In this review, we summarize in vitro and in vivo findings related to neurogenesis, neuroregeneration and neurodegeneration in the context of Alzheimer’s and Parkinson’s diseases.

More and more researchers use single-cell RNA sequencing (scRNA-seq) technology to characterize the transcriptional map at the single-cell level. They use it to study the heterogeneity of complex tissues, transcriptome dynamics, and the diversity of unknown organisms. However, there are generally lots of technical and biological noises in the scRNA-seq data since the randomness of gene expression patterns. These data are often characterized by high-dimension, sparsity, large number of "dropout" values, and affected by batch effects. A large number of "dropout" values in scRNA-seq data seriously conceal the important relationship between genes and hinder the downstream analysis. Therefore, the imputation of dropout values of scRNA-seq data is particularly important. We classify, analyze and compare the current advanced scRNA-seq data imputation methods from different angles. Through the comparison and analysis of the principle, advantages and disadvantages of the algorithm, it can provide suggestions for the selection of imputation methods for specific problems and diverse data, and have basic research significance for the downstream function analysis of data. [ABSTRACT FROM AUTHOR]

The importance of T cells in controlling SARS-CoV-2 infections has been demonstrated widely, but insights into the quality of these responses are still limited due to technical challenges. Indeed, understanding the functionality of the T-cell receptor (TCR) repertoire of a polyclonal antigen-specific population still requires the tedious work of T-cell cloning or TCR re-expression and subsequent characterization. In this work, we show that it is possible to discriminate highly functional and bystander TCRs based on gene signatures of T-cell activation induced by recent peptide stimulation. SARS-CoV-2-specific TCRs previously identified by cytokine release after peptide restimulation and subsequent single-cell RNA sequencing were re-expressed via CRISPR-Cas9-mediated gene editing into a Jurkat-based reporter cell line system suitable for high-throughput screening. We could observe differences in SARS-CoV-2 epitope recognition as well as a wide range of functional avidities. By correlating these in vitro TCR engineered functional data with the transcriptomic profiles of the corresponding TCR-expressing parental T cells, we could validate that gene signatures of recent T-cell activation accurately identify and predict truly SARS-CoV-2-specific TCRs. In summary, this work paves the way for alternative approaches useful for the functional analysis of global antigen-specific TCR repertoires with largely improved throughput. [ABSTRACT FROM AUTHOR]

Quantifying gene expression in individual cells can substantially improve our understanding about complex genetically engineered cell products such as chimeric antigen receptor (CAR) T cells. Here we designed a single-cell RNA sequencing (scRNA-seq) approach to monitor the delivery of a CD19-CAR gene via lentiviral vectors (LVs), i.e., the conventional vesicular stomatitis virus (VSV)-LV and the CD8-targeted CD8-LV. LV-exposed human donor peripheral blood mononuclear cells (PBMCs) were evaluated for a panel of 400 immune response-related genes including LV-specific probes. The resulting data revealed a trimodal expression for the CAR and CD8A, demanding a careful distribution-based identification of CAR T cells and CD8+ lymphocytes in scRNA-seq analysis. The fraction of T cells expressing high CAR levels was in concordance with flow cytometry results. More than 97% of the cells hit by CD8-LV expressed the CD8A gene. Remarkably, the majority of the potential off-target cells were in fact on-target cells, resulting in a target cell selectivity of more than 99%. Beyond that, differential gene expression analysis revealed the upregulation of restriction factors in CAR-negative cells, thus explaining their protection from CAR gene transfer. In summary, we provide a workflow and subsetting approach for scRNA-seq enabling reliable distinction between transduced and untransduced cells during CAR T cell generation.

Abstract Background Epidemiological data associate high levels of combustion-derived particulate matter (PM) with deleterious respiratory outcomes, but the mechanism underlying those outcomes remains elusive. It has been acknowledged by the World Health Organization that PM exposure contributes to more than 4.2 million all-cause mortalities worldwide each year. Current literature demonstrates that PM exacerbates respiratory diseases, impairs lung function, results in chronic respiratory illnesses, and is associated with increased mortality. The proposed mechanisms revolve around oxidative stress and inflammation promoting pulmonary physiological remodeling. However, our previous data found that PM is capable of inducing T helper cell 17 (Th17) immune responses via aryl hydrocarbon receptor (Ahr) activation, which was associated with neutrophilic invasion characteristic of steroid insensitive asthma. Methods In the present study, we utilized a combination of microarray and single cell RNA sequencing data to analyze the immunological landscape in mouse lungs following acute exposure to combustion derived particulate matter. Results We present data that suggest epithelial cells produce specific cytokines in the aryl hydrocarbon receptor (Ahr) pathway that inform dendritic cells to initiate the production of pathogenic T helper (eTh17) cells. Using single-cell RNA sequencing analysis, we observed that upon exposure epithelial cells acquire a transcriptomic profile indicative of increased Il-17 signaling, Ahr activation, Egfr signaling, and T cell receptor and co-stimulatory signaling pathways. Epithelial cells further showed, Ahr activation is brought on by Ahr/ARNT nuclear translocation and activation of tyrosine kinase c-src, Egfr, and subsequently Erk1/2 pathways. Conclusions Collectively, our data corroborates that PM initiates an eTh17 specific inflammatory response causing neutrophilic asthma through pathways in epithelial, dendritic, and T cells that promote eTh17 differentiation during initial PM exposure.

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.

Duck reovirus (DRV) is a universal waterfowl virus that causes significant economic losses in the duck industry. However, the role of the host innate immune response of the Bursa of Fabricius to DRV infection is largely unknown. In the present study, we constructed a single-cell resolution transcriptomic atlas of the Bursa of Fabricius of Cairna moschata after infection with HN10 (a novel DRV). Ten cell-type marker genes were used to annotate the cell type, indicating a high degree of cell heterogeneity in the Bursa of Fabricius. Most of the innate and adaptive immune system-related genes were highly expressed in T cells, B cells, neutrophils, macrophages, and DCs. In the Bursa of Fabricius, the proportions of DCs and macrophages were largely increased by HN10 infection at 14 d, suggesting that DCs and macrophages play important roles in the long-term viral response. Notably, a number of innate and adaptive immune system-related genes were highly expressed at 24 h after HN10 infection, indicating that the Bursa of Fabricius has a very strong immune function even in the early developmental stage. In the immune system, the NOD-like receptor signaling pathway and RIG-I-like receptor signaling pathway were significantly activated at the early stage of HN10 infection, while the Toll-like receptor signaling pathway was significantly activated at the late stage. Enrichment analysis suggested that different immune signaling pathways play roles in specific developmental stages. Our data provide an opportunity to reveal the immune response to DRV infection at the single-cell level. [ABSTRACT FROM AUTHOR]

• Significant differences in ICM expression on TILs between HPVpos and HPVneg OPSCC. • Different immune-therapeutic strategies may be needed managing HPVpos OPSCC. • Combined therapies might be needed to address distinct ICM expressions based on HPV. A substantial proportion of head and neck squamous cell carcinoma (HNSCC), particularly oropharyngeal squamous cell carcinoma (OPSCC), is associated with human papillomavirus (HPV), resulting in distinct molecular phenotypes. In this study, we investigated differential immune checkpoint molecule (ICM) expression by HPV status using RNA sequencing data to identify additional ICM targets that may complement anti-PD1 antibodies. RNA sequencing was performed on 51 OPSCC cases and validated using the TCGA HNSCC dataset. Unsupervised clustering and differential gene expression analyses in R were conducted based on HPV status. Additionally, a published single-cell RNA sequencing (scRNA) dataset of tumor-infiltrating lymphocytes (TIL) and peripheral immune cells (PBMC) (GSE139324) was analyzed with a Seurat pipeline grouped by HPV status. Our study identified a significant upregulation of all examined ICM in HPV-positive OPSCC through bulk RNA sequencing, validated by the TCGA cohort. Unsupervised clustering revealed a strong association between HPV-positive/-negative and high/low ICM expression cases respectively, indicating overlap between ICM and HPV status. In scRNA analysis, CD27, PD-1, OX-40, and BTLA were significantly more highly expressed on TILs of HPV-positive OPSCC. Conversely, VSIR was increased in PBMC and TILs of HPV-negative OPSCC, while LAG3 expression on PBMC was reduced in HPV-negative OPSCC. Our study unveils the intricate interplay of ICMs in OPSCC, emphasizing the necessity for personalized therapeutic approaches based on HPV status and immune profiles. The identified ICMs, including PD1, CD27, and CTLA4, are promising candidates for further investigation and may enhance immunotherapeutic interventions in the HPV-dependent treatment strategies for OPSCC. [ABSTRACT FROM AUTHOR]

Molecular profiles of neurons influence neural development and function but bridging the gap between genes, circuits, and behavior has been very difficult. Here we used single cell RNAseq to generate a complete gene expression atlas of the Drosophila larval central nervous system composed of 131,077 single cells across three developmental stages (1 h, 24 h and 48 h after hatching). We identify 67 distinct cell clusters based on the patterns of gene expression. These include 31 functional mature larval neuron clusters, 1 ring gland cluster, 8 glial clusters, 6 neural precursor clusters, and 13 developing immature adult neuron clusters. Some clusters are present across all stages of larval development, while others are stage specific (such as developing adult neurons). We identify genes that are differentially expressed in each cluster, as well as genes that are differentially expressed at distinct stages of larval life. These differentially expressed genes provide promising candidates for regulating the function of specific neuronal and glial types in the larval nervous system, or the specification and differentiation of adult neurons. The cell transcriptome Atlas of the Drosophila larval nervous system is a valuable resource for developmental biology and systems neuroscience and provides a basis for elucidating how genes regulate neural development and function. [ABSTRACT FROM AUTHOR]

Diagnosis and longitudinal monitoring of neurological diseases are limited by the poor specificity and limited resolution of currently available techniques. Analysis of circulating cells in cerebrospinal fluid (CSF) has emerged as a promising strategy for the diagnosis, molecular characterization, and monitoring of neurological disease. In comparison to bulk sequencing analysis, single-cell sequencing studies can provide novel insights into rare cell populations and uncover heterogeneity in gene expression at a single-cell resolution, which has several implications for understanding disease pathology and treatment. Parallel development of standardized biofluid collection protocols, pre-processing strategies, reliable single-cell isolation strategies, downstream genomic analysis, and robust computational analysis is paramount for comprehensive single-cell sequencing analysis. Here we perform a comprehensive review of studies focusing on single-cell sequencing of cells in the CSF of patients with oncological or non-oncological diseases of the central nervous system. [ABSTRACT FROM AUTHOR]

Long non-coding RNAs (lncRNAs) have attracted attention for their potential roles in modulating keratinocyte differentiation and inflammatory response; however, for many identified skin-expressing lncRNAs, there is no comprehensive characterization regarding their biological roles. In addition, the reported expression profiles for lncRNAs can be ambiguous due to their low-expressing nature. The objective of this review is to utilize large scale genomic data to characterize the prominent skin-expressing lncRNAs, aiming to provide additional insights for their potential roles in the pathology of inflammatory skin of psoriasis and atopic dermatitis by integrating in vitro and in vivo data. We highlighted the different skin-expressing lncRNAs, including H19 , which is significantly down-regulated in lesional skin of AD/psoriasis and upon cytokine stimulation in keratinocytes; it is also negatively correlated with CYP1A1 (r = -0.75, p = 8 × 10−73), a gene involved in drug metabolism and skin barrier homeostasis, in keratinocytes. In addition, SPRR2C , a potential regulator that modulates IL-22 stimulation, was upregulated in both atopic dermatitis and psoriasis lesional skin and was also downstream of the IL-17A and IL-17 + TNF signaling in keratinocytes. Using scRNAseq, we further revealed the cell type specificity of lncRNAs, including basal-expressing nature of H19 in the epidermis. Interestingly, instead of having cell type specific expression profile, we found few lncRNAs that are express across different cell types in skin, including MALAT1 , NEAT1 , and GAS5. While lncRNAs in general have lower expression, our results combining in vitro and in vivo experimental data demonstrate how some of these lncRNAs can play mediator roles in the cytokine-stimulated pathway. [ABSTRACT FROM AUTHOR]

Osteosarcoma is a type of highly aggressive bone tumor arising from primitive cells of mesenchymal origin in adults and is associated with a high rate of tumor relapse. However, there is an urgent need to clarify the molecular mechanisms underlying osteosarcoma development. The present study performed integrated bioinformatics analysis in a single-cell RNA sequencing dataset and explored the potential interactive signaling pathways associated with osteosarcoma development. Single-cell transcriptomic analysis of osteosarcoma tissues was performed by using the Seurat R package, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of differentially expressed genes was performed by using the clusterProfiler R package, and the cell–cell interaction analysis was performed by using the CellPhoneDB package. Our results showed that 11 clustered cell types were identified across 11 osteosarcoma tissues, with cell types including "osteoblastic", "myeloid", "osteoblastic_proli", "osteoclast", and "tumor-infiltrating lymphocytes (TILs)" as the main types. The DEGs between different cell types from primary, metastatic, and recurrent osteosarcomas were mainly enriched in the GO terms including "negative regulation of hydrolase activity", "regulation of peptidase activity", "regulation of binding", "negative regulation of proteolysis", and "negative regulation of peptidase activity" and in the KEGG pathways including "transcriptional misregulation in cancer", "cellular senescence", "apoptosis", "FoxO signaling pathway", "cell cycle", "NF-kappa B signaling pathway", "p53 signaling pathway", "pentose phosphate pathway", and "protein export". For the cell–cell communication network analysis, the different interaction profiles between cell types were detected among primary, metastatic, and recurrent osteosarcomas. Further exploration of the KEGG pathway revealed that these ligand/receptor interactions may be associated with the NF-κB signaling pathway and its interacted mediators. In conclusion, the present study for the first time explored the scRNA-seq dataset in osteosarcoma, and our results revealed the 11 clustered cell types and demonstrated the novel cell–cell interactions among different cell types in primary, metastatic, and recurrent osteosarcomas. The NF-κB signaling pathway may play a key role in regulating the TME of osteosarcoma. The present study may provide new insights into understanding the molecular mechanisms of osteosarcoma pathophysiology. [ABSTRACT FROM AUTHOR]

Long non-coding RNAs (lncRNAs) have attracted attention for their potential roles in modulating keratinocyte differentiation and inflammatory response; however, for many identified skin-expressing lncRNAs, there is no comprehensive characterization regarding their biological roles. In addition, the reported expression profiles for lncRNAs can be ambiguous due to their low-expressing nature. The objective of this review is to utilize large scale genomic data to characterize the prominent skin-expressing lncRNAs, aiming to provide additional insights for their potential roles in the pathology of inflammatory skin of psoriasis and atopic dermatitis by integrating in vitro and in vivo data. We highlighted the different skin-expressing lncRNAs, including H19, which is significantly down-regulated in lesional skin of AD/psoriasis and upon cytokine stimulation in keratinocytes; it is also negatively correlated with CYP1A1 (r = -0.75, p = 8 × 10−73), a gene involved in drug metabolism and skin barrier homeostasis, in keratinocytes. In addition, SPRR2C, a potential regulator that modulates IL-22 stimulation, was upregulated in both atopic dermatitis and psoriasis lesional skin and was also downstream of the IL-17A and IL-17 + TNF signaling in keratinocytes. Using scRNAseq, we further revealed the cell type specificity of lncRNAs, including basal-expressing nature of H19 in the epidermis. Interestingly, instead of having cell type specific expression profile, we found few lncRNAs that are express across different cell types in skin, including MALAT1, NEAT1, and GAS5. While lncRNAs in general have lower expression, our results combining in vitro and in vivo experimental data demonstrate how some of these lncRNAs can play mediator roles in the cytokine-stimulated pathway.

Osteosarcoma is a type of highly aggressive bone tumor arising from primitive cells of mesenchymal origin in adults and is associated with a high rate of tumor relapse. However, there is an urgent need to clarify the molecular mechanisms underlying osteosarcoma development. The present study performed integrated bioinformatics analysis in a single-cell RNA sequencing dataset and explored the potential interactive signaling pathways associated with osteosarcoma development. Single-cell transcriptomic analysis of osteosarcoma tissues was performed by using the Seurat R package, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of differentially expressed genes was performed by using the clusterProfiler R package, and the cell–cell interaction analysis was performed by using the CellPhoneDB package. Our results showed that 11 clustered cell types were identified across 11 osteosarcoma tissues, with cell types including “osteoblastic”, “myeloid”, “osteoblastic_proli”, “osteoclast”, and “tumor-infiltrating lymphocytes (TILs)” as the main types. The DEGs between different cell types from primary, metastatic, and recurrent osteosarcomas were mainly enriched in the GO terms including “negative regulation of hydrolase activity”, “regulation of peptidase activity”, “regulation of binding”, “negative regulation of proteolysis”, and “negative regulation of peptidase activity” and in the KEGG pathways including “transcriptional misregulation in cancer”, “cellular senescence”, “apoptosis”, “FoxO signaling pathway”, “cell cycle”, “NF-kappa B signaling pathway”, “p53 signaling pathway”, “pentose phosphate pathway”, and “protein export”. For the cell–cell communication network analysis, the different interaction profiles between cell types were detected among primary, metastatic, and recurrent osteosarcomas. Further exploration of the KEGG pathway revealed that these ligand/receptor interactions may be associated with the NF-κB signaling pathway and its interacted mediators. In conclusion, the present study for the first time explored the scRNA-seq dataset in osteosarcoma, and our results revealed the 11 clustered cell types and demonstrated the novel cell–cell interactions among different cell types in primary, metastatic, and recurrent osteosarcomas. The NF-κB signaling pathway may play a key role in regulating the TME of osteosarcoma. The present study may provide new insights into understanding the molecular mechanisms of osteosarcoma pathophysiology.

Background: Epidemiological data associate high levels of combustion-derived particulate matter (PM) with deleterious respiratory outcomes, but the mechanism underlying those outcomes remains elusive. It has been acknowledged by the World Health Organization that PM exposure contributes to more than 4.2 million all-cause mortalities worldwide each year. Current literature demonstrates that PM exacerbates respiratory diseases, impairs lung function, results in chronic respiratory illnesses, and is associated with increased mortality. The proposed mechanisms revolve around oxidative stress and inflammation promoting pulmonary physiological remodeling. However, our previous data found that PM is capable of inducing T helper cell 17 (Th17) immune responses via aryl hydrocarbon receptor (Ahr) activation, which was associated with neutrophilic invasion characteristic of steroid insensitive asthma.Methods: In the present study, we utilized a combination of microarray and single cell RNA sequencing data to analyze the immunological landscape in mouse lungs following acute exposure to combustion derived particulate matter.Results: We present data that suggest epithelial cells produce specific cytokines in the aryl hydrocarbon receptor (Ahr) pathway that inform dendritic cells to initiate the production of pathogenic T helper (eTh17) cells. Using single-cell RNA sequencing analysis, we observed that upon exposure epithelial cells acquire a transcriptomic profile indicative of increased Il-17 signaling, Ahr activation, Egfr signaling, and T cell receptor and co-stimulatory signaling pathways. Epithelial cells further showed, Ahr activation is brought on by Ahr/ARNT nuclear translocation and activation of tyrosine kinase c-src, Egfr, and subsequently Erk1/2 pathways.Conclusions: Collectively, our data corroborates that PM initiates an eTh17 specific inflammatory response causing neutrophilic asthma through pathways in epithelial, dendritic, and T cells that promote eTh17 differentiation during initial PM exposure. [ABSTRACT FROM AUTHOR]

More and more researchers use single-cell RNA sequencing (scRNA-seq) technology to characterize the transcriptional map at the single-cell level. They use it to study the heterogeneity of complex tissues, transcriptome dynamics, and the diversity of unknown organisms. However, there are generally lots of technical and biological noises in the scRNA-seq data since the randomness of gene expression patterns. These data are often characterized by high-dimension, sparsity, large number of “dropout” values, and affected by batch effects. A large number of “dropout” values in scRNA-seq data seriously conceal the important relationship between genes and hinder the downstream analysis. Therefore, the imputation of dropout values of scRNA-seq data is particularly important. We classify, analyze and compare the current advanced scRNA-seq data imputation methods from different angles. Through the comparison and analysis of the principle, advantages and disadvantages of the algorithm, it can provide suggestions for the selection of imputation methods for specific problems and diverse data, and have basic research significance for the downstream function analysis of data.

It is well-established that anterior pituitary contains multiple endocrine cell populations, and each of them can secrete one/two hormone(s) to regulate vital physiological processes of vertebrates. However, the gene expression profiles of each pituitary cell population remains poorly characterized in most vertebrate groups. Here we analyzed the transcriptome of each cell population in adult chicken anterior pituitaries using single-cell RNA sequencing technology. The results showed that: (1) four out of five known endocrine cell clusters have been identified and designated as the lactotrophs, thyrotrophs, corticotrophs, and gonadotrophs, respectively. Somatotrophs were not analyzed in the current study. Each cell cluster can express at least one known endocrine hormone, and novel marker genes (e.g., CD24 and HSPB1 in lactotrophs, NPBWR2 and NDRG1 in corticotrophs; DIO2 and SOUL in thyrotrophs, C5H11ORF96 and HPGDS in gonadotrophs) are identified. Interestingly, gonadotrophs were shown to abundantly express five peptide hormones: FSH, LH, GRP, CART and RLN3; (2) four non-endocrine/secretory cell types, including endothelial cells (expressing IGFBP7 and CFD) and folliculo-stellate cells (FS-cells, expressing S100A6 and S100A10), were identified in chicken anterior pituitaries. Among them, FS-cells can express many growth factors, peptides (e.g., WNT5A, HBEGF, Activins, VEGFC, NPY, and BMP4), and progenitor/stem cell-associated genes (e.g., Notch signaling components, CDH1), implying that the FS-cell cluster may act as a paracrine/autocrine signaling center and enrich pituitary progenitor/stem cells; (3) sexually dimorphic expression of many genes were identified in most cell clusters, including gonadotrophs and lactotrophs. Taken together, our data provides a bird’s-eye view on the diverse aspects of anterior pituitaries, including cell composition, heterogeneity, cell-to-cell communication, and gene expression profiles, which facilitates our comprehensive understanding of vertebrate pituitary biology.

Alcohol-associated liver disease (ALD) has been recognized as the most common cause of advanced liver disease worldwide, though mechanisms of pathogenesis remain incompletely understood. The X-linked inhibitor of apoptosis (XIAP) protein was originally described as an anti-apoptotic protein that directly binds and inhibits caspases-3, 7, and 9. Here, we investigated the function of XIAP in hepatocytes in vitro using gain and loss-of-function approaches. We noted an XIAP-dependent increase in caspase activation as well as increased inflammatory markers and pro-inflammatory EV release from hepatocytes in vitro. Primary hepatocytes (PMH) from XiapAlb.Cre and XiaploxP mice exhibited higher cell death but surprisingly, lower expression of inflammation markers. Conditioned media from these isolated Xiap deleted PMH further decrease inflammation in bone marrow-derived macrophages. Also, interestingly, when administered an ethanol plus Fas-agonist-Jo2 model and an ethanol plus CCl4 model, these animals failed to develop an exacerbated disease phenotype in vivo. Of note, neither XiapAlb.Cre nor XiapAAV8.Cre mice presented with aggravated liver injury, hepatocyte apoptosis, liver steatosis, or fibrosis. Since therapeutics targeting XIAP are currently in clinical trials and caspase-induced death is very important for development of ALD, we sought to explore the potential basis of this unexpected lack of effect. We utilized scRNA-seq and spatially reconstructed hepatocyte transcriptome data from human liver tissue and observed that XIAP was significantly zonated, along with its endogenous inhibitor second mitochondria-derived activator of caspases (SMAC) in periportal region. This contrasted with pericentral zonation of other IAPs including cIAP1 and Apollon as well as caspases 3, 7, and 9. Thus providing a potential explanation for compensation of the effect of Xiap deletion by other IAPs. In conclusion, our findings implicate a potential zonallydependent role for SMAC that prevented development of a phenotype in XIAP knockout mice in ALD models. Targeting SMAC may also be important in addition to current efforts of targeting XIAP in treatment of ALD.

Single cell RNA-seq data, like data from other sequencing technology, contain systematic technical noise. Such noise results from a combined effect of unequal efficiencies in the capturing and counting of mRNA molecules, such as extraction/amplification efficiency and sequencing depth. We show that such technical effects are not only cell-specific, but also affect genes differently, thus a simple cell-wise size factor adjustment may not be sufficient. We present a non-linear normalization approach that provides a cell- and gene-specific normalization factor for each gene in each cell. We show that the proposed normalization method (implemented in “SC2P" package) reduces more technical variation than competing methods, without reducing biological variation. When technical effects such as sequencing depths are not balanced between cell populations, SC2P normalization also removes the bias due to uneven technical noise. This method is applicable to scRNA-seq experiments that do not use unique molecular identifier (UMI) thus retain amplification biases.

The importance of T cells in controlling SARS-CoV-2 infections has been demonstrated widely, but insights into the quality of these responses are still limited due to technical challenges. Indeed, understanding the functionality of the T-cell receptor (TCR) repertoire of a polyclonal antigen-specific population still requires the tedious work of T-cell cloning or TCR re-expression and subsequent characterization. In this work, we show that it is possible to discriminate highly functional and bystander TCRs based on gene signatures of T-cell activation induced by recent peptide stimulation. SARS-CoV-2-specific TCRs previously identified by cytokine release after peptide restimulation and subsequent single-cell RNA sequencing were re-expressed via CRISPR-Cas9-mediated gene editing into a Jurkat-based reporter cell line system suitable for high-throughput screening. We could observe differences in SARS-CoV-2 epitope recognition as well as a wide range of functional avidities. By correlating these in vitro TCR engineered functional data with the transcriptomic profiles of the corresponding TCR-expressing parental T cells, we could validate that gene signatures of recent T-cell activation accurately identify and predict truly SARS-CoV-2-specific TCRs. In summary, this work paves the way for alternative approaches useful for the functional analysis of global antigen-specific TCR repertoires with largely improved throughput.

It is well-established that anterior pituitary contains multiple endocrine cell populations, and each of them can secrete one/two hormone(s) to regulate vital physiological processes of vertebrates. However, the gene expression profiles of each pituitary cell population remains poorly characterized in most vertebrate groups. Here we analyzed the transcriptome of each cell population in adult chicken anterior pituitaries using single-cell RNA sequencing technology. The results showed that: (1) four out of five known endocrine cell clusters have been identified and designated as the lactotrophs, thyrotrophs, corticotrophs, and gonadotrophs, respectively. Somatotrophs were not analyzed in the current study. Each cell cluster can express at least one known endocrine hormone, and novel marker genes (e.g., CD24 and HSPB1 in lactotrophs, NPBWR2 and NDRG1 in corticotrophs; DIO2 and SOUL in thyrotrophs, C5H11ORF96 and HPGDS in gonadotrophs) are identified. Interestingly, gonadotrophs were shown to abundantly express five peptide hormones: FSH, LH, GRP, CART and RLN3 ; (2) four non-endocrine/secretory cell types, including endothelial cells (expressing IGFBP7 and CFD) and folliculo-stellate cells (FS-cells, expressing S100A6 and S100A10), were identified in chicken anterior pituitaries. Among them, FS-cells can express many growth factors, peptides (e.g., WNT5A, HBEGF , Activins, VEGFC, NPY , and BMP4), and progenitor/stem cell-associated genes (e.g. , Notch signaling components, CDH1), implying that the FS-cell cluster may act as a paracrine/autocrine signaling center and enrich pituitary progenitor/stem cells; (3) sexually dimorphic expression of many genes were identified in most cell clusters, including gonadotrophs and lactotrophs. Taken together, our data provides a bird's-eye view on the diverse aspects of anterior pituitaries, including cell composition, heterogeneity, cell-to-cell communication, and gene expression profiles, which facilitates our comprehensive understanding of vertebrate pituitary biology. [ABSTRACT FROM AUTHOR]

Alcohol-associated liver disease (ALD) has been recognized as the most common cause of advanced liver disease worldwide, though mechanisms of pathogenesis remain incompletely understood. The X-linked inhibitor of apoptosis (XIAP) protein was originally described as an anti-apoptotic protein that directly binds and inhibits caspases-3, 7, and 9. Here, we investigated the function of XIAP in hepatocytes in vitro using gain and loss-of-function approaches. We noted an XIAP-dependent increase in caspase activation as well as increased inflammatory markers and pro-inflammatory EV release from hepatocytes in vitro. Primary hepatocytes (PMH) from Xiap Alb.Cre and Xiap loxP mice exhibited higher cell death but surprisingly, lower expression of inflammation markers. Conditioned media from these isolated Xiap deleted PMH further decrease inflammation in bone marrow-derived macrophages. Also, interestingly, when administered an ethanol plus Fas-agonist-Jo2 model and an ethanol plus CCl4 model, these animals failed to develop an exacerbated disease phenotype in vivo. Of note, neither Xiap Alb . Cre nor Xiap AAV8.Cre mice presented with aggravated liver injury, hepatocyte apoptosis, liver steatosis, or fibrosis. Since therapeutics targeting XIAP are currently in clinical trials and caspase-induced death is very important for development of ALD, we sought to explore the potential basis of this unexpected lack of effect. We utilized scRNA-seq and spatially reconstructed hepatocyte transcriptome data from human liver tissue and observed that XIAP was significantly zonated, along with its endogenous inhibitor second mitochondria-derived activator of caspases (SMAC) in periportal region. This contrasted with pericentral zonation of other IAPs including cIAP1 and Apollon as well as caspases 3, 7, and 9. Thus providing a potential explanation for compensation of the effect of Xiap deletion by other IAPs. In conclusion, our findings implicate a potential zonallydependent role for SMAC that prevented development of a phenotype in XIAP knockout mice in ALD models. Targeting SMAC may also be important in addition to current efforts of targeting XIAP in treatment of ALD. [ABSTRACT FROM AUTHOR]

Single cell RNA-seq data, like data from other sequencing technology, contain systematic technical noise. Such noise results from a combined effect of unequal efficiencies in the capturing and counting of mRNA molecules, such as extraction/amplification efficiency and sequencing depth. We show that such technical effects are not only cell-specific, but also affect genes differently, thus a simple cell-wise size factor adjustment may not be sufficient. We present a non-linear normalization approach that provides a cell- and gene-specific normalization factor for each gene in each cell. We show that the proposed normalization method (implemented in "SC2P" package) reduces more technical variation than competing methods, without reducing biological variation. When technical effects such as sequencing depths are not balanced between cell populations, SC2P normalization also removes the bias due to uneven technical noise. This method is applicable to scRNA-seq experiments that do not use unique molecular identifier (UMI) thus retain amplification biases. [ABSTRACT FROM AUTHOR]

Coronavirus disease‐2019 (COVID‐19) is a global pandemic and caused by severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2), which has resulted in millions of deaths worldwide. Reports denote SARS‐CoV‐2 uses angiotensin‐converting enzyme 2 (ACE2), transmembrane serine protease 2 (TMPRSS2) as its primary entry point into the host cell. However, understanding the biology behind this viral replication, disease mechanism and drug discovery efforts are limited due to the lack of a suitable experimental model. Here, we used single‐cell RNA sequencing data of human organoids to analyze expressions of ACE2 and TMPRSS2, in addition to an array of RNA receptors to examine their role in SARS‐CoV‐2 pathogenesis. ACE2 is abundant in all organoids, except the prostate and brain, and TMPRSS2 is omnipresent. Innate immune pathways are upregulated in ACE2(+) cells of all organoids, except the lungs. Besides this, the expression of low‐density lipoprotein receptor is highly enriched in ACE2(+) cells in intestinal, lung, and retinal organoids, with the highest expression in lung organoids. Collectively, this study demonstrates that the organoids can be used as an experimental platform to explore this novel virus disease mechanism and for drug development. [ABSTRACT FROM AUTHOR]

Idiopathic pulmonary fibrosis (IPF) is a progressive and fatal lung disease with limited therapeutic options. The current model suggests that chronic or repetitive “micro-injuries” of the alveolar epithelium lead to activation and proliferation of fibroblasts and excessive extracellular matrix (ECM) deposition. Disruption of alveolar type II (ATII) epithelial cell homeostasis and the characteristics of mesenchymal cell populations in IPF have received particular attention in recent years. Emerging data from single cell RNA sequencing (scRNAseq) analysis shed novel light on alterations in ATII cell progenitor dysfunction and the diversity of mesenchymal cells within the fibrotic lung. Within this minireview, we summarize the data from most recent human scRNAseq studies. We aim to collate the current knowledge on cellular plasticity and heterogeneity in the development and progression of IPF, effects of drug treatment on transcriptional changes. Finally, we provide a brief outlook on future challenges and promises for large scale sequencing studies in the development of novel therapeutics for IPF.

Sciatic nerve crush injury triggers sterile inflammation within the distal nerve and axotomized dorsal root ganglia (DRGs). Granulocytes and pro-inflammatory Ly6Chigh monocytes infiltrate the nerve first and rapidly give way to Ly6Cnegative inflammation-resolving macrophages. In axotomized DRGs, few hematogenous leukocytes are detected and resident macrophages acquire a ramified morphology. Single-cell RNA-sequencing of injured sciatic nerve identifies five macrophage subpopulations, repair Schwann cells, and mesenchymal precursor cells. Macrophages at the nerve crush site are molecularly distinct from macrophages associated with Wallerian degeneration. In the injured nerve, macrophages ‘eat’ apoptotic leukocytes, a process called efferocytosis, and thereby promote an anti-inflammatory milieu. Myeloid cells in the injured nerve, but not axotomized DRGs, strongly express receptors for the cytokine GM-CSF. In GM-CSF-deficient (Csf2-/-) mice, inflammation resolution is delayed and conditioning-lesion-induced regeneration of DRG neuron central axons is abolished. Thus, carefully orchestrated inflammation resolution in the nerve is required for conditioning-lesion-induced neurorepair.

Myeloid cell heterogeneity remains poorly studied in breast cancer, and particularly in premalignancy. Here, we used single cell RNA sequencing to characterize macrophage diversity in mouse pre-invasive lesions as compared to lesions undergoing localized invasion. Several subpopulations of macrophages with transcriptionally distinct profiles were identified, two of which resembled macrophages in the steady state. While all subpopulations expressed tumor-promoting genes, many of the populations expressed pro-inflammatory genes, differing from reports in tumor-associated macrophages. Gene profiles of the myeloid cells were similar between early and late stages of premalignancy, although expansion of some subpopulations occurred. These results unravel macrophage heterogeneity in early progression and may provide insight into early intervention strategies that target macrophages.

Myeloid cell heterogeneity remains poorly studied in breast cancer, and particularly in premalignancy. Here, we used single cell RNA sequencing to characterize macrophage diversity in mouse pre-invasive lesions as compared to lesions undergoing localized invasion. Several subpopulations of macrophages with transcriptionally distinct profiles were identified, two of which resembled macrophages in the steady state. While all subpopulations expressed tumor-promoting genes, many of the populations expressed pro-inflammatory genes, differing from reports in tumor-associated macrophages. Gene profiles of the myeloid cells were similar between early and late stages of premalignancy, although expansion of some subpopulations occurred. These results unravel macrophage heterogeneity in early progression and may provide insight into early intervention strategies that target macrophages. [ABSTRACT FROM AUTHOR]

Single-cell RNAsequencing (scRNA-seq) technologies have enabled the large-scale whole-transcriptome profiling of each individual single cell in a cell population. A core analysis of the scRNA-seq transcriptome profiles is to cluster the single cells to reveal cell subtypes and infer cell lineages based on the relations among the cells. This article reviews the machine learning and statistical methods for clustering scRNA-seq transcriptomes developed in the past few years. The review focuses on how conventional clustering techniques such as hierarchical clustering, graph-based clustering, mixture models, |$k$| -means, ensemble learning, neural networks and density-based clustering are modified or customized to tackle the unique challenges in scRNA-seq data analysis, such as the dropout of low-expression genes, low and uneven read coverage of transcripts, highly variable total mRNAs from single cells and ambiguous cell markers in the presence of technical biases and irrelevant confounding biological variations. We review how cell-specific normalization, the imputation of dropouts and dimension reduction methods can be applied with new statistical or optimization strategies to improve the clustering of single cells. We will also introduce those more advanced approaches to cluster scRNA-seq transcriptomes in time series data and multiple cell populations and to detect rare cell types. Several software packages developed to support the cluster analysis of scRNA-seq data are also reviewed and experimentally compared to evaluate their performance and efficiency. Finally, we conclude with useful observations and possible future directions in scRNA-seq data analytics. Availability All the source code and data are available at https://github.com/kuanglab/single-cell-review. [ABSTRACT FROM AUTHOR]

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., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Background: 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., Methods: 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)., Results: 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., Conclusions: 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., (© 2023 Wiley Periodicals LLC.)