Genome-Wide Association Study of Treatment-Resistant Depression: Shared Biology With Metabolic Traits.

OBJECTIVE: Treatment-resistant depression (TRD) occurs in roughly one-third of all individuals with major depressive disorder (MDD). Although research has suggested a significant common variant genetic component of liability to TRD, with heritability estimated at 8% when compared with non-treatment-resistant MDD, no replicated genetic loci have been identified, and the genetic architecture of TRD remains unclear. A key barrier to this work has been the paucity of adequately powered cohorts for investigation, largely because of the challenge in prospectively investigating this phenotype. The objective of this study was to perform a well-powered genetic study of TRD. METHODS: Using receipt of electroconvulsive therapy (ECT) as a surrogate for TRD, the authors applied standard machine learning methods to electronic health record data to derive predicted probabilities of receiving ECT. These probabilities were then applied as a quantitative trait in a genome-wide association study of 154,433 genotyped patients across four large biobanks. RESULTS: Heritability estimates ranged from 2% to 4.2%, and significant genetic overlap was observed with cognition, attention deficit hyperactivity disorder, schizophrenia, alcohol and smoking traits, and body mass index. Two genome-wide significant loci were identified, both previously implicated in metabolic traits, suggesting shared biology and potential pharmacological implications. CONCLUSIONS: This work provides support for the utility of estimation of disease probability for genomic investigation and provides insights into the genetic architecture and biology of TRD.

Polygenic Risk for Type 2 Diabetes in African Americans.

African Americans (AAs) have been underrepresented in polygenic risk score (PRS) studies. Here, we integrated genome-wide data from multiple observational studies on type 2 diabetes (T2D), encompassing a total of 101,987 AAs, to train and optimize an AA-focused T2D PRS (PRSAA), using a Bayesian polygenic modeling method. We further tested the score in three independent studies with a total of 7,275 AAs and compared the PRSAA with other published scores. Results show that a 1-SD increase in the PRSAA was associated with 40-60% increase in the odds of T2D (odds ratio [OR] 1.60, 95% CI 1.37-1.88; OR 1.40, 95% CI 1.16-1.70; and OR 1.45, 95% CI 1.30-1.62) across three testing cohorts. These models captured 1.0-2.6% of the variance (R2) in T2D on the liability scale. The positive predictive values for three calculated score thresholds (the top 2%, 5%, and 10%) ranged from 14 to 35%. The PRSAA, in general, performed similarly to existing T2D PRS. The need remains for larger data sets to continue to evaluate the utility of within-ancestry scores in the AA population.

Whole-exome sequencing in UK Biobank reveals rare genetic architecture for depression.

Nearly two hundred common-variant depression risk loci have been identified by genome-wide association studies (GWAS). However, the impact of rare coding variants on depression remains poorly understood. Here, we present whole-exome sequencing analyses of depression with seven different definitions based on survey, questionnaire, and electronic health records in 320,356 UK Biobank participants. We showed that the burden of rare damaging coding variants in loss-of-function intolerant genes is significantly associated with risk of depression with various definitions. We compared the rare and common genetic architecture across depression definitions by genetic correlation and showed different genetic relationships between definitions across common and rare variants. In addition, we demonstrated that the effects of rare damaging coding variant burden and polygenic risk score on depression risk are additive. The gene set burden analyses revealed overlapping rare genetic variant components with developmental disorder, autism, and schizophrenia. Our study provides insights into the contribution of rare coding variants, separately and in conjunction with common variants, on depression with various definitions and their genetic relationships with neurodevelopmental disorders.

Regional analysis of posterior corneal elevation after laser refractive surgeries for correction of myopia of different degrees.

