Multiplexed chromosome conformation capture sequencing for rapid genome-scale high-resolution detection of long-range chromatin interactions

A non-syndromic orofacial cleft risk locus links tRNA splicing defects to neural crest cell pathologies

Orofacial clefts are the most common form of congenital craniofacial malformation worldwide. The etiology of these birth defects is multifactorial, involving genetic and environmental factors. However, in most cases, the underlying causes remain unexplained, precluding a molecular understanding of disease mechanisms. Here, we integrated genome-wide association data, targeted resequencing of case and control cohorts, tissue- and cell-type-specific epigenomic profiling, and genome architecture analyses to molecularly dissect a genomic locus associated with an increased risk of non-syndromic orofacial cleft. We found that common and rare risk variants associated with orofacial cleft intersect with an enhancer (e2p24.2) that is active in human embryonic craniofacial tissue. We mapped e2p24.2 long-range interactions to a topologically associated domain harboring MYCN, DDX1, and CYRIA. We found that MYCN and DDX1, but not CYRIA, are required during craniofacial development in chicken embryos. We investigated the role of DDX1, a key component of the tRNA splicing complex, in cranial neural crest cells (cNCCs). The loss of DDX1 in cNCCs resulted in the accumulation of unspliced tRNA fragments, depletion of mature intron-containing tRNAs, and ribosome stalling at codons decoded by these tRNAs. This was accompanied by defects in both global protein synthesis and cNCC migration. We further showed that the induction of tRNA fragments is sufficient to disrupt craniofacial development. Together, these results uncovered a molecular mechanism in which impaired tRNA splicing affects cNCCs and craniofacial development and positioned MYCN, DDX1, and tRNA processing defects as risk factors in the pathogenesis of orofacial clefts.

An eRNA transcription checkpoint for diverse signal-dependent enhancer activation programs

The evidence that signal- and ligand-dependent pathways function by activating regulatory enhancer programs suggests that a 'checkpoint' strategy may underline activation of many diversely regulated enhancers. Here we report a molecular mechanism common to several acute signal- and ligand-dependent enhancer activation programs based on release of a shared enhancer RNA (eRNA) transcription checkpoint. It requires recruitment of a DNA-dependent protein kinase catalytic subunit (DNA-PKcs)-phosphorylated RING finger repressor (Krüppel-associated box)-associated protein 1 (KAP1) as a modulator, inhibiting its association with 7SK and E3 small ubiquitin-like modifier (SUMO) ligase activity on the CDK9 subunit of positive transcription elongation factor b (P-TEFb). This facilitates formation of an activated P-TEFb complex, licensing eRNA elongation. Overcoming this checkpoint for signal-dependent enhancer activation occurs in diverse pathways, including estrogen receptor-α, NF-κB-regulated proinflammatory stimulation, androgen receptor and neuronal depolarization. Therefore, a specific strategy required to convert a basal state enhancer P-TEFb complex to an active state to release a conserved checkpoint is apparently employed by several functionally important signal-regulated regulatory enhancers to implement the instructions of the endocrine and paracrine system.

Epigenetic editing at individual age-associated CpGs affects the genome-wide epigenetic aging landscape

Aging is reflected by genome-wide DNA methylation changes, which form the basis of epigenetic clocks, but it is largely unclear how these epigenetic modifications are regulated and whether they directly affect the aging process. In this study, we performed epigenetic editing at age-associated CpG sites to explore the consequences of interfering with epigenetic clocks. CRISPR-guided editing targeted at individual age-related CpGs evoked genome-wide bystander effects, which were highly reproducible and enriched at other age-associated regions. 4C-sequencing at age-associated sites revealed increased interactions with bystander modifications and other age-related CpGs. Subsequently, we multiplexed epigenetic editing in human T cells and mesenchymal stromal cells at five genomic regions that become either hypermethylated or hypomethylated upon aging. While targeted methylation seemed more stable at age-hypermethylated sites, both approaches induced bystander modifications at CpGs with the highest correlations with chronological age. Notably, these effects were simultaneously observed at CpGs that gain and lose methylation with age. Our results demonstrate that epigenetic editing can extensively modulate the epigenetic aging network and interfere with epigenetic clocks.

In vivo deletion of a GWAS-identified Myb distal enhancer acts on Myb expression, globin switching, and clinical erythroid parameters in β-thalassemia

Genome-wide association studies (GWAS) have identified numerous genetic variants linked to human diseases, mostly located in non-coding regions of the genome, particularly in putative enhancers. However, functional assessment of the non-coding GWAS variants has progressed at slow pace, since the functions of the vast majority of genomic enhancers have not been defined, impeding interpretation of disease-susceptibility variants. The HBS1L-MYB intergenic region harbors multiple SNPs associated with clinical erythroid parameters, including fetal hemoglobin levels, a feature impacting disease severity of beta-hemoglobinopathies such as sickle cell anemia and beta-thalassemia. HBS1L-MYB variants cluster in the vicinity of several MYB enhancers, altering MYB expression and globin switching. We and others have highlighted the conserved human MYB - 84kb enhancer, known as the - 81kb enhancer in the mouse, as likely candidate linked to these traits. We report here the generation of a Myb - 81kb enhancer knock-out mouse model, and shed light for the first time on its impact on steady state erythropoiesis and in beta-thalassemia in vivo.

TGF-β signaling controls neural crest developmental plasticity via SMAD2/3

The neural crest is a highly plastic stem cell population that represents an exception to the germ layer theory. Despite being of ectodermal origin, cranial neural crest cells can differentiate into skeletal derivatives typically formed by mesoderm. Here, we report that SMAD2/3-mediated transforming growth factor β (TGF-β) signaling enhances neural crest developmental potential in the chicken embryo. Our results show that TGF-β signaling modulates neural crest axial identity and directly controls the gene circuits that support skeletal differentiation. Cooperation between TGF-β and low levels of WNT signaling in the embryonic head activates cranial-specific cis-regulatory elements. Activation of TGF-β signaling reprogrammed trunk neural crest cells into adopting an anterior identity and led to the development of an improved protocol for the generation of human cranial neural crest cells. Our findings indicate TGF-β signaling is required for the specification of cranial neural crest cells, endowing them with the potential to give rise to the craniofacial skeleton.

ABCA1-Super Enhancer RNA Promotes Cholesterol Efflux, Reduces Macrophage-Mediated Inflammation and Atherosclerosis

We describe a previously uncharacterized ATP-binding cassette A1 super enhancer RNA (ABCA1-seRNA)-mediated cholesterol efflux. In addition, it promoted macrophage inflammatory cytokine release, and was causally correlated with coronary artery disease severity. Mechanistically, ABCA1-seRNA upregulated cholesterol efflux by interacting with mediator complex subunit 23 and recruiting retinoid X receptor-alpha and liver X receptor-alpha to promote ABCA1 transcription in a cis manner. Meanwhile, ABCA1-seRNA induced P65 ubiquitination degradation, and thereby repressed the macrophage inflammatory response. Consistently, overexpression of ABCA1-seRNA in ApoE-/- mice decreased plasma lipids, cytokines, and atherosclerotic plaques. These findings indicate ABCA1-seRNA is a critical epigenetic regulator that maintains cholesterol homeostasis and modulates inflammation, thus promising a therapeutic target for atherosclerotic cardiovascular diseases.

Promoter-centred chromatin interactions associated with EVI1 expression in EVI1+3q- myeloid leukaemia cells

