Single-cell genomics and regulatory networks for 388 human brains

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, multi-omics datasets into a resource comprising >2.8M 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 >550K cell-type-specific regulatory elements and >1.4M 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.

Alzheimer's disease transcriptional landscape in ex-vivo human microglia

Microglia are resident immune cells of the brain and are implicated in the etiology of Alzheimer's Disease (AD) and other diseases. Yet the cellular and molecular processes regulating their function throughout the course of the disease are poorly understood. Here, we present the transcriptional landscape of primary microglia from 189 human postmortem brains, including 58 healthy aging individuals and 131 with a range of disease phenotypes, including 63 patients representing the full spectrum of clinical and pathological severity of AD. We identified transcriptional changes associated with multiple AD phenotypes, capturing the severity of dementia and neuropathological lesions. Transcript-level analyses identified additional genes with heterogeneous isoform usage and AD phenotypes. We identified changes in gene-gene coordination in AD, dysregulation of co-expression modules, and disease subtypes with distinct gene expression. Taken together, these data further our understanding of the key role of microglia in AD biology and nominate candidates for therapeutic intervention.

Activity-Dependent Transcriptional Program in NGN2+ Neurons Enriched for Genetic Risk for Brain-Related Disorders

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.

Convergence of the dysregulated regulome in schizophrenia with polygenic risk and evolutionarily constrained enhancers

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.

Transcriptional correlates of cocaine-associated learning in striatal ARC ensembles

Learned associations between the rewarding effects of drugs and the context in which they are experienced underlie context-induced relapse. Previous work demonstrates the importance of sparse neuronal populations - called neuronal ensembles - in associative learning and cocaine seeking, but it remains unknown whether the encoding vs. retrieval of cocaine-associated memories involves similar or distinct mechanisms of ensemble activation and reactivation in nucleus accumbens (NAc). We use ArcCreER T2 mice to establish that mostly distinct NAc ensembles are recruited by initial vs. repeated exposures to cocaine, which are then differentially reactivated and exert distinct effects during cocaine-related memory retrieval. Single-nuclei RNA-sequencing of these ensembles demonstrates predominant recruitment of D1 medium spiny neurons and identifies transcriptional properties that are selective to cocaine-recruited NAc neurons and could explain distinct excitability features. These findings fundamentally advance our understanding of how cocaine drives pathological memory formation during repeated exposures.

QuadST: A Powerful and Robust Approach for Identifying Cell-Cell Interaction-Changed Genes on Spatially Resolved Transcriptomics

Spatially resolved transcriptomics (SRT) have enabled profiling spatial organization of cells and their transcriptome in situ. Various analytical methods have been developed to uncover cell-cell interaction processes using SRT data. To improve upon existing efforts, we developed a novel statistical framework called QuadST for the robust and powerful identification of interaction-changed genes (ICGs) for cell-type-pair specific interactions on a single-cell SRT dataset. QuadST is motivated by the idea that in the presence of cell-cell interaction, gene expression level can vary with cell-cell distance between cell type pairs, which can be particularly pronounced within and in the vicinity of cell-cell interaction distance. Specifically, QuadST infers ICGs in a specific cell type pair's interaction based on a quantile regression model, which allows us to assess the strength of distance-expression association across entire distance quantiles conditioned on gene expression level. To identify ICGs, QuadST performs a hypothesis testing with an empirically estimated FDR, whose upper bound is determined by the ratio of cumulative associations at symmetrically smaller and larger distance quantiles simultaneously across all genes. Simulation studies illustrate that QuadST provides consistent FDR control and better power performance than other compared methods. Its application on SRT datasets profiled from mouse brains demonstrates that QuadST can identify ICGs presumed to play a role in specific cell type pair interactions (e.g., synaptic pathway genes among excitatory neuron cell interactions). These results suggest that QuadST can be a useful tool to discover genes and regulatory processes involved in specific cell type pair interactions.

Long-read RNA-seq atlas of novel microglia isoforms elucidates disease-associated genetic regulation of splicing

Microglia, the innate immune cells of the central nervous system, have been genetically implicated in multiple neurodegenerative diseases. We previously mapped the genetic regulation of gene expression and mRNA splicing in human microglia, identifying several loci where common genetic variants in microglia-specific regulatory elements explain disease risk loci identified by GWAS. However, identifying genetic effects on splicing has been challenging due to the use of short sequencing reads to identify causal isoforms. Here we present the isoform-centric microglia genomic atlas (isoMiGA) which leverages the power of long-read RNA-seq to identify 35,879 novel microglia isoforms. We show that the novel microglia isoforms are involved in stimulation response and brain region specificity. We then quantified the expression of both known and novel isoforms in a multi-ethnic meta-analysis of 555 human microglia short-read RNA-seq samples from 391 donors, the largest to date, and found associations with genetic risk loci in Alzheimer's disease and Parkinson's disease. We nominate several loci that may act through complex changes in isoform and splice site usage.

Lineage specific 3D genome structure in the adult human brain and neurodevelopmental changes in the chromatin interactome

The human brain is a complex organ comprised of distinct cell types, and the contribution of the 3D genome to lineage specific gene expression remains poorly understood. To decipher cell type specific genome architecture, and characterize fine scale changes in the chromatin interactome across neural development, we compared the 3D genome of the human fetal cortical plate to that of neurons and glia isolated from the adult prefrontal cortex. We found that neurons have weaker genome compartmentalization compared to glia, but stronger TADs, which emerge during fetal development. Furthermore, relative to glia, the neuronal genome shifts more strongly towards repressive compartments. Neurons have differential TAD boundaries that are proximal to active promoters involved in neurodevelopmental processes. CRISPRi on CNTNAP2 in hIPSC-derived neurons reveals that transcriptional inactivation correlates with loss of insulation at the differential boundary. Finally, re-wiring of chromatin loops during neural development is associated with transcriptional and functional changes. Importantly, differential loops in the fetal cortex are associated with autism GWAS loci, suggesting a neuropsychiatric disease mechanism affecting the chromatin interactome. Furthermore, neural development involves gaining enhancer-promoter loops that upregulate genes that control synaptic activity. Altogether, our study provides multi-scale insights on the 3D genome in the human brain.

