Three-Dimensional Genome Organization in Breast and Gynecological Cancers: How Chromatin Folding Influences Tumorigenic Transcriptional Programs

A growing body of research on the transcriptome and cancer genome has demonstrated that many gynecological tumor-specific gene mutations are located in cis-regulatory elements. Through chromosomal looping, cis-regulatory elements interact which each other to control gene expression by bringing distant regulatory elements, such as enhancers and insulators, into close proximity with promoters. It is well known that chromatin connections may be disrupted in cancer cells, promoting transcriptional dysregulation and the expression of abnormal tumor suppressor genes and oncogenes. In this review, we examine the roles of alterations in 3D chromatin interactions. This includes changes in CTCF protein function, cancer-risk single nucleotide polymorphisms, viral integration, and hormonal response as part of the mechanisms that lead to the acquisition of enhancers or super-enhancers. The translocation of existing enhancers, as well as enhancer loss or acquisition of insulator elements that interact with gene promoters, is also revised. Remarkably, similar processes that modify 3D chromatin contacts in gene promoters may also influence the expression of non-coding RNAs, such as long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), which have emerged as key regulators of gene expression in a variety of cancers, including gynecological malignancies.

Dynamic spreading of chromatin-mediated gene silencing and reactivation between neighboring genes in single cells

In mammalian cells genes that are in close proximity are coupled transcriptionally: silencing or activating one gene can affect its neighbors. Understanding these dynamics is important for natural processes, such as heterochromatin spreading during development and aging, and when designing synthetic gene regulation. Here, we systematically dissect this process in single cells by recruiting and releasing repressive chromatin regulators at dual-gene synthetic reporters, and measuring how fast gene silencing and reactivation spread as a function of intergenic distance and configuration of insulator elements. We find that silencing by KRAB, associated with histone methylation, spreads between two genes within hours, with a time delay that increases with distance. This fast KRAB-mediated spreading is not blocked by the classical cHS4 insulators. Silencing by histone deacetylase HDAC4 of the upstream gene can also lead to downstream gene silencing, but with a days-long delay that does not change with distance. This slower silencing can sometimes be stopped by insulators. Gene reactivation of neighboring genes is also coupled, with strong promoters and insulators determining the order of reactivation. We propose a new model of multi-gene regulation, where both gene silencing and gene reactivation can act at a distance, allowing for coordinated dynamics via chromatin regulator recruitment.

Sequential in cis mutagenesis in vivo reveals various functions for CTCF sites at the mouse HoxD cluster

Mammalian Hox gene clusters contain a range of CTCF binding sites. In addition to their importance in organizing a TAD border, which isolates the most posterior genes from the rest of the cluster, the positions and orientations of these sites suggest that CTCF may be instrumental in the selection of various subsets of contiguous genes, which are targets of distinct remote enhancers located in the flanking regulatory landscapes. We examined this possibility by producing an allelic series of cumulative in cis mutations in these sites, up to the abrogation of CTCF binding in the five sites located on one side of the TAD border. In the most impactful alleles, the global chromatin architecture of the locus was modified, yet not drastically, illustrating that CTCF sites located on one side of a strong TAD border are sufficient to organize at least part of this insulation. Spatial colinearity in the expression of these genes along the major body axis was nevertheless maintained, despite abnormal expression boundaries. In contrast, strong effects were scored in the selection of target genes responding to particular enhancers, leading to the misregulation of Hoxd genes in specific structures. Altogether, while most enhancer–promoter interactions can occur in the absence of this series of CTCF sites, the binding of CTCF in the Hox cluster is required to properly transform a rather unprecise process into a highly discriminative mechanism of interactions, which is translated into various patterns of transcription accompanied by the distinctive chromatin topology found at this locus. Our allelic series also allowed us to reveal the distinct functional contributions for CTCF sites within this Hox cluster, some acting as insulator elements, others being necessary to anchor or stabilize enhancer–promoter interactions, and some doing both, whereas they all together contribute to the formation of a TAD border. This variety of tasks may explain the amazing evolutionary conservation in the distribution of these sites among paralogous Hox clusters or between various vertebrates.

Sequential in-cis mutagenesis in vivo reveals various functions for CTCF sites at the mouse HoxD cluster

Mammalian Hox gene clusters contain a range of CTCF binding sites. In addition to their importance in organizing a TAD border, which isolates the most posterior genes from the rest of the cluster, the positions and orientations of these sites suggest that CTCF may be instrumental in the selection of various subsets of contiguous genes, which are targets of distinct remote enhancers located in the flanking regulatory landscapes. We examined this possibility by producing an allelic series of cumulative in-cis mutations in these sites, up to the abrogation of CTCF binding in the five sites located on one side of the TAD border. In the most impactful alleles, the global chromatin architecture of the locus was modified, yet not drastically, illustrating that CTCF sites located on one side of a strong TAD border are sufficient to organize at least part of this insulation. Spatial colinearity in the expression of these genes along the major body axis was nevertheless maintained, despite abnormal expression boundaries. In contrast, strong effects were scored in the selection of target genes responding to particular enhancers, leading to the mis-regulation of Hoxd genes in specific structures. Altogether, while most enhancer-promoter interactions can occur in the absence of this series of CTCF sites, it seems that the binding of CTCF in the Hox cluster is required to properly transform a rather unprecise process into a highly discriminative mechanism of interactions, which is translated into various patterns of transcription accompanied by the distinctive chromatin topology found at this locus. Our allelic series also allowed us to reveal the distinct functional contributions for CTCF sites within this Hox cluster, some acting as insulator elements, others being necessary to anchor or stabilize enhancer-promoter interactions and some doing both, whereas all together contribute to the formation of a TAD border. This variety of tasks may explain the amazing evolutionary conservation in the distribution of these sites amongst paralogous Hox clusters or between various vertebrates.

Boundaries support specific long-distance interactions between enhancers and promoters in Drosophila Bithorax complex

Drosophila bithorax complex (BX-C) is one of the best model systems for studying the role of boundaries (insulators) in gene regulation. Expression of three homeotic genes, Ubx, abd-A, and Abd-B, is orchestrated by nine parasegment-specific regulatory domains. These domains are flanked by boundary elements, which function to block crosstalk between adjacent domains, ensuring that they can act autonomously. Paradoxically, seven of the BX-C regulatory domains are separated from their gene target by at least one boundary, and must “jump over” the intervening boundaries. To understand the jumping mechanism, the Mcp boundary was replaced with Fab-7 and Fab-8. Mcp is located between the iab-4 and iab-5 domains, and defines the border between the set of regulatory domains controlling abd-A and Abd-B. When Mcp is replaced by Fab-7 or Fab-8, they direct the iab-4 domain (which regulates abd-A) to inappropriately activate Abd-B in abdominal segment A4. For the Fab-8 replacement, ectopic induction was only observed when it was inserted in the same orientation as the endogenous Fab-8 boundary. A similar orientation dependence for bypass activity was observed when Fab-7 was replaced by Fab-8. Thus, boundaries perform two opposite functions in the context of BX-C – they block crosstalk between neighboring regulatory domains, but at the same time actively facilitate long distance communication between the regulatory domains and their respective target genes.Author SummaryDrosophila bithorax complex (BX-C) is one of a few examples demonstrating in vivo role of boundary/insulator elements in organization of independent chromatin domains. BX-C contains three HOX genes, whose parasegment-specific pattern is controlled by cis-regulatory domains flanked by boundary/insulator elements. Since the boundaries ensure autonomy of adjacent domains, the presence of these elements poses a paradox: how do the domains bypass the intervening boundaries and contact their proper regulatory targets? According to the textbook model, BX-C regulatory domains are able to bypass boundaries because they harbor special promoter targeting sequences. However, contrary to this model, we show here that the boundaries themselves play an active role in directing regulatory domains to their appropriate HOX gene promoter.

The remote transcriptional control of Hox genes

Since the discovery by Ed Lewis that the order of Hox genes on the chromosome reflects the partitioning of their patterning function along the anterior-posterior axis of the developing fruit fly embryo, extensive efforts have been dedicated to uncovering the regulatory events underlying the collinear expression of Hox genes. These studies have revealed various aspects of Hox regulation, including short-range and long-range transcriptional enhancers, insulator elements and non-coding RNAs. With the development of technologies allowing for high resolution probing of chromatin architecture, notably Chromosome Conformation Capture (3C)-based techniques, a clear relationship is emerging between long-range regulation of Hox genes and the three-dimensional organization of the genome. Here, we provide an overview of these studies and in particular we discuss the functional relevance of genome compartmentalization, CTCF- mediated insulation and the Polycomb Repressive Complexes in the remote control of Hox genes.

Convergence of topological domain boundaries, insulators, and polytene interbands revealed by high-resolution mapping of chromatin contacts in the early Drosophila melanogaster embryo

High-throughput assays of three-dimensional interactions of chromosomes have shed considerable light on the structure of animal chromatin. Despite this progress, the precise physical nature of observed structures and the forces that govern their establishment remain poorly understood. Here we present high resolution Hi-C data from early Drosophila embryos. We demonstrate that boundaries between topological domains of various sizes map to DNA elements that resemble classical insulator elements: short genomic regions sensitive to DNase digestion that are strongly bound by known insulator proteins and are frequently located between divergent promoters. Further, we show a striking correspondence between these elements and the locations of mapped polytene interband regions. We believe it is likely this relationship between insulators, topological boundaries, and polytene interbands extends across the genome, and we therefore propose a model in which decompaction of boundary-insulator-interband regions drives the organization of interphase chromosomes by creating stable physical separation between adjacent domains.

A chromatin extension model for insulator function based on comparison of high-resolution chromatin conformation capture and polytene banding maps

Insulator proteins bind to specific genomic loci and have been shown to play a role in partitioning genomes into independent domains of gene expression and chromatin structure. Despite decades of study, the mechanism by which insulators establish these domains remains elusive. Here, we use genome-wide chromatin conformation capture (Hi-C) to generate a high-resolution map of spatial interactions of chromatin from Drosophila melanogaster embryos. We show that from the earliest stages of development the genome is divided into distinct topologically associated domains (TADs), that we can map the boundaries between TADs to sub-kilobase resolution, and that these boundaries correspond to 500-2000 bp insulator elements. Comparing this map with a detailed assessment of the banding pattern of a region of a polytene chromosome, we show that these insulator elements correspond to low density polytene interbands that divide compacted bands, which correspond to TADs. It has been previously shown that polytene interbands have low packing ratios allowing the conversion of small genomic distances (in base pairs) into a large physical distances. We therefore suggest a simple mechanism for insulator function whereby insulators increase the physical space between adjacent domains via the unpacking and extension of intervening chromatin. This model provides an intuitive explanation for known features of insulators, including the ability to block enhancer-promoter interactions, limit the spread of heterochromatin, and organize the structural features of interphase chromosomes.

Motif Disruption Domains Lead To Cancer Gene Expression Rewiring

Somatic mutations accumulate in non-coding regions of the genome during tumorigenesis, but their functional characterization presents a challenge. Somatic non-coding mutations rarely overlap among patients, which necessitates large sample sizes to detect associations. We analysed somatic mutations called from whole-genome sequencing (WGS) and RNA sequencing (RNAseq) from 3000 tumors across the Pan-Cancer Analysis of Whole Genomes to identify and functionally characterize mutation accumulation and its impact on gene dysregulation in cancer. We identified 1.5 million motif disruption domains (MDDs) across 40 cancer types, which we characterized as pan-cancer targets for recurrent mutation accumulation. These MDDs deregulate gene expression in cancer-specific and pan-cancer patterns by disrupting transcription factor binding sites in regulatory and insulator elements. Disruption is most recurrent across individuals at MDDs in conserved open chromatin, revealing potential drivers. This accumulation of somatic variants targeting regulatory and structural elements in MDDs generates gene expression dysregulation during tumorigenesis.

Targeted insertion in well-characterized Drosophila cell lines using φC31 integrase

We describe an adaptation of φC31 integrase-mediated targeted cassette exchange for use in Drosophila cell lines. Single copies of an attP-bounded docking platform carrying a GFP-expression marker, with and without insulator elements flanking the attP sites, were inserted by P-element transformation into the Kc167 and Sg4 cell lines; each of the resulting docking site lines carries a single mapped copy of one of the docking platforms. Vectors for targeted substitution contain a cloning cassette flanked by attB sites. Targeted substitution occurs by integrase-mediated substitution between the attP sites (integrated) and the attB sites (vector). We describe procedures for isolating cells carrying the substitutions and for eliminating the products of secondary off-target events. We demonstrated the technology by integrating a cassette containing a Cu++-inducible mCherry marker, and we report the expression properties of those lines. When compared with clonal lines made by traditional transformation methods, which lead to the illegitimate insertion of tandem arrays, targeted insertion lines give more uniform expression, lower basal expression and higher induction ratios. Targeted substitution, though intricate, affords results that should greatly improve comparative expression assays ??? a major emphasis of cell-based studies.