- Open Access
Longitudinal analysis of healthy colon establishes aspirin as a suppressor of cancer-related epigenetic aging
Clinical Epigenetics volume 12, Article number: 164 (2020)
Colon cancer (CC) is the third most common cancer worldwide, highlighting the importance of developing effective prevention strategies. Accumulating evidence supports that aspirin use reduces CC incidence. We reported previously that aspirin suppresses age-associated and CC-relevant DNA methylation (DNAm) in healthy colon. Here we addressed the aspirin’s effectiveness in longitudinal cohort.
We measured genome-wide DNAm in 124 healthy normal mucosa samples taken at baseline (time point 1, t1) and after 10-years follow-up (time point 2, t2) from a longitudinal female screening cohort. We investigated the time-dependent methylation drift in aspirin users and nonusers using multivariable regression and related the modulatory effect of aspirin to colonic epigenome-aging and CC.
Over time, compared to nonusers, long-term (≥ 2 years) aspirin users showed less hypermethylated CpGs (proximal: 17% vs. 87%; distal: 16% vs. 70%) and more hypomethylated CpGs (proximal: 83% vs. 13%; distal: 84% vs. 30%). Overall, users showed 2% (P = 0.02) less mean methylation levels than nonusers in proximal colon and displayed repressed methylation age (mAge). Methylation loss in users occurred at several CC-specific tumor suppressors that gained methylation in nonusers. Methylation loss in users effected genes involved in immune system and inflammation, while methylation gain in nonusers effected genes involved in metabolism.
This is the first longitudinal study demonstrating effectiveness of aspirin-use in suppression of age-related and CC-relevant hypermethylation in the normal colon. These findings provide a rationale for future studies to evaluate loci that may serve as markers to identify individuals that will benefit most from aspirin and hence increase its efficiency in CC prevention and therapy.
Colon cancer (CC) is the third most common cancer worldwide and the second most common in terms of cancer-related mortality . Age is the strongest risk factor for sporadic CC ; however, the risk varies among individuals, even within the same age group, which may reflect heterogeneity in exposure to environmental factor. Indeed, several observational studies have shown that diet and lifestyle are associated with increase CC risk [2,3,4]. Thus, primary prevention by modification of environmental and lifestyle exposures and the use of chemoprevention are important strategies to reduce further CC incidence. Among chemopreventive agents, long-term use of aspirin is associated with decreased risk of CC incidence, mortality and cancer recurrence after therapy [5,6,7]. Understanding aspirin’s mechanism of action is important to increase its benefit-risk ratio, given the known hazards of its use [8,9,10]. While the inhibition of cyclooxygenases (COX) is well characterized as an underlying cause , several COX-independent mechanisms of aspirin, such as suppression of cellular proliferation and induction of apoptosis, inhibition of WNT/β-catenine and NF-kB signaling have been suggested to play a role [12,13,14]; however, the impact of aspirin use on epigenome as potential mechanism is poorly studied.
Epigenetic modifications, particularly DNA methylation (DNAm) changes, are recognized as one of the most common molecular alterations in human tumors, including CC . CC-subtype-specific DNAm are readily detectable in healthy normal mucosa [16,17,18], before the occurrence of visible precursor lesions, which ultimately may progress to CC . The stability of DNAm on as molecular marker and its effect on gene expression  facilitates its clinical use in early cancer detection and makes it a potential target to predict treatment outcome and patient’s response to therapy.
We have previously shown in a cross-sectional study that aspirin modulates age-associated DNAm changes in the colon and, thereby, CC risk . Understanding how aspirin controls the stability of DNAm in aging tissue is of both biological and clinical importance. Hence, the present study investigated the effect of aspirin use on DNAm in a longitudinal screening female cohort and its effect on human biological aging in normal colon using three established epigenetic clocks; Horvath’s , Hannum’s  and PhenoAge . Finally, we assessed the functional link between suppressive gene methylation by aspirin use and their link to CC risk.
To assess global changes in the DNAm over time, we subjected 124 healthy normal mucosa samples of screening females to Illumina infinium MethylationEPIC BeadChip. The major clinical characteristics of the participants are summarized in Tables 1 and 2. Principal component analysis (PCA) on methylation profiles separated proximal from distal colon samples (28% variation), aspirin users (U) from nonusers (Nu) and/or t1 from t2 (8% variation) (Fig. 1a). Multivariable linear regression model was used to measure significant differentially methylated CpGs in aspirin users (U-dmCpGs) and in nonusers (Nu-dmCpGs) from t1 to t2. We identified more hypomethylated U-dmCpGs in aspirin users and more hypermethylated Nu-dmCpGs in nonusers (Fig. 1b, Additional file 1: Figure S1a). At U-dmCpGs, median methylation change (hyper and hypo) was generally similar in users and nonusers, while at Nu-dmCpGs, it was lower in users than nonusers (Fig. 1c). A methylation index (MI), designating the percentage of mean methylation across all dmCpGs, was calculated for each individual. In proximal colon, over time, the MI was increased in nonusers (P = 0.01), but not in users, which at t2 showed 2% (P = 0.02) decreased MI in users vs nonusers. In distal colon, the MI was decreased independently of aspirin use (Fig. 1d).
Enrichment analyses showed that in both colon locations dmCpGs were generally under-represented at CpG-island (CGI) promoters. Hypermethylated U-dmCpGs, however, were enriched at promoters regardless of the CGI status and hypomethylated U-dmCpGs were enriched at weakly transcribed regions according to their chromatin state defined by ChromHMM (Fig. 2a). The functional significance of genes associated with dmCpGs was investigated by overlapping them with tumor suppressor genes (TSG) and oncogenes known to be aberrantly expressed in CC [24, 25]. Notably, genes associated with hypomethylated U-dmCpGs and hypermethylated Nu-dmCpGs involved several TSG (Additional file 1: Figure S1b). Among them were genes that regulate DNA and histone methylation, e.g., (KDM8, SIRT1/2, CHD5). These findings suggest a widespread modulatory effect of aspirin on the colonic epigenome (Fig. 2b). Genes associated with U-dmCpGs (hyper and hypo) were allocated in 13 pathways involved in actin cytoskeleton, immune system (estrogen signaling) and inflammation (PI3K-Akt and cAMP signaling). Genes associated with Nu-dmCpGs (hyper and hypo) were enriched in 10 pathways primarily involved in metabolism (Fig. 2c). Genes associating with both U- and Nu-dmCpGs were enriched in 7 pathways, including platelet activation and chemokines signaling pathways (Fig. 2c).
Next, we asked whether the long-term aspirin use has an impact on methylation age (mAge). We utilized previously established three mAge estimators; Hannum (71 CpGs), Horvath (353 CpGs) and PhenoAge (513 CpGs). Overall, the median methylation change on CpGs corresponding to three mAge estimators was significantly reduced in users (Hannum: P = 0.0001; Horvath: P = 0.0001; PhenoAge: P = 3e−06) versus nonusers in proximal colon (Fig. 3a). No difference was found in distal colon. Further, there was moderate to high correlation between chronological age and mAge for all three estimators (r = 0.6–0.8); however, the Horvath based mAge was substantially lower than chronological age (Fig. 3b). To examine if the biological age varies between aspirin users and nonusers, we used age acceleration residuals (AAR). The residuals were obtained from regressing mAge on chronological age of all 124 samples; hence, they adjust for the chronological age bias between aspirin users and nonusers. Over time, in aspirin users, we found significant deceleration of Horvath mAge (P = 0.05) and PhenoAge (P = 0.02), and the deceleration trend was the same for Hannum age (P = 0.056) (Fig. 3c). Taken together, over time, aspirin users show low methylation levels and display decelerated mAge in proximal colon.
We conducted a 10-year follow-up study of our previous female screening cohort, finding the evidence that long-term aspirin use suppresses DNAm even after 10 years of its initial use and is associated with decelerated biological aging and CC risk.
The repressed DNAm was not limited to colon but was also found on aging-associated CpGs derived from multi-tissue and blood, indicating aspirin’s beneficial role beyond tissue type. We recognize that there are baseline DNAm differences between aspirin users and nonusers at t1 that may reflect different age and/or disease risks and reasons for aspirin use; however, baseline DNAm levels and age was corrected in our analysis. The consistency of our findings across Horvath , Hannum  and PhenoAge  defined age-associated CpGs implies that the observed suppression of DNAm in aspirin users is unlikely an artifact of the current data set. Our findings are in accordance with previous studies suggesting that epigenetic modifications are involved in the chemopreventive actions of aspirin. Recently, an association between aspirin use and BRCA1 promoter methylation that affects the mortality after breast cancer has been shown . Similarly, in mouse inflammation model, aspirin has been shown to reduced tumor multiplicity by suppressing histone deacetylases (HDAC), with concomitant increase of H3K27ac in the promoter regions of pro-inflammatory genes .
Previously in a cross-sectional screening female’s cohort, we found that aspirin represses rate of DNAm change in both colon locations. However, in the present longitudinal study, we found overall repressed MI in proximal but not in distal colon. The location-specific effect of aspirin was consistent with the repressed mAge in users in proximal colon for all three epigenetic clocks, but not in distal colon. These results are in line with previous reports showing daily aspirin intake reduces CC incidence and mortality and have a greater effect on CC of the proximal colon than on the distal colon . Further evaluation for colon location dependence is needed with large sample size.
We found decelerated mAge using PhenoAge clock in aspirin users, although the same trend was found for Hannum’s and Horvath’s clock as well. While the accelerated mAge has been identified as a potential risk factor for age-related diseases and all-cause mortality [21, 23], the decelerated aging is expected to reduce disease incidence at any given age. The effect of aspirin on decelerated mAge in the present study is consistent with previous report showing increased lifespan of male mice upon aspirin use . The mechanism(s) that drive deviation of mAge from chronological age are independent of cell proliferation, as these clocks operate in nonproliferation tissue (e.g., brain) as well as proliferating tissue (e.g., peripheral blood leukocytes) [21, 22].
There are several limitations to the current study. First is its small sample size. Second is the observational design of the study cannot imply that observed suppression in DNAm was caused by aspirin use. However, we adjusted for the potential confounders and still have large effects. Third is one important but unresolved question that is why aspirin users have more polyps in total versus nonusers. Interestingly, most of the polyps developed from t1 to t2 in the distal but not in the proximal colon, which is in line with our molecular data showing that the protective effect of aspirin was found predominantly in the proximal colon. This relationship and general significance for broader population needs to be tested in follow-up studies.
Taken together, we show that methylation age in aspirin nonusers increases with increase in chronological age; however, in aspirin users, methylation age acceleration is suppressed (Fig. 4). Our data show that many CC-specific chromatin modulators gain methylation in nonusers over time (such as CHD5, PADI2, KDM8) but lose methylation in aspirin users, suggesting that the aspirin may reverse the hypermethylation of epigenetic modifiers and TSGs in colon. Investigating the effect of aspirin as an epigenetic modulator may help to better characterize the functional mechanisms underlying its anticancer activity.
In conclusion, this is the first longitudinal study demonstrating a suppression of DNAm and epigenetic aging in healthy colon over time by aspirin use. While the study cohort is small and the observation period is short (10 years) in the context of a human lifespan, this pilot study has a strength of its longitudinal design. These findings provide a rationale for future longitudinal studies with large sample size to improve our understanding on epigenetic response to aspirin, optimize its potential benefit in the suppression of age-associated DNAm increase, CC risk, and hence increase its efficiency in CC prevention and therapy.
Study participants and data collection
We contacted 89 healthy females included in our previous study . Colonoscopy could be performed in 31 females with 58 subjects being lost for various reasons (death, no response to invitation, discourage by family doctor for colonoscopy). Total 124 normal mucosa biopsies from cecum (proximal) and sigmoid (distal) colon were obtained at baseline (t1) and after 10 years (t2; range: 10–10.6 years) from 31 screening females (four biopsies per subject). Information on aspirin use was obtained using a self-administered questionnaire with assistance of a study nurse. Aspirin use was determined as never versus (vs.) long-term use (≥ twice per week for ≥ 2 years) at either the baseline (t1) or follow-up (t2). Of the 31 individuals, 17 were aspirin users at baseline and 20 at follow-up (three subjects switched from nonuser at t1 to user at t2). The term user was not restricted to the use of aspirin period immediately prior to the sampling at both time points. Height (cm) and weight (kg) were self-reported, and BMI was calculated (kg/m2). Polyp status was determined as no versus yes (any type) at baseline (t1) and at follow-up (t2). Polyp type was classified as tubular adenoma versus serrated lesions according to WHO classification .
DNA isolation and bisulfite conversion
All endoscopies were performed by using high-definition colonoscopy including narrow-band imaging. Biopsies were taken from macroscopic normal mucosa without microscopically assessment of the collected specimens. Females with symptoms or history of inflammatory bowel disease were excluded from this study. Genomic DNA was extracted from freshly collected biopsies of normal mucosa stored in RNAlater using Quick-DNA Miniprep kit (Zymo Research) including RNAse A treatment according to the manufacturer protocol. Five hundred ng of genomic DNA from each sample was bisulfite-treated using the EZ DNA Methylation kit (Zymo Research, Orange, CA), according to the manufacturer's protocol and stored at − 80 °C until further use.
Genome-Wide DNA methylation profiling
Sodium bisulfite-treated DNA of 124 biopsies was subjected to measure global DNAm by Illumina infinium MethylationEPIC BeadChip . Samples were randomly distributed across 21 chips (8 samples per chip) to avoid confounding batch effects. Illumina’s GenomeStudio software V2011.1 was used to extract the raw signal intensities of 865,859 probes. Downstream analysis was conducted using annotation developed previously mapping to hg38 or otherwise hg19 as indicated. Raw IDAT files were preprocessed using R/Bioconductor package minfi (version 1.32.0) [32, 33]. Samples were normalized using “preprocessFunnorm” including background subtraction and dye-bias normalization. Methylation levels were computed as Beta (β) values (Meth/[Meth + Unmeth + 100]) and as logit-transformed M value (log2 [Meth/Unmeth]). To avoid extreme M values, we use beta threshold of 1e-07 to make sure beta was neither 0 nor 1, before taken the logit. Probes where at least one sample had a detection P value > 0.01 and the probes recommended previously  using hg38 genome as general-purpose masking probes were filtered out. After filtering, 753,630 probes were kept for further analysis. M value were used to identify differentially methylated CpGs (dmCpGs) using bioconductor limma  package.
To derive an estimate of DNA methylation Age (mAge), we used three different algorithms: the blood based Hannum clock , multi-tissue based Horvath clock  and blood based PhenoAge  that differentiate morbidity and mortality risk among same-age individuals were used to determine epigenetic aging. Of the requisite 71 CpGs in Hannum-CpGs, 6 are missing from the EPIC and 4 were filtered out at QC. Of the requisite 353 CpGs in Horvath-CpGs, 17 are missing from EPIC and the 10 were filtered out during data quality control (QC). Of the requisite 513 CpGs in PhenoAge clock, 18 were filtered out during data QC. The methylation levels on the remaining 61 of Hannum-CpGs, 326 Horvath-CpGs and 495 PhenoAge CpGs were utilized to compare median methylation differences in Fig. 3a. The mAge and age acceleration residuals (AAR) for three algorithms were calculated using the software (https://dnamage.genetics.ucla.edu/new). The raw beta values without masking probes based on a detection P value were submitted. This resulted in submission of 64 of Hannum-CpGs, 336 Horvath-CpGs and 513 PhenoAge CpGs. The age calculator automatically imputes missing data across platform based on k-nearest neighbor imputation.
Pathway analysis and annotation
Pathway analysis was done using Kyoto Encyclopedia of Genes and Genomes (KEGG). Gene annotation was obtained from IlluminaHumanMethylationEPICanno.ilm10b4.hg19, and number of CpGs per gene was controlled. For enrichment, all genes present on BeadChip were used as background and pathways < 0.05 adjusted P value were considered. For functional annotation, TSGene , ONGene  and EpiFactors  database were used. Detailed annotation of probes on gene regulatory elements (promoters, enhancers and weak transcripts) were done using ENCODE  as described previously . The CpG Island (CGI) were defined as in UCSC.
Normalized methylation levels were analyzed using multivariable linear regression model between t2 versus t1 adjusted for body mass index (BMI), polyps, age and batch-effect. The paired nature of samples was taken into account by adding subjects in the model matrix. First, model matrix was designed to compare aspirin users to nonusers for each time point and colon location, and then, contrasts were made to address the effect of aspirin in user and nonuser from t1 to t2. Differentially methylated CpGs (dmCpGs) were identified using empirical Bayes moderated t statistics and associated Benjamini–Hochberg (BH) false discovery rate (FDR) adjusted P values < 0.05 were used as cutoff. Shapiro–Wilk test was used for normality testing for Fig. 1d and 3a. For methylation index (MI) analysis in Fig. 1d, paired Wilcoxon signed rank test was used. Wilcoxon rank sum tests were used to determine significance of median methylation differences between groups at each epigenetic clock in Fig. 3a. For enrichment analysis, P values were calculated using Fisher’s test. All analyses were performed using R (version 4.0.2) and R-Studio (version 1.1.463).
Availability of data and materials
The DNA methylation data generated in this study are available through NCBI Gene Expression Omnibus (GEO) https://www.ncbi.nlm.nih.gov repository under the Accession Number GSE142257.
Differentially methylated CpGs in aspirin users
Differentially methylated CpGs in nonusers
Tumors suppressor genes
Time point 1
Time point 2
Keum N, Giovannucci E. Global burden of colorectal cancer: emerging trends, risk factors and prevention strategies. Nat Rev Gastroenterol Hepatol. 2019;16(12):713–32. https://doi.org/10.1038/s41575-019-0189-8.
Dekker E, Tanis PJ, Vleugels JLA, Kasi PM, Wallace MB. Colorectal cancer. Lancet. 2019;394(10207):1467–80. https://doi.org/10.1016/S0140-6736(19)32319-0.
Conti L, Del Corno M, Gessani S. Revisiting the impact of lifestyle on colorectal cancer risk in a gender perspective. Crit Rev Oncol Hematol. 2020;145:102834. https://doi.org/10.1016/j.critrevonc.2019.102834.
Aleksandrova K, Pischon T, Jenab M, Bueno-de-Mesquita HB, Fedirko V, Norat T, et al. Combined impact of healthy lifestyle factors on colorectal cancer: a large European cohort study. BMC Med. 2014;12:168. https://doi.org/10.1186/s12916-014-0168-4.
Katona BW, Weiss JM. Chemoprevention of colorectal cancer. Gastroenterology. 2019. https://doi.org/10.1053/j.gastro.2019.06.047.
Drew DA, Cao Y, Chan AT. Aspirin and colorectal cancer: the promise of precision chemoprevention. Nat Rev Cancer. 2016;16(3):173–86. https://doi.org/10.1038/nrc.2016.4.
Qiao Y, Yang T, Gan Y, Li W, Wang C, Gong Y, et al. Associations between aspirin use and the risk of cancers: a meta-analysis of observational studies. BMC Cancer. 2018;18(1):288. https://doi.org/10.1186/s12885-018-4156-5.
Huang ES, Strate LL, Ho WW, Lee SS, Chan AT. Long-term use of aspirin and the risk of gastrointestinal bleeding. Am J Med. 2011;124(5):426–33. https://doi.org/10.1016/j.amjmed.2010.12.022.
Thun MJ, Jacobs EJ, Patrono C. The role of aspirin in cancer prevention. Nat Rev Clin Oncol. 2012;9(5):259–67. https://doi.org/10.1038/nrclinonc.2011.199.
Avivi D, Moshkowitz M, Detering E, Arber N. The role of low-dose aspirin in the prevention of colorectal cancer. Expert Opin Ther Targets. 2012;16(Suppl 1):S51–62. https://doi.org/10.1517/14728222.2011.647810.
Alfonso L, Ai G, Spitale RC, Bhat GJ. Molecular targets of aspirin and cancer prevention. Br J Cancer. 2014;111(1):61–7. https://doi.org/10.1038/bjc.2014.271.
Bos CL, Kodach LL, van den Brink GR, Diks SH, van Santen MM, Richel DJ, et al. Effect of aspirin on the Wnt/beta-catenin pathway is mediated via protein phosphatase 2A. Oncogene. 2006;25(49):6447–56. https://doi.org/10.1038/sj.onc.1209658.
Yin MJ, Yamamoto Y, Gaynor RB. The anti-inflammatory agents aspirin and salicylate inhibit the activity of I(kappa)B kinase-beta. Nature. 1998;396(6706):77–80. https://doi.org/10.1038/23948.
Shimura T, Toden S, Komarova NL, Boland C, Wodarz D, Goel A. A comprehensive in vivo and mathematic modeling-based kinetic characterization for aspirin-induced chemoprevention in colorectal cancer. Carcinogenesis. 2020. https://doi.org/10.1093/carcin/bgz195.
Ehrlich M. DNA methylation in cancer: too much, but also too little. Oncogene. 2002;21(35):5400–13. https://doi.org/10.1038/sj.onc.1205651.
Menigatti M, Truninger K, Gebbers JO, Marbet U, Marra G, Schar P. Normal colorectal mucosa exhibits sex- and segment-specific susceptibility to DNA methylation at the hMLH1 and MGMT promoters. Oncogene. 2009;28(6):899–909. https://doi.org/10.1038/onc.2008.444.
Noreen F, Roosli M, Gaj P, Pietrzak J, Weis S, Urfer P, et al. Modulation of age- and cancer-associated DNA methylation change in the healthy colon by aspirin and lifestyle. J Natl Cancer Inst. 2014. https://doi.org/10.1093/jnci/dju161.
Wang T, Maden SK, Luebeck GE, Li CI, Newcomb PA, Ulrich CM, et al. Dysfunctional epigenetic aging of the normal colon and colorectal cancer risk. Clin Epigenetics. 2020;12(1):5. https://doi.org/10.1186/s13148-019-0801-3.
Noreen F, Kung T, Tornillo L, Parker H, Silva M, Weis S, et al. DNA methylation instability by BRAF-mediated TET silencing and lifestyle-exposure divides colon cancer pathways. Clin Epigenetics. 2019;11(1):196. https://doi.org/10.1186/s13148-019-0791-1.
Baylin SB, Jones PA. Epigenetic determinants of cancer. Cold Spring Harb Perspect Biol. 2016. https://doi.org/10.1101/cshperspect.a019505.
Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013;14(10):R115. https://doi.org/10.1186/gb-2013-14-10-r115.
Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013;49(2):359–67. https://doi.org/10.1016/j.molcel.2012.10.016.
Levine ME, Lu AT, Quach A, Chen BH, Assimes TL, Bandinelli S, et al. An epigenetic biomarker of aging for lifespan and healthspan. Aging (Albany NY). 2018;10(4):573–91. https://doi.org/10.18632/aging.101414.
Liu Y, Sun J, Zhao M. ONGene: a literature-based database for human oncogenes. J Genet Genomics. 2017;44(2):119–21. https://doi.org/10.1016/j.jgg.2016.12.004.
Zhao M, Sun J, Zhao Z. TSGene: a web resource for tumor suppressor genes. Nucleic Acids Res. 2013;41(Database issue):D970–6. https://doi.org/10.1093/nar/gks937.
Wang T, McCullough LE, White AJ, Bradshaw PT, Xu X, Cho YH, et al. Prediagnosis aspirin use, DNA methylation, and mortality after breast cancer: a population-based study. Cancer. 2019;125(21):3836–44. https://doi.org/10.1002/cncr.32364.
Guo Y, Liu Y, Zhang C, Su ZY, Li W, Huang MT, et al. The epigenetic effects of aspirin: the modification of histone H3 lysine 27 acetylation in the prevention of colon carcinogenesis in azoxymethane- and dextran sulfate sodium-treated CF-1 mice. Carcinogenesis. 2016;37(6):616–24. https://doi.org/10.1093/carcin/bgw042.
Rothwell PM, Wilson M, Elwin CE, Norrving B, Algra A, Warlow CP, et al. Long-term effect of aspirin on colorectal cancer incidence and mortality: 20-year follow-up of five randomised trials. Lancet. 2010;376(9754):1741–50. https://doi.org/10.1016/S0140-6736(10)61543-7.
Strong R, Miller RA, Astle CM, Floyd RA, Flurkey K, Hensley KL, et al. Nordihydroguaiaretic acid and aspirin increase lifespan of genetically heterogeneous male mice. Aging Cell. 2008;7(5):641–50. https://doi.org/10.1111/j.1474-9726.2008.00414.x.
Nagtegaal ID, Odze RD, Klimstra D, Paradis V, Rugge M, Schirmacher P, et al. The 2019 WHO classification of tumours of the digestive system. Histopathology. 2019. https://doi.org/10.1111/his.13975.
Pidsley R, Zotenko E, Peters TJ, Lawrence MG, Risbridger GP, Molloy P, et al. Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling. Genome Biol. 2016;17(1):208. https://doi.org/10.1186/s13059-016-1066-1.
Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, et al. Minfi: a flexible and comprehensive bioconductor package for the analysis of infinium DNA methylation microarrays. Bioinformatics. 2014;30(10):1363–9. https://doi.org/10.1093/bioinformatics/btu049.
Fortin JP, Triche TJ Jr, Hansen KD. Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi. Bioinformatics. 2017;33(4):558–60. https://doi.org/10.1093/bioinformatics/btw691.
Zhou W, Laird PW, Shen H. Comprehensive characterization, annotation and innovative use of Infinium DNA methylation BeadChip probes. Nucleic Acids Res. 2017;45(4):e22. https://doi.org/10.1093/nar/gkw967.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47. https://doi.org/10.1093/nar/gkv007.
Medvedeva YA, Lennartsson A, Ehsani R, Kulakovskiy IV, Vorontsov IE, Panahandeh P, et al. EpiFactors: a comprehensive database of human epigenetic factors and complexes. Database (Oxford). 2015;2015:bav067. https://doi.org/10.1093/database/bav067.
Hoffman MM, Ernst J, Wilder SP, Kundaje A, Harris RS, Libbrecht M, et al. Integrative annotation of chromatin elements from ENCODE data. Nucleic Acids Res. 2013;41(2):827–41. https://doi.org/10.1093/nar/gks1284.
We thank Robert Ivanek for his support and critical reviewing of bioinformatics analyses.
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
Ethics approval and consent to participate
The study was approved by the Ethics Committee of the Medical Center for Postgraduate Education, Warsaw, Poland (No. 25/2006/1/2015).
Consent for publication
All authors give consent for the publication of the manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. a) Overlap of all significant differentially methylated CpGs found in aspirin users and nonusers. The non-overlapping CpGs were then defined as aspirin user specific differentially methylated CpGs (U-dmCpGs) or nonuser specific differentially methylated CpGs (Nu-dmCpGs). b) Overlap of genes affected by Nu-dmCpGs and U-dmCpGs with known differential expressed colon cancer specific tumor suppressor genes (TSGs) and oncogenes. Shown are separate overlaps for proximal and distal colon.
About this article
Cite this article
Noreen, F., Chaber-Ciopinska, A., Regula, J. et al. Longitudinal analysis of healthy colon establishes aspirin as a suppressor of cancer-related epigenetic aging. Clin Epigenet 12, 164 (2020). https://doi.org/10.1186/s13148-020-00956-9
- Colon cancer
- DNA methylation