BrainScope: interactive visual exploration of the spatial and temporal human brain transcriptome

Spatiotemporal expression patterns of genes coding for plasmalemmal chloride transporters and channels in neurological diseases

Neuronal voltage changes which are dependent on chloride transporters and channels are involved in forming neural functions during early development and maintaining their stability until adulthood. The intracellular chloride concentration maintains a steady state, which is delicately regulated by various genes coding for chloride transporters and channels (GClTC) on the plasmalemma; however, the synergistic effect of these genes in central nervous system disorders remains unclear. In this study, we first defined 10 gene clusters with similar temporal expression patterns, and identified 41 GClTC related to brain developmental process. Then, we found 4 clusters containing 22 GClTC were enriched for the neuronal functions. The GClTC from different clusters presented distinct cell type preferences and anatomical heterogeneity. We also observed strong correlations between clustered genes and diseases, most of which were nervous system disorders. Finally, we found that one of the most well-known GClTC, SLC12A2, had a more profound effect on glial cell-related diseases than on neuron-related diseases, which was in accordance with our observation that SLC12A2 was mainly expressed in oligodendrocytes during brain development. Our findings provide a more comprehensive understanding of the temporal and spatial expression characteristics of GClTC, which can help us understand the complex roles of GClTC in the development of the healthy human brain and the etiology of brain disorders.

Identifying multicellular spatiotemporal organization of cells with SpaceFlow

One major challenge in analyzing spatial transcriptomic datasets is to simultaneously incorporate the cell transcriptome similarity and their spatial locations. Here, we introduce SpaceFlow, which generates spatially-consistent low-dimensional embeddings by incorporating both expression similarity and spatial information using spatially regularized deep graph networks. Based on the embedding, we introduce a pseudo-Spatiotemporal Map that integrates the pseudotime concept with spatial locations of the cells to unravel spatiotemporal patterns of cells. By comparing with multiple existing methods on several spatial transcriptomic datasets at both spot and single-cell resolutions, SpaceFlow is shown to produce a robust domain segmentation and identify biologically meaningful spatiotemporal patterns. Applications of SpaceFlow reveal evolving lineage in heart developmental data and tumor-immune interactions in human breast cancer data. Our study provides a flexible deep learning framework to incorporate spatiotemporal information in analyzing spatial transcriptomic data.

Identification of Cysteine Protease Inhibitor CST2 as a Potential Biomarker for Colorectal Cancer

Additional biomarkers for the development and progression of colorectal cancer (CRC) remain to be identified. Hence, the current study aimed to identify potential diagnostic markers for CRC. Analyses of cysteine protease inhibitor [cystatins (CSTs)] expression in CRC samples and its correlation with cancer stage or survival in patients with CRC demonstrated that CRC tissues had greater CST1 and CST2 mRNA expression compared to noncancerous adjacent tissues, while higher CST2 mRNA expression in CRC tissues was correlated with advanced stages and disease-free survival in patients with CRC, encouraging further exploration on the role of CST2 in CRC. Through an online database search and tissue microarray (TMA), we confirmed that CRC samples had higher CST2 expression compared to noncancerous adjacent tissue or normal colorectal tissues at both the mRNA and protein levels. TMA also revealed that colorectal adenoma, CRC, and metastatic CRC tissues exhibited a significantly increased CST2 protein expression. Accordingly, survival analysis demonstrated that the increase in CST2 protein expression was correlated with shorter overall survival of patients with CRC. Moreover, our results found a significant upregulation of CST2 in multiple cancer tissues. Taken together, these findings suggest the potential role of CST2 as a diagnostic and prognostic biomarker for CRC.

Molecular characterization of the stress network in individuals at risk for schizophrenia

The biological mechanisms underlying inter-individual differences in human stress reactivity remain poorly understood. We aimed to identify the molecular underpinning of aberrant neural stress sensitivity in individuals at risk for schizophrenia. Linking mRNA expression data from the Allen Human Brain Atlas to task-based fMRI revealed 201 differentially expressed genes in cortex-specific brain regions differentially activated by stress in individuals with low (healthy siblings of schizophrenia patients) or high (healthy controls) stress sensitivity. These genes are associated with stress-related psychiatric disorders (e.g. schizophrenia and anxiety) and include markers for specific neuronal populations (e.g. ADCYAP1, GABRB1, SSTR1, and TNFRSF12A), neurotransmitter receptors (e.g. GRIN3A, SSTR1, GABRB1, and HTR1E), and signaling factors that interact with the corticosteroid receptor and hypothalamic-pituitary-adrenal axis (e.g. ADCYAP1, IGSF11, and PKIA). Overall, the identified genes potentially underlie altered stress reactivity in individuals at risk for schizophrenia and other psychiatric disorders and play a role in mounting an adaptive stress response in at-risk individuals, making them potentially druggable targets for stress-related diseases.

