MSF: Modulated Sub-graph Finder [version 3; peer review: 1 approved, 2 approved with reservations]

High throughput techniques such as RNA-seq or microarray analysis have proven to be invaluable for the characterizing of global transcriptional gene activity changes due to external stimuli or diseases. Differential gene expression analysis (DGEA) is the first step in the course of data interpretation, typically producing lists of dozens to thousands of differentially expressed genes. To further guide the interpretation of these lists, different pathway analysis approaches have been developed. These tools typically rely on the classification of genes into sets of genes, such as pathways, based on the interactions between the genes and their function in a common biological process. Regardless of technical differences, these methods do not properly account for cross talk between different pathways and most of the methods rely on binary separation into differentially expressed gene and unaffected genes based on an arbitrarily set p-value cut-off. To overcome this limitation, we developed a novel approach to identify concertedly modulated sub-graphs in the global cell signaling network, based on the DGEA results of all genes tested. To this end, expression patterns of genes are integrated according to the topology of their interactions and allow potentially to read the flow of information and identify the effectors. The described software, named Modulated Sub-graph Finder (MSF) is freely available at https://github.com/Modulated-Subgraph-Finder/MSF

RNASeq Based Transcriptional Profiling of Pseudomonas aeruginosa PA14 after Short- and Long-Term Anoxic Cultivation in Synthetic Cystic Fibrosis Sputum Medium

The opportunistic human pathogen Pseudomonas aeruginosa can thrive under microaerophilic to anaerobic conditions in the lungs of cystic fibrosis patients. RNA$_{Seq}$ based comparative RNA profiling of the clinical isolate PA14 cultured in synthetic cystic fibrosis medium was performed after planktonic growth (OD$_{600}$ = 2.0; P), 30 min after shift to anaerobiosis (A-30) and after anaerobic biofilm growth for 96h (B-96) with the aim to reveal differentially regulated functions impacting on sustained anoxic biofilm formation as well as on tolerance towards different antibiotics. Most notably, functions involved in sulfur metabolism were found to be up-regulated in B-96 cells when compared to A-30 cells. Based on the transcriptome studies a set of transposon mutants were screened, which revealed novel functions involved in anoxic biofilm growth.In addition, these studies revealed a decreased and an increased abundance of the oprD and the mexCD-oprJ operon transcripts, respectively, in B-96 cells, which may explain their increased tolerance towards meropenem and to antibiotics that are expelled by the MexCD-OprD efflux pump. The OprI protein has been implicated as a target for cationic antimicrobial peptides, such as SMAP-29. The transcriptome and subsequent Northern-blot analyses showed that the abundance of the oprI transcript encoding the OprI protein is strongly decreased in B-96 cells. However, follow up studies revealed that the susceptibility of a constructed PA14ΔoprI mutant towards SMAP-29 was indistinguishable from the parental wild-type strain, which questions OprI as a target for this antimicrobial peptide in strain PA14

Comparative genomics of Czech vaccine strains of Bordetella pertussis.

Bordetella pertussis is a strictly human pathogen causing the respiratory infectious disease called whooping cough or pertussis. B. pertussis adaptation to acellular pertussis vaccine pressure has been repeatedly highlighted, but recent data indicate that adaptation of circulating strains started already in the era of the whole cell pertussis vaccine (wP) use. We sequenced the genomes of five B. pertussis wP vaccine strains isolated in the former Czechoslovakia in the pre-wP (1954-1957) and early wP (1958-1965) eras, when only limited population travel into and out of the country was possible. Four isolates exhibit a similar genome organization and form a distinct phylogenetic cluster with a geographic signature. The fifth strain is rather distinct, both in genome organization and SNP-based phylogeny. Surprisingly, despite isolation of this strain before 1966, its closest sequenced relative appears to be a recent isolate from the US. On the genome content level, the five vaccine strains contained both new and already described regions of difference. One of the new regions contains duplicated genes potentially associated with transport across the membrane. The prevalence of this region in recent isolates indicates that its spread might be associated with selective advantage leading to increased strain fitness

Negative Control of RpoS Synthesis by the sRNA ReaL in Pseudomonas aeruginosa

Pseudomonas aeruginosa (Pae) is an opportunistic human pathogen, able to resist host defense mechanisms and antibiotic treatment. In Pae, the master regulator of stress responses RpoS (σS) is involved in the regulation of quorum sensing and several virulence genes. Here, we report that the sRNA ReaL translationally silences rpoS mRNA, which results in a decrease of the RpoS levels. Our studies indicated that ReaL base-pairs with the Shine-Dalgarno region of rpoS mRNA. These studies are underlined by a highly similar transcription profile of a rpoS deletion mutant and a reaL over-expressing strain

