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.     

Familial Dyslexia in a Large Swedish Family: A Whole Genome Linkage Scan

There is a compelling body of evidence that developmental dyslexia runs in families and seems to be highly inheritable. Several investigations during the last two decades have shown possible locations of genes that might be involved in dyslexia, including regions of chromosomes 1, 2, 3, 6, 11, 13, 15 and 18. In addition, six candidate genes (KIAA0319, DYX1C1, DCDC2, ROBO1, MRPL19 and C2ORF3) seem to be related to dyslexia. The present study carried out a whole genome scan in a six-generation pedigree. In addition to literacy skills the assessment included cognitive skills and records concerning the history of reading and writing ability. Thirty-five percent were regarded as dyslexic in the family. A linkage analysis using both a quantitative and a qualitative approach has been performed. No evidence was obtained to support the hypothesis that the transmission of dyslexia in this pedigree is due to a highly penetrant major gene, and previous linkage findings were not replicated; however, power in this small study was not adequate to confirm linkage of genes with small to moderate effects. The results were discussed in relation to diagnostic procedures and sample characteristics.     

SNP Variations in the 7q33 Region Containing DGKI are Associated with Dyslexia in the Finnish and German Populations

Four genes, DYX1C1, ROBO1, DCDC2 and KIAA0319 have been studied both genetically and functionally as candidate genes for developmental dyslexia, a common learning disability in children. The identification of novel genes is crucial to better understand the molecular pathways affected in dyslectic individuals. Here, we report results from a fine-mapping approach involving linkage and association analysis in Finnish and German dyslexic cohorts. We restrict a candidate region to 0.3 Mb on chromosome 7q33. This region harbours the gene diacylglycerol kinase, iota (DGKI) which contains overlapping haplotypes associated with dyslexia in both Finnish and German sample sets.     

Family-based association study of DYX1C1 variants in autism

DYX1C1: was recently identified as a candidate gene for developmental dyslexia, which is characterized by an unexpected difficulty in learning to read and write despite adequate intelligence, motivation, and education. It will be important to clarify, whether the phenotype caused by DYX1C1 extends to other language-related or comorbid disorders. Impaired language development is one of the essential features in autism. Therefore, we analyzed the allelic distribution of the DYX1C1 gene by family-based association method in 100 Finnish autism families. No evidence for association was observed with any intragenic marker or with haplotypes constructed from alleles of several adjacent markers. No evidence for deviated allelic diversity was either observed: the frequency of expected dyslexia risk haplotype was comparable to its frequency in Finnish controls. Thus it seems unlikely that DYX1C1 gene would be involved in the genetic etiology of autism in Finnish patients.     

DYX1C1 functions in neuronal migration in developing neocortex

Rodent homologues of two candidate dyslexia susceptibility genes, Kiaa0319 and Dcdc2, have been shown to play roles in neuronal migration in developing cerebral neocortex. This functional role is consistent with the hypothesis that dyslexia susceptibility is increased by interference with normal neural development. In this study we report that in utero RNA interference against the rat homolog of another candidate dyslexia susceptibility gene, DYX1C1, disrupts neuronal migration in developing neocortex. The disruption of migration can be rescued by concurrent overexpression of DYX1C1, indicating that the impairment is not due to off-target effects. Transfection of C- and N-terminal truncations of DYX1C1 shows that the C-terminal TPR domains determine DYX1C1 intracellular localization to cytoplasm and nucleus. RNAi rescue experiments using truncated versions of DYX1C1 further indicate that the C-terminus of DYX1C1 is necessary and sufficient to DYX1C1's function in migration. In conclusion, DYX1C1, similar to two other candidate dyslexia susceptibility genes, functions in neuronal migration in rat neocortex.     

Haplotypes spanning SPEC2, PDZ-GEF2 and ACSL6 genes are associated with schizophrenia