DeepGAMI: deep biologically guided auxiliary learning for multimodal integration and imputation to improve genotype-phenotype prediction

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.

A public resource of single cell transcriptomes and multiscale networks from persons with and without Alzheimer's disease

The emergence of technologies that can support high-throughput profiling of single cell transcriptomes offers to revolutionize the study of brain tissue from persons with and without Alzheimer's disease (AD). Integration of these data with additional complementary multiomics data such as genetics, proteomics and clinical data provides powerful opportunities to link observed cell subpopulations and molecular network features within a broader disease-relevant context. We report here single nucleus RNA sequencing (snRNA-seq) profiles generated from superior frontal gyrus cortical tissue samples from 101 exceptionally well characterized, aged subjects from the Banner Brain and Body Donation Program in combination with whole genome sequences. We report findings that link common AD risk variants with CR1 expression in oligodendrocytes as well as alterations in peripheral hematological lab parameters, with these observations replicated in an independent, prospective cohort study of ageing and dementia. We also observed an AD-associated CD83(+) microglial subtype with unique molecular networks that encompass many known regulators of AD-relevant microglial biology, and which are associated with immunoglobulin IgG4 production in the transverse colon. These findings illustrate the power of multi-tissue molecular profiling to contextualize snRNA-seq brain transcriptomics and reveal novel disease biology. The transcriptomic, genetic, phenotypic, and network data resources described within this study are available for access and utilization by the scientific community.

Multi-omic profiling of the developing human cerebral cortex at the single-cell level

The cellular complexity of the human brain is established via dynamic changes in gene expression throughout development that is mediated, in part, by the spatiotemporal activity of cis-regulatory elements (CREs). We simultaneously profiled gene expression and chromatin accessibility in 45,549 cortical nuclei across six broad developmental time points from fetus to adult. We identified cell type-specific domains in which chromatin accessibility is highly correlated with gene expression. Differentiation pseudotime trajectory analysis indicates that chromatin accessibility at CREs precedes transcription and that dynamic changes in chromatin structure play a critical role in neuronal lineage commitment. In addition, we mapped cell type-specific and temporally specific genetic loci implicated in neuropsychiatric traits, including schizophrenia and bipolar disorder. Together, our results describe the complex regulation of cell composition at critical stages in lineage determination and shed light on the impact of spatiotemporal alterations in gene expression on neuropsychiatric disease.

The neuronal chromatin landscape in adult schizophrenia brains is linked to early fetal development

Non-coding variants increase risk of neuropsychiatric disease. However, our understanding of the cell-type specific role of the non-coding genome in disease is incomplete. We performed population scale (N=1,393) chromatin accessibility profiling of neurons and non-neurons from two neocortical brain regions: the anterior cingulate cortex and dorsolateral prefrontal cortex. Across both regions, we observed notable differences in neuronal chromatin accessibility between schizophrenia cases and controls. A per-sample disease pseudotime was positively associated with genetic liability for schizophrenia. Organizing chromatin into cis- and trans-regulatory domains, identified a prominent neuronal trans-regulatory domain (TRD1) active in immature glutamatergic neurons during fetal development. Polygenic risk score analysis using genetic variants within chromatin accessibility of TRD1 successfully predicted susceptibility to schizophrenia in the Million Veteran Program cohort. Overall, we present the most extensive resource to date of chromatin accessibility in the human cortex, yielding insights into the cell-type specific etiology of schizophrenia.

The neuronal chromatin landscape in adult schizophrenia brains is linked to early fetal development

Non-coding variants increase risk of neuropsychiatric disease. However, our understanding of the cell-type specific role of the non-coding genome in disease is incomplete. We performed population scale (N=1,393) chromatin accessibility profiling of neurons and non-neurons from two neocortical brain regions: the anterior cingulate cortex and dorsolateral prefrontal cortex. Across both regions, we observed notable differences in neuronal chromatin accessibility between schizophrenia cases and controls. A per-sample disease pseudotime was positively associated with genetic liability for schizophrenia. Organizing chromatin into cis- and trans-regulatory domains, identified a prominent neuronal trans-regulatory domain (TRD1) active in immature glutamatergic neurons during fetal development. Polygenic risk score analysis using genetic variants within chromatin accessibility of TRD1 successfully predicted susceptibility to schizophrenia in the Million Veteran Program cohort. Overall, we present the most extensive resource to date of chromatin accessibility in the human cortex, yielding insights into the cell-type specific etiology of schizophrenia.

Neuroimmune Mechanisms of Opioid Use Disorder and Recovery: Translatability to Human Studies, and Future Research Directions

Opioid use disorder (OUD) is a major current cause of morbidity and mortality. Long-term exposure to short-acting opioids (MOP-r agonists such as heroin or fentanyl) results in complex pathophysiological changes to neuroimmune and neuroinflammatory functions, affected in part by peripheral mechanisms (e.g., cytokines in blood), and by neuroendocrine systems such as the hypothalamic-pituitary-adrenal (HPA) stress axis. There are important findings from preclinical models, but their role in the trajectory and outcomes of OUD in humans is not well understood. The goal of this narrative review is to examine available data on immune and inflammatory functions in persons with OUD, and to identify major areas for future research. Peripheral blood biomarker studies revealed a pro-inflammatory state in persons with OUD in withdrawal or early abstinence, consistent with available postmortem brain studies (which show glial activation) and diffusion tensor imaging studies (indicating white matter disruptions), with gradual abstinence-associated recovery. The mechanistic roles of these neuroimmune and neuroinflammatory changes in the trajectory of OUD (including recovery and medication management) cannot be examined practically with postmortem data. Collection of longitudinal data in larger-scale human cohorts would allow examination of these mechanisms associated with OUD stage and progression. Given the heterogeneity in presentation of OUD, a precision medicine approach integrating multi-omic peripheral biomarkers and comprehensive phenotyping, including neuroimaging, can be beneficial in risk stratification, and individually optimized selection of interventions for individuals who will benefit, and assessments under refractory therapy.

