Transcription factor EBF1 is essential for the maintenance of B cell identity and prevention of alternative fates in committed cells

The significance of PAX5 in Merkel cell carcinoma

Merkel cell carcinoma (MCC) is a highly malignant skin cancer that expresses epithelial-, neuroendocrine-, and lymphoid-associated genes. Here, we focused on B-cell differentiation, which is characterised by the coexpression of PAX5 and TdT. PAX5 is the master regulator of B-cell commitment and is expressed in 65% of MCC cases. Yet little is known about the underlying molecular biology of the frequently reported PAX5 expression in MCC. Multi-omics analyses, including RNA next-generation sequencing, RT-qPCR, immunohistochemistry, and western blotting, were performed to assess PAX5 expression in MCC. Differential DNA methylation analysis at 61,043 PAX5 binding sites in enhancer and promoter elements was performed to detect differences between n = 14 MCC tissues and n = 91 various normal B-cell populations. RNA analysis revealed full-length PAX5 expression in MCC at the transcriptional level using both PAX5 transcription start sites. PAX5 protein expression was found in 40 of 41 MCCs and six out of seven MCC cell lines. DNA methylation array analysis revealed 1,084 hypermethylated loci of enhancer and promoter elements located in PAX5 binding sites in MCC. Of these, 702 loci were associated with 257 genes that are not expressed. The PAX5-associated regulatory elements of these 257 genes were enriched for interferon regulatory factor 4 (IRF4) and SPi-proto-oncogene (SPI1) binding motifs. Neither IRF4 or SPI1 could be detected in MCC on RNA or the protein level. Thus, because of the absence of these transcription factors, we conclude that full-length PAX5 alone cannot induce B-cell differentiation. © 2025 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Synergistic and antagonistic activities of IRF8 and FOS enhancer pairs during an immune-cell fate switch

Cell fate instructive genes tend to be regulated by large clusters of enhancers. Whether and how individual enhancers within such clusters cooperate in regulating gene expression is poorly understood. We have previously developed a computational method, SEGCOND, which identifies hubs that we termed Putative Transcriptional Condensates (PTCs), consisting of enhancer clusters and associated target genes. Here, we use SEGCOND to identify PTCs in a CEBPA-induced B-cell-to-macrophage transdifferentiation system. We find that PTCs are enriched for highly expressed, lineage-restricted genes and associate with BRD4, a component of transcriptional condensates. Further, we performed single and combinatorial deletions of enhancers within two PTCs active during induced transdifferentiation, harboring IRF8 and FOS. Two enhancers within the IRF8 PTC were found to provide a backup mechanism when combined, safeguarding IRF8 expression and efficient transdifferentiation. Unexpectedly, two individual enhancers within the FOS PTC antagonize each other on day 1 of transdifferentiation, delaying the conversion of B-cells into macrophages and reducing FOS expression, while on day 7, they cooperate to increase FOS levels induced cells. Our results reveal complex, differentiation-stage-specific interactions between individual enhancers within enhancer clusters.

Active repression of cell fate plasticity by PROX1 safeguards hepatocyte identity and prevents liver tumorigenesis

Cell fate plasticity enables development, yet unlocked plasticity is a cancer hallmark. While transcription master regulators induce lineage-specific genes to restrict plasticity, it remains unclear whether plasticity is actively suppressed by lineage-specific repressors. Here we computationally predict so-called safeguard repressors for 18 cell types that block phenotypic plasticity lifelong. We validated hepatocyte-specific candidates using reprogramming, revealing that prospero homeobox protein 1 (PROX1) enhanced hepatocyte identity by direct repression of alternative fate master regulators. In mice, Prox1 was required for efficient hepatocyte regeneration after injury and was sufficient to prevent liver tumorigenesis. In line with patient data, Prox1 depletion caused hepatocyte fate loss in vivo and enabled the transition of hepatocellular carcinoma to cholangiocarcinoma. Conversely, overexpression promoted cholangiocarcinoma to hepatocellular carcinoma transdifferentiation. Our findings provide evidence for PROX1 as a hepatocyte-specific safeguard and support a model where cell-type-specific repressors actively suppress plasticity throughout life to safeguard lineage identity and thus prevent disease.

Bone Marrow Niche in Cardiometabolic Disease: Mechanisms and Therapeutic Potential

Cardiovascular and cardiometabolic diseases are leading causes of morbidity and mortality worldwide, driven in part by chronic inflammation. Emerging research suggests that the bone marrow microenvironment, or marrow niche, plays a critical role in both immune system regulation and disease progression. The bone marrow niche is essential for maintaining hematopoietic stem cells (HSCs) and orchestrating hematopoiesis. Under normal conditions, this niche ensures a return to immune homeostasis after acute stress. However, in the setting of inflammatory conditions such as those seen in cardiometabolic diseases, it becomes dysregulated, leading to enhanced myelopoiesis and immune activation. This review explores the reciprocal relationship between the bone marrow niche and cardiometabolic diseases, highlighting how alterations in the niche contribute to disease development and progression. The niche regulates HSCs through complex interactions with stromal cells, endothelial cells, and signaling molecules. However, in the setting of chronic diseases such as hypertension, atherosclerosis, and diabetes, inflammatory signals disrupt the balance between HSC self-renewal and differentiation, promoting the excessive production of proinflammatory myeloid cells that exacerbate the disease. Key mechanisms discussed include the effects of hyperlipidemia, hyperglycemia, and sympathetic nervous system activation on HSC proliferation and differentiation. Furthermore, the review emphasizes the role of epigenetic modifications and metabolic reprogramming in creating trained immunity, a phenomenon whereby HSCs acquire long-term proinflammatory characteristics that sustain disease states. Finally, we explore therapeutic strategies aimed at targeting the bone marrow niche to mitigate chronic inflammation and its sequelae. Novel interventions that modulate hematopoiesis and restore niche homeostasis hold promise for the treatment of cardiometabolic diseases. By interrupting the vicious cycle of inflammation and marrow dysregulation, such therapies may offer new avenues for reducing cardiovascular risk and improving patient outcomes.

YY1 knockout in pro-B cells impairs lineage commitment, enabling unusual hematopoietic lineage plasticity

During B-cell development, cells progress through multiple developmental stages, with the pro-B-cell stage defining commitment to the B-cell lineage. YY1 is a ubiquitous transcription factor that is capable of both activation and repression functions. We found here that knockout of YY1 at the pro-B-cell stage eliminates B lineage commitment. YY1 knockout pro-B cells can generate T lineage cells in vitro using the OP9-DL4 feeder system and in vivo after injection into sublethally irradiated Rag1-/- mice. These T lineage-like cells lose their B lineage transcript profile and gain a T-cell lineage profile. Single-cell RNA-seq experiments showed that as YY1 knockout pro-B cells transition into T lineage cells in vitro, various cell clusters adopt transcript profiles representing a multiplicity of hematopoietic lineages, indicating unusual lineage plasticity. In addition, YY1 KO pro-B cells in vivo can give rise to other hematopoietic lineages in vivo. Evaluation of RNA-seq, scRNA-seq, ChIP-seq, and scATAC-seq data indicates that YY1 controls numerous chromatin-modifying proteins leading to increased accessibility of alternative lineage genes in YY1 knockout pro-B cells. Given the ubiquitous nature of YY1 and its dual activation and repression functions, YY1 may regulate commitment in multiple cell lineages.

In vivo CRISPR/Cas9-mediated screen reveals a critical function of TFDP1 and E2F4 transcription factors in hematopoiesis

