Global DNA hypomethylation of colorectal tumours detected in tissue and liquid biopsies may be related to decreased methyl-donor content
BMC Cancer volume 22, Article number: 605 (2022)
Hypomethylation of long interspersed nuclear element 1 (LINE-1) is characteristic of various cancer types, including colorectal cancer (CRC). Malfunction of several factors or alteration of methyl-donor molecules’ (folic acid and S-adenosylmethionine) availability can contribute to DNA methylation changes. Detection of epigenetic alterations in liquid biopsies can assist in the early recognition of CRC. Following the investigations of a Hungarian colon tissue sample set, our goal was to examine the LINE-1 methylation of blood samples along the colorectal adenoma-carcinoma sequence and in inflammatory bowel disease. Moreover, we aimed to explore the possible underlying mechanisms of global DNA hypomethylation formation on a multi-level aspect.
LINE-1 methylation of colon tissue (n = 183) and plasma (n = 48) samples of healthy controls and patients with colorectal tumours were examined with bisulfite pyrosequencing. To investigate mRNA expression, microarray analysis results were reanalysed in silico (n = 60). Immunohistochemistry staining was used to validate DNA methyltransferases (DNMTs) and folate receptor beta (FOLR2) expression along with the determination of methyl-donor molecules’ in situ level (n = 40).
Significantly decreased LINE-1 methylation level was observed in line with cancer progression both in tissue (adenoma: 72.7 ± 4.8%, and CRC: 69.7 ± 7.6% vs. normal: 77.5 ± 1.7%, p ≤ 0.01) and liquid biopsies (adenoma: 80.0 ± 1.7%, and CRC: 79.8 ± 1.3% vs. normal: 82.0 ± 2.0%, p ≤ 0.01). However, no significant changes were recognized in inflammatory bowel disease cases. According to in silico analysis of microarray data, altered mRNA levels of several DNA methylation-related enzymes were detected in tumours vs. healthy biopsies, namely one-carbon metabolism-related genes—which met our analysing criteria—showed upregulation, while FOLR2 was downregulated. Using immunohistochemistry, DNMTs, and FOLR2 expression were confirmed. Moreover, significantly diminished folic acid and S-adenosylmethionine levels were observed in parallel with decreasing 5-methylcytosine staining in tumours compared to normal adjacent to tumour tissues (p ≤ 0.05).
Our results suggest that LINE-1 hypomethylation may have a distinguishing value in precancerous stages compared to healthy samples in liquid biopsies. Furthermore, the reduction of global DNA methylation level could be linked to reduced methyl-donor availability with the contribution of decreased FOLR2 expression.
Colorectal cancer (CRC) is the third most commonly diagnosed cancer type and possesses the second highest mortality worldwide, as reported by the 2020 GLOBOCAN survey .
The vast majority of sporadic CRCs develop from healthy epithelium through benign colorectal adenoma (AD) stages into CRCs . CRC progression can also be associated with the presence of long-lasting inflammatory bowel diseases (IBD), including ulcerative colitis and Crohn’s disease [3, 4]. The prevalence of IBD is increasing continuously worldwide; thus, it is necessary to investigate the molecular background of this disease more thoroughly .
In most CRC cases, the main molecular driving force in tumour development is chromosomal instability, where extensive DNA hypomethylation can be detected as an early molecular event, contributing to the formation of chromosomal instability [6,7,8]. As a consequence of global DNA hypomethylation, mobile genetic elements such as long interspersed nuclear element 1 (LINE-1) can be released from the inhibition induced by DNA methylation [9,10,11]. As copies of LINE-1 retrotransposon compose about 17% of the human genome and their methylation level correlates well with genome-wide DNA methylation [12, 13], LINE-1 bisulfite pyrosequencing is utilised for global DNA methylation estimation. Minimally invasive diagnostic tests, such as the detection of cancer-specific molecular markers in cell-free DNA (cfDNA) fraction by blood sampling, can contribute to the early recognition of cancer formation and thereby can increase the survival rate [14,15,16,17,18,19,20,21,22]. Due to this fact and the prognostic value of LINE-1 methylation status in various cancer types , it is important to explore the possible clinical potential of LINE-1 hypomethylation detection.
DNA methyltransferases, demethylases, chromatin remodelling enzymes, transcription factors, non-coding RNAs, and methyl-donor molecules are necessary for maintaining a cell-specific DNA methylation pattern . Malfunctions related to any of these factors can contribute to abnormal global DNA hypomethylation . In methylation processes, the methyl group originates from the methyl-donor S-adenosylmethionine (SAM) molecule, which is produced by the methionine-cycle as part of one-carbon metabolism . One-carbon metabolism consists of several closely connected pathways centered around folic acid (FA), which is responsible for carrying single carbon units, like methyl groups .
Several studies have focused on the influence of methyl-donor molecule intake on CRC development risk revealing inverse correlations [26,27,28,29,30]. However, results regarding the effect of methyl-donors on global DNA- and LINE-1 methylation are controversial. No LINE-1 methylation alteration was detected, ensuing folate supplementation in the case of normal colon tissue [31, 32], but a slight positive association was observed between serum folate level and global DNA methylation according to Pufulete et al. . Furthermore, a trend of association was noticed between LINE-1 methylation level and folate intake in studies about patients possessing colorectal tumours with DNA hypomethylation [34, 35]. It is also important to note that folate supply could have a different effect on pre-established neoplasms than on normal tissue, as it can contribute to the progression of the former and the maintenance of the latter . As global DNA hypomethylation is closely connected to chromosomal instability along with CRC development, and folic acid intake may have a protective effect before neoplastic transformation, the investigation of the possible underlying mechanisms—including aberrant methyl-donor availability—is crucial.
The present study aims to comprehensively evaluate LINE-1 methylation on a Hungarian sample cohort along the colorectal adenoma-carcinoma sequence using LINE-1 bisulfite-sequencing as a basis for our further examinations. Moreover, we involved liquid biopsy samples to analyse the potential diagnostic value of LINE-1 hypomethylation, including the investigation of patients with IBD due to its high incidence rate and risk of CRC development. Finally, we aimed to explore the possible molecular background behind the emergence of global DNA hypomethylation during colorectal tumour progression on a multi-level aspect. Our purpose was to combine microarray analysis and immunohistochemistry staining to discover possible connections between DNA methylation-related enzymes’ mRNA expression profile, methyl-donor (FA and SAM) content, and global DNA hypomethylation.
Material and methods
A total of 183 fresh frozen (FFT), 40 formalin-fixed paraffin-embedded tissue (FFPET) biopsies, and 48 blood samples taken in K3EDTA Vacuette tubes (Greiner Bio-One Gmbh) were collected at the Department of Internal Medicine and Hematology and the 1st Department of Pathology and Experimental Cancer Research, Semmelweis University, Budapest, Hungary.
The FFT sample cohort contained 45 healthy (N), 23 normal adjacent to colorectal adenoma tissue (AD-NAT), 25 normal adjacent to colorectal carcinoma tissue (CRC-NAT), 37 AD, 38 CRC, and 15 IBD biopsy samples. AD samples included different histological subtypes, namely tubular (TA) and tubulovillous (TVA) biopsies, while the CRC group contained specimens with multiple Astler-Coller modified Dukes’ A-B (early) and Astler-Coller modified Dukes’ C-D (late) stages. NAT biopsies were taken 10–12 cm distance from the lesions. FFPET biopsies included 20 AD and 20 CRC sections containing transitional zones. Moreover, plasma samples of 10 healthy, 14 AD, 13 CRC, and 11 IBD patients were also involved in this study. Exclusion criteria were acute health problems, cancerous diseases apart from CRC, radiotherapy, and chemotherapy. The clinicopathological and demographic data of the investigated patients are summarized in Additional file 1 (Supplemental material Table 1, 2, 3).
The study has been approved by the local ethics committee (Regional and Institutional Committee of Science and Research Ethics; TUKEB Nr: 14,383–2/2017/EKU), and all applied methods were performed in accordance with the relevant guidelines and regulations. Written informed consent was obtained from all the patients prior to sample collection.
DNA isolation from fresh frozen tissue and plasma samples
Genomic DNA was isolated from FFT biopsies with High Pure PCR Template Preparation Kit (Roche) according to the manufacturer’s instruction with overnight proteinase K digestion.
Plasma fractions were separated from whole blood samples by two centrifugation steps (1350 rcf, 12 min). CfDNA isolation was carried out from 3.5 ml plasma samples with Quick-cfDNA™ Serum & Plasma Kit (Zymo Research) according to the manufacturer's recommendation. DNA samples were stored at -20 °C until further examinations. Concentration and purity (OD260/280, OD230/280) of genomic DNA were measured by NanoDrop ND-1000 Spectrophotometer (Thermo Fisher Scientific), while cfDNA quantification was performed with Qubit 1.0 fluorometer using Qubit dsDNA HS Assay Kit (Thermo Fisher Scientific).
LINE-1 bisulfite pyrosequencing
After DNA isolation, 500 ng DNA of each tissue sample and approximately 20 ng cfDNA of each plasma sample were bisulfite-converted using EZ DNA Methylation-Direct Kit (Zymo Research) according to the manufacturer's protocol. A 146 basepair-long region of LINE-1 sequence was amplified by bisulfite-specific PCR using Pyromark Q24 CpG LINE-1 Kit (Qiagen) on Mastercycler EP Gradient S PCR machine (Eppendorf, Hamburg, Germany) with the following thermocycling program: PCR activation for 15 min on 90 °C; 45 cycles of 30 s denaturation on 94 °C, 30 s annealing on 50 °C and 30 s extension on 72 °C; 10 min final extension on 72 °C. The specificity of PCR was verified by gel electrophoresis using 2% agarose gel. Amplicons were prepared and washed on PyroMark Q24 Vacuum Workstation (Qiagen) and pyrosequenced on Pyromark Q24 (Qiagen) instrument using Pyromark Q24 CpG LINE-1 Kit (Qiagen) and PyroMark Gold Q24 Reagents (Qiagen). Methylation of three cytosine-guanine dinucleotides (positions 318, 321, and 328 of LINE-1, GenBank accession number: X58075) was calculated as the percentage of cytosine nucleotides relative to the sum of cytosine and thymine nucleotides by Pyromark Q24 software v2.0.6. The average methylation level of the three cytosine-guanine dinucleotides was assessed as the LINE-1 methylation status of each sample. All pyrosequencing and analysis steps were performed according to the manufacturer's description.
Statistical significance (p ≤ 0.05) and receiver operating characteristic (ROC) curve analysis were calculated using Prism8 software (GraphPad). In multiple comparison tests, data with non-normal distribution were analysed by Kruskal–Wallis and Dunn's multiple comparison tests, while in the case of normal distribution, ANOVA followed by Brown Forsythe and Welch ANOVA tests were applied. In pairwise comparisons, paired t-test or Wilcoxon matched-pairs signed rank tests were performed depending on normal or non-normal data distribution, respectively.
In silico gene expression analysis
Human Transcriptome Array (HTA) 2.0 RNA microarray (Affymetrix) was performed by Kalmár et al. according to the manufacturer’s recommendations to detect the mRNA expression profile of 20 colorectal normal, 20 AD, and 20 CRC samples . The dataset can be found on Gene Expression Omnibus (GEO) data repository (GEO ID GSE100179) . Processed data were reanalysed with Transcriptome Analysis Console 4.0 software (TAC 4.0, Affymetrix). Criteria of our in silico analysis were the following: adjusted p-value ≤ 0.05; fold change (FC) ≥|1.4|.
Immunohistochemistry staining was applied on 5 µm sections obtained from 20 AD and 20 CRC FFPET samples. After deparaffinization, heat-induced epitope retrieval was performed. Tissue sections were treated with 3% hydrogen peroxide solution for 30 min. Nonspecific binding sites were blocked with 1% BSA solution for 60 min at room temperature. According to the manufacturers’ protocol, the following primary antibodies were used: anti-DNMT1 (AB19905, Abcam), anti-DNMT3a (ab16704, Abcam), anti-DNMT3b (ab122932, Abnova), and anti-FOLR2 (ab228643, Abcam, 1:500), which is specific to folate receptor beta (FR-beta) encoded by FOLR2. In the case of FA, SAM, and 5-methylcytosine (5mC) staining, sections were firstly treated with 0.1% Tween-TBS for 10 min after 30 min of 3% hydrogen peroxide incubation. Antigen retrieval was performed with 2 N HCl solution for 15 min at room temperature. The pH was normalized using 100 mM TBS (pH 8.3). Nonspecific binding sites were blocked with 1% BSA for 30 min. Primary antibodies were applied according to the manufacturers’ recommendation (anti-folic acid (PAB0092, Abnova, 1:500), anti-S-adenosylmethionine (Clone 118–18, 6942–25, BioVision, 1:200), anti-5-methylcytosine (ab214727, Abcam, 1:500).
Secondary antibodies were utilised for 1 h at room temperature: anti-rabbit antibody (BA1054, Boster), anti-rat antibody (18–4818-82, Invitrogen), and anti-mouse antibody (62–6520, Invitrogen). The visualisation was performed using DAB Chromogen (K3468, DAKO) with HISTOLS-DAB Substrate (30,014. K, Histopathology Ltd.), and nuclei were stained with hematoxylin.
Slides were digitised with Pannoramic Confocal instrument and were analysed using Pannoramic Viewer digital microscope (v: 1.15.3, 3DHISTECH). The intensity was classified from 0 to + 3, where 0 means negative, while + 1 is a weak, + 2 is an intermediate, and + 3 is a strong immunohistochemical signal. To evaluate the staining, modified Quick-score (Q-score) method was applied (Formula: Q-score = IxP, P: percentage of the positive cells, I: staining intensity of these cells, maximum value: 300 (3 × 100) [38, 39]). Statistical significance (p ≤ 0.05) was assessed by paired t-test, or Wilcoxon matched-pairs signed rank test depending on normal or non-normal data distribution, respectively (Prism8 software, GraphPad).
LINE-1 methylation in tissue samples and plasma fraction
Significant decrease of LINE-1 methylation was found in AD (72.7 ± 4.8%) and in CRC (69.7 ± 7.6%) samples compared to N (77.5 ± 1.7%) colon and IBD (77.1 ± 1.9%) tissue specimens according to LINE-1 bisulfite sequencing (p ≤ 0.01). However, there was no significant methylation alteration in IBD biopsies compared to N (Fig. 1/A).
AD-NAT and CRC-NAT specimens (76.0 ± 2.1%, 76.4 ± 2.0%, respectively) showed significantly reduced methylation levels compared to the samples of healthy individuals (77.5 ± 1.7%) (p ≤ 0.05) (Additional file 2/A). Furthermore, the LINE-1 methylation levels of AD (72.8 ± 4.4%) and CRC (71.2 ± 6.8%) tissue samples decreased significantly compared to their matched NAT pairs (76.0 ± 2.1%, 76.4 ± 2.0%, respectively) (p ≤ 0.05) (Additional file 2/B).
Significantly lower LINE-1 methylation was found in TA (73.8 ± 4.1%) and TVA (69.8 ± 4.5%) in comparison with N samples (77.5 ± 1.7%) (p ≤ 0.05), and a decreasing trend was observed in TVA (69.8 ± 4.5%) vs. TA (73.8 ± 4.1%) (Additional file 2/C, left). LINE-1 methylation levels of both early (71.5 ± 4.9%) and late (68.5 ± 8.7%) carcinomas were significantly reduced compared to N specimens (p ≤ 0.05) (Additional file 2/C, right). A decline was also noticed in late cancer stages (68.5 ± 8.7%) in comparison with early carcinoma (71.5 ± 4.9%) biopsies.
According to our observations, AD (80.0 ± 1.7%) and CRC (79.8 ± 1.3%) plasma specimens showed significantly decreased methylation status compared to N (82.0 ± 2.0%, p ≤ 0.01), while no significant alteration was observed in IBD (81.1 ± 0.8%) in comparison with N samples (Fig. 1/B, left). Receiver operating characteristic (ROC) curve analysis was applied in AD liquid biopsy samples to investigate the potential of LINE-1 methylation as a biomarker that already indicates the precancerous AD stages. At the cut-off value 80.0%, LINE-1 methylation percentage distinguished AD from N with 66.7% sensitivity and 90.0% specificity (AUC: 0.8, p < 0.05, 95% CI 0.6–1.0) (Fig. 1/B, right).
In silico gene expression analysis of genes related to DNA methylation
GeneChip™ HTA 2.0 results of our research group were reanalysed in silico to investigate the mRNA expression profile of folate receptors, one-carbon cycle-, and DNA methylation-related enzymes .
Only candidates showing significant expression changes (p ≤ 0.05) in both AD and CRC tissue samples compared to N were illustrated in Fig. 2. Five genes involved in one-carbon metabolism (GART, MTHFD1, MTHFD2, ATIC, and SHMT2) were found to be overexpressed in tumorous samples in comparison with healthy controls, while FOLR2 owned significantly reduced gene expression alterations with -1.45 and -1.53 FC values in AD vs. N and CRC vs. N comparisons (p ≤ 0.05), respectively (Fig. 2.).
DNMT1 showed significant mRNA expression level elevation with 1.51 FC value in CRC vs. N comparison. APOBEC1 also possessed increased transcription with 1.88 FC value in AD samples compared to N specimens (p ≤ 0.05). APOBEC3B had significantly decreased transcript level (FC = -1.49) in AD vs. healthy biopsies (p ≤ 0.05). mRNA level alterations of DNMT1, APOBEC1, APOBEC3B are represented on Additional file 3.
Validation of DNMTs and FOLR2 GeneChip™ Human Transcriptome Array 2.0 results
Gene expression results of DNMTs and FOLR2 were validated with immunohistochemistry staining
DNMT1 labelling showed significantly elevated intensity in CRC epithelial cells (p ≤ 0.05), while no expression changes were detected in stromal cells of CRC compared to NAT. No significant alterations could be observed in DNMT3A and DNMT3B immunostaining between the analysed areas (Additional file 4).
Significantly lower (p ≤ 0.01) FOLR2-encoded FR-beta expression was detected in epithelial and stromal cells of AD and epithelial cells of CRC specimens compared to their NAT areas (Additional file 5). In CRC stromal cells, a moderate reduction of staining intensity was noticed in comparison with NAT regions.
Analyses of methyl-donor content in colon tissue
Immunohistochemistry staining was applied to investigate in situ methyl-donor (FA and SAM) content and examine its relation to the 5mC level. Significantly reduced intensity of FA (Fig. 3/A, B), SAM (Fig. 3/C, D), and 5mC (Fig. 3/E, F) was observed in epithelial cells of both AD and CRC vs. NAT areas (p ≤ 0.01). Slightly diminished labelling was noticed in stromal cells except for 5mC staining in AD (Fig. 3/E) (p ≤ 0.01), and CRC (Fig. 3/F) (p ≤ 0.01), along with FA labelling in CRC sections (Fig. 3/B) (p ≤ 0.05), where also significantly lower intensity was detected.
Significant LINE-1 hypomethylation has been known in several cancers compared to healthy tissue specimens . In the present study, we focused specifically on tissue samples collected from Hungarian individuals, revealing significant LINE-1 hypomethylation in both AD and CRC biopsies compared to N, paired AD-NAT, and CRC-NAT tissues. Moreover, our measurements of AD-NAT, CRC-NAT, and N colonic biopsies showed significantly lower LINE-1 methylation in NAT samples vs. healthy controls, which also agreed with the recent scientific literature [8, 41, 42].
Significant LINE-1 hypomethylation was noticed in TA and TVA samples compared to N, and tendentiously lower methylation levels were found in TVA specimens compared to TA biopsies. In contrast with the latter result, no LINE-1 methylation alteration was detected in TVA vs. TA comparison, according to the observations of Jiang et al. . Our result can be explained by the fact that ADs with more villous features have a higher risk of cancerous transformation , and global DNA methylation is continuously decreasing with CRC progression .
Significant reduction of LINE-1 methylation was detected in the samples of patients with early and late CRC stages compared to N specimens, while moderately diminished methylation level was observed within the CRC group (Dukes C-D vs. Dukes A-B). This is consistent with the fact that global DNA methylation level reduces during CRC development [45, 46].
The molecular examination of IBD is crucial due to its increasing incidence and its relation to the risk of CRC development [5, 47, 48]. Glória et al. found global DNA hypomethylation in rectal mucosa of ulcerative colitis patients compared to those who had negative colonoscopy . As a novelty of the present study, no significant LINE-1 methylation changes were observed in the IBD vs. N comparison. Nevertheless, AD and CRC tissue samples showed significantly lower LINE-1 methylation levels than IBD. Similar to our observations, significantly declined LINE-1 methylation level was identified in gastric cancer compared to paired NAT samples, but no significant alteration was noticed between healthy gastric, NAT, and chronic gastritis tissue specimens .
The research of Nagai et al. was the first to describe LINE-1 hypomethylation as a possible diagnostic biomarker for early CRC detection using plasma samples . As a novelty of our study, significant LINE-1 hypomethylation was also found in AD besides CRC samples compared to healthy controls in cfDNA fraction. Similar to our findings in inflamed tissue samples, no significant methylation alteration was detected in IBD vs. N plasma samples. Based on these results and the literature data mentioned above, LINE-1 hypomethylation can be considered a characteristic feature of tumour formation rather than inflammation, which can be advantageous in the development of a cancer-specific screening procedure based on LINE-1 hypomethylation. Since reduced LINE-1 methylation is characteristic in various cancer types, it may not serve as a specific diagnostic marker for CRC identification, but it can be useful in a personalized, minimally invasive monitoring method to observe AD formation, CRC progression, or remission . The low number of investigated plasma samples is the limitation of this study; however, our results may expand the opportunities for monitoring colorectal tumour development with a more comfortable method for patients.
The phenomenon of LINE-1 hypomethylation is widely examined, but its development has not been fully understood yet. Thus, the expression alterations of DNA methylation-related enzymes and the in situ level of methyl-donor molecules were examined. Genes that met our investigating criteria during microarray analysis are discussed below. Upregulation of DNA demethylase enzymes can result in DNA hypomethylation . Hence, the transcript level of TETs and APOBECs encoding genes were analysed [53,54,55,56], from which only APOBEC1 and APOBEC3B had significant alterations. mRNA level of APOBEC1 increased, and APOBEC3B diminished significantly in benign compared to normal specimens. As a higher transcript level of a catalytically inactive, shorter APOBEC1 splice-variant has been observed in colorectal cancer before [57, 58], the elevated mRNA expression level also detected in our study may not have an effect on the global DNA methylation level. Regarding APOBEC3B, declined mRNA expression was noticed in CRC vs. N samples by Burns et al. as well . In normal cells, APOBEC3 enzymes can inhibit the retrotransposition of LINE-1 sequences in a manner not coupled to DNA methylation [60, 61]. Consequently, the reduced level of APOBEC3B in colorectal cancer may contribute to the accumulation of LINE-1 copies not through lessened DNA demethylase activity but through decreased restriction of retrotransposons’ mRNA.
SAM producing methionine-cycle is part of the one-carbon metabolism ; therefore, the associated genes were also analysed in this study. Candidates of one-carbon cycle with mRNA expression alterations that met our criteria (SHMT2, MTHFD1, MTHFD2, ATIC, and GART) showed significant upregulation in AD and CRC vs. N cases. Overexpression of one-carbon metabolism-related enzymes is observed in many cancer types, including colon cancer, which can be explained by the elevated nucleotide demand of the highly proliferating cells .
FR-beta and DNMT enzymes are key factors of the DNA methylation process, as FR-beta provides a unidirectional transport of folate to the cells, while DNMTs are the main catalysts of DNA methylation [63, 64]. Hence, in our study, these enzymes were selected for validation with immunohistochemistry staining. Elevated DNMT1 and no significantly altered DNMT3a, -3b expressions were verified by HTA 2.0 microarray analysis and immunohistochemistry. DNMT1 is responsible for the maintenance of global DNA methylation, while DNMT3a and 3b are liable for local, promoter-specific methylation patterns [65, 66]. Increased DNMT1 expression can be related to a regulatory feedback mechanism, as DNA hypomethylation and elevated DNMT1 mRNA levels were detected in response to inhibited DNA methylation [67,68,69].
Regarding FOLR2, we found reduced mRNA expression in polypoid and cancerous lesions compared to healthy tissues. Immunohistochemistry staining also confirmed decreased FR-beta expression in AD and CRC vs. NAT regions. Similar to our findings, weak staining intensity of tumorous cells was noticed by de Boer et al. . The downregulation of FOLR2 may result in decreased FA and SAM levels, which can impact DNA methylation pattern .
Immunohistochemistry staining was applied to examine the in situ tissue presence of methyl-donor molecules. Evaluation of in situ methyl-donor content in colorectal tissue with immunolabelling technique in parallel with 5mC staining was one of our novelties. FA and SAM labelling showed decreased intensity in adenomatous and cancerous areas compared to NAT regions in line with 5mC staining, which also reduced along the normal-adenoma-carcinoma sequence. Similarly, the study published by Alonso-Aperte et al. revealed decreased folate status in colonic neoplastic mucosa compared to normal and also reduced S-adenosylmethinone/S-adenosylhomocysteine ratio, which means diminished methylation potential . Our observations about decreased SAM levels in tumorous areas can be explained by the fact that folate provides single carbon moieties for one-carbon metabolism, so the lower amount of this molecule can lead to decreased SAM level and lessened global DNA methylation [72, 73]. Even though recent findings regarding the effect of folate deficiency on DNA methylation and cancer progression are controversial , our results support the fact that low methyl-donor content may contribute to the formation of global DNA hypomethylation and CRC development.
In summary, we observed significant LINE-1 hypomethylation along the colorectal adenoma-to-carcinoma continuum in both tissue and peripheral blood specimens in a Hungarian colon sample cohort. Plasma specimens are emphasised since intraindividual monitoring of LINE-1 hypomethylation may serve as an approach to provide assistance in the observation of AD formation, CRC progression, or remission. Additionally, LINE-1 methylation did not show a significant alteration in IBD. Hence, it might offer a higher screening value.
Our results on the gene expression changes of methylation- and one-carbon metabolism-related enzymes could not explain the decrease of global DNA methylation level; however, they may have complex regulatory mechanisms in the background. We presume that the diminished expression of FOLR2 in colorectal tumours can contribute to the declined methyl-donor availability, which may lead to the process of global DNA hypomethylation.
Availability of data and materials
The microarray dataset analysed during the current study are available in the Gene Expression Omnibus (GEO) data repository (accession number: GSE100179).
Apolipoprotein B mRNA editing enzyme catalytic subunit 1
Apolipoprotein B mRNA editing enzyme catalytic subunit 3B
5-aminoimidazole-4-carboxamide ribonucleotide formyltransferase/IMP cyclohydrolase
DNA methyltransferase 1
DNA methyltransferase 3a
DNA methyltransferase 3b
Fresh frozen tissue
Formalin-fixed paraffin-embedded tissue
- FOLR2 :
Folate receptor beta
- FR-beta :
Folate receptor beta
Phosphoribosylglycinamide formyltransferas, phosphoribosylglycinamide synthetase, phosphoribosylaminoimidazole synthetasee
Gene Expression Omnibus
Human Transcriptome Array
- IBD :
Inflammatory bowel disease
- LINE-1 :
Long interspersed nuclear element 1
Methylenetetrahydrofolate dehydrogenase, cyclohydrolase and formyltetrahydrofolate synthetase 1
Methylenetetrahydrofolate dehydrogenase (NADP + Dependent) 2, methenyltetrahydrofolate cyclohydrolase
- NAT :
Normal adjacent to tumour tissue
Receiver operating characteristic
- SAM :
Serine Hydroxymethyltransferase 2
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49.
Fearon ER, Vogelstein B. A genetic model for colorectal tumorigenesis. Cell. 1990;61(5):759–67.
Greenstein AJ, Sachar DB, Smith H, Janowitz HD, Aufses AH. Patterns of neoplasia in Crohn’s disease and ulcerative colitis. Cancer. 1980;46(2):403–7.
Gillen CD, Walmsley RS, Prior P, Andrews HA, Allan RN. Ulcerative colitis and Crohn’s disease: a comparison of the colorectal cancer risk in extensive colitis. Gut. 1994;35(11):1590.
Ng SC, Shi HY, Hamidi N, Underwood FE, Tang W, Benchimol EI, et al. Worldwide incidence and prevalence of inflammatory bowel disease in the 21st century: a systematic review of population-based studies. Lancet. 2018;390(10114):2769–78.
Rodriguez J, Frigola J, Vendrell E, Risques R-A, Fraga MF, Morales C, et al. Chromosomal instability correlates with genome-wide DNA demethylation in human primary colorectal cancers. Cancer Res. 2006;66(17):8462–9468.
Sahnane N, Magnoli F, Bernasconi B, Tibiletti MG, Romualdi C, Pedroni M, et al. Aberrant DNA methylation profiles of inherited and sporadic colorectal cancer. Clin Epigenetics. 2015;7:131.
Suter CM, Martin DI, Ward RL. Hypomethylation of L1 retrotransposons in colorectal cancer and adjacent normal tissue. Int J Colorectal Dis. 2004;19(2):95–101.
Feinberg AP, Tycko B. The history of cancer epigenetics. Nat Rev Cancer. 2004;4(2):143–53.
Ehrlich M. DNA methylation in cancer: too much, but also too little. Oncogene. 2002;21(35):5400–13.
Rodić N, Burns KH. Long interspersed element-1 (LINE-1): passenger or driver in human neoplasms? PLoS Genet. 2013;9(3): e1003402.
Phokaew C, Kowudtitham S, Subbalekha K, Shuangshoti S, Mutirangura A. LINE-1 methylation patterns of different loci in normal and cancerous cells. Nucleic Acids Res. 2008;36(17):5704–12.
Cordaux R, Batzer MA. The impact of retrotransposons on human genome evolution. Nat Rev Genet. 2009;10(10):691–703.
Grützmann R, Molnar B, Pilarsky C, Habermann JK, Schlag PM, Saeger HD, et al. Sensitive Detection of Colorectal Cancer in Peripheral Blood by Septin 9 DNA Methylation Assay. PLoS ONE. 2008;3(11):3759.
Zandvakili I, Lazaridis KN. Cell-free DNA testing: future applications in gastroenterology and hepatology. Therap Adv Gastroenterol. 2019;12:1756284819841896.
Pathak AK, Bhutani M, Kumar S, Mohan A, Guleria R. Circulating cell-free DNA in plasma/serum of lung cancer patients as a potential screening and prognostic tool. Clin Chem. 2006;52(10):1833–42.
Shen SY, Singhania R, Fehringer G, Chakravarthy A, Roehrl MHA, Chadwick D, et al. Sensitive tumour detection and classification using plasma cell-free DNA methylomes. Nature. 2018;563(7732):579–83.
Nuzzo PV, Berchuck JE, Korthauer K, Spisak S, Nassar AH, Abou Alaiwi S, et al. Detection of renal cell carcinoma using plasma and urine cell-free DNA methylomes. Nat Med. 2020;26(7):1041–3.
Qin Y, Wu CW, Taylor WR, Sawas T, Burger KN, Mahoney DW, et al. Discovery, Validation, and Application of Novel Methylated DNA Markers for Detection of Esophageal Cancer in Plasma. Clin Cancer Res. 2019;25(24):7396–404.
Gormally E, Caboux E, Vineis P, Hainaut P. Circulating free DNA in plasma or serum as biomarker of carcinogenesis: Practical aspects and biological significance. Mutation Research/Reviews in Mutation Research. 2007;635(2–3):105–17.
Terry MB, Delgado-Cruzata L, Vin-Raviv N, Wu HC, Santella RM. DNA methylation in white blood cells. Epigenetics. 2011;6(7):828–37.
Kitkumthorn N, Tuangsintanakul T, Rattanatanyong P, Tiwawech D, Mutirangura A. LINE-1 methylation in the peripheral blood mononuclear cells of cancer patients. Clin Chim Acta. 2012;413(9–10):869–74.
Xiao-Jie L, Hui-Ying X, Qi X, Jiang X, Shi-Jie M. LINE-1 in cancer: multifaceted functions and potential clinical implications. Genet Med. 2016;18(5):431–9.
Hoffmann MJ, Schulz WA. Causes and consequences of DNA hypomethylation in human cancer. Biochem Cell Biol. 2005;83(3):296–321.
Gregory S, Ducker JDR. One-Carbon Metabolism in Health and Disease. Cell Metab. 2017;25(1):27.
Burr NE, Hull MA, Subramanian V. Folic Acid Supplementation May Reduce Colorectal Cancer Risk in Patients With Inflammatory Bowel Disease: A Systematic Review and Meta-Analysis. J Clin Gastroenterol. 2017;51(3):247.
Pufulete M, Al-Ghnaniem R, Leather AJM, Appleby P, Gout S, Terry C, et al. Folate status, genomic DNA hypomethylation, and risk of colorectal adenoma and cancer: a case control study. Gastroenterology. 2003;124(5):1240–8.
Giovannucci E, Stampfer MJ, Colditz GA, Rimm EB, Trichopoulos D, Rosner BA, et al. Folate, Methionine, and Alcohol Intake and Risk of Colorectal Adenoma. J Natl Cancer Inst. 1993;85(11):875–83.
Su LJ, Arab L. Nutritional status of folate and colon cancer risk: evidence from NHANES I epidemiologic follow-up study. Ann Epidemiol. 2001;11(1):65–72.
Fuchs CS, Willett WC, Colditz GA, Hunter DJ, Stampfer MJ, Speizer FE, et al. The influence of folate and multivitamin use on the familial risk of colon cancer in women. Cancer Epidemiol Biomarkers Prev. 2002;11(3):227–34.
Hübner U, Geisel J, Kirsch SH, Kruse V, Bodis M, Klein C, et al. Effect of 1 year B and D vitamin supplementation on LINE-1 repetitive element methylation in older subjects. Clin Chem Lab Med. 2013;51(3):649–55.
Figueiredo JC, Grau MV, Wallace K, Levine AJ, Shen L, Hamdan R, et al. Global DNA hypomethylation (LINE-1) in the normal colon and lifestyle characteristics and dietary and genetic factors. Cancer Epidemiol Biomarkers Prev. 2009;18(4):1041–9.
Pufulete M, Al-Ghnaniem R, Rennie JA, Appleby P, Harris N, Gout S, et al. Influence of folate status on genomic DNA methylation in colonic mucosa of subjects without colorectal adenoma or cancer. Br J Cancer. 2005;92(5):838–42.
Cravo ML, Pinto AG, Chaves P, Cruz JA, Lage P, Leitão CN, et al. Effect of folate supplementation on DNA methylation of rectal mucosa in patients with colonic adenomas: correlation with nutrient intake. Clin Nutr. 1998;17(2):45–9.
Schernhammer ES, Giovannucci E, Kawasaki T, Rosner B, Fuchs CS, Ogino S. Dietary folate, alcohol and B vitamins in relation to LINE-1 hypomethylation in colon cancer. Gut. 2010;59(6):794–9.
Kim Y-I. Role of Folate in Colon Cancer Development and Progression. J Nutr. 2003;133(11):3731S-3739S.
Kalmár A, Nagy ZB, Galamb O, Csabai I, Bodor A, Wichmann B, et al. Genome-wide expression profiling in colorectal cancer focusing on lncRNAs in the adenoma-carcinoma transition. BMC Cancer. 2019;19(1):1059.
Parris TZ, Aziz L, Kovács A, Hajizadeh S, Nemes S, Semaan M, et al. Clinical relevance of breast cancer-related genes as potential biomarkers for oral squamous cell carcinoma. BMC Cancer. 2014;14:324.
Valcz G, Galamb O, Krenács T, Spisák S, Kalmár A, Patai ÁV, et al. Exosomes in colorectal carcinoma formation: ALIX under the magnifying glass. Mod Pathol. 2016;29(8):928–38.
Chalitchagorn K, Shuangshoti S, Hourpai N, Kongruttanachok N, Tangkijvanich P, Thong-ngam D, et al. Distinctive pattern of LINE-1 methylation level in normal tissues and the association with carcinogenesis. Oncogene. 2004;23(54):8841–6.
Estécio MRH, Gharibyan V, Shen L, Ibrahim AEK, Doshi K, He R, et al. LINE-1 hypomethylation in cancer is highly variable and inversely correlated with microsatellite instability. PLoS ONE. 2007;2(5): e399.
Jiang A, Sutherland S, Buckingham L, Melson J. Correlation of LINE-1 hypomethylation with size and dysplasia in colorectal tubular adenomas, and risk for synchronous and metachronous CRC. J Clin Oncol. 2020;38(15_suppl):e16083–e16083.
Jiang AC, Buckingham L, Barbanera W, Korang AY, Bishehsari F, Melson J. Erratum to: LINE-1 is preferentially hypomethylated within adenomatous polyps in the presence of synchronous colorectal cancer. Clin Epigenetics. 2017;9:76.
Muto T, Bussey HJR, Morson BC. The evolution of cancer of the colon and rectum. Cancer. 1975;36(6):2251–70.
Wilson AS, Power BE, Molloy PL. DNA hypomethylation and human diseases. Biochimica et Biophysica Acta (BBA) - Reviews on Cancer. 2007;1775(1):138–62.
Sunami E, de Maat M, Vu A, Turner RR, Hoon DSB. LINE-1 hypomethylation during primary colon cancer progression. PLoS ONE. 2011;6(4): e18884.
Lakatos P-L, Lakatos L. Risk for colorectal cancer in ulcerative colitis: changes, causes and management strategies. World J Gastroenterol. 2008;14(25):3937–47.
Bernstein CN, Blanchard JF, Kliewer E, Wajda A. Cancer risk in patients with inflammatory bowel disease. Cancer. 2001;91(4):854–62.
Glória L, Cravo M, Pinto A, de Sousa LS, Chaves P, Leitão CN, et al. DNA hypomethylation and proliferative activity are increased in the rectal mucosa of patients with long-standing ulcerative colitis. Cancer. 1996;78(11):2300–6.
Kupcinskas J, Steponaitiene R, Langner C, Smailyte G, Skieceviciene J, Kupcinskas L, et al. LINE-1 hypomethylation is not a common event in preneoplastic stages of gastric carcinogenesis. Sci Rep. 2017;7(1):4828.
Nagai Y, Sunami E, Yamamoto Y, Hata K, Okada S, Murono K, et al. LINE-1 hypomethylation status of circulating cell-free DNA in plasma as a biomarker for colorectal cancer. Oncotarget. 2017;8(7):11906.
Rai K, Sarkar S, Broadbent TJ, Voas M, Grossmann KF, Nadauld LD, et al. DNA demethylase activity maintains intestinal cells in an undifferentiated state following loss of APC. Cell. 2010;142(6):930–42.
Conticello SG. The AID/APOBEC family of nucleic acid mutators. Genome Biol. 2008;9(6):229.
Tahiliani M, Koh KP, Shen Y, Pastor WA, Bandukwala H, Brudno Y, et al. Conversion of 5-Methylcytosine to 5-Hydroxymethylcytosine in Mammalian DNA by MLL Partner TET1. Science. 2009;324(5929):930–5.
Ito S, Shen L, Dai Q, Wu SC, Collins LB, Swenberg JA, et al. Tet Proteins Can Convert 5-Methylcytosine to 5-Formylcytosine and 5-Carboxylcytosine. Science. 2011;333(6047):1300–3.
He Y-F, Li B-Z, Li Z, Liu P, Wang Y, Tang Q, et al. Tet-Mediated Formation of 5-Carboxylcytosine and Its Excision by TDG in Mammalian DNA. Science. 2011;333(6047):1303–7.
Lee RM, Hirano K, Anant S, Baunoch D, Davidson NO. An alternatively spliced form of apobec-1 messenger RNA is overexpressed in human colon cancer. Gastroenterology. 1998;115(5):1096–103.
Greeve J, Lellek H, Apostel F, Hundoegger K, Barialai A, Kirsten R, et al. Absence of APOBEC-1 mediated mRNA editing in human carcinomas. Oncogene. 1999;18(46):6357–66.
Burns MB, Temiz NA, Harris RS. Evidence for APOBEC3B mutagenesis in multiple human cancers. Nat Genet. 2013;45(9):977–83.
Stenglein MD, Harris RS. APOBEC3B and APOBEC3F inhibit L1 retrotransposition by a DNA deamination-independent mechanism. J Biol Chem. 2006;281(25):16837–41.
Kinomoto M, Kanno T, Shimura M, Ishizaka Y, Kojima A, Kurata T, et al. All APOBEC3 family proteins differentially inhibit LINE-1 retrotransposition. Nucleic Acids Res. 2007;35(9):2955–64.
Mehrmohamadi M, Liu X, Shestov AA, Locasale JW. Characterization of the usage of the serine metabolic network in human cancer. Cell Rep. 2014;9(4):1507–19.
Henderson GI, Perez T, Schenker S, Mackins J, Antony AC. Maternal-to-fetal transfer of 5-methyltetrahydrofolate by the perfused human placental cotyledon: evidence for a concentrative role by placental folate receptors in fetal folate delivery. J Lab Clin Med. 1995;126(2):184–203.
Hermann A, Gowher H, Jeltsch A. Biochemistry and biology of mammalian DNA methyltransferases. Cell Mol Life Sci. 2004;61(19–20):2571–87.
Robert M-F, Morin S, Beaulieu N, Gauthier F, Chute IC, Barsalou A, et al. DNMT1 is required to maintain CpG methylation and aberrant gene silencing in human cancer cells. Nat Genet. 2003;33(1):61–5.
Hervouet E, Vallette FM, Cartron P-F. Dnmt3/transcription factor interactions as crucial players in targeted DNA methylation. Epigenetics. 2009;4(7):487–99.
Yang J, Deng C, Hemati N, Hanash SM, Richardson BC. Effect of mitogenic stimulation and DNA methylation on human T cell DNA methyltransferase expression and activity. J Immunol. 1997;159(3):1303–9.
Christman JK, Sheikhnejad G, Dizik M, Abileah S, Wainfan E. Reversibility of changes in nucleic acid methylation and gene expression induced in rat liver by severe dietary methyl deficiency. Carcinogenesis. 1993;14(4):551–7.
Slack A, Cervoni N, Pinard M, Szyf M. Feedback regulation of DNA methyltransferase gene expression by methylation. Eur J Biochem. 1999;264(1):191–9.
de Boer E, Crane LMA, van Oosten M, van der Vegt B, van der Sluis T, Kooijman P, et al. Folate Receptor-Beta Has Limited Value for Fluorescent Imaging in Ovarian, Breast and Colorectal Cancer. PLoS ONE. 2015;10(8): e0135012.
Alonso-Aperte E, González MP, Póo-Prieto R, Varela-Moreiras G. Folate status and S-adenosylmethionine/S-adenosylhomocysteine ratio in colorectal adenocarcinoma in humans. Eur J Clin Nutr. 2008;62(2):295–8.
Abbasi IHR, Abbasi F, Wang L, Abd El Hack ME, Swelum AA, Hao R, et al. Folate promotes S-adenosyl methionine reactions and the microbial methylation cycle and boosts ruminants production and reproduction. AMB Expres. 2018;8(1):65.
Sibani S, Melnyk S, Pogribny IP, Wang W, Hiou-Tim F, Deng L, et al. Studies of methionine cycle intermediates (SAM, SAH), DNA methylation and the impact of folate deficiency on tumor numbers in Min mice. Carcinogenesis. 2002;23(1):61–5.
Mahmoud AM, Ali MM. Methyl Donor Micronutrients that Modify DNA Methylation and Cancer Outcome. Nutrients. 2019;11(3):608.
We kindly thank Gitta Szabó for her statistical advices, Hajnalka Kovács and Anikó Szigeti for sample collection, and William J Kothalawala for the helpful suggestions. Moreover, we are grateful to Kitti Erdő-Bonyár English language editor, and William J Kothalawala English expert from Semmelweis University, for the proofreading of our manuscript.
Open access funding provided by Semmelweis University. This research was funded by the National Research, Development and Innovation Office, Hungary (grant number 2020–1.1.2-PIACI-KFI-2021–00298). The study was also financed by the Higher Education Institutional Excellence Program of the Ministry of Human Capacities in Hungary, within the framework of the molecular biology thematic program of the Semmelweis University.
Ethics approval and consent to participate
The study has been approved by the local ethics committee (Regional and Institutional Committee of Science and Research Ethics; TUKEB Nr: 14383–2/2017/EKU), all applied methods were performed in accordance with the relevant guidelines and regulations. Written informed consent was obtained from all the patients prior to sample collection.
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Summary of LINE-1 methylation levels, clinicopathological and demographic data of the investigated patients
LINE-1 methylation in NAT, tumour and paired NAT specimens, AD subtypes, and CRC stages. Significant LINE-1 methylation reduction was noticed in the NAT vs. N comparison (A) and in tumour specimens compared to their paired NAT tissue biopsies (B) (*p≤0.05). In TA, TVA vs. N comparison, significantly lower LINE-1 methylation was found (*p≤0.05), and a decreasing trend was observed in TVA compared to TA biopsies (C, left). Also, significant LINE-1 hypomethylation was detected in both early (Astler-Coller modified Dukes’ A-B) and late (Astler-Coller modified Dukes’ C-D) carcinomas compared to N specimens (C, right) (*p≤0.05). N: healthy, AD NAT: normal adjacent to adenoma tissue, CRC NAT: normal adjacent to carcinoma tissue, AD: colorectal adenoma, CRC: colorectal carcinoma, TA: tubular adenoma, TVA: tubulovillous adenoma
mRNA expression alterations of DNMT1, APOBEC1, APOBEC3B genes in tumorous compared to normal samples. On the heatmap, distinct colours are coupled to different expression intensity values: green - low, black - intermediate, red - high intensity. Each row represents different genes, and each column indicates the investigated sample types. DNMT1 mRNA level increased significantly only in CRC vs. N comparison, while APOBEC1 significantly elevated, and APOBEC3B significantly decreased only in AD samples compared to healthy controls (p≤0.05).
Expression changes of DNMT enzymes in CRC tissue sections. A) Significant increase of DNMT1 labelling in the epithelial cells (*p≤0.05) and no altered expression in stromal cells were detected in the cancerous area compared to NAT. There were no significant changes in DNMT3a (B) and DNMT3b (C) staining in CRC vs. NAT comparison. Expression differences were illustrated with the Q-score method (below). Scale bars on the top: 1000µm, on the bottom: 20µm. NAT: normal adjacent to tumour tissue, CRC: colorectal carcinoma, ep: epithelial cells, st: stromal cells, Q-score: Quick-score method
Immunohistochemistry of FR-beta encoded by FOLR2 in AD and CRC samples. Significant reduction of FR-beta intensity was found in epithelial and stromal cells in the AD (A) and CRC (B) vs. NAT comparison (**p≤0.01) except for CRC stromal cells, where a slight decrease was observed (B). Histological changes are represented with the Q-score method on the right. Scale bars on the left: 50µm, on the right: 20µm. NAT: normal adjacent to tumour tissue, AD: colorectal adenoma, CRC: colorectal carcinoma, ep: epithelial cells, st: stromal cells, Q score: Quick-score method
About this article
Cite this article
Szigeti, K.A., Kalmár, A., Galamb, O. et al. Global DNA hypomethylation of colorectal tumours detected in tissue and liquid biopsies may be related to decreased methyl-donor content. BMC Cancer 22, 605 (2022). https://doi.org/10.1186/s12885-022-09659-1
- Colorectal cancer
- DNA methylation
- Liquid biopsy
- Folic acid