Search

Background: Orthology is a central tenet of comparative genomics and ortholog identification is instrumental to protein function prediction. Major advances have been made to determine orthology relations among a set of homologous proteins. However, they depend on the comparison of individual sequences and do not take into account divergent orthologs., Results: We have developed an iterative orthology prediction method, Ortho-Profile, that uses reciprocal best hits at the level of sequence profiles to infer orthology. It increases ortholog detection by 20% compared to sequence-to-sequence comparisons. Ortho-Profile predicts 598 human orthologs of mitochondrial proteins from Saccharomyces cerevisiae and Schizosaccharomyces pombe with 94% accuracy. Of these, 181 were not known to localize to mitochondria in mammals. Among the predictions of the Ortho-Profile method are 11 human cytochrome c oxidase (COX) assembly proteins that are implicated in mitochondrial function and disease. Their co-expression patterns, experimentally verified subcellular localization, and co-purification with human COX-associated proteins support these predictions. For the human gene C12orf62, the ortholog of S. cerevisiae COX14, we specifically confirm its role in negative regulation of the translation of cytochrome c oxidase., Conclusions: Divergent homologs can often only be detected by comparing sequence profiles and profile-based hidden Markov models. The Ortho-Profile method takes advantage of these techniques in the quest for orthologs.

Dendritic cells (DC) take up pathogens through phagocytosis and process them into protein and lipid fragments for presentation to T cells. So far, the proteome of the human DC phagosome, a detrimental compartment for antigen processing and presentation as well as for DC activation, remains largely uncharacterized. Here we have analyzed the protein composition of phagosomes from human monocyte-derived DC. For LC-MS/MS analysis we purified phagosomes from DC using latex beads targeted to DC-SIGN, and quantified proteins using a label-free method. We used organellar enrichment ranking (OER) to select proteins with a high potential to be relevant for phagosome function. The method compares phagosome protein abundance with protein abundance in whole DC. Phagosome enrichment indicates specific recruitment to the phagosome rather than co-purification or passive incorporation. Using OER we extracted the most enriched proteins that we further complemented with functionally associated proteins to define a set of 90 phagosomal proteins that included many proteins with established relevance on DC phagosomes as well as high potential novel candidates. We already experimentally confirmed phagosomal recruitment of Galectin-9, which has not been previously associated with phagocytosis, to both bead and pathogen containing phagosomes, suggesting a role for Galectin-9 in DC phagocytosis., (Copyright © 2011 Elsevier B.V. All rights reserved.)

We report a fragmented mitochondrial network and swollen and irregularly shaped mitochondria with partial to complete loss of the cristae in fibroblasts of a patient with a novel TMEM70 gene deletion, which could be completely restored by complementation of the TMEM70 genetic defect. Comparative genomics analysis predicted the topology of TMEM70 in the inner mitochondrial membrane, which could be confirmed by immunogold labeling experiments, and showed that the TMEM70 gene is not restricted to higher multi-cellular eukaryotes. This study demonstrates that the role of complex V in mitochondrial cristae morphology applies to human mitochondrial disease pathology., (Copyright © 2011 Elsevier B.V. and Mitochondria Research Society. All rights reserved. All rights reserved.)

Motivation: The intensification of DNA sequencing will increasingly unveil uncharacterized species with potential alternative genetic codes. A total of 0.65% of the DNA sequences currently in Genbank encode their proteins with a variant genetic code, and these exceptions occur in many unrelated taxa., Results: We introduce FACIL (Fast and Accurate genetic Code Inference and Logo), a fast and reliable tool to evaluate nucleic acid sequences for their genetic code that detects alternative codes even in species distantly related to known organisms. To illustrate this, we apply FACIL to a set of mitochondrial genomic contigs of Globobulimina pseudospinescens. This foraminifer does not have any sequenced close relative in the databases, yet we infer its alternative genetic code with high confidence values. Results are intuitively visualized in a Genetic Code Logo., Availability and Implementation: FACIL is available as a web-based service at http://www.cmbi.ru.nl/FACIL/ and as a stand-alone program.

