Genome sequence and gene compaction of the eukaryote parasite Encephalitozoon cuniculi

ATP/ADP translocases: a common feature of obligate intracellular amoebal symbionts related to Chlamydiae and Rickettsiae

ATP/ADP translocases catalyze the highly specific transport of ATP across a membrane in an exchange mode with ADP. Such unique transport proteins are employed by plant plastids and have among the prokaryotes so far only been identified in few obligate intracellular bacteria belonging to the Chlamydiales and the Rickettsiales. In this study, 12 phylogenetically diverse bacterial endosymbionts of free-living amoebae and paramecia were screened for the presence of genes encoding ATP/ADP transport proteins. The occurrence of ATP/ADP translocase genes was found to be restricted to endosymbionts related to rickettsiae and chlamydiae. We showed that the ATP/ADP transport protein of the Parachlamydia-related endosymbiont of Acanthamoeba sp. strain UWE25, a recently identified relative of the important human pathogens Chlamydia trachomatis and Chlamydophila pneumoniae, is functional when expressed in the heterologous host Escherichia coli and demonstrated the presence of transcripts during the chlamydial developmental cycle. These findings indicate that the interaction between Parachlamydia-related endosymbionts and their amoeba hosts concerns energy parasitism similar to the interaction between pathogenic chlamydiae and their human host cells. Phylogenetic analysis of all known ATP/ADP translocases indicated that the genes encoding ATP/ADP translocases originated from a chlamydial ancestor and were, after an ancient gene duplication, transferred horizontally to rickettsiae and plants.

Guanosine Distribution and Oxidation Resistance in Eight Eukaryotic Genomes

Reactive oxygen species that attack DNA are continuously generated in living cells. Both the guanosine (G) mole fraction and its distribution should affect the stability of genomes and their parts to oxidation. At a lesser G content, genomes should be more oxidation resistant or "ennobled". Oxidant scavenging by G's in nonessential parts of introns and intergenic domains should decrease G oxidation in the essential exons. To determine whether genomes are indeed ennobled and whether oxidant-scavenging domains exist in genomes, the relative rates of guanosine oxidation in average exons, introns, and intergenic domains were estimated. Comparison among genomes indicated that average exons are ennobled in the genomes of Caenorhabditis (worm), Arabidopsis (plant), Saccharomyces (yeast), Schizosaccharomyces (yeast), and Plasmodium (malaria parasite), and that average introns and intergenic domains are ennobled in these genomes and in the genome of Drosophila (fly). The exon oxidation rates estimated for these genomes were less than the rate for the hypothetical "standard" genome, with a 0.25 mole fraction of uniformly distributed G. For Plasmodium the rate was half of that estimated for the standard genome. Average exons were not ennobled in the human or fly genomes; their G distributions were comparable to that in the standard genome. Instead, their exons were situated between introns and intergenic domains that could protect them by oxidant scavenging, the G's of their introns and intergenic domains outnumbering those of their exons 50-fold in humans and 4-fold in flies. The G distribution in the Encephalitozoon (parasite) genome was not protective relative to that of the standard genome.

The positive and negative regulation of Tn10 transposition by IHF is mediated by structurally asymmetric transposon arms

The Tn10 transpososome has symmetrical components on either side: there are two transposon ends each of which has binding sites for a monomer of transposase and an IHF heterodimer. The DNA bending activity of IHF stimulates assembly of an intermediate with tightly folded transposon ends in which transposase has additional 'subterminal' DNA contacts, located distal to the IHF site. These subterminal contacts are required to activate later steps in the reaction. Quantitative hydroxyl radical footprinting and gel retardation unfolding experiments show that the transpososome is fundamentally asymmetric, despite having identical components on either side. Major differences between the transposon ends define alpha and beta sides of the complex. IHF can dissociate from the transposon arm on the beta side of the complex in the absence of metal ion. However, IHF is locked onto the alpha side of the complex, probably by the subterminal transposase contacts, until released by a metal ion-dependent conformational change. Later in the reaction, IHF inhibits target interactions. Using a very short transposon arm, target interactions are demonstrated at a saturating IHF concentration. This suggests that inhibition of target interactions is due to steric hindrance of the target binding site by a single IHF-folded transposon arm.

From a single double helix to paired double helices and back

The propagation of our genomes during cell proliferation depends on the movement of sister DNA molecules produced by DNA replication to opposite sides of the cell before it divides. This feat is achieved by microtubules in eukaryotic cells but it has long remained a mystery how cells ensure that sister DNAs attach to microtubules with opposite orientations, known as amphitelic attachment. It is currently thought that sister chromatid cohesion has a crucial role. By resisting the forces exerted by microtubules, sister chromatid cohesion gives rise to tension that is thought essential for stabilizing kinetochore-microtubule attachments. Efficient amphitelic attachment is therefore achieved by an error correction mechanism that selectively eliminates connections that do not give rise to tension. Cohesion between sister chromatids is mediated by a multisubunit complex called cohesin which forms a gigantic ring structure. It has been proposed that sister DNAs are held together owing to their becoming entrapped within a single cohesin ring. Cohesion between sister chromatids is destroyed at the metaphase to anaphase transition by proteolytic cleavage of cohesin's Scc1 subunit by a thiol protease called separase, which severs the ring and thereby releases sister DNAs.

A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes

We examined functional and evolutionary patterns in the recently constructed set of 5,873 clusters of predicted orthologs from seven eukaryotic genomes. The analysis reveals a conserved core of largely essential eukaryotic genes as well as major diversification and innovation associated with evolution of eukaryotic genomes.

Background

Sequencing the genomes of multiple, taxonomically diverse eukaryotes enables in-depth comparative-genomic analysis which is expected to help in reconstructing ancestral eukaryotic genomes and major events in eukaryotic evolution and in making functional predictions for currently uncharacterized conserved genes.

Results

We examined functional and evolutionary patterns in the recently constructed set of 5,873 clusters of predicted orthologs (eukaryotic orthologous groups or KOGs) from seven eukaryotic genomes: Caenorhabditis elegans, Drosophila melanogaster, Homo sapiens, Arabidopsis thaliana, Saccharomyces cerevisiae, Schizosaccharomyces pombe and Encephalitozoon cuniculi. Conservation of KOGs through the phyletic range of eukaryotes strongly correlates with their functions and with the effect of gene knockout on the organism's viability. The approximately 40% of KOGs that are represented in six or seven species are enriched in proteins responsible for housekeeping functions, particularly translation and RNA processing. These conserved KOGs are often essential for survival and might approximate the minimal set of essential eukaryotic genes. The 131 single-member, pan-eukaryotic KOGs we identified were examined in detail. For around 20 that remained uncharacterized, functions were predicted by in-depth sequence analysis and examination of genomic context. Nearly all these proteins are subunits of known or predicted multiprotein complexes, in agreement with the balance hypothesis of evolution of gene copy number. Other KOGs show a variety of phyletic patterns, which points to major contributions of lineage-specific gene loss and the 'invention' of genes new to eukaryotic evolution. Examination of the sets of KOGs lost in individual lineages reveals co-elimination of functionally connected genes. Parsimonious scenarios of eukaryotic genome evolution and gene sets for ancestral eukaryotic forms were reconstructed. The gene set of the last common ancestor of the crown group consists of 3,413 KOGs and largely includes proteins involved in genome replication and expression, and central metabolism. Only 44% of the KOGs, mostly from the reconstructed gene set of the last common ancestor of the crown group, have detectable homologs in prokaryotes; the remainder apparently evolved via duplication with divergence and invention of new genes.

Conclusions

