Single-cell RNA Sequencing Studies Research Articles

Integrative analysis of COL6A3 in lupus nephritis: insights from single-cell transcriptomics and proteomics

IntroductionLupus nephritis (LN), a severe complication of systemic lupus erythematosus (SLE), presents significant challenges in patient management and treatment outcomes. The identification of novel LN-related biomarkers and therapeutic targets is critical to enhancing treatment outcomes and prognosis for patients.MethodsIn this study, we analyzed single-cell expression data from LN (n=21) and healthy controls (n=3). A total of 143 differentially expressed genes were identified between the LN and control groups. Then, proteomics analysis of LN patients (n=9) and control (SLE patients without LN, n=11) revealed 55 differentially expressed genes among patients with LN and control group. We further utilizes protein-protein interaction network and functional enrichment analyses to elucidate the pivotal role of COL6A3 in key signaling pathways. Its diagnostic value is evaluate through its correlation with disease progression and renal function metrics, as well as Receiver Operating Characteristic Curve (ROC) analysis. Additionally, immunohistochemistry and qPCR experiments were performed to validate the expression of COL6A3 in LN.ResultsBy comparison of single-cell and proteomics data, we discovered that COL6A3 is significantly upregulated, highlighting it as a critical biomarker of LN. Our findings emphasize the substantial involvement of COL6A3 in the pathogenesis of LN, particularly noting its expression in mesangial cells. Through comprehensive protein-protein interaction network and functional enrichment analyses, we uncovered the pivotal role of COL6A3 in key signaling pathways including integrin-mediated signaling pathways, collagen-activated signaling pathways, and ECM-receptor interaction, suggesting potential therapeutic targets. The diagnostic utility is confirmed by its correlation with disease progression and renal function metrics of the glomerular filtration rate. ROC analysis further validates the diagnostic value of COL6A3, with the area under the ROC values of 0.879 in the in-house cohort, and 0.802 and 0.915 in tubular and glomerular external cohort samples, respectively. Furthermore, immunohistochemistry and qPCR experiments were consistent with those obtained from the single-cell RNA sequencing and proteomics studies. DiscussionThese results proved that COL6A3 is a promising biomarker and therapeutic target, advancing personalized medicine strategies for LN.

Transcriptomic imputation of genetic risk variants uncovers novel whole-blood biomarkers of Parkinson’s disease

Blood-based gene expression signatures could potentially be used as biomarkers for PD. However, it is unclear whether genetically-regulated transcriptomic signatures can provide novel gene candidates for use as PD biomarkers. We leveraged on the Genotype-Tissue Expression (GTEx) database to impute whole-blood transcriptomic expression using summary statistics of three large-scale PD GWAS. A random forest classifier was used with the consensus whole-blood imputed gene signature (IGS) to discriminate between cases and controls. Outcome measures included Area under the Curve (AUC) of Receiver Operating Characteristic (ROC) Curve. We demonstrated that the IGS (n = 37 genes) is conserved across PD GWAS studies and brain tissues. IGS discriminated between cases and controls in an independent whole-blood RNA-sequencing study (1176 PD, 254 prodromal, and 860 healthy controls) with mean AUC and accuracy of 64.8% and 69.4% for PD cohort, and 78.8% and 74% for prodromal cohort. PATL2 was the top-performing imputed gene in both PD and prodromal PD cohorts, whose classifier performance varied with biological sex (higher performance for males and females in the PD and prodromal PD, respectively). Single-cell RNA-sequencing studies (scRNA-seq) of healthy humans and PD patients found PATL2 to be enriched in terminal effector CD8+ and cytotoxic CD4+ cells, whose proportions are both increased in PD patients. We demonstrated the utility of GWAS transcriptomic imputation in identifying novel whole-blood transcriptomic signatures which could be leveraged upon for PD biomarker derivation. We identified PATL2 as a potential biomarker in both clinical and prodromic PD.

Highlight of 2023: Advances in mast cells.

In this article for the Highlights of 2023 Series, we consider the growing understanding of mast cell heterogeneity and interactions that has developed from single cell RNA sequencing studies. We also discuss novel concepts concerning mast cell interactions with the central nervous system and evidence for their role in host defense against SARS-CoV-2 infection.

Cellular plasticity in the breast cancer ecosystem.

