Inhibitory input directs astrocyte morphogenesis through glial GABABR.

Communication between neurons and glia has an important role in establishing and maintaining higher-order brain function1. Astrocytes are endowed with complex morphologies, placing their peripheral processes in close proximity to neuronal synapses and directly contributing to their regulation of brain circuits2-4. Recent studies have shown that excitatory neuronal activity promotes oligodendrocyte differentiation5-7; whether inhibitory neurotransmission regulates astrocyte morphogenesis during development is unclear. Here we show that inhibitory neuron activity is necessary and sufficient for astrocyte morphogenesis. We found that input from inhibitory neurons functions through astrocytic GABAB receptor (GABABR) and that its deletion in astrocytes results in a loss of morphological complexity across a host of brain regions and disruption of circuit function. Expression of GABABR in developing astrocytes is regulated in a region-specific manner by SOX9 or NFIA and deletion of these transcription factors results in region-specific defects in astrocyte morphogenesis, which is conferred by interactions with transcription factors exhibiting region-restricted patterns of expression. Together, our studies identify input from inhibitory neurons and astrocytic GABABR as universal regulators of morphogenesis, while further revealing a combinatorial code of region-specific transcriptional dependencies for astrocyte development that is intertwined with activity-dependent processes.

Remote neuronal activity drives glioma progression through SEMA4F.

The tumour microenvironment plays an essential role in malignancy, and neurons have emerged as a key component of the tumour microenvironment that promotes tumourigenesis across a host of cancers1,2. Recent studies on glioblastoma (GBM) highlight bidirectional signalling between tumours and neurons that propagates a vicious cycle of proliferation, synaptic integration and brain hyperactivity3-8; however, the identity of neuronal subtypes and tumour subpopulations driving this phenomenon is incompletely understood. Here we show that callosal projection neurons located in the hemisphere contralateral to primary GBM tumours promote progression and widespread infiltration. Using this platform to examine GBM infiltration, we identified an activity-dependent infiltrating population present at the leading edge of mouse and human tumours that is enriched for axon guidance genes. High-throughput, in vivo screening of these genes identified SEMA4F as a key regulator of tumourigenesis and activity-dependent progression. Furthermore, SEMA4F promotes the activity-dependent infiltrating population and propagates bidirectional signalling with neurons by remodelling tumour-adjacent synapses towards brain network hyperactivity. Collectively our studies demonstrate that subsets of neurons in locations remote to primary GBM promote malignant progression, and also show new mechanisms of glioma progression that are regulated by neuronal activity.

Leveraging single-cell sequencing to classify and characterize tumor subgroups in bulk RNA-sequencing data.

PURPOSE: Accurate classification of cancer subgroups is essential for precision medicine, tailoring treatments to individual patients based on their cancer subtypes. In recent years, advances in high-throughput sequencing technologies have enabled the generation of large-scale transcriptomic data from cancer samples. These data have provided opportunities for developing computational methods that can improve cancer subtyping and enable better personalized treatment strategies. METHODS: Here in this study, we evaluated different feature selection schemes in the context of meningioma classification. To integrate interpretable features from the bulk (n = 77 samples) and single-cell profiling (∼ 10 K cells), we developed an algorithm named CLIPPR which combines the top-performing single-cell models, RNA-inferred copy number variation (CNV) signals, and the initial bulk model to create a meta-model. RESULTS: While the scheme relying solely on bulk transcriptomic data showed good classification accuracy, it exhibited confusion between malignant and benign molecular classes in approximately ∼ 8% of meningioma samples. In contrast, models trained on features learned from meningioma single-cell data accurately resolved the sub-groups confused by bulk-transcriptomic data but showed limited overall accuracy. CLIPPR showed superior overall accuracy and resolved benign-malignant confusion as validated on n = 789 bulk meningioma samples gathered from multiple institutions. Finally, we showed the generalizability of our algorithm using our in-house single-cell (∼ 200 K cells) and bulk TCGA glioma data (n = 711 samples). CONCLUSION: Overall, our algorithm CLIPPR synergizes the resolution of single-cell data with the depth of bulk sequencing and enables improved cancer sub-group diagnoses and insights into their biology.

