The Interplay between Chromatin and Transcription Factor Networks during B Cell Development: Who Pulls the Trigger First?

UBXN3B is crucial for B lymphopoiesis

BACKGROUND

The ubiquitin regulatory X (UBX) domain-containing proteins (UBXNs) are putative adaptors for ubiquitin ligases and valosin-containing protein; however, their in vivo physiological functions remain poorly characterised. We recently showed that UBXN3B is essential for activating innate immunity to DNA viruses and controlling DNA/RNA virus infection. Herein, we investigate its role in adaptive immunity.

METHODS

We evaluated the antibody responses to multiple viruses and pathogenesis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza in tamoxifen-inducible global and constitutive B cell-specific Ubxn3b knockout mice; quantified various immune populations, B lineage progenitors/precursors, B cell receptor (BCR) signalling and apoptosis by flow cytometry, immunoblotting and immunofluorescence microscopy. We also performed bone marrow transfer, single-cell and bulk RNA sequencing.

FINDINGS

Both global and B cell-specific Ubxn3b knockout mice present a marked reduction in small precursor B-II (>60%), immature (>70%) and mature B (>95%) cell numbers. Transfer of wildtype bone marrow to irradiated global Ubxn3b knockouts restores normal B lymphopoiesis, while reverse transplantation does not. The mature B population shrinks rapidly with apoptosis and higher pro and activated caspase-3 protein levels were observed following induction of Ubxn3b knockout. Mechanistically, Ubxn3b deficiency leads to impaired pre-BCR signalling and cell cycle arrest. Ubxn3b knockout mice are highly vulnerable to respiratory viruses, with increased viral loads and prolonged immunopathology in the lung, and reduced production of virus-specific IgM/IgG.

INTERPRETATION

UBXN3B is essential for B lymphopoiesis by maintaining constitutive pre-BCR signalling and cell survival in a cell-intrinsic manner.

FUNDING

United States National Institutes of Health grants, R01AI132526 and R21AI155820.

Signaling pathways and regulation of gene expression in hematopoietic cells

Cellular functions are regulated by signal transduction pathway networks consisting of protein-modifying enzymes that control the activity of many downstream proteins. Protein kinases and phosphatases regulate gene expression by reversible phosphorylation of transcriptional factors, which are their direct substrates. Casein kinase II (CK2) is a serine/threonine kinase that phosphorylates a large number of proteins that have critical roles in cellular proliferation, metabolism and survival. Altered function of CK2 has been associated with malignant transformation, immunological disorders and other types of diseases. Protein phosphatase 1 (PP1) is a serine/threonine phosphatase, which regulates the phosphorylation status of many proteins that are essential for cellular functions. IKAROS is a DNA-binding protein, which functions as a regulator of gene transcription in hematopoietic cells. CK2 directly phosphorylates IKAROS at multiple phosphosites which determines IKAROS activity as a regulator of gene expression. PP1 binds to IKAROS via the PP1-consensus recognition site and dephosphorylates serine/threonine residues that are phosphorylated by CK2. Thus, the interplay between CK2 and PP1 signaling pathways have opposing effects on the phosphorylation status of their mutual substrate - IKAROS. This review summarizes the effects of CK2 and PP1 on IKAROS role in regulation of gene expression and its function as a tumor suppressor in leukemia.

Genome-wide analysis of DNaseI hypersensitivity unveils open chromatin associated with histone H3 modifications after areca nut with lime exposure

Research over the years revealed that precocious anaphase, securin overexpression, and genome instability in both target and nontarget cells are significantly associated with the increased risk of areca nut (AN) and lime-induced oral, esophageal, and gastric cancers. Further, hyperphosphorylation of Rb and histone H3 epigenetic modifications both globally and in the promoter region of the securin gene were demonstrated after AN + lime exposure. This study aims whether the extract of raw AN + lime relaxes chromatin structure which further facilitates the histone H3 epigenetic modifications during the initial phase of carcinogenesis. Three groups of mice (10 in each group) were used. The treated group consumed 1 mg/day/mice of AN extract with lime ad libitum in the drinking water for 60 days. The dose was increased by 1 mg every 60 days. Isolated nuclei were digested with DNaseI and 2 kb and below DNA was eluted from the agarose gel, purified and PCR amplified by using securin and GAPDH primers. Securin and E2F1 expression, pRb phosphorylation, and histone epigenetic modifications were analyzed by immunohistochemistry. The number of DNA fragments within 2 kb in size after DNaseI treatment was higher significantly in AN + lime exposed tissue samples than in the untreated one. The PCR result showed that the number of fragments bearing securin gene promoter and GAPDH gene was significantly higher in AN + lime exposed DNaseI-treated samples. Immunohistochemistry data revealed increased Rb hyperphosphorylation, upregulation of E2F1, and securin in the AN + lime-treated samples. Increased trimethylation of histone H3 lysine 4 and acetylation of H3 lysine 9 and 18 were observed globally in the treated samples. Therefore, the results of this study have led to the hypothesis that AN + lime exposure relaxes the chromatin, changes the epigenetic landscape, and deregulates the Rb–E2F1 circuit which might be involved in the upregulation of securin and some other proto-oncogenes that might play an important role in the initial phases of AN + lime mediated carcinogenesis.

The Pleiotropy of PAX5 Gene Products and Function

