Deli Hong

Chromatin interaction analysis reveals changes in small chromosome and telomere clustering between epithelial and breast cancer cells

BackgroundHigher-order chromatin structure is often perturbed in cancer and other pathological states. Although several genetic and epigenetic differences have been charted between normal and breast cancer tissues, changes in higher-order chromatin organization during tumorigenesis have not been fully explored. To probe the differences in higher-order chromatin structure between mammary epithelial and breast cancer cells, we performed Hi-C analysis on MCF-10A mammary epithelial and MCF-7 breast cancer cell lines.ResultsOur studies reveal that the small, gene-rich chromosomes chr16 through chr22 in the MCF-7 breast cancer genome display decreased interaction frequency with each other compared to the inter-chromosomal interaction frequency in the MCF-10A epithelial cells. Interestingly, this finding is associated with a higher occurrence of open compartments on chr16–22 in MCF-7 cells. Pathway analysis of the MCF-7 up-regulated genes located in altered compartment regions on chr16–22 reveals pathways related to repression of WNT signaling. There are also differences in intra-chromosomal interactions between the cell lines; telomeric and sub-telomeric regions in the MCF-10A cells display more frequent interactions than are observed in the MCF-7 cells.ConclusionsWe show evidence of an intricate relationship between chromosomal organization and gene expression between epithelial and breast cancer cells. Importantly, this work provides a genome-wide view of higher-order chromatin dynamics and a resource for studying higher-order chromatin interactions in two cell lines commonly used to study the progression of breast cancer.

hsa-mir-30c promotes the invasive phenotype of metastatic breast cancer cells by targeting NOV/CCN3.

BackgroundFor treatment and prevention of metastatic disease, one of the premier challenges is the identification of pathways and proteins to target for clinical intervention. Micro RNAs (miRNAs) are short, non-coding RNAs, which regulate cellular activities by either mRNA degradation or translational inhibition. Our studies focused on the invasive properties of hsa-mir30c based on its high expression in MDA-MB-231 metastatic cells and our bioinformatic analysis of the Cancer Genome Atlas that identified aberrant hsa-mir-30c to be associated with poor survival.MethodsContributions of hsa-mir-30c to breast cancer cell invasion were examined by Matrigel invasion transwell assays following modulation of hsa-mir-30c or hsa-mir-30c* levels in MDA-MB-231 cells. hsa-mir-30c in silico predicted targets linked to cell invasion were screened for targeting by hsa-mir-30c in metastatic breast cancer cells by RT-qPCR. The contribution to invasion by a target of hsa-mir-30c, Nephroblastoma overexpressed (NOV), was characterized by siRNA and invasion assays. Significant effects were determined using Student’s T-tests with Welch’s correction for unequal variance.ResultsMCF-7 and MDA-MB-231 cells were used as models of poorly invasive and late-stage metastatic disease, respectively. By modulating the levels of hsa-mir-30c in these cells, we observed concomitant changes in breast cancer cell invasiveness. From predicted targets of hsa-mir-30c that were related to cellular migration and invasion, NOV/CCN3 was identified as a novel target of hsa-mir-30c. Depleting NOV by siRNA caused a significant increase in the invasiveness of MDA-MB-231 cells is a regulatory protein associated with the extracellular matrix.ConclusionsNOV/CCN3 expression, which protects cells from invasion, is known in patient tumors to inversely correlate with advanced breast cancer and metastasis. This study has identified a novel target of hsa-mir-30c, NOV, which is an inhibitor of the invasiveness of metastatic breast cancer cells. Thus, hsa-mir-30c-mediated inhibition of NOV levels promotes the invasive phenotype of MDA-MB-231 cells and significantly, the miR-30/NOV pathways is independent of RUNX2, a known target of hsa-mir-30c that promotes osteolytic disease in metastatic breast cancer cells. Our findings allow for mechanistic insight into the clinical observation of poor survival of patients with elevated hsa-mir-30c levels, which can be considered for miRNA-based translational studies.

