Genetic Influence on Social Support: A Twin Study.

Social support is often considered an environmental factor affecting health, especially in aging populations. However, its genetic underpinnings suggest a more complex origin. This study investigates the heritability of social support through applying a threshold model on data of a large adult sample of twins (N = 8019) from the Netherlands Twin Register, collected between 2009 and 2011. The study employed the Duke - UNC Functional Social Support Questionnaire to assess social support quality. Our analysis revealed genetic contributions to social support, with heritability estimated at 37%, without a contribution of shared environment and no differences between men and women in heritability. The study's results underscore the complexity of social support as a trait influenced by genetic and environmental factors, challenging the notion that it is solely an environmental construct.

Development of a Diabetes Dietary Quality Index: Reproducibility and Associations with Measures of Insulin Resistance, Beta Cell Function, and Hyperglycemia.

AIMS: Various dietary risk factors for type 2 diabetes have been identified. A short assessment of dietary patterns related to the risk for type 2 diabetes mellitus may be relevant in clinical practice given the largely preventable nature of the disease. The aim of this study was to investigate the reproducibility of a short food frequency questionnaire based on available knowledge of diabetes-related healthy diets. In addition, we aimed to investigate whether a Diabetes Dietary Quality Index based on this questionnaire was related to metabolic risk factors, including measures of beta cell function and insulin sensitivity. METHODS: A short food frequency questionnaire was composed by selecting fourteen questions (representing eight dietary factors) from existing food frequency questionnaires on the basis of their reported relationship with diabetes risk. Healthy participants (N = 176) from a Dutch family study completed the questionnaire and a subgroup (N = 123) completed the questionnaire twice. Reproducible items from the short questionnaire were combined into an index. The association between the Diabetes Dietary Quality index and metabolic risk factors was investigated using multiple linear regression analysis. Measures of beta cell function and insulin sensitivity were derived from a mixed meal test and an euglycemic-hyperinsulinemic and modified hyperglycemic clamp test. RESULTS: Our results show that this new short food frequency questionnaire is reliable (Intraclass Correlations ranged between 0.5 and 0.9). A higher Diabetes Dietary Quality index score was associated with lower 2 h post-meal glucose (ß -0.02, SE 0.006, p < 0.05), HbA1c (ß -0.07, SE 0.02, p < 0.05), total cholesterol, (ß -0.02, SE 0.07, p < 0.05), LDL cholesterol, (ß -0.19, SE 0.07, p < 0.05), fasting (ß -0.4, SE 0.16, p < 0.05) and post-load insulin, (ß -3.9, SE 1.40, p < 0.05) concentrations and the incremental AUC of glucose during MMT (ß -1.9, SE 0.97, p < 0.05). The scores obtained for the oral glucose insulin sensitivity-derived mixed meal test were higher in subjects who scored higher on the Diabetes Dietary Quality index (ß 0.89, 0.39, p < 0.05). In contrast, we found no significant associations between the Diabetes Dietary Quality index and clamp measures of beta cell function. CONCLUSIONS: We identified a questionnaire-derived Diabetes Dietary Quality index that was reproducible and inversely associated with a number of type 2 diabetes mellitus and metabolic risk factors, like 2 h post-meal glucose, Hba1c and LDL, and total cholesterol. Once relative validity has been established, the Diabetes Dietary Quality index could be used by health care professionals to identify individuals with diets adversely related to development of type 2 diabetes.

Genomic analysis of intracranial and subcortical brain volumes yields polygenic scores accounting for variation across ancestries.

Subcortical brain structures are involved in developmental, psychiatric and neurological disorders. Here we performed genome-wide association studies meta-analyses of intracranial and nine subcortical brain volumes (brainstem, caudate nucleus, putamen, hippocampus, globus pallidus, thalamus, nucleus accumbens, amygdala and the ventral diencephalon) in 74,898 participants of European ancestry. We identified 254 independent loci associated with these brain volumes, explaining up to 35% of phenotypic variance. We observed gene expression in specific neural cell types across differentiation time points, including genes involved in intracellular signaling and brain aging-related processes. Polygenic scores for brain volumes showed predictive ability when applied to individuals of diverse ancestries. We observed causal genetic effects of brain volumes with Parkinson's disease and attention-deficit/hyperactivity disorder. Findings implicate specific gene expression patterns in brain development and genetic variants in comorbid neuropsychiatric disorders, which could point to a brain substrate and region of action for risk genes implicated in brain diseases.

