Exploring polygenic-environment and residual-environment interactions for depressive symptoms within the UK Biobank

Substantial advances have been made in identifying genetic contributions to depression, but little is known about how the effect of genes can be modulated by the environment, creating a gene-environment interaction. Using multivariate reaction norm models (MRNMs) within the UK Biobank (N = 61294-91644), we investigate whether the polygenic and residual variance components of depressive symptoms are modulated by 17 a priori selected covariate traits-12 environmental variables and 5 biomarkers. MRNMs, a mixed-effects modelling approach, provide unbiased polygenic-covariate interaction estimates for a quantitative trait by controlling for outcome-covariate correlations and residual-covariate interactions. A continuous depressive symptom variable was the outcome in 17 MRNMs-one for each covariate trait. Each MRNM had a fixed-effects model (fixed effects included the covariate trait, demographic variables, and principal components) and a random effects model (where polygenic-covariate and residual-covariate interactions are modelled). Of the 17 selected covariates, 11 significantly modulate deviations in depressive symptoms through the modelled interactions, but no single interaction explains a large proportion of phenotypic variation. Results are dominated by residual-covariate interactions, suggesting that covariate traits (including neuroticism, childhood trauma, and BMI) typically interact with unmodelled variables, rather than a genome-wide polygenic component, to influence depressive symptoms. Only average sleep duration has a polygenic-covariate interaction explaining a demonstrably nonzero proportion of the variability in depressive symptoms. This effect is small, accounting for only 1.22% (95% confidence interval: [0.54, 1.89]) of variation. The presence of an interaction highlights a specific focus for intervention, but the negative results here indicate a limited contribution from polygenic-environment interactions.

The genetic case for cardiorespiratory fitness as a clinical vital sign and the routine prescription of physical activity in healthcare

BACKGROUND

Cardiorespiratory fitness (CRF) and physical activity (PA) are well-established predictors of morbidity and all-cause mortality. However, CRF is not routinely measured and PA not routinely prescribed as part of standard healthcare. The American Heart Association (AHA) recently presented a scientific case for the inclusion of CRF as a clinical vital sign based on epidemiological and clinical observation. Here, we leverage genetic data in the UK Biobank (UKB) to strengthen the case for CRF as a vital sign and make a case for the prescription of PA.

METHODS

We derived two CRF measures from the heart rate data collected during a submaximal cycle ramp test: CRF-vo2max, an estimate of the participants' maximum volume of oxygen uptake, per kilogram of body weight, per minute; and CRF-slope, an estimate of the rate of increase of heart rate during exercise. Average PA over a 7-day period was derived from a wrist-worn activity tracker. After quality control, 70,783 participants had data on the two derived CRF measures, and 89,683 had PA data. We performed genome-wide association study (GWAS) analyses by sex, and post-GWAS techniques to understand genetic architecture of the traits and prioritise functional genes for follow-up.

RESULTS

We found strong evidence that genetic variants associated with CRF and PA influenced genetic expression in a relatively small set of genes in the heart, artery, lung, skeletal muscle and adipose tissue. These functionally relevant genes were enriched among genes known to be associated with coronary artery disease (CAD), type 2 diabetes (T2D) and Alzheimer's disease (three of the top 10 causes of death in high-income countries) as well as Parkinson's disease, pulmonary fibrosis, and blood pressure, heart rate, and respiratory phenotypes. Genetic variation associated with lower CRF and PA was also correlated with several disease risk factors (including greater body mass index, body fat and multiple obesity phenotypes); a typical T2D profile (including higher insulin resistance, higher fasting glucose, impaired beta-cell function, hyperglycaemia, hypertriglyceridemia); increased risk for CAD and T2D; and a shorter lifespan.

CONCLUSIONS

Genetics supports three decades of evidence for the inclusion of CRF as a clinical vital sign. Given the genetic, clinical and epidemiological evidence linking CRF and PA to increased morbidity and mortality, regular measurement of CRF as a marker of health and routine prescription of PA could be a prudent strategy to support public health.

Genetic and clinical characteristics of treatment-resistant depression using primary care records in two UK cohorts

