JMJD6 orchestrates a transcriptional program in favor of endocrine resistance in ER+ breast cancer cells.

High expression of Jumonji domain containing protein 6 (JMJD6) is strongly associated with poor prognosis in estrogen receptor positive (ER+) breast cancer. We overexpressed JMJD6 in MCF7 cells (JOE cells) and performed RNA-seq analysis. 76% of differentially expressed genes (DEGs) overlapped with ER target genes. Pathway analysis revealed that JMJD6 upregulated a larger subset of genes related to cell proliferation as compared to ER. Interestingly, JOE cells showed a decrease in ER target gene expression prompting us to check ER levels. Indeed, JOE cells showed a significant decrease in both ESR1 and ER levels and JMJD6 siRNA transfection increased the expression of both. Additionally, JOE cells showed increased RET and ERK1 expression, events associated with resistance to endocrine therapy. Accordingly, JOE cells displayed lower sensitivity and survived better at higher doses of 4-hydroxy tamoxifen (Tam) as compared to parental MCF-7 cells. Conversely, LTED-I and TAM R that resist Tam induced death, showed high expression of JMJD6. Further, JMJD6 siRNA treatment decreased growth and improved Tam sensitivity in TAM R. Comparison of JOE DEGs with known Tam signature genes showed a substantial overlap. Overall, these data suggest that blocking ER alone in patients may not eradicate proliferation of JMJD6 expressing ER+ cells and JMJD6 may predispose and sustain endocrine therapy resistance. We propose that immunostaining for JMJD6 could be developed as a potential marker for predicting endocrine therapy resistance. Further, antagonizing JMJD6 action in women expressing higher amounts of this protein, may offer a greater clinical benefit than endocrine therapy.

Both EZH2 and JMJD6 regulate cell cycle genes in breast cancer.

BACKGROUND: Strong evidences support the critical role of Jumonji domain containing 6 (JMJD6) in progression of breast cancer. Here we explore potential partners that coregulate gene expression, to understand additional pathways that are activated by higher amounts of JMJD6. METHODS: We used Gene Set Enrichment Analysis (GSEA) data to identify factors that display gene expression similar to cells treated with JMJD6 siRNA. Using chromatin immunoprecipitations (ChIP) against genomic regions that bind JMJD6 identified by in house and public database Encyclopaedia of DNA Elements (ENCODE), we confirmed JMJD6 occupancy by ChIP PCR. We tested the association of co-regulated genes with patient prognosis using The Cancer Genome Atlas (TCGA) datasets. RESULTS: JMJD6 profiles overlapped with those of Enhancer of Zeste homolog 2 (EZH2) and together they appear to co-regulate a unique cassette of genes in both ER+ and ER- cells. 496 genes including aurora kinases, which are currently being tested as novel therapeutic targets in breast cancer were co-regulated in MDA MB 231 cells. JMJD6 and EZH2 neither inter-regulated nor physically interacted with one another. Since both proteins are chromatin modulators, we performed ChIP linked PCR analysis and show that JMJD6 bound in the neighbourhood of co-regulated genes, though EZH2 data did not show any peaks within 100 kb of these sites. Alignment of binding site sequences suggested that atleast two types of binding partners could offer their DNA binding properties to enrich JMJD6 at regulatory sites. In clinical samples, JMJD6 and EZH2 expression significantly correlated in both normal and tumor samples, however the strongest correlation was observed in triple-negative breast cancer (TNBC) subtype. Co-expression of JMJD6 and EZH2 imposed poorer prognosis in breast cancer. CONCLUSIONS: JMJD6 and EZH2 regulate the same crucial cell cycle regulatory and therapeutic targets but their mechanisms appear to be independent of each other. Blocking of a single molecule may not axe cell proliferation completely and blocking both JMJD6 and EZH2 simultaneously may be more effective in breast cancer patients.

Pathways to Endocrine Therapy Resistance in Breast Cancer.

