The PsychENCODE project

Multi-omics in stress and health research: study designs that will drive the field forward

Despite decades of stress research, there still exist substantial gaps in our understanding of how social, environmental, and biological factors interact and combine with developmental stressor exposures, cognitive appraisals of stressors, and psychosocial coping processes to shape individuals' stress reactivity, health, and disease risk. Relatively new biological profiling approaches, called multi-omics, are helping address these issues by enabling researchers to quantify thousands of molecules from a single blood or tissue sample, thus providing a panoramic snapshot of the molecular processes occurring in an organism from a systems perspective. In this review, we summarize two types of research designs for which multi-omics approaches are best suited, and describe how these approaches can help advance our understanding of stress processes and the development, prevention, and treatment of stress-related pathologies. We first discuss incorporating multi-omics approaches into theory-rich, intensive longitudinal study designs to characterize, in high-resolution, the transition to stress-related multisystem dysfunction and disease throughout development. Next, we discuss how multi-omics approaches should be incorporated into intervention research to better understand the transition from stress-related dysfunction back to health, which can help inform novel precision medicine approaches to managing stress and fostering biopsychosocial resilience. Throughout, we provide concrete recommendations for types of studies that will help advance stress research, and translate multi-omics data into better health and health care.

ChIP-DIP: A multiplexed method for mapping hundreds of proteins to DNA uncovers diverse regulatory elements controlling gene expression

Gene expression is controlled by the dynamic localization of thousands of distinct regulatory proteins to precise regions of DNA. Understanding this cell-type specific process has been a goal of molecular biology for decades yet remains challenging because most current DNA-protein mapping methods study one protein at a time. To overcome this, we developed ChIP-DIP (ChIP Done In Parallel), a split-pool based method that enables simultaneous, genome-wide mapping of hundreds of diverse regulatory proteins in a single experiment. We demonstrate that ChIP-DIP generates highly accurate maps for all classes of DNA-associated proteins, including histone modifications, chromatin regulators, transcription factors, and RNA Polymerases. Using these data, we explore quantitative combinations of protein localization on genomic DNA to define distinct classes of regulatory elements and their functional activity. Our data demonstrate that ChIP-DIP enables the generation of 'consortium level', context-specific protein localization maps within any molecular biology lab.

Challenges and opportunities to computationally deconvolve heterogeneous tissue with varying cell sizes using single-cell RNA-sequencing datasets

Deconvolution of cell mixtures in "bulk" transcriptomic samples from homogenate human tissue is important for understanding disease pathologies. However, several experimental and computational challenges impede transcriptomics-based deconvolution approaches using single-cell/nucleus RNA-seq reference atlases. Cells from the brain and blood have substantially different sizes, total mRNA, and transcriptional activities, and existing approaches may quantify total mRNA instead of cell type proportions. Further, standards are lacking for the use of cell reference atlases and integrative analyses of single-cell and spatial transcriptomics data. We discuss how to approach these key challenges with orthogonal "gold standard" datasets for evaluating deconvolution methods.

Towards a biologically annotated brain connectome

The brain is a network of interleaved neural circuits. In modern connectomics, brain connectivity is typically encoded as a network of nodes and edges, abstracting away the rich biological detail of local neuronal populations. Yet biological annotations for network nodes - such as gene expression, cytoarchitecture, neurotransmitter receptors or intrinsic dynamics - can be readily measured and overlaid on network models. Here we review how connectomes can be represented and analysed as annotated networks. Annotated connectomes allow us to reconceptualize architectural features of networks and to relate the connection patterns of brain regions to their underlying biology. Emerging work demonstrates that annotated connectomes help to make more veridical models of brain network formation, neural dynamics and disease propagation. Finally, annotations can be used to infer entirely new inter-regional relationships and to construct new types of network that complement existing connectome representations. In summary, biologically annotated connectomes offer a compelling way to study neural wiring in concert with local biological features.

Artificial intelligence for dementia genetics and omics

Genetics and omics studies of Alzheimer's disease and other dementia subtypes enhance our understanding of underlying mechanisms and pathways that can be targeted. We identified key remaining challenges: First, can we enhance genetic studies to address missing heritability? Can we identify reproducible omics signatures that differentiate between dementia subtypes? Can high-dimensional omics data identify improved biomarkers? How can genetics inform our understanding of causal status of dementia risk factors? And which biological processes are altered by dementia-related genetic variation? Artificial intelligence (AI) and machine learning approaches give us powerful new tools in helping us to tackle these challenges, and we review possible solutions and examples of best practice. However, their limitations also need to be considered, as well as the need for coordinated multidisciplinary research and diverse deeply phenotyped cohorts. Ultimately AI approaches improve our ability to interrogate genetics and omics data for precision dementia medicine. HIGHLIGHTS: We have identified five key challenges in dementia genetics and omics studies. AI can enable detection of undiscovered patterns in dementia genetics and omics data. Enhanced and more diverse genetics and omics datasets are still needed. Multidisciplinary collaborative efforts using AI can boost dementia research.

A region-resolved proteomic map of the human brain enabled by high-throughput proteomics

