Mitochondrial Type Iron-Sulfur Cluster Assembly in the Amitochondriate Eukaryotes Trichomonas vaginalis and Giardia intestinalis, as Indicated by the Phylogeny of IscS

Jan Tachezy, Lidya B. Sánchez and Miklós Müller

Department of Parasitology, Faculty of Science, Charles University, Prague, Czech Republic;
Rockefeller University, New York, NY, USA


    Abstract
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
Pyridoxal-5'-phosphate-dependent cysteine desulfurase (IscS) is an essential enzyme in the assembly of FeS clusters in bacteria as well as in the mitochondria of eukaryotes. Although FeS proteins are particularly important for the energy metabolism of amitochondrial anaerobic eukaryotes, there is no information about FeS cluster formation in these organisms. We identified and sequenced two IscS homologs of Trichomonas vaginalis (TviscS-1 and TviscS-2) and one of Giardia intestinalis (GiiscS). TviscS-1, TviscS-2, and GiiscS possess the typical conserved regions implicated in cysteine desulfurase activity. N-termini of TviscS-1 and TviscS-2 possess eight amino acid extensions, which resemble the N-terminal presequences that target proteins to hydrogenosomes in trichomonads. No presequence was evident in GiiscS from Giardia, an organism that apparently lacks hydrogenosmes or mitochondria. Phylogenetic analysis showed a close relationship among all eukaryotic IscS genes including those of amitochondriates. IscS of proteobacteria formed a sister group to the eukaryotic clade, suggesting that isc-related genes were present in the proteobacterial endosymbiotic ancestor of mitochondria and hydrogenosomes. NifS genes of nitrogen-fixing bacteria, which are IscS homologs required for specific formation of FeS clusters in nitrogenase, formed a more distant group. The phylogeny indicates the presence of a common mechanism for FeS cluster formation in mitochondriates as well as in amitochondriate eukaryotes. Furthermore, the analyses support a common origin of Trichomonas hydrogenosomes and mitochondria, as well as secondary loss of mitochondrion/hydrogenosome-like organelles in Giardia.


    Introduction
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
Amitochondriate eukaryotes can be divided into two metabolic types (Martin and Müller 1998Citation ; Müller 1998Citation ). Type I organisms such as Giardia and Entamoeba lack organelles involved in core energy metabolism, while type II organisms (trichomonads, some ciliates, and chytrid fungi) harbor a double-membrane limited organelle, the hydrogenosome (Müller 1993Citation ; Hackstein et al. 1999Citation ; Kulda 1999Citation ). The hydrogenosome is the site of the FeS protein-mediated metabolism of pyruvate and the formation of molecular hydrogen, which is accompanied by substrate-level phosphorylation ATP synthesis. In type I amitochondriates, the FeS protein-dependent pyruvate metabolism takes place in the cytosol (Reeves 1984Citation ; Ellis et al. 1993Citation ). The evolutionary origin of amitochondriate eukaryotes is much debated. Although they were long regarded as ancestral premitochondrial lineages (Cavalier-Smith 1987Citation ), recent evidence suggests that they had experienced the endosymbiotic event leading to the establishment of the mitochondrion (Embley and Hirt 1998Citation ; Martin and Müller 1998Citation ; Roger 1999Citation ). In type II organisms, hydrogenosomes are regarded as descendants of a common endosymbiont that evolved to either mitochondria or hydrogenosomes (Martin and Müller 1998Citation ; Müller 1997Citation ). A common origin of the two organelles is supported by a number of similarities in their structure, function, and biogenesis (Johnson, Lahti, and Bradley 1993Citation ; Benchimol, Johnson, and De Souza 1996Citation ; Bui, Bradley, and Johnson 1996Citation ; Bradley et al. 1997Citation ; Dyall et al. 2000Citation ), as well as by phylogenetic analysis of several hydrogenosomal metabolic enzymes (Länge, Rozario, and Müller 1994Citation ; Hrdy and Müller 1995a, 1995bCitation ) and heat shock proteins (Müller 1997Citation ; Embley and Hirt 1998Citation ). Although neither mitochondria nor hydrogenosomes have been found in type I organisms, genes of probable mitochondrial origin have been identified in Giardia (Roger et al. 1998Citation ) and Entamoeba (Clark and Roger 1995Citation ). Moreover, a putative mitochondrial "remnant," the mitosome (Tovar, Fischer, and Clark 1999Citation ) or crypton (Mai et al. 1999Citation ), has recently been detected in Entamoeba (Müller 2000Citation ). These considerations led us to the hypothesis that FeS proteins operating in the energy metabolism of mitochondrial as well as secondarily amitochondrial organisms were present in the common ancestral organelle. If so, a common mechanism of FeS cluster assembly may operate in amitochondriate and in mitochondrial eukaryotes.

