Predicting ADHD in alcohol dependence using polygenic risk scores for ADHD

Attention deficit hyperactivity disorder (ADHD) is a neurodevelopmental disorder with a high degree of comorbidity, including substance misuse. We aimed to assess whether ADHD polygenic risk scores (PRS) could predict ADHD diagnosis in alcohol dependence (AD). ADHD PRS were generated for 1223 AD subjects with ADHD diagnosis information and 1818 healthy controls. ADHD PRS distributions were compared to evaluate the differences between healthy controls and AD cases with and without ADHD. We found increased ADHD PRS means in the AD cohort with ADHD (mean 0.30, standard deviation (SD) 0.92; p = 3.9 × 10-6); and without ADHD (mean - 0.00, SD 1.00; p = 5.2 × 10-5) compared to the healthy control subjects (mean - 0.17, SD 0.99). The ADHD PRS means differed within the AD group with a higher ADHD PRS mean in those with ADHD, odds ratio (OR) 1.34, confidence interval (CI) 1.10 to 1.65; p = 0.002. This study showed a positive relationship between ADHD PRS and risk of ADHD in individuals with co-occurring AD indicating that ADHD PRS may have utility in identifying individuals that are at a higher or lower risk of ADHD. Further larger studies need to be conducted to confirm the reliability of the results before ADHD PRS can be considered as a robust biomarker for diagnosis.

GWAS meta-analysis reveals key risk loci in essential tremor pathogenesis

Essential tremor (ET) is a prevalent neurological disorder with a largely unknown underlying biology. In this genome-wide association study meta-analysis, comprising 16,480 ET cases and 1,936,173 controls from seven datasets, we identify 12 sequence variants at 11 loci. Evaluating mRNA expression, splicing, plasma protein levels, and coding effects, we highlight seven putative causal genes at these loci, including CA3 and CPLX1. CA3 encodes Carbonic Anhydrase III and carbonic anhydrase inhibitors have been shown to decrease tremors. CPLX1, encoding Complexin-1, regulates neurotransmitter release. Through gene-set enrichment analysis, we identify a significant association with specific cell types, including dopaminergic and GABAergic neurons, as well as biological processes like Rho GTPase signaling. Genetic correlation analyses reveals a positive association between ET and Parkinson's disease, depression, and anxiety-related phenotypes. This research uncovers risk loci, enhancing our knowledge of the complex genetics of this common but poorly understood disorder, and highlights CA3 and CPLX1 as potential therapeutic targets.

Multi-polygenic scores in psychiatry: From disorder specific to transdiagnostic perspectives

The dense co-occurrence of psychiatric disorders questions the categorical classification tradition and motivates efforts to establish dimensional constructs with neurobiological foundations that transcend diagnostic boundaries. In this study, we examined the genetic liability for eight major psychiatric disorder phenotypes under both a disorder-specific and a transdiagnostic framework. The study sample (n = 513) was deeply phenotyped, consisting of 452 patients from tertiary care with mood disorders, anxiety disorders (ANX), attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorders, and/or substance use disorders (SUD) and 61 unaffected comparison individuals. We computed subject-specific polygenic risk score (PRS) profiles and assessed their associations with psychiatric diagnoses, comorbidity status, as well as cross-disorder behavioral dimensions derived from a rich battery of psychopathology assessments. High PRSs for depression were unselectively associated with the diagnosis of SUD, ADHD, ANX, and mood disorders (p < 1e-4). In the dimensional approach, four distinct functional domains were uncovered, namely the negative valence, social, cognitive, and regulatory systems, closely matching the major functional domains proposed by the Research Domain Criteria (RDoC) framework. Critically, the genetic predisposition for depression was selectively reflected in the functional aspect of negative valence systems (R2  = 0.041, p = 5e-4) but not others. This study adds evidence to the ongoing discussion about the misalignment between current psychiatric nosology and the underlying psychiatric genetic etiology and underscores the effectiveness of the dimensional approach in both the functional characterization of psychiatric patients and the delineation of the genetic liability for psychiatric disorders.

Rare variants with large effects provide functional insights into the pathology of migraine subtypes, with and without aura

Migraine is a complex neurovascular disease with a range of severity and symptoms, yet mostly studied as one phenotype in genome-wide association studies (GWAS). Here we combine large GWAS datasets from six European populations to study the main migraine subtypes, migraine with aura (MA) and migraine without aura (MO). We identified four new MA-associated variants (in PRRT2, PALMD, ABO and LRRK2) and classified 13 MO-associated variants. Rare variants with large effects highlight three genes. A rare frameshift variant in brain-expressed PRRT2 confers large risk of MA and epilepsy, but not MO. A burden test of rare loss-of-function variants in SCN11A, encoding a neuron-expressed sodium channel with a key role in pain sensation, shows strong protection against migraine. Finally, a rare variant with cis-regulatory effects on KCNK5 confers large protection against migraine and brain aneurysms. Our findings offer new insights with therapeutic potential into the complex biology of migraine and its subtypes.

Depression pathophysiology, risk prediction of recurrence and comorbid psychiatric disorders using genome-wide analyses

Depression is a common psychiatric disorder and a leading cause of disability worldwide. Here we conducted a genome-wide association study meta-analysis of six datasets, including >1.3 million individuals (371,184 with depression) and identified 243 risk loci. Overall, 64 loci were new, including genes encoding glutamate and GABA receptors, which are targets for antidepressant drugs. Intersection with functional genomics data prioritized likely causal genes and revealed new enrichment of prenatal GABAergic neurons, astrocytes and oligodendrocyte lineages. We found depression to be highly polygenic, with ~11,700 variants explaining 90% of the single-nucleotide polymorphism heritability, estimating that >95% of risk variants for other psychiatric disorders (anxiety, schizophrenia, bipolar disorder and attention deficit hyperactivity disorder) were influencing depression risk when both concordant and discordant variants were considered, and nearly all depression risk variants influenced educational attainment. Additionally, depression genetic risk was associated with impaired complex cognition domains. We dissected the genetic and clinical heterogeneity, revealing distinct polygenic architectures across subgroups of depression and demonstrating significantly increased absolute risks for recurrence and psychiatric comorbidity among cases of depression with the highest polygenic burden, with considerable sex differences. The risks were up to 5- and 32-fold higher than cases with the lowest polygenic burden and the background population, respectively. These results deepen the understanding of the biology underlying depression, its disease progression and inform precision medicine approaches to treatment.

Deficit of homozygosity among 1.52 million individuals and genetic causes of recessive lethality

