Disease Markers
Volume 2015 (2015), Article ID 241301, 13 pages
http://dx.doi.org/10.1155/2015/241301

Research Article

Identification of Novel Epigenetic Markers of Prostate Cancer by NotI-Microarray Analysis

Alexey A. Dmitriev,1,2 Eugenia E. Rosenberg,³ George S. Krasnov,¹ Ganna V. Gerashchenko,³ Vasily V. Gordiyuk,³ Tatiana V. Pavlova,⁴ Anna V. Kudryavtseva,¹ Artemy D. Beniaminov,¹ Anastasia A. Belova,¹ Yuriy N. Bondarenko,⁵ Rostislav O. Danilets,⁵ Alexander I. Glukhov,⁶ Aleksandr G. Kondratov,³ Andrey Alexeyenko,⁷ Boris Y. Alekseev,² George Klein,⁴ Vera N. Senchenko,¹ and Vladimir I. Kashuba3,4

¹Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, Moscow 119991, Russia
²P.A. Herzen Moscow Cancer Research Institute, Ministry of Healthcare of the Russian Federation, Moscow 125284, Russia
³Institute of Molecular Biology and Genetics, National Academy of Sciences of Ukraine, Kiev 03680, Ukraine
⁴Department of Microbiology, Tumor and Cell Biology, Karolinska Institute, 17177 Stockholm, Sweden
⁵Institute of Urology, National Academy of Medical Sciences of Ukraine, Kiev 04053, Ukraine
⁶Department of Molecular Biology, Kurchatov NBIC Centre NRC “Kurchatov Institute”, Moscow 123182, Russia
⁷Bioinformatics Infrastructure for Life Sciences, Science for Life Laboratory, Karolinska Institute, 17177 Stockholm, Sweden

Received 9 February 2015; Revised 11 July 2015; Accepted 14 July 2015

Academic Editor: Stamatios Theocharis

Copyright © 2015 Alexey A. Dmitriev et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A significant need for reliable and accurate cancer diagnostics and prognosis compels the search for novel biomarkers that would be able to discriminate between indolent and aggressive tumors at the early stages of disease. The aim of this work was identification of potential diagnostic biomarkers for characterization of different types of prostate tumors. NotI-microarrays with 180 clones associated with chromosome 3 genes/loci were applied to determine genetic and epigenetic alterations in 33 prostate tumors. For 88 clones, aberrations were detected in more than 10% of tumors. The major types of alterations were DNA methylation and/or deletions. Frequent methylation of the discovered loci was confirmed by bisulfite sequencing on selective sampling of genes: FGF12, GATA2, and LMCD1. Three genes (BHLHE40, BCL6, and ITGA9) were tested for expression level alterations using qPCR, and downregulation associated with hypermethylation was shown in the majority of tumors. Based on these data, we proposed the set of potential biomarkers for detection of prostate cancer and discrimination between prostate tumors with different malignancy and aggressiveness: BHLHE40, FOXP1, LOC285205, ITGA9, CTDSPL, FGF12, LOC440944/SETD5, VHL, CLCN2, OSBPL10/ZNF860, LMCD1, FAM19A4, CAND2, MAP4, KY, and LRRC58. Moreover, we probabilistically estimated putative functional relations between the genes within each set using the network enrichment analysis.

1. Introduction

Prostate cancer is the second most common cancer type and the fifth leading cause of death from cancer among men worldwide [1]. According to the International Agency for Research on Cancer data morbidity and mortality rates of prostate cancer in 2012 were 1.1 million (15%) and 307 thousands (6.6%), respectively. Prostate cancer is heterogeneous disease divided into several stages. Hyperactivation of androgen receptor pathway (AR) plays a central role in the tumor development. Initially most of prostate malignant tumors are dependent on external androgens in blood. Androgen deprivation approaches are effective for suppression of prostate cancer progression at this stage [2]. However, several mechanisms, synthesis of endogenous androgen, AR mutations and induction of AR splice isoforms, and others, makes tumor independent of presence of external androgens, and this leads to treatment failure [3]. Moreover, some of tumors proceed neuroendocrine differentiation, which makes their growth completely independent of AR pathway [4]. The choice of therapy for an individual is dependent on a number of factors but it is well recognized that different therapies may work equally well. Conversely, many patients will fail a particular treatment despite apparently favorable disease characteristics [5]. The effectiveness of radical prostatectomy for locally advanced prostate cancer is also controversial and is a focus of debate [6, 7]. The identification of markers of tumor aggressiveness is highly demanded [8–10].

The first protein biomarkers proposed for the diagnosis of prostate cancer were serum prostatic acid phosphatase encoded by ACPP gene [11] and prostate-specific antigen encoded by KLK3 gene [12]. However, due to insufficient specificity and sensitivity of prostatic acid phosphatase [13], prostate-specific antigen (PSA) was the gold standard for prostate cancer diagnostics in clinic for a long time. Recent advancement in molecular biology has led to reporting of many novel biomarkers. Moreover, the PSA-based screening nowadays are associated with more harms than benefits to patients due to excessive false positive rate, which may result in subsequent wrong diagnosis and overtreatment [14]. In 2012, the US Food and Drug Administration approved prostate cancer associated 3 (PSA3) as a novel marker of prostate cancer to use when the combined results of repeated biopsy, PSA concentration, and digital rectal examination are controversial. PCA3 together with the prostate health index (phi) showed high prognostic ability [15]. However, approved clinical biomarkers for discrimination between aggressive and nonaggressive types of prostate cancer are still lacking [16].

It is known that chromosomes undergo different types of genetic and epigenetic changes in carcinogenesis. DNA methylation affects CpG-rich regions referred to as CpG-islands [17]. DNA hypermethylation in cancer cells is associated with gene silencing. Most recurrently, the silencing is observed in tumor suppressor genes [18]. On the other hand, the DNA hypomethylation is associated with activation of gene expression, which can affect normally silenced oncogenes. Therefore, this common mechanism of epigenetic aberration could be helpful for searching potential markers of prostate cancer progression [19, 20].

Prostate cancer DNA methylation profiling demonstrated high prognostic value [21–23]. The use of a DNA methylation-based biomarker for prostate cancer is appealing for several reasons: the high stability of DNA, ease of analysis with the current techniques available, and the ability to assess the biomarker in body fluids such as blood, urine, and saliva use [24]. Methylation of several genes was shown to significantly correlate with disease recurrence. Among them are well-known tumors suppressor RASSF1 [22]; cell surface glycoprotein CD44 participating cells interactions, adhesion, and migration; cyclooxygenase PTGS2 [23]; transcription factor PITX2 [25]. However, none of these markers has reached clinical use [24]. The most known potential prostate cancer methylation biomarker is glutathione-S-transferase P1 (GSTP1) gene, which encodes an enzyme required for detoxification and protection of DNA from oxidants and electrophilic metabolites. GSTP1 gene methylation status was shown to be able to discriminate between prostate cancer and benign hyperplasia and predict disease recurrence [26, 27]. Interestingly, GSTP1 methylation in adjacent nonneoplastic tissues also correlated with clinical outcome [28].