EVI1 expression is associated with poor prognosis in myeloid leukaemia, which can result from Chr.3q alterations that juxtapose enhancers to induce EVI1 expression via long-range chromatin interactions. More often, however, EVI1 expression occurs unrelated to 3q alterations, and it remained unclear if, in these cases, EVI1 expression is similarly caused by aberrant enhancer activation. Here, we report that, in EVI1+3q- myeloid leukaemia cells, the EVI1 promoter interacts via long-range chromatin interactions with promoters of distally located, active genes, rather than with enhancer elements. Unlike in 3q+ cells, EVI1 expression and long-range interactions appear to not depend on CTCF/cohesin, though EVI1+3q- cells utilise an EVI1 promoter-proximal site to enhance its expression that is also involved in CTCF-mediated looping in 3q+ cells. Long-range interactions in 3q- cells connect EVI1 to promoters of multiple genes, whose transcription correlates with EVI1 in EVI1+3q- cell lines, suggesting a shared mechanism of transcriptional regulation. In line with this, CRISPR interference-induced silencing of two of these sites minimally, but consistently reduced EVI1 expression. Together, we provide novel evidence of features associated with EVI1 expression in 3q- leukaemia and consolidate the view that EVI1 in 3q- leukaemia is largely promoter-driven, potentially involving long-distance promoter clustering.

DeSUMOylation of chromatin-bound proteins limits the rapid transcriptional reprogramming induced by daunorubicin in acute myeloid leukemias

Genotoxicants have been used for decades as front-line therapies against cancer on the basis of their DNA-damaging actions. However, some of their non-DNA-damaging effects are also instrumental for killing dividing cells. We report here that the anthracycline Daunorubicin (DNR), one of the main drugs used to treat Acute Myeloid Leukemia (AML), induces rapid (3 h) and broad transcriptional changes in AML cells. The regulated genes are particularly enriched in genes controlling cell proliferation and death, as well as inflammation and immunity. These transcriptional changes are preceded by DNR-dependent deSUMOylation of chromatin proteins, in particular at active promoters and enhancers. Surprisingly, inhibition of SUMOylation with ML-792 (SUMO E1 inhibitor), dampens DNR-induced transcriptional reprogramming. Quantitative proteomics shows that the proteins deSUMOylated in response to DNR are mostly transcription factors, transcriptional co-regulators and chromatin organizers. Among them, the CCCTC-binding factor CTCF is highly enriched at SUMO-binding sites found in cis-regulatory regions. This is notably the case at the promoter of the DNR-induced NFKB2 gene. DNR leads to a reconfiguration of chromatin loops engaging CTCF- and SUMO-bound NFKB2 promoter with a distal cis-regulatory region and inhibition of SUMOylation with ML-792 prevents these changes.

Application of the 3C Method to Study the Developmental Genes in Drosophila Larvae

A transition from one developmental stage to another is accompanied by activation of developmental programs and corresponding gene ensembles. Changes in the spatial conformation of the corresponding loci are associated with this activation and can be investigated with the help of the Chromosome Conformation Capture (3C) methodology. Application of 3C to specific developmental stages is a sophisticated task. Here, we describe the use of the 3C method to study the spatial organization of developmental loci in Drosophila larvae. We critically analyzed the existing protocols and offered our own solutions and the optimized protocol to overcome limitations. To demonstrate the efficiency of our procedure, we studied the spatial organization of the developmental locus Dad in 3rd instar Drosophila larvae. Differences in locus conformation were found between embryonic cells and living wild-type larvae. We also observed the establishment of novel regulatory interactions in the presence of an adjacent transgene upon activation of its expression in larvae. Our work fills the gap in the application of the 3C method to Drosophila larvae and provides a useful guide for establishing 3C on an animal model.

LncRNA EPR-induced METTL7A1 modulates target gene translation

EPR is a long non-coding RNA (lncRNA) that controls cell proliferation in mammary gland cells by regulating gene transcription. Here, we report on Mettl7a1 as a direct target of EPR. We show that EPR induces Mettl7a1 transcription by rewiring three-dimensional chromatin interactions at the Mettl7a1 locus. Our data indicate that METTL7A1 contributes to EPR-dependent inhibition of TGF-β signaling. METTL7A1 is absent in tumorigenic murine mammary gland cells and its human ortholog (METTL7A) is downregulated in breast cancers. Importantly, re-expression of METTL7A1 in 4T1 tumorigenic cells attenuates their transformation potential, with the putative methyltransferase activity of METTL7A1 being dispensable for its biological functions. We found that METTL7A1 localizes in the cytoplasm whereby it interacts with factors implicated in the early steps of mRNA translation, associates with ribosomes, and affects the levels of target proteins without altering mRNA abundance. Overall, our data indicates that METTL7A1-a transcriptional target of EPR-modulates translation of select transcripts.

JMJD3 intrinsically disordered region links the 3D-genome structure to TGFβ-dependent transcription activation

Enhancers are key regulatory elements that govern gene expression programs in response to developmental signals. However, how multiple enhancers arrange in the 3D-space to control the activation of a specific promoter remains unclear. To address this question, we exploited our previously characterized TGFβ-response model, the neural stem cells, focusing on a ~374 kb locus where enhancers abound. Our 4C-seq experiments reveal that the TGFβ pathway drives the assembly of an enhancer-cluster and precise gene activation. We discover that the TGFβ pathway coactivator JMJD3 is essential to maintain these structures. Using live-cell imaging techniques, we demonstrate that an intrinsically disordered region contained in JMJD3 is involved in the formation of phase-separated biomolecular condensates, which are found in the enhancer-cluster. Overall, in this work we uncover novel functions for the coactivator JMJD3, and we shed light on the relationships between the 3D-conformation of the chromatin and the TGFβ-driven response during mammalian neurogenesis.

Here the authors demonstrate that TGFβ drives multi-enhancer contacts and ultimately gene activation during neuronal commitment, and that this requires the intrinsically disordered region (IDR) of the histone demethylase JMJD3 likely through its role in promoting phase-separated biomolecular condensates.

How Genomes Emerge, Function, and Evolve: Living Systems Emergence-Genotype-Phenotype-Multilism-Genome/Systems Ecology

What holds together the world in its innermost, what life is, how it emerges, functions, and evolves, has not only been an epic matter of endless romantic sunset poetry and philosophy, but also manifests explicitly in its perhaps most central organization unit-genomes. Their 3D architecture and dynamics, including the interaction networks of regulatory elements, obviously co-evolved as inseparable systems allowing the physical storage, expression, and replication of genetic information. Since we were able to fill finally the much-debated centennial gaps in their 3D architecture and dynamics, now entire new perspectives open beyond epigenetics reaching as far as a general understanding of living systems: besides the previously known DNA double helix and nucleosome structure, the latter compact into a chromatin quasi-fibre folded into stable loops forming stable multi-loop aggregates/rosettes connected by linkers, creating hence the again already known chromosome arms and entire chromosomes forming the cell nucleus. Instantly and for the first time this leads now to a consistent and cross-proven systems statistical mechanics genomics framework elucidating genome intrinsic function and regulation including various components. It balances stability/flexibility ensuring genome integrity, enabling expression/regulation of genetic information, as well as genome replication/spread. Furthermore, genotype and phenotype are multiplisticly entangled being evolutionarily the outcome of both Darwinian natural selection and Lamarckian self-referenced manipulation-all embedded in even broader genome ecology (autopoietic) i(!)n- and environmental scopes. This allows formulating new meta-level functional semantics of genomics, i.e. notions as communication of genes, genomes, and information networks, architectural and dynamic spaces for creativity and innovation, or genomes as central geno-/phenotype entanglements. Beyond and most fundamentally, the paradoxical-seeming local equilibrium substance stability in its entity though far from a universal heat-death-like equilibrium is solved, and system irreversibility, time directionality, and thus the emergence of existence are clarified. Consequently, real deep understandings of genomes, life, and complex systems in general appear in evolutionary perspectives as well as from systems analyses, via system damage/disease (its repair/cure and manipulation) as far as the understanding of extraterrestrial life, the de novo creation and thus artificial life, and even the raison d'etre.

Simulation of Different Three-Dimensional Models of Whole Interphase Nuclei Compared to Experiments - A Consistent Scale-Bridging Simulation Framework for Genome Organization