Genotypes causing pregnancy loss and perinatal mortality are depleted among living individuals and are therefore difficult to find. To explore genetic causes of recessive lethality, we searched for sequence variants with deficit of homozygosity among 1.52 million individuals from six European populations. In this study, we identified 25 genes harboring protein-altering sequence variants with a strong deficit of homozygosity (10% or less of predicted homozygotes). Sequence variants in 12 of the genes cause Mendelian disease under a recessive mode of inheritance, two under a dominant mode, but variants in the remaining 11 have not been reported to cause disease. Sequence variants with a strong deficit of homozygosity are over-represented among genes essential for growth of human cell lines and genes orthologous to mouse genes known to affect viability. The function of these genes gives insight into the genetics of intrauterine lethality. We also identified 1077 genes with homozygous predicted loss-of-function genotypes not previously described, bringing the total set of genes completely knocked out in humans to 4785.

DanMAC5: a browser of aggregated sequence variants from 8,671 whole genome sequenced Danish individuals

OBJECTIVES

Allele counts of sequence variants obtained by whole genome sequencing (WGS) often play a central role in interpreting the results of genetic and genomic research. However, such variant counts are not readily available for individuals in the Danish population. Here, we present a dataset with allele counts for sequence variants (single nucleotide variants (SNVs) and indels) identified from WGS of 8,671 (5,418 females) individuals from the Danish population. The data resource is based on WGS data from three independent research projects aimed at assessing genetic risk factors for cardiovascular, psychiatric, and headache disorders. To enable the sharing of information on sequence variation in Danish individuals, we created summarized statistics on allele counts from anonymized data and made them available through the European Genome-phenome Archive (EGA, https://identifiers.org/ega.

DATASET

EGAD00001009756 ) and in a dedicated browser, DanMAC5 (available at www.danmac5.dk ). The summary level data and the DanMAC5 browser provide insight into the allelic spectrum of sequence variants segregating in the Danish population, which is important in variant interpretation.

DATA DESCRIPTION

Three WGS datasets with an average coverage of 30x were processed independently using the same quality control pipeline. Subsequently, we summarized, filtered, and merged allele counts to create a high-quality summary level dataset of sequence variants.

Genome-wide analyses of ADHD identify 27 risk loci, refine the genetic architecture and implicate several cognitive domains

Attention-deficit hyperactivity disorder (ADHD) is a prevalent neurodevelopmental disorder with a major genetic component. Here, we present a genome-wide association study meta-analysis of ADHD comprising 38,691 individuals with ADHD and 186,843 controls. We identified 27 genome-wide significant loci, highlighting 76 potential risk genes enriched among genes expressed particularly in early brain development. Overall, ADHD genetic risk was associated with several brain-specific neuronal subtypes and midbrain dopaminergic neurons. In exome-sequencing data from 17,896 individuals, we identified an increased load of rare protein-truncating variants in ADHD for a set of risk genes enriched with probable causal common variants, potentially implicating SORCS3 in ADHD by both common and rare variants. Bivariate Gaussian mixture modeling estimated that 84-98% of ADHD-influencing variants are shared with other psychiatric disorders. In addition, common-variant ADHD risk was associated with impaired complex cognition such as verbal reasoning and a range of executive functions, including attention.

A genome-wide meta-analysis identifies 50 genetic loci associated with carpal tunnel syndrome

Carpal tunnel syndrome (CTS) is the most common entrapment neuropathy and has a largely unknown underlying biology. In a genome-wide association study of CTS (48,843 cases and 1,190,837 controls), we found 53 sequence variants at 50 loci associated with the syndrome. The most significant association is with a missense variant (p.Glu366Lys) in SERPINA1 that protects against CTS (P = 2.9 × 10^-24, OR = 0.76). Through various functional analyses, we conclude that at least 22 genes mediate CTS risk and highlight the role of 19 CTS variants in the biology of the extracellular matrix. We show that the genetic component to the risk is higher in bilateral/recurrent/persistent cases than nonrecurrent/nonpersistent cases. Anthropometric traits including height and BMI are genetically correlated with CTS, in addition to early hormonal-replacement therapy, osteoarthritis, and restlessness. Our findings suggest that the components of the extracellular matrix play a key role in the pathogenesis of CTS.

Rare SLC13A1 variants associate with intervertebral disc disorder highlighting role of sulfate in disc pathology

Back pain is a common and debilitating disorder with largely unknown underlying biology. Here we report a genome-wide association study of back pain using diagnoses assigned in clinical practice; dorsalgia (119,100 cases, 909,847 controls) and intervertebral disc disorder (IDD) (58,854 cases, 922,958 controls). We identify 41 variants at 33 loci. The most significant association (OR_IDD = 0.92, P = 1.6 × 10⁻³⁹; OR_dorsalgia = 0.92, P = 7.2 × 10⁻¹⁵) is with a 3’UTR variant (rs1871452-T) in CHST3, encoding a sulfotransferase enzyme expressed in intervertebral discs. The largest effects on IDD are conferred by rare (MAF = 0.07 − 0.32%) loss-of-function (LoF) variants in SLC13A1, encoding a sodium-sulfate co-transporter (LoF burden OR = 1.44, P = 3.1 × 10⁻¹¹); variants that also associate with reduced serum sulfate. Genes implicated by this study are involved in cartilage and bone biology, as well as neurological and inflammatory processes.

Little is known about the biology of back pain, a leading cause of disability. Here the authors report 30 new back pain loci, implicating genes involved in cartilage/bone biology, as well as neurological and inflammatory processes.

Effects of copy number variations on brain structure and risk for psychiatric illness: Large-scale studies from the ENIGMA working groups on CNVs

The Enhancing NeuroImaging Genetics through Meta-Analysis copy number variant (ENIGMA-CNV) and 22q11.2 Deletion Syndrome Working Groups (22q-ENIGMA WGs) were created to gain insight into the involvement of genetic factors in human brain development and related cognitive, psychiatric and behavioral manifestations. To that end, the ENIGMA-CNV WG has collated CNV and magnetic resonance imaging (MRI) data from ~49,000 individuals across 38 global research sites, yielding one of the largest studies to date on the effects of CNVs on brain structures in the general population. The 22q-ENIGMA WG includes 12 international research centers that assessed over 533 individuals with a confirmed 22q11.2 deletion syndrome, 40 with 22q11.2 duplications, and 333 typically developing controls, creating the largest-ever 22q11.2 CNV neuroimaging data set. In this review, we outline the ENIGMA infrastructure and procedures for multi-site analysis of CNVs and MRI data. So far, ENIGMA has identified effects of the 22q11.2, 16p11.2 distal, 15q11.2, and 1q21.1 distal CNVs on subcortical and cortical brain structures. Each CNV is associated with differences in cognitive, neurodevelopmental and neuropsychiatric traits, with characteristic patterns of brain structural abnormalities. Evidence of gene-dosage effects on distinct brain regions also emerged, providing further insight into genotype-phenotype relationships. Taken together, these results offer a more comprehensive picture of molecular mechanisms involved in typical and atypical brain development. This "genotype-first" approach also contributes to our understanding of the etiopathogenesis of brain disorders. Finally, we outline future directions to better understand effects of CNVs on brain structure and behavior.