PAX5, a member of the Paired Box (PAX) transcription factor family, is an essential factor for B-lineage identity during lymphoid differentiation. Mechanistically, PAX5 controls gene expression profiles, which are pivotal to cellular processes such as viability, proliferation, and differentiation. Given its crucial function in B-cell development, PAX5 aberrant expression also correlates with hallmark cancer processes leading to hematological and other types of cancer lesions. Despite the well-established association of PAX5 in the development, maintenance, and progression of cancer disease, the use of PAX5 as a cancer biomarker or therapeutic target has yet to be implemented. This may be partly due to the assortment of PAX5 expressed products, which layers the complexity of their function and role in various regulatory networks and biological processes. In this review, we provide an overview of the reported data describing PAX5 products, their regulation, and function in cellular processes, cellular biology, and neoplasm.

A network of transcription factors governs the dynamics of NODAL/Activin transcriptional responses

SMAD2, an effector of the NODAL/Activin signalling pathway, regulates developmental processes by sensing distinct chromatin states and interacting with different transcriptional partners. However, the network of factors that controls SMAD2 chromatin binding and shapes its transcriptional programme over time is poorly characterised. Here, we combine ATAC-seq with computational footprinting to identify temporal changes in chromatin accessibility and transcription factor activity upon NODAL/Activin signalling. We show that SMAD2 binding induces chromatin opening genome wide. We discover footprints for FOXI3, FOXO3 and ZIC3 at the SMAD2-bound enhancers of the early response genes, Pmepa1 and Wnt3, respectively, and demonstrate their functionality. Finally, we determine a mechanism by which NODAL/Activin signalling induces delayed gene expression, by uncovering a self-enabling transcriptional cascade whereby activated SMADs, together with ZIC3, induce the expression of Wnt3. The resultant activated WNT pathway then acts together with the NODAL/Activin pathway to regulate expression of delayed target genes in prolonged NODAL/Activin signalling conditions. This article has an associated First Person interview with the first author of the paper.

Genomic profile of radiation-induced early-onset mouse B-cell lymphoma recapitulates features of Philadelphia chromosome-like acute lymphoblastic leukemia in humans

Epidemiological studies have revealed a radiation-related increase in the risk of developing acute lymphoblastic leukemia (ALL). Our recent study revealed early induction and increased risk of precursor B-cell (pB) lymphomas in mice after radiation exposure. However, the genomic landscape of radiation-induced B-cell lymphomas remains unclear. To identify the relevant genetic alterations in mice, whole-exome sequencing was performed on both early-onset and late-onset B-cell lymphomas that developed spontaneously or after gamma-irradiation. In addition to multiple driver mutations, the data revealed that interstitial deletion of chromosome 4, including Pax5, and missense mutations in Jak3 are unique genomic alterations in radiation-induced, early-onset B-cell lymphomas. RNA sequencing revealed a pB-cell-type gene-expression profile with no involvement of known fusion genes for human ALLs in the early-onset B-cell lymphomas. Activation of Jak3/Stat5 signaling in early-onset B-cell lymphomas was validated using western capillary electrophoresis. Those features were similar to those of Philadelphia chromosome-like ALL. Our data suggest a critical role for Pax5 loss-of-function mutations in initiating B-cell leukemogenesis coupled with activation of Jak3/Stat5 signaling as a basis for the rapid development of radiation-induced pB-ALL. These molecular signatures for radiation-induced cancers will inform both risk assessment and potential targeted therapies for pB-ALL.

ATP citrate lyase controls hematopoietic stem cell fate and supports bone marrow regeneration

In order to support bone marrow regeneration after myeloablation, hematopoietic stem cells (HSCs) actively divide to provide both stem and progenitor cells. However, the mechanisms regulating HSC function and cell fate choice during hematopoietic recovery remain unclear. We herein provide novel insights into HSC regulation during regeneration by focusing on mitochondrial metabolism and ATP citrate lyase (ACLY). After 5-fluorouracil-induced myeloablation, HSCs highly expressing endothelial protein C receptor (EPCR^high ) were enriched within the stem cell fraction at the expense of more proliferative EPCR^Low HSCs. These EPCR^High HSCs were initially more primitive than EPCR^Low HSCs and enabled stem cell expansion by enhancing histone acetylation, due to increased activity of ACLY in the early phase of hematopoietic regeneration. In the late phase of recovery, HSCs enhanced differentiation potential by increasing the accessibility of cis-regulatory elements in progenitor cell-related genes, such as CD48. In conditions of reduced mitochondrial metabolism and ACLY activity, these HSCs maintained stem cell phenotypes, while ACLY-dependent histone acetylation promoted differentiation into CD48⁺ progenitor cells. Collectively, these results indicate that the dynamic control of ACLY-dependent metabolism and epigenetic alterations is essential for HSC regulation during hematopoietic regeneration.

Tcf1 and Lef1 provide constant supervision to mature CD8+ T cell identity and function by organizing genomic architecture

T cell identity is established during thymic development, but how it is maintained in the periphery remains unknown. Here we show that ablating Tcf1 and Lef1 transcription factors in mature CD8⁺ T cells aberrantly induces genes from non-T cell lineages. Using high-throughput chromosome-conformation-capture sequencing, we demonstrate that Tcf1/Lef1 are important for maintaining three-dimensional genome organization at multiple scales in CD8⁺ T cells. Comprehensive network analyses coupled with genome-wide profiling of chromatin accessibility and Tcf1 occupancy show the direct impact of Tcf1/Lef1 on the T cell genome is to promote formation of extensively interconnected hubs through enforcing chromatin interaction and accessibility. The integrative mechanisms utilized by Tcf1/Lef1 underlie activation of T cell identity genes and repression of non-T lineage genes, conferring fine control of various T cell functionalities. These findings suggest that Tcf1/Lef1 control global genome organization and help form intricate chromatin-interacting hubs to facilitate promoter-enhancer/silencer contact, hence providing constant supervision of CD8⁺ T cell identity and function.

