Evidence for specificity of polygenic contributions to attainment in English, maths and science during adolescence

How well one does at school is predictive of a wide range of important cognitive, socioeconomic, and health outcomes. The last few years have shown marked advancement in our understanding of the genetic contributions to, and correlations with, academic attainment. However, there exists a gap in our understanding of the specificity of genetic associations with performance in academic subjects during adolescence, a critical developmental period. To address this, the Avon Longitudinal Study of Parents and Children was used to conduct genome-wide association studies of standardised national English (N = 5983), maths (N = 6017) and science (N = 6089) tests. High SNP-based heritabilities (h2SNP) for all subjects were found (41-53%). Further, h2SNP for maths and science remained after removing shared variance between subjects or IQ (N = 3197-5895). One genome-wide significant single nucleotide polymorphism (rs952964, p = 4.86 × 10-8) and four gene-level associations with science attainment (MEF2C, BRINP1, S100A1 and S100A13) were identified. Rs952964 remained significant after removing the variance shared between academic subjects. The findings highlight the benefits of using environmentally homogeneous samples for genetic analyses and indicate that finer-grained phenotyping will help build more specific biological models of variance in learning processes and abilities.

RNA sequencing of identical twins discordant for autism reveals blood-based signatures implicating immune and transcriptional dysregulation

Background

A gap exists in our mechanistic understanding of how genetic and environmental risk factors converge at the molecular level to result in the emergence of autism symptoms. We compared blood-based gene expression signatures in identical twins concordant and discordant for autism spectrum condition (ASC) to differentiate genetic and environmentally driven transcription differences, and establish convergent evidence for biological mechanisms involved in ASC.

Methods

Genome-wide gene expression data were generated using RNA-seq on whole blood samples taken from 16 pairs of monozygotic (MZ) twins and seven twin pair members (39 individuals in total), who had been assessed for ASC and autism traits at age 12. Differential expression (DE) analyses were performed between (a) affected and unaffected subjects (N = 36) and (b) within discordant ASC MZ twin pairs (total N = 11) to identify environmental-driven DE. Gene set enrichment and pathway testing was performed on DE gene lists. Finally, an integrative analysis using DNA methylation data aimed to identify genes with consistent evidence for altered regulation in cis.

Results

In the discordant twin analysis, three genes showed evidence for DE at FDR < 10%: IGHG4, EVI2A and SNORD15B. In the case-control analysis, four DE genes were identified at FDR < 10% including IGHG4, PRR13P5, DEPDC1B, and ZNF501. We find enrichment for DE of genes curated in the SFARI human gene database. Pathways showing evidence of enrichment included those related to immune cell signalling and immune response, transcriptional control and cell cycle/proliferation. Integrative methylomic and transcriptomic analysis identified a number of genes showing suggestive evidence for cis dysregulation.

Limitations

Identical twins stably discordant for ASC are rare, and as such the sample size was limited and constrained to the use of peripheral blood tissue for transcriptomic and methylomic profiling. Given these primary limitations, we focused on transcript-level analysis.

Conclusions

Using a cohort of ASC discordant and concordant MZ twins, we add to the growing body of transcriptomic-based evidence for an immune-based component in the molecular aetiology of ASC. Whilst the sample size was limited, the study demonstrates the utility of the discordant MZ twin design combined with multi-omics integration for maximising the potential to identify disease-associated molecular signals.

Genome-Wide Association Study of Latent Cognitive Measures in Adolescence: Genetic Overlap With Intelligence and Education

Individual differences in executive functions (EF) are heritable and predictive of academic attainment (AA). However, little is known about genetic contributions to EFs or their genetic relationship with AA and intelligence. We conducted genome-wide association analyses for processing speed (PS) and the latent EF measures of working memory (WM) and inhibitory control (IC) in 4,611 adolescents from the Avon Longitudinal Study of Parents and Children. While no loci reached genome-wide significance, common genetic variants explained 30% of the variance in WM and 19% in PS. In contrast, we failed to find common genetic contributions to IC. Finally, we examined shared genetic effects between EFs and general intelligence, AA and ADHD. We identified significant genetic correlations between WM, intelligence, and AA. A more specific pattern was observed for PS, with modest genetic overlap with intelligence. Together these findings highlight diversity in the genetic contributions to specific cognitive functions and their genetic relationship with educational and psychiatric outcomes.

Estimation of a significance threshold for epigenome-wide association studies