The KOG analysis reveals a conserved core of largely essential eukaryotic genes as well as major diversification and innovation associated with evolution of eukaryotic genomes. The results provide quantitative support for major trends of eukaryotic evolution noticed previously at the qualitative level and a basis for detailed reconstruction of evolution of eukaryotic genomes and biology of ancestral forms.

Identification and comparative analysis of components from the signal recognition particle in protozoa and fungi

Background

The signal recognition particle (SRP) is a ribonucleoprotein complex responsible for targeting proteins to the ER membrane. The SRP of metazoans is well characterized and composed of an RNA molecule and six polypeptides. The particle is organized into the S and Alu domains. The Alu domain has a translational arrest function and consists of the SRP9 and SRP14 proteins bound to the terminal regions of the SRP RNA. So far, our understanding of the SRP and its evolution in lower eukaryotes such as protozoa and yeasts has been limited. However, genome sequences of such organisms have recently become available, and we have now analyzed this information with respect to genes encoding SRP components.

Results

A number of SRP RNA and SRP protein genes were identified by an analysis of genomes of protozoa and fungi. The sequences and secondary structures of the Alu portion of the RNA were found to be highly variable. Furthermore, proteins SRP9/14 appeared to be absent in certain species. Comparative analysis of the SRP RNAs from different Saccharomyces species resulted in models which contain features shared between all SRP RNAs, but also a new secondary structure element in SRP RNA helix 5. Protein SRP21, previously thought to be present only in Saccharomyces, was shown to be a constituent of additional fungal genomes. Furthermore, SRP21 was found to be related to metazoan and plant SRP9, suggesting that the two proteins are functionally related.

Conclusions

Analysis of a number of not previously annotated SRP components show that the SRP Alu domain is subject to a more rapid evolution than the other parts of the molecule. For instance, the RNA portion is highly variable and the protein SRP9 seems to have evolved into the SRP21 protein in fungi. In addition, we identified a secondary structure element in the Sacccharomyces RNA that has been inserted close to the Alu region. Together, these results provide important clues as to the structure, function and evolution of SRP.

MAPping the Eukaryotic Tree of Life: Structure, Function, and Evolution of the MAP215⧸Dis1 Family of Microtubule-Associated Proteins

The MAP215/Dis1 family of proteins is an evolutionarily ancient family of microtubule-associated proteins, with characterized members in all major kingdoms of eukaryotes, including fungi (Stu2 in S. cerevisiae, Dis1 and Alp14 in S. pombe), Dictyostelium (DdCP224), plants (Mor1 in A. thaliana and TMBP200 in N. tabaccum), and animals (Zyg9 in C. elegans, Msps in Drosophila, XMAP215 in Xenopus, and ch-TOG in humans). All MAP215/Dis1 proteins (with the exception of those in plants) localize to microtubule-organizing centers (MTOCs), including spindle pole bodies in yeast and centrosomes in animals, and all bind to microtubules in vitro and?or in vivo. Diverse roles in regulating microtubule assembly and organization have been proposed for individual family members, and a substantial body of evidence suggests that MAP215/Dis1-related proteins play critical roles in the assembly and function of the meiotic/mitotic spindles and/or cell division. An extensive search of public databases (including both EST and genome databases) identified partial sequences predicted to encode more than three dozen new members of the MAP215/Dis1 family, including putative MAP215/Dis1-related proteins in Giardia lamblia and four other protists, sixteen additional species of fungi, six plants, and twelve animals. The structure and function of MAP215/Dis1 proteins are discussed in relation to the evolution of this ancient family of microtubule-associated proteins.

The universally conserved HCF101 protein is involved in assembly of [4Fe-4S]-cluster-containing complexes in Arabidopsis thaliana chloroplasts

The seedling-lethal nuclear Arabidopsis hcf101 (high chlorophyll fluorescence) mutant is impaired in photosynthesis and complemented by the wild-type HCF101 cDNA. Photosystem I (PSI) activity is abolished, and PSI core complexes fail to accumulate in hcf101, whereas levels of other thylakoid membrane proteins are unaffected. Northern and in vivo labelling analyses as well as studies on polysome loading show that PSI transcript levels and translation rates of proteins, which belong to PSI, are normal in hcf101. PSI-specific fluorescence at 77 K is shifted from 735 to 728 nm in hcf101, indicating that exitons cannot efficiently be transferred to the PSI reaction centre, whereby the PSI antenna is almost unaffected. Mutant plants not only fail to accumulate mature PSI, which contains three [4Fe-4S]clusters (FSCs), but also are characterized by reduced levels of the soluble FSC-containing complex ferredoxin-thioredoxin reductase (FTR) in the stroma. Inhibited FTR maturation is not a secondary effect stemming from lack of PSI because the mutant hcf145, which also lacks PSI, accumulates FTR at normal levels. Levels of the [2Fe-2S] cluster-containing soluble and membrane proteins, ferredoxin and PetC, respectively, were unchanged in hcf101 plants. These data suggest a specific role of HCF101 in FSC biogenesis. HCF101 is plastid localized and belongs to an ancient and universally conserved family of P-loop ATPases previously designated as the 'MRP' (metGrelated protein) family. The function identified for HCF101 suggests a new designation, FSC, for this family.

Glutathione, Altruistic Metabolite in Fungi

