The Young Netherlands Twin Register (YNTR): longitudinal twin and family studies in over 70,000 children

Genome-wide analysis of DNA methylation in buccal cells: a study of monozygotic twins and mQTLs

BACKGROUND

DNA methylation arrays are widely used in epigenome-wide association studies and methylation quantitative trait locus (mQTL) studies. Here, we performed the first genome-wide analysis of monozygotic (MZ) twin correlations and mQTLs on data obtained with the Illumina MethylationEPIC BeadChip (EPIC array) and compared the performance of the EPIC array to the Illumina HumanMethylation450 BeadChip (HM450 array) for buccal-derived DNA.

RESULTS

Good-quality EPIC data were obtained for 102 buccal-derived DNA samples from 49 MZ twin pairs (mean age = 7.5 years, range = 1-10). Differences between MZ twins in the cellular content of buccal swabs were a major driver for differences in their DNA methylation profiles, highlighting the importance to adjust for cellular composition in DNA methylation studies of buccal-derived DNA. After adjusting for cellular composition, the genome-wide mean correlation (r) between MZ twins was 0.21 for the EPIC array, and cis mQTL analysis in 84 twins identified 1,296,323 significant associations (FDR 5%), encompassing 33,749 methylation sites and 616,029 genetic variants. MZ twin correlations were slightly larger (p < 2.2 × 10-16) for novel EPIC probes (N = 383,066, mean r = 0.22) compared to probes that are also present on HM450 (N = 406,822, mean r = 0.20). In line with this observation, a larger percentage of novel EPIC probes was associated with genetic variants (novel EPIC probes with significant mQTL 4.7%, HM450 probes with mQTL 3.9%, p < 2.2 × 10-16). Methylation sites with a large MZ correlation and sites associated with mQTLs were most strongly enriched in epithelial cell DNase I hypersensitive sites (DHSs), enhancers, and histone mark H3K4me3.

CONCLUSIONS

We conclude that the contribution of familial factors to individual differences in DNA methylation and the effect of mQTLs are larger for novel EPIC probes, especially those within regulatory elements connected to active regions specific to the investigated tissue.

Extending Causality Tests with Genetic Instruments: An Integration of Mendelian Randomization with the Classical Twin Design

Although experimental studies are regarded as the method of choice for determining causal influences, these are not always practical or ethical to answer vital questions in health and social research (e.g., one cannot assign individuals to a "childhood trauma condition" in studying the causal effects of childhood trauma on depression). Key to solving such questions are observational studies. Mendelian Randomization (MR) is an influential method to establish causality in observational studies. MR uses genetic variants to test causal relationships between exposures/risk factors and outcomes such as physical or mental health. Yet, individual genetic variants have small effects, and so, when used as instrumental variables, render MR liable to weak instrument bias. Polygenic scores have the advantage of larger effects, but may be characterized by horizontal pleiotropy, which violates a central assumption of MR. We developed the MR-DoC twin model by integrating MR with the Direction of Causation twin model. This model allows us to test pleiotropy directly. We considered the issue of parameter identification, and given identification, we conducted extensive power calculations. MR-DoC allows one to test causal hypotheses and to obtain unbiased estimates of the causal effect given pleiotropic instruments, while controlling for genetic and environmental influences common to the outcome and exposure. Furthermore, the approach allows one to employ strong instrumental variables in the form of polygenic scores, guarding against weak instrument bias, and increasing the power to detect causal effects of exposures on potential outcomes. Beside allowing to test pleiotropy directly, incorporating in MR data collected from relatives provide additional within-family data that resolve additional assumptions like random mating, the absence of the gene-environment interaction/covariance, no dyadic effects. Our approach will enhance and extend MR's range of applications, and increase the value of the large cohorts collected at twin/family registries as they correctly detect causation and estimate effect sizes even in the presence of pleiotropy.

Unraveling the Genetic and Environmental Relationship Between Well-Being and Depressive Symptoms Throughout the Lifespan

Whether well-being and depressive symptoms can be considered as two sides of the same coin is widely debated. The aim of this study was to gain insight into the etiology of the association between well-being and depressive symptoms across the lifespan. In a large twin-design, including data from 43,427 twins between age 7 and 99, we estimated the association between well-being and depressive symptoms throughout the lifespan and assessed genetic and environmental contributions to the observed overlap. For both well-being (range 31-47%) and depressive symptoms (range 49-61%), genetic factors explained a substantial part of the phenotypic variance across the lifespan. Phenotypic correlations between well-being and depressive symptoms across ages ranged from -0.34 in childhood to -0.49 in adulthood. In children, genetic effects explained 49% of the phenotypic correlation while in adolescents and young adults, genetic effects explained 60-77% of the phenotypic correlations. Moderate to high genetic correlations (ranging from -0.59 to -0.66) were observed in adolescence and adulthood, while in childhood environmental correlations were substantial but genetic correlations small. Our results suggest that in childhood genetic and environmental effects are about equally important in explaining the relationship between well-being and depressive symptoms. From adolescence onwards, the role of genetic effects increases compared to environmental effects. These results provided more insights into the etiological underpinnings of well-being and depressive symptoms, possibly allowing to articulate better strategies for health promotion and resource allocation in the future.

Cohort profile: the Geoscience and Health Cohort Consortium (GECCO) in the Netherlands

PURPOSE

In the Netherlands, a great variety of objectively measured geo-data is available, but these data are scattered and measured at varying spatial and temporal scales. The centralisation of these geo-data and the linkage of these data to individual-level data from longitudinal cohort studies enable large-scale epidemiological research on the impact of the environment on public health in the Netherlands. In the Geoscience and Health Cohort Consortium (GECCO), six large-scale and ongoing cohort studies have been enriched with a variety of existing geo-data. Here, we introduce GECCO by describing: (1) the phenotypes of the involved cohort studies, (2) the collected geo-data and their sources, (3) the methodology that was used to link the collected geo-data to individual cohort studies, (4) the similarity of commonly used geo-data between our consortium and the nationwide situation in the Netherlands and (5) the distribution of geo-data within our consortium.

PARTICIPANTS

GECCO includes participants from six prospective cohort studies (eg, 44 657 respondents (18-100 years) in 2006) and it covers all municipalities in the Netherlands. Using postal code information of the participants, geo-data on the address-level, postal code-level as well as neighbourhood-level could be linked to individual-level cohort data.

FINDINGS TO DATE

The geo-data could be successfully linked to almost all respondents of all cohort studies, with successful data-linkage rates ranging from 97.1% to 100.0% between cohort studies. The results show variability in geo-data within and across cohorts. GECCO increases power of analyses, provides opportunities for cross-checking and replication, ensures sufficient geographical variation in environmental determinants and allows for nuanced analyses on specific subgroups.

FUTURE PLANS

GECCO offers unique opportunities for (longitudinal) studies on the complex relationships between the environment and health outcomes. For example, GECCO will be used for further research on environmental determinants of physical/psychosocial functioning and lifestyle behaviours.

Dopaminergic Genetic Variants and Voluntary Externally Paced Exercise Behavior

PURPOSE

Most candidate gene studies on the neurobiology of voluntary exercise behavior have focused on the dopaminergic signaling pathway and its role in the mesolimbic reward system. We hypothesized that dopaminergic candidate genes may influence exercise behavior through additional effects on executive functioning and that these effects are only detected when the types of exercise activity are taken into account.

METHODS

Data on voluntary exercise behavior and at least one single-nucleotide polymorphism/variable number of tandem repeat (VNTR) were available for 12,929 participants of the Netherlands Twin Registry. Exercise activity was classified as externally paced if a high level of executive function skill was required. The total volume of voluntary exercise (minutes per week) as well as the volume specifically spent on externally paced activities were tested for association with nine functional dopaminergic polymorphisms (DRD1: rs265981, DRD2/ANKK1: rs1800497, DRD3: rs6280, DRD4: VNTR 48 bp, DRD5: VNTR 130-166 bp, DBH: rs2519152, DAT1: VNTR 40 bp, COMT: rs4680, MAOA: VNTR 30 bp), a polygenic score (PGS) based on nine alleles leading to lower dopamine responsiveness, and a PGS based on three alleles associated with both higher reward sensitivity and better executive functioning (DRD2/ANKK1: "G" allele, COMT: Met allele, DAT1: 440-bp allele).

