Chromatin accessibility during human first-trimester neurodevelopment.

The human brain develops through a tightly organized cascade of patterning events, induced by transcription factor expression and changes in chromatin accessibility. Although gene expression across the developing brain has been described at single-cell resolution1, similar atlases of chromatin accessibility have been primarily focused on the forebrain2-4. Here we describe chromatin accessibility and paired gene expression across the entire developing human brain during the first trimester (6-13 weeks after conception). We defined 135 clusters and used multiomic measurements to link candidate cis-regulatory elements to gene expression. The number of accessible regions increased both with age and along neuronal differentiation. Using a convolutional neural network, we identified putative functional transcription factor-binding sites in enhancers characterizing neuronal subtypes. We applied this model to cis-regulatory elements linked to ESRRB to elucidate its activation mechanism in the Purkinje cell lineage. Finally, by linking disease-associated single nucleotide polymorphisms to cis-regulatory elements, we validated putative pathogenic mechanisms in several diseases and identified midbrain-derived GABAergic neurons as being the most vulnerable to major depressive disorder-related mutations. Our findings provide a more detailed view of key gene regulatory mechanisms underlying the emergence of brain cell types during the first trimester and a comprehensive reference for future studies related to human neurodevelopment.

Estimating Disorder Probability Based on Polygenic Prediction Using the BPC Approach.

Polygenic Scores (PGSs) summarize an individual's genetic propensity for a given trait in a single value, based on SNP effect sizes derived from Genome-Wide Association Study (GWAS) results. Methods have been developed that apply Bayesian approaches to improve the prediction accuracy of PGSs through optimization of estimated effect sizes. While these methods are generally well-calibrated for continuous traits (implying the predicted values are on average equal to the true trait values), they are not well-calibrated for binary disorder traits in ascertained samples. This is a problem because well-calibrated PGSs are needed to reliably compute the absolute disorder probability for an individual to facilitate future clinical implementation. Here we introduce the Bayesian polygenic score Probability Conversion (BPC) approach, which computes an individual's predicted disorder probability using GWAS summary statistics, an existing Bayesian PGS method (e.g. PRScs, SBayesR), the individual's genotype data, and a prior disorder probability. The BPC approach transforms the PGS to its underlying liability scale, computes the variances of the PGS in cases and controls, and applies Bayes' Theorem to compute the absolute disorder probability; it is practical in its application as it does not require a tuning dataset with both genotype and phenotype data. We applied the BPC approach to extensive simulated data and empirical data of nine disorders. The BPC approach yielded well-calibrated results that were consistently better than the results of another recently published approach.

Simplifying causal gene identification in GWAS loci.

Genome-wide association studies (GWAS) help to identify disease-linked genetic variants, but pinpointing the most likely causal genes in GWAS loci remains challenging. Existing GWAS gene prioritization tools are powerful, but often use complex black box models trained on datasets containing unaddressed biases. Here we present CALDERA, a gene prioritization tool that achieves similar or better performance than state-of-the-art methods, but uses just 12 features and a simple logistic regression model with L1 regularization. We use a data-driven approach to construct a truth set of causal genes in 406 GWAS loci and correct for potential confounders. We demonstrate that CALDERA is well-calibrated in external datasets and prioritizes genes with expected properties, such as being mutation-intolerant (OR = 1.751 for pLI > 90%, P = 8.45×10-3). CALDERA facilitates the prioritization of potentially causal genes in GWAS loci and may help identify novel genetics-driven drug targets.

Evaluation of the Psychometric Properties of the Social Communication Questionnaire in Rural Kenya.