Genomic analysis of intracranial and subcortical brain volumes yields polygenic scores accounting for variation across ancestries.

Subcortical brain structures are involved in developmental, psychiatric and neurological disorders. We performed GWAS meta-analyses of intracranial and nine subcortical brain volumes (brainstem, caudate nucleus, putamen, hippocampus, globus pallidus, thalamus, nucleus accumbens, amygdala and, for the first time, the ventral diencephalon) in 74,898 participants of European ancestry. We identified 254 independent loci associated with these brain volumes, explaining up to 35% of phenotypic variance. We observed gene expression in specific neural cell types across differentiation time points, including genes involved in intracellular signalling and brain ageing-related processes. Polygenic scores for brain volumes showed predictive ability when applied to individuals of diverse ancestries. We observed causal genetic effects of brain volumes with Parkinson's disease and ADHD. Findings implicate specific gene expression patterns in brain development and genetic variants in comorbid neuropsychiatric disorders, which could point to a brain substrate and region of action for risk genes implicated in brain diseases.

Maternal Characteristics in Natural and Medically Assisted Reproduction Dizygotic Twin Pregnancies.

Previous studies have shown that mothers of naturally conceived dizygotic (DZ) twins tend to be taller, older, and smoke more than mothers of naturally conceived monozygotic (MZ) twin and mothers of singletons. Here, we investigate whether mothers of naturally conceived DZ twins differ from mothers who conceived their DZ twins after medically assisted reproduction (MAR) in eight maternal traits related to fertility based on observational survey data. We include data from 33,648 mothers from the Netherlands Twin Register (NTR) and 1660 mothers of twins from the Norwegian Mother, Father and Child Cohort Study (MoBA). We contrast mothers of naturally conceived DZ twins with mothers of MAR DZ twins. Next, we further segment the MAR group into mothers who underwent hormonal induction of ovulation but not in vitro fertilization (IVF) and those who IVF twins, comparing them both to each other and against the mothers of naturally conceived DZ twins. Mothers of naturally conceived DZ twins smoke more often, differ in body composition, have a higher maternal age and have more offspring before the twins than mothers of MZ twins. Compared to MAR DZ twin mothers, mothers of naturally conceived DZ twins have fewer miscarriages, lower maternal age and increased height, more offspring and are more often smokers. BMI before the twin pregnancy is similar in both natural and MAR DZ twin mothers. Mothers who received hormonal induction of ovulation (OI) have a lower maternal age, fewer miscarriages, and a higher number of offspring before their twin pregnancy than twin mothers who received IVF and/or intracytoplasmic sperm injection (ICSI) treatments. Our study shows that twin mothers are a heterogenous group and the differences between twin mothers should be taken into account in epidemiological and genetic research that includes twins.

Genome-wide meta-analysis of ascertainment and symptom structures of major depression in case-enriched and community cohorts.

BACKGROUND: Diagnostic criteria for major depressive disorder allow for heterogeneous symptom profiles but genetic analysis of major depressive symptoms has the potential to identify clinical and etiological subtypes. There are several challenges to integrating symptom data from genetically informative cohorts, such as sample size differences between clinical and community cohorts and various patterns of missing data. METHODS: We conducted genome-wide association studies of major depressive symptoms in three cohorts that were enriched for participants with a diagnosis of depression (Psychiatric Genomics Consortium, Australian Genetics of Depression Study, Generation Scotland) and three community cohorts who were not recruited on the basis of diagnosis (Avon Longitudinal Study of Parents and Children, Estonian Biobank, and UK Biobank). We fit a series of confirmatory factor models with factors that accounted for how symptom data was sampled and then compared alternative models with different symptom factors. RESULTS: The best fitting model had a distinct factor for Appetite/Weight symptoms and an additional measurement factor that accounted for the skip-structure in community cohorts (use of Depression and Anhedonia as gating symptoms). CONCLUSION: The results show the importance of assessing the directionality of symptoms (such as hypersomnia versus insomnia) and of accounting for study and measurement design when meta-analyzing genetic association data.

NMR metabolomics-guided DNA methylation mortality predictors.

