Mg2+是生物体内广泛存在的二价金属阳离子,其在植物的各种生理生化过程中起重要作用[1-2]。Mg2+具有较大的水合半径、小的离子半径和高的电荷密度[1],这些物理结构决定了其转运蛋白的特异性。在高等植物中,Mg2+的转运主要由镁离子转运蛋白（magnesium transporter/mitochondrial RNA splicing 2, MGT/MRS2）和Mg2+/ H+交换体（magnesium-proton exchanger, MHX）介导[2⇓⇓-5]。其中,MGT基因最早在鼠伤寒沙门氏菌（Salmonella typhimurium）中被发现,缺失Mg2+转运系统（即CorA、MgtA和MgtB）的突变株MM281无法在低于10 mmol/L Mg2+的N-Minimal培养基上生长[6-7]。至今MGT家族成员已在拟南芥（Arabidopsis thaliana）、水稻（Oryza sativa）等植物中进行了较为深入的研究[8⇓-10]。在拟南芥的10个MGT中,有9个基因的编码蛋白具有Mg2+转运活性,其中,AtMGT1、AtMGT2和AtMGT10为高亲和型,AtMGT3、AtMGT6、AtMGT7和AtMGT9为低亲和型,AtMGT5为双亲和型[11];对atmgt2、atmgt3、atmgt7等突变体的研究表明,其胞质中的Mg2+浓度明显低于野生型[7,12]。在水稻的9个家族成员中,OsMGT1的上调表达可明显增加细胞中的Mg2+浓度,从而增强植株的耐铝胁迫能力[10]。结构分析显示,MGT蛋白序列的C-末端通常含有两个保守的疏水跨膜结构域（TM1和TM2）及位于TM1末端的三肽基序GMN（Gly-Met-Asn）[13]。
番木瓜（Carica papaya L.）是一种隶属于十字花目番木瓜科的多年生常绿果树,含多种维生素、氨基酸、钾、钙、镁、铁等,被誉为“水果之王”[14]。本研究基于番木瓜的基因组和转录组数据,对一个果实高丰度表达基因CpMGT1进行克隆、生物信息学及表达模式分析,并在此基础上利用沙门氏菌突变株MM281进行功能互补,以期为今后的遗传改良奠定基础。
1 材料与方法
1.1 材料
1.1.1 植物材料
供试番木瓜品系为中白,种植于中国热带农业科学院文昌试验基地,树龄为7个月。选取长势基本一致的植株分成3组,每组3个单株,花授粉后即进行标记,最后于同一天早上8：00之前采集授粉后10、60、120 d的果实及根、茎、叶、花等其他组织,授粉后120 d的果实用报纸包裹后置于25℃培养室中,分别放置6 h、2 d和6 d,每样品3次生物学重复,样品采集后立即放入液氮中速冻,并置于-80℃保存备用。
1.1.2 菌株及载体
大肠杆菌DH5α由本实验室保存;沙门氏菌突变株MM281和质粒pTrc99A由中国热带农业科学院橡胶研究所阳江华副研究员惠赠。
1.1.3 主要试剂
反转录试剂盒PrimeScript II 1st Strand cDNA Synthesis Kit和高保真酶PrimerSTAR Max DNA Polymerase购自宝日医生物技术（北京）有限公司;胶回收试剂盒Gel Extraction Kit和质粒提取试剂盒Plasmid DNA Mini KitⅠ购自Omega Bio-Tek;各类限制性内切酶均购自美国NEB;重组克隆试剂盒ClonExpress II One Step Cloning Kit购自诺唯赞;其余试剂均为国产分析纯。
1.2 方法
1.2.1 总RNA的提取及cDNA的合成
总RNA的提取参照房平昌[15]的方法,主要流程如下：取样品2~3 g,在液氮下用研钵研磨成细粉末;加入10 mL 65℃预热的CTAB提取液和200 μL β-巯基乙醇,涡旋剧烈混匀,65℃水浴5~20 min,中间振荡数次;待冷却到室温后加入等体积的氯仿,振荡混匀,4℃ 15 000×g离心10 min;将上清液转至另一新管中,并加入等体积的水饱和酚和氯仿,充分振荡,4℃ 15 000×g离心10 min（视情况重复上述步骤1~2次）;将上清吸出后,加入总体积1/3的8 mol/L LiCl,-20℃冰箱放置12 h,再次离心,取上清,75%乙醇洗涤沉淀2次,沉淀用50~100 μL DEPC水溶解。质检合格后参照TAKARA反转录试剂盒说明书合成cDNA,-20℃保存待用。
1.2.2 基因克隆与载体构建
通过搜索番木瓜的基因组和转录组数据（https://www.ncbi.nlm.nih.gov/）获得CpMGT1的全长转录本序列,并据此设计如表1所示的3对基因特异性引物。以不同组织混合来源的cDNA作为模板,采用巢式PCR法进行扩增,首轮反应体系为PrimeSTAR Max premix 2×12.5 μL、cDNA 1 μL、CpMGT1F/ CpMGT1R各1 μL、ddH2O 10.5 μL;扩增程序为95℃预变性1 min,98℃变性10 s、53℃退火25 s、72℃延伸1.5 min（35个循环）,最后72℃延伸5 min。PCR产物经电泳检测后切胶回收,稀释50倍作为第二轮PCR扩增的模板,反应体系为PrimeSTAR Max premix 2×50 μL、第一轮PCR产物为cDNA 1 μL;CpMGT1-99F/CpMGT1- 99R各1 μL、ddH2O 22 μL;扩增程序为95℃预变性1 min,98℃变性11 s、54℃退火25 s、72℃延伸1.5 min（25个循环）,最后72℃延伸5 min。表达载体pTrc99A使用限制性核酸内切酶EcoRⅠ（10 U/μL）和PstⅠ（10 U/μL）进行酶切,切胶回收后与第二轮PCR获得的目的基因片段按1∶3比例混匀,并用同源重组酶Exnase II 37℃水浴连接2 h,然后采用冻融法[16]转入DH5α大肠杆菌感受态细胞,菌落PCR验证后选取阳性克隆进行测序分析。
表1 本研究所用引物及其序列Tab. 1 Primer sequence used in this study |
| 引物名称Primer name | 引物序列（5°-3°）Primer sequence (5°-3°) | 目的Usage |
|---|---|---|
| CpM1-115SF | CTCACACAACTCTTCCATCTCT | 全长编码区克隆 |
| CpM1-273SR | CGAGGTCCTGAGAATTACATAC | |
| CpM1-1685AF | CTAATTCCTTAATAGCTTGACTGC | |
| CpM1-1758AR | TTTTGAGCTAGACCAACTGATC | |
| CpMGT1-1F | GTAGTTGCTATTTTTGGGGTGG | 荧光定量分析 |
| CpMGT1-1R | ACCAAACAAACGCACAGAATATG | |
| CpMGT1-2F | AGCTGATCCAGTTTGAGCTG | |
| CpMGT1-2R | CTCCAGTTATTATGAGGACCCATT | |
| CpEIEFF | AGGCAGGCAAGAGAAGAT | |
| CpEIEFR | TTCATACCGAGTAGCGATTC | |
| CpM1-99-hoF | ATTTCACACAGGAAACAGACCATGGCGGACCTCAAAGAGC | |
| CpM1-99-hoR | CTCTCATCCGCCAAAACAGCCTCACAGCGGCATAAGTCTTCTA |
1.2.3 生物信息学分析
用在线软件Expasy（https://web.expasy.org）分析理论分子量和等电点;用NCBI CDD（https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi）分析保守结构域;用在线软件WoLF PSORT（https://www.genscript.com/wolf-psort.html）分析亚细胞定位;用MEGA 6.0软件中的MUSCLE进行序列比对,并用邻接法构建进化树。
1.2.4 基因的表达分析
1.2.5 沙门氏菌功能互补分析
利用化学转化法将重组质粒pTrc99A-CpMGT1转入MM281感受态中,选取阳性单克隆摇菌,阳性对照为pTrc99A- AtMGT1,阴性对照为MM281突变株及转pTrc99A空载的MM281突变株。OD600摇至0.6后,加入1/103 IPTG,将菌液继续摇至1.0,并用N-Minimal培养基洗涤2~3次,分别稀释到5个浓度梯度。准备4个不同Mg2+浓度的N-Minimal固体培养皿,将处理好的菌液呈梯度式点至培养皿,倒置培养2~3 d后,观察长势并拍照。
2 结果与分析
2.1 基因克隆
如图1所示,PCR扩增成功获得1条约1300 bp的目的条带,而pTrc99A的双酶切获得1条约4100 bp的目的条带,将其孵育后构建原核表达载体pTrc99A-CpMGT1。
2.2 生物信息学分析
序列分析表明,CpMGT1的编码区（CDS）为1332 bp（NCBI登录号为XM_022050807）,预测编码443 aa,理论分子量（MW）为50.43 kDa,等电点（pI）为5.12。CDD分析和序列比对显示,CpMGT1含有保守的Mrs2_Mfm1p-like（cd12823）结构域;在C端含有2个疏水跨膜区（即TM1和TM2）,且TM1包含高度保守的GMN基序（图2）。WoLF PSORT分析显示,CpMGT1定位在细胞膜。
图2 CpMGT1与拟南芥和水稻中同源蛋白的多序列比对保守的CorA结构域、C端的疏水跨膜区TM1和TM2用下划线标识;GMN基序用线框标识。 Fig. 2 Multiple sequence alignment of CpMGT1 with homologs in Arabidopsis and rice The conserved CorA domain and two hydrophobic transmembrane regions (i.e. TM1 and TM2) were underlined, whereas the GMN motif at the TM1 were boxed. |
2.3 进化分析
为揭示CpMGT1的进化特征,将其与已报道的拟南芥和水稻MGT蛋白构建进化树（图3）。结果表明,这些MGT被聚为5支,其中,CpMGT1与拟南芥的AtMGT1（NP_001322724）、AtMGT2（NP_001077545）、AtMGT3（NP_001118259）,水稻的OsMRS2-1（XP_015643498）、OsMRS2-9（XP_015633967）聚在一起。相比而言,CpMGT1与AtMGT1和AtMGT2的相似性更高,分别为83.33%和86.04%。
2.4 基因表达分析
2.5 功能互补分析
为鉴定CpMGT1的Mg2+转运活性,将构建的重组质粒pTrc99A-CpMGT1转入MM281突变株中。功能互补实验显示,实验组和对照组在10 mmol/L Mg2+的平板上均可正常生长。然而,当培养基中的Mg2+浓度降低到5 mmol/L时,阴性对照未见生长,而实验组与阳性对照长势相当。随着Mg2+浓度的进一步降低,实验组和阳性对照的生长虽然都呈现减弱趋势,但是趋势相当,这表明CpMGT1是与AtMGT1活性相当的高亲和型Mg2+转运蛋白（图5）。
3 讨论
本研究首次报道了一个隶属于MGT基因家族的番木瓜Mg2+转运蛋白基因CpMGT1,该基因编码443 aa、分子量为50.43 kDa、等电点为5.12,细胞膜定位。更重要的是,蛋白含有该家族特有Mrs2_Mfm1p-like或CorA（PFAM登陆号PF01544）结构域及位于C端疏水跨膜区TM1上的GMN基序。早前的研究表明,GMN基序为Mg2+转运所必需[4,13]。根据基因表达分析结果,CpMGT1总体为一组成型表达基因,其中在根、茎和果实中的表达量较高,暗示该基因可能在多个部位发挥Mg2+转运功能。根据进化分析结果,CpMGT1与拟南芥中AtMGT1和AtMGT2的亲缘关系较近,相似性高达80%以上,这不仅暗示它们之间的直系同源关系及类似的生物学功能,同时,该结果跟番木瓜与拟南芥同属十字花目的分类学关系是一致的[19]。前人的研究结果表明,AtMGT1和AtMGT2均可互补沙门氏菌突变株MM281,且二者均编码高亲和型Mg2+转运蛋白[7,12]。AtMGT1定位于细胞质膜,具有较高的组织特异性,主要在根部和表皮中表达;基因表达受Al3+诱导,其过量表达可增强转基因植株对铝毒的抗性[27⇓-29]。AtMGT2定位于液泡膜,主要在幼嫩叶片中发挥Mg2+转运功能[12],并与AtMGT3共同维持叶肉细胞中液泡的Mg2+平衡[20]。此外,有研究显示,AtMGT2还能互补Mg2+转运缺陷型酵母突变株CM66,并显著提高酵母的生长效率[19]。与进化结果相应的是,CpMGT1的异源过表达可互补MM281突变株,并显著提高工程菌在低Mg2+浓度下的生长速率,这初步证实CpMGT1编码一个高亲和的Mg2+转运蛋白,且其活性与AtMGT1相当。虽然如此,CpMGT1具体是如何调控Mg2+平衡的分子机制还有待深入研究。
综上,本研究完成了番木瓜CpMGT1基因的克隆及序列特征、进化关系、表达模式和Mg2+转运活性分析。这些结果不仅为进一步揭示番木瓜MGT基因在不同组织特别是在果实中Mg2+的积累机制奠定了坚实的基础,也为今后的品种改良提供了宝贵的资源。