Transcription Elongation Can Affect Genome 3D Structure

An Alternatively Translated Isoform of PPARG Suggests AF-1 Domain Inhibition as an Insulin Sensitization Target

ARTICLE HIGHLIGHTS

Genetic screens were performed across PPARG to study how disruptive mutations across the full coding sequence affect function. An alternative translational start site in PPARG generates a truncated isoform, peroxisome proliferator-activated receptor γ (PPARγ) M135, which lacks the N-terminal activation function 1 (AF-1) domain and shows increased agonist-induced transactivation of target genes. In human carriers of rare PPARG variants, AF-1 domain-disrupting genetic variants increase agonist-induced PPARγ activity and decrease metabolic syndrome severity. Targeting the AF-1 domain is a potential therapeutic strategy for insulin sensitization.

Recipes and ingredients for deep learning models of 3D genome folding

Three-dimensional genome folding plays roles in gene regulation and disease. In this review, we compare and contrast recent deep learning models for predicting genome contact maps. We survey preprocessing, architecture, training, evaluation, and interpretation methods, highlighting the capabilities and limitations of different models. In each area, we highlight challenges, opportunities, and potential future directions for genome-folding models.

Human long noncoding RNA VILMIR is induced by major respiratory viral infections and modulates the host interferon response

Long noncoding RNAs (lncRNAs) are a newer class of noncoding transcripts identified as key regulators of biological processes. Here, we aimed to identify novel lncRNA targets that play critical roles in major human respiratory viral infections by systematically mining large-scale transcriptomic data sets. Using bulk RNA-sequencing (RNA-seq) analysis, we identified a previously uncharacterized lncRNA, named virus-inducible lncRNA modulator of interferon response (VILMIR), that was consistently upregulated after in vitro influenza infection across multiple human epithelial cell lines and influenza A virus subtypes. VILMIR was also upregulated after severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and respiratory syncytial virus (RSV) infections in vitro. We experimentally confirmed the response of VILMIR to influenza infection and interferon-beta (IFN-β) treatment in the A549 human epithelial cell line and found the expression of VILMIR was robustly induced by IFN-β treatment in a dose- and time-specific manner. Single-cell RNA-seq analysis of bronchoalveolar lavage fluid samples from coronavirus disease 2019 (COVID-19) patients uncovered that VILMIR was upregulated across various cell types, including at least five immune cells. The upregulation of VILMIR in immune cells was further confirmed in the human T cell and monocyte cell lines, SUP-T1 and THP-1, after IFN-β treatment. Finally, we found that knockdown of VILMIR expression reduced the magnitude of host transcriptional responses to both IFN-β treatment and influenza A virus infection in A549 cells. Together, our results show that VILMIR is a novel interferon-stimulated gene (ISG) that regulates the host interferon response and may be a potential therapeutic target for human respiratory viral infections upon further mechanistic investigation.IMPORTANCEIdentifying host factors that regulate the immune response to human respiratory viral infection is critical to developing new therapeutics. Human long noncoding RNAs (lncRNAs) have been found to play key regulatory roles during biological processes; however, the majority of lncRNA functions within the host antiviral response remain unknown. In this study, we identified that a previously uncharacterized lncRNA, virus-inducible lncRNA modulator of interferon response (VILMIR), is upregulated after major respiratory viral infections including influenza, severe acute respiratory syndrome coronavirus 2, and respiratory syncytial virus. We demonstrated that VILMIR is an interferon-stimulated gene that is upregulated after interferon-beta (IFN-β) in several human cell types. We also found that knockdown of VILMIR reduced the magnitude of host transcriptional responses to IFN-β treatment and influenza A infection in human epithelial cells. Our results reveal that VILMIR regulates the host interferon response and may present a new therapeutic target during human respiratory viral infections.

Bridging spatial and temporal scales of developmental gene regulation

The development of multicellular organisms relies on the precise coordination of molecular events across multiple spatial and temporal scales. Understanding how information flows from molecular interactions to cellular processes and tissue organization during development is crucial for explaining the remarkable reproducibility of complex organisms. This review explores how chromatin-encoded information is transduced from localized transcriptional events to global gene expression patterns, highlighting the challenge of bridging these scales. We discuss recent experimental findings and theoretical frameworks, emphasizing polymer physics as a tool for describing the relationship between chromatin structure and dynamics across scales. By integrating these perspectives, we aim to clarify how gene regulation is coordinated across levels of biological organization and suggest strategies for future experimental approaches.

A comprehensive review and benchmark of differential analysis tools for Hi-C data

MOTIVATION

The 3D organization of the genome plays a crucial role in various biological processes. Hi-C technology is widely used to investigate chromosome structures by quantifying 3D proximity between genomic regions. While numerous computational tools exist for detecting differences in Hi-C data between conditions, a comprehensive review and benchmark comparing their effectiveness is lacking.

RESULTS

This study offers a comprehensive review and benchmark of 10 generic tools for differential analysis of Hi-C matrices at the interaction count level. The benchmark assesses the statistical methods, usability, and performance (in terms of precision and power) of these tools, using both real and simulated Hi-C data. Results reveal a striking variability in performance among the tools, highlighting the substantial impact of preprocessing filters and the difficulty all tools encounter in effectively controlling the false discovery rate across varying resolutions and chromosome sizes.

AVAILABILITY

The complete benchmark is available at https://forgemia.inra.fr/scales/replication-chrocodiff using processed data deposited at https://doi.org/10.57745/LR0W9R.

CONTACT

nathalie.vialaneix@inrae.fr.

Reconstruction of diploid higher-order human 3D genome interactions from noisy Pore-C data using Dip3D

Differential high-order chromatin interactions between homologous chromosomes affect many biological processes. Traditional chromatin conformation capture genome analysis methods mainly identify two-way interactions and cannot provide comprehensive haplotype information, especially for low-heterozygosity organisms such as human. Here, we present a pipeline of methods to delineate diploid high-order chromatin interactions from noisy Pore-C outputs. We trained a previously published single-nucleotide variant (SNV)-calling deep learning model, Clair3, on Pore-C data to achieve superior SNV calling, applied a filtering strategy to tag reads for haplotypes and established a haplotype imputation strategy for high-order concatemers. Learning the haplotype characteristics of high-order concatemers from high-heterozygosity mouse allowed us to devise a progressive haplotype imputation strategy, which improved the haplotype-informative Pore-C contact rate 14.1-fold to 76% in the HG001 cell line. Overall, the diploid three-dimensional (3D) genome interactions we derived using Dip3D surpassed conventional methods in noise reduction and contact distribution uniformity, with better haplotype-informative contact density and genomic coverage rates. Dip3D identified previously unresolved haplotype high-order interactions, in addition to an understanding of their relationship with allele-specific expression, such as in X-chromosome inactivation. These results lead us to conclude that Dip3D is a robust pipeline for the high-quality reconstruction of diploid high-order 3D genome interactions.

lncRNA CARINH regulates expression and function of innate immune transcription factor IRF1 in macrophages

The discovery of long non-coding RNAs (lncRNAs) has provided a new perspective on the centrality of RNA in gene regulation and genome organization. Here, we screened for lncRNAs with putative functions in the host response to single-stranded RNA respiratory viruses. We identify CARINH as a conserved cis-acting lncRNA up-regulated in three respiratory diseases to control the expression of its antisense gene IRF1, a key transcriptional regulator of the antiviral response. CARINH and IRF1 are coordinately increased in the circulation of patients infected with human metapneumovirus, influenza A virus, or SARS-CoV-2, and in macrophages in response to viral infection or TLR3 agonist treatment. Targeted depletion of CARINH or its mouse ortholog Carinh in macrophages reduces the expression of IRF1/Irf1 and their associated target gene networks, increasing susceptibility to viral infection. Accordingly, CRISPR-mediated deletion of Carinh in mice reduces antiviral immunity, increasing viral burden upon sublethal challenge with influenza A virus. Together, these findings identify a conserved role of lncRNA CARINH in coordinating interferon-stimulated genes and antiviral immune responses.

Polymer models of chromatin organization in virally infected cells

Genome architecture is closely tied to essential biological functions, yet a complete understanding of the mechanisms governing DNA folding remains a significant challenge. Theoretical models based on polymer physics have been applied to decipher the complexity of chromatin architecture and uncover the physical processes shaping its structure. Importantly, recent findings suggest that certain viruses can alter the 3D organization of the host genome. In this review, we highlight recent advances in the development of polymer models used to study how chromatin 3D structure within a cell re-organizes following viral infection, with a particular emphasis on the SARS-CoV-2 virus, capable of altering genome organization of the host cell at different scales, including A/B compartments, TADs and gene-enhancer regulatory contacts.

Synchronized long-read genome, methylome, epigenome and transcriptome profiling resolve a Mendelian condition

Resolving the molecular basis of a Mendelian condition remains challenging owing to the diverse mechanisms by which genetic variants cause disease. To address this, we developed a synchronized long-read genome, methylome, epigenome and transcriptome sequencing approach, which enables accurate single-nucleotide, insertion-deletion and structural variant calling and diploid de novo genome assembly. This permits the simultaneous elucidation of haplotype-resolved CpG methylation, chromatin accessibility and full-length transcript information in a single long-read sequencing run. Application of this approach to an Undiagnosed Diseases Network participant with a chromosome X;13-balanced translocation of uncertain significance revealed that this translocation disrupted the functioning of four separate genes (NBEA, PDK3, MAB21L1 and RB1) previously associated with single-gene Mendelian conditions. Notably, the function of each gene was disrupted via a distinct mechanism that required integration of the four 'omes' to resolve. These included fusion transcript formation, enhancer adoption, transcriptional readthrough silencing and inappropriate X-chromosome inactivation of autosomal genes. Overall, this highlights the utility of synchronized long-read multi-omic profiling for mechanistically resolving complex phenotypes.

Evolution of 3D Chromatin Folding

Studies examining the evolution of genomes have focused mainly on sequence conservation. However, the inner working of a cell implies tightly regulated crosstalk between complex gene networks controlled by small dispersed regulatory elements of physically contacting DNA regions. How these different levels of chromatin organization crosstalk in different species underpins the potential for genome evolutionary plasticity. We review the evolution of chromatin organization across the Animal Tree of Life. We introduce general aspects of the mode and tempo of genome evolution to later explore the multiple layers of genome organization. We argue that both genome and chromosome size modulate patterns of chromatin folding and that chromatin interactions facilitate the formation of lineage-specific chromosomal reorganizations, especially in germ cells. Overall, analyzing the mechanistic forces involved in the maintenance of chromatin structure and function of the germ line is critical for understanding genome evolution, maintenance, and inheritance.

Viral oncogene EBNALP regulates YY1 DNA binding and alters host 3D genome organization

The Epstein-Barr virus (EBV) nuclear antigen leader protein (EBNALP) is essential for the immortalization of naive B lymphocytes (NBLs). However, the mechanisms remain elusive. To understand EBNALP's role in B-cell transformation, we compare NBLs infected with wild-type EBV and an EBNALP-null mutant EBV using multi-omics techniques. EBNALP inactivation alters enhancer-promoter interactions, resulting in decreased CCND2 and increased CASP1 and BCL2L11 expression. Mechanistically, EBNALP interacts with and colocalizes with the looping factor YY1. Depletion of EBNALP reduces YY1 DNA-binding and enhancer-promoter interactions, similar to effects observed with YY1 depletion. Furthermore, EBNALP colocalizes with DPF2, a protein that binds to H3K14ac and H4K16ac. CRISPR depletion of DPF2 reduces both EBNALP and YY1 DNA binding, suggesting that the DPF2/EBNALP complex may tether YY1 to DNA to increase enhancer-promoter interactions. EBNALP inactivation also increases enhancer-promoter interactions at the CASP1 and BCL2L11 loci, along with elevated DPF2 and YY1 binding and DNA accessibility. Our data suggest that EBNALP regulates YY1 to rewire the host genome, which might facilitate naive B-cell transformation.

RTP4 restricts influenza A virus infection by targeting the viral NS1 protein

The influenza A virus evades the host innate immune response to establish infection by inhibiting RIG-I activation through its nonstructural protein 1 (NS1). Here, we reported that receptor-transporting protein 4 (RTP4), an interferon-stimulated gene (ISG), targets NS1 to inhibit influenza A virus infection. Depletion of RTP4 significantly increased influenza A virus multiplication, while NS1-deficient viruses were unaffected. Mechanistically, RTP4 interacts with NS1 in an RNA-dependent manner and sequesters it from the TRIM25-RIG-I complex, thereby restoring TRIM25-mediated RIG-I K63-linked ubiquitination and subsequent activation of IRF3. Antiviral activity of RTP4 requires the evolutionarily conserved CXXC motifs and an H149 residue in the zinc finger domain, mutations of which disrupted RTP4-NS1 interaction and abrogated the ability of RTP4 to rescue RIG-I-mediated signaling. Collectively, our findings provided insights into the mechanism by which an ISG restricts influenza A virus replication by reactivating host antiviral signaling.

Bridging spatial and temporal scales of developmental gene regulation

The development of multicellular organisms relies on the precise coordination of molecular events across multiple spatial and temporal scales. Understanding how information flows from molecular interactions to cellular processes and tissue organization during development is crucial for explaining the remarkable reproducibility of complex organisms. This review explores how chromatin-encoded information is transduced from localized transcriptional events to global gene expression patterns, highlighting the challenge of bridging these scales. We discuss recent experimental findings and theoretical frameworks, emphasizing polymer physics as a tool for describing the relationship between chromatin structure and dynamics across scales. By integrating these perspectives, we aim to clarify how gene regulation is coordinated across levels of biological organization and suggest strategies for future experimental approaches.

Genetic coupling of enhancer activity and connectivity in gene expression control

Gene enhancers often form long-range contacts with promoters, but it remains unclear if the activity of enhancers and their chromosomal contacts are mediated by the same DNA sequences and recruited factors. Here, we study the effects of expression quantitative trait loci (eQTLs) on enhancer activity and promoter contacts in primary monocytes isolated from 34 male individuals. Using eQTL-Capture Hi-C and a Bayesian approach considering both intra- and inter-individual variation, we initially detect 19 eQTLs associated with enhancer-eGene promoter contacts, most of which also associate with enhancer accessibility and activity. Capitalising on these shared effects, we devise a multi-modality Bayesian strategy, identifying 629 "trimodal QTLs" jointly associated with enhancer accessibility, eGene promoter contact, and gene expression. Causal mediation analysis and CRISPR interference reveal causal relationships between these three modalities. Many detected QTLs overlap disease susceptibility loci and influence the predicted binding of myeloid transcription factors, including SPI1, GABPB and STAT3. Additionally, a variant associated with PCK2 promoter contact directly disrupts a CTCF binding motif and impacts promoter insulation from downstream enhancers. Jointly, our findings suggest an inherent genetic coupling of enhancer activity and connectivity in gene expression control relevant to human disease and highlight the regulatory role of genetically determined chromatin boundaries.

The α-globin super-enhancer acts in an orientation-dependent manner

Individual enhancers are defined as short genomic regulatory elements, bound by transcription factors, and able to activate cell-specific gene expression at a distance, in an orientation-independent manner. Within mammalian genomes, enhancer-like elements may be found individually or within clusters referred to as locus control regions or super-enhancers (SEs). While these behave similarly to individual enhancers with respect to cell specificity, distribution and distance, their orientation-dependence has not been formally tested. Here, using the α-globin locus as a model, we show that while an individual enhancer works in an orientation-independent manner, the direction of activity of a SE changes with its orientation. When the SE is inverted within its normal chromosomal context, expression of its normal targets, the α-globin genes, is severely reduced and the normally silent genes lying upstream of the α-globin locus are upregulated. These findings add to our understanding of enhancer-promoter specificity that precisely activate transcription.

Pervasive RNA-binding protein enrichment on TAD boundaries regulates TAD organization

Mammalian genome is hierarchically organized by CTCF and cohesin through loop extrusion mechanism to facilitate the organization of topologically associating domains (TADs). Mounting evidence suggests additional factors/mechanisms exist to orchestrate TAD formation and maintenance. In this study, we investigate the potential role of RNA-binding proteins (RBPs) in TAD organization. By integrated analyses of global RBP binding and 3D genome mapping profiles from both K562 and HepG2 cells, our study unveils the prevalent enrichment of RBPs on TAD boundaries and define boundary-associated RBPs (baRBPs). We found that baRBP binding is correlated with enhanced TAD insulation strength and in a CTCF-independent manner. Moreover, baRBP binding is associated with nascent promoter transcription. Additional experimental testing was performed using RBFox2 as a paradigm. Knockdown of RBFox2 in K562 cells causes mild TAD reorganization. Moreover, RBFox2 enrichment on TAD boundaries is a conserved phenomenon in C2C12 myoblast (MB) cells. RBFox2 is downregulated and its bound boundaries are remodeled during MB differentiation into myotubes. Finally, transcriptional inhibition indeed decreases RBFox2 binding and disrupts TAD boundary insulation. Altogether, our findings demonstrate that RBPs can play an active role in modulating TAD organization through co-transcriptional association and synergistic actions with nascent promoter transcripts.

An interferon-stimulated long non-coding RNA USP30-AS1 as an immune modulator in influenza A virus infection

Long non-coding RNAs (lncRNAs) are essential components of innate immunity, maintaining the functionality of immune systems that control virus infection. However, how lncRNAs engage immune responses during influenza A virus (IAV) infection remains unclear. Here, we show that lncRNA USP30-AS1 is up-regulated by infection of multiple different IAV subtypes and is required for tuning inflammatory and antiviral response in IAV infection. Genetically inactivation of USP30-AS1 enhances viral protein synthesis and viral growth. USP30-AS1 is an interferon-stimulated gene, and the induction of USP30-AS1 can be achieved by JAK-STAT mediated signaling activation. The immune regulation of USP30-AS1 is independent of its proximal protein-coding gene USP30. In IAV infection, deletion of USP30-AS1 unleashes high systemic inflammatory responses involving a broad range of pro-inflammatory factors, suggesting USP30-AS1 as a critical modulator of immune responses in IAV infection. Furthermore, we established a database providing well-annotated host gene expression profiles IAV infection or immune stimulation.

Stiff extracellular matrix drives the differentiation of mesenchymal stem cells toward osteogenesis by the multiscale 3D genome reorganization

Extracellular matrix (ECM) stiffness is a major driver of stem cell fate. However, the involvement of the three-dimensional (3D) genomic reorganization in response to ECM stiffness remains unclear. Here, we generated comprehensive 3D chromatin landscapes of mesenchymal stem cells (MSCs) exposed to various ECM stiffness. We found that there were more long-range chromatin interactions, but less compartment A in MSCs cultured on stiff ECM than those cultured on soft ECM. However, the switch from compartment B in MSCs cultured on soft ECM to compartment A in MSCs cultured on stiff ECM included genes encoding proteins primarily enriched in cytoskeleton organization. At the topologically associating domains (TADs) level, stiff ECM tends to have merged TADs on soft ECM. These merged TADs on stiff ECM include upregulated genes encoding proteins enriched in osteogenesis, such as SP1, ETS1, and DCHS1, which were validated by quantitative real-time polymerase chain reaction and found to be consistent with the increase of alkaline phosphatase staining. Knockdown of SP1 or ETS1 led to the downregulation of osteogenic marker genes, including COL1A1, RUNX2, ALP, and OCN in MSCs cultured on stiff ECM. Our study provides an important insight into the stiff ECM-mediated promotion of MSC differentiation towards osteogenesis, emphasizing the influence of mechanical cues on the reorganization of 3D genome architecture and stem cell fate.

The three-dimensional genome drives the evolution of asymmetric gene duplicates via enhancer capture-divergence

Previous evolutionary models of duplicate gene evolution have overlooked the pivotal role of genome architecture. Here, we show that proximity-based regulatory recruitment by distally duplicated genes is an efficient mechanism for modulating tissue-specific production of preexisting proteins. By leveraging genomic asymmetries, we performed a coexpression analysis on Drosophila melanogaster tissue data to show the generality of enhancer capture-divergence (ECD) as a significant evolutionary driver of asymmetric, distally duplicated genes. We use the recently evolved gene HP6/Umbrea as an example of the ECD process. By assaying genome-wide chromosomal conformations in multiple Drosophila species, we show that HP6/Umbrea was inserted near a preexisting, long-distance three-dimensional genomic interaction. We then use this data to identify a newly found enhancer (FLEE1), buried within the coding region of the highly conserved, essential gene MFS18, that likely neofunctionalized HP6/Umbrea. Last, we demonstrate ancestral transcriptional coregulation of HP6/Umbrea's future insertion site, illustrating how enhancer capture provides a highly evolvable, one-step solution to Ohno's dilemma.

Node features of chromosome structure networks and their connections to genome annotation

The 3D conformations of chromosomes can encode biological significance, and the implications of such structures have been increasingly appreciated recently. Certain chromosome structural features, such as A/B compartmentalization, are frequently extracted from Hi-C pairwise genome contact information (physical association between different regions of the genome) and compared with linear annotations of the genome, such as histone modifications and lamina association. We investigate how additional properties of chromosome structure can be deduced using an abstract graph representation of the contact heatmap, and describe specific network properties that can have a strong connection with some of these biological annotations. We constructed chromosome structure networks (CSNs) from bulk Hi-C data and calculated a set of site-resolved (node-based) network properties. These properties are useful for characterizing certain aspects of chromosomal structure. We examined the ability of network properties to differentiate several scenarios, such as haploid vs diploid cells, partially inverted nuclei vs conventional architecture, depletion of chromosome architectural proteins, and structural changes during cell development. We also examined the connection between network properties and a series of other linear annotations, such as histone modifications and chromatin states including poised promoter and enhancer labels. We found that semi-local network properties exhibit greater capability in characterizing genome annotations compared to diffusive or ultra-local node features. For example, the local square clustering coefficient can be a strong classifier of lamina-associated domains. We demonstrated that network properties can be useful for highlighting large-scale chromosome structure differences that emerge in different biological situations.

The first chicken oocyte nucleus whole transcriptomic profile defines the spectrum of maternal mRNA and non-coding RNA genes transcribed by the lampbrush chromosomes

Lampbrush chromosomes, with their unusually high rate of nascent RNA synthesis, provide a valuable model for studying mechanisms of global transcriptome up-regulation. Here, we obtained a whole-genomic profile of transcription along the entire length of all lampbrush chromosomes in the chicken karyotype. With nuclear RNA-seq, we obtained information about a wider set of transcripts, including long non-coding RNAs retained in the nucleus and stable intronic sequence RNAs. For a number of protein-coding genes, we visualized their nascent transcripts on the lateral loops of lampbrush chromosomes by RNA-FISH. The set of genes transcribed on the lampbrush chromosomes is required for basic cellular processes and is characterized by a broad expression pattern. We also present the first high-throughput transcriptome characterization of miRNAs and piRNAs in chicken oocytes at the lampbrush chromosome stage. Major targets of predicted piRNAs include CR1 and long terminal repeat (LTR) containing retrotransposable elements. Transcription of tandem repeat arrays was demonstrated by alignment against the whole telomere-to-telomere chromosome assemblies. We show that transcription of telomere-derived RNAs is initiated at adjacent LTR elements. We conclude that hypertranscription on the lateral loops of giant lampbrush chromosomes is required for synthesizing large amounts of transferred to the embryo maternal RNA for thousands of genes.

Epstein-Barr virus induces host shutoff extensively via BGLF5-independent mechanisms

Epstein-Barr virus (EBV) is a ubiquitous oncogenic virus associated with multiple cancers and autoimmune diseases. Unlike most herpesviruses, EBV reactivation from latency occurs asymptomatically, allowing it to spread efficiently to other hosts. However, available models are limited by the inefficient and asynchronous reactivation from latency into lytic replication. To address this problem, we develop a dual-fluorescent lytic reporter (DFLR) EBV that specifically labels cells in the early and late stages of replication. Using lymphoblastoid cell lines transformed by DFLR EBV as a model for EBV reactivation in B cells, we observe extensive reprogramming of the host cell transcriptome during lytic cycle progression. This includes widespread shutoff of host gene expression and disruption of mRNA processing. Unexpectedly, host shutoff remains extensive even in cells infected with DFLR EBV deleted for the BGLF5 nuclease. These findings implicate BGLF5-independent mechanisms as the primary drivers of host transcriptome remodeling during EBV lytic replication.

A comparative roadmap of PIWI-interacting RNAs across seven species reveals insights into de novo piRNA-precursor formation in mammals

PIWI-interacting RNAs (piRNAs) play a crucial role in safeguarding genome integrity by silencing mobile genetic elements. From flies to humans, piRNAs originate from long single-stranded precursors encoded by genomic piRNA clusters. How piRNA clusters form to adapt to genomic invaders and evolve to maintain protection remain key outstanding questions. Here, we generate a roadmap of piRNA clusters across seven species that highlights both similarities and variations. In mammals, we identify transcriptional readthrough as a mechanism to generate piRNAs from transposon insertions (piCs) downstream of genes (DoG). Together with the well-known stress-dependent DoG transcripts, our findings suggest a molecular mechanism for the formation of piRNA clusters in response to retroviral invasion. Finally, we identify a class of dynamic piRNA clusters in humans, underscoring unique features of human germ cell biology. Our results advance the understanding of conserved principles and species-specific variations in piRNA biology and provide tools for future studies.

Distinct negative-sense RNA viruses induce a common set of transcripts encoding proteins forming an extensive network

The large group of negative-strand RNA viruses (NSVs) comprises many important pathogens. To identify conserved patterns in host responses, we systematically compared changes in the cellular RNA levels after infection of human hepatoma cells with nine different NSVs of different virulence degrees. RNA sequencing experiments indicated that the amount of viral RNA in host cells correlates with the number of differentially expressed host cell transcripts. Time-resolved differential gene expression analysis revealed a common set of 178 RNAs that are regulated by all NSVs analyzed. A newly developed open access web application allows downloads and visualizations of all gene expression comparisons for individual viruses over time or between several viruses. Most of the genes included in the core set of commonly differentially expressed genes (DEGs) encode proteins that serve as membrane receptors, signaling proteins and regulators of transcription. They mainly function in signal transduction and control immunity, metabolism, and cell survival. One hundred sixty-five of the DEGs encode host proteins from which 47 have already been linked to the regulation of viral infections in previous studies and 89 proteins form a complex interaction network that may function as a core hub to control NSV infections.IMPORTANCEThe infection of cells with negative-strand RNA viruses leads to the differential expression of many host cell RNAs. The differential spectrum of virus-regulated RNAs reflects a large variety of events including anti-viral responses, cell remodeling, and cell damage. Here, these virus-specific differences and similarities in the regulated RNAs were measured in a highly standardized model. A newly developed app allows interested scientists a wide range of comparisons and visualizations.

Genome structural dynamics: insights from Gaussian network analysis of Hi-C data

Characterization of the spatiotemporal properties of the chromatin is essential to gaining insights into the physical bases of gene co-expression, transcriptional regulation and epigenetic modifications. The Gaussian network model (GNM) has proven in recent work to serve as a useful tool for modeling chromatin structural dynamics, using as input high-throughput chromosome conformation capture data. We focus here on the exploration of the collective dynamics of chromosomal structures at hierarchical levels of resolution, from single gene loci to topologically associating domains or entire chromosomes. The GNM permits us to identify long-range interactions between gene loci, shedding light on the role of cross-correlations between distal regions of the chromosomes in regulating gene expression. Notably, GNM analysis performed across diverse cell lines highlights the conservation of the global/cooperative movements of the chromatin across different types of cells. Variations driven by localized couplings between genomic loci, on the other hand, underlie cell differentiation, underscoring the significance of the four-dimensional properties of the genome in defining cellular identity. Finally, we demonstrate the close relation between the cell type-dependent mobility profiles of gene loci and their gene expression patterns, providing a clear demonstration of the role of chromosomal 4D features in defining cell-specific differential expression of genes.

Trained immunity is regulated by T cell-induced CD40-TRAF6 signaling

Trained immunity is characterized by histone modifications and metabolic changes in innate immune cells following exposure to inflammatory signals, leading to heightened responsiveness to secondary stimuli. Although our understanding of the molecular regulation of trained immunity has increased, the role of adaptive immune cells herein remains largely unknown. Here, we show that T cells modulate trained immunity via cluster of differentiation 40-tissue necrosis factor receptor-associated factor 6 (CD40-TRAF6) signaling. CD40-TRAF6 inhibition modulates functional, transcriptomic, and metabolic reprogramming and modifies histone 3 lysine 4 trimethylation associated with trained immunity. Besides in vitro studies, we reveal that single-nucleotide polymorphisms in the proximity of CD40 are linked to trained immunity responses in vivo and that combining CD40-TRAF6 inhibition with cytotoxic T lymphocyte antigen 4-immunoglobulin (CTLA4-Ig)-mediated co-stimulatory blockade induces long-term graft acceptance in a murine heart transplantation model. Combined, our results reveal that trained immunity is modulated by CD40-TRAF6 signaling between myeloid and adaptive immune cells and that this can be leveraged for therapeutic purposes.

STAG2 mutations reshape the cohesin-structured spatial chromatin architecture to drive gene regulation in acute myeloid leukemia

Cohesin shapes the chromatin architecture, including enhancer-promoter interactions. Its components, especially STAG2, but not its paralog STAG1, are frequently mutated in myeloid malignancies. To elucidate the underlying mechanisms of leukemogenesis, we comprehensively characterized genetic, transcriptional, and chromatin conformational changes in acute myeloid leukemia (AML) patient samples. Specific loci displayed altered cohesin occupancy, gene expression, and local chromatin activation, which were not compensated by the remaining STAG1-cohesin. These changes could be linked to disrupted spatial chromatin looping in cohesin-mutated AMLs. Complementary depletion of STAG2 or STAG1 in primary human hematopoietic progenitors (HSPCs) revealed effects resembling STAG2-mutant AML-specific changes following STAG2 knockdown, not invoked by the depletion of STAG1. STAG2-deficient HSPCs displayed impaired differentiation capacity and maintained HSPC-like gene expression. This work establishes STAG2 as a key regulator of chromatin contacts, gene expression, and differentiation in the hematopoietic system and identifies candidate target genes that may be implicated in human leukemogenesis.

Artificially inserted strong promoter containing multiple G-quadruplexes induces long-range chromatin modification

Although the role of G-quadruplex (G4) DNA structures has been suggested in chromosomal looping this was not tested directly. Here, to test causal function, an array of G4s, or control sequence that does not form G4s, were inserted within chromatin in cells. In vivo G4 formation of the inserted G4 sequence array, and not the control sequence, was confirmed using G4-selective antibody. Compared to the control insert, we observed a remarkable increase in the number of 3D chromatin looping interactions from the inserted G4 array. This was evident within the immediate topologically associated domain (TAD) and throughout the genome. Locally, recruitment of enhancer histone marks and the transcriptional coactivator p300/Acetylated-p300 increased in the G4-array, but not in the control insertion. Resulting promoter-enhancer interactions and gene activation were clear up to 5 Mb away from the insertion site. Together, these show the causal role of G4s in enhancer function and long-range chromatin interactions. Mechanisms of 3D topology are primarily based on DNA-bound architectural proteins that induce/stabilize long-range interactions. Involvement of the underlying intrinsic DNA sequence/structure in 3D looping shown here therefore throws new light on how long-range chromosomal interactions might be induced or maintained.

Multi-omic profiling of pathogen-stimulated primary immune cells

We performed long-read transcriptome and proteome profiling of pathogen-stimulated peripheral blood mononuclear cells (PBMCs) from healthy donors to discover new transcript and protein isoforms expressed during immune responses to diverse pathogens. Long-read transcriptome profiling reveals novel sequences and isoform switching induced upon pathogen stimulation, including transcripts that are difficult to detect using traditional short-read sequencing. Widespread loss of intron retention occurs as a common result of all pathogen stimulations. We highlight novel transcripts of NFKB1 and CASP1 that may indicate novel immunological mechanisms. RNA expression differences did not result in differences in the amounts of secreted proteins. Clustering analysis of secreted proteins revealed a correlation between chemokine (receptor) expression on the RNA and protein levels in C. albicans- and poly(I:C)-stimulated PBMCs. Isoform aware long-read sequencing of pathogen-stimulated immune cells highlights the potential of these methods to identify novel transcripts, revealing a more complex transcriptome landscape than previously appreciated.

Cell-type-specific loops linked to RNA polymerase II elongation in human neural differentiation

DNA is folded into higher-order structures that shape and are shaped by genome function. The role of long-range loops in the establishment of new gene expression patterns during cell fate transitions remains poorly understood. Here, we investigate the link between cell-specific loops and RNA polymerase II (RNA Pol II) during neural lineage commitment. We find thousands of loops decommissioned or gained de novo upon differentiation of human induced pluripotent stem cells (hiPSCs) to neural progenitor cells (NPCs) and post-mitotic neurons. During hiPSC-to-NPC and NPC-to-neuron transitions, genes changing from RNA Pol II initiation to elongation are >4-fold more likely to anchor cell-specific loops than repressed genes. Elongated genes exhibit significant mRNA upregulation when connected in cell-specific promoter-enhancer loops but not invariant promoter-enhancer loops or promoter-promoter loops or when unlooped. Genes transitioning from repression to RNA Pol II initiation exhibit a slight mRNA increase independent of loop status. Our data link cell-specific loops and robust RNA Pol II-mediated elongation during neural cell fate transitions.

RUNX1::ETO and CBFβ::MYH11 converge on aberrant activation of BCAT1 to confer a therapeutic vulnerability in core-binding factor-acute myeloid leukaemia

Effectively targeting transcription factors in therapeutic interventions remains challenging, especially in core-binding factor-acute myeloid leukaemia (CBF-AML) characterized by RUNX1::ETO and CBFβ::MYH11 fusions. However, recent studies have drawn attention towards aberrant amino acid metabolisms as actionable therapeutic targets. Here, by integrating the expression profile and genetic makeup in AML cohort, we found higher BCAT1 expression in CBF-AML patients compared with other subtypes. Metabolic profiling revealed that high BCAT1 expression led to reprogrammed branch amino acid metabolism in CBF-AML and was associated with sphingolipid pathway relating to the fitness of leukaemia cells, supported by transcriptomic profiling. Mechanistically, we demonstrated in cell lines and primary patient samples that BCAT1 was directly activated by RUNX1::ETO and CBFβ::MYH11 fusion proteins similarly in a RUNX1-dependent manner through rewiring chromatin conformation at the BCAT1 gene locus. Furthermore, BCAT1 inhibition resulted in blunted cell cycle, enhanced apoptosis and myeloid differentiation of CBF-AML cells in vitro, and alleviated leukaemia burden and prolonged survival in vivo. Importantly, pharmacological inhibition of BCAT1 using the specific inhibitor Gabapentin demonstrated therapeutic effects, as evidenced by delayed leukaemia progression and improved survival in vivo. In conclusion, our study uncovers BCAT1 as a genetic vulnerability and a promising targeted therapeutic opportunity for CBF-AML.

Dynamic Changes in Histone Modifications Are Associated with Differential Chromatin Interactions

Eukaryotic genomes are organized into chromatin domains through long-range chromatin interactions which are mediated by the binding of architectural proteins, such as CTCF and cohesin, and histone modifications. Based on the published Hi-C and ChIP-seq datasets in human monocyte-derived macrophages, we identified 206 and 127 differential chromatin interactions (DCIs) that were not located within transcription readthrough regions in influenza A virus- and interferon β-treated cells, respectively, and found that the binding positions of CTCF and RAD21 within more than half of the DCI sites did not change. However, five histone modifications, H3K4me3, H3K27ac, H3K36me3, H3K9me3, and H3K27me3, showed significantly more dramatic changes than CTCF and RAD21 within the DCI sites. For H3K4me3, H3K27ac, H3K36me3, and H3K27me3, significantly more dramatic changes were observed outside than within the DCI sites. We further applied a motif scanning approach to discover proteins that might correlate with changes in histone modifications and chromatin interactions and found that PRDM9, ZNF384, and STAT2 frequently bound to DNA sequences corresponding to 1 kb genomic intervals with gains or losses of a histone modification within the DCI sites. This study explores the dynamic regulation of chromatin interactions and extends the current knowledge of the relationship between histone modifications and chromatin interactions.

IRF2BP2 counteracts the ATF7/JDP2 AP-1 heterodimer to prevent inflammatory overactivation in acute myeloid leukemia (AML) cells

Acute myeloid leukemia (AML) is a hematological malignancy characterized by abnormal proliferation and accumulation of immature myeloid cells in the bone marrow. Inflammation plays a crucial role in AML progression, but excessive activation of cell-intrinsic inflammatory pathways can also trigger cell death. IRF2BP2 is a chromatin regulator implicated in AML pathogenesis, although its precise role in this disease is not fully understood. In this study, we demonstrate that IRF2BP2 interacts with the AP-1 heterodimer ATF7/JDP2, which is involved in activating inflammatory pathways in AML cells. We show that IRF2BP2 is recruited by the ATF7/JDP2 dimer to chromatin and counteracts its gene-activating function. Loss of IRF2BP2 leads to overactivation of inflammatory pathways, resulting in strongly reduced proliferation. Our research indicates that a precise equilibrium between activating and repressive transcriptional mechanisms creates a pro-oncogenic inflammatory environment in AML cells. The ATF7/JDP2-IRF2BP2 regulatory axis is likely a key regulator of this process and may, therefore, represent a promising therapeutic vulnerability for AML. Thus, our study provides new insights into the molecular mechanisms underlying AML pathogenesis and identifies a potential therapeutic target for AML treatment.

Downstream-of-gene (DoG) transcripts contribute to an imbalance in the cancer cell transcriptome

Downstream-of-gene (DoG) transcripts are an emerging class of noncoding RNAs. However, it remains largely unknown how DoG RNA production is regulated and whether alterations in DoG RNA signatures exist in major cancers. Here, through transcriptomic analyses of matched tumors and nonneoplastic tissues and cancer cell lines, we reveal a comprehensive catalog of DoG RNA signatures. Through separate lines of evidence, we support the biological importance of DoG RNAs in carcinogenesis. First, we show tissue-specific and stage-specific differential expression of DoG RNAs in tumors versus paired normal tissues with their respective host genes involved in tumor-promoting versus tumor-suppressor pathways. Second, we identify that differential DoG RNA expression is associated with poor patient survival. Third, we identify that DoG RNA induction is a consequence of treating colon cancer cells with the topoisomerase I (TOP1) poison camptothecin and following TOP1 depletion. Our results underlie the significance of DoG RNAs and TOP1-dependent regulation of DoG RNAs in diversifying and modulating the cancer transcriptome.

A continuum of zinc finger transcription factor retention on native chromatin underlies dynamic genome organization

Transcription factor (TF) residence on chromatin translates into quantitative transcriptional or structural outcomes on genome. Commonly used formaldehyde crosslinking fixes TF-DNA interactions cumulatively and compromises the measured occupancy level. Here we mapped the occupancy level of global or individual zinc finger TFs like CTCF and MAZ, in the form of highly resolved footprints, on native chromatin. By incorporating reinforcing perturbation conditions, we established S-score, a quantitative metric to proxy the continuum of CTCF or MAZ retention across different motifs on native chromatin. The native chromatin-retained CTCF sites harbor sequence features within CTCF motifs better explained by S-score than the metrics obtained from other crosslinking or native assays. CTCF retention on native chromatin correlates with local SUMOylation level, and anti-correlates with transcriptional activity. The S-score successfully delineates the otherwise-masked differential stability of chromatin structures mediated by CTCF, or by MAZ independent of CTCF. Overall, our study established a paradigm continuum of TF retention across binding sites on native chromatin, explaining the dynamic genome organization.

Transcription regulates the spatio-temporal dynamics of genes through micro-compartmentalization

Although our understanding of the involvement of heterochromatin architectural factors in shaping nuclear organization is improving, there is still ongoing debate regarding the role of active genes in this process. In this study, we utilize publicly-available Micro-C data from mouse embryonic stem cells to investigate the relationship between gene transcription and 3D gene folding. Our analysis uncovers a nonmonotonic - globally positive - correlation between intragenic contact density and Pol II occupancy, independent of cohesin-based loop extrusion. Through the development of a biophysical model integrating the role of transcription dynamics within a polymer model of chromosome organization, we demonstrate that Pol II-mediated attractive interactions with limited valency between transcribed regions yield quantitative predictions consistent with chromosome-conformation-capture and live-imaging experiments. Our work provides compelling evidence that transcriptional activity shapes the 4D genome through Pol II-mediated micro-compartmentalization.

Extensive long-range polycomb interactions and weak compartmentalization are hallmarks of human neuronal 3D genome

Chromatin architecture regulates gene expression and shapes cellular identity, particularly in neuronal cells. Specifically, polycomb group (PcG) proteins enable establishment and maintenance of neuronal cell type by reorganizing chromatin into repressive domains that limit the expression of fate-determining genes and sustain distinct gene expression patterns in neurons. Here, we map the 3D genome architecture in neuronal and non-neuronal cells isolated from the Wernicke's area of four human brains and comprehensively analyze neuron-specific aspects of chromatin organization. We find that genome segregation into active and inactive compartments is greatly reduced in neurons compared to other brain cells. Furthermore, neuronal Hi-C maps reveal strong long-range interactions, forming a specific network of PcG-mediated contacts in neurons that is nearly absent in other brain cells. These interacting loci contain developmental transcription factors with repressed expression in neurons and other mature brain cells. But only in neurons, they are rich in bivalent promoters occupied by H3K4me3 histone modification together with H3K27me3, which points to a possible functional role of PcG contacts in neurons. Importantly, other layers of chromatin organization also exhibit a distinct structure in neurons, characterized by an increase in short-range interactions and a decrease in long-range ones.

Single-molecule imaging reveals a direct role of CTCF's zinc fingers in SA interaction and cluster-dependent RNA recruitment

CTCF is a zinc finger protein associated with transcription regulation that also acts as a barrier factor for topologically associated domains (TADs) generated by cohesin via loop extrusion. These processes require different properties of CTCF-DNA interaction, and it is still unclear how CTCF's structural features may modulate its diverse roles. Here, we employ single-molecule imaging to study both full-length CTCF and truncation mutants. We show that CTCF enriches at CTCF binding sites (CBSs), displaying a longer lifetime than observed previously. We demonstrate that the zinc finger domains mediate CTCF clustering and that clustering enables RNA recruitment, possibly creating a scaffold for interaction with RNA-binding proteins like cohesin's subunit SA. We further reveal a direct recruitment and an increase of SA residence time by CTCF bound at CBSs, suggesting that CTCF-SA interactions are crucial for cohesin stability on chromatin at TAD borders. Furthermore, we establish a single-molecule T7 transcription assay and show that although a transcribing polymerase can remove CTCF from CBSs, transcription is impaired. Our study shows that context-dependent nucleic acid binding determines the multifaceted CTCF roles in genome organization and transcription regulation.

A CD36-dependent non-canonical lipid metabolism program promotes immune escape and resistance to hypomethylating agent therapy in AML

Environmental lipids are essential for fueling tumor energetics, but whether these exogenous lipids transported into cancer cells facilitate immune escape remains unclear. Here, we find that CD36, a transporter for exogenous lipids, promotes acute myeloid leukemia (AML) immune evasion. We show that, separately from its established role in lipid oxidation, CD36 on AML cells senses oxidized low-density lipoprotein (OxLDL) to prime the TLR4-LYN-MYD88-nuclear factor κB (NF-κB) pathway, and exogenous palmitate transfer via CD36 further potentiates this innate immune pathway by supporting ZDHHC6-mediated MYD88 palmitoylation. Subsequently, NF-κB drives the expression of immunosuppressive genes that inhibit anti-tumor T cell responses. Notably, high-fat-diet or hypomethylating agent decitabine treatment boosts the immunosuppressive potential of AML cells by hijacking CD36-dependent innate immune signaling, leading to a dampened therapeutic effect. This work is of translational interest because lipid restriction by US Food and Drug Administration (FDA)-approved lipid-lowering statin drugs improves the efficacy of decitabine therapy by weakening leukemic CD36-mediated immunosuppression.

Multi-omics analysis reveals the dynamic interplay between Vero host chromatin structure and function during vaccinia virus infection

The genome folds into complex configurations and structures thought to profoundly impact its function. The intricacies of this dynamic structure-function relationship are not well understood particularly in the context of viral infection. To unravel this interplay, here we provide a comprehensive investigation of simultaneous host chromatin structural (via Hi-C and ATAC-seq) and functional changes (via RNA-seq) in response to vaccinia virus infection. Over time, infection significantly impacts global and local chromatin structure by increasing long-range intra-chromosomal interactions and B compartmentalization and by decreasing chromatin accessibility and inter-chromosomal interactions. Local accessibility changes are independent of broad-scale chromatin compartment exchange (~12% of the genome), underscoring potential independent mechanisms for global and local chromatin reorganization. While infection structurally condenses the host genome, there is nearly equal bidirectional differential gene expression. Despite global weakening of intra-TAD interactions, functional changes including downregulated immunity genes are associated with alterations in local accessibility and loop domain restructuring. Therefore, chromatin accessibility and local structure profiling provide impactful predictions for host responses and may improve development of efficacious anti-viral counter measures including the optimization of vaccine design.

Attraction and disruption: how loop extrusion and compartmentalisation shape the nuclear genome

Chromatin loops, which bring two distal loci of the same chromosome into close physical proximity, are the ubiquitous units of the three-dimensional genome. Recent advances in understanding the spatial organisation of chromatin suggest that several distinct mechanisms control chromatin interactions, such as loop extrusion by cohesin complexes, compartmentalisation by phase separation, direct protein-protein interactions and others. Here, we review different types of chromatin loops and highlight the factors and processes involved in their regulation. We discuss how loop extrusion and compartmentalisation shape chromatin interactions and how these two processes can either positively or negatively influence each other.

Mechanistic drivers of chromatin organization into compartments

The human genome is not just a simple string of DNA, it is a complex and dynamic entity intricately folded within the cell's nucleus. This three-dimensional organization of chromatin, the combination of DNA and proteins in the nucleus, is crucial for many biological processes and has been prominently studied for its intricate relationship to gene expression. Indeed, the transcriptional machinery does not operate in isolation but interacts intimately with the folded chromatin structure. Techniques for chromatin conformation capture, including genome-wide sequencing approaches, have revealed key organizational features of chromatin, such as the formation of loops by CCCTC-binding factor (CTCF) and the division of loci into chromatin compartments. While much of the recent research and reviews have focused on CTCF loops, we discuss several new revelations that have emerged concerning chromatin compartments, with a particular focus on what is known about mechanistic drivers of compartmentalization. These insights challenge the traditional views of chromatin organization and reveal the complexity behind the formation and maintenance of chromatin compartments.

Long-read sequencing for 29 immune cell subsets reveals disease-linked isoforms

Alternative splicing events are a major causal mechanism for complex traits, but they have been understudied due to the limitation of short-read sequencing. Here, we generate a full-length isoform annotation of human immune cells from an individual by long-read sequencing for 29 cell subsets. This contains a number of unannotated transcripts and isoforms such as a read-through transcript of TOMM40-APOE in the Alzheimer's disease locus. We profile characteristics of isoforms and show that repetitive elements significantly explain the diversity of unannotated isoforms, providing insight into the human genome evolution. In addition, some of the isoforms are expressed in a cell-type specific manner, whose alternative 3'-UTRs usage contributes to their specificity. Further, we identify disease-associated isoforms by isoform switch analysis and by integration of several quantitative trait loci analyses with genome-wide association study data. Our findings will promote the elucidation of the mechanism of complex diseases via alternative splicing.

TAD border deletion at the Kit locus causes tissue-specific ectopic activation of a neighboring gene

Topologically associated domains (TADs) restrict promoter-enhancer interactions, thereby maintaining the spatiotemporal pattern of gene activity. However, rearrangements of the TADs boundaries do not always lead to significant changes in the activity pattern. Here, we investigated the consequences of the TAD boundaries deletion on the expression of developmentally important genes encoding tyrosine kinase receptors: Kit, Kdr, Pdgfra. We used genome editing in mice to delete the TADs boundaries at the Kit locus and characterized chromatin folding and gene expression in pure cultures of fibroblasts, mast cells, and melanocytes. We found that although Kit is highly active in both mast cells and melanocytes, deletion of the TAD boundary between the Kit and Kdr genes results in ectopic activation only in melanocytes. Thus, the epigenetic landscape, namely the mutual arrangement of enhancers and actively transcribing genes, is important for predicting the consequences of the TAD boundaries removal. We also found that mice without a TAD border between the Kit and Kdr genes have a phenotypic manifestation of the mutation - a lighter coloration. Thus, the data obtained shed light on the principles of interaction between the 3D chromatin organization and epigenetic marks in the regulation of gene activity.

Suv39h-catalyzed H3K9me3 is critical for euchromatic genome organization and the maintenance of gene transcription

H3K9me3-dependent heterochromatin is critical for the silencing of repeat-rich pericentromeric regions and also has key roles in repressing lineage-inappropriate protein-coding genes in differentiation and development. Here, we investigate the molecular consequences of heterochromatin loss in cells deficient in both SUV39H1 and SUV39H2 (Suv39DKO), the major mammalian histone methyltransferase enzymes that catalyze heterochromatic H3K9me3 deposition. We reveal a paradoxical repression of protein-coding genes in Suv39DKO cells, with these differentially expressed genes principally in euchromatic (Tn5-accessible, H3K4me3- and H3K27ac-marked) rather than heterochromatic (H3K9me3-marked) or polycomb (H3K27me3-marked) regions. Examination of the three-dimensional (3D) nucleome reveals that transcriptomic dysregulation occurs in euchromatic regions close to the nuclear periphery in 3D space. Moreover, this transcriptomic dysregulation is highly correlated with altered 3D genome organization in Suv39DKO cells. Together, our results suggest that the nuclear lamina-tethering of Suv39-dependent H3K9me3 domains provides an essential scaffold to support euchromatic genome organization and the maintenance of gene transcription for healthy cellular function.

Multiscale modelling of chromatin 4D organization in SARS-CoV-2 infected cells

SARS-CoV-2 can re-structure chromatin organization and alter the epigenomic landscape of the host genome, but the mechanisms that produce such changes remain unclear. Here, we use polymer physics to investigate how the chromatin of the host genome is re-organized upon infection with SARS-CoV-2. We show that re-structuring of A/B compartments can be explained by a re-modulation of intra-compartment homo-typic affinities, which leads to the weakening of A-A interactions and the enhancement of A-B mixing. At the TAD level, re-arrangements are physically described by a reduction in the loop extrusion activity coupled with an alteration of chromatin phase-separation properties, resulting in more intermingling between different TADs and a spread in space of the TADs themselves. In addition, the architecture of loci relevant to the antiviral interferon response, such as DDX58 or IFIT, becomes more variable within the 3D single-molecule population of the infected model, suggesting that viral infection leads to a loss of chromatin structural specificity. Analysing the time trajectories of pairwise gene-enhancer and higher-order contacts reveals that this variability derives from increased fluctuations in the chromatin dynamics of infected cells. This suggests that SARS-CoV-2 alters gene regulation by impacting the stability of the contact network in time.

LATS1 controls CTCF chromatin occupancy and hormonal response of 3D-grown breast cancer cells

The cancer epigenome has been studied in cells cultured in two-dimensional (2D) monolayers, but recent studies highlight the impact of the extracellular matrix and the three-dimensional (3D) environment on multiple cellular functions. Here, we report the physical, biochemical, and genomic differences between T47D breast cancer cells cultured in 2D and as 3D spheroids. Cells within 3D spheroids exhibit a rounder nucleus with less accessible, more compacted chromatin, as well as altered expression of ~2000 genes, the majority of which become repressed. Hi-C analysis reveals that cells in 3D are enriched for regions belonging to the B compartment, have decreased chromatin-bound CTCF and increased fusion of topologically associating domains (TADs). Upregulation of the Hippo pathway in 3D spheroids results in the activation of the LATS1 kinase, which promotes phosphorylation and displacement of CTCF from DNA, thereby likely causing the observed TAD fusions. 3D cells show higher chromatin binding of progesterone receptor (PR), leading to an increase in the number of hormone-regulated genes. This effect is in part mediated by LATS1 activation, which favors cytoplasmic retention of YAP and CTCF removal.

Transcription regulation by long non-coding RNAs: mechanisms and disease relevance

Long non-coding RNAs (lncRNAs) outnumber protein-coding transcripts, but their functions remain largely unknown. In this Review, we discuss the emerging roles of lncRNAs in the control of gene transcription. Some of the best characterized lncRNAs have essential transcription cis-regulatory functions that cannot be easily accomplished by DNA-interacting transcription factors, such as XIST, which controls X-chromosome inactivation, or imprinted lncRNAs that direct allele-specific repression. A growing number of lncRNA transcription units, including CHASERR, PVT1 and HASTER (also known as HNF1A-AS1) act as transcription-stabilizing elements that fine-tune the activity of dosage-sensitive genes that encode transcription factors. Genetic experiments have shown that defects in such transcription stabilizers often cause severe phenotypes. Other lncRNAs, such as lincRNA-p21 (also known as Trp53cor1) and Maenli (Gm29348) contribute to local activation of gene transcription, whereas distinct lncRNAs influence gene transcription in trans. We discuss findings of lncRNAs that elicit a function through either activation of their transcription, transcript elongation and processing or the lncRNA molecule itself. We also discuss emerging evidence of lncRNA involvement in human diseases, and their potential as therapeutic targets.

Advances in the multimodal analysis of the 3D chromatin structure and gene regulation

Recent studies have demonstrated that the three-dimensional conformation of the chromatin plays a crucial role in gene regulation, with aberrations potentially leading to various diseases. Advanced methodologies have revealed a link between the chromatin conformation and biological function. This review divides these methodologies into sequencing-based and imaging-based methodologies, tracing their development over time. We particularly highlight innovative techniques that facilitate the simultaneous mapping of RNAs, histone modifications, and proteins within the context of the 3D architecture of chromatin. This multimodal integration substantially improves our ability to establish a robust connection between the spatial arrangement of molecular components in the nucleus and their functional roles. Achieving a comprehensive understanding of gene regulation requires capturing diverse data modalities within individual cells, enabling the direct inference of functional relationships between these components. In this context, imaging-based technologies have emerged as an especially promising approach for gathering spatial information across multiple components in the same cell.

The influencers' era: how the environment shapes chromatin in 3D

Environment-epigenome interactions are emerging as contributors to disease risk and health outcomes. In fact, organisms outside of the laboratory are constantly exposed to environmental changes that can influence chromatin regulation at multiple levels, potentially impacting on genome function. In this review, we will summarize recent findings on how major external cues impact on 3D chromatin organization in different experimental systems. We will describe environment-induced 3D genome alterations ranging from chromatin accessibility to the spatial distribution of the genome and discuss their role in regulating gene expression.