Visualizing profile-profile alignment: pairwise HMM logos

Genome-wide characterization, evolution, structure, and expression analysis of the F-box genes in Caenorhabditis

BACKGROUND

F-box proteins represent a diverse class of adaptor proteins of the ubiquitin-proteasome system (UPS) that play critical roles in the cell cycle, signal transduction, and immune response by removing or modifying cellular regulators. Among closely related organisms of the Caenorhabditis genus, remarkable divergence in F-box gene copy numbers was caused by sizeable species-specific expansion and contraction. Although F-box gene number expansion plays a vital role in shaping genomic diversity, little is known about molecular evolutionary mechanisms responsible for substantial differences in gene number of F-box genes and their functional diversification in Caenorhabditis. Here, we performed a comprehensive evolution and underlying mechanism analysis of F-box genes in five species of Caenorhabditis genus, including C. brenneri, C. briggsae, C. elegans, C. japonica, and C. remanei.

RESULTS

Herein, we identified and characterized 594, 192, 377, 39, 1426 F-box homologs encoding putative F-box proteins in the genome of C. brenneri, C. briggsae, C. elegans, C. japonica, and C. remanei, respectively. Our work suggested that extensive species-specific tandem duplication followed by a small amount of gene loss was the primary mechanism responsible for F-box gene number divergence in Caenorhabditis genus. After F-box gene duplication events occurred, multiple mechanisms have contributed to gene structure divergence, including exon/intron gain/loss, exonization/pseudoexonization, exon/intron boundaries alteration, exon splits, and intron elongation by tandem repeats. Based on high-throughput RNA sequencing data analysis, we proposed that F-box gene functions have diversified by sub-functionalization through highly divergent stage-specific expression patterns in Caenorhabditis species.

CONCLUSIONS

Massive species-specific tandem duplications and occasional gene loss drove the rapid evolution of the F-box gene family in Caenorhabditis, leading to complex gene structural variation and diversified functions affecting growth and development within and among Caenorhabditis species. In summary, our findings outline the evolution of F-box genes in the Caenorhabditis genome and lay the foundation for future functional studies.

MCAW-DB: A glycan profile database capturing the ambiguity of glycan recognition patterns

Glycan-binding protein (GBP) interaction experiments, such as glycan microarrays, are often used to understand glycan recognition patterns. However, oftentimes the interpretation of glycan array experimental data makes it difficult to identify discrete GBP binding patterns due to their ambiguity. It is known that lectins, for example, are non-specific in their binding affinities; the same lectin can bind to different monosaccharides or even different glycan structures. In bioinformatics, several tools to mine the data generated from these sorts of experiments have been developed. These tools take a library of predefined motifs, which are commonly-found glycan patterns such as sialyl-Lewis X, and attempt to identify the motif(s) that are specific to the GBP being analyzed. In our previous work, as opposed to using predefined motifs, we developed the Multiple Carbohydrate Alignment with Weights (MCAW) tool to visualize the state of the glycans being recognized by the GBP under analysis. We previously reported on the effectiveness of our tool and algorithm by analyzing several glycan array datasets from the Consortium of Functional Glycomics (CFG). In this work, we report on our analysis of 1081 data sets which we collected from the CFG, the results of which we have made publicly and freely available as a database called MCAW-DB. We introduce this database, its usage and describe several analysis results. We show how MCAW-DB can be used to analyze glycan-binding patterns of GBPs amidst their ambiguity. For example, the visualization of glycan-binding patterns in MCAW-DB show how they correlate with the concentrations of the samples used in the array experiments. Using MCAW-DB, the patterns of glycans found to bind to various GBP-glycan binding proteins are visualized, indicating the binding "environment" of the glycans. Thus, the ambiguity of glycan recognition is numerically represented, along with the patterns of monosaccharides surrounding the binding region. The profiles in MCAW-DB could potentially be used as predictors of affinity of unknown or novel glycans to particular GBPs by comparing how well they match the existing profiles for those GBPs. Moreover, as the glycan profiles of diseased tissues become available, glycan alignments could also be used to identify glycan biomarkers unique to that tissue. Databases of these alignments may be of great use for drug discovery.

A Toxin-Antitoxin System VapBC15 from Synechocystis sp. PCC 6803 Shows Distinct Regulatory Features

