Enhanced MYC association with the NuA4 histone acetyltransferase complex mediated by the adenovirus E1A N-terminal domain activates a subset of MYC target genes highly expressed in cancer cells

ABSTRACT

The proto-oncogene MYC is a transcription factor over-expressed in many cancers and required for cell survival. Its function is regulated by histone acetyltransferase (HAT) complexes, such as the GCN5 complex and the NuA4/Tip60 complex. However, the roles of the HAT complexes during MYC function in cancer have not been well characterized. We recently showed that adenovirus E1A and its N-terminal 80 aa region, E1A 1-80, interact with the NuA4 complex, through the E1A TRRAP-targeting (ET) domain, and enhance MYC association with the NuA4 complex. We show here that the ET domain mainly targets the MYC-NuA4 complex. By global gene expression analysis using E1A 1-80 and deletion mutants, we have identified a panel of genes activated by targeting the MYC-NuA4 complex and notably enriched for genes involved in ribosome biogenesis and gene expression. A second panel of genes is activated by E1A 1-80 targeting of both the MYC-NuA4 complex and p300, and is enriched for genes involved in DNA replication and cell cycle processes. Both panels of genes are highly expressed in cancer cells. Since the ET domain is essential for E1A-mediated cellular transformation, our results suggest that MYC and the NuA4 complex function cooperatively in cell transformation and cancer.

INTRODUCTION

Adenovirus E1A 243R (E1A) is a viral oncoprotein that transforms cells by targeting several key cellular factors including p300, TRRAP, and Rb [1, 2]. E1A targets TRRAP via an N-terminal domain; however, the functional consequence of E1A-TRRAP association remains unclear. E1A 1-80, the E1A N-terminal 80 aa region, represses transcription of some viral and cellular promoters in vivo and in vitro [3-6]. E1A 1-80 targets p300 and TRRAP [7]. By proteomic and co-immunoprecipitation analyses, we recently showed that E1A 1-80 and the full-length E1A interact with the NuA4 complex and enhance its association with MYC [7]. This activity requires the E1A N-terminal TRRAP-targeting (ET) domain [7]. The ET domain is capable of enhancing MYC association with the NuA4 complex when fused to the N-terminus of MYC [8]. Since association with TRRAP is essential for transformation by E1A and by MYC [9-11], we hypothesize that MYC association with the NuA4 complex may be increased during cellular proliferation and transformation induced by MYC and by E1A. The proto-oncogene MYC is over-expressed in many cancers and required by essentially all cancers for survival [12-14]. MYC over-expression promotes tumorigenesis [14-16], and transient down-regulation of MYC below a certain level causes tumor regression in animal models [17]. Histone acetyltransferase (HAT) complexes, such as the NuA4 complex and the GCN5 complex, participate in transcriptional activation by MYC [18-21]. However, how the HAT complexes participate in MYC function during cancer development remains uncharacterized. The human NuA4 complex, the largest HAT complex containing up to 20 subunits [22], has a core HAT enzyme Tip60, and is involved in chromatin remodeling, gene activation, and DNA damage repair [22-24]. Proteomic analysis of MYC suggests that MYC associates with the NuA4 complex in mouse embryonic stem cells [21]. MYC and some components of the NuA4 complex are found to co-occupy a set of promoters and co-activate those genes (“MYC-module” genes) in mouse embryonic stem cells to help define stem cell characteristics [21]. Significant over-expression of a similar set of genes in human cancers seems to be correlated with tumor invasiveness, suggesting the importance of MYC association with the NuA4 complex in cancer [21].

In this report, we performed RNA-seq analysis with human foreskin fibroblast cells (HS68) expressing E1A 1-80 and mutants lacking a functional ET domain or p300-targeting domain, and identified two panels of genes that are activated by E1A 1-80 targeting of the MYC-NuA4 complex or by targeting of both the MYC-NuA4 complex and p300. Expression of both panels of genes is increased in three cancer cell lines examined compared to HS68 cells, suggesting that expression of both panels of genes in cancer may involve MYC association with the NuA4 complex.

RESULTS AND DISCUSSION

