Jessica Becker

Genome-wide association study identifies two susceptibility loci for nonsyndromic cleft lip with or without cleft palate

We conducted a genome-wide association study for nonsyndromic cleft lip with or without cleft palate (NSCL/P) in 401 affected individuals and 1,323 controls, with replication in an independent sample of 793 NSCL/P triads. We report two new loci associated with NSCL/P at 17q22 (rs227731, combined P = 1.07 × 10−8, relative risk in homozygotes = 1.84, 95% CI 1.34–2.53) and 10q25.3 (rs7078160, combined P = 1.92 × 10−8, relative risk in homozygotes = 2.17, 95% CI 1.32–3.56).

Genome-wide meta-analyses of nonsyndromic cleft lip with or without cleft palate identify six new risk loci

We have conducted the first meta-analyses for nonsyndromic cleft lip with or without cleft palate (NSCL/P) using data from the two largest genome-wide association studies published to date. We confirmed associations with all previously identified loci and identified six additional susceptibility regions (1p36, 2p21, 3p11.1, 8q21.3, 13q31.1 and 15q22). Analysis of phenotypic variability identified the first specific genetic risk factor for NSCLP (nonsyndromic cleft lip plus palate) (rs8001641; PNSCLP = 6.51 × 10−11; homozygote relative risk = 2.41, 95% confidence interval (CI) 1.84–3.16).

First genome-wide association scan on neurophysiological endophenotypes points to trans-regulation effects on SLC2A3 in dyslexic children

Dyslexia is one of the most common learning disorders affecting about 5% of all school-aged children. It has been shown that event-related potential measurements reveal differences between dyslexic children and age-matched controls. This holds particularly true for mismatch negativity (MMN), which reflects automatic speech deviance processing and is altered in dyslexic children. We performed a whole-genome association analysis in 200 dyslexic children, focusing on MMN measurements. We identified rs4234898, a marker located on chromosome 4q32.1, to be significantly associated with the late MMN component. This association could be replicated in an independent second sample of 186 dyslexic children, reaching genome-wide significance in the combined sample (P=5.14e−08). We also found an association between the late MMN component and a two-marker haplotype of rs4234898 and rs11100040, one of its neighboring single nucleotide polymorphisms (SNPs). In the combined sample, this marker combination withstands correction for multiple testing (P=6.71e−08). Both SNPs lie in a region devoid of any protein-coding genes; however, they both show significant association with mRNA-expression levels of SLC2A3 on chromosome 12, the predominant facilitative glucose transporter in neurons. Our results suggest a possible trans-regulation effect on SLC2A3, which might lead to glucose deficits in dyslexic children and could explain their attenuated MMN in passive listening tasks.

Widespread non-additive and interaction effects within HLA loci modulate the risk of autoimmune diseases.

Human leukocyte antigen (HLA) genes confer substantial risk for autoimmune diseases on a log-additive scale. Here we speculated that differences in autoantigen-binding repertoires between a heterozygotes two expressed HLA variants might result in additional non-additive risk effects. We tested the non-additive disease contributions of classical HLA alleles in patients and matched controls for five common autoimmune diseases: rheumatoid arthritis (ncases = 5,337), type 1 diabetes (T1D; ncases = 5,567), psoriasis vulgaris (ncases = 3,089), idiopathic achalasia (ncases = 727) and celiac disease (ncases = 11,115). In four of the five diseases, we observed highly significant, non-additive dominance effects (rheumatoid arthritis, P = 2.5 × 10−12; T1D, P = 2.4 × 10−10; psoriasis, P = 5.9 × 10−6; celiac disease, P = 1.2 × 10−87). In three of these diseases, the non-additive dominance effects were explained by interactions between specific classical HLA alleles (rheumatoid arthritis, P = 1.8 × 10−3; T1D, P = 8.6 × 10−27; celiac disease, P = 6.0 × 10−100). These interactions generally increased disease risk and explained moderate but significant fractions of phenotypic variance (rheumatoid arthritis, 1.4%; T1D, 4.0%; celiac disease, 4.1%) beyond a simple additive model.

Genetic risk factors for nonsyndromic cleft lip with or without cleft palate in a Mesoamerican population: Evidence for IRF6 and variants at 8q24 and 10q25

INTRODUCTION Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is one of the most common of all birth defects. NSCL/P has a multifactorial etiology that includes both genetic and environmental factors. The IRF6 gene and three further susceptibility loci at 8q24, 10q25, and 17q22, which were identified by a recent genome-wide association scan (GWAS), are confirmed genetic risk factors for NSCL/P in patients of European descent. METHODS A case-control association study was performed to investigate whether these four risk loci contribute to NSCL/P in a Mesoamerican population using four single nucleotide polymorphisms to represent IRF6 and the three novel susceptibility loci. A total of 149 NSCL/P patients and 303 controls of Mayan origin were included. RESULTS Single marker analysis revealed a significant association between NSCL/P and risk variants in IRF6 and the 8q24 and 10q25 loci. In contrast to previous findings, the association at the 8q24 locus was driven solely by homozygote carriers of the risk allele. This suggests that this locus might act in a recessive manner in the Mayan population. No evidence for association was found at the 17q22 locus. This may have been attributable to the limited power of the sample. CONCLUSION These results suggest that IRF6 and the 10q25 and 8q24 loci confer a risk for the development of NSCL/P in persons of Mayan origin.

Follow-up of loci from the International Genomics of Alzheimer's Disease Project identifies TRIP4 as a novel susceptibility gene