Chromosome 5q22-33 is a region where studies have repeatedly found evidence for linkage to schizophrenia. In this report, we took a stepwise approach to systematically map this region in the Irish Study of High Density Schizophrenia Families (ISHDSF, 267 families, 1337 subjects) sample. We typed 289 SNPs in the critical interval of 8 million basepairs and found a 758 kb interval coding for the SPEC2/PDZ-GEF2/ACSL6 genes to be associated with the disease. Using sex and genotype-conditioned transmission disequilibrium test analyses, we found that 19 of the 24 typed markers were associated with the disease and the associations were sex-specific. We replicated these findings with an Irish case-control sample (657 cases and 414 controls), an Irish parent-proband trio sample (187 families, 564 subjects), a German nuclear family sample (211 families, 751 subjects) and a Pittsburgh nuclear family sample (247 families, 729 subjects). In all four samples, we replicated the sex-specific associations at the levels of both individual markers and haplotypes using sex- and genotype-conditioned analyses. Three risk haplotypes were identified in the five samples, and each haplotype was found in at least two samples. Consistent with the discovery of multiple estrogen-response elements in this region, our data showed that the impact of these haplotypes on risk for schizophrenia differed in males and females. From these data, we concluded that haplotypes underlying the SPEC2/PDZ-GEF2/ACSL6 region are associated with schizophrenia. However, due to the extended high LD in this region, we were unable to distinguish whether the association signals came from one or more of these genes.     

Haplotype patterns in cancer-related genes with long-range linkage disequilibrium: no evidence of association with breast cancer or positive selection

The average length of linkage disequilibrium (LD) blocks in European populations is about 22 kb. In this study, we have selected 20 genes with LD blocks larger than 60 kb (with a median length of 88 kb) from a total of 121 cancer-related genes. We observed limited haplotype diversity, with an average of three haplotypes per gene accounting for more than 90% of the diversity, two of these being a Yin-Yang pair in 95% of the LD blocks. The mean frequency of the most common haplotype in the Spanish population was just below 50%, similar to those for the HapMap CEU and African samples, but lower than the 60% observed in Asian samples. Genes involved in the regulation of nucleobases and nucleic acid metabolism were overrepresented among these 20 genes with long LD blocks (eight genes ATM, BRCA1, BRCA2, ERCC6, MLH1, MSH3, RAD54B and XRCC4) relative to the other 101 cancer-related genes studied (P=1.23 x 10(-6)). The ancestral haplotype was observed at a frequency greater than 3 in 67% of the genes either in the Spanish or one of the HapMap sampled populations. When observed, the ancestral haplotype had an average 15% frequency in the Spanish sample, less than half that observed in Asian and African samples. The Spanish Yin-Yang haplotype pair represented over 35% of haplotypes in African samples and over 65% in non-African samples. We detected differences in SNP frequencies between populations for five genes (ALDH2, APC, PIK3CB, RB1 and XRCC4, all with Fst>0.4); however, these genes did not show evidence of positive selection. Finally, we found no evidence that the haplotypes formed by SNPs in the 20 genes are associated with breast cancer.     

A comprehensive haplotype analysis of CYP19 and breast cancer risk: the Multiethnic Cohort

The CYP19 gene encodes for aromatase (P450arom), a key steroidogenic enzyme that catalyzes the final step of estrogen biosynthesis. Apart from rare mutations in CYP19 which result in severe phenotypes associated with estrogen insufficiency, little is known about whether common variation in CYP19 is associated with risk of hormone-related diseases. In this study, we employed a haplotype-based approach to search for common disease-associated variants in this candidate breast cancer susceptibility gene among African-American, Hawaiian, Japanese, Latina and White women in the Multiethnic Cohort Study (MEC). We utilized 74 densely spaced single-nucleotide polymorphisms (SNPs) (one every approximately 2.6 kb) spanning 189.4 kb of the CYP19 locus to characterize linkage disequilibrium (LD) and haplotype patterns among 69-70 individuals from each ethnic population. We detected four regions of strong LD (blocks 1-4) that were quite closely conserved across populations. Within each block there was a limited diversity of common haplotypes (5 to 10 with a frequency >/=5%) and most haplotypes were observed to be shared across populations. Twenty-five haplotype-tagging SNPs (htSNPs) were selected to predict the common haplotypes with high probability (average Rh2=0.92) and genotyped in a breast cancer case-control study in the MEC (cases, n=1355; controls, n=2580). We first performed global tests for differences in risk according to the common haplotypes and observed significant haplotype-effects in block 2 [P=0.01; haplotypes 2b (OR=1.23; 95% CI, 1.07-1.40), 2d (OR=1.28; 95% CI, 1.01-1.62)]. We also found a common long-range haplotype comprised of block-specific haplotypes 2b and 3c to be associated with increased risk of breast cancer (haplotype 2b-3c: OR=1.31; 95% CI, 1.11-1.54). Our findings suggest the hypothesis that women with the long-range CYP19 haplotype 2b-3c may be carriers of a predisposing breast cancer susceptibility allele.     

