Luc Jouneau

A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis

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.

IL-35-producing B cells are critical regulators of immunity during autoimmune and infectious diseases

B lymphocytes have critical roles as positive and negative regulators of immunity. Their inhibitory function has been associated primarily with interleukin 10 (IL-10) because B-cell-derived IL-10 can protect against autoimmune disease and increase susceptibility to pathogens. Here we identify IL-35-producing B cells as 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 markedly improved resistance to infection with the intracellular bacterial pathogen Salmonella enterica serovar Typhimurium as shown by their superior containment of the bacterial growth and their prolonged survival after primary infection, and upon secondary challenge, 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 increased function of B cells as antigen-presenting cells (APCs). 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 (also known as Prdm1). During EAE, CD138+ plasma cells were also the main source of B-cell-derived IL-35 and IL-10. Collectively, our data show the importance of IL-35-producing B cells in regulation of immunity and highlight IL-35 production by B cells as a potential therapeutic target for autoimmune and infectious diseases. This study reveals 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.

Maternal Diets Trigger Sex-Specific Divergent Trajectories of Gene Expression and Epigenetic Systems in Mouse Placenta

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.

Existence of CD8α-Like Dendritic Cells with a Conserved Functional Specialization and a Common Molecular Signature in Distant Mammalian Species

The mouse lymphoid organ-resident CD8α+ dendritic cell (DC) subset is specialized in Ag presentation to CD8+ T cells. Recent evidence shows that mouse nonlymphoid tissue CD103+ DCs and human blood DC Ag 3+ DCs share similarities with CD8α+ DCs. We address here whether the organization of DC subsets is conserved across mammals in terms of gene expression signatures, phenotypic characteristics, and functional specialization, independently of the tissue of origin. We study the DC subsets that migrate from the skin in the ovine species that, like all domestic animals, belongs to the Laurasiatheria, a distinct phylogenetic clade from the supraprimates (human/mouse). We demonstrate that the minor sheep CD26+ skin lymph DC subset shares significant transcriptomic similarities with mouse CD8α+ and human blood DC Ag 3+ DCs. This allowed the identification of a common set of phenotypic characteristics for CD8α-like DCs in the three mammalian species (i.e., SIRPlo, CADM1hi, CLEC9Ahi, CD205hi, XCR1hi). Compared to CD26− DCs, the sheep CD26+ DCs show 1) potent stimulation of allogeneic naive CD8+ T cells with high selective induction of the Ifnγ and Il22 genes; 2) dominant efficacy in activating specific CD8+ T cells against exogenous soluble Ag; and 3) selective expression of functional pathways associated with high capacity for Ag cross-presentation. Our results unravel a unifying definition of the CD8α+-like DCs across mammalian species and identify molecular candidates that could be used for the design of vaccines applying to mammals in general.

Teleost fish mount complex clonal IgM and IgT responses in spleen upon systemic viral infection

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.

Naive and primed murine pluripotent stem cells have distinct miRNA expression profiles.

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.

Origin and Evolution of TRIM Proteins: New Insights from the Complete TRIM Repertoire of Zebrafish and Pufferfish

Tripartite motif proteins (TRIM) constitute a large family of proteins containing a RING-Bbox-Coiled Coil motif followed by different C-terminal domains. Involved in ubiquitination, TRIM proteins participate in many cellular processes including antiviral immunity. The TRIM family is ancient and has been greatly diversified in vertebrates and especially in fish. We analyzed the complete sets of trim genes of the large zebrafish genome and of the compact pufferfish genome. Both contain three large multigene subsets - adding the hsl5/trim35-like genes (hltr) to the ftr and the btr that we previously described - all containing a B30.2 domain that evolved under positive selection. These subsets are conserved among teleosts. By contrast, most human trim genes of the other classes have only one or two orthologues in fish. Loss or gain of C-terminal exons generated proteins with different domain organizations; either by the deletion of the ancestral domain or, remarkably, by the acquisition of a new C-terminal domain. Our survey of fish trim genes in fish identifies subsets with different evolutionary dynamics. trims encoding RBCC-B30.2 proteins show the same evolutionary trends in fish and tetrapods: they evolve fast, often under positive selection, and they duplicate to create multigenic families. We could identify new combinations of domains, which epitomize how new trim classes appear by domain insertion or exon shuffling. Notably, we found that a cyclophilin-A domain replaces the B30.2 domain of a zebrafish fintrim gene, as reported in the macaque and owl monkey antiretroviral TRIM5α. Finally, trim genes encoding RBCC-B30.2 proteins are preferentially located in the vicinity of MHC or MHC gene paralogues, which suggests that such trim genes may have been part of the ancestral MHC.

Enhanced or Reduced Fetal Growth Induced by Embryo Transfer into Smaller or Larger Breeds Alters Post-Natal Growth and Metabolism in Pre-Weaning Horses

In equids, placentation is diffuse and nutrient supply to the fetus is determined by uterine size. This correlates with maternal size and affects intra-uterine development and subsequent post-natal growth, as well as insulin sensitivity in the newborn. Long-term effects remain to be described. In this study, fetal growth was enhanced or restricted through ET using pony (P), saddlebred (S) and draft (D) horses. Control P-P (n = 21) and S-S (n = 28) pregnancies were obtained by AI. Enhanced and restricted pregnancies were obtained by transferring P or S embryos into D mares (P-D, n = 6 and S-D, n = 8) or S embryos into P mares (S-P, n = 6), respectively. Control and experimental foals were raised by their dams and recipient mothers, respectively. Weight gain, growth hormones and glucose homeostasis were investigated in the foals from birth to weaning. Fetal growth was enhanced in P-D and these foals remained consistently heavier, with reduced T3 concentrations until weaning compared to P-P. P-D had lower fasting glucose from days 30 to 200 and higher insulin secretion than P-P after IVGTT on day 3. Euglycemic clamps in the immediate post-weaning period revealed no difference in insulin sensitivity between P-D and P-P. Fetal growth was restricted in S-P and these foals remained consistently lighter until weaning compared to S-D, with elevated T3 concentrations in the newborn compared to S-S. S-P exhibited higher fasting glycemia than S-S and S-D from days 30 to 200. They had higher maximum increment in plasma glucose than S-D after IVGTT on day 3 and clamps on day 200 demonstrated higher insulin sensitivity compared to S-D. Neither the restricted nor the enhanced fetal environment affected IGF-1 concentrations. Thus, enhanced and restricted fetal and post-natal environments had combined effects that persisted until weaning. They induced different adaptive responses in post-natal glucose metabolism: an early insulin-resistance was induced in enhanced P-D, while S-P developed increased insulin sensitivity.

Induced pluripotent stem cells derived from rabbits exhibit some characteristics of naïve pluripotency

Summary Not much is known about the molecular and functional features of pluripotent stem cells (PSCs) in rabbits. To address this, we derived and characterized 2 types of rabbit PSCs from the same breed of New Zealand White rabbits: 4 lines of embryonic stem cells (rbESCs), and 3 lines of induced PSCs (rbiPSCs) that were obtained by reprogramming adult skin fibroblasts. All cell lines required fibroblast growth factor 2 for their growth and proliferation. All rbESC lines showed molecular and functional properties typically associated with primed pluripotency. The cell cycle of rbESCs had a prolonged G1 phase and a DNA damage checkpoint before entry into the S phase, which are the 2 features typically associated with the somatic cell cycle. In contrast, the rbiPSC lines exhibited some characteristics of naïve pluripotency, including resistance to single-cell dissociation by trypsin, robust activity of the distal enhancer of the mouse Oct4 gene, and expression of naïve pluripotency-specific genes, as defined in rodents. According to gene expression profiles, rbiPSCs were closer to the rabbit inner cell mass (ICM) than rbESCs. Furthermore, rbiPSCs were capable of colonizing the ICM after aggregation with morulas. Therefore, we propose that rbiPSCs self-renew in an intermediate state between naïve and primed pluripotency, which represents a key step toward the generation of bona fide naïve PSC lines in rabbits.

Culicoides Midge Bites Modulate the Host Response and Impact on Bluetongue Virus Infection in Sheep

Many haematophagous insects produce factors that help their blood meal and coincidently favor pathogen transmission. However nothing is known about the ability of Culicoides midges to interfere with the infectivity of the viruses they transmit. Among these, Bluetongue Virus (BTV) induces a hemorrhagic fever- type disease and its recent emergence in Europe had a major economical impact. We observed that needle inoculation of BTV8 in the site of uninfected C. nubeculosus feeding reduced viraemia and clinical disease intensity compared to plain needle inoculation. The sheep that developed the highest local inflammatory reaction had the lowest viral load, suggesting that the inflammatory response to midge bites may participate in the individual sensitivity to BTV viraemia development. Conversely compared to needle inoculation, inoculation of BTV8 by infected C. nubeculosus bites promoted viraemia and clinical symptom expression, in association with delayed IFN- induced gene expression and retarded neutralizing antibody responses. The effects of uninfected and infected midge bites on BTV viraemia and on the host response indicate that BTV transmission by infected midges is the most reliable experimental method to study the physio-pathological events relevant to a natural infection and to pertinent vaccine evaluation in the target species. It also leads the way to identify the promoting viral infectivity factors of infected Culicoides in order to possibly develop new control strategies against BTV and other Culicoides transmitted viruses.