In spite of the importance of FeS proteins for all living cells, little is known of how and where FeS clusters are synthesized in vivo and which proteins are involved in their insertion into the apoproteins. The best-characterized enzyme participating in this process is a pyridoxal-5'-phosphate-dependent cysteine desulfurase which catalyzes the formation of L-alanine and elemental sulfur by using L-cysteine as substrate. Initially, the enzyme was described as NifS in Azotobacter vinelandii (Zheng et al. 1993Citation ), in which it provides sulfur for FeS cluster formation in nitrogenase (Zheng and Dean 1994Citation ). Later, a NifS homolog designated IscS (iron-sulfur cluster) was found in A. vinelandii, as well as in a number of non–nitrogen-fixing bacteria. It has been proposed that IscS plays a general role in the formation of FeS clusters or repair of FeS proteins with a housekeeping function (Zheng et al. 1998Citation ). More recently, a major role of IscS in de novo FeS cluster synthesis has been demonstrated using an iscS deletion strain of Escherichia coli (Schwartz et al. 2000Citation ). Importantly, IscS homologs have been identified in the genomes of diverse eukaryotes (Arabidopsis, Caenorhabditis, Drosophila, Homo, Mus, Saccharomyces), suggesting a general role for IscS in FeS cluster formation. In mice (Nakai et al. 1998Citation ) and yeast (Strain et al. 1998Citation ), IscS is localized in mitochondria. In human cells, the IscS homologs are targeted either to mitochondria or to the cytosol and nucleus (Land and Rouault 1998Citation ). Mutation in an iscS-like gene in yeast (NFS1) caused reduction in the activities of the mitochondrial FeS proteins, aconitase and succinate dehydrogenase (Strain et al. 1998Citation ). According to Kispal et al. (1999)Citation and Lill and Kispal (2000)Citation , mitochondria also play a crucial role in the FeS cluster formation of extramitochondrial FeS proteins. In addition, IscS homologs have been found to mediate several other functions that are independent of FeS cluster assembly but require IscS as a sulfur donor. Thus, IscS homologs are involved in biosynthesis of thiamin (Lauhon and Kambampati 2000Citation ), NAD (Sun and Setlow 1993Citation ), 4-thiouridine (Kambampati and Lauhon 1999Citation ), and molybdopterin (Amrani et al. 2000Citation ). Finally, IscS/NifS homologs mediate release of elemental selenium from L-selenocysteine (Mihara et al. 2000Citation ), and they may participate in tRNA splicing (Kolman and Söll 1993Citation ).

IscS has not been reported in amitochondriate eukaryotes, although FeS proteins are of particular importance for these organisms (Müller 1998Citation ). Here we report the identification of genes encoding IscS in the type II organism Trichomonas vaginalis and the type I organism Giardia intestinalis. Phylogenetic analysis indicates the presence of a common mechanism for FeS cluster formation in mitochondria and hydrogenosomes, as well as in organisms that secondarily lost the mitochondrion/hydrogenosome-like organelles.


    Materials and Methods
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
Organisms and Genomic DNA Preparation
Trichomonas vaginalis strain NIH-C1 (ATCC 30001) and G. intestinalis strain WB, clone 6 (ATCC 300957), were used. Genomic DNA was isolated from T. vaginalis using a guanidium thiocyanate procedure (Wang and Wang 1985Citation ) and from G. intestinalis using a Blood & Culture DNA kit (Qiagen, Chatsworth, Calif.).

