RNA-seq analysis of amygdala tissue reveals characteristic expression profiles in schizophrenia

Connecting genomic results for psychiatric disorders to human brain cell types and regions reveals convergence with functional connectivity

Identifying cell types and brain regions critical for psychiatric disorders and brain traits is essential for targeted neurobiological research. By integrating genomic insights from genome-wide association studies with a comprehensive single-cell transcriptomic atlas of the adult human brain, we prioritized specific neuronal clusters significantly enriched for the SNP-heritabilities for schizophrenia, bipolar disorder, and major depressive disorder along with intelligence, education, and neuroticism. Extrapolation of cell-type results to brain regions reveals the whole-brain impact of schizophrenia genetic risk, with subregions in the hippocampus and amygdala exhibiting the most significant enrichment of SNP-heritability. Using functional MRI connectivity, we further confirmed the significance of the central and lateral amygdala, hippocampal body, and prefrontal cortex in distinguishing schizophrenia cases from controls. Our findings underscore the value of single-cell transcriptomics in understanding the polygenicity of psychiatric disorders and suggest a promising alignment of genomic, transcriptomic, and brain imaging modalities for identifying common biological targets.

MK-801-exposure induces increased translation efficiency and mRNA hyperacetylation of Grin2a in the mouse prefrontal cortex

Acute exposure to MK-801, the non-competitive N-methyl-D-aspartate (NMDA) receptor antagonist, induces schizophrenia-like behavioural changes in juvenile male mice. However, the effects of acute MK-801 exposure on brain gene expression at the translation level remain unclear. Here, we conducted ribosome profiling analysis on the prefrontal cortex (PFC) of acute MK-801-exposed juvenile male mice. We found 357 differentially translated genes, with the N4-acetylcytidine (ac4C) consensus motif enriched in the transcripts with increased translation efficiency. Acetylated RNA immunoprecipitation sequencing revealed 148 differentially acetylated peaks, of which 121 were hyperacetylated, and 27 were hypoacetylated. Genes harbouring these peaks were enriched in pathways related to axon guidance, Hedgehog signalling pathway, neuron differentiation, and memory. Grin2a encodes an NMDA receptor subunit NMDAR2A, and its human orthologue is a strong susceptibility gene for schizophrenia. Grin2a mRNA was hyperacetylated and exhibited significantly increased translation efficiency. NMDAR2A protein level was increased in MK-801-exposed PFC. Pretreatment of Remodelin, an inhibitor of N-acetyltransferase 10, returned the NMDAR2A protein levels to normal and partially reversed schizophrenia-like behaviours of MK-801-exposed mice, shedding light on the possible role of mRNA acetylation in the aetiology of schizophrenia.

Semaphorin 6D tunes amygdalar circuits for emotional, metabolic, and inflammatory outputs

Regulated neural-metabolic-inflammatory responses are essential for maintaining physiological homeostasis. However, the molecular machinery that coordinates neural, metabolic, and inflammatory responses is largely unknown. Here, we show that semaphorin 6D (SEMA6D) coordinates anxiogenic, metabolic, and inflammatory outputs from the amygdala by maintaining synaptic homeostasis. Using genome-wide approaches, we identify SEMA6D as a pleiotropic gene for both psychiatric and metabolic traits in human. Sema6d deficiency increases anxiety in mice. When fed a high-fat diet, Sema6d-/- mice display attenuated obesity and enhanced myelopoiesis compared with control mice due to higher sympathetic activity via the β3-adrenergic receptor. Genetic manipulation and spatial and single-nucleus transcriptomics reveal that SEMA6D in amygdalar interneurons is responsible for regulating anxiogenic and autonomic responses. Mechanistically, SEMA6D is required for synaptic maturation and γ-aminobutyric acid transmission. These results demonstrate that SEMA6D is important for the normal functioning of the neural circuits in the amygdala, coupling emotional, metabolic, and inflammatory responses.

The impact of sex on gene expression in the brain of schizophrenic patients: a systematic review and meta-analysis of transcriptomic studies

BACKGROUND

Schizophrenia is a severe neuropsychiatric disorder characterized by altered perception, mood, and behavior that profoundly impacts patients and society despite its relatively low prevalence. Sex-based differences have been described in schizophrenia epidemiology, symptomatology and outcomes. Different studies explored the impact of schizophrenia in the brain transcriptome, however we lack a consensus transcriptomic profile that considers sex and differentiates specific cerebral regions.

METHODS

We performed a systematic review on bulk RNA-sequencing studies of post-mortem brain samples. Then, we fulfilled differential expression analysis on each study and summarized their results with regions-specific meta-analyses (prefrontal cortex and hippocampus) and a global all-studies meta-analysis. Finally, we used the consensus transcriptomic profiles to functionally characterize the impact of schizophrenia in males and females by protein-protein interaction networks, enriched biological processes and dysregulated transcription factors.

RESULTS

We discovered the sex-based dysregulation of 265 genes in the prefrontal cortex, 1.414 genes in the hippocampus and 66 genes in the all-studies meta-analyses. The functional characterization of these gene sets unveiled increased processes related to immune response functions in the prefrontal cortex in male and the hippocampus in female schizophrenia patients and the overexpression of genes related to neurotransmission and synapses in the prefrontal cortex of female schizophrenia patients. Considering a meta-analysis of all brain regions available, we encountered the relative overexpression of genes related to synaptic plasticity and transmission in females and the overexpression of genes involved in organizing genetic information and protein folding in male schizophrenia patients. The protein-protein interaction networks and transcription factors activity analyses supported these sex-based profiles.

CONCLUSIONS

Our results report multiple sex-based transcriptomic alterations in specific brain regions of schizophrenia patients, which provides new insight into the role of sex in schizophrenia. Moreover, we unveil a partial overlapping of inflammatory processes in the prefrontal cortex of males and the hippocampus of females.

Validation of Enhancer Regions in Primary Human Neural Progenitor Cells using Capture STARR-seq

Genome-wide association studies (GWAS) and expression analyses implicate noncoding regulatory regions as harboring risk factors for psychiatric disease, but functional characterization of these regions remains limited. We performed capture STARR-sequencing of over 78,000 candidate regions to identify active enhancers in primary human neural progenitor cells (phNPCs). We selected candidate regions by integrating data from NPCs, prefrontal cortex, developmental timepoints, and GWAS. Over 8,000 regions demonstrated enhancer activity in the phNPCs, and we linked these regions to over 2,200 predicted target genes. These genes are involved in neuronal and psychiatric disease-associated pathways, including dopaminergic synapse, axon guidance, and schizophrenia. We functionally validated a subset of these enhancers using mutation STARR-sequencing and CRISPR deletions, demonstrating the effects of genetic variation on enhancer activity and enhancer deletion on gene expression. Overall, we identified thousands of highly active enhancers and functionally validated a subset of these enhancers, improving our understanding of regulatory networks underlying brain function and disease.

Schizophrenia endothelial cells exhibit higher permeability and altered angiogenesis patterns in patient-derived organoids

