Genome-wide association study identifies 74 loci associated with educational attainment

Apparent latent structure within the UK Biobank sample has implications for epidemiological analysis

Large studies use genotype data to discover genetic contributions to complex traits and infer relationships between those traits. Co-incident geographical variation in genotypes and health traits can bias these analyses. Here we show that single genetic variants and genetic scores composed of multiple variants are associated with birth location within UK Biobank and that geographic structure in genotype data cannot be accounted for using routine adjustment for study centre and principal components derived from genotype data. We find that major health outcomes appear geographically structured and that coincident structure in health outcomes and genotype data can yield biased associations. Understanding and accounting for this phenomenon will be important when making inference from genotype data in large studies.

Multivariate genome-wide analyses of the well-being spectrum

We introduce two novel methods for multivariate genome-wide-association meta-analysis (GWAMA) of related traits that correct for sample overlap. A broad range of simulation scenarios supports the added value of our multivariate methods relative to univariate GWAMA. We applied the novel methods to life satisfaction, positive affect, neuroticism, and depressive symptoms, collectively referred to as the well-being spectrum (N_obs = 2,370,390), and found 304 significant independent signals. Our multivariate approaches resulted in a 26% increase in the number of independent signals relative to the four univariate GWAMAs and in an ~57% increase in the predictive power of polygenic risk scores. Supporting transcriptome- and methylome-wide analyses (TWAS and MWAS, respectively) uncovered an additional 17 and 75 independent loci, respectively. Bioinformatic analyses, based on gene expression in brain tissues and cells, showed that genes differentially expressed in the subiculum and GABAergic interneurons are enriched in their effect on the well-being spectrum.

Father Absence and Accelerated Reproductive Development in Non-Hispanic White Women in the United States

Girls who experience father absence in childhood also experience accelerated reproductive development in comparison with peers with present fathers. One hypothesis advanced to explain this empirical pattern is genetic confounding, wherein gene-environment correlation (rGE) causes a spurious relationship between father absence and reproductive timing. We test this hypothesis by constructing polygenic scores for age at menarche and first birth using recently available genome-wide association study results and molecular genetic data on a sample of non-Hispanic white females from the National Longitudinal Study of Adolescent to Adult Health. We find that young women's accelerated menarche polygenic scores are unrelated to their exposure to father absence. In contrast, polygenic scores for earlier age at first birth tend to be higher in young women raised in homes with absent fathers. Nevertheless, father absence and the polygenic scores independently and additively predict reproductive timing. We find no evidence in support of the rGE hypothesis for accelerated menarche and only limited evidence in support of the rGE hypothesis for earlier age at first birth.

Successful Memory Aging

For more than 50 years, psychologists, gerontologists, and, more recently, neuroscientists have considered the possibility of successful aging. How to define successful aging remains debated, but well-preserved age-sensitive cognitive functions, like episodic memory, is an often-suggested criterion. Evidence for successful memory aging comes from cross-sectional and longitudinal studies showing that some older individuals display high and stable levels of performance. Successful memory aging may be accomplished via multiple paths. One path is through brain maintenance, or relative lack of age-related brain pathology. Through another path, successful memory aging can be accomplished despite brain pathology by means of efficient compensatory and strategic processes. Genetic, epigenetic, and lifestyle factors influence memory aging via both paths. Some of these factors can be promoted throughout the life course, which, at the individual as well as the societal level, can positively impact successful memory aging.

Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder

Attention deficit/hyperactivity disorder (ADHD) is a highly heritable childhood behavioral disorder affecting 5% of children and 2.5% of adults. Common genetic variants contribute substantially to ADHD susceptibility, but no variants have been robustly associated with ADHD. We report a genome-wide association meta-analysis of 20,183 individuals diagnosed with ADHD and 35,191 controls that identifies variants surpassing genome-wide significance in 12 independent loci, finding important new information about the underlying biology of ADHD. Associations are enriched in evolutionarily constrained genomic regions and loss-of-function intolerant genes and around brain-expressed regulatory marks. Analyses of three replication studies: a cohort of individuals diagnosed with ADHD, a self-reported ADHD sample and a meta-analysis of quantitative measures of ADHD symptoms in the population, support these findings while highlighting study-specific differences on genetic overlap with educational attainment. Strong concordance with GWAS of quantitative population measures of ADHD symptoms supports that clinical diagnosis of ADHD is an extreme expression of continuous heritable traits.

Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder

Attention deficit/hyperactivity disorder (ADHD) is a highly heritable childhood behavioral disorder affecting 5% of children and 2.5% of adults. Common genetic variants contribute substantially to ADHD susceptibility, but no variants have been robustly associated with ADHD. We report a genome-wide association meta-analysis of 20,183 individuals diagnosed with ADHD and 35,191 controls that identifies variants surpassing genome-wide significance in 12 independent loci, finding important new information about the underlying biology of ADHD. Associations are enriched in evolutionarily constrained genomic regions and loss-of-function intolerant genes and around brain-expressed regulatory marks. Analyses of three replication studies: a cohort of individuals diagnosed with ADHD, a self-reported ADHD sample and a meta-analysis of quantitative measures of ADHD symptoms in the population, support these findings while highlighting study-specific differences on genetic overlap with educational attainment. Strong concordance with GWAS of quantitative population measures of ADHD symptoms supports that clinical diagnosis of ADHD is an extreme expression of continuous heritable traits.

A combined analysis of genetically correlated traits identifies 187 loci and a role for neurogenesis and myelination in intelligence

Intelligence, or general cognitive function, is phenotypically and genetically correlated with many traits, including a wide range of physical, and mental health variables. Education is strongly genetically correlated with intelligence (r_g = 0.70). We used these findings as foundations for our use of a novel approach-multi-trait analysis of genome-wide association studies (MTAG; Turley et al. 2017)-to combine two large genome-wide association studies (GWASs) of education and intelligence, increasing statistical power and resulting in the largest GWAS of intelligence yet reported. Our study had four goals: first, to facilitate the discovery of new genetic loci associated with intelligence; second, to add to our understanding of the biology of intelligence differences; third, to examine whether combining genetically correlated traits in this way produces results consistent with the primary phenotype of intelligence; and, finally, to test how well this new meta-analytic data sample on intelligence predicts phenotypic intelligence in an independent sample. By combining datasets using MTAG, our functional sample size increased from 199,242 participants to 248,482. We found 187 independent loci associated with intelligence, implicating 538 genes, using both SNP-based and gene-based GWAS. We found evidence that neurogenesis and myelination-as well as genes expressed in the synapse, and those involved in the regulation of the nervous system-may explain some of the biological differences in intelligence. The results of our combined analysis demonstrated the same pattern of genetic correlations as those from previous GWASs of intelligence, providing support for the meta-analysis of these genetically-related phenotypes.

New considerations for hiPSC-based models of neuropsychiatric disorders