Exploring the secrets of brain transcriptional regulation: developing methodologies, recent significant findings, and perspectives

Exploring and revealing the secret of the function of the human brain has been the dream of mankind and science. Delineating brain transcriptional regulation has been extremely challenging, but recent technological advances have facilitated a deeper investigation of molecular processes in the brain. Tracing the molecular regulatory mechanisms of different gene expression profiles in the brain is divergent and has made it possible to connect spatial and temporal variations in gene expression to distributed properties of brain structure and function. Here, we review the molecular diversity of the brain among rodents, non-human primates and humans. We also discuss the molecular mechanism of non-coding DNA/RNA at the transcriptional/post-transcriptional level based on recent technical advances to highlight an improved understanding of the complex transcriptional network in the brain. Spatiotemporal and single-cell transcriptomics have attempted to gain novel insight into the development and evolution of the brain as well as the progression of human diseases. Although it is clear that the field is developing and challenges remain to be resolved, the impressive recent progress provides a solid foundation to better understand the brain and evidence-based recommendations for the diagnosis and treatment of brain diseases.

Cellular resolution in clinical MALDI mass spectrometry imaging: the latest advancements and current challenges

Mass spectrometry (MS) is the workhorse of metabolomics, proteomics and lipidomics. Mass spectrometry imaging (MSI), its extension to spatially resolved analysis of tissues, is a powerful tool for visualizing molecular information within the histological context of tissue. This review summarizes recent developments in MSI and highlights current challenges that remain to achieve molecular imaging at the cellular level of clinical specimens. We focus on matrix-assisted laser desorption/ionization (MALDI)-MSI. We discuss the current status of each of the analysis steps and remaining challenges to reach the desired level of cellular imaging. Currently, analyte delocalization and degradation, matrix crystal size, laser focus restrictions and detector sensitivity are factors that are limiting spatial resolution. New sample preparation devices and laser optic systems are being developed to push the boundaries of these limitations. Furthermore, we review the processing of cellular MSI data and images, and the systematic integration of these data in the light of available algorithms and databases. We discuss roadblocks in the data analysis pipeline and show how technology from other fields can be used to overcome these. Finally, we conclude with curative and community efforts that are needed to enable contextualization of the information obtained.

Visualizing and Interpreting Single-Cell Gene Expression Datasets with Similarity Weighted Nonnegative Embedding

High-throughput single-cell gene expression profiling has enabled the definition of new cell types and developmental trajectories. Visualizing these datasets is crucial to biological interpretation, and a popular method is t-stochastic neighbor embedding (t-SNE), which visualizes local patterns well but distorts global structure, such as distances between clusters. We developed similarity weighted nonnegative embedding (SWNE), which enhances interpretation of datasets by embedding the genes and factors that separate cell states on the visualization alongside the cells and maintains fidelity when visualizing local and global structure for both developmental trajectories and discrete cell types. SWNE uses nonnegative matrix factorization to decompose the gene expression matrix into biologically relevant factors; embeds the cells, genes, and factors in a 2D visualization; and uses a similarity matrix to smooth the embeddings. We demonstrate SWNE on single-cell RNA-seq data from hematopoietic progenitors and human brain cells. SWNE is available as an R package at github.com/yanwu2014/swne.

Analyzing the similarity of samples and genes by MG-PCC algorithm, t-SNE-SS and t-SNE-SG maps

BACKGROUND

For analyzing these gene expression data sets under different samples, clustering and visualizing samples and genes are important methods. However, it is difficult to integrate clustering and visualizing techniques when the similarities of samples and genes are defined by PCC(Person correlation coefficient) measure.

RESULTS

