Evidence for the additivity of rare and common variant burden throughout the spectrum of intellectual disability.

Intellectual disability (ID) is a common disorder, yet there is a wide spectrum of impairment from mild to profoundly affected individuals. Mild ID is seen as the low extreme of the general distribution of intelligence, while severe ID is often seen as a monogenic disorder caused by rare, pathogenic, highly penetrant variants. To investigate the genetic factors influencing mild and severe ID, we evaluated rare and common variation in the Northern Finland Intellectual Disability cohort (n = 1096 ID patients), a cohort with a high percentage of mild ID (n = 550) and from a population bottleneck enriched in rare, damaging variation. Despite this enrichment, we found only a small percentage of ID was due to recessive Finnish-enriched variants (0.5%). A larger proportion was linked to dominant variation, with a significant burden of rare, damaging variation in both mild and severe ID. This rare variant burden was enriched in more severe ID (p = 2.4e-4), patients without a relative with ID (p = 4.76e-4), and in those with features associated with monogenic disorders. We also found a significant burden of common variants associated with decreased cognitive function, with no difference between mild and more severe ID. When we included common and rare variants in a joint model, the rare and common variants had additive effects in both mild and severe ID. A multimodel inference approach also found that common and rare variants together best explained ID status (ΔAIC = 16.8, ΔBIC = 10.2). Overall, we report evidence for the additivity of rare and common variant burden throughout the spectrum of intellectual disability.

Genome-wide screen of otosclerosis in population biobanks: 27 loci and shared associations with skeletal structure.

Otosclerosis is one of the most common causes of conductive hearing loss, affecting 0.3% of the population. It typically presents in adulthood and half of the patients have a positive family history. The pathophysiology of otosclerosis is poorly understood. A previous genome-wide association study (GWAS) identified a single association locus in an intronic region of RELN. Here, we report a meta-analysis of GWAS studies of otosclerosis in three population-based biobanks comprising 3504 cases and 861,198 controls. We identify 23 novel risk loci (p < 5 × 10-8) and report an association in RELN and three previously reported candidate gene or linkage regions (TGFB1, MEPE, and OTSC7). We demonstrate developmental stage-dependent immunostaining patterns of MEPE and RUNX2 in mouse otic capsules. In most association loci, the nearest protein-coding genes are implicated in bone remodelling, mineralization or severe skeletal disorders. We highlight multiple genes involved in transforming growth factor beta signalling for follow-up studies.

Large-scale sequencing identifies multiple genes and rare variants associated with Crohn's disease susceptibility.

Genome-wide association studies (GWASs) have identified hundreds of loci associated with Crohn's disease (CD). However, as with all complex diseases, robust identification of the genes dysregulated by noncoding variants typically driving GWAS discoveries has been challenging. Here, to complement GWASs and better define actionable biological targets, we analyzed sequence data from more than 30,000 patients with CD and 80,000 population controls. We directly implicate ten genes in general onset CD for the first time to our knowledge via association to coding variation, four of which lie within established CD GWAS loci. In nine instances, a single coding variant is significantly associated, and in the tenth, ATG4C, we see additionally a significantly increased burden of very rare coding variants in CD cases. In addition to reiterating the central role of innate and adaptive immune cells as well as autophagy in CD pathogenesis, these newly associated genes highlight the emerging role of mesenchymal cells in the development and maintenance of intestinal inflammation.

The 22q11.2 region regulates presynaptic gene-products linked to schizophrenia.

It is unclear how the 22q11.2 deletion predisposes to psychiatric disease. To study this, we generated induced pluripotent stem cells from deletion carriers and controls and utilized CRISPR/Cas9 to introduce the heterozygous deletion into a control cell line. Here, we show that upon differentiation into neural progenitor cells, the deletion acted in trans to alter the abundance of transcripts associated with risk for neurodevelopmental disorders including autism. In excitatory neurons, altered transcripts encoded presynaptic factors and were associated with genetic risk for schizophrenia, including common and rare variants. To understand how the deletion contributed to these changes, we defined the minimal protein-protein interaction network that best explains gene expression alterations. We found that many genes in 22q11.2 interact in presynaptic, proteasome, and JUN/FOS transcriptional pathways. Our findings suggest that the 22q11.2 deletion impacts genes that may converge with psychiatric risk loci to influence disease manifestation in each deletion carrier.

Rare coding variants in ten genes confer substantial risk for schizophrenia.

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

