A high-resolution transcriptomic and spatial atlas of cell types in the whole mouse brain.

The mammalian brain consists of millions to billions of cells that are organized into many cell types with specific spatial distribution patterns and structural and functional properties1-3. Here we report a comprehensive and high-resolution transcriptomic and spatial cell-type atlas for the whole adult mouse brain. The cell-type atlas was created by combining a single-cell RNA-sequencing (scRNA-seq) dataset of around 7 million cells profiled (approximately 4.0 million cells passing quality control), and a spatial transcriptomic dataset of approximately 4.3 million cells using multiplexed error-robust fluorescence in situ hybridization (MERFISH). The atlas is hierarchically organized into 4 nested levels of classification: 34 classes, 338 subclasses, 1,201 supertypes and 5,322 clusters. We present an online platform, Allen Brain Cell Atlas, to visualize the mouse whole-brain cell-type atlas along with the single-cell RNA-sequencing and MERFISH datasets. We systematically analysed the neuronal and non-neuronal cell types across the brain and identified a high degree of correspondence between transcriptomic identity and spatial specificity for each cell type. The results reveal unique features of cell-type organization in different brain regions-in particular, a dichotomy between the dorsal and ventral parts of the brain. The dorsal part contains relatively fewer yet highly divergent neuronal types, whereas the ventral part contains more numerous neuronal types that are more closely related to each other. Our study also uncovered extraordinary diversity and heterogeneity in neurotransmitter and neuropeptide expression and co-expression patterns in different cell types. Finally, we found that transcription factors are major determinants of cell-type classification and identified a combinatorial transcription factor code that defines cell types across all parts of the brain. The whole mouse brain transcriptomic and spatial cell-type atlas establishes a benchmark reference atlas and a foundational resource for integrative investigations of cellular and circuit function, development and evolution of the mammalian brain.

Uveal melanoma immunogenomics predict immunotherapy resistance and susceptibility.

Immune checkpoint inhibition has shown success in treating metastatic cutaneous melanoma but has limited efficacy against metastatic uveal melanoma, a rare variant arising from the immune privileged eye. To better understand this resistance, we comprehensively profile 100 human uveal melanoma metastases using clinicogenomics, transcriptomics, and tumor infiltrating lymphocyte potency assessment. We find that over half of these metastases harbor tumor infiltrating lymphocytes with potent autologous tumor specificity, despite low mutational burden and resistance to prior immunotherapies. However, we observe strikingly low intratumoral T cell receptor clonality within the tumor microenvironment even after prior immunotherapies. To harness these quiescent tumor infiltrating lymphocytes, we develop a transcriptomic biomarker to enable in vivo identification and ex vivo liberation to counter their growth suppression. Finally, we demonstrate that adoptive transfer of these transcriptomically selected tumor infiltrating lymphocytes can promote tumor immunity in patients with metastatic uveal melanoma when other immunotherapies are incapable.

The transcriptomic and spatial organization of telencephalic GABAergic neuronal types.

The telencephalon of the mammalian brain comprises multiple regions and circuit pathways that play adaptive and integrative roles in a variety of brain functions. There is a wide array of GABAergic neurons in the telencephalon; they play a multitude of circuit functions, and dysfunction of these neurons has been implicated in diverse brain disorders. In this study, we conducted a systematic and in-depth analysis of the transcriptomic and spatial organization of GABAergic neuronal types in all regions of the mouse telencephalon and their developmental origins. This was accomplished by utilizing 611,423 single-cell transcriptomes from the comprehensive and high-resolution transcriptomic and spatial cell type atlas for the adult whole mouse brain we have generated, supplemented with an additional single-cell RNA-sequencing dataset containing 99,438 high-quality single-cell transcriptomes collected from the pre- and postnatal developing mouse brain. We present a hierarchically organized adult telencephalic GABAergic neuronal cell type taxonomy of 7 classes, 52 subclasses, 284 supertypes, and 1,051 clusters, as well as a corresponding developmental taxonomy of 450 clusters across different ages. Detailed charting efforts reveal extraordinary complexity where relationships among cell types reflect both spatial locations and developmental origins. Transcriptomically and developmentally related cell types can often be found in distant and diverse brain regions indicating that long-distance migration and dispersion is a common characteristic of nearly all classes of telencephalic GABAergic neurons. Additionally, we find various spatial dimensions of both discrete and continuous variations among related cell types that are correlated with gene expression gradients. Lastly, we find that cortical, striatal and some pallidal GABAergic neurons undergo extensive postnatal diversification, whereas septal and most pallidal GABAergic neuronal types emerge simultaneously during the embryonic stage with limited postnatal diversification. Overall, the telencephalic GABAergic cell type taxonomy can serve as a foundational reference for molecular, structural and functional studies of cell types and circuits by the entire community.

