HAK/KUP/KT (High-affinity K+ transporters/K+ uptake permeases/K+ transporters) is the largest potassium transporter family in plants, and plays pivotal roles in K+ uptake and transport, as well as biotic and abiotic stress responses. However, our understanding of the gene family in barley (Hordeum vulgare L.) is quite limited.
In the present study, we identified 27 barley HAK/KUP/KT genes (hereafter called HvHAKs) through a genome-wide analysis. These HvHAKs were unevenly distributed on seven chromosomes, and could be phylogenetically classified into four clusters. All HvHAK protein sequences possessed the conserved motifs and domains. However, the substantial difference existed among HAK members in cis-acting elements and tissue expression patterns. Wheat had the most orthologous genes to barley HAKs, followed by Brachypodium distachyon, rice and maize. In addition, six barley HAK genes were selected to investigate their expression profiling in response to three abiotic stresses by qRT-PCR, and their expression levels were all up-regulated under salt, hyperosmotic and potassium deficiency treatments.
Twenty seven HAK genes (HvHAKs) were identified in barley, and they differ in tissue expression patterns and responses to salt stress, drought stress and potassium deficiency.
Potassium (K) is the second most abundant mineral nutrient in plants, comprising 2–10% of plant dry weight . K is not only essential for plants to maintain normal physiological and biochemical processes such as stomatal movement, photosynthesis, osmoregulation, protein synthesis, enzyme activation, but also involved in the responses to biotic and abiotic stresses [2, 3]. K in soil exists in several different states/pools [4, 5], and available K for plants is commonly quite low due to the strong K adsorption by 2:1 silicate minerals . In plant cells, the cytosolic K+ concentration is relatively constant, ranging from 60 to 200 mM, while vacuolar K+ concentration is more variable, ranging from 20 mM (in K+-deficient plants) to 500 mM (in specialized guard cells) [1, 6, 7]. In addition, apoplastic K+ concentration also shows a large variation, ranging from 10 to 200 mM or even up to 500 mM . In short, K concentration in plants differs greatly between species, tissues, and cell organs.
As sessile organisms, plants have evolved the efficient K+ transport systems to maintain optimal growth under the conditions of highly variable K+ levels . In plants, K+ is first taken up by roots, and then translocated to aerial parts and distributed within cells into different compartments . All these processes are involved in K+ transmembrane transport, which is achieved mainly by K+ channels and transporters. So far, K+ uptake systems are better identified and characterized in Arabidopsis thaliana than in any other plant species. Currently it is well documented that non-selective cation channels (NSCC) probably play the major role in K+ uptake at high concentrations (> 10 mM), while the low-affinity inward rectifying channel AtAKT1 dominates K+ uptake at intermediate concentrations (1 mM). Under lower external K+ concentrations (100 μM), high-affinity potassium transporter AtHAK5, and together with AtAKT1, contributes to K+ acquisition. AtHAK5 is the only transporter capable of taking up K+ under extremely low concentrations (< 10 μM) . Thus, it is imperative to characterize the function and expression pattern of high-affinity potassium transporters.
There are three families of K+ transporters in plants: HAK/KUP/KT (High-affinity K+ transporters/K+ uptake permeases/K+ transporters), HKT (High-affinity K+ transporters) and CPA (monovalent cation/proton antiporters) . HKT transporters function as K+/Na+ co-transporters or Na+ transporters . CPA family can be further divided into NHX (Na+/H+ exchangers), CHX (Cation/H+ exchangers) and KEA (K+ efflux antiporters) subfamilies , and only several members (viz. AtCHX17, AtKEA4, AtKEA5 and AtKEA6) have been proved to function in K+ transport [12, 13]. The HAK/KUP/KT members in bacteria and fungi were named as KUP (K+ uptake permease) and HAK (High-affinity K+ transporter), respectively [14, 15]. Diverse acronyms were assigned when the family members in plants were identified, viz. HAK for barley , KT (K+ transporter) or KUP for Arabidopsis [17–19], so composite names HAK/KUP/KT or KT/HAK/KUP or KT/KUP/HAK were widely used for the whole transporter family in plants [10, 20, 21]. HAK/KUP/KT (HAK hereafter) is the largest K+ transporter family, and its members are predicted as K+/H+ symporters . HAK genes are absent in animal cells, but present in all known plant genomes , and first identified in barley  and Arabidopsis [17–19], and later also in rice [23, 24], wheat  and maize . However, the mechanisms of HAK transporters in mediating K+ uptake and transport are still completely unclear.
HAK transporters play the diverse roles in K+ uptake and translocation, salt and drought stress response, as well as morphological development of root and shoot . In Arabidopsis, AtHAK5 is a major K+ deprivation-induced high-affinity K+ uptake transporter, and is the only system being able to take up K+ at external concentrations below 10 μM [9, 27]. AtKUP7 is responsible for K+ uptake at high K+ concentrations, and may be also involved in K+ transport into xylem sap, affecting K+ translocation from roots to shoots . In rice, the expression of OsHAK1 is up-regulated by K+ deficiency in roots. Knockout of OsHAK1 reduces total K+ uptake by approximately 80 and 65% at 0.05–0.1 mM and 1 mM of K+ concentrations, respectively, while overexpression of OsHAK1 increases K+ uptake and K+/Na+ ratio . OsHAK5 is expressed highly in root epidermis, stele and vascular tissues, and plays major roles not only in K+ acquisition, but also in K+ upward translocation from roots to shoots under low external K+ conditions . K+ homeostasis and K+/Na+ balance are severely disturbed when plants are exposed to salt stress, as high Na+ concentrations suppress the expression of HAK genes at low K+ concentrations [30, 31]. However, AtHAK5 also plays a major role in maintaining high-affinity K+ acquisition and plant growth under salt stress . OsHAK1 is up-regulated by salt stress under normal K+ supply condition, but down-regulated at low K+ concentrations . The net K+ uptake rate of oshak1 mutant was almost completely blocked by salt stress when K+ concentration was below 0.05 mM, demonstrating that OsHAK1 plays pivotal role in enhancing salt tolerance in rice . Overexpression of OsHAK5 also increases shoot K+/Na+ ratio and salt tolerance in rice . Moreover, OsHAK1 knockout mutants reduce drought tolerance and growth at both vegetative and reproductive stages, while the over-expressed lines enhance drought tolerance and increase grain yield by 35% relative to wild type under drought condition . In addition to K+ uptake and translocation, and stress tolerance, some HAK genes are also associated with tissue (root and shoot) morphology. It has been reported that AtHAK5 is the downstream target of ARF2 (Auxin response factor 2) and its overexpression results in longer primary roots under low K+ conditions . AtKUP4/TRH1 (Tiny root hair 1) is required for polar localization of PIN1 (An auxin efflux transporter) in root apex, and for gravitropic responses, as well as normal initiation and formation of root hairs [34–37]. AtKUP2/SHY3 (Short hypocotyl 3) mutation affects cell expansion and leads to developmental defects in shoots . In short, the versatile roles of HAK genes have been recognized, including K+ transport and abiotic stress response, but little has been known about the relevant molecular mechanisms.
Barley (Hordeum vulgare) ranks the fourth largest cereal crop worldwide with multiple use as animal feed, human food and brewing material . In comparison with other cereal crops including rice, wheat and maize, barley is well known for its higher salt tolerance, which is largely attributed to its more K+ uptake in roots, more Na+ exclusion and vacuolar sequestration [40–42]. Although genome-wide identification of HAK gene family has been accomplished in rice [24, 43], wheat , maize , cassava , peach , pear , and willow , the similar research has not been done in barley. The release of barley genome data  and recently improved annotated reference genome assembly  allow it possible to perform genome-wide identification of HAK gene family in barley. In the present study, we identified 27 HAK genes (HvHAKs) from barley genome, and analyzed their phylogenetic relationships, conserved motifs and domains, gene structure, cis-acting elements, syntenic relationships, tissue expression patterns and expression profiling in responses to salt stress, osmotic stress and potassium deficiency.
Identification of HvHAKs in barley
A total of 27 HvHAK genes were identified (Table 1; Additional file 1) using barley genome data  and recently improved annotated reference genome assembly . The 27 HvHAK genes were unevenly distributed on the seven chromosomes, with chromosome 2 containing eight genes, chromosomes 3 and 7 each containing four genes, chromosomes 1, 5 and 6 each containing three genes, and chromosome 4 containing two genes (Additional file 2). HAK genes in barley were named in the order of their locations on chromosomes (viz. HvHAK1 to HvHAK27), and no tandem duplication event was observed (Table 1; Additional file 2). The length of HvHAK transporters varied from 724 to 875 amino acids (aa), and the number of transmembrane segments ranged from 10 to 14, with 11–12 (66.7%) being the most abundant (Table 1). The isoelectric points (pI) of HvHAK proteins ranged from 5.65 to 8.99, with only 6 members being below 7 (Table 1). The theoretical molecular weights (MW) of HvHAK proteins were in the range of 79.93 to 98.04 kDa. All the HvHAK transporters contained a typical “k_trans” domain (PF02705), which is specific to HAK potassium transporter family members (Table 1). In addition, all HvHAKs were localized on plasma membrane (Table 1).
Phylogenetic analysis of HvHAK transporters
The protein sequences of 13 AtKUPs in Arabidopsis , 27 OsHAKs in rice , 27 ZmHAKs in maize , and together with 27 HvHAKs in barley, were aligned with MAFFT (https://www.ebi.ac.uk/Tools/msa/mafft/) and used to construct phylogenetic tree by maximum-likelihood (ML) method (Fig. 1). Based on the classification criteria of HAK transporters in Arabidopsis, rice and maize, the identified 27 barley HvHAKs could be classified into four clusters, and each of them could be further subdivided into sub-cluster A and B (Fig. 1). Clusters I and II each contained nine transporters, with IA, IB, IIA and IIB each containing two (HvHAK7 and 23), seven (HvHAK3, 8, 9, 12, 15, 17 and 21), four (HvHAK2, 5, 14 and 27) and five members (HvHAK4, 16, 22, 24 and 26), respectively (Fig. 1). Cluster III consisted of six transporters, with three members in both IIIA (HvHAK1, 10 and 20) and IIIB (HvHAK6, 11 and 18) (Fig. 1). Only three HvHAK transporters were categorized into cluster IV, with one member (HvHAK19) in IVA and two members (HvHAK13 and 25) in IVB (Fig. 1).
Motif, domain and gene structure analyses of HvHAKs
The conserved motifs and domains of HvHAK protein sequence and structures of HvHAK genes were further analyzed (Fig. 2). Phylogenetic analysis confirmed the classification of 27 HvHAKs (Figs. 1 and 2a). Ten conserved motifs were identified in all 27 HvHAK sequences (Fig. 2b). According to NCBI conserved domain database, there are eight types of members for “K_trans superfamily” (cl15781), and three of them were identified in HvHAK protein sequences (Fig. 2b). The conserved domains of three, two and six HvHAK proteins were specifically identified as “K_trans”, “PLN00151” and “kup”, respectively, while those of the other 16 HvHAK sequences were identified as “K_trans superfamily” because of relatively low specificity (Fig. 2b). Cluster I and III contained “kup” and “K_trans superfamily”, “PLN00151” and “K_trans superfamily”, respectively; while cluster II and IV were solely comprised of “K_trans superfamily” and “K_trans”, respectively (Fig. 2b). Moreover, gene structures of barley HvHAKs were highly variable, containing 3–10 exons and 2–9 introns (Table 1; Fig. 2c).
Totally there were 14 (51.9%) HvHAK genes which each contained eight introns, and there were 4, 3 and 2 genes which contained nine, six and seven introns, respectively; and other 4 genes possessed five, four, three and two introns, respectively (Fig. 2c). Intron numbers also differed among the different clusters (Fig. 2c); cluster I was the most variable, with six genes each containing two, three, four, six, eight and nine introns, respectively; three genes in cluster II had six, eight and nine introns, respectively; and two genes in both cluster III and IV had eight and nine, five and seven introns, respectively (Fig. 2c). These introns could be classified into three types based on phases: phase 0, phase 1 and phase 2 . Interestingly, phases 0 and 2 introns were detected, while phase 1 introns were not found in barley HvHAK genes (Fig. 2c). The phylogenetically most closely related members are prone to have the similar exon/intron structure. For instance, gene pairs HvHAK6 and HvHAK11, HvHAK1 and HvHAK10, HvHAK5 and HvHAK14, and HvHAK8 and HvHAK17 contained the same exon/intron numbers, respectively, and the exon configuration of each gene pair was almost the same regardless of the difference in intron length (Fig. 2c). However, the gene pairs of HvHAK4 and HvHAK16, and HvHAK7 and HvHAK23 showed the difference in both intron number and exon configuration (Fig. 2c).
Identification of cis-acting elements in HvHAK genes
A 2 kb sequence in the upstream of HvHAK coding sequences was used for identification of cis-acting element. In total, 52 cis-acting elements were identified, and categorized into 6 types according to their functional annotations (Fig. 3; Additional file 3). Light responsiveness-related elements were the most abundant (22, 42.3%), followed by hormone response-related (9, 17.3%), development/tissue specificity-related (8, 15.4%), promoter/enhancer element-related (7, 13.5%) and stress-related (5, 9.6%) ones (Fig. 3). Only one kind of elements was identified as circadian control (Fig. 3). CAAT-box and TATA-box from promoter/enhancer element, which are binding sites for RNA polymerase and responsible for transcription efficiency, were ubiquitously identified in all HvHAKs, indicating that they might play the pivotal roles in controlling the initiation and level of HvHAKs expression. Moreover, ABRE, as well as CGTCA-motif and TGACG-motif (from hormone response type), which are involved in abscisic acid (ABA) and methyl jasmonate (MeJA) responsiveness, respectively, were widely present in HvHAKs (Fig. 3). G-box in light responsiveness was also extensively distributed in promoter regions of HvHAKs (Fig. 3). These results suggest that HvHAKs might play critical roles in phytohormone metabolism and environmental responses.
Synteny analysis of HvHAKs in four monocotyledons
The syntenic relationships of HvHAKs in four monocotyledonous plants (Triticum aestivum, Brachypodium distachyon, Oryza sativa and Zea mays) were investigated (Fig. 4; Additional file 4). Totally 73 pairs of orthologous genes were identified between barley and wheat, comprising 25 (92.6%; except HvHAK3 and HvHAK9) HvHAKs and 72 TaHAKs (Fig. 4). Twenty-two pairs of orthologous genes were observed between barley and Brachypodium distachyon, barley and rice, and barley and maize, respectively, including 19 (70.4%) HvHAKs and 21 BdHAKs, 17 (63.0%) HvHAKs and 21 OsHAKs, 13 (48.1%) HvHAKs and 20 ZmHAKs, respectively (Fig. 4). Orthologous genes of six barley HvHAKs (HvHAK1, 2, 15, 17, 21 and 23) were not detected in Brachypodium distachyon, rice and maize genomes, although found in wheat genome. Totally, 22 HvHAKs each had three orthologous genes in wheat, while only HvHAK24 had three orthologous genes in maize. Three HvHAKs, five HvHAKs and seven HvHAKs each had two orthologous genes in Brachypodium distachyon, rice and maize genomes, respectively. One TaHAK and one BdHAK were orthologous to two HvHAKs (HvHAK4 and HvHAK16), and one OsHAK and one ZmHAK were orthologous to the other two HvHAKs (HvHAK5 and HvHAK14). Besides, another ZmHAK was also orthologous to HvHAK12 and HvHAK13. These results indicate that HvHAKs in barley are phylogenetically closer to TaHAKs in wheat than to BdHAKs in Brachypodium distachyon, OsHAKs in rice and ZmHAKs in maize.
Tissue expression patterns of HvHAKs
Tissue expression patterns of HvHAKs were investigated in 15 tissues based on transcriptomic data from Barlex database (Additional file 5). Hierarchical cluster analysis revealed that 27 HvHAKs could be divided into two groups according to their expression levels (Fig. 5). One group comprised 14 HvHAKs, and of them 13 genes were from cluster II and cluster III, with only one gene (HvHAK8) from cluster I (Fig. 5). Genes in this group were expressed ubiquitously in all 15 tissues with relatively high levels. The other group consisted of 13 genes, with eight from cluster I, two from cluster II and three genes from cluster IV (Fig. 5). The expression levels of HvHAKs in this group were comparatively lower than that of the genes belonging to the former group, with eight genes being not detected in certain specific tissues (Fig. 5). Moreover, the expression of HvHAK2, HvHAK3 and HvHAK21 was undetectable in seedling roots, while detected in the roots at 28 d after heading (Fig. 5). Furthermore, HvHAK13 was not expressed, and HvHAK15 and HvHAK21 were mildly expressed in developing grains (5 d after heading). However, HvHAK13 was mildly expressed, but HvHAK15 and HvHAK21 were not expressed at 15 d after heading (Fig. 5). These results indicate that the expression of HvHAKs is not only tissue-dependent, but also varies with development phase.
Expression patterns of HvHAKs responding to abiotic stresses
The effects of salt stress (200 mM NaCl), drought stress (20% PEG8000) and potassium deficiency (0.01 mM K+) on barley growth were investigated. All three abiotic stress treatments inhibited shoot growth of barley seedlings, and resulted in abnormal root phenotypes (Additional file 6).
Under salt stress, all HvHAKs examined were up-regulated, but their expression patterns differed (Fig. 6). HvHAK5, HvHAK6, HvHAK18 and HvHAK27 were rapidly induced after one-hour salt treatment, and displayed continuous up-regulation, while HvHAK8 was not induced until salt treatment for 3 h (Fig. 6a). Besides, HvHAK17 was first down-regulated after treatment (1–6 h), and was up-regulated after one-day salt treatment (Fig. 6a). In addition, HvHAK6, HvHAK18 and HvHAK27 showed significantly higher up-regulation after long-term (3–6 d) salt treatment than after short-term (1 h to 1 d) treatment (Fig. 6a).
Under osmotic stress, HvHAK5, HvHAK6, and HvHAK27 were continuously induced during the whole experiment period (Fig. 6b). However, HvHAK8, HvHAK17 and HvHAK18 remained little change after one-hour treatment and even down-regulated after three-hour treatment, while they were up-regulated after treatment for 6 h and 1 d (Fig. 6b). Similarly, HvHAK5, HvHAK6 and HvHAK8 had significantly higher expression after long-term (3 d and 6 d) treatment than after short-term (< 1 d) treatment (Fig. 6b).
All HvHAK genes examined were up-regulated under K deficiency, except HvHAK6, which was down-regulated at 1 d after treatment (Fig. 6c). Besides, the expression of HvHAK18 remained relatively stable during the whole period of treatment (Fig. 6c).
So far, the HAK/KUP/KT potassium transporter family has been identified and analyzed in many plant species, including Arabidopsis [50, 51], rice [23, 24, 43], wheat , maize  and cassava . In the current study, we identified 27 HvHAK genes in barley (Table 1) using data of barley genome  and the improved annotated reference genome assembly , and then investigated their phylogenetic relationships, sequence characteristics, cis-acting elements, syntenic relationships, tissue expression patterns and expression profiling in response to abiotic stress.
Nomenclature and classification of HvHAKs
HAK genes in maize and wheat were named according to their phylogenetic relationships with OsHAKs [25, 26]. The first identified HAK gene in barley was named as HvHAK1 due to its high homology to Escherichia coli Kup and Schwanniomyces occidentalis HAK1, and another gene with only partial length (incomplete at 5′ end) was named as HvHAK2 . Rubio et al.  isolated three new cDNAs of HAK genes in barley, and found one of them had more than 90% similarity to HvHAK1. Consequently, they named it as HvHAK1B (HvHAK1 was designated as HvHAK1A), and the other two genes as HvHAK3 and HvHAK4, respectively. HvHAK1/HvHAK1A and HvHAK1B were phylogenetically most closely related to OsHAK1, while HvHAK2, HvHAK3 and HvHAK4 showed the closest phylogenetic relationships with OsHAK7, OsHAK23 and OsHAK2, respectively (Table 1; Fig. 1). Therefore, naming HvHAKs based on their phylogenetic relationships with OsHAKs seems not to be straightforward and understandable. Thus, 27 HvHAKs were designated as HvHAK1 to HvHAK27 according to their chromosomal locations (Table 1).
OsHAKs were divided into four clusters on the basis of phylogenetic analysis, with cluster I, II, II and IV having eight, nine, six and four genes, respectively . Barley HAK genes can be also classified into four clusters according to the classification criteria of OsHAKs in rice, with cluster I, II, III and IV having nine, nine, six and three genes, respectively (Fig. 1). Notably, OsHAK13 was subdivided into cluster IIB previously , while it was subdivided into cluster IIA in this study (Fig. 1). This inconsistency might be caused by the difference in the methods of protein sequence alignment and phylogenetic construction. Gupta et al.  adopted ClustalX in the protein sequence alignment of OsHAKs, while MAFFT was used in this study, which was widely recommended for its higher accuracy in multiple sequence alignment . In addition, neighbor joining (NJ) was used in the phylogenetic construction of OsHAKs . NJ method is widely used due to its high computational efficiency, however, it cannot account for the high variances of large distance estimates and is also sensitive to the gaps in the sequence alignment . Maximum likelihood, which was adopted in this study, has an obvious advantage over NJ method in understanding the sequence evolution process .
Similarities and differences of HAK genes between barley and other plants
In this study, 27 HAK genes were identified in barley, basically similar to those found in rice , maize  and Brachypodium distachyon , although these plants show a huge difference in genome sizes, ranging from ~ 355 Mb of Brachypodium distachyon to ~ 5.3 Gb of barley. In addition, the chromosome distribution of 27 HvHAK genes are uneven. Similarly, ZmHAKs were unevenly distributed on 10 chromosomes . However, the chromosome distribution of OsHAKs are almost even . The duplication events (including segmental duplication and tandem duplication) could be observed in both OsHAKs and ZmHAKs [26, 43], but not detected in HvHAKs (Additional file 2). The length in amino acid sequences of HvHAKs, OsHAKs and ZmHAKs were similar, ranging from 724 to 875, 697 to 877 and 642 to 921 aa, respectively (Table 1) [24, 26]. Introns of HvHAK genes varied from 2 to 9, being similar to those of ZmHAK and OsHAK genes, with an exception of OsHAK22, which had only one intron. Interestingly, the HAK genes with 8 introns were predominant in barley (14, 51.9%), rice (12, 44.4%) and maize (16, 59.3%) [24, 26]. Intron phase determines whether exons are targeted for alternative splicing. Generally, exons flanked by same-phase introns are possible for alternative splicing, while those flanked by the different-phase introns cannot be alternatively spliced . On the whole the most introns of HAK genes in rice and maize belonged to phase 0 and phase 1, with only four and three introns in phase 2 for OsHAKs and ZmHAKs, respectively [24, 26]. However, introns of HvHAKs all belong to phase 0 or phase 2, and no introns were observed at phase 1 (Fig. 2), suggesting that HvHAKs had the different expression over the phase from OsHAKs and ZmHAKs. The numbers of transmembrane segments (TMS) for OsHAKs, ZmHAKs and HvHAKs were in the range of 11–15, 9–14, and 10–14, respectively (Table 1) [24, 26]. The HAK genes with 13 TMS (11, 40.7%) and 11 TMS (9, 33.3%) were dominant in OsHAKs and ZmHAKs, while those with 11 and 12 TMS (18, 66.7%) were the most abundant in HvHAKs (Table 1) [24, 26].
All HAK proteins in rice, maize and wheat contain a conserved domain, “K_trans” (PF02705) [24–26]. In the present study, a conserved domain of “K_trans superfamily” (cl15781) was also identified in all HvHAK transporters (Fig. 2b). In addition, 1 conserved motifs were detected, and they were almost evenly distributed along HvHAK sequences, similar to motifs in TaHAKs . For the cis-acting elements related to light responsiveness, they were the most abundant in HvHAKs, while quite fewer in OsHAKs and ZmHAKs [24, 26]. On the other hand, Ca2+-responsive elements could be detected in the most OsHAKs and ZmHAKs, but not observed in HvHAKs (Fig. 3). It may be concluded that there is a great difference in cis-acting elements among OsHAKs, ZmHAKs and HvHAKs.
Based on the classification criteria of HAK genes in rice and maize, HvHAKs can be also divided into four clusters (Fig. 1). Synteny analysis identified 73 HAK orthologous gene pairs (comprising 25 HvHAKs and 72 TaHAKs) in barley and wheat genomes, much more than those between barley and Brachypodium distachyon genomes, and between rice and maize genomes, respectively (Fig. 4; Additional file 4). Interestingly, each of the 22 HvHAK genes was orthologous to three TaHAK genes, which might be attributed to the fact that wheat is an allohexaploid composed of three distinct ancestral genomes, viz. A, B and D . If this explanation is reasonable, the phenomena that HvHAK15 and HvHAK16 were orthologous to two (lacking orthologous gene on genome D) and one TaHAKs (lacking orthologous genes on genomes A and B), respectively, can be described to gene loss events during evolution (Additional file 4). In addition, HvHAK4 was orthologous to four TaHAK genes on two chromosomes, suggesting that both gene loss and genomic recombination might have occurred (Additional file 4). Notably, according to the synteny analysis in this study, at least 72 HAK genes should be present in wheat genome, while only 56 TaHAK genes were identified in a previous research . The inconsistency can be attributed to the imperfect assembly and annotation of wheat genome, as well as the workflow of gene family identification. Based on the barley genome assembly released in 2017 , only 24 HAK genes were identified in barley genome (three genes with incomplete K+ transporter domains were arbitrarily discarded), while using the recently released barley reference genome assembly , we identified 27 HvHAK genes. On the other hand, 22 HAK orthologous gene pairs were identified between barley and other three plants examined in this study (Fig. 4; Additional file 4). In addition, 19, 17 and 13 HvHAKs were orthologous to 21 BdHAKs, 21 OsHAKs and 20 ZmHAKs, respectively, being consistent with the phylogenetic relationships between barley and these plants (Fig. 4; Additional file 4). No orthologous gene was found for HvHAK3 and HvHAK9 in other three genomes, indicating that these two genes might originate from ectopic duplication of HvHAK2 and HvHAK8, respectively (Fig. 1) .
Expression of HvHAKs in different tissues and in response to abiotic stresses
Analysis of tissue expression pattern revealed that HvHAK genes could be divided into two groups (Fig. 5). One group, consisting of 14 genes, was constitutively and highly expressed in all 15 tissues, while another group, consisting of 13 genes, only expressed in the specific tissues with lower levels (Fig. 5). In the previous studies, expression of HAK genes in Arabidopsis, rice and wheat also showed the great difference in the tissues and levels [24, 25, 58].
Salt, osmotic (drought) and K deficiency stresses had the dramatic effect on the expression of the examined six HvHAK genes (Fig. 6). Maintaining efficient K+ uptake is prerequisite for K+/Na+ homeostasis and salt tolerance when plants are exposed to salt stress . In this study, the expression of the six HvHAKs was up-regulated under salt stress (Fig. 6a), contributing to high K+/Na+ ratio and salt tolerance. The similar results were also observed for OsHAKs in rice and for MeKUPs in cassava [44, 60]. It is well recognized that plants respond to drought stress via both ABA-dependent and ABA-independent pathways [61, 62]. It was reported that increasing cytosolic ion concentrations through absorbing and accumulating K+ is a relatively fast and cost-effective way for plants to enhance intracellular osmotic adjustment and drought tolerance [63, 64]. In the present study, the expression of HvHAK genes were up-regulated in a time-dependent manner when plants were exposed to osmotic stress (drought) (Fig. 6b), being consistent with the previous reports on TaHAKs in wheat and MeKUPs in cassava [25, 44]. In general, K+ transporter genes are up-regulated under K deficiency conditions [25, 60, 65]. In this study, all examined HvHAK genes were up-regulated under K deficiency treatment (Fig. 6c). It was reported that the expression of OsHAK1 was induced by K deficiency or salt stress, leading to increased K+ uptake and K+/Na+ ratio in roots . In this study, the homolog genes of OsHAK1 in barley, HvHAK8 (alternative name as HvHAK1 or HvHAK1A) and HvHAK17 (alternative name as HvHAK1B), were also up-regulated under K deficiency or salt stress (Table 1; Figs. 1 and 6). In addition, the expression level of HvHAKs varied dramatically with the time of stress exposure (Fig. 6). For example, HvHAK5 was slightly up-regulated (1.2–2.4 folds) at 1 h to 1 d after osmotic stress, while it was up-regulated significantly (5.8–6.6 folds) at 3–6 d after the treatment (Fig. 6b). Moreover, three genes (HvHAK2, 3 and 21) were highly expressed in roots at 28 d after heading, but remained little change in seedling roots (Fig. 5). The same was true for HvHAK13, HvHAK15 and HvHAK21 in the developing grains (5 DAP and 15 DAP) (Fig. 5). Obviously, the response patterns of HvHAK genes to abiotic stress are dual-phases or even multi-phases. Thus, the suitable time of sampling should be taken into consideration when expression patterns of these HvHAK genes are analyzed and compared.
In the current study, 27 HAK genes (HvHAKs) were identified in barley. The expression of these HvHAKs vary with the plant tissues and the time of exposed abiotic stress. Although all HvHAKs could be induced by salt, drought and K deficiency stresses, their response patterns and magnitudes showed the great variation. The obtained results should be helpful for us to make the comprehensive understanding of HAK family in barley.
Sequence, chromosomal location and duplication analyses of HvHAKs
Conserved motifs of HvHAKs were identified using MEME program (http://meme-suite.org/tools/meme)  with following parameters: classic motif discovery mode, zero or one occurrence per sequence (zoops), motif number was set to 10. Conserved domains were analyzed using NCBI Conserved Domain Database (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi). Motifs and conserved domains of HvHAKs as well as structures and chromosomal locations of HvHAKs were visualized using TBtools . Gene duplication events were analyzed using TBtools following parameters described by Tombuloglu [69, 71].
Plant materials, growth conditions and abiotic stress treatments
Barley (Hordeum vulgare cv. Morex. From Professor Rugen Xu’s lab, Yangzhou University) seeds were sterilized with 10% commercial NaClO for 15 min and rinsed with tap water for 30 min . Sterilized seeds were germinated in basic salt medium (BSM, 0.5 mM KCl + 0.1 mM CaCl2) for 2 d, and then BSM was changed to one-fifth Hoagland solution for another 4 d in growth room with a photoperiod of 14/10 h, light intensity of 200 ± 25 μmol·m− 2·s− 1, temperature of 23/18 °C (day/night) and relative humidity of 60%. Barley seedlings were grown for 6 d and then subjected to salt stress (200 mM NaCl) , osmotic stress (20% PEG8000)  and potassium deficiency (0.01 mM K+)  in background of 1/5 Hoagland solution. The seedlings growing in one-fifth Hoagland solution were set as control. The solutions were renewed every 2 days. After treatments for 1 h, 3 h, 6 h, 1 d, 3 d, and 6 d, roots of barley seedlings under both control and abiotic stress conditions were sampled for RNA extraction and qRT-PCR. All samples were collected in three replicates.
Total RNA was extracted using MiniBEST Plant RNA Extraction Kit (9769, TaKaRa, Japan) following the manufacturer’s instructions. The cDNA was synthesized from total RNA (1 μg) using PrimeScript RT Master Mix (RR036A, TaKaRa, Japan) and was used as templates for qRT-PCR amplification. qRT-PCR amplification was performed with LightCycler 480 II (Roche, Basel, Switzerland) using iTaq Universal SYBR Green Supermix (1,725,124, Bio-Rad, USA). The relative gene expression was calculated based on the 2−△△CT method using actin as the internal standard . The primer sequences were listed in the Additional file 7.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
Auxin response factor 2
Monovalent cation/proton antiporter
High-affinity K+ transporter
High-affinity K+ transporter
K+ efflux antiporter
K+ uptake permease
Non-selective cation channel
Short hypocotyl 3
Tiny root hair 1
Anschütz U, Becker D, Shabala S. Going beyond nutrition: regulation of potassium homoeostasis as a common denominator of plant adaptive responses to environment. J Plant Physiol. 2014;171(9):670–87. https://doi.org/10.1016/j.jplph.2014.01.009.
Hasanuzzaman M, Bhuyan MHMB, Nahar K, Hossain MS, Al Mahmud J, Hossen MS, et al. Potassium: a vital regulator of plant responses and tolerance to abiotic stresses. Agronomy. 2018;8(3). https://doi.org/10.3390/agronomy8030031.
Fageria NK. The use of nutrients in crop plants. 1st ed. Boca Raton: CRC Press; 2009.
Nieves-Cordones M, Alemán F, Martínez V, Rubio F. K+ uptake in plant roots. The systems involved, their regulation and parallels in other organisms. J Plant Physiol. 2014;171(9):688–95. https://doi.org/10.1016/j.jplph.2013.09.021.
Cellier F, Conéjéro G, Ricaud L, Doan TL, Lepetit M, Gosti F, et al. Characterization of AtCHX17, a member of the cation/H+ exchangers, CHX family, from Arabidopsis thaliana suggests a role in K+ homeostasis. Plant J. 2004;39(6):834–46. https://doi.org/10.1111/j.1365-313X.2004.02177.x.
Schleyer M, Bakker EP. Nucleotide sequence and 3′-end deletion studies indicate that the K(+)-uptake protein kup from Escherichia coli is composed of a hydrophobic core linked to a large and partially essential hydrophilic C terminus. J Bacteriol. 1993;175(21):6925–31. https://doi.org/10.1128/jb.175.21.6925-6931.1993.
Santa-María GE, Rubio F, Dubcovsky J, Rodríguez-Navarro A. The HAK1 gene of barley is a member of a large gene family and encodes a high-affinity potassium transporter. Plant Cell. 1997;9(12):2281–9. https://doi.org/10.1105/tpc.9.12.2281.
Yang T, Zhang S, Hu Y, Wu F, Hu Q, Chen G, et al. The role of a potassium transporter OsHAK5 in potassium acquisition and transport from roots to shoots in rice at low potassium supply levels. Plant Physiol. 2014;166(2):945–59. https://doi.org/10.1104/pp.114.246520.
Bañuelos MA, Garciadeblas B, Cubero B, Rodríguez-Navarro A. Inventory and functional characterization of the HAK potassium transporters of rice. Plant Physiol. 2002;130(2):784–95. https://doi.org/10.1104/pp.007781.
Gupta M, Qiu X, Wang L, Xie W, Zhang C, Xiong L, et al. KT/HAK/KUP potassium transporters gene family and their whole-life cycle expression profile in rice (Oryza sativa). Mol Gen Genomics. 2008;280(5):437–52. https://doi.org/10.1007/s00438-008-0377-7.
Cheng X, Liu X, Mao W, Zhang X, Chen S, Zhan K, et al. Genome-wide identification and analysis of HAK/KUP/KT potassium transporters gene family in wheat (Triticum aestivum L.). Int J Mol Sci. 2018;19(12). https://doi.org/10.3390/ijms19123969.
Zhang Z, Zhang J, Chen Y, Li R, Wang H, Wei J. Genome-wide analysis and identification of HAK potassium transporter gene family in maize (Zea mays L.). Mol Biol Rep. 2012;39(8):8465–73. https://doi.org/10.1007/s11033-012-1700-2.
Gierth M. The potassium transporter AtHAK5 functions in K+ deprivation-induced high-affinity K+ uptake and AKT1 K+ channel contribution to K+ uptake kinetics in Arabidopsis roots. Plant Physiol. 2005;137(3):1105–14. https://doi.org/10.1104/pp.104.057216.
Han M, Wu W, Wu WH, Wang Y. Potassium transporter KUP7 is involved in K+ acquisition and translocation in Arabidopsis root under K+-limited conditions. Mol Plant. 2016;9(3):437–46. https://doi.org/10.1016/j.molp.2016.01.012.
Chen G, Hu Q, Luo L, Yang T, Zhang S, Hu Y, et al. Rice potassium transporter OsHAK1 is essential for maintaining potassium-mediated growth and functions in salt tolerance over low and high potassium concentration ranges. Plant Cell Environ. 2015;38(12):2747–65. https://doi.org/10.1111/pce.12585.
Nieves-Cordones M, Miller AJ, Alemán F, Martínez V, Rubio F. A putative role for the plasma membrane potential in the control of the expression of the gene encoding the tomato high-affinity potassium transporter HAK5. Plant Mol Biol. 2008;68(6):521–32. https://doi.org/10.1007/s11103-008-9388-3.
Nieves-Cordones M, Alemán F, Martínez V, Rubio F. The Arabidopsis thaliana HAK5 K+ transporter is required for plant growth and K+ acquisition from low K+ solutions under saline conditions. Mol Plant. 2010;3(2):326–33. https://doi.org/10.1093/mp/ssp102.
Zhao S, Zhang ML, Ma TL, Wang Y. Phosphorylation of ARF2 relieves its repression of transcription of the K+ transporter gene HAK5 in response to low potassium stress. Plant Cell. 2016;28(12):3005–19. https://doi.org/10.1105/tpc.16.00684.
Rigas S, Ditengou FA, Ljung K, Daras G, Tietz O, Palme K, et al. Root gravitropism and root hair development constitute coupled developmental responses regulated by auxin homeostasis in the Arabidopsis root apex. New Phytol. 2013;197(4):1130–41. https://doi.org/10.1111/nph.12092.
Rigas S, Debrosses G, Haralampidis K, Vicente-Agullo F, Feldmann KA, Grabov A, et al. TRH1 encodes a potassium transporter required for tip growth in Arabidopsis root hairs. Plant Cell. 2001;13(1):139–51. https://doi.org/10.1105/tpc.13.1.139.
Cai K, Chen X, Han Z, Wu X, Zhang S, Li Q, et al. Screening of worldwide barley collection for drought tolerance: the assessment of various physiological measures as the selection criteria. Front Plant Sci. 2020;11:1159. https://doi.org/10.3389/fpls.2020.01159.
Yang Z, Gao Q, Sun C, Li W, Gu S, Xu C. Molecular evolution and functional divergence of HAK potassium transporter gene family in rice (Oryza sativa L.). J Genet Genom. 2009;36(3):161–72. https://doi.org/10.1016/S1673-8527(08)60103-4.
Ou W, Mao X, Huang C, Tie W, Yan Y, Ding Z, et al. Genome-wide identification and expression analysis of the KUP family under abiotic stress in cassava (Manihot esculenta Crantz). Front Physiol. 2018;9:17. https://doi.org/10.3389/fphys.2018.00017.
Song ZZ, Ma RJ, Yu ML. Genome-wide analysis and identification of KT/HAK/KUP potassium transporter gene family in peach (Prunus persica). Genet Mol Res. 2015;14(1):774–87. https://doi.org/10.4238/2015.January.30.21.
Li Y, Peng L, Xie C, Shi X, Dong C, Shen Q, et al. Genome-wide identification, characterization, and expression analyses of the HAK/KUP/KT potassium transporter gene family reveals their involvement in K+ deficient and abiotic stress responses in pear rootstock seedlings. Plant Growth Regul. 2018;85(2):187–98. https://doi.org/10.1007/s10725-018-0382-8.
Liang M, Liang M, Gao Y, Gao Y, Mao T, Mao T, et al. Characterization and expression of KT/HAK/KUP transporter family genes in willow under potassium deficiency, drought, and salt stresses. Biomed Res Int. 2020;2020:2690760.
Mascher M, Gundlach H, Himmelbach A, Beier S, Twardziok SO, Wicker T, et al. A chromosome conformation capture ordered sequence of the barley genome. Nature. 2017;544(7651):427–33. https://doi.org/10.1038/nature22043.
Nieves-Cordones M, Ródenas R, Chavanieu A, Rivero RM, Martinez V, Gaillard I, et al. Uneven HAK/KUP/KT protein diversity among angiosperms: species distribution and perspectives. Front Plant Sci. 2016;7:127.
Cuin TA, Bose J, Stefano G, Jha D, Tester M, Mancuso S, et al. Assessing the role of root plasma membrane and tonoplast Na+/H+ exchangers in salinity tolerance in wheat: in planta quantification methods. Plant Cell Environ. 2011;34(6):947–61. https://doi.org/10.1111/j.1365-3040.2011.02296.x.
Shabala S, Demidchik V, Shabala L, Cuin TA, Smith SJ, Miller AJ, et al. Extracellular Ca2+ ameliorates NaCl-induced K+ loss from Arabidopsis root and leaf cells by controlling plasma membrane K+-permeable channels. Plant Physiol. 2006;141(4):1653–65. https://doi.org/10.1104/pp.106.082388.
Cai K, Gao H, Wu X, Zhang S, Han Z, Chen X, et al. The ability to regulate transmembrane potassium transport in root is critical for drought tolerance in barley. Int J Mol Sci. 2019;20(17):4111. https://doi.org/10.3390/ijms20174111.
Savojardo C, Martelli PL, Fariselli P, Profiti G, Casadio R. BUSCA: an integrative web server to predict subcellular localization of proteins. Nucleic Acids Res. 2018;46(W1):W459–66. https://doi.org/10.1093/nar/gky320.
Madeira F, Park YM, Lee J, Buso N, Gur T, Madhusoodanan N, et al. The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 2019;47(W1):W636–41. https://doi.org/10.1093/nar/gkz268.
Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59(3):307–21. https://doi.org/10.1093/sysbio/syq010.
Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202. https://doi.org/10.1016/j.molp.2020.06.009.
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, et al. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37(Web Server):W202–8. https://doi.org/10.1093/nar/gkp335.
Shen Q, Fu L, Dai F, Jiang L, Zhang G, Wu D. Multi-omics analysis reveals molecular mechanisms of shoot adaption to salt stress in Tibetan wild barley. BMC Genomics. 2016;17(1):1–15. https://doi.org/10.1186/s12864-016-3242-9.
Ye Z, Zeng J, Li X, Zeng F, Zhang G. Physiological characterizations of three barley genotypes in response to low potassium stress. Acta Physiol Plant. 2017;39(1):1–10. https://doi.org/10.1007/s11738-016-2300-x.
We thank PhD candidate Yuanyuan Wang from Zhejiang University for her comments and suggestions on phylogenetic analyses.
This study was supported by National Natural Science Foundation of China (31620103912, 31671678), China Agriculture Research System (CARS-05) and Jiangsu Collaborative Innovation Centre for Modern Crop Production (JCIC-MCP). The funders had no roles in study design, data collection and analysis, and result interpretation and manuscript writing.
Authors and Affiliations
Institute of Crop and Nuclear Technology Utilization, Zhejiang Academy of Agricultural Sciences, Hangzhou, 310021, China
Kangfeng Cai & Junmei Wang
Institute of Crop Science, Zhejiang University, Hangzhou, 310058, China
KC and GZ conceived the project. KC designed and conducted the experiments. KC performed the main data analysis. KC, FZ, JW and GZ contributed to result interpretation. KC wrote the manuscript. GZ revised the manuscript. All authors read and approved the final manuscript.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
Cai, K., Zeng, F., Wang, J. et al. Identification and characterization of HAK/KUP/KT potassium transporter gene family in barley and their expression under abiotic stress.
BMC Genomics22, 317 (2021). https://doi.org/10.1186/s12864-021-07633-y