NCOR2 is a novel candidate gene for migraine-epilepsy phenotype.

HYPOTHESIS: To identify genetic factors predisposing to migraine-epilepsy phenotype utilizing a multi-generational family with known linkage to chr12q24.2-q24.3. METHODS: We used single nucleotide polymorphism (SNP) genotyping and next-generation sequencing technologies to perform linkage, haplotype, and variant analyses in an extended Finnish migraine-epilepsy family (n = 120). In addition, we used a large genome-wide association study (GWAS) dataset of migraine and two biobank studies, UK Biobank and FinnGen, to test whether variants within the susceptibility region associate with migraine or epilepsy related phenotypes in a population setting. RESULTS: The family showed the highest evidence of linkage (LOD 3.42) between rs7966411 and epilepsy. The haplotype shared among 12 out of 13 epilepsy patients in the family covers almost the entire NCOR2 and co-localizes with one of the risk loci of the recent GWAS on migraine. The haplotype harbors nine low-frequency variants with potential regulatory functions. Three of them, in addition to two common variants, show nominal associations with neurological disorders in either UK Biobank or FinnGen. CONCLUSION: We provide several independent lines of evidence supporting association between migraine-epilepsy phenotype and NCOR2. Our study suggests that NCOR2 may have a role in both migraine and epilepsy and thus would provide evidence for shared pathophysiology underlying these two diseases.

A novel variant in SMG9 causes intellectual disability, confirming a role for nonsense-mediated decay components in neurocognitive development.

Biallelic loss-of-function variants in the SMG9 gene, encoding a regulatory subunit of the mRNA nonsense-mediated decay (NMD) machinery, are reported to cause heart and brain malformation syndrome. Here we report five patients from three unrelated families with intellectual disability (ID) and a novel pathogenic SMG9 c.551 T > C p.(Val184Ala) homozygous missense variant, identified using exome sequencing. Sanger sequencing confirmed recessive segregation in each family. SMG9 c.551T > C p.(Val184Ala) is most likely an autozygous variant identical by descent. Characteristic clinical findings in patients were mild to moderate ID, intention tremor, pyramidal signs, dyspraxia, and ocular manifestations. We used RNA sequencing of patients and age- and sex-matched healthy controls to assess the effect of the variant. RNA sequencing revealed that the SMG9 c.551T > C variant did not affect the splicing or expression level of SMG9 gene products, and allele-specific expression analysis did not provide evidence that the nonsense mRNA-induced NMD was affected. Differential gene expression analysis identified prevalent upregulation of genes in patients, including the genes SMOX, OSBP2, GPX3, and ZNF155. These findings suggest that normal SMG9 function may be involved in transcriptional regulation without affecting nonsense mRNA-induced NMD. In conclusion, we demonstrate that the SMG9 c.551T > C missense variant causes a neurodevelopmental disorder and impacts gene expression. NMD components have roles beyond aberrant mRNA degradation that are crucial for neurocognitive development.

Polygenic risk provides biological validity for the ICHD-3 criteria among Finnish migraine families.

BACKGROUND: Migraine is diagnosed using the extensively field-tested International Classification of Headache Disorders (ICHD-3) consensus criteria derived by the International Headache Society. To evaluate the criteria in respect to a measurable biomarker, we studied the relationship between the main ICHD-3 criteria and the polygenic risk score, a measure of common variant burden in migraine. METHODS: We used linear mixed models to study the correlation of ICHD-3 diagnostic criteria, underlying symptoms, and main diagnoses with the polygenic risk score of migraine in a cohort of 8602 individuals from the Finnish Migraine Genome Project. RESULTS: Main diagnostic categories and all underlying diagnostic criteria formed a consistent continuum along the increasing polygenic burden. Polygenic risk was associated with the heterogeneous clinical picture starting from the non-migraine headache (mean 0.07; 95% CI 0.02-0.12; p = 0.008 compared to the non-headache group), to probable migraine (mean 0.13; 95% CI 0.08-0.18; p < 0.001), migraine headache (mean 0.17; 95% CI 0.14-0.21; p < 0.001) and migraine with typical visual aura (mean 0.29; 95% CI 0.26-0.33; p < 0.001), all the way to the hemiplegic aura (mean 0.37; 95% CI 0.31-0.43; p < 0.001). All individual ICHD-3 symptoms and the total number of reported symptoms, a surrogate of migraine complexity, demonstrated a clear inclination with an increasing polygenic risk. CONCLUSIONS: The complex migraine phenotype progressively follows the polygenic burden from individuals with no headache to non-migrainous headache and up to patients with attacks manifesting all the features of the ICHD-3 headache and aura. Results provide further biological support for the ICHD-3 diagnostic criteria.