BACKGROUND: 1H-NMR metabolomics and DNA methylation in blood are widely known biomarkers predicting age-related physiological decline and mortality yet exert mutually independent mortality and frailty signals. METHODS: Leveraging multi-omics data in four Dutch population studies (N = 5238, ∼40% of which male) we investigated whether the mortality signal captured by 1H-NMR metabolomics could guide the construction of DNA methylation-based mortality predictors. FINDINGS: We trained DNA methylation-based surrogates for 64 metabolomic analytes and found that analytes marking inflammation, fluid balance, or HDL/VLDL metabolism could be accurately reconstructed using DNA-methylation assays. Interestingly, a previously reported multi-analyte score indicating mortality risk (MetaboHealth) could also be accurately reconstructed. Sixteen of our derived surrogates, including the MetaboHealth surrogate, showed significant associations with mortality, independent of relevant covariates. INTERPRETATION: The addition of our metabolic analyte-derived surrogates to the well-established epigenetic clock GrimAge demonstrates that our surrogates potentially represent valuable mortality signal. FUNDING: BBMRI-NL, X-omics, VOILA, Medical Delta, NWO, ERC.

Epigenetic Signatures of Asthma: A Comprehensive Study of DNA Methylation and Clinical Markers.

Background: Asthma, a complex respiratory disease, presents with inflammatory symptoms in the lungs, blood, and other tissues. We investigated the relationship between DNA methylation and 35 clinical markers of asthma. The Illumina Infinium EPIC v1 methylation array was used to evaluate 742,442 CpGs in whole blood samples from 319 participants. They were part of the Netherlands Twin Register from families with at least one member suffering from severe asthma. Repeat blood samples were taken after 10 years from 182 of these individuals. Principal component analysis (PCA) on the clinical markers yielded ten principal components (PCs) that explained 92.8% of the total variance. We performed epigenome-wide association studies (EWAS) for each of the ten PCs correcting for familial structure and other covariates. Results: 221 unique CpGs reached genome-wide significance at timepoint 1 (T1) after Bonferroni correction. PC7 accounted for the majority of associations (204), which correlated with loadings of eosinophil counts and immunoglobulin levels. Enrichment analysis via the EWAS Atlas identified 190 of these CpGs to be previously identified in EWASs of asthma and asthma-related traits. Proximity assessment to previously identified SNPs associated with asthma identified 17 unique SNPs within 1 MB of two of the 221 CpGs. EWAS in 182 individuals with epigenetic data at a second timepoint (T2) identified 49 significant CpGs. EWAS Atlas enrichment analysis indicated that 4 of the 49 were previously associated with asthma or asthma-related traits. Comparing the estimates of all the significant associations identified across the two time points (271 in total) yielded a correlation of 0.81. Conclusion: We identified 270 unique CpGs that were associated with PC scores generated from 35 clinical markers of asthma, either cross-sectionally or 10 years later. A strong correlation was present between effect sizes at the 2 timepoints. Most associations were identified for PC7, which captured blood eosinophil counts and immunoglobulin levels and many of these CpGs have previous associations in earlier studies of asthma and asthma-related traits. The results point to using this robust DNA methylation profile as a new, stable biomarker for asthma.

A Bivariate Twin Study of Lifetime cannabis Initiation and Lifetime Regular Tobacco Smoking Across Three Different Countries.

Regular cigarette smoking and cannabis consumption are strongly positively related to each other, yet few studies explore their underlying variation and covariation. We evaluated the genetic and environmental decomposition of variance and covariance of these two traits in twin data from three countries with different social norms and legislation. Data from the Netherlands Twin Register, FinnTwin12/16, and the Minnesota Center for Twin Family Research (total N = 21,617) were analyzed in bivariate threshold models of lifetime regular smoking initiation (RSI) and lifetime cannabis initiation (CI). We ran unstratified models and models stratified by sex and country. Prevalence of RSI was lowest in the Netherlands and prevalence of CI was highest in Minnesota. In the unstratified model, genetic (A) and common environmental factors (C) contributed substantially to the liabilities of RSI (A = 0.47, C = 0.34) and CI (A = 0.28, C = 0.51). The two liabilities were significantly phenotypically (rP = 0.56), genetically (rA = 0.74), and environmentally correlated in the unstratified model (rC = 0.47and rE = 0.48, representing correlations between common and unique environmental factors). The magnitude of phenotypic correlation between liabilities varied by country but not sex (Minnesota rP ~ 0.70, Netherlands rP ~ 0.59, Finland rP ~ 0.45). Comparisons of decomposed correlations could not be reliably tested in the stratified models. The prevalence and association of RSI and CI vary by sex and country. These two behaviors are correlated because there is genetic and environmental overlap between their underlying latent liabilities. There is heterogeneity in the genetic architecture of these traits across country.