Here, for rare samples of gene expression data sets, we use MG-PCC (mini-groups that are defined by PCC) algorithm to divide them into mini-groups, and use t-SNE-SSP maps to display these mini-groups, where the idea of MG-PCC algorithm is that the nearest neighbors should be in the same mini-groups, t-SNE-SSP map is selected from a series of t-SNE(t-statistic Stochastic Neighbor Embedding) maps of standardized samples, and these t-SNE maps have different perplexity parameter. Moreover, for PCC clusters of mass genes, they are displayed by t-SNE-SGI map, where t-SNE-SGI map is selected from a series of t-SNE maps of standardized genes, and these t-SNE maps have different initialization dimensions. Here, t-SNE-SSP and t-SNE-SGI maps are selected by A-value, where A-value is modeled from areas of clustering projections, and t-SNE-SSP and t-SNE-SGI maps are such t-SNE map that has the smallest A-value.

CONCLUSIONS

From the analysis of cancer gene expression data sets, we demonstrate that MG-PCC algorithm is able to put tumor and normal samples into their respective mini-groups, and t-SNE-SSP(or t-SNE-SGI) maps are able to display the relationships between mini-groups(or PCC clusters) clearly. Furthermore, t-SNE-SS(m)(or t-SNE-SG(n)) maps are able to construct independent tree diagrams of the nearest sample(or gene) neighbors, where each tree diagram is corresponding to a mini-group of samples(or genes).

Natriuretic Peptides in Post-mortem Brain Tissue and Cerebrospinal Fluid of Non-demented Humans and Alzheimer's Disease Patients

Animal studies suggest the involvement of natriuretic peptides (NP) in several brain functions that are known to be disturbed during Alzheimer's disease (AD). However, it remains unclear whether such findings extend to humans. In this study, we aimed to: (1) map the gene expression and localization of NP and their receptors (NPR) in human post-mortem brain tissue; (2) compare the relative amounts of NP and NPR between the brain tissue of AD patients and non-demented controls, and (3) compare the relative amounts of NP between the cerebrospinal fluid (CSF) of AD patients and non-demented controls. Using the publicly available Allen Human Brain Atlas dataset, we mapped the gene expression of NP and NPR in healthy humans. Using immunohistochemistry, we visualized the localization of NP and NPR in the frontal cortex of AD patients (n = 10, mean age 85.8 ± 6.2 years) and non-demented controls (mean age = 80.2 ± 9.1 years). Using Western blotting and ELISA, we quantified the relative amounts of NP and NPR in the brain tissue and CSF of these AD patients and non-demented controls. Our results showed that NP and NPR genes were ubiquitously expressed throughout the brain in healthy humans. NP and NPR were present in various cellular structures including in neurons, astrocyte-like structures, and cerebral vessels in both AD patients and non-demented controls. Furthermore, we found higher amounts of NPR type-A in the brain of AD patients (p = 0.045) and lower amounts of NP type-B in the CSF of AD patients (p = 0.029). In conclusion, this study shows the abundance of NP and NPR in the brain of humans suggesting involvement of NP in various brain functions. In addition, our findings suggest alterations of NP levels in the brain of AD patients. The role of NP in the development and progression of AD remains to be elucidated.

A structural equation model for imaging genetics using spatial transcriptomics

Imaging genetics deals with relationships between genetic variation and imaging variables, often in a disease context. The complex relationships between brain volumes and genetic variants have been explored with both dimension reduction methods and model-based approaches. However, these models usually do not make use of the extensive knowledge of the spatio-anatomical patterns of gene activity. We present a method for integrating genetic markers (single nucleotide polymorphisms) and imaging features, which is based on a causal model and, at the same time, uses the power of dimension reduction. We use structural equation models to find latent variables that explain brain volume changes in a disease context, and which are in turn affected by genetic variants. We make use of publicly available spatial transcriptome data from the Allen Human Brain Atlas to specify the model structure, which reduces noise and improves interpretability. The model is tested in a simulation setting and applied on a case study of the Alzheimer’s Disease Neuroimaging Initiative.

Gene set enrichment analysis and expression pattern exploration implicate an involvement of neurodevelopmental processes in bipolar disorder

BACKGROUND

Bipolar disorder (BD) is a common and highly heritable disorder of mood. Genome-wide association studies (GWAS) have identified several independent susceptibility loci. In order to extract more biological information from GWAS data, multi-locus approaches represent powerful tools since they utilize knowledge about biological processes to integrate functional sets of genes at strongly to moderately associated loci.

METHODS