Multi-omic atlas of the parahippocampal gyrus in Alzheimer's disease

Alzheimer's disease (AD) is the most common form of dementia worldwide, with a projection of 151 million cases by 2050. Previous genetic studies have identified three main genes associated with early-onset familial Alzheimer's disease, however this subtype accounts for less than 5% of total cases. Next-generation sequencing has been well established and holds great promise to assist in the development of novel therapeutics as well as biomarkers to prevent or slow the progression of this devastating disease. Here we present a public resource of functional genomic data from the parahippocampal gyrus of 201 postmortem control, mild cognitively impaired (MCI) and AD individuals from the Mount Sinai brain bank, of which whole-genome sequencing (WGS), and bulk RNA sequencing (RNA-seq) were previously published. The genomic data include bulk proteomics and DNA methylation, as well as cell-type-specific RNA-seq and assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) data. We have performed extensive preprocessing and quality control, allowing the research community to access and utilize this public resource available on the Synapse platform at https://doi.org/10.7303/syn51180043.2 .

Genome-wide analysis of a model-derived binge eating disorder phenotype identifies risk loci and implicates iron metabolism

Binge eating disorder (BED) is the most common eating disorder, yet its genetic architecture remains largely unknown. Studying BED is challenging because it is often comorbid with obesity, a common and highly polygenic trait, and it is underdiagnosed in biobank data sets. To address this limitation, we apply a supervised machine-learning approach (using 822 cases of individuals diagnosed with BED) to estimate the probability of each individual having BED based on electronic medical records from the Million Veteran Program. We perform a genome-wide association study of individuals of African (n = 77,574) and European (n = 285,138) ancestry while controlling for body mass index to identify three independent loci near the HFE, MCHR2 and LRP11 genes and suggest APOE as a risk gene for BED. We identify shared heritability between BED and several neuropsychiatric traits, and implicate iron metabolism in the pathophysiology of BED. Overall, our findings provide insights into the genetics underlying BED and suggest directions for future translational research.

The stability of the myelinating oligodendrocyte transcriptome is regulated by the nuclear lamina

Oligodendrocytes are specialized cells that insulate and support axons with their myelin membrane, allowing proper brain function. Here, we identify lamin A/C (LMNA/C) as essential for transcriptional and functional stability of myelinating oligodendrocytes. We show that LMNA/C levels increase with differentiation of progenitors and that loss of Lmna in differentiated oligodendrocytes profoundly alters their chromatin accessibility and transcriptional signature. Lmna deletion in myelinating glia is compatible with normal developmental myelination. However, altered chromatin accessibility is detected in fully differentiated oligodendrocytes together with increased expression of progenitor genes and decreased levels of lipid-related transcription factors and inner mitochondrial membrane transcripts. These changes are accompanied by altered brain metabolism, lower levels of myelin-related lipids, and altered mitochondrial structure in oligodendrocytes, thereby resulting in myelin thinning and the development of a progressively worsening motor phenotype. Overall, our data identify LMNA/C as essential for maintaining the transcriptional and functional stability of myelinating oligodendrocytes.

Genome-wide methylomic regulation of multiscale gene networks in Alzheimer's disease

INTRODUCTION

Recent studies revealed the association of abnormal methylomic changes with Alzheimer's disease (AD) but there is a lack of systematic study of the impact of methylomic alterations over the molecular networks underlying AD.

METHODS

We profiled genome-wide methylomic variations in the parahippocampal gyrus from 201 post mortem control, mild cognitive impaired, and AD brains.

RESULTS

We identified 270 distinct differentially methylated regions (DMRs) associated with AD. We quantified the impact of these DMRs on each gene and each protein as well as gene and protein co-expression networks. DNA methylation had a profound impact on both AD-associated gene/protein modules and their key regulators. We further integrated the matched multi-omics data to show the impact of DNA methylation on chromatin accessibility, which further modulates gene and protein expression.

DISCUSSION

The quantified impact of DNA methylation on gene and protein networks underlying AD identified potential upstream epigenetic regulators of AD.

HIGHLIGHTS

A cohort of DNA methylation data in the parahippocampal gyrus was developed from 201 post mortem control, mild cognitive impaired, and Alzheimer's disease (AD) brains. Two hundred seventy distinct differentially methylated regions (DMRs) were found to be associated with AD compared to normal control. A metric was developed to quantify methylation impact on each gene and each protein. DNA methylation was found to have a profound impact on not only the AD-associated gene modules but also key regulators of the gene and protein networks. Key findings were validated in an independent multi-omics cohort in AD. The impact of DNA methylation on chromatin accessibility was also investigated by integrating the matched methylomic, epigenomic, transcriptomic, and proteomic data.

Dopaminoceptive D1 and D2 neurons in ventral hippocampus arbitrate approach and avoidance in anxiety

The hippocampus 1-7, as well as dopamine circuits 8-11, coordinate decision-making in anxiety-eliciting situations. Yet, little is known about how dopamine modulates hippocampal representations of emotionally-salient stimuli to inform appropriate resolution of approach versus avoidance conflicts. We here study dopaminoceptive neurons in mouse ventral hippocampus (vHipp), molecularly distinguished by their expression of dopamine D1 or D2 receptors. We show that these neurons are transcriptionally distinct and topographically organized across vHipp subfields and cell types. In the ventral subiculum where they are enriched, both D1 and D2 neurons are recruited during anxiogenic exploration, yet with distinct profiles related to investigation and behavioral selection. In turn, they mediate opposite approach/avoidance responses, and are differentially modulated by dopaminergic transmission in that region. Together, these results suggest that vHipp dopamine dynamics gate exploratory behaviors under contextual uncertainty, implicating dopaminoception in the complex computation engaged in vHipp to govern emotional states.

Interleukin-3 coordinates glial-peripheral immune crosstalk to incite multiple sclerosis

