Runx1 shapes the chromatin landscape via a cascade of direct and indirect targets

Shared and distinct interactions of type 1 and type 2 Epstein-Barr Nuclear Antigen 2 with the human genome

BACKGROUND

There are two major genetic types of Epstein-Barr Virus (EBV): type 1 (EBV-1) and type 2 (EBV-2). EBV functions by manipulating gene expression in host B cells, using virus-encoded gene regulatory proteins including Epstein-Barr Nuclear Antigen 2 (EBNA2). While type 1 EBNA2 is known to interact with human transcription factors (hTFs) such as RBPJ, EBF1, and SPI1 (PU.1), type 2 EBNA2 shares only ~ 50% amino acid identity with type 1 and thus may have distinct binding partners, human genome binding locations, and functions.

RESULTS

In this study, we examined genome-wide EBNA2 binding in EBV-1 and EBV-2 transformed human B cells to identify shared and unique EBNA2 interactions with the human genome, revealing thousands of type-specific EBNA2 ChIP-seq peaks. Computational predictions based on hTF motifs and subsequent ChIP-seq experiments revealed that both type 1 and 2 EBNA2 co-occupy the genome with SPI1 and AP-1 (BATF and JUNB) hTFs. However, type 1 EBNA2 showed preferential co-occupancy with EBF1, and type 2 EBNA2 preferred RBPJ. These differences in hTF co-occupancy revealed possible mechanisms underlying type-specific gene expression of known EBNA2 human target genes: MYC (shared), CXCR7 (type 1 specific), and CD21 (type 2 specific). Both type 1 and 2 EBNA2 binding events were enriched at systemic lupus erythematosus (SLE) and multiple sclerosis (MS) risk loci, while primary biliary cholangitis (PBC) risk loci were specifically enriched for type 2 peaks.

CONCLUSIONS

This study reveals extensive type-specific EBNA2 interactions with the human genome, possible differences in EBNA2 interaction partners, and a possible new role for type 2 EBNA2 in autoimmune disorders. Our results highlight the importance of considering EBV type in the control of human gene expression and disease-related investigations.

Cooperative Response to Endocardial Notch Reveals Interaction With Hippo Pathway

BACKGROUND

The endocardium is a crucial signaling center for cardiac valve development and maturation. Genetic analysis has identified several human endocardial genes whose inactivation leads to bicuspid aortic valve formation and calcific aortic valve disease, but knowledge is very limited about the role played in valve development and disease by noncoding endocardial regulatory regions and upstream factors.

METHODS

We manipulated Notch signaling in mouse embryonic endocardial cells by short-term and long-term coculture with OP9 stromal cells expressing Notch ligands and inhibition of Notch activity. We examined the transcriptional profile and chromatin accessibility landscape for each condition, integrated transcriptomic, transcription factor occupancy, chromatin accessibility, and proteomic datasets. We generated in vitro and in vivo models with CRISPR-Cas9-edited deletions of various noncoding regulatory elements and validated their regulatory potential.

RESULTS

We identified primary and secondary transcriptional responses to Notch ligands in the mouse embryonic endocardium, and a NOTCH-dependent transcriptional signature in valve development and disease. By defining the changes in the chromatin accessibility landscape and integrating with the landscape in developing mouse endocardium and adult human valves, we identify potential noncoding regulatory elements, validated selected candidates, propose interacting cofactors, and define the timeframe of their regulatory activity. Additionally, we found cooperative transcriptional repression with Hippo pathway by inhibiting nuclear Yap (Yes-associated protein) activity in the endocardium during cardiac valve development.

CONCLUSIONS

Sequential Notch-dependent transcriptional regulation in the embryonic endocardium involves multiple factors. Notch activates certain noncoding elements through these factors and simultaneously suppresses elements that could hinder cardiac valve development and homeostasis. Biorxviv: https://www.biorxiv.org/content/10.1101/2023.03.23.533882v1.full.

The unique function of Runx1 in skeletal muscle differentiation and regeneration is mediated by an ETS interaction domain

