EBF1 and PAX5 control pro-B cell expansion via opposing regulation of the Myc gene

DNA methylation regulates B cell activation via repressing Pax5 expression in teleost

In mammals, the transcription factor Pax5 is a key regulator of B cell development and maturation and specifically expressed in naive/mature B cells but repressed upon B cell activation. Despite the long-standing proposal that Pax5 repression is essential for proper B cell activation, the underlying mechanisms remain largely elusive. In this study, we used a teleost model to elucidate the mechanisms governing Pax5 repression during B cell activation. Treatment with lipopolysaccharide (LPS) and chitosan oligosaccharide (COS) significantly enhanced the antibody secreting ability and phagocytic capacity of IgM+ B cells in large yellow croaker (Larimichthys crocea), coinciding with upregulated expression of activation-related genes, such as Bcl6, Blimp1, and sIgM, and downregulated expression of Pax5. Intriguingly, two CpG islands were identified within the promoter region of Pax5. Both CpG islands exhibited hypomethylation in naive/mature B cells, while CpG island1 was specifically transited into hypermethylation upon B cell activation. Furthermore, treatment with DNA methylation inhibitor 5-aza-2'-deoxycytidine (AZA) prevented the hypermethylation of CpG island1, and concomitantly impaired the downregulation of Pax5 and activation of B cells. Finally, through in vitro methylation experiments, we demonstrated that DNA methylation exerts an inhibitory effect on promoter activities of Pax5. Taken together, our findings unveil a novel mechanism underlying Pax5 repression during B cell activation, thus promoting the understanding of B cell activation process.

EBF1, PAX5, and MYC: regulation on B cell development and association with hematologic neoplasms

During lymphocyte development, a diverse repertoire of lymphocyte antigen receptors is produced to battle against pathogens, which is the basis of adaptive immunity. The diversity of the lymphocyte antigen receptors arises primarily from recombination-activated gene (RAG) protein-mediated V(D)J rearrangement in early lymphocytes. Furthermore, transcription factors (TFs), such as early B cell factor 1 (EBF1), paired box gene 5 (PAX5), and proto-oncogene myelocytomatosis oncogene (MYC), play critical roles in regulating recombination and maintaining normal B cell development. Therefore, the aberrant expression of these TFs may lead to hematologic neoplasms.

Effect of BAFF blockade on the B cell receptor repertoire and transcriptome in a mouse model of systemic lupus erythematosus

Introduction

Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease. Anti-B-cell-activating factor (BAFF) therapy effectively depletes B cells and reduces SLE disease activity. This research aimed to evaluate the effect of BAFF blockade on B cell receptor (BCR) repertoire and gene expression.

Methods

Through next-generation sequencing, we analyzed gene expression and BCR repertoire in MRL/lpr mice that received long-term anti-BAFF therapy. Based on gene expression profiles, we predicted the relative proportion of immune cells using ImmuCellAI-mouse, validating our predictions via flow cytometry and FluoroSpot.

Results

The loss of BCR repertoire diversity and richness, along with increased clonality and differential frequency distribution of the immunoglobulin heavy chain variable (IGHV) segment gene usage, were observed in BAFF-blockade mice. Meanwhile, the distribution of complementarity-determining region 3 (CDR3) length and CDR3 amino acid usage remained unaffected. BAFF blockade resulted in extensive changes in gene expression, particularly that of genes related to B cells and immunoglobulins. Besides, the tumor necrosis factor (TNF)-α responses and interferon (IFN)-α/γ were downregulated, consistent with the decrease in IFN-γ and TNF-α serum levels following anti-BAFF therapy. In addition, BAFF blockade significantly reduced B cell subpopulations and plasmacytoid dendritic cells, and caused the depletion of antibody-secreting cells.

Discussion

Our comparative BCR repertoire and transcriptome analyses of MRL/lpr mice subjected to BAFF blockade provide innovative insights into the molecular pathophysiology of SLE.

Intra- and interchromosomal contact mapping reveals the Igh locus has extensive conformational heterogeneity and interacts with B-lineage genes

