Spatiotemporal Epigenetic Control of the Histone Gene Chromatin Landscape during the Cell Cycle.

Higher-order genomic organization supports the activation of histone genes in response to cell cycle regulatory cues that epigenetically mediates stringent control of transcription at the G1/S-phase transition. Histone locus bodies (HLBs) are dynamic, non-membranous, phase-separated nuclear domains where the regulatory machinery for histone gene expression is organized and assembled to support spatiotemporal epigenetic control of histone genes. HLBs provide molecular hubs that support synthesis and processing of DNA replication-dependent histone mRNAs. These regulatory microenvironments support long-range genomic interactions among non-contiguous histone genes within a single topologically associating domain (TAD). HLBs respond to activation of the cyclin E/CDK2/NPAT/HINFP pathway at the G1/S transition. HINFP and its coactivator NPAT form a complex within HLBs that controls histone mRNA transcription to support histone protein synthesis and packaging of newly replicated DNA. Loss of HINFP compromises H4 gene expression and chromatin formation, which may result in DNA damage and impede cell cycle progression. HLBs provide a paradigm for higher-order genomic organization of a subnuclear domain that executes an obligatory cell cycle-controlled function in response to cyclin E/CDK2 signaling. Understanding the coordinately and spatiotemporally organized regulatory programs in focally defined nuclear domains provides insight into molecular infrastructure for responsiveness to cell signaling pathways that mediate biological control of growth, differentiation phenotype, and are compromised in cancer.

Identification of tRNA-derived small RNA (tsRNA) responsive to the tumor suppressor, RUNX1, in breast cancer.

Despite recent advances in targeted therapies, the molecular mechanisms driving breast cancer initiation, progression, and metastasis are minimally understood. Growing evidence indicate that transfer RNA (tRNA)-derived small RNAs (tsRNA) contribute to biological control and aberrations associated with cancer development and progression. The runt-related transcription factor 1 (RUNX1) transcription factor is a tumor suppressor in the mammary epithelium whereas RUNX1 downregulation is functionally associated with breast cancer initiation and progression. We identified four tsRNA (ts-19, ts-29, ts-46, and ts-112) that are selectively responsive to expression of the RUNX1 tumor suppressor. Our finding that ts-112 and RUNX1 anticorrelate in normal-like mammary epithelial and breast cancer lines is consistent with tumor-related activity of ts-112 and tumor suppressor activity of RUNX1. Inhibition of ts-112 in MCF10CA1a aggressive breast cancer cells significantly reduced proliferation. Ectopic expression of a ts-112 mimic in normal-like mammary epithelial MCF10A cells significantly increased proliferation. These findings support an oncogenic potential for ts-112. Moreover, RUNX1 may repress ts-112 to prevent overactive proliferation in breast epithelial cells to augment its established roles in maintaining the mammary epithelium.

RUNX1 and RUNX2 transcription factors function in opposing roles to regulate breast cancer stem cells.

Breast cancer stem cells (BCSCs) are competent to initiate tumor formation and growth and refractory to conventional therapies. Consequently BCSCs are implicated in tumor recurrence. Many signaling cascades associated with BCSCs are critical for epithelial-to-mesenchymal transition (EMT). We developed a model system to mechanistically examine BCSCs in basal-like breast cancer using MCF10AT1 FACS sorted for CD24 (negative/low in BCSCs) and CD44 (positive/high in BCSCs). Ingenuity Pathway Analysis comparing RNA-seq on the CD24-/low versus CD24+/high MCF10AT1 indicates that the top activated upstream regulators include TWIST1, TGFß1, OCT4, and other factors known to be increased in BCSCs and during EMT. The top inhibited upstream regulators include ESR1, TP63, and FAS. Consistent with our results, many genes previously demonstrated to be regulated by RUNX factors are altered in BCSCs. The RUNX2 interaction network is the top significant pathway altered between CD24-/low and CD24+/high MCF10AT1. RUNX1 is higher in expression at the RNA level than RUNX2. RUNX3 is not expressed. While, human-specific quantitative polymerase chain reaction primers demonstrate that RUNX1 and CDH1 decrease in human MCF10CA1a cells that have grown tumors within the murine mammary fat pad microenvironment, RUNX2 and VIM increase. Treatment with an inhibitor of RUNX binding to CBFß for 5 days followed by a 7-day recovery period results in EMT suggesting that loss of RUNX1, rather than increase in RUNX2, is a driver of EMT in early stage breast cancer. Increased understanding of RUNX regulation on BCSCs and EMT will provide novel insight into therapeutic strategies to prevent recurrence.