Activation of MYC target genes by E1A 1-80 correlates with enhanced MYC association with the NuA4 complex

We have previously shown that both the full length E1A and its N-terminal 80 aa (E1A 1-80) enhance MYC association with the NuA4 complex [7]. Further, we have shown that E1A 1-80 with deletion of aa 2-11, which inactivates p300 targeting remains capable of interaction with the NuA4 complex, whereas deletion of aa 26-35, which inactivates the E1A TRRAP-targeting (ET) domain abolishes E1A 1-80 interaction with the NuA4 complex (for a summary, see Figure 1A). These data suggest that the ET domain interacts with TRRAP in the NuA4 complex to promote its association with MYC (Figure 1B). Because E1A 1-80 enhances MYC association with the NuA4 complex through the ET domain, it is reasonable to expect that TRRAP-targeting helps assemble a MYC-NuA4 complex that activates a subset of MYC target genes. To examine this possibility, we designed a strategy for examination of MYC target gene activation using contact-inhibited human fibroblast cells (HS68). In culture, HS68 cells become contact-inhibited after forming a confluent monolayer, possibly resulting in down-regulation of genes important for cell proliferation. However, they resume growth and proliferation after re-plating and have been used to examine MYC regulation of gene expression [8].

HS68 cells were infected with Ad viral expression vectors for E1A 1-80FH (Flag-HA tag on the C-terminus of E1A 1-80), and two deletion mutants, ∆2-11 and ∆26-35 which do not interact with p300 or the MYC-NuA4 complex, respectively (Figure 1A) [7]. To examine potential function of E1A 1-80 through MYC, RT-qPCR analysis was performed with gene-specific primers for three selected genes that are known to be activated by MYC: MCT1 (monocarboxylate transporter 1), NCL (nucleolin), and NOP2 (nucleolar protein 2 homologue) [8]. As shown in Figure 1C, both E1A 1-80FH and ∆2-11 activated all three genes, while ∆26-35 was defective. Thus, activation of these MYC target genes is correlated with a functional ET domain which promotes MYC association with the NuA4 complex (Figure 1A). Consistent with previous results [25], expression of E1A 1-80FH or the two deletion mutants repressed cellular MYC mRNA (Figure 1D - MYC). Since E1A 1-80FH repressed MYC mRNA more strongly than either ∆2-11 or ∆26-35, it is possible that targeting of the MYC-NuA4 complex and targeting of p300 are both involved in this repression. However, statistical analysis suggested that ∆2-11 was significantly more effective than ∆26-35 in repression of MYC. Both E1A 1-80FH and ∆2-11 significantly activated TP53 (Figure 1D - TP53), consistent with the observation that over-activation of MYC usually leads to TP53 activation [26, 27]. In contrast, ∆26-35 did not activate TP53, possibly because it lacks a functional ET domain to enhance MYC association with the NuA4 complex. Analysis of the mRNA for E1A 1-80FH and its two deletion mutants showed that all forms of E1A 1-80FH are expressed to comparable levels (Figure 1D - E1A).

The ET domain targets mainly the NuA4 complex

Since the ET domain interacts with TRRAP that is present in at least three HAT complexes [22], it remains unclear as to which HAT complex the ET domain interacts with. Therefore, we constructed an ET-MYC fusion protein that efficiently associates with components of the NuA4 complex by co-IP and Western blot analysis [8]. Since the NuA4 complex contains up to 20 subunits, we examined ET-MYC interaction with the NuA4 complex more thoroughly by a proteomic approach. Lentivirus vectors were used to express i) FH-MYC to detect interaction with the NuA4 complex, ii) ET-MYC to detect enhanced association with the NuA4 complex, and iii) Flag-Tip60 to detect components of the NuA4 complex as control. Cell lysates were immunoprecipitated with Flag antibody beads, and bound protein complexes eluted with Flag peptide and subjected to liquid chromatography-tandem MS analysis and associated proteins identified. As shown in Table 1, this analysis shows that ET-MYC associates with 14 known components of the NuA4 complex, while FH-MYC associates with only three. ET-MYC associated proteins overlap with those associated with Tip60 (control for the NuA4 complex). Thus, the ET-MYC complex closely resembles a MYC-NuA4 complex. Importantly, from this analysis and our previous proteomic analysis [7], the ET domain appears to associate only with the NuA4 complex, and not the GCN5 or PCAF complex that also contains TRRAP [22], suggesting a key role for the NuA4 complex in the function of the ET domain.

