Michael Nothnagel

Sequence variants in IL10, ARPC2 and multiple other loci contribute to ulcerative colitis susceptibility

Inflammatory bowel disease (IBD) typically manifests as either ulcerative colitis (UC) or Crohns disease (CD). Systematic identification of susceptibility genes for IBD has thus far focused mainly on CD, and little is known about the genetic architecture of UC. Here we report a genome-wide association study with 440,794 SNPs genotyped in 1,167 individuals with UC and 777 healthy controls. Twenty of the most significantly associated SNPs were tested for replication in three independent European case-control panels comprising a total of 1,855 individuals with UC and 3,091 controls. Among the four consistently replicated markers, SNP rs3024505 immediately flanking the IL10 (interleukin 10) gene on chromosome 1q32.1 showed the most significant association in the combined verification samples (P = 1.35 × 10−12; OR = 1.46 (1.31–1.62)). The other markers were located in ARPC2 and in the HLA-BTNL2 region. Association between rs3024505 and CD (1,848 cases, 1,804 controls) was weak (P = 0.013; OR = 1.17 (1.01–1.34)). IL10 is an immunosuppressive cytokine that has long been proposed to influence IBD pathophysiology. Our findings strongly suggest that defective IL10 function is central to the pathogenesis of the UC subtype of IBD.

15q13.3 microdeletions increase risk of idiopathic generalized epilepsy

We identified 15q13.3 microdeletions encompassing the CHRNA7 gene in 12 of 1,223 individuals with idiopathic generalized epilepsy (IGE), which were not detected in 3,699 controls (joint P = 5.32 × 10−8). Most deletion carriers showed common IGE syndromes without other features previously associated with 15q13.3 microdeletions, such as intellectual disability, autism or schizophrenia. Our results indicate that 15q13.3 microdeletions constitute the most prevalent risk factor for common epilepsies identified to date.

Correlation between Genetic and Geographic Structure in Europe

Understanding the genetic structure of the European population is important, not only from a historical perspective, but also for the appropriate design and interpretation of genetic epidemiological studies. Previous population genetic analyses with autosomal markers in Europe either had a wide geographic but narrow genomic coverage [1, 2], or vice versa [3-6]. We therefore investigated Affymetrix GeneChip 500K genotype data from 2,514 individuals belonging to 23 different subpopulations, widely spread over Europe. Although we found only a low level of genetic differentiation between subpopulations, the existing differences were characterized by a strong continent-wide correlation between geographic and genetic distance. Furthermore, mean heterozygosity was larger, and mean linkage disequilibrium smaller, in southern as compared to northern Europe. Both parameters clearly showed a clinal distribution that provided evidence for a spatial continuity of genetic diversity in Europe. Our comprehensive genetic data are thus compatible with expectations based upon European population history, including the hypotheses of a south-north expansion and/or a larger effective population size in southern than in northern Europe. By including the widely used CEPH from Utah (CEU) samples into our analysis, we could show that these individuals represent northern and western Europeans reasonably well, thereby confirming their assumed regional ancestry.

Power and sample size calculations for case-control genetic association tests when errors are present: application to single nucleotide polymorphisms.

The purpose of this work is to quantify the effects that errors in genotyping have on power and the sample size necessary to maintain constant asymptotic Type I and Type II error rates (SSN) for case-control genetic association studies between a disease phenotype and a di-allelic marker locus, for example a single nucleotide polymorphism (SNP) locus. We consider the effects of three published models of genotyping errors on the chi-square test for independence in the 2 × 3 table. After specifying genotype frequencies for the marker locus conditional on disease status and error model in both a genetic model-based and a genetic model-free framework, we compute the asymptotic power to detect association through specification of the test’s non-centrality parameter. This parameter determines the functional dependence of SSN on the genotyping error rates. Additionally, we study the dependence of SSN on linkage disequilibrium (LD), marker allele frequencies, and genotyping error rates for a dominant disease model. Increased genotyping error rate requires a larger SSN. Every 1% increase in sum of genotyping error rates requires that both case and control SSN be increased by 2–8%, with the extent of increase dependent upon the error model. For the dominant disease model, SSN is a nonlinear function of LD and genotyping error rate, with greater SSN for lower LD and higher genotyping error rate. The combination of lower LD and higher genotyping error rates requires a larger SSN than the sum of the SSN for the lower LD and for the higher genotyping error rate.

