PlaScope: a targeted approach to assess the plasmidome from genome assemblies at the species level.

Plasmid prediction may be of great interest when studying bacteria of medical importance such as Enterobacteriaceae as well as Staphylococcus aureus or Enterococcus. Indeed, many resistance and virulence genes are located on such replicons with major impact in terms of pathogenicity and spreading capacities. Beyond strain outbreak, plasmid outbreaks have been reported in particular for some extended-spectrum beta-lactamase- or carbapenemase-producing Enterobacteriaceae. Several tools are now available to explore the 'plasmidome' from whole-genome sequences with various approaches, but none of them are able to combine high sensitivity and specificity. With this in mind, we developed PlaScope, a targeted approach to recover plasmidic sequences in genome assemblies at the species or genus level. Based on Centrifuge, a metagenomic classifier, and a custom database containing complete sequences of chromosomes and plasmids from various curated databases, PlaScope classifies contigs from an assembly according to their predicted location. Compared to other plasmid classifiers, PlasFlow and cBar, it achieves better recall (0.87), specificity (0.99), precision (0.96) and accuracy (0.98) on a dataset of 70 genomes of Escherichia coli containing plasmids. In a second part, we identified 20 of the 21 chromosomal integrations of the extended-spectrum beta-lactamase coding gene in a clinical dataset of E. coli strains. In addition, we predicted virulence gene and operon locations in agreement with the literature. We also built a database for Klebsiella and correctly assigned the location for the majority of resistance genes from a collection of 12 Klebsiella pneumoniae strains. Similar approaches could also be developed for other well-characterized bacteria.

MiRImpact, a new bioinformatic method using complete microRNA expression profiles to assess their overall influence on the activity of intracellular molecular pathways.

MicroRNAs (miRs) are short noncoding RNA molecules that regulate expression of target mRNAs. Many published sources provide information about miRs and their targets. However, bioinformatic tools elucidating higher level impact of the established total miR profiles, are still largely missing. Recently, we developed a method termed OncoFinder enabling quantification of the activities of intracellular molecular pathways basing on gene expression data. Here we propose a new technique, MiRImpact, which enables to link miR expression data with its estimated outcome on the regulation of molecular pathways, like signaling, metabolic, cytoskeleton rearrangement, and DNA repair pathways. MiRImpact uses OncoFinder rationale for pathway activity calculations, with the major distinctions that (i) it deals with the concentrations of miRs--known regulators of gene products participating in molecular pathways, and (ii) miRs are considered as negative regulators of target molecules, if other is not specified. MiRImpact operates with 2 types of databases: for molecular targets of miRs and for gene products participating in molecular pathways. We applied MiRImpact to compare regulation of human bladder cancer-specific signaling pathways at the levels of mRNA and miR expression. We took 2 most complete alternative databases of experimentally validated miR targets--miRTarBase and DianaTarBase, and an OncoFinder database featuring 2725 gene products and 271 signaling pathways. We showed that the impact of miRs is orthogonal to pathway regulation at the mRNA level, which stresses the importance of studying posttranscriptional regulation of gene expression. We also report characteristic set of miR and mRNA regulation features linked with bladder cancer.