Probe Preparation, Cloning, and Screening of the Genomic Library
To obtain probes for screening a T. vaginalis genomic library, two pairs of degenerate primers, a GC-rich one and an AT-rich one, were designed based on the conserved regions of IscS/NifS sequences in GenBank (National Center for Biomedical Information): EIIFTSGATE (GC-rich sense: 5'-GARATYATYTTCTCVTCHGGHGCHACHGAR-3'; AT-rich sense: 5'-GAAATWATWTTYACWWSWGGWGCWACWGAA) and HKIH/YGPKGV/IG (GC-rich antisense: 5'-CCRAYDCCYTTTGGDCCRTRRATYTTRTG-3'; AT-rich antisense: 5'-CCWAYWCCYTT TGGWCCRTRWATYTTRTG-3'). Corresponding fragments were amplified by PCR, purified with a gel extraction kit (Qiagen) and cloned into pCR 2.1 vector (TA cloning kit, Invitrogen). The inserts were excised from the vector, gel-purified, and labeled by means of a Random Primers DNA Labeling System (GIBCO/BRL) with {alpha}-[32P]dATP. These probes were used for screening a genomic DNA library in {lambda} ZAP II vector (Stratagene). The sequences of positive clones were determined for both strands by primer walking.

Nucleotide sequences of Escherichia coli and Saccharomyces cerevisiae IscS/NifS homolog genes were used to search the Giardia lamblia genome sequence database (http://www.mbl.edu/baypaul/Giardia-HTML/index2.html; McArthur et al. 2000Citation ) with the BLAST program. Clones Ai0824 and Ki1686 contained sequences homologous to the N- and C-terminal ends of the bacterial and eukaryotic homologs. Based on the nucleotide sequence of these clones, we designed a pair of oligonuclueotide primers (sense: 5'-GATGACGAGCGTGCAAGGAAAGCTC-3'; antisense: 5'-GGTGACTACATGCGGATGCTCAGCC-3') located in the 5' and 3' untranslated regions of the putative G. intestinalis nifS homolog gene, respectively. PCR reactions, utilizing these oligonucleotides and G. intestinalis genomic DNA as template, amplified a 2.4-kb fragment that was purified, cloned into the pCR2.1 vector (Invitrogen), and sequenced.

Sequence Alignment
Nucleotide and protein database searches were performed at the National Center for Biomedical Information using the BLAST program (Altschul et al. 1997Citation ). Sequences were extracted from databases using the BlastAli program (http://www.joern-lewin.de/). The IscS sequences of T. vaginalis and G. intestinalis were aligned to sequences from 64 taxa using ClustalX (Thompson et al. 2000Citation ). The alignment was further edited visually with the use of the ED program of MUST (Philippe 1993Citation ). The alignment of all 67 taxa resulted in 231 shared amino acid positions, while an alignment of 21 selected taxa consisted of 362 shared amino acid positions. The T. vaginalis TviscS-1, T. vaginalis TviscS-2, and G. intestinalis GiiscS sequences have been submitted to GenBank under accession numbers AF321005, AF321006, and AF311744, respectively.

Phylogenetic Analysis
Phylogenetic relationships were analyzed by means of the Neighbor-Joining (NJ) and Maximum-Parsimony (MP) methods using PHYLIP, version 3.6 (Felsenstein 1989Citation ), and by the Maximum-Likelihood (ML) method using the PROTML program in MOLPHY, version 2.3 (Adachi and Hasegawa 1996Citation ). The ML tree was constructed by local rearrangement of an NJ tree using the Jones-Taylor-Thornton model of amino acid substitutions with the F-option (JTT-F) to account for amino acid frequencies in the data set. User-defined trees were analyzed to compare alternative topologies (Kishino and Hasegawa 1989Citation ). The local bootstrap proportion value was calculated for each internal branch of the ML tree using a local rearrangement option of the PROTML program. Bootstrap support for distance and parsimony analyses were based on 100 resampled data sets using SEQBOOT, PHYLIP, version 3.6.


    Results and Discussion
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
Analysis of T. vaginalis and G. intestinalis iscS Genes and Putative Translation Products
Sequences for two complete iscS genes from T. vaginalis (TviscS-1 and TviscS-2) and for one gene from G. intestinalis (GiiscS) were obtained. The T. vaginalis genomic DNA library was screened with two probes derived from PCR products. The products were amplified using the A+T-rich or the G+C-rich degenerate primer pairs. Each pair amplified a distinct DNA fragment of about 4.4 kb. The two products displayed only 60% nucleotide sequence identity. Both fragments were identified as IscS/NifS homologs by a BLAST search. These products used as probes recognized separate sets of positive clones of the genomic library containing either TviscS-1 or TviscS-2 genes. No cross-reactivity between these clones was observed. The isolated clones contained complete putative open reading frames (ORFs) without intron-like sequences. The G+C content of TviscS-1 was rather low (39.7%), and that of TviscS-2 was higher (46.7%). The TviscS-1 and TviscS-2 coded for proteins 385 and 411 amino acids in length, respectively. Predicted molecular mass and an isoelectric point for the putative products for TviscS-1 were 41.9 and 8.2, while they were 44.8 and 8.2 for TviscS-2, respectively. The G. intestinalis gene isolated contained a 52% G+C-rich ORF coding for a polypeptide of 433 amino acids with a predicted molecular mass of 47.7 and an isoelectric point of 6.3.

Amino acid sequences deduced from the Trichomonas and Giardia genes were compared with IscS/NifS homologs from 64 species, including bacteria, fungi, plants, invertebrates, and vertebrates (alignment available on request from J.T.). An alignment of selected sequences including that of eubacterial IscS from A. vinelandii and mitochondrial sequence from S. cerevisiae is shown in figure 1 . Both Trichomonas and Giardia sequences contained all conserved regions proposed to mediate the cysteine desulfurase activity in IscS/NifS-like proteins: (1) His111 (numbered according to TviscS-1), which is involved in initial deprotonation of the substrate (Kaiser et al. 2000Citation ); (2) the pyridoxal-5'-phosphate-binding site with the Schiff base forming Lys222 residue and Asp187 and Gln190, which bind the pyridine nitrogen and the phenolate oxygen of PLP, respectively, and residues involved in forming an additional six hydrogen bonds anchoring the phosphate group: Thr82, His221, Ser/Thr219, and Thr250 (Zheng et al. 1993Citation ); and (3) the substrate-binding site including Cys371, which provides a reactive cysteinyl residue (Zheng et al. 1994Citation ), as well as Arg362, Asn162, and Asn 41, which anchor the cysteine with a salt bridge and hydrogen bond (Kaiser et al. 2000Citation ).



View larger version (64K):
[in this window]
[in a new window]
 
Fig. 1.—Sequence alignment of putative Trichomonas vaginalis (TviscS-1, TviscS-2) and Giardia intestinalis (GiiscS) IscS proteins with eubacterial (Azotobacter vinelandii, Av) and mitochondrial (Saccharomyces cerevisiae, Sc) homologs. Conserved lysine and other residues involved in PLP binding are indicated by closed () and open circles ({circ}), respectively. Invariable cysteine is highlighted by a closed square ({block}), while other residues involved in substrate binding are indicated by open squares ({square}). An arrow indicates the conserved histidine involved in the substrate deprotonation. Putative N-terminal presequences are underlined. Cystein signatures of eukaryotic IscS and C-terminal conserved residues typical for eukaryotic/eubacterial IscS are boxed. ("*" indicates fully conserved residue, ":" indicates conserved "strong" groups, and "." indicates conserved "weaker" groups according to ClustalX). For accession numbers see figure 3

 
Importantly, the alignment revealed short N-terminal extensions of the trichomonad IscS sequences which were not present at the N-terminus of IscS in Giardia or in the eubacterial sequences. The PSORT II program (http://psort.nibb.ac.jp/), designed for prediction of cleavage sites for mitochondrial presequences, recognized the SRS/YF motif with a characteristic arginine at position -2 relative to the cleavage site at the TviscS-2 N-terminus (fig. 2 ). The TviscS-2 extension was typically serine-rich, consisting of eight amino acids, which started with leucine. It resembled N-terminal leader sequences found in proteins targeted to the Trichomonas hydrogenosome and can be assumed to have the same function (fig. 2 ). Interestingly, TviscS-1 also started at the same position and contained an eight-amino-acid N-terminal extension; however, it did not contain the typical -2 arginine, and PSORT II did not recognize the consensus cleavage sequence. Nevertheless, a function of TviscS-1 presequence in organelle targeting cannot be ruled out, as a presequence without arginine has recently been reported for the hydrogenosomal membrane protein Hmp31 (Dyall et al. 2000Citation ). No N-terminal extension of IscS was expected in Giardia, as this organism contains neither mitochondria nor hydrogenosomes. Lack of the N-terminal extension was previously reported in Giardia Hsp60, another mitochondrial-type protein recognized in this organism (Roger 1998Citation ).



View larger version (24K):
[in this window]
[in a new window]
 
Fig. 2.—Comparison of N-terminal amino acid sequences of TviscS-1 and TviscS-2 with leader sequences known to target hydrogenosomal proteins to the organelle. The cleavage sites recognized by PSORT program (http://psort.nibb.c.jp/) are indicated by thin arrows, and those determined by N-terminal sequencing (Dyall et al. 2000Citation ) are indicated by thick arrows. Accession numbers of compared genes are as follows: adenylate kinase, P49983; ferredoxin, P21149; heat shock protein 60, Q95058; malic enzyme A, AAA92714; succinylCoA synthase {alpha}SCS1, P53399; succinylCoA synthase ßSCS, Q03184; pyruvate : ferredoxin oxidoreductase A, AAA85494; hydrogenosomal membrane protein Hmp31, AF216971

 
The presence of all key elements required for cysteine desulfurase activity in the deduced proteins suggests that IscS homologs are involved in Fe-S cluster assembly in T. vaginalis and G. intestinalis. The presence of the conserved N-terminal leader sequence in TviscS-2 suggests that in trichomonads the IscS-like proteins are targeted to hydrogenosomes, which are the likely sites of Fe-S cluster assembly. In Giardia, the FeS clusters are possibly assembled in cytosol.

A C-terminal sequence signature differentiates proteobacterial and eukaryotic IscS from homologs in all other organisms. TviscS-2 and GiiscS also contain this signature, which consists of 20 to 21 amino acids with consensus sequence SPL(W/Y)(E/D)(M/L)X(K/Q)XG(I/V)D(L/I)XX(I/V)XWXXX (fig. 1 ). NifS genes of nitrogen-fixing bacteria also possess a similar C-terminal extension that starts with the Ser-Pro motif, but the subsequent sequence is not conserved. The only sequence from a eukaryote that lacks this signature is TviscS-1. The lack of this signature and the "atypical" N-terminal extension might indicate that Tvisc-1 is a pseudogene or that its product has a different function or localization. Eukaryotic IscS is distinguished from all other organisms, including proteobacteria, by the invariable Cys113 in the substrate deprotonation region. This residue is present in both trichomonad iscS products as well as in GiiscS. Prokaryotes possess Ala, Ser, or Gly at this position. Interestingly, Giardia IscS possessed two unique highly hydrophilic inserts, Thr137-Glu145 and Glu300-Ser321, which are not present in any of the 66 other species. It will be of interest to determine the function of such inserts, which might be associated with a specific localization of the gene product in Giardia.

Our sequencing data provide a solid basis for the prediction of the function of the products of iscS genes in amitochondriate eukaryotes. However, further studies are required to confirm their physiological function and cell localization.

Phylogenetic Analysis of IscS/NifS-like Sequences
In all global phylogenetic reconstructions, IscS/NifS-like homologs formed two distinct groups that were previously designated groups I and II (Mihara et al. 1997Citation ). The IscS sequences of Trichomonas and Giardia and those of the mitochondrial homologs in other eukaryotes formed a single clade (group I) with a high bootstrap value (99%) using the local rearrangement option of the PROTML program (fig. 3 ). Within this clade, Trichomonas and Giardia formed a subgroup together with Plasmodium falciparum and Arabidopsis thaliana. The second eukaryotic subgroup consisted of metazoan IscS, and the third group comprised homologs in fungal mitochondria. The {alpha}-proteobacterium Rickettsia prowazeki, often considered a close relative to the mitochondrial ancestor, clustered together with metazoan mitochondrial IscS (Andersson and Kurland 1999Citation ).



View larger version (47K):
[in this window]
[in a new window]
 
Fig. 3.—Phylogenetic tree based on analysis of IscS/NifS amino acid sequences. The tree was constructed using the maximum-likelihood method with the Jones-Taylor-Thornton model of amino acid substitution (JTT-F). Local bootstrap values calculated by the PROTML program are attached to the internal nodes. The cell localization is indicated by M (mitochondria), H (hydrogenosomes), and C (cytosol). * These sequence data were produced by the Plasmodium falciparum Sequencing Group at the Sanger Centre and can be obtained from ftp://ftp.sanger.ac.uk/pub/pathogens/malaria2/unfinished_ORFS/orf_dna_sequences

 
The IscS of proteobacteria formed a sister group to eukaryotes. It is known that nitrogen-fixing bacteria possess two separate gene clusters that are involved in FeS cluster assembly: products of the nif gene cluster, including NifS, participate in nitrogenase-specific FeS cluster formation (Zheng and Dean 1994Citation ), and isc gene products, including IscS, which have a more general function in the assembly of other Fe-S proteins (Zheng et al. 1998Citation ). Our analysis suggests that isc-related genes were present in the proteobacterial endosymbiotic ancestor of mitochondria and hydrogenosomes. The more specialized NifS sequences formed a distinct group with 99% local bootstrap support in the PROTML analysis which contained exclusively nitrogen-fixing proteobacteria.

Additional IscS/NifS-like homologs of eukaryotic organisms (Mus musculus and A. thaliana) also appeared in another two distinct clades. The mouse homolog was located in a heterogeneous eubacterial group (fig. 3 ). This gene codes for the cytosolic pyridoxal-5'-phosphate-dependent selenocysteine lyase, which resembles NifS in primary structure as well as in catalytic function (Mihara et al. 2000Citation ). If the partial sequence of a human counterpart was also included in the analysis (data not shown), it was a sister group to the mouse sequence. The function of bacterial gene products of this clade has not been studied except in Synechocystis sp. (S74526 corresponds to sll0704 in CyanoBase http://www.kauza.or.jp/cyano/ investigated by Kato et al. [2000]Citation ) and Bacillus subtilis (AAA21613 corresponds to the nifS-like gene according to Sun and Setlow [1993]Citation ). The product of the Synechocystis gene showed selenocysteine lyase activity, although it also acted on L-cysteine sulfinic acid and other substrates (Kato et al. 2000Citation ). The B. subtilis gene product has been suggested to participate in NAD biosynthesis (Sun and Setlow 1993Citation ). Thus, it is likely that other members of this clade also have biochemical functions that are different from those of NifS and IscS proteins. The topology of A. thaliana genes was of particular interest. While one gene was related to the subtree of genes coding for mitochondrial IscS in protists, a second gene was placed in group II together with Synechocystis S76601. Since cyanobacteria share a common ancestor with plastids, we analyzed the second A. thaliana sequence for its possible subcellular localization with the PSORT program. The analysis gave the highest score for a chloroplast stroma localization (certainty = 0.501) of the gene product. This analysis suggests that the second IscS homolog of A. thaliana may operate in the chloroplast. However, it is difficult to predict its possible function. Group II consists of the most divergent NifS/IscS homologs. Function has been established only for products of two E. coli genes, which encode cysteine sulfinate desulfinase (F65063; Mihara et al. 1997Citation ), and for selenocysteine lyase (H64925; Fujii et al. 2000Citation ). Thus, the A. thaliana IscS homolog, as well as other members of this heterogenous bacterial group, might have functions different from FeS cluster formation, and their distance from genes of group I might reflect different evolutionary pressures. In any case, further biochemical studies on the functions of group II members are required.

The global gene tree showed that amitochondrial and mitochondrial IscSs share a common eubacterial ancestor, suggesting a common biosynthetic mechanism for FeS proteins. The tree also showed that both eukaryotes and bacteria possess several paralogous or orthologous IscS/NifS-like genes. The topology of these gene trees possibly reflects their specialized function or cell localization more than their large-scale phylogenetic relationships. Therefore, in a subsequent analysis we restricted the data set to eukaryotic and proteobacterial IscSs. The trees constructed by the ML, MP, and NJ methods confirmed the close relationship between amitochondrial and mitochondrial IscS with high bootstrap support (fig. 4 ). The robustness of the relationship within the eukaryotic group was further assessed through the analysis of alternative tree topologies. We defined six branches on the tree: (1) amitochondriates and Plasmodium, (2) Arabidopsis, (3) Fungi, (4) Metazoa, (5) Rickettsia, and (6) Proteobacteria. Evaluation of the 105 alternative trees confirmed a common ancestry of genes from amitochondriate and mitochondriate organisms. However, the branching order within the eukaryotic clade was not resolved, as several alternative positions for eukaryotic subtrees with comparable significance were found (table 1 ). We further suspected that the subtree of protists could be affected by long-branch attraction in spite of high bootstrap support. Thus, the positions of the Trichomonas and Giardia genes were tested using a data set from which we removed the most divergent sequences of A. thaliana and P. falciparum. The analysis clearly showed the instability of the protist group. Although NJ and MP reconstruction placed Trichomonas and Giardia as sister taxa with bootstrap support of 85% and 87%, respectively, ML constraint analysis gave comparable support to several alternative hypotheses. In the best tree topology and in four other tree topologies with {Delta}ln L <1 SE from the best tree, Trichomonas and Giardia were not placed as sister taxa (table 1 ). Nevertheless, in all alternative trees, both amitochondriates were part of the eukaryotic clade.



View larger version (31K):
[in this window]
[in a new window]
 
Fig. 4.—Tree of eukaryotic-proteobacterial IscS constructed with the Maximum-Likelihood (ML) method (Jones-Taylor-Thornton model). Bootstrap values (local bootstrap proportion) of the PROTML analysis are shown at the considered nodes. For the node of the main interest, the bootstrap values calculated by the ML, Maximum Parsimony (MP), and Neighbor-Joining (NJ) methods are given in a box

 

View this table:
[in this window]
[in a new window]
 
Table 1 Maximum-Likelihood Analysis of Alternative Tree Topologies Using the JTT-F Model of Amino Acid Substitution

 
Our results indicate that a common mechanism mediates FeS cluster assembly in mitochondriate and amitochondriate eukaryotes, even though different sets of FeS proteins function in these organisms (Müller 1998Citation ). Both mitochondriate and amitochondriate eukaryotes require FeS proteins for single-electron transport processes associated with ATP production linked to the utilization of oxygen (mitochondriates) or not (type I and II amitochondriates). In mitochondria, FeS proteins are involved in the respiratory chain (subunits of complexes I, II, and III) and in the citric acid cycle (aconitase), while low–redox-potential FeS proteins such as pyruvate : ferredoxin oxidoreductase and ferredoxin are present in the hydrogenosomes of type I and in the cytosol of type II amitochondriates. Hydrogenosomes also contain hydrogenase. The presence of a common machinery responsible for FeS cluster formation and incorporation into "aerobic" and "ananerobic" types of apoproteins is supported by (1) the identification of the genes encoding IscS, a key member of machinery involved in FeS cluster assembly, in T. vaginalis and G. intestinalis, and (2) the placement of IscS of these amitochondrial organisms within the eukaryotic clade in phylogenetic analysis. In addition, a partial sequence of IscU, another protein involved in FeS cluster formation, has been identified in the database of the Giardia genome project (http://www.mbl.edu/baypaul/Giardia-HTML/index2.html).

A common origin of FeS cluster formation in mitochondriate and amitochondriate eukaryotes could be explained by the recently proposed hydrogen hypothesis of eukaryotic origin (Martin and Müller 1998Citation ). The hypothesis assumes that all eukaryotes, including contemporary amitochondrial organisms, once harbored the mitochondrion/hydrogenosome-like organelle derived from a proteobacterial endosymbiont. The ancestral endosymbiont is viewed as a facultatively anaerobic proteobacterium which possessed both anaerobic and aerobic metabolic machineries for electron transport–linked ATP production, including "aerobic" and "anaerobic" types of FeS proteins. A possible scenario is that the "anaerobic" set of FeS proteins was preserved in hydogenosomes, whereas the "aerobic" set was preserved in mitochondria. Both organelles inherited the common mechanism of the FeS cluster assembly. This scenario is supported by the close phylogenetic relationship between eukaryotic and proteobacterial IscS proteins. Our hypothesis is also congruent with a common origin of mitochondria and Trichomonas hydrogenosomes, as well as a secondary loss of the mitochondrion/hydrogenosome-like organelle in Giardia. We cannot, however, rule out alternative explanations for the origin of hydrogenosomes and biochemistry of amitochondrial eukaryotes, including the mechanism of FeS cluster assembly. Indeed, independent lateral gene transfers (Doolittle 1998Citation ) or preservation of certain biochemical pathways from an anaerobic past of eukaryotic evolution might be involved. Nevertheless, comparative analysis of mechanisms responsible for the formation of FeS clusters, which are considered to be among the most ancient biologically active metal cofactors (Cammack 1996Citation ), appears to be a promising tool for tracing eukaryotic history.


    Acknowledgements
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
We thank Dr. Tetsuo Hashimoto for discussions and critical reading of the manuscript and Dr. Gang Wu for encouragement and help. We also thank Dr. Mitchell L. Sogin and Dr. Andrew G. McArthur for the use of the Giardia genome database. Oligonucleotide synthesis and DNA sequencing were performed at the Protein/DNA Technology Center at Rockefeller University. This work was supported by National Institutes of Health grant AI11942 to M.M. and by a grant from the Grant Agency of the Czech Republic to J.T. (204/00/1561). A major part of this work was performed during J.T.'s stay at the New York laboratory, supported by the NATO Science Fellowships Program.


    Footnotes
 
Geoffrey McFadden, Reviewing Editor

1 Abbreviations: ML, maximum likelihood; MP, maximum parsimony; NJ, neighbor-joining; PLP, pyridoxal-5'-phosphate. Back

2 Keywords: IscS iron-sulfur cluster hydrogenosome Trichomonas vaginalis Giardia intestinalis Back

3 Address for correspondence and reprints: Jan Tachezy, Department of Parasitology, Faculty of Science, Charles University, Vinic 7, Prague 128 44, Czech Republic. tachezy{at}natur.cuni.cz . Back


    References
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 

    Adachi J., M. Hasegawa, 1996 MOLPHY version 2.3: programs for molecular phylogenetics based on maximum likelihood Comput. Sci. Monogr 28:1-150

    Altschul S. F., T. L. Madden, A. A. Schaffer, J. Zhang, Z. Zhang, W. Miller, D. J. Lipman, 1997 Gapped BLAST and PSI-BLAST: a new generation of protein database search programs Nucleic Acids Res 25:3389-3402[Abstract/Free Full Text]

    Amrani L., J. Primus, A. Glatigny, L. Arcangeli, C. Scazzocchio, V. Finnerty, 2000 Comparison of the sequences of the Aspergillus nidulans hxB and Drosophila melanogaster ma-I genes with nifS from Azotobacter vinelandii suggests a mechanism for the insertion of the terminal sulphur atom in the molybdopterin cofactor Mol. Microbiol 38:114-125[ISI][Medline]

    Andersson S. G., C. G. Kurland, 1999 Origins of mitochondria and hydrogenosomes Curr. Opin. Microbiol 2:535-541[ISI][Medline]

    Benchimol M., P. J. Johnson, W. De Souza, 1996 Morphogenesis of the hydrogenosome: an ultrastructural study Biol. Cell 87:197-205[ISI][Medline]

    Bradley P. J., C. J. Lahti, E. Plümper, P. J. Johnson, 1997 Targeting and translocation of proteins into the hydrogenosome of the protist Trichomonas: similarities with mitochondrial protein import EMBO J 16:3484-3493[Abstract/Free Full Text]

    Bui E. T., P. J. Bradley, P. J. Johnson, 1996 A common evolutionary origin for mitochondria and hydrogenosomes Proc. Natl. Acad. Sci. USA 93:9651-9656[Abstract/Free Full Text]

    Cammack R., 1996 Iron and sulfur in the origin and evolution of biological energy conversion systems Pp. 43–69 in H. Baltscheffsky, ed. Origin and evolution of biological energy conversion. VCH, New York

    Cavalier-Smith T., 1987 The origin of eukaryotic and archaebacterial cells Ann. N.Y. Acad. Sci 503:17-54[ISI][Medline]

    Clark C. G., A. J. Roger, 1995 Direct evidence for secondary loss of mitochondria in Entamoeba histolytica Proc. Natl. Acad. Sci. USA 92:6518-6521[Abstract]

    Doolittle W. F., 1998 You are what you eat: a gene transfer ratchet could account for bacterial genes in eukaryotic nuclear genomes Trends Genet 14:307-311[ISI][Medline]

    Dyall S. D., C. M. Koehler, M. G. Delgadillo-Correa, P. J. Bradley, E. Plümper, D. Leuenberger, C. W. Turck, P. J. Johnson, 2000 Presence of a member of the mitochondrial carrier family in hydrogenosomes: conservation of membrane-targeting pathways between hydrogenosomes and mitochondria Mol. Cell. Biol 20:2488-2497[Abstract/Free Full Text]

    Ellis J. E., R. Williams, D. Cole, R. Cammack, D. Lloyd, 1993 Electron transport components of the parasitic protozoon Giardia lamblia FEBS Lett 325:196-200[ISI][Medline]

    Embley T. M., R. P. Hirt, 1998 Early branching eukaryotes? Curr. Opin. Genet. Dev 8:624-629[ISI][Medline]

    Felsenstein J., 1989 PHYLIP—phylogeny inference package (version 3.2) Cladistics 5:164-166

    Fujii T., M. Maeda, H. Mihara, T. Kurihara, N. Esaki, Y. Hata, 2000 Structure of a NifS homologue: X-ray structure analysis of CsdB, an Escherichia coli counterpart of mammalian selenocysteine lyase Biochemistry 39:1263-1273[ISI][Medline]

    Hackstein J. H., A. Akhmanova, B. Boxma, H. R. Harhangi, F. G. Voncken, 1999 Hydrogenosomes: eukaryotic adaptations to anaerobic environments Trends Microbiol 7:441-447[ISI][Medline]

    Hrd I., M. Müller, 1995a Primary structure and eubacterial relationships of the pyruvate : ferredoxin oxidoreductase of the amitochondriate eukaryote Trichomonas vaginalis J. Mol. Evol 41:388-396[ISI][Medline]

    ———. 1995b Primary structure of the hydrogenosomal malic enzyme of Trichomonas vaginalis and its relationship to homologous enzymes J. Eukaryot. Microbiol 42:593-603[ISI][Medline]

    Johnson P. J., C. J. Lahti, P. J. Bradley, 1993 Biogenesis of the hydrogenosome in the anaerobic protist Trichomonas vaginalis J. Parasitol 79:664-670[ISI][Medline]

    Kaiser J. T., T. Clausen, G. P. Bourenkow, H. D. Bartunik, S. Steinbacher, R. Huber, 2000 Crystal structure of a NifS-like protein from Thermotoga maritima: implications for iron sulphur cluster assembly J. Mol. Biol 297:451-464[ISI][Medline]

    Kambampati R., C. T. Lauhon, 1999 IscS is a sulfurtransferase for the in vitro biosynthesis of 4-thiouridine in Escherichia coli tRNA Biochemistry 38:16561-16568[ISI][Medline]

    Kato S.-I., H. Mihara, T. Kurihara, T. Yoshimura, N. Esaki, 2000 Gene cloning, purification, and characterization of two cyanobacterial NifS homologs driving iron-sulfur cluster formation Biosci. Biotechnol. Biochem 64:2412-2419[ISI][Medline]

    Kishino H., M. Hasegawa, 1989 Evaluation of the maximum likelihood estimate of the evolutionary tree topology from DNA sequence data, and the branching order in hominoidea J. Mol. Evol 29:170-179[ISI][Medline]

    Kishino H., T. Miyata, M. Hasegawa, 1990 Maximum likelihood inference of protein phylogeny and the origin of chloroplast J. Mol. Evol 30:151-160

    Kispal G., P. Csere, C. Prohl, R. Lill, 1999 The mitochondrial proteins Atm1p and Nfs1p are essential for biogenesis of cytosolic Fe/S proteins EMBO J 18:3981-3989[Abstract/Free Full Text]

    Kolman C., D. Söll, 1993 SPL1-1, a Saccharomyces cerevisiae mutation affecting tRNA splicing J. Bacteriol 175:1433-1442[Abstract]

    Kulda J., 1999 Trichomonads, hydrogenosomes and drug resistance Int. J. Parasitol 29:199-212[ISI][Medline]

    Land T., T. A. Rouault, 1998 Targeting of a human iron-sulfur cluster assembly enzyme, nifs, to different subcellular compartments is regulated through alternative AUG utilization Mol. Cell 2:807-815[ISI][Medline]

    Länge S., C. Rozario, M. Müller, 1994 Primary structure of the hydrogenosomal adenylate kinase of Trichomonas vaginalis and its phylogenetic relationships Mol. Biochem. Parasitol 66:297-308[ISI][Medline]

    Lauhon C., R. Kambampati, 2000 The iscS gene in Escherichia coli is required for the biosynthesis of 4-thiouridine, thiamin, and NAD J. Biol. Chem 275:20096-20103[Abstract/Free Full Text]

    Lill R., G. Kispal, 2000 Maturation of cellular Fe-S proteins: an essential function of mitochondria Trends Biochem. Sci 25:352-356[ISI][Medline]

    McArthur A. G., H. G. Morrison, J. E. Nixon, et al. (15 co-authors) 2000 The Giardia genome project database FEMS Microbiol Lett 15:271-273

    Mai Z., S. Ghosh, M. Frisardi, B. Rosenthal, R. Rogers, J. Samuelson, 1999 Hsp60 is targeted to a cryptic mitochondrion-derived organelle (‘crypton’) in the microaerophilic protozoan parasite Entamoeba histolytica Mol. Cell. Biol 19:2198-2205[Abstract/Free Full Text]

    Martin W., M. Müller, 1998 The hydrogen hypothesis for the first eukaryote Nature 392:37-41[ISI][Medline]

    Mihara H., T. Kurihara, T. Watanabe, T. Yoshimura, N. Esaki, 2000 cDNA cloning, purification, and characterization of mouse liver selenocysteine lyase. Candidate for selenium delivery protein in selenoprotein synthesis J. Biol. Chem 275:6195-6200[Abstract/Free Full Text]

    Mihara H., T. Kurihara, T. Yoshimura, K. Soda, N. Esaki, 1997 Cysteine sulfinate desulfinase, a NIFS-like protein of Escherichia coli with selenocysteine lyase and cysteine desulfurase activities. Gene cloning, purification, and characterization of a novel pyridoxal enzyme J. Biol. Chem 272:22417-22424[Abstract/Free Full Text]

    Müller M., 1993 The hydrogenosome J. Gen. Microbiol 139:2879-2889[ISI][Medline]

    ———. 1997 Evolutionary origins of trichomonad hydrogenosomes Parasitol. Today 13:166-167[ISI]

    ———. 1998 Enzymes and compartmentation of core energy metabolism of anaerobic protists—a special case in eukaryotic evolution? Pp. 109–131 in G. H. Coombs, K. Vickerman, M. A. Sleigh, and A. Warren, eds. Evolutionary relationship among protozoa. Kluwer, Dordrecht, The Netherlands

    ———. 2000 A mitochondrion in Entamoeba histolytica? Parasitol. Today 16:368-369[ISI][Medline]

    Nakai Y., Y. Yoshihara, H. Hayashi, H. Kagamiyama, 1998 cDNA cloning and characterization of mouse nifS-like protein, m-Nfs1: mitochondrial localization of eukaryotic NifS-like proteins FEBS Lett 433:143-148[ISI][Medline]

    Philippe H., 1993 MUST, a computer package of management utilities for sequences and trees Nucleic Acids Res 21:5264-5272[Abstract]

    Reeves R. E., 1984 Metabolism of Entamoeba histolytica Schaudinn, 1903 Adv. Parasitol 23:105-142[ISI][Medline]

    Roger A. J., 1999 Reconstructing early events in eukaryotic evolution Am. Nat 154:S146-S163[ISI][Medline]

    Roger A. J., S. G. Svärd, J. Tovar, C. G. Clark, M. W. Smith, F. D. Gillin, M. L. Sogin, 1998 A mitochondrial-like chaperonin 60 gene in Giardia lamblia: evidence that diplomonads once harbored an endosymbiont related to the progenitor of mitochondria Proc. Natl. Acad. Sci. USA 95:229-234[Abstract/Free Full Text]

    Schwartz C. J., O. Djaman, J. A. Imlay, P. J. Kiley, 2000 The cysteine desulfurase, IscS, has a major role in in vivo Fe-S cluster formation in Escherichia coli Proc. Natl. Acad. Sci. USA 97:9009-9014[Abstract/Free Full Text]

    Strain J., C. R. Lorenz, J. Bode, S. Garland, G. A. Smolen, D. T. Ta, L. E. Vickery, V. C. Culotta, 1998 Suppressors of superoxide dismutase (SOD1) deficiency in Saccharomyces cerevisiae. Identification of proteins predicted to mediate iron-sulfur cluster assembly J. Biol. Chem 273:31138-31144[Abstract/Free Full Text]

    Sun D., P. Setlow, 1993 Cloning, nucleotide sequence, and regulation of the Bacillus subtilis nadB gene and a nifS-like gene, both of which are essential for NAD biosynthesis J. Bacteriol 175:1423-1432[Abstract]

    Thompson J. D., T. J. Gibson, F. Plewniak, F. Jeanmougin, D. G. Higgins, 2000 The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools Nucleic Acids Res 25:4876-4882[Abstract/Free Full Text]

    Tovar J., A. Fischer, C. G. Clark, 1999 The mitosome, a novel organelle related to mitochondria in the amitochondrial parasite Entamoeba histolytica Mol. Microbiol 32:1013-1021[ISI][Medline]

    Wang A. L., C. C. Wang, 1985 Isolation and characterization of DNA from Tritrichomonas foetus and Trichomonas vaginalis Mol. Biochem. Parasitol 14:323-335[ISI][Medline]

    Zheng L., V. L. Cash, D. H. Flint, D. R. Dean, 1998 Assembly of iron-sulfur clusters. Identification of an iscSUA-hscBA-fdx gene cluster from Azotobacter vinelandii J. Biol. Chem 273:13264-13272[Abstract/Free Full Text]

    Zheng L., D. R. Dean, 1994 Catalytic formation of a nitrogenase iron-sulfur cluster J. Biol. Chem 269:18723-18726[Abstract/Free Full Text]

    Zheng L., R. H. White, V. L. Cash, D. R. Dean, 1994 Mechanism for the desulfurization of L-cysteine catalyzed by the nifS gene product Biochemistry 33:4714-4720[ISI][Medline]

    Zheng L., R. H. White, V. L. Cash, R. F. Jack, D. R. Dean, 1993 Cysteine desulfurase activity indicates a role for NIFS in metallocluster biosynthesis Proc. Natl. Acad. Sci. USA 90:2754-2758[Abstract]

Accepted for publication June 15, 2001.