How CD8⁺ T cell identity is maintained after exit from the thymus is not fully established. Here the authors use multiomics approaches including Hi-C to show that Tcf1 and Lef1 prevent aberrant expression of lineage-inappropriate genes by organizing three-dimensional genomic architecture in CD8⁺ T cells.

Chromatin Remodeling of Colorectal Cancer Liver Metastasis is Mediated by an HGF-PU.1-DPP4 Axis

Colorectal cancer (CRC) metastasizes mainly to the liver, which accounts for the majority of CRC-related deaths. Here it is shown that metastatic cells undergo specific chromatin remodeling in the liver. Hepatic growth factor (HGF) induces phosphorylation of PU.1, a pioneer factor, which in turn binds and opens chromatin regions of downstream effector genes. PU.1 increases histone acetylation at the DPP4 locus. Precise epigenetic silencing by CRISPR/dCas9KRAB or CRISPR/dCas9HDAC revealed that individual PU.1-remodeled regulatory elements collectively modulate DPP4 expression and liver metastasis growth. Genetic silencing or pharmacological inhibition of each factor along this chromatin remodeling axis strongly suppressed liver metastasis. Therefore, microenvironment-induced epimutation is an important mechanism for metastatic tumor cells to grow in their new niche. This study presents a potential strategy to target chromatin remodeling in metastatic cancer and the promise of repurposing drugs to treat metastasis.

Activated B Cells and Plasma Cells Are Resistant to Radiation Therapy

PURPOSE

B cells play a key role in outcomes of cancer patients and responses to checkpoint blockade immunotherapies. However, the effect of radiation therapy on B cell populations is poorly understood. Here we characterize the effects of radiation on the development, survival, and phenotype of physiological B-cell subsets.

METHODS AND MATERIALS

Naïve and immunized tumor bearing and nontumor bearing mice were treated with large-field or focal stereotactic radiation and distinct B-cell subsets of varying developmental stages were analyzed by flow cytometry and real-time reverse transcription polymerase chain reaction.

RESULTS

We first report that focal stereotactic radiation is highly superior to large-field radiation at inducing tumor infiltration of B cells, CD8⁺ T cells, and macrophages. We observed that radiation affects B cell development in the bone marrow, increasing frequencies of early pro-B cells and late pro-B cells while inducing upregulation of programmed cell death protein 1. We then demonstrate that class switched B cells and plasma cells are highly resistant to radiation therapy compared with naïve B cells and upregulate activation markers programmed cell death 1 ligand 2 and major histocompatibility complex class II) after radiation. Mechanistically, radiation upregulates Xbp1 and Bcl6 in plasma cells, conferring radioresistance. Furthermore, using an immunization approach, we demonstrate that radiation enhances activation-induced cytidine deaminase mediated class switching and somatic hypermutation in primed B cells.

CONCLUSIONS

These data demonstrate that stereotactic radiation is superior to large-field radiation at inducing infiltration of immune cells into tumors and that plasma cells and class switched B cells are highly resistant to radiation therapy. These results represent the most comprehensive analysis of the effects of radiation on B cells to date and identify novel mechanisms by which radiation modulates immune cells within the tumor microenvironment.

An Essential Role of UBXN3B in B Lymphopoiesis

Hematopoiesis is finely regulated to enable timely production of the right numbers and types of mature immune cells to maintain tissue homeostasis. Dysregulated hematopoiesis may compromise antiviral immunity and/or exacerbate immunopathogenesis. Herein, we report an essential role of UBXN3B in maintenance of hematopoietic homeostasis and restriction of immunopathogenesis during respiratory viral infection. Ubxn3b deficient ( Ubxn3b ^-/- ) mice are highly vulnerable to SARS-CoV-2 and influenza A infection, characterized by more severe lung immunopathology, lower virus-specific IgG, significantly fewer B cells, but more myeloid cells than Ubxn3b ^+/+ littermates. This aberrant immune compartmentalization is recapitulated in uninfected Ubxn3b ^-/- mice. Mechanistically, UBXN3B controls precursor B-I (pre-BI) transition to pre-BII and subsequent proliferation in a cell-intrinsic manner, by maintaining BLNK protein stability and pre-BCR signaling. These results reveal an essential role of UBXN3B for the early stage of B cell development.

Transcriptional Reprogramming and Constitutive PD-L1 Expression in Melanoma Are Associated with Dedifferentiation and Activation of Interferon and Tumour Necrosis Factor Signalling Pathways

Simple Summary

Melanoma, an aggressive form of skin cancer, is frequently associated with drug resistance in the advanced stages. For instance, frequently resistance is observed in sequential treatment of melanoma with targeted therapy and immunotherapy. In this research, the authors investigated whether potential transcriptional mechanisms and pathways associated with PD-L1 protein expression could underlie targeted therapy drug resistance in melanoma. The authors found a PD-L1 expression transcriptional pattern underlies resistance to targeted therapy in a subgroup of melanomas. These melanomas were markedly dedifferentiated, as compared to melanomas that were not drug resistant. Understanding changes in transcription and molecular pathways that lead to drug resistance could allow researchers to develop interventions to prevent drug resistance from occurring in melanoma, which could also be relevant to other cancer types.

Abstract

