Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity

STUDY QUESTION

Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins?

SUMMARY ANSWER

We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3.

WHAT IS KNOWN ALREADY

The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures.

STUDY DESIGN, SIZE, DURATION

We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls).

PARTICIPANTS/MATERIALS, SETTING, METHODS

Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation.

MAIN RESULTS AND THE ROLE OF CHANCE

This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle.

LARGE SCALE DATA

The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/).

LIMITATIONS, REASONS FOR CAUTION

Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility.

WIDER IMPLICATIONS OF THE FINDINGS

About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility.

STUDY FUNDING/COMPETING INTEREST(S)

Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Skłodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests.

TRIAL REGISTRATION NUMBER

N/A.

A genome-wide association study of total child psychiatric problems scores

Substantial genetic correlations have been reported across psychiatric disorders and numerous cross-disorder genetic variants have been detected. To identify the genetic variants underlying general psychopathology in childhood, we performed a genome-wide association study using a total psychiatric problem score. We analyzed 6,844,199 common SNPs in 38,418 school-aged children from 20 population-based cohorts participating in the EAGLE consortium. The SNP heritability of total psychiatric problems was 5.4% (SE = 0.01) and two loci reached genome-wide significance: rs10767094 and rs202005905. We also observed an association of SBF2, a gene associated with neuroticism in previous GWAS, with total psychiatric problems. The genetic effects underlying the total score were shared with common psychiatric disorders only (attention-deficit/hyperactivity disorder, anxiety, depression, insomnia) (rG > 0.49), but not with autism or the less common adult disorders (schizophrenia, bipolar disorder, or eating disorders) (rG < 0.01). Importantly, the total psychiatric problem score also showed at least a moderate genetic correlation with intelligence, educational attainment, wellbeing, smoking, and body fat (rG > 0.29). The results suggest that many common genetic variants are associated with childhood psychiatric symptoms and related phenotypes in general instead of with specific symptoms. Further research is needed to establish causality and pleiotropic mechanisms between related traits.

Genetic sensitivity analysis: Adjusting for genetic confounding in epidemiological associations

Associations between exposures and outcomes reported in epidemiological studies are typically unadjusted for genetic confounding. We propose a two-stage approach for estimating the degree to which such observed associations can be explained by genetic confounding. First, we assess attenuation of exposure effects in regressions controlling for increasingly powerful polygenic scores. Second, we use structural equation models to estimate genetic confounding using heritability estimates derived from both SNP-based and twin-based studies. We examine associations between maternal education and three developmental outcomes - child educational achievement, Body Mass Index, and Attention Deficit Hyperactivity Disorder. Polygenic scores explain between 14.3% and 23.0% of the original associations, while analyses under SNP- and twin-based heritability scenarios indicate that observed associations could be almost entirely explained by genetic confounding. Thus, caution is needed when interpreting associations from non-genetically informed epidemiology studies. Our approach, akin to a genetically informed sensitivity analysis can be applied widely.

Genetic Meta-Analysis of Twin Birth Weight Shows High Genetic Correlation with Singleton Birth Weight

