The Genetic Architecture of the Human Corpus Callosum and its Subregions.

The corpus callosum (CC) is the largest set of white matter fibers connecting the two hemispheres of the brain. In humans, it is essential for coordinating sensorimotor responses, performing associative/executive functions, and representing information in multiple dimensions. Understanding which genetic variants underpin corpus callosum morphometry, and their shared influence on cortical structure and susceptibility to neuropsychiatric disorders, can provide molecular insights into the CC's role in mediating cortical development and its contribution to neuropsychiatric disease. To characterize the morphometry of the midsagittal corpus callosum, we developed a publicly available artificial intelligence based tool to extract, parcellate, and calculate its total and regional area and thickness. Using the UK Biobank (UKB) and the Adolescent Brain Cognitive Development study (ABCD), we extracted measures of midsagittal corpus callosum morphometry and performed a genome-wide association study (GWAS) meta-analysis of European participants (combined N = 46,685). We then examined evidence for generalization to the non-European participants of the UKB and ABCD cohorts (combined N = 7,040). Post-GWAS analyses implicate prenatal intracellular organization and cell growth patterns, and high heritability in regions of open chromatin, suggesting transcriptional activity regulation in early development. Results suggest programmed cell death mediated by the immune system drives the thinning of the posterior body and isthmus. Global and local genetic overlap, along with causal genetic liability, between the corpus callosum, cerebral cortex, and neuropsychiatric disorders such as attention-deficit/hyperactivity and bipolar disorders were identified. These results provide insight into variability of corpus callosum development, its genetic influence on the cerebral cortex, and biological mechanisms related to neuropsychiatric dysfunction.

Genome-wide meta-analyses of non-response to antidepressants identify novel loci and potential drugs.

Antidepressants exhibit a considerable variation in efficacy, and increasing evidence suggests that individual genetics contribute to antidepressant treatment response. Here, we combined data on antidepressant non-response measured using rating scales for depressive symptoms, questionnaires of treatment effect, and data from electronic health records, to increase statistical power to detect genomic loci associated with non-response to antidepressants in a total sample of 135,471 individuals prescribed antidepressants. We performed genome-wide association meta-analyses, leave-one-out polygenic prediction, and bioinformatics analyses for genetically informed drug prioritization. We identified two novel loci associated with non-response to antidepressants and showed significant polygenic prediction in independent samples. In addition, we investigated drugs that target proteins likely involved in mechanisms underlying antidepressant non-response, and shortlisted drugs that warrant further replication and validation of their potential to reduce depressive symptoms in individuals who do not respond to first-line antidepressant medications. These results suggest that meta-analyses of GWAS utilizing real-world measures of treatment outcomes can increase sample sizes to improve the discovery of variants associated with non-response to antidepressants.

Genetic analyses point to alterations in immune-related pathways underpinning the association between psychiatric disorders and COVID-19.

Current literature suggests that people with psychiatric disorders have a higher risk of COVID-19 infection and a worse prognosis of the disease. We aimed to study the genetic contribution to these associations across seven psychiatric disorders as well as a general psychopathology factor (P-factor) and determine whether these are unique or shared across psychiatric disorders using statistical genetic techniques. Using the largest available genome-wide association studies (GWAS), we found a significant genetic overlap between depression, ADHD, PTSD, and the P-factor with both COVID-19 infection and hospitalization, and between anxiety and COVID-19 hospitalization. We used pairwise GWAS to examine this overlap on a fine-grained scale and identified specific regions of the genome shared between several psychiatric disorders, the P-factor, and COVID-19. Gene-based analysis in these genomic regions suggested possible links with immune-related pathways such as thyroid homeostasis, inflammation, and stress response. Finally, we show preliminary evidence for causal associations between depression, ADHD, PTSD, and the P-factor, and higher COVID-19 infection and hospitalization using Mendelian Randomization and Latent Causal Variable methods. Our results support the hypothesis that the relationship between psychiatric disorders and COVID-19 risk is likely due to shared alterations in immune-related pathways and is not a result of environmental factors alone, shedding light on potentially viable therapeutic targets.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum.