Glial cells and central nervous system (CNS)-infiltrating leukocytes contribute to multiple sclerosis (MS). However, the networks that govern crosstalk among these ontologically distinct populations remain unclear. Here, we show that, in mice and humans, CNS-resident astrocytes and infiltrating CD44hiCD4+ T cells generated interleukin-3 (IL-3), while microglia and recruited myeloid cells expressed interleukin-3 receptor-ɑ (IL-3Rɑ). Astrocytic and T cell IL-3 elicited an immune migratory and chemotactic program by IL-3Rɑ+ myeloid cells that enhanced CNS immune cell infiltration, exacerbating MS and its preclinical model. Multiregional snRNA-seq of human CNS tissue revealed the appearance of IL3RA-expressing myeloid cells with chemotactic programming in MS plaques. IL3RA expression by plaque myeloid cells and IL-3 amount in the cerebrospinal fluid predicted myeloid and T cell abundance in the CNS and correlated with MS severity. Our findings establish IL-3:IL-3RA as a glial-peripheral immune network that prompts immune cell recruitment to the CNS and worsens MS.

Depression pathophysiology, risk prediction of recurrence and comorbid psychiatric disorders using genome-wide analyses

Depression is a common psychiatric disorder and a leading cause of disability worldwide. Here we conducted a genome-wide association study meta-analysis of six datasets, including >1.3 million individuals (371,184 with depression) and identified 243 risk loci. Overall, 64 loci were new, including genes encoding glutamate and GABA receptors, which are targets for antidepressant drugs. Intersection with functional genomics data prioritized likely causal genes and revealed new enrichment of prenatal GABAergic neurons, astrocytes and oligodendrocyte lineages. We found depression to be highly polygenic, with ~11,700 variants explaining 90% of the single-nucleotide polymorphism heritability, estimating that >95% of risk variants for other psychiatric disorders (anxiety, schizophrenia, bipolar disorder and attention deficit hyperactivity disorder) were influencing depression risk when both concordant and discordant variants were considered, and nearly all depression risk variants influenced educational attainment. Additionally, depression genetic risk was associated with impaired complex cognition domains. We dissected the genetic and clinical heterogeneity, revealing distinct polygenic architectures across subgroups of depression and demonstrating significantly increased absolute risks for recurrence and psychiatric comorbidity among cases of depression with the highest polygenic burden, with considerable sex differences. The risks were up to 5- and 32-fold higher than cases with the lowest polygenic burden and the background population, respectively. These results deepen the understanding of the biology underlying depression, its disease progression and inform precision medicine approaches to treatment.

Diversity and Scale: Genetic Architecture of 2,068 Traits in the VA Million Veteran Program

Genome-wide association studies (GWAS) have underrepresented individuals from non-European populations, impeding progress in characterizing the genetic architecture and consequences of health and disease traits. To address this, we present a population-stratified phenome-wide GWAS followed by a multi-population meta-analysis for 2,068 traits derived from electronic health records of 635,969 participants in the Million Veteran Program (MVP), a longitudinal cohort study of diverse U.S. Veterans genetically similar to the respective African (121,177), Admixed American (59,048), East Asian (6,702), and European (449,042) superpopulations defined by the 1000 Genomes Project. We identified 38,270 independent variants associating with one or more traits at experiment-wide P < 4.6 × 10 - 11 significance; fine-mapping 6,318 signals identified from 613 traits to single-variant resolution. Among these, a third (2,069) of the associations were found only among participants genetically similar to non-European reference populations, demonstrating the importance of expanding diversity in genetic studies. Our work provides a comprehensive atlas of phenome-wide genetic associations for future studies dissecting the architecture of complex traits in diverse populations.

Basal forebrain cholinergic neurons are vulnerable in a mouse model of Down syndrome and their molecular fingerprint is rescued by maternal choline supplementation

Basal forebrain cholinergic neuron (BFCN) degeneration is a hallmark of Down syndrome (DS) and Alzheimer's disease (AD). Current therapeutics in these disorders have been unsuccessful in slowing disease progression, likely due to poorly understood complex pathological interactions and dysregulated pathways. The Ts65Dn trisomic mouse model recapitulates both cognitive and morphological deficits of DS and AD, including BFCN degeneration and has shown lifelong behavioral changes due to maternal choline supplementation (MCS). To test the impact of MCS on trisomic BFCNs, we performed laser capture microdissection to individually isolate choline acetyltransferase-immunopositive neurons in Ts65Dn and disomic littermates, in conjunction with MCS at the onset of BFCN degeneration. We utilized single population RNA sequencing (RNA-seq) to interrogate transcriptomic changes within medial septal nucleus (MSN) BFCNs. Leveraging multiple bioinformatic analysis programs on differentially expressed genes (DEGs) by genotype and diet, we identified key canonical pathways and altered physiological functions within Ts65Dn MSN BFCNs, which were attenuated by MCS in trisomic offspring, including the cholinergic, glutamatergic and GABAergic pathways. We linked differential gene expression bioinformatically to multiple neurological functions, including motor dysfunction/movement disorder, early onset neurological disease, ataxia and cognitive impairment via Ingenuity Pathway Analysis. DEGs within these identified pathways may underlie aberrant behavior in the DS mice, with MCS attenuating the underlying gene expression changes. We propose MCS ameliorates aberrant BFCN gene expression within the septohippocampal circuit of trisomic mice through normalization of principally the cholinergic, glutamatergic, and GABAergic signaling pathways, resulting in attenuation of underlying neurological disease functions.

ZBTB7A regulates MDD-specific chromatin signatures and astrocyte-mediated stress vulnerability in orbitofrontal cortex

Hyperexcitability in the orbitofrontal cortex (OFC) is a key clinical feature of anhedonic domains of Major Depressive Disorder (MDD). However, the cellular and molecular substrates underlying this dysfunction remain unknown. Here, cell-population-specific chromatin accessibility profiling in human OFC unexpectedly mapped genetic risk for MDD exclusively to non-neuronal cells, and transcriptomic analyses revealed significant glial dysregulation in this region. Characterization of MDD-specific cis-regulatory elements identified ZBTB7A - a transcriptional regulator of astrocyte reactivity - as an important mediator of MDD-specific chromatin accessibility and gene expression. Genetic manipulations in mouse OFC demonstrated that astrocytic Zbtb7a is both necessary and sufficient to promote behavioral deficits, cell-type-specific transcriptional and chromatin profiles, and OFC neuronal hyperexcitability induced by chronic stress - a major risk factor for MDD. These data thus highlight a critical role for OFC astrocytes in stress vulnerability and pinpoint ZBTB7A as a key dysregulated factor in MDD that mediates maladaptive astrocytic functions driving OFC hyperexcitability.