Schizophrenia (SCZ) is a complex neurodevelopmental disorder characterized by the manifestation of psychiatric symptoms in early adulthood. While many research avenues into the origins of SCZ during brain development have been explored, the contribution of endothelial/vascular dysfunction to the disease remains largely elusive. To model the neuropathology of SCZ during early critical periods of brain development, we utilized patient-derived induced pluripotent stem cells (iPSCs) to generate 3D cerebral organoids and define cell-specific signatures of disease. Single-cell RNA sequencing revealed that while SCZ organoids were similar in their macromolecular diversity to organoids generated from healthy controls (CTRL), SCZ organoids exhibited a higher percentage of endothelial cells when normalized to total cell numbers. Additionally, when compared to CTRL, differential gene expression analysis revealed a significant enrichment in genes that function in vessel formation, vascular regulation, and inflammatory response in SCZ endothelial cells. In line with these findings, data from 23 donors demonstrated that PECAM1+ microvascular vessel-like structures were increased in length and number in SCZ organoids in comparison to CTRL organoids. Furthermore, we report that patient-derived endothelial cells displayed higher paracellular permeability, implicating elevated vascular activity. Collectively, our data identified altered gene expression patterns, vessel-like structural changes, and enhanced permeability of endothelial cells in patient-derived models of SCZ. Hence, brain microvascular cells could play a role in the etiology of SCZ by modulating the permeability of the developing blood brain barrier (BBB).

Connecting genomic results for psychiatric disorders to human brain cell types and regions reveals convergence with functional connectivity

Understanding the temporal and spatial brain locations etiological for psychiatric disorders is essential for targeted neurobiological research. Integration of genomic insights from genome-wide association studies with single-cell transcriptomics is a powerful approach although past efforts have necessarily relied on mouse atlases. Leveraging a comprehensive atlas of the adult human brain, we prioritized cell types via the enrichment of SNP-heritabilities for brain diseases, disorders, and traits, progressing from individual cell types to brain regions. Our findings highlight specific neuronal clusters significantly enriched for the SNP-heritabilities for schizophrenia, bipolar disorder, and major depressive disorder along with intelligence, education, and neuroticism. Extrapolation of cell-type results to brain regions reveals important patterns for schizophrenia with distinct subregions in the hippocampus and amygdala exhibiting the highest significance. Cerebral cortical regions display similar enrichments despite the known prefrontal dysfunction in those with schizophrenia highlighting the importance of subcortical connectivity. Using functional MRI connectivity from cases with schizophrenia and neurotypical controls, we identified brain networks that distinguished cases from controls that also confirmed involvement of the central and lateral amygdala, hippocampal body, and prefrontal cortex. Our findings underscore the value of single-cell transcriptomics in decoding the polygenicity of psychiatric disorders and offer a promising convergence of genomic, transcriptomic, and brain imaging modalities toward common biological targets.

The human lncRNA GOMAFU suppresses neuronal interferon response pathways affected in neuropsychiatric diseases

Long noncoding RNAs (lncRNAs) play multifaceted roles in regulating brain gene networks. LncRNA abnormalities are thought to underlie the complex etiology of numerous neuropsychiatric disorders. One example is the human lncRNA gene GOMAFU, which is found dysregulated in schizophrenia (SCZ) postmortem brains and harbors genetic variants that contribute to the risk of SCZ. However, transcriptome-wide biological pathways regulated by GOMAFU have not been determined. How GOMAFU dysregulation contributes to SCZ pathogenesis remains elusive. Here we report that GOMAFU is a novel suppressor of human neuronal interferon (IFN) response pathways that are hyperactive in the postmortem SCZ brains. We analyzed recently released transcriptomic profiling datasets in clinically relevant brain areas derived from multiple SCZ cohorts and found brain region-specific dysregulation of GOMAFU. Using CRISPR-Cas9 to delete the GOMAFU promoter in a human neural progenitor cell model, we identified transcriptomic alterations caused by GOMAFU deficiency in pathways commonly affected in postmortem brains of SCZ and autism spectrum disorder (ASD), with the most striking effects on upregulation of numerous genes underlying IFN signaling. In addition, expression levels of GOMAFU target genes in the IFN pathway are differentially affected in SCZ brain regions and negatively associated with GOMAFU alterations. Furthermore, acute exposure to IFN-γ causes a rapid decline of GOMAFU and activation of a subclass of GOMAFU targets in stress and immune response pathways that are affected in SCZ brains, which form a highly interactive molecular network. Together, our studies unveiled the first evidence of lncRNA-governed neuronal response pathways to IFN challenge and suggest that GOMAFU dysregulation may mediate environmental risks and contribute to etiological neuroinflammatory responses by brain neurons of neuropsychiatric diseases.

MyBrain-Seq: A Pipeline for MiRNA-Seq Data Analysis in Neuropsychiatric Disorders

High-throughput sequencing of small RNA molecules such as microRNAs (miRNAs) has become a widely used approach for studying gene expression and regulation. However, analyzing miRNA-Seq data can be challenging because it requires multiple steps, from quality control and preprocessing to differential expression and pathway-enrichment analyses, with many tools and databases available for each step. Furthermore, reproducibility of the analysis pipeline is crucial to ensure that the results are accurate and reliable. Here, we present myBrain-Seq, a comprehensive and reproducible pipeline for analyzing miRNA-Seq data that incorporates miRNA-specific solutions at each step of the analysis. The pipeline was designed to be flexible and user-friendly, allowing researchers with different levels of expertise to perform the analysis in a standardized and reproducible manner, using the most common and widely used tools for each step. In this work, we describe the implementation of myBrain-Seq and demonstrate its capacity to consistently and reproducibly identify differentially expressed miRNAs and enriched pathways by applying it to a real case study in which we compared schizophrenia patients who responded to medication with treatment-resistant schizophrenia patients to obtain a 16-miRNA treatment-resistant schizophrenia profile.

Dysregulation of Long Intergenic Non-Coding RNA Expression in the Schizophrenia Brain

BACKGROUND

Transcriptomic studies of the brains of schizophrenia (SZ) patients have produced abundant but largely inconsistent findings about the disorders pathophysiology. These inconsistencies might stem not only from the heterogeneous nature of the disorder, but also from the unbalanced focus on particular cortical regions and protein-coding genes. Compared to protein-coding transcripts, long intergenic non-coding RNA (lincRNA) display substantially greater brain region and disease response specificity, positioning them as prospective indicators of SZ-associated alterations. Further, a growing understanding of the systemic character of the disorder calls for a more systematic screening involving multiple diverse brain regions.

AIM

We aimed to identify and interpret alterations of the lincRNA expression profiles in SZ by examining the transcriptomes of 35 brain regions.

METHODS

We measured the transcriptome of 35 brain regions dissected from eight adult brain specimens, four SZ patients, and four healthy controls, using high-throughput RNA sequencing. Analysis of these data yielded 861 annotated human lincRNAs passing the detection threshold.

RESULTS

