Strong associations between neural tube defects (NTDs) and monozygotic (MZ) twinning have long been noted, and it has been suggested that NTD cases who do not present as MZ twins may be the survivors of MZ twinning events. We have recently shown that MZ twins carry a strong, distinctive DNA methylation signature and have developed an algorithm based on genomewide DNA methylation array data that distinguishes MZ twins from dizygotic twins and other relatives at well above chance level. We have applied this algorithm to published methylation data from five fetal tissues (placental chorionic villi, kidney, spinal cord, brain and muscle) collected from spina bifida cases (n = 22), anencephalic cases (n = 15) and controls (n = 19). We see no difference in signature between cases and controls, providing no support for a common etiological role of MZ twinning in NTDs. The strong associations therefore continue to await elucidation.

There are research questions whose answers require record linkage of multiple databases that may be characterized by limited options for full data sharing. For this purpose, the Open Data Infrastructure for Social Science and Economic Innovations (ODISSEI) consortium has supported the development of the ODISSEI Secure Supercomputer (OSSC) platform that allows researchers to link cohort data to data from Statistics Netherlands and run large-scale analyses in a high-performance computing (HPC) environment. Here, we report a successful record linkage genomewide association (GWA) study on expenditure for total health, mental health, primary and hospital care, and medication. Record linkage for genotype data from 16,726 participants from the Netherlands Twin Register (NTR) with data from Statistics Netherlands was accomplished in the secure OSSC platform, followed by gene-based tests and estimation of total and single nucleotide polymorphism (SNP)-based heritability. The total heritability of expenditure ranged between 29.4% (SE 0.8) and 37.5% (SE 0.8), but GWA analyses did not identify SNPs or genes that were genomewide significantly associated with health care expenditure. SNP-based heritability was between 0.0% (SE 3.5) and 5.4% (SE 4.0) and was different from zero for mental health care and primary care expenditure. We conclude that successfully linking genotype data to administrative health care expenditure data from Statistics Netherlands is feasible and demonstrates a series of analyses on health care expenditure. The OSSC platform offers secure possibilities for analyzing linked data in large scale and realizing sample sizes required for GWA studies, providing invaluable opportunities to answer many new research questions.

The COllaborative project of Development of Anthropometrical measures in Twins (CODATwins) project is a large international collaborative effort to analyze individual-level phenotype data from twins in multiple cohorts from different environments. The main objective is to study factors that modify genetic and environmental variation of height, body mass index (BMI, kg/m2) and size at birth, and additionally to address other research questions such as long-term consequences of birth size. The project started in 2013 and is open to all twin projects in the world having height and weight measures on twins with information on zygosity. Thus far, 54 twin projects from 24 countries have provided individual-level data. The CODATwins database includes 489,981 twin individuals (228,635 complete twin pairs). Since many twin cohorts have collected longitudinal data, there is a total of 1,049,785 height and weight observations. For many cohorts, we also have information on birth weight and length, own smoking behavior and own or parental education. We found that the heritability estimates of height and BMI systematically changed from infancy to old age. Remarkably, only minor differences in the heritability estimates were found across cultural–geographic regions, measurement time and birth cohort for height and BMI. In addition to genetic epidemiological studies, we looked at associations of height and BMI with education, birth weight and smoking status. Within-family analyses examined differences within same-sex and opposite-sex dizygotic twins in birth size and later development. The CODATwins project demonstrates the feasibility and value of international collaboration to address gene-by-exposure interactions that require large sample sizes and address the effects of different exposures across time, geographical regions and socioeconomic status.

Whether monozygotic (MZ) and dizygotic (DZ) twins differ from each other in a variety of phenotypes is important for genetic twin modeling and for inferences made from twin studies in general. We analyzed whether there were differences in individual, maternal and paternal education between MZ and DZ twins in a large pooled dataset. Information was gathered on individual education for 218,362 adult twins from 27 twin cohorts (53% females; 39% MZ twins), and on maternal and paternal education for 147,315 and 143,056 twins respectively, from 28 twin cohorts (52% females; 38% MZ twins). Together, we had information on individual or parental education from 42 twin cohorts representing 19 countries. The original education classifications were transformed to education years and analyzed using linear regression models. Overall, MZ males had 0.26 (95% CI [0.21, 0.31]) years and MZ females 0.17 (95% CI [0.12, 0.21]) years longer education than DZ twins. The zygosity difference became smaller in more recent birth cohorts for both males and females. Parental education was somewhat longer for fathers of DZ twins in cohorts born in 1990–1999 (0.16 years, 95% CI [0.08, 0.25]) and 2000 or later (0.11 years, 95% CI [0.00, 0.22]), compared with fathers of MZ twins. The results show that the years of both individual and parental education are largely similar in MZ and DZ twins. We suggest that the socio-economic differences between MZ and DZ twins are so small that inferences based upon genetic modeling of twin data are not affected.

