Independent studies have shown that candidate genes for dyslexia and specific language impairment (SLI) impact upon reading/language-specific traits in the general population. To further explore the effect of disorder-associated genes on cognitive functions, we investigated whether they play a role in broader cognitive traits. We tested a panel of dyslexia and SLI genetic risk factors for association with two measures of general cognitive abilities, or IQ, (verbal and non-verbal) in the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort (N>5,000). Only the MRPL19/C2ORF3 locus showed statistically significant association (minimum P = 0.00009) which was further supported by independent replications following analysis in four other cohorts. In addition, a fifth independent sample showed association between the MRPL19/C2ORF3 locus and white matter structure in the posterior part of the corpus callosum and cingulum, connecting large parts of the cortex in the parietal, occipital and temporal lobes. These findings suggest that this locus, originally identified as being associated with dyslexia, is likely to harbour genetic variants associated with general cognitive abilities by influencing white matter structure in localised neuronal regions.

Recent advances in the field of language-related disorders have led to the identification of candidate genes for specific language impairment (SLI) and dyslexia. Replication studies have been conducted in independent samples including population-based cohorts, which can be characterised for a large number of relevant cognitive measures. The availability of a wide range of phenotypes allows us to not only identify the most suitable traits for replication of genetic association but also to refine the associated cognitive trait. In addition, it is possible to test for pleiotropic effects across multiple phenotypes which could explain the extensive comorbidity observed across SLI, dyslexia and other neurodevelopmental disorders. The availability of genome-wide genotype data for such cohorts will facilitate this kind of analysis but important issues, such as multiple test corrections, have to be taken into account considering that small effect sizes are expected to underlie such associations.

Recent advances in the field of language-related disorders have led to the identification of candidate genes for specific language impairment (SLI) and dyslexia. Replication studies have been conducted in independent samples including population-based cohorts, which can be characterised for a large number of relevant cognitive measures. The availability of a wide range of phenotypes allows us to not only identify the most suitable traits for replication of genetic association but also to refine the associated cognitive trait. In addition, it is possible to test for pleiotropic effects across multiple phenotypes which could explain the extensive comorbidity observed across SLI, dyslexia and other neurodevelopmental disorders. The availability of genome-wide genotype data for such cohorts will facilitate this kind of analysis but important issues, such as multiple test corrections, have to be taken into account considering that small effect sizes are expected to underlie such associations.

Background

Several susceptibility genes have been proposed for dyslexia (reading disability; RD) and specific language impairment (SLI). RD and SLI show comorbidity, but it is unclear whether a common genetic component is shared.

Methods

We have investigated whether candidate genes for RD and SLI affect specific cognitive traits or have broad effect on cognition. We have analyzed common risk variants within RD (MRPL19/C2ORF3, KIAA0319, and DCDC2) and language impairment (CMIP and ATP2C2) candidate loci in the Avon Longitudinal Study of Parents and Children cohort (n = 3725), representing children born in southwest England in the early 1990s.

Results

We detected associations between reading skills and KIAA0319, DCDC2, and CMIP. We show that DCDC2 is specifically associated with RD, whereas variants in CMIP and KIAA0319 are associated with reading skills across the ability range. The strongest associations were restricted to single-word reading and spelling measures, suggesting that these genes do not extend their effect to other reading and language-related skills. Inclusion of individuals with comorbidity tends to strengthen these associations. Our data do not support MRPL19/C2ORF3 as a locus involved in reading abilities nor CMIP/ATP2C2 as genes regulating language skills.

Conclusions

We provide further support for the role of KIAA0319 and DCDC2 in contributing to reading abilities and novel evidence that the language-disorder candidate gene CMIP is also implicated in reading processes. Additionally, we present novel data to evaluate the prevalence and comorbidity of RD and SLI, and we recommend not excluding individuals with comorbid RD and SLI when designing genetic association studies for RD.