Of the 861 detected lincRNA, 135 showed significant region-dependent expression alterations in SZ (two-way ANOVA, BH-adjusted p 0.05) and 37 additionally showed significant differential expression between HC and SZ individuals in at least one region (post hoc Tukey test, p 0.05). For these 37 differentially expressed lincRNAs (DELs), 88% of the differences occurred in a cluster of brain regions containing axon-rich brain regions and cerebellum. Functional annotation of the DEL targets further revealed stark enrichment in neurons and synaptic transmission terms and pathways.

CONCLUSION

Our study highlights the utility of a systematic brain transcriptome analysis relying on the expression profiles measured across multiple brain regions and singles out white matter regions as a prospective target for further SZ research.

Unraveling Psychiatric Disorders through Neural Single-Cell Transcriptomics Approaches

The development of single-cell and single-nucleus transcriptome technologies is enabling the unraveling of the molecular and cellular heterogeneity of psychiatric disorders. The complexity of the brain and the relationships between different brain regions can be better understood through the classification of individual cell populations based on their molecular markers and transcriptomic features. Analysis of these unique cell types can explain their involvement in the pathology of psychiatric disorders. Recent studies in both human and animal models have emphasized the importance of transcriptome analysis of neuronal cells in psychiatric disorders but also revealed critical roles for non-neuronal cells, such as oligodendrocytes and microglia. In this review, we update current findings on the brain transcriptome and explore molecular studies addressing transcriptomic alterations identified in human and animal models in depression and stress, neurodegenerative disorders (Parkinson's and Alzheimer's disease), schizophrenia, opioid use disorder, and alcohol and psychostimulant abuse. We also comment on potential future directions in single-cell and single-nucleus studies.

Current State of Modeling Human Psychiatric Disorders Using Zebrafish

Psychiatric disorders are highly prevalent brain pathologies that represent an urgent, unmet biomedical problem. Since reliable clinical diagnoses are essential for the treatment of psychiatric disorders, their animal models with robust, relevant behavioral and physiological endpoints become necessary. Zebrafish (Danio rerio) display well-defined, complex behaviors in major neurobehavioral domains which are evolutionarily conserved and strikingly parallel to those seen in rodents and humans. Although zebrafish are increasingly often used to model psychiatric disorders, there are also multiple challenges with such models as well. The field may therefore benefit from a balanced, disease-oriented discussion that considers the clinical prevalence, the pathological complexity, and societal importance of the disorders in question, and the extent of its detalization in zebrafish central nervous system (CNS) studies. Here, we critically discuss the use of zebrafish for modeling human psychiatric disorders in general, and highlight the topics for further in-depth consideration, in order to foster and (re)focus translational biological neuroscience research utilizing zebrafish. Recent developments in molecular biology research utilizing this model species have also been summarized here, collectively calling for a wider use of zebrafish in translational CNS disease modeling.

Overexpression of transmembrane TNFα in brain endothelial cells induces schizophrenia-relevant behaviors

Upregulation of genes and coexpression networks related to immune function and inflammation have been repeatedly reported in the brain of individuals with schizophrenia. However, a causal relationship between the abnormal immune/inflammation-related gene expression and schizophrenia has not been determined. We conducted co-expression networks using publicly available RNA-seq data from prefrontal cortex (PFC) and hippocampus (HP) of 64 individuals with schizophrenia and 64 unaffected controls from the SMRI tissue collections. We identified proinflammatory cytokine, transmembrane tumor necrosis factor-α (tmTNFα), as a potential regulator in the module of co-expressed genes that we find related to the immune/inflammation response in endothelial cells (ECs) and/or microglia of the brain of individuals with schizophrenia. The immune/inflammation-related modules associated with schizophrenia and the TNF signaling pathway that regulate the network were replicated in an independent cohort of brain samples from 68 individuals with schizophrenia and 135 unaffected controls. To investigate the association between the overexpression of tmTNFα in brain ECs and schizophrenia-like behaviors, we induced short-term overexpression of the uncleavable form of (uc)-tmTNFα in ECs of mouse brain for 7 weeks. We found schizophrenia-relevant behavioral deficits in these mice, including cognitive impairment, abnormal sensorimotor gating, and sensitization to methamphetamine (METH) induced locomotor activity and METH-induced neurotransmitter levels. These uc-tmTNFα effects were mediated by TNF receptor2 (TNFR2) and induced activation of TNFR2 signaling in astrocytes and neurons. A neuronal module including neurotransmitter signaling pathways was down-regulated in the brain of mice by the short-term overexpression of the gene, while an immune/inflammation-related module was up-regulated in the brain of mice after long-term expression of 22 weeks. Our results indicate that tmTNFα may play a direct role in regulating neurotransmitter signaling pathways that contribute to the clinical features of schizophrenia.

Infections, Inflammation, and Psychiatric Illness: Review of Postmortem Evidence

While there is an abundance of epidemiological evidence implicating infectious agents in the etiology of severe mental illnesses, postmortem studies have not yet detected an increased incidence of microbial nucleic acid or proteins in the brains of people with mental illness. Nevertheless, abnormally expressed immune and inflammatory markers have consistently been found in the postmortem brain of patients with schizophrenia and mood disorders. Some of these abnormalities may be the result of an infection in utero or early in life that not only impacted the developing immune system but also the developing neurons of the brain. Some of the immune markers that are consistently found to be upregulated in schizophrenia implicate a possible viral infection and the blood brain barrier in the etiology and neuropathology of the disorder.

The gene expression patterns as surrogate indices of pH in the brain

Hydrogen ion (H⁺) is one of the most potent intrinsic neuromodulators in the brain in terms of concentration. Changes in H⁺ concentration, expressed as pH, are thought to be associated with various biological processes, such as gene expression, in the brain. Accumulating evidence suggests that decreased brain pH is a common feature of several neuropsychiatric disorders, including schizophrenia, bipolar disorder, autism spectrum disorder, and Alzheimer's disease. However, it remains unclear whether gene expression patterns can be used as surrogates for pH changes in the brain. In this study, we performed meta-analyses using publicly available gene expression datasets to profile the expression patterns of pH-associated genes, whose expression levels were correlated with brain pH, in human patients and mouse models of major central nervous system (CNS) diseases, as well as in mouse cell-type datasets. Comprehensive analysis of 281 human datasets from 11 CNS disorders revealed that gene expression associated with decreased pH was over-represented in disorders including schizophrenia, bipolar disorder, autism spectrum disorders, Alzheimer's disease, Huntington's disease, Parkinson's disease, and brain tumors. Expression patterns of pH-associated genes in mouse models of neurodegenerative disease showed a common time course trend toward lower pH over time. Furthermore, cell type analysis identified astrocytes as the cell type with the most acidity-related gene expression, consistent with previous experimental measurements showing a lower intracellular pH in astrocytes than in neurons. These results suggest that the expression pattern of pH-associated genes may be a surrogate for the state- and trait-related changes in pH in brain cells. Altered expression of pH-associated genes may serve as a novel molecular mechanism for a more complete understanding of the transdiagnostic pathophysiology of neuropsychiatric and neurodegenerative disorders.

Two subtypes of schizophrenia identified by an individual-level atypical pattern of tensor-based morphometric measurement

