Search

Large-scale genome-wide association studies of schizophrenia have uncovered hundreds of associated loci but with extremely limited representation of African diaspora populations. We surveyed electronic health records of 200,000 individuals of African ancestry in the Million Veteran and All of Us Research Programs, and, coupled with genotype-level data from four case-control studies, realized a combined sample size of 13,012 affected and 54,266 unaffected persons. Three genome-wide significant signals - near PLXNA4 , PMAIP1 , and TRPA1 - are the first to be independently identified in populations of predominantly African ancestry. Joint analyses of African, European, and East Asian ancestries across 86,981 cases and 303,771 controls, yielded 376 distinct autosomal loci, which were refined to 708 putatively causal variants via multi-ancestry fine-mapping. Utilizing single-cell functional genomic data from human brain tissue and two complementary approaches, transcriptome-wide association studies and enhancer-promoter contact mapping, we identified a consensus set of 94 genes across ancestries and pinpointed the specific cell types in which they act. We identified reproducible associations of schizophrenia polygenic risk scores with schizophrenia diagnoses and a range of other mental and physical health problems. Our study addresses a longstanding gap in the generalizability of research findings for schizophrenia across ancestral populations, underlining shared biological underpinnings of schizophrenia across global populations in the presence of broadly divergent risk allele frequencies.

We elucidated the molecular fingerprint of vulnerable excitatory neurons within select cortical lamina of individuals with Down syndrome (DS) for mechanistic understanding and therapeutic potential that also informs Alzheimer's disease (AD) pathophysiology. Frontal cortex (BA9) layer III (L3) and layer V (L5) pyramidal neurons were microisolated from postmortem human DS and age- and sex-matched controls (CTR) to interrogate differentially expressed genes (DEGs) and key biological pathways relevant to neurodegenerative programs. We identified > 2300 DEGs exhibiting convergent dysregulation of gene expression in both L3 and L5 pyramidal neurons in individuals with DS versus CTR subjects. DEGs included over 100 triplicated human chromosome 21 genes in L3 and L5 neurons, demonstrating a trisomic neuronal karyotype in both laminae. In addition, thousands of other DEGs were identified, indicating gene dysregulation is not limited to trisomic genes in the aged DS brain, which we postulate is relevant to AD pathobiology. Convergent L3 and L5 DEGs highlighted pertinent biological pathways and identified key pathway-associated targets likely underlying corticocortical neurodegeneration and related cognitive decline in individuals with DS. Select key DEGs were interrogated as potential hub genes driving dysregulation, namely the triplicated DEGs amyloid precursor protein (APP) and superoxide dismutase 1 (SOD1), along with key signaling DEGs including mitogen activated protein kinase 1 and 3 (MAPK1, MAPK3) and calcium calmodulin dependent protein kinase II alpha (CAMK2A), among others. Hub DEGs determined from multiple pathway analyses identified potential therapeutic candidates for amelioration of cortical neuron dysfunction and cognitive decline in DS with translational relevance to AD., (© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Nucleotide variants in cell type-specific gene regulatory elements in the human brain are risk factors for human disease. We measured chromatin accessibility in 1932 aliquots of sorted neurons and non-neurons from 616 human postmortem brains and identified 34,539 open chromatin regions with chromatin accessibility quantitative trait loci (caQTLs). Only 10.4% of caQTLs are shared between neurons and non-neurons, which supports cell type-specific genetic regulation of the brain regulome. Incorporating allele-specific chromatin accessibility improves statistical fine-mapping and refines molecular mechanisms that underlie disease risk. Using massively parallel reporter assays in induced excitatory neurons, we screened 19,893 brain QTLs and identified the functional impact of 476 regulatory variants. Combined, this comprehensive resource captures variation in the human brain regulome and provides insights into disease etiology.

Sample-wise deconvolution methods estimate cell-type proportions and gene expressions in bulk tissue samples, yet their performance and biological applications remain unexplored, particularly in human brain transcriptomic data. Here, nine deconvolution methods were evaluated with sample-matched data from bulk tissue RNA sequencing (RNA-seq), single-cell/nuclei (sc/sn) RNA-seq, and immunohistochemistry. A total of 1,130,767 nuclei per cells from 149 adult postmortem brains and 72 organoid samples were used. The results showed the best performance of dtangle for estimating cell proportions and bMIND for estimating sample-wise cell-type gene expressions. For eight brain cell types, 25,273 cell-type eQTLs were identified with deconvoluted expressions (decon-eQTLs). The results showed that decon-eQTLs explained more schizophrenia GWAS heritability than bulk tissue or single-cell eQTLs did alone. Differential gene expressions associated with Alzheimer's disease, schizophrenia, and brain development were also examined using the deconvoluted data. Our findings, which were replicated in bulk tissue and single-cell data, provided insights into the biological applications of deconvoluted data in multiple brain disorders.