Multigene approaches perform single-gene prognosis and diagnosis methods [29, 30]. The data obtained in our laboratory [31–33] and other groups [34] suggest that aberrations in the human chromosome 3 frequently accompany the formation of tumors of the epithelial origin. A number of regions in the short arm of chromosome 3 are often either deleted or methylated in cancer genomes. The relevance of these alternative types of activity can be confirmed by a recent meta-analysis of cancer genomics data that clearly demonstrated that the same cancer driver genes might experience alterations via different molecular mechanisms, including methylation and copy number changes [35]. These observations suggest their roles as tumor suppressor genes [36, 37] and warrant evaluation as potential cancer markers.

The aim of this research was to identify a set of novel potential markers for high accuracy early detection of prostate tumors and discrimination between adenoma and carcinoma, aggressive and nonaggressive cases.

We screened genetic (deletions) and epigenetic (DNA methylation) changes in prostate biopsy samples with different pathologies using the NotI-microarrays (NMA) earlier developed by Dr. Eugene R. Zabarovsky [38]. The NMA results were verified for several genes by bisulfite genomic sequencing to confirm that DNA methylation is a frequent inactivation mechanism in prostate tumors. Moreover, expression downregulation was shown for three hypermethylated genes in the majority of examined samples using quantitative PCR (qPCR). Finally, we developed novel sets of potential biomarkers for detection and discrimination among prostate tumors with different pathomorphological characteristics.

2. Materials and Methods

2.1. Sample Collection

Prostate biopsy samples with different pathologies were collected from patients of the Institute of Urology (Kiev, Ukraine). In total, there were 33 samples, among which 15 samples of prostate adenoma, 14 of prostate carcinoma with Gleason score ≤7, and 4 of prostate aggressive carcinoma with Gleason score >7. Gleason score is a sum of grades from two most representative biopsy specimens, which are microscopically examined by pathologist. A pool from four normal prostate inflammation samples was used as a reference. Three normal prostate inflammation samples obtained later were used to verify the suitability of the pool as a reference (no valid aberrations were observed according to NMA). All patients gave written informed consent. The samples were collected in accordance with the Declaration of Helsinki and the guidelines issued by the Ethic Committee of the Institute of Urology, National Academy of Medical Sciences of Ukraine. The Ethics Committee of the Institute of Urology (Kiev, Ukraine) specifically approved this study.

2.2. DNA and RNA Extraction and cDNA Synthesis

The DNA was isolated using phenol-chloroform extraction according to the Maniatis protocol [39]. Total RNA extraction was done with RNeasy Mini kit (Qiagen, Germany) according to the manufacturer’s recommendations. RNA quality was monitored with absorbance spectra (NanoDrop Technologies Inc., USA) and the RNA integrity number (RIN; Agilent Technologies, USA). cDNA was synthesized using Maxima Reverse Transcriptase (Thermo Fisher Scientific, USA) and random primers.

2.3. NotI-Microarray

For the NotI-microarray we used 180 clones of the human chromosome 3. The NotI-probes were prepared as described previously [40–42].

Briefly, the hybridization of coupled NotI-enriched samples was done at 42°C for 15 h in a Lucidea Base device (Amersham Pharmacia Biotech) according to the manufacturer’s instructions. The microarrays were scanned in GenePix 4000A. The results were processed with GenePix Pro 6.0 software (Amersham Pharmacia Biotech). Then the data were analyzed using our program NIMAN [43].

2.4. Bisulfite Genomic Sequencing

The bisulfite conversion of DNA was performed using EZ DNA Methylation Kit (Zymo Research) according to the manufacturer’s instructions. Primers for PCR for FGF12 were as follows:

 forward 5′-ACATTTTCTCCTTAGGACCAAGGGAAG-3′; reverse 5′-CTGCAGCCTCCTCAAATTTTAGCACTGC-3′.

After amplification of the bisulfite treated DNA, the PCR products were cloned and used for automated sequencing or were sequenced directly (ABI Prism 3100-Avant Genetic Analyzer, Applied Biosystems).

Quantitative PCR. Expression of BHLHE40, BCL6, and ITGA9 genes was evaluated using commercial sets of primers and probes (Applied Biosystems, USA) and 7500 Real-Time PCR System (Applied Biosystems). Each reaction was repeated three times. QPCR data were analyzed using three reference genes: GAPDH, ACTB, and B2M. Relative quantification of gene expression was performed as described earlier [44, 45]. At least 2-fold mRNA level alterations were considered significant.

2.5. Statistical Analysis

Fisher’s exact test and criteria were used for analysis of methylation and/or deletion frequencies in groups of prostate tumors with different pathomorphological characteristics. Cases with value below 0.05 were considered statistically significant. Sets of markers for identification and discrimination of prostate tumors were developed using the support vector machine [46]. Sensitivity, specificity, and accuracy of the sets were calculated as the proportion of true positive results, true negative results, and true positive plus true negative results, respectively. Gini coefficient was used to evaluate the predictive power of the developed sets [47]. All statistical procedures were performed using our NIMAN software [43].

2.6. Network Enrichment Analysis

We probabilistically estimated putative functional relations of the obtained gene sets using the methodology of network enrichment analysis (NEA) as described in [48]. The biological network connectivity between genes of the novel lists and genes of known KEGG pathways was quantified as total numbers of links (edges) found in the global interaction network that connected any genes of the novel list to any genes of a given KEGG pathway. Similarly, we quantified the enrichment within the novel lists by counting any links between any gene pairs of each list. In this latter analysis we also utilized the indirect link mode by counting shared network neighbors in such genes pairs. In this analysis we utilized the global network of functional coupling FunCoup [49] which included all edges with confidence score higher than 0.5 with addition of all known links from the curated databases KEGG [50], PhosphoSite, CORUM [51], MSigDB [52], and HTRIdb [53]. This procedure gave a union network of 974,427 functional links between 19,031 distinct HUPO gene symbols.

3. Results and Discussion

The genomic DNA from 15 adenoma and 18 carcinoma prostate biopsy samples was analyzed by NotI-microarrays. A pool of samples from 4 patients with nonmalignant prostate inflammation was used as a reference sample.

The restriction site NotI (GCGGCCGC) is frequently located in CpG-islands within promoter regions and is thus the sensitivity to NotI digestion that reflects the methylation status of the surrounding genes/loci [54]. In the present work, we studied 180 genomic loci (clones) associated with genes from chromosome 3. We detected changes in more than 10% of tumors for 88 out of 180 clones. The most frequent type of changes in all prostate samples (33) was heterozygous deletion/methylation (1163 cases), followed by homozygous deletion/methylation (461 cases), and only 24 cases of amplification/demethylation. The 50 genes/loci with the highest percent of changes in prostate adenoma and cancer are shown in Figure 1.