Breast cancers with positive expression of Estrogen Receptor (ER+) are treated with anti-hormone/endocrine therapy which targets the activity of the receptor, the half-life of the receptor or the availability of estrogen. This has significantly decreased mortality in women with ER+ breast cancer, however, about 25-30% of treated women run the risk or recurrence due to either intrinsic or acquired resistance to endocrine therapies. While ER itself is a predictor of response to such therapies, there exists a need to find more biomarkers and novel targets to treat resistant tumors. In this review, we summarize the known mechanisms and describe the ability of genomics in unraveling rare mutations and gene rearrangements that may impact the development of resistance and therefore treatment of ER+ breast cancer in the near future.

JMJD6 induces HOTAIR, an oncogenic lincRNA, by physically interacting with its proximal promoter.

Using microarray analysis, we found that HOX transcript antisense intergenic RNA (HOTAIR) is up-regulated by Jumonji domain containing-6 (JMJD6), a bifunctional lysyl hydroxylase and arginine demethylase. In breast cancer, both JMJD6 and HOTAIR RNAs increase tumor growth and associate with poor prognosis but no molecular relationship between them is known. We show that overexpression of JMJD6 increased HOTAIR expression and JMJD6 siRNAs suppressed it in ER+ MCF-7, triple negative MDA-MB-231 and non-breast cancer HEK 293 cells. Therefore, JMJD6 regulates HOTAIR independent of ER status. Using various deletion constructs spanning (-1874 to +50) of the HOTAIR promoter, we identified pHP216 (-216 to +50 bp) as the smallest construct that retained maximal JMJD6 responsiveness. In ChIP assays, JMJD6 bound this region suggesting that JMJD6 may be directly recruited to the HOTAIR promoter. Mutant JMJD6H187A that is devoid of enzymatic activity could bind this site but failed to induce transcription. ChIP and electromobility shift assays identified a JMJD6 interaction region from (-123 to -103 bp) within the HOTAIR promoter. In tumor samples but not normal breast tissue, the expression of JMJD6 linearly correlated with HOTAIR suggesting that JMJD6-mediated up-regulation may occur specifically in tumors. Further, concurrent high expression of both genes correlated with poor survival when individual expression of either gene showed no significant association in TCGA datasets. We propose that high JMJD6 expression may achieve higher levels of HOTAIR in breast tumors. Further, since high levels of HOTAIR promote metastasis and death, blocking JMJD6 may be useful in preventing such events.

Whole-genome reconstruction and mutational signatures in gastric cancer.

BACKGROUND: Gastric cancer is the second highest cause of global cancer mortality. To explore the complete repertoire of somatic alterations in gastric cancer, we combined massively parallel short read and DNA paired-end tag sequencing to present the first whole-genome analysis of two gastric adenocarcinomas, one with chromosomal instability and the other with microsatellite instability. RESULTS: Integrative analysis and de novo assemblies revealed the architecture of a wild-type KRAS amplification, a common driver event in gastric cancer. We discovered three distinct mutational signatures in gastric cancer--against a genome-wide backdrop of oxidative and microsatellite instability-related mutational signatures, we identified the first exome-specific mutational signature. Further characterization of the impact of these signatures by combining sequencing data from 40 complete gastric cancer exomes and targeted screening of an additional 94 independent gastric tumors uncovered ACVR2A, RPL22 and LMAN1 as recurrently mutated genes in microsatellite instability-positive gastric cancer and PAPPA as a recurrently mutated gene in TP53 wild-type gastric cancer. CONCLUSIONS: These results highlight how whole-genome cancer sequencing can uncover information relevant to tissue-specific carcinogenesis that would otherwise be missed from exome-sequencing data.

JMJD6 is a driver of cellular proliferation and motility and a marker of poor prognosis in breast cancer.