RUNX1 contributes to higher-order chromatin organization and gene regulation in breast cancer cells

RUNX1 is a transcription factor functioning both as an oncogene and a tumor suppressor in breast cancer. RUNX1 alters chromatin structure in cooperation with chromatin modifier and remodeling enzymes. In this study, we examined the relationship between RUNX1-mediated transcription and genome organization. We characterized genome-wide RUNX1 localization and performed RNA-seq and Hi-C in RUNX1-depleted and control MCF-7 breast cancer cells. RNA-seq analysis showed that RUNX1 depletion led to up-regulation of genes associated with chromatin structure and down-regulation of genes related to extracellular matrix biology, as well as NEAT1 and MALAT1 lncRNAs. Our ChIP-Seq analysis supports a prominent role for RUNX1 in transcriptional activation. About 30% of all RUNX1 binding sites were intergenic, indicating diverse roles in promoter and enhancer regulation and suggesting additional functions for RUNX1. Hi-C analysis of RUNX1-depleted cells demonstrated that overall three-dimensional genome organization is largely intact, but indicated enhanced association of RUNX1 near Topologically Associating Domain (TAD) boundaries and alterations in long-range interactions. These results suggest an architectural role for RUNX1 in fine-tuning local interactions rather than in global organization. Our results provide novel insight into RUNX1-mediated perturbations of higher-order genome organization that are functionally linked with RUNX1-dependent compromised gene expression in breast cancer cells.

MicroRNA-378-mediated suppression of Runx1 alleviates the aggressive phenotype of triple-negative MDA-MB-231 human breast cancer cells.

The Runx1 transcription factor, known for its essential role in normal hematopoiesis, was reported in limited studies to be mutated or associated with human breast tumor tissues. Runx1 increases concomitantly with disease progression in the MMTV-PyMT transgenic mouse model of breast cancer. Compelling questions relate to mechanisms that regulate Runx1 expression in breast cancer. Here, we tested the hypothesis that dysregulation of Runx1-targeting microRNAs (miRNAs) allows for pathologic increase of Runx1 during breast cancer progression. Microarray profiling of the MMTV-PyMT model revealed significant downregulation of numerous miRNAs predicted to target Runx1. One of these, miR-378, was inversely correlated with Runx1 expression during breast cancer progression in mice and in human breast cancer cell lines MCF7 and triple-negative MDA-MB-231 that represent early- and late-stage diseases, respectively. MiR-378 is nearly absent in MDA-MB-231 cells. Luciferase reporter assays revealed that miR-378 binds the Runx1 3′ untranslated region (3′UTR) and inhibits Runx1 expression. Functionally, we demonstrated that ectopic expression of miR-378 in MDA-MB-231 cells inhibited Runx1 and suppressed migration and invasion, while inhibition of miR-378 in MCF7 cells increased Runx1 levels and cell migration. Depletion of Runx1 in late-stage breast cancer cells resulted in increased expression of both the miR-378 host gene PPARGC1B and pre-miR-378, suggesting a feedback loop. Taken together, our study identifies a novel and clinically relevant mechanism for regulation of Runx1 in breast cancer that is mediated by a PPARGC1B-miR-378-Runx1 regulatory pathway. Our results highlight the translational potential of miRNA replacement therapy for inhibiting Runx1 in breast cancer.

Runx1 stabilizes the mammary epithelial cell phenotype and prevents epithelial to mesenchymal transition