PURPOSE: To evaluate regional changes in the posterior corneal elevation after three laser refractive surgeries for correction of myopia of different degrees. SETTINGS AND DESIGN: Retrospective, comparative, and non-randomized study. METHODS: Two hundred patients (200 eyes) who underwent laser epithelial keratoplasty (LASEK), femtosecond-assisted laser in-situ keratomileusis (FS-LASIK), and small-incision lenticule extraction (SMILE) were included in this study. According to preoperative spherical equivalent (SE), each surgical group was divided into two refractive subgroups: low-to-moderate myopia (LM group) and high myopia (H group). The posterior corneal elevation from Pentacam Scheimpflug tomography was analyzed preoperatively and at 1 month, 3 months, 6 months, and 12 months postoperatively. Three subregions of the posterior cornea were divided in this study as the central, paracentral, and peripheral regions. STATISTICAL ANALYSIS USED: Generalized Estimating Equations (GEE). RESULTS: For all three surgical groups, similar changing trends were seen in the two refractive subgroups. H group presented a larger changing magnitude than the LM group in FS-LASIK over time (P < 0.05), whereas no significant difference was noted in the two refractive subgroups of LASEK or SMILE (P > 0.05). At 12 months postoperatively, the central posterior corneal elevation returned to the preoperative level in LASEK (P > 0.05) but shifted forward significantly in FS-LASIK and SMILE (P < 0.05). CONCLUSION: Different posterior corneal regions respond differently to corneal refractive surgeries. LASEK, FS-LASIK, and SMILE demonstrate different trends in the regional changes in posterior corneal elevation. The corneal shape seems more stable in LASEK than in FS-LASIK and SMILE.

Selection, optimization and validation of ten chronic disease polygenic risk scores for clinical implementation in diverse US populations.

Polygenic risk scores (PRSs) have improved in predictive performance, but several challenges remain to be addressed before PRSs can be implemented in the clinic, including reduced predictive performance of PRSs in diverse populations, and the interpretation and communication of genetic results to both providers and patients. To address these challenges, the National Human Genome Research Institute-funded Electronic Medical Records and Genomics (eMERGE) Network has developed a framework and pipeline for return of a PRS-based genome-informed risk assessment to 25,000 diverse adults and children as part of a clinical study. From an initial list of 23 conditions, ten were selected for implementation based on PRS performance, medical actionability and potential clinical utility, including cardiometabolic diseases and cancer. Standardized metrics were considered in the selection process, with additional consideration given to strength of evidence in African and Hispanic populations. We then developed a pipeline for clinical PRS implementation (score transfer to a clinical laboratory, validation and verification of score performance), and used genetic ancestry to calibrate PRS mean and variance, utilizing genetically diverse data from 13,475 participants of the All of Us Research Program cohort to train and test model parameters. Finally, we created a framework for regulatory compliance and developed a PRS clinical report for return to providers and for inclusion in an additional genome-informed risk assessment. The initial experience from eMERGE can inform the approach needed to implement PRS-based testing in diverse clinical settings.

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.

Polygenic risk for suicide attempt is associated with lifetime suicide attempt in US soldiers independent of parental risk.

BACKGROUND: Suicide is a leading cause of death worldwide. Whereas some studies have suggested that a direct measure of common genetic liability for suicide attempts (SA), captured by a polygenic risk score for SA (SA-PRS), explains risk independent of parental history, further confirmation would be useful. Even more unsettled is the extent to which SA-PRS is associated with lifetime non-suicidal self-injury (NSSI). METHODS: We used summary statistics from the largest available GWAS study of SA to generate SA-PRS for two non-overlapping cohorts of soldiers of European ancestry. These were tested in multivariable models that included parental major depressive disorder (MDD) and parental SA. RESULTS: In the first cohort, 417 (6.3 %) of 6573 soldiers reported lifetime SA and 1195 (18.2 %) reported lifetime NSSI. In a multivariable model that included parental history of MDD and parental history of SA, SA-PRS remained significantly associated with lifetime SA [aOR = 1.26, 95%CI:1.13-1.39, p < 0.001] per standardized unit SA-PRS]. In the second cohort, 204 (4.2 %) of 4900 soldiers reported lifetime SA, and 299 (6.1 %) reported lifetime NSSI. In a multivariable model that included parental history of MDD and parental history of SA, SA-PRS remained significantly associated with lifetime SA [aOR = 1.20, 95%CI:1.04-1.38, p = 0.014]. A combined analysis of both cohorts yielded similar results. In neither cohort or in the combined analysis was SA-PRS significantly associated with NSSI. CONCLUSIONS: PRS for SA conveys information about likelihood of lifetime SA (but not NSSI, demonstrating specificity), independent of self-reported parental history of MDD and parental history of SA. LIMITATIONS: At present, the magnitude of effects is small and would not be immediately useful for clinical decision-making or risk-stratified prevention initiatives, but this may be expected to improve with further iterations. Also critical will be the extension of these findings to more diverse populations.

