Single-cell RNA transcriptomic analyses of tumor microenvironment of ovarian metastasis in gastric cancer.

PURPOSE: Ovarian metastasis of gastric cancer (GC), commonly referred to as Krukenberg tumors, leads to a poor prognosis. However, the cause of metastasis remains unknown. Here, we present an integrated single-cell RNA sequencing (scRNA-Seq) analysis of the immunological microenvironment of two paired clinical specimens with ovarian metastasis of GC. METHODS: scRNA-Seq was performed to determine the immunological microenvironment in ovarian metastasis of gastric cancer. CellChat was employed to analyze cell-cell communications across different cell types. Functional enrichment analysis was done by enrichKEGG in clusterProfiler. GEPIA2 was used to assess the influence of certain genes and gene signatures on prognosis. RESULTS: The ovarian metastasis tissues exhibit a heterogenous immunological microenvironment compared to the primary tumors. Exhaustion of T and B cells is observed in the ovarian metastasis tissues. Compared to the paired adjacent non-tumoral and primary tumors, the ratio of endothelial cells and fibroblasts is high in the ovarian metastasis tissues. Compared to primary ovarian cancers, we identify a specific group of tumor-associated fibroblasts with MFAP4 and CAPNS1 expression in the ovarian metastatic tissues of GC. We further define metastasis-related-endothelial and metastasis-related-fibroblast signatures and indicate that patients with these high signature scores have a poor prognosis. In addition, the ovarian metastasis tissue has a lower level of intercellular communications compared to the primary tumor. CONCLUSION: Our findings reveal the immunological microenvironment of ovarian metastasis of gastric cancer and will promote the discovery of new therapeutic strategies for ovarian metastasis in gastric cancer.

Multimodal profiling of biostabilized human skin modules reveals a coordinated ecosystem response to injected mRNA-1273 COVID-19 vaccine.

BACKGROUND: The field of drug development is witnessing a remarkable surge in the development of innovative strategies. There is a need to develop technological platforms capable of generating human data prior to progressing to clinical trials. METHODS: Here we introduce a new flexible solution designed for the comprehensive monitoring of the natural human skin ecosystem's response to immunogenic drugs over time. Based on unique bioengineering to preserve surgical resections in a long survival state, it allows for the first time a comprehensive analysis of resident immune cells response at both organ and single-cell levels. RESULTS: Upon injection of the mRNA-1273 COVID-19 vaccine, we characterized precise sequential molecular events triggered upon detection of the exogenous substance. The vaccine consistently targets DC/macrophages and mast cells, regardless of the administration route, while promoting specific cell-cell communications in surrounding immune cell subsets. CONCLUSION: Given its direct translational relevance, this approach provides a multiscale vision of genuine human tissue immunity that could pave the way toward the development of new vaccination and drug development strategies.

Identification of a Novel ECM Remodeling Macrophage Subset in AKI to CKD Transition by Integrative Spatial and Single-Cell Analysis.

The transition from acute kidney injury (AKI) to chronic kidney disease (CKD) is a critical clinical issue. Although previous studies have suggested macrophages as a key player in promoting inflammation and fibrosis during this transition, the heterogeneity and dynamic characterization of macrophages are still poorly understood. Here, we used integrated single-cell RNA sequencing and spatial transcriptomic to characterize the spatiotemporal heterogeneity of macrophages in murine AKI-to-CKD model of unilateral ischemia-reperfusion injury. A marked increase in macrophage infiltration at day 1 was followed by a second peak at day 14 post AKI. Spatiotemporal profiling revealed that injured tubules and macrophages co-localized early after AKI, whereas in late chronic stages had spatial proximity to fibroblasts. Further pseudotime analysis revealed two distinct lineages of macrophages in this transition: renal resident macrophages differentiated into the pro-repair subsets, whereas infiltrating monocyte-derived macrophages contributed to chronic inflammation and fibrosis. A novel macrophage subset, extracellular matrix remodeling-associated macrophages (EAMs) originating from monocytes, linked to renal fibrogenesis and communicated with fibroblasts via insulin-like growth factors (IGF) signalling. In sum, our study identified the spatiotemporal dynamics of macrophage heterogeneity with a unique subset of EAMs in AKI-to-CKD transition, which could be a potential therapeutic target for preventing CKD development.

