Cross-National and Cross-Generational Evidence That Educational Attainment May Slow the Pace of Aging in European-Descent Individuals

OBJECTIVES

Individuals with more education are at lower risk of developing multiple, different age-related diseases than their less-educated peers. A reason for this might be that individuals with more education age slower. There are 2 complications in testing this hypothesis. First, there exists no definitive measure of biological aging. Second, shared genetic factors contribute toward both lower educational attainment and the development of age-related diseases. Here, we tested whether the protective effect of educational attainment was associated with the pace of aging after accounting for genetic factors.

METHODS

We examined data from 5 studies together totaling almost 17,000 individuals with European ancestry born in different countries during different historical periods, ranging in age from 16 to 98 years old. To assess the pace of aging, we used DunedinPACE, a DNA methylation algorithm that reflects an individual's rate of aging and predicts age-related decline and Alzheimer's disease and related disorders. To assess genetic factors related to education, we created a polygenic score based on the results of a genome-wide association study of educational attainment.

RESULTS

Across the 5 studies, and across the life span, higher educational attainment was associated with a slower pace of aging even after accounting for genetic factors (meta-analysis effect size = -0.20; 95% confidence interval [CI]: -0.30 to -0.10; p = .006). Further, this effect persisted after taking into account tobacco smoking (meta-analysis effect size = -0.13; 95% CI: -0.21 to -0.05; p = .01).

DISCUSSION

These results indicate that higher levels of education have positive effects on the pace of aging, and that the benefits can be realized irrespective of individuals' genetics.

Opioid medication use and blood DNA methylation: epigenome-wide association meta-analysis

Aim: To identify differential methylation related to prescribed opioid use. Methods: This study examined whether blood DNA methylation, measured using Illumina arrays, differs by recent opioid medication use in four population-based cohorts. We meta-analyzed results (282 users; 10,560 nonusers) using inverse-variance weighting. Results: Differential methylation (false discovery rate <0.05) was observed at six CpGs annotated to the following genes: KIAA0226, CPLX2, TDRP, RNF38, TTC23 and GPR179. Integrative epigenomic analyses linked implicated loci to regulatory elements in blood and/or brain. Additionally, 74 CpGs were differentially methylated in males or females. Methylation at significant CpGs correlated with gene expression in blood and/or brain. Conclusion: This study identified DNA methylation related to opioid medication use in general populations. The results could inform the development of blood methylation biomarkers of opioid use.

Inequalities in healthcare disruptions during the COVID-19 pandemic: evidence from 12 UK population-based longitudinal studies

OBJECTIVES

We investigated associations between multiple sociodemographic characteristics (sex, age, occupational social class, education and ethnicity) and self-reported healthcare disruptions during the early stages of the COVID-19 pandemic.

DESIGN

Coordinated analysis of prospective population surveys.

SETTING

Community-dwelling participants in the UK between April 2020 and January 2021.

PARTICIPANTS

Over 68 000 participants from 12 longitudinal studies.

OUTCOMES

Self-reported healthcare disruption to medication access, procedures and appointments.

RESULTS

Prevalence of healthcare disruption varied substantially across studies: between 6% and 32% reported any disruption, with 1%-10% experiencing disruptions in medication, 1%-17% experiencing disruption in procedures and 4%-28% experiencing disruption in clinical appointments. Females (OR 1.27; 95% CI 1.15 to 1.40; I²=54%), older persons (eg, OR 1.39; 95% CI 1.13 to 1.72; I²=77% for 65-75 years vs 45-54 years) and ethnic minorities (excluding white minorities) (OR 1.19; 95% CI 1.05 to 1.35; I²=0% vs white) were more likely to report healthcare disruptions. Those in a more disadvantaged social class were also more likely to report healthcare disruptions (eg, OR 1.17; 95% CI 1.08 to 1.27; I²=0% for manual/routine vs managerial/professional), but no clear differences were observed by education. We did not find evidence that these associations differed by shielding status.

CONCLUSIONS

Healthcare disruptions during the COVID-19 pandemic could contribute to the maintenance or widening of existing health inequalities.

Epigenetic scores for the circulating proteome as tools for disease prediction

