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.

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.

Super-enhancer hijacking drives ectopic expression of hedgehog pathway ligands in meningiomas.

Hedgehog signaling mediates embryologic development of the central nervous system and other tissues and is frequently hijacked by neoplasia to facilitate uncontrolled cellular proliferation. Meningiomas, the most common primary brain tumor, exhibit Hedgehog signaling activation in 6.5% of cases, triggered by recurrent mutations in pathway mediators such as SMO. In this study, we find 35.6% of meningiomas that lack previously known drivers acquired various types of somatic structural variations affecting chromosomes 2q35 and 7q36.3. These cases exhibit ectopic expression of Hedgehog ligands, IHH and SHH, respectively, resulting in Hedgehog signaling activation. Recurrent tandem duplications involving IHH permit de novo chromatin interactions between super-enhancers within DIRC3 and a locus containing IHH. Our work expands the landscape of meningioma molecular drivers and demonstrates enhancer hijacking of Hedgehog ligands as a route to activate this pathway  in neoplasia.

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.

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.

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.

Remote neuronal activity drives glioma infiltration via Sema4f.

The tumor microenvironment (TME) plays an essential role in malignancy and neurons have emerged as a key component of the TME that promotes tumorigenesis across a host of cancers. Recent studies on glioblastoma (GBM) highlight bi-directional signaling between tumors and neurons that propagates a vicious cycle of proliferation, synaptic integration, and brain hyperactivity; however, the identity of neuronal subtypes and tumor subpopulations driving this phenomenon are incompletely understood. Here we show that callosal projection neurons located in the hemisphere contralateral to primary GBM tumors 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 tumors that is enriched for axon guidance genes. High-throughput, in vivo screening of these genes identified Sema4F as a key regulator of tumorigenesis and activity-dependent infiltration. Furthermore, Sema4F promotes the activity-dependent infiltrating population and propagates bi-directional signaling with neurons by remodeling tumor adjacent synapses towards brain network hyperactivity. Collectively, our studies demonstrate that subsets of neurons in locations remote to primary GBM promote malignant progression, while revealing new mechanisms of tumor infiltration that are regulated by neuronal activity.

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.

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.

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 .

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.

PPIL4 is essential for brain angiogenesis and implicated in intracranial aneurysms in humans.

Intracranial aneurysm (IA) rupture leads to subarachnoid hemorrhage, a sudden-onset disease that often causes death or severe disability. Although genome-wide association studies have identified common genetic variants that increase IA risk moderately, the contribution of variants with large effect remains poorly defined. Using whole-exome sequencing, we identified significant enrichment of rare, deleterious mutations in PPIL4, encoding peptidyl-prolyl cis-trans isomerase-like 4, in both familial and index IA cases. Ppil4 depletion in vertebrate models causes intracerebral hemorrhage, defects in cerebrovascular morphology and impaired Wnt signaling. Wild-type, but not IA-mutant, PPIL4 potentiates Wnt signaling by binding JMJD6, a known angiogenesis regulator and Wnt activator. These findings identify a novel PPIL4-dependent Wnt signaling mechanism involved in brain-specific angiogenesis and maintenance of cerebrovascular integrity and implicate PPIL4 gene mutations in the pathogenesis of IA.

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.

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.

Integrated genomic analyses of de novo pathways underlying atypical meningiomas.

This corrects the article DOI: 10.1038/ncomms14433.

ALPK3 gene mutation in a patient with congenital cardiomyopathy and dysmorphic features.

Primary cardiomyopathy is one of the most common inherited cardiac diseases and harbors significant phenotypic and genetic heterogeneity. Because of this, genetic testing has become standard in treatment of this disease group. Indeed, in recent years, next-generation DNA sequencing has found broad applications in medicine, both as a routine diagnostic tool for genetic disorders and as a high-throughput discovery tool for identifying novel disease-causing genes. We describe a male infant with primary dilated cardiomyopathy who was diagnosed using intrauterine echocardiography and found to progress to hypertrophic cardiomyopathy after birth. This proband was born to a nonconsanguineous family with a past history of a male fetus that died because of cardiac abnormalities at 30 wk of gestation. Using whole-exome sequencing, a novel homozygous frameshift mutation (c.2018delC; p.Gln675SerfsX30) in ALPK3 was identified and confirmed with Sanger sequencing. Heterozygous family members were normal with echocardiographic examination. To date, only two studies have reported homozygous pathogenic variants of ALPK3, with a total of seven affected individuals with cardiomyopathy from four unrelated consanguineous families. We include a discussion of the patient's phenotypic features and a review of relevant literature findings.

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.

Longitudinal analysis of treatment-induced genomic alterations in gliomas.

BACKGROUND: Glioblastoma multiforme (GBM) constitutes nearly half of all malignant brain tumors and has a median survival of 15 months. The standard treatment for these lesions includes maximal resection, radiotherapy, and chemotherapy; however, individual tumors display immense variability in their response to these approaches. Genomic techniques such as whole-exome sequencing (WES) provide an opportunity to understand the molecular basis of this variability. METHODS: Here, we report WES-guided treatment of a patient with a primary GBM and two subsequent recurrences, demonstrating the dynamic nature of treatment-induced molecular changes and their implications for clinical decision-making. We also analyze the Yale-Glioma cohort, composed of 110 whole exome- or whole genome-sequenced tumor-normal pairs, to assess the frequency of genomic events found in the presented case. RESULTS: Our longitudinal analysis revealed how the genomic profile evolved under the pressure of therapy. Specifically targeted approaches eradicated treatment-sensitive clones while enriching for resistant ones, generated due to chromothripsis, which we show to be a frequent event in GBMs based on our extended analysis of 110 gliomas in the Yale-Glioma cohort. Despite chromothripsis and the later acquired mismatch-repair deficiency, genomics-guided personalized treatment extended survival to over 5 years. Interestingly, the case displayed a favorable response to immune checkpoint inhibition after acquiring mismatch repair deficiency. CONCLUSIONS: Our study demonstrates the importance of longitudinal genomic profiling to adjust to the dynamic nature of treatment-induced molecular changes to improve the outcomes of precision therapies.

Integrated genomic analyses of de novo pathways underlying atypical meningiomas.

Meningiomas are mostly benign brain tumours, with a potential for becoming atypical or malignant. On the basis of comprehensive genomic, transcriptomic and epigenomic analyses, we compared benign meningiomas to atypical ones. Here, we show that the majority of primary (de novo) atypical meningiomas display loss of NF2, which co-occurs either with genomic instability or recurrent SMARCB1 mutations. These tumours harbour increased H3K27me3 signal and a hypermethylated phenotype, mainly occupying the polycomb repressive complex 2 (PRC2) binding sites in human embryonic stem cells, thereby phenocopying a more primitive cellular state. Consistent with this observation, atypical meningiomas exhibit upregulation of EZH2, the catalytic subunit of the PRC2 complex, as well as the E2F2 and FOXM1 transcriptional networks. Importantly, these primary atypical meningiomas do not harbour TERT promoter mutations, which have been reported in atypical tumours that progressed from benign ones. Our results establish the genomic landscape of primary atypical meningiomas and potential therapeutic targets.