Phylogenetic Analysis of K+ Transporters in Bryophytes, Lycophytes, and Flowering Plants Indicates a Specialization of Vascular Plants

As heritage from early evolution, potassium (K+) is absolutely necessary for all living cells. It plays significant roles as stabilizer in metabolism and is important for enzyme activation, stabilization of protein synthesis, and neutralization of negative charges on cellular molecules as proteins and nucleic acids. Land plants even enlarged this spectrum of K+ utilization after having gone ashore, despite the fact that K+ is far less available in their new oligotrophic habitats than in sea water. Inevitably, plant cells had to improve and to develop unique transport systems for K+ accumulation and distribution. In the past two decades a manifold of K+ transporters from flowering plants has been identified at the molecular level. The recently published genome of the fern ally Selaginella moellendorffii now helps in providing a better understanding on the molecular changes involved in the colonization of land and the development of the vasculature and the seeds. In this article we present an inventory of K+ transporters of this lycophyte and pigeonhole them together with their relatives from the moss Physcomitrella patens, the monocotyledon Oryza sativa, and two dicotyledonous species, the herbaceous plant Arabidopsis thaliana, and the tree Populus trichocarpa. Interestingly, the transition of green plants from an aqueous to a dry environment coincides with a dramatic reduction in the diversity of voltage-gated potassium channels followed by a diversification on the basis of one surviving K+ channel class. The first appearance of K+ release (Kout) channels in S. moellendorffii that were shown in Arabidopsis to be involved in xylem loading and guard cell closure coincides with the specialization of vascular plants and may indicate an important adaptive step.


