Phenotypic and Genetic Associations Between Preschool Fine Motor Skills and Later Neurodevelopment, Psychopathology, and Educational Achievement.

BACKGROUND: Fine motor skills are heritable and comprise important milestones in development, and some evidence suggests that impairments in fine motor skills are associated with neurodevelopmental conditions, psychiatric disorders, and poor educational outcomes. METHODS: In a preregistered study of 9625 preschool children from TEDS (Twins Early Development Study), fine motor assessments (drawing, block building, folding, and questionnaires) were conducted at 2, 3, and 4 years of age. A cross-age fine motor score was derived using principal component analysis. Multivariate regression analysis was used to examine the relationships between the fine motor score and neurodevelopmental traits, psychopathology, and educational outcomes at 3 later ages (7-8, 12, and 16 years) and cross-age psychopathology composite scores. Polygenic scores (PGSs) were created for attention-deficit/hyperactivity disorder (ADHD), autism, schizophrenia, anxiety, major depressive disorder, obsessive-compulsive disorder, and years of education. We ran single-PGS models and a multi-PGS model. RESULTS: Fine motor skills were negatively associated with neurodevelopmental traits and psychopathology across childhood and adolescence and positively associated with educational achievement in adolescence (ß = 0.25, p < .001). Superior fine motor skills were associated with a higher years-of-education PGS (ß = 0.07, p < .001), a lower ADHD PGS (ß = -0.04, p = .011), and a higher anxiety PGS (ß = 0.03, p = .040). Similarly, the multi-PGS model retained the PGSs for years of education (ß = 0.07), ADHD (ß = -0.03), and anxiety (ß = 0.01). A non-preregistered analysis in an independent preschool sample replicated the ADHD PGS association, but not the years of education or anxiety PGS associations. CONCLUSIONS: Fine motor skills are linked genetically and phenotypically to later neurodevelopment, psychopathology, and educational outcomes. Future work should investigate the mechanisms that underlie the role of fine motor development in later outcomes.

Development of the pupillary light reflex from 9 to 24 months: association with common autism spectrum disorder (ASD) genetic liability and 3-year ASD diagnosis.

BACKGROUND: Although autism spectrum disorder (ASD) is heritable, the mechanisms through which genes contribute to symptom emergence remain unclear. Investigating candidate intermediate phenotypes such as the pupillary light reflex (PLR) prospectively from early in development could bridge genotype and behavioural phenotype. METHODS: Using eye tracking, we longitudinally measured the PLR at 9, 14 and 24 months in a sample of infants (N = 264) enriched for a family history of ASD; 27 infants received an ASD diagnosis at 3 years. We examined the 9- to 24-month developmental trajectories of PLR constriction latency (onset; ms) and amplitude (%) and explored their relation to categorical 3-year ASD outcome, polygenic liability for ASD and dimensional 3-year social affect (SA) and repetitive/restrictive behaviour (RRB) traits. Polygenic scores for ASD (PGSASD ) were calculated for 190 infants. RESULTS: While infants showed a decrease in latency between 9 and 14 months, higher PGSASD was associated with a smaller decrease in latency in the first year (ß = -.16, 95% CI = -0.31, -0.002); infants with later ASD showed a significantly steeper decrease in latency (a putative 'catch-up') between 14 and 24 months relative to those with other outcomes (typical: ß = .54, 95% CI = 0.08, 0.99; other: ß = .53, 95% CI = 0.02, 1.04). Latency development did not associate with later dimensional variation in ASD-related traits. In contrast, change in amplitude was not related to categorical ASD or genetics, but decreasing 9- to 14-month amplitude was associated with higher SA (ß = .08, 95% CI = 0.01, 0.14) and RRB (ß = .05, 95% CI = 0.004, 0.11) traits. CONCLUSIONS: These findings corroborate PLR development as possible intermediate phenotypes being linked to both genetic liability and phenotypic outcomes. Future work should incorporate alternative measures (e.g. functionally informed structural and genetic measures) to test whether distinct neural mechanisms underpin PLR alterations.

Internalising and externalising in early adolescence predict later executive function, not the other way around: a cross-lagged panel analysis.

