Genome-wide characterization of circulating metabolic biomarkers.

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.

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.

Normal range CAG repeat size variations in the HTT gene are associated with an adverse lipoprotein profile partially mediated by body mass index.

Tandem cytosine-adenine-guanine (CAG) repeat sizes of 36 or more in the huntingtin gene (HTT) cause Huntington's disease (HD). Apart from neuropsychiatric complications, the disease is also accompanied by metabolic dysregulation and weight loss, which contribute to a progressive functional decline. Recent studies also reported an association between repeats below the pathogenic threshold (<36) for HD and body mass index (BMI), suggesting that HTT repeat sizes in the non-pathogenic range are associated with metabolic dysregulation. In this study, we hypothesized that HTT repeat sizes < 36 are associated with metabolite levels, possibly mediated through reduced BMI. We pooled data from three European cohorts (n = 10 228) with genotyped HTT CAG repeat size and metabolomic measurements. All 145 metabolites were measured on the same targeted platform in all studies. Multilevel mixed-effects analysis using the CAG repeat size in HTT identified 67 repeat size metabolite associations. Overall, the metabolomic profile associated with larger CAG repeat sizes in HTT were unfavorable-similar to those of higher risk of coronary artery disease and type 2 diabetes-and included elevated levels of amino acids, fatty acids, low-density lipoprotein (LDL)-, very low-density lipoprotein- and intermediate density lipoprotein (IDL)-related metabolites while with decreased levels of very large high-density lipoprotein (HDL)-related metabolites. Furthermore, the associations of 50 metabolites, in particular, specific very large HDL-related metabolites, were mediated by lower BMI. However, no mediation effect was found for 17 metabolites related to LDL and IDL. In conclusion, our findings indicate that large non-pathogenic CAG repeat sizes in HTT are associated with an unfavorable metabolomic profile despite their association with a lower BMI.

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.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

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.

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

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

