ALBERT

Description: Motivation: Random sampling of the solution space has emerged as a popular tool to explore and infer properties of large metabolic networks. However, conventional sampling approaches commonly used do not eliminate thermodynamically unfeasible loops. Results: In order to overcome this limitation, we developed an efficient sampling algorithm called loopless Artificially Centered Hit-and-Run on a Box (ll-ACHRB). This algorithm is inspired by the Hit-and-Run on a Box algorithm for uniform sampling from general regions, but employs the directions of choice approach of Artificially Centered Hit-and-Run. A novel strategy for generating feasible warmup points improved both sampling efficiency and mixing. ll-ACHRB shows overall better performance than current strategies to generate feasible flux samples across several models. Furthermore, we demonstrate that a failure to eliminate unfeasible loops greatly affects sample statistics, in particular the correlation structure. Finally, we discuss recommendations for the interpretation of sampling results and possible algorithmic improvements. Availability and implementation: Source code for MATLAB and OCTAVE including examples are freely available for download at http://www.aibn.uq.edu.au/cssb-resources under Software. Optimization runs can use Gurobi Optimizer (by default if available) or GLPK (included with the algorithm). Contact: lars.nielsen@uq.edu.au Supplementary information: Supplementary data are available at Bioinformatics online.

Description: The Gram-negative proteobacteria genus Burkholderia encompasses multiple bacterial species that are pathogenic to humans and other vertebrates. Two pathogenic species of interest within this genus are Burkholderia pseudomallei (Bpm) and the B. cepacia complex (Bcc); the former is the causative agent of melioidosis in humans and other mammals, and the latter is associated with pneumonia in immunocompromised patients. One understudied and shared characteristic of these two pathogenic groups is their ability to persist and establish chronic infection within the host. In this review, we will explore the depth of knowledge about chronic infections caused by persistent Bpm and Bcc. We examine the host risk factors and immune responses associated with more severe chronic infections. We also discuss host adaptation and phenotypes associated with persistent Burkholderia species. Lastly, we survey how other intracellular bacteria associated with chronic infections are combatted and explore possible future applications to target Burkholderia . Our goal is to highlight understudied areas that should be addressed for a more thorough understanding of chronic Burkholderia infections and how to combat them.

Description: Macrophages are phagocytic cells that play a key role maintaining the homeostasis of many tissues. Their function is essential for controlling and eradicating infecting mycobacteria. Human monocytic cell lines such as THP-1 and U937 have provided interesting insights into how mycobacteria subvert the host cell response. However, immortalized cell lines could bring some disadvantages. Here we compare the response of THP-1 and U937 cell lines with human monocyte-derived macrophages (hMDMs) to determine functional differences during infection with different mycobacterial phenotypes (virulent Mycobacterium tuberculosis H37Rv and Mycobacterium bovis , and attenuated M. bovis BCG). The findings of this study indicate that the U937 cell line displays a significantly lower phagocytic capacity than hMDMs and THP-1 macrophages, regardless of the mycobacterial strain. In all cell models, interferon- activation leads to up-regulation of interleukin-12 and nitrite production. However, the phorbol 12-myristate 13-acetate (PMA)-induced differentiation of U937 and THP-1 cell lines induces a significant tumor necrosis factor-α production in resting macrophages. However, this state of activation has no effect on the control of intracellular growth of mycobacteria. Moreover, U937 cells show more discrepancies with hMDM than THP-1. This study demonstrates that THP-1 macrophages exhibit closer functional similarities to hMDMs in response to mycobacterial infection, regardless of the strain.

Description: Urinary tract infections (UTIs) caused by uropathogenic Escherichia coli (UPEC) are common infections in humans. Despite the substantial healthcare cost represented by these infections, the human immune response associated with the infection immediately following the onset of symptoms in patients remains largely undefined. We performed a prospective study aimed at defining the milieu of urinary cytokines in adult inpatients in the 24–48 h period immediately following hospital admission for acute cystitis due to UPEC. Urine samples, analyzed using 27-target multiplex protein assays, were used to generate immune profiles for patients and compared to age- and gender-matched healthy controls. The levels of multiple pro-inflammatory cytokines were significantly elevated in urine as a result of infection, an observation consistent with prior findings in murine models and clinical literature. We also identified significant responses for several novel factors not previously associated with the human response to UTI, including Interleukin (IL)-4, IL-7, IL-9, IL-17A, eotaxin, Granulocyte-macrophage colony-stimulating factor (GM-CSF) and several growth factors. These data establish crucial parallels between the human immune response to UPEC and murine model UTI studies, and emphasize the complex but poorly defined nature of the human immune response to UPEC, particularly in the immediate period following the onset of symptoms for acute cystitis.

Description: Staphylococcus aureus is a ubiquitous, versatile and dangerous pathogen. It colonizes over 30% of the human population, and is one of the leading causes of death by an infectious agent. During S. aureus colonization and invasion, leukocytes are recruited to the site of infection. To combat S. aureus , leukocytes generate an arsenal of reactive species including superoxide, hydrogen peroxide, nitric oxide and hypohalous acids that modify and inactivate cellular macromolecules, resulting in growth defects or death. When S. aureus colonization cannot be cleared by the immune system, antibiotic treatment is necessary and can be effective. Yet, this organism quickly gains resistance to each new antibiotic it encounters. Therefore, it is in the interest of human health to acquire a deeper understanding of how S. aureus evades killing by the immune system. Advances in this field will have implications for the design of future S. aureus treatments that complement and assist the host immune response. In that regard, this review focuses on how S. aureus avoids host-generated oxidative stress, and discusses the mechanisms used by S. aureus to survive oxidative damage including antioxidants, direct repair of damaged proteins, sensing oxidant stress and transcriptional changes. This review will elucidate areas for studies to identify and validate future antimicrobial targets.

Description: Celastrol is a bioactive compound derived from traditional Chinese medicinal herbs of the Celastraceae family. Celastrol is known to possess anti-inflammatory and anti-oxidant activities. Our studies have highlighted the immunomodulatory attributes of celastrol in adjuvant-induced arthritis (AA), an experimental model of human rheumatoid arthritis (RA). RA is an autoimmune disease characterized by chronic inflammation of the synovial lining of the joints, leading eventually to tissue damage and deformities. Identification of the molecular targets of celastrol such as the NF-B pathway, MAPK pathway, JAK/STAT pathway and RANKL/OPG pathway has unraveled its strategic checkpoints in controlling arthritic inflammation and tissue damage in AA. The pathological events that are targeted and rectified by celastrol include increased production of pro-inflammatory cytokines; an imbalance between pathogenic T helper 17 and regulatory T cells; enhanced production of chemokines coupled with increased migration of immune cells into the joints; and increased release of mediators of osteoclastic bone damage. Accordingly, celastrol is a promising candidate for further testing in the clinic for RA therapy. Furthermore, the results of other preclinical studies suggest that celastrol might also be beneficial for the treatment of a few other autoimmune diseases besides arthritis.

Description: Results: Here, we present a comprehensive analysis on the reproducibility of computational characterization of genomic variants using high throughput sequencing data. We reanalyzed the same datasets twice, using the same tools with the same parameters, where we only altered the order of reads in the input (i.e. FASTQ file). Reshuffling caused the reads from repetitive regions being mapped to different locations in the second alignment, and we observed similar results when we only applied a scatter/gather approach for read mapping—without prior shuffling. Our results show that, some of the most common variation discovery algorithms do not handle the ambiguous read mappings accurately when random locations are selected. In addition, we also observed that even when the exact same alignment is used, the GATK HaplotypeCaller generates slightly different call sets, which we pinpoint to the variant filtration step. We conclude that, algorithms at each step of genomic variation discovery and characterization need to treat ambiguous mappings in a deterministic fashion to ensure full replication of results. Availability and Implementation: Code, scripts and the generated VCF files are available at DOI:10.5281/zenodo.32611. Contact: calkan@cs.bilkent.edu.tr Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Animals from worms and insects to birds and mammals show distinct body plans; however, the embryonic development of diverse body plans with tissues and organs within is controlled by a surprisingly few signaling pathways. It is well recognized that combinatorial use of and dynamic interactions among signaling pathways follow specific logic to control complex and accurate developmental signaling and patterning, but it remains elusive what such logic is, or even, what it looks like. Results: We have developed a computational model for Drosophila eye development with innovated methods to reveal how interactions among multiple pathways control the dynamically generated hexagonal array of R8 cells. We obtained two novel findings. First, the coupling between the long-range inductive signals produced by the proneural Hh signaling and the short-range restrictive signals produced by the antineural Notch and EGFR signaling is essential for generating accurately spaced R8s. Second, the spatiotemporal orders of key signaling events reveal a robust pattern of lateral inhibition conducted by Ato-coordinated Notch and EGFR signaling to collectively determine R8 patterning. This pattern, stipulating the orders of signaling and comparable to the protocols of communication, may help decipher the well-appreciated but poorly defined logic of developmental signaling. Availability and implementation: The model is available upon request. Contact: hao.zhu@ymail.com Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: The growing amount of regulatory data from the ENCODE, Roadmap Epigenomics and other consortia provides a wealth of opportunities to investigate the functional impact of single nucleotide polymorphisms (SNPs). Yet, given the large number of regulatory datasets, researchers are posed with a challenge of how to efficiently utilize them to interpret the functional impact of SNP sets. Results: We developed the GenomeRunner web server to automate systematic statistical analysis of SNP sets within a regulatory context. Besides defining the functional impact of SNP sets, GenomeRunner implements novel regulatory similarity/differential analyses, and cell type-specific regulatory enrichment analysis. Validated against literature- and disease ontology-based approaches, analysis of 39 disease/trait-associated SNP sets demonstrated that the functional impact of SNP sets corresponds to known disease relationships. We identified a group of autoimmune diseases with SNPs distinctly enriched in the enhancers of T helper cell subpopulations, and demonstrated relevant cell type-specificity of the functional impact of other SNP sets. In summary, we show how systematic analysis of genomic data within a regulatory context can help interpreting the functional impact of SNP sets. Availability and Implementation: GenomeRunner web server is freely available at http://www.integrativegenomics.org/ . Contact: mikhail.dozmorov@gmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Adverse drug reactions (ADRs) are a central consideration during drug development. Here we present a machine learning classifier to prioritize ADRs for approved drugs and pre-clinical small-molecule compounds by combining chemical structure (CS) and gene expression (GE) features. The GE data is from the Library of Integrated Network-based Cellular Signatures (LINCS) L1000 dataset that measured changes in GE before and after treatment of human cells with over 20 000 small-molecule compounds including most of the FDA-approved drugs. Using various benchmarking methods, we show that the integration of GE data with the CS of the drugs can significantly improve the predictability of ADRs. Moreover, transforming GE features to enrichment vectors of biological terms further improves the predictive capability of the classifiers. The most predictive biological-term features can assist in understanding the drug mechanisms of action. Finally, we applied the classifier to all 〉20 000 small-molecules profiled, and developed a web portal for browsing and searching predictive small-molecule/ADR connections. Availability and Implementation: The interface for the adverse event predictions for the 〉20 000 LINCS compounds is available at http://maayanlab.net/SEP-L1000/ . Contact: avi.maayan@mssm.edu Supplementary information : Supplementary data are available at Bioinformatics online.