Lunatic Fringe-GFP Marks Lamina-Specific Astrocytes That Regulate Sensory Processing.

Astrocytes are the most abundant glial cell in the brain and perform a wide range of tasks that support neuronal function and circuit activities. There is emerging evidence that astrocytes exhibit molecular and cellular heterogeneity; however, whether distinct subpopulations perform these diverse roles remains poorly defined. Here we show that the Lunatic Fringe-GFP (Lfng-GFP) bacteria artificial chromosome mouse line from both sexes specifically labels astrocyte populations within lamina III and IV of the dorsal spinal cord. Transcriptional profiling of Lfng-GFP+ astrocytes revealed unique molecular profiles, featuring an enriched expression of Notch- and Wnt- pathway components. Leveraging CRE-DOG viral tools, we ablated Lfng-GFP+ astrocytes, which decreased neuronal activity in lamina III and IV and impaired mechanosensation associated with light touch. Together, our findings identify Lfng-GFP+ astrocytes as a unique subpopulation that occupies a distinct anatomic location in the spinal cord and directly contributes to neuronal function and sensory responses.SIGNIFICANCE STATEMENT Astrocytes are the most abundant glial cell in the CNS, and their interactions with neurons are essential for brain function. However, understanding the functional diversity of astrocytes has been hindered because of the lack of reporters that mark subpopulations and genetic tools for accessing them. We discovered that the Lfng-GFP reporter mouse labels a laminae-specific subpopulation of astrocytes in the dorsal spinal cord and that ablation of these astrocytes reduces glutamatergic synapses. Further analysis revealed that these astrocytes have a role in maintaining sensory-processing circuity related to light touch.

CXCR4 expression is associated with proneural-to-mesenchymal transition in glioblastoma.

Glioblastoma (GBM) is the most common primary intracranial malignant tumor and consists of three molecular subtypes: proneural (PN), mesenchymal (MES) and classical (CL). Transition between PN to MES subtypes (PMT) is the glioma analog of the epithelial-mesenchymal transition (EMT) in carcinomas and is associated with resistance to therapy. CXCR4 signaling increases the expression of MES genes in glioma cell lines and promotes EMT in other cancers. RNA sequencing (RNAseq) data of PN GBMs in The Cancer Genome Atlas (TCGA) and secondary high-grade gliomas (HGGs) from an internal cohort were examined for correlation between CXCR4 expression and survival as well as expression of MES markers. Publicly available single-cell RNA sequencing (scRNAseq) data was analyzed for cell type specific CXCR4 expression. These results were validated in a genetic mouse model of PN GBM. Higher CXCR4 expression was associated with significantly reduced survival and increased expression of MES markers in TCGA and internal cohorts. CXCR4 was expressed in immune and tumor cells based on scRNAseq analysis. Higher CXCR4 expression within tumor cells on scRNAseq was associated with increased MES phenotype, suggesting a cell-autonomous effect. In a genetically engineered mouse model, tumors induced with CXCR4 exhibited a mesenchymal phenotype and shortened survival. These results suggest that CXCR4 signaling promotes PMT and shortens survival in GBM and highlights its inhibition as a potential therapeutic strategy.

Integrated clinical genomic analysis reveals xenobiotic metabolic genes are downregulated in meningiomas of current smokers.

