Search

Abstract Background Arsenic 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). Results Microarrays 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. Conclusion Gene 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.

Abstract Background elicobacter 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. Results A 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 cagPAI 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 vacAs2m2 allele and lacks the genes encoding the major virulence factors (absence of cagPAI, babB, babC, sabB, and homB). 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. Conclusion These isolates are deprived of the main H. pylori virulence factors characterized previously, but are nonetheless associated with gastric neoplasia.

Abstract Background Both the speciation and toxicity of arsenic are affected by bacterial transformations, i.e. oxidation, reduction or methylation. These transformations have a major impact on environmental contamination and more particularly on arsenic contamination of drinking water. Herminiimonas arsenicoxydans has been isolated from an arsenic- contaminated environment and has developed various mechanisms for coping with arsenic, including the oxidation of As(III) to As(V) as a detoxification mechanism. Results In the present study, a differential transcriptome analysis was used to identify genes, including arsenite oxidase encoding genes, involved in the response of H. arsenicoxydans to As(III). To get insight into the molecular mechanisms of this enzyme activity, a Tn5 transposon mutagenesis was performed. Transposon insertions resulting in a lack of arsenite oxidase activity disrupted aoxR and aoxS genes, showing that the aox operon transcription is regulated by the AoxRS two-component system. Remarkably, transposon insertions were also identified in rpoN coding for the alternative N sigma factor (σ54) of RNA polymerase and in dnaJ coding for the Hsp70 co-chaperone. Western blotting with anti-AoxB antibodies and quantitative RT-PCR experiments allowed us to demonstrate that the rpoN and dnaJ gene products are involved in the control of arsenite oxidase gene expression. Finally, the transcriptional start site of the aoxAB operon was determined using rapid amplification of cDNA ends (RACE) and a putative -12/-24 σ54-dependent promoter motif was identified upstream of aoxAB coding sequences. Conclusion These results reveal the existence of novel molecular regulatory processes governing arsenite oxidase expression in H. arsenicoxydans. These data are summarized in a model that functionally integrates arsenite oxidation in the adaptive response to As(III) in this microorganism.

Abstract Background A number of molecular tools have been developed to monitor the emergence and spread of anti-malarial drug resistance to Plasmodium falciparum. One of the major obstacles to the wider implementation of these tools is the absence of practical methods enabling high throughput analysis. Here a new Zip-code array is described, called FlexiChip, linked to a dedicated software program, which largely overcomes this problem. Methods Previously published microarray probes detecting single-nucleotide polymorphisms (SNP) associated with parasite resistance to anti-malarial drugs (ResMalChip) were adapted for a universal microarray FlexiChip format. To evaluate the overall sensitivity of the FlexiChip package (microarray + software), the results of FlexiChip were compared to ResMalChip microarray, using the same extension probes and with the same PCR products. In both cases, sequence results were used as gold standard to calculate sensitivity and specificity. FlexiChip results obtained with a set of field isolates were then compared to those assessed in an independent reference laboratory. Results The FlexiChip package gave results identical to the ResMalChip results in 92.7% of samples (kappa coefficient 0.8491, with a standard error 0.021) and had a sensitivity of 95.88% and a specificity of 97.68% compared to the sequencing as the reference method. Moreover the method performed well compared to the results obtained in the reference laboratories, with 99.7% of identical results (kappa coefficient 0.9923, S.E. 0.0523). Conclusion Microarrays could be employed to monitor P. falciparum drug resistance markers with greater cost effectiveness and the possibility for high throughput analysis. The FlexiChip package is a promising tool for use in poor resource settings of malaria endemic countries.

Abstract Background In man, infection by the Gram-negative enteropathogen Yersinia pseudotuberculosis is usually limited to the terminal ileum. However, in immunocompromised patients, the microorganism may disseminate from the digestive tract and thus cause a systemic infection with septicemia. Results To gain insight into the metabolic pathways and virulence factors expressed by the bacterium at the blood stage of pseudotuberculosis, we compared the overall gene transcription patterns (the transcriptome) of bacterial cells cultured in either human plasma or Luria-Bertani medium. The most marked plasma-triggered metabolic consequence in Y. pseudotuberculosis was the switch to high glucose consumption, which is reminiscent of the acetogenic pathway (known as "glucose overflow") in Escherichia coli. However, upregulation of the glyoxylate shunt enzymes suggests that (in contrast to E. coli) acetate may be further metabolized in Y. pseudotuberculosis. Our data also indicate that the bloodstream environment can regulate major virulence genes (positively or negatively); the yadA adhesin gene and most of the transcriptional units of the pYV-encoded type III secretion apparatus were found to be upregulated, whereas transcription of the pH6 antigen locus was strongly repressed. Conclusion Our results suggest that plasma growth of Y. pseudotuberculosis is responsible for major transcriptional regulatory events and prompts key metabolic reorientations within the bacterium, which may in turn have an impact on virulence.

Abstract Background Translation 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. Results The 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. Conclusion The 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.

