TIGRFAMs: a protein family resource for the functional identification of proteins

A widespread picture of theStreptococcus thermophilus proteome by cell lysate fractionation and gel-based/gel-free approaches

Among the group of lactic acid bacteria, Streptococcus thermophilus has found a wide application in industrial processes used for the manufacture of dairy products. Taking advantage of different proteome extraction and subfractionation protocols, bacterial cytosolic and membrane proteins were isolated and resolved by independent gel-free and gel-based separation procedures. Whole cytosolic fraction and its acid, basic and low molecular mass protein components were separated by different resolutive 2-DE and tricine 1-DE gels and identified by MALDI-TOF PMF and/or microLC-ESI-IT-MS/MS. Membrane proteins were resolved by 2-DE and SDS-PAGE gels and similarly identified by PMF and TMS analysis. In parallel, whole extract was trypsinized and resulting peptides were identified by shotgun 2-D LC-ESI-IT-MS/MS analysis. Using this combined approach, expression products corresponding to 458 different genes were identified, which cover almost a third of the predicted vegetative proteome. Relative protein concentration and hydrophobicity affected protein detection. Broad recognition was obtained for enzymes involved in carbohydrate, fatty acid, amino acid and nucleotide metabolism, replication, transcription, translation, cell wall synthesis, as well as for proteins affecting bacterial functions important for industrial applications, i.e. milk sugar import and exopolysaccharide biosynthesis. By providing detailed reference electrophoretic/chromatographic maps to be used in future comparative proteomic investigations on bacteria grown under various experimental conditions or on different bacterial strains, our results will favour dedicated studies on S. thermophilus metabolism and its regulation or on detection of biomarkers for selection of optimal strains for industrial applications.

Complete genome sequence of the prototype lactic acid bacterium Lactococcus lactis subsp. cremoris MG1363

Lactococcus lactis is of great importance for the nutrition of hundreds of millions of people worldwide. This paper describes the genome sequence of Lactococcus lactis subsp. cremoris MG1363, the lactococcal strain most intensively studied throughout the world. The 2,529,478-bp genome contains 81 pseudogenes and encodes 2,436 proteins. Of the 530 unique proteins, 47 belong to the COG (clusters of orthologous groups) functional category "carbohydrate metabolism and transport," by far the largest category of novel proteins in comparison with L. lactis subsp. lactis IL1403. Nearly one-fifth of the 71 insertion elements are concentrated in a specific 56-kb region. This integration hot-spot region carries genes that are typically associated with lactococcal plasmids and a repeat sequence specifically found on plasmids and in the "lateral gene transfer hot spot" in the genome of Streptococcus thermophilus. Although the parent of L. lactis MG1363 was used to demonstrate lysogeny in Lactococcus, L. lactis MG1363 carries four remnant/satellite phages and two apparently complete prophages. The availability of the L. lactis MG1363 genome sequence will reinforce its status as the prototype among lactic acid bacteria through facilitation of further applied and fundamental research.

MBGD: a platform for microbial comparative genomics based on the automated construction of orthologous groups

The microbial genome database for comparative analysis (MBGD) is a comprehensive platform for microbial comparative genomics. The central function of MBGD is to create orthologous groups among multiple genomes from precomputed all-against-all similarity relationships using the DomClust algorithm. The database now contains >300 published genomes and the number continues to grow. For researchers who are interested in ongoing genome projects, we have now started a new service called ‘My MBGD,’ which allows users to add their own genome sequences to MBGD for the purpose of identifying orthologs among both the new and the existing genomes. Furthermore, in order to make available the rapidly accumulating information on closely related genome sequences, we enhanced the interface for pairwise genome comparisons using the CGAT interface, which allows users to see nucleotide sequence alignments of non-coding as well as coding regions. MBGD is available at .

Genome sequence of Avery's virulent serotype 2 strain D39 of Streptococcus pneumoniae and comparison with that of unencapsulated laboratory strain R6

Streptococcus pneumoniae (pneumococcus) is a leading human respiratory pathogen that causes a variety of serious mucosal and invasive diseases. D39 is an historically important serotype 2 strain that was used in experiments by Avery and coworkers to demonstrate that DNA is the genetic material. Although isolated nearly a century ago, D39 remains extremely virulent in murine infection models and is perhaps the strain used most frequently in current studies of pneumococcal pathogenesis. To date, the complete genome sequences have been reported for only two S. pneumoniae strains: TIGR4, a recent serotype 4 clinical isolate, and laboratory strain R6, an avirulent, unencapsulated derivative of strain D39. We report here the genome sequences and new annotation of two different isolates of strain D39 and the corrected sequence of strain R6. Comparisons of these three related sequences allowed deduction of the likely sequence of the D39 progenitor and mutations that arose in each isolate. Despite its numerous repeated sequences and IS elements, the serotype 2 genome has remained remarkably stable during cultivation, and one of the D39 isolates contains only five relatively minor mutations compared to the deduced D39 progenitor. In contrast, laboratory strain R6 contains 71 single-base-pair changes, six deletions, and four insertions and has lost the cryptic pDP1 plasmid compared to the D39 progenitor strain. Many of these mutations are in or affect the expression of genes that play important roles in regulation, metabolism, and virulence. The nature of the mutations that arose spontaneously in these three strains, the relative global transcription patterns determined by microarray analyses, and the implications of the D39 genome sequences to studies of pneumococcal physiology and pathogenesis are presented and discussed.

Evolution of sensory complexity recorded in a myxobacterial genome

Myxobacteria are single-celled, but social, eubacterial predators. Upon starvation they build multicellular fruiting bodies using a developmental program that progressively changes the pattern of cell movement and the repertoire of genes expressed. Development terminates with spore differentiation and is coordinated by both diffusible and cell-bound signals. The growth and development of Myxococcus xanthus is regulated by the integration of multiple signals from outside the cells with physiological signals from within. A collection of M. xanthus cells behaves, in many respects, like a multicellular organism. For these reasons M. xanthus offers unparalleled access to a regulatory network that controls development and that organizes cell movement on surfaces. The genome of M. xanthus is large (9.14 Mb), considerably larger than the other sequenced delta-proteobacteria. We suggest that gene duplication and divergence were major contributors to genomic expansion from its progenitor. More than 1,500 duplications specific to the myxobacterial lineage were identified, representing >15% of the total genes. Genes were not duplicated at random; rather, genes for cell-cell signaling, small molecule sensing, and integrative transcription control were amplified selectively. Families of genes encoding the production of secondary metabolites are overrepresented in the genome but may have been received by horizontal gene transfer and are likely to be important for predation.

Sequence analysis and organization of the Neodiprion abietis nucleopolyhedrovirus genome