Difficulties in parsing the multiaspect heterogeneity of schizophrenia (SCZ) based on current nosology highlight the need to subtype SCZ using objective biomarkers. Here, utilizing a large-scale multisite SCZ dataset, we identified and validated 2 neuroanatomical subtypes with individual-level abnormal patterns of the tensor-based morphometric measurement. Remarkably, compared with subtype 1, which showed moderate deficits of some subcortical nuclei and an enlarged striatum and cerebellum, subtype 2, which showed cerebellar atrophy and more severe subcortical nuclei atrophy, had a higher subscale score of negative symptoms, which is considered to be a core aspect of SCZ and is associated with functional outcome. Moreover, with the neuroimaging-clinic association analysis, we explored the detailed relationship between the heterogeneity of clinical symptoms and the heterogeneous abnormal neuroanatomical patterns with respect to the 2 subtypes. And the neuroimaging-transcription association analysis highlighted several potential heterogeneous biological factors that may underlie the subtypes. Our work provided an effective framework for investigating the heterogeneity of SCZ from multilevel aspects and may provide new insights for precision psychiatry.

Unfolding the cascade of SERPINA3: Inflammation to cancer

SERine Protease INhibitor clade A member 3 (SERPINA3), a member of the SERine-Protease INhibitor (SERPIN) superfamily, principally works as a protease inhibitor in maintaining cellular homeostasis. It is a matricellular acute-phase glycoprotein that appears to be the sole nuclear-binding secretory serpin. Several studies have emerged in recent years demonstrating its link to cancer and disease biology. SERPINA3 seems to have cancer- and compartment-specific biological functions, acting either as a tumour promoter or suppressor in different cancers. However, the localization, mechanism of action and the effectors of SERPINA3 in physiological and pathological scenarios remain obscure. Our review aims to consolidate the current evidence of SERPINA3 in various cancers, highlighting its association with the cancer hallmarks and ratifying its status as an emerging cancer biomarker. The elucidation of SERPINA3-mediated cancer progression and its targeting might shed light on the realm of cancer therapeutics.

Schizophrenia is defined by cell-specific neuropathology and multiple neurodevelopmental mechanisms in patient-derived cerebral organoids

Due to an inability to ethically access developing human brain tissue as well as identify prospective cases, early-arising neurodevelopmental and cell-specific signatures of Schizophrenia (Scz) have remained unknown and thus undefined. To overcome these challenges, we utilized patient-derived induced pluripotent stem cells (iPSCs) to generate 3D cerebral organoids to model neuropathology of Scz during this critical period. We discovered that Scz organoids exhibited ventricular neuropathology resulting in altered progenitor survival and disrupted neurogenesis. This ultimately yielded fewer neurons within developing cortical fields of Scz organoids. Single-cell sequencing revealed that Scz progenitors were specifically depleted of neuronal programming factors leading to a remodeling of cell-lineages, altered differentiation trajectories, and distorted cortical cell-type diversity. While Scz organoids were similar in their macromolecular diversity to organoids generated from healthy controls (Ctrls), four GWAS factors (PTN, COMT, PLCL1, and PODXL) and peptide fragments belonging to the POU-domain transcription factor family (e.g., POU3F2/BRN2) were altered. This revealed that Scz organoids principally differed not in their proteomic diversity, but specifically in their total quantity of disease and neurodevelopmental factors at the molecular level. Single-cell sequencing subsequently identified cell-type specific alterations in neuronal programming factors as well as a developmental switch in neurotrophic growth factor expression, indicating that Scz neuropathology can be encoded on a cell-type-by-cell-type basis. Furthermore, single-cell sequencing also specifically replicated the depletion of BRN2 (POU3F2) and PTN in both Scz progenitors and neurons. Subsequently, in two mechanistic rescue experiments we identified that the transcription factor BRN2 and growth factor PTN operate as mechanistic substrates of neurogenesis and cellular survival, respectively, in Scz organoids. Collectively, our work suggests that multiple mechanisms of Scz exist in patient-derived organoids, and that these disparate mechanisms converge upon primordial brain developmental pathways such as neuronal differentiation, survival, and growth factor support, which may amalgamate to elevate intrinsic risk of Scz.

Brain transcriptomic profiling reveals common alterations across neurodegenerative and psychiatric disorders

Neurodegenerative and neuropsychiatric disorders (ND-NPs) are multifactorial, polygenic and complex behavioral phenotypes caused by brain abnormalities. Large-scale collaborative efforts have tried to identify the genetic architecture of these conditions. However, the specific and shared underlying molecular pathobiology of brain illnesses is not clear. Here, we examine transcriptome-wide characterization of eight conditions, using a total of 2,633 post-mortem brain samples from patients with Alzheimer’s disease (AD), Parkinson’s disease (PD), Progressive Supranuclear Palsy (PSP), Pathological Aging (PA), Autism Spectrum Disorder (ASD), Schizophrenia (Scz), Major Depressive Disorder (MDD), and Bipolar Disorder (BP)–in comparison with 2,078 brain samples from matched control subjects.

Similar transcriptome alterations were observed between NDs and NPs with the top correlations obtained between Scz-BP, ASD-PD, AD-PD, and Scz-ASD. Region-specific comparisons also revealed shared transcriptome alterations in frontal and temporal lobes across NPs and NDs. Co-expression network analysis identified coordinated dysregulations of cell-type-specific modules across NDs and NPs. This study provides a transcriptomic framework to understand the molecular alterations of NPs and NDs through their shared- and specific gene expression in the brain.

16p11 Duplication Disrupts Hippocampal-Orbitofrontal-Amygdala Connectivity, Revealing a Neural Circuit Endophenotype for Schizophrenia

Chromosome 16p11.2 duplications dramatically increase risk for schizophrenia, but the mechanisms remain largely unknown. Here, we show that mice with an equivalent genetic mutation (16p11.2 duplication mice) exhibit impaired hippocampal-orbitofrontal and hippocampal-amygdala functional connectivity. Expression of schizophrenia-relevant GABAergic cell markers (parvalbumin and calbindin) is selectively decreased in orbitofrontal cortex, while somatostatin expression is decreased in lateral amygdala. When 16p11.2 duplication mice are tested in cognitive tasks dependent on hippocampal-orbitofrontal connectivity, performance is impaired in an 8-arm maze "N-back" working memory task and in a touchscreen continuous performance task. Consistent with hippocampal-amygdala dysconnectivity, deficits in ethologically relevant social behaviors are also observed. Overall, the cellular/molecular, brain network, and behavioral alterations markedly mirror those observed in schizophrenia patients. Moreover, the data suggest that 16p11.2 duplications selectively impact hippocampal-amygdaloid-orbitofrontal circuitry, supporting emerging ideas that dysfunction in this network is a core element of schizophrenia and defining a neural circuit endophenotype for the disease.

Thalamic, Amygdalar, and hippocampal nuclei morphology and their trajectories in first episode psychosis: A preliminary longitudinal study✰

