Marie-Agnès Dillies

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.

Anti-inflammatory effect of Lactobacillus casei on Shigella-infected human intestinal epithelial cells.

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α.

Virulence strategies for infecting phagocytes deduced from the in vivo transcriptional program of Legionella pneumophila

Adaptation 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 (σ28) deletion mutant identified genes coregulated with the flagellar regulon, including GGDEF/EAL regulators and factors that promote host cell entry and survival.

The Global Regulator CodY Regulates Toxin Gene Expression in Bacillus anthracis and Is Required for Full Virulence

ABSTRACT In gram-positive bacteria, CodY is an important regulator of genes whose expression changes upon nutrient limitation and acts as a repressor of virulence gene expression in some pathogenic species. Here, we report the role of CodY in Bacillus anthracis, the etiologic agent of anthrax. Disruption of codY completely abolished virulence in a toxinogenic, noncapsulated strain, indicating that the activity of CodY is required for full virulence of B. anthracis. Global transcriptome analysis of a codY mutant and the parental strain revealed extensive differences. These differences could reflect direct control for some genes, as suggested by the presence of CodY binding sequences in their promoter regions, or indirect effects via the CodY-dependent control of other regulatory proteins or metabolic rearrangements in the codY mutant strain. The differences included reduced expression of the anthrax toxin genes in the mutant strain, which was confirmed by lacZ reporter fusions and immunoblotting. The accumulation of the global virulence regulator AtxA protein was strongly reduced in the mutant strain. However, in agreement with the microarray data, expression of atxA, as measured using an atxA-lacZ transcriptional fusion and by assaying atxA mRNA, was not significantly affected in the codY mutant. An atxA-lacZ translational fusion was also unaffected. Overexpression of atxA restored toxin component synthesis in the codY mutant strain. These results suggest that CodY controls toxin gene expression by regulating AtxA accumulation posttranslationally.

Identification and characterization of the genetic changes responsible for the characteristic smooth-to-rough morphotype alterations of clinically persistent Mycobacterium abscessus

Mycobacterium 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.

From array-based hybridization of Helicobacter pylori isolates to the complete genome sequence of an isolate associated with MALT lymphoma

Backgroundelicobacter pylori infection is associated with several gastro-duodenal inflammatory diseases of various levels of severity. To determine whether certain combinations of genetic markers can be used to predict the clinical source of the infection, we analyzed well documented and geographically homogenous clinical isolates using a comparative genomics approach.ResultsA set of 254 H. pylori genes was used to perform array-based comparative genomic hybridization among 120 French H. pylori strains associated with chronic gastritis (n = 33), duodenal ulcers (n = 27), intestinal metaplasia (n = 17) or gastric extra-nodal marginal zone B-cell MALT lymphoma (n = 43). Hierarchical cluster analyses of the DNA hybridization values allowed us to identify a homogeneous subpopulation of strains that clustered exclusively with cag PAI minus MALT lymphoma isolates. The genome sequence of B38, a representative of this MALT lymphoma strain-cluster, was completed, fully annotated, and compared with the six previously released H. pylori genomes (i.e. J99, 26695, HPAG1, P12, G27 and Shi470). B38 has the smallest H. pylori genome described thus far (1,576,758 base pairs containing 1,528 CDSs); it contains the vacA s2m2 allele and lacks the genes encoding the major virulence factors (absence of cag PAI, bab B, bab C, sab B, and hom B). Comparative genomics led to the identification of very few sequences that are unique to the B38 strain (9 intact CDSs and 7 pseudogenes). Pair-wise genomic synteny comparisons between B38 and the 6 H. pylori sequenced genomes revealed an almost complete co-linearity, never seen before between the genomes of strain Shi470 (a Peruvian isolate) and B38.ConclusionThese isolates are deprived of the main H. pylori virulence factors characterized previously, but are nonetheless associated with gastric neoplasia.

Dynamic RNA profiling in Plasmodium falciparum synchronized blood stages exposed to lethal doses of artesunate

BackgroundTranslation of the genome sequence of Plasmodium sp. into biologically relevant information relies on high through-put genomics technology which includes transcriptome analysis. However, few studies to date have used this powerful approach to explore transcriptome alterations of P. falciparum parasites exposed to antimalarial drugs.ResultsThe rapid action of artesunate allowed us to study dynamic changes of the parasite transcriptome in synchronous parasite cultures exposed to the drug for 90 minutes and 3 hours. Developmentally regulated genes were filtered out, leaving 398 genes which presented altered transcript levels reflecting drug-exposure. Few genes related to metabolic pathways, most encoded chaperones, transporters, kinases, Zn-finger proteins, transcription activating proteins, proteins involved in proteasome degradation, in oxidative stress and in cell cycle regulation. A positive bias was observed for over-expressed genes presenting a subtelomeric location, allelic polymorphism and encoding proteins with potential export sequences, which often belonged to subtelomeric multi-gene families. This pointed to the mobilization of processes shaping the interface between the parasite and its environment. In parallel, pathways were engaged which could lead to parasite death, such as interference with purine/pyrimidine metabolism, the mitochondrial electron transport chain, proteasome-dependent protein degradation or the integrity of the food vacuole.ConclusionThe high proportion of over-expressed genes encoding proteins exported from the parasite highlight the importance of extra-parasitic compartments as fields for exploration in drug research which, to date, has mostly focused on the parasite itself rather than on its intra and extra erythrocytic environment. Further work is needed to clarify which transcriptome alterations observed reflect a specific response to overcome artesunate toxicity or more general perturbations on the path to cellular death.

