Search

Somatic evolution of cells within the body is well known to lead to cancers. However, spread of somatic mutations within a tissue over time may also contribute to the pathogenesis of non-neoplastic diseases. Recent years have seen the publication of many studies aiming to characterize somatic evolution in healthy tissues. A logical next step is to extend such work to diseased conditions. As our understanding of the interplay between somatic mutations and non-neoplastic disease grows, opportunities for the joint study of germline and somatic variants will present themselves. Here, we present our thoughts on the utility of somatic mutations for understanding both the causes and consequences of common complex disease and the challenges that remain for the joint study of the soma and germline.

Funder: Wellcome Clinical PhD fellowship, Funder: Wellcome PhD Studentship, Funder: Jean Shank/Pathological Society Intermediate Fellowship, Mutation accumulation in somatic cells contributes to cancer development and is proposed as a cause of aging. DNA polymerases Pol ε and Pol δ replicate DNA during cell division. However, in some cancers, defective proofreading due to acquired POLE/POLD1 exonuclease domain mutations causes markedly elevated somatic mutation burdens with distinctive mutational signatures. Germline POLE/POLD1 mutations cause familial cancer predisposition. Here, we sequenced normal tissue and tumor DNA from individuals with germline POLE/POLD1 mutations. Increased mutation burdens with characteristic mutational signatures were found in normal adult somatic cell types, during early embryogenesis and in sperm. Thus human physiology can tolerate ubiquitously elevated mutation burdens. Except for increased cancer risk, individuals with germline POLE/POLD1 mutations do not exhibit overt features of premature aging. These results do not support a model in which all features of aging are attributable to widespread cell malfunction directly resulting from somatic mutation burdens accrued during life.

Somatic mutations are hypothesised to play a role in many non-neoplastic diseases. These diseases may also shape the somatic mutation landscape of affected tissues after onset. We performed whole-exome sequencing of 1182 microbiopsies dissected from lesional and non-lesional epidermis from 111 patients with psoriasis, a chronic inflammatory disease of the skin, to search for evidence that somatic mutations in keratinocytes may influence the disease process and to characterise the effects of the disease on the mutation landscape of the epidermis. We show that psoriasis is associated with increased mutation burden of the cell-intrinsic signatures SBS1 and SBS5 but not of UV-light, which remains the dominant mutagen in psoriatic skin. Despite the hyperproliferation of keratinocytes that characterises psoriasis, lesional skin remains highly polyclonal, showing no evidence of spread of clones carrying potentially pathogenic mutations. We find that the selection forces operating in the epidermis remain mostly unchanged in psoriasis and the mutational landscape continues to be dominated by clones carrying mutations in genes recurrently mutated in normal squamous epithelia. There is evidence of positive selection in previously reported driver genes, NOTCH1, NOTCH2, TP53, FAT1 and PPM1D and we also identify four driver genes (GXYLT1, CHEK2, ZFP36L2 and EEF1A1), that have not been previously described in studies of normal skin but which we hypothesise are selected for in squamous epithelium irrespective of disease status. We describe the mutagenic effects of psoralens, a class of chemicals previously found in some sunscreens and which remain a part of a common photochemotherapy treatment for psoriasis (psoralens and UV-A, PUVA). Psoralens leave a distinct mutational signature in the genomes of exposed cells that is tightly linked with transcription, showing evidence of both transcription-coupled repair and transcription-coupled damage. These results suggest that somatic mutations in keratinocytes are unlikely to influence the pathogenesis of psoriasis and that while psoriasis has only modest effect on the mutation landscape of the skin, PUVA treatment has the potential to exert a unique and larger effects.

Microsatellites are polymorphic tracts of short tandem repeats (STRs) with one to six base-pair (bp) motifs and are some of the most polymorphic markers in the genome. Using 6,084 Icelandic parent-offspring trios we estimate 63.7 (95% CI: 61.9–65.4) microsatellite de novo mutations (mDNMs) per offspring per generation. Paternal mDNMs occur at longer repeats, while maternal mDNMs affect more bp. mDNMs increase by 0.97 (95% CI: 0.90–1.04) and 0.31 (95% CI: 0.25–0.37) per year of father‘s and mother‘s age at conception, respectively. We found two independent coding variants associated with an increased number of mDNMs transmitted to offspring; A missense variant (allele frequency (AF) = 1.9%) in MSH2, a mismatch repair gene, increases transmitted mDNMs from both parents (effect: 13.1 paternal and 7.8 maternal mDNMs). A synonymous variant (AF = 20.3%) in NEIL2, a DNA damage repair gene, increases paternally transmitted mDNMs (effect: 4.4 mDNMs). Thus, the microsatellite mutation rate in humans is in part under genetic control.