Birth weight (BW) is an important predictor of newborn survival and health and has associations with many adult health outcomes, including cardio-metabolic disorders, autoimmune diseases, and mental health. On average, twins have a lower BW than singletons as a result of a different pattern of fetal growth and shorter gestational duration. Therefore, investigations into the genetics of BW often exclude data from twins, leading to a reduction in sample size and remaining ambiguities concerning the genetic contribution to BW in twins. In this study, we carried out a genome-wide association meta-analysis of BW in 42 212 twin individuals and found a positive correlation of beta values (Pearson's r = 0.66, 95% confidence interval [CI]: 0.47-0.77) with 150 previously reported genome-wide significant variants for singleton BW. We identified strong positive genetic correlations between BW in twins and numerous anthropometric traits, most notably with BW in singletons (genetic correlation [rg] = 0.92, 95% CI: 0.66-1.18). Genetic correlations of BW in twins with a series of health-related traits closely resembled those previously observed for BW in singletons. Polygenic scores constructed from a genome-wide association study on BW in UK Biobank demonstrated strong predictive power in a target sample of Dutch twins and singletons. Together, our results indicate that a similar genetic architecture underlies BW in twins and singletons and that future genome-wide studies might benefit from including data from large twin registers.

Integrated Polygenic Tool Substantially Enhances Coronary Artery Disease Prediction

Supplemental Digital Content is available in the text.

Background:

There is considerable interest in whether genetic data can be used to improve standard cardiovascular disease risk calculators, as the latter are routinely used in clinical practice to manage preventative treatment.

Methods:

Using the UK Biobank resource, we developed our own polygenic risk score for coronary artery disease (CAD). We used an additional 60 000 UK Biobank individuals to develop an integrated risk tool (IRT) that combined our polygenic risk score with established risk tools (either the American Heart Association/American College of Cardiology pooled cohort equations [PCE] or UK QRISK3), and we tested our IRT in an additional, independent set of 186 451 UK Biobank individuals.

Results:

The novel CAD polygenic risk score shows superior predictive power for CAD events, compared with other published polygenic risk scores, and is largely uncorrelated with PCE and QRISK3. When combined with PCE into an IRT, it has superior predictive accuracy. Overall, 10.4% of incident CAD cases were misclassified as low risk by PCE and correctly classified as high risk by the IRT, compared with 4.4% misclassified by the IRT and correctly classified by PCE. The overall net reclassification improvement for the IRT was 5.9% (95% CI, 4.7–7.0). When individuals were stratified into age-by-sex subgroups, the improvement was larger for all subgroups (range, 8.3%–15.4%), with the best performance in 40- to 54-year-old men (15.4% [95% CI, 11.6–19.3]). Comparable results were found using a different risk tool (QRISK3) and also a broader definition of cardiovascular disease. Use of the IRT is estimated to avoid up to 12 000 deaths in the United States over a 5-year period.

Conclusions:

An IRT that includes polygenic risk outperforms current risk stratification tools and offers greater opportunity for early interventions. Given the plummeting costs of genetic tests, future iterations of CAD risk tools would be enhanced with the addition of a person’s polygenic risk.

Harmonizing behavioral outcomes across studies, raters, and countries: application to the genetic analysis of aggression in the ACTION Consortium

BACKGROUND

Aggression in children has genetic and environmental causes. Studies of aggression can pool existing datasets to include more complex models of social effects. Such analyses require large datasets with harmonized outcome measures. Here, we made use of a reference panel for phenotype data to harmonize multiple aggression measures in school-aged children to jointly analyze data from five large twin cohorts.

METHODS

Individual level aggression data on 86,559 children (42,468 twin pairs) were available in five European twin cohorts measured by different instruments. A phenotypic reference panel was collected which enabled a model-based phenotype harmonization approach. A bi-factor integration model in the integrative data analysis framework was developed to model aggression across studies while adjusting for rater, age, and sex. Finally, harmonized aggression scores were analyzed to estimate contributions of genes, environment, and social interaction to aggression. The large sample size allowed adequate power to test for sibling interaction effects, with unique dynamics permitted for opposite-sex twins.

RESULTS

The best-fitting model found a high level of overall heritability of aggression (~60%). Different heritability rates of aggression across sex were marginally significant, with heritability estimates in boys of ~64% and ~58% in girls. Sibling interaction effects were only significant in the opposite-sex twin pairs: the interaction effect of males on their female co-twin differed from the effect of females on their male co-twin. An aggressive female had a positive effect on male co-twin aggression, whereas more aggression in males had a negative influence on a female co-twin.

CONCLUSIONS

Opposite-sex twins displayed unique social dynamics of aggressive behaviors in a joint analysis of a large, multinational dataset. The integrative data analysis framework, applied in combination with a reference panel, has the potential to elucidate broad, generalizable results in the investigation of common psychological traits in children.

Aggressive behaviour in childhood and adolescence: the role of smoking during pregnancy, evidence from four twin cohorts in the EU-ACTION consortium

BACKGROUND

Maternal smoking during pregnancy (MSDP) has been linked to offspring's externalizing problems. It has been argued that socio-demographic factors (e.g. maternal age and education), co-occurring environmental risk factors, or pleiotropic genetic effects may account for the association between MSDP and later outcomes. This study provides a comprehensive investigation of the association between MSDP and a single harmonized component of externalizing: aggressive behaviour, measured throughout childhood and adolescence.

METHODS

Data came from four prospective twin cohorts - Twins Early Development Study, Netherlands Twin Register, Childhood and Adolescent Twin Study of Sweden, and FinnTwin12 study - who collaborate in the EU-ACTION consortium. Data from 30 708 unrelated individuals were analysed. Based on item level data, a harmonized measure of aggression was created at ages 9-10; 12; 14-15 and 16-18.

RESULTS

MSDP predicted aggression in childhood and adolescence. A meta-analysis across the four samples found the independent effect of MSDP to be 0.4% (r = 0.066), this remained consistent when analyses were performed separately by sex. All other perinatal factors combined explained 1.1% of the variance in aggression across all ages and samples (r = 0.112). Paternal smoking and aggressive parenting strategies did not account for the MSDP-aggression association, consistent with the hypothesis of a small direct link between MSDP and aggression.

CONCLUSIONS

Perinatal factors, including MSDP, account for a small portion of the variance in aggression in childhood and adolescence. Later experiences may play a greater role in shaping adolescents' aggressive behaviour.

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.

Genome-wide association meta-analysis in 269,867 individuals identifies new genetic and functional links to intelligence

Intelligence is highly heritable1 and a major determinant of human health and well-being2. Recent genome-wide meta-analyses have identified 24 genomic loci linked to variation in intelligence3-7, but much about its genetic underpinnings remains to be discovered. Here, we present a large-scale genetic association study of intelligence (n = 269,867), identifying 205 associated genomic loci (190 new) and 1,016 genes (939 new) via positional mapping, expression quantitative trait locus (eQTL) mapping, chromatin interaction mapping, and gene-based association analysis. We find enrichment of genetic effects in conserved and coding regions and associations with 146 nonsynonymous exonic variants. Associated genes are strongly expressed in the brain, specifically in striatal medium spiny neurons and hippocampal pyramidal neurons. Gene set analyses implicate pathways related to nervous system development and synaptic structure. We confirm previous strong genetic correlations with multiple health-related outcomes, and Mendelian randomization analysis results suggest protective effects of intelligence for Alzheimer's disease and ADHD and bidirectional causation with pleiotropic effects for schizophrenia. These results are a major step forward in understanding the neurobiology of cognitive function as well as genetically related neurological and psychiatric disorders.

Genetics of co-developing conduct and emotional problems during childhood and adolescence

Common genetic influences offer a partial explanation for comorbidity between different psychiatric disorders1-3. However, the genetics underlying co-development-the cross-domain co-occurrence of patterns of change over time-of psychiatric symptoms during childhood and adolescence has not been well explored. Here, we show genetic influence on joint symptom trajectories of parent-reported conduct and emotional problems (overall N = 15,082) across development (4-16 years) using both twin- and genome-wide polygenic score analyses (genotyped N = 2,610). Specifically, we found seven joint symptom trajectories, including two characterized by jointly stable and jointly increasing symptoms of conduct and emotional problems, respectively (7.3% of the sample, collectively). Twin modelling analyses revealed substantial genetic influence on trajectories (heritability estimates range of 0.41-0.78). Furthermore, individuals' risk of being classified in the most symptomatic trajectory classes was significantly predicted by polygenic scores for years-of-education-associated alleles and depressive symptoms-associated alleles. Complementary analyses of child self-reported symptoms across late childhood and early adolescence yielded broadly similar results. Taken together, our results indicate that genetic factors are involved in the co-development of conduct and emotional problems across childhood and adolescence, and that individuals with co-developing symptoms across multiple domains may represent a clinical subgroup characterized by increased levels of genetic risk.

A genome-wide association study for extremely high intelligence

We used a case-control genome-wide association (GWA) design with cases consisting of 1238 individuals from the top 0.0003 (~170 mean IQ) of the population distribution of intelligence and 8172 unselected population-based controls. The single-nucleotide polymorphism heritability for the extreme IQ trait was 0.33 (0.02), which is the highest so far for a cognitive phenotype, and significant genome-wide genetic correlations of 0.78 were observed with educational attainment and 0.86 with population IQ. Three variants in locus ADAM12 achieved genome-wide significance, although they did not replicate with published GWA analyses of normal-range IQ or educational attainment. A genome-wide polygenic score constructed from the GWA results accounted for 1.6% of the variance of intelligence in the normal range in an unselected sample of 3414 individuals, which is comparable to the variance explained by GWA studies of intelligence with substantially larger sample sizes. The gene family plexins, members of which are mutated in several monogenic neurodevelopmental disorders, was significantly enriched for associations with high IQ. This study shows the utility of extreme trait selection for genetic study of intelligence and suggests that extremely high intelligence is continuous genetically with normal-range intelligence in the population.

Multi-polygenic score approach to trait prediction

A primary goal of polygenic scores, which aggregate the effects of thousands of trait-associated DNA variants discovered in genome-wide association studies (GWASs), is to estimate individual-specific genetic propensities and predict outcomes. This is typically achieved using a single polygenic score, but here we use a multi-polygenic score (MPS) approach to increase predictive power by exploiting the joint power of multiple discovery GWASs, without assumptions about the relationships among predictors. We used summary statistics of 81 well-powered GWASs of cognitive, medical and anthropometric traits to predict three core developmental outcomes in our independent target sample: educational achievement, body mass index (BMI) and general cognitive ability. We used regularized regression with repeated cross-validation to select from and estimate contributions of 81 polygenic scores in a UK representative sample of 6710 unrelated adolescents. The MPS approach predicted 10.9% variance in educational achievement, 4.8% in general cognitive ability and 5.4% in BMI in an independent test set, predicting 1.1%, 1.1%, and 1.6% more variance than the best single-score predictions. As other relevant GWA analyses are reported, they can be incorporated in MPS models to maximize phenotype prediction. The MPS approach should be useful in research with modest sample sizes to investigate developmental, multivariate and gene-environment interplay issues and, eventually, in clinical settings to predict and prevent problems using personalized interventions.

Individual and shared effects of social environment and polygenic risk scores on adolescent body mass index

Juvenile obesity is associated with adverse health outcomes. Understanding genetic and environmental influences on body mass index (BMI) during adolescence could inform interventions. We investigated independent and interactive effects of parenting, socioeconomic status (SES) and polygenic risk on BMI pre-adolescence, and on the rate of change in BMI across adolescence. Genome-wide genotype data, BMI and child perceptions of parental warmth and punitive discipline were available at 11 years old, and parental SES was available from birth on 3,414 unrelated participants. Linear models were used to test the effects of social environment and polygenic risk on pre-adolescent BMI. Change in BMI across adolescence was assessed in a subset (N = 1943). Sex-specific effects were assessed. Higher genetic risk was associated with increased BMI pre-adolescence and across adolescence (p < 0.00417, corrected for multiple tests). Negative parenting was not significantly associated with either phenotype, but lower SES was associated with increased BMI pre-adolescence. No interactions passed correction for multiple testing. Polygenic risk scores from adult GWAS meta-analyses are associated with BMI in juveniles, suggesting a stable genetic component. Pre-adolescent BMI was associated with social environment, but parental style has, at most, a small effect.

Genetic influence on social outcomes during and after the Soviet era in Estonia

The etiology of individual differences in educational attainment and occupational status includes genetic as well as environmental factors1-5 and can change as societies change3,6,7. The extent of genetic influence on these social outcomes can be viewed as an index of success in achieving meritocratic values of equality of opportunity by rewarding talent and hard work, which are to a large extent influenced by genetic factors, rather than rewarding environmentally driven privilege. To the extent that the end of the Soviet Union and the independence of Estonia led to an increase in meritocratic selection of individuals in education and occupation, genetic influence should be higher in the post-Soviet era than in the Soviet era. Here we confirmed this hypothesis: DNA differences (single-nucleotide polymorphisms, SNPs) explained twice as much variance in educational attainment and occupational status in the post-Soviet era compared to the Soviet era in both polygenic score analyses and SNP heritability analyses of 12 500 Estonians. This is the first demonstration of a change in the extent of genetic influence in the same population following a massive and abrupt social change - in this case, the shift from a communist to a capitalist society.

Differences in exam performance between pupils attending selective and non-selective schools mirror the genetic differences between them

On average, students attending selective schools outperform their non-selective counterparts in national exams. These differences are often attributed to value added by the school, as well as factors schools use to select pupils, including ability, achievement and, in cases where schools charge tuition fees or are located in affluent areas, socioeconomic status. However, the possible role of DNA differences between students of different schools types has not yet been considered. We used a UK-representative sample of 4814 genotyped students to investigate exam performance at age 16 and genetic differences between students in three school types: state-funded, non-selective schools ('non-selective'), state-funded, selective schools ('grammar') and private schools, which are selective ('private'). We created a genome-wide polygenic score (GPS) derived from a genome-wide association study of years of education (EduYears). We found substantial mean genetic differences between students of different school types: students in non-selective schools had lower EduYears GPS compared to those in grammar (d = 0.41) and private schools (d = 0.37). Three times as many students in the top EduYears GPS decile went to a selective school compared to the bottom decile. These results were mirrored in the exam differences between school types. However, once we controlled for factors involved in pupil selection, there were no significant genetic differences between school types, and the variance in exam scores at age 16 explained by school type dropped from 7% to <1%. These results show that genetic and exam differences between school types are primarily due to the heritable characteristics involved in pupil admission.

Predicting educational achievement from DNA

This corrects the article DOI: 10.1038/mp.2016.107.

Genetics of depressive symptoms in adolescence

BACKGROUND

Despite many attempts to understand the genetic architecture of depression, little progress has been made. The majority of these studies, however, have been carried out in adults and do not account for the potential influence of development.

METHODS

The Avon Longitudinal Study of Parents and Children (ALSPAC) is a longitudinal pregnancy cohort which recruited participants between April 1991 and December 1992. Analyses were replicated in two independent European cohorts. Genome-wide complex trait analysis (GCTA) software was used to investigate SNP-heritability (h2SNP) of depression across adolescence, the role of puberty was investigated by stratifying these estimates according to pubertal onset. Genome-wide association studies were performed to identify genetic variants associated with depression at different stages of development.

RESULTS

Heritability was estimated between the ages of 11 and 18 with sample sizes ranging from 3289 to 5480. Heritability was low with an apparent peak was found at age 13 (h2 = 0.17, p = 0.006). Confidence intervals around these estimates suggest an upper-bound to h2SNP of around 30%. A variant located on chromosome 7 was found to be associated with depressive symptoms at age 13 in ALSPAC (rs138191010: β = 0.142, p = 2.51 × 10-8), although this was not replicated.

CONCLUSIONS

Although power is a potential limitation, the observed patterns provide interesting hypotheses surrounding the heritability of depression at different developmental stages. We found substantially lower estimates for depressive symptoms at age 11 (0.07) compared to those previously estimated in adults (0.21). We also found a peak in heritability at age 13. These findings suggest environmental factors are likely to be more important in the aetiology of depressive symptoms in early adolescence than in adulthood.

Genetic Influence on Intergenerational Educational Attainment

Using twin (6,105 twin pairs) and genomic (5,825 unrelated individuals taken from the twin sample) analyses, we tested for genetic influences on the parent-offspring correspondence in educational attainment. Genetics accounted for nearly half of the variance in intergenerational educational attainment. A genomewide polygenic score (GPS) for years of education was also associated with intergenerational educational attainment: The highest and lowest GPS means were found for offspring in stably educated families (i.e., who had taken A Levels and had a university-educated parent; M = 0.43, SD = 0.97) and stably uneducated families (i.e., who had not taken A Levels and had no university-educated parent; M = -0.19, SD = 0.97). The average GPSs fell in between for children who were upwardly mobile (i.e., who had taken A Levels but had no university-educated parent; M = 0.05, SD = 0.96) and children who were downwardly mobile (i.e., who had not taken A Levels but had a university-educated parent; M = 0.28, SD = 1.03). Genetic influences on intergenerational educational attainment can be viewed as an index of equality of educational opportunity.

Genome-Wide Polygenic Scores Predict Reading Performance Throughout the School Years

It is now possible to create individual-specific genetic scores, called genome-wide polygenic scores (GPS). We used a GPS for years of education (EduYears) to predict reading performance assessed at UK National Curriculum Key Stages 1 (age 7), 2 (age 12) and 3 (age 14) and on reading tests administered at ages 7 and 12 in a UK sample of 5,825 unrelated individuals. EduYears GPS accounts for up to 5% of the variance in reading performance at age 14. GPS predictions remained significant after accounting for general cognitive ability and family socioeconomic status. Reading performance of children in the lowest and highest 12.5% of the EduYears GPS distribution differed by a mean growth in reading ability of approximately two school years. It seems certain that polygenic scores will be used to predict strengths and weaknesses in education.

Paternal Age Alters Social Development in Offspring

OBJECTIVE

Advanced paternal age (APA) at conception has been linked with autism and schizophrenia in offspring, neurodevelopmental disorders that affect social functioning. The current study explored the effects of paternal age on social development in the general population.

METHOD

We used multilevel growth modeling to investigate APA effects on socioemotional development from early childhood until adolescence, as measured by the Strengths and Difficulties Questionnaire (SDQ) in the Twins Early Development Study (TEDS) sample. We also investigated genetic and environmental underpinnings of the paternal age effects on development, using the Additive genetics, Common environment, unique Environment (ACE) and gene-environment (GxE) models.

RESULTS

In the general population, both very young and advanced paternal ages were associated with altered trajectory of social development (intercept: p = .01; slope: p = .03). No other behavioral domain was affected by either young or advanced age at fatherhood, suggesting specificity of paternal age effects. Increased importance of genetic factors in social development was recorded in the offspring of older but not very young fathers, suggesting distinct underpinnings of the paternal age effects at these two extremes.

CONCLUSION

Our findings highlight that the APA-related deficits that lead to autism and schizophrenia are likely continuously distributed in the population.

Predicting educational achievement from DNA

A genome-wide polygenic score (GPS), derived from a 2013 genome-wide association study (N=127,000), explained 2% of the variance in total years of education (EduYears). In a follow-up study (N=329,000), a new EduYears GPS explains up to 4%. Here, we tested the association between this latest EduYears GPS and educational achievement scores at ages 7, 12 and 16 in an independent sample of 5825 UK individuals. We found that EduYears GPS explained greater amounts of variance in educational achievement over time, up to 9% at age 16, accounting for 15% of the heritable variance. This is the strongest GPS prediction to date for quantitative behavioral traits. Individuals in the highest and lowest GPS septiles differed by a whole school grade at age 16. Furthermore, EduYears GPS was associated with general cognitive ability (~3.5%) and family socioeconomic status (~7%). There was no evidence of an interaction between EduYears GPS and family socioeconomic status on educational achievement or on general cognitive ability. These results are a harbinger of future widespread use of GPS to predict genetic risk and resilience in the social and behavioral sciences.

Fine mapping genetic associations between the HLA region and extremely high intelligence

General cognitive ability (intelligence) is one of the most heritable behavioural traits and most predictive of socially important outcomes and health. We hypothesized that some of the missing heritability of IQ might lie hidden in the human leukocyte antigen (HLA) region, which plays a critical role in many diseases and traits but is not well tagged in conventional GWAS. Using a uniquely powered design, we investigated whether fine-mapping of the HLA region could narrow the missing heritability gap. Our case-control design included 1,393 cases with extremely high intelligence scores (top 0.0003 of the population equivalent to IQ > 147) and 3,253 unselected population controls. We imputed variants in 200 genes across the HLA region, one SNP (rs444921) reached our criterion for study-wide significance. SNP-based heritability of the HLA variants was small and not significant (h2 = 0.3%, SE = 0.2%). A polygenic score from the case-control genetic association analysis of SNPs in the HLA region did not significantly predict individual differences in intelligence in an independent unselected sample. We conclude that although genetic variation in the HLA region is important to the aetiology of many disorders, it does not appear to be hiding much of the missing heritability of intelligence.

A Genome-Wide Association Meta-Analysis of Attention-Deficit/Hyperactivity Disorder Symptoms in Population-Based Pediatric Cohorts

OBJECTIVE

The aims of this study were to elucidate the influence of common genetic variants on childhood attention-deficit/hyperactivity disorder (ADHD) symptoms, to identify genetic variants that explain its high heritability, and to investigate the genetic overlap of ADHD symptom scores with ADHD diagnosis.

METHOD

Within the EArly Genetics and Lifecourse Epidemiology (EAGLE) consortium, genome-wide single nucleotide polymorphisms (SNPs) and ADHD symptom scores were available for 17,666 children (<13 years of age) from nine population-based cohorts. SNP-based heritability was estimated in data from the three largest cohorts. Meta-analysis based on genome-wide association (GWA) analyses with SNPs was followed by gene-based association tests, and the overlap in results with a meta-analysis in the Psychiatric Genomics Consortium (PGC) case-control ADHD study was investigated.

RESULTS

SNP-based heritability ranged from 5% to 34%, indicating that variation in common genetic variants influences ADHD symptom scores. The meta-analysis did not detect genome-wide significant SNPs, but three genes, lying close to each other with SNPs in high linkage disequilibrium (LD), showed a gene-wide significant association (p values between 1.46 × 10(-6) and 2.66 × 10(-6)). One gene, WASL, is involved in neuronal development. Both SNP- and gene-based analyses indicated overlap with the PGC meta-analysis results with the genetic correlation estimated at 0.96.

CONCLUSION

The SNP-based heritability for ADHD symptom scores indicates a polygenic architecture, and genes involved in neurite outgrowth are possibly involved. Continuous and dichotomous measures of ADHD appear to assess a genetically common phenotype. A next step is to combine data from population-based and case-control cohorts in genetic association studies to increase sample size and to improve statistical power for identifying genetic variants.

Phenome-wide analysis of genome-wide polygenic scores

Genome-wide polygenic scores (GPS), which aggregate the effects of thousands of DNA variants from genome-wide association studies (GWAS), have the potential to make genetic predictions for individuals. We conducted a systematic investigation of associations between GPS and many behavioral traits, the behavioral phenome. For 3152 unrelated 16-year-old individuals representative of the United Kingdom, we created 13 GPS from the largest GWAS for psychiatric disorders (for example, schizophrenia, depression and dementia) and cognitive traits (for example, intelligence, educational attainment and intracranial volume). The behavioral phenome included 50 traits from the domains of psychopathology, personality, cognitive abilities and educational achievement. We examined phenome-wide profiles of associations for the entire distribution of each GPS and for the extremes of the GPS distributions. The cognitive GPS yielded stronger predictive power than the psychiatric GPS in our UK-representative sample of adolescents. For example, education GPS explained variation in adolescents' behavior problems (~0.6%) and in educational achievement (~2%) but psychiatric GPS were associated with neither. Despite the modest effect sizes of current GPS, quantile analyses illustrate the ability to stratify individuals by GPS and opportunities for research. For example, the highest and lowest septiles for the education GPS yielded a 0.5 s.d. difference in mean math grade and a 0.25 s.d. difference in mean behavior problems. We discuss the usefulness and limitations of GPS based on adult GWAS to predict genetic propensities earlier in development.

Longitudinal heritability of childhood aggression

The genetic and environmental contributions to the variation and longitudinal stability in childhood aggressive behavior were assessed in two large twin cohorts, the Netherlands Twin Register (NTR), and the Twins Early Development Study (TEDS; United Kingdom). In NTR, maternal ratings on aggression from the Child Behavior Checklist (CBCL) were available for 10,765 twin pairs at age 7, for 8,557 twin pairs at age 9/10, and for 7,176 twin pairs at age 12. In TEDS, parental ratings of conduct disorder from the Strength and Difficulty Questionnaire (SDQ) were available for 6,897 twin pairs at age 7, for 3,028 twin pairs at age 9 and for 5,716 twin pairs at age 12. In both studies, stability and heritability of aggressive behavioral problems was high. Heritability was on average somewhat, but significantly, lower in TEDS (around 60%) than in NTR (between 50% and 80%) and sex differences were slightly larger in the NTR sample. In both studies, the influence of shared environment was similar: in boys shared environment explained around 20% of the variation in aggression across all ages while in girls its influence was absent around age 7 and only came into play at later ages. Longitudinal genetic correlations were the main reason for stability of aggressive behavior. Individual differences in CBCL-Aggressive Behavior and SDQ-Conduct disorder throughout childhood are driven by a comparable but significantly different genetic architecture. © 2016 Wiley Periodicals, Inc.

Childhood gene-environment interactions and age-dependent effects of genetic variants associated with refractive error and myopia: The CREAM Consortium

Myopia, currently at epidemic levels in East Asia, is a leading cause of untreatable visual impairment. Genome-wide association studies (GWAS) in adults have identified 39 loci associated with refractive error and myopia. Here, the age-of-onset of association between genetic variants at these 39 loci and refractive error was investigated in 5200 children assessed longitudinally across ages 7-15 years, along with gene-environment interactions involving the major environmental risk-factors, nearwork and time outdoors. Specific variants could be categorized as showing evidence of: (a) early-onset effects remaining stable through childhood, (b) early-onset effects that progressed further with increasing age, or (c) onset later in childhood (N = 10, 5 and 11 variants, respectively). A genetic risk score (GRS) for all 39 variants explained 0.6% (P = 6.6E-08) and 2.3% (P = 6.9E-21) of the variance in refractive error at ages 7 and 15, respectively, supporting increased effects from these genetic variants at older ages. Replication in multi-ancestry samples (combined N = 5599) yielded evidence of childhood onset for 6 of 12 variants present in both Asians and Europeans. There was no indication that variant or GRS effects altered depending on time outdoors, however 5 variants showed nominal evidence of interactions with nearwork (top variant, rs7829127 in ZMAT4; P = 6.3E-04).

Assortative Mating-A Missing Piece in the Jigsaw of Psychiatric Genetics

The topic of assortative (nonrandom) mating might seem esoteric or even salacious. For example, in lectures you have to point out to students that random mating is not about promiscuity. In this issue of JAMA Psychiatry, Nordsletten and colleagues report the first general population study to date of assortative mating for psychiatric disorders, which may help to solve 3 puzzles in psychiatric genetics: Why are psychiatric disorders so highly heritable when they are associated with reduced fecundity? Why are some psychiatric disorders so much more highly heritable than others? Why is there so much genetic comorbidity across psychiatric disorders?

Genetic link between family socioeconomic status and children's educational achievement estimated from genome-wide SNPs

One of the best predictors of children's educational achievement is their family's socioeconomic status (SES), but the degree to which this association is genetically mediated remains unclear. For 3000 UK-representative unrelated children we found that genome-wide single-nucleotide polymorphisms could explain a third of the variance of scores on an age-16 UK national examination of educational achievement and half of the correlation between their scores and family SES. Moreover, genome-wide polygenic scores based on a previously published genome-wide association meta-analysis of total number of years in education accounted for ~3.0% variance in educational achievement and ~2.5% in family SES. This study provides the first molecular evidence for substantial genetic influence on differences in children's educational achievement and its association with family SES.

Thinking positively: The genetics of high intelligence

High intelligence (general cognitive ability) is fundamental to the human capital that drives societies in the information age. Understanding the origins of this intellectual capital is important for government policy, for neuroscience, and for genetics. For genetics, a key question is whether the genetic causes of high intelligence are qualitatively or quantitatively different from the normal distribution of intelligence. We report results from a sibling and twin study of high intelligence and its links with the normal distribution. We identified 360,000 sibling pairs and 9000 twin pairs from 3 million 18-year-old males with cognitive assessments administered as part of conscription to military service in Sweden between 1968 and 2010. We found that high intelligence is familial, heritable, and caused by the same genetic and environmental factors responsible for the normal distribution of intelligence. High intelligence is a good candidate for "positive genetics" - going beyond the negative effects of DNA sequence variation on disease and disorders to consider the positive end of the distribution of genetic effects.

Genetic influence on family socioeconomic status and children's intelligence

Environmental measures used widely in the behavioral sciences show nearly as much genetic influence as behavioral measures, a critical finding for interpreting associations between environmental factors and children's development. This research depends on the twin method that compares monozygotic and dizygotic twins, but key aspects of children's environment such as socioeconomic status (SES) cannot be investigated in twin studies because they are the same for children growing up together in a family. Here, using a new technique applied to DNA from 3000 unrelated children, we show significant genetic influence on family SES, and on its association with children's IQ at ages 7 and 12. In addition to demonstrating the ability to investigate genetic influence on between-family environmental measures, our results emphasize the need to consider genetics in research and policy on family SES and its association with children's IQ.

Strong genetic influence on a UK nationwide test of educational achievement at the end of compulsory education at age 16

We have previously shown that individual differences in educational achievement are highly heritable in the early and middle school years in the UK. The objective of the present study was to investigate whether similarly high heritability is found at the end of compulsory education (age 16) for the UK-wide examination, called the General Certificate of Secondary Education (GCSE). In a national twin sample of 11,117 16-year-olds, heritability was substantial for overall GCSE performance for compulsory core subjects (58%) as well as for each of them individually: English (52%), mathematics (55%) and science (58%). In contrast, the overall effects of shared environment, which includes all family and school influences shared by members of twin pairs growing up in the same family and attending the same school, accounts for about 36% of the variance of mean GCSE scores. The significance of these findings is that individual differences in educational achievement at the end of compulsory education are not primarily an index of the quality of teachers or schools: much more of the variance of GCSE scores can be attributed to genetics than to school or family environment. We suggest a model of education that recognizes the important role of genetics. Rather than a passive model of schooling as instruction (instruere, ‘to build in’), we propose an active model of education (educare, ‘to bring out’) in which children create their own educational experiences in part on the basis of their genetic propensities, which supports the trend towards personalized learning.