The complexity and heterogeneity of schizophrenia have hindered mechanistic elucidation and the development of more effective therapies. Here, we performed single-cell dissection of schizophrenia-associated transcriptomic changes in the human prefrontal cortex across 140 individuals in two independent cohorts. Excitatory neurons were the most affected cell group, with transcriptional changes converging on neurodevelopment and synapse-related molecular pathways. Transcriptional alterations included known genetic risk factors, suggesting convergence of rare and common genomic variants on neuronal population-specific alterations in schizophrenia. Based on the magnitude of schizophrenia-associated transcriptional change, we identified two populations of individuals with schizophrenia marked by expression of specific excitatory and inhibitory neuronal cell states. This single-cell atlas links transcriptomic changes to etiological genetic risk factors, contextualizing established knowledge within the human cortical cytoarchitecture and facilitating mechanistic understanding of schizophrenia pathophysiology and heterogeneity.

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Enhancers play an essential role in the etiology of schizophrenia; however, the dysregulation of enhancer activity and its impact on the regulome in schizophrenia remains understudied. To address this gap in our knowledge, we assessed enhancer and gene expression in 1,382 brain samples comprising cases with schizophrenia and unaffected controls. Dysregulation of enhancer expression was concordant with changes in gene expression, and was more closely associated with schizophrenia polygenic risk, suggesting that enhancer dysregulation is proximal to the genetic etiology of the disease. Modeling the shared variance of cis-coordinated genes and enhancers revealed a gene regulatory program that was highly associated with genetic vulnerability to schizophrenia. By integrating coordinated factors with evolutionary constraints, we found that enhancers acquired during human evolution are more likely to regulate genes that are implicated in neuropsychiatric disorders and, thus, hold potential as therapeutic targets. Our analysis provides a systematic view of regulome dysregulation in schizophrenia and highlights its convergence with schizophrenia polygenic risk and human-gained enhancers., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

Background: Converging evidence from large-scale genetic and postmortem studies highlights the role of aberrant neurotransmission and genetic regulation in brain-related disorders. However, identifying neuronal activity-regulated transcriptional programs in the human brain and understanding how changes contribute to disease remain challenging., Methods: To better understand how the activity-dependent regulome contributes to risk for brain-related disorders, we profiled the transcriptomic and epigenomic changes following neuronal depolarization in human induced pluripotent stem cell-derived glutamatergic neurons (NGN2) from 6 patients with schizophrenia and 5 control participants., Results: Multiomic data integration associated global patterns of chromatin accessibility with gene expression and identified enhancer-promoter interactions in glutamatergic neurons. Within 1 hour of potassium chloride-induced depolarization, independent of diagnosis, glutamatergic neurons displayed substantial activity-dependent changes in the expression of genes regulating synaptic function. Depolarization-induced changes in the regulome revealed significant heritability enrichment for schizophrenia and Parkinson's disease, adding to mounting evidence that sequence variation within activation-dependent regulatory elements contributes to the genetic risk for brain-related disorders. Gene coexpression network analysis elucidated interactions among activity-dependent and disease-associated genes and pointed to a key driver (NAV3) that interacted with multiple genes involved in axon guidance., Conclusions: Overall, we demonstrated that deciphering the activity-dependent regulome in glutamatergic neurons reveals novel targets for advanced diagnosis and therapy., (Copyright © 2023 Society of Biological Psychiatry. All rights reserved.)

Background: Individuals with Down syndrome (DS) have intellectual disability and develop Alzheimer's disease (AD) pathology during midlife, particularly in the hippocampal component of the medial temporal lobe memory circuit. However, molecular and cellular mechanisms underlying selective vulnerability of hippocampal CA1 neurons remains a major knowledge gap during DS/AD onset. This is compounded by evidence showing spatial (e.g., deep versus superficial) localization of pyramidal neurons (PNs) has profound effects on activity and innervation within the CA1 region., Objective: We investigated whether there is a spatial profiling difference in CA1 PNs in an aged female DS/AD mouse model. We posit dysfunction may be dependent on spatial localization and innervation patterns within discrete CA1 subfields., Methods: Laser capture microdissection was performed on trisomic CA1 PNs in an established mouse model of DS/AD compared to disomic controls, isolating the entire CA1 pyramidal neuron layer and sublayer microisolations of deep and superficial PNs from the distal CA1 (CA1a) region., Results: RNA sequencing and bioinformatic inquiry revealed dysregulation of CA1 PNs based on spatial location and innervation patterns. The entire CA1 region displayed the most differentially expressed genes (DEGs) in trisomic mice reflecting innate DS vulnerability, while trisomic CA1a deep PNs exhibited fewer but more physiologically relevant DEGs, as evidenced by bioinformatic inquiry., Conclusions: CA1a deep neurons displayed numerous DEGs linked to cognitive functions whereas CA1a superficial neurons, with approximately equal numbers of DEGs, were not linked to pathways of dysregulation, suggesting the spatial location of vulnerable CA1 PNs plays an important role in circuit dissolution.

