This Article Full Text (PDF) Supplemental Data Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via ISI Web of Science (18) Citing Articles via Google Scholar Google Scholar Articles by Hanikenne, M. Articles by Baurain, D. Search for Related Content PubMed PubMed Citation Articles by Hanikenne, M. Articles by Baurain, D. Agricola Articles by Hanikenne, M. Articles by Baurain, D.

UPDATE ON METAL TRANSPORTERS IN UNICELLULAR ALGAE

A Comparative Inventory of Metal Transporters in the Green Alga Chlamydomonas reinhardtii and the Red Alga Cyanidioschizon merolae^1,[w]

Metal Homeostasis Group, Max Planck Institute for Plant Molecular Physiology, 14476 Golm, Germany (M.H., U.K.); and Algology, Mycology and Experimental Systematics, Department of Life Sciences, University of Liege, Sart Tilman, 4000 Liege, Belgium (V.D., D.B.)

As in all organisms, metal cations are crucial for nutrition in plants. Several metals, such as copper, iron, zinc, and manganese, act as important cofactors for many enzymes and are essential for both mitochondrial and chloroplast functions. However, when supplied in excess, these essential cations can become toxic, like heavy metals with no generally established function, such as cadmium, lead, or mercury. To maintain micronutrient metal homeostasis and to cope with the deleterious effects of nonessential heavy metals, plants have developed a complex network of metal uptake, chelation, trafficking, and storage processes. Metal transporters are required to maintain metal homeostasis and thus constitute important components of this network (Clemens, 2001; Hall and Williams, 2003).

In recent years, a number of membrane transport protein families have been implicated in metal homeostasis in plants. These include the cation diffusion facilitators (CDF), the Zrt-, Irt-like proteins (ZIP), the cation exchangers (CAX), the copper transporters (COPT), the heavy-metal P-type ATPases (HMA), the natural resistance-associated macrophage proteins (NRAMP), and the ATP-binding cassette (ABC) transporters (Williams et al., 2000; Maser et al., 2001; Cobbett et al., 2003; Hall and Williams, 2003). These transporters are encoded by multigene families. For example, 15 ZIP genes, 12 metal tolerance protein (MTP) genes, and 8 HMA genes are present in the Arabidopsis (Arabidopsis thaliana) genome (Cobbett et al., 2003; Delhaize et al., 2003). However, the transport specificities, patterns of expression, or subcellular localizations of metal transport proteins are still largely unknown. To further understand plant metal homeostasis, it will be necessary to elucidate the contribution of each of these transporters to the uptake, trafficking, and storage of essential metals, as well as to the detoxification of toxic heavy metals.

In this article, we present an overview of our current knowledge of the metal transport function in metal homeostasis and tolerance in eukaryotes, with a special emphasis on plants. We also provide a timely inventory of putative metal transporters in two unicellular algal models, the green alga Chlamydomonas reinhardtii and the red alga Cyanidioschizon merolae. These new data should facilitate functional genomics and molecular analysis of metal homeostasis and tolerance in photosynthetic organisms. Moreover, the comparison of metal transporters from species belonging to the red and green algae with those of the land plant Arabidopsis, as well as with their human and yeast homologs, allows some light to be shed on the molecular evolution of metal homeostasis and tolerance systems.

	TWO UNICELLULAR ALGAL MODEL SYSTEMS

TOP TWO UNICELLULAR ALGAL MODEL... A MINING STRATEGY SUITED... THE CDF FAMILY THE ZIP FAMILY THE CAX FAMILY THE COPT FAMILY THE HMA SUBFAMILY OF... THE ABC TRANSPORTER FAMILY IRON TRANSPORTERS METHODOLOGICAL INSIGHTS PHYLOGENETIC INSIGHTS CONCLUSION LITERATURE CITED

 
The use of a simple model organism is often helpful to dissect the functions and the interactions of various transport systems. The yeasts Saccharomyces cerevisiae and, to a lesser extent, Schizosaccharomyces pombe have been successfully used to elucidate the molecular basis of cellular metal homeostasis (Clemens and Simm, 2003; De Freitas et al., 2003; Eide, 2003). However, investigation of the essential chloroplast metal metabolism is not accessible in yeast models.

The green alga Chlamydomonas reinhardtii Dangeard is a well-known model of a photosynthetic cell. This unicellular eukaryote is widely used for studies of a number of physiological processes, such as photosynthesis, respiration, nitrogen assimilation, flagella motility, and basal body function (Rochaix et al., 1998; Harris, 2001; Silflow and Lefebvre, 2001). Its short life cycle is easily controlled, making this organism a powerful tool for genetic analysis. The past few years have seen a dramatic development of the molecular technologies applicable to C. reinhardtii (Lefebvre and Silflow, 1999; Grossman, 2000; Fuhrmann, 2002). In the near future, the completion of the Chlamydomonas genome project, which includes molecular mapping, whole-genome sequencing, and extended expressed sequence tag (EST) generation programs (Grossman et al., 2003; Kathir et al., 2003; Shrager et al., 2003), will make C. reinhardtii an even more accessible model.

C. reinhardtii has also been used in the past to study metal tolerance (for review, see Hanikenne, 2003) and homeostasis, and especially iron and copper metabolism in the chloroplast (e.g. Merchant, 1998; Moseley et al., 2002a, 2002b). In a recent review, we recommended the use of C. reinhardtii as a photosynthetic eukaryotic model to study metal homeostasis and tolerance (Hanikenne, 2003).

The red alga Cyanidioschyzon merolae De Luca, Taddei, and Varano is a blue-green-colored unicellular alga that is found in sulfur-rich acidic hot springs. Its ultrastructure and other cytological and biochemical features suggest that it is one of the most primitive algae. C. merolae has no rigid cell wall and contains one nucleus, one mitochondrion, one chloroplast with a centrally located nucleoid, one Golgi body, and one microbody (Seckbach, 1994; Seckbach and Ott, 1994; Matsuzaki et al., 2004). Belonging to the extremophilic Cyanidiales, it shares the same photosynthetic pigments and ecological niche as Cyanidium caldarium and Galdieria sulphuraria, but differs by being smaller, structurally simpler, and dividing only by binary fission (no formation of endospores; Ott and Seckbach, 1994; Gross, 2000; Ciniglia et al., 2004). C. merolae was extensively used in studies of mitochondrial and chloroplast divisions, which can be highly synchronized by light/dark cycles (Kuroiwa, 1998). Its nuclear genome is one of the smallest known in photosynthetic eukaryotes. It has just been completely characterized (Matsuzaki et al., 2004), with the mitochondrial (Ohta et al., 1998) and plastid genomes (Ohta et al., 2003) also being available.