A suite of enhancer AAVs and transgenic mouse lines for genetic access to cortical cell types.

The mammalian cortex is comprised of cells with different morphological, physiological, and molecular properties that can be classified according to shared properties into cell types. Defining the contribution of each cell type to the computational and cognitive processes that are guided by the cortex is essential for understanding its function in health and disease. We use transcriptomic and epigenomic cortical cell type taxonomies from mice and humans to define marker genes and enhancers, and to build genetic tools for cortical cell types. Here, we present a large toolkit for selective targeting of cortical populations, including mouse transgenic lines and recombinant adeno-associated virus (AAV) vectors containing genomic enhancers. We report evaluation of fifteen new transgenic driver lines and over 680 different enhancer AAVs covering all major subclasses of cortical cells, with many achieving a high degree of specificity, comparable with existing transgenic lines. We find that the transgenic lines based on marker genes can provide exceptional specificity and completeness of cell type labeling, but frequently require generation of a triple-transgenic cross for best usability/specificity. On the other hand, enhancer AAVs are easy to screen and use, and can be easily modified to express diverse cargo, such as recombinases. However, their use depends on many factors, such as viral titer and route of administration. The tools reported here as well as the scaled process of tool creation provide an unprecedented resource that should enable diverse experimental strategies towards understanding mammalian cortex and brain function.

Analysis of LMNB1 duplications in autosomal dominant leukodystrophy provides insights into duplication mechanisms and allele-specific expression.

Autosomal dominant leukodystrophy (ADLD) is an adult onset demyelinating disorder that is caused by duplications of the lamin B1 (LMNB1) gene. However, as only a few cases have been analyzed in detail, the mechanisms underlying LMNB1 duplications are unclear. We report the detailed molecular analysis of the largest collection of ADLD families studied, to date. We have identified the minimal duplicated region necessary for the disease, defined all the duplication junctions at the nucleotide level and identified the first inverted LMNB1 duplication. We have demonstrated that the duplications are not recurrent; patients with identical duplications share the same haplotype, likely inherited from a common founder and that the duplications originated from intrachromosomal events. The duplication junction sequences indicated that nonhomologous end joining or replication-based mechanisms such fork stalling and template switching or microhomology-mediated break induced repair are likely to be involved. LMNB1 expression was increased in patients' fibroblasts both at mRNA and protein levels and the three LMNB1 alleles in ADLD patients show equal expression, suggesting that regulatory regions are maintained within the rearranged segment. These results have allowed us to elucidate duplication mechanisms and provide insights into allele-specific LMNB1 expression levels.

Mitotic CDK1 and 4E-BP1 I: Loss of 4E-BP1 serine 82 phosphorylation promotes proliferative polycystic disease and lymphoma in aged or sublethally irradiated mice.

