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.

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.

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.

Plasmablasts derive from CD23- activated B cells after the extinction of IL-4/STAT6 signaling and IRF4 induction.

The terminal differentiation of B cells into antibody-secreting cells (ASCs) is a critical component of adaptive immune responses. However, it is a very sensitive process, and dysfunctions lead to a variety of lymphoproliferative neoplasias including germinal center-derived lymphomas. To better characterize the late genomic events that drive the ASC differentiation of human primary naive B cells, we used our in vitro differentiation system and a combination of RNA sequencing and Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC sequencing). We discovered 2 mechanisms that drive human terminal B-cell differentiation. First, after an initial response to interleukin-4 (IL-4), cells that were committed to an ASC fate downregulated the CD23 marker and IL-4 signaling, whereas cells that maintained IL-4 signaling did not differentiate. Second, human CD23- cells also increased IRF4 protein to levels required for ASC differentiation, but they did that independently of the ubiquitin-mediated degradation process previously described in mice. Finally, we showed that CD23- cells carried the imprint of their previous activated B-cell status, were precursors of plasmablasts, and had a phenotype similar to that of in vivo preplasmablasts. Altogether, our results provide an unprecedented genomic characterization of the fate decision between activated B cells and plasmablasts, which provides new insights into the pathological mechanisms that drive lymphoma biology.

Locus suicide recombination actively occurs on the functionally rearranged IgH allele in B-cells from inflamed human lymphoid tissues.

B-cell activation yields abundant cell death in parallel to clonal amplification and remodeling of immunoglobulin (Ig) genes by activation-induced deaminase (AID). AID promotes affinity maturation of Ig variable regions and class switch recombination (CSR) in mature B lymphocytes. In the IgH locus, these processes are under control of the 3' regulatory region (3'RR) super-enhancer, a region demonstrated in the mouse to be both transcribed and itself targeted by AID-mediated recombination. Alternatively to CSR, IgH deletions joining Sµ to "like-switch" DNA repeats that flank the 3' super-enhancer can thus accomplish so-called "locus suicide recombination" (LSR) in mouse B-cells. Using an optimized LSR-seq high throughput method, we now show that AID-mediated LSR is evolutionarily conserved and also actively occurs in humans, providing an activation-induced cell death pathway in multiple conditions of B-cell activation. LSR either focuses on the functional IgH allele or is bi-allelic, and its signature is mainly detected when LSR is ongoing while it vanishes from fully differentiated plasma cells or from "resting" blood memory B-cells. Highly diversified breakpoints are distributed either within the upstream (3'RR1) or downstream (3'RR2) copies of the IgH 3' super-enhancer and all conditions activating CSR in vitro also seem to trigger LSR although TLR ligation appeared the most efficient. Molecular analysis of breakpoints and junctions confirms that LSR is AID-dependent and reveals junctional sequences somehow similar to CSR junctions but with increased usage of microhomologies.

The hydroxymethylome of multiple myeloma identifies FAM72D as a 1q21 marker linked to proliferation.

Cell identity relies on the cross-talk between genetics and epigenetics and their impact on gene expression. Oxidation of 5-methylcytosine (5mC) into 5-hydroxymethylcytosine (5hmC) is the first step of an active DNA demethylation process occurring mainly at enhancers and gene bodies and, as such, participates in processes governing cell identity in normal and pathological conditions. Although genetic alterations are well documented in multiple myeloma (MM), epigenetic alterations associated with this disease have not yet been thoroughly analyzed. To gain insight into the biology of MM, genome-wide 5hmC profiles were obtained and showed that regions enriched in this modified base overlap with MM enhancers and super enhancers and are close to highly expressed genes. Through the definition of a MM-specific 5hmC signature, we identified FAM72D as a poor prognostic gene located on 1q21, a region amplified in high risk myeloma. We further uncovered that FAM72D functions as part of the FOXM1 transcription factor network controlling cell proliferation and survival and we evidenced an increased sensitivity of cells expressing high levels of FOXM1 and FAM72 to epigenetic drugs targeting histone deacetylases and DNA methyltransferases.

Purkinje cells and Bergmann glia are primary targets of the TRα1 thyroid hormone receptor during mouse cerebellum postnatal development.

Thyroid hormone is necessary for normal development of the central nervous system, as shown by the severe mental retardation syndrome affecting hypothyroid patients with low levels of active thyroid hormone. The postnatal defects observed in hypothyroid mouse cerebellum are recapitulated in mice heterozygous for a dominant-negative mutation of Thra, the gene encoding the ubiquitous TRα1 receptor. Using CRE/loxP-mediated conditional expression approach, we found that this mutation primarily alters the differentiation of Purkinje cells and Bergmann glia, two cerebellum-specific cell types. These primary defects indirectly affect cerebellum development in a global manner. Notably, the inward migration and terminal differentiation of granule cell precursors is impaired. Therefore, despite the broad distribution of its receptors, thyroid hormone targets few cell types that exert a predominant role in the network of cellular interactions that govern normal cerebellum maturation.

Genome-wide analysis of thyroid hormone receptors shared and specific functions in neural cells.

TRα1 and TRß1, the two main thyroid hormone receptors in mammals, are transcription factors that share similar properties. However, their respective functions are very different. This functional divergence might be explained in two ways: it can reflect different expression patterns or result from different intrinsic properties of the receptors. We tested this second hypothesis by comparing the repertoires of 3,3',5-triiodo-L-thyronine (T3)-responsive genes of two neural cell lines, expressing either TRα1 or TRß1. Using transcriptome analysis, we found that a substantial fraction of the T3 target genes display a marked preference for one of the two receptors. So when placed alone in identical situations, the two receptors have different repertoires of target genes. Chromatin occupancy analysis, performed at a genome-wide scale, revealed that TRα1 and TRß1 cistromes were also different. However, receptor-selective regulation of T3 target genes did not result from receptor-selective chromatin occupancy of their promoter regions. We conclude that modification of TRα1 and TRß1 intrinsic properties contributes in a large part to the divergent evolution of the receptors' function, at least during neurodevelopment.

Thyroid hormone triggers the developmental loss of axonal regenerative capacity via thyroid hormone receptor α1 and krüppel-like factor 9 in Purkinje cells.

Neurons in the CNS of higher vertebrates lose their ability to regenerate their axons at a stage of development that coincides with peak circulating thyroid hormone (T(3)) levels. Here, we examined whether this peak in T(3) is involved in the loss of axonal regenerative capacity in Purkinje cells (PCs). This event occurs at the end of the first postnatal week in mice. Using organotypic culture, we found that the loss of axon regenerative capacity was triggered prematurely by early exposure of mouse PCs to T(3), whereas it was delayed in the absence of T(3). Analysis of mutant mice showed that this effect was mainly mediated by the T(3) receptor α1. Using gain- and loss-of-function approaches, we also showed that Krüppel-like factor 9 was a key mediator of this effect of T(3). These results indicate that the sudden physiological increase in T(3) during development is involved in the onset of the loss of axon regenerative capacity in PCs. This loss of regenerative capacity might be part of the general program triggered by T(3) throughout the body, which adapts the animal to its postnatal environment.

Direct and indirect consequences on gene expression of a thyroid hormone receptor alpha 1 mutation restricted to Sertoli cells.

Thyroid hormone is required for the timely transition of Sertoli cells from proliferative to differentiating and maturing. This transition takes place during a critical developmental period in mammals, which in mice is the first post-natal week. In order to identify the underlying molecular mechanisms of this differentiation process, we used Cre/loxP technology to selectively block the function of the thyroid hormone receptor TRα1 in Sertoli cells. We then used RNA-seq to analyze the changes in gene expression induced in the post-natal testis. This differential analysis provides genetic clues to the initial testicular defects resulting from disrupted thyroid hormone signaling, and suggests that Sertoli cells influence germ cells soon after their birth.

Severe impairment of cerebellum development in mice expressing a dominant-negative mutation inactivating thyroid hormone receptor alpha1 isoform.

Thyroid hormone deficiency is known to deeply affect cerebellum post-natal development. We present here a detailed analysis of the phenotype of a recently generated mouse model, expressing a dominant-negative TRα1 mutation. Although hormonal level is not affected, the cerebellum of these mice displays profound alterations in neuronal and glial differentiation, which are reminiscent of congenital hypothyroidism, indicating a predominant function of this receptor isoform in normal cerebellum development. Some of the observed effects might result from the cell autonomous action of the mutation, while others are more likely to result from a reduction in neurotrophic factor production.

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.

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.

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.

Integrative Analysis of Cell Crosstalk within Follicular Lymphoma Cell Niche: Towards a Definition of the FL Supportive Synapse.

Follicular lymphoma (FL), the most frequent indolent non-Hodgkin's B cell lymphoma, is considered as a prototypical centrocyte-derived lymphoma, dependent on a specific microenvironment mimicking the normal germinal center (GC). In agreement, several FL genetic alterations affect the crosstalk between malignant B cells and surrounding cells, including stromal cells and follicular helper T cells (Tfh). In our study, we sought to deconvolute this complex FL supportive synapse by comparing the transcriptomic profiles of GC B cells, Tfh, and stromal cells, isolated from normal versus FL tissues, in order to identify tumor-specific pathways. In particular, we highlighted a high expression of IL-6 and IL-7 in FL B cells that could favor the activation of FL Tfh overexpressing IFNG, able in turn to stimulate FL B cells without triggering MHC (major histocompatibility) class II expression. Moreover, the glycoprotein clusterin was found up-regulated in FL stromal cells and could promote FL B cell adhesion. Finally, besides its expression on Tfh, CD200 was found overexpressed on tumor B cells and could contribute to the induction of the immunosuppressive enzyme indoleamine-2,3 dioxygenase by CD200R-expressing dendritic cells. Altogether our findings led us to outline the contribution of major signals provided by the FL microenvironment and their interactions with malignant FL B cells.

Cloning and analysis of Nkx6.3 during CNS and gastrointestinal development.

Members of the Nkx family of homeodomain proteins are involved in a variety of developmental processes such as cell fate determination in the CNS and in the pancreas. Here we describe the cloning and developmental expression pattern of Nkx6.3, a new member of the Nkx6 subfamily of homeodomain proteins. Nkx6.3 is expressed in the developing CNS and gastro-intestinal tract. In contrast to Nkx6.1 and Nkx6.2 that are broadly expressed in ventral positions of the developing CNS, Nkx6.3 shows a remarkably selective expression in a subpopulation of differentiating V2 neurons at caudal hindbrain levels. The expression of Nkx6.3 at this level depends on the activity of other Nkx6 proteins. In the gut, Nkx6.3 is expressed in duodenal and glandular stomach endoderm and at the end of gestation Nkx6.3 became restricted to the base of the gastric units in the glandular stomach. The expression of Nkx6.3 overlapped with the expression of Nkx6.2 both in the CNS and in the gut. Transient Nkx6.2 expression was also detected in the developing pancreas. However, analysis of Nkx6.2(-/-) mice did not display any obvious aberrations of pancreatic or stomach development.

Cell-Cycle-Dependent Reconfiguration of the DNA Methylome during Terminal Differentiation of Human B Cells into Plasma Cells.

Molecular mechanisms underlying terminal differentiation of B cells into plasma cells are major determinants of adaptive immunity but remain only partially understood. Here we present the transcriptional and epigenomic landscapes of cell subsets arising from activation of human naive B cells and differentiation into plasmablasts. Cell proliferation of activated B cells was linked to a slight decrease in DNA methylation levels, but followed by a committal step in which an S phase-synchronized differentiation switch was associated with an extensive DNA demethylation and local acquisition of 5-hydroxymethylcytosine at enhancers and genes related to plasma cell identity. Downregulation of both TGF-?1/SMAD3 signaling and p53 pathway supported this final step, allowing the emergence of a CD23-negative subpopulation in transition from B cells to plasma cells. Remarkably, hydroxymethylation of PRDM1, a gene essential for plasma cell fate, was coupled to progression in S phase, revealing an intricate connection among cell cycle, DNA (hydroxy)methylation, and cell fate determination.

From hindbrain segmentation to breathing after birth: developmental patterning in rhombomeres 3 and 4.

Respiration is a rhythmic motor behavior that appears in the fetus and acquires a vital importance at birth. It is generated within central pattern-generating neuronal networks of the hindbrain. This region of the brain is of particular interest since it is the most understood part with respect to the cellular and molecular mechanisms that underlie its development. Hox paralogs and Hox-regulating genes kreisler/mafB and Krox20 are required for the normal formation of rhombomeres in vertebrate embryos. From studies of rhombomeres r3 and r4, the authors review mechanisms whereby these developmental genes may govern the early embryonic development of para-facial neuronal networks and specify patterns of motor activities operating throughout life. A model whereby the regional identity of progenitor cells can be abnormally specified in r3 and r4 after a mutation of these genes is proposed. Novel neuronal circuits may develop from some of these misspecified progenitors while others are eliminated, eventually affecting respiration and survival after birth.

Early development of respiratory rhythm generation in mouse and chick.

We are investigating neuronal circuits resulting from conservative developmental mechanisms orchestrating the segmentation of the vertebrates hindbrain into compartments called rhombomeres (r). Segmentation transcription factors Hoxa1, Krox20 and kreisler are expressed in the future rhombomeres r4-r5, r3 and r5, r5-r6, respectively. In mice, the in vivo and in vitro analysis of neuronal groups after inactivation of these three genes revealed distinct postnatal respiratory phenotypes associated with defects of central respiratory controls resulting from deletion, neoformation or reconfiguration of modular circuits. In chick and mice, we have found neuronal rhythm generators that conform to the rhombomeric anatomical pattern as early as at the end of the segmentation. By isolating chick hindbrain segments in vitro, we have also identified rhombomeric motifs allowing the formation or deletion of a specific (GABAergic) rhythm-promoting module. Therefore, primordial rhombomeric organization of the hindbrain seems to determine a modular organization of the rhythmogenic network, thereby influencing later function of brainstem respiratory control networks.