Association analysis between an epigenetic risk score and blood pressure.

Background: Epigenome-wide association studies have revealed multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases. Results: We generated an alcohol consumption epigenetic risk score (ERS) based on previously reported 144 alcohol-associated CpGs and examined the association of the ERS with systolic blood pressure (SBP), diastolic blood pressure (DBP), and hypertension (HTN) in 3,898 Framingham Heart Study (FHS) participants. We found an association of alcohol intake with the ERS in the meta-analysis with 0.09 units higher ERS per drink consumed per day (p < 0.0001). Cross-sectional analyses in FHS revealed that a one-unit increment of the ERS was associated with 1.93 mm Hg higher SBP (p = 4.64E-07), 0.68 mm Hg higher DBP (p = 0.006), and an odds ratio of 1.78 for HTN (p < 2E-16). Meta-analysis of the cross-sectional association of the ERS with BP traits in eight independent external cohorts (n = 11,544) showed similar relationships with blood pressure levels, i.e., a one-unit increase in ERS was associated with 0.74 (p = 0.002) and 0.50 (p = 0.0006) mm Hg higher SBP and DBP, but could not confirm the association with hypertension. Longitudinal analyses in FHS (n = 3,260) and five independent external cohorts (n = 4,021) showed that the baseline ERS was not associated with a change in blood pressure over time or with incident HTN. Conclusions: Our findings provide proof-of-concept that utilizing an ERS is a useful approach to capture the recent health consequences of lifestyle behaviors such as alcohol consumption.

Longitudinal change in mitochondrial heteroplasmy exhibits positive selection for deleterious variants.

A common feature of human aging is the acquisition of somatic mutations, and mitochondria are particularly prone to mutation due to their inefficient DNA repair and close proximity to reactive oxygen species, leading to a state of mitochondrial DNA heteroplasmy1,2. Cross-sectional studies have demonstrated that detection of heteroplasmy increases with participant age3, a phenomenon that has been attributed to genetic drift4-7. In this first large-scale longitudinal study, we measured heteroplasmy in two prospective cohorts (combined n=1405) at two timepoints (mean time between visits, 8.6 years), demonstrating that deleterious heteroplasmies were more likely to increase in variant allele fraction (VAF). We further demonstrated that increase in VAF was associated with increased risk of overall mortality. These results challenge the claim that somatic mtDNA mutations arise mainly due to genetic drift, instead demonstrating positive selection for predicted deleterious mutations at the cellular level, despite an negative impact on overall mortality.

Weight-Bearing Physical Activity, Lower-Limb Muscle Mass, and Risk of Knee Osteoarthritis.

Importance: It has been demonstrated that total physical activity is not associated with risk of osteoarthritis. However, the association of different types of physical activity with incident knee osteoarthritis remains unclear. Objective: To determine whether weight-bearing recreational physical activities are associated with increased risk of incident knee osteoarthritis. Design, Setting, and Participants: This prospective cohort study used data from the Rotterdam Study (1996 to 2009), including participants with knee x-ray measurements at baseline and follow-up examinations. Participants with knee osteoarthritis at baseline were excluded. Residents aged 45 years and older of the Ommoord district in the city of Rotterdam in The Netherlands were invited to join the Rotterdam Study (78% response rate). Analysis was conducted in June 2023. Exposure: Total, weight-bearing, and non-weight-bearing recreational physical activities collected by questionnaires at baseline. Main Outcomes and Measures: Incident radiographic knee osteoarthritis measured by knee x-ray was the primary outcome, and incident symptomatic knee osteoarthritis defined by x-ray and knee pain questionnaire was the secondary outcome. The association of different types of recreational physical activity with radiographic knee osteoarthritis was examined using logistic regression within generalized estimating equation framework after adjusting for potential confounders. A prespecified stratification analysis was planned on the basis of lower-limb muscle mass index (LMI) tertiles, measured by dual-energy x-ray absorptiometry. Results: A total of 5003 individuals (2804 women [56.0%]; mean [SD] age, 64.5 [7.9] years) were included. The knee osteoarthritis incident rate was 8.4% (793 of 9483 knees) for a mean (SD) follow-up time of 6.33 (2.46) years. Higher weight-bearing activity was associated with increased odds of incident knee osteoarthritis (odds ratio [OR], 1.22; 95% CI, 1.10-1.35; P < .001), but non-weight-bearing activity was not (OR, 1.04; 95% CI, 0.95-1.15; P = .37). In the analysis stratified by LMI tertiles, the association of weight-bearing activity with incident osteoarthritis was found only among 431 patients in the lowest LMI tertile (OR, 1.53; 95% CI, 1.15-2.04; P = .003), but not among patients in the middle or high LMI tertile. Conclusions and Relevance: The findings of this study suggest that weight-bearing activity is associated with incident knee osteoarthritis in people with low levels of lower-limb muscle mass, which might be a promising avenue for tailored advice for physical activity.

