Dynamic Methylome Modification associated with mutational signatures in ageing and etiology of disease

Reversible epigenetic changes within the loci of genes that regulate critical cell processes have recently emerged as important biomarkers of disease pathology. It is then natural to consider the consequences for population health risk of such epigenetic changes during the aging process. Specifically, the interplay between dynamic methylation changes that accompany aging and mutations that accrue in an individuals genome over time need further investigation. The current study investigated the role of dynamic methylation acting together with gene variants in an individual over time to gain insight into the evolving epigenome-genome interplay that affects biochemical pathways controlling physiological processes during aging. We completed a whole-genome methylation and variant analysis in a non-smoking Zoroastrian-Parsi individual, collecting two samples, 12 years apart (at 53 and 65 years respectively) (ZPMetG-Hv2a-1A (old, t0), ZPMetG-Hv2a-1B (recent, t0+12)) and analysing them using a GridION Nanopore sequencer at 13X genome coverage overall. We further identified the single nucleotide variants (SNVs) and indels in known CpG islands by employing the Genome Analysis Tool Kit (GATK) and MuTect2 variant-caller pipeline with the GRCh37 (patch 13) human genome as reference. We found 5258 disease-relevant genes that had been differentially methylated in this individual over 12 years. Employing the GATK pipeline, we found 24,948 genes, corresponding to 4,58,148 variants, specific to ZPMetG-Hv2a-1B, indicating the presence of variants that had accrued over time. A fraction of the gene variants (242/24948) occurred within the CpG regions that were differentially methylated, with 67/247 exactly coincident with a CpG site. Our analysis yielded a critical cluster of 10 genes that were each significantly methylated and had variants at the CpG site or the {+/-}4 bp CpG region window. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment network analysis as well as Reactome and STRING analysis of gene-specific variants indicated an impact on biological processes regulating the immune system, disease networks implicated in cancer and neurodegenerative diseases, and transcriptional control of processes regulating cellular senescence and longevity. Additional analysis of mutational signatures indicated a majority of C>T transitions followed by T>C transitions in the more recent sample, ZPMetG-Hv2a-1B. Our current study provides additional insight into the aging methylome over time and the interplay between different methylation and gene variants in the etiology of disease.