Differential transcriptional responses to Ebola and Marburg virus infection in bat and human cells

The unprecedented outbreak of Ebola in West Africa resulted in over 28,000 cases and 11,000 deaths, underlining the need for a better understanding of the biology of this highly pathogenic virus to develop specific counter strategies. Two filoviruses, the Ebola and Marburg viruses, result in a severe and often fatal infection in humans. However, bats are natural hosts and survive filovirus infections without obvious symptoms. The molecular basis of this striking difference in the response to filovirus infections is not well understood. We report a systematic overview of differentially expressed genes, activity motifs and pathways in human and bat cells infected with the Ebola and Marburg viruses, and we demonstrate that the replication of filoviruses is more rapid in human cells than in bat cells. We also found that the most strongly regulated genes upon filovirus infection are chemokine ligands and transcription factors. We observed a strong induction of the JAK/STAT pathway, of several genes encoding inhibitors of MAP kinases (DUSP genes) and of PPP1R15A, which is involved in ER stress-induced cell death. We used comparative transcriptomics to provide a data resource that can be used to identify cellular responses that might allow bats to survive filovirus infections.Additional co-authors: Andreas J. Gruber, Franziska Hufsky, Henrike Indrischek, Sabina Kanton, Jörg Linde, Nelly Mostajo, Roman Ochsenreiter, Konstantin Riege, Lorena Rivarola-Duarte, Abdullah H. Sahyoun, Sita J. Saunders, Stefan E. Seemann, Andrea Tanzer, Bertram Vogel, Michael T. Wolfinger, Rolf Backofen, Jan Gorodkin, Ivo Grosse, Ivo Hofacker, Steve Hoffmann, Christoph Kaleta, Peter F. Stadler, Stephan Becker, and Manja Marz

Transcriptional profiling of human macrophages during infection with <i>Bordetella pertussis</i>

Bordetella pertussis, a strictly human re-emerging pathogen and the causative agent of whooping cough, exploits a broad variety of virulence factors to establish efficient infection. Here, we used RNA sequencing to analyse the changes in gene expression profiles of human THP-1 macrophages resulting from B. pertussis infection. In parallel, we attempted to determine the changes in intracellular B. pertussis-specific transcriptomic profiles resulting from interaction with macrophages. Our analysis revealed that global gene expression profiles in THP-1 macrophages are extensively rewired 6 h post-infection. Among the highly expressed genes, we identified those encoding cytokines, chemokines, and transcription regulators involved in the induction of the M1 and M2 macrophage polarization programmes. Notably, several host genes involved in the control of apoptosis and inflammation which are known to be hijacked by intracellular bacterial pathogens were overexpressed upon infection. Furthermore, in silico analyses identified large temporal changes in expression of specific gene subsets involved in signalling and metabolic pathways. Despite limited numbers of the bacterial reads, we observed reduced expression of majority of virulence factors and upregulation of several transcriptional regulators during infection suggesting that intracellular B. pertussis cells switch from virulent to avirulent phase and actively adapt to intracellular environment, respectively.Facultad de Ciencias ExactasCentro de Investigación y Desarrollo en Fermentaciones Industriale

Modelling Translation Initiation under the Influence of sRNA

Bacterial small non-coding RNA (sRNA) plays an important role in post-transcriptional gene regulation. Although the number of annotated sRNA is steadily increasing, their functional characterization is still lagging behind. Various computational strategies for finding sRNA–mRNA interactions, and thus putative sRNA targets, were developed. Most of them suffer from a high false positive rate. Here, we present a qualitative model to simulate the effect of an sRNA on the translation initiation of a potential target. Information about the ribosome–mRNA interaction, sRNA–mRNA interaction and expression information from deep sequencing experiments is integrated to calculate the change in translation initiation complex formation, as a proxy for translational activity. This model can be used to post-evaluate predicted targets, hence condensing the list of potential targets. We show that our translation initiation model, under the influence of an sRNA, can successfully simulate thirteen out of fifteen tested sRNA–mRNA interactions in a qualitative manner. To show the gain in specificity, we applied our method to a target search for the Escherichia coli sRNA RyhB. Compared with simple target prediction without post-evaluation, we reduce the number of targets to less than one fourth potential targets, considerably reducing the burden of experimental validation