Contribution of rare and common variants to intellectual disability in a sub-isolate of Northern Finland.

The contribution of de novo variants in severe intellectual disability (ID) has been extensively studied whereas the genetics of mild ID has been less characterized. To elucidate the genetics of milder ID we studied 442 ID patients enriched for mild ID (>50%) from a population isolate of Finland. Using exome sequencing, we show that rare damaging variants in known ID genes are observed significantly more often in severe (27%) than in mild ID (13%) patients. We further observe a significant enrichment of functional variants in genes not yet associated with ID (OR: 2.1). We show that a common variant polygenic risk significantly contributes to ID. The heritability explained by polygenic risk score is the highest for educational attainment (EDU) in mild ID (2.2%) but lower for more severe ID (0.6%). Finally, we identify a Finland enriched homozygote variant in the CRADD ID associated gene.

Quantifying the Impact of Rare and Ultra-rare Coding Variation across the Phenotypic Spectrum.

There is a limited understanding about the impact of rare protein-truncating variants across multiple phenotypes. We explore the impact of this class of variants on 13 quantitative traits and 10 diseases using whole-exome sequencing data from 100,296 individuals. Protein-truncating variants in genes intolerant to this class of mutations increased risk of autism, schizophrenia, bipolar disorder, intellectual disability, and ADHD. In individuals without these disorders, there was an association with shorter height, lower education, increased hospitalization, and reduced age at enrollment. Gene sets implicated from GWASs did not show a significant protein-truncating variants burden beyond what was captured by established Mendelian genes. In conclusion, we provide a thorough investigation of the impact of rare deleterious coding variants on complex traits, suggesting widespread pleiotropic risk.

Common Variant Burden Contributes to the Familial Aggregation of Migraine in 1,589 Families.

Complex traits, including migraine, often aggregate in families, but the underlying genetic architecture behind this is not well understood. The aggregation could be explained by rare, penetrant variants that segregate according to Mendelian inheritance or by the sufficient polygenic accumulation of common variants, each with an individually small effect, or a combination of the two hypotheses. In 8,319 individuals across 1,589 migraine families, we calculated migraine polygenic risk scores (PRS) and found a significantly higher common variant burden in familial cases (n = 5,317, OR = 1.76, 95% CI = 1.71-1.81, p = 1.7 × 10-109) compared to population cases from the FINRISK cohort (n = 1,101, OR = 1.32, 95% CI = 1.25-1.38, p = 7.2 × 10-17). The PRS explained 1.6% of the phenotypic variance in the population cases and 3.5% in the familial cases (including 2.9% for migraine without aura, 5.5% for migraine with typical aura, and 8.2% for hemiplegic migraine). The results demonstrate a significant contribution of common polygenic variation to the familial aggregation of migraine.

The contribution of CACNA1A, ATP1A2 and SCN1A mutations in hemiplegic migraine: A clinical and genetic study in Finnish migraine families.

Objective To study the position of hemiplegic migraine in the clinical spectrum of migraine with aura and to reveal the importance of CACNA1A, ATP1A2 and SCN1A in the development of hemiplegic migraine in Finnish migraine families. Methods The International Classification of Headache Disorders 3rd edition criteria were used to determine clinical characteristics and occurrence of hemiplegic migraine, based on detailed questionnaires, in a Finnish migraine family collection consisting of 9087 subjects. Involvement of CACNA1A, ATP1A2 and SCN1A was studied using whole exome sequencing data from 293 patients with hemiplegic migraine. Results Overall, hemiplegic migraine patients reported clinically more severe headache and aura episodes than non-hemiplegic migraine with aura patients. We identified two mutations, c.1816G>A (p.Ala606Thr) and c.1148G>A (p.Arg383His), in ATP1A2 and one mutation, c.1994C>T (p.Thr665Met) in CACNA1A. Conclusions The results highlight hemiplegic migraine as a clinically and genetically heterogeneous disease. Hemiplegic migraine patients do not form a clearly separate group with distinct symptoms, but rather have an extreme phenotype in the migraine with aura continuum. We have shown that mutations in CACNA1A, ATP1A2 and SCN1A are not the major cause of the disease in Finnish hemiplegic migraine patients, suggesting that there are additional genetic factors contributing to the phenotype.