E1A 1-80 modulates global gene expression by targeting the MYC-NuA4 complex and p300 separately and cooperatively

E1A 1-80 contains the domains of the full-length E1A that target the HAT enzyme p300 and the MYC-NuA4 HAT complex (Figure 1A). To identify genes activated by E1A 1-80 targeting of the these two HAT activities, RNA-seq analysis was performed with mRNAs from HS68 cells infected with adenovirus vectors expressing control lacZ, E1A 1-80FH, and the two deletion mutants, ∆2-11 and ∆26-35. Importantly, the ∆2-11 mutant remains capable of interaction with TRRAP and therefore the NuA4 complex (see Figure 1A) in the absence of potential interference from p300-targeting, and thus may help reveal genes activated by the enhanced MYC association with the NuA4 complex. The ∆26-35, on the other hand, may regulate gene expression by interaction with p300 as one of the major cellular targets for E1A.

RNA-seq analysis revealed more than 2000 genes activated over 35% by E1A 1-80FH (with a P value ≤ 0.05). Three panels of genes are derived from these E1A 1-80-activated genes (Figure 2A): i) those that overlap with genes activated by ∆2-11 and require targeting of the MYC-NuA4 complex (MYC-NuA4 panel, or MNA4 panel); ii) those that overlap with genes activated by ∆26-35 and require targeting of p300 (p300 panel), and iii) genes that are only weakly activated by both ∆2-11 and ∆26-35 and therefore require targeting of both the MYC-NuA4 complex and p300 for activation (MYC-NuA4 and p300 panel, or MNP300 panel). Since E1A 1-80 region has been shown to interact with factors in addition to p300 and the NuA4 complex (which contains TRRAP, p400 and other subunits) [2], it remains to be determined if any of these additional factors contributes to the activation of these >2000 identified genes.

To understand the function of these different panels of genes, PANTHER analysis for gene enrichment in biological processes was performed (see Figure 2A legend). This analysis revealed that the MNA4 panel genes (440 genes total) are strongly enriched for genes involved in gene expression and ribosome biogenesis (Figure 2A and Figure 3), while the p300 panel genes (129 genes total) do not enrich for any identifiable biological processes. The MNP300 panel genes (352 genes total) are most significantly enriched in genes involved in DNA replication and cell cycle (Figure 2A and Figure 3).

To confirm activation of the identified genes by E1A 1-80, selected MNA4 panel genes involved in ribosome biogenesis/gene expression and MNP300 panel genes involved in DNA replication/cell cycle were examined by RT-qPCR. As shown (Figure 2B), the selected MNA4 panel genes involved in ribosome biogenesis/gene expression are activated efficiently by both E1A 1-80FH and ∆2-11 but not by ∆26-35. These results confirm that E1A 1-80 activates the MNA4 panel genes through the ET domain by targeting the NuA4 complex and enhancing its association with MYC (Figure 1B). The selected MNP300 panel genes involved in DNA replication/cell cycle are activated by E1A 1-80FH much more efficiently than by ∆2-11 or by ∆26-35 (Figure 2C), confirming that activation of the MNP300 panel genes involves targeting of both the MYC-NuA4 complex and p300.

We next analyzed the average activation of all the MNA4 panel genes involved in ribosome biogenesis/gene expression (157 total out of the MNA4 panel), utilizing the data from the RNA-seq analysis. This analysis showed that both E1A 1-80FH and ∆2-11 activate this group of genes on average by > 60%, whereas ∆26-35 is defective in activation (Figure 2D). In contrast, the MNP300 panel genes involved in DNA replication/cell cycle (82 total out of the MNP300 panel) were activated by E1A 1-80FH (by ~ 65%), but poorly by both ∆2-11 and ∆26-35 (by < 10%, Figure 2E), confirming that both the MYC-NuA4 complex and p300 are involved in activation of the MNP300 panel genes.