Of 30 baculovirus genomes that have been sequenced to date, the only nonlepidopteran baculoviruses include the dipteran Culex nigripalpus nucleopolyhedrovirus and two hymenopteran nucleopolyhedroviruses that infect the sawflies Neodiprion lecontei (NeleNPV) and Neodiprion sertifer (NeseNPV). This study provides a complete sequence and genome analysis of the nucleopolyhedrovirus that infects the balsam fir sawfly Neodiprion abietis (Hymenoptera, Symphyta, Diprionidae). The N. abietis nucleopolyhedrovirus (NeabNPV) is 84,264 bp in size, with a G+C content of 33.5%, and contains 93 predicted open reading frames (ORFs). Eleven predicted ORFs are unique to this baculovirus, 10 ORFs have a putative sequence homologue in the NeleNPV genome but not the NeseNPV genome, and 1 ORF (neab53) has a putative sequence homologue in the NeseNPV genome but not the NeleNPV genome. Specific repeat sequences are coincident with major genome rearrangements that distinguish NeabNPV and NeleNPV. Genes associated with these repeat regions encode a common amino acid motif, suggesting that they are a family of repeated contiguous gene clusters. Lepidopteran baculoviruses, similarly, have a family of repeated genes called the bro gene family. However, there is no significant sequence similarity between the NeabNPV and bro genes. Homologues of early-expressed genes such as ie-1 and lef-3 were absent in NeabNPV, as they are in the previously sequenced hymenopteran baculoviruses. Analyses of ORF upstream sequences identified potential temporally distinct genes on the basis of putative promoter elements.

Exopolysaccharide-associated protein sorting in environmental organisms: the PEP-CTERM/EpsH system. Application of a novel phylogenetic profiling heuristic

Background

Protein translocation to the proper cellular destination may be guided by various classes of sorting signals recognizable in the primary sequence. Detection in some genomes, but not others, may reveal sorting system components by comparison of the phylogenetic profile of the class of sorting signal to that of various protein families.

Results

We describe a short C-terminal homology domain, sporadically distributed in bacteria, with several key characteristics of protein sorting signals. The domain includes a near-invariant motif Pro-Glu-Pro (PEP). This possible recognition or processing site is followed by a predicted transmembrane helix and a cluster rich in basic amino acids. We designate this domain PEP-CTERM. It tends to occur multiple times in a genome if it occurs at all, with a median count of eight instances; Verrucomicrobium spinosum has sixty-five. PEP-CTERM-containing proteins generally contain an N-terminal signal peptide and exhibit high diversity and little homology to known proteins. All bacteria with PEP-CTERM have both an outer membrane and exopolysaccharide (EPS) production genes. By a simple heuristic for screening phylogenetic profiles in the absence of pre-formed protein families, we discovered that a homolog of the membrane protein EpsH (exopolysaccharide locus protein H) occurs in a species when PEP-CTERM domains are found. The EpsH family contains invariant residues consistent with a transpeptidase function. Most PEP-CTERM proteins are encoded by single-gene operons preceded by large intergenic regions. In the Proteobacteria, most of these upstream regions share a DNA sequence, a probable cis-regulatory site that contains a sigma-54 binding motif. The phylogenetic profile for this DNA sequence exactly matches that of three proteins: a sigma-54-interacting response regulator (PrsR), a transmembrane histidine kinase (PrsK), and a TPR protein (PrsT).

Conclusion

These findings are consistent with the hypothesis that PEP-CTERM and EpsH form a protein export sorting system, analogous to the LPXTG/sortase system of Gram-positive bacteria, and correlated to EPS expression. It occurs preferentially in bacteria from sediments, soils, and biofilms. The novel method that led to these findings, partial phylogenetic profiling, requires neither global sequence clustering nor arbitrary similarity cutoffs and appears to be a rapid, effective alternative to other profiling methods.

Role of topoisomerase IIbeta in the expression of developmentally regulated genes

Mice lacking topoisomerase IIbeta (TopIIbeta) are known to exhibit a perinatal death phenotype. In the current study, transcription profiles of the brains of wild-type and top2beta knockout mouse embryos were generated. Surprisingly, only a small number (1 to 4%) of genes were affected in top2beta knockout embryos. However, the expression of nearly 30% of developmentally regulated genes was either up- or down-regulated. By contrast, the expression of genes encoding general cell growth functions and early differentiation markers was not affected, suggesting that TopIIbeta is not required for early differentiation programming but is specifically required for the expression of developmentally regulated genes at later stages of differentiation. Consistent with this notion, immunohistochemical analysis of brain sections showed that TopIIbeta and histone deacetylase 2, a known TopIIbeta-interacting protein, were preferentially expressed in neurons which are in their later stages of differentiation. Chromatin immunoprecipitation analysis of the developing brains revealed TopIIbeta binding to the 5' region of a number of TopIIbeta-sensitive genes. Further studies of a TopIIbeta-sensitive gene, Kcnd2, revealed the presence of TopIIbeta in the transcription unit with major binding near the promoter region. Together, these results support a role of TopIIbeta in activation/repression of developmentally regulated genes at late stages of neuronal differentiation.

Three-color cDNA microarrays with prehybridization quality control yield gene expression data comparable to that of commercial platforms

Despite their lower cost and high content flexibility, a limitation of in-house-prepared arrays has been their susceptibility to quality control (QC) issues and lack of QC standards across laboratories. Therefore, we developed a novel three-color array system that allows prehybridization QC as well as the Matarray software to facilitate acquisition of accurate gene expression data. In this study, we compared performance of our rat cDNA array to the Affymetrix RG-U34A and Agilent G4130A arrays using 2,824 UniGenes represented on all three arrays. Before data filtering, poor interplatform agreement was observed; however, after data filtering, differentially expressed UniGenes exhibited correlation coefficients of 0.91, 0.88, and 0.92 between the Affymetrix vs. Agilent, Affymetrix vs. cDNA, and Agilent vs. cDNA arrays, respectively. The Affymetrix, Agilent, and cDNA arrays agreed well with quantitative RT-PCR conducted on 42 UniGenes, yielding correlation coefficients of 0.90, 0.90, and 0.96, respectively. Each platform underestimated ratios relative to quantitative RT-PCR, possessing respective slopes of 0.86 ( R2 = 0.81), 0.65 ( R2 = 0.81), and 0.70 ( R2 = 0.92). Overall, these data show that the combination of our novel technical and analytic approaches yield an accurate platform for functional genomics that is concordant with commercial discovery arrays in terms of identifying regulated genes and pathways.

Bacterial genome adaptation to niches: divergence of the potential virulence genes in three Burkholderia species of different survival strategies