Background: Genotypes are strongly associated with disease phenotypes, particularly in brain disorders. However, the molecular and cellular mechanisms behind this association remain elusive. With emerging multimodal data for these mechanisms, machine learning methods can be applied for phenotype prediction at different scales, but due to the black-box nature of machine learning, integrating these modalities and interpreting biological mechanisms can be challenging. Additionally, the partial availability of these multimodal data presents a challenge in developing these predictive models., Method: To address these challenges, we developed DeepGAMI, an interpretable neural network model to improve genotype-phenotype prediction from multimodal data. DeepGAMI leverages functional genomic information, such as eQTLs and gene regulation, to guide neural network connections. Additionally, it includes an auxiliary learning layer for cross-modal imputation allowing the imputation of latent features of missing modalities and thus predicting phenotypes from a single modality. Finally, DeepGAMI uses integrated gradient to prioritize multimodal features for various phenotypes., Results: We applied DeepGAMI to several multimodal datasets including genotype and bulk and cell-type gene expression data in brain diseases, and gene expression and electrophysiology data of mouse neuronal cells. Using cross-validation and independent validation, DeepGAMI outperformed existing methods for classifying disease types, and cellular and clinical phenotypes, even using single modalities (e.g., AUC score of 0.79 for Schizophrenia and 0.73 for cognitive impairment in Alzheimer's disease)., Conclusion: We demonstrated that DeepGAMI improves phenotype prediction and prioritizes phenotypic features and networks in multiple multimodal datasets in complex brains and brain diseases. Also, it prioritized disease-associated variants, genes, and regulatory networks linked to different phenotypes, providing novel insights into the interpretation of gene regulatory mechanisms. DeepGAMI is open-source and available for general use., (© 2023. The Author(s).)

The fading and quenching of fluorescence intensity has been a major problem in the use of fluorescein isothiocyanate (FITC) for immunofluorescence cytochemical techniques, especially with laser confocal microscopy. The companion article by Longin et al. provided an empirical approach to overcoming this problem. The present commentary highlights the significance of the Longin et al. article when it was published and its continued relevance today.

Attention-deficit hyperactivity disorder (ADHD) is a prevalent neurodevelopmental disorder with a major genetic component. Here, we present a genome-wide association study meta-analysis of ADHD comprising 38,691 individuals with ADHD and 186,843 controls. We identified 27 genome-wide significant loci, highlighting 76 potential risk genes enriched among genes expressed particularly in early brain development. Overall, ADHD genetic risk was associated with several brain-specific neuronal subtypes and midbrain dopaminergic neurons. In exome-sequencing data from 17,896 individuals, we identified an increased load of rare protein-truncating variants in ADHD for a set of risk genes enriched with probable causal common variants, potentially implicating SORCS3 in ADHD by both common and rare variants. Bivariate Gaussian mixture modeling estimated that 84-98% of ADHD-influencing variants are shared with other psychiatric disorders. In addition, common-variant ADHD risk was associated with impaired complex cognition such as verbal reasoning and a range of executive functions, including attention., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Post-acute sequelae of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection are debilitating, clinically heterogeneous and of unknown molecular etiology. A transcriptome-wide investigation was performed in 165 acutely infected hospitalized individuals who were followed clinically into the post-acute period. Distinct gene expression signatures of post-acute sequelae were already present in whole blood during acute infection, with innate and adaptive immune cells implicated in different symptoms. Two clusters of sequelae exhibited divergent plasma-cell-associated gene expression patterns. In one cluster, sequelae associated with higher expression of immunoglobulin-related genes in an anti-spike antibody titer-dependent manner. In the other, sequelae associated independently of these titers with lower expression of immunoglobulin-related genes, indicating lower non-specific antibody production in individuals with these sequelae. This relationship between lower total immunoglobulins and sequelae was validated in an external cohort. Altogether, multiple etiologies of post-acute sequelae were already detectable during SARS-CoV-2 infection, directly linking these sequelae with the acute host response to the virus and providing early insights into their development., (© 2022. The Author(s).)

Identification of risk variants for neuropsychiatric diseases within enhancers underscores the importance of understanding population-level variation in enhancer function in the human brain. Besides regulating tissue-specific and cell-type-specific transcription of target genes, enhancers themselves can be transcribed. By jointly analyzing large-scale cell-type-specific transcriptome and regulome data, we cataloged 30,795 neuronal and 23,265 non-neuronal candidate transcribed enhancers. Examination of the transcriptome in 1,382 brain samples identified robust expression of transcribed enhancers. We explored gene-enhancer coordination and found that enhancer-linked genes are strongly implicated in neuropsychiatric disease. We identified expression quantitative trait loci (eQTLs) for both genes and enhancers and found that enhancer eQTLs mediate a substantial fraction of neuropsychiatric trait heritability. Inclusion of enhancer eQTLs in transcriptome-wide association studies enhanced functional interpretation of disease loci. Overall, our study characterizes the gene-enhancer regulome and genetic mechanisms in the human cortex in both healthy and diseased states., (© 2022. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

To characterize the dysregulation of chromatin accessibility in Alzheimer's disease (AD), we generated 636 ATAC-seq libraries from neuronal and nonneuronal nuclei isolated from the superior temporal gyrus and entorhinal cortex of 153 AD cases and 56 controls. By analyzing a total of ~20 billion read pairs, we expanded the repertoire of known open chromatin regions (OCRs) in the human brain and identified cell-type-specific enhancer-promoter interactions. We show that interindividual variability in OCRs can be leveraged to identify cis-regulatory domains (CRDs) that capture the three-dimensional structure of the genome (3D genome). We identified AD-associated effects on chromatin accessibility, the 3D genome and transcription factor (TF) regulatory networks. For one of the most AD-perturbed TFs, USF2, we validated its regulatory effect on lysosomal genes. Overall, we applied a systematic approach to understanding the role of the 3D genome in AD. We provide all data as an online resource for widespread community-based analysis., (© 2022. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Recent efforts have identified genetic loci that are associated with coronavirus disease 2019 (COVID-19) infection rates and disease outcome severity. Translating these genetic findings into druggable genes that reduce COVID-19 host susceptibility is a critical next step. Using a translational genomics approach that integrates COVID-19 genetic susceptibility variants, multi-tissue genetically regulated gene expression (GReX), and perturbagen signatures, we identified IL10RB as the top candidate gene target for COVID-19 host susceptibility. In a series of validation steps, we show that predicted GReX upregulation of IL10RB and higher IL10RB expression in COVID-19 patient blood is associated with worse COVID-19 outcomes and that in vitro IL10RB overexpression is associated with increased viral load and activation of disease-relevant molecular pathways., (© 2022. The Author(s).)

Microglia are brain myeloid cells that play a critical role in neuroimmunity and the etiology of Alzheimer's disease (AD), yet our understanding of how the genetic regulatory landscape controls microglial function and contributes to AD is limited. Here, we performed transcriptome and chromatin accessibility profiling in primary human microglia from 150 donors to identify genetically driven variation and cell-specific enhancer-promoter (E-P) interactions. Integrative fine-mapping analysis identified putative regulatory mechanisms for 21 AD risk loci, of which 18 were refined to a single gene, including 3 new candidate risk genes (KCNN4, FIBP and LRRC25). Transcription factor regulatory networks captured AD risk variation and identified SPI1 as a key putative regulator of microglia expression and AD risk. This comprehensive resource capturing variation in the human microglia regulome provides insights into the etiology of neurodegenerative disease., (© 2022. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Stressors of different natures induce activation of the hypothalamic-pituitary-adrenal (HPA) axis at different magnitudes. Moreover, the HPA axis response to repeated exposure is usually distinct from that elicited by a single session. Paradoxical sleep deprivation (PSD) augments ACTH and corticosterone (CORT) levels, but the nature of this stimulus is not yet defined. The purpose of the present study was to qualitatively compare the stress response of animals submitted to PSD to that of rats exposed once or four times to cold, as a physiological stress, movement restraint (RST) as a mixed stressor and predator odour (PRED) as the psychological stressor, whilst animals were submitted for 1 or 4 days to PSD and respective control groups. None of the stressors altered corticotropin releasing factor immunoreactivity in the paraventricular nucleus of the hypothalamus (PVN), median eminence (ME) or central amygdala, compared to control groups, whereas vasopressin immunoreactivity in PSD animals was decreased in the PVN and increased in the ME, indicating augmented activity of this system. ACTH levels were higher after repeated stress or prolonged PSD than after single- or 1 day-exposure and control groups, whereas the CORT response was habituated by repeated stress, but not by 4-days PSD. This dissociation resulted in changes in the CORT : ACTH ratio, with repeated cold and RST decreasing the ratio compared to single exposure, but no change was seen in PRED and PSD groups. Comparing the magnitude and pattern of pituitary-adrenal response to the different stressors, PSD-induced responses were closer to that shown by PRED-exposed rats. In contrast, the hypothalamic response of PSD-exposed rats was unique, inasmuch as this was the only stressor which increased the activity of the vasopressin system. In conclusion, we propose that the pituitary-adrenal response to PSD is similar to that induced by a psychological stressor., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Moraes, Machado, Koban, Hoffman and Suchecki.)

The anatomy and morphology of gonadotropin-releasing hormone (GnRH) neurons makes them both a joy and a challenge to investigate. They are a highly unique population of neurons given their developmental migration into the brain from the olfactory placode, their relatively small number, their largely scattered distribution within the rostral forebrain, and, in some species, their highly varied individual anatomical characteristics. These unique features have posed technological hurdles to overcome and promoted fertile ground for the establishment and use of creative approaches. Historical and more contemporary discoveries defining GnRH neuron anatomy remain critical in shaping and challenging our views of GnRH neuron function in the regulation of reproductive function. We begin this review with a historical overview of anatomical discoveries and developing methodologies that have shaped our understanding of the reproductive axis. We then highlight significant discoveries across specific groups of mammalian species to address some of the important comparative aspects of GnRH neuroanatomy. Lastly, we touch on unresolved questions and opportunities for future neuroanatomical research on this fascinating and important population of neurons., (© 2022 The Authors. Journal of Neuroendocrinology published by John Wiley & Sons Ltd on behalf of British Society for Neuroendocrinology.)

Chromosomal organization, scaling from the 147-base pair (bp) nucleosome to megabase-ranging domains encompassing multiple transcriptional units, including heritability loci for psychiatric traits, remains largely unexplored in the human brain. In this study, we constructed promoter- and enhancer-enriched nucleosomal histone modification landscapes for adult prefrontal cortex from H3-lysine 27 acetylation and H3-lysine 4 trimethylation profiles, generated from 388 controls and 351 individuals diagnosed with schizophrenia (SCZ) or bipolar disorder (BD) (n = 739). We mapped thousands of cis-regulatory domains (CRDs), revealing fine-grained, 10 4 -10 6 -bp chromosomal organization, firmly integrated into Hi-C topologically associating domain stratification by open/repressive chromosomal environments and nuclear topography. Large clusters of hyper-acetylated CRDs were enriched for SCZ heritability, with prominent representation of regulatory sequences governing fetal development and glutamatergic neuron signaling. Therefore, SCZ and BD brains show coordinated dysregulation of risk-associated regulatory sequences assembled into kilobase- to megabase-scaling chromosomal domains., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

While large-scale, genome-wide association studies (GWAS) have identified hundreds of loci associated with brain-related traits, identification of the variants, genes and molecular mechanisms underlying these traits remains challenging. Integration of GWAS with expression quantitative trait loci (eQTLs) and identification of shared genetic architecture have been widely adopted to nominate genes and candidate causal variants. However, this approach is limited by sample size, statistical power and linkage disequilibrium. We developed the multivariate multiple QTL approach and performed a large-scale, multi-ancestry eQTL meta-analysis to increase power and fine-mapping resolution. Analysis of 3,983 RNA-sequenced samples from 2,119 donors, including 474 non-European individuals, yielded an effective sample size of 3,154. Joint statistical fine-mapping of eQTL and GWAS identified 329 variant-trait pairs for 24 brain-related traits driven by 204 unique candidate causal variants for 189 unique genes. This integrative analysis identifies candidate causal variants and elucidates potential regulatory mechanisms for genes underlying schizophrenia, bipolar disorder and Alzheimer's disease., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Background: While schizophrenia differs between males and females in the age of onset, symptomatology, and disease course, the molecular mechanisms underlying these differences remain uncharacterized., Methods: To address questions about the sex-specific effects of schizophrenia, we performed a large-scale transcriptome analysis of RNA sequencing data from 437 controls and 341 cases from two distinct cohorts from the CommonMind Consortium., Results: Analysis across the cohorts identified a reproducible gene expression signature of schizophrenia that was highly concordant with previous work. Differential expression across sex was reproducible across cohorts and identified X- and Y-linked genes, as well as those involved in dosage compensation. Intriguingly, the sex expression signature was also enriched for genes involved in neurexin family protein binding and synaptic organization. Differential expression analysis testing a sex-by-diagnosis interaction effect did not identify any genome-wide signature after multiple testing corrections. Gene coexpression network analysis was performed to reduce dimensionality from thousands of genes to dozens of modules and elucidate interactions among genes. We found enrichment of coexpression modules for sex-by-diagnosis differential expression signatures, which were highly reproducible across the two cohorts and involved a number of diverse pathways, including neural nucleus development, neuron projection morphogenesis, and regulation of neural precursor cell proliferation., Conclusions: Overall, our results indicate that the effect size of sex differences in schizophrenia gene expression signatures is small and underscore the challenge of identifying robust sex-by-diagnosis signatures, which will require future analyses in larger cohorts., (Copyright © 2021 Society of Biological Psychiatry. All rights reserved.)

Multisystem inflammatory syndrome in children (MIS-C) presents with fever, inflammation and pathology of multiple organs in individuals under 21 years of age in the weeks following severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Although an autoimmune pathogenesis has been proposed, the genes, pathways and cell types causal to this new disease remain unknown. Here we perform RNA sequencing of blood from patients with MIS-C and controls to find disease-associated genes clustered in a co-expression module annotated to CD56 dim CD57 + natural killer (NK) cells and exhausted CD8 + T cells. A similar transcriptome signature is replicated in an independent cohort of Kawasaki disease (KD), the related condition after which MIS-C was initially named. Probing a probabilistic causal network previously constructed from over 1,000 blood transcriptomes both validates the structure of this module and reveals nine key regulators, including TBX21, a central coordinator of exhausted CD8 + T cell differentiation. Together, this unbiased, transcriptome-wide survey implicates downregulation of NK cells and cytotoxic T cell exhaustion in the pathogenesis of MIS-C., (© 2021. The Author(s).)

Background: Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, has been associated with neurological and neuropsychiatric illness in many individuals. We sought to further our understanding of the relationship between brain tropism, neuro-inflammation, and host immune response in acute COVID-19 cases., Methods: Three brain regions (dorsolateral prefrontal cortex, medulla oblongata, and choroid plexus) from 5 patients with severe COVID-19 and 4 controls were examined. The presence of the virus was assessed by western blot against viral spike protein, as well as viral transcriptome analysis covering > 99% of SARS-CoV-2 genome and all potential serotypes. Droplet-based single-nucleus RNA sequencing (snRNA-seq) was performed in the same samples to examine the impact of COVID-19 on transcription in individual cells of the brain., Results: Quantification of viral spike S1 protein and viral transcripts did not detect SARS-CoV-2 in the postmortem brain tissue. However, analysis of 68,557 single-nucleus transcriptomes from three distinct regions of the brain identified an increased proportion of stromal cells, monocytes, and macrophages in the choroid plexus of COVID-19 patients. Furthermore, differential gene expression, pseudo-temporal trajectory, and gene regulatory network analyses revealed transcriptional changes in the cortical microglia associated with a range of biological processes, including cellular activation, mobility, and phagocytosis., Conclusions: Despite the absence of detectable SARS-CoV-2 in the brain at the time of death, the findings suggest significant and persistent neuroinflammation in patients with acute COVID-19., (© 2021. The Author(s).)

Summary: Large-scale transcriptome studies with multiple samples per individual are widely used to study disease biology. Yet, current methods for differential expression are inadequate for cross-individual testing for these repeated measures designs. Most problematic, we observe across multiple datasets that current methods can give reproducible false-positive findings that are driven by genetic regulation of gene expression, yet are unrelated to the trait of interest. Here, we introduce a statistical software package, dream, that increases power, controls the false positive rate, enables multiple types of hypothesis tests, and integrates with standard workflows. In 12 analyses in 6 independent datasets, dream yields biological insight not found with existing software while addressing the issue of reproducible false-positive findings., Availability and Implementation: Dream is available within the variancePartition Bioconductor package at http://bioconductor.org/packages/variancePartition., Contact: gabriel.hoffman@mssm.edu., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Although clinical and laboratory data have long been used to guide medical practice, this information is rarely integrated with multi-omic data to identify endotypes. We present Merged Affinity Network Association Clustering (MANAclust), a coding-free, automated pipeline enabling integration of categorical and numeric data spanning clinical and multi-omic profiles for unsupervised clustering to identify disease subsets. Using simulations and real-world data from The Cancer Genome Atlas, we demonstrate that MANAclust's feature selection algorithms are accurate and outperform competitors. We also apply MANAclust to a clinically and multi-omically phenotyped asthma cohort. MANAclust identifies clinically and molecularly distinct clusters, including heterogeneous groups of "healthy controls" and viral and allergy-driven subsets of asthmatic subjects. We also find that subjects with similar clinical presentations have disparate molecular profiles, highlighting the need for additional testing to uncover asthma endotypes. This work facilitates data-driven personalized medicine through integration of clinical parameters with multi-omics. MANAclust is freely available at https://bitbucket.org/scottyler892/manaclust/src/master/., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

The availability of high-quality RNA-sequencing and genotyping data of post-mortem brain collections from consortia such as CommonMind Consortium (CMC) and the Accelerating Medicines Partnership for Alzheimer's Disease (AMP-AD) Consortium enable the generation of a large-scale brain cis-eQTL meta-analysis. Here we generate cerebral cortical eQTL from 1433 samples available from four cohorts (identifying >4.1 million significant eQTL for >18,000 genes), as well as cerebellar eQTL from 261 samples (identifying 874,836 significant eQTL for >10,000 genes). We find substantially improved power in the meta-analysis over individual cohort analyses, particularly in comparison to the Genotype-Tissue Expression (GTEx) Project eQTL. Additionally, we observed differences in eQTL patterns between cerebral and cerebellar brain regions. We provide these brain eQTL as a resource for use by the research community. As a proof of principle for their utility, we apply a colocalization analysis to identify genes underlying the GWAS association peaks for schizophrenia and identify a potentially novel gene colocalization with lncRNA RP11-677M14.2 (posterior probability of colocalization 0.975).

Background: Phelan-McDermid syndrome (PMS) is a rare genetic disorder with high risk of autism spectrum disorder (ASD), intellectual disability, and language delay, and is caused by 22q13.3 deletions or mutations in the SHANK3 gene. To date, the molecular and pathway changes resulting from SHANK3 haploinsufficiency in PMS remain poorly understood. Uncovering these mechanisms is critical for understanding pathobiology of PMS and, ultimately, for the development of new therapeutic interventions., Methods: We developed human-induced pluripotent stem cell (hiPSC)-based models of PMS by reprogramming peripheral blood samples from individuals with PMS (n = 7) and their unaffected siblings (n = 6). For each participant, up to three hiPSC clones were generated and differentiated into induced neural progenitor cells (hiPSC-NPCs; n = 39) and induced forebrain neurons (hiPSC-neurons; n = 41). Genome-wide RNA-sequencing was applied to explore transcriptional differences between PMS probands and unaffected siblings., Results: Transcriptome analyses identified 391 differentially expressed genes (DEGs) in hiPSC-NPCs and 82 DEGs in hiPSC-neurons, when comparing cells from PMS probands and unaffected siblings (FDR < 5%). Genes under-expressed in PMS were implicated in Wnt signaling, embryonic development, and protein translation, while over-expressed genes were enriched for pre- and postsynaptic density genes, regulation of synaptic plasticity, and G-protein-gated potassium channel activity. Gene co-expression network analysis identified two modules in hiPSC-neurons that were over-expressed in PMS, implicating postsynaptic signaling and GDP binding, and both modules harbored a significant enrichment of genetic risk loci for developmental delay and intellectual disability. Finally, PMS-associated genes were integrated with other ASD hiPSC transcriptome findings and several points of convergence were identified, indicating altered Wnt signaling and extracellular matrix., Limitations: Given the rarity of the condition, we could not carry out experimental validation in independent biological samples. In addition, functional and morphological phenotypes caused by loss of SHANK3 were not characterized here., Conclusions: This is the largest human neural sample analyzed in PMS. Genome-wide RNA-sequencing in hiPSC-derived neural cells from individuals with PMS revealed both shared and distinct transcriptional signatures across hiPSC-NPCs and hiPSC-neurons, including many genes implicated in risk for ASD, as well as specific neurobiological pathways, including the Wnt pathway.

Structural variants (SVs) contribute to many disorders, yet, functionally annotating them remains a major challenge. Here, we integrate SVs with RNA-sequencing from human post-mortem brains to quantify their dosage and regulatory effects. We show that genic and regulatory SVs exist at significantly lower frequencies than intergenic SVs. Functional impact of copy number variants (CNVs) stems from both the proportion of genic and regulatory content altered and loss-of-function intolerance of the gene. We train a linear model to predict expression effects of rare CNVs and use it to annotate regulatory disruption of CNVs from 14,891 independent genome-sequenced individuals. Pathogenic deletions implicated in neurodevelopmental disorders show significantly more extreme regulatory disruption scores and if rank ordered would be prioritized higher than using frequency or length alone. This work shows the deleteriousness of regulatory SVs, particularly those altering CTCF sites and provides a simple approach for functionally annotating the regulatory consequences of CNVs.

To reveal post-traumatic stress disorder (PTSD) genetic risk influences on tissue-specific gene expression, we use brain and non-brain transcriptomic imputation. We impute genetically regulated gene expression (GReX) in 29,539 PTSD cases and 166,145 controls from 70 ancestry-specific cohorts and identify 18 significant GReX-PTSD associations corresponding to specific tissue-gene pairs. The results suggest substantial genetic heterogeneity based on ancestry, cohort type (military versus civilian), and sex. Two study-wide significant PTSD associations are identified in European and military European cohorts; ZNF140 is predicted to be upregulated in whole blood, and SNRNP35 is predicted to be downregulated in dorsolateral prefrontal cortex, respectively. In peripheral leukocytes from 175 marines, the observed PTSD differential gene expression correlates with the predicted differences for these individuals, and deployment stress produces glucocorticoid-regulated expression changes that include downregulation of both ZNF140 and SNRNP35. SNRNP35 knockdown in cells validates its functional role in U12-intron splicing. Finally, exogenous glucocorticoids in mice downregulate prefrontal Snrnp35 expression., Competing Interests: Declaration of Interests J.W.S. is an unpaid member of the Bipolar/Depression Research Community Advisory Panel of 23andMe. D.J.S. has received research grants and/or consultancy honoraria from Biocodex, Lundbeck, and Sun. R.Y. is a co-inventor of the following patent: “Genes associated with posttraumatic-stress disorder, WO 2010029176 A1.” I.L. has been a consultant for ARMGO Pharmaceuticals, Sunovion Pharmaceuticals, and Trimaran Pharma. K.J.R. has received consulting income from Alkermes and is on scientific advisory boards for Janssen, Verily, and Resilience Therapeutics. He has also received sponsored research support from Takeda and Brainsway. N.P.D. has held a part-time paid position at Cohen Veteran Bioscience. He has also been a consultant for Sunovion Pharmaceuticals., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

Motivation: Identifying correlated epigenetic features and finding differences in correlation between individuals with disease compared to controls can give novel insight into disease biology. This framework has been successful in analysis of gene expression data, but application to epigenetic data has been limited by the computational cost, lack of scalable software and lack of robust statistical tests., Results: Decorate, differential epigenetic correlation test, identifies correlated epigenetic features and finds clusters of features that are differentially correlated between two or more subsets of the data. The software scales to genome-wide datasets of epigenetic assays on hundreds of individuals. We apply decorate to four large-scale datasets of DNA methylation, ATAC-seq and histone modification ChIP-seq., Availability and Implementation: decorate R package is available from https://github.com/GabrielHoffman/decorate., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

The scientific community has searched for years for ways of examining neuronal tissue to track neural activity with reliable anatomical markers for stimulated neuronal activity. Existing studies that focused on hypothalamic systems offer a few options but do not always compare approaches or validate them for dependence on cell firing, leaving the reader uncertain of the benefits and limitations of each method. Thus, in this article, potential markers will be presented and, where possible, placed into perspective in terms of when and how these methods pertain to hypothalamic function. An example of each approach is included. In reviewing the approaches, one is guided through how neurons work, the consequences of their stimulation, and then the potential markers that could be applied to hypothalamic systems are discussed. Approaches will use features of neuronal glucose utilization, water/oxygen movement, changes in neuron-glial interactions, receptor translocation, cytoskeletal changes, stimulus-synthesis coupling that includes expression of the heteronuclear or mature mRNA for transmitters or the enzymes that make them, and changes in transcription factors (immediate early gene products, precursor buildup, use of promoter-driven surrogate proteins, and induced expression of added transmitters. This article includes discussion of methodological limitations and the power of combining approaches to understand neuronal function. © 2020 American Physiological Society. Compr Physiol 10:549-575, 2020., (Copyright © 2020 American Physiological Society. All rights reserved.)

NRXN1 undergoes extensive alternative splicing, and non-recurrent heterozygous deletions in NRXN1 are strongly associated with neuropsychiatric disorders. We establish that human induced pluripotent stem cell (hiPSC)-derived neurons well represent the diversity of NRXN1α alternative splicing observed in the human brain, cataloguing 123 high-confidence in-frame human NRXN1α isoforms. Patient-derived NRXN1 +/- hiPSC-neurons show a greater than twofold reduction in half of the wild-type NRXN1α isoforms and express dozens of novel isoforms from the mutant allele. Reduced neuronal activity in patient-derived NRXN1 +/- hiPSC-neurons is ameliorated by overexpression of individual control isoforms in a genotype-dependent manner, whereas individual mutant isoforms decrease neuronal activity levels in control hiPSC-neurons. In a genotype-dependent manner, the phenotypic impact of patient-specific NRXN1 +/- mutations can occur through a reduction in wild-type NRXN1α isoform levels as well as the presence of mutant NRXN1α isoforms.