Genome-wide Polygenic Risk Scores Predict Risk of Glioma and Molecular Subtypes.

Background: Polygenic risk scores (PRS) aggregate the contribution of many risk variants to provide a personalized genetic susceptibility profile. Since sample sizes of glioma genome-wide association studies (GWAS) remain modest, there is a need to find efficient ways of capturing genetic risk factors using available germline data. Methods: We developed a novel PRS (PRS-CS) that uses continuous shrinkage priors to model the joint effects of over 1 million polymorphisms on disease risk and compared it to an approach (PRS-CT) that selects a limited set of independent variants that reach genome-wide significance (P<5×10-8). PRS models were trained using GWAS results stratified by histological (10,346 cases, 14,687 controls) and molecular subtype (2,632 cases, 2,445 controls), and validated in two independent cohorts. Results: PRS-CS was consistently more predictive than PRS-CT across glioma subtypes with an average increase in explained variance (R2) of 21%. Improvements were particularly pronounced for glioblastoma tumors, with PRS-CS yielding larger effect sizes (odds ratio (OR)=1.93, P=2.0×10-54 vs. OR=1.83, P=9.4×10-50) and higher explained variance (R2=2.82% vs. R2=2.56%). Individuals in the 95th percentile of the PRS-CS distribution had a 3-fold higher lifetime absolute risk of IDH mutant (0.63%) and IDH wildtype (0.76%) glioma relative to individuals with average PRS. PRS-CS also showed high classification accuracy for IDH mutation status among cases (AUC=0.895). Conclusions: Our novel genome-wide PRS may improve the identification of high-risk individuals and help distinguish between prognostic glioma subtypes, increasing the potential clinical utility of germline genetics in glioma patient management.

Shared genetic architectures of educational attainment in East Asian and European populations.

Educational attainment (EduYears), a heritable trait often used as a proxy for cognitive ability, is associated with various health and social outcomes. Previous genome-wide association studies (GWASs) on EduYears have been focused on samples of European (EUR) genetic ancestries. Here we present the first large-scale GWAS of EduYears in people of East Asian (EAS) ancestry (n = 176,400) and conduct a cross-ancestry meta-analysis with EduYears GWAS in people of EUR ancestry (n = 766,345). EduYears showed a high genetic correlation and power-adjusted transferability ratio between EAS and EUR. We also found similar functional enrichment, gene expression enrichment and cross-trait genetic correlations between two populations. Cross-ancestry fine-mapping identified refined credible sets with a higher posterior inclusion probability than single population fine-mapping. Polygenic prediction analysis in four independent EAS and EUR cohorts demonstrated transferability between populations. Our study supports the need for further research on diverse ancestries to increase our understanding of the genetic basis of educational attainment.

Principles and methods for transferring polygenic risk scores across global populations.

Polygenic risk scores (PRSs) summarize the genetic predisposition of a complex human trait or disease and may become a valuable tool for advancing precision medicine. However, PRSs that are developed in populations of predominantly European genetic ancestries can increase health disparities due to poor predictive performance in individuals of diverse and complex genetic ancestries. We describe genetic and modifiable risk factors that limit the transferability of PRSs across populations and review the strengths and weaknesses of existing PRS construction methods for diverse ancestries. Developing PRSs that benefit global populations in research and clinical settings provides an opportunity for innovation and is essential for health equity.

