BackgroundThe ability to predict the disease course of individuals with major depressive disorder (MDD) is essential for optimal treatment planning. Here, we used a data-driven machine learning... Show moreBackgroundThe ability to predict the disease course of individuals with major depressive disorder (MDD) is essential for optimal treatment planning. Here, we used a data-driven machine learning approach to assess the predictive value of different sets of biological data (whole-blood proteomics, lipid metabolomics, transcriptomics, genetics), both separately and added to clinical baseline variables, for the longitudinal prediction of 2-year remission status in MDD at the individual-subject level.MethodsPrediction models were trained and cross-validated in a sample of 643 patients with current MDD (2-year remission n = 325) and subsequently tested for performance in 161 individuals with MDD (2-year remission n = 82).ResultsProteomics data showed the best unimodal data predictions (area under the receiver operating characteristic curve = 0.68). Adding proteomic to clinical data at baseline significantly improved 2-year MDD remission predictions (area under the receiver operating characteristic curve = 0.63 vs. 0.78, p = .013), while the addition of other omics data to clinical data did not yield significantly improved model performance. Feature importance and enrichment analysis revealed that proteomic analytes were involved in inflammatory response and lipid metabolism, with fibrinogen levels showing the highest variable importance, followed by symptom severity. Machine learning models outperformed psychiatrists’ ability to predict 2-year remission status (balanced accuracy = 71% vs. 55%).ConclusionsThis study showed the added predictive value of combining proteomic data, but not other omics data, with clinical data for the prediction of 2-year remission status in MDD. Our results reveal a novel multimodal signature of 2-year MDD remission status that shows clinical potential for individual MDD disease course predictions from baseline measurements. Show less
Kluiver, H. de; Jansen, R.; Penninx, B.W.J.H.; Giltay, E.J.; Schoevers, R.A.; Milaneschi, Y. 2023
Depression shows a metabolomic signature overlapping with that of cardiometabolic conditions. Whether this signature is linked to specific depression profiles remains undetermined. Previous... Show moreDepression shows a metabolomic signature overlapping with that of cardiometabolic conditions. Whether this signature is linked to specific depression profiles remains undetermined. Previous research suggested that metabolic alterations cluster more consistently with depressive symptoms of the atypical spectrum related to energy alterations, such as hyperphagia, weight gain, hypersomnia, fatigue and leaden paralysis. We characterized the metabolomic signature of an “atypical/energy-related” symptom (AES) profile and evaluated its specificity and consistency. Fifty-one metabolites measured using the Nightingale platform in 2876 participants from the Netherlands Study of Depression and Anxiety were analyzed. An ‘AES profile’ score was based on five items of the Inventory of Depressive Symptomatology (IDS) questionnaire. The AES profile was significantly associated with 31 metabolites including higher glycoprotein acetyls (β = 0.13, p = 1.35*10-12), isoleucine (β = 0.13, p = 1.45*10-10), very-low-density lipoproteins cholesterol (β = 0.11, p = 6.19*10-9) and saturated fatty acid levels (β = 0.09, p = 3.68*10-10), and lower high-density lipoproteins cholesterol (β = −0.07, p = 1.14*10-4). The metabolites were not significantly associated with a summary score of all other IDS items not included in the AES profile. Twenty-five AES-metabolites associations were internally replicated using data from the same subjects (N = 2015) collected at 6-year follow-up. We identified a specific metabolomic signature—commonly linked to cardiometabolic disorders—associated with a depression profile characterized by atypical, energy-related symptoms. The specific clustering of a metabolomic signature with a clinical profile identifies a more homogenous subgroup of depressed patients at higher cardiometabolic risk, and may represent a valuable target for interventions aiming at reducing depression’s detrimental impact on health. Show less
Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of... Show moreIdentifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European ancestry, we identified 43 genomic loci associated with either number of children ever born (NEB) or childlessness. These loci span diverse aspects of reproductive biology, including puberty timing, age at first birth, sex hormone regulation, endometriosis and age at menopause. Missense variants in ARHGAP27 were associated with higher NEB but shorter reproductive lifespan, suggesting a trade-off at this locus between reproductive ageing and intensity. Other genes implicated by coding variants include PIK3IP1, ZFP82 and LRP4, and our results suggest a new role for the melanocortin 1 receptor (MC1R) in reproductive biology. As NEB is one component of evolutionary fitness, our identified associations indicate loci under present-day natural selection. Integration with data from historical selection scans highlighted an allele in the FADS1/2 gene locus that has been under selection for thousands of years and remains so today. Collectively, our findings demonstrate that a broad range of biological mechanisms contribute to reproductive success.Mathieson et al. carried out a genome-wide association study of reproductive success (number of children born) in humans, revealing the importance of diverse neuro-endocrine and behavioural factors. Show less
Schaick, G. van; Hajjouti, N. el; Nicolardi, S.; Hartog, J. den; Jansen, R.; Hoeven, R. van der; ... ; Dominguez Vega, E. 2022
Xylanases are of great value in various industries, including paper, food, and biorefinery. Due to their biotechnological production, these enzymes can contain a variety of post-translational... Show moreXylanases are of great value in various industries, including paper, food, and biorefinery. Due to their biotechnological production, these enzymes can contain a variety of post-translational modifications, which may have a profound effect on protein function. Understanding the structure-function relationship can guide the development of products with optimal performance. We have developed a workflow for the structural and functional characterization of an endo-1,4-beta-xylanase (ENDO-I) produced by Aspergillus niger with and without applying thermal stress. This workflow relies on orthogonal native separation techniques to resolve proteoforms. Mass spectrometry and activity assays of separated proteoforms permitted the establishment of structure-function relationships. The separation conditions were focus on balancing efficient separation and protein functionality. We employed size exclusion chromatography (SEC) to separate ENDO-I from other co-expressed proteins. Charge variants were investigated with ion exchange chromatography (IEX) and revealed the presence of low abundant glycated variants in the temperature-stressed material. To obtain better insights into the effect on glycation on function, we enriched for these species using boronate affinity chromatography (BAC). The activity measurements showed lower activity of glycated species compared to the non-modified enzyme. Altogether, this workflow allowed in-depth structural and functional characterization of ENDO-I proteoforms. Show less
For a long time it has been thought that habitation and landscape organisation only changed significantly from the Roman Period onwards. However, many developments were already started long before... Show moreFor a long time it has been thought that habitation and landscape organisation only changed significantly from the Roman Period onwards. However, many developments were already started long before Julius Caesar's Roman armies arrived in the southern Netherlands. The Iron Age landscapes were ordered and structured, contrasting with the still open Bronze Age landscapes. Iron Age people inhabited the same places for generations. At the same time they structured their immediate environment and surroundings resulting in a sustainable organisation and arrangement of the landscape.Recent excavations and (micro-)regional archaeological studies into habitation and landscape organisation, among others in the north-eastern region of the province Noord-Brabant, show that relicts from the past strongly dictated the organisation and structuring of later landscapes. The past in the past formed a guideline (dutch: leidraad) for later (Iron Age) inhabitants.The past can also be a guideline for the design, protection and preservation of contemporary landscapes. This aligns with a trend in which archaeologists are explicitly seeking the connection with present society. Therefore this book ends with a plea for a transition of the Dutch archaeological system in which living heritage can also be a guideline for the present. Show less
Trait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and trans... Show moreTrait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and trans-expression quantitative trait locus (eQTL) analyses using blood-derived expression from 31,684 individuals through the eQTLGen Consortium. We detected cis-eQTL for 88% of genes, and these were replicable in numerous tissues. Distal trans-eQTL (detected for 37% of 10,317 trait-associated variants tested) showed lower replication rates, partially due to low replication power and confounding by cell type composition. However, replication analyses in single-cell RNA-seq data prioritized intracellular trans-eQTL. Trans-eQTL exerted their effects via several mechanisms, primarily through regulation by transcription factors. Expression of 13% of the genes correlated with polygenic scores for 1,263 phenotypes, pinpointing potential drivers for those traits. In summary, this work represents a large eQTL resource, and its results serve as a starting point for in-depth interpretation of complex phenotypes.Analyses of expression profiles from whole blood of 31,684 individuals identify cis-expression quantitative trait loci (eQTL) effects for 88% of genes and trans-eQTL effects for 37% of trait-associated variants. Show less
Lagou, V.; Magi, R.; Hottenga, J.J.; Grallert, H.; Perry, J.R.B.; Bouatia-Naji, N.; ... ; Meta-Analyses of Glucose and 2021
Differences between sexes contribute to variation in the levels of fasting glucose and insulin. Epidemiological studies established a higher prevalence of impaired fasting glucose in men and... Show moreDifferences between sexes contribute to variation in the levels of fasting glucose and insulin. Epidemiological studies established a higher prevalence of impaired fasting glucose in men and impaired glucose tolerance in women, however, the genetic component underlying this phenomenon is not established. We assess sex-dimorphic (73,089/50,404 women and 67,506/47,806 men) and sex-combined (151,188/105,056 individuals) fasting glucose/fasting insulin genetic effects via genome-wide association study meta-analyses in individuals of European descent without diabetes. Here we report sex dimorphism in allelic effects on fasting insulin at IRS1 and ZNF12 loci, the latter showing higher RNA expression in whole blood in women compared to men. We also observe sex-homogeneous effects on fasting glucose at seven novel loci. Fasting insulin in women shows stronger genetic correlations than in men with waist-to-hip ratio and anorexia nervosa. Furthermore, waist-to-hip ratio is causally related to insulin resistance in women, but not in men. These results position dissection of metabolic and glycemic health sex dimorphism as a steppingstone for understanding differences in genetic effects between women and men in related phenotypes. Show less
Background DNA methylation is a key epigenetic modification in human development and disease, yet there is limited understanding of its highly coordinated regulation. Here, we identify 818 genes... Show moreBackground DNA methylation is a key epigenetic modification in human development and disease, yet there is limited understanding of its highly coordinated regulation. Here, we identify 818 genes that affect DNA methylation patterns in blood using large-scale population genomics data. Results By employing genetic instruments as causal anchors, we establish directed associations between gene expression and distant DNA methylation levels, while ensuring specificity of the associations by correcting for linkage disequilibrium and pleiotropy among neighboring genes. The identified genes are enriched for transcription factors, of which many consistently increased or decreased DNA methylation levels at multiple CpG sites. In addition, we show that a substantial number of transcription factors affected DNA methylation at their experimentally determined binding sites. We also observe genes encoding proteins with heterogenous functions that have widespread effects on DNA methylation, e.g.,NFKBIE,CDCA7(L), andNLRC5, and for several examples, we suggest plausible mechanisms underlying their effect on DNA methylation. Conclusion We report hundreds of genes that affect DNA methylation and provide key insights in the principles underlying epigenetic regulation. Show less
Objectives: The present study examined associations between immunometabolic characteristics (IMCs) and depressive symptom profiles (DSPs) in probands with lifetime diagnoses of depression and/or... Show moreObjectives: The present study examined associations between immunometabolic characteristics (IMCs) and depressive symptom profiles (DSPs) in probands with lifetime diagnoses of depression and/or anxiety disorders and their siblings. Methods: Data were from the Netherlands Study of Depression and Anxiety, comprising 256 probands with lifetime diagnoses of depression and/or anxiety and their 380 siblings. Measured IMCs included blood pressure, waist circumference, and levels of glucose, triglycerides, HDL cholesterol, CRP, TNF-alpha and IL-6. DSPs included mood, cognitive, somatic and atypical-like profiles. We cross-sectionally examined whether DSPs were associated with IMCs within probands and within siblings, and whether DSPs were associated with IMCs between probands and siblings. Results: Within probands and within siblings, higher BMI and waist circumference were associated with higher somatic and atypical-like profiles. Other IMCs (IL-6, glucose and HDL cholesterol) were significantly related to DSPs either within probands or within siblings. DSPs and IMCs were not associated between probands and siblings. Conclusions: The results suggest that there is a familial component for each trait, but no common familial factors for the association between DSPs and IMCs. Alternative mechanisms, such as direct causal effects or non-shared environmental risk factors, may better fit these results. Show less
Insights into individual differences in gene expression and its heritability (h(2)) can help in understanding pathways from DNA to phenotype. We estimated the heritability of gene expression of 52... Show moreInsights into individual differences in gene expression and its heritability (h(2)) can help in understanding pathways from DNA to phenotype. We estimated the heritability of gene expression of 52,844 genes measured in whole blood in the largest twin RNA-Seq sample to date (1497 individuals including 459 monozygotic twin pairs and 150 dizygotic twin pairs) from classical twin modeling and identity-by-state-based approaches. We estimated for each gene h(total)(2), composed of cis-heritability (h(cis)(2), the variance explained by single nucleotide polymorphisms in the cis-window of the gene), and trans-heritability (h(res)(2), the residual variance explained by all other genome-wide variants). Mean h(total)(2) was 0.26, which was significantly higher than heritability estimates earlier found in a microarray-based study using largely overlapping (>60%) RNA samples (mean h(2) = 0.14, p = 6.15 x 10(-258)). Mean h(cis)(2) was 0.06 and strongly correlated with beta of the top cis expression quantitative loci (eQTL, rho = 0.76, p < 10(-308)) and with estimates from earlier RNA-Seq-based studies. Mean h(res)(2) was 0.20 and correlated with the beta of the corresponding trans-eQTL (rho = 0.04, p < 1.89 x 10(-3)) and was significantly higher for genes involved in cytokine-cytokine interactions (p = 4.22 x 10(-15)), many other immune system pathways, and genes identified in genome-wide association studies for various traits including behavioral disorders and cancer. This study provides a thorough characterization of cis- and trans-h(2) estimates of gene expression, which is of value for interpretation of GWAS and gene expression studies. Show less
Windhorst, R.; Alpaslan, M.; Andrews, S.; Ashcraft, T.; Broadhurst, T.; Coe, D.; ... ; Zitrin, A. 2019
The river area Maaskant and adjacent sand area of Oss, located ‘between’ the current course of the river Meuse and the city Oss, are among the most intensively researched regions in the Netherlands... Show moreThe river area Maaskant and adjacent sand area of Oss, located ‘between’ the current course of the river Meuse and the city Oss, are among the most intensively researched regions in the Netherlands. Extensive archaeological and palynological research provides ample opportunities for an interregional research of the occupation and vegetation history of both areas. This article describes the intertwinement between the Holocene river area and the adjacent Pleistocene sandy soils, to eventually get a first insight of the relation(s) between the inhabitants of both regions in late prehistoric and Early Roman period (3000 BC – 250 AD). Show less
Kaal, S.E.J.; Husson, O.; Dartel, F. van; Hermans, K.; Jansen, R.; Manten-Horst, E.; ... ; Graaf, W.T.A. van der 2018