Temporal dynamics and genomic programming of plasma cell fates.

Affinity-matured plasma cells (PCs) of varying lifespans are generated through a germinal center (GC) response. The developmental dynamics and genomic programs of antigen-specific PC precursors remain to be elucidated. Here, using a model antigen in mice, we demonstrate biphasic generation of PC precursors, with those generating long-lived bone marrow PCs preferentially produced in the late phase of GC response. Clonal tracing using single-cell RNA sequencing and B cell antigen receptor sequencing in spleen and bone marrow compartments, coupled with adoptive transfer experiments, reveals a new PC transition state that gives rise to functionally competent PC precursors. The latter undergo clonal expansion, dependent on inducible expression of TIGIT. We propose a model for the proliferation and programming of precursors of long-lived PCs, based on extended antigen encounters in the GC.

An immunophenotype-coupled transcriptomic atlas of human hematopoietic progenitors.

Analysis of the human hematopoietic progenitor compartment is being transformed by single-cell multimodal approaches. Cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) enables coupled surface protein and transcriptome profiling, thereby revealing genomic programs underlying progenitor states. To perform CITE-seq systematically on primary human bone marrow cells, we used titrations with 266 CITE-seq antibodies (antibody-derived tags) and machine learning to optimize a panel of 132 antibodies. Multimodal analysis resolved >80 stem, progenitor, immune, stromal and transitional cells defined by distinctive surface markers and transcriptomes. This dataset enables flow cytometry solutions for in silico-predicted cell states and identifies dozens of cell surface markers consistently detected across donors spanning race and sex. Finally, aligning annotations from this atlas, we nominate normal marrow equivalents for acute myeloid leukemia stem cell populations that differ in clinical response. This atlas serves as an advanced digital resource for hematopoietic progenitor analyses in human health and disease.

ThPOK is a critical multifaceted regulator of myeloid lineage development.

The transcription factor ThPOK (encoded by Zbtb7b) is well known for its role as a master regulator of CD4 lineage commitment in the thymus. Here, we report an unexpected and critical role of ThPOK as a multifaceted regulator of myeloid lineage commitment, differentiation and maturation. Using reporter and knockout mouse models combined with single-cell RNA-sequencing, progenitor transfer and colony assays, we show that ThPOK controls monocyte-dendritic cell versus granulocyte lineage production during homeostatic differentiation, and serves as a brake for neutrophil maturation in granulocyte lineage-specified cells through transcriptional regulation of lineage-specific transcription factors and RNA via altered messenger RNA splicing to reprogram intron retention.

Coupled analysis of transcriptome and BCR mutations reveals role of OXPHOS in affinity maturation.

Antigen-activated B cells diversify variable regions of B cell antigen receptors by somatic hypermutation in germinal centers (GCs). The positive selection of GC B cells that acquire high-affinity mutations enables antibody affinity maturation. In spite of considerable progress, the genomic states underlying this process remain to be elucidated. Single-cell RNA sequencing and topic modeling revealed increased expression of the oxidative phosphorylation (OXPHOS) module in GC B cells undergoing mitoses. Coupled analysis of somatic hypermutation in immunoglobulin heavy chain (Igh) variable gene regions showed that GC B cells acquiring higher-affinity mutations had further elevated expression of OXPHOS genes. Deletion of mitochondrial Cox10 in GC B cells resulted in reduced cell division and impaired positive selection. Correspondingly, augmentation of OXPHOS activity with oltipraz promoted affinity maturation. We propose that elevated OXPHOS activity promotes B cell clonal expansion and also positive selection by tuning cell division times.

Ubiquitination of hnRNPA1 by TRAF6 links chronic innate immune signaling with myelodysplasia.