4E-BP1 is a tumor suppressor regulating cap-dependent translation that is in turn controlled by mechanistic target of rapamycin (mTOR) or cyclin-dependent kinase 1 (CDK1) phosphorylation. 4E-BP1 serine 82 (S82) is phosphorylated by CDK1, but not mTOR, and the consequences of this mitosis-specific phosphorylation are unknown. Knock-in mice were generated with a single 4E-BP1 S82 alanine (S82A) substitution leaving other phosphorylation sites intact. S82A mice were fertile and exhibited no gross developmental or behavioral abnormalities, but the homozygotes developed diffuse and severe polycystic liver and kidney disease with aging, and lymphoid malignancies after irradiation. Sublethal irradiation caused immature T-cell lymphoma only in S82A mice while S82A homozygous mice have normal T-cell hematopoiesis before irradiation. Whole genome sequencing identified PTEN mutations in S82A lymphoma and impaired PTEN expression was verified in S82A lymphomas derived cell lines. Our study suggests that the absence of 4E-BP1S82 phosphorylation, a subtle change in 4E-BP1 phosphorylation, might predispose to polycystic proliferative disease and lymphoma under certain stressful circumstances, such as aging and irradiation.

Cell-type specific molecular signatures of aging revealed in a brain-wide transcriptomic cell-type atlas.

Biological aging can be defined as a gradual loss of homeostasis across various aspects of molecular and cellular function. Aging is a complex and dynamic process which influences distinct cell types in a myriad of ways. The cellular architecture of the mammalian brain is heterogeneous and diverse, making it challenging to identify precise areas and cell types of the brain that are more susceptible to aging than others. Here, we present a high-resolution single-cell RNA sequencing dataset containing ~1.2 million high-quality single-cell transcriptomic profiles of brain cells from young adult and aged mice across both sexes, including areas spanning the forebrain, midbrain, and hindbrain. We find age-associated gene expression signatures across nearly all 130+ neuronal and non-neuronal cell subclasses we identified. We detect the greatest gene expression changes in non-neuronal cell types, suggesting that different cell types in the brain vary in their susceptibility to aging. We identify specific, age-enriched clusters within specific glial, vascular, and immune cell types from both cortical and subcortical regions of the brain, and specific gene expression changes associated with cell senescence, inflammation, decrease in new myelination, and decreased vasculature integrity. We also identify genes with expression changes across multiple cell subclasses, pointing to certain mechanisms of aging that may occur across wide regions or broad cell types of the brain. Finally, we discover the greatest gene expression changes in cell types localized to the third ventricle of the hypothalamus, including tanycytes, ependymal cells, and Tbx3+ neurons found in the arcuate nucleus that are part of the neuronal circuits regulating food intake and energy homeostasis. These findings suggest that the area surrounding the third ventricle in the hypothalamus may be a hub for aging in the mouse brain. Overall, we reveal a dynamic landscape of cell-type-specific transcriptomic changes in the brain associated with normal aging that will serve as a foundation for the investigation of functional changes in the aging process and the interaction of aging and diseases.

A high-resolution transcriptomic and spatial atlas of cell types in the whole mouse brain.

The mammalian brain is composed of millions to billions of cells that are organized into numerous cell types with specific spatial distribution patterns and structural and functional properties. An essential step towards understanding brain function is to obtain a parts list, i.e., a catalog of cell types, of the brain. Here, we report a comprehensive and high-resolution transcriptomic and spatial cell type atlas for the whole adult mouse brain. The cell type atlas was created based on the combination of two single-cell-level, whole-brain-scale datasets: a single-cell RNA-sequencing (scRNA-seq) dataset of ~7 million cells profiled, and a spatially resolved transcriptomic dataset of ~4.3 million cells using MERFISH. The atlas is hierarchically organized into five nested levels of classification: 7 divisions, 32 classes, 306 subclasses, 1,045 supertypes and 5,200 clusters. We systematically analyzed the neuronal, non-neuronal, and immature neuronal cell types across the brain and identified a high degree of correspondence between transcriptomic identity and spatial specificity for each cell type. The results reveal unique features of cell type organization in different brain regions, in particular, a dichotomy between the dorsal and ventral parts of the brain: the dorsal part contains relatively fewer yet highly divergent neuronal types, whereas the ventral part contains more numerous neuronal types that are more closely related to each other. We also systematically characterized cell-type specific expression of neurotransmitters, neuropeptides, and transcription factors. The study uncovered extraordinary diversity and heterogeneity in neurotransmitter and neuropeptide expression and co-expression patterns in different cell types across the brain, suggesting they mediate a myriad of modes of intercellular communications. Finally, we found that transcription factors are major determinants of cell type classification in the adult mouse brain and identified a combinatorial transcription factor code that defines cell types across all parts of the brain. The whole-mouse-brain transcriptomic and spatial cell type atlas establishes a benchmark reference atlas and a foundational resource for deep and integrative investigations of cell type and circuit function, development, and evolution of the mammalian brain.

