Recent advances in spatially variable gene detection in spatial transcriptomics

With the emergence of advanced spatial transcriptomic technologies, there has been a surge in research papers dedicated to analyzing spatial transcriptomics data, resulting in significant contributions to our understanding of biology. The initial stage of downstream analysis of spatial transcriptomic data has centered on identifying spatially variable genes (SVGs) or genes expressed with specific spatial patterns across the tissue. SVG detection is an important task since many downstream analyses depend on these selected SVGs. Over the past few years, a plethora of new methods have been proposed for the detection of SVGs, accompanied by numerous innovative concepts and discussions. This article provides a selective review of methods and their practical implementations, offering valuable insights into the current literature in this field.

Aqueous Macrophages Contribute to Conserved CCL2 and CXCL10 Gradients in Uveitis

Purpose

Uveitis is a heterogenous group of inflammatory eye disease for which current cytokine-targeted immune therapies are effective for only a subset of patients. We hypothesized that despite pathophysiologic nuances that differentiate individual disease states, all forms of eye inflammation might share common mechanisms for immune cell recruitment. Identifying these mechanisms is critical for developing novel, broadly acting therapeutic strategies.

Design

Experimental study.

Subjects

Biospecimens from patients with active or inactive uveitis and healthy controls.

Methods

Protein concentration and single cell gene expression were assessed in aqueous fluid biopsies and plasma samples from deidentified patients with uveitis or healthy controls.

Main Outcome Measures

The concentration of 31 inflammatory proteins was measured in all aqueous samples, as well as plasma samples from patients with active uveitis. Chemokine and cytokine ligand and receptor expression were assessed in individual cell types from aqueous biopsies obtained from patients with active uveitis.

Results

We identified 6 chemokines that were both elevated in active uveitis compared with controls and enriched in aqueous compared with plasma during active uveitis (C-C motif chemokine ligand [CCL]2, C-X-C motif chemokine ligand [CXCL]10, CXCL9, CXCL8, CCL3, and CCL14), forming potential gradients for migration of immune cells from the blood to the eye. Of these, CCL2 and CXCL10 were consistently enriched in the aqueous of all patients in our cohort, as well as in a larger cohort of patients from a previously published study. These data suggest that CCL2 and CXCL10 are key mediators in immune cell migration to the eye during uveitis. Next, single cell RNA sequencing suggested that macrophages contribute to aqueous enrichment of CCL2 and CXCL10 during human uveitis. Finally, using chemokine ligand and receptor expression mapping, we identified a broad signaling network for macrophage-derived CCL2 and CXCL10 in human uveitis.

Conclusions

These data suggest that ocular macrophages may play a central role, via CCL2 and CXCL10 production, in recruiting inflammatory cells to the eye in patients with uveitis.

Financial Disclosures

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

The landscape of immune dysregulation in pediatric sepsis at a single-cell resolution

Recognizing immune dysregulation as a hallmark of sepsis pathophysiology, leukocytes have attracted major attention of investigation. While adult and pediatric sepsis are clinically distinct, their immunological delineation remains limited. Single cell technologies facilitated the characterization of immune signatures. We tackled to delineate immunological profiles of pediatric sepsis at a single-cell level by analyzing blood samples from six septic children, at both acute and recovery phases, and four healthy children. 16 single-cell transcriptomic datasets were analyzed and compared to adult sepsis dataset. We showed a unique shift in neutrophil subpopulations and functions between acute and recovery phases, along with the regulatory role of resistin. Neutrophil signatures were comparable between adult and pediatric sepsis. Innate-like CD4 T cells were predominantly and uniquely observed in acute phase of pediatric sepsis. Our study serves as a rich source of information about the phenotypic diversity and trajectory of circulating immune cells during pediatric sepsis.

Single-Cell Transcriptomic Analysis of Epstein-Barr Virus-Associated Hemophagocytic Lymphohistiocytosis

Epstein-Barr virus (EBV) infection can lead to infectious mononucleosis (EBV-IM) and, more rarely, EBV-associated hemophagocytic lymphohistiocytosis (EBV-HLH), which is characterized by a life-threatening hyperinflammatory cytokine storm with immune dysregulation. Interferon-gamma (IFNγ) has been identified as a critical mediator for primary HLH; however, the detailed role of IFNγ and other cytokines in EBV-HLH is not fully understood. In this study, we used single-cell RNA sequencing to characterize the immune landscape of EBV-HLH and compared it with EBV-IM. Three pediatric patients with EBV-HLH with different backgrounds, one with X-linked lymphoproliferative syndrome type 1 (XLP1), two with chronic active EBV disease (CAEBV), and two patients with EBV-IM were enrolled. The TUBA1B + STMN1 + CD8 + T cell cluster, a responsive proliferating cluster with rich mRNA detection, was explicitly observed in EBV-IM, and the upregulation of SH2D1A-the gene responsible for XLP1-was localized in this cluster. This proliferative cluster was scarcely observed in EBV-HLH cases. In EBV-HLH cases with CAEBV, upregulation of LAG3 was observed in EBV-infected cells, which may be associated with an impaired response by CD8 + T cells. Additionally, genes involved in type I interferon (IFN) signaling were commonly upregulated in each cell fraction of EBV-HLH, and activation of type II IFN signaling was observed in CD4 + T cells, natural killer cells, and monocytes but not in CD8 + T cells in EBV-HLH. In conclusion, impaired responsive proliferation of CD8 + T cells and upregulation of type I IFN signaling were commonly observed in EBV-HLH cases, regardless of the patients' background, indicating the key features of EBV-HLH.

919 granules improve postpartum depression through the regulation of abnormal peripheral blood IL-1β

Postpartum depression (PPD) has a significant impact on the physical and mental health of mothers, potentially leading to symptoms such as low mood, fatigue, and decreased appetite. It may also affect the healthy growth of the infant. The onset of PPD is closely related to abnormalities in inflammation and the immune system. PPD patients exhibit abnormalities in the proportion of peripheral blood immune cells, along with an increase in pro-inflammatory cytokines. Excessive pro-inflammatory cytokines in peripheral blood can disrupt the blood-brain barrier (BBB) by activating astrocytes and reducing transendothelial electrical resistance (TEER), allowing peripheral immune cells or cytokines to enter the brain and trigger inflammation, ultimately leading to the onset of depression. In addition, PPD lacks safe and effective treatment medications. In this study, we collected peripheral blood from both healthy postpartum women and those with PPD, conducted single cell RNA sequencing (scRNA-seq), and used an in-house analytical tool scSTAR to reveal that PPD patients exhibit elevated proportions of peripheral blood cDC2 and Proliferation B cells, which are significantly correlated with IL-1β. Additionally, animal experiments were designed to validate that 919 granules can improve PPD by modulating the levels of peripheral blood IL-1β, providing a potential therapeutic mechanism for PPD treatment.

Bone marrow CD8+ Trm cells induced by IL-15 and CD16+ monocytes contribute to HSPC destruction in human severe aplastic anemia

Idiopathic severe aplastic anemia (SAA) is a disease of bone marrow failure caused by T-cell-induced destruction of hematopoietic stem and progenitor cells (HSPCs), however the mechanism remains unclear. We performed single-cell RNA sequencing of PBMCs and BMMCs from SAA patients and healthy donors and identified a CD8+ T cell subset with a tissue residency phenotype (Trm) in bone marrow that exhibit high IFN-γ and FasL expression and have a higher ability to induce apoptosis in HSPCs in vitro through FasL expression. CD8+ Trm cells were induced by IL-15 presented by IL-15Rα on monocytes, especially CD16+ monocytes, which were increased in SAA patients. CD16+ monocytes contributed to IL-15-induced CD38+CXCR6+ pre-Trm differentiation into CD8+ Trm cells, which can be inhibited by the CD38 inhibitor 78c. Our results demonstrate that IL-15-induced CD8+ Trm cells are pathogenic cells that mediate HSPC destruction in SAA patients and are therapeutic targets for future treatments.

Split Pool Ligation-based Single-cell Transcriptome sequencing (SPLiT-seq) data processing pipeline comparison

BACKGROUND

Single-cell sequencing techniques are revolutionizing every field of biology by providing the ability to measure the abundance of biological molecules at a single-cell resolution. Although single-cell sequencing approaches have been developed for several molecular modalities, single-cell transcriptome sequencing is the most prevalent and widely applied technique. SPLiT-seq (split-pool ligation-based transcriptome sequencing) is one of these single-cell transcriptome techniques that applies a unique combinatorial-barcoding approach by splitting and pooling cells into multi-well plates containing barcodes. This unique approach required the development of dedicated computational tools to preprocess the data and extract the count matrices. Here we compare eight bioinformatic pipelines (alevin-fry splitp, LR-splitpipe, SCSit, splitpipe, splitpipeline, SPLiTseq-demultiplex, STARsolo and zUMI) that have been developed to process SPLiT-seq data. We provide an overview of the tools, their computational performance, functionality and impact on downstream processing of the single-cell data, which vary greatly depending on the tool used.

RESULTS

We show that STARsolo, splitpipe and alevin-fry splitp can all handle large amount of data within reasonable time. In contrast, the other five pipelines are slow when handling large datasets. When using smaller dataset, cell barcode results are similar with the exception of SPLiTseq-demultiplex and splitpipeline. LR-splitpipe that is originally designed for processing long-read sequencing data is the slowest of all pipelines. Alevin-fry produced different down-stream results that are difficult to interpret. STARsolo functions nearly identical to splitpipe and produce results that are highly similar to each other. However, STARsolo lacks the function to collapse random hexamer reads for which some additional coding is required.