Treatment-resistant depression (TRD) is a major contributor to the disability caused by major depressive disorder (MDD). Primary care electronic health records provide an easily accessible approach to investigate TRD clinical and genetic characteristics. MDD defined from primary care records in UK Biobank (UKB) and EXCEED studies was compared with other measures of depression and tested for association with MDD polygenic risk score (PRS). Using prescribing records, TRD was defined from at least two switches between antidepressant drugs, each prescribed for at least 6 weeks. Clinical-demographic characteristics, SNP-based heritability (h²_SNP) and genetic overlap with psychiatric and non-psychiatric traits were compared in TRD and non-TRD MDD cases. In 230,096 and 8926 UKB and EXCEED participants with primary care data, respectively, the prevalence of MDD was 8.7% and 14.2%, of which 13.2% and 13.5% was TRD, respectively. In both cohorts, MDD defined from primary care records was strongly associated with MDD PRS, and in UKB it showed overlap of 71-88% with other MDD definitions. In UKB, TRD vs healthy controls and non-TRD vs healthy controls h²_SNP was comparable (0.25 [SE = 0.04] and 0.19 [SE = 0.02], respectively). TRD vs non-TRD was positively associated with the PRS of attention deficit hyperactivity disorder, with lower socio-economic status, obesity, higher neuroticism and other unfavourable clinical characteristics. This study demonstrated that MDD and TRD can be reliably defined using primary care records and provides the first large scale population assessment of the genetic, clinical and demographic characteristics of TRD.

Multiple measures of depression to enhance validity of major depressive disorder in the UK Biobank

BACKGROUND

The UK Biobank contains data with varying degrees of reliability and completeness for assessing depression. A third of participants completed a Mental Health Questionnaire (MHQ) containing the gold-standard Composite International Diagnostic Interview (CIDI) criteria for assessing mental health disorders.

AIMS

To investigate whether multiple observations of depression from sources other than the MHQ can enhance the validity of major depressive disorder (MDD).

METHOD

In participants who did not complete the MHQ, we calculated the number of other depression measures endorsed, for example from hospital episode statistics and interview data. We compared cases defined this way with CIDI-defined cases for several estimates: the variance explained by polygenic risk scores (PRS), area under the curve attributable to PRS, single nucleotide polymorphisms (SNPs)-based heritability and genetic correlations with summary statistics from the Psychiatric Genomics Consortium MDD genome-wide association study.

RESULTS

The strength of the genetic contribution increased with the number of measures endorsed. For example, SNP-based heritability increased from 7% in participants who endorsed only one measure of depression, to 21% in those who endorsed four or five measures of depression. The strength of the genetic contribution to cases defined by at least two measures approximated that for CIDI-defined cases. Most genetic correlations between UK Biobank and the Psychiatric Genomics Consortium MDD study exceeded 0.7, but there was variability between pairwise comparisons.

CONCLUSIONS

Multiple measures of depression can serve as a reliable approximation for case status where the CIDI measure is not available, indicating sample size can be optimised using the entire suite of UK Biobank data.

Genome-wide association study of MRI markers of cerebral small vessel disease in 42,310 participants

Cerebral small vessel disease is a major cause of stroke and dementia, but its genetic basis is incompletely understood. We perform a genetic study of three MRI markers of the disease in UK Biobank imaging data and other sources: white matter hyperintensities (N = 42,310), fractional anisotropy (N = 17,663) and mean diffusivity (N = 17,467). Our aim is to better understand the disease pathophysiology. Across the three traits, we identify 31 loci, of which 21 were previously unreported. We perform a transcriptome-wide association study to identify associations with gene expression in relevant tissues, identifying 66 associated genes across the three traits. This genetic study provides insights into the understanding of the biological mechanisms underlying small vessel disease.

Genetic and Inflammatory Biomarkers Classify Small Intestine Inflammation in Asymptomatic First-degree Relatives of Patients With Crohn's Disease

BACKGROUND & AIMS

Relatives of individuals with Crohn's disease (CD) carry CD-associated genetic variants and are often exposed to environmental factors that increase their risk for this disease. We aimed to estimate the utility of genotype, smoking status, family history, and biomarkers can calculate risk in asymptomatic first-degree relatives of patients with CD.

METHODS