The development of human-induced pluripotent stem cells (hiPSCs) has made possible patient-specific modeling across the spectrum of human disease. Here, we discuss recent advances in psychiatric genomics and post-mortem studies that provide critical insights concerning cell-type composition and sample size that should be considered when designing hiPSC-based studies of complex genetic disease. We review recent hiPSC-based models of SZ, in light of our new understanding of critical power limitations in the design of hiPSC-based studies of complex genetic disorders. Three possible solutions are a movement towards genetically stratified cohorts of rare variant patients, application of CRISPR technologies to engineer isogenic neural cells to study the impact of common variants, and integration of advanced genetics and hiPSC-based datasets in future studies. Overall, we emphasize that to advance the reproducibility and relevance of hiPSC-based studies, stem cell biologists must contemplate statistical and biological considerations that are already well accepted in the field of genetics. We conclude with a discussion of the hypothesis of biological convergence of disease-through molecular, cellular, circuit, and patient level phenotypes-and how this might emerge through hiPSC-based studies.

Comprehensive functional genomic resource and integrative model for the human brain

Despite progress in defining genetic risk for psychiatric disorders, their molecular mechanisms remain elusive. Addressing this, the PsychENCODE Consortium has generated a comprehensive online resource for the adult brain across 1866 individuals. The PsychENCODE resource contains ~79,000 brain-active enhancers, sets of Hi-C linkages, and topologically associating domains; single-cell expression profiles for many cell types; expression quantitative-trait loci (QTLs); and further QTLs associated with chromatin, splicing, and cell-type proportions. Integration shows that varying cell-type proportions largely account for the cross-population variation in expression (with >88% reconstruction accuracy). It also allows building of a gene regulatory network, linking genome-wide association study variants to genes (e.g., 321 for schizophrenia). We embed this network into an interpretable deep-learning model, which improves disease prediction by ~6-fold versus polygenic risk scores and identifies key genes and pathways in psychiatric disorders.

Transcriptome-wide isoform-level dysregulation in ASD, schizophrenia, and bipolar disorder

Most genetic risk for psychiatric disease lies in regulatory regions, implicating pathogenic dysregulation of gene expression and splicing. However, comprehensive assessments of transcriptomic organization in diseased brains are limited. In this work, we integrated genotypes and RNA sequencing in brain samples from 1695 individuals with autism spectrum disorder (ASD), schizophrenia, and bipolar disorder, as well as controls. More than 25% of the transcriptome exhibits differential splicing or expression, with isoform-level changes capturing the largest disease effects and genetic enrichments. Coexpression networks isolate disease-specific neuronal alterations, as well as microglial, astrocyte, and interferon-response modules defining previously unidentified neural-immune mechanisms. We integrated genetic and genomic data to perform a transcriptome-wide association study, prioritizing disease loci likely mediated by cis effects on brain expression. This transcriptome-wide characterization of the molecular pathology across three major psychiatric disorders provides a comprehensive resource for mechanistic insight and therapeutic development.

GWAS identifies 14 loci for device-measured physical activity and sleep duration

Physical activity and sleep duration are established risk factors for many diseases, but their aetiology is poorly understood, partly due to relying on self-reported evidence. Here we report a genome-wide association study (GWAS) of device-measured physical activity and sleep duration in 91,105 UK Biobank participants, finding 14 significant loci (7 novel). These loci account for 0.06% of activity and 0.39% of sleep duration variation. Genome-wide estimates of ~ 15% phenotypic variation indicate high polygenicity. Heritability is higher in women than men for overall activity (23 vs. 20%, p = 1.5 × 10-4) and sedentary behaviours (18 vs. 15%, p = 9.7 × 10-4). Heritability partitioning, enrichment and pathway analyses indicate the central nervous system plays a role in activity behaviours. Two-sample Mendelian randomisation suggests that increased activity might causally lower diastolic blood pressure (beta mmHg/SD: -0.91, SE = 0.18, p = 8.2 × 10-7), and odds of hypertension (Odds ratio/SD: 0.84, SE = 0.03, p = 4.9 × 10-8). Our results advocate the value of physical activity for reducing blood pressure.

Regulatory and sequence evolution in response to selection for improved associative learning ability in Nasonia vitripennis

Background

Selection acts on the phenotype, yet only the genotype is inherited. While both the phenotypic and genotypic response to short-term selection can be measured, the link between these is a major unsolved problem in evolutionary biology, in particular for complex behavioural phenotypes.

Results

Here we characterize the genomic and the transcriptomic basis of associative learning ability in the parasitic wasp Nasonia vitripennis and use gene network analysis to link the two. We artificially selected for improved associative learning ability in four independent pairs of lines and identified signatures of selection across the genome. Allele frequency diverged consistently between the selected and control lines in 118 single nucleotide polymorphisms (SNPs), clustering in 51 distinct genomic regions containing 128 genes. The majority of SNPs were found in regulatory regions, suggesting a potential role for gene expression evolution. We therefore sequenced the transcriptomes of selected and control lines and identified 36 consistently differentially expressed transcripts with large changes in expression. None of the differentially expressed genes also showed sequence divergence as a result of selection. Instead, gene network analysis showed many of the genes with consistent allele frequency differences and all of the differentially expressed genes to cluster in a single co-expression network. At a functional level, both genomic and transcriptomic analyses implicated members of gene networks known to be involved in neural plasticity and cognitive processes.

Conclusions

Taken together, our results reveal how specific cognitive abilities can readily respond to selection via a complex interplay between regulatory and sequence evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5310-9) contains supplementary material, which is available to authorized users.

SumHer better estimates the SNP heritability of complex traits from summary statistics

We present SumHer, software for estimating confounding bias, SNP heritability, enrichments of heritability and genetic correlations using summary statistics from genome-wide association studies. The key difference between SumHer and the existing software LD Score Regression (LDSC) is that SumHer allows the user to specify the heritability model. We apply SumHer to results from 24 large-scale association studies (average sample size 121,000) using our recommended heritability model. We show that these studies tended to substantially over-correct for confounding, and as a result the number of genome-wide significant loci was under-reported by about a quarter. We also estimate enrichments for 24 categories of SNPs defined by functional annotations. A previous study using LDSC reported that conserved regions were 13-fold enriched, and found a further six categories with above threefold enrichment. By contrast, our analysis using SumHer finds that none of the categories have enrichment above twofold. SumHer provides an improved understanding of the genetic architecture of complex traits, which enables more efficient analysis of future genetic data.

Genomic analysis of family data reveals additional genetic effects on intelligence and personality