Thirty novel sequence variants impacting human intracranial volume

Intracranial volume, measured through magnetic resonance imaging and/or estimated from head circumference, is heritable and correlates with cognitive traits and several neurological disorders. We performed a genome-wide association study meta-analysis of intracranial volume (n = 79 174) and found 64 associating sequence variants explaining 5.0% of its variance. We used coding variation, transcript and protein levels, to uncover 12 genes likely mediating the effect of these variants, including GLI3 and CDK6 that affect cranial synostosis and microcephaly, respectively. Intracranial volume correlates genetically with volumes of cortical and sub-cortical regions, cognition, learning, neonatal and neurological traits. Parkinson's disease cases have greater and attention deficit hyperactivity disorder cases smaller intracranial volume than controls. Our Mendelian randomization studies indicate that intracranial volume associated variants either increase the risk of Parkinson's disease and decrease the risk of attention deficit hyperactivity disorder and neuroticism or correlate closely with a confounder.

Cohort profile: Copenhagen Hospital Biobank - Cardiovascular Disease Cohort (CHB-CVDC): Construction of a large-scale genetic cohort to facilitate a better understanding of heart diseases

PURPOSE

The aim of Copenhagen Hospital Biobank-Cardiovascular Disease Cohort (CHB-CVDC) is to establish a cohort that can accelerate our understanding of CVD initiation and progression by jointly studying genetics, diagnoses, treatments and risk factors.

PARTICIPANTS

The CHB-CVDC is a large genomic cohort of patients with CVD. CHB-CVDC currently includes 96 308 patients. The cohort is part of CHB initiated in 2009 in the Capital Region of Denmark. CHB is continuously growing with ~40 000 samples/year. Patients in CHB were included in CHB-CVDC if they were above 18 years of age and assigned at least one cardiovascular diagnosis. Additionally, up-to 110 000 blood donors can be analysed jointly with CHB-CVDC. Linkage with the Danish National Health Registries, Electronic Patient Records, and Clinical Quality Databases allow up-to 41 years of medical history. All individuals are genotyped using the Infinium Global Screening Array from Illumina and imputed using a reference panel consisting of whole-genome sequence data from 8429 Danes along with 7146 samples from North-Western Europe. Currently, 39 539 of the patients are deceased.

FINDINGS TO DATE

Here, we demonstrate the utility of the cohort by showing concordant effects between known variants and selected CVDs, that is, >93% concordance for coronary artery disease, atrial fibrillation, heart failure and cholesterol measurements and 85% concordance for hypertension. Furthermore, we evaluated multiple study designs and the validity of using Danish blood donors as part of CHB-CVDC. Lastly, CHB-CVDC has already made major contributions to studies of sick sinus syndrome and the role of phytosterols in development of atherosclerosis.

FUTURE PLANS

In addition to genetics, electronic patient records, national socioeconomic and health registries extensively characterise each patient in CHB-CVDC and provides a promising framework for improved understanding of risk and protective variants. We aim to include other measurable biomarkers for example, proteins in CHB-CVDC making it a platform for multiomics cardiovascular studies.

A genome-wide meta-analysis uncovers six sequence variants conferring risk of vertigo

Vertigo is the leading symptom of vestibular disorders and a major risk factor for falls. In a genome-wide association study of vertigo (Ncases = 48,072, Ncontrols = 894,541), we uncovered an association with six common sequence variants in individuals of European ancestry, including missense variants in ZNF91, OTOG, OTOGL, and TECTA, and a cis-eQTL for ARMC9. The association of variants in ZNF91, OTOGL, and OTOP1 was driven by an association with benign paroxysmal positional vertigo. Using previous reports of sequence variants associating with age-related hearing impairment and motion sickness, we found eight additional variants that associate with vertigo. Although disorders of the auditory and the vestibular system may co-occur, none of the six genome-wide significant vertigo variants were associated with hearing loss and only one was associated with age-related hearing impairment. Our results uncovered sequence variants associating with vertigo in a genome-wide association study and implicated genes with known roles in inner ear development, maintenance, and disease.

Analysis of Diffusion Tensor Imaging Data From the UK Biobank Confirms Dosage Effect of 15q11.2 Copy Number Variation on White Matter and Shows Association With Cognition

BACKGROUND

Copy number variations at the 15q11.2 BP1-BP2 locus are present in 0.5%-1.0% of the population, and the deletion is associated with several neurodevelopmental disorders. Previously, we showed a reciprocal effect of 15q11.2 copy number variation on fractional anisotropy, with widespread increases in deletion carriers. We aim to expand these findings using a larger sample of participants (N = 29,166) and higher resolution imaging and by examining the implications for cognitive performance.

METHODS

Diffusion tensor imaging measures from participants with no neurological or psychiatric diagnoses were obtained from the UK Biobank database. We compared 15q11.2 BP1-BP2 deletion (n = 102) and duplication (n = 113) carriers to a large cohort of control individuals with no neuropsychiatric copy number variants (n = 28,951). Additionally, we assessed how changes in white matter mediated the association between carrier status and cognitive performance.

RESULTS

Deletion carriers showed increases in fractional anisotropy in the internal capsule and cingulum and decreases in the posterior thalamic radiation compared with both duplication carriers and control subjects (who had intermediate values). Compared with control subjects, deletion carriers had lower scores across cognitive tasks, which were partly influenced by white matter. Reduced fractional anisotropy in the posterior thalamic radiation partially contributed to worse cognitive performance in deletion carriers.

CONCLUSIONS

These results, together with our previous findings, provide convergent evidence for an effect of 15q11.2 BP1-BP2 on white matter microstructure, this being more pronounced in deletion carriers. Additionally, changes in white matter were found to partially mediate cognitive ability in deletion carriers, providing a link between white matter changes in 15q11.2 BP1-BP2 carriers and cognitive function.

Identification of genetic overlap and novel risk loci for attention-deficit/hyperactivity disorder and bipolar disorder