To produce a diverse antibody repertoire, immunoglobulin heavy-chain (Igh) loci undergo large-scale alterations in structure to facilitate juxtaposition and recombination of spatially separated variable (VH), diversity (DH), and joining (JH) genes. These chromosomal alterations are poorly understood. Uncovering their patterns shows how chromosome dynamics underpins antibody diversity. Using tiled Capture Hi-C, we produce a comprehensive map of chromatin interactions throughout the 2.8-Mb Igh locus in progenitor B cells. We find that the Igh locus folds into semi-rigid subdomains and undergoes flexible looping of the VH genes to its 3' end, reconciling two views of locus organization. Deconvolution of single Igh locus conformations using polymer simulations identifies thousands of different structures. This heterogeneity may underpin the diversity of V(D)J recombination events. All three immunoglobulin loci also participate in a highly specific, developmentally regulated network of interchromosomal interactions with genes encoding B cell-lineage factors. This suggests a model of interchromosomal coordination of B cell development.

Human and mouse early B cell development: So similar but so different

Early B cell development in the bone marrow ensures the replenishment of the peripheral B cell pool. Immature B cells continuously develop from hematopoietic stem cells, in a process guided by an intricate network of transcription factors as well as chemokine and cytokine signals. Humans and mice possess somewhat similar regulatory mechanisms of B lymphopoiesis. The continuous discovery of monogenetic defects that impact early B cell development in humans substantiates the similarities and differences with B cell development in mice. These differences become relevant when targeted therapeutic approaches are used in patients; therefore, predicting potential immunological adverse events is crucial. In this review, we have provided a phenotypical classification of human and murine early progenitors and B cell stages, based on surface and intracellular protein expression. Further, we have critically compared the role of key transcription factors (Ikaros, E2A, EBF1, PAX5, and Aiolos) and chemo- or cytokine signals (FLT3, c-kit, IL-7R, and CXCR4) during homeostatic and aberrant B lymphopoiesis in both humans and mice.

Transcription factor networks link B-lymphocyte development and malignant transformation in leukemia

Rapid advances in genomics have opened unprecedented possibilities to explore the mutational landscapes in malignant diseases, such as B-cell acute lymphoblastic leukemia (B-ALL). This disease is manifested as a severe defect in the production of normal blood cells due to the uncontrolled expansion of transformed B-lymphocyte progenitors in the bone marrow. Even though classical genetics identified translocations of transcription factor-coding genes in B-ALL, the extent of the targeting of regulatory networks in malignant transformation was not evident until the emergence of large-scale genomic analyses. There is now evidence that many B-ALL cases present with mutations in genes that encode transcription factors with critical roles in normal B-lymphocyte development. These include PAX5, IKZF1, EBF1, and TCF3, all of which are targeted by translocations or, more commonly, partial inactivation in cases of B-ALL. Even though there is support for the notion that germline polymorphisms in the PAX5 and IKZF1 genes predispose for B-ALL, the majority of leukemias present with somatic mutations in transcription factor-encoding genes. These genetic aberrations are often found in combination with mutations in genes that encode components of the pre-B-cell receptor or the IL-7/TSLP signaling pathways, all of which are important for early B-cell development. This review provides an overview of our current understanding of the molecular interplay that occurs between transcription factors and signaling events during normal and malignant B-lymphocyte development.

Profiling of miRNAs and their interfering targets in peripheral blood mononuclear cells from patients with chronic myeloid leukaemia

Introduction

MicroRNAs may be implicated in the acquisition of drug resistance in chronic myeloid leukemia as they regulate the expression of not only BCR-ABL1 but also genes associated with the activation of drug transfer proteins or essential signaling pathways.

Methods