The MNA4 panel genes are authentic MYC target genes in HS68 cells

The 440 MNA4 panel genes (Figure 2A) were compared with the 1037 genes activated by MYC by > 20% under similar cell culture conditions [8] and 276 (63%) are found to overlap, suggesting that enhanced MYC association with the NuA4 complex promoted by E1A 1-80FH (and ∆2-11) activates a large portion of authentic MYC target genes. Other genes of the MNA4 panel may be authentic MYC target genes with lower levels of activation by MYC. In contrast, the 352 MNP300 panel genes (Figure 2A) have 21 (6%) overlapping with the same 1037 MYC activated genes, consistent with the possibility that activation of the MNP300 panel genes requires targeting of both the MYC-NuA4 complex and p300, i.e., targeting of either the MYC-NuA4 complex or p300 is insufficient for activation. Importantly, out of the 54 ribosome biogenesis genes from the MNA4 panel, 50 (93%) are activated by MYC in HS68 cells [8], consistent with their activation by MYC association with the NuA4 complex.

The MNA4 panel and MNP300 panel genes are highly expressed in cancer cells and candidates for important function in cancer

To examine the potential roles of the MYC-NuA4 complex in cancer, we first compared the expression of the MNA4 panel and MNP300 panel genes in three cancer cell lines, HeLa (cervical cancer), MB231 (breast cancer), and U2OS (osteosarcoma), relative to HS68 cells. RNA was isolated and RT-qPCR performed for i) selected MNA4 panel genes involved in ribosome biogenesis/gene expression and ii) selected MNP300 panel genes involved in DNA replication/cell cycle. As shown in Figure 4, most of these selected genes are highly expressed in the three cancer cell lines to significant degrees relative to HS68 cells, with the MNP300 panel genes (Figure 4B) over-expressed at higher levels than the MNA4 panel genes (Figure 4A). Importantly, MYC is also highly expressed in the three cancer cell lines (Figure 4A). In contrast, Tip60 and p300, both with important HAT activities, are not over-expressed in the cancer cell lines (Figure 4A). PolB from the MNP300 panel was also not over-expressed in the cancer cell lines (Figure 4B). The high expression of the genes from the MNA4 panel and MNP300 panel in cancer cell lines examined are suggestive of important roles in cancer.

Since activation of the MNP300 panel genes in HS68 cells by E1A 1-80FH requires both the MYC-NuA4 complex and p300 (Figure 2), it is likely that the MYC-NuA4 complex normally cooperates with other transcription factors for efficient activation of these genes. The lack of a strong activation of these genes by MYC in HS68 cells [8] suggests that such cooperation may not exist in normal cells. However, the high expression of the MNP300 panel genes in three cancer cell lines examined (Figure 4A) suggest that the MYC-NuA4 complex and other transcription factors may cooperate in cancer cells.

The MNA4 panel of genes is most significantly enriched for genes involved in ribosome biogenesis/gene expression. The same conclusion was drawn from global gene expression analysis with the ET-MYC fusion protein [8], which identified an E1A-activated MYC-NuA4 panel (EMNA panel) of genes (262 genes total). Comparing the MNA4 panel from this study with the EMNA panel reveals 78 genes in common which are also significantly enriched for genes involved in ribosome biogenesis (unpublished results). Genes involved in ribosome biogenesis from the MNA4 panel may increase a basic cellular function to accommodate the need of cancer cells for increased protein synthesis. Bioinformatic search for consensus MYC binding motif CACGTG has revealed predominantly ribosome biogenesis genes in human and fly genomes [28]. Importantly, increased ribosome biogenesis function has been suggested to directly contribute to cancer development [29, 30]. Further, mice with haploinsufficient MYC have reduced ribosome biogenesis/protein translation, and a longer lifespan [31], suggesting the essential role of MYC for ribosome biogenesis. Thus, it would appear that MYC association with the NuA4 complex and its activation of genes involved in ribosome biogenesis may contribute to cancer development.