CONCLUSION

Our comprehensive comparative analysis aids users in selecting the most suitable analysis tool for efficient SPLiT-seq data processing, while also detailing the specific prerequisites for each of these pipelines. From the available pipelines, we recommend splitpipe or STARSolo for SPLiT-seq data analysis.

Single cell analyses of cancer cells identified two regulatorily and functionally distinct categories in differentially expressed genes among tumor subclones

To explore the feature of cancer cells and tumor subclones, we analyzed 101,065 single-cell transcriptomes from 12 colorectal cancer (CRC) patients and 92 single cell genomes from one of these patients. We found cancer cells, endothelial cells and stromal cells in tumor tissue expressed much more genes and had stronger cell-cell interactions than their counterparts in normal tissue. We identified copy number variations (CNVs) in each cancer cell and found correlation between gene copy number and expression level in cancer cells at single cell resolution. Analysis of tumor subclones inferred by CNVs showed accumulation of mutations in each tumor subclone along lineage trajectories. We found differentially expressed genes (DEGs) between tumor subclones had two populations: DEGCNV and DEGreg. DEGCNV, showing high CNV-expression correlation and whose expression differences depend on the differences of CNV level, enriched in housekeeping genes and cell adhesion associated genes. DEGreg, showing low CNV-expression correlation and mainly in low CNV variation regions and regions without CNVs, enriched in cytokine signaling genes. Furthermore, cell-cell communication analyses showed that DEGCNV tends to involve in cell-cell contact while DEGreg tends to involve in secreted signaling, which further support that DEGCNV and DEGreg are two regulatorily and functionally distinct categories.

Deciphering the role of SMARCA4 in cardiac disorders: Insights from single-cell studies on dilated cardiomyopathy and coronary heart disease

BACKGROUND

Dilated cardiomyopathy (DCM) and coronary heart disease (CHD) stand as two of the foremost causes of mortality. However, the comprehensive comprehension of the regulatory mechanisms governing DCM and CHD remains limited, particularly from the vantage point of single-cell transcriptional analysis.

METHOD

We used the GSE121893 dataset from the GEO database, analyzing single-cell expressions with tools like DropletUtils, Seurat, and Monocle. We also utilized the GSVA package for comparing gene roles in DCM and CHD, Finally, we conducted qRT-PCR and Western blot analyses to measure the expression levels of SMARCA4, Col1A1, Col3A1 and α-SMA, and the role of SMARCA4 on fibroblasts were explored by EdU and Transwell assay.

RESULTS

Our analysis identified six cell types in heart tissue, with fibroblasts showing the most interaction with other cells. DEGs in fibroblasts were linked to muscle development and morphogenesis. Pseudotime analysis revealed the dynamics of fibroblast changes in both the normal and disease groups and many transcription factors (TFs) potentially involved in this process. Among these TFs, SMARCA4 which was translated into protein BRG1, showed the most significantly difference. In vivo experiments have demonstrated that SMARCA4 indeed promoted fibroblasts proliferation and migration.

CONCLUSION

This study provides a clearer understanding of cell-type dynamics in heart diseases, emphasizing the role of fibroblasts and the significance of SMARCA4 in their function. Our results offer insights into the cellular mechanisms underlying DCM and CHD, potentially guiding future therapeutic strategies.

Single-cell analysis reveals the transcriptional alterations in the submandibular glands of aged mice

OBJECTIVES

Aging-related salivary gland changes, such as lymphocyte infiltration and acinar cell loss decrease saliva secretion, thereby affecting quality of life. The precise molecular mechanisms underlying these changes remain unclear.

METHODS

We here performed single-cell RNA sequencing to clarify gene expression changes in each cell type comprising the submandibular glands (SMGs) of adult and aged mice.

RESULTS

The proportion of acinar cells decreased in various epithelial clusters annotated with cell type-specific marker genes. Expression levels of the cellular senescence markers, Cdkn2a/p16 and Cdkn1a/p21, were increased in the basal and striated ducts of aged SMGs relative to their levels in those of adult SMGs. In contrast, senescence-associated secretory phenotype-related genes, except transforming growth factor-β, exhibited little change in expression in aged SMGs relative to adult SMGs.

CONCLUSIONS

Gene Ontology analysis revealed increased expression levels of genes encoding major histocompatibility complex (MHC) class I components in the ductal component cells of aged SMGs. MHC class I expression may thus be associated with salivary gland aging.

Mapping the cellular landscape of Atlantic salmon head kidney by single cell and single nucleus transcriptomics

Single-cell transcriptomics is the current gold standard for global gene expression profiling, not only in mammals and model species, but also in non-model fish species. This is a rapidly expanding field, creating a deeper understanding of tissue heterogeneity and the distinct functions of individual cells, making it possible to explore the complexities of immunology and gene expression on a highly resolved level. In this study, we compared two single cell transcriptomic approaches to investigate cellular heterogeneity within the head kidney of healthy farmed Atlantic salmon (Salmo salar). We compared 14,149 cell transcriptomes assayed by single cell RNA-seq (scRNA-seq) with 18,067 nuclei transcriptomes captured by single nucleus RNA-Seq (snRNA-seq). Both approaches detected eight major cell populations in common: granulocytes, heamatopoietic stem cells, erythrocytes, mononuclear phagocytes, thrombocytes, B cells, NK-like cells, and T cells. Four additional cell types, endothelial, epithelial, interrenal, and mesenchymal cells, were detected in the snRNA-seq dataset, but appeared to be lost during preparation of the single cell suspension submitted for scRNA-seq library generation. We identified additional heterogeneity and subpopulations within the B cells, T cells, and endothelial cells, and revealed developmental trajectories of heamatopoietic stem cells into differentiated granulocyte and mononuclear phagocyte populations. Gene expression profiles of B cell subtypes revealed distinct IgM and IgT-skewed resting B cell lineages and provided insights into the regulation of B cell lymphopoiesis. The analysis revealed eleven T cell sub-populations, displaying a level of T cell heterogeneity in salmon head kidney comparable to that observed in mammals, including distinct subsets of cd4/cd8-negative T cells, such as tcrγ positive, progenitor-like, and cytotoxic cells. Although snRNA-seq and scRNA-seq were both useful to resolve cell type-specific expression in the Atlantic salmon head kidney, the snRNA-seq pipeline was overall more robust in identifying several cell types and subpopulations. While scRNA-seq displayed higher levels of ribosomal and mitochondrial genes, snRNA-seq captured more transcription factor genes. However, only scRNA-seq-generated data was useful for cell trajectory inference within the myeloid lineage. In conclusion, this study systematically outlines the relative merits of scRNA-seq and snRNA-seq in Atlantic salmon, enhances understanding of teleost immune cell lineages, and provides a comprehensive list of markers for identifying major cell populations in the head kidney with significant immune relevance.

A review on gene regulatory network reconstruction algorithms based on single cell RNA sequencing

BACKGROUND

Understanding gene regulatory networks (GRNs) is essential for unraveling the molecular mechanisms governing cellular behavior. With the advent of high-throughput transcriptome measurement technology, researchers have aimed to reverse engineer the biological systems, extracting gene regulatory rules from their outputs, which represented by gene expression data. Bulk RNA sequencing, a widely used method for measuring gene expression, has been employed for GRN reconstruction. However, it falls short in capturing dynamic changes in gene expression at the level of individual cells since it averages gene expression across mixed cell populations.

OBJECTIVE

In this review, we provide an overview of 15 GRN reconstruction tools and discuss their respective strengths and limitations, particularly in the context of single cell RNA sequencing (scRNA-seq).

METHODS

Recent advancements in scRNA-seq break new ground of GRN reconstruction. They offer snapshots of the individual cell transcriptomes and capturing dynamic changes. We emphasize how these technological breakthroughs have enhanced GRN reconstruction.

CONCLUSION

GRN reconstructors can be classified based on their requirement for cellular trajectory, which represents a dynamical cellular process including differentiation, aging, or disease progression. Benchmarking studies support the superiority of GRN reconstructors that do not require trajectory analysis in identifying regulator-target relationships. However, methods equipped with trajectory analysis demonstrate better performance in identifying key regulatory factors. In conclusion, researchers should select a suitable GRN reconstructor based on their specific research objectives.

Comparative analysis of peripheral blood immunoinflammatory landscapes in patients with acute cholangitis and its secondary septic shock using single-cell RNA sequencing

BACKGROUND

Acute cholangitis (AC) is a key pathogeny of septic shock, which has a high mortality rate. AC has significant clinical heterogeneity, but no study has analyzed the discrepancies in immunoresponsiveness between AC and its secondary septic shock. The immune inflammatory responses play a critical role in the development of septic shock.

METHODS

We performed single-cell RNA sequencing (scRNA-seq) to analyze the differences of immunocytes in immunoresponse and inflammation between the early stages of AC (A1, A2, and A3) and its secondary septic shock (B1, B2, and B3).

RESULTS