To understand the impact of specifically expressed miRNAs in chronic myeloid leukemia and their target genes, we collected peripheral blood mononuclear cells (PBMC) from patients diagnosed with chronic myeloid leukemia (CML) and healthy donors to determine whole miRNA expression by small RNA sequencing and screened out 31 differentially expressed microRNAs (DE-miRNAs) with high expression. With the utilization of miRNA set enrichment analysis tools, we present here a comprehensive analysis of the relevance of DE-miRNAs to disease and biological function. Furthermore, the literature-based miRNA-target gene database was used to analyze the overall target genes of the DE-miRNAs and to define their associated biological responses. We further integrated DE-miRNA target genes to identify CML miRNA targeted gene signature singscore (CMTGSS) and used gene-set enrichment analysis (GSEA) to analyze the correlation between CMTGSS and Hallmark gene-sets in PBMC samples from clinical CML patients. Finally, the association of CMTGSS stratification with multiple CML cell lineage gene sets was validated in PBMC samples from CML patients using GSEA.

Results

Although individual miRNAs have been reported to have varying degrees of impact on CML, overall, our results show that abnormally upregulated miRNAs are associated with apoptosis and aberrantly downregulated miRNAs are associated with cell cycle. The clinical database shows that our defined DE-miRNAs are associated with the prognosis of CML patients. CMTGSS-based stratification analysis presented a tendency for miRNAs to affect cell differentiation in the blood microenvironment.

Conclusion

Collectively, this study defined differentially expressed miRNAs by miRNA sequencing from clinical samples and comprehensively analyzed the biological functions of the differential miRNAs in association with the target genes. The analysis of the enrichment of specific myeloid differentiated cells and immune cells also suggests the magnitude and potential targets of differentially expressed miRNAs in the clinical setting. It helps us to make links between the different results obtained from the multi-faceted studies to provide more potential research directions.

An intrinsically interpretable neural network architecture for sequence-to-function learning

MOTIVATION

Sequence-based deep learning approaches have been shown to predict a multitude of functional genomic readouts, including regions of open chromatin and RNA expression of genes. However, a major limitation of current methods is that model interpretation relies on computationally demanding post hoc analyses, and even then, one can often not explain the internal mechanics of highly parameterized models. Here, we introduce a deep learning architecture called totally interpretable sequence-to-function model (tiSFM). tiSFM improves upon the performance of standard multilayer convolutional models while using fewer parameters. Additionally, while tiSFM is itself technically a multilayer neural network, internal model parameters are intrinsically interpretable in terms of relevant sequence motifs.

RESULTS

We analyze published open chromatin measurements across hematopoietic lineage cell-types and demonstrate that tiSFM outperforms a state-of-the-art convolutional neural network model custom-tailored to this dataset. We also show that it correctly identifies context-specific activities of transcription factors with known roles in hematopoietic differentiation, including Pax5 and Ebf1 for B-cells, and Rorc for innate lymphoid cells. tiSFM's model parameters have biologically meaningful interpretations, and we show the utility of our approach on a complex task of predicting the change in epigenetic state as a function of developmental transition.

AVAILABILITY AND IMPLEMENTATION

The source code, including scripts for the analysis of key findings, can be found at https://github.com/boooooogey/ATAConv, implemented in Python.

PAX5 fusion genes in acute lymphoblastic leukemia: A literature review

Acute lymphoblastic leukemia (ALL) is a common cancer affecting children worldwide. The development of ALL is driven by several genes, some of which can be targeted for treatment by inhibiting gene fusions. PAX5 is frequently mutated in ALL and is involved in chromosomal rearrangements and translocations. Mutations in PAX5 interact with other genes, such as ETV6 and FOXP1, which influence B-cell development. PAX5/ETV6 has been observed in both B-ALL patients and a mouse model. The interaction between PAX5 and FOXP1 negatively suppresses the Pax5 gene in B-ALL patients. Additionally, ELN and PML genes have been found to fuse with PAX5, leading to adverse effects on B-cell differentiation. ELN-PAX5 interaction results in the decreased expression of LEF1, MB1, and BLNK, while PML-PAX5 is critical in the early stages of leukemia. PAX5 fusion genes prevent the transcription of the PAX5 gene, making it an essential target gene for the study of leukemia progression and the diagnosis of B-ALL.

An intrinsically interpretable neural network architecture for sequence to function learning

Motivation

