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

Plasma cells (PC) 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 drives the early stages of PC differentiation. We combined single cell (sc) RNA-seq and single cell ATAC-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 the AP-1 transcription factor family, that may represent key transcriptional nodes involved in PCD. Integration of transcriptomic and epigenetic data at the single cell level revealed that a population of preplasmablasts already undergone epigenetic remodeling related to PC profile together with UPR 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 plasma cell maturation and biological functions. These data thus provide critical insights into epigenetic- and transcriptional-mediated reprogramming events that sustain PC differentiation.

Hepatosplenic T-cell lymphoma displays an original oyster-shell cytological pattern and a distinct genomic profile from that of gamma-delta T-cell large granular lymphocytic leukemia.

Not available.

Regulus infers signed regulatory relations from few samples' information using discretization and likelihood constraints.

MOTIVATION: Transcriptional regulation is performed by transcription factors (TF) binding to DNA in context-dependent regulatory regions and determines the activation or inhibition of gene expression. Current methods of transcriptional regulatory circuits inference, based on one or all of TF, regions and genes activity measurements require a large number of samples for ranking the candidate TF-gene regulation relations and rarely predict whether they are activations or inhibitions. We hypothesize that transcriptional regulatory circuits can be inferred from fewer samples by (1) fully integrating information on TF binding, gene expression and regulatory regions accessibility, (2) reducing data complexity and (3) using biology-based likelihood constraints to determine the global consistency between a candidate TF-gene relation and patterns of genes expressions and region activations, as well as qualify regulations as activations or inhibitions. RESULTS: We introduce Regulus, a method which computes TF-gene relations from gene expressions, regulatory region activities and TF binding sites data, together with the genomic locations of all entities. After aggregating gene expressions and region activities into patterns, data are integrated into a RDF (Resource Description Framework) endpoint. A dedicated SPARQL (SPARQL Protocol and RDF Query Language) query retrieves all potential relations between expressed TF and genes involving active regulatory regions. These TF-region-gene relations are then filtered using biological likelihood constraints allowing to qualify them as activation or inhibition. Regulus provides signed relations consistent with public databases and, when applied to biological data, identifies both known and potential new regulators. Regulus is devoted to context-specific transcriptional circuits inference in human settings where samples are scarce and cell populations are closely related, using discretization into patterns and likelihood reasoning to decipher the most robust regulatory relations.

Five new F10 variants in hereditary factor x deficiency detected by high-throughput sequencing.

INTRODUCTION: Factor X deficiency is a rare inherited bleeding disorder. To date, 181 variants are reported in the recently updated F10-gene variant database. AIM: This study aimed to describe new F10 variants. METHOD: The F10 gene was analysed in 16 consecutive families with FX deficiency by a targeted high-throughput sequencing approach, including F10, F9, F8 genes, and 78 genes dedicated to haematological malignancies. RESULTS: We identified 19 variants (17 missense, one nonsense and one frameshift) and two copy number variations. Two patients presenting a combined FVII-FX deficiency showed a loss of one F10 gene copy (del13q34) associated with a missense variant on the remaining allele, leading to a FX:C significantly lower than the FVII:C level and explaining their unusual bleeding history. We reported five novel variants. Three missense variants (p.Glu22Val affecting the signal peptide cleavage site, p.Cys342Tyr removing the disulphide bond between the FX heavy and light chains, and p.Val385Met located in FX peptidase S1 domain) were detected at compound heterozygosis status in three patients with severe bleeding symptoms and FX:C level below 10 IU/dL. Two truncating variants p.Tyr279* and p.Thr434Aspfs*13 leading to an altered FX protein were found at heterozygous state in two patients with mild bleeding history. CONCLUSION: This study showed the feasibility and the interest of high-throughput sequencing approach for rare bleeding disorders, enabling the report of F10 gene screening in a 3-weeks delay, suitable for clinical use. The description of five new variants may contribute to a better understanding of the phenotype-genotype correlation in FX deficiency.

Final step of B-cell differentiation into plasmablasts; the right time to activate plasma cell PIM2 kinase.

The differentiation of B cells into antibody-secreting plasma cells is a complex process that involves extensive changes in morphology, lifespan, and cellular metabolism to support the high rates of antibody production. During the final stage of differentiation, B cells undergo significant expansion of their endoplasmic reticulum and mitochondria, which induces cellular stress and may lead to cell death in absence of effective inhibition of the apoptotic pathway. These changes are tightly regulated at transcriptional and epigenetic levels, as well as at post-translational level, with protein modifications playing a critical role in the process of cellular modification and adaptation. Our recent research has highlighted the pivotal role of the serine/threonine kinase PIM2 in B cell differentiation, from commitment stage to plasmablast and maintenance of expression in mature plasma cells. PIM2 has been shown to promote cell cycle progression during the final stage of differentiation and to inhibit Caspase 3 activation, raising the threshold for apoptosis. In this review, we examine the key molecular mechanisms controlled by PIM2 that contribute to plasma cell development and maintenance.

The Negative Influence of Baseline Cell-free DNA on Long-term Survival in DLBCL Depends on Frontline Treatment Intensity.

PURPOSE: This study aims to investigate the relationship between the intensity of the initial treatment given to patients with de novo diffuse large B-cell lymphoma (DLBCL) and the impact of their baseline cell-free DNA (cfDNA) levels on their long-term survival. EXPERIMENTAL DESIGN: The GOELAMS 075 randomized clinical trial compared rituximab plus cyclophosphamide, doxorubicin, vincristine and prednisone (R-CHOP) with high-dose R-chemotherapy plus autologous stem cell transplantation (R-HDT) for patients aged ≤60. An interim PET assessment was used to refer patients for salvage therapy. With a median follow-up of more than 5.8 years, we analyzed the effects of the treatment arm, salvage therapy, and cfDNA level at diagnosis on overall survival (OS). RESULTS: In a representative group of 123 patients, a high cfDNA concentration (>55 ng/mL) at diagnosis was associated with poor clinical prognostic factors and constituted a prognostic marker, independently of the age-adjusted International Prognostic Index. A cfDNA level above a threshold value of 55 ng/mL at diagnosis was associated with significantly worse OS. In an intention-to-treat analysis, high-cfDNA R-CHOP patients (but not high-cfDNA R-HDT patients) had worse OS [HR (95% confidence interval), 3.99 (1.98-10.74); P = 0.006]. In patients with high cfDNA levels, salvage therapy and transplantation were associated with a significantly higher OS rate. Among 50 patients with complete response 6 months after the end of treatment, for 11 of 24 R-CHOP patients, the cfDNA did not fall back to normal values. CONCLUSIONS: In this randomized clinical trial, intensive regimens mitigated the negative influence of high cfDNA levels in de novo DLBCL, relative to R-CHOP.

Human anti-smallpox long-lived memory B cells are defined by dynamic interactions in the splenic niche and long-lasting germinal center imprinting.

Memory B cells (MBCs) can persist for a lifetime, but the mechanisms that allow their long-term survival remain poorly understood. Here, we isolated and analyzed human splenic smallpox/vaccinia protein B5-specific MBCs in individuals who were vaccinated more than 40 years ago. Only a handful of clones persisted over such an extended period, and they displayed limited intra-clonal diversity with signs of extensive affinity-based selection. These long-lived MBCs appeared enriched in a CD21hiCD20hi IgG+ splenic B cell subset displaying a marginal-zone-like NOTCH/MYC-driven signature, but they did not harbor a unique longevity-associated transcriptional or metabolic profile. Finally, the telomeres of B5-specific, long-lived MBCs were longer than those in patient-paired naive B cells in all the samples analyzed. Overall, these results imply that separate mechanisms such as early telomere elongation, affinity selection during the contraction phase, and access to a specific niche contribute to ensuring the functional longevity of MBCs.

Improving reusability along the data life cycle: a regulatory circuits case study.

BACKGROUND: In life sciences, there has been a long-standing effort of standardization and integration of reference datasets and databases. Despite these efforts, many studies data are provided using specific and non-standard formats. This hampers the capacity to reuse the studies data in other pipelines, the capacity to reuse the pipelines results in other studies, and the capacity to enrich the data with additional information. The Regulatory Circuits project is one of the largest efforts for integrating human cell genomics data to predict tissue-specific transcription factor-genes interaction networks. In spite of its success, it exhibits the usual shortcomings limiting its update, its reuse (as a whole or partially), and its extension with new data samples. To address these limitations, the resource has previously been integrated in an RDF triplestore so that TF-gene interaction networks could be generated with two SPARQL queries. However, this triplestore did not store the computed networks and did not integrate metadata about tissues and samples, therefore limiting the reuse of this dataset. In particular, it does not enable to reuse only a portion of Regulatory Circuits if a study focuses on a subset of the tissues, nor to combine the samples described in the datasets with samples from other studies. Overall, these limitations advocate for the design of a complete, flexible and reusable representation of the Regulatory Circuits dataset based on Semantic Web technologies. RESULTS: We provide a modular RDF representation of the Regulatory Circuits, called Linked Extended Regulatory Circuits (LERC). It consists in (i) descriptions of biological and experimental context mapped to the references databases, (ii) annotations about TF-gene interactions at the sample level for 808 samples, (iii) annotations about TF-gene interactions at the tissue level for 394 tissues, (iv) metadata connecting the knowledge graphs cited above. LERC is based on a modular organisation into 1,205 RDF named graphs for representing the biological data, the sample-specific and the tissue-specific networks, and the corresponding metadata. In total it contains 3,910,794,050 triples and is available as a SPARQL endpoint. CONCLUSION: The flexible and modular architecture of LERC supports biologically-relevant SPARQL queries. It allows an easy and fast querying of the resources related to the initial Regulatory Circuits datasets and facilitates its reuse in other studies. ASSOCIATED WEBSITE: https://regulatorycircuits-lod.genouest.org.

PIM2 kinase has a pivotal role in plasmablast generation and plasma cell survival, opening up novel treatment options in myeloma.

The differentiation of B cells into plasmablasts (PBs) and then plasma cells (PCs) is associated with extensive cell reprogramming and new cell functions. By using specific inhibition strategies (including a novel morpholino RNA antisense approach), we found that early, sustained upregulation of the proviral integrations of Moloney virus 2 (PIM2) kinase is a pivotal event during human B-cell in vitro differentiation and then continues in mature normal and malignant PCs in the bone marrow. In particular, PIM2 sustained the G1/S transition by acting on CDC25A and p27Kip1 and hindering caspase 3-driven apoptosis through BAD phosphorylation and cytoplasmic stabilization of p21Cip1. In PCs, interleukin-6 triggered PIM2 expression, resulting in antiapoptotic effects on which malignant PCs were particularly dependent. In multiple myeloma, pan-PIM and myeloid cell leukemia-1 (MCL1) inhibitors displayed synergistic activity. Our results highlight a cell-autonomous function that links kinase activity to the newly acquired secretion ability of the PBs and the adaptability observed in both normal and malignant PCs. These findings should finally prompt the reconsideration of PIM2 as a therapeutic target in multiple myeloma.

Natural and human-driven selection of a single non-coding body size variant in ancient and modern canids.

Domestic dogs (Canis lupus familiaris) are the most variable-sized mammalian species on Earth, displaying a 40-fold size difference between breeds.1 Although dogs of variable size are found in the archeological record,2-4 the most dramatic shifts in body size are the result of selection over the last two centuries, as dog breeders selected and propagated phenotypic extremes within closed breeding populations.5 Analyses of over 200 domestic breeds have identified approximately 20 body size genes regulating insulin processing, fatty acid metabolism, TGFß signaling, and skeletal formation.6-10 Of these, insulin-like growth factor 1 (IGF1) predominates, controlling approximately 15% of body size variation between breeds.8 The identification of a functional mutation associated with IGF1 has thus far proven elusive.6,10,11 Here, to identify and elucidate the role of an ancestral IGF1 allele in the propagation of modern canids, we analyzed 1,431 genome sequences from 13 species, including both ancient and modern canids, thus allowing us to define the evolutionary history of both ancestral and derived alleles at this locus. We identified a single variant in an antisense long non-coding RNA (IGF1-AS) that interacts with the IGF1 gene, creating a duplex. While the derived mutation predominates in both modern gray wolves and large domestic breeds, the ancestral allele, which predisposes to small size, was common in small-sized breeds and smaller wild canids. Our analyses demonstrate that this major regulator of canid body size nearly vanished in Pleistocene wolves, before its recent resurgence resulting from human-imposed selection for small-sized breed dogs.

Committed Human CD23-Negative Light-Zone Germinal Center B Cells Delineate Transcriptional Program Supporting Plasma Cell Differentiation.

B cell affinity maturation occurs in the germinal center (GC). Light-zone (LZ) GC B cells (BGC-cells) interact with follicular dendritic cells (FDCs) and compete for the limited, sequential help from T follicular helper cells needed to escape from apoptosis and complete their differentiation. The highest-affinity LZ BGC-cells enter the cell cycle and differentiate into PCs, following a dramatic epigenetic reorganization that induces transcriptome changes in general and the expression of the PRDM1 gene in particular. Human PC precursors are characterized by the loss of IL-4/STAT6 signaling and the absence of CD23 expression. Here, we studied the fate of human LZ BGC-cells as a function of their CD23 expression. We first showed that CD23 expression was restricted to the GC LZ, where it was primarily expressed by FDCs; less than 10% of tonsil LZ BGC-cells were positive. Sorted LZ BGC-cells left in culture and stimulated upregulated CD23 expression but were unable to differentiate into PCs - in contrast to cells that did not upregulate CD23 expression. An in-depth analysis (including single-cell gene expression) showed that stimulated CD23-negative LZ BGC-cells differentiated into plasmablasts and time course of gene expression changes delineates the transcriptional program that sustains PC differentiation. In particular, we identified a B cell proliferation signature supported by a transient MYC gene expression. Overall, the CD23 marker might be of value in answering questions about the differentiation of normal BGC-cells and allowed us to propose an instructive LZ BGC-cells maturation and fate model.

Follicular lymphoma triggers phenotypic and functional remodeling of the human lymphoid stromal cell landscape.

Lymphoid stromal cells (LSCs) are essential organizers of immune responses. We analyzed tonsillar tissue by combining flow cytometry, in situ imaging, RNA sequencing, and functional assays, defining three distinct human LSC subsets. The integrin CD49a designated perivascular stromal cells exhibiting features of local committed LSC precursors and segregated cytokine and chemokine-producing fibroblastic reticular cells (FRCs) supporting B and T cell survival. The follicular dendritic cell transcriptional profile reflected active responses to B cell and non-B cell stimuli. We therefore examined the effect of B cell stimuli on LSCs in follicular lymphoma (FL). FL B cells interacted primarily with CD49a+ FRCs. Transcriptional analyses revealed LSC reprogramming in situ downstream of the cytokines tumor necrosis factor (TNF) and transforming growth factor ß (TGF-ß), including increased expression of the chemokines CCL19 and CCL21. Our findings define human LSC populations in healthy tissue and reveal bidirectional crosstalk between LSCs and malignant B cells that may present a targetable axis in lymphoma.

Early Emergence of Adaptive Mechanisms Sustaining Ig Production: Application to Antibody Therapy.

Antibody therapy, where artificially-produced immunoglobulins (Ig) are used to treat pathological conditions such as auto-immune diseases and cancers, is a very innovative and competitive field. Although substantial efforts have been made in recent years to obtain specific and efficient antibodies, there is still room for improvement especially when considering a precise tissular targeting or increasing antigen affinity. A better understanding of the cellular and molecular steps of terminal B cell differentiation, in which an antigen-activated B cell becomes an antibody secreting cell, may improve antibody therapy. In this review, we use our recently published data about human B cell differentiation, to show that the mechanisms necessary to adapt a metamorphosing B cell to its new secretory function appear quite early in the differentiation process i.e., at the pre-plasmablast stage. After characterizing the molecular pathways appearing at this stage, we will focus on recent findings about two main processes involved in antibody production: unfolded protein response (UPR) and endoplasmic reticulum (ER) stress. We'll show that many genes coding for factors involved in UPR and ER stress are induced at the pre-plasmablast stage, sustaining our hypothesis. Finally, we propose to use this recently acquired knowledge to improve productivity of industrialized therapeutic antibodies.

Extracellular vesicles shed by follicular lymphoma B cells promote polarization of the bone marrow stromal cell niche.

