Examining the Impact of Imputation Errors on Fine-Mapping Using DNA Methylation QTL as a Model Trait

Genetic variants disrupting DNA methylation at CpG dinucleotides (CpG-SNP) provide a set of known causal variants to serve as models for testing fine-mapping methodology. We use 1716 CpG-SNPs to test three fine-mapping approaches (BIMBAM, BSLMM, and the J-test), assessing the impact of imputation errors and the choice of reference panel by using both whole-genome sequence (WGS), and genotype array data on the same individuals (n=1166). The choice of imputation reference panel had a strong effect on imputation accuracy, with the 1000 Genomes Phase 3 (1000G) reference panel (n=2504 from 26 populations) giving a mean non-reference discordance rate between imputed and sequenced genotypes of 3.2% compared to 1.6% when using the Haplotype Reference Consortium (HRC) reference panel (n=32470 Europeans). These imputation errors impacted on whether the CpG-SNP was included in the 95% credible set, with a difference of ∼ 23% and ∼ 7% between the WGS and the 1000G and HRC imputed datasets respectively. All of the fine-mapping methods failed to reach the expected 95% coverage of the CpG-SNP. This is attributed to secondary cis genetic effects that are unable to be statistically separated from the CpG-SNP, and through a masking mechanism where the effect of the methylation disrupting allele at the CpG-SNP is hidden by the effect of a nearby SNP that has strong LD with the CpG-SNP. The reduced accuracy in fine-mapping a known causal variant in a low level biological trait with imputed genetic data has implications for the study of higher order complex traits and disease

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

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

