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.

BackgroundSeveral R packages exist for the detection of differentially expressed genes from RNA-Seq data. The analysis process includes three main steps, namely normalization, dispersion estimation and test for differential expression. Quality control steps along this process are recommended but not mandatory, and failing to check the characteristics of the dataset may lead to spurious results. In addition, normalization methods and statistical models are not exchangeable across the packages without adequate transformations the users are often not aware of. Thus, dedicated analysis pipelines are needed to include systematic quality control steps and prevent errors from misusing the proposed methods.ResultsSARTools is an R pipeline for differential analysis of RNA-Seq count data. It can handle designs involving two or more conditions of a single biological factor with or without a blocking factor (such as a batch effect or a sample pairing). It is based on DESeq2 and edgeR and is composed of an R package and two R script templates (for DESeq2 and edgeR respectively). Tuning a small number of parameters and executing one of the R scripts, users have access to the full results of the analysis, including lists of differentially expressed genes and a HTML report that (i) displays diagnostic plots for quality control and model hypotheses checking and (ii) keeps track of the whole analysis process, parameter values and versions of the R packages used.ConclusionsSARTools provides systematic quality controls of the dataset as well as diagnostic plots that help to tune the model parameters. It gives access to the main parameters of DESeq2 and edgeR and prevents untrained users from misusing some functionalities of both packages. By keeping track of all the parameters of the analysis process it fits the requirements of reproducible research.

SummaryAdaptation to the host environment and exploitation of host cell functions are critical to the success of intracellular pathogens. Here, insight to these virulence mechanisms was obtained for the first time from the transcriptional program of the human pathogen Legionella pneumophila during infection of its natural host, Acanthamoeba castellanii . The biphasic life cycle of L. pneumophila was reflected by a major shift in gene expression from replicative to transmissive phase, concerning nearly half of the genes predicted in the genome. However, three different L. pneumophila strains showed similar in vivo gene expression patterns, indicating that common regulatory mechanisms govern the Legionella life cycle, despite the plasticity of its genome. During the replicative phase, in addition to components of aerobic metabolism and amino acid catabolism, the Entner-Doudoroff pathway, a NADPH producing mechanism used for sugar and/or gluconate assimilation, was expressed, suggesting for the first time that intracellular L. pneumophila may also scavenge host carbohydrates as nutrients and not only proteins. Identification of genes only upregulated in vivo but not in vitro , may explain higher virulence of in vivo grown L. pneumophila. Late in the life cycle, L. pneumophila upregulates genes predicted to promote transmission and manipulation of a new host cell, therewith priming it for the next attack. These including substrates of the Dot/Icm secretion system, other factors associated previously with invasion and virulence, the motility and the type IV pilus machineries, and > 90 proteins not characterized so far. Analysis of a fliA ( s 28 ) deletion mutant identified genes coregulated with the flagellar regulon, including GGDEF/EAL regulators and factors that promote host cell entry and survival.

Listeria monocytogenes is responsible for gastroenteritis in healthy individuals and for a severe invasive disease in immunocompromised patients. Among the three identified L. monocytogenes evolutionary lineages, lineage I strains are overrepresented in epidemic listeriosis outbreaks, but the mechanisms underlying the higher virulence potential of strains of this lineage remain elusive. Here, we demonstrate that Listeriolysin S (LLS), a virulence factor only present in a subset of lineage I strains, is a bacteriocin highly expressed in the intestine of orally infected mice that alters the host intestinal microbiota and promotes intestinal colonization by L. monocytogenes, as well as deeper organ infection. To our knowledge, these results therefore identify LLS as the first bacteriocin described in L. monocytogenes and associate modulation of host microbiota by L. monocytogenes epidemic strains to increased virulence.T he gram-positive bacterium Listeria monocytogenes is a facultative intracellular pathogen that causes foodborne infections in humans and animals. Upon consumption of contaminated food, L. monocytogenes reaches the intestinal lumen, crosses the intestinal barrier, and disseminates within the host. The clinical manifestations of listeriosis vary from a mild, self-limiting gastroenteritis to severe intestinal and systemic infections, with a fatality rate estimated to 20-30% of infected individuals (1). Host gut microbiota plays a critical role in resistance against colonization by invading pathogens within the intestine (2). The mechanisms of L. monocytogenes to compete with the host microbiota to survive in the intestine remain unknown.During the last decades, the majority of Listeria studies in bacterial pathophysiology, cell biology, and immunology compared three pathogenic strains from lineage II: EGD, EGD-e, and 10403S (3). Interestingly, major listeriosis epidemics have been preferentially associated to L. monocytogenes clonal groups belonging to the evolutionary lineage I and, more specifically, to serotype 4b (4, 5), but the molecular mechanisms that contribute to the higher virulence potential of these bacterial strains have not been identified yet.Bacteriocins are bacterially synthesized proteinaceous substances that target and inhibit the growth of closely related bacteria, allowing competition in diverse ecological niches, including the digestive tract (6, 7). Production of these antimicrobial peptides is widespread among bacterial species, and such production is made possible by biosynthetic machineries present in the genome, plasmids, or conjugative transposons (7). A conserved biosynthetic gene cluster for the production of bacteriocins displaying thiazole and oxazole heterocycles was discovered in 2008 in six microbial phyla (8). These gene clusters encode a toxin precursor and all indispensable proteins for toxin maturation in a mode similar to that associated with the bacteriocin microcin B17 (8). This gene cluster in L .monocytogenes was only present in a subset of lineage I strains...