The assembly of mitochondrial respiratory chain complex IV (cytochrome c oxidase) involves the coordinated action of several assembly chaperones. In Saccharomyces cerevisiae, at least 30 different assembly chaperones have been identified. To date, pathogenic mutations leading to a mitochondrial disorder have been identified in only seven of the corresponding human genes. One of the genes for which the relevance to human pathology is unknown is C2orf64, an ortholog of the S. cerevisiae gene PET191. This gene has previously been shown to be a complex IV assembly factor in yeast, although its exact role is still unknown. Previous research in a large cohort of complex IV deficient patients did not support an etiological role of C2orf64 in complex IV deficiency. In this report, a homozygous mutation in C2orf64 is described in two siblings affected by fatal neonatal cardiomyopathy. Pathogenicity of the mutation is supported by the results of a complementation experiment, showing that complex IV activity can be fully restored by retroviral transduction of wild-type C2orf64 in patient-derived fibroblasts. Detailed analysis of complex IV assembly intermediates in patient fibroblasts by 2D-BN PAGE revealed the accumulation of a small assembly intermediate containing subunit COX1 but not the COX2, COX4, or COX5b subunits, indicating that C2orf64 is involved in an early step of the complex IV assembly process. The results of this study demonstrate that C2orf64 is essential for human complex IV assembly and that C2orf64 mutational analysis should be considered for complex IV deficient patients, in particular those with hypertrophic cardiomyopathy., (Copyright © 2011 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Complex I (NADH:ubiquinone oxidoreductase) is the first and largest protein complex of the oxidative phosphorylation. Crystal structures have elucidated the positions of most subunits of bacterial evolutionary origin in the complex, but the positions of the eukaryotic subunits are unknown. Based on the analysis of sequence conservation we propose intra-molecular disulfide bridges and the inter-membrane space localization of three Cx(9)C-containing subunits in human: NDUFS5, NDUFB7 and NDUFA8. We experimentally confirm the localization of the latter two, while our data are consistent with disulfide bridges in NDUFA8. We propose these subunits stabilize the membrane domain of complex I., (Copyright © 2011 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.)

Although the origin of mitochondria from the endosymbiosis of an α-proteobacterium is well established, the nature of the host cell, the metabolic complexity of the endosymbiont and the subsequent evolution of the proto-mitochondrion into all its current appearances are still the subject of discovery and sometimes debate. Here we review what has been inferred about the original composition and subsequent evolution of the mitochondrial proteome and essential mitochondrial systems. The evolutionary mosaic that currently constitutes mitochondrial proteomes contains (i) endosymbiotic proteins (15-45%), (ii) proteins without detectable orthologs outside the eukaryotic lineage (40%), and (iii) proteins that are derived from non-proteobacterial Bacteria, Bacteriophages and Archaea (15%, specifically multiple tRNA-modification proteins). Protein complexes are of endosymbiotic origin, but have greatly expanded with novel eukaryotic proteins; in contrast to mitochondrial enzymes that are both of proteobacterial and non-proteobacterial origin. This disparity is consistent with the complexity hypothesis, which argues that proteins that are a part of large, multi-subunit complexes are unlikely to undergo horizontal gene transfer. We observe that they neither change their subcellular compartments in the course of evolution, even when their genes do., (Copyright © 2010 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Mitochondria are an essential organelle, not only to the human cell, but to all eukaryotic life. This essentiality is reflected in the large number of mutations in genes encoding mitochondrial proteins that lead to disease. Aside from their relevance to disease, mitochondria are, given their endosymbiotic origin, very interesting from an evolutionary point of view. Here, in the year that marks the bicentenary of Darwin's birth and the 150th anniversary of the publication of "On the origin of species" we review approaches that implicitly or explicitly use evolutionary analyses to find new genes involved in mitochondrial disease and to predict their function and involvement in pathways. We show how the phenotypic spectrum of mitochondrial disease is linked to the evolutionary origin of mitochondrial proteins, how combinations of evolutionary data and genomics data have been used to predict the mitochondrial proteome and functional links between the mitochondrial proteins and how the evolution of the mitochondrial proteome has been used to predict new mitochondrial disease genes. For the latter we review and reanalyze the eukaryotic evolution of the NADH:ubiquinone oxidoreductase (complex I) and the proteins involved in its assembly.

Background: Mitochondria are highly complex, membrane-enclosed organelles that are essential to the eukaryotic cell. The experimental elucidation of organellar proteomes combined with the sequencing of complete genomes allows us to trace the evolution of the mitochondrial proteome., Results: We present a systematic analysis of the evolution of mitochondria via gene duplication in the human lineage. The most common duplications are intra-mitochondrial, in which the ancestral gene and the daughter genes encode mitochondrial proteins. These duplications significantly expanded carbohydrate metabolism, the protein import machinery and the calcium regulation of mitochondrial activity. The second most prevalent duplication, inter-compartmental, extended the catalytic as well as the RNA processing repertoire by the novel mitochondrial localization of the protein encoded by one of the daughter genes. Evaluation of the phylogenetic distribution of N-terminal targeting signals suggests a prompt gain of the novel localization after inter-compartmental duplication. Relocalized duplicates are more often expressed in a tissue-specific manner relative to intra-mitochondrial duplicates and mitochondrial proteins in general. In a number of cases, inter-compartmental duplications can be observed in parallel in yeast and human lineages leading to the convergent evolution of subcellular compartments., Conclusions: One-to-one human-yeast orthologs are typically restricted to their ancestral subcellular localization. Gene duplication relaxes this constraint on the cellular location, allowing nascent proteins to be relocalized to other compartments. We estimate that the mitochondrial proteome expanded at least 50% since the common ancestor of human and yeast.

Background: Thanks to recent high coverage mass-spectrometry studies and reconstructed protein complexes, we are now in an unprecedented position to study the evolution of biological systems. Gene duplications, known to be a major source of innovation in evolution, can now be readily examined in the context of protein complexes., Results: We observe that paralogs operating in the same complex fulfill different roles: mRNA dosage increase for more than a hundred cytosolic ribosomal proteins, mutually exclusive participation of at least 54 paralogs resulting in alternative forms of complexes, and 24 proteins contributing to bona fide structural growth. Inspection of paralogous proteins participating in two independent complexes shows that an ancient, pre-duplication protein functioned in both multi-protein assemblies and a gene duplication event allowed the respective copies to specialize and split their roles., Conclusion: Variants with conditionally assembled, paralogous subunits likely have played a role in yeast's adaptation to anaerobic conditions. In a number of cases the gene duplication has given rise to one duplicate that is no longer part of a protein complex and shows an accelerated rate of evolution. Such genes could provide the raw material for the evolution of new functions.

Coding of multiple proteins by overlapping reading frames is not a feature one would associate with eukaryotic genes. Indeed, codependency between codons of overlapping protein-coding regions imposes a unique set of evolutionary constraints, making it a costly arrangement. Yet in cases of tightly coexpressed interacting proteins, dual coding may be advantageous. Here we show that although dual coding is nearly impossible by chance, a number of human transcripts contain overlapping coding regions. Using newly developed statistical techniques, we identified 40 candidate genes with evolutionarily conserved overlapping coding regions. Because our approach is conservative, we expect mammals to possess more dual-coding genes. Our results emphasize that the skepticism surrounding eukaryotic dual coding is unwarranted: rather than being artifacts, overlapping reading frames are often hallmarks of fascinating biology.

Motivation: The accumulation of genome sequences will only accelerate in the coming years. We aim to use this abundance of data to improve the quality of genomic alignments and devise a method which is capable of detecting regions evolving under weak or no evolutionary constraints., Results: We describe a genome alignment program AuberGene, which explores the idea of transitivity of local alignments. Assessment of the program was done based on a 2 Mbp genomic region containing the CFTR gene of 13 species. In this region, we can identify 53% of human sequence sharing common ancestry with mouse, as compared with 44% found using the usual pairwise alignment. Between human and tetraodon 93 orthologous exons are found, as compared with 77 detected by the pairwise human-tetraodon comparison. AuberGene allows the user to (1) identify distant, previously undetected, conserved orthogonal regions such as ORFs or regulatory regions; (2) identify neutrally evolving regions in related species which are often overlooked by other alignment programs; (3) recognize false orthologous genomic regions. The increased sensitivity of the method is not obtained at the cost of reduced specificity. Our results suggest that, over the CFTR region, human shares 10% more sequence with mouse than previously thought ( approximately 50%, instead of 40% found with the pairwise alignment).

Motivation: Internal repeats in coding sequences correspond to structural and functional units of proteins. Moreover, duplication of fragments of coding sequences is known to be a mechanism to facilitate evolution. Identification of repeats is crucial to shed light on the function and structure of proteins, and explain their evolutionary past. The task is difficult because during the course of evolution many repeats diverged beyond recognition., Results: We introduce a new method TRUST, for ab initio determination of internal repeats in proteins. It provides an improvement in prediction quality as compared to alternative state-of-the-art methods. The increased sensitivity and accuracy of the method is achieved by exploiting the concept of transitivity of alignments. Starting from significant local suboptimal alignments, the application of transitivity allows us to (1) identify distant repeat homologues for which no alignments were found; (2) gain confidence about consistently well-aligned regions; and (3) recognize and reduce the contribution of non-homologous repeats. This re-assessment step enables us to derive a virtually noise-free profile representing a generalized repeat with high fidelity. We also obtained superior specificity by employing rigid statistical testing for self-sequence and profile-sequence alignments. Assessment was done using a database of repeat annotations based on structural superpositioning. The results show that TRUST is a useful and reliable tool for mining tandem and non-tandem repeats in protein sequence databases, capable of predicting multiple repeat types with varying intervening segments within a single sequence., Availability: The TRUST server (together with the source code) is available at http://ibivu.cs.vu.nl/programs/trustwww

We present a model of contextual alignment of biological sequences. It is an extension of the classical alignment, in which we assume that the cost of a substitution depends on the surrounding symbols. In this model the cost of transforming one sequence into another depends on the order of editing operations. We present efficient algorithms for calculating this cost, as well as reconstructing (the representation of) all the orders of operations which yield this optimal cost. A precise characterization of the families of linear orders which can emerge this way is given.