Children can be reliably diagnosed with autism as early as 3 years of age, and early interventions are initiated. There is often a significant gap between the age of onset of symptoms (2-3 years) and diagnosis (8-10 years) in Africa. We conducted a study to validate the Social Communication Questionnaire (SCQ) as a screening instrument in a rural setting in Kenya. The study was conducted along the Kenyan Coast. Study participants included 172 children with a neurodevelopmental disorder (NDD) diagnosis (84 of which were autism) and 112 controls. Internal consistency was evaluated through the use of Cronbach's alpha, confirmatory factor analysis (CFA) with maximum likelihood procedure to assess the conceptual model for the SCQ. Additionally, the sensitivity and specificity of cut-off scores using ROC analysis and item difficulties and discrimination quality using an IRT framework were also assessed. Factor analysis revealed an adequate fitting model for the three-factor DSM-IV-TR (root mean squared error of approximation (RMSEA) = 0.050; Comparative Fit Index (CFI) = 0.974; Tucker-Lewis Index (TLI) = 0.973) and two-factor DSM-5 factor structure (RMSEA = 0.050; CFI = 0.972; TLI = 0.974). The reliability coefficient alphas for the whole group for all items (Cronbach's α = 0.90) and all three domains (Cronbach's α = 0.68-0.84) were acceptable to excellent. The recommended cut-off score of 15 yielded 72% sensitivity and 100% specificity in the ASD group compared to the typically developing group. We provide early evidence of the adequate factor structure and good internal consistency of the SCQ. We also note that the recommended cut-off yielded sufficient predictive validity.

Identifying drug targets for schizophrenia through gene prioritization.

Background: Schizophrenia genome-wide association studies (GWASes) have identified >250 significant loci and prioritized >100 disease-related genes. However, gene prioritization efforts have mostly been restricted to locus-based methods that ignore information from the rest of the genome. Methods: To more accurately characterize genes involved in schizophrenia etiology, we applied a combination of highly-predictive tools to a published GWAS of 67,390 schizophrenia cases and 94,015 controls. We combined both locus-based methods (fine-mapped coding variants, distance to GWAS signals) and genome-wide methods (PoPS, MAGMA, ultra-rare coding variant burden tests). To validate our findings, we compared them with previous prioritization efforts, known neurodevelopmental genes, and results from the PsyOPS tool. Results: We prioritized 62 schizophrenia genes, 41 of which were also highlighted by our validation methods. In addition to DRD2, the principal target of antipsychotics, we prioritized 9 genes that are targeted by approved or investigational drugs. These included drugs targeting glutamatergic receptors (GRIN2A and GRM3), calcium channels (CACNA1C and CACNB2), and GABAB receptor (GABBR2). These also included genes in loci that are shared with an addiction GWAS (e.g. PDE4B and VRK2). Conclusions: We curated a high-quality list of 62 genes that likely play a role in the development of schizophrenia. Developing or repurposing drugs that target these genes may lead to a new generation of schizophrenia therapies. Rodent models of addiction more closely resemble the human disorder than rodent models of schizophrenia. As such, genes prioritized for both disorders could be explored in rodent addiction models, potentially facilitating drug development.

Improved functional mapping of complex trait heritability with GSA-MiXeR implicates biologically specific gene sets.

While genome-wide association studies are increasingly successful in discovering genomic loci associated with complex human traits and disorders, the biological interpretation of these findings remains challenging. Here we developed the GSA-MiXeR analytical tool for gene set analysis (GSA), which fits a model for the heritability of individual genes, accounting for linkage disequilibrium across variants and allowing the quantification of partitioned heritability and fold enrichment for small gene sets. We validated the method using extensive simulations and sensitivity analyses. When applied to a diverse selection of complex traits and disorders, including schizophrenia, GSA-MiXeR prioritizes gene sets with greater biological specificity compared to standard GSA approaches, implicating voltage-gated calcium channel function and dopaminergic signaling for schizophrenia. Such biologically relevant gene sets, often with fewer than ten genes, are more likely to provide insights into the pathobiology of complex diseases and highlight potential drug targets.

Distinguishing different psychiatric disorders using DDx-PRS.

Despite great progress on methods for case-control polygenic prediction (e.g. schizophrenia vs. control), there remains an unmet need for a method that genetically distinguishes clinically related disorders (e.g. schizophrenia (SCZ) vs. bipolar disorder (BIP) vs. depression (MDD) vs. control); such a method could have important clinical value, especially at disorder onset when differential diagnosis can be challenging. Here, we introduce a method, Differential Diagnosis-Polygenic Risk Score (DDx-PRS), that jointly estimates posterior probabilities of each possible diagnostic category (e.g. SCZ=50%, BIP=25%, MDD=15%, control=10%) by modeling variance/covariance structure across disorders, leveraging case-control polygenic risk scores (PRS) for each disorder (computed using existing methods) and prior clinical probabilities for each diagnostic category. DDx-PRS uses only summary-level training data and does not use tuning data, facilitating implementation in clinical settings. In simulations, DDx-PRS was well-calibrated (whereas a simpler approach that analyzes each disorder marginally was poorly calibrated), and effective in distinguishing each diagnostic category vs. the rest. We then applied DDx-PRS to Psychiatric Genomics Consortium SCZ/BIP/MDD/control data, including summary-level training data from 3 case-control GWAS ( N =41,917-173,140 cases; total N =1,048,683) and held-out test data from different cohorts with equal numbers of each diagnostic category (total N =11,460). DDx-PRS was well-calibrated and well-powered relative to these training sample sizes, attaining AUCs of 0.66 for SCZ vs. rest, 0.64 for BIP vs. rest, 0.59 for MDD vs. rest, and 0.68 for control vs. rest. DDx-PRS produced comparable results to methods that leverage tuning data, confirming that DDx-PRS is an effective method. True diagnosis probabilities in top deciles of predicted diagnosis probabilities were considerably larger than prior baseline probabilities, particularly in projections to larger training sample sizes, implying considerable potential for clinical utility under certain circumstances. In conclusion, DDx-PRS is an effective method for distinguishing clinically related disorders.

How Real-World Data Can Facilitate the Development of Precision Medicine Treatment in Psychiatry.

Precision medicine has the ambition to improve treatment response and clinical outcomes through patient stratification and holds great potential for the treatment of mental disorders. However, several important factors are needed to transform current practice into a precision psychiatry framework. Most important are 1) the generation of accessible large real-world training and test data including genomic data integrated from multiple sources, 2) the development and validation of advanced analytical tools for stratification and prediction, and 3) the development of clinically useful management platforms for patient monitoring that can be integrated into health care systems in real-life settings. This narrative review summarizes strategies for obtaining the key elements-well-powered samples from large biobanks integrated with electronic health records and health registry data using novel artificial intelligence algorithms-to predict outcomes in severe mental disorders and translate these models into clinical management and treatment approaches. Key elements are massive mental health data and novel artificial intelligence algorithms. For the clinical translation of these strategies, we discuss a precision medicine platform for improved management of mental disorders. We use cases to illustrate how precision medicine interventions could be brought into psychiatry to improve the clinical outcomes of mental disorders.