We recruited 480 healthy first-degree relatives (full siblings, offspring or parents) of patients with CD through the Guy's and St Thomas' NHS Foundation Trust and from members of Crohn's and Colitis, United Kingdom. DNA samples were genotyped using the Immunochip. We calculated a risk score for 454 participants, based on 72 genetic variants associated with CD, family history, and smoking history. Participants were assigned to highest and lowest risk score quartiles. We assessed pre-symptomatic inflammation by capsule endoscopy and measured 22 markers of inflammation in stool and serum samples (reference standard). Two machine-learning classifiers (elastic net and random forest) were used to assess the ability of the risk factors and biomarkers to identify participants with small intestinal inflammation in the same dataset.

RESULTS

The machine-learning classifiers identified participants with pre-symptomatic intestinal inflammation: elastic net (area under the curve, 0.80; 95% CI, 0.62-0.98) and random forest (area under the curve, 0.87; 95% CI, 0.75-1.00). The elastic net method identified 3 variables that can be used to calculate odds for intestinal inflammation: combined family history of CD (odds ratio, 1.31), genetic risk score (odds ratio, 1.14), and fecal calprotectin (odds ratio, 1.04). These same 3 variables were among the 5 factors associated with intestinal inflammation in the random forest model.

CONCLUSION

Using machine learning classifiers, we found that genetic variants associated with CD, family history, and fecal calprotectin together identify individuals with pre-symptomatic intestinal inflammation who are therefore at risk for CD. A tool for detecting people at risk for CD before they develop symptoms would help identify the individuals most likely to benefit from early intervention.

Classical Human Leukocyte Antigen Alleles and C4 Haplotypes Are Not Significantly Associated With Depression

BACKGROUND

The prevalence of depression is higher in individuals with autoimmune diseases, but the mechanisms underlying the observed comorbidities are unknown. Shared genetic etiology is a plausible explanation for the overlap, and in this study we tested whether genetic variation in the major histocompatibility complex (MHC), which is associated with risk for autoimmune diseases, is also associated with risk for depression.

METHODS

We fine-mapped the classical MHC (chr6: 29.6-33.1 Mb), imputing 216 human leukocyte antigen (HLA) alleles and 4 complement component 4 (C4) haplotypes in studies from the Psychiatric Genomics Consortium Major Depressive Disorder Working Group and the UK Biobank. The total sample size was 45,149 depression cases and 86,698 controls. We tested for association between depression status and imputed MHC variants, applying both a region-wide significance threshold (3.9 × 10-6) and a candidate threshold (1.6 × 10-4).

RESULTS

No HLA alleles or C4 haplotypes were associated with depression at the region-wide threshold. HLA-B*08:01 was associated with modest protection for depression at the candidate threshold for testing in HLA genes in the meta-analysis (odds ratio = 0.98, 95% confidence interval = 0.97-0.99).

CONCLUSIONS

We found no evidence that an increased risk for depression was conferred by HLA alleles, which play a major role in the genetic susceptibility to autoimmune diseases, or C4 haplotypes, which are strongly associated with schizophrenia. These results suggest that any HLA or C4 variants associated with depression either are rare or have very modest effect sizes.

Genetic correlations of psychiatric traits with body composition and glycemic traits are sex- and age-dependent

Body composition is often altered in psychiatric disorders. Using genome-wide common genetic variation data, we calculate sex-specific genetic correlations amongst body fat %, fat mass, fat-free mass, physical activity, glycemic traits and 17 psychiatric traits (up to N = 217,568). Two patterns emerge: (1) anorexia nervosa, schizophrenia, obsessive-compulsive disorder, and education years are negatively genetically correlated with body fat % and fat-free mass, whereas (2) attention-deficit/hyperactivity disorder (ADHD), alcohol dependence, insomnia, and heavy smoking are positively correlated. Anorexia nervosa shows a stronger genetic correlation with body fat % in females, whereas education years is more strongly correlated with fat mass in males. Education years and ADHD show genetic overlap with childhood obesity. Mendelian randomization identifies schizophrenia, anorexia nervosa, and higher education as causal for decreased fat mass, with higher body fat % possibly being a causal risk factor for ADHD and heavy smoking. These results suggest new possibilities for targeted preventive strategies.

