Back

Amplification of the epigenetic (gestational) age acceleration signal

Lee, Y.; Bohlin, J.

2021-02-02 bioinformatics
10.1101/2021.02.02.429312 bioRxiv
Show abstract

BackgroundEpigenetic (gestational) age acceleration (E(G)AA) is associated with environmental exposures and health outcomes in humans. However, E(G)AA is the residual term from a regression of epigenetic age (outcome) on chronological (gestational) age (predictor) and therefore strongly obscured by noise from multiple sources. Here, we propose a simple procedure, based on regression, principal component analysis (PCA), and the Lasso, that amplifies E(G)AA signals. More specifically, we first regress given (gestational) age against each CpG used for epigenetic (gestational) age prediction. The CpGs are typically taken from one of several epigenetic clocks available. PCA is subsequently performed on the resulting matrix of residual vectors for each CpG as it projects the E(G)AA signal onto perpendicular principal components (PCs), thereby separating signal from noise. Finally, we use the Lasso to select PCs associated with an outcome of interest. We apply our method to previous studies: EAA in patients with Downs syndrome and Werners syndrome and EGAA of newborns exposed to prenatal smoking as well as associations with maternal BMI. ResultsThe extracted EAA components computed using our proposed procedure revealed a significant association with Downs syndrome (PB<0.05, Bonferroni adjusted for multiple testing) as well as for Werners Syndrome (PB<0.05). For EGAA we find a significant association with maternal prenatal smoking (PB<0.05, also Bonferroni adjusted) and maternal BMI (PB<0.05). Additionally, by examining the loadings of the PCs of interest, and contrary to residual EGAA, our method can identify implicated CpGs. ConclusionsOur findings suggest that our proposed procedure leads to a remarkable amplification of the E(G)AA signal. Furthermore, our method reveals that E(G)AA is a composite signal that can be driven by multiple independent factors.

Matching journals

The top 1 journal accounts for 50% of the predicted probability mass.

50% of probability mass above

"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.