Differential diagnosis between childhood onset attention-deficit/hyperactivity disorder (ADHD) and bipolar disorder (BD) remains a challenge, mainly due to overlapping symptoms and high rates of comorbidity. Despite this, genetic correlation reported for these disorders is low and non-significant. Here we aimed to better characterize the genetic architecture of these disorders utilizing recent large genome-wide association studies (GWAS). We analyzed independent GWAS summary statistics for ADHD (19,099 cases and 34,194 controls) and BD (20,352 cases and 31,358 controls) applying the conditional/conjunctional false discovery rate (condFDR/conjFDR) statistical framework that increases the power to detect novel phenotype-specific and shared loci by leveraging the combined power of two GWAS. We observed cross-trait polygenic enrichment for ADHD conditioned on associations with BD, and vice versa. Leveraging this enrichment, we identified 19 novel ADHD risk loci and 40 novel BD risk loci at condFDR <0.05. Further, we identified five loci jointly associated with ADHD and BD (conjFDR < 0.05). Interestingly, these five loci show concordant directions of effect for ADHD and BD. These results highlight a shared underlying genetic risk for ADHD and BD which may help to explain the high comorbidity rates and difficulties in differentiating between ADHD and BD in the clinic. Improving our understanding of the underlying genetic architecture of these disorders may aid in the development of novel stratification tools to help reduce these diagnostic difficulties.

Germline variants at SOHLH2 influence multiple myeloma risk

Multiple myeloma (MM) is caused by the uncontrolled, clonal expansion of plasma cells. While there is epidemiological evidence for inherited susceptibility, the molecular basis remains incompletely understood. We report a genome-wide association study totalling 5,320 cases and 422,289 controls from four Nordic populations, and find a novel MM risk variant at SOHLH2 at 13q13.3 (risk allele frequency = 3.5%; odds ratio = 1.38; P = 2.2 × 10-14). This gene encodes a transcription factor involved in gametogenesis that is normally only weakly expressed in plasma cells. The association is represented by 14 variants in linkage disequilibrium. Among these, rs75712673 maps to a genomic region with open chromatin in plasma cells, and upregulates SOHLH2 in this cell type. Moreover, rs75712673 influences transcriptional activity in luciferase assays, and shows a chromatin looping interaction with the SOHLH2 promoter. Our work provides novel insight into MM susceptibility.

1q21.1 distal copy number variants are associated with cerebral and cognitive alterations in humans

Low-frequency 1q21.1 distal deletion and duplication copy number variant (CNV) carriers are predisposed to multiple neurodevelopmental disorders, including schizophrenia, autism and intellectual disability. Human carriers display a high prevalence of micro- and macrocephaly in deletion and duplication carriers, respectively. The underlying brain structural diversity remains largely unknown. We systematically called CNVs in 38 cohorts from the large-scale ENIGMA-CNV collaboration and the UK Biobank and identified 28 1q21.1 distal deletion and 22 duplication carriers and 37,088 non-carriers (48% male) derived from 15 distinct magnetic resonance imaging scanner sites. With standardized methods, we compared subcortical and cortical brain measures (all) and cognitive performance (UK Biobank only) between carrier groups also testing for mediation of brain structure on cognition. We identified positive dosage effects of copy number on intracranial volume (ICV) and total cortical surface area, with the largest effects in frontal and cingulate cortices, and negative dosage effects on caudate and hippocampal volumes. The carriers displayed distinct cognitive deficit profiles in cognitive tasks from the UK Biobank with intermediate decreases in duplication carriers and somewhat larger in deletion carriers-the latter potentially mediated by ICV or cortical surface area. These results shed light on pathobiological mechanisms of neurodevelopmental disorders, by demonstrating gene dose effect on specific brain structures and effect on cognitive function.

Identification of Genetic Loci Shared Between Attention-Deficit/Hyperactivity Disorder, Intelligence, and Educational Attainment

BACKGROUND

Attention-deficit/hyperactivity disorder (ADHD) is a neurodevelopmental disorder that is consistently associated with lower levels of educational attainment. A recent large genome-wide association study identified common gene variants associated with ADHD, but most of the genetic architecture remains unknown.

METHODS

We analyzed independent genome-wide association study summary statistics for ADHD (19,099 cases and 34,194 controls), educational attainment (N = 842,499), and general intelligence (N = 269,867) using a conditional/conjunctional false discovery rate (FDR) statistical framework that increases power of discovery by conditioning the FDR on overlapping associations. The genetic variants identified were characterized in terms of function, expression, and biological processes.

RESULTS

We identified 58 linkage disequilibrium-independent ADHD-associated loci (conditional FDR < 0.01), of which 30 were shared between ADHD and educational attainment or general intelligence (conjunctional FDR < 0.01) and 46 were novel risk loci for ADHD.

CONCLUSIONS

These results expand on previous genetic and epidemiological studies and support the hypothesis of a shared genetic basis between these phenotypes. Although the clinical utility of the identified loci remains to be determined, they can be used as resources to guide future studies aiming to disentangle the complex etiologies of ADHD, educational attainment, and general intelligence.

Association of Copy Number Variation of the 15q11.2 BP1-BP2 Region With Cortical and Subcortical Morphology and Cognition

Importance

Recurrent microdeletions and duplications in the genomic region 15q11.2 between breakpoints 1 (BP1) and 2 (BP2) are associated with neurodevelopmental disorders. These structural variants are present in 0.5% to 1.0% of the population, making 15q11.2 BP1-BP2 the site of the most prevalent known pathogenic copy number variation (CNV). It is unknown to what extent this CNV influences brain structure and affects cognitive abilities.

Objective

To determine the association of the 15q11.2 BP1-BP2 deletion and duplication CNVs with cortical and subcortical brain morphology and cognitive task performance.

Design, Setting, and Participants

In this genetic association study, T1-weighted brain magnetic resonance imaging were combined with genetic data from the ENIGMA-CNV consortium and the UK Biobank, with a replication cohort from Iceland. In total, 203 deletion carriers, 45 247 noncarriers, and 306 duplication carriers were included. Data were collected from August 2015 to April 2019, and data were analyzed from September 2018 to September 2019.

Main Outcomes and Measures

The associations of the CNV with global and regional measures of surface area and cortical thickness as well as subcortical volumes were investigated, correcting for age, age2, sex, scanner, and intracranial volume. Additionally, measures of cognitive ability were analyzed in the full UK Biobank cohort.

Results