Pedigree-based analyses of intelligence have reported that genetic differences account for 50-80% of the phenotypic variation. For personality traits these effects are smaller, with 34-48% of the variance being explained by genetic differences. However, molecular genetic studies using unrelated individuals typically report a heritability estimate of around 30% for intelligence and between 0 and 15% for personality variables. Pedigree-based estimates and molecular genetic estimates may differ because current genotyping platforms are poor at tagging causal variants, variants with low minor allele frequency, copy number variants, and structural variants. Using ~20,000 individuals in the Generation Scotland family cohort genotyped for ~700,000 single-nucleotide polymorphisms (SNPs), we exploit the high levels of linkage disequilibrium (LD) found in members of the same family to quantify the total effect of genetic variants that are not tagged in GWAS of unrelated individuals. In our models, genetic variants in low LD with genotyped SNPs explain over half of the genetic variance in intelligence, education, and neuroticism. By capturing these additional genetic effects our models closely approximate the heritability estimates from twin studies for intelligence and education, but not for neuroticism and extraversion. We then replicated our finding using imputed molecular genetic data from unrelated individuals to show that ~50% of differences in intelligence, and ~40% of the differences in education, can be explained by genetic effects when a larger number of rare SNPs are included. From an evolutionary genetic perspective, a substantial contribution of rare genetic variants to individual differences in intelligence, and education is consistent with mutation-selection balance.

Imprint of assortative mating on the human genome

Preference for mates with similar phenotypes; that is, assortative mating, is widely observed in humans1-5 and has evolutionary consequences6-8. Under Fisher's classical theory6, assortative mating is predicted to induce a signature in the genome at trait-associated loci that can be detected and quantified. Here, we develop and apply a method to quantify assortative mating on a specific trait by estimating the correlation (θ) between genetic predictors of the trait from single nucleotide polymorphisms on odd- versus even-numbered chromosomes. We show by theory and simulation that the effect of assortative mating can be quantified in the presence of population stratification. We applied this approach to 32 complex traits and diseases using single nucleotide polymorphism data from ~400,000 unrelated individuals of European ancestry. We found significant evidence of assortative mating for height (θ = 3.2%) and educational attainment (θ = 2.7%), both of which were consistent with theoretical predictions. Overall, our results imply that assortative mating involves multiple traits and affects the genomic architecture of loci that are associated with these traits, and that the consequence of mate choice can be detected from a random sample of genomes.

Atrial fibrillation genetic risk differentiates cardioembolic stroke from other stroke subtypes

OBJECTIVE

We sought to assess whether genetic risk factors for atrial fibrillation (AF) can explain cardioembolic stroke risk.

METHODS

We evaluated genetic correlations between a previous genetic study of AF and AF in the presence of cardioembolic stroke using genome-wide genotypes from the Stroke Genetics Network (N = 3,190 AF cases, 3,000 cardioembolic stroke cases, and 28,026 referents). We tested whether a previously validated AF polygenic risk score (PRS) associated with cardioembolic and other stroke subtypes after accounting for AF clinical risk factors.

RESULTS

We observed a strong correlation between previously reported genetic risk for AF, AF in the presence of stroke, and cardioembolic stroke (Pearson r = 0.77 and 0.76, respectively, across SNPs with p < 4.4 × 10-4 in the previous AF meta-analysis). An AF PRS, adjusted for clinical AF risk factors, was associated with cardioembolic stroke (odds ratio [OR] per SD = 1.40, p = 1.45 × 10-48), explaining ∼20% of the heritable component of cardioembolic stroke risk. The AF PRS was also associated with stroke of undetermined cause (OR per SD = 1.07, p = 0.004), but no other primary stroke subtypes (all p > 0.1).

CONCLUSIONS

Genetic risk of AF is associated with cardioembolic stroke, independent of clinical risk factors. Studies are warranted to determine whether AF genetic risk can serve as a biomarker for strokes caused by AF.

Fetal Origins of Mental Disorders? An Answer Based on Mendelian Randomization

The Barker hypothesis states that low birth weight (BW) is associated with higher risk of adult onset diseases, including mental disorders like schizophrenia, major depressive disorder (MDD), and attention deficit hyperactivity disorder (ADHD). The main criticism of this hypothesis is that evidence for it comes from observational studies. Specifically, observational evidence does not suffice for inferring causality, because the associations might reflect the effects of confounders. Mendelian randomization (MR) - a novel method that tests causality on the basis of genetic data - creates the unprecedented opportunity to probe the causality in the association between BW and mental disorders in observation studies. We used MR and summary statistics from recent large genome-wide association studies to test whether the association between BW and MDD, schizophrenia and ADHD is causal. We employed the inverse variance weighted (IVW) method in conjunction with several other approaches that are robust to possible assumption violations. MR-Egger was used to rule out horizontal pleiotropy. IVW showed that the association between BW and MDD, schizophrenia and ADHD is not causal (all p > .05). The results of all the other MR methods were similar and highly consistent. MR-Egger provided no evidence for pleiotropic effects biasing the estimates of the effects of BW on MDD (intercept = -0.004, SE = 0.005, p = .372), schizophrenia (intercept = 0.003, SE = 0.01, p = .769), or ADHD (intercept = 0.009, SE = 0.01, p = .357). Based on the current evidence, we refute the Barker hypothesis concerning the fetal origins of adult mental disorders. The discrepancy between our results and the results from observational studies may be explained by the effects of confounders in the observational studies, or by the existence of a small causal effect not detected in our study due to weak instruments. Our power analyses suggested that the upper bound for a potential causal effect of BW on mental disorders would likely not exceed an odds ratio of 1.2.

Genetic Overlap between General Cognitive Function and Schizophrenia: A Review of Cognitive GWASs

General cognitive (intelligence) function is substantially heritable, and is a major determinant of economic and health-related life outcomes. Cognitive impairments and intelligence decline are core features of schizophrenia which are evident before the onset of the illness. Genetic overlaps between cognitive impairments and the vulnerability for the illness have been suggested. Here, we review the literature on recent large-scale genome-wide association studies (GWASs) of general cognitive function and correlations between cognitive function and genetic susceptibility to schizophrenia. In the last decade, large-scale GWASs (n > 30,000) of general cognitive function and schizophrenia have demonstrated that substantial proportions of the heritability of the cognitive function and schizophrenia are explained by a polygenic component consisting of many common genetic variants with small effects. To date, GWASs have identified more than 100 loci linked to general cognitive function and 108 loci linked to schizophrenia. These genetic variants are mostly intronic or intergenic. Genes identified around these genetic variants are densely expressed in brain tissues. Schizophrenia-related genetic risks are consistently correlated with lower general cognitive function (rg = -0.20) and higher educational attainment (rg = 0.08). Cognitive functions are associated with many of the socioeconomic and health-related outcomes. Current treatment strategies largely fail to improve cognitive impairments of schizophrenia. Therefore, further study is needed to understand the molecular mechanisms underlying both cognition and schizophrenia.

Are Bigger Brains Smarter? Evidence From a Large-Scale Preregistered Study

A positive relationship between brain volume and intelligence has been suspected since the 19th century, and empirical studies seem to support this hypothesis. However, this claim is controversial because of concerns about publication bias and the lack of systematic control for critical confounding factors (e.g., height, population structure). We conducted a preregistered study of the relationship between brain volume and cognitive performance using a new sample of adults from the United Kingdom that is about 70% larger than the combined samples of all previous investigations on this subject ( N = 13,608). Our analyses systematically controlled for sex, age, height, socioeconomic status, and population structure, and our analyses were free of publication bias. We found a robust association between total brain volume and fluid intelligence ( r = .19), which is consistent with previous findings in the literature after controlling for measurement quality of intelligence in our data. We also found a positive relationship between total brain volume and educational attainment ( r = .12). These relationships were mainly driven by gray matter (rather than white matter or fluid volume), and effect sizes were similar for both sexes and across age groups.

