Cell-type-specific cis-eQTLs in eight human brain cell types identify novel risk genes for psychiatric and neurological disorders

Revealing brain cell-stratified causality through dissecting causal variants according to their cell-type-specific effects on gene expression

The human brain has been implicated in the pathogenesis of several complex diseases. Taking advantage of single-cell techniques, genome-wide association studies (GWAS) have taken it a step further and revealed brain cell-type-specific functions for disease loci. However, genetic causal associations inferred by Mendelian randomization (MR) studies usually include all instrumental variables from GWAS, which hampers the understanding of cell-specific causality. Here, we developed an analytical framework, Cell-Stratified MR (csMR), to investigate cell-stratified causality through colocalizing GWAS signals with single-cell eQTL from different brain cells. By applying to obesity-related traits, our results demonstrate the cell-type-specific effects of GWAS variants on gene expression, and indicate the benefits of csMR to identify cell-type-specific causal effect that is often hidden from bulk analyses. We also found csMR valuable to reveal distinct causal pathways between different obesity indicators. These findings suggest the value of our approach to prioritize target cells for extending genetic causation studies.

Bayesian colocalization of GWAS and eQTL signals reveals cell type-specific genes and regulatory variants for susceptibility to subtypes of ischemic stroke

A colocalization analysis of genome-wide association study (GWAS) signals and expression quantitative trait loci (eQTL) was conducted to pinpoint target genes and their regulatory nucleotide variants for subtypes of ischemic stroke. We utilized GWAS data from prominent meta-analysis consortia (MEGASTROKE and GIGASTROKE) and single-cell eQTL data in brain and blood tissues to enhance accuracy and minimize noise inherent in bulk RNA-seq. Employing Bayesian colocalization methods, we identified ten shared loci between GWAS and eQTL signals, targeting five eGenes. Specifically, RAPH1 and ICA1L were discovered for small vessel stroke (SVS), whereas SCYL3, CAV1, and CAV2 were for cardioembolic stroke (CS). However, no findings have been made for large artery stroke. The exploration and subsequent functional analysis of causal variants within the colocalized regions revealed their regulatory roles, particularly as enhancer variants (e.g., rs144505847 and rs72932755 targeting ICA1L; rs629234 targeting SCYL3; rs3807989 targeting CAV1 and CAV2). Notably, our study unveiled that all eQTL for CS were identified in oligodendrocytes, while those for SVS were across excitatory neurons, astrocytes, and oligodendrocyte precursor cells. This underscores the heterogeneous tissue-specific genetic factors by subtypes of ischemic stroke. The study emphasizes the need for intensive research efforts to discover causative genes and variants, unravelling the cell type-specific genetic architecture of ischemic stroke subtypes. This knowledge is crucial for advancing our understanding of the underlying pathophysiology and paving the way for precision neurology applications.

Omics approaches to investigate the pathogenesis of suicide

Suicide is the second leading cause of death in U.S. adolescents and young adults, and generally associated with a psychiatric disorder. Suicidal behavior has a complex etiology and pathogenesis. Moderate heritability suggests genetic causes. Associations between childhood and recent life adversity indicate contributions from epigenetic factors. Genomic contributions to suicide pathogenesis remain largely unknown. This paper is based on a workshop held to design strategies to identify molecular drivers of suicide neurobiology that would be putative new treatment targets. The panel determined that, while bulk tissue studies provide comprehensive information, single-nucleus approaches identifying cell-type specific changes are needed. While single nuclei techniques lack information on cytoplasm, processes, spines, and synapses, spatial multiomic technologies on intact tissue detect cell alterations specific to brain tissue layers and subregions. Because suicide has genetic and environmental drivers, multiomic approaches combining cell-type specific epigenome, transcriptome, and proteome provide a more complete picture of pathogenesis. To determine the direction of effect of suicide risk gene variants on RNA and protein expression, and how these interact with epigenetic marks, single nuclei and spatial multiomics quantitative trait loci maps should be integrated with whole genome sequencing and genome-wide association databases. The workshop concluded with the recommendation for the formation of an international suicide biology consortium that will bring together brain banks and investigators with expertise in cutting-edge omics technologies to delineate the biology of suicide and identify novel potential treatment targets to be tested in cellular and animal models for drug and biomarkers discovery, to guide suicide prevention.

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, 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.

Systems biology dissection of PTSD and MDD across brain regions, cell types, and blood

The molecular pathology of stress-related disorders remains elusive. Our brain multiregion, multiomic study of posttraumatic stress disorder (PTSD) and major depressive disorder (MDD) included the central nucleus of the amygdala, hippocampal dentate gyrus, and medial prefrontal cortex (mPFC). Genes and exons within the mPFC carried most disease signals replicated across two independent cohorts. Pathways pointed to immune function, neuronal and synaptic regulation, and stress hormones. Multiomic factor and gene network analyses provided the underlying genomic structure. Single nucleus RNA sequencing in dorsolateral PFC revealed dysregulated (stress-related) signals in neuronal and non-neuronal cell types. Analyses of brain-blood intersections in >50,000 UK Biobank participants were conducted along with fine-mapping of the results of PTSD and MDD genome-wide association studies to distinguish risk from disease processes. Our data suggest shared and distinct molecular pathology in both disorders and propose potential therapeutic targets and biomarkers.

Brain cell-type shifts in Alzheimer's disease, autism, and schizophrenia interrogated using methylomics and genetics

Few neuropsychiatric disorders have replicable biomarkers, prompting high-resolution and large-scale molecular studies. However, we still lack consensus on a more foundational question: whether quantitative shifts in cell types-the functional unit of life-contribute to neuropsychiatric disorders. Leveraging advances in human brain single-cell methylomics, we deconvolve seven major cell types using bulk DNA methylation profiling across 1270 postmortem brains, including from individuals diagnosed with Alzheimer's disease, schizophrenia, and autism. We observe and replicate cell-type compositional shifts for Alzheimer's disease (endothelial cell loss), autism (increased microglia), and schizophrenia (decreased oligodendrocytes), and find age- and sex-related changes. Multiple layers of evidence indicate that endothelial cell loss contributes to Alzheimer's disease, with comparable effect size to APOE genotype among older people. Genome-wide association identified five genetic loci related to cell-type composition, involving plausible genes for the neurovascular unit (P2RX5 and TRPV3) and excitatory neurons (DPY30 and MEMO1). These results implicate specific cell-type shifts in the pathophysiology of neuropsychiatric disorders.

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

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.

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

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.