Epigenome-wide association studies (EWAS) are designed to characterise population-level epigenetic differences across the genome and link them to disease. Most commonly, they assess DNA-methylation status at cytosine-guanine dinucleotide (CpG) sites, using platforms such as the Illumina 450k array that profile a subset of CpGs genome wide. An important challenge in the context of EWAS is determining a significance threshold for declaring a CpG site as differentially methylated, taking multiple testing into account. We used a permutation method to estimate a significance threshold specifically for the 450k array and a simulation extrapolation approach to estimate a genome-wide threshold. These methods were applied to five different EWAS datasets derived from a variety of populations and tissue types. We obtained an estimate of α=2.4×10-7 for the 450k array, and a genome-wide estimate of α=3.6×10-8. We further demonstrate the importance of these results by showing that previously recommended sample sizes for EWAS should be adjusted upwards, requiring samples between ∼10% and ∼20% larger in order to maintain type-1 errors at the desired level.

Tissue-specific patterns of allelically-skewed DNA methylation

While DNA methylation is usually thought to be symmetrical across both alleles, there are some notable exceptions. Genomic imprinting and X chromosome inactivation are two well-studied sources of allele-specific methylation (ASM), but recent research has indicated a more complex pattern in which genotypic variation can be associated with allelically-skewed DNA methylation in cis. Given the known heterogeneity of DNA methylation across tissues and cell types we explored inter- and intra-individual variation in ASM across several regions of the human brain and whole blood from multiple individuals. Consistent with previous studies, we find widespread ASM with > 4% of the ∼220,000 loci interrogated showing evidence of allelically-skewed DNA methylation. We identify ASM flanking known imprinted regions, and show that ASM sites are enriched in DNase I hypersensitivity sites and often located in an extended genomic context of intermediate DNA methylation. We also detect examples of genotype-driven ASM, some of which are tissue-specific. These findings contribute to our understanding of the nature of differential DNA methylation across tissues and have important implications for genetic studies of complex disease. As a resource to the community, ASM patterns across each of the tissues studied are available in a searchable online database: http://epigenetics.essex.ac.uk/ASMBrainBlood.

Common variation near ROBO2 is associated with expressive vocabulary in infancy

Twin studies suggest that expressive vocabulary at ~24 months is modestly heritable. However, the genes influencing this early linguistic phenotype are unknown. Here we conduct a genome-wide screen and follow-up study of expressive vocabulary in toddlers of European descent from up to four studies of the EArly Genetics and Lifecourse Epidemiology consortium, analysing an early (15-18 months, 'one-word stage', N(Total) = 8,889) and a later (24-30 months, 'two-word stage', N(Total)=10,819) phase of language acquisition. For the early phase, one single-nucleotide polymorphism (rs7642482) at 3p12.3 near ROBO2, encoding a conserved axon-binding receptor, reaches the genome-wide significance level (P=1.3 × 10(-8)) in the combined sample. This association links language-related common genetic variation in the general population to a potential autism susceptibility locus and a linkage region for dyslexia, speech-sound disorder and reading. The contribution of common genetic influences is, although modest, supported by genome-wide complex trait analysis (meta-GCTA h(2)(15-18-months) = 0.13, meta-GCTA h(2)(24-30-months) = 0.14) and in concordance with additional twin analysis (5,733 pairs of European descent, h(2)(24-months) = 0.20).

Are genetic risk factors for psychosis also associated with dimension-specific psychotic experiences in adolescence?

