油莎豆块茎表达CeTIP4;1基因的克隆、亚细胞定位及表达分析

关键词：油料作物；虎坚果；块茎；液泡膜内在蛋白；亚细胞定位；表达模式

中图分类号：S565.9 文献标志码：A

液泡膜内在蛋白（tonoplast intrinsic protein，TIP）隶属于水通道蛋白（aquaporin， AQP）家族，因其定位在液泡膜而得名[1]。TIP 最先在模式植物拟南芥中被鉴定，拟南芥γ-TIP1（AtTIP1;1）及其同源基因在蟾蜍卵母细胞和酵母中异源表达通常表现出高效的水分转运活性，高于或与细胞膜定位的质膜内在蛋白（plasma membrane intrinsicprotein， PIP）相当[2-3]。在体外，有些TIP 还被证实可以转运甘油、尿素、硼酸、NH3、H2O2 等中性小分子[4-10]。拟南芥（Arabidopsis thaliana）基因组含有10 个TIP 基因，并根据进化关系分为TIP1～TIP5 等5 个亚类，其中，除TIP4 和TIP5仅有1 个成员外，其他亚类含有2～3 个成员不等[11]。晶体结构显示，与PIP 类似，TIP 含有6 个由5个亲水环（LA–E）连接的跨膜α 螺旋（TM1–6）、2 个由LB 和LE 折向膜内形成的半螺旋以及位于半螺旋顶端的NPA 基序[10， 12]，2 个NPA 基序与位于TM2、TM5 和LE 上的4 个保守残基（即ar/R选择性滤器：H2、H5、LE1 和LE2）共同决定了孔道的大小和底物的特异性[13-14]；此外，拟南芥AtTIP2;1 位于LC 的H131 为NH3 的高效转运所必须[10]。与PIP 相比，拟南芥TIP 的ar/R 选择性滤器变异较大，TIP1 为H-I-A-V，TIP2、TIP3 和TIP4为H-I-G/A-R，TIP5 为N-V-G-C，孔径介于0.67～1.56 Å [15]。研究表明，不同TIP 亚类定位于不同类型的液泡，如TIP1 主要定位于裂解型液泡，TIP2 位于蛋白储藏型液泡，TIP3 定位于种子的蛋白储藏型液泡，并在种子的萌发过程中逐步被TIP1 所替代[16]。TIP4 也在裂解型液泡被发现，主要参与水分、尿素和甘油的转运[5， 9]。

油莎豆（Cyperus esculentus L.），又名虎坚果，是一种隶属于禾本目莎草科莎草属的新型草本油料作物[17-20]。不同于传统油料作物在种子中积累油脂，油莎豆不开花或开花不结实，但它可在地下块茎中积累高达35%的油脂，这使其成为研究营养组织油脂合成与调控的理想模型[21-23]。油莎豆块茎起源于地下匍匐茎，其发育过程主要分为起始、膨大和成熟等3 个阶段[21， 24]。通过前期研究，团队发现海南地区块茎的发育非常迅速，自起始后35 d 即可成熟，成熟块茎的含水量约为45%，而在此之前一般维持在85%左右[24]，表明水分平衡在块茎发育和代谢中的重要作用。前期研究显示，在新鲜成熟块茎的蛋白组中鉴定到多个AQP 基因，如CePIP1;1、CePIP2;1、CeTIP1;1和CeTIP2;1[19， 25-28]。其中，CeTIP1;1 隶属于TIP1亚类，在块茎、叶片、叶鞘、根、匍匐茎和芽尖等主要油莎豆组织中均高水平表达，且其编码蛋白定位在烟草叶片的液泡[27]。本研究报道另一个块茎表达的TIP 基因CeTIP4;1，该基因源于团队先前构建的油莎豆全长转录组[29]，其与AtTIP4;1同源。研究显示，CeTIP4;1 的表达模式及其编码蛋白的亚细胞定位与CeTIP1;1 明显不同，这些研究结果将有助于深入解析油莎豆块茎的水分平衡机制。

1 材料与方法

1.1 材料

本研究所用油莎豆品系为热研3 号，植株种植于中国热带农业科学院文昌实验基地，不同发育时期块茎的采集详见文献[24] ； 本氏烟草（Nicotiana benthamiana）种植于中国热带农业科学院热带生物技术研究所温室，具体生长条件详见文献[30] ； 大肠杆菌DH5α 和根癌农杆菌GV3101（内含pSoup-P19）感受态、亚细胞定位载体pNC-Cam1304-SubN 和pNC-HbPIP2;3-RFP[31]均由本实验室保存；各类酶、试剂盒及生化试剂详见文献[17]。

