Netherlands Twin Register: From Twins to Twin Families

Metabolomics profiling of Type D personality traits

OBJECTIVE

Type D (Distressed) personality combines negative affectivity (NA) and social inhibition (SI) and is associated with an increased risk of cardiometabolic diseases. Here, we examined the association of Type D traits with 230 (predominantly) lipid metabolites and metabolite ratios.

METHODS

Four Dutch cohorts were included, comprising 10,834 individuals. Type D personality traits were measured by self-report questionnaires. A proton nuclear magnetic resonance (NMR) metabolomics platform provided 149 absolute measures (98 belonging to lipoprotein subclasses) and 81 derived ratios. For all, linear regression analyses were performed within each cohort, followed by random-effects meta-analyses. A per-measure FDR q-value<0.05 was set as a study-wise significant association.

RESULTS

SI was significantly associated with a lower omega-3 fatty acids to total fatty acids (FAw3.FA%) ratio, and a lower free cholesterol to total lipids ratio in very small VLDL (XS.VLDL.FC%). FAw3.FA% was also associated to NA (no study-wise significance though). NA showed a suggestive replication (p-value<.05) of the previous reported associations with depression for 5 out of 18 metabolites from the same metabolomics platform: triglycerides in HDL, serum total triglycerides, VLDL cholesterol, mean diameter for VLDL particles and VLDL triglycerides.

CONCLUSIONS

In this large meta-analysis, SI was associated with omega-3 fatty acids to total fatty acids ratio, which is suggestive of lower omega-3 fatty acid intake. Only some metabolite biomarkers showed tentative links to Type D and NA. In sum, it seems that there are no major alterations in lipid metabolism associated with Type D traits.

Longitudinal Assessment of Retinal Microvasculature in Preclinical Alzheimer's Disease

Purpose

To investigate if changes in vessel density (VD) and the foveal avascular zone (FAZ) occur in the preclinical phase of Alzheimer's disease (pAD) over time.

Methods

Optical coherence tomography angiography (OCTA) was used to image VD and FAZ at baseline and for a follow-up period of 2 years. Positron emission tomography (PET) was used to determine the amyloid beta (Aβ) status of participants.

Results

The VD and FAZ of 148 participants (54% female) were analyzed at baseline and follow-up (mean time between measurements, 2.24 ± 0.35 years). The mean age of the participants was 68.3 ± 6.0 years at baseline and 70.3 ± 5.9 years at follow-up. Participants were divided into three groups: control group, participants who had negative Aβ status at both measurements (Aβ-, n = 116); converter group, participants who transitioned from negative to positive between baseline and follow-up (Aβ-+, n = 18); and participants who were consistently positive at both visits (Aβ++, n = 14). The VD of both Aβ+ groups demonstrated non-significant increases over time in both macula and optic nerve head (ONH) regions. The Aβ- group was found to be significantly higher in both ONH and macular regions. The VD of the Aβ++ group was significantly higher in the macula inner and outer rings compared to the Aβ-+ and Aβ- groups. No significant change was found in FAZ values over time.

Conclusions

Alterations in VD seem to manifest already in pAD, exhibiting distinct variations between the ONH and macula. Further longitudinal studies with a longer follow-up design and known amyloid pathology should be undertaken to validate these observations.

Post-GWAS multiomic functional investigation of the TNIP1 locus in Alzheimer's disease highlights a potential role for GPX3

INTRODUCTION

Recent genome-wide association studies (GWAS) have reported a genetic association with Alzheimer's disease (AD) at the TNIP1/GPX3 locus, but the mechanism is unclear.

METHODS

We used cerebrospinal fluid (CSF) proteomics data to test (n = 137) and replicate (n = 446) the association of glutathione peroxidase 3 (GPX3) with CSF biomarkers (including amyloid and tau) and the GWAS-implicated variants (rs34294852 and rs871269).

RESULTS

CSF GPX3 levels decreased with amyloid and tau positivity (analysis of variance P = 1.5 × 10-5) and higher CSF phosphorylated tau (p-tau) levels (P = 9.28 × 10-7). The rs34294852 minor allele was associated with decreased GPX3 (P = 0.041). The replication cohort found associations of GPX3 with amyloid and tau positivity (P = 2.56 × 10-6) and CSF p-tau levels (P = 4.38 × 10-9).

DISCUSSION

These results suggest variants in the TNIP1 locus may affect the oxidative stress response in AD via altered GPX3 levels.

HIGHLIGHTS

Cerebrospinal fluid (CSF) glutathione peroxidase 3 (GPX3) levels decreased with amyloid and tau positivity and higher CSF phosphorylated tau. The minor allele of rs34294852 was associated with lower CSF GPX3. levels when also controlling for amyloid and tau category. GPX3 transcript levels in the prefrontal cortex were lower in Alzheimer's disease than controls. rs34294852 is an expression quantitative trait locus for GPX3 in blood, neutrophils, and microglia.

Heritability of adult picky eating in the Netherlands

Adult picky eating (APE), the rejection of familiar and unfamiliar foods leading to a diet with limited variety, is an understudied phenomenon which can have both physical and psychological negative consequences. The aetiology of individual differences in APE is understudied, although there is reason to believe that it is partly heritable. Therefore, we aimed to estimate the heritability of APE with data from the Netherlands Twin Register (n = 8016) with classical genetic structural equation modelling. In order to use these data, we firstly investigated whether a Food Preference Questionnaire (FPQ) could measure APE with a pre-registered prestudy. Adult participants (n = 414) filled in online questionnaires, including a FPQ and measures related to APE. Spearman's rho correlation quantified the relationship between different elements of the Dutch FPQ and different scores on measures of APE. Results of the prestudy showed that the mean liking score on the FPQ could be used to measure APE (ρ > .50). This measure was then used in the main study to estimate the heritability of APE. Results showed that broad-sense heritability for APE is 49 % (additive genetic effects 14 % (95 % CI [00, 38]) + dominance genetic effects 35 % (95 % CI [11, 52]), while the remaining variance is explained by unique environmental factors. Future studies may focus on uncovering the specific genetic and unique environmental factors that play a role in APE.

Contributions of amyloid beta and cerebral small vessel disease in clinical decline

INTRODUCTION

We assessed whether co-morbid small vessel disease (SVD) has clinical predictive value in preclinical or prodromal Alzheimer's disease.

METHODS

In 1090 non-demented participants (65.4 ± 10.7 years) SVD was assessed with magnetic resonance imaging and amyloid beta (Aβ) with lumbar puncture and/or positron emission tomography scan (mean follow-up for cognitive function 3.1 ± 2.4 years).

RESULTS