INTRODUCTION: We developed an analytic strategy that correlates gene expression and clinical outcomes as a means to identify novel candidate oncogenes operative in breast cancer. This analysis, followed by functional characterization, resulted in the identification of Jumonji Domain Containing 6 (JMJD6) protein as a novel driver of oncogenic properties in breast cancer. METHODS: Through microarray informatics, Cox proportional hazards regression was used to analyze the correlation between gene expression and distant metastasis-free survival (DMFS) of patients in 14 independent breast cancer cohorts. JMJD6 emerged as a top candidate gene robustly associated with poor patient survival. Immunohistochemistry, siRNA-mediated silencing, and forced overexpression of JMJD6 in cell-based assays elucidated molecular mechanisms of JMJD6 action in breast cancer progression and shed light on the clinical breast cancer subtypes relevant to JMJD6 action. RESULTS: JMJD6 was expressed at highest levels in tumors associated with worse outcomes, including ER- and basal-like, Claudin-low, Her2-enriched, and ER+ Luminal B tumors. High nuclear JMJD6 protein was associated with ER negativity, advanced grade, and poor differentiation in tissue microarrays. Separation of ER+/LN- patients that received endocrine monotherapy indicated that JMJD6 is predictive of poor outcome in treatment-specific subgroups. In breast cancer cell lines, loss of JMJD6 consistently resulted in suppressed proliferation but not apoptosis, whereas forced stable overexpression increased growth. In addition, knockdown of JMJD6 in invasive cell lines, such as MDA-MB231, decreased motility and invasion, whereas overexpression in MCF-7 cells slightly promoted motility but did not confer invasive growth. Microarray analysis showed that the most significant transcriptional changes occurred in cell-proliferation genes and genes of the TGF-ß tumor-suppressor pathway. High proliferation was characterized by constitutively high cyclin E protein levels. The inverse relation of JMJD6 expression with TGF-ß2 could be extrapolated to the breast cancer cohorts, suggesting that JMJD6 may affect similar pathways in primary breast cancer. CONCLUSIONS: JMJD6 is a novel biomarker of tumor aggressiveness with functional implications in breast cancer growth and migration.

Combined genomic and phenotype screening reveals secretory factor SPINK1 as an invasion and survival factor associated with patient prognosis in breast cancer.

Secretory factors that drive cancer progression are attractive immunotherapeutic targets. We used a whole-genome data-mining approach on multiple cohorts of breast tumours annotated for clinical outcomes to discover such factors. We identified Serine protease inhibitor Kazal-type 1 (SPINK1) to be associated with poor survival in estrogen receptor-positive (ER+) cases. Immunohistochemistry showed that SPINK1 was absent in normal breast, present in early and advanced tumours, and its expression correlated with poor survival in ER+ tumours. In ER- cases, the prognostic effect did not reach statistical significance. Forced expression and/or exposure to recombinant SPINK1 induced invasiveness without affecting cell proliferation. However, down-regulation of SPINK1 resulted in cell death. Further, SPINK1 overexpressing cells were resistant to drug-induced apoptosis due to reduced caspase-3 levels and high expression of Bcl2 and phospho-Bcl2 proteins. Intriguingly, these anti-apoptotic effects of SPINK1 were abrogated by mutations of its protease inhibition domain. Thus, SPINK1 affects multiple aggressive properties in breast cancer: survival, invasiveness and chemoresistance. Because SPINK1 effects are abrogated by neutralizing antibodies, we suggest that SPINK1 is a viable potential therapeutic target in breast cancer.

An oestrogen-receptor-alpha-bound human chromatin interactome.

Genomes are organized into high-level three-dimensional structures, and DNA elements separated by long genomic distances can in principle interact functionally. Many transcription factors bind to regulatory DNA elements distant from gene promoters. Although distal binding sites have been shown to regulate transcription by long-range chromatin interactions at a few loci, chromatin interactions and their impact on transcription regulation have not been investigated in a genome-wide manner. Here we describe the development of a new strategy, chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) for the de novo detection of global chromatin interactions, with which we have comprehensively mapped the chromatin interaction network bound by oestrogen receptor alpha (ER-alpha) in the human genome. We found that most high-confidence remote ER-alpha-binding sites are anchored at gene promoters through long-range chromatin interactions, suggesting that ER-alpha functions by extensive chromatin looping to bring genes together for coordinated transcriptional regulation. We propose that chromatin interactions constitute a primary mechanism for regulating transcription in mammalian genomes.