Description: Motivation: Environmental dissemination of antibiotic resistance genes (ARGs) has become an increasing concern for public health. Metagenomics approaches can effectively detect broad profiles of ARGs in environmental samples; however, the detection and subsequent classification of ARG-like sequences are time consuming and have been severe obstacles in employing metagenomic methods. We sought to accelerate quantification of ARGs in metagenomic data from environmental samples. Results: A Structured ARG reference database (SARG) was constructed by integrating ARDB and CARD, the two most commonly used databases. SARG was curated to remove redundant sequences and optimized to facilitate query sequence identification by similarity. A database with a hierarchical structure (type-subtype-reference sequence) was then constructed to facilitate classification (assigning ARG-like sequence to type, subtype and reference sequence) of sequences identified through similarity search. Utilizing SARG and a previously proposed hybrid functional gene annotation pipeline, we developed an online pipeline called ARGs-OAP for fast annotation and classification of ARG-like sequences from metagenomic data. We also evaluated and proposed a set of criteria important for efficiently conducting metagenomic analysis of ARGs using ARGs-OAP. Availability and Implementation: Perl script for ARGs-OAP can be downloaded from https://github.com/biofuture/Ublastx_stageone . ARGs-OAP can be accessed through http://smile.hku.hk/SARGs . Contact: zhangt@hku.hk or tiedjej@msu.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : Visualizing genomic data in chromosomal context can help detecting errors in data processing and may suggest new hypotheses to be tested. Here, we report a new tool for displaying large and diverse genomic data along chromosomes. The software is implemented in R so that visualization can be easily integrated with its numerous packages for processing genomic data. It supports simultaneous visualization of multiple tracks of data. Large genomic regions such as QTLs or synteny tracts may be shown along histograms of number of genes, genetic variants, or any other type of genomic element. Tracks can also contain values for continuous or categorical variables and the user can choose among points, connected lines, colored segments, or histograms for representing data. chromPlot takes data from tables in data.frame in GRanges formats. The information necessary to draw chromosomes for mouse and human is included with the package. For other organisms, chromPlot can read Gap and cytoBandIdeo tables from the UCSC Genome Browser. We present common use cases here, and a full tutorial is included as the package’s vignette. Availability and Implementation: chromPlot is distributed under a GLP2 licence at http://www.bioconductor.org . Contact: raverdugo@u.uchile.cl Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : The most important features of error correction tools for sequencing data are accuracy, memory efficiency and fast runtime. The previous version of BLESS was highly memory-efficient and accurate, but it was too slow to handle reads from large genomes. We have developed a new version of BLESS to improve runtime and accuracy while maintaining a small memory usage. The new version, called BLESS 2, has an error correction algorithm that is more accurate than BLESS, and the algorithm has been parallelized using hybrid MPI and OpenMP programming. BLESS 2 was compared with five top-performing tools, and it was found to be the fastest when it was executed on two computing nodes using MPI, with each node containing twelve cores. Also, BLESS 2 showed at least 11% higher gain while retaining the memory efficiency of the previous version for large genomes. Availability and implementation: Freely available at https://sourceforge.net/projects/bless-ec Contact: dchen@illinois.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Accurate segmentation of brain electron microscopy (EM) images is a critical step in dense circuit reconstruction. Although deep neural networks (DNNs) have been widely used in a number of applications in computer vision, most of these models that proved to be effective on image classification tasks cannot be applied directly to EM image segmentation, due to the different objectives of these tasks. As a result, it is desirable to develop an optimized architecture that uses the full power of DNNs and tailored specifically for EM image segmentation. Results: In this work, we proposed a novel design of DNNs for this task. We trained a pixel classifier that operates on raw pixel intensities with no preprocessing to generate probability values for each pixel being a membrane or not. Although the use of neural networks in image segmentation is not completely new, we developed novel insights and model architectures that allow us to achieve superior performance on EM image segmentation tasks. Our submission based on these insights to the 2D EM Image Segmentation Challenge achieved the best performance consistently across all the three evaluation metrics. This challenge is still ongoing and the results in this paper are as of June 5, 2015. Availability and Implementation : https://github.com/ahmed-fakhry/dive Contact : sji@eecs.wsu.edu

Description: : Hilbert curves enable high-resolution visualization of genomic data on a chromosome- or genome-wide scale. Here we present the HilbertCurve package that provides an easy-to-use interface for mapping genomic data to Hilbert curves. The package transforms the curve as a virtual axis, thereby hiding the details of the curve construction from the user. HilbertCurve supports multiple-layer overlay that makes it a powerful tool to correlate the spatial distribution of multiple feature types. Availability and implementation: The HilbertCurve package and documentation are freely available from the Bioconductor project: http://www.bioconductor.org/packages/devel/bioc/html/HilbertCurve.html Contact: m.schlesner@dkfz.de Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: The sequences among subgenomes in a polyploid species have high similarity, making it difficult to design genome-specific primers for sequence analysis. Results: We present GSP, a web-based platform to design genome-specific primers that distinguish subgenome sequences in a polyploid genome. GSP uses BLAST to extract homeologous sequences of the subgenomes in existing databases, performs a multiple sequence alignment, and design primers based on sequence variants in the alignment. An interactive primers diagram, a sequence alignment viewer and a virtual electrophoresis are displayed as parts of the primer design result. GSP also designs specific primers from multiple sequences uploaded by users. Availability and implementation: GSP is a user-friendly and efficient web platform freely accessible at http://probes.pw.usda.gov/GSP . Source code and command-line application are available at https://github.com/bioinfogenome/GSP . Contacts: yong.gu@ars.usda.gov or devin.coleman-derr@ars.usda.gov Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : The prediction of protein–protein complexes from the structures of unbound components is a challenging and powerful strategy to decipher the mechanism of many essential biological processes. We present a user-friendly protein–protein docking server based on an improved version of FRODOCK that includes a complementary knowledge-based potential. The web interface provides a very effective tool to explore and select protein–protein models and interactively screen them against experimental distance constraints. The competitive success rates and efficiency achieved allow the retrieval of reliable potential protein–protein binding conformations that can be further refined with more computationally demanding strategies. Availability and Implementation: The server is free and open to all users with no login requirement at http://frodock.chaconlab.org Contact: pablo@chaconlab.org Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : XIBD performs pairwise relatedness mapping on the X chromosome using dense single nucleotide polymorphism (SNP) data from either SNP chips or next generation sequencing data. It correctly accounts for the difference in chromosomal numbers between males and females and estimates global relatedness as well as regions of the genome that are identical by descent (IBD). XIBD also generates novel graphical summaries of all pairwise IBD tracts for a cohort making it very useful for disease locus mapping. Availability and implementation: XIBD is written in R/Rcpp and executed from shell scripts that are freely available from http://bioinf.wehi.edu.au/software/XIBD along with accompanying reference datasets. Contact: henden.l@wehi.edu.au Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : Several approaches to the region-based association analysis of quantitative traits have recently been developed and successively applied. However, no software package has been developed that implements all of these approaches for either independent or structured samples. Here we introduce FREGAT (Family REGional Association Tests), an R package that can handle family and population samples and implements a wide range of region-based association methods including burden tests, functional linear models, and kernel machine-based regression. FREGAT can be used in genome/exome-wide region-based association studies of quantitative traits and candidate gene analysis. FREGAT offers many useful options to empower its users and increase the effectiveness and applicability of region-based association analysis. Availability and Implementation: https://cran.r-project.org/web/packages/FREGAT/index.html Supplementary Information: Supplementary data are available at Bioinformatics Online. Contact: belon@bionet.nsc.ru

Description: Motivation: There is a growing need in bioinformatics for easy-to-use software implementations of algorithms that are usable across platforms. At the same time, reproducibility of computational results is critical and often a challenge due to source code changes over time and dependencies. Results: The approach introduced in this paper addresses both of these needs with AlgoRun, a dedicated packaging system for implemented algorithms, using Docker technology. Implemented algorithms, packaged with AlgoRun, can be executed through a user-friendly interface directly from a web browser or via a standardized RESTful web API to allow easy integration into more complex workflows. The packaged algorithm includes the entire software execution environment, thereby eliminating the common problem of software dependencies and the irreproducibility of computations over time. AlgoRun-packaged algorithms can be published on http://algorun.org , a centralized searchable directory to find existing AlgoRun-packaged algorithms. Availability and implementation: AlgoRun is available at http://algorun.org and the source code under GPL license is available at https://github.com/algorun Contact: laubenbacher@uchc.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Description: [ Pasteurella ] pneumotropica is a ubiquitous bacterium frequently isolated from laboratory rodents. Although this bacterium causes various diseases in immunosuppressed animals, little is known about major virulence factors and their roles in pathogenicity. To identify virulence factors, we sequenced the genome of [ P. ] pneumotropica biotype Heyl strain ATCC 12555, and compared the resulting non-contiguous draft genome sequence with the genome of biotype Jawetz strain ATCC 35149. Among a large number of genes encoding virulence-associated factors in both strains, four genes encoding for YadA-like proteins, which are known virulence factors that function in host cell adherence and invasion in many pathogens. In this study, we assessed YadA distribution and biological activity as an example of one of virulence–associated factor shared, with biotype Jawetz and Heyl. More than half of mouse isolates were found to have at least one of these genes; whereas, the majority of rat isolates did not. Autoagglutination activity, and ability to bind to mouse collagen type IV and mouse fibroblast cells, was significantly higher in YadA-positive than YadA-negative strains. To conclude, we identified a large number of candidate genes predicted to influence [ P. ] pneumotropica pathogenesis.

Description: Antibiotic resistance of bacteria growing in biofilms compared to their planktonic counterparts enhances the difficulty to eradicate biofilm-associated infections. In the last decade, combination antibiotic therapy has emerged as an attractive strategy for treating biofilm infections, even if in most of tolerant biofilms the optimal combinations are still unknown. In this study, an antimicrobial cationic polyacrylamide was used in combination with daptomycin or moxifloxacin against mature biofilms of Staphylococcus aureus clinical isolates to examine a possible improvement of the antibiofilm activity of the two antibiotics. The polymer did not have an effect on moxifloxacin but significantly increased the antibiofilm efficacy of daptomycin. These findings are presumably related to the different mechanism of action of the two drugs. In summary, our data highlighted the ability of polycations to increase daptomycin antibiofilm activity providing a potential strategy to eradicate biofilms in industrial or medical settings.

Description: Motivation: Genome browsers that support fast navigation through vast datasets and provide interactive visual analytics functions can help scientists achieve deeper insight into biological systems. Toward this end, we developed Integrated Genome Browser (IGB), a highly configurable, interactive and fast open source desktop genome browser. Results: Here we describe multiple updates to IGB, including all-new capabilities to display and interact with data from high-throughput sequencing experiments. To demonstrate, we describe example visualizations and analyses of datasets from RNA-Seq, ChIP-Seq and bisulfite sequencing experiments. Understanding results from genome-scale experiments requires viewing the data in the context of reference genome annotations and other related datasets. To facilitate this, we enhanced IGB’s ability to consume data from diverse sources, including Galaxy, Distributed Annotation and IGB-specific Quickload servers. To support future visualization needs as new genome-scale assays enter wide use, we transformed the IGB codebase into a modular, extensible platform for developers to create and deploy all-new visualizations of genomic data. Availability and implementation: IGB is open source and is freely available from http://bioviz.org/igb . Contact: aloraine@uncc.edu

Description: Motivation: Single Molecule Real-Time (SMRT) sequencing technology and Oxford Nanopore technologies (ONT) produce reads over 10 kb in length, which have enabled high-quality genome assembly at an affordable cost. However, at present, long reads have an error rate as high as 10–15%. Complex and computationally intensive pipelines are required to assemble such reads. Results: We present a new mapper, minimap and a de novo assembler, miniasm, for efficiently mapping and assembling SMRT and ONT reads without an error correction stage. They can often assemble a sequencing run of bacterial data into a single contig in a few minutes, and assemble 45-fold Caenorhabditis elegans data in 9 min, orders of magnitude faster than the existing pipelines, though the consensus sequence error rate is as high as raw reads. We also introduce a pairwise read mapping format and a graphical fragment assembly format, and demonstrate the interoperability between ours and current tools. Availability and implementation: https://github.com/lh3/minimap and https://github.com/lh3/miniasm Contact: hengli@broadinstitute.org Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Alternative splicing represents a prime mechanism of post-transcriptional gene regulation whose misregulation is associated with a broad range of human diseases. Despite the vast availability of transcriptome data from different cell types and diseases, bioinformatics-based surveys of alternative splicing patterns remain a major challenge due to limited availability of analytical tools that combine high accuracy and rapidity. Results: We describe here a novel junction-centric method, jSplice, that enables de novo extraction of alternative splicing events from RNA-sequencing data with high accuracy, reliability and speed. Application to clear cell renal carcinoma (ccRCC) cell lines and 65 ccRCC patients revealed experimentally validatable alternative splicing changes and signatures able to prognosticate ccRCC outcome. In the aggregate, our results propose jSplice as a key analytic tool for the derivation of cell context-dependent alternative splicing patterns from large-scale RNA-sequencing datasets. Availability and implementation: jSplice is a standalone Python application freely available at http://www.mhs.biol.ethz.ch/research/krek/jsplice . Contact: wilhelm.krek@biol.ethz.ch Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Single-cell RNA-sequencing technology allows detection of gene expression at the single-cell level. One typical feature of the data is a bimodality in the cellular distribution even for highly expressed genes, primarily caused by a proportion of non-expressing cells. The standard and the over-dispersed gamma-Poisson models that are commonly used in bulk-cell RNA-sequencing are not able to capture this property. Results: We introduce a beta-Poisson mixture model that can capture the bimodality of the single-cell gene expression distribution. We further integrate the model into the generalized linear model framework in order to perform differential expression analyses. The whole analytical procedure is called BPSC. The results from several real single-cell RNA-seq datasets indicate that ~90% of the transcripts are well characterized by the beta-Poisson model; the model-fit from BPSC is better than the fit of the standard gamma-Poisson model in 〉 80% of the transcripts. Moreover, in differential expression analyses of simulated and real datasets, BPSC performs well against edgeR, a conventional method widely used in bulk-cell RNA-sequencing data, and against scde and MAST, two recent methods specifically designed for single-cell RNA-seq data. Availability and Implementation: An R package BPSC for model fitting and differential expression analyses of single-cell RNA-seq data is available under GPL-3 license at https://github.com/nghiavtr/BPSC . Contact: yudi.pawitan@ki.se or mattias.rantalainen@ki.se Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Biological network querying is a problem requiring a considerable computational effort to be solved. Given a target and a query network, it aims to find occurrences of the query in the target by considering topological and node similarities (i.e. mismatches between nodes, edges, or node labels). Querying tools that deal with similarities are crucial in biological network analysis because they provide meaningful results also in case of noisy data. In addition, as the size of available networks increases steadily, existing algorithms and tools are becoming unsuitable. This is rising new challenges for the design of more efficient and accurate solutions. Results: This paper presents APPAGATO , a stochastic and parallel algorithm to find approximate occurrences of a query network in biological networks. APPAGATO handles node, edge and node label mismatches. Thanks to its randomic and parallel nature, it applies to large networks and, compared with existing tools, it provides higher performance as well as statistically significant more accurate results. Tests have been performed on protein–protein interaction networks annotated with synthetic and real gene ontology terms. Case studies have been done by querying protein complexes among different species and tissues. Availability and implementation: APPAGATO has been developed on top of CUDA-C ++ Toolkit 7.0 framework. The software is available online http://profs.sci.univr.it/~bombieri/APPAGATO . Contact: rosalba.giugno@univr.it Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : We introduce SharpViSu, an interactive open-source software with a graphical user interface, which allows performing processing steps for localization data in an integrated manner. This includes common features and new tools such as correction of chromatic aberrations, drift correction based on iterative cross-correlation calculations, selection of localization events, reconstruction of 2D and 3D datasets in different representations, estimation of resolution by Fourier ring correlation, clustering analysis based on Voronoi diagrams and Ripley’s functions. SharpViSu is optimized to work with eventlist tables exported from most popular localization software. We show applications of these on single and double-labelled super-resolution data. Availability and implementation: SharpViSu is available as open source code and as compiled stand-alone application under https://github.com/andronovl/SharpViSu . Contact: klaholz@igbmc.fr Supplementary information: Supplementary data are available at Bioinformatics online.