Figure 1: Pattern of DNA alterations in prostate tumor samples. NotI-microarray data. Horizontally: 33 prostate samples (15 adenomas, 14 nonaggressive carcinomas, and 4 aggressive carcinomas). : Gleason grading system score. Vertically: gene-linked 50 NotI-sites arranged by methylation/deletion frequency (from 82% to 33%). Green squares: methylation/deletion, red: amplification/demethylation, yellow: unchanged, and white: no info.

According to PubMed database search all genes/loci with the highest percent of changes could be divided into three groups: (1) previously not shown to be associated with carcinogenesis: HMGB1L5, LRRC58, GPR149, DZIP1L, C3orf77, and NUDT16; (2) known to be involved in nonprostate cancer: LOC285205, KY, BHLHE40, ROPN1/KALRN, BCL6, PLCL2, ITGA9, CTDSPL (RBSP3), GORASP1/TTC21A, FSTL1, ABHD5/C3orf77, IQSEC1, CLASP2, GNAI2, NEK11, FBLN2, SOX2, MINA, CHCHD6, WNT7A, LOC285375, FGF12, NKIRAS1/RPL15, CGGBP1, PPP2R3A, SOX14, ZIC4, RAP2B, RPL32/IQSEC1, RRP9/PARP3, PPM1M, KBTBD8, FGD5, CMTM8, NBEAL2, TMEM45A, LRRC3B, PDZRN3, USP19, and EPHB1; (3) previously shown to be associated with prostate cancer: FOXP1, MANF, GATA2, ALDH1L1, and EPHB3.

The 50 NotI-sites (gene/locus) with the highest percent of DNA methylation/deletion in prostate adenoma and cancer are shown in Table 1.

Table 1: Methylation/deletion frequencies for 50 NotI-sites (gene/locus) with the highest percent of changes in prostate adenoma and cancer.

No
NotI-site
Gene/locus
Localization
Methylation/deletion frequency

Adenoma
Cancer ()
Cancer ()

1
NR5-IH18RS
BHLHE40
3p26.1
80% (14/15)
79% (11/14)
100% (4/4)

2
NR1-XM13C
IQSEC1
3p25.2
73% (11/15)
79% (11/14)
100% (4/4)

3
NR1-AK24R
BCL6
3q27
73% (11/15)
79% (11/14)
100% (4/4)

4
NL1-BA6R
FOXP1
3p14.1
67% (10/15)
79% (11/14)
100% (4/4)

5
NR1-KJ5R (C)
FBLN2
3p25.1
73% (11/15)
64% (9/14)
100% (4/4)

6
NR1-EP7RS
CLASP2
3p22.3
80% (12/15)
57% (8/14)
100% (4/4)

7
NL1Z216R (D)
MANF
3p21.1
80% (12/15)
64% (9/14)
75% (3/4)

8
NL1-GK21R (C)
ROPN1/KALRN
3q13.3
60% (9/15)
79% (11/14)
100% (4/4)

9
NL4-AP18R (C)
PLCL2
3p24.3
53% (8/15)
86% (12/14)
75% (3/4)

10
NL1A401R (D)
ITGA9
3p21.3
60% (9/15)
71% (10/14)
100% (4/4)

11
NLJ-003RD
CTDSPL (RBSP3)
3p21.3
73% (11/15)
64% (9/14)
75% (3/4)

12
NL3003R (U)
GORASP1/TTC21A
3p22-p21.33
67% (10/15)
71% (10/14)
75% (3/4)

13
NL3-CI2R (C)
LOC285205
3q13.12
73% (11/15)
57% (8/14)
100% (4/4)

14
NR5-FG18R (C)
FSTL1
3q13.33
73% (11/15)
57% (8/14)
100% (4/4)

15
NR1-WD21R (C)
NEK11/NUDT16
3q22.1
60% (9/15)
64% (9/14)
100% (4/4)

16
NL3A001R (D)
GNAI2
3p21.31
53% (8/15)
64% (9/14)
100% (4/4)

17
NR1-AN24RS
ABHD5/C3orf77
3p21
53% (8/15)
71% (10/14)
50% (2/4)

18
NR1-PD1R
ZIC4
3q24
47% (7/15)
64% (9/14)
100% (4/4)

19
NL1-VJ14R (C)
LOC285375
3p25.1
60% (9/15)
43% (6/14)
100% (4/4)

20
NL4-BC8R (C)
ALDH1L1
3q21.3
33% (5/15)
79% (11/14)
75% (3/4)

21
NL1-YJ5R (C)
С3orf46/CHCHD6
3q21.3
60% (9/15)
50% (7/14)
75% (3/4)

22
NL1-GC10C
HMGB1L5 (Pseudo)
3p24
67% (10/15)
43% (6/14)
50% (2/4)

23
NR5-FK16RS
MINA
3q11.2
40% (6/15)
57% (8/14)
100% (4/4)

24
NL1-FK10R (C)
PPP2R3A
3q22.1
33% (5/15)
64% (9/14)
100% (4/4)

25
NL1-ZD4R
SOX2
3q26.3-q27
40% (6/15)
57% (8/14)
100% (4/4)

26
NR1-NH1R (C)
FGF12
3q28
53% (8/15)
64% (9/14)
25% (1/4)

27
NR5-FK11R (C)
CMTM8
3p22.3
67% (10/15)
50% (7/14)
0% (0/4)

28
NL4-BH3R (C)
GATA2
3q21.3
20% (3/15)
71% (10/14)
100% (4/4)

29
NL4-BK12R (C)
WNT7A
3p25
40% (6/15)
43% (6/14)
100% (4/4)

30
NR1-WD23R (C)
LRRC58
3q13.33
53% (8/15)
29% (4/14)
100% (4/4)

31
NL4-BI4RS
RAP2B
3q25.2
47% (7/15)
57% (8/14)
25% (1/4)

32
NL1-ZP13R (C)
KBTBD8
3p14
47% (7/15)
50% (7/14)
25% (1/4)

33
NR1-WE11RS
CGGBP1
3p12-p11.1
33% (5/15)
43% (6/14)
100% (4/4)

34
HSJ4-AB7R (C)
RPL32/IQSEC1
3p25.2
27% (4/15)
43% (6/14)
100% (4/4)

35
NL4-DP2RS
FGD5
3p25.1
33% (5/15)
36% (5/14)
100% (4/4)

36
NL1-CJ4R (C)
NKIRAS1/RPL15
3p24.2
40% (6/15)
50% (7/14)
25% (1/4)

37
NL6-II3R
USP19
3p21.31
27% (4/15)
43% (6/14)
100% (4/4)

38
NL1268R (P65D)
TMEM45A
3q12.2
27% (4/15)
50% (7/14)
75% (3/4)

39
NR1-WJ2RS
SOX14
3q22-q23
27% (4/15)
57% (8/14)
50% (2/4)

40
NL1-VC9R (C)
GPR149
3q25.2
47% (7/15)
43% (6/14)
25% (1/4)

41
NR1-WH9R (C)
RRP9/PARP3
3p21.2
47% (7/15)
29% (4/14)
50% (2/4)