Cellular DNA damage caused by reactive oxygen species is repaired by the base excision repair (BER) pathway which includes the DNA glycosylase MUTYH. Inherited biallelic MUTYH mutations cause predisposition to colorectal adenomas and carcinoma. However, the mechanistic progression from germline MUTYH mutations to MUTYH-Associated Polyposis (MAP) is incompletely understood. Here, we sequenced normal tissue DNAs from 10 individuals with MAP. Somatic base substitution mutation rates in intestinal epithelial cells were elevated 2 to 5-fold in all individuals, except for one showing a 33-fold increase, and were also increased in other tissues. The increased mutation burdens were of multiple mutational signatures characterised by C>A changes. Different mutation rates and signatures between individuals were likely due to different MUTYH mutations or additional inherited mutations in other BER pathway genes. The elevated base substitution rate in normal cells likely accounts for the predisposition to neoplasia in MAP. Despite ubiquitously elevated mutation rates, individuals with MAP do not display overt evidence of premature ageing. Thus, accumulation of somatic mutations may not be sufficient to cause the global organismal functional decline of ageing.Summary

Somatic mutations drive the development of cancer and may contribute to ageing and other diseases1,2. Despite their importance, the difficulty of detecting mutations that are only present in single cells or small clones has limited our knowledge of somatic mutagenesis to a minority of tissues. Here, to overcome these limitations, we developed nanorate sequencing (NanoSeq), a duplex sequencing protocol with error rates of less than five errors per billion base pairs in single DNA molecules from cell populations. This rate is two orders of magnitude lower than typical somatic mutation loads, enabling the study of somatic mutations in any tissue independently of clonality. We used this single-molecule sensitivity to study somatic mutations in non-dividing cells across several tissues, comparing stem cells to differentiated cells and studying mutagenesis in the absence of cell division. Differentiated cells in blood and colon displayed remarkably similar mutation loads and signatures to their corresponding stem cells, despite mature blood cells having undergone considerably more divisions. We then characterized the mutational landscape of post-mitotic neurons and polyclonal smooth muscle, confirming that neurons accumulate somatic mutations at a constant rate throughout life without cell division, with similar rates to mitotically active tissues. Together, our results suggest that mutational processes that are independent of cell division are important contributors to somatic mutagenesis. We anticipate that the ability to reliably detect mutations in single DNA molecules could transform our understanding of somatic mutagenesis and enable non-invasive studies on large-scale cohorts. NanoSeq is used to detect mutations in single DNA molecules and analyses show that mutational processes that are independent of cell division are important contributors to somatic mutagenesis.

Mutation accumulation over time in normal somatic cells contributes to cancer development and is proposed as a cause of ageing. DNA polymerases Pol ε and Pol δ replicate DNA with high fidelity during normal cell divisions. However, in some cancers defective proofreading due to acquired mutations in the exonuclease domains of POLE or POLD1 causes markedly elevated somatic mutation burdens with distinctive mutational signatures. POLE and POLD1 exonuclease domain mutations also cause familial cancer predisposition when inherited through the germline. Here, we sequenced normal tissue DNA from individuals with germline POLE or POLD1 exonuclease domain mutations. Increased mutation burdens with characteristic mutational signatures were found to varying extents in all normal adult somatic cell types examined, during early embryogenesis and in sperm. Mutation burdens were further markedly elevated in neoplasms from these individuals. Thus human physiology is able to tolerate ubiquitously elevated mutation burdens. Indeed, with the exception of early onset cancer, individuals with germline POLE and POLD1 exonuclease domain mutations are not reported to show abnormal phenotypic features, including those of premature ageing. The results, therefore, do not support a simple model in which all features of ageing are attributable to widespread cell malfunction directly resulting from somatic mutation burdens accrued during life.

Item does not contain fulltext 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.