Toll-like receptor (TLR) activation contributes to premalignant hematologic conditions, such as myelodysplastic syndromes (MDS). TRAF6, a TLR effector with ubiquitin (Ub) ligase activity, is overexpressed in MDS hematopoietic stem/progenitor cells (HSPCs). We found that TRAF6 overexpression in mouse HSPC results in impaired hematopoiesis and bone marrow failure. Using a global Ub screen, we identified hnRNPA1, an RNA-binding protein and auxiliary splicing factor, as a substrate of TRAF6. TRAF6 ubiquitination of hnRNPA1 regulated alternative splicing of Arhgap1, which resulted in activation of the GTP-binding Rho family protein Cdc42 and accounted for hematopoietic defects in TRAF6-expressing HSPCs. These results implicate Ub signaling in coordinating RNA processing by TLR pathways during an immune response and in premalignant hematologic diseases, such as MDS.

Combinatorial Single-Cell Analyses of Granulocyte-Monocyte Progenitor Heterogeneity Reveals an Early Uni-potent Neutrophil Progenitor.

Granulocyte-monocyte progenitors (GMPs) have been previously defined for their potential to generate various myeloid progenies such as neutrophils and monocytes. Although studies have proposed lineage heterogeneity within GMPs, it is unclear if committed progenitors already exist among these progenitors and how they may behave differently during inflammation. By combining single-cell transcriptomic and proteomic analyses, we identified the early committed progenitor within the GMPs responsible for the strict production of neutrophils, which we designate as proNeu1. Our dissection of the GMP hierarchy led us to further identify a previously unknown intermediate proNeu2 population. Similar populations could be detected in human samples. proNeu1s, but not proNeu2s, selectively expanded during the early phase of sepsis at the expense of monocytes. Collectively, our findings help shape the neutrophil maturation trajectory roadmap and challenge the current definition of GMPs.

In situ mapping identifies distinct vascular niches for myelopoiesis.

In contrast to nearly all other tissues, the anatomy of cell differentiation in the bone marrow remains unknown. This is owing to a lack of strategies for examining myelopoiesis-the differentiation of myeloid progenitors into a large variety of innate immune cells-in situ in the bone marrow. Such strategies are required to understand differentiation and lineage-commitment decisions, and to define how spatial organizing cues inform tissue function. Here we develop approaches for imaging myelopoiesis in mice, and generate atlases showing the differentiation of granulocytes, monocytes and dendritic cells. The generation of granulocytes and dendritic cells-monocytes localizes to different blood-vessel structures known as sinusoids, and displays lineage-specific spatial and clonal architectures. Acute systemic infection with Listeria monocytogenes induces lineage-specific progenitor clusters to undergo increased self-renewal of progenitors, but the different lineages remain spatially separated. Monocyte-dendritic cell progenitors (MDPs) map with nonclassical monocytes and conventional dendritic cells; these localize to a subset of blood vessels expressing a major regulator of myelopoiesis, colony-stimulating factor 1 (CSF1, also known as M-CSF)1. Specific deletion of Csf1 in endothelium disrupts the architecture around MDPs and their localization to sinusoids. Subsequently, there are fewer MDPs and their ability to differentiate is reduced, leading to a loss of nonclassical monocytes and dendritic cells during both homeostasis and infection. These data indicate that local cues produced by distinct blood vessels are responsible for the spatial organization of definitive blood cell differentiation.

Ontogeny and function of the circadian clock in intestinal organoids.

Circadian rhythms regulate diverse aspects of gastrointestinal physiology ranging from the composition of microbiota to motility. However, development of the intestinal circadian clock and detailed mechanisms regulating circadian physiology of the intestine remain largely unknown. In this report, we show that both pluripotent stem cell-derived human intestinal organoids engrafted into mice and patient-derived human intestinal enteroids possess circadian rhythms and demonstrate circadian phase-dependent necrotic cell death responses to Clostridium difficile toxin B (TcdB). Intriguingly, mouse and human enteroids demonstrate anti-phasic necrotic cell death responses to TcdB. RNA-Seq analysis shows that ~3-10% of the detectable transcripts are rhythmically expressed in mouse and human enteroids. Remarkably, we observe anti-phasic gene expression of Rac1, a small GTPase directly inactivated by TcdB, between mouse and human enteroids, and disruption of Rac1 abolishes clock-dependent necrotic cell death responses. Our findings uncover robust functions of circadian rhythms regulating clock-controlled genes in both mouse and human enteroids governing organism-specific, circadian phase-dependent necrotic cell death responses, and lay a foundation for human organ- and disease-specific investigation of clock functions using human organoids for translational applications.