Background: Advances in high-throughput technologies have originated an ever-increasing availability of omics datasets. The integration of multiple heterogeneous data sources is currently an issue for biology and bioinformatics. Multiple kernel learning (MKL) has shown to be a flexible and valid approach to consider the diverse nature of multi-omics inputs, despite being an underused tool in genomic data mining., Results: We provide novel MKL approaches based on different kernel fusion strategies. To learn from the meta-kernel of input kernels, we adapted unsupervised integration algorithms for supervised tasks with support vector machines. We also tested deep learning architectures for kernel fusion and classification. The results show that MKL-based models can outperform more complex, state-of-the-art, supervised multi-omics integrative approaches., Conclusion: Multiple kernel learning offers a natural framework for predictive models in multi-omics data. It proved to provide a fast and reliable solution that can compete with and outperform more complex architectures. Our results offer a direction for bio-data mining research, biomarker discovery and further development of methods for heterogeneous data integration., Competing Interests: Declarations. Ethics approval and consent to participate: Not applicable. Consent for publication: Not applicable. Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Background: Kernel methods have been proven to be a powerful tool for the integration and analysis of high-throughput technologies generated data. Kernels offer a nonlinear version of any linear algorithm solely based on dot products. The kernelized version of principal component analysis is a valid nonlinear alternative to tackle the nonlinearity of biological sample spaces. This paper proposes a novel methodology to obtain a data-driven feature importance based on the kernel PCA representation of the data., Results: The proposed method, kernel PCA Interpretable Gradient (KPCA-IG), provides a data-driven feature importance that is computationally fast and based solely on linear algebra calculations. It has been compared with existing methods on three benchmark datasets. The accuracy obtained using KPCA-IG selected features is equal to or greater than the other methods' average. Also, the computational complexity required demonstrates the high efficiency of the method. An exhaustive literature search has been conducted on the selected genes from a publicly available Hepatocellular carcinoma dataset to validate the retained features from a biological point of view. The results once again remark on the appropriateness of the computed ranking., Conclusions: The black-box nature of kernel PCA needs new methods to interpret the original features. Our proposed methodology KPCA-IG proved to be a valid alternative to select influential variables in high-dimensional high-throughput datasets, potentially unravelling new biological and medical biomarkers., (© 2023. The Author(s).)

Entamoeba is a genus of Amoebozoa that includes the intestine-colonizing pathogenic species Entamoeba histolytica . To understand the basis of gene regulation in E. histolytica from an evolutionary perspective, we have profiled the transcriptomes of its closely related species E. dispar , E. moshkovskii and E. invadens. Genome-wide identification of transcription start sites (TSS) and polyadenylation sites (PAS) revealed the similarities and differences of their gene regulatory sequences. In particular, we found the widespread initiation of antisense transcription from within the gene coding sequences is a common feature among all Entamoeba species. Interestingly, we observed the enrichment of antisense transcription in genes involved in several processes that are common to species infecting the human intestine, e.g., the metabolism of phospholipids. These results suggest a potentially conserved and compact gene regulatory system in Entamoeba .

Background: Epstein-Barr virus (EBV) targets B-cells where it establishes a latent infection. EBV can transform B-cells in vitro and is recognized as an oncogenic virus, especially in the setting of immune compromise. Indeed, immunodeficient patients may fail to control chronic EBV infection, leading to the development EBV-driven lymphoid malignancies. Ataxia telangiectasia (AT) is a primary immune deficiency caused by mutations in the ATM gene, involved in the repair of double-strand breaks. Patients with AT are at high risk of developing cancers, mostly B-cell lymphoid malignancies, most of which being EBV-related. Aside from immune deficiency secondary to AT, loss of ATM function could also hinder the control of the virus within B-cells, favoring lymphomagenesis in AT patients., Results: We used RNA sequencing on lymphoblastoid cell lines derived from patients with AT and healthy donors to analyze and compare both cellular and viral gene expression. We found numerous deregulated signaling pathways involving transcription, translation, oncogenesis and immune regulation. Specifically, the translational defect was confirmed in vitro, suggesting that the pathogenesis of AT may also involve a ribosomal defect. Concomitant analysis of viral gene expression did not reveal significant differential gene expression, however, analysis of EBV interactome suggests that the viral latency genes EBNA-3A, EBNA-3C and LMP1 may be disrupted in LCL from AT patients., Conclusion: Our data support the notion that ATM deficiency deregulates cellular gene expression possibly disrupting interactions with EBV latent genes, promoting the oncogenic potential of the virus. These preliminary findings provide a new step towards the understanding of EBV regulation and of AT pathogenesis.

Lassa virus (LASV) is endemic in West Africa and induces a viral hemorrhagic fever (VHF) with up to 30% lethality among clinical cases. The mechanisms involved in control of Lassa fever or, in contrast, the ensuing catastrophic illness and death are poorly understood. We used the cynomolgus monkey model to reproduce the human disease with asymptomatic to mild or fatal disease. After initial replication at the inoculation site, LASV reached the secondary lymphoid organs. LASV did not spread further in nonfatal disease and was rapidly controlled by balanced innate and T-cell responses. Systemic viral dissemination occurred during severe disease. Massive replication, a cytokine/chemokine storm, defective T-cell responses, and multiorgan failure were observed. Clinical, biological, immunological, and transcriptomic parameters resembled those observed during septic-shock syndrome, suggesting that similar pathogenesis is induced during Lassa fever. The outcome appears to be determined early, as differentially expressed genes in PBMCs were associated with fatal and non-fatal Lassa fever outcome very early after infection. These results provide a full characterization and important insights into Lassa fever pathogenesis and could help to develop early diagnostic tools.

Natural antisense transcripts (NAT) have been reported in prokaryotes and eukaryotes. While the functions of most reported NATs remain unknown, their potentials in regulating the transcription of their counterparts have been speculated. Entamoeba histolytica, which is a unicellular eukaryotic parasite, has a compact protein-coding genome with very short intronic and intergenic regions. The regulatory mechanisms of gene expression in this compact genome are under-described. In this study, by genome-wide mapping of RNA-Seq data in the genome of E. histolytica, we show that a substantial fraction of its protein-coding genes (28%) has significant transcription on their opposite strand (i.e. NAT). Intriguingly, we found the location of transcription start sites or polyadenylation sites of NAT are determined by the specific motifs encoded on the opposite strand of the gene coding sequences, thereby providing a compact regulatory system for gene transcription. Moreover, we demonstrated that NATs are globally up-regulated under various environmental conditions including temperature stress and pathogenicity. While NATs do not appear to be consequences of spurious transcription, they may play a role in regulating gene expression in E. histolytica, a hypothesis which needs to be tested.

