浙江大学学报（农业与生命科学版）

浙江大学学报（农业与生命科学版）, 2021, 47(3): 303-313 doi: 10.3785/j.issn.1008-9209.2020.08.251

作物栽培与生理

甘蓝型油菜钾离子转运载体HAK/KUP/KT家族的全基因组鉴定与分析

朱乐,,, 赵鑫泽, 蒋立希,,

浙江大学农业与生物技术学院作物科学研究所/浙江省作物种质资源重点实验室，杭州 310058

Genome-wide identification and analysis of potassium ion transporter HAK/KUP/KT family in rapeseed (Brassica napus L.)

ZHU Le,,, ZHAO Xinze, JIANG Lixi,,

Institute of Crop Science, College of Agriculture and Biotechnology, Zhejiang University/Zhejiang Key Laboratory of Crop Gene Resources, Hangzhou 310058, China

通讯作者: 蒋立希（https://orcid.org/0000-0002-8579-0763），Tel：+86-571-88982905，E-mail：jianglx@zju.edu.cn蒋立希（https://orcid.org/0000-0002-8579-0763），Tel：+86-571-88982905，E-mail：jianglx@zju.edu.cn

收稿日期: 2020-08-25   接受日期: 2020-12-28   网络出版日期: 2021-07-05

基金资助: 国家自然科学基金.  31961143008

Received: 2020-08-25   Accepted: 2020-12-28   Online: 2021-07-05

作者简介 About authors

朱乐（https://orcid.org/0000-0001-8410-3717），E-mail：11616028@zju.edu.cn , E-mail：11616028@zju.edu.cn

摘要

钾离子转运载体HAK/KUP/KT（high-affinity K＋/K＋ uptake permease/K＋ transporter）家族成员在植物对钾（K）的吸收和转运以及植物生长发育、耐盐性和渗透势的调节中起着重要作用。本研究以甘蓝型油菜（Brassica napus L.）为材料，共鉴定出45个HAK/KUP/KT基因，分别命名为BnHAK1～BnHAK45，并对油菜中的HAK/KUP/KT基因家族进行基因结构、进化和表达以及染色体定位分析。结果表明，甘蓝型油菜中的HAK/KUP/KT基因家族在染色体上呈现不均匀分布，其中在A01号染色体上分布最多，有5个BnHAKs基因。根据不同物种HAK/KUP/KT蛋白构建系统进化树，将其分为2个组群，其中双子叶植物中的拟南芥、甘蓝、白菜和甘蓝型油菜基因的保守性更强，亲缘关系更近。在表达量分析中，尽管HAK/KUP/KT基因家族的表达模式各不相同，但是总体而言，它们在根中的表达量最高，在角果中的表达量最低，这说明不同HAK/KUP/KT基因家族成员的分工是不同的。BnHAK4、BnHAK9基因在叶中特异性高表达，BnHAK11、BnHAK18和BnHAK23在根中特异性高表达，BnHAK7和BnHAK17在花中的表达量最高，说明该基因家族的表达具有组织特异性。本研究为进一步揭示甘蓝型油菜HAK/KUP/KT基因家族的功能机制提供了相应的参考依据。

关键词： 甘蓝型油菜 ; 钾离子转运载体HAK/KUP/KT ; 基因家族分析

Abstract

The high-affinity K＋/K＋ uptake permease/K＋ transporter (HAK/KUP/KT) family members play important roles in the absorption and transport of potassium (K), regulations of plant growth and development, modulation/acquisition of salt tolerance, and osmotic potential regulation/control. However, this gene family has not been thoroughly studied in Brassica napus. In this study, forty-five BnHAK genes (BnHAK1-BnHAK45) were identified in the B. napus reference genome, and bioinformatic methods were used to resolve the gene structures, evolutionary analysis, expression patterns and chromosome position in different tissues of rape. The results showed that the members of the HAK/KUP/KT gene family were unevenly distributed on the chromosomes, with the most of them located on chromosome A01 which had five BnHAK genes. Phylogenetic analysis of HAK/KUP/KT proteins from multiple species showed that they could be divided into two groups. Among them, Arabidopsis, B. oleracea, B. rapa and B. napus had stronger gene conservation and closer genetic relationships. Although the expression patterns of the HAK/KUP/KT genes were different, in general, they exhibited highest and lowest expression levels in the roots and siliques, respectively, and this showed that there was a functional division of different HAK/KUP/KT gene family members. The expression levels of BnHAK4 and BnHAK9 were highest in the leaves, whilst the expression levels of BnHAK11, BnHAK18, and BnHAK23 were highest in the roots, and BnHAK7 and BnHAK17 were mostly expressed in the flowers, which indicated that the expression of this gene family had tissue specificity. This study provides a reference for further revealing the functional mechanism of HAK/KUP/KT gene family in B. napus.

Keywords： Brassica napus ; potassium ion transporter ; HAK/KUP/KT ; gene family analysis

PDF (3485KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

朱乐, 赵鑫泽, 蒋立希. 甘蓝型油菜钾离子转运载体HAK/KUP/KT家族的全基因组鉴定与分析. 浙江大学学报（农业与生命科学版）[J]. 2021, 47(3): 303-313 doi:10.3785/j.issn.1008-9209.2020.08.251

ZHU Le, ZHAO Xinze, JIANG Lixi. Genome-wide identification and analysis of potassium ion transporter HAK/KUP/KT family in rapeseed (Brassica napus L.). Journal of Zhejiang University（Agriculture & Life Sciences）[J]. 2021, 47(3): 303-313 doi:10.3785/j.issn.1008-9209.2020.08.251

钾（K）是植物内一种重要的大量营养元素，占植物干物质量的10%。钾参与多种生化过程，包括蛋白质合成、碳水化合物代谢和酶激活，以及许多生理过程，如气孔调节和光合作用[1]。此外，钾还可以增强作物对各种生物和非生物胁迫的抵抗力[2-3]。

植物对K＋的需求量很大而植物所处的环境又非常复杂，使得K＋转运系统十分庞大，以至于植物从土壤中吸收的K＋可以通过不同的组织进行再分配[4-5]。几乎在所有的植物组织中，K＋转运系统都是一个大的基因家族，编码多种功能的蛋白质，K＋转运载体家族基因通常被分为4大家族：HAK/KUP/KT、Trk/HKT、KEA和CHX家族。通过与细菌KUP和真菌HAKs的序列进行相似性比较，在拟南芥中鉴定出第1个植物HAK/KUP/KT基因[6-7]。目前，HAK/KUP/KT家族数目庞大，对植物的生长发育调节起着重要作用。植物根系对钾的吸收及其在植物体内的转运是由许多K＋通道和转运子共同完成的[8-9]。K＋通道在外部高K＋浓度（＞0.5 mmol/L）下起作用，被认为是低亲和转运系统；K＋转运体在外部低K＋浓度（＜0.2 mmol/L）下也可以发挥作用，属于高亲和转运系统；但是，随着K＋通道和转运体数目的增加，这2类系统之间的分界线变得越来越模糊。

植物中存在多种HAK/KUP/KT转运蛋白，它们在K的吸收和转运以及植物生长发育、耐盐性和渗透势的调节中发挥着不同的作用[10]。在拟南芥中，K＋转运蛋白AtHAK5和K＋通道蛋白AtAKT1主要在根细胞中表达，促进根细胞吸收外界环境中的K＋[11]。在正常K＋浓度条件下，在拟南芥根系中可检测到AtHAK5转录本，而在芽中检测不到；但在K＋缺乏的情况下，AtHAK5转录本在根系中的表达水平升高，而且在芽中也能检测到AtHAK5转录本[12]。在水稻中，AtHAK5和AtAKT1的功能是由其同源基因OsAKT5和OsHAK1实现的。研究发现，水稻木质部汁液中的K＋浓度以及OsHAK1突变体的K＋净输出率明显低于野生型植物，特别是低K＋处理（0.1 mmol/L）的植株[13]。K＋转运蛋白OsHAK5可提高K＋吸收亲和力，在K＋浓度高于OsHAK1蛋白浓度的情况下，OsHAK5基因在根部维管组织的木质部薄壁组织和韧皮部中强烈表达，在缺K条件下尤为明显，表明OsHAK5可能参与了根与芽之间的K分配[14]。OsHAK21具有K＋转运活性，但不直接参与K＋吸收。OsHAK21主要在中柱的根维管系统中表达，可能在根与芽的长距离K＋分配中起作用[15]。

植物体内K＋/Na＋比值的增大可增强其耐盐性。OsHAK5的过表达增大了植株的K＋/Na＋比值，从而提高了植株的耐盐性。相反，OsHAK5基因敲除突变体的K＋/Na＋比值降低，结果转基因植株对盐胁迫的敏感性增加[14]。OsHAK21在盐胁迫反应中也发挥着重要作用。在高盐胁迫条件下，OsHAK21的表达量显著增加，特别是在盐胁迫4 h后，OsHAK21在根中的表达水平增加了600倍。在高盐胁迫下，OsHAK21突变体在地上部和根部积累的K＋量均低于野生型（对照）。此外，HAK/KUP/KT转运蛋白也在抗旱性中发挥作用。OsHAK1在水稻营养期和生殖期的过表达均增强了水稻的抗旱性，与对照相比，OsHAK1过表达幼苗的脂质过氧化水平降低，脯氨酸积累量增加，抗氧化酶活性提高[15]。在干旱条件下，OsHAK1过表达植株的产量比野生型植株高出35%[16]。

上述研究表明，HAK/KUP/KT家族成员对作物的生长发育具有潜在价值，可以提高作物对干旱和盐胁迫的耐受力。同时，随着我国缺K土壤面积日益扩大，科学合理地施用K肥是提高K肥利用效率的有效措施[17-18]，而筛选K高效基因型作物品种是缓解K肥需求压力的另一有效途径。近年来，越来越多的植物HAK/KUP/KT家族被研究，如梨[19]、甘蔗[20]、辣椒[21]、小麦[22]、大豆[23]和水稻[24]等。油菜是世界上最重要的油料作物之一，其中甘蓝型冬油菜生育期长，对K素需求量大，但是目前关于甘蓝型油菜（Brassica napus L.）HAK/KUP/KT基因家族系统的鉴定和分析研究未见报道。

本研究利用已经发布的甘蓝型油菜基因组数据，对油菜中HAK/KUP/KT类K转运蛋白家族进行系统鉴定，并对HAK/KUP/KT家族成员的基因特征、系统发育树、蛋白质结构域、基因结构、染色体位置和顺式作用元件等进行了分析，为深入了解HAK/KUP/KT基因家族的功能及培育K高效油菜品种提供理论依据。

1 材料与方法

1.1 甘蓝型油菜HAK/KUP/KT家族成员的鉴定与分析

甘蓝型油菜的蛋白质序列来源于已公布的数据库（http://www.genoscope.cns.fr/brassicanapus/）；分别从拟南芥的数据库[The Arabidopsis Information Resource (TAIR), https://www.arabidopsis.org/]和水稻的数据库（Rice Genome Annotation Project, http://rice.plantbiology.msu.edu）中下载拟南芥和水稻的蛋白质序列；从Phytozome数据库（https://phytozome.jgi.doe.gov/pz/portal.html）中下载甘蓝、白菜和大麦的蛋白质序列；随后，以其他物种的HAK/KUP/KT蛋白质序列为查询序列，与甘蓝型油菜的蛋白质序列进行Blastp比对，筛选阈值小于10－10、相似性大于50%的蛋白质序列；再在Pfam数据库（http://pfam.xfam.org/）中下载HAK/KUP/KT基因家族的隐马尔可夫模型（PF02705）。利用在线分析软件HMMER（https://www.ebi.ac.uk/Tools/hmmer/）和SMART（http://smart.embl-heidelberg.de/）分别对蛋白质序列的保守结构域进行鉴定，最终确定甘蓝型油菜特异的HAK/KUP/KT家族成员。通过ExPASy网站（https://www.expasy.org/）的在线分析软件对蛋白质的分子质量、等电点等特征信息进行分析。

1.2 甘蓝型油菜HAK/KUP/KT蛋白家族的系统分析和进化树构建

对甘蓝型油菜中HAK/KUP/KT蛋白家族进行系统进化分析。将从其他物种和甘蓝型油菜中鉴定出来的HAK/KUP/KT蛋白序列利用Clustal X进行多序列比对，并将比对结果利用MEGA X软件进行格式转换对齐，然后利用邻接法绘制系统发育树，其中自展值（bootstrap）设置为1 000[25]。

1.3 甘蓝型油菜HAK/KUP/KT基因的染色体定位与分析

根据甘蓝型油菜HAK/KUP/KT基因序列信息，获得每个基因在油菜染色体上的位置信息，并利用MapChart软件（https://www.wur.nl/en/show/Mapchart.html）进行染色体定位和作图[26]。

1.4 甘蓝型油菜HAK/KUP/KT蛋白序列的基序分析

利用MEME在线工具（http://meme-suite.org/)对甘蓝型油菜的HAK/KUP/KT蛋白序列进行基序（motif）分析，基序个数设置为10，其他参数为默认数值。

1.5 甘蓝型油菜HAK/KUP/KT基因结构分析

从甘蓝型油菜基因组数据库中获取BnHAKs基因的编码序列（coding sequence, CDS）和非编码区（untranslated region, UTR）信息，利用GSBS（http://gsds.cbi.pku.edu.cn/index.php）和TBtools（https://github.com/CJ-Chen/TBtools）对基因的内含子和外显子进行可视化展示。

1.6 甘蓝型油菜HAK/KUP/KT基因顺式作用元件分析

从甘蓝型油菜基因组数据库中提取BnHAKs基因的启动子序列（转录起始位置为上游1 500 bp序列处）。顺式作用元件的预测通过PLANTCARE在线网站（http://bioinformatics.psb.ugent.be/webtools/plantcare/html/）进行。

1.7 甘蓝型油菜HAK/KUP/KT基因表达模式分析

从公开的RNA-seq数据库[27]下载甘蓝型油菜根、茎、叶、萼片、果实和花等6个组织的表达谱数据文件，进行标准均一化处理后，使用MeV 4.9软件绘制基因表达热图。

1.8 实时荧光定量聚合酶链式反应（real‑time quantitative polymerase chain reaction, RT-qPCR）分析

利用快速通用植物RNA提取试剂盒[华越洋生物（北京）科技有限公司]提取样品RNA。第1链cDNA利用GoScriptTM反转录系统（GoScriptTM Reverse Transcription Mix）、Oligo （dT）[普洛麦格（北京）生物技术有限公司]合成。20 μL RT-qPCR扩增体系：2 μL cDNA模板，10 μmol/L上、下游引物各 1 μL（引物信息如表1所示）、10 μL Promega Eastep qPCR主混合物和6 μL ddH2O。PCR程序：50 ℃预变性600 s；95 ℃变性180 s；95 ℃变性15 s，60 ℃退火和延伸60 s，共45个循环；从60 ℃升温到95 ℃，进行熔解曲线分析。每个样品设置3个生物学重复，以油菜ACTIN7为内参基因，各基因相对表达量使用2－△△CT法计算得出，并且使用软件Excel 2016进行数据统计和分析。

表1   RT-qPCR引物序列

Table 1  Primer sequences for RT-qPCR

引物名称

Primer name

引物序列（5´→3´）

Primer sequence (5´→3´)

HAK1-FCAACAATCCAAAGATGAGCAGCGAC
HAK1-RGCTCCGATGAAGAGGAACCAAATG
HAK4-FCCTCCTCTTAGCTTACATGGGCCA
HAK4-RAACCTTGGGAAGCAGCCGAGA
HAK9-FGATGGGATTCTCACACCAGCTATGT
HAK9-RAGACGATGTAATACGGGTTCAGTGC
HAK13-FAGCCTCTTCCTTTTCATCCGCC
HAK13-RCTGGCTCGACCCAGTCACCG
HAK18-FGAAGCATGTGGGATTTGGACCA
HAK18-RGAGAGCTCCAATGATGTCTTCAGGG
HAK37-FAACAATCCAAAGATGAGCAGCGAC
HAK37-RGATAGGTGCAAAGAGCCAGCCTACT
ACTIN7-FGCTGACCGTATGAGCAAAG
ACTIN7-RAAGATGGATGGACCCGAC

新窗口打开| 下载CSV


2 结果与分析

2.1 甘蓝型油菜HAK/KUP/KT家族成员的鉴定与分析结果

使用Blastp方法对已报道的拟南芥、水稻与甘蓝型油菜的HAK/KUP/KT蛋白序列进行比对，将所得的序列比对Pfam数据库，初步确定了45个编码HAK/KUP/KT蛋白的基因。根据这些基因在染色体上的分布位置，分别将其命名为BnHAK1～BnHAK45（附表1，http://www.zjujournals.com/agr/CN/10.3785/j.issn.1008-9209.2020.08.251）。对鉴定到的45个BnHAKs蛋白的理化性质进行分析，发现甘蓝型油菜中BnHAKs蛋白的氨基酸长度在204～1 526个氨基酸残基之间，蛋白质分子质量为22.61～171.55 kDa，等电点范围为4.88～9.63，大部分BnHAKs编码碱性蛋白。

2.2 甘蓝型油菜HAK/KUP/KT蛋白家族的系统进化分析

为充分研究甘蓝型油菜HAK/KUP/KT基因家族的进化关系，分别利用甘蓝（23个）、白菜（19个）、拟南芥（13个）、水稻（27个）和大麦（25个）的HAK/KUP/KT蛋白序列与甘蓝型油菜的45条BnHAKs蛋白序列构建系统发育树。

结果（图1）表明，HAK/KUP/KT家族分成2大类，双子叶植物油菜HAK/KUP/KT家族成员与双子叶植物拟南芥、甘蓝和白菜的亲缘关系较近，而与单子叶植物水稻和大麦的亲缘关系较远。在同一亚家族内的蛋白序列保守度较高，而不同亚家族间的差异较大，在同一亚家族内、不同亚组之间也存在一定差异。亚家族Ⅱb的甘蓝和单子叶植物水稻聚集在一起，说明该基因家族出现得较早，随着物种进化，功能比较保守；甘蓝型油菜（B. napus L.，2n=AACC=38）是由甘蓝（B. oleracea L.，2n=CC=18）和白菜（B. rapa L.，2n=AA=20）通过天然种间杂交而形成的异源四倍体。甘蓝型油菜和甘蓝处于同一姊妹分支的比例高于白菜，表明甘蓝型油菜HAK/KUP/KT家族成员与甘蓝的亲缘关系更近，HAK/KUP/KT家族基因可能更多的来自A染色体，与前人对甘蓝型油菜的全基因组分析结果[28]一致。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   甘蓝型油菜与其他植物HAK/KUP/KT蛋白的系统进化分析

Fig. 1   Phylogenetic analysis of HAK/KUP/KT proteins in B. napus and other plants


2.3 甘蓝型油菜HAK/KUP/KT基因染色体定位结果

从HAK/KUP/KT在甘蓝型油菜基因组上的分布来看，45个BnHAKs基因分布在油菜的1 922条染色体上（图2）。其中：A01号染色体上分布得最多，有5个BnHAKs基因；A03、C04和C07这3条染色体上分布的其次，均有4个BnHAKs基因；A02、A05、A08、C01、C02、C03、C06、C09、Cnn_random和C01_random这10条染色体上均分布2个BnHAKs基因；其余的8条染色体上均分布1个。本研究中的45个甘蓝型油菜HAK/KUP/KT基因没有形成基因簇分布，说明甘蓝型油菜HAK/KUP/KT基因家族无法通过串联重复来实现扩增，因此，HAK/KUP/KT基因家族的扩增可能主要为分散复制和片段复制的形式。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   甘蓝型油菜HAK/KUP/KT基因家族的染色体定位

Fig. 2   Chromosome location of HAK/KUP/KT gene family in B. napus


2.4 甘蓝型油菜HAK/KUP/KT基因家族的结构域分析

根据系统发育分析，可以将甘蓝型油菜HAK/KUP/KT基因家族分成4个亚家族（图3A）。分析表明：BnHAKs基因序列中至少含有4个外显子和3个内含子，至多含有20个外显子和19个内含子。BnHAK10和BnHAK23都只含有4个外显子和3个内含子，是所有已鉴定的BnHAKs基因中含外显子和内含子最少的基因；BnHAK42含有20个外显子和19个内含子，是所有已鉴定的BnHAKs基因中外显子和内含子数量最多的基因。甘蓝型油菜HAK/KUP/KT基因家族含有9和10个外显子的基因各有10个（22.2%），含有8和9个内含子的基因数量最多。在同一个组群中，外显子和内含子的数目也不相同（图3B）。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   甘蓝型油菜HAK/KUP/KT基因的系统发育、基因结构和蛋白结构域分析

A. BnHAKs蛋白系统发育树分析；B. BnHAKs基因结构分析；C. BnHAKs蛋白的基序分析；D. 10种基序中的氨基酸保守性分析。UTR：非编码区；CDS：编码序列。

Fig. 3   Phylogenetic analysis and structure of HAK/KUP/KT gene and its architecture of conserved protein motifs from B. napus

A. Phylogenetic tree of BnHAKs; B. Exon/intron organization of BnHAKs; C. Motif composition of BnHAKs; D. Conserve amino acids in ten motifs. UTR: Untranslated region; CDS: Coding sequence.


对每个组群进行详细分析发现，在组群Ⅰ的2个亚组群中，只有BnHAK23和BnHAK19蛋白基序缺少较多；组群Ⅰ的基序与其他3组基序整齐度一致；组群Ⅱ的4个亚组群中BnHAK3、BnHAK5、BnHAK10、BnHAK20、BnHAK23、BnHAK38和BnHAK43蛋白基序缺少较多，特别是Ⅱb亚组群的BnHAK3蛋白基序只有2个，保守性较差。甘蓝型油菜BnHAKs蛋白结构域基序整体保守性较好，只有个别蛋白结构域的基序保守性差（图3C）。

通过MEME在线工具对甘蓝型油菜BnHAKs蛋白序列进行分析，共鉴定出10个保守基序（图3D），大多数的保守基序位于K＋超级家族结构域的序列中。总的来说，这些基序几乎均匀分布，在所有组群中，基序1、2、3、4、5、6、8、9和10都是保守的，只有1个基因因缺少特定基序而例外。

2.5 甘蓝型油菜HAK/KUP/KT基因家族顺式作用元件分析结果

为了进一步研究HAK/KUP/KT基因家族在非生物胁迫反应中的调控机制，对45条甘蓝型油菜HAK/KUP/KT家族基因的上游1 500 bp序列进行顺式作用元件分析。结果（图4）显示了11种非生物胁迫响应元件，包括低温应答元件（LTR）、赤霉素应答元件（TATC-box）、防御和胁迫应答元件（TC-rich repeats）、脱落酸应答元件（ABRE）、干旱应答元件（MBS）、茉莉酸甲酯（MeJA）应答元件（CGTCA-motif和TGACG-motif）、水杨酸应答元件（SARE和TCA-element）以及生长素应答元件（AuxRE和TGA-element）。上述结果表明，多数顺式元件参与激素调控和抗性相关的反应，所以我们推测HAK/KUP/KT基因家族在调控植物的抗性和生长发育的过程中发挥着重要作用。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   甘蓝型油菜HAK/KUP/KT基因家族顺势作用元件分析

Fig. 4   Analysis of the Cis-elements of HAK/KUP/KT gene family in B. napus


2.6 甘蓝型油菜HAK/KUP/KT基因家族在不同组织中的表达分析

利用公开的RNA-seq数据库预测HAK/KUP/KT基因家族在甘蓝型油菜根、茎、叶、萼片、果实和花中的表达量并绘制热图。在被检测的45个基因中，42个基因的数据来自甘蓝型油菜已发表的RNA序列数据库，而其他3个基因（BnHAK20、BnHAK26和BnHAK38）的信息丢失，其在6个组织中的表达量都设置为0。结果发现，42个BnHAKs基因的表达模式各不相同，总体而言，它们在根中的表达量最高，在角果中的表达量最低（图5A和附表1，http://www.zjujournals.com/agr/CN/10.3785/j.issn.1008-9209.2020.08.251）。

在组群Ⅰ的2个亚组群中，除BnHAK19、BnHAK24和BnHAK34之外，其余基因在6个组织中均有表达；组群Ⅱ的4个亚组群的大部分基因在不同组织中都有较高的表达水平，其中BnHAK41、BnHAK7和BnHAK8都在花中有较高的表达丰度，而BnHAK14只在根和茎中高表达，在其他组织中的表达量较低。

为了验证基因表达趋势是否正确，利用RT-qPCR分析了油菜根和叶中的基因表达情况，包括BnHAK1、BnHAK4、BnHAK9、BnHAK13、BnHAK18和BnHAK37。以甘蓝型油菜‘中双11’的第3片真叶的RNA为模板，基因表达结果如图5B所示：HAK/KUP/KT基因家族的表达具有组织特异性，其中BnHAK1在叶中的表达量远高于在根中的表达量，而BnHAK37在根中的表达丰度特别高，与转录表达结果一致。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   甘蓝型油菜HAK/KUP/KT基因家族在不同组织中的表达量

A. 转录表达分析；B. 荧光定量PCR表达分析。

Fig. 5   Expression levels of HAK/KUP/KT gene family in different tissues of B. napus

A. Transcriptional expression analysis; B. RT-qPCR expression analysis.


3 讨论

油菜是我国重要的油料作物，2017年油菜的种植面积约665.3万hm2，总产量132.7亿kg，占全国油料作物年产总量的38.2%，是当前我国食用植物油的主要来源。我国油菜种植面积和总产均居世界第一[29]，其产量和品质受K元素的影响较大。K是植物生存的必需元素，而K＋必须从细胞外摄取以满足植物自身的生理需要，因此，植物对K的吸收及转运也一直是植物生理研究的重点。

本研究对甘蓝型油菜基因组测序，共鉴定出45个HAK/KUP/KT基因，分别命名为BnHAK1～BnHAK45。对甘蓝型油菜和其他物种的HAK/KUP/KT的系统进化分析发现，HAK/KUP/KT家族可以分为2类，其中：第Ⅰ类包含的成员最少；第Ⅱ类成员最多，可分为4个亚组，所有物种均有涉及。尽管甘蓝型油菜是由甘蓝和白菜通过天然种间杂交而形成的异源四倍体，但是甘蓝型油菜和甘蓝在同一姊妹分支的比例高于白菜，表明甘蓝型油菜HAK/KUP/KT基因家族成员与甘蓝的亲缘关系更近。大多数的BnHAKs旁系同源基因具有相似的基序（motif），表明它们之间的功能具有一定的保守性；也有部分旁系同源基因具有不同的基序，这部分基因可能在进化过程中产生了功能分化。利用MEME分析共鉴定出10个基序，其中大部分被预测为K＋转运域的组成部分。激素信号传导途径及逆境胁迫响应调控相关的顺式元件在HAK/KUP/KT基因启动子区域大量存在，暗示着HAK/KUP/KT基因家族在多种激素信号转导途径与逆境胁迫响应的调控过程中可能发挥重要作用。

表达模式对基因功能的研究具有重要作用，BnHAKs在不同组织中的表达模式具有显著差异。对甘蓝型油菜6个组织中的HAK/KUP/KT基因家族成员分析发现，这45个成员具有不同的时空表达特异性，如BnHAK4、BnHAK9在叶中特异性高表达，显示出二者对叶发育的重要性；BnHAK7和BnHAK17在花中的表达量最高，表明它们在花发育中具有重要作用；BnHAK11、BnHAK18和BnHAK23在根中特异性高表达，表明它们在根发育中具有重要作用。BnHAK2、BnHAK5、BnHAK12、BnHAK22、BnHAK32、BnHAK33、BnHAK36、BnHAK39和BnHAK45在大部分组织中的表达量均较高，而BnHAK3、BnHAK19、BnHAK24、BnHAK25、BnHAK29、BnHAK34和BnHAK43在大部分组织中的表达量均较低，这说明不同HAK/KUP/KT基因家族成员的分工是不同的。

4 结论

本研究鉴定出45个甘蓝型油菜HAK/KUP/KT基因，并对其蛋白质理化性质、染色体定位、基因结构和蛋白功能结构域进行了分析，预测了BnHAKs启动子区域响应多种激素和胁迫的调控元件，检测了HAK/KUP/KT基因在甘蓝型油菜6个组织中的表达模式，为进一步揭示甘蓝型油菜HAK/KUP/KT基因家族的功能机制奠定了基础。

参考文献

HASANUZZAMAN M, BHUYAN M H M B, NAHAR K, et al.

Potassium: a vital regulator of plant responses and tolerance to abiotic stresses

Agronomy, 2018,8(3):31. DOI:10.3390/agronomy8030031

[本文引用: 1]

AMRUTHA R N, SEKHAR P N, VARSHNEY R K, et al.

Genome-wide analysis and identification of genes related to potassium transporter families in rice (Oryza sativa L

.). Plant Science, 2007,172(4):708-721. DOI:10.1016/j.plantsci.2006.11.019

[本文引用: 1]

ASHLEY M K, GRANT M, GRABOY A.

Plant responses to potassium deficiencies: a role for potassium transport proteins

Journal of Experimental Botany, 2006,57(2):425-436. DOI:10.1093/jxb/erj034

[本文引用: 1]

WANG Y, WU W H.

Regulation of potassium transport and signaling in plants

Science Direct, 2015,37(4):33. DOI:10.1016/j.pbi.2017.06.006

[本文引用: 1]

WANG M, ZHENG Q S, SHEN Q R, et al.

The critical role of potassium in plant stress response

International Journal of Molecular Sciences, 2013,14(4):7370-7390. DOI:10.3390/ijms14047370

[本文引用: 1]

AHN S J, SHIN R, SCHACHTMAN D P.

Expression of KT/KUP genes in Arabidopsis and the role of root hairs in K＋ uptake

Plant Physiology, 2004,134(3):1135-1145. DOI:10.1104/pp.103.034660.

[本文引用: 1]

GUPTA M, QIU X H, WANG L, et al.

KT/HAK/KUP potassium transporters gene family and their whole-life cycle expression profile in rice (Oryza sativa)

Molecular Genetics and Genomics, 2008,280(5):437-452. DOI:10.1007/s00438-008-0377-7

[本文引用: 1]

VÉRY A A, NIEVES-CORDONES M, DALY M, et al.

Molecular biology of K＋ transport across the plant cell membrane: What do we learn from comparison between plant species?

Journal of Plant Physiology, 2014,171(9):748-769. DOI:10.1016/j.jplph.2014.01.011

[本文引用: 1]

WANG Y, WU W H.

Potassium transport and signaling in higher plants

Annual Review of Plant Biology, 2013,64(1):451-476. DOI:10.1146/annurev-arplant-050312-120153

[本文引用: 1]

李学文,游西龙,王艳.

钾离子转运载体HAK/KUP/KT家族参与植物耐盐性的研究进展

.植物科学学报,2019,37(1):101-108. DOI:10.11913/PSJ.2095-0837.2019.10101

[本文引用: 1]

LI X W, YOU X L, WANG Y.

Research progress of HAK/KUP/KT potassium transporter family in plant response to salt stress

Plant Science Journal, 2019,37(1):101-108. (in Chinese with English abstract)

DOI:10.11913/PSJ.2095-0837.2019.10101      [本文引用: 1]

NIEVES-CORDONES M, MARTÍNEZ V, BENITO B, et al.

Comparison between Arabidopsis and rice for main pathways of K＋ and Na＋ uptake by roots

Frontiers in Plant Science, 2016,7:992. DOI:10.3389/fpls.2016.00992

[本文引用: 1]

RUBIO F, SANTA-MARÍA G E, RODRÍGUEZ-NAVARRO A.

Cloning of Arabidopsis and barley cDNAs encoding HAK potassium transporters in root and shoot cells

Physiologia Plantarum, 2000,109(1):34-43. DOI:10.1034/j.1399-3054.2000.100106.x

[本文引用: 1]

CHEN G, HU Q D, LUO L, 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 & Environment, 2015,38(12):2747-2765. DOI:10.1111/pce.12585

[本文引用: 1]

YANG T Y, ZHANG S, HU Y B, 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 Physiology, 2014,166(2):945-959. DOI:10.1104/pp.114.246520

[本文引用: 2]

SHEN Y, SHEN L K, SHEN Z X, et al.

The potassium transporter OsHAK21 functions in the maintenance of ion homeostasis and tolerance to salt stress in rice

Plant, Cell & Environment, 2015, 38(12):2766-2779. DOI:10.1111/pce.12586

[本文引用: 2]

CHEN G, LIU C L, GAO Z Y, et al.

OsHAK1, a high-affinity potassium transporter, positively regulates responses to drought stress in rice

Frontiers in Plant Science, 2017,8:1885. DOI:10.3389/fpls.2017.01885

[本文引用: 1]

高祥照,马文奇,崔勇,等.

我国耕地土壤养分变化与肥料投入状况

.植物营养与肥料学报,2000,6(4):363-369.

[本文引用: 1]

GAO X Z, MA W Q, CUI Y, et al.

Changes of soil nutrient contents and input of nutrients in arable of China

Plant Nutrition and Fertilizers Science, 2000,6(4):363-369. (in Chinese with English abstract)

[本文引用: 1]

易九红,刘爱玉.

作物钾效率基因型差异及缺钾反应

.作物研究,2007,21(5):536-540. DOI:10.16848/j.cnki.issn.1001-5280.2007.s1.011

[本文引用: 1]

YI J H, LIU A Y.

Genotype difference of potassium efficiency and potassium deficiency reaction in crops

Crop Research, 2007,21(5):536-540. (in Chinese)

DOI:10.16848/j.cnki.issn.1001-5280.2007.s1.011      [本文引用: 1]

WANG Y Z, LÜ J H, CHEN D, et al.

Genome-wide identification, evolution, and expression analysis of the KT/HAK/KUP family in pear

Genome, 2018,61(10):755-765. DOI:10.1139/gen-2017-0254

[本文引用: 1]

FENG X M, WANG Y J, ZHANG N N, et al.

Genome-wide systematic characterization of the HAK/KUP/KT gene family and its expression profile during plant growth and in response to low-K＋ stress in Saccharum

BMC Plant Biology, 2020,20:20. DOI:10.1186/s12870-019-2227-7

[本文引用: 1]

司伟娜,孙毅,李晓玉.

栽培型与野生型辣椒HAK/KUP/KT基因家族的进化分析

.安徽农业大学学报,2018,45(4):753-761. DOI:10.13610/j.cnki.1672-352x.20180825.021

[本文引用: 1]

SI W N, SUN Y, LI X Y.

Evolutionary analysis of HAK/KUP/KT gene family in cultivated and wild peppers

Journal of Anhui Agricultural University, 2018,45(4):753-761. (in Chinese with English abstract)

DOI:10.13610/j.cnki.1672-352x.20180825.021      [本文引用: 1]

CHENG X Y, LIU X D, MAO W W, et al.

Genome-wide identification and analysis of HAK/KUP/KT potassium transporters gene family in wheat (Triticum aestivum L

.). International Journal of Molecular Sciences, 2018,19(12):3969. DOI:10.3390/ijms19123969

[本文引用: 1]

晁毛妮,温青玉,张晋玉,等.

大豆KUP/HAK/KT钾转运体基因家族的鉴定与表达分析

.西北植物学报,2017,37(2):239-249. DOI:10.7606/j.issn.1000-4025.2017.02.0239

[本文引用: 1]

CHAO M N, WEN Q Y, ZHANG J Y, et al.

Identification and expression analysis of KUP/HAK/KT potassium transporter gene family in soybean [Glycine mac (L

.) Merr.]. Acta Botanica Boreali-Occidentalia Sinica, 2017,37(2):239-249. (in Chinese with English abstract)

DOI:10.7606/j.issn.1000-4025.2017.02.0239      [本文引用: 1]

OKADA T, YAMANE S, YAMAGUCHI M, et al.

Characterization of rice KT/HAK/KUP potassium trans-porters and K＋ uptake by HAK1 from Oryza sativa

Plant Biotechnology, 2018,35(2):101-111. DOI:10.5511/plantbiotechnology.18.0308a

[本文引用: 1]

LARKIN M A, BLACKSHIELDS G, BROWN N P, et al.

Clustal W and Clustal X version 2

.0. Bioinformatics, 2007,23(21):2947-2948. DOI:10.1093/bioinformatics/btm404

[本文引用: 1]

赵艳,瓮巧云,马海莲,等.

谷子ARF基因家族的鉴定与生物信息学分析

.植物遗传资源学报,2016,17(3):547-554. DOI:10.13430/j.cnki.jpgr.2016.03.022

[本文引用: 1]

ZHAO Y, WENG Q Y, MA H L, et al.

Genome-wide identification and bioinformatics analysis of ARF gene family in foxtail millet Setaria italica

Journal of Plant Genetic Resources, 2016,17(3):547-554. (in Chinese with English abstract)

DOI:10.13430/j.cnki.jpgr.2016.03.022      [本文引用: 1]

SUN F M, FAN G Y, HU Q, et al.

The high-quality genome of Brassica napus cultivar ‘ZS11’ reveals the introgression history in semi-winter morphotype

The Plant Journal, 2017,92(3):452-468. DOI:10.1111/tpj.13669

[本文引用: 1]

张毅,伍向苹,张芳,等.

基于基因组数据解析中国油菜品种演化历程及方向

.中国油料作物学报,2020,42(3):325-333. DOI:10.19802/j.issn.1007-9084.2019310

[本文引用: 1]

ZHANG Y, WU X P, ZHANG F, et al.

Analysis of evolutionary history and future development of Chinese rapeseed varieties based on genomic data

Chinese Journal of Oil Crop Sciences, 2020,42(3):325-333. (in Chinese with English abstract)

DOI:10.19802/j.issn.1007-9084.2019310      [本文引用: 1]

苏益,蔺万煌,马立英.

植物高亲和性钾离子转运系统

.安徽农业科学,2010,38(13):6646-6648.

[本文引用: 1]

SU Y, LIN W H, MA L Y.

High-affinity potassium transport system in plants

Journal of Anhui Agricultural Sciences, 2010,38(13):6646-6648. (in Chinese with English abstract)

[本文引用: 1]

/

〈 〉