Developmental changes in the brain networks involved in emotion regulation are thought to contribute to vulnerability to mental health problems during adolescence. Executive control is often viewed as allowing top-down regulation of emotional responses. However, while associations between executive control and mental health are commonly observed in both clinical and non-clinical populations, the direction of these associations remains unclear. Low, or immature, cognitive control could limit emotion regulation. Reversely, high emotionality could impede cognitive functioning. The scarcity of longitudinal studies testing for bi-directional effects, particularly in adolescence, has made it difficult to draw conclusions. This study analysed data from 1,445 participants of a longitudinal cohort in a cross-lagged panel design to understand bi-directional longitudinal associations between executive function and emotional behaviours across adolescence. Executive function was assessed using experimental working memory and inhibitory control tasks, emotional behaviours through parental report of internalising and externalising behaviours. Cross-sectional associations were replicated. Controlling for cross-sectional associations, early executive functions were not found to predict later emotional behaviours. Instead, early emotional behaviours predicted later executive function, with the strongest link observed between early externalising and later working memory. These results suggest that emotional well-being may affect the maturation of executive function during adolescence.

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

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

Look duration at the face as a developmental endophenotype: elucidating pathways to autism and ADHD.

Identifying developmental endophenotypes on the pathway between genetics and behavior is critical to uncovering the mechanisms underlying neurodevelopmental conditions. In this proof-of-principle study, we explored whether early disruptions in visual attention are a unique or shared candidate endophenotype of autism spectrum disorder (ASD) and attention-deficit/hyperactivity disorder (ADHD). We calculated the duration of the longest look (i.e., peak look) to faces in an array-based eye-tracking task for 335 14-month-old infants with and without first-degree relatives with ASD and/or ADHD. We leveraged parent-report and genotype data available for a proportion of these infants to evaluate the relation of looking behavior to familial (n = 285) and genetic liability (using polygenic scores, n = 185) as well as ASD and ADHD-relevant temperament traits at 2 years of age (shyness and inhibitory control, respectively, n = 272) and ASD and ADHD clinical traits at 6 years of age (n = 94).Results showed that longer peak looks at the face were associated with elevated polygenic scores for ADHD (ß = 0.078, p = .023), but not ASD (ß = 0.002, p = .944), and with elevated ADHD traits in mid-childhood (F(1,88) = 6.401, p = .013, $\eta _p^2$=0.068; ASD: F (1,88) = 3.218, p = .076), but not in toddlerhood (ps > 0.2). This pattern of results did not emerge when considering mean peak look duration across face and nonface stimuli. Thus, alterations in attention to faces during spontaneous visual exploration may be more consistent with a developmental endophenotype of ADHD than ASD. Our work shows that dissecting paths to neurodevelopmental conditions requires longitudinal data incorporating polygenic contribution, early neurocognitive function, and clinical phenotypic variation.

Leveraging epigenetics to examine differences in developmental trajectories of social attention: A proof-of-principle study of DNA methylation in infants with older siblings with autism.

Preliminary evidence suggests that changes in DNA methylation, a widely studied epigenetic mechanism, contribute to the etiology of Autism Spectrum Disorder (ASD). However, data is primarily derived from post-mortem brain samples or peripheral tissue from adults. Deep-phenotyped longitudinal infant cohorts are essential to understand how epigenetic modifications relate to early developmental trajectories and emergence of ASD symptoms. We present a proof-of-principle study designed to evaluate the potential of prospective epigenetic studies of infant siblings of children with ASD. Illumina genome-wide 450 K DNA methylation data from buccal swabs was generated for 63 male infants at multiple time-points from 8 months to 2 years of age (total N = 107 samples). 11 of those infants received a diagnosis of ASD at 3 years. We conducted a series of analyses to characterize DNA methylation signatures associated with categorical outcome and neurocognitive measures from parent-report questionnaire, eye-tracking and electro-encephalography. Effects observed across the entire genome (epigenome-wide association analyses) suggest that collecting DNA methylation samples within infant-sibling designs allows for the detection of meaningful signals with smaller sample sizes than previously estimated. Mapping networks of co-methylated probes associated with neural correlates of social attention implicated enrichment of pathways involved in brain development. Longitudinal modelling found covariation between phenotypic traits and DNA methylation levels in the proximity of genes previously associated with cognitive development, although larger samples and more complete datasets are needed to obtain generalizable results. In conclusion, assessment of DNA methylation profiles at multiple time-points in infant-sibling designs is a promising avenue to comprehend developmental origins and mechanisms of ASD.

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

