The epigenetic modifier JMJD6 is amplified in mammary tumors and cooperates with c-Myc to enhance cellular transformation, tumor progression, and metastasis
© Aprelikova et al. 2016
Received: 28 January 2016
Accepted: 6 April 2016
Published: 14 April 2016
Oncogene overexpression in primary cells often triggers the induction of a cellular safeguard response promoting senescence or apoptosis. Secondary cooperating genetic events are generally required for oncogene-induced tumorigenesis to overcome these biologic obstacles. We employed comparative genomic hybridization for eight genetically engineered mouse models of mammary cancer to identify loci that might harbor genes that enhance oncogene-induced tumorigenesis.
Unlike many other mammary tumor models, the MMTV-Myc tumors displayed few copy number variants except for amplification of distal mouse chromosome 11 in 80 % of the tumors (syntenic to human 17q23-qter often amplified in human breast cancer). Analyses of candidate genes located in this region identified JMJD6 as an epigenetic regulatory gene that cooperates with Myc to enhance tumorigenesis. It suppresses Myc-induced apoptosis under varying stress conditions through inhibition of p19ARF messenger RNA (mRNA) and protein, leading to reduced levels of p53. JMJD6 binds to the p19ARF promoter and exerts its inhibitory function through demethylation of H4R3me2a. JMJD6 overexpression in MMTV-Myc cell lines increases tumor burden, induces EMT, and greatly enhances tumor metastasis. Importantly, we demonstrate that co-expression of high levels of JMJD6 and Myc is associated with poor prognosis for human ER+ breast cancer patients.
A novel epigenetic mechanism has been identified for how JMJD6 cooperates with Myc during oncogenic transformation. Combined high expression of Myc and JMJD6 confers a more aggressive phenotype in mouse and human tumors. Given the pleiotropic pro-tumorigenic activities of JMJD6, it may be useful as a prognostic factor and a therapeutic target for Myc-driven mammary tumorigenesis.
Tumorigenesis is a multistep process involving the accumulation of genetic aberrations. Initiating oncogenic alterations promote the selection of additional genetic and epigenetic changes that favor transformation, tumor growth, and metastases. Initiating oncogenic alterations in untransformed cells in most, if not all cases, trigger cellular safeguard mechanisms that induce cellular senescence (e.g., Ras) or apoptosis (e.g., Myc, E2F1, or E1A). Therefore, oncogene-induced tumorigenesis generally requires cooperating genetic events to overcome these safeguard mechanisms.
It has been firmly established that many genetically engineered mouse models (GEMs) of mammary cancer are valuable systems to dissect tumorigenic pathways that may involve multiple genetic aberrations. Many models have been designed to mimic human breast cancer by either overexpression of known oncogenes (Myc, Ras, Wnt, PyMT, Erbb2) or deletion of tumor suppressor genes (BRCA1/2, p53, and Rb) [1, 2]. For example, several mouse models of breast cancer, including MMTV-PyMT or WAP-Myc, express markers associated with human luminal-type breast cancer . Murine models (C3(1)-Tag or BRCA1 deficiency together with p53 mutation) allow tumor development with a characteristic basal-like phenotype [1, 3–5]. On the other hand, some other models show mixed features and high levels of heterogeneity . Unlike human breast cancers, many mouse models of mammary cancer are based upon the induction of a specific oncogenic event through the overexpression of relevant oncogenes or inactivation of tumor suppressor genes, leading to the evolution of oncogene-specific secondary pathways.
A predominant mechanism leading to additional genetic alterations required for progression of tumorigenesis appears to be related to changes in genome copy number variants (CNV). A recent study demonstrated that 22 % of the haploid genome in breast cancer is affected by chromosome rearrangements , thus indicating that CNVs are the major contributor to the accumulation of additional genetic changes during tumor progression.
Array comparative genomic hybridization (CGH) has been a powerful tool to identify chromosomal regions that may harbor amplified oncogenes or deleted tumor suppressor genes. These techniques combined with gene expression analysis, revealed a significant correlation between human and multiple GEM models of mammary gland tumors [3, 7, 8]. Recently 662 regions of chromosomal aberrations conserved between human and mouse breast cancers were identified . These studies allow not only the identification of novel drivers of tumorigenesis but also find supporting genetic alterations necessary for promoting cellular transformation and tumor progression.
We used array CGH to identify CNVs in eight genetically engineered mouse models of mammary cancer to find recurrent CNVs that potentially could harbor critical genes enhancing oncogene-induced transformation. We determined that the distal region of chromosome 11 was amplified in 80 % of MMTV-Myc-driven mouse mammary gland tumors. This locus is syntenic to human chromosome 17q23-qter, a region that is often amplified in human breast cancers. We hypothesized that gene(s) located in this locus are critical for Myc-induced tumor development and progression. Myc is commonly amplified in many cancers of different origins. Importantly, Myc is overexpressed in 25–30 % of all breast cancer cases [9–12].
Further functional analyses of several candidate genes overexpressed in the mouse 11q locus revealed that the epigenetic modifier JMJD6 is able to inhibit Myc-induced apoptosis, which is critical for tumor progression. While Myc-triggered cell death may involve several pathways, the predominant responder to aberrant induction of Myc in primary cells is p19ARF and upregulation of p53. We found that JMJD6 represses p19ARF, at least in part, by demethylation of Arg3 of histone H4 associated within the p19ARF promoter. Therefore, JMJD6 amplification may cooperate with Myc to enhance neoplastic transformation of primary epithelial cells. Additionally, when overexpressed in MMTV-Myc-driven tumor cells that lack the 11q amplification, JMJD6 induces EMT, increases cell migration and invasion in vitro, and stimulates tumor growth in vivo. Most importantly, JMJD6 dramatically increased lung metastatic colonization of these otherwise non-metastatic cells. Bioinformatics analyses of human breast cancer tumors revealed a significant decrease in survival of patients with ER+ tumors when Myc and JMJD6 were highly expressed together as compared to high Myc expression alone.
Identification of JMJD6 as a gene that cooperates with Myc to enhance tumorigenesis could provide a novel therapeutic target for breast and other cancers where Myc is an essential driver of tumorigenesis, since to date no successful therapies directly targeting Myc have been developed.
Identification of a chromosome 11 amplicon as the major genomic alteration in MMTV-Myc mammary tumors
Identification of genes overexpressed in the chromosome 11 amplicon in MMTV-Myc mammary tumors
To identify potential candidate cell lines derived from MMTV-Myc tumors for further in vitro studies, we performed CGH analysis of several cell lines and found that the Myc83 cell line harbors the chromosome 11 amplicon, as previously observed , while the 88CT1 cell line has relatively few CNVs without amplification of the chromosome 11 locus (Additional file 1: Figure S2A). We also found increased expression levels of selected genes (JMJD6, Tnrc6c, and Ube2o) by RT-qPCR in the Myc83 cells compared to 88CT1 cells, consistent with the status of the chromosome 11 amplification (Additional file 1: Figure S2B).
Identification of JMJD6 as a gene that suppresses Myc-induced apoptosis
Since Myc expression increases apoptosis in primary cells, which is a major response preventing full transformation of cells by Myc alone, we tested the hypothesis that the increased expression of candidate genes would suppress Myc-induced cell death, whereas depletion of any of the seven selected genes would increase Myc-induced cell death in vitro. Myc-induced apoptosis in many cases requires an intact p53 pathway. Sequence analysis confirmed the wild-type status of p53 in both Myc83 and 88CT1 cells. However, etoposide treatment of these cell lines revealed that Myc83 responded to etoposide with a robust increase in p53 and p21 protein levels, while 88CT1 had a much lower expression of p53 (possibly because of the amplification of MDM2 in these cells, as determined by CGH analysis). Therefore, for the primary apoptosis screen, we chose Myc83 cells that express high constitutive levels of c-Myc and contain the chromosome 11 amplification. The limitation of this model is that one cannot conclude that cell death is Myc-dependent or whether the tested genes have a more general effect on cell viability.
The anti-apoptosis effect of JMJD6 is dependent upon JMJD6 enzymatic activity
JMJD6 is an iron- and 2-oxoglutarate-dependent dioxygenase with the capability to hydroxylate lysine residues in histone and non-histone proteins [26, 29, 30]. Lysine hydroxylation of RNA splicing factors results in production of differentially spliced pre-messenger RNA (mRNA) molecules , while hydroxylation of histone proteins may result in changes in transcriptional regulation of targeted gene expression . In order to determine whether the enzymatic activity of JMJD6 is necessary for the inhibition of Myc-induced cell death, we mutated His187 to Ala in the iron-coordinating center of the enzyme that has previously been showed to inactivate JMJD6 and demonstrated that mutant JMJD6 does not suppress cell death (Fig. 4e). Western blotting confirmed that the levels of mutated protein expression were similar to the expression of the wild-type JMJD6 (Fig. 4f).
In order to understand mechanisms by which JMJD6 inhibits Myc-induced cell death, we first analyzed its effect on Myc protein levels and cellular localization. Additional file 1: Figure S5A shows that MycER™ or endogenous levels of c-Myc were not compromised by overexpression of JMJD6. Also, in the presence of 4-hydroxytamoxifen, ectopically expressed MycER™ re-localizes to the cell nucleus similar to control cells without JMJD6 overexpression (Additional file 1: Figure S6B).
Identification of mechanisms contributing to JMJD6 inhibition of Myc-induced apoptosis
Our finding that JMJD6 may reduce p19ARF stimulation by Myc led us to determine whether the blunting effect of JMJD6 on Myc-induced cell death is p19ARF dependent. We overexpressed p19ARF protein in cells with JMJD6 and the inducible MycERTM and determined levels of cell death after starving cells in serum-free conditions. Indeed, increased levels of p19ARF protein rescued Myc-induced apoptosis in cells with high JMJD6 expression, demonstrating that p19ARF suppression by JMJD6 contributes to the alleviation of cell death (Fig. 6c, d).
To further confirm that high expression of JMJD6 would correlate with decreased expression of p19ARF in vivo, we analyzed expression of these two genes in mouse MMTV-Myc mammary gland tumors. This experiment also showed robust negative correlation between JMJD6 and p19ARF (Additional file 1: Figure S7).
JMJD6 reduces p19ARF expression through histone H4 modifications on the p19ARF promoter region
JMJD6 is involved in post-translational modifications of non-histone and histone proteins [29, 30, 32, 36]. However, when JMJD6 regulates histone modifications, it modifies mRNA levels of the target gene. To discriminate which of the pathways is involved in p19ARF silencing by JMJD6, we studied p19ARF transcript levels in the presence of high JMJD6 expression. We clearly observed that JMJD6 overexpression significantly reduces mRNA levels of p19ARF (Fig. 6e), while three different JMJD6 shRNA constructs significantly increased p19ARF mRNA when expressed in NMuMG cells (Fig. 6f).
Since endogenous mouse JMJD6 is poorly recognized by the antibody available for ChIP studies, we decided to test human HeLa cells to analyze whether endogenous JMJD6 in these cells is found on the ARF promoter. We tested the human p14ARF promoter with two primer sets using a human JMJD6 antibody for ChIP and found higher levels of JMJD6 associated with the ARF promoter compared to immunoprecipitation with a non-specific IgG (Fig. 7d). In summary, we found that in non-tumorigenic mouse epithelial cells, JMJD6 may collaborate with c-Myc to initiate tumor formation by suppressing Myc-induced cell death by inhibiting transcriptional expression of the p19ARF tumor suppressor protein.
JMJD6 promotes Myc-induced tumor formation
Interestingly, we demonstrated that 88CT1 cells established from MMTV mice have a compromised p53 pathway, possibly resulting from amplification of MDM2 (gain of 8 copies) as determined by array CGH. We found that the reduction in the percentage of apoptotic cells in 88CT1 cells overexpressing JMJD6 was associated with elevated levels of the Bcl2 anti-apoptotic family members Bcl-xl and Bcl-w (Fig. 8e) that were previously shown to inhibit Myc-induced cell death [41–43].
JMJD6 enhances cell migration, invasion, and metastases
MMTV-Myc mammary gland tumors usually produce well-differentiated tumors with a relatively long latency and relatively few metastases [12, 44, 45], implying that additional genetic alterations are required for a more aggressive phenotype. However, several observations raise the possibility that elevated expression of Myc protein may inhibit cellular migration, invasion, and metastasis formation in mouse xenograft models of breast cancer . We therefore investigated whether JMJD6 might enhance migration, invasion, and metastases of MMTV-Myc mammary tumor cells.
High expression of JMJD6 in human breast tumors is associated with a worse prognosis of Myc-high tumors
Our array CNV results of 8 GEM models of mammary cancer demonstrated that there are common genomic amplifications and losses in multiple models as well as alterations that are more pronounced in specific models. These data, consistent with prior reports, demonstrates that the evolution of genomic alterations during tumorigenesis in these models appears to depend upon the genetic event that initiates transformation. We focused on the MMTV-Myc model since Myc is overexpressed in about 30 % of human breast cancers and this model displays relatively few CNVs except for a prominent amplification of mouse distal chromosome 11.
We identified several candidate genes within this locus that could potentially cooperate with Myc and enhance tumorigenesis. Functional analyses determined that the epigenetic modifier JMJD6 cooperates with Myc to augment oncogenic transformation of normal mammary epithelial cells as well as to increase the aggressiveness and metastatic potential of Myc-driven tumors. Increased expression of JMJD6 resulted from the in vivo selection for the amplification of the distal region of chromosome 11 in Myc-overexpressing tumors. The chromosome 11 amplified locus is not exclusively specific for Myc-driven tumors but also has been found in mouse mammary gland tumors with overexpression of Erbb2, PyMT, or BRCA1 mutations (this study and [8, 13–16, 23]), which implicates JMJD6 as potentially augmenting tumorigenesis in other oncogenic pathways.
Myc activation in normal cells is not capable of promoting neoplastic transformation because of the overwhelming induction of cell death triggered by Myc protein expression above normal physiological levels [47–51]. The mechanisms of Myc-induced apoptosis have been studied very extensively, and involve multiple pathways of Myc-dependent cell death and proliferation. First, similar to other oncogenes, supraphysiological levels of Myc induce p19ARF (p14ARF in human) that sequesters MDM2 E3 ubiquitin ligase which leads to stabilization of p53 [33, 35]. Consistent with this finding, loss of p19ARF or p53 allows Myc to immortalize primary cells and synergize with Myc-induced tumorigenesis [52–54]. There is also evidence that ARF can limit Myc activity by direct protein-protein interaction independent of p53 and MDM2 [55, 56]. Therefore, ARF mediates its tumor suppressor activity through p53-dependent and -independent mechanisms to prevent Myc-stimulated proliferation and cellular transformation. In this study, we found that JMJD6 decreases Myc-induced apoptosis, at least in part, by suppressing p19ARF expression. It would be interesting to analyze the cooperation of JMJD6 with other oncogenes whose action is also restricted by ARF.
There are multiple functions assigned to the JMJD6 enzyme. It was first discovered as a histone arginine demethylase  and a lysine hydroxylase of histone and non-histone proteins [29, 30]. The function of JMJD6 as an epigenetic modifier of histone proteins has been widely debated in the literature. The unique property of JMJD6 is its ability to function as an arginine demethylase. Currently, it is the only known arginine demethylase that is capable of targeting both histone and non-histone proteins. In vitro studies have demonstrated the removal of monomethyl groups from histone arginine residues as well as symmetric and asymmetric methyl groups [32, 57]. That implies that JMJD6 may function as a transcriptional activator and transcriptional repressor, as H4R3me2s (symmetric) is a repressive mark [37, 38] and H4R3me2a (asymmetric) methylation is a mark of active chromatin [37, 39, 40]. Consistent with this notion, diminished expression of p19ARF in cells with high JMJD6 levels presented in this study was accompanied by a decrease in histone H4R3me2a modification (a marker of active chromatin) in the p19 promoter. This was associated with reduced p53 expression and reduced apoptosis.
Consequences of lysine hydroxylation of non-histone proteins may be variable. For example, when JMJD6 hydroxylates splice factor U2AF65 [30, 58], it changes VegfR1 mRNA splicing in endothelial cells in favor of the angiogenesis promoting long form of the protein [31, 59]. Another example of non-histone protein modification by JMJD6 has recently identified its ability to hydroxylate Lys382 in p53 protein in the colorectal cancer cell line HCT116 . This modification does not affect p53 protein stability but rather inhibits its transcriptional activity. This is, however, not the case for normal mouse epithelial cells used in this study, as we clearly see decreased p53 at the protein level (Fig. 6). We are currently exploring the possibility that JMJD6 also induces p53 hydroxylation in non-tumorigenic mouse and human cells.
Several publications identified Brd4 as a binding partner of JMJD6 in different cell types [32, 60, 61]. Brd4 is classified as a histone code “reader” that functions by recognition of acetylated histones and recruiting epigenetic modifiers to the target gene . Assembly of the Brd4-JMJD6 complex in distant enhancers of target genes stimulates histone H4R3me2(s) demethylation as well as cap removal from 7SKsnRNA, promoting transcriptional activation regulated by anti-pause enhancers . Interestingly, Myc protein was also implicated in the transcriptional pause release of active genes . When Myc levels were increased, similar to Myc-amplified tumors, Myc also occupied distant enhancers of actively transcribed genes  which raises the question of whether JMJD6 amplification in MMTV-Myc tumors is necessary for the overall increase of transcriptional activation by Myc. Further experimentation will be required to fully understand the extent of Myc-JMJD6 cooperation.
Elevated levels of JMJD6 protein and mRNA were reported for several tumor types, breast, lung, and colon cancers [25–27]. For colon cancer, elevated JMJD6 expression positively correlated with depth of invasion, lymph node metastasis, and advanced tumor node metastasis stage . In agreement with this, we observed a causal role of JMJD6 in dramatically promoting metastasis in Myc-driven tumors (Fig. 9).
Myc is a proto-oncogene that triggers tumorigenesis by increasing cellular proliferation. However, paradoxically, it has been reported that in breast cancer cells, it can act as a metastasis suppressor, decreasing cell invasion in vitro and reducing metastatic burden in mice in vivo . In our experiments, overexpression of JMJD6 in MMTV-Myc tumor cells was able to dramatically increase the number of metastatic nodules in lungs of mice after tail vein injection, as well as invasion and migration of these cells in vitro. In addition, analysis of a large cohort of breast cancer patients (more than 2000) revealed that high JMJD6 expression in human breast cancer with high Myc results in a lower overall survival presumably due to a higher metastatic burden (Fig. 10). This is consistent with a previous report that demonstrated that high JMJD6 was associated with a worse prognosis in ER+ breast cancer . ERα function is regulated in part by methylation of R260, which is required for estrogen-induced complex formation of ERα with Src and PI3K and activation of the Akt pathway . JMJD6 demethylates this R260me2a after treatment with estrogen , which implies a compromised, non-genomic function of ERα. However, it remains unclear how JMJD6 expression leads to a worse prognosis in ER+ patients.
Our studies have revealed important cooperativity between JMJD6 and the Myc proto-oncogene and possibly other types of oncogene-driven breast cancers. These results suggest that JMJD6 protein expression may be used as a prognostic biomarker in future studies. Additionally, since Myc is overexpressed in many breast cancers but not successfully targeted by drugs, inhibiting JMJD6 function may provide a means of inhibiting Myc-dependent tumorigenesis.
DNA copy number
Genomic DNA was extracted from tumors or spleen (as reference) from 8 GEM models of mammary cancer using DNeasy Blood and Tissue kit (Qiagen, USA). Five hundred nanograms of reference or test sample DNA was labeled with Cy3 and Cy5 using Enzo CGH labeling kit (Enzo Life Sciences, USA) according to the manufacturer’s instructions. The labeled DNAs were purified with Amicon Ultra-0.5, 30Kd filters and hybridized to Agilent 44K CGH arrays for 44 h at 65 °C in a rotating microarray hybridization chamber and then washed with buffer 1 for 5 min at room temperature and buffer 2 for 1 min at 37 °C. After a brief rinse in acetonitrile and Agilent stabilizer (Agilent 5185–5979), the arrays were scanned in an array holder covered with the ozone barrier using the Agilent DNA microarray scanner G2539A. The array features were extracted using the Agilent Feature Extraction v10.7. The copy numbers were analyzed using Agilent Genomic Workbench 5.0 program and Nexus Copy Number software (BioDiscovery, USA).
Data have been deposited to the Gene Expression Omnibus (accession number GSE75331). The reviewer access url is http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=yxclqmqynzadxgx&acc=GSE75331.
Gene expression microarray
Total RNA (1 μg) was reverse transcribed with a T7-oligo(dT) primer and labeled with biotin using the Affymetrix One Cycle Target Labeling kit (Affymetrix, USA), following the manufacturer’s protocol. RNA was then labeled and hybridized to the mouse genome 430A 2.0 GeneChip (Affymetrix, USA) and scanned on an Affymetrix GeneChip scanner 3000. Data were collected using Affymetrix GCOS software. Data analysis was described in Zhu et al., 2011 . Data have been deposited to the Gene Expression Omnibus (accession number GSE23938).
NMuMG epithelial cells were purchased from ATCC. MMTV-Myc mammary gland tumor cell lines Myc83 and 88CT1 were obtained from Dr. R. Dickson (Georgetown University). Cells were cultured in DMEM (high glucose) tissue culture media, supplemented with 10 % fetal bovine serum (Life Technologies, USA) and penicillin/streptomycin. Cultures were incubated in humidified air containing 5 % CO2 at 37 °C. Cells in nutrient deprivation experiments were cultured in DMEM glucose-free, glutamine-free media (Life Technologies, USA), supplemented with 10 % dialyzed FBS (dialyzed against 3 changes of 100 volumes of PBS) . The remaining components (glucose or glutamine) were added separately. 4-hydroxytamoxifen was purchased from Sigma and used at a final concentration of 150 nM. Etoposide (Sigma, USA) was added for 24 h at a final concentration of 100 μM.
The lentiviral vectors pLenti6.2 expressing LacZ-V5 or JMJD6-V5 were described by Lee et al.  and were a kind gift of Dr. K. Desai (National Institute of Biomedical Genomics, India). The pre-designed lentiviral shRNA vectors for knock-down of JMJD6, TK1, Tnrc6c, Sumo2, FBF1, Birc5, and Ube2o were cloned into MISSION vector with puromycin resistance and purchased from Sigma. PBabePuro-MycERTM used for Myc induction in cells overexpressing JMJD6 or LacZ was a gift from Dr. G. Evan (UCSF), and pBabeHygro-MycERTM used with shRNA knock-down experiments was a gift from Dr. Y. Lazebnik (Cold Spring Harbor laboratory). Retroviral vector pBabe-p19ARF was a gift from Dr. P. Johnson (NCI).
Real-time quantitative PCR
Total RNA was isolated from cells by TRIzol extraction (Life Technology). One microgram of total RNA was reverse transcribed in a 50-μl reaction using TaqMan Reverse Transcription Reagents (Life Technologies, USA) according to manufacturer’s instructions. Five microliters of the resultant reaction was PCR amplified in a total volume of 20 μl for 40 cycles using an Applied Biosystems 7500 instrument. TaqMan probes were purchased from Life Technologies. Mouse Tbp (TATA binding protein) was used as a standard for normalization. The gene expression for p19ARF was also validated with the specific pre-designed probe Mm.PT.58.8388138 (Integrated DNA Technologies, USA) with SYBR Green qPCR master mix (Bio-Rad, USA). All reactions were performed in triplicate.
For immunoprecipitation with anti-H4R3me2a antibody (Active Motif, USA) cells were fixed with 1 % formaldehyde (Sigma, USA) for 15 min at room temperature. For immunoprecipitation with anti-JMJD6 antibody (Millipore, USA), cells were first fixed with 2 mM disuccinimidyl glutarate (DSG) (ProteoChem, USA) for 45 min at room temperature, followed by two washes with cold PBS and additional fixation with 1 % formaldehyde for 15 min. ChIP assays were performed using a ChIP-IT High Sensitivity Kit from Active Motif (USA). Fixation was stopped by addition of 1/20 volume of stop buffer, and aliquots of 4 × 106 cells were lysed in 400 μl of SDS-lysis buffer. DNA shearing was performed for a total of 16 min at high power setting in BIORUPTOR water bath sonicator (Diagenode, USA). Immunoprecipitation was performed with 1.4 μg of anti-JMJD6 or 10 μg of anti-H4R3me2a antibody or 40 μl of V5-agarose (Sigma) at 4 °C overnight. DNA protein complexes were collected with protein G agarose beads and washed according to manufacturer’s instructions. After elution and reversal of crosslinking, samples were treated with RNase A and proteinase K and purified with QIAQuick PCR columns. Five microliters of the eluted DNA was used for quantitative PCR analysis with SYBR Green qPCR master mix (Bio-Rad, USA). For mouse p19ARF, primers used  were forward primer 5′-TTCCAGGCCTTGCCATCTTCCTAT-3′ and reverse primer 5′-TGGTCTGGCTGCAGTAAAGTAGCA-3′. For human p14ARF, we used two pairs of primers: distal forward 5′-CATTTCTGAGGAAGGGCTACTT-3′ and distal reverse 5′-GCCACCTTTCTGTCTAGTATGG-3′, proximal forward 5′-TCGCCAAGACAACCATTCTAC-3′ and proximal reverse 5′-CGCTTCTTCCTCTTTCCTCTTC-3′. For amplification of mouse p16 promoter, the primers were forward 5′-CAGATTGCCCTCCGATGACTTC-3′ and reverse 5′-TGGACCCGCACAGCAAAGAAGT-3′. Primers for β-globin used as a positive control were provided by Millipore together with ChIP quality anti-JMJD6 antibody.
Cell cytotoxicity assay
Cells were seeded at 10,000 cells per well in 96-well non-transparent plates. The next day, cells were treated with 150 nM 4-hydroxytamoxifen (4-OHT) to induce Myc or with equal volume of EtOH as a control. Twenty-four hour later, cells were washed twice with PBS and the media was changed for “stress media” that included either serum-free media, glucose-free media, glutamine-free media, or 100 μM etoposide. After 18–20 hours the percentage of dead cells relative to the total amount of cells was measured using the luminescent CytoTox-Glo cytotoxicity assay (Promega, USA). Briefly, 50 μl of CytoTox-Glo cytotoxicity reagent was added to each well, mixed and incubated for 15 min at room temperature. After measuring the luminescent signal, 50 μl of lysis reagent containing digitonin was added for another 15 min to lyse the remaining viable cells. After incubation, the second luminescence measurement was taken and the percentage of dead cells was calculated by dividing the first luminescence by second after subtraction of the background signal for media only. Each experiment was performed in triplicates and repeated 2–3 times.
The CytoTox-Glo assay measures the number of dead cells irrespective of how cells died. To confirm that majority of cell die from apoptosis after Myc induction, we repeated the experiments, performing FACS analysis after staining the cells with FITC-Annexin V and propidium iodide using the Apoptosis Detection Kit (BD Pharmingen, USA). The data are shown in Additional file 1: Figure S9.
Migration and invasion assay
2 × 104 cells expressing JMJD6 or LacZ control were placed into the top compartment of transwell Boyden chambers (8 μm, Corning, USA) in serum-free media. The lower compartment contained complete media with 10 % FBS. For the invasion assay, we used the same setting with Matrigel™-covered chambers prepared according to the manufacturer’s instructions (Corning, USA). After 24 h, migrated cells were fixed and stained with the Diff-Quik Stain Set (Dade Behring Inc., USA). The cells that migrated through the pores of the membrane were photographed and quantified using ImageJ software.
Western blot analysis
Protein lysates were prepared in RIPA buffer and 20–50 μg of protein were fractionated in a 10 % SDS-PAGE gel and transferred to a PVDF membrane (Immobilon-P, Millipore). The membrane was blocked with 5 % milk (Bio-Rad Laboratories), probed with antibodies, followed by washing and incubation with HRP-conjugated secondary antibody (Pierce) and developed using the WesternBright ECLTM reagent (Advansta, USA). Antibodies used were anti-V5 (Life Technologies, USA), anti-p53 (Becton Dickinson, USA), anti-Myc, anti-p53-Ser15, anti-p21 (Santa Cruz Biotechnology, USA), anti-Caspase3, anti-PARP, anti-Bmi1, anti-Bcl2, anti-Bcl-xL, anti-Snail, anti-Slug (Cell Signaling, USA), anti-Bcl-W (GeneTex, USA), anti-p19ARF, anti-JMJD6 (Abcam, USA), anti-vimentin, and anti-β-actin (Millipore, USA).
Cells were trypsinized and washed twice with PBS, and 106 cells in 100-μl PBS were surgically implanted into the #4 mammary fat pad of FVB/N mice (Charles River, USA). One week later, staples were removed and tumor volume was measured weekly with a caliper over 1 month. For the lung colonization assay, cells were prepared the same way and 5 × 105 cells in 100-μl PBS were injected by tail vein into FVB/N mice. Three weeks later, mice were euthanized, lungs were formalin fixed and paraffin embedded, and sections were stained with hematoxylin-eosin solution.
All GEM models were described by Zhu et al. .
All mice were treated in accordance with the Guide for the Care and Use of Laboratory Animals (NIH publication no. 86–23, 1985) under an animal protocol (LCBG-063) approved by the IACUC of the National Cancer Institute (NCI).
Lungs from mice were formalin fixed and embedded in paraffin, and three 5-μm step sections were obtained for each lung. Sections were deparaffinized by standard treatment with xylene and ethanol. The following procedure was performed using the ApopTag Peroxidase In Situ Apoptosis Detection Kit (Millipore, USA) according to the manufacturer’s instructions. Briefly, the tissue was pretreated with proteinase K (20 μg/ml), endogenous peroxidase was quenched by 3 % hydrogen peroxide for 5 min, and slides were incubated with TdT enzyme for 1 h at 37 °C. After washing, slides were incubated with anti-digoxigenin conjugate for 30 min at room temperature, washed and developed with peroxidase substrate. Slides were counterstained with methyl green and after dehydration mounted under glass coverslips using Permount media (Fisher Scientific, USA). Three photographs were taken from each slide and analyzed using ImageJ software.
Cells expressing JMJD6 or LacZ were grown on cover slips, fixed with 4 % paraformaldehyde for 15 min, and permeabilized with 0.2 % Triton X-100 for 4 min, and immunostaining was performed by 1-h incubation with anti-V5 antibody (Life Technologies, USA) followed by 30-min incubation with the secondary goat anti-mouse antibody, conjugated with Alexa Fluor 568 (Life Technologies, USA). Nuclei were counterstained with DAPI and images were taken with Nikon Eclipse E800 fluorescent microscope (Nikon).
Statistical analysis was performed with GraphPad Prism5 software that shows mean values with 95 % confidence interval. The p values were calculated using Student’s t test (two tailed). The TaqMan data were analyzed using ΔΔCt method and presented as mean fold change ±SEM. The Pearson correlation coefficient for mouse mammary gland tumors was determined to assess the correlation of JMJD6 and p19ARF expressions.
The expression data for JMJD6 and Myc in human mammary gland tumors were obtained from the METABRIC . For each gene, the binary expression values (high or low) were defined by dichotomizing continuous expression values using the median as the threshold. We used the two genes JMJD6 and Myc jointly to define four patient sub-groups: JMJD6 high/Myc low, JMJD6 low/Myc low, JMJD6 high/Myc high, and JMJD6 low/Myc high. To examine the effect of JMJD6, we compared the Kaplan-Meier curves between the two sub-groups JMJD6 high/Myc low versus JMJD6 low/Myc low. We also compared JMJD6 high/Myc high versus JMJD6 low/Myc high. p values were obtained using log-rank test.
This work was supported by the Intramural Research Program of the NIH, CCR, and NCI. We would like to thank Dr. Evan and Dr. Lazebnik for pBabe-MycERTM puro and hygro plasmids, Dr. Desai for pLV-JMJD6 and LacZ plasmids, Dr. Johnson for pBabe-p19ARF plasmid, Dr. Liu for depositing data into Gene Expression Omnibus, and members of Green lab for helpful discussions. We are grateful to Dr. Kent Hunter for critical reading of the manuscript. We would also like to thank the Array Core Facility of the National Cancer Institute for performing array CGH analysis.
- Green JE, Shibata MA, Yoshidome K, Liu ML, Jorcyk C, Anver MR, et al. The C3(1)/SV40 T-antigen transgenic mouse model of mammary cancer: ductal epithelial cell targeting with multistage progression to carcinoma. Oncogene. 2000;19:1020–7.View ArticlePubMedGoogle Scholar
- Hennighausen L. Mouse models for breast cancer. Oncogene. 2000;19:966–7.View ArticlePubMedGoogle Scholar
- Herschkowitz JI, Simin K, Weigman VJ, Mikaelian I, Usary J, Hu Z, et al. Identification of conserved gene expression features between murine mammary carcinoma models and human breast tumors. Genome Biol. 2007;8:R76.View ArticlePubMedPubMed CentralGoogle Scholar
- Foulkes WD, Stefansson IM, Chappuis PO, Begin LR, Goffin JR, Wong N, et al. Germline BRCA1 mutations and a basal epithelial phenotype in breast cancer. J Natl Cancer Inst. 2003;95:1482–5.View ArticlePubMedGoogle Scholar
- Sorlie T, Tibshirani R, Parker J, Hastie T, Marron JS, Nobel A, et al. Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci U S A. 2003;100:8418–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Inaki K, Menghi F, Woo XY, Wagner JP, Jacques PE, Lee YF, et al. Systems consequences of amplicon formation in human breast cancer. Genome Res. 2014;24:1559–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Andrechek ER, Cardiff RD, Chang JT, Gatza ML, Acharya CR, Potti A, et al. Genetic heterogeneity of Myc-induced mammary tumors reflecting diverse phenotypes including metastatic potential. Proc Natl Acad Sci U S A. 2009;106:16387–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Silva GO, He X, Parker JS, Gatza ML, Carey LA, Hou JP, et al. Cross-species DNA copy number analyses identifies multiple 1q21-q23 subtype-specific driver genes for breast cancer. Breast Cancer Res Treat. 2015;152:347–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Alitalo K, Schwab M. Oncogene amplification in tumor cells. Adv Cancer Res. 1986;47:235–81.View ArticlePubMedGoogle Scholar
- Liao DJ, Dickson RB. c-Myc in breast cancer. Endocr Relat Cancer. 2000;7:143–64.View ArticlePubMedGoogle Scholar
- Nesbit CE, Tersak JM, Prochownik EV. Myc oncogenes and human neoplastic disease. Oncogene. 1999;18:3004–16.View ArticlePubMedGoogle Scholar
- Shiu RP, Watson PH, Dubik D. c-Myc oncogene expression in estrogen-dependent and -independent breast cancer. Clin Chem. 1993;39:353–5.PubMedGoogle Scholar
- Hodgson JG, Malek T, Bornstein S, Hariono S, Ginzinger DG, Muller WJ, et al. Copy number aberrations in mouse breast tumors reveal loci and genes important in tumorigenic receptor tyrosine kinase signaling. Cancer Res. 2005;65:9695–704.View ArticlePubMedGoogle Scholar
- McCormack SJ, Weaver Z, Deming S, Natarajan G, Torri J, Johnson MD, et al. Myc/p53 interactions in transgenic mouse mammary development, tumorigenesis and chromosomal instability. Oncogene. 1998;16:2755–66.View ArticlePubMedGoogle Scholar
- Montagna C, Lyu MS, Hunter K, Lukes L, Lowther W, Reppert T, et al. The septin 9 (MSF) gene is amplified and overexpressed in mouse mammary gland adenocarcinomas and human breast cancer cell lines. Cancer Res. 2003;63:2179–87.PubMedGoogle Scholar
- Weaver Z, Montagna C, Xu X, Howard T, Gadina M, Brodie SG, et al. Mammary tumors in mice conditionally mutant for Brca1 exhibit gross genomic instability and centrosome amplification yet display a recurring distribution of genomic imbalances that is similar to human breast cancer. Oncogene. 2002;21:5097–107.View ArticlePubMedGoogle Scholar
- Liu ML, Shibata MA, Von Lintig FC, Wang W, Cassenaer S, Boss GR, et al. Haploid loss of Ki-ras delays mammary tumor progression in C3 (1)/SV40 Tag transgenic mice. Oncogene. 2001;20:2044–9.View ArticlePubMedGoogle Scholar
- Liu ML, Von Lintig FC, Liyanage M, Shibata MA, Jorcyk CL, Ried T, et al. Amplification of Ki-ras and elevation of MAP kinase activity during mammary tumor progression in C3(1)/SV40 Tag transgenic mice. Oncogene. 1998;17:2403–11.View ArticlePubMedGoogle Scholar
- Kallioniemi A, Kallioniemi OP, Piper J, Tanner M, Stokke T, Chen L, et al. Detection and mapping of amplified DNA sequences in breast cancer by comparative genomic hybridization. Proc Natl Acad Sci U S A. 1994;91:2156–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Somiari SB, Shriver CD, He J, Parikh K, Jordan R, Hooke J, et al. Global search for chromosomal abnormalities in infiltrating ductal carcinoma of the breast using array-comparative genomic hybridization. Cancer Genet Cytogenet. 2004;155:108–18.View ArticlePubMedGoogle Scholar
- Toffoli S, Bar I, Abdel-Sater F, Delree P, Hilbert P, Cavallin F, et al. Identification by array comparative genomic hybridization of a new amplicon on chromosome 17q highly recurrent in BRCA1 mutated triple negative breast cancer. Breast Cancer Res. 2014;16:466.View ArticlePubMedPubMed CentralGoogle Scholar
- Waldman FM, DeVries S, Chew KL, Moore 2nd DH, Kerlikowske K, Ljung BM. Chromosomal alterations in ductal carcinomas in situ and their in situ recurrences. J Natl Cancer Inst. 2000;92:313–20.View ArticlePubMedGoogle Scholar
- Weaver ZA, McCormack SJ, Liyanage M, du Manoir S, Coleman A, Schrock E, et al. A recurring pattern of chromosomal aberrations in mammary gland tumors of MMTV-cMyc transgenic mice. Genes Chromosomes Cancer. 1999;25:251–60.View ArticlePubMedGoogle Scholar
- Littlewood TD, Hancock DC, Danielian PS, Parker MG, Evan GI. A modified oestrogen receptor ligand-binding domain as an improved switch for the regulation of heterologous proteins. Nucleic Acids Res. 1995;23:1686–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee YF, Miller LD, Chan XB, Black MA, Pang B, Ong CW, et al. JMJD6 is a driver of cellular proliferation and motility and a marker of poor prognosis in breast cancer. Breast Cancer Res. 2012;14:R85.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang F, He L, Huangyang P, Liang J, Si W, Yan R, et al. JMJD6 promotes colon carcinogenesis through negative regulation of p53 by hydroxylation. PLoS Biol. 2014;12, e1001819.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang J, Ni SS, Zhao WL, Dong XC, Wang JL. High expression of JMJD6 predicts unfavorable survival in lung adenocarcinoma. Tumour Biol. 2013;34:2397–401.View ArticlePubMedGoogle Scholar
- Poulard C, Rambaud J, Lavergne E, Jacquemetton J, Renoir JM, Tredan O, et al. Role of JMJD6 in breast tumourigenesis. PLoS One. 2015;10, e0126181.View ArticlePubMedPubMed CentralGoogle Scholar
- Unoki M, Masuda A, Dohmae N, Arita K, Yoshimatsu M, Iwai Y, et al. Lysyl 5-hydroxylation, a novel histone modification, by Jumonji domain containing 6 (JMJD6). J Biol Chem. 2013;288:6053–62.View ArticlePubMedPubMed CentralGoogle Scholar
- Webby CJ, Wolf A, Gromak N, Dreger M, Kramer H, Kessler B, et al. Jmjd6 catalyses lysyl-hydroxylation of U2AF65, a protein associated with RNA splicing. Science. 2009;325:90–3.View ArticlePubMedGoogle Scholar
- Boeckel JN, Guarani V, Koyanagi M, Roexe T, Lengeling A, Schermuly RT, et al. Jumonji domain-containing protein 6 (Jmjd6) is required for angiogenic sprouting and regulates splicing of VEGF-receptor 1. Proc Natl Acad Sci U S A. 2011;108:3276–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu W, Ma Q, Wong K, Li W, Ohgi K, Zhang J, et al. Brd4 and JMJD6-associated anti-pause enhancers in regulation of transcriptional pause release. Cell. 2013;155:1581–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen D, Kon N, Zhong J, Zhang P, Yu L, Gu W. Differential effects on ARF stability by normal versus oncogenic levels of c-Myc expression. Mol Cell. 2013;51:46–56.View ArticlePubMedGoogle Scholar
- Zhang Y, Xiong Y, Yarbrough WG. ARF promotes MDM2 degradation and stabilizes p53: ARF-INK4a locus deletion impairs both the Rb and p53 tumor suppression pathways. Cell. 1998;92:725–34.View ArticlePubMedGoogle Scholar
- Zindy F, Eischen CM, Randle DH, Kamijo T, Cleveland JL, Sherr CJ, et al. Myc signaling via the ARF tumor suppressor regulates p53-dependent apoptosis and immortalization. Genes Dev. 1998;12:2424–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Chang B, Chen Y, Zhao Y, Bruick RK. JMJD6 is a histone arginine demethylase. Science. 2007;318:444–7.View ArticlePubMedGoogle Scholar
- Wysocka J, Allis CD, Coonrod S. Histone arginine methylation and its dynamic regulation. Front Biosci. 2006;11:344–55.View ArticlePubMedGoogle Scholar
- Xu X, Hoang S, Mayo MW, Bekiranov S. Application of machine learning methods to histone methylation ChIP-Seq data reveals H4R3me2 globally represses gene expression. BMC Bioinf. 2010;11:396.Google Scholar
- Li X, Hu X, Patel B, Zhou Z, Liang S, Ybarra R, et al. H4R3 methylation facilitates beta-globin transcription by regulating histone acetyltransferase binding and H3 acetylation. Blood. 2010;115:2028–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang Y, Lu Y, Espejo A, Wu J, Xu W, Liang S, et al. TDRD3 is an effector molecule for arginine-methylated histone marks. Mol Cell. 2010;40:1016–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Beverly LJ, Varmus HE. Myc-induced myeloid leukemogenesis is accelerated by all six members of the antiapoptotic BCL family. Oncogene. 2009;28:1274–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Eischen CM, Packham G, Nip J, Fee BE, Hiebert SW, Zambetti GP, et al. Bcl-2 is an apoptotic target suppressed by both c-Myc and E2F-1. Oncogene. 2001;20:6983–93.View ArticlePubMedGoogle Scholar
- Shiio Y, Suh KS, Lee H, Yuspa SH, Eisenman RN, Aebersold R. Quantitative proteomic analysis of Myc-induced apoptosis: a direct role for Myc induction of the mitochondrial chloride ion channel, mtCLIC/CLIC4. J Biol Chem. 2006;281:2750–6.View ArticlePubMedGoogle Scholar
- D’Cruz CM, Gunther EJ, Boxer RB, Hartman JL, Sintasath L, Moody SE, et al. c-Myc induces mammary tumorigenesis by means of a preferred pathway involving spontaneous Kras2 mutations. Nat Med. 2001;7:235–9.View ArticlePubMedGoogle Scholar
- Stewart TA, Pattengale PK, Leder P. Spontaneous mammary adenocarcinomas in transgenic mice that carry and express MTV/Myc fusion genes. Cell. 1984;38:627–37.View ArticlePubMedGoogle Scholar
- Liu H, Radisky DC, Yang D, Xu R, Radisky ES, Bissell MJ, et al. Myc suppresses cancer metastasis by direct transcriptional silencing of alphav and beta3 integrin subunits. Nat Cell Biol. 2012;14:567–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Evan GI, Wyllie AH, Gilbert CS, Littlewood TD, Land H, Brooks M, et al. Induction of apoptosis in fibroblasts by c-Myc protein. Cell. 1992;69:119–28.View ArticlePubMedGoogle Scholar
- Felsher DW, Bishop JM. Transient excess of Myc activity can elicit genomic instability and tumorigenesis. Proc Natl Acad Sci U S A. 1999;96:3940–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Gibson AW, Cheng T, Johnston RN. Apoptosis induced by c-Myc overexpression is dependent on growth conditions. Exp Cell Res. 1995;218:351–8.View ArticlePubMedGoogle Scholar
- Hoffman B, Liebermann DA. Apoptotic signaling by c-Myc. Oncogene. 2008;27:6462–72.View ArticlePubMedGoogle Scholar
- Nilsson JA, Cleveland JL. Myc pathways provoking cell suicide and cancer. Oncogene. 2003;22:9007–21.View ArticlePubMedGoogle Scholar
- Eischen CM, Weber JD, Roussel MF, Sherr CJ, Cleveland JL. Disruption of the ARF-MDM2-p53 tumor suppressor pathway in Myc-induced lymphomagenesis. Genes Dev. 1999;13:2658–69.View ArticlePubMedPubMed CentralGoogle Scholar
- Hermeking H, Eick D. Mediation of c-Myc-induced apoptosis by p53. Science. 1994;265:2091–3.View ArticlePubMedGoogle Scholar
- Schmitt CA, McCurrach ME, de Stanchina E, Wallace-Brodeur RR, Lowe SW. INK4a/ARF mutations accelerate lymphomagenesis and promote chemoresistance by disabling p53. Genes Dev. 1999;13:2670–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Datta A, Nag A, Pan W, Hay N, Gartel AL, Colamonici O, et al. Myc-ARF (alternate reading frame) interaction inhibits the functions of Myc. J Biol Chem. 2004;279:36698–707.View ArticlePubMedGoogle Scholar
- Qi Y, Gregory MA, Li Z, Brousal JP, West K, Hann SR. p19ARF directly and differentially controls the functions of c-Myc independently of p53. Nature. 2004;431:712–7.View ArticlePubMedGoogle Scholar
- Poulard C, Rambaud J, Hussein N, Corbo L, Le Romancer M. JMJD6 regulates ERalpha methylation on arginine. PLoS One. 2014;9, e87982.View ArticlePubMedPubMed CentralGoogle Scholar
- Heim A, Grimm C, Muller U, Haussler S, Mackeen MM, Merl J, et al. Jumonji domain containing protein 6 (JMJD6) modulates splicing and specifically interacts with arginine-serine-rich (RS) domains of SR- and SR-like proteins. Nucleic Acids Res. 2014;42:7833–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Qi JW, Qin TT, Xu LX, Zhang K, Yang GL, Li J, et al. TNFSF15 inhibits vasculogenesis by regulating relative levels of membrane-bound and soluble isoforms of VEGF receptor 1. Proc Natl Acad Sci U S A. 2013;110:13863–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Alsarraj J, Faraji F, Geiger TR, Mattaini KR, Williams M, Wu J, et al. BRD4 short isoform interacts with RRP1B, SIPA1 and components of the LINC complex at the inner face of the nuclear membrane. PLoS One. 2013;8, e80746.View ArticlePubMedPubMed CentralGoogle Scholar
- Rahman S, Sowa ME, Ottinger M, Smith JA, Shi Y, Harper JW, et al. The Brd4 extraterminal domain confers transcription activation independent of pTEFb by recruiting multiple proteins, including NSD3. Mol Cell Biol. 2011;31:2641–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Shi J, Vakoc CR. The mechanisms behind the therapeutic activity of BET bromodomain inhibition. Mol Cell. 2014;54:728–36.View ArticlePubMedGoogle Scholar
- Rahl PB, Lin CY, Seila AC, Flynn RA, McCuine S, Burge CB, et al. c-Myc regulates transcriptional pause release. Cell. 2010;141:432–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Lin CY, Loven J, Rahl PB, Paranal RM, Burge CB, Bradner JE, et al. Transcriptional amplification in tumor cells with elevated c-Myc. Cell. 2012;151:56–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Le Romancer M, Treilleux I, Leconte N, Robin-Lespinasse Y, Sentis S, Bouchekioua-Bouzaghou K, et al. Regulation of estrogen rapid signaling through arginine methylation by PRMT1. Mol Cell. 2008;31:212–21.View ArticlePubMedGoogle Scholar
- Zhu M, Yi M, Kim CH, Deng C, Li Y, Medina D, et al. Integrated miRNA and mRNA expression profiling of mouse mammary tumor models identifies miRNA signatures associated with mammary tumor lineage. Genome Biol. 2011;12:R77.View ArticlePubMedPubMed CentralGoogle Scholar
- Yuneva M, Zamboni N, Oefner P, Sachidanandam R, Lazebnik Y. Deficiency in glutamine but not glucose induces Myc-dependent apoptosis in human cells. J Cell Biol. 2007;178:93–105.View ArticlePubMedPubMed CentralGoogle Scholar
- Zheng Y, Zhao YD, Gibbons M, Abramova T, Chu PY, Ash JD, et al. Tgfbeta signaling directly induces ARF promoter remodeling by a mechanism involving Smads 2/3 and p38 MAPK. J Biol Chem. 2010;285:35654–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Curtis C, Shah SP, Chin SF, Turashvili G, Rueda OM, Dunning MJ, et al. The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature. 2012;486:346–52.PubMedPubMed CentralGoogle Scholar