Open Access

Aberrant promoter methylation status is associated with upregulation of the E2F4 gene in breast cancer

  • Authors:
    • Farman Ullah Farman
    • Farhan Haq
    • Noor Muhammad
    • Nawab Ali
    • Hazir Rahman
    • Muhammad Saeed
  • View Affiliations

  • Published online on: March 29, 2018     https://doi.org/10.3892/ol.2018.8382
  • Pages: 8461-8469
  • Copyright: © Farman et al. This is an open access article distributed under the terms of Creative Commons Attribution License.

Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

E2F4 is an important basal transcription factor with the potential to promote tumor growth. Its upregulation in various types of cancer has been linked to numerous genetic factors; however, the nature of the involvement of epigenetic mechanisms, including DNA methylation, remains elusive. In the present study, E2F4 expression profiles were determined in 100 paired breast tumor and control samples, through RT‑qPCR using the SYBR® green method. Furthermore, the E2F4 promoter methylation status in each of these samples was assessed using methylation specific PCR, in order to evaluate its impact on gene expression. A two‑fold increase in E2F4 gene expression was observed in the breast tumors compared with in their respective controls (P=0.022); of these tumors, ~72% were under‑methylated. The change in methylation status was also significantly higher (P<0.001) in the tumor samples. Methylation status was negatively correlated (r=‑30) with E2F4 expression profiles, indicating that a decrease in methylation may promote higher expression of E2F4. The two study cohorts (>45 and ≤45 years) had comparable methylation profiles, though they had significantly decreased methylation status compared with controls. Various histo‑pathological types also have different methylation profiles, indicating the presence of a tissue specific methylation signature. The results of the present study demonstrated that E2F4 methylation status can have a notable influence on its expression, and that it may have prognostic value in breast carcinogenesis.

Introduction

Cell cycle regulation is critical to normal growth and division. Any deregulation of activity in the cell cycle regulatory proteins may proceed to uncontrolled division and tumor formation. Cell cycle checkpoints, particularly the G1-to-S checkpoint, serve crucial roles in the initiation and progression of cancer, due to the availability of checkpoint regulatory proteins, including E2F basal transcription factors and retinoblastoma (Rb) proteins (1).

The E2F transcription factors constitute a superfamily of basal elongation factors, which is comprised of eight members (2). Based on the regulatory roles that these proteins have, they are broadly categorized as activators (E2F1-E2F3) and repressors (E2F4-E2F8) (35). The association of pocket proteins (pRb, Rbl1/p107, Rbl2/p130) with E2F proteins abrogates their access to target promoters, hence hindering entry to the G1-to-S (6,7). The release of E2Fs from these pocket proteins is necessary for the transcription of genes that are required to cross this checkpoint (8,9). It has been previously reported that the regulatory roles of E2F1-E2F3 in differentiating cells may switch those of repressor proteins (10). The expression profiles of E2F1-E2F3 are strictly regulated in a cell cycle stage-specific manner; however, E2F4 and E2F5 have been demonstrated to maintain constant expression levels at all stages of the cell cycle (3), which indicates that their functions are not limited to any particular stage (11).

Studies researching the involvement of E2F4 in tumorigenesis first reported that a mutated form of this gene was identifiable in several types of cancer, including gastrointestinal and prostate cancer (1214). Elevated expression levels of E2F4 have been associated with gastric cancer (1), solid osteosarcoma (15) and breast tumor progression (16); whereas E2F4-deficient mice exhibited severe defects in hematopoietic development (17). The E2F4 deficient cells have also been identified to easily undergo apoptosis upon DNA damage, possibly due to the absence of E2F4/DP1/p130 complexes (18). Therefore, it can be concluded that the elevated E2F4 gene levels can be attributed to frequent genetic alterations (13) and abnormal splicing (19). The epigenetic mechanisms, including promoter DNA methylation and histone modifications influencing E2F4 expression, remain elusive; although, the involvement of epigenetic mechanisms in cancer and metastasis, specifically in regulating the expression patterns of various genes, is well known (20,21).

Promoter DNA methylation can produce effects similar to gene mutation in order to achieve gain or loss of function of certain genes, leading to carcinogenesis (22). Furthermore, aberrant promoter methylation has been identified in various types of tumor cells, and thus has been theorized to be a biomarker for the early prognosis of cancer (23). In previous studies, aberrant methylation was identified at the pRb (24), Rbl2/p130 (25) and E2F5 promoters (26), and was significantly associated with their expression profiles in breast as well as head and neck cancer.

The present study was designed to investigate the expression levels of the E2F4 mRNA transcript and its promoter methylation status in breast tumor tissues, as well as in their adjacent normal control (ANCT) samples. This allowed an evaluation of whether promoter methylation is a decisive factor for E2F4 gene transcription. Associations between promoter methylation status and with various clinico-pathological outcomes of breast tumors were also established to evaluate the prognostic value of E2F4 promoter methylation in breast carcinogenesis within a clinical setting.

Materials and methods

Recruitment of patients with breast cancer, and specimen collection

