KAT6A, a Chromatin Modifier from the 8p11-p12 Amplicon is a Candidate Oncogene in Luminal Breast Cancer

The chromosome 8p11-p12 amplicon is present in 12% to 15% of breast cancers, resulting in an increase in copy number and expression of several chromatin modifiers in these tumors, including KAT6A. Previous analyses in SUM-52 breast cancer cells showed amplification and overexpression of KAT6A, and subsequent RNAi screening identified KAT6A as a potential driving oncogene. KAT6A is a histone acetyltransferase previously identified as a fusion partner with CREB binding protein in acute myeloid leukemia. Knockdown of KAT6A in SUM-52 cells, a luminal breast cancer cell line harboring the amplicon, resulted in reduced growth rate compared to non-silencing controls and profound loss of clonogenic capacity both in mono-layer and in soft agar. The normal cell line MCF10A, however, did not exhibit slower growth with knockdown of KAT6A. SUM-52 cells with KAT6A knockdown formed fewer mammospheres in culture compared to controls, suggesting a possible role for KAT6A in self-renewal. Previous data from our laboratory identified FGFR2 as a driving oncogene in SUM-52 cells. The colony forming efficiency of SUM-52 KAT6A knockdown cells in the presence of FGFR inhibition was significantly reduced compared to cells with KAT6A knockdown only. These data suggest that KAT6A may be a novel oncogene in breast cancers bearing the 8p11-p12 amplicon. While there are other putative oncogenes in the amplicon, the identification of KAT6A as a driving oncogene suggests that chromatin-modifying enzymes are a key class of oncogenes in these cancers, and play an important role in the selection of this amplicon in luminal B breast cancers.