The three-dimensional architecture of chromosomes, their arrangement, and dynamics within cell nuclei are still subject of debate. Obviously, the function of genomes-the storage, replication, and transcription of genetic information-has closely coevolved with this architecture and its dynamics, and hence are closely connected. In this work a scale-bridging framework investigates how of the 30 nm chromatin fibre organizes into chromosomes including their arrangement and morphology in the simulation of whole nuclei. Therefore, mainly two different topologies were simulated with corresponding parameter variations and comparing them to experiments: The Multi-Loop-Subcompartment (MLS) model, in which (stable) small loops form (stable) rosettes, connected by chromatin linkers, and the Random-Walk/Giant-Loop (RW/GL) model, in which large loops are attached to a flexible non-protein backbone, were simulated for various loop and linker sizes. The 30 nm chromatin fibre was modelled as a polymer chain with stretching, bending and excluded volume interactions. A spherical boundary potential simulated the confinement to nuclei with different radii. Simulated annealing and Brownian Dynamics methods were applied in a four-step decondensation procedure to generate from metaphase decondensated interphase configurations at thermodynamical equilibrium. Both the MLS and the RW/GL models form chromosome territories, with different morphologies: The MLS rosettes result in distinct subchromosomal domains visible in electron and confocal laser scanning microscopic images. In contrast, the big RW/GL loops lead to a mostly homogeneous chromatin distribution. Even small changes of the model parameters induced significant rearrangements of the chromatin morphology. The low overlap of chromosomes, arms, and subchromosomal domains observed in experiments agrees only with the MLS model. The chromatin density distribution in CLSM image stacks reveals a bimodal behaviour in agreement with recent experiments. Combination of these results with a variety of (spatial distance) measurements favour an MLS like model with loops and linkers of 63 to 126 kbp. The predicted large spaces between the chromatin fibres allow typically sized biological molecules to reach nearly every location in the nucleus by moderately obstructed diffusion and is in disagreement with the much simplified assumption that defined channels between territories for molecular transport as in the Interchromosomal Domain (ICD) hypothesis exist and are necessary for transport. All this is also in agreement with recent selective high-resolution chromosome interaction capture (T2C) experiments, the scaling behaviour of the DNA sequence, the dynamics of the chromatin fibre, the diffusion of molecules, and other measurements. Also all other chromosome topologies can in principle be excluded. In summary, polymer simulations of whole nuclei compared to experimental data not only clearly favour only a stable loop aggregate/rosette like genome architecture whose local topology is tightly connected to the global morphology and dynamics of the cell nucleus and hence can be used for understanding genome organization also in respect to diagnosis and treatment. This is in agreement with and also leads to a general novel framework of genome emergence, function, and evolution.

TNF-α-producing macrophages determine subtype identity and prognosis via AP1 enhancer reprogramming in pancreatic cancer

Large-scale genomic profiling of pancreatic cancer (PDAC) has revealed two distinct subtypes: 'classical' and 'basal-like'. Their variable coexistence within the stromal immune microenvironment is linked to differential prognosis; however, the extent to which these neoplastic subtypes shape the stromal immune landscape and impact clinical outcome remains unclear. By combining preclinical models, patient-derived xenografts, as well as FACS-sorted PDAC patient biopsies, we show that the basal-like neoplastic state is sustained via BRD4-mediated cJUN/AP1 expression, which induces CCL2 to recruit tumor necrosis factor (TNF)-α-secreting macrophages. TNF-α+ macrophages force classical neoplastic cells into an aggressive phenotypic state via lineage reprogramming. Integration of ATAC-, ChIP- and RNA-seq data revealed distinct JUNB/AP1 (classical) and cJUN/AP1 (basal-like)-driven regulation of PDAC subtype identity. Pharmacological inhibition of BRD4 led to suppression of the BRD4-cJUN-CCL2-TNF-α axis, restoration of classical subtype identity and a favorable prognosis. Hence, patient-tailored therapy for a cJUNhigh/TNF-αhigh subtype is paramount in overcoming highly inflamed and aggressive PDAC states.

Translational Applications of Linear and Circular Long Noncoding RNAs in Endometriosis

Endometriosis is a chronic gynecologic disease that negatively affects the quality of life of many women. Unfortunately, endometriosis does not have a cure. The current medical treatments involve hormonal manipulation with unwanted side effects and high recurrence rates after stopping the medication. Sadly, a definitive diagnosis for endometriosis requires invasive surgical procedures, with the risk of complications, additional surgeries in the future, and a high rate of recurrence. Both improved therapies and noninvasive diagnostic tests are needed. The unique molecular features of endometriosis have been studied at the coding gene level. While the molecular components of endometriosis at the small RNA level have been studied extensively, other noncoding RNAs, such as long intergenic noncoding RNAs and the more recently discovered subset of long noncoding RNAs called circular RNAs, have been studied more limitedly. This review describes the molecular formation of long noncoding and the unique circumstances of the formation of circular long noncoding RNAs, their expression and function in endometriosis, and promising preclinical studies. Continued translational research on long noncoding RNAs, including the more stable circular long noncoding RNAs, may lead to improved therapeutic and diagnostic opportunities.

Investigation of the Basic Steps in the Chromosome Conformation Capture Procedure

A constellation of chromosome conformation capture methods (С-methods) are an important tool for biochemical analysis of the spatial interactions between DNA regions that are separated in the primary sequence. All these methods are based on the long sequence of basic steps of treating cells, nuclei, chromatin, and finally DNA, thus representing a significant technical challenge. Here, we present an in-depth study of the basic steps in the chromatin conformation capture procedure (3С), which was performed using Drosophila Schneider 2 cells as a model. We investigated the steps of cell lysis, nuclei washing, nucleoplasm extraction, chromatin treatment with SDS/Triton X-100, restriction enzyme digestion, chromatin ligation, reversion of cross-links, DNA extraction, treatment of a 3C library with RNases, and purification of the 3C library. Several options were studied, and optimal conditions were found. Our work contributes to the understanding of the 3C basic steps and provides a useful guide to the 3C procedure.

Targeted chromatin conformation analysis identifies novel distal neural enhancers of ZEB2 in pluripotent stem cell differentiation

The transcription factor zinc finger E-box binding protein 2 (ZEB2) controls embryonic and adult cell fate decisions and cellular maturation in many stem/progenitor cell types. Defects in these processes in specific cell types underlie several aspects of Mowat–Wilson syndrome (MOWS), which is caused by ZEB2 haplo-insufficiency. Human ZEB2, like mouse Zeb2, is located on chromosome 2 downstream of a ±3.5 Mb-long gene-desert, lacking any protein-coding gene. Using temporal targeted chromatin capture (T2C), we show major chromatin structural changes based on mapping in-cis proximities between the ZEB2 promoter and this gene desert during neural differentiation of human-induced pluripotent stem cells, including at early neuroprogenitor cell (NPC)/rosette state, where ZEB2 mRNA levels increase significantly. Combining T2C with histone-3 acetylation mapping, we identified three novel candidate enhancers about 500 kb upstream of the ZEB2 transcription start site. Functional luciferase-based assays in heterologous cells and NPCs reveal co-operation between these three enhancers. This study is the first to document in-cis Regulatory Elements located in ZEB2’s gene desert. The results further show the usability of T2C for future studies of ZEB2 REs in differentiation and maturation of multiple cell types and the molecular characterization of newly identified MOWS patients that lack mutations in ZEB2 protein-coding exons.

Enhancer release and retargeting activates disease-susceptibility genes

The functional engagement between an enhancer and its target promoter ensures precise gene transcription1. Understanding the basis of promoter choice by enhancers has important implications for health and disease. Here we report that functional loss of a preferred promoter can release its partner enhancer to loop to and activate an alternative promoter (or alternative promoters) in the neighbourhood. We refer to this target-switching process as 'enhancer release and retargeting'. Genetic deletion, motif perturbation or mutation, and dCas9-mediated CTCF tethering reveal that promoter choice by an enhancer can be determined by the binding of CTCF at promoters, in a cohesin-dependent manner-consistent with a model of 'enhancer scanning' inside the contact domain. Promoter-associated CTCF shows a lower affinity than that at chromatin domain boundaries and often lacks a preferred motif orientation or a partnering CTCF at the cognate enhancer, suggesting properties distinct from boundary CTCF. Analyses of cancer mutations, data from the GTEx project and risk loci from genome-wide association studies, together with a focused CRISPR interference screen, reveal that enhancer release and retargeting represents an overlooked mechanism that underlies the activation of disease-susceptibility genes, as exemplified by a risk locus for Parkinson's disease (NUCKS1-RAB7L1) and three loci associated with cancer (CLPTM1L-TERT, ZCCHC7-PAX5 and PVT1-MYC).