Efficient differential expression analysis of large-scale single cell transcriptomics data using dreamlet

Advances in single-cell and -nucleus transcriptomics have enabled generation of increasingly large-scale datasets from hundreds of subjects and millions of cells. These studies promise to give unprecedented insight into the cell type specific biology of human disease. Yet performing differential expression analyses across subjects remains difficult due to challenges in statistical modeling of these complex studies and scaling analyses to large datasets. Our open-source R package dreamlet (DiseaseNeurogenomics.github.io/dreamlet) uses a pseudobulk approach based on precision-weighted linear mixed models to identify genes differentially expressed with traits across subjects for each cell cluster. Designed for data from large cohorts, dreamlet is substantially faster and uses less memory than existing workflows, while supporting complex statistical models and controlling the false positive rate. We demonstrate computational and statistical performance on published datasets, and a novel dataset of 1.4M single nuclei from postmortem brains of 150 Alzheimer's disease cases and 149 controls.

Applications of artificial intelligence-machine learning for detection of stress: a critical overview

Psychological distress is a major contributor to human physiology and pathophysiology, and it has been linked to several conditions, such as auto-immune diseases, metabolic syndrome, sleep disorders, and suicidal thoughts and inclination. Therefore, early detection and management of chronic stress is crucial for the prevention of several diseases. Artificial intelligence (AI) and Machine Learning (ML) have promoted a paradigm shift in several areas of biomedicine including diagnosis, monitoring, and prognosis of disease. Here, our review aims to present some of the AI and ML applications for solving biomedical issues related to psychological stress. We provide several lines of evidence from previous studies highlighting that AI and ML have been able to predict stress and detect the brain normal states vs. abnormal states (notably, in post-traumatic stress disorder (PTSD)) with accuracy around 90%. Of note, AI/ML-driven technology applied to identify ubiquitously present stress exposure may not reach its full potential, unless future analytics focus on detecting prolonged distress through such technology instead of merely assessing stress exposure. Moving forward, we propose that a new subcategory of AI methods called Swarm Intelligence (SI) can be used towards detecting stress and PTSD. SI involves ensemble learning techniques to efficiently solve a complex problem, such as stress detection, and it offers particular strength in clinical settings, such as privacy preservation. We posit that AI and ML approaches will be beneficial for the medical and patient community when applied to predict and assess stress levels. Last, we encourage additional research to bring AI and ML into the standard clinical practice for diagnostics in the not-too-distant future.

Efficient differential expression analysis of large-scale single cell transcriptomics data using dreamlet

Advances in single-cell and -nucleus transcriptomics have enabled generation of increasingly large-scale datasets from hundreds of subjects and millions of cells. These studies promise to give unprecedented insight into the cell type specific biology of human disease. Yet performing differential expression analyses across subjects remains difficult due to challenges in statistical modeling of these complex studies and scaling analyses to large datasets. Our open-source R package dreamlet (DiseaseNeurogenomics.github.io/dreamlet) uses a pseudobulk approach based on precision-weighted linear mixed models to identify genes differentially expressed with traits across subjects for each cell cluster. Designed for data from large cohorts, dreamlet is substantially faster and uses less memory than existing workflows, while supporting complex statistical models and controlling the false positive rate. We demonstrate computational and statistical performance on published datasets, and a novel dataset of 1.4M single nuclei from postmortem brains of 150 Alzheimer's disease cases and 149 controls.

Evaluating performance and applications of sample-wise cell deconvolution methods on human brain transcriptomic data

Sample-wise deconvolution methods have been developed to estimate cell-type proportions and gene expressions in bulk-tissue samples. However, the performance of these methods and their biological applications has not been evaluated, particularly on human brain transcriptomic data. Here, nine deconvolution methods were evaluated with sample-matched data from bulk-tissue RNAseq, single-cell/nuclei (sc/sn) RNAseq, and immunohistochemistry. A total of 1,130,767 nuclei/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 expression. 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 alone. Differential gene expression associated with multiple phenotypes were also examined using the deconvoluted data. Our findings, which were replicated in bulk-tissue RNAseq and sc/snRNAseq data, provided new insights into the biological applications of deconvoluted data.

Genetic regulation of cell-type specific chromatin accessibility shapes the etiology of brain diseases

Nucleotide variants in cell type-specific gene regulatory elements in the human brain are major risk factors of human disease. We measured chromatin accessibility in sorted neurons and glia from 1,932 samples of human postmortem brain and identified 34,539 open chromatin regions with chromatin accessibility quantitative trait loci (caQTL). Only 10.4% of caQTL are shared between neurons and glia, supporting the cell type specificity of genetic regulation of the brain regulome. Incorporating allele specific chromatin accessibility improves statistical fine-mapping and refines molecular mechanisms underlying disease risk. Using massively parallel reporter assays in induced excitatory neurons, we screened 19,893 brain QTLs, identifying the functional impact of 476 regulatory variants. Combined, this comprehensive resource captures variation in the human brain regulome and provides novel insights into brain disease etiology.

African ancestry GWAS of dementia in a large military cohort identifies significant risk loci