Association analysis between an epigenetic alcohol risk score and blood pressure.

Background: Epigenome-wide association studies have revealed multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases. Results: We generated an alcohol consumption epigenetic risk score (ERS) based on previously reported 144 alcohol-associated CpGs and examined the association of the ERS with systolic blood pressure (SBP), diastolic blood pressure (DBP), and hypertension (HTN) in 3,898 Framingham Heart Study (FHS) participants. We found an association of alcohol intake with the ERS in the meta-analysis with 0.09 units higher ERS per drink consumed per day (p < 0.0001). Cross-sectional analyses in FHS revealed that a one-unit increment of the ERS was associated with 1.93 mm Hg higher SBP (p = 4.64E-07), 0.68 mm Hg higher DBP (p = 0.006), and an odds ratio of 1.78 for HTN (p < 2E-16). Meta-analysis of the cross-sectional association of the ERS with BP traits in eight independent external cohorts (n = 11,544) showed similar relationships with blood pressure levels, i.e., a one-unit increase in ERS was associated with 0.74 (p = 0.002) and 0.50 (p = 0.0006) mm Hg higher SBP and DBP, but could not confirm the association with hypertension. Longitudinal analyses in FHS (n = 3,260) and five independent external cohorts (n = 4,021) showed that the baseline ERS was not associated with a change in blood pressure over time or with incident HTN. Conclusions: Our findings provide proof-of-concept that utilizing an ERS is a useful approach to capture the recent health consequences of lifestyle behaviors such as alcohol consumption.

Advanced Glycation End Products (AGEs) in Diet and Skin in Relation to Stool Microbiota: The Rotterdam Study.

BACKGROUND: Advanced glycation end products (AGEs) are involved in age-related diseases, but the interaction of gut microbiota with dietary AGEs (dAGEs) and tissue AGEs in the population is unknown. OBJECTIVE: Our objective was to investigate the association of dietary and tissue AGEs with gut microbiota in the population-based Rotterdam Study, using skin AGEs as a marker for tissue accumulation and stool microbiota as a surrogate for gut microbiota. DESIGN: Dietary intake of three AGEs (dAGEs), namely carboxymethyl-lysine (CML), N-(5-hydro-5-methyl-4-imidazolon-2-yl)-ornithine (MGH1), and carboxyethyl-lysine (CEL), was quantified at baseline from food frequency questionnaires. Following up after a median of 5.7 years, skin AGEs were measured using skin autofluorescence (SAF), and stool microbiota samples were sequenced (16S rRNA) to measure microbial composition (including alpha-diversity, beta-dissimilarity, and taxonomic abundances) as well as predict microbial metabolic pathways. Associations of both dAGEs and SAF with microbial measures were investigated using multiple linear regression models in 1052 and 718 participants, respectively. RESULTS: dAGEs and SAF were not associated with either the alpha-diversity or beta-dissimilarity of the stool microbiota. After multiple-testing correction, dAGEs were not associated with any of the 188 genera tested, but were nominally inversely associated with the abundance of Barnesiella, Colidextribacter, Oscillospiraceae UCG-005, and Terrisporobacter, in addition to being positively associated with Coprococcus, Dorea, and Blautia. A higher abundance of Lactobacillus was associated with a higher SAF, along with several nominally significantly associated genera. dAGEs and SAF were nominally associated with several microbial pathways, but none were statistically significant after multiple-testing correction. CONCLUSIONS: Our findings did not solidify a link between habitual dAGEs, skin AGEs, and overall stool microbiota composition. Nominally significant associations with several genera and functional pathways suggested a potential interaction between gut microbiota and AGE metabolism, but validation is required. Future studies are warranted, to investigate whether gut microbiota modifies the potential impact of dAGEs on health.