Pitfalls in genetic testing: the story of missed SCN1A mutations.

BACKGROUND: Sanger sequencing, still the standard technique for genetic testing in most diagnostic laboratories and until recently widely used in research, is gradually being complemented by next-generation sequencing (NGS). No single mutation detection technique is however perfect in identifying all mutations. Therefore, we wondered to what extent inconsistencies between Sanger sequencing and NGS affect the molecular diagnosis of patients. Since mutations in SCN1A, the major gene implicated in epilepsy, are found in the majority of Dravet syndrome (DS) patients, we focused on missed SCN1A mutations. METHODS: We sent out a survey to 16 genetic centers performing SCN1A testing. RESULTS: We collected data on 28 mutations initially missed using Sanger sequencing. All patients were falsely reported as SCN1A mutation-negative, both due to technical limitations and human errors. CONCLUSION: We illustrate the pitfalls of Sanger sequencing and most importantly provide evidence that SCN1A mutations are an even more frequent cause of DS than already anticipated.

Genome-wide Polygenic Burden of Rare Deleterious Variants in Sudden Unexpected Death in Epilepsy.

Sudden unexpected death in epilepsy (SUDEP) represents the most severe degree of the spectrum of epilepsy severity and is the commonest cause of epilepsy-related premature mortality. The precise pathophysiology and the genetic architecture of SUDEP remain elusive. Aiming to elucidate the genetic basis of SUDEP, we analysed rare, protein-changing variants from whole-exome sequences of 18 people who died of SUDEP, 87 living people with epilepsy and 1479 non-epilepsy disease controls. Association analysis revealed a significantly increased genome-wide polygenic burden per individual in the SUDEP cohort when compared to epilepsy (P = 5.7 × 10(- 3)) and non-epilepsy disease controls (P = 1.2 × 10(- 3)). The polygenic burden was driven both by the number of variants per individual, and over-representation of variants likely to be deleterious in the SUDEP cohort. As determined by this study, more than a thousand genes contribute to the observed polygenic burden within the framework of this study. Subsequent gene-based association analysis revealed five possible candidate genes significantly associated with SUDEP or epilepsy, but no one single gene emerges as common to the SUDEP cases. Our findings provide further evidence for a genetic susceptibility to SUDEP, and suggest an extensive polygenic contribution to SUDEP causation. Thus, an overall increased burden of deleterious variants in a highly polygenic background might be important in rendering a given individual more susceptible to SUDEP. Our findings suggest that exome sequencing in people with epilepsy might eventually contribute to generating SUDEP risk estimates, promoting stratified medicine in epilepsy, with the eventual aim of reducing an individual patient's risk of SUDEP.

Genetic analysis for a shared biological basis between migraine and coronary artery disease.

OBJECTIVE: To apply genetic analysis of genome-wide association data to study the extent and nature of a shared biological basis between migraine and coronary artery disease (CAD). METHODS: Four separate methods for cross-phenotype genetic analysis were applied on data from 2 large-scale genome-wide association studies of migraine (19,981 cases, 56,667 controls) and CAD (21,076 cases, 63,014 controls). The first 2 methods quantified the extent of overlapping risk variants and assessed the load of CAD risk loci in migraineurs. Genomic regions of shared risk were then identified by analysis of covariance patterns between the 2 phenotypes and by querying known genome-wide significant loci. RESULTS: We found a significant overlap of genetic risk loci for migraine and CAD. When stratified by migraine subtype, this was limited to migraine without aura, and the overlap was protective in that patients with migraine had a lower load of CAD risk alleles than controls. Genes indicated by 16 shared risk loci point to mechanisms with potential roles in migraine pathogenesis and CAD, including endothelial dysfunction (PHACTR1) and insulin homeostasis (GIP). CONCLUSIONS: The results suggest that shared biological processes contribute to risk of migraine and CAD, but surprisingly this commonality is restricted to migraine without aura and the impact is in opposite directions. Understanding the mechanisms underlying these processes and their opposite relationship to migraine and CAD may improve our understanding of both disorders.

Distribution and medical impact of loss-of-function variants in the Finnish founder population.