Melanoma is the most aggressive type of skin cancer, with increasing incidence worldwide. Advances in targeted therapy and immunotherapy have improved the survival of melanoma patients experiencing recurrent disease, but unfortunately treatment resistance frequently reduces patient survival. Resistance to targeted therapy is associated with transcriptomic changes and has also been shown to be accompanied by increased expression of programmed death ligand 1 (PD-L1), a potent inhibitor of immune response. Intrinsic upregulation of PD-L1 is associated with genome-wide DNA hypomethylation and widespread alterations in gene expression in melanoma cell lines. However, an in-depth analysis of the transcriptomic landscape of melanoma cells with intrinsically upregulated PD-L1 expression is lacking. To determine the transcriptomic landscape of intrinsically upregulated PD-L1 expression in melanoma, we investigated transcriptomes in melanomas with constitutive versus inducible PD-L1 expression (referred to as PD-L1_CON and PD-L1_IND). RNA-Seq analysis was performed on seven PD-L1_CON melanoma cell lines and ten melanoma cell lines with low inducible PD-L1_IND expression. We observed that PD-L1_CON melanoma cells had a reprogrammed transcriptome with a characteristic pattern of dedifferentiated gene expression, together with active interferon (IFN) and tumour necrosis factor (TNF) signalling pathways. Furthermore, we identified key transcription factors that were also differentially expressed in PD-L1_CON versus PD-L1_IND melanoma cell lines. Overall, our studies describe transcriptomic reprogramming of melanomas with PD-L1_CON expression.

Using high-throughput multi-omics data to investigate structural balance in elementary gene regulatory network motifs

MOTIVATION

The simultaneous availability of ATAC-seq and RNA-seq experiments allows to obtain a more in-depth knowledge on the regulatory mechanisms occurring in gene regulatory networks. In this article, we highlight and analyze two novel aspects that leverage on the possibility of pairing RNA-seq and ATAC-seq data. Namely we investigate the causality of the relationships between transcription factors, chromatin and target genes and the internal consistency between the two omics, here measured in terms of structural balance in the sample correlations along elementary length-3 cycles.

RESULTS

We propose a framework that uses the a priori knowledge on the data to infer elementary causal regulatory motifs (namely chains and forks) in the network. It is based on the notions of conditional independence and partial correlation, and can be applied to both longitudinal and non-longitudinal data. Our analysis highlights a strong connection between the causal regulatory motifs that are selected by the data and the structural balance of the underlying sample correlation graphs: strikingly, >97% of the selected regulatory motifs belong to a balanced subgraph. This result shows that internal consistency, as measured by structural balance, is close to a necessary condition for 3-node regulatory motifs to satisfy causality rules.

AVAILABILITY AND IMPLEMENTATION

The analysis was carried out in MATLAB and the code can be found at https://github.com/albertozenere/Multi-omics-elementary-regulatory-motifs.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Enhancing B-Cell Malignancies-On Repurposing Enhancer Activity towards Cancer

B-cell lymphomas and leukemias derive from B cells at various stages of maturation and are the 6th most common cancer-related cause of death. While the role of several oncogenes and tumor suppressors in the pathogenesis of B-cell neoplasms was established, recent research indicated the involvement of non-coding, regulatory sequences. Enhancers are DNA elements controlling gene expression in a cell type- and developmental stage-specific manner. They ensure proper differentiation and maturation of B cells, resulting in production of high affinity antibodies. However, the activity of enhancers can be redirected, setting B cells on the path towards cancer. In this review we discuss different mechanisms through which enhancers are exploited in malignant B cells, from the well-studied translocations juxtaposing oncogenes to immunoglobulin loci, through enhancer dysregulation by sequence variants and mutations, to enhancer hijacking by viruses. We also highlight the potential of therapeutic targeting of enhancers as a direction for future investigation.

Role of Tet2 in Regulating Adaptive and Innate Immunity

Accumulated evidence indicates that epigenetic modifications play central roles in gene expression regulation and participate in developing many autoimmune and autoinflammatory diseases. Mechanistically, epigenetic modifications act as a bridge between environmental and cellular factors and susceptibility genes. DNA methylation is a critical epigenetic modification that is regulated by ten-eleven translocation (TET) enzymes. Accumulating evidence has revealed that TET family proteins function as gene regulators and antitumor drug targets mainly because of their ability to oxidize 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC). Recently, the effect of Tet2, an essential TET protein, on the development of autoimmune diseases has been explored. In this review, we summarize the current understanding of Tet2 in immune response regulation, clarify the mechanisms of Tet2 in B and T cell differentiation and function, and discuss the opposing effects of Tet2 on inflammatory gene expression in the immune system to provide new potential therapeutic targets for related diseases.

Identification of novel, clonally stable, somatic mutations targeting transcription factors PAX5 and NKX2-3, the epigenetic regulator LRIF1, and BRAF in a case of atypical B-cell chronic lymphocytic leukemia harboring a t(14;18)(q32;q21)

