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.

Alternative splicing induced by bacterial pore-forming toxins sharpens CIRBP-mediated cell response to Listeria infection.

Cell autonomous responses to intracellular bacteria largely depend on reorganization of gene expression. To gain isoform-level resolution of these modes of regulation, we combined long- and short-read transcriptomic analyses of the response of intestinal epithelial cells to infection by the foodborne pathogen Listeria monocytogenes. Among the most striking isoform-based types of regulation, expression of the cellular stress response regulator CIRBP (cold-inducible RNA-binding protein) and of several SRSFs (serine/arginine-rich splicing factors) switched from canonical transcripts to nonsense-mediated decay-sensitive isoforms by inclusion of 'poison exons'. We showed that damage to host cell membranes caused by bacterial pore-forming toxins (listeriolysin O, perfringolysin, streptolysin or aerolysin) led to the dephosphorylation of SRSFs via the inhibition of the kinase activity of CLK1, thereby driving CIRBP alternative splicing. CIRBP isoform usage was found to have consequences on infection, since selective repression of canonical CIRBP reduced intracellular bacterial load while that of the poison exon-containing isoform exacerbated it. Consistently, CIRBP-bound mRNAs were shifted towards stress-relevant transcripts in infected cells, with increased mRNA levels or reduced translation efficiency for some targets. Our results thus generalize the alternative splicing of CIRBP and SRSFs as a common response to biotic or abiotic stresses by extending its relevance to the context of bacterial infection.

Cis-acting variation is common across regulatory layers but is often buffered during embryonic development.

Precise patterns of gene expression are driven by interactions between transcription factors, regulatory DNA sequences, and chromatin. How DNA mutations affecting any one of these regulatory "layers" are buffered or propagated to gene expression remains unclear. To address this, we quantified allele-specific changes in chromatin accessibility, histone modifications, and gene expression in F1 embryos generated from eight Drosophila crosses at three embryonic stages, yielding a comprehensive data set of 240 samples spanning multiple regulatory layers. Genetic variation (allelic imbalance) impacts gene expression more frequently than chromatin features, with metabolic and environmental response genes being most often affected. Allelic imbalance in cis-regulatory elements (enhancers) is common and highly heritable, yet its functional impact does not generally propagate to gene expression. When it does, genetic variation impacts RNA levels through two alternative mechanisms involving either H3K4me3 or chromatin accessibility and H3K27ac. Changes in RNA are more predictive of variation in H3K4me3 than vice versa, suggesting a role for H3K4me3 downstream from transcription. The impact of a substantial proportion of genetic variation is consistent across embryonic stages, with 50% of allelic imbalanced features at one stage being also imbalanced at subsequent developmental stages. Crucially, buffering, as well as the magnitude and evolutionary impact of genetic variants, is influenced by regulatory complexity (i.e., number of enhancers regulating a gene), with transcription factors being most robust to cis-acting, but most influenced by trans-acting, variation.

Deciphering and modelling the TGF-ß signalling interplays specifying the dorsal-ventral axis of the sea urchin embryo.

During sea urchin development, secretion of Nodal and BMP2/4 ligands and their antagonists Lefty and Chordin from a ventral organiser region specifies the ventral and dorsal territories. This process relies on a complex interplay between the Nodal and BMP pathways through numerous regulatory circuits. To decipher the interplay between these pathways, we used a combination of treatments with recombinant Nodal and BMP2/4 proteins and a computational modelling approach. We assembled a logical model focusing on cell responses to signalling inputs along the dorsal-ventral axis, which was extended to cover ligand diffusion and enable multicellular simulations. Our model simulations accurately recapitulate gene expression in wild-type embryos, accounting for the specification of ventral ectoderm, ciliary band and dorsal ectoderm. Our model simulations further recapitulate various morphant phenotypes, reveal a dominance of the BMP pathway over the Nodal pathway and stress the crucial impact of the rate of Smad activation in dorsal-ventral patterning. These results emphasise the key role of the mutual antagonism between the Nodal and BMP2/4 pathways in driving early dorsal-ventral patterning of the sea urchin embryo.

RSAT 2022: regulatory sequence analysis tools.

RSAT (Regulatory Sequence Analysis Tools) enables the detection and the analysis of cis-regulatory elements in genomic sequences. This software suite performs (i) de novo motif discovery (including from genome-wide datasets like ChIP-seq/ATAC-seq) (ii) genomic sequences scanning with known motifs, (iii) motif analysis (quality assessment, comparisons and clustering), (iv) analysis of regulatory variations and (v) comparative genomics. RSAT comprises 50 tools. Six public Web servers (including a teaching server) are offered to meet the needs of different biological communities. RSAT philosophy and originality are: (i) a multi-modal access depending on the user needs, through web forms, command-line for local installation and programmatic web services, (ii) a support for virtually any genome (animals, bacteria, plants, totalizing over 10 000 genomes directly accessible). Since the 2018 NAR Web Software Issue, we have developed a large REST API, extended the support for additional genomes and external motif collections, enhanced some tools and Web forms, and developed a novel tool that builds or refine gene regulatory networks using motif scanning (network-interactions). The RSAT website provides extensive documentation, tutorials and published protocols. RSAT code is under open-source license and now hosted in GitHub. RSAT is available at http://www.rsat.eu/.

Synthetic STARR-seq reveals how DNA shape and sequence modulate transcriptional output and noise.

The binding of transcription factors to short recognition sequences plays a pivotal role in controlling the expression of genes. The sequence and shape characteristics of binding sites influence DNA binding specificity and have also been implicated in modulating the activity of transcription factors downstream of binding. To quantitatively assess the transcriptional activity of tens of thousands of designed synthetic sites in parallel, we developed a synthetic version of STARR-seq (synSTARR-seq). We used the approach to systematically analyze how variations in the recognition sequence of the glucocorticoid receptor (GR) affect transcriptional regulation. Our approach resulted in the identification of a novel highly active functional GR binding sequence and revealed that sequence variation both within and flanking GR's core binding site can modulate GR activity without apparent changes in DNA binding affinity. Notably, we found that the sequence composition of variants with similar activity profiles was highly diverse. In contrast, groups of variants with similar activity profiles showed specific DNA shape characteristics indicating that DNA shape may be a better predictor of activity than DNA sequence. Finally, using single cell experiments with individual enhancer variants, we obtained clues indicating that the architecture of the response element can independently tune expression mean and cell-to cell variability in gene expression (noise). Together, our studies establish synSTARR as a powerful method to systematically study how DNA sequence and shape modulate transcriptional output and noise.

RSAT 2018: regulatory sequence analysis tools 20th anniversary.

RSAT (Regulatory Sequence Analysis Tools) is a suite of modular tools for the detection and the analysis of cis-regulatory elements in genome sequences. Its main applications are (i) motif discovery, including from genome-wide datasets like ChIP-seq/ATAC-seq, (ii) motif scanning, (iii) motif analysis (quality assessment, comparisons and clustering), (iv) analysis of regulatory variations, (v) comparative genomics. Six public servers jointly support 10 000 genomes from all kingdoms. Six novel or refactored programs have been added since the 2015 NAR Web Software Issue, including updated programs to analyse regulatory variants (retrieve-variation-seq, variation-scan, convert-variations), along with tools to extract sequences from a list of coordinates (retrieve-seq-bed), to select motifs from motif collections (retrieve-matrix), and to extract orthologs based on Ensembl Compara (get-orthologs-compara). Three use cases illustrate the integration of new and refactored tools to the suite. This Anniversary update gives a 20-year perspective on the software suite. RSAT is well-documented and available through Web sites, SOAP/WSDL (Simple Object Access Protocol/Web Services Description Language) web services, virtual machines and stand-alone programs at http://www.rsat.eu/.

Krox20 hindbrain regulation incorporates multiple modes of cooperation between cis-acting elements.

Developmental genes can harbour multiple transcriptional enhancers that act simultaneously or in succession to achieve robust and precise spatiotemporal expression. However, the mechanisms underlying cooperation between cis-acting elements are poorly documented, notably in vertebrates. The mouse gene Krox20 encodes a transcription factor required for the specification of two segments (rhombomeres) of the developing hindbrain. In rhombomere 3, Krox20 is subject to direct positive feedback governed by an autoregulatory enhancer, element A. In contrast, a second enhancer, element C, distant by 70 kb, is active from the initiation of transcription independent of the presence of the KROX20 protein. Here, using both enhancer knock-outs and investigations of chromatin organisation, we show that element C possesses a dual activity: besides its classical enhancer function, it is also permanently required in cis to potentiate the autoregulatory activity of element A, by increasing its chromatin accessibility. This work uncovers a novel, asymmetrical, long-range mode of cooperation between cis-acting elements that might be essential to avoid promiscuous activation of positive autoregulatory elements.

Logical modeling of lymphoid and myeloid cell specification and transdifferentiation.

Blood cells are derived from a common set of hematopoietic stem cells, which differentiate into more specific progenitors of the myeloid and lymphoid lineages, ultimately leading to differentiated cells. This developmental process is controlled by a complex regulatory network involving cytokines and their receptors, transcription factors, and chromatin remodelers. Using public data and data from our own molecular genetic experiments (quantitative PCR, Western blot, EMSA) or genome-wide assays (RNA-sequencing, ChIP-sequencing), we have assembled a comprehensive regulatory network encompassing the main transcription factors and signaling components involved in myeloid and lymphoid development. Focusing on B-cell and macrophage development, we defined a qualitative dynamical model recapitulating cytokine-induced differentiation of common progenitors, the effect of various reported gene knockdowns, and the reprogramming of pre-B cells into macrophages induced by the ectopic expression of specific transcription factors. The resulting network model can be used as a template for the integration of new hematopoietic differentiation and transdifferentiation data to foster our understanding of lymphoid/myeloid cell-fate decisions.

RSAT matrix-clustering: dynamic exploration and redundancy reduction of transcription factor binding motif collections.

Transcription factor (TF) databases contain multitudes of binding motifs (TFBMs) from various sources, from which non-redundant collections are derived by manual curation. The advent of high-throughput methods stimulated the production of novel collections with increasing numbers of motifs. Meta-databases, built by merging these collections, contain redundant versions, because available tools are not suited to automatically identify and explore biologically relevant clusters among thousands of motifs. Motif discovery from genome-scale data sets (e.g. ChIP-seq) also produces redundant motifs, hampering the interpretation of results. We present matrix-clustering, a versatile tool that clusters similar TFBMs into multiple trees, and automatically creates non-redundant TFBM collections. A feature unique to matrix-clustering is its dynamic visualisation of aligned TFBMs, and its capability to simultaneously treat multiple collections from various sources. We demonstrate that matrix-clustering considerably simplifies the interpretation of combined results from multiple motif discovery tools, and highlights biologically relevant variations of similar motifs. We also ran a large-scale application to cluster ∼11 000 motifs from 24 entire databases, showing that matrix-clustering correctly groups motifs belonging to the same TF families, and drastically reduced motif redundancy. matrix-clustering is integrated within the RSAT suite (http://rsat.eu/), accessible through a user-friendly web interface or command-line for its integration in pipelines.

Role of the chromatin landscape and sequence in determining cell type-specific genomic glucocorticoid receptor binding and gene regulation.

The genomic loci bound by the glucocorticoid receptor (GR), a hormone-activated transcription factor, show little overlap between cell types. To study the role of chromatin and sequence in specifying where GR binds, we used Bayesian modeling within the universe of accessible chromatin. Taken together, our results uncovered that although GR preferentially binds accessible chromatin, its binding is biased against accessible chromatin located at promoter regions. This bias can only be explained partially by the presence of fewer GR recognition sequences, arguing for the existence of additional mechanisms that interfere with GR binding at promoters. Therefore, we tested the role of H3K9ac, the chromatin feature with the strongest negative association with GR binding, but found that this correlation does not reflect a causative link. Finally, we find a higher percentage of promoter-proximal GR binding for genes regulated by GR across cell types than for cell type-specific target genes. Given that GR almost exclusively binds accessible chromatin, we propose that cell type-specific regulation by GR preferentially occurs via distal enhancers, whose chromatin accessibility is typically cell type-specific, whereas ubiquitous target gene regulation is more likely to result from binding to promoter regions, which are often accessible regardless of cell type examined.

ChIP-exo signal associated with DNA-binding motifs provides insight into the genomic binding of the glucocorticoid receptor and cooperating transcription factors.

The classical DNA recognition sequence of the glucocorticoid receptor (GR) appears to be present at only a fraction of bound genomic regions. To identify sequences responsible for recruitment of this transcription factor (TF) to individual loci, we turned to the high-resolution ChIP-exo approach. We exploited this signal by determining footprint profiles of TF binding at single-base-pair resolution using ExoProfiler, a computational pipeline based on DNA binding motifs. When applied to our GR and the few available public ChIP-exo data sets, we find that ChIP-exo footprints are protein- and recognition sequence-specific signatures of genomic TF association. Furthermore, we show that ChIP-exo captures information about TFs other than the one directly targeted by the antibody in the ChIP procedure. Consequently, the shape of the ChIP-exo footprint can be used to discriminate between direct and indirect (tethering to other DNA-bound proteins) DNA association of GR. Together, our findings indicate that the absence of classical recognition sequences can be explained by direct GR binding to a broader spectrum of sequences than previously known, either as a homodimer or as a heterodimer binding together with a member of the ETS or TEAD families of TFs, or alternatively by indirect recruitment via FOX or STAT proteins. ChIP-exo footprints also bring structural insights and locate DNA:protein cross-link points that are compatible with crystal structures of the studied TFs. Overall, our generically applicable footprint-based approach uncovers new structural and functional insights into the diverse ways of genomic cooperation and association of TFs.

Identification and characterization of DNA sequences that prevent glucocorticoid receptor binding to nearby response elements.

Out of the myriad of potential DNA binding sites of the glucocorticoid receptor (GR) found in the human genome, only a cell-type specific minority is actually bound, indicating that the presence of a recognition sequence alone is insufficient to specify where GR binds. Cooperative interactions with other transcription factors (TFs) are known to contribute to binding specificity. Here, we reasoned that sequence signals preventing GR recruitment to certain loci provide an alternative means to confer specificity. Motif analyses uncovered candidate Negative Regulatory Sequences (NRSs) that interfere with genomic GR binding. Subsequent functional analyses demonstrated that NRSs indeed prevent GR binding to nearby response elements. We show that NRS activity is conserved across species, found in most tissues and that they also interfere with the genomic binding of other TFs. Interestingly, the effects of NRSs appear not to be a simple consequence of changes in chromatin accessibility. Instead, we find that NRSs interact with proteins found at sub-nuclear structures called paraspeckles and that these proteins might mediate the repressive effects of NRSs. Together, our studies suggest that the joint influence of positive and negative sequence signals partition the genome into regions where GR can bind and those where it cannot.

Histone Chaperone SSRP1 is Essential for Wnt Signaling Pathway Activity During Osteoblast Differentiation.

Cellular differentiation is accompanied by dramatic changes in chromatin structure which direct the activation of lineage-specific transcriptional programs. Structure-specific recognition protein-1 (SSRP1) is a histone chaperone which is important for chromatin-associated processes such as transcription, DNA replication and repair. Since the function of SSRP1 during cell differentiation remains unclear, we investigated its potential role in controlling lineage determination. Depletion of SSRP1 in human mesenchymal stem cells elicited lineage-specific effects by increasing expression of adipocyte-specific genes and decreasing the expression of osteoblast-specific genes. Consistent with a role in controlling lineage specification, transcriptome-wide RNA-sequencing following SSRP1 depletion and the induction of osteoblast differentiation revealed a specific decrease in the expression of genes involved in biological processes related to osteoblast differentiation. Importantly, we observed a specific downregulation of target genes of the canonical Wnt signaling pathway, which was accompanied by decreased nuclear localization of active ß-catenin. Together our data uncover a previously unknown role for SSRP1 in promoting the activation of the Wnt signaling pathway activity during cellular differentiation. Stem Cells 2016;34:1369-1376.

RSAT 2015: Regulatory Sequence Analysis Tools.

RSAT (Regulatory Sequence Analysis Tools) is a modular software suite for the analysis of cis-regulatory elements in genome sequences. Its main applications are (i) motif discovery, appropriate to genome-wide data sets like ChIP-seq, (ii) transcription factor binding motif analysis (quality assessment, comparisons and clustering), (iii) comparative genomics and (iv) analysis of regulatory variations. Nine new programs have been added to the 43 described in the 2011 NAR Web Software Issue, including a tool to extract sequences from a list of coordinates (fetch-sequences from UCSC), novel programs dedicated to the analysis of regulatory variants from GWAS or population genomics (retrieve-variation-seq and variation-scan), a program to cluster motifs and visualize the similarities as trees (matrix-clustering). To deal with the drastic increase of sequenced genomes, RSAT public sites have been reorganized into taxon-specific servers. The suite is well-documented with tutorials and published protocols. The software suite is available through Web sites, SOAP/WSDL Web services, virtual machines and stand-alone programs at http://www.rsat.eu/.

A naturally occurring insertion of a single amino acid rewires transcriptional regulation by glucocorticoid receptor isoforms.

In addition to guiding proteins to defined genomic loci, DNA can act as an allosteric ligand that influences protein structure and activity. Here we compared genome-wide binding, transcriptional regulation, and, using NMR, the conformation of two glucocorticoid receptor (GR) isoforms that differ by a single amino acid insertion in the lever arm, a domain that adopts DNA sequence-specific conformations. We show that these isoforms differentially regulate gene expression levels through two mechanisms: differential DNA binding and altered communication between GR domains. Our studies suggest a versatile role for DNA in both modulating GR activity and also in directing the use of GR isoforms. We propose that the lever arm is a "fulcrum" for bidirectional allosteric signaling, conferring conformational changes in the DNA reading head that influence DNA sequence selectivity, as well as conferring changes in the dimerization domain that connect functionally with remote regulatory surfaces, thereby influencing which genes are regulated and the magnitude of their regulation.

RSAT peak-motifs: motif analysis in full-size ChIP-seq datasets.

ChIP-seq is increasingly used to characterize transcription factor binding and chromatin marks at a genomic scale. Various tools are now available to extract binding motifs from peak data sets. However, most approaches are only available as command-line programs, or via a website but with size restrictions. We present peak-motifs, a computational pipeline that discovers motifs in peak sequences, compares them with databases, exports putative binding sites for visualization in the UCSC genome browser and generates an extensive report suited for both naive and expert users. It relies on time- and memory-efficient algorithms enabling the treatment of several thousand peaks within minutes. Regarding time efficiency, peak-motifs outperforms all comparable tools by several orders of magnitude. We demonstrate its accuracy by analyzing data sets ranging from 4000 to 1,28,000 peaks for 12 embryonic stem cell-specific transcription factors. In all cases, the program finds the expected motifs and returns additional motifs potentially bound by cofactors. We further apply peak-motifs to discover tissue-specific motifs in peak collections for the p300 transcriptional co-activator. To our knowledge, peak-motifs is the only tool that performs a complete motif analysis and offers a user-friendly web interface without any restriction on sequence size or number of peaks.

Theoretical and empirical quality assessment of transcription factor-binding motifs.

Position-specific scoring matrices (PSSMs) are routinely used to predict transcription factor (TF)-binding sites in genome sequences. However, their reliability to predict novel binding sites can be far from optimum, due to the use of a small number of training sites or the inappropriate choice of parameters when building the matrix or when scanning sequences with it. Measures of matrix quality such as E-value and information content rely on theoretical models, and may fail in the context of full genome sequences. We propose a method, implemented in the program 'matrix-quality', that combines theoretical and empirical score distributions to assess reliability of PSSMs for predicting TF-binding sites. We applied 'matrix-quality' to estimate the predictive capacity of matrices for bacterial, yeast and mouse TFs. The evaluation of matrices from RegulonDB revealed some poorly predictive motifs, and allowed us to quantify the improvements obtained by applying multi-genome motif discovery. Interestingly, the method reveals differences between global and specific regulators. It also highlights the enrichment of binding sites in sequence sets obtained from high-throughput ChIP-chip (bacterial and yeast TFs), and ChIP-seq and experiments (mouse TFs). The method presented here has many applications, including: selecting reliable motifs before scanning sequences; improving motif collections in TFs databases; evaluating motifs discovered using high-throughput data sets.

RSAT 2011: regulatory sequence analysis tools.

RSAT (Regulatory Sequence Analysis Tools) comprises a wide collection of modular tools for the detection of cis-regulatory elements in genome sequences. Thirteen new programs have been added to the 30 described in the 2008 NAR Web Software Issue, including an automated sequence retrieval from EnsEMBL (retrieve-ensembl-seq), two novel motif discovery algorithms (oligo-diff and info-gibbs), a 100-times faster version of matrix-scan enabling the scanning of genome-scale sequence sets, and a series of facilities for random model generation and statistical evaluation (random-genome-fragments, random-motifs, random-sites, implant-sites, sequence-probability, permute-matrix). Our most recent work also focused on motif comparison (compare-matrices) and evaluation of motif quality (matrix-quality) by combining theoretical and empirical measures to assess the predictive capability of position-specific scoring matrices. To process large collections of peak sequences obtained from ChIP-seq or related technologies, RSAT provides a new program (peak-motifs) that combines several efficient motif discovery algorithms to predict transcription factor binding motifs, match them against motif databases and predict their binding sites. Availability (web site, stand-alone programs and SOAP/WSDL (Simple Object Access Protocol/Web Services Description Language) web services): http://rsat.ulb.ac.be/rsat/.

Logical modelling of in vitro differentiation of human monocytes into dendritic cells unravels novel transcriptional regulatory interactions.

Dendritic cells (DCs) are the major specialized antigen-presenting cells, thereby connecting innate and adaptive immunity. Because of their role in establishing adaptive immunity, they constitute promising targets for immunotherapy. Monocytes can differentiate into DCs in vitro in the presence of colony-stimulating factor 2 (CSF2) and interleukin 4 (IL4), activating four signalling pathways (MAPK, JAK/STAT, NFKB and PI3K). However, the downstream transcriptional programme responsible for DC differentiation from monocytes (moDCs) remains unknown. By analysing the scientific literature on moDC differentiation, we established a preliminary logical model that helped us identify missing information regarding the activation of genes responsible for this differentiation, including missing targets for key transcription factors (TFs). Using ChIP-seq and RNA-seq data from the Blueprint consortium, we defined active and inactive promoters, together with differentially expressed genes in monocytes, moDCs and macrophages, which correspond to an alternative cell fate. We then used this functional genomic information to predict novel targets for previously identified TFs. By integrating this information, we refined our model and recapitulated the main established facts regarding moDC differentiation. Prospectively, the resulting model should be useful to develop novel immunotherapies targeting moDCs.