Genome-wide characterization and expression analysis of the HAK gene family in response to abiotic stresses in Medicago

The high-affinity K+ transporter (HAK) family plays a vital role in K+ uptake and transport as well as in salt and drought stress responses. In the present study, we identified 22 HAK genes in each Medicago truncatula and Medicago sativa genome. Phylogenetic analysis suggested that these HAK proteins could be divided into four clades, and the members of the same subgroup share similar gene structure and conserved motifs. Many cis-acting elements related with defense and stress were found in their promoter region. In addition, gene expression profiles analyzed with genechip and transcriptome data showed that these HAK genes exhibited distinct expression pattern in different tissues, and in response to salt and drought treatments. Furthermore, co-expression analysis showed that 6 homologous HAK hub gene pairs involved in direct network interactions. RT-qPCR verified that the expression level of six HAK gene pairs was induced by NaCl and mannitol treatment to different extents. In particular, MtHK2/7/12 from M. truncatula and MsHAK2/6/7 from M. sativa were highly induced. The expression level of MsHAK1/2/11 determined by RT-qPCR showed significantly positive correlation with transcriptome data. In conclusion, our study shows that HAK genes play a key role in response to various abiotic stresses in Medicago, and the highly inducible candidate HAK genes could be used for further functional studies and molecular breeding in Medicago.