Diagnosis of B-cell chronic lymphocytic leukemia (B-CLL) is usually straightforward, involving clinical, immunophenotypic (Matutes score), and (immuno)genetic analyses (to refine patient prognosis for treatment). CLL cases with atypical presentation (e.g., Matutes ≤ 3) are also encountered, and for these diseases, biology and prognostic impact are less clear. Here we report the genomic characterization of a case of atypical B-CLL in a 70-yr-old male patient; B-CLL cells showed a Matutes score of 3, chromosomal translocation t(14;18)(q32;q21) (BCL2/IGH), mutated IGHV, deletion 17p, and mutations in BCL2, NOTCH1 (subclonal), and TP53 (subclonal). Quite strikingly, a novel PAX5 mutation that was predicted to be loss of function was also seen. Exome sequencing identified, in addition, a potentially actionable BRAF mutation, together with novel somatic mutations affecting the homeobox transcription factor NKX2-3, known to control B-lymphocyte development and homing, and the epigenetic regulator LRIF1, which is implicated in chromatin compaction and gene silencing. Neither NKX2-3 nor LRIF1 mutations, predicted to be loss of function, have previously been reported in B-CLL. Sequencing confirmed the presence of these mutations together with BCL2, NOTCH1, and BRAF mutations, with the t(14;18)(q32;q21) translocation, in the initial diagnostic sample obtained 12 yr prior. This is suggestive of a role for these novel mutations in B-CLL initiation and stable clonal evolution, including upon treatment withdrawal. This case extends the spectrum of atypical B-CLL with t(14;18)(q32;q21) and highlights the value of more global precision genomics for patient follow-up and treatment in these patients.

Tet2 and Tet3 in B cells are required to repress CD86 and prevent autoimmunity

A contribution of epigenetic modifications to B cell tolerance has been proposed but not directly tested. Here we report that deficiency of ten-eleven translocation (Tet) DNA demethylase family members Tet2 and Tet3 in B cells led to hyperactivation of B and T cells, autoantibody production and lupus-like disease in mice. Mechanistically, in the absence of Tet2 and Tet3, downregulation of CD86, which normally occurs following chronic exposure of self-reactive B cells to self-antigen, did not take place. The importance of dysregulated CD86 expression in Tet2- and Tet3-deficient B cells was further demonstrated by the restriction, albeit not complete, on aberrant T and B cell activation following anti-CD86 blockade. Tet2- and Tet3-deficient B cells had decreased accumulation of histone deacetylase 1 (HDAC1) and HDAC2 at the Cd86 locus. Thus, our findings suggest that Tet2- and Tet3-mediated chromatin modification participates in repression of CD86 on chronically stimulated self-reactive B cells, which contributes, at least in part, to preventing autoimmunity.

Mechanisms governing the pioneering and redistribution capabilities of the non-classical pioneer PU.1

Establishing gene regulatory networks during differentiation or reprogramming requires master or pioneer transcription factors (TFs) such as PU.1, a prototype master TF of hematopoietic lineage differentiation. To systematically determine molecular features that control its activity, here we analyze DNA-binding in vitro and genome-wide in vivo across different cell types with native or ectopic PU.1 expression. Although PU.1, in contrast to classical pioneer factors, is unable to access nucleosomal target sites in vitro, ectopic induction of PU.1 leads to the extensive remodeling of chromatin and redistribution of partner TFs. De novo chromatin access, stable binding, and redistribution of partner TFs both require PU.1's N-terminal acidic activation domain and its ability to recruit SWI/SNF remodeling complexes, suggesting that the latter may collect and distribute co-associated TFs in conjunction with the non-classical pioneer TF PU.1.

Challenges and Opportunities for Consistent Classification of Human B Cell and Plasma Cell Populations

The increasingly recognized role of different types of B cells and plasma cells in protective and pathogenic immune responses combined with technological advances have generated a plethora of information regarding the heterogeneity of this human immune compartment. Unfortunately, the lack of a consistent classification of human B cells also creates significant imprecision on the adjudication of different phenotypes to well-defined populations. Additional confusion in the field stems from: the use of non-discriminatory, overlapping markers to define some populations, the extrapolation of mouse concepts to humans, and the assignation of functional significance to populations often defined by insufficient surface markers. In this review, we shall discuss the current understanding of human B cell heterogeneity and define major parental populations and associated subsets while discussing their functional significance. We shall also identify current challenges and opportunities. It stands to reason that a unified approach will not only permit comparison of separate studies but also improve our ability to define deviations from normative values and to create a clean framework for the identification, functional significance, and disease association with new populations.

Cellular signaling and epigenetic regulation of gene expression in leukemia

Alterations in normal regulation of gene expression is one of the key features of hematopoietic malignancies. In order to gain insight into the mechanisms that regulate gene expression in these diseases, we dissected the role of the Ikaros protein in leukemia. Ikaros is a DNA-binding, zinc finger protein that functions as a transcriptional regulator and a tumor suppressor in leukemia. The use of ChIP-seq, RNA-seq, and ATAC-seq—coupled with functional experiments—revealed that Ikaros regulates both the global epigenomic landscape and epigenetic signature at promoter regions of its target genes. Casein kinase II (CK2), an oncogenic kinase that is overexpressed in leukemia, directly phosphorylates Ikaros at multiple, evolutionarily-conserved residues. Phosphorylation of Ikaros impairs the protein's ability to regulate both the transcription of its target genes and global epigenetic landscape in leukemia. Treatment of leukemia cells with a specific inhibitor of CK2 restores Ikaros function, resulting in cytotoxicity of leukemia cells. Here, we review the mechanisms through which the CK2-Ikaros signaling axis regulates the global epigenomic landscape and expression of genes that control cellular proliferation in leukemia.

A framework for identifying dysregulated chromatin regulators as master regulators in human cancer

MOTIVATION

Chromatin regulators (CRs) are frequently dysregulated to reprogram the epigenetic landscape of the cancer genome. However, the underpinnings of the dysregulation of CRs and their downstream effectors remain to be elucidated.