Epigenetic and Metabolomic Biomarkers for Biological Age: A Comparative Analysis of Mortality and Frailty Risk.

Biological age captures a person's age-related risk of unfavorable outcomes using biophysiological information. Multivariate biological age measures include frailty scores and molecular biomarkers. These measures are often studied in isolation, but here we present a large-scale study comparing them. In 2 prospective cohorts (n = 3 222), we compared epigenetic (DNAm Horvath, DNAm Hannum, DNAm Lin, DNAm epiTOC, DNAm PhenoAge, DNAm DunedinPoAm, DNAm GrimAge, and DNAm Zhang) and metabolomic-based (MetaboAge and MetaboHealth) biomarkers in reflection of biological age, as represented by 5 frailty measures and overall mortality. Biomarkers trained on outcomes with biophysiological and/or mortality information outperformed age-trained biomarkers in frailty reflection and mortality prediction. DNAm GrimAge and MetaboHealth, trained on mortality, showed the strongest association with these outcomes. The associations of DNAm GrimAge and MetaboHealth with frailty and mortality were independent of each other and of the frailty score mimicking clinical geriatric assessment. Epigenetic, metabolomic, and clinical biological age markers seem to capture different aspects of aging. These findings suggest that mortality-trained molecular markers may provide novel phenotype reflecting biological age and strengthen current clinical geriatric health and well-being assessment.

Higher testosterone is associated with open-angle glaucoma in women: a genetic predisposition?

BACKGROUND: Testosterone may be a possible modifiable risk factor for open-angle glaucoma (OAG) and intraocular pressure (IOP), but evidence has been scarce and conflicting. In this study we evaluated the association of testosterone and its genetic predisposition with incident (i) OAG, IOP, retinal nerve fiber layer (RNFL), and ganglion cell-inner plexiform layer (GCL +). METHODS: Participants aged 45-100 years were derived from the prospective, population-based Rotterdam Study. Ophthalmic examinations and serum testosterone measurements (including bioavailable and free testosterone) were performed from 1991 onwards. Follow-up took place every 4-5 years. A total of 187 out of 7898 participants were diagnosed with incident (i) OAG during follow-up. Genotyping was performed in 165 glaucoma cases and 6708 controls. We calculated sex-specific weighted genetic risk scores (GRS) for total and bioavailable testosterone. Associations with iOAG were analyzed using multivariable logistic regression. Associations with IOP, RNFL, and GCL + were analyzed with multivariable linear regression. Analyses were stratified on sex and adjusted for at least age, body mass index, and follow-up duration. RESULTS: In men, testosterone was not associated with iOAG. However, the GRS for higher total testosterone was associated with an increased iOAG risk (odds ratio [OR] with 95% confidence interval [95% CI]: 2.48 [1.18; 5.22], per unit). In women, higher values of bioavailable testosterone (2.05 [1.00; 4.18] per nmol/L) and free testosterone (1.79 [1.00; 3.20] per ng/dL) were significantly associated with increased risk of iOAG. Moreover, the GRS for higher bioavailable testosterone was associated with an increased iOAG risk (2.48 [1.09; 5.65], per unit). Higher bioavailable and free testosterone were adversely associated with IOP (0.58 [0.05; 1.10] per nmol/L and 0.47 [0.04; 0.90] per ng/dL). Higher total testosterone was inversely associated with peripapillary RNFL and GCL + (Beta [95% CI]: - 3.54 [- 7.02; - 0.06] per nmol/L and - 2.18 [- 4.11; - 0.25] per nmol/L, respectively). CONCLUSIONS: In women, higher testosterone levels increased the risk of iOAG. Both IOP-dependent and IOP-independent mechanisms may underlie this association. Managing testosterone levels may be particularly relevant for the prevention of neurodegeneration in the eye. Future research should confirm these findings.

Dietary and supplemental intake of vitamins C and E is associated with altered DNA methylation in an epigenome-wide association study meta-analysis.