Docosahexaenoic Acid Inhibits Helicobacter pylori Growth In Vitro and Mice Gastric Mucosa Colonization

H. pylori drug-resistant strains and non-compliance to therapy are the major causes of H. pylori eradication failure. For some bacterial species it has been demonstrated that fatty acids have a growth inhibitory effect. Our main aim was to assess the ability of docosahexaenoic acid (DHA) to inhibit H. pylori growth both in vitro and in a mouse model. The effectiveness of standard therapy (ST) in combination with DHA on H. pylori eradication and recurrence prevention success was also investigated. The effects of DHA on H. pylori growth were analyzed in an in vitro dose-response study and n in vivo model. We analized the ability of H. pylori to colonize mice gastric mucosa following DHA, ST or a combination of both treatments. Our data demonstrate that DHA decreases H. pylori growth in vitro in a dose-dependent manner. Furthermore, DHA inhibits H. pylori gastric colonization in vivo as well as decreases mouse gastric mucosa inflammation. Addition of DHA to ST was also associated with lower H. pylori infection recurrence in the mouse model. In conclusion, DHA is an inhibitor of H. pylori growth and its ability to colonize mouse stomach. DHA treatment is also associated with a lower recurrence of H. pylori infection in combination with ST. These observations pave the way to consider DHA as an adjunct agent in H. pylori eradication treatment.

Temporal transcriptomic response during arsenic stress in Herminiimonas arsenicoxydans.

BackgroundArsenic is present in numerous ecosystems and microorganisms have developed various mechanisms to live in such hostile environments. Herminiimonas arsenicoxydans, a bacterium isolated from arsenic contaminated sludge, has acquired remarkable capabilities to cope with arsenic. In particular our previous studies have suggested the existence of a temporal induction of arsenite oxidase, a key enzyme in arsenic metabolism, in the presence of As(III).ResultsMicroarrays were designed to compare gene transcription profiles under a temporal As(III) exposure. Transcriptome kinetic analysis demonstrated the existence of two phases in arsenic response. The expression of approximatively 14% of the whole genome was significantly affected by an As(III) early stress and 4% by an As(III) late exposure. The early response was characterized by arsenic resistance, oxidative stress, chaperone synthesis and sulfur metabolism. The late response was characterized by arsenic metabolism and associated mechanisms such as phosphate transport and motility. The major metabolic changes were confirmed by chemical, transcriptional, physiological and biochemical experiments. These early and late responses were defined as general stress response and specific response to As(III), respectively.ConclusionGene expression patterns suggest that the exposure to As(III) induces an acute response to rapidly minimize the immediate effects of As(III). Upon a longer arsenic exposure, a broad metabolic response was induced. These data allowed to propose for the first time a kinetic model of the As(III) response in bacteria.

Endoplasmic Reticulum Stress-Sensing Mechanism Is Activated in Entamoeba histolytica upon Treatment with Nitric Oxide

The Endoplasmic Reticulum stores calcium and is a site of protein synthesis and modification. Changes in ER homeostasis lead to stress responses with an activation of the unfolded protein response (UPR). The Entamoeba histolytica endomembrane system is simple compared to those of higher eukaryotes, as a canonical ER is not observed. During amoebiasis, an infection of the human intestine and liver by E. histolytica, nitric oxide (NO) triggers an apoptotic-like event preceded by an impairment of energy production and a loss of important parasite pathogenic features. We address the question of how this ancient eukaryote responds to stress induced by immune components (i.e. NO) and whether stress leads to ER changes and subsequently to an UPR. Gene expression analysis suggested that NO triggers stress responses marked by (i) dramatic up-regulation of hsp genes although a bona fide UPR is absent; (ii) induction of DNA repair and redox gene expression and iii) up-regulation of glycolysis-related gene expression. Enzymology approaches demonstrate that NO directly inhibits glycolysis and enhance cysteine synthase activity. Using live imaging and confocal microscopy we found that NO dramatically provokes extensive ER fragmentation. ER fission in E. histolytica appears as a protective response against stress, as it has been recently proposed for neuron self-defense during neurologic disorders. Chronic ER stress is also involved in metabolic diseases including diabetes, where NO production reduces ER calcium levels and activates cell death. Our data highlighted unique cellular responses of interest to understand the mechanisms of parasite death during amoebiasis.