In summary, our results show that genes up-regulated by the MYC-NuA4 complex alone and together with p300 are highly expressed in cancer cells and are potential candidates for important functions in cancer.

MATERIALS AND METHODS

Expression constructs and proteomic analysis

Adenoviral expression constructs for E1A 1-80FH, ∆2-11 and ∆26-35, and purification of adenoviruses are described earlier [7]. Lentiviral expression constructs for FH-MYC, ET-MYC, and Flag-Tip60 are described earlier [7, 8]. Generation of lentiviruses and proteomic analyses are under conditions described [7].

Cell culture and viral infection

Human HS68, HeLa, MB231, and U2OS are from ATCC, and cultured in DMEM (Life Technologies) supplemented with 10% fetal bovine serum and 50 U/ml of Pen/Strep. For RNA-seq analysis, HS68 cells were cultured to confluence and maintained for three days. Cells were re-plated at 1.2 x 10⁶ cells/T75 flask, and infected with 40 PFU/cell of the adenoviral expression vector in triplicates for 18 h. RNA was prepared as described previously [25] .

RNA-seq analysis

RNA samples in triplicates were subjected to PolyA selection and RNA-seq analysis at the Washington University Genome Technology Access Center as described earlier [8].

PANTHER gene enrichment analysis

PANTHER (version 12) gene enrichment analysis [32-34] was performed online (http://www.geneontology.org). Biological processes with Bonferroni corrected (http://www.geneontology.org) P values of smaller than a set value were shown in Figure 3.

ACKNOWLEDGMENTS

We thank Michael Zhao for help with managing RNA-seq files, and the Genome Technology Access Center in the Department of Genetics at Washington University School of Medicine for help with RNA-seq analysis. The Center is partially supported by NCI Cancer Center Support Grant #P30 CA91842 to the Siteman Cancer Center and by ICTS/CTSA Grant# UL1TR000448. This publication is solely the responsibility of the authors and does not necessarily represent the official view of NCRR or NIH.

CONFLICTS OF INTEREST

The authors declare that they have no conflicts of interest with the contents of this article.

FUNDING

MG is supported by Lifetime Research Career Award (RCA #AI-04739) from the NIH. Saint Louis University Cancer Research Fund has provided support for the project.

In Memorium

Dr. Maurice Green (b. May 5, 1926) passed away on Tuesday, December 5 – only a few days after notification that this manuscript had been accepted for publication in Genes & Cancer. This is the 308th publication originating from his rich and productive six-decade career. Maurice’s contributions to virology, molecular virology, and tumor biology are significant, and too extensive to even briefly detail here. Interested readers are referred to his Wikipedia entry at https://en.wikipedia.org/wiki/Maurice_Green_(virologist). He will be deeply missed as a mentor, colleague, and friend.

• 1. Berk AJ. Recent lessons in gene expression. cell cycle control, and cell biology from adenovirus. Oncogene. 2005; 24(52): 7673-85. [PubMed]
• 2. Peter Pelka, Jailal NG Ablack, Gregory J Fonseca, Ahmed F Yousef, and Joe S Mymryk. Intrinsic structural disorder in adenovirus E1A: a viral molecular hub linking multiple diverse processes. J Virol. 2008; 82(15): 7252-63. [PubMed] https://doi.org/10.1128/JVI.00104-08.
• 3. Green M, Panesar NK, and Loewenstein PM. The transcription-repression domain of the adenovirus E1A oncoprotein targets p300 at the promoter. Oncogene. 2008; 27(32): 4446-55. [PubMed]
• 4. Song CZ, Loewenstein PM, and Green M. Repression in vitro, by human adenovirus E1A protein domains, of basal or Tat-activated transcription of the human immunodeficiency virus type 1 long terminal repeat. J Virol. 1995; 69(5): 2907-11. [PubMed]
• 5. Song CZ, Loewenstein PM, Toth PM, Tang Q, Nishikawa A, and Green M. The adenovirus E1A repression domain disrupts the interaction between the TATA binding protein and the TATA box in a manner reversible by TFIIB. Mol Cell Biol. 1997; 17(4): 2186-93. [PubMed]
• 6. Loewenstein PM, Wu SY, Chiang CM, and Green M. The adenovirus E1A N-terminal repression domain represses transcription from a chromatin template in vitro. Virology. 2012; 428(1): 70-5. [PubMed] https://doi.org/10.1016/j.virol.2012.03.021.
• 7. Zhao LJ, Loewenstein PM, and Green M. Ad E1A 243R oncoprotein promotes association of proto-oncogene product MYC with the NuA4/Tip60 complex via the E1A N-terminal repression domain. Virology. 2016; 499: 178184. [PubMed] https://doi.org/10.1016/j.virol.2012.03.021.
• 8. Zhao LJ, Loewenstein PM, and Green M. Adenovirus E1A TRRAP-targeting domain-mediated enhancement of MYC association with the NuA4 complex activates a panel of MYC target genes enriched for gene expression and ribosome biogenesis. Virology. 2017; 512: 172-179. [PubMed] https://doi.org/10.1016/j.virol.2017.08.010.
• 9. Deleu L, Shellard S, Alevizopoulos K, Amati B, and Land H. Recruitment of TRRAP required for oncogenic transformation by E1A. Oncogene. 2001; 20(57): 8270-5. [PubMed]
• 10. Park J, Kunjibettu S, McMahon SB, and Cole MD.The ATM-related domain of TRRAP is required for histone acetyltransferase recruitment and Myc-dependent oncogenesis. Genes Dev. 2001; 15(13): 1619-24.. 2001; 20(57): 8270-5. [PubMed] https://doi.org/10.1101/gad.900101.
• 11. McMahon SB, Van Buskirk HA, Dugan KA, Copeland TD, and Cole MD. The novel ATM-related protein TRRAP is an essential cofactor for the c-Myc and E2F oncoproteins. Cell. 1998; 94(3): 363-74. [PubMed]
• 12. Tansey WP. Mammalian MYC Proteins and Cancer. New Journal of Science. 2014; 2014: 1-27. [PubMed] https://doi.org/10.1158/2159-8290.CD-15-0507.
• 13. Dang CV. MYC on the path to cancer. Cell. 2012; 149(1): 22-35. [PubMed] https://doi.org/10.1016/j.cell.2012.03.003.
• 14. Meyer N and LZ Penn. Reflecting on 25 years with MYC. Nat Rev Cancer. 2008; 8(12): 976-90. [PubMed]
• 15. Shachaf CM, Gentles AJ, Elchuri S, Sahoo D, Soen Y, Sharpe O, Perez OD, Chang M, Mitchel D, Robinson WH, Dill D, Nolan GP, Plevritis SK, et al. Genomic and proteomic analysis reveals a threshold level of MYC required for tumor maintenance. Cancer Res. 2008; 68(13): 5132-42. [PubMed] https://doi.org/10.1158/0008-5472.CAN-07-6192.
• 16. Murphy DJ, Junttila MR, Pouyet L, Karnezis A, Shchors K, Bui DA, Brown-Swigart L, Johnson L, and Evan GI. Distinct thresholds govern Myc’s biological output in vivo. Cancer Cell. 2008; 14(6): 447-57. [PubMed] https://doi.org/10.1016/j.ccr.2008.10.018.
• 17. Jain M, Arvanitis C, Chu K, Dewey W, Leonhardt E, Trinh M, Sundberg CD, Bishop JM, and Felsher DW. Sustained loss of a neoplastic phenotype by brief inactivation of MYC. Science. 2002; 297(5578): 102-4. [PubMed]
• 18. Liu X, Tesfai J, Evrard YA, Dent SY, and Martinez E. c-Myc transformation domain recruits the human STAGA complex and requires TRRAP and GCN5 acetylase activity for transcription activation. J Biol Chem. 2003; 278(22): 20405-12. [PubMed] https://doi.org/10.1074/jbc.M211795200.
• 19. McMahon SB, Wood MA, and Cole MD. The essential cofactor TRRAP recruits the histone acetyltransferase hGCN5 to c-Myc. Mol Cell Biol. 2000; 20(2): 556-62. [PubMed]
• 20. Frank SR, Parisi T, Taubert S, Fernandez P, Fuchs M, Chan HM, Livingston DM, and Amati B. MYC recruits the TIP60 histone acetyltransferase complex to chromatin. EMBO Rep. 2003; 4(6): 575-80. [PubMed] https://doi.org/10.1038/sj.embor.embor861.
• 21. Kim J, Woo AJ, Chu J, Snow JW, Fujiwara Y, Kim CG, Cantor AB, and Orkin SH. A Myc network accounts for similarities between embryonic stem and cancer cell transcription programs. Cell. 2010; 143(2): 313-24. [PubMed] https://doi.org/10.1016/j.cell.2010.09.010.
• 22. Lee KK and Workman JL. Histone acetyltransferase complexes: one size doesn’t fit all. Nat Rev Mol Cell Biol. 2007; 8(4): 284-95. [PubMed]
• 23. Doyon Y, Selleck W, Lane WS, Tan S, and Cote J. Structural and functional conservation of the NuA4 histone acetyltransferase complex from yeast to humans. Mol Cell Biol. 2004; 24(5): 1884-96. [PubMed] https://doi.org/10.1128/MCB.24.5.1884-1896.2004.
• 24. Squatrito M, Gorrini C, and Amati B. Tip60 in DNA damage response and growth control: many tricks in one HAT. Trends Cell Biol. 2006; 16(9): 433-42. [PubMed]
• 25. Zhao LJ, Loewenstein PM, and Green M. The adenoviral E1A N-terminal domain represses MYC transcription in human cancer cells by targeting both p300 and TRRAP and inhibiting MYC promoter acetylation of H3K18 and H4K16. Genes Cancer. 2016; 7(3-4): 98-109. [PubMed] https://doi.org/10.18632/genesandcancer.99.
• 26. Lindstrom MS and Wiman KG. Myc and E2F1 induce p53 through p14ARF-independent mechanisms in human fibroblasts. Oncogene. 2003; 22(32): 4993-5005. [PubMed]
• 27. Hong S, Pusapati RV, Powers JT, and Johnson DG. Oncogenes and the DNA damage response: Myc and E2F1 engage the ATM signaling pathway to activate p53 and induce apoptosis. Cell Cycle. 2006; 5(8): 801-3. [PubMed]
• 28. Brown SJ, Cole MD, and Erives AJ. Evolution of the holozoan ribosome biogenesis regulon. BMC Genomics. 2008; 9: 442. [PubMed] https://doi.org/10.1186/1471-2164-9-442.
• 29. Derenzini M, Montanaro L, and Trere D. Ribosome biogenesis and cancer. Acta Histochem. 2017; 119(3): 190197. [PubMed] https://doi.org/10.1186/1471-2164-9-442.
• 30. Montanaro L, Trere D, and Derenzini M. Nucleolus, ribosomes, and cancer. Am J Pathol. 2008; 173(2): 301-10. [PubMed] https://doi.org/10.2353/ajpath.2008.070752.
• 31. Hofmann JW, Zhao X, De Cecco M, Peterson AL, Pagliaroli L, Manivannan J, Hubbard GB, Ikeno Y, Zhang Y, Feng B, Li X, Serre T, Qi W, et al. Reduced expression of MYC increases longevity and enhances healthspan. Cell. 2015; 160(3): 477-88. [PubMed] https://doi.org/10.1016/j.cell.2014.12.016.
• 32. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000; 25(1): 25-9. [PubMed] https://doi.org/10.1038/75556.
• 33. Mi H, Muruganujan A, Casagrande JT, and Thomas PD. Large-scale gene function analysis with the PANTHER classification system. Nat Protoc. 2013; 8(8): 1551-66. [PubMed] https://doi.org/10.1038/nprot.2013.092.
• 34. Mi H, Huang X, Muruganujan A, Tang H, Mills C, Kang D, and Thomas PD. PANTHER version 11: expanded annotation data from Gene Ontology and Reactome pathways, and data analysis tool enhancements. Nucleic Acids Res. 2017; 45(D1): D183-D189. [PubMed] https://doi.org/10.1093/nar/gkw1138.