Genome-Wide Association Analysis in Primary Sclerosing Cholangitis

BACKGROUND & AIMS We aimed to characterize the genetic susceptibility to primary sclerosing cholangitis (PSC) by means of a genome-wide association analysis of single nucleotide polymorphism (SNP) markers. METHODS A total of 443,816 SNPs on the Affymetrix SNP Array 5.0 (Affymetrix, Santa Clara, CA) were genotyped in 285 Norwegian PSC patients and 298 healthy controls. Associations detected in this discovery panel were re-examined in independent case-control panels from Scandinavia (137 PSC cases and 368 controls), Belgium/The Netherlands (229 PSC cases and 735 controls), and Germany (400 cases and 1832 controls). RESULTS The strongest associations were detected near HLA-B at chromosome 6p21 (rs3099844: odds ratio [OR], 4.8; 95% confidence interval [CI], 3.6-6.5; P = 2.6 x 10(-26); and rs2844559: OR, 4.7; 95% CI, 3.5-6.4; P = 4.2 x 10(-26) in the discovery panel). Outside the HLA complex, rs9524260 at chromosome 13q31 showed significant associations in 3 of 4 study panels. Lentiviral silencing of glypican 6, encoded at this locus, led to the up-regulation of proinflammatory markers in a cholangiocyte cell line. Of 15 established ulcerative colitis susceptibility loci, significant replication was obtained at chromosomes 2q35 and 3p21 (rs12612347: OR, 1.26; 95% CI, 1.06-1.50; and rs3197999: OR, 1.22; 95% CI, 1.02-1.47, respectively), with circumstantial evidence supporting the G-protein-coupled bile acid receptor 1 and macrophage-stimulating 1, respectively, as the likely disease genes. CONCLUSIONS Strong HLA associations and a subset of genes involved in bile homeostasis and other inflammatory conditions constitute key components of the genetic architecture of PSC.

Mutations in GRIN2A cause idiopathic focal epilepsy with rolandic spikes

Idiopathic focal epilepsy (IFE) with rolandic spikes is the most common childhood epilepsy, comprising a phenotypic spectrum from rolandic epilepsy (also benign epilepsy with centrotemporal spikes, BECTS) to atypical benign partial epilepsy (ABPE), Landau-Kleffner syndrome (LKS) and epileptic encephalopathy with continuous spike and waves during slow-wave sleep (CSWS). The genetic basis is largely unknown. We detected new heterozygous mutations in GRIN2A in 27 of 359 affected individuals from 2 independent cohorts with IFE (7.5%; P = 4.83 × 10−18, Fishers exact test). Mutations occurred significantly more frequently in the more severe phenotypes, with mutation detection rates ranging from 12/245 (4.9%) in individuals with BECTS to 9/51 (17.6%) in individuals with CSWS (P = 0.009, Cochran-Armitage test for trend). In addition, exon-disrupting microdeletions were found in 3 of 286 individuals (1.0%; P = 0.004, Fishers exact test). These results establish alterations of the gene encoding the NMDA receptor NR2A subunit as a major genetic risk factor for IFE.

Identification of a shared genetic susceptibility locus for coronary heart disease and periodontitis.

Recent studies indicate a mutual epidemiological relationship between coronary heart disease (CHD) and periodontitis. Both diseases are associated with similar risk factors and are characterized by a chronic inflammatory process. In a candidate-gene association study, we identify an association of a genetic susceptibility locus shared by both diseases. We confirm the known association of two neighboring linkage disequilibrium regions on human chromosome 9p21.3 with CHD and show the additional strong association of these loci with the risk of aggressive periodontitis. For the lead SNP of the main associated linkage disequilibrium region, rs1333048, the odds ratio of the autosomal-recessive mode of inheritance is 1.99 (95% confidence interval 1.33–2.94; P = 6.9×10−4) for generalized aggressive periodontitis, and 1.72 (1.06–2.76; P = 2.6×10−2) for localized aggressive periodontitis. The two associated linkage disequilibrium regions map to the sequence of the large antisense noncoding RNA ANRIL, which partly overlaps regulatory and coding sequences of CDKN2A/CDKN2B. A closely located diabetes-associated variant was independent of the CHD and periodontitis risk haplotypes. Our study demonstrates that CHD and periodontitis are genetically related by at least one susceptibility locus, which is possibly involved in ANRIL activity and independent of diabetes associated risk variants within this region. Elucidation of the interplay of ANRIL transcript variants and their involvement in increased susceptibility to the interactive diseases CHD and periodontitis promises new insight into the underlying shared pathogenic mechanisms of these complex common diseases.

