UNASSIGNED: Metabolic liver diseases, including alcohol- and non-alcoholic fatty liver diseases (ALD/NAFLD), are characterized by inflammation and decreased ability to prevent infections. Patients with severe alcohol-associated hepatitis (sAH) are particularly susceptible to infections while undergoing treatment with steroids. Understanding the immunological mechanisms for these responses is critical to managing the treatment of patients with metabolic liver diseases. Cytotoxic NK cells and CD8 T cells, using cytolytic granules, serve an important immunological role by killing infected cells, including monocytes. However, patients with sAH have dysfunctional NK cells, which cannot kill target cells, though the mechanism is unknown.
UNASSIGNED: We performed an exploratory study using single-cell RNA-seq (scRNA-seq) (n = 4) and multi-panel intracellular flow cytometry (n = 7-8 for all patient groups) on PBMCs isolated from patients with sAH and healthy controls (HC).
UNASSIGNED: ScRNA-seq revealed receptors in NK cells and CD8 T cells required for cytotoxic cell recognition of activated monocytes were downregulated in patients with sAH compared to healthy controls. Granulysin was the most downregulated gene in both NK cells and effector CD8 T cells. In NK cells from HC, expression of granulysin, perforin, and granzymes A and B was highly correlated; however, in sAH, these genes lost coordinate expression, indicative of dysfunctional cytolytic granule formation. Finally, the expression of cytolytic granule proteins in NK cells was decreased from sAH, indicating reduced cytolytic granules.
UNASSIGNED: Together, these results suggest a loss of cytotoxic cell function in PBMCs from sAH that may contribute to a decreased ability to communicate with other immune cells, such as monocytes, and prevent the killing of infected cells, thus increasing the risk of infection.

Our understanding of human fetal cerebellum development during the late second trimester, a critical period for the generation of astrocytes, oligodendrocytes, and unipolar brush cells (UBCs), remains limited. Here, we performed single-cell RNA sequencing (scRNA-seq) in human fetal cerebellum samples from gestational weeks (GWs) 18-25. We find that proliferating UBC progenitors distribute in the subventricular zone of the rhombic lip (RLSVZ) near white matter (WM), forming a layer structure. We also delineate two trajectories from astrogenic radial glia (ARGs) to Bergmann glial progenitors (BGPs) and recognize oligodendrogenic radial glia (ORGs) as one source of primitive oligodendrocyte progenitor cells (PriOPCs). Additionally, our scRNA-seq analysis of the trisomy 21 fetal cerebellum at this stage reveals abnormal upregulated genes in pathways such as the cell adhesion pathway and focal adhesion pathway, which potentially promote neuronal differentiation. Overall, our research provides valuable insights into normal and abnormal development of the human fetal cerebellum.

OBJECTIVE: Osteosarcoma, a highly malignant primary bone tumor primarily affecting adolescents, frequently develops resistance to initial chemotherapy, leading to metastasis and limited treatment options. Our study aims to uncover novel therapeutic targets for metastatic and recurrent osteosarcoma.
METHODS: In this study, we proved the potential of modulating the YAP1-regulated glutamine metabolic pathway to augment the response of OS to DFMO. We initially employed single-cell transcriptomic data to gauge the activation level of polyamine metabolism in MTAP-deleted OS patients. This was further substantiated by transcriptome sequencing data from recurrent and non-recurrent patient tissues, confirming the activation of polyamine metabolism in progressive OS. Through high-throughput drug screening, we pinpointed CIL56, a YAP1 inhibitor, as a promising candidate for a combined therapeutic strategy with DFMO. In vivo, we utilized PDX and CDX models to validate the therapeutic efficacy of this drug combination. In vitro, we conducted western blot analysis, qPCR analysis, immunofluorescence staining, and PuMA experiments to monitor alterations in molecular expression, distribution, and tumor metastasis capability. We employed CCK-8 and colony formation assays to assess the proliferative capacity of cells in the experimental group. We used flow cytometry and reactive oxygen probes to observe changes in ROS and glutamine metabolism within the cells. Finally, we applied RNA-seq in tandem with metabolomics to identify metabolic alterations in OS cells treated with a DFMO and CIL56 combination. This enabled us to intervene and validate the role of the YAP1-mediated glutamine metabolic pathway in DFMO resistance.
RESULTS: Through single-cell RNA-seq data analysis, we pinpointed a subset of late-stage OS cells with significantly upregulated polyamine metabolism. This upregulation was further substantiated by transcriptomic profiling of recurrent and non-recurrent OS tissues. High-throughput drug screening revealed a promising combination strategy involving DFMO and CIL56. DFMO treatment curbs the phosphorylation of YAP1 protein in OS cells, promoting nuclear entry and initiating the YAP1-mediated glutamine metabolic pathway. This reduces intracellular ROS levels, countering DFMO\'s anticancer effect. The therapeutic efficacy of DFMO can be amplified both in vivo and in vitro by combining it with the YAP1 inhibitor CIL56 or the glutaminase inhibitor CB-839. This underscores the significant potential of targeting the YAP1-mediated glutamine metabolic pathway to enhance efficacy of DFMO.
CONCLUSIONS: Our findings elucidate YAP1-mediated glutamine metabolism as a crucial bypass mechanism against DFMO, following the inhibition of polyamine metabolism. Our study provides valuable insights into the potential role of DFMO in an \"One-two Punch\" therapy of metastatic and recurrent osteosarcoma.