The complex interplay between genetically diverse tumor cells and their microenvironment significantly influences cancer progression and therapeutic responses. This review highlights recent findings on cellular plasticity and heterogeneity within the breast cancer ecosystem, focusing on the roles of cancer-associated fibroblasts (CAFs) and tumor-associated macrophages (TAMs). We discuss evidence suggesting that breast cancer cells exhibit phenotypic plasticity driven by both intrinsic genetic factors and external microenvironmental cues, impacting treatment responses and disease recurrence. Moreover, single-cell RNA sequencing studies reveal diverse subtypes of CAFs and TAMs, each with distinct functional gene expression programs and spatial organization within the tumor microenvironment. Understanding the hierarchical relationships and niche cues governing cellular phenotypes offers new opportunities for targeted therapeutic interventions. By elucidating the organizational principles of the tumor ecosystem, future therapies may target phenotypic states or entire cellular niches, advancing precision medicine approaches in breast cancer treatment.

Autoimmunity-associated allele of tyrosine phosphatase gene PTPN22 enhances anti-viral immunity.

The 1858C>T allele of the tyrosine phosphatase PTPN22 is present in 5-10% of the North American population and is strongly associated with numerous autoimmune diseases. Although research has been done to define how this allele potentiates autoimmunity, the influence PTPN22 and its pro-autoimmune allele has in anti-viral immunity remains poorly defined. Here, we use single cell RNA-sequencing and functional studies to interrogate the impact of this pro-autoimmune allele on anti-viral immunity during Lymphocytic Choriomeningitis Virus clone 13 (LCMV-cl13) infection. Mice homozygous for this allele (PEP-619WW) clear the LCMV-cl13 virus whereas wildtype (PEP-WT) mice cannot. This is associated with enhanced anti-viral CD4 T cell responses and a more immunostimulatory CD8α- cDC phenotype. Adoptive transfer studies demonstrated that PEP-619WW enhanced anti-viral CD4 T cell function through virus-specific CD4 T cell intrinsic and extrinsic mechanisms. Taken together, our data show that the pro-autoimmune allele of Ptpn22 drives a beneficial anti-viral immune response thereby preventing what is normally a chronic virus infection.

Disease-associated astrocyte epigenetic memory promotes CNS pathology.

Disease-associated astrocyte subsets contribute to the pathology of neurologic diseases, including multiple sclerosis and experimental autoimmune encephalomyelitis1-8 (EAE), an experimental model for multiple sclerosis. However, little is known about the stability of these astrocyte subsets and their ability to integrate past stimulation events. Here we report the identification of an epigenetically controlled memory astrocyte subset that exhibits exacerbated pro-inflammatory responses upon rechallenge. Specifically, using a combination of single-cell RNA sequencing, assay for transposase-accessible chromatin with sequencing, chromatin immunoprecipitation with sequencing, focused interrogation of cells by nucleic acid detection and sequencing, and cell-specific in vivo CRISPR-Cas9-based genetic perturbation studies we established that astrocyte memory is controlled by the metabolic enzyme ATP-citrate lyase (ACLY), which produces acetyl coenzyme A (acetyl-CoA) that is used by histone acetyltransferase p300 to control chromatin accessibility. The number of ACLY+p300+ memory astrocytes is increased in acute and chronic EAE models, and their genetic inactivation ameliorated EAE. We also detected the pro-inflammatory memory phenotype in human astrocytes in vitro; single-cell RNA sequencing and immunohistochemistry studies detected increased numbers of ACLY+p300+ astrocytes in chronic multiple sclerosis lesions. In summary, these studies define an epigenetically controlled memory astrocyte subset that promotes CNS pathology in EAE and, potentially, multiple sclerosis. These findings may guide novel therapeutic approaches for multiple sclerosis and other neurologic diseases.

Revealing the molecular landscape of human placenta: a systematic review and meta-analysis of single-cell RNA sequencing studies.