While genome wide association studies (GWASs) of Alzheimer's Disease (AD) in European (EUR) ancestry cohorts have identified approximately 83 potentially independent AD risk loci, progress in non-European populations has lagged. In this study, data from the Million Veteran Program (MVP), a biobank which includes genetic data from more than 650,000 US Veteran participants, was used to examine dementia genetics in an African descent (AFR) cohort. A GWAS of Alzheimer's disease and related dementias (ADRD), an expanded AD phenotype including dementias such as vascular and non-specific dementia that included 4012 cases and 18,435 controls age 60+ in AFR MVP participants was performed. A proxy dementia GWAS based on survey-reported parental AD or dementia (n = 4385 maternal cases, 2256 paternal cases, and 45,970 controls) was also performed. These two GWASs were meta-analyzed, and then subsequently compared and meta-analyzed with the results from a previous AFR AD GWAS from the Alzheimer's Disease Genetics Consortium (ADGC). A meta-analysis of common variants across the MVP ADRD and proxy GWASs yielded GWAS significant associations in the region of APOE (p = 2.48 × 10-101), in ROBO1 (rs11919682, p = 1.63 × 10-8), and RNA RP11-340A13.2 (rs148433063, p = 8.56 × 10-9). The MVP/ADGC meta-analysis yielded additional significant SNPs near known AD risk genes TREM2 (rs73427293, p = 2.95 × 10-9), CD2AP (rs7738720, p = 1.14 × 10-9), and ABCA7 (rs73505251, p = 3.26 × 10-10), although the peak variants observed in these genes differed from those previously reported in EUR and AFR cohorts. Of the genes in or near suggestive or genome-wide significant associated variants, nine (CDA, SH2D5, DCBLD1, EML6, GOPC, ABCA7, ROS1, TMCO4, and TREM2) were differentially expressed in the brains of AD cases and controls. This represents the largest AFR GWAS of AD and dementia, finding non-APOE GWAS-significant common SNPs associated with dementia. Increasing representation of AFR participants is an important priority in genetic studies and may lead to increased insight into AD pathophysiology and reduce health disparities.

Schizophrenia risk conferred by rare protein-truncating variants is conserved across diverse human populations

Schizophrenia (SCZ) is a chronic mental illness and among the most debilitating conditions encountered in medical practice. A recent landmark SCZ study of the protein-coding regions of the genome identified a causal role for ten genes and a concentration of rare variant signals in evolutionarily constrained genes1. This recent study-and most other large-scale human genetics studies-was mainly composed of individuals of European (EUR) ancestry, and the generalizability of the findings in non-EUR populations remains unclear. To address this gap, we designed a custom sequencing panel of 161 genes selected based on the current knowledge of SCZ genetics and sequenced a new cohort of 11,580 SCZ cases and 10,555 controls of diverse ancestries. Replicating earlier work, we found that cases carried a significantly higher burden of rare protein-truncating variants (PTVs) among evolutionarily constrained genes (odds ratio = 1.48; P = 5.4 × 10-6). In meta-analyses with existing datasets totaling up to 35,828 cases and 107,877 controls, this excess burden was largely consistent across five ancestral populations. Two genes (SRRM2 and AKAP11) were newly implicated as SCZ risk genes, and one gene (PCLO) was identified as shared by individuals with SCZ and those with autism. Overall, our results lend robust support to the rare allelic spectrum of the genetic architecture of SCZ being conserved across diverse human populations.

Genome-wide analyses of ADHD identify 27 risk loci, refine the genetic architecture and implicate several cognitive domains

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.

Detecting and Adjusting for Hidden Biases due to Phenotype Misclassification in Genome-Wide Association Studies

With the advent of healthcare-based genotyped biobanks, genome-wide association studies (GWAS) leverage larger sample sizes, incorporate patients with diverse ancestries and introduce noisier phenotypic definitions. Yet the extent and impact of phenotypic misclassification on large-scale datasets is not currently well understood due to a lack of statistical methods to estimate relevant parameters from empirical data. Here, we develop a statistical method and scalable software, PheMED, Phenotypic Measurement of Effective Dilution, to quantify phenotypic misclassification across GWAS using only summary statistics. We illustrate how the parameters estimated by PheMED relate to the negative and positive predictive value of the labeled phenotype, compared to ground truth, and how misclassification of the phenotype yields diluted effect-sizes of variant-phenotype associations. Furthermore, we apply our methodology to detect multiple instances of statistically significant dilution in real-world data. We demonstrate how effective dilution biases downstream GWAS replication and heritability analyses despite utilizing current best practices, and provide a dilution-aware meta-analysis approach that outperforms existing methods. Consequently, we anticipate that PheMED will be a valuable tool for researchers to address phenotypic data quality issues both within and across cohorts.

Genome-wide association study of school grades identifies genetic overlap between language ability, psychopathology and creativity

Cognitive functions of individuals with psychiatric disorders differ from that of the general population. Such cognitive differences often manifest early in life as differential school performance and have a strong genetic basis. Here we measured genetic predictors of school performance in 30,982 individuals in English, Danish and mathematics via a genome-wide association study (GWAS) and studied their relationship with risk for six major psychiatric disorders. When decomposing the school performance into math and language-specific performances, we observed phenotypically and genetically a strong negative correlation between math performance and risk for most psychiatric disorders. But language performance correlated positively with risk for certain disorders, especially schizophrenia, which we replicate in an independent sample (n = 4547). We also found that the genetic variants relating to increased risk for schizophrenia and better language performance are overrepresented in individuals involved in creative professions (n = 2953) compared to the general population (n = 164,622). The findings together suggest that language ability, creativity and psychopathology might stem from overlapping genetic roots.

Molecular states during acute COVID-19 reveal distinct etiologies of long-term sequelae

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.

Rescue of deficits by Brwd1 copy number restoration in the Ts65Dn mouse model of Down syndrome

With an incidence of ~1 in 800 births, Down syndrome (DS) is the most common chromosomal condition linked to intellectual disability worldwide. While the genetic basis of DS has been identified as a triplication of chromosome 21 (HSA21), the genes encoded from HSA21 that directly contribute to cognitive deficits remain incompletely understood. Here, we found that the HSA21-encoded chromatin effector, BRWD1, was upregulated in neurons derived from iPS cells from an individual with Down syndrome and brain of trisomic mice. We showed that selective copy number restoration of Brwd1 in trisomic animals rescued deficits in hippocampal LTP, cognition and gene expression. We demonstrated that Brwd1 tightly binds the BAF chromatin remodeling complex, and that increased Brwd1 expression promotes BAF genomic mistargeting. Importantly, Brwd1 renormalization rescued aberrant BAF localization, along with associated changes in chromatin accessibility and gene expression. These findings establish BRWD1 as a key epigenomic mediator of normal neurodevelopment and an important contributor to DS-related phenotypes.