The conserved Runt-related (RUNX) transcription factor family are well-known master regulators of developmental and regenerative processes. Runx1 and Runx2 are both expressed in satellite cells (SC) and skeletal myotubes. Conditional deletion of Runx1 in adult SC negatively impacted self-renewal and impaired skeletal muscle maintenance. Runx1- deficient SC retain Runx2 expression but cannot support muscle regeneration in response to injury. To determine the unique molecular functions of Runx1 that cannot be compensated by Runx2 we deleted Runx1 in C2C12 that retain Runx2 expression and established that myoblasts differentiation was blocked in vitro due in part to ectopic expression of Mef2c, a target repressed by Runx1 . Structure-function analysis demonstrated that the Ets-interacting MID/EID region of Runx1, absent from Runx2, is critical to regulating myoblasts proliferation, differentiation, and fusion. Analysis of in-house and published ChIP-seq datasets from Runx1 (T-cells, muscle) versus Runx2 (preosteoblasts) dependent tissue identified enrichment for a Ets:Runx composite site in Runx1 -dependent tissues. Comparing ATACseq datasets from WT and Runx1KO C2C12 cells showed that the Ets:Runx composite motif was enriched in peaks open exclusively in WT cells compared to peaks unique to Runx1KO cells. Thus, engagement of a set of targets by the RUNX1/ETS complex define the non-redundant functions of Runx1 .

Runt-related Transcription Factors and Gene Regulatory Mechanisms in Skeletal Development and Diseases

PURPOSE OF REVIEW

Runt-related transcription factors (RUNX) play critical roles in skeletal development, metabolism, and diseases. In mammals, three RUNX members, namely RUNX1, RUNX2, and RUNX3, play distinct and redundant roles, although RUNX2 is a dominant factor in skeletal development and several skeletal diseases. This review is to provide an overview of the current understanding of RUNX-mediated transcriptional regulation in different skeletal cell types.

RECENT FINDINGS

Advances in chromatin immunoprecipitation and next-generation sequencing (ChIP-seq) have revealed genome-wide RUNX-mediated gene regulatory mechanisms, including their association with cis-regulatory elements and putative target genes. Further studies with genome-wide analysis and biochemical assays have shed light on RUNX-mediated pioneering action and involvements of RUNX2 in lipid-lipid phase separation. Emerging multi-layered mechanisms of RUNX-mediated gene regulations help us better understanding of skeletal development and diseases, which also provides clues to think how genome-wide studies can help develop therapeutic strategies for skeletal diseases.

Transcription factor genetics and biology in predisposition to bone marrow failure and hematological malignancy

Transcription factors (TFs) play a critical role as key mediators of a multitude of developmental pathways, with highly regulated and tightly organized networks crucial for determining both the timing and pattern of tissue development. TFs can act as master regulators of both primitive and definitive hematopoiesis, tightly controlling the behavior of hematopoietic stem and progenitor cells (HSPCs). These networks control the functional regulation of HSPCs including self-renewal, proliferation, and differentiation dynamics, which are essential to normal hematopoiesis. Defining the key players and dynamics of these hematopoietic transcriptional networks is essential to understanding both normal hematopoiesis and how genetic aberrations in TFs and their networks can predispose to hematopoietic disease including bone marrow failure (BMF) and hematological malignancy (HM). Despite their multifaceted and complex involvement in hematological development, advances in genetic screening along with elegant multi-omics and model system studies are shedding light on how hematopoietic TFs interact and network to achieve normal cell fates and their role in disease etiology. This review focuses on TFs which predispose to BMF and HM, identifies potential novel candidate predisposing TF genes, and examines putative biological mechanisms leading to these phenotypes. A better understanding of the genetics and molecular biology of hematopoietic TFs, as well as identifying novel genes and genetic variants predisposing to BMF and HM, will accelerate the development of preventative strategies, improve clinical management and counseling, and help define targeted treatments for these diseases.

Emerging RUNX2-Mediated Gene Regulatory Mechanisms Consisting of Multi-Layered Regulatory Networks in Skeletal Development