Ingileif Jonsdottir, Bjorn Nilsson, Kari Stefansson and colleagues perform a genome-wide association study for immunoglobulin levels in Icelandic and Swedish cohorts. They find 38 new variants associated with IgA, IgG, IgM or composite immunoglobulin traits and identify candidate genes underlying the regulation of immunoglobulin levels. Immunoglobulins are the effector molecules of the adaptive humoral immune system. In a genome-wide association study of 19,219 individuals, we found 38 new variants and replicated 5 known variants associating with IgA, IgG or IgM levels or with composite immunoglobulin traits, accounted for by 32 loci. Variants at these loci also affect the risk of autoimmune diseases and blood malignancies and influence blood cell development. Notable associations include a rare variant at RUNX3 decreasing IgA levels by shifting isoform proportions (rs188468174[C>T]: P = 8.3 × 10−55, β = −0.90 s.d.), a rare in-frame deletion in FCGR2B abolishing IgG binding to the encoded receptor (p.Asn106del: P = 4.2 × 10−8, β = 1.03 s.d.), four IGH locus variants influencing class switching, and ten new associations with the HLA region. Our results provide new insight into the regulation of humoral immunity.

We use polygenic risk scores (PRSs) for schizophrenia (SCZ) and bipolar disorder (BPD) to predict smoking, and addiction to nicotine, alcohol or drugs in individuals not diagnosed with psychotic disorders. Using PRSs for 144 609 subjects, including 10 036 individuals admitted for in-patient addiction treatment and 35 754 smokers, we find that diagnoses of various substance use disorders and smoking associate strongly with PRSs for SCZ (P = 5.3 × 10-50 -1.4 × 10-6 ) and BPD (P = 1.7 × 10-9 -1.9 × 10-3 ), showing shared genetic etiology between psychosis and addiction. Using standardized scores for SCZ and BPD scaled to a unit increase doubling the risk of the corresponding disorder, the odds ratios for alcohol and substance use disorders range from 1.19 to 1.31 for the SCZ-PRS, and from 1.07 to 1.29 for the BPD-PRS. Furthermore, we show that as regular smoking becomes more stigmatized and less prevalent, these biological risk factors gain importance as determinants of the behavior.

Background Under normal physiological conditions, colonic crypts accrue ~40 somatic mutations for every year of life. That somatic mutations contribute to the development of cancer is well established, but their patterns, burden and functional consequences in diseases other than cancer have not been extensively studied and our understanding of the effects of chronic inflammation on the mutational profile and clonal structure of the colon is limited. Here, we investigated how the recurrent cycles of inflammation, ulceration and regeneration seen in IBD impact the mutational and clonal structure of intestinal epithelia. Methods We isolated and whole-genome sequenced ~400 individual colonic crypts from 46 IBD patients and compared these to ~400 crypts from 41 non-IBD controls. We compared the mutation burden, mutational signature exposure, clonal structure and cancer driver mutation landscape in crypts from actively and previously inflamed regions with crypts dissected from controls. Results We estimated the base substitution rate of affected colonic epithelial cells to be doubled after IBD onset. This change was primarily driven by acceleration of mutational processes ubiquitously observed in normal colon (Figure 1), and we did not detect an IBD-specific mutational process. In contrast to the normal colon, where clonal expansions outside the confines of the crypt are rare, we observed widespread millimetre-scale clonal expansions, even in the absence of mutations in KRAS, TP53 and APC (Figure 2). We discovered that mutations in ARID1A, PIGR and ZC3H12A, and genes in the interleukin 17 and Toll-like receptor pathways, were under positive selection in colonic crypts from IBD patients (Figure 3). With the exception of ARID1A, these genes and pathways have not been previously associated with cancer risk. A previously published mouse model of ZC3H12A suggests that LoF mutations in this gene may facilitate healing of affected mucosa without promoting tumorigenesis. This could make the encoded protein an attractive drug target. The observed enrichment of mutations in PIGR and IL17 and TLR pathways suggests that somatic mutations may initiate, maintain or perpetuate IBD pathogenesis through disruption of microbe-epithelial homeostasis. Conclusion Our results provide new insights into the consequences of chronic intestinal inflammation on the mutational profile and clonal structure of colonic epithelia. We identify the mutagens driving the increase in mutation burden and mutations which are under positive selection in the context of inflammation. Our results suggest that studying somatic mutations in the colon can reveal putative drug targets and pathogenic mechanisms for IBD.