BACKGROUND

Two closely related species Burkholderia mallei (Bm) and Burkholderia pseudomallei (Bp) are serious human health hazards and are potential bio-warfare agents, whereas another closely related species Burkholderia thailandensis (Bt) is a non-pathogenic saprophyte. To investigate the genomic factors resulting in such a dramatic difference, we first identified the Bm genes responsive to the mouse environment, and then examined the divergence of these genes in Bp and Bt.

RESULTS

The genes down-expressed, which largely encode cell growth-related proteins, are conserved well in all three species, whereas those up-expressed, which include potential virulence genes, are less well conserved or absent notably in Bt. However, a substantial number of up-expressed genes is still conserved in Bt. Bm and Bp further diverged from each other in a small number of genes resulting from unit number changes in simple sequence repeats (ssr) in the homologs.

CONCLUSION

Our data suggest that divergent evolution of a small set of genes, rather than acquisition or loss of pathogenic islands, is associated with the development of different life styles in these bacteria of similar genomic contents. Further divergence between Bm and Bp mediated by ssr changes may reflect different adaptive processes of Bm and Bp fine-tuning into their host environments.

Defining genes in the genome of the hyperthermophilic archaeon Pyrococcus furiosus: implications for all microbial genomes

The original genome annotation of the hyperthermophilic archaeon Pyrococcus furiosus contained 2,065 open reading frames (ORFs). The genome was subsequently automatically annotated in two public databases by the Institute for Genomic Research (TIGR) and the National Center for Biotechnology Information (NCBI). Remarkably, more than 500 of the originally annotated ORFs differ in size in the two databases, many very significantly. For example, more than 170 of the predicted proteins differ at their N termini by more than 25 amino acids. Similar discrepancies were observed in the TIGR and NCBI databases with the other archaeal and bacterial genomes examined. In addition, the two databases contain 60 (NCBI) and 221 (TIGR) ORFs not present in the original annotation of P. furiosus. In the present study we have experimentally assessed the validity of 88 previously unannotated ORFs. Transcriptional analyses showed that 11 of 61 ORFs examined were expressed in P. furiosus when grown at either 95 or 72 degrees C. In addition, 7 of 54 ORFs examined yielded heat-stable recombinant proteins when they were expressed in Escherichia coli, although only one of the seven ORFs was expressed in P. furiosus under the growth conditions tested. It is concluded that the P. furiosus genome contains at least 17 ORFs not previously recognized in the original annotation. This study serves to highlight the discrepancies in the public databases and the problems of accurately defining the number and sizes of ORFs within any microbial genome.

Insights into genome plasticity and pathogenicity of the plant pathogenic bacterium Xanthomonas campestris pv. vesicatoria revealed by the complete genome sequence

The gram-negative plant-pathogenic bacterium Xanthomonas campestris pv. vesicatoria is the causative agent of bacterial spot disease in pepper and tomato plants, which leads to economically important yield losses. This pathosystem has become a well-established model for studying bacterial infection strategies. Here, we present the whole-genome sequence of the pepper-pathogenic Xanthomonas campestris pv. vesicatoria strain 85-10, which comprises a 5.17-Mb circular chromosome and four plasmids. The genome has a high G+C content (64.75%) and signatures of extensive genome plasticity. Whole-genome comparisons revealed a gene order similar to both Xanthomonas axonopodis pv. citri and Xanthomonas campestris pv. campestris and a structure completely different from Xanthomonas oryzae pv. oryzae. A total of 548 coding sequences (12.2%) are unique to X. campestris pv. vesicatoria. In addition to a type III secretion system, which is essential for pathogenicity, the genome of strain 85-10 encodes all other types of protein secretion systems described so far in gram-negative bacteria. Remarkably, one of the putative type IV secretion systems encoded on the largest plasmid is similar to the Icm/Dot systems of the human pathogens Legionella pneumophila and Coxiella burnetii. Comparisons with other completely sequenced plant pathogens predicted six novel type III effector proteins and several other virulence factors, including adhesins, cell wall-degrading enzymes, and extracellular polysaccharides.

The R1 resistance gene cluster contains three groups of independently evolving, type I R1 homologues and shows substantial structural variation among haplotypes of Solanum demissum

Cultivated and wild potatoes contain a major disease-resistance cluster on the short arm of chromosome V, including the R1 resistance (R) gene against potato late blight. To explore the functional and evolutionary significance of clustering in the generation of novel disease-resistance genes, we constructed three approximately 1 Mb physical maps in the R1 gene region, one for each of the three genomes (haplotypes) of allohexaploid Solanum demissum, the wild potato progenitor of the R1 locus. Totals of 691, 919 and 559 kb were sequenced for each haplotype, and three distinct resistance-gene families were identified, one homologous to the potato R1 gene and two others homologous to either the Prf or the Bs4 R-gene of tomato. The regions with R1 homologues are highly divergent among the three haplotypes, in contrast to the conserved flanking non-resistance gene regions. The R1 locus shows dramatic variation in overall length and R1 homologue number among the three haplotypes. Sequence comparisons of the R1 homologues show that they form three distinct clades in a distance tree. Frequent sequence exchanges were detected among R1 homologues within each clade, but not among those in different clades. These frequent sequence exchanges homogenized the intron sequences of homologues within each clade, but did not homogenize the coding sequences. Our results suggest that the R1 homologues represent three independent groups of fast-evolving type I resistance genes, characterized by chimeric structures resulting from frequent sequence exchanges among group members. Such genes were first identified among clustered RGC2 genes in lettuce, where they were distinguished from slow-evolving type II R-genes. Our findings at the R1 locus in S. demissum may indicate that a common or similar mechanism underlies the previously reported differentiation of type I and type II R-genes and the differentiation of type I R-genes into distinct groups, identified here.

A structure-based method for identifying DNA-binding proteins and their sites of DNA-interaction

A classification model of a DNA-binding protein chain was created based on identification of alpha helices within the chain likely to bind to DNA. Using the model, all chains in the Protein Data Bank were classified. For many of the chains classified with high confidence, previous documentation for DNA-binding was found, yet no sequence homology to the structures used to train the model was detected. The result indicates that the chain model can be used to supplement sequence based methods for annotating the function of DNA-binding. Four new candidates for DNA-binding were found, including two structures solved through structural genomics efforts. For each of the candidate structures, possible sites of DNA-binding are indicated by listing the residue ranges of alpha helices likely to interact with DNA.

Defective expression of Tamm-Horsfall protein/uromodulin in COX-2-deficient mice increases their susceptibility to urinary tract infections

