Fabian Amman

Dual RNA-seq unveils noncoding RNA functions in host–pathogen interactions

Bacteria express many small RNAs for which the regulatory roles in pathogenesis have remained poorly understood due to a paucity of robust phenotypes in standard virulence assays. Here we use a generic ‘dual RNA-seq’ approach to profile RNA expression simultaneously in pathogen and host during Salmonella enterica serovar Typhimurium infection and reveal the molecular impact of bacterial riboregulators. We identify a PhoP-activated small RNA, PinT, which upon bacterial internalization temporally controls the expression of both invasion-associated effectors and virulence genes required for intracellular survival. This riboregulatory activity causes pervasive changes in coding and noncoding transcripts of the host. Interspecies correlation analysis links PinT to host cell JAK–STAT signalling, and we identify infection-specific alterations in multiple long noncoding RNAs. Our study provides a paradigm for a sensitive RNA-based analysis of intracellular bacterial pathogens and their hosts without physical separation, as well as a new discovery route for hidden functions of pathogen genes.

Archaeal Signal Transduction: Impact of Protein Phosphatase Deletions on Cell Size, Motility, and Energy Metabolism in Sulfolobus acidocaldarius

In this study, the in vitro and in vivo functions of the only two identified protein phosphatases, Saci-PTP and Saci-PP2A, in the crenarchaeal model organism Sulfolobus acidocaldarius were investigated. Biochemical characterization revealed that Saci-PTP is a dual-specific phosphatase (against pSer/pThr and pTyr), whereas Saci-PP2A exhibited specific pSer/pThr activity and inhibition by okadaic acid. Deletion of saci_pp2a resulted in pronounced alterations in growth, cell shape and cell size, which could be partially complemented. Transcriptome analysis of the three strains (Δsaci_ptp, Δsaci_pp2a and the MW001 parental strain) revealed 155 genes that were differentially expressed in the deletion mutants, and showed significant changes in expression of genes encoding the archaella (archaeal motility structure), components of the respiratory chain and transcriptional regulators. Phosphoproteome studies revealed 801 unique phosphoproteins in total, with an increase in identified phosphopeptides in the deletion mutants. Proteins from most functional categories were affected by phosphorylation, including components of the motility system, the respiratory chain, and regulatory proteins. In the saci_pp2a deletion mutant the up-regulation at the transcript level, as well as the observed phosphorylation pattern, resembled starvation stress responses. Hypermotility was also observed in the saci_pp2a deletion mutant. The results highlight the importance of protein phosphorylation in regulating essential cellular processes in the crenarchaeon S. acidocaldarius.

TSSAR: TSS annotation regime for dRNA-seq data