Genome-wide association study of circulating interleukin 6 levels identifies novel loci.

Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67 428 (ndiscovery = 52 654 and nreplication = 14 774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (Pcombined = 1.8 × 10-11), HLA-DRB1/DRB5 rs660895 on Chr6p21 (Pcombined = 1.5 × 10-10) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (Pcombined = 1.2 × 10-122). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology.

Genome-wide Association Analysis in Humans Links Nucleotide Metabolism to Leukocyte Telomere Length.

Leukocyte telomere length (LTL) is a heritable biomarker of genomic aging. In this study, we perform a genome-wide meta-analysis of LTL by pooling densely genotyped and imputed association results across large-scale European-descent studies including up to 78,592 individuals. We identify 49 genomic regions at a false dicovery rate (FDR) < 0.05 threshold and prioritize genes at 31, with five highlighting nucleotide metabolism as an important regulator of LTL. We report six genome-wide significant loci in or near SENP7, MOB1B, CARMIL1, PRRC2A, TERF2, and RFWD3, and our results support recently identified PARP1, POT1, ATM, and MPHOSPH6 loci. Phenome-wide analyses in >350,000 UK Biobank participants suggest that genetically shorter telomere length increases the risk of hypothyroidism and decreases the risk of thyroid cancer, lymphoma, and a range of proliferative conditions. Our results replicate previously reported associations with increased risk of coronary artery disease and lower risk for multiple cancer types. Our findings substantially expand current knowledge on genes that regulate LTL and their impact on human health and disease.

Depression, cardiometabolic disease, and their co-occurrence after childhood maltreatment: an individual participant data meta-analysis including over 200,000 participants.

BACKGROUND: Childhood maltreatment is associated with depression and cardiometabolic disease in adulthood. However, the relationships with these two diseases have so far only been evaluated in different samples and with different methodology. Thus, it remains unknown how the effect sizes magnitudes for depression and cardiometabolic disease compare with each other and whether childhood maltreatment is especially associated with the co-occurrence ("comorbidity") of depression and cardiometabolic disease. This pooled analysis examined the association of childhood maltreatment with depression, cardiometabolic disease, and their comorbidity in adulthood. METHODS: We carried out an individual participant data meta-analysis on 13 international observational studies (N = 217,929). Childhood maltreatment comprised self-reports of physical, emotional, and/or sexual abuse before 18 years. Presence of depression was established with clinical interviews or validated symptom scales and presence of cardiometabolic disease with self-reported diagnoses. In included studies, binomial and multinomial logistic regressions estimated sociodemographic-adjusted associations of childhood maltreatment with depression, cardiometabolic disease, and their comorbidity. We then additionally adjusted these associations for lifestyle factors (smoking status, alcohol consumption, and physical activity). Finally, random-effects models were used to pool these estimates across studies and examined differences in associations across sex and maltreatment types. RESULTS: Childhood maltreatment was associated with progressively higher odds of cardiometabolic disease without depression (OR [95% CI] = 1.27 [1.18; 1.37]), depression without cardiometabolic disease (OR [95% CI] = 2.68 [2.39; 3.00]), and comorbidity between both conditions (OR [95% CI] = 3.04 [2.51; 3.68]) in adulthood. Post hoc analyses showed that the association with comorbidity was stronger than with either disease alone, and the association with depression was stronger than with cardiometabolic disease. Associations remained significant after additionally adjusting for lifestyle factors, and were present in both males and females, and for all maltreatment types. CONCLUSIONS: This meta-analysis revealed that adults with a history of childhood maltreatment suffer more often from depression and cardiometabolic disease than their non-exposed peers. These adults are also three times more likely to have comorbid depression and cardiometabolic disease. Childhood maltreatment may therefore be a clinically relevant indicator connecting poor mental and somatic health. Future research should investigate the potential benefits of early intervention in individuals with a history of maltreatment on their distal mental and somatic health (PROSPERO CRD42021239288).

Dissection of depression heterogeneity using proteomic clusters.

BACKGROUND: The search for relevant biomarkers of major depressive disorder (MDD) is challenged by heterogeneity; biological alterations may vary in patients expressing different symptom profiles. Moreover, most research considers a limited number of biomarkers, which may not be adequate for tagging complex network-level mechanisms. Here we studied clusters of proteins and examined their relation with MDD and individual depressive symptoms. METHODS: The sample consisted of 1621 subjects from the Netherlands Study of Depression and Anxiety (NESDA). MDD diagnoses were based on DSM-IV criteria and the Inventory of Depressive Symptomatology questionnaire measured endorsement of 30 symptoms. Serum protein levels were detected using a multi-analyte platform (171 analytes, immunoassay, Myriad RBM DiscoveryMAP 250+). Proteomic clusters were computed using weighted correlation network analysis (WGCNA). RESULTS: Six proteomic clusters were identified, of which one was nominally significantly associated with current MDD (p = 9.62E-03, Bonferroni adj. p = 0.057). This cluster contained 21 analytes and was enriched with pathways involved in inflammation and metabolism [including C-reactive protein (CRP), leptin and insulin]. At the individual symptom level, this proteomic cluster was associated with ten symptoms, among which were five atypical, energy-related symptoms. After correcting for several health and lifestyle covariates, hypersomnia, increased appetite, panic and weight gain remained significantly associated with the cluster. CONCLUSIONS: Our findings support the idea that alterations in a network of proteins involved in inflammatory and metabolic processes are present in MDD, but these alterations map predominantly to clinical symptoms reflecting an imbalance between energy intake and expenditure.

Metabolomics dissection of depression heterogeneity and related cardiometabolic risk.

BACKGROUND: A recent hypothesis postulates the existence of an 'immune-metabolic depression' (IMD) dimension characterized by metabolic dysregulations. Combining data on metabolomics and depressive symptoms, we aimed to identify depressions associated with an increased risk of adverse metabolic alterations. METHOD: Clustering data were from 1094 individuals with major depressive disorder in the last 6 months and measures of 149 metabolites from a 1H-NMR platform and 30 depressive symptoms (IDS-SR30). Canonical correlation analyses (CCA) were used to identify main independent metabolite-symptom axes of variance. Then, for the replication, we examined the association of the identified dimensions with metabolites from the same platform and cardiometabolic diseases in an independent population-based cohort (n = 6572). RESULTS: CCA identified an overall depression dimension and a dimension resembling IMD, in which symptoms such as sleeping too much, increased appetite, and low energy level had higher relative loading. In the independent sample, the overall depression dimension was associated with lower cardiometabolic risk, such as (i.e. per s.d.) HOMA-1B -0.06 (95% CI -0.09 - -0.04), and visceral adipose tissue -0.10 cm2 (95% CI -0.14 - -0.07). In contrast, the IMD dimension was associated with well-known cardiometabolic diseases such as higher visceral adipose tissue 0.08 cm2 (95% CI 0.04-0.12), HOMA-1B 0.06 (95% CI 0.04-0.09), and lower HDL-cholesterol levels -0.03 mmol/L (95% CI -0.05 - -0.01). CONCLUSIONS: Combining metabolomics and clinical symptoms we identified a replicable depression dimension associated with adverse metabolic alterations, in line with the IMD hypothesis. Patients with IMD may be at higher cardiometabolic risk and may benefit from specific treatment targeting underlying metabolic dysregulations.

A network analysis of depressive symptoms and metabolomics.

BACKGROUND: Depression is associated with metabolic alterations including lipid dysregulation, whereby associations may vary across individual symptoms. Evaluating these associations using a network perspective yields a more complete insight than single outcome-single predictor models. METHODS: We used data from the Netherlands Study of Depression and Anxiety (N = 2498) and leveraged networks capturing associations between 30 depressive symptoms (Inventory of Depressive Symptomatology) and 46 metabolites. Analyses involved 4 steps: creating a network with Mixed Graphical Models; calculating centrality measures; bootstrapping for stability testing; validating central, stable associations by extra covariate-adjustment; and validation using another data wave collected 6 years later. RESULTS: The network yielded 28 symptom-metabolite associations. There were 15 highly-central variables (8 symptoms, 7 metabolites), and 3 stable links involving the symptoms Low energy (fatigue), and Hypersomnia. Specifically, fatigue showed consistent associations with higher mean diameter for VLDL particles and lower estimated degree of (fatty acid) unsaturation. These remained present after adjustment for lifestyle and health-related factors and using another data wave. CONCLUSIONS: The somatic symptoms Fatigue and Hypersomnia and cholesterol and fatty acid measures showed central, stable, and consistent relationships in our network. The present analyses showed how metabolic alterations are more consistently linked to specific symptom profiles.

The association between adiposity and atypical energy-related symptoms of depression: A role for metabolic dysregulations.

BACKGROUND: Adiposity has been shown to be linked with atypical energy-related symptoms (AES) of depression. We used genomics to separate the effect of adiposity from that of metabolic dysregulations to examine whether the link between obesity and AES is dependent on the presence of metabolic dysregulations. METHOD: Data were from NEO (n = 5734 individuals) and NESDA (n = 2238 individuals) cohorts, in which the Inventory of Depressive Symptomatology (IDS-SR30) was assessed. AES profile was based on four symptoms: increased appetite, increased weight, low energy level, and leaden paralysis. We estimated associations between AES and two genetic risk scores (GRS) indexing increasing total body fat with (metabolically unhealthy adiposity, GRS-MUA) and without (metabolically healthy adiposity, GRS-MHA) metabolic dysregulations. RESULTS: We validated that both GRS-MUA and GRS-MHA were associated with higher total body fat in NEO study, but divergently associated with biomarkers of metabolic health (e.g., fasting glucose and HDL-cholesterol) in both cohorts. In the pooled results, per standard deviation, GRS-MUA was specifically associated with a higher AES score (ß = 0.03, 95%CI: 0.01; 0.05), while there was no association between GRS-MHA and AES (ß = -0.01, 95%CI: -0.03; 0.01). CONCLUSION: These results suggest that the established link between adiposity and AES profile emerges in the presence of metabolic dysregulations, which may represent the connecting substrate between the two conditions.

Inflammation and depression in young people: a systematic review and proposed inflammatory pathways.

Depression onset peaks during adolescence and young adulthood. Current treatments are only moderately effective, driving the search for novel pathophysiological mechanisms underlying youth depression. Inflammatory dysregulation has been shown in adults with depression, however, less is known about inflammation in youth depression. This systematic review identified 109 studies examining the association between inflammation and youth depression and showed subtle evidence for inflammatory dysregulation in youth depression. Longitudinal studies support the bidirectional association between inflammation and depression in youth. We hypothesise multiple inflammatory pathways contributing to depression. More research is needed on anti-inflammatory treatments, potentially tailored to individual symptom profiles.

The phenotypic and genotypic association of grip strength with frailty, physical performance and functional limitations over time in older adults.

OBJECTIVES: To replicate the phenotypic associations of grip strength with frailty, physical performance and functional limitations in older adults for longer follow-up periods and to examine whether these associations are due to shared genetic factors. METHODS: In total 2,262 participants 55 years and older with follow-up data up to 23 years (Nobservations = 8,262) from the Longitudinal Aging Study Amsterdam were included. Weighted polygenic risk scores for grip strength (PRS-GS) were built using the genome-wide meta-analysis results from UK Biobank as reference. Grip strength was measured two times on each hand using a dynamometer. Frailty index (FI) and frailty phenotype were operationalised following standard procedures. Performance tests included a timed walk test, a repeated chair stands test and put on-take off cardigan test. Functional limitations were assessed using a questionnaire with six items. RESULTS: Higher grip strength was phenotypically associated with lower FI (b = -0.013, 95% CI (-0.016, -0.009)), better physical performance (b = 0.040, 95% CI (0.026, 0.054)) and less functional limitations (OR = 0.965, 95% CI (0.954, 0.977)) over time for follow-up periods up to 23 years. However, PRS-GS was not associated with any of the traits. CONCLUSION: The phenotypic associations between grip strength, frailty, physical performance and functional limitations were replicated for follow-up periods up to 23 years. However, the associations between the traits could not be explained by shared genetics potentially indicating a more relevant involvement of non-genetic factors.

Weighing psychosocial factors in relatives for the risk of psychopathology: a study of patients with depressive and anxiety disorders and their siblings.

PURPOSE: Siblings of probands with depressive and anxiety disorders are at increased risk for psychopathology, but little is known about how risk factors operate within families to increase psychopathology for siblings. We examined the additional impact of psychosocial risk factors in probands-on top of or in combination with those in siblings-on depressive/anxious psychopathology in siblings. METHODS: The sample included 636 participants (Mage = 49.7; 62.4% female) from 256 families, each including a proband with lifetime depressive and/or anxiety disorders and their sibling(s) (N = 380 proband-sibling pairs). Sixteen psychosocial risk factors were tested. In siblings, depressive and anxiety disorders were determined with standardized psychiatric interviews; symptom severity was measured using self-report questionnaires. Analyses were performed with mixed-effects models accounting for familial structure. RESULTS: In siblings, various psychosocial risk factors (female gender, low income, childhood trauma, poor parental bonding, being single, smoking, hazardous alcohol use) were associated with higher symptomatology and likelihood of disorder. The presence of the same risk factor in probands was independently associated (low income, being single) with higher symptomatology in siblings or moderated (low education, childhood trauma, hazardous alcohol use)-by reducing its strength-the association between the risk factor and symptomatology in siblings. There was no additional impact of risk factors in probands on likelihood of disorder in siblings. CONCLUSION: Our findings demonstrate the importance of weighing psychosocial risk factors within a family context, as it may provide relevant information on the risk of affective psychopathology for individuals.

Familial risk for depressive and anxiety disorders: associations with genetic, clinical, and psychosocial vulnerabilities.

BACKGROUND: In research and clinical practice, familial risk for depression and anxiety is often constructed as a simple Yes/No dichotomous family history (FH) indicator. However, this measure may not fully capture the liability to these conditions. This study investigated whether a continuous familial loading score (FLS), incorporating family- and disorder-specific characteristics (e.g. family size, prevalence of depression/anxiety), (i) is associated with a polygenic risk score (PRS) for major depression and with clinical/psychosocial vulnerabilities and (ii) still captures variation in clinical/psychosocial vulnerabilities after information on FH has been taken into account. METHODS: Data came from 1425 participants with lifetime depression and/or anxiety from the Netherlands Study of Depression and Anxiety. The Family Tree Inventory was used to determine FLS/FH indicators for depression and/or anxiety. RESULTS: Persons with higher FLS had higher PRS for major depression, more severe depression and anxiety symptoms, higher disease burden, younger age of onset, and more neuroticism, rumination, and childhood trauma. Among these variables, FH was not associated with PRS, severity of symptoms, and neuroticism. After regression out the effect of FH from the FLS, the resulting residualized measure of FLS was still associated with severity of symptoms of depression and anxiety, rumination, and childhood trauma. CONCLUSIONS: Familial risk for depression and anxiety deserves clinical attention due to its associated genetic vulnerability and more unfavorable disease profile, and seems to be better captured by a continuous score that incorporates family- and disorder-specific characteristics than by a dichotomous FH measure.

Metabolomic and inflammatory signatures of symptom dimensions in major depression.

BACKGROUND: Major depressive disorder (MDD) is a highly heterogenous disease, both in terms of clinical profiles and pathobiological alterations. Recently, immunometabolic dysregulations were shown to be correlated with atypical, energy-related symptoms but less so with the Melancholic or Anxious distress symptom dimensions of depression in The Netherlands Study of Depression and Anxiety (NESDA) study. In this study, we aimed to replicate these immunometabolic associations and to characterize the metabolomic correlates of each of the three MDD dimensions. METHODS: Using three clinical rating scales, Melancholic, and Anxious distress, and Immunometabolic (IMD) dimensions were characterized in 158 patients who participated in the Predictors of Remission to Individual and Combined Treatments (PReDICT) study and from whom plasma and serum samples were available. The NESDA-defined inflammatory index, a composite measure of interleukin-6 and C-reactive protein, was measured from pre-treatment plasma samples and a metabolomic profile was defined using serum samples analyzed on three metabolomics platforms targeting fatty acids and complex lipids, amino acids, acylcarnitines, and gut microbiome-derived metabolites among other metabolites of central metabolism. RESULTS: The IMD clinical dimension and the inflammatory index were positively correlated (r = 0.19, p = 0.019) after controlling for age, sex, and body mass index, whereas the Melancholic and Anxious distress dimensions were not, replicating the previous NESDA findings. The three symptom dimensions had distinct metabolomic signatures using both univariate and set enrichment statistics. IMD severity correlated mainly with gut-derived metabolites and a few acylcarnitines and long chain saturated free fatty acids. Melancholia severity was significantly correlated with several phosphatidylcholines, primarily the ether-linked variety, lysophosphatidylcholines, as well as several amino acids. Anxious distress severity correlated with several medium and long chain free fatty acids, both saturated and polyunsaturated ones, sphingomyelins, as well as several amino acids and bile acids. CONCLUSION: The IMD dimension of depression appears reliably associated with markers of inflammation. Metabolomics provides powerful tools to inform about depression heterogeneity and molecular mechanisms related to clinical dimensions in MDD, which include a link to gut microbiome and lipids implicated in membrane structure and function.