Identification of shared and differentiating genetic architecture for autism spectrum disorder, attention-deficit hyperactivity disorder and case subgroups

Attention-deficit hyperactivity disorder (ADHD) and autism spectrum disorder (ASD) are highly heritable neurodevelopmental conditions, with considerable overlap in their genetic etiology. We dissected their shared and distinct genetic etiology by cross-disorder analyses of large datasets. We identified seven loci shared by the disorders and five loci differentiating them. All five differentiating loci showed opposite allelic directions in the two disorders and significant associations with other traits, including educational attainment, neuroticism and regional brain volume. Integration with brain transcriptome data enabled us to identify and prioritize several significantly associated genes. The shared genomic fraction contributing to both disorders was strongly correlated with other psychiatric phenotypes, whereas the differentiating portion was correlated most strongly with cognitive traits. Additional analyses revealed that individuals diagnosed with both ASD and ADHD were double-loaded with genetic predispositions for both disorders and showed distinctive patterns of genetic association with other traits compared with the ASD-only and ADHD-only subgroups. These results provide insights into the biological foundation of the development of one or both conditions and of the factors driving psychopathology discriminatively toward either ADHD or ASD.

Population-level variation in enhancer expression identifies disease mechanisms in the human brain

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.

Impact of schizophrenia GWAS loci converge onto distinct pathways in cortical interneurons vs glutamatergic neurons during development

Remarkable advances have been made in schizophrenia (SCZ) GWAS, but gleaning biological insight from these loci is challenging. Genetic influences on gene expression (e.g., eQTLs) are cell type-specific, but most studies that attempt to clarify GWAS loci's influence on gene expression have employed tissues with mixed cell compositions that can obscure cell-specific effects. Furthermore, enriched SCZ heritability in the fetal brain underscores the need to study the impact of SCZ risk loci in specific developing neurons. MGE-derived cortical interneurons (cINs) are consistently affected in SCZ brains and show enriched SCZ heritability in human fetal brains. We identified SCZ GWAS risk genes that are dysregulated in iPSC-derived homogeneous populations of developing SCZ cINs. These SCZ GWAS loci differential expression (DE) genes converge on the PKC pathway. Their disruption results in PKC hyperactivity in developing cINs, leading to arborization deficits. We show that the fine-mapped GWAS locus in the ATP2A2 gene of the PKC pathway harbors enhancer marks by ATACseq and ChIPseq, and regulates ATP2A2 expression. We also generated developing glutamatergic neurons (GNs), another population with enriched SCZ heritability, and confirmed their functionality after transplantation into the mouse brain. Then, we identified SCZ GWAS risk genes that are dysregulated in developing SCZ GNs. GN-specific SCZ GWAS loci DE genes converge on the ion transporter pathway, distinct from those for cINs. Disruption of the pathway gene CACNA1D resulted in deficits of Ca2+ currents in developing GNs, suggesting compromised neuronal function by GWAS loci pathway deficits during development. This study allows us to identify cell type-specific and developmental stage-specific mechanisms of SCZ risk gene function, and may aid in identifying mechanism-based novel therapeutic targets.

The three-dimensional landscape of cortical chromatin accessibility in Alzheimer's disease

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.

The Neuroepigenome: Implications of Chemical and Physical Modifications of Genomic DNA in Schizophrenia

Schizophrenia is a chronic mental illness with a substantial genetic component. To unfold the complex etiology of schizophrenia, it is important to understand the interplay between genetic and nongenetic factors. Genetic factors involve variation in the DNA sequences of protein-coding genes, which directly contribute to phenotypic traits, and variation in noncoding sequences, which comprise 98% of the genome and contain DNA elements known to play a role in regulating gene expression. The epigenome refers to the chemical modifications on both DNA and the structural proteins that package DNA into the nucleus, which together regulate gene expression in specific cell types, conditions, and developmental stages. The dynamic nature of the epigenome makes it an ideal tool to investigate the relationship between inherited genetic mutations associated with schizophrenia and altered gene regulation throughout the course of brain development. In this review, we focus on the current understanding of the role of epigenetic marks and their three-dimensional nuclear organization in the developmental trajectory of distinct brain cell types to decipher the complex gene regulatory mechanisms that are disrupted in schizophrenia.

Penetrance and Pleiotropy of Polygenic Risk Scores for Schizophrenia, Bipolar Disorder, and Depression Among Adults in the US Veterans Affairs Health Care System

Importance

Serious mental illnesses, including schizophrenia, bipolar disorder, and depression, are heritable, highly multifactorial disorders and major causes of disability worldwide.

Objective

To benchmark the penetrance of current neuropsychiatric polygenic risk scores (PRSs) in the Veterans Health Administration health care system and to explore associations between PRS and broad categories of human disease via phenome-wide association studies.

Design, Setting, and Participants

Extensive Veterans Health Administration's electronic health records were assessed from October 1999 to January 2021, and an embedded cohort of 9378 individuals with confirmed diagnoses of schizophrenia or bipolar 1 disorder were found. The performance of schizophrenia, bipolar disorder, and major depression PRSs were compared in participants of African or European ancestry in the Million Veteran Program (approximately 400 000 individuals), and associations between PRSs and 1650 disease categories based on ICD-9/10 billing codes were explored. Last, genomic structural equation modeling was applied to derive novel PRSs indexing common and disorder-specific genetic factors. Analysis took place from January 2021 to January 2022.

Main Outcomes and Measures

Diagnoses based on in-person structured clinical interviews were compared with ICD-9/10 billing codes. PRSs were constructed using summary statistics from genome-wide association studies of schizophrenia, bipolar disorder, and major depression.

Results

