Integrative single-cell chromatin and transcriptome analysis of human plasma cell differentiation.

ABSTRACT: Plasma cells (PCs) are highly specialized cells representing the end stage of B-cell differentiation. We have shown that PC differentiation can be reproduced in vitro using elaborate culture systems. The molecular changes occurring during PC differentiation are recapitulated in this in vitro differentiation model. However, a major challenge exists to decipher the spatiotemporal epigenetic and transcriptional programs that drive the early stages of PC differentiation. We combined single cell (sc) RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin with high throughput sequencing (scATAC-seq) to decipher the trajectories involved in PC differentiation. ScRNA-seq experiments revealed a strong heterogeneity of the preplasmablastic and plasmablastic stages. Among genes that were commonly identified using scATAC-seq and scRNA-seq, we identified several transcription factors with significant stage specific potential importance in PC differentiation. Interestingly, differentially accessible peaks characterizing the preplasmablastic stage were enriched in motifs of BATF3, FOS and BATF, belonging to activating protein 1 (AP-1) transcription factor family that may represent key transcriptional nodes involved in PC differentiation. Integration of transcriptomic and epigenetic data at the single cell level revealed that a population of preplasmablasts had already undergone epigenetic remodeling related to PC profile together with unfolded protein response activation and are committed to differentiate in PC. These results and the supporting data generated with our in vitro PC differentiation model provide a unique resource for the identification of molecular circuits that are crucial for early and mature PC maturation and biological functions. These data thus provide critical insights into epigenetic- and transcription-mediated reprogramming events that sustain PC differentiation.

ShIVA: a user-friendly and interactive interface giving biologists control over their single-cell RNA-seq data.

Single-cell technologies have revolutionised biological research and applications. As they continue to evolve with multi-omics and spatial resolution, analysing single-cell datasets is becoming increasingly complex. For biologists lacking expert data analysis resources, the problem is even more crucial, even for the simplest single-cell transcriptomics datasets. We propose ShIVA, an interface for the analysis of single-cell RNA-seq and CITE-seq data specifically dedicated to biologists. Intuitive, iterative and documented by video tutorials, ShIVA allows biologists to follow a robust and reproducible analysis process, mostly based on the Seurat v4 R package, to fully explore and quantify their dataset, to produce useful figures and tables and to export their work to allow more complex analyses performed by experts.

Germinal center-dependent and -independent immune responses of tumor-infiltrating B cells in human cancers.

B cells play essential roles in immunity, mainly through the production of high affinity plasma cells (PCs) and memory B (Bmem) cells. The affinity maturation and differentiation of B cells rely on the integration of B-cell receptor (BCR) intrinsic and extrinsic signals provided by antigen binding and the microenvironment, respectively. In recent years, tumor infiltrating B (TIL-B) cells and PCs (TIL-PCs) have been revealed as important players in antitumor responses in human cancers, but their interplay and dynamics remain largely unknown. In lymphoid organs, B-cell responses involve both germinal center (GC)-dependent and GC-independent pathways for Bmem cell and PC production. Affinity maturation of BCR repertoires occurs in GC reactions with specific spatiotemporal dynamics of signal integration by B cells. In general, the reactivation of high-affinity Bmem cells by antigens triggers GC-independent production of large numbers of PC without BCR rediversification. Understanding B-cell dynamics in immune responses requires the integration of multiple tools and readouts such as single-cell phenotyping and RNA-seq, in situ analyses, BCR repertoire analysis, BCR specificity and affinity assays, and functional tests. Here, we review how those tools have recently been applied to study TIL-B cells and TIL-PC in different types of solid tumors. We assessed the published evidence for different models of TIL-B-cell dynamics involving GC-dependent or GC-independent local responses and the resulting production of antigen-specific PCs. Altogether, we highlight the need for more integrative B-cell immunology studies to rationally investigate TIL-B cells as a leverage for antitumor therapies.

Novel mouse models based on intersectional genetics to identify and characterize plasmacytoid dendritic cells.