Skeletal development is tightly coordinated by chondrocytes and osteoblasts, which are derived from skeletal progenitors, and distinct cell-type gene regulatory programs underlie the specification and differentiation of cells. Runt-related transcription factor 2 (Runx2) is essential to chondrocyte hypertrophy and osteoblast differentiation. Genetic studies have revealed the biological functions of Runx2 and its involvement in skeletal genetic diseases. Meanwhile, molecular biology has provided a framework for our understanding of RUNX2-mediated transactivation at a limited number of cis-regulatory elements. Furthermore, studies using next-generation sequencing (NGS) have provided information on RUNX2-mediated gene regulation at the genome level and novel insights into the multiple layers of gene regulatory mechanisms, including the modes of action of RUNX2, chromatin accessibility, the concept of pioneer factors and phase separation, and three-dimensional chromatin organization. In this review, I summarize the emerging RUNX2-mediated regulatory mechanism from a multi-layer perspective and discuss future perspectives for applications in the treatment of skeletal diseases.

The Landscape of Secondary Genetic Rearrangements in Pediatric Patients with B-Cell Acute Lymphoblastic Leukemia with t(12;21)

The most frequent chromosomal rearrangement in childhood B-cell acute lymphoblastic leukemia (B-ALL) is translocation t(12;21)(p13;q22). It results in the fusion of the ETV6::RUNX1 gene, which is active in the regulation of multiple crucial cellular pathways. Recent studies hypothesize that many translocations are influenced by RAG-initiated deletions, as well as defects in the RAS and NRAS pathways. According to a "two-hit" model for the molecular pathogenesis of pediatric ETV6::RUNX1-positive B-ALL, the t(12;21) translocation requires leukemia-causing secondary mutations. Patients with ETV6::RUNX1 express up to 60 different aberrations, which highlights the heterogeneity of this B-ALL subtype and is reflected in differences in patient response to treatment and chances of relapse. Most studies of secondary genetic changes have concentrated on deletions of the normal, non-rearranged ETV6 allele. Other predominant structural changes included deletions of chromosomes 6q and 9p, loss of entire chromosomes X, 8, and 13, duplications of chromosome 4q, or trisomy of chromosomes 21 and 16, but the impact of these changes on overall survival remains unclarified. An equally genetically diverse group is the recently identified new B-ALL subtype ETV6::RUNX1-like ALL. In our review, we provide a comprehensive description of recurrent secondary mutations in pediatric B-ALL with t(12;21) to emphasize the value of investigating detailed molecular mechanisms in ETV6::RUNX1-positive B-ALL, both for our understanding of the etiology of the disease and for future clinical advances in patient treatment and management.

Global DNA Methylation Analysis of Cancer-Associated Fibroblasts Reveals Extensive Epigenetic Rewiring Linked with RUNX1 Upregulation in Breast Cancer Stroma

Cancer cells recruit and rewire normal fibroblasts in their microenvironment to become protumorigenic cancer-associated fibroblasts (CAF). These CAFs are genomically stable, yet their transcriptional programs are distinct from those of their normal counterparts. Transcriptional regulation plays a major role in this reprogramming, but the extent to which epigenetic modifications of DNA also contribute to the rewiring of CAF transcription is not clear. Here we address this question by dissecting the epigenetic landscape of breast CAFs. Applying tagmentation-based whole-genome bisulfite sequencing in a mouse model of breast cancer, we found that fibroblasts undergo massive DNA methylation changes as they transition into CAFs. Transcriptional and epigenetic analyses revealed RUNX1 as a potential mediator of this process and identified a RUNX1-dependent stromal gene signature. Coculture and mouse models showed that both RUNX1 and its stromal signature are induced as normal fibroblasts transition into CAFs. In breast cancer patients, RUNX1 was upregulated in CAFs, and expression of the RUNX1 signature was associated with poor disease outcome, highlighting the relevance of these findings to human disease. This work presents a comprehensive genome-wide map of DNA methylation in CAFs and reveals a previously unknown facet of the dynamic plasticity of the stroma.

SIGNIFICANCE

The first genome-wide map of DNA methylation in breast cancer-associated fibroblasts unravels a previously unknown facet of the dynamic plasticity of the stroma, with far-reaching therapeutic implications.