Summary Inflammatory bowel disease (IBD) is a chronic inflammatory disease associated with increased risk of gastrointestinal cancers. We whole-genome sequenced 446 colonic crypts from 46 IBD patients and compared these to 412 crypts from 41 non-IBD controls from our previous publication on the mutation landscape of the normal colon. The average mutation rate of affected colonic epithelial cells is 2.4-fold that of healthy colon, and this increase is mostly driven by acceleration of mutational processes ubiquitously observed in normal colon. In contrast to the normal colon, where clonal expansions outside the confines of the crypt are rare, we observed widespread millimeter-scale clonal expansions. We discovered non-synonymous mutations in ARID1A, FBXW7, PIGR, ZC3H12A, and genes in the interleukin 17 and Toll-like receptor pathways, under positive selection in IBD. These results suggest distinct selection mechanisms in the colitis-affected colon and that somatic mutations potentially play a causal role in IBD pathogenesis., Graphical Abstract, Highlights • IBD-affected colons accrue substitutions and indels 2.4 and 7 times faster than normal • 17 signatures of mutational processes, including treatment • Millimeter-scale clonal expansions late in molecular time • Distinct mechanisms of positive selection of mutations in immune-related genes, Whole-genome sequencing of inflammatory bowel disease patient samples allows insight into mutational burdens and processes associated with disease, including putative driver mutations positively selected in the diseased colon.

The colorectal adenoma-carcinoma sequence has provided a paradigmatic framework for understanding the successive somatic genetic changes and consequent clonal expansions that lead to cancer1. However, our understanding of the earliest phases of colorectal neoplastic changes-which may occur in morphologically normal tissue-is comparatively limited, as for most cancer types. Here we use whole-genome sequencing to analyse hundreds of normal crypts from 42 individuals. Signatures of multiple mutational processes were revealed; some of these were ubiquitous and continuous, whereas others were only found in some individuals, in some crypts or during certain periods of life. Probable driver mutations were present in around 1% of normal colorectal crypts in middle-aged individuals, indicating that adenomas and carcinomas are rare outcomes of a pervasive process of neoplastic change across morphologically normal colorectal epithelium. Colorectal cancers exhibit substantially increased mutational burdens relative to normal cells. Sequencing normal colorectal cells provides quantitative insights into the genomic and clonal evolution of cancer.

These authors contributed equally: Lara Kular, Yun Liu, Ingrid Kockum, Andrew P. Feinberg, Maja Jagodic., The human leukocyte antigen (HLA) haplotype DRB1*15:01 is the major risk factor for multiple sclerosis (MS). Here, we find that DRB1*15:01 is hypomethylated and predominantly expressed in monocytes among carriers of DRB1*15:01. A differentially methylated region (DMR) encompassing HLA-DRB1 exon 2 is particularly affected and displays methylation-sensitive regulatory properties in vitro. Causal inference and Mendelian randomization provide evidence that HLA variants mediate risk for MS via changes in the HLA-DRB1 DMR that modify HLA-DRB1 expression. Meta-analysis of 14,259 cases and 171,347 controls confirms that these variants confer risk from DRB1*15:01 and also identifies a protective variant (rs9267649, p, This work was supported by grants from the Swedish Research Council, the Swedish Association for Persons with Neurological Disabilities, the Swedish Brain Foundation, the Swedish MS Foundation, Petrus and Augusta Hedlunds Foundation, the Swedish AFA Insurance, Knut and Alice Wallenberg Foundation, the Stockholm County Council (ALF project), and AstraZeneca (AstraZeneca-Science for Life Laboratory collaboration). A.P.F. received grant support from NIH (DP1 ES022579). L.K. was supported by fellowship from the Margaretha af Ugglas Foundation. D.G.-C. and J.T. were supported by EU FP7 306000 STATegra. Y.L. was supported in part by the National Natural Science Foundation (grant no. 31471212). We acknowledge the International Multiple Sclerosis Genetics Consortium (IMSGC) for providing SNP genotypes used in this study and BEA core facility (Karolinska Institutet) for processing 450K array data on monocytes. The computations were performed on resources provided by SNIC through Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX). The Norwegian MS Registry and Biobank and the MS research group in Oslo are acknowledged for access to data from Norwegian MS patients. B.H. was supported by the German research foundation in the framework of the Collaborative research group TR128, the German MS competence network, and the EU project MultipleMS. This work was supported by the German Ministry for Education and Research (BMBF) as part of the “German Competence Network Multiple Sclerosis” (KKNMS) (grant nos. 01GI0916 and 01GI0917). F.Z., H.W., and R.G. were supported by the German research foundation in the framework of the Collaborative research group TR128, the German MS competence network. The KORA study was initiated and financed by the Helmholtz Zentrum München-German Research Center for Environmental Health, which is funded by the BMBF and by the State of Bavaria. Furthermore, KORA research was supported within the Munich Center of Health Sciences (MC-Health), Ludwig-Maximilians-Universität, as part of LMUinnovativ. The collection of probands in the Heinz Nixdorf RECALL Study (HNR) (PIs: K.-H. Jöckel, R. Erbel) was supported by the Heinz Nixdorf Foundation. The genotyping of HNR probands was financed through a grant of the BMBF to M. M. Nöthen. The Dortmund Health Study was supported by the German Migraine and Headache Society (DMKG) and unrestricted grants of equal share from Almirall, AstraZeneca, Berlin-Chemie, Boehringer, Boots Healthcare, GlaxoSmithKline (GSK), Janssen-Cilag, McNeil Pharma, Merck Sharp & Dohme (MSD), and Pfizer to the University of Münster. Blood collection was done through funds from the Institute of Epidemiology and Social Medicine, University of Münster (K. Berger and J. Wellmann), genotyping was supported by the BMBF (grant no. 01ER0816). SHIP is part of the Community Medicine Research Network of the University Medicine Greifswald, Germany (www.community-medicine.de), which was initiated and funded by the BMBF (grants no. 01ZZ9603, 01ZZ0103, and 01ZZ0403), the Ministry of Cultural Affairs and the Social Ministry of the Federal State of Mecklenburg-West Pomerania; genome-wide data have been supported by the BMBF (grant no. 03ZIK012). The FoCus study was supported by the BMBF (grant no. 0315540A).