Chromosome territories and the global regulation of the genome.

Spatial positioning is a fundamental principle governing nuclear processes. Chromatin is organized as a hierarchy from nucleosomes to Mbp chromatin domains (CD) or topologically associating domains (TADs) to higher level compartments culminating in chromosome territories (CT). Microscopic and sequencing techniques have substantiated chromatin organization as a critical factor regulating gene expression. For example, enhancers loop back to interact with their target genes almost exclusively within TADs, distally located coregulated genes reposition into common transcription factories upon activation, and Mbp CDs exhibit dynamic motion and configurational changes in vivo. A longstanding question in the nucleus field is whether an interactive nuclear matrix provides a direct link between structure and function. The findings of nonrandom radial positioning of CT within the nucleus suggest the possibility of preferential interaction patterns among populations of CT. Sequential labeling up to 10 CT followed by application of computer imaging and geometric graph mining algorithms revealed cell-type specific interchromosomal networks (ICN) of CT that are altered during the cell cycle, differentiation, and cancer progression. It is proposed that the ICN correlate with the global level of genome regulation. These approaches also demonstrated that the large scale 3-D topology of CT is specific for each CT. The cell-type specific proximity of certain chromosomal regions in normal cells may explain the propensity of distinct translocations in cancer subtypes. Understanding how genes are dysregulated upon disruption of the normal "wiring" of the nucleus by translocations, deletions, and amplifications that are hallmarks of cancer, should enable more targeted therapeutic strategies.

RUNX1-dependent mechanisms in biological control and dysregulation in cancer.

The RUNX1 transcription factor has recently been shown to be obligatory for normal development. RUNX1 controls the expression of genes essential for proper development in many cell lineages and tissues including blood, bone, cartilage, hair follicles, and mammary glands. Compromised RUNX1 regulation is associated with many cancers. In this review, we highlight evidence for RUNX1 control in both invertebrate and mammalian development and recent novel findings of perturbed RUNX1 control in breast cancer that has implications for other solid tumors. As RUNX1 is essential for definitive hematopoiesis, RUNX1 mutations in hematopoietic lineage cells have been implicated in the etiology of several leukemias. Studies of solid tumors have revealed a context-dependent function for RUNX1 either as an oncogene or a tumor suppressor. These RUNX1 functions have been reported for breast, prostate, lung, and skin cancers that are related to cancer subtypes and different stages of tumor development. Growing evidence suggests that RUNX1 suppresses aggressiveness in most breast cancer subtypes particularly in the early stage of tumorigenesis. Several studies have identified RUNX1 suppression of the breast cancer epithelial-to-mesenchymal transition. Most recently, RUNX1 repression of cancer stem cells and tumorsphere formation was reported for breast cancer. It is anticipated that these new discoveries of the context-dependent diversity of RUNX1 functions will lead to innovative therapeutic strategies for the intervention of cancer and other abnormalities of normal tissues.

SMARCA4 regulates gene expression and higher-order chromatin structure in proliferating mammary epithelial cells.