Mice lacking a functional cyclooxygenase-2 (COX-2) gene develop abnormal kidneys that contain hypoplastic glomeruli and reduced proximal tubular mass, and they often die of renal failure. A comparison of kidney-specific gene expression between wild-type and COX-2-deficient mice by cDNA microarrays revealed that although more than 500 mRNAs were differentially expressed between the two strains of mice depending on their ages, the genes encoding pre-pro-epidermal growth factor (pre-pro-EGF) and Tamm-Horsfall protein (THP)/uromodulin were aberrantly expressed in the kidneys of COX-2 −/− mice at all stages of their development. Downregulation of EGF could potentially affect renal development, and THP/uromodulin gene has been implicated in abnormal kidney development and end-stage renal failure in humans. We assessed in detail mechanism of defective THP/uromodulin gene expression and its potential consequences in COX-2-deficient mice. Consistent with the microarray data, the steady-state levels of THP/uromodulin mRNA were severely reduced in the COX-2 −/− kidney. Furthermore, reduced expression of renal THP/uromodulin, as assessed by Western blot and immunohistological methods, was closely corroborated by a corresponding decline in the urinary secretion of THP/uromodulin in COX-2 −/− mice. Finally, we demonstrate that the bladders of COX-2 −/− mice, in contrast to those of the wild-type mice, are highly susceptible to colonization by uropathogenic Escherichia coli.

Detection and genotyping of Mycobacterium species from clinical isolates and specimens by oligonucleotide array

Identification of pathogenic Mycobacterium species is important for a successful diagnosis of mycobacteriosis. The purpose of this study was to develop an oligonucleotide array which could detect and differentiate mycobacteria to the species level by using the internal transcribed spacer (ITS) sequence. Using a genus-specific probe and 20 species-specific probes including two M. avium-intracellulare complex (MAC)-specific probes, we have developed an ITS-based oligonucleotide array for the rapid and reliable detection and discrimination of M. tuberculosis, MAC, M. fortuitum, M. chelonae, M. abscessus, M. kansasii, M. gordonae, M. scrofulaceum, M. szulgai, M. vaccae, M. xenopi, M. terrae, M. flavescens, M. smegmatis, M. malmoense, M. simiae, M. marinum, M. ulcerans, M. gastri, and M. leprae. All mycobacteria were hybridized with a genus-specific probe (PAN-03) for detection of the genus Mycobacterium. Mycobacterial species were expected to show a unique hybridization pattern with species-specific probes, except for M. marinum and M. ulcerans, which were not differentiated by ITS-based probe. Among the species-specific probes, two kinds of species-specific probes were designed for MAC in which there were many subspecies. The performance of the oligonucleotide array assay was demonstrated by using 46 reference strains, 149 clinical isolates, and 155 clinical specimens. The complete procedure (DNA extraction, PCR, DNA hybridization, and scanning) was carried out in 4.5 h. Our results indicated that the oligonucleotide array is useful for the identification and discrimination of mycobacteria from clinical isolates and specimens in an ordinary clinical laboratory.

PROTEIN FAMILIES AND THEIR EVOLUTION—A STRUCTURAL PERSPECTIVE

We can now assign about two thirds of the sequences from completed genomes to as few as 1400 domain families for which structures are known and thus more ancient evolutionary relationships established. About 200 of these domain families are common to all kingdoms of life and account for nearly 50% of domain structure annotations in the genomes. Some of these domain families have been very extensively duplicated within a genome and combined with different domain partners giving rise to different multidomain proteins. The ways in which these domain combinations evolve tend to be specific to the organism so that less than 15% of the protein families found within a genome appear to be common to all kingdoms of life. Recent analyses of completed genomes, exploiting the structural data, have revealed the extent to which duplication of these domains and modifications of their functions can expand the functional repertoire of the organism, contributing to increasing complexity.

The TetR family of transcriptional repressors

We have developed a general profile for the proteins of the TetR family of repressors. The stretch that best defines the profile of this family is made up of 47 amino acid residues that correspond to the helix-turn-helix DNA binding motif and adjacent regions in the three-dimensional structures of TetR, QacR, CprB, and EthR, four family members for which the function and three-dimensional structure are known. We have detected a set of 2,353 nonredundant proteins belonging to this family by screening genome and protein databases with the TetR profile. Proteins of the TetR family have been found in 115 genera of gram-positive, alpha-, beta-, and gamma-proteobacteria, cyanobacteria, and archaea. The set of genes they regulate is known for 85 out of the 2,353 members of the family. These proteins are involved in the transcriptional control of multidrug efflux pumps, pathways for the biosynthesis of antibiotics, response to osmotic stress and toxic chemicals, control of catabolic pathways, differentiation processes, and pathogenicity. The regulatory network in which the family member is involved can be simple, as in TetR (i.e., TetR bound to the target operator represses tetA transcription and is released in the presence of tetracycline), or more complex, involving a series of regulatory cascades in which either the expression of the TetR family member is modulated by another regulator or the TetR family member triggers a cell response to react to environmental insults. Based on what has been learned from the cocrystals of TetR and QacR with their target operators and from their three-dimensional structures in the absence and in the presence of ligands, and based on multialignment analyses of the conserved stretch of 47 amino acids in the 2,353 TetR family members, two groups of residues have been identified. One group includes highly conserved positions involved in the proper orientation of the helix-turn-helix motif and hence seems to play a structural role. The other set of less conserved residues are involved in establishing contacts with the phosphate backbone and target bases in the operator. Information related to the TetR family of regulators has been updated in a database that can be accessed at www.bactregulators.org.

Patellamide A and C biosynthesis by a microcin-like pathway in Prochloron didemni, the cyanobacterial symbiont of Lissoclinum patella

Prochloron spp. are obligate cyanobacterial symbionts of many didemnid family ascidians. It has been proposed that the cyclic peptides of the patellamide class found in didemnid extracts are synthesized by Prochloron spp., but studies in which host and symbiont cells are separated and chemically analyzed to identify the biosynthetic source have yielded inconclusive results. As part of the Prochloron didemni sequencing project, we identified patellamide biosynthetic genes and confirmed their function by heterologous expression of the whole pathway in Escherichia coli. The primary sequence of patellamides A and C is encoded on a single ORF that resembles a precursor peptide. We propose that this prepatellamide is heterocyclized to form thiazole and oxazoline rings, and the peptide is cleaved to yield the two cyclic patellamides, A and C. This work represents the full sequencing and functional expression of a marine natural-product pathway from an obligate symbiont. In addition, a related cluster was identified in Trichodesmium erythraeum IMS101, an important bloom-forming cyanobacterium.

Meta-DP: domain prediction meta-server