BACKGROUND: Pancreatic fibrosis is an early diagnostic feature of the common inherited disorder cystic fibrosis (CF). Many people with CF (pwCF) are pancreatic insufficient from birth and the replacement of acinar tissue with cystic lesions and fibrosis is a progressive phenotype that may later lead to diabetes. Little is known about the initiating events in the fibrotic process though it may be a sequela of inflammation in the pancreatic ducts resulting from loss of CFTR impairing normal fluid secretion. Here we use a sheep model of CF (CFTR-/-) to examine the evolution of pancreatic disease through gestation.
METHODS: Fetal pancreas was collected at six time points from 50-days of gestation through to term, which is equivalent to ~ 13 weeks to term in human. RNA was extracted from tissue for bulk RNA-seq and single cells were prepared from 80-day, 120-day and term samples for scRNA-seq. Data were validated by immunochemistry.
RESULTS: Transcriptomic evidence from bulk RNA-seq showed alterations in the CFTR-/- pancreas by 65-days of gestation, which are accompanied by marked pathological changes by 80-days of gestation. These include a fibrotic response, confirmed by immunostaining for COL1A1, αSMA and SPARC, together with acinar loss. Moreover, using scRNA-seq we identify a unique cell population that is significantly overrepresented in the CFTR-/- animals at 80- and 120-days gestation, as are stellate cells at term.
CONCLUSIONS: The transcriptomic changes and cellular imbalance that we observe likely have pivotal roles in the evolution of CF pancreatic disease and may provide therapeutic opportunities to delay or prevent pancreatic destruction in CF.

We describe a time-resolved nascent single-cell RNA sequencing (RNA-seq) approach that measures gene-specific transcriptional noise and the fraction of active genes in S. cerevisiae. Most genes are expressed with near-constitutive behavior, while a subset of genes show high mRNA variance suggestive of transcription bursting. Transcriptional noise is highest in the cofactor/coactivator-redundant (CR) gene class (dependent on both SAGA and TFIID) and strongest in TATA-containing CR genes. Using this approach, we also find that histone gene transcription switches from a low-level, low-noise constitutive mode during M and M/G1 to an activated state in S phase that shows both an increase in the fraction of active promoters and a switch to a noisy and bursty transcription mode. Rapid depletion of cofactors SAGA and MED Tail indicates that both factors play an important role in stimulating the fraction of active promoters at CR genes, with a more modest role in transcriptional noise.

Abdominal aortic aneurysm (AAA) is characterized by permanent luminal expansion and a high mortality rate due to aortic rupture. Despite the identification of abnormalities in the mevalonate pathway (MVA) in many diseases, including cardiovascular diseases, the potential impact of this pathway on AAA remains unclear. This study aims to investigate whether the expression of the MVA-related enzyme is altered during the progression of angiotensin II (Ang II) -induced AAA.Ang II 28D and Ang II 5D groups were continuously perfused with Ang II for 28 days and 5 days, respectively, and the Sham group was perfused with saline. The general and remodeling characteristics of AAA were determined by biochemical and histological analysis. Alteration of MVA-related enzyme expressions was revealed by western blot and single-cell RNA sequencing (scRNA-seq).The continuous Ang II infusion for 28 days showed significant aorta expansion and arterial remodeling. Although the arterial diameter slightly increased, the aneurysm formation was not found in Ang II induction for 5 days. MVA-related enzyme expression and activation of small GTP-binding proteins were significantly increased after Ang II-induced. As verified by scRNA-seq, the key enzyme gene expression was also higher in Ang II 28D. Similarly, it was detected that the expression levels of the above enzymes and the activity of small G proteins were elevated in the early stage of AAA as induced by Ang II infusion for 5 days.Continuous Ang II infusion-induced abdominal aortic expansion and arterial remodeling were accompanied by altered expression of key enzymes in the MVA.