The packaging of DNA into chromatin plays an important role in transcriptional regulation and nuclear processes. Brahma-related gene-1 SMARCA4 (also known as BRG1), the essential ATPase subunit of the mammalian SWI/SNF chromatin remodeling complex, uses the energy from ATP hydrolysis to disrupt nucleosomes at target regions. Although the transcriptional role of SMARCA4 at gene promoters is well-studied, less is known about its role in higher-order genome organization. SMARCA4 knockdown in human mammary epithelial MCF-10A cells resulted in 176 up-regulated genes, including many related to lipid and calcium metabolism, and 1292 down-regulated genes, some of which encode extracellular matrix (ECM) components that can exert mechanical forces and affect nuclear structure. ChIP-seq analysis of SMARCA4 localization and SMARCA4-bound super-enhancers demonstrated extensive binding at intergenic regions. Furthermore, Hi-C analysis showed extensive SMARCA4-mediated alterations in higher-order genome organization at multiple resolutions. First, SMARCA4 knockdown resulted in clustering of intra- and inter-subtelomeric regions, demonstrating a novel role for SMARCA4 in telomere organization. SMARCA4 binding was enriched at topologically associating domain (TAD) boundaries, and SMARCA4 knockdown resulted in weakening of TAD boundary strength. Taken together, these findings provide a dynamic view of SMARCA4-dependent changes in higher-order chromatin organization and gene expression, identifying SMARCA4 as a novel component of chromatin organization.

Epithelial-to-mesenchymal transition and cancer stem cells contribute to breast cancer heterogeneity.

Breast cancer is the most common cancer in women, and accounts for ~30% of new cancer cases and 15% of cancer-related deaths. Tumor relapse and metastasis are primary factors contributing to breast cancer-related deaths. Therefore, the challenge for breast cancer treatment is to sustain remission. A driving force behind tumor relapse is breast cancer heterogeneity (both intertumor, between different patients, and intratumor, within the same tumor). Understanding breast cancer heterogeneity is necessary to develop preventive interventions and targeted therapies. A recently emerging concept is that intratumor heterogeneity is driven by cancer stem cells (CSCs) that are capable of giving rise to a multitude of different cells within a tumor. Studies have highlighted linkage of CSC formation with epithelial-to-mesenchymal transition (EMT). In this review, we summarize the current understanding of breast cancer heterogeneity, links between EMT and CSCs, regulation of EMT by Runx transcription factors, and potential therapeutic strategies targeting these processes.

Selective expression of long non-coding RNAs in a breast cancer cell progression model.

Long non-coding RNAs (lncRNAs) are acknowledged as regulators of cancer biology and pathology. Our goal was to perform a stringent profiling of breast cancer cell lines that represent disease progression. We used the MCF-10 series, which includes the normal-like MCF-10A, HRAS-transformed MCF-10AT1 (pre-malignant), and MCF-10CA1a (malignant) cells, to perform transcriptome wide sequencing. From these data, we have identified 346 lncRNAs with dysregulated expression across the progression series. By comparing lncRNAs from these datasets to those from an additional set of cell lines that represent different disease stages and subtypes, MCF-7 (early stage, luminal), and MDA-MB-231 (late stage, basal), 61 lncRNAs that are associated with breast cancer progression were identified. Querying breast cancer patient data from The Cancer Genome Atlas, we selected a lncRNA, IGF-like family member 2 antisense RNA 1 (IGFL2-AS1), of potential clinical relevance for functional characterization. Among the 61 lncRNAs, IGFL2-AS1 was the most significantly decreased. Our results indicate that this lncRNA plays a role in downregulating its nearest neighbor, IGFL1, and affects migration of breast cancer cells. Furthermore, the lncRNAs we identified provide a valuable resource to mechanistically and clinically understand the contribution of lncRNAs in breast cancer progression.

Higher order genomic organization and regulatory compartmentalization for cell cycle control at the G1/S-phase transition.

Fidelity of histone gene regulation, and ultimately of histone protein biosynthesis, is obligatory for packaging of newly replicated DNA into chromatin. Control of histone gene expression within the 3-dimensional context of nuclear organization is reflected by two well documented observations. DNA replication-dependent histone mRNAs are synthesized at specialized subnuclear domains designated histone locus bodies (HLBs), in response to activation of the growth factor dependent Cyclin E/CDK2/HINFP/NPAT pathway at the G1/S transition in mammalian cells. Complete loss of the histone gene regulatory factors HINFP or NPAT disrupts HLB integrity that is necessary for coordinate control of DNA replication and histone gene transcription. Here we review the molecular histone-related requirements for G1/S-phase progression during the cell cycle. Recently developed experimental strategies, now enable us to explore mechanisms involved in dynamic control of histone gene expression in the context of the temporal (cell cycle) and spatial (HLBs) remodeling of the histone gene loci.

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.

Large-scale probabilistic 3D organization of human chromosome territories.

There is growing evidence that chromosome territories (CT) have a probabilistic non-random arrangement within the cell nucleus of mammalian cells including radial positioning and preferred patterns of interchromosomal interactions that are cell-type specific. While it is generally assumed that the three-dimensional (3D) arrangement of genes within the CT is linked to genomic regulation, the degree of non-random organization of individual CT remains unclear. As a first step to elucidating the global 3D organization (topology) of individual CT, we performed multi-color fluorescence in situ hybridization using six probes extending across each chromosome in human WI38 lung fibroblasts. Six CT were selected ranging in size and gene density (1, 4, 12, 17, 18 and X). In-house computational geometric algorithms were applied to measure the 3D distances between every combination of probes and to elucidate data-mined structural patterns. Our findings demonstrate a high degree of non-random arrangement of individual CT that vary from chromosome to chromosome and display distinct changes during the cell cycle. Application of a classic, well-defined data mining and pattern recognition approach termed the 'k-means' generated 3D models for the best fit arrangement of each chromosome. These predicted models correlated well with the detailed distance measurements and analysis. We propose that the unique 3D topology of each CT and characteristic changes during the cell cycle provide the structural framework for the global gene expression programs of the individual chromosomes.

C-ing the Genome: A Compendium of Chromosome Conformation Capture Methods to Study Higher-Order Chromatin Organization.

Three-dimensional organization of the chromatin has important roles in transcription, replication, DNA repair, and pathologic events such as translocations. There are two fundamental ways to study higher-order chromatin organization: microscopic and molecular approaches. In this review, we briefly introduce the molecular approaches, focusing on chromosome conformation capture or "3C" technology and its derivatives, which can be used to probe chromatin folding at resolutions beyond that provided by microscopy techniques. We further discuss the different types of data generated by the 3C-based methods and how they can be used to answer distinct biological questions.

Chromosomes at Work: Organization of Chromosome Territories in the Interphase Nucleus.

The organization of interphase chromosomes in chromosome territories (CTs) was first proposed more than one hundred years ago. The introduction of increasingly sophisticated microscopic and molecular techniques, now provide complementary strategies for studying CTs in greater depth than ever before. Here we provide an overview of these strategies and how they are being used to elucidate CT interactions and the role of these dynamically regulated, nuclear-structure building blocks in directly supporting nuclear function in a physiologically responsive manner.

Wide-scale alterations in interchromosomal organization in breast cancer cells: defining a network of interacting chromosomes.