A genome-wide association study confirms APOE as the major gene influencing survival in long-lived individuals

We conducted a case-control genome-wide association study (GWAS) of human longevity, comparing 664,472 autosomal SNPs in 763 long-lived individuals (LLI; mean age: 99.7 years) and 1085 controls (mean age: 60.2 years) from Germany. Only one association, namely that of SNP rs4420638 near the APOC1 gene, achieved genome-wide significance (allele-based P=1.8×10(-10)). However, logistic regression analysis revealed that this association, which was replicated in an independent German sample, is fully explicable by linkage disequilibrium with the APOE allele ɛ4, the only variant hitherto established as a major genetic determinant of survival into old age. Our GWAS failed to identify any additional autosomal susceptibility genes. One explanation for this lack of success in our study would be that GWAS provide only limited statistical power for a polygenic phenotype with loci of small effect such as human longevity. A recent GWAS in Dutch LLI independently confirmed the APOE-longevity association, thus strengthening the conclusion that this locus is a very, if not the most, important genetic factor influencing longevity.

Wnt Signaling and Dupuytren's Disease

BACKGROUND Dupuytrens disease is a benign fibromatosis of the hands and fingers that leads to flexion contractures. We hypothesized that multiple genetic and environmental factors influence susceptibility to this disease and sought to identify susceptibility genes to better understand its pathogenesis. METHODS We conducted a genomewide association study of 960 Dutch persons with Dupuytrens disease and 3117 controls (the discovery set) to test for association between the disease and genetic markers. We tested the 35 single-nucleotide polymorphisms (SNPs) most strongly associated with Dupuytrens disease (P<1×10(-4)) in the discovery set in three additional, independent case series comprising a total of 1365 affected persons and 8445 controls from Germany, the United Kingdom, and The Netherlands. RESULTS Initially, we observed a significant genomewide association between Dupuytrens disease and 8 SNPs at three loci. Tests of replication and joint analysis of all data from 2325 patients with Dupuytrens disease and 11,562 controls yielded an association with 11 SNPs from nine different loci (P<5.0×10(-8)). Six of these loci contain genes known to be involved in the Wnt-signaling pathway: WNT4 (rs7524102) (P=2.8×10(-9); odds ratio, 1.28), SFRP4 (rs16879765) (P=5.6×10(-39); odds ratio, 1.98), WNT2 (rs4730775) (P=3.0×10(-8); odds ratio, 0.83), RSPO2 (rs611744) (P=7.9×10(-15); odds ratio, 0.75), SULF1 (rs2912522) (P=2.0×10(-13); odds ratio, 0.72), and WNT7B (rs6519955) (P=3.2×10(-33); odds ratio, 1.54). CONCLUSIONS This study implicates nine different loci involved in genetic susceptibility to Dupuytrens disease. The fact that six of these nine loci harbor genes encoding proteins in the Wnt-signaling pathway suggests that aberrations in this pathway are key to the process of fibromatosis in Dupuytrens disease.

De novo loss- or gain-of-function mutations in KCNA2 cause epileptic encephalopathy.

Epileptic encephalopathies are a phenotypically and genetically heterogeneous group of severe epilepsies accompanied by intellectual disability and other neurodevelopmental features. Using next-generation sequencing, we identified four different de novo mutations in KCNA2, encoding the potassium channel KV1.2, in six isolated patients with epileptic encephalopathy (one mutation recurred three times independently). Four individuals presented with febrile and multiple afebrile, often focal seizure types, multifocal epileptiform discharges strongly activated by sleep, mild to moderate intellectual disability, delayed speech development and sometimes ataxia. Functional studies of the two mutations associated with this phenotype showed almost complete loss of function with a dominant-negative effect. Two further individuals presented with a different and more severe epileptic encephalopathy phenotype. They carried mutations inducing a drastic gain-of-function effect leading to permanently open channels. These results establish KCNA2 as a new gene involved in human neurodevelopmental disorders through two different mechanisms, predicting either hyperexcitability or electrical silencing of KV1.2-expressing neurons.