Large-scale cross-trait genetic analysis highlights shared genetic backgrounds of autoimmune diseases.

Disorders associated with the immune system burden multiple organs, although the shared biology exists across the diseases. Preceding family-based studies reveal that immune diseases are heritable to varying degrees, providing the basis for immunogenomics. The recent cost reduction in genetic analysis intensively promotes biobank-scale studies and the development of frameworks for statistical genetics. The accumulating multi-layer omics data, including genome-wide association studies (GWAS) and RNA-sequencing at single-cell resolution, enable us to dissect the genetic backgrounds of immune-related disorders. Although autoimmune and allergic diseases are generally categorized into different disease categories, epidemiological studies reveal the high incidence of autoimmune and allergic disease complications, suggesting the shared genetics and biology between the disease categories. Biobank resources and consortia cover multiple immune-related disorders to accumulate phenome-wide associations of genetic variants and enhance researchers to analyze the shared and heterogeneous genetic backgrounds. The emerging post-GWAS and integrative multi-omics analyses provide genetic and biological insights into the multicategorical disease associations.

Hepatic BMAL1 and HIF1α regulate a time-dependent hypoxic response and prevent hepatopulmonary-like syndrome.

The transcriptional response to hypoxia is temporally regulated, yet the molecular underpinnings and physiological implications are unknown. We examined the roles of hepatic Bmal1 and Hif1α in the circadian response to hypoxia in mice. We found that the majority of the transcriptional response to hypoxia is dependent on either Bmal1 or Hif1α, through shared and distinct roles that are daytime determined. We further show that hypoxia-inducible factor (HIF)1α accumulation upon hypoxia is temporally regulated and Bmal1 dependent. Unexpectedly, mice lacking both hepatic Bmal1 and Hif1α are hypoxemic and exhibit increased mortality upon hypoxic exposure in a daytime-dependent manner. These mice display mild liver dysfunction with pulmonary vasodilation likely due to extracellular signaling regulated kinase (ERK) activation, endothelial nitric oxide synthase, and nitric oxide accumulation in lungs, suggestive of hepatopulmonary syndrome. Our findings indicate that hepatic BMAL1 and HIF1α are key time-dependent regulators of the hypoxic response and can provide molecular insights into the pathophysiology of hepatopulmonary syndrome.

A human fetal cerebellar map of the late second trimester reveals developmental molecular characteristics and abnormality in trisomy 21.

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.

Transcriptional noise, gene activation, and roles of SAGA and Mediator Tail measured using nucleotide recoding single-cell RNA-seq.

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.

A Low Expression of NRF2 Enhances Oxidative Stress and Autophagy in Myofibroblasts, Promoting Progression of Chronic Obstructive Pulmonary Disease.

BACKGROUND: The inflammation phenotypes are often closely related to oxidative stress and autophagy pathway activation, which could be developed as a treatment target. AIMS: The aim of this study was to explore the underlying mechanism of inflammation in chronic obstructive pulmonary disease (COPD). METHODS: The lung tissue single-cell RNA-seq (scRNA-seq) dataset of GSE171541 was downloaded from the Gene Expression Omnibus (GEO) database. The marker genes were obtained from the CellMarker database. "Seurat" and "harmony" R packages were used for single-cell profiling analysis. Then, the "AUCell" R package was employed to calculate the reactive oxygen species (ROS) and autophagy pathway scores. Gene regulation network analysis was performed by applying the "SCENIC" package, followed by conducting correlation analysis with Spearman's rank correlation method. The cigarettes were used to develop a traumatic model in mice, and the expression of relevant genes was measured by qRT-PCR. RESULTS: The scRNA-seq analysis classified 12 cell subgroups in which the contractility of myofibroblasts was closely associated with the progression of COPD. Further analysis showed that ROS and autophagy pathways were significantly activated in myofibroblasts and that the nuclear factor erythroid 2-related factor 2 (NRF2) and its mediated oxidative stress pathway were inhibited in myofibroblasts. In addition, the downregulated NRF2 gene was negatively correlated with the expression of autophagy and ROS activation. In the traumatic mice model, NRF2 was downregulated in COPD mice but further elevated in the COPD+NRF2 mice group. Interestingly, the mRNA levels of Kelchlike ECH-associated protein 1 (Keap1), NADPH oxidase (NOX), and Cathepsin B (CTSB) were upregulated in COPD group in comparison to the control group but they were downregulated by NRF2. These results suggested that low-expressed NFR2 promoted autophagy and ROS pathway activation in myofibroblasts for COPD progression. CONCLUSION: We identified a cell myofibroblast cluster closely associated with COPD progression using the scRNA-seq analysis. The downregulated NFR2, as a key risk factor, mediated myofibroblast death by activating the oxidative stress and autophagy pathway for COPD progression.