School Performance and Educational Attainment as Early-Life Predictors of Age-Related Memory Decline: Protective Influences in Later-Born Cohorts

Objectives

Evidence is accumulating that early-life characteristics and experiences contribute significantly to differences in cognitive aging. This study investigated whether school performance at age 12 predicted late-life level and rate of memory change over 15–25 years, and whether its potential protective influence on memory change was mediated by educational attainment or income.

Methods

Latent growth curve models were fitted to 15–25 year longitudinal memory data from a population-based sample, stratified on age cohorts (n = 227, born 1909–1935; n = 301, born 1938–1954).

Results

A latent-level school grade variable significantly predicted both memory level and slope in later-born cohorts. Higher grades were associated with higher level and reduced decline, measured between ages 45 and 70 years, on average. In the earlier-born cohorts, grades predicted memory level, but not slope, measured between ages 66 and 81 years. Follow-up analyses indicated that the protective influence of higher school grades in later-born cohorts was partially mediated by educational attainment, but independent of income.

Discussion

The results suggest that higher childhood school performance is protective against age-related cognitive decline in younger or later-born cohorts, for which further education has been more accessible. Education may exert such influence through increased cognitive reserve or more well-informed health- and lifestyle decisions.

NFIB Haploinsufficiency Is Associated with Intellectual Disability and Macrocephaly

The nuclear factor I (NFI) family of transcription factors play an important role in normal development of multiple organs. Three NFI family members are highly expressed in the brain, and deletions or sequence variants in two of these, NFIA and NFIX, have been associated with intellectual disability (ID) and brain malformations. NFIB, however, has not previously been implicated in human disease. Here, we present a cohort of 18 individuals with mild ID and behavioral issues who are haploinsufficient for NFIB. Ten individuals harbored overlapping microdeletions of the chromosomal 9p23-p22.2 region, ranging in size from 225 kb to 4.3 Mb. Five additional subjects had point sequence variations creating a premature termination codon, and three subjects harbored single-nucleotide variations resulting in an inactive protein as determined using an in vitro reporter assay. All individuals presented with additional variable neurodevelopmental phenotypes, including muscular hypotonia, motor and speech delay, attention deficit disorder, autism spectrum disorder, and behavioral abnormalities. While structural brain anomalies, including dysgenesis of corpus callosum, were variable, individuals most frequently presented with macrocephaly. To determine whether macrocephaly could be a functional consequence of NFIB disruption, we analyzed a cortex-specific Nfib conditional knockout mouse model, which is postnatally viable. Utilizing magnetic resonance imaging and histology, we demonstrate that Nfib conditional knockout mice have enlargement of the cerebral cortex but preservation of overall brain structure and interhemispheric connectivity. Based on our findings, we propose that haploinsufficiency of NFIB causes ID with macrocephaly.

Genome-wide association analysis links multiple psychiatric liability genes to oscillatory brain activity

Oscillatory activity is crucial for information processing in the brain, and has a long history as a biomarker for psychopathology. Variation in oscillatory activity is highly heritable, but current understanding of specific genetic influences remains limited. We performed the largest genome-wide association study to date of oscillatory power during eyes-closed resting electroencephalogram (EEG) across a range of frequencies (delta 1-3.75 Hz, theta 4-7.75 Hz, alpha 8-12.75 Hz, and beta 13-30 Hz) in 8,425 subjects. Additionally, we performed KGG positional gene-based analysis and brain-expression analyses. GABRA2-a known genetic marker for alcohol use disorder and epilepsy-significantly affected beta power, consistent with the known relation between GABAA interneuron activity and beta oscillations. Tissue-specific SNP-based imputation of gene-expression levels based on the GTEx database revealed that hippocampal GABRA2 expression may mediate this effect. Twenty-four genes at 3p21.1 were significant for alpha power (FDR q < .05). SNPs in this region were linked to expression of GLYCTK in hippocampal tissue, and GNL3 and ITIH4 in the frontal cortex-genes that were previously implicated in schizophrenia and bipolar disorder. In sum, we identified several novel genetic variants associated with oscillatory brain activity; furthermore, we replicated and advanced understanding of previously known genes associated with psychopathology (i.e., schizophrenia and alcohol use disorders). Importantly, these psychopathological liability genes affect brain functioning, linking the genes' expression to specific cortical/subcortical brain regions.

Polygenic Risk Scores, School Achievement, and Risk for Schizophrenia: A Danish Population-Based Study

BACKGROUND

Studies have suggested that poor school achievement is associated with increased risk of schizophrenia; however, the possible genetic contribution to this association is unknown. We investigated the possible effect of the polygenic risk score (PRS) for schizophrenia (PRS_SCZ) and for educational attainment (PRS_EDU) on the association between school performance and later schizophrenia.

METHODS

We conducted a case-cohort study on a Danish population-based sample born from 1987 to 1995 comprising 1470 individuals with schizophrenia and 7318 subcohort noncases. Genome-wide data, school performance, and family psychiatric and socioeconomic background information were obtained from national registers and neonatal biobanks. PRS_SCZ and PRS_EDU were calculated using discovery effect size estimates from a meta-analysis of 34,600 cases and 45,968 controls and 293,723 individuals.

RESULTS

Higher PRS_SCZ increased the risk (incidence rate ratio [IRR]: 1.28; 95% confidence interval [CI], 1.19-1.36), whereas higher PRS_EDU decreased the risk of schizophrenia (IRR, 0.87; 95% CI, 0.82-0.92) per standard deviation. Not completing primary school and receiving low school marks were associated with increased risk of schizophrenia (IRR, 2.92; 95% CI, 2.37-3.60; and IRR, 1.58; 95% CI, 1.27-1.97, respectively), which was not confounded by PRS_SCZ or PRS_EDU. Adjusting for social factors and parental psychiatric history, effects of not completing primary school and receiving low school marks were attenuated by up to 25% (IRR, 2.19; 95% CI, 1.75-2.73; and IRR, 1.39; 95% CI, 1.11-1.75, respectively). Increasing PRS_EDU correlated with better school performance (p < .01; R² = 7.6%). PRS_SCZ and PRS_EDU was significantly negatively correlated (r = -.31, p < .01).

CONCLUSIONS

The current PRS did not account for the observed association between primary school performance and risk of schizophrenia.

The genetics of university success