Protein biomarkers have been identified across many age-related morbidities. However, characterising epigenetic influences could further inform disease predictions. Here, we leverage epigenome-wide data to study links between the DNA methylation (DNAm) signatures of the circulating proteome and incident diseases. Using data from four cohorts, we trained and tested epigenetic scores (EpiScores) for 953 plasma proteins, identifying 109 scores that explained between 1% and 58% of the variance in protein levels after adjusting for known protein quantitative trait loci (pQTL) genetic effects. By projecting these EpiScores into an independent sample (Generation Scotland; n = 9537) and relating them to incident morbidities over a follow-up of 14 years, we uncovered 137 EpiScore-disease associations. These associations were largely independent of immune cell proportions, common lifestyle and health factors, and biological aging. Notably, we found that our diabetes-associated EpiScores highlighted previous top biomarker associations from proteome-wide assessments of diabetes. These EpiScores for protein levels can therefore be a valuable resource for disease prediction and risk stratification.

A proteomic survival predictor for COVID-19 patients in intensive care

Global healthcare systems are challenged by the COVID-19 pandemic. There is a need to optimize allocation of treatment and resources in intensive care, as clinically established risk assessments such as SOFA and APACHE II scores show only limited performance for predicting the survival of severely ill COVID-19 patients. Additional tools are also needed to monitor treatment, including experimental therapies in clinical trials. Comprehensively capturing human physiology, we speculated that proteomics in combination with new data-driven analysis strategies could produce a new generation of prognostic discriminators. We studied two independent cohorts of patients with severe COVID-19 who required intensive care and invasive mechanical ventilation. SOFA score, Charlson comorbidity index, and APACHE II score showed limited performance in predicting the COVID-19 outcome. Instead, the quantification of 321 plasma protein groups at 349 timepoints in 50 critically ill patients receiving invasive mechanical ventilation revealed 14 proteins that showed trajectories different between survivors and non-survivors. A predictor trained on proteomic measurements obtained at the first time point at maximum treatment level (i.e. WHO grade 7), which was weeks before the outcome, achieved accurate classification of survivors (AUROC 0.81). We tested the established predictor on an independent validation cohort (AUROC 1.0). The majority of proteins with high relevance in the prediction model belong to the coagulation system and complement cascade. Our study demonstrates that plasma proteomics can give rise to prognostic predictors substantially outperforming current prognostic markers in intensive care.

Early life predictors of late life cerebral small vessel disease in four prospective cohort studies

Development of cerebral small vessel disease, a major cause of stroke and dementia, may be influenced by early life factors. It is unclear whether these relationships are independent of each other, of adult socio-economic status or of vascular risk factor exposures. We examined associations between factors from birth (ponderal index, birth weight), childhood (IQ, education, socio-economic status), adult small vessel disease, and brain volumes, using data from four prospective cohort studies: STratifying Resilience And Depression Longitudinally (STRADL) (n = 1080; mean age = 59 years); the Dutch Famine Birth Cohort (n = 118; mean age = 68 years); the Lothian Birth Cohort 1936 (LBC1936; n = 617; mean age = 73 years), and the Simpson's cohort (n = 110; mean age = 78 years). We analysed each small vessel disease feature individually and summed to give a total small vessel disease score (range 1-4) in each cohort separately, then in meta-analysis, adjusted for vascular risk factors and adult socio-economic status. Higher birth weight was associated with fewer lacunes [odds ratio (OR) per 100 g = 0.93, 95% confidence interval (CI) = 0.88 to 0.99], fewer infarcts (OR = 0.94, 95% CI = 0.89 to 0.99), and fewer perivascular spaces (OR = 0.95, 95% CI = 0.91 to 0.99). Higher childhood IQ was associated with lower white matter hyperintensity burden (OR per IQ point = 0.99, 95% CI 0.98 to 0.998), fewer infarcts (OR = 0.98, 95% CI = 0.97 to 0.998), fewer lacunes (OR = 0.98, 95% CI = 0.97 to 0.999), and lower total small vessel disease burden (OR = 0.98, 95% CI = 0.96 to 0.999). Low education was associated with more microbleeds (OR = 1.90, 95% CI = 1.33 to 2.72) and lower total brain volume (mean difference = -178.86 cm3, 95% CI = -325.07 to -32.66). Low childhood socio-economic status was associated with fewer lacunes (OR = 0.62, 95% CI = 0.40 to 0.95). Early life factors are associated with worse small vessel disease in later life, independent of each other, vascular risk factors and adult socio-economic status. Risk for small vessel disease may originate in early life and provide a mechanistic link between early life factors and risk of stroke and dementia. Policies investing in early child development may improve lifelong brain health and contribute to the prevention of dementia and stroke in older age.

Genome-wide methylation data improves dissection of the effect of smoking on body mass index