Meta-DP, a domain prediction meta-server provides a simple interface to predict domains in a given protein sequence using a number of domain prediction methods. The Meta-DP is a convenient resource because through accessing a single site, users automatically obtain the results of the various domain prediction methods along with a consensus prediction. The Meta-DP is currently coupled to 10 domain prediction servers and can be extended to include any number of methods. Meta-DP can thus become a centralized repository of available methods. Meta-DP was also used to evaluate the performance of 13 domain prediction methods in the context of CAFASP-DP.

AVAILABILITY

The Meta-DP server is freely available at http://meta-dp.bioinformatics.buffalo.edu and the CAFASP-DP evaluation results are available at http://cafasp4.bioinformatics.buffalo.edu/dp/update.html

CONTACT

hkaur@bioinformatics.buffalo.edu

SUPPLEMENTARY INFORMATION

Available at http://cafasp4.bioinformatics.buffalo.edu/dp/update.html.

Genomic characterization of non-O1, non-O139 Vibrio cholerae reveals genes for a type III secretion system

Non-O1, non-O139 Vibrio cholerae can cause gastroenteritis and extraintestinal infections, but, unlike O1 and O139 strains of V. cholerae, little is known about the virulence gene content of non-O1, non-O139 strains and their phylogenetic relationship to other pathogenic V. cholerae. Comparative genomic microarray analysis of four pathogenic non-O1, non-O139 strains indicates that these strains are quite divergent from O1 and O139 strains. Genomic sequence analysis of a non-O1, non-O139 strain (AM-19226) that appeared particularly pathogenic in experimental animals suggests that this strain carries a type III secretion system (TTSS) that is related to the TTSS2 gene cluster found in a pandemic clone of Vibrio parahaemolyticus. The genes for this V. cholerae TTSS system appear to be present in many clinical and environmental non-O1, non-O139 strains, including at least one clone that is globally distributed. We hypothesize that the TTSS present in some pathogenic strains of non-O1, non-O139 V. cholerae may be involved in the virulence and environmental fitness of these strains.

The Universal Protein Resource (UniProt)