Morphological changes of the limbic system associated with acute and chronic low-back pain: A UK biobank imaging study.

BACKGROUND: Low back pain (LBP) is a major public health issue that influences physical and emotional factors integral to the limbic system. This study aims to investigate the association between LBP and brain morphometry alterations as the duration of LBP increases (acute vs. chronic). METHODS: We used the UK Biobank data to investigate the morphological features of the limbic system in acute LBP (N = 115), chronic LBP (N = 243) and controls (N = 358), and tried to replicate our findings with an independent dataset composed of 45 acute LBP participants evaluated at different timepoints throughout 1 year from the OpenPain database. RESULTS: We found that in comparison with chronic LBP and pain-free controls, acute LBP was associated with increased volumes of the nucleus accumbens, amygdala, hippocampus, and thalamus, and increased grey matter volumes in the hippocampus and posterior cingulate gyrus. In the replication cohort, we found non-significantly larger hippocampus and thalamus volumes in the 3-month visit (acute LBP) compared to the 1-year visit (chronic LBP), with similar effect sizes as the UK Biobank dataset. CONCLUSIONS: Our results suggest that acute LBP is associated with dramatic morphometric increases in the limbic system and mesolimbic pathway, which may reflect an active brain response and self-regulation in the early stage of LBP. SIGNIFICANCE: Our study suggests that LBP in the acute phase is associated with the brain morphometric changes (increase) in some limbic areas, indicating that the acute phase of LBP may represent a crucial stage of self-regulation and active response to the disease's onset.

Beyond the Global Brain Differences: Intraindividual Variability Differences in 1q21.1 Distal and 15q11.2 BP1-BP2 Deletion Carriers.

BACKGROUND: Carriers of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants exhibit regional and global brain differences compared with noncarriers. However, interpreting regional differences is challenging if a global difference drives the regional brain differences. Intraindividual variability measures can be used to test for regional differences beyond global differences in brain structure. METHODS: Magnetic resonance imaging data were used to obtain regional brain values for 1q21.1 distal deletion (n = 30) and duplication (n = 27) and 15q11.2 BP1-BP2 deletion (n = 170) and duplication (n = 243) carriers and matched noncarriers (n = 2350). Regional intra-deviation scores, i.e., the standardized difference between an individual's regional difference and global difference, were used to test for regional differences that diverge from the global difference. RESULTS: For the 1q21.1 distal deletion carriers, cortical surface area for regions in the medial visual cortex, posterior cingulate, and temporal pole differed less and regions in the prefrontal and superior temporal cortex differed more than the global difference in cortical surface area. For the 15q11.2 BP1-BP2 deletion carriers, cortical thickness in regions in the medial visual cortex, auditory cortex, and temporal pole differed less and the prefrontal and somatosensory cortex differed more than the global difference in cortical thickness. CONCLUSIONS: We find evidence for regional effects beyond differences in global brain measures in 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants. The results provide new insight into brain profiling of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants, with the potential to increase understanding of the mechanisms involved in altered neurodevelopment.

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.

Analysis across Taiwan Biobank, Biobank Japan, and UK Biobank identifies hundreds of novel loci for 36 quantitative traits.

Genome-wide association studies (GWASs) have identified tens of thousands of genetic loci associated with human complex traits. However, the majority of GWASs were conducted in individuals of European ancestries. Failure to capture global genetic diversity has limited genomic discovery and has impeded equitable delivery of genomic knowledge to diverse populations. Here we report findings from 102,900 individuals across 36 human quantitative traits in the Taiwan Biobank (TWB), a major biobank effort that broadens the population diversity of genetic studies in East Asia. We identified 968 novel genetic loci, pinpointed novel causal variants through statistical fine-mapping, compared the genetic architecture across TWB, Biobank Japan, and UK Biobank, and evaluated the utility of cross-phenotype, cross-population polygenic risk scores in disease risk prediction. These results demonstrated the potential to advance discovery through diversifying GWAS populations and provided insights into the common genetic basis of human complex traits in East Asia.

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.