Integrating single cell expression quantitative trait loci summary statistics to understand complex trait risk genes

Transcriptome-wide association study (TWAS) is a popular approach to dissect the functional consequence of disease associated non-coding variants. Most existing TWAS use bulk tissues and may not have the resolution to reveal cell-type specific target genes. Single-cell expression quantitative trait loci (sc-eQTL) datasets are emerging. The largest bulk- and sc-eQTL datasets are most conveniently available as summary statistics, but have not been broadly utilized in TWAS. Here, we present a new method EXPRESSO (EXpression PREdiction with Summary Statistics Only), to analyze sc-eQTL summary statistics, which also integrates 3D genomic data and epigenomic annotation to prioritize causal variants. EXPRESSO substantially improves existing methods. We apply EXPRESSO to analyze multi-ancestry GWAS datasets for 14 autoimmune diseases. EXPRESSO uniquely identifies 958 novel gene x trait associations, which is 26% more than the second-best method. Among them, 492 are unique to cell type level analysis and missed by TWAS using whole blood. We also develop a cell type aware drug repurposing pipeline, which leverages EXPRESSO results to identify drug compounds that can reverse disease gene expressions in relevant cell types. Our results point to multiple drugs with therapeutic potentials, including metformin for type 1 diabetes, and vitamin K for ulcerative colitis.

Asian-European differentiation of schizophrenia-associated genes driven by admixture and natural selection

The European-centered genome-wide association studies of schizophrenia (SCZ) may not be well applied to non-European populations. We analyzed 1,592 reported SCZ-associated genes using the public genome data and found an overall higher Asian-European differentiation on the SCZ-associated variants than at the genome-wide level. Notable examples included 15 missense variants, a regulatory variant SLC5A10-rs1624825, and a damaging variant TSPAN18-rs1001292. Independent local adaptations in recent 25,000 years, after the Asian-European divergence, could have contributed to such genetic differentiation, as were identified at a missense mutation LTN1-rs57646126-A in Asians, and a non-risk allele ZSWIM6-rs72761442-G in Europeans. Altai-Neanderthal-derived alleles may have opposite effects on SCZ susceptibility between ancestries. Furthermore, adaptive introgression was detected on the non-risk haplotype at 1q21.2 in Europeans, while in Asians it was observed on the SCZ risk haplotype at 3p21.31 which is also potentially ultra-violet protective. This study emphasizes the importance of including more representative Asian samples in future SCZ studies.

Efficient and accurate mixed model association tool for single-cell eQTL analysis

Understanding the genetic basis of gene expression can help us understand the molecular underpinnings of human traits and disease. Expression quantitative trait locus (eQTL) mapping can help in studying this relationship but have been shown to be very cell-type specific, motivating the use of single-cell RNA sequencing and single-cell eQTLs to obtain a more granular view of genetic regulation. Current methods for single-cell eQTL mapping either rely on the "pseudobulk" approach and traditional pipelines for bulk transcriptomics or do not scale well to large datasets. Here, we propose SAIGE-QTL, a robust and scalable tool that can directly map eQTLs using single-cell profiles without needing aggregation at the pseudobulk level. Additionally, SAIGE-QTL allows for testing the effects of less frequent/rare genetic variation through set-based tests, which is traditionally excluded from eQTL mapping studies. We evaluate the performance of SAIGE-QTL on both real and simulated data and demonstrate the improved power for eQTL mapping over existing pipelines.

Altered Expression of PDE4 Genes in Schizophrenia: Insights from a Brain and Blood Sample Meta-Analysis and iPSC-Derived Neurons

Schizophrenia symptomatology includes negative symptoms and cognitive impairment. Several studies have linked schizophrenia with the PDE4 family of enzymes due to their genetic association and function in cognitive processes such as long-term potentiation. We conducted a systematic gene expression meta-analysis of four PDE4 genes (PDE4A-D) in 10 brain sample datasets (437 samples) and three blood sample datasets (300 samples). Subsequently, we measured mRNA levels in iPSC-derived hippocampal dentate gyrus neurons generated from fibroblasts of three groups: healthy controls, healthy monozygotic twins (MZ), and their MZ siblings with schizophrenia. We found downregulation of PDE4B in brain tissues, further validated by independent data of the CommonMind consortium (515 samples). Interestingly, the downregulation signal was present in a subgroup of the patients, while the others showed no differential expression or even upregulation. Notably, PDE4A, PDE4B, and PDE4D exhibited upregulation in iPSC-derived neurons compared to healthy controls, whereas in blood samples, PDE4B was found to be upregulated while PDE4A was downregulated. While the precise mechanism and direction of altered PDE4 expression necessitate further investigation, the observed multilevel differential expression across the brain, blood, and iPSC-derived neurons compellingly suggests the involvement of PDE4 genes in the pathophysiology of schizophrenia.

From genetic associations to genes: methods, applications, and challenges

Genome-wide association studies (GWASs) have identified numerous genetic loci associated with human traits and diseases. However, pinpointing the causal genes remains a challenge, which impedes the translation of GWAS findings into biological insights and medical applications. In this review, we provide an in-depth overview of the methods and technologies used for prioritizing genes from GWAS loci, including gene-based association tests, integrative analysis of GWAS and molecular quantitative trait loci (xQTL) data, linking GWAS variants to target genes through enhancer-gene connection maps, and network-based prioritization. We also outline strategies for generating context-dependent xQTL data and their applications in gene prioritization. We further highlight the potential of gene prioritization in drug repurposing. Lastly, we discuss future challenges and opportunities in this field.

Genome-wide association analyses identify 95 risk loci and provide insights into the neurobiology of post-traumatic stress disorder

Post-traumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 new). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (for example, GRIA1, GRM8 and CACNA1E), developmental, axon guidance and transcription factors (for example, FOXP2, EFNA5 and DCC), synaptic structure and function genes (for example, PCLO, NCAM1 and PDE4B) and endocrine or immune regulators (for example, ESR1, TRAF3 and TANK). Additional top genes influence stress, immune, fear and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Single-nucleus and spatial transcriptomic profiling of human temporal cortex and white matter reveals novel associations with AD pathology