Granulocyte-Monocyte Progenitors and Monocyte-Dendritic Cell Progenitors Independently Produce Functionally Distinct Monocytes.

Granulocyte-monocyte progenitors (GMPs) and monocyte-dendritic cell progenitors (MDPs) produce monocytes during homeostasis and in response to increased demand during infection. Both progenitor populations are thought to derive from common myeloid progenitors (CMPs), and a hierarchical relationship (CMP-GMP-MDP-monocyte) is presumed to underlie monocyte differentiation. Here, however, we demonstrate that mouse MDPs arose from CMPs independently of GMPs, and that GMPs and MDPs produced monocytes via similar but distinct monocyte-committed progenitors. GMPs and MDPs yielded classical (Ly6Chi) monocytes with gene expression signatures that were defined by their origins and impacted their function. GMPs produced a subset of "neutrophil-like" monocytes, whereas MDPs gave rise to a subset of monocytes that yielded monocyte-derived dendritic cells. GMPs and MDPs were also independently mobilized to produce specific combinations of myeloid cell types following the injection of microbial components. Thus, the balance of GMP and MDP differentiation shapes the myeloid cell repertoire during homeostasis and following infection.

Mouse models of neutropenia reveal progenitor-stage-specific defects.

Advances in genetics and sequencing have identified a plethora of disease-associated and disease-causing genetic alterations. To determine causality between genetics and disease, accurate models for molecular dissection are required; however, the rapid expansion of transcriptional populations identified through single-cell analyses presents a major challenge for accurate comparisons between mutant and wild-type cells. Here we generate mouse models of human severe congenital neutropenia (SCN) using patient-derived mutations in the GFI1 transcription factor. To determine the effects of SCN mutations, we generated single-cell references for granulopoietic genomic states with linked epitopes1, aligned mutant cells to their wild-type equivalents and identified differentially expressed genes and epigenetic loci. We find that GFI1-target genes are altered sequentially, as cells go through successive states of differentiation. These insights facilitated the genetic rescue of granulocytic specification but not post-commitment defects in innate immune effector function, and underscore the importance of evaluating the effects of mutations and therapy within each relevant cell state.

Retinoid X receptor promotes hematopoietic stem cell fitness and quiescence and preserves hematopoietic homeostasis.

Hematopoietic stem cells (HSCs) balance self-renewal and differentiation to maintain hematopoietic fitness throughout life. In steady-state conditions, HSC exhaustion is prevented by the maintenance of most HSCs in a quiescent state, with cells entering the cell cycle only occasionally. HSC quiescence is regulated by retinoid and fatty-acid ligands of transcriptional factors of the nuclear retinoid X receptor (RXR) family. Herein, we show that dual deficiency for hematopoietic RXRα and RXRß induces HSC exhaustion, myeloid cell/megakaryocyte differentiation, and myeloproliferative-like disease. RXRα and RXRß maintain HSC quiescence, survival, and chromatin compaction; moreover, transcriptome changes in RXRα;RXRß-deficient HSCs include premature acquisition of an aging-like HSC signature, MYC pathway upregulation, and RNA intron retention. Fitness loss and associated RNA transcriptome and splicing alterations in RXRα;RXRß-deficient HSCs are prevented by Myc haploinsufficiency. Our study reveals the critical importance of RXRs for the maintenance of HSC fitness and their protection from premature aging.

Author Correction: Single-cell analysis of mixed-lineage states leading to a binary cell fate choice.

In this Letter, the first name of author Virendra K. Chaudhri was incorrectly spelled 'Viren'; author Meenakshi Venkatasubramanian should also be associated with 'Department of Electrical Engineering and Computer Science, University of Cincinnati, Cincinnati, Ohio 45221, USA'; authors Bruce J. Aronow, Nathan Salomonis, Harinder Singh and H. Leighton Grimes should also be associated with 'Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati, Ohio 45229, USA'. The Letter has not been corrected online.

Bromodomain inhibition overcomes treatment resistance in distinct molecular subtypes of melanoma.