Cell-type-specific Alzheimer's disease polygenic risk scores are associated with distinct disease processes in Alzheimer's disease.

Many of the Alzheimer's disease (AD) risk genes are specifically expressed in microglia and astrocytes, but how and when the genetic risk localizing to these cell types contributes to AD pathophysiology remains unclear. Here, we derive cell-type-specific AD polygenic risk scores (ADPRS) from two extensively characterized datasets and uncover the impact of cell-type-specific genetic risk on AD endophenotypes. In an autopsy dataset spanning all stages of AD (n = 1457), the astrocytic ADPRS affected diffuse and neuritic plaques (amyloid-ß), while microglial ADPRS affected neuritic plaques, microglial activation, neurofibrillary tangles (tau), and cognitive decline. In an independent neuroimaging dataset of cognitively unimpaired elderly (n = 2921), astrocytic ADPRS was associated with amyloid-ß, and microglial ADPRS was associated with amyloid-ß and tau, connecting cell-type-specific genetic risk with AD pathology even before symptom onset. Together, our study provides human genetic evidence implicating multiple glial cell types in AD pathophysiology, starting from the preclinical stage.

Shared Genetic Risk in the Association of Screen Time With Psychiatric Problems in Children.

Importance: Children's exposure to screen time has been associated with poor mental health outcomes, yet the role of genetic factors remains largely unknown. Objective: To assess the extent of genetic confounding in the associations between screen time and attention problems or internalizing problems in preadolescent children. Design, Setting, and Participants: This cohort study analyzed data obtained between 2016 and 2019 from the Adolescent Brain Cognitive Development Study at 21 sites in the US. The sample included children aged 9 to 11 years of genetically assigned European ancestry with self-reported screen time. Data were analyzed between November 2021 and September 2023. Exposure: Child-reported daily screen time (in hours) was ascertained from questionnaires completed by the children at baseline. Main Outcomes and Measures: Child psychiatric problems, specifically attention and internalizing problems, were measured with the parent-completed Achenbach Child Behavior Checklist at the 1-year follow-up. Genetic sensitivity analyses model (Gsens) was used, which incorporated polygenic risk scores (PRSs) of both exposure and outcomes as well as either single-nucleotide variant (SNV; formerly single-nucleotide polymorphism)-based heritability or twin-based heritability to estimate genetic confounding. Results: The 4262 children in the sample included 2269 males (53.2%) with a mean (SD) age of 9.9 (0.6) years. Child screen time was associated with attention problems (ß = 0.10 SD; 95% CI, 0.07-0.13 SD) and internalizing problems (ß = 0.03 SD; 95% CI, 0.003-0.06 SD). The television time PRS was associated with child screen time (ß = 0.18 SD; 95% CI, 0.14-0.23 SD), the attention-deficit/hyperactivity disorder PRS was associated with attention problems (ß = 0.13 SD; 95% CI, 0.10-0.16 SD), and the depression PRS was associated with internalizing problems (ß = 0.10 SD; 95% CI, 0.07-0.13 SD). These PRSs were associated with cross-traits, suggesting genetic confounding. Estimates using PRSs and SNV-based heritability showed that genetic confounding accounted for most of the association between child screen time and attention problems and for 42.7% of the association between child screen time and internalizing problems. When PRSs and twin-based heritability estimates were used, genetic confounding fully explained both associations. Conclusions and Relevance: Results of this study suggest that genetic confounding may explain a substantial part of the associations between child screen time and psychiatric problems. Genetic confounding should be considered in sociobehavioral studies of modifiable factors for youth mental health.