A meta-analysis of publicly available summary statistics on multiple sclerosis combined with three Nordic multiple sclerosis cohorts (21,079 cases, 371,198 controls) revealed seven sequence variants associating with multiple sclerosis, not reported previously. Using polygenic risk scores based on public summary statistics of variants outside the major histocompatibility complex region we quantified genetic overlap between common autoimmune diseases in Icelanders and identified disease clusters characterized by autoantibody presence/absence. As multiple sclerosis-polygenic risk scores captures the risk of primary biliary cirrhosis and vice versa (P = 1.6 × 10−7, 4.3 × 10−9) we used primary biliary cirrhosis as a proxy-phenotype for multiple sclerosis, the idea being that variants conferring risk of primary biliary cirrhosis have a prior probability of conferring risk of multiple sclerosis. We tested 255 variants forming the primary biliary cirrhosis-polygenic risk score and found seven multiple sclerosis-associating variants not correlated with any previously established multiple sclerosis variants. Most of the variants discovered are close to or within immune-related genes. One is a low-frequency missense variant in TYK2, another is a missense variant in MTHFR that reduces the function of the encoded enzyme affecting methionine metabolism, reported to be dysregulated in multiple sclerosis brain., We are grateful to all the participants in the study. We thank Theodora Baldursdottir and Ottar M. Bergmann for validating the clinical information on PBC in Iceland. We acknowledge that this work would not have been possible without the important contribution of the staff at recruitment centers and genotyping and informatics facilities in deCODE Genetics and the MS research group at Oslo University Hospital and The Norwegian MS Registry and Biobank, Bergen Norway. We also thank the International MS Genetics Consortium and the Wellcome Trust Case Control Consortium for collaboration in the genotyping of the Norwegian samples. The MS research at Karolinska Institutet, as a basis for the Swedish MS case-control material has received grant support from the Swedish Research Council, The Knut and Alice Wallenberg Foundation, The AFA Foundation and the Swedish Brain Foundation. We are grateful for the access to the GWAS and summary statistics for the autoimmune diseases in Immunobase and IBD Genetics.

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., The authors are grateful to the participants and we thank the research nurses and staff at the Krokhals recruitment centre. We also thank the staff at deCODE genetics core facilities for their important contribution to this work. The research leading to these results has received support from the European Community's Seventh Framework Programme under the Marie Curie Industry-Academia Partnership and Pathways (PsychDPC, GA 286213), the Innovative Medicines Initiative Joint Undertaking under grant agreement 115008 and in part by the National Institute for Health Research (NIHR) Biomedical Research Centre for Mental Health at South London and Maudsley NHS Foundation Trust and King's College London.