Follicular lymphoma (FL) originates in the lymph nodes (LNs) and infiltrates bone marrow (BM) early in the course of the disease. BM FL B cells are characterized by a lower cytological grade, decreased proliferation, and a specific phenotypic and subclonal profile. Mesenchymal stromal cells (MSCs) obtained from FL BM display a specific gene expression profile (GEP), including enrichment for a lymphoid stromal cell signature, and an increased capacity to sustain FL B-cell growth. However, the mechanisms triggering the formation of the medullar FL permissive stromal niche have not been identified. In the current work, we demonstrate that FL B cells produce extracellular vesicles (EVs) that can be internalized by BM-MSCs, making them more efficient to support FL B-cell survival and quiescence. Accordingly, EVs purified from FL BM plasma activate transforming growth factor ß-dependent and independent pathways in BM-MSCs and modify their GEP, triggering an upregulation of factors classically associated with hematopoietic stem cell niche, including CXCL12 and angiopoietin-1. Moreover, we provide the first characterization of BM FL B-cell GEP, allowing the definition of the landscape of molecular interactions they could engage with EV-primed BM-MSCs. This work identifies FL-derived EVs as putative mediators of BM stroma polarization and supports further investigation of their clinical interest for targeting the crosstalk between BM-MSCs and malignant B cells.

Lenalidomide triggers T-cell effector functions in vivo in patients with follicular lymphoma.

The immunomodulatory drug lenalidomide is used in patients with follicular lymphoma (FL) with the aim of stimulating T-cell antitumor immune response. However, little is known about the effects of lenalidomide on T-cell biology in vivo in patients with FL. We thus undertook an extensive longitudinal immunologic study, including phenotypic, transcriptomic, and functional analyses, on 44 first-line and 27 relapsed/refractory patients enrolled in the GALEN trial (Obinutuzumab Combined With Lenalidomide for Relapsed or Refractory Follicular B-Cell Lymphoma) to test the efficacy of lenalidomide and obinutuzumab combination in patients with FL. Lenalidomide rapidly and transiently induced an activated T-cell phenotype, including HLA-DR, Tim-3, CD137, and programmed cell death protein 1 (PD-1) upregulation. Furthermore, sequential RNA-sequencing of sorted PD-1+ and PD-1- T-cell subsets revealed that lenalidomide triggered a strong enrichment for several gene signatures related to effector memory T-cell features, including proliferation, antigen receptor signaling, and immune synapse restoration; all were validated at the phenotypic level and with ex vivo functional assays. Correlative analyses pinpointed a negative clinical impact of high effector T-cell and regulatory T-cell percentages before and during treatment. Our findings bring new insight in lenalidomide mechanisms of action at work in vivo and will fuel a new rationale for the design of combination therapies.

Functional characterization of PD1+TIM3+ tumor-infiltrating T cells in DLBCL and effects of PD1 or TIM3 blockade.

In diffuse large B-cell lymphoma (DLBCL), tumor-infiltrating T lymphocytes (TILs) are involved in therapeutic responses. However, tumor-specific TILs can be dysfunctional, with impaired effector functions. Various mechanisms are involved in this exhaustion, and the increased expression of programmed cell death receptor 1 (PD1) and TIM3 on dysfunctional cells suggests their involvement. However, conflicting data have been published regarding their expression or coexpression in DLBCL. We evaluated the presence and phenotype of CD4+ and CD8+ TILs in freshly collected tumor tissues in DLBCL and compared the results with those in follicular lymphoma, classical Hodgkin lymphoma, and nonmalignant reactive lymphadenopathy. We found that TILs expressing both PD1 and TIM3 were expanded in DLBCL, particularly in the activated B cell-like subgroup. Isolated PD1+TIM3+ TILs exhibited a transcriptomic signature related to T-cell exhaustion associated with a reduction in cytokine production, both compromising the antitumor immune response. However, these cells expressed high levels of cytotoxic molecules. In line with this, stimulated PD1+TIM3+ TILs from DLBCL patients exhibited reduced proliferation and impaired secretion of interferon-γ, but these functions were restored by the blockade of PD1 or TIM3. In summary, the PD1+TIM3+ TIL population is expanded and exhausted in DLBCL but can be reinvigorated with appropriate therapies.