Regulation of mRNAs encoding MMP-9 and MMP-2, and their inhibitors TIMP-1 and TIMP-2 by androgens in the rat ventral prostate.

The relative activities of matrix metalloproteinases (MMPs) and their natural inhibitors (tissue inhibitors of matrix metalloproteinases, TIMPs) determine the extent of matrix degradation in any tissue. Their identification and characterization is key towards understanding remodeling of the prostate in the context of both castration induced atrophy and tumor invasion and metastasis. Although the expression of MMPs and TIMPs in prostate tumors has been reported, their regulation by androgens has not been studied. Here, we show that androgen ablation by castration increases the steady state mRNA levels of MMP-9, MMP-2, TIMP-1 and TIMP-2. Blockade of the androgen receptor using flutamide, however, has differential effects on the steady state mRNA expressions of these genes. We also show that both castration and flutamide treatment cause enhanced expression of a high molecular weight gelatinolytic activity in the rat ventral prostate (RVP). Actinomycin D does not affect the increase in steady state mRNA levels of MMP-9 and TIMP-1. Furthermore we show that actinomycin D alone enhances the steady state mRNA and protein levels of these genes. Using RNA gel shift assay with 3'-UTR of TIMP-1, we show that an RNA binding protein is induced following castration. Taken together our data suggest that the induction of MMP-2, MMP-9, TIMP-1 and TIMP-2 mRNAs post-castration could be at least in part due to post-transcriptional stabilization.

Smad3 deficiency alters key structural elements of the extracellular matrix and mechanotransduction of wound closure.

The loss of TGFbeta or its downstream mediator, Smad3, key players in tissue repair, accelerates closure of incisional wounds in mice. In contrast, we now report that excisional ear wounds in mice lacking Smad3 enlarge compared with wild-type controls resulting from changes in extracellular matrix molecules, which alter the mechanotransduction properties of these wounds. Specifically, levels of elastin and glycosoaminoglycans are increased, collagen fibers are more compactly organized, and matrix modulators like integrins, TGFbeta1, and matrix metalloproteinases (MMPs) are altered both basally and after wounding in Smad3 knockout mice. Mechanical testing of dorsal skin correlates these changes in matrix composition with functional parameters, specifically an increased elastic modulus, suggesting an imbalance of tissue forces. We propose that the altered mechanical elastic properties translate into a persistent retractile force that is opposed by decreased wound contractile forces contributing to the enlarging ear wound in Smad3 knockout mice. These studies highlight a previously undescribed role for Smad3 in the mechanotransduction of matrix unsupported ear wound closure.

Changes in gene expression during the development of mammary tumors in MMTV-Wnt-1 transgenic mice.

BACKGROUND: In human breast cancer normal mammary cells typically develop into hyperplasia, ductal carcinoma in situ, invasive cancer, and metastasis. The changes in gene expression associated with this stepwise progression are unclear. Mice transgenic for mouse mammary tumor virus (MMTV)-Wnt-1 exhibit discrete steps of mammary tumorigenesis, including hyperplasia, invasive ductal carcinoma, and distant metastasis. These mice might therefore be useful models for discovering changes in gene expression during cancer development. RESULTS: We used cDNA microarrays to determine the expression profiles of five normal mammary glands, seven hyperplastic mammary glands and 23 mammary tumors from MMTV-Wnt-1 transgenic mice, and 12 mammary tumors from MMTV-Neu transgenic mice. Adipose tissues were used to control for fat cells in the vicinity of the mammary glands. In these analyses, we found that the progression of normal virgin mammary glands to hyperplastic tissues and to mammary tumors is accompanied by differences in the expression of several hundred genes at each step. Some of these differences appear to be unique to the effects of Wnt signaling; others seem to be common to tumors induced by both Neu and Wnt-1 oncogenes. CONCLUSION: We described gene-expression patterns associated with breast-cancer development in mice, and identified genes that may be significant targets for oncogenic events. The expression data developed provide a resource for illuminating the molecular mechanisms involved in breast cancer development, especially through the identification of genes that are critical in cancer initiation and progression.