Stellate cells are in utero markers of pancreatic disease in cystic fibrosis.

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. CONCLUSION: 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.

Disrupting YAP1-mediated glutamine metabolism induces synthetic lethality alongside ODC1 inhibition in osteosarcoma.

PURPOSE: 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. CONCLUSION: 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.

Single-cell RNA-seq reveals heterogeneity in metastatic renal cell carcinoma and effect of anti-angiogenesis therapy in the pancreas metastatic lesion.

Metastatic clear cell renal cell carcinoma has heterogenous tumor microenvironment (TME). Among the metastatic lesions, pancreas metastasis is rare and controversy in treatment approaches. Here, extensive primary and metastatic lesion samples were included by single-cell RNA-seq to decipher the distinct metastasis TME. The hypoxic and inflammatory TME of pancreas metastasis was decoded in this study, and the activation of PAX8-myc signaling, and metabolic reprogramming were observed. The active components including endothelial cells, fibroblasts and T cells were profiled. Meanwhile, we also evaluated the effect of anti-angiogenesis treatment in the pancreas metastasis patient. The potential mechanisms of pancreatic tropism, instability of genome, and the response of immunotherapy were also discussed in this work. Taken together, our findings suggest a clue to the heterogeneity in metastasis TME and provide evidence for the treatment of pancreas metastasis in renal cell carcinoma patients.

Combining Bulk and Single Cell RNA-Sequencing Data to Identify Hub Genes of Fibroblasts in Dilated Cardiomyopathy.

Background: Dilated cardiomyopathy (DCM) is the second leading cause of heart failure, with intricate pathophysiological underpinnings. In order to shed fresh light on the mechanistic research of DCM, we combined bulk RNA-seq and single-cell RNA-seq (scRNA-seq) data to examine significant cells and genes implicated in the disease. Methods: This analysis employed publicly accessible bulk RNA-seq and scRNA-seq DCM datasets. The scRNA-seq data underwent normalization, principal component, and t-distribution stochastic neighbor embedding analysis. Cell-to-cell communication networks and activity analysis were conducted using CellChat. Utilizing enrichment analysis, the marker genes' role in the active cells was evaluated. After screening by limma software and weighted gene co-expression network analysis, the differentially expressed genes (DEGs) served as hub genes. Furthermore, these hub genes were subjected to immunological studies, transcription factor expression, and gene set enrichment. Lastly, the expression of the four hub genes and their connection to DCM were verified using the rat models. Results: Fibroblasts and monocytes were chosen as hub cells from among the eight identified cell clusters; their marker genes intersected with DEGs to yield six hub genes. In addition, the six hub genes and the essential module genes intersected to yield four essential genes (ASPN, SFRP4, LUM, and FRZB) that were connected to the Wnt signaling pathway and highly expressed in fibroblast. The four hub DEGs had an expression pattern in the DCM rat model experiment results that was in line with the findings of the bioinformatics study. Additionally, there was a strong correlation between decreased cardiac function and the up-regulation of ASPN, SFRP4, LUM, and FRZB. Conclusion: Ultimately, bulk RNA-seq and scRNA-seq data identified fibroblasts and monocytes as the main cell types implicated in DCM. The highly expressed genes ASPN, FRZB, LUM, and SFRP4 in fibroblasts may aid in the mechanistic investigation of DCM.

Erratum: Extensive patient-to-patient single nucleus transcriptome heterogeneity in pheochromocytomas and paragangliomas.

[This corrects the article DOI: 10.3389/fonc.2022.965168.].

Luteinizing Hormone Receptor Mutation (LHRN316S) Causes Abnormal Follicular Development Revealed by Follicle Single-Cell Analysis and CRISPR/Cas9.