Hematopoiesis is a continuous process of blood cell production driven by hematopoietic stem and progenitor cells (HSPCs) in the bone marrow. Proliferation and differentiation of HSPCs are regulated by complex transcriptional networks. In order to identify transcription factors with key roles in HSPC-mediated hematopoietic reconstitution, we developed an efficient and robust CRISPR/Cas9-based in vivo genetic screen. Using this experimental system, we identified the TFDP1 transcription factor to be essential for HSPC proliferation and post-transplant hematopoiesis. We further discovered that E2F4, an E2F transcription factor, serves as a binding partner of TFDP1 and is required for HSPC proliferation. Deletion of TFDP1 caused downregulation of genes associated with the cell cycle, with around 50% of these genes being identified as direct targets of TFDP1 and E2F4. Thus, our study expands the transcriptional network governing hematopoietic development through an in vivo CRISPR/Cas9-based genetic screen and identifies TFDP1/E2F4 as positive regulators of cell cycle genes in HSPCs.

Keratinocytes drive the epithelial hyperplasia key to sea lice resistance in coho salmon

BACKGROUND

Salmonid species have followed markedly divergent evolutionary trajectories in their interactions with sea lice. While sea lice parasitism poses significant economic, environmental, and animal welfare challenges for Atlantic salmon (Salmo salar) aquaculture, coho salmon (Oncorhynchus kisutch) exhibit near-complete resistance to sea lice, achieved through a potent epithelial hyperplasia response leading to rapid louse detachment. The molecular mechanisms underlying these divergent responses to sea lice are unknown.

RESULTS

We characterized the cellular and molecular responses of Atlantic salmon and coho salmon to sea lice using single-nuclei RNA sequencing. Juvenile fish were exposed to copepodid sea lice (Lepeophtheirus salmonis), and lice-attached pelvic fin and skin samples were collected 12 h, 24 h, 36 h, 48 h, and 60 h after exposure, along with control samples. Comparative analysis of control and treatment samples revealed an immune and wound-healing response that was common to both species, but attenuated in Atlantic salmon, potentially reflecting greater sea louse immunomodulation. Our results revealed unique but complementary roles of three layers of keratinocytes in the epithelial hyperplasia response leading to rapid sea lice rejection in coho salmon. Our results suggest that basal keratinocytes direct the expansion and mobility of intermediate and, especially, superficial keratinocytes, which eventually encapsulate the parasite.

CONCLUSIONS

Our results highlight the key role of keratinocytes in coho salmon's sea lice resistance and the diverged biological response of the two salmonid host species when interacting with this parasite. This study has identified key pathways and candidate genes that could be manipulated using various biotechnological solutions to improve Atlantic salmon sea lice resistance.

Unusual lineage plasticity revealed by YY1 knockout in pro-B cells

During B cell development, cells progress through multiple developmental stages with the pro-B cell stage defining commitment to the B cell lineage. YY1 is a ubiquitous transcription factor that is capable of both activation and repression functions. We find here that knockout of YY1 at the pro-B cell stage eliminates B lineage commitment. YY1 knockout pro-B cells can generate T lineage cells in vitro using the OP9- DL4 feeder system, as well as in vivo after injection into sub-lethally irradiated Rag1 -/- mice. These T lineage-like cells lose their B lineage transcript profile and gain a T cell lineage profile. Single cell-RNA-seq experiments showed that as YY1 knockout pro-B cells transition into T lineage cells, various cell clusters adopt transcript profiles representing a multiplicity of hematopoietic lineages indicating unusual lineage plasticity. Given the ubiquitous nature of YY1 and its dual activation and repression functions, YY1 likely regulates commitment in multiple cell lineages.

Progression of Notch signaling regulation of B cells under radiation exposure

With the continuous development of nuclear technology, the radiation exposure caused by radiation therapy is a serious health hazard. It is of great significance to further develop effective radiation countermeasures. B cells easily succumb to irradiation exposure along with immunosuppressive response. The approach to ameliorate radiation-induced B cell damage is rarely studied, implying that the underlying mechanisms of B cell damage after exposure are eager to be revealed. Recent studies suggest that Notch signaling plays an important role in B cell-mediated immune response. Notch signaling is a critical regulator for B cells to maintain immune function. Although accumulating studies reported that Notch signaling contributes to the functionality of hematopoietic stem cells and T cells, its role in B cells is scarcely appreciated. Presently, we discussed the regulation of Notch signaling on B cells under radiation exposure to provide a scientific basis to prevent radiation-induced B cell damage.

CaSSiDI: novel single-cell "Cluster Similarity Scoring and Distinction Index" reveals critical functions for PirB and context-dependent Cebpb repression

PirB is an inhibitory cell surface receptor particularly prominent on myeloid cells. PirB curtails the phenotypes of activated macrophages during inflammation or tumorigenesis, but its functions in macrophage homeostasis are obscure. To elucidate PirB-related functions in macrophages at steady-state, we generated and compared single-cell RNA-sequencing (scRNAseq) datasets obtained from myeloid cell subsets of wild type (WT) and PirB-deficient knockout (PirB KO) mice. To facilitate this analysis, we developed a novel approach to clustering parameter optimization called "Cluster Similarity Scoring and Distinction Index" (CaSSiDI). We demonstrate that CaSSiDI is an adaptable computational framework that facilitates tandem analysis of two scRNAseq datasets by optimizing clustering parameters. We further show that CaSSiDI offers more advantages than a standard Seurat analysis because it allows direct comparison of two or more independently clustered datasets, thereby alleviating the need for batch-correction while identifying the most similar and different clusters. Using CaSSiDI, we found that PirB is a novel regulator of Cebpb expression that controls the generation of Ly6Clo patrolling monocytes and the expansion properties of peritoneal macrophages. PirB's effect on Cebpb is tissue-specific since it was not observed in splenic red pulp macrophages (RPMs). However, CaSSiDI revealed a segregation of the WT RPM population into a CD68loIrf8+ "neuronal-primed" subset and an CD68hiFtl1+ "iron-loaded" subset. Our results establish the utility of CaSSiDI for single-cell assay analyses and the determination of optimal clustering parameters. Our application of CaSSiDI in this study has revealed previously unknown roles for PirB in myeloid cell populations. In particular, we have discovered homeostatic functions for PirB that are related to Cebpb expression in distinct macrophage subsets.

Differential carbonic anhydrase activities control EBV-induced B-cell transformation and lytic cycle reactivation

Epstein-Barr virus (EBV) contributes to ~1% of all human cancers including several B-cell neoplasms. A characteristic feature of EBV life cycle is its ability to transform metabolically quiescent B-lymphocytes into hyperproliferating B-cell blasts with the establishment of viral latency, while intermittent lytic cycle induction is necessary for the production of progeny virus. Our RNA-Seq analyses of both latently infected naïve B-lymphocytes and transformed B-lymphocytes upon lytic cycle replication indicate a contrasting expression pattern of a membrane-associated carbonic anhydrase isoform CA9, an essential component for maintaining cell acid-base homeostasis. We show that while CA9 expression is transcriptionally activated during latent infection model, lytic cycle replication restrains its expression. Pharmacological inhibition of CA-activity using specific inhibitors retards EBV induced B-cell transformation, inhibits B-cells outgrowth and colony formation ability of transformed B-lymphocytes through lowering the intracellular pH, induction of cell apoptosis and facilitating degradation of CA9 transcripts. Reanalyses of ChIP-Seq data along with utilization of EBNA2 knockout virus, ectopic expression of EBNA2 and sh-RNA mediated knockdown of CA9 expression we further demonstrate that EBNA2 mediated CA9 transcriptional activation is essential for EBV latently infected B-cell survival. In contrast, during lytic cycle reactivation CA9 expression is transcriptionally suppressed by the key EBV lytic cycle transactivator, BZLF1 through its transactivation domain. Overall, our study highlights the dynamic alterations of CA9 expression and its activity in regulating pH homeostasis act as one of the major drivers for EBV induced B-cell transformation and subsequent B-cell lymphomagenesis.