RESULTS

Here, we designed an integrated framework based on multi-omics data to identify candidate master regulatory CRs affected by genomic alterations across eight cancer types in The Cancer Genome Atlas. Most of them showed consistent activated or repressed (i.e. oncogenic or tumor-suppressive) roles in cancer initiation and progression. In order to further explore the insight mechanism of the dysregulated CRs, we developed an R package ModReg based on differential connectivity to identify CRs as modulators of transcription factors (TFs) involved in tumorigenesis. Our analysis revealed that the connectivity between TFs and their target genes (TGs) tended to be disrupted in the patients who had a high expression of oncogenic CRs or low-expression of tumor-suppressive CRs. As a proof-of-principle study, 14 (82.4%) of the top-ranked 17 driver CRs in liver cancer were able to be validated by literature mining or experiments including shRNA knockdown and dCas9-based epigenetic editing. Moreover, we confirmed that CR SIRT7 physically interacted with TF NFE2L2, and positively modulated the transcriptional program of NFE2L2 by affecting ∼64% of its TGs.

AVAILABILITY AND IMPLEMENTATION

ModReg is freely accessible at http://cis.hku.hk/software/ModReg.tar.gz.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Ikaros tumor suppressor function includes induction of active enhancers and super-enhancers along with pioneering activity

Ikaros encodes a transcription factor that functions as a tumor suppressor in T-cell acute lymphoblastic leukemia (T-ALL). The mechanisms through which Ikaros regulates gene expression and cellular proliferation in T-ALL are unknown. Re-introduction of Ikaros into Ikaros-null T-ALL cells resulted in cessation of cellular proliferation and induction of T-cell differentiation. We performed dynamic, global, epigenomic and gene expression analyses to determine the mechanisms of Ikaros tumor suppressor activity. Our results identified novel Ikaros functions in the epigenetic regulation of gene expression: Ikaros directly regulates de novo formation and depletion of enhancers, de novo formation of active enhancers and activation of poised enhancers; Ikaros directly induces the formation of super-enhancers; and Ikaros demonstrates pioneering activity by directly regulating chromatin accessibility. Dynamic analyses demonstrate the long-lasting effects of Ikaros DNA binding on enhancer activation, de novo formation of enhancers and super-enhancers, and chromatin accessibility. Our results establish that Ikaros’ tumor suppressor function occurs via global regulation of the enhancer and super-enhancer landscape and through pioneering activity. Expression analysis identified a large number of novel signaling pathways that are directly regulated by Ikaros and Ikaros-induced enhancers, and that are responsible for the cessation of proliferation and induction of T-cell differentiation in T-ALL cells.

The E3 ubiquitin ligase Itch is required for B-cell development

The E3 ubiquitin ligase Itch interacts with Foxo1 and targets it for ubiquitination and degradation during follicular helper T-cell differentiation, whereas the transcription factor Foxo1 plays a critical role in B-cell development. Thus, we proposed that Itch mediates B-cell differentiation. Unexpectedly, we found that Itch deficiency downregulated Foxo1 expression in B cells. Itch cKO (conditional knock out in B cells) mice had fewer pro-B cells in the bone marrow, more small resting IgM⁻IgD⁻B cells in the periphery, and lower B-cell numbers in the lymph nodes through decreased Foxo1-mediated IL-7Rα, RAG, and CD62L expression, respectively. Importantly, Itch deficiency reduced Foxo1 mRNA expression by up-regulating JunB-mediated miR-182. Finally, Foxo1 negatively regulated JunB expression by up-regulating Itch. Thus, we have identified a novel regulatory axis between Itch and Foxo1 in B cells, suggesting that Itch is essential for B-cell development.

BRWD1 orchestrates epigenetic landscape of late B lymphopoiesis

Transcription factor (TF) networks determine cell fate in hematopoiesis. However, how TFs cooperate with other regulatory mechanisms to instruct transcription remains poorly understood. Here we show that in small pre-B cells, the lineage restricted epigenetic reader BRWD1 closes early development enhancers and opens the enhancers of late B lymphopoiesis to TF binding. BRWD1 regulates over 7000 genes to repress proliferative and induce differentiation programs. However, BRWD1 does not regulate the expression of TFs required for B lymphopoiesis. Hypogammaglobulinemia patients with BRWD1 mutations have B-cell transcriptional profiles and enhancer landscapes similar to those observed in Brwd1^-/- mice. These data indicate that, in both mice and humans, BRWD1 is a master orchestrator of enhancer accessibility that cooperates with TF networks to drive late B-cell development.

B-cell development is tightly regulated by transcription programs that are coordinated by transcription factors (TF) and locus accessibility. Here the authors show that, in mice and humans, the epigenetic reader BRWD1 inhibits and promotes the accessibility of enhancers for early and late B lymphopoiesis, respectively.

Comprehensive profiling of transcriptional networks specific for lactogenic differentiation of HC11 mammary epithelial stem-like cells