Hematopoietic and Chronic Myeloid Leukemia Stem Cells: Multi-Stability versus Lineage Restriction

There is compelling evidence to support the view that the cell-of-origin for chronic myeloid leukemia is a hematopoietic stem cell. Unlike normal hematopoietic stem cells, the progeny of the leukemia stem cells are predominantly neutrophils during the disease chronic phase and there is a mild anemia. The hallmark oncogene for chronic myeloid leukemia is the BCR-ABLp210 fusion gene. Various studies have excluded a role for BCR-ABLp210 expression in maintaining the population of leukemia stem cells. Studies of BCR-ABLp210 expression in embryonal stem cells that were differentiated into hematopoietic stem cells and of the expression in transgenic mice have revealed that BCR-ABLp210 is able to veer hematopoietic stem and progenitor cells towards a myeloid fate. For the transgenic mice, global changes to the epigenetic landscape were observed. In chronic myeloid leukemia, the ability of the leukemia stem cells to choose from the many fates that are available to normal hematopoietic stem cells appears to be deregulated by BCR-ABLp210 and changes to the epigenome are also important. Even so, we still do not have a precise picture as to why neutrophils are abundantly produced in chronic myeloid leukemia.

Runx2 regulates chromatin accessibility to direct the osteoblast program at neonatal stages

The transcriptional regulator Runx2 (runt-related transcription factor 2) has essential but distinct roles in osteoblasts and chondrocytes in skeletal development. However, Runx2-mediated regulatory mechanisms underlying the distinctive programming of osteoblasts and chondrocytes are not well understood. Here, we perform an integrative analysis to investigate Runx2-DNA binding and chromatin accessibility ex vivo using neonatal osteoblasts and chondrocytes. We find that Runx2 engages with cell-type-distinct chromatin-accessible regions, potentially interacting with different combinations of transcriptional regulators, forming cell-type-specific hotspots, and potentiating chromatin accessibility. Genetic analysis and direct cellular reprogramming studies suggest that Runx2 is essential for establishment of chromatin accessibility in osteoblasts. Functional enhancer studies identify an Sp7 distal enhancer driven by Runx2-dependent binding and osteoblast-specific chromatin accessibility, contributing to normal osteoblast differentiation. Our findings provide a framework for understanding the regulatory landscape encompassing Runx2-mediated and cell-type-distinct enhancer networks that underlie the specification of osteoblasts.

Hojo et al. investigate the gene-regulatory landscape underlying specification of skeletal cell types in neonatal mice. Runx2, an osteoblast determinant, engages with cell-type-distinct chromatin-accessible regions and is essential for establishment of chromatin accessibility in osteoblasts. The study provides insights into enhancer networks in skeletal development.

RUNX regulated immune-associated genes predicts prognosis in breast cancer

Background: Breast cancer is the most common malignant tumor in women. RUNX family has been involved in the regulation of different carcinogenic processes and signaling pathways with cancer, which is closely related to immunity and prognosis of various tumors, and also plays an important role in the development and prognosis of breast cancer.

Methods: We discovered the expression of RUNX family through GEPIA Dataset and then evaluated the relationship between RUNX family and immune-related genes and the prognosis of breast cancer through analyzing TCGA database. A prognostic model was established and verified via cox proportional hazards regression model using R packages. We evaluated the accuracy of the prognostic model by Kaplan-Meier curves and receiver operating characteristic (ROC) curves. Additionally, we obtained the relationship between the RUNX family and immune infiltration by TIMER database. Finally, the dual luciferase reporter assay was used to verify the regulation of RUNX3 on potential target genes ULBP2 and TRDV1, and the effects of ULBP2 and TRDV1 on the growth of breast cancer cells were explored by CCK-8, colony formation and wound healing assays.