Variation in obesity-related traits has a genetic basis with heritabilities between 40 and 70%. While the global obesity pandemic is usually associated with environmental changes related to lifestyle and socioeconomic changes, most genetic studies do not include all relevant environmental covariates, so the genetic contribution to variation in obesity-related traits cannot be accurately assessed. Some studies have described interactions between a few individual genes linked to obesity and environmental variables but there is no agreement on their total contribution to differences between individuals. Here we compared self-reported smoking data and a methylation-based proxy to explore the effect of smoking and genome-by-smoking interactions on obesity related traits from a genome-wide perspective to estimate the amount of variance they explain. Our results indicate that exploiting omic measures can improve models for complex traits such as obesity and can be used as a substitute for, or jointly with, environmental records to better understand causes of disease.

Using symptom-based case predictions to identify host genetic factors that contribute to COVID-19 susceptibility

Epidemiological and genetic studies on COVID-19 are currently hindered by inconsistent and limited testing policies to confirm SARS-CoV-2 infection. Recently, it was shown that it is possible to predict COVID-19 cases using cross-sectional self-reported disease-related symptoms. Here, we demonstrate that this COVID-19 prediction model has reasonable and consistent performance across multiple independent cohorts and that our attempt to improve upon this model did not result in improved predictions. Using the existing COVID-19 prediction model, we then conducted a GWAS on the predicted phenotype using a total of 1,865 predicted cases and 29,174 controls. While we did not find any common, large-effect variants that reached genome-wide significance, we do observe suggestive genetic associations at two SNPs (rs11844522, p = 1.9x10-7; rs5798227, p = 2.2x10-7). Explorative analyses furthermore suggest that genetic variants associated with other viral infectious diseases do not overlap with COVID-19 susceptibility and that severity of COVID-19 may have a different genetic architecture compared to COVID-19 susceptibility. This study represents a first effort that uses a symptom-based predicted phenotype as a proxy for COVID-19 in our pursuit of understanding the genetic susceptibility of the disease. We conclude that the inclusion of symptom-based predicted cases could be a useful strategy in a scenario of limited testing, either during the current COVID-19 pandemic or any future viral outbreak.

Cohort profile: DOLORisk Dundee: a longitudinal study of chronic neuropathic pain

PURPOSE

Neuropathic pain is a common disorder of the somatosensory system that affects 7%-10% of the general population. The disorder places a large social and economic burden on patients as well as healthcare services. However, not everyone with a relevant underlying aetiology develops corresponding pain. DOLORisk Dundee, a European Union-funded cohort, part of the multicentre DOLORisk consortium, was set up to increase current understanding of this variation in onset. In particular, the cohort will allow exploration of psychosocial, clinical and genetic predictors of neuropathic pain onset.

PARTICIPANTS

DOLORisk Dundee has been constructed by rephenotyping two pre-existing Scottish population cohorts for neuropathic pain using a standardised 'core' study protocol: Genetics of Diabetes Audit and Research in Tayside Scotland (GoDARTS) (n=5236) consisting of predominantly type 2 diabetics from the Tayside region, and Generation Scotland: Scottish Family Health Study (GS:SFHS; n=20 221). Rephenotyping was conducted in two phases: a baseline postal survey and a combined postal and online follow-up survey. DOLORisk Dundee consists of 9155 participants (GoDARTS=1915; GS:SFHS=7240) who responded to the baseline survey, of which 6338 (69.2%; GoDARTS=1046; GS:SFHS=5292) also responded to the follow-up survey (18 months later).

FINDINGS TO DATE

At baseline, the proportion of those with chronic neuropathic pain (Douleur Neuropathique en 4 Questions questionnaire score ≥3, duration ≥3 months) was 30.5% in GoDARTS and 14.2% in Generation Scotland. Electronic record linkage enables large scale genetic association studies to be conducted and risk models have been constructed for neuropathic pain.

FUTURE PLANS

