Genome-wide Association Studies in Ancestrally Diverse Populations: Opportunities, Methods, Pitfalls, and Recommendations.

Genome-wide association studies (GWASs) have focused primarily on populations of European descent, but it is essential that diverse populations become better represented. Increasing diversity among study participants will advance our understanding of genetic architecture in all populations and ensure that genetic research is broadly applicable. To facilitate and promote research in multi-ancestry and admixed cohorts, we outline key methodological considerations and highlight opportunities, challenges, solutions, and areas in need of development. Despite the perception that analyzing genetic data from diverse populations is difficult, it is scientifically and ethically imperative, and there is an expanding analytical toolbox to do it well.

Upregulated GIRK2 Counteracts Ethanol-Induced Changes in Excitability and Respiration in Human Neurons.

Genome-wide association studies (GWAS) of electroencephalographic endophenotypes for alcohol use disorder (AUD) has identified noncoding polymorphisms within the KCNJ6 gene. KCNJ6 encodes GIRK2, a subunit of a G-protein-coupled inwardly rectifying potassium channel that regulates neuronal excitability. We studied the effect of upregulating KCNJ6 using an isogenic approach with human glutamatergic neurons derived from induced pluripotent stem cells (male and female donors). Using multielectrode arrays, population calcium imaging, single-cell patch-clamp electrophysiology, and mitochondrial stress tests, we find that elevated GIRK2 acts in concert with 7-21 d of ethanol exposure to inhibit neuronal activity, to counteract ethanol-induced increases in glutamate response, and to promote an increase intrinsic excitability. Furthermore, elevated GIRK2 prevented ethanol-induced changes in basal and activity-dependent mitochondrial respiration. These data support a role for GIRK2 in mitigating the effects of ethanol and a previously unknown connection to mitochondrial function in human glutamatergic neurons.

Examining associations between genetic and neural risk for externalizing behaviors in adolescence and early adulthood.

BACKGROUND: Researchers have identified genetic and neural risk factors for externalizing behaviors. However, it has not yet been determined if genetic liability is conferred in part through associations with more proximal neurophysiological risk markers. METHODS: Participants from the Collaborative Study on the Genetics of Alcoholism, a large, family-based study of alcohol use disorders were genotyped and polygenic scores for externalizing (EXT PGS) were calculated. Associations with target P3 amplitude from a visual oddball task (P3) and broad endorsement of externalizing behaviors (indexed via self-report of alcohol and cannabis use, and antisocial behavior) were assessed in participants of European (EA; N = 2851) and African ancestry (AA; N = 1402). Analyses were also stratified by age (adolescents, age 12-17 and young adults, age 18-32). RESULTS: The EXT PGS was significantly associated with higher levels of externalizing behaviors among EA adolescents and young adults as well as AA young adults. P3 was inversely associated with externalizing behaviors among EA young adults. EXT PGS was not significantly associated with P3 amplitude and therefore, there was no evidence that P3 amplitude indirectly accounted for the association between EXT PGS and externalizing behaviors. CONCLUSIONS: Both the EXT PGS and P3 amplitude were significantly associated with externalizing behaviors among EA young adults. However, these associations with externalizing behaviors appear to be independent of each other, suggesting that they may index different facets of externalizing.

Alcohol milestones and internalizing, externalizing, and executive function: longitudinal and polygenic score associations.

BACKGROUND: Although the link between alcohol involvement and behavioral phenotypes (e.g. impulsivity, negative affect, executive function [EF]) is well-established, the directionality of these associations, specificity to stages of alcohol involvement, and extent of shared genetic liability remain unclear. We estimate longitudinal associations between transitions among alcohol milestones, behavioral phenotypes, and indices of genetic risk. METHODS: Data came from the Collaborative Study on the Genetics of Alcoholism (n = 3681; ages 11-36). Alcohol transitions (first: drink, intoxication, alcohol use disorder [AUD] symptom, AUD diagnosis), internalizing, and externalizing phenotypes came from the Semi-Structured Assessment for the Genetics of Alcoholism. EF was measured with the Tower of London and Visual Span Tasks. Polygenic scores (PGS) were computed for alcohol-related and behavioral phenotypes. Cox models estimated associations among PGS, behavior, and alcohol milestones. RESULTS: Externalizing phenotypes (e.g. conduct disorder symptoms) were associated with future initiation and drinking problems (hazard ratio (HR)⩾1.16). Internalizing (e.g. social anxiety) was associated with hazards for progression from first drink to severe AUD (HR⩾1.55). Initiation and AUD were associated with increased hazards for later depressive symptoms and suicidal ideation (HR⩾1.38), and initiation was associated with increased hazards for future conduct symptoms (HR = 1.60). EF was not associated with alcohol transitions. Drinks per week PGS was linked with increased hazards for alcohol transitions (HR⩾1.06). Problematic alcohol use PGS increased hazards for suicidal ideation (HR = 1.20). CONCLUSIONS: Behavioral markers of addiction vulnerability precede and follow alcohol transitions, highlighting dynamic, bidirectional relationships between behavior and emerging addiction.