42
NR1-NC7RS
PPM1M
3p21.2
40% (6/15)
43% (6/14)
25% (1/4)

43
NR5-IG2R (C)
KY
3q22.2
60% (9/15)
14% (2/14)
50% (2/4)

44
NR1-WL7R (C)
DZIP1L
3q22.3
53% (8/15)
36% (5/14)
0% (0/4)

45
NR1-NM7R (C)
B3GALNT1
3q25
27% (4/15)
43% (6/14)
75% (3/4)

46
NL3-CA11RS
LRRC3B
3p24
27% (4/15)
29% (4/14)
100% (4/4)

47
NL6-AF21R (C)
PDZRN3
3p13
40% (6/15)
29% (4/14)
50% (2/4)

48
NL1A079R (D)
EPHB1
3q21-q23
27% (4/15)
36% (5/14)
75% (3/4)

49
NL3A006R (D)
NBEAL2
3p21.31
20% (3/15)
36% (5/14)
75% (3/4)

50
NR1-WB21R (C)
EPHB3
3q21-qter
40% (6/15)
36% (5/14)
0% (0/4)

Note. : Gleason grading system score.

According to our recent results, some members of the second and third groups (CTDSPL, ALDH1L1, LRRC3B, IQSEC1, FOXP1, GNAI2, EPHB1, WNT7A, etc.) manifested frequent methylation/deletion in tumors of other localizations, such as colorectal cancer [55], nonsmall cell lung cancer [43], cervical cancer [56], clear cell renal cell carcinoma [57, 58], and high-grade serous ovarian cancer [59]. It may indicate potential tumor suppressor role of these genes/loci and needs further investigation for prostate cancer.

In order to validate the results of NotI-microarrays we performed bisulfite sequencing for 4 genes: FGF12, GATA2, LMCD1, and TESSP2 (control). The fibroblast growth factor (FGF) family exhibits a broad range of mitogenic and cell survival activities. FGF12 was shown to be significantly methylated in breast and colorectal cancer [60, 61]. On the other hand, it was found to be overexpressed in lung squamous cell carcinoma [62]. Thus, its possible tumor suppressor role needs further investigation. We selected a few samples (, = 4–9) of prostate cancer with genetic/epigenetic changes (according to NMA data) for bisulfite sequencing. The results of bisulfite sequencing for FGF12 gene promoter region in the sample with Gleason Score 9 are shown in Figure 2. Amplified and cloned region contained 34 CG-pairs including 2 of NotI-site. NotI-site methylation was observed in 6 out of 11 clones (55%) and associated with dense methylation of the entire sequenced region. The other 11 prostate cancer samples showed 40–80% rate of NotI-site methylation. High density of methylation in the examined samples (30–70%) was also observed for GATA2 and LMCD1 genes and almost no methylation (<10%) was detected for TESSP2, which was used as a negative control. Thus, the results of the bisulfite sequencing were in concordance with the NotI-microarray data and confirmed that methylation is a frequent event in prostate carcinomas.

(a)

(b)

Figure 2: Bisulfite sequencing data for the promoter region of FGF12 gene in prostate carcinoma with Gleason score = 9. Thirty-four CG-pairs (a) are shown in bold and grey. Primers for bisulfite sequencing (a) are in bold and italic. NotI-site is underlined (a). In the table (b) methylated (grey squares) and unmethylated (white squares) CG-pairs are shown in eleven sequenced clones. CG-pairs that correspond to NotI-site ((7) and (8)) are in bold and underlined. Crossed squares: no data.

Expression level of three genes with high methylation/deletion frequency (BHLHE40, BCL6, and ITGA9) was evaluated in 11 prostate cancer samples using qPCR. All genes showed on the average 3-fold downregulation in the majority of cases (Figure 3). In the rest cases, the retention of mRNA level was observed. In 5 of 11 samples, expression decrease was observed for all three genes simultaneously. Thus, high frequency of methylation/deletion was associated with expression downregulation for BHLHE40, BCL6, and ITGA9 genes.

Figure 3: Relative mRNA level of BHLHE40, BCL6, and ITGA9 genes in 11 prostate cancer samples. QPCR data. Grey dashed lines represent 2-fold alteration interval.

Using our NotI-microarray data, we attempted to construct prediction models for detection of prostate cancer and discrimination between aggressive and nonaggressive subtypes. To detect prostate tumors, a set of six markers could be proposed: BHLHE40, FOXP1, LOC285205, ITGA9, CTDSPL (RBSP3), and FGF12. For a sample to be classified as a prostate tumor, we required DNA methylation and/or deletion to be detected in two or more of these marker genes. Under this condition, both the sensitivity and the specificity of this set were 94% for the examined sampling. The Gini coefficient was in the range 0.92–0.99 (Figure 4(a)). However, these values need further refinement using additional samplings.

(a)

(b)

(c)

Figure 4: Lift diagrams for 3 prediction models: detection of prostate tumors (a), isolation of aggressive prostate cancer (b), and discrimination between nonaggressive prostate cancer and adenoma (c). Solid black line: ideal lift curve; solid grey line: lift curve under the most favorable conditions; dashed black line: lift curve under the least favorable conditions; the diagonal line corresponds to a random guess.

Gene BHLHE40 (region 3p26) encodes a transcription factor of the basic helix-loop-helix family. It is a candidate tumor suppressor gene [63]. Gene ITGA9 (3p21.3) encodes alpha 9 integrin. Gene RBSP3 (3p21.3) encodes CTD small phosphatase-like protein involved in activation of retinoblastoma protein (Rb). Until recently, little was known about functions of locus LOC285205. According to the latest HGNC annotation, it encodes LINC00636, the long intergenic nonprotein coding RNA 636 (HGNC Acc. 27702). We noted that its surrounding regions contain a potential promoter site to at least four transcription factors (CUX1, GATA1, POU2F1, and POU3F2), each of which also has binding sites in the promoter region of BHLHE40. Although these promoters are located 35–45 thousands base pairs upstream of LINC00636, there are no alternative gene coding regions closer than 100 thousands base pairs, except another noncoding RNA LINC00635. The deletion in 3p14 region containing FOXP1 (forkhead box P1 protein) gene is associated with TMPRSS2-ERG fusion events which are very common in prostate cancer genomes. This deletion might suggest the tumor suppressor role of FOXP1 [64].

By comparing aggressive prostate carcinoma cases (Gleason score > 7) against adenomas and nonaggressive (Gleason score ≤ 7) carcinomas, we identified 9 genes/loci with statistically significant association (Table 2).

Table 2: Methylation/deletion frequencies in aggressive prostate cancer versus adenoma plus nonaggressive cancer.

Gene/locus
Methylation/deletion frequency, %
value

Aggressive cancer
Adenoma and nonaggressive cancer

LOC440944/SETD5
100 (4/4)
10 (3/29)
<0.001

OSBPL10/ZNF860
100 (4/4)
7 (2/29)
<0.001

CLCN2
100 (4/4)
7 (2/29)
<0.001