EBF1 Limits the Numbers of Cochlear Hair and Supporting Cells and Forms the Scala Tympani and Spiral Limbus during Inner Ear Development

Early B-cell factor 1 (EBF1) is a basic helix-loop-helix transcription factor essential for the differentiation of various tissues. Our single-cell RNA sequencing data suggest that Ebf1 is expressed in the sensory epithelium of the mouse inner ear. Here, we found that the murine Ebf1 gene and its protein are expressed in the prosensory domain of the inner ear, medial region of the cochlear duct floor, otic mesenchyme, and cochleovestibular ganglion. Ebf1 deletion in mice results in incomplete formation of the spiral limbus and scala tympani, increased number of cells in the organ of Corti and Kölliker's organ, and aberrant course of the spiral ganglion axons. Ebf1 deletion in the mouse cochlear epithelia caused the proliferation of SOX2-positive cochlear cells at E13.5, indicating that EBF1 suppresses the proliferation of the prosensory domain and cells of Kölliker's organ to facilitate the development of appropriate numbers of hair and supporting cells. Furthermore, mice with deletion of cochlear epithelium-specific Ebf1 showed poor postnatal hearing function. Our results suggest that Ebf1 is essential for normal auditory function in mammals.

Early B-Cell Factor 1: An Archetype for a Lineage-Restricted Transcription Factor Linking Development to Disease

The development of highly specialized blood cells from hematopoietic stem cells (HSCs) in the bone marrow (BM) is dependent upon a stringently orchestrated network of stage- and lineage-restricted transcription factors (TFs). Thus, the same stem cell can give rise to various types of differentiated blood cells. One of the key regulators of B-lymphocyte development is early B-cell factor 1 (EBF1). This TF belongs to a small, but evolutionary conserved, family of proteins that harbor a Zn-coordinating motif and an IPT/TIG (immunoglobulin-like, plexins, transcription factors/transcription factor immunoglobulin) domain, creating a unique DNA-binding domain (DBD). EBF proteins play critical roles in diverse developmental processes, including body segmentation in the Drosophila melanogaster embryo, and retina formation in mice. While several EBF family members are expressed in neuronal cells, adipocytes, and BM stroma cells, only B-lymphoid cells express EBF1. In the absence of EBF1, hematopoietic progenitor cells (HPCs) fail to activate the B-lineage program. This has been attributed to the ability of EBF1 to act as a pioneering factor with the ability to remodel chromatin, thereby creating a B-lymphoid-specific epigenetic landscape. Conditional inactivation of the Ebf1 gene in B-lineage cells has revealed additional functions of this protein in relation to the control of proliferation and apoptosis. This may explain why EBF1 is frequently targeted by mutations in human leukemia cases. This chapter provides an overview of the biochemical and functional properties of the EBF family proteins, with a focus on the roles of EBF1 in normal and malignant B-lymphocyte development.

Single-cell sequencing reveals the evolution of immune molecules across multiple vertebrate species

INTRODUCTION

Both innate and adaptive immune system undergo evolution from low to high vertebrates. Due to the limitation of conventional approaches in identifying broader spectrum of immune cells and molecules from various vertebrates, it remains unclear how immune molecules evolve among vertebrates.

OBJECTIVES

Here, we utilized carry out comparative transcriptome analysis in various immune cells across seven vertebrate species.

METHODS

Single-cell RNA sequencing (scRNA-seq).

RESULTS

We uncovered both conserved and species-specific profiling of gene expression in innate and adaptive immunity. Macrophages exhibited highly-diversified genes and developed sophisticated molecular signaling networks along with evolution, indicating effective and versatile functions in higher species. In contrast, B cells conservatively evolved with less differentially-expressed genes in analyzed species. Interestingly, T cells represented a dominant immune cell populations in all species and unique T cell populations were identified in zebrafish and pig. We also revealed compensatory TCR cascade components utilized by different species. Inter-species comparison of core gene programs demonstrated mouse species has the highest similarity in immune transcriptomes to human.

CONCLUSIONS

Therefore, our comparative study reveals gene transcription characteristics across multiple vertebrate species during the evolution of immune system, providing insights for species-specific immunity as well as the translation of animal studies to human physiology and disease.

The transcription factor EBF1 non-cell-autonomously regulates cardiac growth and differentiation

Reciprocal interactions between non-myocytes and cardiomyocytes regulate cardiac growth and differentiation. Here, we report that the transcription factor Ebf1 is highly expressed in non-myocytes and potently regulates heart development. Ebf1-deficient hearts display myocardial hypercellularity and reduced cardiomyocyte size, ventricular conduction system hypoplasia, and conduction system disease. Growth abnormalities in Ebf1 knockout hearts are observed as early as embryonic day 13.5. Transcriptional profiling of Ebf1-deficient embryonic cardiac non-myocytes demonstrates dysregulation of Polycomb repressive complex 2 targets, and ATAC-Seq reveals altered chromatin accessibility near many of these same genes. Gene set enrichment analysis of differentially expressed genes in cardiomyocytes isolated from E13.5 hearts of wild-type and mutant mice reveals significant enrichment of MYC targets and, consistent with this finding, we observe increased abundance of MYC in mutant hearts. EBF1-deficient non-myocytes, but not wild-type non-myocytes, are sufficient to induce excessive accumulation of MYC in co-cultured wild-type cardiomyocytes. Finally, we demonstrate that BMP signaling induces Ebf1 expression in embryonic heart cultures and controls a gene program enriched in EBF1 targets. These data reveal a previously unreported non-cell-autonomous pathway controlling cardiac growth and differentiation.

Chromatin accessibility profiling of targeted cell populations with laser capture microdissection coupled to ATAC-seq

Spatially resolved omics technologies reveal context-dependent cellular regulatory networks in tissues of interest. Beyond transcriptome analysis, information on epigenetic traits and chromatin accessibility can provide further insights on gene regulation in health and disease. Nevertheless, compared to the enormous advancements in spatial transcriptomics technologies, the field of spatial epigenomics is much younger and still underexplored. In this study, we report laser capture microdissection coupled to ATAC-seq (LCM-ATAC-seq) applied to fresh frozen samples for the spatial characterization of chromatin accessibility. We first demonstrate the efficient use of LCM coupled to in situ tagmentation and evaluate its performance as a function of cell number, microdissected areas, and tissue type. Further, we demonstrate its use for the targeted chromatin accessibility analysis of discrete contiguous or scattered cell populations in tissues via single-nuclei capture based on immunostaining for specific cellular markers.

Pre-pregnancy gene expression signatures are associated with subsequent improvement/worsening of rheumatoid arthritis during pregnancy

BACKGROUND

While many women with rheumatoid arthritis (RA) improve during pregnancy and others worsen, there are no biomarkers to predict this improvement or worsening. In our unique RA pregnancy cohort that includes a pre-pregnancy baseline, we have examined pre-pregnancy gene co-expression networks to identify differences between women with RA who subsequently improve during pregnancy and those who worsen.