BackgroundDifferential RNA sequencing (dRNA-seq) is a high-throughput screening technique designed to examine the architecture of bacterial operons in general and the precise position of transcription start sites (TSS) in particular. Hitherto, dRNA-seq data were analyzed by visualizing the sequencing reads mapped to the reference genome and manually annotating reliable positions. This is very labor intensive and, due to the subjectivity, biased.ResultsHere, we present TSSAR, a tool for automated de novo TSS annotation from dRNA-seq data that respects the statistics of dRNA-seq libraries. TSSAR uses the premise that the number of sequencing reads starting at a certain genomic position within a transcriptional active region follows a Poisson distribution with a parameter that depends on the local strength of expression. The differences of two dRNA-seq library counts thus follow a Skellam distribution. This provides a statistical basis to identify significantly enriched primary transcripts.We assessed the performance by analyzing a publicly available dRNA-seq data set using TSSAR and two simple approaches that utilize user-defined score cutoffs. We evaluated the power of reproducing the manual TSS annotation. Furthermore, the same data set was used to reproduce 74 experimentally validated TSS in H. pylori from reliable techniques such as RACE or primer extension. Both analyses showed that TSSAR outperforms the static cutoff-dependent approaches.ConclusionsHaving an automated and efficient tool for analyzing dRNA-seq data facilitates the use of the dRNA-seq technique and promotes its application to more sophisticated analysis. For instance, monitoring the plasticity and dynamics of the transcriptomal architecture triggered by different stimuli and growth conditions becomes possible.The main asset of a novel tool for dRNA-seq analysis that reaches out to a broad user community is usability. As such, we provide TSSAR both as intuitive RESTful Web service (http://rna.tbi.univie.ac.at/TSSAR) together with a set of post-processing and analysis tools, as well as a stand-alone version for use in high-throughput dRNA-seq data analysis pipelines.

Fast accessibility-based prediction of RNA–RNA interactions

MOTIVATION Currently, the best RNA-RNA interaction prediction tools are based on approaches that consider both the inter- and intramolecular interactions of hybridizing RNAs. While accurate, these methods are too slow and memory-hungry to be employed in genome-wide RNA target scans. Alternative methods neglecting intramolecular structures are fast enough for genome-wide applications, but are too inaccurate to be of much practical use. RESULTS A new approach for RNA-RNA interaction was developed, with a prediction accuracy that is similar to that of algorithms that explicitly consider intramolecular structures, but running at least three orders of magnitude faster than RNAup. This is achieved by using a combination of precomputed accessibility profiles with an approximate energy model. This approach is implemented in the new version of RNAplex. The software also provides a variant using multiple sequences alignments as input, resulting in a further increase in specificity. AVAILABILITY RNAplex is available at www.bioinf.uni-leipzig.de/Software/RNAplex.

2D Meets 4G: G-Quadruplexes in RNA Secondary Structure Prediction

G-quadruplexes are abundant locally stable structural elements in nucleic acids. The combinatorial theory of RNA structures and the dynamic programming algorithms for RNA secondary structure prediction are extended here to incorporate G-quadruplexes using a simple but plausible energy model. With preliminary energy parameters, we find that the overwhelming majority of putative quadruplex-forming sequences in the human genome are likely to fold into canonical secondary structures instead. Stable G-quadruplexes are strongly enriched, however, in the 5Ê1UTR of protein coding mRNAs.

Animal snoRNAs and scaRNAs with exceptional structures

The overwhelming majority of small nucleolar RNAs (snoRNAs) fall into two clearly defined classes characterized by distinctive secondary structures and sequence motifs. A small group of diverse ncRNAs, however, shares the hallmarks of one or both classes of snoRNAs but differs substantially from the norm in some respects. Here, we compile the available information on these exceptional cases, conduct a thorough homology search throughout the available metazoan genomes, provide improved and expanded alignments, and investigate the evolutionary histories of these ncRNA families as well as their mutual relationships.

Impact of Hfq on the Bacillus subtilis Transcriptome

The RNA chaperone Hfq acts as a central player in post-transcriptional gene regulation in several Gram-negative Bacteria, whereas comparatively little is known about its role in Gram-positive Bacteria. Here, we studied the function of Hfq in Bacillus subtilis, and show that it confers a survival advantage. A comparative transcriptome analysis revealed mRNAs with a differential abundance that are governed by the ResD-ResE system required for aerobic and anaerobic respiration. Expression of resD was found to be up-regulated in the hfq− strain. Furthermore, several genes of the GerE and ComK regulons were de-regulated in the hfq− background. Surprisingly, only six out of >100 known and predicted small RNAs (sRNAs) showed altered abundance in the absence of Hfq. Moreover, Hfq positively affected the transcript abundance of genes encoding type I toxin-antitoxin systems. Taken the moderate effect on sRNA levels and mRNAs together, it seems rather unlikely that Hfq plays a central role in RNA transactions in Bacillus subtilis.

ViennaNGS: A toolbox for building efficient next- generation sequencing analysis pipelines

Recent achievements in next-generation sequencing (NGS) technologies lead to a high demand for reuseable software components to easily compile customized analysis workflows for big genomics data. We present ViennaNGS, an integrated collection of Perl modules focused on building efficient pipelines for NGS data processing. It comes with functionality for extracting and converting features from common NGS file formats, computation and evaluation of read mapping statistics, as well as normalization of RNA abundance. Moreover, ViennaNGS provides software components for identification and characterization of splice junctions from RNA-seq data, parsing and condensing sequence motif data, automated construction of Assembly and Track Hubs for the UCSC genome browser, as well as wrapper routines for a set of commonly used NGS command line tools.

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.

Transcriptional profiling of Bordetella pertussis reveals requirement of RNA chaperone Hfq for Type III secretion system functionality

Bordetella pertussis, the causative agent of human whooping cough (pertussis) produces a complex array of virulence factors in order to establish efficient infection in the host. The RNA chaperone Hfq and small regulatory RNAs are key players in posttranscriptional regulation in bacteria and have been shown to play an essential role in virulence of a broad spectrum of bacterial pathogens. This study represents the first attempt to characterize the Hfq regulon of the human pathogen B. pertussis under laboratory conditions as well as upon passage in the host and indicates that loss of Hfq has a profound effect on gene expression in B. pertussis. Comparative transcriptional profiling revealed that Hfq is required for expression of several virulence factors in B. pertussis cells including the Type III secretion system (T3SS). In striking contrast to the wt strain, T3SS did not become operational in the hfq mutant passaged either through mice or macrophages thereby proving that Hfq is required for the functionality of the B. pertussis T3SS. Likewise, expression of virulence factors vag8 and tcfA encoding autotransporter and tracheal colonization factor, respectively, was strongly reduced in the hfq mutant. Importantly, for the first time we demonstrate that B. pertussis T3SS can be activated upon contact with macrophage cells in vitro.