The BLAST programs are widely used tools for searching protein and DNA databases for sequence similarities. For protein comparisons, a variety of definitional, algorithmic and statistical refinements described here permits the execution time of the BLAST programs to be decreased substantially while enhancing their sensitivity to weak similarities. A new criterion for triggering the extension of word hits, combined with a new heuristic for generating gapped alignments, yields a gapped BLAST program that runs at approximately three times the speed of the original. In addition, a method is introduced for automatically combining statistically significant alignments produced by BLAST into a position-specific score matrix, and searching the database using this matrix. The resulting Position-Specific Iterated BLAST (PSIBLAST) program runs at approximately the same speed per iteration as gapped BLAST, but in many cases is much more sensitive to weak but biologically relevant sequence similarities. PSI-BLAST is used to uncover several new and interesting members of the BRCT superfamily.

/pdf/gapped-blast-and-psi-blast-a-new-generation-of-protein-2ybl67x5yx.pdf

Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.

Genomic sequencing has made it clear that a large fraction of the genes specifying the core biological functions are shared by all eukaryotes. Knowledge of the biological role of such shared proteins in one organism can often be transferred to other organisms. The goal of the Gene Ontology Consortium is to produce a dynamic, controlled vocabulary that can be applied to all eukaryotes even as knowledge of gene and protein roles in cells is accumulating and changing. To this end, three independent ontologies accessible on the World-Wide Web (http://www.geneontology.org) are being constructed: biological process, molecular function and cellular component.

/pdf/gene-ontology-tool-for-the-unification-of-biology-ppsktrcye3.pdf

Gene Ontology: tool for the unification of biology

We present the latest version of the Molecular Evolutionary Genetics Analysis (Mega) software, which contains many sophisticated methods and tools for phylogenomics and phylomedicine. In this major upgrade, Mega has been optimized for use on 64-bit computing systems for analyzing larger datasets. Researchers can now explore and analyze tens of thousands of sequences in Mega The new version also provides an advanced wizard for building timetrees and includes a new functionality to automatically predict gene duplication events in gene family trees. The 64-bit Mega is made available in two interfaces: graphical and command line. The graphical user interface (GUI) is a native Microsoft Windows application that can also be used on Mac OS X. The command line Mega is available as native applications for Windows, Linux, and Mac OS X. They are intended for use in high-throughput and scripted analysis. Both versions are available from www.megasoftware.net free of charge.

/pdf/mega7-molecular-evolutionary-genetics-analysis-version-7-0-1asv9rsrcf.pdf

MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets

MicroRNAs: Genomics, Biogenesis, Mechanism, and Function

We report a major update of the MAFFT multiple sequence alignment program. This version has several new features, including options for adding unaligned sequences into an existing alignment, adjustment of direction in nucleotide alignment, constrained alignment and parallel processing, which were implemented after the previous major update. This report shows actual examples to explain how these features work, alone and in combination. Some examples incorrectly aligned by MAFFT are also shown to clarify its limitations. We discuss how to avoid misalignments, and our ongoing efforts to overcome such limitations.

https://www.science-open.com/document_file/00d79e6c-98e4-425c-9fd1-3a85d5fdefbc/PubMedCentral/00d79e6c-98e4-425c-9fd1-3a85d5fdefbc.pdf

MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability

Fragile-X mental retardation is caused by loss of function of a single gene encoding the Fragile-X mental retardation protein, FMRP, an RNA-binding protein that harbors two KH-type and one RGG-type RNA-binding domains. Previous studies identified intramolecular G-quartet RNAs as high-affinity targets for the RGG box, but the relationship of RNA binding to FMRP function and mental retardation remains unclear. One severely affected patient harbors a missense mutation (I304N) within the second KH domain (KH2), and some evidence suggests this domain may be involved in the proposed role of FMRP in translational regulation. We now identify the RNA target for the KH2 domain as a sequence-specific element within a complex tertiary structure termed the FMRP kissing complex. We demonstrate that the association of FMRP with brain polyribosomes is abrogated by competition with the FMRP kissing complex RNA, but not by high-affinity G-quartet RNAs. We conclude that mental retardation associated with the I304N mutation, and likely the Fragile-X syndrome more generally, may relate to a crucial role for RNAs harboring the kissing complex motif as targets for FMRP translational regulation.

/pdf/kissing-complex-rnas-mediate-interaction-between-the-fragile-2njc9ohpay.pdf

Kissing complex RNAs mediate interaction between the Fragile-X mental retardation protein KH2 domain and brain polyribosomes

A Tree Viewer (ATV) is a Java tool for the display and manipulation of annotated phylogenetic trees. It can be utilized both as a standalone application and as an applet in a web browser.

/pdf/atv-display-and-manipulation-of-annotated-phylogenetic-trees-2qnf4amlrq.pdf

ATV: display and manipulation of annotated phylogenetic trees

Many sequence alignment programs use the BLOSUM62 score matrix to score pairs of aligned residues. Where did BLOSUM62 come from?

Where did the BLOSUM62 alignment score matrix come from

RNA secondary structure prediction methods based on probabilistic modeling can be developed using stochastic context-free grammars (SCFGs). Such methods can readily combine different sources of information that can be expressed probabilistically, such as an evolutionary model of comparative RNA sequence analysis and a biophysical model of structure plausibility. However, the number of free parameters in an integrated model for consensus RNA structure prediction can become untenable if the underlying SCFG design is too complex. Thus a key question is, what small, simple SCFG designs perform best for RNA secondary structure prediction? Nine different small SCFGs were implemented to explore the tradeoffs between model complexity and prediction accuracy. Each model was tested for single sequence structure prediction accuracy on a benchmark set of RNA secondary structures. Four SCFG designs had prediction accuracies near the performance of current energy minimization programs. One of these designs, introduced by Knudsen and Hein in their PFOLD algorithm, has only 21 free parameters and is significantly simpler than the others.

/pdf/evaluation-of-several-lightweight-stochastic-context-free-2o14fufn9m.pdf

Evaluation of several lightweight stochastic context-free grammars for RNA secondary structure prediction.

A new method, R-scape, tests whether observed sequence covariation supports a conserved secondary structure in RNA. The program finds no evidence for previously proposed conserved secondary structures in several lncRNAs. Many functional RNAs have an evolutionarily conserved secondary structure. Conservation of RNA base pairing induces pairwise covariations in sequence alignments. We developed a computational method, R-scape (RNA Structural Covariation Above Phylogenetic Expectation), that quantitatively tests whether covariation analysis supports the presence of a conserved RNA secondary structure. R-scape analysis finds no statistically significant support for proposed secondary structures of the long noncoding RNAs HOTAIR, SRA, and Xist.

/pdf/a-statistical-test-for-conserved-rna-structure-shows-lack-of-15q7dht4q5.pdf

Sean R. Eddy

Papers

Kissing complex RNAs mediate interaction between the Fragile-X mental retardation protein KH2 domain and brain polyribosomes

ATV: display and manipulation of annotated phylogenetic trees

Where did the BLOSUM62 alignment score matrix come from

Evaluation of several lightweight stochastic context-free grammars for RNA secondary structure prediction.

A statistical test for conserved RNA structure shows lack of evidence for structure in lncRNAs