This study has identified seven cell types, including T cells, B cells, plasma cells, neutrophils, monocytes, platelets and erythrocytes. We mainly focused on neutrophils, monocytes, and T cells. Neutrophil subpopulation analysis indicated that neutrophil progenitors (proNeus) were identified in neutrophil subsets. Compared with patients suffering from AC, the gene phenotypes of proNeus (ELANE, AZU1, MPO, and PRTN3) were significantly upregulated in septic shock. The differentiation direction of neutrophil subsets in peripheral blood mononuclear cells (PBMCs) was determined; Moreover, the proNeus in septic shock presented a state of "expansion", with upregulation of neutrophil degranulation and downregulation of monocyte and T cell proliferation. Neutrophils-7 (CCL5, RPL23A, RPL13, RPS19 and RPS18) were mainly involved in the regulation of cellular functions. The neutrophils-7 subpopulation in septic shock were in a state of "exhaustion", and its biological functions showed the characteristics of weakening neutrophil migration and phagocytosis, etc., which maked infection difficult to control and aggravated the development of septic shock. Analysis of monocyte and T cell subpopulations showed that the expression genes and biological functions of subpopulations were closely related to immunoinflammatory regulation. In addition, CCL3 - CCR1, CXCL1 - CXCR2 and other ligand-receptors were highly expressed in neutrophils and monocytes, enhancing interactions between immune cells.

CONCLUSION

ScRNA-seq revealed significant differences in immune cells between AC and its secondary septic shock, which were primarily manifested in the cellular numbers, differentially expressed genes, functions of cellular subsets, differentiation trajectories, cell-cell interactions and so on. We identified many subsets of neutrophil, T cell and monocyte were associated with inflammation and immunosuppression induced by septic shock. These provided a reference for accurately evaluating the pathological severity of patients with AC and discovering the targets for therapy.

A SELECTIVE REVIEW OF RECENT DEVELOPMENTS IN SPATIALLY VARIABLE GENE DETECTION FOR SPATIAL TRANSCRIPTOMICS

With the emergence of advanced spatial transcriptomic technologies, there has been a surge in research papers dedicated to analyzing spatial transcriptomics data, resulting in significant contributions to our understanding of biology. The initial stage of downstream analysis of spatial transcriptomic data has centered on identifying spatially variable genes (SVGs) or genes expressed with specific spatial patterns across the tissue. SVG detection is an important task since many downstream analyses depend on these selected SVGs. Over the past few years, a plethora of new methods have been proposed for the detection of SVGs, accompanied by numerous innovative concepts and discussions. This article provides a selective review of methods and their practical implementations, offering valuable insights into the current literature in this field.

A new and effective two-step clustering approach for single cell RNA sequencing data

BACKGROUND

The rapid devolvement of single cell RNA sequencing (scRNA-seq) technology leads to huge amounts of scRNA-seq data, which greatly advance the research of many biomedical fields involving tissue heterogeneity, pathogenesis of disease and drug resistance etc. One major task in scRNA-seq data analysis is to cluster cells in terms of their expression characteristics. Up to now, a number of methods have been proposed to infer cell clusters, yet there is still much space to improve their performance.

RESULTS

In this paper, we develop a new two-step clustering approach to effectively cluster scRNA-seq data, which is called TSC - the abbreviation of Two-Step Clustering. Particularly, by dividing all cells into two types: core cells (those possibly lying around the centers of clusters) and non-core cells (those locating in the boundary areas of clusters), we first clusters the core cells by hierarchical clustering (the first step) and then assigns the non-core cells to the corresponding nearest clusters (the second step). Extensive experiments on 12 real scRNA-seq datasets show that TSC outperforms the state of the art methods.

CONCLUSION

TSC is an effective clustering method due to its two-steps clustering strategy, and it is a useful tool for scRNA-seq data analysis.

A genomic perspective of the aging human and mouse lung with a focus on immune response and cellular senescence

BACKGROUND

The aging lung is a complex process and influenced by various stressors, especially airborne pathogens and xenobiotics. Additionally, a lifetime exposure to antigens results in structural and functional changes of the lung; yet an understanding of the cell type specific responses remains elusive. To gain insight into age-related changes in lung function and inflammaging, we evaluated 89 mouse and 414 individual human lung genomic data sets with a focus on genes mechanistically linked to extracellular matrix (ECM), cellular senescence, immune response and pulmonary surfactant, and we interrogated single cell RNAseq data to fingerprint cell type specific changes.

RESULTS

We identified 117 and 68 mouse and human genes linked to ECM remodeling which accounted for 46% and 27%, respectively of all ECM coding genes. Furthermore, we identified 73 and 31 mouse and human genes linked to cellular senescence, and the majority code for the senescence associated secretory phenotype. These cytokines, chemokines and growth factors are primarily secreted by macrophages and fibroblasts. Single-cell RNAseq data confirmed age-related induced expression of marker genes of macrophages, neutrophil, eosinophil, dendritic, NK-, CD4+, CD8+-T and B cells in the lung of aged mice. This included the highly significant regulation of 20 genes coding for the CD3-T-cell receptor complex. Conversely, for the human lung we primarily observed macrophage and CD4+ and CD8+ marker genes as changed with age. Additionally, we noted an age-related induced expression of marker genes for mouse basal, ciliated, club and goblet cells, while for the human lung, fibroblasts and myofibroblasts marker genes increased with age. Therefore, we infer a change in cellular activity of these cell types with age. Furthermore, we identified predominantly repressed expression of surfactant coding genes, especially the surfactant transporter Abca3, thus highlighting remodeling of surfactant lipids with implications for the production of inflammatory lipids and immune response.

CONCLUSION

We report the genomic landscape of the aging lung and provide a rationale for its growing stiffness and age-related inflammation. By comparing the mouse and human pulmonary genome, we identified important differences between the two species and highlight the complex interplay of inflammaging, senescence and the link to ECM remodeling in healthy but aged individuals.

Systematic Analysis of Tumor Stem Cell-related Gene Characteristics to Predict the PD-L1 Immunotherapy and Prognosis of Gastric Cancer

AIMS

We aimed to develop a prognostic model with stemness-correlated genes to evaluate prognosis and immunotherapy responsiveness in gastric cancer (GC).

BACKGROUND

Tumor stemness is related to intratumoral heterogeneity, immunosuppression, and anti-tumor resistance. We developed a prognostic model with stemness-correlated genes to evaluate prognosis and immunotherapy responsiveness in GC.

OBJECTIVE

We aimed to develop a prognostic model with stemness-correlated genes to evaluate prognosis and immunotherapy responsiveness in GC.

METHODS

We downloaded single-cell RNA sequencing (scRNA-seq) data of GC patients from the Gene-Expression Omnibus (GEO) database and screened GC stemness- related genes using CytoTRACE. We characterized the association of tumor stemness with immune checkpoint blockade (ICB) and immunity. Thereafter, a 9-stemness signature-based prognostic model was developed using weighted gene co-expression network analysis (WGCNA), univariate Cox regression analysis, and the least absolute shrinkage and selection operator (LASSO) regression analysis. The model predictive value was evaluated with a nomogram.

RESULTS

Early GC patients had significantly higher levels of stemness. The stemness score showed a negative relationship to tumor immune dysfunction and exclusion (TIDE) score and immune infiltration, especially T cells and B cells. A stemness-based signature based on 9 genes (ERCC6L, IQCC, NKAPD1, BLMH, SLC25A15, MRPL4, VPS35, SUMO3, and CINP) was constructed with good performance in prognosis prediction, and its robustness was validated in GSE26942 cohort. Additionally, nomogram and risk score exhibited the most powerful ability for prognosis prediction. High-risk patients exhibited a tendency to develop immune escape and low response to PD-L1 immunotherapy.

CONCLUSION

We developed a stemness-based gene signature for prognosis prediction with accuracy and reliability. This signature also helps clinical decision-making of immunotherapy for GC patients.

Group 1 metabotropic glutamate receptor expression defines a T cell memory population during chronic Toxoplasma infection that enhances IFN-gamma and perforin production in the CNS