Results: We screened out RUNX family-regulated immune-related genes associated with the prognosis of breast cancer. These predictors included PSME2, ULBP2, IL-18, TSLP, NPR3, TRDV1. Then a prognosis-related risk score model was built using the independent risk factors to provide a clinically appropriate method predicting the overall survival (OS) probability of the patients with breast cancer. In addition, a further research was made on the functions of high risk immune gene ULBP2 and low risk immune gene TRDV1 which regulated by RUNX3, the results showed that down-regulation of ULBP2 suppressed breast cancer cell proliferation and TRDV1 had the opposite functions. The prognostic model we constructed could promote the development of prognostic, and was associated with lower immune infiltration.

Conclusion: The expression of RUNX family was closely related to the prognosis of breast cancer. At the same time, RUNX family could modulate the functions of immune-related genes, and affect the development and prognosis of breast cancer. These immune-related genes regulated by RUNX family could be promising prognostic biomarkers and therapeutic targets in breast cancer.

Plasmacytoid Dendritic Cells, a Novel Target in Myeloid Neoplasms

Plasmacytoid dendritic cells (pDC) are the main type I interferon producing cells in humans and are able to modulate innate and adaptive immune responses. Tumor infiltration by plasmacytoid dendritic cells is already well described and is associated with poor outcomes in cancers due to the tolerogenic activity of pDC. In hematological diseases, Blastic Plasmacytoid Dendritic Cells Neoplasm (BPDCN), aggressive leukemia derived from pDCs, is well described, but little is known about tumor infiltration by mature pDC described in Myeloid Neoplasms (MN). Recently, mature pDC proliferation (MPDCP) has been described as a differential diagnosis of BPDCN associated with acute myeloid leukemia (pDC-AML), myelodysplastic syndrome (pDC-MDS) and chronic myelomonocytic leukemia (pDC-CMML). Tumor cells are myeloid blasts and/or mature myeloid cells from related myeloid disorders and pDC derived from a clonal proliferation. The poor prognosis associated with MPDCP requires a better understanding of pDC biology, MN oncogenesis and immune response. This review provides a comprehensive overview about the biological aspects of pDCs, the description of pDC proliferation in MN, and an insight into putative therapies in pDC-AML regarding personalized medicine.

Leveraging transcriptome and epigenome landscapes to infer regulatory networks during the onset of sexual maturation

Background

Despite sexual development being ubiquitous to vertebrates, the molecular mechanisms underpinning this fundamental transition remain largely undocumented in many organisms. We designed a time course experiment that successfully sampled the period when Atlantic salmon commence their trajectory towards sexual maturation.

Results

Through deep RNA sequencing, we discovered key genes and pathways associated with maturation in the pituitary-ovarian axis. Analyzing DNA methylomes revealed a bias towards hypermethylation in ovary that implicated maturation-related genes. Co-analysis of DNA methylome and gene expression changes revealed chromatin remodeling genes and key transcription factors were both significantly hypermethylated and upregulated in the ovary during the onset of maturation. We also observed changes in chromatin state landscapes that were strongly correlated with fundamental remodeling of gene expression in liver. Finally, a multiomic integrated analysis revealed regulatory networks and identified hub genes including TRIM25 gene (encoding the estrogen-responsive finger protein) as a putative key regulator in the pituitary that underwent a 60-fold change in connectivity during the transition to maturation.

Conclusion

The study successfully documented transcriptome and epigenome changes that involved key genes and pathways acting in the pituitary – ovarian axis. Using a Systems Biology approach, we identified hub genes and their associated networks deemed crucial for onset of maturation. The results provide a comprehensive view of the spatiotemporal changes involved in a complex trait and opens the door to future efforts aiming to manipulate puberty in an economically important aquaculture species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08514-8.

Epigenetic and transcriptional dysregulation in CD4+ T cells in patients with atopic dermatitis