RESULTS

No association with total exercise volume or externally paced exercise volume was found for individual alleles or the nine-allele PGS. The volume of externally paced exercise behavior was significantly associated with the reward and executive function congruent PGS. This association was driven by the DAT1 440-bp and COMT Met allele, which acted as increaser alleles for externally paced exercise behavior.

CONCLUSIONS

Taking into account the types of exercise activity may increase the success of identifying genetic variants and unraveling the neurobiology of voluntary exercise behavior.

Genetic and Environmental Influences on Self-Control: Assessing Self-Control with the ASEBA Self-Control Scale

This study used a theoretically-derived set of items of the Achenbach System of Empirically Based Assessment to develop the Achenbach Self-Control Scale (ASCS) for 7-16 year olds. Using a large dataset of over 20,000 children, who are enrolled in the Netherlands Twin Register, we demonstrated the psychometric properties of the ASCS for parent-, self- and teacher-report by examining internal and criterion validity, and inter-rater and test-retest reliability. We found associations between the ASCS and measures of well-being, educational achievement, and substance use. Next, we applied the classical twin design to estimate the genetic and environmental contributions to self-control. Genetic influences accounted for 64-75% of the variance in self-control based on parent- and teacher-report (age 7-12), and for 47-49% of the variance in self-control based on self-report (age 12-16), with the remaining variance accounted by non-shared environmental influences. In conclusion, we developed a validated and accessible self-control scale, and show that genetic influences explain a majority of the individual differences in self-control across youth aged 7-16 years.

Association between structural brain network efficiency and intelligence increases during adolescence

Adolescence represents an important period during which considerable changes in the brain take place, including increases in integrity of white matter bundles, and increasing efficiency of the structural brain network. A more efficient structural brain network has been associated with higher intelligence. Whether development of structural network efficiency is related to intelligence, and if so to which extent genetic and environmental influences are implicated in their association, is not known. In a longitudinal study, we mapped FA-weighted efficiency of the structural brain network in 310 twins and their older siblings at an average age of 10, 13, and 18 years. Age-trajectories of global and local FA-weighted efficiency were related to intelligence. Contributions of genes and environment were estimated using structural equation modeling. Efficiency of brain networks changed in a non-linear fashion from childhood to early adulthood, increasing between 10 and 13 years, and leveling off between 13 and 18 years. Adolescents with higher intelligence had higher global and local network efficiency. The dependency of FA-weighted global efficiency on IQ increased during adolescence (rph =0.007 at age 10; 0.23 at age 18). Global efficiency was significantly heritable during adolescence (47% at age 18). The genetic correlation between intelligence and global and local efficiency increased with age; genes explained up to 87% of the observed correlation at age 18. In conclusion, the brain's structural network differentiates depending on IQ during adolescence, and is under increasing influence of genes that are also associated with intelligence as it develops from late childhood to adulthood.

Genetic and environmental influences on conduct and antisocial personality problems in childhood, adolescence, and adulthood

Conduct problems in children and adolescents can predict antisocial personality disorder and related problems, such as crime and conviction. We sought an explanation for such predictions by performing a genetic longitudinal analysis. We estimated the effects of genetic, shared environmental, and unique environmental factors on variation in conduct problems measured at childhood and adolescence and antisocial personality problems measured at adulthood and on the covariation across ages. We also tested whether these estimates differed by sex. Longitudinal data were collected in the Netherlands Twin Register over a period of 27 years. Age appropriate and comparable measures of conduct and antisocial personality problems, assessed with the Achenbach System of Empirically Based Assessment, were available for 9783 9-10-year-old, 6839 13-18-year-old, and 7909 19-65-year-old twin pairs, respectively; 5114 twins have two or more assessments. At all ages, men scored higher than women. There were no sex differences in the estimates of the genetic and environmental influences. During childhood, genetic and environmental factors shared by children in families explained 43 and 44% of the variance of conduct problems, with the remaining variance due to unique environment. During adolescence and adulthood, genetic and unique environmental factors equally explained the variation. Longitudinal correlations across age varied between 0.20 and 0.38 and were mainly due to stable genetic factors. We conclude that shared environment is mainly of importance during childhood, while genetic factors contribute to variation in conduct and antisocial personality problems at all ages, and also underlie its stability over age.

Bullying and Victimization: The Effect of Close Companionship

Peer bullying and victimization are a widespread phenomenon among school-age children and can have detrimental effects on the development of children. To examine whether having a close companion during childhood increases or decreases risk of victimization and bullying, this study compared twins to singleton children. A large group of twins (n = 9,909) were included who were compared to their related non-twin siblings (n = 1,534) aged 7-12 from the Netherlands Twin Register, thus creating optimal matching between twins and non-twins. Bullying and victimization were each based on a four-item scale filled out by their teachers. Prevalence rates for either bullying or victimization did not differ between twins and singletons. In total, in the past couple of months, 36% of children bullied peers moderately to severely, and 35% suffered moderately to severely from victimization. Boys were more likely to bully and were more prone to becoming a victim than girls. The most notable finding is that female twin pairs placed together in the same classroom did not bully more often, but were victimized less often, thus pointing to a protective effect of having a close companion in the classroom.

The Genetic Overlap Between Hair and Eye Color

We identified the genetic variants for eye color by Genome-Wide Association Study (GWAS) in a Dutch Caucasian family-based population sample and examined the genetic correlation between hair and eye color using data from unrelated participants from the Netherlands Twin Register. With the Genome-wide Complex Trait Analysis software package, we found strong genetic correlations between various combinations of hair and eye colors. The strongest positive correlations were found for blue eyes with blond hair (0.87) and brown eyes with dark hair (0.71), whereas blue eyes with dark hair and brown eyes with blond hair showed the strongest negative correlations (-0.64 and -0.94, respectively). Red hair with green/hazel eyes showed the weakest correlation (-0.14). All analyses were corrected for age and sex, and we explored the effects of correcting for principal components (PCs) that represent ancestry and describe the genetic stratification of the Netherlands. When including the first three PCs as covariates, the genetic correlations between the phenotypes disappeared. This is not unexpected since hair and eye colors strongly indicate the ancestry of an individual. This makes it difficult to separate the effects of population stratification and the true genetic effects of variants on these particular phenotypes.

Genetic Overlap Between Schizophrenia and Developmental Psychopathology: Longitudinal and Multivariate Polygenic Risk Prediction of Common Psychiatric Traits During Development

Background

Several nonpsychotic psychiatric disorders in childhood and adolescence can precede the onset of schizophrenia, but the etiology of this relationship remains unclear. We investigated to what extent the association between schizophrenia and psychiatric disorders in childhood is explained by correlated genetic risk factors.

Methods

Polygenic risk scores (PRS), reflecting an individual's genetic risk for schizophrenia, were constructed for 2588 children from the Netherlands Twin Register (NTR) and 6127 from the Avon Longitudinal Study of Parents And Children (ALSPAC). The associations between schizophrenia PRS and measures of anxiety, depression, attention deficit hyperactivity disorder (ADHD), and oppositional defiant disorder/conduct disorder (ODD/CD) were estimated at age 7, 10, 12/13, and 15 years in the 2 cohorts. Results were then meta-analyzed, and a meta-regression analysis was performed to test differences in effects sizes over, age and disorders.

Results

Schizophrenia PRS were associated with childhood and adolescent psychopathology. Meta-regression analysis showed differences in the associations over disorders, with the strongest association with childhood and adolescent depression and a weaker association for ODD/CD at age 7. The associations increased with age and this increase was steepest for ADHD and ODD/CD. Genetic correlations varied between 0.10 and 0.25.

Conclusion

By optimally using longitudinal data across diagnoses in a multivariate meta-analysis this study sheds light on the development of childhood disorders into severe adult psychiatric disorders. The results are consistent with a common genetic etiology of schizophrenia and developmental psychopathology as well as with a stronger shared genetic etiology between schizophrenia and adolescent onset psychopathology.

