Single-cell analysis identifies dynamic gene expression networks that govern B cell development and transformation.

Integration of external signals and B-lymphoid transcription factor activities organise B cell lineage commitment through alternating cycles of proliferation and differentiation, producing a diverse repertoire of mature B cells. We use single-cell transcriptomics/proteomics to identify differentially expressed gene networks across B cell development and correlate these networks with subtypes of B cell leukemia. Here we show unique transcriptional signatures that refine the pre-B cell expansion stages into pre-BCR-dependent and pre-BCR-independent proliferative phases. These changes correlate with reciprocal changes in expression of the transcription factor EBF1 and the RNA binding protein YBX3, that are defining features of the pre-BCR-dependent stage. Using pseudotime analysis, we further characterize the expression kinetics of different biological modalities across B cell development, including transcription factors, cytokines, chemokines, and their associated receptors. Our findings demonstrate the underlying heterogeneity of developing B cells and characterise developmental nodes linked to B cell transformation.

Ku70 suppresses alternative end joining in G1-arrested progenitor B cells.

Classical nonhomologous end joining (C-NHEJ) repairs DNA double-strand breaks (DSBs) throughout interphase but predominates in G1 phase when homologous recombination is unavailable. Complexes containing the Ku70/80 ("Ku") and XRCC4/ligase IV (Lig4) core C-NHEJ factors are required, respectively, for sensing and joining DSBs. While XRCC4/Lig4 are absolutely required for joining RAG1/2 endonuclease ("RAG")-initiated DSBs during V(D)J recombination in G1-phase progenitor lymphocytes, cycling cells deficient for XRCC4/Lig4 also can join chromosomal DSBs by alternative end-joining (A-EJ) pathways. Restriction of V(D)J recombination by XRCC4/Lig4-mediated joining has been attributed to RAG shepherding V(D)J DSBs exclusively into the C-NHEJ pathway. Here, we report that A-EJ of DSB ends generated by RAG1/2, Cas9:gRNA, and Zinc finger endonucleases in Lig4-deficient G1-arrested progenitor B cell lines is suppressed by Ku. Thus, while diverse DSBs remain largely as free broken ends in Lig4-deficient G1-arrested progenitor B cells, deletion of Ku70 increases DSB rejoining and translocation levels to those observed in Ku70-deficient counterparts. Correspondingly, while RAG-initiated V(D)J DSB joining is abrogated in Lig4-deficient G1-arrested progenitor B cell lines, joining of RAG-generated DSBs in Ku70-deficient and Ku70/Lig4 double-deficient lines occurs through a translocation-like A-EJ mechanism. Thus, in G1-arrested, Lig4-deficient progenitor B cells are functionally end-joining suppressed due to Ku-dependent blockage of A-EJ, potentially in association with G1-phase down-regulation of Lig1. Finally, we suggest that differential impacts of Ku deficiency versus Lig4 deficiency on V(D)J recombination, neuronal apoptosis, and embryonic development results from Ku-mediated inhibition of A-EJ in the G1 cell cycle phase in Lig4-deficient developing lymphocyte and neuronal cells.

Notch2-mediated plasticity between marginal zone and follicular B cells.

Follicular B (FoB) and marginal zone B (MZB) cells are functionally and spatially distinct mature B cell populations in the spleen, originating from a Notch2-dependent fate decision after splenic influx of immature transitional B cells. In the B cell follicle, a Notch2-signal is provided by DLL-1-expressing fibroblasts. However, it is unclear whether FoB cells, which are in close contact with these DLL-1 expressing fibroblasts, can also differentiate to MZB cells if they receive a Notch2-signal. Here, we show induced Notch2IC-expression in FoB cells re-programs mature FoB cells into bona fide MZB cells as is evident from the surface phenotype, localization, immunological function and transcriptome of these cells. Furthermore, the lineage conversion from FoB to MZB cells occurs in immunocompetent wildtype mice. These findings demonstrate plasticity between mature FoB and MZB cells that can be driven by a singular signaling event, the activation of Notch2.

