JBrowse Jupyter: a Python interface to JBrowse 2.

MOTIVATION: JBrowse Jupyter is a package that aims to close the gap between Python programming and genomic visualization. Web-based genome browsers are routinely used for publishing and inspecting genome annotations. Historically they have been deployed at the end of bioinformatics pipelines, typically decoupled from the analysis itself. However, emerging technologies such as Jupyter notebooks enable a more rapid iterative cycle of development, analysis and visualization. RESULTS: We have developed a package that provides a Python interface to JBrowse 2's suite of embeddable components, including the primary Linear Genome View. The package enables users to quickly set up, launch and customize JBrowse views from Jupyter notebooks. In addition, users can share their data via Google's Colab notebooks, providing reproducible interactive views. AVAILABILITY AND IMPLEMENTATION: JBrowse Jupyter is released under the Apache License and is available for download on PyPI. Source code and demos are available on GitHub at https://github.com/GMOD/jbrowse-jupyter.

Apollo: Democratizing genome annotation.

Genome annotation is the process of identifying the location and function of a genome's encoded features. Improving the biological accuracy of annotation is a complex and iterative process requiring researchers to review and incorporate multiple sources of information such as transcriptome alignments, predictive models based on sequence profiles, and comparisons to features found in related organisms. Because rapidly decreasing costs are enabling an ever-growing number of scientists to incorporate sequencing as a routine laboratory technique, there is widespread demand for tools that can assist in the deliberative analytical review of genomic information. To this end, we present Apollo, an open source software package that enables researchers to efficiently inspect and refine the precise structure and role of genomic features in a graphical browser-based platform. Some of Apollo's newer user interface features include support for real-time collaboration, allowing distributed users to simultaneously edit the same encoded features while also instantly seeing the updates made by other researchers on the same region in a manner similar to Google Docs. Its technical architecture enables Apollo to be integrated into multiple existing genomic analysis pipelines and heterogeneous laboratory workflow platforms. Finally, we consider the implications that Apollo and related applications may have on how the results of genome research are published and made accessible.

Historian: accurate reconstruction of ancestral sequences and evolutionary rates.

Motivation: Reconstruction of ancestral sequence histories, and estimation of parameters like indel rates, are improved by using explicit evolutionary models and summing over uncertain alignments. The previous best tool for this purpose (according to simulation benchmarks) was ProtPal, but this tool was too slow for practical use. Results: Historian combines an efficient reimplementation of the ProtPal algorithm with performance-improving heuristics from other alignment tools. Simulation results on fidelity of rate estimation via ancestral reconstruction, along with evaluations on the structurally informed alignment dataset BAliBase 3.0, recommend Historian over other alignment tools for evolutionary applications. Availability and Implementation: Historian is available at https://github.com/evoldoers/historian under the Creative Commons Attribution 3.0 US license. Contact: ihholmes+historian@gmail.com.

BioMake: a GNU make-compatible utility for declarative workflow management.

MOTIVATION: The Unix 'make' program is widely used in bioinformatics pipelines, but suffers from problems that limit its application to large analysis datasets. These include reliance on file modification times to determine whether a target is stale, lack of support for parallel execution on clusters, and restricted flexibility to extend the underlying logic program. RESULTS: We present BioMake, a make-like utility that is compatible with most features of GNU Make and adds support for popular cluster-based job-queue engines, MD5 signatures as an alternative to timestamps, and logic programming extensions in Prolog. AVAILABILITY AND IMPLEMENTATION: BioMake is available for MacOSX and Linux systems from https://github.com/evoldoers/biomake under the BSD3 license. The only dependency is SWI-Prolog (version 7), available from http://www.swi-prolog.org/. CONTACT: ihholmes + biomake@gmail.com or cmungall + biomake@gmail.com. SUPPLEMENTARY INFORMATION: Feature table comparing BioMake to similar tools. Supplementary data are available at Bioinformatics online.

Solving the master equation for Indels.

BACKGROUND: Despite the long-anticipated possibility of putting sequence alignment on the same footing as statistical phylogenetics, theorists have struggled to develop time-dependent evolutionary models for indels that are as tractable as the analogous models for substitution events. MAIN TEXT: This paper discusses progress in the area of insertion-deletion models, in view of recent work by Ezawa (BMC Bioinformatics 17:304, 2016); (BMC Bioinformatics 17:397, 2016); (BMC Bioinformatics 17:457, 2016) on the calculation of time-dependent gap length distributions in pairwise alignments, and current approaches for extending these approaches from ancestor-descendant pairs to phylogenetic trees. CONCLUSIONS: While approximations that use finite-state machines (Pair HMMs and transducers) currently represent the most practical approach to problems such as sequence alignment and phylogeny, more rigorous approaches that work directly with the matrix exponential of the underlying continuous-time Markov chain also show promise, especially in view of recent advances.

RNAcentral: A vision for an international database of RNA sequences.