With increasing significance of developmental programming effects associated with placental dysfunction, more investigations are devoted to improving the characterization and understanding of placental signatures in health and disease. The placenta is a transitory but dynamic organ adapting to the shifting demands of fetal development and available resources of the maternal supply throughout pregnancy. Trophoblasts (cytotrophoblasts, syncytiotrophoblasts, and extravillous trophoblasts) are placental-specific cell types responsible for the main placental exchanges and adaptations. Transcriptomic studies with single-cell resolution have led to advances in understanding the placenta's role in health and disease. These studies, however, often show discrepancies in characterization of the different placental cell types. We aim to review the knowledge regarding placental structure and function gained from the use of single-cell RNA sequencing (scRNAseq), followed by comparing cell-type-specific genes, highlighting their similarities and differences. Moreover, we intend to identify consensus marker genes for the various trophoblast cell types across studies. Finally, we will discuss the contributions and potential applications of scRNAseq in studying pregnancy-related diseases. We conducted a comprehensive systematic literature review to identify different cell types and their functions at the human maternal-fetal interface, focusing on all original scRNAseq studies on placentas published before March 2023 and published reviews (total of 28 studies identified) using PubMed search. Our approach involved curating cell types and subtypes that had previously been defined using scRNAseq and comparing the genes used as markers or identified as potential new markers. Next, we reanalyzed expression matrices from the six available scRNAseq raw datasets with cell annotations (four from first trimester and two at term), using Wilcoxon rank-sum tests to compare gene expression among studies and annotate trophoblast cell markers in both first trimester and term placentas. Furthermore, we integrated scRNAseq raw data available from 18 healthy first trimester and nine term placentas, and performed clustering and differential gene expression analysis. We further compared markers obtained with the analysis of annotated and raw datasets with the literature to obtain a common signature gene list for major placental cell types. Variations in the sampling site, gestational age, fetal sex, and subsequent sequencing and analysis methods were observed between the studies. Although their proportions varied, the three trophoblast types were consistently identified across all scRNAseq studies, unlike other non-trophoblast cell types. Notably, no marker genes were shared by all studies for any of the investigated cell types. Moreover, most of the newly defined markers in one study were not observed in other studies. These discrepancies were confirmed by our analysis on trophoblast cell types, where hundreds of potential marker genes were identified in each study but with little overlap across studies. From 35461 and 23378 cells of high quality in the first trimester and term placentas, respectively, we obtained major placental cell types, including perivascular cells that previously had not been identified in the first trimester. Importantly, our meta-analysis provides marker genes for major placental cell types based on our extensive curation. This review and meta-analysis emphasizes the need for establishing a consensus for annotating placental cell types from scRNAseq data. The marker genes identified here can be deployed for defining human placental cell types, thereby facilitating and improving the reproducibility of trophoblast cell annotation.

Abstract P35: High-grade B-cell lymphomas with 11q aberrations show a dark zone expression profile similar to Burkitt lymphoma but with additional extension towards intermediate germinal center B-cells

Abstract P35: High-grade B-cell lymphomas with 11q aberrations show a dark zone expression profile similar to Burkitt lymphoma but with additional extension towards intermediate germinal center B-cells

Cardiac and perivascular myofibroblasts, matrifibrocytes, and immune fibrocytes in hypertension; commonalities and differences with other cardiovascular diseases.

Hypertension is a major cause of cardiovascular diseases such as myocardial infarction and stroke. Cardiovascular fibrosis occurs with hypertension and contributes to vascular resistance, aortic stiffness, and cardiac hypertrophy. However, the molecular mechanisms leading to fibroblast activation in hypertension remain largely unknown. There are two types of fibrosis: replacement fibrosis and reactive fibrosis. Replacement fibrosis occurs in response to the loss of viable tissue to form a scar. Reactive fibrosis occurs in response to an increase in mechanical and neurohormonal stress. Although both types of fibrosis are considered adaptive processes, they become maladaptive when the tissue loss is too large, or the stress persists. Myofibroblasts represent a subpopulation of activated fibroblasts that have gained contractile function to promote wound healing. Therefore, myofibroblasts are a critical cell type that promotes replacement fibrosis. Although myofibroblasts were recognized as the fibroblasts participating in reactive fibrosis, recent experimental evidence indicated there are distinct fibroblast populations in cardiovascular reactive fibrosis. Accordingly, we will discuss the updated definition of fibroblast subpopulations, the regulatory mechanisms, and their potential roles in cardiovascular pathophysiology utilizing new knowledge from various lineage tracing and single-cell RNA sequencing studies. Among the fibroblast subpopulations, we will highlight the novel roles of matrifibrocytes and immune fibrocytes in cardiovascular fibrosis including experimental models of hypertension, pressure overload, myocardial infarction, atherosclerosis, aortic aneurysm, and nephrosclerosis. Exploration into the molecular mechanisms involved in the differentiation and activation of those fibroblast subpopulations may lead to novel treatments for end-organ damage associated with hypertension and other cardiovascular diseases.