A total of 100 pairs of tumor and ANCT tissue samples were collected and stored in RNA Later™ stabilization solution (cat. no. AM7024; Thermo Fisher Scientific, Inc., Waltham, MA, USA) at the time of surgery from patients with breast cancer at Lady Reading Hospital (Khyber Pakhtunkhwa, Pakistan) between June 2012 and August 2015. All the patients were female. The ANCT sections were taken from an area ~2 cm from the cancer lesion site, as selected by a histo-pathologist. The patients with a history of other infectious, metabolic or familial diseases were excluded from the present study. Furthermore, prior approval was obtained from the Ethical Review Committee (approval no. CIIT-09-10-14) of the COMSTAS Institute of Information Technology (Islamabad, Pakistan) and the collaborating hospital aforementioned for the initiation of this project. In addition, informed and written consent was received from all patients enrolled in the present study, prior to sample collection.

mRNA extraction and real-time qPCR analysis

RNA extraction from all tissue samples (tumor and ANCT) was performed using the standard TRIzol® method (27). The extracted mRNA was converted to cDNA using a cDNA synthesis kit (Thermo Fisher Scientific, Inc.). The primers used for mRNA expression analysis were as follows: E2F4 forward, 5′-TCAGAAATCTTTGATCCCAC-3′ and reverse, 5′-AGATATAATCGTGGTCTCCC-3′; β-actin (internal control) forward, 5′-CACTCTTCCAGCCTTCCTTC-3′ and reverse, 5′-TGATCTCCTTCTGCATCGTG-3′. PCR amplification was performed using SYBR®-Green method (28) through Step One Plus Real-Time PCR system (Applied Biosystems; Thermo Fisher Scientific, Inc.). The thermocycler conditions were as follows: Initial denaturation at 95°C for 10 min, followed by 40 cycles of 95°C for 15 sec, 56°C for 60 sec and 72°C for 20 sec; and a final extension at 72°C for 1 min. mRNA expression of the E2F4 gene was calculated using the 2−ΔΔCq method (29).

Extraction and bisulfite modification of DNA

DNA from tumor and ANCT tissues was extracted using the standard phenol-chloroform method (30), and confirmed by 2% agarose gel electrophoresis. Genomic DNA (~2–5 µg) was bisulfite modified using an EpiJET™ Bisulfite Conversion kit (Thermo Fisher Scientific, Inc.), as per the manufacturer's instructions. Bisulfite treatment was used to convert non-methylated cytosine ‘C’ bases to thymine ‘T’ bases; however, the methylated cytosine ‘C’ was left unchanged.

Primer designing for methylation analysis