A Developmentally-Informative Genome-wide Association Study of Alcohol Use Frequency.

Contemporary genome-wide association study (GWAS) methods typically do not account for variability in genetic effects throughout development. We applied genomic structural equation modeling to combine developmentally-informative phenotype data and GWAS to create polygenic scores (PGS) for alcohol use frequency that are specific to developmental stage. Longitudinal cohort studies targeted for gene-identification analyses include the Collaborative Study on the Genetics of Alcoholism (adolescence n = 1,118, early adulthood n = 2,762, adulthood n = 5,255), the National Longitudinal Study of Adolescent to Adult Health (adolescence n = 3,089, early adulthood n = 3,993, adulthood n = 5,149), and the Avon Longitudinal Study of Parents and Children (ALSPAC; adolescence n = 5,382, early adulthood n = 3,613). PGS validation analyses were conducted in the COGA sample using an alternate version of the discovery analysis with COGA removed. Results suggest that genetic liability for alcohol use frequency in adolescence may be distinct from genetic liability for alcohol use frequency later in developmental periods. The age-specific PGS predicts an increase of 4 drinking days per year per PGS standard deviation when modeled separately from the common factor PGS in adulthood. The current work was underpowered at all steps of the analysis plan. Though small sample sizes and low statistical power limit the substantive conclusions that can be drawn regarding these research questions, this work provides a foundation for future genetic studies of developmental variability in the genetic underpinnings of alcohol use behaviors and genetically-informed, age-matched phenotype prediction.

RNA alternative splicing impacts the risk for alcohol use disorder.

Alcohol use disorder (AUD) is a complex genetic disorder characterized by problems arising from excessive alcohol consumption. Identifying functional genetic variations that contribute to risk for AUD is a major goal. Alternative splicing of RNA mediates the flow of genetic information from DNA to gene expression and expands proteome diversity. We asked whether alternative splicing could be a risk factor for AUD. Herein, we used a Mendelian randomization (MR)-based approach to identify skipped exons (the predominant splicing event in brain) that contribute to AUD risk. Genotypes and RNA-seq data from the CommonMind Consortium were used as the training dataset to develop predictive models linking individual genotypes to exon skipping in the prefrontal cortex. We applied these models to data from the Collaborative Studies on Genetics of Alcoholism to examine the association between the imputed cis-regulated splicing outcome and the AUD-related traits. We identified 27 exon skipping events that were predicted to affect AUD risk; six of these were replicated in the Australian Twin-family Study of Alcohol Use Disorder. Their host genes are DRC1, ELOVL7, LINC00665, NSUN4, SRRM2 and TBC1D5. The genes downstream of these splicing events are enriched in neuroimmune pathways. The MR-inferred impacts of the ELOVL7 skipped exon on AUD risk was further supported in four additional large-scale genome-wide association studies. Additionally, this exon contributed to changes of gray matter volumes in multiple brain regions, including the visual cortex known to be involved in AUD. In conclusion, this study provides strong evidence that RNA alternative splicing impacts the susceptibility to AUD and adds new information on AUD-relevant genes and pathways. Our framework is also applicable to other types of splicing events and to other complex genetic disorders.

Genomic risk for post-traumatic stress disorder in families densely affected with alcohol use disorders.