The effects of parental education on exercise behavior in childhood and youth: a study in Dutch and Finnish twins

Twin studies have estimated the relative contribution of genes and the environment to variance in exercise behavior and it is known that parental education positively affects exercise levels. This study investigates the role of parental education as a potential modifier of variance in exercise behavior from age 7 to 18 years. The study is based on large datasets from the Netherlands Twin Register (NTR: N = 24 874 twins; surveys around the ages of 7, 10, 12, 14, 16 and 18 years) and two Finnish twin cohorts (FinnTwin12: N = 4399; 12, 14 and 17 years; FinnTwin16: N = 4648; 16, 17 and 18 years). Regular participation in moderate-to-vigorous exercise activities during leisure time was assessed by survey. Parental education was dichotomized ("both parents with a low education" vs "at least one parent with a high education"). The mean in exercise behavior tended to be higher and the variance tended to be lower in children of high educated parents. Evidence for gene-by-environment interaction was weak. To develop successful interventions that specifically target children of low educated parents, the mechanisms causing the mean and variance differences between the two groups should be better understood.

Sum Scores in Twin Growth Curve Models: Practicality Versus Bias

To study behavioral or psychiatric phenotypes, multiple indices of the behavior or disorder are often collected that are thought to best reflect the phenotype. Combining these items into a single score (e.g. a sum score) is a simple and practical approach for modeling such data, but this simplicity can come at a cost in longitudinal studies, where the relevance of individual items often changes as a function of age. Such changes violate the assumptions of longitudinal measurement invariance (MI), and this violation has the potential to obfuscate the interpretation of the results of latent growth models fit to sum scores. The objectives of this study are (1) to investigate the extent to which violations of longitudinal MI lead to bias in parameter estimates of the average growth curve trajectory, and (2) whether absence of MI affects estimates of the heritability of these growth curve parameters. To this end, we analytically derive the bias in the estimated means and variances of the latent growth factors fit to sum scores when the assumption of longitudinal MI is violated. This bias is further quantified via Monte Carlo simulation, and is illustrated in an empirical analysis of aggression in children aged 3-12 years. These analyses show that measurement non-invariance across age can indeed bias growth curve mean and variance estimates, and our quantification of this bias permits researchers to weigh the costs of using a simple sum score in longitudinal studies. Simulation results indicate that the genetic variance decomposition of growth factors is, however, not biased due to measurement non-invariance across age, provided the phenotype is measurement invariant across birth-order and zygosity in twins.

Genetic transmission of reading ability

Reading is the processing of written language. Family resemblance for reading (dis)ability might be due to transmission of a genetic liability or due to family environment, including cultural transmission from parents to offspring. Familial-risk studies exploring neurobehavioral precursors for dyslexia and twin studies can only speak to some of these issues, but a combined twin-family study can resolve the nature of the transmitted risk. Word-reading fluency scores of 1100 participants from 431 families (with twins, siblings and their parents) were analyzed to estimate genetic and environmental sources of variance, and to test the presence of assortative mating and cultural transmission. Results show that variation in reading ability is mainly caused by additive and non-additive genetic factors (64%). The substantial assortative mating (r_{father-mother}=0.38) has scientific and clinical implications. We conclude that parents and offspring tend to resemble each other for genetic reasons, and not due to cultural transmission.

The etiology of autistic traits in preschoolers: a population-based twin study

BACKGROUND

Autism Spectrum Disorders (ASD) are highly heritable, but the exact etiological mechanisms underlying the condition are still unclear.

METHODS

Using a multiple rater twin design in a large sample of general population preschool twins, this study aimed to (a) estimate the contribution of genetic and environmental factors to autistic traits, controlling for the possible effects of rater bias, (b) to explore possible sex differences in etiology and (c) to investigate the discordance in autistic traits in monozygotic and same-sex dizygotic twin pairs. The Netherlands Twin Register collected maternal and paternal ratings on autistic traits from a general population of 38,798 three-year-old twins. Autistic traits were assessed with the DSM-oriented Pervasive Developmental Problems scale of the Child Behavior Check List for preschoolers (1½-5 years).

RESULTS

Mother and fathers showed high agreement in their assessment of autistic traits (r = .60-.66). Differences between children in autistic traits were largely accounted for by genetic effects (boys: 78% and girls: 83%). Environmental effects that are unique to a child also played a modest role. Environmental effects shared by children growing up in the same family were negligible, once rater bias was controlled for. While the prevalence for clinical ASD is higher in boys than in girls, this study did not find evidence for striking differences in the etiology of autistic traits across the sexes. Even though the heritability was high, 29% of MZ twin pairs were discordant for high autistic traits (clinical range vs. normal development), suggesting that despite high genetic risk, environmental factors might lead to resilience, unaffected status in the context of genetic risk, in some children.

CONCLUSIONS

It is important to focus future research on risk factors that might interplay with a genetic disposition for ASD, but also on protective factors that make a difference in the lives of children at genetic risk.

Disentangling Heterogeneity of Childhood Disruptive Behavior Problems Into Dimensions and Subgroups

OBJECTIVE

Irritable and oppositional behaviors are increasingly considered as distinct dimensions of oppositional defiant disorder. However, few studies have explored this multidimensionality across the broader spectrum of disruptive behavior problems (DBPs). This study examined the presence of dimensions and distinct subgroups of childhood DBPs, and the cross-sectional and longitudinal associations between these dimensions.

METHOD

Using factor mixture models (FMMs), the presence of dimensions and subgroups of DBPs was assessed in the Generation R Study at ages 6 (n = 6,209) and 10 (n = 4,724) years. Replications were performed in two population-based cohorts (Netherlands Twin Registry, n = 4,402, and Swedish Twin Study of Child and Adolescent Development, n = 1,089) and a clinical sample (n = 1,933). We used cross-lagged modeling in the Generation R Study to assess cross-sectional and longitudinal associations between dimensions. DBPs were assessed using mother-reported responses to the Child Behavior Checklist.

RESULTS

Empirically obtained dimensions of DBPs were oppositional behavior (age 6 years), disobedient behavior, rule-breaking behavior (age 10 years), physical aggression, and irritability (both ages). FMMs suggested that one-class solutions had the best model fit for all dimensions in all three population-based cohorts. Similar results were obtained in the clinical sample. All three dimensions, including irritability, predicted subsequent physical aggression (range, 0.08-0.16).

CONCLUSION

This study showed that childhood DBPs should be regarded as a multidimensional phenotype rather than comprising distinct subgroups. Incorporating multidimensionality will improve diagnostic accuracy and refine treatment. Future studies need to address the biological validity of the DBP dimensions observed in this study; herein lies an important opportunity for neuroimaging and genetic measures.

Associations between subjective well-being and subcortical brain volumes

To study the underpinnings of individual differences in subjective well-being (SWB), we tested for associations of SWB with subcortical brain volumes in a dataset of 724 twins and siblings. For significant SWB-brain associations we probed for causal pathways using Mendelian Randomization (MR) and estimated genetic and environmental contributions from twin modeling. Another independent measure of genetic correlation was obtained from linkage disequilibrium (LD) score regression on published genome-wide association summary statistics. Our results indicated associations of SWB with hippocampal volumes but not with volumes of the basal ganglia, thalamus, amygdala, or nucleus accumbens. The SWB-hippocampus relations were nonlinear and characterized by lower SWB in subjects with relatively smaller hippocampal volumes compared to subjects with medium and higher hippocampal volumes. MR provided no evidence for an SWB to hippocampal volume or hippocampal volume to SWB pathway. This was in line with twin modeling and LD-score regression results which indicated non-significant genetic correlations. We conclude that low SWB is associated with smaller hippocampal volume, but that genes are not very important in this relationship. Instead other etiological factors, such as exposure to stress and stress hormones, may exert detrimental effects on SWB and the hippocampus to bring about the observed association.

Heritability of the affective response to exercise and its correlation to exercise behavior

OBJECTIVES

Individual differences in adolescent exercise behavior are strongly influenced by genetic factors. The affective response to exercise is a potential source of these genetic influences. To test its role in the motivation to exercise, we estimated the heritability of the affective responses during and after exercise and the overlap with the genetic factors influencing regular voluntary exercise behavior.