Psychosis has been hypothesised to be a continuously distributed quantitative phenotype and disorders such as schizophrenia and bipolar disorder represent its extreme manifestations. Evidence suggests that common genetic variants play an important role in liability to both schizophrenia and bipolar disorder. Here we tested the hypothesis that these common variants would also influence psychotic experiences measured dimensionally in adolescents in the general population. Our aim was to test whether schizophrenia and bipolar disorder polygenic risk scores (PRS), as well as specific single nucleotide polymorphisms (SNPs) previously identified as risk variants for schizophrenia, were associated with adolescent dimension-specific psychotic experiences. Self-reported Paranoia, Hallucinations, Cognitive Disorganisation, Grandiosity, Anhedonia, and Parent-rated Negative Symptoms, as measured by the Specific Psychotic Experiences Questionnaire (SPEQ), were assessed in a community sample of 2,152 16-year-olds. Polygenic risk scores were calculated using estimates of the log of odds ratios from the Psychiatric Genomics Consortium GWAS stage-1 mega-analysis of schizophrenia and bipolar disorder. The polygenic risk analyses yielded no significant associations between schizophrenia and bipolar disorder PRS and the SPEQ measures. The analyses on the 28 individual SNPs previously associated with schizophrenia found that two SNPs in TCF4 returned a significant association with the SPEQ Paranoia dimension, rs17512836 (p-value = 2.57×10⁻⁴) and rs9960767 (p-value = 6.23×10⁻⁴). Replication in an independent sample of 16-year-olds (N = 3,427) assessed using the Psychotic-Like Symptoms Questionnaire (PLIKS-Q), a composite measure of multiple positive psychotic experiences, failed to yield significant results. Future research with PRS derived from larger samples, as well as larger adolescent validation samples, would improve the predictive power to test these hypotheses further. The challenges of relating adult clinical diagnostic constructs such as schizophrenia to adolescent psychotic experiences at a genetic level are discussed.

Genome-wide association study of receptive language ability of 12-year-olds

PURPOSE

Researchers have previously shown that individual differences in measures of receptive language ability at age 12 are highly heritable. In the current study, the authors attempted to identify some of the genes responsible for the heritability of receptive language ability using a genome-wide association approach.

METHOD

The authors administered 4 Internet-based measures of receptive language (vocabulary, semantics, syntax, and pragmatics) to a sample of 2,329 twelve-year-olds for whom DNA and genome-wide genotyping were available. Nearly 700,000 single-nucleotide polymorphisms (SNPs) and 1 million imputed SNPs were included in a genome-wide association analysis of receptive language composite scores.

RESULTS

No SNP associations met the demanding criterion of genome-wide significance that corrects for multiple testing across the genome ( p < 5 × 10 -8). The strongest SNP association did not replicate in an additional sample of 2,639 twelve-year-olds.

CONCLUSIONS

These results indicate that individual differences in receptive language ability in the general population do not reflect common genetic variants that account for more than 3% of the phenotypic variance. The search for genetic variants associated with language skill will require larger samples and additional methods to identify and functionally characterize the full spectrum of risk variants.

Genetics of callous-unemotional behavior in children

Callous-unemotional behavior (CU) is currently under consideration as a subtyping index for conduct disorder diagnosis. Twin studies routinely estimate the heritability of CU as greater than 50%. It is now possible to estimate genetic influence using DNA alone from samples of unrelated individuals, not relying on the assumptions of the twin method. Here we use this new DNA method (implemented in a software package called Genome-wide Complex Trait Analysis, GCTA) for the first time to estimate genetic influence on CU. We also report the first genome-wide association (GWA) study of CU as a quantitative trait. We compare these DNA results to those from twin analyses using the same measure and the same community sample of 2,930 children rated by their teachers at ages 7, 9 and 12. GCTA estimates of heritability were near zero, even though twin analysis of CU in this sample confirmed the high heritability of CU reported in the literature, and even though GCTA estimates of heritability were substantial for cognitive and anthropological traits in this sample. No significant associations were found in GWA analysis, which, like GCTA, only detects additive effects of common DNA variants. The phrase ‘missing heritability’ was coined to refer to the gap between variance associated with DNA variants identified in GWA studies versus twin study heritability. However, GCTA heritability, not twin study heritability, is the ceiling for GWA studies because both GCTA and GWA are limited to the overall additive effects of common DNA variants, whereas twin studies are not. This GCTA ceiling is very low for CU in our study, despite its high twin study heritability estimate. The gap between GCTA and twin study heritabilities will make it challenging to identify genes responsible for the heritability of CU.

First genome-wide association study on anxiety-related behaviours in childhood

BACKGROUND

Twin studies have shown that anxiety in a general population sample of children involves both domain-general and trait-specific genetic effects. For this reason, in an attempt to identify genes responsible for these effects, we investigated domain-general and trait-specific genetic associations in the first genome-wide association (GWA) study on anxiety-related behaviours (ARBs) in childhood.

METHODS

The sample included 2810 7-year-olds drawn from the Twins Early Development Study (TEDS) with data available for parent-rated anxiety and genome-wide DNA markers. The measure was the Anxiety-Related Behaviours Questionnaire (ARBQ), which assesses four anxiety traits and also yields a general anxiety composite. Affymetrix GeneChip 6.0 DNA arrays were used to genotype nearly 700,000 single-nucleotide polymorphisms (SNPs), and IMPUTE v2 was used to impute more than 1 million SNPs. Several GWA associations from this discovery sample were followed up in another TEDS sample of 4804 children. In addition, Genome-wide Complex Trait Analysis (GCTA) was used on the discovery sample, to estimate the total amount of variance in ARBs that can be accounted for by SNPs on the array.