The thalamus, amygdala, and hippocampus play important pathophysiologic roles in psychosis. Few studies have prospectively examined subcortical nuclei in relation to predicting clinical outcomes after a first-episode of psychosis (FEP). Here, we examined volumetric differences and trajectories among subcortical nuclei in FEP patients and their associations with illness severity. Clinical and brain volume measures were collected using a 1.5T MRI scanner and processed using FreeSurfer 6.0 from a prospective study of antipsychotic-naïve FEP patients of FEP-schizophrenia (FEP-SZ) (baseline, n = 38; follow-up, n = 17), FEP non-schizophrenia (FEP-NSZ) (baseline, n = 23; follow-up, n = 13), and healthy controls (HCs) (baseline, n = 47; follow-up, n = 29). Compared to FEP-NSZ and HCs, FEP-SZ had significantly smaller thalamic anterior nuclei volume at baseline. Longitudinally, FEP-SZ showed a positive rate of change in the amygdala compared to controls or FEP-NSZ, as well as in the basal, central and accessory basal nuclei compared to FEP-NSZ. Enlargement in the thalamic anterior nuclei predicted a worsening in overall psychosis symptoms. Baseline thalamic anterior nuclei alterations further specify key subcortical regions associated with FEP-SZ pathophysiology. Longitudinally, anterior nuclei volume enlargement may signal symptomatic worsening. The amygdala and thalamus structures may show diagnostic differences between schizophrenia and non-schizophrenia psychoses, while the thalamus changes may reflect disease or treatment related changes in clinical outcome.

Brain-immune interactions in neuropsychiatric disorders: Lessons from transcriptome studies for molecular targeting

Understanding the pathophysiological mechanisms of neuropsychiatric disorders has been a challenging quest for neurobiologists. Recent years have witnessed enormous technological advances in the field of neuroimmunology, blurring boundaries between the central nervous system and the periphery. Consequently, the discipline has expanded to cover interactions between the nervous and immune systems in health and diseases. The complex interplay between the peripheral and central immune pathways in neuropsychiatric disorders has recently been documented in various studies, but the genetic determinants remain elusive. Recent transcriptome studies have identified dysregulated genes involved in peripheral immune cell activation, blood-brain barrier integrity, glial cell activation, and synaptic plasticity in major depressive disorder, bipolar disorder, autism spectrum disorder, and schizophrenia. Herein, the key transcriptomic techniques applied in investigating differentially expressed genes and pathways responsible for altered brain-immune interactions in neuropsychiatric disorders are discussed. The application of transcriptomics that can aid in identifying molecular targets in various neuropsychiatric disorders is highlighted.

Concordance of Immune-Related Markers in Lymphocytes and Prefrontal Cortex in Schizophrenia

Schizophrenia is a severe neuropsychiatric disorder associated with a wide array of transcriptomic and neurobiochemical changes. Genome-wide transcriptomic profiling conducted in postmortem brain have provided novel insights into the pathophysiology of this disorder, and identified biological processes including immune/inflammatory-related responses, metabolic, endocrine, and synaptic function. However, few studies have investigated whether similar changes are present in peripheral tissue. Here, we used RNA-sequencing to characterize transcriptomic profiles of lymphocytes in 18 nonpsychotic controls and 19 individuals with schizophrenia. We identified 2819 differentially expressed transcripts (P nominal < .05) in the schizophrenia group when compared to controls. Bioinformatic analyses conducted on a subset of 293 genes (P nominal < .01 and |log2 FC| > 0.5) highlighted immune/inflammatory responses as key biological processes in our dataset. Differentially expressed genes in lymphocytes were highly enriched in gene expression profiles associated with cortex layer 5a and immune cells. Thus, we investigated whether the changes in transcripts levels observed in lymphocytes could also be detected in the prefrontal cortex (PFC, BA10) in a second replication cohort of schizophrenia subjects. Remarkably, mRNA levels detected in the PFC and lymphocytes were in strong agreement, and measurements obtained using RNA-sequencing positively correlated with data obtained by reverse transcriptase-quantitative polymerase chain reaction analysis. Collectively, our work supports a role for immune dysfunction in the pathogenesis of schizophrenia and suggests that peripheral markers can be used as accessible surrogates to investigate putative central nervous system disruptions.

nMAGMA: a network-enhanced method for inferring risk genes from GWAS summary statistics and its application to schizophrenia

MOTIVATION

Annotating genetic variants from summary statistics of genome-wide association studies (GWAS) is crucial for predicting risk genes of various disorders. The multimarker analysis of genomic annotation (MAGMA) is one of the most popular tools for this purpose, where MAGMA aggregates signals of single nucleotide polymorphisms (SNPs) to their nearby genes. In biology, SNPs may also affect genes that are far away in the genome, thus missed by MAGMA. Although different upgrades of MAGMA have been proposed to extend gene-wise variant annotations with more information (e.g. Hi-C or eQTL), the regulatory relationships among genes and the tissue specificity of signals have not been taken into account.

RESULTS

We propose a new approach, namely network-enhanced MAGMA (nMAGMA), for gene-wise annotation of variants from GWAS summary statistics. Compared with MAGMA and H-MAGMA, nMAGMA significantly extends the lists of genes that can be annotated to SNPs by integrating local signals, long-range regulation signals (i.e. interactions between distal DNA elements), and tissue-specific gene networks. When applied to schizophrenia (SCZ), nMAGMA is able to detect more risk genes (217% more than MAGMA and 57% more than H-MAGMA) that are involved in SCZ compared with MAGMA and H-MAGMA, and more of nMAGMA results can be validated with known SCZ risk genes. Some disease-related functions (e.g. the ATPase pathway in Cortex) are also uncovered in nMAGMA but not in MAGMA or H-MAGMA. Moreover, nMAGMA provides tissue-specific risk signals, which are useful for understanding disorders with multitissue origins.

Increased RNA editing in maternal immune activation model of neurodevelopmental disease

The etiology of major neurodevelopmental disorders such as schizophrenia and autism is unclear, with evidence supporting a combination of genetic factors and environmental insults, including viral infection during pregnancy. Here we utilized a mouse model of maternal immune activation (MIA) with the viral mimic PolyI:C infection during early gestation. We investigated the transcriptional changes in the brains of mouse fetuses following MIA during the prenatal period, and evaluated the behavioral and biochemical changes in the adult brain. The results reveal an increase in RNA editing levels and dysregulation in brain development-related gene pathways in the fetal brains of MIA mice. These MIA-induced brain editing changes are not observed in adulthood, although MIA-induced behavioral deficits are observed. Taken together, our findings suggest that MIA induces transient dysregulation of RNA editing at a critical time in brain development.

Reelin Supplementation Into the Hippocampus Rescues Abnormal Behavior in a Mouse Model of Neurodevelopmental Disorders