Differential gene expression in peripheral blood mononuclear cells from children immunized with inactivated influenza vaccine.

The human immune response to inactivated influenza vaccine is dynamic and impacted by age and preexisting immunity. Our goal was to identify postvaccination transcriptomic changes in peripheral blood mononuclear cells from children. Blood samples were obtained before and at 3 or 7 days postvaccination with 2016-2017 quadrivalent inactivated influenza vaccine and RNA sequencing was performed. There were 1,466 differentially expressed genes (DEGs) for the Day 0-Day 3 group and 513 DEGs for the Day 0-Day 7 group. Thirty-three genes were common between the two groups. The majority of the transcriptomic changes at Day 3 represented innate inflammation and apoptosis pathways. Day 7 DEGs were characterized by activation of cellular processes, including the regulation of cytoskeleton, junctions, and metabolism, and increased expression of immunoglobulin genes. DEGs at Day 3 were compared between older and younger children revealing increased inflammatory gene expression in the older group. Vaccine history in the year prior to the study was characterized by robust DEGs at Day 3 with decreased phagosome and dendritic cell maturation in those who had been vaccinated in the previous year. PBMC responses to inactivated influenza vaccination in children differed significantly by the timing of sampling, patient age, and vaccine history. These data provide insight into the expected molecular pathways to be temporally altered by influenza vaccination in children.

A Novel Small Molecule Targets Androgen Receptor and Its Splice Variants in Castration-Resistant Prostate Cancer.

Reactivation of androgen receptor (AR) appears to be the major mechanism driving the resistance of castration-resistant prostate cancer (CRPC) to second-generation antiandrogens and involves AR overexpression, AR mutation, and/or expression of AR splice variants lacking ligand-binding domain. There is a need for novel small molecules targeting AR, particularly those also targeting AR splice variants such as ARv7. A high-throughput/high-content screen was previously reported that led to the discovery of a novel lead compound, 2-(((3,5-dimethylisoxazol-4-yl)methyl)thio)-1-(4-(2,3-dimethylphenyl)piperazin-1-yl)ethan-1-one (IMTPPE), capable of inhibiting nuclear AR level and activity in CRPC cells, including those resistant to enzalutamide. A novel analogue of IMTPPE, JJ-450, has been investigated with evidence for its direct and specific inhibition of AR transcriptional activity via a pulldown assay and RNA-sequencing analysis, PSA-based luciferase, qPCR, and chromatin immunoprecipitation assays, and xenograft tumor model 22Rv1. JJ-450 blocks AR recruitment to androgen-responsive elements and suppresses AR target gene expression. JJ-450 also inhibits ARv7 transcriptional activity and its target gene expression. Importantly, JJ-450 suppresses the growth of CRPC tumor xenografts, including ARv7-expressing 22Rv1. Collectively, these findings suggest JJ-450 represents a new class of AR antagonists with therapeutic potential for CRPC, including those resistant to enzalutamide.

Immunogenomic correlates of response to cetuximab monotherapy in head and neck squamous cell carcinoma.

BACKGROUND: Mechanisms of resistance to immune-modulating cancer treatments are poorly understood. Using a novel cohort of patients with head and neck squamous cell carcinoma (HNSCC), we investigated mechanisms of immune escape from epidermal growth factor receptor-specific monoclonal antibody (mAb) therapy. METHODS: HNSCC tumors (n = 20) from a prospective trial of neoadjuvant cetuximab monotherapy underwent whole-exome sequencing. Expression of killer-cell immunoglobulin-like receptor (KIR) and human leukocyte antigen-C (HLA-C) and the effect of KIR blockade were assessed in HNSCC cell lines. RESULTS: Nonresponders to cetuximab had an increased rate of mutations in HLA-C compared to responders and HNSCC tumors (n = 528) in The Cancer Genome Atlas (P < 0.00001). In vitro, cetuximab-activated natural killer (NK) cells induced upregulation of HLA-C on HNSCC cells (P < 0.01) via interferon gamma. Treatment of NK cells with the anti-KIR mAb lirilumab increased killing of HNSCC cells (P < 0.001). CONCLUSIONS: Alterations in HLA-C may provide a mechanism of immune evasion through disruption of NK activation.