RESULTS

No SNP associations met the demanding criterion of genome-wide significance that corrects for multiple testing across the genome (p<5×10(-8)). Attempts to replicate the top associations did not yield significant results. In contrast to the substantial twin study estimates of heritability which ranged from 0.50 (0.03) to 0.61 (0.01), the GCTA estimates of phenotypic variance accounted for by the SNPs were much lower 0.01 (0.11) to 0.19 (0.12).

CONCLUSIONS

Taken together, these GWAS and GCTA results suggest that anxiety--similar to height, weight and intelligence--is affected by many genetic variants of small effect, but unlike these other prototypical polygenic traits, genetic influence on anxiety is not well tagged by common SNPs.

Common DNA markers can account for more than half of the genetic influence on cognitive abilities

For nearly a century, twin and adoption studies have yielded substantial estimates of heritability for cognitive abilities, although it has proved difficult for genomewide-association studies to identify the genetic variants that account for this heritability (i.e., the missing-heritability problem). However, a new approach, genomewide complex-trait analysis (GCTA), forgoes the identification of individual variants to estimate the total heritability captured by common DNA markers on genotyping arrays. In the same sample of 3,154 pairs of 12-year-old twins, we directly compared twin-study heritability estimates for cognitive abilities (language, verbal, nonverbal, and general) with GCTA estimates captured by 1.7 million DNA markers. We found that DNA markers tagged by the array accounted for .66 of the estimated heritability, reaffirming that cognitive abilities are heritable. Larger sample sizes alone will be sufficient to identify many of the genetic variants that influence cognitive abilities.

Allele-specific methylation in the human genome: implications for genetic studies of complex disease

Across the genome, outside of a small number of known imprinted genes and regions subject to X-inactivation in females, DNA methylation at CpG dinucleotides is often assumed to be complementary across both alleles in a diploid cell. However, recent findings suggest the reality is more complex, with the discovery that allele-specific methylation (ASM) is a common feature across the genome. A key observation is that the majority of ASM is associated with genetic variation in cis, although a noticeable proportion is also non-cis in nature and mediated, for example, by parental origin. ASM appears to be both quantitative, characterized by subtle skewing of DNA methylation between alleles, and heterogeneous, varying across tissues and between individuals. These findings have important implications for complex disease genetics; whilst cis-mediated ASM provides a functional consequence for non-coding genetic variation, heterogeneous and quantitative ASM complicates the identification of disease-associated loci. We propose that non-cis ASM could contribute toward the 'missing heritability' of complex diseases, rendering certain loci hemizygous and masking the direct association between genotype and phenotype. We suggest that the interpretation of results from genomewide association studies can be improved by the incorporation of epi-allelic information, and that in order to fully understand the extent and consequence of ASM in the human genome, a comprehensive sequencing-based analysis of allelic methylation patterns across tissues and individuals is required.

Gene set enrichment; a problem of pathways

Gene Set Enrichment (GSE) is a computational technique which determines whether a priori defined set of genes show statistically significant differential expression between two phenotypes. Currently, the gene sets used for GSE are derived from annotation or pathway databases, which often contain computationally based and unrepresentative data. Here, we propose a novel approach for the generation of comprehensive and biologically derived gene sets, deriving sets through the application of machine learning techniques to gene expression data. These gene sets can be produced for specific tissues, developmental stages or environments. They provide a powerful and functionally meaningful way in which to mine genomewide association and next generation sequencing data in order to identify disease-associated variants and pathways.

In search of genes associated with risk for psychopathic tendencies in children: a two-stage genome-wide association study of pooled DNA

BACKGROUND

Quantitative genetic data from our group indicates that antisocial behaviour (AB) is strongly heritable when coupled with psychopathic, callous-unemotional (CU) personality traits. We have also demonstrated that the genetic influences for AB and CU overlap considerably. We conducted a genome-wide association scan that capitalises on these findings in an attempt to identify quantitative trait loci (QTLs) that may increase risk for psychopathic tendencies (AB+/CU+).

METHODS