Recent genome-wide association studies (GWAS) have identified genetic markers of post-traumatic stress disorder (PTSD) in civilian and military populations. However, studies have yet to examine the genetics of PTSD while factoring in risk for alcohol dependence, which commonly co-occur. We examined genome-wide associations for DSM-IV PTSD among 4,978 trauma-exposed participants (31% with alcohol dependence, 50% female, 30% African ancestry) from the Collaborative Study on the Genetics of Alcoholism (COGA). We also examined associations of polygenic risk scores (PRS) derived from the Psychiatric Genomics Consortium (PGC)-PTSD Freeze 2 (N = 3533) and Million Veterans Program GWAS of PTSD (N = 5200) with PTSD and substance dependence in COGA, and moderating effects of sex and alcohol dependence. 7.3% of COGA participants met criteria for PTSD, with higher rates in females (10.1%) and those with alcohol dependence (12.3%). No independent loci met genome-wide significance in the PTSD meta-analysis of European (EA) and African ancestry (AA) participants. The PGC-PTSD PRS was associated with increased risk for PTSD (B = 0.126, p < 0.001), alcohol dependence (B = 0.231, p < 0.001), and cocaine dependence (B = 0.086, p < 0.01) in EA individuals. A significant interaction was observed, such that EA individuals with alcohol dependence and higher polygenic risk for PTSD were more likely to have PTSD (B = 0.090, p < 0.01) than those without alcohol dependence. These results further support the importance of examining substance dependence, specifically alcohol dependence, and PTSD together when investigating genetic influence on these disorders.

Genetic nurture effects for alcohol use disorder.

We tested whether aspects of the childhood/adolescent home environment mediate genetic risk for alcohol problems within families across generations. Parental relationship discord and parental divorce were the focal environments examined. The sample included participants of European ancestry (N = 4806, 51% female) and African ancestry (N = 1960, 52% female) from the high-risk Collaborative Study on the Genetics of Alcoholism. Alcohol outcomes in the child generation included lifetime criterion counts for DSM-5 Alcohol Use Disorder (AUD), lifetime maximum drinks in 24 h, age at initiation of regular drinking, and age at first alcohol intoxication. Predictors in the parent generation included relationship discord, divorce, alcohol measures parallel to those in the child generation, and polygenic scores for alcohol problems. Parental polygenic scores were partitioned into alleles that were transmitted and non-transmitted to the child. The results from structural equation models were consistent with genetic nurture effects in European ancestry families. Exposure to parental relationship discord and parental divorce mediated, in part, the transmission of genetic risk for alcohol problems from parents to children to predict earlier ages regular drinking (ßindirect = -0.018 [-0.026, -0.011]) and intoxication (ßindirect = -0.015 [-0.023, -0.008]), greater lifetime maximum drinks (ßindirect = 0.006 [0.002, 0.01]) and more lifetime AUD criteria (ßindirect = 0.011 [0.006, 0.016]). In contrast, there was no evidence that parental alleles had indirect effects on offspring alcohol outcomes via parental relationship discord or divorce in the smaller number of families of African ancestry. In conclusion, parents transmit genetic risk for alcohol problems to their children not only directly, but also indirectly via genetically influenced aspects of the home environment. Further investigation of genetic nurture in non-European samples is needed.

Alcohol reverses the effects of KCNJ6 (GIRK2) noncoding variants on excitability of human glutamatergic neurons.

Synonymous and noncoding single nucleotide polymorphisms (SNPs) in the KCNJ6 gene, encoding G protein-gated inwardly rectifying potassium channel subunit 2 (GIRK2), have been linked with increased electroencephalographic frontal theta event-related oscillations (ERO) in subjects diagnosed with alcohol use disorder (AUD). To identify molecular and cellular mechanisms while retaining the appropriate genetic background, we generated induced excitatory glutamatergic neurons (iN) from iPSCs derived from four AUD-diagnosed subjects with KCNJ6 variants ("Affected: AF") and four control subjects without variants ("Unaffected: UN"). Neurons were analyzed for changes in gene expression, morphology, excitability and physiological properties. Single-cell RNA sequencing suggests that KCNJ6 AF variant neurons have altered patterns of synaptic transmission and cell projection morphogenesis. Results confirm that AF neurons express lower levels of GIRK2, have greater neurite area, and elevated excitability. Interestingly, exposure to intoxicating concentrations of ethanol induces GIRK2 expression and reverses functional effects in AF neurons. Ectopic overexpression of GIRK2 alone mimics the effect of ethanol to normalize induced excitability. We conclude that KCNJ6 variants decrease GIRK2 expression and increase excitability and that this effect can be minimized or reduced with ethanol.