Transcriptionally inactive hepatitis B virus episome DNA preferentially resides in the vicinity of chromosome 19 in 3D host genome upon infection

The hepatitis B virus (HBV) infects 257 million people worldwide. HBV infection requires establishment and persistence of covalently closed circular (ccc) DNA, a viral episome, in nucleus. Here, we study cccDNA spatial localization in the 3D host genome by using chromosome conformation capture-based sequencing analysis and fluorescence in situ hybridization (FISH). We show that transcriptionally inactive cccDNA is not randomly distributed in host nucleus. Rather, it is preferentially accumulated at specialized areas, including regions close to chromosome 19 (chr.19). Activation of the cccDNA is apparently associated with its re-localization, from a pre-established heterochromatin hub formed by 5 regions of chr.19 to transcriptionally active regions formed by chr.19 and nearby chromosomes including chr.16, 17, 20, and 22. This active versus inactive positioning at discrete regions of the host genome is primarily controlled by the viral HBx protein and by host factors including the structural maintenance of chromosomes protein 5/6 (SMC5/6) complex.

DNA methylation changes during long-term in vitro cell culture are caused by epigenetic drift

Culture expansion of primary cells evokes highly reproducible DNA methylation (DNAm) changes. We have identified CG dinucleotides (CpGs) that become continuously hyper- or hypomethylated during long-term culture of mesenchymal stem cells (MSCs) and other cell types. Bisulfite barcoded amplicon sequencing (BBA-seq) demonstrated that DNAm patterns of neighboring CpGs become more complex without evidence of continuous pattern development and without association to oligoclonal subpopulations. Circularized chromatin conformation capture (4C) revealed reproducible changes in nuclear organization between early and late passages, while there was no enriched interaction with other genomic regions that also harbor culture-associated DNAm changes. Chromatin immunoprecipitation of CTCF did not show significant differences during long-term culture of MSCs, however culture-associated hypermethylation was enriched at CTCF binding sites and hypomethylated CpGs were devoid of CTCF. Taken together, our results support the notion that DNAm changes during culture-expansion are not directly regulated by a targeted mechanism but rather resemble epigenetic drift.

Julia Franzen et al. investigate if changes in DNA methylation at specific genetic loci during cell culture expansion are due to a specific mechanism or gradual deregulation of an epigenetic state. Their results suggest that changes in CpG methylation are due to indirect epigenetic drift, rather than a consequence of targeting by DNA methyltransferases.

Insulator foci distance correlates with cellular and nuclear morphology in early Drosophila embryos

The three-dimensional (3D) organization of the genome is highly dynamic, changing during development and varying across different tissues and cell types. Recent studies indicate that these changes alter regulatory interactions, leading to changes in gene expression. Despite its importance, the mechanisms that influence genomic organization remain poorly understood. We have previously identified a network of chromatin boundary elements, or insulators, in the Drosophila Antennapedia homeotic complex (ANT-C). These genomic elements interact with one another to tether chromatin loops that could block or promote enhancer-promoter interactions. To understand the function of these insulators, we assessed their interactions by measuring their 3D nuclear distance in developing animal tissues. Our data suggest that the ANT-C Hox complex might be in a folded or looped configuration rather than in a random or extended form. The architecture of the ANT-C complex, as read out by the pair-wise distance between insulators, undergoes a strong compression during late embryogenesis, coinciding with the reduction of cell and nuclear diameters due to continued cell divisions in post-cleavage cells. Our results suggest that genomic architecture and gene regulation may be influenced by cellular morphology and movement during development.

miRNA-independent function of long noncoding pri-miRNA loci

Among the large, diverse set of mammalian long noncoding RNAs (lncRNAs), long noncoding primary microRNAs (lnc-pri-miRNAs) are those that host miRNAs. Whether lnc-pri-miRNA loci have important biological function independent of their cognate miRNAs is poorly understood. From a genome-scale lncRNA screen, lnc-pri-miRNA loci were enriched for function in cell proliferation, and in glioblastoma (i.e., GBM) cells with DGCR8 or DROSHA knockdown, lnc-pri-miRNA screen hits still regulated cell growth. To molecularly dissect the function of a lnc-pri-miRNA locus, we studied LOC646329 (also known as MIR29HG), which hosts the miR-29a/b1 cluster. In GBM cells, LOC646329 knockdown reduced miR-29a/b1 levels, and these cells exhibited decreased growth. However, genetic deletion of the miR-29a/b1 cluster (LOC646329-miR29Δ) did not decrease cell growth, while knockdown of LOC646329-miR29Δ transcripts reduced cell proliferation. The miR-29a/b1-independent activity of LOC646329 corresponded to enhancer-like activation of a neighboring oncogene (MKLN1), regulating cell propagation. The LOC646329 locus interacts with the MKLN1 promoter, and antisense oligonucleotide knockdown of the lncRNA disrupts these interactions and reduces the enhancer-like activity. More broadly, analysis of genome-wide data from multiple human cell types showed that lnc-pri-miRNA loci are significantly enriched for DNA looping interactions with gene promoters as well as genomic and epigenetic characteristics of transcriptional enhancers. Functional studies of additional lnc-pri-miRNA loci demonstrated cognate miRNA-independent enhancer-like activity. Together, these data demonstrate that lnc-pri-miRNA loci can regulate cell biology via both miRNA-dependent and miRNA-independent mechanisms.

HOXBLINC long non-coding RNA activation promotes leukemogenesis in NPM1-mutant acute myeloid leukemia

Nucleophosmin (NPM1) is the most commonly mutated gene in acute myeloid leukemia (AML) resulting in aberrant cytoplasmic translocation of the encoded nucleolar protein (NPM1c⁺). NPM1c⁺ maintains a unique leukemic gene expression program, characterized by activation of HOXA/B clusters and MEIS1 oncogene to facilitate leukemogenesis. However, the mechanisms by which NPM1c⁺ controls such gene expression patterns to promote leukemogenesis remain largely unknown. Here, we show that the activation of HOXBLINC, a HOXB locus-associated long non-coding RNA (lncRNA), is a critical downstream mediator of NPM1c⁺-associated leukemic transcription program and leukemogenesis. HOXBLINC loss attenuates NPM1c⁺-driven leukemogenesis by rectifying the signature of NPM1c⁺ leukemic transcription programs. Furthermore, overexpression of HoxBlinc (HoxBlincTg) in mice enhances HSC self-renewal and expands myelopoiesis, leading to the development of AML-like disease, reminiscent of the phenotypes seen in the Npm1 mutant knock-in (Npm1^c/+) mice. HoxBlincTg and Npm1^c/+ HSPCs share significantly overlapped transcriptome and chromatin structure. Mechanistically, HoxBlinc binds to the promoter regions of NPM1c⁺ signature genes to control their activation in HoxBlincTg HSPCs, via MLL1 recruitment and promoter H3K4me3 modification. Our study reveals that HOXBLINC lncRNA activation plays an essential oncogenic role in NPM1c⁺ leukemia. HOXBLINC and its partner MLL1 are potential therapeutic targets for NPM1c⁺ AML.

Nucleophosmin (NPM1) gene mutation induces a specific gene expression program leading to acute myeloid leukaemia. Here, the authors show that mutant NPM1 activates a HOXB locus-associated long non-coding RNA which is essential for its associated oncogenic transcriptional program and leukaemia development.

Interplay between FLI-1 and the LDB1 complex in murine erythroleukemia cells and during megakaryopoiesis