Gene expression profiling identifies a unique androgen-mediated inflammatory/immune signature and a PTEN (phosphatase and tensin homolog deleted on chromosome 10)-mediated apoptotic response specific to the rat ventral prostate.

Understanding androgen regulation of gene expression is critical for deciphering mechanisms responsible for the transition from androgen-responsive (AR) to androgen-independent (AI) prostate cancer (PCa). To identify genes differentially regulated by androgens in each prostate lobe, the rat castration model was used. Microarray analysis was performed to compare dorsolateral (DLP) and ventral prostate (VP) samples from sham-castrated, castrated, and testosterone-replenished castrated rats. Our data demonstrate that, after castration, the VP and the DLP differed in the number of genes with altered expression (1496 in VP vs. 256 in DLP) and the nature of pathways modulated. Gene signatures related to apoptosis and immune response specific to the ventral prostate were identified. Microarray and RT-PCR analyses demonstrated the androgen repression of IGF binding protein-3 and -5, CCAAT-enhancer binding protein-delta, and phosphatase and tensin homolog deleted on chromosome 10 (PTEN) genes, previously implicated in apoptosis. We show that PTEN protein was increased only in the luminal epithelial cells of the VP, suggesting that it may be a key mediator of VP apoptosis in the absence of androgens. The castration-induced immune/inflammatory gene cluster observed specifically in the VP included IL-15 and IL-18. Immunostaining of the VP, but not the DLP, showed an influx of T cells, macrophages, and mast cells, suggesting that these cells may be the source of the immune signature genes. Interestingly, IL-18 was localized mainly to the basal epithelial cells and the infiltrating macrophages in the regressing VP, whereas IL-15 was induced in the luminal epithelium. The VP castration model exhibits immune cell infiltration and loss of PTEN that is often observed in progressive PCa, thereby making this model useful for further delineation of androgen-regulated gene expression with relevance to PCa.

Identification of a novel putative non-selenocysteine containing phospholipid hydroperoxide glutathione peroxidase (NPGPx) essential for alleviating oxidative stress generated from polyunsaturated fatty acids in breast cancer cells.

A dramatic reduction in the expression of a novel phospholipid hydroperoxide glutathione peroxidase (PHGPx), which incorporates cysteine instead of selenocysteine in the conserved catalytic motif was observed in a microarray analysis using cDNAs amplified from mRNA of Brca1-null mouse embryonic fibroblasts. This non-selenocysteine PHGPx named NPGPx is a cytoplasmic protein with molecular mass of approximately 22 kDa and has little detectable glutathione peroxidase activity in vitro. Ectopic expression of NPGPx in Brca1-null cells that were sensitive to oxidative stress induced by hydrogen peroxide conferred a similar resistance level to that of the wild-type cells, suggesting the importance of this protein in reducing oxidative stress. Expression of NPGPx was found in many tissues, including developing mammary gland. However, the majority of breast cancer cell lines studied (11 of 12) expressed very low or undetectable levels of NPGPx irrespective of BRCA1 status. Re-expression of NPGPx in breast cancer lines, MCF-7 and HCC1937, which have very little or no endogenous NPGPx, induced resistance to eicosapentaenoic acid (an omega-3 type of polyunsaturated fatty acid)-mediated cell death. Conversely, inhibition of the expression of NPGPx by the specific small interfering RNA in HS578T breast cancer cells that originally express substantial amounts of endogenous NPGPx increased their sensitivity to eicosapentaenoic acid-mediated cell death. Thus, NPGPx plays an essential role in breast cancer cells in alleviating oxidative stress generated from polyunsaturated fatty acid metabolism.

Application of gene expression profiling for validating models of human breast cancer.