We conducted gene set enrichment analyses (GSEA) using 2.3 million single-nucleotide polymorphisms, 397 Reactome pathways and 24,025 patients with BD and controls. RNA expression of implicated individual genes and gene sets were examined in post-mortem brains across lifespan.

RESULTS

Two pathways showed a significant enrichment after correction for multiple comparisons in the GSEA: GRB2 events in ERBB2 signaling, for which 6 of 21 genes were BD associated (PFDR = 0.0377), and NCAM signaling for neurite out-growth, for which 11 out of 62 genes were BD associated (PFDR = 0.0451). Most pathway genes showed peaks of RNA co-expression during fetal development and infancy and mapped to neocortical areas and parts of the limbic system.

LIMITATIONS

Pathway associations were technically reproduced by two methods, although they were not formally replicated in independent samples. Gene expression was explored in controls but not in patients.

CONCLUSIONS

Pathway analysis in large GWAS data of BD and follow-up of gene expression patterns in healthy brains provide support for an involvement of neurodevelopmental processes in the etiology of this neuropsychiatric disease. Future studies are required to further evaluate the relevance of the implicated genes on pathway functioning and clinical aspects of BD.

Interactive Visual Exploration of 3D Mass Spectrometry Imaging Data Using Hierarchical Stochastic Neighbor Embedding Reveals Spatiomolecular Structures at Full Data Resolution

Technological advances in mass spectrometry imaging (MSI) have contributed to growing interest in 3D MSI. However, the large size of 3D MSI data sets has made their efficient analysis and visualization and the identification of informative molecular patterns computationally challenging. Hierarchical stochastic neighbor embedding (HSNE), a nonlinear dimensionality reduction technique that aims at finding hierarchical and multiscale representations of large data sets, is a recent development that enables the analysis of millions of data points, with manageable time and memory complexities. We demonstrate that HSNE can be used to analyze large 3D MSI data sets at full mass spectral and spatial resolution. To benchmark the technique as well as demonstrate its broad applicability, we have analyzed a number of publicly available 3D MSI data sets, recorded from various biological systems and spanning different mass-spectrometry ionization techniques. We demonstrate that HSNE is able to rapidly identify regions of interest within these large high-dimensionality data sets as well as aid the identification of molecular ions that characterize these regions of interest; furthermore, through clearly separating measurement artifacts, the HSNE analysis exhibits a degree of robustness to measurement batch effects, spatially correlated noise, and mass spectral misalignment.

Co-expression Patterns between ATN1 and ATXN2 Coincide with Brain Regions Affected in Huntington's Disease

Cytosine-adenine-guanine (CAG) repeat expansions in the coding regions of nine polyglutamine (polyQ) genes (HTT, ATXN1, ATXN2, ATXN3, CACNA1A, ATXN7, ATN1, AR, and TBP) are the cause of several neurodegenerative diseases including Huntington’s disease (HD), six different spinocerebellar ataxias (SCAs), dentatorubral-pallidoluysian atrophy, and spinobulbar muscular atrophy. The expanded CAG repeat length in the causative gene is negatively related to the age-at-onset (AAO) of clinical symptoms. In addition to the expanded CAG repeat length in the causative gene, the normal CAG repeats in the other polyQ genes can affect the AAO, suggesting functional interactions between the polyQ genes. However, there is no detailed assessment of the relationships among polyQ genes in pathologically relevant brain regions. We used gene co-expression analysis to study the functional relationships among polyQ genes in different brain regions using the Allen Human Brain Atlas (AHBA), a spatial map of gene expression in the healthy brain. We constructed co-expression networks for seven anatomical brain structures, as well as a region showing a specific pattern of atrophy in HD patients detected by magnetic resonance imaging (MRI) of the brain. In this HD-associated region, we found that ATN1 and ATXN2 were co-expressed and shared co-expression partners which were enriched for DNA repair genes. We observed a similar co-expression pattern in the frontal lobe, parietal lobe, and striatum in which this relation was most pronounced. Given that the co-expression patterns for these anatomical structures were similar to those for the HD-associated region, our results suggest that their disruption is likely involved in HD pathology. Moreover, ATN1 and ATXN2 also shared many co-expressed genes with HTT, the causative gene of HD, across the brain. Although this triangular relationship among these three polyQ genes may also be dysregulated in other polyQ diseases, stronger co-expression patterns between ATN1 and ATXN2 observed in the HD-associated region, especially in the striatum, may be more specific to HD.