University success, which includes enrolment in and achievement at university, as well as quality of the university, have all been linked to later earnings, health and wellbeing. However, little is known about the causes and correlates of differences in university-level outcomes. Capitalizing on both quantitative and molecular genetic data, we perform the first genetically sensitive investigation of university success with a UK-representative sample of 3,000 genotyped individuals and 3,000 twin pairs. Twin analyses indicate substantial additive genetic influence on university entrance exam achievement (57%), university enrolment (51%), university quality (57%) and university achievement (46%). We find that environmental effects tend to be non-shared, although the shared environment is substantial for university enrolment. Furthermore, using multivariate twin analysis, we show moderate to high genetic correlations between university success variables (0.27-0.76). Analyses using DNA alone also support genetic influence on university success. Indeed, a genome-wide polygenic score, derived from a 2016 genome-wide association study of years of education, predicts up to 5% of the variance in each university success variable. These findings suggest young adults select and modify their educational experiences in part based on their genetic propensities and highlight the potential for DNA-based predictions of real-world outcomes, which will continue to increase in predictive power.

Education can reduce health differences related to genetic risk of obesity

This work investigates whether genetic makeup moderates the effects of education on health. Low statistical power and endogenous measures of environment have been obstacles to the credible estimation of such gene-by-environment interactions. We overcome these obstacles by combining a natural experiment that generated variation in secondary education with polygenic scores for a quarter-million individuals. The additional schooling affected body size, lung function, and blood pressure in middle age. The improvements in body size and lung function were larger for individuals with high genetic predisposition to obesity. As a result, education reduced the gap in unhealthy body size between those in the top and bottom terciles of genetic risk of obesity from 20 to 6 percentage points.

Genetik der allgemeinen kognitiven Fähigkeit

Intelligenz ist eines der bestuntersuchten Konstrukte der empirischen Verhaltenswissenschaften und stellt eine allgemeine geistige Kapazität dar, die unter anderem die Fähigkeit zum schlussfolgernden Denken, zum Lösen neuartiger Probleme, zum abstrakten Denken sowie zum schnellen Lernen umfasst. Diese kognitiven Fähigkeiten spielen eine große Rolle in der Erklärung und Vorhersage individueller Unterschiede in zentralen Bereichen des gesellschaftlichen Lebens, wie Schul- und Bildungserfolg, Berufserfolg, sozioökonomischer Status und Gesundheitsverhalten. Verhaltensgenetische Studien zeigen konsistent, dass genetische Einflüsse einen substanziellen Beitrag zur Erklärung individueller Unterschiede leisten, die über 60 % der Intelligenzunterschiede im Erwachsenenalter erklären. In den letzten Jahren konnten in großen genomweiten Assoziationsstudien mit häufigen genetischen Varianten Hunderte mit Intelligenz assoziierte Loci identifiziert werden sowie über 1300 assoziierte Gene mit differentieller Expression überwiegend im Gehirn. Mehrere Signalwege waren angereichert, vor allen für Neurogenese, Regulation der Entwicklung des Nervensystems sowie der synaptischen Struktur und Aktivität. Die Mehrzahl der assoziierten Loci betraf regulatorische Regionen und interessanterweise lag die Hälfte intronisch. Von den über 1300 Genen überlappen nur 9,2 % mit solchen, die mit monogenen neurokognitiven Störungen assoziiert sind. Insgesamt bestätigen die Befunde ein polygenes Modell Tausender additiver Faktoren, wobei die einzelnen Loci eine sehr geringe Effektstärke aufweisen. Insgesamt erklären die jetzigen Befunde ca. 10 % der Gesamtvarianz des Merkmals. Diese Ergebnisse sind ein wichtiger Ausgangspunkt für zukünftige Forschung sowohl in der Genetik als auch den Verhaltenswissenschaften.

A Further Comment on 'Large-Scale Cognitive GWAS Meta-Analysis Reveals Tissue-Specific Neural Expression and Potential Nootropic Drug Targets' by Lam et al

Lam et al. (2018) respond to a commentary of their paper entitled ‘Large-Scale Cognitive GWAS Meta-Analysis Reveals Tissue-Specific Neural Expression and Potential Nootropic Drug Targets’ Lam et al. (2017). While Lam et al. (2018) have now provided the recommended quality control metrics for their paper, problems remain. Specifically, Lam et al. (2018) do not dispute that the results of their multi-trait analysis of genome-wide association study (MTAG) analysis has produced a phenotype with a genetic correlation of one with three measures of education, but do claim the associations found are specific to the trait of cognitive ability. In this brief paper, it is empirically demonstrated that the phenotype derived by Lam et al. (2017) is more genetically similar to education than cognitive ability. In addition, it is shown that of the genome-wide significant loci identified by Lam et al. (2017) are loci that are associated with education rather than with cognitive ability.

Phenotype-Specific Enrichment of Mendelian Disorder Genes near GWAS Regions across 62 Complex Traits

Although recent studies provide evidence for a common genetic basis between complex traits and Mendelian disorders, a thorough quantification of their overlap in a phenotype-specific manner remains elusive. Here, we have quantified the overlap of genes identified through large-scale genome-wide association studies (GWASs) for 62 complex traits and diseases with genes containing mutations known to cause 20 broad categories of Mendelian disorders. We identified a significant enrichment of genes linked to phenotypically matched Mendelian disorders in GWAS gene sets; of the total 1,240 comparisons, a higher proportion of phenotypically matched or related pairs (n = 50 of 92 [54%]) than phenotypically unmatched pairs (n = 27 of 1,148 [2%]) demonstrated significant overlap, confirming a phenotype-specific enrichment pattern. Further, we observed elevated GWAS effect sizes near genes linked to phenotypically matched Mendelian disorders. Finally, we report examples of GWAS variants localized at the transcription start site or physically interacting with the promoters of genes linked to phenotypically matched Mendelian disorders. Our results are consistent with the hypothesis that genes that are disrupted in Mendelian disorders are dysregulated by non-coding variants in complex traits and demonstrate how leveraging findings from related Mendelian disorders and functional genomic datasets can prioritize genes that are putatively dysregulated by local and distal non-coding GWAS variants.

A genetic perspective on the relationship between eudaimonic -and hedonic well-being

Whether hedonism or eudaimonia are two distinguishable forms of well-being is a topic of ongoing debate. To shed light on the relation between the two, large-scale available molecular genetic data were leveraged to gain more insight into the genetic architecture of the overlap between hedonic and eudaimonic well-being. Hence, we conducted the first genome-wide association studies (GWAS) of eudaimonic well-being (N = ~108 K) and linked it to a GWAS of hedonic well-being (N = ~222 K). We identified the first two genome-wide significant independent loci for eudaimonic well-being and six independent loci for hedonic well-being. Joint analyses revealed a moderate phenotypic correlation (r = 0.53) and a high genetic correlation (rg = 0.78) between eudaimonic and hedonic well-being. This indicates that the genetic etiology of hedonic and eudaimonic well-being is substantially shared, with divergent (environmental) factors contributing to their phenotypic divergence. Loci regulating expression showed significant enrichment in the brain cortex, brain cerebellum, frontal cortex, as well as the cerebellar hemisphere for eudaimonic well-being. No significant enrichment for hedonic well-being is observed, although brain tissues were top ranked. Genetic correlations patterns with a range of positive and negative related phenotypes were largely similar for hedonic -and eudaimonic well-being. Our results reveal a large overlap between the genes that influence hedonism and the genes that influence eudaimonia.

