SWATH分析玉米蛋白质丰度、转录水平和组织多样性揭示其发育调控规律

上一篇 / 下一篇 2018-04-27 17:25:38/ 个人分类：蛋白质组

文献：An integrated analysis of protein abundance, transcript. level and tissue diversity to reveal developmental regulation of maize.Journal of Proteome Research（IF4.268）

关键技术：SWATH

SWATH分析玉米蛋白质丰度、转录水平和组织多样性揭示发育调控规律

研究背景

　　采用SWATH-MS技术定量分析四种玉米组织中的蛋白质：未成熟雌穗（V7阶段），未成熟雄穗（V8阶段），授粉后20天（20_DAP）幼胚和14日龄幼苗根，并进行蛋白组＆转录组关联分析，鉴定组织特异性高表达基因和蛋白质，了解组织结构和器官发生的调节机制。

研究内容

1）数据分析和蛋白质特性

　　合并来自四个组织样本的蛋白组学数据，作者通过SWATH共获得646,658个图谱。用不同的置信度阈值和错误发现率（FDR）来评估SWATH-MS的检测能力。发现。一定程度上较低阈值可鉴定到更多独特肽和蛋白质（图1A，1B）。在置信度0.85和FDR 0.05下，检测的多肽和蛋白质数量与使用更严格的阈值标准—置信度0.90和FDR 0.02相比，检测的多肽和蛋白数量稍有增加（图1A，1B），但差异蛋白质数量变化不大（图1D）。因此，最终选择置信度0.85和FDR 0.05阈值，因为该阈值能检测更多的蛋白质来扩增蛋白质组数据，而不会显著增加差异蛋白数量。共有117,184个unique谱图匹配到10,606个unique肽段，最后共鉴定到4,551个蛋白质（图1E）。蛋白质丰度在重复实验中重复性很好，相关系数在0.84-0.90之间（图S1）。在雌穗，雄穗，幼胚和幼根中分别鉴定到3916、3707、3702和2871种蛋白质，在这四种组织中共鉴定到2269种蛋白质（图1E）。大多数蛋白质（70％）被不少于2个unique谱图覆盖，主要由10至25个氨基酸组成（图S2A，S2B）。大约64％的蛋白质显示出> 5％的蛋白序列覆盖率度，84％的蛋白质的分子质量> 20kDa。

2）蛋白质组&转录组关联分析

　　基于Doreen Ware实验室的转录组数据，及qTeller中下载的RNA-seq数据，作者对蛋白质组和转录组数据进行了关联分析。先将每个组织中的蛋白质丰度与相应的mRNA水平进行比较。然后过滤掉那些低mRNA水平（RPKM<1）的蛋白质，总共保留了4,314种（94.8％）蛋白，其中雌穗3554个，雄穗3404个，幼胚3417个和幼根2370个（图2A）。四个组织中共鉴定到2045个蛋白，表明大量检测到的蛋白存在于各个组织中。此外，还鉴定到了组织特异性高度表达的蛋白，其中幼胚253个，幼根181个，雌穗123个和雄穗43个（图2A）。其次，尽管蛋白质丰度与组织中的mRNA水平呈正相关，但Pearson相关系数很低，从0.35（幼根）到0.43（幼胚）（图2B），这表明转录水平并不总是与蛋白质丰度一一对应。

　　去掉组织特异性高度表达的蛋白质，以变化倍数> 2倍和P值<0.01，FDR <0.05为筛选条件，共有3714个差异蛋白至少在两个组织中共同鉴定到。结果表明，在雌穗和雄穗之间鉴定到227个差异蛋白，在雌穗和幼根之间鉴定到799个差异蛋白（图2D-E）。差异蛋白质的数量反映了各组织形态和功能的差异。

　　为了探究蛋白质与转录本丰度相关性较低的原因，作者将这些差异蛋白质分成如下三个亚组：i）正相关亚组（PC）：mRNA水平与蛋白质丰度正相关；ii）低相关性亚组（LC）：编码差异蛋白的转录本与蛋白质丰度相关但不显著；iii）负相关亚组（NC）：mRNA水平与蛋白质丰度呈负相关。作者发现大约50％的差异蛋白质可以归类到PC亚组中（图2E），这表明过半的差异蛋白质的差异丰度主要由这些基因的表达差异所决定。然而，大约42％的差异蛋白，从31％到62％不等，归类到到LC亚组（图2E）。这种现象通常是由于RNA转录，加工和转换以及蛋白质翻译和转换的所导致的。4.9％至15％蛋白质归类到NC亚组（图2E），这表明转录组不能完全解释在组织或器官的不同功能和形态的分子基础上的蛋白质组，由此突出了蛋白质组学研究的重要性。