Plasmacytoid dendritic cells (pDCs) are the main source of type I interferon (IFN-I) during viral infections. Their other functions are debated, due to a lack of tools to identify and target them in vivo without affecting pDC-like cells and transitional DCs (tDCs), which harbor overlapping phenotypes and transcriptomes but a higher efficacy for T cell activation. In the present report, we present a reporter mouse, pDC-Tom, designed through intersectional genetics based on unique Siglech and Pacsin1 coexpression in pDCs. The pDC-Tom mice specifically tagged pDCs and, on breeding with Zbtb46GFP mice, enabled transcriptomic profiling of all splenic DC types, unraveling diverging activation of pDC-like cells versus tDCs during a viral infection. The pDC-Tom mice also revealed initially similar but later divergent microanatomical relocation of splenic IFN+ versus IFN- pDCs during infection. The mouse models and specific gene modules we report here will be useful to delineate the physiological functions of pDCs versus other DC types.

Neuroblasts contribute to oligodendrocytes generation upon demyelination in the adult mouse brain.

After demyelinating insult, the neuronal progenitors of the adult mouse sub-ventricular zone (SVZ) called neuroblasts convert into oligodendrocytes that participate to the remyelination process. We use this rare example of spontaneous fate conversion to identify the molecular mechanisms governing these processes. Using in vivo cell lineage and single cell RNA-sequencing, we demonstrate that SVZ neuroblasts fate conversion proceeds through formation of a non-proliferating transient cellular state co-expressing markers of both neuronal and oligodendrocyte identities. Transition between the two identities starts immediately after demyelination and occurs gradually, by a stepwise upregulation/downregulation of key TFs and chromatin modifiers. Each step of this fate conversion involves fine adjustments of the transcription and translation machineries as well as tight regulation of metabolism and migratory behaviors. Together, these data constitute the first in-depth analysis of a spontaneous cell fate conversion in the adult mammalian CNS.

Neuropilin-1 cooperates with PD-1 in CD8+ T cells predicting outcomes in melanoma patients treated with anti-PD1.

Targeting immune checkpoints, such as Programmed cell Death 1 (PD1), has improved survival in cancer patients by restoring antitumor immune responses. Most patients, however, relapse or are refractory to immune checkpoint blocking therapies. Neuropilin-1 (NRP1) is a transmembrane glycoprotein required for nervous system and angiogenesis embryonic development, also expressed in immune cells. We hypothesized that NRP1 could be an immune checkpoint co-receptor modulating CD8+ T cells activity in the context of the antitumor immune response. Here, we show that NRP1 is recruited in the cytolytic synapse of PD1+CD8+ T cells, cooperates and enhances PD-1 activity. In mice, CD8+ T cells specific deletion of Nrp1 improves anti-PD1 antibody antitumor immune responses. Likewise, in human metastatic melanoma, the expression of NRP1 in tumor infiltrating CD8+ T cells predicts poor outcome of patients treated with anti-PD1. NRP1 is a promising target to overcome resistance to anti-PD1 therapies.

Viral infection engenders bona fide and bystander subsets of lung-resident memory B cells through a permissive mechanism.

Lung-resident memory B cells (MBCs) provide localized protection against reinfection in respiratory airways. Currently, the biology of these cells remains largely unexplored. Here, we combined influenza and SARS-CoV-2 infection with fluorescent-reporter mice to identify MBCs regardless of antigen specificity. We found that two main transcriptionally distinct subsets of MBCs colonized the lung peribronchial niche after infection. These subsets arose from different progenitors and were both class switched, somatically mutated, and intrinsically biased in their differentiation fate toward plasma cells. Combined analysis of antigen specificity and B cell receptor repertoire segregated these subsets into "bona fide" virus-specific MBCs and "bystander" MBCs with no apparent specificity for eliciting viruses generated through an alternative permissive process. Thus, diverse transcriptional programs in MBCs are not linked to specific effector fates but rather to divergent strategies of the immune system to simultaneously provide rapid protection from reinfection while diversifying the initial B cell repertoire.