Within the brain, a pro-inflammatory response is essential to prevent clinical disease due to Toxoplasma gondii reactivation. Infection in the immunocompromised leads to lethal Toxoplasmic encephalitis while in the immunocompetent, there is persistent low-grade inflammation which is devoid of clinical symptoms. This signifies that there is a well-balanced and regulated inflammatory response to T. gondii in the brain. T cells are the dominant immune cells that prevent clinical disease, and this is mediated through the secretion of effector molecules such as perforins and IFN-γ. The presence of cognate antigen, the expression of survival cytokines, and the alteration of the epigenetic landscape drive the development of memory T cells. However, specific extrinsic signals that promote the formation and maintenance of memory T cells within tissue are poorly understood. During chronic infection, there is an increase in extracellular glutamate that, due to its function as an excitatory neurotransmitter, is normally tightly controlled in the CNS. Here we demonstrate that CD8+ T cells from the T. gondii-infected brain parenchyma are enriched for metabotropic glutamate receptors (mGluR's). Characterization studies determined that mGluR+ expression by CD8+ T cells defines a distinct memory population at the transcriptional and protein level. Finally, using receptor antagonists and agonists we demonstrate mGluR signaling is required for optimal CD8+ T cell production of the effector cytokine IFNγ. This work suggests that glutamate is an important environmental signal of inflammation that promotes T cell function. Understanding glutamate's influence on T cells in the brain can provide insights into the mechanisms that govern protective immunity against CNS-infiltrating pathogens and neuroinflammation.

Exploring the global immune landscape of peripheral blood mononuclear cells in H5N6-infected patient with single-cell transcriptomics

BACKGROUND

Avian influenza viruses (AIV), particularly H5N6, have risen in infection frequency, prompting major concerns. Single-cell RNA sequencing (scRNA-seq) can illustrate the immune cell landscape present in the peripheral circulation of influenza H5N6-infected individuals at the single-cell level. This study attempted to employ scRNA-seq technology to map the potentially hidden single cell landscape of influenza H5N6.

METHODS

High-quality transcriptomes were generated from scRNA-seq data of peripheral blood mononuclear cells (PBMCs), which were taken from a critically-ill child diagnosed with H5N6 avian influenza infection and one healthy control donor. Cluster analysis was then performed on the scRNA-seq data to identify the different cell types. The pathways, pseudotime developmental trajectories and gene regulatory networks involved in different cell subpopulations were also explored.

RESULTS

In total, 3,248 single cell transcriptomes were captured by scRNA-seq from PBMC of the child infected with H5N6 avian influenza and the healthy control donor and further identified seven immune microenvironment cell types. In addition, a subsequent subpopulation analysis of innate lymphoid cells (ILC) and CD4+ T cells revealed that subpopulations of ILC and CD4+ T cells were involved in cytokine and inflammation-related pathways and had significant involvement in the biological processes of oxidative stress and cell death.

CONCLUSION

In conclusion, characterizing the overall immune cell composition of H5N6-infected individuals by assessing the immune cell landscape in the peripheral circulation of H5N6 avian influenza-infected and healthy control donors at single-cell resolution provides key information for understanding H5N6 pathogenesis.

Exploring oncogenes for renal clear cell carcinoma based on G protein-coupled receptor-associated genes

G protein-coupled receptors (GPCRs) are a class of receptors on cell membranes that regulate various biological processes in cells, such as cell proliferation, differentiation, migration, apoptosis, and metabolism, by interacting with G proteins. However, the role of G protein-coupled receptors in predicting the prognosis of renal clear cell carcinoma is still unknown. The transcriptome data and clinical profiles of renal clear cell carcinoma patients, were downloaded from TCGA databases, and the validation group data were downloaded from number GSE167573, including 63 tumor samples and 14 normal samples. Single-cell RNA sequencing data were downloaded from the GEO database, No. GSE152938 and selected samples were used for GSEA enrichment analysis, WGCNA subgroup analysis, single-cell data analysis, and mutation analysis to explore the role of G protein-coupled receptor-related genes in the diagnosis and prognosis of renal clear cell carcinoma and to verify their reliability with cellular experiments. Finally, this study establishes a disease model based on G protein-coupled receptor-related genes, which may help to propose targeted therapeutic regimens in different strata of renal cell carcinoma patients.Author names: Please confirm if the author names are presented accurately and in the correct sequence (given name, middle name/initial, family name). Author: Given name [Lisa Jia] Last name [Tran].It's ok!

β-catenin/CBP activation of mTORC1 signaling promotes partial epithelial-mesenchymal states in head and neck cancer

Head and neck cancers, which include oral squamous cell carcinoma (OSCC) as a major subsite, exhibit cellular plasticity that includes features of an epithelial-mesenchymal transition (EMT), referred to as partial-EMT (p-EMT). To identify molecular mechanisms contributing to OSCC plasticity, we performed a multiphase analysis of single cell RNA sequencing (scRNAseq) data from human OSCC. This included a multiresolution characterization of cancer cell subgroups to identify pathways and cell states that are heterogeneously represented, followed by casual inference analysis to elucidate activating and inhibitory relationships between these pathways and cell states. This approach revealed signaling networks associated with hierarchical cell state transitions, which notably included an association between β-catenin-driven CREB-binding protein (CBP) activity and mTORC1 signaling. This network was associated with subpopulations of cancer cells that were enriched for markers of the p-EMT state and poor patient survival. Functional analyses revealed that β-catenin/CBP induced mTORC1 activity in part through the transcriptional regulation of a raptor-interacting protein, chaperonin containing TCP1 subunit 5 (CCT5). Inhibition of β-catenin-CBP activity through the use of the orally active small molecule, E7386, reduced the expression of CCT5 and mTORC1 activity in vitro, and inhibited p-EMT-associated markers and tumor development in a murine model of OSCC. Our study highlights the use of multiresolution network analyses of scRNAseq data to identify targetable signals for therapeutic benefit, thus defining an underappreciated association between β-catenin/CBP and mTORC1 signaling in head and neck cancer plasticity.

Deciphering the heterogeneity and immunosuppressive function of regulatory T cells in osteosarcoma using single-cell RNA transcriptome

Osteosarcoma (OS) is a highly invasive malignant neoplasm with poor prognosis. The tumor microenvironment (TME) plays an essential role in the occurrence and development of OS. Regulatory T cells (Tregs) are known to facilitate immunosuppression, tumor progression, invasion, and metastasis. However, the effect of Tregs in the TME of OS remains unclear. In this study, single-cell RNA sequencing (scRNA-seq) data was used to identify Tregs and various other cell clusters in the TME of OS. Gene set variation analysis (GSVA) was used to investigate the signaling pathways in Tregs from OS and adjacent tissues. The CellChat and iTALK packages were used to analyze cellular communication. In addition, a prognostic model was established based on the Tregs-specific genes using bulk RNA-seq from the TARGET database, and it was verified using a Gene Expression Omnibus dataset. The pRRophetic package was used to predict drug sensitivity. Immunohistochemistry was used to verify the expression of candidate genes in OS. Based on the above methods, we showed that the OS samples were highly infiltrated with Tregs. GSVA revealed that oxidative phosphorylation, angiogenesis and mammalian target of rapamycin complex 1 (mTORC1) were highly activated in Tregs from OS compared with those from adjacent tissues. Using cellular communication analysis, we found that Tregs interacted with osteoblastic, endothelial, and myeloid cells via C-X-C motif chemokine ligand (CXCL) signaling; particularly, they strongly affected the expression of C-X-C motif chemokine receptor 4 (CXCR4) and interacted with other cell clusters through CXCL12/transforming growth factor β1 (TGFB1) to collectively enable tumor growth and progression. Subsequently, two Tregs-specific genes-CD320 and MAF-were screened through univariate, least absolute shrinkage and selection operator regression (LASSO) and multivariate analysis to construct a prognostic model, which showed excellent prognostic accuracy in two independent cohorts. In addition, drug sensitivity analysis demonstrated that OS patients at high Tregs risk were sensitive to sunitinib, sorafenib, and axitinib. We also used immunohistochemistry to validate that CD320 and MAF were significantly upregulated in OS tissues compared with adjacent tissues. Overall, this study reveals the heterogeneity of Tregs in the OS TME, providing new insights into the invasion and treatment of this cancer.

Age-dependent transcriptomic profiles of leukocytes in pediatric population

Immunity at birth is considered immature. Following birth, our immune function is considered to grow and reach maturation over time. To obtain granular information of leukocyte functions and transcriptomic profiles in pediatric cohort, we examined leukocyte profiles in infants, preschool and school children using single cell RNA sequencing of their peripheral blood mononuclear cells (PBMCs). Monocytes and natural killer (NK) cells showed immaturity in infants. Their innate and adaptive immunity was developed by preschool age. Adaptive immune cells showed different maturation patterns. CD4, CD8 naïve T cells and plasma cells continued to mature untill school age. In CD8 naïve T cells, innate immunity was upregulated in infants, in support of our knowledge that they manifests more innate cell-like phenotype soon after birth. Many signaling pathways have been differentially up- and/or down-regulated in infants, preschool and school children. Their contribution to the development of the immune system needs to be delineated.

Abnormal accumulation of extracellular vesicles in hippocampal dystrophic axons and regulation by the primary cilia in Alzheimer's disease

Dystrophic neurites (DNs) are abnormal axons and dendrites that are swollen or deformed in various neuropathological conditions. In Alzheimer's disease (AD), DNs play a crucial role in impairing neuronal communication and function, and they may also contribute to the accumulation and spread of amyloid beta (Aβ) in the brain of AD patients. However, it is still a challenge to understand the DNs of specific neurons that are vulnerable to Aβ in the pathogenesis of AD. To shed light on the development of radiating DNs, we examined enriched dystrophic hippocampal axons in a mouse model of AD using a three-dimensional rendering of projecting neurons. We employed the anterograde spread of adeno-associated virus (AAV)1 and conducted proteomic analysis of synaptic compartments obtained from hippocampo-septal regions. Our findings revealed that DNs were formed due to synaptic loss at the axon terminals caused by the accumulation of extracellular vesicle (EV). Abnormal EV-mediated transport and exocytosis were identified in association with primary cilia, indicating their involvement in the accumulation of EVs at presynaptic terminals. To further address the regulation of DNs by primary cilia, we conducted knockdown of the Ift88 gene in hippocampal neurons, which impaired EV-mediated secretion of Aβ and promoted accumulation of axonal spheroids. Using single-cell RNA sequencing, we identified the septal projecting hippocampal somatostatin neurons (SOM) as selectively vulnerable to Aβ with primary cilia dysfunction and vesicle accumulation. Our study suggests that DNs in AD are initiated by the ectopic accumulation of EVs at the neuronal axon terminals, which is affected by neuronal primary cilia.

A novel signature of the ligand and receptor genes associated with disulfidoptosis for prediction of prognosis, immunologic therapy responses in hepatocellular carcinoma

Backgroud

We aimed to explore the prognostic features of ligand and receptor genes associated with disulfidoptosis in hepatocellular carcinoma (HCC) and establish a risk signature utilizing these genes to predict the prognosis of HCC patients.

Methods

We used scRNA-seq data from GSE166635 to differentiate malignant cells from normal cells using "copykat".The study thoroughly examined the disparities in disulfidoptosis scores and the associated gene expressions between malignant and normal cells.We identified key ligand and receptor genes that are specific to HCC cells.Subsequently, Correlation analysis was conducted to ascertain the ligand and receptor genes associated with disulfidoptosis.We performed univariate Cox regression analysis to identify prognostic ligand and receptor genes associated with disulfidoptosis.We employed LASSO to construct a risk signature using prognostic ligand and receptor genes associated with disulfidoptosis.Lastly, we developed a nomogram model that integrates the risk signature with clinicopathological characteristics.

Results

Malignant cells displayed a marked increase in disulfidoptosis scores and the expression of associated genes compared to normal cells.We identified 110 receptor and ligand genes significantly associated with disulfidoptosis, and narrowed them down to create a risk signature comprising eight genes.Multivariate analysis confirmed the risk signature as an independent prognostic factor for HCC and validated its predictive value for immunotherapy outcomes.A novel nomogram was developed, incorporating stage information and the risk signature derived from disulfidoptosis-related receptor and ligand genes, demonstrating excellent predictive accuracy and reliability in HCC prognosis prediction.

Conclusion

Risk signatures based on disulfidoptosis-associated ligand and receptor genes can effectively predict HCC prognosis and may inform immunotherapy strategies.

B cell depletion therapy does not resolve chronic active multiple sclerosis lesions

BACKGROUND

Chronic active lesions (CAL) in multiple sclerosis (MS) have been observed even in patients taking high-efficacy disease-modifying therapy, including B-cell depletion. Given that CAL are a major determinant of clinical progression, including progression independent of relapse activity (PIRA), understanding the predicted activity and real-world effects of targeting specific lymphocyte populations is critical for designing next-generation treatments to mitigate chronic inflammation in MS.

METHODS

We analyzed published lymphocyte single-cell transcriptomes from MS lesions and bioinformatically predicted the effects of depleting lymphocyte subpopulations (including CD20 B-cells) from CAL via gene-regulatory-network machine-learning analysis. Motivated by the results, we performed in vivo MRI assessment of PRL changes in 72 adults with MS, 46 treated with anti-CD20 antibodies and 26 untreated, over ∼2 years.

FINDINGS

Although only 4.3% of lymphocytes in CAL were CD20 B-cells, their depletion is predicted to affect microglial genes involved in iron/heme metabolism, hypoxia, and antigen presentation. In vivo, tracking 202 PRL (150 treated) and 175 non-PRL (124 treated), none of the treated paramagnetic rims disappeared at follow-up, nor was there a treatment effect on PRL for lesion volume, magnetic susceptibility, or T1 time. PIRA occurred in 20% of treated patients, more frequently in those with ≥4 PRL (p = 0.027).

INTERPRETATION

Despite predicted effects on microglia-mediated inflammatory networks in CAL and iron metabolism, anti-CD20 therapies do not fully resolve PRL after 2-year MRI follow up. Limited tissue turnover of B-cells, inefficient passage of anti-CD20 antibodies across the blood-brain-barrier, and a paucity of B-cells in CAL could explain our findings.

FUNDING

Intramural Research Program of NINDS, NIH; NINDS grants R01NS082347 and R01NS082347; Dr. Miriam and Sheldon G. Adelson Medical Research Foundation; Cariplo Foundation (grant #1677), FRRB Early Career Award (grant #1750327); Fund for Scientific Research (FNRS).

Time-resolved single-cell RNAseq profiling identifies a novel Fabp5+ subpopulation of inflammatory myeloid cells with delayed cytotoxic profile in chronic spinal cord injury

Traumatic spinal cord injuries (SCI) are a group of highly debilitating pathologies affecting thousands annually, and adversely affecting quality of life. Currently, no fully restorative therapies exist, and SCI still results in significant personal, societal and financial burdens. Inflammation plays a major role in the evolution of SCI, with myeloid cells, including bone marrow derived macrophages (BMDMs) and microglia (MG) being primary drivers of both early secondary pathogenesis and delayed wound healing events. The precise role of myeloid cell subsets is unclear as upon crossing the blood-spinal cord barrier, infiltrating bone marrow derived macrophages (BMDMs) may take on the morphology of resident microglia, and upregulate canonical microglia markers, thus making the two populations difficult to distinguish. Here, we used time-resolved scRNAseq and transgenic fate-mapping to chart the transcriptional profiles of tissue-resident and -infiltrating myeloid cells in a mouse model of thoracic contusion SCI. Our work identifies a novel subpopulation of foam cell-like inflammatory myeloid cells with increased expression of Fatty Acid Binding Protein 5 (Fabp5) and comprise both tissue-resident and -infiltrating cells. Fabp5+ inflammatory myeloid cells display a delayed cytotoxic profile that is predominant at the lesion epicentre and extends into the chronic phase of SCI.

Integration of single-cell regulon atlas and multi-omics data for prognostic stratification and personalized treatment prediction in human lung adenocarcinoma

Transcriptional programs are often dysregulated in cancers. A comprehensive investigation of potential regulons is critical to the understanding of tumorigeneses. We first constructed the regulatory networks from single-cell RNA sequencing data in human lung adenocarcinoma (LUAD). We next introduce LPRI (Lung Cancer Prognostic Regulon Index), a precision oncology framework to identify new biomarkers associated with prognosis by leveraging the single cell regulon atlas and bulk RNA sequencing or microarray datasets. We confirmed that LPRI could be a robust biomarker to guide prognosis stratification across lung adenocarcinoma cohorts. Finally, a multi-omics data analysis to characterize molecular alterations associated with LPRI was performed from The Cancer Genome Atlas (TCGA) dataset. Our study provides a comprehensive chart of regulons in LUAD. Additionally, LPRI will be used to help prognostic prediction and developing personalized treatment for future studies.

scQCEA: a framework for annotation and quality control report of single-cell RNA-sequencing data

BACKGROUND

Systematic description of library quality and sequencing performance of single-cell RNA sequencing (scRNA-seq) data is imperative for subsequent downstream modules, including re-pooling libraries. While several packages have been developed to visualise quality control (QC) metrics for scRNA-seq data, they do not include expression-based QC to discriminate between true variation and background noise.

RESULTS

We present scQCEA (acronym of the single-cell RNA sequencing Quality Control and Enrichment Analysis), an R package to generate reports of process optimisation metrics for comparing sets of samples and visual evaluation of quality scores. scQCEA can import data from 10X or other single-cell platforms and includes functions for generating an interactive report of QC metrics for multi-omics data. In addition, scQCEA provides automated cell type annotation on scRNA-seq data using differential gene expression patterns for expression-based quality control. We provide a repository of reference gene sets, including 2348 marker genes, which are exclusively expressed in 95 human and mouse cell types. Using scRNA-seq data from 56 gene expressions and V(D)J T cell replicates, we show how scQCEA can be applied for the visual evaluation of quality scores for sets of samples. In addition, we use the summary of QC measures from 342 human and mouse shallow-sequenced gene expression profiles to specify optimal sequencing requirements to run a cell-type enrichment analysis function.

CONCLUSIONS

The open-source R tool will allow examining biases and outliers over biological and technical measures, and objective selection of optimal cluster numbers before downstream analysis. scQCEA is available at https://isarnassiri.github.io/scQCEA/ as an R package. Full documentation, including an example, is provided on the package website.

Single cell RNA sequencing reveals mesenchymal heterogeneity and critical functions of Cd271 in tooth development

BACKGROUND

Accumulating evidence suggests that the maxillary process, to which cranial crest cells migrate, is essential to tooth development. Emerging studies indicate that Cd271 plays an essential role in odontogenesis. However, the underlying mechanisms have yet to be elucidated.

AIM

To establish the functionally heterogeneous population in the maxillary process, elucidate the effects of Cd271 deficiency on gene expression differences.

METHODS

p75NTR knockout (Cd271^-/-) mice (from American Jackson laboratory) were used to collect the maxillofacial process tissue of p75NTR knockout mice, and the wild-type maxillofacial process of the same pregnant mouse wild was used as control. After single cell suspension, the cDNA was prepared by loading the single cell suspension into the 10x Genomics Chromium system to be sequenced by NovaSeq6000 sequencing system. Finally, the sequencing data in Fastq format were obtained. The FastQC software is used to evaluate the quality of data and CellRanger analyzed the data. The gene expression matrix is read by R software, and Seurat is used to control and standardize the data, reduce the dimension and cluster. We search for marker genes for subgroup annotation by consulting literature and database; explore the effect of p75NTR knockout on mesenchymal stem cells (MSCs) gene expression and cell proportion by cell subgrouping, differential gene analysis, enrichment analysis and protein-protein interaction network analysis; understand the interaction between MSCs cells and the differentiation trajectory and gene change characteristics of p75NTR knockout MSCs by cell communication analysis and pseudo-time analysis. Last we verified the findings single cell sequencing in vitro.

RESULTS

We identified 21 cell clusters, and we re-clustered these into three subclusters. Importantly, we revealed the cell–cell communication networks between clusters. We clarified that Cd271 was significantly associated with the regulation of mineralization.

CONCLUSION

This study provides comprehensive mechanistic insights into the maxillary- process-derived MSCs and demonstrates that Cd271 is significantly associated with the odontogenesis in mesenchymal populations.

Microtubule stabilization targeting regenerative chondrocyte cluster for cartilage regeneration

Purpose: Chondrocytes (CHs) in cartilage undergo several detrimental events during the development of osteoarthritis (OA). However, the mechanism underlying CHs regeneration involved in pathogenesis is largely unknown. The aim of this study was to explore the underlying mechanism of regeneration of CHs involved in the pathological condition and the potential therapeutic strategies of cartilage repair. Methods and Materials: CHs were isolated from human cartilage in different OA stages and the high-resolution cellular architecture of human osteoarthritis was examined by applying single-cell RNA sequencing. The analysis of gene differential expression and gene set enrichment was utilized to reveal the relationship of cartilage regeneration and microtubule stabilization. Microtubule destabilizer (nocodazole) and microtubule stabilizer (docetaxel) treated-human primary CHs and rats cartilage defect model were used to investing the effects and downstream signaling pathway of microtubule stabilization on cartilage regeneration. Results: CHs subpopulations were identified on the basis of their gene markers and the data indicated an imbalance caused by an increase in the degeneration and disruption of CHs regeneration in OA samples. Interestingly, the CHs subpopulation namely CHI3L1⁺ CHs, was characterized by the cell regenerative capacity, stem cell potency and the activated microtubule (MT) process. Furthermore, the data indicated that MT stabilization was effective in promoting cartilage regeneration in rats with cartilage injury model by inhibiting YAP activity. Conclusion: These findings lead to a new understanding of CHs regeneration in the OA pathophysiology context and suggest that MT stabilization is a promising therapeutic target for OA and cartilage injury.

Single-cell RNA sequencing in cornea research: Insights into limbal stem cells and their niche regulation

The corneal epithelium is composed of stratified squamous epithelial cells on the outer surface of the eye, which acts as a protective barrier and is critical for clear and stable vision. Its continuous renewal or wound healing depends on the proliferation and differentiation of limbal stem cells (LSCs), a cell population that resides at the limbus in a highly regulated niche. Dysfunction of LSCs or their niche can cause limbal stem cell deficiency, a disease that is manifested by failed epithelial wound healing or even blindness. Nevertheless, compared to stem cells in other tissues, little is known about the LSCs and their niche. With the advent of single-cell RNA sequencing, our understanding of LSC characteristics and their microenvironment has grown considerably. In this review, we summarized the current findings from single-cell studies in the field of cornea research and focused on important advancements driven by this technology, including the heterogeneity of the LSC population, novel LSC markers and regulation of the LSC niche, which will provide a reference for clinical issues such as corneal epithelial wound healing, ocular surface reconstruction and interventions for related diseases.

Single-cell RNA seq identifies Plg-RKT-PLG as signals inducing phenotypic transformation of scar-associated macrophage in liver fibrosis

Hepatic macrophages play a central role in liver fibrosis. Scar-associated macrophages (SAMs), a recently identified subgroup of macrophages, play an important role in this process. However, the mechanism by which SAMs transform during liver fibrosis is still unclear. In this study, we aimed to characterize SAMs and elucidate the underlying mechanism of SAM transformation. Bile duct ligation (BDL) and carbon tetrachloride (CCl₄) were used to induce mouse liver fibrosis. Non-parenchymal cells were isolated from normal/fibrotic livers and were analyzed using single cell RNA sequencing (scRNA-seq) or mass cytometry (CyTOF). The glucan-encapsulated siRNA particles (siRNA-GeRPs) was employed to perform macrophage selective gene knockdown. The results of scRNA-seq and CyTOF revealed that SAMs, which derived from bone marrow-derived macrophages (BMMs), accumulated in mouse fibrotic livers. Further analysis showed that SAMs highly expressed genes related to fibrosis, indicating the pro-fibrotic functions of SAMs. Moreover, plasminogen receptor Plg-R_KT was highly expressed by SAMs, suggesting the role of Plg-R_KT and plasminogen (PLG) in SAM transformation. In vitro, PLG-treated BMMs transformed into SAMs and expressed SAM functional genes. Knockdown of Plg-R_KT blocked the effects of PLG. In vivo, selective knockdown of Plg-R_KT in intrahepatic macrophages of BDL- and CCl₄-treated mice reduced the number of SAMs and alleviated BDL- and CCl₄-induced liver fibrosis, suggesting that Plg-R_KT-PLG played an important role in liver fibrosis by mediating SAM transformation. Our findings reveal that SAMs are crucial participants in liver fibrosis. Inhibition of SAM transformation by blocking Plg-R_KT might be a potential therapeutic target for liver fibrosis.

Alemtuzumab treatment exemplifies discordant immune effects of blood and cerebrospinal fluid in multiple sclerosis

BACKGROUND AND OBJECTIVES

Immune responses in the central nervous system (CNS) are highly compartmentalized and cerebrospinal fluid (CSF) in particular often reflects CNS pathology better than peripheral blood. While CSF leukocytes are known to be distinct from blood, the immediate effects of peripheral leukocyte depletion on CSF leukocytes have not been studied in humans.

METHODS

We here analyzed CSF and blood from two relapsing-remitting multiple sclerosis (RRMS) patients early after peripheral leukocyte depletion with the anti-CD52 antibody alemtuzumab compared to untreated RRMS and control patients using single cell RNA-sequencing.

RESULTS

As expected for alemtuzumab, most leukocyte lineages including T cells were synchronously depleted from CSF and blood, while - surprisingly - pDCs were maintained in CSF but depleted from blood by alemtuzumab. Transcriptionally, genes associated with migration were elevated only in the CSF after alemtuzumab. Predicted cellular interactions indicated a central role of pDCs and enhanced migration signaling in the CSF after alemtuzumab.

DISCUSSION

The CSF and blood compartments are thus partially uncoupled, emphasizing that the CNS is only partially accessible even for treatments profoundly affecting the blood.

Multi-modal profiling of peripheral blood cells across the human lifespan reveals distinct immune cell signatures of aging and longevity

BACKGROUND

Age-related changes in immune cell composition and functionality are associated with multimorbidity and mortality. However, many centenarians delay the onset of aging-related disease suggesting the presence of elite immunity that remains highly functional at extreme old age.

METHODS

To identify immune-specific patterns of aging and extreme human longevity, we analyzed novel single cell profiles from the peripheral blood mononuclear cells (PBMCs) of a random sample of 7 centenarians (mean age 106) and publicly available single cell RNA-sequencing (scRNA-seq) datasets that included an additional 7 centenarians as well as 52 people at younger ages (20-89 years).

FINDINGS

The analysis confirmed known shifts in the ratio of lymphocytes to myeloid cells, and noncytotoxic to cytotoxic cell distributions with aging, but also identified significant shifts from CD4+ T cell to B cell populations in centenarians suggesting a history of exposure to natural and environmental immunogens. We validated several of these findings using flow cytometry analysis of the same samples. Our transcriptional analysis identified cell type signatures specific to exceptional longevity that included genes with age-related changes (e.g., increased expression of STK17A, a gene known to be involved in DNA damage response) as well as genes expressed uniquely in centenarians' PBMCs (e.g., S100A4, part of the S100 protein family studied in age-related disease and connected to longevity and metabolic regulation).

INTERPRETATION

Collectively, these data suggest that centenarians harbor unique, highly functional immune systems that have successfully adapted to a history of insults allowing for the achievement of exceptional longevity.

FUNDING

TK, SM, PS, GM, SA, TP are supported by NIH-NIAUH2AG064704 and U19AG023122. MM and PS are supported by NIHNIA Pepper center: P30 AG031679-10. This project is supported by the Flow Cytometry Core Facility at BUSM. FCCF is funded by the NIH Instrumentation grant: S10 OD021587.

The heterogeneous cellular landscape of atherosclerosis: Implications for future research and therapies. A collaborative review from the EAS young fellows

Single cell technologies, lineage tracing mouse models and advanced imaging techniques unequivocally improved the resolution of the cellular landscape of atherosclerosis. Although the discovery of the heterogeneous nature of the cellular plaque architecture has undoubtedly improved our understanding of the specific cellular states in atherosclerosis progression, it also adds more complexity to current and future research and will change how we approach future drug development. In this review, we will discuss how the revolution of new single cell technologies allowed us to map the cellular networks in the plaque, but we will also address current (technological) limitations that confine us to identify the cellular drivers of the disease and to pinpoint a specific cell state, cell subset or cell surface antigen as new candidate drug target for atherosclerosis.

Identification of diagnostic molecules and potential therapeutic agents for atopic dermatitis by single-cell RNA sequencing combined with a systematic computing framework that integrates network pharmacology

Atopic dermatitis (AD) is composed of highly flexible cellular participants. To better understand its pathobiology and molecular regulation mechanisms, it is necessary to combine single-cell RNA sequencing (scRNA-seq) with new computing frameworks or specific technologies, which may contribute to the development of better treatments for AD. The scRNA-seq data of GSE180885 and bulk RNA-seq data of GSE193309 were obtained from Gene Expression Omnibus (GEO) database, and the scRNA-seq data was analyzed by Seurat package to identify the cell types in AD. The genes related to the activity of AD topical drugs were obtained from the ChEMBL database, which provided a variety of bioactivity data such as multiple drugs and targets. AD-related genes were obtained from DisGeNET and CTD databases synthesizing human disease-related genes; the intersection of AD-related genes from these three sources with differentially expressed genes (DEGs) between non-diseased AD and normal human skin (NHS) samples and differential cell type marker genes was taken. The proximity analysis of drug gene network was performed based on the gene with the largest area of receiver operating characteristic (ROC) curve. Ten distinct cell types of AD and NHS were identified, except for phagocytes cells. Three hub genes, F10 and CALCRL and CTSB, were obtained. The area under the curve of ROC based on CTSB expression was the largest, which was 60.15%. By binding drug CTSB-related gene interaction network, we identified 145 potential drugs. Among them, the score of DB07045 and CTSB docking was the lowest, and molecular docking and molecular dynamics (MD) simulation confirmed the close and stable binding of DB07045 and cathepsin B. This work identified diagnostic molecules and potential therapeutic drugs of AD by scRNA-seq combined with a systematic computing framework of network pharmacology, which may provide valuable clues for drug design.

Single cell molecular alterations reveal target cells and pathways of conditioned fear memory

OBJECTIVES

Recent evidence indicates that hippocampus is important for conditioned fear memory (CFM). Though few studies consider the roles of various cell types' contribution to such a process, as well as the accompanying transcriptome changes during this process. The purpose of this study was to explore the transcriptional regulatory genes and the targeted cells that are altered by CFM reconsolidation.

METHODS

A fear conditioning experiment was established on adult male C57 mice, after day 3 tone-cued CFM reconsolidation test, hippocampus cells were dissociated. Using single cell RNA sequencing (scRNA-seq) technique, alterations of transcriptional genes expression were detected and cell cluster analysis were performed and compared with those in sham group.

RESULTS

Seven non-neuronal and eight neuronal cell clusters (including four known neurons and four newly identified neuronal subtypes) has been explored. Among them, CA subtype 1 has characteristic gene markers of Ttr and Ptgds, which is speculated to be the outcome of acute stress and promotes the production of CFM. The results of KEGG pathway enrichment indicate the differences in the expression of certain molecular protein functional subunits in long-term potentiation (LTP) pathway between two types of neurons (DG and CA1) and astrocytes, thus providing a new transcriptional perspective for the role of hippocampus in the CFM reconsolidation. More importantly, the correlation between the reconsolidation of CFM and neurodegenerative diseases-linked genes is substantiated by the results from cell-cell interactions and KEGG pathway enrichment. Further analysis shows that the reconsolidation of CFM inhibits the risk-factor genes App and ApoE in Alzheimer's Disease (AD) and activates the protective gene Lrp1.

CONCLUSIONS

This study reports the transcriptional genes expression changes of hippocampal cells driven by CFM, which confirm the involvement of LTP pathway and suggest the possibility of CFM-like behavior in preventing AD. However, the current research is limited to normal C57 mice, and further studies on AD model mice are needed to prove this preliminary conclusion.

The heterogeneity of fibrosis and angiogenesis in endometriosis revealed by single-cell RNA-sequencing

Fibrosis, angiogenesis and chronic inflammation are the inherent characteristics of endometriosis (EMS). The cellular heterogeneity of ectopic and non-ectopic endometrium by single-cell RNA-sequencing (scRNA-seq)at secretory phase without the disturbance of hormone drugs hasn't been explored so far. In this study, scRNA-seq was adopted to explore the properties of ectopic endometrium (ECE), eutopic endometrium (EUE) and normal endometrium (NOE) at secretory phase. We found that (i) The proportion of myofibroblasts, pericytes, endothelial cells and macrophages in ECE overwhelms that of non-ectopic tissues (EUE and NOE), and Myofibro.C2 was the predominant myofibroblast sub-cluster in ECE. (ii) Myofibroblasts were mainly fibroblast-to-myofibroblast transdifferentiation (FMT) and pericytes were endothelial cell-dependent differentiation in ECE. (iii) Both myofibroblasts and pericytes had a low differentiation potential. (iv) The increased inflammation score, deceased NK cells, T cell exhaustion score and antigen-presenting capacity in ECE confirmed the inflammatory properties and immunodeficiency of ECE. These findings suggested that myofibroblasts, pericytes and macrophages may be the potential targets for anti-fibrotic, anti-angiogenic and anti-inflammatory therapy of EMS.

Myeloid cell heterogeneity in the tumor microenvironment and therapeutic implications for childhood central nervous system (CNS) tumors

Central nervous system (CNS) tumors are the most common type of solid tumors in children and the leading cause of cancer deaths in ages 0-14. Recent advances in the field of tumor biology and immunology have underscored the disparate nature of these distinct CNS tumor types. In this review, we briefly introduce pediatric CNS tumors and discuss various components of the TME, with a particular focus on myeloid cells. Although most studies regarding myeloid cells have been done on adult CNS tumors and animal models, we discuss the role of myeloid cell heterogeneity in pediatric CNS tumors and describe how these cells may contribute to tumorigenesis and treatment response. In addition, we present studies within the last 5 years that highlight human CNS tumors, the utility of various murine CNS tumor models, and the latest multi-dimensional tools that can be leveraged to investigate myeloid cell infiltration in young adults and children diagnosed with select CNS tumors.

Guidance on Processing the 10x Genomics Single Cell Gene Expression Assay

The demand for technologies that allow the study of gene expression at single cell resolution continues to increase. One such assay was launched in 2016 by the US-based company 10x Genomics Inc. Utilizing the power of the single cell on a large scale (Zheng et al. Nat Commun 8:14049, 2017)-capturing thousands of cells at once-has shaped life sciences ever since and allowed researchers to discover new insights within their respective fields of study such as oncology, neurobiology, and immunology (among others). Obtaining high-data quality is the key to being able to make these meaningful discoveries, which in turn is directly linked to the quality of the initial cell (or nuclei) suspension that is used to load the 10x Genomics Chromium Single Cell Gene Expression assay. A successful workflow relies on a cell suspension which is fully dissociated, extremely clean, and of high viability. While the workflow itself has been detailed elsewhere (De Simone et al. Methods Mol Biol 1979:87-110, 2019), in this chapter we will focus on the importance of the quality of the initial cell suspension, as well as common mistakes that can occur while running a Single Cell Gene Expression assay. The descriptions of these tips and tricks refer to the current version of the 10x Genomics User Guide (Chromium Single Cell 3' Reagent Kits User Guide (v3.1 Chemistry Dual Index). https://support10xgenomicscom/single-cell-gene-expression/index/doc/user-guide-chromium-single-cell-3-reagent-kits-user-guide-v31-chemistry-dual-index) which can be downloaded from the Support section on the 10x Genomics website (10x Genomics website. https://www10xgenomicscom). These documents and user guides are continuously improved and updated; hence, it is important to regularly check the company's website for the most recent version.

Ensuring Quality Cell Input for Single Cell Sequencing Experiments by Viability and Singlet Enrichment Using Cell Sorting

As with every sensitive analysis technology, the golden principle of "input quality rules output quality" also applies to single cell sequencing methods. Given the sensitivity of the current methods in single cell sequencing and the minuscule amounts of RNA present within a single cell, any extrinsic source of variability should be reduced by ensuring a homogenous input right at the start. Not every tissue is as readily handled as a single cell suspension like blood and most tissues will have to undergo digestions to free the cells from their spatial organization to undergo single cell transcriptomics workflows. This chapter provides working protocols for two simple, but very precise and powerful methods to ensure only the most viable cells are introduced into single cell assays.

Strategy of Patient-Specific Therapeutics in Cardiovascular Disease Through Single-Cell RNA Sequencing

Recently, single cell RNA sequencing (scRNA-seq) technology has enabled the discovery of novel or rare subtypes of cells and their characteristics. This technique has advanced unprecedented biomedical research by enabling the profiling and analysis of the transcriptomes of single cells at high resolution and throughput. Thus, scRNA-seq has contributed to recent advances in cardiovascular research by the generation of cell atlases of heart and blood vessels and the elucidation of mechanisms involved in cardiovascular development and diseases. This review summarizes the overall workflow of the scRNA-seq technique itself and key findings in the cardiovascular development and diseases based on the previous studies. In particular, we focused on how the single-cell sequencing technology can be utilized in clinical field and precision medicine to treat specific diseases.

Transcriptomic mapping of the metzincin landscape in human trophoblasts

The metzincin family of metalloproteases coordinates tissue developmental processes through regulation of growth factor availability, receptor signaling, and cell-cell/cell-matrix adhesion. While roles for select metzincins in controlling trophoblast functions in human placental development have been described, a comprehensive understanding of metzincin dynamics during trophoblast differentiation is lacking. To address this knowledge gap, single cell transcriptomic datasets derived from first trimester chorionic villi and decidua were used to decipher metzincin expression profiles and kinetics in diverse cell types within the utero-placental interface. Further, specific protease-substrate interactions within progenitor trophoblasts were examined to better define the progenitor niche. Within the uterine-placental compartment, 43 metzincin proteases were expressed across 15 cell-type clusters. Metzincin subgroups expressed in placental trophoblasts, placental mesenchymal cells, uterine stromal, and immune cells included multiple matrix metalloproteases (MMPs), a disintegrin and metalloproteases (ADAMs), a disintegrin and metalloproteases with thrombospondin repeats (ADAMTSs), pappalysins, and astacins. Within the trophoblast compartment, eight distinct trophoblasts states were identified: four cytotrophoblast (CTB), one syncytiotrophoblast precursor (SCTp), two column CTB (cCTB), and one extravillous trophoblast (EVT). Within these states 7 MMP, 8 ADAM, 4 ADAMTS, 2 pappalysin, and 3 astacin proteases were expressed. Cell trajectory modeling shows that expression of most (19/24) metzincins increase during EVT differentiation, though expression of select metalloproteases increase along the villous pathway. Eleven metzincins (ADAM10, -17, MMP14, -15, -19, -23B, ADAMTS1, -6, -19, TLL-1, -2) showed enrichment within CTB progenitors, and analysis of metzincin-substrate interactions identified ∼150 substrates and binding partners, including FBN2 as an ADAMTS6-specific substrate. Together, this work characterizes the metzincin landscape in human first trimester trophoblasts and establishes insight into the roles specific proteases perform within distinct trophoblast niches and across trophoblast differentiation. This resource serves as a guide for future investigations into the roles of metzincin proteases in human placental development.

BISC: accurate inference of transcriptional bursting kinetics from single-cell transcriptomic data

Gene expression in mammalian cells is inherently stochastic and mRNAs are synthesized in discrete bursts. Single-cell transcriptomics provides an unprecedented opportunity to explore the transcriptome-wide kinetics of transcriptional bursting. However, current analysis methods provide limited accuracy in bursting inference due to substantial noise inherent to single-cell transcriptomic data. In this study, we developed BISC, a Bayesian method for inferring bursting parameters from single cell transcriptomic data. Based on a beta-gamma-Poisson model, BISC modeled the mean-variance dependency to achieve accurate estimation of bursting parameters from noisy data. Evaluation based on both simulation and real intron sequential RNA fluorescence in situ hybridization data showed improved accuracy and reliability of BISC over existing methods, especially for genes with low expression values. Further application of BISC found bursting frequency but not bursting size was strongly associated with gene expression regulation. Moreover, our analysis provided new mechanistic insights into the functional role of enhancer and superenhancer by modulating both bursting frequency and size. BISC also formulated a downstream framework to identify differential bursting (in frequency and size separately) genes in samples under different conditions. Applying to multiple datasets (a mouse embryonic cell and fibroblast dataset, a human immune cell dataset and a human pancreatic cell dataset), BISC identified known cell-type signature genes that were missed by differential expression analysis, providing additional insights in understanding the cell-specific stochastic gene transcription. Applying to datasets of human lung and colon cancers, BISC successfully detected tumor signature genes based on alterations in bursting kinetics, which illustrates its value in understanding disease development regarding transcriptional bursting. Collectively, BISC provides a new tool for accurately inferring bursting kinetics and detecting differential bursting genes. This study also produced new insights in the role of transcriptional bursting in regulating gene expression, cell identity and tumor progression.

Stable iPSC-derived NKX2-1+ lung bud tip progenitor organoids give rise to airway and alveolar cell types

Bud tip progenitors (BTPs) in the developing lung give rise to all epithelial cell types found in the airways and alveoli. This work aimed to develop an iPSC organoid model enriched with NKX2-1+ BTP-like cells. Building on previous studies, we optimized a directed differentiation paradigm to generate spheroids with more robust NKX2-1 expression. Spheroids were expanded into organoids that possessed NKX2-1+/CPM+ BTP-like cells, which increased in number over time. Single cell RNA-sequencing analysis revealed a high degree of transcriptional similarity between induced BTPs (iBTPs) and in vivo BTPs. Using FACS, iBTPs were purified and expanded as induced bud tip progenitor organoids (iBTOs), which maintained an enriched population of bud tip progenitors. When iBTOs were directed to differentiate into airway or alveolar cell types using well-established methods, they gave rise to organoids composed of organized airway or alveolar epithelium, respectively. Collectively, iBTOs are transcriptionally and functionally similar to in vivo BTPs, providing an important model for studying human lung development and differentiation.

Immune features of COVID-19 convalescent individuals revealed by a single-cell RNA sequencing

It remains unclear whether immune responses following natural infection can be sustained or potentially prove critical for long-term immune protection against SARS-CoV-2 reinfection. Here, we systematically mapped the phenotypic landscape of SARS-CoV-2-specific immune responses in peripheral blood samples of convalescent patients with COVID-19 by single-cell RNA sequencing. The relative percentage of the CD8 + effector memory subset was increased in both convalescent moderate and severe cases, but NKT-CD160 and marginal zone B clusters were decreased. Innate immune responses were attenuated reflected by decreased expression of genes involved in interferon-gamma, leukocyte migration and neutrophil mediated immune response in convalescent COVID-19 patients. Functions of T cell were strengthened in convalescent COVID-19 patients by clear endorsement of increased expression of genes involved in biological processes of regulation of T cell activation, differentiation and cell-cell adhesion. In addition, T cell mediated immune responses were enhanced with remarkable clonal expansions of TCR and increased transition of CD4 + effector memory and CD8 + effector-GNLY in severe subjects. B cell immune responses displayed complicated and dualfunctions during convalescence of COVID-19, providing a novel mechanism that B cell activation was observed especially in moderate while humoral immune response was weakened. Interestingly, HLA class I genes displayed downregulation while HLA class II genes upregulation in both T and B cell subsets in convalescent individuals. Our results showed that innate immunity was declined but SARS-CoV-2-specific T cell responses were retained even strengthened whereas complicated and dualfunctions of B cells, including declined humoral immunity were presented at several months following infections.

Notch signaling determines cell-fate specification of the two main types of vomeronasal neurons of rodents

The ability of terrestrial vertebrates to find food and mating partners, and to avoid predators, relies on the detection of chemosensory information. Semiochemicals responsible for social and sexual behaviors are detected by chemosensory neurons of the vomeronasal organ (VNO), which transmits information to the accessory olfactory bulb. The vomeronasal sensory epithelium of most mammalian species contains a uniform vomeronasal system; however, rodents and marsupials have developed a more complex binary vomeronasal system, containing vomeronasal sensory neurons (VSNs) expressing receptors of either the V1R or V2R family. In rodents, V1R/apical and V2R/basal VSNs originate from a common pool of progenitors. Using single cell RNA-sequencing, we identified differential expression of Notch1 receptor and Dll4 ligand between the neuronal precursors at the VSN differentiation dichotomy. Our experiments show that Notch signaling is required for effective differentiation of V2R/basal VSNs. In fact, Notch1 loss of function in neuronal progenitors diverts them to the V1R/apical fate, whereas Notch1 gain of function redirects precursors to V2R/basal. Our results indicate that Notch signaling plays a pivotal role in triggering the binary differentiation dichotomy in the VNO of rodents.

Visualisation of tumour cells in bone in vivo at single-cell resolution

The skeleton is a common site for the establishment of distant metastases. Once cancers occupy bone, the prognosis is poor as disease recurrence and visceral spread is imminent. Understanding the pathways and cellular interactions, which regulate tumour cell seeding, dormancy and growth in bone, is pertinent to improving outcomes for patients with advanced cancers. Advances in imaging techniques have facilitated the development of the concept that the behavior of bone marrow resident cells dictates the fate of tumour cells upon arrival in bone. This review summarises recent findings achieved through intravital imaging. It highlights the importance of developing both longitudinal static and acute dynamic data to develop our understanding of tumour cell engraftment, dormancy, activation and the subsequent establishment of metastases. We also describe how imaging techniques have developed our knowledge of the elements that make up the complex bone microenvironment which tumour cells interact with to survive and grow. We also discuss how through combining these imaging insights with single cell RNA sequencing data, we are entering a new era of research which has the power to define the cell-cell interactions which control tumour cell growth in bone.

Culture shock: microglial heterogeneity, activation, and disrupted single-cell microglial networks in vitro

BACKGROUND

Microglia, the resident immune cells of the brain, play a critical role in numerous diseases, but are a minority cell type and difficult to genetically manipulate in vivo with viral vectors and other approaches. Primary cultures allow a more controlled setting to investigate these cells, but morphological and transcriptional changes upon removal from their normal brain environment raise many caveats from in vitro studies.

METHODS

To investigate whether cultured microglia recapitulate in vivo microglial signatures, we used single-cell RNA sequencing (scRNAseq) to compare microglia freshly isolated from the brain to primary microglial cultures. We performed cell population discovery, differential expression analysis, and gene co-expression module analysis to compare signatures between in vitro and in vivo microglia. We constructed causal predictive network models of transcriptional regulators from the scRNAseq data and identified a set of potential key drivers of the cultured phenotype. To validate this network analysis, we knocked down two of these key drivers, C1qc and Prdx1, in primary cultured microglia and quantified changes in microglial activation markers.

RESULTS

We found that, although often assumed to be a relatively homogenous population of cells in culture, in vitro microglia are a highly heterogeneous population consisting of distinct subpopulations of cells with transcriptional profiles reminiscent of macrophages and monocytes, and are marked by transcriptional programs active in neurodegeneration and other disease states. We found that microglia in vitro presented transcriptional activation of a set of "culture shock genes" not found in freshly isolated microglia, characterized by strong upregulation of disease-associated genes including Apoe, Lyz2, and Spp1, and downregulation of homeostatic microglial markers, including Cx3cr1, P2ry12, and Tmem119. Finally, we found that cultured microglia prominently alter their transcriptional machinery modulated by key drivers from the homeostatic to activated phenotype. Knockdown of one of these drivers, C1qc, resulted in downregulation of microglial activation genes Lpl, Lyz2, and Ccl4.

CONCLUSIONS

Overall, our data suggest that when removed from their in vivo home environment, microglia suffer a severe case of "culture shock", drastically modulating their transcriptional regulatory network state from homeostatic to activated through upregulation of modules of culture-specific genes. Consequently, cultured microglia behave as a disparate cell type that does not recapitulate the homeostatic signatures of microglia in vivo. Finally, our predictive network model discovered potential key drivers that may convert activated microglia back to their homeostatic state, allowing for more accurate representation of in vivo states in culture. Knockdown of key driver C1qc partially attenuated microglial activation in vitro, despite C1qc being only weakly upregulated in culture. This suggests that even genes that are not strongly differentially expressed across treatments or preparations may drive downstream transcriptional changes in culture.