Pubertal timing varies considerably and is associated with later health outcomes. We performed multi-ancestry genetic analyses on ~800,000 women, identifying 1,080 signals for age at menarche. Collectively, these explained 11% of trait variance in an independent sample. Women at the top and bottom 1% of polygenic risk exhibited ~11 and ~14-fold higher risks of delayed and precocious puberty, respectively. We identified several genes harboring rare loss-of-function variants in ~200,000 women, including variants in ZNF483, which abolished the impact of polygenic risk. Variant-to-gene mapping approaches and mouse gonadotropin-releasing hormone neuron RNA sequencing implicated 665 genes, including an uncharacterized G-protein-coupled receptor, GPR83, which amplified the signaling of MC3R, a key nutritional sensor. Shared signals with menopause timing at genes involved in DNA damage response suggest that the ovarian reserve might signal centrally to trigger puberty. We also highlight body size-dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Clozapine Efficacy and Adverse Drug Reactions Among a Nationwide Study of 1021 Australians Prescribed Clozapine: The ClozaGene Study.

BACKGROUND AND HYPOTHESIS: The ClozaGene Study is a nationwide cohort of adults who have been treated with clozapine. While clozapine is indicated in the management of treatment-resistant schizophrenia, it is associated with a considerable adverse drug reaction (ADR) burden, and not all patients achieve adequate symptomatic response. The current study focuses on self-reported experiences of clozapine use and response, clozapine-associated ADRs, and mental health comorbidity. STUDY DESIGN: A total of 1021 participants (41.0% female; aged 46.2 ±â€…10.6 years [range 18-66]) were recruited via a mail-out based on prescriptions for clozapine. Participants completed a self-report questionnaire. STUDY RESULTS: Most participants (90.1%, n = 912) were living with schizophrenia while 41.5% reported a lifetime diagnosis of depression, 15.6% bipolar disorder, and 8.1% schizoaffective disorder. Clozapine was currently prescribed to 944 (92.5%) participants and 37.8% of these participants self-reported currently taking additional antipsychotic medication. Nearly 3 quarters of participants living with schizophrenia reported that clozapine helped control their schizophrenia symptoms moderately to very well. The most commonly reported ADRs were sialorrhea (80.3%), weight gain (71.0%), constipation (56.9%), and sedation (52.8%). The prevalence of clozapine cessation due to clozapine-induced myocarditis and neutropenia was 1% and 0.4%, respectively. CONCLUSIONS: Our findings highlight the high rate of psychotic and metabolic symptoms and ADRs among adults prescribed clozapine in the general Australian population. Future genomic analyses will focus on identifying genetic variants influencing clozapine treatment response and side effects.

Heritability of cerebellar subregion volumes in adolescent and young adult twins.

Twin studies have found gross cerebellar volume to be highly heritable. However, whether fine-grained regional volumes within the cerebellum are similarly heritable is still being determined. Anatomical MRI scans from two independent datasets (QTIM: Queensland Twin IMaging, N = 798, mean age 22.1 years; QTAB: Queensland Twin Adolescent Brain, N = 396, mean age 11.3 years) were combined with an optimised and automated cerebellum parcellation algorithm to segment and measure 28 cerebellar regions. We show that the heritability of regional volumetric measures varies widely across the cerebellum ( h 2 $$ {h}^2 $$ 47%-91%). Additionally, the good to excellent test-retest reliability for a subsample of QTIM participants suggests that non-genetic variance in cerebellar volumes is due primarily to unique environmental influences rather than measurement error. We also show a consistent pattern of strong associations between the volumes of homologous left and right hemisphere regions. Associations were predominantly driven by genetic effects shared between lobules, with only sparse contributions from environmental effects. These findings are consistent with similar studies of the cerebrum and provide a first approximation of the upper bound of heritability detectable by genome-wide association studies.