Accurate Genomic Prediction of Human Height

We construct genomic predictors for heritable but extremely complex human quantitative traits (height, heel bone density, and educational attainment) using modern methods in high dimensional statistics (i.e., machine learning). The constructed predictors explain, respectively, ∼40, 20, and 9% of total variance for the three traits, in data not used for training. For example, predicted heights correlate ∼0.65 with actual height; actual heights of most individuals in validation samples are within a few centimeters of the prediction. The proportion of variance explained for height is comparable to the estimated common SNP heritability from genome-wide complex trait analysis (GCTA), and seems to be close to its asymptotic value (i.e., as sample size goes to infinity), suggesting that we have captured most of the heritability for SNPs. Thus, our results close the gap between prediction R-squared and common SNP heritability. The ∼20k activated SNPs in our height predictor reveal the genetic architecture of human height, at least for common variants. Our primary dataset is the UK Biobank cohort, comprised of almost 500k individual genotypes with multiple phenotypes. We also use other datasets and SNPs found in earlier genome-wide association studies (GWAS) for out-of-sample validation of our results.

Detecting genome-wide directional effects of transcription factor binding on polygenic disease risk

Biological interpretation of genome-wide association study data frequently involves assessing whether SNPs linked to a biological process, for example, binding of a transcription factor, show unsigned enrichment for disease signal. However, signed annotations quantifying whether each SNP allele promotes or hinders the biological process can enable stronger statements about disease mechanism. We introduce a method, signed linkage disequilibrium profile regression, for detecting genome-wide directional effects of signed functional annotations on disease risk. We validate the method via simulations and application to molecular quantitative trait loci in blood, recovering known transcriptional regulators. We apply the method to expression quantitative trait loci in 48 Genotype-Tissue Expression tissues, identifying 651 transcription factor-tissue associations including 30 with robust evidence of tissue specificity. We apply the method to 46 diseases and complex traits (average n = 290 K), identifying 77 annotation-trait associations representing 12 independent transcription factor-trait associations, and characterize the underlying transcriptional programs using gene-set enrichment analyses. Our results implicate new causal disease genes and new disease mechanisms.

Does education protect against depression? Evidence from the Young Finns Study using Mendelian randomization

Using participants (N = 1733) drawn from the nationally representative longitudinal Young Finns Study (YFS) we estimate the effect of education on depressive symptoms. In 2007, when the participants were between 30 and 45 years old, they reported their depressive symptoms using a revised version of Beck's Depression Inventory. Education was measured using register information on the highest completed level of education in 2007, which was converted to years of education. To identify a causal relationship between education and depressive symptoms we use an instrumental variables approach (Mendelian randomization, MR) with a genetic risk score as an instrument for years of education. The genetic risk score was based on 74 genetic variants, which were associated with years of education in a genome-wide association study (GWAS). Because the genetic variants are randomly assigned at conception, they induce exogenous variation in years of education and thus identify a causal effect if the assumptions of the MR approach are met. In Ordinary Least Squares (OLS) estimation years of education in 2007 were negatively associated with depressive symptoms in 2007 (b = -0.027, 95% Confidence Interval (CI) = -0.040, -0.015). However, the results based on Mendelian randomization suggested that the effect is not causal (b = 0.017; 95% CI = -0.144, 0.178). This indicates that omitted variables correlated with education and depression may bias the linear regression coefficients and exogenous variation in education caused by differences in genetic make-up does not seem to protect against depressive symptoms.

Testing the validity of value-added measures of educational progress with genetic data

Value-added measures of educational progress have been used by education researchers and policy-makers to assess the performance of teachers and schools, contributing to performance-related pay and position in school league tables. They are designed to control for all underlying differences between pupils and should therefore provide unbiased measures of school and teacher influence on pupil progress, however, their effectiveness has been questioned. We exploit genetic data from a UK birth cohort to investigate how successfully value-added measures control for genetic differences between pupils. We use raw value-added, contextual value-added (which additionally controls for background characteristics) and teacher-reported value-added measures built from data at ages 11, 14 and 16. Sample sizes for analyses range from 4,600 to 6,518. Our findings demonstrate that genetic differences between pupils explain little variation in raw value-added measures but explain up to 20% of the variation in contextual value-added measures (95% CI = 6.06% to 35.71%). Value-added measures built from teacher-rated ability have a greater proportion of variance explained by genetic differences between pupils, with 36.3% of their cross-sectional variation being statistically accounted for by genetics (95% CI = 22.8% to 49.8%). By contrast, a far greater proportion of variance is explained by genetic differences for raw test scores at each age of at least 47.3% (95% CI: 35.9 to 58.7). These findings provide evidence that value-added measures of educational progress can be influenced by genetic differences between pupils, and therefore may provide a biased measure of school and teacher performance. We include a glossary of genetic terms for educational researchers interested in the use of genetic data in educational research.

Understanding the effects of education through the lens of biology

Early educational interventions aim to close gaps in achievement levels between children. However, early interventions do not eliminate individual differences in populations and the effects of early interventions often fade-out over time, despite changes of the mean of the population immediately following the intervention. Here, we discuss biological factors that help to better understand why early educational interventions do not eliminate achievement gaps. Children experience and respond to educational interventions differently. These stable individual differences are a consequence of biological mechanisms that support the interplay between genetic predispositions and the embedding of experience into our biology. Accordingly, we argue that it is not plausible to conceptualize the goals of educational interventions as both a shifting of the mean and a narrowing of the distribution of a particular measure of educational attainment assumed to be of utmost importance (such as a standardized test score). Instead of aiming to equalize the performance of students, the key goal of educational interventions should be to maximize potential at the individual level and consider a kaleidoscope of educational outcomes across which individuals vary. Additionally, in place of employing short-term interventions in the hope of achieving long-term gains, educational interventions need to be sustained throughout development and their long-term, rather than short-term, efficacy be evaluated. In summary, this paper highlights how biological research is valuable for driving a re-evaluation of how educational success across development can be conceptualized and thus what policy implications may be drawn.

Genome-wide association study of primary open-angle glaucoma in continental and admixed African populations

Primary open angle glaucoma (POAG) is a complex disease with a major genetic contribution. Its prevalence varies greatly among ethnic groups, and is up to five times more frequent in black African populations compared to Europeans. So far, worldwide efforts to elucidate the genetic complexity of POAG in African populations has been limited. We conducted a genome-wide association study in 1113 POAG cases and 1826 controls from Tanzanian, South African and African American study samples. Apart from confirming evidence of association at TXNRD2 (rs16984299; OR[T] 1.20; P = 0.003), we found that a genetic risk score combining the effects of the 15 previously reported POAG loci was significantly associated with POAG in our samples (OR 1.56; 95% CI 1.26-1.93; P = 4.79 × 10-5). By genome-wide association testing we identified a novel candidate locus, rs141186647, harboring EXOC4 (OR[A] 0.48; P = 3.75 × 10-8), a gene transcribing a component of the exocyst complex involved in vesicle transport. The low frequency and high degree of genetic heterogeneity at this region hampered validation of this finding in predominantly West-African replication sets. Our results suggest that established genetic risk factors play a role in African POAG, however, they do not explain the higher disease load. The high heterogeneity within Africans remains a challenge to identify the genetic commonalities for POAG in this ethnicity, and demands studies of extremely large size.

