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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Toxicogenomic analysis of the ability of brominated flame retardants TBBPA and BDE-209 to disrupt thyroid hormone signaling in neural cells.

Brominated flame retardants are suspected to act as disruptors of thyroid hormone signaling. This raises the concern that they might affect children's cognitive functions by influencing thyroid hormone signaling in the developing brain. We present here an in vitro analysis of the ability of the most common compounds, tetrabromobisphenol A (TBBPA) and BDE-209, to alter thyroid hormone response based on a model neural cell line and genome-wide analysis of gene expression.

Deciphering direct and indirect influence of thyroid hormone with mouse genetics.

T3, the active form of thyroid hormone, binds nuclear receptors that regulate the transcription of a large number of genes in many cell types. Unraveling the direct and indirect effect of this hormonal stimulation, and establishing links between these molecular events and the developmental and physiological functions of the hormone, is a major challenge. New mouse genetics tools, notably those based on Cre/loxP technology, are suitable to perform a multiscale analysis of T3 signaling and achieve this task.

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.

Genome-wide search reveals the existence of a limited number of thyroid hormone receptor alpha target genes in cerebellar neurons.

Thyroid hormone (T3) has a major influence on cerebellum post-natal development. The major phenotypic landmark of exposure to low levels of T3 during development (hypothyroidism) in the cerebellum is the retarded inward migration of the most numerous cell type, granular neurons. In order to identify the direct genetic regulation exerted by T3 on cerebellar neurons and their precursors, we used microarray RNA hybridization to perform a time course analysis of T3 induced gene expression in primary cultures of cerebellar neuronal cell. These experiments suggest that we identified a small set of genes which are directly regulated, both in vivo and in vitro, during cerebellum post-natal development. These modest changes suggest that T3 does not acts directly on granular neurons and mainly indirectly influences the cellular interactions taking place during development.

A bimodal influence of thyroid hormone on cerebellum oligodendrocyte differentiation.

Thyroid hormone (T(3)) can trigger a massive differentiation of cultured oligodendrocytes precursor cells (OPC) by binding the nuclear T(3) receptor α1 (TRα1). Whether this reflects a physiological function of TRα1 remains uncertain. Using a recently generated mouse model, in which CRE/loxP recombination is used to block its function, we show that TRα1 acts at two levels for the in vivo differentiation of OPC in mouse cerebellum. At the early postnatal stage, it promotes the secretion of several neurotrophic factors by acting in Purkinje neurons and astrocytes, defining an environment suitable for OPC differentiation. At later stages, TRα1 acts in a cell-autonomous manner to ensure the complete arrest of OPC proliferation. These data explain contradictory observations made on various models and outline the importance of T(3) signaling both for synchronizing postnatal neurodevelopment and restraining OPC proliferation in adult brain.