Michel G. Nivard

The Young Netherlands Twin Register (YNTR): Longitudinal Twin and Family Studies in Over 70,000 Children

The Netherlands Twin Register (NTR) began in 1987 with data collection in twins and their families, including families with newborn twins and triplets. Twenty-five years later, the NTR has collected at least one survey for 70,784 children, born after 1985. For the majority of twins, longitudinal data collection has been done by age-specific surveys. Shortly after giving birth, mothers receive a first survey with items on pregnancy and birth. At age 2, a survey on growth and achievement of milestones is sent. At ages 3, 7, 9/10, and 12 parents and teachers receive a series of surveys that are targeted at the development of emotional and behavior problems. From age 14 years onward, adolescent twins and their siblings report on their behavior problems, health, and lifestyle. When the twins are 18 years and older, parents are also invited to take part in survey studies. In sub-groups of different ages, in-depth phenotyping was done for IQ, electroencephalography , MRI, growth, hormones, neuropsychological assessments, and cardiovascular measures. DNA and biological samples have also been collected and large numbers of twin pairs and parents have been genotyped for zygosity by either micro-satellites or sets of short nucleotide polymorphisms and repeat polymorphisms in candidate genes. Subject recruitment and data collection is still ongoing and the longitudinal database is growing. Data collection by record linkage in the Netherlands is beginning and we expect these combined longitudinal data to provide increased insights into the genetic etiology of development of mental and physical health in children and adolescents.

Genetic and environmental influences interact with age and sex in shaping the human methylome

The methylome is subject to genetic and environmental effects. Their impact may depend on sex and age, resulting in sex- and age-related physiological variation and disease susceptibility. Here we estimate the total heritability of DNA methylation levels in whole blood and estimate the variance explained by common single nucleotide polymorphisms at 411,169 sites in 2,603 individuals from twin families, to establish a catalogue of between-individual variation in DNA methylation. Heritability estimates vary across the genome (mean=19%) and interaction analyses reveal thousands of sites with sex-specific heritability as well as sites where the environmental variance increases with age. Integration with previously published data illustrates the impact of genome and environment across the lifespan at methylation sites associated with metabolic traits, smoking and ageing. These findings demonstrate that our catalogue holds valuable information on locations in the genome where methylation variation between people may reflect disease-relevant environmental exposures or genetic variation.

Genetic and Environmental Stability in Attention Problems Across the Lifespan: Evidence From the Netherlands Twin Register

OBJECTIVE To review findings on attention-deficit/hyperactivity disorder and attention problems (AP) in children, adolescents, and adults, as established in the database of the Netherlands Twin Register and increase the understanding of stability in AP across the lifespan as a function of genetic and environmental influences. METHOD A longitudinal model was fitted on Netherlands Twin Register AP scores from 44,607 child (<12-year-old), adolescent (12- to 18-year-old), and adult (>18-year-old) twins. RESULTS Mean AP showed a downward trend with age. Age-to-age correlations ranged from 0.33 (50-≥60 years old) to 0.73 (10-12 years old). Stability in individual differences in AP was due to genetic and environmental factors, and change was due primarily to environmental factors. Nonadditive genetic influences were present from childhood to adulthood. Total genetic variance decreased slightly throughout aging, whereas environmental variance increased substantially with the switch from maternal to self-ratings at 12 years of age. As a result, heritability coefficients decreased from 0.70 to 0.74 in childhood (maternal ratings) to 0.51 to 0.56 in adolescence (self-ratings), and 0.40 to 0.54 in adulthood (self-ratings). In childhood, male subjects scored higher than female subjects. After the rater switch at 12 years of age, female subjects tended to score higher than male subjects. CONCLUSIONS Stability of AP is the result of genetic and environmental stability. The decrease in estimated heritability at 12 years of age is due to an increase in occasion-specific environmental variance and likely reflects a methodologic effect. Because environmental influences have lasting effects on AP, their early detection is crucial.

Meta-analysis of genome-wide association studies of anxiety disorders