Approximately 90% of humans are right-handed. Handedness is a heritable trait, yet the genetic basis is not well understood. Here we report a genome-wide association study for a quantitative measure of relative hand skill in individuals with dyslexia [reading disability (RD)]. The most highly associated marker, rs11855415 (P = 4.7 × 10−7), is located within PCSK6. Two independent cohorts with RD show the same trend, with the minor allele conferring greater relative right-hand skill. Meta-analysis of all three RD samples is genome-wide significant (n = 744, P = 2.0 × 10−8). Conversely, in the general population (n = 2666), we observe a trend towards reduced laterality of hand skill for the minor allele (P = 0.0020). These results provide molecular evidence that cerebral asymmetry and dyslexia are linked. Furthermore, PCSK6 is a protease that cleaves the left–right axis determining protein NODAL. Functional studies of PCSK6 promise insights into mechanisms underlying cerebral lateralization and dyslexia.

The number of significant genetic associations with common complex traits is constantly increasing. However, most of these associations have not been understood at molecular level. One of the mechanisms mediating the effect of DNA variants on phenotypes is gene expression, which has been shown to be particularly relevant for complex traits1.

This method tests in a cellular context the effect of specific DNA sequences on gene expression. The principle is to measure the relative abundance of transcripts arising from the two alleles of a gene, analysing cells which carry one copy of the DNA sequences associated with disease (the risk variants)2,3. Therefore, the cells used for this method should meet two fundamental genotypic requirements: they have to be heterozygous both for DNA risk variants and for DNA markers, typically coding polymorphisms, which can distinguish transcripts based on their chromosomal origin (Figure 1). DNA risk variants and DNA markers do not need to have the same allele frequency but the phase (haplotypic) relationship of the genetic markers needs to be understood. It is also important to choose cell types which express the gene of interest. This protocol refers specifically to the procedure adopted to extract nucleic acids from fibroblasts but the method is equally applicable to other cells types including primary cells.

DNA and RNA are extracted from the selected cell lines and cDNA is generated. DNA and cDNA are analysed with a primer extension assay, designed to target the coding DNA markers4. The primer extension assay is carried out using the MassARRAY (Sequenom)5 platform according to the manufacturer's specifications. Primer extension products are then analysed by matrix-assisted laser desorption/ionization time of-flight mass spectrometry (MALDI-TOF/MS). Because the selected markers are heterozygous they will generate two peaks on the MS profiles. The area of each peak is proportional to the transcript abundance and can be measured with a function of the MassARRAY Typer software to generate an allelic ratio (allele 1: allele 2) calculation. The allelic ratio obtained for cDNA is normalized using that measured from genomic DNA, where the allelic ratio is expected to be 1:1 to correct for technical artifacts. Markers with a normalised allelic ratio significantly different to 1 indicate that the amount of transcript generated from the two chromosomes in the same cell is different, suggesting that the DNA variants associated with the phenotype have an effect on gene expression. Experimental controls should be used to confirm the results.

Background

Six independent studies have identified linkage to chromosome 18 for developmental dyslexia or general reading ability. Until now, no candidate genes have been identified to explain this linkage. Here, we set out to identify the gene(s) conferring susceptibility by a two stage strategy of linkage and association analysis.

Methodology/Principal Findings

Linkage analysis: 264 UK families and 155 US families each containing at least one child diagnosed with dyslexia were genotyped with a dense set of microsatellite markers on chromosome 18. Association analysis: Using a discovery sample of 187 UK families, nearly 3000 SNPs were genotyped across the chromosome 18 dyslexia susceptibility candidate region. Following association analysis, the top ranking SNPs were then genotyped in the remaining samples. The linkage analysis revealed a broad signal that spans approximately 40 Mb from 18p11.2 to 18q12.2. Following the association analysis and subsequent replication attempts, we observed consistent association with the same SNPs in three genes; melanocortin 5 receptor (MC5R), dymeclin (DYM) and neural precursor cell expressed, developmentally down-regulated 4-like (NEDD4L).

Conclusions

Along with already published biological evidence, MC5R, DYM and NEDD4L make attractive candidates for dyslexia susceptibility genes. However, further replication and functional studies are still required.

Background