Exome sequencing studies in complex diseases are challenged by the allelic heterogeneity, large number and modest effect sizes of associated variants on disease risk and the presence of large numbers of neutral variants, even in phenotypically relevant genes. Isolated populations with recent bottlenecks offer advantages for studying rare variants in complex diseases as they have deleterious variants that are present at higher frequencies as well as a substantial reduction in rare neutral variation. To explore the potential of the Finnish founder population for studying low-frequency (0.5-5%) variants in complex diseases, we compared exome sequence data on 3,000 Finns to the same number of non-Finnish Europeans and discovered that, despite having fewer variable sites overall, the average Finn has more low-frequency loss-of-function variants and complete gene knockouts. We then used several well-characterized Finnish population cohorts to study the phenotypic effects of 83 enriched loss-of-function variants across 60 phenotypes in 36,262 Finns. Using a deep set of quantitative traits collected on these cohorts, we show 5 associations (p<5×10⁻8) including splice variants in LPA that lowered plasma lipoprotein(a) levels (P = 1.5×10⁻¹¹7). Through accessing the national medical records of these participants, we evaluate the LPA finding via Mendelian randomization and confirm that these splice variants confer protection from cardiovascular disease (OR = 0.84, P = 3×10⁻4), demonstrating for the first time the correlation between very low levels of LPA in humans with potential therapeutic implications for cardiovascular diseases. More generally, this study articulates substantial advantages for studying the role of rare variation in complex phenotypes in founder populations like the Finns and by combining a unique population genetic history with data from large population cohorts and centralized research access to National Health Registers.

Genome-wide meta-analysis identifies new susceptibility loci for migraine.

Migraine is the most common brain disorder, affecting approximately 14% of the adult population, but its molecular mechanisms are poorly understood. We report the results of a meta-analysis across 29 genome-wide association studies, including a total of 23,285 individuals with migraine (cases) and 95,425 population-matched controls. We identified 12 loci associated with migraine susceptibility (P<5×10(-8)). Five loci are new: near AJAP1 at 1p36, near TSPAN2 at 1p13, within FHL5 at 6q16, within C7orf10 at 7p14 and near MMP16 at 8q21. Three of these loci were identified in disease subgroup analyses. Brain tissue expression quantitative trait locus analysis suggests potential functional candidate genes at four loci: APOA1BP, TBC1D7, FUT9, STAT6 and ATP5B.

Genome-wide association analysis identifies susceptibility loci for migraine without aura.

Migraine without aura is the most common form of migraine, characterized by recurrent disabling headache and associated autonomic symptoms. To identify common genetic variants associated with this migraine type, we analyzed genome-wide association data of 2,326 clinic-based German and Dutch individuals with migraine without aura and 4,580 population-matched controls. We selected SNPs from 12 loci with 2 or more SNPs associated with P values of <1 × 10(-5) for replication testing in 2,508 individuals with migraine without aura and 2,652 controls. SNPs at two of these loci showed convincing replication: at 1q22 (in MEF2D; replication P = 4.9 × 10(-4); combined P = 7.06 × 10(-11)) and at 3p24 (near TGFBR2; replication P = 1.0 × 10(-4); combined P = 1.17 × 10(-9)). In addition, SNPs at the PHACTR1 and ASTN2 loci showed suggestive evidence of replication (P = 0.01; combined P = 3.20 × 10(-8) and P = 0.02; combined P = 3.86 × 10(-8), respectively). We also replicated associations at two previously reported migraine loci in or near TRPM8 and LRP1. This study identifies the first susceptibility loci for migraine without aura, thereby expanding our knowledge of this debilitating neurological disorder.

Metabonomic, transcriptomic, and genomic variation of a population cohort.

Comprehensive characterization of human tissues promises novel insights into the biological architecture of human diseases and traits. We assessed metabonomic, transcriptomic, and genomic variation for a large population-based cohort from the capital region of Finland. Network analyses identified a set of highly correlated genes, the lipid-leukocyte (LL) module, as having a prominent role in over 80 serum metabolites (of 134 measures quantified), including lipoprotein subclasses, lipids, and amino acids. Concurrent association with immune response markers suggested the LL module as a possible link between inflammation, metabolism, and adiposity. Further, genomic variation was used to generate a directed network and infer LL module's largely reactive nature to metabolites. Finally, gene co-expression in circulating leukocytes was shown to be dependent on serum metabolite concentrations, providing evidence for the hypothesis that the coherence of molecular networks themselves is conditional on environmental factors. These findings show the importance and opportunity of systematic molecular investigation of human population samples. To facilitate and encourage this investigation, the metabonomic, transcriptomic, and genomic data used in this study have been made available as a resource for the research community.