Ivica Tamas

Automatic clustering of orthologs and inparalogs shared by multiple proteomes

MOTIVATION The complete sequencing of many genomes has made it possible to identify orthologous genes descending from a common ancestor. However, reconstruction of evolutionary history over long time periods faces many challenges due to gene duplications and losses. Identification of orthologous groups shared by multiple proteomes therefore becomes a clustering problem in which an optimal compromise between conflicting evidences needs to be found. RESULTS Here we present a new proteome-scale analysis program called MultiParanoid that can automatically find orthology relationships between proteins in multiple proteomes. The software is an extension of the InParanoid program that identifies orthologs and inparalogs in pairwise proteome comparisons. MultiParanoid applies a clustering algorithm to merge multiple pairwise ortholog groups from InParanoid into multi-species ortholog groups. To avoid outparalogs in the same cluster, MultiParanoid only combines species that share the same last ancestor. To validate the clustering technique, we compared the results to a reference set obtained by manual phylogenetic analysis. We further compared the results to ortholog groups in KOGs and OrthoMCL, which revealed that MultiParanoid produces substantially fewer outparalogs than these resources. AVAILABILITY MultiParanoid is a freely available standalone program that enables efficient orthology analysis much needed in the post-genomic era. A web-based service providing access to the original datasets, the resulting groups of orthologs, and the source code of the program can be found at http://multiparanoid.cgb.ki.se.

Endosymbiont gene functions impaired and rescued by polymerase infidelity at poly(A) tracts

Among host-dependent bacteria that have evolved by extreme reductive genome evolution, long-term bacterial endosymbionts of insects have the smallest (160–790 kb) and most A + T-rich (>70%) bacterial genomes known to date. These genomes are riddled with poly(A) tracts, and 5–50% of genes contain tracts of 10 As or more. Here, we demonstrate transcriptional slippage at poly(A) tracts within genes of Buchnera aphidicola associated with aphids and Blochmannia pennsylvanicus associated with ants. Several tracts contain single frameshift deletions; these apparent pseudogenes showed patterns of constraint consistent with purifying selection on the encoded proteins. Transcriptional slippage yielded a heterogeneous population of transcripts with variable numbers of As in the tract. Across several frameshifted genes, including B. aphidicola cell wall biosynthesis genes and a B. pennsylvanicus histidine biosynthesis gene, 12–50% of transcripts contained corrected reading frames that could potentially yield full-length proteins. In situ immunostaining confirmed the production of the cell wall biosynthetic enzyme UDP-N-acetylmuramyl pentapeptide synthase encoded by the frameshifted murF gene. Simulation studies indicated an overrepresentation of poly(A) tracts in endosymbiont genomes relative to other A + T-rich bacterial genomes. Polymerase infidelity at poly(A) tracts rescues the functionality of genes with frameshift mutations and, conversely, reduces the efficiency of expression for in-frame genes carrying poly(A) regions. These features of homopolymeric tracts could be exploited to manipulate gene expression in small synthetic genomes.

Sequencing of the Francisella tularensis strain Schu 4 genome reveals the shikimate and purine metabolic pathways, targets for the construction of a rationally attenuated auxotrophic vaccine.

Francisella tularensis is the etiological agent of tularemia, a serious disease in several Northern hemisphere countries. The organism has fastidious growth requirements and is very poorly understood at the genetic and molecular levels. Given the lack of data on this organism, we undertook the sample sequencing of its genome. A random library of DNA fragments from a highly virulent strain (Schu 4) of F. tularensis was constructed and the nucleotide sequences of 13,904 cloned fragments were determined and assembled into 353 contigs. A total of 1.83 Mb of nucleotide sequence was obtained that had a G+C content of 33.2%. Genes located on plasmids pOM1 and pNFL10, which had been previously isolated from low virulence strains of F. tularensis, were absent but all of the other known F. tularensis genes were represented in the assembled data. F. tularensis Schu4 was able to grow in the absence of aromatic amino acids and orthologues of genes which could encode enzymes in the shikimate pathway in other bacteria were identified in the assembled data. Genes that could encode all of the enzymes in the purine biosynthetic and most of the en- zymes in the purine salvage pathways were also identified. This data will be used to develop defined rationally attenuated mutants of F. tularensis, which could be used as replacements for the existing genetically undefined live vaccine strain.

Preliminary analysis and annotation of the partial genome sequence of Francisella tularensis strain Schu 4.

Preliminary analysis and annotation of the partial genome sequence of Francisella tularensis strain Schu 4.

Mutualists and parasites: how to paint yourself into a (metabolic) corner

Eukaryotes have developed an elaborate series of interactions with bacteria that enter their bodies and/or cells. Genome evolution of symbiotic and parasitic bacteria multiplying inside eukaryotic cells results in both convergent and divergent changes. The genome sequences of the symbiotic bacteria of aphids, Buchnera aphidicola, and the parasitic bacteria of body louse and humans, Rickettsia prowazekii, provide insights into these processes. Convergent genome characteristics include reduction in genome sizes and lowered G+C content values. Divergent evolution was recorded for amino acid and cell wall biosynthetic genes. The presence of pseudogenes in both genomes provides examples of recent gene inactivation events and offers clues to the process of genome deterioration and host‐cell adaptation.

Determinative value of a portion of the nifH sequence for the genera Nostoc and Anabaena (Cyanobacteria)

The taxonomic positions of Nostoc and Anabaena strains are currently disputed. We selected three Nostoc and Anabaena strains, using the classic criteria of morphology and life cycle. DNA sequences of a part of the nifH gene were determined from these strains and aligned with homologous sequences from 10 other Nostoc/Anabaena strains in the public databases. Phylogenetic reconstructions were carried out to test the consistency of the taxonomic placement of these strains. The phylogenetic trees do not separate these strains into distinct groups. Our results are in agreement with other molecular-based phylogenies that also fail to differentiate the Nostoc-Anabaena groups. The data suggest that the currently recognized genera Nostoc and Anabaena may in fact belong within a single, broadly defined genus.

Fisher: a program for the detection of H/ACA snoRNAs using MFE secondary structure prediction and comparative genomics -- assessment and update.

BackgroundThe H/ACA family of small nucleolar RNAs (snoRNAs) plays a central role in guiding the pseudouridylation of ribosomal RNA (rRNA). In an effort to systematically identify the complete set of rRNA-modifying H/ACA snoRNAs from the genome sequence of the budding yeast, Saccharomyces cerevisiae, we developed a program – Fisher – and previously presented several candidate snoRNAs based on our analysis [1].FindingsIn this report, we provide a brief update of this work, which was aborted after the publication of experimentally-identified snoRNAs [2] identical to candidates we had identified bioinformatically using Fisher. Our motivation for revisiting this work is to report on the status of the candidate snoRNAs described in [1], and secondly, to report that a modified version of Fisher together with the available multiple yeast genome sequences was able to correctly identify several H/ACA snoRNAs for modification sites not identified by the snoGPS program [3]. While we are no longer developing Fisher, we briefly consider the merits of the Fisher algorithm relative to snoGPS, which may be of use for workers considering pursuing a similar search strategy for the identification of small RNAs. The modified source code for Fisher is made available as supplementary material.ConclusionOur results confirm the validity of using minimum free energy (MFE) secondary structure prediction to guide comparative genomic screening for RNA families with few sequence constraints.