The cohort is being maintained by an access committee, through which collaborations are encouraged. Details of how to do this will be available on the study website (http://dolorisk.eu/). Further follow-up surveys of the cohort are planned and funding applications are being prepared to this effect. This will be conducted in harmony with similar pain rephenotyping of UK Biobank.

Investigating genetic links between grapheme-colour synaesthesia and neuropsychiatric traits

Synaesthesia is a neurological phenomenon affecting perception, where triggering stimuli (e.g. letters and numbers) elicit unusual secondary sensory experiences (e.g. colours). Family-based studies point to a role for genetic factors in the development of this trait. However, the contributions of common genomic variation to synaesthesia have not yet been investigated. Here, we present the SynGenes cohort, the largest genotyped collection of unrelated people with grapheme-colour synaesthesia (n = 723). Synaesthesia has been associated with a range of other neuropsychological traits, including enhanced memory and mental imagery, as well as greater sensory sensitivity. Motivated by the prior literature on putative trait overlaps, we investigated polygenic scores derived from published genome-wide scans of schizophrenia and autism spectrum disorder (ASD), comparing our SynGenes cohort to 2181 non-synaesthetic controls. We found a very slight association between schizophrenia polygenic scores and synaesthesia (Nagelkerke's R2 = 0.0047, empirical p = 0.0027) and no significant association for scores related to ASD (Nagelkerke's R2 = 0.00092, empirical p = 0.54) or body mass index (R2 = 0.00058, empirical p = 0.60), included as a negative control. As sample sizes for studying common genomic variation continue to increase, genetic investigations of the kind reported here may yield novel insights into the shared biology between synaesthesia and other traits, to complement findings from neuropsychology and brain imaging. This article is part of a discussion meeting issue 'Bridging senses: novel insights from synaesthesia'.

Improved precision of epigenetic clock estimates across tissues and its implication for biological ageing

BACKGROUND

DNA methylation changes with age. Chronological age predictors built from DNA methylation are termed 'epigenetic clocks'. The deviation of predicted age from the actual age ('age acceleration residual', AAR) has been reported to be associated with death. However, it is currently unclear how a better prediction of chronological age affects such association.

METHODS

In this study, we build multiple predictors based on training DNA methylation samples selected from 13,661 samples (13,402 from blood and 259 from saliva). We use the Lothian Birth Cohorts of 1921 (LBC1921) and 1936 (LBC1936) to examine whether the association between AAR (from these predictors) and death is affected by (1) improving prediction accuracy of an age predictor as its training sample size increases (from 335 to 12,710) and (2) additionally correcting for confounders (i.e., cellular compositions). In addition, we investigated the performance of our predictor in non-blood tissues.

RESULTS

We found that in principle, a near-perfect age predictor could be developed when the training sample size is sufficiently large. The association between AAR and mortality attenuates as prediction accuracy increases. AAR from our best predictor (based on Elastic Net, https://github.com/qzhang314/DNAm-based-age-predictor ) exhibits no association with mortality in both LBC1921 (hazard ratio = 1.08, 95% CI 0.91-1.27) and LBC1936 (hazard ratio = 1.00, 95% CI 0.79-1.28). Predictors based on small sample size are prone to confounding by cellular compositions relative to those from large sample size. We observed comparable performance of our predictor in non-blood tissues with a multi-tissue-based predictor.

CONCLUSIONS

This study indicates that the epigenetic clock can be improved by increasing the training sample size and that its association with mortality attenuates with increased prediction of chronological age.

Genome-wide meta-analysis of 158,000 individuals of European ancestry identifies three loci associated with chronic back pain

Back pain is the #1 cause of years lived with disability worldwide, yet surprisingly little is known regarding the biology underlying this symptom. We conducted a genome-wide association study (GWAS) meta-analysis of chronic back pain (CBP). Adults of European ancestry were included from 15 cohorts in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and from the UK Biobank interim data release. CBP cases were defined as those reporting back pain present for ≥3-6 months; non-cases were included as comparisons ("controls"). Each cohort conducted genotyping using commercially available arrays followed by imputation. GWAS used logistic regression models with additive genetic effects, adjusting for age, sex, study-specific covariates, and population substructure. The threshold for genome-wide significance in the fixed-effect inverse-variance weighted meta-analysis was p<5×10(-8). Suggestive (p<5×10(-7)) and genome-wide significant (p<5×10(-8)) variants were carried forward for replication or further investigation in the remaining UK Biobank participants not included in the discovery sample. The discovery sample comprised 158,025 individuals, including 29,531 CBP cases. A genome-wide significant association was found for the intronic variant rs12310519 in SOX5 (OR 1.08, p = 7.2×10(-10)). This was subsequently replicated in 283,752 UK Biobank participants not included in the discovery sample, including 50,915 cases (OR 1.06, p = 5.3×10(-11)), and exceeded genome-wide significance in joint meta-analysis (OR 1.07, p = 4.5×10(-19)). We found suggestive associations at three other loci in the discovery sample, two of which exceeded genome-wide significance in joint meta-analysis: an intergenic variant, rs7833174, located between CCDC26 and GSDMC (OR 1.05, p = 4.4×10(-13)), and an intronic variant, rs4384683, in DCC (OR 0.97, p = 2.4×10(-10)). In this first reported meta-analysis of GWAS for CBP, we identified and replicated a genetic locus associated with CBP (SOX5). We also identified 2 other loci that reached genome-wide significance in a 2-stage joint meta-analysis (CCDC26/GSDMC and DCC).