A polygenic resilience score moderates the genetic risk for schizophrenia: Replication in 18,090 cases and 28,114 controls from the Psychiatric Genomics Consortium.

Identifying heritable factors that moderate the genetic risk for schizophrenia (SCZ) could help clarify why some individuals remain unaffected despite having relatively high genetic liability. Previously, we developed a framework to mine genome-wide association (GWAS) data for common genetic variants that protect high-risk unaffected individuals from SCZ, leading to derivation of the first-ever "polygenic resilience score" for SCZ (resilient controls n = 3786; polygenic risk score-matched SCZ cases n = 18,619). Here, we performed a replication study to verify the moderating effect of our polygenic resilience score on SCZ risk (OR = 1.09, p = 4.03 × 10-5 ) using newly released GWAS data from 23 independent case-control studies collated by the Psychiatric Genomics Consortium (PGC) (resilient controls n = 2821; polygenic risk score-matched SCZ cases n = 5150). Additionally, we sought to optimize our polygenic resilience-scoring formula to improve subsequent modeling of resilience to SCZ and other complex disorders. We found significant replication of the polygenic resilience score, and found that strict pruning of SNPs based on linkage disequilibrium to known risk SNPs and their linked loci optimizes the performance of the polygenic resilience score.

Differences in genetic correlations between posttraumatic stress disorder and alcohol-related problems phenotypes compared to alcohol consumption-related phenotypes.

BACKGROUND: Posttraumatic Stress Disorder (PTSD) tends to co-occur with greater alcohol consumption as well as alcohol use disorder (AUD). However, it is unknown whether the same etiologic factors that underlie PTSD-alcohol-related problems comorbidity also contribute to PTSD- alcohol consumption. METHODS: We used summary statistics from large-scale genome-wide association studies (GWAS) of European-ancestry (EA) and African-ancestry (AA) participants to estimate genetic correlations between PTSD and a range of alcohol consumption-related and alcohol-related problems phenotypes. RESULTS: In EAs, there were positive genetic correlations between PTSD phenotypes and alcohol-related problems phenotypes (e.g. Alcohol Use Disorders Identification Test (AUDIT) problem score) (rGs: 0.132-0.533, all FDR adjusted p < 0.05). However, the genetic correlations between PTSD phenotypes and alcohol consumption -related phenotypes (e.g. drinks per week) were negatively associated or non-significant (rGs: -0.417 to -0.042, FDR adjusted p: <0.05-NS). For AAs, the direction of correlations was sometimes consistent and sometimes inconsistent with that in EAs, and the ranges were larger (rGs for alcohol-related problems: -0.275 to 0.266, FDR adjusted p: NS, alcohol consumption-related: 0.145-0.699, FDR adjusted p: NS). CONCLUSIONS: These findings illustrate that the genetic associations between consumption and problem alcohol phenotypes and PTSD differ in both strength and direction. Thus, the genetic factors that may lead someone to develop PTSD and high levels of alcohol consumption are not the same as those that lead someone to develop PTSD and alcohol-related problems. Discussion around needing improved methods to better estimate heritabilities and genetic correlations in diverse and admixed ancestry samples is provided.

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

BACKGROUND: Alcohol use disorder (AUD) and schizophrenia (SCZ) frequently co-occur, and large-scale genome-wide association studies (GWAS) have identified significant genetic correlations between these disorders. METHODS: We used the largest published GWAS for AUD (total cases = 77 822) and SCZ (total cases = 46 827) to identify genetic variants that influence both disorders (with either the same or opposite direction of effect) and those that are disorder specific. RESULTS: We identified 55 independent genome-wide significant single nucleotide polymorphisms with the same direction of effect on AUD and SCZ, 8 with robust effects in opposite directions, and 98 with disorder-specific effects. We also found evidence for 12 genes whose pleiotropic associations with AUD and SCZ are consistent with mediation via gene expression in the prefrontal cortex. The genetic covariance between AUD and SCZ was concentrated in genomic regions functional in brain tissues (p = 0.001). CONCLUSIONS: Our findings provide further evidence that SCZ shares meaningful genetic overlap with AUD.