In the majority of schizophrenia patients, chronic atypical antipsychotic administration produces a significant reduction in or even complete remission of psychotic symptoms such as hallucinations and delusions. However, these drugs are not effective in improving cognitive and emotional deficits in patients with schizophrenia. Atypical antipsychotic drugs have a high affinity for the dopamine D₂ receptor, and a modest affinity for the serotonin 5-HT_2A receptor. The cognitive and emotional deficits in schizophrenia are thought to involve neural networks beyond the classical dopaminergic mesolimbic pathway, however, including serotonergic systems. For example, mutations in the RELN gene, which encodes Reelin, an extracellular matrix protein involved in neural development and synaptic plasticity, are associated with neurodevelopmental disorders such as schizophrenia and autism spectrum disorder. Furthermore, hippocampal Reelin levels are down-regulated in the brains of both schizophrenic patients and in rodent models of schizophrenia. In the present study, we investigated the effect of Reelin microinjection into the mouse hippocampus on behavioral phenotypes to evaluate the role of Reelin in neurodevelopmental disorders and to test a therapeutic approach that extends beyond classical monoamine targets. To model the cognitive and emotional deficits, as well as histological decreases in Reelin-positive cell numbers and hippocampal synaptoporin distribution, a synaptic vesicle protein, offspring that were prenatally exposed to maternal immune activation were used. Microinjections of recombinant Reelin protein into the hippocampus rescued impairments in object memory and anxiety-like behavior and recruited synaptoporin in the hippocampus in offspring exposed to antenatal inflammation. These results suggest that Reelin supplementation has the potential to treat cognitive and emotional impairments, as well as synaptic disturbances, in patients with neurodevelopmental disorders such as schizophrenia.

Gene expression and response prediction to amisulpride in the OPTiMiSE first episode psychoses

A fundamental shortcoming in the current treatment of schizophrenia is the lack of valid criteria to predict who will respond to antipsychotic treatment. The identification of blood-based biological markers of the therapeutic response would enable clinicians to identify the subgroup of patients in whom conventional antipsychotic treatment is ineffective and offer alternative treatments. As part of the Optimisation of Treatment and Management of Schizophrenia in Europe (OPTiMiSE) programme, we conducted an RNA-Seq analysis on 188 subjects with first episode psychosis, all of whom were subsequently treated with amisulpride for 4 weeks. We compared gene expression on total RNA from patients' blood before and after treatment and identified 32 genes for which the expression changed after treatment in good responders only. These findings were replicated in an independent sample of 24 patients with first episode psychosis. Six genes showed a significant difference in expression level between good and poor responders before starting treatment, allowing to predict treatment outcome with a predictive value of 93.8% when combined with clinical features. Collectively, these findings identified new mechanisms to explain symptom improvement after amisulpride medication and highlight the potential of combining gene expression profiling with clinical data to predict treatment response in first episode psychoses.

Lasting and Sex-Dependent Impact of Maternal Immune Activation on Molecular Pathways of the Amygdala

The prolonged and sex-dependent impact of maternal immune activation (MIA) during gestation on the molecular pathways of the amygdala, a brain region that influences social, emotional, and other behaviors, is only partially understood. To address this gap, we investigated the effects of viral-elicited MIA during gestation on the amygdala transcriptome of pigs, a species of high molecular and developmental homology to humans. Gene expression levels were measured using RNA-Seq on the amygdala for 3-week-old female and male offspring from MIA and control groups. Among the 403 genes that exhibited significant MIA effect, a prevalence of differentially expressed genes annotated to the neuroactive ligand-receptor pathway, glutamatergic functions, neuropeptide systems, and cilium morphogenesis were uncovered. Genes in these categories included corticotropin-releasing hormone receptor 2, glutamate metabotropic receptor 4, glycoprotein hormones, alpha polypeptide, parathyroid hormone 1 receptor, vasointestinal peptide receptor 2, neurotensin, proenkephalin, and gastrin-releasing peptide. These categories and genes have been associated with the MIA-related human neurodevelopmental disorders, including schizophrenia and autism spectrum disorders. Gene network reconstruction highlighted differential vulnerability to MIA effects between sexes. Our results advance the understanding necessary for the development of multifactorial therapies targeting immune modulation and neurochemical dysfunction that can ameliorate the effects of MIA on offspring behavior later in life.

Regulation of Janus Kinase 2 by an Inflammatory Bowel Disease Causal Non-coding Single Nucleotide Polymorphism

BACKGROUND AND AIMS

Among the >240 genetic loci described to date which confer susceptibility to inflammatory bowel disease, a small subset have been fine-mapped to an individual, non-coding single nucleotide polymorphism [SNP]. To illustrate a model mechanism by which a presumed-causal non-coding SNP can function, we analysed rs1887428, located in the promoter region of the Janus kinase 2 [JAK2] gene.

METHODS

We utilized comparative affinity purification-mass spectrometry, DNA-protein binding assays, CRISPR/Cas9 genome editing, transcriptome sequencing and methylome quantitative trait locus methods to characterize the role of this SNP.

RESULTS

We determined that the risk allele of rs1887428 is bound by the transcription factor [TF] RBPJ, while the protective allele is bound by the homeobox TF CUX1. While rs188748 only has a very modest influence on JAK2 expression, this effect was amplified downstream through the expression of pathway member STAT5B and epigenetic modification of the JAK2 locus.

CONCLUSION

Despite the absence of a consensus TF-binding motif or expression quantitative trait locus, we have used improved methods to characterize a putatively causal SNP to yield insight into inflammatory bowel disease mechanisms.

PODCAST

This article has an associated podcast which can be accessed at https://academic.oup.com/ecco-jcc/pages/podcast.

Molecular alterations in the medial temporal lobe in schizophrenia

The medial temporal lobe (MTL) and its individual structures have been extensively implicated in schizophrenia pathophysiology, with considerable efforts aimed at identifying structural and functional differences in this brain region. The major structures of the MTL for which prominent differences have been revealed include the hippocampus, the amygdala and the superior temporal gyrus (STG). The different functions of each of these regions have been comprehensively characterized, and likely contribute differently to schizophrenia. While neuroimaging studies provide an essential framework for understanding the role of these MTL structures in various aspects of the disease, ongoing efforts have sought to employ molecular measurements in order to elucidate the biology underlying these macroscopic differences. This review provides a summary of the molecular findings in three major MTL structures, and discusses convergent findings in cellular architecture and inter-and intra-cellular networks. The findings of this effort have uncovered cell-type, network and gene-level specificity largely unique to each brain region, indicating distinct molecular origins of disease etiology. Future studies should test the functional implications of these molecular changes at the circuit level, and leverage new advances in sequencing technology to further refine our understanding of the differential contribution of MTL structures to schizophrenia.

Involvement of hemoglobins in the pathophysiology of Alzheimer's disease

Hemoglobins (Hbs) are heme-containing proteins binding oxygen, carbon monoxide, and nitric oxide. While erythrocytes are the most well-known location of Hbs, Hbs also exist in neurons, glia and oligodendroglia and they are primarily localized in the inner mitochondrial membrane of neurons with likely roles in cellular respiration and buffering protons. Recently, studies have suggested links between hypoxia and neurodegenerative disorders such as Alzheimer Disease (AD) and furthermore suggested involvement of Hbs in the pathogenesis of AD. While cellular immunohistochemical studies on AD brains have observed reduced levels of Hb in the cytoplasm of pre-tangle and tangle-bearing neurons, other studies on homogenates of AD brain samples observed increased Hb levels. This potential discrepancy may result from differential presence and function of intracellular versus extracellular Hbs. Intracellular Hbs may protect neurons against hypoxia and hyperoxia. On the other hand, extracellular free Hb and its degradation products may trigger inflammatory immune and oxidative reactions against neural macromolecules and/or damage the blood-brain barrier. Therefore, biological processes leading to reduction of Hb transcription (including clinically silent Hb mutations) may influence intra-erythrocytic and neural Hbs, and reduce the transport of oxygen, carbon monoxide and nitric oxide which may be involved in the (patho)physiology of neurodegenerative disorders such as AD. Agents such as erythropoietin, which stimulate both erythropoiesis, reduce eryptosis and induce intracellular neural Hbs may exert multiple beneficial effects on the onset and course of AD. Thus, evidence accumulates for a role of Hbs in the central nervous system while Hbs deserve more attention as possible candidate molecules involved in AD.