The development of mammary gland as a lactogenic tissue is a highly coordinated multistep process. The epithelial cells of lactiferous tubules undergo profound changes during the developmental window of puberty, pregnancy, and lactation. Several hormones including estrogen, progesterone, glucocorticoids and prolactin act in concert, and orchestrate the development of mammary gland. Understanding the gene regulatory networks that coordinate proliferation and differentiation of HC11 Mammary Epithelial stem-like Cells (MEC) under the influence of lactogenic hormones is critical for elucidating the mechanism of lactogenesis in detail. In this study, we analyzed transcriptome profiles of undifferentiated MEC (normal) and compared them with Murine Embryonic Stem Cells (ESC) using next-generation mRNA sequencing. Further, we analyzed the transcriptome output during lactogenic differentiation of MEC following treatment with glucocorticoids (primed state) and both glucocorticoids and prolactin together (prolactin state). We established stage-specific gene regulatory networks in ESC and MEC (normal, priming and prolactin states). We validated the top up-and downregulated genes in each stage of differentiation of MEC by RT-PCR and found that they are comparable with that of RNA-seq data. HC11 MEC display decreased expression of Pou5f1 and Sox2, which is crucial for the differentiation of MEC, which otherwise ensure pluripotency to ESC. Cited4 is induced during priming and is involved in milk secretion. MEC upon exposure to both glucocorticoids and prolactin undergo terminal differentiation, which is associated with the expression of several genes, including Xbp1 and Cbp that are required for cell growth and differentiation. Our study also identified differential expression of transcription factors and epigenetic regulators in each stage of lactogenic differentiation. We also analyzed the transcriptome data for the pathways that are selectively activated during lactogenic differentiation. Further, we found that selective expression of chromatin modulators (Dnmt3l, Chd9) in response to glucocorticoids suggests a highly coordinated stage-specific lactogenic differentiation of MEC.

Regulation of the Development and Function of B Cells by ZBTB Transcription Factors

The large ZBTB family comprises a diverse group of transcriptional factors. Several ZBTB proteins have emerged as critical factors that regulate the lineage commitment, differentiation, and function of lymphoid cells as well as many other developmental events. For instance, dysfunctions of ZBTB20 or ZBTB24 have been linked to multisystem failures in humans. Within the B-cell lineage, BCL6, ZBTB7A, ZBTB17, and ZBTB1 regulate the development/differentiation of B cells in both bone marrow and peripheral lymphoid organs, while ZBTB20 and ZBTB32 seem to mainly impact the maintenance of terminal plasma cells. Given the importance of B cells in the prevention and treatment of infectious or autoimmune disorders, we herein summarize the roles of seven ZBTB family members (BCL6, ZBTB7A, ZBTB17, ZBTB20, ZBTB32, ZBTB1, and ZBTB24) in the development, differentiation, and function of B cells as well as the underlying molecular mechanisms.

Charting the dynamic epigenome during B-cell development

The epigenetic landscape undergoes a widespread modulation during embryonic development and cell differentiation. Within the hematopoietic system, B cells are perhaps the cell lineage with a more dynamic DNA methylome during their maturation process, which involves approximately one third of all the CpG sites of the genome. Although each B-cell maturation step displays its own DNA methylation fingerprint, the DNA methylome is more extensively modified in particular maturation transitions. These changes are gradually accumulated in specific chromatin environments as cell differentiation progresses and reflect different features and functional states of B cells. Promoters and enhancers of B-cell transcription factors acquire activation-related epigenetic marks and are sequentially expressed in particular maturation windows. These transcription factors further reconfigure the epigenetic marks and activity state of their target sites to regulate the expression of genes related to B-cell functions. Together with this observation, extensive DNA methylation changes in areas outside gene regulatory elements such as hypomethylation of heterochromatic regions and hypermethylation of CpG-rich regions, also take place in mature B cells, which intriguingly have been described as hallmarks of cancer. This process starts in germinal center B cells, a highly proliferative cell type, and becomes particularly apparent in long-lived cells such as memory and plasma cells. Overall, the characterization of the DNA methylome during B-cell differentiation not only provides insights into the complex epigenetic network of regulatory elements that mediate the maturation process but also suggests that late B cells also passively accumulate epigenetic changes related to cell proliferation and longevity.

PBX1 as Pioneer Factor: A Case Still Open

Pioneer factors are proteins that can recognize their target sites in barely accessible chromatin and initiate a cascade of events that allows for later transcriptional activation of the respective genes. Pioneer factors are therefore particularly well-suited to initiate cell fate changes. To date, only a small number of pioneer factors have been identified and studied in depth, such as FOXD3/FOXA1, OCT4, or SOX2. Interestingly, several recent studies reported that the PBC transcription factor PBX1 can access transcriptionally inactive genomic loci. Here, we summarize the evidence linking PBX1 with transcriptional pioneer functions, suggest potential mechanisms involved and discuss open questions to be resolved.

Feature-based classification of human transcription factors into hypothetical sub-classes related to regulatory function

Background

Transcription factors are key proteins in the regulation of gene transcription. An important step in this process is the opening of chromatin in order to make genomic regions available for transcription. Data on DNase I hypersensitivity has previously been used to label a subset of transcription factors as Pioneers, Settlers and Migrants to describe their potential role in this process. These labels represent an interesting hypothesis on gene regulation and possibly a useful approach for data analysis, and therefore we wanted to expand the set of labeled transcription factors to include as many known factors as possible. We have used a well-annotated dataset of 1175 transcription factors as input to supervised machine learning methods, using the subset with previously assigned labels as training set. We then used the final classifier to label the additional transcription factors according to their potential role as Pioneers, Settlers and Migrants. The full set of labeled transcription factors was used to investigate associated properties and functions of each class, including an analysis of interaction data for transcription factors based on DNA co-binding and protein-protein interactions. We also used the assigned labels to analyze a previously published set of gene lists associated with a time course experiment on cell differentiation.

Results