Transcription factors are key players in a broad range of cellular processes such as cell-fate decision. Understanding how they act to control these processes is of critical importance for therapy purposes. FLI-1 controls several hematopoietic lineage differentiation including megakaryopoiesis and erythropoiesis. Its aberrant expression is often observed in cancer and is associated with poor prognosis. We showed that FLI-1 interacts with the LDB1 complex, which also plays critical roles in erythropoiesis and megakaryopoiesis. In this study, we aimed to unravel how FLI-1 and the LDB1 complex act together in murine erythroleukemia cells and in megakaryocyte. Combining omics techniques, we show that FLI-1 enables the recruitment of the LDB1 complex to regulatory sequences of megakaryocytic genes and to enhancers. We show as well for the first time that FLI-1 is able to modulate the 3D chromatin organization by promoting chromatin looping between enhancers and promoters most likely through the LDB1 complex.

•

FLI-1 is important for the recruitment of the LDB1 complex

•

FLI-1 is important for chromatin looping

•

FLI-1 and the LDB1 complex co-regulate megakaryopoiesis

Ultrasound May Suppress Tumor Growth, Inhibit Inflammation, and Establish Tolerogenesis by Remodeling Innatome via Pathways of ROS, Immune Checkpoints, Cytokines, and Trained Immunity/Tolerance

BACKGROUND

The immune mechanisms underlying low-intensity ultrasound- (LIUS-) mediated suppression of inflammation and tumorigenesis remain poorly determined.

METHODS

