Sensitive and error-tolerant annotation of protein-coding DNA with BATH

We present BATH, a tool for highly sensitive annotation of protein-coding DNA based on direct alignment of that DNA to a database of protein sequences or profile hidden Markov models (pHMMs). BATH is built on top of the HMMER3 code base, and simplifies the annotation workflow for pHMM-based annotation by providing a straightforward input interface and easy-to-interpret output. BATH also introduces novel frameshift-aware algorithms to detect frameshift-inducing nucleotide insertions and deletions (indels). BATH matches the accuracy of HMMER3 for annotation of sequences containing no errors, and produces superior accuracy to all tested tools for annotation of sequences containing nucleotide indels. These results suggest that BATH should be used when high annotation sensitivity is required, particularly when frameshift errors are expected to interrupt protein-coding regions, as is true with long read sequencing data and in the context of pseudogenes.

IsoSplitter: identification and characterization of alternative splicing sites without a reference genome

Long-read transcriptome sequencing is designed to sequence full-length RNA molecules and advantageous for identifying alternative splice isoforms; however, in the absence of a reference genome, it is difficult to accurately locate splice sites, because of the diversity of patterns of alternative splicing (AS). Based on long-read transcriptome data we developed a versatile tool, IsoSplitter, to reverse-trace and validate AS gene "split-sites" with the following features: (1) IsoSplitter initially invokes a modified SIM4 program to find transcript split-sites; (2) each split-site is then quantified, to reveal transcript diversity, and putative isoforms are grouped into gene clusters; (3) an optional step for aligning short-reads is provided, to validate split-sites by identifying unique junction reads, and revealing and quantifying tissue-specific alternative splice isoforms. We tested IsoSplitter AS prediction using datasets from multiple model and non-model plant species, and showed that IsoSplitter pipeline is efficient to handle different transcriptomes with high accuracy. Furthermore, we evaluated the IsoSplitter pipeline compared with that of the splice junction identification tools, Program to Assemble Spliced Alignments (PASA-software needs a reference genome for AS identification) and AStrap, using data from the model plant Arabidopsis thaliana. We found that, IsoSplitter determined more than twice as many AS events than AStrap analysis; and 94.13% of the IsoSplitter predicted AS events were also identified by the PASA analysis. Starting from a simple sequence file, IsoSplitter is an assembly-free tool for identification and characterization of AS. IsoSplitter is developed and implemented in Python 3.5 using the Linux platform and is freely available at https://github.com/Hengfu-Yin/IsoSplitter.

ULTRA: A Model Based Tool to Detect Tandem Repeats

In biological sequences, tandem repeats consist of tens to hundreds of residues of a repeated pattern, such as atgatgatgatgatg ('atg' repeated), often the result of replication slippage. Over time, these repeats decay so that the original sharp pattern of repetition is somewhat obscured, but even degenerate repeats pose a problem for sequence annotation: when two sequences both contain shared patterns of similar repetition, the result can be a false signal of sequence homology. We describe an implementation of a new hidden Markov model for detecting tandem repeats that shows substantially improved sensitivity to labeling decayed repetitive regions, presents low and reliable false annotation rates across a wide range of sequence composition, and produces scores that follow a stable distribution. On typical genomic sequence, the time and memory requirements of the resulting tool (ULTRA) are competitive with the most heavily used tool for repeat masking (TRF). ULTRA is released under an open source license and lays the groundwork for inclusion of the model in sequence alignment tools and annotation pipelines.

The hypertrophic cardiomyopathy myosin mutation R453C alters ATP binding and hydrolysis of human cardiac β-myosin

The human hypertrophic cardiomyopathy mutation R453C results in one of the more severe forms of the myopathy. Arg-453 is found in a conserved surface loop of the upper 50-kDa domain of the myosin motor domain and lies between the nucleotide binding pocket and the actin binding site. It connects to the cardiomyopathy loop via a long α-helix, helix O, and to Switch-2 via the fifth strand of the central β-sheet. The mutation is, therefore, in a position to perturb a wide range of myosin molecular activities. We report here the first detailed biochemical kinetic analysis of the motor domain of the human β-cardiac myosin carrying the R453C mutation. A recent report of the same mutation (Sommese, R. F., Sung, J., Nag, S., Sutton, S., Deacon, J. C., Choe, E., Leinwand, L. A., Ruppel, K., and Spudich, J. A. (2013) Proc. Natl. Acad. Sci. U.S.A. 110, 12607-12612) found reduced ATPase and in vitro motility but increased force production using an optical trap. Surprisingly, our results show that the mutation alters few biochemical kinetic parameters significantly. The exceptions are the rate constants for ATP binding to the motor domain (reduced by 35%) and the ATP hydrolysis step/recovery stroke (slowed 3-fold), which could be the rate-limiting step for the ATPase cycle. Effects of the mutation on the recovery stroke are consistent with a perturbation of Switch-2 closure, which is required for the recovery stroke and the subsequent ATP hydrolysis.

The superfast human extraocular myosin is kinetically distinct from the fast skeletal IIa, IIb, and IId isoforms

Humans express five distinct myosin isoforms in the sarcomeres of adult striated muscle (fast IIa, IId, the slow/cardiac isoform I/β, the cardiac specific isoform α, and the specialized extraocular muscle isoform). An additional isoform, IIb, is present in the genome but is not normally expressed in healthy human muscles. Muscle fibers expressing each isoform have distinct characteristics including shortening velocity. Defining the properties of the isoforms in detail has been limited by the availability of pure samples of the individual proteins. Here we study purified recombinant human myosin motor domains expressed in mouse C₂C₁₂ muscle cells. The results of kinetic analysis show that among the closely related adult skeletal isoforms, the affinity of ADP for actin·myosin (K_AD) is the characteristic that most readily distinguishes the isoforms. The three fast muscle myosins have K_AD values of 118, 80, and 55 μm for IId, IIa, and IIb, respectively, which follows the speed in motility assays from fastest to slowest. Extraocular muscle is unusually fast with a far weaker K_AD = 352 μm. Sequence comparisons and homology modeling of the structures identify a few key areas of sequence that may define the differences between the isoforms, including a region of the upper 50-kDa domain important in signaling between the nucleotide pocket and the actin-binding site.