Of 45 756 included individuals, the mean (SD) age was 55.8 (18.3) years, and 23 754 (51.9%) were female. Compared with noncarriers, deletion carriers had a lower surface area (Cohen d = -0.41; SE, 0.08; P = 4.9 × 10-8), thicker cortex (Cohen d = 0.36; SE, 0.07; P = 1.3 × 10-7), and a smaller nucleus accumbens (Cohen d = -0.27; SE, 0.07; P = 7.3 × 10-5). There was also a significant negative dose response on cortical thickness (β = -0.24; SE, 0.05; P = 6.8 × 10-7). Regional cortical analyses showed a localization of the effects to the frontal, cingulate, and parietal lobes. Further, cognitive ability was lower for deletion carriers compared with noncarriers on 5 of 7 tasks.

Conclusions and Relevance

These findings, from the largest CNV neuroimaging study to date, provide evidence that 15q11.2 BP1-BP2 structural variation is associated with brain morphology and cognition, with deletion carriers being particularly affected. The pattern of results fits with known molecular functions of genes in the 15q11.2 BP1-BP2 region and suggests involvement of these genes in neuronal plasticity. These neurobiological effects likely contribute to the association of this CNV with neurodevelopmental disorders.

Dose response of the 16p11.2 distal copy number variant on intracranial volume and basal ganglia

Carriers of large recurrent copy number variants (CNVs) have a higher risk of developing neurodevelopmental disorders. The 16p11.2 distal CNV predisposes carriers to e.g., autism spectrum disorder and schizophrenia. We compared subcortical brain volumes of 12 16p11.2 distal deletion and 12 duplication carriers to 6882 non-carriers from the large-scale brain Magnetic Resonance Imaging collaboration, ENIGMA-CNV. After stringent CNV calling procedures, and standardized FreeSurfer image analysis, we found negative dose-response associations with copy number on intracranial volume and on regional caudate, pallidum and putamen volumes (β = -0.71 to -1.37; P < 0.0005). In an independent sample, consistent results were obtained, with significant effects in the pallidum (β = -0.95, P = 0.0042). The two data sets combined showed significant negative dose-response for the accumbens, caudate, pallidum, putamen and ICV (P = 0.0032, 8.9 × 10-6, 1.7 × 10-9, 3.5 × 10-12 and 1.0 × 10-4, respectively). Full scale IQ was lower in both deletion and duplication carriers compared to non-carriers. This is the first brain MRI study of the impact of the 16p11.2 distal CNV, and we demonstrate a specific effect on subcortical brain structures, suggesting a neuropathological pattern underlying the neurodevelopmental syndromes.

Brain age prediction using deep learning uncovers associated sequence variants

Machine learning algorithms can be trained to estimate age from brain structural MRI. The difference between an individual’s predicted and chronological age, predicted age difference (PAD), is a phenotype of relevance to aging and brain disease. Here, we present a new deep learning approach to predict brain age from a T1-weighted MRI. The method was trained on a dataset of healthy Icelanders and tested on two datasets, IXI and UK Biobank, utilizing transfer learning to improve accuracy on new sites. A genome-wide association study (GWAS) of PAD in the UK Biobank data (discovery set: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N=12378$$\end{document}N=12378, replication set: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N=4456$$\end{document}N=4456) yielded two sequence variants, rs1452628-T (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta =-0.08$$\end{document}β=−0.08, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P=1.15\times{10}^{-9}$$\end{document}P=1.15×10−9) and rs2435204-G (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta =0.102$$\end{document}β=0.102, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P=9.73\times 1{0}^{-12}$$\end{document}P=9.73×10−12). The former is near KCNK2 and correlates with reduced sulcal width, whereas the latter correlates with reduced white matter surface area and tags a well-known inversion at 17q21.31 (H2).

Machine learning algorithms can be trained to estimate age from brain structural MRI. Here, the authors introduce a new deep-learning-based age prediction approach, and then carry out a GWAS of the difference between predicted and chronological age, revealing two associated variants.

Attention-deficit hyperactivity disorder shares copy number variant risk with schizophrenia and autism spectrum disorder

Attention-deficit/hyperactivity disorder (ADHD) is a highly heritable common childhood-onset neurodevelopmental disorder. Some rare copy number variations (CNVs) affect multiple neurodevelopmental disorders such as intellectual disability, autism spectrum disorders (ASD), schizophrenia and ADHD. The aim of this study is to determine to what extent ADHD shares high risk CNV alleles with schizophrenia and ASD. We compiled 19 neuropsychiatric CNVs and test 14, with sufficient power, for association with ADHD in Icelandic and Norwegian samples. Eight associate with ADHD; deletions at 2p16.3 (NRXN1), 15q11.2, 15q13.3 (BP4 & BP4.5-BP5) and 22q11.21, and duplications at 1q21.1 distal, 16p11.2 proximal, 16p13.11 and 22q11.21. Six of the CNVs have not been associated with ADHD before. As a group, the 19 CNVs associate with ADHD (OR = 2.43, P = 1.6 × 10^-21), even when comorbid ASD and schizophrenia are excluded from the sample. These results highlight the pleiotropic effect of the neuropsychiatric CNVs and add evidence for ADHD, ASD and schizophrenia being related neurodevelopmental disorders rather than distinct entities.

Reciprocal White Matter Changes Associated With Copy Number Variation at 15q11.2 BP1-BP2: A Diffusion Tensor Imaging Study

BACKGROUND

The 15q11.2 BP1-BP2 cytogenetic region has been associated with learning and motor delays, autism, and schizophrenia. This region includes a gene that codes for the cytoplasmic FMR1 interacting protein 1 (CYFIP1). The CYFIP1 protein is involved in actin cytoskeletal dynamics and interacts with the fragile X mental retardation protein. Absence of fragile X mental retardation protein causes fragile X syndrome. Because abnormal white matter microstructure has been reported in both fragile X syndrome and psychiatric disorders, we looked at the impact of 15q11.2 BP1-BP2 dosage on white matter microstructure.

METHODS

Combining a brain-wide voxel-based approach and a regional-based analysis, we analyzed diffusion tensor imaging data from healthy individuals with the deletion (n = 30), healthy individuals with the reciprocal duplication (n = 27), and IQ-matched control subjects with no large copy number variants (n = 19), recruited from a large genotyped population sample.

RESULTS

We found global mirror effects (deletion > control > duplication) on fractional anisotropy. The deletion group showed widespread increased fractional anisotropy when compared with duplication. Regional analyses revealed a greater effect size in the posterior limb of the internal capsule and a tendency for decreased fractional anisotropy in duplication.

CONCLUSIONS