Epigenetic prediction of complex traits and death

BACKGROUND

Genome-wide DNA methylation (DNAm) profiling has allowed for the development of molecular predictors for a multitude of traits and diseases. Such predictors may be more accurate than the self-reported phenotypes and could have clinical applications.

RESULTS

Here, penalized regression models are used to develop DNAm predictors for ten modifiable health and lifestyle factors in a cohort of 5087 individuals. Using an independent test cohort comprising 895 individuals, the proportion of phenotypic variance explained in each trait is examined for DNAm-based and genetic predictors. Receiver operator characteristic curves are generated to investigate the predictive performance of DNAm-based predictors, using dichotomized phenotypes. The relationship between DNAm scores and all-cause mortality (n = 212 events) is assessed via Cox proportional hazards models. DNAm predictors for smoking, alcohol, education, and waist-to-hip ratio are shown to predict mortality in multivariate models. The predictors show moderate discrimination of obesity, alcohol consumption, and HDL cholesterol. There is excellent discrimination of current smoking status, poorer discrimination of college-educated individuals and those with high total cholesterol, LDL with remnant cholesterol, and total:HDL cholesterol ratios.

CONCLUSIONS

DNAm predictors correlate with lifestyle factors that are associated with health and mortality. They may supplement DNAm-based predictors of age to identify the lifestyle profiles of individuals and predict disease risk.

Common genetic variants contribute to risk of rare severe neurodevelopmental disorders

There are thousands of rare human disorders that are caused by single deleterious, protein-coding genetic variants¹. However, patients with the same genetic defect can have different clinical presentations^2-4, and some individuals who carry known disease-causing variants can appear unaffected⁵. Here, to understand what explains these differences, we study a cohort of 6,987 children assessed by clinical geneticists to have severe neurodevelopmental disorders such as global developmental delay and autism, often in combination with abnormalities of other organ systems. Although the genetic causes of these neurodevelopmental disorders are expected to be almost entirely monogenic, we show that 7.7% of variance in risk is attributable to inherited common genetic variation. We replicated this genome-wide common variant burden by showing, in an independent sample of 728 trios (comprising a child plus both parents) from the same cohort, that this burden is over-transmitted from parents to children with neurodevelopmental disorders. Our common-variant signal is significantly positively correlated with genetic predisposition to lower educational attainment, decreased intelligence and risk of schizophrenia. We found that common-variant risk was not significantly different between individuals with and without a known protein-coding diagnostic variant, which suggests that common-variant risk affects patients both with and without a monogenic diagnosis. In addition, previously published common-variant scores for autism, height, birth weight and intracranial volume were all correlated with these traits within our cohort, which suggests that phenotypic expression in individuals with monogenic disorders is affected by the same variants as in the general population. Our results demonstrate that common genetic variation affects both overall risk and clinical presentation in neurodevelopmental disorders that are typically considered to be monogenic.

Missing compared to what? Revisiting heritability, genes and culture

Standard models for the determination of phenotypes from genes are grounded in simple assumptions that are inherent in the modern evolutionary synthesis (MES), which was developed in the 1930s, 1940s and 1950s. The MES was framed in the context of Mendelian genetic transmission enhanced by the Fisherian view of the way discretely inherited genes determine continuously quantitative phenotypes. The statistical models that are used to estimate and interpret genetic contributions to human phenotypes-including behavioural traits-are constructed within the framework of the MES. Variance analysis constitutes the main tool and is used under this framework to characterize genetic inheritance, and hence determination of phenotypes. In this essay, we show that cultural inheritance, when incorporated into models for the determination of phenotypes, can sharply reduce estimates of the genetic contribution to these phenotypes. Recognition of the importance of non-genetic transmission of many human traits is becoming ever more necessary to prevent regression to the debates of the 1970s and 1980s concerning policies based on genetic determination of complex human phenotypes.This article is part of the theme issue 'Bridging cultural gaps: interdisciplinary studies in human cultural evolution'.

The accuracy of LD Score regression as an estimator of confounding and genetic correlations in genome-wide association studies

To infer that a single-nucleotide polymorphism (SNP) either affects a phenotype or is linkage disequilibrium with a causal site, we must have some assurance that any SNP-phenotype correlation is not the result of confounding with environmental variables that also affect the trait. In this study, we study the properties of linkage disequilibrium (LD) Score regression, a recently developed method for using summary statistics from genome-wide association studies to ensure that confounding does not inflate the number of false positives. We do not treat the effects of genetic variation as a random variable and thus are able to obtain results about the unbiasedness of this method. We demonstrate that LD Score regression can produce estimates of confounding at null SNPs that are unbiased or conservative under fairly general conditions. This robustness holds in the case of the parent genotype affecting the offspring phenotype through some environmental mechanism, despite the resulting correlation over SNPs between LD Scores and the degree of confounding. Additionally, we demonstrate that LD Score regression can produce reasonably robust estimates of the genetic correlation, even when its estimates of the genetic covariance and the two univariate heritabilities are substantially biased.

Orderliness predicts academic performance: behavioural analysis on campus lifestyle

Quantitative understanding of relationships between students' behavioural patterns and academic performances is a significant step towards personalized education. In contrast to previous studies that were mainly based on questionnaire surveys, recent literature suggests that unobtrusive digital data bring us unprecedented opportunities to study students' lifestyles in the campus. In this paper, we collect behavioural records from undergraduate students' (N = 18 960) smart cards and propose two high-level behavioural characters, orderliness and diligence. The former is a novel entropy-based metric that measures the regularity of campus daily life, which is estimated here based on temporal records of taking showers and having meals. Empirical analyses on such large-scale unobtrusive behavioural data demonstrate that academic performance (GPA) is significantly correlated with orderliness. Furthermore, we show that orderliness is an important feature to predict academic performance, which improves the prediction accuracy even in the presence of students' diligence. Based on these analyses, education administrators could quantitatively understand the major factors leading to excellent or poor performance, detect undesirable abnormal behaviours in time and thus implement effective interventions to better guide students' campus lives at an early stage when necessary.

Systematic Mendelian randomization framework elucidates hundreds of CpG sites which may mediate the influence of genetic variants on disease

