We describe MUSCLE, a new computer program for creating multiple alignments of protein sequences. Elements of the algorithm include fast distance estimation using kmer counting, progressive alignment using a new profile function we call the logexpectation score, and refinement using treedependent restricted partitioning. The speed and accuracy of MUSCLE are compared with T-Coffee, MAFFT and CLUSTALW on four test sets of reference alignments: BAliBASE, SABmark, SMART and a new benchmark, PREFAB. MUSCLE achieves the highest, or joint highest, rank in accuracy on each of these sets. Without refinement, MUSCLE achieves average accuracy statistically indistinguishable from T-Coffee and MAFFT, and is the fastest of the tested methods for large numbers of sequences, aligning 5000 sequences of average length 350 in 7 min on a current desktop computer. The MUSCLE program, source code and PREFAB test data are freely available at http://www.drive5. com/muscle.

MUSCLE: multiple sequence alignment with high accuracy and high throughput

The human genome holds an extraordinary trove of information about human development, physiology, medicine and evolution. Here we report the results of an international collaboration to produce and make freely available a draft sequence of the human genome. We also present an initial analysis of the data, describing some of the insights that can be gleaned from the sequence.

/pdf/initial-sequencing-and-analysis-of-the-human-genome-5dapk1rsx1.pdf

Initial sequencing and analysis of the human genome.

抗原变异可使得多种致病微生物易于逃避宿主免疫应答。表达在感染红细胞表面的恶性疟原虫红细胞表面蛋白1（PfPMP1）与感染红细胞、内皮细胞、树突状细胞以及胎盘的单个或多个受体作用，在黏附及免疫逃避中起关键的作用。每个单倍体基因组var基因家族编码约60种成员，通过启动转录不同的var基因变异体为抗原变异提供了分子基础。

疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A

Pfam is a widely used database of protein families and domains. This article describes a set of major updates that we have implemented in the latest release (version 24.0). The most important change is that we now use HMMER3, the latest version of the popular profile hidden Markov model package. This software is approximately 100 times faster than HMMER2 and is more sensitive due to the routine use of the forward algorithm. The move to HMMER3 has necessitated numerous changes to Pfam that are described in detail. Pfam release 24.0 contains 11,912 families, of which a large number have been significantly updated during the past two years. Pfam is available via servers in the UK (http://pfam.sanger.ac.uk/), the USA (http://pfam.janelia.org/) and Sweden (http://pfam.sbc.su.se/).

/pdf/the-pfam-protein-families-database-3t3lguk949.pdf

The Pfam protein families database

Sequence similarity searching is a very important bioinformatics task. While Basic Local Alignment Search Tool (BLAST) outperforms exact methods through its use of heuristics, the speed of the current BLAST software is suboptimal for very long queries or database sequences. There are also some shortcomings in the user-interface of the current command-line applications. We describe features and improvements of rewritten BLAST software and introduce new command-line applications. Long query sequences are broken into chunks for processing, in some cases leading to dramatically shorter run times. For long database sequences, it is possible to retrieve only the relevant parts of the sequence, reducing CPU time and memory usage for searches of short queries against databases of contigs or chromosomes. The program can now retrieve masking information for database sequences from the BLAST databases. A new modular software library can now access subject sequence data from arbitrary data sources. We introduce several new features, including strategy files that allow a user to save and reuse their favorite set of options. The strategy files can be uploaded to and downloaded from the NCBI BLAST web site. The new BLAST command-line applications, compared to the current BLAST tools, demonstrate substantial speed improvements for long queries as well as chromosome length database sequences. We have also improved the user interface of the command-line applications.

/pdf/blast-architecture-and-applications-2udguz40ns.pdf

BLAST+: architecture and applications.

Insights from the architecture of the bacterial transcription apparatus

The ubiquitin (Ub) system is centered on conjugation and deconjugation of Ub and Ub-like (Ubls) proteins by a system of ligases and peptidases, respectively. Ub/Ubls contain the beta-grasp fold, also found in numerous proteins with biochemically distinct roles unrelated to the conventional Ub-system. The beta-GF underwent an early radiation spawning at least seven clades prior to the divergence of extant organisms from their last universal common ancestor, first emerging in the context of translation-related RNA-interactions and subsequently exploding to occupy various functional niches. Most beta-GF diversification occurred in prokaryotes, with the Ubl clade showing dramatic expansion in the eukaryotes. Diversification of Ubl families in eukaryotes played a major role in emergence of characteristic eukaryotic cellular sub-structures and systems. Recent comparative genomics studies indicate precursors of the eukaryotic Ub-system emerged in prokaryotes. The simplest of these combine an Ubl and an E1-like enzyme in metabolic pathways. Sampylation in archaea and Urmylation in eukaryotes appear to represent recruitment of such systems as simple protein-tagging apparatuses. However, other prokaryotic systems incorporated further components and mirror the eukaryotic condition in possessing an E2, a RING-type E3 or both of these components. Additionally, prokaryotes have evolved conjugation systems independent of Ub ligases, such as the Pup system.

The natural history of ubiquitin and ubiquitin-related domains.

Unlike classical 2-oxoglutarate and iron-dependent dioxygenases, which include several nucleic acid modifiers, the structurally similar jumonji-related dioxygenase superfamily was only known to catalyze peptide modifications. Using comparative genomics methods, we predict that a family of jumonji-related enzymes catalyzes wybutosine hydroxylation/peroxidation at position 37 of eukaryotic tRNAPhe. Identification of this enzyme raised questions regarding the emergence of protein- and nucleic acid-modifying activities among jumonji-related domains. We addressed these with a natural classification of DSBH domains and reconstructed the precursor of the dioxygenases as a sugar-binding domain. This precursor gave rise to sugar epimerases and metal-binding sugar isomerases. The sugar isomerase active site was exapted for catalysis of oxygenation, with a radiation of these enzymes in bacteria, probably due to impetus from the primary oxygenation event in Earth’s history. 2-Oxoglutarate-dependent versions appear to have further expanded with rise of the tricarboxylic acid cycle. We identify previously under-appreciated aspects of their active site and multiple independent innovations of 2-oxoacid-binding basic residues among these superfamilies. We show that double-stranded β-helix dioxygenases diversified extensively in biosynthesis and modification of halogenated siderophores, antibiotics, peptide secondary metabolites and glycine-rich collagen-like proteins in bacteria. Jumonji-related domains diversified into three distinct lineages in bacterial secondary metabolism systems and these were precursors of the three major clades of eukaryotic enzymes. The specificity of wybutosine hydroxylase/peroxidase probably relates to the structural similarity of the modified moiety to the ancestral amino acid substrate of this superfamily.

/pdf/origin-and-evolution-of-peptide-modifying-dioxygenases-and-ba5a1stu48.pdf

Origin and evolution of peptide-modifying dioxygenases and identification of the wybutosine hydroxylase/hydroperoxidase

Origin of multicellular eukaryotes - insights from proteome comparisons.

The availability of completely sequenced genomes and the wealth of literature on gene regulation have enabled researchers to model the transcriptional regulation system of some organisms in the form of a network. In order to reconstruct such networks in non-model organisms, three principal approaches have been taken. First, one can transfer the interactions between homologous components from a model organism to the organism of interest. Second, microarray experiments can be used to detect patterns in gene expression that stem from regulatory interactions. Finally, knowledge of experimentally characterized transcription factor binding sites can be used to analyze the promoter sequences in a genome in order to identify potential binding sites. In this chapter, we will focus in detail on the first approach and describe methods to reconstruct and analyze the transcriptional regulatory networks of uncharacterized organisms by using a known regulatory network as a template.

L. Aravind

Papers

Insights from the architecture of the bacterial transcription apparatus

The natural history of ubiquitin and ubiquitin-related domains.

Origin and evolution of peptide-modifying dioxygenases and identification of the wybutosine hydroxylase/hydroperoxidase

Origin of multicellular eukaryotes - insights from proteome comparisons.

Methods to reconstruct and compare transcriptional regulatory networks.