Brain-age prediction: Systematic evaluation of site effects, and sample age range and size.

Structural neuroimaging data have been used to compute an estimate of the biological age of the brain (brain-age) which has been associated with other biologically and behaviorally meaningful measures of brain development and aging. The ongoing research interest in brain-age has highlighted the need for robust and publicly available brain-age models pre-trained on data from large samples of healthy individuals. To address this need we have previously released a developmental brain-age model. Here we expand this work to develop, empirically validate, and disseminate a pre-trained brain-age model to cover most of the human lifespan. To achieve this, we selected the best-performing model after systematically examining the impact of seven site harmonization strategies, age range, and sample size on brain-age prediction in a discovery sample of brain morphometric measures from 35,683 healthy individuals (age range: 5-90 years; 53.59% female). The pre-trained models were tested for cross-dataset generalizability in an independent sample comprising 2101 healthy individuals (age range: 8-80 years; 55.35% female) and for longitudinal consistency in a further sample comprising 377 healthy individuals (age range: 9-25 years; 49.87% female). This empirical examination yielded the following findings: (1) the accuracy of age prediction from morphometry data was higher when no site harmonization was applied; (2) dividing the discovery sample into two age-bins (5-40 and 40-90 years) provided a better balance between model accuracy and explained age variance than other alternatives; (3) model accuracy for brain-age prediction plateaued at a sample size exceeding 1600 participants. These findings have been incorporated into CentileBrain (https://centilebrain.org/#/brainAGE2), an open-science, web-based platform for individualized neuroimaging metrics.

Overlapping genetic influences between creativity and borderline personality symptoms in a large Dutch sample.

Creativity and mental disorders are sometimes seen as intertwined, but research is still unclear on whether, how much, and why. Here we explore the potential role of shared genetic factors behind creativity and symptoms of borderline personality disorder (BPD, characterized by mood swings and randomness of thoughts). Data were collected from 6745 twins (2378 complete pairs) by the Netherlands Twin Register on BPD scores (PAI-BOR questionnaire) and working in a creative profession (proxy for creativity). First, we tested whether there is an association between BPD symptoms and creative professions. Results confirmed that individuals scoring higher on the BPD spectrum are more likely to have a creative profession (Cohen's d = 0.16). Next, we modeled how much of this association reflects underlying genetic and/or environmental correlations-by using a bivariate classical twin design. We found that creativity and BPD were each influenced by genetic factors (heritability = 0.45 for BPD and 0.67 for creativity) and that these traits are genetically correlated rG = 0.17. Environmental influences were not correlated. This is evidence for a common genetic mechanism between borderline personality scores and creativity which may reflect causal effects and shed light on mechanisms.

Twin-based Mendelian Randomization Analyses Highlight Smoking's Effects on Blood DNA Methylation, with Putative Reverse Causation.

Epigenome-wide association studies (EWAS) aim to identify differentially methylated loci associated with complex traits and disorders. EWAS of cigarette smoking shows some of the most widespread DNA methylation (DNAm) associations in blood. However, traditional EWAS cannot differentiate between causation and confounding, leading to ambiguity in etiological interpretations. Here, we apply an integrated approach combining Mendelian Randomization and twin-based Direction-of-Causation analyses (MR-DoC) to examine causality underlying smoking-associated blood DNAm changes in the Netherlands Twin Register (N=2577). Evidence across models suggests that current smoking's causal effects on DNAm likely drive many of the previous EWAS findings, implicating functional pathways relevant to several adverse health outcomes of smoking, including hemopoiesis, cell- and neuro-development, and immune regulation. Additionally, we find evidence of potential reverse causal influences at some DNAm sites, with 17 of these sites enriched for gene regulatory functional elements in the brain. The top three sites with evidence of DNAm's effects on smoking annotate to genes involved in G protein-coupled receptor signaling (GNG7, RGS3) and innate immune response (SLC15A4), elucidating potential biological risk factors for smoking. This study highlights the utility of integrating genotypic and DNAm measures in twin cohorts to clarify the causal relationships between health behaviors and blood DNAm.

Exposome-Wide Association Study of Body Mass Index Using a Novel Meta-Analytical Approach for Random Forest Models.

BACKGROUND: Overweight and obesity impose a considerable individual and social burden, and the urban environments might encompass factors that contribute to obesity. Nevertheless, there is a scarcity of research that takes into account the simultaneous interaction of multiple environmental factors. OBJECTIVES: Our objective was to perform an exposome-wide association study of body mass index (BMI) in a multicohort setting of 15 studies. METHODS: Studies were affiliated with the Dutch Geoscience and Health Cohort Consortium (GECCO), had different population sizes (688-141,825), and covered the entire Netherlands. Ten studies contained general population samples, others focused on specific populations including people with diabetes or impaired hearing. BMI was calculated from self-reported or measured height and weight. Associations with 69 residential neighborhood environmental factors (air pollution, noise, temperature, neighborhood socioeconomic and demographic factors, food environment, drivability, and walkability) were explored. Random forest (RF) regression addressed potential nonlinear and nonadditive associations. In the absence of formal methods for multimodel inference for RF, a rank aggregation-based meta-analytic strategy was used to summarize the results across the studies. RESULTS: Six exposures were associated with BMI: five indicating neighborhood economic or social environments (average home values, percentage of high-income residents, average income, livability score, share of single residents) and one indicating the physical activity environment (walkability in 5-km buffer area). Living in high-income neighborhoods and neighborhoods with higher livability scores was associated with lower BMI. Nonlinear associations were observed with neighborhood home values in all studies. Lower neighborhood home values were associated with higher BMI scores but only for values up to €300,000. The directions of associations were less consistent for walkability and share of single residents. DISCUSSION: Rank aggregation made it possible to flexibly combine the results from various studies, although between-study heterogeneity could not be estimated quantitatively based on RF models. Neighborhood social, economic, and physical environments had the strongest associations with BMI. https://doi.org/10.1289/EHP13393.

Genome-wide analyses identify 21 infertility loci and over 400 reproductive hormone loci across the allele frequency spectrum.

Genome-wide association studies (GWASs) may help inform treatments for infertility, whose causes remain unknown in many cases. Here we present GWAS meta-analyses across six cohorts for male and female infertility in up to 41,200 cases and 687,005 controls. We identified 21 genetic risk loci for infertility (P≤5E-08), of which 12 have not been reported for any reproductive condition. We found positive genetic correlations between endometriosis and all-cause female infertility (rg=0.585, P=8.98E-14), and between polycystic ovary syndrome and anovulatory infertility (rg=0.403, P=2.16E-03). The evolutionary persistence of female infertility-risk alleles in EBAG9 may be explained by recent directional selection. We additionally identified up to 269 genetic loci associated with follicle-stimulating hormone (FSH), luteinising hormone, oestradiol, and testosterone through sex-specific GWAS meta-analyses (N=6,095-246,862). While hormone-associated variants near FSHB and ARL14EP colocalised with signals for anovulatory infertility, we found no rg between female infertility and reproductive hormones (P>0.05). Exome sequencing analyses in the UK Biobank (N=197,340) revealed that women carrying testosterone-lowering rare variants in GPC2 were at higher risk of infertility (OR=2.63, P=1.25E-03). Taken together, our results suggest that while individual genes associated with hormone regulation may be relevant for fertility, there is limited genetic evidence for correlation between reproductive hormones and infertility at the population level. We provide the first comprehensive view of the genetic architecture of infertility across multiple diagnostic criteria in men and women, and characterise its relationship to other health conditions.

Enhanced resolution profiling in twins reveals differential methylation signatures of type 2 diabetes with links to its complications.

BACKGROUND: Type 2 diabetes (T2D) susceptibility is influenced by genetic and environmental factors. Previous findings suggest DNA methylation as a potential mechanism in T2D pathogenesis and progression. METHODS: We profiled DNA methylation in 248 blood samples from participants of European ancestry from 7 twin cohorts using a methylation sequencing platform targeting regulatory genomic regions encompassing 2,048,698 CpG sites. FINDINGS: We find and replicate 3 previously unreported T2D differentially methylated CpG positions (T2D-DMPs) at FDR 5% in RGL3, NGB and OTX2, and 20 signals at FDR 25%, of which 14 replicated. Integrating genetic variation and T2D-discordant monozygotic twin analyses, we identify both genetic-based and genetic-independent T2D-DMPs. The signals annotate to genes with established GWAS and EWAS links to T2D and its complications, including blood pressure (RGL3) and eye disease (OTX2). INTERPRETATION: The results help to improve our understanding of T2D disease pathogenesis and progression and may provide biomarkers for its complications. FUNDING: Funding acknowledgements for each cohort can be found in the Supplementary Note.

Genetic Insights into Externalizing and Internalizing Traits through Integration of the Research Domain Criteria and Hierarchical Taxonomy of Psychopathology.

Background: There is considerable comorbidity between externalizing (EXT) and internalizing (INT) psychopathology. Understanding the shared genetic underpinnings of these spectra is crucial for advancing knowledge of their biological bases and potential health impacts, and for informing empirical models like the Research Domain Criteria (RDoC) and Hierarchical Taxonomy of Psychopathology (HiTOP). Methods: We conducted a multivariate genome-wide association study (GWAS) of EXT and INT psychopathology by applying genomic structural equation modeling to summary statistics from 16 EXT and INT traits in European-ancestry individuals (n = 16,400 to 1,074,629). Downstream analyses explored associations across RDoC units of analysis (i.e., genes, molecules, cells, circuits, physiology, and behaviors). Results: The GWAS identified 409 lead single nucleotide polymorphisms (SNPs) for EXT, 85 for INT, and 256 for EXT+INT (i.e., shared) traits. Bivariate causal mixture models estimated that nearly all EXT and INT causal variants overlapped, despite a genetic correlation of 0.37 (SE = 0.02). Drug repurposing analyses identified potential therapeutic targets, including perturbagens affecting dopamine and serotonin pathways. EXT genes had enriched expression in GABAergic, cortical, and hippocampal neurons, while INT genes were more narrowly linked to GABAergic neurons. EXT+INT liability was associated with reduced grey matter volumes in the amygdala and subcallosal cortex. Conclusions: These findings reveal both genetic overlap and distinct molecular and neurobiological pathways underlying EXT and INT psychopathology. By integrating genomic insights with the RDoC and HiTOP frameworks, this study advances our understanding of the mechanisms driving these dimensions of psychopathology.

Breast Cancer Polygenic Risk Score Validation and Effects of Variable Imputation.

Breast cancer (BC) is a complex disease affecting one in eight women in the USA. Advances in population genomics have led to the development of polygenic risk scores (PRSs) with the potential to augment current risk models, but replication is often limited. We evaluated 2 robust PRSs with 313 and 3820 SNPs and the effects of multiple genotype imputation replications in BC cases and control populations. Biological samples from BC cases and cancer-free controls were drawn from three European ancestry cohorts. Genotyping on the Illumina Global Screening Array was followed by stringent quality control measures and 20 genotype imputation replications. A total of 468 unrelated cases and 4337 controls were scored, revealing significant differences in mean PRS percentiles between cases and controls (p < 0.001) for both SNP sets (313-SNP PRS: 52.81 and 48.07; 3820-SNP PRS: 55.45 and 49.81), with receiver operating characteristic curve analysis showing area under the curve values of 0.596 and 0.603 for the 313-SNP and 3820-SNP PRS, respectively. PRS fluctuations (from ~2-3% up to 9%) emerged across imputation iterations. Our study robustly reaffirms the predictive capacity of PRSs for BC by replicating their performance in an independent BC population and showcases the need to average imputed scores for reliable outcomes.

Genomics Research of Lifetime Depression in the Netherlands: The BIObanks Netherlands Internet Collaboration (BIONIC) Project.

In this cohort profile article we describe the lifetime major depressive disorder (MDD) database that has been established as part of the BIObanks Netherlands Internet Collaboration (BIONIC). Across the Netherlands we collected data on Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) lifetime MDD diagnosis in 132,850 Dutch individuals. Currently, N = 66,684 of these also have genomewide single nucleotide polymorphism (SNP) data. We initiated this project because the complex genetic basis of MDD requires large population-wide studies with uniform in-depth phenotyping. For standardized phenotyping we developed the LIDAS (LIfetime Depression Assessment Survey), which then was used to measure MDD in 11 Dutch cohorts. Data from these cohorts were combined with diagnostic interview depression data from 5 clinical cohorts to create a dataset of N = 29,650 lifetime MDD cases (22%) meeting DSM-5 criteria and 94,300 screened controls. In addition, genomewide genotype data from the cohorts were assembled into a genomewide association study (GWAS) dataset of N = 66,684 Dutch individuals (25.3% cases). Phenotype data include DSM-5-based MDD diagnoses, sociodemographic variables, information on lifestyle and BMI, characteristics of depressive symptoms and episodes, and psychiatric diagnosis and treatment history. We describe the establishment and harmonization of the BIONIC phenotype and GWAS datasets and provide an overview of the available information and sample characteristics. Our next step is the GWAS of lifetime MDD in the Netherlands, with future plans including fine-grained genetic analyses of depression characteristics, international collaborations and multi-omics studies.

Genome-Wide Association Study of Obsessive-Compulsive Symptoms including 33,943 individuals from the general population.

While 1-2% of individuals meet the criteria for a clinical diagnosis of obsessive-compulsive disorder (OCD), many more (~13-38%) experience subclinical obsessive-compulsive symptoms (OCS) during their life. To characterize the genetic underpinnings of OCS and its genetic relationship to OCD, we conducted the largest genome-wide association study (GWAS) meta-analysis of parent- or self-reported OCS to date (N = 33,943 with complete phenotypic and genome-wide data), combining the results from seven large-scale population-based cohorts from Sweden, the Netherlands, England, and Canada (including six twin cohorts and one cohort of unrelated individuals). We found no genome-wide significant associations at the single-nucleotide polymorphism (SNP) or gene-level, but a polygenic risk score (PRS) based on the OCD GWAS previously published by the Psychiatric Genetics Consortium (PGC-OCD) was significantly associated with OCS (Pfixed = 3.06 × 10-5). Also, one curated gene set (Mootha Gluconeogenesis) reached Bonferroni-corrected significance (Ngenes = 28, Beta = 0.79, SE = 0.16, Pbon = 0.008). Expression of genes in this set is high at sites of insulin mediated glucose disposal. Dysregulated insulin signaling in the etiology of OCS has been suggested by a previous study describing a genetic overlap of OCS with insulin signaling-related traits in children and adolescents. We report a SNP heritability of 4.1% (P = 0.0044) in the meta-analyzed GWAS, and heritability estimates based on the twin cohorts of 33-43%. Genetic correlation analysis showed that OCS were most strongly associated with OCD (rG = 0.72, p = 0.0007) among all tested psychiatric disorders (N = 11). Of all 97 tested phenotypes, 24 showed a significant genetic correlation with OCS, and 66 traits showed concordant directions of effect with OCS and OCD. OCS have a significant polygenic contribution and share genetic risk with diagnosed OCD, supporting the hypothesis that OCD represents the extreme end of widely distributed OCS in the population.

Normative modelling of brain morphometry across the lifespan with CentileBrain: algorithm benchmarking and model optimisation.

The value of normative models in research and clinical practice relies on their robustness and a systematic comparison of different modelling algorithms and parameters; however, this has not been done to date. We aimed to identify the optimal approach for normative modelling of brain morphometric data through systematic empirical benchmarking, by quantifying the accuracy of different algorithms and identifying parameters that optimised model performance. We developed this framework with regional morphometric data from 37 407 healthy individuals (53% female and 47% male; aged 3-90 years) from 87 datasets from Europe, Australia, the USA, South Africa, and east Asia following a comparative evaluation of eight algorithms and multiple covariate combinations pertaining to image acquisition and quality, parcellation software versions, global neuroimaging measures, and longitudinal stability. The multivariate fractional polynomial regression (MFPR) emerged as the preferred algorithm, optimised with non-linear polynomials for age and linear effects of global measures as covariates. The MFPR models showed excellent accuracy across the lifespan and within distinct age-bins and longitudinal stability over a 2-year period. The performance of all MFPR models plateaued at sample sizes exceeding 3000 study participants. This model can inform about the biological and behavioural implications of deviations from typical age-related neuroanatomical changes and support future study designs. The model and scripts described here are freely available through CentileBrain.