Cell Type-specific Differences Research Articles

Cell type specific transcriptomic differences in depression show similar patterns between males and females but implicate distinct cell types and genes

Major depressive disorder (MDD) is a common, heterogenous, and potentially serious psychiatric illness. Diverse brain cell types have been implicated in MDD etiology. Significant sexual differences exist in MDD clinical presentation and outcome, and recent evidence suggests different molecular bases for male and female MDD. We evaluated over 160,000 nuclei from 71 female and male donors, leveraging new and pre-existing single-nucleus RNA-sequencing data from the dorsolateral prefrontal cortex. Cell type specific transcriptome-wide threshold-free MDD-associated gene expression patterns were similar between the sexes, but significant differentially expressed genes (DEGs) diverged. Among 7 broad cell types and 41 clusters evaluated, microglia and parvalbumin interneurons contributed the most DEGs in females, while deep layer excitatory neurons, astrocytes, and oligodendrocyte precursors were the major contributors in males. Further, the Mic1 cluster with 38% of female DEGs and the ExN10_L46 cluster with 53% of male DEGs, stood out in the meta-analysis of both sexes.

The transcriptional control of the VEGFA-VEGFR1 (FLT1) axis in alternatively polarized murine and human macrophages.

Macrophages significantly contribute to the regulation of vessel formation under physiological and pathological conditions. Although the angiogenesis-regulating role of alternatively polarized macrophages is quite controversial, a growing number of evidence shows that they can participate in the later phases of angiogenesis, including vessel sprouting and remodeling or regression. However, the epigenetic and transcriptional regulatory mechanisms controlling this angiogenesis-modulating program are not fully understood. Here we show that IL-4 can coordinately regulate the VEGFA-VEGFR1 (FLT1) axis via simultaneously inhibiting the proangiogenic Vegfa and inducing the antiangiogenic Flt1 expression in murine bone marrow-derived macrophages, which leads to the attenuated proangiogenic activity of alternatively polarized macrophages. The IL-4-activated STAT6 and IL-4-STAT6 signaling pathway-induced EGR2 transcription factors play a direct role in the transcriptional regulation of the Vegfa-Flt1 axis. We demonstrated that this phenomenon is not restricted to the murine bone marrow-derived macrophages, but can also be observed in different murine tissue-resident macrophages ex vivo and parasites-elicited macrophages in vivo with minor cell type-specific differences. Furthermore, IL-4 exposure can modulate the hypoxic response of genes in both murine and human macrophages leading to a blunted Vegfa/VEGFA and synergistically induced Flt1/FLT1 expression. Our findings establish that the IL-4-activated epigenetic and transcriptional program can determine angiogenesis-regulating properties in alternatively polarized macrophages under normoxic and hypoxic conditions.

The Knowns and Unknowns of Herpesvirus Nuclear Egress.

Nuclear egress of herpesvirus capsids across the intact nuclear envelope is an exceptional vesicle-mediated nucleocytoplasmic translocation resulting in the delivery of herpesvirus capsids into the cytosol. Budding of the (nucleo)capsid at and scission from the inner nuclear membrane (INM) is mediated by the viral nuclear egress complex (NEC) resulting in a transiently enveloped virus particle in the perinuclear space followed by fusion of the primary envelope with the outer nuclear membrane (ONM). The dimeric NEC oligomerizes into a honeycomb-shaped coat underlining the INM to induce membrane curvature and scission. Mutational analyses complemented structural data defining functionally important regions. Questions remain, including where and when the NEC is formed and how membrane curvature is mediated, vesicle formation is regulated, and directionality is secured. The composition of the primary enveloped virion and the machinery mediating fusion of the primary envelope with the ONM is still debated. While NEC-mediated budding apparently follows a highly conserved mechanism, species and/or cell type-specific differences complicate understanding of later steps.

Abstract 4560: Assay-ready tissue-engineered cancer microspheres for high-throughput screening

