During the last 3 years, a number of approaches for the normalization of RNA sequencing data have emerged in the literature, differing both in the type of bias adjustment and in the statistical strategy adopted. However, as data continue to accumulate, there has been no clear consensus on the appropriate normalization method to be used or the impact of a chosen method on the downstream analysis. In this work, we focus on a comprehensive comparison of seven recently proposed normalization methods for the differential analysis of RNA-seq data, with an emphasis on the use of varied real and simulated datasets involving different species and experimental designs to represent data characteristics commonly observed in practice. Based on this comparison study, we propose practical recommendations on the appropriate normalization method to be used and its impact on the differential analysis of RNA-seq data.

SUMMARY B lymphocytes have critical roles as positive and negative regulators of immunity. Their inhibitory function has so far been associated primarily with interleukin (IL)-10 because B cell-derived IL-10 can protect against autoimmune disease and increase susceptibility to pathogens1,2. Here, we identify IL-35-producing B cells as novel key players in the negative regulation of immunity. Mice in which only B cells did not express IL-35 lost their ability to recover from the T cell-mediated demyelinating autoimmune disease experimental autoimmune encephalomyelitis (EAE). In contrast, these mice displayed a strikingly improved resistance to infection with the intracellular bacterial pathogen Salmonella typhimurium, as shown by their superior containment of the bacterial growth and their prolonged survival both after primary infection, and upon secondary challenge after vaccination, compared to control mice. The increased immunity found in mice lacking IL-35 production by B cells was associated with a higher activation of macrophages and inflammatory T cells, as well as an enhanced stimulatory function of B cells as antigen-presenting cells (APC). During Salmonella infection IL-35- and IL-10-producing B cells corresponded to two largely distinct sets of surface-IgM+CD138hiTACI+CXCR4+CD1dintTim1int plasma cells expressing the transcription factor Blimp1. During EAE CD138+ plasma cells were also the major source of B cell-derived IL-35 and IL-10. Collectively, our data unravel the importance of IL-35-producing B cells in regulation of immunity, and highlight IL-35 production by B cells as a novel therapeutic target for autoimmune and infectious diseases. More generally, this study emphasizes the central role of activated B cells, particularly plasma cells, and their production of cytokines in the regulation of immune responses in health and disease.

SummaryB lymphocytes can suppress immunity through interleukin (IL)-10 production in infectious, autoimmune, and malignant diseases. Here, we have identified a natural plasma cell subset that distinctively expresses the inhibitory receptor LAG-3 and mediates this function in vivo. These plasma cells also express the inhibitory receptors CD200, PD-L1, and PD-L2. They develop from various B cell subsets in a B cell receptor (BCR)-dependent manner independently of microbiota in naive mice. After challenge they upregulate IL-10 expression via a Toll-like receptor-driven mechanism within hours and without proliferating. This function is associated with a unique transcriptome and epigenome, including the lowest amount of DNA methylation at the Il10 locus compared to other B cell subsets. Their augmented accumulation in naive mutant mice with increased BCR signaling correlates with the inhibition of memory T cell formation and vaccine efficacy after challenge. These natural regulatory plasma cells may be of broad relevance for disease intervention.

Males and females responses to gestational overnutrition set the stage for subsequent sex-specific differences in adult onset non communicable diseases. Placenta, as a widely recognized programming agent, contibutes to the underlying processes. According to our previous findings, a high-fat diet during gestation triggers sex-specific epigenetic alterations within CpG and throughout the genome, together with the deregulation of clusters of imprinted genes. We further investigated the impact of diet and sex on placental histology, transcriptomic and epigenetic signatures in mice. Both basal gene expression and response to maternal high-fat diet were sexually dimorphic in whole placentas. Numerous genes showed sexually dimorphic expression, but only 11 genes regardless of the diet. In line with the key role of genes belonging to the sex chromosomes, 3 of these genes were Y-specific and 3 were X-specific. Amongst all the genes that were differentially expressed under a high-fat diet, only 16 genes were consistently affected in both males and females. The differences were not only quantitative but remarkably qualitative. The biological functions and networks of genes dysregulated differed markedly between the sexes. Seven genes of the epigenetic machinery were dysregulated, due to effects of diet, sex or both, including the Y- and X-linked histone demethylase paralogues Kdm5c and Kdm5d, which could mark differently male and female epigenomes. The DNA methyltransferase cofactor Dnmt3l gene expression was affected, reminiscent of our previous observation of changes in global DNA methylation. Overall, this striking sexual dimorphism of programming trajectories impose a considerable revision of the current dietary interventions protocols.

Material

Upon infection, B-lymphocytes expressing antibodies specific for the intruding pathogen develop clonal responses triggered by pathogen recognition via the B-cell receptor. The constant region of antibodies produced by such responding clones dictates their functional properties. In teleost fish, the clonal structure of B-cell responses and the respective contribution of the three isotypes IgM, IgD and IgT remain unknown. The expression of IgM and IgT are mutually exclusive, leading to the existence of two B-cell subsets expressing either both IgM and IgD or only IgT. Here, we undertook a comprehensive analysis of the variable heavy chain (VH) domain repertoires of the IgM, IgD and IgT in spleen of homozygous isogenic rainbow trout (Onchorhynchus mykiss) before, and after challenge with a rhabdovirus, the Viral Hemorrhagic Septicemia Virus (VHSV), using CDR3-length spectratyping and pyrosequencing of immunoglobulin (Ig) transcripts. In healthy fish, we observed distinct repertoires for IgM, IgD and IgT, respectively, with a few amplified μ and τ junctions, suggesting the presence of IgM- and IgT-secreting cells in the spleen. In infected animals, we detected complex and highly diverse IgM responses involving all VH subgroups, and dominated by a few large public and private clones. A lower number of robust clonal responses involving only a few VH were detected for the mucosal IgT, indicating that both IgM+ and IgT+ spleen B cells responded to systemic infection but at different degrees. In contrast, the IgD response to the infection was faint. Although fish IgD and IgT present different structural features and evolutionary origin compared to mammalian IgD and IgA, respectively, their implication in the B-cell response evokes these mouse and human counterparts. Thus, it appears that the general properties of antibody responses were already in place in common ancestors of fish and mammals, and were globally conserved during evolution with possible functional convergences.

Over the last years, the microRNA (miRNA) pathway has emerged as a key component of the regulatory network of pluripotency. Although clearly distinct states of pluripotency have been described in vivo and ex vivo, differences in miRNA expression profiles associated with the developmental modulation of pluripotency have not been extensively studied so far. Here, we performed deep sequencing to profile miRNA expression in naive (embryonic stem cell [ESC]) and primed (epiblast stem cell [EpiSC]) pluripotent stem cells derived from mouse embryos of identical genetic background. We developed a graphical representation method allowing the rapid identification of miRNAs with an atypical profile including mirtrons, a small nucleolar RNA (snoRNA)-derived miRNA, and miRNAs whose biogenesis may differ between ESC and EpiSC. Comparison of mature miRNA profiles revealed that ESCs and EpiSCs exhibit very different miRNA signatures with one third of miRNAs being differentially expressed between the two cell types. Notably, differential expression of several clusters, including miR290-295, miR17-92, miR302/367, and a large repetitive cluster on chromosome 2, was observed. Our analysis also showed that differentiation priming of EpiSC compared to ESC is evidenced by changes in miRNA expression. These dynamic changes in miRNAs signature are likely to reflect both redundant and specific roles of miRNAs in the fine-tuning of pluripotency during development.

Whereas major histocompatibility class-1 (MH1) proteins present peptides to T cells displaying a large T-cell receptor (TR) repertoire, MH1Like proteins, such as CD1D and MR1, present glycolipids and microbial riboflavin precursor derivatives, respectively, to T cells expressing invariant TR-α (iTRA) chains. The groove of such MH1Like, as well as iTRA chains used by mucosal-associated invariant T (MAIT) and natural killer T (NKT) cells, respectively, may result from a coevolution under particular selection pressures. Herein, we investigated the evolutionary patterns of the iTRA of MAIT and NKT cells and restricting MH1Like proteins: MR1 appeared 170 Mya and is highly conserved across mammals, evolving more slowly than other MH1Like. It has been pseudogenized or independently lost three times in carnivores, the armadillo, and lagomorphs. The corresponding TRAV1 gene also evolved slowly and harbors highly conserved complementarity determining regions 1 and 2. TRAV1 is absent exclusively from species in which MR1 is lacking, suggesting that its loss released the purifying selection on MR1. In the rabbit, which has very few NKT and no MAIT cells, a previously unrecognized iTRA was identified by sequencing leukocyte RNA. This iTRA uses TRAV41, which is highly conserved across several groups of mammals. A rabbit MH1Like gene was found that appeared with mammals and is highly conserved. It was independently lost in a few groups in which MR1 is present, like primates and Muridae, illustrating compensatory emergences of new MH1Like/Invariant T-cell combinations during evolution. Deciphering their role is warranted to search similar effector functions in humans.MHC | TCR | MAIT | evolution | mammals