According to many phylogenetic studies, red algae and green plants (including green algae) are considered to be sister groups (Cavalier-Smith, 1982; Bhattacharya and Medlin, 1995; Moreira et al., 2000). In particular, red and green plastids are thought to be derived from the same primary endosymbiotic event, the former retaining more cyanobacterial features than the latter (Nozaki et al., 2003; Ohta et al., 2003). Along with fossil evidence (Butterfield, 2000), this supports the idea that red algae are among the most primitive eukaryotes still living today. Red algae are also worth being studied because their plastid seems to be at the origin of the plastids of various algal groups of fundamental importance in marine primary production (diatoms, dinophytes, and brown algae) through a process of secondary endosymbiosis (Moreira and Philippe, 2001; Yoon et al., 2002).

	A MINING STRATEGY SUITED FOR MULTIGENE FAMILIES

TOP TWO UNICELLULAR ALGAL MODEL... A MINING STRATEGY SUITED... THE CDF FAMILY THE ZIP FAMILY THE CAX FAMILY THE COPT FAMILY THE HMA SUBFAMILY OF... THE ABC TRANSPORTER FAMILY IRON TRANSPORTERS METHODOLOGICAL INSIGHTS PHYLOGENETIC INSIGHTS CONCLUSION LITERATURE CITED

 
We used a high-throughput semiautomated approach to identify metal transporters in the recently released genomic sequences of C. reinhardtii (U.S. Department of Energy Joint Genome Institute [JGI; http://www.jgi.doe.gov]) and C. merolae (Matsuzaki et al., 2004). Briefly, numerous BLAST searches were conducted on both algal genomes, using all members of the human, yeast, and Arabidopsis CDF, ZIP, CAX, COPT, HMA, yellow-stripe 1-like (YSL), Fe transporter (FTR), NRAMP, and iron-regulated 1 (IREG1)-like (sub)families, as well as multidrug resistance-associated protein (MRP) and ABC transporter of the mitochondria (ATM)/heavy-metal tolerance (HMT) subfamilies of ABC transporters (see Supplemental Table I). This custom strategy allowed us to identify, respectively, 41 and 25 putative metal transporters in C. reinhardtii and C. merolae. Their distribution among the different transporter families, compared to the three reference organisms, is given in Table I. The identified proteins were first subjected to topology and targeting predictions, then used to mine available EST databases (Tables II and III; Supplemental Tables II–VII). Finally, evolutionary relationships within each transporter family were examined using robust phylogenetic trees (Figs. 1–7 ). A detailed description of the whole procedure is given in Supplemental Figure 1.

View larger version (16K): [in this window] [in a new window]   Figure 1. Phylogenetic tree of the CDF family. The tree is the strict consensus of the five most parsimonious trees (8,358 steps) for which branch lengths were estimated by maximum likelihood (ML) using the VT model of amino acid substitution (log L = –37,883.42). Quartet puzzling support values above 50% from an independent ML reconstruction (10,000 puzzling steps) are printed in bold type and bootstrap percentages above 50% from an independent distance reconstruction (1,000 replicates) are printed in italic type (see Supplemental Fig. 1 for methodological details). Protein accession numbers from organisms other than C. reinhardtii and C. merolae are given in Supplemental Table I, except for AhMTP1 (CAD89013, TgMTP1 (AAS67026, PtdMTP1 (AAR23528, ShMTP1 (AAO38707, and SpZHF1 (NP_593645). Ah, Arabidopsis halleri; Ptd, Populus balsamifera subsp. trichocarpa x Populus deltoids; Sh, Stylosanthes hamata; Sp, Schizosaccharomyces pombe; Tg, Thlaspi goesingense.

View larger version (26K): [in this window] [in a new window]   Figure 2. Phylogenetic tree of the ZIP family. The tree is the most parsimonious tree (10,168 steps) for which branch lengths were estimated by maximum likelihood (ML) using the WAG model of amino acid substitution (log L = –43,192.29). Quartet puzzling support values above 50% from an independent ML reconstruction (25,000 puzzling steps) are printed in bold type and bootstrap percentages above 50% from an independent distance reconstruction (1,000 replicates) are printed in italic type (see Supplemental Fig. 1 for methodological details). CrZIP13 and CrZIP14 positions were ignored in previous steps because they were subject to long-branch attraction artifacts in both the distance and ML reconstructions. The inset is enlarged twice. Protein accession numbers from organisms other than C. reinhardtii and C. merolae are given in Supplemental Table I, except for TcZNT1 (AAF61374, TcZNT2 (AAK69429, and MxGUFA (CAA50380. Mx, Myxococcus xanthus; Tc, Thlaspi caerulescens.

View larger version (8K): [in this window] [in a new window]   Figure 3. Phylogenetic tree of the CAX family. The tree is the most parsimonious tree (4,207 steps) for which branch lengths were estimated by maximum likelihood (ML) using the VT model of amino acid substitution (log L = –19,012.92). Quartet puzzling support values above 50% from an independent ML reconstruction (1,000 puzzling steps) are printed in bold type and bootstrap percentages above 50% from an independent distance reconstruction (1,000 replicates) are printed in italic type (see Supplemental Fig. 1 for methodological details). Protein accession numbers from organisms other than C. reinhardtii and C. merolae are given in Supplemental Table I.

View larger version (14K): [in this window] [in a new window]   Figure 4. Phylogenetic tree of the HMA subfamily of P-type ATPases. The tree is the most parsimonious tree (9,500 steps) for which branch lengths were estimated by maximum likelihood (ML) using the WAG model of amino acid substitution (log L = –44,546.19). Quartet puzzling support values above 50% from an independent ML reconstruction (1,000 puzzling steps) are printed in bold type and bootstrap percentages above 50% from an independent distance reconstruction (1,000 replicates) are printed in italic type (see Supplemental Fig. 1 for methodological details). Proteins having the uncharacteristic Ser/Pro/Cys motif in the sixth transmembrane domain are marked by an asterisk. Protein accession numbers from organisms other than C. reinhardtii and C. merolae are given in Supplemental Table I, except for SaCadA (BAC57487, EcZntA (NP_417926), and SynCoaT (S77012). Ec, Escherichia coli; Sa, Staphylococcus aureus; Syn, Synechocystis sp.

View larger version (18K): [in this window] [in a new window]   Figure 5. Phylogenetic tree of the MRP subfamily of ABC transporters. The tree is the most parsimonious tree (25,015 steps) for which branch lengths were estimated by maximum likelihood (ML) using the VT model of amino acid substitution (log L = –113,214.98). Quartet puzzling support values above 50% from an independent ML reconstruction (10,000 puzzling steps) are printed in bold type and bootstrap percentages above 50% from an independent distance reconstruction (1,000 replicates) are printed in italic type (see Supplemental Fig. 1 for methodological details). Protein accession numbers from organisms other than C. reinhardtii and C. merolae are given in Supplemental Table I. CrMRP1 (AAL35383 has been named previously by Im and Grossman (2002).

View larger version (22K): [in this window] [in a new window]   Figure 6. Phylogenetic tree of the ATM/HMT subfamily of ABC transporters. The tree is the most parsimonious tree (3,950 steps) for which branch lengths were estimated by maximum likelihood (ML) using the VT model of amino acid substitution (log L = –20,342.67). Quartet puzzling support values above 50% from an independent ML reconstruction (1,000 puzzling steps) are printed in bold type and bootstrap percentages above 50% from an independent distance reconstruction (1,000 replicates) are printed in italic type (see Supplemental Fig. 1 for methodological details). Protein accession numbers from organisms other than C. reinhardtii and C. merolae are given in Supplemental Table I. CrCds1 (AAQ19847 has been named previously by Hanikenne et al. (2005).

View larger version (12K): [in this window] [in a new window]   Figure 7. Phylogenetic tree of the NRAMP family. The tree is the strict consensus of the two most parsimonious trees (4,055 steps) for which branch lengths were estimated by maximum likelihood (ML) using the VT model of amino acid substitution (log L = –20,415.07). Quartet puzzling support values above 50% from an independent ML reconstruction (1,000 puzzling steps) are printed in bold type and bootstrap percentages above 50% from an independent distance reconstruction (1,000 replicates) are printed in italic type (see Supplemental Fig. 1 for methodological details). Protein accession numbers from organisms other than C. reinhardtii and C. merolae are given in Supplemental Table I, except for NpNRAMP1 (ZP_00111957), EcNRAMP1 (NP_416893), StNRAMP1 (NP_461349), and SaNRAMP1 (YP_040492). Ec, Escherichia coli; Np, Nostoc punctiforme; Sa, Staphylococcus aureus; St, Salmonella typhymurium.

	THE CDF FAMILY

TOP TWO UNICELLULAR ALGAL MODEL... A MINING STRATEGY SUITED... THE CDF FAMILY THE ZIP FAMILY THE CAX FAMILY THE COPT FAMILY THE HMA SUBFAMILY OF... THE ABC TRANSPORTER FAMILY IRON TRANSPORTERS METHODOLOGICAL INSIGHTS PHYLOGENETIC INSIGHTS CONCLUSION LITERATURE CITED

We have identified, respectively, five and three MTPs in the genome sequences of C. reinhardtii and C. merolae (Tables II and III). CrMTP1 and CmMTP1 are related to zinc transporters of higher plants, humans, and, more distantly, yeasts (Fig. 1 ). AtMTP1 (formerly ZAT) confers zinc resistance to Arabidopsis when overexpressed (van der Zaal et al., 1999). Homologs of AtMTP1 have been identified in several species, including poplar and the metal hyperaccumulators Arabidopsis halleri and Thlaspi goesingense, and were shown to be involved in zinc homeostasis (Persans et al., 2001; Blaudez et al., 2003; Becher et al., 2004; Drager et al., 2004; Kim et al., 2004). In A. halleri, two AhMTP1 genes, which are particularly highly expressed compared to Arabidopsis, cosegregate with zinc tolerance. The AhMTP1 protein is localized in the vacuolar membrane and likely to mediate the accumulation of zinc in the vacuole (Becher et al., 2004; Drager et al., 2004). In S. cerevisiae, the ScZRC1 and ScCOT1 proteins are located in the vacuolar membrane and export zinc, cobalt, and possibly cadmium from the cytoplasm into the vacuole, playing an essential role in the tolerance to these metals (Gaither and Eide, 2001). In S. pombe, SpZHF1 is an endoplasmic reticulum (ER)-localized protein and is responsible for zinc storage in this compartment. The disruption of the gene results in zinc and cobalt sensitivity but tolerance to cadmium and nickel. SpZHF1 is also required for growth under zinc deficiency (Borrelly et al., 2002; Clemens et al., 2002). Nine MTPs (HsZnT1–9) are found in humans, with HsZnT1 to 7 being responsible for zinc transport in different cell types and organs (for review, see Palmiter and Huang, 2004). Both CrMTP1 and CmMTP1 are predicted to localize to the vacuolar membrane (Tables II and III), similar to their plant and yeast homologs. Interestingly, the three ESTs corresponding to CrMTP1 were only identified in the stress II cDNA library (Supplemental Table IV) prepared from cells exposed to cadmium (Shrager et al., 2003), therefore suggesting that CrMTP1 might be involved cadmium detoxification.

CrMTP2 to 4 and CmMTP2 cluster with AtMTP8 to 11 and with ShMTP1, a subgroup of putative manganese transporters. The ShMTP1 protein of the manganese-tolerant legume Stylosanthes hamata confers manganese resistance to yeast and Arabidopsis when overexpressed. It is localized in the plant vacuolar membrane and was proposed to be involved in manganese sequestration in this organelle (Delhaize et al., 2003). By analogy, AtMTP8 to 11 may have similar substrate specificity, although no experimental data are available. Our phylogenetic analysis suggests that the diversification of this subfamily in C. reinhardtii and Arabidopsis evolved independently through gene duplications. In that respect, it is interesting to note that CrMTP2 to 4 are found on the same genome scaffold very close to each other (Supplemental Table II). Given the important need of manganese to maintain photosynthesis, it is possible that at least one of these proteins is located in the chloroplast, although our prediction analysis only localized these proteins in the vacuole or the mitochondria (Supplemental Tables VI and VII).

CrMTP5 groups with AtMTP7 and HsZnT9, two proteins of unknown function, whereas the remaining CmMTP3 is distantly related to a few other CDF proteins. It is, however, worth mentioning that, similar to ScMSC2, AtMTP12, and HsZnT5 (Fig. 1), CmMTP3 possesses a long N-terminal extension with 6 predicted additional transmembrane domains that share no homology with other CDFs. HsZnT5 is believed to act in the sequestration of zinc in the Golgi complex of different tissues (Palmiter and Huang, 2004), and ScMSC2 is localized in the ER and is essential to maintain proper zinc homeostasis in this compartment (Ellis et al., 2004). Our analysis suggests that CmMTP3 is targeted to the secretory pathway as well (Table III).

Finally, neither C. reinhardtii nor C. merolae seem to possess a homolog of the yeast ScMMT1 and ScMMT2 (Fig. 1), two mitochondrial proteins believed to play a role in iron homeostasis (Li and Kaplan, 1997).

	THE ZIP FAMILY

TOP TWO UNICELLULAR ALGAL MODEL... A MINING STRATEGY SUITED... THE CDF FAMILY THE ZIP FAMILY THE CAX FAMILY THE COPT FAMILY THE HMA SUBFAMILY OF... THE ABC TRANSPORTER FAMILY IRON TRANSPORTERS METHODOLOGICAL INSIGHTS PHYLOGENETIC INSIGHTS CONCLUSION LITERATURE CITED

 
The ZIP protein family forms another ubiquitous transporter family involved in metal homeostasis, generally mediating the influx of metal cations, like zinc, iron, cadmium, or manganese, from outside the cell or from a subcellular compartment into the cytoplasm. ZIPs are predicted to have eight transmembrane domains, with the N and C termini being extracytoplasmic. As a common feature, ZIPs possess a long cytoplasmic loop (the so-called variable region) of variable length and sequence between transmembrane domains III and IV. The variable region very often contains a probable metal-binding His-rich domain. Transmembrane domains IV and V are amphipathic and believed to form a polar cavity required for the cation metal transport, while the loop between transmembrane domains II and III could be the site of initial binding of the substrate (Guerinot and Eide, 1999; Guerinot, 2000; Gaither and Eide, 2001).

Based on sequence conservation, Gaither and Eide (2001) have classified the ZIP proteins into four subfamilies. Subfamily I includes ZIPs from plants and fungi, while subfamily II contains nematode and mammalian proteins. The GUFA subfamily contains proteins of mainly unknown function present in both prokaryotes and eukaryotes, and the LIV1 group is only found in eukaryotes.

We have identified, respectively, 14 and 4 ZIPs in the genomes of C. reinhardtii and C. merolae and, altogether, our findings suggest that these algal ZIPs are mainly localized in the vacuolar or the plasma membrane (Tables II and III).

CrZIP1 to 5 and CmZIP1 cluster with subfamily I, which contains most of the Arabidopsis and S. cerevisae ZIPs (Fig. 2 ). Interestingly, the proteins of each of the four species form distinct groups within the subfamily. The protein diversification probably occurred independently in the different groups through duplications of an ancestral gene. In that respect, it is noteworthy that CrZIP4 and 5 are found in tandem on the same genome scaffold (Supplemental Table II). The algal proteins seem to branch shortly before the ScZRT1 and ScZRT2 proteins, the yeast high- and low-affinity zinc uptake systems, respectively, located in the plasma membrane (Zhao and Eide, 1996a, 1996b).

A few Arabidopsis subfamily I ZIPs have been characterized through yeast complementation and expression analyses. AtZIP1 to 4 play a role in cellular zinc uptake, AtZIP1, 3, and 4 being induced under zinc-limiting conditions at the transcriptional level (Guerinot, 2000; Gaither and Eide, 2001; Hall and Williams, 2003). In the zinc/cadmium hyperaccumulator Thlaspi caerulescens, TcZNT1 and TcZNT2, two ZIP transporters, are probably involved in zinc hyperaccumulation (Lasat et al., 2000; Pence et al., 2000; Assuncao et al., 2001), while in the zinc hyperaccumulator A. halleri, AhZIP6 and AhZIP9 are constitutively highly expressed compared to Arabidopsis, respectively, in shoots and in roots (Becher et al., 2004; Weber et al., 2004). In Arabidopsis, the AtIRT1 and AtIRT2 proteins are expressed in roots under iron-limiting conditions. In addition to iron, AtIRT1 also transports zinc, manganese, and cadmium, whereas AtIRT2 is only responsible for zinc uptake (Guerinot, 2000; Gaither and Eide, 2001; Hall and Williams, 2003). In Arabidopsis, AtIRT1 represents the major iron uptake system in roots. None of the CrZIPs identified in this article appears to be clearly related to this plant iron uptake system. Our findings are in agreement with the results of Herbik et al. (2002) and La Fontaine et al. (2002), revealing the occurrence in C. reinhardtii of an iron assimilation pathway related to the high-affinity iron uptake pathway of S. cerevisiae (see also below).

CrZIP6 belongs to ZIP subfamily II, together with CmZIP2 and HsZIP1 to 3 (Fig. 2). While HsZIP3 has not been characterized, HsZIP1 and 2 are involved in zinc uptake across the plasma membrane (Eide, 2004). Surprisingly, C. reinhardtii possesses six GUFA-related proteins (CrZIP7–12), whereas only one is found in C. merolae (CmZIP3/GUFA), yeast (ScZRT3), Arabidopsis (AtGUFA), and humans (HsZIP11; Fig. 2). The distant relationship of CrZIP7 and 12 with the GUFA subfamily is only supported by neighbor-joining analysis (data not shown). The GUFA subfamily has been named after the GUFA protein of the bacteria Myxococcus xanthus, a protein of unknown function identified by genome sequencing (McGowan et al., 1993). The yeast protein ScZRT3 is the only GUFA protein that has been functionally characterized. It is localized in the vacuolar membrane and involved in the remobilization of stored zinc under zinc-deficiency conditions (MacDiarmid et al., 2000).

Finally, a LIV1-like protein is found in C. merolae but not in C. reinhardtii. CmZIP4/LIV1 contains a highly conserved motif [(H,E) E(L,F) P(H,Q,A) E(L,I,V,M)(G,S) D(F,L,V)(M,A,V,G)XL(L,I,V), defined by Taylor and Nicholson (2003) as a signature for human LIV1 proteins] in the transmembrane domain V, which is specific to this subfamily of ZIPs. Initially implicated in human metastatic breast cancer, the LIV1 proteins have recently been shown to be capable of zinc transport (Taylor et al., 2003, 2004).

It should be mentioned that the CrZIP13 and 14 proteins are only distantly related to the ZIP family. Due to difficulties in aligning them properly with other ZIPs, their positions in the tree vary according to the phylogenetic method used (data not shown). These proteins may be based on inaccurate gene models and should be considered as possibly distant, putative ZIPs.

Nonetheless, our findings surprisingly reveal the presence of a very high number of putative zinc transporters of the ZIP family in the C. reinhardtii genome compared to yeast and C. merolae. The unicellular C. reinhardtii possesses as many ZIPs as the complex multicellular eukaryotes Arabidopsis and humans, although the ZIP diversification occurred in different subfamilies in each of these organisms (within subfamily I for Arabidopsis, LIV1 for humans, and GUFA and subfamily I for C. reinhardtii; Fig. 2) from possibly four ancestral genes still present in C. merolae.

	THE CAX FAMILY

TOP TWO UNICELLULAR ALGAL MODEL... A MINING STRATEGY SUITED... THE CDF FAMILY THE ZIP FAMILY THE CAX FAMILY THE COPT FAMILY THE HMA SUBFAMILY OF... THE ABC TRANSPORTER FAMILY IRON TRANSPORTERS METHODOLOGICAL INSIGHTS PHYLOGENETIC INSIGHTS CONCLUSION LITERATURE CITED

 
The CAX proteins are divalent cation/H⁺ antiporters generally containing 10 to 14 transmembrane domains. The AtCAX1 and AtCAX2 proteins have been identified by functional complementation of a S. cerevisiae mutant defective in vacuolar calcium accumulation. Both having 11 putative transmembrane domains, AtCAX1 is a vacuolar high-affinity Ca²⁺/H⁺ antiporter, while AtCAX2 displays a lower affinity for Ca²⁺ (Hirschi et al., 1996) and has been proposed to transport Mn²⁺ and Cd²⁺ across the tonoplast of plant cells (Hirschi et al., 2000). Nine other CAX proteins have been identified in the Arabidopsis genome, but have not been functionally characterized yet (Maser et al., 2001). An additional member of this family is the AtMHX1 protein, which is an H⁺-coupled antiporter that can transport Mg²⁺ and Zn²⁺ across the vacuole membrane (Shaul et al., 1999).

We have found, respectively, four and two CAXs in the genomes of C. reinhardtii and C. merolae (Tables II and III). CrCAX1 to 2 and CmCAX1 are related to the yeast ScVCX1 and AtCAX1 to 6 (Fig. 3 ). CrCAX3 to 4 and CmCAX2 are spread among the other CAXs. Maser et al. (2001) highlighted the difficulty of determining transport properties from sequence information within this family. Experimental data will be necessary to establish whether any of the newly identified CAXs plays a role in transition metal homeostasis in algae. To support this idea, it is worth mentioning that CrCAX1 is expressed in the stress II and stress III EST libraries (Supplemental Table IV), both of which were prepared from metal-stressed cells (Shrager et al., 2003).

Finally, it is interesting to note that neither C. reinhardtii nor C. merolae seem to possess a close AtMHX1 homolog, suggesting that this function is lacking in algae. In Arabidopsis, AtMHX1 is mainly expressed in vascular tissues and believed to participate in Mg²⁺ and Zn²⁺ partitioning between plant organs (Shaul et al., 1999). This function would have evolved associated with the emergence of land plants. Alternatively, it might be performed by a more distant member of the CAX family.

	THE COPT FAMILY

TOP TWO UNICELLULAR ALGAL MODEL... A MINING STRATEGY SUITED... THE CDF FAMILY THE ZIP FAMILY THE CAX FAMILY THE COPT FAMILY THE HMA SUBFAMILY OF... THE ABC TRANSPORTER FAMILY IRON TRANSPORTERS METHODOLOGICAL INSIGHTS PHYLOGENETIC INSIGHTS CONCLUSION LITERATURE CITED

 
The COPT proteins form a eukaryotic family of copper transporters (Eide, 1998). COPTs share common topology features, including (1) the presence of three transmembrane domains with an extracellular N terminus and a cytoplasmic C terminus; (2) the presence of two Met-rich regions in the N terminus that may act as a copper scavenger in the extracellular milieu; and (3) an additional Met-rich region in transmembrane domain II, probably involved in copper coordination during transmembrane transport (Petris, 2004). In S. cerevisiae, high-affinity copper uptake, like iron uptake (see below), requires plasma membrane reductases to reduce Cu(II) to Cu(I), the reduction step being mediated by the same reductases that are involved in iron uptake. Cu(I) is then transported by the redundant high-affinity copper transporters ScCTR1 and ScCTR3, via a Cu(I)/2 K(I) antiport mechanism that is highly specific for Cu(I) over other metal ions (Eide, 1998). S. cerevisiae possesses an additional CTR-related protein (ScCTR2). Located in the vacuole, the function of ScCTR2 remains unknown (Eide, 1998; Petris, 2004). Two CTR homologs are found in humans, but only HsCTR1 has been functionally characterized. It is a plasma membrane protein involved in copper uptake in various cell types (Zhou and Gitschier, 1997; Petris, 2004). Sancenon et al. (2003) have recently identified 5 COPT proteins (AtCOPT1 to 5) in Arabidopsis and showed that AtCOPT1, 2, 3, and 5 restore copper uptake in a ctr1 ctr3 yeast mutant. Antisense AtCOPT1 transgenic lines display reduced copper uptake together with increased root length and defects in pollen development (Sancenon et al., 2004). AtCOPT1 thus plays an important role in copper acquisition in Arabidopsis.

We have identified a single COPT protein in the genomes of both C. reinhardtii and C. merolae (Tables II and III). CrCOPT1 is more related to the Arabidopsis COPTs than to the human and yeast transporters, while CmCOPT1 is only distantly related to the other COPTs (data not shown). It should also be mentioned that CmCOPT1 is substantially larger than the plant and human COPTs. CrCOPT1 might be located in the plasma membrane (Supplemental Table VI), while CmCOPT1 is predicted to be localized in the vacuolar membrane (Supplemental Table VII).

	THE HMA SUBFAMILY OF P-TYPE ATPASES

TOP TWO UNICELLULAR ALGAL MODEL... A MINING STRATEGY SUITED... THE CDF FAMILY THE ZIP FAMILY THE CAX FAMILY THE COPT FAMILY THE HMA SUBFAMILY OF... THE ABC TRANSPORTER FAMILY IRON TRANSPORTERS METHODOLOGICAL INSIGHTS PHYLOGENETIC INSIGHTS CONCLUSION LITERATURE CITED

 
P-type ATPases are transporters characterized by the formation of a phosphorylated intermediate in the reaction cycle. These proteins typically contain 8 to 12 transmembrane domains and a large cytoplasmic loop, including ATP-binding and phosphorylation sites. P-type ATPases transport a broad range of small cations, and possibly phospholipids, and have been classified into five subfamilies according to their predicted substrate specificities and phylogenies. The type 1_B subfamily proteins (HMAs or CPx ATPases) are involved in heavy-metal transport. HMAs possess eight transmembrane domains, the sixth of which contains a conserved Cys-Pro-Cys/His/Ser motif (CPx motif) believed to be involved in metal cation translocation across the membrane. Identified in a wide range of organisms, HMAs can be divided into two main groups with different substrate specificities (monovalent Cu⁺/Ag⁺ cations or divalent Zn²⁺/Co²⁺/Cd²⁺/Pb²⁺ cations; Axelsen and Palmgren, 2001; Cobbett et al., 2003).

We have identified three and two HMAs in the genomes of C. reinhardtii and C. merolae, respectively (Tables II and III). CrHMA1 and CmHMA1 are related to the divalent cation transporters of bacteria (EcZntA, SaCadA) and Arabidopsis (AtHMA1 to 4), especially to AtHMA1 (Fig. 4 ). EcZntA and SaCadA are involved in zinc or cadmium excretion and tolerance in Escherichia coli and Staphylococcus aureus, respectively (Nies, 2003). AtHMA2 and AtHMA4 probably play a role in zinc translocation in the plant (Mills et al., 2003; Hussain et al., 2004), while functional studies in yeast suggest that AtHMA3 is a cadmium/lead transporter (Gravot et al., 2004). In addition, AhHMA3 and TcHMA4 are constitutively highly expressed in the zinc hyperaccumulator A. halleri and the zinc/cadmium hyperaccumulator T. caerulescens, respectively, when compared to Arabidopsis. These proteins could be at least partly responsible for the metal hyperaccumulation and hypertolerance in these two species (Becher et al., 2004; Bernard et al., 2004). So far, nothing is known about the function of AtHMA1. According to phylogenetic analysis, it belongs to the divalent cation transporter subfamily of P_1B-ATPases but is only distantly related to AtHMA2 to 4 (Cobbett et al., 2003). Interestingly, like AtHMA1, CrHMA1 and CmHMA1 have an uncharacteristic Ser/Pro/Cys motif in the sixth transmembrane domain instead of the common Cys-Pro-Cys/His/Ser motif. This property is shared with the cobalt transporter CoaT of the cyanobacterium Synechocystis and might influence the substrate specificity of the transporter. On that basis, Cobbett et al. (2003) suggested that AtHMA1 might be a cobalt transporter.

Although widespread in prokaryotes, the occurrence of these divalent cation-transporting HMAs in eukaryotes is apparently limited to photosynthetic organisms (Fig. 4). Cobbett et al. (2003) suggested that these proteins might have evolved from one or more, probably two, horizontal gene transfers from the prokaryotic endosymbiont believed to be at the origin of the plastid. The AtHMA2 to 4 genes apparently result from the recent repeated duplication of an ancestral gene (Cobbett et al., 2003). The AtHMA2 to 4 subgroup is absent in C. merolae and C. reinhardtii, indicating that this function might have been lost in both algae.

CrHMA2 and 3 and CmHMA2 cluster with the monovalent cation transporters of the P_1B-ATPases and are most probably copper transporters (Fig. 4). In particular, CrHMA2 and CmHMA2 are closely related to yeast (ScCCC2), plant (AtRAN1), and human (HsATP7A and B) transporters known to deliver copper to proteins in the trans-Golgi network (Cobbett et al., 2003; Hall and Williams, 2003). The CrHMA2 protein corresponds to the EST 1021007C06, which has been previously identified, and is believed to be involved, together with the chaperone ATX1, in copper delivery to the secretory pathway (La Fontaine et al., 2002).

CrHMA3 clusters with AtPAA1 and AtHMA6 (Fig. 4). AtPAA1 is a chloroplastic protein mediating the transport of copper into the chloroplast. Arabidopsis paa mutants have a lower chloroplast copper content, accumulate reduced levels of holoplastocyanin, and display reduced chloroplastic copper/zinc superoxide dismutase activity (Shikanai et al., 2003). By analogy, CrHMA3 could play a role in copper delivery to the chloroplast. However, this protein is predicted to localize to the vacuolar membrane (Table II). In some cyanobacteria and green algae, plastocyanin (involved in the transfer of photosynthetic electrons from cytochrome b₆/f to PSI) may be functionally replaced by a soluble c-type (c-552 or c-553) cytochrome. Plastocyanin is expressed under copper-replete conditions and replaced by cytochrome c₆ when copper availability is limited (Kerfeld and Krogmann, 1998; Merchant, 1998). Whereas higher plants only express plastocyanin, all examined red algae (along with Chrysophyceae and Xanthophyceae) seem to use the c-type cytochrome alternative (Sandmann et al., 1983; Price et al., 1991). Accordingly, C. merolae also lacks plastocyanin (data not shown) and its plastid genome contains a petJ gene encoding a c-553 cytochrome (Ohta et al., 2003). The fact that photosynthesis is not based on plastocyanin in C. merolae may explain the lack of an AtPAA1-like function in this organism (Fig. 4).

	THE ABC TRANSPORTER FAMILY

TOP TWO UNICELLULAR ALGAL MODEL... A MINING STRATEGY SUITED... THE CDF FAMILY THE ZIP FAMILY THE CAX FAMILY THE COPT FAMILY THE HMA SUBFAMILY OF... THE ABC TRANSPORTER FAMILY IRON TRANSPORTERS METHODOLOGICAL INSIGHTS PHYLOGENETIC INSIGHTS CONCLUSION LITERATURE CITED

 
ABC transporters are ubiquitous transporters involved in a large number of physiological processes. This family is one of the largest protein families with 29, 128, and 48 members in S. cerevisiae, Arabidopsis, and humans, respectively. Typical ABC transporters (the so-called full-size transporters) possess two conserved nucleotide-binding folds responsible for ATP hydrolysis alternating with two highly hydrophobic domains (containing 4–6 transmembrane spans) that specify the substrates to be transported. The half-size ABC transporters possess a single copy of each domain and are assumed to function as homo- or heterodimers (Holland et al., 2003).

Based on structural similarities, ABC transporters can be classified in several subfamilies (Decottignies and Goffeau, 1997; Sanchez-Fernandez et al., 2001; Dean et al., 2003). Among them, members of only two subfamilies (MRP and ATM/HMT) have been implicated in metal transport (Table I).

MRPs are full-size ABC transporters mainly acting as glutathione S-conjugate pumps (Rea et al., 1998). In S. cerevisiae, the ScYCF1 protein ensures the transport of bis(glutathionato) cadmium complexes from the cytoplasm into the vacuole. The lack of the transporter determines hypersensitivity to cadmium, arsenate, and mercury (Szczypka et al., 1994; Li et al., 1997; Ghosh et al., 1999; Gueldry et al., 2003). A second vacuolar yeast MRP, named ScBTP1, is also involved in cadmium tolerance, but its role is only marginal compared to ScYCF1 (Klein et al., 2002; Sharma et al., 2002). The overexpression of ScYCF1 in Arabidopsis enhances the tolerance of the plant to lead and cadmium (Song et al., 2003). Among 14 MRP proteins (AtMRP1–14) identified in Arabidopsis (Sanchez-Fernandez et al., 2001; Kolukisaoglu et al., 2002; Martinoia et al., 2002), only AtMRP3 has been implicated in cadmium detoxification and transport. The AtMRP3 gene expression is up-regulated in response to cadmium and this induction is apparently correlated with the accumulation of cadmium in the plant organs (Bovet et al., 2003). Moreover, AtMRP3 complements the cadmium sensitivity of a yeast ycf1 mutant (Tommasini et al., 1998). None of the 4 other cloned and characterized AtMRPs (AtMRP1, 2, 4, and 5) or the 12 HsMRPs (or HsABCC) are involved in metal detoxification, but in the transport of a wide range of substrates (for review, see Dean et al., 2003; Rea et al., 2003).

We have found seven and two MRPs in the genome sequences of C. reinhardtii and C. merolae, respectively (Tables II and III). Although it is obviously a MRP transporter, protein model 155613 is too small (428 amino acid residues) compared to the other members of the family. With a TBLASTN search, we have identified additional ABC transporter domains in close proximity to the gene model on scaffold 142 (Table II; Supplemental Table II). Using three gene prediction software programs (GreenGenie, GeneMark, and GENSCAN), we determined an alternative gene model (encoding CrMRP6, 1,524 amino acid residues), which, albeit imperfect, probably much better reflects the real structure of the gene (data not shown; available upon request). The cloning of the corresponding cDNA will nevertheless be necessary to fully elucidate the gene structure.

CrMRP1 and 2 group with human (HsABCC4, 5, 11, and 12), plant (AtMRP11 and 15), and yeast (ScYOR1) proteins (Fig. 5 ), all of which lack the N-terminal domain characteristic of the members of the family (Dean et al., 2003; Rea et al., 2003). CrMRP1 is a chloroplastic protein that might be involved in bicarbonate uptake and is regulated by light intensity and CO₂ level (Im and Grossman, 2002). The other members of the MRP subfamily identified here, CrMRP3 to 7 and CmMRP1 and 2, are related to AtMRP11 and HsABCC10 (Fig. 5). Based on sequence homology, these algal transporters are probably glutathione S-conjugate pumps. However, to understand the function of these proteins in drug and/or metal tolerance, a detailed analysis (including the determination of the membrane localization and substrate specificity) will be required.

ATM/HMTs are half-size transporters located either in the mitochondrial or vacuolar membranes. Mitochondrial transporters (HsABCB6 and 7, AtATM1–3, and ScATM1) are involved in the export of iron/sulfur clusters from the mitochondrial matrix to the cytoplasm (Kispal et al., 1997; Csere et al., 1998; Allikmets et al., 1999; Kispal et al., 1999; Mitsuhashi et al., 2000; Kushnir et al., 2001). We have recently shown that a C. reinhardtii homolog of these mitochondrial transporters, named CrCds1, plays an essential role in cadmium tolerance, possibly by exporting cadmium out of the mitochondria, thereby protecting the mitochondrial function from cadmium toxicity (Hanikenne et al., 2001, 2005). The ability of mitochondrial ABC transporters to act in cadmium detoxification could be a new property among green organisms (Hanikenne et al., 2005). Out of 7 ESTs corresponding to the Cds1 gene (Supplemental Table IV), 3 were identified in the stress II cDNA library prepared from cells incubated in the presence of cadmium (Shrager et al., 2003), suggesting that the gene might be regulated by the metal. This was indeed demonstrated by northern-blot analysis, showing that Cds1 is strongly induced in the presence of cadmium (Hanikenne et al., 2005). This highlights that the careful analysis of EST data might be used to generate testable hypotheses concerning gene regulation.

The vacuolar SpHMT1 of the fission yeast S. pombe is involved in the transport of cadmium-phytochelatin complexes from the cytoplasm into the vacuole. A mutant strain lacking this transporter is unable to accumulate high-M_r cadmium-phytochelatin complexes in the vacuoles and displays hypersensitivity to cadmium (Ortiz et al., 1992, 1995). Although the existence of such a transport function has been biochemically confirmed in plant vacuolar membranes (Vögeli-Lange and Wagner, 1990; Salt and Rauser, 1995), no functionally related SpHMT1 analog has been identified yet in the Arabidopsis genome. The Arabidopsis proteins most highly similar to SpHMT1 are the ATM transporters (Sanchez-Fernandez et al., 2001; Martinoia et al., 2002).

We have identified two and three ATM/HMTs in the genomes of C. reinhardtii and C. merolae, respectively (Tables II and III). Intriguingly, we could not find any protein model corresponding to CrCds1 in the JGI sequence data but have found by TBLASTN that it falls within scaffold 122 (Table II; Supplemental Table II).

Two subclusters can be distinguished within the ATM/HMT subfamily (Fig. 6 ). On the one hand, subcluster I, which contains CrATM/HMT-2, CmATM/HMT-1 and -2, only includes mitochondrial transporters (AtATM1–3, ScATM1, HsABCB7) possessing 5 or 6 conserved transmembrane-spanning regions. Since they share structural and targeting predictions with ScATM1 (Supplemental Tables VI and VII), the 3 algal proteins are likely functional homologs of the yeast transporter and might have a role in iron homeostasis of mitochondria. On the other hand, subcluster II, which contains CrCds1, CrATM/HMT-3, and CmATM/HMT-3, includes both mitochondrial (HsABCB6) and vacuolar (SpHMT1) transporters possessing 5 additional transmembrane-spanning segments at the N-terminal end of the protein. Phytochelatins are the main intracellular chelators for cadmium and accumulate in the vacuole in C. reinhardtii (Howe and Merchant, 1992; Hu et al., 2001). Therefore, a functional analog of SpHMT1 should be present in the alga and might possibly be fulfilled by CrATM/HMT-3. We will have to wait for experimental evidence supporting these hypotheses. Although a putative phytochelatin synthase gene (protein model CMI111C) is present in the genome of C. merolae, the occurrence of phytochelatins in this organism is not known. It is thus difficult to speculate on the function of CmATM/HMT-3.

	IRON TRANSPORTERS

TOP TWO UNICELLULAR ALGAL MODEL... A MINING STRATEGY SUITED... THE CDF FAMILY THE ZIP FAMILY THE CAX FAMILY THE COPT FAMILY THE HMA SUBFAMILY OF... THE ABC TRANSPORTER FAMILY IRON TRANSPORTERS METHODOLOGICAL INSIGHTS PHYLOGENETIC INSIGHTS CONCLUSION LITERATURE CITED

 
Two strategies are known for iron uptake in higher plants. Strategy I (or reduction strategy) occurs in all plants, except graminaceous monocots, and involves 3 steps: (1) soil acidification by H⁺ ATPases to solubilize iron; (2) reduction of ferric iron [Fe(III)] by plasma membrane ferric chelate reductases; and (3) uptake of ferrous iron [Fe(II)] by AtIRT1, a member of the ZIP family (Guerinot and Yi, 1994; Connolly and Guerinot, 2002). We have already mentioned above that neither C. reinhardtii nor C. merolae seem to possess any homologs of IRT1 (Fig. 2). Strategy II (or chelation strategy) is found in graminaceous monocots and involves the release of Fe(III) chelators called phytosiderophores, coupled to the induction of a specific Fe(III) chelate transporter (Guerinot and Yi, 1994; Connolly and Guerinot, 2002). In maize, this transporter, ZmYS1, has recently been cloned (Curie et al., 2001). ZmYS1 belongs to the oligopeptide transporter family and is induced under iron deficiency both at transcriptional and translational levels. It has been shown to cotransport protons and phytosiderophore- or nicotianamine-chelated iron. It is interesting to mention that, although ZmYS1 displays a broad metal specificity (iron, zinc, nickel, copper, and cobalt) in Xenopus oocytes and yeast, it seems to be regulated only by iron availability in planta (Roberts et al., 2004; Schaaf et al., 2004). Moreover, the rice OsYSL2 has been shown to transport iron and manganese chelated by nicotianamine (Koike et al., 2004). Phytosiderophores are not synthesized in Arabidopsis, a strategy I species, but eight YS1 homologs (AtYSL1–8) are found in the plant genome. The AtYSL2 protein has been shown to transport iron and copper chelated by nicotianamine, a precursor of phytosiderophores (DiDonato et al., 2004).

To our knowledge, phytosiderophores or nicotianamine are not found in algae. Nevertheless, we searched for YS1-like proteins in the C. reinhardtii and C. merolae genome sequence but could not identify any homolog. This function may thus have evolved after the emergence of land plants. As soon as genome sequence data become available, it would be interesting to analyze green algae that belong to the evolutionary lineage of higher plants (Charophyta, sensu lato; van den Hoek et al., 1995) in this respect.

In S. cerevisiae, the first step of iron uptake involves, as for the plant strategy I, the reduction of Fe(III) to Fe(II) by ferric reductases. Then occurs high-affinity iron uptake mediated by (1) a multicopper oxidase (ScFET3) that reoxidizes Fe(II) to Fe(III) and (2) an iron permease (ScFTR1) transporting Fe(III) into the cell. ScFET3 and ScFTR1 are induced under iron deficiency and the corresponding proteins form a complex at the plasma membrane. In iron-sufficient conditions, iron uptake is driven by the low-affinity transporter ScFET4 (Radisky and Kaplan, 1999).

Recently, La Fontaine et al. (2002) identified genes encoding a multicopper oxidase (CrFOX1) and an iron permease (CrFTR1) in the C. reinhardtii EST databases and showing, in a photosynthetic organism, the occurrence of an iron assimilation pathway analogous to the high-affinity iron uptake pathway of S. cerevisiae (for review, see Hanikenne, 2003). Although no additional FTR1 homolog was found in the C. reinhardtii genome, we have identified four FTR proteins (CmFTR1–4) in C. merolae (Tables II and III). CmFTR1 and 2 are almost identical in sequence and are highly represented in the EST databases (Table III; Supplemental Table V). Both are predicted to localize in the plasma membrane like ScFTR1 and CrFTR1 (Tables II and III). Interestingly, we could not find any multicopper oxidase homolog in C. merolae (data not shown). Living in an iron- and sulfur-rich milieu with low copper availability (see below), C. merolae might have an alternative Fe(II) oxidase associated with the FTRs and its iron uptake might be independent on copper, in contrast to S. cerevisiae. CmFTR3 and 4 are predicted to reside in the mitochondria (Supplemental Table VII). This putative function is not found in yeast and C. reinhardtii and might represent an independent evolution in C. merolae. Alternatively, such a function might have been lost in other studied unicellular eukaryotes.

Members of the NRAMP family are also known to transport Fe(II) ions. These transporters use the transmembrane proton gradient to facilitate transport of divalent cations and iron in particular. NRAMPs are ubiquitous proteins that possess common structural features, including the presence of 12 transmembrane domains, 2 conserved His residues in transmembrane domain 6, and a transport motif in the intracellular loop between transmembrane domains 8 and 9 (Forbes and Gros, 2001; Mackenzie and Hediger, 2004). In humans, HsNRAMP1 is required for natural resistance to intracellular pathogenic bacteria by transporting Mn(II) and Fe(II) across the phagosomal membranes in the macrophages, thus restricting metal availability for pathogens, whereas HsNRAMP2 (or HsDMT1) plays a role in the uptake of dietary Fe(II) at the apical surface of enterocytes (Forbes and Gros, 2001; Mackenzie and Hediger, 2004