SummaryMycobacterium abscessus is an emerging pathogen that is increasingly recognized as a relevant cause of human lung infection in cystic fibrosis patients. This highly antibiotic-resistant mycobacterium is an exception within the rapidly growing mycobacteria, which are mainly saprophytic and non-pathogenic organisms. M. abscessus manifests as either a smooth (S) or a rough (R) colony morphotype, which is of clinical importance as R morphotypes are associated with more severe and persistent infections. To better understand the molecular mechanisms behind the S/R alterations, we analysed S and R variants of three isogenic M. abscessus S/R pairs using an unbiased approach involving genome and transcriptome analyses, transcriptional fusions and integrating constructs. This revealed different small insertions, deletions (indels) or single nucleotide polymorphisms within the non-ribosomal peptide synthase gene cluster mps1-mps2-gap or mmpl4b in the three R variants, consistent with the transcriptional differences identified within this genomic locus that is implicated in the synthesis and transport of Glyco-Peptido-Lipids (GPL). In contrast to previous reports, the identification of clearly defined genetic lesions responsible for the loss of GPL-production or transport makes a frequent switching back-and-forth between smooth and rough morphologies in M. abscessus highly unlikely, which is important for our understanding of persistent M. abscessus infections.

Shigella invades the human intestinal mucosa, thus causing bacillary dysentery, an acute recto-colitis responsible for lethal complications, mostly in infants and toddlers. Conversely, commensal bacteria live in a mutualistic relationship with the intestinal mucosa that is characterized by homeostatic control of innate responses, thereby contributing to tolerance to the flora. Cross-talk established between commensals and the intestinal epithelium mediate this active process, the mechanisms of which remain largely uncharacterized. Probiotics such as Lactobacillus casei belong to a subclass of these commensals that modulate mucosal innate responses and possibly display anti-inflammatory properties. We analyzed whether L. casei could attenuate the pro-inflammatory signaling induced by Shigella flexneri after invasion of the epithelial lining. Cultured epithelial cells were infected with L. casei, followed by a challenge with S. flexneri. Using macroarray DNA chips, we observed that L. casei down-regulated the transcription of a number of genes encoding pro-inflammatory effectors such as cytokines and chemokines and adherence molecules induced by invasive S. flexneri. This resulted in an anti-inflammatory effect that appeared mediated by the inhibition of the NF-κB pathway, particularly through stabilization of I-κBα. In a time-course experiment using GeneChip hybridization analysis, the expression of many genes involved in ubiquitination and proteasome processes were modulated during L. casei treatment. Thus, L. casei has developed a sophisticated means to maintain intestinal homeostasis through a process that involves manipulation of the ubiquitin/proteasome pathway upstream of I-κBα.

The RIG-I-like receptors (RLRs) play a major role in sensing RNA virus infection to initiate and modulate antiviral immunity. They interact with particular viral RNAs, most of them being still unknown. To decipher the viral RNA signature on RLRs during viral infection, we tagged RLRs (RIG-I, MDA5, LGP2) and applied tagged protein affinity purification followed by next-generation sequencing (NGS) of associated RNA molecules. Two viruses with negative- and positive-sense RNA genome were used: measles (MV) and chikungunya (CHIKV). NGS analysis revealed that distinct regions of MV genome were specifically recognized by distinct RLRs: RIG-I recognized defective interfering genomes, whereas MDA5 and LGP2 specifically bound MV nucleoprotein-coding region. During CHIKV infection, RIG-I associated specifically to the 3’ untranslated region of viral genome. This study provides the first comparative view of the viral RNA ligands for RIG-I, MDA5 and LGP2 in the presence of infection.DOI: http://dx.doi.org/10.7554/eLife.11275.001

The RpoS/σS sigma subunit of RNA polymerase (RNAP) controls a global adaptive response that allows many Gram-negative bacteria to survive starvation and various stresses. σS also contributes to biofilm formation and virulence of the food-borne pathogen Salmonella enterica serovar Typhimurium (S. Typhimurium). In this study, we used directional RNA-sequencing and complementary assays to explore the σS-dependent transcriptome of S. Typhimurium during late stationary phase in rich medium. This study confirms the large regulatory scope of σS and provides insights into the physiological functions of σS in Salmonella. Extensive regulation by σS of genes involved in metabolism and membrane composition, and down-regulation of the respiratory chain functions, were important features of the σS effects on gene transcription that might confer fitness advantages to bacterial cells and/or populations under starving conditions. As an example, we show that arginine catabolism confers a competitive fitness advantage in stationary phase. This study also provides a firm basis for future studies to address molecular mechanisms of indirect regulation of gene expression by σS. Importantly, the σS-controlled downstream network includes small RNAs that might endow σS with post-transcriptional regulatory functions. Of these, four (RyhB-1/RyhB-2, SdsR, SraL) were known to be controlled by σS and deletion of the sdsR locus had a competitive fitness cost in stationary phase. The σS-dependent control of seven additional sRNAs was confirmed in Northern experiments. These findings will inspire future studies to investigate molecular mechanisms and the physiological impact of post-transcriptional regulation by σS.