The interchromosomal spatial positionings of a subset of human chromosomes was examined in the human breast cell line MCF10A (10A) and its malignant counterpart MCF10CA1a (CA1a). The nine chromosomes selected (#1, 4, 11, 12, 15, 16, 18, 21 and X) cover a wide range in size and gene density and compose ∼40% of the total human genome. Radial positioning of the chromosome territories (CT) was size dependent with certain of the CT more peripheral in CA1a. Each CT was in close proximity (interaction) with a similar number of other CT except the inactive CTXi. It had lower levels of interchromosomal partners in 10A which increased strikingly in CA1a. Major alterations from 10A to CA1a were detected in the pairwise interaction profiles which were subdivided into five types of altered interaction profiles: overall increase, overall decrease, switching from 1 to ≥2, vice versa or no change. A global data mining program termed the chromatic median calculated the most probable overall association network for the entire subset of CT. This interchromosomal network was drastically altered in CA1a with only 1 of 20 shared connections. We conclude that CT undergo multiple and preferred interactions with other CT in the cell nucleus and form preferred-albeit probabilistic-interchromosomal networks. This network of interactions is highly altered in malignant human breast cells. It is intriguing to consider the relationship of these alterations to the corresponding changes in the gene expression program of these malignant cancer cells.

Non-Random Patterns in the Distribution of NOR-Bearing Chromosome Territories in Human Fibroblasts: A Network Model of Interactions.

We present a 3-D mapping in WI38 human diploid fibroblast cells of chromosome territories (CT) 13,14,15,21, and 22, which contain the nucleolar organizing regions (NOR) and participate in the formation of nucleoli. The nuclear radial positioning of NOR-CT correlated with the size of chromosomes with smaller CT more interior. A high frequency of pairwise associations between NOR-CT ranging from 52% (CT13-21) to 82% (CT15-21) was detected as well as a triplet arrangement of CT15-21-22 (72%). The associations of homologous CT were significantly lower (24-36%). Moreover, singular contacts between CT13-14 or CT13-22 were found in the majority of cells, while CT13-15 or CT13-21 predominantly exhibited multiple interactions. In cells with multiple nucleoli, one of the nucleoli (termed "dominant") always associated with a higher number of CT. Moreover, certain CT pairs more frequently contributed to the same nucleolus than to others. This nonrandom pattern suggests that a large number of the NOR-chromosomes are poised in close proximity during the postmitotic nucleolar recovery and through their NORs may contribute to the formation of the same nucleolus. A global data mining program termed the chromatic median determined the most probable interchromosomal arrangement of the entire NOR-CT population. This interactive network model was significantly above randomized simulation and was composed of 13 connections among the NOR-CT. We conclude that the NOR-CT form a global interactive network in the cell nucleus that may be a fundamental feature for the regulation of nucleolar and other genomic functions.

Gene density and chromosome territory shape.

Despite decades of study of chromosome territories (CT) in the interphase nucleus of mammalian cells, our understanding of the global shape and 3-D organization of the individual CT remains very limited. Past microscopic analysis of CT suggested that while many of the CT appear to be very regular ellipsoid-like shapes, there were also those with more irregular shapes. We have undertaken a comprehensive analysis to determine the degree of shape regularity of different CT. To be representative of the whole human genome, 12 different CT (~41 % of the genome) were selected that ranged from the largest (CT 1) to the smallest (CT 21) in size and from the highest (CT 19) to lowest (CT Y) in gene density. Using both visual inspection and algorithms that measure the degree of shape ellipticity and regularity, we demonstrate a strong inverse correlation between the degree of regular CT shape and gene density for those CT that are most gene-rich (19, 17, 11) and gene-poor (18, 13, Y). CT more intermediate in gene density showed a strong negative correlation with shape regularity, but not with ellipticity. An even more striking correlation between gene density and CT shape was determined for the nucleolar-associated NOR-CT. Correspondingly, striking differences in shape between the X active and inactive CT implied that aside from gene density, the overall global level of gene transcription on individual CT is also an important determinant of chromosome territory shape.

Cell type specific alterations in interchromosomal networks across the cell cycle.

The interchromosomal organization of a subset of human chromosomes (#1, 4, 11, 12, 16, 17, and 18) was examined in G1 and S phase of human WI38 lung fibroblast and MCF10A breast epithelial cells. Radial positioning of the chromosome territories (CTs) was independent of gene density, but size dependent. While no changes in radial positioning during the cell cycle were detected, there were stage-specific differences between cell types. Each CT was in close proximity (interaction) with a similar number of other CT except the gene rich CT17 which had significantly more interactions. Furthermore, CT17 was a member of the highest pairwise CT combinations with multiple interactions. Major differences were detected in the pairwise interaction profiles of MCF10A versus WI38 including cell cycle alterations from G1 to S. These alterations in interaction profiles were subdivided into five types: overall increase, overall decrease, switching from 1 to ≥2 interactions, vice versa, or no change. A global data mining program termed the chromatic median determined the most probable overall association network for the entire subset of CT. This probabilistic interchromosomal network was nearly completely different between the two cell lines. It was also strikingly altered across the cell cycle in MCF10A, but only slightly in WI38. We conclude that CT undergo multiple and preferred interactions with other CT in the nucleus and form preferred -albeit probabilistic- interchromosomal networks. This network of interactions is altered across the cell cycle and between cell types. It is intriguing to consider the relationship of these alterations to the corresponding changes in the gene expression program across the cell cycle and in different cell types.

Spatiotemporal higher-order chromatin landscape of human histone gene clusters at histone locus bodies during the cell cycle in breast cancer progression.

Human Histone Locus Bodies (HLBs) are nuclear subdomains comprised of clustered histone genes that are coordinately regulated throughout the cell cycle. We addressed temporal-spatial higher-order genome organization for time-dependent chromatin remodeling at HLBs that supports control of cell proliferation. Proximity distances of specific genomic contacts within histone gene clusters exhibit subtle changes during the G1 phase in MCF10 breast cancer progression model cell lines. This approach directly demonstrates that the two principal histone gene regulatory proteins, HINFP (H4 gene regulator) and NPAT, localize at chromatin loop anchor-points, denoted by CTCF binding, supporting the stringent requirement for histone biosynthesis to package newly replicated DNA as chromatin. We identified a novel enhancer region located âˆ¼ 2 MB distal to histone gene sub-clusters on chromosome 6 that consistently makes genomic contacts with HLB chromatin and is bound by NPAT. During G1 progression the first DNA loops form between one of three histone gene sub-clusters bound by HINFP and the distal enhancer region. Our findings are consistent with a model that the HINFP/NPAT complex controls the formation and dynamic remodeling of higher-order genomic organization of histone gene clusters at HLBs in early to late G1 phase to support transcription of histone mRNAs in S phase.

Epigenetic-Mediated Regulation of Gene Expression for Biological Control and Cancer: Fidelity of Mechanisms Governing the Cell Cycle.

The cell cycle is governed by stringent epigenetic mechanisms that, in response to intrinsic and extrinsic regulatory cues, support fidelity of DNA replication and cell division. We will focus on (1) the complex and interdependent processes that are obligatory for control of proliferation and compromised in cancer, (2) epigenetic and topological domains that are associated with distinct phases of the cell cycle that may be altered in cancer initiation and progression, and (3) the requirement for mitotic bookmarking to maintain intranuclear localization of transcriptional regulatory machinery to reinforce cell identity throughout the cell cycle to prevent malignant transformation.

Epigenetic-Mediated Regulation of Gene Expression for Biological Control and Cancer: Cell and Tissue Structure, Function, and Phenotype.

Epigenetic gene regulatory mechanisms play a central role in the biological control of cell and tissue structure, function, and phenotype. Identification of epigenetic dysregulation in cancer provides mechanistic into tumor initiation and progression and may prove valuable for a variety of clinical applications. We present an overview of epigenetically driven mechanisms that are obligatory for physiological regulation and parameters of epigenetic control that are modified in tumor cells. The interrelationship between nuclear structure and function is not mutually exclusive but synergistic. We explore concepts influencing the maintenance of chromatin structures, including phase separation, recognition signals, factors that mediate enhancer-promoter looping, and insulation and how these are altered during the cell cycle and in cancer. Understanding how these processes are altered in cancer provides a potential for advancing capabilities for the diagnosis and identification of novel therapeutic targets.