Smoller, Jordan W.

ABSTRACT Amidst the opioid crisis, understanding the genetic basis of opioid use disorder (OUD) is crucial for identifying biological mechanisms and intervention points. However, genome‐wide association studies (GWASs) have been hampered by inadequate sample sizes and often the use of control populations not assessed for prior opioid exposure. Because opioid exposure is a prerequisite for the development of OUD, consideration of exposure history in controls is important. Electronic health record data (EHR) paired with genomic information allow a broader sampling of patients with OUD and exposed controls. We leveraged data across two healthcare systems to evaluate the impact of using controls not screened for opioid exposure (‘generic’) versus minimally opioid‐exposed control (‘exposed’). First, at the phenotypic level, we conducted phenome‐wide association studies (PheWAS) to compare the medical comorbidity profiles of OUD cases when using generic versus exposed controls. While PheWAS results for OUD‐related comorbidities were more pronounced when using the generic group, 83% of the disease associations were overlapping and of similar effect sizes. Second, at the genetic level, we conducted GWAS (cases vs. generic; cases vs. exposed) and assessed differences in genetic correlations and degrees of phenotypic misclassification. Genetic results were concordant across control groups based on heritability (generic: 0.16 ± 0.07 vs. 0.10 ± 0.07), associations with the coding OPRM1 variant rs1799971 ( p generic = 8.83E‐03 vs. p exposed = 1.83E‐02) and genetic correlations with prior OUD GWAS ( r g‐generic = 0.83 ± 0.26 vs. r g‐exposed = 0.78 ± 0.27). Although GWASs were limited by sample size ( N generic = 6269, N exposed = 6365), compared to an independent OUD GWAS ( N = 425 944), the dilution value for the two GWAS was not different from 1, suggesting no major impact of phenotypic misclassification. This study represents the first effort to enhance OUD genetic research through optimization of control definitions using EHR data. Generic controls ascertained within the US health systems, where exposure to prescription opioids is high, offer a practical alternative for genetic studies of OUD.

ABSTRACT Bipolar disorder’s (BD) clinical heterogeneity has an unresolved genetic basis. We meta-analyzed genome-wide association studies (GWAS) of 16 BD subphenotypes in 226,032 individuals from 57 cohorts (38,022 cases); 10 advanced to multivariate and multi-trait analyses. Four factors (compulsive, psychotic, dysregulated, internalizing) explained 82.8% of shared genetic variance. BD1 and BD2 loaded on distinct factors despite a high genetic correlation; 87.0% of common-factor loci were significant in neither subtype. Unipolar mania aligned with psychosis over internalizing, and was distinguishable from BD1, and rapid cycling showed heritable cross-domain liability. We identified 356 risk loci, 158 novel, including the first univariate-GWAS associations for psychosis, unipolar mania, rapid cycling and schizoaffective disorder—and 249 credible genes (89 high-confidence), 12 with approved-drug or clinical-phase annotations. Cell-type association showed a midbrain dopaminergic–GABAergic gradient along the psychotic factor. BD’s genetic architecture appears hierarchical—a general liability resolving into dimensions of course and comorbidity beyond diagnostic subtypes.

Abstract Current studies regarding the secondary use of electronic health records (EHR) predominantly rely on domain expertise and existing medical knowledge. Though significant efforts have been devoted to investigating the application of machine learning algorithms in the EHR, efficient and powerful representation of patients is needed to unleash the potential of discovering new medical patterns underlying the EHR. Here, we present an unsupervised method for embedding high-dimensional EHR data at the patient level, aimed at characterizing patient heterogeneity in complex diseases and identifying new disease patterns associated with clinical outcome disparities. Inspired by the architecture of modern language models—specifically transformers with attention mechanisms, we use patient diagnosis and procedure codes as vocabularies and treat each patient as a sentence to perform the patient embedding. We applied this approach to 34,851 unique medical codes across 1,046,649 longitudinal patient events, including 102,739 patients from the electronic Medical Records and GEnomics (eMERGE) Network. The resulting patient vectors demonstrated excellent performance in predicting future disease events (median AUROC = 0.87 within one year) and bulk phenotyping (median AUROC = 0.84). We then illustrated the utility of these patient vectors in revealing heterogeneous comorbidity patterns, exemplified by disease subtypes in colorectal cancer and systemic lupus erythematosus, and capturing distinct longitudinal disease trajectories. External validation using EHR data from the University of Washington confirmed robust model performance, with median AUROCs of 0.83 and 0.84 for bulk phenotyping tasks and disease onset prediction, respectively. Importantly, the model reproduced the clustering results of disease subtypes identified in the eMERGE cohort and uncovered variations in overall mortality among these subtypes. Together, these results underscore the potential of representation learning in EHRs to enhance patient characterization and associated clinical outcomes, thereby advancing disease forecasting and facilitating personalized medicine.