Alzheimer's disease (AD) is a neurodegenerative disorder with complex pathological manifestations and is the leading cause of cognitive decline and dementia in elderly individuals. A major goal in AD research is to identify new therapeutic pathways by studying the molecular and cellular changes in the disease, either downstream or upstream of the pathological hallmarks. In this study, we present a comprehensive investigation of cellular heterogeneity from the temporal cortex region of 40 individuals, comprising healthy donors and individuals with differing tau and amyloid burden. Using single-nucleus transcriptome analysis of 430,271 nuclei from both gray and white matter of these individuals, we identified cell type-specific subclusters in both neuronal and glial cell types with varying degrees of association with AD pathology. In particular, these associations are present in layer specific glutamatergic (excitatory) neuronal types, along with GABAergic (inhibitory) neurons and glial subtypes. These associations were observed in early as well as late pathological progression. We extended this analysis by performing multiplexed in situ hybridization using the CARTANA platform, capturing 155 genes in 13 individuals with varying levels of tau pathology. By modeling the spatial distribution of these genes and their associations with the pathology, we not only replicated key findings from our snRNA data analysis, but also identified a set of cell type-specific genes that show selective enrichment or depletion near pathological inclusions. Together, our findings allow us to prioritize specific cell types and pathways for targeted interventions at various stages of pathological progression in AD.

A novel multi-omics mendelian randomization method for gene set enrichment and its application to psychiatric disorders

Genome-wide association studies (GWAS) of psychiatric disorders (PD) yield numerous loci with significant signals, but often do not implicate specific genes. Because GWAS risk loci are enriched in expression/protein/methylation quantitative loci (e/p/mQTL, hereafter xQTL), transcriptome/proteome/methylome-wide association studies (T/P/MWAS, hereafter XWAS) that integrate xQTL and GWAS information, can link GWAS signals to effects on specific genes. To further increase detection power, gene signals are aggregated within relevant gene sets (GS) by performing gene set enrichment (GSE) analyses. Often GSE methods test for enrichment of "signal" genes in curated GS while overlooking their linkage disequilibrium (LD) structure, allowing for the possibility of increased false positive rates. Moreover, no GSE tool uses xQTL information to perform mendelian randomization (MR) analysis. To make causal inference on association between PD and GS, we develop a novel MR GSE (MR-GSE) procedure. First, we generate a "synthetic" GWAS for each MSigDB GS by aggregating summary statistics for x-level (mRNA, protein or DNA methylation (DNAm) levels) from the largest xQTL studies available) of genes in a GS. Second, we use synthetic GS GWAS as exposure in a generalized summary-data-based-MR analysis of complex trait outcomes. We applied MR-GSE to GWAS of nine important PD. When applied to the underpowered opioid use disorder GWAS, none of the four analyses yielded any signals, which suggests a good control of false positive rates. For other PD, MR-GSE greatly increased the detection of GO terms signals (2,594) when compared to the commonly used (non-MR) GSE method (286). Some of the findings might be easier to adapt for treatment, e.g., our analyses suggest modest positive effects for supplementation with certain vitamins and/or omega-3 for schizophrenia, bipolar and major depression disorder patients. Similar to other MR methods, when applying MR-GSE researchers should be mindful of the confounding effects of horizontal pleiotropy on statistical inference.

Investigating the role of common cis-regulatory variants in modifying penetrance of putatively damaging, inherited variants in severe neurodevelopmental disorders

Recent work has revealed an important role for rare, incompletely penetrant inherited coding variants in neurodevelopmental disorders (NDDs). Additionally, we have previously shown that common variants contribute to risk for rare NDDs. Here, we investigate whether common variants exert their effects by modifying gene expression, using multi-cis-expression quantitative trait loci (cis-eQTL) prediction models. We first performed a transcriptome-wide association study for NDDs using 6987 probands from the Deciphering Developmental Disorders (DDD) study and 9720 controls, and found one gene, RAB2A, that passed multiple testing correction (p = 6.7 × 10-7). We then investigated whether cis-eQTLs modify the penetrance of putatively damaging, rare coding variants inherited by NDD probands from their unaffected parents in a set of 1700 trios. We found no evidence that unaffected parents transmitting putatively damaging coding variants had higher genetically-predicted expression of the variant-harboring gene than their child. In probands carrying putatively damaging variants in constrained genes, the genetically-predicted expression of these genes in blood was lower than in controls (p = 2.7 × 10-3). However, results for proband-control comparisons were inconsistent across different sets of genes, variant filters and tissues. We find limited evidence that common cis-eQTLs modify penetrance of rare coding variants in a large cohort of NDD probands.

Multimodal Nature of the Single-cell Primate Brain Atlas: Morphology, Transcriptome, Electrophysiology, and Connectivity

Primates exhibit complex brain structures that augment cognitive function. The neocortex fulfills high-cognitive functions through billions of connected neurons. These neurons have distinct transcriptomic, morphological, and electrophysiological properties, and their connectivity principles vary. These features endow the primate brain atlas with a multimodal nature. The recent integration of next-generation sequencing with modified patch-clamp techniques is revolutionizing the way to census the primate neocortex, enabling a multimodal neuronal atlas to be established in great detail: (1) single-cell/single-nucleus RNA-seq technology establishes high-throughput transcriptomic references, covering all major transcriptomic cell types; (2) patch-seq links the morphological and electrophysiological features to the transcriptomic reference; (3) multicell patch-clamp delineates the principles of local connectivity. Here, we review the applications of these technologies in the primate neocortex and discuss the current advances and tentative gaps for a comprehensive understanding of the primate neocortex.

Cell-type-specific and disease-associated expression quantitative trait loci in the human lung

Common genetic variants confer substantial risk for chronic lung diseases, including pulmonary fibrosis. Defining the genetic control of gene expression in a cell-type-specific and context-dependent manner is critical for understanding the mechanisms through which genetic variation influences complex traits and disease pathobiology. To this end, we performed single-cell RNA sequencing of lung tissue from 66 individuals with pulmonary fibrosis and 48 unaffected donors. Using a pseudobulk approach, we mapped expression quantitative trait loci (eQTLs) across 38 cell types, observing both shared and cell-type-specific regulatory effects. Furthermore, we identified disease interaction eQTLs and demonstrated that this class of associations is more likely to be cell-type-specific and linked to cellular dysregulation in pulmonary fibrosis. Finally, we connected lung disease risk variants to their regulatory targets in disease-relevant cell types. These results indicate that cellular context determines the impact of genetic variation on gene expression and implicates context-specific eQTLs as key regulators of lung homeostasis and disease.

Cell biology of Parkinson's disease: Mechanisms of synaptic, lysosomal, and mitochondrial dysfunction