Description: The ICH S6(R1) recommendations on safety evaluation of biotherapeutics have led to uncertainty in determining what would constitute a cause for concern that would require genotoxicity testing. A Health and Environmental Sciences Institute’s Genetic Toxicology Technical Committee Workgroup was formed to review the current practice of genotoxicity assessment of peptide/protein-related biotherapeutics. There are a number of properties of peptide/protein-related biotherapeutics that distinguish such products from traditional ‘small molecule’ drugs and need to be taken into consideration when assessing whether genotoxicity testing may be warranted and if so, how to do it appropriately. Case examples were provided by participating companies and decision trees were elaborated to determine whether and when genotoxicity evaluation is needed for peptides containing natural amino acids, non-natural amino acids and other chemical entities and for unconjugated and conjugated proteins. From a scientific point of view, there is no reason for testing peptides containing exclusively natural amino acids irrespective of the manufacturing process. If non-natural amino acids, organic linkers and other non-linker chemical components have already been tested for genotoxicity, there is no need to re-evaluate them when used in different peptide/protein-related biotherapeutics. Unless the peptides have been modified to be able to enter the cells, it is generally more appropriate to evaluate the peptides containing the non-natural amino acids and other non-linker chemical moieties in vivo where the cleavage products can be formed. For linkers, it is important to determine if exposure to reactive forms are likely to occur and from which origin. When the linkers are anticipated to be potential mutagenic impurities they should be evaluated according to ICH M7. If linkers are expected to be catabolic products, it is recommended to test the entire conjugate in vivo , as this would ensure that the relevant ‘free’ linker forms stemming from in vivo catabolism are tested.

Description: Genome sequences that contain tandem repeats of guanine can form stable four-stranded structures known as G-quadruplex, or G4 DNA. While the molecular mechanisms are not fully defined, such guanine-rich loci are prone to mutagenesis and recombination. Various repair pathways function to reduce the potential for genome instability by correcting base damage and replication errors; however, it is not yet fully defined how well these processes function at G4 DNA. One frequent form of base damage occurs from cytidine deamination, resulting in deoxyuracil and UG mismatches. In duplex and single-stranded DNA, uracil bases are recognised and excised by uracil glycosylases. Here, we tested the efficiency of uracil glycosylase activity in vitro on uracil bases located directly adjacent to guanine repeats and G4 DNA. We show that uracil excision by bacterial UDG and human hUNG2 is reduced at uracils positioned directly 5' or 3' of a guanine tetrad. Control reactions using oligonucleotides disrupted for G4 formation or reaction conditions that do not favour G4 formation resulted in full uracil excision activity. Based on these in vitro results, we suggest that folding of guanine-rich DNA into G4 DNA results in a DNA conformation that is resistant to uracil glycosylase-initiated repair and this has the potential to increase the risk of instability at guanine repeats in the genome.