INTRODUCTION
The absolute requirement for K + in all living cells was already fixed from the cradle of evolution in the sea. Among all the cations that were present in the marine environment K + was utilized by cells as the major cation for essential functions as maintaining electroneutrality and osmotic equilibrium. Further evolutionary steps in the cellular K + -rich environment then employed K + as regulator of protein activities being essential for several biochemical processes. The interactions of potassium with these proteins depend on the unique electrochemical properties of K + ions, i.e., the topology of their electrical charge-density. These features cannot or only incompletely be mimicked by Na + or by any other cation because they all differ from K + in their electron shell configuration and consequently also in the arrangement of the surrounding hydration shell. K + thus became indispensably necessary for living cells; a dependency also inherited to Embryophyta, where K + can contribute up to 10% of the dry mass (Leigh and Wyn Jones, 1984). Terrestrial plants even developed new functions for K + such as turgor-driven processes like stomatal movement, phototropism, gravitropism, and cell elongation (Ashley et al., 2006;Rodriguez-Navarro and Rubio, 2006;Amtmann and Armengaud, 2009;Amtmann and Blatt, 2009;Maathuis, 2009;Szczerba et al., 2009). Embryophyta need to survive in oligotrophic environments where K + is present at much lower concentrations than in sea water; the potassium concentration in normal soil solution (10-100 µM) is considerably variable and about three to four orders of magnitude lower than in the plant. Therefore, not only for potassium homeostasis (the maintenance of a dynamic equilibrium in the cellular K + concentration) but also for K + uptake from the environment and its distribution throughout the organism, plants have to invest energy and need a set of specialized transporter proteins.
Pioneering work by Epstein et al. (1963) proposed that K + uptake from soil into plant cells is mediated by two mechanisms that take advantage of the electrical gradient and/or the proton motive force established by H + -ATPases. One was characterized as a high-affinity system (mechanism I), showing apparent affinities in the range of ∼20 µM, that can transport also Na + when K + is not present. The other one (mechanism II) showed a much lower affinity and provided an increasing contribution from >200 µM to mM external K + concentrations. During the last two decades a variety of potassium-permeable transmembrane transport systems -potentially underlying these two components -were identified at the molecular level. They were classified into five major gene families (Maser et al., 2001(Maser et al., , 2002bSentenac, 2002, 2003;Lebaudy et al., 2007): (i) voltage-gated K + channels, (ii) non-voltage-gated (tandem-pore) K + (TPK) channels, (iii) high-affinity K + transporters of the HAK type, (iv) highaffinity K + transporters of the HKT type, and (v) cation-proton antiporters (CPAs). K + channels likely underlie the experimentally observed low affinity component in plants, whereas HAK transporters contribute to the high-affinity K + uptake component. HKT transporters are responsible for a K + -dependent Na + component (Rodriguez-Navarro and Rubio, 2006). However, this is not a generalized strict separation. Channels contribute to the highaffinity K + uptake component and transporters might contribute under certain conditions also to low affinity transport.
Here we took advantage from the recently published genome of Selaginella moellendorffii (Banks et al., 2011) and prepared an inventory of K + transporters in this lycophyte. We focused especially on transporters of the HAK and HKT type and on K + channels. For information on other potentially K + -permeable transporters such as KEA or CHX belonging to the class of monovalent CPAs we refer to a recent excellent review especially dedicated to these proteins (Chanroj et al., 2012). We are comparing the results from S. moellendorffii with those from the moss Physcomitrella patens, the monocotyledon Oryza sativa, as well as with those from two dicotyledonous species, the Brassicaceae Arabidopsis thaliana, and the tree Populus trichocarpa. Voltage-gated K + channels are also compared with those from Chlorophyta.

POTASSIUM TRANSPORTERS OF THE HAK TYPE
The high-affinity K + (HAK) transporter gene family -also called KT or KUP transporter family -is an ancient large family with members in Bacteria, Archaea, Fungi, Amoebozoa, and probably also in some species of Animalia (Grabov, 2007;Benito et al., 2011). Initially HAK genes have been deduced from plants by their similarity to K + uptake permeases (KUP) from E. coli (Schleyer and Bakker, 1993) and high-affinity K + transporters (HAK) from fungi (Banuelos et al., 1995;Quintero and Blatt, 1997;Santa-Maria et al., 1997;Fu and Luan, 1998;Kim et al., 1998). Several members of this family were shown to function as K + uptake transporters in plants especially when the external potassium concentration was in the low µM range (Gierth et al., 2005;Aleman et al., 2011) indicating that HAK transporters are involved in high-affinity K + uptake. Interestingly, all plant genomes analyzed so far contain genes encoding HAK transporters, while in Bacteria, Archaea, and Fungi they were found only in a subset of species (Grabov, 2007;Benito et al., 2011). The HAK family is the largest family of potential K + transporters in plants and members of this family are expressed in nearly all tested plant tissues suggesting that HAK transporters have a general function in K + supply (Banuelos et al., 2002).
To date, the topology of HAK transporters has neither been determined experimentally nor by in silico predictions; nevertheless, hydropathy profiles of these proteins suggest about 12 putative transmembrane segments and a long hydrophilic COOH-terminal region. In genome-wide screenings, HAK transporter proteins can be pinpointed by the presence of several conserved consensus motifs (see Materials and Methods). Our screenings identified 13 HAKs in Arabidopsis, 27 in rice, 22 in poplar, 18 in P. patens, and 11 in S. moellendorffii (Table 1; see Yang et al., 2009, for comparison). It is likely that the genome of P. trichocarpa contains more genes coding for HAK transporters, because especially in the screening of poplar we discarded partial sequences resulting from pre-mature gene annotation. Phylogenetic analyses allowed subdividing them into six independent groups ( Figure 1A; see Rubio et al., 2000, for initial grouping-into Groups I-IV) and revealed that the last recent common ancestor of all embryophytes had two HAK transporters ( Figure 1B). One of these diverged into current Group II, whereas the other was duplicated at least three times before the origin of tracheophytes. Two early duplication events got lost in the lineage leading to tracheophytes and led to P. patens-specific gene family amplifications (Groups V and VI). HAK transporters in S. moellendorffii spread over the other clades (Groups I-IV).
HAK transporters may not only mediate transport across the plasma membrane. Transient expression of the Ory-sat-HAK10::GFP fusion protein in living onion epidermal cells targeted this protein (from Group II) to the tonoplast (Banuelos et al., 2002); and Ara-tha-KUP/HAK/KT12 (from Group I) was found in the chloroplast proteome (Kleffmann et al., 2004;Peltier et al., 2004). Additionally, HAK transporters may not exclusively transport K + . Phy-pat-HAK1 and Ara-tha-HAK5, for instance, were reported to be permeable also for Cs + (Garciadeblas et al., 2007;Qi et al., 2008); and Phy-pat-HAK13 belonging to Group IV was recently characterized as a high-affinity Na + uptake transporter . We therefore propose for the closely related Sel-moe-HAK1 and Sel-moe-HAK11 from S. moellendorffii similar sodium-transport features. This phylogenetic divergence may indicate a -so far underexplored -diversity of HAK transporters in fine-tuned function of K + uptake and re-distribution, cellular expression, and/or sub-cellular targeting.

POTASSIUM TRANSPORTERS OF THE HKT TYPE
HKTs in plants belong to a family of monovalent cation transporters comprising also the fungal TRKs (K + transporters) and bacterial KtrABs (Na + -dependent K + transporter), for instance (Corratge-Faillie et al., 2010). Proteins of this family share a common structure of four TM-P-TM motifs (every two transmembrane α-helices are connected by ∼30 aa-long pore-forming Frontiers in Plant Science | Plant Evolution and Development Ory-sat-HAK2

POPTR_0010s11100
Pop-tri-HAK22 P segments), which might have evolved from an ancestor related to the bacterial KscA K + channel of Streptomyces lividans (Durell and Guy, 1999; Figure 2A). The plant HKT family comprises transporters that mediate Na + uptake in roots or in other plant organs. They accumulate Na + from the soil and recirculate it throughout the plant. There are two types of plant HKT transporters that can be distinguished by the amino acid sequence of the selectivity filter (the narrowest part of the permeation pathways that selects one ion species over others) of the first TM-P-TM motif: (i) S-S-M and (ii) [T,S,I]-G-L. Two HKTs from O. sativa and A. thaliana belonging to the first type have been well characterized in planta as Na + uptake transporters (Uozumi et al., 2000;Rus et al., 2001;Maser et al., 2002a,c;Berthomieu et al., 2003;Garciadeblas et al., 2003; www.frontiersin.org

FIGURE 1 | Continued Evolutionary relationships among HAK transporters in land plants. (B)
Reconciliation analysis of HAK transporters. The last recent common ancestor of all embryophytes had two HAK transporters. One of these diverged into current Group II, whereas the other was duplicated at least three times before the origin of tracheophytes. Two of these duplications got lost in the lineage leading to tracheophytes forming P. patens-specific groups. Red "D"s at branching points indicate predicted gene duplications.
Gray branches indicate gene losses. Sunarpi et al., 2005;Horie et al., 2007;Xue et al., 2011). The function of the second type has not been studied in plants. Nonetheless, two members of this group, from barley and wheat, mediate K + or Na + uniport or Na + -K + symport -depending on the protein expression level -when heterologously expressed in yeast cells (Haro et al., 2005;Banuelos et al., 2008). Functional expression of this type of transporters in Xenopus oocytes produced similar results with slight variations regarding K + versus Na + permeability, symport activity, and permeability to divalent cations (Rubio et al., 1995;Gassman et al., 1996;Jabnoune et al., 2009;Lan et al., 2010;Horie et al., 2011a;Oomen et al., 2012).
In genome-wide screenings, proteins of the HKT type can be pinpointed by the presence of several conserved consensus motifs (see Materials and Methods). Our screenings identified one HKTcoding gene in Arabidopsis, seven in rice, one in poplar, one in P. patens, and six in S. moellendorffii ( Table 2). Phylogenetic analyses grouped all of them into a single group of orthologs ( Figure 2B) indicating that the most recent common ancestor of all embryophytes comprised a single protein of the HKT type. P. patens has a single extant representative (Phy-pat-HKT1), whereas in tracheophytes several duplication events occurred in different lineages ( Figure 2C). The obviously independent multiplication of HKT-coding genes in rice and S. moellendorffii may be correlated with the affinity of these vascular plants to moisture environments. Probably, a larger variety of Na + /K + transporters provides some advantage for better adaptation. Initially, the HKT family has been partitioned into the two subfamilies one and two, and transporter nomenclature was adjusted accordingly of the type "species HKT subfamily; No" (Platten et al., 2006). Subfamily one gathers transporters with the S-S-M signature in the selectivity filter of the first TM-P-TM motif. Our analysis now reveals that this subfamily division emerged in land plants only after the separation of Lycopodiophyta. Thus, the proposed unified nomenclature cannot be applied to all plant HKT genes. The rules fail, for instance, for HKTs from S. moellendorffii and P. patens (see also Haro et al., 2010).
At the functional level, the six HKTs of S. moellendorffii very likely share properties of the orthologs from other species. They may thus be implicated in K + /Na + recirculation in this vascular plant and could contribute not only to K + transport but in first line to desalination and Na + detoxification.

VOLTAGE-INDEPENDENT K + CHANNELS
Potassium channels play important roles in many physiological aspects of higher plants such as osmoregulation, turgor-driven movements, and ion uptake. It is estimated that K + channels can contribute to more than 50% of the nutritional K + uptake www.frontiersin.org under most field conditions (Spalding et al., 1999;Amtmann and Blatt, 2009). In angiosperms there are two large groups of K + channels: voltage-gated channels, the activity of which is regulated by the transmembrane voltage (Dreyer and Blatt, 2009), and non-voltage-gated K + channels. Non-voltage-gated K + channels form the class of tandem-pore K + (TPK) channels. Functional TPK channels are proposed to form dimers consisting of two identical subunits (Maitrejean et al., 2011). Each subunit is characterized by a structure with four transmembrane domains and two pore-forming loops between the first and second and the third and forth membrane-spanning domain, respectively ( Figure 3A; Voelker et al., 2010). Searching for proteins with the characteristic pore-forming region, in the genome of S. moellendorffii four genes coding for TPK channel subunits could be identified (Table 3). Together with the six TPKs from Arabidopsis (Ara-tha-TPK1-5 and Ara-tha-KCO3), three from rice, ten from poplar, and three from P. patens they could be classed into two groups of orthologs ( Figure 3B). This implies that the ancestor of land plants had already two of these genes. A deeper phylogenetic analysis revealed several duplication events in the two groups, both species-specific and at higher levels ( Figure 3C). A remarkable example in this context is KCO3 from A. thaliana. This subunit lacks the first of the two pore loops and was originally considered as founder of a separate channel family with structural features (TM-P-TM; two transmembrane α-helices, and pore-forming P segment) similar to the simplest class of K + channels from bacteria and animals. It became evident, however, that Ara-tha-KCO3 developed through a very recent evolutionary event involving gene duplication of the Ara-tha-TPK2 gene followed by partial deletion (Marcel et al., 2010;Voelker et al., 2010). And indeed, in line with this concept, neither the genome of S. moellendorffii nor that of P. patens appears to contain genes coding for K + channels of the TM-P-TM type.
TPK channels in plants were reported to be targeted to the vacuolar membrane (Czempinski et al., 2002;Voelker et al., 2006;Latz et al., 2007;Dunkel et al., 2008;Isayenkov et al., 2011a,b). The exception is Ara-tha-TPK4 which has been reported to be targeted also to the plasma membrane (Becker et al., 2004). However, orthologs of Ara-tha-TPK4 were only found in the genus Arabidopsis so far (i.e., A. thaliana and A. lyrata; Voelker et al., 2010) but not in other plant species indicating a rather recent evolutionary event in channel specialization. Therefore, we have justified reasons to hypothesize that TPKs in bryophytes and lycophytes are vacuolar K + channels. Their physiological role, however, remains as speculative as that of TPKs in other plants (Voelker et al., 2010).

VOLTAGE-GATED K + CHANNELS OF THE SHAKER -TYPE
Voltage-gated K + channels are tetrameric proteins built of four α-subunits. One subunit shows usually a structure with six transmembrane domains and one pore loop (5TM-P-TM). The first four transmembrane domains fold into the voltage-sensor module, and the pore loop together with the fifth and sixth transmembrane domains establishes the permeation pathway module ( Figure 4A; Dreyer and Blatt, 2009).
Voltage-gated potassium channels in angiosperms are targeted to the plasma membrane and could normally be grouped into the class of Shaker-like K + channels that subdivides into four functional subgroups: (a) Inward-rectifying (K in ) channels open at membrane hyperpolarization and are responsible for K + uptake. (b) Silent (K silent ) channel subunits assemble with K in subunits and modulate K + uptake channel properties. (c) Weakly rectifying (K weak ) channels are specialized K in channels that show a bi-modal gating behavior. They appear to play a special role in the energy household of vascular tissues. (d) Outward-rectifying (K out ) channel subunits open at depolarizing voltages and mediate K + release, e.g., during xylem loading or stomata closure (see Dreyer and Uozumi, 2011, for a contemporary review). Our screening strategy based on the characteristic pore-forming region identified nine -already known -genes coding for Shaker-like channels in Arabidopsis, eleven in rice, and eleven in poplar. Genes coding for Shaker-like K + channels were also identified in the moss P. patens and in the fern ally S. moellendorffii (Table 4). However, whereas

FIGURE 4 | Continued channel tetramer in top view (right). (B) Evolutionary relationships among
Shaker -like channels in land plants. Extensive functional analyses identified inward-rectifying (K in ) channels, outward-rectifying (K out ) channels, weakly rectifying (K weak ) channels, and silent (K silent ) channel subunits that assemble with K in subunits and modulate K + uptake channel properties. (C) Reconciliation analysis of Shaker -like K + channels. The common ancestor of land plants had a single Shaker -like K + channel. Since then several amplifications have occurred: a bryophyte specific amplification following the split between the lineage leading to P. patens and the tracheophytes, and duplications in the tracheophyte lineage. The common ancestor of tracheophytes had two genes coding for Shaker -like K + channel subunits. One of these got lost in the lineage leading to S. moellendorffii after the split of angiosperms. Red "D"s at branching points indicate predicted gene duplications. Gray branches indicate gene losses. For Pop-tri-Kc06 only partial sequence information was available (indicated by an asterisk).
four K in -like channels were identified in P. patens, despite a very careful screening strategy no such gene could be found in S. moellendorffii. Instead there, a gene coding for a K out channel subunit was discovered as the only Shaker-like channel (Figure 4). This channel comprises all the essential structural features that were shown in the Ara-tha-SKOR K out channel to be responsible for a unique K + sensing property (Johansson et al., 2006). K out channels open upon depolarization but additionally adjust their gating to the prevailing concentration of K + outside. As a consequence, they open only at voltages positive of the K + equilibrium voltage, E K , when the electrochemical driving force is directed outward and so ensure K + efflux regardless of the extracellular K + concentration. This ability to adapt channel gating to the cation concentration outside guarantees an efficient K + release during xylem loading and stomatal closure, for instance, even under varying external K + (from 10 nM to 100 mM; Blatt, 1988;Schroeder, 1988;Wegner and de Boer, 1997;Gaymard et al., 1998;Ache et al., 2000). From analogy we may postulate that the presence of a K out channel in the vascular plant S. moellendorffii and its absence in the non-vascular plant P. patens is correlated with the important evolutionary step of vascularization. In contrast, it is rather difficult to find an explanation for the loss of the K in /K weak /K silent channel branch in S. moellendorffii.

OTHER TYPES OF VOLTAGE-GATED K + CHANNELS IN ALGAE, BRYOPHYTES, AND LYCOPHYTES
In addition to Shaker-like channels the genomes of both, S. moellendorffii and P. patens, contain members of another class of putatively voltage-gated potassium channels (Table 5). These channels show some similarity with large conductance Ca 2+ -activated K + channels ("big K" = BK channels), a channel type that is widely present in animals (including humans) but absent in flowering plants, for instance. BK channels activate in response to membrane depolarization and binding of intracellular Ca 2+ and Mg 2+ (Latorre et al., 2010). These channels are built of αand β-subunits, where -as in Shaker-like channels -four α-subunits form the per se functional permeation pathway-establishing unit and the β-subunits just modulate and fine-tune channel properties. In contrast to Shaker-like channels the BK channel protein consists of seven (instead of six) transmembrane domains (6TM-P-TM structure) that lead to an exoplasmic N-terminus. Also BK-like channel

Pp1s2_170V6
Phy-pat-AKT4 S. moellendorffii 453399 Sel-moe-SmORK α-subunits from P. patens and S. moellendorffii show a larger Nterminal region compared to plant Shaker-like channels. It might thus be speculated that also these proteins fold into a 6TM-P-TM structure instead of the 5TM-P-TM Shaker-like topology. The functional properties and the physiological roles of plant BK-like channels are unknown. In mammalian tissues, BK channels serve as a negative-feedback mechanism for excitatory events that lead to increases in calcium concentration or membrane depolarization. In this way, they play a key role, for instance, in regulating the contractile tone in vascular smooth muscle cells or help to terminate the action potential and thus modulate secretion in chromaffin cells. It might be speculated that -at least in S. moellendorffii -the two BK-like channels could compensate for the absent Shaker-like K in channels in carrying out functions in K + uptake and distribution. To assess the evolutionary origin of the non-Shaker-like channels we screened the genomes of the green algae Chlamydomonas reinhardtii, Coccomyxa sp.C-169, Micromonas sp. RCC299, Ostreococcus tauri, and Volvox carteri for voltage-gated K + channels. Despite the fact that this transporter class in algae exhibits a huge structural diversity comprising also homologs of plant Shaker-like www.frontiersin.org channels (Figure 5), a clear trace leading to BK-like channels in S. moellendorffii or P. patens could not be identified. The two most similar channels from Volvox and Chlamydomonas share an identity of 16-19% over a stretch of ∼500 amino acids. In comparison, a BLAST search at NCBI 1 limited to a query coverage of >45% resulted as best hit outside Animalia, S. moellendorffii, or P. patens in a voltage-gated K + channel from Phytophthora infestans (XM_002998337) with ∼28% identity over a stretch of ∼500 amino acids. Unfortunately, from all these results we cannot resolve unequivocally the origin of BK-like channels in S. moellendorffii and P. patens. Neither can we exclude the possibility that bryophytes and lycophytes may have acquired these K + channel genes from Fungi or Protozoa. However, our data ( Figure 5) 1 http://blast.ncbi.nlm.nih.gov/Blast.cgi FIGURE 5 | Voltage-gated K + channels in Chlorophyta, Bryophyta, and Lycophyta. In contrast to land plants, voltage-gated K + channels in algae show a large structural diversity. The functional variety of K + channels in higher plants (Figure 4; Shaker -like K + channels) developed from only one of these channel types.
clearly indicate that the diversity of voltage-gated K + channels observable in Chlorophyta collapsed contemporaneously with the transition of green plants from an aqueous to a dry environment. In higher plants the subsequent functional diversification into K in /K out /K weak /K silent (Figure 4) took place on the basis of only one surviving channel class.

SUMMARY
In the haploid genome of the spike moss S. moellendorffii we identified 1 homolog of voltage-gated outward-rectifying K + release channels, 4 homologs of voltage-independent tandem pore K + channels, 2 homologs with some similarity to large conductance Ca 2+ -activated K + (BK) channels, 11 homologs of transporters of the HAK type, as well as 6 homologs of the HKT type ( Table 6). On the basis of phylogenetic analyses, detailed functional properties can be predicted for a few of them. Most probable is that Selmoe-SmORK forms voltage-gated K + release channels involved in stomatal closure and/or in K + loading into the vascular bundles.

GENOME-WIDE SEARCH FOR K + TRANSPORTERS
Putative K + transporters were identified using the conceptual proteomes of A. thaliana (TAIR10 Genome release), O. sativa spp. Indica, P. trichocarpa, P. patens, and S. moellendorffii (Phytozome v6.0) and the algae genomes Coc-comyxa_C169, Micromonas RCC299, O. tauri (v2, v3, and v4 respectively 2 ), V. carteri, and C. reinhardtii (Phytozome v8.0) by screening with different transporter class-specific protein motifs: three motifs for K + channels ( (1) (Rice et al., 2000). Additionally, results were checked against BLAST searches in the five genomes using known transporters of different classes from Arabidopsis and rice as templates. In order to eliminate false-positives the resulting raw-data were curated in a semi-automatic way. In a first step sequences with a length <70% of the average length between the outermost motifs in the corresponding Arabidopsis transporters were discarded. Subsequently, the remaining n protein sequences of each transporter type of each species were pairwise aligned using ClustalW2 3 . From the resulting n(n−1)/2 pairs those with a score of <20 and of 100 (identical sequences) were removed. The residual pairs fragmented the sequences into distinct groups. That group with the highest similarity to the corresponding Arabidopsis transporters was selected for further analyses.
To verify whether the screening for K + channels in S. moellendorffii was exhaustive, its genome was screened in the six-frame Table 6 | Summary -Molecular toolkit for K + uptake and re-distribution in S. moellendorffii.
11 Genes coding for transporters of the HAK type 6 Genes coding for transporters of the HKT type 4 Genes coding for TPK subunits that dimerize into non-voltage-gated K + channels 1 Gene coding for a subunit that homotetramerizes into voltage-gated K + release (K out ) channels 2 Genes coding for subunits that tetramerize into channels with some similarity to animal large conductance Ca 2+ -activated K + ("big K" = BK) channels translations using the program SIXPACK from EMBOSS (Rice et al., 2000) for the presence of the K + -selectivity filter motif G-Y-G in ORFs. Following a positive hit, the closer environment of the GYG was inspected manually for further characteristic sequence features allowing categorizing the peptide to be part of a K + channel. As a result -besides the K + channels obtained already in the first screening -only the K out channel SmORK could be identified in addition.

PHYLOGENETIC ANALYSES
Sequences from each family were aligned using MAFFT (Katoh and Toh, 2010), and alignments were filtered using GBlocks (Castresana, 2000) in order to eliminate regions of low quality. Briefly, the minimum number of sequences for a conserved position was half the number of sequences, the minimum number of sequences for a flanking position was half the number of sequences, the maximum length of contiguous non-conserved positions was 20, and the minimum length of a block was two, positions with gaps were not treated differently from other position. Evolutionary relationships were inferred by Maximum Likelihood using RAxML and 1000 bootstrap replicates (Stamatakis, 2006). The evolutionary model used for phylogenetic analyses was inferred using ProtTest (Darriba et al., 2011). For two pore channels and HAK transporters the model was LG + Γ, for HKTs and Shaker-like channels it was JTT + Γ. In order to root and resolve the gene trees we performed a gene tree-species tree reconciliation analysis using the species tree from Lang et al. (2010;TreeBase 10409). Reconciliation analysis was carried out in Notung 2.6 (Chen et al., 2000;Vernot et al., 2007). To get an idea of the phylogenetic structure of the other voltage-gated K + channels displayed in Figure 5, the sequences were hierarchically clustered based on pairwise identities between two sequences using UPGMA (Unweighted Pair Group Method with Arithmetic Mean). UPGMA analyses were carried out in MAFFT 4 .

ACKNOWLEDGMENTS
This work was supported by grants from the Spanish Ministerio de Economía y Competitividad to Ingo Dreyer and Alonso Rodríguez-Navarro (BFU2011-28815; AGL2007-61705), a Marie Curie Career Integration Grant to Ingo Dreyer (FP7-PEOPLE-2011-CIG No. 303674 -Regopoc), as well as by a Marie-Curie Cofund fellowship to Judith Lucia Gomez-Porras. Kamil Sklodowski is a recipient of a doctoral fellowship from the Max-Planck Research School "Primary Metabolism and Plant Growth."

SUPPLEMENTARY MATERIAL
The alignments used for generating the phylogenetic trees presented in this study are available online as Supplementary Material.