The analysis showed that the classification of transcription factors with respect to their potential role in chromatin opening largely was determined by how they bind to DNA. Each subclass of transcription factors was enriched for properties that seemed to characterize the subclass relative to its role in gene regulation, with very general functions for Pioneers, whereas Migrants to a larger extent were associated with specific processes. Further analysis showed that the expanded classification is a useful resource for analyzing other datasets on transcription factors with respect to their potential role in gene regulation. The analysis of transcription factor interaction data showed complementary differences between the subclasses, where transcription factors labeled as Pioneers often interact with other transcription factors through DNA co-binding, whereas Migrants to a larger extent use protein-protein interactions. The analysis of time course data on cell differentiation indicated a shift in the regulatory program associated with Pioneer-like transcription factors during differentiation.

Conclusions

The expanded classification is an interesting resource for analyzing data on gene regulation, as illustrated here on transcription factor interaction data and data from a time course experiment. The potential regulatory function of transcription factors seems largely to be determined by how they bind DNA, but is also influenced by how they interact with each other through cooperativity and protein-protein interactions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1349-2) contains supplementary material, which is available to authorized users.

Distinct Genetic Networks Orchestrate the Emergence of Specific Waves of Fetal and Adult B-1 and B-2 Development

B cell development is often depicted as a linear process initiating in the fetus and continuing postnatally. Using a PU.1 hypomorphic mouse model, we found that B-1 and B-2 lymphopoiesis occurred in distinct fetal and adult waves differentially dependent on the Sfpi1 14 kB upstream regulatory element. The initial wave of fetal B-1 development was absent in PU.1 hypomorphic mice, while subsequent fetal and adult waves emerged. In contrast, B-2 lymphopoiesis occurred in distinct fetal and adult waves. Whole-transcriptome profiling of fetal and adult B cell progenitors supported the existence of three waves of B-1 and two waves of B-2 development and revealed that the network of transcription factors governing B lineage specification and commitment was highly divergent between B-1 and B-2 progenitors. These findings support the view that the B-1 and B-2 lineages are distinct and provide a genetic basis for layering of immune system development.

Zfp521 promotes B-cell viability and cyclin D1 gene expression in a B cell culture system

•

Knockdown of Zfp521 in BCL1 cell culture reduces viability and promotes apoptosis.

•

Genes expressed in B cells are down-regulated in cells with Zfp521 knockdown.

•

Cyclin D1 expression is increased in mouse tumors with Zfp521 over-expression.

Leukemia arises due to the dysregulated proliferation of hematopoietic progenitor cells. Errors in the multi-step commitment process result in excessive numbers of immature lymphocytes, causing malignant disease. Genes involved in the differentiation of lymphocytes are often associated with leukemia. One such gene, Zfp521, has been found to cause B-cell leukemia in mice when over-expressed. The role of Zfp521 in B-cell differentiation, and the mechanisms by which it leads to leukemic transformation, are unclear. In this study we report that Zfp521 knockdown causes apoptosis in a B-cell culture system and promotes down-regulation of genes acting at late stages of B-cell differentiation. We identify Pax5 and cyclin D1 as Zfp521 target genes, and suggest that excessive B-cell proliferation observed in mice with retroviral insertions near the Zfp521 gene is due to an up-regulation of cyclin D1 in B-cells. Overall, these results suggest links between dysregulated Zfp521 and B-cell survival.

ZNF423 and ZNF521: EBF1 Antagonists of Potential Relevance in B-Lymphoid Malignancies

The development of the B-lymphoid cell lineage is tightly controlled by the concerted action of a network of transcriptional and epigenetic regulators. EBF1, a central component of this network, is essential for B-lymphoid specification and commitment as well as for the maintenance of the B-cell identity. Genetic alterations causing loss of function of these B-lymphopoiesis regulators have been implicated in the pathogenesis of B-lymphoid malignancies, with particular regard to B-cell acute lymphoblastic leukaemias (B-ALLs), where their presence is frequently detected. The activity of the B-cell regulatory network may also be disrupted by the aberrant expression of inhibitory molecules. In particular, two multi-zinc finger transcription cofactors named ZNF423 and ZNF521 have been characterised as potent inhibitors of EBF1 and are emerging as potentially relevant contributors to the development of B-cell leukaemias. Here we will briefly review the current knowledge of these factors and discuss the importance of their functional cross talk with EBF1 in the development of B-cell malignancies.

Enhancer repertoires are reshaped independently of early priming and heterochromatin dynamics during B cell differentiation

A widely accepted model posits that activation of enhancers during differentiation goes through a priming step prior to lineage commitment. To investigate the chronology of enhancer repertoire establishment during hematopoiesis, we monitored epigenome dynamics during three developmental stages representing hematopoietic stem cells, B-cell progenitors and mature B-cells. We find that only a minority of enhancers primed in stem cells or progenitors become active at later stages. Furthermore, most enhancers active in differentiated cells were not primed in earlier stages. Thus, the enhancer repertoire is reshaped dynamically during B-cell differentiation and enhancer priming in early stages does not appear to be an obligate step for enhancer activation. Furthermore, our data reveal that heterochromatin and Polycomb-mediated silencing have only a minor contribution in shaping enhancer repertoires during cell differentiation. Together, our data revisit the prevalent model about epigenetic reprogramming during hematopoiesis and give insights into the formation of gene regulatory networks.

Enhancers in differentiated haematopoietic cells are generally believed to be primed prior to lineage commitment. Here, the authors show that early priming and Polycomb group mediated silencing have minor roles in shaping the enhancer repertoire in differentiated B cells and that most active enhancers are generated de novo.