IL-33-mediated IL-13 secretion by ST2+ Tregs controls inflammation after lung injury.

Acute respiratory distress syndrome is an often fatal disease that develops after acute lung injury and trauma. How released tissue damage signals, or alarmins, orchestrate early inflammatory events is poorly understood. Herein we reveal that IL-33, an alarmin sequestered in the lung epithelium, is required to limit inflammation after injury due to an unappreciated capacity to mediate Foxp3+ Treg control of local cytokines and myeloid populations. Specifically, Il33-/- mice are more susceptible to lung damage-associated morbidity and mortality that is typified by augmented levels of the proinflammatory cytokines and Ly6Chi monocytes in the bronchoalveolar lavage fluid. Local delivery of IL-33 at the time of injury is protective but requires the presence of Treg cells. IL-33 stimulates both mouse and human Tregs to secrete IL-13. Using Foxp3Cre × Il4/Il13fl/fl mice, we show that Treg expression of IL-13 is required to prevent mortality after acute lung injury by controlling local levels of G-CSF, IL-6, and MCP-1 and inhibiting accumulation of Ly6Chi monocytes. Our study identifies a regulatory mechanism involving IL-33 and Treg secretion of IL-13 in response to tissue damage that is instrumental in limiting local inflammatory responses and may shape the myeloid compartment after lung injury.

EAF2 and p53 Co-Regulate STAT3 Activation in Prostate Cancer.

The tumor suppressor genes EAF2 and p53 are frequently dysregulated in prostate cancers. Recently, we reported that concurrent p53 nuclear staining and EAF2 downregulation were associated with high Gleason score. Combined loss of EAF2 and p53 in a murine model induced prostate tumors, and concurrent knockdown of EAF2 and p53 in prostate cancer cells enhanced proliferation and migration, further suggesting that EAF2 and p53 could functionally interact in the suppression of prostate tumorigenesis. Here, RNA-seq analyses identified differentially regulated genes in response to concurrent knockdown of p53 and EAF2. Several of these genes were associated with the STAT3 signaling pathway, and this was verified by significantly increased p-STAT3 immunostaining in the Eaf2-/-p53-/- mouse prostate. STAT3 knockdown abrogated the stimulation of C4-2 cell proliferation by concurrent knockdown of EAF2 and p53. Furthermore, immunostaining of p-STAT3 was increased in human prostate cancer specimens with EAF2 downregulation and/or p53 nuclear staining. Our findings suggest that simultaneous inactivation of EAF2 and p53 can act to activate STAT3 and drive prostate tumorigenesis.

Research Resource: The Dexamethasone Transcriptome in Hypothalamic Embryonic Neural Stem Cells.

