A phenome-wide association and Mendelian randomisation study of alcohol use variants in a diverse cohort comprising over 3 million individuals.

BACKGROUND: Alcohol consumption is associated with numerous negative social and health outcomes. These associations may be direct consequences of drinking, or they may reflect common genetic factors that influence both alcohol consumption and other outcomes. METHODS: We performed exploratory phenome-wide association studies (PheWAS) of three of the best studied protective single nucleotide polymorphisms (SNPs) in genes encoding ethanol metabolising enzymes (ADH1B: rs1229984-T, rs2066702-A; ADH1C: rs698-T) using up to 1109 health outcomes across 28 phenotypic categories (e.g., substance-use, mental health, sleep, immune, cardiovascular, metabolic) from a diverse 23andMe cohort, including European (N ≤ 2,619,939), Latin American (N ≤ 446,646) and African American (N ≤ 146,776) populations to uncover new and perhaps unexpected associations. These SNPs have been consistently implicated by both candidate gene studies and genome-wide association studies of alcohol-related behaviours but have not been investigated in detail for other relevant phenotypes in a hypothesis-free approach in such a large cohort of multiple ancestries. To provide insight into potential causal effects of alcohol consumption on the outcomes significant in the PheWAS, we performed univariable two-sample and one-sample Mendelian randomisation (MR) analyses. FINDINGS: The minor allele rs1229984-T, which is protective against alcohol behaviours, showed the highest number of PheWAS associations across the three cohorts (N = 232, European; N = 29, Latin American; N = 7, African American). rs1229984-T influenced multiple domains of health. We replicated associations with alcohol-related behaviours, mental and sleep conditions, and cardio-metabolic health. We also found associations with understudied traits related to neurological (migraines, epilepsy), immune (allergies), musculoskeletal (fibromyalgia), and reproductive health (preeclampsia). MR analyses identified evidence of causal effects of alcohol consumption on liability for 35 of these outcomes in the European cohort. INTERPRETATION: Our work demonstrates that polymorphisms in genes encoding alcohol metabolising enzymes affect multiple domains of health beyond alcohol-related behaviours. Understanding the underlying mechanisms of these effects could have implications for treatments and preventative medicine. FUNDING: MVJ, NCK, SBB, SSR and AAP were supported by T32IR5226 and 28IR-0070. SSR was also supported by NIDA DP1DA054394. NCK and RBC were also supported by R25MH081482. ASH was supported by funds from NIAAA K01AA030083. JLMO was supported by VA 1IK2CX002095. JLMO and JJMM were also supported by NIDA R21DA050160. JJMM was also supported by the Kavli Postdoctoral Award for Academic Diversity. EGA was supported by K01MH121659 from the NIMH/NIH, the Caroline Wiess Law Fund for Research in Molecular Medicine and the ARCO Foundation Young Teacher-Investigator Fund at Baylor College of Medicine. MSA was supported by the Instituto de Salud Carlos III and co-funded by the European Union Found: Fondo Social Europeo Plus (FSE+) (P19/01224, PI22/00464 and CP22/00128).

Latin American Trans-ancestry INitiative for OCD genomics (LATINO): Study protocol.