Investigation of Dyslexia and SLI Risk Variants in Reading- and Language-Impaired Subjects

Dyslexia (or reading disability) and specific language impairment (or SLI) are common childhood disorders that show considerable co-morbidity and diagnostic overlaps and have been suggested to share some genetic aetiology. Recently, genetic risk variants have been identified for SLI and dyslexia enabling the direct evaluation of possible shared genetic influences between these disorders. In this study we investigate the role of variants in these genes (namely MRPL19/C20RF3, ROBO1, DCDC2, KIAA0319, DYX1C1, CNTNAP2, ATP2C2 and CMIP) in the aetiology of SLI and dyslexia. We perform case–control and quantitative association analyses using measures of oral and written language skills in samples of SLI and dyslexic families and cases. We replicate association between KIAA0319 and DCDC2 and dyslexia and provide evidence to support a role for KIAA0319 in oral language ability. In addition, we find association between reading-related measures and variants in CNTNAP2 and CMIP in the SLI families.     

A Cohort of Balanced Reciprocal Translocations Associated with Dyslexia: Identification of Two Putative Candidate Genes at DYX1

Dyslexia is one of the most common neurodevelopmental disorders where likely many genes are involved in the pathogenesis. So far six candidate dyslexia genes have been proposed, and two of these were identified by rare chromosomal translocations in affected individuals. By systematic re-examination of all translocation carriers in Denmark, we have identified 16 different translocations associated with dyslexia. In four families, where the translocation co-segregated with the phenotype, one of the breakpoints concurred (at the cytogenetic level) with either a known dyslexia linkage region—at 15q21 (DYX1), 2p13 (DYX3) and 1p36 (DYX8)—or an unpublished linkage region at 19q13. As a first exploitation of this unique cohort, we identify three novel candidate dyslexia genes, ZNF280D and TCF12 at 15q21, and PDE7B at 6q23.3, by molecular mapping of the familial translocation with the 15q21 breakpoint.     

Evaluation of candidate genes for DYX1 and DYX2 in families with dyslexia.

Dyslexia is a common heterogeneous disorder with a significant genetic component. Multiple studies have replicated the evidence for linkage between variously defined phenotypes of dyslexia and chromosomal regions on 15q21 (DYX1) and 6p22.2 (DYX2). Based on association studies and the possibility for functional significance of several polymorphisms, candidate genes responsible for the observed linkage signal have been proposed-DYX1C1 for 15q21, and KIAA0319 and DCDC2 for 6p22.2. We investigated the evidence for contribution of these candidate genes to dyslexia in our sample of multigenerational families. Our previous quantitative linkage analyses in this dataset provided supportive evidence for linkage of dyslexia to the locus on chromosome 15, but not to the locus on chromosome 6. In the current study, we used probands from 191 families for a case control analysis, and proband-parent trios for family-based TDT analyses. The observation of weak evidence for transmission disequilibrium for one of the two studied polymorphisms in DYX1C1 suggests involvement of this gene in dyslexia in our dataset. We did not find evidence for the association of KIAA0319 or DCDC2 alleles to dyslexia in our sample. We observed a slight tendency for an intronic deletion in DCDC2 to be associated with worse performance on some quantitative measures of dyslexia in the probands in our sample, but not in their parents.     

Mapping for dyslexia and related cognitive trait loci provides strong evidence for further risk genes on chromosome 6p21

In a genome‐wide linkage scan, we aimed at mapping risk loci for dyslexia in the German population. Our sample comprised 1,030 individuals from 246 dyslexia families which were recruited through a single‐proband sib pair study design and a detailed assessment of dyslexia and related cognitive traits. We found evidence for a major dyslexia locus on chromosome 6p21. The cognitive trait rapid naming (objects/colors) produced a genome‐wide significant LOD score of 5.87 (P = 1.00 × 10^?7) and the implicated 6p‐risk region spans around 10 Mb. Although our finding maps close to DYX2, where the dyslexia candidate genes DCDC2 and KIAA0319 have already been identified, our data point to the presence of an additional risk gene in this region and are highlighting the impact of 6p21 in dyslexia and related cognitive traits. © 2010 Wiley‐Liss, Inc.     

A new gene (DYX3) for dyslexia is located on chromosome 2