To follow-up loci discovered by the International Genomics of Alzheimer’s Disease Project, we attempted independent replication of 19 single nucleotide polymorphisms (SNPs) in a large Spanish sample (Fundació ACE data set; 1808 patients and 2564 controls). Our results corroborate association with four SNPs located in the genes INPP5D, MEF2C, ZCWPW1 and FERMT2, respectively. Of these, ZCWPW1 was the only SNP to withstand correction for multiple testing (P=0.000655). Furthermore, we identify TRIP4 (rs74615166) as a novel genome-wide significant locus for Alzheimer’s disease risk (odds ratio=1.31; confidence interval 95% (1.19–1.44); P=9.74 × 10−9).

Common variants in the HLA-DQ region confer susceptibility to idiopathic achalasia

Idiopathic achalasia is characterized by a failure of the lower esophageal sphincter to relax due to a loss of neurons in the myenteric plexus. This ultimately leads to massive dilatation and an irreversibly impaired megaesophagus. We performed a genetic association study in 1,068 achalasia cases and 4,242 controls and fine-mapped a strong MHC association signal by imputing classical HLA haplotypes and amino acid polymorphisms. An eight-residue insertion at position 227–234 in the cytoplasmic tail of HLA-DQβ1 (encoded by HLA-DQB1*05:03 and HLA-DQB1*06:01) confers the strongest risk for achalasia (P = 1.73 × 10−19). In addition, two amino acid substitutions in the extracellular domain of HLA-DQα1 at position 41 (lysine encoded by HLA-DQA1*01:03; P = 5.60 × 10−10) and of HLA-DQβ1 at position 45 (glutamic acid encoded by HLA-DQB1*03:01 and HLA-DQB1*03:04; P = 1.20 × 10−9) independently confer achalasia risk. Our study implies that immune-mediated processes are involved in the pathophysiology of achalasia.

Characterizing the genetic basis of innate immune response in TLR4-activated human monocytes

Toll-like receptors (TLRs) play a key role in innate immunity. Apart from their function in host defense, dysregulation in TLR signalling can confer risk to autoimmune diseases, septic shock or cancer. Here we report genetic variants and transcripts that are active only during TLR signalling and contribute to interindividual differences in immune response. Comparing unstimulated versus TLR4-stimulated monocytes reveals 1,471 expression quantitative trait loci (eQTLs) that are unique to TLR4 stimulation. Among these we find functional SNPs for the expression of NEU4, CCL14, CBX3 and IRF5 on TLR4 activation. Furthermore, we show that SNPs conferring risk to primary biliary cirrhosis (PBC), inflammatory bowel disease (IBD) and celiac disease are immune response eQTLs for PDGFB and IL18R1. Thus, PDGFB and IL18R1 represent plausible candidates for studying the pathophysiology of these disorders in the context of TLR4 activation. In summary, this study presents novel insights into the genetic basis of the innate immune response and exemplifies the value of eQTL studies in the context of exogenous cell stimulation.

Genetic analysis of dyslexia candidate genes in the European cross-linguistic NeuroDys cohort

Dyslexia is one of the most common childhood disorders with a prevalence of around 5–10% in school-age children. Although an important genetic component is known to have a role in the aetiology of dyslexia, we are far from understanding the molecular mechanisms leading to the disorder. Several candidate genes have been implicated in dyslexia, including DYX1C1, DCDC2, KIAA0319, and the MRPL19/C2ORF3 locus, each with reports of both positive and no replications. We generated a European cross-linguistic sample of school-age children – the NeuroDys cohort – that includes more than 900 individuals with dyslexia, sampled with homogenous inclusion criteria across eight European countries, and a comparable number of controls. Here, we describe association analysis of the dyslexia candidate genes/locus in the NeuroDys cohort. We performed both case–control and quantitative association analyses of single markers and haplotypes previously reported to be dyslexia-associated. Although we observed association signals in samples from single countries, we did not find any marker or haplotype that was significantly associated with either case–control status or quantitative measurements of word-reading or spelling in the meta-analysis of all eight countries combined. Like in other neurocognitive disorders, our findings underline the need for larger sample sizes to validate possibly weak genetic effects.

A systematic eQTL study of cis-trans epistasis in 210 HapMap individuals

We aimed at identifying transcripts whose expression is regulated by a SNP–SNP interaction. Out of 47 294 expression phenotypes we used 3107 transcripts that survived an extensive quality control and 86 613 linkage disequilibrium-pruned SNP markers that have been genotyped in 210 individuals. For each transcript we defined cis-SNPs, tested them for epistasis with all trans-SNPs, and corrected all observed cis–trans-regulated expression effects for multiple testing. We determined that the expression of about 15% of all included transcripts is regulated by a significant two-locus interaction, which is more than expected (P=2.86 × 10−144). Our findings suggest further that cis-markers with so called ‘marginal effects’ are more likely to be involved in two-locus gene regulation than expected (P=8.27 × 10−05), although the majority of interacting cis-markers showed no one-locus regulation. Furthermore, we found evidence that gene-mediated trans-effects are not a major source of epistasis, as no enrichment of genes has been found in close vicinity of trans-SNPs. In addition, our data support the notion that neither chromosomal regions nor cellular processes are enriched in epistatic interactions. Finally, some of the cis–trans regulated genes have been found in genome-wide association studies, which might be interesting for follow-up studies of the corresponding disorders. In summary, our results provide novel insights into the complex genome-transcriptome regulation.