Computational Methods for Analysis of DNA Methylation Data, Volume II

Introduction: Highly active antiretroviral therapy (HAART) helps improve some measures of accelerated epigenetic aging in persons living with HIV (PLWH), but its overall impact on the epigenome is not fully understood.

Methods: In this study, we analyzed the DNA methylation profiles of PLWH (n = 187) shortly before and approximately 2–3 years after they started HAART, as well as matched seronegative (SN) controls (n = 187), taken at two time intervals. Our aim was to identify specific CpGs and biologic pathways associated with HIV infection and initiation of HAART. Additionally, we attempted to identify epigenetic changes associated with HAART initiation that were independent of HIV-associated changes, using matched HIV seronegative (SN) controls (matched on age, hepatitis C status, and interval between visits) to identify CpGs that did not differ between PLWH and SN pre-HAART but were significantly associated with HAART initiation while being unrelated to HIV viral load. Epigenome-wide association studies (EWAS) on >850,000 CpG sites were performed using pre- and post-HAART samples from PLWH. The results were then annotated using the Genomic Regions Enrichment of Annotations Tool (GREAT).

Results: When only pre- and post-HAART visits in PLWH were compared, gene ontologies related to immune function and diseases related to immune function were significant, though with less significance for PLWH with detectable HIV viral loads (>50 copies/mL) at the post-HAART visit. To specifically elucidate the effects of HAART separately from HIV-induced methylation changes, we performed EWAS of HAART while also controlling for HIV viral load, and found gene ontologies associated with transplant rejection, transplant-related diseases, and other immunologic signatures. Additionally, we performed a more focused analysis that examined CpGs reaching genome-wide significance (p < 1 × 10⁻⁷) from the viral load-controlled EWAS that did not differ between all PLWH and matched SN controls pre-HAART. These CpGs were found to be near genes that play a role in retroviral drug metabolism, diffuse large B cell lymphoma proliferation, and gastric cancer metastasis.

Discussion: Overall, this study provides insight into potential biological functions associated with DNA methylation changes induced by HAART initiation in persons living with HIV.

Introduction: DNA methylation clocks presents advantageous characteristics with respect to the ambitious goal of identifying very early markers of disease, based on the concept that accelerated ageing is a reliable predictor in this sense.

Methods: Such tools, being epigenomic based, are expected to be conditioned by sex and tissue specificities, and this work is about quantifying this dependency as well as that from the regression model and the size of the training set.

Results: Our quantitative results indicate that elastic-net penalization is the best performing strategy, and better so when—unsurprisingly—the data set is bigger; sex does not appear to condition clocks performances and tissue specific clocks appear to perform better than generic blood clocks. Finally, when considering all trained clocks, we identified a subset of genes that, to the best of our knowledge, have not been presented yet and might deserve further investigation: CPT1A, MMP15, SHROOM3, SLIT3, and SYNGR.

Conclusion: These factual starting points can be useful for the future medical translation of clocks and in particular in the debate between multi-tissue clocks, generally trained on a large majority of blood samples, and tissue-specific clocks.