Comprehensive Sex-Stratified Genetic Analysis of 28 Blood Biomarkers and Depression Reveals a Significant Association between Depression and Low Levels of Total Protein in Females.

Introduction: Major depression (MD) is more common amongst women than men, and MD episodes have been associated with fluctuations in reproductive hormones amongst women. To investigate biological underpinnings of heterogeneity in MD, the associations between depression, stratified by sex and including perinatal depression (PND), and blood biomarkers, using UK Biobank (UKB) data, were evaluated, and extended to include the association of depression with biomarker polygenic scores (PGS), generated as proxy for each biomarker. Method: Using female (N = 39,761) and male (N = 38,821) UKB participants, lifetime MD and PND were tested for association with 28 blood biomarkers. A GWAS was conducted for each biomarker and genetic correlations with depression subgroups were estimated. Using independent data from the Australian Genetics of Depression Study, PGS were constructed for each biomarker, and tested for association with depression status (n [female cases/controls] = 9,006/6,442; n [male cases/controls] = 3,106/6,222). Regions of significant local genetic correlation between depression subgroups and biomarkers highlighted by the PGS analysis were identified. Results: Depression in females was significantly associated with levels of twelve biomarkers, including total protein (OR = 0.90, CI = [0.86, 0.94], p = 3.9 × 10-6) and vitamin D (OR = 0.94, CI = [0.90, 0.97], p = 2.6 × 10-4), and PND with five biomarker levels, also including total protein (OR = 0.88, CI = [0.81, 0.96], p = 4.7 × 10-3). Depression in males was significantly associated with levels of eleven biomarkers. In the independent Australian Genetics of Depression Study, PGS analysis found significant associations for female depression and PND with total protein (female depression: OR = 0.93, CI = [0.88, 0.98], p = 3.6 × 10-3; PND: OR = 0.91, CI = [0.86, 0.96], p = 1.1 × 10-3), as well as with vitamin D (female depression: OR = 0.93, CI = [0.89, 0.97], p = 2.0 × 10-3; PND: OR = 0.92, CI = [0.87, 0.97], p = 1.4 × 10-3). The male depression sample did not report any significant results, and the point estimate of total protein (OR = 0.98, CI = [0.92-1.04], p = 4.7 × 10-1) did not indicate any association. Local genetic correlation analysis highlighted significant genetic correlation between PND and total protein, located in 5q13.3 (rG = 0.68, CI = [0.33, 1.0], p = 3.6 × 10-4). Discussion and Conclusion: Multiple lines of evidence from genetic analysis highlight an association between total serum protein levels and depression in females. Further research involving prospective measurement of total protein and depressive symptoms is warranted.

Patterns of stressful life events and polygenic scores for five mental disorders and neuroticism among adults with depression.