Parkinson's disease (PD) is a growing cause of disability worldwide and there is a critical need for the development of disease-modifying therapies to slow or stop disease progression. Recent advances in characterizing the genetics of PD have expanded our understanding of the cell biology of this disorder. Mitochondrial oxidative stress, defects in synaptic function, and impaired lysosomal activity have been shown to be linked in PD, resulting in a pathogenic feedback cycle involving the accumulation of toxic oxidized dopamine and alpha-synuclein. In this review, we will highlight recent data on a subset of PD-linked genes which have key roles in these pathways and the pathogenic cycle. We will furthermore discuss findings highlighting the importance of dynamic mitochondria-lysosome contact sites that mediate direct inter-organelle cross-talk in the pathogenesis of PD and related disorders.

Single-cell multiplex chromatin and RNA interactions in ageing human brain

Dynamically organized chromatin complexes often involve multiplex chromatin interactions and sometimes chromatin-associated RNA1-3. Chromatin complex compositions change during cellular differentiation and ageing, and are expected to be highly heterogeneous among terminally differentiated single cells4-7. Here we introduce the multinucleic acid interaction mapping in single cells (MUSIC) technique for concurrent profiling of multiplex chromatin interactions, gene expression and RNA-chromatin associations within individual nuclei. When applied to 14 human frontal cortex samples from older donors, MUSIC delineated diverse cortical cell types and states. We observed that nuclei exhibiting fewer short-range chromatin interactions were correlated with both an 'older' transcriptomic signature and Alzheimer's disease pathology. Furthermore, the cell type exhibiting chromatin contacts between cis expression quantitative trait loci and a promoter tends to be that in which these cis expression quantitative trait loci specifically affect the expression of their target gene. In addition, female cortical cells exhibit highly heterogeneous interactions between XIST non-coding RNA and chromosome X, along with diverse spatial organizations of the X chromosomes. MUSIC presents a potent tool for exploration of chromatin architecture and transcription at cellular resolution in complex tissues.

Cell subtype-specific effects of genetic variation in the Alzheimer's disease brain

The relationship between genetic variation and gene expression in brain cell types and subtypes remains understudied. Here, we generated single-nucleus RNA sequencing data from the neocortex of 424 individuals of advanced age; we assessed the effect of genetic variants on RNA expression in cis (cis-expression quantitative trait loci) for seven cell types and 64 cell subtypes using 1.5 million transcriptomes. This effort identified 10,004 eGenes at the cell type level and 8,099 eGenes at the cell subtype level. Many eGenes are only detected within cell subtypes. A new variant influences APOE expression only in microglia and is associated with greater cerebral amyloid angiopathy but not Alzheimer's disease pathology, after adjusting for APOEε4, providing mechanistic insights into both pathologies. Furthermore, only a TMEM106B variant affects the proportion of cell subtypes. Integration of these results with genome-wide association studies highlighted the targeted cell type and probable causal gene within Alzheimer's disease, schizophrenia, educational attainment and Parkinson's disease loci.

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.

SingleQ: a comprehensive database of single-cell expression quantitative trait loci (sc-eQTLs) cross human tissues

Mapping of expression quantitative trait loci (eQTLs) and other molecular QTLs can help characterize the modes of action of disease-associated genetic variants. However, current eQTL databases present data from bulk RNA-seq approaches, which cannot shed light on the cell type- and environment-specific regulation of disease-associated genetic variants. Here, we introduce our Single-cell eQTL Interactive Database which collects single-cell eQTL (sc-eQTL) datasets and provides online visualization of sc-eQTLs across different cell types in a user-friendly manner. Although sc-eQTL mapping is still in its early stage, our database curates the most comprehensive summary statistics of sc-eQTLs published to date. sc-eQTL studies have revolutionized our understanding of gene regulation in specific cellular contexts, and we anticipate that our database will further accelerate the research of functional genomics. Database URL: http://www.sqraolab.com/scqtl.

Deep sequencing of proteotoxicity modifier genes uncovers a Presenilin-2/beta-amyloid-actin genetic risk module shared among alpha-synucleinopathies

Whether neurodegenerative diseases linked to misfolding of the same protein share genetic risk drivers or whether different protein-aggregation pathologies in neurodegeneration are mechanistically related remains uncertain. Conventional genetic analyses are underpowered to address these questions. Through careful selection of patients based on protein aggregation phenotype (rather than clinical diagnosis) we can increase statistical power to detect associated variants in a targeted set of genes that modify proteotoxicities. Genetic modifiers of alpha-synuclein (ɑS) and beta-amyloid (Aβ) cytotoxicity in yeast are enriched in risk factors for Parkinson's disease (PD) and Alzheimer's disease (AD), respectively. Here, along with known AD/PD risk genes, we deeply sequenced exomes of 430 ɑS/Aβ modifier genes in patients across alpha-synucleinopathies (PD, Lewy body dementia and multiple system atrophy). Beyond known PD genes GBA1 and LRRK2, rare variants AD genes (CD33, CR1 and PSEN2) and Aβ toxicity modifiers involved in RhoA/actin cytoskeleton regulation (ARGHEF1, ARHGEF28, MICAL3, PASK, PKN2, PSEN2) were shared risk factors across synucleinopathies. Actin pathology occurred in iPSC synucleinopathy models and RhoA downregulation exacerbated ɑS pathology. Even in sporadic PD, the expression of these genes was altered across CNS cell types. Genome-wide CRISPR screens revealed the essentiality of PSEN2 in both human cortical and dopaminergic neurons, and PSEN2 mutation carriers exhibited diffuse brainstem and cortical synucleinopathy independent of AD pathology. PSEN2 contributes to a common-risk signal in PD GWAS and regulates ɑS expression in neurons. Our results identify convergent mechanisms across synucleinopathies, some shared with AD.

Genetic contribution to the comorbidity between attention-deficit/hyperactivity disorder and substance use disorders