Therapy of BRAF-mutant melanoma with selective inhibitors of BRAF (BRAFi) and MEK (MEKi) represents a major clinical advance but acquired resistance to therapy has emerged as a key obstacle. To date, no clinical approaches successfully resensitize to BRAF/MEK inhibition. Here, we develop a therapeutic strategy for melanoma using bromosporine, a bromodomain inhibitor. Bromosporine (bromo) monotherapy produced significant anti-tumor effects against established melanoma cell lines and patient-derived xenografts (PDXs). Combinatorial therapy involving bromosporine and cobimetinib (bromo/cobi) showed synergistic anti-tumor effects in multiple BRAFi-resistant PDX models. The bromo/cobi combination was superior in vivo to standard BRAFi/MEKi therapy in the treatment-naive BRAF-mutant setting and to MEKi alone in the setting of immunotherapy-resistant NRAS- and NF1-mutant melanoma. RNA sequencing of xenografts treated with bromo/cobi revealed profound down-regulation of genes critical to cell division and mitotic progression. Bromo/cobi treatment resulted in marked DNA damage and cell-cycle arrest, resulting in induction of apoptosis. These studies introduce bromodomain inhibition, alone or combined with agents targeting the mitogen activated protein kinase pathway, as a rational therapeutic approach for melanoma refractory to standard targeted or immunotherapeutic approaches.

CellDrift: inferring perturbation responses in temporally sampled single-cell data.