Thirty-nine percent had neither Aβ nor SVD (A-V-), 21% had SVD only (A-V+), 23% Aβ only (A+V-), and 17% had both (A+V+). Pooled cohort linear mixed model analyses demonstrated that compared to A-V- (reference), A+V- had a faster rate of cognitive decline. Co-morbid SVD (A+V+) did not further increase rate of decline. Cox regression showed that dementia risk was modestly increased in A-V+ (hazard ratio [95% confidence interval: 1.8 [1.0-3.2]) and most strongly in A+ groups. Also, mortality risk was increased in A+ groups.

DISCUSSION

In non-demented persons Aβ was predictive of cognitive decline, dementia, and mortality. SVD modestly predicts dementia in A-, but did not increase deleterious effects in A+.

HIGHLIGHTS

Amyloid beta (Aβ; A) was predictive for cognitive decline, dementia, and mortality. Small vessel disease (SVD) had no additional deleterious effects in A+. SVD modestly predicted dementia in A-. Aβ should be assessed even when magnetic resonance imaging indicates vascular cognitive impairment.

Trans-ancestry epigenome-wide association meta-analysis of DNA methylation with lifetime cannabis use

Cannabis is widely used worldwide, yet its links to health outcomes are not fully understood. DNA methylation can serve as a mediator to link environmental exposures to health outcomes. We conducted an epigenome-wide association study (EWAS) of peripheral blood-based DNA methylation and lifetime cannabis use (ever vs. never) in a meta-analysis including 9436 participants (7795 European and 1641 African ancestry) from seven cohorts. Accounting for effects of cigarette smoking, our trans-ancestry EWAS meta-analysis revealed four CpG sites significantly associated with lifetime cannabis use at a false discovery rate of 0.05 ( p < 5.85 × 10 - 7 ) : cg22572071 near gene ADGRF1, cg15280358 in ADAM12, cg00813162 in ACTN1, and cg01101459 near LINC01132. Additionally, our EWAS analysis in participants who never smoked cigarettes identified another epigenome-wide significant CpG site, cg14237301 annotated to APOBR. We used a leave-one-out approach to evaluate methylation scores constructed as a weighted sum of the significant CpGs. The best model can explain 3.79% of the variance in lifetime cannabis use. These findings unravel the DNA methylation changes associated with lifetime cannabis use that are independent of cigarette smoking and may serve as a starting point for further research on the mechanisms through which cannabis exposure impacts health outcomes.

The relationship of school performance with self-control and grit is strongly genetic and weakly causal

The non-cognitive skills self-control and grit are often considered predictors of school performance, but whether this relationship is causal remains unclear. We investigated the causality of this association using a twin design. Specifically, we evaluated the direct impact of self-control and grit on school performance, while controlling for genetic or environmental influences common to all three traits (i.e., confounding). Teachers of 4891 Dutch 12-year-old twin pairs (of which 3837 were complete pairs) completed a survey about school performance (school grades), self-control (ASEBA self-control scale), and the perseverance aspect of grit. Our analysis aimed to determine the direct impact of self-control and grit on school performance, while simultaneously controlling for genetic or environmental confounding. Establishing the regression relationship corrected for confounding supports the interpretation of the regression relationship as causal. In all analyses, we corrected for sex, rater bias of the teachers, and parental socioeconomic status. Initially, in the standard regression, self-control, and grit explained 28.4% of the school performance variance. However, allowing for genetic confounding (due to genetic pleiotropy) revealed that most of this association could be attributed to genetic influences that the three traits share. In the presence of genetic pleiotropy, the phenotypic regression of school performance on self-control and grit accounted for only 4.4% (i.e., the effect size association with the causal hypothesis). In conclusion, self-control and grit predict school performance primarily due to genetic pleiotropy, with a much smaller causal effect (R2 = 4.4%). This suggests that interventions targeting self-control and grit alone may yield limited improvements in school performance.

Quantitative trait loci mapping of circulating metabolites in cerebrospinal fluid to uncover biological mechanisms involved in brain-related phenotypes

Genomic studies of molecular traits have provided mechanistic insights into complex disease, though these lag behind for brain-related traits due to the inaccessibility of brain tissue. We leveraged cerebrospinal fluid (CSF) to study neurobiological mechanisms in vivo , measuring 5,543 CSF metabolites, the largest panel in CSF to date, in 977 individuals of European ancestry. Individuals originated from two separate cohorts including cognitively healthy subjects (n=490) and a well-characterized memory clinic sample, the Amsterdam Dementia Cohort (ADC, n=487). We performed metabolite quantitative trait loci (mQTL) mapping on CSF metabolomics and found 126 significant mQTLs, representing 65 unique CSF metabolites across 51 independent loci. To better understand the role of CSF mQTLs in brain-related disorders, we performed a metabolome-wide association study (MWAS), identifying 40 associations between CSF metabolites and brain traits. Similarly, over 90% of significant mQTLs demonstrated colocalized associations with brain-specific gene expression, unveiling potential neurobiological pathways.

Self-control and grit are associated with school performance mainly because of shared genetic effects

Background

By combining the classical twin design with regression analysis, we investigated the role of two non-cognitive factors, self-control and grit, in the prediction of school performance. We did so at the phenotypic, genetic, and environmental level.

Methods

Teachers filled out a survey on the twins' school performance (school grades for reading, literacy, and math), self-control (ASEBA self-control scale), and grit (the perseverance aspect) for 4891 Dutch 12-years-old twin pairs (3837 pairs with data for both and 1054 pairs with data for one of the twins). We employed regression analyses to first assess the contributions of self-control and grit to school performance at the phenotypic level, and next at the genetic and environmental level, while correcting for rater (teacher) effects, parental SES, and sex.

Results

Higher SES was associated with better school performance, self-control, and grit. On average, girls had more self-control and grit than boys. Corrected for sex, SES, and teacher rater effects, genetic factors accounted for 74%, 69%, and 58% of the phenotypic variance of school performance, self-control, and grit, respectively. Phenotypically, self-control and grit explained 28.3% of the variance in school performance. We found that this phenotypic result largely reflected genetic influences.

Conclusions

Children who have better self-control and are grittier tend to do better in school. Individual differences in these three traits are not correlated because of shared environmental influences, but mainly because of shared genetic factors.

Novel insight into the aetiology of rheumatoid arthritis gained by a cross-tissue transcriptome-wide association study

BACKGROUND

Although genome-wide association studies (GWASs) have identified more than 100 loci associated with rheumatoid arthritis (RA) susceptibility, the causal genes and biological mechanisms remain largely unknown.

METHODS

A cross-tissue transcriptome-wide association study (TWAS) using the unified test for molecular signaturestool was performed to integrate GWAS summary statistics from 58 284 individuals (14 361 RA cases and 43 923 controls) with gene-expression matrix in the Genotype-Tissue Expression project. Subsequently, a single tissue by using FUSION software was conducted to validate the significant associations. We also compared the TWAS with different gene-based methodologies, including Summary Data Based Mendelian Randomization (SMR) and Multimarker Analysis of Genomic Annotation (MAGMA). Further in silico analyses (conditional and joint analysis, differential expression analysis and gene-set enrichment analysis) were used to deepen our understanding of genetic architecture and comorbidity aetiology of RA.

RESULTS

We identified a total of 47 significant candidate genes for RA in both cross-tissue and single-tissue test after multiple testing correction, of which 40 TWAS-identified genes were verified by SMR or MAGMA. Among them, 13 genes were situated outside of previously reported significant loci by RA GWAS. Both TWAS-based and MAGMA-based enrichment analyses illustrated the shared genetic determinants among autoimmune thyroid disease, asthma, type I diabetes mellitus and RA.

CONCLUSION

Our study unveils 13 new candidate genes whose predicted expression is associated with risk of RA, providing new insights into the underlying genetic architecture of RA.

Genome- and transcriptome-wide association studies show that pulmonary embolism is associated with bone-forming proteins

BACKGROUND

Pulmonary embolism (PE) is a leading cause of death in stroke patients and a severe health burden worldwide. There is a pressing need to understand the mechanisms by which it occurs and to identify at-risk patients efficiently and accurately.

OBJECTIVES

The aim of this paper was to analyze the genetic correlation between PE and human plasma proteins through genome-wide association study (GWAS) with transcriptome-wide association study (TWAS), in combination with mRNA expression profiling at three levels: DNA, RNA, and protein.

METHODS

First, based on data from GWAS in European populations, we performed a linkage disequilibrium score regression (LDSC) analysis of plasma proteins and PE in 3,283 individuals and additionally analyzed the genetic association between PE and fracture. Then, we performed a TWAS on PE GWAS data using skeletal muscle and blood for gene expression references. Finally, we validated the genetic correlation between PE and human plasma proteins by co-matching the genes encoding the identified proteins and those identified using TWAS with the differentially expressed genes obtained from mRNA expression profiling of PE (Figure1).

RESULTS

We identified five plasma proteins associated with PE, including hydroxycarboxylic acid receptor 2, defensin 118, and bone morphogenetic protein (BMP) 7, as well as a relationship between PE and fracture. Comparison of genes encoding these proteins with genes obtained from TWAS and then with differentially expressed genes obtained from PE mRNA expression profiling revealed that PE was highly correlated with the BMP family of genes.

Multivariate Genetic Structure of Externalizing Behavior and Structural Brain Development in a Longitudinal Adolescent Twin Sample

Externalizing behavior in its more extreme form is often considered a problem to the individual, their families, teachers, and society as a whole. Several brain structures have been linked to externalizing behavior and such associations may arise if the (co)development of externalizing behavior and brain structures share the same genetic and/or environmental factor(s). We assessed externalizing behavior with the Child Behavior Checklist and Youth Self Report, and the brain volumes and white matter integrity (fractional anisotropy [FA] and mean diffusivity [MD]) with magnetic resonance imaging in the BrainSCALE cohort, which consisted of twins and their older siblings from 112 families measured longitudinally at ages 10, 13, and 18 years for the twins. Genetic covariance modeling based on the classical twin design, extended to also include siblings of twins, showed that genes influence externalizing behavior and changes therein (h² up to 88%). More pronounced externalizing behavior was associated with higher FA (observed correlation r_ph up to +0.20) and lower MD (r_ph up to −0.20), with sizeable genetic correlations (FA r_a up to +0.42; MD r_a up to −0.33). The cortical gray matter (CGM; r_ph up to −0.20) and cerebral white matter (CWM; r_ph up to +0.20) volume were phenotypically but not genetically associated with externalizing behavior. These results suggest a potential mediating role for global brain structures in the display of externalizing behavior during adolescence that are both partially explained by the influence of the same genetic factor.

Subgroup Identification and Regression Analysis of Clustered and Heterogeneous Interval-Censored Data

Clustered and heterogeneous interval-censored data occur in many fields such as medical studies. For example, in a migraine study with the Netherlands Twin Registry, the information including time to diagnosis of migraine and gender was collected for 3975 monozygotic and dizygotic twins. Since each study subject is observed only at discrete and periodic follow-up time points, the failure times of interest (i.e., the time when the individual first had a migraine) are known only to belong to certain intervals and hence are interval-censored. Furthermore, these twins come from different genetic backgrounds and may be associated with differential risks for developing migraines. For simultaneous subgroup identification and regression analysis of such data, we propose a latent Cox model where the number of subgroups is not assumed a priori but rather data-driven estimated. The nonparametric maximum likelihood method and an EM algorithm with monotone ascent property are also developed for estimating the model parameters. Simulation studies are conducted to assess the finite sample performance of the proposed estimation procedure. We further illustrate the proposed methodologies by an empirical analysis of migraine data.

Genetic association of primary nonresponse to anti-TNFα therapy in patients with inflammatory bowel disease

OBJECTIVES

Primary nonresponse (PNR) to antitumor necrosis factor-α (TNFα) biologics is a serious concern in patients with inflammatory bowel disease (IBD). We aimed to identify the genetic variants associated with PNR.

PATIENTS AND METHODS

Patients were recruited from outpatient GI clinics and PNR was determined using both clinical and endoscopic findings. A case-control genome-wide association study was performed in 589 IBD patients and associations were replicated in an independent cohort of 293 patients. Effect of the associated variant on gene expression and TNFα secretion was assessed by cell-based assays. Pleiotropic effects were investigated by Phenome-wide association study (PheWAS).

RESULTS

We identified rs34767465 as associated with PNR to anti-TNFα therapy (odds ratio: 2.07, 95% CI, 1.46-2.94, P = 2.43 × 10-7, [replication odds ratio: 1.8, 95% CI, 1.04-3.16, P = 0.03]). rs34767465 is a multiple-tissue expression quantitative trait loci for FAM114A2. Using RNA-sequencing and protein quantification from HapMap lymphoblastoid cell lines (LCLs), we found a significant decrease in FAM114A2 mRNA and protein expression in both heterozygous and homozygous genotypes when compared to wild type LCLs. TNFα secretion was significantly higher in THP-1 cells [differentiated into macrophages] with FAM114A2 knockdown versus controls. Immunoblotting experiments showed that depletion of FAM114A2 impaired autophagy-related pathway genes suggesting autophagy-mediated TNFα secretion as a potential mechanism. PheWAS showed rs34767465 was associated with comorbid conditions found in IBD patients (derangement of joints [P = 3.7 × 10-4], pigmentary iris degeneration [P = 5.9 × 10-4], diverticulum of esophagus [P = 7 × 10-4]).

CONCLUSIONS

We identified a variant rs34767465 associated with PNR to anti-TNFα biologics, which increases TNFα secretion through mechanism related to autophagy. rs34767465 may also explain the comorbidities associated with IBD.

Transcriptome-wide association study identifies multiple genes and pathways associated with thyroid function

Thyroid dysfunction is a common endocrine disease measured by thyroid-stimulating hormone (TSH) level. Although more than 70 genetic loci associated with TSH have been reported through genome-wide association studies (GWASs), the variants can only explain a small fraction of the thyroid function heritability. To identify novel candidate genes for thyroid function, we conducted the first large-scale transcriptome-wide association study (TWAS) for thyroid function using GWAS-summary data for TSH levels in up to 119 715 individuals combined with pre-computed gene expression weights of six panels from four tissue types. The candidate genes identified by TWAS were further validated by TWAS replication and gene expression profiles. We identified 74 conditionally independent genes significantly associated with thyroid function, such as PDE8B (P = 1.67 × 10-282), PDE10A (P = 7.61 × 10-119), NR3C2 (P = 1.50 × 10-92), and CAPZB (P = 3.13 × 10-79). After TWAS replication using UKBB datasets, 26 genes were replicated for significant associations with thyroid-relevant diseases/traits. Among them, 16 gene were causal for their associations to thyroid-relevant diseases/traits and further validated in differential expression analyses, including two novel genes (MFSD6 and RBM47) that did not implicate in previous GWASs. Enrichment analyses detected several pathways associated with thyroid function, such as the cAMP signaling pathway (P = 7.27 × 10-4), hemostasis (P = 3.74 × 10-4), and platelet activation, signaling, and aggregation (P = 9.98 × 10-4). Our study identified multiple candidate genes and pathways associated with thyroid function, providing novel clues for revealing the genetic mechanisms of thyroid function and disease.

Fat metabolism is associated with telomere length in six population-based studies

Telomeres are repetitive DNA sequences located at the end of chromosomes, which are associated to biological aging, cardiovascular disease, cancer, and mortality. Lipid and fatty acid metabolism have been associated with telomere shortening. We have conducted an in-depth study investigating the association of metabolic biomarkers with telomere length (LTL). We performed an association analysis of 226 metabolic biomarkers with LTL using data from 11 775 individuals from six independent population-based cohorts (BBMRI-NL consortium). Metabolic biomarkers include lipoprotein lipids and subclasses, fatty acids, amino acids, glycolysis measures and ketone bodies. LTL was measured by quantitative polymerase chain reaction or FlowFISH. Linear regression analysis was performed adjusting for age, sex, lipid-lowering medication and cohort-specific covariates (model 1) and additionally for body mass index (BMI) and smoking (model 2), followed by inverse variance-weighted meta-analyses (significance threshold pmeta = 6.5x10-4). We identified four metabolic biomarkers positively associated with LTL, including two cholesterol to lipid ratios in small VLDL (S-VLDL-C % and S-VLDL-ce %) and two omega-6 fatty acid ratios (FAw6/FA and LA/FA). After additionally adjusting for BMI and smoking, these metabolic biomarkers remained associated with LTL with similar effect estimates. In addition, cholesterol esters in very small VLDL (XS-VLDL-ce) became significantly associated with LTL (p = 3.6x10-4). We replicated the association of FAw6/FA with LTL in an independent dataset of 7845 individuals (p = 1.9x10-4). To conclude, we identified multiple metabolic biomarkers involved in lipid and fatty acid metabolism that may be involved in LTL biology. Longitudinal studies are needed to exclude reversed causation.

Identical twins carry a persistent epigenetic signature of early genome programming

Monozygotic (MZ) twins and higher-order multiples arise when a zygote splits during pre-implantation stages of development. The mechanisms underpinning this event have remained a mystery. Because MZ twinning rarely runs in families, the leading hypothesis is that it occurs at random. Here, we show that MZ twinning is strongly associated with a stable DNA methylation signature in adult somatic tissues. This signature spans regions near telomeres and centromeres, Polycomb-repressed regions and heterochromatin, genes involved in cell-adhesion, WNT signaling, cell fate, and putative human metastable epialleles. Our study also demonstrates a never-anticipated corollary: because identical twins keep a lifelong molecular signature, we can retrospectively diagnose if a person was conceived as monozygotic twin.

Continuity of Genetic Risk for Aggressive Behavior Across the Life-Course

We test whether genetic influences that explain individual differences in aggression in early life also explain individual differences across the life-course. In two cohorts from The Netherlands (N = 13,471) and Australia (N = 5628), polygenic scores (PGSs) were computed based on a genome-wide meta-analysis of childhood/adolescence aggression. In a novel analytic approach, we ran a mixed effects model for each age (Netherlands: 12-70 years, Australia: 16-73 years), with observations at the focus age weighted as 1, and decaying weights for ages further away. We call this approach a 'rolling weights' model. In The Netherlands, the estimated effect of the PGS was relatively similar from age 12 to age 41, and decreased from age 41-70. In Australia, there was a peak in the effect of the PGS around age 40 years. These results are a first indication from a molecular genetics perspective that genetic influences on aggressive behavior that are expressed in childhood continue to play a role later in life.

The Link between Type 2 Diabetes Mellitus and the Polymorphisms of Glutathione-Metabolizing Genes Suggests a New Hypothesis Explaining Disease Initiation and Progression

The present study investigated whether type 2 diabetes (T2D) is associated with polymorphisms of genes encoding glutathione-metabolizing enzymes such as glutathione synthetase (GSS) and gamma-glutamyl transferase 7 (GGT7). A total of 3198 unrelated Russian subjects including 1572 T2D patients and 1626 healthy subjects were enrolled. Single nucleotide polymorphisms (SNPs) of the GSS and GGT7 genes were genotyped using the MassArray-4 system. We found that the GSS and GGT7 gene polymorphisms alone and in combinations are associated with T2D risk regardless of sex, age, and body mass index, as well as correlated with plasma glutathione, hydrogen peroxide, and fasting blood glucose levels. Polymorphisms of GSS (rs13041792) and GGT7 (rs6119534 and rs11546155) genes were associated with the tissue-specific expression of genes involved in unfolded protein response and the regulation of proteostasis. Transcriptome-wide association analysis has shown that the pancreatic expression of some of these genes such as EDEM2, MYH7B, MAP1LC3A, and CPNE1 is linked to the genetic risk of T2D. A comprehensive analysis of the data allowed proposing a new hypothesis for the etiology of type 2 diabetes that endogenous glutathione deficiency might be a key condition responsible for the impaired folding of proinsulin which triggered an unfolded protein response, ultimately leading to beta-cell apoptosis and disease development.

Longitudinal retinal layer changes in preclinical Alzheimer's disease

Purpose

Several studies found reduced retinal thickness on optical coherence tomography (OCT) in Alzheimer’s disease (AD), even in preclinical stages, labelling this technique of interest as biomarker. In this study, we examine retinal thickness changes in preclinical AD, as defined by cognitively normal individuals with amyloid‐beta (Aβ) on positron emission tomography (PET).

Methods

For this monocentre study, 145 cognitively healthy monozygotic twins aged ≥ 60 were included from the Netherlands Twin Register taking part in the EMIF‐AD PreclinAD study. At baseline, participants underwent [¹⁸F] flutemetamol PET that was visually rated for cortical Aβ. Binding potential was calculated as continuous measure for Aβ. Optical coherence tomography (OCT) was performed at baseline and after 22 months to assess changes in total and individual inner retinal layer thickness in the macular region (ETDRS circles) and peripapillary retinal nerve fibre layer thickness. Differences in rate of change between amyloid‐beta positive and negative individuals and associations between binding potential and change in retinal thickness were evaluated.

Results

Sixteen participants (11%) were positive for Aβ. Change in retinal thickness did not differ in any region between Aβ+ and Aβ− individuals. A positive association between binding potential and change in inner plexiform layer thickness was observed in the inner macular ring (beta = 1.708, CI = 0.575 to 2.841, p = 0.003).

Conclusion

Aβ+ individuals did not differ in rate of change of any retinal layer compared to controls, but higher binding potential at baseline was associated with less IPL thinning over time. Optical coherence tomography (OCT) as a longitudinal screening tool for preclinical AD seems limited, but IPL changes offer leads for further research.

DNA methylation signatures of aggression and closely related constructs: A meta-analysis of epigenome-wide studies across the lifespan

DNA methylation profiles of aggressive behavior may capture lifetime cumulative effects of genetic, stochastic, and environmental influences associated with aggression. Here, we report the first large meta-analysis of epigenome-wide association studies (EWAS) of aggressive behavior (N = 15,324 participants). In peripheral blood samples of 14,434 participants from 18 cohorts with mean ages ranging from 7 to 68 years, 13 methylation sites were significantly associated with aggression (alpha = 1.2 × 10-7; Bonferroni correction). In cord blood samples of 2425 children from five cohorts with aggression assessed at mean ages ranging from 4 to 7 years, 83% of these sites showed the same direction of association with childhood aggression (r = 0.74, p = 0.006) but no epigenome-wide significant sites were found. Top-sites (48 at a false discovery rate of 5% in the peripheral blood meta-analysis or in a combined meta-analysis of peripheral blood and cord blood) have been associated with chemical exposures, smoking, cognition, metabolic traits, and genetic variation (mQTLs). Three genes whose expression levels were associated with top-sites were previously linked to schizophrenia and general risk tolerance. At six CpGs, DNA methylation variation in blood mirrors variation in the brain. On average 44% (range = 3-82%) of the aggression-methylation association was explained by current and former smoking and BMI. These findings point at loci that are sensitive to chemical exposures with potential implications for neuronal functions. We hope these results to be a starting point for studies leading to applications as peripheral biomarkers and to reveal causal relationships with aggression and related traits.

Onset of Preclinical Alzheimer Disease in Monozygotic Twins

Objective

The present work was undertaken to study the genetic contribution to the start of Alzheimer's disease (AD) with amyloid and tau biomarkers in cognitively intact older identical twins.

Methods

We studied in 96 monozygotic twin‐pairs relationships between amyloid‐beta (Aβ) aggregation as measured by the Aβ1–42/1–40 ratio in cerebrospinal fluid (CSF; n = 126) and positron emission tomography (PET, n = 194), and CSF markers for Aβ production (beta‐secretase 1, Aβ1–40, and Aβ1–38) and CSF tau. Associations among markers were tested with generalized estimating equations including a random effect for twin status, adjusted for age, gender, and apolipoprotein E ε4 genotype. We used twin analyses to determine relative contributions of genetic and/or environmental factors to AD pathophysiological processes.

Results

Twenty‐seven individuals (14%) had an abnormal amyloid PET, and 14 twin‐pairs (15%) showed discordant amyloid PET scans. Within twin‐pairs, Aβ production markers and total‐tau (t‐tau) levels strongly correlated (r range = 0.73–0.86, all p < 0.0001), and Aβ aggregation markers and 181‐phosphorylated‐tau (p‐tau) levels correlated moderately strongly (r range = 0.50–0.64, all p < 0.0001). Cross‐twin cross‐trait analysis showed that Aβ1–38 in one twin correlated with Aβ1–42/1–40 ratios, and t‐tau and p‐tau levels in their cotwins (r range = −0.28 to 0.58, all p < .007). Within‐pair differences in Aβ production markers related to differences in tau levels (r range = 0.49–0.61, all p < 0.0001). Twin discordance analyses suggest that Aβ production and tau levels show coordinated increases in very early AD.

Interpretation

Our results suggest a substantial genetic/shared environmental background contributes to both Aβ and tau increases, suggesting that modulation of environmental risk factors may aid in delaying the onset of AD pathophysiological processes. ANN NEUROL 2021;89:987–1000

Obesity and atypical depression symptoms: findings from Mendelian randomization in two European cohorts

Studies considering the causal role of body mass index (BMI) for the predisposition of major depressive disorder (MDD) based on a Mendelian Randomization (MR) approach have shown contradictory results. These inconsistent findings may be attributable to the heterogeneity of MDD; in fact, several studies have documented associations between BMI and mainly the atypical subtype of MDD. Using a MR approach, we investigated the potential causal role of obesity in both the atypical subtype and its five specific symptoms assessed according to the Statistical Manual of Mental Disorders (DSM), in two large European cohorts, CoLaus|PsyCoLaus (n = 3350, 1461 cases and 1889 controls) and NESDA|NTR (n = 4139, 1182 cases and 2957 controls). We first tested general obesity measured by BMI and then the body fat distribution measured by waist-to-hip ratio (WHR). Results suggested that BMI is potentially causally related to the symptom increase in appetite, for which inverse variance weighted, simple median and weighted median MR regression estimated slopes were 0.68 (SE = 0.23, p = 0.004), 0.77 (SE = 0.37, p = 0.036), and 1.11 (SE = 0.39, p = 0.004). No causal effect of BMI or WHR was found on the risk of the atypical subtype or for any of the other atypical symptoms. Our findings show that higher obesity is likely causal for the specific symptom of increase in appetite in depressed participants and reiterate the need to study depression at the granular level of its symptoms to further elucidate potential causal relationships and gain additional insight into its biological underpinnings.

The Speed of Development of Adolescent Brain Age Depends on Sex and Is Genetically Determined

Children and adolescents show high variability in brain development. Brain age-the estimated biological age of an individual brain-can be used to index developmental stage. In a longitudinal sample of adolescents (age 9-23 years), including monozygotic and dizygotic twins and their siblings, structural magnetic resonance imaging scans (N = 673) at 3 time points were acquired. Using brain morphology data of different types and at different spatial scales, brain age predictors were trained and validated. Differences in brain age between males and females were assessed and the heritability of individual variation in brain age gaps was calculated. On average, females were ahead of males by at most 1 year, but similar aging patterns were found for both sexes. The difference between brain age and chronological age was heritable, as was the change in brain age gap over time. In conclusion, females and males show similar developmental ("aging") patterns but, on average, females pass through this development earlier. Reliable brain age predictors may be used to detect (extreme) deviations in developmental state of the brain early, possibly indicating aberrant development as a sign of risk of neurodevelopmental disorders.

Evolution, genes, and inter‐disciplinary personality research

Most commentaries welcomed an evolutionary genetic approach to personality, but several raised concerns about our integrative model. In response, we clarify the scientific status of evolutionary genetic theory and explain the plausibility and value of our evolutionary genetic model of personality, despite some shortcomings with the currently available theories and data. We also have a closer look at mate choice for personality traits, point to promising ways to assess evolutionarily relevant environmental factors and defend higher‐order personality domains and the g‐factor as the best units for evolutionary genetic analyses. Finally, we discuss which extensions of and alternatives to our model appear most fruitful, and end with a call for more inter‐disciplinary personality research grounded in evolutionary theory. Copyright © 2007 John Wiley & Sons, Ltd.

The Heritability of Type D Personality by an Extended Twin-Pedigree Analysis in the Netherlands Twin Register

Type D (Distressed) personality combines negative affectivity (NA) and social inhibition (SI) and is associated with an increased risk of cardiovascular disease. We aimed to (1) validate a new proxy based on the Achenbach System of Empirically Based Assessment (ASEBA) for Type D personality and its NA and SI subcomponents and (2) estimate the heritability of the Type D proxy in an extended twin-pedigree design in the Netherlands Twin Register (NTR). Proxies for the dichotomous Type D classification, and continuous NA, SI, and NAxSI (the continuous measure of Type D) scales were created based on 12 ASEBA items for 30,433 NTR participants (16,449 twins and 13,984 relatives from 11,106 pedigrees) and sources of variation were analyzed in the ‘Mendel’ software package. We estimated additive and non-additive genetic variance components, shared household and unique environmental variance components and ran bivariate models to estimate the genetic and non-genetic covariance between NA and SI. The Type D proxy showed good reliability and construct validity. The best fitting genetic model included additive and non-additive genetic effects with broad-sense heritabilities for NA, SI and NAxSI estimated at 49%, 50% and 49%, respectively. Household effects showed small contributions (4–9%) to the total phenotypic variation. The genetic correlation between NA and SI was .66 (reflecting both additive and non-additive genetic components). Thus, Type D personality and its NA and SI subcomponents are heritable, with a shared genetic basis for the two subcomponents.

The building blocks of social competence: Contributions of the Consortium of Individual Development

Social competence refers to the ability to engage in meaningful interactions with others. It is a crucial skill potentially malleable to interventions. Nevertheless, it remains difficult to select which children, which periods in a child's life, and which underlying skills form optimal targets for interventions. Development of social competence is complex to characterize because (a) it is by nature context- dependent; (b) it is subserved by multiple relevant processes that develop at different times in a child's life; and (c) over the years multiple, possibly conflicting, ways have been coined to index a child's social competence. The current paper elaborates upon a theoretical model of social competence developed by Rose-Krasnor (Rose- Krasnor, 1997; Rose-Krasnor and Denham, 2009), and it makes concrete how underlying skills and the variety of contexts of social interaction are both relevant dimensions of social competence that might change over development. It then illustrates how the cohorts and work packages in the Consortium on Individual Development each provide empirical contributions necessary for testing this model on the development of social competence.

Associations of Brain Pathology Cognitive and Physical Markers With Age in Cognitively Normal Individuals Aged 60-102 Years

The prevalence of brain pathologies increases with age and cognitive and physical functions worsen over the lifetime. It is unclear whether these processes show a similar increase with age. We studied the association of markers for brain pathology cognitive and physical functions with age in 288 cognitively normal individuals aged 60-102 years selected from the cross-sectional EMIF-AD PreclinAD and 90+ Study at the Amsterdam UMC. An abnormal score was consistent with a score below the 5th percentile in the 60- to 70-year-old individuals. Prevalence of abnormal scores was estimated using Generalized Estimating Equations (GEE) models. The prevalence of abnormal handgrip strength, the Digit Symbol Substitution Test, and hippocampal volume showed the fastest increase with age and abnormal MMSE score, muscle mass, and amyloid aggregation the lowest. The increase in prevalence of abnormal markers was partly dependent on sex, level of education, and amyloid aggregation. We did not find a consistent pattern in which markers of brain pathology cognitive and physical processes became abnormal with age.

Causes of Variation in Food Preference in the Netherlands

Our current society is characterized by an increased availability of industrially processed foods with high salt, fat and sugar content. How is it that some people prefer these unhealthy foods while others prefer more healthy foods? It is suggested that both genetic and environmental factors play a role. The aim of this study was to (1) identify food preference clusters in the largest twin-family study into food preference to date and (2) determine the relative contribution of genetic and environmental factors to individual differences in food preference in the Netherlands. Principal component analysis was performed to identify the preference clusters by using data on food liking/disliking from 16,541 adult multiples and their family members. To estimate the heritability of food preference, the data of 7833 twins were used in structural equation models. We identified seven food preference clusters (Meat, Fish, Fruits, Vegetables, Savory snacks, Sweet snacks and Spices) and one cluster with Drinks. Broad-sense heritability (additive [A] + dominant [D] genetic factors) for these clusters varied between .36 and .60. Dominant genetic effects were found for the clusters Fruit, Fish (males only) and Spices. Quantitative sex differences were found for Meat, Fish and Savory snacks and Drinks. To conclude, our study convincingly showed that genetic factors play a significant role in food preference. A next important step is to identify these genes because genetic vulnerability for food preference is expected to be linked to actual food consumption and different diet-related disorders.

Illicit drug use and the genetic overlap with Cannabis use

BACKGROUND

The use of illicit substances is correlated, meaning that individuals who use one illicit substance are more likely to also use another illicit substance. This association could (partly) be explained by overlapping genetic factors. Genetic overlap may indicate a common underlying genetic predisposition, or can be the result of a causal association.

METHODS

Polygenic scores for lifetime cannabis use were generated in a sample of Dutch participants (N = 8348). We tested the association of a PGS for cannabis use with ecstasy, stimulants and a broad category of illicit drug use. To explore the nature of the relationship: (1) these analyses were repeated separately in cannabis users and non-users and (2) monozogytic twin pairs discordant for cannabis use were compared on their drug use.

RESULTS

The lifetime prevalence was 24.8 % for cannabis, 6.2 % for ecstasy, 6.5 % for stimulants and 7.1 % for any illicit drug use. Significant, positive associations were found between PGS for cannabis use with ecstasy use, stimulants and any illicit drug use. These associations seemed to be stronger in cannabis users compared to non-users for both ecstasy and stimulant use, but only in people born after 1968 and not significant after correction for multiple testing. The discordant twin pair analyses suggested that cannabis use could play a causal role in drug use.

CONCLUSIONS

The genetic liability underlying cannabis use significantly explained variability in ecstasy, stimulant and any illicit drug use. Further research should further explore the underlying mechanism to understand the nature of the association.

Exploring the Temporal Relation between Body Mass Index and Corticosteroid Metabolite Excretion in Childhood

Childhood obesity is associated with alterations in hypothalamus–pituitary–adrenal (HPA) axis activity. However, it is unknown whether these alterations are a cause or a consequence of obesity. This study aimed to explore the temporal relationship between cortisol production and metabolism, and body mass index (BMI). This prospective follow-up study included 218 children (of whom 50% were male), born between 1995 and 1996, who were assessed at the ages of 9, 12 and 17 years. Morning urine samples were collected for assessment of cortisol metabolites by gas chromatography-tandem mass spectrometry, enabling the calculation of cortisol metabolite excretion rate and cortisol metabolic pathways. A cross-lagged regression model was used to determine whether BMI at various ages during childhood predicted later cortisol production and metabolism parameters, or vice versa. The cross-lagged regression coefficients showed that BMI positively predicted cortisol metabolite excretion (p = 0.03), and not vice versa (p = 0.33). In addition, BMI predicted the later balance of 11β-hydroxysteroid dehydrogenase (HSD) activities (p = 0.07), and not vice versa (p = 0.55). Finally, cytochrome P450 3A4 activity positively predicted later BMI (p = 0.01). Our study suggests that changes in BMI across the normal range predict alterations in HPA axis activity. Therefore, the alterations in HPA axis activity as observed in earlier studies among children with obesity may be a consequence rather than a cause of increased BMI.

Integration of transcriptome-wide association study and messenger RNA expression profile to identify genes associated with osteoarthritis

Aims

Osteoarthritis (OA) is the most prevalent joint disease. However, the specific and definitive genetic mechanisms of OA are still unclear.

Methods

Tissue-related transcriptome-wide association studies (TWAS) of hip OA and knee OA were performed utilizing the genome-wide association study (GWAS) data of hip OA and knee OA (including 2,396 hospital-diagnosed hip OA patients versus 9,593 controls, and 4,462 hospital-diagnosed knee OA patients versus 17,885 controls) and gene expression reference to skeletal muscle and blood. The OA-associated genes identified by TWAS were further compared with the differentially expressed genes detected by the messenger RNA (mRNA) expression profiles of hip OA and knee OA. Functional enrichment and annotation analysis of identified genes was performed by the DAVID and FUMAGWAS tools.

Results

We detected 33 common genes, eight common gene ontology (GO) terms, and one common pathway for hip OA, such as calcium and integrin-binding protein 1 (CIB1) (P_TWAS = 0.025, FC_mRNA = -1.575 for skeletal muscle), adrenomedullin (ADM) (P_TWAS = 0.022, FC_mRNA = -4.644 for blood), Golgi apparatus (P_TWAS <0.001, P_mRNA = 0.012 for blood), and phosphatidylinositol 3' -kinase-protein kinase B (PI3K-Akt) signalling pathway (P_TWAS = 0.033, P_mRNA = 0.005 for blood). For knee OA, we detected 24 common genes, eight common GO terms, and two common pathways, such as histocompatibility complex, class II, DR beta 1 (HLA-DRB1) (P_TWAS = 0.040, FC_mRNA = 4.062 for skeletal muscle), Follistatin-like 1 (FSTL1) (P_TWAS = 0.048, FC_mRNA = 3.000 for blood), cytoplasm (P_TWAS < 0.001, P_mRNA = 0.005 for blood), and complement and coagulation cascades (P_TWAS = 0.017, P_mRNA = 0.001 for skeletal muscle).

Conclusion

We identified a group of OA-associated genes and pathways, providing novel clues for understanding the genetic mechanism of OA.

Cite this article:Bone Joint Res. 2020;9(3):130–138.

Ocular biomarkers for cognitive impairment in nonagenarians; a prospective cross-sectional study

Background

Ocular imaging receives much attention as a source of potential biomarkers for dementia. In the present study, we analyze these ocular biomarkers in cognitively impaired and healthy participants in a population aged over 90 years (= nonagenarian), and elucidate the effects of age on these biomarkers.

Methods

For this prospective cross-sectional study, we included individuals from the EMIF-AD 90+ study, consisting of a cognitively healthy (N = 67) and cognitively impaired group (N = 33), and the EMIF-AD PreclinAD study, consisting of cognitively healthy controls aged ≥60 (N = 198). Participants underwent Optical Coherence Tomography (OCT) and fundus photography of both eyes. OCT was used to asses total and individual inner retinal layer thickness in the macular region (Early Treatment Diabetic Retinopathy Study circles) as well as peripapillary retinal nerve fiber layer thickness, fundus images were analyzed with Singapore I Vessel Assessment to obtain 7 retinal vascular parameters. Values for both eyes were averaged. Differences in ocular biomarkers between the 2 nonagenarian groups were analyzed using linear regression, differences between the individual nonagenarian groups and controls were analyzed using generalized estimating equations.

Results

Ocular biomarkers did not differ between the healthy and cognitively impaired nonagenarian groups. 19 out of 22 ocular biomarkers assessed in this study differed between either nonagenarian group and the younger controls.

Conclusion

The ocular biomarkers assessed in this study were not associated with cognitive impairment in nonagenarians, making their use as a screening tool for dementing disorders in this group limited. However, ocular biomarkers were significantly associated with chronological age, which were very similar to those ascribed to occur in Alzheimer’s Disease.

Predicting Loneliness from Where and What People Do

The many devastating mental health outcomes associated with chronic loneliness is the motivation behind research into examining personal and demographic characteristics of the lonely. The present study sought to examine the connection of where people live (degree of urbanization) and what people do (leisure activities) with self-report of loneliness in a large sample (N = 8356) of unrelated Dutch adults. Information regarding where people live and what they do in their leisure time was entered into a regression analysis for self-reported loneliness. The overall regression was significant and accounted for 2.8% of the loneliness scale scores. Significant independent predictors for loneliness were living in heavily urbanized areas and engaging in fewer social activities. People who went sightseeing or to amusement parks/zoos or who participated in clubs reported being less lonely. Spending time using a computer predicted higher self-report loneliness scores. Consistent with previous research, after controlling for other variables, gender was not a significant predictor of loneliness but both a younger age and a curvilinear or U-shaped curve of age predicted loneliness (the younger and the much older). The results suggest that meaningful interpersonal interactions may result in lower feelings of loneliness.

Genetic Vulnerability for Smoking and Cannabis Use: Associations With E-Cigarette and Water Pipe Use

INTRODUCTION

Cigarette smoking and cannabis use are heritable traits and share, at least in part, a common genetic substrate. In recent years, the prevalence of alternative methods of nicotine intakes, such as electronic cigarette (e-cigarette) and water pipe use, has risen substantially. We tested whether the genetic vulnerability underlying cigarettes smoking and cannabis use explained variability in e-cigarette and water pipe use phenotypes, as these vaping methods are alternatives for smoking tobacco cigarettes and joints.

METHODS

On the basis of the summary statistics of the International Cannabis Consortium and the Tobacco and Genetics Consortium, we generated polygenic risk scores (PRSs) for smoking and cannabis use traits, and used these to predict e-cigarette and water pipe use phenotypes in a sample of 5025 individuals from the Netherlands Twin Register.

RESULTS

PRSs for cigarettes per day were positively associated with lifetime e-cigarette use and early initiation of water pipe use, but only in ex-smokers (odds ratio = 1.43, R2 = 1.56%, p = .011) and never cigarette smokers (odds ratio = 1.35, R2 = 1.60%, p = .013) respectively.

CONCLUSIONS

Most associations of PRSs for cigarette smoking and cannabis use with e-cigarette and water pipe use were not significant, potentially due to a lack of power. The significant associations between genetic liability to smoking heaviness with e-cigarette and water pipe phenotypes are in line with studies indicating a common genetic background for substance-use phenotypes. These associations emerged only in nonsmokers, and future studies should investigate the nature of this observation.

IMPLICATIONS

Our study showed that genetic vulnerability to smoking heaviness is associated with lifetime e-cigarette use and age at initiation of water pipe use. This finding has implications for the current debate on whether alternative smoking methods, such as usage of vaping devices, predispose to smoking initiation and related behaviors.

Urinary Amine and Organic Acid Metabolites Evaluated as Markers for Childhood Aggression: The ACTION Biomarker Study

Biomarkers are of interest as potential diagnostic and predictive instruments in personalized medicine. We present the first urinary metabolomics biomarker study of childhood aggression. We aim to examine the association of urinary metabolites and neurotransmitter ratios involved in key metabolic and neurotransmitter pathways in a large cohort of twins (N = 1,347) and clinic-referred children (N = 183) with an average age of 9.7 years. This study is part of ACTION (Aggression in Children: Unraveling gene-environment interplay to inform Treatment and InterventiON strategies), in which we developed a standardized protocol for large-scale collection of urine samples in children. Our analytical design consisted of three phases: a discovery phase in twins scoring low or high on aggression (N = 783); a replication phase in twin pairs discordant for aggression (N = 378); and a validation phase in clinical cases and matched twin controls (N = 367). In the discovery phase, 6 biomarkers were significantly associated with childhood aggression, of which the association of O-phosphoserine (β = 0.36; SE = 0.09; p = 0.004), and gamma-L-glutamyl-L-alanine (β = 0.32; SE = 0.09; p = 0.01) remained significant after multiple testing. Although non-significant, the directions of effect were congruent between the discovery and replication analyses for six biomarkers and two neurotransmitter ratios and the concentrations of 6 amines differed between low and high aggressive twins. In the validation analyses, the top biomarkers and neurotransmitter ratios, with congruent directions of effect, showed no significant associations with childhood aggression. We find suggestive evidence for associations of childhood aggression with metabolic dysregulation of neurotransmission, oxidative stress, and energy metabolism. Although replication is required, our findings provide starting points to investigate causal and pleiotropic effects of these dysregulations on childhood aggression.

Long-Term Stability of Cortisol Production and Metabolism Throughout Adolescence: Longitudinal Twin Study

Life-course experiences have been postulated to program hypothalamus-pituitary-adrenal (HPA) axis activity, suggesting that HPA axis activity is, at least partially, stable over time. Yet, there is paucity of data on the long-term stability of cortisol production and metabolism. We performed a prospective follow-up study in twins recruited from a nationwide register to estimate the stability of cortisol production and metabolism over time, and the contribution of genetic and environmental factors to this stability. In total, 218 healthy mono- and dizygotic twins were included. At the ages of 9, 12 and 17 years, morning urine samples were collected for assessment (by gas chromatography-tandem mass spectrometry) of cortisol metabolites, enabling the calculation of cortisol metabolite excretion rate and cortisol metabolism activity. Our results showed a low stability for both cortisol metabolite excretion rate (with correlations <.20) and cortisol metabolism activity indices (with correlations of .25 to .46 between 9 and 12 years, -.02 to .15 between 12 and 17 years and .09 to .28 between 9 and 17 years). Because of the low stability over time, genetic and environmental contributions to this stability were difficult to assess, although it seemed to be mostly determined by genetic factors. The low stability in both cortisol production and metabolism between ages 9 and 17 years reflects the dynamic nature of the HPA axis.

Metabolomics Profile in Depression: A Pooled Analysis of 230 Metabolic Markers in 5283 Cases With Depression and 10,145 Controls

BACKGROUND

Depression has been associated with metabolic alterations, which adversely impact cardiometabolic health. Here, a comprehensive set of metabolic markers, predominantly lipids, was compared between depressed and nondepressed persons.

METHODS

Nine Dutch cohorts were included, comprising 10,145 control subjects and 5283 persons with depression, established with diagnostic interviews or questionnaires. A proton nuclear magnetic resonance metabolomics platform provided 230 metabolite measures: 51 lipids, fatty acids, and low-molecular-weight metabolites; 98 lipid composition and particle concentration measures of lipoprotein subclasses; and 81 lipid and fatty acids ratios. For each metabolite measure, logistic regression analyses adjusted for gender, age, smoking, fasting status, and lipid-modifying medication were performed within cohort, followed by random-effects meta-analyses.

RESULTS

Of the 51 lipids, fatty acids, and low-molecular-weight metabolites, 21 were significantly related to depression (false discovery rate q < .05). Higher levels of apolipoprotein B, very-low-density lipoprotein cholesterol, triglycerides, diglycerides, total and monounsaturated fatty acids, fatty acid chain length, glycoprotein acetyls, tyrosine, and isoleucine and lower levels of high-density lipoprotein cholesterol, acetate, and apolipoprotein A1 were associated with increased odds of depression. Analyses of lipid composition indicators confirmed a shift toward less high-density lipoprotein and more very-low-density lipoprotein and triglyceride particles in depression. Associations appeared generally consistent across gender, age, and body mass index strata and across cohorts with depressive diagnoses versus symptoms.

CONCLUSIONS

This large-scale meta-analysis indicates a clear distinctive profile of circulating lipid metabolites associated with depression, potentially opening new prevention or treatment avenues for depression and its associated cardiometabolic comorbidity.

Optical coherence tomography angiography in preclinical Alzheimer's disease

BACKGROUND/AIMS

As a protrusion from the brain, the retina might reflect the status of the brain. Previous studies showed a decrease in vessel density and foveal avascular zone (FAZ) enlargement on optical coherence tomography angiography (OCTA) in individuals suffering from Alzheimer's disease (AD). This study aims to assess whether such changes are already present in preclinical stages of AD, in a population of monozygotic (MZ) twins.

METHODS

124 cognitively healthy individuals (MZ twins, ages 60-93 years) underwent [18F]flutemetamol amyloid positron emission tomography (PET) scanning and OCTA. PET scans were visually rated for cortical amyloid-beta (Aβ) positivity. Parametric global cortical non-displaceable binding potential (BPND) was used as a continuous measure for Aβ aggregation. FAZ size and vessel densities for the inner and outer ring of the macular ETDRS grid and in a 3-6 mm ring around the optic nerve head (ONH) were measured.OCTA measures were associated with visual Aβ score, BPND and amyloid load estimated by twin concordance on visual Aβ score. Twin correlations were estimated as a measure of maximum heritability of OCTA measures.

RESULTS

13 of 124 participants were Aβ+. Aβ+ individuals had significantly higher vessel density than Aβ- individuals in all regions but did not differ in FAZ size. Twin analyses showed a positive association between and vessel densities in all regions. BPND tended to be associated with higher vessel density in the inner ring. Twin correlations were moderate/high for all OCTA parameters except vessel density around the ONH, which correlated weakly.

CONCLUSION

Retinal vessel density was higher in individuals with preclinical AD.

Heritability of Cortisol Production and Metabolism Throughout Adolescence

CONTEXT

Inter-individual differences in cortisol production and metabolism emerge with age and may be explained by genetic factors.

OBJECTIVE

To estimate the relative contributions of genetic and environmental factors to inter-individual differences in cortisol production and metabolism throughout adolescence.

DESIGN

Prospective follow-up study of twins.

SETTING

Nationwide register.

PARTICIPANTS

218 mono- and dizygotic twins (N = 109 pairs) born between 1995 amd 1996, recruited from the Netherlands Twin Register. Cortisol metabolites were determined in 213, 169, and 160 urine samples at the ages of 9, 12, and 17, respectively.

MAIN OUTCOME MEASURES

The total contribution of genetic factors (broad-sense heritability) and shared and unshared environmental influences to inter-individual differences in cortisol production and activities of 5α-reductase, 5β-reductase, and 11β-hydroxysteroid dehydrogenases and cytochrome P450 3A4.

RESULTS

For cortisol production rate at the ages of 9, 12, and 17, broad-sense heritability was estimated as 42%, 30%, and 0%, respectively, and the remainder of the variance was explained by unshared environmental factors. For cortisol metabolism indices, the following heritability was observed: for the A-ring reductases (5α-and 5β-reductases), broad-sense heritability increased with age (to >50%), while for the other indices (renal 11β-HSD2, global 11β-HSD, and CYP3A4), the contribution of genetic factors was highest (68%, 18%, and 67%, respectively) at age 12.

CONCLUSIONS

The contribution of genetic factors to inter-individual differences in cortisol production decreased between 12 and 17y, indicative of a predominant role of individual circumstances. For cortisol metabolism, distinct patterns of genetic and environmental influences were observed, with heritability that either increased with age or peaked at age 12y.

Integration of epidemiologic, pharmacologic, genetic and gut microbiome data in a drug-metabolite atlas

Progress in high-throughput metabolic profiling provides unprecedented opportunities to obtain insights into the effects of drugs on human metabolism. The Biobanking BioMolecular Research Infrastructure of the Netherlands has constructed an atlas of drug-metabolite associations for 87 commonly prescribed drugs and 150 clinically relevant plasma-based metabolites assessed by proton nuclear magnetic resonance. The atlas includes a meta-analysis of ten cohorts (18,873 persons) and uncovers 1,071 drug-metabolite associations after evaluation of confounders including co-treatment. We show that the effect estimates of statins on metabolites from the cross-sectional study are comparable to those from intervention and genetic observational studies. Further data integration links proton pump inhibitors to circulating metabolites, liver function, hepatic steatosis and the gut microbiome. Our atlas provides a tool for targeted experimental pharmaceutical research and clinical trials to improve drug efficacy, safety and repurposing. We provide a web-based resource for visualization of the atlas (http://bbmri.researchlumc.nl/atlas/).

Heritability estimates for 361 blood metabolites across 40 genome-wide association studies

Metabolomics examines the small molecules involved in cellular metabolism. Approximately 50% of total phenotypic differences in metabolite levels is due to genetic variance, but heritability estimates differ across metabolite classes. We perform a review of all genome-wide association and (exome-) sequencing studies published between November 2008 and October 2018, and identify >800 class-specific metabolite loci associated with metabolite levels. In a twin-family cohort (N = 5117), these metabolite loci are leveraged to simultaneously estimate total heritability (h2total), and the proportion of heritability captured by known metabolite loci (h2Metabolite-hits) for 309 lipids and 52 organic acids. Our study reveals significant differences in h2Metabolite-hits among different classes of lipids and organic acids. Furthermore, phosphatidylcholines with a high degree of unsaturation have higher h2Metabolite-hits estimates than phosphatidylcholines with low degrees of unsaturation. This study highlights the importance of common genetic variants for metabolite levels, and elucidates the genetic architecture of metabolite classes.

Retinal layer thickness in preclinical Alzheimer's disease

Purpose

There is urgent need for non‐invasive diagnostic biomarkers in the preclinical phase of Alzheimer's Disease (AD). Several studies suggest that retinal thickness is reduced in AD. Here, we aim to test the diagnostic value of retinal thickness in preclinical AD, as defined by cognitively normal individuals with amyloid pathology on PET.

Methods

One hundred and sixty five cognitively healthy monozygotic twins aged ≥ 60 were included from the Netherlands Twin Register taking part in the European Medical Information Framework for Alzheimer's Disease PreclinAD study. Participants underwent [¹⁸F] flutemetamol PET that was visually rated for presence or absence of cortical amyloid beta (Aβ). Binding potential (BP_ND) was calculated as continuous measure for Aβ. Spectral Domain OCT was used to asses total and individual inner retinal layer thickness in the macular region (ETDRS circles) as well as peripapillary retinal nerve fibre layer (pRNFL) thickness. Differences between Aβ+ and Aβ− individuals and associations between BP_ND and retinal thickness were analyzed.

Results

No differences were found in retinal layer thickness in the macula or pRNFL between Aβ+ and Aβ− individuals. A positive associations between BP_ND and macular total retinal thickness was observed in the inner ring (p = 0.018), but this was not statistically significant after correction for multiple testing (p = 0.144). Brain/eye parameters had moderate to high intra‐twin correlations (p < 0.001) except visual rating score of Aβ, which did not correlate (r = 0.21, p = 0.068).

Conclusion

Variation in retinal thickness likely reflects genetic differences between individuals, but cannot discriminate between healthy and preclinical AD cases, making its use as biomarker in these early stages limited.

The Netherlands Twin Register: Longitudinal Research Based on Twin and Twin-Family Designs

The Netherlands Twin Register (NTR) is a national register in which twins, multiples and their parents, siblings, spouses and other family members participate. Here we describe the NTR resources that were created from more than 30 years of data collections; the development and maintenance of the newly developed database systems, and the possibilities these resources create for future research. Since the early 1980s, the NTR has enrolled around 120,000 twins and a roughly equal number of their relatives. The majority of twin families have participated in survey studies, and subsamples took part in biomaterial collection (e.g., DNA) and dedicated projects, for example, for neuropsychological, biomarker and behavioral traits. The recruitment into the NTR is all inclusive without any restrictions on enrollment. These resources - the longitudinal phenotyping, the extended pedigree structures and the multigeneration genotyping - allow for future twin-family research that will contribute to gene discovery, causality modeling, and studies of genetic and cultural inheritance.

The German Twin Family Panel (TwinLife)

The German Twin Family Panel (TwinLife) is a German longitudinal study of monozygotic and dizygotic same-sex twin pairs and their families that was designed to investigate the development of social inequalities over the life course. The study covers an observation period from approximately 2014 to 2023. The target population of the sample are reared-together twins of four different age cohorts that were born in 2009/2010 (cohort 1), in 2003/2004 (cohort 2), in 1997/1998 (cohort 3) and between 1990 and 1993 (cohort 4). In the first wave, the study included data on 4097 twin families. Families were recruited in all parts of Germany so that the sample comprises the whole range of the educational, occupational and income structure. As of 2019, two face-to-face, at-home interviews and two telephone interviews have been conducted. Data from the first home and telephone interviews are already available free of charge as a scientific use-file from the GESIS data archive. This report aims to provide an overview of the study sample and design as well as constructs that are unique in TwinLife in comparison with previous twin studies - such as an assessment of cognitive abilities or information based on the children's medical records and report cards. In addition, major findings based on the data already released are displayed, and future directions of the study are presented and discussed.

Genetic Similarity Assessment of Twin-Family Populations by Custom-Designed Genotyping Array

Twin registries often take part in large collaborative projects and are major contributors to genome-wide association (GWA) meta-analysis studies. In this article, we describe genotyping of twin-family populations from Australia, the Midwestern USA (Avera Twin Register), the Netherlands (Netherlands Twin Register), as well as a sample of mothers of twins from Nigeria to assess the extent, if any, of genetic differences between them. Genotyping in all cohorts was done using a custom-designed Illumina Global Screening Array (GSA), optimized to improve imputation quality for population-specific GWA studies. We investigated the degree of genetic similarity between the populations using several measures of population variation with genotype data generated from the GSA. Visualization of principal component analysis (PCA) revealed that the Australian, Dutch and Midwestern American populations exhibit negligible interpopulation stratification when compared to each other, to a reference European population and to globally distant populations. Estimations of fixation indices (FST values) between the Australian, Midwestern American and Netherlands populations suggest minimal genetic differentiation compared to the estimates between each population and a genetically distinct cohort (i.e., samples from Nigeria genotyped on GSA). Thus, results from this study demonstrate that genotype data from the Australian, Dutch and Midwestern American twin-family populations can be reasonably combined for joint-genetic analysis.

Genetic and environmental influences on functional connectivity within and between canonical cortical resting-state networks throughout adolescent development in boys and girls

The human brain is active during rest and hierarchically organized into intrinsic functional networks. These functional networks are largely established early in development, with reports of a shift from a local to more distributed organization during childhood and adolescence. It remains unknown to what extent genetic and environmental influences on functional connectivity change throughout adolescent development. We measured functional connectivity within and between eight cortical networks in a longitudinal resting-state fMRI study of adolescent twins and their older siblings on two occasions (mean ages 13 and 18 years). We modelled the reliability for these inherently noisy and head-motion sensitive measurements by analyzing data from split-half sessions. Functional connectivity between resting-state networks decreased with age whereas functional connectivity within resting-state networks generally increased with age, independent of general cognitive functioning. Sex effects were sparse, with stronger functional connectivity in the default mode network for girls compared to boys, and stronger functional connectivity in the salience network for boys compared to girls. Heritability explained up to 53% of the variation in functional connectivity within and between resting-state networks, and common environment explained up to 33%. Genetic influences on functional connectivity remained stable during adolescent development. In conclusion, longitudinal age-related changes in functional connectivity within and between cortical resting-state networks are subtle but wide-spread throughout adolescence. Genes play a considerable role in explaining individual variation in functional connectivity with mostly stable influences throughout adolescence.

Integrating transcriptome-wide association study and mRNA expression profiling identifies novel genes associated with bone mineral density

To scan novel candidate genes associated with osteoporosis, a two-stage transcriptome-wide association study (TWAS) of bone mineral density (BMD) was conducted. The BMD-associated genes identified by TWAS were then compared with the gene expression profiling of BMD in bone cells, B cells, and mesenchymal stem cells. We identified multiple candidate genes and gene ontology (GO) terms associated with BMD.

INTRODUCTION

Osteoporosis (OP) is a metabolic bone disease characterized by decrease in BMD. Our objective is to scan novel candidate genes associated with OP.

METHODS

A transcriptome-wide association study (TWAS) was performed by integrating the genome-wide association study (GWAS) summary of bone mineral density (BMD) with two pre-computed mRNA expression weights of peripheral blood and muscle skeleton. Then, another independent GWAS data of BMD was used to verify the discovery results. The BMD-associated genes identified between discovery and replicate TWAS were further subjected to gene ontology (GO) analysis implemented by DAVID. Finally, the BMD-associated genes and GO terms were further compared with the mRNA expression profiling results of BMD to detect the common genes and GO terms shared by both DNA-level TWAS and mRNA expression profile analysis.

RESULTS

TWAS identified 95 common genes with permutation P value < 0.05 for peripheral blood and muscle skeleton, such as TMTC4 in muscle skeleton and DDX17 in peripheral blood. Further comparing the genes detected by discovery-replicate TWAS with the differentially expressed genes identified by mRNA expression profiling of OP patients found 18 overlapped genes, such as MUL1 in muscle skeleton and SPTBN1 in peripheral blood. GO analysis of the genes identified by discovery-replicate TWAS detected 12 BMD-associated GO terms, such as negative regulation of cell growth and regulation of glycogen catabolic process. Further comparing the GO results of discovery-replicate TWAS and mRNA expression profiles found 6 overlapped GO terms, such as membrane and cell adhesion.

CONCLUSION

Our study identified multiple candidate genes and GO terms for BMD, providing novel clues for understanding the genetic mechanism of OP.

Men Have a Stronger Monocyte-Derived Cytokine Production Response upon Stimulation with the Gram-Negative Stimulus Lipopolysaccharide than Women: A Pooled Analysis Including 15 Study Populations

The incidence of bacterial infections and sepsis, as well as the mortality risk from sepsis, is sex specific. These clinical findings have been attributed to sex differences in immune responsiveness. The aim of the present study was to investigate sex differences in monocyte-derived cytokine production response upon stimulation with the gram-negative stimulus lipopolysaccharide (LPS) using cytokine data from 15 study populations. Individual data on ex vivo cytokine production response upon stimulation with LPS in whole blood were available for 4,020 subjects originating from these 15 study populations, either from the general population or from patient populations with specific diseases. Men had a stronger cytokine production response than women to LPS for tumour necrosis factor-α, interleukin (IL)-6, IL-12, IL-1β, IL-1RA, and IL-10, but not for interferon-γ. The granulocyte-macrophage colony-stimulating factor production response was lower in men than in women. These sex differences were independent of chronological age. As men had higher monocyte concentrations, we normalized the cytokine production responses for monocyte concentration. After normalization, the sex differences in cytokine production response to LPS disappeared, except for IL-10, for which the production response was lower in men than in women. A sex-based approach to interpreting immune responsiveness is crucial.

Retinal and Cerebral Microvasculopathy: Relationships and Their Genetic Contributions

Purpose

Retinal microvasculopathy may reflect small vessel disease in the brain. Here we test the relationships between retinal vascular parameters and small vessel disease, the influence of cardiovascular risk factors on these relationships, and their common genetic background in a monozygotic twin cohort.

Methods

We selected 134 cognitively healthy individuals (67 monozygotic twin pairs) aged ≥60 years from the Netherlands Twin Register for the EMIF-AD PreclinAD study. We measured seven retinal vascular parameters averaged over both eyes using fundus images analyzed with Singapore I Vessel Assessment. Small vessel disease was assessed on MRI by a volumetric measurement of periventricular and deep white matter hyperintensities. We calculated associations between RVPs and WMH, estimated intratwin pair correlations, and performed twin-specific analyses on relationships of interest.

Results

Deep white matter hyperintensities volume was positively associated with retinal tortuosity in veins (P = 0.004) and fractal dimension in arteries (P = 0.001) and veins (P = 0.032), periventricular white matter hyperintensities volume was positively associated with retinal venous width (P = 0.028). Intratwin pair correlations were moderate to high for all small vessel disease/retinal vascular parameter variables (r = 0.49-0.87, P < 0.001). Cross-twin cross-trait analyses showed that retinal venous tortuosity of twin 1 could predict deep white matter hyperintensities volume of the co-twin (r = 0.23, P = 0.030). Within twin-pair differences for retinal venous tortuosity were associated with within twin-pair differences in deep white matter hyperintensities volume (r = 0.39, P = 0.001).

Conclusions

Retinal arterial fractal dimension and venous tortuosity have associations with deep white matter hyperintensities volume. Twin-specific analyses suggest that retinal venous tortuosity and deep white matter hyperintensities volume have a common etiology driven by both shared genetic factors and unique environmental factors, supporting the robustness of this relationship.