PRSS42/MYL3
100 (4/4)
0 (0/29)
<0.001

VHL
75 (3/4)
0 (0/29)
<0.001

BBX
100 (4/4)
17 (5/29)
0.003

LMCD1
100 (4/4)
21 (6/29)
0.005

CMTM6
100 (4/4)
21 (6/29)
0.005

FAM19A4
100 (4/4)
21 (6/29)
0.005

The promising set of 5 markers enabling isolation of aggressive cases of prostate cancer could be proposed: LOC440944/SETD5, VHL, CLCN2, OSBPL10/ZNF860, and LMCD1. A sample would be recognized as aggressive given methylation and/or deletion detected in three or more of these markers (100% sensitivity and 97% specificity). The Gini coefficient was in the range 0.93–1.00 (Figure 4(b)).

Gene VHL (von Hippel-Lindau tumor suppressor gene; 3p25.3) is frequently lost in clear-cell renal carcinomas [65]. Increased levels of VHL induced apoptosis in prostate cells [66]. CLCN2 (3q27-q28) encodes the chloride channel 2 protein. Little is known about its role in cancerogenesis, although this gene was suggested as a novel drug target for tumor inhibition in malignant glioma cells [67]. Gene LMCD1 (3p26-p24) encodes LIM and cysteine-rich domains protein 1 (Dyxin). This gene is suggested as a potential oncogene in hepatocellular carcinoma [68]; however, its role in prostate cancer is not known. Locus LOC440944/SETD5 is situated in 3p25.3 region and probably is a noncoding RNA (NCBI Gene ID: 440944) with unknown functions. Gene OSBPL10 (3p22.3) encodes oxysterol binding protein-like protein 10. It is a member of OSBP family and plays a key role in the maintenance of the cholesterol balance [69]. Gene ZNF860 encodes zinc finger protein 860. Our network analysis demonstrated that five genes of this set were functionally related to pathways which enable tissue formation and intercellular communications: “adherens junction,” “focal adhesion,” and “tight junction.” We could see that these statistically significant relations (NEA FDR < 0.01; Figure 5(a)) affected the subnetwork at multiple points and via different molecular mechanisms (e.g., VHL is a tumor suppressor, LMCD1 is a transcription factor, and CMTM6 is a regulator of cytokine signaling). Hence the role of these genes/loci in cancer and, particularly in prostate cancer, warrants further investigation.

(a)

(b)

Figure 5: NEA: network analysis of the sets of potential markers in regard of individual genes and KEGG pathways. (a) Detailed view of network connections between genes of Table (yellow) and genes of three (partially overlapping) KEGG pathways “adherens junction,” “focal adhesion,” and “tight junction” (red). Many genes of these pathways also had experimentally verified binding sited to LMCD1 (HTRIdb database) which are not shown here for the sake of simplicity. The network enrichment of the red genes as a whole set, on the one hand, against the genes of Table , on the other hand, was probabilistically evaluated and is part of the more general figure (b). (b) Generalized view on functional relations between novel sets of potential markers (yellow) and KEGG pathways, based on the network enrichment analysis. Numeric edge labels denote the number of individual gene-gene links behind each relation. Arrow opacity reflects statistical confidence of relations (although each relation is based on at least 10 individual gene-gene links and has false discovery rate <0.01). The shades of red and brown reflect the cumulative connectivity of the KEGG pathways in the global network.

The comparison of nonaggressive prostate carcinomas () with prostate adenomas permitted us to identify six genes/loci with statistically significant association of methylation/deletion events (): CAND2, GATA2, FAM19A4, KY, ALDH1L1, and MAP4 (Table 3).

Table 3: Methylation/deletion frequencies in prostate adenoma versus nonaggressive cancer.

Gene
Methylation/deletion frequency, %
value

Adenoma
Cancer

CAND2
47 (7/15)
0 (0/14)
0.006

GATA2
20 (3/15)
71 (10/14)
0.009

FAM19A4
40 (6/15)
0 (0/14)
0.017

KY
60 (9/15)
14 (2/14)
0.021

ALDH1L1
33 (5/15)
79 (11/14)
0.025

MAP4
33 (5/15)
0 (0/14)
0.042