Substantial efforts are underway to deepen our understanding of human brain morphology, structure, and function using high-resolution imaging as well as high-content molecular profiling technologies. The current work adds to these approaches by providing a comprehensive and quantitative protein expression map of 13 anatomically distinct brain regions covering more than 11,000 proteins. This was enabled by the optimization, characterization, and implementation of a high-sensitivity and high-throughput microflow liquid chromatography timsTOF tandem mass spectrometry system (LC-MS/MS) capable of analyzing more than 2,000 consecutive samples prepared from formalin-fixed paraffin embedded (FFPE) material. Analysis of this proteomic resource highlighted brain region-enriched protein expression patterns and functional protein classes, protein localization differences between brain regions and individual markers for specific areas. To facilitate access to and ease further mining of the data by the scientific community, all data can be explored online in a purpose-built R Shiny app (https://brain-region-atlas.proteomics.ls.tum.de).

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.

The molecular genetic landscape of human brain size variation

Human brain size changes dynamically through early development, peaks in adolescence, and varies up to 2-fold among adults. However, the molecular genetic underpinnings of interindividual variation in brain size remain unknown. Here, we leveraged postmortem brain RNA sequencing and measurements of brain weight (BW) in 2,531 individuals across three independent datasets to identify 928 genome-wide significant associations with BW. Genes associated with higher or lower BW showed distinct neurodevelopmental trajectories and spatial patterns that mapped onto functional and cellular axes of brain organization. Expression of BW genes was predictive of interspecies differences in brain size, and bioinformatic annotation revealed enrichment for neurogenesis and cell-cell communication. Genome-wide, transcriptome-wide, and phenome-wide association analyses linked BW gene sets to neuroimaging measurements of brain size and brain-related clinical traits. Cumulatively, these results represent a major step toward delineating the molecular pathways underlying human brain size variation in health and disease.

The Imageable Genome

Understanding human disease on a molecular level, and translating this understanding into targeted diagnostics and therapies are central tenets of molecular medicine1. Realizing this doctrine requires an efficient adaptation of molecular discoveries into the clinic. We present an approach to facilitate this process by describing the Imageable Genome, the part of the human genome whose expression can be assessed via molecular imaging. Using a deep learning-based hybrid human-AI pipeline, we bridge individual genes and their relevance in human diseases with specific molecular imaging methods. Cross-referencing the Imageable Genome with RNA-seq data from over 60,000 individuals reveals diagnostic, prognostic and predictive imageable genes for a wide variety of major human diseases. Having both the critical size and focus to be altered in its expression during the development and progression of any human disease, the Imageable Genome will generate new imaging tools that improve the understanding, diagnosis and management of human diseases.

The impact of common variants on gene expression in the human brain: from RNA to protein to schizophrenia risk

Background

The impact of genetic variants on gene expression has been intensely studied at the transcription level, yielding in valuable insights into the association between genes and the risk of complex disorders, such as schizophrenia (SCZ). However, the downstream impact of these variants and the molecular mechanisms connecting transcription variation to disease risk are not well understood.

Results

We quantitated ribosome occupancy in prefrontal cortex samples of the BrainGVEX cohort. Together with transcriptomics and proteomics data from the same cohort, we performed cis-Quantitative Trait Locus (QTL) mapping and identified 3,253 expression QTLs (eQTLs), 1,344 ribosome occupancy QTLs (rQTLs), and 657 protein QTLs (pQTLs) out of 7,458 genes quantitated in all three omics types from 185 samples. Of the eQTLs identified, only 34% have their effects propagated to the protein level. Further analysis on the effect size of prefrontal cortex eQTLs identified from an independent dataset showed clear post-transcriptional attenuation of eQTL effects. To investigate the biological relevance of the attenuated eQTLs, we identified 70 expression-specific QTLs (esQTLs), 51 ribosome-occupancy-specific QTLs (rsQTLs), and 107 protein-specific QTLs (psQTLs). Five of these omics-specific QTLs showed strong colocalization with SCZ GWAS signals, three of them are esQTLs. The limited number of GWAS colocalization discoveries from omics-specific QTLs and the apparent prevalence of eQTL attenuation prompted us to take a complementary approach to investigate the functional relevance of attenuated eQTLs. Using S-PrediXcan we identified 74 SCZ risk genes, 34% of which were novel, and 67% of these risk genes were replicated in a MR-Egger test. Notably, 52 out of 74 risk genes were identified using eQTL data and 70% of these SCZ-risk-gene-driving eQTLs show little to no evidence of driving corresponding variations at the protein level.

Conclusion

The effect of eQTLs on gene expression in the prefrontal cortex is commonly attenuated post-transcriptionally. Many of the attenuated eQTLs still correlate with SCZ GWAS signal. Further investigation is needed to elucidate a mechanistic link between attenuated eQTLs and SCZ disease risk.

Integrative analyses highlight functional regulatory variants associated with neuropsychiatric diseases

Noncoding variants of presumed regulatory function contribute to the heritability of neuropsychiatric disease. A total of 2,221 noncoding variants connected to risk for ten neuropsychiatric disorders, including autism spectrum disorder, attention deficit hyperactivity disorder, bipolar disorder, borderline personality disorder, major depression, generalized anxiety disorder, panic disorder, post-traumatic stress disorder, obsessive-compulsive disorder and schizophrenia, were studied in developing human neural cells. Integrating epigenomic and transcriptomic data with massively parallel reporter assays identified differentially-active single-nucleotide variants (daSNVs) in specific neural cell types. Expression-gene mapping, network analyses and chromatin looping nominated candidate disease-relevant target genes modulated by these daSNVs. Follow-up integration of daSNV gene editing with clinical cohort analyses suggested that magnesium transport dysfunction may increase neuropsychiatric disease risk and indicated that common genetic pathomechanisms may mediate specific symptoms that are shared across multiple neuropsychiatric diseases.

A single-cell multi-omic atlas spanning the adult rhesus macaque brain

Cataloging the diverse cellular architecture of the primate brain is crucial for understanding cognition, behavior, and disease in humans. Here, we generated a brain-wide single-cell multimodal molecular atlas of the rhesus macaque brain. Together, we profiled 2.58 M transcriptomes and 1.59 M epigenomes from single nuclei sampled from 30 regions across the adult brain. Cell composition differed extensively across the brain, revealing cellular signatures of region-specific functions. We also identified 1.19 M candidate regulatory elements, many previously unidentified, allowing us to explore the landscape of cis-regulatory grammar and neurological disease risk in a cell type-specific manner. Altogether, this multi-omic atlas provides an open resource for investigating the evolution of the human brain and identifying novel targets for disease interventions.

Genome-wide Analysis of Motor Progression in Parkinson Disease

Background and Objectives

The genetic basis of Parkinson disease (PD) motor progression is largely unknown. Previous studies of the genetics of PD progression have included small cohorts and shown a limited overlap with genetic PD risk factors from case-control studies. Here, we have studied genomic variation associated with PD motor severity and early-stage progression in large longitudinal cohorts to help to define the biology of PD progression and potential new drug targets.

Methods

We performed a GWAS meta-analysis of early PD motor severity and progression up to 3 years from study entry. We used linear mixed-effect models with additive effects, corrected for age at diagnosis, sex, and the first 5 genetic principal components to assess variability in axial, limb, and total Movement Disorder Society-Unified Parkinson's Disease Rating Scale (MDS-UPDRS) III scores.

Results

We included 3,572 unrelated European ancestry patients with PD from 5 observational cohorts and 1 drug trial. The average AAO was 62.6 years (SD = 9.83), and 63% of participants were male. We found an average increase in the total MDS-UPDRS III score of 2.3 points/year. We identified an association between PD axial motor progression and variation at the GJA5 locus at 1q12 (β = -0.25, SE = 0.04, p = 3.4e-10). Exploration of the regulation of gene expression in the region (cis-expression quantitative trait loci [eQTL] analysis) showed that the lead variant was associated with expression of ACP6, a lysophosphatidic acid phosphatase that regulates mitochondrial lipid biosynthesis (cis-eQTL p-values in blood and brain RNA expression data sets: <10-14 in eQTLGen and 10-7 in PsychEncode).

Discussion

Our study highlights the potential role of mitochondrial lipid homeostasis in the progression of PD, which may be important in establishing new drug targets that might modify disease progression.

Epigenomic dissection of Alzheimer's disease pinpoints causal variants and reveals epigenome erosion

Recent work has identified dozens of non-coding loci for Alzheimer's disease (AD) risk, but their mechanisms and AD transcriptional regulatory circuitry are poorly understood. Here, we profile epigenomic and transcriptomic landscapes of 850,000 nuclei from prefrontal cortexes of 92 individuals with and without AD to build a map of the brain regulome, including epigenomic profiles, transcriptional regulators, co-accessibility modules, and peak-to-gene links in a cell-type-specific manner. We develop methods for multimodal integration and detecting regulatory modules using peak-to-gene linking. We show AD risk loci are enriched in microglial enhancers and for specific TFs including SPI1, ELF2, and RUNX1. We detect 9,628 cell-type-specific ATAC-QTL loci, which we integrate alongside peak-to-gene links to prioritize AD variant regulatory circuits. We report differential accessibility of regulatory modules in late AD in glia and in early AD in neurons. Strikingly, late-stage AD brains show global epigenome dysregulation indicative of epigenome erosion and cell identity loss.

Genetic insights into human cortical organization and development through genome-wide analyses of 2,347 neuroimaging phenotypes

Our understanding of the genetics of the human cerebral cortex is limited both in terms of the diversity and the anatomical granularity of brain structural phenotypes. Here we conducted a genome-wide association meta-analysis of 13 structural and diffusion magnetic resonance imaging-derived cortical phenotypes, measured globally and at 180 bilaterally averaged regions in 36,663 individuals and identified 4,349 experiment-wide significant loci. These phenotypes include cortical thickness, surface area, gray matter volume, measures of folding, neurite density and water diffusion. We identified four genetic latent structures and causal relationships between surface area and some measures of cortical folding. These latent structures partly relate to different underlying gene expression trajectories during development and are enriched for different cell types. We also identified differential enrichment for neurodevelopmental and constrained genes and demonstrate that common genetic variants associated with cortical expansion are associated with cephalic disorders. Finally, we identified complex interphenotype and inter-regional genetic relationships among the 13 phenotypes, reflecting the developmental differences among them. Together, these analyses identify distinct genetic organizational principles of the cortex and their correlates with neurodevelopment.

Whole-genome sequencing analysis of suicide deaths integrating brain-regulatory eQTLs data to identify risk loci and genes

Recent large-scale genome-wide association studies (GWAS) have started to identify potential genetic risk loci associated with risk of suicide; however, a large portion of suicide-associated genetic factors affecting gene expression remain elusive. Dysregulated gene expression, not assessed by GWAS, may play a significant role in increasing the risk of suicide death. We performed the first comprehensive genomic association analysis prioritizing brain expression quantitative trait loci (eQTLs) within regulatory regions in suicide deaths from the Utah Suicide Genetic Risk Study (USGRS). 440,324 brain-regulatory eQTLs were obtained by integrating brain eQTLs, histone modification ChIP-seq, ATAC-seq, DNase-seq, and Hi-C results from publicly available data. Subsequent genomic analyses were conducted in whole-genome sequencing (WGS) data from 986 suicide deaths of non-Finnish European (NFE) ancestry and 415 ancestrally matched controls. Additional independent USGRS suicide deaths with genotyping array data (n = 4657) and controls from the Genome Aggregation Database were explored for WGS result replication. One significant eQTL locus, rs926308 (p = 3.24e-06), was identified. The rs926308-T is associated with lower expression of RFPL3S, a gene important for neocortex development and implicated in arousal. Gene-based analyses performed using Sherlock Bayesian statistical integrative analysis also detected 20 genes with expression changes that may contribute to suicide risk. From analyzing publicly available transcriptomic data, ten of these genes have previous evidence of differential expression in suicide death or in psychiatric disorders that may be associated with suicide, including schizophrenia and autism (ZNF501, ZNF502, CNN3, IGF1R, KLHL36, NBL1, PDCD6IP, SNX19, BCAP29, and ARSA). Electronic health records (EHR) data was further merged to evaluate if there were clinically relevant subsets of suicide deaths associated with genetic variants. In summary, our study identified one risk locus and ten genes associated with suicide risk via gene expression, providing new insight into possible genetic and molecular mechanisms leading to suicide.

Human forebrain organoid-based multi-omics analyses of PCCB as a schizophrenia associated gene linked to GABAergic pathways

Identifying genes whose expression is associated with schizophrenia (SCZ) risk by transcriptome-wide association studies (TWAS) facilitates downstream experimental studies. Here, we integrated multiple published datasets of TWAS, gene coexpression, and differential gene expression analysis to prioritize SCZ candidate genes for functional study. Convergent evidence prioritized Propionyl-CoA Carboxylase Subunit Beta (PCCB), a nuclear-encoded mitochondrial gene, as an SCZ risk gene. However, the PCCB's contribution to SCZ risk has not been investigated before. Using dual luciferase reporter assay, we identified that SCZ-associated SNPs rs6791142 and rs35874192, two eQTL SNPs for PCCB, showed differential allelic effects on transcriptional activities. PCCB knockdown in human forebrain organoids (hFOs) followed by RNA sequencing analysis revealed dysregulation of genes enriched with multiple neuronal functions including gamma-aminobutyric acid (GABA)-ergic synapse. The metabolomic and mitochondrial function analyses confirmed the decreased GABA levels resulted from inhibited tricarboxylic acid cycle in PCCB knockdown hFOs. Multielectrode array recording analysis showed that PCCB knockdown in hFOs resulted into SCZ-related phenotypes including hyper-neuroactivities and decreased synchronization of neural network. In summary, this study utilized hFOs-based multi-omics analyses and revealed that PCCB downregulation may contribute to SCZ risk through regulating GABAergic pathways, highlighting the mitochondrial function in SCZ.

Pathogenic Mis-splicing of CPEB4 in Schizophrenia

BACKGROUND

Schizophrenia (SCZ) is caused by an interplay of polygenic risk and environmental factors, which may alter regulators of gene expression leading to pathogenic misexpression of SCZ risk genes. The CPEB family of RNA-binding proteins (CPEB1-4) regulates translation of target RNAs (approximately 40% of overall genes). We previously identified CPEB4 as a key dysregulated translational regulator in autism spectrum disorder (ASD) because its neuronal-specific microexon (exon 4) is mis-spliced in ASD brains, causing underexpression of numerous ASD risk genes. The genetic factors and pathogenic mechanisms shared between SCZ and ASD led us to hypothesize CPEB4 mis-splicing in SCZ leading to underexpression of multiple SCZ-related genes.

METHODS

We performed MAGMA-enrichment analysis on Psychiatric Genomics Consortium genome-wide association study data and analyzed RNA sequencing data from the PsychENCODE Consortium. Reverse transcriptase polymerase chain reaction and Western blot were performed on postmortem brain tissue, and the presence/absence of antipsychotics was assessed through toxicological analysis. Finally, mice with mild overexpression of exon 4-lacking CPEB4 (CPEB4Δ4) were generated and analyzed biochemically and behaviorally.

RESULTS

First, we found enrichment of SCZ-associated genes for CPEB4-binder transcripts. We also found decreased usage of CPEB4 microexon in SCZ probands, which was correlated with decreased protein levels of CPEB4-target SCZ-associated genes only in antipsychotic-free individuals. Interestingly, differentially expressed genes fit those reported for SCZ, specifically in the SCZ probands with decreased CPEB4-microexon inclusion. Finally, we demonstrated that mice with mild overexpression of CPEB4Δ4 showed decreased protein levels of CPEB4-target SCZ genes and SCZ-linked behaviors.

CONCLUSIONS

We identified aberrant CPEB4 splicing and downstream misexpression of SCZ risk genes as a novel etiological mechanism in SCZ.

Autistic traits in myotonic dystrophy type 1 due to MBNL inhibition and RNA mis-splicing

Tandem repeat expansions are enriched in autism spectrum disorder, including CTG expansion in the DMPK gene that underlines myotonic muscular dystrophy type 1. Although the clinical connection of autism to myotonic dystrophy is corroborated, the molecular links remained unknown. Here, we show a mechanistic path of autism via repeat expansion in myotonic dystrophy. We found that inhibition of muscleblind-like (MBNL) splicing factors by expanded CUG RNAs alerts the splicing of autism-risk genes during brain development especially a class of autism-relevant microexons. To provide in vivo evidence that the CTG expansion and MBNL inhibition axis leads to the presentation of autistic traits, we demonstrate that CTG expansion and MBNL-null mouse models recapitulate autism-relevant mis-splicing profiles and demonstrate social deficits. Our findings indicate that DMPK CTG expansion-associated autism arises from developmental mis-splicing. Understanding this pathomechanistic connection provides an opportunity for greater in-depth investigations of mechanistic threads in autism.

Identification of drug repurposing candidates for the treatment of anxiety: A genetic approach

Anxiety disorders are a group of prevalent and heritable neuropsychiatric diseases. We previously conducted a genome-wide association study (GWAS) which identified genomic loci associated with anxiety; however, the biological consequences underlying the genetic associations are largely unknown. Integrating GWAS and functional genomic data may improve our understanding of the genetic effects on intermediate molecular phenotypes such as gene expression. This can provide an opportunity for the discovery of drug targets for anxiety via drug repurposing. We used the GWAS summary statistics to determine putative causal genes for anxiety using MAGMA and colocalization analyses. A transcriptome-wide association study was conducted to identify genes with differential genetically regulated levels of gene expression in human brain tissue. The genes were integrated with a large drug-gene expression database (Connectivity Map), discovering compounds that are predicted to "normalise" anxiety-associated expression changes. The study identified 64 putative causal genes associated with anxiety (35 genes upregulated; 29 genes downregulated). Drug mechanisms adrenergic receptor agonists, sigma receptor agonists, and glutamate receptor agonists gene targets were enriched in anxiety-associated genetic signal and exhibited an opposing effect on the anxiety-associated gene expression signature. The significance of the project demonstrated genetic links for novel drug candidates to potentially advance anxiety therapeutics.

New "drugs and targets" in the GWAS era of bipolar disorder

OBJECTIVE

Due to the phenotypic heterogeneity and etiological complexity of bipolar disorder (BD), many patients do not respond well to the current medications, and developing novel effective treatment is necessary. Whether any BD genome-wide association study (GWAS) risk genes were targets of existing drugs or novel drugs that can be repurposed in the clinical treatment of BD is a hot topic in the GWAS era of BD.

METHODS

A list of 425 protein-coding BD risk genes was distilled through the BD GWAS, and 4479 protein-coding druggable targets were retrieved from the druggable genome. The overlapped genes/targets were subjected to further analyses in DrugBank, Pharos, and DGIdb datasets in terms of their FDA status, mechanism of action and primary indication, to identify their potential for repurposing.

RESULTS

We identified 58 BD GWAS risk genes grouped as the druggable targets, and several genes were given higher priority. These BD risk genes were targets of antipsychotics, antidepressants, antiepileptics, calcium channel antagonists, as well as anxiolytics and analgesics, either existing clinically-approved drugs for BD or the drugs than can be repurposed for treatment of BD in the future. Those genes were also likely relevant to BD pathophysiology, as many of them encode ion channel, ion transporter or neurotransmitter receptor, or the mice manipulating those genes are likely to mimic the phenotypes manifest in BD patients.

CONCLUSIONS

This study identifies several targets that may facilitate the discovery of novel treatments in BD, and implies the value of conducting GWAS into clinical translation.

Single-Cell RNA Sequencing and Its Applications in the Study of Psychiatric Disorders

Neuroscience is currently one of the most challenging research fields owing to the enormous complexity of the mammalian nervous system. We are yet to understand precise transcriptional programs that govern cell fate during neurodevelopment, resolve the connectome of the mammalian brain, and determine the etiology of various neurodegenerative and psychiatric disorders. Technological advances in the past decade, notably single-cell RNA sequencing, have enabled huge progress in our understanding of such features. Our current knowledge of the transcriptome is largely derived from bulk RNA sequencing, which reveals only the average gene expression of millions of cells, potentially missing out on minor transcriptome differences between cells detectable only at single-cell resolution. Since 2009, several single-cell RNA sequencing techniques have emerged that enable the accurate classification of neuronal and glial cell subtypes beyond classical molecular markers and electrophysiological features and allow the identification of previously unknown cell types. Furthermore, it enables the interrogation of molecular and disease-relevant mechanisms and offers further possibilities for the discovery of new drug targets and disease biomarkers. This review intends to familiarize the reader with the main single-cell RNA sequencing techniques developed throughout the past decade and discusses their application in the fields of brain cell taxonomy, neurodevelopment, and psychiatric disorders.

Integrating GWAS and proteome data to identify novel drug targets for MU

Mouth ulcers have been associated with numerous loci in genome wide association studies (GWAS). Nonetheless, it remains unclear what mechanisms are involved in the pathogenesis of mouth ulcers at these loci, as well as what the most effective ulcer drugs are. Thus, we aimed to screen hub genes responsible for mouth ulcer pathogenesis. We conducted an imputed/in-silico proteome-wide association study to discover candidate genes that impact the development of mouth ulcers and affect the expression and concentration of associated proteins in the bloodstream. The integrative analysis revealed that 35 genes play a significant role in the development of mouth ulcers, both in terms of their protein and transcriptional levels. Following this analysis, the researchers identified 6 key genes, namely BTN3A3, IL12B, BPI, FAM213A, PLXNB2, and IL22RA2, which were related to the onset of mouth ulcers. By combining with multidimensional data, six genes were found to correlate with mouth ulcer pathogenesis, which can be useful for further biological and therapeutic research.

Cell type deconvolution of bulk blood RNA-Seq to reveal biological insights of neuropsychiatric disorders

Genome-wide association studies (GWAS) have uncovered susceptibility loci associated with psychiatric disorders like bipolar disorder (BP) and schizophrenia (SCZ). However, most of these loci are in non-coding regions of the genome with unknown causal mechanisms of the link between genetic variation and disease risk. Expression quantitative trait loci (eQTL) analysis of bulk tissue is a common approach to decipher underlying mechanisms, though this can obscure cell-type specific signals thus masking trait-relevant mechanisms. While 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 for individuals with BP and SCZ this study estimated cell type proportions and their relation with disease status and medication. We found between 2,875 and 4,629 eGenes for each cell type, including 1,211 eGenes that are not found using bulk expression alone. We performed a colocalization test between cell type eQTLs and various traits and identified hundreds of associations between cell type eQTLs and GWAS loci that are not detected in bulk eQTLs. Finally, we investigated the effects of lithium use on cell type expression regulation and found examples of genes that are differentially regulated dependent on lithium use. Our study suggests that computational methods can be applied to large bulk RNA-Seq datasets of non-brain tissue to identify disease-relevant, cell type specific biology of psychiatric disorders and psychiatric medication.

Challenges and opportunities to computationally deconvolve heterogeneous tissue with varying cell sizes using single cell RNA-sequencing datasets

Deconvolution of cell mixtures in "bulk" transcriptomic samples from homogenate human tissue is important for understanding the pathologies of diseases. However, several experimental and computational challenges remain in developing and implementing transcriptomics-based deconvolution approaches, especially those using a single cell/nuclei RNA-seq reference atlas, which are becoming rapidly available across many tissues. Notably, deconvolution algorithms are frequently developed using samples from tissues with similar cell sizes. However, brain tissue or immune cell populations have cell types with substantially different cell sizes, total mRNA expression, and transcriptional activity. When existing deconvolution approaches are applied to these tissues, these systematic differences in cell sizes and transcriptomic activity confound accurate cell proportion estimates and instead may quantify total mRNA content. Furthermore, there is a lack of standard reference atlases and computational approaches to facilitate integrative analyses, including not only bulk and single cell/nuclei RNA-seq data, but also new data modalities from spatial -omic or imaging approaches. New multi-assay datasets need to be collected with orthogonal data types generated from the same tissue block and the same individual, to serve as a "gold standard" for evaluating new and existing deconvolution methods. Below, we discuss these key challenges and how they can be addressed with the acquisition of new datasets and approaches to analysis.

Characterization of Gene Regulatory Elements in Human Fetal Cortical Development: Enhancing Our Understanding of Neurodevelopmental Disorders and Evolution

The neocortex is the region that most distinguishes human brain from other mammals and primates [Annu Rev Genet. 2021 Nov;55(1):555-81]. Studying the development of human cortex is important in understanding the evolutionary changes occurring in humans relative to other primates, as well as in elucidating mechanisms underlying neurodevelopmental disorders. Cortical development is a highly regulated process, spatially and temporally coordinated by expression of essential transcriptional factors in response to signaling pathways [Neuron. 2019 Sep;103(6):980-1004]. Enhancers are the most well-understood cis-acting, non-protein-coding regulatory elements that regulate gene expression [Nat Rev Genet. 2014 Apr;15(4):272-86]. Importantly, given the conservation of both DNA sequence and molecular function of the majority of proteins across mammals [Genome Res. 2003 Dec;13(12):2507-18], enhancers [Science. 2015 Mar;347(6226):1155-9], which are far more divergent at the sequence level, likely account for the phenotypes that distinguish the human brain by changing the regulation of gene expression. In this review, we will revisit the conceptual framework of gene regulation during human brain development, as well as the evolution of technologies to study transcriptional regulation, with recent advances in genome biology that open a window allowing us to systematically characterize cis-regulatory elements in developing human brain [Hum Mol Genet. 2022 Oct;31(R1):R84-96]. We provide an update on work to characterize the suite of all enhancers in the developing human brain and the implications for understanding neuropsychiatric disorders. Finally, we discuss emerging therapeutic ideas that utilize our emerging knowledge of enhancer function.

Widespread RNA hypoediting in schizophrenia and its relevance to mitochondrial function

RNA editing, the endogenous modification of nucleic acids, is known to be altered in genes with important neurological function in schizophrenia (SCZ). However, the global profile and molecular functions of disease-associated RNA editing remain unclear. Here, we analyzed RNA editing in postmortem brains of four SCZ cohorts and uncovered a significant and reproducible trend of hypoediting in patients of European descent. We report a set of SCZ-associated editing sites via WGCNA analysis, shared across cohorts. Using massively parallel reporter assays and bioinformatic analyses, we observed that differential 3' untranslated region (3'UTR) editing sites affecting host gene expression were enriched for mitochondrial processes. Furthermore, we characterized the impact of two recoding sites in the mitofusin 1 (MFN1) gene and showed their functional relevance to mitochondrial fusion and cellular apoptosis. Our study reveals a global reduction of editing in SCZ and a compelling link between editing and mitochondrial function in the disease.

Proprotein convertase subtilisin/kexin type 9 (PCSK9) in the central nervous system

The gene encoding proprotein convertase subtilisin/kexin type 9 (PCSK9) and its protein product have been widely studied for their role in cholesterol and lipid metabolism. PCSK9 increases the rate of metabolic degradation of low-density lipoprotein receptors, preventing the diffusion of low-density lipoprotein (LDL) from plasma into cells and contributes to high lipoprotein-bound cholesterol levels in the plasma. While most research has focused on the regulation and disease relevance of PCSK9 to the cardiovascular system and lipid metabolism, there is a growing body of evidence that PCSK9 plays a crucial role in pathogenic processes in other organ systems, including the central nervous system. PCSK9's impact on the brain is not yet fully understood, though several recent studies have sought to illuminate its impact on various neurodegenerative and psychiatric disorders, as well as its connection with ischemic stroke. Cerebral PCSK9 expression is low but is highly upregulated during disease states. Among others, PCSK9 is known to play a role in neurogenesis, neural cell differentiation, central LDL receptor metabolism, neural cell apoptosis, neuroinflammation, Alzheimer's Disease, Alcohol Use Disorder, and stroke. The PCSK9 gene contains several polymorphisms, including both gain-of-function and loss-of-function mutations which profoundly impact normal PCSK9 signaling and cholesterol metabolism. Gain-of-function mutations lead to persistent hypercholesterolemia and poor health outcomes, while loss-of-function mutations generally lead to hypocholesterolemia and may serve as a protective factor against diseases of the liver, cardiovascular system, and central nervous system. Recent genomic studies have sought to identify the end-organ effects of such mutations and continue to identify evidence of a much broader role for PCSK9 in extrahepatic organ systems. Despite this, there remain large gaps in our understanding of PCSK9, its regulation, and its effects on disease risk outside the liver. This review, which incorporates data from a wide range of scientific disciplines and experimental paradigms, is intended to describe PCSK9's role in the central nervous system as it relates to cerebral disease and neuropsychiatric disorders, and to examine the clinical potential of PCSK9 inhibitors and genetic variation in the PCSK9 gene on disease outcomes, including neurological and neuropsychiatric disease.

Human forebrain organoids-based multi-omics analyses reveal PCCB's regulation on GABAergic system contributing to schizophrenia

Identifying genes whose expression is associated with schizophrenia (SCZ) risk by transcriptome-wide association studies (TWAS) facilitates downstream experimental studies. Here, we integrated multiple published datasets of TWAS (including FUSION, PrediXcan, summary-data-based Mendelian randomization (SMR), joint-tissue imputation approach with Mendelian randomization (MR-JTI)), gene coexpression, and differential gene expression analysis to prioritize SCZ candidate genes for functional study. Convergent evidence prioritized Propionyl-CoA Carboxylase Subunit Beta ( PCCB ), a nuclear-encoded mitochondrial gene, as an SCZ risk gene. However, the PCCB ’s contribution to SCZ risk has not been investigated before. Using dual luciferase reporter assay, we identified that SCZ-associated SNP rs35874192, an eQTL SNP for PCCB , showed differential allelic effects on transcriptional activities. PCCB knockdown in human forebrain organoids (hFOs) followed by RNA-seq revealed dysregulation of genes enriched with multiple neuronal functions including gamma-aminobutyric acid (GABA)-ergic synapse, as well as genes dysregulated in postmortem brains of SCZ patients or in cerebral organoids derived from SCZ patients. The metabolomic and mitochondrial function analyses confirmed the deceased GABA levels resulted from reduced tricarboxylic acid cycle in PCCB knockdown hFOs. Multielectrode array recording analysis showed that PCCB knockdown in hFOs resulted into SCZ-related phenotypes including hyper-neuroactivities and decreased synchronization of neural network. In summary, this study utilized hFOs-based multi-omics data and revealed that PCCB downregulation may contribute to SCZ risk through regulating GABAergic system, highlighting the mitochondrial function in SCZ.

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.

The Foundational Data Initiative for Parkinson Disease: Enabling efficient translation from genetic maps to mechanism

The Foundational Data Initiative for Parkinson Disease (FOUNDIN-PD) is an international collaboration producing fundamental resources for Parkinson disease (PD). FOUNDIN-PD generated a multi-layered molecular dataset in a cohort of induced pluripotent stem cell (iPSC) lines differentiated to dopaminergic (DA) neurons, a major affected cell type in PD. The lines were derived from the Parkinson's Progression Markers Initiative study, which included participants with PD carrying monogenic PD variants, variants with intermediate effects, and variants identified by genome-wide association studies and unaffected individuals. We generated genetic, epigenetic, regulatory, transcriptomic, and longitudinal cellular imaging data from iPSC-derived DA neurons to understand molecular relationships between disease-associated genetic variation and proximate molecular events. These data reveal that iPSC-derived DA neurons provide a valuable cellular context and foundational atlas for modeling PD genetic risk. We have integrated these data into a FOUNDIN-PD data browser as a resource for understanding the molecular pathogenesis of PD.

Multivariate genome-wide association meta-analysis of over 1 million subjects identifies loci underlying multiple substance use disorders

Genetic liability to substance use disorders can be parsed into loci that confer general or substance-specific addiction risk. We report a multivariate genome-wide association meta-analysis that disaggregates general and substance-specific loci for published summary statistics of problematic alcohol use, problematic tobacco use, cannabis use disorder, and opioid use disorder in a sample of 1,025,550 individuals of European descent and 92,630 individuals of African descent. Nineteen independent SNPs were genome-wide significant (P < 5e-8) for the general addiction risk factor (addiction-rf), which showed high polygenicity. Across ancestries, PDE4B was significant (among other genes), suggesting dopamine regulation as a cross-substance vulnerability. An addiction-rf polygenic risk score was associated with substance use disorders, psychopathologies, somatic conditions, and environments associated with the onset of addictions. Substance-specific loci (9 for alcohol, 32 for tobacco, 5 for cannabis, 1 for opioids) included metabolic and receptor genes. These findings provide insight into genetic risk loci for substance use disorders that could be leveraged as treatment targets.

Investigation of convergent and divergent genetic influences underlying schizophrenia and alcohol use disorder

BACKGROUND

Alcohol use disorder (AUD) and schizophrenia (SCZ) frequently co-occur, and large-scale genome-wide association studies (GWAS) have identified significant genetic correlations between these disorders.

METHODS

We used the largest published GWAS for AUD (total cases = 77 822) and SCZ (total cases = 46 827) to identify genetic variants that influence both disorders (with either the same or opposite direction of effect) and those that are disorder specific.

RESULTS

We identified 55 independent genome-wide significant single nucleotide polymorphisms with the same direction of effect on AUD and SCZ, 8 with robust effects in opposite directions, and 98 with disorder-specific effects. We also found evidence for 12 genes whose pleiotropic associations with AUD and SCZ are consistent with mediation via gene expression in the prefrontal cortex. The genetic covariance between AUD and SCZ was concentrated in genomic regions functional in brain tissues (p = 0.001).

CONCLUSIONS

Our findings provide further evidence that SCZ shares meaningful genetic overlap with AUD.

Multi-Omics Studies in Historically Excluded Populations: The Road to Equity

Over the past few decades, genomewide association studies (GWASs) have identified the specific genetics variants contributing to many complex diseases by testing millions of genetic variations across the human genome against a variety of phenotypes. However, GWASs are limited in their ability to uncover mechanistic insight given that most significant associations are found in non-coding region of the genome. Furthermore, the lack of diversity in studies has stymied the advance of precision medicine for many historically excluded populations. In this review, we summarize most popular multi-omics approaches (genomics, transcriptomics, proteomics, and metabolomics) related to precision medicine and highlight if diverse populations have been included and how their findings have advance biological understanding of disease and drug response. New methods that incorporate local ancestry have been to improve the power of GWASs for admixed populations (such as African Americans and Latinx). Because most signals from GWAS are in the non-coding region, other machine learning and omics approaches have been developed to identify the potential causative single-nucleotide polymorphisms and genes that explain these phenotypes. These include polygenic risk scores, expression quantitative trait locus mapping, and transcriptome-wide association studies. Analogous protein methods, such as proteins quantitative trait locus mapping, proteome-wide association studies, and metabolomic approaches provide insight into the consequences of genetic variation on protein abundance. Whereas, integrated multi-omics studies have improved our understanding of the mechanisms for genetic association, we still lack the datasets and cohorts for historically excluded populations to provide equity in precision medicine and pharmacogenomics.

Toward Best Practices for Imaging Transcriptomics of the Human Brain

Modern brainwide transcriptional atlases provide unprecedented opportunities for investigating the molecular correlates of brain organization, as quantified using noninvasive neuroimaging. However, integrating neuroimaging data with transcriptomic measures is not straightforward, and careful consideration is required to make valid inferences. In this article, we review recent work exploring how various methodological choices affect 3 main phases of imaging transcriptomic analyses, including 1) processing of transcriptional atlas data; 2) relating transcriptional measures to independently derived neuroimaging phenotypes; and 3) evaluating the functional implications of identified associations through gene enrichment analyses. Our aim is to facilitate the development of standardized and reproducible approaches for this rapidly growing field. We identify sources of methodological variability, key choices that can affect findings, and considerations for mitigating false positive and/or spurious results. Finally, we provide an overview of freely available open-source toolboxes implementing current best-practice procedures across all 3 analysis phases.

BOMA, a machine-learning framework for comparative gene expression analysis across brains and organoids

Our machine-learning framework, brain and organoid manifold alignment (BOMA), first performs a global alignment of developmental gene expression data between brains and organoids. It then applies manifold learning to locally refine the alignment, revealing conserved and specific developmental trajectories across brains and organoids. Using BOMA, we found that human cortical organoids better align with certain brain cortical regions than with other non-cortical regions, implying organoid-preserved developmental gene expression programs specific to brain regions. Additionally, our alignment of non-human primate and human brains reveals highly conserved gene expression around birth. Also, we integrated and analyzed developmental single-cell RNA sequencing (scRNA-seq) data of human brains and organoids, showing conserved and specific cell trajectories and clusters. Further identification of expressed genes of such clusters and enrichment analyses reveal brain- or organoid-specific developmental functions and pathways. Finally, we experimentally validated important specific expressed genes through the use of immunofluorescence. BOMA is open-source available as a web tool for community use.

Ultra-low-coverage genome-wide association study-insights into gestational age using 17,844 embryo samples with preimplantation genetic testing

BACKGROUND

Very low-coverage (0.1 to 1×) whole genome sequencing (WGS) has become a promising and affordable approach to discover genomic variants of human populations for genome-wide association study (GWAS). To support genetic screening using preimplantation genetic testing (PGT) in a large population, the sequencing coverage goes below 0.1× to an ultra-low level. However, the feasibility and effectiveness of ultra-low-coverage WGS (ulcWGS) for GWAS remains undetermined.

METHODS

We built a pipeline to carry out analysis of ulcWGS data for GWAS. To examine its effectiveness, we benchmarked the accuracy of genotype imputation at the combination of different coverages below 0.1× and sample sizes from 2000 to 16,000, using 17,844 embryo PGT samples with approximately 0.04× average coverage and the standard Chinese sample HG005 with known genotypes. We then applied the imputed genotypes of 1744 transferred embryos who have gestational ages and complete follow-up records to GWAS.

RESULTS

The accuracy of genotype imputation under ultra-low coverage can be improved by increasing the sample size and applying a set of filters. From 1744 born embryos, we identified 11 genomic risk loci associated with gestational ages and 166 genes mapped to these loci according to positional, expression quantitative trait locus, and chromatin interaction strategies. Among these mapped genes, CRHBP, ICAM1, and OXTR were more frequently reported as preterm birth related. By joint analysis of gene expression data from previous studies, we constructed interrelationships of mainly CRHBP, ICAM1, PLAGL1, DNMT1, CNTLN, DKK1, and EGR2 with preterm birth, infant disease, and breast cancer.

CONCLUSIONS

This study not only demonstrates that ulcWGS could achieve relatively high accuracy of adequate genotype imputation and is capable of GWAS, but also provides insights into the associations between gestational age and genetic variations of the fetal embryos from Chinese population.

Neuroimmune transcriptome changes in patient brains of psychiatric and neurological disorders

Neuroinflammation has been implicated in multiple brain disorders but the extent and the magnitude of change in immune-related genes (IRGs) across distinct brain disorders has not been directly compared. In this study, 1275 IRGs were curated and their expression changes investigated in 2467 postmortem brains of controls and patients with six major brain disorders, including schizophrenia (SCZ), bipolar disorder (BD), autism spectrum disorder (ASD), major depressive disorder (MDD), Alzheimer's disease (AD), and Parkinson's disease (PD). There were 865 IRGs present across all microarray and RNA-seq datasets. More than 60% of the IRGs had significantly altered expression in at least one of the six disorders. The differentially expressed immune-related genes (dIRGs) shared across disorders were mainly related to innate immunity. Moreover, sex, tissue, and putative cell type were systematically evaluated for immune alterations in different neuropsychiatric disorders. Co-expression networks revealed that transcripts of the neuroimmune systems interacted with neuronal-systems, both of which contribute to the pathology of brain disorders. However, only a few genes with expression changes were also identified as containing risk variants in genome-wide association studies. The transcriptome alterations at gene and network levels may clarify the immune-related pathophysiology and help to better define neuropsychiatric and neurological disorders.

Multi-omics data integration methods and their applications in psychiatric disorders

To study mental illness and health, in the past researchers have often broken down their complexity into individual subsystems (e.g., genomics, transcriptomics, proteomics, clinical data) and explored the components independently. Technological advancements and decreasing costs of high throughput sequencing has led to an unprecedented increase in data generation. Furthermore, over the years it has become increasingly clear that these subsystems do not act in isolation but instead interact with each other to drive mental illness and health. Consequently, individual subsystems are now analysed jointly to promote a holistic understanding of the underlying biological complexity of health and disease. Complementing the increasing data availability, current research is geared towards developing novel methods that can efficiently combine the information rich multi-omics data to discover biologically meaningful biomarkers for diagnosis, treatment, and prognosis. However, clinical translation of the research is still challenging. In this review, we summarise conventional and state-of-the-art statistical and machine learning approaches for discovery of biomarker, diagnosis, as well as outcome and treatment response prediction through integrating multi-omics and clinical data. In addition, we describe the role of biological model systems and in silico multi-omics model designs in clinical translation of psychiatric research from bench to bedside. Finally, we discuss the current challenges and explore the application of multi-omics integration in future psychiatric research. The review provides a structured overview and latest updates in the field of multi-omics in psychiatry.

MicroRNA-eQTLs in the developing human neocortex link miR-4707-3p expression to brain size

Expression quantitative trait loci (eQTL) data have proven important for linking non-coding loci to protein-coding genes. But eQTL studies rarely measure microRNAs (miRNAs), small non-coding RNAs known to play a role in human brain development and neurogenesis. Here, we performed small-RNA sequencing across 212 mid-gestation human neocortical tissue samples, measured 907 expressed miRNAs, discovering 111 of which were novel, and identified 85 local-miRNA-eQTLs. Colocalization of miRNA-eQTLs with GWAS summary statistics yielded one robust colocalization of miR-4707-3p expression with educational attainment and brain size phenotypes, where the miRNA expression increasing allele was associated with decreased brain size. Exogenous expression of miR-4707-3p in primary human neural progenitor cells decreased expression of predicted targets and increased cell proliferation, indicating miR-4707-3p modulates progenitor gene regulation and cell fate decisions. Integrating miRNA-eQTLs with existing GWAS yielded evidence of a miRNA that may influence human brain size and function via modulation of neocortical brain development.

Integrative Post-Genome-Wide Association Study Analyses Relevant to Psychiatric Disorders: Imputing Transcriptome and Proteome Signals

Background

The genome-wide association study (GWAS) is a common tool to identify genetic variants associated with complex traits, including psychiatric disorders (PDs). However, post-GWAS analyses are needed to extend the statistical inference to biologically relevant entities, e.g., genes, proteins, and pathways. To achieve this goal, researchers developed methods that incorporate biologically relevant intermediate molecular phenotypes, such as gene expression and protein abundance, which are posited to mediate the variant-trait association. Transcriptome-wide association study (TWAS) and proteome-wide association study (PWAS) are commonly used methods to test the association between these molecular mediators and the trait.

Summary

In this review, we discuss the most recent developments in TWAS and PWAS. These methods integrate existing "omic" information with the GWAS summary statistics for trait(s) of interest. Specifically, they impute transcript/protein data and test the association between imputed gene expression/protein level with phenotype of interest by using (i) GWAS summary statistics and (ii) reference transcriptomic/proteomic/genomic datasets. TWAS and PWAS are suitable as analysis tools for (i) primary association scan and (ii) fine-mapping to identify potentially causal genes for PDs.

Key Messages

As post-GWAS analyses, TWAS and PWAS have the potential to highlight causal genes for PDs. These prioritized genes could indicate targets for the development of novel drug therapies. For researchers attempting such analyses, we recommend Mendelian randomization tools that use GWAS statistics for both trait and reference datasets, e.g., summary Mendelian randomization (SMR). We base our recommendation on (i) being able to use the same tool for both TWAS and PWAS, (ii) not requiring the pre-computed weights (and thus easier to update for larger reference datasets), and (iii) most larger transcriptome reference datasets are publicly available and easy to transform into a compatible format for SMR analysis.

The evolutionary origin of psychosis

Imagination, the driving force of creativity, and primary psychosis are human-specific, since we do not observe behaviors in other species that would convincingly suggest they possess the same traits. Both these traits have been linked to the function of the prefrontal cortex, which is the most evolutionarily novel region of the human brain. A number of evolutionarily novel genetic and epigenetic changes that determine the human brain-specific structure and function have been discovered in recent years. Among them are genomic loci subjected to increased rates of single nucleotide substitutions in humans, called human accelerated regions. These mostly regulatory regions are involved in brain development and sometimes contain genetic variants that confer a risk for schizophrenia. On the other hand, neuroimaging data suggest that mind wandering and related phenomena (as a proxy of imagination) are in many ways similar to rapid eye movement dreaming, a function also present in non-human species. Furthermore, both functions are similar to psychosis in several ways: for example, the same brain areas are activated both in dreams and visual hallucinations. In the present Perspective we hypothesize that imagination is an evolutionary adaptation of dreaming, while primary psychosis results from deficient control by higher-order brain areas over imagination. In the light of this, human accelerated regions might be one of the key drivers in evolution of human imagination and the pathogenesis of psychotic disorders.

Genomic regulatory sequences in the pathogenesis of bipolar disorder

The lifetime prevalence of bipolar disorder is estimated to be about 2%. Epigenetics defines regulatory mechanisms that determine relatively stable patterns of gene expression by controlling all key steps, from DNA to messenger RNA to protein. This Mini Review highlights recent discoveries of modified epigenetic control resulting from genetic variants associated with bipolar disorder in genome-wide association studies. The revealed epigenetic abnormalities implicate gene transcription and post-transcriptional regulation. In the light of these discoveries, the Mini Review focuses on the genes PACS1, MCHR1, DCLK3, HAPLN4, LMAN2L, TMEM258, GNL3, LRRC57, CACNA1C, CACNA1D, and NOVA2 and their potential biological role in the pathogenesis of bipolar disorder. Molecular mechanisms under control of these genes do not translate into a unified picture and substantially more research is needed to fill the gaps in knowledge and to solve current limitations in prognosis and treatment of bipolar disorder. In conclusion, the genetic and functional studies confirm the complex nature of bipolar disorder and indicate future research directions to explore possible targeted treatment options, eventually working toward a personalized approach.

DeepPerVar: a multi-modal deep learning framework for functional interpretation of genetic variants in personal genome

MOTIVATION

Understanding the functional consequence of genetic variants, especially the non-coding ones, is important but particularly challenging. Genome-wide association studies (GWAS) or quantitative trait locus analyses may be subject to limited statistical power and linkage disequilibrium, and thus are less optimal to pinpoint the causal variants. Moreover, most existing machine-learning approaches, which exploit the functional annotations to interpret and prioritize putative causal variants, cannot accommodate the heterogeneity of personal genetic variations and traits in a population study, targeting a specific disease.

RESULTS

By leveraging paired whole-genome sequencing data and epigenetic functional assays in a population study, we propose a multi-modal deep learning framework to predict genome-wide quantitative epigenetic signals by considering both personal genetic variations and traits. The proposed approach can further evaluate the functional consequence of non-coding variants on an individual level by quantifying the allelic difference of predicted epigenetic signals. By applying the approach to the ROSMAP cohort studying Alzheimer's disease (AD), we demonstrate that the proposed approach can accurately predict quantitative genome-wide epigenetic signals and in key genomic regions of AD causal genes, learn canonical motifs reported to regulate gene expression of AD causal genes, improve the partitioning heritability analysis and prioritize putative causal variants in a GWAS risk locus. Finally, we release the proposed deep learning model as a stand-alone Python toolkit and a web server.

AVAILABILITY AND IMPLEMENTATION

https://github.com/lichen-lab/DeepPerVar.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-wide association study of REM sleep behavior disorder identifies polygenic risk and brain expression effects

Rapid-eye movement (REM) sleep behavior disorder (RBD), enactment of dreams during REM sleep, is an early clinical symptom of alpha-synucleinopathies and defines a more severe subtype. The genetic background of RBD and its underlying mechanisms are not well understood. Here, we perform a genome-wide association study of RBD, identifying five RBD risk loci near SNCA, GBA, TMEM175, INPP5F, and SCARB2. Expression analyses highlight SNCA-AS1 and potentially SCARB2 differential expression in different brain regions in RBD, with SNCA-AS1 further supported by colocalization analyses. Polygenic risk score, pathway analysis, and genetic correlations provide further insights into RBD genetics, highlighting RBD as a unique alpha-synucleinopathy subpopulation that will allow future early intervention.

Illuminating links between cis-regulators and trans-acting variants in the human prefrontal cortex

BACKGROUND

Neuropsychiatric disorders afflict a large portion of the global population and constitute a significant source of disability worldwide. Although Genome-wide Association Studies (GWAS) have identified many disorder-associated variants, the underlying regulatory mechanisms linking them to disorders remain elusive, especially those involving distant genomic elements. Expression quantitative trait loci (eQTLs) constitute a powerful means of providing this missing link. However, most eQTL studies in human brains have focused exclusively on cis-eQTLs, which link variants to nearby genes (i.e., those within 1 Mb of a variant). A complete understanding of disease etiology requires a clearer understanding of trans-regulatory mechanisms, which, in turn, entails a detailed analysis of the relationships between variants and expression changes in distant genes.

METHODS

By leveraging large datasets from the PsychENCODE consortium, we conducted a genome-wide survey of trans-eQTLs in the human dorsolateral prefrontal cortex. We also performed colocalization and mediation analyses to identify mediators in trans-regulation and use trans-eQTLs to link GWAS loci to schizophrenia risk genes.

RESULTS

We identified ~80,000 candidate trans-eQTLs (at FDR<0.25) that influence the expression of ~10K target genes (i.e., "trans-eGenes"). We found that many variants associated with these candidate trans-eQTLs overlap with known cis-eQTLs. Moreover, for >60% of these variants (by colocalization), the cis-eQTL's target gene acts as a mediator for the trans-eQTL SNP's effect on the trans-eGene, highlighting examples of cis-mediation as essential for trans-regulation. Furthermore, many of these colocalized variants fall into a discernable pattern wherein cis-eQTL's target is a transcription factor or RNA-binding protein, which, in turn, targets the gene associated with the candidate trans-eQTL. Finally, we show that trans-regulatory mechanisms provide valuable insights into psychiatric disorders: beyond what had been possible using only cis-eQTLs, we link an additional 23 GWAS loci and 90 risk genes (using colocalization between candidate trans-eQTLs and schizophrenia GWAS loci).

CONCLUSIONS

We demonstrate that the transcriptional architecture of the human brain is orchestrated by both cis- and trans-regulatory variants and found that trans-eQTLs provide insights into brain-disease biology.

Efficacy difference of antipsychotics in Alzheimer's disease and schizophrenia: explained with network efficiency and pathway analysis methods

Approximately 50% of Alzheimer's disease (AD) patients will develop psychotic symptoms and these patients will experience severe rapid cognitive decline compared with those without psychosis (AD-P). Currently, no medication has been approved by the Food and Drug Administration for AD with psychosis (AD+P) specifically, although atypical antipsychotics are widely used in clinical practice. These drugs have demonstrated modest efficacy in managing psychosis in individuals with AD, with an increased frequency of adverse events, including excess mortality. We compared the differences between the genetic variations/genes associated with AD+P and schizophrenia from existing Genome-Wide Association Study and differentially expressed genes (DEGs). We also constructed disease-specific protein-protein interaction networks for AD+P and schizophrenia. Network efficiency was then calculated to characterize the topological structures of these two networks. The efficiency of antipsychotics in these two networks was calculated. A weight adjustment based on binding affinity to drug targets was later applied to refine our results, and 2013 and 2123 genes were identified as related to AD+P and schizophrenia, respectively, with only 115 genes shared. Antipsychotics showed a significantly lower efficiency in the AD+P network than in the schizophrenia network (P < 0.001) indicating that antipsychotics may have less impact in AD+P than in schizophrenia. AD+P may be caused by mechanisms distinct from those in schizophrenia which result in a decreased efficacy of antipsychotics in AD+P. In addition, the network analysis methods provided quantitative explanations of the lower efficacy of antipsychotics in AD+P.

Association between the LRP1B and APOE loci in the development of Parkinson's disease dementia

Parkinson's disease is one of the most common age-related neurodegenerative disorders. Although predominantly a motor disorder, cognitive impairment and dementia are important features of Parkinson's disease, particularly in the later stages of the disease. However, the rate of cognitive decline varies among Parkinson's disease patients, and the genetic basis for this heterogeneity is incompletely understood. To explore the genetic factors associated with rate of progression to Parkinson's disease dementia, we performed a genome-wide survival meta-analysis of 3,923 clinically diagnosed Parkinson's disease cases of European ancestry from four longitudinal cohorts. In total, 6.7% of individuals with Parkinson's disease developed dementia during study follow-up, on average 4.4 ± 2.4 years from disease diagnosis. We have identified the APOE ε4 allele as a major risk factor for the conversion to Parkinson's disease dementia [hazards ratio = 2.41 (1.94-3.00), P = 2.32 × 10-15], as well as a new locus within the ApoE and APP receptor LRP1B gene [hazards ratio = 3.23 (2.17-4.81), P = 7.07 × 10-09]. In a candidate gene analysis, GBA variants were also identified to be associated with higher risk of progression to dementia [hazards ratio = 2.02 (1.21-3.32), P = 0.007]. CSF biomarker analysis also implicated the amyloid pathway in Parkinson's disease dementia, with significantly reduced levels of amyloid β42 (P = 0.0012) in Parkinson's disease dementia compared to Parkinson's disease without dementia. These results identify a new candidate gene associated with faster conversion to dementia in Parkinson's disease and suggest that amyloid-targeting therapy may have a role in preventing Parkinson's disease dementia.