These results show a reciprocal effect of 15q11.2 BP1-BP2 on white matter microstructure, suggesting that reciprocal chromosomal imbalances may lead to opposite changes in brain structure. Findings in the deletion overlap with previous white matter differences reported in fragile X syndrome patients, suggesting common pathogenic mechanisms derived from disruptions of cytoplasmic CYFIP1-fragile X mental retardation protein complexes. Our data begin to identify specific components of the 15q11.2 BP1-BP2 phenotype and neurobiological mechanisms of potential relevance to the increased risk for disorder.

Identification of common genetic risk variants for autism spectrum disorder

Autism spectrum disorder (ASD) is a highly heritable and heterogeneous group of neurodevelopmental phenotypes diagnosed in more than 1% of children. Common genetic variants contribute substantially to ASD susceptibility, but to date no individual variants have been robustly associated with ASD. With a marked sample-size increase from a unique Danish population resource, we report a genome-wide association meta-analysis of 18,381 individuals with ASD and 27,969 controls that identified five genome-wide-significant loci. Leveraging GWAS results from three phenotypes with significantly overlapping genetic architectures (schizophrenia, major depression, and educational attainment), we identified seven additional loci shared with other traits at equally strict significance levels. Dissecting the polygenic architecture, we found both quantitative and qualitative polygenic heterogeneity across ASD subtypes. These results highlight biological insights, particularly relating to neuronal function and corticogenesis, and establish that GWAS performed at scale will be much more productive in the near term in ASD.

Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder

Attention deficit/hyperactivity disorder (ADHD) is a highly heritable childhood behavioral disorder affecting 5% of children and 2.5% of adults. Common genetic variants contribute substantially to ADHD susceptibility, but no variants have been robustly associated with ADHD. We report a genome-wide association meta-analysis of 20,183 individuals diagnosed with ADHD and 35,191 controls that identifies variants surpassing genome-wide significance in 12 independent loci, finding important new information about the underlying biology of ADHD. Associations are enriched in evolutionarily constrained genomic regions and loss-of-function intolerant genes and around brain-expressed regulatory marks. Analyses of three replication studies: a cohort of individuals diagnosed with ADHD, a self-reported ADHD sample and a meta-analysis of quantitative measures of ADHD symptoms in the population, support these findings while highlighting study-specific differences on genetic overlap with educational attainment. Strong concordance with GWAS of quantitative population measures of ADHD symptoms supports that clinical diagnosis of ADHD is an extreme expression of continuous heritable traits.

Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder

Attention deficit/hyperactivity disorder (ADHD) is a highly heritable childhood behavioral disorder affecting 5% of children and 2.5% of adults. Common genetic variants contribute substantially to ADHD susceptibility, but no variants have been robustly associated with ADHD. We report a genome-wide association meta-analysis of 20,183 individuals diagnosed with ADHD and 35,191 controls that identifies variants surpassing genome-wide significance in 12 independent loci, finding important new information about the underlying biology of ADHD. Associations are enriched in evolutionarily constrained genomic regions and loss-of-function intolerant genes and around brain-expressed regulatory marks. Analyses of three replication studies: a cohort of individuals diagnosed with ADHD, a self-reported ADHD sample and a meta-analysis of quantitative measures of ADHD symptoms in the population, support these findings while highlighting study-specific differences on genetic overlap with educational attainment. Strong concordance with GWAS of quantitative population measures of ADHD symptoms supports that clinical diagnosis of ADHD is an extreme expression of continuous heritable traits.

MAP1B mutations cause intellectual disability and extensive white matter deficit

Discovery of coding variants in genes that confer risk of neurodevelopmental disorders is an important step towards understanding the pathophysiology of these disorders. Whole-genome sequencing of 31,463 Icelanders uncovers a frameshift variant (E712KfsTer10) in microtubule-associated protein 1B (MAP1B) that associates with ID/low IQ in a large pedigree (genome-wide corrected P = 0.022). Additional stop-gain variants in MAP1B (E1032Ter and R1664Ter) validate the association with ID and IQ. Carriers have 24% less white matter (WM) volume (β = −2.1SD, P = 5.1 × 10⁻⁸), 47% less corpus callosum (CC) volume (β = −2.4SD, P = 5.5 × 10⁻¹⁰) and lower brain-wide fractional anisotropy (P = 6.7 × 10⁻⁴). In summary, we show that loss of MAP1B function affects general cognitive ability through a profound, brain-wide WM deficit with likely disordered or compromised axons.

Intellectual disability (ID) is characterized by an intelligence quotient of below 70 and impaired adaptive skills. Here, analyzing whole genome sequences from 31,463 Icelanders, Walters et al. identify variants in MAP1B associated with ID and extensive brain-wide white matter deficits.

Rare and Common Variants Conferring Risk of Tooth Agenesis

We present association results from a large genome-wide association study of tooth agenesis (TA) as well as selective TA, including 1,944 subjects with congenitally missing teeth, excluding third molars, and 338,554 controls, all of European ancestry. We also tested the association of previously identified risk variants, for timing of tooth eruption and orofacial clefts, with TA. We report associations between TA and 9 novel risk variants. Five of these variants associate with selective TA, including a variant conferring risk of orofacial clefts. These results contribute to a deeper understanding of the genetic architecture of tooth development and disease. The few variants previously associated with TA were uncovered through candidate gene studies guided by mouse knockouts. Knowing the etiology and clinical features of TA is important for planning oral rehabilitation that often involves an interdisciplinary approach.

Reproductive fitness and genetic risk of psychiatric disorders in the general population

The persistence of common, heritable psychiatric disorders that reduce reproductive fitness is an evolutionary paradox. Here, we investigate the selection pressures on sequence variants that predispose to schizophrenia, autism, bipolar disorder, major depression and attention deficit hyperactivity disorder (ADHD) using genomic data from 150,656 Icelanders, excluding those diagnosed with these psychiatric diseases. Polygenic risk of autism and ADHD is associated with number of children. Higher polygenic risk of autism is associated with fewer children and older age at first child whereas higher polygenic risk of ADHD is associated with having more children. We find no evidence for a selective advantage of a high polygenic risk of schizophrenia or bipolar disorder. Rare copy-number variants conferring moderate to high risk of psychiatric illness are associated with having fewer children and are under stronger negative selection pressure than common sequence variants.

Why genetic variants that confer risk for psychiatric disorders persist in the genome is an evolutionary conundrum. Here, Mullins et al. report association of polygenic risk for autism with having fewer children and polygenic risk for ADHD with higher reproductive fitness.

Sequence variant at 4q25 near PITX2 associates with appendicitis