We characterized the genetic architecture of the attention-deficit hyperactivity disorder-substance use disorder (ADHD-SUD) relationship by investigating genetic correlation, causality, pleiotropy, and common polygenic risk. Summary statistics from genome-wide association studies (GWAS) were used to investigate ADHD (Neff = 51,568), cannabis use disorder (CanUD, Neff = 161,053), opioid use disorder (OUD, Neff = 57,120), problematic alcohol use (PAU, Neff = 502,272), and problematic tobacco use (PTU, Neff = 97,836). ADHD, CanUD, and OUD GWAS meta-analyses included cohorts with case definitions based on different diagnostic criteria. PAU GWAS combined information related to alcohol use disorder, alcohol dependence, and the items related to alcohol problematic consequences assessed by the alcohol use disorders identification test. PTU GWAS was generated a multi-trait analysis including information regarding Fagerström Test for Nicotine Dependence and cigarettes per day. Linkage disequilibrium score regression analyses indicated positive genetic correlation with CanUD, OUD, PAU, and PTU. Genomic structural equation modeling showed that these genetic correlations were related to two latent factors: one including ADHD, CanUD, and PTU and the other with OUD and PAU. The evidence of a causal effect of PAU and PTU on ADHD was stronger than the reverse in the two-sample Mendelian randomization analysis. Conversely, similar strength of evidence was found between ADHD and CanUD. CADM2 rs62250713 was a pleiotropic SNP between ADHD and all SUDs. We found seven, one, and twenty-eight pleiotropic variants between ADHD and CanUD, PAU, and PTU, respectively. Finally, OUD, CanUD, and PAU PRS were associated with increased odds of ADHD. Our findings demonstrated the contribution of multiple pleiotropic mechanisms to the comorbidity between ADHD and SUDs.

Machine learning nominates the inositol pathway and novel genes in Parkinson's disease

There are 78 loci associated with Parkinson's disease in the most recent genome-wide association study (GWAS), yet the specific genes driving these associations are mostly unknown. Herein, we aimed to nominate the top candidate gene from each Parkinson's disease locus and identify variants and pathways potentially involved in Parkinson's disease. We trained a machine learning model to predict Parkinson's disease-associated genes from GWAS loci using genomic, transcriptomic and epigenomic data from brain tissues and dopaminergic neurons. We nominated candidate genes in each locus and identified novel pathways potentially involved in Parkinson's disease, such as the inositol phosphate biosynthetic pathway (INPP5F, IP6K2, ITPKB and PPIP5K2). Specific common coding variants in SPNS1 and MLX may be involved in Parkinson's disease, and burden tests of rare variants further support that CNIP3, LSM7, NUCKS1 and the polyol/inositol phosphate biosynthetic pathway are associated with the disease. Functional studies are needed to further analyse the involvements of these genes and pathways in Parkinson's disease.

Genomic insights into the comorbidity between type 2 diabetes and schizophrenia

Multimorbidity represents an increasingly important public health challenge with far-reaching implications for health management and policy. Mental health and metabolic diseases have a well-established epidemiological association. In this study, we investigate the genetic intersection between type 2 diabetes and schizophrenia. We use Mendelian randomization to examine potential causal relationships between the two conditions and related endophenotypes. We report no compelling evidence that type 2 diabetes genetic liability potentially causally influences schizophrenia risk and vice versa. Our findings show that increased body mass index (BMI) has a protective effect against schizophrenia, in contrast to the well-known risk-increasing effect of BMI on type 2 diabetes risk. We identify evidence of colocalization of association signals for these two conditions at 11 genomic loci, six of which have opposing directions of effect for type 2 diabetes and schizophrenia. To elucidate these colocalizing signals, we integrate multi-omics data from bulk and single-cell gene expression studies, along with functional information. We identify putative effector genes and find that they are enriched for homeostasis and lipid-related pathways. We also highlight drug repurposing opportunities including N-methyl-D-aspartate (NMDA) receptor antagonists. Our findings provide insights into shared biological mechanisms for type 2 diabetes and schizophrenia, highlighting common factors that influence the risk of the two conditions in opposite directions and shedding light on the complex nature of this comorbidity.

Brain tissue- and cell type-specific eQTL Mendelian randomization reveals efficacy of FADS1 and FADS2 on cognitive function

Epidemiological studies suggested an association between omega-3 fatty acids and cognitive function. However, the causal role of the fatty acid desaturase (FADS) gene, which play a key role in regulating omega-3 fatty acids biosynthesis, on cognitive function is unclear. Hence, we used two-sample Mendelian randomization (MR) to estimate the gene-specific causal effect of omega-3 fatty acids (N = 114,999) on cognitive function (N = 300,486). Tissue- and cell type-specific effects of FADS1/FADS2 expression on cognitive function were estimated using brain tissue cis-expression quantitative trait loci (cis-eQTL) datasets (GTEx, N ≤ 209; MetaBrain, N ≤ 8,613) and single cell cis-eQTL data (N = 373), respectively. These causal effects were further evaluated in whole blood cis-eQTL data (N ≤ 31,684). A series of sensitivity analyses were conducted to validate MR assumptions. Leave-one-out MR showed a FADS gene-specific effect of omega-3 fatty acids on cognitive function [β = -1.3 × 10-2, 95% confidence interval (CI) (-2.2 × 10-2, -5 × 10-3), P = 2 × 10-3]. Tissue-specific MR showed an effect of increased FADS1 expression in cerebellar hemisphere and FADS2 expression in nucleus accumbens basal ganglia on maintaining cognitive function, while decreased FADS1 expression in nine brain tissues on maintaining cognitive function [colocalization probability (PP.H4) ranged from 71.7% to 100.0%]. Cell type-specific MR showed decreased FADS1/FADS2 expression in oligodendrocyte was associated with maintaining cognitive function (PP.H4 = 82.3%, respectively). Increased FADS1/FADS2 expression in whole blood showed an effect on cognitive function maintenance (PP.H4 = 86.6% and 88.4%, respectively). This study revealed putative causal effect of FADS1/FADS2 expression in brain tissues and blood on cognitive function. These findings provided evidence to prioritize FADS gene as potential target gene for maintenance of cognitive function.

Cell-type deconvolution of bulk-blood RNA-seq reveals biological insights into neuropsychiatric disorders

Genome-wide association studies (GWASs) have uncovered susceptibility loci associated with psychiatric disorders such as bipolar disorder (BP) and schizophrenia (SCZ). However, most of these loci are in non-coding regions of the genome, and the causal mechanisms of the link between genetic variation and disease risk is unknown. Expression quantitative trait locus (eQTL) analysis of bulk tissue is a common approach used for deciphering underlying mechanisms, although this can obscure cell-type-specific signals and thus mask trait-relevant mechanisms. Although single-cell sequencing can be prohibitively expensive in large cohorts, computationally inferred cell-type proportions and cell-type gene expression estimates have the potential to overcome these problems and advance mechanistic studies. Using bulk RNA-seq from 1,730 samples derived from whole blood in a cohort ascertained from individuals with BP and SCZ, this study estimated cell-type proportions and their relation with disease status and medication. For each cell type, we found between 2,875 and 4,629 eGenes (genes with an associated eQTL), including 1,211 that are not found on the basis of bulk expression alone. We performed a colocalization test between cell-type eQTLs and various traits and identified hundreds of associations that occur between cell-type eQTLs and GWASs but that are not detected in bulk eQTLs. Finally, we investigated the effects of lithium use on the regulation of cell-type expression loci and found examples of genes that are differentially regulated according to lithium use. Our study suggests that applying computational methods to large bulk RNA-seq datasets of non-brain tissue can identify disease-relevant, cell-type-specific biology of psychiatric disorders and psychiatric medication.