Atopic dermatitis (AD) is one of the most common skin disorders among children. Disease etiology involves genetic and environmental factors, with 29 independent AD risk loci enriched for risk allele-dependent gene expression in the skin and CD4+ T cell compartments. We investigated the potential epigenetic mechanisms responsible for the genetic susceptibility of CD4+ T cells. To understand the differences in gene regulatory activity in peripheral blood T cells in AD, we measured chromatin accessibility (an assay based on transposase-accessible chromatin sequencing, ATAC-seq), nuclear factor kappa B subunit 1 (NFKB1) binding (chromatin immunoprecipitation with sequencing, ChIP-seq), and gene expression levels (RNA-seq) in stimulated CD4+ T cells from subjects with active moderate-to-severe AD, as well as in age-matched non-allergic controls. Open chromatin regions in stimulated CD4+ T cells were highly enriched for AD genetic risk variants, with almost half of the AD risk loci overlapping AD-dependent ATAC-seq peaks. AD-specific open chromatin regions were strongly enriched for NF-κB DNA-binding motifs. ChIP-seq identified hundreds of NFKB1-occupied genomic loci that were AD- or control-specific. As expected, the AD-specific ChIP-seq peaks were strongly enriched for NF-κB DNA-binding motifs. Surprisingly, control-specific NFKB1 ChIP-seq peaks were not enriched for NFKB1 motifs, but instead contained motifs for other classes of human transcription factors, suggesting a mechanism involving altered indirect NFKB1 binding. Using DNA sequencing data, we identified 63 instances of altered genotype-dependent chromatin accessibility at 36 AD risk variant loci (30% of AD risk loci) that might lead to genotype-dependent gene expression. Based on these findings, we propose that CD4+ T cells respond to stimulation in an AD-specific manner, resulting in disease- and genotype-dependent chromatin accessibility alterations involving NFKB1 binding.

In the spotlight: the role of TGFβ signalling in haematopoietic stem and progenitor cell emergence

Haematopoietic stem and progenitor cells (HSPCs) sustain haematopoiesis by generating precise numbers of mature blood cells throughout the lifetime of an individual. In vertebrates, HSPCs arise during embryonic development from a specialised endothelial cell population, the haemogenic endothelium (HE). Signalling by the Transforming Growth Factor β (TGFβ) pathway is key to regulate haematopoiesis in the adult bone marrow, but evidence for a role in the formation of HSPCs has only recently started to emerge. In this review, we examine recent work in various model systems that demonstrate a key role for TGFβ signalling in HSPC emergence from the HE. The current evidence underpins two seemingly contradictory views of TGFβ function: as a negative regulator of HSPCs by limiting haematopoietic output from HE, and as a positive regulator, by programming the HE towards the haematopoietic fate. Understanding how to modulate the requirement for TGFβ signalling in HSC emergence may have critical implications for the generation of these cells in vitro for therapeutic use.

Pioneer transcription factors are associated with the modulation of DNA methylation patterns across cancers

Methylation of cytosines on DNA is a prominent modification associated with gene expression regulation. Aberrant DNA methylation patterns have recurrently been linked to dysregulation of the regulatory program in cancer cells. To shed light on the underlying molecular mechanism driving this process, we hypothesised that aberrant methylation patterns could be controlled by the binding of specific transcription factors (TFs) across cancer types. By combining DNA methylation arrays and gene expression data with TF binding sites (TFBSs), we explored the interplay between TF binding and DNA methylation in 19 cancer types. We performed emQTL (expression–methylation quantitative trait loci) analyses independently in each cancer type and identified 13 TFs whose expression levels are correlated with local DNA methylation patterns around their binding sites in at least 2 cancer types. The 13 TFs are mainly associated with local demethylation and are enriched for pioneer function, suggesting a specific role for these TFs in modulating chromatin structure and transcription in cancer patients. Furthermore, we confirmed that de novo methylation is precluded across cancers at CpGs lying in genomic regions enriched for TF binding signatures associated with SP1, CTCF, NRF1, GABPA, KLF9, and/or YY1. The modulation of DNA methylation associated with TF binding was observed at cis-regulatory regions controlling immune- and cancer-associated pathways, corroborating that the emQTL signals were derived from both cancer and tumor-infiltrating cells. As a case example, we experimentally confirmed that FOXA1 knock-down is associated with higher methylation in regions bound by FOXA1 in breast cancer MCF-7 cells. Finally, we reported physical interactions between FOXA1 with TET1 and TET2 both in an in vitro setup and in vivo at physiological levels in MCF-7 cells, adding further support for FOXA1 attracting TET1 and TET2 to induce local demethylation in cancer cells.