BACKGROUND: Dietary intake of antioxidants such as vitamins C and E protect against oxidative stress, and may also be associated with altered DNA methylation patterns. METHODS: We meta-analysed epigenome-wide association study (EWAS) results from 11,866 participants across eight population-based cohorts to evaluate the association between self-reported dietary and supplemental intake of vitamins C and E with DNA methylation. EWAS were adjusted for age, sex, BMI, caloric intake, blood cell type proportion, smoking status, alcohol consumption, and technical covariates. Significant results of the meta-analysis were subsequently evaluated in gene set enrichment analysis (GSEA) and expression quantitative trait methylation (eQTM) analysis. RESULTS: In meta-analysis, methylation at 4,656 CpG sites was significantly associated with vitamin C intake at FDR ≤ 0.05. The most significant CpG sites associated with vitamin C (at FDR ≤ 0.01) were enriched for pathways associated with systems development and cell signalling in GSEA, and were associated with downstream expression of genes enriched in the immune response in eQTM analysis. Furthermore, methylation at 160 CpG sites was significantly associated with vitamin E intake at FDR ≤ 0.05, but GSEA and eQTM analysis of the top most significant CpG sites associated with vitamin E did not identify significant enrichment of any biological pathways investigated. CONCLUSIONS: We identified significant associations of many CpG sites with vitamin C and E intake, and our results suggest that vitamin C intake may be associated with systems development and the immune response.

A functional genomics approach reveals suggestive quantitative trait loci associated with combined TLR4 and BCP crystal-induced inflammation and osteoarthritis.

OBJECTIVE: Basic calcium phosphate (BCP) crystals can activate the NLRP3 inflammasome and are potentially involved in the pathogenesis of osteoarthritis (OA). In order to elucidate relevant inflammatory mechanisms in OA, we used a functional genomics approach to assess genetic variation influencing BCP crystal-induced cytokine production. METHOD: Peripheral blood mononuclear cells (PBMCs) were isolated from healthy volunteers who were previously genotyped and stimulated with BCP crystals and/or lipopolysaccharide (LPS) after which cytokines release was assessed. Cytokine quantitative trait locus (cQTL) mapping was performed. For in vitro validation of the cQTL located in anoctamin 3 (ANO3), PBMCs were incubated with Tamoxifen and Benzbromarone prior to stimulation. Additionally, we performed co-localisation analysis of our top cQTLs with the most recent OA meta-analysis of genome-wide association studies (GWAS). RESULTS: We observed that BCP crystals and LPS synergistically induce IL-1ß in human PBMCs. cQTL analysis revealed several suggestive loci influencing cytokine release upon stimulation, among which are quantitative trait locus annotated to ANO3 and GLIS3. As functional validation, anoctamin inhibitors reduced IL-1ß release in PBMCs after stimulation. Co-localisation analysis showed that the GLIS3 locus was shared between LPS/BCP crystal-induced IL-1ß and genetic association with Knee OA. CONCLUSIONS: We identified and functionally validated a new locus, ANO3, associated with LPS/BCP crystal-induced inflammation in PBMCs. Moreover, the cQTL in the GLIS3 locus co-localises with the previously found locus associated with Knee OA, suggesting that this Knee OA locus might be explained through an inflammatory mechanism. These results form a basis for further exploration of inflammatory mechanisms in OA.

Younger facial looks are associate with a lower likelihood of several age-related morbidities in the middle-aged to elderly.

BACKGROUND: Looking older for one's chronological age is associated with a higher mortality rate. Yet it remains unclear how perceived facial age relates to morbidity and the degree to which facial ageing reflects systemic ageing of the human body. OBJECTIVES: To investigate the association between ΔPA and age-related morbidities of different organ systems, where ΔPA represents the difference between perceived age (PA) and chronological age. METHODS: We performed a cross-sectional analysis on data from the Rotterdam Study, a population-based cohort study in the Netherlands. High-resolution facial photographs of 2679 men and women aged 51.5-87.8 years of European descent were used to assess PA. PA was estimated and scored in 5-year categories using these photographs by a panel of men and women who were blinded for chronological age and medical history. A linear mixed model was used to generate the mean PAs. The difference between the mean PA and chronological age was calculated (ΔPA), where a higher (positive) ΔPA means that the person looks younger for their age and a lower (negative) ΔPA that the person looks older. ΔPA was tested as a continuous variable for association with ageing-related morbidities including cardiovascular, pulmonary, ophthalmological, neurocognitive, renal, skeletal and auditory morbidities in separate regression analyses, adjusted for age and sex (model 1) and additionally for body mass index, smoking and sun exposure (model 2). RESULTS: We observed 5-year higher ΔPA (i.e. looking younger by 5 years for one's age) to be associated with less osteoporosis [odds ratio (OR) 0.76, 95% confidence interval (CI) 0.62-0.93], less chronic obstructive pulmonary disease (OR 0.85, 95% CI 0.77-0.95), less age-related hearing loss (model 2; B = -0.76, 95% CI -1.35 to -0.17) and fewer cataracts (OR 0.84, 95% CI 0.73-0.97), but with better global cognitive functioning (g-factor; model 2; B = 0.07, 95% CI 0.04-0.10). CONCLUSIONS: PA is associated with multiple morbidities and better cognitive function, suggesting that systemic ageing and cognitive ageing are, to an extent, externally visible in the human face.

