Search

Abstract Existing metabolomic clocks exhibit deficiencies in capturing the heterogeneous aging rates among individuals with the same chronological age. Yet, the modifiable and non‐modifiable factors in metabolomic aging have not been systematically studied. Here, a new aging measure—MetaboAgeMort—is developed using metabolomic profiles from 239,291 UK Biobank participants for 10‐year all‐cause mortality prediction. The MetaboAgeMort showed significant associations with all‐cause mortality, cause‐specific mortality, and diverse incident diseases. Adding MetaboAgeMort to a conventional risk factors model improved the predictive ability of 10‐year mortality. A total of 99 modifiable factors across seven categories are identified for MetaboAgeMort. Among these, 16 factors representing pulmonary function, body composition, socioeconomic status, dietary quality, smoking status, alcohol intake, and disease status showed quantitatively stronger associations. The genetic analyses revealed 99 genomic risk loci and 271 genes associated with MetaboAgeMort. The tissue‐enrichment analysis showed significant enrichment in liver. While the external validation of the MetaboAgeMort is required, this study illuminates heterogeneous metabolomic aging across the same age, providing avenues for identifying high‐risk individuals, developing anti‐aging therapies, and personalizing interventions, thus promoting healthy aging and longevity.

Abstract Background The gut-derived metabolite Trimethylamine N-oxide (TMAO) and its precursors - betaine, carnitine, choline, and deoxycarnitine – have been associated with an increased risk of cardiovascular disease, but their relation to cognition, neuroimaging markers, and dementia remains uncertain. Methods In the population-based Rotterdam Study, we used multivariable regression models to study the associations between plasma TMAO, its precursors, and cognition in 3,143 participants. Subsequently, we examined their link to structural brain MRI markers in 2,047 participants, with a partial validation in the Leiden Longevity Study (n = 318). Among 2,517 participants, we assessed the risk of incident dementia using multivariable Cox proportional hazard models. Following this, we stratified the longitudinal associations by medication use and sex, after which we conducted a sensitivity analysis for individuals with impaired renal function. Results Overall, plasma TMAO was not associated with cognition, neuroimaging markers or incident dementia. Instead, higher plasma choline was significantly associated with poor cognition (adjusted mean difference: -0.170 [95% confidence interval (CI) -0.297;-0.043]), brain atrophy and more markers of cerebral small vessel disease, such as white matter hyperintensity volume (0.237 [95% CI: 0.076;0.397]). By contrast, higher carnitine concurred with lower white matter hyperintensity volume (-0.177 [95% CI: -0.343;-0.010]). Only among individuals with impaired renal function, TMAO appeared to increase risk of dementia (hazard ratio (HR): 1.73 [95% CI: 1.16;2.60]). No notable differences were observed in stratified analyses. Conclusions Plasma choline, as opposed to TMAO, was found to be associated with cognitive decline, brain atrophy, and markers of cerebral small vessel disease. These findings illustrate the complexity of relationships between TMAO and its precursors, and emphasize the need for concurrent study to elucidate gut-brain mechanisms.

Background: Aging triggers intricate physiological changes, particularly in whole-body fat-free mass (FFM) and handgrip strength, affecting overall health and independence. Despite existing research, the broader significance of how muscle health is affected by the intricate interplay of lifestyle factors simultaneously during aging needs more exploration. This study aims to examine how nutrition, physical activity, and sleep impact on FFM and handgrip strength in middle-aged men and women, facilitating future personalized recommendations for preserving muscle health. Methods: The cross-sectional analysis of the UK Biobank involved 45,984 individuals (54 % women) aged 40–70 years with a complete dataset. Multiple linear regression explored determinants of FFM and handgrip strength, considering traditional, socio-demographics, medication use and smoking as covariates, with sex and age (younger and older than 55 years) stratifications. Results: In older men and women, higher physical activity beneficially affect both FFM (respectively Β = 3.36 × 10−3, p-value = 1.66 × 10−3; Β = 2.52 × 10−3, p-value = 3.57 × 10−4) and handgrip strength (Β = 6.05 × 10−3, p-value = 7.99 × 10−5, Β = 8.98 × 10−3, p-value = 2.95 × 10−15). Similar results were found in fiber intake for FFM in older men and women (respectively B = 3.00 × 10−2, p-value = 2.76 × 10−5; B = 2.68 × 10−2, p-value = 1.78 × 10−9) and handgrip strength (Β = 3.27 × 10−2, p-value = 1.40 × 10−3; Β = 3.12 × 10−2, p-value = 1.34 × 10−5). Other lifestyle factors influence FFM and handgrip strength differently. Key determinants influencing handgrip strength included higher protein intake, lower water intake, higher alcohol intake, and extended sleep duration whereas mainly higher water intake is associated with higher FFM. Conclusions: In both men and women, the main factors associated with FFM and handgrip strength are physical activity and fiber intake, which may underlie a connection between gut and muscle health. Given the observed complexity of muscle health in the age and sex strata, further longitudinal research is needed to provide personalized lifestyle recommendations.

