- Open Access
Factor interaction analysis for chromosome 8 and DNA methylation alterations highlights innate immune response suppression and cytoskeletal changes in prostate cancer
Molecular Cancer volume 6, Article number: 14 (2007)
Alterations of chromosome 8 and hypomethylation of LINE-1 retrotransposons are common alterations in advanced prostate carcinoma. In a former study including many metastatic cases, they strongly correlated with each other. To elucidate a possible interaction between the two alterations, we investigated their relationship in less advanced prostate cancers.
In 50 primary tumor tissues, no correlation was observed between chromosome 8 alterations determined by comparative genomic hybridization and LINE-1 hypomethylation measured by Southern blot hybridization. The discrepancy towards the former study, which had been dominated by advanced stage cases, suggests that both alterations converge and interact during prostate cancer progression. Therefore, interaction analysis was performed on microarray-based expression profiles of cancers harboring both alterations, only one, or none. Application of a novel bioinformatic method identified Gene Ontology (GO) groups related to innate immunity, cytoskeletal organization and cell adhesion as common targets of both alterations. Many genes targeted by their interaction were involved in type I and II interferon signaling and several were functionally related to hereditary prostate cancer genes. In addition, the interaction appeared to influence a switch in the expression pattern of EPB41L genes encoding 4.1 cytoskeleton proteins. Real-time RT-PCR revealed GADD45A, MX1, EPB41L3/DAL1, and FBLN1 as generally downregulated in prostate cancer, whereas HOXB13 and EPB41L4B were upregulated. TLR3 was downregulated in a subset of the cases and associated with recurrence. Downregulation of EPB41L3, but not of GADD45A, was associated with promoter hypermethylation, which was detected in 79% of carcinoma samples.
Alterations of chromosome 8 and DNA hypomethylation in prostate cancer probably do not cause each other, but converge during progression. The present analysis implicates their interaction in innate immune response suppression and cytoskeletal changes during prostate cancer progression. The study thus highlights novel mechanisms in prostate cancer progression and identifies novel candidate genes for diagnostic and therapeutic purposes. In particular, TLR3 expression might be useful for prostate cancer prognosis and EPB41L3 hypermethylation for its detection.
Up to 40% of all elderly men may harbor prostate carcinomas, less than 20% develop symptomatic disease, which in about 3% becomes the cause of death. Great advances have been made in prostate cancer detection due to refined PSA-based assays, imaging and histopathology. In their wake, two questions have gained importance, i.e. which tumors represent clinically significant disease and how tumors having progressed locally or metastasized can be identified and appropriately treated. Molecular research is thus challenged to provide a reliable classification of prostate cancers and to identify targets for novel therapies in those tumors no longer containable by surgery, irradiation, and anti-androgenic therapy.
Among such aggressive prostate cancers, molecular alterations are not uniform. Instead, these tumors contain different combinations of genetic and epigenetic aberrations that each appear to influence the biological processes crucial for the cancer phenotype [1–3]. It is therefore necessary to disentangle the influence of individual alterations by factor analysis approaches. Through defining the biological consequences of particular recurrent molecular alterations – considered as factors – it may become possible to identify subgroups of prostate carcinoma whose behavior is determined by them. These factors could then be used to predict the prognosis in each particular subgroup. In addition, elucidation of the mechanisms by which a factor exerts its effects on prostate cancer progression would provide specific targets for therapy within the respective subgroup. This approach does not imply that factors act independently of each other. Indeed, we suggest here that two molecular alterations previously identified as being associated with more aggressive prostate cancers appear to interact in a synergistic fashion during tumor progression.
The first factor, alteration of chromosome 8, is found in up to 50% of prostate cancers [1–3]. A large number of publications have concordantly reported that chromosome 8 alterations are significantly associated with various established histopathological indicators of poor prognosis or directly with clinical outcome. Curiously, 8p losses or 8q gains both seem to have similar clinical effects (reviewed in ). One reason for this unusual relationship is that 8q gain is often mechanistically linked to 8p loss through isochromosome 8q formation . Moreover, it has been difficult to trace the effect of chromosome 8 alterations to any single gene . Interestingly, chromosome 8 has also emerged as the site of a hereditary prostate cancer gene in genome-wide searches. A good candidate is MSR1 at 8p23.1. Inherited mutations at this locus are presumed to modify immune responses during carcinogenesis . Thus, the importance of chromosome 8 alterations as a factor in prostate cancer progression is undisputed, but the question of which biological processes they influence is open.
The second factor is hypomethylation of LINE-1 retrotransposons, which constitute approximately 18% of the human genome. LINE-1 hypomethylation is part of a broader process, 'genome-wide' hypomethylation that affects not only repeat sequences, but also single-copy genes inactivated by DNA methylation in normal adult tissues . While hypomethylation occurs at early stages in some cancer types, it is associated with tumor progression in others. Conceivably, this association is brought about by increased chromosomal instability as a consequence of repetitive sequence hypomethylation and by reactivation of specific genes favoring tumor cell survival and adaptation during invasion and metastasis . However, the details of the relationship between genome-wide hypomethylation and cancer progression are poorly understood. In prostate cancer, LINE-1 hypomethylation is detected in only a subset of localized cancers, but is highly prevalent in high-stage and metastatic cases [8, 9]. In contrast, hypermethylation of selected single copy genes occurs consistently at early stages of carcinogenesis [9, 10].
The two factors, chromosome 8 alterations and LINE-1 hypomethylation, may not be independent of each other. In a previous study of 55 prostate carcinomas including many advanced stage and recurrent cases, we found LINE-1 hypomethylation to be highly significantly associated with chromosomal instability , in keeping with expectations from other cancer types and animal models [6, 7]. Unexpectedly, we found a particularly close relationship between hypomethylation and chromosome 8 alterations. This raised the question whether this association was (i) mechanistic, i.e. a consequence of chromosome 8 destabilization by hypomethylation or, conversely, deregulation of DNA methylation as a consequence of chromosome 8 changes, (ii) the consequence of a common factor causing both alterations, or (iii) an indication of a functional interaction between both factors that promotes prostate cancer progression. Here, we will present evidence that explanation (iii) is most likely. Our analysis revealed that the interaction of the two factors appears to relate to several biological processes already implicated in prostate cancer progression, prominently suppression of immune responses and changes in the cortical actin cytoskeleton and the extracellular matrix.
Comparative genomic hybridization and DNA methylation analysis
Our previous study had revealed a highly significant association between LINE-1 hypomethylation and chromosome 8 alterations, particularly in tumors with distant metastases or recurring after anti-androgenic therapy . To further elucidate this association, we investigated whether it also existed in less advanced cancers. From a previous study of DNA methylation alterations in prostate carcinoma , a subset of 50 primary tumors without distant metastases were selected by the criteria of availability of high-quality DNA and RNA and complete follow-up data (Table 1). Almost all tumor tissues presented hypermethylation in several genes typically hypermethylated in prostate cancer (Table 1).
For all cases, comparative genomic hybridization (CGH) analysis was performed using the same DNA aliquots employed in methylation analysis (Table 1). Overall, 17 tumors revealed typical changes in chromosome 8 dosage. Loss of 8p was detected in 6 cases and gain of 8q in 13 cases. Three tumors exhibiting nontypical 8p gain or 8q loss were not categorized as 'chromosome 8 altered'.
Significant LINE-1 hypomethylation (>4%) was found in 17 tumor specimens (Table 1). Only five of these harbored 8p losses or 8q gains, 12 tumors each contained either hypomethylation or chromosome 8 alterations, and in 21 cases neither alteration was detected. Thus, in the present series no association was found between LINE-1 hypomethylation and chromosome 8 alterations. The previously measured association therefore likely reflects a convergence in late stage prostate cancers, i.e. both molecular alterations may be selected for during progression of prostate cancers towards an aggressive phenotype. This co-selection is most straightforwardly explained by a synergistic influence of the two factors on biological processes important for tumor progression.
Microarray expression analysis
To identify which biological processes might be targeted by the supposed synergism, expression profiling was performed. From those cases for which ample high quality RNA was available, 24 tumors were selected exhibiting only chromosome 8 alterations (8p loss, 8q gain, or both), only LINE-1 hypomethylation, chromosome 8 alterations as well as LINE-1 hypomethylation ('both' group), or neither alteration (Table 1). The specimens lacking both changes were selected to harbor a maximum of other chromosomal alterations. Total RNA from each case was hybridized to Affymetrix U133A oligonucleotide arrays, without pre-amplification.
This experimental setup allows three comparisons, i.e. (i) cancers with or without chromosome 8 alterations, (ii) cancers with or without LINE-1 hypomethylation, and (iii) an interaction analysis of both factors. In the interaction analysis, we considered those genes as significant, whose changes in expression were significantly greater in 'both' cancers than the sum of the changes in cancers with one alteration, each compared to the cancers without either alteration (see Methods). Moreover, since we were interested in the identification of biological processes targeted by the presumed interaction, all comparisons were performed for Gene Ontology (GO) 'biological process' groups (see Methods). In this analysis, a Gene Ontology group is highlighted if it is enriched, i.e. if it contains statistically significantly more differentially expressed genes than expected by chance.
Intriguingly, few differences were discovered in the comparisons based upon differential methylation or upon the presence of chromosome 8 alterations (Figure 1A). In contrast, the interaction analysis revealed a number of highly significant differences (Figure 1A). Initially, the classic algorithm (see Methods) for testing enrichment was used. The majority of the significant terms obtained with this algorithm were related to 'immune response' and 'response towards wounding' (areas 1. and 2. in Figure 1B). A prominent unrelated GO group was #30865 (area 7. in Figure 1B). Individual genes significant within GO groups are listed in Additional file 1 (further information on Figures 1B and 1C can be found in Additional file 3).
Expression of immune response-related genes
Several genes individually significant in the interaction analysis were selected for closer analysis. The microarray results for these and some related other genes are shown in Figure 2. Figure 3 presents a comparison of their expression analyzed by quantitative RT-PCR in 47 cancer samples (including almost all used in the microarray analysis) and 13 morphologically normal prostate tissues obtained from distant sites of cancer-carrying prostates. This latter analysis allows to determine whether these genes were differentially expressed not only among prostate cancer groups, but also between prostate cancers and noncancerous tissue in general. As a quality control, HPN (Hepsin) expression was first determined in all these samples and was strongly increased in cancerous over benign tissues (Figure 3A), as expected [14–17]. As another quality control, DNA from the benign tissues used exhibited neither GSTP1 nor RARB2 hypermethylation, excluding contamination with tumor [9, 10].
Like Hepsin, HOXB13 expression was overall significantly elevated in cancer tissue (Figure 3B), although the increase was highly variable. This variation was obviously related to the presence of chromosome 8 alterations and LINE-1 hypomethylation (Figure 2A). No other classical homeobox gene was significant, including the closest HOXB13 paralogs, HOXC13 and HOXD13 (data not shown). The upregulation of HOXB13, but also its variability are in keeping with previous studies [17, 19].
In contrast, BCCIP was identified as significant (p = 0.0036) in the interaction analysis, with maximum expression in the 'both' group (Figure 2B), but its expression did not differ significantly between cancer and noncancerous tissues overall (Figure 3C). BCCIP interacts with BRCA2. Comparisons of BRCA1, BRCA2, and BARD1 did not yield differences (data not shown), but another gene interacting with BRCA2, HMG20B, was significant in the interaction analysis (Figure 2C). Altered expression of HMG20B has previously been reported in prostate cancer .
GADD45A, which is also induced by BRCA-dependent pathways, was significantly decreased in prostate cancers compared to noncancerous tissues (Figure 3D), with the lowest levels in 'both' cancers (Figure 2D). Expression of GADD45B or GADD45G was not significantly different between the tumor groups.
Several genes related to interferon signaling were downregulated in the 'both' cancers, most significantly, the prototypic type I interferon response gene MX1 ( also named MXA) (Figure 2E). MX2 (MXB) behaved similarly (Figure 2F). MX1 was also downregulated in cancer tissues overall (Figure 3E). A significant decrease in the 'both' group was also observed for CCL5 encoding a chemokine enhancing interferon responses (Figure 2G). Expression of IFNB1 (encoding IFNβ), but not of IFNG or IFNA genes, was significantly lowered according to the interaction analysis (Figure 2H), and so was accordingly IRF3 (Figure 2I) encoding a major transcriptional activator of IFNB1. Another important transcription factor for type I interferons, STAT3, showed minor decreases (Figure 2J). Interestingly, selected TLRs were downregulated in 'both' cancers, especially TLR2 (Figure 2K) and TLR3 (Figure 2L). TLR3 expression was decreased in many, but not all prostate cancers compared to normal tissues (Figure 3F).
The expression differences in immune-response genes might potentially be due to altered proportions of innate immune cells in the different tumor groups. We therefore analyzed, whether several typical marker genes of NK cells, mast cells, macrophages, dendritic cells, and granulocytes were differentially expressed in the 'both' cancers [see Additional File 2]. Only a single marker, ITGAM (more commonly known as MAC-1) showed a significant, but small difference between the groups. From this analysis, the differences in the proportions of innate immune cells between the tumor groups appear to be small and cannot account for the overall observed expression differences.
Expression of cytoskeleton and extracellular matrix genes
Since the classic algorithm for detecting enriched GO terms treats each GO term independently, related biological terms often appear simultaneously among the top scoring significant GO terms. In the more sophisticated weight algorithm, genes annotated to a GO term receive weights based on the significance of neighboring GO terms (see Methods). This approach identifies local dependencies between related GO terms and highlights those terms that receive a more significant enrichment score than all their neighbors in the GO graph. Upon application of the weight algorithm the most significant GO terms are spread over additional areas in the GO graph (compare Figure 1B and 1C), identifying a larger variety of biological processes (Table 2), in addition to those related to immune response. Most significantly, the weight algorithm highlights another set of GO groups concerning the interacting processes of cortical actin cytoskeletal organization and cell adhesion (areas 3. and 7. in Figure 1C). In fact, the significance of GO groups 'regulation of embryonic development' and 'positive regulation of signal transduction' (areas 6. and 8. in Figure 1C) also derives predominantly from genes related to cell adhesion [see Additional file 1].
The conspicuous GO group #30865 comprises mostly genes encoding band 4.1 proteins, officially termed EPB41L. Most family members were down-regulated in the 'both' cancers, especially EPB41L3 (Figure 2M) encoding 4.1B. Its expression was highly significantly downregulated in all prostate cancers compared to normal tissues (Figure 3G). In contrast, a more distant paralog, EPB41L4B encoding EHM2, was strongly upregulated in the 'both' group (Figure 2N), and in cancers overall (Figure 3H). Expression of EPB41L3 and EPB41L4B correlated inversely with each other (Pearson -0.574, p < 0.0001).
A related group of genes identified by interaction analysis encodes components of basement membranes and extracellular matrix. Genes encoding laminin subunits (Figure 2O) and spondins (Figure 2P) have previously been identified as differentially expressed in prostate cancer [14–17, 19]. In addition, we observed FBLN1 encoding the basement membrane component Fibulin-1 to be strongly downregulated in 'both' cancers (Figure 2Q), and in cancers overall (Figure 3I).
In addition to GO groups related to immune response and cell adhesion or whose significance derives from expression differences in genes belong to these larger categories, interaction analysis by the weight algorithm revealed several others (Figure 1C, Table 2), which were not investigated further here. Interestingly, many genes significant in these groups have previously been implicated in prostate cancer. For instance, GO group #43122 contains genes related to NFκB regulation. Constitutive activation of NFκB is widespread in prostate cancer and likely associated with progression (reviewed in ). Increased activity of Gα proteins and PKA (GO #9187) is thought to enhance androgen receptor activity  and influence hedgehog signaling . Expression of monoamine oxidase responsible for the significance of GO group #42135 correlates closely with Gleason score and prostate cancer prognosis . Phosphatic acid phosphatase 2a responsible for the significance of the GO group #8354 was previously identified in a screen for genes differentially expressed in prostate cancer .
For all genes in Fig. 3 the relationship of expression to tumor stage, lymph node involvement, Gleason score, and biochemical recurrence was investigated. After adjustment for multiple testing, significant associations were obtained for two relations. HPN expression was associated with Gleason score, and TLR3 expression was associated with recurrence. As for several other genes, altered expression of TLR3 tended furthermore to be associated with increased tumor stage.
DNA methylation analysis of novel prostate cancer genes
Downregulation of EPB41L3 in other carcinomas is caused by allelic loss at 18p11.3 or promoter hypermethylation [25, 26]. In the present study, eleven cancers showed loss, but only 4 gains at chromosome 18p (Table 1). Bisulfite sequencing revealed EPB41L3 promoter hypermethylation in prostate cancer tissues and cell lines, low methylation in benign prostate tissues and none in leukocytes (Figure 4A). Methylation-specific PCR (Figure 4B) detected EPB41L3 hypermethylation in 79% of the cancer tissues. In prostate carcinoma cell lines, with the exception of PC3, EPB41L3 expression was undetectable even by highly sensitive real-time RT-PCR. Accordingly, the EPB41L3 promoter was strongly methylated in all cell lines, but less so in PC3 (Figure 4A). Combined treatment with a DNA methyltransferase inhibitor, 5-aza-dC, and a histone deacetylase inhibitor, SAHA, induced EPB41L3 expression to detectable levels in Du145 and 22Rv1, while expression in PC3 and LNCaP remained unchanged (data not shown). In the same experiment, 5-aza-dC strongly induced CTCFL and SAHA induced CDKN1A/p21CIP1 expression, as described ().
Treatment with a combination of 5-aza-dC and SAHA increased expression of GADD45A in two prostate cancer cell lines, PC3 and DU145. Induction was slight in 22Rv1, and the inhibitors were inefficacious in LNCaP (Figure 4C). In spite of the response to 5-aza-dC, bisulfite sequencing revealed the GADD45A promoter to be unmethylated in prostate cancer cell lines and prostate cancer tissues with low expression (Figure 4D).
The starting point of the present investigation was our previous observation that chromosome 8 alterations and LINE-1 hypomethylation concurred in many prostate cancers . Such a relationship could be mechanistic in the sense that one alteration causes the other. For instance, chromosome 1 alterations may sometimes be a consequence of hypomethylation of CpG-rich satellites enriched in its juxtacentromeric heterochromatin . Clearly, the fact that the highly significant association was not replicated in the present study makes mechanistic relationships of this kind unlikely. The major difference in the former study was a high proportion of metastatic and recurrent cancers, which exhibited LINE-1 hypomethylation regularly and chromosome 8 alterations often. Therefore, a likely explanation for the discrepancy is that LINE-1 hypomethylation and chromosome 8 alterations both contribute to prostate cancer progression, but are brought about by independent causes.
The fact that they concur in advanced cancers suggests a synergism between the two factors that could be exerted in two ways. Each alteration could target independent biological processes contributing to tumor progression, e.g. one might promote cell proliferation, while the other blocks apoptosis. Alternatively, both factors could target the same process(es). The expression profiling experiment in the present study was designed to distinguish between these alternatives and to identify relevant biological processes. The results clearly indicate that the two alterations act upon the same processes. Very few genes were significantly different between the tumors with or without LINE-1 hypomethylation or between those with or without chromosome 8 alterations, whereas interaction analysis identified several clearly circumscribed biological processes. One prominent process can be roughly categorized as suppression of immune responses, a second one as altered cytoskeletal organization and cell adhesion. Both processes are considered generally important in carcinogenesis, especially during tumor progression. Of note, these processes were identified by the interaction analysis as particularly important in the group of cancers with both alterations. The analysis of individual genes by quantitative real-time PCR revealed that many (e.g. MX1, EPB41L3), but not all (notably BCCIP and TLR3) of the genes significant in the interaction analysis were also overexpressed generally in prostate cancers compared to benign tissues.
More specifically, most individual genes differentially expressed in the immune-response GO groups relate to the activation of innate immune responses following viral infection. A key step in this response is the production of type I interferons by infected cells. These exert antiproliferative effects themselves and initiate adaptive immune responses. Antiproliferative and proapoptotic factors from T-cells, including IFNγ, then help to eliminate infected cells. The same sequence is thought to be involved in immunosurveillance against tumors caused by viruses and other carcinogens . In particular, loss of responsiveness to IFNγ has previously been linked to prostate cancer progression [30, 31]. This result is strengthened by our finding that genes associated with interferon responses are down-regulated in prostate cancer compared to normal tissue overall and even more so in a subgroup of cancers harboring molecular alterations typical of advanced stage tumors. Importantly, however, the genes identified as downregulated in our analysis include several that are involved in the earlier phase mediated by type I interferons, especially the prototypic type I interferon response genes MX1 and MX2. Of note, MX1 has a direct antiproliferative effect , as does GADD45α . IFNB1 itself and IRF3 encoding a key transcription factor for its synthesis appeared downregulated in the group of prostate cancers with both alterations, too. In addition, expression of certain toll-like receptors including TLR2 and TLR3 that recognize 'alien' structures was diminished, at least in this subgroup. Downregulation of TLR3, which detects dsRNA of viral origin and aberrantly methylated endogenous RNAs  is apparent from a previous prostate cancer microarray study, but not was not followed up there . Our findings thus suggest that the immune response to some prostate cancers may fail at at an early step because the cancers are not recognized. The finding that TLR3 down-regulation occurred in a subset of cancers and was associated with recurrence may mean that some cancers are more successful than others in escaping detection by the immune system.
Our findings are particularly intriguing when considered together with current knowledge on hereditary prostate cancer genes. The best characterized HPC1 gene encodes RNaseL, a protein involved in response to dsRNA produced during viral infections or certain other aberrant RNAs . Another hereditary prostate cancer gene candidate is MSR1 at 8p23 likewise involved in innate immune responses . Polymorphisms in TLR genes have also been implicated as predisposition factors [36, 37]. Intriguingly, two genes (BCCIP, HMG20B) appeared in our analysis that interact with BRCA2, another, albeit weaker candidate for a hereditary prostate cancer gene . GADD45A found downregulated in prostate cancers in our study is a well-established BRCA1 target. BRCA1 and IFNγ together act upon the RNaseL pathway .
The preponderance of immune response genes among hereditary prostate cancer genes has been interpreted as indicating a role for an aberrant immune response towards a viral agent in prostate carcinogenesis [1, 35]. Obviously, our data are in line with that interpretation. Of note, the cancers displaying the most pronounced down-regulation of immune response-related genes harbored significant hypomethylation of LINE-1 retrotransposons and chromosome 8 alterations. Diminished methylation may facilitate LINE-1 reexpression . LINE-1 hypomethylation in cancer parallels that of HERV endogenous retroviruses [39, 40] and of ALU (SINE) sequences . Activation of endogenous retroelements occurs in response to various types of genotoxic stress and may also take place during infection by exogenous viruses (reviewed in ). Indeed, processed HERV transcripts have been observed in prostate cancer . HERV proteins are recognized as autoantigens in cancer and autoimmune diseases [41, 43]. ALU transcripts influence the dsRNA-dependent protein kinase PKR, which regulates RNaseL activation, and TLRs . Conceivably, therefore, what appears as an aberrant response to an exogenous viral agent in prostate carcinogenesis might reflect to some extent the activation of endogenous retroelements. Interestingly, the recent completion of the chromosome 8 sequence has highlighted an unusual concentration of genes involved in the regulation of the immune response on this particular chromosome . Therefore, a speculative interpretation of our findings is that increased expression of endogenous retroelements in prostate cancer cells with hypomethylated genomes would contribute to the activation of immune responses, but is tolerated as a consequence of chromosome 8 alterations. We suggest this as a working hypothesis for further studies.
Similarly, our study implies that altered expression patterns of cytoskeletal and extracellular matrix proteins in prostate carcinoma are linked to alterations of chromosome 8 and LINE-1 hypomethylation. A novel finding in this context is the involvement of 4.1 proteins in prostate cancer. These proteins contain a phosphatidyl inositol phosphate-binding domain and connect a variety of transmembrane proteins to the actin cytoskeleton, thereby organizing cell polarity and motility. The two members of the family investigated here more closely are already implicated in other cancers. EPB41L3 was proposed as a potential tumor suppressor in lung cancer  and was shown to be down-regulated by promoter hypermethylation and allelic loss in renal cell carcinoma . We report here that hypermethylation of EPB41L3 is also prevalent in prostate cancer. The 18p losses detected by CGH in 11/50 cases, could contribute to down-regulation. EHM2 encoded by EPB41L4B was identified as overexpressed in metastatic melanoma cells . We show here that its overexpression is common in prostate cancer. This is the first explicit study on this gene, but our findings are in line with significant differences in previous microarray data [14, 15]. Taken together, the results suggest a shift in the pattern of 4.1 proteins associated with prostate cancer progression whose biological and clinical implications deserve further investigation.
Altered expression of Fibulin-1, which interacts with laminins reported to be downregulated in several previous studies [14–17], was previously reported only in gynaecological cancers [47, 48]. In these cancers, FBLN1 becomes upregulated during progression. In contrast and therefore surprisingly, in prostate cancer FBLN1 appears to become generally downregulated. The decreases in laminin and Fibulin-1 expression may be related to the dissolution of the basement membrane in prostate cancer tissues. Like the 4.1 proteins, Fibulin-1 influences cell motility and polarity. As overexpression and abnormal localization in breast cancer can lead to recognition of Fibulin-1 as an autoantigen , its decreased expression in prostate cancers may also lower their immunogenicity.
Promoter hypermethylation underlies GADD45A downregulation in breast cancer . In prostate cancer altered methylation at a more distant, but unfortunately unspecified site in the gene has been reported . We found the actual promoter unmethylated in prostate cancer tissues and cell lines. Induction of GADD45A by SAHA treatment in prostate cancer cell lines could mean that downregulation is associated with histone deacetylation and altered chromatin structure. The additive effect of 5-aza-dC treatment might reflect the methylation at the distant site .
LINE-1 hypomethylation and chromosome 8 alterations are commonly associated in highly advanced prostate cancers. The present study suggests that both alterations do not cause each other. Rather, they appear to act as converging and even synergistic factors contributing to prostate cancer progression. Interaction analysis identified suppression of innate immune responses and cytoskeletal and extracellular matrix changes as common targets. This implies that candidate tumor genes on chromosome 8 as well as DNA hypomethylation should be considered for their influence on these processes in prostate cancer. Intriguingly, many hereditary prostate cancers genes are also involved in regulation of innate immune responses. The data stress the importance of altered cell adhesion and cytoskeletal organization in prostate cancer and specifically implicate changes in 4.1 protein expression in this process. The 4.1 protein encoding genes and several others newly identified here ought to be investigated for their usefulness in prostate cancer detection and classification through expression or methylation analyses.
From a previous study of DNA methylation alterations in a series of 113 prostate carcinomas , a subset of 50 specimens (Table 1) was selected according to DNA availability, RNA quality and complete follow-up (median period 62 months). Cases with distant metastases at the time of surgery were excluded. Twelve cases had lymph node metastases. Benign tissues were taken from distant locations of cancer-carrying prostates as described . The study was approved by the HHU medical faculty ethics committee.
The prostate carcinoma cell lines LNCaP, 22RV1, PC3, and DU145 were cultured and treated with epigenetic inhibitors, 5-aza-2'-deoxycytidine (5-aza-dC; Sigma, Taufkirchen, Germany), suberoylanilide hydroxamic acid (SAHA, Biomol, Hamburg, Germany) as described . 5-aza-dC was supplied at 2 μM every 24 h for 3 days and suberoylanilide hydroxamic acid at 2 μM for 2 days.
DNA and RNA extraction
Comparative genomic hybridization
CGH was performed as described , using the same aliquots as for DNA methylation analysis.
RNA microarray analysis
High-quality total RNA (5 μg each) was converted to double-stranded cDNA before in vitro transcription with biotinylated dNTPs using the GeneChip® Expression 3' Amplification One-Cycle Target Labeling Kit (Affymetrix, Santa Clara, USA). The resulting biotinylated cRNA was fragmented and hybridized to HG-U133A (Affymetrix) microarrays according to the manufacturer's specifications. After hybridization, the microarrays were washed, stained with streptavidin/phycoerythrin conjugate and biotinylated antibody, and scanned according to the manufacturer's recommendations using an Affymetrix GeneChip® Scanner 3000. Images were processed using GeneChip® Operating Software (GCOS, Version 1.3, Affymetrix) and total intensity normalization was applied by normalizing all arrays to an average signal level of 500 counts.
Bioinformatic analysis of microarray data
For every gene separately, a multivariate linear model for predicting logarithmic gene expression was estimated. Two main effects for alteration of chromosome 8 and for hypomethylation of LINE-1 retrotransposons and an interaction effect for joint occurrence of these two factors were included in the model. Significance values for every gene and for the three effects, respectively, were obtained by applying the function lm for fitting linear models as implemented in the R programming language R .
GO group scoring
Two methods for scoring the significance of enrichment of a list of differentially expressed genes with genes belonging to a Gene Ontology (GO) group were applied, namely the classic and the weight algorithm . In both algorithms, genes are first ranked according to a score that quantifies the amount of differential expression. Corresponding p-values are then adjusted for multiple testing according to the false discovery rate (fdr) method. All genes with adjusted p-values p < 0.05 are included in the list of significant genes. The significance of a GO term is then obtained by comparing the observed number of genes in the respective GO group that are members of this list with the expected number calculated in a model that assumes independence between the GO group and the list of significant genes. In the classic algorithm, this comparison is based on Fisher's exact test which is applied for each GO term independently. Due to the graph topology of the Gene Ontology, pairs of terms with parent-child relationships automatically receive similar p-values and thus appear simultaneously among the most significant GO terms. In the weight algorithm, genes annotated to a GO term receive weights based on the significance of neighboring GO terms . This approach identifies local dependencies between related GO terms and highlights those terms that receive a more significant score than all their neighbors. The algorithms for GO group scoring were implemented in the R programming language . The results were obtained using R version 2.3.0 and the libraries provided by the Bioconductor project, version 1.8 (released on April 27th, 2006).
Quantitative real-time RT-PCR was performed using RNA from the identical tissue specimens on an ABI 7900 instrument using commercially available primers and probes specific for the respective mRNAs (Applied Biosystems, Weiterstadt, Germany). Each run was standardized using a dilution series of a strongly expressing cell line or normal tissue. Experimental variation for each sample was below 10%. All results were expressed relative to TBP used as a reference gene.
DNA methylation analysis
Bisulfite sequencing DNA methylation analysis was performed as described [9, 11] using the following novel primers for EPB41L3 DAL1fwd 5'-GTAATAGGGGG(T/C)GGGGGGAATAG-3', DAL1rev 5'-AACCCCCTC(A/G)CAATCCCCCACTC-3' for GADD45A: GADD45Afwd 5'-TTAGTGGTTGGTAGGTAGTGGTT-3', GADD45Arev 5'-CCTCCAAAATCATATTACAAACTAC-3'. Methylation-specific PCR was performed as described  using the novel primer pairs for EPB41L3 DAL1US 5'-TTTGTGTATTGTTGTTGAGGAGTG-3' and DAL1UAS 5'-CACAATCCCCCACTCCAAAAAACA-3' to detect unmethylated sequences or DAL1MS 5'-TTGCGTATCGTCGTCGTCGAGGACG-3' and DAL1MAS 5'-CGCAATCCCCCACTCCGAAAAACG-3' to detect methylated sequences at 61°C and 64°C annealing temperature, respectively.
DeMarzo AM, Nelson WG, Isaacs WB, Epstein JI: Pathological and molecular aspects of prostate cancer. Lancet. 2003, 361: 955-964. 10.1016/S0140-6736(03)12779-1
Dong JT: Prevalent mutations in prostate cancer. J Cell Biochem. 2006, 97: 433-447. 10.1002/jcb.20696
Schulz WA, Hatina J: Epigenetics of prostate cancer: beyond DNA methylation. J Cell Mol Med. 2006, 10: 100-125.
Kindich R, Florl AR, Kamradt J, Lehmann J, Müller M, Wullich B, Schulz WA: Relationship of NKX3.1 and MYC gene copy number ratio and DNA hypomethylation to prostate carcinoma stage. Eur Urol. 2006, 49: 169-175. 10.1016/j.eururo.2005.09.012
Macoska JA, Beheshti B, Rhim JS, Hukku B, Lehr J, Pienta KJ, Squire JA: Genetic characterization of immortalized human prostate epithelial cell cultures. Evidence for structural rearrangements of chromosome 8 and i(8q) chromosome formation in primary tumor-derived cells. Cancer Genet Cytogenet. 2000, 120: 50-57. 10.1016/S0165-4608(99)00248-4
Ehrlich M: DNA methylation in cancer: too much, but also too little. Oncogene. 2002, 21: 5400-5413. 10.1038/sj.onc.1205651
Hoffmann MJ, Schulz WA: Causes and consequences of DNA hypomethylation in human cancer. Biochem Cell Biol. 2005, 83: 296-321. 10.1139/o05-036
Schulz WA, Elo JP, Florl AR, Pennanen S, Santourlidis S, Engers R, Buchardt M, Seifert HH, Visakorpi T: Genomewide DNA hypomethylation is associated with alterations on chromosome 8 in prostate carcinoma. Genes Chromosomes Cancer. 2002, 35: 58-65. 10.1002/gcc.10092
Florl AR, Steinhoff C, Müller M, Seifert HH, Hader C, Engers R, Ackermann R, Schulz WA: Coordinate hypermethylation at specific genes in prostate carcinoma precedes LINE-1 hypomethylation. Br J Cancer. 2004, 91: 985-994.
Yegnasubramanian S, Kowalski J, Gonzalgo ML, Zahurak M, Piantadosi S, Walsh PC, Bova GS, De Marzo AM, Isaacs WB, Nelson WG: Hypermethylation of CpG islands in primary and metastatic human prostate cancer. Cancer Res. 2004, 64: 1975-1986. 10.1158/0008-5472.CAN-03-3972
Hoffmann MJ, Müller M, Engers R, Schulz WA: Epigenetic control of CTCFL/BORIS and OCT4 in urogenital malignancies. Biochem Pharmacol. 2006, 72: 1577-1588. 10.1016/j.bcp.2006.06.020
Jung V, Kindich R, Kamradt J, Jung M, Müller M, Schulz WA, Engers R, Unteregger G, Stöckle M, Zimmermann R, Wullich B: Genomic and expression analysis of the 3q25-q26 amplification unit reveals TLOC1/SEC62 as a probable target gene in prostate cancer. Mol Cancer Res. 2006, 4: 169-176. 10.1158/1541-7786.MCR-05-0165
Alexa A, Rahnenführer J, Lengauer T: Improved scoring of functional groups from gene expression data by decorrelating GO graph structure. Bioinformatics. 2006, 22: 1600-1607. 10.1093/bioinformatics/btl140
Dhanasekaran SM, Barrette TR, Ghosh D, Shah R, Varambally S, Kurachi K, Pienta KJ, Rubin MA, Chinnaiyan AM: Delineation of prognostic biomarkers in prostate cancer. Nature. 2001, 412: 822-826. 10.1038/35090585
Luo J, Duggan DJ, Chen Y, Sauvageot J, Ewing CM, Bittner ML, Trent JM, Isaacs WB: Human prostate cancer and benign prostatic hyperplasia: molecular dissection by gene expression profiling. Cancer Res. 2001, 61: 4683-4688.
Welsh JB, Sapinoso LM, Su AI, Kern SG, Wang-Rodriguez J, Moskaluk CA, Frierson HF, Hampton GM: Analysis of gene expression identifies candidate markers and pharmacological targets in prostate cancer. Cancer Res. 2001, 61: 5974-5978.
Singh D, Febbo PG, Ross K, Jackson DG, Manola J, Ladd C, Tamayo P, Renshaw AA, D'Amico AV, Richie JP, Lander ES, Loda M, Kantoff PW, Golub TR, Sellers WR: Gene expression correlates of clinical prostate cancer behavior. Cancer Cell. 2002, 1: 203-209. 10.1016/S1535-6108(02)00030-2
Wang C, McCarty IM, Balazs L, Li Y, Steiner MS: Immunohistological detection of BRAF25 in human prostate tumor and cancer specimens. Biochem Biophys Res Commun. 2002, 295: 136-141. 10.1016/S0006-291X(02)00625-3
Edwards S, Campbell C, Flohr P, Shipley J, Giddings I, Te-Poele R, Dodson A, Foster C, Clark J, Jhavar S, Kovacs G, Cooper CS: Expression analysis onto microarrays of randomly selected cDNA clones highlights HOXB13 as a marker of human prostate cancer. Br J Cancer. 2005, 92: 376-381. 10.1038/sj.bjc.6602261
Lessard L, Mes-Masson AM, Saad F: NFêB: un nouveau marqueur kappable de prédire l'évolution du cancer de la prostate. Bull Cancer. 2006, 93: 891-899.
Kasbohm EA, Guo R, Yowell CW, Bagchi G, Kelly P, Arora P, Casey PJ, Daaka Y: Androgen receptor activation by G(s) signaling in prostate cancer cells. J Biol Chem. 2005, 280: 11583-11589. 10.1074/jbc.M414423200
Jia J, Tong C, Wang B, Luo L, Jiang J: Hedgehog signalling activity of Smoothened requires phosphorylation by protein kinase A and casein kinase I. Nature. 2004, 432: 1045-1050. 10.1038/nature03179
True L, Coleman I, Hawley S, Huang CY, Gifford D, Coleman R, Beer TM, Gelmann E, Datta M, Mostaghel E, Knudsen B, Lange P, Vessela R, Lin D, Hood L, Nelson PS: A molecular corelate to the Gleason grading system for prostate adenocarcinoma. Proc Natl Acad Sci USA. 2006, 103: 10991-10996. 10.1073/pnas.0603678103
Porkka KP, Visakorpi T: Detection of differentially expressed genes in prostate cancer by combining suppression subtractive hybridization and cDNA library array. J Pathol. 2001, 193: 73-79. 10.1002/1096-9896(2000)9999:9999<::AID-PATH751>3.0.CO;2-Y
Tran YK, Bogler O, Gorse KM, Wieland I, Green MR, Newsham IF: A novel member of the NF2/ERM/4.1 superfamily with growth suppressing properties in lung cancer. Cancer Res. 1999, 59: 35-43.
Yamada D, Kikuchi S, Williams YN, Sakurai-Yageta M, Masuda M, Maruyama T, Tomita K, Gutmann DH, Kakizoe T, Kitamura T, Kanai Y, Murakami Y: Promoter hypermethylation of the potential tumor suppressor DAL-1/4.1B gene in renal clear cell carcinoma. Int J Cancer. 2006, 118: 916-923. 10.1002/ijc.21450
Wang W, Huper G, Guo Y, Murphy SK, Olson JA, Marks JR: Analysis of methylation-sensitive transcriptome identifies GADD45a as a frequently methylated gene in breast cancer. Oncogene. 2005, 24: 2705-2714. 10.1038/sj.onc.1208464
Ji W, Hernandez R, Zhang XY, Qu GZ, Frady A, Varela M, Ehrlich M: DNA demethylation and pericentromeric rearrangements of chromosome 1. Mutat Res. 1997, 379: 33-41.
Dunn GP, Old LJ, Schreiber RD: The immunobiology of cancer immunosurveillance and immunoediting. Immunity. 2004, 21: 137-148. 10.1016/j.immuni.2004.07.017
Shou J, Soriano R, Hayward SW, Cunha GR, Williams PM, Gao WQ: Expression profiling of a human cell line model of prostatic cancer reveals a direct involvement of interferon signaling in prostate tumor progression. Proc Natl Acad Sci USA. 2002, 99: 2830-2835. 10.1073/pnas.052705299
Nagano K, Masters JR, Akpan A, Yang A, Corless S, Wood C, Hastie C, Zvelebil M, Cramer R, Naaby-Hansen S: Differential protein synthesis and expression levels in normal and neoplastic human prostate cells and their regulation by type I and II interferons. Oncogene. 2004, 23: 1693-1703. 10.1038/sj.onc.1207297
Numajiri A, Mibayashi M, Nagata K: Stimulus-dependent and domain-dependent cell death acceleration by an IFN-inducible protein, human MxA. J Interferon Cytokine Res. 2006, 26: 214-219. 10.1089/jir.2006.26.214
Mak SK, Kultz D: Gadd45 proteins induce G2/M arrest and modulate apoptosis in kidney cells exposed to hyperosmotic stress. J Biol Chem. 2004, 279: 39075-39084. 10.1074/jbc.M406643200
Kariko K, Buckstein M, Ni H, Weissman D: Suppression of RNA recognition by Toll-like receptors: the impact of nucleoside modification and the evolutionary origin of RNA. Immunity. 2005, 23: 165-175. 10.1016/j.immuni.2005.06.008
Silverman RH: Implications for RNase L in prostate cancer biology. Biochemistry. 2003, 42: 1805-1812. 10.1021/bi027147i
Zheng SL, Augustsson-Balter K, Chang B, Hedelin M, Li L, Adami HO, Bensen J, Li G, Johnasson JE, Turner AR, Adams TS, Meyers DA, Isaacs WB, Xu J, Gronberg H: Sequence variants of toll-like receptor 4 are associated with prostate cancer risk: results from the CAncer Prostate in Sweden Study. Cancer Res. 2004, 64: 2918-2922. 10.1158/0008-5472.CAN-03-3280
Sun J, Wiklund F, Hsu FC, Balter K, Zheng SL, Johansson JE, Chang B, Liu W, Li T, Turner AR, Li L, Li G, Adami HO, Isaacs WB, Xu J, Gronberg H: Interactions of sequence variants in interleukin-1 receptor-associated kinase4 and the toll-like receptor 6-1-10 gene cluster increase prostate cancer risk. Cancer Epidemiol Biomarkers Prev. 2006, 15: 480-485. 10.1158/1055-9965.EPI-05-0645
Mullan PB, Hosey AM, Buckley NE, Quinn JE, Kennedy RD, Johnston PG, Harkin DP: The 2, 5 oligoadenylate synthetase/RNaseL pathway is a novel effector of BRCA1 and interferon γ induced apoptosis. Oncogene. 2005, 24: 5492-5501. 10.1038/sj.onc.1208698
Florl AR, Löwer R, Schmitz-Dräger BJ, Schulz WA: DNA methylation and expression of LINE-1 and HERV-K provirus sequences in urothelial and renal cell carcinomas. Br J Cancer. 1999, 80: 1312-1321. 10.1038/sj.bjc.6690524
Weisenberger DJ, Campan M, Long TI, Kim M, Woods C, Fiala E, Ehrlich M, Laird PW: Analysis of repetitive element DNA methylation by MethyLight. Nucleic Acids Res. 2005, 33: 6823-6836. 10.1093/nar/gki987
Schulz WA, Steinhoff C, Florl AR: Methylation of endogenous human retroelements in health and disease. Curr Topics Microbiol Immunol. 2006, 310: 211-250.
Wang-Johanning F, Frost AR, Jian B, Azerou R, Lu DW, Chen DT, Johanning GL: Detecting the expression of human endogenous retrovirus E envelope transcripts in human prostate adenocarcinoma. Cancer. 2003, 98: 187-197. 10.1002/cncr.11451
Bannert N, Kurth R: Retroelements and the human genome: new perspectives on an old relation. Proc Natl Acad Sci USA. 2004, 101 (Suppl 2): 14572-14579. 10.1073/pnas.0404838101
Li TH, Schmid CW: Differential stress induction of individual Alu loci: implications for transcription and retrotransposition. Gene. 2001, 276: 135-141. 10.1016/S0378-1119(01)00637-0
Nusbaum C, Mikkelsen TS, Zody MC, Asakawa S, Taudien S, Garber M, Kodira CD, Schueler MG, Shimizu A, Whittaker CA, Chang JL, Cuomo CA, Dewar K, FitzGerald MG, Yang X, Allen NR, Anderson S, Asakawa T, Blechschmidt K, Bloom T, Borowsky ML, Butler J, Cook A, Corum B, DeArellano K, DeCaprio D, Dooley KT, Dorris L, Engels R, Glockner G, Hafez N, Hagopian DS, Hall JL, Ishikawa SK, Jaffe DB, Kamat A, Kudoh J, Lehmann R, Lokitsang T, Macdonald P, Major JE, Matthews CD, Mauceli E, Menzel U, Mihalev AH, Minoshima S, Murayama Y, Naylor JW, Nicol R, Nguyen C, O'Leary SB, O'Neill K, Parker SC, Polley A, Raymond CK, Reichwald K, Rodriguez J, Sasaki T, Schilhabel M, Siddiqui R, Smith CL, Sneddon TP, Talamas JA, Tenzin P, Topham K, Venkataraman V, Wen G, Yamazaki S, Young SK, Zeng Q, Zimmer AR, Rosenthal A, Birren BW, Platzer M, Shimizu N, Lander ES: DNA sequence and analysis of human chromosome 8. Nature. 2006, 439: 331-335. 10.1038/nature04406
Shimizu K, Nagamachi Y, Tani M, Kimura K, Shiroishi T, Wakana S, Yokota J: Molecular cloning of a novel NF2/ERM/4.1 superfamily gene, ehm2, that is expressed in high-metastatic K1735 murine melanoma cells. Genomics. 2000, 65: 113-120. 10.1006/geno.2000.6154
Moll F, Katsaros D, Lazennec G, Hellio N, Roger P, Giacalone PL, Chalbos D, Maudelonde T, Rochefort H, Pujol P: Estrogen induction and overexpression of fibulin-1C mRNA in ovarian cancer cells. Oncogene. 2002, 21: 1097-1107. 10.1038/sj.onc.1205171
Pupa SM, Argraves WS, Forti S, Casalini P, Berno V, Agresti R, Aiello P, Invernizzi A, Baldassari P, Twal WO, Mortarini R, Anichini A, Menard S: Immunological and pathobiological roles of fibulin-1 in breast cancer. Oncogene. 2004, 23: 2153-2160. 10.1038/sj.onc.1207323
Wang Y, Yu Q, Cho AH, Rondeau G, Welsh J, Adamson E, Mercola D, McClelland M: Survey of differentially methylated promoters in prostate cancer cell lines. Neoplasia. 2005, 7: 748-760. 10.1593/neo.05289
The R Project for Statistical Computing., [http://www.r-project.org]
We are grateful to Dr. Gernot Röder (Institute of Oncological Chemistry, Heinrich Heine University) for carrying out the microarray experiment and to Dr. Jiri Hatina (supported by the German Academic Exchange Service) for helpful comments.
Financial support was provided by the Deutsche Krebshilfe (70–3193 Schu I, to WAS and BW) and the BMBF (grant No. 01GR0453 to AA and JR). The work at the Max-Planck-Institute for Informatics was performed in the context of the BioSapiens Network of Excellence (EU contract no. LSHG-CT-2003-503265). The funding agencies did not exert influence on decisions in the course of the experimental study, the writing of the manuscript, or the decision to submit the manuscript for publication.
The author(s) declare that they have no competing interests.
WAS, BW, and JR conceived the study. WAS and JR drafted the manuscript. VJ performed the CGH analysis evaluating it together with BW. CH, MJH, MY, SF, AW, and ARF carried out and evaluated the expression and methylation analyses with guidance by WAS. AA, TL, and JR performed the bioinformatic analyses. MM contributed the clinical data and according statistical evaluations. RE controlled the tissue sampling and histopathological data. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table of genes significant in the 20 most significant GO groups. Genes annotated to the 20 most significant GO groups that are significantly differentially expressed (p < 0.01) between prostate cancer groups according to interaction analysis (DOC 386 KB)
Additional file 2: Box plot representation of microarray analysis results for genes encoding markers of innate immune-response cell types. In each graph, the expression values (log-scale) for the respective genes are depicted for the cancers with neither hypomethylation of LINE-1 retrotransposons nor alteration of chromosome 8 (null), hypomethylation only (hypo), alteration of chromosome 8 only (chrom8) and both alterations (chrom8 AND hypo). The p-values refer to the result of the interaction analysis (see Methods for details). KLRD1 and ITGA2 gene products are characteristic of NK cells, CPA3 of mast cells, CD163 is relatively characteristic of macrophages (the more characteristic marker CD68 was not represented on the HGU133A microarray), CD83 and ADAM19 are typical for dendritic cells, and ITGAM encodes a marker of granulocytes. (PDF 93 KB)
About this article
Cite this article
Schulz, W.A., Alexa, A., Jung, V. et al. Factor interaction analysis for chromosome 8 and DNA methylation alterations highlights innate immune response suppression and cytoskeletal changes in prostate cancer. Mol Cancer 6, 14 (2007). https://doi.org/10.1186/1476-4598-6-14
- Prostate Cancer
- Gene Ontology
- Interaction Analysis
- Promoter Hypermethylation