Figure 5(b) displays that this small set was highly enriched with functional links to a number of metabolic pathways related to glycan and keratin synthesis as well as to the cytokine signaling pathway. A detailed view (not shown) revealed that these enrichment patterns were enabled mostly by the presence of binding sites for the transcription factor GATA2 in the promoters of respective genes. The binding sites were predicted from sequence (http://genome.ucsc.edu/) and validated experimentally (HTRIdb database). To a much lesser extent, the enrichment was enabled by the presence of functional links of KY and MAP4 in the FunCoup network [49].

The resulting set also comprised five markers: FAM19A4, CAND2, MAP4, KY, and LRRC58. A sample would be recognized as cancerous if methylation and/or deletion were found in less than two of these markers, which yielded sensitivity 93% and specificity 73%. Gini coefficient was in the range 0.61–0.87 (Figure 4(c)).

Gene FAM19A4 (3p14.1) encodes chemokine-like protein TAFA4. The role of this gene in prostate cancer is not clear, but methylation of FAM19A4 was found in cervical [70] and breast tumors [71]. It may indicate that FAM19A4 is potential tumor suppressor gene. Gene CAND2 (3p25.2) encodes cullin-associated and neddylation-dissociated 2 protein. It may be a transcription regulator and play a role in the assembly of ubiquitin ligase complexes [72]. Its role in cancerogenesis needs investigation. Gene MAP4 (3p21) encodes microtubule-associated protein 4, which stabilizes microtubules and controls their dynamics in mitosis. Its functions in prostate cancer have yet to be fully elucidated. Gene KY (3q22.2) encodes kyphoscoliosis peptidase. Its function in prostate cancer is not clear. Gene LRRC58 located in 3q13.33 region and encodes leucine-rich repeat-containing protein 58. There is a few information about this gene.

In summary, the suggested set of 16 markers (BHLHE40, FOXP1, LOC285205, ITGA9, CTDSPL, FGF12, LOC440944/SETD5, VHL, CLCN2, OSBPL10/ZNF860, LMCD1, FAM19A4, CAND2, MAP4, KY, and LRRC58) allowed for discriminating/diagnosing the majority of prostate tumor cases with accuracy more than 83% for examined sampling (Table 4). This is one of the possible classifiers that could be constructed from the data. We selected the set based on the reliability, biological interpretation, and maximized statistical significance given the current data. The prediction power of the developed sets should be further validated on additional collections of prostate biopsy samples.

Table 4: Detection and discrimination of aggressive/nonaggressive prostate cancer and adenoma using the set of 16 selected markers.

Use
Set of markers

Detection of aggressive/nonaggressive cancer and adenoma
BHLHE40, FOXP1, LOC285205, ITGA9, CTDSPL (RBSP3), FGF12  
Sp = 94% Sn = 94% Ac = 94%

Discrimination between aggressive cancer and nonaggressive cancer/adenoma
LOC440944/SETD5, VHL, CLCN2, OSBPL10/ZNF860, LMCD1  
Sp = 97% Sn = 100% Ac = 97%

Discrimination between nonaggressive cancer and adenoma
FAM19A4, CAND2, MAP4, KY, LRRC58  
Sp = 73% Sn = 93% Ac = 83%

Note. Sp/Sn/Ac: specificity/sensitivity/accuracy of the set.

4. Conclusions

Our study shows that alterations on chromosome 3 often accompany formation of prostate tumors. Fifty genes with frequent (>30%) methylation/deletion aberrations in prostate tumors with different pathomorphological characteristics were revealed using NotI-microarray technology. Frequent methylation of FGF12, GATA2, and LMCD1 genes was confirmed by bisulfite sequencing. Hypermethylation of BHLHE40, BCL6, and ITGA9 genes was associated with their downregulation according to qPCR analysis. Six of fifty genes were not previously known to associate with cancer (HMGB1L5, LRRC58, GPR149, DZIP1L, C3orf77, and NUDT16). Further analysis of alterations in signaling pathways involving these genes is of high interest. According to our and other authors’ works seventeen genes/loci (PLCL2, IQSEC1, ZIC4, ALDH1L1, WNT7A, KY, PPP2R3A, GATA2, LOC285205, NKIRAS1, ITGA9, CGGBP1, FOXP1, GORASP1, NBEAL2, RBSP3, and LRRC3B) are associated with cancers of other localizations.

We investigated functional coherence of the presented gene sets using the network enrichment analysis. Beyond the informative functional relations we employed the ability of NEA to evaluate internal consistency of a given set. We established that the discovered sets were enriched in connection with each other, that is, within the respective sets (NEA FDR < 0.05 in each case). In other words, each of these observations was highly unlikely to be made by chance, that is, in a random set of genes or in a random network (given the genes and the network of same topological properties).

These results approved NotI-microarray technology as a powerful method for screening of epigenetic and genetic alterations in prostate cancer. The final set of 16 promising markers for detection of prostate tumors and discrimination between prostate adenoma and carcinomas with different aggressiveness is suggested for further studies and refinement.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

Alexey A. Dmitriev and Eugenia E. Rosenberg contributed equally to this work.

Acknowledgments

This work was supported by Grant 0110U004744 from National Academy of Sciences of Ukraine; Grants 14-04-31978 mol_a, 13-04-01885 a, and 15-04-08731 a from the Russian Foundation for Basic Research; the grant from RAS Presidium Program “Molecular and Cellular Biology.” Part of this work (qPCR analysis) was performed using the equipment of EIMB RAS “Genome” center (http://www.eimb.ru/RUSSIAN_NEW/INSTITUTE/ccu_genome_c.php) under the financial support by the Ministry of Education and Science of the Russian Federation (Contract 14.621.21.0001, project’s unique identifier RFMEFI62114X0001). The support by BILS (Bioinformatics Infrastructure for Life Sciences in Sweden) is gratefully acknowledged.

References

J. Ferlay, I. Soerjomataram, M. Ervik et al., GLOBOCAN 2012 v1.0, Cancer Incidence and Mortality Worldwide: IARC CancerBase No. 11, International Agency for Research on Cancer, Lyon, France, 2012, http://globocan.iarc.fr/.

H. Ahmadi and S. Daneshmand, “Androgen deprivation therapy for prostate cancer: long-term safety and patient outcomes,” Patient Related Outcome Measures, vol. 5, pp. 63–70, 2014.

M. Katsogiannou, H. Ziouziou, S. Karaki, C. Andrieu, M. H. de Villeneuve, and P. Rocchi, “The hallmarks of castration-resistant prostate cancers,” Cancer Treatment Reviews, vol. 41, no. 7, pp. 588–597, 2015.

A. D. Grigore, E. Ben-Jacob, and M. C. Farach-Carson, “Prostate cancer and neuroendocrine differentiation: more neuronal, less endocrine?” Frontiers in Oncology, vol. 5, article 37, 2015.

V. J. Gnanapragasam, “Molecular markers to guide primary radical treatment selection in localized prostate cancer,” Expert Review of Molecular Diagnostics, vol. 14, no. 7, pp. 871–881, 2014.

S. M. Nazim and F. Abbas, “Role of surgery in locally advanced prostate cancer,” Pakistan Journal of Medical Sciences, vol. 31, no. 3, pp. 710–716, 2015.

E. F. Faria, B. F. Chapin, R. L. Muller, R. D. Machado, R. B. Reis, and S. F. Matin, “Radical prostatectomy for locally advanced prostate cancer: current status,” Urology, vol. 86, no. 1, pp. 10–15, 2015.

J. Sullivan, R. Kopp, K. Stratton et al., “An analysis of the association between prostate cancer risk loci, PSA levels, disease aggressiveness and disease-specific mortality,” British Journal of Cancer, vol. 113, no. 1, pp. 166–172, 2015.

M. Danilewicz, O. Stasikowska-Kanicka, and M. Wagrowska-Danilewicz, “Augmented immunoexpression of survivin correlates with parameters of aggressiveness in prostate cancer,” Polish Journal of Pathology, vol. 66, no. 1, pp. 44–48, 2015.

S. I. Berndt, Z. Wang, M. Yeager et al., “Two susceptibility loci identified for prostate cancer aggressiveness,” Nature Communications, vol. 6, p. 6889, 2015.

A. B. Gutman and E. B. Gutman, “An ‘acid’ phosphatase occurring in the serum of patients with metastasizing carcinoma of the prostate gland,” The Journal of Clinical Investigation, vol. 17, no. 4, pp. 473–478, 1938.

T. Bhavsar, P. McCue, and R. Birbe, “Molecular diagnosis of prostate cancer: are we up to age?” Seminars in Oncology, vol. 40, no. 3, pp. 259–275, 2013.

H. Y. Kong and J. Byun, “Emerging roles of human prostatic acid phosphatase,” Biomolecules and Therapeutics, vol. 21, no. 1, pp. 10–20, 2013.

R. Chou, J. M. Croswell, T. Dana et al., “Screening for prostate cancer: a review of the evidence for the U.S. Preventive services task force,” Annals of Internal Medicine, vol. 155, no. 11, pp. 762–771, 2011.

M. Ferro, D. Bruzzese, S. Perdonà et al., “Prostate Health Index (Phi) and Prostate Cancer Antigen 3 (PCA3) significantly improve prostate cancer detection at initial biopsy in a total PSA range of 2–10 ng/ml,” PLoS ONE, vol. 8, no. 7, Article ID e67687, 2013.

E. E. Rosenberg, T. Y. Prudnikova, G. V. Gerashchenko, E. V. Grigorieva, and V. I. Kashuba, “Search for genes—potential markers of aggressiveness and metastasis for human prostate cancer,” Biopolymers and Cell, vol. 29, no. 6, pp. 499–505, 2013.

C. A. Iacobuzio-Donahue, “Epigenetic changes in cancer,” Annual Review of Pathology: Mechanisms of Disease, vol. 4, pp. 229–249, 2009.

R. Singal, L. Ferdinand, I. M. Reis, and J. J. Schlesselman, “Methylation of multiple genes in prostate cancer and the relationship with clinicopathological features of disease,” Oncology Reports, vol. 12, no. 3, pp. 631–637, 2004.

J. Pan, Y. Chen, C. Mo et al., “Association of DSC3 mRNA down-regulation in prostate cancer with promoter hypermethylation and poor prognosis,” PLoS ONE, vol. 9, no. 3, Article ID e92815, 2014.

J. Sandoval, L. Peiró-Chova, F. V. Pallardó, and J. L. García-Giménez, “Epigenetic biomarkers in laboratory diagnostics: emerging approaches and opportunities,” Expert Review of Molecular Diagnostics, vol. 13, no. 5, pp. 457–471, 2013.

S. H. Strand, T. F. Orntoft, and K. D. Sorensen, “Prognostic DNA methylation markers for prostate cancer,” International Journal of Molecular Sciences, vol. 15, no. 9, pp. 16544–16576, 2014.

K. Daniunaite, S. Jarmalaite, N. Kalinauskaite et al., “Prognostic value of RASSF1 promoter methylation in prostate cancer,” Journal of Urology, vol. 192, no. 6, pp. 1849–1855, 2014.

K. Woodson, K. J. O'Reilly, D. E. Ward et al., “CD44 and PTGS2 methylation are independent prognostic markers for biochemical recurrence among prostate cancer patients with clinically localized disease,” Epigenetics, vol. 1, no. 4, pp. 183–186, 2006.

K. Chiam, C. Ricciardelli, and T. Bianco-Miotto, “Epigenetic biomarkers in prostate cancer: current and future uses,” Cancer Letters, vol. 342, no. 2, pp. 248–256, 2014.

G. Weiss, S. Cottrell, J. Distler et al., “DNA methylation of the PITX2 gene promoter region is a strong independent prognostic marker of biochemical recurrence in patients with prostate cancer after radical prostatectomy,” The Journal of Urology, vol. 181, no. 4, pp. 1678–1685, 2009.

L. Maldonado, M. Brait, M. Loyo et al., “GSTP1 promoter methylation is associated with recurrence in early stage prostate cancer,” Journal of Urology, vol. 192, no. 5, pp. 1542–1548, 2014.

W. Zhang, H. Jiao, X. Zhang et al., “Correlation between the expression of DNMT1, and GSTP1 and APC, and the methylation status of GSTP1 and APC in association with their clinical significance in prostate cancer,” Molecular Medicine Reports, vol. 12, no. 1, pp. 141–146, 2015.

L. Richiardi, V. Fiano, C. Grasso et al., “Methylation of APC and GSTP1 in non-neoplastic tissue adjacent to prostate tumour and mortality from prostate cancer,” PLoS ONE, vol. 8, no. 7, Article ID e68162, 2013.

W. Li, R. Wang, Z. Yan, L. Bai, and Z. Sun, “High accordance in prognosis prediction of colorectal cancer across independent datasets by multi-gene module expression profiles,” PLoS ONE, vol. 7, no. 3, Article ID e33653, 2012.

R. Tsunashima, Y. Naoi, N. Kagara et al., “Construction of multi-gene classifier for prediction of response to and prognosis after neoadjuvant chemotherapy for estrogen receptor positive breast cancers,” Cancer Letters, vol. 365, no. 2, pp. 166–173, 2015.

G. V. Gerashchenko, O. O. Bogatyrova, E. E. Rudenko et al., “Genetic and epigenetic changes of NKIRAS1 gene in human renal cell carcinomas,” Experimental Oncology, vol. 32, no. 2, pp. 71–75, 2010.

V. V. Gordiyuk, “Genetic and epigenetic changes of genes on chromosome 3 in human urogenital tumors,” Biopolymers and Cell, vol. 27, no. 1, pp. 25–35, 2011.

V. I. Kashuba, I. Y. Skrypkina, D. V. Saraev et al., “Identification of changes in gene loci potentially associated with cervical cancer using NotI microarrays,” Ukrain'skyi Biokhimichnyi Zhurnal, vol. 78, no. 2, pp. 113–120, 2006.

L. B. Hesson, W. N. Cooper, and F. Latif, “Evaluation of the 3p21.3 tumour-suppressor gene cluster,” Oncogene, vol. 26, no. 52, pp. 7283–7301, 2007.

G. Ciriello, M. L. Miller, B. A. Aksoy, Y. Senbabaoglu, N. Schultz, and C. Sander, “Emerging landscape of oncogenic signatures across human cancers,” Nature Genetics, vol. 45, no. 10, pp. 1127–1133, 2013.

L. Ji, J. D. Minna, and J. A. Roth, “3p21.3 tumor suppressor cluster: prospects for translational applications,” Future Oncology, vol. 1, no. 1, pp. 79–92, 2005.

A. J. H. Smith, J. Xian, M. Richardson, K. A. Johnstone, and P. Rabbitts, “Cre-loxP chromosome engineering of a targeted deletion in the mouse corresponding to the 3p21.3 region of homozygous loss in human tumours,” Oncogene, vol. 21, no. 29, pp. 4521–4529, 2002.

J. Li, A. Protopopov, F. Wang et al., “NotI subtraction and NotI-specific microarrays to detect copy number and methylation changes in whole genomes,” Proceedings of the National Academy of Sciences of the United States of America, vol. 99, no. 16, pp. 10724–10729, 2002.

T. Maniatis, E. F. Fritsch, and J. Sambrook, Molecular Cloning: A Laboratory Manual, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA, 1982.

T. L. Dunwell, L. B. Hesson, T. Pavlova et al., “Epigenetic analysis of childhood acute lymphoblastic leukemia,” Epigenetics, vol. 4, no. 3, pp. 185–193, 2009.

V. V. Gordiyuk, G. V. Gerashchenko, I. Y. Skrypkina et al., “Identification of chromosome 3 epigenetic and genetic abnormalities and gene expression changes in ovarian cancer,” Biopolymers and Cell, vol. 24, no. 4, pp. 323–332, 2008.

K. Haraldson, V. I. Kashuba, A. A. Dmitriev et al., “LRRC3B gene is frequently epigenetically inactivated in several epithelial malignancies and inhibits cell growth and replication,” Biochimie, vol. 94, no. 5, pp. 1151–1157, 2012.

A. A. Dmitriev, V. I. Kashuba, K. Haraldson et al., “Genetic and epigenetic analysis of non-small cell lung cancer with NotI-microarrays,” Epigenetics, vol. 7, no. 5, pp. 502–513, 2012.

V. N. Senchenko, G. S. Krasnov, A. A. Dmitriev et al., “Differential expression of CHL1 gene during development of major human cancers,” PLoS ONE, vol. 6, no. 3, Article ID e15612, 2011.

G. S. Krasnov, N. Y. Oparina, A. A. Dmitriev et al., “RPN1, a new reference gene for quantitative data normalization in lung and kidney cancer,” Molecular Biology, vol. 45, no. 2, pp. 211–220, 2011.

C. Campbell and Y. Ying, Learning with Support Vector Machines, Morgan & Claypool Publishers, San Rafael, Calif, USA, 2011.

R. Nisbet, J. Elder, and G. Miner, Handbook of Statistical Analysis and Data Mining Applications, Academic Press, Burlington, Mass, USA, 2009.

A. Alexeyenko, W. Lee, M. Pernemalm et al., “Network enrichment analysis: extension of gene-set enrichment analysis to gene networks,” BMC Bioinformatics, vol. 13, no. 1, article 226, 2012.

A. Alexeyenko and E. L. L. Sonnhammer, “Global networks of functional coupling in eukaryotes from comprehensive data integration,” Genome Research, vol. 19, no. 6, pp. 1107–1116, 2009.

M. Kanehisa, S. Goto, S. Kawashima, and A. Nakaya, “The KEGG databases at GenomeNet,” Nucleic Acids Research, vol. 30, no. 1, pp. 42–46, 2002.

A. Ruepp, B. Waegele, M. Lechner et al., “CORUM: the comprehensive resource of mammalian protein complexes—2009,” Nucleic Acids Research, vol. 38, pp. D497–D501, 2010.

A. Subramanian, P. Tamayo, V. K. Mootha et al., “Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 43, pp. 15545–15550, 2005.

L. A. Bovolenta, M. L. Acencio, and N. Lemke, “HTRIdb: an open-access database for experimentally verified human transcriptional regulation interactions,” BMC Genomics, vol. 13, article 405, 2012.

E. R. Zabarovsky, R. Gizatullin, R. M. Podowski et al., “NotI clones in the analysis of the human genome,” Nucleic Acids Research, vol. 28, no. 7, pp. 1635–1639, 2000.

G. V. Gerashchenko, V. V. Gordiyuk, I. Y. Skrypkina et al., “Screening of epigenetic and genetic disturbances of human chromosome 3 genes in colorectal cancer,” Ukrainskiĭ Biokhimicheskiĭ Zhurna, vol. 81, no. 4, pp. 81–87, 2009.

V. N. Senchenko, N. P. Kisseljova, T. A. Ivanova et al., “Novel tumor suppressor candidates on chromosome 3 revealed by NotI-microarrays in cervical cancer,” Epigenetics, vol. 8, no. 4, pp. 409–420, 2013.

A. A. Dmitriev, E. E. Rudenko, A. V. Kudryavtseva et al., “Epigenetic alterations of chromosome 3 revealed by noti-microarrays in clear cell renal cell carcinoma,” BioMed Research International, vol. 2014, Article ID 735292, 9 pages, 2014.

A. G. Kondratov, S. M. Kvasha, L. A. Stoliar et al., “Alterations of the WNT7A gene in clear cell renal cell carcinomas,” PLoS ONE, vol. 7, no. 10, Article ID e47012, 2012.

V. Kashuba, A. A. Dmitriev, G. S. Krasnov et al., “Noti microarrays: novel epigenetic markers for early detection and prognosis of high grade serous ovarian cancer,” International Journal of Molecular Sciences, vol. 13, no. 10, pp. 13352–13377, 2012.

H. Li, Y. Du, D. Zhang et al., “Identification of novel DNA methylation markers in colorectal cancer using MIRA-based microarrays,” Oncology Reports, vol. 28, no. 1, pp. 99–104, 2012.

Z.-Q. Lian, Q. Wang, W.-P. Li, A.-Q. Zhang, and L. Wu, “Screening of significantly hypermethylated genes in breast cancer using microarray-based methylated-CpG island recovery assay and identification of their expression levels,” International Journal of Oncology, vol. 41, no. 2, pp. 629–638, 2012.

J. U. Kang, S. H. Koo, K. C. Kwon, J. W. Park, and J. M. Kim, “Identification of novel candidate target genes, including EPHB3, MASP1 and SST at 3q26.2-q29 in squamous cell carcinoma of the lung,” BMC Cancer, vol. 9, article 237, 2009.

T. V. Pavlova, V. I. Kashuba, O. V. Muravenko et al., “Use of NotI microarrays in analysis of epigenetic and structural changes in epithelial tumor genomes by the example of human chromosome 3,” Molecular Biology, vol. 43, no. 2, pp. 313–320, 2009.

B. S. Taylor, N. Schultz, H. Hieronymus et al., “Integrative genomic profiling of human prostate cancer,” Cancer Cell, vol. 18, no. 1, pp. 11–22, 2010.

P. Bastola, Y. Stratton, E. Kellner et al., “Folliculin contributes to VHL tumor suppressing activity in renal cancer through regulation of autophagy,” PLoS ONE, vol. 8, no. 7, Article ID e70030, 2013.

M. A. Chetram, D. A. Bethea, V. A. Odero-Marah, A. S. Don-Salu-Hewage, K. J. Jones, and C. V. Hinton, “ROS-mediated activation of AKT induces apoptosis via pVHL in prostate cancer cells,” Molecular and Cellular Biochemistry, vol. 376, no. 1-2, pp. 63–71, 2013.

X.-Y. Yang, X.-G. Lai, Y. Zhang, J.-M. Pei, A.-G. Yang, and S.-S. Zhou, “siRNA-mediated silencing of ClC-2 gene inhibits proliferation of human U-87 glioma cells,” Ai Zheng, vol. 25, no. 7, pp. 805–810, 2006.

C.-Y. Chang, S.-C. Lin, W.-H. Su, C.-M. Ho, and Y.-S. Jou, “Somatic LMCD1 mutations promoted cell migration and tumor metastasis in hepatocellular carcinoma,” Oncogene, vol. 31, no. 21, pp. 2640–2652, 2012.

M. Lehto, S. Laitinen, G. Chinetti et al., “The OSBP-related protein family in humans,” The Journal of Lipid Research, vol. 42, no. 8, pp. 1203–1213, 2001.

R. D. M. Steenbergen, M. Ongenaert, S. Snellenberg et al., “Methylation-specific digital karyotyping of HPV16E6E7-expressing human keratinocytes identifies novel methylation events in cervical carcinogenesis,” Journal of Pathology, vol. 231, no. 1, pp. 53–62, 2013.

K. Miyamoto, T. Fukutomi, S. Akashi-Tanaka et al., “Identification of 20 genes aberrantly methylated in human breast cancers,” International Journal of Cancer, vol. 116, no. 3, pp. 407–414, 2005.

S. Shiraishi, C. Zhou, T. Aoki et al., “TBP-interacting protein 120B (TIP120B)/cullin-associated and neddylation-dissociated 2 (CAND2) inhibits SCF-dependent ubiquitination of myogenin and accelerates myogenic differentiation,” Journal of Biological Chemistry, vol. 282, no. 12, pp. 9017–9028, 2007.