BACKGROUND: Breast cancer (BC) is a heterogeneous disease, with the ductal subtype exhibiting significant cellular diversity that influences prognosis and response to treatment. Single-cell RNA sequencing data from the GEO database were utilized in this study to investigate the underlying mechanisms of cellular heterogeneity and to identify potential prognostic markers and therapeutic targets.
METHODS: Bioinformatics analysis was conducted using R packages to analyze the single-cell sequencing data. The presence of highly variable genes and differences in malignant potency within the same BC samples were examined. Differential gene expression and biological function between Type 1 and Type 2 ductal epithelial cells were identified. Lasso regression and Cox proportional hazards regression analyses were employed to identify genes associated with patient prognosis. Experimental validation was performed in vitro and in vivo to confirm the functional relevance of the identified genes.
RESULTS: The analysis revealed notable heterogeneity among BC cells, with the presence of highly variable genes and differences in malignant behavior within the same samples. Significant disparities in gene expression and biological function were identified between Type 1 and Type 2 ductal epithelial cells. Through regression analyses, CYP24A1 and TFPI2 were identified as pivotal genes associated with patient prognosis. Kaplan-Meier curves demonstrated their prognostic significance, and experimental validation confirmed their inhibitory effects on malignant behaviors of ductal BC cells.
CONCLUSIONS: This study highlights the cellular heterogeneity in ductal subtype breast cancer and delineates the differential gene expressions and biological functions between Type 1 and Type 2 ductal epithelial cells. The genes CYP24A1 and TFPI2 emerged as promising prognostic markers and therapeutic targets, exhibiting inhibitory effects on BC cell malignancy in vitro and in vivo. These findings offer the potential for improved BC management and the development of targeted treatment strategies.

Our laboratory is interested in investigating the maturation process of zebrafish thrombocytes, which are functional equivalents to human platelets. We have adopted the zebrafish model to gain insights into mammalian platelet production, or thrombopoiesis. Notably, zebrafish exhibit two distinct populations of thrombocytes in their circulating blood: young and mature thrombocytes. This observation is intriguing because maturation appears to occur in circulation, yet the precise mechanisms governing this maturation remain elusive. Our goal is to understand the mechanisms underlying thrombocyte maturation by conducting single-cell RNA sequencing (scRNA-Seq) on young and mature thrombocytes, analyzing these transcriptomes to identify genes specific to each thrombocyte population, and elucidating the role of these genes in the maturation process, by quantifying thrombocyte numbers after the piggyback knockdown of each of these genes. In this chapter, we present a comprehensive, step-by-step protocol detailing the multifaceted methodology involved in understanding thrombocyte maturation, which encompasses the collection of zebrafish blood, the separation of young and mature thrombocytes using flow cytometry, scRNA-Seq analysis of these distinct thrombocyte populations, identification of genes specific to young and mature thrombocytes, and subsequent validation through gene knockdown techniques.

Single-cell transcriptomics allows unbiased characterization of cell heterogeneity in a sample by profiling gene expression at single-cell level. These profiles capture snapshots of transient or steady states in dynamic processes, such as cell cycle, activation, or differentiation, which can be computationally ordered into a \"flip-book\" of cell development using trajectory inference methods. However, prediction of more complex topology structures, such as multifurcations or trees, remains challenging. In this chapter, we present two user-friendly protocols for inferring tree-shaped single-cell trajectories and pseudotime from single-cell transcriptomics data with Totem. Totem is a trajectory inference method that offers flexibility in inferring both nonlinear and linear trajectories and usability by avoiding the cumbersome fine-tuning of parameters. The QuickStart protocol provides a simple and practical example, whereas the GuidedStart protocol details the analysis step-by-step. Both protocols are demonstrated using a case study of human bone marrow CD34+ cells, allowing the study of the branching of three lineages: erythroid, lymphoid, and myeloid. All the analyses can be fully reproduced in Linux, macOS, and Windows operating systems (amd64 architecture) with >8 Gb of RAM using the provided docker image distributed with notebooks, scripts, and data in Docker Hub (elolab/repro-totem-ti). These materials are shared online under open-source license at https://elolab.github.io/Totem-protocol .

This chapter shows applying the Asymmetric Within-Sample Transformation to single-cell RNA-Seq data matched with a previous dropout imputation. The asymmetric transformation is a special winsorization that flattens low-expressed intensities and preserves highly expressed gene levels. Before a standard hierarchical clustering algorithm, an intermediate step removes noninformative genes according to a threshold applied to a per-gene entropy estimate. Following the clustering, a time-intensive algorithm is shown to uncover the molecular features associated with each cluster. This step implements a resampling algorithm to generate a random baseline to measure up/downregulated significant genes. To this aim, we adopt a GLM model as implemented in DESeq2 package. We render the results in graphical mode. While the tools are standard heat maps, we introduce some data scaling to clarify the results\' reliability.