Sequence-based deep learning approaches have been shown to predict a multitude of functional genomic readouts, including regions of open chromatin and RNA expression of genes. However, a major limitation of current methods is that model interpretation relies on computationally demanding post-hoc analyses, and even then, we often cannot explain the internal mechanics of highly parameterized models. Here, we introduce a deep learning architecture called tiSFM (totally interpretable sequence to function model). tiSFM improves upon the performance of standard multi-layer convolutional models while using fewer parameters. Additionally, while tiSFM is itself technically a multi-layer neural network, internal model parameters are intrinsically interpretable in terms of relevant sequence motifs.

Results

tiSFM's model architecture makes use of convolutions with a fixed set of kernel weights representing known transcription factor (TF) binding site motifs. We analyze published open chromatin measurements across hematopoietic lineage cell-types and demonstrate that tiSFM outperforms a state- of-the-art convolutional neural network model custom-tailored to this dataset. We also show that it correctly identifies context specific activities of transcription factors with known roles in hematopoietic differentiation, including Pax5 and Ebf1 for B-cells, and Rorc for innate lymphoid cells. tiSFM's model parameters have biologically meaningful interpretations, and we show the utility of our approach on a complex task of predicting the change in epigenetic state as a function of developmental transition.

Availability and implementation

The source code, including scripts for the analysis of key findings, can be found at https://github.com/boooooogey/ATAConv , implemented in Python.

Contact

atb44@pitt.edu.

EBF1 primes B-lymphoid enhancers and limits the myeloid bias in murine multipotent progenitors

Hematopoietic stem cells (HSCs) and multipotent progenitors (MPPs) generate all cells of the blood system. Despite their multipotency, MPPs display poorly understood lineage bias. Here, we examine whether lineage-specifying transcription factors, such as the B-lineage determinant EBF1, regulate lineage preference in early progenitors. We detect low-level EBF1 expression in myeloid-biased MPP3 and lymphoid-biased MPP4 cells, coinciding with expression of the myeloid determinant C/EBPα. Hematopoietic deletion of Ebf1 results in enhanced myelopoiesis and reduced HSC repopulation capacity. Ebf1-deficient MPP3 and MPP4 cells exhibit an augmented myeloid differentiation potential and a transcriptome with an enriched C/EBPα signature. Correspondingly, EBF1 binds the Cebpa enhancer, and the deficiency and overexpression of Ebf1 in MPP3 and MPP4 cells lead to an up- and downregulation of Cebpa expression, respectively. In addition, EBF1 primes the chromatin of B-lymphoid enhancers specifically in MPP3 cells. Thus, our study implicates EBF1 in regulating myeloid/lymphoid fate bias in MPPs by constraining C/EBPα-driven myelopoiesis and priming the B-lymphoid fate.

PAX5 alterations in B-cell acute lymphoblastic leukemia

PAX5, a master regulator of B cell development and maintenance, is one of the most common targets of genetic alterations in B-cell acute lymphoblastic leukemia (B-ALL). PAX5 alterations consist of copy number variations (whole gene, partial, or intragenic), translocations, and point mutations, with distinct distribution across B-ALL subtypes. The multifaceted functional impacts such as haploinsufficiency and gain-of-function of PAX5 depending on specific variants have been described, thereby the connection between the blockage of B cell development and the malignant transformation of normal B cells has been established. In this review, we provide the recent advances in understanding the function of PAX5 in orchestrating the development of both normal and malignant B cells over the past decade, with a focus on the PAX5 alterations shown as the initiating or driver events in B-ALL. Recent large-scale genomic analyses of B-ALL have identified multiple novel subtypes driven by PAX5 genetic lesions, such as the one defined by a distinct gene expression profile and PAX5 P80R mutation, which is an exemplar leukemia entity driven by a missense mutation. Although altered PAX5 is shared as a driver in B-ALL, disparate disease phenotypes and clinical outcomes among the patients indicate further heterogeneity of the underlying mechanisms and disturbed gene regulation networks along the disease development. In-depth mechanistic studies in human B-ALL and animal models have demonstrated high penetrance of PAX5 variants alone or concomitant with other genetic lesions in driving B-cell malignancy, indicating the altered PAX5 and deregulated genes may serve as potential therapeutic targets in certain B-ALL cases.

Termination of CD40L co-stimulation promotes human B cell differentiation into antibody-secreting cells