Genome-wide association study in individuals of European and African ancestry and multi-trait analysis of opioid use disorder identifies 19 independent genome-wide significant risk loci.

Despite the large toll of opioid use disorder (OUD), genome-wide association studies (GWAS) of OUD to date have yielded few susceptibility loci. We performed a large-scale GWAS of OUD in individuals of European (EUR) and African (AFR) ancestry, optimizing genetic informativeness by performing MTAG (Multi-trait analysis of GWAS) with genetically correlated substance use disorders (SUDs). Meta-analysis included seven cohorts: the Million Veteran Program, Psychiatric Genomics Consortium, iPSYCH, FinnGen, Partners Biobank, BioVU, and Yale-Penn 3, resulting in a total N = 639,063 (Ncases = 20,686;Neffective = 77,026) across ancestries. OUD cases were defined as having a lifetime OUD diagnosis, and controls as anyone not known to meet OUD criteria. We estimated SNP-heritability (h2SNP) and genetic correlations (rg). Based on genetic correlation, we performed MTAG on OUD, alcohol use disorder (AUD), and cannabis use disorder (CanUD). A leave-one-out polygenic risk score (PRS) analysis was performed to compare OUD and OUD-MTAG PRS as predictors of OUD case status in Yale-Penn 3. The EUR meta-analysis identified three genome-wide significant (GWS; p ≤ 5 × 10-8) lead SNPs-one at FURIN (rs11372849; p = 9.54 × 10-10) and two OPRM1 variants (rs1799971, p = 4.92 × 10-09; rs79704991, p = 1.11 × 10-08; r2 = 0.02). Rs1799971 (p = 4.91 × 10-08) and another OPRM1 variant (rs9478500; p = 1.95 × 10-08; r2 = 0.03) were identified in the cross-ancestry meta-analysis. Estimated h2SNP was 12.75%, with strong rg with CanUD (rg = 0.82; p = 1.14 × 10-47) and AUD (rg = 0.77; p = 6.36 × 10-78). The OUD-MTAG resulted in a GWAS Nequivalent = 128,748 and 18 independent GWS loci, some mapping to genes or gene regions that have previously been associated with psychiatric or addiction phenotypes. The OUD-MTAG PRS accounted for 3.81% of OUD variance (beta = 0.61;s.e. = 0.066; p = 2.00 × 10-16) compared to 2.41% (beta = 0.45; s.e. = 0.058; p = 2.90 × 10-13) explained by the OUD PRS. The current study identified OUD variant associations at OPRM1, single variant associations with FURIN, and 18 GWS associations in the OUD-MTAG. The genetic architecture of OUD is likely influenced by both OUD-specific loci and loci shared across SUDs.

Rare variants implicate NMDA receptor signaling and cerebellar gene networks in risk for bipolar disorder.

Bipolar disorder is an often-severe mental health condition characterized by alternation between extreme mood states of mania and depression. Despite strong heritability and the recent identification of 64 common variant risk loci of small effect, pathophysiological mechanisms remain unknown. Here, we analyzed genome sequences from 41 multiply-affected pedigrees and identified variants in 741 genes with nominally significant linkage or association with bipolar disorder. These 741 genes overlapped known risk genes for neurodevelopmental disorders and clustered within gene networks enriched for synaptic and nuclear functions. The top variant in this analysis - prioritized by statistical association, predicted deleteriousness, and network centrality - was a missense variant in the gene encoding D-amino acid oxidase (DAOG131V). Heterologous expression of DAOG131V in human cells resulted in decreased DAO protein abundance and enzymatic activity. In a knock-in mouse model of DAOG131, DaoG130V/+, we similarly found decreased DAO protein abundance in hindbrain regions, as well as enhanced stress susceptibility and blunted behavioral responses to pharmacological inhibition of N-methyl-D-aspartate receptors (NMDARs). RNA sequencing of cerebellar tissue revealed that DaoG130V resulted in decreased expression of two gene networks that are enriched for synaptic functions and for genes expressed, respectively, in Purkinje neurons or granule neurons. These gene networks were also down-regulated in the cerebellum of patients with bipolar disorder compared to healthy controls and were enriched for additional rare variants associated with bipolar disorder risk. These findings implicate dysregulation of NMDAR signaling and of gene expression in cerebellar neurons in bipolar disorder pathophysiology and provide insight into its genetic architecture.