Genomics of body fat percentage may contribute to sex bias in anorexia nervosa

Anorexia nervosa (AN) occurs nine times more often in females than in males. Although environmental factors likely play a role, the reasons for this imbalanced sex ratio remain unresolved. AN displays high genetic correlations with anthropometric and metabolic traits. Given sex differences in body composition, we investigated the possible metabolic underpinnings of female propensity for AN. We conducted sex-specific GWAS in a healthy and medication-free subsample of the UK Biobank (n = 155,961), identifying 77 genome-wide significant loci associated with body fat percentage (BF%) and 174 with fat-free mass (FFM). Partitioned heritability analysis showed an enrichment for central nervous tissue-associated genes for BF%, which was more prominent in females than males. Genetic correlations of BF% and FFM with the largest GWAS of AN by the Psychiatric Genomics Consortium were estimated to explore shared genomics. The genetic correlations of BF%male and BF%female with AN differed significantly from each other (p < .0001, δ = -0.17), suggesting that the female preponderance in AN may, in part, be explained by sex-specific anthropometric and metabolic genetic factors increasing liability to AN.

Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa

Characterized primarily by a low body-mass index, anorexia nervosa is a complex and serious illness1, affecting 0.9-4% of women and 0.3% of men2-4, with twin-based heritability estimates of 50-60%5. Mortality rates are higher than those in other psychiatric disorders6, and outcomes are unacceptably poor7. Here we combine data from the Anorexia Nervosa Genetics Initiative (ANGI)8,9 and the Eating Disorders Working Group of the Psychiatric Genomics Consortium (PGC-ED) and conduct a genome-wide association study of 16,992 cases of anorexia nervosa and 55,525 controls, identifying eight significant loci. The genetic architecture of anorexia nervosa mirrors its clinical presentation, showing significant genetic correlations with psychiatric disorders, physical activity, and metabolic (including glycemic), lipid and anthropometric traits, independent of the effects of common variants associated with body-mass index. These results further encourage a reconceptualization of anorexia nervosa as a metabo-psychiatric disorder. Elucidating the metabolic component is a critical direction for future research, and paying attention to both psychiatric and metabolic components may be key to improving outcomes.

ukbtools: An R package to manage and query UK Biobank data

INTRODUCTION

The UK Biobank (UKB) is a resource that includes detailed health-related data on about 500,000 individuals and is available to the research community. However, several obstacles limit immediate analysis of the data: data files vary in format, may be very large, and have numerical codes for column names.

RESULTS

ukbtools removes all the upfront data wrangling required to get a single dataset for statistical analysis. All associated data files are merged into a single dataset with descriptive column names. The package also provides tools to assist in quality control by exploring the primary demographics of subsets of participants; query of disease diagnoses for one or more individuals, and estimating disease frequency relative to a reference variable; and to retrieve genetic metadata.

CONCLUSION

Having a dataset with meaningful variable names, a set of UKB-specific exploratory data analysis tools, disease query functions, and a set of helper functions to explore and write genetic metadata to file, will rapidly enable UKB users to undertake their research.

Genetic fine mapping of systemic lupus erythematosus MHC associations in Europeans and African Americans

Genetic variation within the major histocompatibility complex (MHC) contributes substantial risk for systemic lupus erythematosus, but high gene density, extreme polymorphism and extensive linkage disequilibrium (LD) have made fine mapping challenging. To address the problem, we compared two association techniques in two ancestrally diverse populations, African Americans (AAs) and Europeans (EURs). We observed a greater number of Human Leucocyte Antigen (HLA) alleles in AA consistent with the elevated level of recombination in this population. In EUR we observed 50 different A-C-B-DRB1-DQA-DQB multilocus haplotype sequences per hundred individuals; in the AA sample, these multilocus haplotypes were twice as common compared to Europeans. We also observed a strong narrow class II signal in AA as opposed to the long-range LD observed in EUR that includes class I alleles. We performed a Bayesian model choice of the classical HLA alleles and a frequentist analysis that combined both single nucleotide polymorphisms (SNPs) and classical HLA alleles. Both analyses converged on a similar subset of risk HLA alleles: in EUR HLA- B*08:01 + B*18:01 + (DRB1*15:01 frequentist only) + DQA*01:02 + DQB*02:01 + DRB3*02 and in AA HLA-C*17:01 + B*08:01 + DRB1*15:03 + (DQA*01:02 frequentist only) + DQA*02:01 + DQA*05:01+ DQA*05:05 + DQB*03:19 + DQB*02:02. We observed two additional independent SNP associations in both populations: EUR rs146903072 and rs501480; AA rs389883 and rs114118665. The DR2 serotype was best explained by DRB1*15:03 + DQA*01:02 in AA and by DRB1*15:01 + DQA*01:02 in EUR. The DR3 serotype was best explained by DQA*05:01 in AA and by DQB*02:01 in EUR. Despite some differences in underlying HLA allele risk models in EUR and AA, SNP signals across the extended MHC showed remarkable similarity and significant concordance in direction of effect for risk-associated variants.