PICALO: principal interaction component analysis for the identification of discrete technical, cell-type, and environmental factors that mediate eQTLs

Expression quantitative trait loci (eQTL) offer insights into the regulatory mechanisms of trait-associated variants, but their effects often rely on contexts that are unknown or unmeasured. We introduce PICALO, a method for hidden variable inference of eQTL contexts. PICALO identifies and disentangles technical from biological context in heterogeneous blood and brain bulk eQTL datasets. These contexts are biologically informative and reproducible, outperforming cell counts or expression-based principal components. Furthermore, we show that RNA quality and cell type proportions interact with thousands of eQTLs. Knowledge of hidden eQTL contexts may aid in the inference of functional mechanisms underlying disease variants.

Microglia regulation of central nervous system myelin health and regeneration

Microglia are resident macrophages of the central nervous system that have key functions in its development, homeostasis and response to damage and infection. Although microglia have been increasingly implicated in contributing to the pathology that underpins neurological dysfunction and disease, they also have crucial roles in neurological homeostasis and regeneration. This includes regulation of the maintenance and regeneration of myelin, the membrane that surrounds neuronal axons, which is required for axonal health and function. Myelin is damaged with normal ageing and in several neurodegenerative diseases, such as multiple sclerosis and Alzheimer disease. Given the lack of approved therapies targeting myelin maintenance or regeneration, it is imperative to understand the mechanisms by which microglia support and restore myelin health to identify potential therapeutic approaches. However, the mechanisms by which microglia regulate myelin loss or integrity are still being uncovered. In this Review, we discuss recent work that reveals the changes in white matter with ageing and neurodegenerative disease, how this relates to microglia dynamics during myelin damage and regeneration, and factors that influence the regenerative functions of microglia.

Disentangling Accelerated Cognitive Decline from the Normal Aging Process and Unraveling Its Genetic Components: A Neuroimaging-Based Deep Learning Approach

Background

The progressive cognitive decline, an integral component of Alzheimer's disease (AD), unfolds in tandem with the natural aging process. Neuroimaging features have demonstrated the capacity to distinguish cognitive decline changes stemming from typical brain aging and AD between different chronological points.

Objective

To disentangle the normal aging effect from the AD-related accelerated cognitive decline and unravel its genetic components using a neuroimaging-based deep learning approach.

Methods

We developed a deep-learning framework based on a dual-loss Siamese ResNet network to extract fine-grained information from the longitudinal structural magnetic resonance imaging (MRI) data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) study. We then conducted genome-wide association studies (GWAS) and post-GWAS analyses to reveal the genetic basis of AD-related accelerated cognitive decline.

Results

We used our model to process data from 1,313 individuals, training it on 414 cognitively normal people and predicting cognitive assessment for all participants. In our analysis of accelerated cognitive decline GWAS, we identified two genome-wide significant loci: APOE locus (chromosome 19 p13.32) and rs144614292 (chromosome 11 p15.1). Variant rs144614292 (G > T) has not been reported in previous AD GWA studies. It is within the intronic region of NELL1, which is expressed in neurons and plays a role in controlling cell growth and differentiation. The cell-type-specific enrichment analysis and functional enrichment of GWAS signals highlighted the microglia and immune-response pathways.

Conclusions

Our deep learning model effectively extracted relevant neuroimaging features and predicted individual cognitive decline. We reported a novel variant (rs144614292) within the NELL1 gene.

Decreased CNNM2 expression in prefrontal cortex affects sensorimotor gating function, cognition, dendritic spine morphogenesis and risk of schizophrenia

Genome-wide association studies (GWASs) have reported multiple single nucleotide polymorphisms (SNPs) associated with schizophrenia, yet the underlying molecular mechanisms are largely unknown. In this study, we aimed to identify schizophrenia relevant genes showing alterations in mRNA and protein expression associated with risk SNPs at the 10q24.32-33 GWAS locus. We carried out the quantitative trait loci (QTL) and summary data-based Mendelian randomization (SMR) analyses, using the PsychENCODE dorsolateral prefrontal cortex (DLPFC) expression QTL (eQTL) database, as well as the ROSMAP and Banner DLPFC protein QTL (pQTL) datasets. The gene CNNM2 (encoding a magnesium transporter) at 10q24.32-33 was identified to be a robust schizophrenia risk gene, and was highly expressed in human neurons according to single cell RNA-seq (scRNA-seq) data. We further revealed that reduced Cnnm2 in the mPFC of mice led to impaired cognition and compromised sensorimotor gating function, and decreased Cnnm2 in primary cortical neurons altered dendritic spine morphogenesis, confirming the link between CNNM2 and endophenotypes of schizophrenia. Proteomics analyses showed that reduced Cnnm2 level changed expression of proteins associated with neuronal structure and function. Together, these results identify a robust gene in the pathogenesis of schizophrenia.

Overview: Research on the Genetic Architecture of the Developing Cerebral Cortex in Norms and Diseases

The human brain is characterized by high cell numbers, diverse cell types with diverse functions, and intricate connectivity with an exceedingly broad surface of the cortex. Human-specific brain development was accomplished by a long timeline for maturation from the prenatal period to the third decade of life. The long timeline makes complicated architecture and circuits of human cerebral cortex possible, and it makes human brain vulnerable to intrinsic and extrinsic insults resulting in the development of variety of neuropsychiatric disorders. Unraveling the molecular and cellular processes underlying human brain development under the elaborate regulation of gene expression in a spatiotemporally specific manner, especially that of the cortex will provide a biological understanding of human cognition and behavior in health and diseases. Global research consortia and the advancing technologies in brain science including functional genomics equipped with emergent neuroinformatics such as single-cell multiomics, novel human models, and high-volume databases with high-throughput computation facilitate the biological understanding of the development of the human brain cortex. Knowing the process of interplay of the genome and the environment in cortex development will lead us to understand the human-specific cognitive function and its individual diversity. Thus, it is worthwhile to overview the recent progress in neurotechnology to foresee further understanding of the human brain and norms and diseases.