Background: A gap exists in our mechanistic understanding of how genetic and environmental risk factors converge at the molecular level to result in the emergence of autism symptoms. We compared blood-based gene expression signatures in identical twins concordant and discordant for autism spectrum condition (ASC) to differentiate genetic and environmentally driven transcription differences, and establish convergent evidence for biological mechanisms involved in ASC. Methods: Genome-wide gene expression data were generated using RNA-seq on whole blood samples taken from 16 pairs of monozygotic (MZ) twins and seven twin pair members (39 individuals in total), who had been assessed for ASC and autism traits at age 12. Differential expression (DE) analyses were performed between (a) affected and unaffected subjects (N = 36) and (b) within discordant ASC MZ twin pairs (total N = 11) to identify environmental-driven DE. Gene set enrichment and pathway testing was performed on DE gene lists. Finally, an integrative analysis using DNA methylation data aimed to identify genes with consistent evidence for altered regulation in cis. Results: In the discordant twin analysis, three genes showed evidence for DE at FDR < 10%: IGHG4, EVI2A and SNORD15B. In the case-control analysis, four DE genes were identified at FDR < 10% including IGHG4, PRR13P5, DEPDC1B, and ZNF501. We find enrichment for DE of genes curated in the SFARI human gene database. Pathways showing evidence of enrichment included those related to immune cell signalling and immune response, transcriptional control and cell cycle/proliferation. Integrative methylomic and transcriptomic analysis identified a number of genes showing suggestive evidence for cis dysregulation. Limitations: Identical twins stably discordant for ASC are rare, and as such the sample size was limited and constrained to the use of peripheral blood tissue for transcriptomic and methylomic profiling. Given these primary limitations, we focused on transcript-level analysis. Conclusions: Using a cohort of ASC discordant and concordant MZ twins, we add to the growing body of transcriptomic-based evidence for an immune-based component in the molecular aetiology of ASC. Whilst the sample size was limited, the study demonstrates the utility of the discordant MZ twin design combined with multi-omics integration for maximising the potential to identify disease-associated molecular signals.

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

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

Estimation of a significance threshold for epigenome-wide association studies.

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

Tissue-specific patterns of allelically-skewed DNA methylation.

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

Heritability and genome-wide analyses of problematic peer relationships during childhood and adolescence.

Peer behaviour plays an important role in the development of social adjustment, though little is known about its genetic architecture. We conducted a twin study combined with a genome-wide complex trait analysis (GCTA) and a genome-wide screen to characterise genetic influences on problematic peer behaviour during childhood and adolescence. This included a series of longitudinal measures (parent-reported Strengths-and-Difficulties Questionnaire) from a UK population-based birth-cohort (ALSPAC, 4-17 years), and a UK twin sample (TEDS, 4-11 years). Longitudinal twin analysis (TEDS; N ≤ 7,366 twin pairs) showed that peer problems in childhood are heritable (4-11 years, 0.60 < twin-h(2) ≤ 0.71) but genetically heterogeneous from age to age (4-11 years, twin-r(g) = 0.30). GCTA (ALSPAC: N ≤ 5,608, TEDS: N ≤ 2,691) provided furthermore little support for the contribution of measured common genetic variants during childhood (4-12 years, 0.02 < GCTA-h(2)(Meta) ≤ 0.11) though these influences become stronger in adolescence (13-17 years, 0.14 < GCTA-h (2)(ALSPAC) ≤ 0.27). A subsequent cross-sectional genome-wide screen in ALSPAC (N ≤ 6,000) focussed on peer problems with the highest GCTA-heritability (10, 13 and 17 years, 0.0002 < GCTA-P ≤ 0.03). Single variant signals (P ≤ 10(-5)) were followed up in TEDS (N ≤ 2835, 9 and 11 years) and, in search for autism quantitative trait loci, explored within two autism samples (AGRE: N Pedigrees = 793; ACC: N Cases = 1,453/N Controls = 7,070). There was, however, no evidence for association in TEDS and little evidence for an overlap with the autistic continuum. In summary, our findings suggest that problematic peer relationships are heritable but genetically complex and heterogeneous from age to age, with an increase in common measurable genetic variation during adolescence.

Common variation near ROBO2 is associated with expressive vocabulary in infancy.

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

Genetic contributions to visuospatial cognition in Williams syndrome: insights from two contrasting partial deletion patients.