Summary: Recent studies have shown that elevated concentrations of unconjugated bilirubin (UCB) may be a protective host factor against the development of noncommunicable diseases (NCDs), whereas low levels of UCB are associated with the opposite effect. The results of this European study, in which 2,489 samples were tested for their UCB concentration using high-performance liquid chromatography (HPLC) and additional data from the MARK-AGE database were used for analysis, provide further evidence that elevated UCB concentrations are linked to a lower risk of developing NCDs and may act as a predictive marker of biological aging as individuals with elevated UCB concentrations showed favorable outcomes in metabolic health and oxidative-stress-related biomarkers. These findings underline the significance of studying individuals with moderate hyperbilirubinemia and investigate UCB routinely, also in the setting of aging, since this condition affects millions of people worldwide but has been underrepresented in clinical research and practice until now.

Abstract Background There is divergence in the rate at which people age. The concept of biological age is postulated to capture this variability, and hence to better represent an individual’s true global physiological state than chronological age. Biological age predictors are often generated based on cross-sectional data, using biochemical or molecular markers as predictor variables. It is assumed that the difference between chronological and predicted biological age is informative of one’s chronological age-independent aging divergence ∆. Methods We investigated the statistical assumptions underlying the most popular cross-sectional biological age predictors, based on multiple linear regression, the Klemera-Doubal method or principal component analysis. We used synthetic and real data to illustrate the consequences if this assumption does not hold. Results The most popular cross-sectional biological age predictors all use the same strong underlying assumption, namely that a candidate marker of aging’s association with chronological age is directly informative of its association with the aging rate ∆. We called this the identical-association assumption and proved that it is untestable in a cross-sectional setting. If this assumption does not hold, weights assigned to candidate markers of aging are uninformative, and no more signal may be captured than if markers would have been assigned weights at random. Conclusions Cross-sectional methods for predicting biological age commonly use the untestable identical-association assumption, which previous literature in the field had never explicitly acknowledged. These methods have inherent limitations and may provide uninformative results, highlighting the importance of researchers exercising caution in the development and interpretation of cross-sectional biological age predictors.