A Genome-Wide Association Study Meta-Analysis of Alpha Angle Suggests Cam-Type Morphology May Be a Specific Feature of Hip Osteoarthritis in Older Adults.

OBJECTIVE: To examine the genetic architecture of cam morphology using alpha angle (AA) as a proxy measure and conduct an AA genome-wide association study (GWAS) followed by Mendelian randomization (MR) to evaluate its causal relationship with hip osteoarthritis (OA). METHODS: Observational analyses examined associations between AA measurements derived from hip dual x-ray absorptiometry (DXA) scans from the UK Biobank study and radiographic hip OA outcomes and subsequent total hip replacement. Following these analyses, an AA GWAS meta-analysis was performed (N = 44,214) using AA measurements previously derived in the Rotterdam Study. Linkage disequilibrium score regression assessed the genetic correlation between AA and hip OA. Genetic associations considered significant (P < 5 × 10-8 ) were used as AA genetic instrument for 2-sample MR analysis. RESULTS: DXA-derived AA showed expected associations between AA and radiographic hip OA (adjusted odds ratio [OR] 1.63 [95% confidence interval (95% CI) 1.58, 1.67]) and between AA and total hip replacement (adjusted hazard ratio 1.45 [95% CI 1.33, 1.59]) in the UK Biobank study cohort. The heritability of AA was 10%, and AA had a moderate genetic correlation with hip OA (rg  = 0.26 [95% CI 0.10, 0.43]). Eight independent genetic signals were associated with AA. Two-sample MR provided weak evidence of causal effects of AA on hip OA risk (inverse variance weighted OR 1.84 [95% CI 1.14, 2.96], P = 0.01). In contrast, genetic predisposition for hip OA had stronger evidence of a causal effect on increased AA (inverse variance weighted ß = 0.09 [95% CI 0.04, 0.13], P = 4.58 × 10-5 ). CONCLUSION: Expected observational associations between AA and related clinical outcomes provided face validity for the DXA-derived AA measurements. Evidence of bidirectional associations between AA and hip OA, particularly for risk of hip OA on AA, suggests that hip shape modeling secondary to a genetic predisposition to hip OA contributes to the well-established relationship between hip OA and cam morphology in older adults.

Epigenome-Wide Association Study Reveals CpG Sites Associated with Thyroid Function and Regulatory Effects on KLF9.

Background: Thyroid hormones play a key role in differentiation and metabolism and are known regulators of gene expression through both genomic and epigenetic processes including DNA methylation. The aim of this study was to examine associations between thyroid hormones and DNA methylation. Methods: We carried out a fixed-effect meta-analysis of epigenome-wide association study (EWAS) of blood DNA methylation sites from 8 cohorts from the ThyroidOmics Consortium, incorporating up to 7073 participants of both European and African ancestry, implementing a discovery and replication stage. Statistical analyses were conducted using normalized beta CpG values as dependent and log-transformed thyrotropin (TSH), free thyroxine, and free triiodothyronine levels, respectively, as independent variable in a linear model. The replicated findings were correlated with gene expression levels in whole blood and tested for causal influence of TSH and free thyroxine by two-sample Mendelian randomization (MR). Results: Epigenome-wide significant associations (p-value <1.1E-7) of three CpGs for free thyroxine, five for free triiodothyronine, and two for TSH concentrations were discovered and replicated (combined p-values = 1.5E-9 to 4.3E-28). The associations included CpG sites annotated to KLF9 (cg00049440) and DOT1L (cg04173586) that overlap with all three traits, consistent with hypothalamic-pituitary-thyroid axis physiology. Significant associations were also found for CpGs in FKBP5 for free thyroxine, and at CSNK1D/LINCO1970 and LRRC8D for free triiodothyronine. MR analyses supported a causal effect of thyroid status on DNA methylation of KLF9. DNA methylation of cg00049440 in KLF9 was inversely correlated with KLF9 gene expression in blood. The CpG at CSNK1D/LINC01970 overlapped with thyroid hormone receptor alpha binding peaks in liver cells. The total additive heritability of the methylation levels of the six significant CpG sites was between 25% and 57%. Significant methylation QTLs were identified for CpGs at KLF9, FKBP5, LRRC8D, and CSNK1D/LINC01970. Conclusions: We report novel associations between TSH, thyroid hormones, and blood-based DNA methylation. This study advances our understanding of thyroid hormone action particularly related to KLF9 and serves as a proof-of-concept that integrations of EWAS with other -omics data can provide a valuable tool for unraveling thyroid hormone signaling in humans by complementing and feeding classical in vitro and animal studies.