Abstract For drug discovery, new in vitro cancer models are needed to obtain more translatable study outcomes in a low-cost and high-throughput manner. For this purpose, 3D cancer spheroids have been established as more effective than 2D models. Current commercial techniques, however, rely heavily on self-aggregation of dissociated cells and cannot replicate key features of the native tumor microenvironment, particularly due to a lack of control over extracellular matrix components and heterogeneity in size and aggregate-forming tendencies. Also, current spheroidal techniques are typically limited to one spheroid per well, therefore providing a narrow range of cell numbers per well, disadvantageous for assay development in drug screening. Here, we overcome these challenges by coupling tissue engineering toolsets with microfluidic technologies to create engineered cancer microspheres and sorting desired numbers of microspheres into assay-ready well-plate format. To form the engineered cancer microspheres, MCF7 (non-metastatic) and MDA-MB-231 (metastatic) breast cancer cells were encapsulated within poly(ethylene glycol)-fibrinogen hydrogels using our previously developed microfluidic platform. Highly uniform cancer microspheres (intra and inter-batch coefficient of variation ≤ 5%) with high cell densities (over 20 × 106 cells/ml) were produced rapidly, which is critical for use in drug testing. The microspheres supported the 3D culture of both breast cancer cell lines over at least 14 days in culture. Encapsulated cells displayed cell type-specific differences in morphology, proliferation, metabolic activity, ultrastructure, and overall microsphere size distribution and bulk stiffness. To prepare assay-ready pre-plated microspheres, a COPAS FP flow cytometer was used for its ability to analyze and sort large sample particles such as tumor spheroids and hydrogel cancer microspheres generated in this study. When using a 96-well plate, the sorting rate varied from 2.5 - 6 microspheres per second, depending on the sample concentration. When sorting a desired number of microspheres per well, the accuracy was greater than 95% as verified visually by microscopy. Viability of sorted microspheres was verified 24 hours post-sort. Shipping conditions were established that maintained cell viability for remote use in drug testing. Methods for compound addition by pinning and imaging were tested and optimized. Using these approaches, the microsphere system was shown to be compatible with an automated liquid handling system for administration of drug compounds; MDA-MB-231 microspheres were distributed in 384 well plates and treated with chemotherapeutic drugs. Expected responses were quantitated using CellTiter-Glo® 3D and detected using automated imaging. Overall, our results demonstrate initial applicability for the tissue-engineered cancer microspheres for drug screening. Citation Format: Yuan Tian, Yongwoon Kim, Wen J. Seeto, Shantanu Pradhan, Dmitriy Minond, Rock Pulak, Elizabeth A. Lipke. Assay-ready tissue-engineered cancer microspheres for high-throughput screening. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 4560.

Distinct subpopulations of D1 medium spiny neurons exhibit unique transcriptional responsiveness to cocaine

Drugs of abuse increase extracellular concentrations of dopamine in the nucleus accumbens (NAc), resulting in transcriptional alterations that drive long-lasting cellular and behavioral adaptations. While decades of research have focused on the transcriptional mechanisms by which drugs of abuse influence neuronal physiology and function, few studies have comprehensively defined NAc cell type heterogeneity in transcriptional responses to drugs of abuse. Here, we used single nucleus RNA-seq (snRNA-seq) to characterize the transcriptome of over 39,000 NAc cells from male and female adult Sprague-Dawley rats following acute or repeated cocaine experience. This dataset identified 16 transcriptionally distinct cell populations, including two populations of medium spiny neurons (MSNs) that express the Drd1 dopamine receptor (D1-MSNs). Critically, while both populations expressed classic marker genes of D1-MSNs, only one population exhibited a robust transcriptional response to cocaine. Validation of population-selective transcripts using RNA in situ hybridization revealed distinct spatial compartmentalization of these D1-MSN populations within the NAc. Finally, analysis of published NAc snRNA-seq datasets from non-human primates and humans demonstrated conservation of MSN subtypes across rat and higher order mammals, and further highlighted cell type-specific transcriptional differences across the NAc and broader striatum. These results highlight the utility in using snRNA-seq to characterize both cell type heterogeneity and cell type-specific responses to cocaine and provides a useful resource for cross-species comparisons of NAc cell composition.

Abstract PD4-06: PD4-06 Obesity-associated changes in transcriptomic profile and immune landscape of primary breast cancer revealed by bulk and single-cell gene expression data