DESIGN

226 twin pairs and 38 siblings completed two submaximal exercise tests on a cycle ergometer and a treadmill and a maximal exercise test on a cycle ergometer. Affective responses were assessed by the Feeling Scale (FS), Borg's Rating of Perceived Exertion (RPE) and the Activation-Deactivation Adjective Checklist (AD ACL).

METHODS

Multivariate structural equation modeling was used to estimate heritability of the affective responses during and after submaximal and maximal exercise and the (genetic) correlation with self-reported regular voluntary exercise behavior over the past year.

RESULTS

Genetic factors explained 15% of the individual differences in FS responses during the cycle ergometer test, as well as 29% and 35% of the individual differences in RPE during the cycle ergometer and treadmill tests, respectively. For the AD ACL scales, heritability estimates ranged from 17% to 37% after submaximal exercise and from 12% to 37% after maximal exercise. Without exception, more positive affective responses were associated with higher amounts of regular exercise activity (.15 < r < .21) and this association was accounted for by an overlap in genetic factors influencing affective responding and exercise behavior.

CONCLUSIONS

We demonstrate low to moderate heritability estimates for the affective response during and after exercise and significant (genetic) associations with regular voluntary exercise behavior. These innate individual differences in the affective responses to exercise should be taken into account in interventions aiming to motivate adolescents to adopt and maintain regular exercise.

Attention Deficit Hyperactivity Disorder Symptoms and Low Educational Achievement: Evidence Supporting A Causal Hypothesis

Attention Deficit Hyperactivity Disorder (ADHD) and educational achievement are negatively associated in children. Here we test the hypothesis that there is a direct causal effect of ADHD on educational achievement. The causal effect is tested in a genetically sensitive design to exclude the possibility of confounding by a third factor (e.g. genetic pleiotropy) and by comparing educational achievement and secondary school career in children with ADHD who take or do not take methylphenidate. Data on ADHD symptoms, educational achievement and methylphenidate usage were available in a primary school sample of ~10,000 12-year-old twins from the Netherlands Twin Register. A substantial group also had longitudinal data at ages 7-12 years. ADHD symptoms were cross-sectionally and longitudinally, associated with lower educational achievement at age 12. More ADHD symptoms predicted a lower-level future secondary school career at age 14-16. In both the cross-sectional and longitudinal analyses, testing the direct causal effect of ADHD on educational achievement, while controlling for genetic and environmental factors, revealed an association between ADHD symptoms and educational achievement independent of genetic and environmental pleiotropy. These findings were confirmed in MZ twin intra-pair differences models, twins with more ADHD symptoms scored lower on educational achievement than their co-twins. Furthermore, children with ADHD medication, scored significantly higher on the educational achievement test than children with ADHD who did not use medication. Taken together, the results are consistent with a direct causal effect of ADHD on educational achievement.

Psychopathology in 7-year-old children: Differences in maternal and paternal ratings and the genetic epidemiology

The assessment of children's psychopathology is often based on parental report. Earlier studies have suggested that rater bias can affect the estimates of genetic, shared environmental and unique environmental influences on differences between children. The availability of a large dataset of maternal as well as paternal ratings of psychopathology in 7-year old children enabled (i) the analysis of informant effects on these assessments, and (ii) to obtain more reliable estimates of the genetic and non-genetic effects. DSM-oriented measures of affective, anxiety, somatic, attention-deficit/hyperactivity, oppositional-defiant, conduct, and obsessive-compulsive problems were rated for 12,310 twin pairs from the Netherlands Twin Register by mothers (N = 12,085) and fathers (N = 8,516). The effects of genetic and non-genetic effects were estimated on the common and rater-specific variance. For all scales, mean scores on maternal ratings exceeded paternal ratings. Parents largely agreed on the ranking of their child's problems (r 0.60-0.75). The heritability was estimated over 55% for maternal and paternal ratings for all scales, except for conduct problems (44-46%). Unbiased shared environmental influences, i.e., on the common variance, were significant for affective (13%), oppositional (13%), and conduct problems (37%). In clinical settings, different cutoffs for (sub)clinical scores could be applied to paternal and maternal ratings of their child's psychopathology. Only for conduct problems, shared environmental and genetic influences explain an equal amount in differences between children. For the other scales, genetic factors explain the majority of the variance, especially for the common part that is free of rater bias. © 2016 The Authors. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics Published by Wiley Periodicals, Inc.

Prevalence of dieting and fear of weight gain across ages: a community sample from adolescents to the elderly

OBJECTIVES

The current study aimed to define the prevalence of dieting and fear of weight gain among men and women across the entire lifespan and identify factors associated with them.

METHODS

Data were available for 31,636 participants (60.2% women; age 13-98 years) from the Netherlands Twin Register. Dieting and fear of weight gain were described by age and sex. Associations with BMI, exercise behavior, urbanization and educational attainment were examined by regression analyses in 19,294 participants.

RESULTS

Dieting was most frequently reported by 35- to 65-year-old women (56.6-63%), and 45- to 65-year-old men (31.7-31.9%). Fear of weight gain was most prevalent in women between 16 and 25 (73.2-74.3%), and in 25- to 55-year-old men (43.2-46.1%). In addition to sex and BMI, dieting and fear of weight gain were associated with each other. Furthermore, fear was associated with the age × sex interaction and educational attainment.

CONCLUSIONS

Dieting and fear of weight gain is common during the entire lifespan for women, but is also endorsed by a substantial number of men. Given the low rate of overweight in young women, the high levels of fear of weight gain are striking.

Heritability of Behavioral Problems in 7-Year Olds Based on Shared and Unique Aspects of Parental Views

In studies of child psychopathology, phenotypes of interest are often obtained by parental ratings. When behavioral ratings are obtained in the context of a twin study, this allows for the decomposition of the phenotypic variance, into a genetic and a non-genetic part. If a phenotype is assessed by a single rater, heritability is based on the child’s behavior as expressed in the presence of that particular rater, whereas heritability based on assessments by multiple raters allows for the estimation of the heritability of the phenotype based on rater agreement, as well as the heritability of the rater specific view of the behavior. The aim of this twin study was to quantify the rater common and rater specific contributions to the variation in children’s behavioral problems. We estimated the heritability of maternal and paternal ratings of the Child Behavior Checklist (CBCL) 6–18 empirical emotional and behavioral problem scales in a large sample of 12,310 7-year old Dutch twin pairs. Between 30 and 59% of variation in the part of the phenotype parents agree upon was explained by genetic effects. Common environmental effects that make children in the same family similar explained less variance, ranging between 0 and 32%. For unique views of their children’s behavioral problems, heritability ranged between 0 and 20% for maternal and between 0 and 22% for paternal views. Between 7 and 24% of the variance was accounted for by common environmental factors specific to mother and father’s views. The proportion of rater shared and rater specific heritability can be translated into genetic correlations between parental views and inform the design and interpretation of results of molecular genetic studies. Genetic correlations were nearly or above 0.7 for all CBCL based psychopathology scales. Such large genetic correlations suggest two practical guidelines for genome-wide association studies (GWAS): when studies have collected data from either fathers or mothers, the shared genetic aetiology in parental ratings indicates that is possible to analyze paternal and maternal assessments in a single GWAS or meta-analysis. Secondly, if a study has collected information from both parents, a gain in statistical power may be realized in GWAS by the simultaneous analysis of the data.

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.

A powerful phenotype for gene-finding studies derived from trajectory analyses of symptoms of anxiety and depression between age seven and 18