Novel Approach for Pancreas Transcriptomics Reveals the Cellular Landscape in Homeostasis and Acute Pancreatitis

Acinar cells produce digestive enzymes that impede transcriptomic characterization of the exocrine pancreas. Thus, single-cell RNA-sequencing studies of the pancreas underrepresent acinar cells relative to histological expectations, and a robust approach to capture pancreatic cell responses in disease states is needed. We sought to innovate a method that overcomes these challenges to accelerate study of the pancreas in health and disease. We leverage FixNCut, a single-cell RNA-sequencing approach in which tissue is reversibly fixed with dithiobis(succinimidyl propionate) before dissociation and single-cell preparation. We apply FixNCut to an established mouse model of acute pancreatitis, validate findings using GeoMx whole transcriptome atlas profiling, and integrate our data with prior studies to compare our method in both mouse and human pancreas datasets. FixNCut achieves unprecedented definition of challenging pancreatic cells, including acinar and immune populations in homeostasis and acute pancreatitis, and identifies changes in all major cell types during injury and recovery. We define the acinar transcriptome during homeostasis and acinar-to-ductal metaplasia and establish a unique gene set to measure deviation from normal acinar identity. We characterize pancreatic immune cells, and analysis of T-cell subsets reveals a polarization of the homeostatic pancreas toward type-2 immunity. We report immune responses during acute pancreatitis and recovery, including early neutrophil infiltration, expansion of dendritic cell subsets, and a substantial shift in the transcriptome of macrophages due to both resident macrophage activation and monocyte infiltration. FixNCut preserves pancreatic transcriptomes to uncover novel cell states during homeostasis and following pancreatitis, establishing a broadly applicable approach and reference atlas for study of pancreas biology and disease.

Macrophage profiling in atherosclerosis: understanding the unstable plaque.

The development and rupture of atherosclerotic plaques is a major contributor to myocardial infarctions and ischemic strokes. The dynamic evolution of the plaque is largely attributed to monocyte/macrophage functions, which respond to various stimuli in the plaque microenvironment. To this end, macrophages play a central role in atherosclerotic lesions through the uptake of oxidized low-density lipoprotein that gets trapped in the artery wall, and the induction of an inflammatory response that can differentially affect the stability of the plaque in men and women. In this environment, macrophages can polarize towards pro-inflammatory M1 or anti-inflammatory M2 phenotypes, which represent the extremes of the polarization spectrum that include Mhem, M(Hb), Mox, and M4 populations. However, this traditional macrophage model paradigm has been redefined to include numerous immune and nonimmune cell clusters based on in-depth unbiased single-cell approaches. The goal of this review is to highlight (1) the phenotypic and functional properties of monocyte subsets in the circulation, and macrophage populations in atherosclerotic plaques, as well as their contribution towards stable or unstable phenotypes in men and women, and (2) single-cell RNA sequencing studies that have advanced our knowledge of immune, particularly macrophage signatures present in the atherosclerotic niche. We discuss the importance of performing high-dimensional approaches to facilitate the development of novel sex-specific immunotherapies that aim to reduce the risk of cardiovascular events.

An advanced comprehensive muti-cell-type-specific model for predicting anti-PD-1 therapeutic effect in melanoma.

Rationale: Immune checkpoint inhibitors targeting the programmed cell death (PD)-1/PD-L1 pathway have promise in patients with advanced melanoma. However, drug resistance usually results in limited patient benefits. Recent single-cell RNA sequencing studies have elucidated that MM patients display distinctive transcriptional features of tumor cells, immune cells and interstitial cells, including loss of antigen presentation function of tumor cells, exhaustion of CD8+T and extracellular matrix secreted by fibroblasts to prevents immune infiltration, which leads to a poor response to immune checkpoint inhibitors (ICIs). However, cell subgroups beneficial to anti-tumor immunity and the model developed by them remain to be further identified. Methods: In this clinical study of neoadjuvant therapy with anti-PD-1 in advanced melanoma, tumor tissues were collected before and after treatment for single-nucleus sequencing, and the results were verified using multicolor immunofluorescence staining and public datasets. Results: This study describes four cell subgroups which are closely associated with the effectiveness of anti-PD-1 treatment. It also describes a cell-cell communication network, in which the interaction of the four cell subgroups contributes to anti-tumor immunity. Furthermore, we discuss a newly developed predictive model based on these four subgroups that holds significant potential for assessing the efficacy of anti-PD-1 treatment. Conclusions: These findings elucidate the primary mechanism of anti-PD-1 resistance and offer guidance for clinical drug administration for melanoma.