油莎豆的基因组和转录组数据分别下载于CNGBdb（https：//db.cngb.org/search/assembly/CNA-0051961/）和NCBI（https：//www.ncbi.nlm.nih.gov/bioproject/703731）数据库；拟南芥和水稻（Oryzasativa）的TIP 蛋白下载于TAIR11（https：//www.arabidopsis.org/）和RGAP7（http：//rice.plantbiology.msu.edu/）数据库；菠菜（Spinacia oleracea）的SoPIP2;1 下载于NCBI（https：//www.ncbi.nlm.nih.gov/Structure/pdb/1Z98）数据库。

1.2 方法

1.2.1 总RNA 的提取和反转录 油莎豆块茎总RNA 的提取及反转录方法参照文献[17]，即用天根植物多糖多酚RNA 提取试剂盒抽提总RNA，并用宝生物反转录试剂盒合成cDNA 第一链。

1.2.2 基因克隆与载体构建 本课题组前期通过对油莎豆的全长转录组[29]进行功能注释，发现1条与AtTIP4;1 同源的转录本，序列全长1098 bp，包含1 个753 bp 的开放读码框。为克隆该基因，采用Primer Premier 5.0 软件设计巢式PCR 引物对CeTIP4;1F/R （ AAACAAAGCAACCAAGACAG/CTTGCATCTCATGATCCAAA）和CeTIP4;1HF/R（ AGTGGTCTCTGTCCAGTCCTATGGCAAAGATTGCATTAGG/GGTCTCAGCAGACCACAAGTTCAAAATTCTTCACCATCTC）。PCR 扩增和载体构建参照文献[17]：首先利用外侧引物CeTIP4;1F/R 分离基因的全长转录本，随后采用内侧引物CeTIP4;1HF/R 扩增基因的编码区序列，最后运用NC 克隆试剂盒将目的片段克隆到亚细胞定位载体pNC-Cam1304-SubN。

1.2.3 序列的生物信息学分析 利用本地BLASTN 程序将上述CeTIP4;1 的转录本比对到油莎豆的基因组， 提取相应的基因序列后用GSDS 2.0（http：//gsds.gao-lab.org/）在线软件分析其基因结构；蛋白的理化特性和保守结构域分析分别使用Protparam（https：//web.expasy.org/protparam/）和CDD（https：//www.ncbi.nlm.nih.gov/cdd）在线软件；亚细胞定位、二级和3D 结构的预测分别使用Plant-mPLoc（http：//www.csbio.sjtu.edu.cn/bioinf/plant-multi/）、SOPMA（https：//npsaprabi.ibcp.fr/cgi-bin/npsa_automat.pl？page=/NPSA/npsa_sopma.html）和SWISS-MODEL（https：//swissmodel.expasy.org/interactive）在线软件；跨膜螺旋区及保守残基的鉴定基于与已结晶AtTIP2;1[10]和SoPIP2;1[12]的序列比对。

1.2.4 多序列比对与进化树的构建 蛋白的多序列比对采用MUSCLE 软件，进化树的构建采用MEGA（v6.06）软件，所用参数为最大似然法、1000 次自举重复。

1.2.5 亚细胞定位分析 以前期报道的pNCHbPIP2;3-RFP[31]作为细胞膜定位的阳性标记，参照文献[30]将重组质粒pNC-Cam1304-CeTIP4;1导入农杆菌，挑选单克隆，用含Kan 和Rif 的LB培养基摇菌，随后与含阳性标记的农杆菌工程菌等量混合，制备浸染液，转化受体，选取4 周龄的烟草叶片，荧光观察采用共聚焦显微镜ZeissLMS880。

1.2.6 基因表达分析 转录谱分析采用StringTie（v2.2.0）软件，即将转录组数据比对到基因组，然后用FKPM（fragments per kilobase of exon permillion fragments mapped）法均一化基因的相对表达水平。qRT-PCR 分析参照文献[20]，以CeUCE2（Fq： ATCATCAAGGAGACCCAGCG， Rq： CTTAGGGGCAGCCATAGGA）作为内参基因，所用基因特异性引物为CeTIP4;1Fq/Rq（GTCCTCACTTTCTCCCTCCTCT/CACTAACCCGACAAGCAACG）。

2 结果与分析

2.1 CeTIP4;1的克隆与序列分析