NF-κB-dependent IRF1 activation programs cDC1 dendritic cells to drive antitumor immunity.

Conventional type 1 dendritic cells (cDC1s) are critical for antitumor immunity. They acquire antigens from dying tumor cells and cross-present them to CD8+ T cells, promoting the expansion of tumor-specific cytotoxic T cells. However, the signaling pathways that govern the antitumor functions of cDC1s in immunogenic tumors are poorly understood. Using single-cell transcriptomics to examine the molecular pathways regulating intratumoral cDC1 maturation, we found nuclear factor κB (NF-κB) and interferon (IFN) pathways to be highly enriched in a subset of functionally mature cDC1s. We identified an NF-κB-dependent and IFN-γ-regulated gene network in cDC1s, including cytokines and chemokines specialized in the recruitment and activation of cytotoxic T cells. By mapping the trajectory of intratumoral cDC1 maturation, we demonstrated the dynamic reprogramming of tumor-infiltrating cDC1s by NF-κB and IFN signaling pathways. This maturation process was perturbed by specific inactivation of either NF-κB or IFN regulatory factor 1 (IRF1) in cDC1s, resulting in impaired expression of IFN-γ-responsive genes and consequently a failure to efficiently recruit and activate antitumoral CD8+ T cells. Last, we demonstrate the relevance of these findings to patients with melanoma, showing that activation of the NF-κB/IRF1 axis in association with cDC1s is linked with improved clinical outcome. The NF-κB/IRF1 axis in cDC1s may therefore represent an important focal point for the development of new diagnostic and therapeutic approaches to improve cancer immunotherapy.

Heterogeneity of germinal center B cells: New insights from single-cell studies.

Upon antigen exposure, activated B cells in antigen-draining lymphoid organs form microanatomical structures, called germinal centers (GCs), where affinity maturation occurs. Within the GC microenvironment, GC B cells undergo proliferation and B cell receptor (BCR) genes somatic hypermutation in the dark zone (DZ), and affinity-based selection in the light zone (LZ). In the current paradigm of GC dynamics, high-affinity LZ B cells may be selected by cognate T- follicular helper cells to either differentiate into plasma cells or memory B cells, or re-enter the DZ and initiate a new round of proliferation and BCR diversification, before migrating back to the LZ. Given the diversity of cell states and potential cell fates that GC B cells may adopt, the two-state DZ-LZ paradigm has been challenged by studies that explored GC B-cell heterogeneity with a variety of single-cell technologies. Here, we review studies and single-cell technologies which have allowed to refine the working model of GC B-cell cellular and molecular heterogeneity during affinity maturation. This review also covers the use of single-cell quantitative data for mathematical modeling of GC reactions, and the application of single-cell genomics to the study of GC-derived malignancies.

Follicular lymphoma dynamics.

Follicular lymphoma (FL) is an indolent yet challenging disease. Despite a generally favorable response to immunochemotherapy regimens, a fraction of patients does not respond or relapses early with unfavorable prognosis. For the vast majority of those who initially respond, relapses will repeatedly occur with increasing refractoriness to available treatments. Addressing the clinical challenges in FL warrants deep understanding of the nature of treatment-resistant FL cells seeding relapses, and of the biological basis of early disease progression. Great progress has been made in the last decade in the description and interrogation of the (epi)genomic landscape of FL cells, of their major dependency to the tumor microenvironment (TME), and of the stepwise lymphomagenesis process, from healthy to subclinical disease and to overt FL. A new picture is emerging, in which an ever-evolving tumor-TME duo sparks a complex and multilayered clonal and functional heterogeneity, blurring the discovery of prognostic biomarkers, patient stratification and reliable designs of risk-adapted treatments. Novel technological approaches allowing to decipher both tumor and TME heterogeneity at the single-cell level are beginning to unravel unsuspected cell dynamics and plasticity of FL cells. The upcoming drawing of a comprehensive functional picture of FL within its ecosystem holds great promise to address the unmet medical needs of this complex lymphoma.