Autism spectrum disorders (ASDs) are characterized by social, communication, and behavioral deficits and complex genetic etiology. A recent study of 517 ASD families implicated DOCK4 by single nucleotide polymorphism (SNP) association and a microdeletion in an affected sibling pair.

Methods

The DOCK4 microdeletion on 7q31.1 was further characterized in this family using QuantiSNP analysis of 1M SNP array data and reverse transcription polymerase chain reaction. Extended family members were tested by polymerase chain reaction amplification of junction fragments. DOCK4 dosage was measured in additional samples using SNP arrays. Since QuantiSNP analysis identified a novel CNTNAP5 microdeletion in the same affected sibling pair, this gene was sequenced in 143 additional ASD families. Further polymerase chain reaction-restriction fragment length polymorphism analysis included 380 ASD cases and suitable control subjects.

Results

The maternally inherited microdeletion encompassed chr7:110,663,978-111,257,682 and led to a DOCK4-IMMP2L fusion transcript. It was also detected in five extended family members with no ASD. However, six of nine individuals with this microdeletion had poor reading ability, which prompted us to screen 606 other dyslexia cases. This led to the identification of a second DOCK4 microdeletion co-segregating with dyslexia. Assessment of genomic background in the original ASD family detected a paternal 2q14.3 microdeletion disrupting CNTNAP5 that was also transmitted to both affected siblings. Analysis of other ASD cohorts revealed four additional rare missense changes in CNTNAP5. No exonic deletions of DOCK4 or CNTNAP5 were seen in 2091 control subjects.

Conclusions

This study highlights two new risk factors for ASD and dyslexia and demonstrates the importance of performing a high-resolution assessment of genomic background, even after detection of a rare and likely damaging microdeletion using a targeted approach.

Specific language impairment (SLI) is a common developmental disorder characterized by difficulties in language acquisition despite otherwise normal development and in the absence of any obvious explanatory factors. We performed a high-density screen of SLI1, a region of chromosome 16q that shows highly significant and consistent linkage to nonword repetition, a measure of phonological short-term memory that is commonly impaired in SLI. Using two independent language-impaired samples, one family-based (211 families) and another selected from a population cohort on the basis of extreme language measures (490 cases), we detected association to two genes in the SLI1 region: that encoding c-maf-inducing protein (CMIP, minP = 5.5 × 10−7 at rs6564903) and that encoding calcium-transporting ATPase, type2C, member2 (ATP2C2, minP = 2.0 × 10−5 at rs11860694). Regression modeling indicated that each of these loci exerts an independent effect upon nonword repetition ability. Despite the consistent findings in language-impaired samples, investigation in a large unselected cohort (n = 3612) did not detect association. We therefore propose that variants in CMIP and ATP2C2 act to modulate phonological short-term memory primarily in the context of language impairment. As such, this investigation supports the hypothesis that some causes of language impairment are distinct from factors that influence normal language variation. This work therefore implicates CMIP and ATP2C2 in the etiology of SLI and provides molecular evidence for the importance of phonological short-term memory in language acquisition.

Specific language impairment (SLI) is a common developmental disorder characterized by difficulties in language acquisition despite otherwise normal development and in the absence of any obvious explanatory factors. We performed a high-density screen of SLI1, a region of chromosome 16q that shows highly significant and consistent linkage to nonword repetition, a measure of phonological short-term memory that is commonly impaired in SLI. Using two independent language-impaired samples, one family-based (211 families) and another selected from a population cohort on the basis of extreme language measures (490 cases), we detected association to two genes in the SLI1 region: that encoding c-maf-inducing protein (CMIP, minP = 5.5 × 10−7 at rs6564903) and that encoding calcium-transporting ATPase, type2C, member2 (ATP2C2, minP = 2.0 × 10−5 at rs11860694). Regression modeling indicated that each of these loci exerts an independent effect upon nonword repetition ability. Despite the consistent findings in language-impaired samples, investigation in a large unselected cohort (n = 3612) did not detect association. We therefore propose that variants in CMIP and ATP2C2 act to modulate phonological short-term memory primarily in the context of language impairment. As such, this investigation supports the hypothesis that some causes of language impairment are distinct from factors that influence normal language variation. This work therefore implicates CMIP and ATP2C2 in the etiology of SLI and provides molecular evidence for the importance of phonological short-term memory in language acquisition.

