Ammonium and Urea Transporter Inventory of the Selaginella and Physcomitrella Genomes

Ammonium and urea are important nitrogen sources for autotrophic organisms. Plant genomes encode several families of specific transporters for these molecules, plus other uptake mechanisms such as aquaporins and ABC transporters. Selaginella and Physcomitrella are representatives of lycophytes and bryophytes, respectively, and the recent completion of their genome sequences provided us with an opportunity for comparative genome studies, with special emphasis on the adaptive processes that accompanied the conquest of dry land and the evolution of a vascular system. Our phylogenetic analysis revealed that the number of genes encoding urea transporters underwent a progressive reduction during evolution, eventually down to a single copy in vascular plants. Conversely, no clear evolutionary pattern was found for ammonium transporters, and their number and distribution in families varies between species. In particular Selaginella, similar to rice, favors the AMT2/MEP family of ammonium transporters over the plant-specific AMT1 type. In comparison, Physcomitrella presents several members belonging to both families.


INTRODUCTION
Of all mineral elements required by plants, nitrogen is quantitatively the most important and thus often growth-limiting for many plants. Nitrogen is found in many organic compounds such as amino acids, nucleic acid, and consequently in proteins, as well as in nucleic acids. In most soils, nitrogen is heterogeneously distributed and found in various forms such as ammonium, nitrate, urea, amino acids, peptides, and even in water-insoluble forms (Jackson et al., 1990). Among the different nitrogen sources, ammonium is often preferred by plants, since its assimilation requires less energy. Moreover, ammonium is the predominant nitrogen form in anoxic soils, such as paddy rice fields, and it is abundantly used as fertilizer.
Due to the impermeability or poor permeability of the lipid bilayer of the plasma membrane to most nutrients including nitrogen compounds, uni and multicellular organisms developed a suite of nitrogen transporters. In 1979, electrophysiological studies in the unicellular algae Chara australis showed a positive inward current across the plasma membrane induced by ammonium addition (Walker et al., 1979). Fifteen years later, the first ammonium transporter, AtAMT1;1 from Arabidopsis, was cloned by yeast complementation (Ninnemann et al., 1994), in parallel with the yeast MEP orthologs (Marini et al., 1994). Since then, many para and orthologs were isolated either by yeast complementation or via homology cloning (Lauter et al., 1996;Gazzarrini et al., 1999;Simon-Rosin et al., 2003;Sonoda et al., 2003). Homologs have been found in bacteria, fungi, algae, plants, and animals (Marini et al., 1997(Marini et al., , 2000Gonzalez-Ballester et al., 2004) and the transport of the charged form of ammonium (NH + 4 ) was confirmed for the plant ammonium transporter AMT1;1 by electrophysiological characterization in Xenopus oocytes (Ludewig et al., 2002;Wood et al., 2006;Loqué et al., 2009).
In both Eubacteria and Archaea, AMT activity is regulated by the interaction with GlnK proteins (Coutts et al., 2002). Importantly, at least the bacterial AMTs, similar to the plant nitrate transporter CHL1 and the yeast glucose and amino acid transporters RGT2, SNF1, and GAP1 have a dual function as sensors and transporters, and thus can be defined as transceptors (Hyde et al., 2007;Ho et al., 2009;Thevelein and Voordeckers, 2009). Therefore, at least the bacterial AMTs link energy and nitrogen status of the cell with ammonium uptake to control ammoniumdependent transcription (Tremblay and Hallenbeck, 2009). Plants also contain GlnK homologs, however these are confined to the chloroplast (Hsieh et al., 1998). In plants, plasma membrane AMT activity is regulated by an allosteric auto-inhibition driven by the conserved cytosolic C-terminal tail, which connects the adjacent members in the trimeric complex (Loqué et al., 2007(Loqué et al., , 2009Neuhauser et al., 2007;Lanquar et al., 2009;Lanquar and Frommer, 2010). An increase in the concentration of extracellular ammonium leads to the phosphorylation of a specific threonine residue www.frontiersin.org in the cytosolic C-terminus of AtAMT1;1, resulting in closure of the pores of the transporter complex (Loqué et al., 2007;Lanquar et al., 2009). Since the C-terminus, and this regulatory threonine in particular, is highly conserved in plant AMT1s, it is conceivable that this feedback mechanism represents a common feature of plant ammonium transporters. We recently demonstrated that the Archaeglobus AMT1 also uses a trans-activation allosteric mechanism involving the conserved C-terminus, indicating that this unique regulatory system was developed early in evolution and has been maintained. This was surprising, because point mutations in the pore region were found that release the transporter from the strict dependence on trans-activation (Loqué et al., 2009). Such point mutations must have occurred during evolution, but apparently have not been fixed, suggesting a strong pressure to retain this regulatory mechanism, potentially to protect against accumulation of toxic levels of ammonium.
Urea is another nitrogen source used by plants that is commonly available in nature. Urea is mainly derived from decomposition of nitrogen-containing molecules and from urine secretion by animals. Additionally, urea is widely used in agriculture as a fertilizer due to its stability, low toxicity and cost-effective production.
Urea transporters are present in most organisms and belong to different classes, some functioning as passive channels, others as secondary active transporters. Interestingly, urea transporters commonly found in animals and bacteria (UT, UreI, Yut) are absent in plant and fungal genomes (Wang et al., 2008). Plant and fungal genomes encode a different family of urea transporters, called DUR3. In 1993 the first plant urea transporter (Dur3) was cloned from Arabidopsis by homology to the yeast transporter DUR3 (Elberry et al., 1993). DUR3 proteins are secondary active transporters, structurally belonging to the sodium solute symporter (SSS) superfamily (Wang et al., 2008). Electrophysiological studies showed that a proton, rather than sodium, is co-transported with urea by AtDUR3 (Liu et al., 2003). In planta characterization showed that DUR3 genes from Arabidopsis and rice are induced by nitrogen deficiency (Liu et al., 2003;Wang et al., 2011). DUR3s have an important role in high affinity urea transport under nitrogen limiting growth conditions (Kojima et al., 2007;Wang et al., 2011). DUR3 is not the only mechanism for urea uptake, since several aquaporins in the TIP and NOD26 families have been shown to facilitate urea import through the plasma membrane (for review, see Kojima et al., 2006). Additionally, ABC transporters are known to transport urea in cyanobacteria (Valladares et al., 2002). Typically, plant genomes contain > 100 ABC transporter homologs; whether any of these transport urea remains an open question.
The genome sequence of the lycophyte Selaginella moellendorffii (Banks et al., 2011) provides a valuable tool for comparative genomic studies in plants. The lineage of this species arose ∼400 millions years ago early after vascular plants had evolved, and shares many primitive traits with bryophytes. This species adds to the increasing list of plants for which genome sequences have become available, including the moss Physcomitrella patens (Rensing et al., 2008).
Here, we compared the ammonium and urea transporter homologs of two dicots (Arabidopsis and poplar), one monocot (rice) a lycophyte (Selaginella), a bryophyte (Physcomitrella), and a unicellular alga (Chlamydomonas reinhardtii) with the aim of gaining more insight into the evolution of these important transporters. For Arabidopsis, poplar, rice and Chlamydomonas, AMT and DUR3 members have already been annotated and characterized elsewhere, although in the case of poplar annotation is still not complete ( Table 1; Loqué and von Wirén, 2004;Cuturier et al., 2006;Fernandez and Galvan, 2007;Wang et al., 2008;Li et al., 2009). Our results indicate a progressive reduction of the number of genes encoding urea transporters during higher plant evolution, possibly suggesting a shift in nutrient form preference of angiosperms for other nitrogen sources. Analysis of number and distribution of ammonium transporters between the AMT1 and AMT2 families revealed that most species, with the exception of Physcomitrella, present a bias for one or the other subfamily (Arabidopsis: AMT1, Selaginella and rice: AMT2).

MATERIALS AND METHODS
Protein sequences of the six members of Arabidopsis AMTs and AtDUR3 were used to retrieve homologs of the ammonium and  (Table 1). For poplar, only the paralogs already annotated as AMT (Cuturier et al., 2006) were further considered in the phylogenic analysis. Sequences were blasted against a translated-nucleotide sequence database in the Ensembl Genome browser (www.ensemblgenomes.org). Searches for the urea transporter classes UT-A, -B, -C retrieved no hits in the plant databases. Within a species, coding and flanking sequences of the paralogs were carefully analyzed to distinguish between gene duplications and alleles. Protein sequences were aligned by ClustalW and poorly conserved N-and C-terminal domains (corresponding to amino acids 1-39 and 475-501 for AtAMT1;1 and 1-27 and 660-704 for AtDUR3) were trimmed (Figures A1 and A2 in Appendix). Phylogenetic trees were created on the basis of these alignments by the Mega5.05 program (Tamura et al., 2011) using a Neighbor-Joining, Poisson-corrected method, with pairwise deletion option and 1000 bootstraps. Additional trees were generated by Maximum-Likelihood method, Poisson-corrected, with 1000 bootstraps (Figures A3 and A4 in Appendix).

RESULTS AND DISCUSSION
Analysis of the S. moellendorffii genome indicates the presence of four ammonium transporters, with little variation between the two haplotypes sequenced (Tables 1 and 2). One AMT (SmAMT1) belongs to the AMT1 family and are closely related to the Physcomitrella AMTs (Figure 1; Figure A3 in Appendix).
The three other Selaginella AMTs fall into the AMT2 subfamily.
The genome of the moss P. patens revealed a higher and more diversified number of ammonium transporters, with six members of the AMT1 subfamily and seven members of the AMT2 subfamily (Table 1). AMT paralogs often differ with respect to affinities for ammonium and regulation. The existence of multiple paralogs may either help to increase transport capacity, ensure flexibility for regulation, and help cover a wider range of ammonium levels. In both AMT1 and AMT2 clades, the members of the Physcomitrella AMTs group closely together, suggesting a recent origin of the paralogs within the group (Figure 1; Figure A3 in Appendix). Three of the AMT1 members (named PpAMT1;1, PpAMT1;2, and PpAMT1;3), highly similar, correspond to recent gene duplications as they occur in tandem at the same genetic locus. The persistence of several paralogs in Physcomitrella is not limited to the AMT family, and in fact it seems to be a common feature derived from genome duplication events (Lang et al., 2005;Rensing et al., 2007).
The unicellular alga C. reinhardtii has eight AMTs, all belonging to the AMT1 family, and sharing only 31-46% identity with AMT1s of Arabidopsis or Selaginella. Two of them (CrAMT-E and -G) also show alternative splicing (Fernandez and Galvan, 2007). Regarding phylogeny, they fall into three subclades, separated from other plant genes. This divergence is emphasized also by the lack of the conserved regulatory threonine in the C-terminus (Loqué et al., 2007). Other residues, usually conserved in AMT1s but missing or misplaced in some Chlamydomonas AMTs are the N-terminal cysteines, which have been shown to be important to promote interactions among monomers in the trimeric complex  (Graff et al., 2011). It will be interesting to study how AMTs in this organism assemble and are regulated. Apart from the AMT genes, the genome of Chlamydomonas also contains homologs of the animal Rhesus ammonia transporters (Rh), that are distantly related to the AMT family (Soupene et al., 2004). Rh genes are however absent in Selaginella, Physcomitrella, and in any plant species sequenced so far. An interesting observation is that the proportion of AMT1 and AMT2 genes in the various species is different ( Table 1). The Arabidopsis genome possesses five AMT1 members and one AMT2, but in rice the proportion is inverted (three AMT1s and nine AMT2s). Similarly, the Selaginella genome is more enriched in AMT2s (three members) and a single AMT1. Similar to poplar, Physcomitrella is more balanced, with six AMT1 and seven AMT2 types; finally, Chlamydomonas is so far the species containing more AMT1 members (eight) and totally lacking the AMT2 type. It is interesting to observe that the C-terminus of the AMT1 family is highly conserved also in Selaginella and Physcomitrella (but not Chlamydomonas) as is the equivalent of threonine-460 in the Cterminus AtAMT1;1, which is involved in ammonium-dependent allosteric regulation (Loqué et al., 2007;Lanquar et al., 2009), suggesting that the regulatory mechanism via phosphorylation of AMT1 members is maintained in all multicellular plants.

UREA TRANSPORTER GENE FAMILY
Analysis of the Selaginella genome suggests the presence of a unique DUR3 homolog (named SmDUR3), similar to the one found in Arabidopsis, poplar, or rice genomes ( Table 1). Whereas vascular plants have a single urea transporter, the moss Physcomitrella possesses two DUR3 homologs (PpDUR3A and B), with 85% identity between each other. Going further back in the evolutionary tree, Chlamydomonas has three different copies, named DUR3A, B and C, more closely related to each other (about 80% identity) than with the DUR3 members of multicellular plants (Figure 2; Figure A4 in Appendix). Other algae, as well as fungi, are known to have additional copies of DUR3 genes (Kakinuma et al., 2008;Morel et al., 2008).
The progressive loss of genes for urea transporters could indicate that the evolution of the vascular system shifted the preference Frontiers in Plant Science | Plant Evolution and Development for other nitrogen sources, such as ammonium, nitrate, and amino acids, all present in the xylem sap of vascular plants (Smirnoff and Stewart, 1985;Atkins, 2000;Schjoerring et al., 2002).

CONCLUSION
The comparative analysis of ammonium transporters revealed that different lineages of plants vary in the number of AMTs belonging to either family, and in many cases they show a bias toward one of the two AMT subfamilies, with some organisms particularly enriched in the plant-specific AMT1 clade (Arabidopsis, Chlamydomonas), while others rely more on the AMT2 family, typical for bacteria and fungi. The physiological and evolutionary relevance of these shifts remains to be determined. The number of urea transporters decreased during evolution, maybe indicating a preference of vascular plants for other nitrogen sources or the evolution of yet uncharacterized urea transporters.

ACKNOWLEDGMENTS
This work was supported by a grant from the National Science Foundation (MCB-1021677) to Wolf B. Frommer.