【目的】从芝麻全基因组中分离鉴定水通道蛋白AQP(aquaporin)家族基因,并进行系统进化关系、连锁群定位、基因结构、跨膜结构域和亚细胞定位预测、保守性氨基酸残基以及青枯雷尔氏菌诱导表达分析,为芝麻AQP的功能研究与利用奠定基础。【方法】通过生物信息学手段,结合芝麻基因组注释信息,鉴定芝麻AQP家族成员序列信息,并用InterPro逐一进行验证。利用ClustalW2对芝麻、拟南芥和水稻的AQP以及马铃薯的XIPs进行多序列比对,用MEGA6.0构建进化树。通过MapInspect和Gene Structure Display Server 2.0进行连锁群定位和基因结构分析。采用ProtParam、WoLF PSORT和TMHMM Server v.2.0在线工具预测芝麻水通道蛋白的分子质量和等电点、亚细胞定位及跨膜结构域。通过芝麻、拟南芥和水稻AQP及马铃薯XIPs的蛋白多序列比对结果推测NPA基序、ar/R滤器及P1—P5的氨基酸残基。利用前期研究获得的转录组测序结果进行青枯雷尔氏菌诱导表达分析,并通过qRT-PCR技术对差异表达较为明显的12个芝麻AQP进行验证。【结果】系统分析鉴定了36个芝麻AQP家族基因,根据多序列比对及系统进化分析将其分为5个亚家族:13个质膜内在蛋白(PIP)、12个液胞膜内在蛋白(TIP)、8个类NOD26膜内在蛋白(NIP)、2个膜内在小分子碱性蛋白(SIP)和1个未知内在蛋白(XIP),其中有34个基因定位在12个连锁群上。同一亚家族成员在基因结构、蛋白序列、亚细胞定位预测及保守性氨基酸残基等方面都较为相似。青枯雷尔氏菌诱导表达分析显示,NIPs、SIPs和XIPs表达无明显变化,部分PIPs和TIPs能够响应青枯菌诱导。SiPIP1;2、SiPIP1;3、SiPIP2;3和SiPIP2;4受青枯菌诱导后表达上调,SiPIP1;3和SiPIP2;3为持续上调,而SiPIP1;2和SiPIP2;4的表达先下调后上调;与之相反,表达下调较为显著的有SiPIP1;4、SiPIP2;1、SiPIP2;6、SiTIP1;1、SiTIP1;3、SiTIP2;1及SiTIP2;2。上述差异�
