A computational tool (H-MAGMA) for improved prediction of brain-disorder risk genes by incorporating brain chromatin interaction profiles

Genome-wide association studies of coffee intake in UK/US participants of European ancestry uncover cohort-specific genetic associations

Coffee is one of the most widely consumed beverages. We performed a genome-wide association study (GWAS) of coffee intake in US-based 23andMe participants (N = 130,153) and identified 7 significant loci, with many replicating in three multi-ancestral cohorts. We examined genetic correlations and performed a phenome-wide association study across hundreds of biomarkers, health, and lifestyle traits, then compared our results to the largest available GWAS of coffee intake from the UK Biobank (UKB; N = 334,659). We observed consistent positive genetic correlations with substance use and obesity in both cohorts. Other genetic correlations were discrepant, including positive genetic correlations between coffee intake and psychiatric illnesses, pain, and gastrointestinal traits in 23andMe that were absent or negative in the UKB, and genetic correlations with cognition that were negative in 23andMe but positive in the UKB. Phenome-wide association study using polygenic scores of coffee intake derived from 23andMe or UKB summary statistics also revealed consistent associations with increased odds of obesity- and red blood cell-related traits, but all other associations were cohort-specific. Our study shows that the genetics of coffee intake associate with substance use and obesity across cohorts, but also that GWAS performed in different populations could capture cultural differences in the relationship between behavior and genetics.

Human genetics and epigenetics of alcohol use disorder

Alcohol use disorder (AUD) is a prominent contributor to global morbidity and mortality. Its complex etiology involves genetics, epigenetics, and environmental factors. We review progress in understanding the genetics and epigenetics of AUD, summarizing the key findings. Advancements in technology over the decades have elevated research from early candidate gene studies to present-day genome-wide scans, unveiling numerous genetic and epigenetic risk factors for AUD. The latest GWAS on more than one million participants identified more than 100 genetic variants, and the largest epigenome-wide association studies (EWAS) in blood and brain samples have revealed tissue-specific epigenetic changes. Downstream analyses revealed enriched pathways, genetic correlations with other traits, transcriptome-wide association in brain tissues, and drug-gene interactions for AUD. We also discuss limitations and future directions, including increasing the power of GWAS and EWAS studies as well as expanding the diversity of populations included in these analyses. Larger samples, novel technologies, and analytic approaches are essential; these include whole-genome sequencing, multiomics, single-cell sequencing, spatial transcriptomics, deep-learning prediction of variant function, and integrated methods for disease risk prediction.

Resolving the three-dimensional interactome of Human Accelerated Regions during human and chimpanzee neurodevelopment

Human Accelerated Regions (HARs) are highly conserved across species but exhibit a significant excess of human-specific sequence changes, suggesting they may have gained novel functions in human evolution. HARs include transcriptional enhancers with human-specific activity and have been implicated in the evolution of the human brain. However, our understanding of how HARs contributed to uniquely human features of the brain is hindered by a lack of insight into the genes and pathways that HARs regulate. It is unclear whether HARs acted by altering the expression of gene targets conserved between HARs and their chimpanzee orthologs or by gaining new gene targets in human, a mechanism termed enhancer hijacking. We generated a high-resolution map of chromatin interactions for 1,590 HARs and their orthologs in human and chimpanzee neural stem cells (NSCs) to comprehensively identify gene targets in both species. HARs and their chimpanzee orthologs targeted a conserved set of 2,963 genes enriched for neurodevelopmental processes including neurogenesis and synaptic transmission. Changes in HAR enhancer activity were correlated with changes in conserved gene target expression. Conserved targets were enriched among genes differentially expressed between human and chimpanzee NSCs or between human and non-human primate developing and adult brain. Species-specific HAR gene targets did not converge on known biological functions and were not significantly enriched among differentially expressed genes, suggesting that HARs did not alter gene expression via enhancer hijacking. HAR gene targets, including differentially expressed targets, also showed cell type-specific expression patterns in the developing human brain, including outer radial glia, which are hypothesized to contribute to human cortical expansion. Our findings support that HARs influenced human brain evolution by altering the expression of conserved gene targets and provide the means to functionally link HARs with novel human brain features.

Microglial function interacts with the environment to affect sex-specific depression risk

There is a two-fold higher incidence of depression in females compared to men with recent studies suggesting a role for microglia in conferring this sex-dependent depression risk. In this study we investigated the nature of this relation. Using GWAS enrichment, gene-set enrichment analysis and Mendelian randomization, we found minimal evidence for a direct relation between genes functionally related to microglia and sex-dependent genetic risk for depression. We then used expression quantitative trait loci and single nucleus RNA-sequencing resources to generate polygenic scores (PGS) representative of individual variation in microglial function in the adult (UK Biobank; N = 54753-72682) and fetal (ALSPAC; N = 1452) periods. The adult microglial PGS moderated the association between BMI (UK Biobank; beta = 0.001, 95 %CI 0.0009 to 0.003, P = 7.74E-6) and financial insecurity (UK Biobank; beta = 0.001, 95 %CI 0.005 to 0.015, P = 2E-4) with depressive symptoms in females. The fetal microglia PGS moderated the association between maternal prenatal depressive symptoms and offspring depressive symptoms at 24 years in females (ALSPAC; beta = 0.04, 95 %CI 0.004 to 0.07, P = 0.03). We found no evidence for an interaction between the microglial PGS and depression risk factors in males. Our results illustrate a role for microglial function in the conferral of sex-dependent depression risk following exposure to a depression risk factor.

The genetic architecture of biological age in nine human organ systems

Investigating the genetic underpinnings of human aging is essential for unraveling the etiology of and developing actionable therapies for chronic diseases. Here, we characterize the genetic architecture of the biological age gap (BAG; the difference between machine learning-predicted age and chronological age) across nine human organ systems in 377,028 participants of European ancestry from the UK Biobank. The BAGs were computed using cross-validated support vector machines, incorporating imaging, physical traits and physiological measures. We identify 393 genomic loci-BAG pairs (P < 5 × 10-8) linked to the brain, eye, cardiovascular, hepatic, immune, metabolic, musculoskeletal, pulmonary and renal systems. Genetic variants associated with the nine BAGs are predominantly specific to the respective organ system (organ specificity) while exerting pleiotropic links with other organ systems (interorgan cross-talk). We find that genetic correlation between the nine BAGs mirrors their phenotypic correlation. Further, a multiorgan causal network established from two-sample Mendelian randomization and latent causal variance models revealed potential causality between chronic diseases (for example, Alzheimer's disease and diabetes), modifiable lifestyle factors (for example, sleep duration and body weight) and multiple BAGs. Our results illustrate the potential for improving human organ health via a multiorgan network, including lifestyle interventions and drug repurposing strategies.

Genome-wide association studies of lifetime and frequency cannabis use in 131,895 individuals

Cannabis is one of the most widely used drugs globally. Decriminalization of cannabis is further increasing cannabis consumption. We performed genome-wide association studies (GWASs) of lifetime (N=131,895) and frequency (N=73,374) of cannabis use. Lifetime cannabis use GWAS identified two loci, one near CADM2 (rs11922956, p=2.40E-11) and another near GRM3 (rs12673181, p=6.90E-09). Frequency of use GWAS identified one locus near CADM2 (rs4856591, p=8.10E-09; r2 =0.76 with rs11922956). Both traits were heritable and genetically correlated with previous GWASs of lifetime use and cannabis use disorder (CUD), as well as other substance use and cognitive traits. Polygenic scores (PGSs) for lifetime and frequency of cannabis use associated cannabis use phenotypes in AllofUs participants. Phenome-wide association study of lifetime cannabis use PGS in a hospital cohort replicated associations with substance use and mood disorders, and uncovered associations with celiac and infectious diseases. This work demonstrates the value of GWASs of CUD transition risk factors.

Navigating centriolar satellites: the role of PCM1 in cellular and organismal processes

Centriolar satellites are ubiquitous membrane-less organelles that play critical roles in numerous cellular and organismal processes. They were initially discovered through electron microscopy as cytoplasmic granules surrounding centrosomes in vertebrate cells. These structures remained enigmatic until the identification of pericentriolar material 1 protein (PCM1) as their molecular marker, which has enabled their in-depth characterization. Recently, centriolar satellites have come into the spotlight due to their links to developmental and neurodegenerative disorders. This review presents a comprehensive summary of the major advances in centriolar satellite biology, with a focus on studies that investigated their biology associated with the essential scaffolding protein PCM1. We begin by exploring the molecular, cellular, and biochemical properties of centriolar satellites, laying the groundwork for a deeper understanding of their functions and mechanisms at both cellular and organismal levels. We then examine the implications of their dysregulation in various diseases, particularly highlighting their emerging roles in neurodegenerative and developmental disorders, as revealed by organismal models of PCM1. We conclude by discussing the current state of knowledge and posing questions about the adaptable nature of these organelles, thereby setting the stage for future research.

The Genesis of Schizophrenia: An Origin Story

Schizophrenia is routinely referred to as a neurodevelopmental disorder, but the role of brain development in a disorder typically diagnosed during early adult life is enigmatic. The authors revisit the neurodevelopmental model of schizophrenia with genomic insights from the most recent schizophrenia clinical genetic association studies, transcriptomic and epigenomic analyses from human postmortem brain studies, and analyses from cellular models that recapitulate neurodevelopment. Emerging insights into schizophrenia genetic risk continue to converge on brain development, particularly stages of early brain development, that may be perturbed to deviate from a typical, normative course, resulting in schizophrenia clinical symptomatology. As the authors explicate, schizophrenia genetic risk is likely dynamic and context dependent, with effects of genetic risk varying spatiotemporally, across the neurodevelopmental continuum. Optimizing therapeutic strategies for the heterogeneous collective of individuals with schizophrenia may likely be guided by leveraging markers of genetic risk and derivative functional insights, well before the emergence of psychosis. Ultimately, rather than a focus on therapeutic intervention during adolescence or adulthood, principles of prediction and prophylaxis in the pre- and perinatal and neonatal stages may best comport with the biology of schizophrenia to address the early-stage perturbations that alter the normative neurodevelopmental trajectory.

Cortical gene expression architecture links healthy neurodevelopment to the imaging, transcriptomics and genetics of autism and schizophrenia

Human brain organization involves the coordinated expression of thousands of genes. For example, the first principal component (C1) of cortical transcription identifies a hierarchy from sensorimotor to association regions. In this study, optimized processing of the Allen Human Brain Atlas revealed two new components of cortical gene expression architecture, C2 and C3, which are distinctively enriched for neuronal, metabolic and immune processes, specific cell types and cytoarchitectonics, and genetic variants associated with intelligence. Using additional datasets (PsychENCODE, Allen Cell Atlas and BrainSpan), we found that C1-C3 represent generalizable transcriptional programs that are coordinated within cells and differentially phased during fetal and postnatal development. Autism spectrum disorder and schizophrenia were specifically associated with C1/C2 and C3, respectively, across neuroimaging, differential expression and genome-wide association studies. Evidence converged especially in support of C3 as a normative transcriptional program for adolescent brain development, which can lead to atypical supragranular cortical connectivity in people at high genetic risk for schizophrenia.

Multi-ancestry meta-analysis of tobacco use disorder identifies 461 potential risk genes and reveals associations with multiple health outcomes

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviours and although strides have been made using genome-wide association studies to identify risk variants, most variants identified have been for nicotine consumption, rather than TUD. Here we leveraged four US biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records) in 653,790 individuals (495,005 European, 114,420 African American and 44,365 Latin American) and data from UK Biobank (ncombined = 898,680). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviours in children and hundreds of medical outcomes, including HIV infection, heart disease and pain. This work furthers our biological understanding of TUD and establishes electronic health records as a source of phenotypic information for studying the genetics of TUD.

Single-cell multi-cohort dissection of the schizophrenia transcriptome

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

Cross-ancestry atlas of gene, isoform, and splicing regulation in the developing human brain

Neuropsychiatric genome-wide association studies (GWASs), including those for autism spectrum disorder and schizophrenia, show strong enrichment for regulatory elements in the developing brain. However, prioritizing risk genes and mechanisms is challenging without a unified regulatory atlas. Across 672 diverse developing human brains, we identified 15,752 genes harboring gene, isoform, and/or splicing quantitative trait loci, mapping 3739 to cellular contexts. Gene expression heritability drops during development, likely reflecting both increasing cellular heterogeneity and the intrinsic properties of neuronal maturation. Isoform-level regulation, particularly in the second trimester, mediated the largest proportion of GWAS heritability. Through colocalization, we prioritized mechanisms for about 60% of GWAS loci across five disorders, exceeding adult brain findings. Finally, we contextualized results within gene and isoform coexpression networks, revealing the comprehensive landscape of transcriptome regulation in development and disease.

The pivotal role of the X-chromosome in the genetic architecture of the human brain

Genes on the X-chromosome are extensively expressed in the human brain. However, little is known for the X-chromosome's impact on the brain anatomy, microstructure, and functional network. We examined 1,045 complex brain imaging traits from 38,529 participants in the UK Biobank. We unveiled potential autosome-X-chromosome interactions, while proposing an atlas outlining dosage compensation (DC) for brain imaging traits. Through extensive association studies, we identified 72 genome-wide significant trait-locus pairs (including 29 new associations) that share genetic architectures with brain-related disorders, notably schizophrenia. Furthermore, we discovered unique sex-specific associations and assessed variations in genetic effects between sexes. Our research offers critical insights into the X-chromosome's role in the human brain, underscoring its contribution to the differences observed in brain structure and functionality between sexes.

Dopamine signaling enriched striatal gene set predicts striatal dopamine synthesis and physiological activity in vivo

The polygenic architecture of schizophrenia implicates several molecular pathways involved in synaptic function. However, it is unclear how polygenic risk funnels through these pathways to translate into syndromic illness. Using tensor decomposition, we analyze gene co-expression in the caudate nucleus, hippocampus, and dorsolateral prefrontal cortex of post-mortem brain samples from 358 individuals. We identify a set of genes predominantly expressed in the caudate nucleus and associated with both clinical state and genetic risk for schizophrenia that shows dopaminergic selectivity. A higher polygenic risk score for schizophrenia parsed by this set of genes predicts greater dopamine synthesis in the striatum and greater striatal activation during reward anticipation. These results translate dopamine-linked genetic risk variation into in vivo neurochemical and hemodynamic phenotypes in the striatum that have long been implicated in the pathophysiology of schizophrenia.

Dissecting Schizophrenia Biology Using Pleiotropy with Cognitive Genomics

Given the increasingly large number of loci discovered by psychiatric GWAS, specification of the key biological pathways underlying these loci has become a priority for the field. We have previously leveraged the pleiotropic genetic relationships between schizophrenia and two cognitive phenotypes (educational attainment and cognitive task performance) to differentiate two subsets of illness-relevant SNPs: (1) those with "concordant" alleles, which are associated with reduced cognitive ability/education and increased schizophrenia risk; and (2) those with "discordant" alleles linked to reduced educational and/or cognitive levels but lower schizophrenia susceptibility. In the present study, we extend our prior work, utilizing larger input GWAS datasets and a more powerful statistical approach to pleiotropic meta-analysis, the Pleiotropic Locus Exploration and Interpretation using Optimal test (PLEIO). Our pleiotropic meta-analysis of schizophrenia and the two cognitive phenotypes revealed 768 significant loci (159 novel). Among these, 347 loci harbored concordant SNPs, 270 encompassed discordant SNPs, and 151 "dual" loci contained concordant and discordant SNPs. Competitive gene-set analysis using MAGMA related concordant SNP loci with neurodevelopmental pathways (e.g., neurogenesis), whereas discordant loci were associated with mature neuronal synaptic functions. These distinctions were also observed in BrainSpan analysis of temporal enrichment patterns across developmental periods, with concordant loci containing more prenatally expressed genes than discordant loci. Dual loci were enriched for genes related to mRNA translation initiation, representing a novel finding in the schizophrenia literature.

A multi-ancestry genetic study of pain intensity in 598,339 veterans

Chronic pain is a common problem, with more than one-fifth of adult Americans reporting pain daily or on most days. It adversely affects the quality of life and imposes substantial personal and economic costs. Efforts to treat chronic pain using opioids had a central role in precipitating the opioid crisis. Despite an estimated heritability of 25-50%, the genetic architecture of chronic pain is not well-characterized, in part because studies have largely been limited to samples of European ancestry. To help address this knowledge gap, we conducted a cross-ancestry meta-analysis of pain intensity in 598,339 participants in the Million Veteran Program, which identified 126 independent genetic loci, 69 of which are new. Pain intensity was genetically correlated with other pain phenotypes, level of substance use and substance use disorders, other psychiatric traits, education level and cognitive traits. Integration of the genome-wide association studies findings with functional genomics data shows enrichment for putatively causal genes (n = 142) and proteins (n = 14) expressed in brain tissues, specifically in GABAergic neurons. Drug repurposing analysis identified anticonvulsants, β-blockers and calcium-channel blockers, among other drug groups, as having potential analgesic effects. Our results provide insights into key molecular contributors to the experience of pain and highlight attractive drug targets.

Loci on chromosome 20 interact with rs16969968 to influence cigarettes per day in European ancestry individuals

BACKGROUND

The understanding of the molecular genetic contributions to smoking is largely limited to the additive effects of individual single nucleotide polymorphisms (SNPs), but the underlying genetic risk is likely to also include dominance, epistatic, and gene-environment interactions.

METHODS

To begin to address this complexity, we attempted to identify genetic interactions between rs16969968, the most replicated SNP associated with smoking quantity, and all SNPs and genes across the genome.

RESULTS

Using the UK Biobank European subsample, we found one SNP, rs1892967, and two genes, PCNA and TMEM230, that showed a significant genome-wide interaction with rs16969968 for log10 CPD and raw CPD, respectively, in a sample of 116 442 individuals who self-reported currently or previously smoking. We extended these analyses to individuals of South Asian descent and meta-analyzed the combined sample of 117 212 individuals of European and South Asian ancestry. We replicated the gene findings in a meta-analysis of five Finnish samples (N=40 140): FinHealth, FINRISK, Finnish Twin Cohort, GeneRISK, and Health-2000-2011.

CONCLUSIONS

To our knowledge, this represents the first reliable epistatic association between single nucleotide polymorphisms for smoking behaviors and provides a novel direction for possible future functional studies related to this interaction. Furthermore, this work demonstrates the feasibility of these analyses by pooling multiple datasets across various ancestries, which may be applied to other top SNPs for smoking and/or other phenotypes.

Genetic and phenotypic similarity across major psychiatric disorders: a systematic review and quantitative assessment

There is widespread overlap across major psychiatric disorders, and this is the case at different levels of observations, from genetic variants to brain structures and function and to symptoms. However, it remains unknown to what extent these commonalities at different levels of observation map onto each other. Here, we systematically review and compare the degree of similarity between psychiatric disorders at all available levels of observation. We searched PubMed and EMBASE between January 1, 2009 and September 8, 2022. We included original studies comparing at least four of the following five diagnostic groups: Schizophrenia, Bipolar Disorder, Major Depressive Disorder, Autism Spectrum Disorder, and Attention Deficit Hyperactivity Disorder, with measures of similarities between all disorder pairs. Data extraction and synthesis were performed by two independent researchers, following the PRISMA guidelines. As main outcome measure, we assessed the Pearson correlation measuring the degree of similarity across disorders pairs between studies and biological levels of observation. We identified 2975 studies, of which 28 were eligible for analysis, featuring similarity measures based on single-nucleotide polymorphisms, gene-based analyses, gene expression, structural and functional connectivity neuroimaging measures. The majority of correlations (88.6%) across disorders between studies, within and between levels of observation, were positive. To identify a consensus ranking of similarities between disorders, we performed a principal component analysis. Its first dimension explained 51.4% (95% CI: 43.2, 65.4) of the variance in disorder similarities across studies and levels of observation. Based on levels of genetic correlation, we estimated the probability of another psychiatric diagnosis in first-degree relatives and showed that they were systematically lower than those observed in population studies. Our findings highlight that genetic and brain factors may underlie a large proportion, but not all of the diagnostic overlaps observed in the clinic.

Clustering Schizophrenia Genes by Their Temporal Expression Patterns Aids Functional Interpretation

BACKGROUND

Schizophrenia is a highly heritable brain disorder with a typical symptom onset in early adulthood. The 2-hit hypothesis posits that schizophrenia results from differential early neurodevelopment, predisposing an individual, followed by a disruption of later brain maturational processes that trigger the onset of symptoms.

STUDY DESIGN

We applied hierarchical clustering to transcription levels of 345 genes previously linked to schizophrenia, derived from cortical tissue samples from 56 donors across the lifespan. We subsequently calculated clustered-specific polygenic risk scores for 743 individuals with schizophrenia and 743 sex- and age-matched healthy controls.

STUDY RESULTS

Clustering revealed a set of 183 genes that was significantly upregulated prenatally and downregulated postnatally and 162 genes that showed the opposite pattern. The prenatally upregulated set of genes was functionally annotated to fundamental cell cycle processes, while the postnatally upregulated set was associated with the immune system and neuronal communication. We found an interaction between the 2 scores; higher prenatal polygenic risk showed a stronger association with schizophrenia diagnosis at higher levels of postnatal polygenic risk. Importantly, this finding was replicated in an independent clinical cohort of 3233 individuals.

CONCLUSIONS

We provide genetics-based evidence that schizophrenia is shaped by disruptions of separable biological processes acting at distinct phases of neurodevelopment. The modeling of genetic risk factors that moderate each other's effect, informed by the timing of their expression, will aid in a better understanding of the development of schizophrenia.

Massively parallel disruption of enhancers active in human neural stem cells

Changes in gene regulation have been linked to the expansion of the human cerebral cortex and to neurodevelopmental disorders, potentially by altering neural progenitor proliferation. However, the effects of genetic variation within regulatory elements on neural progenitors remain obscure. We use sgRNA-Cas9 screens in human neural stem cells (hNSCs) to disrupt 10,674 genes and 26,385 conserved regions in 2,227 enhancers active in the developing human cortex and determine effects on proliferation. Genes with proliferation phenotypes are associated with neurodevelopmental disorders and show biased expression in specific fetal human brain neural progenitor populations. Although enhancer disruptions overall have weaker effects than gene disruptions, we identify enhancer disruptions that severely alter hNSC self-renewal. Disruptions in human accelerated regions, implicated in human brain evolution, also alter proliferation. Integrating proliferation phenotypes with chromatin interactions reveals regulatory relationships between enhancers and their target genes contributing to neurogenesis and potentially to human cortical evolution.

SATB2 organizes the 3D genome architecture of cognition in cortical neurons

The DNA-binding protein SATB2 is genetically linked to human intelligence. We studied its influence on the three-dimensional (3D) epigenome by mapping chromatin interactions and accessibility in control versus SATB2-deficient cortical neurons. We find that SATB2 affects the chromatin looping between enhancers and promoters of neuronal-activity-regulated genes, thus influencing their expression. It also alters A/B compartments, topologically associating domains, and frequently interacting regions. Genes linked to SATB2-dependent 3D genome changes are implicated in highly specialized neuronal functions and contribute to cognitive ability and risk for neuropsychiatric and neurodevelopmental disorders. Non-coding DNA regions with a SATB2-dependent structure are enriched for common variants associated with educational attainment, intelligence, and schizophrenia. Our data establish SATB2 as a cell-type-specific 3D genome modulator, which operates both independently and in cooperation with CCCTC-binding factor (CTCF) to set up the chromatin landscape of pyramidal neurons for cognitive processes.

Multi-ancestry genome-wide association study of major depression aids locus discovery, fine mapping, gene prioritization and causal inference

Most genome-wide association studies (GWAS) of major depression (MD) have been conducted in samples of European ancestry. Here we report a multi-ancestry GWAS of MD, adding data from 21 cohorts with 88,316 MD cases and 902,757 controls to previously reported data. This analysis used a range of measures to define MD and included samples of African (36% of effective sample size), East Asian (26%) and South Asian (6%) ancestry and Hispanic/Latin American participants (32%). The multi-ancestry GWAS identified 53 significantly associated novel loci. For loci from GWAS in European ancestry samples, fewer than expected were transferable to other ancestry groups. Fine mapping benefited from additional sample diversity. A transcriptome-wide association study identified 205 significantly associated novel genes. These findings suggest that, for MD, increasing ancestral and global diversity in genetic studies may be particularly important to ensure discovery of core genes and inform about transferability of findings.

Advances in the understanding of the pathophysiology of schizophrenia and bipolar disorder through induced pluripotent stem cell models

The pathophysiology of schizophrenia and bipolar disorder involves a complex interaction between genetic and environmental factors that begins in the early stages of neurodevelopment. Recent advancements in the field of induced pluripotent stem cells (iPSCs) offer a promising tool for understanding the neurobiological alterations involved in these disorders and, potentially, for developing new treatment options. In this review, we summarize the results of iPSC-based research on schizophrenia and bipolar disorder, showing disturbances in neurodevelopmental processes, imbalance in glutamatergic-GABAergic transmission and neuromorphological alterations. The limitations of the reviewed literature are also highlighted, particularly the methodological heterogeneity of the studies, the limited number of studies developing iPSC models of both diseases simultaneously, and the lack of in-depth clinical characterization of the included samples. Further studies are needed to advance knowledge on the common and disease-specific pathophysiological features of schizophrenia and bipolar disorder and to promote the development of new treatment options.

Multi-ancestry study of the genetics of problematic alcohol use in over 1 million individuals

Problematic alcohol use (PAU), a trait that combines alcohol use disorder and alcohol-related problems assessed with a questionnaire, is a leading cause of death and morbidity worldwide. Here we conducted a large cross-ancestry meta-analysis of PAU in 1,079,947 individuals (European, N = 903,147; African, N = 122,571; Latin American, N = 38,962; East Asian, N = 13,551; and South Asian, N = 1,716 ancestries). We observed a high degree of cross-ancestral similarity in the genetic architecture of PAU and identified 110 independent risk variants in within- and cross-ancestry analyses. Cross-ancestry fine mapping improved the identification of likely causal variants. Prioritizing genes through gene expression and chromatin interaction in brain tissues identified multiple genes associated with PAU. We identified existing medications for potential pharmacological studies by a computational drug repurposing analysis. Cross-ancestry polygenic risk scores showed better performance of association in independent samples than single-ancestry polygenic risk scores. Genetic correlations between PAU and other traits were observed in multiple ancestries, with other substance use traits having the highest correlations. This study advances our knowledge of the genetic etiology of PAU, and these findings may bring possible clinical applicability of genetics insights-together with neuroscience, biology and data science-closer.

Multi-omics cannot replace sample size in genome-wide association studies

The integration of multi-omics information (e.g., epigenetics and transcriptomics) can be useful for interpreting findings from genome-wide association studies (GWAS). It has been suggested that multi-omics could circumvent or greatly reduce the need to increase GWAS sample sizes for novel variant discovery. We tested whether incorporating multi-omics information in earlier and smaller-sized GWAS boosts true-positive discovery of genes that were later revealed by larger GWAS of the same/similar traits. We applied 10 different analytic approaches to integrating multi-omics data from 12 sources (e.g., Genotype-Tissue Expression project) to test whether earlier and smaller GWAS of 4 brain-related traits (alcohol use disorder/problematic alcohol use, major depression/depression, schizophrenia, and intracranial volume/brain volume) could detect genes that were revealed by a later and larger GWAS. Multi-omics data did not reliably identify novel genes in earlier less-powered GWAS (PPV <0.2; 80% false-positive associations). Machine learning predictions marginally increased the number of identified novel genes, correctly identifying 1-8 additional genes, but only for well-powered early GWAS of highly heritable traits (i.e., intracranial volume and schizophrenia). Although multi-omics, particularly positional mapping (i.e., fastBAT, MAGMA, and H-MAGMA), can help to prioritize genes within genome-wide significant loci (PPVs = 0.5-1.0) and translate them into information about disease biology, it does not reliably increase novel gene discovery in brain-related GWAS. To increase power for discovery of novel genes and loci, increasing sample size is required.

The genetic relationships between brain structure and schizophrenia

Genetic risks for schizophrenia are theoretically mediated by genetic effects on brain structure but it has been unclear which genes are associated with both schizophrenia and cortical phenotypes. We accessed genome-wide association studies (GWAS) of schizophrenia (N = 69,369 cases; 236,642 controls), and of three magnetic resonance imaging (MRI) metrics (surface area, cortical thickness, neurite density index) measured at 180 cortical areas (N = 36,843, UK Biobank). Using Hi-C-coupled MAGMA, 61 genes were significantly associated with both schizophrenia and one or more MRI metrics. Whole genome analysis with partial least squares demonstrated significant genetic covariation between schizophrenia and area or thickness of most cortical regions. Genetic similarity between cortical areas was strongly coupled to their phenotypic covariance, and genetic covariation between schizophrenia and brain phenotypes was strongest in the hubs of structural covariance networks. Pleiotropically associated genes were enriched for neurodevelopmental processes and positionally concentrated in chromosomes 3p21, 17q21 and 11p11. Mendelian randomization analysis indicated that genetically determined variation in a posterior cingulate cortical area could be causal for schizophrenia. Parallel analyses of GWAS on bipolar disorder, Alzheimer's disease and height showed that pleiotropic association with MRI metrics was stronger for schizophrenia compared to other disorders.

Systematic investigation of allelic regulatory activity of schizophrenia-associated common variants

Genome-wide association studies (GWASs) have successfully identified 145 genomic regions that contribute to schizophrenia risk, but linkage disequilibrium makes it challenging to discern causal variants. We performed a massively parallel reporter assay (MPRA) on 5,173 fine-mapped schizophrenia GWAS variants in primary human neural progenitors and identified 439 variants with allelic regulatory effects (MPRA-positive variants). Transcription factor binding had modest predictive power, while fine-map posterior probability, enhancer overlap, and evolutionary conservation failed to predict MPRA-positive variants. Furthermore, 64% of MPRA-positive variants did not exhibit expressive quantitative trait loci signature, suggesting that MPRA could identify yet unexplored variants with regulatory potentials. To predict the combinatorial effect of MPRA-positive variants on gene regulation, we propose an accessibility-by-contact model that combines MPRA-measured allelic activity with neuronal chromatin architecture.

Dopamine and schizophrenia from bench to bedside: Discovery of a striatal co-expression risk gene set that predicts in vivo measures of striatal function

Schizophrenia (SCZ) is characterized by a polygenic risk architecture implicating diverse molecular pathways important for synaptic function. However, how polygenic risk funnels through these pathways to translate into syndromic illness is unanswered. To evaluate biologically meaningful pathways of risk, we used tensor decomposition to characterize gene co-expression in post-mortem brain (of neurotypicals: N=154; patients with SCZ: N=84; and GTEX samples N=120) from caudate nucleus (CN), hippocampus (HP), and dorsolateral prefrontal cortex (DLPFC). We identified a CN-predominant gene set showing dopaminergic selectivity that was enriched for genes associated with clinical state and for genes associated with SCZ risk. Parsing polygenic risk score for SCZ based on this specific gene set (parsed-PRS), we found that greater pathway-specific SCZ risk predicted greater in vivo striatal dopamine synthesis capacity measured by [ 18 F]-FDOPA PET in three independent cohorts of neurotypicals and patients (total N=235) and greater fMRI striatal activation during reward anticipation in two additional independent neurotypical cohorts (total N=141). These results reveal a 'bench to bedside' translation of dopamine-linked genetic risk variation in driving in vivo striatal neurochemical and hemodynamic phenotypes that have long been implicated in the pathophysiology of SCZ.

Multi-ancestry meta-analysis of tobacco use disorder prioritizes novel candidate risk genes and reveals associations with numerous health outcomes

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviors, and although strides have been made using genome-wide association studies (GWAS) to identify risk variants, the majority of variants identified have been for nicotine consumption, rather than TUD. We leveraged five biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records, EHR) in 898,680 individuals (739,895 European, 114,420 African American, 44,365 Latin American). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviors in children, and hundreds of medical outcomes, including HIV infection, heart disease, and pain. This work furthers our biological understanding of TUD and establishes EHR as a source of phenotypic information for studying the genetics of TUD.

Multiple genes in a single GWAS risk locus synergistically mediate aberrant synaptic development and function in human neurons

The mechanistic tie between genome-wide association study (GWAS)-implicated risk variants and disease-relevant cellular phenotypes remains largely unknown. Here, using human induced pluripotent stem cell (hiPSC)-derived neurons as a neurodevelopmental model, we identify multiple schizophrenia (SZ) risk variants that display allele-specific open chromatin (ASoC) and are likely to be functional. Editing the strongest ASoC SNP, rs2027349, near vacuolar protein sorting 45 homolog (VPS45) alters the expression of VPS45, lncRNA AC244033.2, and a distal gene, C1orf54. Notably, the transcriptomic changes in neurons are associated with SZ and other neuropsychiatric disorders. Neurons carrying the risk allele exhibit increased dendritic complexity and hyperactivity. Interestingly, individual/combinatorial gene knockdown shows that these genes alter cellular phenotypes in a non-additive synergistic manner. Our study reveals that multiple genes at a single GWAS risk locus mediate a compound effect on neural function, providing a mechanistic link between a non-coding risk variant and disease-related cellular phenotypes.

Genome-Wide Association Studies of Coffee Intake in UK/US Participants of European Ancestry Uncover Gene-Cohort Influences

Coffee is one of the most widely consumed beverages. We performed a genome-wide association study (GWAS) of coffee intake in US-based 23andMe participants (N=130,153) and identified 7 significant loci, with many replicating in three multi-ancestral cohorts. We examined genetic correlations and performed a phenome-wide association study across thousands of biomarkers and health and lifestyle traits, then compared our results to the largest available GWAS of coffee intake from UK Biobank (UKB; N=334,659). The results of these two GWAS were highly discrepant. We observed positive genetic correlations between coffee intake and psychiatric illnesses, pain, and gastrointestinal traits in 23andMe that were absent or negative in UKB. Genetic correlations with cognition were negative in 23andMe but positive in UKB. The only consistent observations were positive genetic correlations with substance use and obesity. Our study shows that GWAS in different cohorts could capture cultural differences in the relationship between behavior and genetics.

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.

Toward Identification of Functional Sequences and Variants in Noncoding DNA

Understanding the noncoding part of the genome, which encodes gene regulation, is necessary to identify genetic mechanisms of disease and translate findings from genome-wide association studies into actionable results for treatments and personalized care. Here we provide an overview of the computational analysis of noncoding regions, starting from gene-regulatory mechanisms and their representation in data. Deep learning methods, when applied to these data, highlight important regulatory sequence elements and predict the functional effects of genetic variants. These and other algorithms are used to predict damaging sequence variants. Finally, we introduce rare-variant association tests that incorporate functional annotations and predictions in order to increase interpretability and statistical power.

Prioritizing genes associated with brain disorders by leveraging enhancer-promoter interactions in diverse neural cells and tissues

BACKGROUND

Prioritizing genes that underlie complex brain disorders poses a considerable challenge. Despite previous studies have found that they shared symptoms and heterogeneity, it remained difficult to systematically identify the risk genes associated with them.

METHODS

By using the CAGE (Cap Analysis of Gene Expression) read alignment files for 439 human cell and tissue types (including primary cells, tissues and cell lines) from FANTOM5 project, we predicted enhancer-promoter interactions (EPIs) of 439 cell and tissue types in human, and examined their reliability. Then we evaluated the genetic heritability of 17 diverse brain disorders and behavioral-cognitive phenotypes in each neural cell type, brain region, and developmental stage. Furthermore, we prioritized genes associated with brain disorders and phenotypes by leveraging the EPIs in each neural cell and tissue type, and analyzed their pleiotropy and functionality for different categories of disorders and phenotypes. Finally, we characterized the spatiotemporal expression dynamics of these associated genes in cells and tissues.

RESULTS

We found that identified EPIs showed activity specificity and network aggregation in cell and tissue types, and enriched TF binding in neural cells played key roles in synaptic plasticity and nerve cell development, i.e., EGR1 and SOX family. We also discovered that most neurological disorders exhibit heritability enrichment in neural stem cells and astrocytes, while psychiatric disorders and behavioral-cognitive phenotypes exhibit enrichment in neurons. Furthermore, our identified genes recapitulated well-known risk genes, which exhibited widespread pleiotropy between psychiatric disorders and behavioral-cognitive phenotypes (i.e., FOXP2), and indicated expression specificity in neural cell types, brain regions, and developmental stages associated with disorders and phenotypes. Importantly, we showed the potential associations of brain disorders with brain regions and developmental stages that have not been well studied.

CONCLUSIONS

Overall, our study characterized the gene-enhancer regulatory networks and genetic mechanisms in the human neural cells and tissues, and illustrated the value of reanalysis of publicly available genomic datasets.

Neurogenetic and multi-omic sources of overlap among sensation seeking, alcohol consumption, and alcohol use disorder

Sensation seeking is bidirectionally associated with levels of alcohol consumption in both adult and adolescent samples and shared neurobiological and genetic influences may in part explain this association. Links between sensation seeking and alcohol use disorder (AUD) may primarily manifest via increased alcohol consumption rather than through direct effects on increasing problems and consequences. Here the overlap between sensation seeking, alcohol consumption, and AUD was examined using multivariate modeling approaches for genome-wide association study (GWAS) summary statistics in conjunction with neurobiologically-informed analyses at multiple levels of investigation. Meta-analytic and genomic structural equation modeling (GenomicSEM) approaches were used to conduct GWAS of sensation seeking, alcohol consumption, and AUD. Resulting summary statistics were used in downstream analyses to examine shared brain tissue enrichment of heritability and genome-wide evidence of overlap (e.g., stratified GenomicSEM, RRHO, genetic correlations with neuroimaging phenotypes) and to identify genomic regions likely contributing to observed genetic overlap across traits (e.g., HMAGMA, LAVA). Across approaches, results supported shared neurogenetic architecture between sensation seeking and alcohol consumption characterized by overlapping enrichment of genes expressed in midbrain and striatal tissues and variants associated with increased cortical surface area. Alcohol consumption and AUD evidenced overlap in relation to variants associated with decreased frontocortical thickness. Finally, genetic mediation models provided evidence of alcohol consumption mediating associations between sensation seeking and AUD. This study extends previous research by examining critical sources of neurogenetic and multi-omic overlap among sensation seeking, alcohol consumption, and AUD which may underlie observed phenotypic associations.

Polygenic influences on the behavioral effects of alcohol withdrawal in a mixed-ancestry population from the collaborative study on the genetics of alcoholism (COGA)

Alcohol withdrawal (AW) is a feature of alcohol use disorder that may occur in up to half of individuals with chronic, heavy alcohol consumption whenever alcohol use is abruptly stopped or significantly reduced. To date, few genes have been robustly associated with AW; this may be partly due to most studies defining AW as a binary construct despite the multiple symptoms and their range in severity from mild to severe. The current study examined the effects of genome-wide loci on a factor score for AW in high risk and community family samples in the Collaborative Study for the Genetics of Alcoholism (COGA). In addition, we tested whether differentially expressed genes associated with alcohol withdrawal in model organisms are enriched in human genome-wide association study (GWAS) effects. Analyses employed roughly equal numbers of males and females (mean age 35, standard deviation = 15; total N = 8009) and included individuals from multiple ancestral backgrounds. Genomic data were imputed to the HRC reference panel and underwent strict quality control procedures using Plink2. Analyses controlled for age, sex, and population stratification effects using ancestral principal components. We found support that AW is a polygenic disease (SNP-heritability = 0.08 [95 % CI = 0.01, 0.15; pedigree-based heritability = 0.12 [0.08,0.16]. We identified five single nucleotide variants that met genomewide significance, some of which have previously been associated with alcohol phenotypes. Gene-level analyses suggest a role for COL19A1 in AW; H-MAGMA analyses implicated 12 genes associated with AW. Cross-species enrichment analyses indicated that variation within genes identified in model organism studies explained <1 % of the phenotypic variability in human AW. Notably, the surrounding regulatory regions of model organism genes explained more variance than expected by chance, indicating that these regulatory regions and gene sets may be important for human AW. Lastly, when comparing the overlap in genes identified from the human GWAS and H-MAGMA analyses with the genes identified from the animal studies, there was modest overlap, indicating some convergence between the methods and organisms.

Unrevealing the shared genetic mechanisms underlying C-reactive protein and schizophrenia

Longitudinal observational studies and Mendelian randomization research have obtained contradictory conclusions regarding the association between C-reactive protein (CRP) level and schizophrenia risk. However, the shared genetic mechanisms underlying CRP and schizophrenia remain poorly understood. Here, we examined the global and local genetic correlations using summary statistics from large-scale genome-wide association studies (GWAS) on CRP level and schizophrenia. Furthermore, we identified their shared genetic variants by applying the conditional false discovery rate approach and performed functional analyses of shared variants to explore the shared genetic mechanisms underlying CRP level and schizophrenia. We found a significant negative genetic correlation at the whole genome level and five significant local genetic correlations between CRP level and schizophrenia. Eight-three shared genetic loci were identified, from which single-nucleotide polymorphism (SNP) presents mixed effects on the increased CRP level and schizophrenia risk. Additionally, we identified 64 and 73 candidate genes that were mapped from SNPs with"concordant effect"(ceSNPs) and"discordant effect"(deSNPs) on the CRP level and schizophrenia risk respectively. Functional analyses revealed that genes mapped from ceSNPs and deSNPs exhibited similar patterns of human brain developmental expression trajectories and biological processes, but differed in expression levels and cell-type-specific enrichment in brain tissues. Our findings demonstrated mixed effects of shared genetic architecture between CRP level and schizophrenia, proving a deeper insight into the shared genetic aetiology underlying the CRP level and schizophrenia.

The functional and evolutionary impacts of human-specific deletions in conserved elements

Conserved genomic sequences disrupted in humans may underlie uniquely human phenotypic traits. We identified and characterized 10,032 human-specific conserved deletions (hCONDELs). These short (average 2.56 base pairs) deletions are enriched for human brain functions across genetic, epigenomic, and transcriptomic datasets. Using massively parallel reporter assays in six cell types, we discovered 800 hCONDELs conferring significant differences in regulatory activity, half of which enhance rather than disrupt regulatory function. We highlight several hCONDELs with putative human-specific effects on brain development, including HDAC5, CPEB4, and PPP2CA. Reverting an hCONDEL to the ancestral sequence alters the expression of LOXL2 and developmental genes involved in myelination and synaptic function. Our data provide a rich resource to investigate the evolutionary mechanisms driving new traits in humans and other species.

Consensus molecular environment of schizophrenia risk genes in coexpression networks shifting across age and brain regions

Schizophrenia is a neurodevelopmental brain disorder whose genetic risk is associated with shifting clinical phenomena across the life span. We investigated the convergence of putative schizophrenia risk genes in brain coexpression networks in postmortem human prefrontal cortex (DLPFC), hippocampus, caudate nucleus, and dentate gyrus granule cells, parsed by specific age periods (total N = 833). The results support an early prefrontal involvement in the biology underlying schizophrenia and reveal a dynamic interplay of regions in which age parsing explains more variance in schizophrenia risk compared to lumping all age periods together. Across multiple data sources and publications, we identify 28 genes that are the most consistently found partners in modules enriched for schizophrenia risk genes in DLPFC; twenty-three are previously unidentified associations with schizophrenia. In iPSC-derived neurons, the relationship of these genes with schizophrenia risk genes is maintained. The genetic architecture of schizophrenia is embedded in shifting coexpression patterns across brain regions and time, potentially underwriting its shifting clinical presentation.

Examining the shared etiology of psychopathology with genome-wide association studies

Genome-wide association studies (GWASs) have ushered in a new era of reproducible discovery in psychiatric genetics. The field has now identified hundreds of common genetic variants that are associated with mental disorders, and many of them influence more than one disorder. By advancing the understanding of causal biology underlying psychopathology, GWAS results are poised to inform the development of novel therapeutics, stratification of at-risk patients, and perhaps even the revision of top-down classification systems in psychiatry. Here, we provide a concise review of GWAS findings with an emphasis on findings that have elucidated the shared genetic etiology of psychopathology, summarizing insights at three levels of analysis: 1) genome-wide architecture; 2) networks, pathways, and gene sets; and 3) individual variants/genes. Three themes emerge from these efforts. First, all psychiatric phenotypes are heritable, highly polygenic, and influenced by many pleiotropic variants with incomplete penetrance. Second, GWAS results highlight the broad etiological roles of neuronal biology, system-wide effects over localized effects, and early neurodevelopment as a critical period. Third, many loci that are robustly associated with multiple forms of psychopathology harbor genes that are involved in synaptic structure and function. Finally, we conclude our review by discussing the implications that GWAS results hold for the field of psychiatry, as well as expected challenges and future directions in the next stage of psychiatric genetics.

Development and Application of Brain Region-Specific Organoids for Investigating Psychiatric Disorders

Human society has been burdened by psychiatric disorders throughout the course of its history. The emergence and rapid advances of human brain organoid technology provide unprecedented opportunities for investigation of potential disease mechanisms and development of targeted or even personalized treatments for various psychiatric disorders. In this review, we summarize recent advances for generating organoids from human pluripotent stem cells to model distinct brain regions and diverse cell types. We also highlight recent progress, discuss limitations, and propose potential improvements in using patient-derived or genetically engineered brain region-specific organoids for investigating various psychiatric disorders.

Integration of epigenetic and genetic profiles identifies multiple sclerosis disease-critical cell types and genes

Genome-wide association studies (GWAS) successfully identified multiple sclerosis (MS) susceptibility variants. Despite this notable progress, understanding the biological context of these associations remains challenging, due in part to the complexity of linking GWAS results to causative genes and cell types. Here, we aimed to address this gap by integrating GWAS data with single-cell and bulk chromatin accessibility data and histone modification profiles from immune and nervous systems. MS-GWAS associations are significantly enriched in regulatory regions of microglia and peripheral immune cell subtypes, especially B cells and monocytes. Cell-specific polygenic risk scores were developed to examine the cumulative impact of the susceptibility genes on MS risk and clinical phenotypes, showing significant associations with risk and brain white matter volume. The findings reveal enrichment of GWAS signals in B cell and monocyte/microglial cell-types, consistent with the known pathology and presumed targets of effective MS therapeutics.

The genetic architecture of pain intensity in a sample of 598,339 U.S. veterans

Chronic pain is a common problem, with more than one-fifth of adult Americans reporting pain daily or on most days. It adversely affects quality of life and imposes substantial personal and economic costs. Efforts to treat chronic pain using opioids played a central role in precipitating the opioid crisis. Despite an estimated heritability of 25-50%, the genetic architecture of chronic pain is not well characterized, in part because studies have largely been limited to samples of European ancestry. To help address this knowledge gap, we conducted a cross-ancestry meta-analysis of pain intensity in 598,339 participants in the Million Veteran Program, which identified 125 independent genetic loci, 82 of which are novel. Pain intensity was genetically correlated with other pain phenotypes, level of substance use and substance use disorders, other psychiatric traits, education level, and cognitive traits. Integration of the GWAS findings with functional genomics data shows enrichment for putatively causal genes (n = 142) and proteins (n = 14) expressed in brain tissues, specifically in GABAergic neurons. Drug repurposing analysis identified anticonvulsants, beta-blockers, and calcium-channel blockers, among other drug groups, as having potential analgesic effects. Our results provide insights into key molecular contributors to the experience of pain and highlight attractive drug targets.

Cross-ancestry, cell-type-informed atlas of gene, isoform, and splicing regulation in the developing human brain

Genomic regulatory elements active in the developing human brain are notably enriched in genetic risk for neuropsychiatric disorders, including autism spectrum disorder (ASD), schizophrenia, and bipolar disorder. However, prioritizing the specific risk genes and candidate molecular mechanisms underlying these genetic enrichments has been hindered by the lack of a single unified large-scale gene regulatory atlas of human brain development. Here, we uniformly process and systematically characterize gene, isoform, and splicing quantitative trait loci (xQTLs) in 672 fetal brain samples from unique subjects across multiple ancestral populations. We identify 15,752 genes harboring a significant xQTL and map 3,739 eQTLs to a specific cellular context. We observe a striking drop in gene expression and splicing heritability as the human brain develops. Isoform-level regulation, particularly in the second trimester, mediates the greatest proportion of heritability across multiple psychiatric GWAS, compared with eQTLs. Via colocalization and TWAS, we prioritize biological mechanisms for ~60% of GWAS loci across five neuropsychiatric disorders, nearly two-fold that observed in the adult brain. Finally, we build a comprehensive set of developmentally regulated gene and isoform co-expression networks capturing unique genetic enrichments across disorders. Together, this work provides a comprehensive view of genetic regulation across human brain development as well as the stage-and cell type-informed mechanistic underpinnings of neuropsychiatric disorders.

Pathway-specific polygenic risk scores correlate with clinical status and Alzheimer's-related biomarkers

Background: APOE is the largest genetic risk factor for sporadic Alzheimer's disease (AD), but there is a substantial polygenic component as well. Polygenic risk scores (PRS) can summarize small effects across the genome but may obscure differential risk associated with different molecular processes and pathways. Variability at the genetic level may contribute to the extensive phenotypic heterogeneity of Alzheimer's disease (AD). Here, we examine polygenic risk impacting specific pathways associated with AD and examined its relationship with clinical status and AD biomarkers of amyloid, tau, and neurodegeneration (A/T/N). Methods: A total of 1,411 participants from the Alzheimer's Disease Neuroimaging Initiative (ADNI) with genotyping data were included. Sets of variants identified from a pathway analysis of AD GWAS summary statistics were combined into clusters based on their assigned pathway. We constructed pathway-specific PRSs for each participant and tested their associations with diagnostic status (AD vs cognitively normal), abnormal levels of amyloid and ptau (positive vs negative), and hippocampal volume. The APOE region was excluded from all PRSs, and analyses controlled for APOE -ε4 carrier status. Results: Thirteen pathway clusters were identified relating to categories such as immune response, amyloid precursor processing, protein localization, lipid transport and binding, tyrosine kinase, and endocytosis. Eight pathway-specific PRSs were significantly associated with AD dementia diagnosis. Amyloid-positivity was associated with endocytosis and fibril formation, response misfolded protein, and regulation protein tyrosine PRSs. Ptau positivity and hippocampal volume were both related to protein localization and mitophagy PRS, and ptau positivity was additionally associated with an immune signaling PRS. A global AD PRS showed stronger associations with diagnosis and all biomarkers compared to pathway PRSs, suggesting a strong synergistic effect of all loci contributing to the global AD PRS. Conclusions: Pathway PRS may contribute to understanding separable disease processes, but do not appear to add significant power for predictive purposes. These findings demonstrate that, although genetic risk for AD is widely distributed, AD-phenotypes may be preferentially associated with risk in specific pathways. Defining genetic risk along multiple dimensions at the individual level may help clarify the etiological heterogeneity in AD.

Examining the biological mechanisms of human mental disorders resulting from gene-environment interdependence using novel functional genomic approaches

We explore how functional genomics approaches that integrate datasets from human and non-human model systems can improve our understanding of the effect of gene-environment interplay on the risk for mental disorders. We start by briefly defining the G-E paradigm and its challenges and then discuss the different levels of regulation of gene expression and the corresponding data existing in humans (genome wide genotyping, transcriptomics, DNA methylation, chromatin modifications, chromosome conformational changes, non-coding RNAs, proteomics and metabolomics), discussing novel approaches to the application of these data in the study of the origins of mental health. Finally, we discuss the multilevel integration of diverse types of data. Advance in the use of functional genomics in the context of a G-E perspective improves the detection of vulnerabilities, informing the development of preventive and therapeutic interventions.

BIGKnock: fine-mapping gene-based associations via knockoff analysis of biobank-scale data

We propose BIGKnock (BIobank-scale Gene-based association test via Knockoffs), a computationally efficient gene-based testing approach for biobank-scale data, that leverages long-range chromatin interaction data, and performs conditional genome-wide testing via knockoffs. BIGKnock can prioritize causal genes over proxy associations at a locus. We apply BIGKnock to the UK Biobank data with 405,296 participants for multiple binary and quantitative traits, and show that relative to conventional gene-based tests, BIGKnock produces smaller sets of significant genes that contain the causal gene(s) with high probability. We further illustrate its ability to pinpoint potential causal genes at [Formula: see text] of the associated loci.

Role of the Gut-Brain Axis in the Shared Genetic Etiology Between Gastrointestinal Tract Diseases and Psychiatric Disorders: A Genome-Wide Pleiotropic Analysis

IMPORTANCE

Comorbidities and genetic correlations between gastrointestinal tract diseases and psychiatric disorders have been widely reported, with the gut-brain axis (GBA) hypothesized as a potential biological basis. However, the degree to which the shared genetic determinants are involved in these associations underlying the GBA is unclear.

OBJECTIVE

To investigate the shared genetic etiology between gastrointestinal tract diseases and psychiatric disorders and to identify shared genomic loci, genes, and pathways.

DESIGN, SETTING, AND PARTICIPANTS

This genome-wide pleiotropic association study using genome-wide association summary statistics from publicly available data sources was performed with various statistical genetic approaches to sequentially investigate the pleiotropic associations from genome-wide single-nucleotide variation (SNV; formerly single-nucleotide polymorphism [SNP]), and gene levels and biological pathways to disentangle the underlying shared genetic etiology between 4 gastrointestinal tract diseases (inflammatory bowel disease, irritable bowel syndrome, peptic ulcer disease, and gastroesophageal reflux disease) and 6 psychiatric disorders (schizophrenia, bipolar disorder, major depressive disorder, attention-deficit/hyperactivity disorder, posttraumatic stress disorder, and anorexia nervosa). Data were collected from March 10, 2021, to August 25, 2021, and analysis was performed from January 8 through May 30, 2022.

MAIN OUTCOMES AND MEASURES

The primary outcomes consisted of a list of genetic loci, genes, and pathways shared between gastrointestinal tract diseases and psychiatric disorders.

RESULTS

Extensive genetic correlations and genetic overlaps were found among 22 of 24 trait pairs. Pleiotropic analysis under a composite null hypothesis identified 2910 significant potential pleiotropic SNVs in 19 trait pairs, with 83 pleiotropic loci and 24 colocalized loci detected. Gene-based analysis found 158 unique candidate pleiotropic genes, which were highly enriched in certain GBA-related phenotypes and tissues, whereas pathway enrichment analysis further highlighted biological pathways primarily involving cell adhesion, synaptic structure and function, and immune cell differentiation. Several identified pleiotropic loci also shared causal variants with gut microbiomes. Mendelian randomization analysis further illustrated vertical pleiotropy across 8 pairwise traits. Notably, many pleiotropic loci were identified for multiple pairwise traits, such as 1q32.1 (INAVA), 19q13.33 (FUT2), 11q23.2 (NCAM1), and 1p32.3 (LRP8).

CONCLUSIONS AND RELEVANCE

These findings suggest that the pleiotropic genetic determinants between gastrointestinal tract diseases and psychiatric disorders are extensively distributed across the genome. These findings not only support the shared genetic basis underlying the GBA but also have important implications for intervention and treatment targets of these diseases simultaneously.

TTLL11 gene is associated with sustained attention performance and brain networks: A genome-wide association study of a healthy Chinese sample

Genetic studies on attention have mainly focused on children with attention-deficit/hyperactivity disorder (ADHD), so little systematic research has been conducted on genetic correlates of attention performance and their potential brain mechanisms among healthy individuals. The current study included a genome-wide association study (GWAS, N = 1145 healthy young adults) aimed to identify genes associated with sustained attention and an imaging genetics study (an independent sample of 483 healthy young adults) to examine any identified genes' influences on brain function. The GWAS found that TTLL11 showed genome-wide significant associations with sustained attention, with rs13298112 as the most significant SNP and the GG homozygotes showing more impulsive but also more focused responses than the A allele carriers. A retrospective examination of previously published ADHD GWAS results confirmed an un-reported, small but statistically significant effect of TTLL11 on ADHD. The imaging genetics study replicated this association and showed that the TTLL11 gene was associated with resting state activity and connectivity of the somatomoter network, and can be predicted by dorsal attention network connectivity. Specifically, the GG homozygotes showed lower brain activity, weaker brain network connectivity, and non-significant brain-attention association compared to the A allele carriers. Expression database showed that expression of this gene is enriched in the brain and that the G allele is associated with lower expression level than the A allele. These results suggest that TTLL11 may play a major role in healthy individuals' attention performance and may also contribute to the etiology of ADHD.

Multi-ancestry study of the genetics of problematic alcohol use in >1 million individuals

Problematic alcohol use (PAU) is a leading cause of death and disability worldwide. To improve our understanding of the genetics of PAU, we conducted a large cross-ancestry meta-analysis of PAU in 1,079,947 individuals. We observed a high degree of cross-ancestral similarity in the genetic architecture of PAU and identified 110 independent risk variants in within- and cross-ancestry analyses. Cross-ancestry fine-mapping improved the identification of likely causal variants. Prioritizing genes through gene expression and/or chromatin interaction in brain tissues identified multiple genes associated with PAU. We identified existing medications for potential pharmacological studies by drug repurposing analysis. Cross-ancestry polygenic risk scores (PRS) showed better performance in independent sample than single-ancestry PRS. Genetic correlations between PAU and other traits were observed in multiple ancestries, with other substance use traits having the highest correlations. The analysis of diverse ancestries contributed significantly to the findings, and fills an important gap in the literature.