The dominant ('general') version of the diathesis-stress theory of depression views stressors and genetic vulnerability as independent risks. In the Australian Genetics of Depression Study (N = 14,146; 75% female), we tested whether polygenic scores (PGS) for major depression, bipolar disorder, schizophrenia, anxiety, ADHD, and neuroticism were associated with reported exposure to 32 childhood, past-year, lifetime, and accumulated stressful life events (SLEs). In false discovery rate-corrected models, the clearest PGS-SLE relationships were for the ADHD- and depression-PGSs, and to a lesser extent, the anxiety- and schizophrenia-PGSs. We describe the associations for childhood and accumulated SLEs, and the 2-3 strongest past-year/lifetime SLE associations. Higher ADHD-PGS was associated with all childhood SLEs (emotional abuse, emotional neglect, physical neglect; ORs = 1.09-1.14; p's < 1.3 × 10-5), more accumulated SLEs, and reported exposure to sudden violent death (OR = 1.23; p = 3.6 × 10-5), legal troubles (OR = 1.15; p = 0.003), and sudden accidental death (OR = 1.14; p = 0.006). Higher depression-PGS was associated with all childhood SLEs (ORs = 1.07-1.12; p's < 0.013), more accumulated SLEs, and severe human suffering (OR = 1.17; p = 0.003), assault with a weapon (OR = 1.12; p = 0.003), and living in unpleasant surroundings (OR = 1.11; p = 0.001). Higher anxiety-PGS was associated with childhood emotional abuse (OR = 1.08; p = 1.6 × 10-4), more accumulated SLEs, and serious accident (OR = 1.23; p = 0.004), physical assault (OR = 1.08; p = 2.2 × 10-4), and transportation accident (OR = 1.07; p = 0.001). Higher schizophrenia-PGS was associated with all childhood SLEs (ORs = 1.12-1.19; p's < 9.3-8), more accumulated SLEs, and severe human suffering (OR = 1.16; p = 0.003). Higher neuroticism-PGS was associated with living in unpleasant surroundings (OR = 1.09; p = 0.007) and major financial troubles (OR = 1.06; p = 0.014). A reversed pattern was seen for the bipolar-PGS, with lower odds of reported physical assault (OR = 0.95; p = 0.014), major financial troubles (OR = 0.93; p = 0.004), and living in unpleasant surroundings (OR = 0.92; p = 0.007). Genetic risk for several mental disorders influences reported exposure to SLEs among adults with moderately severe, recurrent depression. Our findings emphasise that stressors and diatheses are inter-dependent and challenge diagnosis and subtyping (e.g., reactive/endogenous) based on life events.

Genetic architecture of the structural connectome.

Myelinated axons form long-range connections that enable rapid communication between distant brain regions, but how genetics governs the strength and organization of these connections remains unclear. We perform genome-wide association studies of 206 structural connectivity measures derived from diffusion magnetic resonance imaging tractography of 26,333 UK Biobank participants, each representing the density of myelinated connections within or between a pair of cortical networks, subcortical structures or cortical hemispheres. We identify 30 independent genome-wide significant variants after Bonferroni correction for the number of measures studied (126 variants at nominal genome-wide significance) implicating genes involved in myelination (SEMA3A), neurite elongation and guidance (NUAK1, STRN, DPYSL2, EPHA3, SEMA3A, HGF, SHTN1), neural cell proliferation and differentiation (GMNC, CELF4, HGF), neuronal migration (CCDC88C), cytoskeletal organization (CTTNBP2, MAPT, DAAM1, MYO16, PLEC), and brain metal transport (SLC39A8). These variants have four broad patterns of spatial association with structural connectivity: some have disproportionately strong associations with corticothalamic connectivity, interhemispheric connectivity, or both, while others are more spatially diffuse. Structural connectivity measures are highly polygenic, with a median of 9.1 percent of common variants estimated to have non-zero effects on each measure, and exhibited signatures of negative selection. Structural connectivity measures have significant genetic correlations with a variety of neuropsychiatric and cognitive traits, indicating that connectivity-altering variants tend to influence brain health and cognitive function. Heritability is enriched in regions with increased chromatin accessibility in adult oligodendrocytes (as well as microglia, inhibitory neurons and astrocytes) and multiple fetal cell types, suggesting that genetic control of structural connectivity is partially mediated by effects on myelination and early brain development. Our results indicate pervasive, pleiotropic, and spatially structured genetic control of white-matter structural connectivity via diverse neurodevelopmental pathways, and support the relevance of this genetic control to healthy brain function.

Lifespan reference curves for harmonizing multi-site regional brain white matter metrics from diffusion MRI.

Age-related white matter (WM) microstructure maturation and decline occur throughout the human lifespan, complementing the process of gray matter development and degeneration. Here, we create normative lifespan reference curves for global and regional WM microstructure by harmonizing diffusion MRI (dMRI)-derived data from ten public datasets (N = 40,898 subjects; age: 3-95 years; 47.6% male). We tested three harmonization methods on regional diffusion tensor imaging (DTI) based fractional anisotropy (FA), a metric of WM microstructure, extracted using the ENIGMA-DTI pipeline. ComBat-GAM harmonization provided multi-study trajectories most consistent with known WM maturation peaks. Lifespan FA reference curves were validated with test-retest data and used to assess the effect of the ApoE4 risk factor for dementia in WM across the lifespan. We found significant associations between ApoE4 and FA in WM regions associated with neurodegenerative disease even in healthy individuals across the lifespan, with regional age-by-genotype interactions. Our lifespan reference curves and tools to harmonize new dMRI data to the curves are publicly available as eHarmonize (https://github.com/ahzhu/eharmonize).

Identifying genetic differences between bipolar disorder and major depression through multiple GWAS.

Background: Accurate diagnosis of bipolar disorder (BD) is difficult in clinical practice, with an average delay between symptom onset and diagnosis of about 7 years. A key reason is that the first manic episode is often preceded by a depressive one, making it difficult to distinguish BD from unipolar major depressive disorder (MDD). Aims: Here, we use genome-wide association analyses (GWAS) to identify differential genetic factors and to develop predictors based on polygenic risk scores that may aid early differential diagnosis. Methods: Based on individual genotypes from case-control cohorts of BD and MDD shared through the Psychiatric Genomics Consortium, we compile case-case-control cohorts, applying a careful merging and quality control procedure. In a resulting cohort of 51,149 individuals (15,532 BD cases, 12,920 MDD cases and 22,697 controls), we perform a variety of GWAS and polygenic risk scores (PRS) analyses. Results: While our GWAS is not well-powered to identify genome-wide significant loci, we find significant SNP-heritability and demonstrate the ability of the resulting PRS to distinguish BD from MDD, including BD cases with depressive onset. We replicate our PRS findings, but not signals of individual loci in an independent Danish cohort (iPSYCH 2015 case-cohort study, N=25,966). We observe strong genetic correlation between our case-case GWAS and that of case-control BD. Conclusions: We find that MDD and BD, including BD with a depressive onset, are genetically distinct. Further, our findings support the hypothesis that Controls - MDD - BD primarily lie on a continuum of genetic risk. Future studies with larger and richer samples will likely yield a better understanding of these findings and enable the development of better genetic predictors distinguishing BD and, importantly, BD with depressive onset from MDD.

Evening Chronotypes With Depression Report Poorer Outcomes of Selective Serotonin Reuptake Inhibitors: A Survey-Based Study of Self-Ratings.

BACKGROUND: Preliminary evidence suggests that evening chronotype is related to poorer efficacy of selective serotonin reuptake inhibitors. It is unknown whether this is specific to particular medications, self-rated chronotype, or efficacy. METHODS: In the Australian Genetics of Depression Study (n = 15,108; 75% women; 18-90 years; 68% with ≥1 other lifetime diagnosis), a survey recorded experiences with 10 antidepressants, and the reduced Morningness-Eveningness Questionnaire was used to estimate chronotype. A chronotype polygenic score was calculated. Age- and sex-adjusted regression models (Bonferroni-corrected) estimated associations among antidepressant variables (how well the antidepressant worked [efficacy], duration of symptom improvement, side effects, discontinuation due to side effects) and self-rated and genetic chronotypes. RESULTS: The chronotype polygenic score explained 4% of the variance in self-rated chronotype (r = 0.21). Higher self-rated eveningness was associated with poorer efficacy of escitalopram (odds ratio [OR] = 1.04; 95% CI, 1.02 to 1.06; p = .000035), citalopram (OR = 1.03; 95% CI, 1.01 to 1.05; p = .004), fluoxetine (OR = 1.03; 95% CI, 1.01 to 1.05; p = .001), sertraline (OR = 1.02; 95% CI, 1.01 to 1.04; p = .0008), and desvenlafaxine (OR = 1.03; 95% CI, 1.01 to 1.05; p = .004), and a profile of increased side effects (80% of those recorded; ORs = 0.93-0.98), with difficulty getting to sleep the most common. Self-rated chronotype was unrelated to duration of improvement or discontinuation. The chronotype polygenic score was only associated with suicidal thoughts and attempted suicide (self-reported). While our measures are imperfect, and not of circadian phase under controlled conditions, the model coefficients suggest that dysregulation of the phenotypic chronotype relative to its genetic proxy drove relationships with antidepressant outcomes. CONCLUSIONS: The idea that variation in circadian factors influences response to antidepressants was supported and encourages exploration of circadian mechanisms of depressive disorders and antidepressant treatments.

Scientific clickbait: Examining media coverage and readability in genome-wide association research.

In the present study, we analyzed a large corpus of English-language online media articles covering genome-wide association studies (GWAS), exemplifying the use of computational methods to study science communication in biological sciences. We analyzed trends in media coverage, readability, themes, and mentions of ethical and social issues, in over 5,000 websites published from 2005 to 2018 from 3,555 GWAS publications on 1,943 different traits, identified via GWAS Catalog using a text-mining approach to inform the discussion about genetic literacy and media coverage. We found that 22.9% of GWAS papers received media attention but most were described in language too complex to be understood by the public. Ethical issues are rarely mentioned and mentions of translation are increasing over time. We predicted media attention based on year of publication, number of genetic associations identified, study sample size, and journal impact factor, using a regression model (r2 = 38.7%). We found that chronotype, educational attainment, alcohol and coffee consumption, sexual orientation, tanning, and hair color received substantially more attention than predicted by the regression model. We also evaluated the prevalence of the clickbait "one gene, one disease" headlines (e.g., "Scientists Say They've Found Gene That Causes Breast Cancer") and found that it is declining. In sum, online media coverage of GWAS should be more accessible, introduce more modern genetics terms, and when appropriate, ELSI should be mentioned. Science communication research can benefit from big data and text-mining techniques which allow us to study trends and changes in coverage trends across thousands of media outlets. Results can be explored interactively in a website we have built for this manuscript: https://jjmorosoli.shinyapps.io/newas/.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

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.

Discontinuation of antidepressant treatment: a retrospective cohort study on more than 20,000 participants.

BACKGROUND: Factors influencing antidepressant treatment discontinuation are poorly understood. In the present study, we aimed to estimate the prevalence of antidepressant treatment discontinuation and identify demographic characteristics, psychiatric comorbidities, and specific side effects associated with treatment discontinuation. METHODS: We leveraged data from the Australian Genetics of Depression Study (AGDS; N = 20,941) to perform a retrospective cohort study on antidepressant treatment discontinuation. Participants were eligible if they were over 18 years of age, had taken antidepressants in the past 4 years, and provided informed consent. RESULTS: Among the ten antidepressants studied, the highest discontinuation rates were observed for Mirtazapine (57.3%) and Amitriptyline (51.6%). Discontinuation rates were comparable across sexes except for Mirtazapine, for which women were more likely to discontinue. The two most common side effects, reduced sexual function and weight gain, were not associated with increased odds of treatment discontinuation. Anxiety, agitation, suicidal thoughts, vomiting, and rashes were associated with higher odds for treatment discontinuation, as were lifetime diagnoses of PTSD, ADHD, and a higher neuroticism score. Educational attainment showed a negative (protective) association with discontinuation across medications. CONCLUSIONS: Our study suggests that not all side effects contribute equally to discontinuation. Common side effects such as reduced sexual function and weight gain may not necessarily increase the risk of treatment discontinuation. Side effects linked to discontinuation can be divided into two groups, psychopathology related and allergy/intolerance.

GWAS of Dizygotic Twinning in an Enlarged Australian Sample of Mothers of DZ Twins.

Female fertility is a complex trait with age-specific changes in spontaneous dizygotic (DZ) twinning and fertility. To elucidate factors regulating female fertility and infertility, we conducted a genome-wide association study (GWAS) on mothers of spontaneous DZ twins (MoDZT) versus controls (3273 cases, 24,009 controls). This is a follow-up study to the Australia/New Zealand (ANZ) component of that previously reported (Mbarek et al., 2016), with a sample size almost twice that of the entire discovery sample meta-analysed in the previous article (and five times the ANZ contribution to that), resulting from newly available additional genotyping and representing a significant increase in power. We compare analyses with and without male controls and show unequivocally that it is better to include male controls who have been screened for recent family history, than to use only female controls. Results from the SNP based GWAS identified four genomewide significant signals, including one novel region, ZFPM1 (Zinc Finger Protein, FOG Family Member 1), on chromosome 16. Previous signals near FSHB (Follicle Stimulating Hormone beta subunit) and SMAD3 (SMAD Family Member 3) were also replicated (Mbarek et al., 2016). We also ran the GWAS with a dominance model that identified a further locus ADRB2 on chr 5. These results have been contributed to the International Twinning Genetics Consortium for inclusion in the next GWAS meta-analysis (Mbarek et al., in press).

Uncovering the complex relationship between balding, testosterone and skin cancers in men.

Male-pattern baldness (MPB) is related to dysregulation of androgens such as testosterone. A previously observed relationship between MPB and skin cancer may be due to greater exposure to ultraviolet radiation or indicate a role for androgenic pathways in the pathogenesis of skin cancers. We dissected this relationship via Mendelian randomization (MR) analyses, using genetic data from recent male-only meta-analyses of cutaneous melanoma (12,232 cases; 20,566 controls) and keratinocyte cancers (KCs) (up to 17,512 cases; >100,000 controls), followed by stratified MR analysis by body-sites. We found strong associations between MPB and the risk of KC, but not with androgens, and multivariable models revealed that this relationship was heavily confounded by MPB single nucleotide polymorphisms involved in pigmentation pathways. Site-stratified MR analyses revealed strong associations between MPB with head and neck squamous cell carcinoma and melanoma, suggesting that sun exposure on the scalp, rather than androgens, is the main driver. Men with less hair covering likely explains, at least in part, the higher incidence of melanoma in men residing in countries with high ambient UV.

Preliminary results from the Australian Genetics of Bipolar Disorder Study: A nation-wide cohort.

OBJECTIVE: The Australian Genetics of Bipolar Disorder Study is a nation-wide cohort of adults living with bipolar disorder. The study aims to detect the relationships between genetic risk, symptom severity, and the lifetime prevalence of bipolar disorder, treatment response and medication side effects, and patterns and costs of health care usage. METHODS: A total of 6682 participants (68.3% female; aged 44.8 ± 13.6 years [range = 18-90]) were recruited in three waves: a nation-wide media campaign, a mail-out based on prescriptions for lithium carbonate and through the Australian Genetics of Depression Study. Participants completed a self-report questionnaire. A total of 4706 (70%) participants provided a saliva sample and were genotyped and 5506 (82%) consented to record linkage of their Pharmaceutical and Medicare Benefits Schedule data. RESULTS: Most participants were living with bipolar I disorder (n = 4068) while 1622 participants were living with bipolar II disorder and 992 with sub-threshold bipolar disorder. The mean age of bipolar disorder diagnosis was 32.7 ± 11.6 years but was younger in bipolar I (p = 2.0E-26) and females (p = 5.7E-23). Excluding depression with onset prior to bipolar disorder diagnosis, 64.5% of participants reported one or more co-occurring psychiatric disorders: most commonly generalised anxiety disorder (43.5%) and posttraumatic stress disorder (20.7%). Adverse drug reactions were common and resulted in discontinuation rates ranging from 33.4% for lithium to 63.0% for carbamazepine. CONCLUSION: Our findings highlight the high rate of comorbidities and adverse drug reactions among adults living with bipolar disorder in the general Australian population. Future genomic analyses focus on identifying genetic variants influencing pharmacotherapy treatment response and side effects.