Appendicitis is one of the most common conditions requiring acute surgery and can pose a threat to the lives of affected individuals. We performed a genome-wide association study of appendicitis in 7,276 Icelandic and 1,139 Dutch cases and large groups of controls. In a combined analysis of the Icelandic and Dutch data, we detected a single signal represented by an intergenic variant rs2129979 [G] close to the gene PITX2 associating with increased risk of appendicitis (OR = 1.15, P = 1.8 × 10-11). We only observe the association in patients diagnosed in adulthood. The marker is close to, but distinct from, a set of markers reported to associate with atrial fibrillation, which have been linked to PITX2. PITX2 has been implicated in determination of right-left symmetry during development. Anomalies in organ arrangement have been linked to increased prevalence of gastrointestinal and intra-abdominal complications, which may explain the effect of rs2129979 on appendicitis risk.

15q11.2 CNV affects cognitive, structural and functional correlates of dyslexia and dyscalculia

Several copy number variants have been associated with neuropsychiatric disorders and these variants have been shown to also influence cognitive abilities in carriers unaffected by psychiatric disorders. Previously, we associated the 15q11.2(BP1-BP2) deletion with specific learning disabilities and a larger corpus callosum. Here we investigate, in a much larger sample, the effect of the 15q11.2(BP1-BP2) deletion on cognitive, structural and functional correlates of dyslexia and dyscalculia. We report that the deletion confers greatest risk of the combined phenotype of dyslexia and dyscalculia. We also show that the deletion associates with a smaller left fusiform gyrus. Moreover, tailored functional magnetic resonance imaging experiments using phonological lexical decision and multiplication verification tasks demonstrate altered activation in the left fusiform and the left angular gyri in carriers. Thus, by using convergent evidence from neuropsychological testing, and structural and functional neuroimaging, we show that the 15q11.2(BP1-BP2) deletion affects cognitive, structural and functional correlates of both dyslexia and dyscalculia.

Epigenetic and genetic components of height regulation

Adult height is a highly heritable trait. Here we identified 31.6 million sequence variants by whole-genome sequencing of 8,453 Icelanders and tested them for association with adult height by imputing them into 88,835 Icelanders. Here we discovered 13 novel height associations by testing four different models including parent-of-origin (|β|=0.4-10.6 cm). The minor alleles of three parent-of-origin signals associate with less height only when inherited from the father and are located within imprinted regions (IGF2-H19 and DLK1-MEG3). We also examined the association of these sequence variants in a set of 12,645 Icelanders with birth length measurements. Two of the novel variants, (IGF2-H19 and TET1), show significant association with both adult height and birth length, indicating a role in early growth regulation. Among the parent-of-origin signals, we observed opposing parental effects raising questions about underlying mechanisms. These findings demonstrate that common variations affect human growth by parental imprinting.

CNVs conferring risk of autism or schizophrenia affect cognition in controls

In a small fraction of patients with schizophrenia or autism, alleles of copy-number variants (CNVs) in their genomes are probably the strongest factors contributing to the pathogenesis of the disease. These CNVs may provide an entry point for investigations into the mechanisms of brain function and dysfunction alike. They are not fully penetrant and offer an opportunity to study their effects separate from that of manifest disease. Here we show in an Icelandic sample that a few of the CNVs clearly alter fecundity (measured as the number of children by age 45). Furthermore, we use various tests of cognitive function to demonstrate that control subjects carrying the CNVs perform at a level that is between that of schizophrenia patients and population controls. The CNVs do not all affect the same cognitive domains, hence the cognitive deficits that drive or accompany the pathogenesis vary from one CNV to another. Controls carrying the chromosome 15q11.2 deletion between breakpoints 1 and 2 (15q11.2(BP1-BP2) deletion) have a history of dyslexia and dyscalculia, even after adjusting for IQ in the analysis, and the CNV only confers modest effects on other cognitive traits. The 15q11.2(BP1-BP2) deletion affects brain structure in a pattern consistent with both that observed during first-episode psychosis in schizophrenia and that of structural correlates in dyslexia.

A common biological basis of obesity and nicotine addiction

Smoking influences body weight such that smokers weigh less than non-smokers and smoking cessation often leads to weight increase. The relationship between body weight and smoking is partly explained by the effect of nicotine on appetite and metabolism. However, the brain reward system is involved in the control of the intake of both food and tobacco. We evaluated the effect of single-nucleotide polymorphisms (SNPs) affecting body mass index (BMI) on smoking behavior, and tested the 32 SNPs identified in a meta-analysis for association with two smoking phenotypes, smoking initiation (SI) and the number of cigarettes smoked per day (CPD) in an Icelandic sample (N=34,216 smokers). Combined according to their effect on BMI, the SNPs correlate with both SI (r=0.019, P=0.00054) and CPD (r=0.032, P=8.0 × 10(-7)). These findings replicate in a second large data set (N=127,274, thereof 76,242 smokers) for both SI (P=1.2 × 10(-5)) and CPD (P=9.3 × 10(-5)). Notably, the variant most strongly associated with BMI (rs1558902-A in FTO) did not associate with smoking behavior. The association with smoking behavior is not due to the effect of the SNPs on BMI. Our results strongly point to a common biological basis of the regulation of our appetite for tobacco and food, and thus the vulnerability to nicotine addiction and obesity.

Rate of de novo mutations and the importance of father's age to disease risk

Mutations generate sequence diversity and provide a substrate for selection. The rate of de novo mutations is therefore of major importance to evolution. Here we conduct a study of genome-wide mutation rates by sequencing the entire genomes of 78 Icelandic parent-offspring trios at high coverage. We show that in our samples, with an average father's age of 29.7, the average de novo mutation rate is 1.20 × 10(-8) per nucleotide per generation. Most notably, the diversity in mutation rate of single nucleotide polymorphisms is dominated by the age of the father at conception of the child. The effect is an increase of about two mutations per year. An exponential model estimates paternal mutations doubling every 16.5 years. After accounting for random Poisson variation, father's age is estimated to explain nearly all of the remaining variation in the de novo mutation counts. These observations shed light on the importance of the father's age on the risk of diseases such as schizophrenia and autism.

A genome-wide association search for type 2 diabetes genes in African Americans