METHODS

Blood samples were collected before pregnancy (T0) from 19 women with RA and 13 healthy women enrolled in our prospective pregnancy cohort. RA improvement/worsening between T0 and 3rd trimester was assessed by changes in the Clinical Disease Activity Index (CDAI). Pre-pregnancy expression profiles were examined by RNA sequencing and differential gene expression analysis. Weighted gene co-expression network analysis (WGCNA) was used to identify co-expression modules correlated with the improvement/worsening of RA during pregnancy and to assess their functional relevance.

RESULTS

Of the 19 women with RA, 14 improved during pregnancy (RAimproved) while 5 worsened (RAworsened). At the T0 baseline, however, the mean CDAI was similar between the two groups. WGCNA identified one co-expression module related to B cell function that was significantly correlated with the worsening of RA during pregnancy and was significantly enriched in genes differentially expressed between the RAimproved and RAworsened groups. A neutrophil-related expression signature was also identified in the RAimproved group at the T0 baseline.

CONCLUSION

The pre-pregnancy gene expression signatures identified represent potential biomarkers to predict the subsequent improvement/worsening of RA during pregnancy, which has important implications for the personalized treatment of RA during pregnancy.

CEBPA phase separation links transcriptional activity and 3D chromatin hubs

Cell identity is orchestrated through an interplay between transcription factor (TF) action and genome architecture. The mechanisms used by TFs to shape three-dimensional (3D) genome organization remain incompletely understood. Here we present evidence that the lineage-instructive TF CEBPA drives extensive chromatin compartment switching and promotes the formation of long-range chromatin hubs during induced B cell-to-macrophage transdifferentiation. Mechanistically, we find that the intrinsically disordered region (IDR) of CEBPA undergoes in vitro phase separation (PS) dependent on aromatic residues. Both overexpressing B cells and native CEBPA-expressing cell types such as primary granulocyte-macrophage progenitors, liver cells, and trophectoderm cells reveal nuclear CEBPA foci and long-range 3D chromatin hubs at CEBPA-bound regions. In short, we show that CEBPA can undergo PS through its IDR, which may underlie in vivo foci formation and suggest a potential role of PS in regulating CEBPA function.

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.

ExplaiNN: interpretable and transparent neural networks for genomics

Deep learning models such as convolutional neural networks (CNNs) excel in genomic tasks but lack interpretability. We introduce ExplaiNN, which combines the expressiveness of CNNs with the interpretability of linear models. ExplaiNN can predict TF binding, chromatin accessibility, and de novo motifs, achieving performance comparable to state-of-the-art methods. Its predictions are transparent, providing global (cell state level) as well as local (individual sequence level) biological insights into the data. ExplaiNN can serve as a plug-and-play platform for pretrained models and annotated position weight matrices. ExplaiNN aims to accelerate the adoption of deep learning in genomic sequence analysis by domain experts.

A single-cell transcriptional timelapse of mouse embryonic development, from gastrula to pup

The house mouse, Mus musculus, is an exceptional model system, combining genetic tractability with close homology to human biology. Gestation in mouse development lasts just under three weeks, a period during which its genome orchestrates the astonishing transformation of a single cell zygote into a free-living pup composed of >500 million cells. Towards a global framework for exploring mammalian development, we applied single cell combinatorial indexing (sci-*) to profile the transcriptional states of 12.4 million nuclei from 83 precisely staged embryos spanning late gastrulation (embryonic day 8 or E8) to birth (postnatal day 0 or P0), with 2-hr temporal resolution during somitogenesis, 6-hr resolution through to birth, and 20-min resolution during the immediate postpartum period. From these data (E8 to P0), we annotate dozens of trajectories and hundreds of cell types and perform deeper analyses of the unfolding of the posterior embryo during somitogenesis as well as the ontogenesis of the kidney, mesenchyme, retina, and early neurons. Finally, we leverage the depth and temporal resolution of these whole embryo snapshots, together with other published data, to construct and curate a rooted tree of cell type relationships that spans mouse development from zygote to pup. Throughout this tree, we systematically nominate sets of transcription factors (TFs) and other genes as candidate drivers of the in vivo differentiation of hundreds of mammalian cell types. Remarkably, the most dramatic shifts in transcriptional state are observed in a restricted set of cell types in the hours immediately following birth, and presumably underlie the massive changes in physiology that must accompany the successful transition of a placental mammal to extrauterine life.

Hematopoietic Stem Cells and the Immune System in Development and Aging

Hematopoietic stem cells (HSCs) support haematopoiesis throughout life and give rise to the whole variety of cells of the immune system. Developing in the early embryo, passing through the precursor stage, and maturing into the first HSCs, they undergo a fairly large number of divisions while maintaining a high regenerative potential due to high repair activity. This potential is greatly reduced in adult HSCs. They go into a state of dormancy and anaerobic metabolism to maintain their stemness throughout life. However, with age, changes occur in the pool of HSCs that negatively affect haematopoiesis and the effectiveness of immunity. Niche aging and accumulation of mutations with age reduces the ability of HSCs to self-renew and changes their differentiation potential. This is accompanied by a decrease in clonal diversity and a disturbance of lymphopoiesis (decrease in the formation of naive T- and B-cells) and the predominance of myeloid haematopoiesis. Aging also affects mature cells, regardless of HSC, therefore, phagocytic activity and the intensity of the oxidative burst decrease, and the efficiency of processing and presentation of antigens by myeloid cells is impaired. Aging cells of innate and adaptive immunity produce factors that form a chronic inflammatory background. All these processes have a serious negative impact on the protective properties of the immune system, increasing inflammation, the risk of developing autoimmune, oncological, and cardiovascular diseases with age. Understanding the mechanisms of reducing the regenerative potential in a comparative analysis of embryonic and aging HSCs, the features of inflammatory aging will allow us to get closer to deciphering the programs for the development, aging, regeneration and rejuvenation of HSCs and the immune system.

miRNAs as cornerstones in chronic lymphocytic leukemia pathogenesis and therapeutic resistance- An emphasis on the interaction of signaling pathways

Chronic lymphocytic leukemia (CLL) accounts for the vast majority of cases of leukemia. Patients of advanced age are more likely to develop the condition, which has a highly varied clinical course. Consideration of illness features and preceding treatment sequence, as well as patient preferences and comorbidities, is necessary for selecting the appropriate treatment for the appropriate patient. Therefore, there is an urgent need for novel biomarkers with high sensitivity and specificity to detect CLL early, monitor CLL patients, select the treatment responders, and reduce ineffective treatment, unwanted side effects, and unnecessary expenses. In both homeostasis and illness, microRNAs (miRNAs/miRs) play a vital role as master regulators of gene expression and, by extension, protein expression. MiRNAs typically reduce the stability of mRNAs, including those encoding genes involved in tumorigenesis processes as cell cycle regulation, inflammation, stress response, angiogenesis, differentiation, apoptosis, and invasion. Due to their unique properties, miRNAs are rapidly being exploited as accurate biomarkers for illness detection, and medicines based on miRNA targets are finding widespread application in clinical practice. Accordingly, the current review serves as a quick primer on CLL and the biogenesis of miRNAs. In addition to providing a brief overview of the miRNAs whose function in the progression of CLL has been established by recent in vitro or in vivo research through articulating the influence of these miRNAs on a wide variety of cellular functions, including increased proliferative potential; support for angiogenesis; cell cycle aberration; evasion of apoptosis; promotion of metastasis; and reduced sensitivity to specific treatments.

IDH2 and TET2 mutations synergize to modulate T Follicular Helper cell functional interaction with the AITL microenvironment

Angioimmunoblastic T cell lymphoma (AITL) is a peripheral T cell lymphoma that originates from T follicular helper (Tfh) cells and exhibits a prominent tumor microenvironment (TME). IDH2 and TET2 mutations co-occur frequently in AITL, but their contribution to tumorigenesis is poorly understood. We developed an AITL mouse model that is driven by Idh2 and Tet2 mutations. Malignant Tfh cells display aberrant transcriptomic and epigenetic programs that impair TCR signaling. Neoplastic Tfh cells bearing combined Idh2 and Tet2 mutations show altered cross-talk with germinal center B cells that promotes B cell clonal expansion while decreasing Fas-FasL interaction and reducing B cell apoptosis. The plasma cell count and angiogenesis are also increased in the Idh2-mutated tumors, implying a major relationship between Idh2 mutation and the characteristic AITL TME. Our mouse model recapitulates several features of human IDH2-mutated AITL and provides a rationale for exploring therapeutic targeting of Tfh-TME cross-talk for AITL patients.

HyperIgE in hypomorphic recombination-activating gene defects

Increased immunogloblulin-E (IgE) levels associated with eosinophilia represent a common finding observed in Omenn syndrome, a severe immunodeficiency caused by decreased V(D)J recombination, leading to restricted T- and B-cell receptor repertoire. V(D)J recombination is initiated by the lymphoid-restricted recombination-activating gene (RAG) recombinases. The lack of RAG proteins causes a block in lymphocyte differentiation, resulting in T^-B^- severe combined immunodeficiency. Conversely, hypomorphic mutations allow the generation of few T and B cells, leading to a spectrum of immunological phenotypes, in which immunodeficiency associates to inflammation, immune dysregulation, and autoimmunity. Elevated IgE levels are frequently observed in hypomorphic RAG patients. Here, we describe the role of RAG genes in lymphocyte differentiation and maintenance of immune tolerance.

Cholinergic control of Th17 cell pathogenicity in experimental autoimmune encephalomyelitis

Experimental autoimmune encephalomyelitis (EAE) is a mouse model of multiple sclerosis (MS) in which Th17 cells have a crucial but unclear function. Here we show that choline acetyltransferase (ChAT), which synthesizes acetylcholine (ACh), is a critical driver of pathogenicity in EAE. Mice with ChAT-deficient Th17 cells resist disease progression and show reduced brain-infiltrating immune cells. ChAT expression in Th17 cells is linked to strong TCR signaling, expression of the transcription factor Bhlhe40, and increased Il2, Il17, Il22, and Il23r mRNA levels. ChAT expression in Th17 cells is independent of IL21r signaling but dampened by TGFβ, implicating ChAT in controlling the dichotomous nature of Th17 cells. Our study establishes a cholinergic program in which ACh signaling primes chronic activation of Th17 cells, and thereby constitutes a pathogenic determinant of EAE. Our work may point to novel targets for therapeutic immunomodulation in MS.

Indirect Mechanisms of Transcription Factor-Mediated Gene Regulation during Cell Fate Changes

Transcription factors (TFs) are the master regulators of cellular identity, capable of driving cell fate transitions including differentiations, reprogramming, and transdifferentiations. Pioneer TFs recognize partial motifs exposed on nucleosomal DNA, allowing for TF-mediated activation of repressed chromatin. Moreover, there is evidence suggesting that certain TFs can repress actively expressed genes either directly through interactions with accessible regulatory elements or indirectly through mechanisms that impact the expression, activity, or localization of other regulatory factors. Recent evidence suggests that during reprogramming, the reprogramming TFs initiate opening of chromatin regions rich in somatic TF motifs that are inaccessible in the initial and final cellular states. It is postulated that analogous to a sponge, these transiently accessible regions "soak up" somatic TFs, hence lowering the initial barriers to cell fate changes. This indirect TF-mediated gene regulation event, which is aptly named the "sponge effect," may play an essential role in the silencing of the somatic transcriptional network during different cellular conversions.

Maintenance of Lineage Identity: Lessons from a B Cell

The maintenance of B cell identity requires active transcriptional control that enforces a B cell-specific program and suppresses alternative lineage genes. Accordingly, disrupting the B cell identity regulatory network compromises B cell function and induces cell fate plasticity by allowing derepression of alternative lineage-specific transcriptional programs. Although the B lineage is incredibly resistant to most differentiating factors, loss of just a single B lineage-specific transcription factor or the forced expression of individual non-B cell lineage transcription factors can radically disrupt B cell maintenance and allow dedifferentiation or transdifferentiation into entirely distinct lineages. B lymphocytes thereby offer an insightful and useful case study of how a specific cell lineage can maintain a stable identity throughout life and how perturbations of a single master regulator can induce cellular plasticity. In this article, we review the regulatory mechanisms that safeguard B cell identity, and we discuss how dysregulation of the B cell maintenance program can drive malignant transformation and enable therapeutic resistance.

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.

Ikaros Regulates microRNA Networks in Acute Lymphoblastic Leukemia

The hematopoietic transcription factor Ikaros (IKZF1) regulates normal B cell development and functions as a tumor suppressor in precursor B cell acute lymphoblastic leukemia (B-ALL). MicroRNAs (miRNAs) are small regulatory RNAs that through post-transcriptional gene regulation play critical roles in intracellular processes including cell growth in cancer. However, the role of Ikaros in the regulation of miRNA expression in developing B cells is unknown. In this study, we examined the Ikaros-regulated miRNA targets using human IKZF1-mutated Ph+ B-ALL cell lines. Inducible expression of wild-type Ikaros (the Ik1 isoform) caused B-ALL growth arrest and exit from the cell cycle. Global miRNA expression analysis revealed a total of 31 miRNAs regulated by IK1, and ChIP-seq analysis showed that Ikaros bound to several Ik1-responsive miRNA genes. Examination of the prognostic significance of miRNA expression in B-ALL indicate that the IK1-regulated miRNAs hsa-miR-26b, hsa-miR-130b and hsa-miR-4649 are significantly associated with outcome in B-ALL. Our findings establish a potential regulatory circuit between the tumor-suppressor Ikaros and the oncogenic miRNA networks in IKZF1-mutated B-ALL. These results indicate that Ikaros regulates the expression of a subset of miRNAs, of which several may contribute to B-ALL growth.

Systems-level identification of key transcription factors in immune cell specification

Transcription factors (TFs) are crucial for regulating cell differentiation during the development of the immune system. However, the key TFs for orchestrating the specification of distinct immune cells are not fully understood. Here, we integrated the transcriptomic and epigenomic measurements in 73 mouse and 61 human primary cell types, respectively, that span the immune cell differentiation pathways. We constructed the cell-type-specific transcriptional regulatory network and assessed the global importance of TFs based on the Taiji framework, which is a method we have previously developed that can infer the global impact of TFs using integrated transcriptomic and epigenetic data. Integrative analysis across cell types revealed putative driver TFs in cell lineage-specific differentiation in both mouse and human systems. We have also identified TF combinations that play important roles in specific developmental stages. Furthermore, we validated the functions of predicted novel TFs in murine CD8+ T cell differentiation and showed the importance of Elf1 and Prdm9 in the effector versus memory T cell fate specification and Kdm2b and Tet3 in promoting differentiation of CD8+ tissue resident memory (Trm) cells, validating the approach. Thus, we have developed a bioinformatic approach that provides a global picture of the regulatory mechanisms that govern cellular differentiation in the immune system and aids the discovery of novel mechanisms in cell fate decisions.

Assessment of DNA methylation in porcine immune cells reveals novel regulatory elements associated with cell-specific gene expression and immune capacity traits

Background

Genetics studies in the porcine immune system have enhanced selection practices for disease resistance phenotypes and increased the efficacy of porcine models in biomedical research; however limited functional annotation of the porcine immunome has hindered progress on both fronts. Among epigenetic mechanisms that regulate gene expression, DNA methylation is the most ubiquitous modification made to the DNA molecule and influences transcription factor binding as well as gene and phenotype expression. Human and mouse DNA methylation studies have improved mapping of regulatory elements in these species, but comparable studies in the pig have been limited in scope.

Results

We performed whole-genome bisulfite sequencing to assess DNA methylation patterns in nine pig immune cell populations: CD21⁺ and CD21⁻ B cells, four T cell fractions (CD4⁺, CD8⁺, CD8⁺CD4⁺, and SWC6γδ⁺), natural killer and myeloid cells, and neutrophils. We identified 54,391 cell differentially methylated regions (cDMRs), and clustering by cDMR methylation rate grouped samples by cell lineage. 32,737 cDMRs were classified as cell lowly methylated regions (cLMRs) in at least one cell type, and cLMRs were broadly enriched in genes and regions of intermediate CpG density. We observed strong correlations between differential methylation and expression across immune cell populations, with cell-specific low methylation disproportionately impacting genes exhibiting enriched gene expression in the same cell type. Motif analysis of cLMRs revealed cell type-specific enrichment of transcription factor binding motifs, indicating that cell-specific methylation patterns may influence accessibility by trans-acting factors. Lastly, cDMRs were enriched for immune capacity GWAS SNPs, and many such overlaps occurred within genes known to influence immune cell development and function (CD8B, NDRG1).

Conclusion

Our DNA methylation data improve functional annotation of the porcine genome through characterization of epigenomic regulatory patterns that contribute to immune cell identity and function, and increase the potential for identifying mechanistic links between genotype and phenotype.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08773-5.

Tnpo3 enables EBF1 function in conditions of antagonistic Notch signaling

Transcription factor EBF1 (early B cell factor 1) acts as a key regulator of B cell specification. The transcriptional network in which EBF1 operates has been extensively studied; however, the regulation of EBF1 function remains poorly defined. By mass spectrometric analysis of proteins associated with endogenous EBF1 in pro-B cells, we identified the nuclear import receptor Transportin-3 (Tnpo3) and found that it interacts with the immunoglobulin-like fold domain of EBF1. We delineated glutamic acid 271 of EBF1 as a critical residue for the association with Tnpo3. EBF1E271A showed normal nuclear localization; however, it had an impaired B cell programming ability in conditions of Notch signaling, as determined by retroviral transduction of Ebf1 -/- progenitors. By RNA-seq analysis of EBF1E271A-expressing progenitors, we found an up-regulation of T lineage determinants and down-regulation of early B genes, although similar chromatin binding of EBF1E271A and EBF1wt was detected in pro-B cells expressing activated Notch1. B lineage-specific inactivation of Tnpo3 in mice resulted in a block of early B cell differentiation, accompanied by a down-regulation of B lineage genes and up-regulation of T and NK lineage genes. Taken together, our observations suggest that Tnpo3 ensures B cell programming by EBF1 in nonpermissive conditions.

Magnitude and breadth of antibody cross-reactivity induced by recombinant influenza hemagglutinin trimer vaccine is enhanced by combination adjuvants

The effects of adjuvants for increasing the immunogenicity of influenza vaccines are well known. However, the effect of adjuvants on increasing the breadth of cross-reactivity is less well understood. In this study we have performed a systematic screen of different toll-like receptor (TLR) agonists, with and without a squalene-in-water emulsion on the immunogenicity of a recombinant trimerized hemagglutinin (HA) vaccine in mice after single-dose administration. Antibody (Ab) cross-reactivity for other variants within and outside the immunizing subtype (homosubtypic and heterosubtypic cross-reactivity, respectively) was assessed using a protein microarray approach. Most adjuvants induced broad IgG profiles, although the response to a combination of CpG, MPLA and AddaVax (termed 'IVAX-1') appeared more quickly and reached a greater magnitude than the other formulations tested. Antigen-specific plasma cell labeling experiments show the components of IVAX-1 are synergistic. This adjuvant preferentially stimulates CD4 T cells to produce Th1>Th2 type (IgG2c>IgG1) antibodies and cytokine responses. Moreover, IVAX-1 induces identical homo- and heterosubtypic IgG and IgA cross-reactivity profiles when administered intranasally. Consistent with these observations, a single-cell transcriptomics analysis demonstrated significant increases in expression of IgG1, IgG2b and IgG2c genes of B cells in H5/IVAX-1 immunized mice relative to naïve mice, as well as significant increases in expression of the IFNγ gene of both CD4 and CD8 T cells. These data support the use of adjuvants for enhancing the breath and durability of antibody responses of influenza virus vaccines.

FOXO Dictates Initiation of B Cell Development and Myeloid Restriction in Common Lymphoid Progenitors

The development of B cells relies on an intricate network of transcription factors critical for developmental progression and lineage commitment. In the B cell developmental trajectory, a temporal switch from predominant Foxo3 to Foxo1 expression occurs at the CLP stage. Utilizing VAV-iCre mediated conditional deletion, we found that the loss of FOXO3 impaired B cell development from LMPP down to B cell precursors, while the loss of FOXO1 impaired B cell commitment and resulted in a complete developmental block at the CD25 negative proB cell stage. Strikingly, the combined loss of FOXO1 and FOXO3 resulted in the failure to restrict the myeloid potential of CLPs and the complete loss of the B cell lineage. This is underpinned by the failure to enforce the early B-lineage gene regulatory circuitry upon a predominantly pre-established open chromatin landscape. Altogether, this demonstrates that FOXO3 and FOXO1 cooperatively govern early lineage restriction and initiation of B-lineage commitment in CLPs.

Lamin B1 deletion in myeloid neoplasms causes nuclear anomaly and altered hematopoietic stem cell function

Abnormal nuclear morphology is a hallmark of malignant cells widely used in cancer diagnosis. Pelger-Huët anomaly (PHA) is a common abnormality of neutrophil nuclear morphology of unknown molecular etiology in myeloid neoplasms (MNs). We show that loss of nuclear lamin B1 (LMNB1) encoded on chromosome 5q, which is frequently deleted in MNs, induces defects in nuclear morphology and human hematopoietic stem cell (HSC) function associated with malignancy. LMNB1 deficiency alters genome organization inducing in vitro and in vivo expansion of HSCs, myeloid-biased differentiation with impaired lymphoid commitment, and genome instability due to defective DNA damage repair. Nuclear dysmorphology of neutrophils in patients with MNs is associated with 5q deletions spanning the LMNB1 locus, and lamin B1 loss is both necessary and sufficient to cause PHA in normal and 5q-deleted neutrophils. LMNB1 loss thus causes acquired PHA and links abnormal nuclear morphology with HSCs and progenitor cell fate determination via genome organization.

Integrating transcription-factor abundance with chromatin accessibility in human erythroid lineage commitment

Master transcription factors (TFs) directly regulate present and future cell states by binding DNA regulatory elements and driving gene-expression programs. Their abundance influences epigenetic priming to different cell fates at the chromatin level, especially in the context of differentiation. In order to link TF protein abundance to changes in TF motif accessibility and open chromatin, we developed InTAC-seq, a method for simultaneous quantification of genome-wide chromatin accessibility and intracellular protein abundance in fixed cells. Our method produces high-quality data and is a cost-effective alternative to single-cell techniques. We showcase our method by purifying bone marrow (BM) progenitor cells based on GATA-1 protein levels and establish high GATA-1-expressing BM cells as both epigenetically and functionally similar to erythroid-committed progenitors.

EBF1 nuclear repositioning instructs chromatin refolding to promote therapy resistance in T leukemic cells

Chromatin misfolding has been implicated in cancer pathogenesis; yet, its role in therapy resistance remains unclear. Here, we systematically integrated sequencing and imaging data to examine the spatial and linear chromatin structures in targeted therapy-sensitive and -resistant human T cell acute lymphoblastic leukemia (T-ALL). We found widespread alterations in successive layers of chromatin organization including spatial compartments, contact domain boundaries, and enhancer positioning upon the emergence of targeted therapy resistance. The reorganization of genome folding structures closely coincides with the restructuring of chromatin activity and redistribution of architectural proteins. Mechanistically, the derepression and repositioning of the B-lineage-determining transcription factor EBF1 from the heterochromatic nuclear envelope to the euchromatic interior instructs widespread genome refolding and promotes therapy resistance in leukemic T cells. Together, our findings suggest that lineage-determining transcription factors can instruct changes in genome topology as a driving force for epigenetic adaptations in targeted therapy resistance.

Noncoding RNAs and Their Response Predictive Value in Azacitidine-treated Patients With Myelodysplastic Syndrome and Acute Myeloid Leukemia With Myelodysplasia-related Changes

BACKGROUND/AIM

Prediction of response to azacitidine (AZA) treatment is an important challenge in hematooncology. In addition to protein coding genes (PCGs), AZA efficiency is influenced by various noncoding RNAs (ncRNAs), including long ncRNAs (lncRNAs), circular RNAs (circRNAs), and transposable elements (TEs).

MATERIALS AND METHODS

RNA sequencing was performed in patients with myelodysplastic syndromes or acute myeloid leukemia before AZA treatment to assess contribution of ncRNAs to AZA mechanisms and propose novel disease prediction biomarkers.

RESULTS

Our analyses showed that lncRNAs had the strongest predictive potential. The combined set of the best predictors included 14 lncRNAs, and only four PCGs, one circRNA, and no TEs. Epigenetic regulation and recombinational repair were suggested as crucial for AZA response, and network modeling defined three deregulated lncRNAs (CTC-482H14.5, RP11-419K12.2, and RP11-736I24.4) associated with these processes.

CONCLUSION

The expression of various ncRNAs can influence the effect of AZA and new ncRNA-based predictive biomarkers can be defined.

CT-FOCS: a novel method for inferring cell type-specific enhancer–promoter maps

Spatiotemporal gene expression patterns are governed to a large extent by the activity of enhancer elements, which engage in physical contacts with their target genes. Identification of enhancer–promoter (EP) links that are functional only in a specific subset of cell types is a key challenge in understanding gene regulation. We introduce CT-FOCS (cell type FOCS), a statistical inference method that uses linear mixed effect models to infer EP links that show marked activity only in a single or a small subset of cell types out of a large panel of probed cell types. Analyzing 808 samples from FANTOM5, covering 472 cell lines, primary cells and tissues, CT-FOCS inferred such EP links more accurately than recent state-of-the-art methods. Furthermore, we show that strictly cell type-specific EP links are very uncommon in the human genome.

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.

A regulatory network of microRNAs confers lineage commitment during early developmental trajectories of B and T lymphocytes

The commitment of hematopoietic multipotent progenitors (MPPs) toward a particular lineage involves activation of cell type-specific genes and silencing of genes that promote alternate cell fates. Although the gene expression programs of early-B and early-T lymphocyte development are mutually exclusive, we show that these cell types exhibit significantly correlated microRNA (miRNA) profiles. However, their corresponding miRNA targetomes are distinct and predominated by transcripts associated with natural killer, dendritic cell, and myeloid lineages, suggesting that miRNAs function in a cell-autonomous manner. The combinatorial expression of miRNAs miR-186-5p, miR-128-3p, and miR-330-5p in MPPs significantly attenuates their myeloid differentiation potential due to repression of myeloid-associated transcripts. Depletion of these miRNAs caused a pronounced de-repression of myeloid lineage targets in differentiating early-B and early-T cells, resulting in a mixed-lineage gene expression pattern. De novo motif analysis combined with an assay of promoter activities indicates that B as well as T lineage determinants drive the expression of these miRNAs in lymphoid lineages. Collectively, we present a paradigm that miRNAs are conserved between developing B and T lymphocytes, yet they target distinct sets of promiscuously expressed lineage-inappropriate genes to suppress the alternate cell-fate options. Thus, our studies provide a comprehensive compendium of miRNAs with functional implications for B and T lymphocyte development.

Tree-Based Co-Clustering Identifies Chromatin Accessibility Patterns Associated With Hematopoietic Lineage Structure

Chromatin accessibility, as measured by ATACseq, varies between hematopoietic cell types in different lineages of the hematopoietic differentiation tree, e.g. T cells vs. B cells, but methods that associate variation in chromatin accessibility to the lineage structure of the differentiation tree are lacking. Using an ATACseq dataset recently published by the ImmGen consortium, we construct associations between chromatin accessibility and hematopoietic cell types using a novel co-clustering approach that accounts for the structure of the hematopoietic, differentiation tree. Under a model in which all loci and cell types within a co-cluster have a shared accessibility state, we show that roughly 80% of cell type associated accessibility variation can be captured through 12 cell type clusters and 20 genomic locus clusters, with the cell type clusters reflecting coherent components of the differentiation tree. Using publicly available ChIPseq datasets, we show that our clustering reflects transcription factor binding patterns with implications for regulation across cell types. We show that traditional methods such as hierarchical and kmeans clusterings lead to cell type clusters that are more dispersed on the tree than our tree-based algorithm. We provide a python package, chromcocluster, that implements the algorithms presented.

The transcriptomic blueprint of molt in rooster using various tissues from Ginkkoridak (Korean long-tailed chicken)

BACKGROUND

Annual molt is a critical stage in the life cycle of birds. Although the most extensively documented aspects of molt are the renewing of plumage and the remodeling of the reproductive tract in laying hens, in chicken, molt deeply affects various tissues and physiological functions. However, with exception of the reproductive tract, the effect of molt on gene expression across the tissues known to be affected by molt has to date never been investigated. The present study aimed to decipher the transcriptomic effects of molt in Ginkkoridak, a Korean long-tailed chicken. Messenger RNA data available across 24 types of tissue samples (9 males) and a combination of mRNA and miRNA data on 10 males and 10 females blood were used.

RESULTS

The impact of molt on gene expression and gene transcript usage appeared to vary substantially across tissues types in terms of histological entities or physiological functions particularly related to nervous system. Blood was the tissue most affected by molt in terms of differentially expressed genes in both sexes, closely followed by meninges, bone marrow and heart. The effect of molt in blood appeared to differ between males and females, with a more than fivefold difference in the number of down-regulated genes between both sexes. The blueprint of molt in roosters appeared to be specific to tissues or group of tissues, with relatively few genes replicating extensively across tissues, excepted for the spliceosome genes (U1, U4) and the ribosomal proteins (RPL21, RPL23). By integrating miRNA and mRNA data, when chickens molt, potential roles of miRNA were discovered such as regulation of neurogenesis, regulation of immunity and development of various organs. Furthermore, reliable candidate biomarkers of molt were found, which are related to cell dynamics, nervous system or immunity, processes or functions that have been shown to be extensively modulated in response to molt.