Single-cell profiling reveals the trajectories of natural killer cell differentiation in bone marrow and a stress signature induced by acute myeloid leukemia.

Natural killer (NK) cells are innate cytotoxic lymphoid cells (ILCs) involved in the killing of infected and tumor cells. Among human and mouse NK cells from the spleen and blood, we previously identified by single-cell RNA sequencing (scRNAseq) two similar major subsets, NK1 and NK2. Using the same technology, we report here the identification, by single-cell RNA sequencing (scRNAseq), of three NK cell subpopulations in human bone marrow. Pseudotime analysis identified a subset of resident CD56bright NK cells, NK0 cells, as the precursor of both circulating CD56dim NK1-like NK cells and CD56bright NK2-like NK cells in human bone marrow and spleen under physiological conditions. Transcriptomic profiles of bone marrow NK cells from patients with acute myeloid leukemia (AML) exhibited stress-induced repression of NK cell effector functions, highlighting the profound impact of this disease on NK cell heterogeneity. Bone marrow NK cells from AML patients exhibited reduced levels of CD160, but the CD160high group had a significantly higher survival rate.

TLR-9 agonist and CD40-targeting vaccination induces HIV-1 envelope-specific B cells with a diversified immunoglobulin repertoire in humanized mice.

The development of HIV-1 vaccines is challenged by the lack of relevant models to accurately induce human B- and T-cell responses in lymphoid organs. In humanized mice reconstituted with human hematopoietic stem cells (hu-mice), human B cell-development and function are impaired and cells fail to efficiently transition from IgM B cells to IgG B cells. Here, we found that CD40-targeted vaccination combined with CpG-B adjuvant overcomes the usual defect of human B-cell switch and maturation in hu-mice. We further dissected hu-B cell responses directed against the HIV-1 Env protein elicited by targeting Env gp140 clade C to the CD40 receptor of antigen-presenting cells. The anti-CD40.Env gp140 vaccine was injected with CpG-B in a homologous prime/boost regimen or as a boost of a NYVAC-KC pox vector encoding Env gp140 clade C. Both regimens elicited Env-specific IgG-switched memory hu-B cells at a greater magnitude in hu-mice primed with NYVAC-KC. Single-cell RNA-seq analysis showed gp140-specific hu-B cells to express polyclonal IgG1 and IgG3 isotypes and a broad Ig VH/VL repertoire, with predominant VH3 family gene usage. These cells exhibited a higher rate of somatic hypermutation than the non-specific IgG+ hu-B-cell counterpart. Both vaccine regimens induced splenic GC-like structures containing hu-B and hu-Tfh-like cells expressing PD-1 and BCL-6. We confirmed in this model that circulating ICOS+ memory hu-Tfh cells correlated with the magnitude of gp140-specific B-cell responses. Finally, the NYVAC-KC heterologous prime led to a more diverse clonal expansion of specific hu-B cells. Thus, this study shows that CD40-targeted vaccination induces human IgG production in hu-mice and provides insights for the development of a CD40-targeting vaccine to prevent HIV-1 infection in humans.

Corrigendum: FB5P-seq: FACS-Based 5-Prime End Single-Cell RNA-seq for Integrative Analysis of Transcriptome and Antigen Receptor Repertoire in B and T Cells.

[This corrects the article DOI: 10.3389/fimmu.2020.00216.].

Corrigendum: FB5P-seq: FACS-Based 5-Prime End Single-Cell RNA-seq for Integrative Analysis of Transcriptome and Antigen Receptor Repertoire in B and T Cells.

[This corrects the article DOI: 10.3389/fimmu.2020.00216.].

Integrative molecular profiling of autoreactive CD4 T cells in autoimmune hepatitis.