Runx1 is a well characterized transcription factor essential for hematopoietic differentiation and Runx1 mutations are the cause of leukemias. Runx1 is highly expressed in normal epithelium of most glands and recently has been associated with solid tumors. Notably, the function of Runx1 in the mammary gland and how it is involved in initiation and progression of breast cancer is still unclear. Here we demonstrate the consequences of Runx1 loss in normal mammary epithelial and breast cancer cells. We first observed that Runx1 is decreased in tumorigenic and metastatic breast cancer cells. We also observed loss of Runx1 expression upon induction of epithelial-mesenchymal transition (EMT) in MCF10A (normal-like) cells. Furthermore depletion of Runx1 in MCF10A cells resulted in striking changes in cell shape, leading to mesenchymal cell morphology. The epithelial phenotype could be restored in breast cancer cells by re-expressing Runx1. Analyses of breast tumors and patient data revealed that low Runx1 expression is associated with poor prognosis and decreased survival. We addressed mechanisms for the function of Runx1 in maintaining the epithelial phenotype and find Runx1 directly regulates E-cadherin; and serves as a downstream transcription factor mediating TGFβ signaling. We also observed through global gene expression profiling of growth factor depleted cells that induction of EMT and loss of Runx1 is associated with activation of TGFβ and WNT pathways. Thus these findings have identified a novel function for Runx1 in sustaining normal epithelial morphology and preventing EMT and suggest Runx1 levels could be a prognostic indicator of tumor progression.

Identifying Nuclear Matrix-Attached DNA Across the Genome.

Experimental approaches to define the relationship between gene expression and nuclear matrix attachment regions (MARs) have given contrasting and method‐specific results. We have developed a next generation sequencing strategy to identify MARs across the human genome (MAR‐Seq). The method is based on crosslinking chromatin to its nuclear matrix attachment sites to minimize changes during biochemical processing. We used this method to compare nuclear matrix organization in MCF‐10A mammary epithelial‐like cells and MDA‐MB‐231 breast cancer cells and evaluated the results in the context of global gene expression (array analysis) and positional enrichment of gene‐regulatory histone modifications (ChIP‐Seq). In the normal‐like cells, nuclear matrix‐attached DNA was enriched in expressed genes, while in the breast cancer cells, it was enriched in non‐expressed genes. In both cell lines, the chromatin modifications that mark transcriptional activation or repression were appropriately associated with gene expression. Using this new MAR‐Seq approach, we provide the first genome‐wide characterization of nuclear matrix attachment in mammalian cells and reveal that the nuclear matrix‐associated genome is highly cell‐context dependent. J. Cell. Physiol. 232: 1295–1305, 2017.

Intranuclear and higher-order chromatin organization of the major histone gene cluster in breast cancer†

Alterations in nuclear morphology are common in cancer progression. However, the degree to which gross morphological abnormalities translate into compromised higher‐order chromatin organization is poorly understood. To explore the functional links between gene expression and chromatin structure in breast cancer, we performed RNA‐seq gene expression analysis on the basal breast cancer progression model based on human MCF10A cells. Positional gene enrichment identified the major histone gene cluster at chromosome 6p22 as one of the most significantly upregulated (and not amplified) clusters of genes from the normal‐like MCF10A to premalignant MCF10AT1 and metastatic MCF10CA1a cells. This cluster is subdivided into three sub‐clusters of histone genes that are organized into hierarchical topologically associating domains (TADs). Interestingly, the sub‐clusters of histone genes are located at TAD boundaries and interact more frequently with each other than the regions in‐between them, suggesting that the histone sub‐clusters form an active chromatin hub. The anchor sites of loops within this hub are occupied by CTCF, a known chromatin organizer. These histone genes are transcribed and processed at a specific sub‐nuclear microenvironment termed the major histone locus body (HLB). While the overall chromatin structure of the major HLB is maintained across breast cancer progression, we detected alterations in its structure that may relate to gene expression. Importantly, breast tumor specimens also exhibit a coordinate pattern of upregulation across the major histone gene cluster. Our results provide a novel insight into the connection between the higher‐order chromatin organization of the major HLB and its regulation during breast cancer progression.

Expression of the IL‐11 Gene in Metastatic Cells Is Supported by Runx2‐Smad and Runx2‐cJun Complexes Induced by TGFβ1

