The Dyslexia Candidate Locus on 2p12 Is Associated with General Cognitive Ability and White Matter Structure

Abstract

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.

About the Authors:

Thomas S. Scerri

Current address: The Walter and Eliza Hall Institute for Medical Research, Melbourne, Australia

Affiliation: Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom

Fahimeh Darki

Affiliation: Department of Neuroscience, Karolinska Institutet, Stockholm, Sweden

Dianne F. Newbury

Affiliation: Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom

Andrew J. O. Whitehouse

Affiliation: Telethon Institute for Child Health Research, Centre for Child Health Research, University of Western Australia, West Perth, Australia

Myriam Peyrard-Janvid

Affiliation: Department of Biosciences and Nutrition at Novum, Karolinska Institutet, Stockholm, Sweden

Hans Matsson

Affiliation: Department of Biosciences and Nutrition at Novum, Karolinska Institutet, Stockholm, Sweden

Qi W. Ang

Affiliation: School of Women's and Infants' Health, University of Western Australia, Crawley, Australia

Craig E. Pennell

Affiliation: School of Women's and Infants' Health, University of Western Australia, Crawley, Australia

Susan Ring

Affiliation: School of Social and Community Medicine, University of Bristol, Bristol, United Kingdom

John Stein

Affiliation: Department of Physiology, University of Oxford, Oxford, United Kingdom

Andrew P. Morris

Affiliation: Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom

Anthony P. Monaco

Affiliation: Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom

Juha Kere

Affiliations Department of Biosciences and Nutrition at Novum, Karolinska Institutet, Stockholm, Sweden, Molecular Medicine Research Programme, University of Helsinki, and Folkhälsan Institute of Genetics, Helsinki, Finland, Science for Life Laboratory, Karolinska Institutet, Stockholm, Sweden

Joel B. Talcott

Affiliation: Aston Brain Centre, School of Life and Health Sciences, Aston University, Birmingham, United Kingdom

Torkel Klingberg

Affiliation: Department of Neuroscience, Karolinska Institutet, Stockholm, Sweden

Silvia Paracchini

* E-mail: [email protected]

Affiliations Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom, School of Medicine, University of St Andrews, St Andrews, United Kingdom

Introduction

Dyslexia (or reading disability, RD) and specific language impairment (SLI) are common neurodevelopmental disorders, reflecting specific deficits in the acquisition of literacy skills and oral language, respectively [1]. For both disorders, a diagnosis is achieved by excluding known causes of the deficits, such as co-occurring sensory or neurological impairment or lack of educational opportunity. RD and SLI are complex traits resulting from the interaction of multiple factors of both genetic and environmental origin, yet the biological underpinnings and cascading cognitive deficits remain poorly understood. Several genes have been proposed as susceptibility candidates for both RD and SLI. The RD candidates include the ROBO1, KIAA0319, DCDC2 and DYX1C1 genes and the MRPL19/C2ORF3 locus [2]. The SLI candidates include the CMIP, ATP2C2 and CNTNAP2 genes [3]. It has been shown that most of the dyslexia risk genes (KIAA0319, DCDC2, DYX1C1 and ROBO1) are involved in cortical development and specifically in neuronal migration [4]. An important unanswered research question is how genes involved in such a general process as cortical development contribute to the risk for specific neurodevelopmental disorders. Several studies have addressed this research question by investigating whether candidate genes implicated in a specific disorder show pleiotropic effects, which could, for example, help explain the high comorbidity observed between SLI, RD and ADHD [5], [6], [7]. To date only modest associations have been reported for shared genetic factors across these separate conditions [8], [9], [10]. Recently we conducted an association study in the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort to explore potential genetic overlaps between SLI and RD [11]. We analysed a panel of single nucleotide polymorphisms (SNPs), previously reported to be associated with either RD or SLI, and tested them for association with literacy and language related phenotypes. We reported highly specific effects for DCDC2, KIAA0319 and CMIP with measures of single-word reading and spelling [11]. These data suggested genetic effects on specific and independent aspects of cognitive function rather than on multiple or more generalised phenotypes.

A direct genetic effect of the impact on broader cognitive abilities for RD and SLI candidate genes has never been tested. General cognitive ability, assessed with standardised intelligence tests, is a highly heritable trait [12], [13], yet very few genes have been identified that impact upon these trait [14]. Most of the reported candidate genes have lacked adequate sample size and replications [15]. A recent genome-wide association study (GWAS) supports the strong heritability of cognitive abilities, but with effects that are spread across a large number of genetic factors and therefore not easily detectable in isolation [16]. In the present study, we used a candidate gene approach to investigate a precise question, namely whether genetic risk factors for RD and SLI have a broader impact on general cognitive function. This hypothesis is supported by the consistent observation of significant correlation between reading abilities and general cognition. We analysed RD and SLI candidate genes for association with general cognitive ability in the ALSPAC cohort and detected a statistically significant association at the chromosome 2p12 dyslexia-associated locus. The effect size was small but reproducible in independent samples. Given previous associations between white matter structure and both language-related phenotypes [17] and IQ [18], we looked for any correlations between genotypes at the MRPL19/C2ORF3 locus and white matter structure. We found associations in the posterior corpus callosum and cingulum, connecting large sections of the parietal, occipital and temporal cortices. The widespread connectivity of this white matter region is consistent with a more general effect on both language and intellectual function.

