Single Nucleotide Polymorphisms Associated with Reading Ability Show Connection to Socio-Economic Outcomes

Behavior Genetics, Jul 2017

Impairments in reading and in language have negative consequences on life outcomes, but it is not known to what extent genetic effects influence this association. We constructed polygenic scores for difficulties with language and learning to read from genome-wide data in ~6,600 children, adolescents and young adults, and tested their association with health, socioeconomic outcomes and brain structure measures collected in adults (maximal N = 111,749). Polygenic risk of reading difficulties was associated with reduced income, educational attainment, self-rated health and verbal-numerical reasoning (p < 0.00055). Polygenic risk of language difficulties predicted income (p = 0.0005). The small effect sizes ranged 0.01–0.03 of a standard deviation, but these will increase as genetic studies for reading ability get larger. Polygenic scores for childhood cognitive ability and educational attainment were correlated with polygenic scores of reading and language (up to 0.09 and 0.05, respectively). But when they were included in the prediction models, the observed associations between polygenic reading and adult outcomes mostly remained. This suggests that the pathway from reading ability to social outcomes is not only via associated polygenic loads for general cognitive function and educational attainment. The presence of non-overlapping genetic effect is indicated by the genetic correlations of around 0.40 (childhood intelligence) and 0.70 (educational attainment) with reading ability. Mendelian randomization approaches will be important to dissociate any causal and moderating effects of reading and related traits on social outcomes.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1007%2Fs10519-017-9859-x.pdf

Single Nucleotide Polymorphisms Associated with Reading Ability Show Connection to Socio-Economic Outcomes