Abnormal interaction between granulosa cells and oocytes causes disordered development of ovarian follicles. However, the interactions between oocytes and cumulus granulosa cells (CGs), oocytes and mural granulosa cells (MGs), and CGs and MGs remain to be fully explored. Using single-cell RNA-sequencing (scRNA-seq), we determined the transcriptional profiles of oocytes, CGs and MGs in antral follicles. Analysis of scRNA-seq data revealed that CGs may regulate follicular development through the BMP15-KITL-KIT-PI3K-ARF6 pathway with elevated expression of luteinizing hormone receptor (LHR). Because internalization of the LHR is regulated by Arf6, we constructed LHRN316S mice by CRISPR/Cas9 to further explore mechanisms of follicular development and novel treatment strategies for female infertility. Ovaries of LHRN316S mice exhibited reduced numbers of corpora lutea and ovulation. The LHRN316S mice had a reduced rate of oocyte maturation in vitro and decreased serum progesterone levels. Mating LHRN316S female mice with ICR wild type male mice revealed that the infertility rate of LHRN316S mice was 21.4% (3/14). Litter sizes from LHRN316S mice were smaller than those from control wild type female mice. The oocytes from LHRN316S mice had an increased rate of maturation in vitro after progesterone administration in vitro. Furthermore, progesterone treated LHRN316S mice produced offspring numbers per litter equivalent to WT mice. These findings provide key insights into cellular interactions in ovarian follicles and provide important clues for infertility treatment.

Benchmarking Algorithms for Gene Set Scoring of Single-cell ATAC-seq Data.

Gene set scoring (GSS) has been routinely conducted for gene expression analysis of bulk or single-cell RNA sequencing (RNA-seq) data, which helps to decipher single-cell heterogeneity and cell type-specific variability by incorporating prior knowledge from functional gene sets. Single-cell assay for transposase accessible chromatin using sequencing (scATAC-seq) is a powerful technique for interrogating single-cell chromatin-based gene regulation, and genes or gene sets with dynamic regulatory potentials can be regarded as cell type-specific markers as if in single-cell RNA-seq (scRNA-seq). However, there are few GSS tools specifically designed for scATAC-seq, and the applicability and performance of RNA-seq GSS tools on scATAC-seq data remain to be investigated. Here, we systematically benchmarked ten GSS tools, including four bulk RNA-seq tools, five scRNA-seq tools, and one scATAC-seq method. First, using matched scATAC-seq and scRNA-seq datasets, we found that the performance of GSS tools on scATAC-seq data was comparable to that on scRNA-seq, suggesting their applicability to scATAC-seq. Then, the performance of different GSS tools was extensively evaluated using up to ten scATAC-seq datasets. Moreover, we evaluated the impact of gene activity conversion, dropout imputation, and gene set collections on the results of GSS. Results show that dropout imputation can significantly promote the performance of almost all GSS tools, while the impact of gene activity conversion methods or gene set collections on GSS performance is more dependent on GSS tools or datasets. Finally, we provided practical guidelines for choosing appropriate preprocessing methods and GSS tools in different application scenarios.

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

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

A systematic overview of single-cell transcriptomics databases, their use cases, and limitations.

Rapid advancements in high-throughput single-cell RNA-seq (scRNA-seq) technologies and experimental protocols have led to the generation of vast amounts of transcriptomic data that populates several online databases and repositories. Here, we systematically examined large-scale scRNA-seq databases, categorizing them based on their scope and purpose such as general, tissue-specific databases, disease-specific databases, cancer-focused databases, and cell type-focused databases. Next, we discuss the technical and methodological challenges associated with curating large-scale scRNA-seq databases, along with current computational solutions. We argue that understanding scRNA-seq databases, including their limitations and assumptions, is crucial for effectively utilizing this data to make robust discoveries and identify novel biological insights. Such platforms can help bridge the gap between computational and wet lab scientists through user-friendly web-based interfaces needed for democratizing access to single-cell data. These platforms would facilitate interdisciplinary research, enabling researchers from various disciplines to collaborate effectively. This review underscores the importance of leveraging computational approaches to unravel the complexities of single-cell data and offers a promising direction for future research in the field.

Unveiling heterogeneity and prognostic markers in ductal breast cancer through single-cell RNA-seq.

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. CONCLUSION: 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.

Unsupervised Single-Cell Clustering with Asymmetric Within-Sample Transformation and Per-Cluster Supervised Features Selection.

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.

Zebrafish Thrombocyte Transcriptome Analysis and Functional Genomics.

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.