Results

We analysed 19 SNPs for association with verbal IQ (VIQ) and non-verbal IQ (performance IQ; PIQ) in the ALSPAC child cohort (Table 1) across the MRPL19/C2ORF3, KIAA0319, DCDC2, ATP2C2 and CMIP genes (Table 2). These markers were selected for previously reported associations with reading- and language-related phenotypes [11]. We detected significant associations (p-values<0.001) with VIQ and the SNPs rs714939 (MRPL12/C2ORF3 locus) and rs6935076 (KIAA0319). The rs714939 SNP was also the only marker that showed a trend for association with PIQ (P = 0.006). A marker in CMIP, rs6564903, also showed a trend of association for VIQ (P = 0.004). Our previous study [11] showed that KIAA0319 and CMIP are associated with a measure of single word reading (READ) [19] in a subset of the ALSPAC sample. Given the high correlation between VIQ and READ (r = 0.571; Table 1) [11], it is plausible that the association we observed here was driven by reading ability. Therefore, we included READ as covariate when analysing the association with VIQ (Table 2). This analysis showed that the associations with rs6935076 (KIAA0319) and rs6564903 (CMIP) became attenuated, whereas rs714939 remained significantly associated with VIQ. Therefore, while the association with the KIAA0319 and CMIP markers may have been driven by the correlation with reading, the associations at the 2p12 locus appeared to be specific to IQ.

[Figure omitted. See PDF.]

Table 1. Details of phenotypic measures.

https://doi.org/10.1371/journal.pone.0050321.t001

[Figure omitted. See PDF.]

Table 2. Single-SNP analyses in the ALSPAC cohort.

https://doi.org/10.1371/journal.pone.0050321.t002

Another SNP at the MRPL19/C2ORF3 locus, rs917235, also showed a trend of association with VIQ (Table 2). For both rs917235 and rs714939, the G alleles are associated with lower performance. These were the same alleles associated with dyslexia in the original report [20]. The most robust association in this previous study was reported for two haplotypes, in different populations, that overlapped for rs917235 and rs714939; specifically, the “GG” haplotype formed by these two markers increased susceptibility to dyslexia. Accordingly, we also tested the rs917235/rs714939 haplotypes for association with VIQ (Table 3) in the ALSPAC cohort. The “GG” haplotype gave the strongest association with lower VIQ scores (P = 0.00009). When READ was used as a covariate, the association became weaker but remained significant (P = 0.00024). Furthermore, the “AA” haplotype showed significant association with higher VIQ scores and with an effect size of similar magnitude.

[Figure omitted. See PDF.]

Table 3. Summary results for rs917235/rs714939 haplotype associations with VIQ in the ALSPAC cohort.

https://doi.org/10.1371/journal.pone.0050321.t003

To replicate these results, we tested the MRPL19/C2ORF3 markers for association with VIQ and PIQ in four independent samples (Table 4). We consistently found a trend of association across the four collections, where either of the two originally associated markers had the G allele associated with lower performance. Associations were particularly evident in the SLI families (rs714939; P = 0.009; PIQ) and in the unrelated cases with dyslexia (rs917235; P = 0.0004; VIQ). The same SNPs did not show association with single word reading measures in previous studies using the same samples [9], [11], [21] excluding the possibility that these associations are driven by the high correlation between VIQ and reading ability. We also re-analysed all of the samples for association with reading using IQ as covariate, but again none of the samples including the ALSPAC cohort showed any associations (data not shown). The unrelated cases with dyslexia was the only sample that showed stronger association with VIQ, while the other samples showed associations mainly with PIQ. The dyslexia families and Raine samples only showed a modest trend of association, but it was in the same allelic direction (rs917235; P = 0.064 and rs714939; P = 0.059, respectively). Haplotype analysis in all replication sets revealed a similar trend, with the “GG” haplotype associated with lower performance, but did not show associations stronger than single SNP analyses (data not shown).

[Figure omitted. See PDF.]

Table 4. P-values from single-SNP analyses for VIQ and PIQ in replication samples.

https://doi.org/10.1371/journal.pone.0050321.t004

We assessed the effect of genetic variants at the MRPL19/C2ORF3 locus on white matter structure in the Swedish cohort. Among seven genotyped markers, only rs917235 showed significant association with variation in white matter volume (Pcorrected = 1.27×10−3 at cluster level with the threshold of P<0.01). Specifically, the G allele, already shown to be associated with lower IQ, was associated with lower white matter volume. This white matter cluster was located bilaterally and confined mainly to the posterior part of the corpus callosum and the cingulum (Figure 1A).

[Figure omitted. See PDF.]

Figure 1. Imaging results.

