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

T cell activation is initiated upon ligand engagement of the T cell receptor (TCR) and costimulatory receptors. The CD28 molecule acts as a major costimulatory receptor in promoting full activation of naive T cells. However, despite extensive studies, why naive T cell activation requires concurrent stimulation of both the TCR and costimulatory receptors remains poorly understood. Here, we explore this issue by analyzing calcium response as a key early signaling event to elicit T cell activation. Experiments using mouse naive CD4+ T cells showed that engagement of the TCR or CD28 with the respective cognate ligand was able to trigger a rise in fluctuating calcium mobilization levels, as shown by the frequency and average response magnitude of the reacting cells compared with basal levels occurred in unstimulated cells. The engagement of both TCR and CD28 enabled a further increase of these two metrics. However, such increases did not sufficiently explain the importance of the CD28 pathways to the functionally relevant calcium responses in T cell activation. Through the autocorrelation analysis of calcium time series data, we found that combined but not separate TCR and CD28 stimulation significantly prolonged the average decay time (τ) of the calcium signal amplitudes determined with the autocorrelation function, compared with its value in unstimulated cells. This increasement of decay time (τ) uniquely characterizes the fluctuating calcium response triggered by concurrent stimulation of TCR and CD28, as it could not be achieved with either stronger TCR stimuli or by co-engaging both TCR and LFA-1, and likely represents an important feature of competent early signaling to provoke efficient T cell activation. Our work has thus provided new insights into the interplay between the TCR and CD28 early signaling pathways critical to trigger naive T cell activation.

Natural killer (NK) cells are innate lymphoid cells (ILCs) involved in antimicrobial and antitumoral responses. Several NK cell subsets have been reported in humans and mice, but their heterogeneity across organs and species remains poorly characterized. We assessed the diversity of human and mouse NK cells by single-cell RNA sequencing on thousands of individual cells isolated from spleen and blood. Unbiased transcriptional clustering revealed two distinct signatures differentiating between splenic and blood NK cells. This analysis at single-cell resolution identified three subpopulations in mouse spleen and four in human spleen, and two subsets each in mouse and human blood. A comparison of transcriptomic profiles within and between species highlighted the similarity of the two major subsets, NK1 and NK2, across organs and species. This unbiased approach provides insight into the biology of NK cells and establishes a rationale for the translation of mouse studies to human physiology and disease.

In the present paper, an experimental feasibility study on the detection of long-range intermolecular interactions through three-dimensional molecular diffusion in solution is performed. This follows recent theoretical and numerical analyses reporting that long-range electrodynamic forces between biomolecules could be identified through deviations from Brownian diffusion. The suggested experimental technique was fluorescence correlation spectroscopy (FCS). By considering two oppositely charged molecular species in aqueous solution, namely, lysozymes and fluorescent dye molecules (Alexa488), the diffusion coefficient of the dyes has been measured for different values of the concentration of lysozyme, that is, for different average distances between the oppositely charged molecules. For our model, long-range interactions are of electrostatic origin, suggesting that their action radius can be varied by changing the ionic strength of the solution. The experimental outcomes clearly prove the detectability of long-range intermolecular interactions by means of the FCS technique. Molecular dynamics simulations provide a clear and unambiguous interpretation of the experimental results.

In innate immune responses, induction of type‐I interferons (IFNs) prevents virus spreading while viral replication is delayed by protein synthesis inhibition. We asked how cells perform these apparently contradictory activities. Using single fibroblast monitoring by flow cytometry and mathematical modeling, we demonstrate that type‐I IFN production is linked to cell's ability to enter dsRNA‐activated PKR‐dependent translational arrest and then overcome this inhibition by decreasing eIF2α phosphorylation through phosphatase 1c cofactor GADD34 (Ppp1r15a) expression. GADD34 expression, shown here to be dependent on the IRF3 transcription factor, is responsible for a biochemical cycle permitting pulse of IFN synthesis to occur in cells undergoing protein synthesis inhibition. Translation arrest is further demonstrated to be key for anti‐viral response by acting synergistically with MAVS activation to amplify TBK1 signaling and IFN‐β mRNA transcription, while GADD34‐dependent protein synthesis recovery contributes to the heterogeneous expression of IFN observed in dsRNA‐activated cells.

One paradigm of random monoallelic gene expression is that of T-cell receptor (TCR) β allelic exclusion in T lymphocytes. However, the dynamics that sustain asymmetric choice in TCR β dual allele usage and the production of TCR β monoallelic expressing T-cells remain poorly understood. Here, we develop a computational model to explore a scheme of TCR β allelic exclusion based on the stochastic initiation of DNA rearrangement [V(D)J recombination] at homologous alleles in T-cell progenitors, and thus account for the genotypic profiles typically associated with allelic exclusion in differentiated T-cells. Disturbances in these dynamics at the level of an individual allele have limited consequences on these pro1les, robust feature of the system that is underscored by our simulations. Our study predicts a biological system in which locus-specific, prime epigenetic allelic activation effects set the stage to both optimize the production of TCR β allelically excluded T-cells and curtail the emergence of their allelically included counterparts.

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.

As you can see below, the referees raised a number of constructive concerns. One important point is that further in vivo and virus infection data is needed to support the physiological relevance of the key findings. Both referees #1 and 2 raise this point. This issue and the other raised points need to be resolved for consideration here. To extend the analysis along the lines indicated by the referees will require some efforts. Therefore before taking the decision in this case I would like to ask you to provide me with a point-by-point response of what you can do within 3-6 months to address the raised concerns. Based upon your response I will then take the decision on your manuscript.

T-cell acute lymphoblastic leukaemias (T-ALL) are aggressive malignant proliferations characterized by high relapse rates and great genetic heterogeneity. TAL1 is amongst the most frequently deregulated oncogenes. Yet, over half of the TAL1(+) cases lack TAL1 lesions, suggesting unrecognized (epi)genetic deregulation mechanisms. Here we show that TAL1 is normally silenced in the T-cell lineage, and that the polycomb H3K27me3-repressive mark is focally diminished in TAL1(+) T-ALLs. Sequencing reveals that >20% of monoallelic TAL1(+) patients without previously known alterations display microinsertions or RAG1/2-mediated episomal reintegration in a single site 5' to TAL1. Using 'allelic-ChIP' and CrispR assays, we demonstrate that such insertions induce a selective switch from H3K27me3 to H3K27ac at the inserted but not the germline allele. We also show that, despite a considerable mechanistic diversity, the mode of oncogenic TAL1 activation, rather than expression levels, impact on clinical outcome. Altogether, these studies establish site-specific epigenetic desilencing as a mechanism of oncogenic activation.

Allelic exclusion of immunoglobulin (Ig) and T-cell receptor (TCR) genes ensures the development of B and T lymphocytes operating under the mode of clonal selection. This phenomenon associates asynchronous V(D)J recombination events at Ig or TCR alleles and inhibitory feedback control. Despite years of intense research, however, the mechanisms that sustain asymmetric choice in random Ig/TCR dual allele usage and the production of Ig/TCR monoallelic expressing B and T lymphocytes remain unclear and open for debate. In this chapter, we first recapitulate the biological evidence that almost from the start appeared to link V(D)J recombination and allelic exclusion. We review the theoretical models previously proposed to explain this connection. Finally, we introduce our own mathematical modeling views based on how the developmental dynamics of individual lymphoid cells combine to sustain allelic exclusion.

The dynamical properties and diffusive behavior of a collection of mutually interacting particles are numerically investigated for two types of long-range interparticle interactions: Coulomb-electrostatic and dipole-electrodynamic. It is shown that when the particles are uniformly distributed throughout the accessible space, the self-diffusion coefficient is always lowered by the considered interparticle interactions, irrespective of their attractive or repulsive character. This fact is also confirmed by a simple model to compute the correction to the Brownian diffusion coefficient due to the interactions among the particles. These interactions are also responsible for the onset of dynamical chaos and an associated chaotic diffusion which still follows an Einstein-Fick-like law for the mean-square displacement as a function of time. Transitional phenomena are observed for Coulomb-electrostatic (repulsive) and dipole-electrodynamic (attractive) interactions considered both separately and in competition. The outcomes reported in this paper clearly indicate a feasible experimental method to probe the activation of resonant electrodynamic interactions among biomolecules.

To perform their specific functional role, B and T lymphocytes, cells of the adaptive immune system of jawed vertebrates, need to express one (and, preferably, only one) form of antigen receptor, i.e. the immunoglobulin or T-cell receptor (TCR), respectively. This end goal depends initially on a series of DNA cis-rearrangement events between randomly chosen units from separate clusters of V, D (at some immunoglobulin and TCR loci) and J gene segments, a biomolecular process collectively referred to as V(D)J recombination. V(D)J recombination takes place in immature T and B cells and relies on the so-called RAG nuclease, a site-specific DNA cleavage apparatus that corresponds to the lymphoid-specific moiety of the VDJ recombinase. At the genome level, this recombinase's mission presents substantial biochemical challenges. These relate to the huge distance between (some of) the gene segments that it eventually rearranges and the need to achieve cell-lineage-restricted and developmentally ordered routines with at times, mono-allelic versus bi-allelic discrimination. The entire process must be completed without any recombination errors, instigators of chromosome instability, translocation and, potentially, tumorigenesis. As expected, such a precisely choreographed and yet potentially risky process demands sophisticated controls; epigenetics demonstrates what is possible when calling upon its many facets. In this vignette, we will recall the evidence that almost from the start appeared to link the two topics, V(D)J recombination and epigenetics, before reviewing the latest advances in our knowledge of this joint venture.

Immunology - though deeply experimental in everyday practice - is also a theoretical discipline. Recent advances in the understanding of innate immunity, how it is triggered and how it shares features that have previously been uniquely ascribed to the adaptive immune system, can contribute to the refinement of the theoretical framework of immunology. In particular, natural killer cells and macrophages are activated by transient modifications, but adapt to long-lasting modifications that occur in the surrounding tissue environment. This process facilitates the maintenance of self-tolerance while permitting efficient immune responses. In this Essay we extend this idea to other components of the immune system and we propose some general principles that lay the foundations for a unifying theory of immunity - the discontinuity theory. According to this theoretical framework, effector immune responses (namely, activated responses that lead to the potential elimination of the target antigen) are induced by an antigenic discontinuity; that is, by the sudden modification of molecular motifs with which immune cells interact.

BACKGROUND:DNA methylation contributes to the regulation of gene expression during development and cellular differentiation. The recently developed Methylated DNA ImmunoPrecipitation (MeDIP) assay allows a comprehensive analysis of this epigenetic mark at the genomic level in normal and disease-derived cells. However, estimating the efficiency of the MeDIP technique is difficult without previous knowledge of the methylation status of a given cell population. Attempts to circumvent this problem have involved the use of in vitro methylated DNA in parallel to the investigated samples. Taking advantage of this stratagem, we sought to improve the sensitivity of the approach and to assess potential biases resulting from DNA amplification and hybridization procedures using MeDIP samples.FINDINGS:We performed MeDIP assays using in vitro methylated DNA, with or without previous DNA amplification, and hybridization to a human promoter array. We observed that CpG content at gene promoters indeed correlates strongly with the MeDIP signal obtained using in vitro methylated DNA, even when lowering significantly the amount of starting material. In analyzing MeDIP products that were subjected to whole genome amplification (WGA), we also revealed a strong bias against CpG-rich promoters during this amplification procedure, which may potentially affect the significance of the resulting data.CONCLUSION:We illustrate the use of in vitro methylated DNA to assess the efficiency and accuracy of MeDIP procedures. We report that efficient and reproducible genome-wide data can be obtained via MeDIP experiments using relatively low amount of starting genomic DNA; and emphasize for the precaution that must be taken in data analysis when an additional DNA amplification step is required.

Allelic exclusion represents a major aspect of TCRβ gene assembly by V(D)J recombination in developing T lymphocytes. Despite recent progress, its comprehension remains problematic when confronted with experimental data. Existing models fall short in terms of incorporating into a unique distribution all the cell subsets emerging from the TCRβ assembly process. To revise this issue, we propose dynamical, continuous-time Markov chain-based modeling whereby essential steps in the biological procedure (D-J and V-DJ rearrangements and feedback inhibition) evolve independently on the two TCRβ alleles in every single cell while displaying random modes of initiation and duration. By selecting parameters via fitting procedures, we demonstrate the capacity of the model to offer accurate fractions of all distinct TCRβ genotypes observed in studies using developing and mature T cells from wild-type or mutant mice. Selected parameters in turn afford relative duration for each given step, hence updating TCRβ recombination distinctive timings. Overall, our dynamical modeling integrating allele independence and noise in recombination and feedback-inhibition events illustrates how the combination of these ingredients alone may enforce allelic exclusion at the TCRβ locus.

Motivation: High-density tiling microarrays are increasingly used in combination with ChIP assays to study transcriptional regulation. To ease the analysis of the large amounts of data generated by this approach, we have developed ChIP-on-chip Analysis Suite (CoCAS), a standalone software suite which implements optimized ChIP-on-chip data normalization, improved peak detection, as well as quality control reports. Our software allows dye swap, replicate correlation and connects easily with genome browsers and other peak detection algorithms. CoCAS can readily be used on the latest generation of Agilent high-density arrays. Also, the implemented peak detection methods are suitable for other datasets, including ChIP-Seq output. Availability: The software is available for download along with a sample dataset at http://www.ciml.univ-mrs.fr/software/ferrier.htm. Contact: ferrier@ciml.univ-mrs.fr; andrau@ciml.univ-mrs.fr; spicuglia@ciml.univ-mrs.fr Supplementary information: Supplementary data are available at Bioinformatics online.

The TCRβ gene enhancer (Eβ) commands TCRβ gene expression through the lifespan of T lymphocytes. Genetic and molecular studies have implied that in early thymocytes, Eβ directs chromatin opening over the Dβ-Jβ-Cβ domains and triggers initial Dβ-Jβ recombination. In mature T cells, Eβ is required for expression of the assembled TCRβ gene. Whether these separate activities rely on distinct Eβ regulatory sequences and involve differing modes of activation is unclear. Using gene targeting in mouse embryonic stem cells, we replaced Eβ by a conserved core fragment (Eβ169). We found that Eβ169-carrying alleles were capable of sustaining β gene expression and the development of mature T cells in homozygous knockin mice. Surprisingly, these procedures and underlying molecular transactions were affected to a wide range of degrees depending on the developmental stage. Early thymocytes barely achieved Dβ-Jβ germline transcription and recombination. In contrast, T cells displayed substantial though heterogeneous levels of VDJ-rearranged TCRβ gene expression. Our results have implications regarding enhancer function in cells of the adaptive immune system and, potentially, TCRβ gene recombination and allelic exclusion.

T cell receptor (TCR) β gene assembly by V(D)J recombination proceeds via successive Dβ-to-Jβ and Vβ-to-DJβ rearrangements. This two-step process is enforced by a constraint, termed beyond (B)12/23, which prohibits direct Vβ-to-Jβ rearrangements. However the B12/23 restriction does not explain the order of TCRβ assembly for which the regulation remains an unresolved issue. The initiation of V(D)J recombination consists of the introduction of single-strand DNA nicks at recombination signal sequences (RSSs) containing a 12 base-pairs spacer. An RSS containing a 23 base-pairs spacer is then captured to form a 12/23 RSSs synapse leading to coupled DNA cleavage. Herein, we probed RSS nicks at the TCRβ locus and found that nicks were only detectable at Dβ-associated RSSs. This pattern implies that Dβ 12RSS and, unexpectedly, Dβ 23RSS initiate V(D)J recombination and capture their respective Vβ or Jβ RSS partner. Using both in vitro and in vivo assays, we further demonstrate that the Dβ1 23RSS impedes cleavage at the adjacent Dβ1 12RSS and consequently Vβ-to-Dβ1 rearrangement first requires the Dβ1 23RSS excision. Altogether, our results provide the molecular explanation to the B12/23 constraint and also uncover a ‘Dβ1 23RSS-mediated’ restriction operating beyond chromatin accessibility, which directs Dβ1 ordered rearrangements.

Motivation: High-density tiling microarrays are increasingly used in combination with ChIP assays to study transcriptional regulation. To ease the analysis of the large amounts of data generated by this approach, we have developed ChIP-on-chip Analysis Suite (CoCAS), a standalone software suite which implements optimized ChIP-onchip data normalization, improved peak detection, as well as quality control reports. Our software allows dye swap, replicate correlation and connects easily with genome browsers and other peak detection algorithms. CoCAS can readily be used on the latest generation of Agilent high-density arrays. Also, the implemented peak detection methods are suitable for other datasets, including ChIP-Seq output. Availability: The software is available for download along with a sample dataset at http://www.ciml.univ-mrs.fr/software/ferrier.htm. Contact: ferrier@ciml.univ-mrs.fr; andrau@ciml.univ-mrs.fr; spicug lia@ciml.univ-mrs.fr Supplementary information: Supplementary data are available at Bioinformatics online.

Since an immunoreceptor tyrosine-based inhibition motif (ITIM) was first identified in the intracytoplasmic domain of Fc gamma RIIB, ITIMs have been found in a large number of inhibitory molecules that were shown to negatively regulate cell activation. Due to their wide tissue distribution and to the variety of their extracellular ligands, ITIM-containing molecules are involved in the control of a large spectrum of biological functions, mostly but not exclusively related to immunity. On the basis of sequence comparison, ITIMs were structurally defined as 6-amino acid sequences containing a tyrosine (Y) with loosely conserved N-terminal (Y-2) and C-terminal (Y+3) residues. Molecular analysis of signaling events demonstrated that when coaggregated with activating receptors, ITIMs are phosphorylated by Src-family tyrosine kinases, which enables them to recruit Src homology 2 domain-containing phosphatases that antagonize activation signals. Because ITIM-dependent negative regulation seems to be a fundamental regulatory mechanism, both in rodents and in humans, and because it can be used either as a target or as a powerful tool in various diseases, we undertook (i) a genome-wide search of potential novel ITIM-containing molecules in humans, mice, frogs, birds, and flies and (ii) a comparative analysis of potential ITIMs in major animal phyla, from mammals to protozoa. We found a surprisingly high number of potential ITIM-containing molecules, having a great diversity of extracellular domains, and being expressed by a variety of immune and non-immune cells. ITIMs could be traced back to the most primitive metazoa. The genes that encode ITIM-containing molecules that belong to the immunoglobulin superfamily or to the C-lectin family seem to derive from a common set of ancestor genes and to have dramatically expanded and diverged in Gnathostomata (from fish to mammals).