Chromosomal coordination and differential structure of asynchronous replicating regions.

Stochastic asynchronous replication timing (AS-RT) is a phenomenon in which the time of replication of each allele is different, and the identity of the early allele varies between cells. By taking advantage of stable clonal pre-B cell populations derived from C57BL6/Castaneous mice, we have mapped the genome-wide AS-RT loci, independently of genetic differences. These regions are characterized by differential chromatin accessibility, mono-allelic expression and include new gene families involved in specifying cell identity. By combining population level mapping with single cell FISH, our data reveal the existence of a novel regulatory program that coordinates a fixed relationship between AS-RT regions on any given chromosome, with some loci set to replicate in a parallel and others set in the anti-parallel orientation. Our results show that AS-RT is a highly regulated epigenetic mark established during early embryogenesis that may be used for facilitating the programming of mono-allelic choice throughout development.

Ikaros antagonizes DNA binding by STAT5 in pre-B cells.

The IKZF1 gene, which encodes the Ikaros transcription factor, is frequently deleted or mutated in patients with B-cell precursor acute lymphoblastic leukemias that express oncogenes, like BCR-ABL, which activate the JAK-STAT5 pathway. Ikaros functionally antagonizes the transcriptional programs downstream of IL-7/STAT5 during B cell development, as well as STAT5 activity in leukemic cells. However, the mechanisms by which Ikaros interferes with STAT5 function is unknown. We studied the genomic distribution of Ikaros and STAT5 on chromatin in a murine pre-B cell line, and found that both proteins colocalize on >60% of STAT5 target regions. Strikingly, Ikaros activity leads to widespread loss of STAT5 binding at most of its genomic targets within two hours of Ikaros induction, suggesting a direct mechanism. Ikaros did not alter the level of total or phosphorylated STAT5 proteins, nor did it associate with STAT5. Using sequences from the Cish, Socs2 and Bcl6 genes that Ikaros and STAT5 target, we show that both proteins bind overlapping sequences at GGAA motifs. Our results demonstrate that Ikaros antagonizes STAT5 DNA binding, in part by competing for common target sequences. Our study has implications for understanding the functions of Ikaros and STAT5 in B cell development and transformation.

The Ig heavy chain protein but not its message controls early B cell development.

Development of progenitor B cells (ProB cells) into precursor B cells (PreB cells) is dictated by immunoglobulin heavy chain checkpoint (IgHCC), where the IgHC encoded by a productively rearranged Igh allele assembles into a PreB cell receptor complex (PreBCR) to generate signals to initiate this transition and suppressing antigen receptor gene recombination, ensuring that only one productive Igh allele is expressed, a phenomenon known as Igh allelic exclusion. In contrast to a productively rearranged Igh allele, the Igh messenger RNA (mRNA) (IgHR) from a nonproductively rearranged Igh allele is degraded by nonsense-mediated decay (NMD). This fact prohibited firm conclusions regarding the contribution of stable IgHR to the molecular and developmental changes associated with the IgHCC. This point was addressed by generating the IghTer5H∆TM mouse model from IghTer5H mice having a premature termination codon at position +5 in leader exon of IghTer5H allele. This prohibited NMD, and the lack of a transmembrane region (∆TM) prevented the formation of any signaling-competent PreBCR complexes that may arise as a result of read-through translation across premature Ter5 stop codon. A highly sensitive sandwich Western blot revealed read-through translation of IghTer5H message, indicating that previous conclusions regarding a role of IgHR in establishing allelic exclusion requires further exploration. As determined by RNA sequencing (RNA-Seq), this low amount of IgHC sufficed to initiate PreB cell markers normally associated with PreBCR signaling. In contrast, the IghTer5H∆TM knock-in allele, which generated stable IgHR but no detectable IgHC, failed to induce PreB development. Our data indicate that the IgHCC is controlled at the level of IgHC and not IgHR expression.

A Prion-like Domain in Transcription Factor EBF1 Promotes Phase Separation and Enables B Cell Programming of Progenitor Chromatin.

Establishment of B-lineage-specific gene expression requires the binding of transcription factors to inaccessible chromatin of progenitors. The transcription factor EBF1 can bind genomic regions prior to the detection of chromatin accessibility in a manner dependent on EBF1's C-terminal domain (CTD) and independent of cooperating transcription factors. Here, we studied the mechanism whereby the CTD enables this pioneering function. The CTD of EBF1 was dispensable for initial chromatin targeting but stabilized occupancy via recruitment of the chromatin remodeler Brg1. We found that the CTD harbors a prion-like domain (PLD) with an ability of liquid-liquid phase separation, which was enhanced by interaction of EBF1 with the RNA-binding protein FUS. Brg1 also partitioned into phase-separated FUS condensates and coincided with EBF1 and FUS foci in pro-B cells. Heterologous PLDs conferred pioneering function on EBF1ΔCTD. Thus, the phase separation ability of EBF1 facilitates Brg1-mediated chromatin opening and the transition of naive progenitor chromatin to B-lineage-committed chromatin.

Wapl repression by Pax5 promotes V gene recombination by Igh loop extrusion.

Nuclear processes, such as V(D)J recombination, are orchestrated by the three-dimensional organization of chromosomes at multiple levels, including compartments1 and topologically associated domains (TADs)2,3 consisting of chromatin loops4. TADs are formed by chromatin-loop extrusion5-7, which depends on the loop-extrusion function of the ring-shaped cohesin complex8-12. Conversely, the cohesin-release factor Wapl13,14 restricts loop extension10,15. The generation of a diverse antibody repertoire, providing humoral immunity to pathogens, requires the participation of all V genes in V(D)J recombination16, which depends on contraction of the 2.8-Mb-long immunoglobulin heavy chain (Igh) locus by Pax517,18. However, how Pax5 controls Igh contraction in pro-B cells remains unknown. Here we demonstrate that locus contraction is caused by loop extrusion across the entire Igh locus. Notably, the expression of Wapl is repressed by Pax5 specifically in pro-B and pre-B cells, facilitating extended loop extrusion by increasing the residence time of cohesin on chromatin. Pax5 mediates the transcriptional repression of Wapl through a single Pax5-binding site by recruiting the polycomb repressive complex 2 to induce bivalent chromatin at the Wapl promoter. Reduced Wapl expression causes global alterations in the chromosome architecture, indicating that the potential to recombine all V genes entails structural changes of the entire genome in pro-B cells.

Structural Insights into Pseudokinase Domains of Receptor Tyrosine Kinases.

Despite their apparent lack of catalytic activity, pseudokinases are essential signaling molecules. Here, we describe the structural and dynamic properties of pseudokinase domains from the Wnt-binding receptor tyrosine kinases (PTK7, ROR1, ROR2, and RYK), which play important roles in development. We determined structures of all pseudokinase domains in this family and found that they share a conserved inactive conformation in their activation loop that resembles the autoinhibited insulin receptor kinase (IRK). They also have inaccessible ATP-binding pockets, occluded by aromatic residues that mimic a cofactor-bound state. Structural comparisons revealed significant domain plasticity and alternative interactions that substitute for absent conserved motifs. The pseudokinases also showed dynamic properties that were strikingly similar to those of IRK. Despite the inaccessible ATP site, screening identified ATP-competitive type-II inhibitors for ROR1. Our results set the stage for an emerging therapeutic modality of "conformational disruptors" to inhibit or modulate non-catalytic functions of pseudokinases deregulated in disease.

TBL1XR1 Mutations Drive Extranodal Lymphoma by Inducing a Pro-tumorigenic Memory Fate.

The most aggressive B cell lymphomas frequently manifest extranodal distribution and carry somatic mutations in the poorly characterized gene TBL1XR1. Here, we show that TBL1XR1 mutations skew the humoral immune response toward generating abnormal immature memory B cells (MB), while impairing plasma cell differentiation. At the molecular level, TBL1XR1 mutants co-opt SMRT/HDAC3 repressor complexes toward binding the MB cell transcription factor (TF) BACH2 at the expense of the germinal center (GC) TF BCL6, leading to pre-memory transcriptional reprogramming and cell-fate bias. Upon antigen recall, TBL1XR1 mutant MB cells fail to differentiate into plasma cells and instead preferentially reenter new GC reactions, providing evidence for a cyclic reentry lymphomagenesis mechanism. Ultimately, TBL1XR1 alterations lead to a striking extranodal immunoblastic lymphoma phenotype that mimics the human disease. Both human and murine lymphomas feature expanded MB-like cell populations, consistent with a MB-cell origin and delineating an unforeseen pathway for malignant transformation of the immune system.

The role of metabolic checkpoint regulators in B cell survival and transformation.

In response to mitogenic stimulation, B cells activate different pro-anabolic signaling pathways such as c-Myc- and mTORC1-dependent networks to satisfy the energetic demands of biomass synthesis and proliferation. In order to preserve viability and function, cell growth cannot progress unchecked and must be adjusted according to the availability of nutrients. Nutrient-sensing proteins such as AMPK antagonize mTORC1 activity in response to starvation. If pro-anabolic signaling pathways are aberrantly activated, B cells may lack the metabolic capacity to accommodate their energetic needs, which can lead to cell death. On the other hand, metabolic hyperactivation is a salient feature of cancer cells, suggesting that mechanisms exist, which allow B cells to cope with metabolic stress. The aim of this review is to discuss how B cells respond to a mismatch between energy supply and demand and what the consequences are of metabolic dysregulation in normal and malignant B cells.

Increased O-GlcNAcylation of c-Myc Promotes Pre-B Cell Proliferation.

O-linked ß-N-acetylglucosamine (O-GlcNAc) modification regulates the activity of hundreds of nucleocytoplasmic proteins involved in a wide variety of cellular processes, such as gene expression, signaling, and cell growth; however, the mechanism underlying the regulation of B cell development and function by O-GlcNAcylation remains largely unknown. Here, we demonstrate that changes in cellular O-GlcNAc levels significantly affected the growth of pre-B cells, which rapidly proliferate to allow expansion of functional clones that express successfully rearranged heavy chains at the pro-B stage during early B cell development. In our study, the overall O-GlcNAc levels in these proliferative pre-B cells, which are linked to the glucose uptake rate, were highly induced when compared with those in pro-B cells. Thus, pharmacologically, genetically, or nutritionally, inhibition of O-GlcNAcylation in pre-B cells markedly downregulated c-Myc expression, resulting in cell cycle arrest via blockade of cyclin expression. Importantly, the population of B cells after the pro-B cell stage in mouse bone marrow was severely impaired by the administration of an O-GlcNAc inhibitor. These results strongly suggest that O-GlcNAcylation-dependent expression of c-Myc represents a new regulatory component of pre-B cell proliferation, as well as a potential therapeutic target for the treatment of pre-B cell-derived leukemia.

An R-Derived FlowSOM Process to Analyze Unsupervised Clustering of Normal and Malignant Human Bone Marrow Classical Flow Cytometry Data.

Multiparameter flow cytometry (MFC) is a powerful and versatile tool to accurately analyze cell subsets, notably to explore normal and pathological hematopoiesis. Yet, mostly supervised subjective strategies are used to identify cell subsets in this complex tissue. In the past few years, the implementation of mass cytometry and the big data generated have led to a blossoming of new software solutions. Their application to classical MFC in hematology is however still seldom reported. Here, we show how one of these new tools, the FlowSOM R solution, can be applied, together with the Kaluza® software, to a new delineation of hematopoietic subsets in normal human bone marrow (BM). We thus combined the unsupervised discrimination of cell subsets provided by FlowSOM and their expert-driven node-by-node assignment to known or new hematopoietic subsets. We also show how this new tool could modify the MFC exploration of hematological malignancies both at diagnosis (Dg) and follow-up (FU). This can be achieved by direct comparison of merged listmodes of reference normal BM, Dg, and FU samples of a representative acute myeloblastic case tested with the same immunophenotyping panel. This provides an immediate unsupervised evaluation of minimal residual disease. © 2019 International Society for Advancement of Cytometry.