During the last decade there has been a great increase in the number of noncoding RNA genes identified, including new classes such as microRNAs and piRNAs. There is also a large growth in the amount of experimental characterization of these RNA components. Despite this growth in information, it is still difficult for researchers to access RNA data, because key data resources for noncoding RNAs have not yet been created. The most pressing omission is the lack of a comprehensive RNA sequence database, much like UniProt, which provides a comprehensive set of protein knowledge. In this article we propose the creation of a new open public resource that we term RNAcentral, which will contain a comprehensive collection of RNA sequences and fill an important gap in the provision of biomedical databases. We envision RNA researchers from all over the world joining a federated RNAcentral network, contributing specialized knowledge and databases. RNAcentral would centralize key data that are currently held across a variety of databases, allowing researchers instant access to a single, unified resource. This resource would facilitate the next generation of RNA research and help drive further discoveries, including those that improve food production and human and animal health. We encourage additional RNA database resources and research groups to join this effort. We aim to obtain international network funding to further this endeavor.

JBrowse 2: a modular genome browser with views of synteny and structural variation.

We present JBrowse 2, a general-purpose genome annotation browser offering enhanced visualization of complex structural variation and evolutionary relationships. It retains core features of JBrowse while adding new views for synteny, dotplots, breakpoints, gene fusions, and whole-genome overviews. It allows users to share sessions, open multiple genomes, and navigate between views. It can be embedded in a web page, used as a standalone application, or run from Jupyter notebooks or R sessions. These improvements are enabled by a ground-up redesign using modern web technology. We describe application functionality, use cases, performance benchmarks, and implementation notes for web administrators and developers.

JBrowse: a next-generation genome browser.

We describe an open source, portable, JavaScript-based genome browser, JBrowse, that can be used to navigate genome annotations over the web. JBrowse helps preserve the user's sense of location by avoiding discontinuous transitions, instead offering smoothly animated panning, zooming, navigation, and track selection. Unlike most existing genome browsers, where the genome is rendered into images on the webserver and the role of the client is restricted to displaying those images, JBrowse distributes work between the server and client and therefore uses significantly less server overhead than previous genome browsers. We report benchmark results empirically comparing server- and client-side rendering strategies, review the architecture and design considerations of JBrowse, and describe a simple wiki plug-in that allows users to upload and share annotation tracks.

Colorstock, SScolor, Ratón: RNA alignment visualization tools.

UNLABELLED: Interactive examination of RNA multiple alignments for covariant mutations is a useful step in non-coding RNA sequence analysis. We present three parallel implementations of an RNA visualization metaphor: Colorstock, a command-line script using ANSI terminal color; SScolor, a Perl script that generates static HTML pages; and Ratón, an AJAX web application generating dynamic HTML. Each tool can be used to color RNA alignments by secondary structure and to visually highlight compensatory mutations in stems. AVAILABILITY: All source code is freely available under the GPL. The source code can be downloaded and a prototype of Ratón can be accessed at http://biowiki.org/RnaAlignmentViewers.

JBrowse: a dynamic web platform for genome visualization and analysis.

BACKGROUND: JBrowse is a fast and full-featured genome browser built with JavaScript and HTML5. It is easily embedded into websites or apps but can also be served as a standalone web page. RESULTS: Overall improvements to speed and scalability are accompanied by specific enhancements that support complex interactive queries on large track sets. Analysis functions can readily be added using the plugin framework; most visual aspects of tracks can also be customized, along with clicks, mouseovers, menus, and popup boxes. JBrowse can also be used to browse local annotation files offline and to generate high-resolution figures for publication. CONCLUSIONS: JBrowse is a mature web application suitable for genome visualization and analysis.

Identification of crucial residues of conformational epitopes on VP2 protein of infectious bursal disease virus by phage display.

A new method for identifying epitopes in viral proteins expressed by filamentous phage has been developed. Filamentous phage fUSE 1 containing the variable region of the VP2 gene of infectious bursal disease virus (IBDV) strain 002-73 was constructed. Neutralizing monoclonal antibodies 17-82 and 33-10 raised against VP2 protein were used to bind phage containing the original variable region of VP2. The phage bound to monoclonal antibodies, were removed by protein G Sepharose and the unbound phage (escape mutants) were isolated for sequencing to locate the mutations. The crucial amino acid residues for conformational neutralizing epitopes recognized by the monoclonal antibodies were located in the first main hydrophilic region (amino acids from 210 to 225) and the central region of the variable region of VP2. The amino acid residues on both ends of the variable region of VP2 affected considerably the binding of monoclonal antibodies. This technique might be useful for selecting escape mutants of phage displaying the original antigenic regions of other viruses to define the crucial amino acid residues of their conformational epitopes, especially viruses that cannot be grown in cell cultures.

Mapping of conformational epitopes on capsid protein VP2 of infectious bursal disease virus by fd-tet phage display.