INTRODUCTION: Meningiomas are the most common primary intracranial tumor. Recently, various genetic classification systems for meningioma have been described. We sought to identify clinical drivers of different molecular changes in meningioma. As such, clinical and genomic consequences of smoking in patients with meningiomas remain unexplored. METHODS: 88 tumor samples were analyzed in this study. Whole exome sequencing (WES) was used to assess somatic mutation burden. RNA sequencing data was used to identify differentially expressed genes (DEG) and genes sets (GSEA). RESULTS: Fifty-seven patients had no history of smoking, twenty-two were past smokers, and nine were current smokers. The clinical data showed no major differences in natural history across smoking status. WES revealed absence of AKT1 mutation rate in current or past smokers compared to non-smokers (p = 0.046). Current smokers had increased mutation rate in NOTCH2 compared to past and never smokers (p < 0.05). Mutational signature from current and past smokers showed disrupted DNA mismatch repair (cosine-similarity = 0.759 and 0.783). DEG analysis revealed the xenobiotic metabolic genes UGT2A1 and UGT2A2 were both significantly downregulated in current smokers compared to past (Log2FC = - 3.97, padj = 0.0347 and Log2FC = - 4.18, padj = 0.0304) and never smokers (Log2FC = - 3.86, padj = 0.0235 and Log2FC = - 4.20, padj = 0.0149). GSEA analysis of current smokers showed downregulation of xenobiotic metabolism and enrichment for G2M checkpoint, E2F targets, and mitotic spindle compared to past and never smokers (FDR < 25% each). CONCLUSION: In this study, we conducted a comparative analysis of meningioma patients based on their smoking history, examining both their clinical trajectories and molecular changes. Meningiomas from current smokers were more likely to harbor NOTCH2 mutations, and AKT1 mutations were absent in current or past smokers. Moreover, both current and past smokers exhibited a mutational signature associated with DNA mismatch repair. Meningiomas from current smokers demonstrate downregulation of xenobiotic metabolic enzymes UGT2A1 and UGT2A2, which are downregulated in other smoking related cancers. Furthermore, current smokers exhibited downregulation xenobiotic metabolic gene sets, as well as enrichment in gene sets related to mitotic spindle, E2F targets, and G2M checkpoint, which are hallmark pathways involved in cell division and DNA replication control. In aggregate, our results demonstrate novel alterations in meningioma molecular biology in response to systemic carcinogens.

XCVATR: detection and characterization of variant impact on the Embeddings of single -cell and bulk RNA-sequencing samples.

BACKGROUND: RNA-sequencing has become a standard tool for analyzing gene activity in bulk samples and at the single-cell level. By increasing sample sizes and cell counts, this technique can uncover substantial information about cellular transcriptional states. Beyond quantification of gene expression, RNA-seq can be used for detecting variants, including single nucleotide polymorphisms, small insertions/deletions, and larger variants, such as copy number variants. Notably, joint analysis of variants with cellular transcriptional states may provide insights into the impact of mutations, especially for complex and heterogeneous samples. However, this analysis is often challenging due to a prohibitively high number of variants and cells, which are difficult to summarize and visualize. Further, there is a dearth of methods that assess and summarize the association between detected variants and cellular transcriptional states. RESULTS: Here, we introduce XCVATR (eXpressed Clusters of Variant Alleles in Transcriptome pRofiles), a method that identifies variants and detects local enrichment of expressed variants within embedding of samples and cells in single-cell and bulk RNA-seq datasets. XCVATR visualizes local "clumps" of small and large-scale variants and searches for patterns of association between each variant and cellular states, as described by the coordinates of cell embedding, which can be computed independently using any type of distance metrics, such as principal component analysis or t-distributed stochastic neighbor embedding. Through simulations and analysis of real datasets, we demonstrate that XCVATR can detect enrichment of expressed variants and provide insight into the transcriptional states of cells and samples. We next sequenced 2 new single cell RNA-seq tumor samples and applied XCVATR. XCVATR revealed subtle differences in CNV impact on tumors. CONCLUSIONS: XCVATR is publicly available to download from https://github.com/harmancilab/XCVATR .

EpiSAFARI: sensitive detection of valleys in epigenetic signals for enhancing annotations of functional elements.