This research was supported by the Australian Research Council (DP160102400), the Australian National Health and Medical Research Council (1078037, 1078901, 1103418, 1107258, 1127440 and 1113400), and the Sylvia & Charles Viertel Charitable Foundation. Riccardo Marioni was supported by Alzheimer’s Research UK Major Project Grant [ARUK-PG2017B-10]. Generation Scotland received core support from the Chief Scientist Office of the Scottish Government Health Directorates [CZD/16/6] and the Scottish Funding Council [HR03006]. Genotyping and DNA methylation profiling of the GS:SFHS samples was carried out by the Genetics Core Laboratory at the Wellcome Trust Clinical Research Facility, Edinburgh, Scotland, and was funded by the Medical Research Council UK and the Wellcome Trust (Wellcome Trust Strategic Award “STratifying Resilience and Depression Longitudinally” ((STRADL) Reference 104036/Z/14/Z).Peer reviewedPublisher PD

Identification of novel risk loci, causal insights, and heritable risk for Parkinson's disease: a meta-analysis of genome-wide association studies

Background Genome-wide association studies (GWAS) in Parkinson's disease have increased the scope of biological knowledge about the disease over the past decade. We aimed to use the largest aggregate of GWAS data to identify novel risk loci and gain further insight into the causes of Parkinson's disease. Methods We did a meta-analysis of 17 datasets from Parkinson's disease GWAS available from European ancestry samples to nominate novel loci for disease risk. These datasets incorporated all available data. We then used these data to estimate heritable risk and develop predictive models of this heritability. We also used large gene expression and methylation resources to examine possible functional consequences as well as tissue, cell type, and biological pathway enrichments for the identified risk factors. Additionally, we examined shared genetic risk between Parkinson's disease and other phenotypes of interest via genetic correlations followed by Mendelian randomisation. Findings Between Oct 1, 2017, and Aug 9, 2018, we analysed 7·8 million single nucleotide polymorphisms in 37 688 cases, 18 618 UK Biobank proxy-cases (ie, individuals who do not have Parkinson's disease but have a first degree relative that does), and 1·4 million controls. We identified 90 independent genome-wide significant risk signals across 78 genomic regions, including 38 novel independent risk signals in 37 loci. These 90 variants explained 16–36% of the heritable risk of Parkinson's disease depending on prevalence. Integrating methylation and expression data within a Mendelian randomisation framework identified putatively associated genes at 70 risk signals underlying GWAS loci for follow-up functional studies. Tissue-specific expression enrichment analyses suggested Parkinson's disease loci were heavily brain-enriched, with specific neuronal cell types being implicated from single cell data. We found significant genetic correlations with brain volumes (false discovery rate-adjusted p=0·0035 for intracranial volume, p=0·024 for putamen volume), smoking status (p=0·024), and educational attainment (p=0·038). Mendelian randomisation between cognitive performance and Parkinson's disease risk showed a robust association (p=8·00 × 10−7). Interpretation These data provide the most comprehensive survey of genetic risk within Parkinson's disease to date, to the best of our knowledge, by revealing many additional Parkinson's disease risk loci, providing a biological context for these risk factors, and showing that a considerable genetic component of this disease remains unidentified. These associations derived from European ancestry datasets will need to be followed-up with more diverse data. Funding The National Institute on Aging at the National Institutes of Health (USA), The Michael J Fox Foundation, and The Parkinson's Foundation (see appendix for full list of funding sources)

The genetic architecture of the human cerebral cortex

The cerebral cortex underlies our complex cognitive capabilities, yet little is known about the specific genetic loci that influence human cortical structure. To identify genetic variants that affect cortical structure, we conducted a genome-wide association meta-analysis of brain magnetic resonance imaging data from 51,665 individuals. We analyzed the surface area and average thickness of the whole cortex and 34 regions with known functional specializations. We identified 199 significant loci and found significant enrichment for loci influencing total surface area within regulatory elements that are active during prenatal cortical development, supporting the radial unit hypothesis. Loci that affect regional surface area cluster near genes in Wnt signaling pathways, which influence progenitor expansion and areal identity. Variation in cortical structure is genetically correlated with cognitive function, Parkinson's disease, insomnia, depression, neuroticism, and attention deficit hyperactivity disorder

The genetic architecture of the human cerebral cortex

The cerebral cortex underlies our complex cognitive capabilities, yet little is known about the specific genetic loci that influence human cortical structure. To identify genetic variants that affect cortical structure, we conducted a genome-wide association meta-analysis of brain magnetic resonance imaging data from 51,665 individuals. We analyzed the surface area and average thickness of the whole cortex and 34 regions with known functional specializations. We identified 199 significant loci and found significant enrichment for loci influencing total surface area within regulatory elements that are active during prenatal cortical development, supporting the radial unit hypothesis. Loci that affect regional surface area cluster near genes in Wnt signaling pathways, which influence progenitor expansion and areal identity. Variation in cortical structure is genetically correlated with cognitive function, Parkinson's disease, insomnia, depression, neuroticism, and attention deficit hyperactivity disorder

Meta-analysis of genome-wide DNA methylation identifies shared associations across neurodegenerative disorders

Background People with neurodegenerative disorders show diverse clinical syndromes, genetic heterogeneity, and distinct brain pathological changes, but studies report overlap between these features. DNA methylation (DNAm) provides a way to explore this overlap and heterogeneity as it is determined by the combined effects of genetic variation and the environment. In this study, we aim to identify shared blood DNAm differences between controls and people with Alzheimer’s disease, amyotrophic lateral sclerosis, and Parkinson’s disease. Results We use a mixed-linear model method (MOMENT) that accounts for the effect of (un)known confounders, to test for the association of each DNAm site with each disorder. While only three probes are found to be genome-wide significant in each MOMENT association analysis of amyotrophic lateral sclerosis and Parkinson’s disease (and none with Alzheimer’s disease), a fixed-effects meta-analysis of the three disorders results in 12 genome-wide significant differentially methylated positions. Predicted immune cell-type proportions are disrupted across all neurodegenerative disorders. Protein inflammatory markers are correlated with profile sum-scores derived from disease-associated immune cell-type proportions in a healthy aging cohort. In contrast, they are not correlated with MOMENT DNAm-derived profile sum-scores, calculated using effect sizes of the 12 differentially methylated positions as weights. Conclusions We identify shared differentially methylated positions in whole blood between neurodegenerative disorders that point to shared pathogenic mechanisms. These shared differentially methylated positions may reflect causes or consequences of disease, but they are unlikely to reflect cell-type proportion differences

Epigenome wide association study of response to methotrexate in early rheumatoid arthritis patients.

AimTo identify differentially methylated positions (DMPs) and regions (DMRs) that predict response to Methotrexate (MTX) in early rheumatoid arthritis (RA) patients.Materials and methodsDNA from baseline peripheral blood mononuclear cells was extracted from 72 RA patients. DNA methylation, quantified using the Infinium MethylationEPIC, was assessed in relation to response to MTX (combination) therapy over the first 3 months.ResultsBaseline DMPs associated with response were identified; including hits previously described in RA. Additionally, 1309 DMR regions were observed. However, none of these findings were genome-wide significant. Likewise, no specific pathways were related to response, nor could we replicate associations with previously identified DMPs.ConclusionNo baseline genome-wide significant differences were identified as biomarker for MTX (combination) therapy response; hence meta-analyses are required

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

OBJECTIVES: The aim of this study was to identify biomarkers for radiographic osteoarthritis severity and progression acting within the inflammation and metabolic pathways. METHODS: For 3,517 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, hand osteoarthritis, 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 osteoarthritis burden, of which 5 stayed significant after multiple testing correction: circulating Cartilage acidic protein 1 (CRTAC1), Cartilage oligomeric matrix protein (COMP), Thrombospondin 4 (THBS4), Interleukin 18 receptor 1 (IL18R1) and Tumor necrosis factor ligand superfamily member 14 (TNFSF14). These proteins were also associated with progression of knee OA, with the exception of IL18R1. The strongest association was found for the level of CRTAC1 with 1 SD increase in protein level resulting in an increase of 0.09 (95%CI : 0.06-0.12) in overall osteoarthritis KLsum-score (p= 2.9x10-8), in the model adjusted for age, sex, BMI and cell counts. This association was also present with severity of osteoarthritis in all three joints, progression of knee osteoarthritis and was independent of BMI. We observed a stronger association for CRTAC1 with osteoarthritis than for the well-known osteoarthritis-biomarker COMP. CONCLUSION: We identified several compelling biomarkers reflecting overall osteoarthritis burden and increased risk for osteoarthritis progression. CRTAC1 was the most compelling and robust biomarker for osteoarthritis severity and progression. Such a biomarker may be used for disease monitoring

Plasma proteomic signature of fatty liver disease:The Rotterdam Study

Background and Aims: Fatty liver disease (FLD) is caused by excess fat in the liver, and its global prevalence exceeds 33%. The role of protein expression on the pathogenesis of FLD and accompanied fibrosis and its potential as a disease biomarker is currently not clear. Hence, we aimed to identify plasma proteomics associated with FLD and fibrosis using population-based data. Approach and Results:Blood samples were collected from 2578 participants from the population-based Rotterdam Study cohort. The proximity extension assay reliably measured plasma levels of 171 cardiometabolic and inflammatory-related proteins (Olink Proteomics). FLD was assessed by ultrasound, and fibrosis by transient elastography. Logistic regression models quantified the association of plasma proteomics with FLD and fibrosis. In addition, we aimed to validate our results in liver organoids. The cross-sectional analysis identified 27 proteins significantly associated with FLD surpassing the Bonferroni-corrected p&lt;2.92×10-4. The strongest association was observed for FGF-21 (β=0.45, p=1.07×10-18) and carboxylesterase 1 (CES1) protein (β=0.66, p=4.91×10-40). Importantly, 15 of the 27 proteins significantly associated with FLD were also associated with liver fibrosis. Finally, consistent with plasma proteomic profiling, we found the expression levels of IL-18 receptor 1 (IL-18R1) and CES1 to be upregulated in an FLD model of 3-dimensional culture human liver organoids. Conclusions: Among the general population, several inflammatory and cardiometabolic plasma proteins were associated with FLD and fibrosis. Particularly, plasma levels of FGF-21, IL-18R1, and CES1 were largely dependent on the presence of FLD and fibrosis and may therefore be important in their pathogenesis.</p