Galamb O, Kalmár A, Barták BK, Patai &V, Leiszter K, Péterfia B, Wichmann B, Valcz G, Veres G, Tulassay Z, Molnár B. Aging related methylation influences the gene expression of key control genes in colorectal cancer and adenoma. World J Gastroenterol 2016; 22(47): 10325-10340 [PMID: 28058013 DOI: 10.3748/wjg.v22.i47.10325]
Corresponding Author of This Article
Orsolya Galamb, PhD, 2nd Department of Internal Medicine, Semmelweis University, Szentkirályi str 46, H-1088 Budapest, Hungary. orsg1@yahoo.com
Research Domain of This Article
Gastroenterology & Hepatology
Article-Type of This Article
Basic Study
Open-Access Policy of This Article
This article is an open-access article which was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Orsolya Galamb, Barnabás Wichmann, Gábor Valcz, Zsolt Tulassay, Béla Molnár, Molecular Medicine Research Group, Hungarian Academy of Sciences, H-1088 Budapest, Hungary
Orsolya Galamb, Alexandra Kalmár, Barbara Kinga Barták, Árpád V Patai, Katalin Leiszter, Bálint Péterfia, 2nd Department of Internal Medicine, Semmelweis University, H-1088 Budapest, Hungary
Gábor Veres, 1st Department of Paediatrics, Semmelweis University, H-1083 Budapest, Hungary
ORCID number: $[AuthorORCIDs]
Author contributions: Galamb O, Kalmár A, Péterfia B and Molnár B designed the study; Patai ÁV, Veres G and Molnár B collected the samples; Galamb O, Patai ÁV, Leiszter K, Valcz G and Veres G contributed to the collection of clinical data and histological analysis of the samples; Galamb O, Kalmár A, Barták BK and Patai ÁV performed the experiments; Galamb O, Kalmár A, Wichmann B and Valcz G analyzed the experimental data; Tulassay Z and Molnár B contributed to the design and critical review of the manuscript, obtained fundings; all authors were involved in writing the paper, made a critical revision of the manuscript for important intellectual content and had final approval of the submitted and published versions.
Supported bythe National Research, Development and Innovation Office, No. KMR-12-1-2012-0216; and the Hungarian Scientific Research Fund, No. OTKA-K111743.
Institutional review board statement: All routine colonic biopsy samples from the patients were taken after informed consent and ethical permission was obtained for participation in the study.
Conflict-of-interest statement: The authors declare that no conflict of interest exists.
Data sharing statement: Additional data are available in a supplementary file.
Open-Access: This article is an open-access article which was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Correspondence to: Orsolya Galamb, PhD, 2nd Department of Internal Medicine, Semmelweis University, Szentkirályi str 46, H-1088 Budapest, Hungary. orsg1@yahoo.com
Telephone: +36-1-2660926 Fax: +36-1-2660816
Received: June 22, 2016 Peer-review started: June 24, 2016 First decision: August 22, 2016 Revised: September 20, 2016 Accepted: November 13, 2016 Article in press: November 13, 2016 Published online: December 21, 2016 Processing time: 180 Days and 3 Hours
Abstract
AIM
To analyze colorectal carcinogenesis and age-related DNA methylation alterations of gene sequences associated with epigenetic clock CpG sites.
METHODS
In silico DNA methylation analysis of 353 epigenetic clock CpG sites published by Steve Horvath was performed using methylation array data for a set of 123 colonic tissue samples [64 colorectal cancer (CRC), 42 adenoma, 17 normal; GEO accession number: GSE48684]. Among the differentially methylated age-related genes, secreted frizzled related protein 1 (SFRP1) promoter methylation was further investigated in colonic tissue from 8 healthy adults, 19 normal children, 20 adenoma and 8 CRC patients using bisulfite-specific PCR followed by methylation-specific high resolution melting (MS-HRM) analysis. mRNA expression of age-related “epigenetic clock” genes was studied using Affymetrix HGU133 Plus2.0 whole transcriptome data of 153 colonic biopsy samples (49 healthy adult, 49 adenoma, 49 CRC, 6 healthy children) (GEO accession numbers: GSE37364, GSE10714, GSE4183, GSE37267). Whole promoter methylation analysis of genes showing inverse DNA methylation-gene expression data was performed on 30 colonic samples using methyl capture sequencing.
RESULTS
Fifty-seven age-related CpG sites including hypermethylated PPP1R16B, SFRP1, SYNE1 and hypomethylated MGP, PIPOX were differentially methylated between CRC and normal tissues (P < 0.05, Δβ≥ 10%). In the adenoma vs normal comparison, 70 CpG sites differed significantly, including hypermethylated DKK3, SDC2, SFRP1, SYNE1 and hypomethylated CEMIP, SPATA18 (P < 0.05, Δβ≥ 10%). In MS-HRM analysis, the SFRP1 promoter region was significantly hypermethylated in CRC (55.0% ± 8.4 %) and adenoma tissue samples (49.9% ± 18.1%) compared to normal adult (5.2% ± 2.7%) and young (2.2% ± 0.7%) colonic tissue (P < 0.0001). DNA methylation of SFRP1 promoter was slightly, but significantly increased in healthy adults compared to normal young samples (P < 0.02). This correlated with significantly increased SFRP1 mRNA levels in children compared to normal adult samples (P < 0.05). In CRC tissue the mRNA expression of 117 age-related genes were changed, while in adenoma samples 102 genes showed differential expression compared with normal colonic tissue (P < 0.05, logFC > 0.5). The change of expression for several genes including SYNE1, CLEC3B, LTBP3 and SFRP1, followed the same pattern in aging and carcinogenesis, though not for all genes (e.g., MGP).
CONCLUSION
Several age-related DNA methylation alterations can be observed during CRC development and progression affecting the mRNA expression of certain CRC- and adenoma-related key control genes.
Core tip: Several age-related DNA methylation alterations could be observed during colorectal cancer (CRC) formation and progression affecting the mRNA expression of certain CRC- and adenoma-related key control genes such as hypermethylated secreted frizzled related protein 1 (SFRP1), spectrin repeat containing nuclear envelope protein 1 and hypomethylated cell migration-inducing protein. For the first time significantly lower SFRP1 methylation levels were demonstrated in colonic tissue from children (under 18 years) compared to healthy adults. The main CRC-associated signal transduction pathways, such as WNT signaling and PI3K/Akt pathways are also influenced during aging.
Citation: Galamb O, Kalmár A, Barták BK, Patai &V, Leiszter K, Péterfia B, Wichmann B, Valcz G, Veres G, Tulassay Z, Molnár B. Aging related methylation influences the gene expression of key control genes in colorectal cancer and adenoma. World J Gastroenterol 2016; 22(47): 10325-10340
DNA methylation alterations in connection with aging include epigenetic drift and epigenetic clock phenomena. Epigenetic drift is defined as the global DNA methylation changes caused by random and environmental individual-specific factors, while the epigenetic clock is defined as a group of progressive age-related epigenetic alterations at specific genomic sites which are common across individuals and occassionally across various tissue types[1,2]. The epigenetic clock concept is an approach to biological age prediction of different tissues based on the DNA methylation status of 353 CpG sites measured using the Illumina Beadchip450K methylation array platform[2].
Although age-related (A type) and cancer-related (C type) DNA methylation are often distinguished, the main age-related disease is cancer and the age of patients is one of the risk factor for carcinogenesis[3]. In human development, following a transient increase in average DNA methylation in early childhood (during the first year of life)[4,5], global hypomethylation is characteristic during aging[6,7]. Similarly global hypomethylation is observed in various types of cancers including colorectal cancer (CRC)[8]. With aging, besides global hypomethylation, local hypermethylation can occur on promoters of certain genes, including tumor suppressor gene promoters in various types of cancers, and many tumor suppressor genes were reported among the age-dependently hypermethylated genes[6]. Among others, promoter hypermethylation of APC[7,9-12], CDKN2A[7,9], ESR1[7,13,14], GATA5[15,16], HPP1[7,15,17], SFRP1[7,12,15,18-23] and SFRP2[7,18-21] genes was reported for colonic tissues during both aging and colorectal carcinogenesis. Although DNA methylation data from adult colonic tissue samples has been expansively published, data for children/young patients are limited.
In this study, we analyzed DNA methylation and/or gene expression changes of genes covered by the 353 epigenetic clock CpG sites[2] for patients of different ages as well as stages in the progression through to CRC in order to study the possible relationship between age-related and cancer-associated epigenetic alterations. Gene expression analysis was performed using colonic tissue samples from healthy children, healthy adults, and patients with adenomas and CRC. Among the differentially methylated/expressed age-related genes, secreted frizzled related protein 1 (SFRP1) promoter methylation was further analyzed in healthy, premalignant and cancerous colonic tissue samples, and to our knowledge this is the first study to also include colonic biopsy specimens from children.
MATERIALS AND METHODS
In silico DNA methylation analysis
The DNA methylation status of 353 age-related CpG sites[2] was analyzed in silico using methylation array data from the Illumina BeadChip450K. Analysis was performed on 123 CRC, adenoma and normal tissue samples available in the NCBI Gene Expression Database database (GEO accession number: GSE48684[24]). Differences between average methylation values of the compared diagnostic groups (Δβ-values) and P values were determined for each CpG site (cg IDs). For statistical evaluation, normal distribution was checked using Kolmogorov-Smirnov test. Hence normal distribution was observed in any cases, Student’s t-test with Benjamini and Hochberg correction was applied for paired group comparisons. Significance criteries were P < 0.05 in all cases.
In silico gene expression analysis
The expression of age-related “epigenetic clock” genes was analyzed using whole transcriptome data from Affymetrix HGU133 Plus2.0. Data was obtained from 153 colonic biopsy samples (49 healthy, 49 adenoma, 49 CRC and 6 healthy children) previously hybridized by our research group (GEO serial accession numbers: GSE37364[25], GSE10714[26], GSE4183[27], GSE37267[28]). Gene expression levels were compared using unpaired Student’s t-test with Benjamini and Hochberg correction (P value of < 0.05 was considered as significant). For gene expression analysis, normal distribution was found using Kolmogorov-Smirnov test, therefore Student’s t-test (in case of differentiation of two groups with equal variances) or Welch-test (in case of differentiation of two groups with unequal variances) and ANOVA (when more than two groups were compared) were applied. For paired comparisons Benjamini and Hochberg correction was applied. In case of ANOVA, Tukey HSD post-test was used in order to find out which group refers to the differentiation if any. Significance criteries were P < 0.05 in any cases. For the logFC calculation, the differences between the averages of groups were considered (abs logFC ≥ 0.5 criteria).
Methyl capture sequencing - in silico data analysis
Whole methylome data from 6 normal adjacent tissue (NAT), 15 adenoma and 9 CRC tissue samples were determined in a previous study using methyl capture sequencing[12]. Using this dataset, the whole promoter methylation status of genes showing an inverse relation between gene expression and DNA methylation was evaluated. Differentially methylated genes were determined as described earlier[12]. For statistical evaluation normal distribution was determined and the applied tests were chosen according to the above-mentioned criteria. Differences with P < 0.05 were considered as significant. Methylation alterations between diagnostic groups were characterized by Δβ-values (the differences of the average β-values of sample groups).
Clinical samples
All patients provided informed consent. Colorectal biopsy samples were obtained during routine endoscopic intervention at the 2nd Department of Internal Medicine and 1st Department of Paediatrics, Semmelweis University, Budapest, Hungary. In total 55 colonic tissue samples (from 19 healthy children (under age of 18 years), 8 healthy adults, 20 patients with adenomas and 8 CRC samples) were tested in SFRP1 methylation-specific high resolution melting (MS-HRM) study (Table 1). Biopsy samples from all adults and 5 of children were stored in RNALater Stabilization Solution (Ambion, ThermoFisher Scientific) at -80 °C until use. Biopsy samples from the same site were immediately fixed in buffered formalin for histological evaluation. For 14 children, only FFPE blocks were available. Histological diagnoses were established by experienced pathologists. Altogether 27 tissue samples (19 from children and 8 from adults) with normal histology (so called healthy normal colonic tissue samples) were involved in SFRP1 MS-HRM study. Children and adults in the study had been referred to the outpatient clinic with rectal bleeding, constipation or chronic abdominal pain. Ileocolonoscopy was part of their diagnostic procedure (exclude organic disease) and the biopsy specimens showed normal macroscopic appearance and histology[28]. The study was conducted according to the Helsinki declaration and approved by the local ethics committee and government authorities (Regional and Institutional Committee of Science and Research Ethics (TUKEB) Nr.: 69/2008, 202/2009 and 23970/2011 Semmelweis University, Budapest, Hungary).
Table 1 Clinical data of samples involved in the high resolution melting study.
Tissue samples were homogenized in 2% sodium dodecyl sulphate, and digested with 4 mg/mL proteinase K for 16 h at 56 °C. Genomic DNA was isolated using the High Pure PCR Template Preparation Kit (Roche Applied Science) according to the manufacturer’s instructions[18]. DNA was eluted in 2 × 100 μL RNase- and DNase-free water and stored at -20 °C. The quantity of the isolated DNA samples was measured by Qubit fluorometer using the Qubit dsDNA HS Assay (Invitrogen, ThermoFisher Scientific).
Bisulfite-specific PCR and MS-HRM experiments for SFRP1 promoter methylation analysis
Bisulfite conversion was performed using the EZ DNA Methylation Direct™ Kit (Zymo Research, Irvine, CA, United States) according to the manufacturer’s instructions. For fresh frozen samples, 500 ng isolated DNA was converted, while for FFPE samples the total recovered DNA after deparaffinization and digestion was bisulfite converted. Bisulfite-specific PCR (BS-PCR) reactions were performed in 15 μL volume using LightCycler 480 Probes Master (Roche Applied Science), LightCycler® 480 ResoLight Dye (Roche), SFRP1 primers at 0.2 μmol/L final concentrations and bisulfite converted DNA (bcDNA) samples (approx, 5 ng bcDNA/well). The sequences of the applied SFRP1 BS-PCR primers were previously described[19]. Real-time PCR amplification was carried out on LightCycler 480 System with the following thermocycling conditions: 95 °C for 10 min, then 95 °C for 30 s, 62 °C with 0.4 °C decreasement/cycle for 30 s, 72 °C for 30 s for 10 touchdown cycles, followed by the amplification at 95 °C for 30 s, 58 °C for 30 s, and 72 °C for 30 s in 50 cycles. For HRM calibration, unmethylated and methylated bisulfite converted control DNA (EpiTect PCR Control DNA Set, Qiagen) were used in different ratios (0%, 10%, 25%, 50%, 75% and 100% methylated controls). HRM analyis was performed according to the following thermal conditions: after denaturation at 95 °C for 1 min, and cooling at 40 °C for 1 min, the reactions were continuously warmed up to 95 °C with a 25 acquisition/°C rate. Raw HRM data were exported and the HRM peak heights of the negative derivative of fluorescence over temperature curves (-(d/dT) Fluorescence) of the biological samples were retrieved at the melting temperatures of the methylated and unmethylated standards. The methylation percent was calculated by the ratio of values at the methylated and unmethylated melting temperatures. For statistical evaluation normal distribution was determined and the applied tests were chosen according to the above-mentioned criteria.
SFRP1 immunohistochemistry
Parallel with our epigenetic examinations, 4 nm thick FFPE samples from healthy children (n = 6) and healthy adults (n = 7) were examined. SFRP1 immunohistochemistry was performed on colonic tissue samples of healthy adults (n = 7; mean age at histology examination: 48 ± 17 years; 5 f/2 m) and of children (mean age: 12 ± 6 years); 3 f/3 m). Histology was diagnosed by an expert pathologist on routinely stained HE sections. Following deparaffinization and rehydration, microwave based antigen retrieval was performed in TRIS EDTA buffer (pH 9.0) (900 W/10 min, then 340 W/40 min). Samples were immunostained with SFRP1 polyclonal antibody (ab4193, Abcam, Cambridge, United Kingdom, 1:80 dilution) with diaminobenzidine - hydrogen peroxidase - chromogen-substrate system (cat#30014.K, HISTOLS-DAB, Histopathology Ltd., Hungary) and were digitalized by Pannoramic 250 Flash II scanner (with Zeiss Plan-Apochromat 20 × objective; 3DHISTECH Ltd, Hungary). Digital slides were semi-quantitatively analyzed with Pannoramic Viewer (ver.:1.15.3; 3DHISTECH) based on Q-score method (scored by multiplying the percentage of positive cells (P) by the intensity (I: +3, +2, +1, 0). Formula: Q = P × I; Maximum = 300). Epithelial and stromal compartments were examined separately, then these scores were summarized (Σ) (Σ Q-score maximum: 600) in order to have comparable results with our whole biopsy methylation analyses.
Statistical analysis
The applied statistical methods are outlined above after the descriptions of molecular and in silico analyses. The statistical review of the study was performed by a biomedical statistician.
RESULTS
Gene ontology of 353 CpG sites of the "epigenetic clock"
The “epigenetic clock” signature includes 353 CpG sites[2] including different genes, gene promoters and other genomic regions such as enhancers, insulators, Polycomb-repressed regions. From the above 353 CpG sites, DNA methylation levels of 193 were positively and of 160 were negatively correlated with age[2]. First we updated the annotation of the CpG sites and assigned official gene symbols according to the newest version of NCBI Gene Database. Approximately 80% of the genes belonging to 353 CpG sites could be classified into functional groups including highly represented transcriptional regulation, translation (15.93%), metabolism (12.36%), development and ontogenesis (8.24%) and transport (8.24%). Approximately 20% of the genes had unknown function. According to the Encode ChromHMM results of nine human cell lines, the majority of the 353 CpG sites (76.2 %) were located in the promoter regions of genes. Also, 57.2% of them were categorized as active promoters (declared if it was found “active promoter” in at least one of the nine analyzed cell lines), while 19% were found to be “weak promoter” (“weak promoter” in at least one of the nine analyzed cell lines). One fourth (23.8%) of the CpG sites was located in non-promoter regions such as enhancers, insulators, transcribed and repressed regions (Supplementary Table 1).
In silico DNA methylation analysis
Analysis of the Illumnina Beadchip450K methylation array data set of Luo et al[24] showed 137 (38.8%) of the epigenetic clock CpG sites were found to be significantly differentially methylated between CRC and normal tissue samples (P < 0.05). Approximately two third of these CpG sites had similar methylation changes in CRC samples as during aging, while one third of these CpG sites showed opposite alterations in CRC tissue as during aging (Supplementary Table 1). Among these, 57 CpG sites showed at least a 10% methylation difference: from the 57 CpG sites 45 were hypermethylated (including ADHFE1/cg08090772, MCAM/cg21096399, DKK3/cg13216057, SFRP1/cg02388150, SYNE1/cg26620959), while 12 CpG sites showed hypomethylation (such as MGP/cg00431549, PIPOX/cg06144905, STRA6/cg00075967, ERG/cg17274064) in CRCs (P < 0.05, Δβ≥ 10%) (Table 2, Supplementary Table 1). In the adenoma vs normal comparison, DNA methylation of 165 CpG sites (46.7%) were significantly altered from which 70 CpG sites showed a ≥ 10% methylation difference: 36 CpG sites were hypermethylated (e.g., SDC2/cg25070637, SFRP1/cg02388150, SYNE1/cg26620959) and 34 showed decreased methylation levels (including CEMIP/cg20828084, SPATA18/cg03103192, STRA6/cg00075967) in adenoma samples (P < 0.05, Δβ≥ 10%) (Table 3). In CRC samples 33 CpG sites were found to be hypermethylated (such as KRT20/cg00091693, STRA6/cg00075967, UROS/cg19346193) and only one (LTBP3/cg08965235) was hypomethylated compared to adenomas (P < 0.05, Δβ≥ 10%) (Table 4). A heatmap for the differentially methylated epigenetic clock CpG sites (P < 0.05, Δβ≥ 10%), with hierarchical cluster analysis results of normal, adenoma and CRC samples is shown in Figure 1.
Figure 1 DNA methylation heatmap of normal, adenoma and colorectal carcinoma samples according to the methylation status of age-related CpG sites.
From the 353 epigenetic clock CpG sites (cg IDs)[2] significantly differentially methylated in CRC vs normal, adenoma vs normal and CRC vs adenoma comparisons were selected and colonic tissue samples (GSE48684[24]) were classified according to their methylation levels. The analyzed samples are illustrated on X axis, significantly altered CpG sites (cg IDs) are represented on Y axis. DNA methylation intensities (β values) are visualized, red shows hypermethylation, while hypomethylation was marked with green color. CRC: Colorectal carcinoma (light green); Ad: Adenoma (dark blue); N: Normal (light blue).
Table 2 Significant DNA methylation alterations of age-related CpG sites in CRC samples compared to normal tissue.
Genes belonging to 353 age-related CpG sites were mapped to 768 Affymerix transcript IDs. In the CRC vs N comparison, 215 “epigenetic clock” genes were found to be significantly differentially expressed (P < 0.05), of which 117 were altered with absolute logFC > 0.5 (70 upregulated such as ERG, MGP, MCAM, CEMIP and 47 downregulated like SFRP1, KRT20, CLEC3B, SYNE1) (Supplementary Table 2A). Expression of 196 “epigenetic clock” genes changed significantly (P < 0.05) in adenoma samples compared to healthy normal controls, 102 with absolute logFC higher than 0.5 (47 overexpressed such as CEMIP, PLK1, CCNF and 55 underexpressed like SFRP1, SDC2, SYNE1) (Supplementary Table 2B). Forty-three genes including MCAM, MGP and AKT3 showed increased expression in CRC compared to adenoma samples, while 17 genes including CEMIP, SPATA18 were downregulated (P < 0.05, absolute logFC > 0.5) (Supplementary Table 2C).
For genes with an inverse relation between gene expression and promoter methylation, the genes with both significant mRNA expression changes with absolute logFC > 0.5 and significant DNA methylation alterations with at least 10% difference were taken into consideration. Based on these criteria, eleven genes, including ERG, MGP, PIPOX, CLEC3B, LTK, SFRP1 and SYNE1 were found to be inversely expressed with the promoter methylation status in CRC compared to normal tissue. Compared to the promoter methylation alterations, the expression of 8 genes, including CEMIP, SPATA18, SDC2, SFRP1 and SYNE1 changed oppositely in AD vs N comparison, while in CRC vs AD tissues 3 genes, namely CEMIP, SPATA18 and SLC28A2 showed this expression pattern. The genes showing an inverse relation between gene expression and DNA methylation in CRC vs normal, AD vs normal and CRC vs AD comparisons are represented in Table 5.
Table 5 Genes showing inverse DNA methylation and gene expression data.
In the comparison of healthy young colonic samples and normal adult tissues, 150 genes showed significantly altered expression from which 94 genes with absolute logFC > 0.5 including overexpressed LTBP3, REEP1, MGP, PLK1, SFRP1, SYNE1 and downregulated PRKG2, PDCD6IP and TMEM56 (P < 0.05) (Supplementary Table 2D). The pattern of expression of several genes including SYNE1, CLEC3B, LTBP3 (Figure 2) and SFRP1 (Figure 3A-C) was similar for increased age as that observed for cancer progression. However, there were some genes such as MGP (Figure 2) which showed similar expression pattern in young people and in cancer patients compared to healthy adult tissue.
Figure 2 Genes showing both age- and carcinogenesis-related expression alterations.SYNE1 (spectrin repeat containing nuclear envelope protein 1), LTBP3 (latent transforming growth factor beta binding protein 3) and CLEC3B (C-type lectin domain family 3 member B) genes were downregulated during the colorectal carcinogenesis, similar decreasing expression was found during aging (significantly higher mRNA levels were detected in young colonic samples than in healthy adult biopsy specimens). MGP (matrix Gla protein) was overexpressed in children and in CRC samples compared to adenoma and healthy adults (P < 0.035), hence its opposite expression was found during aging and colorectal carcinogenesis. X axis shows the analyzed sample groups, the normalized mRNA expression can be seen on Y axis. Red dots indicate the normalized mRNA expression values, boxplots represent the medians and standard deviations. Ch: Children; N: Normal; Ad: Adenoma; CRC: Colorectal cancer.
Figure 3 SFRP1 mRNA expression and promoter DNA methylation alterations during aging and in different stages of colorectal carcinogenesis.SFRP1 mRNA expression was significantly downregulated in adenoma and CRC samples compared to normal controls in case of all three Affymetrix probeset IDs representing SFRP1 [202035_s_at: P < 0.05 (A); 202036_s_at: P < 0.0003 (B); 202037_s_at: P < 0.005 (C)]. In colonic biopsy samples of healthy young patients, higher SFRP1 mRNA levels could be measured than in normal adults samples, this overexpression was proven to be significant in two of three transcript IDs [202035_s_at: P < 0.05 (A); 202037_s_at: P < 0.005 (C)]. SFRP1 promoter region was significantly hypermethylated in CRC and adenoma tissue samples compared to normal adult and young colonic tissue (P < 0.0001) (D). In pairwise comparison, DNA methylation of SFRP1 promoter was slightly, but significantly increased in healthy adults compared to normal young samples (P < 0.02) (E). The analyzed sample groups are illustrated on X axis, the normalized mRNA expression (A, B, C) and percentage of SFRP1 promoter methylation (D, E) are represented on Y axis. Red dots indicate the normalized mRNA expression values (A, B, C) and the normalized DNA methylation percentage values (D, E), respectively. Boxplots represent the medians and standard deviations. Ch: Children; N: Normal; Ad: Adenoma; CRC: Colorectal cancer.
Whole promoter methylation status of genes showing inverse relationship between gene expression and DNA methylation
The DNA methylation status of whole promoter regions of genes showing an inverse relation between gene expression and DNA methylation (Table 5) was determined using methyl capture sequencing data of 6 NAT, 15 adenoma and 9 CRC tissue samples[12]. In the CRC vs N/NAT comparison, similar DNA methylation alterations (such as hypomethylated AKT3, MGP promoters and hypermethylated PPP1R16B, SFRP1, SYNE1 promoters) were observed in the promoter regions of 7 of the 11 inversely expressed genes (Table 6). Between adenoma and normal samples, promoter regions of 7 of the 8 inversely expressed genes showed DNA methylation differences (e.g., hypomethylated CEMIP/KIAA1199, SPATA18 promoters and hypermethylated DKK3, SDC2, SFRP1, SYNE1 promoters) (Table 6) as detected in case of age-related CpG sites. In CRC samples compared to adenoma tissue, significant hypermethylation of CEMIP/KIAA199 promoter could be demonstrated in “epigenetic clock” CpG sites and whole promoter methylation analyses (Table 6).
Table 6 Whole promoter methylation status of genes with inverse age-related CpG site DNA methylation and gene expression data.
SFRP1 promoter methylation analysis in healthy children, healthy adult, adenoma and CRC tissues
Based on the gene expression analysis results, SFRP1 was found to be overexpressed in normal adult samples compared to adenoma and CRC biopsy specimens, and in healthy young patients even higher SFRP1 mRNA levels could be detected than in normal adult samples (Figure 3A-C). As SFRP1 is proven to be a methylation-regulated gene with literature data regarding its age-related DNA methylation alterations, hence it was chosen for detailed methylation analysis of normal, premalignant and cancerous colonic specimens including tissue samples from healthy children.
SFRP1 promoter sequences were highly methylated in CRC samples (average methylation% = 55.0% ± 8.4%) and in adenoma tissue (49.9% ± 18.1 %), while low methylation levels could be measured in colonic tissues of both healthy adults (5.2% ± 2.7%) and children (2.2% ± 0.7%). Significant considerable hypermethylation was found in SFRP1 promoter both between CRC and adult normal and between CRC and healthy children colonic tissue samples (P < 0.0001) (Figure 3D). In the healthy adult vs healthy children comparison, significant, but moderate DNA methylation alterations were detected: in adults higher DNA methylation levels were found in the analyzed region of SFRP1 promoter (P = 0.017) (Figure 3E).
SFRP1 protein expression in colonic tissue samples of healthy normal children and adults
In healthy children samples the epithelial layer showed strong (representative scoring values: +3 and +2), diffuse SFRP1 expression both in cytoplasmic and nuclear region (Q-score: 226.67 ± 17.51), whereas the stromal cells showed heterogeneous protein expression (scoring values: form +3 to 0; Q-score: 176.66 ± 18.61; ΣQ-score: 403.33 ± 22.51; Figure 4A). Among stromal cells subepithelial fibroblast and several immune cells showed strong cytoplasmic and/or nuclear SFRP1 expression. Not significantly, but remarkably lower (representative scoring value: +2) SFRP1 protein expression was detectable both in epithelial (Q-score: 202.14 ± 24.12) and stromal component (140.71 ± 41.47; Σ Q-score: 351.42 ± 68.66; Figure 4B) of adult persons (P values: 0.063, 0.073 and 0.105 respectively).
Figure 4 SFRP1 protein expression in colonic tissue samples of healthy normal children and adults.
Strong/moderate cytoplasmic and nuclear SFRP1 expression both in epithelial and stromal compartments of healthy children (A) and healthy adult (B) samples. Digital microscopic images; magnification × 40; scale: 50 μm.
DISCUSSION
Cancer is considered a primary age-related disease[3], and therefore age-related molecular changes including epigenetic alterations such as epigenetic drift and epigenetic clock[1,2] necessarily show relationship with carcinogenesis-associated differences. Besides global hypomethylation, local, genomic site specific hypermethylation principally in the promoter regions of tumor suppressor genes can occur during both processes[6-8].
In this study, the potential correspondence between age-related and CRC-associated DNA methylation changes was studied using the 353 epigenetic clock CpG sites published by Horvath[2] as a model for age-related DNA methylation changes. With the analysis of methyl capture sequencing and Illumina BeadChip450K methylation array data, the methylation status of age-related CpG sites and genes was determined during CRC development and progression, and the relevant mRNA expression changes were also evaluated. Among the differentially methylated/expressed age-related genes, SFRP1 promoter methylation was further analyzed in healthy, premalignant and cancerous colonic tissue samples, including biopsy specimens from young children.
Similarly to previous findings[29], DNA methylation alterations in a considerable proportion of age-related CpG sites/gene promoters (approximately 40%) were observed in samples representing different stages of CRC formation and progression. Approximately two third of these CpG sites had similar DNA methylation alterations in CRC compared to normal tissue samples as during aging. When the effect of DNA methylation of epigenetic clock genes was studied, whole promoter methylation was also observed in addition to the analysis of DNA methylation status of representing age-related CpG sites.
In accordance with our results, hypermethylation of several genes belonging to aging-associated CpG sites such as SFRP1[7,12,15,18-23], TFAP2E[30], TBX5[31], GNAS[32], DKK3[18,33], DGKI[34], SYNE1[35-37], SDC2[38,39], ADHFE1[40-42] was observed in tissue and/or blood samples of CRC patients. SFRP1 tumor suppressor protein with a putative Wnt-binding site impedes the frizzled ligand - Wnt receptor interaction. Its reduced expression caused by promoter hypermethylation can lead to constitutive activation of Wnt pathway which is best characterized signaling pathway in CRC pathogenesis[7,12,15,18-23,43]. Worthley et al[15] showed strong positive correlation between SFRP1 methylation and age (Spearmen’s rank P = 0.72, P < 0.0001) on a set of 166 CRC tissue samples from adults [median age: 61.2 years (22.8-89.2 years)]. In this study, the increase of SFRP1 methylation during the aging was also observed, moreover to our best knowledge, we provide the first evidence of significantly lower SFRP1 methylation levels in children (under 18 years) compared to healthy adult colonic tissues. Preliminary results of methylation-sensitive restriction enzyme methylation array analysis of our reseach groups suggested SFRP1 hypomethylation in young colonic tissue samples, though the high standard deviations of methylation percentages and the low samples size limited our conclusions[18]. In accordance with the DNA methylation data, elevated mRNA and protein levels could be detected in colonic tissues of normal children compared to adults.
Dickkopf Wnt signaling pathway inhibitors including DKK3 are also frequent targets of epigenetic silencing in gastrointestinal tumors promoting carcinogenesis by loss of/reduction their expression[33]. Hypermethylated syndecan 2 (SDC2) is a potential biomarker for early CRC detection both in serum and tissue, although the gene silencing effect of elevated promoter methylation is not unambiguous according to the literature data[44], moreover some observations support its tumorigenic activity in CRC[45]. Hypermethylation of spectrin repeat containing nuclear envelope protein 1 (SYNE1) suggests its tumor suppressor function in CRC[37], which was detected not only in CRC tissue samples[35,36], but it appears to be a promising marker for blood-based CRC detection[37]. Alcohol dehydrogenase, iron containing 1 (ADHFE1) promoter hypermethylation was found to be associated with CRC differentiation[41], furthermore it is involved in cell proliferation induction by alcohol in colon carcinoma cells[46].
In the case of some genes like AKT3, CEMIP and DDR1, promoter hypomethylation was observed in different types of cancers such as breast cancer, lung cancer and CRC[47-49]. PI3K/Akt pathway is thought to be the most commonly activated intracellular signaling pathway in human malignancies[50]. AKT kinases including AKT3 are remarkable contributors to malignant diseases as they are involved in the regulation of cell proliferation, growth and survival[50,51]. Hypomethylation and overexpression of the cell migration-inducing protein (CEMIP/KIAA1199) was previously described in CRC[49,52]. In our study, hypomethylation of this gene could be detected mainly in adenoma samples, however a slight but significant decrease in methylation level was observed in CRC samples compared to normal controls. In accordance with the methylation data, strong upregulation of CEMIP mRNAs was shown both in adenoma and CRC samples with higher expression values in adenoma tissue[53,54]. Due to its robust overexpression at mRNA and also at protein levels, CEMIP is considered as a candidate prognostic marker for CRC and a potential therapeutic target[55]. CEMIP facilitates colon cancer cell proliferation via enhancing Wnt signaling[49] and promotes tumor growth[55] and cancer dissemination under hypoxic conditions[56].
In conclusion, our results regarding DNA methylation alterations of age-related, epigenetic clock genes during colorectal carcinogenesis supports the concept that aging is one of the main factors predisposing cancer including CRC. Several age-related DNA methylation alterations could be observed during development and progression of CRC affecting the mRNA expression of certain CRC- and adenoma-related key control genes. The main CRC-associated signal transduction pathways, such as WNT signaling and PI3K/Akt pathways are also influenced during aging.
ACKNOWLEDGMENTS
We thank hereby the help of Gabriella Kónyáné Farkas with the immunohistochemical analysis. Furthermore we thank Theo deVos PhD for his careful language assistance.
COMMENTS
Background
Cancer is considered a primary age-related disease, and therefore age-related molecular changes including epigenetic alterations such as epigenetic drift and epigenetic clock necessarily show relationship with carcinogenesis-associated differences. Besides global hypomethylation, local, genomic site specific hypermethylation principally in the promoter regions of tumor suppressor genes can occur during both processes.
Research frontiers
Several age-related DNA methylation alterations could be observed during colorectal cancer (CRC) formation and progression affecting the mRNA expression of certain CRC- and adenoma-related key control genes such as hypermethylated secreted frizzled related protein 1 (SFRP1), spectrin repeat containing nuclear envelope protein 1 and hypomethylated cell migration-inducing protein.
Innovations and breakthroughs
For the first time significantly lower SFRP1 methylation levels were demonstrated in colonic tissue from children (under 18 years) compared to healthy adults. The main CRC-associated signal transduction pathways, such as WNT signaling and PI3K/Akt pathways are also influenced during aging.
Peer-review
In this paper, the authors analyzed the methylation and expression levels of 353 age-related “epigenetic clock” genes in colonic tissue samples. They identified many differentially methylated and/or differentially expressed genes. Among these genes, the DNA methylation and mRNA levels of SFRP1 was further analyzed. This is an interesting work using a large number of data and the results may be useful in related field.
Footnotes
Manuscript source: Invited manuscript
Specialty type: Gastroenterology and hepatology
Country of origin: Hungary
Peer-review report classification
Grade A (Excellent): 0
Grade B (Very good): B
Grade C (Good): C, C, C, C
Grade D (Fair): 0
Grade E (Poor): 0
P- Reviewer: Essani K, Lakatos PL, M’Koma A, Xiao Y, Lakatos PL S- Editor: Yu J L- Editor: A E- Editor: Wang CH
Galamb O, Kalmár A, Péterfia B, Csabai I, Bodor A, Ribli D, Krenács T, Patai ÁV, Wichmann B, Barták BK, Tóth K, Valcz G, Spisák S, Tulassay Z, Molnár B. Aberrant DNA methylation of WNT pathway genes in the development and progression of CIMP-negative colorectal cancer.Epigenetics. 2016;11:588-602.
[RCA] [PubMed] [DOI] [Full Text][Cited by in Crossref: 57][Cited by in RCA: 58][Article Influence: 6.4][Reference Citation Analysis (0)]
Patai ÁV, Valcz G, Hollósi P, Kalmár A, Péterfia B, Patai Á, Wichmann B, Spisák S, Barták BK, Leiszter K, Tóth K, Sipos F, Kovalszky I, Péter Z, Miheller P, Tulassay Z, Molnár B. Comprehensive DNA Methylation Analysis Reveals a Common Ten-Gene Methylation Signature in Colorectal Adenomas and Carcinomas.PLoS One. 2015;10:e0133836.
[RCA] [PubMed] [DOI] [Full Text] [Full Text (PDF)][Cited by in Crossref: 39][Cited by in RCA: 40][Article Influence: 4.0][Reference Citation Analysis (0)]
Toyota M, Ahuja N, Ohe-Toyota M, Herman JG, Baylin SB, Issa JP. CpG island methylator phenotype in colorectal cancer.Proc Natl Acad Sci USA. 1999;96:8681-8686.
[PubMed] [DOI]
Beggs AD, Dilworth MP, Domingo E, Midgley R, Kerr D, Tomlinson IP, Middleton GW. Methylation changes in the TFAP2E promoter region are associated with BRAF mutation and poorer overall & amp; disease free survival in colorectal cancer.Oncoscience. 2015;2:508-516.
[PubMed] [DOI]