MOTIVATION: Functional genomics experiments generate genomewide signal profiles that are dense information sources for annotating the regulatory elements. These profiles measure epigenetic activity at the nucleotide resolution and they exhibit distinctive patterns as they fluctuate along the genome. Most notable of these patterns are the valley patterns that are prevalently observed in assays such as ChIP Sequencing and bisulfite sequencing. The genomic positions of valleys pinpoint locations of cis-regulatory elements such as enhancers and insulators. Systematic identification of the valleys provides novel information for delineating the annotation of regulatory elements. Nevertheless, the valleys are not reported by majority of the analysis pipelines. RESULTS: We describe EpiSAFARI, a computational method for sensitive detection of valleys from diverse types of epigenetic profiles. EpiSAFARI employs a novel smoothing method for decreasing noise in signal profiles and accounts for technical factors such as sparse signals, mappability and nucleotide content. In performance comparisons, EpiSAFARI performs favorably in terms of accuracy. The histone modification valleys detected by EpiSAFARI exhibit high conservation, transcription factor binding and they are enriched in nascent transcription. In addition, the large clusters of histone valleys are found to be enriched at the promoters of the developmentally associated genes. Differential histone valleys exhibit concordance with differential DNase signal at cell line specific valleys. DNA methylation valleys exhibit elevated conservation and high transcription factor binding. Specifically, we observed enriched binding of transcription factors associated with chromatin structure around methyl-valleys. AVAILABILITY AND IMPLEMENTATION: EpiSAFARI is publicly available at https://github.com/harmancilab/EpiSAFARI. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

METAP1 mutation is a novel candidate for autosomal recessive intellectual disability.

Intellectual disability (ID) is a genetic and clinically heterogeneous common disease and underlying molecular pathogenesis can frequently not be identified by whole-exome/genome testing. Here, we report four siblings born to a consanguineous union who presented with intellectual disability and discuss the METAP1 pathway as a novel etiology of ID. Genomic analyses demonstrated that patients harbor a novel homozygous nonsense mutation in the gene METAP1. METAP1 codes for methionine aminopeptidase 1 (MetAP1) which oversees the co-translational excision of the first methionine remnants in eukaryotes. The loss-of-function mutations to this gene may result in a defect in the translation of many essential proteins within a cell. Improper neuronal function resulting from this loss of essential proteins could lead to neurologic impairment and ID.

Identification of novel fusion transcripts in meningioma.

INTRODUCTION: Meningiomas are the most common primary intracranial tumor. Recent next generation sequencing analyses have elaborated the molecular drivers of this disease. We aimed to identify and characterize novel fusion genes in meningiomas. METHODS: We performed a secondary analysis of our RNA sequencing data of 145 primary meningioma from 140 patients to detect fusion genes. Semi-quantitative rt-PCR was performed to confirm transcription of the fusion genes in the original tumors. Whole exome sequencing was performed to identify copy number variations within each tumor sample. Comparative RNA seq analysis was performed to assess the clonality of the fusion constructs within the tumor. RESULTS: We detected six fusion events (NOTCH3-SETBP1, NF2-SPATA13, SLC6A3-AGBL3, PHF19-FOXP2 in two patients, and ITPK1-FBP2) in five out of 145 tumor samples. All but one event (NF2-SPATA13) led to extremely short reading frames, making these events de facto null alleles. Three of the five patients had a history of childhood radiation. Four out of six fusion events were detected in expression type C tumors, which represent the most aggressive meningioma. We validated the presence of the RNA transcripts in the tumor tissue by semi-quantitative RT PCR. All but the two PHF19-FOXP2 fusions demonstrated high degrees of clonality. CONCLUSIONS: Fusion genes occur infrequently in meningiomas and are more likely to be found in tumors with greater degree of genomic instability (expression type C) or in patients with history of cranial irradiation.

Exome analysis of the evolutionary path of hepatocellular adenoma-carcinoma transition, vascular invasion and brain dissemination.