Introduction A significant step in breast cancer progression is activation of oncogenes via gene amplification and overexpression [1]. The chromosome 8p11-p12 amplicon, containing approximately 55 genes, is present in 12% to 15% of breast cancers and is correlated with poor prognosis in primary breast tumors. Amplification of 8p11-p12 is also correlated with histologic grade, increased Ki-67 proliferation index, and decreased 5-year metastasis-free survival [2][3][4][5][6][7]. Due to its relevance in breast cancer, many studies have been aimed at characterizing this amplification and identifying the driving oncogenes in this region. Over the past several years, our laboratory and several others have analyzed the 8p11-p12 amplicon to identify possible driver oncogenes and determine the clinical relevance of these gene amplification events in breast cancer. These studies have resulted in the identification of a number of genes that play a role in breast cancer when the amplicon is present, including LSM-1, C8orf4 (TC-1), RAB11FIP1, WHSC1L1, ERLIN2, PROSC, PPAPDC, DDHD2 and others [2,6,[8][9][10][11][12][13]. Gelsi-Boyer et al. demonstrated that the amplicon can be sub-divided into four distinct regions that can be amplified independently of each www.neoplasia.com Volume 16 Number 8 August 2014 pp. 644-655 644 other, and this partly explains the large number of candidate oncogenes identified to date from this region [3]. Recently, our lab has taken a genome-scale shRNA screening strategy to identify additional driving oncogenes in a panel of SUM breast cancer cell lines, with and without 8p11-p12 amplification, for which copy number and expression data are available. KAT6A was found to be a consistent hit in the shRNA screen and is significantly amplified and overexpressed in SUM-52 cells.
KAT6A, a gene commonly found in the 8p11-p12 amplicon, was first identified in 1996 as part of a chromosomal translocation t(8;16)(p11;p13) with CREB-binding protein in a subtype of acute myeloid leukemia (AML) [14]. Since then, additional translocations involving KAT6A in AML have been discovered, including translocations with the HATs p300 and TIF2 [15]. The KAT6A-TIF2 fusion is capable of immortalizing myeloid progenitor cells and inducing acute myeloid leukemia upon injection into irradiated mice [16]. Characterization of KAT6A in these gene rearrangements provides compelling evidence that KAT6A is an oncogene in leukemogenesis, much like c-myc, which is well documented in various chromosomal translocations promoting several forms of leukemia, and is also amplified in breast cancer.
In this study, we provide evidence that KAT6A functions as a novel oncogene in luminal breast cancers bearing the 8p11-p12 amplicon. We identified KAT6A as a top shRNA screen hit among genes both amplified and overexpressed in the SUM-52 breast cancer cell line. The influence of KAT6A knockdown on proliferation and viability of SUM-52 cells, and other cell lines with the 8p11-p12 amplicon was validated in separate experiments involving several short hairpin RNAs (shRNAs). Finally, data obtained with primary human breast cancers coupled with TCGA data support a role for KAT6A as a breast cancer oncogene.

Reagents
All chemicals were purchased from Sigma Aldrich unless otherwise noted. Antibodies against Actin and anti-acetyl Histone 3 were purchased from Millipore. The total-Histone 3 antibody was purchased from Abcam, and the KAT6A antibody from Novus Biologicals. The FGFR inhibitor PD170374 was purchased from Selleckchem. Expression Arrest GIPZ lentiviral KAT6A shRNA constructs and packaging mix were purchased from Open Biosystems. Primers used for RT-PCR reactions were ordered from IDT. The SUM breast cancer cell lines were maintained in Hams F-12 cell culture medium as described previously [17][18][19].

Cellecta Screening Method
Virus pools expressing shRNA constructs were prepared according to the Cellecta Pooled Lentiviral shRNA Libraries User Manual protocol (www.cellecta.com). HEK 293 T cells were transfected with each of the three Cellecta library plasmid DNA pools (Human Modules 1-3) and the Cellecta Ready-to-Use Packaging Mix (Cat #CPCP-K2A). For each module, virus was titered and used to transduce 5 × 10 7 SUM-52 cells at an MOI of 0.5 in the presence of 5 μg/ml polybrene. Following transduction, cells were cultured for 3 days to allow expression of the resistance marker. Non-transduced cells were eliminated from the culture by the addition of 2 μg/ml puromycin to the growth media. Three days after the addition of puromycin, cells were trypsinized and one-half of the total population was harvested for genomic DNA preparation. This DNA served as the reference time point DNA. The remaining cells were plated and grown for approximately 7 population doublings before being harvested for genomic DNA preparation. Genomic DNA was prepared by phenol:chloroform extraction according to the Cellecta Pooled Lentiviral shRNA Libraries User Manual protocol.
See Supplementary Methods for detailed descriptions of the shRNA screen and statistical analyses.

RT-PCR
For RT-PCR reactions, total RNA was isolated using the Qiagen RNeasy Plus Mini Kit. RNA was converted into cDNA using the iScript Advanced cDNA Synthesis Kit (BioRad). Following cDNA conversion, RT-PCR reactions were performed in 96-well whitewalled plates using FastStart Universal SYBR Green Master mix (Roche). Reactions were performed on an Eppendorf Mastercycler Realplex2 under the following parameters: 1 cycle for 10 minutes at 95°C, 1 cycle for 15 seconds at 95°C, and 40 cycles at 53°C for 30 seconds. Primer sets included a control GAPDH primer set and a primer set specific to the KAT6A transcript.

Growth and Colony-Forming Assays
For growth assays, cells were plated in triplicate in 6-well plates at a density of 35,000 cells/well. Nuclei were harvested and cell counts obtained using a Beckman Coulter Counter as described previously [7].
For colony-forming assays, cells were plated in triplicate in 6-well plates at varying cell densities. Colonies formed for two weeks and then were fixed with 1 mL/well of 3.7% paraformaldehyde for 20 minutes at RT. Colonies were stained with 1 mL/well of 0.2% crystal violet for 15 minutes at RT and then de-stained with dH 2 O. For colony-forming assays including treatment with the FGFR inhibitor PD170374, cells were plated as stated above and treated with 0.1 μM of drug for 3 consecutive days. Colonies formed for an additional 10-14 days following drug treatment and were stained with crystal violet as described. Colony counts were generated using the GelCount System from Oxford Optronix.

Mammosphere Formation
Cells were plated in triplicate for primary mammosphere formation in low-attachment 6-well plates. Spheres formed for one week, at which point two wells were transferred to standard 6-well plates in media supplemented with FBS. Attached spheres were fixed with 1 mL/well of 3.7% paraformaldehyde for 20 minutes at RT and stained with 1 mL/well of 0.2% crystal violet for 15 minutes at RT. For secondary mammosphere formation, spheres from one well of the primary mammosphere assay were harvested and centrifuged at 2500 RPM for 5 minutes. Spheres were re-suspended in 1 mL of trypsin and incubated at RT for 5 minutes. The trypsin was diluted with media and centrifuged as before. The cell pellet was washed with 5 mL of PBS and centrifuged, and the final cell pellet was resuspended in 1 mL of 1 × Hanks buffer solution. Cells were re-plated in triplicate in 6-well low attachment plates. Spheres were allowed to form for 2 weeks and then transferred and stained with 0.2% crystal violet as described above. Images of primary and secondary spheres were taken using an EVOS Imaging System (Life Technologies). Mammosphere counts for both primary and secondary assays were generated using the GelCount System from Oxford Optronix. RNA sequencing on an Illumina HiScanSQ. RNA integrity was verified on an Agilent 2100 TapeStation (Agilent Technologies, Palo Alto, CA). 100-200 ng of total RNA was used to prepare RNA-Seq libraries using the TruSeq RNA Sample Prep kit following the protocol described by the manufacturer (Illumina, San Diego, CA). Sequencing was performed on an Illumina HiScanSQ. Samples were demultiplexed using CASAVA (Illumina, San Diego, CA). Fastq files were used to map reads to the human genome (hg19, UCSC) utilizing Tophat2 with default settings. Bam files were imported into Partek Genomics Suite (St. Louis, MO) for analysis. Samples were filtered for ≥ 10 reads. Differential gene expression was performed by comparison of biological duplicate samples in Genomics Suite using advanced analysis of variance (ANOVA). Differential gene express values were then filtered by ± 1.5 fold change and by unadjusted P ≤ .05.

Rationale for Studying the Oncogenic Potential of KAT6A
In an effort to comprehensively identify genes that control cell proliferation and survival in breast cancer cells, our laboratory is performing genome-scale, RNAi-based negative selection screens in our panel of SUM breast cancer cell lines. In this screen, the SUM-52 breast cancer cell line was transduced with a library of 82,000 lentiviral vectors that express shRNAs targeting 15,377 cellular genes with a minimum of 4 shRNAs per gene. Cells transduced with this library were harvested at either day 3 after infection or after approximately 7 population doublings, and the abundance of shRNAs at both time points was determined by PCR amplification, purification and sequencing of shRNA-associated barcode sequences. Fold depletion values were calculated by comparing the abundance of each shRNA at the final time point with the abundance at the initial time point. Fold depletion scores following statistical analysis for all of the shRNAs targeting each gene were used to identify genes that are necessary for the growth/survival of the SUM-52 breast cancer cell line. Integration of this list with the list of genes that are significantly copy number amplified in this cell line identified a small number of genes that are amplified, overexpressed and necessary for growth/survival of the cell line. KAT6A was found to be a consistent hit in this screen, as well as an additional screen using an early generation lower complexity library, and is significantly amplified and overexpressed in SUM-52 cells ( Figure 1A).
Based on the RNAi screen data, we performed microarray analysis to examine the expression level of KAT6A in a small panel of primary breast cancers available in our lab that potentially carry the 8p11-p12 amplicon. As seen in Figure 1B, several primary tumor samples from this panel exhibited high levels of KAT6A gene expression. Additionally, tumors with the highest expression of KAT6A also had a KAT6A gene amplification. Likewise, a series of SUM breast cancer cell lines were analyzed by microarray, and KAT6A expression was highest in cells harboring the 8p11-p12 amplicon, particularly SUM-52 cells ( Figure 1C). The SUM-52 cell line was generated from a patient with luminal B breast cancer who presented with an ER-positive primary tumor. The data in Figure 1D shows that KAT6A protein expression was highest in SUM-52 cells compared to normal MCF10A cells, SUM-149 and SUM-225 cell lines. Whereas SUM-225 cells have two small regions of 8p11 amplification, KAT6A is not amplified in these cells, and SUM-149 cells do not harbor the 8p11-p12 amplicon. As expected, KAT6A protein expression in SUM-225 cells was barely detectable; however, SUM-149 cells, although lacking the 8p11-p12 amplicon, express KAT6A protein at levels similar to SUM-52 cells. Because KAT6A is known as a histone acetyltransferase (HAT) with HAT activity on histone 3, we examined the acetylation status of histone 3 in cells with amplification of KAT6A. Total histone 3 acetylation was highest in SUM-52 cells compared to MCF10A, SUM-225, and SUM-149 cells ( Figure 1D); however, the difference in histone 3 acetylation was most interesting between SUM-52 and SUM-149 cells. Although both cell lines express comparable levels of KAT6A protein, histone 3 acetylation was only increased in SUM-52 cells.
To explore further the relevance of KAT6A in breast cancer, we examined TCGA data derived from 1,000 primary breast cancers, and found that KAT6A expression is altered in 22% of breast cancer cases. Additionally, there is a strong correlation between copy number and expression at the mRNA level ( Figure 1E). Kaplan-Meier analysis of more than 4000 breast cancer patients (http://kmplot.com/analysis/) showed that high level KAT6A expression (highest quartile) was associated with a statistically significant reduction in overall survival ( Figure 1F). Interestingly, the Kaplan-Meier curves diverge after approximately 100 months of follow-up, suggesting that KAT6A plays a role in late recurrences, an important clinical problem in luminal breast cancer. These data support the hypothesis that KAT6A functions as a driving oncogene in breast cancer, and warrant further investigation to validate its activity as a transforming oncogene.

KAT6A Knockdown Reduces Growth in 8p11-p12 Amplicon-Bearing Breast Cancer Cells
To validate the shRNA screen results, KAT6A knockdown studies were performed with several cell lines. Lentiviral shRNA vectors were prepared by transfection of 293FT cells with five shRNA plasmids targeting different KAT6A sequences. Initially, only SUM-52 and MCF10A cells were transduced with the lentiviral shRNAs. Following lentiviral infection and puromycin selection, RNA was harvested to confirm knockdown of KAT6A by RT-PCR. The results showed that two shRNAs yielded the best knockdown, sh1-566 (sh1) and sh3-895 (sh3), compared to a non-silencing shRNA sequence, (NS) (Supplementary Figure 1). To measure the effects of reduced KAT6A expression on the growth of the SUM-52 and SUM-225 cell lines, the cells were plated in monolayer in 6-well plates, and their growth was assessed by cell counts on days 1, 4, and 7. Knockdown of KAT6A at the message and protein level in SUM-52 cells (Figure 2A) resulted in 50% reduction in cell growth ( Figure 2C). KAT6A is not part of the 8p11-p12 amplicon in SUM-225 cells, and only a limited growth response to reduced KAT6A expression was observed ( Figure 2D). Infection of SUM-149 cells, which lack the 8p11-p12 amplicon but still express the protein, did not affect proliferation ( Figure 2E). Additionally, knockdown of KAT6A in MCF10A cells ( Figure 2B) showed no effect on proliferation ( Figure 2F). These data suggest that the effects of KAT6A on breast cancer cell proliferation are restricted to cells in which the gene is amplified and overexpressed.

KAT6A Reduces the Clonogenic Capacity of SUM-52 Cells
Since reduction of KAT6A expression yielded growth inhibition in SUM-52 cells, we next addressed the colony-forming efficiency of SUM-52 cells after knockdown of KAT6A. As shown in Figure 3A, SUM-52 cells with depleted KAT6A levels exhibited a dramatic reduction in colony-forming ability compared to control cells. Moreover, colonies that did grow after KAT6A knockdown were not only smaller in size, but exhibited an abnormal morphology ( Figure 3A and B). Cells with KAT6A knockdown produced small, scattered colonies composed of cells showing morphological signs of pyknosis. By contrast, non-silencing control cells produced tightlypacked colonies with cells displaying a normal appearance. Figure 3C shows the size distribution of colonies derived from KAT6A knockdown cells. The largest colonies from the knockdown cells were approximately 300 μm in diameter, whereas colonies from control cells were greater than 1000 μm in diameter. The Figure 1. Rationale for studying the oncogenic potential of KAT6A. A. Fold depletion scores from the genome-wide shRNA screen were used to generate probability density plots. The black plot was generated using depletion scores for all genes, the blue plot was generated using depletion scores for commonly essential genes, the turquoise plot was generated using depletion scores for the negative control gene set, and the green plot is a normal distribution whose mean was determined by the median of the fold depletion scores for all genes. The vertical red line indicates the 95th percentile of the normal distribution that was used as a cutoff for determining hits in the screen. colony-forming efficiency of MCF10A cells, however, was unaffected by knockdown of KAT6A ( Figure 3D and E). In a separate experiment, we examined the role of KAT6A knockdown in anchorageindependent colony formation. KAT6A knockdown also inhibited the soft agar colony-forming ability of SUM-52 cells, as KAT6A knockdown cells formed fewer colonies compared to the nonsilencing control ( Supplementary Figure 2A and B). To extend these observations to other 8p11-p12 harboring cancer cell lines in which KAT6A is amplified and overexpressed, the ovarian cancer cell lines OAW28 and 59 M were assayed for colony-forming efficiency following KAT6A knockdown (Supplementary Figure 3A and B). These additional cell lines were identified using the TCGA Bioportal and the Broad Institute's Cancer Cell Line Encyclopedia, which indicated a GISTIC copy number score of 1.563 for 59 M cells and  Neoplasia Vol. 16, No. 8, 2014 KAT6A is a Candidate Oncogene in Luminal Breast Cancer Turner-Ivey et al.
1.665 for OAW28 cells. As seen in Figure 3F and G, both OAW28 and 59 M cells, respectively, responded dramatically to knockdown of KAT6A, with both cell lines displaying a reduction in colony forming ability. Based on these results, KAT6A functions in regulating the clonogenic capacity of not only SUM-52 cells, but also of other cancer cell lines with the 8p11-p12 amplicon and KAT6A amplification/overexpression.

KAT6A Plays a Role in Breast Cancer Stem Cell Activity and Differentiation
In light of the effects of KAT6A knockdown on the growth and colony-forming efficiency of SUM-52 cells, the role of KAT6A in breast cancer stem cell activity and self-renewal was examined in a primary mammosphere assay. As shown in Figure 4A, approximately 2% of non-silencing control SUM-52 cells formed mammospheres. By contrast, less than 1% of KAT6A knockdown cells formed spheres.
To confirm and extend these results, a secondary mammosphere assay was performed. The non-silencing control cells formed 4-fold more mammospheres than in the primary assay, suggesting that the cells plated for secondary mammosphere assay contained a higher percentage of cancer-initiating stem cells. Consistent with results in the primary mammosphere assay, non-silencing control cells again formed more mammospheres than KAT6A knockdown cells in the secondary mammosphere assay ( Figure 4B). Images of mammospheres revealed striking differences not only between non-silencing and KAT6A knockdown spheres, but also between primary and secondary mammosphere assays. Mammospheres derived from control cells displayed classical sphere morphology. Although KAT6A knockdown spheres appeared similar, they were smaller in size compared to the control spheres ( Figure 4C). Allowing secondary mammosphere cultures to grow over an extended period revealed a remarkable difference in morphology between the non-silencing and KAT6A knockdown mammospheres. Mammospheres derived from non-silencing control cells began to form three-dimensional duct-like structures after two weeks in culture ( Figure 4D), which were never observed in mammosphere cultures derived from KAT6A knockdown cells. These results indicate that SUM-52 cells can form complex three dimensional structures after extended mammosphere culture, and that this capacity is lost upon KAT6A knockdown.
In many breast cancer cell types, expression of surface markers is a measure of stemness and self-renewal potential. In particular, the CD44 + /CD24 -/low expression pattern is associated with the breast cancer stem cell phenotype in some cell lines and tumors. Because knockdown of KAT6A resulted in reduced mammosphere formation, we investigated the expression levels of CD24 and CD44 in KAT6A knockdown and control cells by flow cytometry. Knockdown of KAT6A resulted in an increase in cells positive for both CD24 and CD44 to a level approximately three times above the non-silencing control. Perhaps more interesting was the increase in CD44 expression seen in KAT6A knockdown cells, as evidence now suggests that CD44 positivity is associated with differentiation [20]. Only 3.1% of non-silencing control cells were positive for CD44, Figure 4. KAT6A plays a role in breast cancer stem cell activity. A. SUM-52 control and KAT6A knockdown cells were plated in triplicate in 6-well low-attachment plates for a primary mammosphere assay. After one week, mammospheres were transferred to standard 6-well plates overnight, stained with 0.2% crystal violet and counted. B. Cells from the primary mammosphere assay were harvested, trypsinized, and re-seeded in a secondary mammosphere assay at 3,000 cells/well. After one week, spheres were processed as above for staining. C. Images of primary mammospheres following one week of growth. D. Images of secondary mammospheres following two weeks of growth.

Influence of KAT6A on Gene Expression in SUM-52 Cells
To gain insight into the role of KAT6A as a breast cancer oncogene, we performed RNA sequencing in SUM-52 KAT6A knockdown and control cells. The heat map in Figure 5A was generated from four replicate experiments using the two different shRNAs described above. Since the replicates from each shRNA clustered together, we pooled the data resulting in four data sets for the bioinformatic analysis of the effect of KAT6A knock-down on gene expression. As can be seen from the figure, KAT6A knockdown had a significant effect on gene expression. We used NIH DAVID to analyze the differentially expressed genes that resulted from KAT6A knockdown, and identified several biological processes, molecular functions, and KEGG pathways associated with this gene set. The biological processes affected by KAT6A knockdown are listed in Table 1, with genes involved in Response to Hormone Stimulus and Gland Development processes as top hits. From this analysis, several well-known genes involved in regulating growth were highlighted. Some of the significantly down-regulated genes in KAT6A knockdown cells included ESR1, IRS1, IGF1R, and ERBB4 (Table 2 and Supplementary Table 1, 2). Interestingly, not only was the adenylate cyclase gene ADCY7 significantly down-regulated in response to KAT6A knockdown, but so were its two regulatory subunits, PRKAR1B and PRKAR2B. A number of interesting transcription factors were also identified using NIH DAVID analysis. As shown in Table 3 (and Supplementary Table 3), GATA3, TBX3, and MYC were all significantly down-regulated upon KAT6A knockdown. KEGG pathway analysis yielded additional insight into the role of KAT6A, showing signaling pathways effected by knockdown of KAT6A. As seen in Table 5, some of the major pathways enriched following KAT6A knockdown included Wnt, PI3K-Akt, Insulin/IGF-1, and ErbB signaling. Finally, merging the list of significantly regulated genes identified by RNA sequencing in KAT6A knockdown cells with hits from the genome-wide shRNA screen in SUM-52 cells identified several overlapping genes, including IGF1R, IRS1, IGF2BP3, and GATA3 ( Table 5). The identification of genes that are both regulated in expression by KAT6A and hits in the shRNA screen provides functional evidence for the importance of these genes and pathways in breast cancer cells with KAT6A amplification. (See Table 4.)

KAT6A Interacts With FGFR2 to Enhance the Tumorigenicity of SUM-52 Cells
Previous work from our laboratory identified FGFR2 as a driving oncogene in SUM-52 cells. We have observed that while FGFR inhibition can inhibit the growth of these cells, this growth arrest is reversible, as washout of the drug results in prompt resumption of cell growth and no significant loss of colony forming ability (unpublished   data and Figure 6A). To address the potential interaction of FGFR2 and KAT6A in SUM-52 cells, the colony-forming ability of SUM-52 non-silencing control cells or SUM-52 KAT6A knockdown cells was examined following treatment with the FGFR inhibitor PD170374. As shown in Figure 6A, exposure of non-silencing control cells to FGFR inhibition did not significantly affect the colony-forming ability of the cells after washout of the drug. As observed previously, the colony-forming efficiency of SUM-52 KAT6A knockdown cells was reduced compared to controls. However, unlike the control cells, the colony-forming efficiency of the KAT6A knockdown cells was decreased further following 72 hours of exposure to PD173074, and this decrease was statistically significant (P b .001). Images of colonies with KAT6A knockdown revealed a striking morphological difference compared to control colonies. As seen in Figure 6B, non-silencing control colonies appeared healthy in size and morphology after treatment and washout of PD170374; however, the KAT6A knockdown colonies that did grow after drug treatment were much smaller in size, with cells showing signs of pyknosis. MCF10A cells, however, were not negatively affected by either KAT6A knockdown or PD170374 treatment (seen in Figure 6C). These results suggest that the KAT6A oncogene interacts with the FGFR2 oncogene, a clear driver of transformed phenotypes in SUM-52 cells, and has an influence on the response of cells to a targeted receptor tyrosine kinase inhibitor.

Discussion
Elucidation of novel candidate oncogenes from the 8p11-p12 amplicon has been the subject of on-going studies in our laboratory and a number of other laboratories over many years. The 8p11-p12 amplicon occurs in approximately 15% of breast cancers, the majority of which are luminal cancers. An interesting feature of this amplicon is the large number of candidate oncogenes already identified in this the region, and the consistent finding that many of these genes are coordinately amplified and overexpressed in breast cancer and lung cancer [2][3][4]8,9,12,13,21]. In the present study, we demonstrate that KAT6A is another candidate oncogene from this region. Gelsi-Boyer et al. provided compelling evidence that the 8p11-p12 amplicon can be sub-divided into four distinct regions, each of which can be amplified independently of the others [3]. We have reported previously on candidate oncogenes from region III of this amplicon, which contains genes such as WHSC1L1, DDHD2, and FGFR1 [5][6][7]10,11,22]. The KAT6A oncogene is in region IV of the amplicon, which is closest to the centromere. KAT6A is co-amplified with region III oncogenes in approximately 7% of cases and amplified independently of this region in approximately 3% of breast cancer cases according to TCGA data. We found KAT6A to be a consistent hit in an RNAi-based negative selection screen in SUM-52 breast cancer cells. Since KAT6A is both amplified and overexpressed in SUM-52 cells, the identification of KAT6A as a positive screen hit was suggestive of its importance as an oncogene in this cell line. Validation experiments using several KAT6A-specific shRNAs confirmed the screen result and demonstrated a significantly reduced growth rate and reduction in clonogenic potential in knockdown cells compared to non-silencing controls. Similar results were obtained using ovarian cancer cell lines that have KAT6A amplification, but normal mammary epithelial cells, or other breast cancer cells without the 8p11-p12 amplicon were not affected by KAT6A knockdown. Analysis of greater than 1000 primary breast tumors in the TCGA data base showed that KAT6A is amplified in 8.7% of invasive breast cancers, and this amplification is associated with overexpression at the mRNA level. Furthermore, overexpression of KAT6A was found to be associated with late recurrence of breast cancer. Thus, both functional genomics data obtained in our lab and data from the TCGA point to KAT6A as an important breast cancer oncogene.
Cancer cells are heterogeneous, and within any population of cancer cells in a cell line, only a relatively small fraction are considered clonogens based on their ability to form continuously expanding colonies when grown at clonal density, either in monolayer or in soft-agar. In the work reported here, we demonstrate a greater than 10-fold decrease in the clonogenic capacity of SUM-52 cells following knockdown of KAT6A. Colony size distribution shows that the majority of KAT6A knockdown colonies were smaller than the non-silencing control colonies. We interpret these results to indicate that reduction of KAT6A levels not only reduced the clonogenic fraction of SUM52 cells, but also reduced the growth rate of colonies that did form from the remaining clonogens.
We were intrigued by the relationship between the dramatic decrease in colony and mammosphere formation in SUM-52 cells upon knockdown of KAT6A, and the genes regulated by KAT6A. IGF1R, IRS1, and ERBB4 were among the genes reduced in expression as a result of KAT6A knockdown, which is consistent with a decrease in proliferation and colony formation by these cells. Interestingly, a recent publication by Heiser et al. showing the efficacy of 77 compounds in more than 50 breast cancer cell lines showed that the SUM-52 cell line was among the most sensitive to IGF1R inhibition [23], a result consistent with our association between KAT6A overexpression and the IGF1R.
Several interesting transcription factors were also down-regulated following knockdown of KAT6A, including MYC, TBX3, and P value indicates the probability of obtaining the number of screen hits in the pathway by chance alone. FDR indicates the global probability corrected for false discovery rate of obtaining the indicated number of pathway hits by chance alone. GATA3. It is known that gene rearrangements involving KAT6A that occur in leukemia results in an inability of cells to differentiate, which is in keeping with the known roles of KAT6A and other MYST family genes in regulation of stem cell self-renewal and differentiation in hematopoiesis [24][25][26]. Thus, one mechanism by which KAT6A may affect breast tumorigenesis is by altering self-renewal capacity and differentiation of cells of the luminal mammary lineage. Several of our results are consistent with this hypothesis. Mammosphere formation in SUM-52 KAT6A knockdown cells was not only significantly reduced, but so too was the ability of the cells to form large, branching structures when cultured for extended periods in low-attachment plates in stem cell medium. The down-regulation of key transcription factors associated with mammary gland development in KAT6A knockdown cells is consistent with this observation. C-MYC is a well-known breast cancer oncogene and mechanisms that drive MYC expression in breast cancer cells play a key role in the growth and survival of breast cancer cells [27][28][29]. GATA3 is a key transcription factor in mammary gland development and plays a role in determining luminal cell fate [30,31], and as we and others have shown, KAT6A and other oncogenes of the 8p11 amplicon are over expressed in luminal breast cancers like the SUM-52 cells. Finally, the finding of KAT6A regulation of TBX3 provides a link between breast cancer biology and leukemogenesis induced by this oncogene. KAT6A regulates expression of TBX1 in hematopoietic progenitor cells [26], while TBX3 has been shown to play an essential role in mammary gland development and ductal branching. Furthermore, TBX3 has been shown to be over expressed in a fraction of breast cancers, notably estrogen receptor positive, luminal breast cancers [32]. TBX3 has also been found to be an estrogen responsive gene in MCF-7 cells, and our results show that KAT6A regulates expression of ESR1 in SUM-52 cells. Thus, the genes regulated by KAT6A in SUM-52 cells are all key genes in the development of breast cancers of the luminal lineage, which further implicates KAT6A as a key driver of luminal breast cancer. KAT6A, also known as monocytic leukemia zinc finger protein (MOZ) or MYST3, belongs to the highly conserved MYST family of histone acetyltransferases, important in post-translational modification of histones and chromatin structure. Currently, the MYST family includes Tip60, MOF, HBO1, MOZ (MYST3), and MORF [22]. The MYST family is known to acetylate lysine residues on histones H4, H3, H2A and H2A variants, and recent evidence suggests that MYST proteins regulate themselves by acetylation. Specifically, KAT6A has acetyltransferase activity towards lysine residues on histones H2B, H3 (K14) and H4 (K5, K8, K12, and K16) [22,26]. Acetylation of histone H3K9 by KAT6A has been shown in vivo [33]. Although KAT6A is widely-known as an acetyltransferase, some HAT-independent functions have been reported. For example, KAT6A interacts with RUNX1 [34,35], leading to transcription of hematopoietic genes. KAT6A also interacts with RUNX2 and plays a role in osteogenesis [36,37] and these functions are independent of HAT activity. KAT6A also interacts with PU.1, important in the development of lymphoid and myeloid lineages, and with MLL, a histone methyl transferase that induces di-and tri-methylation of Histone H3 on lysine 4 [26]. ASH2L, another amplified and over expressed gene from the 8p11 amplicon has been shown to be the histone methylase in the MLL complex. The H3K4 di-and tri-methyl marks are associated with active gene expression, and KAT6A-MLL interactions have been associated with expression of homeobox genes in hematopoietic progenitor cells. Interestingly, two homeobox genes, HOXA11 and HOXA13 were regulated by KAT6A in SUM-52 cells. Thus, in breast cancer cells with the 8p11 amplicon, KAT6A may cooperate with other epigenome modifying oncogenes to aberrantly drive expression of a number of important cancer genes.
The significance of KAT6A in breast cancer is perhaps best reflected in the Kaplan-Meier plot shown in Figure 1F, which demonstrates a significant reduction in overall survival of patients with ER positive breast cancer and KAT6A amplification and overexpression. Not only is there a statistically significant divergence in mortality among patients with KAT6A overexpression, but this affect occurs following approximately 100 months of survival. This suggests that KAT6A plays a role in late recurrence of breast cancer, an issue of unique importance in ER-positive breast cancer. The results reported here suggest that amplification/overexpression of KAT6A could be a useful biomarker for predicting the risk of late recurrence, which could aid in the clinical decisions regarding the duration of anti-estrogen therapy. It is well-accepted that cancer is a disease driven by genetic abnormalities; however, it is becoming clear that alterations in epigenetic pathways play a role in oncogenesis. The identification of KAT6A as a driving oncogene in breast cancers with 8p11-p12 amplification is consistent with the emerging field of cancer epigenetics. Future studies will explore the role of KAT6A in breast cancer self-renewal and mechanisms governing clonogenicity. While there are other putative oncogenes in the amplicon, KAT6A has emerged as an oncogene that plays an important role in the selection of this amplicon in luminal breast cancers.