Comprehensive Integration of Multiple Single-Cell Transcriptomic Data Sets Defines Distinct Cell Populations and Their Phenotypic Changes in Murine Atherosclerosis.

The application of single-cell transcriptomic (single-cell RNA sequencing) analysis to the study of atherosclerosis has provided unique insights into the molecular and genetic mechanisms that mediate disease risk and pathophysiology. However, nonstandardized methodologies and relatively high costs associated with the technique have limited the size and replication of existing data sets and created disparate or contradictory findings that have fostered misunderstanding and controversy. To address these uncertainties, we have performed a conservative integration of multiple published single-cell RNA sequencing data sets into a single meta-analysis, performed extended analysis of native resident vascular cells, and used in situ hybridization to map the disease anatomic location of the identified cluster cells. To investigate the transdifferentiation of smooth muscle cells to macrophage phenotype, we have developed a classifying algorithm based on the quantification of reporter transgene expression. The reporter gene expression tool indicates that within the experimental limits of the examined studies, transdifferentiation of smooth muscle cell to the macrophage lineage is extremely rare. Validated transition smooth muscle cell phenotypes were defined by clustering, and the location of these cells was mapped to lesion anatomy with in situ hybridization. We have also characterized 5 endothelial cell phenotypes and linked these cellular species to different vascular structures and functions. Finally, we have identified a transcriptomically unique cellular phenotype that constitutes the aortic valve. Taken together, these analyses resolve a number of outstanding issues related to differing results reported with vascular disease single-cell RNA sequencing studies, and significantly extend our understanding of the role of resident vascular cells in anatomy and disease.

Contribution of collagen XIII to lung function and development of pulmonary fibrosis

BackgroundCollagen XIII is a transmembrane collagen associated with neuromuscular junction development, and in humans its deficiency results in congenital myasthenic syndrome type 19 (CMS19), which leads to breathing difficulties. CMS19...

Natural Killer Lymphocytes Mediate Renal Fibrosis Due to Acute Cardiorenal Syndrome.

The acute kidney injury (AKI) to chronic kidney disease (CKD) transition presents an opportunity for intervention to prevent CKD. Our lab developed a novel murine model of AKI-CKD transition, cardiac arrest and cardiopulmonary resuscitation (CA/CPR), in which all animals develop CKD at 7 weeks. The purpose of this study was to identify potential immune drivers of fibrosis after CA/CPR. Cardiac arrest was induced by potassium chloride and mice were resuscitated with chest compressions and epinephrine. The kidney immune landscape after CA/CPR was profiled using 11-color flow cytometry analysis and immunofluorescence. Immune cell-derived mediators of fibrosis were identified by analyzing data from three previously published single cell or single nuclear RNA sequencing (sc/snRNASeq) studies. NRK49F fibroblasts were treated with granzyme A (GzA) in vitro, then cell proliferation was quantified using 5-ethynyl-2'-deoxyuridine. GzA was pharmacologically inhibited both in vitro and in vivo. Immune cells infiltrated the kidney after CA/CPR, consisting primarily of innate immune cells, including monocytes/macrophages, neutrophils, and natural killer (NK) cells. NK cell infiltration immediately preceded mesenchymal cell expansion, which occurred starting 7 days after CA/CPR. Immune cells colocalized with mesenchymal cells, accumulating in areas of fibrosis. Analysis of previously published sc/snRNASeq data revealed GzA as a potential mediator of immune to mesenchymal communication. GzA administration to fibroblasts in vitro induced cell growth and proliferation. Pharmacologic blockade of GzA signaling in vivo attenuated fibrosis and improved renal function after CA/CPR. Renal inflammation occurs during cardiorenal syndrome, which correlates with mesenchymal cell expansion. GzA, produced by NK cells, presents a novel therapeutic target to prevent the transition to CKD after AKI.