Hepatocellular adenoma (HCA) is a rare benign liver tumor, predominantly seen in young women. Its major complications are malignant transformation, spontaneous hemorrhage, and rupture. We describe a case of a young female with no underlying liver disease who presented with acute abdominal pain and was found to have a 17cm heterogeneous mass in the left lobe of the liver. She underwent left hepatectomy and pathology revealed a 14cm moderately differentiated hepatocellular carcinoma (HCC) arising in a shell of a HCA. At that time, vascular invasion was already present. She rapidly developed recurrent multifocal hepatic lesions and subsequent spread to the brain, leading to her death 18months after surgery. To investigate the underlying genetic events occurring during hepatocellular adenoma-carcinoma transition and extra-hepatic dissemination, we performed whole exome sequencing of DNA isolated from peripheral blood leucocytes, HCA, HCC, tumor thrombus and brain metastasis. Our data show a step-wise addition of somatic mutations and copy number variations with disease progression, suggesting a linear tumor evolution, which is supported by clonality analysis. Specifically, using a model based clustering of somatic mutations, one single founding clone arising in the HCA, which included catenin beta 1 (CTNNB1) and IL6ST driver mutations, was identified and displayed an increasing clonality rate in HCC, tumor thrombus and brain metastasis. Our data highlight the feasibility of performing whole exome capture, sequencing and analysis using formalin-fixed paraffin-embedded (FFPE) samples, and we describe the first genomic longitudinal study of hepatocellular adenoma-carcinoma transition, vascular invasion and brain metastasis with detailed clinicopathologic annotation.

A patient with a novel homozygous missense mutation in FTO and concomitant nonsense mutation in CETP.

The fat mass and obesity associated (FTO) gene has previously been associated with a variety of diseases and conditions, notably obesity, acute coronary syndrome and metabolic syndrome. Reports describing mutations in FTO as well as in FTO animal models have further demonstrated a role for FTO in the development of the brain and other organs. Here, we describe a patient born of consanguineous union who presented with microcephaly, developmental delay, behavioral abnormalities, dysmorphic facial features, hypotonia and other various phenotypic abnormalities. Whole-exome sequencing revealed a novel homozygous missense mutation in FTO and a nonsense mutation in the cholesteryl ester transfer protein (CETP). Exome copy number variation analysis revealed no disease-causing large duplications or deletions within coding regions. Patient's, her parents' and non-related control' fibroblasts were analyzed for morphologic defects, abnormal proliferation, apoptosis and transcriptome profile. We have shown that FTO is located in the nucleus of cells from each tested sample. Western blot analysis demonstrated no changes in patient FTO. Quantitative (qPCR) analysis revealed slightly decreased levels of FTO expression in patient cells compared with controls. No morphological or proliferation differences between the patient and control fibroblasts were observed. There is still much to be learned about the molecular mechanisms by which mutations in FTO contribute to such severe phenotypes.

Paediatric hepatocellular carcinoma due to somatic CTNNB1 and NFE2L2 mutations in the setting of inherited bi-allelic ABCB11 mutations.

Hepatocellular carcinoma (HCC) rarely occurs in childhood. We describe a patient with new onset of pruritus at 8 months of age who at 17 months of age was found to have a 2.5 cm HCC. To delineate the possible genetic basis of this tumour, we performed whole exome sequencing (WES) of the germline DNA and identified two novel predictably deleterious missense mutations in ABCB11, encoding bile salt export pump (BSEP), confirmed in the parental DNA as bi-allelic and inherited. Although inherited ABCB11 mutations have previously been linked to HCC in a small number of cases, the molecular mechanisms of hepatocellular carcinogenesis in ABCB11 disease are unknown. WES of the HCC tissue uncovered somatic driver mutations in the beta-catenin (CTNNB1) and nuclear-factor-erythroid-2-related-factor-2 (NFE2L2) genes. Moreover, clonality analysis predicted that the CTNNB1 mutation was clonal and occurred earlier during carcinogenesis, whereas the NFE2L2 mutation was acquired later. Interestingly, background liver parenchyma showed no inflammation or fibrosis and BSEP expression was preserved. This is the first study to identify somatic CTNNB1 and NFE2L2 mutations in early childhood arisen in the setting of inherited bi-allelic ABCB11 mutations. Rapid WES analysis expedited this child's diagnosis and treatment, and likely improved her prognosis.

Integrated electrophysiological and genomic profiles of single cells reveal spiking tumor cells in human glioma.

