Author: Davide Piffer. Email: firstname.lastname@example.org
A recent GWAS has examined the additive genetic variance accounting for variation in general cognitive function or fluid g (Davies et al., 2015)These were assessed using a battery of information-processing tests including memory, block design, matrix reasoning, reaction time, letter-number sequencing (Davies et al. , 2015).
Since my use of an educational attainment GWAS has been criticized for being affected by environmental variables and for not being strictly an intelligence measure, I decided to see if I could replicate this result on an independent sample and using different measures, hopefully tapping into a more “culture-free” construct, such as fluid g. The typical reaction to using educational attainment is that it could be influenced by environmental variables correlated to genetic variation (see for example the comments by this reviewer: http://openpsych.net/forum/showthread.php?tid=262&pid=3737#pid3737).
13 SNPs with genome-wide significance (p<5*10-8) were identified (Davies et al., 2015). 10 hits (i.e. the allele with a positive effect on the phenotype) were derived and 3 were ancestral alleles. Table 1 reports the average frequency of the 13 SNPs for the 26 populations in 1000 Genomes and the frequency of the top 10 SNPs with an effect on years on education from Rietveld et al. (2013) . The correlation between the two polygenic scores (e.g. average population frequency of GWAS hits) is very high: r= 0.964. Their correlation to population IQ is also substantial: r= 0.817 and 0.715 for Davies et al, 2015 and Rietveld et al, 2013, respectively.
Table 1. Average frequency of intelligence (fluid g) and education (years of education)-increasing alleles from two independent GWAS.
|Population||Davies et al, 2015. Top 13 SNPs||Rietveld et al., 2013. Top 10 SNPs||IQ|
|Gujarati Indian, Tx||0.391||0.403|
|Indian Telegu, UK||0.280||0.293|
|Mende, Sierra Leone||0.311||0.355||64|
|Mexican in L.A.||0.358||0.370||88|
|Sri Lankan, UK||0.308||0.323||79|
As overrepresentation of derived alleles among GWAS hits is a potential counfound (due to different frequencies of derived alleles among population caused by drift and bottlenecks or GWAS artifacts: see my previous posts for an explanation), a baseline frequency of derived alleles (DAF) was estimated using the 693 SNPs significant for human stature in the largest GWAS to date (Wood et al, 2014).
A multiple regression was ran with population IQ and the two variables (baseline DAF and polygenic score) was ran for the two GWAS hits.
Table 2. Standardized beta coefficients. DAF= derived allele frequency. DP (derived alleles with positive effect on the trait).
|Baseline DAF||Davies DP|
|Rietveld et al., 2013||0.406||0.464|
|Davies et al, 2015||0.307||0.587|
Both polygenic scores emerged as better predictors than baseline DAF. A DAF-calibrated score was calculated by subtracting baseline DAF from the frequency of derived hits. This likely represents selection signal on derived alleles as it controls for evolutionary dynamics such as random drift and population bottlenecks. Since the two population-level polygenic scores were highly correlated (r= 0.953), an average score was computed and is reported in table 3, ranked in descending order. This score is highly correlated to the average of the two polygenic scores obtained using all the SNPs (table 1), r= 0.971. However, the correlation with population IQ is slightly lower, at r= 0.687.
Table 3. DAF-calibrated polygenic scores for derived alleles and average polygenic score. Ranked in descending order. DAF= derived allele frequency.
|Population||DAF-free Derived hits. Rietveld et al, 2013||DAF-free Derived hits. Davies et al, 2013||Average|
|Gujarati Indian, Tx||0.033||0.038||0.036|
|Mende, Sierra Leone||0.021||0.038||0.029|
|Mexican in L.A.||-0.007||0.000||-0.003|
|Sri Lankan, UK||-0.057||-0.042||-0.050|
|Indian Telegu, UK||-0.090||-0.064||-0.077|
We can see that genetic variants increasing fluid intelligence and educational attainment are highly correlated at the population-level, suggesting two things: 1) there are common selection pressures on the two phenotypes or 2) educational attainment is a good proxy for g and the SNPs found by Rietveld et al., 2013 are actually g-related (as was suggested by their replication on g in a sub-sample). The findings in the present study debunk two criticisms of my work: 1) That the observed allele frequency differences were “specific” to educational attainment and not really about intelligence and 2) that derived allele differences caused by GWAS artifacts or random drift could mediate the effects. I showed that the observed effects are not due to different baseline derived allele frequencies, thus ruling this out as a possible confound. A discrepancy with IQ estimates is that East Asians lag behind Europeans and that South Asians and Hispanics don’t perform better than sub-Saharan Africans, a finding that is difficult to explain at present.
Again, we observe a tendency for derived alleles (human-specific mutations or not shared with non-human primates) to be overrepresented among the most significant intelligence GWAS hits, confirming the prediction stemming from the evolutionary fact that intelligence has dramatically increased during human evolution.
Davies, G., Armstrong, N., Bis, J. C., et al. (2015). Genetic contributions to variation in general cognitive function: a meta-analysis of genome-wide association studies in the CHARGE consortium (N=53949).Molecular Psychiatry, 20:183-192. doi: 10.1038/mp.2014.188
Rietveld, C.A., Medland, S.E., Derringer, J., Yang, J., Esko, T., Martin, N.W., et al. (2013). GWAS of 126,559 individuals identifies genetic variants associated with educational attainment. Science, 340, 1467-1471. doi: http://doi.org/10.1126/science.1235488
Wood AR, Esko T, Yang J,et al.: Defining the role of common variation in the genomic and biological architecture of adult human height. Nat Genet. 2014; 46(11): 1173–86.