Abstract Background Dietary intake of n-3 polyunsaturated fatty acids (PUFA) may have a protective effect on the development of cardiovascular diseases, diabetes, depression and cancer, while a high intake of n-6 PUFA was often reported to be associated with inflammation-related traits. The effect of PUFAs on health outcomes might be mediated by DNA methylation (DNAm). The aim of our study is to identify the impact of PUFA intake on DNAm in the Cooperative Health Research in the Region of Augsburg (KORA) FF4 cohort and the Leiden Longevity Study (LLS). Results DNA methylation levels were measured in whole blood from the population-based KORA FF4 study (N = 1354) and LLS (N = 448), using the Illumina MethylationEPIC BeadChip and Illumina HumanMethylation450 array, respectively. We assessed associations between DNAm and intake of eight and four PUFAs in KORA and LLS, respectively. Where possible, results were meta-analyzed. Below the Bonferroni correction threshold (p

Abstract To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Plant-based diets (PBD) may offer various health benefits and contribute to a sustainable way of life, but, if not planned correctly, may also confer risks, e.g., by focusing on plant foods with low nutrient density, such as foods primarily consisting of refined carbohydrates. A plant-based diet index (PDI) differentiating between a healthful, unhealthful, and overall PBD, offers a promising approach to standardize and compare studies and integrate results. In this review we (1) summarize current evidence on the PDI and disease risk of relevance to public health, (2) discuss the methodology of the PDI and how it can be sensibly applied in further studies and (3) indicate areas with a lack of knowledge, such as vulnerable populations. In summary, our amalgamation shows, that adherence to a healthier plant-based diet is associated with an 8–68% lower risk for metabolic risk factors, diabetes, and cardiovascular disease, while adherence to an unhealthier plant-based diet is associated with a 10–63% higher risk. Although differences in calculation methods and underlying diet patterns between populations should be accounted for, the PDI can be a useful tool to assess adherence to different plant-based diet patterns and their association with health outcomes in cohort studies across cultures.

Abstract Trimethylamine N‐oxide (TMAO) is a circulating microbiome‐derived metabolite implicated in the development of atherosclerosis and cardiovascular disease (CVD). We investigated whether plasma levels of TMAO, its precursors (betaine, carnitine, deoxycarnitine, choline), and TMAO‐to‐precursor ratios are associated with clinical outcomes, including CVD and mortality. This was followed by an in‐depth analysis of their genetic, gut microbial, and dietary determinants. The analyses were conducted in five Dutch prospective cohort studies including 7834 individuals. To further investigate association results, Mendelian Randomization (MR) was also explored. We found only plasma choline levels (hazard ratio [HR] 1.17, [95% CI 1.07; 1.28]) and not TMAO to be associated with CVD risk. Our association analyses uncovered 10 genome‐wide significant loci, including novel genomic regions for betaine (6p21.1, 6q25.3), choline (2q34, 5q31.1), and deoxycarnitine (10q21.2, 11p14.2) comprising several metabolic gene associations, for example, CPS1 or PEMT. Furthermore, our analyses uncovered 68 gut microbiota associations, mainly related to TMAO‐to‐precursors ratios and the Ruminococcaceae family, and 16 associations of food groups and metabolites including fish‐TMAO, meat‐carnitine, and plant‐based food‐betaine associations. No significant association was identified by the MR approach. Our analyses provide novel insights into the TMAO pathway, its determinants, and pathophysiological impact on the general population.

Abstract Background B vitamins such as folate (B9), B6, and B12 are key in one carbon metabolism, which generates methyl donors for DNA methylation. Several studies have linked differential methylation to self-reported intakes of folate and B12, but these estimates can be imprecise, while metabolomic biomarkers can offer an objective assessment of dietary intakes. We explored blood metabolomic biomarkers of folate and vitamins B6 and B12, to carry out epigenome-wide analyses across up to three European cohorts. Associations between self-reported habitual daily B vitamin intakes and 756 metabolites (Metabolon Inc.) were assessed in serum samples from 1064 UK participants from the TwinsUK cohort. The identified B vitamin metabolomic biomarkers were then used in epigenome-wide association tests with fasting blood DNA methylation levels at 430,768 sites from the Infinium HumanMethylation450 BeadChip in blood samples from 2182 European participants from the TwinsUK and KORA cohorts. Candidate signals were explored for metabolite associations with gene expression levels in a subset of the TwinsUK sample (n = 297). Metabolomic biomarker epigenetic associations were also compared with epigenetic associations of self-reported habitual B vitamin intakes in samples from 2294 European participants. Results Eighteen metabolites were associated with B vitamin intakes after correction for multiple testing (Bonferroni-adj. p

Abstract Globally, the lifespan of populations increases but the healthspan is lagging behind. Previous research showed that survival into extreme ages (longevity) clusters in families as illustrated by the increasing lifespan of study participants with each additional long-lived family member. Here we investigate whether the healthspan in such families follows a similar quantitative pattern using three-generational data from two databases, LLS (Netherlands), and SEDD (Sweden). We study healthspan in 2143 families containing index persons with 26 follow-up years and two ancestral generations, comprising 17,539 persons. Our results provide strong evidence that an increasing number of long-lived ancestors associates with up to a decade of healthspan extension. Further evidence indicates that members of long-lived families have a delayed onset of medication use, multimorbidity and, in mid-life, healthier metabolomic profiles than their partners. We conclude that both lifespan and healthspan are quantitatively linked to ancestral longevity, making family data invaluable to identify protective mechanisms of multimorbidity.

Summary: Background: Smoking impacts DNA methylation, but data are lacking on smoking-related differential methylation by sex or dietary intake, recent smoking cessation (

Introduction: Endothelial dysfunction has various detrimental effects on the brain that may lead to cognitive decline. However, data on the direct association between endothelial dysfunction and cognitive performance in humans are scarce. The aim of the present study was to assess the cross-sectional association between plasma biomarkers of endothelial dysfunction and cognitive performance. Methods: Data from 9,344 individuals without dementia from seven Dutch cohorts participating in the Netherlands Consortium of Dementia Cohorts (NCDC) were included: Doetinchem Cohort Study, Rotterdam Study, The Maastricht Study, Leiden Longevity Study, the EMIF-AD 90+ cohort, Longitudinal Aging Study Amsterdam, and SMART-MR Study (mean age range: 57-93 years, 18-58% women). Plasma biomarkers of endothelial dysfunction included intercellular adhesion molecule-1 (sICAM-1), soluble vascular cell adhesion molecule-1 (sVCAM-1) and soluble E-selectin (sE-selectin), which were measured in EDTA plasma samples with commercially available ELISA kits (MSD and R&D Systems Europe), or the Olink Cardiometabolic 96-plex panels. Subsequently, a standardized endothelial dysfunction composite score was calculated. Cognitive performance was measured on six domains: executive function, processing speed, immediate and delayed memory, attention, and language. A two-step individual participant data (IPD) meta-analysis was performed: linear regressions were run in the individual cohorts. Resulting standardized regression coefficients were subsequently pooled using random-effects models. Results: Pooled analyses showed an association between a higher endothelial dysfunction composite score and worse executive function (β = -0.04; 95% CI -0.06, -0.02), slower processing speed (β = -0.03; 95% CI -0.05, -0.01), worse immediate memory (β = -0.03; 95% CI -0.04, -0.01), worse delayed memory (β = -0.03; 95% CI -0.05, -0.00), and lower attention scores (β = -0.02; 95% CI -0.04, -0.00) (Figure 1). These associations were consistent across cohorts and independent of demographics, lifestyle and cardiovascular risk factors. No association was found between endothelial dysfunction and language (β = -0.02; 95% CI -0.05, 0.01). Discussion: Endothelial dysfunction was consistently associated with worse cognitive functioning across all domains, except for language. These findings support the hypothesis that endothelial dysfunction may contribute to cognitive impairment and may identify a new prevention target.Future research is needed to assess any longitudinal association between endothelial dysfunction and cognitive decline.

Circulating lipids can influence immune cell function, which could have implications for inflammatory diseases such as rheumatoid arthritis. Here, the authors use Mendelian randomization to identify genes whose expression is influenced by triglyceride levels in blood, implicating genes involved in lipid metabolism and allergic response.

Abstract Population-scale expression profiling studies can provide valuable insights into biological and disease-underlying mechanisms. The availability of phenotypic traits is essential for studying clinical effects. Therefore, missing, incomplete, or inaccurate phenotypic information can make analyses challenging and prevent RNA-seq or other omics data to be reused. A possible solution are predictors that infer clinical or behavioral phenotypic traits from molecular data. While such predictors have been developed based on different omics data types and are being applied in various studies, metabolomics-based surrogates are less commonly used than predictors based on DNA methylation profiles.In this study, we inferred 17 traits, including diabetes status and exposure to lipid medication, using previously trained metabolomic predictors. We evaluated whether these metabolomic surrogates can be used as an alternative to reported information for studying the respective phenotypes using expression profiling data of four population cohorts. For the majority of the 17 traits, the metabolomic surrogates performed similarly to the reported phenotypes in terms of effect sizes, number of significant associations, replication rates, and significantly enriched pathways.The application of metabolomics-derived surrogate outcomes opens new possibilities for reuse of multi-omics data sets. In studies where availability of clinical metadata is limited, missing or incomplete information can be complemented by these surrogates, thereby increasing the size of available data sets. Additionally, the availability of such surrogates could be used to correct for potential biological confounding. In the future, it would be interesting to further investigate the use of molecular predictors across different omics types and cohorts.

1H-NMR metabolomics data is increasingly used to track health and disease. Nightingale Health, a major supplier of 1H-NMR metabolomics, has recently updated the quantification strategy to further align with clinical standards. Such updates, however, might influence backward replicability, particularly affecting studies with repeated measures. Using data from BBMRI-NL consortium (~28,000 samples from 28 cohorts), we compared Nightingale data, originally released in 2014 and 2016, with a re-quantified version released in 2020, of which both versions were based on the same NMR spectra. Apart from two discontinued and twenty-three new analytes, we generally observe a high concordance between quantification versions with 73 out of 222 (33%) analytes showing a mean ρ > 0.9 across all cohorts. Conversely, five analytes consistently showed lower Spearman’s correlations (ρ < 0.7) between versions, namely acetoacetate, LDL-L, saturated fatty acids, S-HDL-C, and sphingomyelins. Furthermore, previously trained multi-analyte scores, such as MetaboAge or MetaboHealth, might be particularly sensitive to platform changes. Whereas MetaboHealth replicated well, the MetaboAge score had to be retrained due to use of discontinued analytes. Notably, both scores in the re-quantified data recapitulated mortality associations observed previously. Concluding, we urge caution in utilizing different platform versions to avoid mixing analytes, having different units, or simply being discontinued.

Abstract Background Epigenetic clocks use DNA methylation (DNAm) levels of specific sets of CpG dinucleotides to accurately predict individual chronological age. A popular application of these clocks is to explore whether the deviation of predicted age from chronological age is associated with disease phenotypes, where this deviation is interpreted as a potential biomarker of biological age. This wide application, however, contrasts with the limited insight in the processes that may drive the running of epigenetic clocks. Results We perform a functional genomics analysis on four epigenetic clocks, including Hannum’s blood predictor and Horvath’s multi-tissue predictor, using blood DNA methylome and transcriptome data from 3132 individuals. The four clocks result in similar predictions of individual chronological age, and their constituting CpGs are correlated in DNAm level and are enriched for similar histone modifications and chromatin states. Interestingly, DNAm levels of CpGs from the clocks are commonly associated with gene expression in trans. The gene sets involved are highly overlapping and enriched for T cell processes. Further analysis of the transcriptome and methylome of sorted blood cell types identifies differences in DNAm between naive and activated T and NK cells as a probable contributor to the clocks. Indeed, within the same donor, the four epigenetic clocks predict naive cells to be up to 40 years younger than activated cells. Conclusions The ability of epigenetic clocks to predict chronological age involves their ability to detect changes in proportions of naive and activated immune blood cells, an established feature of immuno-senescence. This finding may contribute to the interpretation of associations between clock-derived measures and age-related health outcomes.

Abstract Background Observational studies suggest interconnections between thyroid status, metabolism, and risk of coronary artery disease (CAD), but causality remains to be proven. The present study aimed to investigate the potential causal relationship between thyroid status and cardiovascular disease and to characterize the metabolomic profile associated with thyroid status. Methods Multi-cohort two-sample Mendelian randomization (MR) was performed utilizing genome-wide significant variants as instruments for standardized thyrotropin (TSH) and free thyroxine (fT4) within the reference range. Associations between TSH and fT4 and metabolic profile were investigated in a two-stage manner: associations between TSH and fT4 and the full panel of 161 metabolomic markers were first assessed hypothesis-free, then directional consistency was assessed through Mendelian randomization, another metabolic profile platform, and in individuals with biochemically defined thyroid dysfunction. Results Circulating TSH was associated with 52/161 metabolomic markers, and fT4 levels were associated with 21/161 metabolomic markers among 9432 euthyroid individuals (median age varied from 23.0 to 75.4 years, 54.5% women). Positive associations between circulating TSH levels and concentrations of very low-density lipoprotein subclasses and components, triglycerides, and triglyceride content of lipoproteins were directionally consistent across the multivariable regression, MR, metabolomic platforms, and for individuals with hypo- and hyperthyroidism. Associations with fT4 levels inversely reflected those observed with TSH. Among 91,810 CAD cases and 656,091 controls of European ancestry, per 1-SD increase of genetically determined TSH concentration risk of CAD increased slightly, but not significantly, with an OR of 1.03 (95% CI 0.99–1.07; p value 0.16), whereas higher genetically determined fT4 levels were not associated with CAD risk (OR 1.00 per SD increase of fT4; 95% CI 0.96–1.04; p value 0.59). Conclusions Lower thyroid status leads to an unfavorable lipid profile and a somewhat increased cardiovascular disease risk.

The mechanisms underlying how monozygotic (or identical) twins arise are yet to be determined. Here, the authors investigate this in an epigenome-wide association study, showing that monozygotic twinning has a characteristic DNA methylation signature in adult somatic tissues.

Abstract Background Aging is a multifactorial process that affects multiple tissues and is characterized by changes in homeostasis over time, leading to increased morbidity. Whole blood gene expression signatures have been associated with aging and have been used to gain information on its biological mechanisms, which are still not fully understood. However, blood is composed of many cell types whose proportions in blood vary with age. As a result, previously observed associations between gene expression levels and aging might be driven by cell type composition rather than intracellular aging mechanisms. To overcome this, previous aging studies already accounted for major cell types, but the possibility that the reported associations are false positives driven by less prevalent cell subtypes remains. Results Here, we compared the regression model from our previous work to an extended model that corrects for 33 additional white blood cell subtypes. Both models were applied to whole blood gene expression data from 3165 individuals belonging to the general population (age range of 18–81 years). We evaluated that the new model is a better fit for the data and it identified fewer genes associated with aging (625, compared to the 2808 of the initial model; P ≤ 2.5⨯10−6). Moreover, 511 genes (~ 18% of the 2808 genes identified by the initial model) were found using both models, indicating that the other previously reported genes could be proxies for less abundant cell types. In particular, functional enrichment of the genes identified by the new model highlighted pathways and GO terms specifically associated with platelet activity. Conclusions We conclude that gene expression analyses in blood strongly benefit from correction for both common and rare blood cell types, and recommend using blood-cell count estimates as standard covariates when studying whole blood gene expression.

Abstract Background Sleep traits are associated with cardiometabolic disease risk, with evidence from Mendelian randomization (MR) suggesting that insomnia symptoms and shorter sleep duration increase coronary artery disease risk. We combined adjusted multivariable regression (AMV) and MR analyses of phenotypes of unfavourable sleep on 113 metabolomic traits to investigate possible biochemical mechanisms linking sleep to cardiovascular disease. Methods We used AMV (N = 17,368) combined with two-sample MR (N = 38,618) to examine effects of self-reported insomnia symptoms, total habitual sleep duration, and chronotype on 113 metabolomic traits. The AMV analyses were conducted on data from 10 cohorts of mostly Europeans, adjusted for age, sex, and body mass index. For the MR analyses, we used summary results from published European-ancestry genome-wide association studies of self-reported sleep traits and of nuclear magnetic resonance (NMR) serum metabolites. We used the inverse-variance weighted (IVW) method and complemented this with sensitivity analyses to assess MR assumptions. Results We found consistent evidence from AMV and MR analyses for associations of usual vs. sometimes/rare/never insomnia symptoms with lower citrate (− 0.08 standard deviation (SD)[95% confidence interval (CI) − 0.12, − 0.03] in AMV and − 0.03SD [− 0.07, − 0.003] in MR), higher glycoprotein acetyls (0.08SD [95% CI 0.03, 0.12] in AMV and 0.06SD [0.03, 0.10) in MR]), lower total very large HDL particles (− 0.04SD [− 0.08, 0.00] in AMV and − 0.05SD [− 0.09, − 0.02] in MR), and lower phospholipids in very large HDL particles (− 0.04SD [− 0.08, 0.002] in AMV and − 0.05SD [− 0.08, − 0.02] in MR). Longer total sleep duration associated with higher creatinine concentrations using both methods (0.02SD per 1 h [0.01, 0.03] in AMV and 0.15SD [0.02, 0.29] in MR) and with isoleucine in MR analyses (0.22SD [0.08, 0.35]). No consistent evidence was observed for effects of chronotype on metabolomic measures. Conclusions Whilst our results suggested that unfavourable sleep traits may not cause widespread metabolic disruption, some notable effects were observed. The evidence for possible effects of insomnia symptoms on glycoprotein acetyls and citrate and longer total sleep duration on creatinine and isoleucine might explain some of the effects, found in MR analyses of these sleep traits on coronary heart disease, which warrant further investigation.

Abstract Background Although human longevity tends to cluster within families, genetic studies on longevity have had limited success in identifying longevity loci. One of the main causes of this limited success is the selection of participants. Studies generally include sporadically long-lived individuals, i.e. individuals with the longevity phenotype but without a genetic predisposition for longevity. The inclusion of these individuals causes phenotype heterogeneity which results in power reduction and bias. A way to avoid sporadically long-lived individuals and reduce sample heterogeneity is to include family history of longevity as selection criterion using a longevity family score. A main challenge when developing family scores are the large differences in family size, because of real differences in sibship sizes or because of missing data. Methods We discussed the statistical properties of two existing longevity family scores: the Family Longevity Selection Score (FLoSS) and the Longevity Relatives Count (LRC) score and we evaluated their performance dealing with differential family size. We proposed a new longevity family score, the mLRC score, an extension of the LRC based on random effects modeling, which is robust for family size and missing values. The performance of the new mLRC as selection tool was evaluated in an intensive simulation study and illustrated in a large real dataset, the Historical Sample of the Netherlands (HSN). Results Empirical scores such as the FLOSS and LRC cannot properly deal with differential family size and missing data. Our simulation study showed that mLRC is not affected by family size and provides more accurate selections of long-lived families. The analysis of 1105 sibships of the Historical Sample of the Netherlands showed that the selection of long-lived individuals based on the mLRC score predicts excess survival in the validation set better than the selection based on the LRC score . Conclusions Model-based score systems such as the mLRC score help to reduce heterogeneity in the selection of long-lived families. The power of future studies into the genetics of longevity can likely be improved and their bias reduced, by selecting long-lived cases using the mLRC.

Abstract 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 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), and NLRC5, 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.

Blood metabolite levels are under the influence of environmental and genetic factors. Here, Hagenbeek et al. perform heritability estimations for metabolite measures and determine the contribution of known metabolite loci to metabolite levels using data from 40 genome-wide association studies.

Variants associated with blood lipid levels may be population-specific. To identify low-frequency variants associated with this phenotype, population-specific reference panels may be used. Here we impute nine large Dutch biobanks (~35,000 samples) with the population-specific reference panel created by the Genome of The Netherlands Project and perform association testing with blood lipid levels. We report the discovery of five novel associations at four loci (P value

Educational attainment is a key behavioural measure in studies of cognitive and physical health, and socioeconomic status. We measured DNA methylation at 410,746 CpGs (N = 4152) and identified 58 CpGs associated with educational attainment at loci characterized by pleiotropic functions shared with neuronal, immune and developmental processes. Associations overlapped with those for smoking behaviour, but remained after accounting for smoking at many CpGs: Effect sizes were on average 28% smaller and genome-wide significant at 11 CpGs after adjusting for smoking and were 62% smaller in never smokers. We examined sources and biological implications of education-related methylation differences, demonstrating correlations with maternal prenatal folate, smoking and air pollution signatures, and associations with gene expression in cis, dynamic methylation in foetal brain, and correlations between blood and brain. Our findings show that the methylome of lower-educated people resembles that of smokers beyond effects of their own smoking behaviour and shows traces of various other exposures.

The presence of circulating microbiome in blood has been reported in both physiological and pathological conditions, although its origins, identities and function remain to be elucidated. This study aimed to investigate the presence of blood microbiome by quantitative real-time PCRs targeting the 16S rRNA gene. To our knowledge, this is the first study in which the circulating microbiome has been analyzed in such a large sample of individuals since the study was carried out on 1285 Randomly recruited Age-Stratified Individuals from the General population (RASIG). The samples came from several different European countries recruited within the EU Project MARK-AGE in which a series of clinical biochemical parameters were determined. The results obtained reveal an association between microbial DNA copy number and geographic origin. By contrast, no gender and age-related difference emerged, thus demonstrating the role of the environment in influencing the above levels independent of age and gender at least until the age of 75. In addition, a significant positive association was found with Free Fatty Acids (FFA) levels, leukocyte count, insulin, and glucose levels. Since these factors play an essential role in both health and disease conditions, their association with the extent of the blood microbiome leads us to consider the blood microbiome as a potential biomarker of human health.

Statins effectively lower LDL cholesterol levels in large studies and the observed interindividual response variability may be partially explained by genetic variation. Here we perform a pharmacogenetic meta-analysis of genome-wide association studies (GWAS) in studies addressing the LDL cholesterol response to statins, including up to 18,596 statin-treated subjects. We validate the most promising signals in a further 22,318 statin recipients and identify two loci, SORT1/CELSR2/PSRC1 and SLCO1B1, not previously identified in GWAS. Moreover, we confirm the previously described associations with APOE and LPA. Our findings advance the understanding of the pharmacogenetic architecture of statin response.

Biomarkers that predict mortality are of interest for clinical as well as research applications. Here, the authors analyze metabolomics data from 44,168 individuals and identify key metabolites independently associated with all-cause mortality risk.

Genome-wide association studies have only revealed a handful of genetic loci for longevity. Here, in a case–control design based on phenotype definitions of individuals surviving at or beyond the age corresponding to the 90th and 99th survival percentile, the authors report two additional loci located in the APOE locus and near GPR78.

Christina M. Lill, who contributed to analysis of data, was inadvertently omitted from the author list in the originally published version of this article. This has now been corrected in both the PDF and HTML versions of the article.

While human lifespan is only moderately heritable, “getting old” runs in families. Here, van den Berg et al. study mortality data from three-generation cohorts to define a threshold for longevity and find that individuals have an increasing survival advantage with each additional relative in the top 10% survivors of their birth cohort.

Biomarkers of low-grade inflammation have been associated with susceptibility to a severe infectious disease course, even when measured prior to disease onset. We investigated whether metabolic biomarkers measured by nuclear magnetic resonance (NMR) spectroscopy could be associated with susceptibility to severe pneumonia (2507 hospitalised or fatal cases) and severe COVID-19 (652 hospitalised cases) in 105,146 generally healthy individuals from UK Biobank, with blood samples collected 2007–2010. The overall signature of metabolic biomarker associations was similar for the risk of severe pneumonia and severe COVID-19. A multi-biomarker score, comprised of 25 proteins, fatty acids, amino acids, and lipids, was associated equally strongly with enhanced susceptibility to severe COVID-19 (odds ratio 2.9 [95%CI 2.1–3.8] for highest vs lowest quintile) and severe pneumonia events occurring 7–11 years after blood sampling (2.6 [1.7–3.9]). However, the risk for severe pneumonia occurring during the first 2 years after blood sampling for people with elevated levels of the multi-biomarker score was over four times higher than for long-term risk (8.0 [4.1–15.6]). If these hypothesis generating findings on increased susceptibility to severe pneumonia during the first few years after blood sampling extend to severe COVID-19, metabolic biomarker profiling could potentially complement existing tools for identifying individuals at high risk. These results provide novel molecular understanding on how metabolic biomarkers reflect the susceptibility to severe COVID-19 and other infections in the general population.

Introduction: High mortality rates of approximately 20% within 1 year after treatment are observed for patients with proximal femoral fractures. This preliminary study explores the prognostic value of a previously constructed mortality risk score based on a set of 14 metabolites for the survival and functional recovery in patients with proximal femoral fractures. Materials and Methods: A prospective observational cohort study was conducted including patients admitted with a proximal femoral fracture. The primary outcome was patient survival, and the recovery of independence in activities of daily living was included as a secondary outcome. The mortality risk score was constructed for each patient and its prognostic value was tested for the whole population. Results: Data was available form 136 patients. The mean age of all patients was 82.1 years, with a median follow-up of 6 months. Within this period, 19.0% of all patients died and 51.1% recovered to their prefracture level of independence. The mortality score was significantly associated with mortality (HR, 2.74; 95% CI, 1.61-4.66; P < 0.001), but showed only a fair prediction accuracy (AUC = 0.68) and a borderline significant comparison of the mortality score tertile groups in survival analyses (P = 0.049). No decisive associations were found in any of the analyses for the functional recovery of patients. Discussion: These findings support the previously determined prognostic value of the mortality risk score. However, the independent prognostic value when adjusted for potential confounding factors is yet to be assessed. Also, a risk score constructed for this specific patient population might achieve higher accuracies for the prediction of survival and functional recovery. Conclusions: A modest prediction accuracy was observed for the mortality risk score in this population. More elaborate studies are needed to validate these findings and develop a tailored model for clinical purposes in this patient population.

Thyroid dysfunction is a common public health problem and associated with cardiovascular co-morbidities. Here, the authors carry out genome-wide meta-analysis for thyroid hormone (TH) levels, hyper- and hypothyroidism and identify SLC17A4 as a TH transporter and AADAT as a TH metabolizing enzyme.

DNA methylation is critically involved in X chromosome inactivation (XCI) and dosage compensation, yet some X-chromosomal genes escape XCI. Here, Lujik et al. identify three autosomal genetic loci that associate with differential DNA methylation near genes that variably escape XCI in females.

Identification of causal drivers behind expression is essential for understanding gene function. Here authors develop a method for the large-scale inference of gene–gene interactions in observational population genomics data and characterize a network of trans-effects for 6600 genes.

Cognitive function is associated with health and important life outcomes. Here, the authors perform a genome-wide association study for general cognitive function in 300,486 individuals and identify genetic loci that implicate neural and cell developmental pathways in this trait.

Abstract Background The well-established association of chronological age with changes in DNA methylation is primarily founded on the analysis of large sets of blood samples, while conclusions regarding tissue-specificity are typically based on small number of samples, tissues and CpGs. Here, we systematically investigate the tissue-specific character of age-related DNA methylation changes at the level of the CpG, functional genomic region and nearest gene in a large dataset. Results We assembled a compendium of public data, encompassing genome-wide DNA methylation data (Illumina 450k array) on 8092 samples from 16 different tissues, including 7 tissues with moderate to high sample numbers (Dataset size range 96–1202, N total = 2858). In the 7 tissues (brain, buccal, liver, kidney, subcutaneous fat, monocytes and T-helper cells), we identified 7850 differentially methylated positions that gained (gain-aDMPs; cut-offs: P bonf ≤ 0.05, effect size ≥ 2%/10 years) and 4,287 that lost DNA methylation with age (loss-aDMPs), 92% of which had not previously been reported for whole blood. The majority of all aDMPs identified occurred in one tissue only (gain-aDMPs: 85.2%; loss-aDMPs: 97.4%), an effect independent of statistical power. This striking tissue-specificity extended to both the functional genomic regions (defined by chromatin state segmentation) and the nearest gene. However, aDMPs did accumulate in regions with the same functional annotation across tissues, namely polycomb-repressed CpG islands for gain-aDMPs and regions marked by active histone modifications for loss-aDMPs. Conclusion Our analysis shows that age-related DNA methylation changes are highly tissue-specific. These results may guide the development of improved tissue-specific markers of chronological and, perhaps, biological age.

Vitamin D deficiency is associated with multiple human pathologic conditions. In a genome-wide association study of 79,366 individuals, Jiang et al. replicate four and identify two new genetic loci for serum levels of 25-hydroxyvitamin D and find evidence for a shared genetic basis with autoimmune diseases.

Variability in human longevity is genetically influenced. Using genetic data of parental lifespan, the authors identify associations at HLA-DQA/DRB1 and LPA and find that genetic variants that increase educational attainment have a positive effect on lifespan whereas increasing BMI negatively affects lifespan.

Animal lifespan is plastic and is regulated by conserved signalling pathways. Here, Tikuet al.show that longevity-enhancing mutations or interventions are associated with reduced nucleolar size in worms, flies, mice and humans, and that nucleolar size can predict life-expectancy in individual worms.

Genome-wide association studies (GWAS) have become a key tool to discover genetic markers for complex traits; however, environmental factors that interact with genes are rarely considered. Here, the authors conduct a GWAS of obesity traits, and find that smoking may alter genetic susceptibilities.

Introduction: Up to 33% and 25% of patients with end-stage hip and knee osteoarthritis (OA) are considered frail by the Groningen Frailty Indicator (GFI). This study aims to assess whether frail patients have lower functional gains after arthroplasty and to assess GFI as a tool to discriminate between good and adverse change score. Materials and Methods: Patients with end-stage hip/knee OA scheduled for arthroplasty were recruited from the Longitudinal Leiden Orthopaedics Outcomes of Osteo-Arthritis Study . Functional outcome was measured as change score on the Hip Osteoarthritis Outcome Score/Knee Osteoarthritis Outcome Score (HOOS/KOOS), by subtracting preoperative score from 1-year postsurgery score and then dichotomized based on a cutoff of 20 points. For each HOOS/KOOS subscale, 3 models were estimated: GFI univariate (model 1), GFI and baseline score (model 2), and baseline score univariate (model 3). A receiver operating characteristic analysis was performed to assess the discriminative ability of each model. Results: Eight hundred five patients with end-stage hip OA (31.4% frail) and 640 patients with end-stage knee OA (25.4% frail) were included. Frail patients were older, had a higher body mass index, had more comorbidities, and lived more often alone. Persons considered frail by GFI had significant lower baseline score; however, except for “function in sports and recreation” and “quality of life,” change scores were similar in frail and nonfrail persons. The discriminatory value of GFI was negligible for all HOOS/KOOS subscales. Baseline score, however, was adequate to discriminate between total knee arthroplasty patients with more or less than twice the minimally clinically important difference on KOOS symptoms subscale (area under the curve = 0.802). Discussion/Conclusion: Although frail patients with OA have lower functioning scores at baseline, the change scores on HOOS/KOOS subscales are similar for both frail and nonfrail patients. Exploring other heath assessements may improve patient-specific outcome prediction.

A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-22613-2

Structural variants (SVs) are prevalent in genomes of the general population. Here, Guryev and The Genome of the Netherlands Consortium describe the reference panel of haplotype-resolved SVs from 769 individuals from 250 Dutch families and show its utility for studying heritable traits.

Differential impact of genetic and environmental influences on DNA methylation may result in sex- and age-related physiological variation and disease susceptibility. By analysing DNA methylome of 2,603 individuals from twin families, here, the authors establish a catalogue of between-individual variation in DNA methylation.