The Universal Protein Resource (UniProt) provides the scientific community with a single, centralized, authoritative resource for protein sequences and functional information. Formed by uniting the Swiss-Prot, TrEMBL and PIR protein database activities, the UniProt consortium produces three layers of protein sequence databases: the UniProt Archive (UniParc), the UniProt Knowledgebase (UniProt) and the UniProt Reference (UniRef) databases. The UniProt Knowledgebase is a comprehensive, fully classified, richly and accurately annotated protein sequence knowledgebase with extensive cross-references. This centrepiece consists of two sections: UniProt/Swiss-Prot, with fully, manually curated entries; and UniProt/TrEMBL, enriched with automated classification and annotation. During 2004, tens of thousands of Knowledgebase records got manually annotated or updated; we introduced a new comment line topic: TOXIC DOSE to store information on the acute toxicity of a toxin; the UniProt keyword list got augmented by additional keywords; we improved the documentation of the keywords and are continuously overhauling and standardizing the annotation of post-translational modifications. Furthermore, we introduced a new documentation file of the strains and their synonyms. Many new database cross-references were introduced and we started to make use of Digital Object Identifiers. We also achieved in collaboration with the Macromolecular Structure Database group at EBI an improved integration with structural databases by residue level mapping of sequences from the Protein Data Bank entries onto corresponding UniProt entries. For convenient sequence searches we provide the UniRef non-redundant sequence databases. The comprehensive UniParc database stores the complete body of publicly available protein sequence data. The UniProt databases can be accessed online (http://www.uniprot.org) or downloaded in several formats (ftp://ftp.uniprot.org/pub). New releases are published every two weeks.

transition priors for protein hidden Markov models: an empirical study towards maximum discrimination

Insertions and deletions in a profile hidden Markov model (HMM) are modeled by transition probabilities between insert, delete and match states. These are estimated by combining observed data and prior probabilities. The transition prior probabilities can be defined either ad hoc or by maximum likelihood (ML) estimation. We show that the choice of transition prior greatly affects the HMM's ability to discriminate between true and false hits. HMM discrimination was measured using the HMMER 2.2 package applied to 373 families from Pfam. We measured the discrimination between true members and noise sequences employing various ML transition priors and also systematically scanned the parameter space of ad hoc transition priors. Our results indicate that ML priors produce far from optimal discrimination, and we present an empirically derived prior that considerably decreases the number of misclassifications compared to ML. Most of the difference stems from the probabilities for exiting a delete state. The ML prior, which is unaware of noise sequences, estimates a delete-to-delete probability that is relatively high and does not penalize noise sequences enough for optimal discrimination.

Cross-talk between activated human NK cells and CD4+ T cells via OX40-OX40 ligand interactions

It is important to understand which molecules are relevant for linking innate and adaptive immune cells. In this study, we show that OX40 ligand is selectively induced on IL-2, IL-12, or IL-15-activated human NK cells following stimulation through NKG2D, the low affinity receptor for IgG (CD16) or killer cell Ig-like receptor 2DS2. CD16-activated NK cells costimulate TCR-induced proliferation, and IFN-gamma produced by autologous CD4+ T cells and this process is dependent upon expression of OX40 ligand and B7 by the activated NK cells. These findings suggest a novel and unexpected link between the natural and specific immune responses, providing direct evidence for cross-talk between human CD4+ T cells and NK receptor-activated NK cells.

Biphenyl and benzoate metabolism in a genomic context: outlining genome-wide metabolic networks in Burkholderia xenovorans LB400

We designed and successfully implemented the use of in situ-synthesized 45-mer oligonucleotide DNA microarrays (XeoChips) for genome-wide expression profiling of Burkholderia xenovorans LB400, which is among the best aerobic polychlorinated biphenyl degraders known so far. We conducted differential gene expression profiling during exponential growth on succinate, benzoate, and biphenyl as sole carbon sources and investigated the transcriptome of early-stationary-phase cells grown on biphenyl. Based on these experiments, we outlined metabolic pathways and summarized other cellular functions in the organism relevant for biphenyl and benzoate degradation. All genes previously identified as being directly involved in biphenyl degradation were up-regulated when cells were grown on biphenyl compared to expression in succinate-grown cells. For benzoate degradation, however, genes for an aerobic coenzyme A activation pathway were up-regulated in biphenyl-grown cells, while the pathway for benzoate degradation via hydroxylation was up-regulated in benzoate-grown cells. The early-stationary-phase biphenyl-grown cells showed similar expression of biphenyl pathway genes, but a surprising up-regulation of C(1) metabolic pathway genes was observed. The microarray results were validated by quantitative reverse transcription PCR with a subset of genes of interest. The XeoChips showed a chip-to-chip variation of 13.9%, compared to the 21.6% variation for spotted oligonucleotide microarrays, which is less variation than that typically reported for PCR product microarrays.

Structural flexibility in the Burkholderia mallei genome

The complete genome sequence of Burkholderia mallei ATCC 23344 provides insight into this highly infectious bacterium's pathogenicity and evolutionary history. B. mallei, the etiologic agent of glanders, has come under renewed scientific investigation as a result of recent concerns about its past and potential future use as a biological weapon. Genome analysis identified a number of putative virulence factors whose function was supported by comparative genome hybridization and expression profiling of the bacterium in hamster liver in vivo. The genome contains numerous insertion sequence elements that have mediated extensive deletions and rearrangements of the genome relative to Burkholderia pseudomallei. The genome also contains a vast number (>12,000) of simple sequence repeats. Variation in simple sequence repeats in key genes can provide a mechanism for generating antigenic variation that may account for the mammalian host's inability to mount a durable adaptive immune response to a B. mallei infection.

Genome Properties: a system for the investigation of prokaryotic genetic content for microbiology, genome annotation and comparative genomics

MOTIVATION

The presence or absence of metabolic pathways and structures provide a context that makes protein annotation far more reliable. Compiling such information across microbial genomes improves the functional classification of proteins and provides a valuable resource for comparative genomics.

RESULTS

We have created a Genome Properties system to present key aspects of prokaryotic biology using standardized computational methods and controlled vocabularies. Properties reflect gene content, phenotype, phylogeny and computational analyses. The results of searches using hidden Markov models allow many properties to be deduced automatically, especially for families of proteins (equivalogs) conserved in function since their last common ancestor. Additional properties are derived from curation, published reports and other forms of evidence. Genome Properties system was applied to 156 complete prokaryotic genomes, and is easily mined to find differences between species, correlations between metabolic features and families of uncharacterized proteins, or relationships among properties.

AVAILABILITY

Genome Properties can be found at http://www.tigr.org/Genome_Properties

SUPPLEMENTARY INFORMATION

http://www.tigr.org/tigr-scripts/CMR2/genome_properties_references.spl.

Impact of human cytomegalovirus latent infection on myeloid progenitor cell gene expression

Herpesviruses establish lifelong latent infections in their hosts. Human cytomegalovirus (CMV) targets a population of bone marrow-derived myeloid lineage progenitor cells that serve as a reservoir for reactivation; however, the mechanisms by which latent CMV infection is maintained are unknown. To gain insights into mechanisms of maintenance and reactivation, we employed microarrays of approximately 26,900 sequence-verified human cDNAs to assess global changes in cellular gene expression during experimental CMV latent infection of granulocyte-macrophage progenitors (GM-Ps). This analysis revealed at least 29 host cell genes whose expression was increased and six whose expression was decreased during CMV latency. These changes in transcript levels appeared to be authentic, judging on the basis of further analysis of a subset by semiquantitative reverse transcription-PCR. This study provides a comprehensive snapshot of changes in host cell gene expression that result from latent infection and suggest that CMV regulates genes that encode proteins involved in immunity and host defense, cell growth, signaling, and transcriptional regulation. The host genes whose expression we found altered are likely to contribute to an environment that sustains latent infection.

The genome sequence of the anaerobic, sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough

Desulfovibrio vulgaris Hildenborough is a model organism for studying the energy metabolism of sulfate-reducing bacteria (SRB) and for understanding the economic impacts of SRB, including biocorrosion of metal infrastructure and bioremediation of toxic metal ions. The 3,570,858 base pair (bp) genome sequence reveals a network of novel c-type cytochromes, connecting multiple periplasmic hydrogenases and formate dehydrogenases, as a key feature of its energy metabolism. The relative arrangement of genes encoding enzymes for energy transduction, together with inferred cellular location of the enzymes, provides a basis for proposing an expansion to the 'hydrogen-cycling' model for increasing energy efficiency in this bacterium. Plasmid-encoded functions include modification of cell surface components, nitrogen fixation and a type-III protein secretion system. This genome sequence represents a substantial step toward the elucidation of pathways for reduction (and bioremediation) of pollutants such as uranium and chromium and offers a new starting point for defining this organism's complex anaerobic respiration.

Comparison of the genome of the oral pathogen Treponema denticola with other spirochete genomes

We present the complete 2,843,201-bp genome sequence of Treponema denticola (ATCC 35405) an oral spirochete associated with periodontal disease. Analysis of the T. denticola genome reveals factors mediating coaggregation, cell signaling, stress protection, and other competitive and cooperative measures, consistent with its pathogenic nature and lifestyle within the mixed-species environment of subgingival dental plaque. Comparisons with previously sequenced spirochete genomes revealed specific factors contributing to differences and similarities in spirochete physiology as well as pathogenic potential. The T. denticola genome is considerably larger in size than the genome of the related syphilis-causing spirochete Treponema pallidum. The differences in gene content appear to be attributable to a combination of three phenomena: genome reduction, lineage-specific expansions, and horizontal gene transfer. Genes lost due to reductive evolution appear to be largely involved in metabolism and transport, whereas some of the genes that have arisen due to lineage-specific expansions are implicated in various pathogenic interactions, and genes acquired via horizontal gene transfer are largely phage-related or of unknown function.

High-throughput molecular pathology in human tissues as a method for driving drug discovery

To facilitate prioritization of potential drug targets, gene expression can be localized to individual cell types in normal and diseased tissues. Given the complexity of molecular physiology and pathology, the creation of large-scale molecular pathology databases collating data obtained from human tissues is a challenging marriage of old and new technologies, particularly when considering the many issues that preclude easy access to substantial quantities of human tissues. Molecular pathology databases are powerful tools and are essential for early-stage drug discovery, enabling informed decisions to be made with respect to scientific direction and follow-up research.

A genomewide oscillation in transcription gates DNA replication and cell cycle

Microarray analysis from a yeast continuous synchrony culture system shows a genomewide oscillation in transcription. Maximums in transcript levels occur at three nearly equally spaced intervals in this approximately 40-min cycle of respiration and reduction. Two temporal clusters (4,679 of 5,329) are maximally expressed during the reductive phase of the cycle, whereas a third cluster (650) is maximally expressed during the respiratory phase. Transcription is organized functionally into redox-state superclusters with genes known to be important in respiration or reduction being synthesized in opposite phases of the cycle. The transcriptional cycle gates synchronous bursts in DNA replication in a constant fraction of the population at 40-min intervals. Restriction of DNA synthesis to the reductive phase of the cycle may be an evolutionarily important mechanism for reducing oxidative damage to DNA during replication.

Structural genomics: computational methods for structure analysis

The success of structural genomics initiatives requires the development and application of tools for structure analysis, prediction, and annotation. In this paper we review recent developments in these areas; specifically structure alignment, the detection of remote homologs and analogs, homology modeling and the use of structures to predict function. We also discuss various rationales for structural genomics initiatives. These include the structure-based clustering of sequence space and genome-wide function assignment. It is also argued that structural genomics can be integrated into more traditional biological research if specific biological questions are included in target selection strategies.

Complete genome sequence of the broad-host-range vibriophage KVP40: comparative genomics of a T4-related bacteriophage

The complete genome sequence of the T4-like, broad-host-range vibriophage KVP40 has been determined. The genome sequence is 244,835 bp, with an overall G+C content of 42.6%. It encodes 386 putative protein-encoding open reading frames (CDSs), 30 tRNAs, 33 T4-like late promoters, and 57 potential rho-independent terminators. Overall, 92.1% of the KVP40 genome is coding, with an average CDS size of 587 bp. While 65% of the CDSs were unique to KVP40 and had no known function, the genome sequence and organization show specific regions of extensive conservation with phage T4. At least 99 KVP40 CDSs have homologs in the T4 genome (Blast alignments of 45 to 68% amino acid similarity). The shared CDSs represent 36% of all T4 CDSs but only 26% of those from KVP40. There is extensive representation of the DNA replication, recombination, and repair enzymes as well as the viral capsid and tail structural genes. KVP40 lacks several T4 enzymes involved in host DNA degradation, appears not to synthesize the modified cytosine (hydroxymethyl glucose) present in T-even phages, and lacks group I introns. KVP40 likely utilizes the T4-type sigma-55 late transcription apparatus, but features of early- or middle-mode transcription were not identified. There are 26 CDSs that have no viral homolog, and many did not necessarily originate from Vibrio spp., suggesting an even broader host range for KVP40. From these latter CDSs, an NAD salvage pathway was inferred that appears to be unique among bacteriophages. Features of the KVP40 genome that distinguish it from T4 are presented, as well as those, such as the replication and virion gene clusters, that are substantially conserved.

Classification schemes for protein structure and function

We examine the structural and functional classifications of the protein universe, providing an overview of the existing classification schemes, their features and inter-relationships. We argue that a unified scheme should be based on a natural classification approach and that more comparative analyses of the present schemes are required both to understand their limitations and to help delimit the number of known protein folds and their corresponding functional roles in cells.

Bioverse: Functional, structural and contextual annotation of proteins and proteomes

Functional annotation is routinely performed for large-scale genomics projects and databases. Researchers working on more specific problems, for instance on an individual pathway or complex, also need to be able to quickly, completely and accurately annotate sequences. The Bioverse sequence annotation server (http://bioverse.compbio.washington.edu) provides a web-based interface to allow users to submit protein sequences to the Bioverse framework. Sequences are functionally and structurally annotated and potential contextual annotations are provided. Researchers can also submit candidate genomes for annotation of all proteins encoded by the genome (proteome).

Complete genome sequence of the Q-fever pathogen Coxiella burnetii

The 1,995,275-bp genome of Coxiella burnetii, Nine Mile phase I RSA493, a highly virulent zoonotic pathogen and category B bioterrorism agent, was sequenced by the random shotgun method. This bacterium is an obligate intracellular acidophile that is highly adapted for life within the eukaryotic phagolysosome. Genome analysis revealed many genes with potential roles in adhesion, invasion, intracellular trafficking, host-cell modulation, and detoxification. A previously uncharacterized 13-member family of ankyrin repeat-containing proteins is implicated in the pathogenesis of this organism. Although the lifestyle and parasitic strategies of C. burnetii resemble that of Rickettsiae and Chlamydiae, their genome architectures differ considerably in terms of presence of mobile elements, extent of genome reduction, metabolic capabilities, and transporter profiles. The presence of 83 pseudogenes displays an ongoing process of gene degradation. Unlike other obligate intracellular bacteria, 32 insertion sequences are found dispersed in the chromosome, indicating some plasticity in the C. burnetii genome. These analyses suggest that the obligate intracellular lifestyle of C. burnetii may be a relatively recent innovation.

The Gene Ontology Annotation (GOA) project: implementation of GO in SWISS-PROT, TrEMBL, and InterPro

Gene Ontology Annotation (GOA) is a project run by the European Bioinformatics Institute (EBI) that aims to provide assignments of terms from the Gene Ontology (GO) resource to gene products in a number of its databases (http://www.ebi.ac.uk/GOA). In the first stage of this project, GO assignments have been applied to a data set representing the complete human proteome by a combination of electronic mappings and manual curation. This vocabulary has also been applied to the nonredundant proteome sets for all other completely sequenced organisms as well as to proteins from a wide range of organisms where the proteome is not yet complete.

Bioinformatics in the post-sequence era

In the past decade, bioinformatics has become an integral part of research and development in the biomedical sciences. Bioinformatics now has an essential role both in deciphering genomic, transcriptomic and proteomic data generated by high-throughput experimental technologies and in organizing information gathered from traditional biology. Sequence-based methods of analyzing individual genes or proteins have been elaborated and expanded, and methods have been developed for analyzing large numbers of genes or proteins simultaneously, such as in the identification of clusters of related genes and networks of interacting proteins. With the complete genome sequences for an increasing number of organisms at hand, bioinformatics is beginning to provide both conceptual bases and practical methods for detecting systemic functional behaviors of the cell and the organism.

Domains, motifs and clusters in the protein universe

The rapid growth of bio-sequence information has resulted in an increasing demand for reliable methods that group proteins. A few databases with curated alignments of protein families have demonstrated that expert-driven repositories can keep up with the data deluge in the genome era. These original resources implicitly identify domain-like modules in proteins. An increasing number of automatic methods have sprouted over the past few years that cluster the protein universe. Many of these implicitly dissect proteins into structural domain-like fragments. In a very coarse-grained evaluation, some of the automatic methods appear to be on par with expert-driven approaches. However, neither automatic nor manual methods are currently entirely up to the challenges of tasks such as target selection in structural genomics. Thus, we urgently need refined and sustained automatic clustering tools.

Automated annotation of microbial proteomes in SWISS-PROT

Large-scale sequencing of prokaryotic genomes demands the automation of certain annotation tasks currently manually performed in the production of the SWISS-PROT protein knowledgebase. The HAMAP project, or 'High-quality Automated and Manual Annotation of microbial Proteomes', aims to integrate manual and automatic annotation methods in order to enhance the speed of the curation process while preserving the quality of the database annotation. Automatic annotation is only applied to entries that belong to manually defined orthologous families and to entries with no identifiable similarities (ORFans). Many checks are enforced in order to prevent the propagation of wrong annotation and to spot problematic cases, which are channelled to manual curation. The results of this annotation are integrated in SWISS-PROT, and a website is provided at http://www.expasy.org/sprot/hamap/.

Protein family classification and functional annotation

With the accelerated accumulation of genomic sequence data, there is a pressing need to develop computational methods and advanced bioinformatics infrastructure for reliable and large-scale protein annotation and biological knowledge discovery. The Protein Information Resource (PIR) provides an integrated public resource of protein informatics to support genomic and proteomic research. PIR produces the Protein Sequence Database of functionally annotated protein sequences. The annotation problems are addressed by a classification-driven and rule-based method with evidence attribution, coupled with an integrated knowledge base system being developed. The approach allows sensitive identification, consistent and rich annotation, and systematic detection of annotation errors, as well as distinction of experimentally verified and computationally predicted features. The knowledge base consists of two new databases, sequence analysis tools, and graphical interfaces. PIR-NREF, a non-redundant reference database, provides a timely and comprehensive collection of all protein sequences, totaling more than 1,000,000 entries. iProClass, an integrated database of protein family, function, and structure information, provides extensive value-added features for about 830,000 proteins with rich links to over 50 molecular databases. This paper describes our approach to protein functional annotation with case studies and examines common identification errors. It also illustrates that data integration in PIR supports exploration of protein relationships and may reveal protein functional associations beyond sequence homology.

Target selection and determination of function in structural genomics

The first crucial step in any structural genomics project is the selection and prioritization of target proteins for structure determination. There may be a number of selection criteria to be satisfied, including that the proteins have novel folds, that they be representatives of large families for which no structure is known, and so on. The better the selection at this stage, the greater is the value of the structures obtained at the end of the experimental process. This value can be further enhanced once the protein structures have been solved if the functions of the given proteins can also be determined. Here we describe the methods used at either end of the experimental process: firstly, sensitive sequence comparison techniques for selecting a high-quality list of target proteins, and secondly the various computational methods that can be applied to the eventual 3D structures to determine the most likely biochemical function of the proteins in question.

The TIGRFAMs database of protein families

TIGRFAMs is a collection of manually curated protein families consisting of hidden Markov models (HMMs), multiple sequence alignments, commentary, Gene Ontology (GO) assignments, literature references and pointers to related TIGRFAMs, Pfam and InterPro models. These models are designed to support both automated and manually curated annotation of genomes. TIGRFAMs contains models of full-length proteins and shorter regions at the levels of superfamilies, subfamilies and equivalogs, where equivalogs are sets of homologous proteins conserved with respect to function since their last common ancestor. The scope of each model is set by raising or lowering cutoff scores and choosing members of the seed alignment to group proteins sharing specific function (equivalog) or more general properties. The overall goal is to provide information with maximum utility for the annotation process. TIGRFAMs is thus complementary to Pfam, whose models typically achieve broad coverage across distant homologs but end at the boundaries of conserved structural domains. The database currently contains over 1600 protein families. TIGRFAMs is available for searching or downloading at www.tigr.org/TIGRFAMs.

The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003

The SWISS-PROT protein knowledgebase (http://www.expasy.org/sprot/ and http://www.ebi.ac.uk/swissprot/) connects amino acid sequences with the current knowledge in the Life Sciences. Each protein entry provides an interdisciplinary overview of relevant information by bringing together experimental results, computed features and sometimes even contradictory conclusions. Detailed expertise that goes beyond the scope of SWISS-PROT is made available via direct links to specialised databases. SWISS-PROT provides annotated entries for all species, but concentrates on the annotation of entries from human (the HPI project) and other model organisms to ensure the presence of high quality annotation for representative members of all protein families. Part of the annotation can be transferred to other family members, as is already done for microbes by the High-quality Automated and Manual Annotation of microbial Proteomes (HAMAP) project. Protein families and groups of proteins are regularly reviewed to keep up with current scientific findings. Complementarily, TrEMBL strives to comprise all protein sequences that are not yet represented in SWISS-PROT, by incorporating a perpetually increasing level of mostly automated annotation. Researchers are welcome to contribute their knowledge to the scientific community by submitting relevant findings to SWISS-PROT at swiss-prot@expasy.org.

MBGD: microbial genome database for comparative analysis

MBGD is a workbench system for comparative analysis of completely sequenced microbial genomes. The central function of MBGD is to create an orthologous gene classification table using precomputed all-against-all similarity relationships among genes in multiple genomes. In MBGD, an automated classification algorithm has been implemented so that users can create their own classification table by specifying a set of organisms and parameters. This feature is especially useful when the user's interest is focused on some taxonomically related organisms. The created classification table is stored into the database and can be explored combining with the data of individual genomes as well as similarity relationships among genomes. Using these data, users can carry out comparative analyses from various points of view, such as phylogenetic pattern analysis, gene order comparison and detailed gene structure comparison. MBGD is accessible at http://mbgd.genome.ad.jp/.

The European Bioinformatics Institute's data resources

As the amount of biological data grows, so does the need for biologists to store and access this information in central repositories in a free and unambiguous manner. The European Bioinformatics Institute (EBI) hosts six core databases, which store information on DNA sequences (EMBL-Bank), protein sequences (SWISS-PROT and TrEMBL), protein structure (MSD), whole genomes (Ensembl) and gene expression (ArrayExpress). But just as a cell would be useless if it couldn't transcribe DNA or translate RNA, our resources would be compromised if each existed in isolation. We have therefore developed a range of tools that not only facilitate the deposition and retrieval of biological information, but also allow users to carry out searches that reflect the interconnectedness of biological information. The EBI's databases and tools are all available on our website at www.ebi.ac.uk.

The InterPro Database, 2003 brings increased coverage and new features

InterPro, an integrated documentation resource of protein families, domains and functional sites, was created in 1999 as a means of amalgamating the major protein signature databases into one comprehensive resource. PROSITE, Pfam, PRINTS, ProDom, SMART and TIGRFAMs have been manually integrated and curated and are available in InterPro for text- and sequence-based searching. The results are provided in a single format that rationalises the results that would be obtained by searching the member databases individually. The latest release of InterPro contains 5629 entries describing 4280 families, 1239 domains, 95 repeats and 15 post-translational modifications. Currently, the combined signatures in InterPro cover more than 74% of all proteins in SWISS-PROT and TrEMBL, an increase of nearly 15% since the inception of InterPro. New features of the database include improved searching capabilities and enhanced graphical user interfaces for visualisation of the data. The database is available via a webserver (http://www.ebi.ac.uk/interpro) and anonymous FTP (ftp://ftp.ebi.ac.uk/pub/databases/interpro).