BACKGROUND: Williams syndrome (WS) is a rare neurodevelopmental disorder arising from a hemizygotic deletion of approximately 27 genes on chromosome 7, at locus 7q11.23. WS is characterised by an uneven cognitive profile, with serious deficits in visuospatial tasks in comparison to relatively proficient performance in some other cognitive domains such as language and face processing. Individuals with partial genetic deletions within the WS critical region (WSCR) have provided insights into the contribution of specific genes to this complex phenotype. However, the combinatorial effects of different genes remain elusive. METHODS: WE REPORT ON VISUOSPATIAL COGNITION IN TWO INDIVIDUALS WITH CONTRASTING PARTIAL DELETIONS IN THE WSCR: one female (HR), aged 11 years 9 months, with haploinsufficiency for 24 of the WS genes (up to GTF2IRD1), and one male (JB), aged 14 years 2 months, with the three most telomeric genes within the WSCR deleted, or partially deleted. RESULTS: Our in-depth phenotyping of the visuospatial domain from table-top psychometric, and small- and large-scale experimental tasks reveal a profile in HR in line with typically developing controls, albeit with some atypical features. These data are contrasted with patient JB's atypical profile of strengths and weaknesses across the visuospatial domain, as well as with more substantial visuospatial deficits in individuals with the full WS deletion. CONCLUSIONS: Our findings point to the contribution of specific genes to spatial processing difficulties associated with WS, highlighting the multifaceted nature of spatial cognition and the divergent effects of genetic deletions within the WSCR on different components of visuospatial ability. The importance of general transcription factors at the telomeric end of the WSCR, and their combinatorial effects on the WS visuospatial phenotype are also discussed.

The correlation between reading and mathematics ability at age twelve has a substantial genetic component.

Dissecting how genetic and environmental influences impact on learning is helpful for maximizing numeracy and literacy. Here we show, using twin and genome-wide analysis, that there is a substantial genetic component to children's ability in reading and mathematics, and estimate that around one half of the observed correlation in these traits is due to shared genetic effects (so-called Generalist Genes). Thus, our results highlight the potential role of the learning environment in contributing to differences in a child's cognitive abilities at age twelve.

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

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

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

PURPOSE: Researchers have previously shown that individual differences in measures of receptive language ability at age 12 are highly heritable. In the current study, the authors attempted to identify some of the genes responsible for the heritability of receptive language ability using a genome-wide association approach. METHOD: The authors administered 4 Internet-based measures of receptive language (vocabulary, semantics, syntax, and pragmatics) to a sample of 2,329 twelve-year-olds for whom DNA and genome-wide genotyping were available. Nearly 700,000 single-nucleotide polymorphisms (SNPs) and 1 million imputed SNPs were included in a genome-wide association analysis of receptive language composite scores. RESULTS: No SNP associations met the demanding criterion of genome-wide significance that corrects for multiple testing across the genome ( p < 5 × 10 -8). The strongest SNP association did not replicate in an additional sample of 2,639 twelve-year-olds. CONCLUSIONS: These results indicate that individual differences in receptive language ability in the general population do not reflect common genetic variants that account for more than 3% of the phenotypic variance. The search for genetic variants associated with language skill will require larger samples and additional methods to identify and functionally characterize the full spectrum of risk variants.

Genetics of callous-unemotional behavior in children.

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

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

BACKGROUND: Twin studies have shown that anxiety in a general population sample of children involves both domain-general and trait-specific genetic effects. For this reason, in an attempt to identify genes responsible for these effects, we investigated domain-general and trait-specific genetic associations in the first genome-wide association (GWA) study on anxiety-related behaviours (ARBs) in childhood. METHODS: The sample included 2810 7-year-olds drawn from the Twins Early Development Study (TEDS) with data available for parent-rated anxiety and genome-wide DNA markers. The measure was the Anxiety-Related Behaviours Questionnaire (ARBQ), which assesses four anxiety traits and also yields a general anxiety composite. Affymetrix GeneChip 6.0 DNA arrays were used to genotype nearly 700,000 single-nucleotide polymorphisms (SNPs), and IMPUTE v2 was used to impute more than 1 million SNPs. Several GWA associations from this discovery sample were followed up in another TEDS sample of 4804 children. In addition, Genome-wide Complex Trait Analysis (GCTA) was used on the discovery sample, to estimate the total amount of variance in ARBs that can be accounted for by SNPs on the array. RESULTS: No SNP associations met the demanding criterion of genome-wide significance that corrects for multiple testing across the genome (p<5×10(-8)). Attempts to replicate the top associations did not yield significant results. In contrast to the substantial twin study estimates of heritability which ranged from 0.50 (0.03) to 0.61 (0.01), the GCTA estimates of phenotypic variance accounted for by the SNPs were much lower 0.01 (0.11) to 0.19 (0.12). CONCLUSIONS: Taken together, these GWAS and GCTA results suggest that anxiety--similar to height, weight and intelligence--is affected by many genetic variants of small effect, but unlike these other prototypical polygenic traits, genetic influence on anxiety is not well tagged by common SNPs.