Mechanisms of Aging Vary with Developmental Stages in the Human Hematopoietic Trajectory

Background : Hematopoietic stem and progenitor cells (HSPCs) undergo functional changes with age. An unresolved issue is whether aging is induced by a shift in intrinsic state in which all HSPCs undergo coordinated changes in functional potential, or by clonal evolution of a senescent population.Identification of senescent cells in human hematopoiesis poses another challenge. In this study, we have identified lineage-specific, senescence-associated pathways at different developmental stages of HSPCs and have developed an “aging signature HSC” for non-primed human HSCs. Methods: HSPCs (CD34+ cells) were isolated from the bone marrow of healthy human subjects (old = >60 years; young = 21 to <35 years). Single-cell RNA-sequencing (scRNAseq) studies and downstream analysis were performed using different R packages for normalization, dimension reduction, clustering, classification, pseudotime ordering, and developmental trajectory analyses. Gene Set Enrichment Analysis, GSEA (version 4.3.2., Broad Institute, Inc) was applied to identify pathways and processes that were coordinately up- or down-regulated with aging. Results: HSPCs (CD34+ cells) were derived from 13 healthy human subjects (old subjects n = 6; young subjects n = 7). There were 26,076 genes identified. The HSPCs were ordered and classified along their developmental trajectory during differentiation processes from (1) hematopoietic stem cells (HSC) to (2) myeloid-lymphoid progenitors (MLP), (3) megakaryocyte-erythrocyte progenitors (MEP), (4) granulocytic-monocytic progenitors (GMP), (5) granulocyte- (G), and (6) monocyte- (Mono) progenitors, (7) dendritic cell (DC), (8) erythroid progenitors (E), as well as various subsets of lymphocytes (Ly1 to Ly4). The results are depicted in Figure 1A. Remarkable are the increases in HSCs (p < 0.05) and in Ly4 (p < 0.001) in old subjects. The Kolmogorov-Smirnov test showed a significant delay in the developmental “pseudotimes” of HSCs (p < 0.001) and of Ly4 (p < 0.05). We then compared the differences in GSEA between old and young subjects in each of the following clusters: (A) whole HSPCs; (B) HSPCs along myeloid pathway (mye-HSPC); (C) HSC-MLP-MEP; (D) HSC. The major findings relevant to aging are summarized in Fig. 1B. Pathways that are significantly increased in the old subjects and throughout all comparisons include: “Formation of senescence associated heterochromatin foci”, “Aging signature HSC”, “DNA damage telomere stress induced senescence”, and “TNFR1 induced NFKB signaling pathway”. “Hypoxia”, “P53 pathway”, “IL2 signaling”, “TNFA signaling via NFKB” were elevated upon enrichment of the primitive HSC population, while other pathways decreased: “Senescence associated secretory phenotype (SASP)”, “DNA methylation”. “Cell cycle” and “DNA repair” are consistently more abundant in young subjects. We analyzed the 512 cells in the HSC cluster with higher resolution and were able to identify a sub-cluster “MPP2A” that is significantly increased in old subjects (4.5% of all HSPCs and 28.5% of HSCs, versus 1.9% and 12.9% respectively in young subjects; p < 0.005). The results are depicted in the highlighted circles in Fig. 1A. This population is quiescent and is highly enriched in “aging signature HSC” in old subjects (p <0.0001). Prominent examples of these signature genes within this list are shown in Fig. 1C). Further along the trajectory, a unique population of B cells “Ly4” was found predominantly in old subjects. The transcriptome profile Ly4 is similar to immature and age-associated B cells (ABCs). This cluster of lymphoid cells may represent the precursors of ABCs described thus far only in peripheral blood, providing another evidence that the aging mechanisms at various differentiation stages are different. Conclusions: In this study, we have defined specific pathways that distinctly separate aging HSCPs from young ones, and shown that the definition of “aging signature” varies with developmental stages. The accumulation a small cluster in the non-primed HSC compartment, with “aging signature HSC” and cell cycle arrest in old subjects indicates that aging of human hematopoiesis is induced by clonal evolution. The abundance of a unique population among lymphoid precursors with expression profiles similar to ABCs demonstrates that lymphoid and myeloid progenitors age differently.