To investigate the utility of longitudinal data in genetic analyses of symptoms of anxiety and depression, we assessed individual differences between age 7 and 18 using growth mixture models, and investigated the genetic and non-genetic factors contributing to the trajectories. Mothers of 7,706 girl and 7,418 boy twins from the Netherlands Twin Register rated the anxious depression scale (SxAnxDep) of the Child Behavior Check List (CBCL) at age 7, 10, and 12 years. Two thousand seven hundred and six girl and 1,856 boy twins completed the Youth Self Report (YSR) at age 14, 16, and 18. While individual trajectories varied considerably, these differences were largely idiosyncratic and could not be grouped into separate latent classes with class-specific average growth curves. The intercept, which reflects the individuals' baseline level of SxAnxDep across time, explained 55-58% of the total phenotypic variance. The slope factor, which captures a common average trend over time, did not explain variance in the phenotype. This finding also underlines the high level of idiosyncrasy of trajectories that lack a common longitudinal structure. The analyses of twin data showed that the random intercept factor of SxAnxDep during childhood and during adolescence is considerably more heritable than the observations at any single age, namely between 60% and 84%. One explanation is that different factors contribute to the level of symptoms of anxiety and depression at any given time point, including temporary events and emotions. When considering baseline stability, these temporary influences average out, with the result of a more reliable and more heritable phenotype. © 2015 Wiley Periodicals, Inc.

Heritability of heart rate recovery and vagal rebound after exercise

Purpose

The prognostic power of heart rate recovery (HRR) after exercise has been well established but the exact origin of individual differences in HRR remains unclear. This study aims to estimate the heritability of HRR and vagal rebound after maximal exercise in adolescents. Furthermore, the role of voluntary regular exercise behavior (EB) in HRR and vagal rebound is tested.

Methods

491 healthy adolescent twins and their siblings were recruited for maximal exercise testing, followed by a standardized cooldown with measurement of the electrocardiogram and respiratory frequency. Immediate and long-term HRR (HRR60 and HRR180) and vagal rebound (heart rate variability in the respiratory frequency range) were assessed 1 and 3 min after exercise. Multivariate twin modeling was used to estimate heritability of all measured variables and to compute the genetic contribution to their covariance.

Results

Heritability of HRR60, HRR180 and immediate and long-term vagal rebound is 60 % (95 % CI: 48–67), 65 % (95 % CI: 54–73), 23 % (95 % CI: 11–35) and 3 % (95 % CI: 0–11), respectively. We find evidence for two separate genetic factors with one factor influencing overall cardiac vagal control, including resting heart rate and respiratory sinus arrhythmia, and a specific factor for cardiac vagal exercise recovery. EB was only modestly associated with resting heart rate (r = −0.27) and HRR (rHRR60 = 0.10; rHRR180 = 0.19) with very high genetic contribution to these associations (88–91 %).

Conclusions

Individual differences in HRR and immediate vagal rebound can to a large extent be explained by genetic factors. These innate cardiac vagal exercise recovery factors partly reflect the effects of heritable differences in EB.

The genetic architecture of body mass index from infancy to adulthood modified by parental education

OBJECTIVE

A higher prevalence of obesity in lower socioeconomic classes is common in Western societies. This study examined the role of gene-environment interactions in the association between parental education and body mass index (BMI) from infancy to the onset of adulthood.

METHODS

Parentally reported BMI from 1 to 13 and self-reported BMI from 14 to 20 years of age were collected in 16,646 complete Dutch twin pairs and analyzed by genetic twin modeling.

RESULTS

At 7 to 8 years of age, children whose parents had middle or low educational levels had more excess weight than the children of more highly educated parents, and the difference increased until 18 to 20 years of age. The major part of the BMI variation was explained by additive genetic factors (a(2)  = 0.55-0.85), but environmental factors common for co-twins also played a significant role, especially from 3 to 7-8 years of age (c(2)  = 0.15-0.29). The genetic variation in BMI was higher in children whose parents had middle or low educational levels compared with children whose parents had a high educational level.

CONCLUSIONS

The interaction between genetic factors and the childhood social environment may contribute to the formation of socioeconomic differences in obesity.

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.

Genetic variants associated with subjective well-being, depressive symptoms, and neuroticism identified through genome-wide analyses

Very few genetic variants have been associated with depression and neuroticism, likely because of limitations on sample size in previous studies. Subjective well-being, a phenotype that is genetically correlated with both of these traits, has not yet been studied with genome-wide data. We conducted genome-wide association studies of three phenotypes: subjective well-being (n = 298,420), depressive symptoms (n = 161,460), and neuroticism (n = 170,911). We identify 3 variants associated with subjective well-being, 2 variants associated with depressive symptoms, and 11 variants associated with neuroticism, including 2 inversion polymorphisms. The two loci associated with depressive symptoms replicate in an independent depression sample. Joint analyses that exploit the high genetic correlations between the phenotypes (|ρ^| ≈ 0.8) strengthen the overall credibility of the findings and allow us to identify additional variants. Across our phenotypes, loci regulating expression in central nervous system and adrenal or pancreas tissues are strongly enriched for association.

Identification of Common Genetic Variants Influencing Spontaneous Dizygotic Twinning and Female Fertility

Spontaneous dizygotic (DZ) twinning occurs in 1%-4% of women, with familial clustering and unknown physiological pathways and genetic origin. DZ twinning might index increased fertility and has distinct health implications for mother and child. We performed a GWAS in 1,980 mothers of spontaneous DZ twins and 12,953 control subjects. Findings were replicated in a large Icelandic cohort and tested for association across a broad range of fertility traits in women. Two SNPs were identified (rs11031006 near FSHB, p = 1.54 × 10(-9), and rs17293443 in SMAD3, p = 1.57 × 10(-8)) and replicated (p = 3 × 10(-3) and p = 1.44 × 10(-4), respectively). Based on ∼90,000 births in Iceland, the risk of a mother delivering twins increased by 18% for each copy of allele rs11031006-G and 9% for rs17293443-C. A higher polygenic risk score (PRS) for DZ twinning, calculated based on the results of the DZ twinning GWAS, was significantly associated with DZ twinning in Iceland (p = 0.001). A higher PRS was also associated with having children (p = 0.01), greater lifetime parity (p = 0.03), and earlier age at first child (p = 0.02). Allele rs11031006-G was associated with higher serum FSH levels, earlier age at menarche, earlier age at first child, higher lifetime parity, lower PCOS risk, and earlier age at menopause. Conversely, rs17293443-C was associated with later age at last child. We identified robust genetic risk variants for DZ twinning: one near FSHB and a second within SMAD3, the product of which plays an important role in gonadal responsiveness to FSH. These loci contribute to crucial aspects of reproductive capacity and health.

Testing Causal Effects of Maternal Smoking During Pregnancy on Offspring's Externalizing and Internalizing Behavior

Maternal smoking during pregnancy (SDP) is associated with increased risk of externalizing and internalizing behaviors in offspring. Two explanations (not mutually exclusive) for this association are direct causal effects of maternal SDP and the effects of genetic and environmental factors common to parents and offspring which increase smoking as well as problem behaviors. Here, we examined the associations between parental SDP and mother rated offspring externalizing and internalizing behaviors (rated by the Child Behavior Checklist/2-3) at age three in a population-based sample of Dutch twins (N = 15,228 pairs). First, as a greater effect of maternal than of paternal SDP is consistent with a causal effect of maternal SDP, we compared the effects of maternal and paternal SDP. Second, as a beneficial effect of quitting smoking before pregnancy is consistent with the causal effect, we compared the effects of SDP in mothers who quit smoking before pregnancy, and mothers who continued to smoke during pregnancy. All mothers were established smokers before their pregnancy. The results indicated a greater effect of maternal SDP, compared to paternal SDP, for externalizing, aggression, overactive and withdrawn behavior. Quitting smoking was associated with less externalizing, overactive behavior, aggression, and oppositional behavior, but had no effect on internalizing, anxious depression, or withdrawn behavior. We conclude that these results are consistent with a causal, but small, effect of smoking on externalizing problems at age 3. The results do not support a causal effect of maternal SDP on internalizing behaviors.

Heritability of compulsive Internet use in adolescents