We used microarray datasets from the NCBI GEO DataSet repository and conducted comprehensive data-mining analyses, where we examined the gene expression of 1376 innate immune regulators (innatome genes (IGs) in cells treated with LIUS.

RESULTS

We made the following findings: (1) LIUS upregulates proinflammatory IGs and downregulates metastasis genes in cancer cells, and LIUS upregulates adaptive immunity pathways but inhibits danger-sensing and inflammation pathways and promote tolerogenic differentiation in bone marrow (BM) cells. (2) LIUS upregulates IGs encoded for proteins localized in the cytoplasm, extracellular space, and others, but downregulates IG proteins localized in nuclear and plasma membranes, and LIUS downregulates phosphatases. (3) LIUS-modulated IGs act partially via several important pathways of reactive oxygen species (ROS), reverse signaling of immune checkpoint receptors B7-H4 and BTNL2, inflammatory cytokines, and static or oscillatory shear stress and heat generation, among which ROS is a dominant mechanism. (4) LIUS upregulates trained immunity enzymes in lymphoma cells and downregulates trained immunity enzymes and presumably establishes trained tolerance in BM cells. (5) LIUS modulates chromatin long-range interactions to differentially regulate IGs expression in cancer cells and noncancer cells.

CONCLUSIONS

Our analysis suggests novel molecular mechanisms that are utilized by LIUS to induce tumor suppression and inflammation inhibition. Our findings may lead to development of new treatment protocols for cancers and chronic inflammation.

A regulatory sub-circuit downstream of Wnt signaling controls developmental transitions in neural crest formation

The process of cell fate commitment involves sequential changes in the gene expression profiles of embryonic progenitors. This is exemplified in the development of the neural crest, a migratory stem cell population derived from the ectoderm of vertebrate embryos. During neural crest formation, cells transition through distinct transcriptional states in a stepwise manner. The mechanisms underpinning these shifts in cell identity are still poorly understood. Here we employ enhancer analysis to identify a genetic sub-circuit that controls developmental transitions in the nascent neural crest. This sub-circuit links Wnt target genes in an incoherent feedforward loop that controls the sequential activation of genes in the neural crest lineage. By examining the cis-regulatory apparatus of Wnt effector gene AXUD1, we found that multipotency factor SP5 directly promotes neural plate border identity, while inhibiting premature expression of specification genes. Our results highlight the importance of repressive interactions in the neural crest gene regulatory network and illustrate how genes activated by the same upstream signal become temporally segregated during progressive fate restriction.

A genetic variant controls interferon-β gene expression in human myeloid cells by preventing C/EBP-β binding on a conserved enhancer

Interferon β (IFN-β) is a cytokine that induces a global antiviral proteome, and regulates the adaptive immune response to infections and tumors. Its effects strongly depend on its level and timing of expression. Therefore, the transcription of its coding gene IFNB1 is strictly controlled. We have previously shown that in mice, the TRIM33 protein restrains Ifnb1 transcription in activated myeloid cells through an upstream inhibitory sequence called ICE. Here, we show that the deregulation of Ifnb1 expression observed in murine Trim33^-/- macrophages correlates with abnormal looping of both ICE and the Ifnb1 gene to a 100 kb downstream region overlapping the Ptplad2/Hacd4 gene. This region is a predicted myeloid super-enhancer in which we could characterize 3 myeloid-specific active enhancers, one of which (E5) increases the response of the Ifnb1 promoter to activation. In humans, the orthologous region contains several single nucleotide polymorphisms (SNPs) known to be associated with decreased expression of IFNB1 in activated monocytes, and loops to the IFNB1 gene. The strongest association is found for the rs12553564 SNP, located in the E5 orthologous region. The minor allele of rs12553564 disrupts a conserved C/EBP-β binding motif, prevents binding of C/EBP-β, and abolishes the activation-induced enhancer activity of E5. Altogether, these results establish a link between a genetic variant preventing binding of a transcription factor and a higher order phenotype, and suggest that the frequent minor allele (around 30% worldwide) might be associated with phenotypes regulated by IFN-β expression in myeloid cells.

Genome-wide association studies identify multiple genetic variants associated with higher order phenotypes. Pinpointing the causative variant and understanding its molecular mode of action is a complex task. Using a murine model of interferon-β transcriptional deregulation, we characterize a super-enhancer controlling Ifnb1 expression in myeloid cells. The most active enhancer of this locus is conserved in humans, but presents a frequent variant found in around 30% of the population worldwide. This variant prevents binding of the C/EBP-β transcription factor, and is associated with decreased expression of IFNB1 in activated monocytes. When mimicked in the murine enhancer, it abolishes its inducible enhancer activity. Our results describe the molecular link between a point mutation and a cellular phenotype that could influence clinical situations.

Defining Essential Enhancers for Pluripotent Stem Cells Using a Features-Oriented CRISPR-Cas9 Screen

cis-regulatory elements (CREs) regulate the expression of genes in their genomic neighborhoods and influence cellular processes such as cell-fate maintenance and differentiation. To date, there remain major gaps in the functional characterization of CREs and the identification of their target genes in the cellular native environment. In this study, we perform a features-oriented CRISPR-utilized systematic (FOCUS) screen of OCT4-bound CREs using CRISPR-Cas9 to identify functional enhancers important for pluripotency maintenance in mESCs. From the initial 235 candidates tested, 16 CREs are identified to be essential stem cell enhancers. Using RNA-seq and genomic 4C-seq, we further uncover a complex network of candidate CREs and their downstream target genes, which supports the growth and self-renewal of mESCs. Notably, an essential enhancer, CRE111, and its target, Lrrc31, form the important switch to modulate the LIF-JAK1-STAT3 signaling pathway.

Global epigenomic analysis of KSHV-infected primary effusion lymphoma identifies functional MYC superenhancers and enhancer RNAs

Enhancers play indispensable roles in cell proliferation and survival through spatiotemporally regulating gene transcription. Active enhancers and superenhancers often produce noncoding enhancer RNAs (eRNAs) that precisely control RNA polymerase II activity. Kaposi's sarcoma-associated herpesvirus (KSHV) is a human oncogenic gamma-2 herpesvirus that causes Kaposi's sarcoma and primary effusion lymphoma (PEL). It is well characterized that KSHV utilizes host epigenetic machineries to control the switch between two lifecycles, latency and lytic replication. However, how KSHV impacts host epigenome at different stages of viral lifecycle is not well understood. Using global run-on sequencing (GRO-seq) and chromatin-immunoprecipitation sequencing (ChIP-seq), we profiled the dynamics of host transcriptional regulatory elements during latency and lytic replication of KSHV-infected PEL cells. This revealed that a number of critical host genes for KSHV latency, including MYC proto-oncogene, were under the control of superenhancers whose activities were globally repressed upon viral reactivation. The eRNA-expressing MYC superenhancers were located downstream of the MYC gene in KSHV-infected PELs and played a key role in MYC expression. RNAi-mediated depletion or dCas9-KRAB CRISPR inhibition of eRNA expression significantly reduced MYC mRNA level in PELs, as did the treatment of an epigenomic drug that globally blocks superenhancer function. Finally, while cellular IRF4 acted upon eRNA expression and superenhancer function for MYC expression during latency, KSHV viral IRF4 repressed cellular IRF4 expression, decreasing MYC expression and thereby, facilitating lytic replication. These results indicate that KSHV acts as an epigenomic driver that modifies host epigenomic status upon reactivation by effectively regulating host enhancer function.

The Myb-MuvB Complex Is Required for YAP-Dependent Transcription of Mitotic Genes

YAP and TAZ, downstream effectors of the Hippo pathway, are important regulators of proliferation. Here, we show that the ability of YAP to activate mitotic gene expression is dependent on the Myb-MuvB (MMB) complex, a master regulator of genes expressed in the G2/M phase of the cell cycle. By carrying out genome-wide expression and binding analyses, we found that YAP promotes binding of the MMB subunit B-MYB to the promoters of mitotic target genes. YAP binds to B-MYB and stimulates B-MYB chromatin association through distal enhancer elements that interact with MMB-regulated promoters through chromatin looping. The cooperation between YAP and B-MYB is critical for YAP-mediated entry into mitosis. Furthermore, the expression of genes coactivated by YAP and B-MYB is associated with poor survival of cancer patients. Our findings provide a molecular mechanism by which YAP and MMB regulate mitotic gene expression and suggest a link between two cancer-relevant signaling pathways.

Exonuclease combinations reduce noises in 3D genomics technologies

Chromosome conformation-capture technologies are widely used in 3D genomics; however, experimentally, such methods have high-noise limitations and, therefore, require significant bioinformatics efforts to extract reliable distal interactions. Miscellaneous undesired linear DNAs, present during proximity-ligation, represent a main noise source, which needs to be minimized or eliminated. In this study, different exonuclease combinations were tested to remove linear DNA fragments from a circularized DNA preparation. This method efficiently removed linear DNAs, raised the proportion of annulation and increased the valid-pairs ratio from ∼40% to ∼80% for enhanced interaction detection in standard Hi-C. This strategy is applicable for development of various 3D genomics technologies, or optimization of Hi-C sequencing efficiency.

Analysis of 4C-seq data: A comparison of methods

The circular chromosome conformation capture technique followed by sequencing (4C-seq) has been used in a number of studies to investigate chromosomal interactions between DNA fragments. Computational pipelines have been developed and published that offer various possibilities of 4C-seq data processing and statistical analysis. Here, we present an overview of four of such pipelines (fourSig, FourCSeq, 4C-ker and w4Cseq) taking into account the most important stages of computations. We provide comparisons of the methods and discuss their advantages and possible weaknesses. We illustrate the results with the use of data obtained for two different species, in a study devoted to vernalization control in Arabidopsis thaliana by the FLOWERING LOCUS C (FLC) gene and to long-range chromatin interactions in mouse embryonic stem cells.

Distinct IL-1α-responsive enhancers promote acute and coordinated changes in chromatin topology in a hierarchical manner

How cytokine-driven changes in chromatin topology are converted into gene regulatory circuits during inflammation still remains unclear. Here, we show that interleukin (IL)-1α induces acute and widespread changes in chromatin accessibility via the TAK1 kinase and NF-κB at regions that are highly enriched for inflammatory disease-relevant SNPs. Two enhancers in the extended chemokine locus on human chromosome 4 regulate the IL-1α-inducible IL8 and CXCL1-3 genes. Both enhancers engage in dynamic spatial interactions with gene promoters in an IL-1α/TAK1-inducible manner. Microdeletions of p65-binding sites in either of the two enhancers impair NF-κB recruitment, suppress activation and biallelic transcription of the IL8/CXCL2 genes, and reshuffle higher-order chromatin interactions as judged by i4C interactome profiles. Notably, these findings support a dominant role of the IL8 "master" enhancer in the regulation of sustained IL-1α signaling, as well as for IL-8 and IL-6 secretion. CRISPR-guided transactivation of the IL8 locus or cross-TAD regulation by TNFα-responsive enhancers in a different model locus supports the existence of complex enhancer hierarchies in response to cytokine stimulation that prime and orchestrate proinflammatory chromatin responses downstream of NF-κB.

Enhancers as regulators of antigen receptor loci three-dimensional chromatin structure

Enhancers are defined as regulatory elements that control transcription in a cell-type and developmental stage-specific manner. They achieve this by physically interacting with their cognate gene promoters. Significantly, these interactions can occur through long genomic distances since enhancers may not be near their cognate promoters. The optimal coordination of enhancer-regulated transcription is essential for the function and identity of the cell. Although great efforts to fully understand the principles of this type of regulation are ongoing, other potential functions of the long-range chromatin interactions (LRCIs) involving enhancers are largely unexplored. We recently uncovered a new role for enhancer elements in determining the three-dimensional (3D) structure of the immunoglobulin kappa (Igκ) light chain receptor locus suggesting a structural function for these DNA elements. This enhancer-mediated locus configuration shapes the resulting Igκ repertoire. We also propose a role for enhancers as critical components of sub-topologically associating domain (subTAD) formation and nuclear spatial localization.

A 16q22.1 variant confers susceptibility to colorectal cancer as a distal regulator of ZFP90

Genome-wide association studies (GWASs) implicate 16q22.1 locus in risk for colorectal cancer (CRC). However, the underlying oncogenic mechanisms remain unknown. Here, through comprehensive filtration, we prioritized rs7198799, a common SNP in the second intron of the CDH1, as the putative causal variant. In addition, we found an association of CRC-risk allele C of rs7198799 with elevated transcript level of biological plausible candidate gene ZFP90 via expression quantitative trait loci analysis. Mechanistically, causal variant rs7198799 resides in an enhancer element and remotely regulate ZFP90 expression by targeting the transcription factor NFATC2. Remarkably, CRISPR/Cas9-guided single-nucleotide editing demonstrated the direct effect of rs7198799 on ZFP90 expression and CRC cellular malignant phenotype. Furthermore, ZFP90 affects several oncogenic pathways, including BMP4, and promotes carcinogenesis in patients and in animal models with ZFP90 specific genetic manipulation. Taken together, these findings reveal a risk SNP-mediated long-range regulation on the NFATC2-ZFP90-BMP4 pathway underlying the initiation of CRC.

The transcription factor PU.1 mediates enhancer-promoter looping that is required for IL-1β eRNA and mRNA transcription in mouse melanoma and macrophage cell lines

The DNA-binding protein PU.1 is a myeloid lineage-determining and pioneering transcription factor due to its ability to bind "closed" genomic sites and maintain "open" chromatin state for myeloid lineage-specific genes. The precise mechanism of PU.1 in cell type-specific programming is yet to be elucidated. The melanoma cell line B16BL6, although it is nonmyeloid lineage, expressed Toll-like receptors and activated the transcription factor NF-κB upon stimulation by the bacterial cell wall component lipopolysaccharide. However, it did not produce cytokines, such as IL-1β mRNA. Ectopic PU.1 expression induced remodeling of a novel distal enhancer (located ∼10 kbp upstream of the IL-1β transcription start site), marked by nucleosome depletion, enhancer-promoter looping, and histone H3 lysine 27 acetylation (H3K27ac). PU.1 induced enhancer-promoter looping and H3K27ac through two distinct PU.1 regions. These PU.1-dependent events were independently required for subsequent signal-dependent and co-dependent events: NF-κB recruitment and further H3K27ac, both of which were required for enhancer RNA (eRNA) transcription. In murine macrophage RAW264.7 cells, these PU.1-dependent events were constitutively established and readily expressed eRNA and subsequently IL-1β mRNA by lipopolysaccharide stimulation. In summary, this study showed a sequence of epigenetic events in programming IL-1β transcription by the distal enhancer priming and eRNA production mediated by PU.1 and the signal-dependent transcription factor NF-κB.

Interplay between FACT subunit SPT16 and TRIM33 can remodel chromatin at macrophage distal regulatory elements

BACKGROUND

Cell type-specific use of cis-acting regulatory elements is mediated by the combinatorial activity of transcription factors involved in lineage determination and maintenance of cell identity. In macrophages, specific transcriptional programs are dictated by the transcription factor PU.1 that primes distal regulatory elements for macrophage identities and makes chromatin competent for activity of stimuli-dependent transcription factors. Although the advances in genome-wide approaches have elucidated the functions of these macrophage-specific distal regulatory elements in transcriptional responses, chromatin structures associated with PU.1 priming and the underlying mechanisms of action of these cis-acting sequences are not characterized.

RESULTS

Here, we show that, in macrophages, FACT subunit SPT16 can bind to positioned nucleosomes directly flanking PU.1-bound sites at previously uncharacterized distal regulatory elements located near genes essential for macrophage development and functions. SPT16 can interact with the transcriptional co-regulator TRIM33 and binds to half of these sites in a TRIM33-dependent manner. Using the Atp1b3 locus as a model, we show that FACT binds to two positioned nucleosomes surrounding a TRIM33/PU.1-bound site in a region, located 35 kb upstream the Atp1b3 TSS, that interact with the Atp1b3 promoter. At this - 35 kb region, TRIM33 deficiency leads to FACT release, loss of the two positioned nucleosomes, RNA Pol II recruitment and bidirectional transcription. These modifications are associated with higher levels of FACT binding at the Atp1b3 promoter, an increase of RNA Pol II recruitment and an increased expression of Atp1b3 in Trim33^-/- macrophages.

CONCLUSIONS

Thus, sequestering of SPT16/FACT by TRIM33 at PU.1-bound distal regions might represent a new regulatory mechanism for RNA Pol II recruitment and transcription output in macrophages.

Male-Specific Long Noncoding RNA TTTY15 Inhibits Non-Small Cell Lung Cancer Proliferation and Metastasis via TBX4

Gender affects cancer susceptibility. Currently, there are only a few studies on Y chromosome-linked long noncoding RNAs (lncRNAs), and the potential association between lncRNAs and cancers in males has not been fully elucidated. Here, we examined the expression of testis-specific transcript Y-linked 15 (TTTY15) in 37 males with non-small cell lung cancer (NSCLC), and performed circular chromosome conformation capture with next-generation sequencing to determine the genomic interaction regions of the TTTY15 gene. Our results showed that the expression levels of TTTY15 were lower in NSCLC tissues. Lower TTTY15 expression levels were associated with Tumor-Node-Metastasis (TNM) stage. A TTTY15 knockdown promoted malignant transformation of NSCLC cells. Based on the bioinformatics analysis of circular chromosome conformation capture data, we found that T-box transcription factor 4 (TBX4) may be a potential target gene of TTTY15. The RNA immunoprecipitation and chromatin immunoprecipitation results showed that TTTY15 may interact with DNA (cytosine-5)-methyltransferase 3A (DNMT3A), and the TTTY15 knockdown increased the binding of DNMT3A to the TBX4 promoter. We concluded that low TTTY15 expression correlates with worse prognosis among patients with NSCLC. TTTY15 promotes TBX4 expression via DNMT3A-mediated regulation. The identification of lncRNAs encoded by male-specific genes may help to identify potential targets for NSCLC therapy.

Experimental Data-Mining Analyses Reveal New Roles of Low-Intensity Ultrasound in Differentiating Cell Death Regulatome in Cancer and Non-cancer Cells via Potential Modulation of Chromatin Long-Range Interactions

Background: The mechanisms underlying low intensity ultrasound (LIUS) mediated suppression of inflammation and tumorigenesis remain poorly determined. Methods: We used microarray datasets from NCBI GEO Dataset databases and conducted a comprehensive data mining analyses, where we studied the gene expression of 299 cell death regulators that regulate 13 different cell death types (cell death regulatome) in cells treated with LIUS. Results: We made the following findings: (1) LIUS exerts a profound effect on the expression of cell death regulatome in cancer cells and non-cancer cells. Of note, LIUS has the tendency to downregulate the gene expression of cell death regulators in non-cancer cells. Most of the cell death regulator genes downregulated by LIUS in non-cancer cells are responsible for mediating inflammatory signaling pathways; (2) LIUS activates different cell death transcription factors in cancer and non-cancer cells. Transcription factors TP-53 and SRF- were induced by LIUS exposure in cancer cells and non-cancer cells, respectively; (3) As two well-accepted mechanisms of LIUS, mild hyperthermia and oscillatory shear stress induce changes in the expression of cell death regulators, therefore, may be responsible for inducing LIUS mediated changes in gene expression patterns of cell death regulators in cells; (4) LIUS exposure may change the redox status of the cells. LIUS may induce more of antioxidant effects in non-cancer cells compared to cancer cells; and (5) The genes modulated by LIUS in cancer cells have distinct chromatin long range interaction (CLRI) patterns to that of non-cancer cells. Conclusions: Our analysis suggests novel molecular mechanisms that may be utilized by LIUS to induce tumor suppression and inflammation inhibition. Our findings may lead to development of new treatment protocols for cancers and chronic inflammation.

GemC1 governs multiciliogenesis through direct interaction with and transcriptional regulation of p73

A distinct combination of transcription factors elicits the acquisition of a specific fate and the initiation of a differentiation program. Multiciliated cells (MCCs) are a specialized type of epithelial cells that possess dozens of motile cilia on their apical surface. Defects in cilia function have been associated with ciliopathies that affect many organs, including brain and airway epithelium. Here we show that the geminin coiled-coil domain-containing protein 1 GemC1 (also known as Lynkeas) regulates the transcriptional activation of p73, a transcription factor central to multiciliogenesis. Moreover, we show that GemC1 acts in a trimeric complex with transcription factor E2F5 and tumor protein p73 (officially known as TP73), and that this complex is important for the activation of the p73 promoter. We also provide in vivo evidence that GemC1 is necessary for p73 expression in different multiciliated epithelia. We further show that GemC1 regulates multiciliogenesis through the control of chromatin organization, and the epigenetic marks/tags of p73 and Foxj 1. Our results highlight novel signaling cues involved in the commitment program of MCCs across species and tissues.This article has an associated First Person interview with the first author of the paper.

Increased acetylation of H3K14 in the genomic regions that encode trained immunity enzymes in lysophosphatidylcholine-activated human aortic endothelial cells - Novel qualification markers for chronic disease risk factors and conditional DAMPs

To test our hypothesis that proatherogenic lysophosphatidylcholine (LPC) upregulates trained immunity pathways (TIPs) in human aortic endothelial cells (HAECs), we conducted an intensive analyses on our RNA-Seq data and histone 3 lysine 14 acetylation (H3K14ac)-CHIP-Seq data, both performed on HAEC treated with LPC. Our analysis revealed that: 1) LPC induces upregulation of three TIPs including glycolysis enzymes (GE), mevalonate enzymes (ME), and acetyl-CoA generating enzymes (ACE); 2) LPC induces upregulation of 29% of 31 histone acetyltransferases, three of which acetylate H3K14; 3) LPC induces H3K14 acetylation (H3K14ac) in the genomic DNA that encodes LPC-induced TIP genes (79%) in comparison to that of in LPC-induced effector genes (43%) including ICAM-1; 4) TIP pathways are significantly different from that of EC activation effectors including adhesion molecule ICAM-1; 5) reactive oxygen species generating enzyme NOX2 deficiency decreases, but antioxidant transcription factor Nrf2 deficiency increases, the expressions of a few TIP genes and EC activation effector genes; and 6) LPC induced TIP genes(81%) favor inter-chromosomal long-range interactions (CLRI, trans-chromatin interaction) while LPC induced effector genes (65%) favor intra-chromosomal CLRIs (cis-chromatin interaction). Our findings demonstrated that proatherogenic lipids upregulate TIPs in HAECs, which are a new category of qualification markers for chronic disease risk factors and conditional DAMPs and potential mechanisms for acute inflammation transition to chronic ones. These novel insights may lead to identifications of new cardiovascular risk factors in upregulating TIPs in cardiovascular cells and novel therapeutic targets for the treatment of metabolic cardiovascular diseases, inflammation, and cancers. (total words: 245).