Exposure to excess glucocorticoids during fetal development has long-lasting physiological and behavioral consequences, although the mechanisms are poorly understood. The impact of prenatal glucocorticoids exposure on stress responses in juvenile and adult offspring implicates the developing hypothalamus as a target of adverse prenatal glucocorticoid action. Therefore, primary cultures of hypothalamic neural-progenitor/stem cells (NPSCs) derived from mouse embryos (embryonic day 14.5) were used to identify the glucocorticoid transcriptome in both males and females. NPSCs were treated with vehicle or the synthetic glucocorticoid dexamethasone (dex; 100nM) for 4 hours and total RNA analyzed using RNA-Sequencing. Bioinformatic analysis demonstrated that primary hypothalamic NPSC cultures expressed relatively high levels of a number of genes regulating stem cell proliferation and hypothalamic progenitor function. Interesting, although these cells express glucocorticoid receptors (GRs), only low levels of sex-steroid receptors are expressed, which suggested that sex-specific differentially regulated genes identified are mediated by genetic and not hormonal influences. We also identified known or novel GR-target coding and noncoding genes that are either regulated equivalently in male and female NPSCs or differential responsiveness in one sex. Using gene ontology analysis, the top functional network identified was cell proliferation and using bromodeoxyuridine (BrdU) incorporation observed a reduction in proliferation of hypothalamic NPSCs after dexamethasone treatment. Our studies provide the first characterization and description of glucocorticoid-regulated pathways in male and female embryonically derived hypothalamic NPSCs and identified GR-target genes during hypothalamic development. These findings may provide insight into potential mechanisms responsible for the long-term consequences of fetal glucocorticoid exposure in adulthood.

TCGA Expedition: A Data Acquisition and Management System for TCGA Data.

BACKGROUND: The Cancer Genome Atlas Project (TCGA) is a National Cancer Institute effort to profile at least 500 cases of 20 different tumor types using genomic platforms and to make these data, both raw and processed, available to all researchers. TCGA data are currently over 1.2 Petabyte in size and include whole genome sequence (WGS), whole exome sequence, methylation, RNA expression, proteomic, and clinical datasets. Publicly accessible TCGA data are released through public portals, but many challenges exist in navigating and using data obtained from these sites. We developed TCGA Expedition to support the research community focused on computational methods for cancer research. Data obtained, versioned, and archived using TCGA Expedition supports command line access at high-performance computing facilities as well as some functionality with third party tools. For a subset of TCGA data collected at University of Pittsburgh, we also re-associate TCGA data with de-identified data from the electronic health records. Here we describe the software as well as the architecture of our repository, methods for loading of TCGA data to multiple platforms, and security and regulatory controls that conform to federal best practices. RESULTS: TCGA Expedition software consists of a set of scripts written in Bash, Python and Java that download, extract, harmonize, version and store all TCGA data and metadata. The software generates a versioned, participant- and sample-centered, local TCGA data directory with metadata structures that directly reference the local data files as well as the original data files. The software supports flexible searches of the data via a web portal, user-centric data tracking tools, and data provenance tools. Using this software, we created a collaborative repository, the Pittsburgh Genome Resource Repository (PGRR) that enabled investigators at our institution to work with all TCGA data formats, and to interrogate these data with analysis pipelines, and associated tools. WGS data are especially challenging for individual investigators to use, due to issues with downloading, storage, and processing; having locally accessible WGS BAM files has proven invaluable. CONCLUSION: Our open-source, freely available TCGA Expedition software can be used to create a local collaborative infrastructure for acquiring, managing, and analyzing TCGA data and other large public datasets.

Needs Assessment for Research Use of High-Throughput Sequencing at a Large Academic Medical Center.

Next Generation Sequencing (NGS) methods are driving profound changes in biomedical research, with a growing impact on patient care. Many academic medical centers are evaluating potential models to prepare for the rapid increase in NGS information needs. This study sought to investigate (1) how and where sequencing data is generated and analyzed, (2) research objectives and goals for NGS, (3) workforce capacity and unmet needs, (4) storage capacity and unmet needs, (5) available and anticipated funding resources, and (6) future challenges. As a precursor to informed decision making at our institution, we undertook a systematic needs assessment of investigators using survey methods. We recruited 331 investigators from over 60 departments and divisions at the University of Pittsburgh Schools of Health Sciences and had 140 respondents, or a 42% response rate. Results suggest that both sequencing and analysis bottlenecks currently exist. Significant educational needs were identified, including both investigator-focused needs, such as selection of NGS methods suitable for specific research objectives, and program-focused needs, such as support for training an analytic workforce. The absence of centralized infrastructure was identified as an important institutional gap. Key principles for organizations managing this change were formulated based on the survey responses. This needs assessment provides an in-depth case study which may be useful to other academic medical centers as they identify and plan for future needs.