Abstract Background: Breast cancer (BC) is one of the cancer types recognized as an obesity-associated disease. Current understandings of molecular mechanisms underlying the BC-obesity connection however largely came from experimental models while systematic investigation of the impact of obesity on BC biology in large patient series is still lacking. The purpose of this study is to discover changes in the transcriptomic profile of primary BC according to patients’ body mass index (BMI). Data and Methods: Bulk and single-cell gene expression data from treatment-naïve primary breast tumors from non-underweight patients were retrieved from the MINDACT trial (NCT00433589; N = 1481) and the pre-treatment cohort of the BioKey trial (NCT03197389, N = 36), respectively. Three categories were considered for BMI: lean, overweight and obese. The main analyses focused on the invasive carcinoma of no special type (NST) estrogen receptor-positive/HER2-negative (ER+/HER2-, N_bulk = 735, N_single-cell(sc) = 10) and NST ER-/HER2- (N_bulk = 118, N_sc = 15) subgroups. The bulk expression data was subjected to differential gene expression analyses according to BMI which was adjusted for menopausal status and tumor grade, then followed by gene set enrichment analyses. Clustering and cluster annotation were performed on the single-cell profiling data before differentially expressed genes according to BMI were identified for each of the present cell types. Results: Obesity-associated differences in the transcriptomic profile of breast tumors, which were subtle but potentially indicative of a biological relationship, were revealed by the bulk data. In both investigated subgroups, tumors from obese patients were shown to be enriched in cell cycle hallmarks. In ER-/HER2- tumors, adiposity further increased MYC signaling. We also observed different obesity-associated changes according to the ER status. Among ER+/HER2- tumors, those from obese patients were enriched in hallmarks related to inflammatory response compared to those from lean patients. In contrast, these hallmarks appeared to be enriched in the ER-/HER2- tumors from lean patients. Our investigation of the single-cell data further revealed shifts in the cell composition of tumor tissue and cell type-specific transcriptomic differences according to BMI which were more pronounced than those detected from the bulk data. ER+/HER2- tumors from obese patients have a higher frequency of immunosuppressive and pro-tumoral cell subpopulations such as dendritic cells (DC) enriched in immunoregulatory molecules (p = .03), LYVE1+ macrophages (p = .02) and myofibroblasts (p = .03) than those from lean patients. Overexpression of Cyclin D1 and CD24 was found in cancer cells in ER+/HER2- tumors from obese patients. A reduction in anti-tumor immune responses was evident with downregulation of multiple interferons in CD8+ and CD4+ T cells as well as B cells. We observed in the ER-/HER2- subgroup increased infiltration of plasmacytoid DC (p = .01), CCL2+ macrophages (p = .01) in tumors from obese versus lean patients, while fibroblasts showed an opposite tendency. Additionally, significant obesity-associated downregulation of major histocompatibility complex (MHC) molecules class I in cancer cells and MHC class II molecules in B cells could be suggestive of deficient antigen presentation and activation of cytotoxic and helper T cells. Conclusion: We highlighted the impact of obesity on the remodeling of tumor and tumor microenvironment which might generally lead to a suppression of anti-tumor immune responses, albeit potentially via diverse axes according to the ER status. Although investigation on a larger cohort is warranted, our current results suggest that obesity-associated transcriptomic changes in BC could be highly cell type-specific, hence we recommend single-cell approaches in addition to spatial multi-omics analysis to further elucidate the interplay between obesity and BC. Citation Format: Ha-Linh Nguyen, Tatjana Geukens, Marion Maetens, Karen Van Baelen, Maxim De Schepper, Coralie Poncet, Mauro Delorenzi, Marjanka K. Schmid, Emiel Rutgers, Laura Van ’t Veer, Martine Piccart, Fatima Cardoso, Giuseppe Viale, Ayse Bassez, Hanne Vos, Patrick Neven, Ines Nevelsteen, Kevin Punie, Hans Wildiers, Giuseppe Floris, Diether Lambrechts, Ann Smeets, Elia Biganzoli, François Richard, Christine Desmedt. PD4-06 Obesity-associated changes in transcriptomic profile and immune landscape of primary breast cancer revealed by bulk and single-cell gene expression data [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr PD4-06.

Single nucleus multiomics identifies ZEB1 and MAFB as candidate regulators of Alzheimer's disease-specific cis-regulatory elements.

Cell type-specific transcriptional differences between brain tissues from donors with Alzheimer's disease (AD) and unaffected controls have been well documented, but few studies have rigorously interrogated the regulatory mechanisms responsible for these alterations. We performed single nucleus multiomics (snRNA-seq plus snATAC-seq) on 105,332 nuclei isolated from cortical tissues from 7 AD and 8 unaffected donors to identify candidate cis-regulatory elements (CREs) involved in AD-associated transcriptional changes. We detected 319,861 significant correlations, or links, between gene expression and cell type-specific transposase accessible regions enriched for active CREs. Among these, 40,831 were unique to AD tissues. Validation experiments confirmed the activity of many regions, including several candidateregulators of APP expression. We identified ZEB1 and MAFB as candidate transcription factorsplaying important roles in AD-specific gene regulation in neurons and microglia, respectively. Microglia links were globally enriched for heritability of AD risk and previously identified active regulatory regions.

Multi–Cell Line Analysis of Lysosomal Proteomes Reveals Unique Features and Novel Lysosomal Proteins

Lysosomes, the main degradative organelles of mammalian cells, play a key role in the regulation of metabolism. It is becoming more and more apparent, that they are highly active, diverse, and involved in a large variety of processes. The essential role of lysosomes is exemplified by the detrimental consequences of their malfunction, which can result in lysosomal storage disorders, neurodegenerative diseases, and cancer. Using lysosome enrichment and mass spectrometry, we investigated the lysosomal proteomes of HEK293, HeLa, HuH-7, SH-SY5Y, MEF, and NIH3T3 cells. We provide evidence on a large scale for cell-type specific differences of lysosomes, showing that levels of distinct lysosomal proteins are highly variable within one cell type, while expression of others is highly conserved across several cell lines. Using differentially stable isotope labeled cells and bimodal distribution analysis, we furthermore identify a high confidence population of lysosomal proteins for each cell line. Multi cell line correlation of these data reveals potential novel lysosomal proteins, and we confirm lysosomal localization for six candidates. All data are available via ProteomeXchange with identifier PXD020600.

Lineage-Mismatched Mitochondrial Replacement in an Inducible Mitochondrial Depletion Model Effectively Restores the Original Proteomic Landscape of Recipient Cells.

In addition to critical roles in bioenergetics, mitochondria are key contributors to the regulation of many other functions in cells, ranging from steroidogenesis to apoptosis. Numerous studies further demonstrate that cell type-specific differences exist in mitochondria, with cells of a given lineage tailoring their endogenous mitochondrial population to suit specific functional needs. These findings, coupled with studies of the therapeutic potential of mitochondrial transplantation, provide a strong impetus to better understand how mitochondria can influence cell function or fate. Here an inducible mitochondrial depletion modelis used to study how cells lacking endogenous mitochondria respond, on a global protein expression level, to transplantation with lineage-mismatched (LM) mitochondria. It is shown that LM mitochondrial transplantation does not alter the proteomic profile in nonmitochondria-depleted recipient cells; however, enforced depletion of endogenous mitochondria results in dramatic changes in the proteomic landscape, which returns to the predepletion state following internalization of LM mitochondria. These data, derived from a cell system that can be rendered free of influence by endogenous mitochondria, indicate that transplantation of mitochondria-even from a source that differs significantly from the recipient cell population, effectively restores a normal proteomic landscape to cells lacking their own mitochondria.

Differential in vitro effects of targeted therapeutics in primary human liver cancer: importance for combined liver cancer

The incidence of primary liver tumors, hepatocellular carcinoma (HCC), intrahepatic cholangiocellular carcinoma (ICC), and combined HCC/ICC (cHCC/CC) is increasing. For ICC, targeted therapy exists only for a small subpopulation of patients, while for HCC, Sorafenib and Lenvatinib are in use. Diagnosis of cHCC/CC is a great challenge and its incidence is underestimated, bearing the risk of unintended non-treatment of ICC. Here, we investigated effects of targeted inhibitors on human ICC cell lines (HUH28, RBE, SSP25), in comparison to extrahepatic (E)CC lines (EGI1, CCC5, TFK1), and HCC/hepatoblastoma cell lines (HEP3B, HUH7, HEPG2). Cells were challenged with: AKT inhibitor MK-2206; multikinase inhibitors Sorafenib, Lenvatinib and Dasatinib; PI3-kinase inhibitors BKM-120, Wortmannin, LY294002, and CAL-101; and mTOR inhibitor Rapamycin. Dosage of the substances was based on the large number of published data of recent years. Proliferation was analyzed daily for four days. All cell lines were highly responsive to MK-2206. Thereby, MK-2206 reduced expression of phospho(p)-AKT in all ICC, ECC, and HCC lines, which mostly corresponded to reduction of p-mTOR, whereas p-ERK1/2 was upregulated in many cases. Lenvatinib showed inhibitory effects on the two HCC cell lines, but not on HEPG2, ICCs and ECCs. Sorafenib inhibited proliferation of all cells, except the ECC line CCC5. However, at reduced dosage, we observed increased cell numbers in some ICC experiments. Dasatinib was highly effective especially in ICC cell lines. Inhibitory effects were observed with all four PI3-kinase inhibitors. However, cell type-specific differences were also evident here. Rapamycin was most effective in the two HCC cell lines. Our studies show that the nine inhibitors differentially target ICC, ECC, and HCC/hepatoblastoma lines. Caution should be taken with Lenvatinib and Sorafenib administration in patients with cHCC/CC as the drugs may have no effects on, or might even stimulate, ICC.

Single-Cell Analysis of Human Clonal Hematopoiesis Identifies Distinct Impact of DNMT3A and TET2 mutations on Hematopoietic Differentiation

Clonal hematopoiesis (CH) increases the risk of hematological malignancies and other adverse outcomes. The most common driver mutations occur in DNMT3A and TET2. However, it is still unclear how mutant DNMT3A and TET2 perturb human stem/progenitor cell homeostasis allowing mutant clones to expand. Studies of human CH have been limited by availability of samples from subjects without co-existing malignancy and high-precision methods for comparing mutant and wild-type (WT) cells in the same individual. To address this gap, we screened 195 bone marrow (BM) samples collected from subjects undergoing hip replacement surgery for CH by targeted DNA sequencing. As DNMT3A and TET2 regulate DNA methylation and gene expression, we employed multiomic single-cell analysis with TARGET-seq to combine single-cell genotyping with transcriptome sequencing on FACS index sorted cells. We developed an improved method, TARGET-seq+, which provides greater sensitivity transcript detection, resulting in a median of 4500 genes/cell with high-sensitivity genotyping (95% of cells genotyped; 7-11% allele dropout rate). We performed TARGET-seq+ on 12,942 BM stem, progenitor and mature cells from 9 CH samples with predicted loss-of-function mutations (4 single mutation cases & 5 multiple mutation cases) and 3 normal controls (Fig. 1). Of these, 87% of cells passed quality controls and were successfully genotyped. Integrated transcriptomic and genotyping analysis showed DNMT3Amut and TET2mut cells occupy a common differentiation trajectory with WT cells. Within DNMT3A-mutant samples, the mutant clone constituted 5-70% of transcriptomic and immunophenotypic long-term hematopoietic stem cells (LT-HSC). There was a further 1.5 to 1.7-fold expansion of DNMT3Amut cells in early erythro-megakaryocytic multipotent progenitors and lympho-myeloid multipotent progenitors (LMPP) compared to LT-HSC. Across most samples, clone size was maintained throughout differentiation with no evidence of differentiation delay or lineage bias. Differential expression analysis revealed subtle cell type-specific differences between WT and mutant cells, with upregulation of megakaryocytic lineage markers in mutant LMPP and downregulation of erythroid lineage markers in mutant erythroid progenitors. TET2mut clones were also identified in LT-HSCs at a clone size of 1-21% but expanded further in multiple downstream progenitors. Maximal expansion (2-12 fold) was seen in the myeloid lineage, but clones also expanded in the erythroid, megakaryocytic and lymphoid lineages (mean clone size 2 to 5-fold larger than in LT-HSC). Comparison of mutant and WT granulocyte-monocyte progenitors (GMP) identified upregulation of myeloid lineage-affiliated transcription factors, (e.g. CEBPD and IRF8) in TET2mut cells. Increased expression of myeloid gene signatures was also seen within TET2mut HSC and LMPP, suggesting myeloid lineage bias arises early in differentiation. Concordantly, TET2mut HSC were shifted towards a more differentiated state in pseudotime and downregulated LT-HSC signatures. In one sample with biallelic TET2 mutations, the TET2 single and double mutant clones were each detected in only 1 out of 109 (0.9%) CD90+CD49f+ LT-HSCs. However, the double mutant clone expanded dramatically in the myeloid lineage, contributing to >50% of cells within GMP, dendritic cell progenitor and neutrophil/monocyte precursor compartments, consistent with frequent detection of biallelic TET2 mutations in myeloid leukemias. Comparison of gene expression between differentiation stage-matched WT cells in CH samples and in age-matched normal controls allowed us to identify the impact of CH clones on WT cells, co-existing within the same BM micro-environment. WT HSCs in DNMT3A and TET2 mutant samples showed enhanced expression of inflammatory, quiescence and chemokine gene signatures compared to WT HSCs in non-CH samples. In summary, we present the first detailed analysis of the impact of DNMT3A and TET2 mutations on hematopoiesis in humans with unperturbed CH. The selective advantage of DNMT3A-mutant clones occurs primarily at the level of HSC and early progenitors, whereas TET2-mutant clones expand not only in HSC but also further throughout differentiation (Fig. 2). There is a non-cell autonomous impact on WT HSCs in CH. Our work provides a foundation for understanding the mechanistic basis of pre-leukemic clonal expansion. Figure 1View largeDownload PPTFigure 1View largeDownload PPT Close modal

Offset Responses in the Auditory Cortex Show Unique History Dependence.

Sensory responses typically vary depending on the recent history of sensory experience. This is essential for processes, including adaptation, efficient coding, and change detection. In the auditory cortex (AC), the short-term history dependence of sound-evoked (onset) responses has been well characterized. Yet many AC neurons also respond to sound terminations, and little is known about the history dependence of these “offset” responses, whether the short-term dynamics of onset and offset responses are correlated, or how these properties are distributed among cell types. Here we presented awake male and female mice with repeating noise burst stimuli while recording single-unit activity from primary AC. We identified parvalbumin and somatostatin interneurons through optotagging, and also separated narrow-spiking from broad-spiking units. We found that offset responses are typically less depressive than onset responses, and this result was robust to a variety of stimulus parameters, controls, measurement types, and selection criteria. Whether a cell's onset response facilitates or depresses does not predict whether its offset response facilitates or depresses. Cell types differed in the dynamics of their onset responses, and in the prevalence, but not the dynamics, of their offset responses. Finally, we clustered cells according to spiking responses and found that response clusters were associated with cell type. Each cluster contained cells of several types, but even within a cluster, cells often showed cell type-specific response dynamics. We conclude that onset and offset responses are differentially influenced by recent sound history, and discuss the implications of this for the encoding of ongoing sound stimuli.SIGNIFICANCE STATEMENT Sensory neuron responses depend on stimulus history. This history dependence is crucial for sensory processing, is precisely controlled at individual synapses and circuits, and is adaptive to the specific requirements of different sensory systems. In the auditory cortex, neurons respond to sound cessation as well as to sound itself, but how history dependence is used along this separate, “offset” information stream is unknown. We show that offset responses are more facilitatory than sound responses, even in neurons where sound responses depress. In contrast to sound onset responses, offset responses are absent in many cells, are relatively homogeneous, and show no cell type-specific differences in history dependence. Offset responses thus show unique response dynamics, suggesting their unique functions.

Abstract P3006: A Human Arterial Cell Atlas

Background: Human vascular diseases are the worldwide leading causes of morbidity and mortality. Nearly all human vascular diseases have arterial segment-specific tropisms despite identical exposures to genetic and environmental risk factors. Understanding the cellular and transcriptomic determinants of arterial identities may hold the key to identifying novel pathophysiology and potential therapies. Methods: To specifically determine arterial site-specific differences independent of inter-individual variation, we have generated a human arterial cellular atlas by simultaneously collecting and analyzing up to 8 arterial sites from multiple healthy transplant donors. We performed single cell transcriptomic analysis on arterial segments to determine the differences in cellular composition and transcriptomic programs. We subsequently integrated human genetic data with cell-type specific transcriptomic differences across vascular beds to identify probable causal cells and causal genes associated with human vascular phenotypes. Results/Conclusions: Single cell transcriptomic analysis of >150,000 cells sequenced at >50,000 reads per cell revealed that the dominant cellular drivers of transcriptomic differences between distinct arterial segments, i.e. determinants of arterial identity, are fibroblasts and smooth muscle cells, not endothelial cells or macrophages. Adult vascular cells transcriptomes from different segments are most influenced by their embryonic origins but not by anatomical proximity. Differentially regulated genes in fibroblast across different vascular beds were particularly enriched for vascular disease associated genetic signals, suggesting a prominent role for these cells in human disease. While the majority of endothelial cells were transcriptionally similar across vascular beds, a rare, previously undescribed, cluster of endothelial cells were identified who expressed segment-specific transcriptomic signatures. Differentially expressed genes in these cells were enriched for vascular disease signals, suggesting a possible role of these rare cells in human disease.

A Survey of the Metabolic Landscape of the Developing Cerebellum at Single-Cell Resolution

The use of cell-culture models to investigate development and disease of the cerebellum is a recent advance, facilitated by the discovery that patterning of precursors is capable of giving rise to cells with specific neuronal identity. Pluripotent stem cell-derived organoids, which exhibit self-organisational characteristics reminiscent of early cerebellar tissue, present a number of challenges including recapitulation of conditions resembling the mature brain. An understanding of the processes driving fetal and postnatal maturation is required to reproduce these conditions in vitro and advance the capability of the system to model adult-onset disease. A key tool for achieving this is single-cell RNA sequencing, which enables visualisation of key transcriptional features of subpopulations comprising tissues. Here, we explore and compare available single-cell RNA sequencing data derived from the developing human cerebellum and its synthetic, in vitro counterpart (stem cell-derived cerebellar organoids). We focus on performing a qualitative assessment of the expression of key metabolic pathway genes, given recent findings exemplifying tissue-specific metabolic activity, including hypoxia and metabolic shifts associated with neuronal expansion. Signatures indicative of known cell type-specific metabolic differences, such as the astrocyte-neuron lactate shuttle and glutamate-glutamine cycle were evident at a transcriptional level. Cerebellar tissue and cerebellar organoids showed a number of behavioural similarities, including HIF1 signalling, which may serve to drive expansion of granule cell progenitors in both settings. We further highlight numerous differences between cultured organoids and native tissue which may provide clarity on the state of metabolic state following differentiation of organoids, providing the future framework to test and further hypotheses regarding promoting maturation. Overall, this analysis provides insight into understanding the state of in vitro models of the cerebellum, a critical factor required for modelling susceptibility of various cell types to cerebellar disease.

Abstract 175: Production of cancer tissue-engineered microspheres for high-throughput screening

Abstract There is a need for new in vitro systems that enable pharmaceutical companies to collect more physiologically-relevant information on drug response in a low-cost and high-throughput manner. For this purpose, three-dimensional (3D) spheroidal models have been established as more effective than two-dimensional models. Current commercial techniques, however, rely heavily on self-aggregation of dissociated cells and are unable to replicate key features of the native tumor microenvironment, particularly due to a lack of control over extracellular matrix components and heterogeneity in shape, size, and aggregate forming tendencies. In this study, we overcome these challenges by coupling tissue engineering toolsets with microfluidics technologies to create engineered cancer microspheres. Specifically, we employ biosynthetic hydrogels composed of conjugated poly(ethylene glycol) (PEG) and fibrinogen protein (PEG-Fb) to create engineered breast and colorectal cancer tissue microspheres for 3D culture, tumorigenic characterization, and examination of potential for high-throughput screening (HTS). MCF7 and MDA-MB-231 cell lines were used to create breast cancer microspheres and the HT29 cell line and cells from a stage II patient-derived xenograft (PDX) were encapsulated to produce colorectal cancer (CRC) microspheres. Using our previously developed microfluidic system, highly uniform cancer microspheres (intra-batch coefficient of variation (CV) ≤ 5%, inter-batch CV < 2%) with high cell densities (>20×106 cells/ml) were produced rapidly, which is critical for use in drug testing. Encapsulated cells maintained high viability and displayed cell type-specific differences in morphology, proliferation, metabolic activity, ultrastructure, and overall microsphere size distribution and bulk stiffness. For PDX CRC microspheres, the percentage of human (70%) and CRC (30%) cells was maintained over time and similar to the original PDX tumor, and the mechanical stiffness also exhibited a similar order of magnitude (103 Pa) to the original tumor. The cancer microsphere system was shown to be compatible with an automated liquid handling system for administration of drug compounds; MDA-MB-231 microspheres were distributed in 384 well plates and treated with staurosporine (1 μM) and doxorubicin (10 μM). Expected responses were quantified using CellTiter-Glo® 3D, demonstrating initial applicability to HTS drug discovery. PDX CRC microspheres were treated with Fluorouracil (5FU) (10 to 500 μM) and displayed a decreasing trend in metabolic activity with increasing drug concentration. Providing a more physiologically relevant tumor microenvironment in a high-throughput and low-cost manner, the PF hydrogel-based cancer microspheres could potentially improve the translational success of drug candidates by providing more accurate in vitro prediction of in vivo drug efficacy. Citation Format: Elizabeth A. Lipke, Wen J. Seeto, Yuan Tian, Mohammadjafar Hashemi, Iman Hassani, Benjamin Anbiah, Nicole L. Habbit, Michael W. Greene, Dmitriy Minond, Shantanu Pradhan. Production of cancer tissue-engineered microspheres for high-throughput screening [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 175.

POS0096 SJÖGREN’S DISEASE AND SYSTEMIC LUPUS ERYTHEMATOSUS DDX6-CXCR5 RISK INTERVALS REVEAL COMMON SNPS WITH FUNCTIONAL SIGNIFICANCE IN IMMUNE AND SALIVARY GLAND CELLS

BackgroundSjögren’s Disease (SjD) and Systemic Lupus Erythematosus (SLE) are autoimmune diseases with several shared characteristics and similar genome-wide significant associations with the DDX6-CXCR5 locus. DDX6 suppresses interferon-stimulated gene expression and CXCR5 regulates T cell functions implicated in autoimmunity.ObjectivesTo identify and characterize functional SNPs in the DDX6-CXCR5 interval.MethodsImmunoChip data from European populations (3785 SLE cases; 1916 SjD cases; 6893 controls) were imputed and SNP-trait associations tested. Bayesian statistics defined a credible SNP set that was refined using bioinformatic analyses (RegulomeDB, Haploreg, ENCODE, promoter capture Hi-C, eQTLs, etc.). Electrophoretic mobility shift assays (EMSAs) and luciferase expression assays were used to test allele-specific SNP function in EBV-transformed B (EBV B) cells, Daudi B cells, Jurkat T cells, THP1 monocytes, and A253 salivary gland cell lines. Chromatin conformation capture with quantitative PCR (3C-qPCR) was used to assess long-range chromatin interactions.ResultsFine mapping of the SjD and SLE associations found similar SNP associations. Bioinformatic analyses identified 5 common SNPs with strong evidence of functionality in immune cell types: rs57494551 in an intron of DDX6, and rs4938572, rs4936443, rs7117261, and rs4938573 in the promoter/enhancer region of DDX6 and CXCR5. EMSAs and luciferase experiments showed cell type-specific differences in protein binding and promoter or enhancer activity, respectively, at each SNP. Risk allele of rs57494551 increased enhancer activity in B cells and A253 cells (p<0.001), but decreased promoter activity in T cells and A253 cells (p<0.01). SNP rs4938572 is an eQTL of DDX6 in T cells, and the risk allele significantly increased protein binding, promoter and enhancer activity in T cells (p<0.01). Risk allele of rs4938572 also increased promoter activity in A253 cells (p<0.001), but had no effect on promoter or enhancer activity in B cells. SNP rs4936443 showed no promoter or enhancer activity in immune cells, but the risk allele showed significant promoter and enhancer (p<0.001) activity in A253 cells. SNP rs7117261 showed decreased enhancer activity in EBV B cells, T cells, and A253 cells (p<0.05) and increased promoter activity in A253 cells (p<0.001). SNP rs4938573 showed decreased promoter activity in EBV B cells, T cell and A253 cells (p<0.05), decreased promoter activity in EBV B cells (p<0.05), and increased enhancer activity in A253 cells (p<0.0001). Overall, A253 cells exhibited more allele-specific effects on promoter and enhancer activity across the five SNPs compared to tested immune cells. In addition to DDX6 and CXCR5, rs57494551 and/or rs4938572 are reported eQTLs for several other genes of interest in the local chromatin regulatory network: IL10RA in T cells, TRAPPC4 in salivary gland and activated macrophages, and long non-coding (lnc)RNA AP002954.1 in T cells and whole blood. 3C-qPCR in EBV B and A253 cells showed that the two regulatory regions carrying rs4938572 or rs57494551 interacted with a region upstream of DDX6 that includes AP002954.1. Hi-C data showed looping between AP002954.1 and the regulatory region carrying rs4938572 and rs57494551 in T cells.ConclusionSjD and SLE share similar genomic architecture across the DDX6-CXCR5 risk interval with several common SNPs showing immune and salivary gland cell type-specific allelic effects on protein binding and/or enhancer/promoter activity. Extensive bioinformatic analyses suggest that the SNPs likely work within the local chromatin regulatory network to regulate cell type-specific expression of several genes on the interval. Ongoing studies will use 3C-qPCR to assess allele-specific chromatin interactions between the SNPs and these genes in different cells types, and CRISPR to determine how the risk alleles alters expression.Disclosure of InterestsMandi M Wiley: None declared, Bhuwan Khatri: None declared, Kandice L Tessneer: None declared, Michelle L Joachims: None declared, Anna M Stolarczyk: None declared, Anna Nagel: None declared, Astrid Rasmussen: None declared, Simon J. Bowman Consultant of: Abbvie, Galapagos, and Novartis in 2020-2021, Lida Radfar: None declared, Roald Omdal: None declared, Marie Wahren-Herlenius: None declared, Blake M Warner: None declared, Torsten Witte: None declared, Roland Jonsson: None declared, Maureen Rischmueller: None declared, Patrick M Gaffney: None declared, Judith A. James: None declared, Lars Ronnblom: None declared, R Hal Scofield: None declared, Xavier Mariette: None declared, Wan Fai Ng: None declared, Kathy Sivils Employee of: current employee of Janssen., Gunnel Nordmark: None declared, Betty Tsao: None declared, Christopher Lessard: None declared

Abstract 107: A Human Arterial Cell Atlas

Background: Human vascular diseases are the worldwide leading causes of morbidity and mortality. Nearly all human vascular diseases have arterial segment-specific tropisms despite identical exposures to genetic and environmental risk factors. Understanding the cellular and transcriptomic determinants of arterial identities may hold the key to identifying novel pathophysiology and potential therapies. Methods: To specifically determine arterial site-specific differences independent of inter-individual variation, we have generated a human arterial cellular atlas by simultaneously collecting and analyzing up to 8 arterial sites from multiple healthy transplant donors. We performed single cell transcriptomic analysis on arterial segments to determine the differences in cellular composition and transcriptomic programs. We subsequently integrated human genetic data with cell-type specific transcriptomic differences across vascular beds to identify probable causal cells and causal genes associated with human vascular phenotypes. Results/Conclusions: Single cell transcriptomic analysis of > 150,000 cells sequenced at > 50,000 reads per cell revealed that the dominant cellular drivers of transcriptomic differences between distinct arterial segments, i.e. determinants of arterial identity, are fibroblasts and smooth muscle cells, not endothelial cells or macrophages. Adult vascular cells from different segments clustered not by anatomical proximity but by embryonic origin. Differentially regulated genes in fibroblast across different vascular beds were particularly enriched for vascular disease associated genetic signals, suggesting a prominent role for these cells in human disease. While the majority of endothelial cells were transcriptionally similar across vascular beds, a rare, previously undescribed, cluster of endothelial cells were identified who expressed segment-specific transcriptomic signatures. Differentially expressed genes in these cells were enriched for vascular disease signals, suggesting a possible role of these rare cells in human disease.

TRAF3 enhances type I interferon (IFN I) receptor signaling in T cells through modulation of the phosphatase PTPN22

Abstract Type I interferons (IFN I) are among the most powerful tools host cells deploy against intracellular pathogens. Their effectiveness is due to both the rapid, directly anti-viral effects of IFN-stimulated gene products, and to the effects of IFN I on responding immune cells. IFN I signaling through its receptor, IFNAR, is tightly regulated at multiple steps in the signaling cascade. Despite the distinct phenotypic changes induced by IFN I in different immune cell types, few studies have addressed the possibility of cell type-specific differences in regulation. In this study, we show that TNFR-associated factor 3 (TRAF3) promotes robust IFNAR signaling specifically in CD4+ T cells by preventing recruitment of the negative regulatory phosphatase PTPN22 to the IFNAR complex. The balance between signals through IFNAR and other differentiating cytokine receptors strongly influences CD4+ T cell differentiation and function during infections. Our work introduces TRAF3 and PTPN22 as newly appreciated regulators key to influencing CD4+ T cell activation by IFN I. Supported by R01 AI123107, T32 AI007260, T32 HL007344, Department of Veterans Affairs

Combined Single Cell Transcriptome and Surface Epitope Profiling Identifies Potential Biomarkers of Psoriatic Arthritis and Facilitates Diagnosis via Machine Learning.

Early diagnosis of psoriatic arthritis (PSA) is important for successful therapeutic intervention but currently remains challenging due, in part, to the scarcity of non-invasive biomarkers. In this study, we performed single cell profiling of transcriptome and cell surface protein expression to compare the peripheral blood immunocyte populations of individuals with PSA, individuals with cutaneous psoriasis (PSO) alone, and healthy individuals. We identified genes and proteins differentially expressed between PSA, PSO, and healthy subjects across 30 immune cell types and observed that some cell types, as well as specific phenotypic subsets of cells, differed in abundance between these cohorts. Cell type-specific gene and protein expression differences between PSA, PSO, and healthy groups, along with 200 previously published genetic risk factors for PSA, were further used to perform machine learning classification, with the best models achieving AUROC ≥ 0.87 when either classifying subjects among the three groups or specifically distinguishing PSA from PSO. Our findings thus expand the repertoire of gene, protein, and cellular biomarkers relevant to PSA and demonstrate the utility of machine learning-based diagnostics for this disease.

Selective Enrichment of Munc13-2 in Presynaptic Active Zones of Hippocampal Pyramidal Cells That Innervate mGluR1α Expressing Interneurons.

Selective distribution of proteins in presynaptic active zones (AZs) is a prerequisite for generating postsynaptic target cell type-specific differences in presynaptic vesicle release probability (Pv) and short-term plasticity, a characteristic feature of cortical pyramidal cells (PCs). In the hippocampus of rodents, somatostatin and mGluR1α expressing interneurons (mGluR1α+ INs) receive small, facilitating excitatory postsynaptic currents (EPSCs) from PCs and express Elfn1 that trans-synaptically recruits mGluR7 into the presynaptic AZ of PC axons. Here we show that Elfn1 also has a role in the selective recruitment of Munc13-2, a synaptic vesicle priming and docking protein, to PC AZs that innervate mGluR1α+ INs. In Elfn1 knock-out mice, unitary EPSCs (uEPSCs) in mGluR1α+ INs have threefold larger amplitudes with less pronounced short-term facilitation, which might be the consequence of the loss of either mGluR7 or Munc13-2 or both. Conditional genetic deletion of Munc13-2 from CA1 PCs results in the loss of Munc13-2, but not mGluR7 from the AZs, and has no effect on the amplitude of uEPSCs and leaves the characteristic short-term facilitation intact at PC to mGluR1α+ IN connection. Our results demonstrate that Munc13-1 alone is capable of imposing low Pv at PC to mGluR1α+ IN synapses and Munc13-2 has yet an unknown role in this synapse.