Over the past decades, Internet use has grown substantially, and it now serves people as a supportive tool that is used regularly and—in large parts of the world—inevitably. Some people develop problematic Internet use, which may lead to addictive behavior and it is becoming important to explore the risk factors for compulsive Internet use. Data were analyzed on compulsive Internet use [with the Compulsive Internet Use Scale (CIUS)] from 5247 monozygotic (MZ) and dizygotic (DZ) adolescent twins registered with the Netherlands Twin Register. The participants form a sample that is informative for genetic analyses, allowing the investigation of the causes of individual differences in compulsive Internet use. The internal consistency of the instrument was high and the 1.6‐year test–retest correlation in a subsample (n = 902) was 0.55. CIUS scores increased slightly with age. Remarkably, gender did not explain variation in CIUS scores, as mean scores on the CIUS were the same in boys and girls. However, the time spent on specific Internet activities differed: boys spent more time on gaming, whereas girls spent more time on social network sites and chatting. The heritability estimates were the same for boys and girls: 48 percent of the individual differences in CIUS score were influenced by genetic factors. The remaining variance (52 percent) was due to environmental influences that were not shared between family members. Because a life without Internet is almost impossible nowadays, it is important to further explore the determinants of compulsive Internet use, including genetic risk factors.

Arithmetic, reading and writing performance has a strong genetic component: A study in primary school children

Even children attending the same primary school and taught by the same teacher differ greatly in their performance. In the Netherlands, performance at the end of primary school determines the enrollment in a particular level of secondary education. Identifying the impact of genes and the environment on individual differences in educational achievement between children is important. The Netherlands Twin Register has collected data on scores of tests used in primary school (ages 6 to 12) to monitor a child's educational progress in four domains, i.e. arithmetic, word reading, reading comprehension and spelling (1058 MZ and 1734 DZ twin pairs), and of a final test (2451 MZ and 4569 DZ twin pairs) in a large Dutch cohort. In general, individual differences in educational achievement were to a large extent due to genes and the influence of the family environment was negligible. Moreover, there is no evidence for gender differences in the underlying etiology.

No Evidence of Causal Effects of Blood Pressure on Cognition in the Population at Large

The large body of literature on the association between blood pressure (BP) and cognitive functioning has yielded mixed results, possibly due to the presence of non-linear effects across age, or because BP affects specific brain areas differently, impacting more on some cognitive skills than on others. If a robust association was detected among BP and specific cognitive tasks, the causal nature of reported associations between BP and cognition could be investigated in twin data, which allow a test of alternative explanations, including genetic pleiotropy. The present study first examines the association between BP and cognition in a sample of 1,140 participants with an age range between 10 and 86 years. Linear and quadratic effects of systolic BP (SBP) and diastolic BP (DBP) on cognitive functioning were examined for 17 tests across five functions. Associations were corrected for effects of sex and linear and quadratic effects of age. Second, to test a causal model, data from 123 monozygotic (MZ) twin pairs were analyzed to test whether cognitive functioning of the twins with the higher BP was different from that of the co-twins with lower BP. Associations between BP and cognitive functioning were absent for the majority of the cognitive tests, with the exception of a lower speed of emotion identification and verbal reasoning in subjects with high diastolic BP. In the MZ twin pair analyses, no effects of BP on cognition were found. We conclude that in the population at large, BP level is not associated with cognitive functioning in a clinically meaningful way.

Twin-sibling study and meta-analysis on the heritability of maximal oxygen consumption

Large individual differences exist in aerobic fitness in childhood and adolescence, but the relative contribution of genetic factors to this variation remains to be established. In a sample of adolescent twins and siblings (n = 479), heart rate (HR) and maximal oxygen uptake (V̇o2max) were recorded during the climax of a graded maximal exercise test. In addition, V̇o2max was predicted in two graded submaximal exercise tests on the cycle ergometer and the treadmill, using extrapolation of the HR/V̇o2 curve to the predicted HRmax. Heritability estimates for measured V̇o2max were 60% in ml/min and 55% for V̇o2max in ml·min(-1)·kg(-1). Phenotypic correlations between measured V̇o2max and predicted V̇o2max from either submaximal treadmill or cycle ergometer tests were modest (0.57 < r < 0.70), in part because of the poor agreement between predicted and actual HRmax. The majority of this correlation was explained by genetic factors; therefore, the submaximal exercise tests still led to very comparable estimates of heritability of V̇o2max. To arrive at a robust estimate for the heritability of V̇o2max in children to young adults, a sample size weighted meta-analysis was performed on all extant twin and sibling studies in this age range. Eight studies, including the current study, were meta-analyzed and resulted in a weighted heritability estimate of 59% (ml/min) and 72% (ml·min(-1)·kg(-1)) for V̇o2max. Taken together, the twin-sibling study and meta-analyses showed that from childhood to early adulthood genetic factors determine more than half of the individual differences in V̇o2max.

Epigenome-Wide Association Study of Wellbeing

Wellbeing (WB) is a major topic of research across several scientific disciplines, partly driven by its strong association with psychological and mental health. Twin-family studies have found that both genotype and environment play an important role in explaining the variance in WB. Epigenetic mechanisms, such as DNA methylation, regulate gene expression, and may mediate genetic and environmental effects on WB. Here, for the first time, we apply an epigenome-wide association study (EWAS) approach to identify differentially methylated sites associated with individual differences in WB. Subjects were part of the longitudinal survey studies of the Netherlands Twin Register (NTR) and participated in the NTR biobank project between 2002 and 2011. WB was assessed by a short inventory that measures satisfaction with life (SAT). DNA methylation was measured in whole blood by the Illumina Infinium HumanMethylation450 BeadChip (HM450k array) and the association between WB and DNA methylation level was tested at 411,169 autosomal sites. Two sites (cg10845147, p = 1.51 * 10−8 and cg01940273, p = 2.34 * 10−8) reached genome-wide significance following Bonferonni correction. Four more sites (cg03329539, p = 2.76* 10−7; cg09716613, p = 3.23 * 10−7; cg04387347, p = 3.95 * 10−7; and cg02290168, p = 5.23 * 10−7) were considered to be genome-wide significant when applying the widely used criterion of a FDR q value < 0.05. Gene ontology (GO) analysis highlighted enrichment of several central nervous system categories among higher-ranking methylation sites. Overall, these results provide a first insight into the epigenetic mechanisms associated with WB and lay the foundations for future work aiming to unravel the biological mechanisms underlying a complex trait like WB.

Development of the brain's structural network efficiency in early adolescence: A longitudinal DTI twin study

The brain is a network and our intelligence depends in part on the efficiency of this network. The network of adolescents differs from that of adults suggesting developmental changes. However, whether the network changes over time at the individual level and, if so, how this relates to intelligence, is unresolved in adolescence. In addition, the influence of genetic factors in the developing network is not known. Therefore, in a longitudinal study of 162 healthy adolescent twins and their siblings (mean age at baseline 9.9 [range 9.0-15.0] years), we mapped local and global structural network efficiency of cerebral fiber pathways (weighted with mean FA and streamline count) and assessed intelligence over a three-year interval. We find that the efficiency of the brain's structural network is highly heritable (locally up to 74%). FA-based local and global efficiency increases during early adolescence. Streamline count based local efficiency both increases and decreases, and global efficiency reorganizes to a net decrease. Local FA-based efficiency was correlated to IQ. Moreover, increases in FA-based network efficiency (global and local) and decreases in streamline count based local efficiency are related to increases in intellectual functioning. Individual changes in intelligence and local FA-based efficiency appear to go hand in hand in frontal and temporal areas. More widespread local decreases in streamline count based efficiency (frontal cingulate and occipital) are correlated with increases in intelligence. We conclude that the teenage brain is a network in progress in which individual differences in maturation relate to level of intellectual functioning.

Smoking During Adolescence as a Risk Factor for Attention Problems

BACKGROUND

Cigarette smoking and attention-deficit/hyperactivity disorder (ADHD) are highly comorbid. One explanation is that individuals with ADHD use cigarettes as "self-medication" to alleviate their attention problems. However, animal studies reported that exposure to nicotine during adolescence influences the developing brain and negatively affects attention. This is the first human study exploring the effects of smoking during adolescence on attention problems.

METHODS

Longitudinal data on smoking and attention problems were available for 1987 adult and 648 adolescent monozygotic twin pairs from the Netherlands Twin Register. Twin pairs were classified as concordant/discordant for smoking and compared on attention problems. Within adult discordant pairs, the difference in attention problems between the smoking and never-smoking twins was first assessed cross-sectionally. In longitudinal analyses, the increase in attention problems from adolescence, when neither twin smoked, to adulthood was compared within discordant pairs. In subgroups with longitudinal data from childhood and adolescence, changes in smoking concordance and subsequent changes in attention problems were explored.