The monocyte–lymphocyte ratio (MLR) is a useful biomarker for disease development, but little is known about the extent to which genetic and environmental factors influence MLR variation. Here, we study the genetic architecture of MLR and determine the influence of demographic and lifestyle factors on MLR in data from a Dutch non-patient twin-family population. Data were obtained in 9,501 individuals from the Netherlands Twin Register. We used regression analyses to determine the effects of age, sex, smoking, and body mass index (BMI) on MLR and its subcomponents. Data on twins, siblings and parents (N = 7,513) were analyzed by genetic structural equation modeling to establish heritability and genome wide single nucleotide polymorphism (SNP) data from a genotyped subsample (N = 5,892) and used to estimate heritability explained by SNPs. SNP and phenotype data were also analyzed in a genome-wide association study to identify the genes involved in MLR. Linkage disequilibrium (LD) score regression and expression quantitative trait loci (eQTL) analyses were performed to further explore the genetic findings. Results showed that age, sex, and age × sex interaction effects were present for MLR and its subcomponents. Variation in MLR was not related to BMI, but smoking was positively associated with MLR. Heritability was estimated at 40% for MLR, 58% for monocyte, and 58% for lymphocyte count. The Genome-wide association study (GWAS) identified a locus on ITGA4 that was associated with MLR and only marginally significantly associated with monocyte count. For monocyte count, additional genetic variants were identified on ITPR3, LPAP1, and IRF8. For lymphocyte count, GWAS provided no significant findings. Taking all measured SNPs together, their effects accounted for 13% of the heritability of MLR, while all known and identified genetic loci explained 1.3% of variation in MLR. eQTL analyses showed that these genetic variants were unlikely to be eQTLs. In conclusion, variation in MLR level in the general population is heritable and influenced by age, sex, and smoking. We identified gene variants in the ITGA4 gene associated with variation in MLR. The significant SNP-heritability indicates that more genetic variants are likely to be involved.

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.

Approximately half of the variation in wellbeing measures overlaps with variation in personality traits. Studies of non-human primate pedigrees and human twins suggest that this is due to common genetic influences. We tested whether personality polygenic scores for the NEO Five-Factor Inventory (NEO-FFI) domains and for item response theory (IRT) derived extraversion and neuroticism scores predict variance in wellbeing measures. Polygenic scores were based on published genome-wide association (GWA) results in over 17,000 individuals for the NEO-FFI and in over 63,000 for the IRT extraversion and neuroticism traits. The NEO-FFI polygenic scores were used to predict life satisfaction in 7 cohorts, positive affect in 12 cohorts, and general wellbeing in 1 cohort (maximal N = 46,508). Meta-analysis of these results showed no significant association between NEO-FFI personality polygenic scores and the wellbeing measures. IRT extraversion and neuroticism polygenic scores were used to predict life satisfaction and positive affect in almost 37,000 individuals from UK Biobank. Significant positive associations (effect sizes <0.05%) were observed between the extraversion polygenic score and wellbeing measures, and a negative association was observed between the polygenic neuroticism score and life satisfaction. Furthermore, using GWA data, genetic correlations of -0.49 and -0.55 were estimated between neuroticism with life satisfaction and positive affect, respectively. The moderate genetic correlation between neuroticism and wellbeing is in line with twin research showing that genetic influences on wellbeing are also shared with other independent personality domains.