(A) Main association of rs917235 with white matter volume. Sagittal view of the significant cluster associated with rs917235 at the MNI coordinates of X = −13, −8, 8, and 13 mm relative to the midline. The background image is the average of all individuals' white matter segmented images. The colour-bar shows the z-scores of the statistical analysis. (B–D) Structural connectivity. (B) A sample of fiber tracking from one individual passing through the region associated with rs917235. Different colours show the direction of pathways (red for left-right, blue for superior-inferior, and green for anterior-posterior). (C) Overlay of tract tracing from 30 randomly selected individuals. The colour-bar is a count for the number of subjects; yellow shows the most probable pathways. (D) Cortical end-points of the white matter pathways showing that the fibers connect right postcentral gyrus, superior parietal lobule, precuneous, occipital cortex and temporal fusiform gyrus to the analogous left regions.

https://doi.org/10.1371/journal.pone.0050321.g001

The significant cluster in the posterior corpus callosum (Figure 1A) was then used as a seed region for fiber tracking (see Figure 1B for tracking in one individual) performed in 30 randomly selected subjects. The pathways from each subject were then overlaid on a common template to identify the most consistent localization of pathways (Figure 1C). Next, we identified the parts of cortex that were closest to the end-points of the white matter tracts. This analysis showed that the pathways passing through the region connected the right postcentral gyrus, superior parietal lobule, precuneous, lateral occipital cortex and fusiform gyrus to analogous areas in the left hemisphere (Figure 1D).

Discussion

We evaluated the effect of candidate genes for RD and SLI on measures of general cognitive ability in the ALSPAC cohort. The MRPL19/C2ORF3 locus showed a statistically significant association with general cognitive abilities (VIQ and PIQ). KIAA0319 and CMIP variants also showed associations, but these genes have previously been found associated with reading ability in the ALSPAC sample [11], hence the associations detected here could be an artefact due to the correlation between VIQ and measures of single word-reading. Indeed, the associations between the KIAA0319 and CMIP markers and VIQ were attenuated when a reading measure was used as a covariate. Conversely, MRPL19/C2ORF3 markers did not show any association with single-word reading in the ALSPAC cohort [11]. Furthermore, this association remained significant when reading ability was used as a covariate and is supported by a consistent trend of association with IQ in four independent samples. The association would not stand a genome-wide correction for significance level, but this would be consistent with a recent GWAS that predicts multiple and small-size factors contributing to IQ [16]. Our candidate gene approach was based on the selection of previously reported disorder-associated genotypes, suggesting a more generalised effect across multiple cognitive traits. Our interpretation is that the chromosome 2 locus has an effect on cognition by influencing neurodevelopment. The effect on specific endophenotypes will depend on ascertainment criteria, methodology for phenotypic assessment and differences in genetic background. All these elements will have a significant role in the outcome of association analyses when sample sizes are relatively small.

To date, only a few common genetic variants have been associated with measures of verbal and nonverbal reasoning. These include variants in the CHRM2, COMT and BDNF genes [14]. GWASs, which have recently had an enormous success in identifying genetic variants contributing to complex traits, had little success in mapping variants associated with cognitive abilities, with no indication of major genetic contributing loci [16], [22], [23]. The generally accepted model proposes that cognitive abilities are influenced by many genes of small effect [24] and are therefore difficult to map in the relatively small-sized samples currently available. Gene group-base analysis might be more effective by testing biological functions determined by multiple genes [25].