Obsessive-compulsive disorder (OCD) is a debilitating psychiatric disorder. Worldwide, its prevalence is ~2% and its etiology is mostly unknown. Identifying biological factors contributing to OCD will elucidate underlying mechanisms and might contribute to improved treatment outcomes. Genomic studies of OCD are beginning to reveal long-sought risk loci, but >95% of the cases currently in analysis are of homogenous European ancestry. If not addressed, this Eurocentric bias will result in OCD genomic findings being more accurate for individuals of European ancestry than other ancestries, thereby contributing to health disparities in potential future applications of genomics. In this study protocol paper, we describe the Latin American Trans-ancestry INitiative for OCD genomics (LATINO, https://www.latinostudy.org). LATINO is a new network of investigators from across Latin America, the United States, and Canada who have begun to collect DNA and clinical data from 5000 richly phenotyped OCD cases of Latin American ancestry in a culturally sensitive and ethical manner. In this project, we will utilize trans-ancestry genomic analyses to accelerate the identification of OCD risk loci, fine-map putative causal variants, and improve the performance of polygenic risk scores in diverse populations. We will also capitalize on rich clinical data to examine the genetics of treatment response, biologically plausible OCD subtypes, and symptom dimensions. Additionally, LATINO will help elucidate the diversity of the clinical presentations of OCD across cultures through various trainings developed and offered in collaboration with Latin American investigators. We believe this study will advance the important goal of global mental health discovery and equity.

Cross-Species Convergence of Brain Transcriptomic and Epigenomic Findings in Posttraumatic Stress Disorder: A Systematic Review.

Introduction: Posttraumatic stress disorder (PTSD) is a complex multifactorial disorder influenced by the interaction of genetic and environmental factors. Analyses of epigenomic and transcriptomic modifications may help to dissect the biological factors underlying the gene-environment interplay in PTSD. To date, most human PTSD epigenetics studies have used peripheral tissue, and these findings have complex and poorly understood relationships to brain alterations. Studies examining brain tissue may help characterize the brain-specific transcriptomic and epigenomic profiles of PTSD. In this review, we compiled and integrated brain-specific molecular findings of PTSD from humans and animals. Methods: A systematic literature search according to the PRISMA criteria was performed to identify transcriptomic and epigenomic studies of PTSD, focusing on brain tissue from human postmortem samples or animal-stress paradigms. Results: Gene- and pathway-level convergence analyses revealed PTSD-dysregulated genes and biological pathways across brain regions and species. A total of 243 genes converged across species, with 17 of them significantly enriched for PTSD. Chemical synaptic transmission and signaling by G-protein-coupled receptors were consistently enriched across omics and species. Discussion: Our findings point out dysregulated genes highly replicated across PTSD studies in humans and animal models and suggest a potential role for the corticotropin-releasing hormone/orexin pathway in PTSD's pathophysiology. Further, we highlight current knowledge gaps and limitations and recommend future directions to address them.

Latin American Trans-ancestry INitiative for OCD genomics (LATINO): Study Protocol.

Obsessive-compulsive disorder (OCD) is a debilitating psychiatric disorder. Worldwide, its prevalence is ~2% and its etiology is mostly unknown. Identifying biological factors contributing to OCD will elucidate underlying mechanisms and might contribute to improved treatment outcomes. Genomic studies of OCD are beginning to reveal long-sought risk loci, but >95% of the cases currently in analysis are of homogenous European ancestry. If not addressed, this Eurocentric bias will result in OCD genomic findings being more accurate for individuals of European ancestry than other ancestries, thereby contributing to health disparities in potential future applications of genomics. In this study protocol paper, we describe the Latin American Trans-ancestry INitiative for OCD genomics (LATINO, www.latinostudy.org). LATINO is a new network of investigators from across Latin America, the United States, and Canada who have begun to collect DNA and clinical data from 5,000 richly-phenotyped OCD cases of Latin American ancestry in a culturally sensitive and ethical manner. In this project, we will utilize trans-ancestry genomic analyses to accelerate the identification of OCD risk loci, fine-map putative causal variants, and improve the performance of polygenic risk scores in diverse populations. We will also capitalize on rich clinical data to examine the genetics of treatment response, biologically plausible OCD subtypes, and symptom dimensions. Additionally, LATINO will help elucidate the diversity of the clinical presentations of OCD across cultures through various trainings developed and offered in collaboration with Latin American investigators. We believe this study will advance the important goal of global mental health discovery and equity.

Rare tandem repeat expansions associate with genes involved in synaptic and neuronal signaling functions in schizophrenia.

Tandem repeat expansions (TREs) are associated with over 60 monogenic disorders and have recently been implicated in complex disorders such as cancer and autism spectrum disorder. The role of TREs in schizophrenia is now emerging. In this study, we have performed a genome-wide investigation of TREs in schizophrenia. Using genome sequence data from 1154 Swedish schizophrenia cases and 934 ancestry-matched population controls, we have detected genome-wide rare (<0.1% population frequency) TREs that have motifs with a length of 2-20 base pairs. We find that the proportion of individuals carrying rare TREs is significantly higher in the schizophrenia group. There is a significantly higher burden of rare TREs in schizophrenia cases than in controls in genic regions, particularly in postsynaptic genes, in genes overlapping brain expression quantitative trait loci, and in brain-expressed genes that are differentially expressed between schizophrenia cases and controls. We demonstrate that TRE-associated genes are more constrained and primarily impact synaptic and neuronal signaling functions. These results have been replicated in an independent Canadian sample that consisted of 252 schizophrenia cases of European ancestry and 222 ancestry-matched controls. Our results support the involvement of rare TREs in schizophrenia etiology.

Recommendations to encourage participation of individuals from diverse backgrounds in psychiatric genetic studies.

We present innovative research practices in psychiatric genetic studies to ensure representation of individuals from diverse ancestry, sex assigned at birth, gender identity, age, body shape and size, and socioeconomic backgrounds. Due to histories of inappropriate and harmful practices against marginalized groups in both psychiatry and genetics, people of certain identities may be hesitant to participate in research studies. Yet their participation is essential to ensure diverse representation, as it is incorrect to assume that the same genetic and environmental factors influence the risk for various psychiatric disorders across all demographic groups. We present approaches developed as part of the Eating Disorders Genetics Initiative (EDGI), a study that required tailored approaches to recruit diverse populations across many countries. Considerations include research priorities and design, recruitment and study branding, transparency, and community investment and ownership. Ensuring representation in participants is costly and funders need to provide adequate support to achieve diversity in recruitment in prime awards, not just as supplemental afterthoughts. The need for diverse samples in genetic studies is critical to minimize the risk of perpetuating health disparities in psychiatry and other health research. Although the EDGI strategies were designed specifically to attract and enroll individuals with eating disorders, our approach is broadly applicable across psychiatry and other fields.

Gene expression changes following chronic antipsychotic exposure in single cells from mouse striatum.

Schizophrenia is an idiopathic psychiatric disorder with a high degree of polygenicity. Evidence from genetics, single-cell transcriptomics, and pharmacological studies suggest an important, but untested, overlap between genes involved in the etiology of schizophrenia and the cellular mechanisms of action of antipsychotics. To directly compare genes with antipsychotic-induced differential expression to genes involved in schizophrenia, we applied single-cell RNA-sequencing to striatal samples from male C57BL/6 J mice chronically exposed to a typical antipsychotic (haloperidol), an atypical antipsychotic (olanzapine), or placebo. We identified differentially expressed genes in three cell populations identified from the single-cell RNA-sequencing (medium spiny neurons [MSNs], microglia, and astrocytes) and applied multiple analysis pipelines to contextualize these findings, including comparison to GWAS results for schizophrenia. In MSNs in particular, differential expression analysis showed that there was a larger share of differentially expressed genes (DEGs) from mice treated with olanzapine compared with haloperidol. DEGs were enriched in loci implicated by genetic studies of schizophrenia, and we highlighted nine genes with convergent evidence. Pathway analyses of gene expression in MSNs highlighted neuron/synapse development, alternative splicing, and mitochondrial function as particularly engaged by antipsychotics. In microglia, we identified pathways involved in microglial activation and inflammation as part of the antipsychotic response. In conclusion, single-cell RNA sequencing may provide important insights into antipsychotic mechanisms of action and links to findings from psychiatric genomic studies.

THUNDER: A reference-free deconvolution method to infer cell type proportions from bulk Hi-C data.

Hi-C data provide population averaged estimates of three-dimensional chromatin contacts across cell types and states in bulk samples. Effective analysis of Hi-C data entails controlling for the potential confounding factor of differential cell type proportions across heterogeneous bulk samples. We propose a novel unsupervised deconvolution method for inferring cell type composition from bulk Hi-C data, the Two-step Hi-c UNsupervised DEconvolution appRoach (THUNDER). We conducted extensive simulations to test THUNDER based on combining two published single-cell Hi-C (scHi-C) datasets. THUNDER more accurately estimates the underlying cell type proportions compared to reference-free methods (e.g., TOAST, and NMF) and is more robust than reference-dependent methods (e.g. MuSiC). We further demonstrate the practical utility of THUNDER to estimate cell type proportions and identify cell-type-specific interactions in Hi-C data from adult human cortex tissue samples. THUNDER will be a useful tool in adjusting for varying cell type composition in population samples, facilitating valid and more powerful downstream analysis such as differential chromatin organization studies. Additionally, THUNDER estimated contact profiles provide a useful exploratory framework to investigate cell-type-specificity of the chromatin interactome while experimental data is still rare.

Bayesian modeling of skewed X inactivation in genetically diverse mice identifies a novel Xce allele associated with copy number changes.

Female mammals are functional mosaics of their parental X-linked gene expression due to X chromosome inactivation (XCI). This process inactivates one copy of the X chromosome in each cell during embryogenesis and that state is maintained clonally through mitosis. In mice, the choice of which parental X chromosome remains active is determined by the X chromosome controlling element (Xce), which has been mapped to a 176-kb candidate interval. A series of functional Xce alleles has been characterized or inferred for classical inbred strains based on biased, or skewed, inactivation of the parental X chromosomes in crosses between strains. To further explore the function structure basis and location of the Xce, we measured allele-specific expression of X-linked genes in a large population of F1 females generated from Collaborative Cross (CC) strains. Using published sequence data and applying a Bayesian "Pólya urn" model of XCI skew, we report two major findings. First, inter-individual variability in XCI suggests mouse epiblasts contain on average 20-30 cells contributing to brain. Second, CC founder strain NOD/ShiLtJ has a novel and unique functional allele, Xceg, that is the weakest in the Xce allelic series. Despite phylogenetic analysis confirming that NOD/ShiLtJ carries a haplotype almost identical to the well-characterized C57BL/6J (Xceb), we observed unexpected patterns of XCI skewing in females carrying the NOD/ShiLtJ haplotype within the Xce. Copy number variation is common at the Xce locus and we conclude that the observed allelic series is a product of independent and recurring duplications shared between weak Xce alleles.

Shared genetic risk between eating disorder- and substance-use-related phenotypes: Evidence from genome-wide association studies.

Eating disorders and substance use disorders frequently co-occur. Twin studies reveal shared genetic variance between liabilities to eating disorders and substance use, with the strongest associations between symptoms of bulimia nervosa and problem alcohol use (genetic correlation [rg ], twin-based = 0.23-0.53). We estimated the genetic correlation between eating disorder and substance use and disorder phenotypes using data from genome-wide association studies (GWAS). Four eating disorder phenotypes (anorexia nervosa [AN], AN with binge eating, AN without binge eating, and a bulimia nervosa factor score), and eight substance-use-related phenotypes (drinks per week, alcohol use disorder [AUD], smoking initiation, current smoking, cigarettes per day, nicotine dependence, cannabis initiation, and cannabis use disorder) from eight studies were included. Significant genetic correlations were adjusted for variants associated with major depressive disorder and schizophrenia. Total study sample sizes per phenotype ranged from ~2400 to ~537 000 individuals. We used linkage disequilibrium score regression to calculate single nucleotide polymorphism-based genetic correlations between eating disorder- and substance-use-related phenotypes. Significant positive genetic associations emerged between AUD and AN (rg = 0.18; false discovery rate q = 0.0006), cannabis initiation and AN (rg = 0.23; q < 0.0001), and cannabis initiation and AN with binge eating (rg = 0.27; q = 0.0016). Conversely, significant negative genetic correlations were observed between three nondiagnostic smoking phenotypes (smoking initiation, current smoking, and cigarettes per day) and AN without binge eating (rgs = -0.19 to -0.23; qs < 0.04). The genetic correlation between AUD and AN was no longer significant after co-varying for major depressive disorder loci. The patterns of association between eating disorder- and substance-use-related phenotypes highlights the potentially complex and substance-specific relationships among these behaviors.

Antipsychotic Behavioral Phenotypes in the Mouse Collaborative Cross Recombinant Inbred Inter-Crosses (RIX).

Schizophrenia is an idiopathic disorder that affects approximately 1% of the human population, and presents with persistent delusions, hallucinations, and disorganized behaviors. Antipsychotics are the standard pharmacological treatment for schizophrenia, but are frequently discontinued by patients due to inefficacy and/or side effects. Chronic treatment with the typical antipsychotic haloperidol causes tardive dyskinesia (TD), which manifests as involuntary and often irreversible orofacial movements in around 30% of patients. Mice treated with haloperidol develop many of the features of TD, including jaw tremors, tongue protrusions, and vacuous chewing movements (VCMs). In this study, we used genetically diverse Collaborative Cross (CC) recombinant inbred inter-cross (RIX) mice to elucidate the genetic basis of antipsychotic-induced adverse drug reactions (ADRs). We performed a battery of behavioral tests in 840 mice from 73 RIX lines (derived from 62 CC strains) treated with haloperidol or placebo in order to monitor the development of ADRs. We used linear mixed models to test for strain and treatment effects. We observed highly significant strain effects for almost all behavioral measurements investigated (P < 0.001). Further, we observed strong strain-by-treatment interactions for most phenotypes, particularly for changes in distance traveled, vertical activity, and extrapyramidal symptoms (EPS). Estimates of overall heritability ranged from 0.21 (change in body weight) to 0.4 (VCMs and change in distance traveled) while the portion attributable to the interactions of treatment and strain ranged from 0.01 (for change in body weight) to 0.15 (for change in EPS). Interestingly, close to 30% of RIX mice exhibited VCMs, a sensitivity to haloperidol exposure, approximately similar to the rate of TD in humans chronically exposed to haloperidol. Understanding the genetic basis for the susceptibility to antipsychotic ADRs may be possible in mouse, and extrapolation to humans could lead to safer therapeutic approaches for schizophrenia.

Robust Hi-C Maps of Enhancer-Promoter Interactions Reveal the Function of Non-coding Genome in Neural Development and Diseases.

Genome-wide mapping of chromatin interactions at high resolution remains experimentally and computationally challenging. Here we used a low-input "easy Hi-C" protocol to map the 3D genome architecture in human neurogenesis and brain tissues and also demonstrated that a rigorous Hi-C bias-correction pipeline (HiCorr) can significantly improve the sensitivity and robustness of Hi-C loop identification at sub-TAD level, especially the enhancer-promoter (E-P) interactions. We used HiCorr to compare the high-resolution maps of chromatin interactions from 10 tissue or cell types with a focus on neurogenesis and brain tissues. We found that dynamic chromatin loops are better hallmarks for cellular differentiation than compartment switching. HiCorr allowed direct observation of cell-type- and differentiation-specific E-P aggregates spanning large neighborhoods, suggesting a mechanism that stabilizes enhancer contacts during development. Interestingly, we concluded that Hi-C loop outperforms eQTL in explaining neurological GWAS results, revealing a unique value of high-resolution 3D genome maps in elucidating the disease etiology.

Increased burden of ultra-rare structural variants localizing to boundaries of topologically associated domains in schizophrenia.

Despite considerable progress in schizophrenia genetics, most findings have been for large rare structural variants and common variants in well-imputed regions with few genes implicated from exome sequencing. Whole genome sequencing (WGS) can potentially provide a more complete enumeration of etiological genetic variation apart from the exome and regions of high linkage disequilibrium. We analyze high-coverage WGS data from 1162 Swedish schizophrenia cases and 936 ancestry-matched population controls. Our main objective is to evaluate the contribution to schizophrenia etiology from a variety of genetic variants accessible to WGS but not by previous technologies. Our results suggest that ultra-rare structural variants that affect the boundaries of topologically associated domains (TADs) increase risk for schizophrenia. Alterations in TAD boundaries may lead to dysregulation of gene expression. Future mechanistic studies will be needed to determine the precise functional effects of these variants on biology.

Correction: Common-variant associations with fragile X syndrome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

LEGAL AND ETHICAL IMPLICATIONS OF CRISPR APPLICATIONS IN PSYCHIATRY.

Gene-environment interactions play a key role in how psychiatric disorders manifest and develop. Psychiatric genetics researchers are making progress in identifying genomic correlates of many disorders. And recently, the field of genetics has given rise to a technology that many claim will revolutionize the biological sciences and propel the field into a transformative phase: the powerful gene-editing tool known as CRISPR-Cas9. This Article illustrates which psychiatric conditions are likely to make an attractive target for CRISPR as the technology evolves and CRISPR therapies becomes a viable tool to manage or prevent disorders in a clinical setting. We examine the potential scientific and clinical challenges of applying CRISPR in the mental health context, along with the regulatory, ethical, and legal issues that might arise as a consequence of these applications.

Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa.

Characterized primarily by a low body-mass index, anorexia nervosa is a complex and serious illness1, affecting 0.9-4% of women and 0.3% of men2-4, with twin-based heritability estimates of 50-60%5. Mortality rates are higher than those in other psychiatric disorders6, and outcomes are unacceptably poor7. Here we combine data from the Anorexia Nervosa Genetics Initiative (ANGI)8,9 and the Eating Disorders Working Group of the Psychiatric Genomics Consortium (PGC-ED) and conduct a genome-wide association study of 16,992 cases of anorexia nervosa and 55,525 controls, identifying eight significant loci. The genetic architecture of anorexia nervosa mirrors its clinical presentation, showing significant genetic correlations with psychiatric disorders, physical activity, and metabolic (including glycemic), lipid and anthropometric traits, independent of the effects of common variants associated with body-mass index. These results further encourage a reconceptualization of anorexia nervosa as a metabo-psychiatric disorder. Elucidating the metabolic component is a critical direction for future research, and paying attention to both psychiatric and metabolic components may be key to improving outcomes.

Non-coding variability at the APOE locus contributes to the Alzheimer's risk.

Alzheimer's disease (AD) is a leading cause of mortality in the elderly. While the coding change of APOE-ε4 is a key risk factor for late-onset AD and has been believed to be the only risk factor in the APOE locus, it does not fully explain the risk effect conferred by the locus. Here, we report the identification of AD causal variants in PVRL2 and APOC1 regions in proximity to APOE and define common risk haplotypes independent of APOE-ε4 coding change. These risk haplotypes are associated with changes of AD-related endophenotypes including cognitive performance, and altered expression of APOE and its nearby genes in the human brain and blood. High-throughput genome-wide chromosome conformation capture analysis further supports the roles of these risk haplotypes in modulating chromatin states and gene expression in the brain. Our findings provide compelling evidence for additional risk factors in the APOE locus that contribute to AD pathogenesis.

Common-variant associations with fragile X syndrome.

Fragile X syndrome is rare but a prominent cause of intellectual disability. It is usually caused by a de novo mutation that occurs on multiple haplotypes and thus would not be expected to be detectible using genome-wide association (GWA). We conducted GWA in 89 male FXS cases and 266 male controls, and detected multiple genome-wide significant signals near FMR1 (odds ratio = 8.10, P = 2.5 × 10-10). These findings withstood robust attempts at falsification. Fine-mapping yielded a minimum P = 1.13 × 10-14, but did not narrow the interval. Comprehensive functional genomic integration did not provide a mechanistic hypothesis. Controls carrying a risk haplotype had significantly longer FMR1 CGG repeats than controls with the protective haplotype (P = 4.75 × 10-5), which may predispose toward increases in CGG number to the premutation range over many generations. This is a salutary reminder of the complexity of even "simple" monogenetic disorders.