Of 707 299 enrolled study participants, 459 667 were genotyped at the time of writing; 84 806 were of broadly African ancestry (mean [SD] age, 58 [12.1] years) and 314 909 were of broadly European ancestry (mean [SD] age, 66.4 [13.5] years). Among 9378 individuals with confirmed diagnoses of schizophrenia or bipolar 1 disorder, 8962 (95.6%) were correctly identified using ICD-9/10 codes (2 or more). Among those of European ancestry, PRSs were robustly associated with having received a diagnosis of schizophrenia (odds ratio [OR], 1.81 [95% CI, 1.76-1.87]; P < 10-257) or bipolar disorder (OR, 1.42 [95% CI, 1.39-1.44]; P < 10-295). Corresponding effect sizes in participants of African ancestry were considerably smaller for schizophrenia (OR, 1.35 [95% CI, 1.29-1.42]; P < 10-38) and bipolar disorder (OR, 1.16 [95% CI, 1.11-1.12]; P < 10-10). Neuropsychiatric PRSs were associated with increased risk for a range of psychiatric and physical health problems.

Conclusions and Relevance

Using diagnoses confirmed by in-person structured clinical interviews and current neuropsychiatric PRSs, the validity of an electronic health records-based phenotyping approach in US veterans was demonstrated, highlighting the potential of PRSs for disentangling biological and mediated pleiotropy.

A translational genomics approach identifies IL10RB as the top candidate gene target for COVID-19 susceptibility

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.

Genetics of the human microglia regulome refines Alzheimer's disease risk loci

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.

Transcriptome-wide association study in UK Biobank Europeans identifies associations with blood cell traits

Previous genome-wide association studies (GWAS) of hematological traits have identified over 10 000 distinct trait-specific risk loci. However, at these loci, the underlying causal mechanisms remain incompletely characterized. To elucidate novel biology and better understand causal mechanisms at known loci, we performed a transcriptome-wide association study (TWAS) of 29 hematological traits in 399 835 UK Biobank (UKB) participants of European ancestry using gene expression prediction models trained from whole blood RNA-seq data in 922 individuals. We discovered 557 gene-trait associations for hematological traits distinct from previously reported GWAS variants in European populations. Among the 557 associations, 301 were available for replication in a cohort of 141 286 participants of European ancestry from the Million Veteran Program. Of these 301 associations, 108 replicated at a strict Bonferroni adjusted threshold ($\alpha$= 0.05/301). Using our TWAS results, we systematically assigned 4261 out of 16 900 previously identified hematological trait GWAS variants to putative target genes. Compared to coloc, our TWAS results show reduced specificity and increased sensitivity in external datasets to assign variants to target genes.

Engagement of vascular early response genes typifies mild cognitive impairment

INTRODUCTION

Molecular responses in the brains of persons with mild cognitive impairment (MCI), the earliest transitional state between normal aging and early Alzheimer's disease (AD), are poorly understood.

METHODS

We examined AD-related neuropathology and transcriptome changes in the neocortex of individuals with MCI relative to controls and temporal responses to the mild hypoxia in mouse brains.

RESULTS

Subsets of vascular early response to hypoxia genes were upregulated in MCI prior to the buildup of AD neuropathology. Early activation of pro-angiogenic hypoxia-inducible factor signaling in response to mild hypoxia was detected in mouse brains similar to those that were altered in MCI. Protracted responses to hypoxia were characterized by activation of phosphoinositide 3-kinase (PI3K)-protein kinase B (Akt)-the mammalian target of rapamycin (mTOR) pathways in brain microvessel isolates.

DISCUSSION

These findings suggest that cerebrovascular remodeling is an important antecedent to the development of dementia and a component of the homeostatic response to reduced oxygen tension in aging prior to the onset of AD.

Altered gene expression and PTSD symptom dimensions in World Trade Center responders

Despite experiencing a significant trauma, only a subset of World Trade Center (WTC) rescue and recovery workers developed posttraumatic stress disorder (PTSD). Identification of biomarkers is critical to the development of targeted interventions for treating disaster responders and potentially preventing the development of PTSD in this population. Analysis of gene expression from these individuals can help in identifying biomarkers of PTSD. We established a well-phenotyped sample of 371 WTC responders, recruited from a longitudinal WTC responder cohort using stratified random sampling, by obtaining blood, self-reported and clinical interview data. Using bulk RNA-sequencing from whole blood, we examined the association between gene expression and WTC-related PTSD symptom severity on (i) highest lifetime Clinician-Administered PTSD Scale (CAPS) score, (ii) past-month CAPS score, and (iii) PTSD symptom dimensions using a 5-factor model of re-experiencing, avoidance, emotional numbing, dysphoric arousal and anxious arousal symptoms. We corrected for sex, age, genotype-derived principal components and surrogate variables. Finally, we performed a meta-analysis with existing PTSD studies (total N = 1016), using case/control status as the predictor and correcting for these variables. We identified 66 genes significantly associated with total highest lifetime CAPS score (FDR-corrected p < 0.05), and 31 genes associated with total past-month CAPS score. Our more granular analyses of PTSD symptom dimensions identified additional genes that did not reach statistical significance in our analyses with total CAPS scores. In particular, we identified 82 genes significantly associated with lifetime anxious arousal symptoms. Several genes significantly associated with multiple PTSD symptom dimensions and total lifetime CAPS score (SERPINA1, RPS6KA1, and STAT3) have been previously associated with PTSD. Geneset enrichment of these findings has identified pathways significant in metabolism, immune signaling, other psychiatric disorders, neurological signaling, and cellular structure. Our meta-analysis revealed 10 genes that reached genome-wide significance, all of which were downregulated in cases compared to controls (CIRBP, TMSB10, FCGRT, CLIC1, RPS6KB2, HNRNPUL1, ALDOA, NACA, ZNF429 and COPE). Additionally, cellular deconvolution highlighted an enrichment in CD4 T cells and eosinophils in responders with PTSD compared to controls. The distinction in significant genes between total lifetime CAPS score and the anxious arousal symptom dimension of PTSD highlights a potential biological difference in the mechanism underlying the heterogeneity of the PTSD phenotype. Future studies should be clear about methods used to analyze PTSD status, as phenotypes based on PTSD symptom dimensions may yield different gene sets than combined CAPS score analysis. Potential biomarkers implicated from our meta-analysis may help improve therapeutic target development for PTSD.