A multivariate genome-wide association study of psycho-cardiometabolic multimorbidity.

Coronary artery disease (CAD), type 2 diabetes (T2D) and depression are among the leading causes of chronic morbidity and mortality worldwide. Epidemiological studies indicate a substantial degree of multimorbidity, which may be explained by shared genetic influences. However, research exploring the presence of pleiotropic variants and genes common to CAD, T2D and depression is lacking. The present study aimed to identify genetic variants with effects on cross-trait liability to psycho-cardiometabolic diseases. We used genomic structural equation modelling to perform a multivariate genome-wide association study of multimorbidity (Neffective = 562,507), using summary statistics from univariate genome-wide association studies for CAD, T2D and major depression. CAD was moderately genetically correlated with T2D (rg = 0.39, P = 2e-34) and weakly correlated with depression (rg = 0.13, P = 3e-6). Depression was weakly correlated with T2D (rg = 0.15, P = 4e-15). The latent multimorbidity factor explained the largest proportion of variance in T2D (45%), followed by CAD (35%) and depression (5%). We identified 11 independent SNPs associated with multimorbidity and 18 putative multimorbidity-associated genes. We observed enrichment in immune and inflammatory pathways. A greater polygenic risk score for multimorbidity in the UK Biobank (N = 306,734) was associated with the co-occurrence of CAD, T2D and depression (OR per standard deviation = 1.91, 95% CI = 1.74-2.10, relative to the healthy group), validating this latent multimorbidity factor. Mendelian randomization analyses suggested potentially causal effects of BMI, body fat percentage, LDL cholesterol, total cholesterol, fasting insulin, income, insomnia, and childhood maltreatment. These findings advance our understanding of multimorbidity suggesting common genetic pathways.

The metabolome-wide signature of major depressive disorder.

Major Depressive Disorder (MDD) is a common, frequently chronic condition characterized by substantial molecular alterations and pathway dysregulations. Single metabolite and targeted metabolomics platforms have revealed several metabolic alterations in depression, including energy metabolism, neurotransmission, and lipid metabolism. More comprehensive coverage of the metabolome is needed to further specify metabolic dysregulations in depression and reveal previously untargeted mechanisms. Here, we measured 820 metabolites using the metabolome-wide Metabolon platform in 2770 subjects from a large Dutch clinical cohort with extensive clinical phenotyping (1101 current MDD, 868 remitted MDD, 801 healthy controls) at baseline, which were repeated in 1805 subjects at 6-year follow up (327 current MDD, 1045 remitted MDD, 433 healthy controls). MDD diagnosis was based on DSM-IV psychiatric interviews. Depression severity was measured with the Inventory of Depressive Symptomatology Self-report. Associations between metabolites and MDD status and depression severity were assessed at baseline and at 6-year follow-up. At baseline, 139 and 126 metabolites were associated with current MDD status and depression severity, respectively, with 79 overlapping metabolites. Adding body mass index and lipid-lowering medication to the models changed results only marginally. Among the overlapping metabolites, 34 were confirmed in internal replication analyses using 6-year follow-up data. Downregulated metabolites were enriched with long-chain monounsaturated (P = 6.7e-07) and saturated (P = 3.2e-05) fatty acids; upregulated metabolites were enriched with lysophospholipids (P = 3.4e-4). Mendelian randomization analyses using genetic instruments for metabolites (N = 14,000) and MDD (N = 800,000) showed that genetically predicted higher levels of the lysophospholipid 1-linoleoyl-GPE (18:2) were associated with greater risk of depression. The identified metabolome-wide profile of depression indicated altered lipid metabolism with downregulation of long-chain fatty acids and upregulation of lysophospholipids, for which causal involvement was suggested using genetic tools. This metabolomics signature offers a window on depression pathophysiology and a potential access point for the development of novel therapeutic approaches.

Discontinuation of psychotropic medication: a synthesis of evidence across medication classes.