Anxiety disorders (ADs), namely generalized AD, panic disorder and phobias, are common, etiologically complex conditions with a partially genetic basis. Despite differing on diagnostic definitions based on clinical presentation, ADs likely represent various expressions of an underlying common diathesis of abnormal regulation of basic threat–response systems. We conducted genome-wide association analyses in nine samples of European ancestry from seven large, independent studies. To identify genetic variants contributing to genetic susceptibility shared across interview-generated DSM-based ADs, we applied two phenotypic approaches: (1) comparisons between categorical AD cases and supernormal controls, and (2) quantitative phenotypic factor scores (FS) derived from a multivariate analysis combining information across the clinical phenotypes. We used logistic and linear regression, respectively, to analyze the association between these phenotypes and genome-wide single nucleotide polymorphisms. Meta-analysis for each phenotype combined results across the nine samples for over 18 000 unrelated individuals. Each meta-analysis identified a different genome-wide significant region, with the following markers showing the strongest association: for case–control contrasts, rs1709393 located in an uncharacterized non-coding RNA locus on chromosomal band 3q12.3 (P=1.65 × 10−8); for FS, rs1067327 within CAMKMT encoding the calmodulin-lysine N-methyltransferase on chromosomal band 2p21 (P=2.86 × 10−9). Independent replication and further exploration of these findings are needed to more fully understand the role of these variants in risk and expression of ADs.

Power in GWAS: lifting the curse of the clinical cut-off

Although genome-wide association studies (GWAS), in general, facilitated important discovery of new biological knowledge about diseases,1, 2, 3 identified variants for psychiatric disorders explain little variation, and insight into the role of genes in highly heritable psychiatric traits remains poor.4, 5 Low statistical power is seen as the main reason for the failure to locate more variants, and has resulted in a call for larger samples. Our present study, however, shows that better use of (available) phenotypic information can also increase power considerably.

Stability in symptoms of anxiety and depression as a function of genotype and environment: a longitudinal twin study from ages 3 to 63 years.

BACKGROUND The influence of genetic factors on major depressive disorder is lower than on other psychiatric disorders. Heritability estimates mainly derive from cross-sectional studies, and knowledge on the longitudinal aetiology of symptoms of anxiety and depression (SxAnxDep) across the lifespan is limited. We aimed to assess phenotypic, genetic and environmental stability in SxAnxDep between ages 3 and 63 years. METHOD We used a cohort-sequential design combining data from 49 524 twins followed from birth to age ⩾20 years, and from adolescence into adulthood. SxAnxDep were assessed repeatedly with a maximum of eight assessments over a 25-year period. Data were ordered in 30 age groups and analysed with longitudinal genetic models. RESULTS Over age, there was a significant increase during adolescence in mean scores with sex differences (women>men) emerging. Heritability was high in childhood and decreased to 30-40% during adulthood. This decrease in heritability was due to an increase in environmental variance. Phenotypic stability was moderate in children (correlations across ages ~0.5) and high in adolescents (r = 0.6), young adults (r = 0.7), and adults (r = 0.8). Longitudinal stability was mostly attributable to genetic factors. During childhood and adolescence there was also significant genetic innovation, which was absent in adults. Environmental effects contributed to short-term stability. CONCLUSIONS The substantial stability in SxAnxDep is mainly due to genetic effects. The importance of environmental effects increases with age and explains the relatively low heritability of depression in adults. The environmental effects are transient, but the contribution to stability increases with age.

A Genome-Wide Approach to Children's Aggressive Behavior: The EAGLE consortium

Individual differences in aggressive behavior emerge in early childhood and predict persisting behavioral problems and disorders. Studies of antisocial and severe aggression in adulthood indicate substantial underlying biology. However, little attention has been given to genome‐wide approaches of aggressive behavior in children. We analyzed data from nine population‐based studies and assessed aggressive behavior using well‐validated parent‐reported questionnaires. This is the largest sample exploring childrens aggressive behavior to date (N = 18,988), with measures in two developmental stages (N = 15,668 early childhood and N = 16,311 middle childhood/early adolescence). First, we estimated the additive genetic variance of childrens aggressive behavior based on genome‐wide SNP information, using genome‐wide complex trait analysis (GCTA). Second, genetic associations within each study were assessed using a quasi‐Poisson regression approach, capturing the highly right‐skewed distribution of aggressive behavior. Third, we performed meta‐analyses of genome‐wide associations for both the total age‐mixed sample and the two developmental stages. Finally, we performed a gene‐based test using the summary statistics of the total sample. GCTA quantified variance tagged by common SNPs (10–54%). The meta‐analysis of the total sample identified one region in chromosome 2 (2p12) at near genome‐wide significance (top SNP rs11126630, P = 5.30 × 10−8). The separate meta‐analyses of the two developmental stages revealed suggestive evidence of association at the same locus. The gene‐based analysis indicated association of variation within AVPR1A with aggressive behavior. We conclude that common variants at 2p12 show suggestive evidence for association with childhood aggression. Replication of these initial findings is needed, and further studies should clarify its biological meaning.