Diverticular disease is characterized by pouches (that is, diverticulae) due to weakness in the bowel wall, which can become infected and inflamed causing diverticulitis, with potentially severe complications. Here, we test 32.4 million sequence variants identified through whole-genome sequencing (WGS) of 15,220 Icelanders for association with diverticular disease (5,426 cases) and its more severe form diverticulitis (2,764 cases). Subsequently, 16 sequence variants are followed up in a diverticular disease sample from Denmark (5,970 cases, 3,020 controls). In the combined Icelandic and Danish data sets we observe significant association of intronic variants in ARHGAP15 (Rho GTPase-activating protein 15; rs4662344-T: P=1.9 × 10−18, odds ratio (OR)=1.23) and COLQ (collagen-like tail subunit of asymmetric acetylcholinesterase; rs7609897-T: P=1.5 × 10−10, OR=0.87) with diverticular disease and in FAM155A (family with sequence similarity 155A; rs67153654-A: P=3.0 × 10−11, OR=0.82) with diverticulitis. These are the first loci shown to associate with diverticular disease in a genome-wide study., A hallmark of diverticular disease is pouches in the bowel wall which can become infected and inflamed, causing the more severe diverticulitis. Here, the authors report the first genome-wide association study on these interconnected conditions and identify ARHGAP15, COLQ and FAM155A as novel risk loci.

The polygenic risk score for schizophrenia (SZ-PRS) is associated with creativity and psychiatric disorders, including cannabis use. We have investigated the selection pressures on sequence variants that predispose to schizophrenia in the Icelandic population and find no evidence of selective advantage of a high SZ-PRS. Rare copy-number variants conferring moderate to high risk of psychiatric illness are, however, associated with having fewer children and are under negative selection. We also examined the correlation between SZ-PRS and the number of relationships producing offspring in 81,713 adult Icelanders without diagnosis of schizophrenia. SZ-PRS is associated with increased likelihood of having offspring with more than one mate (P = 7.9 x 10-12). Thus, while common variants that increase schizophrenia risk may not be under strong selection they influence mating behaviour in adults who do not have the disorder.

Only a few common variants in the sequence of the genome have been shown to impact cognitive traits. Here we demonstrate that polygenic scores of educational attainment predict specific aspects of childhood cognition, as measured with IQ. Recently, three sequence variants were shown to associate with educational attainment, a confluence phenotype of genetic and environmental factors contributing to academic success. We show that one of these variants associating with educational attainment, rs4851266-T, also associates with Verbal IQ in dyslexic children (P = 4.3 × 10−4, β = 0.16 s.d.). The effect of 0.16 s.d. corresponds to 1.4 IQ points for heterozygotes and 2.8 IQ points for homozygotes. We verified this association in independent samples consisting of adults (P = 8.3 × 10−5, β = 0.12 s.d., combined P = 2.2 x 10−7, β = 0.14 s.d.). Childhood cognition is unlikely to be affected by education attained later in life, and the variant explains a greater fraction of the variance in verbal IQ than in educational attainment (0.7% vs 0.12%,. P = 1.0 × 10−5).

Multiple sclerosis (MS) is a chronic inflammatory, likely autoimmune disease of the central nervous system with a combination of genetic and environmental risk factors, among which Epstein-Barr virus (EBV) infection is a strong suspect. We have previously identified increased autoantibody levels toward the chloride-channel protein Anoctamin 2 (ANO2) in MS. Here, IgG antibody reactivity toward ANO2 and EBV nuclear antigen 1 (EBNA1) was measured using bead-based multiplex serology in plasma samples from 8,746 MS cases and 7,228 controls. We detected increased anti-ANO2 antibody levels in MS (P = 3.5 × 10(−36)) with 14.6% of cases and 7.8% of controls being ANO2 seropositive (odds ratio [OR] = 1.6; 95% confidence intervals [95%CI]: 1.5 to 1.8). The MS risk increase in ANO2-seropositive individuals was dramatic when also exposed to 3 known risk factors for MS: HLA-DRB1*15:01 carriage, absence of HLA-A*02:01, and high anti-EBNA1 antibody levels (OR = 24.9; 95%CI: 17.9 to 34.8). Reciprocal blocking experiments with ANO2 and EBNA1 peptides demonstrated antibody cross-reactivity, mapping to ANO2 [aa 140 to 149] and EBNA1 [aa 431 to 440]. HLA gene region was associated with anti-ANO2 antibody levels and HLA-DRB1*04:01 haplotype was negatively associated with ANO2 seropositivity (OR = 0.6; 95%CI: 0.5 to 0.7). Anti-ANO2 antibody levels were not increased in patients from 3 other inflammatory disease cohorts. The HLA influence and the fact that specific IgG production usually needs T cell help provides indirect evidence for a T cell ANO2 autoreactivity in MS. We propose a hypothesis where immune reactivity toward EBNA1 through molecular mimicry with ANO2 contributes to the etiopathogenesis of MS.