Increasing Upstream Chromatin Long-Range Interactions May Favor Induction of Circular RNAs in LysoPC-Activated Human Aortic Endothelial Cells

Circular RNAs (circRNAs) are non-coding RNAs that form covalently closed continuous loops, and act as gene regulators in physiological and disease conditions. To test our hypothesis that proatherogenic lipid lysophosphatidylcholine (LPC) induce a set of circRNAs in human aortic endothelial cell (HAEC) activation, we performed circRNA analysis by searching our RNA-Seq data from LPC-activated HAECs, and found: (1) LPC induces significant modulation of 77 newly characterized cirRNAs, among which 47 circRNAs (61%) are upregulated; (2) 34 (72%) out of 47 upregulated circRNAs are upregulated when the corresponding mRNAs are downregulated, suggesting that the majority of circRNAs are upregulated presumably via LPC-induced “abnormal splicing” when the canonical splicing for generation of corresponding mRNAs is suppressed; (3) Upregulation of 47 circRNAs is temporally associated with mRNAs-mediated LPC-upregulated cholesterol synthesis-SREBP2 pathway and LPC-downregulated TGF-β pathway; (4) Increase in upstream chromatin long-range interaction sites to circRNA related genes is associated with preferred circRNA generation over canonical splicing for mRNAs, suggesting that shifting chromatin long-range interaction sites from downstream to upstream may promote induction of a list of circRNAs in lysoPC-activated HAECs; (5) Six significantly changed circRNAs may have sponge functions for miRNAs; and (6) 74% significantly changed circRNAs contain open reading frames, suggesting that putative short proteins may interfere with the protein interaction-based signaling. Our findings have demonstrated for the first time that a new set of LPC-induced circRNAs may contribute to homeostasis in LPC-induced HAEC activation. These novel insights may lead to identifications of new therapeutic targets for treating metabolic cardiovascular diseases, inflammations, and cancers.

