The oxidative demethylase ALKBH3 marks hyperactive gene promoters in human cancer cells
Genome Medicine volume 7, Article number: 66 (2015)
The oxidative DNA demethylase ALKBH3 targets single-stranded DNA (ssDNA) in order to perform DNA alkylation damage repair. ALKBH3 becomes upregulated during tumorigenesis and is necessary for proliferation. However, the underlying molecular mechanism remains to be understood.
To further elucidate the function of ALKBH3 in cancer, we performed ChIP-seq to investigate the genomic binding pattern of endogenous ALKBH3 in PC3 prostate cancer cells coupled with microarray experiments to examine the expression effects of ALKBH3 depletion.
We demonstrate that ALKBH3 binds to transcription associated locations, such as places of promoter-proximal paused RNA polymerase II and enhancers. Strikingly, ALKBH3 strongly binds to the transcription initiation sites of a small number of highly active gene promoters. These promoters are characterized by high levels of transcriptional regulators, including transcription factors, the Mediator complex, cohesin, histone modifiers, and active histone marks. Gene expression analysis showed that ALKBH3 does not directly influence the transcription of its target genes, but its depletion induces an upregulation of ALKBH3 non-bound inflammatory genes.
The genomic binding pattern of ALKBH3 revealed a putative novel hyperactive promoter type. Further, we propose that ALKBH3 is an intrinsic DNA repair protein that suppresses transcription associated DNA damage at highly expressed genes and thereby plays a role to maintain genomic integrity in ALKBH3-overexpressing cancer cells. These results raise the possibility that ALKBH3 may be a potential target for inhibiting cancer progression.
Genomic DNA is continuously subjected to various harmful insults, such as UV light, ionizing radiation, or nucleic-acid modifying compounds, resulting in thousands of DNA alterations in each cell every day . Such lesions can lead to DNA damage, which in turn favors mutagenesis, carcinogenesis, inflammation, and aging [2–5]. Accordingly, cells have multiple mechanisms to reverse damaging DNA modifications. In particular, DNA alkylation, a process of methylating specific nucleic acids, often requires repair to maintain genomic integrity. Alkylating agents are found ubiquitously in the environment, but DNA can also be alkylated as a natural by-product of cellular metabolism [6, 7]. For example, the universal methyl donor S-adenosylmethionine non-enzymatically methylates DNA [8, 9]. Alkylating agents preferentially attack single-stranded DNA (ssDNA) in the genome due to its higher accessibility [10–13], and some DNA modifications such as 1-methyladenine (1-meA) and 3-methylcytosine (3-meC) are primarily generated in ssDNA, because these positions are shielded in double-stranded DNA (dsDNA) .
DNA alkylation can be removed by base-excision repair (BER), direct reversal by methylguanine methyltransferase (MGMT), and dealkylation via the AlkB family [6, 7]. The AlkB enzymes belong to a large family of non-heme Fe(II) and 2-oxoglutarate-dependent dioxygenases, which catalyze numerous biological reactions, such as proline hydroxylation and histone demethylation . AlkB was original discovered in E. coli, where it demethylates 1-methyladenine (1meA) and 3-methylcytosine (3meC) by oxidation of the N-linked methyl moiety. This reaction creates an unstable methyl-iminium intermediate that spontaneously hydrolyzes into formaldehyde and the non-alkylated base [15–17]. In mammalians, at least nine ALKB family members are known (ALKBH1-8 and FTO). DNA damage dealkylation reactions are mainly catalyzed by ALKBH2 and ALKBH3 . Notably, ALKBH2 preferentially demethylates dsDNA while ALKBH3 demethylates ssDNA and RNA substrates, modified by 3-meC or 1-meA [16, 18, 19]. Since these modifications are predominantly generated in ssDNA and RNA, it has been proposed that the ALKBH3 repair function could be linked to transcription [6, 20]. Incomplete removal of DNA alkylation leads to DNA damage, resulting in cell cycle arrest, inflammation and apoptosis [3–5, 21]. Induction of DNA alkylation by chemotherapeutic agents is a common strategy in cancer treatment to prevent cancer cells from dividing and proliferating . Enzymes that facilitate DNA alkylation damage repair, such as ALKBH2 and ALKBH3, can contribute to resistance to this treatment and insights into their molecular function could provide the basis for developing more efficient cancer therapies .
Recently, we described the cooperativity of ALKBH3 and the ASCC3 DNA helicase complex to promote DNA alkylation damage repair in various cancer cells. ALKBH3 knockdown causes elevated levels of 3-meC accompanied by increased DNA damage response (DDR) and reduced cell proliferation . However, the mechanisms of in vivo genomic targeting of ALKBH3 are not yet fully understood.
Herein, using chromatin immunoprecipitation experiments followed by massively parallel sequencing analysis (ChIP-seq) we find that in PC3 prostate cancer cells ALKBH3 binding is enriched at transcription associated genomic loci, where ssDNA is accessible. Specifically, we find ALKBH3 bound at active gene promoters, enhancers, and regions with putative quadruplex DNA. Unexpectedly, ALKBH3 binds strongly to the initiation sites of some particularly highly expressed gene promoters. Interestingly, these promoters are bound by an unusually large number of transcriptional regulators, indicating a highly regulated ‘hyperactive’ promoter class. However, we find that loss of ALKBH3 does not directly affect expression of ALKBH3 occupied genes, suggesting a transcription unrelated function of ALKBH3. Instead, upon ALKBH3 knockdown we observe an increased expression of genes involved in inflammatory pathways, which could be a downstream effect of elevated DNA damage after ALKBH3 depletion [24, 25]. The genomic localization of ALKBH3 at transcription-related loci raises the possibility that ALKBH3 could have a role in suppressing transcription-associated DNA damage to preserve the genomic integrity.
Cell culture and viral transduction
U2OS, 293 T, NCI-H23, and PC3 cells were obtained from the American Type Culture Collection (ATCC) and maintained as previously described . ShRNAs constructs, preparation of viruses and cell transduction have been described previously . Cells infected with lentiviral shRNAs were selected after infection with puromycin (1 μg/mL) for at least 48 h.
Rabbit anti-ALKBH3 antibodies were obtained from Millipore (Catalog #09-882).
U2OS and PC3 cells were used for IF and were maintained as described previously . The cells were infected with the indicated lentiviral shRNAs and plated onto coverslips. Cells were then fixed with PBS (ph 7.4) containing 3.2 % paraform for 20 min, washed extensively with IF wash buffer (1X PBS containing 0.5 % NP-40 and 0.02 % NaN3), then incubated with blocking buffer (IF wash buffer with 10 % fetal bovine serum), and finally stained with anti-rabbit-ALKBH3 (Millipore) diluted in blocking buffer. Secondary antibodies (goat anti-rabbit Alexa Fluor 488) were from Millipore.
RNAi experiments, microarray analysis, and qRT-PCR
Lentiviral shRNA against ALKBH3, ALKBH2, and GFP as control were used as described before . Briefly, lentiviral production was carried out in 293 T cells and target cells were infected for 48 h followed by selection with puromycin for 48 h or 96 h. RNA extracted from ALKBH3 knockdown and control cells were sent to Transcriptome Analysis Laboratory (TAL, University Medical Center, Göttingen) for expression analysis. Microarrays were done using the ‘Low RNA Input linear Amplification Kit Plus, One Color’ protocol (Agilent Technologies, Cat. No.: 5188–5339) and the Agilent RNA Spike-In Kit for One color (Agilent Technologies, Cat. No.: 5188–5282) following the manufacturer’s standard protocol. Global gene expression analysis was applied using the Human Gene Expression 4x44K v2 Microarray Kit (Agilent Technologies, Cat. No.: G4845A). 200 ng of total RNA from each sample from PC3 cells were used as a starting material to prepare cDNA. The hybridizations were performed for 17 h at 10 rpm and 65 °C in the Hybridization Oven (Agilent). Washing and staining of the arrays were done according to the manufacturer’s recommendation. Cy3 intensities were detected by one-color scanning using an Agilent DNA microarray scanner (G2505B). Intensity data were extracted using Agilent’s Feature Extraction (FE) software (version 10.5.3.1) including a quality control based on internal controls using Agilent’s protocol GE1_107_Sep09. A subset of genes from microarray results was verified by qRT-PCR (Fig. 3d). For RT-PCR, cells were infected with lentivirus shRNA against ALKBH3 or GFP and selected with puromycin. qRT-PCR was performed using the One-Step Sybr No Rox Kit from Bioline on a CFX384 Real-Time System (Bio-Rad). Primers are shown in Table S1 (Additional file 1).
ChIP and ChIP-Seq
Chromatin was prepared from PC3 cells as previously described , except that it was sonicated to 200 base pairs (bp). Chromatin was incubated with total of 6 μg specific antibodies against ALKBH3 (Millipore) or IgG (Abcam) overnight at 4 °C and then mixed with 50 % slurry protein A beads (Millipore) for 1 h at 4 °C. The beads were washed extensively, de-cross-linked at 65 °C, and treated with RNase A and proteinase K. Samples were next subjected to phenol-chloroform extraction and precipitated with ice-cold ethanol. Libraries were constructed of 50 ng of ChIP’d DNA following Illumina’s® protocol and sequenced using an Illumina Genome Analyzer before being further analyzed bioinformatically. Analysis of DNA via qRT-PCR was performed as described above, with gene-specific primers (Additional file 1: Table S2).
Microarrays were normalized and analyzed using the limma package for Bioconductor [27, 28]. ChIP-Seq data (Additional file 1: Table S3) were mapped to human genome hg19 using bowtie version 1.0 , allowing one mismatch (n = 1) and maximal three possible alignments (m = 3). All subsequent analyses of ChIP-Seq data were performed using the Cistrome platform [30, 31] (Galaxy Code 2014.5.5). If possible, data were also directly uploaded into Cistrome from the GEO database. For Promoter definition, RefSeq genes were downloaded from the UCSC Genome Browser. After removal of duplicates with identical transcription start site, 31,296 promoters, including genes with alternative transcription start sites, were used for analysis. Promoters with weakly bound ALKBH3 were identified using the k-means clustering function within the heatmap feature in Cistrome. Enhancers were defined as overlapping peaks (called by MACS with P value 1e-05) of H3K4me1 and H3K27ac (from LNCaP cells), which do not overlap with promoters (−1,000/+1,000). ETS transcription factor bindings sites were called using MACS with a cutoff P value of 1e-05. Promoters overlapping with an ETS binding site were considered as ETS factor bound promoters. Predicted G4 DNA sites were downloaded from  and converted to hg19. Only sites that are not at promoter or enhancers sites were used for analysis. CpG islands and Transcription factor binding sites (TfbsClusteredV3) were downloaded from the UCSC table browser. A promoter transcription factor binding event has been defined as an overlap of a clustered transcription factor binding site with a promoter site. For TATA-box analysis, conserved transcription factor binding sites (tfbsConsSites) were downloaded from UCSC browser. Conserved TBP bindings sites (V$TBP_01) overlapping with promoter sites were considered as TATA-Box. The counting of the ChIP-Seq tags at each promoter was done using a custom R script for Bioconductor.
ChIP-Seq and Microarray data are available at the GEO repository with the accession numbers GSE57568 and GSE57591.
The significance of the data was either calculated by Cistrome, via unpaired Student’s t-tests, or has been evaluated using hypergeometric probability tests.
ALKBH3 occupies ubiquitously expressed promoters in PC3 cells
Our previous work showed that ALKBH3 depletion in PC3 prostate cancer cells increases global 3-meC levels and induces H2A.X phosphorylation as well as 53BP1 foci formation , demonstrating the occurrence of systemic DNA damage. To gain further genome-wide insights into the DNA repair function of ALKBH3, we performed ChIP-seq experiments in PC3 cells using specific anti-ALKBH3 antibodies. First, we performed Model-based Analysis for ChIP-Seq (MACS) (cutoff P value: 1e-04) within the Cistrome platform [30, 31] to identify 423 high confidence ALKBH3 binding sites. ALKBH3 predominantly occupies promoters and regions downstream of the transcription start site (TSS) (5′-UTR), but is excluded from introns and other genomic regions (Fig. 1a). Out of the 423 peaks, 354 are present at known promoters. Further analysis using Cistrome revealed that ALKBH3 is weakly bound but enriched at more than 4,000 additional promoters. We segregated ALKBH3 bound promoters into two groups (I and II) depending on the ALKBH3 binding strength (Figs. 1b, 4a). At these promoters ALKBH3 localizes around the TSS (Fig. 1c). We confirmed the specificity of the antibody via immunofluorescence and ChIP experiments in ALKBH3 knockdown cells (Fig. 1d-f, and Additional file 1: Figure S1a and b). Interestingly, the ALKBH3 gene promoter is bound by the ALKBH3 protein (Fig. 1f, left panel) and belongs to the promoter group I with strong ALKBH3 enrichment.
In order to examine the general properties of ALKBH3 bound promoters, we compared ALKBH3 binding with publicly available datasets for histone marks, bound proteins, and other features (Additional file 1: Table S3). To examine ALKBH3’s role in cancer we preferentially used data from prostate (for example, LNCaP) or other human cancer cell lines. ALKBH3 bound promoters have elevated gene expression (Fig. 2a), are enriched for CpG islands (Fig. 2b), are bound by multiple transcription factors (TFs) (Fig. 2c), and are enriched for the active histone mark H3K4me3 (Fig. 2e) and RNA Polymerase II (Fig. 2f). In addition, ALKBH3 bound promoters are depleted for the TATA-Box and have highly ordered nucleosome positioning (Fig. 2d and g). These properties are typically found at ubiquitously expressed, house-keeping promoters  and gene ontology analysis using DAVID  confirmed a strong enrichment of genes involved in general cellular processes, such as translation, RNA processing, and cell metabolism (Fig. 2h). A motif search identified the ETS transcription factor binding motif to be most significantly enriched (Fig. 2i), suggesting that genes occupied by ALKBH3 are possibly activated by one or several ETS factors (Additional file 1: Figure S2), which are often over-expressed in prostate cancer [35, 36]. However, other TFs binding motifs are highly enriched as well (Fig. 2i), supporting the general conclusion that ALKBH3 bound promoters are strongly regulated and transcriptionally active.
Depletion of ALKBH3 induces an inflammatory response
Previous work for ALKBH1, ALKBH2, and ALKBH4 suggested that oxidative DNA demethylases could be directly involved in gene regulation [37–39]. To address if ALKBH3 affects target gene transcription, we knocked down ALKBH3 and extracted mRNA for microarray analysis at two different time points after selection (48 h and 96 h). We consistently identified 150 upregulated and 58 downregulated genes (Fig. 3a). However, unlike ALKBH1  or ALKBH2  depletion, genes bound by ALKBH3 did not show any significant expression level changes upon ALKBH3 knockdown (Fig. 3b).
To gain deeper insight into which genes are affected by ALKBH3 knockdown, we performed gene ontology (GO) analysis of the upregulated genes and found that many genes are involved in the innate/inflammatory immune response (Fig. 3a, c, and e). We confirmed, by qRT-PCR, the specific upregulation of several inflammatory genes after ALKBH3 but not ALKBH2 knockdown (Fig. 3d). Interferon stimulated genes (ISGs) are strongest induced at 48 h, but their expression declines at 96 h after ALKBH3 knockdown, indicating a negative feedback loop could exist (Fig. 3e).
Since most inflammatory genes are not directly bound by ALKBH3, we speculate that the inflammatory response may be an indirect consequence of DNA damage accumulation in the ALKBH3 deficient cells , because DNA damage triggers the activation of inflammatory and other pathways [3–5, 25, 40]. Interestingly, a similar phenomenon has been described in HeLa cells for the knockdown of the ALKBH3 interacting  ASCC3 DNA helicase . ALKBH3 knockdown also causes an upregulation of inflammatory genes in the non-small cell lung adenocarcinoma cell line NCI-H23 (Additional file 1: Figure S1c), suggesting that ablation of ALKBH3-dependent DNA repair mechanisms induces inflammatory pathways in multiple cancer cell lines.
Together these data suggest that ALKBH3 does not function to regulate transcriptional activity. Instead, the induction of inflammatory genes supports a potential role of ALKBH3 at its genomic targets to remove DNA alkylation adducts, such as 3-meC or 1-meA [16, 18, 19, 24], in order to prevent DNA damage.
ALKBH3 occupies places with ssDNA
Next we wished to determine if ALKBH3 was bound to specific DNA regions. Previous work showed that for DNA demethylation ALKBH3 prefers to demethylate ssDNA over dsDNA substrates [16, 18, 19]. This raises the possibility that for its DNA repair function, ALKBH3 is recruited to genomic places with ssDNA. At promoters, ssDNA is created during transcription initiation and at paused RNA polymerases II complexes. We asked whether ALKBH3 occupancy could be associated with those features. The ALKBH3 binding profile of group I promoters is characterized by a very strong enrichment slightly upstream of the transcription start site, while the group II promoter profile peaks mainly downstream of the TSS and shows only a moderate ALKBH3 binding level (Fig. 4a). Since the two distinct positions and binding strengths might reflect different ALKBH3 targets, we analyzed these two groups further separately. Interestingly, group I promoters have higher transcriptional activity compared to group II promoters (Additional file 1: Figure S3e and f).
The group II promoter profile is similar to the RNA polymerase II profile - both peak at +60 bp and possess a ‘shoulder’ upstream of the transcription start site (Fig. 4a and b). This correlation suggests that ALKBH3 recruitment could be dependent on the formation of single-stranded transcription bubbles at places with paused RNA polymerase II. To investigate this possibility, we divided all promoters (excluding group I promoters) into five different classes according their transcription levels. RNA polymerase II occupancy and ALKBH3 binding showed an almost linear correlation (Fig. 4b). These findings are consistent with the idea that ALKBH3 might get recruited to the transcription bubble of paused RNA polymerase II, analogous to what has been hypothesized before for AID (Activation Induced Cytidine Deaminase) .
Interestingly, the ALKBH3 group I promoter profile does not overlap significantly with RNA polymerase II. Instead, the ALKBH3 binding profile culminates slightly upstream of the TSS, at the site of transcription initiation (Fig. 4a). We speculated that at group I promoters ALKBH3 might target the single stranded DNA bubble established upon transcription initiation. DNA unwinding and formation of ssDNA during initiation depends on the DNA helicase XPB . Investigation of published ChIP-Seq data for XPB from HT1080 cells  revealed that XPB also has stronger binding at group I relative to group II promoters. Furthermore, XPB is enriched upstream of the TSS at group I but not at group II promoters, similar to ALKBH3 (Fig. 4c). These observations imply a stronger activity of XPB at the initiation site of group I promoters, which may lead to an elevated occurrence of ssDNA. We hypothesize that this increased ssDNA level might be the basis for the increased ALKBH3 recruitment. Notably, HT1080 cells are fibrosarcoma cells, but despite the different origin of PC3 and HT1080 cells, we still see this correlation of ALKBH3 and XPB binding. We reasoned that this promoter type may be present in a broad range of cell types, and we therefore analyzed these promoters bioinformatically using additional datasets (see below).
Single-stranded DNA occurs not only at promoters, but also at other places in the genome. During transcription elongation, the transcription bubble moves along the gene. Interestingly, highly expressed genes, which have elevated transcription elongation, have elevated ALKBH3 levels in the gene body (Additional file 1: Figure S4a). We also observed a mild ALKBH3 enrichment at enhancers where RNA polymerase II transcription has been reported  (Additional file 1: Figure S4b and d). Recently the role of G-quadruplex DNA (G4 DNA) during genome-wide regulatory processes was described . G4 DNA contains an accumulation of guanines that leads to the formation of a stable DNA quadruplex and the occurrence of ssDNA, in particular on the reverse strand . We found ALKBH3 is mildly enriched at places with potential G4 DNA (Additional file 1: Figure S4c and d).
Together these findings support the hypothesis that ALKBH3 is recruited to places in the genome with accessible ssDNA (summarized in Fig. 4d), in agreement with ALKBH3’s repair function at ssDNA [6, 16, 18, 19].
Group I promoters are a putative novel hyperactive promoter class
The binding of ALKBH3 at the transcription initiation site of group I promoters (Fig. 4a) prompted us to examine this promoter group in more detail. Interestingly, comparison of group I and group II promoters revealed differences in TF binding frequency (Additional file 1: Figure S3h). We found that group II promoters have on average 75 TF binding events while group I promoters have on average 112 TF binding events (50 % more). ETS transcription factors display a more pronounced difference (80 % more) (Additional file 1: Figure S2c). This indicates a clustering of transcription factors at these promoters, which is commonly associated with cohesin, the Mediator complex and other features [47, 48]. To further characterize the group I promoters, we performed a comprehensive analysis of numerous transcription associated features (Fig. 5a, Additional file 1: Figure S5). A large proportion of the investigated features are enriched at group I relative to group II promoters, indicating non-typical regulation (Fig. 5a). TAF1, a major subunit of TFIID, is one the most enriched factor (50 % more). In contrast, TBP (TATA-Box binding protein), another subunit of TFIID, is not significantly enriched, suggesting that at those promoters a TBP-independent recruitment of TFIID may occur more often than on other promoters . Further, the helicase CHD7, but not CHD1, CHD2, and CHD4, is strongly enriched (30 % more), raising the possibility that CHD7 could play a role in unwinding the DNA, in addition to XPB. Most active histone marks are mildly enriched, while the repressive histone mark H3K27me3 and its methyltransferase EZH2 are depleted. Interestingly, the histone modification with the strongest increase at group I promoters is H3K122 acetylation (40 % more), which has been demonstrated to enhance nucleosome eviction (Fig. 5a and b) . Since H3K122 acetylation is mediated by the histone acetyltransferase p300, and we also found a significant increase of p300, we speculate that acetylation of H3K122 by p300 and subsequent histone eviction could be of particular importance at these promoters . This idea is further supported by a reduced nucleosome occupancy and increased DNase I hypersensitivity at the initiation site of group I promoters (Additional file 1: Figure S3c and d). In addition to histone acetyltransferases, enzymes that regulate histone methylation, such as KMT2D (deposits H3K4me3) and PHF8 (removes H3K9me1/2, H3K27me2, and H4K20me1) are enriched as well (Fig. 5a and b).
Taken together, these data suggest that numerous activating factors are highly enriched at group I promoters, indicating a massive regulation, which consequently leads to higher transcriptional activity (‘hyperactive’) (Fig. 5c and d; Additional file 1: Figure S3e and f). A list of genes with this promoter type in PC3 cells is presented in Table S4 (Additional file 1). The presence of ALKBH3 at some of those gene promoters was confirmed in NCI-H23 lung cancer cells (Additional file 1: Figure S1d), suggesting that the genomic binding pattern of ALKBH3 as well as the occurrence of hyperactive promoters is similar in ALKBH3 over-expressing cancer cells.
Previous work suggested a pivotal role of ALKBH3 in multiple cancer types, such as prostate [24, 51], pancreatic , urothelial , non-small-cell lung , papillary thyroid , and brain  cancer. However, the cellular role of ALKBH3 is not yet fully understood. In order to gain further insights into the function of ALKBH3 particular in cancer, we applied genome wide approaches using PC3 prostate cancer cells as a model.
We initially performed genome localization studies to identify genomic regions bound by ALKBH3. We found that ALKBH3 preferentially occupies locations where ssDNA is occurring, such as promoters, enhancers, and G4 DNA (Fig. 4). This finding suggests that ALKBH3 is directed to ssDNA regions. In the future it will be exciting to determine whether this correlation is reflective of a regulated process and what the precise mechanism for ALKBH3 recruitment could be. Does this recruitment depend on ssDNA formation?
ALKBH3 was most enriched at the initiation site of a small number of highly expressed genes (Fig. 4a). Further investigation of these promoters led to the hypothesis that they are a putative novel ‘hyperactive’ subgroup of ubiquitously expressed gene promoters (Fig. 5).
Characterization of the ALKBH3 bound promoters was carried out using publicly available data from different research groups and different cell types. Despite use of this relatively heterogeneous dataset, we were still able to detect a correlation between strong ALKBH3 binding at group I promoters and enrichment of factors involved in gene activation. This study highlights how publicly available datasets can be utilized to develop novel hypotheses, while creating these data ab initio would neither be timely nor financially feasible.
Since we see this correlation of highest ALKBH3 binding with enriched binding of other factors across multiple cancer cell types we conclude that these hyperactive promoters exist in a relatively cell type independent manner. However, we cannot exclude that a transition between group I and group II promoters can occur and that a group II promoter might get promoted to a hyperactive promoter, and vice versa, once the level of activating factors at the promoter exceeds or falls below a certain threshold, respectively. Since we hardly see promoters where the ALKBH3 binding is in an intermediate state between group I and group II promoters, we hypothesize that these transitions – if they take place – are rare. It also remains to be determined whether these hyperactive promoters are restricted to highly proliferative cancer cells or if they occur in all cell types.
ALKBH3 binding to group I promoters is about eight-fold stronger than to group II promoters. We hypothesize that an elevated formation or an increased accessibility of ssDNA at the initiation site of group I promoters is the source for the increased ALKBH3 binding (Fig. 4d, left panel). We speculate that an enhanced presence of DNA helicases, such as the TFIIH helicases XPB and XPD, CHD7, or the ALKBH3 interacting ASCC3 helicase could lead to increased exposed ssDNA and therefore facilitate ALKBH3 recruitment [24, 57, 58]. The formation of ssDNA could also rely on bending of DNA upon binding of transcription factors, such as ETS factors and YY1 [59, 60]. It is also possible that group I promoters do not have increased accessible ssDNA and that the recruitment of ALKBH3 is facilitated by some unknown mechanism that does not depend on ssDNA.
We performed microarray analysis to elucidate whether ALKBH3 affects transcription of its target genes (Fig. 3). We observed no significant changes of the transcription of the ALKBH3 bound genes, suggesting that ALKBH3 binding does not directly regulate transcription. However, since ALKBH3 bound genes are highly expressed, subtle changes caused by ALKBH3 depletion might not be detected via microarray. Thus, our results do not explicitly exclude the possibility that ALKBH3 plays a similar expression regulating role as described for ALKBH1 and ALKBH2 [38, 39].
Knockdown of ALKBH3 in PC3 cells leads to an induction of inflammatory response gene expression (Fig. 3), which might be a consequence of elevated 3-meC levels and DNA damage after ALKBH3 depletion [3, 4, 24, 25, 61]. Most cancer cells, including PC3 cells, have rapid proliferation and accordingly elevated transcriptional activity . One consequence of elevated transcription is increased sensitivity of DNA for DNA alkylation damage, including 3-meC [6, 63–65]. The ALKBH3 genomic binding profile suggests that ALKBH3 is an intrinsic DNA repair protein that removes DNA alkylation that might occur naturally during transcription [63, 66]. ALKBH3 upregulation in cancer could be an important step during cancerogenesis to achieve an increased proliferation rate while maintaining genomic integrity (Fig. 6). If true, this would suggest that ALKBH3 inhibition could potentially slow cancer progression .
Since ALKBH3 demethylates RNA in addition to ssDNA [16, 18], it is possible that ALKBH3 could also have a role in the demethylation of primary RNA transcripts at its genomic targets . Whether this could be of physiological relevance will be of interest to be determined in the future.
Alkylating agents that methylate DNA and disrupt genomic integrity of fast proliferating cells are widely used in cancer treatment . Due to its localization at transcription-associated loci, we suggest that ALKBH3 is an intrinsic DNA repair protein that suppresses transcription associated DNA alkylation damage at highly expressed genes. The over-expression of ALKBH3 in cancer cells might facilitate alkylation damage resistance during cancer treatment and therefore raises the possibility of ALKBH3 as a potential anticancer target in the future [23, 67].
The genome-wide binding pattern of ALKBH3 revealed a strong binding to the initiation sites of a small number of highly active promoters. We hypothesize that these promoters are a new class of ubiquitously expressed promoters, which may have a specific initiation process, allowing ALKBH3 to ‘mark’ these promoters. Follow-up investigation of those hyperactive promoters in vitro, in vivo, and via bioinformatics will help to better understand their mechanisms as well as their role in transcription regulation and during cancer progression.
Lindahl T, Wood RD. Quality control by DNA repair. Science. 1999;286:1897–905.
Vijg J, Suh Y. Genome instability and aging. Annu Rev Physiol. 2013;75:645–68.
Gasser S, Raulet D. The DNA damage response, immunity and cancer. Semin Cancer Biol. 2006;16:344–7.
Chatzinikolaou G, Karakasilioti I, Garinis GA. DNA damage and innate immunity: links and trade-offs. Trends Immunol. 2014;35:429–35.
Janssens S, Tschopp J. Signals from within: the DNA-damage-induced NF-kappaB response. Cell Death Differ. 2006;13:773–84.
Sedgwick B. Repairing DNA-methylation damage. Nat Rev Mol Cell Biol. 2004;5:148–57.
Yi C, He C. DNA repair by reversal of DNA damage. Cold Spring Harb Perspect Biol. 2013;5:a012575.
Rydberg B, Lindahl T. Nonenzymatic methylation of DNA by the intracellular methyl group donor S-adenosyl-L-methionine is a potentially mutagenic reaction. EMBO J. 1982;1:211–6.
Barrows LR, Magee PN. Nonenzymatic methylation of DNA by S-adenosylmethionine in vitro. Carcinogenesis. 1982;3:349–51.
Roberts SA, Sterling J, Thompson C, Harris S, Mav D, Shah R, et al. Clustered mutations in yeast and in human cancers can arise from damaged long single-strand DNA regions. Mol Cell. 2012;46:424–35.
Yang Y, Sterling J, Storici F, Resnick MA, Gordenin DA. Hypermutability of damaged single-strand DNA formed at double-strand breaks and uncapped telomeres in yeast Saccharomyces cerevisiae. PLoS Genet. 2008;4:e1000264.
Chan K, Sterling JF, Roberts SA, Bhagwat AS, Resnick MA, Gordenin DA. Base damage within single-strand DNA underlies in vivo hypermutability induced by a ubiquitous environmental agent. PLoS Genet. 2012;8:e1003149.
Yang Y, Gordenin DA, Resnick MA. A single-strand specific lesion drives MMS-induced hyper-mutability at a double-strand break in yeast. DNA Repair (Amst). 2010;9:914–21.
Mosammaparast N, Shi Y. Reversal of histone methylation: biochemical and molecular mechanisms of histone demethylases. Annu Rev Biochem. 2010;79:155–79.
Delaney JC, Essigmann JM. Mutagenesis, genotoxicity, and repair of 1-methyladenine, 3-alkylcytosines, 1-methylguanine, and 3-methylthymine in alkB Escherichia coli. Proc Natl Acad Sci U S A. 2004;101:14051–6.
Falnes PO, Bjoras M, Aas PA, Sundheim O, Seeberg E. Substrate specificities of bacterial and human AlkB proteins. Nucleic Acids Res. 2004;32:3456–61.
Trewick SC, Henshaw TF, Hausinger RP, Lindahl T, Sedgwick B. Oxidative demethylation by Escherichia coli AlkB directly reverts DNA base damage. Nature. 2002;419:174–8.
Aas PA, Otterlei M, Falnes PO, Vagbo CB, Skorpen F, Akbari M, et al. Human and bacterial oxidative demethylases repair alkylation damage in both RNA and DNA. Nature. 2003;421:859–63.
Duncan T, Trewick SC, Koivisto P, Bates PA, Lindahl T, Sedgwick B. Reversal of DNA alkylation damage by two human dioxygenases. Proc Natl Acad Sci U S A. 2002;99:16660–5.
Dinglay S, Trewick SC, Lindahl T, Sedgwick B. Defective processing of methylated single-stranded DNA by E coli AlkB mutants. Genes Dev. 2000;14:2097–105.
Sperka T, Wang J, Rudolph KL. DNA damage checkpoints in stem cells, ageing and cancer. Nat Rev Mol Cell Biol. 2012;13:579–90.
Fu D, Calvo JA, Samson LD. Balancing repair and tolerance of DNA damage caused by alkylating agents. Nat Rev Cancer. 2012;12:104–20.
Helleday T, Petermann E, Lundin C, Hodgson B, Sharma RA. DNA repair pathways as targets for cancer therapy. Nat Rev Cancer. 2008;8:193–204.
Dango S, Mosammaparast N, Sowa ME, Xiong LJ, Wu F, Park K, et al. DNA unwinding by ASCC3 helicase is coupled to ALKBH3-dependent DNA alkylation repair and cancer cell proliferation. Mol Cell. 2011;44:373–84.
Hartlova A, Erttmann SF, Raffi FA, Schmalz AM, Resch U, Anugula S, et al. DNA damage primes the type I interferon system via the cytosolic DNA sensor STING to promote anti-microbial innate immunity. Immunity. 2015;42:332–43.
Lan F, Bayliss PE, Rinn JL, Whetstine JR, Wang JK, Chen S, et al. A histone H3 lysine 27 demethylase regulates animal posterior development. Nature. 2007;449:689–94.
Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
Liu T, Ortiz JA, Taing L, Meyer CA, Lee B, Zhang Y, et al. Cistrome: an integrative platform for transcriptional regulation studies. Genome Biol. 2011;12:R83.
Liefke R, Borggrefe T. Probing the epigenetic status at Notch target genes. Methods Mol Biol. 2014;1187:255–76.
Du Z, Zhao Y, Li N. Genome-wide colonization of gene regulatory elements by G4 DNA motifs. Nucleic Acids Res. 2009;37:6784–98.
Lenhard B, Sandelin A, Carninci P. Metazoan promoters: emerging characteristics and insights into transcriptional regulation. Nat Rev Genet. 2012;13:233–45.
Da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
Sharrocks AD. The ETS-domain transcription factor family. Nat Rev Mol Cell Biol. 2001;2:827–37.
Kumar-Sinha C, Tomlins SA, Chinnaiyan AM. Recurrent gene fusions in prostate cancer. Nat Rev Cancer. 2008;8:497–511.
Bjornstad LG, Meza TJ, Otterlei M, Olafsrud SM, Meza-Zepeda LA, Falnes PO. Human ALKBH4 interacts with proteins associated with transcription. PLoS One. 2012;7:e49045.
Li P, Gao S, Wang L, Yu F, Li J, Wang C, et al. ABH2 couples regulation of ribosomal DNA transcription with DNA alkylation repair. Cell Rep. 2013;4:817–29.
Ougland R, Lando D, Jonson I, Dahl JA, Moen MN, Nordstrand LM, et al. ALKBH1 is a histone H2A dioxygenase involved in neural differentiation. Stem Cells. 2012;30:2672–82.
Brzostek-Racine S, Gordon C, Van Scoy S, Reich NC. The DNA damage response induces IFN. J Immunol. 2011;187:5336–45.
Li J, Ding SC, Cho H, Chung BC, Gale Jr M, Chanda SK, et al. A short hairpin RNA screen of interferon-stimulated genes identifies a novel negative regulator of the cellular antiviral response. MBio. 2013;4:e00385–13.
Chaudhuri J, Tian M, Khuong C, Chua K, Pinaud E, Alt FW. Transcription-targeted DNA deamination by the AID antibody diversification enzyme. Nature. 2003;422:726–30.
Tirode F, Busso D, Coin F, Egly JM. Reconstitution of the transcription factor TFIIH: assignment of functions for the three enzymatic subunits, XPB, XPD, and cdk7. Mol Cell. 1999;3:87–95.
Gray LT, Vallur AC, Eddy J, Maizels N. G quadruplexes are genomewide targets of transcriptional helicases XPB and XPD. Nat Chem Biol. 2014;10:313–8.
De Santa F, Barozzi I, Mietton F, Ghisletti S, Polletti S, Tusi BK, et al. A large fraction of extragenic RNA pol II transcription sites overlap enhancers. PLoS Biol. 2010;8:e1000384.
Maizels N, Gray LT. The G4 genome. PLoS Genet. 2013;9:e1003468.
Yan J, Enge M, Whitington T, Dave K, Liu J, Sur I, et al. Transcription factor binding in human cells occurs in dense clusters formed around cohesin anchor sites. Cell. 2013;154:801–13.
Borggrefe T, Yue X. Interactions between subunits of the Mediator complex with gene-specific transcription factors. Semin Cell Dev Biol. 2011;22:759–68.
Wieczorek E, Brand M, Jacq X, Tora L. Function of TAF(II)-containing complex without TBP in transcription by RNA polymerase II. Nature. 1998;393:187–91.
Tropberger P, Pott S, Keller C, Kamieniarz-Gdula K, Caron M, Richter F, et al. Regulation of transcription through acetylation of H3K122 on the lateral surface of the histone octamer. Cell. 2013;152:859–72.
Koike K, Ueda Y, Hase H, Kitae K, Fusamae Y, Masai S, et al. anti-tumor effect of AlkB homolog 3 knockdown in hormone- independent prostate cancer cells. Curr Cancer Drug Targets. 2012;12:847–56.
Yamato I, Sho M, Shimada K, Hotta K, Ueda Y, Yasuda S, et al. PCA-1/ALKBH3 contributes to pancreatic cancer by supporting apoptotic resistance and angiogenesis. Cancer Res. 2012;72:4829–39.
Shimada K, Fujii T, Tsujikawa K, Anai S, Fujimoto K, Konishi N. ALKBH3 contributes to survival and angiogenesis of human urothelial carcinoma cells through NADPH oxidase and tweak/Fn14/VEGF signals. Clin Cancer Res. 2012;18:5247–55.
Tasaki M, Shimada K, Kimura H, Tsujikawa K, Konishi N. ALKBH3, a human AlkB homologue, contributes to cell survival in human non-small-cell lung cancer. Br J Cancer. 2011;104:700–6.
Neta G, Brenner AV, Sturgis EM, Pfeiffer RM, Hutchinson AA, Aschebrook-Kilfoy B, et al. Common genetic variants related to genomic integrity and risk of papillary thyroid cancer. Carcinogenesis. 2011;32:1231–7.
Cetica V, Genitori L, Giunti L, Sanzo M, Bernini G, Massimino M, et al. Pediatric brain tumors: mutations of two dioxygenases (hABH2 and hABH3) that directly repair alkylation damage. J Neurooncol. 2009;94:195–201.
Jung DJ, Sung HS, Goo YW, Lee HM, Park OK, Jung SY, et al. Novel transcription coactivator complex containing activating signal cointegrator 1. Mol Cell Biol. 2002;22:5203–11.
Kim HJ, Yi JY, Sung HS, Moore DD, Jhun BH, Lee YC, et al. Activating signal cointegrator 1, a novel transcription coactivator of nuclear receptors, and its cytosolic localization under conditions of serum deprivation. Mol Cell Biol. 1999;19:6323–32.
Pio F, Kodandapani R, Ni CZ, Shepard W, Klemsz M, McKercher SR, et al. New insights on DNA recognition by ets proteins from the crystal structure of the PU.1 ETS domain-DNA complex. J Biol Chem. 1996;271:23329–37.
Natesan S, Gilman MZ. DNA bending and orientation-dependent function of YY1 in the c-fos promoter. Genes Dev. 1993;7:2497–509.
Muralidharan S, Mandrekar P. Cellular stress response and innate immune signaling: integrating pathways in host defense and inflammation. J Leukoc Biol. 2013;94:1167–84.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144:646–74.
Kim N, Jinks-Robertson S. Transcription as a source of genome instability. Nat Rev Genet. 2012;13:204–14.
Fix D, Canugovi C, Bhagwat AS. Transcription increases methylmethane sulfonate-induced mutations in alkB strains of Escherichia coli. DNA Repair (Amst). 2008;7:1289–97.
Wrzesinski M, Nieminuszczy J, Sikora A, Mielecki D, Chojnacka A, Kozlowski M, et al. Contribution of transcription-coupled DNA repair to MMS-induced mutagenesis in E. coli strains deficient in functional AlkB protein. Mutat Res. 2010;688:19–27.
Aguilera A, Garcia-Muse T. R loops: from transcription byproducts to threats to genome stability. Mol Cell. 2012;46:115–24.
Nakao S, Mabuchi M, Shimizu T, Itoh Y, Takeuchi Y, Ueda M, et al. Design and synthesis of prostate cancer antigen-1 (PCA-1/ALKBH3) inhibitors as anti-prostate cancer drugs. Bioorg Med Chem Lett. 2014;24:1071–4.
LeRoy G, Chepelev I, DiMaggio PA, Blanco MA, Zee BM, Zhao K, et al. Proteogenomic characterization and mapping of nucleosomes decoded by Brd and HP1 proteins. Genome Biol. 2012;13:R68.
Asangani IA, Dommeti VL, Wang X, Malik R, Cieslik M, Yang R, et al. Therapeutic targeting of BET bromodomain proteins in castration-resistant prostate cancer. Nature. 2014;510:278–82.
Galbraith MD, Allen MA, Bensard CL, Wang X, Schwinn MK, Qin B, et al. HIF1A employs CDK8-mediator to stimulate RNAPII elongation in response to hypoxia. Cell. 2013;153:1327–39.
He HH, Meyer CA, Chen MW, Jordan VC, Brown M, Liu XS. Differential DNase I hypersensitivity reveals factor-dependent chromatin dynamics. Genome Res. 2012;22:1015–25.
Yu J, Mani RS, Cao Q, Brenner CJ, Cao X, Wang X, et al. An integrated network of androgen receptor, polycomb, and TMPRSS2-ERG gene fusions in prostate cancer progression. Cancer Cell. 2010;17:443–54.
Hollenhorst PC, Ferris MW, Hull MA, Chae H, Kim S, Graves BJ. Oncogenic ETS proteins mimic activated RAS/MAPK signaling in prostate cells. Genes Dev. 2011;25:2147–57.
Chen Y, Chi P, Rockowitz S, Iaquinta PJ, Shamu T, Shukla S, et al. ETS factors reprogram the androgen receptor cistrome and prime prostate tumorigenesis in response to PTEN loss. Nat Med. 2013;19:1023–9.
Xu K, Wu ZJ, Groner AC, He HH, Cai C, Lis RT, et al. EZH2 oncogenic activity in castration-resistant prostate cancer cells is Polycomb-independent. Science. 2012;338:1465–9.
Wang D, Garcia-Bassets I, Benner C, Li W, Su X, Zhou Y, et al. Reprogramming transcription by distinct classes of enhancers functionally defined by eRNA. Nature. 2011;474:390–4.
Lin CY, Loven J, Rahl PB, Paranal RM, Burge CB, Bradner JE, et al. Transcriptional amplification in tumor cells with elevated c-Myc. Cell. 2012;151:56–67.
Hu D, Gao X, Morgan MA, Herz HM, Smith ER, Shilatifard A. The MLL3/MLL4 branches of the COMPASS family function as major histone H3K4 monomethylases at enhancers. Mol Cell Biol. 2013;33:4745–54.
Tan PY, Chang CW, Chng KR, Wansa KD, Sung WK, Cheung E. Integration of regulatory networks by NKX3-1 promotes androgen-dependent prostate cancer survival. Mol Cell Biol. 2012;32:399–414.
ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489:57–74.
Kundaje A, Kyriazopoulou-Panagiotopoulou S, Libbrecht M, Smith CL, Raha D, Winters EE. Ubiquitous heterogeneity and asymmetry of the chromatin environment at regulatory elements. Genome Res. 2012;22:1735–47.
Plotnik JP, Budka JA, Ferris MW, Hollenhorst PC. ETS1 is a genome-wide effector of RAS/ERK signaling in epithelial cells. Nucleic Acids Res. 2014;42:11928–40.
We thank Drs. Andres Blanco, Tilman Borggrefe, Eric L. Greer, and Yang Shi for critical reading the manuscript. The authors would like to thank all peer reviewers for their constructive comments that helped improve the manuscript. This work has been supported from German Research Foundation (DFG, DA 1213/4-1 to SD, LI 2057/1-1 to RL) and by the Open Access Publication Funds of the Göttingen University.
The authors declare that they have no competing interests.
SD, RL, and IWJ designed the study and interpreted the data. IWJ performed most experiments. SD, JG, GSR, and MG performed experiments and analyzed data. RL performed the bioinformatics analyses. FW contributed to the bioinformatics analysis. SD and RL drafted the manuscript. All authors read and approved the final manuscript.
Robert Liefke and Indra M. Windhof-Jaidhauser contributed equally to this work.
Tables S1 to S4 and Figures S1 to S5. Table S1. List of the used qRT-PCR primers for cDNA. Table S2. List of the used qRT-PCR primers for ChIP. Table S3. List of public data used for bioinformatics analysis [44, 47, 50, 68–82]. Table S4. List of hyperactive promoters in PC3 cells. Figure S1. ALKBH3 Immunofluorescence in PC3 cells, ChIP at additional ALKBH3 targets in PC3 cells, induction of inflammatory genes in NCI-H23 cells, ALKBH3 ChIP in NCI-H23 cells. Figure S2. Bioinformatic analysis of relationship of ETS transcription factors, ALKBH3 and gene expression. Figure S3. A more detailed analysis of promoters such as in Fig. 2. Figure S4. Shows enrichment of ALKBH3 in the gene body, at enhancers and places of G4 DNA. Figure S5. Shows whisker blots of the counted ChIP-Seq tags at promoters from all analyzed features in Fig. 5a.
About this article
Cite this article
Liefke, R., Windhof-Jaidhauser, I.M., Gaedcke, J. et al. The oxidative demethylase ALKBH3 marks hyperactive gene promoters in human cancer cells. Genome Med 7, 66 (2015). https://doi.org/10.1186/s13073-015-0180-0
- Transcription Factor Binding Site
- Active Histone Mark
- Transcription Bubble
- Cluster Transcription Factor Binding Site
- ALKB Family