The intestinal microbiota modulates host physiology and gene expression via mechanisms that are not fully understood. Here we examine whether host epitranscriptomic marks are affected by the gut microbiota. We use methylated RNA-immunoprecipitation and sequencing (MeRIP-seq) to identify N6-methyladenosine (m 6 A) modifications in mRNA of mice carrying conventional, modified, or no microbiota. We find that variations in the gut microbiota correlate with m 6 A modifications in the cecum, and to a lesser extent in the liver, affecting pathways related to metabolism, inflammation and antimicrobial responses. We analyze expression levels of several known writer and eraser enzymes, and find that the methyltransferase Mettl16 is downregulated in absence of a microbiota, and one of its target mRNAs, encoding S-adenosylmethionine synthase Mat2a, is less methylated. We furthermore show that Akkermansia muciniphila and Lactobacillus plantarum affect specific m 6 A modifications in mono-associated mice. Our results highlight epitranscriptomic modifications as an additional level of interaction between commensal bacteria and their host.

Lassa fever is a major threat in Western Africa. The large number of people living at risk for this disease calls for the development of a vaccine against Lassa virus (LASV). We generated live-attenuated LASV vaccines based on measles virus and Mopeia virus platforms and expressing different LASV antigens, with the aim to develop a vaccine able to protect after a single shot. We compared the efficacy of these vaccines against LASV in cynomolgus monkeys. The vaccines were well tolerated and protected the animals from LASV infection and disease after a single immunization but with varying efficacy. Analysis of the immune responses showed that complete protection was associated with robust secondary T cell and antibody responses against LASV. Transcriptomic and proteomic analyses showed an early activation of innate immunity and T cell priming after immunization with the most effective vaccines, with changes detectable as early as 2 days after immunization. The most efficacious vaccine candidate, a measles vector simultaneously expressing LASV glycoprotein and nucleoprotein, has been selected for further clinical evaluation., (Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Spontaneous control of human immunodeficiency virus (HIV) is generally associated with an enhanced capacity of CD8 + T cells to eliminate infected CD4 + T cells, but the molecular characteristics of these highly functional CD8 + T cells are largely unknown. In the present study, using single-cell analysis, it was shown that HIV-specific, central memory CD8 + T cells from spontaneous HIV controllers (HICs) and antiretrovirally treated non-controllers have opposing transcriptomic profiles. Genes linked to effector functions and survival are upregulated in cells from HICs. In contrast, genes associated with activation, exhaustion and glycolysis are upregulated in cells from non-controllers. It was shown that HIV-specific CD8 + T cells from non-controllers are largely glucose dependent, whereas those from HICs have more diverse metabolic resources that enhance both their survival potential and their capacity to develop anti-HIV effector functions. The functional efficiency of the HIV-specific CD8 + T cell response in HICs is thus engraved in their memory population and related to their metabolic programme. Metabolic reprogramming in vitro through interleukin-15 treatment abrogated the glucose dependency and enhanced the antiviral potency of HIV-specific CD8 + T cells from non-controllers.

Lassa virus (LASV) causes a viral haemorrhagic fever in humans and is a major public health concern in West Africa. An efficient immune response to LASV appears to rely on type I interferon (IFN-I) production and T-cell activation. We evaluated the response of plasmacytoid dendritic cells (pDC) to LASV, as they are an important and early source of IFN-I. We compared the response of primary human pDCs to LASV and Mopeia virus (MOPV), which is very closely related to LASV, but non-pathogenic. We showed that pDCs are not productively infected by either MOPV or LASV, but produce IFN-I. However, the activation of pDCs was more robust in response to MOPV than LASV. In vivo, pDC activation may support the control of viral replication through IFN-I production, but also improve the induction of a global immune response. Therefore, pDC activation could play a role in the control of LASV infection.

Background: Infectious diseases are still a leading cause of death and, with the emergence of drug resistance, pose a great threat to human health. New drugs and strategies are thus urgently needed to improve treatment efficacy and limit drug-associated side effects. Nanotechnology-based drug delivery systems are promising approaches, offering hope in the fight against drug resistant bacteria. However, how nanocarriers influence the response of innate immune cells to bacterial infection is mostly unknown., Results: Here, we used Mycobacterium tuberculosis as a model of bacterial infection to examine the impact of mannose functionalization of chitosan nanocarriers (CS-NCs) on the human macrophage response. Both ungrafted and grafted CS-NCs were similarly internalized by macrophages, via an actin cytoskeleton-dependent process. Although tri-mannose ligands did not modify the capacity of CS-NCs to escape lysosomal degradation, they profoundly remodeled the response of M. tuberculosis-infected macrophages. mRNA sequencing showed nearly 900 genes to be differentially expressed due to tri-mannose grafting. Unexpectedly, the set of modulated genes was enriched for pathways involved in cell metabolism, particularly oxidative phosphorylation and sugar metabolism., Conclusions: The ability to modulate cell metabolism by grafting ligands at the surface of nanoparticles may thus be a promising strategy to reprogram immune cells and improve the efficacy of encapsulated drugs.

Lassa virus (LASV) is responsible for a viral hemorrhagic fever in humans and the death of 3,000 to 5,000 people every year. The immune response to LASV is poorly understood, but type I interferon (IFN-I) and T-cell responses appear to be critical for the host. We studied the response of myeloid dendritic cells (mDC) to LASV, as mDCs are involved in both IFN-I production and T-cell activation. We compared the response of primary human mDCs to LASV and Mopeia virus (MOPV), which is similar to LASV, but non-pathogenic. We showed that mDCs produced substantial amounts of IFN-I in response to both LASV and MOPV. However, only MOPV-infected mDCs were able to activate T cells. More surprisingly, coculture with T cells completely inhibited the activation of LASV-infected mDCs. These differences between LASV and MOPV were mostly due to the LASV nucleoprotein, which has major immunosuppressive properties, but the glycoprotein was also involved. Overall, these results suggest that mDCs may be important for the global response to LASV and play a role in the outcome of Lassa fever., Competing Interests: The authors have declared that no competing interests exist.

Background: The SOS response is an almost ubiquitous response of cells to genotoxic stresses. The full complement of genes in the SOS regulon for Vibrio species has only been addressed through bioinformatic analyses predicting LexA binding box consensus and in vitro validation. Here, we perform whole transcriptome sequencing from Vibrio cholerae treated with mitomycin C as an SOS inducer to characterize the SOS regulon and other pathways affected by this treatment., Results: Comprehensive transcriptional profiling allowed us to define the full landscape of promoters and transcripts active in V. cholerae. We performed extensive transcription start site (TSS) mapping as well as detection/quantification of the coding and non-coding RNA (ncRNA) repertoire in strain N16961. To improve TSS detection, we developed a new technique to treat RNA extracted from cells grown in various conditions. This allowed for identification of 3078 TSSs with an average 5'UTR of 116 nucleotides, and peak distribution between 16 and 64 nucleotides; as well as 629 ncRNAs. Mitomycin C treatment induced transcription of 737 genes and 28 ncRNAs at least 2 fold, while it repressed 231 genes and 17 ncRNAs. Data analysis revealed that in addition to the core genes known to integrate the SOS regulon, several metabolic pathways were induced. This study allowed for expansion of the Vibrio SOS regulon, as twelve genes (ubiEJB, tatABC, smpA, cep, VC0091, VC1190, VC1369-1370) were found to be co-induced with their adjacent canonical SOS regulon gene(s), through transcriptional read-through. Characterization of UV and mitomycin C susceptibility for mutants of these newly identified SOS regulon genes and other highly induced genes and ncRNAs confirmed their role in DNA damage rescue and protection., Conclusions: We show that genotoxic stress induces a pervasive transcriptional response, affecting almost 20% of the V. cholerae genes. We also demonstrate that the SOS regulon is larger than previously known, and its syntenic organization is conserved among Vibrio species. Furthermore, this specific co-localization is found in other γ-proteobacteria for genes recN-smpA and rmuC-tatABC, suggesting SOS regulon conservation in this phylum. Finally, we comment on the limitations of widespread NGS approaches for identification of all RNA species in bacteria.

Dengue virus (DENV) causes more human infections than any other mosquito-borne virus. The current lack of antiviral strategies has prompted genome-wide screens for host genes that are required for DENV infectivity. Earlier transcriptomic studies that identified DENV host factors in the primary vector Aedes aegypti used inbred laboratory colonies and/or pools of mosquitoes that erase individual variation. Here, we performed transcriptome sequencing on individual midguts in a field-derived Ae. aegypti population to identify new candidate host factors modulating DENV replication. We analyzed the transcriptomic data using an approach that accounts for individual co-variation between viral RNA load and gene expression. This approach generates a prediction about the agonist or antagonist effect of candidate genes on DENV replication based on the sign of the correlation between gene expression and viral RNA load. Using this method, we identified 39 candidate genes that went undetected by conventional pairwise comparison of gene expression levels between DENV-infected midguts and uninfected controls. Only four candidate genes were detected by both methods, emphasizing their complementarity. We demonstrated the value of our approach by functional validation of a candidate agonist gene encoding a sterol regulatory element-binding protein (SREBP), which was identified by correlation analysis but not by pairwise comparison. We confirmed that SREBP promotes DENV infection in the midgut by RNAi-mediated gene knockdown in vivo. We suggest that our approach for transcriptomic analysis can empower genome-wide screens for potential agonist or antagonist factors by leveraging inter-individual variation in gene expression. More generally, this method is applicable to a wide range of phenotypic traits displaying inter-individual variation.

Bacteria of the genus Thiomonas are found ubiquitously in arsenic contaminated waters such as acid mine drainage (AMD), where they contribute to the precipitation and the natural bioremediation of arsenic. In these environments, these bacteria have developed a large range of resistance strategies among which the capacity to form particular biofilm structures. The biofilm formation is one of the most ubiquitous adaptive response observed in prokaryotes to various stresses, such as those induced in the presence of toxic compounds. This study focused on the process of biofilm formation in three Thiomonas strains (CB1, CB2 and CB3) isolated from the same AMD. The results obtained here show that these bacteria are all capable of forming biofilms, but the architecture and the kinetics of formation of these biofilms differ depending on whether arsenite is present in the environment and from one strain to another. Indeed, two strains favoured biofilm formation, whereas one favoured motility in the presence of arsenite. To identify the underlying mechanisms, the patterns of expression of some genes possibly involved in the process of biofilm formation were investigated in Thiomonas sp. CB2 in the presence and absence of arsenite, using a transcriptomic approach (RNA-seq). The findings obtained here shed interesting light on how the formation of biofilms, and the motility processes contribute to the adaptation of Thiomonas strains to extreme environments., (© 2016 The Authors. Microbial Biotechnology published by John Wiley & Sons Ltd and Society for Applied Microbiology.)

Motivation: With the continued improvement of requisite mass spectrometers and UHPLC systems, Hydrogen/Deuterium eXchange Mass Spectrometry (HDX-MS) workflows are rapidly evolving towards the investigation of more challenging biological systems, including large protein complexes and membrane proteins. The analysis of such extensive systems results in very large HDX-MS datasets for which specific analysis tools are required to speed up data validation and interpretation., Results: We introduce a web application and a new R-package named 'MEMHDX' to help users analyze, validate and visualize large HDX-MS datasets. MEMHDX is composed of two elements. A statistical tool aids in the validation of the results by applying a mixed-effects model for each peptide, in each experimental condition, and at each time point, taking into account the time dependency of the HDX reaction and number of independent replicates. Two adjusted P-values are generated per peptide, one for the 'Change in dynamics' and one for the 'Magnitude of ΔD', and are used to classify the data by means of a 'Logit' representation. A user-friendly interface developed with Shiny by RStudio facilitates the use of the package. This interactive tool allows the user to easily and rapidly validate, visualize and compare the relative deuterium incorporation on the amino acid sequence and 3D structure, providing both spatial and temporal information., Availability and Implementation: MEMHDX is freely available as a web tool at the project home page http://memhdx.c3bi.pasteur.fr CONTACT: marie-agnes.dillies@pasteur.fr or sebastien.brier@pasteur.frSupplementary information: Supplementary data is available at Bioinformatics online., (© The Author 2016. Published by Oxford University Press.)

Amoebiasis is a human infectious disease due to the amoeba parasite Entamoeba histolytica. The disease appears in only 20% of the infections. Diversity in phenotypes may occur within the same infectious strain in the gut; for instance, parasites can be commensal (in the intestinal lumen) or pathogenic (inside the tissue). The degree of pathogenesis of clinical isolates varies greatly. These findings raise the hypothesis that genetic derivation may account for amoebic diverse phenotypes. The main goal of this study was to analyse gene expression changes of a single virulent amoebic strain in different environmental contexts where it exhibit different degrees of virulence, namely isolated from humans and maintained through animal liver passages, in contact with the human colon and short or prolonged in vitro culture. The study reveals major transcriptome changes in virulent parasites upon contact with human colon explants, including genes related to sugar metabolism, cytoskeleton rearrangement, stress responses and DNA repair. Furthermore, in long-term cultured parasites, drastic changes in gene expression for proteins with functions for proteasome and tRNA activities were found. Globally we conclude that rapid changes in gene expression rather than genetic derivation can sustain the invasive phenotype of a single virulent isolate of E. histolytica.

The fungal cell wall is a rigid structure because of fibrillar and branched β-(1,3)-glucan linked to chitin. Softening of the cell wall is an essential phenomenon during fungal morphogenesis, wherein rigid cell wall structures are cleaved by glycosylhydrolases. During the search for glycosylhydrolases acting on β-(1,3)-glucan, we identified seven genes in the Aspergillus fumigatus genome coding for potential endo-β-(1,3)-glucanase. ENG1 (previously characterized and named ENGL1, Mouyna et al., ), belongs to the Glycoside-Hydrolase 81 (GH81) family, while ENG2 to ENG7, to GH16 family. ENG1 and four GH16 genes (ENG2-5) were expressed in the resting conidia as well as during germination, suggesting an essential role during A. fumigatus morphogenesis. Here, we report the effect of sequential deletion of AfENG2-5 (GH16) followed by AfENG1 (GH81) deletion in the Δeng2,3,4,5 mutant. The Δeng1,2,3,4,5 mutant showed conidial defects, with linear chains of conidia unable to separate while the germination rate was not affected. These results show, for the first time in a filamentous fungus, that endo β-(1,3)-glucanases are essential for proper conidial cell wall assembly and thus segregation of conidia during conidiation., (© 2016 John Wiley & Sons Ltd.)

Uropathogenic Escherichia coli (UPEC) is the leading cause of urinary tract infections (UTIs) worldwide, causing over 150 million clinical cases annually. There is currently no specific treatment addressing the asymptomatic carriage in the gut of UPEC before they initiate UTIs. This study investigates the efficacy of virulent bacteriophages to decrease carriage of gut pathogens. Three virulent bacteriophages infecting an antibiotic-resistant UPEC strain were isolated and characterized both in vitro and in vivo. A new experimental murine model of gut carriage of E. coli was elaborated and the impact of virulent bacteriophages on colonization levels and microbiota diversity was assessed. A single dose of a cocktail of the three bacteriophages led to a sharp decrease in E. coli levels throughout the gut. We also observed that microbiota diversity was much less affected by bacteriophages than by antibiotics. Therefore, virulent bacteriophages can efficiently target UPEC strains residing in the gut, with potentially profound public health and economic impacts. These results open a new area with the possibility to manipulate specifically the microbiota using virulent bacteriophages, which could have broad applications in many gut-related disorders/diseases and beyond., (© 2016 Society for Applied Microbiology and John Wiley & Sons Ltd.)

Background: Several 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., Results: SARTools 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., Conclusions: SARTools 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.

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.

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.

BAH domain-containing protein 1 (BAHD1) is involved in heterochromatin formation and gene repression in human cells. BAHD1 also localizes to the inactive X chromosome (Xi), but the functional significance of this targeting is unknown. So far, research on this protein has been hampered by its low endogenous abundance and its role in epigenetic regulation remains poorly explored. In this work, we used whole-genome bisulfite sequencing (BS-seq) to compare the DNA methylation profile of HEK293 cells expressing low levels of BAHD1 (HEK-CT) to that of isogenic cells stably overexpressing BAHD1 (HEK-BAHD1). We show that increasing BAHD1 levels induces de novo DNA methylation on autosomes and a marked hypomethylation on the X chromosome (chrX). We identified 91,358 regions that have different methylation patterns in HEK-BAHD1 compared to HEK-CT cells (termed "BAHD1-DMRs"), of which 83,850 mapped on autosomes and 7508 on the X chromosome (chrX). Autosomal BAHD1-DMRs were predominantly hypermethylated and located to satellites, interspersed repeats, and intergenic regions. In contrast, BAHD1-DMRs on chrX were mainly hypomethylated and located to gene bodies and enhancers. We further found that BAHD1-DMRs display a higher-order organization by being clustered within large chromosomal domains. Half of these "BAHD1-Associated differentially methylated Domains" (BADs) overlapped with lamina-associated domains (LADs). Based on these results, we propose that BAHD1-mediated heterochromatin formation is linked to DNA methylation and may play a role in the spatial architecture of the genome.

Background and Aims: In cystic fibrosis (CF), Pseudomonas aeruginosa is not eradicated from the lower respiratory tract and is associated with epithelial inflammation that eventually causes tissue damage. To identify the molecular determinants of an effective response to P. aeruginosa infection, we performed a transcriptomic analysis of primary human bronchial epithelial cells from healthy donors (CTRL) 2, 4, and 6 h after induced P. aeruginosa infection. Compared to noninfected cells, infected cells showed changes in gene activity, which were most marked 6 h postinfection and usually consisted in upregulation., Results: By comparing for each time point of infection, the transcriptomic response of epithelial cells from CF patients and healthy donors, we identified 851, 638, 667, and 980 differentially expressed genes 0, 2, 4, and 6 h postinfection, respectively. Gene selection followed by bioinformatic analysis showed that most of the differentially expressed genes, either up- or downregulated, were in the protein-binding and catalytic gene-ontology categories. Finally, we established that the protein products of the genes exhibiting the greatest differential upregulation (CSF2, CCL2, TNF, CSF3, MMP1, and MMP10) between CF patients and CTRL were produced in higher amounts by infected cells from CF patients versus CTRL., Conclusions: The differentially expressed genes in CF patients may constitute a signature for a detrimental inflammatory response and for an inefficient P. aeruginosa host-cell response.

Acid mine drainage (AMD) is a highly toxic environment for most living organisms due to the presence of many lethal elements including arsenic (As). Thiomonas (Tm.) bacteria are found ubiquitously in AMD and can withstand these extreme conditions, in part because they are able to oxidize arsenite. In order to further improve our knowledge concerning the adaptive capacities of these bacteria, we sequenced and assembled the genome of six isolates derived from the Carnoulès AMD, and compared them to the genomes of Tm. arsenitoxydans 3As (isolated from the same site) and Tm. intermedia K12 (isolated from a sewage pipe). A detailed analysis of the Tm. sp. CB2 genome revealed various rearrangements had occurred in comparison to what was observed in 3As and K12 and over 20 genomic islands (GEIs) were found in each of these three genomes. We performed a detailed comparison of the two arsenic-related islands found in CB2, carrying the genes required for arsenite oxidation and As resistance, with those found in K12, 3As, and five other Thiomonas strains also isolated from Carnoulès (CB1, CB3, CB6, ACO3 and ACO7). Our results suggest that these arsenic-related islands have evolved differentially in these closely related Thiomonas strains, leading to divergent capacities to survive in As rich environments.

Background: Streptococcus agalactiae, or Group B Streptococcus, is a leading cause of neonatal infections and an increasing cause of infections in adults with underlying diseases. In an effort to reconstruct the transcriptional networks involved in S. agalactiae physiology and pathogenesis, we performed an extensive and robust characterization of its transcriptome through a combination of differential RNA-sequencing in eight different growth conditions or genetic backgrounds and strand-specific RNA-sequencing., Results: Our study identified 1,210 transcription start sites (TSSs) and 655 transcript ends as well as 39 riboswitches and cis-regulatory regions, 39 cis-antisense non-coding RNAs and 47 small RNAs potentially acting in trans. Among these putative regulatory RNAs, ten were differentially expressed in response to an acid stress and two riboswitches sensed directly or indirectly the pH modification. Strikingly, 15% of the TSSs identified were associated with the incorporation of pseudo-templated nucleotides, showing that reiterative transcription is a pervasive process in S. agalactiae. In particular, 40% of the TSSs upstream genes involved in nucleotide metabolism show reiterative transcription potentially regulating gene expression, as exemplified for pyrG and thyA encoding the CTP synthase and the thymidylate synthase respectively., Conclusions: This comprehensive map of the transcriptome at the single nucleotide resolution led to the discovery of new regulatory mechanisms in S. agalactiae. It also provides the basis for in depth analyses of transcriptional networks in S. agalactiae and of the regulatory role of reiterative transcription following variations of intra-cellular nucleotide pools.

The RpoS/σ(S) sigma subunit of RNA polymerase (RNAP) activates transcription of stationary phase genes in many Gram-negative bacteria and controls adaptive functions, including stress resistance, biofilm formation and virulence. In this study, we address an important but poorly understood aspect of σ(S)-dependent control, that of a repressor. Negative regulation by σ(S) has been proposed to result largely from competition between σ(S) and other σ factors for binding to a limited amount of core RNAP (E). To assess whether σ(S) binding to E alone results in significant downregulation of gene expression by other σ factors, we characterized an rpoS mutant of Salmonella enterica serovar Typhimurium producing a σ(S) protein proficient for Eσ(S) complex formation but deficient in promoter DNA binding. Genome expression profiling and physiological assays revealed that this mutant was defective for negative regulation, indicating that gene repression by σ(S) requires its binding to DNA. Although the mechanisms of repression by σ(S) are likely specific to individual genes and environmental conditions, the study of transcription downregulation of the succinate dehydrogenase operon suggests that σ competition at the promoter DNA level plays an important role in gene repression by Eσ(S)., (© The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The mechanisms underlying the heterogeneity of clinical malaria remain largely unknown. We hypothesized that differential gene expression contributes to phenotypic variation of parasites which results in a specific interaction with the host, leading to different clinical features of malaria. In this study, we analyzed the transcriptomes of isolates obtained from asymptomatic carriers and patients with uncomplicated or cerebral malaria. We also investigated the transcriptomes of 3D7 clone and 3D7-Lib that expresses severe malaria associated-variant surface antigen. Our findings revealed a specific up-regulation of genes involved in pathogenesis, adhesion to host cell, and erythrocyte aggregation in parasites from patients with cerebral malaria and 3D7-Lib, compared to parasites from asymptomatic carriers and 3D7, respectively. However, we did not find any significant difference between the transcriptomes of parasites from cerebral malaria and uncomplicated malaria, suggesting similar transcriptomic pattern in these two parasite populations. The difference between isolates from asymptomatic children and cerebral malaria concerned genes coding for exported proteins, Maurer's cleft proteins, transcriptional factor proteins, proteins implicated in protein transport, as well as Plasmodium conserved and hypothetical proteins. Interestingly, UPs A1, A2, A3 and UPs B1 of var genes were predominantly found in cerebral malaria-associated isolates and those containing architectural domains of DC4, DC5, DC13 and their neighboring rif genes in 3D7-lib. Therefore, more investigations are needed to analyze the effective role of these genes during malaria infection to provide with new knowledge on malaria pathology. In addition, concomitant regulation of genes within the chromosomal neighborhood suggests a common mechanism of gene regulation in P. falciparum.

Limited research has been conducted on the role of transcriptional regulators in relation to virulence in Leptospira interrogans, the etiological agent of leptospirosis. Here, we identify an L. interrogans locus that encodes a sensor protein, an anti-sigma factor antagonist, and two genes encoding proteins of unknown function. Transposon insertion into the gene encoding the sensor protein led to dampened transcription of the other 3 genes in this locus. This lb139 insertion mutant (the lb139(-) mutant) displayed attenuated virulence in the hamster model of infection and reduced motility in vitro. Whole-transcriptome analyses using RNA sequencing revealed the downregulation of 115 genes and the upregulation of 28 genes, with an overrepresentation of gene products functioning in motility and signal transduction and numerous gene products with unknown functions, predicted to be localized to the extracellular space. Another significant finding encompassed suppressed expression of the majority of the genes previously demonstrated to be upregulated at physiological osmolarity, including the sphingomyelinase C precursor Sph2 and LigB. We provide insight into a possible requirement for transcriptional regulation as it relates to leptospiral virulence and suggest various biological processes that are affected due to the loss of native expression of this genetic locus.

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.

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.

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., (© 2013 John Wiley & Sons Ltd.)

Archaeal viruses display unusually high genetic and morphological diversity. Studies of these viruses proved to be instrumental for the expansion of knowledge on viral diversity and evolution. The Sulfolobus islandicus rod-shaped virus 2 (SIRV2) is a model to study virus-host interactions in Archaea. It is a lytic virus that exploits a unique egress mechanism based on the formation of remarkable pyramidal structures on the host cell envelope. Using whole-transcriptome sequencing, we present here a global map defining host and viral gene expression during the infection cycle of SIRV2 in its hyperthermophilic host S. islandicus LAL14/1. This information was used, in combination with a yeast two-hybrid analysis of SIRV2 protein interactions, to advance current understanding of viral gene functions. As a consequence of SIRV2 infection, transcription of more than one-third of S. islandicus genes was differentially regulated. While expression of genes involved in cell division decreased, those genes playing a role in antiviral defense were activated on a large scale. Expression of genes belonging to toxin-antitoxin and clustered regularly interspaced short palindromic repeat (CRISPR)-Cas systems was specifically pronounced. The observed different degree of activation of various CRISPR-Cas systems highlights the specialized functions they perform. The information on individual gene expression and activation of antiviral defense systems is expected to aid future studies aimed at detailed understanding of the functions and interplay of these systems in vivo.

RNA viruses exhibit small-sized genomes encoding few proteins, but still establish complex networks of protein-protein and RNA-protein interactions within a cell to achieve efficient replication and spreading. Deciphering these interactions is essential to reach a comprehensive understanding of the viral infection process. To study RNA-protein complexes directly in infected cells, we developed a new approach based on recombinant viruses expressing tagged viral proteins that were purified together with their specific RNA partners. High-throughput sequencing was then used to identify these RNA molecules. As a proof of principle, this method was applied to measles virus nucleoprotein (MV-N). It revealed that in addition to full-length genomes, MV-N specifically interacted with a unique population of 5' copy-back defective interfering RNA genomes that we characterized. Such RNA molecules were able to induce strong activation of interferon-stimulated response element promoter preferentially via the cytoplasmic pattern recognition receptor RIG-I protein, demonstrating their biological functionality. Thus, this method provides a new platform to explore biologically active RNA-protein networks that viruses establish within infected cells.

Alternative splicing and polyadenylation were observed pervasively in eukaryotic messenger RNAs. These alternative isoforms could either be consequences of physiological regulation or stochastic noise of RNA processing. To quantify the extent of stochastic noise in splicing and polyadenylation, we analyzed the alternative usage of splicing and polyadenylation sites in Entamoeba histolytica using RNA-Seq. First, we identified a large number of rarely spliced alternative junctions and then showed that the occurrence of these alternative splicing events is correlated with splicing site sequence, occurrence of constitutive splicing events and messenger RNA abundance. Our results implied the majority of these alternative splicing events are likely to be stochastic error of splicing machineries, and we estimated the corresponding error rates. Second, we observed extensive microheterogeneity of polyadenylation cleavage sites, and the extent of such microheterogeneity is correlated with the occurrence of constitutive cleavage events, suggesting most of such microheterogeneity is likely to be stochastic. Overall, we only observed a small fraction of alternative splicing and polyadenylation isoforms that are unlikely to be solely stochastic, implying the functional relevance of alternative splicing and polyadenylation in E. histolytica is limited. Lastly, we revised the gene models and annotated their 3'UTR in AmoebaDB, providing valuable resources to the community.

Arsenic is widespread in the environment and its presence is a result of natural or anthropogenic activities. Microbes have developed different mechanisms to deal with toxic compounds such as arsenic and this is to resist or metabolize the compound. Here, we present the first reference set of genomic, transcriptomic and proteomic data of an Alphaproteobacterium isolated from an arsenic-containing goldmine: Rhizobium sp. NT-26. Although phylogenetically related to the plant-associated bacteria, this organism has lost the major colonizing capabilities needed for symbiosis with legumes. In contrast, the genome of Rhizobium sp. NT-26 comprises a megaplasmid containing the various genes, which enable it to metabolize arsenite. Remarkably, although the genes required for arsenite oxidation and flagellar motility/biofilm formation are carried by the megaplasmid and the chromosome, respectively, a coordinate regulation of these two mechanisms was observed. Taken together, these processes illustrate the impact environmental pressure can have on the evolution of bacterial genomes, improving the fitness of bacterial strains by the acquisition of novel functions.

Entamoeba histolytica is the pathogenic amoeba responsible for amoebiasis, an infectious disease targeting human tissues. Amoebiasis arises when virulent trophozoites start to destroy the muco-epithelial barrier by first crossing the mucus, then killing host cells, triggering inflammation and subsequently causing dysentery. The main goal of this study was to analyse pathophysiology and gene expression changes related to virulent (i.e. HM1:IMSS) and non-virulent (i.e. Rahman) strains when they are in contact with the human colon. Transcriptome comparisons between the two strains, both in culture conditions and upon contact with human colon explants, provide a global view of gene expression changes that might contribute to the observed phenotypic differences. The most remarkable feature of the virulent phenotype resides in the up-regulation of genes implicated in carbohydrate metabolism and processing of glycosylated residues. Consequently, inhibition of gene expression by RNA interference of a glycoside hydrolase (β-amylase absent from humans) abolishes mucus depletion and tissue invasion by HM1:IMSS. In summary, our data suggest a potential role of carbohydrate metabolism in colon invasion by virulent E. histolytica.

Unlabelled: We developed a modular and scalable framework called Eoulsan, based on the Hadoop implementation of the MapReduce algorithm dedicated to high-throughput sequencing data analysis. Eoulsan allows users to easily set up a cloud computing cluster and automate the analysis of several samples at once using various software solutions available. Our tests with Amazon Web Services demonstrated that the computation cost is linear with the number of instances booked as is the running time with the increasing amounts of data., Availability and Implementation: Eoulsan is implemented in Java, supported on Linux systems and distributed under the LGPL License at: http://transcriptome.ens.fr/eoulsan/

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.

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.

Although Yersinia pestis and Yersinia pseudotuberculosis are genetically very similar (97% nucleotide sequence identity for most of the chromosomal genes), they exhibit very different patterns of infection. Y. pestis causes plague which is usually fatal in the absence of treatment, whereas Y. pseudotuberculosis generally triggers non-life-threatening intestinal symptoms. This drastic difference in pathogenicity may result from the acquisition of a few species-specific genes, but also from differences in their transcriptional regulation networks. In this study, we performed an in silico comparative whole-genome transcriptome analysis of Y. pestis and Y. pseudotuberculosis grown in parallel under 8 distinct conditions to determine whether they exhibit differences in their regulatory networks. In this analysis, 304 genes common to both species were found to display significant inter-species differences in transcriptional levels, with 91% of them being more expressed in Y. pestis. Remarkably, 3 major virulence determinants conserved in the 2 species (the pYV virulence plasmid, the High Pathogenicity Island, and the ail locus) were among the genes more expressed in Y. pestis. Furthermore, the induction at 37°C of pYV-borne genes was considerably greater in Y. pestis than in Y. pseudotuberculosis. Conversely, the rovA transcriptional regulator gene was more transcribed in Y. pseudotuberculosis. We also performed a clustering analysis of the transcriptome data of both Y. pestis and Y. pseudotuberculosis, which allowed to group genes according to their expression profiles. This analysis identified groups of genes with unknown functions which, based on regulation patterns similar to those of known virulence genes, are potential new virulence determinants in Y. pestis. In conclusion, this is the first comparative analysis at the whole-genome level of the transcription profiles of Y. pestis and Y. pseudotuberculosis. Our results suggest that the higher pathogenicity of the plague bacillus may not only result from the acquisition of new genetic material, but also from a higher expression level of common crucial virulence genes. This in silico analysis thus opens new avenues for investigating Y. pestis gain of pathogenicity and new potential virulence factors., (Copyright © 2010 Elsevier GmbH. All rights reserved.)

Background: Arsenic 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)., Results: Microarrays 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., Conclusion: Gene 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.