Genetic variation at 16q24.2 is associated with small vessel stroke

Objective

Genome‐wide association studies (GWAS) have been successful at identifying associations with stroke and stroke subtypes, but have not yet identified any associations solely with small vessel stroke (SVS). SVS comprises one quarter of all ischemic stroke and is a major manifestation of cerebral small vessel disease, the primary cause of vascular cognitive impairment. Studies across neurological traits have shown that younger‐onset cases have an increased genetic burden. We leveraged this increased genetic burden by performing an age‐at‐onset informed GWAS meta‐analysis, including a large younger‐onset SVS population, to identify novel associations with stroke.

Methods

We used a three‐stage age‐at‐onset informed GWAS to identify novel genetic variants associated with stroke. On identifying a novel locus associated with SVS, we assessed its influence on other small vessel disease phenotypes, as well as on messenger RNA (mRNA) expression of nearby genes, and on DNA methylation of nearby CpG sites in whole blood and in the fetal brain.

Results

We identified an association with SVS in 4,203 cases and 50,728 controls on chromosome 16q24.2 (odds ratio [OR; 95% confidence interval {CI}] = 1.16 [1.10–1.22]; p = 3.2 × 10⁻⁹). The lead single‐nucleotide polymorphism (rs12445022) was also associated with cerebral white matter hyperintensities (OR [95% CI] = 1.10 [1.05–1.16]; p = 5.3 × 10⁻⁵; N = 3,670), but not intracerebral hemorrhage (OR [95% CI] = 0.97 [0.84–1.12]; p = 0.71; 1,545 cases, 1,481 controls). rs12445022 is associated with mRNA expression of ZCCHC14 in arterial tissues (p = 9.4 × 10⁻⁷) and DNA methylation at probe cg16596957 in whole blood (p = 5.3 × 10⁻⁶).

Interpretation

16q24.2 is associated with SVS. Associations of the locus with expression of ZCCHC14 and DNA methylation suggest the locus acts through changes to regulatory elements. Ann Neurol 2017;81:383–394

Trajectories of Callous-Unemotional Traits in Childhood Predict Different Forms of Peer Victimization in Adolescence

Callous-unemotional (CU) traits (e.g., lack of empathy and guilt) differentiate a group of children at particularly high risk for engaging in aggressive behavior, notably bullying. However, little is known about whether youths with CU traits are at risk for being victimized by their peers. We examined the associations between trajectories of CU traits in childhood (between 7 and 12 years old) and peer victimization in adolescence (14 years old). The participants were drawn from the Twins Early Development Study, a longitudinal population-based study of twins born in England and in Wales. The trajectories of CU traits (i.e., stable high, increasing, decreasing and stable low) were identified through general growth mixture modeling. Four forms of peer victimization were considered: physical victimization, verbal victimization, social manipulation, and attacks on property. We found that youths with stable high levels, increasing levels, and decreasing levels of CU traits in childhood had higher levels of physical victimization in adolescence, not explained by other predictors at age 7 (e.g., conduct problems). Youths with increasing levels of CU traits, compared with the ones with stable low levels, also had higher levels of verbal victimization, social manipulation, and attacks on property. Our findings highlight the importance of distinct trajectories of CU traits in accounting for the experience of different forms of peer victimization. Youths with CU traits may benefit from bullying prevention programs, as they are likely to be the targets of peer victimization.