Teacher ratings at age 7 were used to screen 8374 twins with available DNA samples for individuals that were high vs. low on both AB and CU. In Stage 1, we screened for allele frequency differences in 642,432 autosomal single-nucleotide polymorphisms (SNPs) using the Affymetrix 6.0 GeneChip with pooled DNA for high-scoring (AB+/CU+) versus low-scoring children (N = approximately 300/group). In Stage 2, we tested the 3000 most strongly associated SNPs from Stage 1 for association in the same direction in a second sample of high- versus low-scoring children from the same twin study (18% co-twins).

RESULTS

Using allele frequencies estimated from pooled DNA, we found suggestive evidence for enrichment of association in the second stage of our two-stage genome-wide association design and focus on reporting the 30 top-ranking SNPs nominally associated with psychopathic tendencies. These SNPs include neurodevelopmental genes such as ROBO2.

CONCLUSIONS

Although none of the SNPs reached genome-wide statistical significance we have generated a list of SNPs that are potentially associated with psychopathic tendencies, which we believe warrant verification and replication in large independent and clinical samples.

A genome-wide association study identifies multiple loci associated with mathematics ability and disability

Numeracy is as important as literacy and exhibits a similar frequency of disability. Although its etiology is relatively poorly understood, quantitative genetic research has demonstrated mathematical ability to be moderately heritable. In this first genome-wide association study (GWAS) of mathematical ability and disability, 10 out of 43 single nucleotide polymorphism (SNP) associations nominated from two high- vs. low-ability (n = 600 10-year-olds each) scans of pooled DNA were validated (P < 0.05) in an individually genotyped sample of (*)2356 individuals spanning the entire distribution of mathematical ability, as assessed by teacher reports and online tests. Although the effects are of the modest sizes now expected for complex traits and require further replication, interesting candidate genes are implicated such as NRCAM which encodes a neuronal cell adhesion molecule. When combined into a set, the 10 SNPs account for 2.9% (F = 56.85; df = 1 and 1881; P = 7.277e-14) of the phenotypic variance. The association is linear across the distribution consistent with a quantitative trait locus (QTL) hypothesis; the third of children in our sample who harbour 10 or more of the 20 risk alleles identified are nearly twice as likely (OR = 1.96; df = 1; P = 3.696e-07) to be in the lowest performing 15% of the distribution. Our results correspond with those of quantitative genetic research in indicating that mathematical ability and disability are influenced by many genes generating small effects across the entire spectrum of ability, implying that more highly powered studies will be needed to detect and replicate these QTL associations.

Allelic skewing of DNA methylation is widespread across the genome

DNA methylation is assumed to be complementary on both alleles across the genome, although there are exceptions, notably in regions subject to genomic imprinting. We present a genome-wide survey of the degree of allelic skewing of DNA methylation with the aim of identifying previously unreported differentially methylated regions (DMRs) associated primarily with genomic imprinting or DNA sequence variation acting in cis. We used SNP microarrays to quantitatively assess allele-specific DNA methylation (ASM) in amplicons covering 7.6% of the human genome following cleavage with a cocktail of methylation-sensitive restriction enzymes (MSREs). Selected findings were verified using bisulfite-mapping and gene-expression analyses, subsequently tested in a second tissue from the same individuals, and replicated in DNA obtained from 30 parent-child trios. Our approach detected clear examples of ASM in the vicinity of known imprinted loci, highlighting the validity of the method. In total, 2,704 (1.5%) of our 183,605 informative and stringently filtered SNPs demonstrate an average relative allele score (RAS) change > or =0.10 following MSRE digestion. In agreement with previous reports, the majority of ASM ( approximately 90%) appears to be cis in nature, and several examples of tissue-specific ASM were identified. Our data show that ASM is a widespread phenomenon, with >35,000 such sites potentially occurring across the genome, and that a spectrum of ASM is likely, with heterogeneity between individuals and across tissues. These findings impact our understanding about the origin of individual phenotypic differences and have implications for genetic studies of complex disease.

Assessing individual differences in genome-wide gene expression in human whole blood: reliability over four hours and stability over 10 months