In tumor cells, two factors are abnormally increased that contribute to metastatic bone disease: Runx2, a transcription factor that promotes expression of metastasis related and osteolytic genes; and IL‐11, a secreted osteolytic cytokine. Here, we addressed a compelling question: Does Runx2 regulate IL‐11 gene expression? We find a positive correlation between Runx2, IL‐11 and TGFβ1, a driver of the vicious cycle of metastatic bone disease, in prostate cancer (PC) cell lines representing early (LNCaP) and late (PC3) stage disease. Further, like Runx2 knockdown, IL‐11 knockdown significantly reduced expression of several osteolytic factors. Modulation of Runx2 expression results in corresponding changes in IL‐11 expression. The IL‐11 gene has Runx2, AP‐1 sites and Smad binding elements located on the IL‐11 promoter. Here, we demonstrated that Runx2‐c‐Jun as well as Runx2‐Smad complexes upregulate IL‐11 expression. Functional studies identified a significant loss of IL‐11 expression in PC3 cells in the presence of the Runx2‐HTY mutant protein, a mutation that disrupts Runx2‐Smad signaling. In response to TGFβ1 and in the presence of Runx2, we observed a 30‐fold induction of IL‐11 expression, accompanied by increased c‐Jun binding to the IL‐11 promoter. Immunoprecipitation and in situ co‐localization studies demonstrated that Runx2 and c‐Jun form nuclear complexes in PC3 cells. Thus, TGFβ1 signaling induces two independent transcriptional pathways ‐ AP‐1 and Runx2. These transcriptional activators converge on IL‐11 as a result of Runx2‐Smad and Runx2‐c‐Jun interactions to amplify IL‐11 gene expression that, together with Runx2, supports the osteolytic pathology of cancer induced bone disease. J. Cell. Biochem. 116: 2098–2108, 2015.

Bivalent Epigenetic Control of Oncofetal Gene Expression in Cancer

ABSTRACT Multiple mechanisms of epigenetic control that include DNA methylation, histone modification, noncoding RNAs, and mitotic gene bookmarking play pivotal roles in stringent gene regulation during lineage commitment and maintenance. Experimental evidence indicates that bivalent chromatin domains, i.e., genome regions that are marked by both H3K4me3 (activating) and H3K27me3 (repressive) histone modifications, are a key property of pluripotent stem cells. Bivalency of developmental genes during the G1 phase of the pluripotent stem cell cycle contributes to cell fate decisions. Recently, some cancer types have been shown to exhibit partial recapitulation of bivalent chromatin modifications that are lost along with pluripotency, suggesting a mechanism by which cancer cells reacquire properties that are characteristic of undifferentiated, multipotent cells. This bivalent epigenetic control of oncofetal gene expression in cancer cells may offer novel insights into the onset and progression of cancer and may provide specific and selective options for diagnosis as well as for therapeutic intervention.

RUNX1 suppresses breast cancer stemness and tumor growth

Recent studies have revealed that mutations in the transcription factor Runx1 are prevalent in breast tumors. Yet, how loss of Runx1 contributes to breast cancer (BCa) remains unresolved. We demonstrate for the first time that Runx1 represses the breast cancer stem cell (BCSC) phenotype and consequently, functions as a tumor suppressor in breast cancer. Runx1 ectopic expression in MCF10AT1 and MCF10CA1a BCa cells reduces (60%) migration, invasion and in vivo tumor growth in mouse mammary fat pad (P<0.05). Runx1 is decreased in BCSCs, and overexpression of Runx1 suppresses tumorsphere formation and reduces the BCSC population. Furthermore, Runx1 inhibits Zeb1 expression, while Runx1 depletion activates Zeb1 and the epithelial-mesenchymal transition. Mechanistically Runx1 functions as a tumor suppressor in breast cancer through repression of cancer stem cell activity. This key regulation of BCSCs by Runx1 may be shared in other epithelial carcinomas, highlighting the importance of Runx1 in solid tumors.