African Americans are disproportionately affected by type 2 diabetes (T2DM) yet few studies have examined T2DM using genome-wide association approaches in this ethnicity. The aim of this study was to identify genes associated with T2DM in the African American population. We performed a Genome Wide Association Study (GWAS) using the Affymetrix 6.0 array in 965 African-American cases with T2DM and end-stage renal disease (T2DM-ESRD) and 1029 population-based controls. The most significant SNPs (n = 550 independent loci) were genotyped in a replication cohort and 122 SNPs (n = 98 independent loci) were further tested through genotyping three additional validation cohorts followed by meta-analysis in all five cohorts totaling 3,132 cases and 3,317 controls. Twelve SNPs had evidence of association in the GWAS (P<0.0071), were directionally consistent in the Replication cohort and were associated with T2DM in subjects without nephropathy (P<0.05). Meta-analysis in all cases and controls revealed a single SNP reaching genome-wide significance (P<2.5×10(-8)). SNP rs7560163 (P = 7.0×10(-9), OR (95% CI) = 0.75 (0.67-0.84)) is located intergenically between RND3 and RBM43. Four additional loci (rs7542900, rs4659485, rs2722769 and rs7107217) were associated with T2DM (P<0.05) and reached more nominal levels of significance (P<2.5×10(-5)) in the overall analysis and may represent novel loci that contribute to T2DM. We have identified novel T2DM-susceptibility variants in the African-American population. Notably, T2DM risk was associated with the major allele and implies an interesting genetic architecture in this population. These results suggest that multiple loci underlie T2DM susceptibility in the African-American population and that these loci are distinct from those identified in other ethnic populations.

Common genetic variants associated with open-angle glaucoma

Open-angle glaucoma (glaucoma) is a major eye disorder characterized by optic disc pathology. Recent genome-wide association studies identified new loci associated with clinically relevant optic disc parameters, such as the optic disc area and vertical cup-disc ratio (VCDR). We examined to what extent these loci are involved in glaucoma. The loci studied include ATOH7, CDC7/TGFBR3 and SALL1 for optic disc area, and CDKN2B, SIX1, SCYL1/LTBP3, CHEK2, ATOH7 and DCLK1 for VCDR. We performed a meta-analysis using data from six independent studies including: the Rotterdam Study (n= 5736), Genetic Research in Isolated Populations combined with Erasmus Rucphen Family study (n= 1750), Amsterdam Glaucoma Study (n= 296) and cohorts from Erlangen and Tübingen (n= 1363), Southampton (n= 702) and deCODE (n= 36 151) resulting in a total of 3161 glaucoma cases and 42 837 controls. Of the eight loci, we found significant evidence (P= 1.41 × 10(-8)) for the association of CDKN2B with glaucoma [odds ratio (OR) for those homozygous for the risk allele: 0.76; 95% confidence interval (CI): 0.70-0.84], for the role of ATOH7 (OR: 1.28; 95% CI: 1.12-1.47) and for SIX1 (OR: 1.20; 95% CI: 1.10-1.31) when adjusting for the number of tested loci. Furthermore, there was a borderline significant association of CDC7/TGFBR3 and SALL1 (both P= 0.04) with glaucoma. In conclusion, we found consistent evidence for three common variants (CDKN2B, ATOH7 and SIX1) significantly associated with glaucoma. These findings may shed new light on the pathophysiological protein pathways leading to glaucoma, and point to pathways involved in the growth and development of the optic nerve.

A rare variant in MYH6 is associated with high risk of sick sinus syndrome

Through complementary application of SNP genotyping, whole-genome sequencing and imputation in 38,384 Icelanders, we have discovered a previously unidentified sick sinus syndrome susceptibility gene, MYH6, encoding the alpha heavy chain subunit of cardiac myosin. A missense variant in this gene, c.2161C>T, results in the conceptual amino acid substitution p.Arg721Trp, has an allelic frequency of 0.38% in Icelanders and associates with sick sinus syndrome with an odds ratio = 12.53 and P = 1.5 × 10⁻²⁹. We show that the lifetime risk of being diagnosed with sick sinus syndrome is around 6% for non-carriers of c.2161C>T but is approximately 50% for carriers of the c.2161C>T variant.

Sequence variants at CYP1A1-CYP1A2 and AHR associate with coffee consumption

Coffee is the most commonly used stimulant and caffeine is its main psychoactive ingredient. The heritability of coffee consumption has been estimated at around 50%. We performed a meta-analysis of four genome-wide association studies of coffee consumption among coffee drinkers from Iceland (n = 2680), The Netherlands (n = 2791), the Sorbs Slavonic population isolate in Germany (n = 771) and the USA (n = 369) using both directly genotyped and imputed single nucleotide polymorphisms (SNPs) (2.5 million SNPs). SNPs at the two most significant loci were also genotyped in a sample set from Iceland (n = 2430) and a Danish sample set consisting of pregnant women (n = 1620). Combining all data, two sequence variants significantly associated with increased coffee consumption: rs2472297-T located between CYP1A1 and CYP1A2 at 15q24 (P = 5.4 · 10(-14)) and rs6968865-T near aryl hydrocarbon receptor (AHR) at 7p21 (P = 2.3 · 10(-11)). An effect of ∼0.2 cups a day per allele was observed for both SNPs. CYP1A2 is the main caffeine metabolizing enzyme and is also involved in drug metabolism. AHR detects xenobiotics, such as polycyclic aryl hydrocarbons found in roasted coffee, and induces transcription of CYP1A1 and CYP1A2. The association of these SNPs with coffee consumption was present in both smokers and non-smokers.

Common variants near CAV1 and CAV2 are associated with primary open-angle glaucoma

We conducted a genome-wide association study for primary open-angle glaucoma (POAG) in 1,263 affected individuals (cases) and 34,877 controls from Iceland. We identified a common sequence variant at 7q31 (rs4236601[A], odds ratio (OR) = 1.36, P = 5.0 × 10⁻¹⁰). We then replicated the association in sample sets of 2,175 POAG cases and 2,064 controls from Sweden, the UK and Australia (combined OR = 1.18, P = 0.0015) and in 299 POAG cases and 580 unaffected controls from Hong Kong and Shantou, China (combined OR = 5.42, P = 0.0021). The risk variant identified here is located close to CAV1 and CAV2, both of which are expressed in the trabecular meshwork and retinal ganglion cells that are involved in the pathogenesis of POAG.

Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis

By combining genome-wide association data from 8,130 individuals with type 2 diabetes (T2D) and 38,987 controls of European descent and following up previously unidentified meta-analysis signals in a further 34,412 cases and 59,925 controls, we identified 12 new T2D association signals with combined P<5x10(-8). These include a second independent signal at the KCNQ1 locus; the first report, to our knowledge, of an X-chromosomal association (near DUSP9); and a further instance of overlap between loci implicated in monogenic and multifactorial forms of diabetes (at HNF1A). The identified loci affect both beta-cell function and insulin action, and, overall, T2D association signals show evidence of enrichment for genes involved in cell cycle regulation. We also show that a high proportion of T2D susceptibility loci harbor independent association signals influencing apparently unrelated complex traits.

New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes.