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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Long-Range Enhancer Interactions Are Prevalent in Mouse Embryonic Stem Cells and Are Reorganized upon Pluripotent State Transition.

Transcriptional enhancers, including super-enhancers (SEs), form physical interactions with promoters to regulate cell-type-specific gene expression. SEs are characterized by high transcription factor occupancy and large domains of active chromatin, and they are commonly assigned to target promoters using computational predictions. How promoter-SE interactions change upon cell state transitions, and whether transcription factors maintain SE interactions, have not been reported. Here, we used promoter-capture Hi-C to identify promoters that interact with SEs in mouse embryonic stem cells (ESCs). We found that SEs form complex, spatial networks in which individual SEs contact multiple promoters, and a rewiring of promoter-SE interactions occurs between pluripotent states. We also show that long-range promoter-SE interactions are more prevalent in ESCs than in epiblast stem cells (EpiSCs) or Nanog-deficient ESCs. We conclude that SEs form cell-type-specific interaction networks that are partly dependent on core transcription factors, thereby providing insights into the gene regulatory organization of pluripotent cells.

Lineage-specific dynamic and pre-established enhancer-promoter contacts cooperate in terminal differentiation.

Chromosome conformation is an important feature of metazoan gene regulation; however, enhancer-promoter contact remodeling during cellular differentiation remains poorly understood. To address this, genome-wide promoter capture Hi-C (CHi-C) was performed during epidermal differentiation. Two classes of enhancer-promoter contacts associated with differentiation-induced genes were identified. The first class ('gained') increased in contact strength during differentiation in concert with enhancer acquisition of the H3K27ac activation mark. The second class ('stable') were pre-established in undifferentiated cells, with enhancers constitutively marked by H3K27ac. The stable class was associated with the canonical conformation regulator cohesin, whereas the gained class was not, implying distinct mechanisms of contact formation and regulation. Analysis of stable enhancers identified a new, essential role for a constitutively expressed, lineage-restricted ETS-family transcription factor, EHF, in epidermal differentiation. Furthermore, neither class of contacts was observed in pluripotent cells, suggesting that lineage-specific chromatin structure is established in tissue progenitor cells and is further remodeled in terminal differentiation.

Global reorganisation of cis-regulatory units upon lineage commitment of human embryonic stem cells.

Long-range cis-regulatory elements such as enhancers coordinate cell-specific transcriptional programmes by engaging in DNA looping interactions with target promoters. Deciphering the interplay between the promoter connectivity and activity of cis-regulatory elements during lineage commitment is crucial for understanding developmental transcriptional control. Here, we use Promoter Capture Hi-C to generate a high-resolution atlas of chromosomal interactions involving ~22,000 gene promoters in human pluripotent and lineage-committed cells, identifying putative target genes for known and predicted enhancer elements. We reveal extensive dynamics of cis-regulatory contacts upon lineage commitment, including the acquisition and loss of promoter interactions. This spatial rewiring occurs preferentially with predicted changes in the activity of cis-regulatory elements and is associated with changes in target gene expression. Our results provide a global and integrated view of promoter interactome dynamics during lineage commitment of human pluripotent cells.

Polycomb repressive complex PRC1 spatially constrains the mouse embryonic stem cell genome.

The Polycomb repressive complexes PRC1 and PRC2 maintain embryonic stem cell (ESC) pluripotency by silencing lineage-specifying developmental regulator genes. Emerging evidence suggests that Polycomb complexes act through controlling spatial genome organization. We show that PRC1 functions as a master regulator of mouse ESC genome architecture by organizing genes in three-dimensional interaction networks. The strongest spatial network is composed of the four Hox gene clusters and early developmental transcription factor genes, the majority of which contact poised enhancers. Removal of Polycomb repression leads to disruption of promoter-promoter contacts in the Hox gene network. In contrast, promoter-enhancer contacts are maintained in the absence of Polycomb repression, with accompanying widespread acquisition of active chromatin signatures at network enhancers and pronounced transcriptional upregulation of network genes. Thus, PRC1 physically constrains developmental transcription factor genes and their enhancers in a silenced but poised spatial network. We propose that the selective release of genes from this spatial network underlies cell fate specification during early embryonic development.