Analyzing DNA methylation patterns in subjects diagnosed with schizophrenia using machine learning methods

Schizophrenia is a common mental disorder with high heritability. It is genetically complex and to date more than a hundred risk loci have been identified. Association of environmental factors and schizophrenia has also been reported, while epigenetic analyses have yielded ambiguous and sometimes conflicting results. Here, we analyzed fresh frozen post-mortem brain tissue from a cohort of 73 subjects diagnosed with schizophrenia and 52 control samples, using the Illumina Infinium HumanMethylation450 Bead Chip, to investigate genome-wide DNA methylation patterns in the two groups. Analysis of differential methylation was performed with the Bioconductor Minfi package and modern machine-learning and visualization techniques, which were shown previously to be successful in detecting and highlighting differentially methylated patterns in case-control studies. In this dataset, however, these methods did not uncover any significant signals discerning the patient group and healthy controls, suggesting that if there are methylation changes associated with schizophrenia, they are heterogeneous and complex with small effect.

Deletion of Dtnbp1 in mice impairs threat memory consolidation and is associated with enhanced inhibitory drive in the amygdala

Schizophrenia is a severe and highly heritable disorder. Dystrobrevin-binding protein 1 (DTNBP1), also known as dysbindin-1, has been implicated in the pathophysiology of schizophrenia. Specifically, dysbindin-1 mRNA and protein expression are decreased in the brains of subjects with this disorder. Mice lacking dysbinidn-1 also display behavioral phenotypes similar to those observed in schizophrenic patients. However, it remains unknown whether deletion of dysbindin-1 impacts functions of the amygdala, a brain region that is critical for emotional processing, which is disrupted in patients with schizophrenia. Here, we show that dysbindin-1 is expressed in both excitatory and inhibitory neurons of the basolateral amygdala (BLA). Deletion of dysbindin-1 in male mice (Dys-/-) impaired cued and context-dependent threat memory, without changes in measures of anxiety. The behavioral deficits observed in Dys-/- mice were associated with perturbations in the BLA, including the enhancement of GABAergic inhibition of pyramidal neurons, increased numbers of parvalbumin interneurons, and morphological abnormalities of dendritic spines on pyramidal neurons. Our findings highlight an important role for dysbindin-1 in the regulation of amygdalar function and indicate that enhanced inhibition of BLA pyramidal neuron activity may contribute to the weakened threat memory expression observed in Dys-/- mice.

Identifying Hub Genes for Heat Tolerance in Water Buffalo (Bubalus bubalis) Using Transcriptome Data

Heat stress has a detrimental effect on the physiological and production performance of buffaloes. Elucidating the underlying mechanisms of heat stress is challenging, therefore identifying candidate genes is urgent and necessary. We evaluated the response of buffaloes (n = 30) to heat stress using the physiological parameters, ELISA indexes, and hematological parameters. We then performed mRNA and microRNA (miRNA) expression profiles analysis between heat tolerant (HT, n = 4) and non-heat tolerant (NHT, n = 4) buffaloes, as well as the specific modules, significant genes, and miRNAs related to the heat tolerance identified using the weighted gene co-expression network analysis (WGCNA). The results indicated that the buffaloes in HT had a significantly lower rectal temperature (RT) and respiratory rate (RR) and displayed a higher plasma heat shock protein (HSP70 and HSP90) and cortisol (COR) levels than those of NHT buffaloes. Differentially expressed analysis revealed a total of 753 differentially expressed genes (DEGs) and 16 differentially expressed miRNAs (DEmiRNAs) were identified between HT and NHT. Using the WGCNA analysis, these DEGs assigned into 5 modules, 4 of which were significantly correlation with the heat stress indexes. Interestingly, 158 DEGs associated with heat tolerance in the turquoise module were identified, 35 of which were found within the protein-protein interaction network. Several hub genes (IL18RAP, IL6R, CCR1, PPBP, IL1B, and IL1R1) were identified that significantly enriched in the Cytokine-cytokine receptor interaction. The findings may help further elucidate the underlying mechanisms of heat tolerance in buffaloes.

Shared co-expression networks in frontal cortex of the normal aged brain and schizophrenia

Previous studies on the brain of people with schizophrenia have identified structural changes and gene expression changes, suggesting that brain aging maybe accelerated in people with schizophrenia. To better characterize gene expression profiles in schizophrenia and in the aged population we constructed co-expression networks using RNA-Seq data from frontal cortex. The first data set analysed was from 62 subjects with schizophrenia and 51 unaffected controls ranging in age from 19 to 63 years. The second separate data set was from normal control individuals ranging in age from 29 to 106 years. In the first data set, we found two co-expression modules significantly associated with schizophrenia. One was a downregulated co-expression module enriched for neuron function related genes and the other was an upregulated immune/inflammation-related module. In the second data set of normal individuals, we found seven co-expression modules significantly correlated with age. A comparison of the co-expression modules from the two data sets revealed a significant consensus in nodes associated with schizophrenia and those associated with normal aging. The results indicate that a co-expression module related to neuronal function is downregulated and an immune/inflammation related co-expression module is upregulated, and associated with cells of the blood vessels, in both schizophrenia and in normal aging. This finding adds further support to the hypothesis that there may be accelerated brain aging in schizophrenia.

Towards precision medicine for pain: diagnostic biomarkers and repurposed drugs

We endeavored to identify objective blood biomarkers for pain, a subjective sensation with a biological basis, using a stepwise discovery, prioritization, validation, and testing in independent cohorts design. We studied psychiatric patients, a high risk group for co-morbid pain disorders and increased perception of pain. For discovery, we used a powerful within-subject longitudinal design. We were successful in identifying blood gene expression biomarkers that were predictive of pain state, and of future emergency department (ED) visits for pain, more so when personalized by gender and diagnosis. MFAP3, which had no prior evidence in the literature for involvement in pain, had the most robust empirical evidence from our discovery and validation steps, and was a strong predictor for pain in the independent cohorts, particularly in females and males with PTSD. Other biomarkers with best overall convergent functional evidence for involvement in pain were GNG7, CNTN1, LY9, CCDC144B, and GBP1. Some of the individual biomarkers identified are targets of existing drugs. Moreover, the biomarker gene expression signatures were used for bioinformatic drug repurposing analyses, yielding leads for possible new drug candidates such as SC-560 (an NSAID), and amoxapine (an antidepressant), as well as natural compounds such as pyridoxine (vitamin B6), cyanocobalamin (vitamin B12), and apigenin (a plant flavonoid). Our work may help mitigate the diagnostic and treatment dilemmas that have contributed to the current opioid epidemic.