Studying the causes and correlates of natural variation in gene expression in healthy populations assumes that individual differences in gene expression can be reliably and stably assessed across time. However, this is yet to be established. We examined 4-hour test-retest reliability and 10 month test-retest stability of individual differences in gene expression in ten 12-year-old children. Blood was collected on four occasions: 10 a.m. and 2 p.m. on Day 1 and 10 months later at 10 a.m. and 2 p.m. Total RNA was hybridized to Affymetrix-U133 plus 2.0 arrays. For each probeset, the correlation across individuals between 10 a.m. and 2 p.m. on Day 1 estimates test-retest reliability. We identified 3,414 variable and abundantly expressed probesets whose 4-hour test-retest reliability exceeded .70, a conventionally accepted level of reliability, which we had 80% power to detect. Of the 3,414 reliable probesets, 1,752 were also significantly reliable 10 months later. We assessed the long-term stability of individual differences in gene expression by correlating the average expression level for each probe-set across the two 4-hour assessments on Day 1 with the average level of each probe-set across the two 4-hour assessments 10 months later. 1,291 (73.7%) of the 1,752 probe-sets that reliably detected individual differences across 4 hours on two occasions, 10 months apart, also stably detected individual differences across 10 months. Heritability, as estimated from the MZ twin intraclass correlations, is twice as high for the 1,752 reliable probesets versus all present probesets on the array (0.68 vs 0.34), and is even higher (0.76) for the 1,291 reliable probesets that are also stable across 10 months. The 1,291 probesets that reliably detect individual differences from a single peripheral blood collection and stably detect individual differences over 10 months are promising targets for research on the causes (e.g., eQTLs) and correlates (e.g., psychopathology) of individual differences in gene expression.

Increasing heritability of BMI and stronger associations with the FTO gene over childhood

The growing evidence of health risks associated with the rise in childhood obesity adds to the urgency of understanding the determinants of BMI. Twin analyses on repeated assessments of BMI in a longitudinal sample of >7,000 children indicated that the genetic influence on BMI becomes progressively stronger, with heritability increasing from 0.48 at age 4 to 0.78 at age 11. In the same large twin sample, the association between a common variant in the FTO gene and BMI increased in parallel with the rise in heritability, going from R(2) < 0.001 at age 4 to R(2) = 0.01 at age 11. These findings suggest that expression of FTO may become stronger throughout childhood. Increases in heritability may also be due to children increasingly selecting environments correlated with their genetic propensities.

Quantitative trait locus association scan of early reading disability and ability using pooled DNA and 100K SNP microarrays in a sample of 5760 children

Quantitative genetic research suggests that reading disability is the quantitative extreme of the same genetic and environmental factors responsible for normal variation in reading ability. This finding warrants a quantitative trait locus (QTL) strategy that compares low versus high extremes of the normal distribution of reading in the search for QTLs associated with variation throughout the distribution. A low reading ability group (N=755) and a high reading group (N=747) were selected from a representative UK sample of 7-year-olds assessed on two measures of reading that we have shown to be highly heritable and highly genetically correlated. The low and high reading ability groups were each divided into 10 independent DNA pools and the 20 pools were assayed on 100 K single nucleotide polymorphism (SNP) microarrays to screen for the largest allele frequency differences between the low and high reading ability groups. Seventy five of these nominated SNPs were individually genotyped in an independent sample of low (N=452) and high (N=452) reading ability children selected from a second sample of 4258 7-year-olds. Nine of the seventy-five SNPs were nominally significant (P<0.05) in the predicted direction. These 9 SNPs and 14 other SNPs showing low versus high allele frequency differences in the predicted direction were genotyped in the rest of the second sample to test the QTL hypothesis. Ten SNPs yielded nominally significant linear associations in the expected direction across the distribution of reading ability. However, none of these SNP associations accounted for more than 0.5% of the variance of reading ability, despite 99% power to detect them. We conclude that QTL effect sizes, even for highly heritable common disorders and quantitative traits such as early reading disability and ability, might be much smaller than previously considered.

Reading and Generalist Genes

Twin-study research suggests that many (but not all) of the same genes contribute to genetic influence on diverse learning abilities and disabilities, a hypothesis called generalist genes. This generalist genes hypothesis was tested using a set of 10 DNA markers (single nucleotide polymorphisms [SNPs]) found to be associated with early reading ability in a study of 4,258 7-year-old children that screened 100,000 SNPs. Using the same sample, we show that this early reading SNP set also correlates with other aspects of literacy, components of mathematics, and more general cognitive abilities. These results provide support for the generalist genes hypothesis. Although the effect size of the current SNP set is small, such SNP sets could eventually be used to predict genetic risk for learning disabilities as well as to prescribe genetically tailored intervention and prevention programs.