RESULTS

Adult twins who ever smoked reported significantly more attention problems than their never-smoking co-twin. Longitudinal analyses showed a larger increase in attention problems from adolescence to adulthood in smoking twins than their never-smoking co-twin (p < .05). In childhood and adolescence, smoking twins had more attention problems than their never-smoking co-twin, whereas scores were similar before smoking was initiated or after both twins started smoking (not significant in all groups).

CONCLUSIONS

Results from this genetically informative study suggest smoking during adolescence leads to higher attention problem scores, lasting into adulthood.

Analysis of Behavioral and Emotional Problems in Children Highlights the Role of Genotype × Environment Interaction

This study tested for Genotype × Environment (G × E) interaction on behavioral and emotional problems in children using new methods that do not require identification of candidate genes or environments, can distinguish between interaction with shared and unique environment, and are insensitive to scale effects. Parental ratings of problem behavior from 14,755 twin pairs (5.3 years, SD = 0.22) indicated G × E interaction on emotional liability, social isolation, aggression, attention problems, dependency, anxiety, and physical coordination. Environmental influences increased in children who were genetically more predisposed to problem behavior, with ~20% of the variance due to G × E interaction (8% for anxiety to 37% for attention problems). Ignoring G × E interaction does not greatly bias heritability estimates, but it does offer a comprehensive model of the etiology for childhood problems.

Epigenome-Wide Association Study of Aggressive Behavior

Aggressive behavior is highly heritable, while environmental influences, particularly early in life, are also important. Epigenetic mechanisms, such as DNA methylation, regulate gene expression throughout development and adulthood, and may mediate genetic and environmental effects on complex traits. We performed an epigenome-wide association study (EWAS) to identify regions in the genome where DNA methylation level is associated with aggressive behavior. Subjects took part in longitudinal survey studies from the Netherlands Twin Register (NTR) and participated in the NTR biobank project between 2004 and 2011 (N = 2,029, mean age at blood sampling = 36.4 years, SD = 12.4, females = 69.2%). Aggressive behavior was rated with the ASEBA Adult Self-Report (ASR). DNA methylation was measured in whole blood by the Illumina HM450k array. The association between aggressive behavior and DNA methylation level at 411,169 autosomal sites was tested. Association analyses in the entire cohort showed top sites at cg01792876 (chr8; 116,684,801, nearest gene = TRPS1, p = 7.6 × 10−7, False discovery rate (FDR) = 0.18) and cg06092953 (chr18; 77,905,699, nearest gene = PARD6G-AS1, p = 9.0 ×10−7, FDR = 0.18). Next, we compared methylation levels in 20 pairs of monozygotic (MZ) twins highly discordant for aggression. Here the top sites were cg21557159 (chr 11; 107,795,699, nearest gene = RAB39, p = 5.7 × 10−6, FDR = 0.99), cg08648367 (chr 19; 51,925,472, nearest gene = SIGLEC10, p = 7.6 × 10−6, FDR = 0.99), and cg14212412 (chr 6; 105,918,992, nearest gene = PREP, p = 8.0 × 10−6, FDR = 0.99). The two top hits based on the entire cohort showed the same direction of effect in discordant MZ pairs (cg01792876, Pdiscordant twins = 0.09 and cg06092953, Pdiscordant twins = 0.24). The other way around, two of the three most significant sites in discordant MZ pairs showed the same direction of effect in the entire cohort (cg08648367, Pentire EWAS = 0.59 and cg14212412, Pentire EWAS = 3.1 × 10−3). Gene ontology analysis highlighted significant enrichment of various central nervous system categories among higher-ranking methylation sites. Higher-ranking methylation sites also showed enrichment for DNase I hypersensitive sites and promoter regions, showing that DNA methylation in peripheral tissues is likely to be associated with aggressive behavior.

Epigenome-Wide Association Study of Tic Disorders

Tic disorders are moderately heritable common psychiatric disorders that can be highly troubling, both in childhood and in adulthood. In this study, we report results obtained in the first epigenome-wide association study (EWAS) of tic disorders. The subjects are participants in surveys at the Netherlands Twin Register (NTR) and the NTR biobank project. Tic disorders were measured with a self-report version of the Yale Global Tic Severity Scale Abbreviated version (YGTSS-ABBR), included in the 8th wave NTR data collection (2008). DNA methylation data consisted of 411,169 autosomal methylation sites assessed by the Illumina Infinium HumanMethylation450 BeadChip Kit (HM450k array). Phenotype and DNA methylation data were available in 1,678 subjects (mean age = 41.5). No probes reached genome-wide significance (p < 1.2 × 10(-7)). The strongest associated probe was cg15583738, located in an intergenic region on chromosome 8 (p = 1.98 × 10(-6)). Several of the top ranking probes (p < 1 × 10(-4)) were in or nearby genes previously associated with neurological disorders (e.g., GABBRI, BLM, and ADAM10), warranting their further investigation in relation to tic disorders. The top significantly enriched gene ontology (GO) terms among higher ranking methylation sites included anatomical structure morphogenesis (GO:0009653, p = 4.6 × 10-(15)) developmental process (GO:0032502, p = 2.96 × 10(-12)), and cellular developmental process (GO:0048869, p = 1.96 × 10(-12)). Overall, these results provide a first insight into the epigenetic mechanisms of tic disorders. This first study assesses the role of DNA methylation in tic disorders, and it lays the foundations for future work aiming to unravel the biological mechanisms underlying the architecture of this disorder.

Differences in Adolescent Physical Fitness: A Multivariate Approach and Meta-analysis

Physical fitness can be defined as a set of components that determine exercise ability and influence performance in sports. This study investigates the genetic and environmental influences on individual differences in explosive leg strength (vertical jump), handgrip strength, balance, and flexibility (sit-and-reach) in 227 healthy monozygotic and dizygotic twin pairs and 38 of their singleton siblings (mean age 17.2 ± 1.2). Heritability estimates were 49 % (95 % CI 35–60 %) for vertical jump, 59 % (95 % CI 46–69 %) for handgrip strength, 38 % (95 % CI 22–52 %) for balance, and 77 % (95 % CI 69–83 %) for flexibility. In addition, a meta-analysis was performed on all twin studies in children, adolescents and young adults reporting heritability estimates for these phenotypes. Fifteen studies, including results from our own study, were meta-analyzed by computing the weighted average heritability. This showed that genetic factors explained most of the variance in vertical jump (62 %; 95 % CI 47–77 %, N = 874), handgrip strength (63 %; 95 % CI 47–73 %, N = 4516) and flexibility (50 %; 95 % CI 38–61 %, N = 1130) in children and young adults. For balance this was 35 % (95 % CI 19–51 %, N = 978). Finally, multivariate modeling showed that the phenotypic correlations between the phenotypes in current study (0.07 < r < 0.27) were mostly driven by genetic factors. It is concluded that genetic factors contribute significantly to the variance in muscle strength, flexibility and balance; factors that may play a key role in the individual differences in adolescent exercise ability and sports performance.

Evidence for Gender-Dependent Genotype by Environment Interaction in Adult Depression

Depression in adults is heritable with about 40 % of the phenotypic variance due to additive genetic effects and the remaining phenotypic variance due to unique (unshared) environmental effects. Common environmental effects shared by family members are rarely found in adults. One possible explanation for this finding is that there is an interaction between genes and the environment which may mask effects of the common environment. To test this hypothesis, we investigated genotype by environment interaction in a large sample of female and male adult twins aged 18-70 years. The anxious depression subscale of the Adult Self Report from the Achenbach System of Empirically Based Assessment (Achenbach and Rescorla in Manual for the ASEBA adult: forms and profiles, 2003) was completed by 13,022 twins who participate in longitudinal studies of the Netherlands Twin Register. In a single group analysis, we found genotype by unique environment interaction, but no genotype by common environment interaction. However, when conditioning on gender, we observed genotype by common environment interaction in men, with larger common environmental variance in men who are genetically less at risk to develop depression. Although the effect size of the interaction is characterized by large uncertainty, the results show that there is at least some variance due to the common environment in adult depression in men.