A systematic review of the sex differences in risk factors for knee osteoarthritis.

OBJECTIVES: Previous systematic reviews focused on the evidence of common risk factors for knee OA (KOA); however, the effect and strength of association between risk factors and KOA might be different between the two sexes. The aim of the present systematic review was to determine the current evidence on sex differences in the association between risk factors and KOA and their prevalence. METHODS: We searched the following electronic bibliographic databases: MEDLINE (PubMed), EMBASE and Web of Science. A methodological quality assessment was conducted independently by two researchers according to an adapted version of the standardized set of criteria known as the Newcastle-Ottawa Quality Assessment Scale (NOS). The NOS, a star system, was converted to three categories of quality. RESULTS: In total, 27 studies reported sex-specific risk estimates on several risk factors for KOA. Out of the 22 longitudinal cohort studies (except one nested case-control), 12 were of good quality and 10 were of fair quality. The five cross-sectional studies consisted of one of good, three of fair and one of poor quality. There was an indication of sex differences in risk factors leading to higher risk of KOA: high BMI, alcohol consumption, atherosclerosis and high vitamin E levels in women, and high physical activity, soft drink consumption and abdominal obesity in men. Knee injury, high blood pressure and low step rate seem to affect both women and men. CONCLUSION: More good quality studies are needed to assess sex differences in risk factors for KOA, especially for symptomatic/clinical OA.

Plasma proteomics identifies CRTAC1 as a biomarker for osteoarthritis severity and progression.

OBJECTIVES: The aim of this study was to identify biomarkers for radiographic OA severity and progression acting within the inflammation and metabolic pathways. METHODS: For 3517 Rotterdam Study participants, 184 plasma protein levels were measured using Olink inflammation and cardiometabolic panels. We studied associations with severity and progression of knee, hip and hand OA and a composite overall OA burden score by multivariable regression models, adjusting for age, sex, cell counts and BMI. RESULTS: We found 18 significantly associated proteins for overall OA burden, of which 5 stayed significant after multiple testing correction: circulating cartilage acidic protein 1 (CRTAC1), cartilage oligomeric matrix protein (COMP), thrombospondin 4, IL-18 receptor 1 (IL-18R1) and TNF ligand superfamily member 14. These proteins were also associated with progression of knee OA, with the exception of IL-18R1. The strongest association was found for the level of CRTAC1, with 1 s.d. increase in protein level resulting in an increase of 0.09 (95% CI 0.06, 0.12) in the overall OA Kellgren-Lawrence sum score (P = 2.9 × 10-8) in the model adjusted for age, sex, BMI and cell counts. This association was also present with the severity of OA in all three joints and progression of knee OA and was independent of BMI. We observed a stronger association for CRTAC1 with OA than for the well-known OA biomarker COMP. CONCLUSION: We identified several compelling biomarkers reflecting the overall OA burden and the increased risk for OA progression. CRTAC1 was the most compelling and robust biomarker for OA severity and progression. Such a biomarker may be used for disease monitoring.

Gut microbiome-wide association study of depressive symptoms.