While classical histopathologic approaches are invaluable in classifying tumors and understanding aspects of cellular interactions, genomic approaches provide a means to molecularly dissect tumorigenesis. The relationship of gene expression to the development of neoplasia remains an area of intensive research. With the advent of large-scale genomic platforms, alterations in gene expression can be related to the morphological development of cancer. The feasibility of using large-scale genomic analysis platforms has dramatically changed the landscape of biological sciences, as cellular processes must be considered in the context of complex networks. Alterations in gene expression must now be understood in a systems approach in which the relationships between genes expression changes are studied by considering the interplay of multiple regulatory networks. Ultimately, such changes must be understood at the protein level. We have begun to apply this technology to determine changes in gene expression that differentiate various types of mammary cancers that arise in mouse models that have been initiated by different genetic alterations. Ultimately, a molecular catalogue of similarities and differences between rodent and human tumors can be created which will serve to validate or credential particular models for specific experimental purposes, such as preclinical testing. These approaches have led to new insights into molecular pathways involved in oncogenesis, new classifications of human breast cancer, and the identification of new genes that may be relevant to understanding and treating human cancer.

Comparative functional analysis of rat TGF-beta1 and Xenopus laevis TGF-beta5 promoters suggest differential regulations.

We have carried out a comparative functional analysis of the rat TGF-beta1 and Xenopus laevis TGF-beta5 promoters across several mammalian and amphibian cell lines. Progressive deletion constructs of both the promoters have been made using a PCR based approach and the basal promoter activities studied in Xenopus tadpole cell line (XTC), Xenopus adult kidney fibroblast cell line (A6), human hepatoma cell line (HepG2), normal rat kidney cell line (NRK), and Chinese hamster ovary cell line (CHO). Data suggests that the basal promoter activity of TGF-beta1 is low as compared to TGF-beta5 promoter in XTC cells but comparable in A6 cells, while TGF-beta5 promoter shows nearly negligible activity as compared to TGF-beta5 promoter in all the tested mammalian cell lines. Moreover, TGF-beta5 promoter is found to be repressed in XTC cells on treatment with TGF-beta5 protein. Thus, the regulation of TGF-beta1 and TGF-beta5 promoters is distinct in amphibian and mammalian species. We therefore suggest that contrary to the suggested functional equivalence of TGF-beta1 and TGF-beta5 proteins, TGF-beta1 and TGF-beta5 genes have distinct functions in their respective species.

Interrogating mouse mammary cancer models: insights from gene expression profiling.

Numerous mouse models for mammary cancer have been developed and characterized based upon their biological, molecular, and histopathological features. In an effort to dissect the molecular anatomy of such models and compare their gene expression profiles to those of human breast cancer, six models representing various oncogenic pathways have been investigated using cDNA microarray technology. Results of these analyses are presented and discussed in the context of technological challenges presented by analyzing data on such a large scale. Further expression profiling coupled with emerging proteomic technologies will more completely define and distinguish mouse models of mammary cancer from each other and provide a comprehensive basis for comparing such models with the human disease they are intended to represent.

Initiating oncogenic event determines gene-expression patterns of human breast cancer models.

Molecular expression profiling of tumors initiated by transgenic overexpression of c-myc, c-neu, c-ha-ras, polyoma middle T antigen (PyMT) or simian virus 40 T/t antigen (T-ag) targeted to the mouse mammary gland have identified both common and oncogene-specific events associated with tumor formation and progression. The tumors shared great similarities in their gene-expression profiles as compared with the normal mammary gland with an induction of cell-cycle regulators, metabolic regulators, zinc finger proteins, and protein tyrosine phosphatases, along with the suppression of some protein tyrosine kinases. Selection and hierarchical clustering of the most variant genes, however, resulted in separating the mouse models into three groups with distinct oncogene-specific patterns of gene expression. Such an identification of targets specified by particular oncogenes may facilitate development of lesion-specific therapeutics and preclinical testing. Moreover, similarities in gene expression between human breast cancers and the mouse models have been identified, thus providing an important component for the validation of transgenic mammary cancer models.