Single Nucleotide Polymorphisms Associated with Reading Ability Show Connection to Socio-Economic Outcomes Michelle Luciano 0 1 2 3 4 5 6 Saskia P. Hagenaars 0 1 2 3 4 5 6 Simon R. Cox 0 1 2 3 4 5 6 William David Hill 0 1 2 3 4 5 6 Gail Davies 0 1 2 3 4 5 6 Sarah E. Harris 0 1 2 3 4 5 6 Ian J. Deary 0 1 2 3 4 5 6 David M. Evans 0 1 2 3 4 5 6 Nicholas G. Martin 0 1 2 3 4 5 6 Margaret J. Wright 0 1 2 3 4 5 6 Timothy C. Bates 0 1 2 3 4 5 6 0 Division of Psychiatry, University of Edinburgh , Edinburgh EH10 5HF , UK 1 Department of Psychology, School of Philosophy, Psychology and Language Sciences, Centre for Cognitive Ageing and Cognitive Epidemiology, The University of Edinburgh , 7 George Square, Edinburgh EH8 9JZ , UK 2 Michelle Luciano 3 Queensland Brain Institute, University of Queensland , Brisbane , Australia 4 Genetic Epidemiology, QIMR Berghofer Medical Research Institute , Brisbane , Australia 5 MRC Integrative Epidemiology Unit, University of Bristol , Bristol , UK 6 Medical Genetics Section, University of Edinburgh Centre for Genomic and Experimental Medicine and MRC Institute of Genetics and Molecular Medicine, Western General Hospital , Crewe Road, Edinburgh EH4 2XU , UK Impairments in reading and in language have negative consequences on life outcomes, but it is not known to what extent genetic effects influence this association. We constructed polygenic scores for difficulties with language and learning to read from genome-wide data in ~6,600 children, adolescents and young adults, and tested their association with health, socioeconomic outcomes and brain structure measures collected in adults (maximal N = 111,749). Polygenic risk of reading difficulties was associated with reduced income, educational attainment, self-rated health and verbal-numerical reasoning (p < 0.00055). Polygenic risk of language difficulties predicted income (p= 0.0005). The small effect sizes ranged 0.01-0.03 of a standard deviation, but these will increase as genetic studies for reading ability get larger. Polygenic scores for childhood cognitive ability and educational attainment were correlated with polygenic scores of reading and language (up to 0.09 and Vol.:(011233456789) Cross-trait linkage disequilibrium regression; MRI; UK Biobank; Polygenic scores; Genetic correlation - 0.05, respectively). But when they were included in the prediction models, the observed associations between polygenic reading and adult outcomes mostly remained. This suggests that the pathway from reading ability to social outcomes is not only via associated polygenic loads for general cognitive function and educational attainment. The presence of non-overlapping genetic effect is indicated by the genetic correlations of around 0.40 (childhood intelligence) and 0.70 (educational attainment) with reading ability. Mendelian randomization approaches will be important to dissociate any causal and moderating effects of reading and related traits on social outcomes. Disabilities in reading and language are often detected late in childhood, and are a major source of disadvantage not only in formal education (Young et  al. 2002; Richardson and Wydell 2003) but beyond (Ritchie and Bates 2013) . For example, lower occupational status (Maughan 1995), 4 6 University of Queensland Diamantina Institute, Translational Research Institute, Brisbane, QLD, Australia greater psychological distress (Boetsch et  al. 1996) , and over-representation in the prison system (Svensson 2011) are all associated with poor acquisition of reading. Whether genes underlie the relationship between childhood reading or language abilities and adult life outcomes is not known because there are no longitudinal family studies which span such a long time. The advent of polygenic scores, based on thousands of small molecular genetic effects influencing a trait, allows unprecedented access to genetic pathways because it does not require data collection of the related traits of interest in the same sample. We utilise this method to test the association between genes influencing reading and language traits in development and measures of educational attainment, income, self-rated health, depression, handedness and brain magnetic resonance imaging (MRI) in mid-to-late adulthood. Variation in reading and language abilities is influenced by genetic and environmental determinants (Defries et  al. 1987; Bates et al. 2007; Kovas et al. 2005; Hayiou-Thomas 2008) , with the same genetic factors appearing to be relevant for ability and disability (Harlaar et al. 2005; Haworth and Plomin 2010; Lind et  al. 2010) . These data, largely from twin samples, have been supplemented by major advances in molecular genetic analysis, revealing common genetic variants (i.e., Single Nucleotide Polymorphism or SNPs) in or near genes (e.g., KIAA0319) affecting neuronal migration and raising risk for developmental communication disorders (Meng et  al. 2005; Scerri et  al. 2011) . The era of low-cost genome-wide association (GWA) studies has further contributed to understanding of complex traits, not only with new candidates, but also as a basis to build polygenic risk scores. This method takes advantage of the fact that, independent of the stringent criteria for significance of evidence for association (which with small samples and small effects will fall well short of the genomewide significance criterion), the individual effect weights for each of the thousands of independent SNPs may be summed to form an estimate of risk. These scores may be used to predict the same trait in fully-independent samples, or as a linking function, to evaluate relationships of the trait with other outcomes (Evans et  al. 2009) . It is this latter function which we exploit here. Genome-wide association studies of reading and language traits indicate that, like most psychological traits, they are genetically complex: that is, communication is influenced by many thousands of gene differences, each of very small individual effect. Even among the candidate genes for reading and language that have received the most support (e.g., DCDC2, KIAA0319), effect sizes do not typically exceed 1% (Lind et  al. 2010; Francks et  al. 2004) . GWA studies of reading and language traits indicate that the majority of novel association loci will have effect sizes less than 0.4% (Meaburn et  al. 2008; Gialluisi et  al. 2014; Luciano et  al. 2013) . Polygenic scores, which combine information from SNP effects on a trait, allow us to study the effects of genes influencing a trait without requiring each individual SNP to be significant. Such a method was effectively used, for instance, by Belsky and colleagues (2016) to study the relationship of a polygenic score for higher educational attainment with social mobility and adult economic status. Their polygenic scores for educational attainment were based on results from a very large GWA study of around 126,000 individuals, although their independent prediction sample included only 1037 participants. In the present study, the sample used to generate the polygenic scores for reading and language is modest, although similar-sized GWA studies have successfully predicted cross-trait variation. For instance, polygenic neuroticism (based on a GWA sample of 6268) predicted depressive symptoms (Luciano et al. 2012) and polygenic autism spectrum disorder (GWA N ~6700) predicted cognitive abilities (Clarke et al. 2016) . Less reliable polygenic scores will limit the variance that can be explained in a trait, but the sample we predict into is extremely large, with up to 111,785 participants from UK Biobank, giving us sufficient power to detect the smaller expected effect sizes. Importantly, the GWA study on which we base our reading and language polygenic scores controlled for the effects of non-verbal IQ, so any associations that we uncover should be independent of the effects of intelligence, which is also known to associate with various adult outcome measures of interest. To confirm this, genetic correlations between reading and language traits with childhood intelligence are estimated. Successful prediction by polygenic reading and language scores indirectly validates the results from the prior GWAS on which the scores are based. The results further complement those from GWAS of broader traits, like IQ or educational attainment, because the polygenic score can be tied more specifically to the narrowly defined and measured trait. Methods Participants Participants were drawn from the UK Biobank (http:// www.ukbiobank.ac.uk), an open resource enabling the study of factors influencing disease in mid- to late (40–69 years) adulthood (Sudlow et  al. 2015) . The baseline survey measures (including questionnaire and biological samples) were gathered online and in specific locations in the UK between 2006 and 2010 on 502,655 community residing individuals. At baseline, measures of income, education, self-rated health, depression, handedness, and cognitive ability (reasoning, reaction time) were relevant to the present study (N varied 18,321–111,749). Between 2014 and 2015, a web-based questionnaire, completed remotely, obtained data for a further two cognitive tests (trail making and symbol digit substitution) on a subsample of 23,757 and 26,914, respectively. MRI brain scans were available for a subset of participants (N = 5455), who were imaged around 4 years after initial recruitment. Only a subset of these passed quality control of MRI data and had genetic array data (maximum of N = 1206). DNA was extracted from blood samples and genotyping performed with either the UK BiLEVE array or the UK Biobank axiom array. Standard quality control procedures, including exclusions for gender mismatch and non-British ancestry, were applied. Further description can be found in (Hagenaars et  al. 2016) . PRS were calculated using the observed genotypes. The recoding from numeric (1, 2) allele coding to standard ACGT format is described in (Hagenaars et al. 2016) . UK Biobank received ethical approval from the Research Ethics Committee (REC reference 11/NW/0382). This study has been completed under UK Biobank application 10,279. Measures Socio-economic Income The household average total income before tax was measured on a five-point scale: <£18,000; £18,000 to £30,999; £31,000 to £51,999; £52,000 to £100,000; >£100 000. ‘Do not know’ (N = 4319) and ‘Prefer not to answer’ (N = 10,553) responses were coded as missing. Education Participants were asked to report their qualifications using six categories, with multiple responses being allowed—College or University Degree; Advanced levelor Advanced Supplementary-levels or equivalent; Ordinary-levels (or General Certificate of Secondary Education or equivalent); Certificate of Secondary Education or equivalent; National Vocational Qualifications or Higher National Diploma or Higher National Certificate or equivalent; and Other professional qualifications e.g. nursing, teaching. A binary variable was constructed for use in analysis: college/university degree versus none. Thirteen logic/reasoning-type questions were given, each with a two-minute time restriction. The sum-score of correct response was used as our dependent variable. In the UK Biobank protocol (http://biobank.ctsu.ox.ac.uk/crystal/ field.cgi?id=20016) this variable is labelled ‘fluid intelligence’ although it includes verbal and numerical content. The variable had reasonable internal reliability (Cronbach alpha coefficient of 0.62) given the mixture of item types and test–retest reliability (r = 0.65) over an average 4-year interval (Hagenaars et al. 2016; Lyall et al. 2016) . Reaction time Reaction time (RT) was measured by the speed with which the participant pressed a button in response to a pair of matching symbols. It was described to the participant as akin to the card game ‘Snap’ (http://biobank.ctsu.ox.ac. uk/crystal/field.cgi?id=20023). Mean response time in milliseconds across correct trials was the dependent variable. The task has shown good internal reliability with a Cronbach’s alpha of 0.85, but lower test–retest reliability (r = 0.57) over a mean 4-year interval in UK Biobank (Hagenaars et al. 2016; Lyall et al. 2016) . Trail making This is a computer adapted version of the Trail making test that has two parts: A and B. Trails A required the participant to quickly and accurately link numbers appearing in circles on the screen in numeric order; they do this using the computer mouse, or touchpad of laptop or tablet. In Trails B, numbers and letters are presented in the circles, and the participant is required to alternate between linking numbers in numeric order and letters in alphabetical order. The time taken to correctly link the circles is the dependent variable, and here we focus on Trails B, which reflects executive functioning. Symbol digit In this processing speed task, participants are shown a key at the top of the screen in which eight symbols are paired with the numbers one to eight. Below this key test items are presented in a series of grids containing symbols which require a matched numeric response using the key. Responses are made by clicking on the number pad below the grid, with participants instructed to complete the grid from left to right as quickly as they can within the specified 2-min time period. The dependent variable is the number of correct items. Participants were asked whether they were right- or lefthanded; ambidextrous individuals were excluded from the analysis. Coding was in the positive direction of left-handers. On a four-point scale, participants were asked to rate their overall health as excellent, good, fair or poor, therefore, higher scores reflected poorer health. The experience of depression, both previous and current, was assessed by asking questions about lifetime experience of minor and major depression as detailed online in the UK Biobank touchscreen questionnaire. Further items from the Patient Health Questionnaire (PHQ: Spitzer et  al. 1999) and history of help-seeking for mental health problems was used to categorise depression into recurrent episode status as used in the present study. Details on the classification procedure are contained in Smith et  al. (2013 ); diagnostic criteria were arrived at by agreement of 12 mental health/ cognition researchers and aligned as much as possible with the International Classification of Diseases (ICD-10) and the American Psychiatric Association’s Diagnostic and Statistical Manual (DSM-IV). The following items were relevant for classification: “Looking back over your life, have you ever had a time when you were feeling depressed or down for at least a whole week?”; “Have you ever had a period of time lasting at least 2 days when you were so irritable that you found yourself shouting at people or starting fights or arguments?”; “How many weeks was the longest period when you were feeling depressed or down?”; “How many periods have you had when you were feeling depressed or down for at least a whole week?”; “Have you ever seen a general practitioner (GP) for nerves, anxiety, tension or depression?”; “Have you ever seen a psychiatrist for nerves, anxiety, tension or depression?” Brain magnetic resonance imaging (MRI) Data were acquired on a single standard Siemens Skyra 3T scanner with a standard Siemens 32-channel RF receiver head coil, with the imaging matrix at a downward angle of 16° from the AC-PC line. All data were acquired, processed, quality checked and distributed by UK Biobank as Imaging Derived Phenotypes (IDPs). Full information is publicly accessible on the UK Biobank website in the form of a Protocol (http://biobank.ctsu. ox.ac.uk/crystal/refer.cgi?id=2367), Brain Imaging Documentation (http://biobank.ctsu.ox.ac.uk/crystal/ refer.cgi?id=1977) alongside further processing pipeline information (Miller 2016) . In short, T1-weighted volumes were acquired using a 3D MPRAGE sagittal sequence at 1  mm3 resolution with a 208 × 256 × 256 field of view (FoV). For diffusion MRI, a spin-echo echoplanar imaging (SE-EPI) sequence with 10 T2-weighted (b ≈ 0  s  mm−2) baseline, 50 b = 1000  s  mm−2 and 50 b = 2000  s  mm−2 diffusion-weighted volumes acquired with 100 distinct diffusion-encoding directions and three times multi-slice acquisition was used. 2  mm isotropic voxels were provided with a field of view of 104 × 104  mm, imaging matrix 52 × 52, 72 slices with slice thickness 2 mm. Volumetric and diffusion MRI and tractography processing details can also be found in Cox et  al. (2016 ). For brain volume, we focus on grey matter, white matter and total brain volume; these have been corrected for head size by using a SIENAX-type procedure which applies a scaling factor to the volumes that had been extracted from the normalisation transform matrix resulting from the affine registration of skull tissue between T1-weighted volume and MNI152 space. For diffusivity, we focus on measures of fractional anisotropy (FA) and mean diffusivity (MD) in available regions of interest that have been previously tied to reading and/or language: the inferior longitudinal fasciculus and superior longitudinal fasciculus in separate left and right hemispheres (Elnakib et  al. 2014; Vandermosten et  al. 2012) . Following correction for gradient distortion (http://github.com/Washington-University/Pipelines), head motion and eddy currents (http://fsl.fmrib.ox.ac. uk/fsl/fslwiki/EDDY), the within-voxel tract orientation structure was modelled using BEDPOSTx, followed by probabilistic tractography (using PROBTRACKx: Behrens et al. 2007). Tract-averaged measures of FA and MD were then derived using automatic tract mapping with the AutoPtx plug-in for FSL (de Groot et  al. 2013) . Participant data was excluded from the present analysis if they reported, at contemporaneous medical interview, any of the following conditions: of dementia, Parkinson’s disease or any other chronic degenerative neurological problem (including demyelinating diseases), brain cancer, brain haemorrhage, brain abscess, aneurysm, cerebral palsy, encephalitis, head injury, nervous system infection, head or neurological injury or trauma, or stroke. Statistical analysis Polygenic scores were based on genome-wide meta-analysis results (Luciano et al. 2013) for two measures of reading (word reading, general reading and spelling component) and one measure of language (non-word repetition). The meta-analysis was based on two samples: the Brisbane Adolescent Twin Study (BATS; N = 1177, age range 12–26 years, mean of 18 years for reading tests and 20 for language test) and the Avon Longitudinal Study of Parents and their Children (ALSPAC; N = 5472, aged 8 years for language test and 9 years for reading test). All reading and language measures were controlled for the effects of nonverbal (performance) IQ in the genome-wide analyses. Polygenic scores for childhood (6–18 years) general intelligence (Benyamin et al. 2014) and adult educational attainment (college degree) (Rietveld et  al. 2013) were further created in UK Biobank based on the largest publicly available GWA studies estimated in respective sample sizes of 17,989 and 101,069. They were used to test whether reading and language polygenic scores effects were independent from the polygenic influences on general cognitive ability, and, as expected, whether they correlated with genetic influences on educational attainment. Using PRSice software (Euesden et al. 2015) , five polygenic scores were created for each variable of interest based on the significance value of each SNP from the GWA metaanalysis: p < 0.01, p < 0.05, p < 0.1, p < 0.5, and p < 1. We focus on the results of polygenic scores based on all SNPs that perform as well as or better than scores based on smaller sets (Dudbridge 2013) , but all results can be found in the online supplementary material. Prior to calculating the scores, exclusions were made for: low minor allele frequency (<0.01) SNPs, and SNPs in linkage disequilibrium (r2 > 0.25) using a clumping method within a 250 kb window. The clumping method preferentially selects SNPs showing the greatest association for the variable of interest. For the word reading polygenic scores, calculations were based on 134,035 SNPs, reading and spelling component polygenic scores on 134,315 SNPs, and non-word repetition polygenic scores on 134,266 SNPs. The regression models for these polygenic scores predicting the range of adult outcome variables included additional covariates: age at survey, sex, genotyping batch and array, assessment centre, and the first ten genetic principal components (to correct for population stratification). False discovery rate correction was applied to the 300 tests (i.e., 3 polygenic reading and language traits at 5 threshold levels × 20 phenotypic traits). For any significant tests, a further analysis including childhood intelligence polygenic scores as a predictor were performed to confirm whether the reading/ language polygenic effects were independent of those influencing general cognitive ability. Pearson’s correlations between the reading and language polygenic scores with cognitive ability and educational attainment polygenic scores were estimated and genetic correlations were derived by a cross-trait linkage disequilibrium score regression method as described in Bulik-Sullivan et al. (2015). This latter method only relies on GWA summary statistics and is not biased by sample overlap. Results Data distributions were checked for normality (log transformation applied to RT) and extreme scores identified and removed [0.05% for log RT (long RTs), 0.23% for symbol digit (zero scores possibly indicating non-comprehension of the task) and 0.27% for trail making (long RTs)]. The descriptive statistics for the adult socio-economic, health, cognitive, and MRI traits in UK Biobank are shown in Table S1 in the Supplemental Material online. The correlations among the reading and language polygenic scores were generally lower at the more restrictive polygenic inclusion levels. At the p < 0.01 SNP inclusion level, the non-word repetition polygenic score showed correlations of ~0.05 with the two reading polygenic scores, increasing to ~0.55 when the SNP inclusion threshold was p < 1. Correlations between the two reading polygenic scores ranged 0.47 at the p < 0.01 SNP inclusion level to 0.77 at the p < 1 SNP inclusion level. In comparison, at the p < 1 SNP inclusion level, the childhood intelligence polygenic score correlated 0.094 with the reading and spelling component polygenic score, 0.063 with the word reading polygenic score, and 0.093 with the non-word repetition polygenic score (all p = 2.2 × 10−16). Educational attainment polygenic scores correlated 0.053 with the reading and spelling component, 0.036 with word reading, and 0.043 with non-word repetition polygenic scores (all p = 2.2 × 10−16). Genetic correlations, estimated by linkage disequilibrium regression, between the reading and language measures with childhood intelligence and educational attainment are shown in Table  1. Despite observed limitations in power, a significant genetic correlation was shown for the reading and spelling component and childhood intelligence (0.40) and for all measures with educational attainment (ranging 0.56 for non-word repetition to 0.78 for word reading). The associations between the three polygenic scores (word reading, reading and spelling component, non-word repetition) and the socio-economic, health, cognitive, and MRI traits are shown in Table  2. Regression betas for Phenotypes Childhood intelligence  Reading and spelling component  Word reading  Non-word repetition Educational attainment  Reading and spelling component  Word reading  Non-word repetition Genetic correlation Heritability Mean χ2 Z-score 0.021 0.131 0.098 The heritability Z-score and the mean χ2 indicate the level of power to detect association where a heritability Z-score of >4 and a mean χ2 > 1.02 being considered well powered (Bulik-Sullivan et al. 2015) The models include adjustment for standard covariates. All nominally significant results (p< 0.05) are in bold typeface, although FDR significance was at a p< 5.55 × 10−4 ILF inferior longitudinal fasciculus, SLF superior longitudinal fasciculus The model includes adjustment for standard covariates A sensitivity analysis including an additional 5 population stratification components as covariates produced no change in Beta values, with minor increases in p-values polygenic scores based on more restrictive SNP inclusion criteria can be found in Table  S2 in the Supplemental Material online. 56 of the 300 tests were significant at a nominal p < 0.05 level; FDR correction supported significant association for 24 of these. For the associations reported in Table  2, the reading and spelling component polygenic score significantly predicted educational attainment, income, self-rated health, and verbal-numerical reasoning. The word reading polygenic score was significantly associated with verbal-numerical reasoning. Associations between polygenic scores and the brain MRI phenotypes were of comparable magnitudes, but were available in a much reduced sample, and we did not reliably detect any significant associations. The results for the most consistently associated variables—income, educational attainment, self-rated health, and verbal-numerical reasoning— at all polygenic score SNP inclusion levels are graphically depicted in Figure S1 in the Supplemental Material online. For the reading and spelling component polygenic associations that withstood FDR correction, further regression analyses which adjusted for childhood intelligence polygenic effects (Table 3), and additionally, educational attainment polygenic effects (Table 4) were performed. The results in Table 3 show that the reading and spelling component polygenic effects remained significant although they were smaller than those due to polygenic childhood intelligence. A sensitivity analysis which adjusted for a further five population stratification components (to confirm stability of the parameters) left these betas unchanged. Additional adjustment for educational attainment showed further weakening of the reading and spelling polygenic associations, with the income association no longer significant. The educational attainment polygenic score explained more variance in the adult outcomes than both the reading ability and childhood intelligence polygenic predictors. We examined how genetic risk for reading and language difficulties, assessed using polygenic scores, may, independent of intelligence, influence consequential adult outcomes such as income and health. The major finding was that, expressed in positive terms, polygenic scores acting to increase reading ability were significantly associated with higher educational attainment, higher incomes and greater self-rated health as well as with higher verbal-numerical reasoning. Low to moderate genetic correlations between reading traits and childhood intelligence were found, but accounting for polygenic childhood intelligence effects did not significantly alter the polygenic reading associations with adult outcomes. These data might reflect a causal role of the pathway from reading to long-term social outcomes, consistent with a report of reading and language skills mediating the genetic effects on these social outcomes (Ritchie and Bates 2013; Ritchie et  al. 2013) . The present study could not address causality, but given that reading and language skills are developed early in life a causal direction from these abilities to later life outcomes is a tenable hypothesis if one assumes no genetic confounding from factors like parental care and socio-economic status. The polygenic reading ability associations were reduced when polygenic educational attainment was controlled for in the models. This suggests that shared genetic influences on reading ability and adult outcomes are, in part, mediated by genetic effects on educational attainment. Indeed, the strong genetic correlation between reading and educational attainment indicate a dependence on common genes. However, some unique predictive power of polygenic reading scores remained for education, self-rated health and, especially, verbal-numerical reasoning in these models. The shared polygenic variance possibly related to genetic effects influencing reading ability that are independent of the polygenic load associated with childhood intelligence, and educational attainment. Ritchie and colleagues (2013) found that higher reading achievement at age seven was associated with higher socio-economic status at age 42 in women independent of its effects on intelligence at age 11, academic motivation at age 16 and years of education (these indirect effects were significant in men and women). Given that the reading polygenic effects we tested were independent of childhood intelligence polygenic effects, the correlation between polygenic reading scores and income could reflect such a direct effect between reading and income or an indirect effect mediated by motivation and years of education [consistent with the moderate genetic correlation observed between reading and educational attainment and also consistent with the large genetic correlation, rg = 0.90, between the number of years spent in education and household income (Hill et al. 2016) ]. A larger effect might be expected for a measure of individual’s income rather than the household measure available in the present study. The association we found between polygenic reading scores and verbal-numerical reasoning might also reflect a direct process given that the effect was present when childhood intelligence polygenic scores were also in the model. This is consistent with nonshared environmental effects on reading ability that lead to variation in later intelligence (Ritchie et al. 2015) , although in this case it manifests as genetic variation. In terms of causal mechanism, one possibility is that a gene-environment correlation is operating such that those who are genetically predisposed to a higher reading ability are able to better engage with stimuli present in the educational setting and indeed seek out more challenging cognitive endeavours that propels their learning and resultant cognitive abilities, and in turn, educational attainment. The moderate to strong genetic correlations between reading and language traits with educational attainment suggest that reading and language skills are contributors to genetic variation in education attainment. That is, additive genetic effects do not act directly on educational attainment but on those factors (such as reading skill) that influence whether one completes college or not. There was no association between polygenic scores with processing speed or executive function, which might suggest that the advantage that genes influencing reading confer are on abilities that are largely dependent on culture (within UK Biobank verbalnumerical reasoning is arguably a measure of crystallized intelligence; Hagenaars et al. 2016) . Alternatively, the null findings for the other tests might be due to their poor psychometric properties (Lyall et al. 2016) . With regard to health, correlational studies in children and young adults have shown that the presence of dyslexia or reading difficulties is associated with greater mood disturbance, like depressive symptoms (Mugnaini et al. 2009) . Our findings show that this association is not due to genetic overlap, at least for any enduring relationship between early reading ability to depression in mid-life and later. It is possible that over time, children with reading and/or language difficulties develop strategies to overcome the associated life challenges that might otherwise trigger depression. Notwithstanding, we find no support for genes influencing reading ability being related to recurrent depression. A systematic review (Dewalt et  al. 2004) found that reading ability was related to hospitalisation, general indices of health, various chronic ailments and knowledge concerning health. Two large studies in the US showed an association between low literacy and poor self-reported health (Baker et  al. 1997; Gazmararian et  al. 1999) ; in one of these, the effect of literacy was greater than years of education (Baker et al. 1997). In our study, the reading and spelling component polygenic score was correlated with self-rated health, in the same direction as previous phenotypic reports. One might speculate that this is again an instance of mediated pleiotropy where reading ability affects one’s health literacy, resultant care, and self-reported health. The genetic correlation between years of education and self-rated health in UK Biobank was rg = 0.59 (Harris et al. 2016) , it is possible that part of this covariance is due to polygenic influences on reading ability. Research shows that brain indices of white matter structural connectivity, particularly in a region including the left arcuate fasciculus and the left corona radiata (Vandermosten et al. 2012) , and of grey and white matter volumes (Elnakib et al. 2014) differ between dyslexics and controls. Furthermore, genetic risk variants for reading impairment have been linked to various imaging variables (Eicher and Gruen 2013) . Despite this, our study found no associations between genes influencing reading and language traits with MRI variables or handedness, potentially linked to brain asymmetry. However, the sample size for the MRI analyses (maximum N of 1206) was much smaller than the other analyses (which were at least 20 times larger); notably, the effect sizes for MRI traits were larger. Characteristics of cortical regions, such as volume, surface area and thickness, may also offer valuable insights but are currently not available in this sample. Future work on focussed regions associated with reading and language should be a priority. Our effect sizes were small with one standard deviation change in reading polygenic scores equal to at most a 0.03 of a standard deviation change in verbal-numerical reasoning. The same effects for polygenic childhood intelligence and education attainment were around three times larger, but this was expected given that verbal-numerical reasoning is an index of general cognitive ability and that educational attainment is largely determined by one’s general cognitive ability. The much larger GWA samples on which the latter two scores were based would have also contributed to their better prediction. That we could detect polygenic reading ability effects for various socioeconomic, cognitive, and health measures, suggests that even relatively small GWA samples can capture genetic signal (albeit weak) for reading ability. Nevertheless, much larger GWA studies of reading and language traits are needed to confirm the size and distribution of individual gene effects and to produce more robust polygenic scores which will predict greater variance. They will also allow sufficient statistical power to confirm genetic correlations with other traits, which our study was underpowered to do. Observationally, the lower genetic correlations between reading ability and childhood intelligence than with adult intelligence (via the proxy, educational attainment) suggest that with development reading ability and intelligence (or at least academic achievement) become coupled. Language ability did not show this difference. Given that the reading and language ability GWA analyses controlled for non-verbal IQ, the genetic correlations could be indicative of direct effects of reading/language on intelligence and educational attainment rather than the reverse or pleiotropy. Our study has demonstrated that the genes influencing reading ability in childhood through to young adulthood affect later adult outcomes related to social status, education and health. These may operate via a causal process from reading ability to adult outcome, via an unmeasured variable (e.g., motivation) affecting both reading ability and adult outcome, or as direct pleiotropic gene effects on reading ability and adult outcome. Mendelian randomization studies may be a way forward to disentangle the unresolved causal relationships reflected in this genetic overlap, although it will be important to model known mediators or use a multi-phenotype approach. Funding This research has been conducted using the UK Biobank Resource. This work was supported by The University of Edinburgh Centre for Cognitive Ageing and Cognitive Epidemiology, part of the cross council Lifelong Health and Wellbeing Initiative (MR/ K026992/1). Funding from the Biotechnology and Biological Sciences Research Council (BBSRC) and Medical Research Council (MRC) is gratefully acknowledged. W.D.H. is supported by a grant from Age UK (Disconnected Mind Project). Author contributions M Luciano and TC Bates conceived the study. UK Biobank data were prepared for use by SP Hagenaars, SR Cox, G Davies, and SE Harris. DM Evans, NG Martin, MJ Wright managed the data collection and/or analysis of the genomewide study. M Luciano and WD Hill performed data analysis and M Luciano drafted the manuscript with all authors providing revisions. All authors approved the final version of the manuscript for submission. Compliance with ethical standards Conflict of interest Michelle Luciano, Saskia P. Hagenaars, Simon R. Cox, William David Hill, Gail Davies, Sarah E. Harris, Ian J. Deary, David M. Evans, Nicholas G. Martin, Margaret J. Wright, and Timothy C. Bates declare that they have no conflict of interest. Ethical approval All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. Informed consent Informed consent was obtained from all individual participants included in the study. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Baker DW , Parker RM , Williams MV , Clark WS , Nurss J ( 1997 ) The relationship of patient reading ability to self-reported health and use of health services . Am J Public Health 87 ( 6 ): 1027 - 1030 Bates TC , Castles A , Luciano M , Wright MJ , Coltheart M , Martin NG ( 2007 ) Genetic and environmental bases of reading and spelling: a unified genetic dual route model . Read Writ 20 ( 1-2 ): 147 - 171 Behrens TE , Berg HJ , Jbabdi S , Rushworth MF , Woolrich MW ( 2007 ) Probabilistic diffusion tractography with multiple fibre orientations: what can we gain? NeuroImage 34 ( 1 ): 144 - 155 Belsky DW , Moffitt TE , Corcoran DL , Domingue B , Harrington H , Hogan S et al ( 2016 ) The genetics of success: how single-nucleotide polymorphisms associated with educational attainment relate to life-course development . Psychol Sci 27 ( 7 ): 957 - 972 Benyamin B , Pourcain B , Davis OS , Davies G , Hansell NK , Brion MJ et al ( 2014 ) Childhood intelligence is heritable, highly polygenic and associated with FNBP1L . Mol Psychiatry 19 ( 2 ): 253 - 258 Boetsch EA , Green PA , Pennington BF ( 1996 ) Psychosocial correlates of dyslexia across the life span . Dev Psychopathol 8 ( 3 ): 539 - 562 Bulik-Sullivan B , Finucane HK , Anttila V , Gusev A , Day FR , Loh PR et al ( 2015 ) An atlas of genetic correlations across human diseases and traits . Nat Genet 47 ( 11 ): 1236 - 1241 Clarke TK , Lupton MK , Fernandez-Pujals AM , Starr J , Davies G , Cox S et  al ( 2016 ) Common polygenic risk for autism spectrum disorder (ASD) is associated with cognitive ability in the general population . Mol Psychiatry 21 ( 3 ): 419 - 425 Cox SR , Ritchie SJ , Tucker-Drob EM , Liewald DC , Hagenaars SP , Davies G et al ( 2016 ) Ageing and brain white matter structure in 3,513 UK Biobank participants . Nat Commun 7 : 13629 . doi: 10 .1038/ncomms13629 de Groot M , Vernooij MW , Klein S , Ikram MA , Vos FM , Smith SM et  al ( 2013 ) Improving alignment in tract-based spatial statistics: evaluation and optimization of image registration . NeuroImage 76 : 400 - 411 Defries JC , Fulker DW , Labuda MC ( 1987 ) Evidence for a genetic etiology in reading-disability of twins . Nature 329 ( 6139 ): 537 - 539 Dewalt DA , Berkman ND , Sheridan S , Lohr KN , Pignone MP ( 2004 ) Literacy and health outcomes: a systematic review of the literature . J General Intern Med 19 ( 12 ): 1228 - 1239 Dudbridge F ( 2013 ) Power and predictive accuracy of polygenic risk Scores . PLoS Genet 9 ( 3 ): e1003348 Eicher JD , Gruen JR ( 2013 ) Imaging-genetics in dyslexia: connecting risk genetic variants to brain neuroimaging and ultimately to reading impairments . Mol Genet Metab 110 ( 3 ): 201 - 212 Elnakib A , Soliman A , Nitzken M , Casanova MF , Gimel'farb G , El-Baz A ( 2014 ) Magnetic resonance imaging findings for dyslexia: a review . J Biomed Nanotechnol 10 ( 10 ): 2778 - 2805 Euesden J , Lewis CM , O'Reilly PF ( 2015 ) PRSice: polygenic risk score software . Bioinformatics 31 (9): 1466 - 1468 Evans DM , Visscher PM , Wray NR ( 2009 ) Harnessing the information contained within genome-wide association studies to improve individual prediction of complex disease risk . Hum Mol Genet 18 ( 18 ): 3525 - 3531 Francks C , Paracchini S , Smith SD , Richardson AJ , Scerri TS , Cardon LR 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 ( 6 ): 1046 - 1058 Gazmararian JA , Baker DW , Williams MV , Parker RM , Scott TL , Green DC et al ( 1999 ) Health literacy among Medicare enrollees in a managed care organization . JAMA 281 ( 6 ): 545 - 551 Gialluisi A , Newbury DF , Wilcutt EG , Olson RK , DeFries JC , Brandler WM et  al ( 2014 ) Genome-wide screening for DNA variants associated with reading and language traits . Genes, Brain, Behav 13 ( 7 ): 686 - 701 Hagenaars SP , Harris SE , Davies G , Hill WD , Liewald DCM , Ritchie SJ et al ( 2016 ) Shared genetic aetiology between cognitive functions and physical and mental health in UK Biobank (N = 112[thinsp]151) and 24 GWAS consortia . Mol Psychiatry . doi: 10 .1038/mp. 2015 .225 Harlaar N , Spinath FM , Dale PS , Plomin R ( 2005 ) Genetic influences on early word recognition abilities and disabilities: a study of 7-year-old twins . J Child Psychol Psychiatry 46 ( 4 ): 373 - 384 Harris SE , Hagenaars SP , Davies G , Hill WD , Liewald DCM , Ritchie SJ et  al ( 2016 ) Molecular genetic contributions to self-rated health . Int J Epidemiol. doi:10 .1093/ije/dyw219 Haworth CMA , Plomin R ( 2010 ) Quantitative genetics in the era of molecular genetics: learning abilities and disabilities as an example . J Am Acad Child Psy 49 ( 8 ): 783 - 793 Hayiou-Thomas ME ( 2008 ) Genetic and environmental influences on early speech, language and literacy development . J Commun Disord 41 ( 5 ): 397 - 408 Hill WD , Hagenaars SP , Marioni RE , Harris SE , Liewald DC , Davies G et  al ( 2016 ) Molecular genetic contributions to social deprivation and household income in UK Biobank . Curr Biol 26 ( 22 ): 3083 - 3089 Kovas Y , Hayiou-Thomas ME , Oliver B , Dale PS , Bishop DVM , Plomin R ( 2005 ) Genetic influences in different aspects of language development: The etiology of language skills in 4.5-year-old twins . Child Dev 76 ( 3 ): 632 - 651 Lind PA , Luciano M , Wright MJ , Montgomery GW , Martin NG , Bates TC ( 2010 ) Dyslexia and DCDC2: normal variation in reading and spelling is associated with DCDC2 polymorphisms in an Australian population sample . Eur J Hum Genet 18 ( 6 ): 668 - 673 Luciano M , Huffman JE , Arias-Vasquez A , Vinkhuyzen AA , Middeldorp CM , Giegling I et  al ( 2012 ) Genome-wide association uncovers shared genetic effects among personality traits and mood states . Am J Med Genet Part B Neuropsychiatr Genet 159B ( 6 ): 684 - 695 Luciano M , Evans DM , Hansell NK , Medland SE , Montgomery GW , Martin NG et  al ( 2013 ) A genome-wide association study for reading and language abilities in two population cohorts . Genes, Brain, Behav 12 ( 6 ): 645 - 652 Lyall DM , Cullen B , Allerhand M , Smith DJ , Mackay D , Evans J et al ( 2016 ) Cognitive test scores in UK Biobank: data reduction in 480,416 participants and longitudinal stability in 20,346 participants . PloS ONE 11 ( 4 ): e0154222 Maughan B ( 1995 ) Annotation: long-term outcomes of developmental reading problems . J Child Psychol Psychiatry 36 ( 3 ): 357 - 371 Meaburn EL , Harlaar N , Craig IW , Schalkwyk LC , Plomin R ( 2008 ) Quantitative trait locus association scan of early reading disability and ability using pooled DNA and 100 K SNP microarrays in a sample of 5760 children . Mol Psychiatry 13 ( 7 ): 729 - 740 Meng HY , Smith SD , Hager K , Held M , Liu J , Olson RK et  al ( 2005 ) DCDC2 is associated with reading disability and modulates neuronal development in the brain . P Natl Acad Sci USA 102 ( 47 ): 17053 - 17058 Miller KL , Alfaro-Almagro F , Bangerter NK , Thomas DL , Yacoub E , Xu J , Griffanti L ( 2016 ) Multimodal population brain imaging in the UK Biobank prospective epidemiological study . Natneurosci 19 ( 11 ): 1523 - 1536 Mugnaini D , Lassi S , La Malfa G , Albertini G ( 2009 ) Internalizing correlates of dyslexia . World J Pediatr 5 ( 4 ): 255 - 264 Richardson JTE , Wydell TN ( 2003 ) The representation and attainment of students with dyslexia in UK higher education . Read Writ 16 ( 5 ): 475 - 503 Rietveld CA , Medland SE , Derringer J , Yang J , Esko T , Martin NW et  al ( 2013 ) GWAS of 126,559 individuals identifies genetic variants associated with educational attainment . Science 340 ( 6139 ): 1467 - 1471 Ritchie SJ , Bates TC ( 2013 ) Enduring links from childhood mathematics and reading achievement to adult socioeconomic status . Psychol Sci 24 ( 7 ): 1301 - 1308 Ritchie SJ , Bates TC , Plomin R ( 2013 ) Influences of reading ability on the development of intelligence: evidence from a longitudinal monozygotic twin-discordant design . Behav Genet 43 ( 6 ): 538 - 538 Ritchie SJ , Bates TC , Plomin R ( 2015 ) Does learning to read improve intelligence? A longitudinal multivariate analysis in identical twins from age 7 to 16 . Child Dev 86 ( 1 ): 23 - 36 Scerri TS , Morris AP , Buckingham LL , Newbury DF , Miller LL , Monaco AP et  al ( 2011 ) DCDC2, KIAA0319 and CMIP are associated with reading-related traits . Biol Psychiat 70 ( 3 ): 237 - 245 Smith DJ , Nicholl BI , Cullen B , Martin D , Ul-Haq Z , Evans J et  al ( 2013 ) Prevalence and characteristics of probable major depression and bipolar disorder within UK biobank: cross-sectional study of 172,751 participants . PloS one 8 ( 11 ): e75362 Spitzer RL , Kroenke K , Williams JB ( 1999 ) Validation and utility of a self-report version of PRIME-MD: the PHQ primary care study. Primary care evaluation of mental disorders. patient health questionnaire . JAMA 282 ( 18 ): 1737 - 1744 Sudlow C , Gallacher J , Allen N , Beral V , Burton P , Danesh J et  al ( 2015 ) UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age . PLoS Med 12 ( 3 ): e1001779 Svensson I ( 2011 ) Reading and writing disabilities among inmates in correctional settings. A Swedish perspective . Learn Individ Differ 21 ( 1 ): 19 - 29 Vandermosten M , Boets B , Wouters J , Ghesquiere P ( 2012 ) A qualitative and quantitative review of diffusion tensor imaging studies in reading and dyslexia . Neurosci Biobehav Rev 36 ( 6 ): 1532 - 1552 Young AR , Beitchman JH , Johnson C , Douglas L , Atkinson L , Escobar M , Wilson B ( 2002 ) Young adult academic outcomes in a longitudinal sample of early identified language impaired and control children . J Child Psychol Psychiatry 43 ( 5 ): 635 - 645


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1007%2Fs10519-017-9859-x.pdf

Michelle Luciano, Saskia P. Hagenaars, Simon R. Cox, William David Hill, Gail Davies, Sarah E. Harris, Ian J. Deary, David M. Evans, Nicholas G. Martin, Margaret J. Wright, Timothy C. Bates. Single Nucleotide Polymorphisms Associated with Reading Ability Show Connection to Socio-Economic Outcomes, Behavior Genetics, 2017, 1-11, DOI: 10.1007/s10519-017-9859-x