The methylation status of the E2F4 promoter was analyzed in the tumor and ANCT tissues from all recruited patients with breast cancer. The target region is located at a site 84 bp upstream and 117 bp downstream from the transcription start site (TSS). Furthermore, methylation-specific primers were designed for the target site, according to National Center for Biotechnology Information Nucleotide database (https://www.ncbi.nlm.nih.gov/nucleotide?cmd=search) and scanned 29 CpG sites within the targeted sequence. The two sets of methylated primers were as follows: Forward, 5′-CGTTACGTTTTTTGGAAGGC-3′ and reverse, 5′-CGTACCGACTTAAAATACCCG-3′. Un-methylated primers were as follows: Forward, 5′-TAGTGTTATGTTTTTTGGAAGGTGT-3′ and reverse, 5′-TTCATACCAACTTAAAATACCCAAA-3′. Promoter methylation status was analyzed using the MethPrimer™ software (version 1.0; http://www.urogene.org/cgi-bin/methprimer/methprimer.cgi; Peking Union Medical Colege Hospital, Beijing, China) available online (31).

Methylation-specific PCR (MSP)

The targeted promoter regions around the TSS of E2F4, amplified from bisulfite converted genomic DNA, was performed using MSP reactions. Bisulfite converted DNA was amplified in a one-step MSP reaction, using two independent sets of primers aforementioned (methylated and un-methylated) and Maxima® Hot Start Taq DNA polymerase (Thermo Fisher Scientific, Inc.). Briefly, bisulfite converted DNA (~50 ng) was amplified with 0.2 µM of each deoxynucleotide triphosphate and primer (both methylated and un-methylated forward and reverse primers) using 1 U Maxima® Hot Start Taq DNA polymerase.

The PCR conditions for E2F4 were as follows: A single cycle at 98°C for 30 sec (initial denaturation); followed by 35 cycles at 98°C for 10 sec (denaturation), 62°C for 30 sec (annealing) and 72°C for 30 sec (extension); and a single cycle at 72°C for 7 min (final extension) using 1 U Maxima® Hot Start Taq DNA polymerase. MSP validity was checked by amplifying CpG-methylated human genomic DNA (Thermo Fisher Scientific, Inc.) as the positive control with PCR water (Solis BioDyne, Tartu, Estonia) used as the negative control; however, for the non-methylation-specific PCR, bisulfite-unconverted human genomic DNA was used as the positive control for all primers specific to the E2F4 gene promoter. All PCR reactions were carried out using a BIOER™ Thermal Cycler 9500 (Hangzhou Bioer Technology Co., Ltd., Binjiang, China).

Gel electrophoresis and quantification

The PCR products from the aforementioned reactions were resolved on a 2% agarose gel and visualized under an UV illuminator (BioDoc Analyzer, Biometra; Analytik Jena, Jena, Germany) using a 100 bp DNA size ladder. To establish the degree of methylation from qualitative MSP data, the change in methylation (Δ meth) value was calculated, as previously described (25,26). The Δ meth value was calculated by subtracting the un-methylation values of a particular sample from its methylation values, when considering un-methylation as the default condition.

Statistical analysis

Statistical analyses were performed using OriginPro 2015 software (OriginLab Corporation, Northampton, MA, USA). The results are expressed as the mean and (±) standard error of the mean to determine the descriptive statistics for qualitative data. For the analysis of transcript profiles, the expression data of the target gene (E2F4) were normalized against those of the internal control gene β-actin. The correlations between various factors were assessed using the Pearson's correlation coefficient. Depending on the experiment, the statistical significance was determined to 95% confidence intervals using Pearson's correlation coefficient, Kruskal-Wallis and Mann Whitney tests and one-way analysis of variance. P<0.05 was considered to indicate a statistically significant difference. Statistical analysis enabled an investigation of the correlations between the methylation frequency and expression levels of the E2F4 gene.

Results

Promoter methylation analysis of E2F4 for patients with breast cancer

The promoter methylation status of the E2F4 gene was analyzed through an MSP approach. Of the breast tumors that were observed, ~71.5% were un-methylated at the E2F4 promoter; whereas ~21.5% of breast tumors were determined to be exclusively methylated at this promoter (Fig. 1 and Table I). Conversely, ANCT control tissues were identified as differentially methylated (~51.5% methylated and ~48.5% un-methylated) at the E2F4 promoter. The Δ meth value also exhibited a remarkable difference between the breast tumor (mean=3,450.76) and the ANCT (mean=722.29) tissue samples (Fig. 2A; Table I). This indicates a statistically significant (P<0.001) decrease in E2F4 promoter methylation among patients with breast cancer, as compared with the controls. These findings highlight a possible role for E2F4 promoter methylation in breast cancer prognosis, via its elevated mRNA expression levels (Tables II and III; Fig. 2A).

Table I.

Statistical analysis of E2F4 mRNA transcript expression and its promoter methylation status in various breast cancer cohorts.

Table I.

Statistical analysis of E2F4 mRNA transcript expression and its promoter methylation status in various breast cancer cohorts.

ANCTTumors


95% CI 95% CI Control vs. tumor





CharacteristicsMeanSEMUpperLowerSkewnessMeanSEMUpperLowerSkewnessP-value
mRNA expression
  Overall1.000.562.15−0.154.891.780.873.57−0.014.090.022a
Methylation status overall
  Un-meth3,398.61745.614,937.471,859.754.594,087.59375.314,862.193,312.980.11 <0.001a
  Meth2,676.32685.994,092.121,260.511.02636.82275.641,205.7267.932.67
  Δ meth722.291,074.182,939.28−1,494.701.933,450.76407.764,292.342,609.19−0.28
Age (years) ∆ meth
  ≤45509.88235.735,754.31−4,734.561.863,589.46600.744,927.992,250.930.680.005a
  >45889.19671.322,339.49−561.12−1.743,341.79572.664,578.942,104.63−0.85
Age (years) of menarche ∆ meth
  ≤12 (early)3,016.731,617.466,540.87−507.422.454,025.39440.894,986.023,064.760.510.010a
  >12 (late)−1,763.351,035.76516.35−4,043.05−0.172,828.25676.964,318.241,338.27−0.09
Menopausal stage ∆ meth
  Pre-menopausal−1,102.75116.641,456.19−3,661.69−0.303,572.11676.345,060.722,083.51−0.060.010a
  Post-menopausal2,406.95167.756,053.72−1,239.832.493,338.75500.264,428.712,248.79−0.90
∆ meth status of histopathological tumor type
  IDC2,018.931,446.915,102.96−1,065.092.214,010.891,484.535,043.652,978.13−0.480.001a
  ILC−2,387.76444.473,010.75−7,786.27−0.031,835.901,018.744,664.27−992.470.73
  DCIS−576.721,638.794,638.67−5,792.13−0.483,228.821,605.715,156.471,301.180.07
Tumor stage ∆ meth
  I−1,227.12855.83587.15−3,041.39−0.433,372.69519.794,474.622,270.78−0.54 <0.001a
  II6,203.493,652.7116,344.66−3,937.692.232,841.46667.114,693.58989.341.34
  III2,633.60278.283,830.961,436.24−1.734,908.641,226.9210,187.65−370.36−0.19

{ label (or @symbol) needed for fn[@id='tfn1-ol-0-0-8382'] } ∆ meth, change in methylation; ANCT, adjacent normal control tissues; un-meth, un-methylation; meth, methylation; SEM, standard error of the mean (±); CI, confidence interval

a P<0.05; IDC, invasive ductal carcinoma; ILC, invasive lobular carcinoma; DCIS, ductal carcinoma in situ.

Table II.

Correlation analysis of the E2F4 gene mRNA transcript expression and promoter methylation for various clinico-pathological factors in patients with breast cancer.

Table II.

Correlation analysis of the E2F4 gene mRNA transcript expression and promoter methylation for various clinico-pathological factors in patients with breast cancer.

ANCTTumors


CharacteristicsPearson's rAdj. R2F valueProb>FPearson's rAdj. R2F valueProb>F
mRNA expression
  Overall−0.07−0.040.120.74−0.240.021.410.25
Methylation status overall
  Un-meth−0.05−0.040.070.790.570.2911.090.003
  Meth−0.10−0.030.240.63−0.06−0.040.070.79
  Δ Meth0.03−0.040.020.890.560.2910.660.003
Age (years) ∆ meth
  ≤45−0.07−0.110.040.850.680.407.700.02
  >450.26−0.010.900.360.460.143.160.10
Age (years) of menarche ∆ meth
  ≤12 (early)−0.25−0.020.830.380.640.367.810.02
  >12 (late)0.400.081.940.190.610.315.930.04
Menopausal stage ∆ meth
  Pre-menopausal0.20−0.060.410.540.780.5816.590.01
  Post-menopausal−0.17−0.060.360.560.320.021.230.29
∆ meth status of histopathological tumor types
  IDC−0.11−0.060.160.700.640.369.570.01
  ILC−0.13−0.310.050.830.490.010.950.40
  DCIS0.22−0.430.100.77−0.220.430.100.77
Tumor stage ∆ meth
  I0.12−0.050.220.640.600.318.040.01
  II−0.36−0.160.460.540.35−0.170.420.56
  III0.860.492.900.330.53−0.430.390.64

[i] ∆ meth, change in methylation; ANCT, adjacent normal controls tissues; un-meth, un-methylation; meth, methylation; IDC, invasive ductal carcinoma; ILC, invasive lobular carcinoma; DCIS, ductal carcinoma in situ; prob, probability.

Table III.

Promoter methylation frequencies of the E2F4 gene among various clinico-pathological parameters in breast cancer.

Table III.

Promoter methylation frequencies of the E2F4 gene among various clinico-pathological parameters in breast cancer.

ANCTTumors


CharacteristicsMeth no. (%)Un-meth no. (%)Meth no. (%)Un-meth no. (%)
Overall103 (51.5)97 (48.5)57 (28.5)143 (71.5)
Age (years)
  ≤4545 (53.57)39 (46.43)28 (33.33)56 (66.67)
  >4558 (50.0)58 (50.0)29 (25.0)87 (75.0)
Age (years) of menarche
  ≤12 (early)53 (43.44)69 (56.56)20 (16.39)102 (83.61)
  >12 (late)50 (64.10)28 (35.6)37 (47.44)41 (52.56)
Menopausal status
  Pre-menopausal80 (62.5)48 (37.5)46 (35.94)82 (64.06)
  Post-menopausal23 (31.94)49 (68.06)11 (15.28)61 (84.72)
Histo-pathological types
  IDC63 (45.65)75 (54.35)35 (25.36)103 (74.64)
  ILC31 (64.58)17 (35.42)19 (39.58)29 (60.42)
  DCIS9 (64.29)5 (35.71)3 (21.43)11 (78.57)
Stages
  I61 (48.41)65 (51.59)40 (31.75)86 (68.25)
  II23 (57.5)17 (42.5)10 (25.0)30 (75.0)
  III19 (55.88)15 (44.12)7 (20.59)27 (79.41)

[i] IDC, invasive ductal carcinoma; ILC, invasive lobular carcinoma; DCIS, ductal carcinoma in situ; ANCT, adjacent normal control tissues; un-meth, un-methylation; meth, methylation.

Promoter methylation status of the E2F4 gene among distinct clinicopathological breast cancer cohorts

The promoter methylation status of the E2F4 gene was analyzed in two age-based study cohorts, consisting of patient's ≤45 and >45 years. The study cohort comprising patients ≤45 years of age, exhibited hypo-methylation in ~66.67% of tumor; whereas, ~46.43% of their ANCT controls were hypo-methylated, indicating that non-cancerous tissues have higher methylation levels at the E2F4 promoter (Table I). Similarly, patients in the >45 years cohort exhibited hypo-methylation in ~75% of tumors, as compared with in ~50% of control samples, which again indicated higher methylation levels in controls. In addition, the Δ meth values between each study cohort were identified to be significantly different (P=0.005). Furthermore, the Pearson's correlation coefficient values for Δ meth revealed a strong positive correlation (r=0.68 for >45 years and r=0.46 for ≤45 years) for each study cohort (Tables IIII and Fig. 2B). These findings indicated an age-independent trend in E2F4 promoter methylation status that may have implications for disease prognosis. Similar methylation profiles were also observed in pre- and post-menopausal study cohorts (Table I and Fig. 2C) with 45 years being the mean age for these cohorts. This indicated that the variable methylation signature in the cohorts may be a consequence of a hormonal imbalance in these patients, revealed a mechanistic insight into the whole process of DNA methylation. Patients with and early (≤12 years) or late (>12 years) age of menarche were observed to have statistically significant (P=0.010) variations in their promoter methylation status (Table I and Fig. 2D); the results demonstrated that there was a strong positive correlation for early (r=0.64) and late (r=0.61) age of menarche patients (Table II). This indicates that early age of menarche is a putative risk that may result in notable methylation variations.

Similarly, promoter methylation status among tumors at different disease stages revealed a gradual reduction in methylation frequency with advanced tumor stage, i.e. from SI to SIII (Table III; Fig. 2E). As the Δ meth was revealed to be statistically significant (P<0.001) among these tumors, it was suggested that promoter methylation had a notable implication for disease progression and prognosis. Furthermore, for the various histopathological types of breast cancer, including invasive ductal carcinoma, invasive lobular carcinoma and ductal carcinoma in situ, the Δ meth value was statistically significant (P=0.001; Tables I and III; Fig. 2F). This demonstrates that methylation may be controlled in a tissue specific manner.

Quantitative PCR (qPCR) mRNA transcript analysis of E2F4

Quantitative mRNA transcript analysis of the E2F4 gene was performed using qPCR with all biopsy samples from the enrolled patients with breast cancer. β-actin was used as the internal standard for gene amplification. It was identified that there was an almost two-fold elevation in E2F4 gene expression in tumor samples when compared with ANCT tissues (Fig. 3); this difference was statistically significant (P=0.022; Table I) and indicated a potential involvement of E2F4 expression in tumor pathogenesis. The upregulation of E2F4 in breast tumor samples was determined to be negatively correlated with its promoter methylation, as per the Pearson's correlation test. The value of Pearson's ‘r’ (Pearson's r=−0.24) indicated a negative correlation with disease outcome, which suggested that a higher E2F4 expression would lead to a poorer outcome. Notably, no correlation (Pearson's r=−0.07) was calculated for the control samples (Table II).

Expression of E2F4 correlates negatively with promoter methylation

The E2F4 promoter methylation status may correspond to its expression; the degree of E2F4 promoter methylation has a strong positive effect on its transcript expression. Statistical analysis revealed a moderately negative correlation between promoter methylation (Δ meth) and gene expression (r=−0.30) in breast tumor tissues, as well as in the ANCT samples (r=−0.32; Table IV). This demonstrates that methylation in these samples may be the causative factor for their gene upregulation. Furthermore, there was a moderately positive association between the methylation status of controls and relative expression of E2F4 gene (r=0.41); whereas, there was no association observed between tumor methylation and their relative expression (r=−0.03), although, a moderately negative correlation (r=−0.34) between the un-methylation status of these samples and their expression was identified. This demonstrates an involvement of the hypo-methylated status of E2F4 in regulating its expression in breast tumor samples (Table IV).

Table IV.

Pearson's correlation between E2F4 gene expression and the promoter methylation status in patients with breast cancer.

Table IV.

Pearson's correlation between E2F4 gene expression and the promoter methylation status in patients with breast cancer.

ANCTTumor


Pearson's correlationRelative expressionUn-methMeth∆ MethRelative expressionUn-methMeth∆ Meth
Control relative expression1−0.090.41a−0.32−0.080.330.81a−0.25
Diseased relative expression−0.080.89a−0.090.67a1−0.34−0.03−0.30
Diseased un-meth0.32−0.160.12−0.19−0.3410.240.76a
Diseased meth0.81a−0.130.67a−0.52a−0.030.241−0.45a
Diseased ∆ meth−0.25−0.06−0.350.18−0.300.76a−0.45a1
Control un-meth−0.091−0.120.77a0.89a−0.16−0.13−0.06
Control meth0.41a−0.121−0.73a−0.090.120.67a−0.35
Control ∆ meth−0.320.77a−0.73a10.67a−0.19−0.52a0.18

a Strong correlations. Un-meth, un-methylation; meth, methylation; Δ meth, change in methylation; ANCT, adjacent normal control tissues.

Discussion

E2F4 being a basal transcription factor serves vital roles in cell cycle regulation, and its deregulated expression may lead to un-controlled growth and tumorigenesis (32). The transcriptional activity of E2F4 requires hetero-dimerization with its DP partner, which selectively binds to the CpG rich DNA sequence TTTC/GC/GCGC/G (33,34), regulating downstream sequences. Furthermore, the transcriptional domain of E2F4 is occupied by members of the RB family of proteins (pRB, p130/RBL2, p107/RBL1) (35,36); therefore, making a cyclin E repressor complex (CERC) (37), which causes the E2F4/DP complex to lose the ability to initiate transcription at the required promoters. As a result, the aberrant expression of E2F4 creates an unbalance in the normal composition of the cell cycle complex, causing uncontrolled cell division and tumor formation (26). High expression of E2F4 has previously been identified in various types of cancer, including gastric (1), prostate (5) and breast cancer (16); however, the key factors regulating E2F4 gene expression remain elusive.

The results of the present study support those of prior studies, which reported upregulation of the E2F4 gene in breast tumor tissues (mean=1.78) compared with in the controls (mean=1.00). Furthermore, the present study demonstrated that promoter methylation status is a primary cause of upregulation. The findings of the present study indicated that the frequency of E2F4 promoter methylation was notably low in the breast tumor tissues (28.5% of all cases), as compared with in the respective ANCT samples (51.5%). Furthermore, the Δ meth value was identified to be significantly different (P<0.001) between the control and tumor tissues. Promoter methylation is an epigenetic phenomenon that is crucial to activating or deactivating the regulatory sequence of a gene, and, therefore, to controlling the gene expression and effecting the normal growth pattern (21,22). This epigenetic phenomenon has already been demonstrated to serve an important role in influencing the expression profiles in a number of growth associated genes, including members of the RB family (pRB, Rbl2) (24,25), as well as other E2Fs (25,26). The CpG rich regions can be modified through methylation, which creates a competitive environment for E2F4 binding in these regions, even for the regulatory elements of E2F4 (38).

The elevated mRNA expression of E2F4 has a moderately negative association with its aberrant methylation level (r=−0.30) in cancerous tissues, and may serve a prognostic role. This indicates that un-methylation is associated with transcribing promoters; however, this must be further evaluated in other cancer types. The hypo-methylated promoters were repeatedly identified in all study cohorts, based on various factors, including patient age, the age of menarche, menopausal status, tumor histo-pathological types, and tumor stages. The decrease in promoter methylation highlights the prognostic involvement of E2F4 in tumor progression, and may be useful as a predictor for the early diagnosis of breast cancer. The hypothesis is supported by observations made during the present study, specifically that well-differentiated (SI) tumors were more methylated (~31%) compared with the poorly differentiated (SIII) tumors (~20%; Table III). Notably, analysis demonstrated that tumors from various histological origins exhibited varying degrees of methylation at the E2F4 promoter; this indicated that the methylation patterns could be a consequence of the activities of various methyl transferases and demethylases present at the tumor origin.

Generally, It is considered that gene silencing through aberrant promoter methylation is a consequence of aging (39). However, the present study identified an age-independent pattern to the E2F4 promoter methylation status. Δ meth values obtained for patients of >45 years and ≤45 of age, indicated that Δ meth values between the tumor and control groups are age-independent, and may have a causative role in tumor progression. However, these patients (≤45 and >45 years) exhibited significant Δ meth values (P=0.005) when compared with their respective controls (Table II). Age- and disease-dependent mutations in CpG-rich regions of E2F4 gene promoter can be a consequence of variable methylation levels (40) in these patients. Hence scanning these methylation hotspots (CpG islands) for single nucleotide polymorphisms may be an interesting avenue to explore. This will also highlight an interesting interplay of genetic and epigenetic factors in determining disease fate and other ageing phenotypes. The E2F4 gene has 10 reported C/T transitions at the targeted chromosomal location (31) flanking over ~205 bp around the E2F4 TSS. This supports the assumption that promoter methylation in this region may have a crucial role in promoting tumor formation. There is currently limited knowledge regarding the disease susceptibility of a specific population harboring these single-nucleotide polymorphisms. Based on these observations, analysis of the C>T transition in a population of patients with breast cancer and other associated diseases may be a notable domain to investigate.

Acknowledgements

All authors cordially acknowledge and are thankful to the participants directly and indirectly involved in the present study, including Dr Mumtaz Khan, Dr Irum Sabir Ali and Dr Ahmad Faraz (PGMI, Lady Reading Hospital, Peshawar, Pakistan). We also give thanks to the supporting staff of Cancer Genetics and Epigenetics Laboratory (Islamabad, Pakistan) for their support.

Funding

No funding was received.

Availability of data and materials

All data generated or analyzed during this study are included in this article.

Authors' contributions

FUF performed data acquisition and specimen analysis. FH performed data analysis and interpretation. NM performed data analysis and manuscript preparation. NA performed data analysis. HR performed data acquisition and interpretation. MS carried out study design, data analysis, interpretation, manuscript preparation and final approval of the version to be published.

Ethics approval and consent to participate

This study was initiated with a prior approval from departmental Ethical Review Committee (approval no. CIIT-09-10-14) of the COMSTAS Institute of Information Technology (Islamabad, Pakistan) and the Lady Reading Hospital (LRH) Peshawar. Informed and written consent was received from all patients enrolled in the present study, prior to sample collection for only data publication.

Consent for publication

All the information regarding patient's detail, are confidential and confined to Cancer Genetics and Epigenetics Laboratory, Department of Biosciences, COMSATS Institute of Information Technology (CIIT), 45550, Islamabad and the collaborating hospital LRH, Peshawar. Informed and written consent was received from all patients enrolled in the present study, prior to sample collection for only data publication.

Competing interests

The authors declare that they have no competing interests.

Glossary

Abbreviations

Abbreviations:

PCR

polymerase chain reaction

MSP

methylation specific PCR

Rb

retinoblastoma

TSS

transcription start site

qPCR

quantitative PCR

ANCT

adjacent normal control

Δ meth

change in methylation

bp

base pairs

References

1 

Zhang X, Ni Z, Duan Z, Xin Z, Wang H, Tan J, Wang G and Li F: Overexpression of E2F mRNAs associated with gastric cancer progression identified by the transcription factor and miRNA co-regulatory network analysis. PLoS One. 10:e01169792015. View Article : Google Scholar : PubMed/NCBI

2 

Chen HZ, Tsai SY and Leone G: Emerging roles of E2Fs in cancer: An exit from cell cycle control. Nat Rev Cancer. 9:785–797. 2009. View Article : Google Scholar : PubMed/NCBI

3 

Lee BK, Bhinge AA and Iyer VR: Wide-ranging functions of E2F4 in transcriptional activation and repression revealed by genome-wide analysis. Nucleic Acids Res. 39:3558–3573. 2011. View Article : Google Scholar : PubMed/NCBI

4 

Popov B, Chang LS and Serikov V: Cell cycle-related transformation of the E2F4-p130 repressor complex. Biochem Biophys Res Commun. 336:762–769. 2005. View Article : Google Scholar : PubMed/NCBI

5 

Crosby ME, Jacobberger J, Gupta D, Macklis RM and Almasan A: E2F4 regulates a stable G2 arrest response to genotoxic stress in prostate carcinoma. Oncogene. 26:1897–1909. 2007. View Article : Google Scholar : PubMed/NCBI

6 

Lee TA and Farnham PJ: Exogenous E2F expression is growth inhibitory before, during, and after cellular transformation. Oncogene. 19:2257–2268. 2000. View Article : Google Scholar : PubMed/NCBI

7 

Mott JL: MicroRNAs involved in tumor suppressor and oncogene pathways: Implications for hepatobiliary neoplasia. Hepatology. 50:630–637. 2009. View Article : Google Scholar : PubMed/NCBI

8 

Bertoli C, Skotheim JM and de Bruin RA: Control of cell cycle transcription during G1 and S phases. Nat Rev Mol Cell Biol. 14:518–528. 2013. View Article : Google Scholar : PubMed/NCBI

9 

Crijns AP, Fehrmann RS, de Jong S, Gerbens F, Meersma GJ, Klip HG, Hollema H, Hofstra RM, te Meerman GJ, de Vries EG and van der Zee AG: Survival-related profile, pathways, and transcription factors in ovarian cancer. PLoS Med. 6:e10000242009. View Article : Google Scholar

10 

Chong JL, Wenzel PL, Saenz-Robles MT, Nair V, Ferrey A, Hagan JP, Gomez YM, Sharma N, Chen HZ, Ouseph M, et al: E2f1-3 switch from activators in progenitor cells to repressors in differentiating cells. Nature. 462:930–934. 2009. View Article : Google Scholar : PubMed/NCBI

11 

Lee EY, Cam H, Ziebold U, Rayman JB, Lees JA and Dynlacht BD: E2F4 loss suppresses tumorigenesis in Rb mutant mice. Cancer Cell. 2:463–472. 2002. View Article : Google Scholar : PubMed/NCBI

12 

Schwemmle S and Pfeifer GP: Genomic structure and mutation screening of the E2F4 gene in human tumors. Int J Cancer. 86:672–677. 2000. View Article : Google Scholar : PubMed/NCBI

13 

Souza RF, Yin J, Smolinski KN, Zou TT, Wang S, Shi YQ, Rhyu MG, Cottrell J, Abraham JM, Biden K, et al: Frequent mutation of the E2F-4 cell cycle gene in primary human gastrointestinal tumors. Cancer Res. 57:2350–2353. 1997.PubMed/NCBI

14 

Khaleel SS, Andrews EH, Ung M, DiRenzo J and Cheng C: E2F4 regulatory program predicts patient survival prognosis in breast cancer. Breast Cancer Res. 16:4862014. View Article : Google Scholar : PubMed/NCBI

15 

Fuchs B, Zhang K, Schabel A, Bolander ME and Sarkar G: Identification of twenty-two candidate markers for human osteogenic sarcoma. Gene. 278:245–252. 2001. View Article : Google Scholar : PubMed/NCBI

16 

Kim EY, Jo JH and Lee MJ: Expressions of E2F4 and E2F2 transcription factors in breast carcinoma. Korean J Pathol. 39:301–306. 2005.

17 

Enos M, Bancos SA, Bushnell I and Crispe IN: E2F4 modulates differentiation and gene expression in hemopoietic progenitor cells during commitment to the lymphoid lineage. J Immunol. 180:3699–3707. 2008. View Article : Google Scholar : PubMed/NCBI

18 

Ma Y, Freeman SN and Cress WD: E2F4 deficiency promotes drug-induced apoptosis. Cancer Biol Ther. 3:1262–1269. 2004. View Article : Google Scholar : PubMed/NCBI

19 

Lassmann S, Weis R, Makowiec F, Roth J, Danciu M, Hopt U and Werner M: Array CGH identifies distinct DNA copy number profiles of oncogenes and tumor suppressor genes in chromosomal- and microsatellite-unstable sporadic colorectal carcinomas. J Mol Med (Berl). 85:293–304. 2007. View Article : Google Scholar : PubMed/NCBI

20 

Sharma S, Kelly TK and Jones PA: Epigenetics in cancer. Carcinogenesis. 31:27–36. 2010. View Article : Google Scholar : PubMed/NCBI

21 

Baxter E, Windloch K, Gannon F and Lee JS: Epigenetic regulation in cancer progression. Cell Biosci. 4:452014. View Article : Google Scholar : PubMed/NCBI

22 

Ng JM and Yu J: Promoter hypermethylation of tumour suppressor genes as potential biomarkers in colorectal cancer. Int J Mol Sci. 16:2472–2496. 2015. View Article : Google Scholar : PubMed/NCBI

23 

Fearon ER: Molecular genetics of colorectal cancer. Annu Rev Pathol. 6:479–507. 2011. View Article : Google Scholar : PubMed/NCBI

24 

Sabir M, Baig RM, Ali K, Mahjabeen I, Saeed M and Kayani MA: Retinoblastoma (RB1) pocket domain mutations and promoter hyper-methylation in head and neck cancer. Cell Oncol (Dordr). 37:203–213. 2014. View Article : Google Scholar : PubMed/NCBI

25 

Ullah F, Khan T, Ali N, Malik FA, Kayani MA, Shah Abbas ST and Saeed M: Promoter methylation status modulate the expression of tumor suppressor (RbL2/p130) gene in breast cancer. PLoS One. 10:e01346872015. View Article : Google Scholar : PubMed/NCBI

26 

Ali A, Ullah F, Sabir AI, Faraz A, Khan M, Shah Ali ST, Ali N and Saeed M: Aberrant promoter methylation at CpG cytosines induce the upregulation of the E2F5 gene in breast cancer. J Breast Cancer. 19:133–141. 2016. View Article : Google Scholar : PubMed/NCBI

27 

Lee JT, Tsan WH and Chow KL: Simple modifications to standard TRIzol® protocol allow high-yield RNA extraction from cells on resorbable materials. J Biomater Nanobiotechnol. 2:41–48. 2011. View Article : Google Scholar

28 

Wang Y, Zhu W and Levy DE: Nuclear and cytoplasmic mRNA quantification by SYBR green based real-time RT-PCR. Methods. 39:356–362. 2006. View Article : Google Scholar : PubMed/NCBI

29 

Livak KJ and Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 25:402–408. 2001. View Article : Google Scholar : PubMed/NCBI

30 

Strauss WM: Preparation of genomic DNA from mammalian tissue. Curr Protoc Mol Biol: Chapter 2. 2001. View Article : Google Scholar

31 

Li LC and Dahiya R: MethPrimer: Designing primers for methylation PCRs. Bioinformatics. 18:1427–1431. 2002. View Article : Google Scholar : PubMed/NCBI

32 

Wang D, Russell JL and Johnson DG: E2F4 and E2F1 have similar proliferative properties but different apoptotic and oncogenic properties in vivo. Mol Cell Biol. 20:3417–3424. 2000. View Article : Google Scholar : PubMed/NCBI

33 

Bandara LR, Buck VM, Zamanian M, Johnston LH and La Thangue NB: Functional synergy between DP-1 and E2F-1 in the cell cycle-regulating transcription factor DRTF1/E2F. EMBO J. 12:4317–4324. 1993.PubMed/NCBI

34 

Nevins JR: E2F: A link between the Rb tumor suppressor protein and viral oncoproteins. Science. 258:424–429. 1992. View Article : Google Scholar : PubMed/NCBI

35 

Beijersbergen RL, Kerkhoven RM, Zhu L, Carlée L, Voorhoeve PM and Bernards R: E2F-4, a new member of the E2F gene family, has oncogenic activity and associates with p107 in vivo. Genes Dev. 8:2680–2690. 1994. View Article : Google Scholar : PubMed/NCBI

36 

Ginsberg D, Vairo G, Chittenden T, Xiao ZX, Xu G, Wydner KL, DeCaprio JA, Lawrence JB and Livingston DM: E2F-4, a new member of the E2F transcription factor family, interacts with p107. Genes Dev. 8:2665–2679. 1994. View Article : Google Scholar : PubMed/NCBI

37 

Le Cam L, Polanowska J, Fabbrizio E, Olivier M, Philips A, Ng Eaton E, Classon M, Geng Y and Sardet C: Timing of cyclin E gene expression depends on the regulated association of a bipartite repressor element with a novel E2F complex. EMBO J. 18:1878–1890. 1999. View Article : Google Scholar : PubMed/NCBI

38 

Zhu WG, Srinivasan K, Dai Z, Duan W, Druhan LJ, Ding H, Yee L, Villalona-Calero MA, Plass C and Otterson GA: Methylation of adjacent CpG sites affects Sp1/Sp3 binding and activity in the p21(Cip1) promoter. Mol Cell Biol. 23:4056–4065. 2003. View Article : Google Scholar : PubMed/NCBI

39 

Almen MS, Nilsson EK, Jacobsson JA, Kalnina I, Klovins J, Fredriksson R and Schiöth HB: Genome-wide analysis reveals DNA methylation markers that vary with both age and obesity. Gene. 548:61–67. 2014. View Article : Google Scholar : PubMed/NCBI

40 

Jin B and Robertson KD: DNA methyltransferases, DNA damage repair, and cancer. Adv Exp Med Biol. 754:3–29. 2013. View Article : Google Scholar : PubMed/NCBI

Related Articles

Journal Cover

June-2018
Volume 15 Issue 6

Print ISSN: 1792-1074
Online ISSN:1792-1082

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Farman FU, Haq F, Muhammad N, Ali N, Rahman H and Saeed M: Aberrant promoter methylation status is associated with upregulation of the E2F4 gene in breast cancer. Oncol Lett 15: 8461-8469, 2018
APA
Farman, F.U., Haq, F., Muhammad, N., Ali, N., Rahman, H., & Saeed, M. (2018). Aberrant promoter methylation status is associated with upregulation of the E2F4 gene in breast cancer. Oncology Letters, 15, 8461-8469. https://doi.org/10.3892/ol.2018.8382
MLA
Farman, F. U., Haq, F., Muhammad, N., Ali, N., Rahman, H., Saeed, M."Aberrant promoter methylation status is associated with upregulation of the E2F4 gene in breast cancer". Oncology Letters 15.6 (2018): 8461-8469.
Chicago
Farman, F. U., Haq, F., Muhammad, N., Ali, N., Rahman, H., Saeed, M."Aberrant promoter methylation status is associated with upregulation of the E2F4 gene in breast cancer". Oncology Letters 15, no. 6 (2018): 8461-8469. https://doi.org/10.3892/ol.2018.8382