Functional analysis of SAP11 gene in cassava
-
摘要:目的
木薯Manihot esculenta Crantz是全球热带地区重要的粮食作物和经济作物,在生长发育过程中极易遭受低温、干旱、盐碱等非生物胁迫而导致减产。胁迫相关蛋白(Stress-associated protein,SAP)是一类新型的A20/AN1锌指蛋白,在模式作物应对多种非生物胁迫过程中发挥重要作用。目前,SAP基因在木薯应对非生物胁迫中的生物学功能尚不明确。本研究旨在分析木薯SAP家族成员的蛋白结构特征和表达模式,以及MeSAP11的互作蛋白,为进一步解析该家族基因在木薯抗逆中的功能提供理论支撑。
方法利用生物信息学技术对木薯SAP家族成员的进化关系、蛋白基序信息以及时空表达模式开展系统分析。同时,通过qRT-PCR研究各基因成员在不同组织中的特异表达以及对不同非生物胁迫的响应。进一步运用酵母双杂交结合高通量测序技术获得与MeSAP11相互作用的蛋白及对应生物学通路。
结果木薯SAP基因家族共6个大类16个成员,该家族成员在木薯根部和叶片中表达量较高,部分家族成员的表达在低温和盐胁迫中显著上调,在干旱、钾饥饿和氮饥饿显著下调。MeSAP11的表达受不同胁迫条件的显著调控,亚细胞定位结果表明MeSAP11蛋白主要定位在细胞核。利用酵母双杂交筛库技术筛选到256个与MeSAP11互作的蛋白,KEGG分析表明这些互作基因主要参与蛋白泛素化降解、内质网蛋白质加工通路等途径,暗示MeSAP11可能通过上述通路发挥功能。
结论木薯SAP 家族大部分成员显著响应低温、干旱、高盐以及缺氮、缺钾胁迫,研究结果为进一步研究MeSAP11在木薯响应非生物胁迫过程中的功能并解析其调控网络奠定了基础。下一步将把MeSAP11基因列为调控非生物逆境变化的候选基因开展深入研究。
Abstract:ObjectiveCassava is an important food and economic crop in tropical regions worldwide, and it is highly susceptible to yield reduction due to non-biological stresses such as low temperature, drought, and salinity during its growth and development. Stress-associated proteins (SAPs) are a novel class of A20/AN1 zinc finger proteins that play important roles in the response of model crops to various non-biological stresses. The biological functions of SAPs in cassava’s response to non-biological stresses are not yet clear. This study aims to analyze the protein structure characteristics and expression patterns of the cassava SAPs, as well as the interacting proteins of MeSAP11, and provide theoretical support for further understanding the function of cassava SAPs in response to abiotic stresses.
MethodBioinformatics techniques were used to systematically analyze the evolutionary relationships, protein motif information, and spatiotemporal expression patterns of the cassava SAPs. Additionally, qRT-PCR was used to study the specific expression of each gene member in different tissues and their response to abiotic stresses. Furthermore, yeast two-hybrid combined with high-throughput sequencing technology was used to identify the proteins interacting with MeSAP11 and their corresponding biological pathways.
ResultThe cassava SAP gene family consisted of six major classes and 16 members. The expression levels of these family members were higher in cassava roots and leaves. The expressions of several family members were significantly up-regulated by low temperature and salt stress, and significantly down-regulated by drought, potassium starvation and nitrogen starvation. The expression of MeSAP11 was significantly regulated under different stress conditions, and subcellular localization results indicated that the MeSAP11 protein was mainly located in the nucleus. Using yeast two-hybrid screening, 256 proteins were identified to interact with MeSAP11, and KEGG analysis indicated that these interacting genes are mainly involved in protein ubiquitination degradation and endoplasmic reticulum protein processing pathways, suggesting that MeSAP11 may function through these pathways.
ConclusionThe majority of the cassava SAPs are significantly regulated by low temperature, drought, high salinity, nitrogen deficiency, and potassium deficiency stresses. These results lay a theoretical fundation for the function exploration of MeSAP11 in cassava responding to abiotic stress, and point out a direction for further study to unravel its regulatory network. MeSAP11 is identified as a candidate gene for further in-depth research on regulating non-biological stress changes.
-
Keywords:
- Cassava /
- Abiotic stress /
- MeSAP gene /
- Protein-protein interaction screening
-
木薯Manihot esculenta Crantz起源于南美洲[1],是全球第六大粮食作物,具有很高的经济价值和产量潜力[2];同时也是一种多用途作物,可用于食品加工和工业原料。随着世界人口不断增长,人们对粮食的需求量急剧增加[3],提高粮食作物产量尤为重要。近年来受多种因素影响,气候变化对作物生产力的影响越来越大[4]。极端温度、干旱和土壤盐碱化是植物经常遇到的不利环境条件[5-6],对植物分子、生化生理、形态建成等产生严重影响,从而影响作物的生长发育以及产量的形成[7]。为适应诸多逆境,植物必须形成相应的调节机制才能保护自己免受其侵袭,比如调整生命周期、改变细胞代谢、调节信号转导途径等[8],这些抵抗机制需要复杂的抗逆应答网络[9]。
SAP (Stress-associated proteins,SAPs)是一类含有A20/AN1锌指结构域的蛋白质,是参与多种非生物胁迫的新型E3泛素连接酶之一[10],可以调节植物中的胁迫信号[11],减少产量损失,在前人研究中已被证明可以赋予植物对多种非生物胁迫的耐受性,并已成为植物非生物胁迫研究的重要对象之一[12-13]。SAP基因在被子植物中高度保守,OsSAP1首先在植物中分离出来,它们的同系基因已在水稻、拟南芥、胡杨、土豆、棉花和芸苔油菜中鉴定出来[14-21]。迄今为止,来自不同物种的几个具有一个A20结构域和一个AN1结构域的SAP基因已被证实参与对多种非生物胁迫的耐受性。在水稻中,OsSAP1的表达是由不同的胁迫(寒冷、干旱、盐、重金属、ABA和损伤)诱导的,OsSAP1的过表达提高了对干旱、寒冷和盐胁迫的耐受性,并增强了其对转基因植物病原体的免疫力[10],过表达OsiSAP8的转基因植株也表现出对盐、干旱和寒冷胁迫的高耐受性[19]。在棉花中,AtSAP5的异位表达通过上调内源性胁迫响应基因的表达来提高棉花对干旱和热胁迫的耐受性[20]。过度表达海滨獐毛AlSAP通过维持转基因水稻的光合作用来增强对寒冷、干旱和盐胁迫的耐受性,并在干旱条件下提高水稻籽粒产量[21]。在拟南芥中过表达毛果杨PtSAP13可以显著提高耐盐性,并且转基因植株种子的萌发率比野生型更高[22]。目前,木薯SAP蛋白在非生物胁迫响应中的功能依旧未知。因此,本研究以模式植物拟南芥为参照,采用同源序列法识别木薯SAP同源基因,并对其进行生物信息学及时空表达分析,以期深入挖掘抗逆相关重要基因,为木薯分子设计育种提供重要的基因资源和理论支撑。
1. 材料与方法
1.1 试验材料
本试验材料为木薯栽培品种‘cv. 60444’,由中国热带农业科学院热带生物技术研究所功能基因组研究实验室提供。选取长势优良的‘cv. 60444’组培苗,对其根、茎、叶、芽、叶柄的组织部位进行收集用于RNA的抽提。再选取生长期和生长势一致的木薯‘cv. 60444’组培苗进行胁迫处理,分别取处理组和对照组的叶片进行液氮速冻后保存于−80 ℃用于RNA的抽提。
RNAprep Pure多糖多酚植物总RNA提取试剂盒(DP441)购自天根生化科技(北京)有限公司;反转录试剂盒(D2639A)和qRT-PCR所用的SYBR Premix Ex Taq Ⅱ Kit (DRR081A)购自莫奈生物科技有限公司;各种引物合成以及载体测序由北京擎科生物科技股份有限公司完成;2× Taq PCR反应体系(KT201-02)购自天根生化科技(北京)有限公司;PCR产物回收纯化试剂盒(DR02)、质粒提取试剂盒(PL03)购自艾德莱生物有限公司;大肠埃希菌DH5α、农杆菌LBA4404和酵母Y2HGold感受态细胞购自上海唯地生物技术有限公司;Clontech酵母双杂交试剂盒(630489)及酵母培养基购自宝生物有限公司;植物表达载体pGAMBIA1302-35S::GFP由中国热带农业科学院热带生物技术研究所改造保存。
1.2 试验方法
1.2.1 叶片总RNA的提取及cDNA的合成
根据RNAprep Pure多糖多酚植物总RNA提取试剂盒说明书,提取木薯总RNA,使用超微量分光光度计NanoDrop 2000检测总RNA的浓度和纯度,参照Fastking gDNA Dispelling RT SuperMix反转录试剂盒(TaKaRa)的说明书进行反转录,得到的cDNA用于后续试验。
1.2.2 MeSAP基因家族成员鉴定及理化性质分析
利用隐马尔可夫模型(Hidden Markov Model,HMM)[23]对其基因组进行搜索并分别通过NCBI-CDD、Pfam和SMART 3个公共数据库鉴定木薯MeSAP基因家族成员。从JGI植物基因组数据库(https://phytozome.jgi.doe. gov/pz/portal.html)[24]获得木薯MeSAP家族成员的基因ID、编码区(CDS)序列和氨基酸序列。通过TB-tools的Protein Paramter Calc功能[25],获得木薯MeSAP家族蛋白的相对分子质量、理论等电点、不稳定系数、蛋白疏水性和脂溶性。
1.2.3 MeSAP基因家族系统进化树分析
从JGI植物基因组数据库下载拟南芥、木薯的SAP家族蛋白序列。以木薯、拟南芥SAP家族蛋白的氨基酸序列为基础,通过MEGA 7.1的NJ (Neighbour-joining)法构建系统发育进化树[26],使用MEGA 7.1进行1 000次迭代bootstrap测试。
1.2.4 MeSAP家族蛋白保守结构分析和氨基酸序列比对
为进一步分析MeSAP的蛋白结构,使用MEME对SAP家族蛋白进行Motif分析,并利用TB-Tools工具进行绘图。从JGI植物基因组数据库下载木薯SAP蛋白序列。利用DNAMAN7软件对木薯SAP蛋白家族序列进行多序列比对。利用TB-tools软件对木薯SAP蛋白家族结构以及Motif基序进行绘制。
1.2.5 MeSAP基因家族的表达模式分析
为鉴定木薯MeSAP基因家族在木薯不同组织部位以及在木薯不同胁迫处理下的表达情况,以2月龄的木薯‘cv. 60444’组培无菌苗为材料,以MS液体培养基处理作为对照组,同时分别进行低温(4 ℃低温培养1 d)、干旱(200 g/L PEG6000溶液模拟干旱处理6 h)、高盐(200 mmol/L NaCl处理6 h)、缺氮(无氮MS液体培养基处理6 d)、缺钾(无钾MS液体培养基处理6 d)胁迫处理,收集对照组和处理组植株的幼嫩叶片和顶芽抽提RNA进行后续试验。为鉴定MeSAP基因家族成员的组织特异性表达特征,分别选取木薯幼苗的根、茎、顶芽、顶芽下方第1—3叶位的叶片及叶柄作为材料抽提RNA。以各材料的cDNA为模板,利用NCBI设计MeSAP基因家族成员的qRT-PCR引物,先通过半定量PCR确定MeSAP基因家族成员在木薯5个组织部位的表达,再通过qRT-PCR分析MeSAP基因家族成员在不同处理下的表达差异。
qRT-PCR用TIANGEN的Fastking gDNA dispelling RT SuperMix (4992251);在StepOne Real-Time PCR系统(Applied Biosystems)上使用TIANGEN的SYBR Premix Ex Taq Ⅱ Kit (fp202-02)对每个cDNA样本进行3次重复PCR;利用MeActin作为内参基因,MeSAP基因家族成员的表达量以2−△△Ct来计算。为分析MeSAP基因家族在不同组织部位的表达情况,通过GraphPad Prism 8工具绘图。使用SPSS 26.0的t检验功能分析数据间的差异显著性,以P<0.01表示差异极显著,P<0.05表示差异显著。
1.2.6 亚细胞定位
将MeSAP11基因cDNA全长片段构建到植物表达载体pGAMBIA1302-35S::GFP上,构建pGAMBIA1302-35S::GFP-MeSAP11植物表达载体,将测序正确的质粒转化至LBA4404农杆菌中,以pGAMBIA1302-35S::GFP载体为对照,经鉴定后的农杆菌按照烟草瞬时转化的方法注射到烟草叶片的下表皮,注射2 d后观察绿色荧光蛋白和mCherry的亚细胞定位。
1.2.7 酵母双杂交筛选MeSAP11互作蛋白
先将MeSAP11基因cDNA全长片段插入酵母表达载体pGBKT7中,将测序正确的质粒和pGADT7共转至酵母Y2H菌株中。分别挑取鉴定后的转化酵母菌株以9 g/L的NaCl溶液重悬,并稀释到10−1~10−4浓度梯度后分别点在SD/DDO和SD/QDO平板上,根据所稀释的菌液是否在SD/QDO培养基上长出菌斑判断MeSAP11是否具有自激活活性。
利用本试验已建好的木薯叶片cDNA文库筛选可能与MeSAP11互作的蛋白,筛选方法按酵母双杂交试剂盒说明书进行。将筛选到的阳性克隆进行菌落PCR测序鉴定。
2. 结果与分析
2.1 MeSAP基因家族成员鉴定及理化性质分析
为了鉴定木薯MeSAP基因家族成员,验证后得到16个木薯SAP基因家族成员,按照MeSAP家族成员染色体编号由小到大进行编序,并分别命名为MeSAP1~16。如表1所示,其编码的SAP家族蛋白由51~293个氨基酸组成,相对分子质量为6 000~32 400,理论等电点除MeSAP2外其余成员均大于7.9,多为碱性蛋白。蛋白不稳定系数小于40表明稳定性好,大于40提示蛋白可能不稳定;本研究发现不稳定蛋白均属于同组成员,其中MeSAP5和MeSAP13稳定性较好。家族成员的疏水指数为38~66,脂溶性系数为−0.8~−0.2,亲水系数大于0,说明16个家族成员均属于疏水性蛋白。
表 1 木薯MeSAP家族蛋白的理化性质Table 1. Physical and chemical properties of MeSAP family protein in cassava名称
Name氨基酸数量
Number of
amino acid相对分子质量
Relative molecular
weight理论等电点
Theoretical
PI不稳定系数
Instability
index蛋白疏水性
Aliphatic
index脂溶性系数
Coefficient of fat
solubilityMeSAP1 136 15 101.02 8.02 46.82 57.50 −0.744 MeSAP2 51 6 004.96 5.40 49.63 65.10 −0.698 MeSAP3 135 14 950.05 8.93 52.51 61.48 −0.583 MeSAP4 179 18 922.35 7.99 34.34 58.94 −0.380 MeSAP5 173 18 308.73 8.68 25.06 58.21 −0.387 MeSAP6 171 18 467.02 8.44 30.59 61.11 −0.480 MeSAP7 179 19 204.52 9.01 44.22 52.91 −0.618 MeSAP8 170 17 892.20 8.48 32.07 60.88 −0.284 MeSAP9 156 17 325.93 8.90 49.73 58.21 −0.440 MeSAP10 188 20 393.91 8.90 40.45 55.59 −0.545 MeSAP11 176 18 929.37 9.13 54.33 60.40 −0.431 MeSAP12 192 21 446.91 9.37 44.90 38.70 −0.782 MeSAP13 173 18 534.17 8.44 26.36 66.53 −0.349 MeSAP14 160 17 700.14 8.75 43.69 57.31 −0.589 MeSAP15 293 32 420.86 8.65 46.82 62.56 −0.597 MeSAP16 192 21 220.33 8.98 37.91 57.97 −0.588 2.2 MeSAP蛋白家族进化树分析
为了更好地了解MeSAP基因家族的起源和演化历程,利用MEGA7.0软件对木薯、拟南芥的SAP蛋白构建系统发育进化树。2个物种的SAP蛋白系统发育树分为6组(图1),A组分别为AtSAP3、AtSAP4、AtSAP6、MeSAP9、MeSAP11、MeSAP14,其中6个成员均包含AN1和A20结构域,从进化树中可以看出木薯MeSAP11与拟南芥当中的AtSAP3亲缘关系最近,在功能上可能存在相似。B组主要有5个蛋白,分别为AtSAP1、AtSAP7、AtSAP9、MeSAP4、MeSAP5,其中MeSAP4和MeSAP5的序列相似性为99%。C组多为MeSAP蛋白,有5种木薯蛋白MeSAP6、MeSAP7、MeSAP8、MeSAP10、MeSAP13和2种拟南芥蛋白AtSAP2、AtSAP5,其中MeSAP6和MeSAP13序列相似性为99%,MeSAP7和MeSAP10序列相似性为100%并聚类到AtSAP5,预示MeSAP7和MeSAP10与AtSAP5的功能可能相似。D组主要有AtSAP10、MeSAP1、MeSAP2、MeSAP3、MeSAP12,MeSAP3与组中其他成员序列相似性均较高。此外E组中仅有1个拟南芥蛋白AtSAP8,表明AtSAP8与其他成员关系较远。F组仅分布2种木薯蛋白MeSAP15和MeSAP16,且都不具有A20和AN1结构域,表明两者功能与其他成员存在一定差异。系统发育关系比较结果表明,聚集在一起或与其他物种聚集在一起的SAP成员可能具有相似的抗胁迫生物学功能。
2.3 MeSAP家族蛋白分析和氨基酸序列比对
Motif分析结果表明,16个MeSAP存在10个保守基序,基序长度位于0~200个氨基酸之间(图2A)。15个家族成员中均存在保守基序Motif 3,其中,除MeSAP2、MeSAP15、MeSAP16外其余家族成员均具有Motif 1和Motif 2保守基序,Motif 1为AN1锌指结构域,Motif 2为A20锌指结构域,可作为转录调控因子参与许多生理过程,并参与基因表达调控(图2B)。与其他蛋白相比,SAP11、SAP9和SAP14在起始位置独特保有Motif 6基序,在木薯中可能具有重要功能。从基因结构来看,在木薯SAP家族的16个基因中,外显子与内含子分布差异较大,SAP2、SAP3和SAP10无内含子(图2C),表明MeSAP家族在功能上可能存在差异。
通对16个MeSAP蛋白序列进行多序列比对(图2D),发现MeSAP家族大部分成员蛋白序列都在N端含有1个由22个氨基酸LCXNX2CGX2GX4MNLCSKC组成的A20结构域,在C端均含有1个由36个氨基酸CX2CX4GLTGXFXCX2CGX2FCX2HRX4HXCX组成的AN1保守结构域,MeSAP2、MeSAP15、MeSAP16蛋白序列不含有A20结构域和AN1保守结构域。多序列比对结果和Motif分析表明MeSAP基因家族成员序列高度保守。
2.4 MeSAP基因家族在木薯不同组织部位中特异性表达
采用qRT-PCR技术检测16个MeSAP基因在木薯根、茎、叶、芽、叶柄5个组织的表达谱,结果表明,MeSAP家族成员在茎中的表达水平均较低,与根组织表达量相比差异显著;部分基因在芽中的表达水平也不明显,如MeSAP3、MeSAP8、MeSAP9、MeSAP12、MeSAP14、MeSAP16;家族成员在根、叶和叶柄中表达较高,说明其可能在木薯营养生长期发挥重要作用;此外,MeSAP3在叶柄中几乎不表达,在芽中的表达也处于较低水平,可能不参与木薯顶部组织生长(图3)。
2.5 胁迫诱导下MeSAP基因家族表达模式
为更好地了解MeSAP基因家族对胁迫的响应,木薯分别经低温4 ℃、PEG模拟干旱、盐、氮饥饿、钾饥饿处理,采用qRT-PCR技术对MeSAP进行表达分析。结果表明,MeSAP1、MeSAP2、MeSAP3、MeSAP4、MeSAP11、MeSAP12在冷害和盐胁迫下显著上调,相反在干旱、钾饥饿或氮饥饿下表现出显著下调,MeSAP16在冷害和盐胁迫下也极显著响应。此外,部分成员MeSAP5、MeSAP8和MeSAP13在各处理条件下具有相对一致的表达模式,而MeSAP14、MeSAP15对氮饥饿响应较为显著(图4)。
2.6 MeSAP11蛋白亚细胞定位
为了进一步研究MeSAP基因家族成员的功能,本文选用受不同非生物逆境显著调控的MeSAP11作为候选基因,构建pGAMBIA1302-35S::GFP-MeSAP11载体,通过瞬时转化烟草叶片的方法对于MeSAP11蛋白的亚细胞定位进行研究。结果(图5)表明,对照组烟草叶片表皮细胞中均分布较强的绿色荧光,而35S::GFP-MeSAP11融合蛋白表达绿色荧光,和35S::H2B-mCherry蛋白共同定位在细胞核内。
2.7 MeSAP11互作蛋白的筛选
为深入解析MeSAP11蛋白的生物学功能,利用酵母双杂交筛库技术,对MeSAP11的互作蛋白进行筛选与鉴定。转录自激活试验表明,MeSAP11蛋白对酵母细胞没有毒性且无转录激活功能(图6)。进一步以pGBKT7-MeSAP11为诱饵,通过酵母双杂交系统从木薯逆境cDNA文库中筛选与MeSAP11互作的蛋白,将生长在四缺培养基上的酵母菌斑进行NGS测序,获得256个潜在的互作蛋白,表2列出了测序丰度较高的17个基因。点对点验证结果表明MeSAP11与HSF24、NAC47、SPL9存在蛋白互作(图6)。
表 2 MeSAP11酵母在cDNA文库的筛选结果Table 2. Screening of MeSAP11 yeast in cDNA library基因编号
Gene number注释
AnnotationManes.13G013400.1 小热休克蛋白HSP20 Manes.11G058600.1 聚泛素3 Polyubiquitin 3 Manes.17G035300.1 泛素样蛋白 Ubiquitin-like proteins Manes.16G032100.1 DNAJ同源家族C成员 DNAJ member C Manes.04G165900.1 GTPase激活蛋白 AGD11 Manes.01G042200.1/
Manes.13G124500.1热休克蛋白热应激转录因子
Heat shock protein HSF24Manes.07G019300.1 聚泛素4 Polyubiquitin 4 Manes.09G144100.1 转录因子 NAC47 Manes.09G032800.1 转录因子 SPL9 Manes.14G148600.1 I 类热休克蛋白
Class I heat shock proteinManes.09G042800.1 未知蛋白 Uncharacterized protein Manes.12G078200.1 WD重复蛋白 WD protein Manes.15G054800.1 转换因子 Translation factor Manes.05G204500.1 一种配子表达的跨膜蛋白HAP8 Manes.13G087500.1 蛋白激酶家族蛋白
Protein kinase family proteinsManes.09G036800.1 泛素样蛋白 Ubiquitin-like proteins 对测序获得的256个互作蛋白进行聚类分析,GO富集分析发现,叶绿体、细胞质、蛋白质结合以及对盐胁迫、热激、光敏和缺水响应相关通路蛋白显著富集(图7A)。KEGG结果发现蛋白泛素化降解、内质网蛋白质加工通路相关基因被显著富集(图7B)。综上,说明MeSAP11可能通过参与蛋白降解途径或与热激蛋白、转录因子协同参与调节蛋白质的合成或降解,以介导植物逆境响应的相关进程。
3. 讨论与结论
SAP蛋白在植物响应非生物胁迫中起着重要作用。水稻中的SAP蛋白OsSAP1已被证明通过调节内源性胁迫相关基因的表达,在过量表达时赋予植物非生物胁迫耐受性[10]。OsiSAP8通过A20和AN1锌指结构域的相互作用,在烟草和水稻中过表达后,使其在种子萌发及幼苗阶段对盐、干旱和冷胁迫的耐受性显著提高[19]。OsiSAP7通过充当泛素连接酶作为ABA和缺水应激信号的负调节因子[27]。在拟南芥中,AtSAP5与转录协同调节因子MBF1C共同作用建立拟南芥热应激的耐受性[11]。AtSAP5在棉花的异位表达通过上调内源性响应基因的表达来提高棉花对干旱和热胁迫的耐受性[21]。AtSAP10在拟南芥中过表达后,拟南芥对镍、锰、锌等重金属和高温胁迫具有很强的耐受性,与野生型相比,生长更加健康,生物量增加数倍,根更长[28]。毛果杨中的PtSAP13在拟南芥中过表达后,在盐胁迫处理下植株表现出比野生型更高的种子萌发率和更好的生长[22]。
通过对木薯SAP基因家族编码的蛋白结构及特点进行生物信息学分析发现,MeSAP蛋白家族也具有A20和AN1锌指结构域,另外SAP基因家族中不含或含有少量的内含子是这一类蛋白的一个重要特点[29]。通过亚细胞定位预测及验证发现本研究中的MeSAP11定位于细胞核,与其他物种中的SAPs相像,绝大部分均定位于细胞核。通过对文库筛选的互作蛋白进行分析发现多数为热休克蛋白、泛素酶以及转录因子等,其中热休克蛋白HSP20基因家族在马铃薯中表现为在热、盐胁迫和干旱胁迫等多种非生物胁迫中均显著上调表达[30]。在植物非生物胁迫响应调控中的功能机制研究中,转录因子一直备受关注,在酵母文库筛选中发现SPL和NAC是植物特有的一类转录因子,在植物抵抗非生物胁迫具有重要的功能,前人研究发现木薯MeSPL9通过负调控茉莉酸、脯氨酸和花青素等干旱逆境响应相关的代谢物关键合成酶基因的表达来影响木薯的抗旱性[31];过表达水稻NAC基因家族成员SNAC1基因能够提高转基因棉花的抗旱和抗盐的能力[32]。
目前,关于MeSAP参与木薯逆境胁迫抵御方面的研究尚不清楚,对木薯SAP家族基因的功能仍未知,本研究利用木薯基因组信息对木薯MeSAP家族进行分析鉴定,并利用生物信息学的方法对其系统进化、理化性质以及结构特征进行深入分析。同时结合qRT-PCR分析了MeSAP在木薯组织中的时空表达以及在逆境胁迫下的表达模式。此外,在木薯MeSAP中确定了1个在木薯根和叶片中优势表达的具有A20和AN1保守结构域的候选基因MeSAP11,通过亚细胞定位发现其主要定位于细胞核。并通过酵母双杂交系统从木薯叶片cDNA文库中筛选到部分可能与MeSAP11互作的候选蛋白,并初步验证了MeSAP11与其中的热激蛋白和胁迫应激方面的转录因子的互作关系。研究结果为进一步探究MeSAP11的功能及分子机制奠定了一定基础。
干旱和热胁迫会显著降低作物产量和质量,基因表达的调节是植物适应胁迫的关键因素[19],提高植物对干旱和热激耐受性的分子策略可以集中在增加特定酶和保护蛋白数量或操纵某些信号的转导途径[33],特定调节蛋白(如转录因子、蛋白激酶或泛素连接酶)的表达改变可用于调节下游基因的表达,以影响对应激的代谢反应[12]。通过前期的试验结果,对于是否确定MeSAP11与互作蛋白共同参与调控木薯的抗旱功能,还需要制备相应的转基因木薯来进一步验证其功能。
-
表 1 木薯MeSAP家族蛋白的理化性质
Table 1 Physical and chemical properties of MeSAP family protein in cassava
名称
Name氨基酸数量
Number of
amino acid相对分子质量
Relative molecular
weight理论等电点
Theoretical
PI不稳定系数
Instability
index蛋白疏水性
Aliphatic
index脂溶性系数
Coefficient of fat
solubilityMeSAP1 136 15 101.02 8.02 46.82 57.50 −0.744 MeSAP2 51 6 004.96 5.40 49.63 65.10 −0.698 MeSAP3 135 14 950.05 8.93 52.51 61.48 −0.583 MeSAP4 179 18 922.35 7.99 34.34 58.94 −0.380 MeSAP5 173 18 308.73 8.68 25.06 58.21 −0.387 MeSAP6 171 18 467.02 8.44 30.59 61.11 −0.480 MeSAP7 179 19 204.52 9.01 44.22 52.91 −0.618 MeSAP8 170 17 892.20 8.48 32.07 60.88 −0.284 MeSAP9 156 17 325.93 8.90 49.73 58.21 −0.440 MeSAP10 188 20 393.91 8.90 40.45 55.59 −0.545 MeSAP11 176 18 929.37 9.13 54.33 60.40 −0.431 MeSAP12 192 21 446.91 9.37 44.90 38.70 −0.782 MeSAP13 173 18 534.17 8.44 26.36 66.53 −0.349 MeSAP14 160 17 700.14 8.75 43.69 57.31 −0.589 MeSAP15 293 32 420.86 8.65 46.82 62.56 −0.597 MeSAP16 192 21 220.33 8.98 37.91 57.97 −0.588 表 2 MeSAP11酵母在cDNA文库的筛选结果
Table 2 Screening of MeSAP11 yeast in cDNA library
基因编号
Gene number注释
AnnotationManes.13G013400.1 小热休克蛋白HSP20 Manes.11G058600.1 聚泛素3 Polyubiquitin 3 Manes.17G035300.1 泛素样蛋白 Ubiquitin-like proteins Manes.16G032100.1 DNAJ同源家族C成员 DNAJ member C Manes.04G165900.1 GTPase激活蛋白 AGD11 Manes.01G042200.1/
Manes.13G124500.1热休克蛋白热应激转录因子
Heat shock protein HSF24Manes.07G019300.1 聚泛素4 Polyubiquitin 4 Manes.09G144100.1 转录因子 NAC47 Manes.09G032800.1 转录因子 SPL9 Manes.14G148600.1 I 类热休克蛋白
Class I heat shock proteinManes.09G042800.1 未知蛋白 Uncharacterized protein Manes.12G078200.1 WD重复蛋白 WD protein Manes.15G054800.1 转换因子 Translation factor Manes.05G204500.1 一种配子表达的跨膜蛋白HAP8 Manes.13G087500.1 蛋白激酶家族蛋白
Protein kinase family proteinsManes.09G036800.1 泛素样蛋白 Ubiquitin-like proteins -
[1] 曹升, 陈江枫, 黄富宇, 等. 广西木薯产业现状分析及其发展建议[J]. 南方农业学报, 2021, 52(6): 1468-1476. doi: 10.3969/j.issn.2095-1191.2021.06.005 [2] WANG W Q, FENG B X, XIAO J F, et al. Cassava genome from a wild ancestor to cultivated varieties[J]. Nature Communications, 2014, 5: 5110. doi: 10.1038/ncomms6110
[3] COBB J N, DECLERCK G, GREENBERG A, et al. Next-generation phenotyping: Requirements and strategies for enhancing our understanding of genotype-phenotype relationships and its relevance to crop improvement[J]. Theoretical and Applied Genetics, 2013, 126(4): 867-887. doi: 10.1007/s00122-013-2066-0
[4] 李智博, 董世满, 李淑霞, 等. 木薯SR45亚家族基因鉴定及表达[J]. 华南农业大学学报, 2022, 43(5): 20-28. doi: 10.7671/j.issn.1001-411X.202111032 [5] GONG Z, XIONG L, SHI H, et al. Plant abiotic stress response and nutrient use efficiency[J]. Science China-Life Sciences, 2020, 63(5): 635-674. doi: 10.1007/s11427-020-1683-x
[6] LANDI S, HAUSMAN J F, GUERRIERO G, et al. Poaceae vs. abiotic stress: Focus on drought and salt stress, recent insights and perspectives[J]. Frontiers in Plant Science, 2017, 8: 1214. doi: 10.3389/fpls.2017.01214
[7] 朱健康, 倪建平. 植物非生物胁迫信号转导及应答[J]. 中国稻米, 2016, 22(6): 52-60. doi: 10.3969/j.issn.1006-8082.2016.06.012 [8] MUKHOPADHYAY A, VIJ S, TYAJI A. Overexpression of a zinc-finger protein gene from rice confers tolerance to cold, dehydration, and salt stress in transgenic tobacco[J]. Proceedings of the National Academy of Sciences of the United States of America, 2004, 101(16): 6309-6314.
[9] 崔江慧, 杨溥原, 常金华. 高粱GRF基因家族鉴定及在非生物胁迫下的表达分析[J]. 中国农业科技导报, 2021, 23(4): 37-46. [10] KOTHARI K S, DANSANA P K, GIRI J, et al. Rice stress associated protein 1 (OsSAP1) interacts with aminotransferase (OsAMTR1) and pathogenesis-related 1a protein (OsSCP) and regulates abiotic stress responses[J]. Frontiers in Plant Science, 2016, 7: 1057.
[11] KANG M, ABDELMAGEED H, LEE S, et al. AtMBP-1, an alternative translation product of LOS2, affects abscisic acid responses and is modulated by the E3 ubiquitin ligase AtSAP5[J]. The Plant Journal, 2013, 76(3): 481-93. doi: 10.1111/tpj.12312
[12] WANG W X, VINOCUR B, ALTMAN A. Plant responses to drought, salinity and extreme temperatures: Towards genetic engineering for stress tolerance[J]. Planta Medica, 2003, 218(1): 1-14.
[13] DANSANA P K, KOTHARI K S, VIJ S, et al. OsiSAP1 overexpression improves water-deficit stress tolerance in transgenic rice by affecting expression of endogenous stress-related genes[J]. Plant Cell Reports, 2014, 33(9): 1425-1440. doi: 10.1007/s00299-014-1626-3
[14] HUANG J, WANG M M, JIANG Y, et al. Expression analysis of rice A20/AN1-type zinc finger genes and characterization of ZFP177 that contributes to temperature stress tolerance[J]. Genes, 2008, 420(2): 135-144.
[15] VIJ S, TYAGI A. Genome-wide analysis of the stress associated protein (SAP) gene family containing A20/AN1 zinc-finger(s) in rice and their phylogenetic relationship with Arabidopsis[J]. Molecular Genetics And Genomics, 2006, 276(6): 565-575. doi: 10.1007/s00438-006-0165-1
[16] JIA H X, LI J B, ZHANG J, et al. Genome-wide survey and expression analysis of the stress-associated protein gene family in desert poplar, Populus euphratica[J]. Tree Genetics & Genomes, 2016, 12(4): 78.
[17] BILLAH S A, KHAN N Z, ALI W, et al. Genome-wide in silico identification and characterization of the stress associated protein (SAP) gene family encoding A20/AN1 zinc-finger proteins in potato (Solanum tuberosum L.)[J]. PLoS One, 2022, 17(8): e0273416. doi: 10.1371/journal.pone.0273416
[18] HE X, XIE S, XIE P, et al. Genome-wide identification of stress-associated proteins (SAP) with A20/AN1 zinc finger domains associated with abiotic stresses responses in Brassica napus[J]. Environmental and Experimental Botany, 2019, 165: 108-119. doi: 10.1016/j.envexpbot.2019.05.007
[19] KANNEGANTI V, GUPTA A K. Overexpression of OsiSAP8, a member of stress associated protein (SAP) gene family of rice confers tolerance to salt, drought and cold stress in transgenic tobacco and rice[J]. Plant Molecular Biology Reporter, 2008, 66(5): 445-462. doi: 10.1007/s11103-007-9284-2
[20] HOZAIN M, ABDELMAGEED H, LEE J, et al. Expression of AtSAP5 in cotton up-regulates putative stress-responsive genes and improves the tolerance to rapidly developing water deficit and moderate heat stress[J]. Journal of Plant Physiology, 2012, 169(13): 1261-1270. doi: 10.1016/j.jplph.2012.04.007
[21] BEN S R, MEYNARD D, BEN R W, et al. The promoter of the AlSAP gene from the halophyte grass Aeluropus littoralis directs a stress-inducible expression pattern in transgenic rice plants[J]. Plant Cell Reports, 2015, 34(10): 1791-1806. doi: 10.1007/s00299-015-1825-6
[22] LI J, SUN P, XIA Y, et al. A stress-associated protein, PtSAP13, from Populus trichocarpa provides tolerance to salt stress[J]. International Journal of Molecular Sciences, 2019, 20(22): 5782. doi: 10.3390/ijms20225782
[23] FUENTES-BEALS C, VALDÉS-JIMÉNEZ A, RIADI G. Hidden Markov Modeling with HMMTeacher[J]. PLoS Computational Biology, 2022, 18(2): e1009703. doi: 10.1371/journal.pcbi.1009703
[24] GOODSTEIN D M, SHU S, HOWSON R, et al. Phytozome: A comparative platform for green plant genomics[J]. Nucleic Acids Research, 2012, 40(D1): D1178-D1186. doi: 10.1093/nar/gkr944
[25] CHEN C J, CHEN H, ZHANG Y, et al. TBtools: An integrative tookit developed for interactive analyses of big biological data[J]. Molecular Plant Pathology, 2020, 13(8): 1194-1202.
[26] KUMAR S, STECHER G, TAMURA K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets[J]. Molecular Biology and Evolution, 2016, 33(7): 1870-1874. doi: 10.1093/molbev/msw054
[27] SHARMA G, GIRI J, TYAGI A K. Rice OsiSAP7 negatively regulates ABA stress signalling and imparts sensitivity to water-deficit stress in Arabidopsis[J]. Plant Science, 2015, 237: 80-92. doi: 10.1016/j.plantsci.2015.05.011
[28] DIXIT A R, DHANKHER O P. A novel stress-associated protein ‘AtSAP10’ from Arabidopsis thaliana confers tolerance to nickel, manganese, zinc, and high temperature stress[J]. PLoS One, 2011, 6(6): e20921. doi: 10.1371/journal.pone.0020921
[29] 张子奇, 李可, 陈银华, 等. 木薯MeSAP13基因的克隆及其抗细菌性枯萎病功能鉴定[J]. 热带作物学报, 2022, 43(10): 1981-1988. [30] ZHAO P, WANG D D, WANG R Q, et al. Genome-wide analysis of the potato Hsp20 gene family: Identification, genomic organization and expression profiles in response to heat stress[J]. BMC Genomics, 2018, 19(1): 61. doi: 10.1186/s12864-018-4443-1
[31] LI S X, CHENG Z H, LI Z B, et al. MeSPL9 attenuates drought resistance by regulating JA signaling and protectant metabolite contents in cassava[J]. Theoretical and Applied Genetics, 2022, 135(3): 817-832. doi: 10.1007/s00122-021-04000-z
[32] LIU G Z, LI X L, JIN S X, et al. Overexpression of rice NAC gene SNAC1 improves drought and salt tolerance by enhancing root development and reducing transpiration rate in transgenic cotton[J]. PLoS One, 2014, 9(1): e86895. doi: 10.1371/journal.pone.0086895
[33] KIM G D, CHO Y H, YOO S D. Regulatory functions of evolutionarily conserved AN1/A20-like Zinc finger family proteins in Arabidopsis stress responses under high temperature[J]. Biochemical and Biophysical Research Communications, 2015, 457(2): 213-220. doi: 10.1016/j.bbrc.2014.12.090
-
期刊类型引用(0)
其他类型引用(1)