0 引言
【研究意义】转运蛋白（transport protein）是膜蛋白的一大类,介导生物膜内外的化学物质以及信号交换。多药物和有毒化合物外排家族（multidrug and toxic compound extrusion family,MATE）,属于次级转运蛋白家族,是生物中5类解毒输出转运蛋白家族之一[1]。近年来,越来越多的遗传学、生物化学和分子生物学研究表明MATE转运蛋白参与植物类黄酮运输。原花青素（proanthocyanidin,PA）是植物酚类物质中普遍存在的组分之一,在植物界广泛存在,呈现出多样的生物和生物化学活性[2,3,4,5,6],并对人体健康有潜在的有益作用[7,8,9,10,11]。荞麦属于蓼科荞麦属植物,是一种具有很高食用价值和医疗保健作用的经济作物,也是一种极具开发潜力的功能性食品原料。荞麦mate的研究为解析MATE转运蛋白家族在荞麦次生代谢产物运输和逆境胁迫中的作用奠定基础。【前人研究进展】MATE转运蛋白首先在副溶血性弧菌和大肠杆菌中被发现,功能为多药物排出,但由于与其他转运蛋白缺乏序列同源性而被认定为新的转运蛋白家族并命名为MATE家族[12,13]。随着研究深入,研究****发现MATE家族成员分布广泛,遍布古细菌、细菌、酵母、动物和植物,且家族直系同源成员较多,如研究发现拟南芥中至少有56个MATE家族基因[14],亚洲棉中有34个,雷蒙德氏棉中有42个[15]。现阶段对于MATE家族大多数成员的功能知之甚少,已研究得到的植物中MATE型转运蛋白的功能涉及异生物外排,次生代谢物包括生物碱和类黄酮的积累、铁（Fe3+）易位、铝（Al3+）解毒和植物激素信号传导等[16]。拟南芥tt12编码一种类黄酮循环必需的、属于MATE家族的多药次生类转运蛋白,其生物学功能是将胞浆内合成的花青素、原花青素等多酚类色素的单体物转运到液泡等亚细胞器官中,然后花青素在酸性条件下显色或原花青素进一步聚合后形成有色色素[17];参与Al3+解毒的MATE转运蛋白首先在高粱（Sorghum bicolor,SbMATE）和大麦（Hordeum vulgare,HvAACT1）中被鉴定;它们的基因表达水平与从根部释放的柠檬酸的量以及在Al3+存在下的根伸长率相关,表明这些MATE转运蛋白的高表达对于Al3+耐受是必需的[18,19,20]。拟南芥Atfrd3编码MATE型转运蛋白负责将柠檬酸盐（Fe的螯合剂）转运到根木质部以高效地转运铁[21]。【本研究切入点】荞麦种子转录组测序中[22],发现mate注释的unigene共有54个。但MATE在荞麦中的研究报道较少,只有LEI等[23]在荞麦中发现2个mate——FeMATE1和FeMATE2,并通过功能分析得到两者与荞麦中Al离子胁迫相关。【拟解决的关键问题】本研究拟通过RACE的方法从荞麦中克隆2条mate全长,并对其功能进行验证,以期为筛选高营养成分含量和抗逆性较强的种质资源提供依据。1 材料与方法
1.1 植物材料
试验所用材料为普通荞麦（Fagopyrum esculentum Moench）,来自国家中期种质库（No.000677）,于2016年9月栽培于中国农业科学院作物科学研究所温室中。采集成株期的根、茎、叶、花苞和花及开花后7、14、21和28 d的籽粒经液氮速冻后于-80℃保存,用于总RNA的提取。1.2 总RNA的提取和cDNA的合成
取液氮中保存的材料约100 mg,采用RNeasy Plant Mini Kit（QIAGEN,USA）提取植物样品总RNA。参照QuantiTect Reverse Transcription Kit（QIAGEN,USA）说明书,以纯化好的2 μg RNA为模板逆转录合成cDNA,保存于-20℃备用。1.3 荞麦mate全长cDNA克隆
用BLAST在线网上工具对荞麦转录组测序已获得的部分mate序列与其他物种已发表的基因序列进行比对,选取保守区域设计引物扩增合成cDNA中间片段,引物序列信息列于表1,扩增产物经电泳回收纯化后连接pMD19-T载体并转化大肠杆菌DH5α感受态细胞,经蓝白斑筛选重组子,挑取阳性克隆送至生工（上海生物工程股份有限公司）测序,并在NCBI数据库中比对分析测序结果。Table 1
表1
表1荞麦mate克隆所用引物序列信息
Table 1Sequence information of primers used for cloning buckwheat mate gene
| 引物名称 Primers name | 引物序列 Primers sequence (5'-3') |
|---|---|
| FeTT12-F | TGGTGGCCGAGGCTAGTGGCTTGG |
| FeTT12-R | TATCATTCCCCACCATATACC |
| FeMATE3-F | GCTGATCCAATTGCTTCTCTAGT |
| FeMATE3-R | ACCAGTCCCAATCCTCCAGAATCCTG |
| FeTT12-cDNA-F | GGACTCTCTTCAACCCACCCCATT |
| FeTT12-cDNA-R | CTCTCTTCAACCCACCCCATTTCT |
| FeMATE3-cDNA-F | ATGGCGGAAAAAGTTGGGTTTATTG |
| FeMATE3-cDNA-R | GACCTTGGAAATTCCTGAGTTAA |
| FeTT12-3-out | gattacgccaagcttTGCTTGCCTTGCTCGTCCACATCT |
| FeTT12-3-inner | gattacgccaagcttTCATGGTGGGTTATCGTCTTGGGG |
| FeTT12-5-out | gattacgccaagcttGCATGACCCCAGAAGAAGCAGAGA |
| FeTT12-5-iner | gattacgccaagcttCAGCCCCAAGACGATAACCCACCA |
| FeMATE3-3-out | gattacgccaagcttTACACCCAGCACAGCAGTACTTGA |
| FeMATE3-3-inner | gattacgccaagcttAGCTATGCAAGGAGTTTTTCGCGG |
| FeMATE3-5-out | gattacgccaagcttACTGTGATATGACATGGGCAATGG |
| FeMATE3-5-inner | gattacgccaagcttTGAAACGTGGAGGCAATAGATCAA |
| FeTT12-RT-F | GAGTGAGGGTGGCGAATG |
| FeTT12-RT-R | TGAGGAGGATGGTGAAGG |
| Histone3-F | GAAAGGCTCCTAGGAAGCAA |
| Histone3-R | GAACAAGCCTCTGGAATGGCA |
新窗口打开
以获得的cDNA中间片段为模板,设计荞麦mate的RACE用特异引物（表1）。根据SMARTer RACE 5′/3′ Kit（TaKaRa, Dalian）说明书进行巢式PCR扩增。利用DNAMAN软件将mate的5′端序列和3′端序列以及中间cDNA序列进行拼接获得荞麦mate全长cDNA序列。参照拼接所得cDNA序列设计引物（表1）,以cDNA为模板扩增出mate全长cDNA序列。
1.4 荞麦mate生物信息学分析
利用DNAMAN软件搜索并翻译mate的ORF区域,得到的氨基酸序列与已发表的其他物种基因编码的氨基酸序列进行同源分析。利用MEGA软件中的NJ法构建荞麦MATE蛋白与其他物种MATE的系统发育树。通过ExPASy中的protparam对MATE蛋白的分子量及等电点进行预测;利用ProtScale对该蛋白的氨基酸序列进行疏水性分析;利用TMHMMv对MATE蛋白的跨膜区域进行预测;利用SOPMA预测MATE蛋白的二级结构。1.5 qRT-PCR分析
依据Fett12 cDNA序列设计引物FeTT12-RT- F/FeTT12-RT-R（表1）,以荞麦根、茎、叶、花和花苞,开花后7、14、21和28 d的种子材料cDNA为模板进行PCR扩增,以检测各部分Fett12的表达水平,以H3为内参基因,每个样品重复3次。1.6 荞麦原花青素（PAs）含量的测定
参照吴楠等[24]DMACA-HCl的方法测定荞麦花苞、花及不同发育阶段种子中可溶的和不可溶的PA含量。取组织冻干样品100 mg,加入10 mL甲醇,超声波萃取1 h,其间每隔20 min旋涡震荡1 min,4 500 r/min离心10 min,上清即为可溶的PAs;用10 mL含有1% HCl的甲醇重悬沉淀,60℃水浴1 h;水浴后,4℃ 4 500 r/min离心10 min,上清即为不溶性的PA。采用分光光度法对原花青素含量进行测定,以儿茶素标品为对照,作标准曲线。分别取200 µL含有原花青素的提取液和2.8 mL的反应液（含有0.1% DMACA和5% HCl的甲醇溶液）加入5 mL离心管中,室温下反应15 min,然后测定643 nm的吸光度,原花青素的含量参照儿茶素的标准曲线进行计算。2 结果
2.1 荞麦mate克隆
利用RACE克隆的方法,从荞麦中获得2条mate全长cDNA序列,分别命名为Fett12和Femate3。核苷酸序列分析表明,Fett12的全长cDNA为1 707 bp,包含69 bp的5'非翻译区,159 bp的3'非翻译区和1 479 bp的开放阅读框;Femate3全长cDNA为2 163 bp,包含1 548 bp的开放阅读框、387 bp的5'非翻译区和228 bp的3'非翻译区。通过ExPASy中的protparam在线分析得到Fett12编码一个含有492个氨基酸残基的假定蛋白,其分子量为53.81 kD,等电点为6.75;Femate3编码一个含有516个氨基酸残基的假定蛋白,分子量为56.12 kD,等电点为6.52。进一步对其疏水性,跨膜区域及二级结构分析预测（图1）,Fett12 编码的假定蛋白二级结构主要以α螺旋为主,占44.51%,含有10个跨膜结构域;Femate3编码的假定蛋白二级结构同样以α螺旋为主,占54.84%,含有9个跨膜结构域。
显示原图|下载原图ZIP|生成PPT图1FeTT12和FeMATE3 跨膜区及二级结构预测
A：FeTT12跨膜区域预测图;B：FeTT12二级结构预测图;C：FeMATE3跨膜区域预测图;D：FeMATE3二级结构预测图
-->Fig. 1FeTT12 and FeMATE3 transmembrane regions and secondary structure predictions
A: FeTT12 transmembrane region prediction; B: FeTT12 secondary structure prediction; C: FeMATE3 transmembrane region prediction; D: FeMATE3 secondary structure prediction
-->
2.2 荞麦FeTT12和FeMATE3的系统进化与同源分析
目前,研究发现植物MATE型转运蛋白的转运底物有烟碱、原花青素、花青素、柠檬酸盐、脱落酸等,并以此发挥不同的功能[1]。本研究选取已发表的其他物种的原花青素MATE转运蛋白、花青素MATE转运蛋白、柠檬酸盐MATE转运蛋白的氨基酸序列与克隆的2条荞麦MATE氨基酸序列构建系统发育树（图2）。从图中可以看出,选取的氨基酸序列按其功能被分成了3组,FeTT12被分到了第一组,暗示其在荞麦中可能参与原花青素的转运与累积。进一步将FeTT12氨基酸序列与其他植物TT12编码的氨基酸序列进行比对（图3）和同源性分析（表2）,得到荞麦FeTT12与其他物种的TT12蛋白的氨基酸序列具有极高的同源性,其中与桑属TT12蛋白同源性最高,为77.3%;与拟南芥TT12蛋白同源性最低,为41.5%。
显示原图|下载原图ZIP|生成PPT图2FeTT12和FeMATE3系统进化树分析
MtTT12：蒺藜苜蓿Medicago truncatula,FJ858726.1;GmTT12：大豆Glycine max,XM_003545107.1;MdTT12：苹果Malus domestica,NP_001280841.1;GhTT12：棉花Gossypium hirsutum,AGW32085.1;AtTT12：拟南芥Arabidopsis thaliana,NM_115765.3;ZmTT12：玉米Zea mays,ACG28719;MnTT12：川桑Morus notabilis,EXC08256.1;VvMATE：葡萄Vitis vinifera,ACN88706.1;MtMATE2：蒺藜苜蓿Medicago truncatula,ADV04045.1;CcMATE：木豆Cajanus cajan,XP 020236867.1;HvMATE：大麦Hordeum vulgare,AMO65654.1;SbMATE：高粱Sorghum bicolor,ABS89149.1;VvMATE：芒果Vigna mungo,ALX81166.1;BnMATE：油菜Brassica napus,NP 001302612.1。下同The same as below
-->Fig. 2FeTT12 and FeMATE3 phylogenetic trees
-->
显示原图|下载原图ZIP|生成PPT图3荞麦Fett12与已公开其他植物tt12编码的氨基酸序列比对
下划线表示跨膜结构域
-->Fig. 3Amino acid sequence alignment of Fett12 gene of buckwheat with the tt12 gene of other plants
Lines denote regions represent transmembrane domains
-->
Table 2
表2
表2FeTT12与其他物种TT12之间的同源分析
Table 2The homology analysis of FeTT12 with TT12 from other species (%)
| 蛋白名称 Protein name | GmTT12 | GhTT12 | MdTT12 | MtTT12 | AtTT12 | MnTT12 | ZmTT12 | FeTT12 |
|---|---|---|---|---|---|---|---|---|
| GmTT12 | – | |||||||
| GhTT12 | 76.80 | – | ||||||
| MdTT12 | 76.80 | 82.60 | – | |||||
| MtTT12 | 83.30 | 80.00 | 76.10 | – | ||||
| AtTT12 | 72.20 | 76.20 | 72.50 | 70.70 | – | |||
| MnTT12 | 41.30 | 41.60 | 41.60 | 41.90 | 40.40 | – | ||
| ZmTT12 | 41.00 | 41.70 | 40.30 | 41.90 | 40.10 | 66.90 | – | |
| FeTT12 | 42.40 | 43.50 | 43.80 | 43.80 | 41.50 | 77.30 | 64.90 | – |
新窗口打开
FeMATE 3被分到了第三组,暗示其功能可能与荞麦的离子胁迫相关。通过与其他物种中转运柠檬酸盐的MATE型蛋白氨基酸多序列比对（图4）及同源分析（表3）,发现其与Al离子胁迫相关的荞麦MATE蛋白FeMATE2的一致性达到96.33%,由此可以说明两者为同一基因序列,但可能由于选用材料的不同,导致存在个别氨基酸残基差异（红色标注）,主要集中在第二个跨膜区与第三个跨膜区之间和羧基端。
显示原图|下载原图ZIP|生成PPT图4FeMATE3氨基酸序列比对图
FeMATE2：荞麦Fagopyrum esculentum;AtMATE1：拟南芥Arabidopsis thaliana,NP_974000;OsFRDL2：水稻Oryza sativa, XP_015615043;红色标注为FeMATE2与FeMATE3有差异的氨基酸残基; 下划线表示跨膜结构域
-->Fig. 4Multiple alignment of FeMATE3
Red is marked as different amino acid residue among FeMATE2 and FeMATE3; Lines denote regions represent transmembrane domains
-->
Table 3
表3
表3FeMATE3同源分析表
Table 3The homology analysis of FeMATE3 (%)
| 蛋白名称 Protein name | FeMATE3 | FeMATE2 | AtMATE1 | OsFRDL2 |
|---|---|---|---|---|
| FeMATE3 | – | |||
| FeMATE2 | 96.33 | – | ||
| AtMATE1 | 69.17 | 69.04 | – | |
| OsFRDL2 | 60.70 | 61.51 | 59.15 | – |
新窗口打开
2.3 荞麦Fett12在不同组织中的表达模式
为进一步分析FeTT12在荞麦各个组织中相对表达量,以荞麦的根、茎、叶、花苞、花和种子为材料进行qRT-PCR分析。结果（图5）表明,Ftt12在荞麦的各个组织中均有表达,其中在叶中的表达量最高。为研究Fett12在荞麦种子不同发育阶段种子中的表达模式,同时也测定了荞麦开花后7、14、21和28 d种子中Fett12的相对表达量,结果表明,伴随着荞麦种子的不断发育成熟,Fett12的相对表达量逐渐增高,在成熟期（28DPA）种子中的表达量最高。
显示原图|下载原图ZIP|生成PPT图5Fett12在荞麦不同组织中的相对表达量
-->Fig. 5Relative expression of Fett12 in different tissues of buckwheat
-->
2.4 荞麦不同组织中PA的含量测定
基于同源分析得到FeTT12可能与原花青素的转运和累积相关,进一步检测荞麦不同组织及不同发育阶段种子中可溶和不可溶的PA的含量。相比其他组织,PA在荞麦花中的含量明显较高;在不同发育阶段种子中的含量模式为随着种子的不断成熟,PA的含量逐渐降低;其中可溶的PA的含量在从14DPA到21DPA阶段的种子材料中下降明显,21DPA到28DPA阶段的种子材料中变化不明显。
显示原图|下载原图ZIP|生成PPT图6荞麦不同组织中可溶和不可溶的PA的含量测定
A：可溶的PA含量;B：不可溶的PA含量。不同小写字母表示差异显著（P<0.05）
-->Fig. 6Determination of PAs (soluble and insoluble) in different groups of buckwheat
A: Soluble PAs content; B: Insoluble PAs content. Different letters indicate significant differences (P<0.05)
-->
3 讨论
3.1 荞麦mate序列特征
研究发现,所有的MATE蛋白没有明显的共有保守序列,但大多数成员都具有大约40%的序列相似性[25]。在结构上,由400—550个氨基酸残基组成,具有12个跨膜螺旋。在荞麦中克隆到的2条MATE蛋白,分别由492和516个氨基酸残基组成;与其余已发表物种的MATE蛋白的序列相似性都达到40%及以上,编码的假定蛋白均为疏水性极强的膜蛋白。通过进一步构建系统发育树及同源分析,发现2条mate被分到了不同的组里,其中一条mate为拟南芥tt12在荞麦中的同源基因,并因此命名为Fett12;Femate3与其他物种中转运底物为柠檬酸盐的mate聚集到一起,并与LEI等[23]在荞麦中发现的Al离子胁迫相关的Femate2的序列相似性达到96%,只有个别氨基酸残基存在差异,主要集中在第二个跨膜区与第三个跨膜区之间,同时这一区段也是不同物种这一同源基因进行氨基酸序列比对时差异最多的区段。结合跨膜区预测结果分析获得,其差异主要存在于胞内区域和跨膜区域且跨膜区域的氨基酸差异之间的极性一致。3.2 FeTT12参与荞麦中原花青素的累积
荞麦是一种高黄酮含量的作物,有着与其他植物物种相似的类黄酮生物合成通路。原花青素作为合成通路中的一员,其在荞麦中的含量,合成与转运等被广泛的研究。PA前体的聚合及随后转化为棕色氧化产物发生在液泡中,PA和花青素也积累在其中[26]。但花青素还原酶（anthocyanidin reductase,ANR）是一种可溶性胞质酶,因此,用于PA生物合成的表儿茶素或其衍生物可能具有外源性,需要一个机制将其从细胞质转运至液泡中。拟南芥tt12编码一个定位于液泡膜的MATE型转运蛋白,tt12突变体的显微分析显示原花青素在拟南芥内皮细胞空泡中的积累强烈减少[17]。进一步研究表明拟南芥TT12在体内功能性运输表儿茶素3′-O-葡萄糖苷作为原花青素合成的前体[27]。本次研究中克隆得到的Fett12为拟南芥tt12在荞麦中的同源基因;同时,与棉花[24]、苹果[28]、桑[29]、葡萄[30]等物种中的tt12同源基因有较高的序列相似度,因此FeTT12有可能在荞麦中发挥转运原花青素前体的功能。随后,研究了Fett12在荞麦不同组织中的表达差异和PA的含量之间的相关性。基因表达模式研究表明其具有组织表达特异性,且在叶中的表达量最高;在PA的含量测定中,发现作为次级代谢产物PA的分布具有明显的组织器官差异,且在荞麦花中的含量明显较高,这也与ZIELIŃSKA等[31]通过测定荞麦茎、叶、花、未成熟与成熟种子的总黄酮及芦丁含量,得到其在荞麦花中的含量最高的结果一致。因之前很多研究表明,PA主要积累在植物种子中[25],因此我们也测定了Fett12在荞麦不同成熟阶段种子中的表达量及PA的含量变化。随着种子的不断成熟Fett12的表达量逐渐增高,PA的含量逐渐降低,表明两者之间存在一定关联。在可溶的PA的含量测定中,其在14DPA到21DPA阶段的种子材料中下降明显,可能与此阶段合成结构复杂的原花青素进而转变为不溶的PA相关。3.3 FeMATE3与荞麦的Al离子胁迫相关
与其他物种如小麦相比,荞麦对Al毒害具有较高的耐受性;此外,它还在叶子中累积高Al而不会表现出毒性症状[32,33,34]。植物MATE的功能研究中,发现MATE型转运蛋白能通过转运柠檬酸盐来解除重金属对植物的毒害。此次在荞麦中克隆到的Femate3在系统进化分析中被归类到柠檬酸盐MATE型转运蛋白,表明其可能在荞麦中Fe异位或Al解毒发挥重要作用。参考LEI等[23]最近发表的荞麦中与AL胁迫相关的2条mate,通过多序列比对,发现FeMATE3与他发表的FeMATE2的氨基酸一致性达到96%,进而可以一定程度说明FeMATE3在荞麦Al解毒方面发挥重要作用。4 结论
从荞麦中成功克隆获得2条mate——Fett12和Femate3,且Fett12与荞麦中原花青素的转运和累积相关。The authors have declared that no competing interests exist.