Discovery of a CD10-negative B-progenitor in human fetal life identifies unique ontogeny-related developmental programs.

Human lymphopoiesis is a dynamic lifelong process that starts in utero 6 weeks postconception. Although fetal B-lymphopoiesis remains poorly defined, it is key to understanding leukemia initiation in early life. Here, we provide a comprehensive analysis of the human fetal B-cell developmental hierarchy. We report the presence in fetal tissues of 2 distinct CD19+ B-progenitors, an adult-type CD10+ve ProB-progenitor and a new CD10-ve PreProB-progenitor, and describe their molecular and functional characteristics. PreProB-progenitors and ProB-progenitors appear early in the first trimester in embryonic liver, followed by a sustained second wave of B-progenitor development in fetal bone marrow (BM), where together they form >40% of the total hematopoietic stem cell/progenitor pool. Almost one-third of fetal B-progenitors are CD10-ve PreProB-progenitors, whereas, by contrast, PreProB-progenitors are almost undetectable (0.53% ± 0.24%) in adult BM. Single-cell transcriptomics and functional assays place fetal PreProB-progenitors upstream of ProB-progenitors, identifying them as the first B-lymphoid-restricted progenitor in human fetal life. Although fetal BM PreProB-progenitors and ProB-progenitors both give rise solely to B-lineage cells, they are transcriptionally distinct. As with their fetal counterparts, adult BM PreProB-progenitors give rise only to B-lineage cells in vitro and express the expected B-lineage gene expression program. However, fetal PreProB-progenitors display a distinct, ontogeny-related gene expression pattern that is not seen in adult PreProB-progenitors, and they share transcriptomic signatures with CD10-ve B-progenitor infant acute lymphoblastic leukemia blast cells. These data identify PreProB-progenitors as the earliest B-lymphoid-restricted progenitor in human fetal life and suggest that this fetal-restricted committed B-progenitor might provide a permissive cellular context for prenatal B-progenitor leukemia initiation.

Deletion of HO-1 blocks development of B lymphocytes in mice.

B lymphocytes, a key cluster of cells composing the immune system, can protect against abnormal biological factors. Heme oxygenase-1 (HO-1) plays important roles in cell proliferation and immune regulation, but its effects on the development and growth of B lymphocytes are still unknown. Herein, the count of B lymphocytes in HO-1 gene knockout (HO-1+/-) mice was significantly lower than that of the HO-1 gene wild-type (HO-1WT) mice. Meanwhile, the cell count of HO-1+/- mice did not recover after irradiation for one week, due to the G0/G1 phase arrest of Pro-B cells and the augmented apoptosis of Pre-B cells. Up-regulation of HO-1 by lentivirus attenuated the Pro-B cell cycle arrest and Pre-B cell apoptosis. To understand the molecular mechanism by which HO-1 knockout blocked B lymphocyte development, protein-to-protein interaction network and Western blot were used. The PI3K/AKT signaling pathway mediated the regulatory effects of HO-1 on B lymphocytes. In conclusion, HO-1 is a crucial transcriptional repressor for B cell development.

The Roles of Bursal Nonapeptide (BP9) on AIV Vaccine Immune Response in Chick Immunization and on Avian Immature B Cell.