Common genetic variation impacts stress response in the brain

To explain why individuals exposed to identical stressors experience divergent clinical outcomes, we determine how molecular encoding of stress modifies genetic risk for brain disorders. Analysis of post-mortem brain (n=304) revealed 8557 stress-interactive expression quantitative trait loci (eQTLs) that dysregulate expression of 915 eGenes in response to stress, and lie in stress-related transcription factor binding sites. Response to stress is robust across experimental paradigms: up to 50% of stress-interactive eGenes validate in glucocorticoid treated hiPSC-derived neurons (n=39 donors). Stress-interactive eGenes show brain region- and cell type-specificity, and, in post-mortem brain, implicate glial and endothelial mechanisms. Stress dysregulates long-term expression of disorder risk genes in a genotype-dependent manner; stress-interactive transcriptomic imputation uncovered 139 novel genes conferring brain disorder risk only in the context of traumatic stress. Molecular stress-encoding explains individualized responses to traumatic stress; incorporating trauma into genomic studies of brain disorders is likely to improve diagnosis, prognosis, and drug discovery.

Mapping the landscape of lineage-specific dynamic regulation of gene expression using single-cell transcriptomics and application to genetics of complex disease

Single-cell transcriptome data can provide insights into how genetic variation influences biological processes involved in human biology and disease. However, the identification of gene-level associations in distinct cell types faces several challenges, including the limited reference resource from population scale studies, data sparsity in single-cell RNA sequencing, and the complex cell-state pattern of expression within individual cell types. Here we develop genetic models of cell type specific and cell state adjusted gene expression in mid-brain neurons in the process of specializing from induced pluripotent stem cells. The resulting framework quantifies the dynamics of the genetic regulation of gene expression and estimates its cell type specificity. As an application, we show that the approach detects known and new genes associated with schizophrenia and enables insights into context-dependent disease mechanisms. We provide a genomic resource from a phenome-wide application of our models to more than 1500 phenotypes from the UK Biobank. Using longitudinal genetically determined expression, we implement a predictive causality framework, evaluating the prediction of future values of a target gene expression using prior values of a putative regulatory gene. Collectively, this work demonstrates the insights that can be gained into the molecular underpinnings of diseases by quantifying the genetic control of gene expression at single-cell resolution.

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.

eQTL studies: from bulk tissues to single cells

An expression quantitative trait locus (eQTL) is a chromosomal region where genetic variants are associated with the expression levels of specific genes that can be both nearby or distant. The identifications of eQTLs for different tissues, cell types, and contexts have led to a better understanding of the dynamic regulations of gene expressions and implications of functional genes and variants for complex traits and diseases. Although most eQTL studies have been performed on data collected from bulk tissues, recent studies have demonstrated the importance of cell-type-specific and context-dependent gene regulations in biological processes and disease mechanisms. In this review, we discuss statistical methods that have been developed to enable the detection of cell-type-specific and context-dependent eQTLs from bulk tissues, purified cell types, and single cells. We also discuss the limitations of the current methods and future research opportunities.

Systems genetics approaches for understanding complex traits with relevance for human disease

Quantitative traits are often complex because of the contribution of many loci, with further complexity added by environmental factors. In medical research, systems genetics is a powerful approach for the study of complex traits, as it integrates intermediate phenotypes, such as RNA, protein, and metabolite levels, to understand molecular and physiological phenotypes linking discrete DNA sequence variation to complex clinical and physiological traits. The primary purpose of this review is to describe some of the resources and tools of systems genetics in humans and rodent models, so that researchers in many areas of biology and medicine can make use of the data.

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.

Genomic insights into the comorbidity between type 2 diabetes and schizophrenia

Multimorbidity represents an increasingly important public health challenge with far-reaching implications for health management and policy. Mental health and metabolic diseases have a well-established epidemiological association. In this study, we investigate the genetic intersection between type 2 diabetes and schizophrenia. We use Mendelian randomization to examine potential causal relationships between the two conditions and related endophenotypes. We report no compelling evidence that type 2 diabetes genetic liability potentially causally influences schizophrenia risk and vice versa. Our findings show that increased body mass index (BMI) has a protective effect against schizophrenia, in contrast to the well-known risk-increasing effect of BMI on type 2 diabetes risk. We identify evidence of colocalization of association signals for these two conditions at 11 genomic loci, six of which have opposing directions of effect for type 2 diabetes and schizophrenia. To elucidate these colocalizing signals, we integrate multi-omics data from bulk and single-cell gene expression studies, along with functional information. We identify high-confidence effector genes and find that they are enriched for homeostasis and lipid-related pathways. We also highlight drug repurposing opportunities including N-methyl-D-aspartate (NMDA) receptor antagonists. Our findings provide insights into shared biological mechanisms for type 2 diabetes and schizophrenia, highlighting common factors that influence the risk of the two conditions in opposite directions and shedding light on the complex nature of this comorbidity.

Interindividual variation in human cortical cell type abundance and expression

Single-cell transcriptomic studies have identified a conserved set of neocortical cell types from small postmortem cohorts. We extended these efforts by assessing cell type variation across 75 adult individuals undergoing epilepsy and tumor surgeries. Nearly all nuclei map to one of 125 robust cell types identified in the middle temporal gyrus. However, we found interindividual variance in abundances and gene expression signatures, particularly in deep-layer glutamatergic neurons and microglia. A minority of donor variance is explainable by age, sex, ancestry, disease state, and cell state. Genomic variation was associated with expression of 150 to 250 genes for most cell types. This characterization of cellular variation provides a baseline for cell typing in health and disease.

Causal inference in drug discovery and development

To discover new drugs is to seek and to prove causality. As an emerging approach leveraging human knowledge and creativity, data, and machine intelligence, causal inference holds the promise of reducing cognitive bias and improving decision-making in drug discovery. Although it has been applied across the value chain, the concepts and practice of causal inference remain obscure to many practitioners. This article offers a nontechnical introduction to causal inference, reviews its recent applications, and discusses opportunities and challenges of adopting the causal language in drug discovery and development.

Multitissue H3K27ac profiling of GTEx samples links epigenomic variation to disease