BACKGROUND & AIMS: In most autoimmune disorders, crosstalk of B cells and CD4 T cells results in the accumulation of autoantibodies. In autoimmune hepatitis (AIH), the presence of anti-soluble liver antigen (SLA) autoantibodies is associated with reduced overall survival, but the associated autoreactive CD4 T cells have not yet been characterised. Herein, we isolated and deeply characterised SLA-specific CD4 T cells in patients with AIH. METHODS: We used brief ex vivo restimulation with overlapping SLA peptides to isolate and phenotype circulating SLA-specific CD4 T cells, and integrative single-cell RNA-seq (scRNA-seq) to characterise their transcriptome and T-cell receptor (TCR) repertoire. Autoreactive TCRs were cloned and used to identify dominant SLA-derived epitopes. SLA-specific CD4 T cells were tracked in peripheral blood through TCR sequencing to identify their phenotypic niche. We further characterised disease-associated peripheral blood T cells by high-content flow cytometry in 42 patients with AIH and 17 controls with non-alcoholic steatohepatitis. RESULTS: Autoreactive SLA-specific CD4 T cells were only detected in patients with anti-SLA autoantibodies and had a memory PD-1+CXCR5-CCR6-CD27+ phenotype. ScRNA-seq revealed their pro-inflammatory/B-helper profile. SLA81-100 and SLA177-204 contain dominant T-cell epitopes. Autoreactive TCR clonotypes were predominantly found in the memory PD-1+CXCR5-CD4 T cells, which were significantly increased in the blood of patients with AIH and supported B-cell differentiation through IL-21. Finally, we identified specific T-cell phenotypes linked to disease activity and IgG level during AIH. CONCLUSIONS: We provide a deep characterisation of rare circulating autoreactive CD4 T cells and identify their peripheral reservoir in AIH. We also propose a specific phenotype of autoreactive T cells related to AIH disease activity, which will be essential to track, delineate, and potentially target these pathogenic cells. LAY SUMMARY: One principal characteristic of autoimmune hepatitis (AIH), like for many other autoimmune diseases, is the accumulation of autoantibodies produced by B lymphocytes following their interaction with autoreactive CD4 T lymphocytes. In this study, we identified and characterised with high resolution these CD4 T cells. This will be essential to track, delineate, and potentially target them during AIH.

The activation trajectory of plasmacytoid dendritic cells in vivo during a viral infection.

Plasmacytoid dendritic cells (pDCs) are a major source of type I interferon (IFN-I). What other functions pDCs exert in vivo during viral infections is controversial, and more studies are needed to understand their orchestration. In the present study, we characterize in depth and link pDC activation states in animals infected by mouse cytomegalovirus by combining Ifnb1 reporter mice with flow cytometry, single-cell RNA sequencing, confocal microscopy and a cognate CD4 T cell activation assay. We show that IFN-I production and T cell activation were performed by the same pDC, but these occurred sequentially in time and in different micro-anatomical locations. In addition, we show that pDC commitment to IFN-I production was marked early on by their downregulation of leukemia inhibitory factor receptor and was promoted by cell-intrinsic tumor necrosis factor signaling. We propose a new model for how individual pDCs are endowed to exert different functions in vivo during a viral infection, in a manner tightly orchestrated in time and space.

FB5P-seq: FACS-Based 5-Prime End Single-Cell RNA-seq for Integrative Analysis of Transcriptome and Antigen Receptor Repertoire in B and T Cells.

Single-cell RNA sequencing (scRNA-seq) allows the identification, characterization, and quantification of cell types in a tissue. When focused on B and T cells of the adaptive immune system, scRNA-seq carries the potential to track the clonal lineage of each analyzed cell through the unique rearranged sequence of its antigen receptor (BCR or TCR, respectively) and link it to the functional state inferred from transcriptome analysis. Here we introduce FB5P-seq, a FACS-based 5'-end scRNA-seq method for cost-effective, integrative analysis of transcriptome and paired BCR or TCR repertoire in phenotypically defined B and T cell subsets. We describe in detail the experimental workflow and provide a robust bioinformatics pipeline for computing gene count matrices and reconstructing repertoire sequences from FB5P-seq data. We further present two applications of FB5P-seq for the analysis of human tonsil B cell subsets and peripheral blood antigen-specific CD4 T cells. We believe that our novel integrative scRNA-seq method will be a valuable option to study rare adaptive immune cell subsets in immunology research.