通过将CeTIP4;1 的转录本比对到基因组，发现该基因位于Scaffold40 反向链的5 229 867～5 231 478 位，基因全长1612 bp，包含2 个内含子，5ʹ和3ʹ UTR 分别为154、191 bp（图1A）。重组质粒pNC-Cam1304-CeTIP4;1 的桑格测序显示，克隆到的序列与基因的对应区域（外显子）完全一致。进一步的序列分析显示，CeTIP4;1 编码区的GC 含量为58.03%，编码250 个氨基酸（aa），理论分子量（MW）为25.57 kDa，等电点（pI）为6.02，总平均疏水指数（GRAVY）为0.884，脂肪族指数（AI）为119.40，不稳定系数（II）为33.70，属于稳定的酸性疏水型蛋白，这与CeTIP1;1、CeTIP2;1 和AtTIP2;1 类似，但不同于SoPIP2;1 的碱性特性（表1）。CDD 分析显示，CeTIP4;1 的第3～235 位残基为保守的MIP 结构域（图1B）。CeTIP4;1 与CeTIP1;1、CeTIP2;1、AtTIP2;1 和SoPIP2;1 的序列相似性分别为67.19%、70.92%、73.81%和44.08%，拥有6 个相对保守的跨膜螺旋（TM1–6）、2 个半螺旋（HB 和HE）以及2 个典型的NPA 基序（图1C）。CeTIP4;1 的ar/R 选择性滤器为H-I-A-R，不同于CeTIP1;1 的H-I-A-V、CeTIP2;1 和AtTIP2;1 的H-I-G-R 以及SoPIP2;1的F-H-T-R；CeTIP4;1 的Froger 位点为T-S-A-YW，与CeTIP1;1 和AtTIP2;1 一致，但不同于CeTIP2;1 的T-S-A-Y-I 和SoPIP2;1 的M-S-A-F-W。与SoPIP2;1 相比，CeTIP4;1、CeTIP1;1 和AtTIP2;1具有较短的N 端，且仅含有对应于SoPIP2;1 的S96 磷酸化位点；CeTIP4;1 的LC 区含有一个对应于AtTIP2;1 的H131 组氨酸，而CeTIP1;1 和SoPIP2;1 的相应位点分别为F 和N（图1C）。与CeTIP1;1、CeTIP2;1AtTIP2;1 和SoPIP2;1 类似，CeTIP4;1 的二级结构以无规则卷曲和α-螺旋为主，分别占36.80%和36.00%（表1）。基于AtTIP2;1的同源建模显示，CeTIP4;1 含有6 个跨膜螺旋，且可形成同源四聚体（图1D）。

2.2 CeTIP4;1 进化分析

为确认CeTIP4;1 的进化关系，本研究将其与已报道的AtTIP[11]和OsTIP[32]构建如图1E 所示的无根进化树。结果显示，这些蛋白被聚为5 组，分别对应于前期报道的TIP1～TIP5[11]，其中，CeTIP4;1与拟南芥和水稻中的TIP4 亚类聚在一起，其与AtTIP4;1、OsTIP4;1、OsTIP4;2 和OsTIP4;3 的序列相似性分别为78.57%、71.71%、69.62%和82.94%。以上结果证实CeTIP4;1 属于AtTIP4;1的直系同源基因，而水稻中的3 个成员可能在其与油莎豆分化之后产生。

2.3 CeTIP4;1 亚细胞定位分析

为揭示CeTIP4;1 所起的功能部位，利用含pNC-Cam1304-CeTIP4;1 的农杆菌工程菌对烟草叶片进行瞬时转化，共聚焦显微镜观察结果如图2 所示。虽然Plant-mPLoc 预测CeTIP4;1 定位在液泡膜，但本研究结果显示其定位在细胞膜，并与细胞膜标记HbPIP2;3-RFP 的红色荧光信号高度重合。

2.4 CeTIP4;1 基因表达分析

如图3A 所示，转录谱分析显示CeTIP4;1 主要在块茎、根、匍匐茎和芽尖中表达，而在幼嫩叶片、叶鞘和成熟叶片中的表达量极低，表现出明显的组织特异性。为进一步揭示基因在块茎发育过程中的表达模式，本研究选取起始期、膨大中期、膨大晚期和成熟期等4 个典型的时期进行qRT-PCR 分析，结果显示，CeTIP4;1 在4 个典型时期的表达量呈先升后降的趋势，峰值出现在膨大中期，最低的是成熟期（图3B）。

3 讨论