Human naïve B cells are notoriously difficult to differentiate into antibody-secreting cells (ASCs) in vitro while maintaining sufficient cell numbers to evaluate the differentiation process. B cells require T follicular helper (T_FH ) cell-derived signals like CD40L and IL-21 during germinal center (GC) responses to undergo differentiation into ASCs. Cognate interactions between B and T_FH cells are transient; after T_FH contact, B cells cycle between GC light and dark zones where T_FH contact is present and absent, respectively. Here, we elucidated that the efficacy of naïve B cells in ACS differentiation is dramatically enhanced by the release of CD40L stimulation. Multiparameter phospho-flow and transcription factor (TF)-flow cytometry revealed that termination of CD40L stimulation downmodulates NF-κB and STAT3 signaling. Furthermore, the termination of CD40 signaling downmodulates C-MYC, while promoting ASC TFs BLIMP1 and XBP-1s. Reduced levels of C-MYC in the differentiating B cells are later associated with crucial downmodulation of the B cell signature TF PAX5 specifically upon the termination of CD40 signaling, resulting in the differentiation of BLIMP1 high expressing cells into ASCs. The data presented here are the first steps to provide further insights how the transient nature of CD40 signaling is in fact needed for efficient human naïve B cell differentiation to ASCs.

EBNA2-EBF1 complexes promote MYC expression and metabolic processes driving S-phase progression of Epstein-Barr virus-infected B cells

Epstein-Barr virus (EBV) is a human tumor virus which preferentially infects resting human B cells. Upon infection in vitro, EBV activates and immortalizes these cells. The viral latent protein EBV nuclear antigen 2 (EBNA2) is essential for B cell activation and immortalization; it targets and binds the cellular and ubiquitously expressed DNA-binding protein CBF1, thereby transactivating a plethora of viral and cellular genes. In addition, EBNA2 uses its N-terminal dimerization (END) domain to bind early B cell factor 1 (EBF1), a pioneer transcription factor specifying the B cell lineage. We found that EBNA2 exploits EBF1 to support key metabolic processes and to foster cell cycle progression of infected B cells in their first cell cycles upon activation. The α1-helix within the END domain was found to promote EBF1 binding. EBV mutants lacking the α1-helix in EBNA2 can infect and activate B cells efficiently, but activated cells fail to complete the early S phase of their initial cell cycle. Expression of MYC, target genes of MYC and E2F, as well as multiple metabolic processes linked to cell cycle progression are impaired in EBVΔα1-infected B cells. Our findings indicate that EBF1 controls B cell activation via EBNA2 and, thus, has a critical role in regulating the cell cycle of EBV-infected B cells. This is a function of EBF1 going beyond its well-known contribution to B cell lineage specification.

Whole-Transcriptome Profiling and circRNA-miRNA-mRNA Regulatory Networks in B-Cell Development

The generation and differentiation of B lymphocytes (B cells) is a flexible process with many critical regulatory factors. Previous studies indicated that non-coding RNAs play multiple roles in the development of lymphocytes. However, little has been known about the circular RNA (circRNA) profiles and their competing endogenous RNA (ceRNA) networks in B-cell development and differentiation. Here, four B-cell subsets were purified from single-cell suspensions of mouse bone marrow. Then RNA sequencing (RNA-Seq) was used to display expression profiles of circRNAs, miRNAs and mRNAs during B-cell differentiation. 175, 203, 219 and 207 circRNAs were specifically expressed in pro-B cells, pre-B cells, immature B cells and mature B cells, respectively. The circRNA-associated ceRNA networks constructed in two sequential stages of B-cell differentiation revealed the potential mechanism of circRNAs in these processes. This study is the first to explore circRNA profiles and circRNA-miRNA-mRNA networks in different B-cell developmental stages of mouse bone marrow, which contribute to further research on their mechanism in B-cell development and differentiation.

Chromodomain helicase DNA-binding 4 (CHD4) regulates early B cell identity and V(D)J recombination

