- Open Access
Knockdown of NAT12/NAA30 reduces tumorigenic features of glioblastoma-initiating cells
Molecular Cancer volume 14, Article number: 160 (2015)
Glioblastoma (GBM) is the most common primary brain malignancy and confers a dismal prognosis. GBMs harbor glioblastoma-initiating cells (GICs) that drive tumorigenesis and contribute to therapeutic resistance and tumor recurrence. Consequently, there is a strong rationale to target this cell population in order to develop new molecular therapies against GBM. Accumulating evidence indicates that Nα-terminal acetyltransferases (NATs), that are dysregulated in numerous human cancers, can serve as therapeutic targets.
Microarrays were used to study the expression of several NATs including NAT12/NAA30 in clinical samples and stem cell cultures. The expression of NAT12/NAA30 was analyzed using qPCR, immunolabeling and western blot. We conducted shRNA-mediated knockdown of NAT12/NAA30 gene in GICs and studied the effects on cell viability, sphere-formation and hypoxia sensitivity. Intracranial transplantation to SCID mice enabled us to investigate the effects of NAT12/NAA30 depletion in vivo. Using microarrays we identified genes and biochemical pathways whose expression was altered upon NAT12/NAA30 down-regulation.
While decreased expression of the distal 3’UTR of NAT12/NAA30 was generally observed in GICs and GBMs, this gene was strongly up-regulated at the protein level in GBM and GICs. The increased protein levels were not caused by increased levels of the steady state mRNA but rather by other mechanisms. Also, shorter 3’UTR of NAT12/NAA30 correlated with poor survival in glioma patients.
As well, we observed previously not described nuclear localization of this typically cytoplasmic protein. When compared to non-silencing controls, cells featuring NAT12/NAA30 knockdown exhibited reduced cell viability, sphere-forming ability, and mitochondrial hypoxia tolerance. Intracranial transplantation showed that knockdown of NAT12/NAA30 resulted in prolonged animal survival.
Microarray analysis of the knockdown cultures showed reduced levels of HIF1α and altered expression of several other genes involved in the hypoxia response. Furthermore, NAT12/NAA30 knockdown correlated with expressional dysregulation of genes involved in the p53 pathway, ribosomal assembly and cell proliferation. Western blot analysis revealed reduction of HIF1α, phospho-MTOR(Ser2448) and higher levels of p53 and GFAP in these cultures.
NAT12/NAA30 plays an important role in growth and survival of GICs possibly by regulating hypoxia response (HIF1α), levels of p-MTOR (Ser2448) and the p53 pathway.
Glioblastoma multiforme (GBM) is the most common and the most aggressive primary brain cancer. The average annual age-adjusted incidence rate is ~3 per 100.000 person-years . Despite multimodal therapy and inclusion of temozolomide, the overall survival remains dismal [2, 3]. This grim scenario mandates a pivotal shift in our approaches to develop new therapy as well as understanding the underlying disease mechanisms. During the last decade the GBM field has shown interest in a subgroup of cells harbored within these tumors with stem cell characteristics [4–6] that forms invasive tumors, similar to the tumor of origin, upon orthotopic xenotransplantation . These cells are commonly referred to as glioblastoma-initiating cells (GICs) and deemed important for the characteristics of GBM. They possess enhanced invasive properties [8, 9], promote tumor angiogenesis [10–12] and are resistant to irradiation  and chemotherapy [14, 15]. Thus, there is a strong rationale for validating potential molecular targets in GICs.
Dysregulation of co- or post-translational protein modifications is a trait of many human cancers. Nα- or N-terminal (Nt-) acetylation is one of the most common covalent modifications of all soluble human proteins  and occurs predominantly co-translationally . N-terminal acetyltransferases (NATs) or Nα-acetyltransferases (NAAs) are arranged in complexes and believed to target ~80-90 % of all soluble human proteins [17, 18]. Six NAT complexes have been identified in humans, NAT-A – NAT-F and each complex consists of one catalytic subunit and auxiliary subunit(s) . With an increasing body of knowledge about their diverse functions, substrates and downstream targets, NATs are emerging as potential targets in several cancers . Knockdown studies targeting NATs in various malignancies including colon , thyroid  and hepatocellular cancers  have shown that depletion of the catalytic subunit causes a less aggressive phenotype with reduced cell proliferation and/or increased apoptosis.
NAT12/NAA30 is the catalytic subunit of the NAT-C complex . Knockdown of each of the NAT-C subunits, led to p53-dependent apoptosis in HeLa and colon carcinoma cell lines with the strongest phenotype observed with depletion of NAT12/NAA30 . The latter study also pointed out mammalian target of rapamycin (MTOR) as a substrate for NAT-C. Another report also suggested TOR as a target of NAT-C activity .
In the present study we investigated the expression of NAT12/NAA30 in GBM tissue samples, GICs, normal brain tissue, and neural stem cells (NSCs) from the adult human brain as well as in a neural fetal cell line (NFCs). Using immunolabeling, we revealed a hitherto undescribed nuclear localization of NAT12/NAA30 protein. To study the function of NAT12/NAA30, we performed gene knockdown using RNA interference (RNAi) technology. Knockdown of NAT12/NAA30 resulted in markedly reduced cell viability and sphere-forming ability of GICs. To study genes and pathways downstream of NAT12/NAA30, we used microarray analysis and western blot. This enabled us to identify several pathways such as p53, ribosomal assembly, hypoxia response and cell proliferation regulated by NAT12/NAA30. Furthermore, we documented a reduction of phospho-MTOR (Ser2448) and increased levels of p53 and glial fibrillary acidic protein (GFAP) in the knockdown cultures. We show that intracranial transplantations into severe combined immunodeficient (SCID) mice of GICs featuring NAT12/NAA30 knockdown, resulted in a significant prolongation of animal survival compared to controls.
Expression of NAT12/NAA30 in brain tissues, GBM, NSCs, GICs and NFCs
To investigate the expression of NAT12/NAA30 in GBM tumor biopsies, normal human brain, NSCs, GICs and NFCs we used microarrays, real-time quantitative reverse-transcription PCR (qPCR), western blots, immunolabeling and public data mining.
Microarray analysis was performed on GIC cultures from seven patients, ten NSC cultures from five patients, two GBM biopsies, two normal brain tissue samples and one NFC culture. NSC cultures were isolated from four brain regions: the subventricular zone (SVZ), hippocampus (HPC), white- and grey matter (WM and GM respectively). Normal tissues were from WM and GM.
NAT12/NAA30 expression (measured with the 3’ terminal ILMN_2128087 reporter) was moderately but significantly higher in NSC compared to GIC cultures (Fig. 1a). It was also higher in the cell cultures compared to tissues (Additional file 1: Figure S1). We also analyzed the expression of several other NATs in the same set of microarrays and found that NAT5, NAT14, NAT10 and NAT9 were among the most abundant NATs in the analyzed cell cultures and tissues (Fig. 1b). NAT12/NAA30 was among the relatively lowly expressed NATs (Fig. 1b). K-means clustering (Additional file 2: Figure S2) showed that NAT12/NAA30 was co-expressed with NAT1, NAT13 and NAT5. Hierarchical clustering with distance matrix using Pearson correlation as a distance measure (Fig. 1c) revealed that NAT12/NAA30 was highly co-expressed with NAT8, NAT8B, NAT1, NAT6, NAT13, and co-expressed with NAT2 and NAT5 (Fig. 1c). Several NATs from this group had significantly lower expression in GIC than in NSC cultures (NAT12, NAT6, NAT8 and NAT8B). Several other NATs such as NAT9, NAT10, NAT14 and NAT15 were expressed to a similar degree in GIC and NSC cultures. NAT8L was the only NAT that was significantly up-regulated in GIC cultures and GBM (Additional file 1: Figure S1).
qPCR analysis was used to validate the microarray results. We used seven sets of oligonucleotides that covered both coding (A-C) and non-coding (D-G) regions of the gene (Fig. 1d-e, Additional file 3) to investigate NAT12/NAA30 expression in seven GIC cultures, two NSC cultures and one NFC culture.
qPCR revealed that the relative expression of the coding regions of NAT12/NAA30 was very variable between different GIC cultures. However, the average expression values were similar in GIC and NSC cultures as shown by qPCR with oligonucleotides corresponding to the coding regions of the gene (Fig. 1d; Oligonucleotide sets A and B, Additional file 4) and Taqman probes corresponding to the exons (Fig. 1f).
Interestingly, when oligonucleotides corresponding to the 3’UTR (C, E, F and G) were used (Fig. 1e), qPCR analysis showed that NAT12/NAA30 expression was reduced in GIC cultures (Fig. 1d). We observed almost complete absence of expression (97 % reduction) in GIC cultures when using oligonucleotide set G, which corresponds to the most distal region of the 3’UTR (Fig. 1d, Additional file 4).
Western blot. NAT12/NAA30 protein expression in NSC and GIC cultures, GBM tissues and NFCs was investigated using western blots (Fig. 1g). We detected a stronger band at around 50 kDa in the GICs although the expected sizes of the two isoforms of the protein were 35 and 39 kDa (Fig. 1g). This might indicate either protein modification or the presence of a different isoform in the GICs. Bands of comparable sizes were detected using another antibody and recombinant protein confirming that the correct bands were detected (results not shown). Western blot analysis revealed that NAT12/NAA30 expression in GIC cultures was significantly higher than in NSC cultures (Fig. 1g) although this was not observed at RNA level (Fig. 1a, d and f).
Survival analysis. Public database mining was done using data from the REMBRANDT database (https://caintegrator.nci.nih.gov/rembrandt/) where 3 probes were used to measure the expression of NAT12/NAA30 covering different parts of the transcript (Fig. 1e; microarray probes ending with “_at”).
Interestingly, the 3’UTR of NAT12/NAA30 seemed to be universally down-regulated in glioma tissue samples. This was in accordance with our data obtained by both microarray analysis and qPCR (Fig. 1a and d). Glioma patients with decreased expression of the distal 3’UTR end of the NAT12/NAA30 transcript (Probe: 225679_at, Fig. 1e) had significantly shorter survival (Additional file 5: Figure S3A), while the increased expression of the coding regions of this gene (Probes: 228322_at, 228321_at, Fig. 1e) did not seem to correlate with patient survival (Additional file 5: Figure S3B-C). For comparison, we also checked the association between the expression levels of other NATs in glioma tissues and patient survival. Increased expression of each of NAT1, NAT2 and NAT10 correlated negatively with patient survival in the group of “all gliomas” patients (Additional file 5: Figure S3D-F).
Immunolabeling. We tested the expression of the NAT12/NAA30 protein in three GIC and three NSC cultures, and in frozen sections from two GBMs and two normal brain tissues (Fig. 2 and Additional file 6: FigureS4). NAT12/NAA30 protein was barely detectable in samples from normal brain (Fig. 2a, Additional file 6: Figure S4A) but it was highly up-regulated in the GBM biopsies (Fig. 2c, Additional file 6: Figure S4B). GBM tissue sections were co-stained with anti-NAT12/NAA30 and anti-nestin antibodies (Fig. 2d). This revealed a cell population expressing both proteins (Fig. 2e). The cells expressing both nestin and NAT12/NAA30 were detected frequently in GICs (Fig. 2h-j). Interestingly, in the GBM sections, NAT12/NAA30 expressing cells were located in the perivascular spaces (Fig. 2f). The same staining pattern was confirmed with another antibody against NAT12/NAA30 (results not shown). The observation that NAT12/NAA30 was highly expressed in cells surrounding blood vessels, was verified by co-staining with anti-CD31 antibody (Fig. 2g).
NAT12/NAA30 was previously detected in the cytoplasm of various eukaryotic cells [23–26]. However, our immunolabeling revealed an hitherto unreported nuclear localization of NAT12/NAA30. The nuclear localization was consistently detected in all tested GIC cultures but was more sporadic in the NSC cultures (Fig. 2k-n, Additional file 6: Figure S4C-F).
In order to study the functional relevance of NAT12/NAA30 in GICs, we established short hairpin RNA (shRNA)-mediated knockdown of NAT12/NAA30. We transduced the GIC culture T65 with lentiviruses harboring two different shRNA sequences and obtained transgenic cell cultures (KD1 and KD2). Non-silencing (NS) shRNA was used as a control. Knockdown efficiency was assessed with qPCR, western blot and immunolabeling to confirm depletion of NAT12/NAA30 transcript and -protein. The gene knockdown efficiency was measured with oligonucleotides corresponding to the coding (A-C) and non-coding (D-G) regions of the gene (Fig. 1e, Fig. 3a-b). Relative expression of NAT12/NAA30 was reduced with 60 ± 14 % in KD1 (Fig. 3a) and 56 ± 3 % in KD2 (Fig. 3b). For this calculation we used average values for two sets of oligonucleotides (B and C) corresponding to the coding regions of the gene (Fig. 1e, Additional file 4). Transcripts corresponding to the 3’UTR were predominantly not down-regulated in knockdown cultures KD1 and KD2 as shown by qPCR using oligonucleotides D-G (Fig. 3a-b).
Western blot showed reduction of NAT12/NAA30 protein levels in both KD1 and KD2 cultures (Fig. 3c). To quantify the efficiency of knockdown at the protein level, we used three western blots and two antibodies. We measured an average reduction of 27 ± 13 % in KD1 and 40 ± 6 % in KD2 (Fig. 3d).
Immunolabeling confirmed reduced protein expression in GICs featuring NAT12/NAA30 gene knockdown compared to both untreated control and NS (Fig. 3e). This was also confirmed with another antibody against NAT12/NAA30 (results not shown).
Knockdown of NAT12/NAA30 in GICs reduces cell viability and sphere-forming ability
We investigated the functional effects of the NAT12/NAA30 knockdown in GICs by measuring cell viability and sphere-forming ability in NS, KD1 and KD2 cultures.
Values obtained by the XTT assay for the knockdown cultures were normalized to the NS control. The number of viable cells was reduced to 58 ± 8 % for KD1 (p = 0.0314; unpaired t-test with Welch’s correction) and 43 ± 8 % (p = 0.0172) for KD2 compared to the NS control. This demonstrated a significant reduction in cell viability of GICs featuring NAT12/NAA30 knockdown (Fig. 4a).
The number of spheres formed from cells seeded at a low density (500 cells per well) was 29 ± 7 (p = 0.0472) in KD1 and 17 ± 2 (p = 0.121) in KD2 compared to 58 ± 10 in the NS control (Fig. 4b-c). Our results thus indicate that the NAT12/NAA30 knockdown reduced sphere-forming ability by 50–71 %. The sphere-size was, however, not affected (results not shown).
Knockdown of NAT12/NAA30 in GICs prolonged survival of intracranially xenografted mice
To investigate the effect of NAT12/NAA30 knockdown on GIC growth and tumor formation in vivo we used a robust orthotopic xenograft model [27, 28]. Cells from the knockdown cultures KD1 and KD2, and the NS control culture were transplanted into the right hemisphere of SCID mice (n = 6 for each group).
Median survival for mice transplanted with NS cells was 125 days versus 189 for KD1 (p = 0.0192; Log-rank test) and 143 for KD2 (p = 0.0385) (Fig. 4d). All 6/6 animals in the NS control group died whereas 2 animals in KD1 group and 1 in KD2 group did not show any sign of distress at the end of the observation period of 200 days. Representative cryosections stained with Hematoxylin and Eosin (HE) from the xenograft brains are shown in Additional file 6: Figure S4G.
Knockdown of NAT12/NAA30 in GICs decreased mitochondrial hypoxia tolerance
The ability to tolerate hypoxic conditions is a well known trait of GICs and it contributes to the tumor cells’ evasion of apoptosis . It has previously been shown that NAT12/NAA30 is involved in regulation of apoptosis . Mitochondrial depolarization is an early step in programmed cell death. To test if down-regulation of NAT12/NAA30 influences mitochondrial properties, such as mitochondrial hypoxia tolerance, we compared the mitochondrial response to acute severe hypoxia in NS and knockdown cultures, KD1 and KD2.
Baseline recordings of mitochondrial membrane potential (ΔΨm) were stable in all cells. Increase in Rh123 fluorescence, indicating mitochondrial membrane depolarization, was seen in all three groups within 3–5 min of hypoxia. Compared to the NS control, there was a steeper increment of the fluorescent signal in KD1 and KD2 (Fig. 4e-f; KD1 and KD2 vs. NS, both p < 0.01, one-way ANOVA). The remaining ΔΨm was released upon addition of the protonophore FCCP. The FCCP response of the control cells was almost twice that of the NAT12/NAA30 knockdown cells (relative fluorescence increase: NS 1.63; KD1 0.92; KD2 0.95) indicating that the remaining mitochondrial membrane potential after hypoxia was higher in the NS control cells than in KD1 and KD2. The difference was, however, not statistically significant (Fig. 4e-f; p = 0.07, two-way ANOVA).
Identification of genes and biological pathways downstream from NAT12/NAA30
To identify genes and biological pathways downstream from NAT12/NAA30, we used microarrays, qPCR and western blot.
Microarray analysis showed that 661 Illumina probes identified transcripts that were differentially regulated (1.5-fold change) in GICs featuring NAT12/NAA30 knockdown (Additional file 7). Of these 464 were recognized and analyzed by the DAVID functional annotation tool [30, 31]. This analysis revealed that the NAT12/NAA30 knockdown had a significant effect (p = 3.00E-05) on the expression of genes coding for ribosomal proteins (Additional file 7). Many of these were directly involved in protein translation (Additional file 7). KEGG-pathway analysis identified three significantly dysregulated pathways: 1) Ribosome (13 genes, p = 5.6E-06), 2) p53 pathway (8 genes, p = 3.1E-03) and 3) Sphingolipid metabolism (6 genes, p = 4.8E-03) (Additional file 8: Figure S5–S7 and Additional file 7).
A more detailed analysis of the p53 pathway revealed that the expression of eight genes in this pathway was differentially regulated in the NAT12/NAA30 knockdown cultures. These genes were involved in apoptosis, cell cycle arrest and DNA repair (Additional file 8: Figure S6 and Additional file 7). Former studies have shown that knockdown of the NAT-C complex leads to reduced proliferation and p53-dependent cell death in human cancer cell lines . Starheim et al. reported increased levels of total p53, phospho-p53 (Ser37) and apoptosis in HeLA cells featuring NAT12/NAA30 knockdown . In our knockdown cultures, western blot showed a slight increment in the total p53 (Fig. 5a), while the levels of the phospho-p53 (Ser37) were not altered (Additional Additional file 9: Figure S8A). We could not detect any significant changes in apoptosis in the knockdown cultures (results not shown).
Ten genes involved in the response to hypoxia (p = 6.80E-03; GOTERM_BP_FAT: response to hypoxia) were differentially regulated in the knockdown cultures (Additional file 10: Figure S9 and Additional file 7). Further analysis of the hypoxia response pathway showed that the hypoxia induced factor 1α (HIF1α), potassium large conductance calcium-activated channel (KCNMA1) and the integrin receptor (ITGA2) (involved in cell survival and growth via the PI3K pathway) were down-regulated in the knockdown cultures. HIF1α is one of the central regulators of transcriptional control in tumor cells’ response to hypoxic conditions. In GBM, it is regarded as a key player in tumor progression and an important factor for the maintenance of GIC phenotype [32, 33]. Therefore, it was of great interest to analyze the expression of HIF1α protein in the NAT12/NAA30 knockdown cultures. Notably, western blot showed that the levels of HIF1α were decreased in the knockdown cultures compared to the controls (Fig. 5d), thus confirming the microarray results. The levels of KCNMA1 protein were, however, not altered in the knockdown cultures (Additional file 9: Figure S8B). Among the up-regulated genes were the growth inhibiting factor metallothionein 3 (MT3), CD24, chemokine receptor type 4 (CXCR4) and chemokine ligand 2 (CCL2) (Additional file 10: Figure S9 and Additional file 7).
Among other genes differentially expressed in the NAT12/NAA30 knockdown cultures, were groups of genes coding for: a) proteins involved in organization of nucleosomes (11 genes, p = 1.1E-04), b) cell proliferation (22 genes, p = 3.50E-03) and c) protein kinase activity (22 genes, p = 1.8E-04) (Additional file 7).
Previous studies have identified MTOR as a potential substrate for the NAT-C complex  and a downstream target of NAT-C activity . Microarray analysis of the GIC cultures featuring NAT12/NAA30 knockdown, showed increased expression of insulin-like growth factor binding protein 3 (IGFBP3), the inhibitor of the IGF1/MTOR pathway as compared to NS controls (Additional file 8: Figure S6). Western blot showed that total MTOR was not changed while the levels of phospho-MTOR (Ser2448) were reduced in the knockdown cultures (Fig. 5a). Western blot analysis of selected important factors in the mTOR/AKT pathway is provided in Additional file 9: Figure S8C.
One of the key features of GICs is the high expression levels of the pluripotency marker nestin (NES) and that they can be differentiated into a neuronal or glial lineage [5, 27]. Upon differentiation the proliferative capacity of these cells is lost while the levels of the neuronal marker β3-tubulin (TUBB3) and the astrocytic marker GFAP increase . As we have shown NAT12/NAA30 knockdown reduced GIC proliferation (Fig. 4a-b). Thus, we wanted to investigate whether knockdown of NAT12/NAA30 altered the expression levels of the aforementioned markers towards a more mature cell type. Interestingly, in both knockdown cultures GFAP was highly up-regulated while the expression of NES and TUBB3 remained unchanged as shown by qPCR analysis (Fig. 5b). However, the expression of nestin at the protein level was reduced in the knockdown cultures (Additional file 9: Figure S8B). Western blot showed a so far undescribed variant of GFAP (30 kDa) highly upregulated in the knockdown cultures (Fig. 5c).
We have tested the expression of stem cell related genes in the NAT12/NAA30 knockdown cultures. Flow cytometric analysis showed that the CD133+ population was depleted in both knockdown cultures compared to the NS control (Additional file 11: Figure S10). We did not detect the expression of CD15 in T65 (results not shown). We also evaluated the expression of NANOG, POU5F1/OCT4, SOX2 and MYC using microarrays. None of these genes were differentially regulated at the mRNA level in the knockdown cultures (Additional file 9: Figure S8C). However, the levels of SOX2 protein were slightly decreased in the two knockdown cultures as shown by western blot (Additional file 9: Figure S8B).
Western blot analysis further revealed that the levels of active β-catenin (active CTNNB1/ABC), phospho-p70 S6 kinase (Thr389), phospho-STAT3 (Ser727), GLI1 and CCND1 and several other proteins were only slightly affected or not altered in the knockdown cultures (Additional file 8: Figure S8A-B).
GBM remains one of the most aggressive human cancers. Accumulating evidence supports the role of GICs in disease progression, tumor cell invasion and the tumor’s notorious resistance to chemo- and radiotherapy. Consequently, new molecular therapies against GBM would also need to target this cell population.
Recently, protein Nt-acetylation has been implicated in cancer development . In the present study we analyzed the expression of several NATs in GIC cultures and GBM tissues showing which of these are co-expressed and thus might be functionally redundant (Fig. 1b-c and Additional file 4: Figure S2). We also show that the increased expression of several NATs in GBM tissues is associated with poor patient survival and that these enzymes may be relevant as new potential therapeutic targets. To further reinforce the validity of our study we used NSC cultures isolated from several different regions of the brain as controls. We identified NATs whose expression was significantly down- or up-regulated in GBM and GICs (Additional file 3: Figure S1).
By detailed expression analysis of NAT12/NAA30 in tissues and cell cultures (GIC and NSC cultures, and NFCs), we detected the presence of alternative transcript(s) that had 3’UTR regions of variable length (Fig. 1d). We have demonstrated that these were differentially regulated in GICs compared to NSCs. Especially the distal 3’UTR of NAT12/NAA30 was strongly down-regulated in GBM and GIC cultures as shown by two independent methods (microarray and qPCR) (Fig. 1a and d). Moreover, public database mining revealed that the decreased expression of the distal 3’UTR region correlated with a shorter survival of glioma patients (Additional file 5: Figure S3A). Several studies have shown that mRNA isoforms with differences in their 3’UTRs have different stability and translational activity [34–36]. In yeast it was found that shorter mRNAs more frequently formed a closed-loop structure that enhanced protein translation . Mayr et al. studied 27 cancer cell lines from different tissues and showed that shortening of the 3’UTR in the tested mRNAs seemed to be important for activation of oncogenesis . Shortening of the 3’UTRs resulted in enhanced translation in some of these cancer cell lines. More importantly, shorter mRNAs had greater stability and produced more protein . Another study reported that shortening of 3’UTRs correlated with poor patient survival in breast and lung cancer . To our knowledge we are the first to report the existence of an analogous mechanism associated with GBM and GICs. In the present work we show that the up-regulation of NAT12/NAA30 at the protein level in GICs is not caused by the increased levels of the steady state mRNA but rather by other mechanisms such as shortening of the 3’UTR and possibly greater stability of mRNA. Other processes involved in the complex regulation of NAT12/NAA30 expression might be alternative splicing, alternative poly-adenylation sites (Stangeland, unpublished), and the presence of regulatory and/or coding sequences within the 3’UTR.
Human NAT12/NAA30 is known to be present in the cytoplasm [23–26]. Our data show that it is also located in the nuclei of GICs and some NSCs (Fig. 2i and k). The transport of proteins from the cytosol to the nucleus is a tightly regulated process facilitated by importins . NAT12/NAA30 contains a nuclear localization signal (RRGYIAMLAVDSKYRRN at 243 according to pSORT II, http://psort.hgc.jp/) that confers a high probability of being imported to the cell nuclei. It has also been reported that 3’UTRs of other genes contain regulatory regions that influence nuclear export and subcellular localization . The observed differential subcellular localization of NAT12/NAA30 in the nucleus is novel and might imply an additional role for NAT12/NAA30. Another member of the NAT family, NAT-D, is known to acetylate histones H2A and H4 . Whether NAT12/NAA30 protein can have a similar role remains to be shown.
Furthermore, evaluating the staining pattern of NAT12/NAA30 in GBM biopsies led to an interesting observation (Fig. 2b-f). NAT12/NAA30 was found to be expressed predominantly in cells surrounding blood vessels. In this perivascular region the GBM cells have a much stronger expression of nestin than in other parts of the tumor specimen. Co-staining with these two markers showed a large degree of overlap in the GBM cells around blood vessels. Recently, it was shown that nestin regulates stemness, cell growth and invasion in GBMs . The authors reported that the overexpression of NESTIN increased cell growth, sphere formation and cell invasion while depletion of NESTIN resulted in decreased expression of stem cell markers. The expression of nestin has been found to be an independent prognostic factor in glioma patients . Since GICs are believed to reside in the perivascular niche of the tumor [10–12], our immunolabeling might indicate that NAT12/NAA30 is predominantly expressed in an immature cell types in GBM.
RNAi-mediated gene silencing of NAT12/NAA30 enabled us to obtain a stable knockdown at both transcript and protein level as determined by qPCR, western blot and immunolabeling (Fig. 3). NATs are evolutionarily highly conserved [17, 22]. The partial knockdown we obtained might be explained by selection mechanisms having enabled GICs to avoid the lethality that would result from total NAT12/NAA30 knockdown. Nevertheless, functional analysis of the knockdown cultures revealed a clear phenotype with restricted cell growth in vitro and in vivo (Fig. 4a-d). Reduced sphere-forming ability in the knockdown cultures was stable across passages. This finding is especially important because sphere formation has been reported to be an independent predictor of clinical outcome .
The number of spheres is assumed to be both an indicator of the aggressiveness of the tumor and the number of tumor initiating cells. The ability to form spheres strongly correlates with tumor growth and survival also in animal models . Our results from intracranial transplantation into SCID mice were in agreement with this as mice transplanted with GICs featuring NAT12/NAA30 gene knockdown survived significantly longer than controls. Altogether 3/11 mice did not form tumors until the end of the observation period in the NAT12/NAA30 knockdown group while all 6/6 mice died from tumors in the control group.
There are few known substrates of NAT-C, but based on the Nt-amino acid composition of target proteins NAT-C can potentially Nt-acetylate up to 14.5 % of all cytoplasmic human proteins . The N-terminus of MTOR protein begins with Met-Leu and is a strong candidate as a direct substrate [23, 24]. In a study in zebrafish TOR was suggested as a downstream target of NAT-C . The authors also demonstrated that overexpression of TOR rescued the effect of depletion of a NAT-C subunit . Knockdown of MTOR is a critical effector of the PIK3-AKT pathway which is dysregulated in many cancers . In GBM AKT-signaling activity is significantly correlated with phosphorylation of MTOR .
In our study, knockdown of NAT12/NAA30 led to decreased levels of phospho-MTOR (Ser2448), while the total MTOR levels remained unchanged (Fig. 5a). Phospho-MTOR (Ser2448) binds to the effector complexes mTORC1 and mTORC2 and is important for mTORC1 activity . Implication of MTOR as a (direct or indirect) target of NAT12/NAA30 was further supported by the expression analysis that showed increased levels of IGFBP3, an inhibitor of IGF1/MTOR pathway, in the cell cultures featuring NAT12/NAA30 knockdown (Additional file 8: Figure S6).
Our results also indicate that NAT12/NAA30 acts via the p53 pathway in GICs. In addition to an increment in p53 protein levels in the knockdown cultures we also provide evidence for dysregulation of the p53 pathway genes. Genes involved in cell-cycle regulation, apoptosis and DNA repair were those that suffered the most substantial expressional alteration (Additional file 8: Figure S6). We also show that NAT12/NAA30 affects a considerable number of genes that regulate cell proliferation and protein kinases.
We investigated the expression of stem cell related genes in NAT12/NAA30 knockdown cultures. Of the tested stem cell markers, we could detect reduced levels of nestin and SOX2 at the protein level. The knockdown of NAT12/NAA30 also reduced the percentage of CD133+ cells. Our data thus indicate the effect of NAT12/NAA30 knockdown on stemness in GICs.
A previous study has shown that mitochondrial proteins are substrates of NAT12/NAA30 . Interestingly, we found that knockdown of NAT12/NAA30 in GICs causes a more abrupt and severe mitochondrial membrane depolarization compared to the NS control cultures when these cells are exposed to hypoxia. This indicates a reduced mitochondrial tolerance to acute hypoxia upon NAT12/NAA30 knockdown. This notion was further supported by microarray analysis showing that the role of NAT12/NAA30 in hypoxia response involves HIF1α and several other genes (Additional file 10: Figure S9). Importantly, western blot analysis confirmed the reduced expression of HIF1α protein (Fig. 5d). Hypoxia induced factors are key mediators in the response of cancer stem cells to hypoxia [32, 51]. Several studies have pointed out that GICs reside in the hypoxic regions of GBMs  and that hypoxia-induced molecular changes regulate and maintain the phenotype of GICs . This might imply that targeting NAT12/NAA30 increases the vulnerability of GICs to hypoxic conditions.
Although N-terminal acetylation complexes have lately received a lot of attention, the functional role of the NAT-C complex in cancer stem cells has not been investigated. In our work, we show that the knockdown of its catalytic subunit, NAT12/NAA30, reduced tumorigenic features of GICs and confirmed that MTOR and HIF1α are downstream targets of NAT12/NAA30 in GICs. We therefore propose that, NAT12/NAA30, might serve as a potential therapeutic target in GBM. Our study further indicates involvement of NAT12/NAA30 in gliogenesis and regulation of the p53 pathway. Specific inhibitors against NATs are under development  and might represent promising candidates for clinical testing.
Materials and methods
Tissue specimens and cell culture
Tissue specimens were harvested from consenting patients after approval by the Norwegian National Committee for Medical Research Ethics. Tumor biopsies were obtained as a part of surgical procedures for treating GBM. Normal brain tissue (SVZ, HPC, GM and WM) was obtained from fresh human temporal lobes surgically resected to treat medically refractory epilepsy. All biopsy specimens were evaluated by neuropathologists.
Tumor biopsies underwent mechanical dissociation and Trypsin-EDTA (Gibco, Life Technologies, NYC, NY, USA) was added for enzymatic dissociation. Subsequently, 2 mg/ml human albumin (Octapharma pharmazeutika produktionges, Vienna, Austria) was used to block the Trypsin effect and the cells were washed in L-15 (Lonza, Basel, Switzerland) before being plated in serum-free neurosphere medium containing 10 ng/ml basic fibroblast growth factor (bFGF) and 20 ng/ml epidermal growth factor (EGF) (both R&D Inc., Minneapolis, MN, USA), B27-supplement (1:50, Invitrogen, Carlsbad, CA, USA), 100 U/ml Penicillin/streptomycin (Lonza), 1 ng/ml Heparin (Leo Pharma, Ballerup, Denmark) and 8 mM Hepes (Lonza) in Dulbecco’s modified essential medium with nutrient mix F-12 and Glutamax (DMEM/F12, Invitrogen) [24, 25]. Dissociated GBM biopsies grown as free-floating spheres in serum free-medium containing mitogens (EGF and bFGF) are highly enriched for GICs [7, 54]. The GIC cultures used in this current work (T65, T08, T59) were analyzed for the expression of stem cell markers by flow cytometry and showed high expression of CD44, CXCR4, CD166 and CD9 while the expression of CD133 and CD15 was variable . Immunolabeling showed that most cells were nestin and Sox2 positive (Additional file 6: Figure S4H-I, Stangeland, B. et al., in press). The cells were cultured in 75 cm2 non-treated flasks (Nunc, Roskilde, Denmark) at a density of 105 cells/ml and supplemented with EGF and bFGF twice a week. When the spheres reached approximately 100 μm in diameter they were dissociated into single cells as previously described . We have previously shown that the GIC cultures T65, T08, TC3 and TC4 formed invasive tumors upon orthotopic transplantation to SCID mice [55, 56].
NFCs (ReNcell VM Human Neural Progenitor Cell Line, SCC008, Merck Millipore, Darmstadt, Germany) were cultured as spheres in serum-free Neurobasal A medium (Gibco) containing B27 (Gibco), 2 mM L-glutamine, 10 ng/ml bFGF, and 20 ng/ml EGF (both from R&D Systems).
Tissue specimens from the adult human brain were dissociated into single cells and cultured according to our FAILSAFE protocol (1 % FBS, 10 ng/ml bFGF and 20 ng/ml TGF α) that ensures robust long-term propagation of multipotent stem cells from the adult human brain .
Microarray analysis and public database mining
Microarray analysis was performed using Illumina gene chip. The data were quantile normalized and analyzed using J-Express (Molmine, Bergen, Norway) analysis software. We used published microarrays  with submission numbers GSE41470 (encompassing GSE41390, and GSE41394), GSE53800 (GSM1301030, GSM1301033, GSM1301042) as well as the microarrays submitted in connection with this work GSE60818 (encompassing GSE60705 and GSE60706). For more information on microarray analysis see Additional file 1: Figure S1 and Additional file 2: Figure S2 and Additional file 3 and Additional file 7.
REMBRANDT: Microarray data from the Repository for Molecular Brain Neoplasia Data (REMBRANDT, National Cancer Institute, 2005, http://rembrandt.nci.nih.gov) were accessed on January the 15th 2014. Hierarchical clustering with distance matrix and the K-means clustering of the NAT genes’ expression were performed using J-Express software (Molmine). For calculation of statistical parameters, we used Graphpad Prism (www.graphpad.com).
Western blot and quantification of protein expression
The cells were homogenized by triturating in Cell Extraction Buffer (Mammalian cell extraction kit, Biovision, Milpitas, CA, USA) and centrifuged through a QIAshredder (Qiagen, Venlo, Netherlands). For HIF1α and the western blots with phospho-proteins, the cells were homogenized in 10 mmol/L Tris–HCl (pH 7.4), 1 % SDS, 10 mmol/L NaF, and 2 mmol/L Na3VO4 as previously described . 20–40 μg of whole protein extract were mixed with loading buffer (NuPAGE, Invitrogen) and loaded onto a 4–12 % gradient Nu-PAGE gel (Invitrogen). Protein gels were blotted onto 0.45 μm PVDF membranes (Invitrogen). Membranes were blocked with 5 % skimmed milk in TBS/0.1 % Tween 20 (TBST) and probed with primary antibodies diluted in the same solution. Primary antibodies from Cell Signaling Technologies (CST, Danvers, MA, USA) were incubated in bovine serum albumin (BSA) according to recommended procedures. For a complete list of the antibodies used in this study see Additional file 12. For a complete list of the antibodies used in this study see Additional file 12. The blots were developed using the Lumiglo Reserve CL Substrate kit, and detected using the Epi Chemi II Darkroom (UVP-Laboratory Products, Upland, CA, USA). The intensities of the protein bands were quantified using Photoshop (Adobe, San Jose, CA, USA), background corrected and normalized to the intensities of the corresponding ACTB bands, so that the relative protein expression (RPE) could be calculated (RPEprotx = PEprotx/PEACT).
RNA isolation and real-time quantitative reverse-transcription PCR (qPCR)
Total RNA was isolated using an RNeasy Mini Kit and QIAshredder (Qiagen). For cDNA synthesis, experimental set up and oligonucleotide design we used the procedure previously described . For expression analysis of NAT12/NAA30 we used seven sets of oligonucleotides covering translated and non-translated regions of the reference sequence NM_001011713.2. For oligonucleotide sequence information see Additional file 3. qPCR was performed on an ABI PRISM 7900HT PCR machine (Applied Biosystems, Life Technologies, Foster City, Ca, USA) using SYBR Premix Ex Taq™ (Takara, Otsu, Japan) or Taqman probes (Applied Biosystems) according to the manufacturer's protocol. For analysis of the NAT12/NAA30 coding regions we used the following Taqman probes: ex2-3 (Hs02340852), ex3-4(Hs02340853) and ex4-5(Hs02340854). Crossing point (CP) values were generated using second-derivative calculation software (SDS2.2). Gene expression levels in GIC cultures were calculated using two house-keeping genes and multiple controls (values obtained for all tested NSCs and NFCs). For analysis of qPCR results we used 2-∆∆ CT-method and REST software .
RNAi-mediated gene knockdown
GIC culture T65 was used to establish stable NAT12/NAA30 knockdown cultures using shRNA1 (RHS4430-99149700; Clone ID V2LHS_180063), shRNA2 (RHS4430-99137454; Clone ID V2LHS_180058) and a non-silencing (NS) shRNA as a control (RHS4346) (all from Open Biosystems, Thermo Scientific, Huntsville, AL, USA). Cell cultures harboring shRNA constructs 1 and 2 are referred to as KD1 and KD2 respectively.
Production of virus in the 293 FT cell line using plasmid DNA and concentration of virus was done as following: Nine μg of plasmid DNA was transfected into the 293FT cell line using Arrest-In transfection reagent (according to the manufacturer’s protocol). Viral supernatants were collected after 48 and 72 h, centrifuged at 3000 rpm for 20 min at 4 °C and filtered through a sterile 0.45 μm low protein binding filter (Sarstedt, Nümbrecht, Germany). The virus was then concentrated in sterile SW28 ultracentrifuge tubes by ultracentrifugation (Beckman Optima™ LE-80 K ultracentrifuge, Fullerton, CA, USA) equipped with a SW-28 rotor at 23,000 rpm for 1.5 h at 4 °C. The pellet was resuspended in 200 μl of DMEM and aliquots of the concentrated virus were stored at −80 °C. GSC cultures (10×104 cells/ well) were then transduced in 24-well plates by adding 10 μl concentrated virus/well and the cells were incubated for 48 h at 37 °C in 5 % CO2. Three to five days after transduction, GFP positive cells were selected by FACS sorting using a FACS Diva cell sorter (Becton Dickinson, Franklin Lakes, NJ, USA) equipped with an argon ion laser, ‘TurboSort Plus’ option, and Diva software (Becton Dickinson). Alternatively, the cells were selected in 2 μg/ml Puromycin (Sigma-Aldrich, St. Louis, MO, USA) for 3–4 weeks before use for functional assays.
Cell viability, apoptosis and sphere-forming assay
For analysis of cell viability we used a colorimetric test based on tetrazolium salt – XTT (Roche Diagnostics, Indianapolis, IN, USA). Briefly, after sphere dissociation single cells were plated in neurosphere medium at a density of 1×104 cells per well in a flat-bottom 96-well plate (Sarstedt). Cells were incubated overnight and XTT-reagents added as recommended by the manufacturer. The colorimetric changes (absorbance) were measured after 24 h at 490 nm using a plate reader (Victor™ X5 Multilabel Plate Reader, Perkin Elmer 2030, Waltham, MA, USA). Five wells were evaluated for each cell culture. All results were presented as a mean of three independent experiments ± standard deviation. P-values were calculated using unpaired t-test with Welch’s correction.
Apoptosis was assessed using a fluorimetric assay which measures the activity of activated caspases (Roche Diagnostics). 1×104 cells were plated per well in a flat-bottom 96-well plate (Sarstedt) and incubated overnight before adding the reagents as recommended by the manufacturer. The fluorimetric changes were measured after 24 h at 485/535 nm using a plate reader (Victor™ X5 Multilabel Plate Reader). Five wells were evaluated for each cell culture.
The sphere-forming assay was done by seeding single cell suspensions containing 500 cells per well in ultra-low attachment flat bottom 96-well plates for 10 days (Sarstedt). Subsequently, the plates were imaged using GelCount™ (Oxford Optronics, Abington, UK). Only spheres >50 μm were taken into consideration and 10 wells were evaluated for each cell culture. The number and the size (average area) of spheres were measured using software supplied by the manufacturers. All results were presented as a mean of five independent experiments ± standard deviation. P-values were calculated using unpaired t-test with Welch’s correction.
Intracranial transplantation to SCID mice
All animal procedures were approved by the National Animal Research Authority. C.B.-17 SCID male mice (8–9 weeks old) were obtained from Taconic (Ejby, Denmark). They were acclimatized for > 1 week and divided into 3 groups: (i) T65 NS cells (n = 6), (ii) T65 KD1 cells (n = 6) and (iii) T65 KD2 cells (n = 6).
T65 GIC cultures were sorted by FACS after viral transduction and cultured for a few passages before transplantation. After tumor sphere dissociation, cells were counted and plated at a low density for overnight incubation. Prior to the inoculation of cells, mice were anesthetized and placed in a stereotactic frame (Kopf Instruments, Tujunga, CA, USA) and ~3×105 incubated cells were inoculated into the right striatum (AP 0 mm, RV 2 mm, ML 2 mm) using a Hamilton syringe and needle (Hamilton, Bonaduz, Switzerland) .
All animals were monitored for adverse effects and examined regularly for general appearance, signs of distress, and weight measurements taken. Mice were sacrificed when they lost >10 % of their bodyweight and developed neurological symptoms. Brains were harvested and further processed for HE-staining as described previously . Results for all groups are presented as median survival in days. Log-rank test was used to calculate the level of significance.
Mitochondrial depolarization assay
After passage, KD1, KD2 and NS cultures were seeded on fibronection-coated cover slips, and cultured in neurosphere medium in separate culture dishes for 1–6 days before experiments.
For each experiment, a cover slip was removed from its cell culture dish, and mounted in the microscopy chamber. Artificial CSF (aCSF) with the fluorescent cationic dye Rhodamine 123 (30 μM) (Invitrogen) was added to the chamber, and cells were loaded with dye for 15 min at 37 °C on the microscope. The chamber was covered with a lid, and subsequently perfused with air and 5 % CO2 for maintenance of appropriate pH. After loading of dye, the chamber was perfused with aCSF at a flow rate of approximately 1 ml/min. Imaging was started after at least 10 min of flow. Healthy looking cells were chosen by the criteria clear cell membrane and intact processes. Live imaging was performed using an Olympus IX81 inverted microscope equipped with a MT20 fluorescence light unit and Olympus xcellance software. Images were captured every 15 s. Rhodamine 123 (Rh123) used at 30 μM concentration accumulates in healthy mitochondria and is released to the cytosol upon depolarization of the mitochondria . This release is shown by increment in fluorescence from the cell. Since all cells expressed GFP, the red spectrum of Rh123 was recorded by using filters suitable for red fluorescent protein.
The cells were exposed to severe hypoxia for 15 min. Hypoxia was obtained by removing most of the oxygen by nitrogen perfusion of the flow solution. The oxygen scavenger Sodium Dithionite (0.75 mM) (Sigma-Aldrich) induced severe hypoxia. Control experiments were done with Sodium Dithionite perfusing the fluid with oxygen to control for any toxic effects other than caused by hypoxia. These control experiments did not show toxic effects on mitochondrial function. The protonophore Carbonyl cyanide 4-(trifluoromethoxy) phenylhydrazone (FCCP) (1 μM) (Sigma-Aldrich) was applied before the end of all experiments. FCCP fully depolarizes mitochondria, and thereby causes release of any remaining dye from the mitochondria. In this way any remaining mitochondrial membrane potential in the cell is revealed by an increase in Rh123 fluorescence.
GIC cultures were plated in 24-well plates (Sarstedt) or chamber slides (Sigma-Aldrich) pretreated with Retronectin 50 μg/ml (Takara) and incubated overnight to facilitate cell adhesion. Cells were fixed with 4 % PFA and washed with PBS. Tissue slides were fixed in 10 % buffered formalin and washed with PBS. Immunolabeling was performed as previously described . Briefly, permeabilization in 0.1 % Triton, was followed by blocking with 5 % BSA and 5 % blocking serum, and incubations firstly with primary antibody in 0.1 % Tween-20 (o/n) and then with secondary antibody (2 h). The protocol was slightly modified for the co-staining of NAT12/NAA30 and CD31 and 0.1 % Triton was replaced with 0.1 % (w/v) Saponin (Sigma-Aldrich). After staining of nuclei with Hoechst 33342 (Invitrogen) the slides were mounted with antifade reagent (ProLong Gold, Invitrogen).
For detection of NAT12/NAA30 protein we used the primary antibody #AV48508 [(Sigma-Aldrich), rabbit 1:200]. For confirmation of the obtained results we used #NBP1-70631 [(Novus Biologicals, Littleton, CO, USA), rabbit 1:200]. Immunofluorescence images shown in the article were obtained with #AV48508. We used primary antibodies against human nestin (NES) [#ab6320 (Abcam, Cambridge, UK), mouse 1:400], CD31 [#ab54211 (Abcam), mouse 1:100] and SOX2 [#AF2018 (R&D), goat 1:100].
Secondary antibodies were: Alexa Fluor 488 [#A-11059 (Invitrogen), donkey anti-mouse, 1:500], Alexa Fluor 555 [#A-31572 (Invitrogen), donkey anti-rabbit, 1:500], Alexa Fluor 594 [#R37119 (Invitrogen), donkey anti-rabbit, 1:500], and Alexa Fluor 594 [#A-11058 (Invitrogen), donkey anti-goat, 1:500].
Spheres from GIC culture T65 were dissociated into single cells, and incubated overnight to recover. Cells were pelleted, and incubated with primary conjugated antibodies against CD133/2 –PE conjugated (Miltenyi biotec) and SSEA1-PE (Becton-Dickinson). Cells were then analyzed using a LSRII flow cytometer (Becton-Dickinson).
Ostrom QT, Gittleman H, Liao P, Rouse C, Chen Y, Dowling J, et al. CBTRUS Statistical Report: Primary Brain and Central Nervous System Tumors Diagnosed in the United States in 2007–2011. Neuro Oncol. 2014;16 suppl 4:iv1–63.
Helseth R, Helseth E, Johannesen TB, Langberg CW, Lote K, Ronning P, et al. Overall survival, prognostic factors, and repeated surgery in a consecutive series of 516 patients with glioblastoma multiforme. Acta Neurol Scand. 2010;122(3):159–67.
Ronning PA, Helseth E, Meling TR, Johannesen TB. A population-based study on the effect of temozolomide in the treatment of glioblastoma multiforme. Neuro Oncol. 2012;14(9):1178–84.
Galli R, Binda E, Orfanelli U, Cipelletti B, Gritti A, De Vitis S, et al. Isolation and characterization of tumorigenic, stem-like neural precursors from human glioblastoma. Cancer Res. 2004;64(19):7011–21.
Singh SK, Clarke ID, Terasaki M, Bonn VE, Hawkins C, Squire J, et al. Identification of a cancer stem cell in human brain tumors. Cancer Res. 2003;63(18):5821–8.
Vescovi AL, Galli R, Reynolds BA. Brain tumour stem cells. Nat Rev Cancer. 2006;6(6):425–36.
Lee J, Kotliarova S, Kotliarov Y, Li A, Su Q, Donin NM, et al. Tumor stem cells derived from glioblastomas cultured in bFGF and EGF more closely mirror the phenotype and genotype of primary tumors than do serum-cultured cell lines. Cancer Cell. 2006;9(5):391–403.
Cheng L, Wu Q, Guryanova OA, Huang Z, Huang Q, Rich JN, et al. Elevated invasive potential of glioblastoma stem cells. Biochem Biophys Res Commun. 2011;406(4):643–8.
Molina JR, Hayashi Y, Stephens C, Georgescu MM. Invasive glioblastoma cells acquire stemness and increased Akt activation. Neoplasia. 2010;12(6):453–63.
Folkins C, Shaked Y, Man S, Tang T, Lee CR, Zhu Z, et al. Glioma tumor stem-like cells promote tumor angiogenesis and vasculogenesis via vascular endothelial growth factor and stromal-derived factor 1. Cancer Res. 2009;69(18):7243–51.
Das S, Marsden PA. Angiogenesis in glioblastoma. N Engl J Med. 2013;369(16):1561–3.
Bao S, Wu Q, Sathornsumetee S, Hao Y, Li Z, Hjelmeland AB, et al. Stem cell-like glioma cells promote tumor angiogenesis through vascular endothelial growth factor. Cancer Res. 2006;66(16):7843–8.
Bao S, Wu Q, McLendon RE, Hao Y, Shi Q, Hjelmeland AB, et al. Glioma stem cells promote radioresistance by preferential activation of the DNA damage response. Nature. 2006;444(7120):756–60.
Chen J, Li Y, Yu TS, McKay RM, Burns DK, Kernie SG, et al. A restricted cell population propagates glioblastoma growth after chemotherapy. Nature. 2012;488(7412):522–6.
Liu G, Yuan X, Zeng Z, Tunici P, Ng H, Abdulkadir IR, et al. Analysis of gene expression and chemoresistance of CD133+ cancer stem cells in glioblastoma. Mol Cancer. 2006;5:67.
Kalvik TV, Arnesen T. Protein N-terminal acetyltransferases in cancer. Oncogene. 2013;32(3):269–76.
Arnesen T, Van Damme P, Polevoda B, Helsens K, Evjenth R, Colaert N, et al. Proteomics analyses reveal the evolutionary conservation and divergence of N-terminal acetyltransferases from yeast and humans. Proc Natl Acad Sci U S A. 2009;106(20):8157–62.
Brown JL, Roberts WK. Evidence that approximately eighty per cent of the soluble proteins from Ehrlich ascites cells are Nalpha-acetylated. J Biol Chem. 1976;251(4):1009–14.
Arnesen T, Gromyko D, Pendino F, Ryningen A, Varhaug JE, Lillehaug JR. Induction of apoptosis in human cells by RNAi-mediated knockdown of hARD1 and NATH, components of the protein N-alpha-acetyltransferase complex. Oncogene. 2006;25(31):4350–60.
Starheim KK, Arnesen T, Gromyko D, Ryningen A, Varhaug JE, Lillehaug JR. Identification of the human N(alpha)-acetyltransferase complex B (hNatB): a complex important for cell-cycle progression. Biochem J. 2008;415(2):325–31.
Ametzazurra A, Larrea E, Civeira MP, Prieto J, Aldabe R. Implication of human N-alpha-acetyltransferase 5 in cellular proliferation and carcinogenesis. Oncogene. 2008;27(58):7296–306.
Polevoda B, Sherman F. Composition and function of the eukaryotic N-terminal acetyltransferase subunits. Biochem Biophys Res Commun. 2003;308(1):1–11.
Starheim KK, Gromyko D, Evjenth R, Ryningen A, Varhaug JE, Lillehaug JR, et al. Knockdown of human N alpha-terminal acetyltransferase complex C leads to p53-dependent apoptosis and aberrant human Arl8b localization. Mol Cell Biol. 2009;29(13):3569–81.
Wenzlau JM, Garl PJ, Simpson P, Stenmark KR, West J, Artinger KB, et al. Embryonic growth-associated protein is one subunit of a novel N-terminal acetyltransferase complex essential for embryonic vascular development. Circ Res. 2006;98(6):846–55. Epub 2006 Feb 16.
Pesaresi P, Gardner NA, Masiero S, Dietzmann A, Eichacker L, Wickner R, et al. Cytoplasmic N-terminal protein acetylation is required for efficient photosynthesis in Arabidopsis. Plant Cell. 2003;15(8):1817–32.
Polevoda B, Brown S, Cardillo TS, Rigby S, Sherman F. Yeast N(alpha)-terminal acetyltransferases are associated with ribosomes. J Cell Biochem. 2008;103(2):492–508.
Varghese M, Olstorn H, Sandberg C, Vik-Mo EO, Noordhuis P, Nister M, et al. A comparison between stem cells from the adult human brain and from brain tumors. Neurosurgery. 2008;63(6):1022–33. discussion 1033–4.
Vik-Mo EO, Sandberg C, Olstorn H, Varghese M, Brandal P, Ramm-Pettersen J, et al. Brain tumor stem cells maintain overall phenotype and tumorigenicity after in vitro culturing in serum-free conditions. Neuro Oncol. 2010;12(12):1220–30.
Heddleston JM, Li Z, Lathia JD, Bao S, Hjelmeland AB, Rich JN. Hypoxia inducible factors in cancer stem cells. Br J Cancer. 2010;102(5):789–95.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
da Huang W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.
Li Z, Bao S, Wu Q, Wang H, Eyler C, Sathornsumetee S, et al. Hypoxia-inducible factors regulate tumorigenic capacity of glioma stem cells. Cancer Cell. 2009;15(6):501–13.
Yang L, Lin C, Wang L, Guo H, Wang X. Hypoxia and hypoxia-inducible factors in glioblastoma multiforme progression and therapeutic implications. Exp Cell Res. 2012;318(19):2417–26.
Miyamoto S, Chiorini JA, Urcelay E, Safer B. Regulation of gene expression for translation initiation factor eIF-2 alpha: importance of the 3' untranslated region. Biochem J. 1996;315(Pt 3):791–8.
Edwalds-Gilbert G, Veraldi KL, Milcarek C. Alternative poly(A) site selection in complex transcription units: means to an end? Nucleic Acids Res. 1997;25(13):2547–61.
Lutz CS. Alternative polyadenylation: a twist on mRNA 3' end formation. ACS Chem Biol. 2008;3(10):609–17.
Amrani N, Ghosh S, Mangus DA, Jacobson A. Translation factors promote the formation of two states of the closed-loop mRNP. Nature. 2008;453(7199):1276–80.
Mayr C, Bartel DP. Widespread shortening of 3'UTRs by alternative cleavage and polyadenylation activates oncogenes in cancer cells. Cell. 2009;138(4):673–84.
Lembo A, Di Cunto F, Provero P. Shortening of 3'UTRs correlates with poor prognosis in breast and lung cancer. PLoS ONE. 2012;7(2), e31129.
Stewart M. Molecular mechanism of the nuclear protein import cycle. Nat Rev Mol Cell Biol. 2007;8(3):195–208.
Moore MJ. From birth to death: the complex lives of eukaryotic mRNAs. Science. 2005;309(5740):1514–8.
Song OK, Wang X, Waterborg JH, Sternglanz R. An Nalpha-acetyltransferase responsible for acetylation of the N-terminal residues of histones H4 and H2A. J Biol Chem. 2003;278(40):38109–12.
Matsuda Y, Ishiwata T, Yoshimura H, Hagio M, Arai T. Inhibition of nestin suppresses stem cell phenotype of glioblastomas through the alteration of post-translational modification of heat shock protein HSPA8/HSC71. Cancer Lett. 2015;357(2):602–11.
Dahlrot RH, Hansen S, Jensen SS, Schroder HD, Hjelmborg J, Kristensen BW. Clinical value of CD133 and nestin in patients with glioma: a population-based study. Int J Clin Exp Pathol. 2014;7(7):3739–51.
Laks DR, Masterman-Smith M, Visnyei K, Angenieux B, Orozco NM, Foran I, et al. Neurosphere formation is an independent predictor of clinical outcome in malignant glioma. Stem Cells. 2009;27(4):980–7.
Suva ML, Riggi N, Janiszewska M, Radovanovic I, Provero P, Stehle JC, et al. EZH2 is essential for glioblastoma cancer stem cell maintenance. Cancer Res. 2009;69(24):9211–8.
Starheim KK, Gevaert K, Arnesen T. Protein N-terminal acetyltransferases: when the start matters. Trends Biochem Sci. 2012;37(4):152–61.
Guertin DA, Sabatini DM. Defining the role of mTOR in cancer. Cancer Cell. 2007;12(1):9–22.
Choe G, Horvath S, Cloughesy TF, Crosby K, Seligson D, Palotie A, et al. Analysis of the phosphatidylinositol 3'-kinase signaling pathway in glioblastoma patients in vivo. Cancer Res. 2003;63(11):2742–6.
Rosner M, Siegel N, Valli A, Fuchs C, Hengstschlager M. mTOR phosphorylated at S2448 binds to raptor and rictor. Amino Acids. 2010;38(1):223–8.
Lund EL, Hog A, Olsen MW, Hansen LT, Engelholm SA, Kristjansen PE. Differential regulation of VEGF, HIF1alpha and angiopoietin-1, −2 and −4 by hypoxia and ionizing radiation in human glioblastoma. Int J Cancer. 2004;108(6):833–8.
Pistollato F, Abbadi S, Rampazzo E, Persano L, Della Puppa A, Frasson C, et al. Intratumoral hypoxic gradient drives stem cells distribution and MGMT expression in glioblastoma. Stem Cells. 2010;28(5):851–62.
Foyn H, Jones JE, Lewallen D, Narawane R, Varhaug JE, Thompson PR, et al. Design, Synthesis, and Kinetic Characterization of Protein N-Terminal Acetyltransferase Inhibitors. ACS Chem Biol, 2013.
Janiszewska M, Suva ML, Riggi N, Houtkooper RH, Auwerx J, Clement-Schatlo V, et al. Imp2 controls oxidative phosphorylation and is crucial for preserving glioblastoma cancer stem cells. Genes Dev. 2012;26(17):1926–44.
Joel M, Mughal AA, Grieg Z, Murrell W, Palmero S, Mikkelsen B, et al. Targeting PBK/TOPK decreases growth and survival of glioma initiating cells in vitro and attenuates tumor growth in vivo. Mol Cancer. 2015;14(1):121.
Vik-Mo EO, Nyakas M, Mikkelsen BV, Moe MC, Due-Tonnesen P, Suso EM, et al. Therapeutic vaccination against autologous cancer stem cells with mRNA-transfected dendritic cells in patients with glioblastoma. Cancer Immunol Immunother. 2013;62(9):1499–509.
Murrell W, Palmero E, Bianco J, Stangeland B, Joel M, Paulson L, et al. Expansion of multipotent stem cells from the adult human brain. PLoS ONE. 2013;8(8), e71334.
Moe IT, Pham TA, Hagelin EM, Ahmed MS, Attramadal H. CCN2 exerts direct cytoprotective actions in adult cardiac myocytes by activation of the PI3-kinase/Akt/GSK-3beta signaling pathway. J Cell Commun Signal. 2013;7(1):31–47.
Stangeland B, Rosenhave EM, Winge P, Berg A, Amundsen SS, Karabeg M, et al. AtMBD8 is involved in control of flowering time in the C24 ecotype of Arabidopsis thaliana. Physiol Plant. 2009;136(1):110–26.
Pfaffl MW, Horgan GW, Dempfle L. Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002;30(9), e36.
Duchen MR, Biscoe TJ. Relative mitochondrial membrane potential and [Ca2+]i in type I cells isolated from the rabbit carotid body. J Physiol. 1992;450:33–61.
Kunke D, Bryja V, Mygland L, Arenas E, Krauss S. Inhibition of canonical Wnt signaling promotes gliogenesis in P0-NSCs. Biochem Biophys Res Commun. 2009;386(4):628–33.
We would like to thank SFI-CAST and the Norwegian Research Council for supporting B. Stangeland and A. Mughal. We would also like to thank Emily Palmero, Sheryl Palmero, Ingunn Ramsnes and Grazyna Babinska/Trygve B. Leergaard for excellent technical help, and Prof. Ansgar Aasen and Dr. Sissel Reinlie for providing excellent working conditions.
The authors declare that they have no competing interests.
AAM participated in the design of the study, carried out cell culture work, immunolabeling, in vivo work, statistical analysis and drafted the manuscript. ZG participated in cell culture work and performed western blots. HS carried out the mitochondrial depolarization experiments and helped to draft the manuscript. AF participated in the work related to the transplantation experiments. ML participated in cell culture work and functional assays. MJ participated in the preparation of plasmids and viral transduction. MSA participated in planning and performing western blots. WM provided the NSC cultures and their microarray data and proofread the manuscript. EVM designed and performed flow cytometry, participated in discussion of individual experiments and helped to draft the manuscript. IAL initiated the study of stem cells, provided the biological material, had the general supervision of the lab, and participated in discussion of individual experiments and edited the manuscript. BS concieved the study, participated in the design of the study, established the shRNA technology, conducted qPCR, microarray and the bioinformatic analysis, participated in cell culture work and western blots and drafted the manuscript. All authors read and approved the final manuscript.
Additional file 1: Figure S1.
Microarray analysis of the selected NAT genes. A, Expression of NAA30/NAT12 in tissues and cells. Expression was measured by microarray analysis of 17 cell cultures (NSCs, NFCs and GICs) and four tissue specimens from normal and GBM tissues. Median expression in cell cultures was 126.8 while in tissues it was 94.29 with p value of 0.0013. B-H, Expression of the selected NATs in NSCs, GICs, normal brain tissues, GBM and NFCs. P values were calculated using the Mann Whitney test. Due to small sample sizes the p values proved reliable only for GICs vs. NSCs and were as follows: p(NAT12)=0.023, p(NAT6)=0.0002, p(NAT14)=0.78, p(NAT13)=0.3231, p(NAT8L)=0.0057, p(NAT8)=0.0006, p(NAT8B)=0.003 and p(NAT10)=0.4331.
Additional file 2: Figure S2.
K-means clustering of the expression data for the selected NAT genes.
Additional file 3: Supplementary File 1.
Oligonucleotides and microarray probes. (PDF 115 kb)
Additional file 4: Supplemetary File 2.
Additional file 5: Figure S3.
Kaplan-Meier survival plot for NAT genes (NAT12, NAT1, NAT2 and NAT10). P values (Log-rank test) signify the differences in patient survival in the selected groups with the following parameters: Up-Regulated: fold ≥ 2; Down-Regulated: fold ≥ 2. Reporter Type: Affymetrix.
Additional file 6: Figure S4.
Protein expression of NAT12/NAA30 detected with an anti-NAT12/NAA30 antibody (red) in biopsy specimens from normal brain (A), in GBM (B), in NSC (C-D) and GIC cultures (E-F). Cryostat sections with HE-staining from xenograft experiments showed that the control cultures (T65 and T65 NS) form visible tumors upon intracranial transplantation in contrast to the NAT12/NAA30 knockdowns (KD1 and KD2) (G). High magnification images showed that tumor cells are present in all groups even though the tumors are not visible in the knockdowns. Protein expression of NESTIN in GIC culture T08 (H) and SOX2 in T65 (I).
Additional file 7: Supplementary File 3.
Functional annotation based on the microarray analysis of NAT12/NAA30 knockdown cultures. (XLSX 201 kb)
Additional file 8: Figure S5-S7.
Figure S5. NAT12/NAA30 knockdown resulted in dysregulation of ribosome assembly as shown by microarray analysis. Using the DAVID functional annotation tool we analyzed the KEGG pathway “Ribosome” where 13 genes (p=5.6E-06) were differentially regulated in KD1 and KD2. See also Supplementary File 3. Dysregulated genes are marked with red asterisks. Figure S6. NAT12/NAA30 knockdown resulted in dysregulation of the p53 pathway as shown by microarray analysis. Using the DAVID functional annotation tool we analyzed the KEGG pathway “p53” where 8 genes (p=3.1E-03) were differentially regulated in KD1 and KD2. See also Supplementary File 3. Dysregulated genes are marked with red asterisks. Figure S7. NAT12/NAA30 knockdown resulted in dysregulation of sphingolipid metabolism as shown by microarray analysis. Using the DAVID functional annotation tool we analyzed the KEGG pathway “Sphingolipid metabolism” where 6 genes (p=4.8E-03) were differentially regulated in KD1 and KD2. See also Supplementary File 3.
Additional file 9: Figure S8.
Expression of the selected genes and proteins in the NAT12/NAA30 knockdown cultures. A, Western blot showing that the protein levels of active β-Catenin (ABC), phospho-p70 S6 kinase (Thr 389), phospho-STAT3 (Ser727) and GLI1 were slightly affected or not altered in the knockdown cultures. The signals obtained with Anti-SHH and Anti-phospho-p53 (Ser15) were below the detection level of western blot (not shown). For a complete list of ABs see Supplementary File 4. B, Western blots showing the expression of neural stem cell markers SOX2 and NESTIN, proteins representing hypoxia response (KCNMA1), cell cycle (CCND1) and ACTB (ACTIN). C, Western blots showing the expression of p-AKT (Ser473), total AKT, p-AKT (Ser308), p-p4E-BP1 and total p4E-BP1. D, Expression of selected stem cell genes SOX2, NANOG, POU5F1/OCT4, MYC and CCND1 on microarrays (n=4 for each cell culture). None of the selected genes were differentially regulated in the NAT12/NAA30 knockdown cultures. The bottom and top of each box indicate the 25th and 75th percentile (the lower and upper quartiles, respectively), and the band near the middle of the box represents the 50th percentile (the median). Interestingly, the knockdown strength was in correlation with the effect on downstream genes. In the KD2 culture that exhibited stronger knockdown effect at the protein level and better reduction in cell viability and number of spheres we could also detect decreased levels of GLI1, nestin and phospho proteins such as phospho-AKT (Ser308) and p-4E-BP1(Thr37/46) (Figure S8 A-C).
Additional file 10: Figure S9.
NAT12/NAA30 knockdown resulted in dysregulation of the hypoxia response pathway as shown by microarray analysis. Asterisks signify p values and indicate level of significance. *=(p≈0.01-0.05), **=(p≈0.001-0.01).
Additional file 11: Figure S10.
NAT12/NAA30 knockdown resulted in reduction of CD133+ cells as shown by flow cytometry.
Additional file 12: Supplementary File 4.
List of antibodies used for western blot. (DOCX 53 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Mughal, A.A., Grieg, Z., Skjellegrind, H. et al. Knockdown of NAT12/NAA30 reduces tumorigenic features of glioblastoma-initiating cells. Mol Cancer 14, 160 (2015). https://doi.org/10.1186/s12943-015-0432-z
- N-terminal acetyltransferases
- NAA30, NAT12
- NSCs and HIF1α