在植物中，细胞的水分平衡主要由细胞膜定位的PIP 和液泡膜定位的TIP 介导[1， 33]。PIP 高度保守，仅存在2 个进化小组，相比而言，TIP高度分化，其在高等植物中至少存在5 个进化小组[3， 11， 32， 34-37]。本研究报道的CeTIP4;1 隶属于TIP4 亚组，其与AtTIP4;1、OsTIP4;1、OsTIP4;2和OsTIP4;3 在蛋白水平的相似性介于69.62%～82.94%，且在进化分析中被聚在一起，推测它们可能具有类似的生物学功能。CeTIP4;1 含有2 个内含子，这与AtTIP4;1 和OsTIP4;1 一样，但不同于仅含有1 个内含子的OsTIP4;2 和OsTIP4;3[11， 32]。蓖麻（ Ricinus communis ）、麻疯树（ Jatrophacurcas）、木薯（Manihot esculenta）、橡胶树（Heveabrasiliensis）和杨树（Populus trichocarpa）等植物中的TIP4 成员均含有2 个内含子，且内含子位置高度保守[11， 34-37]，这极可能是OsTIP4;2 和OsTIP4;3 的第2 内含子发生了丢失， 而非CeTIP4;1 和OsTIP4;1 获得了新的内含子。

CeTIP4;1 含有1 个保守的MIP 结构域，其中包含6 个跨膜螺旋、2 个半螺旋以及2 个典型的NPA 基序，且具有与AtTIP4;1 和OsTIP4;3 完全一样的ar/R 选择性滤器（即H-I-A-R），其孔径大小与AtTIP1;1 相当（0.72 Å vs 0.67 Å）[15]，推测CeTIP4;1 可能具有AtTIP1;1 类似的水分转运活性。事实上，AtTIP4;1 和OsTIP4;1 均已被证实具有高效的水分转运活性[5， 9]。相比于SoPIP2;1[12]，CeTIP4;1 与AtTIP2;1[10]的结构比较接近，且其LC 区含有1 个与NH3 转运所必须的H 残基[10]，推测CeTIP4;1 可能具有高效的NH3 转运活性。此外， 研究还显示AtTIP4;1 可以转运尿素[5]，OsTIP4;1 可以转运甘油[9]。

在拟南芥原生质中的亚细胞定位分析显示AtTIP4;1 主要定位在液泡[5]，而本研究结果显示CeTIP4;1 定位在烟草叶片的细胞膜，其中的原由还有待进一步研究。可能的原因是：（1）TIP4 的亚细胞定位与组织类型有关，其中包括液泡的类型与数量，前人用的是原生质[5]，而本研究用的是幼嫩的烟草叶片；（2）TIP4 原本定位在液泡，但可以与PIP 互作从而定位到细胞膜，因为有研究显示AtTIP1;2、AtTIP2;1 和AtTIP3;1 均可与AtPIP2;1 互作[38]；（3）CeTIP4;1 本身就定位在细胞膜，但这与其pI<7 的酸性特性不太相符。

作为主要块茎表达的TIP1 基因，CeTIP4;1同时被发现在根、匍匐茎、芽尖中也有较高水平的表达，却很少在光合组织叶片和叶鞘中表达，这与CeTIP1;1 和CeTIP2;1 的组成型高水平表达不同[27-28]。在块茎发育过程中， CeTIP4;1 与CeTIP1;1 和CeTIP2;1 均表现出钟型趋势，但峰值出现不同，CeTIP4;1 和CeTIP2;1 较早，出现在膨大中期，而CeTIP1;1 出现在膨大晚期；此外，CeTIP4;1 和CeTIP2;1 在成熟期的表达量显著低于起始期，而CeTIP1;1 在起始期和成熟期的表达量差异不显著[27-28]。这表明CeTIP4;1 与CeTIP2;1类似，主要在块茎发育前期起作用，同时也解释了CeTIP4;1 未能在成熟块茎的蛋白组中被鉴定到的原因[19]。与CeTIP1;1 相比，CeTIP4;1 在块茎发育过程中的表达模式与含水量趋势更为吻合，因为在本研究分析的4 个时期中，起始期和膨大中期的含水量约为85%，膨大晚期为75%，成熟期为45%[24]。

综上，本研究首次报道了1个油莎豆TIP4 基因CeTIP4;1，该基因含有2 个保守的内含子，为AtTIP4;1 的直系同源基因；CeTIP4;1 定位在细胞膜，具有AtTIP2;1 类似的理化特性和结构特征；CeTIP4;1 的表达具有明显的组织特异性，且其在块茎发育过程中的表达模式与含水量趋势基本吻合。本研究结果为深入解析油莎豆块茎的水分平衡机制及遗传改良奠定坚实的基础。