Conformational epitopes on VP2 protein of infectious bursal disease virus (IBDV) were mapped using fd-tet phage display. A gene-targeted phage display library was made using DNA fragments ranging approximately from 80 to 400 bp of the hypervariable region of the VP2 gene of IBDV strain 002-73, as neutralizing monoclonal antibodies against the VP2 protein recognize VP2 conformation-dependent epitopes within the hypervariable region. The phages were selected using immobilized monoclonal antibodies. Epitopes on five phages selected with monoclonal antibody 17-82 were located between amino acids 211 and 344. A constructed phage containing amino acids from 204 to 344 strongly reacted with monoclonal antibodies. Compared to that of the constructed phage, the binding of monoclonal antibodies to the five selected phages was dramatically reduced when several amino acids at either terminus or both termini were absent. The binding of a phage, with conversion of the first hydrophilic region into a hydrophobic region as a result of a chance frameshift mutation from amino acids 214 to 225, dropped sharply. It indicates that conformational epitopes may be up to 423 bp long and the commonly suggested fragments of 50-300 bp for making gene-targeted phage display libraries are not long enough to cover the conformational epitopes. This technique can be used to identify the minimum length of the conformational epitopes for developing recombinant vaccines and specific diagnostic tests.

Web Apollo: a web-based genomic annotation editing platform.

Web Apollo is the first instantaneous, collaborative genomic annotation editor available on the web. One of the natural consequences following from current advances in sequencing technology is that there are more and more researchers sequencing new genomes. These researchers require tools to describe the functional features of their newly sequenced genomes. With Web Apollo researchers can use any of the common browsers (for example, Chrome or Firefox) to jointly analyze and precisely describe the features of a genome in real time, whether they are in the same room or working from opposite sides of the world.

Setting up the JBrowse genome browser.

JBrowse is a Web-based tool for visualizing genomic data. Unlike most other Web-base genome browsers, JBrowse exploits the capabilities of the user's Web browser to make scrolling and zooming fast and smooth. It supports the browsers used by almost all Internet users, and is relatively simple to install. JBrowse can utilize multiple types of data in a variety of common genomic data formats, including genomic feature data in bioperl databases, GFF files, BED files, and quantitative data in wiggle files. This unit describes how to obtain the JBrowse software, set it up on a Linux or Mac OS X computer running as a Web server, and incorporate genome annotation data from multiple sources into JBrowse. After completing the protocols described in this unit, the reader will have a Web site that other users can visit to browse the genomic data.

Tools for simulating evolution of aligned genomic regions with integrated parameter estimation.

Controlled simulations of genome evolution are useful for benchmarking tools. However, many simulators lack extensibility and cannot measure parameters directly from data. These issues are addressed by three new open-source programs: GSIMULATOR (for neutrally evolving DNA), SIMGRAM (for generic structured features) and SIMGENOME (for syntenic genome blocks). Each offers algorithms for parameter measurement and reconstruction of ancestral sequence. All three tools out-perform the leading neutral DNA simulator (DAWG) in benchmarks. The programs are available at http://biowiki.org/SimulationTools.

Rotaviruses interact with alpha4beta7 and alpha4beta1 integrins by binding the same integrin domains as natural ligands.

Group A rotaviruses are major intestinal pathogens that express potential alpha4beta1 and alpha4beta7 integrin ligand sequences Leu-Asp-Val and Leu-Asp-Ile in their outer capsid protein VP7, and Ile-Asp-Ala in their spike protein VP4. Monkey rotavirus SA11 can use recombinant alpha4beta1 as a cellular receptor. In this study a new potential alpha4beta1, alpha4beta7 and alpha9beta1 integrin ligand sequence, Tyr-Gly-Leu, was identified in VP4. It was shown that several human and monkey rotaviruses bound alpha4beta1 and alpha4beta7, but not alpha9beta1. Binding to alpha4beta1 mediated the infectivity and growth of monkey rotaviruses, and binding to alpha4beta7 mediated their infectivity. A porcine rotavirus interacted with alpha4 integrins at a post-binding stage to facilitate infection. Activation of alpha4beta1 increased rotavirus infectivity. Cellular treatment with peptides containing the alpha4 integrin ligand sequences Tyr-Gly-Leu and Ile-Asp-Ala eliminated virus binding to alpha4 integrins and infectivity. In contrast, rotavirus recognition of alpha4 integrins was unaffected by a peptide containing the sequence Leu-Asp-Val or by a mutation in the VP7 Leu-Asp-Val sequence. VP4 involvement in rotavirus recognition of alpha4beta1 was demonstrated with rotavirus reassortants. Swapping and point mutagenesis of alpha4 surface loops showed that rotaviruses required the same alpha4 residues and domains for binding as the natural alpha4 integrin ligands: mucosal addressin cell adhesion molecule-1, fibronectin and vascular cell adhesion molecule-1. Several rotaviruses are able to use alpha4beta7 and alpha4beta1 for cell binding or entry, through the recognition of the same alpha4-subunit domains as natural alpha4 ligands.