Genetic variants associated with complex traits are primarily noncoding, and their effects on gene-regulatory activity remain largely uncharacterized. To address this, we profile epigenomic variation of histone mark H3K27ac across 387 brain, heart, muscle and lung samples from Genotype-Tissue Expression (GTEx). We annotate 282 k active regulatory elements (AREs) with tissue-specific activity patterns. We identify 2,436 sex-biased AREs and 5,397 genetically influenced AREs associated with 130 k genetic variants (haQTLs) across tissues. We integrate genetic and epigenomic variation to provide mechanistic insights for disease-associated loci from 55 genome-wide association studies (GWAS), by revealing candidate tissues of action, driver SNPs and impacted AREs. Lastly, we build ARE-gene linking scores based on genetics (gLink scores) and demonstrate their unique ability to prioritize SNP-ARE-gene circuits. Overall, our epigenomic datasets, computational integration and mechanistic predictions provide valuable resources and important insights for understanding the molecular basis of human diseases/traits such as schizophrenia.

Single-Nucleus Transcriptome Profiling of Dorsolateral Prefrontal Cortex: Mechanistic Roles for Neuronal Gene Expression, Including the 17q21.31 Locus, in PTSD Stress Response

OBJECTIVE

Multidisciplinary studies of posttraumatic stress disorder (PTSD) and major depressive disorder (MDD) implicate the dorsolateral prefrontal cortex (DLPFC) in disease risk and pathophysiology. Postmortem brain studies have relied on bulk-tissue RNA sequencing (RNA-seq), but single-cell RNA-seq is needed to dissect cell-type-specific mechanisms. The authors conducted the first single-nucleus RNA-seq postmortem brain study in PTSD to elucidate disease transcriptomic pathology with cell-type-specific resolution.

METHOD

Profiling of 32 DLPFC samples from 11 individuals with PTSD, 10 with MDD, and 11 control subjects was conducted (∼415K nuclei; >13K cells per sample). A replication sample included 15 DLPFC samples (∼160K nuclei; >11K cells per sample).

RESULTS

Differential gene expression analyses identified significant single-nucleus RNA-seq differentially expressed genes (snDEGs) in excitatory (EX) and inhibitory (IN) neurons and astrocytes, but not in other cell types or bulk tissue. MDD samples had more false discovery rate-corrected significant snDEGs, and PTSD samples had a greater replication rate. In EX and IN neurons, biological pathways that were differentially enriched in PTSD compared with MDD included glucocorticoid signaling. Furthermore, glucocorticoid signaling in induced pluripotent stem cell (iPSC)-derived cortical neurons demonstrated greater relevance in PTSD and opposite direction of regulation compared with MDD, especially in EX neurons. Many snDEGs were from the 17q21.31 locus and are particularly interesting given causal roles in disease pathogenesis and DLPFC-based neuroimaging (PTSD: ARL17B, LINC02210-CRHR1, and LRRC37A2; MDD: LRRC37A and LRP4), while others were regulated by glucocorticoids in iPSC-derived neurons (PTSD: SLC16A6, TAF1C; MDD: CDH3).

CONCLUSIONS

The study findings point to cell-type-specific mechanisms of brain stress response in PTSD and MDD, highlighting the importance of examining cell-type-specific gene expression and indicating promising novel biomarkers and therapeutic targets.

Disentangling accelerated cognitive decline from the normal aging process and unraveling its genetic components: A neuroimaging-based deep learning approach

Background

The progressive cognitive decline that is an integral component of AD unfolds in tandem with the natural aging process. Neuroimaging features have demonstrated the capacity to distinguish cognitive decline changes stemming from typical brain aging and Alzheimer's disease between different chronological points.

Methods

We developed a deep-learning framework based on dual-loss Siamese ResNet network to extract fine-grained information from the longitudinal structural magnetic resonance imaging (MRI) data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) study. We then conducted genome-wide association studies (GWAS) and post-GWAS analyses to reveal the genetic basis of AD-related accelerated cognitive decline.

Results

We used our model to process data from 1,313 individuals, training it on 414 cognitively normal people and predicting cognitive assessment for all participants. In our analysis of accelerated cognitive decline GWAS, we identified two genome-wide significant loci: APOE locus (chromosome 19 p13.32) and rs144614292 (chromosome 11 p15.1). Variant rs144614292 (G>T) has not been reported in previous AD GWA studies. It is within the intronic region of NELL1, which is expressed in neuron and plays a role in controlling cell growth and differentiation. In addition, MUC7 and PROL1/OPRPNon chromosome 4 were significant at the gene level. The cell-type-specific enrichment analysis and functional enrichment of GWAS signals highlighted the microglia and immune-response pathways. Furthermore, we found that the cognitive decline slope GWAS was positively correlated with previous AD GWAS.

Conclusion

Our deep learning model was demonstrated effective on extracting relevant neuroimaging features and predicting individual cognitive decline. We reported a novel variant (rs144614292) within the NELL1 gene. Our approach has the potential to disentangle accelerated cognitive decline from the normal aging process and to determine its related genetic factors, leveraging opportunities for early intervention.

Discovery of 95 PTSD loci provides insight into genetic architecture and neurobiology of trauma and stress-related disorders

Posttraumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 novel). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (e.g., GRIA1, GRM8, CACNA1E ), developmental, axon guidance, and transcription factors (e.g., FOXP2, EFNA5, DCC ), synaptic structure and function genes (e.g., PCLO, NCAM1, PDE4B ), and endocrine or immune regulators (e.g., ESR1, TRAF3, TANK ). Additional top genes influence stress, immune, fear, and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk post-mortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA-editing and alternative polyadenylation (APA), within a cell-type-specific population of human neural progenitors and neurons. More RNA-editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting genetically mediated post-transcriptional regulation during brain development lead to differences in brain function.

Methods and Insights from Single-Cell Expression Quantitative Trait Loci

Recent advancements in single-cell technologies have enabled expression quantitative trait locus (eQTL) analysis across many individuals at single-cell resolution. Compared with bulk RNA sequencing, which averages gene expression across cell types and cell states, single-cell assays capture the transcriptional states of individual cells, including fine-grained, transient, and difficult-to-isolate populations at unprecedented scale and resolution. Single-cell eQTL (sc-eQTL) mapping can identify context-dependent eQTLs that vary with cell states, including some that colocalize with disease variants identified in genome-wide association studies. By uncovering the precise contexts in which these eQTLs act, single-cell approaches can unveil previously hidden regulatory effects and pinpoint important cell states underlying molecular mechanisms of disease. Here, we present an overview of recently deployed experimental designs in sc-eQTL studies. In the process, we consider the influence of study design choices such as cohort, cell states, and ex vivo perturbations. We then discuss current methodologies, modeling approaches, and technical challenges as well as future opportunities and applications.