Glutathione (GSH; gamma-L-glutamyl-L-cysteinyl-glycine), a non-protein thiol with a very low redox potential (E'0 = 240 mV for thiol-disulfide exchange), is present in high concentration up to 10 mM in yeasts and filamentous fungi. GSH is concerned with basic cellular functions as well as the maintenance of mitochondrial structure, membrane integrity, and in cell differentiation and development. GSH plays key roles in the response to several stress situations in fungi. For example, GSH is an important antioxidant molecule, which reacts non-enzymatically with a series of reactive oxygen species. In addition, the response to oxidative stress also involves GSH biosynthesis enzymes, NADPH-dependent GSH-regenerating reductase, glutathione S-transferase along with peroxide-eliminating glutathione peroxidase and glutaredoxins. Some components of the GSH-dependent antioxidative defence system confer resistance against heat shock and osmotic stress. Formation of protein-SSG mixed disulfides results in protection against desiccation-induced oxidative injuries in lichens. Intracellular GSH and GSH-derived phytochelatins hinder the progression of heavy metal-initiated cell injuries by chelating and sequestering the metal ions themselves and/or by eliminating reactive oxygen species. In fungi, GSH is mobilized to ensure cellular maintenance under sulfur or nitrogen starvation. Moreover, adaptation to carbon deprivation stress results in an increased tolerance to oxidative stress, which involves the induction of GSH-dependent elements of the antioxidant defence system. GSH-dependent detoxification processes concern the elimination of toxic endogenous metabolites, such as excess formaldehyde produced during the growth of the methylotrophic yeasts, by formaldehyde dehydrogenase and methylglyoxal, a by-product of glycolysis, by the glyoxalase pathway. Detoxification of xenobiotics, such as halogenated aromatic and alkylating agents, relies on glutathione S-transferases. In yeast, these enzymes may participate in the elimination of toxic intermediates that accumulate in stationary phase and/or act in a similar fashion as heat shock proteins. GSH S-conjugates may also form in a glutathione S-transferases-independent way, e.g. through chemical reaction between GSH and the antifugal agent Thiram. GSH-dependent detoxification of penicillin side-chain precursors was shown in Penicillium sp. GSH controls aging and autolysis in several fungal species, and possesses an anti-apoptotic feature.

MitoP2, an integrated database on mitochondrial proteins in yeast and man

The aim of the MitoP2 database (http://ihg.gsf.de/mitop2) is to provide a comprehensive list of mitochondrial proteins of yeast and man. Based on the current literature we created an annotated reference set of yeast and human proteins. In addition, data sets relevant to the study of the mitochondrial proteome are integrated and accessible via search tools and links. They include computational predictions of signalling sequences, and summarize results from proteome mapping, mutant screening, expression profiling, protein-protein interaction and cellular sublocalization studies. For each individual approach, specificity and sensitivity for allocating mitochondrial proteins was calculated. By providing the evidence for mitochondrial candidate proteins the MitoP2 database lends itself to the genetic characterization of human mitochondriopathies.

Gene loss, protein sequence divergence, gene dispensability, expression level, and interactivity are correlated in eukaryotic evolution

Lineage-specific gene loss, to a large extent, accounts for the differences in gene repertoires between genomes, particularly among eukaryotes. We derived a parsimonious scenario of gene losses for eukaryotic orthologous groups (KOGs) from seven complete eukaryotic genomes. The scenario involves substantial gene loss in fungi, nematodes, and insects. Based on this evolutionary scenario and estimates of the divergence times between major eukaryotic phyla, we introduce a numerical measure, the propensity for gene loss (PGL). We explore the connection among the propensity of a gene to be lost in evolution (PGL value), protein sequence divergence, the effect of gene knockout on fitness, the number of protein-protein interactions, and expression level for the genes in KOGs. Significant correlations between PGL and each of these variables were detected. Genes that have a lower propensity to be lost in eukaryotic evolution accumulate fewer substitutions in their protein sequences and tend to be essential for the organism viability, tend to be highly expressed, and have many interaction partners. The dependence between PGL and gene dispensability and interactivity is much stronger than that for sequence evolution rate. Thus, propensity of a gene to be lost during evolution seems to be a direct reflection of its biological importance.

Plasmodium Biology

The highly A+T rich genomes of human and rodent malarial parasites offer unprecedented glimpses of a lineage that is distinct from other model organisms. Plasmodium is distinguished by the presence of numerous low complexity inserts within globular domains of proteins. It displays several peculiarities in its transcription apparatus, and its DNA repair system appears to favor a certain innate level of mutability. Plasmodium possesses many cell surface molecules with "animal-like" adhesion modules. Potential genetic footprints of the ancestral eukaryotic algal precursor of the apicoplast are also detectable in its genome.

Rearrangements and chromosomal evolution

Comparisons of the genome sequences of related species suggests varying patterns of chromosomal rearrangements in different evolutionary lineages. In this review, I focus on the quantitative characterization of rearrangement processes and discuss specific inventories that have been compiled to date. Of particular interest are the statistical distribution of the lengths of inverted or locally transposed chromosome fragments (notably very short ones), inhomogeneities in susceptibility to evolutionary breakpoints in chromosomal regions, the relative importance of genome doubling in the history of multicellular eukaryotes, and of lateral transfer versus gene gain and loss in prokaryotes. These developments provide challenges to computational biologists to refine, revise and scale up mathematical models and algorithms for analyzing genome rearrangements.

The trans-membrane protein p25 forms highly specialized domains that regulate membrane composition and dynamics

Trans-membrane proteins of the p24 family are abundant, oligomeric proteins predominantly found in cis-Golgi membranes. They are not easily studied in vivo and their functions are controversial. We found that p25 can be targeted to the plasma membrane after inactivation of its canonical KKXX motif (KK to SS, p25SS), and that p25SS causes the co-transport of other p24 proteins beyond the Golgi complex, indicating that wild-type p25 plays a crucial role in retaining p24 proteins in cis-Golgi membranes. We then made use of these observations to study the intrinsic properties of these proteins, when present in a different membrane context. At the cell surface, the p25SS mutant segregates away from both the transferrin receptor and markers of lipid rafts, which are enriched in cholesterol and glycosphingolipids. This suggests that p25SS localizes to, or contributes to form, specialized membrane domains, presumably corresponding to oligomers of p25SS and other p24 proteins. Once at the cell surface, p25SS is endocytosed, together with other p24 proteins, and eventually accumulates in late endosomes, where it remains confined to well-defined membrane regions visible by electron microscopy. We find that this p25SS accumulation causes a concomitant accumulation of cholesterol in late endosomes, and an inhibition of their motility – two processes that are functionally linked. Yet, the p25SS-rich regions themselves seem to exclude not only Lamp1 but also accumulated cholesterol. One may envision that p25SS accumulation, by excluding cholesterol from oligomers, eventually overloads neighboring late endosomal membranes with cholesterol beyond their capacity (see Discussion). In any case, our data show that p25 and presumably other p24 proteins are endowed with the intrinsic capacity to form highly specialized domains that control membrane composition and dynamics. We propose that p25 and other p24 proteins control the fidelity of membrane transport by maintaining cholesterol-poor membranes in the Golgi complex.

Cryptosporidium parvum Cpn60 targets a relict organelle

Chaperonin 60 (Cpn60) is a well-established marker protein for eukaryotic mitochondria and plastids. In order to determine whether the small double-membrane-bounded organelle posterior to the nucleus in the apicomplexan Cryptosporidium parvum is a mitochondrion, the Cpn60 gene of C. parvum sporozoites ( CpCpn60) was analyzed and antibodies were generated for localization of the peptide. Sequence and phylogenetic analyses indicated that CpCpn60 is a mitochondrial isotype and that antibodies against it localize to the rough endoplasmic reticulum-enveloped remnant organelle of C. parvum sporozoites. These data show this organelle is of mitochondrial origin.

Comparative genomics, minimal gene-sets and the last universal common ancestor

Comparative genomics, using computational and experimental methods, enables the identification of a minimal set of genes that is necessary and sufficient for sustaining a functional cell. For most essential cellular functions, two or more unrelated or distantly related proteins have evolved; only about 60 proteins, primarily those involved in translation, are common to all cellular life. The reconstruction of ancestral life-forms is based on the principle of evolutionary parsimony, but the size and composition of the reconstructed ancestral gene-repertoires depend on relative rates of gene loss and horizontal gene-transfer. The present estimate suggests a simple last universal common ancestor with only 500-600 genes.

Mitochondrial remnant organelles of Giardia function in iron-sulphur protein maturation

Giardia intestinalis (syn. lamblia) is one of the most widespread intestinal protozoan pathogens worldwide, causing hundreds of thousands of cases of diarrhoea each year. Giardia is a member of the diplomonads, often described as an ancient protist group whose primitive nature is suggested by the lack of typical eukaryotic organelles (for example, mitochondria, peroxisomes), the presence of a poorly developed endomembrane system and by their early branching in a number of gene phylogenies. The discovery of nuclear genes of putative mitochondrial ancestry in Giardia and the recent identification of mitochondrial remnant organelles in amitochondrial protists such as Entamoeba histolytica and Trachipleistophora hominis suggest that the eukaryotic amitochondrial state is not a primitive condition but is rather the result of reductive evolution. Using an in vitro protein reconstitution assay and specific antibodies against IscS and IscU--two mitochondrial marker proteins involved in iron-sulphur cluster biosynthesis--here we demonstrate that Giardia contains mitochondrial remnant organelles (mitosomes) bounded by double membranes that function in iron-sulphur protein maturation. Our results indicate that Giardia is not primitively amitochondrial and that it has retained a functional organelle derived from the original mitochondrial endosymbiont.

Molecular phylogeny of the kelch-repeat superfamily reveals an expansion of BTB/kelch proteins in animals

Background

The kelch motif is an ancient and evolutionarily-widespread sequence motif of 44–56 amino acids in length. It occurs as five to seven repeats that form a β-propeller tertiary structure. Over 28 kelch-repeat proteins have been sequenced and functionally characterised from diverse organisms spanning from viruses, plants and fungi to mammals and it is evident from expressed sequence tag, domain and genome databases that many additional hypothetical proteins contain kelch-repeats. In general, kelch-repeat β-propellers are involved in protein-protein interactions, however the modest sequence identity between kelch motifs, the diversity of domain architectures, and the partial information on this protein family in any single species, all present difficulties to developing a coherent view of the kelch-repeat domain and the kelch-repeat protein superfamily. To understand the complexity of this superfamily of proteins, we have analysed by bioinformatics the complement of kelch-repeat proteins encoded in the human genome and have made comparisons to the kelch-repeat proteins encoded in other sequenced genomes.

Results

We identified 71 kelch-repeat proteins encoded in the human genome, whereas 5 or 8 members were identified in yeasts and around 18 in C. elegans, D. melanogaster and A. gambiae. Multiple domain architectures were identified in each organism, including previously unrecognised forms. The vast majority of kelch-repeat domains are predicted to form six-bladed β-propellers. The most prevalent domain architecture in the metazoan animal genomes studied was the BTB/kelch domain organisation and we uncovered 3 subgroups of human BTB/kelch proteins. Sequence analysis of the kelch-repeat domains of the most robustly-related subgroups identified differences in β-propeller organisation that could provide direction for experimental study of protein-binding characteristics.

Conclusion

The kelch-repeat superfamily constitutes a distinct and evolutionarily-widespread family of β-propeller domain-containing proteins. Expansion of the family during the evolution of multicellular animals is mainly accounted for by a major expansion of the BTB/kelch domain architecture. BTB/kelch proteins constitute 72 % of the kelch-repeat superfamily of H. sapiens and form three subgroups, one of which appears the most-conserved during evolution. Distinctions in propeller blade organisation between subgroups 1 and 2 were identified that could provide new direction for biochemical and functional studies of novel kelch-repeat proteins.

The COG database: an updated version includes eukaryotes

BACKGROUND

The availability of multiple, essentially complete genome sequences of prokaryotes and eukaryotes spurred both the demand and the opportunity for the construction of an evolutionary classification of genes from these genomes. Such a classification system based on orthologous relationships between genes appears to be a natural framework for comparative genomics and should facilitate both functional annotation of genomes and large-scale evolutionary studies.

RESULTS

We describe here a major update of the previously developed system for delineation of Clusters of Orthologous Groups of proteins (COGs) from the sequenced genomes of prokaryotes and unicellular eukaryotes and the construction of clusters of predicted orthologs for 7 eukaryotic genomes, which we named KOGs after eukaryotic orthologous groups. The COG collection currently consists of 138,458 proteins, which form 4873 COGs and comprise 75% of the 185,505 (predicted) proteins encoded in 66 genomes of unicellular organisms. The eukaryotic orthologous groups (KOGs) include proteins from 7 eukaryotic genomes: three animals (the nematode Caenorhabditis elegans, the fruit fly Drosophila melanogaster and Homo sapiens), one plant, Arabidopsis thaliana, two fungi (Saccharomyces cerevisiae and Schizosaccharomyces pombe), and the intracellular microsporidian parasite Encephalitozoon cuniculi. The current KOG set consists of 4852 clusters of orthologs, which include 59,838 proteins, or approximately 54% of the analyzed eukaryotic 110,655 gene products. Compared to the coverage of the prokaryotic genomes with COGs, a considerably smaller fraction of eukaryotic genes could be included into the KOGs; addition of new eukaryotic genomes is expected to result in substantial increase in the coverage of eukaryotic genomes with KOGs. Examination of the phyletic patterns of KOGs reveals a conserved core represented in all analyzed species and consisting of approximately 20% of the KOG set. This conserved portion of the KOG set is much greater than the ubiquitous portion of the COG set (approximately 1% of the COGs). In part, this difference is probably due to the small number of included eukaryotic genomes, but it could also reflect the relative compactness of eukaryotes as a clade and the greater evolutionary stability of eukaryotic genomes.

CONCLUSION

The updated collection of orthologous protein sets for prokaryotes and eukaryotes is expected to be a useful platform for functional annotation of newly sequenced genomes, including those of complex eukaryotes, and genome-wide evolutionary studies.

Intraspecies genotype variability of the microsporidian parasite Encephalitozoon hellem

Seven isolates of Encephalitozoon hellem from human immunodeficiency virus-positive patients were genotyped through a series of markers: the internal transcribed spacer (ITS) of ribosomal DNA, the polar tube protein (PTP) gene, and two intergenic spacers (IGS-TH and IGS-HZ) whose polymorphism is newly reported. The genome markers were all analyzed at three levels: PCR amplification followed by polyacrylamide gel electrophoresis, single-strand conformation analysis (SSCA), and DNA sequencing. The polymorphisms detected involve insertions/deletions and point mutations. SSCA can distinguish any pair of sequences, even those differing by a single base pair. The different isolates studied fit into the previously described ITS genotype 1A, except one which seems to be a 2A derivative variant (2D). When PTP and the new markers IGS-TH and IGS-HZ were analyzed, most of the isolates displayed different genotypes, demonstrating that E. hellem has a strong intraspecies variability. A set of markers such as those used here may be very useful in genotyping of clinical samples and in the assessment of epidemiological relationships among E. hellem strains.

Glycogen Phosphorylase Sequences from the Amitochondriate Protists, Trichomonas vaginalis, Mastigamoeba balamuthi, Entamoeba histolytica and Giardia intestinalis1

Glycogen phosphorylase genes or messages from four amitochondriate eukaryotes, Trichomonas vaginalis, Mastigamoeba balamuthi, Entamoeba histolytica (two genes) and Giardia intestinalis, have been isolated and sequenced. The sequences of the amitochondriate protist enzymes appear to share a most recent common ancestor. The clade containing these sequences is closest to that of another protist, the slime mold (Dictyostelium discoideum), and is more closely related to fungal and plant phosphorylases than to mammalian and eubacterial homologs. Structure-based amino acid alignment shows conservation of the residues and domains involved in catalysis and allosteric regulation by glucose 6-phosphate but high divergence at domains involved in phosphorylation-dependent regulation and AMP binding in fungi and animals. Protist phosphorylases, as their prokaryotic and plant counterparts, are probably not regulated by phosphorylation.

Integrated mapping, chromosomal sequencing and sequence analysis of Cryptosporidium parvum

The apicomplexan Cryptosporidium parvum is one of the most prevalent protozoan parasites of humans. We report the physical mapping of the genome of the Iowa isolate, sequencing and analysis of chromosome 6, and approximately 0.9 Mbp of sequence sampled from the remainder of the genome. To construct a robust physical map, we devised a novel and general strategy, enabling accurate placement of clones regardless of clone artefacts. Analysis reveals a compact genome, unusually rich in membrane proteins. As in Plasmodium falciparum, the mean size of the predicted proteins is larger than that in other sequenced eukaryotes. We find several predicted proteins of interest as potential therapeutic targets, including one exhibiting similarity to the chloroquine resistance protein of Plasmodium. Coding sequence analysis argues against the conventional phylogenetic position of Cryptosporidium and supports an earlier suggestion that this genus arose from an early branching within the Apicomplexa. In agreement with this, we find no significant synteny and surprisingly little protein similarity with Plasmodium. Finally, we find two unusual and abundant repeats throughout the genome. Among sequenced genomes, one motif is abundant only in C. parvum, whereas the other is shared with (but has previously gone unnoticed in) all known genomes of the Coccidia and Haemosporida. These motifs appear to be unique in their structure, distribution and sequences.

Homology modeling and calculation of the cobalt cluster charges of the Encephazlitozoon cuniculi methionine aminopeptidase, a potential target for drug design

Fumagillin is a potent anti-angiogenic drug used in cancer treatments. It is also one of the few molecules active against the Enterocytozoon and Encephalitozoon parasites responsible for various clinical syndromes in HIV-infected or immunosuppressive treated patients. Its toxicity, however, makes desirable the design of more specific molecules. The fumagillin target, as anti-angiogenic agent, is the methionine aminopeptidase, an ubiquitous metallo-enzyme responsible for the removing of the N-terminal methionine in nascent proteins. By analogy, it has been proposed that this enzyme could also be the target in the parasites. As a first approach to verify this and to determine if it would be possible to design a more specific derivative, we have built a homology model of the E. cuniculi aminopeptidase. The charges of the two cobalt ions present in the active site and of the side-chains involved in their binding were computed using ab-initio methods. A preliminary comparison of the interactions of the fumagillin and of a related compound, the TNP-470, with both the human and the parasitic enzymes strongly support the hypothesis that the parasitic aminopeptidase is indeed the target of the fumagillin. It also suggests that the TNP-470 interact identically with both enzymes while there could be small differences in case of the fumagillin.

Assessing functional divergence in EF-1alpha and its paralogs in eukaryotes and archaebacteria

A number of methods have recently been published that use phylogenetic information extracted from large multiple sequence alignments to detect sites that have changed properties in related protein families. In this study we use such methods to assess functional divergence between eukaryotic EF-1alpha (eEF-1alpha), archaebacterial EF-1alpha (aEF-1alpha) and two eukaryote-specific EF-1alpha paralogs-eukaryotic release factor 3 (eRF3) and Hsp70 subfamily B suppressor 1 (HBS1). Overall, the evolutionary modes of aEF-1alpha, HBS1 and eRF3 appear to significantly differ from that of eEF-1alpha. However, functionally divergent (FD) sites detected between aEF-1alpha and eEF-1alpha only weakly overlap with sites implicated as putative EF-1beta or aminoacyl-tRNA (aa-tRNA) binding residues in EF-1alpha, as expected based on the shared ancestral primary translational functions of these two orthologs. In contrast, FD sites detected between eEF-1alpha and its paralogs significantly overlap with the putative EF-1beta and/or aa-tRNA binding sites in EF-1alpha. In eRF3 and HBS1, these sites appear to be released from functional constraints, indicating that they bind neither eEF-1beta nor aa-tRNA. These results are consistent with experimental observations that eRF3 does not bind to aa-tRNA, but do not support the 'EF-1alpha-like' function recently proposed for HBS1. We re-assess the available genetic data for HBS1 in light of our analyses, and propose that this protein may function in stop codon-independent peptide release.

Incorporation of iron into Tritrichomonas foetus cell compartments reveals ferredoxin as a major iron-binding protein in hydrogenosomes

The intracellular transport of iron and its incorporation into organelles are poorly understood processes in eukaryotes and virtually unknown in parasitic protists. The transport of iron is of particular interest in trichomonads, which possess hydrogenosomes instead of mitochondria. The metabolic functions of hydrogenosomes, which contain a specific set of FeS proteins, entirely depend on iron acquisition. In this work the incorporation of iron into the cattle parasite Tritrichomonas foetus was monitored. Iron was efficiently taken up from (59)Fe-nitrilotriacetic acid and accumulated in the cytosol (88.9 %) and hydrogenosomes (4.7 % of the total radioactivity). Using atomic absorption spectrophotometry, an unusually high steady-state iron concentration in hydrogenosomes was determined [54.4+/-1.1 nmol Fe (mg protein)(-1)]. The concentration of iron in the cytosol was 13.4+/-0.5 nmol Fe (mg protein)(-1). Qualitative analysis of incorporated iron was performed using native gradient PAGE. The majority of the (59)Fe in the cytosol appeared as the labile-iron pool, which represents weakly bound iron associated with compounds of molecular mass ranging from 5000 to 30000 Da. Ferritin was not observed in Tt. foetus, nor in two other anaerobic protists, Entamoeba histolytica and Giardia intestinalis. Analysis of Tt. foetus hydrogenosomes showed at least nine iron-binding compounds, which were absent in metronidazole-resistant mutants. The major iron-binding compound was identified as [2Fe-2S] ferredoxin of the adrenodoxin type.

400000 nematode ESTs on the Net

The parasitic nematode expressed sequence tag (EST) project, a collaboration between University of Edinburgh and the Wellcome Trust Sanger Institute in the UK and the Genome Sequencing Center, St Louis, MO, USA, is currently generating sequence information from >30 different species of nematode. Over 400000 nematode ESTs are now available and at least another 130000 are planned. Here, an update is provided on the status of the project and describes the database tools being developed to disseminate these data.

Phylogenomics: intersection of evolution and genomics

Much has been gained from genomic and evolutionary studies of species. Combining the perspectives of these different approaches suggests that an integrated phylogenomic approach will be beneficial.

Molecular characterization of Encephalitozoon intestinalis (Microspora) replication kinetics in a murine intestinal cell line

Microsporidia are obligate intracellular pathogens of invertebrate and vertebrate animals. Most human infections are caused by Enterocytozoon bieneusi or Encephalitozoon intestinalis, and result in chronic diarrhea. In order to determine the signals involved in microsporidial spore activation and invasion, kinetics of in vitro E. intestinalis replication were defined using real-time quantitative PCR. Segments of small subunit ribosomal RNA and polar tube protein 2 genes of E. intestinalis were used to quantify parasite gene copy number following infection in murine colon carcinoma cells. Parasite DNA was detectable in small but significant amounts within host cells as early as 4 h postinoculation, genome replication was completed by 36 h, and parasite progeny were released into the supernatant beginning 72 h postinoculation. Heat-treating spores did not prevent transfer of parasite DNA into cells, but did inhibit parasite replication. Treating cell cultures with albendazole suppressed but did not completely inhibit parasite replication. These results confirm observations that E. intestinalis completes its life cycle within the turnover time of its target host cells; invasion into susceptible host cells occurs independently of spore viability; and real-time quantitative PCR is a sensitive and reproducible method with which to monitor microsporidial infection under varying treatments or conditions.

Analysis of four human microsporidian isolates by MALDI-TOF mass spectrometry

Spores of four species of microsporidia isolated from humans were analyzed by matrix-assisted laser desorption ionization time of flight mass spectrometry (MALDI-TOF MS) and specific biomarkers were found for each. The microsporidia analyzed included three species, Encephalitozoon cuniculi, Encephalitozoon hellem, and Encephalitozoon intestinalis and the fourth organism is the recently described Brachiola algerae. Whole spores, spore shells, and soluble fractions were applied directly to the MALDI target without further purification steps. MALDI-TOF MS analysis of both whole spores and soluble fractions of the four isolates revealed a group of unique, characteristic, and reproducible spectral markers in the mass range of 2,000-8,000 Da. Statistical analysis of the averaged centroided masses uncovered two distinct sets of unique peptides or biomarkers, one originated from whole spores and the other from soluble fractions, that can differentiate the four microsporidian species studied. MALDI-TOF MS analysis of whole organisms is a rapid, sensitive, and specific option to characterize microsporidian isolates and has the potential for several applications in parasitology.

Phylogeny and taxonomic revision of plastid-containing euglenophytes based on SSU rDNA sequence comparisons and synapomorphic signatures in the SSU rRNA secondary structure

Sequence comparisons and a revised classification of the Euglenophyceae were based on 92 new SSU rDNA sequences obtained from strains of Euglena, Astasia, Phacus, Trachelomonas, Colacium, Cryptoglena, Lepocinclis, Eutreptia, Eutreptiella and Tetreutreptia. Sequence data also provided molecular signatures for taxa from genus to class level in the SSU rRNA secondary structure, revealed by a novel approach (search for non-homoplasious synapomorphies) and used for taxonomic diagnoses. Photosynthetic euglenoids and secondary heterotrophs formed a clade, designated as Euglenophyceae (emend.) with two orders: Euglenales and Eutreptiales. The mostly marine Eutreptiales (Eutreptia, Eutreptiella; not Distigma) comprised taxa with two or four emergent flagella (the quadriflagellate Tetreutreptia was integrated within Eutreptiella). The Euglenales (freshwater genera with < or = one emergent flagellum) formed nine clades and two individual branches (single strains); however, only two clades were congruent with traditional genera: Trachelomonas (incl. Strombomonas) and Colacium. Euglena was polyphyletic and diverged into four independent clades (intermixed with Astasia, Khawkinea and Lepocinclis) and two individual branches (e.g. E. polymorpha). Phacus was also subdivided into Phacus s. str. and two combined lineages (mixed with Lepocinclis spp. or Cryptoglena). In consequence, Euglena (s. str.), Phacus and other genera were emended and one lineage (mixed Phacus/Lepocinclis-clade) was recognized as the previously neglected genus Monomorphina Mereschkowsky (1877). The sister clade of Phacus s. str. (mixed Euglena/Lepocinclis-clade) was identified as Lepocinclis Perty (emended).

Congruent evidence from alpha-tubulin and beta-tubulin gene phylogenies for a zygomycete origin of microsporidia

The origin of microsporidia and the evolutionary relationships among the major lineages of fungi have been examined by molecular phylogeny using alpha-tubulin and beta-tubulin. Chytrids, basidiomycetes, ascomycetes, and microsporidia were all recovered with high support, and the zygomycetes were consistently paraphyletic. The microsporidia were found to branch within zygomycetes, and showed relationships with members of the Entomophthorales and Zoopagales. This provides support for the microsporidia having evolved from within the fungi, however, the tubulin genes are difficult to interpret unambiguously since fungal and microsporidian tubulins are very divergent. Rapid evolutionary rates a characteristic of practically all microsporidian genes studied, so determining their evolutionary history will never be completely free of such difficulties. While the tubulin phylogenies do not provide a decisive conclusion, they do further narrow the probable origin of microsporidia to a zygomycete-like ancestor.

Mitochondrial connection to the origin of the eukaryotic cell

Phylogenetic evidence is presented that primitively amitochondriate eukaryotes containing the nucleus, cytoskeleton, and endomembrane system may have never existed. Instead, the primary host for the mitochondrial progenitor may have been a chimeric prokaryote, created by fusion between an archaebacterium and a eubacterium, in which eubacterial energy metabolism (glycolysis and fermentation) was retained. A Rickettsia-like intracellular symbiont, suggested to be the last common ancestor of the family Rickettsiaceae and mitochondria, may have penetrated such a host (pro-eukaryote), surrounded by a single membrane, due to tightly membrane-associated phospholipase activity, as do present-day rickettsiae. The relatively rapid evolutionary conversion of the invader into an organelle may have occurred in a safe milieu via numerous, often dramatic, changes involving both partners, which resulted in successful coupling of the host glycolysis and the symbiont respiration. Establishment of a potent energy-generating organelle made it possible, through rapid dramatic changes, to develop genuine eukaryotic elements. Such sequential, or converging, global events could fill the gap between prokaryotes and eukaryotes known as major evolutionary discontinuity.

Characterization of aldehyde ferredoxin oxidoreductase gene defective mutant in Magnetospirillum magneticum AMB-1

A non-magnetic mutant of Magnetospirillum magneticum AMB-1, designated as NMA21, was generated by mini-Tn5 transposon mutagenesis to identify genes involved in bacterial magnetic particle (BMP) synthesis. Alignment of the DNA sequences flanking the transposon allowed the isolation of an open reading frame (ORF2) within an operon consisting of five genes. The amino acid sequence of ORF2 showed homology with tungsten-containing aldehyde ferredoxin oxidoreductase (AOR) from Pyrococcus furiosus (48% identity and 64% similarity), which functions for aldehyde oxidation. AOR was found to be expressed under microaerobic conditions and localized in the cytoplasm of AMB-1. Iron uptake and growth of NMA21 were lower than wild type. Transmission electron microscopy (TEM) of NMA21 revealed that no BMPs were completely synthesized, but polyhydroxybutyrate (PHB)-like granules were persistently produced. These results indicate that AOR may contribute to ferric iron reduction during BMP synthesis in M. magneticum AMB-1 under microaerobic respiration.

Phylogenetic relationships of non-mitochondrial nucleotide transport proteins in bacteria and eukaryotes

Current knowledge about the nucleotide metabolism of intracellular bacteria is very limited. Here we report on the identification of nucleotide transport proteins (NTT) of two obligate endoparasites, Caedibacter caryophila and Holospora obtusa, both alpha-proteobacteria, which reside in the vegetative macronucleus of Paramecium caudatum. For comparative studies, we also identified the first nucleotide transporter in chloroplasts of a red alga, i.e. Galdieria sulphuraria, and further homologs in plant chloroplasts. Heterologous expression of the NTT proteins from C. caryophila, H. obtusa, and G. sulphuraria in Escherichia coli demonstrate that the nucleotide influx mediated by these transporters is specific for ATP and ADP. The NTT proteins of C. caryophila and H. obtusa exhibit substantial sequence identity with their counterparts in chloroplasts and intracellular bacterial pathogens of humans, but not with the nucleotide transport system of mitochondria. Comprehensive phylogenetic analyses of bacterial and chloroplast NTT proteins showed that homologs in chloroplasts from plants, and green, red, stramenopile and glaucocystophyte algae are monophyletic. In contrast, the evolutionary relationships of the bacterial counterparts appear highly complex. In the presented phylogeny, NTT proteins of C. caryophila and H. obtusa are only distantly related to one another, although these two taxa are close relatives in 16S rRNA trees. The tree topology indicates that some bacterial NTT paralogs have arisen by gene duplications and others by horizontal transfer.

Molecular evolution meets the genomics revolution

Changes in technology in the past decade have had such an impact on the way that molecular evolution research is done that it is difficult now to imagine working in a world without genomics or the Internet. In 1992, GenBank was less than a hundredth of its current size and was updated every three months on a huge spool of tape. Homology searches took 30 minutes and rarely found a hit. Now it is difficult to find sequences with only a few homologs to use as examples for teaching bioinformatics. For molecular evolution researchers, the genomics revolution has showered us with raw data and the information revolution has given us the wherewithal to analyze it. In broad terms, the most significant outcome from these changes has been our newfound ability to examine the evolution of genomes as a whole, enabling us to infer genome-wide evolutionary patterns and to identify subsets of genes whose evolution has been in some way atypical.

Conserved N-terminal motifs of telomerase reverse transcriptase required for ribonucleoprotein assembly in vivo

Telomerase is a ribonucleoprotein (RNP) reverse transcriptase responsible for the maintenance of one strand of the telomere terminal repeats. The key protein subunit of the telomerase complex, known as TERT, possesses reverse transcriptase (RT)-like motifs that directly mediate nucleotide addition. The RT motifs are located in the C-terminal region of the polypeptide. Sequence alignments also revealed the existence of four conserved motifs (named GQ, CP, QFP, and T) in the N-terminal region of TERT. The GQ motif of yeast TERT has been demonstrated previously to be essential for telomerase catalysis and may participate in RNP formation. In this report, we show that substitution of conserved residues in the CP, QFP, and T motifs of yeast TERT also impairs both telomere maintenance and telomerase activity, thus confirming the validity of the sequence alignment. The extent of telomere shortening correlates with the extent of reduction in the level of telomerase activity, TERT protein, and TERT-associated TLC1 RNA. Overexpression of the mutant proteins does not result in telomere shortening, implying that assembly rather than catalytic function was affected. This notion was further supported by comparing the efficiency of RNP formation in the wild type and the overexpression strains. Taken together, our results show that three of the four N-terminal motifs are required for efficient telomerase RNP formation in vivo but not for the enzymatic function of telomerase. We also show that the majority of telomerase-associated TLC1 RNA has a more upstream 3' end than previously reported, consistent with additional processing events during RNP maturation.

Fish microsporidia: fine structural diversity and phylogeny

Structural diversity of fish microsporidian life cycle stages and of the host-parasite interface is reviewed. In the infected cell of the fish host, microsporidia may either cause serious degradation of the cytoplasm and demise of the cell, or they may elicit host cell hypertrophy, producing a parasite-hypertrophic host cell complex, the xenoma. The structure of the xenoma and of its cell wall may differ according to the genus of the parasite, and seems to express properties of the parasite rather than those of the host. In merogony, the parasite cell surface interacts with the host cell in diverse ways, the most conspicuous being the production of thick envelopes of different types. Sporogony stages reveal different types of walls or membranes encasing the sporoblasts and later the spores and these envelopes may be of host or parasite origin. Nucleospora differs from all other fish microsporidia by its unique process of sporogony. Except for the formation of conspicuous xenomas, there are no essentially different structures in fish-infecting microsporidia compared with microsporidia from other hosts. Although the structures associated with the development of fish microsporidia cannot be attributed importance in tracing the phylogeny, they are relevant for practical determination and assessing the relation to the host. The possibility of the existence of an intermediate host is discussed. Higher-level classification of Microsporidia is briefly discussed and structure and evolutionary rates in microsporidian rDNA are reviewed. Discussion of rDNA molecular phylogeny of fish-infecting microsporidia is followed by classification of these parasites. Most form a rather cohesive clade. Outside this clade is the genus Nucleospora, separated at least at the level of Order. Within the main clade, however, there are six species infecting hosts other than fish. Based on data available for analysis, a tentative classification of fish-infecting microsporidia into five groups is proposed. Morphologically defined groups represent families, others are referred to as clades. Group 1, represented by family Pleistophoridae, includes Pleistophora, Ovipleistophora and Heterosporis; Vavraia and Trachipleistophora infect non-fish hosts. Group 2, represented by family Glugeidae, is restricted to genus Glugea and Tuzetia weidneri from crustaceans. Group 3 comprises three clades: Loma and a hyperparasitic microsporidian from a myxosporean; Ichthyosporidium and Pseudoloma clade and the Loma acerinae clade. For the latter species a new genus has to be established. Group 4 contains two families, Spragueidae with the genus Spraguea and Tetramicridae with genera Microgemma and Tetramicra, and the Kabatana and Microsporidium seriolae clade. Group 5 is represented by the family Enterocytozoonidae with the genus Nucleospora and mammal-infecting genus Enterocytozoon.

The homeobox genes of Encephalitozoon cuniculi (Microsporidia) reveal a putative mating-type locus

Homeobox genes have been found in animals, fungi, and plants. Recently, the complete genomic sequence of Encephalitozoon cuniculi has become available and it was shown that this Microsporida species is related to fungi. Given this close relatedness to fungi, I have searched the genome of E. cuniculifor homeobox genes. There are 12 homeobox genes as well as one STE12 orthologue. The large number of homeobox genes when compared to the annotations, which only list one, suggests that possibly other smaller genes may have been overlooked in the published analysis and E. cuniculi may have more than 1,997 genes. Sequence and phylogenetic analyses show that the E. cuniculi homeobox genes also fall into two distinct groups, with two TALE and ten typical homeobox genes being present. Like in the mating-type locus of yeast and other fungi, one TALE and one typical homeobox are found in close proximity of each other on the chromosome, suggesting that Microsporidia also contain a mating-type locus.

Mobile genetic elements colonizing the genomes of metazoan parasites

A substantial fraction of the genome of most eukaryotes, including those of metazoan parasites, is predicted to comprise repetitive sequences. Mobile genetic elements (MGEs) will make up much of these repetitive sequences, particularly the interspersed sequences. This article reviews information on MGEs that have colonized the genomes of metazoan parasites (i.e. parasites of parasites). Helminth and mosquito genomes, in particular, are compared with those of better-understood model organisms. MGEs from the genomes of metazoan parasites can be expected to have practical uses in transgenesis and epidemiological studies.

Mitochondria and hydrogenosomes are two forms of the same fundamental organelle

Published data suggest that hydrogenosomes, organelles found in diverse anaerobic eukaryotes that make energy and hydrogen, were once mitochondria. As hydrogenosomes generally lack a genome, the conversion is probably one way. The sources of the key hydrogenosomal enzymes, pyruvate : ferredoxin oxidoreductase (PFO) and hydrogenase, are not resolved by current phylogenetic analyses, but it is likely that both were present at an early stage of eukaryotic evolution. Once thought to be restricted to a few unusual anaerobic eukaryotes, the proteins are intimately integrated into the fabric of diverse eukaryotic cells, where they are targeted to different cell compartments, and not just hydrogenosomes. There is no evidence supporting the view that PFO and hydrogenase originated from the mitochondrial endosymbiont, as posited by the hydrogen hypothesis for eukaryogenesis. Other organelles derived from mitochondria have now been described in anaerobic and parasitic microbial eukaryotes, including species that were once thought to have diverged before the mitochondrial symbiosis. It thus seems possible that all eukaryotes may eventually be shown to contain an organelle of mitochondrial ancestry, to which different types of biochemistry can be targeted. It remains to be seen if, despite their obvious differences, this family of organelles shares a common function of importance for the eukaryotic cell, other than energy production, that might provide the underlying selection pressure for organelle retention.

On the origins of cells: a hypothesis for the evolutionary transitions from abiotic geochemistry to chemoautotrophic prokaryotes, and from prokaryotes to nucleated cells

All life is organized as cells. Physical compartmentation from the environment and self-organization of self-contained redox reactions are the most conserved attributes of living things, hence inorganic matter with such attributes would be life's most likely forebear. We propose that life evolved in structured iron monosulphide precipitates in a seepage site hydrothermal mound at a redox, pH and temperature gradient between sulphide-rich hydrothermal fluid and iron(II)-containing waters of the Hadean ocean floor. The naturally arising, three-dimensional compartmentation observed within fossilized seepage-site metal sulphide precipitates indicates that these inorganic compartments were the precursors of cell walls and membranes found in free-living prokaryotes. The known capability of FeS and NiS to catalyse the synthesis of the acetyl-methylsulphide from carbon monoxide and methylsulphide, constituents of hydrothermal fluid, indicates that pre-biotic syntheses occurred at the inner surfaces of these metal-sulphide-walled compartments, which furthermore restrained reacted products from diffusion into the ocean, providing sufficient concentrations of reactants to forge the transition from geochemistry to biochemistry. The chemistry of what is known as the RNA-world could have taken place within these naturally forming, catalyticwalled compartments to give rise to replicating systems. Sufficient concentrations of precursors to support replication would have been synthesized in situ geochemically and biogeochemically, with FeS (and NiS) centres playing the central catalytic role. The universal ancestor we infer was not a free-living cell, but rather was confined to the naturally chemiosmotic, FeS compartments within which the synthesis of its constituents occurred. The first free-living cells are suggested to have been eubacterial and archaebacterial chemoautotrophs that emerged more than 3.8 Gyr ago from their inorganic confines. We propose that the emergence of these prokaryotic lineages from inorganic confines occurred independently, facilitated by the independent origins of membrane-lipid biosynthesis: isoprenoid ether membranes in the archaebacterial and fatty acid ester membranes in the eubacterial lineage. The eukaryotes, all of which are ancestrally heterotrophs and possess eubacterial lipids, are suggested to have arisen ca. 2 Gyr ago through symbiosis involving an autotrophic archaebacterial host and a heterotrophic eubacterial symbiont, the common ancestor of mitochondria and hydrogenosomes. The attributes shared by all prokaryotes are viewed as inheritances from their confined universal ancestor. The attributes that distinguish eubacteria and archaebacteria, yet are uniform within the groups, are viewed as relics of their phase of differentiation after divergence from the non-free-living universal ancestor and before the origin of the free-living chemoautotrophic lifestyle. The attributes shared by eukaryotes with eubacteria and archaebacteria, respectively, are viewed as inheritances via symbiosis. The attributes unique to eukaryotes are viewed as inventions specific to their lineage. The origin of the eukaryotic endomembrane system and nuclear membrane are suggested to be the fortuitous result of the expression of genes for eubacterial membrane lipid synthesis by an archaebacterial genetic apparatus in a compartment that was not fully prepared to accommodate such compounds, resulting in vesicles of eubacterial lipids that accumulated in the cytosol around their site of synthesis. Under these premises, the most ancient divide in the living world is that between eubacteria and archaebacteria, yet the steepest evolutionary grade is that between prokaryotes and eukaryotes.

Genomic reduction and evolution of novel genetic membranes and protein-targeting machinery in eukaryote-eukaryote chimaeras (meta-algae)

Chloroplasts originated just once, from cyanobacteria enslaved by a biciliate protozoan to form the plant kingdom (green plants, red and glaucophyte algae), but subsequently, were laterally transferred to other lineages to form eukaryote-eukaryote chimaeras or meta-algae. This process of secondary symbiogenesis (permanent merger of two phylogenetically distinct eukaryote cells) has left remarkable traces of its evolutionary role in the more complex topology of the membranes surrounding all non-plant (meta-algal) chloroplasts. It took place twice, soon after green and red algae diverged over 550 Myr ago to form two independent major branches of the eukaryotic tree (chromalveolates and cabozoa), comprising both meta-algae and numerous secondarily non-photosynthetic lineages. In both cases, enslavement probably began by evolving a novel targeting of endomembrane vesicles to the perialgal vacuole to implant host porter proteins for extracting photosynthate. Chromalveolates arose by such enslavement of a unicellular red alga and evolution of chlorophyll c to form the kingdom Chromista and protozoan infrakingdom Alveolata, which diverged from the ancestral chromalveolate chimaera. Cabozoa arose when the common ancestor of euglenoids and cercozoan chlorarachnean algae enslaved a tetraphyte green alga with chlorophyll a and b. I suggest that in cabozoa the endomembrane vesicles originally budded from the Golgi, whereas in chromalveolates they budded from the endoplasmic reticulum (ER) independently of Golgi-targeted vesicles, presenting a potentially novel target for drugs against alveolate Sporozoa such as malaria parasites and Toxoplasma. These hypothetical ER-derived vesicles mediated fusion of the perialgal vacuole and rough ER (RER) in the ancestral chromist, placing the former red alga within the RER lumen. Subsequently, this chimaera diverged to form cryptomonads, which retained the red algal nucleus as a nucleomorph (NM) with approximately 464 protein-coding genes (30 encoding plastid proteins) and a red or blue phycobiliprotein antenna pigment, and the chromobiotes (heterokonts and haptophytes), which lost phycobilins and evolved the brown carotenoid fucoxanthin that colours brown seaweeds, diatoms and haptophytes. Chromobiotes transferred the 30 genes to the nucleus and lost the NM genome and nuclear-pore complexes, but retained its membrane as the periplastid reticulum (PPR), putatively the phospholipid factory of the periplastid space (former algal cytoplasm), as did the ancestral alveolate independently. The chlorarachnean NM has three minute chromosomes bearing approximately 300 genes riddled with pygmy introns. I propose that the periplastid membrane (PPM, the former algal plasma membrane) of chromalveolates, and possibly chlorarachneans, grows by fusion of vesicles emanating from the NM envelope or PPR. Dinoflagellates and euglenoids independently lost the PPM and PPR (after diverging from Sporozoa and chlorarachneans, respectively) and evolved triple chloroplast envelopes comprising the original plant double envelope and an extra outermost membrane, the EM, derived from the perialgal vacuole. In all metaalgae most chloroplast proteins are coded by nuclear genes and enter the chloroplast by using bipartite targeting sequences--an upstream signal sequence for entering the ER and a downstream chloroplast transit sequence. I present a new theory for the four-fold diversification of the chloroplast OM protein translocon following its insertion into the PPM to facilitate protein translocation across it (of both periplastid and plastid proteins). I discuss evidence from genome sequencing and other sources on the contrasting modes of protein targeting, cellular integration, and evolution of these two major lineages of eukaryote "cells within cells". They also provide powerful evidence for natural selection's effectiveness in eliminating most functionless DNA and therefore of a universally useful non-genic function for nuclear non-coding DNA, i.e. most DNA in the biosphere, and dramatic examples of genomic reduction. I briefly argue that chloroplast replacement in dinoflagellates, which happened at least twice, may have been evolutionarily easier than secondary symbiogenesis because parts of the chromalveolate protein-targeting machinery could have helped enslave the foreign plastids.

On the origin of mitochondria: a genomics perspective

The availability of complete genome sequence data from both bacteria and eukaryotes provides information about the contribution of bacterial genes to the origin and evolution of mitochondria. Phylogenetic analyses based on genes located in the mitochondrial genome indicate that these genes originated from within the alpha-proteobacteria. A number of ancestral bacterial genes have also been transferred from the mitochondrial to the nuclear genome, as evidenced by the presence of orthologous genes in the mitochondrial genome in some species and in the nuclear genome of other species. However, a multitude of mitochondrial proteins encoded in the nucleus display no homology to bacterial proteins, indicating that these originated within the eukaryotic cell subsequent to the acquisition of the endosymbiont. An analysis of the expression patterns of yeast nuclear genes coding for mitochondrial proteins has shown that genes predicted to be of eukaryotic origin are mainly translated on polysomes that are free in the cytosol whereas those of putative bacterial origin are translated on polysomes attached to the mitochondrion. The strong relationship with alpha-proteobacterial genes observed for some mitochondrial genes, combined with the lack of such a relationship for others, indicates that the modern mitochondrial proteome is the product of both reductive and expansive processes.

Microsporidia: biology and evolution of highly reduced intracellular parasites

Microsporidia are a large group of microbial eukaryotes composed exclusively of obligate intracellular parasites of other eukaryotes. Almost 150 years of microsporidian research has led to a basic understanding of many aspects of microsporidian biology, especially their unique and highly specialized mode of infection, where the parasite enters its host through a projectile tube that is expelled at high velocity. Molecular biology and genomic studies on microsporidia have also drawn attention to many other unusual features, including a unique core carbon metabolism and genomes in the size range of bacteria. These seemingly simple parasites were once thought to be the most primitive eukaryotes; however, we now know from molecular phylogeny that they are highly specialized fungi. The fungal nature of microsporidia indicates that microsporidia have undergone severe selective reduction permeating every level of their biology: From cell structures to metabolism, and from genomics to gene structure, microsporidia are reduced.