Description: Galectin-4 is a member of the galectin family which consists of 15 galactoside-binding proteins. Previously, galectin-4 has been shown to have a role in cancer progression and metastasis and it is found upregulated in many solid tumours, including colorectal cancer (CRC). Recently, the role in the metastatic process was suggested to be via promoting cancer cells to adhere to blood vascular endothelium. In the present study, the regulatory region of LGALS4 (galectin-4) in seven colon cell lines was investigated with respect to genetic variation that could be linked to expression levels and therefore a tumourigenic effect. Interestingly, qRT-PCR and sequencing results revealed that galectin-4 upregulation is associated with SNPs rs116896264 and rs73933062. By use of luciferase reporter- and pull-down assays, we confirmed the association between the gene upregulation and the two SNPs. Also, using pull-down assay followed by mass spectrometry, we found that the presence rs116896264 and rs73933062 is changing transcription factors binding sites. In order to assess the frequencies of the two SNPs among colon cancer patients and healthy individuals, we genotyped 75 colon cancer patients, 12 patients with adenomatous polyposis and 17 patients with ulcerative colitis and we performed data mining in the 1000 genomes databank. We found the two SNPs co-occuring in 21% of 75 CRC patients, 0 out of 12 patients of adenomatous polyposis, and 6 out of 17 patients (35%) with ulcerative colitis. Both in the patient samples and in the 1000 genomes project, the two SNPs were found to co-occur whenever present (D' = 1).

Description: The dose effect between nucleoplasmic bridges (NPB) and relatively low doses of ionising radiation remains unknown. Accordingly, this study investigated the NPB frequencies in human peripheral blood lymphocytes exposed to low-dose 60 Co -rays. Complex anomalies, including fused nuclei (FUS), horse-shoe nuclei (HS) and circular nuclei (CIR), which possibly originated from multiple NPBs, were also scored. Human peripheral blood samples were collected from three healthy males and irradiated with 0–1 and 0–0.4 Gy 60 Co -rays. A cytokinesis-block micronucleus cytome assay was then conducted to analyse NPB, PFHC (NPB plus three complex nuclear anomalies) and micronucleus (MN) in binucleated cells. All dose–response curves followed the linear model for both NPB frequency and PFHC cell frequency. The dose–response curves between NPB frequency and absorbed dose at 0–1 and 0–0.4 Gy were y = 0.0037 x + 0.0005 ( R 2 = 0.979, P 〈 0.05) and y = 0.0043 x + 0.0004 ( R 2 = 0.941, P 〈 0.05), respectively. The dose–response curves between PFHC cell frequency and absorbed dose at 0–1 and 0–0.4 Gy were y = 0.0044 x + 0.0007 ( R 2 = 0.982, P 〈 0.05) and y = 0.0059 x + 0.0005 ( R 2 = 0.969, P 〈 0.05), respectively. The statistical significance of differences between the irradiated groups (0–0.4 Gy) and background levels of NPB, PFHC and MN were also analysed. The lowest analysable doses of NPB, PFHC and MN were 0.12, 0.08 and 0.08 Gy, respectively. In conclusion, NPBs and PFHC positively correlated with the absorbed radiation at a relatively low dose.

Description: The aim of the study was to investigate how coadministration of resveratrol (RSV) at different time after the start of irradiation influences the frequency of micronuclei (MN) in reticulocytes of bone marrow and peripheral blood, and if the RSV supplementation after termination of irradiation may influence the recovery process of damaged cells. Coadministration of RSV with 1-day delay after 1 Gy irradiation significantly decreased the levels of MN in bone marrow and in peripheral blood, whereas with 1-week delay, only in bone marrow reticulocytes. Above combined treatment did not improve the process of recovery. RSV supplementation with 1-day delay relatively to 0.5 Gy irradiation, significantly decreased the frequencies of MN, especially after coadministration with 28mg/kg bw of RSV. Coadministration of RSV since eighth day did not influence the frequencies of MN compared to irradiated cells. The recovery process in the presence of RSV proceeded faster. Supplementation of RSV following initiation of irradiation is beneficial in case of irradiation with lower doses. RSV should be supplemented as soon as possible. Supplementation of RSV after termination of irradiation significantly speed up the recovery. Current results confirmed the ability of RSV to mitigate the effect of irradiation.

Description: Various naturally occurring stilbene-like compounds that are related to resveratrol (RSV) possess some of the beneficial effects of the parent molecule and provide even further benefits. Therefore, a series of methoxylated analogues of RSV were prepared with the aim of increasing antitumour and proapoptotic activity. In a previous article, we studied two methoxy-derivatives, pterostilbene (PTERO) and trimethoxystilbene (TRIMETHOXY), in which the first was formed by the substitution of two hydroxyl groups with two methoxy groups ( trans -3,5-dimethoxy-4'-hydroxystilbene) and the second was formed by the replacement of all three OH groups with methoxy groups ( trans -3,5,4'-trimethoxystilbene). Both methoxy-derivatives showed stronger antioxidant activity when compared with RSV. In the present article, we focused on the analysis of the ability of RSV and its two methoxylated derivatives to protect proliferating non-tumoural cells from the damage induced by ionising radiation (IR). First we showed that the methoxy derivatives, contrary to their parental compound, are unable to affect topoisomerase enzyme and consequently are not clastogenic per se . Second we showed that both PTERO and TRIMETHOXY more efficiently reduce the chromosome damage induced by IR. Furthermore, TRIMETHOXY, but not PTERO, causes a delay in cell proliferation, particularly in mitosis progression increasing the number of cells in metaphase at the expense of prophases and ana/telophases.

Description: α-, β- and -asarone are naturally occurring phenylpropenes that occur in different plant families, mainly in Aristolochiaceae , Acoraceae and Lauraceae. Plants containing asarones are used as flavouring ingredients in alcoholic beverages (bitters), traditional phytomedicines and the rhizome of e.g. Acorus calamus is used to prepare tea. Although α- and β-asarone show a potential in the treatment of several diseases, previous studies have shown carcinogenicity in rodents (duodenum, liver). However, the mechanism of action remained unclear. Studies on the mutagenicity of propenylic α- and β-asarone are inconsistent and data on carcinogenicity and genotoxicity of allylic -asarone are lacking completely. Thus, the present study determined the mutagenicity of the three asarone isomers using the Ames fluctuation assay with and without exogenous metabolic activation (S9 mix) in the standard Salmonella typhimurium strains TA98 and TA100. A concentration dependent increase in mutagenicity could be verified for α- and β-asarone in strain TA100 in the presence of rat liver homogenate. The side-chain epoxides of α- and β-asarone, major metabolites formed in liver microsomes, caused mutations in TA100, supporting the hypothesis that epoxidation of the side chain plays a key role in mutagenicity of the propenylic alkenylbenzenes. The allylic -asarone, not undergoing detectable side-chain epoxidation in liver microsomes, was supposed to be activated via side-chain hydroxylation and further sulphonation, a typical pathway for other allylic alkenylbenzenes like estragole or methyleugenol. However, neither y-asarone nor 1'-OH--asarone showed any mutagenic effect even in the human SULT-expressing Salmonella strains (TA100-hSULT1A1 and TA100-hSULT1C2), while 1'-OH-methyleugenol used as a positive control was mutagenic under these conditions. These results indicate that the propenylic asarones are genotoxic via metabolic formation of side-chain epoxides while the side-chain hydroxylation/sulphonation pathway is either not operative or does not lead to mutagenicity with the allylic -asarone.

Description: Prior to the downstream development of chemical substances, including pharmaceuticals and cosmetics, their influence on the genetic apparatus has to be tested. Several in vitro and in vivo assays have been developed to test for genotoxicity. In a first tier, a battery of two to three in vitro tests is recommended to cover mutagenicity, clastogenicity and aneugenicity as main endpoints. This regulatory in vitro test battery is known to have a high sensitivity, which is at the expense of the specificity. The high number of false positive in vitro results leads to excessive in vivo follow-up studies. In the case of cosmetics it may even induce the ban of the particular compound since in Europe the use of experimental animals is no longer allowed for cosmetics. In this article, an alternative approach to derisk a misleading positive Ames test is explored. Hereto we first tested the performance of five existing computational tools to predict the potential mutagenicity of a data set of 132 cosmetic compounds with a known genotoxicity profile. Furthermore, we present, as a proof-of-principle, a strategy in which a combination of computational tools and mechanistic information derived from in vitro transcriptomics analyses is used to derisk a misleading positive Ames test result. Our data shows that this strategy may represent a valuable tool in a weight-of-evidence approach to further evaluate a positive outcome in an Ames test.

Description: Pioglitazone (PTZ) is an oral antidiabetic agent whose anti-cancer properties have been described recently. Since PTZ increases the production of reactive oxygen species in mammalian cells, the aim of current study was to evaluate the cytotoxic, mutagenic and recombinogenic effects of PTZ using respectively the in vitro mitotic index assay and the in vitro mammalian cell micronucleus test in human peripheral lymphocytes, and the in vivo homozygotization assay in Aspergillus nidulans , which detects the loss of heterozygosity due to somatic recombination. Although the lowest PTZ concentrations (4–36 μM) did not show any significant rise in the micronucleus production, the higher PTZ concentration (108 μM) produced a statistically higher number of micronuclei than the negative control and significantly altered the cell-proliferation kinetics, demonstrating the mutagenic and antiproliferative effects of PTZ, respectively. The recombinogenic activity of PTZ, demonstrated here for the first time, was observed at the highest tested concentration (400 μM) through the homozygotization index rates significantly different from the negative control. Taken together, our results show that PTZ is genotoxic at a concentration higher than the therapeutic plasma concentration. This PTZ genotoxicity may be a potential benefit to its previously described antitumour activity.

Description: Environmental pollutants are complex mixtures in which metals are ubiquitous. Metal mixtures of arsenic, cadmium and lead are present in the occupational environment and generate health effects such as cardiovascular, renal and cancer diseases. Cell transformation induced by metal mixtures that depend on reactive oxygen species (ROS) generation, cell viability maintenance and avoidance of senescence was previously reported by our group. The aim of the present study was to explore the role of a Obg-like ATPase1 (OLA1) in the cell transformation of BALB/c 3T3 A31-1-1 clonal cells induced by a metal mixture (2 µM NaAsO 2 , 2 µM CdCl 2 and 5 µM Pb(C 2 H 3 O 2 ) 2 . 3H 2 O) through ROS generation. The interest in OLA1 is justified because this protein has been proposed to be a negative regulator of the cellular antioxidant response. Small interfering RNA (siRNA) was used to knockdown OLA1 before the initiation stage of the transformation assay. We evaluated (ROS) and OLA1 protein expression throughout the initiation and promotion stages of transformation. OLA1 knockdown modulated metal mixture-induced cell transformation more strongly when the metal mixture was an initiator stimulus than when it was a promoter. The ability of the metal mixture to initiate cell transformation was diminished by OLA1 knockdown, an effect that depended on intracellular ROS levels. The effect of OLA1 was synergistic with N -Acetyl- l -cysteine (NAC) co-treatment. Oxidative stress-associated transcription factors Egr1 and Smad were also down-regulated by the OLA1 knockdown, contributing to the rescue of metal mixture cell transformation.

Description: G-quadruplexes (G4) are highly stable tetra-stranded DNA secondary structures known to mediate gene regulation and to trigger genomic instability events during replication. G4 structural stability can be affected by DNA methylation and oxidation modifications; thus nutrients such as folate that have the ability to alter these processes could potentially modify the genomic occurrence of G4 elements. Hela cells were cultured in a range of folate concentrations or in the presence or absence of 5-aza-2'-deoxycytidine, a DNA-methyltransferase inhibitor. G4 structures were then quantified by immunofluorescence using an automated quantitative imaging system. G4 frequency in Hela cells and nuclei area mean were increased in 20nM folate medium compared with 2000nM folate, as well as in the presence of 5-aza-2'-deoxycytidine when compared to cells non-exposed to 5-aza-2'-deoxycytidine. These changes were exacerbated when pyridostatin, a G4 stabilising ligand, was added to the culture medium. G4 intensity in Hela cells cultured in deficient folate condition with pyridostatin was highly correlated with DNA damage as measured by H2AX immunofluorescence ( r = 0.71). This study showed for the first time that cellular G4 balance is modifiable by low folate concentrations and that these changes may occur as a consequence of DNA hypomethylation. Although the exact mechanism by which these changes occur is unclear, these findings establish the possibility that nutrients could be utilised as a tool for sustaining genome integrity by modifying G4 frequency at a cellular level.

Description: Metal oxide nanoparticles (NPs), including zinc oxide (ZnO) NPs have shown success for use as vehicles for drug delivery and targeting gene delivery in many diseases like cancer. Current anticancer chemotherapeutics fail to effectively differentiate between cancerous and normal cells. There is an urgent need to develop novel drug delivery system that can better target cancer cells while sparing normal cells and tissues. Particularly, ZnO NPs exhibit a high degree of cancer cell selectivity and induce cell death, oxidative stress, interference with the cell cycle progression and genotoxicity in cancerous cells. In this scenario, effective cellular uptake of NP seems to be crucial, which is shown to be affected by cell cycle progression. In the present study, the cytotoxic potential of ZnO NPs and the effect of different cell cycle phases on the uptake of ZnO NPs were examined in A431 cells. It is shown that the ZnO NPs led to cell death and reactive oxygen species generation and were able to induce cell cycle arrest in S and G 2 /M phase with the higher uptake in G 2 /M phase compared with other phases.

Description: Malondialdehyde (MDA), a biomarker of lipid peroxidation and oxidative stress, is a mutagenic and carcinogenic compound that can react with DNA to form several types of DNA adducts including the deoxyguanosine adduct (M 1 dG). The aim of this cross-sectional study was to evaluate the association between individual dietary and lifestyle habits and M 1 dG levels, measured in peripheral leukocytes in a large representative sample of the general population of Florence City (Italy). Selected anthropometric measurements, detailed information on dietary and lifestyle habits and blood samples were available for 313 adults of the Florence City Sample enrolled in the frame of European Prospective Investigation into Cancer and nutrition (EPIC) study. A multivariate regression analysis adjusted for selected individual characteristics possibly related to M 1 dG levels (sex, age, BMI, smoke, physical activity level, education level, total caloric intake and a Mediterranean dietary score) was performed to estimate the association between these parameters and M 1 dG levels. M 1 dG levels were significantly higher in women ( P = 0.014) and lower in moderately active or active subjects ( P = 0.037).We also found a significant inverse association with the Modified Mediterranean dietary score ( P for trend = 0.049), particularly evident for the highest categories of adherence. Our results indicate that M 1 dG levels can be modulated by selected individual characteristics such as gender, physical activity and a Mediterranean dietary pattern.

Description: : A gene tree-species tree reconciliation explains the evolution of a gene tree within the species tree given a model of gene-family evolution. We describe ecceTERA, a program that implements a generic parsimony reconciliation algorithm, which accounts for gene duplication, loss and transfer (DTL) as well as speciation, involving sampled and unsampled lineages, within undated, fully dated or partially dated species trees. The ecceTERA reconciliation model and algorithm generalize or improve upon most published DTL parsimony algorithms for binary species trees and binary gene trees. Moreover, ecceTERA can estimate accurate species-tree aware gene trees using amalgamation. Availability and implementation : ecceTERA is freely available under http://mbb.univ-montp2.fr/MBB/download_sources/16__ecceTERA and can be run online at http://mbb.univ-montp2.fr/MBB/subsection/softExec.php?soft=eccetera . Contact: celine.scornavacca@umontpellier.fr Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : The popularity of using NMR spectroscopy in metabolomics and natural products has driven the development of an array of NMR spectral analysis tools and databases. Particularly, web applications are well used recently because they are platform-independent and easy to extend through reusable web components. Currently available web applications provide the analysis of NMR spectra. However, they still lack the necessary processing and interactive visualization functionalities. To overcome these limitations, we present NMRPro, a web component that can be easily incorporated into current web applications, enabling easy-to-use online interactive processing and visualization. NMRPro integrates server-side processing with client-side interactive visualization through three parts: a python package to efficiently process large NMR datasets on the server-side, a Django App managing server-client interaction, and SpecdrawJS for client-side interactive visualization. Availability and implementation: Demo and installation instructions are available at http://mamitsukalab.org/tools/nmrpro/ Contact: mohamed@kuicr.kyoto-u.ac.jp Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Streptococcus mutans is a Gram-positive bacterium involved in development to caries, the most common infectious disease of our time. Streptococcus mutans interacts with other microbes, like the fungus Candida albicans and both are commonly isolated from patients with caries. Since the role of C. albicans in caries remains unknown, our aim was to unravel this using an in vitro dual-species cariogenic oral biofilm model. Biofilms were grown for 24–72 h on glass cover slips or hydroxyapatite (HA) disks to mimic the surface of teeth. Medium pH, lactic acid production capacity and calcium release from HA disks were determined. All 24-h biofilms had external pH values below the critical pH of 5.5 where enamel dissolves. In contrast, 72-h dual-species biofilms had significantly higher pH (above the critical pH) and consequently decreased calcium release compared to single-species S. mutans biofilms. Counter intuitively, lactic acid production and growth of S. mutans were increased in 72-h dual-species biofilms. Candida albicans modulates the pH in dual-species biofilms to values above the critical pH where enamel dissolves. Our results suggest that C. albicans is not by definition a cariogenic microorganism; it could prevent caries by actively increasing pH preventing mineral loss.

Description: Emergence of drug resistance during visceral leishmaniasis (VL) is a major obstacle imposed during successful therapy. An effective vaccine strategy against this disease is therefore necessary. Our present study exploited the SLA (soluble leishmanial antigen) and PGN (peptidoglycan) stimulated bone marrow-derived dendritic cells (DCs) as a suitable vaccine candidate during experimental VL. SLA-PGN-stimulated DCs showed a significant decrease in hepatic and splenic parasite burden, which were associated with increased production of nitric oxide and pro-inflammatory cytokines such as IL-12, IFN- and IL-17. Elevated level of IL-17 was accompanied with the generation of more Th17 cells. Further studies on DC provided the evidence that these SLA-PGN-stimulated DCs played an important role in providing necessary cytokines such as IL-6, IL-23 and TGF-β for the generation of Th17 cells. Interestingly, inhibition of protein kinase C-β (PKCβ) in DCs led to decreased production of Th17 polarizing cytokines, causing reduction of the Th17 population size. Altogether, our finding highlighted the important role of DC-based PKCβ in regulation of the function and generation of Th17 cells.

Description: TLR2 heterodimers with TLR1 or TLR6 recognize distinct pathogen-associated molecules such as tri- and di-acylated lipopeptides. The activated TLR2 heterodimers recruit Toll-IL-1R domain- (TIR-) containing adapter proteins, TIRAP and MyD88, through the receptor TIR domains. Molecular recognition mechanisms responsible for agonist-driven, TIR domain-mediated receptor–adapter interactions as well as the structure of resultant signaling complexes remain unknown. We previously reported that the cell-permeable peptide derived from helix D of TLR2 TIR (2R9) specifically binds TIRAP in vitro and in cells and thereby inhibits TIRAP-dependent TLR signaling. This study demonstrates that cell-permeable peptides from D helix of TLR1 or TLR6, peptides 1R9 and 6R9 respectively, inhibit signaling mediated by cognate TLR2 co-receptors. Interestingly, 1R9 and 6R9 bind different TLR2 adapters, as they selectively bind MyD88 and TIRAP TIR, respectively. Both peptides block the agonist-induced co-immunoprecipitation (co-IP) of TLR2 with TIRAP or MyD88, but not TLR2 co-IP with co-receptors. Our data suggest that D helices of TLR1 and TLR6 TIR domains are adapter recruitment sites in both co-receptors; yet the sites recruit different adapters. The D helix in TLR1 is the MyD88 docking site, whereas in TLR6 this site recruits TIRAP.

Description: Vibrio parahaemolyticus is a Gram-negative, motile, nonspore-forming pathogen that causes foodborne illness associated with the consumption of contaminated seafoods. Although many cases of foodborne outbreaks caused by V. parahaemolyticus have been reported, the genomes of only five strains have been completely sequenced and analyzed using bioinformatics. In order to characterize overall virulence factors and pathogenesis of V. parahaemolyticus associated with foodborne outbreak in South Korea, a new strain FORC_008 was isolated from flounder fish and its genome was completely sequenced. The genomic analysis revealed that the genome of FORC_008 consists of two circular DNA chromosomes of 3266 132 bp (chromosome I) and 1772 036 bp (chromosome II) with a GC content of 45.36% and 45.53%, respectively. The entire genome contains 4494 predicted open reading frames, 129 tRNAs and 31 rRNA genes. While the strain FORC_008 does not have genes encoding thermostable direct hemolysin (TDH) and TDH-related hemolysin (TRH), its genome encodes many other virulence factors including hemolysins, pathogenesis-associated secretion systems and iron acquisition systems, suggesting that it may be a potential pathogen. This report provides an extended understanding on V. parahaemolyticus in genomic level and would be helpful for rapid detection, epidemiological investigation and prevention of foodborne outbreak in South Korea.

Description: Candida glabrata is an opportunistic pathogen that is considered the second most common cause of candidiasis after Candida albicans . Many characteristics of its mechanisms of pathogenicity remain unknown. Recent studies have focused on determining the events that underlie interactions between C. glabrata and immune cells, but the relationship between this yeast and osteoblasts has not been studied in detail. The aim of this study was to determine the mechanisms of interaction between human osteoblasts and C. glabrata , and to identify the roles played by some of the molecules that are produced by these cells in response to infection. We show that C. glabrata adheres to and is internalized by human osteoblasts. Adhesion is independent of opsonization, and internalization depends on the rearrangement of the actin cytoskeleton. We show that C. glabrata survives and replicates in osteoblasts and that this intracellular behavior is related to the level of production of nitric oxide and reactive oxygen species. Opsonized C. glabrata stimulates the production of IL-6, IL-8 and MCP-1 cytokines. Adhesion and internalization of the pathogen and the innate immune response of osteoblasts require viable C. glabrata . These results suggest that C. glabrata modulates immunological mechanisms in osteoblasts to survive inside the cell.

Description: Biofilms are complex communities of microorganisms embedded in an extracellular matrix and adherent to a surface. The development was described as a four-stage process leading to the formation of a mature biofilm which was resistant to immune system and antibiotic actions. In bone and joint infections (BJIs), the formation of biofilms is a leading cause of treatment failure. Here we study the capacity of 11 antibiotics commonly used in the treatment of BJIs to inhibit the biofilm formation on 29 clinical Staphylococcus aureus isolates by a new test called Antibiofilmogram ® . The minimal inhibitory concentration (MIC) and biofilm MIC (bMIC) were determined in vitro and showed similar values for clindamycin, fusidic acid, linezolid and rifampin. Reversely, daptomycin, fosfomycin, gentamicin and ofloxacin showed a bMIC distribution different from MIC with bMIC above breakpoint. Finally, cloxacillin, teicoplanin and vancomycin revealed an intermediate bMIC distribution with a strain-dependent pattern. A murine in vivo model of catheter-associated S. aureus infection was made and showed a significant reduction, but not total prevention, of catheter colonization with cloxacillin at bMIC, and no or limited reduction with cloxacillin at MIC. Antibiofilmogram ® could be of great interest after surgical operations on contaminated prostheses and after bacteremia in order to prevent the colonization of the device.

Description: Many prokaryotes utilize type IV secretion systems (T4SSs) to translocate substrates (e.g. nucleoprotein, DNA, protein) across the cell envelope, and/or to elaborate surface structures (i.e. pili or adhesins). Among eight distinct T4SS classes, P-T4SSs are typified by the Agrobacterium tumefaciens vir T4SS, which is comprised of 12 scaffold components (VirB1–VirB11, VirD4). While most P-T4SSs include all 12 Vir proteins, some differ from the vir archetype by either containing additional scaffold components not analogous to Vir proteins or lacking one or more of the Vir proteins. In a special case, the Rickettsiales vir homolog ( rvh ) P-T4SS comprises unprecedented gene family expansion. rvh contains three families of gene duplications ( rvhB9 , rvhB8 , rvhB4 ): RvhB9,8,4-I are conserved relative to equivalents in other P-T4SSs, while RvhB9,8,4-II have evolved atypical features that deviate substantially from other homologs. Furthermore, rvh contains five VirB6-like genes ( rvhB6a – e ), which are tandemly arrayed and contain large N- and C-terminal extensions. Our work herein focuses on the complexity underpinned by rvh gene family expansion. Furthermore, we describe an RvhB10 insertion, which occurs in a region that forms the T4SS pore. The significance of these curious properties to rvh structure and function is evaluated, shedding light on a highly complex T4SS.

Description: Enteroaggregative Escherichia coli (EAEC) causes acute and persistent diarrhea among children, HIV-infected patients, and travelers to developing countries. We have searched for 18 genes-encoding virulence factors associated with aggregative adherence, dispersion, biofilm, toxins, serine protease autotransporters of Enterobacteriaceae (SPATEs) and siderophores, analyzed in 172 well-characterized EAEC strains ( aggR + ) isolated from stool samples of 97 children with diarrhea and 75 healthy controls from a passive surveillance diarrhea cohort study in Peru. Eighty-one different genetic profiles were identified, 37 were found only associated with diarrhea and 25 with control samples. The most frequent genetic profile was aggC + aatA + aap + shf + fyuA + , present in 19 strains, including diarrhea and controls. The profile set1A + set1B + pic + was associated with diarrhea ( P 〈 0.05). Of all genes evaluated, the most frequent were aatA (CVD 342) present in 159 strains (92.4%) and fyuA in 157 (91.3%). When EAEC strains were analyzed as a single pathogen (excluding co-infections), only pic was associated with diarrhea ( P 〈 0.05) and with prolonged diarrhea (diarrhea ≥ 7 days) ( P 〈 0.05). In summary, this is the first report on the prevalence of a large set of EAEC virulence genes and its association with diarrhea in Peruvian children. More studies are needed to elucidate the exact role of each virulence factor.

Description: Natural killer T (NKT) cells play a critical role in the host's innate immune response. CD1d-mediated presentation of glycolipid antigens to NKT cells has been established; however, the mechanisms by which NKT cells recognize infected or cancerous cells remain unclear. 5 ' -AMP activated protein kinase (AMPK) is a master regulator of lipogenic pathways. We hypothesized that activation of AMPK during infection and malignancy could alter the repertoire of antigens presented by CD1d and serve as a danger signal to NKT cells. In this study, we examined the effect of alterations in metabolism on CD1d-mediated antigen presentation to NKT cells and found that an infection with lymphocytic choriomeningitis virus rapidly increased CD1d-mediated antigen presentation. Hypoxia inducible factors (HIF) enhance T-cell effector functions during infection, therefore antigen presenting cells pretreated with pharmacological agents that inhibit glycolysis, induce HIF and activate AMPK were assessed for their ability to induce NKT-cell responses. Pretreatment with 2-deoxyglucose, cobalt chloride, AICAR and metformin significantly enhanced CD1d-mediated NKT-cell activation. In addition, NKT cells preferentially respond to malignant B cells and B-cell lymphomas express HIF-1α. These data suggest that targeting cellular metabolism may serve as a novel means of inducing innate immune responses.

Description: The active subunit (S1) of pertussis toxin (PT), a major virulence factor of Bordetella pertussis , ADP-ribosylates G i proteins in the mammalian cell cytosol to inhibit GPCR signaling. The intracellular pathway of PT includes endocytosis and retrograde transport to the trans-Golgi network (TGN) and endoplasmic reticulum (ER). Subsequent translocation of S1 to the cytosol is presumably preceded by dissociation from the holotoxin. In vitro , such dissociation is stimulated by interaction of PT with ATP. To investigate the role of this interaction in cellular events, we engineered a form of PT (PTDM) with changes to two amino acids involved in the interaction with ATP. PTDM was reduced in (1) binding to ATP, (2) dissociability by interaction with ATP, (3) in vitro enzymatic activity and (4) cellular ADP-ribosylation activity. In cells treated with PTDM carrying target sequences for organelle-specific modifications, normal transport to the TGN and ER occurred, but N-glycosylation patterns of the S1 and S4 subunits were consistent with an inability of PTDM to dissociate in the ER. These results indicate a requirement for interaction with ATP for PT dissociation in the ER and cellular activity. They also indicate that the retrograde transport route is the cellular intoxication pathway for PT.

Description: Sphingosine 1-phosphate (S1P) is a sphingosine containing lipid intermediate obtained from ceramide. S1P is known to be an important signaling molecule and plays multiple roles in the context of immunity. This lysophospholipid binds and activates G-protein-coupled receptors (GPCRs) known as S1P receptors 1–5 (S1P 1–5 ). Once activated, these GPCRs mediate signaling that can lead to alterations in cell proliferation, survival or migration, and can also have other effects such as promoting angiogenesis. In this review, we will present evidence demonstrating a role for S1P in lymphocyte migration, inflammation and infection, as well as in cancer. The therapeutic potential of targeting S1P receptors, kinases and lyase will also be discussed.

Description: Increasing morbidity and mortality from Clostridium difficile infection (CDI) present an enormous challenge to healthcare systems. Clostridium difficile express type IV pili (T4P), but their function remains unclear. Many chronic and recurrent bacterial infections result from biofilms, surface-associated bacterial communities embedded in an extracellular matrix. CDI may be biofilm mediated; T4P are important for biofilm formation in a number of organisms. We evaluate the role of T4P in C. difficile biofilm formation using RNA sequencing, mutagenesis and complementation of the gene encoding the major pilin pilA1 , and microscopy. RNA sequencing demonstrates that, in comparison to other growth phenotypes, C. difficile growing in a biofilm has a distinct RNA expression profile, with significant differences in T4P gene expression. Microscopy of T4P-expressing and T4P-deficient strains suggests that T4P play an important role in early biofilm formation. A non-piliated pilA1 mutant forms an initial biofilm of significantly reduced mass and thickness in comparison to the wild type. Complementation of the pilA1 mutant strain leads to formation of a biofilm which resembles the wild-type biofilm. These findings suggest that T4P play an important role in early biofilm formation. Novel strategies for confronting biofilm infections are emerging; our data suggest that similar strategies should be investigated in CDI.

Description: Motivation: DNA methylation is an important epigenetic modification related to a variety of diseases including cancers. We focus on the methylation data from Illumina’s Infinium HumanMethylation450 BeadChip. One of the key issues of methylation analysis is to detect the differential methylation sites between case and control groups. Previous approaches describe data with simple summary statistics or kernel function, and then use statistical tests to determine the difference. However, a summary statistics-based approach cannot capture complicated underlying structure, and a kernel function-based approach lacks interpretability of results. Results: We propose a novel method D 3 M, for detection of differential distribution of methylation, based on distribution-valued data. Our method can detect the differences in high-order moments, such as shapes of underlying distributions in methylation profiles, based on the Wasserstein metric. We test the significance of the difference between case and control groups and provide an interpretable summary of the results. The simulation results show that the proposed method achieves promising accuracy and shows favorable results compared with previous methods. Glioblastoma multiforme and lower grade glioma data from The Cancer Genome Atlas show that our method supports recent biological advances and suggests new insights. Availability and Implementation: R implemented code is freely available from https://github.com/ymatts/D3M/ . Contact: ymatsui@med.nagoya-u.ac.jp or shimamura@med.nagoya-u.ac.jp Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Similarity-based methods have been widely used in order to infer the properties of genes and gene products containing little or no experimental annotation. New approaches that overcome the limitations of methods that rely solely upon sequence similarity are attracting increased attention. One of these novel approaches is to use the organization of the structural domains in proteins. Results: We propose a method for the automatic annotation of protein sequences in the UniProt Knowledgebase (UniProtKB) by comparing their domain architectures, classifying proteins based on the similarities and propagating functional annotation. The performance of this method was measured through a cross-validation analysis using the Gene Ontology (GO) annotation of a sub-set of UniProtKB/Swiss-Prot. The results demonstrate the effectiveness of this approach in detecting functional similarity with an average F-score: 0.85. We applied the method on nearly 55.3 million uncharacterized proteins in UniProtKB/TrEMBL resulted in 44 818 178 GO term predictions for 12 172 114 proteins. 22% of these predictions were for 2 812 016 previously non-annotated protein entries indicating the significance of the value added by this approach. Availability and implementation: The results of the method are available at: ftp://ftp.ebi.ac.uk/pub/contrib/martin/DAAC/ . Contact: tdogan@ebi.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Species identification and quantification are common tasks in metagenomics and pathogen detection studies. The most recent techniques are built on mapping the sequenced reads against a reference database (e.g. whole genomes, marker genes, proteins) followed by application-dependent analysis steps. Although these methods have been proven to be useful in many scenarios, there is still room for improvement in species and strain level detection, mainly for low abundant organisms. Results: We propose a new method: DUDes, a reference-based taxonomic profiler that introduces a novel top-down approach to analyze metagenomic Next-generation sequencing (NGS) samples. Rather than predicting an organism presence in the sample based only on relative abundances, DUDes first identifies possible candidates by comparing the strength of the read mapping in each node of the taxonomic tree in an iterative manner. Instead of using the lowest common ancestor we propose a new approach: the deepest uncommon descendent. We showed in experiments that DUDes works for single and multiple organisms and can identify low abundant taxonomic groups with high precision. Availability and Implementation: DUDes is open source and it is available at http://sf.net/p/dudes Supplementary information: Supplementary data are available at Bioinformatics online. Contact: renardB@rki.de

Description: Motivation: Moonlighting proteins (MPs) show multiple cellular functions within a single polypeptide chain. To understand the overall landscape of their functional diversity, it is important to establish a computational method that can identify MPs on a genome scale. Previously, we have systematically characterized MPs using functional and omics-scale information. In this work, we develop a computational prediction model for automatic identification of MPs using a diverse range of protein association information. Results: We incorporated a diverse range of protein association information to extract characteristic features of MPs, which range from gene ontology (GO), protein–protein interactions, gene expression, phylogenetic profiles, genetic interactions and network-based graph properties to protein structural properties, i.e. intrinsically disordered regions in the protein chain. Then, we used machine learning classifiers using the broad feature space for predicting MPs. Because many known MPs lack some proteomic features, we developed an imputation technique to fill such missing features. Results on the control dataset show that MPs can be predicted with over 98% accuracy when GO terms are available. Furthermore, using only the omics-based features the method can still identify MPs with over 75% accuracy. Last, we applied the method on three genomes: Saccharomyces cerevisiae , Caenorhabditis elegans and Homo sapiens , and found that about 2–10% of proteins in the genomes are potential MPs. Availability and Implementation: Code available at http://kiharalab.org/MPprediction Contact: dkihara@purdue.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Design of protein–protein interaction (PPI) inhibitors is a major challenge in Structural Bioinformatics. Peptides, especially short ones (5–15 amino acid long), are natural candidates for inhibition of protein–protein complexes due to several attractive features such as high structural compatibility with the protein binding site (mimicking the surface of one of the proteins), small size and the ability to form strong hotspot binding connections with the protein surface. Efficient rational peptide design is still a major challenge in computer aided drug design, due to the huge space of possible sequences, which is exponential in the length of the peptide, and the high flexibility of peptide conformations. Results: In this article we present PinaColada, a novel computational method for the design of peptide inhibitors for protein–protein interactions. We employ a version of the ant colony optimization heuristic, which is used to explore the exponential space ( 20n ) of length n peptide sequences, in combination with our fast robotics motivated PepCrawler algorithm, which explores the conformational space for each candidate sequence. PinaColada is being run in parallel, on a DELL PowerEdge 2.8 GHZ computer with 20 cores and 256 GB memory, and takes up to 24 h to design a peptide of 5-15 amino acids length. Availability and implementation: An online server available at: http://bioinfo3d.cs.tau.ac.il/PinaColada/. Contact: danielza@post.tau.ac.il ; wolfson@tau.ac.il

Description: Motivation: Whole-genome low-coverage sequencing has been combined with linkage-disequilibrium (LD)-based genotype refinement to accurately and cost-effectively infer genotypes in large cohorts of individuals. Most genotype refinement methods are based on hidden Markov models, which are accurate but computationally expensive. We introduce an algorithm that models LD using a simple multivariate Gaussian distribution. The key feature of our algorithm is its speed. Results: Our method is hundreds of times faster than other methods on the same data set and its scaling behaviour is linear in the number of samples. We demonstrate the performance of the method on both low- and high-coverage samples. Availability and implementation: The source code is available at https://github.com/illumina/marvin Contact: rarthur@illumina.com Supplementary information: Supplementary data are available at Bioinformatics online .

Description: Motivation : T-cell epitopes serve as molecular keys to initiate adaptive immune responses. Identification of T-cell epitopes is also a key step in rational vaccine design. Most available methods are driven by informatics and are critically dependent on experimentally obtained training data. Analysis of a training set from Immune Epitope Database (IEDB) for several alleles indicates that the sampling of the peptide space is extremely sparse covering a tiny fraction of the possible nonamer space, and also heavily skewed, thus restricting the range of epitope prediction. Results : We present a new epitope prediction method that has four distinct computational modules: (i) structural modelling, estimating statistical pair-potentials and constraint derivation, (ii) implicit modelling and interaction profiling, (iii) feature representation and binding affinity prediction and (iv) use of graphical models to extract peptide sequence signatures to predict epitopes for HLA class I alleles. Conclusions : HLaffy is a novel and efficient epitope prediction method that predicts epitopes for any Class-1 HLA allele, by estimating the binding strengths of peptide-HLA complexes which is achieved through learning pair-potentials important for peptide binding. It relies on the strength of the mechanistic understanding of peptide-HLA recognition and provides an estimate of the total ligand space for each allele. The performance of HLaffy is seen to be superior to the currently available methods. Availability and implementation : The method is made accessible through a webserver http://proline.biochem.iisc.ernet.in/HLaffy . Contact : nchandra@biochem.iisc.ernet.in Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: In vitro and in vivo cell proliferation is often studied using the dye carboxyfluorescein succinimidyl ester (CFSE). The CFSE time-series data provide information about the proliferation history of populations of cells. While the experimental procedures are well established and widely used, the analysis of CFSE time-series data is still challenging. Many available analysis tools do not account for cell age and employ optimization methods that are inefficient (or even unreliable). Results: We present a new model-based analysis method for CFSE time-series data. This method uses a flexible description of proliferating cell populations, namely, a division-, age- and label-structured population model. Efficient maximum likelihood and Bayesian estimation algorithms are introduced to infer the model parameters and their uncertainties. These methods exploit the forward sensitivity equations of the underlying partial differential equation model for efficient and accurate gradient calculation, thereby improving computational efficiency and reliability compared with alternative approaches and accelerating uncertainty analysis. The performance of the method is assessed by studying a dataset for immune cell proliferation. This revealed the importance of different factors on the proliferation rates of individual cells. Among others, the predominate effect of cell age on the division rate is found, which was not revealed by available computational methods. Availability and implementation: The MATLAB source code implementing the models and algorithms is available from http://janhasenauer.github.io/ShAPE-DALSP/ . Contact: jan.hasenauer@helmholtz-muenchen.de Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: The challenges of successfully applying causal inference methods include: (i) satisfying underlying assumptions, (ii) limitations in data/models accommodated by the software and (iii) low power of common multiple testing approaches. Results: The causal inference test (CIT) is based on hypothesis testing rather than estimation, allowing the testable assumptions to be evaluated in the determination of statistical significance. A user-friendly software package provides P -values and optionally permutation-based FDR estimates ( q -values) for potential mediators. It can handle single and multiple binary and continuous instrumental variables, binary or continuous outcome variables and adjustment covariates. Also, the permutation-based FDR option provides a non-parametric implementation. Conclusion: Simulation studies demonstrate the validity of the cit package and show a substantial advantage of permutation-based FDR over other common multiple testing strategies. Availability and implementation: The cit open-source R package is freely available from the CRAN website ( https://cran.r-project.org/web/packages/cit/index.html ) with embedded C ++ code that utilizes the GNU Scientific Library, also freely available ( http://www.gnu.org/software/gsl/ ). Contact: joshua.millstein@usc.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: The vast majority of the many thousands of disease-associated single nucleotide polymorphisms (SNPs) lie in the non-coding part of the genome. They are likely to affect regulatory elements, such as enhancers and promoters, rather than the function of a protein. To understand the molecular mechanisms underlying genetic diseases, it is therefore increasingly important to study the effect of a SNP on nearby molecular traits such as chromatin or transcription factor binding. Results: We developed SNPhood , a user-friendly Bioconductor R package to investigate, quantify and visualise the local epigenetic neighbourhood of a set of SNPs in terms of chromatin marks or TF binding sites using data from NGS experiments. Availability and implementation: SNPhood is publicly available and maintained as an R Bioconductor package at http://bioconductor.org/packages/SNPhood/ . Contact: judith.zaugg@embl.de Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Versatile and efficient variant calling tools are needed to analyze large scale sequencing datasets. In particular, identification of copy number changes remains a challenging task due to their complexity, susceptibility to sequencing biases, variation in coverage data and dependence on genome-wide sample properties, such as tumor polyploidy or polyclonality in cancer samples. Results: We have developed a new tool, Canvas, for identification of copy number changes from diverse sequencing experiments including whole-genome matched tumor-normal and single-sample normal re-sequencing, as well as whole-exome matched and unmatched tumor-normal studies. In addition to variant calling, Canvas infers genome-wide parameters such as cancer ploidy, purity and heterogeneity. It provides fast and easy-to-run workflows that can scale to thousands of samples and can be easily incorporated into variant calling pipelines. Availability and Implementation: Canvas is distributed under an open source license and can be downloaded from https://github.com/Illumina/canvas . Contact: eroller@illumina.com Supplementary information: Supplementary data are available at Bioinformatics online.

Description: p ileup.js is a new browser-based genome viewer. It is designed to facilitate the investigation of evidence for genomic variants within larger web applications. It takes advantage of recent developments in the JavaScript ecosystem to provide a modular, reliable and easily embedded library. Availability and implementation: The code and documentation for pileup.js is publicly available at https://github.com/hammerlab/pileup.js under the Apache 2.0 license. Contact : correspondence@hammerlab.org

Description: Motivation: We present an update to the pathway enrichment analysis tool ‘Pathway Analysis by Randomization Incorporating Structure (PARIS)’ that determines aggregated association signals generated from genome-wide association study results. Pathway-based analyses highlight biological pathways associated with phenotypes. PARIS uses a unique permutation strategy to evaluate the genomic structure of interrogated pathways, through permutation testing of genomic features, thus eliminating many of the over-testing concerns arising with other pathway analysis approaches. Results: We have updated PARIS to incorporate expanded pathway definitions through the incorporation of new expert knowledge from multiple database sources, through customized user provided pathways, and other improvements in user flexibility and functionality. Availability and implementation: PARIS is freely available to all users at https://ritchielab.psu.edu/software/paris-download . Contact: jnc43@case.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : Nucleotide Similarity Scanner (NSimScan) is specialized for searching massive DNA databases for distant similarities. Its targeted applications include phylogenomics, comparative and functional studies of non-coding sequences, contamination detection, etc. NSimScan outperforms industry standard tools in combined sensitivity, accuracy and speed, operating at sensitivity similar to BLAST, accuracy of ssearch and speed of MegaBLAST. Availability and implementation: NSimScan is available at https://github.com/abadona/qsimscan as a part of QSimScan package. It is implemented in C ++, distributed under MIT license and supported on Linux, OS X and Windows (with cygwin). Contact: dkaznadzey@yahoo.com Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : We present TreeDom, a web tool for graphically analysing the evolutionary history of domains in multi-domain proteins. Individual domains on the same protein chain may have distinct evolutionary histories, which is important to grasp in order to understand protein function. For instance, it may be important to know whether a domain was duplicated recently or long ago, to know the origin of inserted domains, or to know the pattern of domain loss within a protein family. TreeDom uses the Pfam database as the source of domain annotations, and displays these on a sequence tree. An advantage of TreeDom is that the user can limit the analysis to N sequences that are most similar to a query, or provide a list of sequence IDs to include. Using the Pfam alignment of the selected sequences, a tree is built and displayed together with the domain architecture of each sequence. Availablility and implementation: http://TreeDom.sbc.su.se Contact: Erik.Sonnhammer@scilifelab.se

Description: : The NCI-60 human tumor cell line panel is an invaluable resource for cancer researchers, providing drug sensitivity, molecular and phenotypic data for a range of cancer types. CellMiner is a web resource that provides tools for the acquisition and analysis of quality-controlled NCI-60 data. CellMiner supports queries of up to 150 drugs or genes, but the output is an Excel file for each drug or gene. This output format makes it difficult for researchers to explore the data from large queries. CellMiner Companion is a web application that facilitates the exploration and visualization of output from CellMiner, further increasing the accessibility of NCI-60 data. Availability and Implementation: The web application is freely accessible at https://pul-bioinformatics.shinyapps.io/CellMinerCompanion . The R source code can be downloaded at https://github.com/pepascuzzi/CellMinerCompanion.git . Contact: ppascuzz@purdue.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : SMeagol is a software tool to simulate highly realistic microscopy data based on spatial systems biology models, in order to facilitate development, validation and optimization of advanced analysis methods for live cell single molecule microscopy data. Availability and implementation: SMeagol runs on Matlab R2014 and later, and uses compiled binaries in C for reaction–diffusion simulations. Documentation, source code and binaries for Mac OS, Windows and Ubuntu Linux can be downloaded from http://smeagol.sourceforge.net . Contact: johan.elf@icm.uu.se Supplementary information : Supplementary data are available at Bioinformatics online.

Description: Bacterial biofilms are implicated in a wide range of implant-based and chronic infections. These infections are often associated with adverse therapeutic outcomes, owing to the decreased antibiotic susceptibility of biofilms compared with their planktonic counterparts. This altered biofilm susceptibility has been attributed to multiple factors, including a reduced antibiotic penetration. Although several studies have addressed the role of penetration barrier in biofilm-associated drug resistance, it remains inconclusive. This study was done to elucidate antibiotic penetration through biofilms formed by Staphylococcus aureus , S. epidermidis , Escherichia coli and Klebsiella pneumoniae , using an agar disk diffusion assay. Penetration capacity of six antimicrobial drugs from different classes (β-lactams, aminoglycosides, tetracyclines, phenicols, fluoroquinolones and glycopeptides) through biofilms formed by standard strains and clinical isolates from catheter-related bloodstream infections (CRBSI) was elucidated by measuring their growth-inhibition zones in lawn cultures on Mueller–Hinton agar, following diffusion of an antibiotic from an overlying disk through their biofilm to the agar medium. Penetration of only select antimicrobials (vancomycin and chloramphenicol) was hindered through biofilms. There was considerable variation in biofilm-permeating capacity depending upon the genus, strain/CRBSI isolate and antibiotic tested. Furthermore, antibiotics failed to kill the biofilm cells independent of penetration, indicating that other factors contributed substantially to biofilm resistance.

Description: Motivation : Photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation (PAR-CLIP) is an experimental method based on next-generation sequencing for identifying the RNA interaction sites of a given protein. The method deliberately inserts T-to-C substitutions at the RNA-protein interaction sites, which provides a second layer of evidence compared with other CLIP methods. However, the experiment includes several sources of noise which cause both low-frequency errors and spurious high-frequency alterations. Therefore, rigorous statistical analysis is required in order to separate true T-to-C base changes, following cross-linking, from noise. So far, most of the existing PAR-CLIP data analysis methods focus on discarding the low-frequency errors and rely on high-frequency substitutions to report binding sites, not taking into account the possibility of high-frequency false positive substitutions. Results : Here, we introduce BMix , a new probabilistic method which explicitly accounts for the sources of noise in PAR-CLIP data and distinguishes cross-link induced T-to-C substitutions from low and high-frequency erroneous alterations. We demonstrate the superior speed and accuracy of our method compared with existing approaches on both simulated and real, publicly available human datasets. Availability and implementation : The model is freely accessible within the BMix toolbox at www.cbg.bsse.ethz.ch/software/BMix , available for Matlab and R. Supplementary information: Supplementary data is available at Bioinformatics online. Contact : niko.beerenwinkel@bsse.ethz.ch

Description: Motivation: Gene networks have become a central tool in the analysis of genomic data but are widely regarded as hard to interpret. This has motivated a great deal of comparative evaluation and research into best practices. We explore the possibility that this may lead to overfitting in the field as a whole. Results: We construct a model of ‘research communities’ sampling from real gene network data and machine learning methods to characterize performance trends. Our analysis reveals an important principle limiting the value of replication, namely that targeting it directly causes ‘easy’ or uninformative replication to dominate analyses. We find that when sampling across network data and algorithms with similar variability, the relationship between replicability and accuracy is positive (Spearman’s correlation, r s ~0.33) but where no such constraint is imposed, the relationship becomes negative for a given gene function ( r s ~ –0.13). We predict factors driving replicability in some prior analyses of gene networks and show that they are unconnected with the correctness of the original result, instead reflecting replicable biases. Without these biases, the original results also vanish replicably. We show these effects can occur quite far upstream in network data and that there is a strong tendency within protein–protein interaction data for highly replicable interactions to be associated with poor quality control. Availability and implementation: Algorithms, network data and a guide to the code available at: https://github.com/wimverleyen/AggregateGeneFunctionPrediction . Contact: jgillis@cshl.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Recent advances in single molecule real-time (SMRT) and nanopore sequencing technologies have enabled high-quality assemblies from long and inaccurate reads. However, these approaches require high coverage by long reads and remain expensive. On the other hand, the inexpensive short reads technologies produce accurate but fragmented assemblies. Thus, a hybrid approach that assembles long reads (with low coverage) and short reads has a potential to generate high-quality assemblies at reduced cost. Results: We describe hybrid SPA des algorithm for assembling short and long reads and benchmark it on a variety of bacterial assembly projects. Our results demonstrate that hybrid SPA des generates accurate assemblies (even in projects with relatively low coverage by long reads) thus reducing the overall cost of genome sequencing. We further present the first complete assembly of a genome from single cells using SMRT reads. Availability and implementation: hybrid SPA des is implemented in C++ as a part of SPAdes genome assembler and is publicly available at http://bioinf.spbau.ru/en/spades Contact: d.antipov@spbu.ru Supplementary information: supplementary data are available at Bioinformatics online.

Description: Motivation: There are numerous examples of RNA–RNA complexes, including microRNA–mRNA and small RNA–mRNA duplexes for regulation of translation, guide RNA interactions with target RNA for post-transcriptional modification and small nuclear RNA duplexes for splicing. Predicting the base pairs formed between two interacting sequences remains difficult, at least in part because of the competition between unimolecular and bimolecular structure. Results: Two algorithms were developed for improved prediction of bimolecular RNA structure that consider the competition between self-structure and bimolecular structure. These algorithms utilize two novel approaches to evaluate accessibility: free energy density minimization and pseudo-energy minimization. Free energy density minimization minimizes the folding free energy change per nucleotide involved in an intermolecular secondary structure. Pseudo-energy minimization (called AccessFold) minimizes the sum of free energy change and a pseudo-free energy penalty for bimolecular pairing of nucleotides that are unlikely to be accessible for bimolecular structure. The pseudo-free energy, derived from unimolecular pairing probabilities, is applied per nucleotide in bimolecular pairs, and this approach is able to predict binding sites that are split by unimolecular structures. A benchmark set of 17 bimolecular RNA structures was assembled to assess structure prediction. Pseudo-energy minimization provides a statistically significant improvement in sensitivity over the method that was found in a benchmark to be the most accurate previously available method, with an improvement from 36.8% to 57.8% in mean sensitivity for base pair prediction. Availability and implementation: Pseudo-energy minimization is available for download as AccessFold, under an open-source license and as part of the RNAstructure package, at: http://rna.urmc.rochester.edu/RNAstructure.html . Contact: david_mathews@urmc.rochester.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Simulating complex evolution scenarios of multiple populations is an important task for answering many basic questions relating to population genomics. Apart from the population samples, the underlying Ancestral Recombinations Graph (ARG) is an additional important means in hypothesis checking and reconstruction studies. Furthermore, complex simulations require a plethora of interdependent parameters making even the scenario-specification highly non-trivial. Results: We present an algorithm SimRA that simulates generic multiple population evolution model with admixture. It is based on random graphs that improve dramatically in time and space requirements of the classical algorithm of single populations. Using the underlying random graphs model, we also derive closed forms of expected values of the ARG characteristics i.e., height of the graph, number of recombinations, number of mutations and population diversity in terms of its defining parameters. This is crucial in aiding the user to specify meaningful parameters for the complex scenario simulations, not through trial-and-error based on raw compute power but intelligent parameter estimation. To the best of our knowledge this is the first time closed form expressions have been computed for the ARG properties. We show that the expected values closely match the empirical values through simulations. Finally, we demonstrate that SimRA produces the ARG in compact forms without compromising any accuracy. We demonstrate the compactness and accuracy through extensive experiments. Availability and implementation : SimRA ( Sim ulation based on R andom graph A lgorithms) source, executable, user manual and sample input-output sets are available for downloading at: https://github.com/ComputationalGenomics/SimRA Contact : parida@us.ibm.com Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation : The Optical Mapping System discovers structural variants and potentiates sequence assembly of genomes via scaffolding and comparisons that globally validate or correct sequence assemblies. Despite its utility, there are few publicly available tools for aligning optical mapping datasets. Results : Here we present software, named ‘Maligner’, for the alignment of both single molecule restriction maps (Rmaps) and in silico restriction maps of sequence contigs to a reference. Maligner provides two modes of alignment: an efficient, sensitive dynamic programming implementation that scales to large eukaryotic genomes, and a faster indexed based implementation for finding alignments with unmatched sites in the reference but not the query. We compare our software to other publicly available tools on Rmap datasets and show that Maligner finds more correct alignments in comparable runtime. Lastly, we introduce the M-Score statistic for normalizing alignment scores across restriction maps and demonstrate its utility for selecting high quality alignments. Availability and implementation : The Maligner software is written in C ++ and is available at https://github.com/LeeMendelowitz/maligner under the GNU General Public License. Contact : mpop@umiacs.umd.edu

Description: Motivation: The detection of genomic structural variation (SV) has advanced tremendously in recent years due to progress in high-throughput sequencing technologies. Novel sequence insertions, insertions without similarity to a human reference genome, have received less attention than other types of SVs due to the computational challenges in their detection from short read sequencing data, which inherently involves de novo assembly. De novo assembly is not only computationally challenging, but also requires high-quality data. Although the reads from a single individual may not always meet this requirement, using reads from multiple individuals can increase power to detect novel insertions. Results: We have developed the program PopIns , which can discover and characterize non-reference insertions of 100 bp or longer on a population scale. In this article, we describe the approach we implemented in PopIns. It takes as input a reads-to-reference alignment, assembles unaligned reads using a standard assembly tool, merges the contigs of different individuals into high-confidence sequences, anchors the merged sequences into the reference genome, and finally genotypes all individuals for the discovered insertions. Our tests on simulated data indicate that the merging step greatly improves the quality and reliability of predicted insertions and that PopIns shows significantly better recall and precision than the recent tool MindTheGap. Preliminary results on a dataset of 305 Icelanders demonstrate the practicality of the new approach. Availability and implementation: The source code of PopIns is available from http://github.com/bkehr/popins . Contact: birte.kehr@decode.is Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: High-throughput sequencing technologies provide access to an increasing number of bacterial genomes. Today, many analyses involve the comparison of biological properties among many strains of a given species, or among species of a particular genus. Tools that can help the microbiologist with these tasks become increasingly important. Results: Insyght is a comparative visualization tool whose core features combine a synchronized navigation across genomic data of multiple organisms with a versatile interoperability between complementary views. In this work, we have greatly increased the scope of the Insyght public dataset by including 2688 complete bacterial genomes available in Ensembl thus vastly improving its phylogenetic coverage. We also report the development of a virtual machine that allows users to easily set up and customize their own local Insyght server. Availability and implementation: http://genome.jouy.inra.fr/Insyght Contact: Thomas.Lacroix@jouy.inra.fr

Description: : Breast cancer is one of the most frequent cancers among women. Extensive studies into the molecular heterogeneity of breast cancer have produced a plethora of molecular subtype classification and prognosis prediction algorithms, as well as numerous gene expression signatures. However, reimplementation of these algorithms is a tedious but important task to enable comparison of existing signatures and classification models between each other and with new models. Here, we present the genefu R/Bioconductor package, a multi-tiered compendium of bioinformatics algorithms and gene signatures for molecular subtyping and prognostication in breast cancer. Availability and implementation: The genefu package is available from Bioconductor. http://www.bioconductor.org/packages/devel/bioc/html/genefu.html . Source code is also available on Github https://github.com/bhklab/genefu . Contact: bhaibeka@uhnresearch.ca Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Aberrant glycosylation is highly associated with cancer progression. The aim of this study was to compare bifucosylated N -glycans in sera obtained from healthy controls and breast cancer patients, with the goal of identifying a potential indicator for monitoring the recurrence and metastasis of breast cancer. A unique structural pattern of bifucosylated N -glycan, with both core and antennary fucosylation, was identified in breast cancer patients. The spectrum of antennary fucosylation was a composite of the standard spectra of Lewis X and H2, indicating a mixture of the two epitopes. Permethylated N -glycans of the glycoproteins extracted from 91 breast cancer patients and 43 healthy controls were detected using linear ion-trap quadrupole-electrospray ionization mass spectrometry, which appeared to be a highly sensitive and useful approach in the detection and identification of N -glycans. To evaluate MS profile data, several statistical tools were applied, including Student's t -test, partial least squares discriminant analysis and receiver-operating characteristic curve. The results showed that the measurement of bifucosylation degree and CEA levels had an improved diagnostic performance compared with that of CEA alone. We compared the potential of bifucosylated N -glycan as an indicator of breast cancer recurrence with the current clinical biomarkers, i.e., CEA, CA 15-3 and CA125. The result revealed that, compared with CEA, CA 15-3 and CA125, the bifucosylation degree of N -glycans could be a more reliable indicator of breast cancer recurrence.

Description: Colitose, also known as 3,6-dideoxy- l -galactose or 3-deoxy- l -fucose, is one of only five naturally occurring 3,6-dideoxyhexoses. Colitose was found in lipopolysaccharide of a number of infectious bacteria, including Escherichia coli O55 & O111 and Vibrio cholera O22 & O139. To date, no colitosyltransferase (ColT) has been characterized, probably due to the inaccessibility of the sugar donor, GDP-colitose. In this study, starting with chemically prepared colitose, 94.6 mg of GDP-colitose was prepared via a facile and efficient one-pot two-enzyme system involving an l -fucokinase/GDP- l -Fuc pyrophosphorylase and an inorganic pyrophosphatase (EcPpA). WbgN, a putative ColT from E. coli O55:H5 was then cloned, overexpressed, purified and biochemically characterized by using GDP-colitose as a sugar donor. Activity assay and structural identification of the synthetic product clearly demonstrated that wbgN encodes an α1,2-ColT. Biophysical study showed that WbgN does not require metal ion, and is highly active at pH 7.5–9.0. In addition, acceptor specificity study indicated that WbgN exclusively recognizes lacto- N -biose (Galβ1,3-GlcNAc). Most interestingly, it was found that WbgN exhibits similar activity toward GDP- l -Fuc ( k cat / K m = 9.2 min –1 mM –1 ) as that toward GDP-colitose ( k cat / K m = 12 min –1 mM –1 ). Finally, taking advantage of this, type 1 H-antigen was successfully synthesized in preparative scale.

Description: Information about specificity and affinity is critical for use of carbohydrate-binding antibodies. Herein, we evaluated eight monoclonal antibodies to the blood group A (BG-A) antigen. Antibodies 87-G, 9A, HE-10, HE-24, HE-193, HE-195, T36 and Z2A were profiled on a glycan microarray to assess specificity, relative affinity and the influence of glycan density on recognition. Our studies highlight several noteworthy recognition properties. First, most antibodies bound GalNAcα1–3Gal and the BG-A trisaccharide nearly as well as larger BG-A oligosaccharides. Second, several antibodies only bound the BG-A trisaccharide when displayed on certain glycan chains. These first two points indicate that the carrier glycan chains primarily influence selectivity, rather than binding strength. Third, binding of some antibodies was highly dependent on glycan density, illustrating the importance of glycan presentation for recognition. Fourth, some antibodies recognized the tumor-associated Tn antigen, and one antibody only bound the variant composed of a GalNAc-alpha-linked to a serine residue. Collectively, these results provide new insights into the recognition properties of anti-BG-A antibodies.

Description: New and effective tuberculosis (TB) vaccines are urgently needed to control pulmonary TB, and in particular to prevent the spread of drug-resistant strains of Mycobacterium tuberculosis . These drug-resistant strains can range from those resistant to first-line drugs to those that are almost impossible to treat. To develop new and effective vaccines for HIV and malaria has been difficult and it is proving to be just as challenging for TB. TB is a complicated disease with a spectrum from apparently controlled latent infection to active clinical disease and so different types of preventive or post-exposure vaccine may be needed. Identifying the most promising vaccine candidates to move into clinical trials is difficult, as we lack biomarker signatures that can predict protective efficacy. There is a risk that the failure of the MVA-85A vaccine to show efficacy when given to previously BCG-vaccinated South African infants will impact on the resources available for the development and trials of other candidate TB vaccines. Continued support for the development of new TB vaccines should remain a priority as an effective vaccine would bring huge public health benefits.

Description: : Precise regulatory control of genes, particularly in eukaryotes, frequently requires the joint action of multiple sequence-specific transcription factors. A cis -regulatory module (CRM) is a genomic locus that is responsible for gene regulation and that contains multiple transcription factor binding sites in close proximity. Given a collection of known transcription factor binding motifs, many bioinformatics methods have been proposed over the past 15 years for identifying within a genomic sequence candidate CRMs consisting of clusters of those motifs. Results: The MCAST algorithm uses a hidden Markov model with a P -value-based scoring scheme to identify candidate CRMs. Here, we introduce a new version of MCAST that offers improved graphical output, a dynamic background model, statistical confidence estimates based on false discovery rate estimation and, most significantly, the ability to predict CRMs while taking into account epigenomic data such as DNase I sensitivity or histone modification data. We demonstrate the validity of MCAST’s statistical confidence estimates and the utility of epigenomic priors in identifying CRMs. Availability and implementation: MCAST is part of the MEME Suite software toolkit. A web server and source code are available at http://meme-suite.org and http://alternate.meme-suite.org . Contact: t.bailey@imb.uq.edu.au or william-noble@uw.edu Supplementary information : Supplementary data are available at Bioinformatics online.

Description: : Pharmacogenomics holds great promise for the development of biomarkers of drug response and the design of new therapeutic options, which are key challenges in precision medicine. However, such data are scattered and lack standards for efficient access and analysis, consequently preventing the realization of the full potential of pharmacogenomics. To address these issues, we implemented PharmacoGx , an easy-to-use, open source package for integrative analysis of multiple pharmacogenomic datasets. We demonstrate the utility of our package in comparing large drug sensitivity datasets, such as the Genomics of Drug Sensitivity in Cancer and the Cancer Cell Line Encyclopedia. Moreover, we show how to use our package to easily perform Connectivity Map analysis. With increasing availability of drug-related data, our package will open new avenues of research for meta-analysis of pharmacogenomic data. Availability and implementation : PharmacoGx is implemented in R and can be easily installed on any system. The package is available from CRAN and its source code is available from GitHub. Contact : bhaibeka@uhnresearch.ca or benjamin.haibe.kains@utoronto.ca Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : In population genetics and phylogeography, haplotype genealogy graphs are important tools for the visualization of population structure based on sequence data. In this type of graph, node sizes are often drawn in proportion to haplotype frequencies and edge lengths represent the minimum number of mutations separating adjacent nodes. I here present Fitchi, a new program that produces publication-ready haplotype genealogy graphs based on the Fitch algorithm. Availability and implementation: http://www.evoinformatics.eu/fitchi.htm Contact : michaelmatschiner@mac.com Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : Genome-wide association studies (GWASs) have successfully identified many sequence variants that are significantly associated with common diseases and traits. Tens of thousands of such trait-associated SNPs have already been cataloged, which we believe form a great resource for genomic research. Recent studies have demonstrated that the collection of trait-associated SNPs can be exploited to indicate whether a given genomic interval or intervals are likely to be functionally connected with certain phenotypes or diseases. Despite this importance, currently, there is no ready-to-use computational tool able to connect genomic intervals to phenotypes. Here, we present traseR , an easy-to-use R Bioconductor package that performs enrichment analyses of trait-associated SNPs in arbitrary genomic intervals with flexible options, including testing method, type of background and inclusion of SNPs in LD. Availability and implementation: The traseR R package preloaded with up-to-date collection of trait-associated SNPs are freely available in Bioconductor Contact: zhaohui.qin@emory.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : Data-dependent acquisition (DDA) is the most common method used to control the acquisition process of shotgun proteomics experiments. While novel DDA approaches have been proposed, their evaluation is made difficult by the need of programmatic control of a mass spectrometer. An alternative is in silico analysis, for which suitable software has been unavailable. To meet this need, we have developed MSAcquisitionSimulator—a collection of C ++ programs for simulating ground truth LC-MS data and the subsequent application of custom DDA algorithms. It provides an opportunity for researchers to test, refine and evaluate novel DDA algorithms prior to implementation on a mass spectrometer. Availability and implementation: The software is freely available from its Github repository http://www.github.com/DennisGoldfarb/MSAcquisitionSimulator/ which contains further documentation and usage instructions. Contact: weiwang@cs.ucla.edu or ben_major@med.unc.edu Supplementary information : Supplementary data are available at Bioinformatics online.

Description: : Coarse grain (CG) models allow long-scale simulations with a much lower computational cost than that of all-atom simulations. However, the absence of atomistic detail impedes the analysis of specific atomic interactions that are determinant in most interesting biomolecular processes. In order to study these phenomena, it is necessary to reconstruct the atomistic structure from the CG representation. This structure can be analyzed by itself or be used as an onset for atomistic molecular dynamics simulations. In this work, we present a computer program that accurately reconstructs the atomistic structure from a CG model for proteins, using a simple geometrical algorithm. Availability and implementation: The software is free and available online at http://www.ic.fcen.uba.ar/cg2aa/cg2aa.py Supplementary information: Supplementary data are available at Bioinformatics online. Contact: lula@qi.fcen.uba.ar

Description: : The Sun Grid Engine (SGE) high-performance computing batch queueing system is commonly used in bioinformatics analysis. Creating re-usable scripts for the SGE is a common challenge. The qsubsec template language and interpreter described here allow researchers to easily create generic template definitions that encapsulate a particular computational job, effectively separating the process logic from the specific run details. At submission time, the generic template is filled in with specific values. This system provides an intermediate level between simple scripting and complete workflow management tools. Availability and implementation: Qsubsec is open-source and is available at https://github.com/alastair-droop/qsubsec . Contact: a.p.droop@leeds.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Description: : The lack of visualization frameworks to guide interpretation and facilitate discovery is a potential bottleneck for precision medicine, systems genetics and other studies. To address this we have developed an interactive, reproducible, web-based prioritization approach that builds on our earlier work. HitWalker2 is highly flexible and can utilize many data types and prioritization methods based upon available data and desired questions, allowing it to be utilized in a diverse range of studies such as cancer, infectious disease and psychiatric disorders. Availability and implementation: Source code is freely available at https://github.com/biodev/HitWalker2 and implemented using Python/Django, Neo4j and Javascript (D3.js and jQuery). We support major open source browsers (e.g. Firefox and Chromium/Chrome). Contact: wilmotb@ohsu.edu Supplementary information: Supplementary data are available at Bioinformatics online. Additional information/instructions are available at https://github.com/biodev/HitWalker2/wiki

Description: Motivation: Whole genome sequencing (WGS) of parent-offspring trios is a powerful approach for identifying disease-associated genes via detecting copy number variations (CNVs). Existing approaches, which detect CNVs for each individual in a trio independently, usually yield low-detection accuracy. Joint modeling approaches leveraging Mendelian transmission within the parent-offspring trio can be an efficient strategy to improve CNV detection accuracy. Results: In this study, we developed TrioCNV, a novel approach for jointly detecting CNVs in parent-offspring trios from WGS data. Using negative binomial regression, we modeled the read depth signal while considering both GC content bias and mappability bias. Moreover, we incorporated the family relationship and used a hidden Markov model to jointly infer CNVs for three samples of a parent-offspring trio. Through application to both simulated data and a trio from 1000 Genomes Project, we showed that TrioCNV achieved superior performance than existing approaches. Availability and implementation: The software TrioCNV implemented using a combination of Java and R is freely available from the website at https://github.com/yongzhuang/TrioCNV . Contact: ydwang@hit.edu.cn Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: Recent advancements in molecular methods have made it possible to capture physical contacts between multiple chromatin fragments. The resulting association matrices provide a noisy estimate for average spatial proximity that can be used to gain insights into the genome organization inside the nucleus. However, extracting topological information from these data is challenging and their integration across resolutions is still poorly addressed. Recent findings suggest that a hierarchical approach could be advantageous for addressing these challenges. Results: We present an algorithmic framework, which is based on hierarchical block matrices (HBMs), for topological analysis and integration of chromosome conformation capture (3C) data. We first describe chromoHBM, an algorithm that compresses high-throughput 3C (HiT-3C) data into topological features that are efficiently summarized with an HBM representation. We suggest that instead of directly combining HiT-3C datasets across resolutions, which is a difficult task, we can integrate their HBM representations, and describe chromoHBM-3C, an algorithm which merges HBMs. Since three-dimensional (3D) reconstruction can also benefit from topological information, we further present chromoHBM-3D, an algorithm which exploits the HBM representation in order to gradually introduce topological constraints to the reconstruction process. We evaluate our approach in light of previous image microscopy findings and epigenetic data, and show that it can relate multiple spatial scales and provide a more complete view of the 3D genome architecture. Availability and implementation: The presented algorithms are available from: https://github.com/yolish/hbm . Contact: ys388@cam.ac.uk or pl219@cam.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Motivation: There are various reasons for rerunning bioinformatics tools and pipelines on sequencing data, including reproducing a past result, validation of a new tool or workflow using a known dataset, or tracking the impact of database changes. For identical results to be achieved, regularly updated reference sequence databases must be versioned and archived. Database administrators have tried to fill the requirements by supplying users with one-off versions of databases, but these are time consuming to set up and are inconsistent across resources. Disk storage and data backup performance has also discouraged maintaining multiple versions of databases since databases such as NCBI nr can consume 50 Gb or more disk space per version, with growth rates that parallel Moore's law. Results: Our end-to-end solution combines our own Kipper software package—a simple key-value large file versioning system—with BioMAJ (software for downloading sequence databases), and Galaxy (a web-based bioinformatics data processing platform). Available versions of databases can be recalled and used by command-line and Galaxy users. The Kipper data store format makes publishing curated FASTA databases convenient since in most cases it can store a range of versions into a file marginally larger than the size of the latest version. Availability and implementation: Kipper v1.0.0 and the Galaxy Versioned Data tool are written in Python and released as free and open source software available at https://github.com/Public-Health-Bioinformatics/kipper and https://github.com/Public-Health-Bioinformatics/versioned_data , respectively; detailed setup instructions can be found at https://github.com/Public-Health-Bioinformatics/versioned_data/blob/master/doc/setup.md Contact : Damion.Dooley@Bccdc.Ca or William.Hsiao@Bccdc.Ca Supplementary information: Supplementary data are available at Bioinformatics online.