Developmental dyslexia is a specific reading disability affecting children and adults who otherwise possess normal intelligence, cognitive skills, and adequate schooling. Difficulties in spelling and reading may persist through adult life. Possible localisations of genes for dyslexia have been reported on chromosomes 15 (DYX1), 6p21.3-23 (DYX2), and 1p over the last 15 years. Only the localisation to 6p21.3-23 has been clearly confirmed and a genome search has not previously been carried out. We have investigated a large Norwegian family in which dyslexia is inherited as an autosomal dominant trait. A genome wide search for linkage with an average 20 cM marker density was initiated in 36 of the 80 family members. The linkage analysis was performed under three different diagnostic models. Linkage analysis in the family identified a region in 2p15-p16 which cosegregated with dyslexia. Maximum lod scores of 3.54, 2.92, and 4.32 for the three different diagnostic models were obtained. These results were confirmed by a non-parametric multipoint GENEHUNTER analysis in which the most likely placement of the gene was in a 4 cM interval between markers D2S2352 and D2S1337. Localisation of a gene for dyslexia to 2p15-16, together with the confirmed linkage to 6p21.3-23, constitute strong evidence for genetic heterogeneity in dyslexia. Since no gene for dyslexia has been isolated, little is known about the molecular processes involved. The isolation and molecular characterisation of this newly reported gene on chromosome 2 (DYX3) and DYX1 will thus provide new and exciting insights into the processes involved in reading and spelling.     

Immunoglobulin constant heavy chain gene polymorphisms of Caucasoids of west European origin.

We have used restriction fragment length polymorphism (RFLP) analysis to investigate the immunoglobulin constant heavy chain (IgCH) loci and the associated locus Dl4Sl, of Caucasoids from South East England and South-West. West Germany, Haplotypes were determined using probes to the Ig heavy chain switch loci S mu and S alpha 1, the IgC gamma 3 and IgC gamma 2 loci as well as the Dl4Sl locus which is 3' of the IgCH loci. The 6.3:1.7 kilobase (kb) Bst EII C gamma 3-C gamma 2 haplotype was the most prevalent in the population from South-East England (frequency 0.364), whilst the 2.3:3.7 kb C gamma 3-C gamma 2 haplotype was the major haplotype in the German population (frequency 0.400). With one exception, the major haplotypes of these two populations differed from the ones previously published for a Caucasoid population from California. This suggests that there may be a major ancestral IgCH haplotype which has been maintained in the population, whilst other haplotypes tend to be specific for a particular group of Caucasoids.     

The structure of the tau haplotype in controls and in progressive supranuclear palsy

The group of neurodegenerative diseases collectively known as tauopathies are characterized by hallmark lesions consisting of fibrillar aggregates of the microtubule-associated protein, tau (MAPT). Mutations of the tau gene (MAPT) are the cause of frontotemporal dementia with parkinsonism linked to chromosome 17, giving tau a central role in the pathogenic process. The chromosomal region containing MAPT has been shown to evolve into two major haplotypes, H1 and H2, which are defined by linkage disequilibrium (LD) between several polymorphisms over the entire MAPT gene. Studies to date suggest a complete absence of recombination between these two haplotypes. The more common haplotype H1 is over-represented in patients with progressive supranuclear palsy (PSP) and corticobasal degeneration. Using single nucleotide polymorphisms, we mapped LD in the regions flanking MAPT and have established the maximum extent of the haplotype block on chromosome 17q21.31 as a region covering approximately 2 Mb. This gene-rich region extends centromerically beyond the corticotrophin releasing hormone receptor 1 gene (CRHR1) to a region of approximately 400 kb, where there is a complete loss of LD. The telomeric end is defined by an approximately 150 kb region just beyond the WNT3 gene. We show that the entire, fully extended H1 haplotype is associated with PSP, which implicates several other genes in addition to MAPT, as candidate pathogenic loci.     

Mutational risk in highly repetitive exon ORF15 of the RPGR multidisease gene is not associated with haplotype background