Pharmacotherapy is an effective treatment modality across psychiatric disorders. Nevertheless, many patients discontinue their medication at some point. Evidence-based guidance for patients, clinicians, and policymakers on rational discontinuation strategies is vital to enable the best, personalized treatment for any given patient. Nonetheless, there is a scarcity of guidelines on discontinuation strategies. In this perspective, we therefore summarize and critically appraise the evidence on discontinuation of six major psychotropic medication classes: antidepressants, antipsychotics, benzodiazepines, mood stabilizers, opioids, and stimulants. For each medication class, a wide range of topics pertaining to each of the following questions are discussed: (1) Who can discontinue (e.g., what are risk factors for relapse?); (2) When to discontinue (e.g., after 1 year or several years of antidepressant use?); and (3) How to discontinue (e.g., what's the efficacy of dose reduction compared to full cessation and interventions to mitigate relapse risk?). We thus highlight how comparing the evidence across medication classes can identify knowledge gaps, which may pave the way for more integrated research on discontinuation.

Psychosocial factors, health behaviors and risk of cancer incidence: Testing interaction and effect modification in an individual participant data meta-analysis.

Depression, anxiety and other psychosocial factors are hypothesized to be involved in cancer development. We examined whether psychosocial factors interact with or modify the effects of health behaviors, such as smoking and alcohol use, in relation to cancer incidence. Two-stage individual participant data meta-analyses were performed based on 22 cohorts of the PSYchosocial factors and CAncer (PSY-CA) study. We examined nine psychosocial factors (depression diagnosis, depression symptoms, anxiety diagnosis, anxiety symptoms, perceived social support, loss events, general distress, neuroticism, relationship status), seven health behaviors/behavior-related factors (smoking, alcohol use, physical activity, body mass index, sedentary behavior, sleep quality, sleep duration) and seven cancer outcomes (overall cancer, smoking-related, alcohol-related, breast, lung, prostate, colorectal). Effects of the psychosocial factor, health behavior and their product term on cancer incidence were estimated using Cox regression. We pooled cohort-specific estimates using multivariate random-effects meta-analyses. Additive and multiplicative interaction/effect modification was examined. This study involved 437,827 participants, 36,961 incident cancer diagnoses, and 4,749,481 person years of follow-up. Out of 744 combinations of psychosocial factors, health behaviors, and cancer outcomes, we found no evidence of interaction. Effect modification was found for some combinations, but there were no clear patterns for any particular factors or outcomes involved. In this first large study to systematically examine potential interaction and effect modification, we found no evidence for psychosocial factors to interact with or modify health behaviors in relation to cancer incidence. The behavioral risk profile for cancer incidence is similar in people with and without psychosocial stress.

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.5 × 10-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.6 × 10-4). We replicated the association of FAw6/FA with LTL in an independent dataset of 7845 individuals (P = 1.9 × 10-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.

Associations of carotid intima media thickness with gene expression in whole blood and genetically predicted gene expression across 48 tissues.

Carotid intima media thickness (cIMT) is a biomarker of subclinical atherosclerosis and a predictor of future cardiovascular events. Identifying associations between gene expression levels and cIMT may provide insight to atherosclerosis etiology. Here, we use two approaches to identify associations between mRNA levels and cIMT: differential gene expression analysis in whole blood and S-PrediXcan. We used microarrays to measure genome-wide whole blood mRNA levels of 5647 European individuals from four studies. We examined the association of mRNA levels with cIMT adjusted for various potential confounders. Significant associations were tested for replication in three studies totaling 3943 participants. Next, we applied S-PrediXcan to summary statistics from a cIMT genome-wide association study (GWAS) of 71 128 individuals to estimate the association between genetically determined mRNA levels and cIMT and replicated these analyses using S-PrediXcan on an independent GWAS on cIMT that included 22 179 individuals from the UK Biobank. mRNA levels of TNFAIP3, CEBPD and METRNL were inversely associated with cIMT, but these associations were not significant in the replication analysis. S-PrediXcan identified associations between cIMT and genetically determined mRNA levels for 36 genes, of which six were significant in the replication analysis, including TLN2, which had not been previously reported for cIMT. There was weak correlation between our results using differential gene expression analysis and S-PrediXcan. Differential expression analysis and S-PrediXcan represent complementary approaches for the discovery of associations between phenotypes and gene expression. Using these approaches, we prioritize TNFAIP3, CEBPD, METRNL and TLN2 as new candidate genes whose differential expression might modulate cIMT.

Long-term glucocorticoids in relation to the metabolic syndrome and cardiovascular disease: A systematic review and meta-analysis.

The striking link of Cushing's syndrome with the metabolic syndrome (MetS) and cardiovascular disease (CVD) suggests that long-term exposure to extremely high cortisol levels catalyzes cardiometabolic deterioration. However, it remained unclear whether the findings from the extreme glucocorticoid overabundance observed in Cushing's syndrome could be translated into more subtle variations in long-term glucocorticoid levels among the general population, for example, due to chronic stress. Here, we performed a systematic review (PROSPERO: CRD42023425541) of evidence regarding the role of subtle variations in long-term biological stress, measured as levels of scalp hair cortisol (HairF) and cortisone (HairE), in the context of MetS and CVD in adults. We also performed a meta-analysis on the cross-sectional difference in HairF levels between individuals with versus without CVD. Seven studies were included regarding MetS, sixteen regarding CVD, and one regarding both. Most articles indicated a strong, consistent cross-sectional association of higher HairF and HairE levels with CVD, which was confirmed by our meta-analysis for HairF (eight studies, SMD = 0.48, 95% confidence intervals [CIs]: 0.16-0.79, p = 0.0095). Moreover, these relationships appear largely independent of standard risk factors. Age seems relevant as the effect seems stronger in younger individuals. Results regarding the associations of HairF and HairE with MetS were inconsistent. Altogether, long-term biological stress, measured as HairF and HairE, is associated with the presence of CVD, and less consistently with MetS. Prospective studies need to evaluate the directionality of this relationship and determine whether HairF and HairE can be used in addition to standard risk factors in predicting future cardiometabolic deterioration.

Features of immunometabolic depression as predictors of antidepressant treatment outcomes: pooled analysis of four clinical trials.

BACKGROUND: Profiling patients on a proposed 'immunometabolic depression' (IMD) dimension, described as a cluster of atypical depressive symptoms related to energy regulation and immunometabolic dysregulations, may optimise personalised treatment. AIMS: To test the hypothesis that baseline IMD features predict poorer treatment outcomes with antidepressants. METHOD: Data on 2551 individuals with depression across the iSPOT-D (n = 967), CO-MED (n = 665), GENDEP (n = 773) and EMBARC (n = 146) clinical trials were used. Predictors included baseline severity of atypical energy-related symptoms (AES), body mass index (BMI) and C-reactive protein levels (CRP, three trials only) separately and aggregated into an IMD index. Mixed models on the primary outcome (change in depressive symptom severity) and logistic regressions on secondary outcomes (response and remission) were conducted for the individual trial data-sets and pooled using random-effects meta-analyses. RESULTS: Although AES severity and BMI did not predict changes in depressive symptom severity, higher baseline CRP predicted smaller reductions in depressive symptoms (n = 376, ßpooled = 0.06, P = 0.049, 95% CI 0.0001-0.12, I2 = 3.61%); this was also found for an IMD index combining these features (n = 372, ßpooled = 0.12, s.e. = 0.12, P = 0.031, 95% CI 0.01-0.22, I2= 23.91%), with a higher - but still small - effect size compared with CRP. Confining analyses to selective serotonin reuptake inhibitor users indicated larger effects of CRP (ßpooled = 0.16) and the IMD index (ßpooled = 0.20). Baseline IMD features, both separately and combined, did not predict response or remission. CONCLUSIONS: Depressive symptoms of people with more IMD features improved less when treated with antidepressants. However, clinical relevance is limited owing to small effect sizes in inconsistent associations. Whether these patients would benefit more from treatments targeting immunometabolic pathways remains to be investigated.

Day-to-day affect fluctuations in adults with childhood trauma history: a two-week ecological momentary assessment study.

BACKGROUND: Childhood trauma (CT) may increase vulnerability to psychopathology through affective dysregulation (greater variability, autocorrelation, and instability of emotional symptoms). However, CT associations with dynamic affect fluctuations while considering differences in mean affect levels across CT status have been understudied. METHODS: 346 adults (age = 49.25 ± 12.55, 67.0% female) from the Netherlands Study of Depression and Anxiety participated in ecological momentary assessment. Positive and negative affect (PA, NA) were measured five times per day for two weeks by electronic diaries. Retrospectively-reported CT included emotional neglect and emotional/physical/sexual abuse. Linear regressions determined associations between CT and affect fluctuations, controlling for age, sex, education, and mean affect levels. RESULTS: Compared to those without CT, individuals with CT reported significantly lower mean PA levels (Cohen's d = -0.620) and higher mean NA levels (d = 0.556) throughout the two weeks. CT was linked to significantly greater PA variability (d = 0.336), NA variability (d = 0.353), and NA autocorrelation (d = 0.308), with strongest effects for individuals reporting higher CT scores. However, these effects were entirely explained by differences in mean affect levels between the CT groups. Findings suggested consistency of results in adults with and without lifetime depressive/anxiety disorders and across CT types, with sexual abuse showing the smallest effects. CONCLUSIONS: Individuals with CT show greater affective dysregulation during the two-week monitoring of emotional symptoms, likely due to their consistently lower PA and higher NA levels. It is essential to consider mean affect level when interpreting the impact of CT on affect dynamics.

Metabolic syndrome after childhood trauma: a 9-year longitudinal analysis.

BACKGROUND: Childhood trauma (CT) has been cross-sectionally associated with metabolic syndrome (MetS), a group of biological risk factors for cardiometabolic disease. Longitudinal studies, while rare, would clarify the development of cardiometabolic dysregulations over time. Therefore, we longitudinally investigated the association of CT with the 9-year course of MetS components. METHODS: Participants (N = 2958) from the Netherlands Study of Depression and Anxiety were assessed four times across 9 years. The CT interview retrospectively assessed childhood emotional neglect and physical, emotional, and sexual abuse. Metabolic outcomes encompassed continuous MetS components (waist circumference, triglycerides, high-density lipoprotein [HDL] cholesterol, blood pressure [BP], and glucose) and count of clinically elevated MetS components. Mixed-effects models estimated sociodemographic- and lifestyle-adjusted longitudinal associations of CT with metabolic outcomes over time. Time interactions evaluated change in these associations. RESULTS: CT was reported by 49% of participants. CT was consistently associated with increased waist (b = 0.32, s.e. = 0.10, p = 0.001), glucose (b = 0.02, s.e. = 0.01, p < 0.001), and count of MetS components (b = 0.04, s.e. = 0.01, p < 0.001); and decreased HDL cholesterol (b = -0.01, s.e.<0.01, p = .020) and systolic BP (b = -0.33, s.e. = 0.13, p = 0.010). These associations were mainly driven by severe CT and unaffected by lifestyle. Only systolic BP showed a CT-by-time interaction, where CT was associated with lower systolic BP initially and with higher systolic BP at the last follow-up. CONCLUSIONS: Over time, adults with CT have overall persistent poorer metabolic outcomes than their non-maltreated peers. Individuals with CT have an increased risk for cardiometabolic disease and may benefit from monitoring and early interventions targeting metabolism.

Identification of latent classes in mood and anxiety disorders and their transitions over time: a follow-up study in the adult general population.

BACKGROUND: Mood and anxiety disorders are heterogeneous conditions with variable course. Knowledge on latent classes and transitions between these classes over time based on longitudinal disorder status information provides insight into clustering of meaningful groups with different disease prognosis. METHODS: Data of all four waves of the Netherlands Mental Health Survey and Incidence Study-2 were used, a representative population-based study of adults (mean duration between two successive waves = 3 years; N at T0 = 6646; T1 = 5303; T2 = 4618; T3 = 4007; this results in a total number of data points: 20 574). Presence of eight mood and anxiety DSM-IV disorders was assessed with the Composite International Diagnostic Interview. Latent class analysis and latent Markov modelling were used. RESULTS: The best fitting model identified four classes: a healthy class (prevalence: 94.1%), depressed-worried class (3.6%; moderate-to-high proportions of mood disorders and generalized anxiety disorder (GAD)), fear class (1.8%; moderate-to-high proportions of panic and phobia disorders) and high comorbidity class (0.6%). In longitudinal analyses over a three-year period, the minority of those in the depressed-worried and high comorbidity class persisted in their class over time (36.5% and 38.4%, respectively), whereas the majority in the fear class did (67.3%). Suggestive of recovery is switching to the healthy class, this was 39.7% in the depressed-worried class, 12.5% in the fear class and 7.0% in the high comorbidity class. CONCLUSIONS: People with panic or phobia disorders have a considerably more persistent and chronic disease course than those with depressive disorders including GAD. Consequently, they could especially benefit from longer-term monitoring and disease management.

Lagged effects of childhood depressive symptoms on adult epigenetic aging.

BACKGROUND: Cross-sectional studies have identified health risks associated with epigenetic aging. However, it is unclear whether these risks make epigenetic clocks 'tick faster' (i.e. accelerate biological aging). The current study examines concurrent and lagged within-person changes of a variety of health risks associated with epigenetic aging. METHODS: Individuals from the Great Smoky Mountains Study were followed from age 9 to 35 years. DNA methylation profiles were assessed from blood, at multiple timepoints (i.e. waves) for each individual. Health risks were psychiatric, lifestyle, and adversity factors. Concurrent (N = 539 individuals; 1029 assessments) and lagged (N = 380 individuals; 760 assessments) analyses were used to determine the link between health risks and epigenetic aging. RESULTS: Concurrent models showed that BMI (r = 0.15, PFDR < 0.01) was significantly correlated to epigenetic aging at the subject-level but not wave-level. Lagged models demonstrated that depressive symptoms (b = 1.67 months per symptom, PFDR = 0.02) in adolescence accelerated epigenetic aging in adulthood, also when models were fully adjusted for BMI, smoking, and cannabis and alcohol use. CONCLUSIONS: Within-persons, changes in health risks were unaccompanied by concurrent changes in epigenetic aging, suggesting that it is unlikely for risks to immediately 'accelerate' epigenetic aging. However, time lagged analyses indicated that depressive symptoms in childhood/adolescence predicted epigenetic aging in adulthood. Together, findings suggest that age-related biological embedding of depressive symptoms is not instant but provides prognostic opportunities. Repeated measurements and longer follow-up times are needed to examine stable and dynamic contributions of childhood experiences to epigenetic aging across the lifespan.

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.

The mediating role of health behaviors in the association between depression, anxiety and cancer incidence: an individual participant data meta-analysis.

BACKGROUND: Although behavioral mechanisms in the association among depression, anxiety, and cancer are plausible, few studies have empirically studied mediation by health behaviors. We aimed to examine the mediating role of several health behaviors in the associations among depression, anxiety, and the incidence of various cancer types (overall, breast, prostate, lung, colorectal, smoking-related, and alcohol-related cancers). METHODS: Two-stage individual participant data meta-analyses were performed based on 18 cohorts within the Psychosocial Factors and Cancer Incidence consortium that had a measure of depression or anxiety (N = 319 613, cancer incidence = 25 803). Health behaviors included smoking, physical inactivity, alcohol use, body mass index (BMI), sedentary behavior, and sleep duration and quality. In stage one, path-specific regression estimates were obtained in each cohort. In stage two, cohort-specific estimates were pooled using random-effects multivariate meta-analysis, and natural indirect effects (i.e. mediating effects) were calculated as hazard ratios (HRs). RESULTS: Smoking (HRs range 1.04-1.10) and physical inactivity (HRs range 1.01-1.02) significantly mediated the associations among depression, anxiety, and lung cancer. Smoking was also a mediator for smoking-related cancers (HRs range 1.03-1.06). There was mediation by health behaviors, especially smoking, physical inactivity, alcohol use, and a higher BMI, in the associations among depression, anxiety, and overall cancer or other types of cancer, but effects were small (HRs generally below 1.01). CONCLUSIONS: Smoking constitutes a mediating pathway linking depression and anxiety to lung cancer and smoking-related cancers. Our findings underline the importance of smoking cessation interventions for persons with depression or anxiety.

Running therapy or antidepressants as treatments for immunometabolic depression in patients with depressive and anxiety Disorders: A secondary analysis of the MOTAR study.

BACKGROUND: Exercise promotes immunometabolic health and is increasingly recognized as an effective depression treatment. Exercise may be beneficial for patients with immunometabolic depression (IMD), who experience inflammatory and metabolic dysregulations and may respond less to antidepressants. This secondary analysis of the MOTAR study compared the effects of running therapy and antidepressants on IMD features among patients with depression and/or anxiety disorder. We additionally assessed whether baseline IMD moderated intervention effects on depression. METHODS: Participants received 16 weeks of group-based running therapy (N = 96) or escitalopram/sertraline (N = 45) in a partially randomized patient preference design. IMD features included atypical, energy-related symptom (AES) severity, inflammation index (CRP, INF-γ, IL-6, TNF-α), metabolic syndrome index, three metabolite principle components (PC) (derived from 73 metabolites) and a composite IMD index. RESULTS: Interventions differed in changes in the metabolic syndrome index (d = 0.59, p = 0.026) and IMD index (d = 0.85, p < 0.001). While running therapy decreased both outcomes, the antidepressant group showed an increased IMD index. Although groups did not differ statistically significant in changes in AES severity, inflammation index, and metabolite PC1, results indicated a consistent trend towards greater improvement with running therapy across these outcomes as well (d = 0.38 to 0.52). Baseline IMD did not moderate intervention effects on depression outcomes. CONCLUSIONS: This study suggests that exercise more effectively targets the IMD dimension than antidepressants. Patients with IMD did not benefit more from running therapy than antidepressants in terms of reductions in depression. Exercise should be considered an alternative or complementary treatment to particularly reduce IMD features in depressed patients. TRIAL REGISTRATION: Trialregister.nl Number of identification: NTR3460.

Methylome-wide association study of anxiety disorders.

Anxiety Disorders (ANX) such as panic disorder, generalized anxiety disorder, and phobias, are highly prevalent conditions that are moderately heritable. Evidence suggests that DNA methylation may play a role, as it is involved in critical adaptations to changing environments. Applying an enrichment-based sequencing approach covering nearly 28 million autosomal CpG sites, we conducted a methylome-wide association study (MWAS) of lifetime ANX in 1132 participants (618 cases/514 controls) from the Netherlands Study of Depression and Anxiety. Using epigenomic deconvolution, we performed MWAS for the main cell types in blood: granulocytes, T-cells, B-cells and monocytes. Cell-type specific analyses identified 280 and 82 methylome-wide significant associations (q-value < 0.1) in monocytes and granulocytes, respectively. Our top finding in monocytes was located in ZNF823 on chromosome 19 (p = 1.38 × 10-10) previously associated with schizophrenia. We observed significant overlap (p < 1 × 10-06) with the same direction of effect in monocytes (210 sites), T-cells (135 sites), and B-cells (727 sites) between this Discovery MWAS signal and a comparable replication dataset from the Great Smoky Mountains Study (N = 433). Overlapping Discovery-Replication MWAS signal was enriched for findings from published GWAS of ANX, major depression, and post-traumatic stress disorder. In monocytes, two specific sites in the FZR1 gene showed significant replication after Bonferroni correction with an additional 15 nominally replicated sites in monocytes and 4 in T-cells. FZR1 regulates neurogenesis in the hippocampus, and its knockout leads to impairments in associative fear memory and long-term potentiation in mice. In the largest and most extensive methylome-wide study of ANX, we identified replicable methylation sites located in genes of potential relevance for brain mechanisms of psychiatric conditions.

The Mood and Resilience in Offspring (MARIO) project: a longitudinal cohort study among offspring of parents with and without a mood disorder.

BACKGROUND: One of the most robust risk factors for developing a mood disorder is having a parent with a mood disorder. Unfortunately, mechanisms explaining the transmission of mood disorders from one generation to the next remain largely elusive. Since timely intervention is associated with a better outcome and prognosis, early detection of intergenerational transmission of mood disorders is of paramount importance. Here, we describe the design of the Mood and Resilience in Offspring (MARIO) cohort study in which we investigate: 1. differences in clinical, biological and environmental (e.g., psychosocial factors, substance use or stressful life events) risk and resilience factors in children of parents with and without mood disorders, and 2. mechanisms of intergenerational transmission of mood disorders via clinical, biological and environmental risk and resilience factors. METHODS: MARIO is an observational, longitudinal cohort study that aims to include 450 offspring of parents with a mood disorder (uni- or bipolar mood disorders) and 100-150 offspring of parents without a mood disorder aged 10-25 years. Power analyses indicate that this sample size is sufficient to detect small to medium sized effects. Offspring are recruited via existing Dutch studies involving patients with a mood disorder and healthy controls, for which detailed clinical, environmental and biological data of the index-parent (i.e., the initially identified parent with or without a mood disorder) is available. Over a period of three years, four assessments will take place, in which extensive clinical, biological and environmental data and data on risk and resilience are collected through e.g., blood sampling, face-to-face interviews, online questionnaires, actigraphy and Experience Sampling Method assessment. For co-parents, information on demographics, mental disorder status and a DNA-sample are collected. DISCUSSION: The MARIO cohort study is a large longitudinal cohort study among offspring of parents with and without mood disorders. A unique aspect is the collection of granular data on clinical, biological and environmental risk and resilience factors in offspring, in addition to available parental data on many similar factors. We aim to investigate the mechanisms underlying intergenerational transmission of mood disorders, which will ultimately lead to better outcomes for offspring at high familial risk.

Association of CYP2D6 and CYP2C19 metabolizer status with switching and discontinuing antidepressant drugs: an exploratory study.

BACKGROUND: Tailoring antidepressant drugs (AD) to patients' genetic drug-metabolism profile is promising. However, literature regarding associations of ADs' treatment effect and/or side effects with drug metabolizing genes CYP2D6 and CYP2C19 has yielded inconsistent results. Therefore, our aim was to longitudinally investigate associations between CYP2D6 (poor, intermediate, and normal) and CYP2C19 (poor, intermediate, normal, and ultrarapid) metabolizer-status, and switching/discontinuing of ADs. Next, we investigated whether the number of perceived side effects differed between metabolizer statuses. METHODS: Data came from the multi-site naturalistic longitudinal cohort Netherlands Study of Depression and Anxiety (NESDA). We selected depression- and/or anxiety patients, who used AD at some point in the course of the 9 years follow-up period (n = 928). Medication use was followed to assess patterns of AD switching/discontinuation over time. CYP2D6 and CYP2C19 alleles were derived using genome-wide data of the NESDA samples and haplotype data from the PharmGKB database. Logistic regression analyses were conducted to investigate the association of metabolizer status with switching/discontinuing ADs. Mann-Whitney U-tests were conducted to compare the number of patient-perceived side effects between metabolizer statuses. RESULTS: No significant associations were observed of CYP metabolizer status with switching/discontinuing ADs, nor with the number of perceived side effects. CONCLUSIONS: We found no evidence for associations between CYP metabolizer statuses and switching/discontinuing AD, nor with side effects of ADs, suggesting that metabolizer status only plays a limited role in switching/discontinuing ADs. Additional studies with larger numbers of PM and UM patients are needed to further determine the potential added value of pharmacogenetics to guide pharmacotherapy.

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.

Educational attainment does not influence brain aging.

Education has been related to various advantageous lifetime outcomes. Here, using longitudinal structural MRI data (4,422 observations), we tested the influential hypothesis that higher education translates into slower rates of brain aging. Cross-sectionally, education was modestly associated with regional cortical volume. However, despite marked mean atrophy in the cortex and hippocampus, education did not influence rates of change. The results were replicated across two independent samples. Our findings challenge the view that higher education slows brain aging.