Depression is one of the most poorly understood diseases due to its elusive pathogenesis. There is an urgency to identify molecular and biological mechanisms underlying depression and the gut microbiome is a novel area of interest. Here we investigate the relation of fecal microbiome diversity and composition with depressive symptoms in 1,054 participants from the Rotterdam Study cohort and validate these findings in the Amsterdam HELIUS cohort in 1,539 subjects. We identify association of thirteen microbial taxa, including genera Eggerthella, Subdoligranulum, Coprococcus, Sellimonas, Lachnoclostridium, Hungatella, Ruminococcaceae (UCG002, UCG003 and UCG005), LachnospiraceaeUCG001, Eubacterium ventriosum and Ruminococcusgauvreauiigroup, and family Ruminococcaceae with depressive symptoms. These bacteria are known to be involved in the synthesis of glutamate, butyrate, serotonin and gamma amino butyric acid (GABA), which are key neurotransmitters for depression. Our study suggests that the gut microbiome composition may play a key role in depression.

Accelerated menopausal changes as human disease model 'FOCUM' for the development of osteoarthritis and other degenerative disorders: protocol for a prospective cohort study.

INTRODUCTION: The incidence of degenerative disorders, including osteoarthritis (OA), increases rapidly in women after menopause. However, the influence of the menopause is still insufficiently investigated due to the slowness of menopausal transition. In this study, a novel human model is used in which it is expected that menopausal-related changes will occur faster. This is the Females discontinuing Oral Contraceptives Use at Menopausal age model. The ultimate aim is to link these changes to OA and other degenerative disorders, including cardiovascular diseases, diabetes, osteoporosis and tendinopathies. METHODS AND ANALYSIS: This is a pilot observational prospective cohort study with 2 years of follow-up. Fifty women aged 50-60 who use oral contraceptive (OC) and have the intention to stop are included. Measurements are performed once before stopping OC, and four times thereafter at 6 weeks, 6 months, 1 year and 2 years. At every time point, a questionnaire is filled in and a sample of blood is drawn. At the first and final time points, a physical examination, hand radiographs and a MRI scan of one knee are performed. The primary OA outcome is progression of the MRI Osteoarthritis Knee Score. Secondary OA outcomes are the development of clinical knee and hand OA, development of knee OA according to the MRI definition, and progression of radiographic features for hand OA. Principal component analysis will be used to assess which changes occur after stopping OC. Univariate and multivariate generalised estimating equation models will be used to test for associations between these components and OA. ETHICS AND DISSEMINATION: The study has been approved by the Medical Ethics Committee of the Erasmus MC University Medical Center Rotterdam (MEC-2019-0592). All participants must give informed consent before data collection. Results will be disseminated in national and international journals. TRIAL REGISTRATION NUMBER: NL70796.078.19.

Risk Assessment for Hip and Knee Osteoarthritis Using Polygenic Risk Scores.

OBJECTIVE: Polygenic risk scores (PRS) allow risk stratification using common single-nucleotide polymorphisms (SNPs), and clinical applications are currently explored for several diseases. This study was undertaken to assess the risk of hip and knee osteoarthritis (OA) using PRS. METHODS: We analyzed 12,732 individuals from a population-based cohort from the Rotterdam Study (n = 11,496), a clinical cohort (Cohort Hip and Cohort Knee [CHECK] study; n = 908), and a high-risk cohort of overweight women (Prevention of Knee OA in Overweight Females [PROOF] study; n = 328), for the association of the PRS with prevalence/incidence of radiographic OA, of clinical OA, and of total hip replacement (THR) or total knee replacement (TKR). The hip PRS and knee PRS contained 44 and 24 independent SNPs, respectively, and were derived from a recent genome-wide association study meta-analysis. Standardized PRS (with Z transformation) were used in all analyses. RESULTS: We found a stronger association of the PRS for clinically defined OA compared to radiographic OA phenotypes, and we observed the highest PRS risk stratification for TKR/THR. The odds ratio (OR) per SD was 1.3 for incident THR (95% confidence interval [95% CI] 1.1-1.5) and 1.6 (95% CI 1.3-1.9) for incident TKR in the Rotterdam Study. The knee PRS was associated with incident clinical knee OA in the CHECK study (OR 1.3 [95% CI 1.1-1.5]), but not for the PROOF study (OR 1.2 [95% CI 0.8-1.7]). The OR for OA increased gradually across the PRS distribution, up to 2.1 (95% CI 1.4-3.2) for individuals with the 10% highest PRS compared to the middle 50% of the PRS distribution. CONCLUSION: Our findings validated the association of PRS across OA definitions. Since OA is becoming frequent and primary prevention is not commonly applicable, PRS-based risk assessment could play a role in OA prevention. However, the utility of PRS is dependent on the setting. Further studies are needed to test the integration of genetic risk assessment in diverse health care settings.