Both rs714939 and rs917235 are within an intergenic region between the FAM176A and MRPL19/C2ORF3 loci. The original study associating this locus with dyslexia showed an effect on gene expression of the co-regulated genes MRPL19 and C2ORF3, and differential expression across a set of brain regions [20]. These genes are transcribed in a head-to-head orientation on chromosome 2p12. Visualisation of MRPL19 and C2ORF3 in the Allen Brain Atlas (a resource of human gene expression data derived from 3 post-mortem males aged 24 to 57 years old at time of death; http://www.brain-map.org/), show that these two genes are most highly expressed in white matter of the corpus callosum and the cingulum when compared to all other regions. This is in contrast to KIAA0319, CMIP and ATP2C2 which show very low expression in white matter, and FAM176A which shows moderate expression. There is no direct evidence that these two SNPs are the causative variants and it is likely that together they tag the effect of a nearby functional factor. The FAM176A gene, also expressed in fetal brain, cannot be ruled out as being influenced by these two SNPs. Interestingly, rs714939 is located in a region of high H3K4Me1 marks (ENCODE ChIP-Seq data, GRCh37/hg19 assembly visualised in the UCSC Genome Browser (http://genome.ucsc.edu/cgi-bin/hgGateway)). H3K4Me1 is a monomethylation of lysin 4 of the H3 histone protein and this modification is often found near regulatory elements. This region of high H3K4Me1 spans 9–10 kb and could be involved in the cis-regulation of any of the three neighbouring genes.

The imaging results provide a neural correlate to the genetic polymorphism. Our interpretation is that the genetic variants contribute to structural variation in a relatively large area of white matter which affects both reasoning and reading. It should be pointed out, however, that this is an association study that does not provide any direction of causality. It cannot be excluded that genes affect behavior, which in turn affect the white matter. Our neuroimaging analysis revealed a significant association between rs917235 and white matter volume of the posterior part of the corpus callosum and cingulum. The G allele, associated with lower IQ in the behavioural analysis, was specifically associated with lower white matter volume in this region. Previous studies have reported correlations between white matter structure and measures of IQ [26], [27], [28], [29], supporting the idea that neural connectivity of the brain is an endophenotype underpinning general intellectual ability [18]. Tract tracing of fibers passing through the correlated white matter region showed that they establish long-range connections with large regions of the parietal and occipital cortices, and a smaller region within temporo-occipital cortex. Intra- and superior parietal have been associated with performance on reasoning tasks [30] and the inferior parietal cortex is important for language-related functions [31]. The thickness of the splenium of the corpus callosum, which connects large parts of the occipital, parietal and temporal lobes, has been previously associated with intelligence [32]. Based on parieto-frontal integration theory of intelligence (P-FIT) [30], the lateral cortex and fusiform gyrus from the temporal lobe are involved in cognitive ability since they participate in visual perception, recognition and imagination. The superior parietal, supramarginal and angular gyri are involved in structural symbolism and abstraction. These regions are also important for language-related functions. The posterior part of the corpus callosum has been previously reported as a white matter region with structural differences between normal and dyslexic readers. [33], [34], [35], [36]. It has been also shown that the posterior part of the corpus callosum is bigger in children with dyslexia rather than in typically developing children [37] The widespread connectivity of the white matter region associated with rs917235 is thus consistent with the previous neuroimaging associations to both language and general cognitive abilities. However, substantially more cortical regions and associated brain functions are likely to underlie general cognitive abilities. The integration of genetics with structural and functional imaging approaches holds potential for further elucidating the effects of genes on the normal and atypical development of cognitive function. With a similar approach we showed that three dyslexia susceptibility genes DYX1C1, DCDC2 and KIAA0319 are associated with white matter volume in distinct but overlapping regions of the left temporo-parietal hemisphere, and that the white matter volume in these regions also correlated with reading ability [38]. A recent study found association between candidate genes for language and reading impairment (FOXP2 and KIAA0319, respectively) and regionally specific brain activations assessed with fMRI [39]. Another study reported two independent SNP associations, both on chromosome 12, with brain volume measures as well as suggestive evidence of an effect on cognitive abilities [40]. Most recently, a study has found a genome-wide significant association between a SNP (rs2298948) within C2ORF3 (called GCFC2 in that study) and hippocampal volume [41]. The hippocampus plays an important role in learning and spatial memory, and this is correlated with hippocampus volume.

The association of cognitive abilities with a genetic locus, originally identified as a candidate for dyslexia susceptibility, is independent from a clinical diagnosis of dyslexia and SLI. The association was statistically significant in the ALSPAC sample, which represents the general population, and we did not observe any associations with reading-related measures in either ALSPAC or the other replication samples, including the ones selected for dyslexia or SLI. It is possible that the original association with developmental dyslexia at the MRPL19/C2ORF3 locus [20] was due to a sampling effect and the high correlation between the verbal component of cognitive abilities and reading. Alternatively, the same genetic variants may have different phenotypic effects when combined with alternative genetic or environmental factors and would become apparent in separate sample subgroups. Multivariate genetic analysis have consistently suggested a correlation of about 0.6 between general cognitive ability and learning [42] but there is less agreement in estimating the effects across the range of the entire phenotypic distribution [43]. This would imply that some factors may have specific effects only at the extremes of the phenotype where disorder diagnosis would apply.

In summary, we report an association of measures of general cognitive abilities with the chromosome 2p12 locus implicated in dyslexia. We show, for the first time, that the same genetic locus is associated with white matter volume in the posterior corpus callosum. Furthermore, fibers throughout this region connected cortical regions involved in both language and general cognitive abilities. Follow-up studies might identify the functional genetic variant(s) and the gene(s) implicated. Such findings will contribute to our understanding of the biological pathways underlying normal and atypical cognition and the possible shared factor(s) mediating general cognitive functions and highly prevalent developmental disorders such as dyslexia.

Materials and Methods

Ethics Statement

Informed written consent was obtained from the parents, with the option for them or their children to withdraw at any time. Ethical approval for the ALSPAC cohort was obtained from the ALSPAC Law and Ethics Committee and the Local Research Ethics Committees. Ethical approval for the SLIC cohort was granted by local ethics committees. Ethical approval for the Oxford/Reading and Aston studies was acquired from the Oxfordshire Psychiatric Research Ethics Committee (OPREC O01.02). For the Raine Study, participant recruitment and all follow-ups of the study families were approved by the Human Ethics Committee at King Edward Memorial Hospital and/or Princess Margaret Hospital for Children in Perth. For the Swedish sample, the study was approved by the Ethics Board of the Karolinska University Hospital (Stockholm).

Initial Sample

The ALSPAC cohort consists of over 14,000 children from the southwest of England that had expected dates of delivery between 1st April 1991 and 31st December 1992 [44]. From age 7 years, all children were annually assessed for a wide range of physical, behavioural, and neuropsychological traits, including reading and language-related measures. DNA is available for approximately 11,000 children. For this study, individuals with a non-white ancestry were excluded and after filtering for missing genotypic or phenotypic data we conducted the analysis in a sample of 5905 children. Cognitive ability was assessed using the Wechsler Intelligence Scales for Children (WISC-III) [45] for both verbal and performance IQ (VIQ and PIQ, respectively; Table 1). The VIQ scale included the subtests for Information, Similarities, Arithmetic, Vocabulary, Comprehension and Digit span. The PIQ scale included the subtests for Picture completion, Coding, Picture arrangement, Block design and Object assembly.

Replication Samples

The SLI Consortium (SLIC) cohort has been described in detail previously [46], [47], [48]. This family-based sample includes approximately 400 individuals from 181 families. The samples were assessed at one of five separate centres across the UK: The Newcomen Centre at Guy's Hospital, London, the Cambridge Language and Speech Project (CLASP [49]), the Child Life and Health Department at the University of Edinburgh [50], the Department of Child Health at the University of Aberdeen and the Manchester Language Study [51], [52]. Cognitive ability of all children in the SLIC sample was assessed using the WISC-III [45] applying the same VIQ and PIQ subtests listed for ALSPAC. In the SLIC sample, the PIQ score (cut-off at PIQ>80) was used to exclude children whose language problems were accompanied by deficits in non-verbal skills.

The dyslexia-based sample has been described previously [53], [54]. It includes 684 siblings from 288 unrelated nuclear families and 282 unrelated cases with dyslexia, recruited through the Dyslexia Research Centre clinics in Oxford and Reading, and the Aston Dyslexia and Development Assessment Centre in Birmingham. VIQ and PIQ were obtained from the BAS similarities and BAS matrices subtests respectively [55]. The similarities sub-scale of the Wechsler Adult Intelligence Scales (WAIS), a measure analogous to the BAS similarities test, was used when age was >17.5 years [45].

The Western Australian Pregnancy Cohort (Raine) study is a longitudinal investigation of 2900 pregnant women and their offspring recruited between 1989 and 1991 [56]. From the original cohort, 2868 children have been followed over two decades. The Raine sample is representative of the larger Australian population (88% Caucasian); only those children with both biological parents of white European origin were included in the current analyses. Verbal and non-verbal ability was assessed at 10 years of age using the Peabody Picture Vocabulary Test–Revised (PPVT-R) [57] and the Raven's Coloured Progressive Matrices (RCPM) [58], respectively. The PPVT-R provides a measure of receptive vocabulary, requiring children to select which of four pictures corresponds to an aurally presented word. Raw scores are converted to a Verbal IQ, standardized for age 2 years and above (based around a mean of 100 and a SD of 15). RCPM is a 36 item multiple choice test that presents a matrix-like arrangement of figural symbols and requires the child to select the missing symbols from a set of six alternatives. Raw scores are converted to percentiles, which provide an indication of performance relative to other children of a similar age. This assessment is standardized for children between 4.9 and 12.0 years of age.

The Swedish sample used for neuroimaging consists of 76 Swedish speaking children and young adults (age range 6 to 25 years) randomly selected from the “Brainchild” study, a longitudinal study of typical development [59], [60]. The participants were from the population register in the town of Nynäshamn, Sweden, and showed no evidence of neurological or psychological disorders. DNA was available for all subjects.

Genotyping and statistical analysis

We analysed 19 SNPs across the MRPL19/C2ORF3, KIAA0319, DCDC2, ATP2C2 and CMIP loci that were recently genotyped in the ALSPAC child cohort [11]. The sample size in the present study is larger because we did not apply an IQ filter as described in the previous analysis. SNPs were genotyped using either Sequenom iPLEX assays according to the manufacturer's instructions or the KBiosciences service using their in-house technology (http://www.kbioscience.co.uk/). The genotyping in the samples selected for dyslexia and SLI was conducted with Sequenom iPLEX assays as part of previous studies [9].

For the Raine study, DNA samples were genotyped on an Illumina 660 Quad Array [61].

For the Swedish sample, seven SNPs (rs3088180, rs4853169, rs917235, rs6732511, rs714939, rs17689640 and rs17689863) located at the MRPL19/C2ORF3 locus were genotyped using the Sequenom iPLEX system as previously described [59].

Quantitative analyses were performed within PLINK (1.07) [62] using additive tests of association. We included two additional phenotypes (VIQ and PIQ) to the multiple testing correction applied in our previous study of the ALSPAC sample [11] which corrected for the analysis of 11 clusters of SNPs showing significant linkage disequilibrium (LD; r2>0.6) and 2 phenotypes. Therefore we corrected here the significance level of P = 0.05 for 44 independent tests (11 SNP clusters analysed for four phenotypes) resulting in P = 0.001. The ALSPAC cohort has been tested previously for other SNPs and phenotypes, therefore these tests should be considered in calculating a significant threshold p-value, or the genome-significant threshold of 5×10−8 should be applied. However, this would be too conservative for the scope of this study which analyses previously reported genetic markers and tests a specific hypothesis rather than conducting an explorative exercise. Family-based cohorts were analysed using QTDT [63].

Image analysis

For the Swedish sample, three-dimensional structural T1-weighted imaging with magnetization-prepared rapid gradient echo sequence (TR = 2300 ms, TE = 2.92 ms, 256×256 mm, 176 sagittal slices and 1 mm3 voxel size) with the field of view of 256×256 mm, 256 slices, and 1 mm3 voxel size was carried out on all the participants and repeated two years later for 69 subjects. White matter segmentation, followed by an alignment technique, was performed on the structural data using the Diffeomorphic Anatomical Registration using Exponentiated Lie algebra (DARTEL) toolbox in SPM (www.fil.ion.ucl.ac.uk/spm/software/spm5). Images were then spatially smoothed with an 8 mm Gaussian kernel. The DARTEL outputs are white matter segmented images which reflect the signal intensity modulated by volume transformations applied to individual images to register them into the MNI template.

Diffusion tensor imaging (DTI) was also acquired using a field of view 230×230 mm, matrix size 128×128 mm, 19 slices with 6.5 mm thickness b-max 1000 s/mm2 in 20 directions. Eddy current and head motion were corrected using affine registration to a reference volume using FSL (www.fmrib.ox.ac.uk/fsl/). The diffusion tensors were then computed for each voxel and the DTI and fractional anisotropy (FA) data were then constructed.

The seven SNPs were entered separately as a main factor in a flexible factorial design second-level SPM analysis. This included both the individual images with and without repeated measures, to assess the variation of white matter with respect to the genetic markers and was corrected for the effect of age, sex, handedness and total white matter volume. Age X gene and gender X gene interaction effects were also included in the model. As a part of this exploratory analysis, the significance level was corrected at the cluster level using non-stationary cluster extent correction [64] for multiple comparisons resulting from searching the entire white matter volume as well as for the seven SNPs (Bonferroni correction, Pcorrected<0.0014 for comparison of searching the entire brain with a threshold of P<0.01, as implemented in the SPM software).

The region showing the significant effect was saved as a binary region of interest (ROI). This ROI was then transformed to each individual's DTI space, to be used as a seed ROI for white matter fiber tracking. Deterministic fiber tracking was applied on 30 randomly selected subjects by starting tractography from the ROI following the principal eigenvector direction using 1 mm steps, considering thresholds of 0.15 for FA values and 30 for angular degree using ExploreDTI v4.7.3. (www.exploredti.com). Computed tracts from all individuals were then averaged across all 30 participants to derive a probabilistic map of the white matter pathways passing through the overlapping areas.

Acknowledgments

We are grateful to all the families and sample donors who took part in this study. We are grateful to Dorothy Bishop for her guidance in handling cognitive phenotypes and useful discussion. SP is a Royal Society University Research Fellow. DFN is an MRC Career Development Fellow and a Junior Research Fellow of St John's College, Oxford. AW is a NHMRC Career Development Fellow. This publication is the work of the authors and SP will serve as guarantor for the contents of this paper.

Author Contributions

Conceived and designed the experiments: CEP SR JS A.P. Monaco JK JBT TK SP. Performed the experiments: TSS FD DN AJOW MP HM SP. Analyzed the data: TSS FD DN AJOW MP HM QWA A.P. Morris JBT SP. Wrote the paper: TSS FD DN AJOW MP HM JK JBT SP.

Citation: Scerri TS, Darki F, Newbury DF, Whitehouse AJO, Peyrard-Janvid M, Matsson H, et al. (2012) The Dyslexia Candidate Locus on 2p12 Is Associated with General Cognitive Ability and White Matter Structure. PLoS ONE 7(11): e50321. https://doi.org/10.1371/journal.pone.0050321

References

1. Pennington BF, Bishop DV (2009) Relations among speech, language, and reading disorders. Annu Rev Psychol 60: 283–306.

2. Scerri TS, Schulte-Korne G (2010) Genetics of developmental dyslexia. Eur Child Adolesc Psychiatry 19: 179–197.

3. Newbury DF, Monaco AP (2010) Genetic advances in the study of speech and language disorders. Neuron 68: 309–320.

4. Galaburda AM, LoTurco J, Ramus F, Fitch RH, Rosen GD (2006) From genes to behavior in developmental dyslexia. Nat Neurosci 9: 1213–1217.

5. Snowling M, Bishop DV, Stothard SE (2000) Is preschool language impairment a risk factor for dyslexia in adolescence? J Child Psychol Psychiatry 41: 587–600.

6. McArthur GM, Hogben JH, Edwards VT, Heath SM, Mengler ED (2000) On the “specifics” of specific reading disability and specific language impairment. J Child Psychol Psychiatry 41: 869–874.

7. Pennington BF (2006) From single to multiple deficit models of developmental disorders. Cognition 101: 385–413.

8. Rice ML, Smith SD, Gayan J (2009) Convergent genetic linkage and associations to language, speech and reading measures in families of probands with Specific Language Impairment. J Neurodev Disord 1: 264–282.

9. Newbury DF, Paracchini S, Scerri TS, Winchester L, Addis L, et al. (2011) Investigation of dyslexia and SLI risk variants in reading and language-impaired subjects. Behav Genet 41: 90–104.

10. Bates TC, Luciano M, Medland SE, Montgomery GW, Wright MJ, et al. (2011) Genetic Variance in a Component of the Language Acquisition Device: ROBO1 Polymorphisms Associated with Phonological Buffer Deficits. Behav Genet 41: 50–57.

11. Scerri TS, Morris AP, Buckingham LL, Newbury DF, Miller LL, et al. (2011) DCDC2, KIAA0319 and CMIP are associated with reading-related traits. Biol Psychiatry 70: 237–245.

12. Bouchard TJ Jr, McGue M (1981) Familial studies of intelligence: a review. Science 212: 1055–1059.

13. Plomin R, Spinath FM (2002) Genetics and general cognitive ability (g). Trends Cogn Sci 6: 169–176.

14. Deary IJ, Johnson W, Houlihan LM (2009) Genetic foundations of human intelligence. Hum Genet 126: 215–232.

15. Chabris CF, Hebert BM, Benjamin DJ, Beauchamp JP, Cesarini D, et al. (2012) Most reported genetic associations with general intelligence are probably false positives. Psychological Science In press.

16. Davies G, Tenesa A, Payton A, Yang J, Harris SE, et al. (2011) Genome-wide association studies establish that human intelligence is highly heritable and polygenic. Mol Psychiatry 16: 996–1005.

17. Klingberg T, Hedehus M, Temple E, Salz T, Gabrieli JD, et al. (2000) Microstructure of temporo-parietal white matter as a basis for reading ability: evidence from diffusion tensor magnetic resonance imaging. Neuron 25: 493–500.

18. Deary IJ, Penke L, Johnson W (2010) The neuroscience of human intelligence differences. Nat Rev Neurosci 11: 201–211.

19. Rust J, Golombok S, Trickey G (1993) WORD: Wechsler Objective Reading Dimensional Manual. Sidcup, UK: Psychological Corporation.

20. Anthoni H, Zucchelli M, Matsson H, Muller-Myhsok B, Fransson I, et al. (2007) A locus on 2p12 containing the co-regulated MRPL19 and C2ORF3 genes is associated to dyslexia. Hum Mol Genet 16: 667–677.

21. Paracchini S, Ang QW, Stanley FJ, Monaco AP, Pennell CE, et al. (2010) Analysis of dyslexia candidate genes in the Raine cohort representing the general Australian population. Genes Brain Behav 10: 158–165.

22. Davis OS, Butcher LM, Docherty SJ, Meaburn EL, Curtis CJ, et al. (2010) A three-stage genome-wide association study of general cognitive ability: hunting the small effects. Behav Genet 40: 759–767.

23. Pan Y, Wang KS, Aragam N (2011) NTM and NR3C2 polymorphisms influencing intelligence: family-based association studies. Prog Neuropsychopharmacol Biol Psychiatry 35: 154–160.

24. Plomin R (1999) Genetics and general cognitive ability. Nature 402: C25–29.

25. Ruano D, Abecasis GR, Glaser B, Lips ES, Cornelisse LN, et al. (2010) Functional gene group analysis reveals a role of synaptic heterotrimeric G proteins in cognitive ability. Am J Hum Genet 86: 113–125.

26. Deary IJ, Bastin ME, Pattie A, Clayden JD, Whalley LJ, et al. (2006) White matter integrity and cognition in childhood and old age. Neurology 66: 505–512.

27. Schmithorst VJ, Wilke M, Dardzinski BJ, Holland SK (2005) Cognitive functions correlate with white matter architecture in a normal pediatric population: a diffusion tensor MRI study. Hum Brain Mapp 26: 139–147.

28. Chiang MC, Barysheva M, Shattuck DW, Lee AD, Madsen SK, et al. (2009) Genetics of brain fiber architecture and intellectual performance. J Neurosci 29: 2212–2224.

29. Fjell AM, Westlye LT, Amlien IK, Walhovd KB (2011) Reduced white matter integrity is related to cognitive instability. J Neurosci 31: 18060–18072.

30. Jung RE, Haier RJ (2007) The Parieto-Frontal Integration Theory (P-FIT) of intelligence: converging neuroimaging evidence. Behav Brain Sci 30: 135–154 discussion 154–187.

31. Richlan F, Kronbichler M, Wimmer H (2011) Meta-analyzing brain dysfunctions in dyslexic children and adults. Neuroimage 56: 1735–1742.

32. Luders E, Thompson PM, Narr KL, Zamanyan A, Chou YY, et al. (2011) The link between callosal thickness and intelligence in healthy children and adolescents. Neuroimage 54: 1823–1830.

33. Vandermosten M, Boets B, Poelmans H, Sunaert S, Wouters J, et al. (2012) A tractography study in dyslexia: neuroanatomic correlates of orthographic, phonological and speech processing. Brain 135: 935–948.

34. Frye RE, Hasan K, Xue L, Strickland D, Malmberg B, et al. (2008) Splenium microstructure is related to two dimensions of reading skill. Neuroreport 19: 1627–1631.

35. Dougherty RF, Ben-Shachar M, Deutsch GK, Hernandez A, Fox GR, et al. (2007) Temporal-callosal pathway diffusivity predicts phonological skills in children. Proc Natl Acad Sci U S A 104: 8556–8561.

36. Odegard TN, Farris EA, Ring J, McColl R, Black J (2009) Brain connectivity in non-reading impaired children and children diagnosed with developmental dyslexia. Neuropsychologia 47: 1972–1977.

37. Hasan KM, Molfese DL, Walimuni IS, Stuebing KK, Papanicolaou AC, et al. (2012) Diffusion tensor quantification and cognitive correlates of the macrostructure and microstructure of the corpus callosum in typically developing and dyslexic children. NMR Biomed

38. Darki F, Peyrard-Janvid M, Matsson H, Kere J, Klingberg T (2012) Three Dyslexia Susceptibility Genes, DYX1C1, DCDC2, and KIAA0319, Affect Temporo-Parietal White Matter Structure. Biol Psychiatry

39. Pinel P, Fauchereau F, Moreno A, Barbot A, Lathrop M, et al. (2012) Genetic Variants of FOXP2 and KIAA0319/TTRAP/THEM2 Locus Are Associated with Altered Brain Activation in Distinct Language-Related Regions. J Neurosci 32: 817–825.

40. Stein JL, Medland SE, Vasquez AA, Hibar DP, Senstad RE, et al. (2012) Identification of common variants associated with human hippocampal and intracranial volumes. Nat Genet

41. Melville SA, Buros J, Parrado AR, Vardarajan B, Logue MW, et al. (2012) Multiple loci influencing hippocampal degeneration identified by genome scan. Ann Neurol 72: 65–75.

42. Plomin R, Kovas Y (2005) Generalist genes and learning disabilities. Psychol Bull 131: 592–617.

43. Haworth CM, Dale PS, Plomin R (2009) Generalist genes and high cognitive abilities. Behav Genet 39: 437–445.

44. Golding J, Pembrey M, Jones R (2001) ALSPAC–the Avon Longitudinal Study of Parents and Children. I. Study methodology. Paediatr Perinat Epidemiol 15: 74–87.

45. Wechsler D, Golombok S, Rust J (1992) WISC-IIIUK: Wechsler Intelligence Scale for Children. Sidcup, UK: The Psycological Corporation.

46. SLIC (2002) A genomewide scan identifies two novel loci involved in specific language impairment. Am J Hum Genet 70: 384–398.

47. SLIC (2004) Highly significant linkage to the SLI1 locus in an expanded sample of individuals affected by specific language impairment. Am J Hum Genet 74: 1225–1238.

48. Falcaro M, Pickles A, Newbury DF, Addis L, Banfield E, et al. (2008) Genetic and phenotypic effects of phonological short-term memory and grammatical morphology in specific language impairment. Genes Brain Behav 7: 393–402.

49. Burden V, Stott CM, Forge J, Goodyer I (1996) The Cambridge Language and Speech Project (CLASP). I. Detection of language difficulties at 36 to 39 months. Dev Med Child Neurol 38: 613–631.

50. Clark A, O'Hare A, Watson J, Cohen W, Cowie H, et al. (2007) Severe receptive language disorder in childhood–familial aspects and long-term outcomes: results from a Scottish study. Arch Dis Child 92: 614–619.

51. Conti-Ramsden G, Botting N (1999) Characteristics of children attending language units in England: a national study of 7-year-olds. Int J Lang Commun Disord 34: 359–366.

52. Conti-Ramsden G, Crutchley A, Botting N (1997) The extent to which psychometric tests differentiate subgroups of children with SLI. J Speech Lang Hear Res 40: 765–777.

53. Francks C, Paracchini S, Smith SD, Richardson AJ, Scerri TS, et al. (2004) A 77-kilobase region of chromosome 6p22.2 is associated with dyslexia in families from the United Kingdom and from the United States. Am J Hum Genet 75: 1046–1058.

54. Scerri TS, Paracchini S, Morris A, MacPhie IL, Talcott J, et al. (2010) Identification of candidate genes for dyslexia susceptibility on chromosome 18. PLoS One 5: e13712.

55. Elliot C, Murray D, Pearson L (1983) British Abilities Scales Windsor, UK: NFER-Nelson.

56. Newnham JP, Evans SF, Michael CA, Stanley FJ, Landau LI (1993) Effects of frequent ultrasound during pregnancy: a randomised controlled trial. Lancet 342: 887–891.

57. Dunn L, Dunn L (1981) Peabody Picture Vocabulary Test-Revised: Manual. Circle Pines, MN: American Guidance Services.

58. Raven JC (1976) Ravens Coloured Progressive Matrices. London: Oxford Psychologists Press.

59. Dumontheil I, Roggeman C, Ziermans T, Peyrard-Janvid M, Matsson H, et al. (2011) Influence of the COMT genotype on working memory and brain activity changes during development. Biol Psychiatry 70: 222–229.

60. Soderqvist S, McNab F, Peyrard-Janvid M, Matsson H, Humphreys K, et al. (2010) The SNAP25 gene is linked to working memory capacity and maturation of the posterior cingulate cortex during childhood. Biol Psychiatry 68: 1120–1125.

61. Whitehouse AJ, Bishop DV, Ang QW, Pennell CE, Fisher SE (2011) CNTNAP2 variants affect early language development in the general population. Genes Brain Behav

62. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.

63. Abecasis GR, Cardon LR, Cookson WO (2000) A general test of association for quantitative traits in nuclear families. Am J Hum Genet 66: 279–292.

64. Hayasaka S, Phan KL, Liberzon I, Worsley KJ, Nichols TE (2004) Nonstationary cluster-size inference with random field and permutation methods. Neuroimage 22: 676–687.