INFERNO: inferring the molecular mechanisms of noncoding genetic variants

The majority of variants identified by genome-wide association studies (GWAS) reside in the noncoding genome, affecting regulatory elements including transcriptional enhancers. However, characterizing their effects requires the integration of GWAS results with context-specific regulatory activity and linkage disequilibrium annotations to identify causal variants underlying noncoding association signals and the regulatory elements, tissue contexts, and target genes they affect. We propose INFERNO, a novel method which integrates hundreds of functional genomics datasets spanning enhancer activity, transcription factor binding sites, and expression quantitative trait loci with GWAS summary statistics. INFERNO includes novel statistical methods to quantify empirical enrichments of tissue-specific enhancer overlap and to identify co-regulatory networks of dysregulated long noncoding RNAs (lncRNAs). We applied INFERNO to two large GWAS studies. For schizophrenia (36,989 cases, 113,075 controls), INFERNO identified putatively causal variants affecting brain enhancers for known schizophrenia-related genes. For inflammatory bowel disease (IBD) (12,882 cases, 21,770 controls), INFERNO found enrichments of immune and digestive enhancers and lncRNAs involved in regulation of the adaptive immune response. In summary, INFERNO comprehensively infers the molecular mechanisms of causal noncoding variants, providing a sensitive hypothesis generation method for post-GWAS analysis. The software is available as an open source pipeline and a web server.

Transcriptomic analysis reveals oxidative phosphorylation activation in an adolescent social isolation rat model

Complex interactions between genetic and environmental factors exert a sustained influence on the pathogenesis of schizophrenia (SCZ). Adolescent social isolation is regarded as a typical paradigm for SCZ. However, the underlying pathological mechanisms are not fully understood. In this study, adolescent Sprague-Dawley (SD) rats were placed in isolation rearing (IR) or social rearing (SR) conditions from postnatal day (PND) 21 to 34 to establish a SCZ disease model and a control model, respectively. Prepulse inhibition (PPI) assays and elevated plus maze tests were performed on PND 56. Next, prefrontal cortex (PFC) tissues were isolated for transcriptomic sequencing and RT-qPCR analyses. The results indicated that adolescent social isolation induced anxious behaviors and disrupted PPIs as well as specific PFC gene expression patterns in adult SD rats. A total of 196 genes were identified as upregulated, and 748 genes were identified as down-regulated in the IR group compared with those in the SR group. Differentially expressed genes (DEGs) were highly enriched in the KEGG pathways associated with the comorbidity of neurological disorder and oxidative phosphorylation (OXPHOS); 26 out of 27 comorbid neurological disorder-associated DEGs overlapped with 31 OXPHOS-associated DEGs. Those 26 overlapping DEGs were all upregulated in the IR group and could easily distinguish the IR group from the SR group; 6 of these DEGs (COX7C, NDUFB11, NDUFA2, NDUFC2, ATP5C1, and COX6A1) were verified by RT-qPCR. Here, we provide a systematic overview of gene expression alterations in adolescent-social-isolation-induced SCZ (ASI-SCZ), which suggests that genes that are associated with the comorbidity of neurological disorders, especially OXPHOS-related genes, contribute to the pathogenesis of ASI-SCZ.

The Long Noncoding RNA Landscape in Amygdala Tissues from Schizophrenia Patients

To date, most transcriptome studies of schizophrenia focus on the analysis of protein-coding genes. Long noncoding RNAs (lncRNAs) are emerging as key tissue-specific regulators of cellular and disease processes. The amygdala brain region has been implicated in the pathophysiology of schizophrenia. We performed unbiased whole transcriptome profiling of amygdala tissues from 22 schizophrenia patients and 24 non-psychiatric controls using RNA-seq. We reconstructed amygdala transcriptome and employed systems biology approaches to annotating the functional roles of lncRNAs. As a result, we identified 839 novel lncRNAs in amygdala. We found in amygdala lncRNAs are more subtype-specific than protein-coding genes. We identified functional modules associated with "synaptic transmission", "ribosome", and "immune responses" which were related to schizophrenia pathophysiology that involved lncRNAs. Integrative functional analyses associating individual lncRNAs with specific pathways and functions further show that amygdala lncRNAs are connected with all of these pathways. Our study presents the first systematic landscape of lncRNAs in amygdala tissue from schizophrenia cases.

SAP97 regulates behavior and expression of schizophrenia risk enriched gene sets in mouse hippocampus

Synapse associated protein of 97KDa (SAP97) belongs to a family of scaffolding proteins, the membrane-associated guanylate kinases (MAGUKs), that are highly enriched in the postsynaptic density of synapses and play an important role in organizing protein complexes necessary for synaptic development and plasticity. The Dlg-MAGUK family of proteins are structurally very similar, and an effort has been made to parse apart the unique function of each Dlg-MAGUK protein by characterization of knockout mice. Knockout mice have been generated and characterized for PSD-95, PSD-93, and SAP102, however SAP97 knockout mice have been impossible to study because the SAP97 null mice die soon after birth due to a craniofacial defect. We studied the transcriptomic and behavioral consequences of a brain-specific conditional knockout of SAP97 (SAP97-cKO). RNA sequencing from hippocampi from control and SAP97-cKO male animals identified 67 SAP97 regulated transcripts. As large-scale genetic studies have implicated MAGUKs in neuropsychiatric disorders such as intellectual disability, autism spectrum disorders, and schizophrenia (SCZ), we analyzed our differentially expressed gene (DEG) set for enrichment of disease risk-associated genes, and found our DEG set to be specifically enriched for SCZ-related genes. Subjecting SAP97-cKO mice to a battery of behavioral tests revealed a subtle male-specific cognitive deficit and female-specific motor deficit, while other behaviors were largely unaffected. These data suggest that loss of SAP97 may have a modest contribution to organismal behavior. The SAP97-cKO mouse serves as a stepping stone for understanding the unique role of SAP97 in biology.

Coexpression network analysis identifies transcriptional modules associated with genomic alterations in neuroblastoma

Neuroblastoma is a highly complex and heterogeneous cancer in children. Acquired genomic alterations including MYCN amplification, 1p deletion and 11q deletion are important risk factors and biomarkers in neuroblastoma. Here, we performed a co-expression-based gene network analysis to study the intrinsic association between specific genomic changes and transcriptome organization. We identified multiple gene coexpression modules which are recurrent in two independent datasets and associated with functional pathways including nervous system development, cell cycle, immune system process and extracellular matrix/space. Our results also indicated that modules involved in nervous system development and cell cycle are highly associated with MYCN amplification and 1p deletion, while modules responding to immune system process are associated with MYCN amplification only. In summary, this integrated analysis provides novel insights into molecular heterogeneity and pathogenesis of neuroblastoma. This article is part of a Special Issue entitled: Accelerating Precision Medicine through Genetic and Genomic Big Data Analysis edited by Yudong Cai & Tao Huang.