Genome-wide association study of lifetime cannabis use based on a large meta-analytic sample of 32 330 subjects from the International Cannabis Consortium

Cannabis is the most widely produced and consumed illicit psychoactive substance worldwide. Occasional cannabis use can progress to frequent use, abuse and dependence with all known adverse physical, psychological and social consequences. Individual differences in cannabis initiation are heritable (40–48%). The International Cannabis Consortium was established with the aim to identify genetic risk variants of cannabis use. We conducted a meta-analysis of genome-wide association data of 13 cohorts (N=32 330) and four replication samples (N=5627). In addition, we performed a gene-based test of association, estimated single-nucleotide polymorphism (SNP)-based heritability and explored the genetic correlation between lifetime cannabis use and cigarette use using LD score regression. No individual SNPs reached genome-wide significance. Nonetheless, gene-based tests identified four genes significantly associated with lifetime cannabis use: NCAM1, CADM2, SCOC and KCNT2. Previous studies reported associations of NCAM1 with cigarette smoking and other substance use, and those of CADM2 with body mass index, processing speed and autism disorders, which are phenotypes previously reported to be associated with cannabis use. Furthermore, we showed that, combined across the genome, all common SNPs explained 13–20% (P<0.001) of the liability of lifetime cannabis use. Finally, there was a strong genetic correlation (rg=0.83; P=1.85 × 10−8) between lifetime cannabis use and lifetime cigarette smoking implying that the SNP effect sizes of the two traits are highly correlated. This is the largest meta-analysis of cannabis GWA studies to date, revealing important new insights into the genetic pathways of lifetime cannabis use. Future functional studies should explore the impact of the identified genes on the biological mechanisms of cannabis use.

A genome-wide association meta-analysis of preschool internalizing problems

OBJECTIVE Preschool internalizing problems (INT) are highly heritable and moderately genetically stable from childhood into adulthood. Gene-finding studies are scarce. In this study, the influence of genome-wide measured single nucleotide polymorphisms (SNPs) was investigated in 3 cohorts (total N = 4,596 children) in which INT was assessed with the same instrument, the Child Behavior Checklist (CBCL). METHOD First, genome-wide association (GWA) results were used for density estimation and genome-wide complex trait analysis (GCTA) to calculate the variance explained by all SNPs. Next, a fixed-effect inverse variance meta-analysis of the 3 GWA analyses was carried out. Finally, the overlap in results with prior GWA studies of childhood and adulthood psychiatric disorders and treatment responses was tested by examining whether SNPs associated with these traits jointly showed a significant signal for INT. RESULTS Genome-wide SNPs explained 13% to 43% of the total variance. This indicates that the genetic architecture of INT mirrors the polygenic model underlying adult psychiatric traits. The meta-analysis did not yield a genome-wide significant signal but was suggestive for the PCSK2 gene located on chromosome 20p12.1. SNPs associated with other psychiatric disorders appeared to be enriched for signals with INT (λ = 1.26, p < .03). CONCLUSION Our study provides evidence that INT is influenced by many common genetic variants, each with a very small effect, and that, even as early as age 3, genetic variants influencing INT overlap with variants that play a role in childhood and adulthood psychiatric disorders.

Ultra-rare disruptive and damaging mutations influence educational attainment in the general population

Disruptive, damaging ultra-rare variants in highly constrained genes are enriched in individuals with neurodevelopmental disorders. In the general population, this class of variants was associated with a decrease in years of education (YOE). This effect was stronger among highly brain-expressed genes and explained more YOE variance than pathogenic copy number variation but less than common variants. Disruptive, damaging ultra-rare variants in highly constrained genes influence the determinants of YOE in the general population.