Integrative analysis of clinical and epigenetic biomarkers of mortality.

DNA methylation (DNAm) has been reported to be associated with many diseases and with mortality. We hypothesized that the integration of DNAm with clinical risk factors would improve mortality prediction. We performed an epigenome-wide association study of whole blood DNAm in relation to mortality in 15 cohorts (n = 15,013). During a mean follow-up of 10 years, there were 4314 deaths from all causes including 1235 cardiovascular disease (CVD) deaths and 868 cancer deaths. Ancestry-stratified meta-analysis of all-cause mortality identified 163 CpGs in European ancestry (EA) and 17 in African ancestry (AA) participants at p < 1 × 10-7 , of which 41 (EA) and 16 (AA) were also associated with CVD death, and 15 (EA) and 9 (AA) with cancer death. We built DNAm-based prediction models for all-cause mortality that predicted mortality risk after adjusting for clinical risk factors. The mortality prediction model trained by integrating DNAm with clinical risk factors showed an improvement in prediction of cancer death with 5% increase in the C-index in a replication cohort, compared with the model including clinical risk factors alone. Mendelian randomization identified 15 putatively causal CpGs in relation to longevity, CVD, or cancer risk. For example, cg06885782 (in KCNQ4) was positively associated with risk for prostate cancer (Beta = 1.2, PMR  = 4.1 × 10-4 ) and negatively associated with longevity (Beta = -1.9, PMR  = 0.02). Pathway analysis revealed that genes associated with mortality-related CpGs are enriched for immune- and cancer-related pathways. We identified replicable DNAm signatures of mortality and demonstrated the potential utility of CpGs as informative biomarkers for prediction of mortality risk.

DNA methylation signature of chronic low-grade inflammation and its role in cardio-respiratory diseases.

We performed a multi-ethnic Epigenome Wide Association study on 22,774 individuals to describe the DNA methylation signature of chronic low-grade inflammation as measured by C-Reactive protein (CRP). We find 1,511 independent differentially methylated loci associated with CRP. These CpG sites show correlation structures across chromosomes, and are primarily situated in euchromatin, depleted in CpG islands. These genomic loci are predominantly situated in transcription factor binding sites and genomic enhancer regions. Mendelian randomization analysis suggests altered CpG methylation is a consequence of increased blood CRP levels. Mediation analysis reveals obesity and smoking as important underlying driving factors for changed CpG methylation. Finally, we find that an activated CpG signature significantly increases the risk for cardiometabolic diseases and COPD.

Shorter periconception maternal telomere length and the risk of congenital cardiac outflow defects in the offspring.

BACKGROUND: Congenital cardiac outflow defects (COD) are the largest group of congenital heart defects, with ventricular septal defect (VSD) as the most prevalent phenotype. Increased maternal age, excessive oxidative stress and inflammation are involved in the pathophysiology of COD and enhance telomere length (TL) shortening. We investigated the association between periconception maternal TL and the risk of having a child with COD. METHODS: From a multicentre case-control trial, 306 case mothers of a child with COD and 424 control mothers of a child without a congenital malformation were selected. Relative TL was measured by qPCR. Multivariable logistic regression was used to compute crude and adjusted odds ratios, per standard deviation decrease, between maternal T/S ratio and COD and VSD risk. Adjustments were made for maternal age. Additional adjustments were made in a second model. RESULTS: Shorter maternal relative TL was significantly associated with an OR of 1.29 (95% CI 1.04-1.61), p = .02, for the risk of VSD in offspring, which remained significant after an adjustment for maternal age (adjOR 1.25(95% CI 1.01-1.55), p = .04). No association between maternal TL and the risk of overall COD in offspring was observed. CONCLUSION: Shorter maternal relative TL is associated with an approximately 1.3-OR for the risk, per SD in relative TL shortening, of VSD in the offspring. These findings need further confirmation in other studies on the predictive value of maternal TL.