TMIC-63. SINGLE CELL SEQUENCING OF RECURRENT GLIOBLASTOMA DEFINES CELLULAR STATES ASSOCIATED WITH THERAPY RESISTANCE

Abstract Glioblastoma (GBM), an aggressive, primary tumor is composed of multiple clonal populations. Resistance to therapy is a hallmark of GBM, making it a deadly disease. Upon therapeutic stress, tumor heterogeneity and glial stem cell plasticity within the primary tumor contributes to reprogramming of the transcriptional activity. This alters the tumor landscape and generates therapy resistant clonal populations. Single cell RNA sequencing studies of primary GBM have defined various tumor cell states namely, the NPC-like (Neuronal precursor-like), OPC-like (oligodendrocyte precursor-like) and MES-like (mesenchymal-like) cellular states. These studies have also defined the immune components and the normal brain cells associated with primary GBM. Yet, there are no studies focused on the single cell analysis of recurrent GBM. To define the tumor landscape at recurrence, we performed single cell RNA sequencing of 180K cells from excised tumor tissue obtained from 11 patients after failure of therapeutic interventions. A corresponding cohort of 240K cells were sequenced from primary tumors obtained from 13 patients. We performed cellular state analyses of tumor cells in both the primary (~70K) and recurrent (~10K) samples. Primary GBM showed higher fractions of AC-like (0.55), compared to NPC-like (0.10), OPC-like (0.20) and MES-like (0.15). On the other hand, recurrent tissues have very low AC-like (0.09) but elevated NPC-like (0.23), OPC-like (0.40) and MES-like (0.27). This data shows a significant decrease in AC-like cells and corresponding increase in MES-like, OPC-like, NPC-like cells after recurrence. Within this dataset we also have matched longitudinal samples for analysis. We are currently performing trajectory analysis to identify genes and transcription programs activated during tumor recurrence. Experimentally, using known markers for the different cellular states we have also performed immunofluorescence analyses to validate the different cell states altered during recurrence. Our efforts to generate comprehensive recurrent GBM atlas will prove to be valuable to the GBM community.

IMMU-11. NLRP3+ TUMOR-ASSOCIATED MYELOID CELLS ARE ENRICHED IN ASTROCYTIC GLIOMAS AND REPRESENT A THERAPEUTIC VULNERABILITY

Abstract BACKGROUND Tumor-associated myeloid cells (TAMs) are known drivers of immunosuppression and glioma progression. Hence novel approaches for repolarizing TAMs from a protumorigenic to an antitumorigenic phenotype is an active area of investigation. While single-cell RNA sequencing (scRNA-seq) studies have characterized glioma TAMs in great detail, their uniqueness or similarity to TAMs in other tumors remains unknown. Cross-comparison of glioma TAMs with TAM subsets from other tumor types will aid in identifying novel TAM-targeted immunotherapies. METHODS We performed scRNA-seq of ~26K myeloid cells from nine patients with high-grade astrocytoma (n=3) and oligodendroglioma (n=6). The scRNA-seq findings were validated using TCGA RNA-seq data, public scRNA-Seq data, and immunohistochemistry (IHC) analysis. RESULTS The TAM clusters in gliomas predominantly exhibited a microglia-specific gene signature, while macrophage and monocyte clusters were also identified. Utilizing canonical TAM markers, we identified three major TAM subtypes, SPP1+ regulatory, FOLR2+ tissue-resident, and NLRP3+ inflammatory, enriched in astrocytomas. A gene set enrichment analysis of the three TAM subtype genes confirmed their enrichment in astrocytomas (IDHMut-non-codel) (P=0.0005; normalized enrichment score [NES]=1.79) in the TCGA bulk RNAseq data from IDH1-mutant LGGs. Investigation of the three TAM subtype markers in public scRNA-seq data confirmed their enrichment in immune/myeloid clusters across multiple glioma types (gangliogliomas, pilocytic astrocytomas, astrocytomas, diffuse midline gliomas, and glioblastomas). Finally, IHC analysis of 14 astrocytomas demonstrated regional variation in NLRP3 and FOLR2 expression in all cases, with a trend toward an inverse association between NLRP3 and FOLR2 expression. CONCLUSIONS 1) Astrocytomas (IDHMut-non-codel) comprise TAM subsets common across multiple tumor types. 2) The FOLR2+ TAMs predict favorable outcomes, while the NLRP3+ TAMs are immunosuppressive and predict adverse outcomes in breast cancer. Given the inverse expression of NLRP3 and FOLR2 in astrocytomas, we hypothesize that NLRP3 TAMs+ are immunosuppressive and a potential immunotherapeutic target in astrocytomas.