INTRODUCTION The cerebral cortex underlies our complex cognitive capabilities. Variations in human cortical surface area and thickness are associated with neurological, psychological, and behavioral traits and can be measured in vivo by magnetic resonance imaging (MRI). Studies in model organisms have identified genes that influence cortical structure, but little is known about common genetic variants that affect human cortical structure. RATIONALE To identify genetic variants associated with human cortical structure at both global and regional levels, we conducted a genome-wide association meta-analysis of brain MRI data from 51,665 individuals across 60 cohorts. We analyzed the surface area and average thickness of the whole cortex and 34 cortical regions with known functional specializations. RESULTS We identified 369 nominally genome-wide significant loci ( P < 5 × 10 −8 ) associated with cortical structure in a discovery sample of 33,992 participants of European ancestry. Of the 360 loci for which replication data were available, 241 loci influencing surface area and 66 influencing thickness remained significant after replication, with 237 loci passing multiple testing correction ( P < 8.3 × 10 −10 ; 187 influencing surface area and 50 influencing thickness). Common genetic variants explained 34% (SE = 3%) of the variation in total surface area and 26% (SE = 2%) in average thickness; surface area and thickness showed a negative genetic correlation ( r G = −0.32, SE = 0.05, P = 6.5 × 10 −12 ), which suggests that genetic influences have opposing effects on surface area and thickness. Bioinformatic analyses showed that total surface area is influenced by genetic variants that alter gene regulatory activity in neural progenitor cells during fetal development. By contrast, average thickness is influenced by active regulatory elements in adult brain samples, which may reflect processes that occur after mid-fetal development, such as myelination, branching, or pruning. When considered together, these results support the radial unit hypothesis that different developmental mechanisms promote surface area expansion and increases in thickness. To identify specific genetic influences on individual cortical regions, we controlled for global measures (total surface area or average thickness) in the regional analyses. After multiple testing correction, we identified 175 loci that influence regional surface area and 46 that influence regional thickness. Loci that affect regional surface area cluster near genes involved in the Wnt signaling pathway, which is known to influence areal identity. We observed significant positive genetic correlations and evidence of bidirectional causation of total surface area with both general cognitive functioning and educational attainment. We found additional positive genetic correlations between total surface area and Parkinson’s disease but did not find evidence of causation. Negative genetic correlations were evident between total surface area and insomnia, attention deficit hyperactivity disorder, depressive symptoms, major depressive disorder, and neuroticism. CONCLUSION This large-scale collaborative work enhances our understanding of the genetic architecture of the human cerebral cortex and its regional patterning. The highly polygenic architecture of the cortex suggests that distinct genes are involved in the development of specific cortical areas. Moreover, we find evidence that brain structure is a key phenotype along the causal pathway that leads from genetic variation to differences in general cognitive function. Identifying genetic influences on human cortical structure. ( A ) Measurement of cortical surface area and thickness from MRI. ( B ) Genomic locations of common genetic variants that influence global and regional cortical structure. ( C ) Our results support the radial unit hypothesis that the expansion of cortical surface area is driven by proliferating neural progenitor cells. ( D ) Cortical surface area shows genetic correlation with psychiatric and cognitive traits. Error bars indicate SE. IMAGE CREDITS: (A) K. COURTNEY; (C) M. R. GLASS

We examined the Institutional Review Board (IRB) process at 9 academic institutions in the electronic Medical Records and Genomics (eMERGE) Network, for proposed electronic health record-based genomic medicine studies, to identify common questions and concerns. Sequencing of 109 disease related genes and genotyping of 14 actionable variants is being performed in ~28,100 participants from the 9 sites. Pathogenic/likely pathogenic variants in actionable genes are being returned to study participants. We examined each site’s research protocols, informed-consent materials, and interactions with IRB staff. Research staff at each site completed questionnaires regarding their IRB interactions. The time to prepare protocols for IRB submission, number of revisions and time to approval ranged from 10–261 days, 0–11, and 11–90 days, respectively. IRB recommendations related to the readability of informed consent materials, specifying the full range of potential risks, providing options for receiving limited results or withdrawal, sharing of information with family members, and establishing the mechanisms to answer participant questions. IRBs reviewing studies that involve the return of results from genomic sequencing have a diverse array of concerns, and anticipating these concerns can help investigators to more effectively engage IRBs.