Clinical, environmental, and genetic risk factors for substance use disorders: characterizing combined effects across multiple cohorts.

Substance use disorders (SUDs) incur serious social and personal costs. The risk for SUDs is complex, with risk factors ranging from social conditions to individual genetic variation. We examined whether models that include a clinical/environmental risk index (CERI) and polygenic scores (PGS) are able to identify individuals at increased risk of SUD in young adulthood across four longitudinal cohorts for a combined sample of N = 15,134. Our analyses included participants of European (NEUR = 12,659) and African (NAFR = 2475) ancestries. SUD outcomes included: (1) alcohol dependence, (2) nicotine dependence; (3) drug dependence, and (4) any substance dependence. In the models containing the PGS and CERI, the CERI was associated with all three outcomes (ORs = 01.37-1.67). PGS for problematic alcohol use, externalizing, and smoking quantity were associated with alcohol dependence, drug dependence, and nicotine dependence, respectively (OR = 1.11-1.33). PGS for problematic alcohol use and externalizing were also associated with any substance dependence (ORs = 1.09-1.18). The full model explained 6-13% of the variance in SUDs. Those in the top 10% of CERI and PGS had relative risk ratios of 3.86-8.04 for each SUD relative to the bottom 90%. Overall, the combined measures of clinical, environmental, and genetic risk demonstrated modest ability to distinguish between affected and unaffected individuals in young adulthood. PGS were significant but added little in addition to the clinical/environmental risk index. Results from our analysis demonstrate there is still considerable work to be done before tools such as these are ready for clinical applications.

Alcohol use polygenic risk score, social support, and alcohol use among European American and African American adults.

Alcohol use is influenced by genetic and environmental factors. We examined the interactive effects between genome-wide polygenic risk scores for alcohol use (alc-PRS) and social support in relation to alcohol use among European American (EA) and African American (AA) adults across sex and developmental stages (emerging adulthood, young adulthood, and middle adulthood). Data were drawn from 4,011 EA and 1,274 AA adults from the Collaborative Study on the Genetics of Alcoholism who were between ages 18-65 and had ever used alcohol. Participants completed the Semi-Structured Assessment for the Genetics of Alcoholism and provided saliva or blood samples for genotyping. Results indicated that social support from friends, but not family, moderated the association between alc-PRS and alcohol use among EAs and AAs (only in middle adulthood for AAs); alc-PRS was associated with higher levels of alcohol use when friend support was low, but not when friend support was high. Associations were similar across sex but differed across developmental stages. Findings support the important role of social support from friends in buffering genetic risk for alcohol use among EA and AA adults and highlight the need to consider developmental changes in the role of social support in relation to alcohol use.

Principal Component Analysis Reduces Collider Bias in Polygenic Score Effect Size Estimation.

In this study, we test principal component analysis (PCA) of measured confounders as a method to reduce collider bias in polygenic association models. We present results from simulations and application of the method in the Collaborative Study of the Genetics of Alcoholism (COGA) sample with a polygenic score for alcohol problems, DSM-5 alcohol use disorder as the target phenotype, and two collider variables: tobacco use and educational attainment. Simulation results suggest that assumptions regarding the correlation structure and availability of measured confounders are complementary, such that meeting one assumption relaxes the other. Application of the method in COGA shows that PC covariates reduce collider bias when tobacco use is used as the collider variable. Application of this method may improve PRS effect size estimation in some cases by reducing the effect of collider bias, making efficient use of data resources that are available in many studies.

Predicting risk for Alcohol Use Disorder using longitudinal data with multimodal biomarkers and family history: a machine learning study.