Type II toxin-antitoxin (TA) systems play important roles in bacterial stress survival by regulating cell growth or death. They are highly abundant in cyanobacteria yet remain poorly characterized. Here, we report the identification and regulation of a putative type II TA system from Synechocystis PCC6803, VapBC15. The VapBC15 system is encoded by the chromosomal operon vapBC15. Exogenous expression of VapC15 dramatically arrested cell growth of Escherichia coli and reduced the numbers of colony-forming units (CFU). The VapC15 toxicity could be which was counteracted neutralized by simultaneous or delayed production of VapB15. Biochemical analysis demonstrated the formation of VapB15-VapC15 complexes by the physical interaction between VapB15 and VapC15. Notably, the VapB15 antitoxin up-regulated the transcription of the vapBC15 operon by directly binding to the promoter region, and the VapC15 toxin abolished the up-regulatory effect by destabilizing the binding. Moreover, VapB15 can be degraded by the proteases Lons and ClpXP2s from Synechocystis PCC6803, thus activating the latent toxicity of VapBC15. These findings suggest that VapBC15 represents a genuine TA system that utilizes a distinct mechanism to regulate toxin activity.

Expanding the repertoire of secretory peptides controlling root development with comparative genome analysis and functional assays

Plant genomes encode numerous small secretory peptides (SSPs) whose functions have yet to be explored. Based on structural features that characterize SSP families known to take part in postembryonic development, this comparative genome analysis resulted in the identification of genes coding for oligopeptides potentially involved in cell-to-cell communication. Because genome annotation based on short sequence homology is difficult, the criteria for the de novo identification and aggregation of conserved SSP sequences were first benchmarked across five reference plant species. The resulting gene families were then extended to 32 genome sequences, including major crops. The global phylogenetic pattern common to the functionally characterized SSP families suggests that their apparition and expansion coincide with that of the land plants. The SSP families can be searched online for members, sequences and consensus (http://bioinformatics.psb.ugent.be/webtools/PlantSSP/). Looking for putative regulators of root development, Arabidopsis thaliana SSP genes were further selected through transcriptome meta-analysis based on their expression at specific stages and in specific cell types in the course of the lateral root formation. As an additional indication that formerly uncharacterized SSPs may control development, this study showed that root growth and branching were altered by the application of synthetic peptides matching conserved SSP motifs, sometimes in very specific ways. The strategy used in the study, combining comparative genomics, transcriptome meta-analysis and peptide functional assays in planta, pinpoints factors potentially involved in non-cell-autonomous regulatory mechanisms. A similar approach can be implemented in different species for the study of a wide range of developmental programmes.

In silico identification of AMPylating enzymes and study of their divergent evolution

AMPylation is a novel post-translational modification (PTM) involving covalent attachment of an AMP moiety to threonine/tyrosine side chains of a protein. AMPylating enzymes belonging to three different families, namely Fic/Doc, GS-ATase and DrrA have been experimentally characterized. Involvement of these novel enzymes in a myriad of biological processes makes them interesting candidates for genome-wide search. We have used SVM and HMM to develop a computational protocol for identification of AMPylation domains and their classification into various functional subfamilies catalyzing AMPylation, deAMPylation, phosphorylation and phosphocholine transfer. Our analysis has not only identified novel PTM catalyzing enzymes among unannotated proteins, but has also revealed how this novel enzyme family has evolved to generate functional diversity by subtle changes in sequence/structures of the proteins. Phylogenetic analysis of Fic/Doc has revealed three new isofunctional subfamilies, thus adding to their functional divergence. Also, frequent occurrence of Fic/Doc proteins on highly mobile and unstable genomic islands indicated their evolution via extensive horizontal gene transfers. On the other hand phylogenetic analyses indicate lateral evolution of GS-ATase family and an early duplication event responsible for AMPylation and deAMPylation activity of GS-ATase. Our analysis also reveals molecular basis of substrate specificity of DrrA proteins.

Structure and mechanism of the phycobiliprotein lyase CpcT

Pigmentation of light-harvesting phycobiliproteins of cyanobacteria requires covalent attachment of open-chain tetrapyrroles, bilins, to the apoproteins. Thioether formation via addition of a cysteine residue to the 3-ethylidene substituent of bilins is mediated by lyases. T-type lyases are responsible for attachment to Cys-155 of phycobiliprotein β-subunits. We present crystal structures of CpcT (All5339) from Nostoc (Anabaena) sp. PCC 7120 and its complex with phycocyanobilin at 1.95 and 2.50 Å resolution, respectively. CpcT forms a dimer and adopts a calyx-shaped β-barrel fold. Although the overall structure of CpcT is largely retained upon chromophore binding, arginine residues at the opening of the binding pocket undergo major rotameric rearrangements anchoring the propionate groups of phycocyanobilin. Based on the structure and mutational analysis, a reaction mechanism is proposed that accounts for chromophore stabilization and regio- and stereospecificity of the addition reaction. At the dimer interface, a loop extending from one subunit partially shields the opening of the phycocyanobilin binding pocket in the other subunit. Deletion of the loop or disruptions of the dimer interface significantly reduce CpcT lyase activity, suggesting functional relevance of the dimer. Dimerization is further enhanced by chromophore binding. The chromophore is largely buried in the dimer, but in the monomer, the 3-ethylidene group is accessible for the apophycobiliprotein, preferentially from the chromophore α-side. Asp-163 and Tyr-65 at the β- and α-face near the E-configured ethylidene group, respectively, support the acid-catalyzed nucleophilic Michael addition of cysteine 155 of the apoprotein to an N-acylimmonium intermediate proposed by Grubmayr and Wagner (Grubmayr, K., and Wagner, U. G. (1988) Monatsh. Chem. 119, 965-983).

DUF581 is plant specific FCS-like zinc finger involved in protein-protein interaction

Zinc fingers are a ubiquitous class of protein domain with considerable variation in structure and function. Zf-FCS is a highly diverged group of C₂-C₂ zinc finger which is present in animals, prokaryotes and viruses, but not in plants. In this study we identified that a plant specific domain of unknown function, DUF581 is a zf-FCS type zinc finger. Based on HMM-HMM comparison and signature motif similarity we named this domain as FCS-Like Zinc finger (FLZ) domain. A genome wide survey identified that FLZ domain containing genes are bryophytic in origin and this gene family is expanded in spermatophytes. Expression analysis of selected FLZ gene family members of A. thaliana identified an overlapping expression pattern suggesting a possible redundancy in their function. Unlike the zf-FCS domain, the FLZ domain found to be highly conserved in sequence and structure. Using a combination of bioinformatic and protein-protein interaction tools, we identified that FLZ domain is involved in protein-protein interaction.

The diversity and functions of choline sulphatases in microorganisms

Choline sulphates have two putative roles in microorganisms: as a reservoir of C, N and S and as osmoprotectants. Although there is no established connection to date regarding the relative distribution of these two functions in microbial communities, this information is crucial in determining the role of choline sulphate in soils, particularly in cultivated soils where S is limiting. Therefore, in order to establish such a connection, the diversity of choline sulphatase (betC) genes was investigated in this study using numerous fully sequenced microbes available in GenBank. Our genomic analyses revealed unequivocally that the betICBA operon is restricted to Rhizobiaceae family members, which live under symbiotic conditions that prevent elemental depletion. Together with the uniform genetic organisation of the betICBA operon in Rhizobiaceae, BetC appears to be both utilised for osmoprotection or S replenishment. In contrast, betC in a wide variety of free-living microbes (including fungi, archaea and bacteria) was found in a cassette encoding only BetC and a choline sulphate transporter, a configuration that appears to be responsible for fulfilling elemental S requirements. Lastly, the relatively high number of BetC sequences available allowed the identification of a specific signature sequence that discriminates between these two functions and also globally defines some conserved motifs in microbial choline sulphatases. Due to the widespread presence of BetC in microbes and the wide repartition of the betC cassette system, the potential importance of choline sulphatase in global S recycling requires further clarification.

Evolutionary history of the TBP-domain superfamily

The TATA binding protein (TBP) is an essential transcription initiation factor in Archaea and Eucarya. Bacteria lack TBP, and instead use sigma factors for transcription initiation. TBP has a symmetric structure comprising two repeated TBP domains. Using sequence, structural and phylogenetic analyses, we examine the distribution and evolutionary history of the TBP domain, a member of the helix-grip fold family. Our analyses reveal a broader distribution than for TBP, with TBP-domains being present across all three domains of life. In contrast to TBP, all other characterized examples of the TBP domain are present as single copies, primarily within multidomain proteins. The presence of the TBP domain in the ubiquitous DNA glycosylases suggests that this fold traces back to the ancestor of all three domains of life. The TBP domain is also found in RNase HIII, and phylogenetic analyses show that RNase HIII has evolved from bacterial RNase HII via TBP-domain fusion. Finally, our comparative genomic screens confirm and extend earlier reports of proteins consisting of a single TBP domain among some Archaea. These monopartite TBP-domain proteins suggest that this domain is functional in its own right, and that the TBP domain could have first evolved as an independent protein, which was later recruited in different contexts.

Regulation of expression and catalytic activity of Escherichia coli RsmG methyltransferase

RsmG is an AdoMet-dependent methyltransferase responsible for the synthesis of m(7)G527 in the 530 loop of bacterial 16S rRNA. This loop is universally conserved, plays a key role in ribosomal accuracy, and is a target for streptomycin binding. Loss of the m(7)G527 modification confers low-level streptomycin resistance and may affect ribosomal functioning. Here, we explore the mechanisms controlling RsmG expression and activity, which may somehow respond to the demand set by the amount of rRNA. We confirm that rsmG is the second member in a bicistronic operon and demonstrate that rsmG also has its own promoter, which appears, in actively growing cells, as a control device to offset both the relatively low stability of RsmG and inhibition of the operon promoter. RsmG levels decrease under conditions that down-regulate rRNA synthesis. However, coordination between rRNA and RsmG expression does not seem to occur at the level of transcription initiation. Instead, it might depend on the activity of an inverted repeated region, located between the rsmG promoter and ribosome binding site, which we show to work as a weak transcriptional terminator. To gain insights into the enzymatic mechanism of RsmG, highly conserved residues were mutated and the abilities of the resulting proteins to confer streptomycin resistance, to modify rRNA, and to bind AdoMet were explored. Our data demonstrate for the first time the critical importance of some residues located in the active site of Escherichia coli RsmG for the m(7)G modification process and suggest a role for them in rRNA binding and catalysis.

Elucidating essential role of conserved carboxysomal protein CcmN reveals common feature of bacterial microcompartment assembly

Bacterial microcompartments are organelles composed of a protein shell that surrounds functionally related proteins. Bioinformatic analysis of sequenced genomes indicates that homologs to shell protein genes are widespread among bacteria and suggests that the shell proteins are capable of encapsulating diverse enzymes. The carboxysome is a bacterial microcompartment that enhances CO(2) fixation in cyanobacteria and some chemoautotrophs by sequestering ribulose-1,5-bisphosphate carboxylase/oxygenase and carbonic anhydrase in the microcompartment shell. Here, we report the in vitro and in vivo characterization of CcmN, a protein of previously unknown function that is absolutely conserved in β-carboxysomal gene clusters. We show that CcmN localizes to the carboxysome and is essential for carboxysome biogenesis. CcmN has two functionally distinct regions separated by a poorly conserved linker. The N-terminal portion of the protein is important for interaction with CcmM and, by extension, ribulose-1,5-bisphosphate carboxylase/oxygenase and the carbonic anhydrase CcaA, whereas the C-terminal peptide is essential for interaction with the carboxysome shell. Deletion of the peptide abolishes carboxysome formation, indicating that its interaction with the shell is an essential step in microcompartment formation. Peptides with similar length and sequence properties to those in CcmN can be bioinformatically detected in a large number of diverse proteins proposed to be encapsulated in functionally distinct microcompartments, suggesting that this peptide and its interaction with its cognate shell proteins are common features of microcompartment assembly.

Endophytic life strategies decoded by genome and transcriptome analyses of the mutualistic root symbiont Piriformospora indica

Recent sequencing projects have provided deep insight into fungal lifestyle-associated genomic adaptations. Here we report on the 25 Mb genome of the mutualistic root symbiont Piriformospora indica (Sebacinales, Basidiomycota) and provide a global characterization of fungal transcriptional responses associated with the colonization of living and dead barley roots. Extensive comparative analysis of the P. indica genome with other Basidiomycota and Ascomycota fungi that have diverse lifestyle strategies identified features typically associated with both, biotrophism and saprotrophism. The tightly controlled expression of the lifestyle-associated gene sets during the onset of the symbiosis, revealed by microarray analysis, argues for a biphasic root colonization strategy of P. indica. This is supported by a cytological study that shows an early biotrophic growth followed by a cell death-associated phase. About 10% of the fungal genes induced during the biotrophic colonization encoded putative small secreted proteins (SSP), including several lectin-like proteins and members of a P. indica-specific gene family (DELD) with a conserved novel seven-amino acids motif at the C-terminus. Similar to effectors found in other filamentous organisms, the occurrence of the DELDs correlated with the presence of transposable elements in gene-poor repeat-rich regions of the genome. This is the first in depth genomic study describing a mutualistic symbiont with a biphasic lifestyle. Our findings provide a significant advance in understanding development of biotrophic plant symbionts and suggest a series of incremental shifts along the continuum from saprotrophy towards biotrophy in the evolution of mycorrhizal association from decomposer fungi.

Phylogenetic distribution and evolutionary history of bacterial DEAD-Box proteins

DEAD-box proteins are found in all domains of life and participate in almost all cellular processes that involve RNA. The presence of DEAD and Helicase_C conserved domains distinguish these proteins. DEAD-box proteins exhibit RNA-dependent ATPase activity in vitro, and several also show RNA helicase activity. In this study, we analyzed the distribution and architecture of DEAD-box proteins among bacterial genomes to gain insight into the evolutionary pathways that have shaped their history. We identified 1,848 unique DEAD-box proteins from 563 bacterial genomes. Bacterial genomes can possess a single copy DEAD-box gene, or up to 12 copies of the gene, such as in Shewanella. The alignment of 1,208 sequences allowed us to perform a robust analysis of the hallmark motifs of DEAD-box proteins and determine the residues that occur at high frequency, some of which were previously overlooked. Bacterial DEAD-box proteins do not generally contain a conserved C-terminal domain, with the exception of some members that possess a DbpA RNA-binding domain (RBD). Phylogenetic analysis showed a separation of DbpA-RBD-containing and DbpA-RBD-lacking sequences and revealed a group of DEAD-box protein genes that expanded mainly in the Proteobacteria. Analysis of DEAD-box proteins from Firmicutes and γ-Proteobacteria, was used to deduce orthologous relationships of the well-studied DEAD-box proteins from Escherichia coli and Bacillus subtilis. These analyses suggest that DbpA-RBD is an ancestral domain that most likely emerged as a specialized domain of the RNA-dependent ATPases. Moreover, these data revealed numerous events of gene family expansion and reduction following speciation.

ZINC-INDUCED FACILITATOR-LIKE family in plants: lineage-specific expansion in monocotyledons and conserved genomic and expression features among rice (Oryza sativa) paralogs

BACKGROUND

Duplications are very common in the evolution of plant genomes, explaining the high number of members in plant gene families. New genes born after duplication can undergo pseudogenization, neofunctionalization or subfunctionalization. Rice is a model for functional genomics research, an important crop for human nutrition and a target for biofortification. Increased zinc and iron content in the rice grain could be achieved by manipulation of metal transporters. Here, we describe the ZINC-INDUCED FACILITATOR-LIKE (ZIFL) gene family in plants, and characterize the genomic structure and expression of rice paralogs, which are highly affected by segmental duplication.

RESULTS

Sequences of sixty-eight ZIFL genes, from nine plant species, were comparatively analyzed. Although related to MSF_1 proteins, ZIFL protein sequences consistently grouped separately. Specific ZIFL sequence signatures were identified. Monocots harbor a larger number of ZIFL genes in their genomes than dicots, probably a result of a lineage-specific expansion. The rice ZIFL paralogs were named OsZIFL1 to OsZIFL13 and characterized. The genomic organization of the rice ZIFL genes seems to be highly influenced by segmental and tandem duplications and concerted evolution, as rice genome contains five highly similar ZIFL gene pairs. Most rice ZIFL promoters are enriched for the core sequence of the Fe-deficiency-related box IDE1. Gene expression analyses of different plant organs, growth stages and treatments, both from our qPCR data and from microarray databases, revealed that the duplicated ZIFL gene pairs are mostly co-expressed. Transcripts of OsZIFL4, OsZIFL5, OsZIFL7, and OsZIFL12 accumulate in response to Zn-excess and Fe-deficiency in roots, two stresses with partially overlapping responses.

CONCLUSIONS

We suggest that ZIFL genes have different evolutionary histories in monocot and dicot lineages. In rice, concerted evolution affected ZIFL duplicated genes, possibly maintaining similar expression patterns between pairs. The enrichment for IDE1 boxes in rice ZIFL gene promoters suggests a role in Zn-excess and Fe-deficiency up-regulation of ZIFL transcripts. Moreover, this is the first description of the ZIFL gene family in plants and the basis for functional studies on this family, which may play important roles in Zn and Fe homeostasis in plants.

Analysis of subgroup C of fungal chitinases containing chitin-binding and LysM modules in the mycoparasite Trichoderma atroviride

Fungi have a plethora of chitinases, which can be phylogenetically divided into three subgroups (A, B and C). Subgroup C (sgC) chitinases are especially interesting due to their multiple carbohydrate-binding modules, but they have not been investigated in detail yet. In this study, we analyzed sgC chitinases in the mycoparasites Trichoderma atroviride and Trichoderma virens. The expression of sgC chitinase genes in T. atroviride was induced during mycoparasitism of the fungal prey Botrytis cinerea, but not Rhizoctonia solani and correspondingly only by fungal cell walls of the former. Interestingly, only few sgC chitinase genes were inducible by chitin, suggesting that non-chitinous cell wall components can act as inducers. In contrast, the transcriptional profile of the most abundantly expressed sgC chitinase gene tac6 indicated a role of the protein in hyphal network formation. This shows that sgC chitinases have diverse functions and are not only involved in the mycoparasitic attack. However, sequence analysis and 3D modelling revealed that TAC6 and also its ortholog in T. virens have potentially detrimental deletions in the substrate-binding site and are thus probably not catalytically active enzymes. Genomic analysis showed that the genes neighboring sgC chitinases often encode proteins that are solely composed of multiple LysM modules, which were induced by similar stimuli as their neighboring sgC chitinase genes. This study provides first insights into fungal sgC chitinases and their associated LysM proteins.

Visualization of multiple alignments, phylogenies and gene family evolution

Software for visualizing sequence alignments and trees are essential tools for life scientists. In this review, we describe the major features and capabilities of a selection of stand-alone and web-based applications useful when investigating the function and evolution of a gene family. These range from simple viewers, to systems that provide sophisticated editing and analysis functions. We conclude with a discussion of the challenges that these tools now face due to the flood of next generation sequence data and the increasingly complex network of bioinformatics information sources.

New insights into granin-derived peptides: evolution and endocrine roles

The granin protein family is composed of two chromogranin and five secretogranin members that are acidic, heat-stable proteins in secretory granules in cells of the nervous and endocrine systems. We report that there is little evidence for evolutionary relationships among the granins except for the chromogranin group. The main granin members, including chromogranin A and B, and secretogranin II are moderately conserved in the vertebrates. Several small bioactive peptides can be generated by proteolysis from those homologous domains existing within the granin precursors, reflecting the conservation of biological activities in different vertebrates. In this context, we focus on reviewing the distribution and function of the major granin-derived peptides, including vasostatin, bovine CgB(1-41) and secretoneurin in vertebrate endocrine systems, especially those associated with growth, glucose metabolism and reproduction.

HHsvm: fast and accurate classification of profile-profile matches identified by HHsearch

MOTIVATION

Recently developed profile-profile methods rival structural comparisons in their ability to detect homology between distantly related proteins. Despite this tremendous progress, many genuine relationships between protein families cannot be recognized as comparisons of their profiles result in scores that are statistically insignificant.

RESULTS

Using known evolutionary relationships among protein superfamilies in SCOP database, support vector machines were trained on four sets of discriminatory features derived from the output of HHsearch. Upon validation, it was shown that the automatic classification of all profile-profile matches was superior to fixed threshold-based annotation in terms of sensitivity and specificity. The effectiveness of this approach was demonstrated by annotating several domains of unknown function from the Pfam database.

AVAILABILITY

Programs and scripts implementing the methods described in this manuscript are freely available from http://hhsvm.dlakiclab.org/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

webPRC: the Profile Comparer for alignment-based searching of public domain databases

Profile–profile methods are well suited to detect remote evolutionary relationships between protein families. Profile Comparer (PRC) is an existing stand-alone program for scoring and aligning hidden Markov models (HMMs), which are based on multiple sequence alignments. Since PRC compares profile HMMs instead of sequences, it can be used to find distant homologues. For this purpose, PRC is used by, for example, the CATH and Pfam-domain databases. As PRC is a profile comparer, it only reports profile HMM alignments and does not produce multiple sequence alignments. We have developed webPRC server, which makes it straightforward to search for distant homologues or similar alignments in a number of domain databases. In addition, it provides the results both as multiple sequence alignments and aligned HMMs. Furthermore, the user can view the domain annotation, evaluate the PRC hits with the Jalview multiple alignment editor and generate logos from the aligned HMMs or the aligned multiple alignments. Thus, this server assists in detecting distant homologues with PRC as well as in evaluating and using the results. The webPRC interface is available at http://www.ibi.vu.nl/programs/prcwww/.

Sequence signatures in envelope protein may determine whether flaviviruses produce hemorrhagic or encephalitic syndromes

We analyzed the envelope proteins in pathogenic flaviviruses to determine whether there are sequence signatures associated with the tendency of viruses to produce hemorrhagic disease (H-viruses) or encephalitis (E-viruses). We found that, at the position corresponding to the glycosylated Asn-67 in dengue virus, asparagine (Asn) occurs in all seven viral species that cause hemorrhagic disease in humans. Furthermore, Asn was extremely rare at position 67 in six flaviviruses that cause encephalitis, being replaced by Asp in four of them. Of the 3,246 sequences from H- and E-viruses, we found that 2,916 sequences (90%) contained Asn in position 67 for H-viruses or Asp in position 67 for E-viruses. The change from Asn-67 that is prevalent in H-viruses to Asp-67 (common in E-viruses) contributes to a stronger electrostatically negative surface in the E-viruses as compared to the H-viruses. These findings should help predicting the disease potential of emerging and re-emerging flaviviruses and understanding the relationship between protein structure and disease outcome.

Profile Comparer: a program for scoring and aligning profile hidden Markov models

Summary: Profile Comparer (PRC) is a stand-alone program for scoring and aligning profile hidden Markov models (HMMs) of protein families. PRC can read models produced by SAM and HMMER, two popular profile HMM packages, as well as PSI-BLAST checkpoint files. This application note provides a brief description of the profile–profile algorithm used by PRC.

Availability: The C source code licensed under the GNU General Public Licence and Linux and Mac OS X binaries can be downloaded from http://supfam.org/PRC.

Contact:martin.madera@gmail.com

Supplementary information:Supplementary data are available at Bioinformatics online.

Evolution of genes and repeats in the Nimrod superfamily

The recently identified Nimrod superfamily is characterized by the presence of a special type of EGF repeat, the NIM repeat, located right after a typical CCXGY/W amino acid motif. On the basis of structural features, nimrod genes can be divided into three types. The proteins encoded by Draper-type genes have an EMI domain at the N-terminal part and only one copy of the NIM motif, followed by a variable number of EGF-like repeats. The products of Nimrod B-type and Nimrod C-type genes (including the eater gene) have different kinds of N-terminal domains, and lack EGF-like repeats but contain a variable number of NIM repeats. Draper and Nimrod C-type (but not Nimrod B-type) proteins carry a transmembrane domain. Several members of the superfamily were claimed to function as receptors in phagocytosis and/or binding of bacteria, which indicates an important role in the cellular immunity and the elimination of apoptotic cells. In this paper, the evolution of the Nimrod superfamily is studied with various methods on the level of genes and repeats. A hypothesis is presented in which the NIM repeat, along with the EMI domain, emerged by structural reorganizations at the end of an EGF-like repeat chain, suggesting a mechanism for the formation of novel types of repeats. The analyses revealed diverse evolutionary patterns in the sequences containing multiple NIM repeats. Although in the Nimrod B and Nimrod C proteins show characteristics of independent evolution, many internal NIM repeats in Eater sequences seem to have undergone concerted evolution. An analysis of the nimrod genes has been performed using phylogenetic and other methods and an evolutionary scenario of the origin and diversification of the Nimrod superfamily is proposed. Our study presents an intriguing example how the evolution of multigene families may contribute to the complexity of the innate immune response.

The solution structure of the C-terminal domain of NfeD reveals a novel membrane-anchored OB-fold

Nodulation formation efficiency D (NfeD) is a member of a class of membrane-anchored ClpP-class proteases. There is a second class of NfeD homologs that lack the ClpP domain. The genes of both NfeD classes usually are part of an operon that also contains a gene for a prokaryotic homolog of stomatin. (Stomatin is a major integral-membrane protein of mammalian erythrocytes.) Such NfeD/stomatin homolog gene pairs are present in more than 290 bacterial and archaeal genomes, and their protein products may be part of the machinery used for quality control of membrane proteins. Herein, we report the structure of the isolated C-terminal domain of PH0471, a Pyrococcus horikoshii NfeD homolog, which lacks the ClpP domain. This C-terminal domain (termed NfeDC) contains a five-strand beta-barrel, which is structurally very similar to the OB-fold (oligosaccharide/oligonucleotide-binding fold) domain. However, there is little sequence similarity between it and previously characterized OB-fold domains. The NfeDC domain lacks the conserved surface residues that are necessary for the binding of an OB-fold domain to DNA/RNA, an ion. Instead, its surface is composed of residues that are uniquely conserved in NfeD homologs and that form the structurally conserved surface turns and beta-bulges. There is also a conserved tryptophan present on the surface. We propose that, in general, NfeDC domains may interact with other spatially proximal membrane proteins and thereby regulate their activities.

The protist, Monosiga brevicollis, has a tyrosine kinase signaling network more elaborate and diverse than found in any known metazoan

Tyrosine kinase signaling has long been considered a hallmark of intercellular communication, unique to multicellular animals. Our genomic analysis of the unicellular choanoflagellate Monosiga brevicollis discovers a remarkable count of 128 tyrosine kinases, 38 tyrosine phosphatases, and 123 phosphotyrosine (pTyr)-binding SH2 proteins, all higher counts than seen in any metazoan. This elaborate signaling network shows little orthology to metazoan counterparts yet displays many innovations reminiscent of metazoans. These include extracellular domains structurally related to those of metazoan receptor kinases, alternative methods for membrane anchoring and phosphotyrosine interaction in cytoplasmic kinases, and domain combinations that link kinases to small GTPase signaling and transcription. These proteins also display a wealth of combinations of known signaling domains. This uniquely divergent and elaborate signaling network illuminates the early evolution of pTyr signaling, explores innovative ways to traverse the cellular signaling circuitry, and shows extensive convergent evolution, highlighting pervasive constraints on pTyr signaling.

An introduction to hidden Markov models

This unit introduces the concept of hidden Markov models in computational biology. It describes them using simple biological examples, requiring as little mathematical knowledge as possible. The unit also presents a brief history of hidden Markov models and an overview of their current applications before concluding with a discussion of their limitations.

SCOOP: a simple method for identification of novel protein superfamily relationships

MOTIVATION

Profile searches of sequence databases are a sensitive way to detect sequence relationships. Sophisticated profile-profile comparison algorithms that have been recently introduced increase search sensitivity even further.

RESULTS

In this article, a simpler approach than profile-profile comparison is presented that has a comparable performance to state-of-the-art tools such as COMPASS, HHsearch and PRC. This approach is called SCOOP (Simple Comparison Of Outputs Program), and is shown to find known relationships between families in the Pfam database as well as detect novel distant relationships between families. Several novel discoveries are presented including the discovery that a domain of unknown function (DUF283) found in Dicer proteins is related to double-stranded RNA-binding domains.

AVAILABILITY

SCOOP is freely available under a GNU GPL license from http://www.sanger.ac.uk/Users/agb/SCOOP/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Subfamily logos: visualization of sequence deviations at alignment positions with high information content

Background

Recognition of relevant sequence deviations can be valuable for elucidating functional differences between protein subfamilies. Interesting residues at highly conserved positions can then be mutated and experimentally analyzed. However, identification of such sites is tedious because automated approaches are scarce.

Results

Subfamily logos visualize subfamily-specific sequence deviations. The display is similar to classical sequence logos but extends into the negative range. Positive, upright characters correspond to residues which are characteristic for the subfamily, negative, upside-down characters to residues typical for the remaining sequences. The symbol height is adjusted to the information content of the alignment position. Residues which are conserved throughout do not appear.

Conclusion

Subfamily logos provide an intuitive display of relevant sequence deviations. The method has proven to be valid using a set of 135 aligned aquaporin sequences in which established subfamily-specific positions were readily identified by the algorithm.

Identification of multiple distinct Snf2 subfamilies with conserved structural motifs

The Snf2 family of helicase-related proteins includes the catalytic subunits of ATP-dependent chromatin remodelling complexes found in all eukaryotes. These act to regulate the structure and dynamic properties of chromatin and so influence a broad range of nuclear processes. We have exploited progress in genome sequencing to assemble a comprehensive catalogue of over 1300 Snf2 family members. Multiple sequence alignment of the helicase-related regions enables 24 distinct subfamilies to be identified, a considerable expansion over earlier surveys. Where information is known, there is a good correlation between biological or biochemical function and these assignments, suggesting Snf2 family motor domains are tuned for specific tasks. Scanning of complete genomes reveals all eukaryotes contain members of multiple subfamilies, whereas they are less common and not ubiquitous in eubacteria or archaea. The large sample of Snf2 proteins enables additional distinguishing conserved sequence blocks within the helicase-like motor to be identified. The establishment of a phylogeny for Snf2 proteins provides an opportunity to make informed assignments of function, and the identification of conserved motifs provides a framework for understanding the mechanisms by which these proteins function.

Pfam: clans, web tools and services

Pfam is a database of protein families that currently contains 7973 entries (release 18.0). A recent development in Pfam has enabled the grouping of related families into clans. Pfam clans are described in detail, together with the new associated web pages. Improvements to the range of Pfam web tools and the first set of Pfam web services that allow programmatic access to the database and associated tools are also presented. Pfam is available on the web in the UK (http://www.sanger.ac.uk/Software/Pfam/), the USA (http://pfam.wustl.edu/), France (http://pfam.jouy.inra.fr/) and Sweden (http://pfam.cgb.ki.se/).