Exon ORF15 is an alternative exon in the retinitis pigmentosa GTPase regulator (RPGR) gene containing a highly repetitive, purine-rich internal region. It constitutes a mutational hot spot giving rise to a group of heterogeneous X-linked retinal disorders. We sought to determine whether non-pathogenic substitutions and sequence length variations in the repetitive sequence have an influence on the risk of pathogenic exon ORF15 mutations. The type and distribution of exon ORF15 polymorphisms were assessed by genotyping 107 healthy German males using standard procedures. Polymorphisms were grouped into haplotypes and their frequencies determined in normal controls and previously analyzed patients with X-linked retinitis pigmentosa (XLRP). In the control group we identified 6 complex variants of the most common, ancestral exon ORF15 allele corresponding to the GenBank reference sequence (accession no. AF286472). Exon ORF15 mutations in XLRP patients were associated with the ancestral allele in 75% of affected cases. Four of the most recent founder haplotypes termed H3, H4, H6 and H7 were not identified in the patient samples. Our analysis and review of polymorphism data from the published literature suggests the presence of common exon ORF15 haplotypes in the European population. While the mutational risk in the RPGR gene appears not to be altered by the haplotype background, exon ORF15 haplotype analysis may be useful for tracing the evolutionary history of RP3-associated diseases.     

Haplotypes in the APOA1-C3-A4-A5 gene cluster affect plasma lipids in both humans and baboons

Genetic studies in non-human primates serve as a potential strategy for identifying genomic intervals where polymorphisms impact upon human disease-related phenotypes. It remains unclear, however, whether independently arising polymorphisms in orthologous regions of non-human primates leads to similar variation in a quantitative trait found in both species. To explore this paradigm, we studied a baboon apolipoprotein gene cluster (APOA1/C3/A4/A5) for which the human gene orthologs have well-established roles in influencing plasma HDL-cholesterol and triglyceride concentrations. Our extensive polymorphism analysis of this 68 kb gene cluster in 96 pedigreed baboons identified several haplotype blocks each with limited diversity, consistent with haplotype findings in humans. To determine whether baboons, like humans, also have particular haplotypes associated with lipid phenotypes, we genotyped 634 well-characterized baboons using 16 haplotype tagging SNPs. Genetic analysis of single SNPs, as well as haplotypes, revealed an association of APOA5 and APOC3 variants with HDL-cholesterol and triglyceride concentrations, respectively. Thus, independent variation in orthologous genomic intervals does associate with similar quantitative lipid traits in both species, supporting the possibility of uncovering human quantitative trait loci genes in a highly controlled non-human primate model.     

Linkage disequilibrium and haplotype diversity in the genes of the renin-angiotensin system: findings from the family blood pressure program

Association studies of candidate genes with complex traits have generally used one or a few single nucleotide polymorphisms (SNPs), although variation in the extent of linkage disequilibrium (LD) within genes markedly influences the sensitivity and precision of association studies. The extent of LD and the underlying haplotype structure for most candidate genes are still unavailable. We sampled 193 blacks (African-Americans) and 160 whites (European-Americans) and estimated the intragenic LD and the haplotype structure in four genes of the renin-angiotensin system. We genotyped 25 SNPs, with all but one of the pairs spaced between 1 and 20 kb, thus providing resolution at small scale. The pattern of LD within a gene was very heterogeneous. Using a robust method to define haplotype blocks, blocks of limited haplotype diversity were identified at each locus; between these blocks, LD was lost owing to the history of recombination events. As anticipated, there was less LD among blacks, the number of haplotypes was substantially larger, and shorter haplotype segments were found, compared with whites. These findings have implications for candidate-gene association studies and indicate that variation between populations of European and African origin in haplotype diversity is characteristic of most genes.     

Sequence variation, linkage disequilibrium and association with Crohn's disease on chromosome 5q31

Chromosome 5q31 contains a cluster of genes involved in immune response, including a 250 kb risk haplotype associated with Crohn's disease (CD) susceptibility. Recently, two functional variants in SLC22A4 and SLC22A5 (L503F and G-207C), encoding the cation transporters OCTN1 and OCTN2, were proposed as causal variants for CD, but with conflicting genetic evidence regarding their contribution. We investigated this locus by resequencing the coding regions of 10 genes in 24 CD cases and deriving a linkage disequilibrium (LD) map of the 27 single nucleotide polymorphisms (SNPs) detected. Ten SNPs representative of the LD groups observed, were tested for CD association. L503F in SLC22A4 was the only nonsynonymous SNP significantly associated with CD (P=0.003), but was not associated with disease in the absence of other markers of the 250 kb risk haplotype. Two other SNPs, rs11242115 in IRF1 and rs17166050 in RAD50, lying outside the 250 kb risk haplotype, also showed CD association (P=0.019 and P=0.0080, respectively). The RAD50 gene contains a locus control region regulating expression of the Th2 cytokine genes at this locus. Other as yet undiscovered SNPs in this region may therefore modulate gene expression and contribute to the risk of CD, and perhaps of other inflammatory phenotypes.