STEM-22. SINGLE-CELL LINEAGE TRACING IN PRIMARY GLIOBLASTOMA REVEALS DISTINCT PROGENITOR SUBTYPES DRIVING INTRATUMORAL HETEROGENEITY

Abstract Glioblastoma (GBM) is the most common type of adult primary brain tumor with nearly a 100% rate of recurrence. One of the key clinical challenges is the vast degree of inter- and intra-tumoral variation in cell types and cell states, which change throughout the course of disease. Single-cell RNA-sequencing (scRNA-seq) studies in GBM samples have illuminated the complex landscape of cell type composition in primary and recurrent GBM, however these studies do not address how this diverse array of cell types developed or which cell populations are principally responsible for generating the heterogeneity and treatment resistance in the tumor. In this study, we implement a scRNA-seq compatible lineage tracing strategy, CellTagging, to track the behavior of individual tumor cells over time. Human GBM samples are dissociated immediately upon resection and infected with the CellTag DNA barcode library. The labeled cells are transplanted onto human cortical organoids, which facilitate the crucial normal-tumor cell interactions that enable GBM growth and invasion. After a period of proliferation, tumor cells are dissociated from the organoid and a fraction are harvested for scRNA-seq for transcriptomic and DNA barcode analysis. The remaining tumor cells are re-transplanted onto new organoids and harvested for scRNA-seq at later timepoints. Clonal relationships are determined by identifying cells with the same DNA barcode profile, and cell type characterizations are assigned via transcriptomic analysis. Transcriptomic profiling demonstrates the existence of multiple subtypes of tumor progenitors across timepoints. Strikingly, clonal analysis indicates that each progenitor subtype gives rise to numerous mature tumor cell types, suggesting multiple distinct lineages that contribute to tumor heterogeneity. Through differential expression analysis, we can predict which progenitors are most likely to drive which cell states. By introducing a treatment paradigm into our lineage tracing protocol, we identify tumor progenitor populations that contribute to treatment resistance.

Organotypic heterogeneity in microvascular endothelial cell responses in sepsis-a molecular treasure trove and pharmacological Gordian knot.

In the last decades, it has become evident that endothelial cells (ECs) in the microvasculature play an important role in the pathophysiology of sepsis-associated multiple organ dysfunction syndrome (MODS). Studies on how ECs orchestrate leukocyte recruitment, control microvascular integrity and permeability, and regulate the haemostatic balance have provided a wealth of knowledge and potential molecular targets that could be considered for pharmacological intervention in sepsis. Yet, this information has not been translated into effective treatments. As MODS affects specific vascular beds, (organotypic) endothelial heterogeneity may be an important contributing factor to this lack of success. On the other hand, given the involvement of ECs in sepsis, this heterogeneity could also be leveraged for therapeutic gain to target specific sites of the vasculature given its full accessibility to drugs. In this review, we describe current knowledge that defines heterogeneity of organ-specific microvascular ECs at the molecular level and elaborate on studies that have reported EC responses across organ systems in sepsis patients and animal models of sepsis. We discuss hypothesis-driven, single-molecule studies that have formed the basis of our understanding of endothelial cell engagement in sepsis pathophysiology, and include recent studies employing high-throughput technologies. The latter deliver comprehensive data sets to describe molecular signatures for organotypic ECs that could lead to new hypotheses and form the foundation for rational pharmacological intervention and biomarker panel development. Particularly results from single cell RNA sequencing and spatial transcriptomics studies are eagerly awaited as they are expected to unveil the full spatiotemporal signature of EC responses to sepsis. With increasing awareness of the existence of distinct sepsis subphenotypes, and the need to develop new drug regimen and companion diagnostics, a better understanding of the molecular pathways exploited by ECs in sepsis pathophysiology will be a cornerstone to halt the detrimental processes that lead to MODS.