3）组织特异性高表达蛋白的功能分类

　　基于蛋白质丰度作者在四种组织中共鉴定了600种组织特异性高表达蛋白质（图3A）。在细胞的整个蛋白质组中，蛋白质的低覆盖率是众所周知的，在所有现有的蛋白质组技术，包括SWATH-MS都有明显的缺陷。这种低覆盖率可能导致对组织特异性高表达蛋白的数量的过高估计。因此，作者检索了来自于qTeller的全长转录组数据和RNA-seq数据中编码600个组织特异性高表达蛋白的基因的表达水平，发现大部分基因（84.3％）在多种组织中都有表达，只有94种蛋白编码基因(15.7%)显示出组织特异性表达模式(图3A)。

　　另外，26种蛋白质在雌穗和雄穗中共表达，称为花序特异性蛋白质。作者利用qRT-PCR交叉验证了18个基因的表达水平，15个为组织特异性表达基因，3个呈现出花序特异性表达模式（图3B，3D）。此外，作者还利用Western blot验证了2个幼胚特异性蛋白（GRMZM2G354013和GRMZM2G054916）（图3C，3E）。

　　作者鉴定了94个组织特异性蛋白，其中幼根60个，幼胚27个，雌穗5个和雄穗2个。幼根特异性蛋白质直接或间接地与氧化还原过程和脱氢有关，包括12个过氧化物酶超家族蛋白，3个还原酶和8个水解酶活性蛋白，表明对氧化胁迫的响应在幼根的发育和生长中起重要作用。在27个幼胚特异性蛋白中，有5种植物脂质转移蛋白/种子储存相关蛋白质，3种氧化还原酶活性蛋白质和2种已知的与光形态相关的蛋白质转录因子（AtPIF1和AtSPT转录因子）。7个花序特异性蛋白与花序结构和花器官的发育高度相关。值得注意的是，有18种蛋白质在雌穗和雄穗中都有很高的丰度，但其mRNA水平在四个组织中都呈现出的低水平或没有被检测到。

4）差异表达蛋白的功能分类

　　GO富集确定差异蛋白功能类别。较其他组织对，雌穗和雄穗间鉴定到的差异蛋白更少，这些差异蛋白主要富集到“细胞分裂素介导的信号调节”，“乙酰辅酶A代谢过程”和“分生组织生长”（图4A），表明细胞分裂素介导的信号在雌穗和雄穗间的细胞分裂活性存在差异。雄穗中UB3蛋白质丰度比雌穗中高1.7倍，表明UB蛋白质积累的增加可能支持雄穗长枝的长出。雌穗和幼胚间的差异蛋白主要富集到“脂质生物合成过程”，“分生组织起始”和“生殖结构发育”条目。雄穗和幼胚间的差异蛋白也富集到这些GO条目（图4B，4C）。差异蛋白中不同GO类别反映了花序和幼胚间生理和发育上的差异。

　　此外，雌穗/雄穗和幼根间的差异蛋白显著富集到“对无机物质的响应”和“有机酸生物合成过程”（图4D，4F）。幼胚和幼根间的差异蛋白主要富集到“有机酸代谢过程”，“对无机物质的反应”和“氧化还原”（p = 2.6 E-12）等。与“脂质生物合成过程”和“脂肪酸生物合成过程”有关的蛋白质在幼胚中显著富集，而涉及“响应无机物质”，“氧化还原”和“醇代谢过程”的蛋白质在幼根中显著富集。

文章小结

　　作者利用SWATH-MS技术在玉米4种组织中共鉴定到4551种蛋白质，与之前iTRAQ鉴定到的蛋白数目相当。基于蛋白质组和转录组关联分析，作者推断出蛋白质丰度与mRNA水平呈正相关，具有弱到中等的相关系数。

　　然而，一些特定组织中与功能或结构相关的一些关键蛋白质在转录水平上被时空调控，在4种组织中不同蛋白表现出不同的表达模式。此外，组织特异性高表达的蛋白和差异蛋白能富集到不同GO类别中。组织特异性高度表达的蛋白质亚组通过交叉实验验证可作为潜在生物标志物。这加深了对玉米组织和器官发育分子机制的理解，并为其生物标志物的发现提供了全面的线索。

备注：本文出于传递更多信息为目的解读该文献内容，不希望转载个人可与我联系，会立即删除处理。