Chorionicity and Heritability Estimates from Twin Studies: The Prenatal Environment of Twins and Their Resemblance Across a Large Number of Traits

There are three types of monozygotic (MZ) twins. MZ twins can either share one chorion and one amnion, each twin can have its own amnion, or MZ twins can—like dizygotic twins—each have their own chorion and amnion. Sharing the same chorion may create a more similar/dissimilar prenatal environment and bias heritability estimates, but most twin studies do not distinguish between these three types of MZ twin pairs. The aim of this paper is to investigate the effect of chorion sharing on the similarity within MZ twin pairs for a large number of traits. Information on chorion status was obtained for the Netherlands twin register (NTR) by linkage to the records from the database of the dutch pathological anatomy national automated archive (PALGA). Record linkage was successful for over 9000 pairs. Effect of chorion type was tested by comparing the within-pair similarity between monochorionic (MC) and dichorionic (DC) MZ twins on 66 traits including weight, height, motor milestones, child problem behaviors, cognitive function, wellbeing and personality. For only 10 traits, within-pair similarity differed between MCMZ and DCMZ pairs. For traits influenced by birth weight (e.g. weight and height in young children) we expected that MC twins would be more discordant. This was found for 5 out of 13 measures. When looking at traits where blood supply is important, we saw MCMZ twins to be more concordant than DCMZ’s for 3 traits. We conclude that the influence on the MZ twin correlation of the intra-uterine prenatal environment, as measured by sharing a chorion type, is small and limited to a few phenotypes. This implies that the assumption of equal prenatal environment of mono- and DC MZ twins, which characterizes the classical twin design, is largely tenable.

Regular exercise behaviour in youth is not related to current body mass index or body mass index at 7-year follow-up

Objective

This population‐based study aimed (1) to test the presence of an association between regular voluntary exercise behaviour (EB) that is performed in leisure time and body mass index (BMI) in youth and (2) to investigate the causal nature of this association using a longitudinal design in genetically informative subjects.

Design and methods

Both EB and BMI were assessed repeatedly over time in 21 458 twin individuals from the Netherlands Twin Register (47.5% male) – first by parental report (ages 7, 10 and 12) and subsequently through self‐report surveys (ages 14, 16 and 18). EB was quantified as weekly metabolic equivalent of task hours.

Results

Correlations over time were higher for BMI than for EB (r ≈ 0.70 vs. r ≈ 0.35) across 12 different follow‐up periods. Cross‐sectionally, regular involvement in EB was not associated with lower BMI in childhood and in genetically identical twin pairs discordant for EB; the exercising twin did not have a lower BMI than the non‐exercising twin. Longitudinally, linear and quadratic relationships between EB and BMI were non‐significant. Changes in EB over time did not induce opposite changes in BMI.

Conclusions

No consistent association between regular EB and BMI was observed from ages 7 to 18 years.

Single Nucleotide Polymorphism Heritability of Behavior Problems in Childhood: Genome-Wide Complex Trait Analysis

OBJECTIVE

Genetic factors contribute to individual differences in behavior problems. In children, genome-wide association studies (GWAS) have yielded the first suggestive results when aiming to identify genetic variants that explain heritability, but the proportion of genetic variance that can be attributed to common single nucleotide polymorphisms (SNPs) remains to be determined, as only a few studies have estimated SNP heritability, with diverging results.

METHOD

Genomic-relationship-matrix restricted maximum likelihood (GREML) as implemented in the software Genome-Wide Complex Trait Analysis (GCTA) was used to estimate SNP heritability (SNP h(2)) for multiple phenotypes within 4 broad domains of children's behavioral problems (attention-deficit/hyperactivity symptoms, internalizing, externalizing, and pervasive developmental problems) and cognitive function. We combined phenotype and genotype data from 2 independent, population-based Dutch cohorts, yielding a total number of 1,495 to 3,175 of 3-, 7-, and 9-year-old children.

RESULTS

Significant SNP heritability estimates were found for attention-deficit/hyperactivity symptoms (SNP h(2) = 0.37-0.71), externalizing problems (SNP h(2) = 0.44), and total problems (SNP h(2) = 0.18), rated by mother or teacher. Sensitivity analyses with exclusion of extreme cases and quantile normalization of the phenotype data decreased SNP h(2) as expected under genetic inheritance, but they remained statistically significant for most phenotypes.

CONCLUSION

We provide evidence of the influence of common SNPs on child behavior problems in an ethnically homogenous sample. These results support the continuation of large GWAS collaborative efforts to unravel the genetic basis of complex child behaviors.

The Influence of Chorion Type on Health Measures at Birth and Dental Development in Australian and Dutch Twins: A Comparative Study

Chorion type may significantly influence the prenatal environment of twins. This study explored the associations between chorion type and gestational age, birth weight, birth length, and the timing of emergence of the first primary tooth in two populations of twins, Australian and Dutch. Additionally, we investigated the relationship between chorion type and birth weight discordance (BWD) in order to determine whether a significant relationship existed between discordance in birth weight and discordance in the timing of emergence of the first primary tooth. The two study samples consisted of 409 Australian twin pairs and 301 Dutch twin pairs, all of European ancestry. Data were collected through a combination of questionnaires and recording charts administered to the parents and through linkage with biological databases. In the Australian sample, monozygotic monochorionic (MZMC) twins experienced the shortest mean gestation time (35 weeks), the lowest mean birth length (46 cm) and the lowest mean birth weight (2.3 kg) compared with other twin groups. For the same variables in the Dutch sample, these trends with MZMC twinning were not observed. Chorion type did not significantly affect the mean timing of emergence of the first primary tooth in either sample. Monochorionicity was found to be significantly associated with BWD in both samples, but there was a significant association between BWD in MZMC twin pairs and timing of emergence of the first primary tooth only in the Australian sample. Results from this study support previous findings that the timing of emergence of the first primary tooth is influenced strongly by genetic factors and is well protected from environmental disturbances.

Estimation of Genetic Relationships Between Individuals Across Cohorts and Platforms: Application to Childhood Height

Combining genotype data across cohorts increases power to estimate the heritability due to common single nucleotide polymorphisms (SNPs), based on analyzing a Genetic Relationship Matrix (GRM). However, the combination of SNP data across multiple cohorts may lead to stratification, when for example, different genotyping platforms are used. In the current study, we address issues of combining SNP data from different cohorts, the Netherlands Twin Register (NTR) and the Generation R (GENR) study. Both cohorts include children of Northern European Dutch background (N = 3102 + 2826, respectively) who were genotyped on different platforms. We explore imputation and phasing as a tool and compare three GRM-building strategies, when data from two cohorts are (1) just combined, (2) pre-combined and cross-platform imputed and (3) cross-platform imputed and post-combined. We test these three strategies with data on childhood height for unrelated individuals (N = 3124, average age 6.7 years) to explore their effect on SNP-heritability estimates and compare results to those obtained from the independent studies. All combination strategies result in SNP-heritability estimates with a standard error smaller than those of the independent studies. We did not observe significant difference in estimates of SNP-heritability based on various cross-platform imputed GRMs. SNP-heritability of childhood height was on average estimated as 0.50 (SE = 0.10). Introducing cohort as a covariate resulted in ≈2 % drop. Principal components (PCs) adjustment resulted in SNP-heritability estimates of about 0.39 (SE = 0.11). Strikingly, we did not find significant difference between cross-platform imputed and combined GRMs. All estimates were significant regardless the use of PCs adjustment. Based on these analyses we conclude that imputation with a reference set helps to increase power to estimate SNP-heritability by combining cohorts of the same ethnicity genotyped on different platforms. However, important factors should be taken into account such as remaining cohort stratification after imputation and/or phenotypic heterogeneity between and within cohorts. Whether one should use imputation, or just combine the genotype data, depends on the number of overlapping SNPs in relation to the total number of genotyped SNPs for both cohorts, and their ability to tag all the genetic variance related to the specific trait of interest.