Numerous genetic association studies have implicated the KIAA0319 gene on human chromosome 6p22 in dyslexia susceptibility. The causative variant(s) remains unknown but may modulate gene expression, given that (1) a dyslexia-associated haplotype has been implicated in the reduced expression of KIAA0319, and (2) the strongest association has been found for the region spanning exon 1 of KIAA0319. Here, we test the hypothesis that variant(s) responsible for reduced KIAA0319 expression resides on the risk haplotype close to the gene's transcription start site. We identified seven single-nucleotide polymorphisms on the risk haplotype immediately upstream of KIAA0319 and determined that three of these are strongly associated with multiple reading-related traits. Using luciferase-expressing constructs containing the KIAA0319 upstream region, we characterized the minimal promoter and additional putative transcriptional regulator regions. This revealed that the minor allele of rs9461045, which shows the strongest association with dyslexia in our sample (max p-value = 0.0001), confers reduced luciferase expression in both neuronal and non-neuronal cell lines. Additionally, we found that the presence of this rs9461045 dyslexia-associated allele creates a nuclear protein-binding site, likely for the transcriptional silencer OCT-1. Knocking down OCT-1 expression in the neuronal cell line SHSY5Y using an siRNA restores KIAA0319 expression from the risk haplotype to nearly that seen from the non-risk haplotype. Our study thus pinpoints a common variant as altering the function of a dyslexia candidate gene and provides an illustrative example of the strategic approach needed to dissect the molecular basis of complex genetic traits.

Author Summary

Dyslexia, or reading disability, is a common disorder caused by both genetic and environmental factors. Genetic studies have implicated a number of genes as candidates for playing a role in dyslexia. We functionally characterized one such gene (KIAA0319) to identify variant(s) that might affect gene expression and contribute to the disorder. We discovered a variant residing outside of the protein-coding region of KIAA0319 that reduces expression of the gene. This variant creates a binding site for the transcription factor OCT-1. Previous studies have shown that OCT-1 binding to a specific DNA sequence upstream of a gene can reduce the expression of that gene. In this case, reduced KIAA0319 expression could lead to improper development of regions of the brain involved in reading ability. This is the first study to identify a functional variant implicated in dyslexia. More broadly, our study illustrates the steps that can be utilized for identifying mutations causing other complex genetic disorders.

The caste system has persisted in Indian Hindu society for around 3,500 years. Like the Y chromosome, caste is defined at birth, and males cannot change their caste. In order to investigate the genetic consequences of this system, we have analysed male-lineage variation in a sample of 227 Indian men of known caste, 141 from the Jaunpur district of Uttar Pradesh and 86 from the rest of India. We typed 131 Y-chromosomal binary markers and 16 microsatellites. We find striking evidence for male substructure: in particular, Brahmins and Kshatriyas (but not other castes) from Jaunpur each show low diversity and the predominance of a single distinct cluster of haplotypes. These findings confirm the genetic isolation and drift within the Jaunpur upper castes, which are likely to result from founder effects and social factors. In the other castes, there may be either larger effective population sizes, or less strict isolation, or both.

We have established the use of a primer extension/mass spectrometry method (the PinPoint assay) for high-throughput SNP genotyping of the human Y chromosome. 118 markers were used to define 116 haplogroups and typing was organised in a hierarchical fashion. Twenty multiplex PCR/primer extension reactions were set up and each sample could be assigned to a haplogroup with only two to five of these multiplex analyses. A single aliquot of one enzyme was found to be sufficient for both PCR and primer extension. We observed 100% accuracy in blind validation tests. The technique thus provides a reliable, cost-effective and automated method for Y genotyping, and the advantages of using a hierarchical strategy can be applied to any DNA segment lacking recombination.