Prior studies have described the complex interplay that exists between glioma cells and neurons, however, the electrophysiological properties endogenous to tumor cells remain obscure. To address this, we employed Patch-sequencing on human glioma specimens and found that one third of patched cells in IDH mutant (IDH mut ) tumors demonstrate properties of both neurons and glia by firing single, short action potentials. To define these hybrid cells (HCs) and discern if they are tumor in origin, we developed a computational tool, Single Cell Rule Association Mining (SCRAM), to annotate each cell individually. SCRAM revealed that HCs represent tumor and non-tumor cells that feature GABAergic neuron and oligodendrocyte precursor cell signatures. These studies are the first to characterize the combined electrophysiological and molecular properties of human glioma cells and describe a new cell type in human glioma with unique electrophysiological and transcriptomic properties that are likely also present in the non-tumor mammalian brain.

Social deprivation induces astrocytic TRPA1-GABA suppression of hippocampal circuits.

Social experience is essential for the development and maintenance of higher-order brain function. Social deprivation results in a host of cognitive deficits, and cellular studies have largely focused on associated neuronal dysregulation; how astrocyte function is impacted by social deprivation is unknown. Here, we show that hippocampal astrocytes from juvenile mice subjected to social isolation exhibit increased Ca2+ activity and global changes in gene expression. We found that the Ca2+ channel TRPA1 is upregulated in astrocytes after social deprivation and astrocyte-specific deletion of TRPA1 reverses the physiological and cognitive deficits associated with social deprivation. Mechanistically, TRPA1 inhibition of hippocampal circuits is mediated by a parallel increase of astrocytic production and release of the inhibitory neurotransmitter GABA after social deprivation. Collectively, our studies reveal how astrocyte function is tuned to social experience and identifies a social-context-specific mechanism by which astrocytic TRPA1 and GABA coordinately suppress hippocampal circuit function.

Inhibitory input directs astrocyte morphogenesis through glial GABA B R.

Communication between neurons and glia plays an important role in establishing and maintaining higher order brain function. Astrocytes are endowed with complex morphologies which places their peripheral processes in close proximity to neuronal synapses and directly contributes to their regulation of brain circuits. Recent studies have shown that excitatory neuronal activity promotes oligodendrocyte differentiation; whether inhibitory neurotransmission regulates astrocyte morphogenesis during development is unknown. Here we show that inhibitory neuron activity is necessary and sufficient for astrocyte morphogenesis. We found that input from inhibitory neurons functions through astrocytic GABA B R and that its deletion in astrocytes results in a loss of morphological complexity across a host of brain regions and disruption of circuit function. Expression of GABA B R in developing astrocytes is regulated in a region-specific manner by SOX9 or NFIA and deletion of these transcription factors results in region-specific defects in astrocyte morphogenesis, which is conferred by interactions with transcription factors exhibiting region-restricted patterns of expression. Together our studies identify input from inhibitory neurons and astrocytic GABA B R as universal regulators of morphogenesis, while further revealing a combinatorial code of region-specific transcriptional dependencies for astrocyte development that is intertwined with activity-dependent processes.

Glioma epileptiform activity and progression are driven by IGSF3-mediated potassium dysregulation.

Seizures are a frequent pathophysiological feature of malignant glioma. Recent studies implicate peritumoral synaptic dysregulation as a driver of brain hyperactivity and tumor progression; however, the molecular mechanisms that govern these phenomena remain elusive. Using scRNA-seq and intraoperative patient ECoG recordings, we show that tumors from seizure patients are enriched for gene signatures regulating synapse formation. Employing a human-to-mouse in vivo functionalization pipeline to screen these genes, we identify IGSF3 as a mediator of glioma progression and dysregulated neural circuitry that manifests as spreading depolarization (SD). Mechanistically, we discover that IGSF3 interacts with Kir4.1 to suppress potassium buffering and found that seizure patients exhibit reduced expression of potassium handlers in proliferating tumor cells. In vivo imaging reveals that dysregulated synaptic activity emanates from the tumor-neuron interface, which we confirm in patients. Our studies reveal that tumor progression and seizures are enabled by ion dyshomeostasis and identify SD as a driver of disease.