Predictive models have succeeded in distinguishing between individuals with Alcohol use Disorder (AUD) and controls. However, predictive models identifying who is prone to develop AUD and the biomarkers indicating a predisposition to AUD are still unclear. Our sample (n = 656) included offspring and non-offspring of European American (EA) and African American (AA) ancestry from the Collaborative Study of the Genetics of Alcoholism (COGA) who were recruited as early as age 12 and were unaffected at first assessment and reassessed years later as AUD (DSM-5) (n = 328) or unaffected (n = 328). Machine learning analysis was performed for 220 EEG measures, 149 alcohol-related single nucleotide polymorphisms (SNPs) from a recent large Genome-wide Association Study (GWAS) of alcohol use/misuse and two family history (mother DSM-5 AUD and father DSM-5 AUD) features using supervised, Linear Support Vector Machine (SVM) classifier to test which features assessed before developing AUD predict those who go on to develop AUD. Age, gender, and ancestry stratified analyses were performed. Results indicate significant and higher accuracy rates for the AA compared with the EA prediction models and a higher model accuracy trend among females compared with males for both ancestries. Combined EEG and SNP features model outperformed models based on only EEG features or only SNP features for both EA and AA samples. This multidimensional superiority was confirmed in a follow-up analysis in the AA age groups (12-15, 16-19, 20-30) and EA age group (16-19). In both ancestry samples, the youngest age group achieved higher accuracy score than the two other older age groups. Maternal AUD increased the model's accuracy in both ancestries' samples. Several discriminative EEG measures and SNPs features were identified, including lower posterior gamma, higher slow wave connectivity (delta, theta, alpha), higher frontal gamma ratio, higher beta correlation in the parietal area, and 5 SNPs: rs4780836, rs2605140, rs11690265, rs692854, and rs13380649. Results highlight the significance of sampling uniformity followed by stratified (e.g., ancestry, gender, developmental period) analysis, and wider selection of features, to generate better prediction scores allowing a more accurate estimation of AUD development.

Allele-specific expression and high-throughput reporter assay reveal functional genetic variants associated with alcohol use disorders.

Genome-wide association studies (GWAS) of complex traits, such as alcohol use disorders (AUD), usually identify variants in non-coding regions and cannot by themselves distinguish whether the associated variants are functional or in linkage disequilibrium with the functional variants. Transcriptome studies can identify genes whose expression differs between alcoholics and controls. To test which variants associated with AUD may cause expression differences, we integrated data from deep RNA-seq and GWAS of four postmortem brain regions from 30 subjects with AUD and 30 controls to analyze allele-specific expression (ASE). We identified 88 genes with differential ASE in subjects with AUD compared to controls. Next, to test one potential mechanism contributing to the differential ASE, we analyzed single nucleotide polymorphisms (SNPs) in the 3' untranslated regions (3'UTR) of these genes. Of the 88 genes with differential ASE, 61 genes contained 437 SNPs in the 3'UTR with at least one heterozygote among the subjects studied. Using a modified PASSPORT-seq (parallel assessment of polymorphisms in miRNA target-sites by sequencing) assay, we identified 25 SNPs that affected RNA levels in a consistent manner in two neuroblastoma cell lines, SH-SY5Y and SK-N-BE(2). Many of these SNPs are in binding sites of miRNAs and RNA-binding proteins, indicating that these SNPs are likely causal variants of AUD-associated differential ASE. In sum, we demonstrate that a combination of computational and experimental approaches provides a powerful strategy to uncover functionally relevant variants associated with the risk for AUD.

A polygenic resilience score moderates the genetic risk for schizophrenia.

Based on the discovery by the Resilience Project (Chen R. et al. Nat Biotechnol 34:531-538, 2016) of rare variants that confer resistance to Mendelian disease, and protective alleles for some complex diseases, we posited the existence of genetic variants that promote resilience to highly heritable polygenic disorders1,0 such as schizophrenia. Resilience has been traditionally viewed as a psychological construct, although our use of the term resilience refers to a different construct that directly relates to the Resilience Project, namely: heritable variation that promotes resistance to disease by reducing the penetrance of risk loci, wherein resilience and risk loci operate orthogonal to one another. In this study, we established a procedure to identify unaffected individuals with relatively high polygenic risk for schizophrenia, and contrasted them with risk-matched schizophrenia cases to generate the first known "polygenic resilience score" that represents the additive contributions to SZ resistance by variants that are distinct from risk loci. The resilience score was derived from data compiled by the Psychiatric Genomics Consortium, and replicated in three independent samples. This work establishes a generalizable framework for finding resilience variants for any complex, heritable disorder.