Cells and tissues respond to perturbations in multiple ways that can be sensitively reflected in the alterations of gene expression. Current approaches to finding and quantifying the effects of perturbations on cell-level responses over time disregard the temporal consistency of identifiable gene programs. To leverage the occurrence of these patterns for perturbation analyses, we developed CellDrift (https://github.com/KANG-BIOINFO/CellDrift), a generalized linear model-based functional data analysis method that is capable of identifying covarying temporal patterns of various cell types in response to perturbations. As compared to several other approaches, CellDrift demonstrated superior performance in the identification of temporally varied perturbation patterns and the ability to impute missing time points. We applied CellDrift to multiple longitudinal datasets, including COVID-19 disease progression and gastrointestinal tract development, and demonstrated its ability to identify specific gene programs associated with sequential biological processes, trajectories and outcomes.

pyInfinityFlow: optimized imputation and analysis of high-dimensional flow cytometry data for millions of cells.

MOTIVATION: While conventional flow cytometry is limited to dozens of markers, new experimental and computational strategies, such as Infinity Flow, allow for the generation and imputation of hundreds of cell surface protein markers in millions of cells. Here, we describe an end-to-end analysis workflow for Infinity Flow data in Python. RESULTS: pyInfinityFlow enables the efficient analysis of millions of cells, without down-sampling, through direct integration with well-established Python packages for single-cell genomics analysis. pyInfinityFlow accurately identifies both common and extremely rare cell populations which are challenging to define from single-cell genomics studies alone. We demonstrate that this workflow can nominate novel markers to design new flow cytometry gating strategies for predicted cell populations. pyInfinityFlow can be extended to diverse cell discovery analyses with flexibility to adapt to diverse Infinity Flow experimental designs. AVAILABILITY AND IMPLEMENTATION: pyInfinityFlow is freely available in GitHub (https://github.com/KyleFerchen/pyInfinityFlow) and on PyPI (https://pypi.org/project/pyInfinityFlow/). Package documentation with tutorials on a test dataset is available by Read the Docs (pyinfinityflow.readthedocs.io). The scripts and data for reproducing the results are available at https://github.com/KyleFerchen/pyInfinityFlow/tree/main/analysis_scripts, along with the raw flow cytometry input data.

miR-876-3p is a tumor suppressor on 9p21 that is inactivated in melanoma and targets ERK.

BACKGROUND: While melanomas commonly harbor losses of 9p21, on which CDKN2A resides, the presence of additional tumor suppressor elements at this locus is incompletely characterized. Here we assess the expression levels and functional role of microRNA-876-3p (miR-876), whose gene also maps to 9p21. METHODS: Expression of miR-876 was assessed in human tissues and cell lines using quantitative miRNA reverse transcriptase polymerase chain reaction (qRT-PCR). MIR876 copy number was determined in The Cancer Genome Atlas (TCGA) melanoma cohort. The consequences of regulation of miR-876 expression were assessed on melanoma cell colony formation, migration, invasion, apoptosis, cell cycle progression, and drug sensitivity in culture, and on in vivo tumor growth in a xenograft model. Genome-wide transcriptomic changes induced by miR-876 overexpression were determined using RNA sequencing (RNA-Seq). RESULTS: miR-876 expression was significantly decreased in primary melanoma samples when compared with nevi, and in human melanoma cell lines when compared with human melanocytes. Analysis of the TCGA cohort revealed deletions in MIR876 in > 50% of melanomas. miR-876 overexpression resulted in decreased melanoma cell colony formation, migration, and invasion, which was accompanied by cell cycle arrest and increased apoptosis. Intra-tumoral injections of miR-876 significantly suppressed melanoma growth in vivo. RNA-Seq analysis of miR-876-treated tumors revealed downregulation of several growth-promoting genes, along with upregulation of tumor suppressor genes, which was confirmed by qRT-PCR analysis. Computational analyses identified MAPK1 (or ERK2) as a possible target of miR-876 action. Overexpression of miR-876 significantly suppressed luciferase expression driven by the MAPK1/ERK2 3' UTR, and resulted in decreased ERK protein expression in melanoma cells. MAPK1/ERK2 cDNA overexpression rescued the effects of miR-876 on melanoma colony formation. miR-876 overexpression sensitized melanoma cells to treatment with the BRAF inhibitor vemurafenib. CONCLUSIONS: These studies identify miR-876 as a distinct tumor suppressor on 9p21 that is inactivated in melanoma and suggest miR-876 loss as an additional mechanism to activate ERK and the mitogen activated protein kinase (MAPK) pathway in melanoma. In addition, they suggest the therapeutic potential of combining miR-876 overexpression with BRAF inhibition as a rational therapeutic strategy for melanoma.

Erythroblastic islands foster granulopoiesis in parallel to terminal erythropoiesis.

The erythroblastic island (EBI), composed of a central macrophage surrounded by maturing erythroblasts, is the erythroid precursor niche. Despite numerous studies, its precise composition is still unclear. Using multispectral imaging flow cytometry, in vitro island reconstitution, and single-cell RNA sequencing of adult mouse bone marrow (BM) EBI-component cells enriched by gradient sedimentation, we present evidence that the CD11b+ cells present in the EBIs are neutrophil precursors specifically associated with BM EBI macrophages, indicating that erythro-(myelo)-blastic islands are a site for terminal granulopoiesis and erythropoiesis. We further demonstrate that the balance between these dominant and terminal differentiation programs is dynamically regulated within this BM niche by pathophysiological states that favor granulopoiesis during anemia of inflammation and favor erythropoiesis after erythropoietin stimulation. Finally, by molecular profiling, we reveal the heterogeneity of EBI macrophages by cellular indexing of transcriptome and epitope sequencing of mouse BM EBIs at baseline and after erythropoietin stimulation in vivo and provide a searchable online viewer of these data characterizing the macrophage subsets serving as hematopoietic niches. Taken together, our findings demonstrate that EBIs serve a dual role as niches for terminal erythropoiesis and granulopoiesis and the central macrophages adapt to optimize production of red blood cells or neutrophils.

Pbx4 limits heart size and fosters arch artery formation by partitioning second heart field progenitors and restricting proliferation.

Vertebrate heart development requires the integration of temporally distinct differentiating progenitors. However, few signals are understood that restrict the size of the later-differentiating outflow tract (OFT). We show that improper specification and proliferation of second heart field (SHF) progenitors in zebrafish lazarus (lzr) mutants, which lack the transcription factor Pbx4, produces enlarged hearts owing to an increase in ventricular and smooth muscle cells. Specifically, Pbx4 initially promotes the partitioning of the SHF into anterior progenitors, which contribute to the OFT, and adjacent endothelial cell progenitors, which contribute to posterior pharyngeal arches. Subsequently, Pbx4 limits SHF progenitor (SHFP) proliferation. Single cell RNA sequencing of nkx2.5+ cells revealed previously unappreciated distinct differentiation states and progenitor subpopulations that normally reside within the SHF and arterial pole of the heart. Specifically, the transcriptional profiles of Pbx4-deficient nkx2.5+ SHFPs are less distinct and display characteristics of normally discrete proliferative progenitor and anterior, differentiated cardiomyocyte populations. Therefore, our data indicate that the generation of proper OFT size and arch arteries requires Pbx-dependent stratification of unique differentiation states to facilitate both homeotic-like transformations and limit progenitor production within the SHF.

DeepImmuno: deep learning-empowered prediction and generation of immunogenic peptides for T-cell immunity.

Cytolytic T-cells play an essential role in the adaptive immune system by seeking out, binding and killing cells that present foreign antigens on their surface. An improved understanding of T-cell immunity will greatly aid in the development of new cancer immunotherapies and vaccines for life-threatening pathogens. Central to the design of such targeted therapies are computational methods to predict non-native peptides to elicit a T-cell response, however, we currently lack accurate immunogenicity inference methods. Another challenge is the ability to accurately simulate immunogenic peptides for specific human leukocyte antigen alleles, for both synthetic biological applications, and to augment real training datasets. Here, we propose a beta-binomial distribution approach to derive peptide immunogenic potential from sequence alone. We conducted systematic benchmarking of five traditional machine learning (ElasticNet, K-nearest neighbors, support vector machine, Random Forest and AdaBoost) and three deep learning models (convolutional neural network (CNN), Residual Net and graph neural network) using three independent prior validated immunogenic peptide collections (dengue virus, cancer neoantigen and SARS-CoV-2). We chose the CNN as the best prediction model, based on its adaptivity for small and large datasets and performance relative to existing methods. In addition to outperforming two highly used immunogenicity prediction algorithms, DeepImmuno-CNN correctly predicts which residues are most important for T-cell antigen recognition and predicts novel impacts of SARS-CoV-2 variants. Our independent generative adversarial network (GAN) approach, DeepImmuno-GAN, was further able to accurately simulate immunogenic peptides with physicochemical properties and immunogenicity predictions similar to that of real antigens. We provide DeepImmuno-CNN as source code and an easy-to-use web interface.

PHIP drives glioblastoma motility and invasion by regulating the focal adhesion complex.

The invasive behavior of glioblastoma is essential to its aggressive potential. Here, we show that pleckstrin homology domain interacting protein (PHIP), acting through effects on the force transduction layer of the focal adhesion complex, drives glioblastoma motility and invasion. Immunofluorescence analysis localized PHIP to the leading edge of glioblastoma cells, together with several focal adhesion proteins: vinculin (VCL), talin 1 (TLN1), integrin beta 1 (ITGB1), as well as phosphorylated forms of paxillin (pPXN) and focal adhesion kinase (pFAK). Confocal microscopy specifically localized PHIP to the force transduction layer, together with TLN1 and VCL. Immunoprecipitation revealed a physical interaction between PHIP and VCL. Targeted suppression of PHIP resulted in significant down-regulation of these focal adhesion proteins, along with zyxin (ZYX), and produced profoundly disorganized stress fibers. Live-cell imaging of glioblastoma cells overexpressing a ZYX-GFP construct demonstrated a role for PHIP in regulating focal adhesion dynamics. PHIP silencing significantly suppressed the migratory and invasive capacity of glioblastoma cells, partially restored following TLN1 or ZYX cDNA overexpression. PHIP knockdown produced substantial suppression of tumor growth upon intracranial implantation, as well as significantly reduced microvessel density and secreted VEGF levels. PHIP copy number was elevated in the classical glioblastoma subtype and correlated with elevated EGFR levels. These results demonstrate PHIP's role in regulating the actin cytoskeleton, focal adhesion dynamics, and tumor cell motility, and identify PHIP as a key driver of glioblastoma migration and invasion.