Genetic Factors Influencing Coagulation Factor XIII B-Subunit Contribute to Risk of Ischemic Stroke

Background and Purpose—

Abnormal coagulation has been implicated in the pathogenesis of ischemic stroke, but how this association is mediated and whether it differs between ischemic stroke subtypes is unknown. We determined the shared genetic risk between 14 coagulation factors and ischemic stroke and its subtypes.

Methods—

Using genome-wide association study results for 14 coagulation factors from the population-based TwinsUK sample (N≈2000 for each factor), meta-analysis results from the METASTROKE consortium ischemic stroke genome-wide association study (12 389 cases, 62 004 controls), and genotype data for 9520 individuals from the WTCCC2 ischemic stroke study (3548 cases, 5972 controls—the largest METASTROKE subsample), we explored shared genetic risk for coagulation and stroke. We performed three analyses: (1) a test for excess concordance (or discordance) in single nucleotide polymorphism effect direction across coagulation and stroke, (2) an estimation of the joint effect of multiple coagulation-associated single nucleotide polymorphisms in stroke, and (3) an evaluation of common genetic risk between coagulation and stroke.

Results—

One coagulation factor, factor XIII subunit B (FXIIIB), showed consistent effects in the concordance analysis, the estimation of polygenic risk, and the validation with genotype data, with associations specific to the cardioembolic stroke subtype. Effect directions for FXIIIB-associated single nucleotide polymorphisms were significantly discordant with cardioembolic disease (smallest P=5.7×10⁻⁰⁴); the joint effect of FXIIIB-associated single nucleotide polymorphisms was significantly predictive of ischemic stroke (smallest P=1.8×10⁻⁰⁴) and the cardioembolic subtype (smallest P=1.7×10⁻⁰⁴). We found substantial negative genetic covariation between FXIIIB and ischemic stroke (rG=−0.71, P=0.01) and the cardioembolic subtype (rG=−0.80, P=0.03).

Conclusions—

Genetic markers associated with low FXIIIB levels increase risk of ischemic stroke cardioembolic subtype.

Why do spatial abilities predict mathematical performance?

Spatial ability predicts performance in mathematics and eventual expertise in science, technology and engineering. Spatial skills have also been shown to rely on neuronal networks partially shared with mathematics. Understanding the nature of this association can inform educational practices and intervention for mathematical underperformance. Using data on two aspects of spatial ability and three domains of mathematical ability from 4174 pairs of 12-year-old twins, we examined the relative genetic and environmental contributions to variation in spatial ability and to its relationship with different aspects of mathematics. Environmental effects explained most of the variation in spatial ability (~70%) and in mathematical ability (~60%) at this age, and the effects were the same for boys and girls. Genetic factors explained about 60% of the observed relationship between spatial ability and mathematics, with a substantial portion of the relationship explained by common environmental influences (26% and 14% by shared and non-shared environments respectively). These findings call for further research aimed at identifying specific environmental mediators of the spatial-mathematics relationship.

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

PURPOSE

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

METHOD

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

RESULTS

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

CONCLUSIONS

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Understanding the science-learning environment: A genetically sensitive approach

Previous studies have shown that environmental influences on school science performance increase in importance from primary to secondary school. Here we assess for the first time the relationship between the science-learning environment and science performance using a genetically sensitive approach to investigate the aetiology of this link. 3000 pairs of 14-year-old twins from the UK Twins Early Development Study reported on their experiences of the science-learning environment and were assessed for their performance in science using a web-based test of scientific enquiry. Multivariate twin analyses were used to investigate the genetic and environmental links between environment and outcome. The most surprising result was that the science-learning environment was almost as heritable (43%) as performance on the science test (50%), and showed negligible shared environmental influence (3%). Genetic links explained most (56%) of the association between learning environment and science outcome, indicating gene–environment correlation.

Chaotic homes and children's disruptive behavior: a longitudinal cross-lagged twin study