Characterizing the phenotypic and genetic structure of psychopathology in UK Biobank.

Mental conditions exhibit a higher-order transdiagnostic factor structure which helps to explain the widespread comorbidity observed in psychopathology. However, the phenotypic and genetic structures of psychopathology may differ, raising questions about the validity and utility of these factors. Here, we study the phenotypic and genetic factor structures of ten psychiatric conditions using UK Biobank and public genomic data. Although the factor structure of psychopathology was generally genetically and phenotypically consistent, conditions related to externalizing (e.g., alcohol use disorder) and compulsivity (e.g., eating disorders) exhibited cross-level disparities in their relationships with other conditions, plausibly due to environmental influences. Domain-level factors, especially thought disorder and internalizing factors, were more informative than a general psychopathology factor in genome-wide association and polygenic index analyses. Collectively, our findings enhance the understanding of comorbidity and shared etiology, highlight the intricate interplay between genes and environment, and offer guidance for psychiatric research using polygenic indices.

Language network lateralization is reflected throughout the macroscale functional organization of cortex.

Hemispheric specialization is a fundamental feature of human brain organization. However, it is not yet clear to what extent the lateralization of specific cognitive processes may be evident throughout the broad functional architecture of cortex. While the majority of people exhibit left-hemispheric language dominance, a substantial minority of the population shows reverse lateralization. Using twin and family data from the Human Connectome Project, we provide evidence that atypical language dominance is associated with global shifts in cortical organization. Individuals with atypical language organization exhibit corresponding hemispheric differences in the macroscale functional gradients that situate discrete large-scale networks along a continuous spectrum, extending from unimodal through association territories. Analyses reveal that both language lateralization and gradient asymmetries are, in part, driven by genetic factors. These findings pave the way for a deeper understanding of the origins and relationships linking population-level variability in hemispheric specialization and global properties of cortical organization.

Transcriptome-Wide Structural Equation Modeling of 13 Major Psychiatric Disorders for Cross-Disorder Risk and Drug Repurposing.

Importance: Psychiatric disorders display high levels of comorbidity and genetic overlap, necessitating multivariate approaches for parsing convergent and divergent psychiatric risk pathways. Identifying gene expression patterns underlying cross-disorder risk also stands to propel drug discovery and repurposing in the face of rising levels of polypharmacy. Objective: To identify gene expression patterns underlying genetic convergence and divergence across psychiatric disorders along with existing pharmacological interventions that target these genes. Design, Setting, and Participants: This genomic study applied a multivariate transcriptomic method, transcriptome-wide structural equation modeling (T-SEM), to investigate gene expression patterns associated with 5 genomic factors indexing shared risk across 13 major psychiatric disorders. Follow-up tests, including overlap with gene sets for other outcomes and phenome-wide association studies, were conducted to better characterize T-SEM results. The Broad Institute Connectivity Map Drug Repurposing Database and Drug-Gene Interaction Database public databases of drug-gene pairs were used to identify drugs that could be repurposed to target genes found to be associated with cross-disorder risk. Data were collected from database inception up to February 20, 2023. Main Outcomes and Measures: Gene expression patterns associated with genomic factors or disorder-specific risk and existing drugs that target these genes. Results: In total, T-SEM identified 466 genes whose expression was significantly associated (z ≥ 5.02) with genomic factors and 36 genes with disorder-specific effects. Most associated genes were found for a thought disorders factor, defined by bipolar disorder and schizophrenia. Several existing pharmacological interventions were identified that could be repurposed to target genes whose expression was associated with the thought disorders factor or a transdiagnostic p factor defined by all 13 disorders. Conclusions and Relevance: The findings from this study shed light on patterns of gene expression associated with genetic overlap and uniqueness across psychiatric disorders. Future versions of the multivariate drug repurposing framework outlined here have the potential to identify novel pharmacological interventions for increasingly common, comorbid psychiatric presentations.