BACKGROUND: Bursa of Fabricius plays the vital functions on B cell development and antibody production in poultry. The bursal-derived peptide plays the essential roles on avian immature B cell development. OBJECTIVES: Here we explored the functions of the recently reported bursal nonapeptide (BP9) on the antibody production and the molecular basis of BP9 on avian immature B cell. METHODS: Chicken were twice immunized with Avian Influenza Virus (AIV) inactivated vaccine plus with BP9 at three dosages, respectively. On two weeks after the second immunization, sera samples were collected from all experimental groups to measure AIV-specific Agglutination Inhibition (HI) antibody titers. Also, on 7th day after the second immunization, spleen lymphocytes were isolated from the immunized chicken to detect the lymphocyte viabilities. DT40 cells were treated with BP9 from 0.02 to 2 µg/mL for 4 and 20h to detect sIgM mRNA levels, and total RNAs from BP9-treated DT40 cells were collected to investigate the gene expression profiles of DT40 cells, and to analyze the enriched pathways and functional biological processes. Finally, nine gene expressions were validated with quantitative PCR (qPCR). RESULTS: Our investigation proved the strong regulatory roles of BP9 on AIV-specific HI antibody titers and lymphocyte viabilities. BP9 promoted sIgM mRNA levels in DT40 cells, and upregulated 598 gene expressions and downregulated 395 gene expressions in DT40 cells with 0.2µg/mL BP9 treatment. Moreover, our findings verified the significantly enriched six pathways and various the biological functional processes of BP9 on avian immature B cell. Also, we found eight signaling pathways in the enriched biological processes of BP9-treated DT40 cells, and the expressions of nine selected genes with qPCR were identical to that of microarray data. CONCLUSION: BP9 promoted the antibody production in the 21-old-day chicken immunization, and stimulated the sIgM expression in DT40 cells. Furthermore, we analyzed the gene expression profile and immune-related biological processes of DT40 cells treated with BP9, which provided some new insights into the mechanism on immature B cell development, and provided important references for adjuvant development on vaccine improvement and clinical application.

The Chimpanzee SIV Envelope Trimer: Structure and Deployment as an HIV Vaccine Template.

Epitope-targeted HIV vaccine design seeks to focus antibody responses to broadly neutralizing antibody (bnAb) sites by sequential immunization. A chimpanzee simian immunodeficiency virus (SIV) envelope (Env) shares a single bnAb site, the variable loop 2 (V2)-apex, with HIV, suggesting its possible utility in an HIV immunization strategy. Here, we generate a chimpanzee SIV Env trimer, MT145K, which displays selective binding to HIV V2-apex bnAbs and precursor versions, but no binding to other HIV specificities. We determine the structure of the MT145K trimer by cryo-EM and show that its architecture is remarkably similar to HIV Env. Immunization of an HIV V2-apex bnAb precursor Ab-expressing knockin mouse with the chimpanzee MT145K trimer induces HIV V2-specific neutralizing responses. Subsequent boosting with an HIV trimer cocktail induces responses that exhibit some virus cross-neutralization. Overall, the chimpanzee MT145K trimer behaves as expected from design both in vitro and in vivo and is an attractive potential component of a sequential immunization regimen to induce V2-apex bnAbs.

Feedforward regulation of Myc coordinates lineage-specific with housekeeping gene expression during B cell progenitor cell differentiation.

The differentiation of self-renewing progenitor cells requires not only the regulation of lineage- and developmental stage-specific genes but also the coordinated adaptation of housekeeping functions from a metabolically active, proliferative state toward quiescence. How metabolic and cell-cycle states are coordinated with the regulation of cell type-specific genes is an important question, because dissociation between differentiation, cell cycle, and metabolic states is a hallmark of cancer. Here, we use a model system to systematically identify key transcriptional regulators of Ikaros-dependent B cell-progenitor differentiation. We find that the coordinated regulation of housekeeping functions and tissue-specific gene expression requires a feedforward circuit whereby Ikaros down-regulates the expression of Myc. Our findings show how coordination between differentiation and housekeeping states can be achieved by interconnected regulators. Similar principles likely coordinate differentiation and housekeeping functions during progenitor cell differentiation in other cell lineages.