The mouse KLF1 Nan variant impairs nuclear condensation and erythroid maturation

Krüppel-like factor 1 (KLF1) is an essential transcription factor for erythroid development, as demonstrated by Klf1 knockout mice which die around E14 due to severe anemia. In humans, >140 KLF1 variants, causing different erythroid phenotypes, have been described. The KLF1 Nan variant, a single amino acid substitution (p.E339D) in the DNA binding domain, causes hemolytic anemia and is dominant over wildtype KLF1. Here we describe the effects of the KLF1 Nan variant during fetal development. We show that Nan embryos have defects in erythroid maturation. RNA-sequencing of the KLF1 Nan fetal liver cells revealed that Exportin 7 (Xpo7) was among the 782 deregulated genes. This nuclear exportin is implicated in terminal erythroid differentiation; in particular it is involved in nuclear condensation. Indeed, KLF1 Nan fetal liver cells had larger nuclei and reduced chromatin condensation. Knockdown of XPO7 in wildtype erythroid cells caused a similar phenotype. We propose that reduced expression of XPO7 is partially responsible for the erythroid defects observed in KLF1 Nan erythroid cells.

Dismissal of RNA Polymerase II Underlies a Large Ligand-Induced Enhancer Decommissioning Program

Nuclear receptors induce both transcriptional activation and repression programs responsible for development, homeostasis, and disease. Here, we report a previously overlooked enhancer decommissioning strategy underlying a large estrogen receptor alpha (ERα)-dependent transcriptional repression program. The unexpected signature for this E₂-induced program resides in indirect recruitment of ERα to a large cohort of pioneer factor basally active FOXA1-bound enhancers that lack cognate ERα DNA-binding elements. Surprisingly, these basally active estrogen-repressed (BAER) enhancers are decommissioned by ERα-dependent recruitment of the histone demethylase KDM2A, functioning independently of its demethylase activity. Rather, KDM2A tethers the E3 ubiquitin-protein ligase NEDD4 to ubiquitylate/dismiss Pol II to abrogate eRNA transcription, with consequent target gene downregulation. Thus, our data reveal that Pol II ubiquitylation/dismissal may serve as a potentially broad strategy utilized by indirectly bound nuclear receptors to abrogate large programs of pioneer factor-mediated, eRNA-producing enhancers.

Identification of rheumatoid arthritis causal genes using functional genomics

Over the past decade, genome-wide association studies have contributed a wealth of knowledge to our understanding of polygenic disorders such as rheumatoid arthritis. As the size of sample cohorts has improved so too have the computational and experimental methods used to robustly define variants associated with disease susceptibility. The challenge now remains to translate these findings into improved understanding of disease aetiology and patient care. Whilst much of the focus of translating the findings of genome-wide association studies has been on global analysis of all variants identified, careful functional study of individual disease susceptibility loci will be required in order to refine our understanding of how individual variants contribute to disease risk. Here, we present the argument behind such an approach and describe some of the novel tools being used to investigate risk loci. This includes the use of chromosomal conformation capture techniques and modifications of the CRISPR-Cas9 system, with several examples of their implementation being described.

The structural and functional roles of CTCF in the regulation of cell type-specific and human disease-associated super-enhancers

BACKGROUND

Super-enhancers play critical roles in cell-type specific gene controls and human disease progression. CCCTC-binding factor (CTCF), a transcriptional repressor that insulates the expression of neighboring genes and is involved in chromatin interactions, is frequently present in the boundary regions of or within super-enhancers. However, the structural and functional roles of CTCF in regulating super-enhancers remain elusive.

OBJECTIVE

To provide a comprehensive review describing the distinct chromatin features and functional roles of CTCF within super-enhancers.

METHODS

This review compares the various tools used to study the three-dimensional (3D) chromatin architecture of super-enhancers; summarizes the chromatin features of CTCF within cell-type specific super-enhancers and their in vivo biological activities, as determined by CRISPR/Cas9 genome editing; and describes the structural and functional activities of CTCF within human disease-associated super-enhancers.

CONCLUSION

This review provides fundamental insights into the regulatory mechanisms of super-enhancers and facilitates studies of tissue-specific developmental processes and human disease progression.

A B-Cell-Specific Enhancer Orchestrates Nuclear Architecture to Generate a Diverse Antigen Receptor Repertoire

The genome is organized into topologically associated domains (TADs) that enclose smaller subTADs. Here, we identify and characterize an enhancer that is located in the middle of the V gene region of the immunoglobulin kappa light chain (Igκ) locus that becomes active preceding the stage at which this locus undergoes V(D)J recombination. This enhancer is a hub of long-range chromatin interactions connecting subTADs in the V gene region with the recombination center at the J genes. Deletion of this element results in a highly altered long-range chromatin interaction pattern across the locus and, importantly, affects individual V gene utilization locus-wide. These results indicate the existence of an enhancer-dependent framework in the Igκ locus and further suggest that the composition of the diverse antibody repertoire is regulated in a subTAD-specific manner. This enhancer thus plays a structural role in orchestrating the proper folding of the Igκ locus in preparation for V(D)J recombination.

HOXA9 Reprograms the Enhancer Landscape to Promote Leukemogenesis

Aberrant expression of HOXA9 is a prominent feature of acute leukemia driven by diverse oncogenes. Here we show that HOXA9 overexpression in myeloid and B progenitor cells leads to significant enhancer reorganizations with prominent emergence of leukemia-specific de novo enhancers. Alterations in the enhancer landscape lead to activation of an ectopic embryonic gene program. We show that HOXA9 functions as a pioneer factor at de novo enhancers and recruits CEBPα and the MLL3/MLL4 complex. Genetic deletion of MLL3/MLL4 blocks histone H3K4 methylation at de novo enhancers and inhibits HOXA9/MEIS1-mediated leukemogenesis in vivo. These results suggest that therapeutic targeting of HOXA9-dependent enhancer reorganization can be an effective therapeutic strategy in acute leukemia with HOXA9 overexpression.