We analyzed birth order differences in means and variances of height and body mass index (BMI) in monozygotic (MZ) and dizygotic (DZ) twins from infancy to old age. The data were derived from the international CODATwins database. The total number of height and BMI measures from 0.5 to 79.5 years of age was 397,466. As expected, first-born twins had greater birth weight than second-born twins. With respect to height, first-born twins were slightly taller than second-born twins in childhood. After adjusting the results for birth weight, the birth order differences decreased and were no longer statistically significant. First-born twins had greater BMI than the second-born twins over childhood and adolescence. After adjusting the results for birth weight, birth order was still associated with BMI until 12 years of age. No interaction effect between birth order and zygosity was found. Only limited evidence was found that birth order influenced variances of height or BMI. The results were similar among boys and girls and also in MZ and DZ twins. Overall, the differences in height and BMI between first- and second-born twins were modest even in early childhood, while adjustment for birth weight reduced the birth order differences but did not remove them for BMI.

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.

Low birth weight (LBW) can have an impact on health outcomes in later life, especially in relation to pre-disposition to metabolic disease. Several studies suggest that LBW resulting from restricted intrauterine growth leaves a footprint on DNA methylation in utero, and this influence likely persists into adulthood. To investigate this further, we performed epigenome-wide association analyses of blood DNA methylation using Infinium HumanMethylation450 BeadChip profiles in 71 adult monozygotic (MZ) twin pairs who were extremely discordant for birth weight. A signal mapping to the IGF1R gene (cg12562232, p = 2.62 × 10−8), was significantly associated with birth weight discordance at a genome-wide false-discovery rate (FDR) of 0.05. We pursued replication in three additional independent datasets of birth weight discordant MZ pairs and observed the same direction of association, but the results were not significant. However, a meta-analysis across the four independent samples, in total 216 birth-weight discordant MZ twin pairs, showed a significant positive association between birth weight and DNA methylation differences at IGF1R (random-effects meta-analysis p = .04), and the effect was particularly pronounced in older twins (random-effects meta-analysis p = .008, 98 older birth-weight discordant MZ twin pairs). The results suggest that severe intra-uterine growth differences (birth weight discordance >20%) are associated with methylation changes in the IGF1R gene in adulthood, independent of genetic effects.

A trend toward greater body size in dizygotic (DZ) than in monozygotic (MZ) twins has been suggested by some but not all studies, and this difference may also vary by age. We analyzed zygosity differences in mean values and variances of height and body mass index (BMI) among male and female twins from infancy to old age. Data were derived from an international database of 54 twin cohorts participating in the COllaborative project of Development of Anthropometrical measures in Twins (CODATwins), and included 842,951 height and BMI measurements from twins aged 1 to 102 years. The results showed that DZ twins were consistently taller than MZ twins, with differences of up to 2.0 cm in childhood and adolescence and up to 0.9 cm in adulthood. Similarly, a greater mean BMI of up to 0.3 kg/m2 in childhood and adolescence and up to 0.2 kg/m2 in adulthood was observed in DZ twins, although the pattern was less consistent. DZ twins presented up to 1.7% greater height and 1.9% greater BMI than MZ twins; these percentage differences were largest in middle and late childhood and decreased with age in both sexes. The variance of height was similar in MZ and DZ twins at most ages. In contrast, the variance of BMI was significantly higher in DZ than in MZ twins, particularly in childhood. In conclusion, DZ twins were generally taller and had greater BMI than MZ twins, but the differences decreased with age in both sexes.

For over 100 years, the genetics of human anthropometric traits has attracted scientific interest. In particular, height and body mass index (BMI, calculated as kg/m2) have been under intensive genetic research. However, it is still largely unknown whether and how heritability estimates vary between human populations. Opportunities to address this question have increased recently because of the establishment of many new twin cohorts and the increasing accumulation of data in established twin cohorts. We started a new research project to analyze systematically (1) the variation of heritability estimates of height, BMI and their trajectories over the life course between birth cohorts, ethnicities and countries, and (2) to study the effects of birth-related factors, education and smoking on these anthropometric traits and whether these effects vary between twin cohorts. We identified 67 twin projects, including both monozygotic (MZ) and dizygotic (DZ) twins, using various sources. We asked for individual level data on height and weight including repeated measurements, birth related traits, background variables, education and smoking. By the end of 2014, 48 projects participated. Together, we have 893,458 height and weight measures (52% females) from 434,723 twin individuals, including 201,192 complete twin pairs (40% monozygotic, 40% same-sex dizygotic and 20% opposite-sex dizygotic) representing 22 countries. This project demonstrates that large-scale international twin studies are feasible and can promote the use of existing data for novel research purposes.