CONCLUSIONS

Our results provide a comprehensive description at the scale of the whole organism deciphering the effects of molt on the transcriptome in chicken. Also, the conclusion of this study can be used as a valuable resource in transcriptome analyses of chicken in the future and provide new insights related to molt.

Rituximab-resistant splenic memory B cells and newly engaged naive B cells fuel relapses in patients with immune thrombocytopenia

Rituximab (RTX), an antibody targeting CD20, is widely used as a first-line therapeutic strategy in B cell-mediated autoimmune diseases. However, a large proportion of patients either do not respond to the treatment or relapse during B cell reconstitution. Here, we characterize the cellular basis responsible for disease relapse in secondary lymphoid organs in humans, taking advantage of the opportunity offered by therapeutic splenectomy in patients with relapsing immune thrombocytopenia. By analyzing the B and plasma cell immunoglobulin gene repertoire at bulk and antigen-specific single-cell level, we demonstrate that relapses are associated with two responses coexisting in germinal centers and involving preexisting mutated memory B cells that survived RTX treatment and naive B cells generated upon reconstitution of the B cell compartment. To identify distinctive characteristics of the memory B cells that escaped RTX-mediated depletion, we analyzed RTX refractory patients who did not respond to treatment at the time of B cell depletion. We identified, by single-cell RNA sequencing (scRNA-seq) analysis, a population of quiescent splenic memory B cells that present a unique, yet reversible, RTX-shaped phenotype characterized by down-modulation of B cell-specific factors and expression of prosurvival genes. Our results clearly demonstrate that these RTX-resistant autoreactive memory B cells reactivate as RTX is cleared and give rise to plasma cells and further germinal center reactions. Their continued surface expression of CD19 makes them efficient targets for current anti-CD19 therapies. This study thus identifies a pathogenic contributor to autoimmune diseases that can be targeted by available therapeutic agents.

ACE2-lentiviral transduction enables mouse SARS-CoV-2 infection and mapping of receptor interactions

SARS-CoV-2 uses the human ACE2 (hACE2) receptor for cell attachment and entry, with mouse ACE2 (mACE2) unable to support infection. Herein we describe an ACE2-lentivirus system and illustrate its utility for in vitro and in vivo SARS-CoV-2 infection models. Transduction of non-permissive cell lines with hACE2 imparted replication competence, and transduction with mACE2 containing N30D, N31K, F83Y and H353K substitutions, to match hACE2, rescued SARS-CoV-2 replication. Intrapulmonary hACE2-lentivirus transduction of C57BL/6J mice permitted significant virus replication in lung epithelium. RNA-Seq and histological analyses illustrated that this model involved an acute inflammatory disease followed by resolution and tissue repair, with a transcriptomic profile similar to that seen in COVID-19 patients. hACE2-lentivirus transduction of IFNAR-/- and IL-28RA-/- mouse lungs was used to illustrate that loss of type I or III interferon responses have no significant effect on virus replication. However, their importance in driving inflammatory responses was illustrated by RNA-Seq analyses. We also demonstrate the utility of the hACE2-lentivirus transduction system for vaccine evaluation in C57BL/6J mice. The ACE2-lentivirus system thus has broad application in SARS-CoV-2 research, providing a tool for both mutagenesis studies and mouse model development.

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

Genes encoding B lineage-restricted transcription factors are frequently mutated in B-lymphoid leukemias, suggesting a close link between normal and malignant B-cell development. One of these transcription factors is early B-cell factor 1 (EBF1), a protein of critical importance for lineage specification and survival of B-lymphoid progenitors. Here, we report that impaired EBF1 function in mouse B-cell progenitors results in reduced expression of Myc. Ectopic expression of MYC partially rescued B-cell expansion in the absence of EBF1 both in vivo and in vitro. Using chromosome conformation analysis in combination with ATAC-sequencing, chromatin immunoprecipitation-sequencing, and reporter gene assays, six EBF1-responsive enhancer elements were identified within the Myc locus. CRISPR-Cas9-mediated targeting of EBF1-binding sites identified one element of key importance for Myc expression and pro-B cell expansion. These data provide evidence that Myc is a direct target of EBF1. Furthermore, chromatin immunoprecipitation-sequencing analysis revealed that several regulatory elements in the Myc locus are targets of PAX5. However, ectopic expression of PAX5 in EBF1-deficient cells inhibits the cell cycle and reduces Myc expression, suggesting that EBF1 and PAX5 act in an opposing manner to regulate Myc levels. This hypothesis is further substantiated by the finding that Pax5 inactivation reduces requirements for EBF1 in pro-B-cell expansion. The binding of EBF1 and PAX5 to regulatory elements in the human MYC gene in a B-cell acute lymphoblastic leukemia cell line indicates that the EBF1:PAX5:MYC regulatory loop is conserved and may control both normal and malignant B-cell development.

The Evolution of Rag Gene Enhancers and Transcription Factor E and Id Proteins in the Adaptive Immune System

Adaptive immunity relies on the V(D)J DNA recombination of immunoglobulin (Ig) and T cell receptor (TCR) genes, which enables the recognition of highly diverse antigens and the elicitation of antigen-specific immune responses. This process is mediated by recombination-activating gene (Rag) 1 and Rag2 (Rag1/2), whose expression is strictly controlled in a cell type-specific manner; the expression of Rag1/2 genes represents a hallmark of lymphoid lineage commitment. Although Rag genes are known to be evolutionally conserved among jawed vertebrates, how Rag genes are regulated by lineage-specific transcription factors (TFs) and how their regulatory system evolved among vertebrates have not been fully elucidated. Here, we reviewed the current body of knowledge concerning the cis-regulatory elements (CREs) of Rag genes and the evolution of the basic helix-loop-helix TF E protein regulating Rag gene CREs, as well as the evolution of the antagonist of this protein, the Id protein. This may help to understand how the adaptive immune system develops along with the evolution of responsible TFs and enhancers.

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.

The Multifaceted Role and Utility of MicroRNAs in Indolent B-Cell Non-Hodgkin Lymphomas

Normal B-cell development is a tightly regulated complex procedure, the deregulation of which can lead to lymphomagenesis. One common group of blood cancers is the B-cell non-Hodgkin lymphomas (NHLs), which can be categorized according to the proliferation and spread rate of cancer cells into indolent and aggressive ones. The most frequent indolent B-cell NHLs are follicular lymphoma and marginal zone lymphoma. MicroRNAs (miRNAs) are small non-coding RNAs that can greatly influence protein expression. Based on the multiple interactions among miRNAs and their targets, complex networks of gene expression regulation emerge, which normally are essential for proper B-cell development. Multiple miRNAs have been associated with B-cell lymphomas, as the deregulation of these complex networks can lead to such pathological states. The aim of the present review is to summarize the existing information regarding the multifaceted role of miRNAs in indolent B-cell NHLs, affecting the main B-cell subpopulations. We attempt to provide insight into their biological function, the complex miRNA-mRNA interactions, and their biomarker utility in these malignancies. Lastly, we address the limitations that hinder the investigation of the role of miRNAs in these lymphomas and discuss ways that these problems could be overcome in the future.