RNA Interactions Are Essential for CTCF-Mediated Genome Organization.

The function of the CCCTC-binding factor (CTCF) in the organization of the genome has become an important area of investigation, but the mechanisms by which CTCF dynamically contributes to genome organization are not clear. We previously discovered that CTCF binds to large numbers of endogenous RNAs, promoting its self-association. In this regard, we now report two independent features that disrupt CTCF association with chromatin: inhibition of transcription and disruption of CTCF-RNA interactions through mutations of 2 of its 11 zinc fingers that are not required for CTCF binding to its cognate DNA site: zinc finger 1 (ZF1) or zinc finger 10 (ZF10). These mutations alter gene expression profiles as CTCF mutants lose their ability to form chromatin loops and thus the ability to insulate chromatin domains and to mediate CTCF long-range genomic interactions. Our results point to the importance of CTCF-mediated RNA interactions as a structural component of genome organization.

Gene expression dynamics during temperature-dependent sex determination in a sea turtle.

Fifty years ago, researchers discovered a link between ambient temperature and the sex of turtle embryos. More recently, significant progress has been made in understanding the influence of temperature on freshwater turtles. However, our understanding of the key genetic factors in other turtle groups, such as sea turtles, remains limited. To address this gap, we conducted RNA-seq analyses on embryonic tissues from the sea olive ridley turtle during the thermosensitive period (stages 21-26) at temperatures known to produce males (26 °C) and females (33 °C). Our findings revealed that incubation temperatures primarily influence genes with broad expression across tissues due to differential cell division rates and later have an effect regulating gonad-specific transcripts. This effect is mostly related to gene activation rather than transcription repression. We performed transcriptome analyses following shifts in incubation temperatures of bi-potential gonads. This approach allowed us to identify genes that respond rapidly and may be closer to the beginning of the temperature-sensing pathway. Notably, we observed swift adaptations in the expression levels of chromatin modifiers JARID2 and KDM6B, as well as the splicing factor SRSF5, and transcription regulators THOC2, DDX3X and CBX3, but little impact in the overall gonad-specific pathways, indicating that temperature-sensing genes may change rapidly but the rewiring of the gonad's developmental fate is complex and resilient. AUTHOR SUMMARY: Sea turtles, one of the most iconic creatures of our oceans, confront a troubling reality of endangerment, a peril magnified by the looming specter of climate change. This climatic shift is gradually increasing the temperature of the nesting beaches thus causing dramatic male/female population biases. Conservation efforts will need genetic and molecular information to reverse the negative effects of climate change on the populations. In this study, we conducted the first transcriptomic analysis of embryonic tissues, including gonads, brain, liver, and mesonephros, in the olive ridley sea turtle during the critical thermosensitive period spanning stages 21-26. We examined both male-producing (26 °C) and female-producing (33 °C) temperatures and found that incubation temperatures influence temperature-sensitive genes that are either expressed globally or specifically associated with the gonads. These findings indicate that incubation temperatures predominantly sway genes with broad expression patterns due to differential cell division rates. This natural process was opted in the gonads to drive sex determination. We also identified genes that are rapidly capable of sensing temperature changes and that could play a role in the activation of the sex determination pathway. Overall, our study sheds light on the intricate interplay between temperature and gene expression during sea turtle development, revealing dynamic changes in the transcriptome and highlighting the involvement of key genetic players in sex determination.

Chromatin Structure from Development to Ageing.

Nuclear structure influences genome architecture, which contributes to determine patterns of gene expression. Global changes in chromatin dynamics are essential during development and differentiation, and are one of the hallmarks of ageing. This chapter describes the molecular dynamics of chromatin structure that occur during development and ageing. In the first part, we introduce general information about the nuclear lamina, the chromatin structure, and the 3D organization of the genome. Next, we detail the molecular hallmarks found during development and ageing, including the role of DNA and histone modifications, 3D genome dynamics, and changes in the nuclear lamina. Within the chapter we discuss the implications that genome structure has on the mechanisms that drive development and ageing, and the physiological consequences when these mechanisms fail.

The pluripotent regulatory circuitry connecting promoters to their long-range interacting elements.

The mammalian genome harbors up to one million regulatory elements often located at great distances from their target genes. Long-range elements control genes through physical contact with promoters and can be recognized by the presence of specific histone modifications and transcription factor binding. Linking regulatory elements to specific promoters genome-wide is currently impeded by the limited resolution of high-throughput chromatin interaction assays. Here we apply a sequence capture approach to enrich Hi-C libraries for >22,000 annotated mouse promoters to identify statistically significant, long-range interactions at restriction fragment resolution, assigning long-range interacting elements to their target genes genome-wide in embryonic stem cells and fetal liver cells. The distal sites contacting active genes are enriched in active histone modifications and transcription factor occupancy, whereas inactive genes contact distal sites with repressive histone marks, demonstrating the regulatory potential of the distal elements identified. Furthermore, we find that coregulated genes cluster nonrandomly in spatial interaction networks correlated with their biological function and expression level. Interestingly, we find the strongest gene clustering in ES cells between transcription factor genes that control key developmental processes in embryogenesis. The results provide the first genome-wide catalog linking gene promoters to their long-range interacting elements and highlight the complex spatial regulatory circuitry controlling mammalian gene expression.

A novel cis-regulatory element regulates αD and αA-globin gene expression in chicken erythroid cells.

Background: Cis-regulatory elements (CREs) play crucial roles in regulating gene expression during erythroid cell differentiation. Genome-wide erythroid-specific CREs have not been characterized in chicken erythroid cells, which is an organism model used to study epigenetic regulation during erythropoiesis. Methods: Analysis of public genome-wide accessibility (ATAC-seq) maps, along with transcription factor (TF) motif analysis, CTCF, and RNA Pol II occupancy, as well as transcriptome analysis in fibroblasts and erythroid HD3 cells, were used to characterize erythroid-specific CREs. An α-globin CRE was identified, and its regulatory activity was validated in vitro and in vivo by luciferase activity and genome-editing assays in HD3 cells, respectively. Additionally, circular chromosome conformation capture (UMI-4C) assays were used to distinguish its role in structuring the α-globin domain in erythroid chicken cells. Results: Erythroid-specific CREs displayed occupancy by erythroid TF binding motifs, CTCF, and RNA Pol II, as well as an association with genes involved in hematopoiesis and cell differentiation. An α-globin CRE, referred to as CRE-2, was identified as exhibiting enhancer activity over αD and αA genes in vitro and in vivo. Induction of terminal erythroid differentiation showed that α-globin CRE-2 is required for the induction of αD and αA. Analysis of TF binding motifs at α-globin CRE-2 shows apparent regulation mediated by GATA-1, YY1, and CTCF binding. Conclusion: Our findings demonstrate that cell-specific CREs constitute a key mechanism that contributes to the fine-tuning gene regulation of erythroid cell differentiation and provide insights into the annotation and characterization of CREs in chicken cells.

Obtention of viable cell suspensions from breast cancer tumor biopsies for 3D chromatin conformation and single-cell transcriptome analysis.

Molecular and cellular characterization of tumors is essential due to the complex and heterogeneous nature of cancer. In recent decades, many bioinformatic tools and experimental techniques have been developed to achieve personalized characterization of tumors. However, sample handling continues to be a major challenge as limitations such as prior treatments before sample acquisition, the amount of tissue obtained, transportation, or the inability to process fresh samples pose a hurdle for experimental strategies that require viable cell suspensions. Here, we present an optimized protocol that allows the recovery of highly viable cell suspensions from breast cancer primary tumor biopsies. Using these cell suspensions we have successfully characterized genome architecture through Hi-C. Also, we have evaluated single-cell gene expression and the tumor cellular microenvironment through single-cell RNAseq. Both technologies are key in the detailed and personalized molecular characterization of tumor samples. The protocol described here is a cost-effective alternative to obtain viable cell suspensions from biopsies simply and efficiently.

An insulator embedded in the chicken α-globin locus regulates chromatin domain configuration and differential gene expression.

Genome organization into transcriptionally active domains denotes one of the first levels of gene expression regulation. Although the chromatin domain concept is generally accepted, only little is known on how domain organization impacts the regulation of differential gene expression. Insulators might hold answers to address this issue as they delimit and organize chromatin domains. We have previously identified a CTCF-dependent insulator with enhancer-blocking activity embedded in the 5' non-coding region of the chicken α-globin domain. Here, we demonstrate that this element, called the αEHS-1.4 insulator, protects a transgene against chromosomal position effects in stably transfected cell lines and transgenic mice. We found that this insulator can create a regulated chromatin environment that coincides with the onset of adult α-globin gene expression. Furthermore, such activity is in part dependent on the in vivo regulated occupancy of CTCF at the αEHS-1.4 element. Insulator function is also regulated by CTCF poly(ADP-ribosyl)ation. Our results suggest that the αEHS-1.4 insulator contributes in organizing the chromatin structure of the α-globin gene domain and prevents activation of adult α-globin gene expression at the erythroblast stage via CTCF.

Chromosome conformation capture technologies as tools to detect structural variations and their repercussion in chromatin 3D configuration.

3D genome organization regulates gene expression in different physiological and pathological contexts. Characterization of chromatin structure at different scales has provided information about how the genome organizes in the nuclear space, from chromosome territories, compartments of euchromatin and heterochromatin, topologically associated domains to punctual chromatin loops between genomic regulatory elements and gene promoters. In recent years, chromosome conformation capture technologies have also been used to characterize structural variations (SVs) de novo in pathological conditions. The study of SVs in cancer, has brought information about transcriptional misregulation that relates directly to the incidence and prognosis of the disease. For example, gene fusions have been discovered arising from chromosomal translocations that upregulate oncogenes expression, and other types of SVs have been described that alter large genomic regions encompassing many genes. However, studying SVs in 2D cannot capture all their regulatory implications in the genome. Recently, several bioinformatic tools have been developed to identify and classify SVs from chromosome conformation capture data and clarify how they impact chromatin structure in 3D, resulting in transcriptional misregulation. Here, we review recent literature concerning bioinformatic tools to characterize SVs from chromosome conformation capture technologies and exemplify their vast potential to rebuild the 3D landscape of genomes in cancer. The study of SVs from the 3D perspective can produce essential information about drivers, molecular targets, and disease evolution.

RNA polymerase II pausing regulates chromatin organization in erythrocytes.

Chicken erythrocytes are nucleated cells often considered to be transcriptionally inactive, although the epigenetic changes and chromatin remodeling that would mediate transcriptional repression and the extent of gene silencing during avian terminal erythroid differentiation are not fully understood. Here, we characterize the changes in gene expression, chromatin accessibility, genome organization and chromatin nuclear disposition during the terminal stages of erythropoiesis in chicken and uncover complex chromatin reorganization at different genomic scales. We observe a robust decrease in transcription in erythrocytes, but a set of genes maintains their expression, including genes involved in RNA polymerase II (Pol II) promoter-proximal pausing. Erythrocytes exhibit a reoriented nuclear architecture, with accessible chromatin positioned towards the nuclear periphery together with the paused RNA Pol II. In erythrocytes, chromatin domains are partially lost genome-wide, except at minidomains retained around paused promoters. Our results suggest that promoter-proximal pausing of RNA Pol II contributes to the transcriptional regulation of the erythroid genome and highlight the role of RNA polymerase in the maintenance of local chromatin organization.

Down-Regulation-Resistant STAT4 Risk Haplotype Contributes to Lupus Nephritis Through CD4+ T Cell Interferon-γ Production.

OBJECTIVE: Variants in STAT4 are associated with systemic lupus erythematosus (SLE) and other autoimmune diseases. We undertook this study to investigate how disease-associated variants affect STAT4 expression, in particular in CD4+ T cells where STAT4 plays an essential role. METHODS: We compared Th1 differentiation between naive CD4+ T cells from healthy donors homozygous for the risk (R/R) or nonrisk (NR/NR) alleles. We analyzed epigenetic marks in STAT4 and evaluated the relevance of its third intron, assessed the consequences of Stat4 overexpression in vivo in mice, and analyzed the effects of the STAT4 genotype in patients with lupus nephritis. RESULTS: Naive CD4+ T cells from NR/NR healthy donors down-regulated STAT4 in response to interleukin-12 (IL-12). In contrast, cells from R/R healthy donors maintained high levels. R/R cells exhibited a higher abundance of transcriptionally active STAT4 and increased interferon-γ production. Accordingly, R/R healthy donors exhibited a stronger induction of local active enhancer marks. Genetic editing confirmed the presence of a negative regulatory region in the STAT4 third intron, where most of the SLE-associated STAT4 single-nucleotide polymorphisms (SNPs) are located. In vivo forced expression demonstrated that increases in Stat4 levels in T cells enhanced glomerulonephritis in mice. Accordingly, the R/R genotype was associated with suboptimal response to treatment and with worse clinical outcomes in patients with proliferative lupus nephritis. CONCLUSION: The SLE-associated STAT4 haplotype correlates with an abnormal IL-12-mediated STAT4 transcriptional regulation. Carriers of the risk variant exhibit exaggerated CD4+ proinflammatory capacities that, in the context of SLE, contribute to more severe disease. R/R patients may benefit from blockade of the IL-12/STAT4 pathway.

Gain of DNA methylation is enhanced in the absence of CTCF at the human retinoblastoma gene promoter.

BACKGROUND: Long-term gene silencing throughout cell division is generally achieved by DNA methylation and other epigenetic processes. Aberrant DNA methylation is now widely recognized to be associated with cancer and other human diseases. Here we addressed the contribution of the multifunctional nuclear factor CTCF to the epigenetic regulation of the human retinoblastoma (Rb) gene promoter in different tumoral cell lines. METHODS: To assess the DNA methylation status of the Rb promoter, genomic DNA from stably transfected human erythroleukemic K562 cells expressing a GFP reporter transgene was transformed with sodium bisulfite, and then PCR-amplified with modified primers and sequenced. Single- and multi-copy integrants with the CTCF binding site mutated were isolated and characterized by Southern blotting. Silenced transgenes were reactivated using 5-aza-2'-deoxycytidine and Trichostatin-A, and their expression was monitored by fluorescent cytometry. Rb gene expression and protein abundance were assessed by RT-PCR and Western blotting in three different glioma cell lines, and DNA methylation of the promoter region was determined by sodium bisulfite sequencing, together with CTCF dissociation and methyl-CpG-binding protein incorporation by chromatin immunoprecipitation assays. RESULTS: We found that the inability of CTCF to bind to the Rb promoter causes a dramatic loss of gene expression and a progressive gain of DNA methylation. CONCLUSIONS: This study indicates that CTCF plays an important role in maintaining the Rb promoter in an optimal chromatin configuration. The absence of CTCF induces a rapid epigenetic silencing through a progressive gain of DNA methylation. Consequently, CTCF can now be seen as one of the epigenetic components that allows the proper configuration of tumor suppressor gene promoters. Its aberrant dissociation can then predispose key genes in cancer cells to acquire DNA methylation and epigenetic silencing.

In-Nucleus Hi-C in Drosophila Cells.

The genome is organized into topologically associating domains (TADs) delimited by boundaries that isolate interactions between domains. In Drosophila, the mechanisms underlying TAD formation and boundaries are still under investigation. The application of the in-nucleus Hi-C method described here helped to dissect the function of architectural protein (AP)-binding sites at TAD boundaries isolating the Notch gene. Genetic modification of domain boundaries that cause loss of APs results in TAD fusion, transcriptional defects, and long-range topological alterations. These results provided evidence demonstrating the contribution of genetic elements to domain boundary formation and gene expression control in Drosophila. Here, the in-nucleus Hi-C method has been described in detail, which provides important checkpoints to assess the quality of the experiment along with the protocol. Also shown are the required numbers of sequencing reads and valid Hi-C pairs to analyze genomic interactions at different genomic scales. CRISPR/Cas9-mediated genetic editing of regulatory elements and high-resolution profiling of genomic interactions using this in-nucleus Hi-C protocol could be a powerful combination for the investigation of the structural function of genetic elements.

The global and promoter-centric 3D genome organization temporally resolved during a circadian cycle.

BACKGROUND: Circadian gene expression is essential for organisms to adjust their physiology and anticipate daily changes in the environment. The molecular mechanisms controlling circadian gene transcription are still under investigation. In particular, how chromatin conformation at different genomic scales and regulatory elements impact rhythmic gene expression has been poorly characterized. RESULTS: Here we measure changes in the spatial chromatin conformation in mouse liver using genome-wide and promoter-capture Hi-C alongside daily oscillations in gene transcription. We find topologically associating domains harboring circadian genes that switch assignments between the transcriptionally active and inactive compartment at different hours of the day, while their boundaries stably maintain their structure over time. To study chromatin contacts of promoters at high resolution over time, we apply promoter capture Hi-C. We find circadian gene promoters displayed a maximal number of chromatin contacts at the time of their peak transcriptional output. Furthermore, circadian genes, as well as contacted and transcribed regulatory elements, reach maximal expression at the same timepoints. Anchor sites of circadian gene promoter loops are enriched in DNA binding sites for liver nuclear receptors and other transcription factors, some exclusively present in either rhythmic or stable contacts. Finally, by comparing the interaction profiles between core clock and output circadian genes, we show that core clock interactomes are more dynamic compared to output circadian genes. CONCLUSION: Our results identify chromatin conformation dynamics at different scales that parallel oscillatory gene expression and characterize the repertoire of regulatory elements that control circadian gene transcription through rhythmic or stable chromatin configurations.

Heterochromatin as an Important Driver of Genome Organization.

Heterochromatin is a constituent of eukaryotic genomes with functions spanning from gene expression silencing to constraining DNA replication and repair. Inside the nucleus, heterochromatin segregates spatially from euchromatin and is localized preferentially toward the nuclear periphery and surrounding the nucleolus. Despite being an abundant nuclear compartment, little is known about how heterochromatin regulates and participates in the mechanisms driving genome organization. Here, we review pioneer and recent evidence that explores the functional role of heterochromatin in the formation of distinct chromatin compartments and how failure of the molecular mechanisms forming heterochromatin leads to disarray of genome conformation and disease.

In situ dissection of domain boundaries affect genome topology and gene transcription in Drosophila.

Chromosomes are organized into high-frequency chromatin interaction domains called topologically associating domains (TADs), which are separated from each other by domain boundaries. The molecular mechanisms responsible for TAD formation are not yet fully understood. In Drosophila, it has been proposed that transcription is fundamental for TAD organization while the participation of genetic sequences bound by architectural proteins (APs) remains controversial. Here, we investigate the contribution of domain boundaries to TAD organization and the regulation of gene expression at the Notch gene locus in Drosophila. We find that deletion of domain boundaries results in TAD fusion and long-range topological defects that are accompanied by loss of APs and RNA Pol II chromatin binding as well as defects in transcription. Together, our results provide compelling evidence of the contribution of discrete genetic sequences bound by APs and RNA Pol II in the partition of the genome into TADs and in the regulation of gene expression in Drosophila.

Sequential chromatin immunoprecipitation protocol: ChIP-reChIP.

Chromatin immunoprecipitation has been widely used to determine the status of histone covalent modifications and also to investigate DNA-protein and protein-protein associations to a particular genomic location in vivo. Generally, DNA regulatory elements nucleate the interaction of several transcription factors in conjunction with ubiquitous and/or tissue-specific cofactors in order to regulate gene transcription. Therefore, it has become relevant to determine the cohabitation of several proteins in a particular developmental stage and cell type. Furthermore, multiple post-translational histone modifications can be analyzed on the same genomic location with the aim of deciphering the combinatorial pattern of histone modifications associated to specific transcriptional stages during cell commitment. Here we describe the ChIP-reChIP assay that represents a direct strategy to determine the in vivo colocalization of proteins interacting or in close contact in a chromatinized template on the basis of double and independent rounds of immunoprecipitations with high-quality ChIP grade antibodies.

RNA proximity sequencing reveals the spatial organization of the transcriptome in the nucleus.

The global, three-dimensional organization of RNA molecules in the nucleus is difficult to determine using existing methods. Here we introduce Proximity RNA-seq, which identifies colocalization preferences for pairs or groups of nascent and fully transcribed RNAs in the nucleus. Proximity RNA-seq is based on massive-throughput RNA barcoding of subnuclear particles in water-in-oil emulsion droplets, followed by cDNA sequencing. Our results show RNAs of varying tissue-specificity of expression, speed of RNA polymerase elongation and extent of alternative splicing positioned at varying distances from nucleoli. The simultaneous detection of multiple RNAs in proximity to each other distinguishes RNA-dense from sparse compartments. Application of Proximity RNA-seq will facilitate study of the spatial organization of transcripts in the nucleus, including non-coding RNAs, and its functional relevance.

Shaping Up the Embryo: The Role of Genome 3D Organization.

The spatial organization of the chromatinized genome inside the cell nucleus impacts genomic function. In transcription, the hierarchical genome structure creates spatial regulatory landscapes, in which modulating elements like enhancers can contact their target genes and activate their expression, as a result of restricting their exploration to a specific topological neighbourhood. Here we describe exciting recent findings obtained through "C" technologies in pluripotent cells and early embryogenesis and emphasize some of the key unanswered questions arising from them.

Promoter Capture Hi-C: High-resolution, Genome-wide Profiling of Promoter Interactions.

The three-dimensional organization of the genome is linked to its function. For example, regulatory elements such as transcriptional enhancers control the spatio-temporal expression of their target genes through physical contact, often bridging considerable (in some cases hundreds of kilobases) genomic distances and bypassing nearby genes. The human genome harbors an estimated one million enhancers, the vast majority of which have unknown gene targets. Assigning distal regulatory regions to their target genes is thus crucial to understand gene expression control. We developed Promoter Capture Hi-C (PCHi-C) to enable the genome-wide detection of distal promoter-interacting regions (PIRs), for all promoters in a single experiment. In PCHi-C, highly complex Hi-C libraries are specifically enriched for promoter sequences through in-solution hybrid selection with thousands of biotinylated RNA baits complementary to the ends of all promoter-containing restriction fragments. The aim is to then pull-down promoter sequences and their frequent interaction partners such as enhancers and other potential regulatory elements. After high-throughput paired-end sequencing, a statistical test is applied to each promoter-ligated restriction fragment to identify significant PIRs at the restriction fragment level. We have used PCHi-C to generate an atlas of long-range promoter interactions in dozens of human and mouse cell types. These promoter interactome maps have contributed to a greater understanding of mammalian gene expression control by assigning putative regulatory regions to their target genes and revealing preferential spatial promoter-promoter interaction networks. This information also has high relevance to understanding human genetic disease and the identification of potential disease genes, by linking non-coding disease-associated sequence variants in or near control sequences to their target genes.