Chaotic home lives are correlated with behavior problems in children. In the study reported here, we tested whether there was a cross-lagged relation between children’s experience of chaos and their disruptive behaviors (conduct problems and hyperactivity-inattention). Using genetically informative models, we then tested for the first time whether the influence of household chaos on disruptive behavior was environmentally mediated and whether genetic influences on children’s disruptive behaviors accounted for the heritability of household chaos. We measured children’s perceptions of household chaos and parents’ ratings of children’s disruptive behavior at ages 9 and 12 in a sample of 6,286 twin pairs from the Twins Early Development Study (TEDS). There was a phenotypic cross-lagged relation between children’s experiences of household chaos and their disruptive behavior. In genetically informative models, we found that the effect of household chaos on subsequent disruptive behavior was environmentally mediated. However, genetic influences on disruptive behavior did not explain why household chaos was heritable.

Socioeconomic status (SES) and children's intelligence (IQ): in a UK-representative sample SES moderates the environmental, not genetic, effect on IQ

Background

The environment can moderate the effect of genes - a phenomenon called gene-environment (GxE) interaction. Several studies have found that socioeconomic status (SES) modifies the heritability of children's intelligence. Among low-SES families, genetic factors have been reported to explain less of the variance in intelligence; the reverse is found for high-SES families. The evidence however is inconsistent. Other studies have reported an effect in the opposite direction (higher heritability in lower SES), or no moderation of the genetic effect on intelligence.

Methods

Using 8716 twin pairs from the Twins Early Development Study (TEDS), we attempted to replicate the reported moderating effect of SES on children's intelligence at ages 2, 3, 4, 7, 9, 10, 12 and 14: i.e., lower heritability in lower-SES families. We used a twin model that allowed for a main effect of SES on intelligence, as well as a moderating effect of SES on the genetic and environmental components of intelligence.

Results

We found greater variance in intelligence in low-SES families, but minimal evidence of GxE interaction across the eight ages. A power calculation indicated that a sample size of about 5000 twin pairs is required to detect moderation of the genetic component of intelligence as small as 0.25, with about 80% power - a difference of 11% to 53% in heritability, in low- (−2 standard deviations, SD) and high-SES (+2 SD) families. With samples at each age of about this size, the present study found no moderation of the genetic effect on intelligence. However, we found the greater variance in low-SES families is due to moderation of the environmental effect – an environment-environment interaction.

Conclusions

In a UK-representative sample, the genetic effect on intelligence is similar in low- and high-SES families. Children's shared experiences appear to explain the greater variation in intelligence in lower SES.

Chaotic homes and school achievement: a twin study

BACKGROUND

Chaotic homes predict poor school performance. Given that it is known that genes affect both children's experience of household chaos and their school achievement, to what extent is the relationship between high levels of noise and environmental confusion in the home, and children's school performance, mediated by heritable child effects? This is the first study to explore the genetic and environmental pathways between household chaos and academic performance.

METHOD

Children's perceptions of family chaos at ages 9 and 12 and their school performance at age 12 were assessed in more than 2,300 twin pairs. The use of child-specific measures in a multivariate genetic analysis made it possible to investigate the genetic and environmental origins of the covariation between children's experience of chaos in the home and their school achievement.

RESULTS

Children's experience of family chaos and their school achievement were significantly correlated in the expected negative direction (r = -.26). As expected, shared environmental factors explained a large proportion (63%) of the association. However, genetic factors accounted for a significant proportion (37%) of the association between children's experience of household chaos and their school performance.

CONCLUSIONS

The association between chaotic homes and poor performance in school, previously assumed to be entirely environmental in origin, is in fact partly genetic. How children's home environment affects their academic achievement is not simply in the direction environment → child → outcome. Instead, genetic factors that influence children's experience of the disordered home environment also affect how well they do at school. The relationship between the child, their environment and their performance at school is complex: both genetic and environmental factors play a role.

The nature (and nurture) of children's perceptions of family chaos

Chaos in the home is a key environment in cognitive and behavioral development. However, we show that children's experience of home chaos is partly genetically mediated. We assessed children's perceptions of household chaos at ages 9 and 12 in 2337 pairs of twins. Using child-specific reports allowed us to use structural equation modeling to explore the genetic and environmental etiology of children's perceptions of chaos. We found that these perceptions are significantly heritable (22%), with the remainder explained by environmental influences. Finding that genes influence children's experience of chaotic environments has far-reaching implications for how we conceptualize the family home and its impact on cognitive and behavioral development.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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