B lymphocytes develop from uncommitted precursors into immunoglobulin (antibody)-producing B cells, a major arm of adaptive immunity. Progression of early progenitors to antibody-expressing cells in the bone marrow is orchestrated by the temporal regulation of different gene programs at discrete developmental stages. A major question concerns how B cells control the accessibility of these genes to transcription factors. Research has implicated nucleosome remodeling ATPases as mediators of chromatin accessibility. Here, we describe studies of chromodomain helicase DNA-binding 4 (CHD4; also known as Mi-2β) in early B cell development. CHD4 comprises multiple domains that function in nucleosome mobilization and histone binding. CHD4 is a key component of Nucleosome Remodeling and Deacetylase, or NuRD (Mi-2) complexes, which assemble with other proteins that mediate transcriptional repression. We review data demonstrating that CHD4 is necessary for B lineage identity: early B lineage progression, proliferation in response to interleukin-7, responses to DNA damage, and cell survival in vivo. CHD4-NuRD is also required for the Ig heavy-chain repertoire by promoting utilization of distal variable (V_H ) gene segments in V(D)J recombination. In conclusion, the regulation of chromatin accessibility by CHD4 is essential for production of antibodies by B cells, which in turn mediate humoral immune responses to pathogens and disease.

PAX5 haploinsufficiency induced CD8+ T cells dysfunction or exhaustion by high expression of immune inhibitory-related molecules

PURPOSE

PAX5 haploinsufficiency promoting tumorigenesis is related to immune escape. But the mechanisms of PAX5 mutations inducing tumor immune escape have not been clarified. Our aim was to study how PAX5 haploinsufficiency influences effector CD8 + T cells in tumor microenvironment.

METHODS

We estimated the proportions of 22 immune cell types and the expressions of immune inhibitory-related molecules based on gene expression profiles (GEPs) from children's B- acute lymphoblastic leukemia(B-ALL) with PAX5 mutations by CIBERSORT, an established algorithm. We constructed the PAX5 haplodeletion A20 cell lines, built allografted A20 tumor models and evaluated the effect of PAX5 haplodeletion on immune inhibitory-related molecules in the tumor microenvironment (TME).

RESULTS

Our results indicated the percentages of T cells in bone marrow of children's B-ALL with PAX5 mutations were not statistically different from that in bone marrow of B-ALL without PAX5 mutations, except for T follicular helper (Tfh) cells. But a variety of up-regulated immune inhibitory-related molecules in bone marrow of children's B- ALL with PAX5 mutations were identified. By different approaches, we found that several immune inhibitory-related molecules of CD8+ T cells in TME of PAX5 haplodeletion clones such as TIM3, NR4A1 and BATF, were increased significantly compared with that of PAX5 wild type control. The IFN-ɤ of CD8+ T cells in TME of PAX5 haplodeletion tumors was decreased significantly compared with that of PAX5 wild type control.

CONCLUSION

Our study showed that PAX5 haploinsufficiency induced CD8+ T cells dysfunction or exhaustion by high expression of TIM3, NR4A1 and BATF in the CD8+ T cells of TME.

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.

B Lymphocyte Specification Is Preceded by Extensive Epigenetic Priming in Multipotent Progenitors

B lymphocyte development is dependent on the interplay between the chromatin landscape and lineage-specific transcription factors. It has been suggested that B lineage commitment is associated with major changes in the nuclear chromatin environment, proposing a critical role for lineage-specific transcription factors in the formation of the epigenetic landscape. In this report, we have used chromosome conformation capture in combination with assay for transposase-accessible chromatin sequencing analysis to enable highly efficient annotation of both proximal and distal transcriptional control elements to genes activated in B lineage specification in mice. A large majority of these genes were annotated to at least one regulatory element with an accessible chromatin configuration in multipotent progenitors. Furthermore, the majority of binding sites for the key regulators of B lineage specification, EBF1 and PAX5, occurred in already accessible regions. EBF1 did, however, cause a dynamic change in assay for transposase-accessible chromatin accessibility and was critical for an increase in distal promoter-enhancer interactions. Our data unravel an extensive epigenetic priming at regulatory elements annotated to lineage-restricted genes and provide insight into the interplay between the epigenetic landscape and transcription factors in cell specification.