We have undertaken a systematic Mendelian randomization (MR) study using methylation quantitative trait loci (meQTL) as genetic instruments to assess the relationship between genetic variation, DNA methylation and 139 complex traits. Using two-sample MR, we identified 1148 associations across 61 traits where genetic variants were associated with both proximal DNA methylation (i.e. cis-meQTL) and complex trait variation (P < 1.39 × 10-08). Joint likelihood mapping provided evidence that the genetic variant which influenced DNA methylation levels for 348 of these associations across 47 traits was also responsible for variation in complex traits. These associations showed a high rate of replication in the BIOS QTL and UK Biobank datasets for 14 selected traits, as 101 of the attempted 128 associations survived multiple testing corrections (P < 3.91 × 10-04). Integrating expression quantitative trait loci (eQTL) data suggested that genetic variants responsible for 306 of the 348 refined meQTL associations also influence gene expression, which indicates a coordinated system of effects that are consistent with causality. CpG sites were enriched for histone mark peaks in tissue types relevant to their associated trait and implicated genes were enriched across relevant biological pathways. Though we are unable to distinguish mediation from horizontal pleiotropy in these analyses, our findings should prove valuable in prioritizing candidate loci where DNA methylation may influence traits and help develop mechanistic insight into the aetiology of complex disease.

The stability of educational achievement across school years is largely explained by genetic factors

Little is known about the etiology of developmental change and continuity in educational achievement. Here, we study achievement from primary school to the end of compulsory education for 6000 twin pairs in the UK-representative Twins Early Development Study sample. Results showed that educational achievement is highly heritable across school years and across subjects studied at school (twin heritability ~60%; SNP heritability ~30%); achievement is highly stable (phenotypic correlations ~0.70 from ages 7 to 16). Twin analyses, applying simplex and common pathway models, showed that genetic factors accounted for most of this stability (70%), even after controlling for intelligence (60%). Shared environmental factors also contributed to the stability, while change was mostly accounted for by individual-specific environmental factors. Polygenic scores, derived from a genome-wide association analysis of adult years of education, also showed stable effects on school achievement. We conclude that the remarkable stability of achievement is largely driven genetically even after accounting for intelligence.

Genetics of self-reported risk-taking behaviour, trans-ethnic consistency and relevance to brain gene expression

Risk-taking behaviour is an important component of several psychiatric disorders, including attention-deficit hyperactivity disorder, schizophrenia and bipolar disorder. Previously, two genetic loci have been associated with self-reported risk taking and significant genetic overlap with psychiatric disorders was identified within a subsample of UK Biobank. Using the white British participants of the full UK Biobank cohort (n = 83,677 risk takers versus 244,662 controls) for our primary analysis, we conducted a genome-wide association study of self-reported risk-taking behaviour. In secondary analyses, we assessed sex-specific effects, trans-ethnic heterogeneity and genetic overlap with psychiatric traits. We also investigated the impact of risk-taking-associated SNPs on both gene expression and structural brain imaging. We identified 10 independent loci for risk-taking behaviour, of which eight were novel and two replicated previous findings. In addition, we found two further sex-specific risk-taking loci. There were strong positive genetic correlations between risk-taking and attention-deficit hyperactivity disorder, bipolar disorder and schizophrenia. Index genetic variants demonstrated effects generally consistent with the discovery analysis in individuals of non-British White, South Asian, African-Caribbean or mixed ethnicity. Polygenic risk scores comprising alleles associated with increased risk taking were associated with lower white matter integrity. Genotype-specific expression pattern analyses highlighted DPYSL5, CGREF1 and C15orf59 as plausible candidate genes. Overall, our findings substantially advance our understanding of the biology of risk-taking behaviour, including the possibility of sex-specific contributions, and reveal consistency across ethnicities. We further highlight several putative novel candidate genes, which may mediate these genetic effects.

Characterizing the Relation Between Expression QTLs and Complex Traits: Exploring the Role of Tissue Specificity

Measurement of gene expression levels and detection of eQTLs (expression quantitative trait loci) are difficult in tissues with limited sample availability, such as the brain. However, eQTL overlap between tissues might be high, which would allow for inference of eQTL functioning in the brain via eQTLs detected in readily accessible tissues, e.g. whole blood. Applying Stratified Linkage Disequilibrium Score Regression (SLDSR), we quantified the enrichment in polygenic signal of blood and brain eQTLs in genome-wide association studies (GWAS) of 11 complex traits. We looked at eQTLs discovered in 44 tissues by the Genotype-Tissue Expression (GTEx) consortium and two other large representative studies, and found no tissue-specific eQTL effects. Next, we integrated the GTEx eQTLs with regions associated with tissue-specific histone modifiers, and interrogated their effect on rheumatoid arthritis and schizophrenia. We observed substantially enriched effects of eQTLs located inside regions bearing modification H3K4me1 on schizophrenia, but not rheumatoid arthritis, and not tissue-specific. Finally, we extracted eQTLs associated with tissue-specific differentially expressed genes and determined their effects on rheumatoid arthritis and schizophrenia, these analysis revealed limited enrichment of eQTLs associated with gene specifically expressed in specific tissues. Our results pointed to strong enrichment of eQTLs in their effect on complex traits, without evidence for tissue-specific effects. Lack of tissue-specificity can be either due to a lack of statistical power or due to the true absence of tissue-specific effects. We conclude that eQTLs are strongly enriched in GWAS signal and that the enrichment is not specific to the eQTL discovery tissue. Until sample sizes for eQTL discovery grow sufficiently large, working with relatively accessible tissues as proxy for eQTL discovery is sensible and restricting lookups for GWAS hits to a specific tissue for which limited samples are available might not be advisable.

Relatedness disequilibrium regression estimates heritability without environmental bias

Heritability measures the proportion of trait variation that is due to genetic inheritance. Measurement of heritability is important in the nature-versus-nurture debate. However, existing estimates of heritability may be biased by environmental effects. Here, we introduce relatedness disequilibrium regression (RDR), a novel method for estimating heritability. RDR avoids most sources of environmental bias by exploiting variation in relatedness due to random Mendelian segregation. We used a sample of 54,888 Icelanders who had both parents genotyped to estimate the heritability of 14 traits, including height (55.4%, s.e. 4.4%) and educational attainment (17.0%, s.e. 9.4%). Our results suggest that some other estimates of heritability may be inflated by environmental effects.

MAP1B mutations cause intellectual disability and extensive white matter deficit

Discovery of coding variants in genes that confer risk of neurodevelopmental disorders is an important step towards understanding the pathophysiology of these disorders. Whole-genome sequencing of 31,463 Icelanders uncovers a frameshift variant (E712KfsTer10) in microtubule-associated protein 1B (MAP1B) that associates with ID/low IQ in a large pedigree (genome-wide corrected P = 0.022). Additional stop-gain variants in MAP1B (E1032Ter and R1664Ter) validate the association with ID and IQ. Carriers have 24% less white matter (WM) volume (β = −2.1SD, P = 5.1 × 10⁻⁸), 47% less corpus callosum (CC) volume (β = −2.4SD, P = 5.5 × 10⁻¹⁰) and lower brain-wide fractional anisotropy (P = 6.7 × 10⁻⁴). In summary, we show that loss of MAP1B function affects general cognitive ability through a profound, brain-wide WM deficit with likely disordered or compromised axons.

Intellectual disability (ID) is characterized by an intelligence quotient of below 70 and impaired adaptive skills. Here, analyzing whole genome sequences from 31,463 Icelanders, Walters et al. identify variants in MAP1B associated with ID and extensive brain-wide white matter deficits.