VPS28通过泛素化信号通路调控奶牛乳腺上皮细胞中乳蛋白的合成

任 宇，燕文荃，丁玥竹，扈孟雪，周杰珑，吴培福，陈粉粉，刘莉莉

(西南林业大学生命科学学院，昆明650224)

牛奶中乳蛋白的含量约为3%～5%，是牛奶的主要成分之一。因其几乎含有人体所需的全部必需氨基酸并易于吸收，满足机体生长和维持日常生命活动的需要，被认为是评价牛奶价值的重要指标之一。如何提高牛奶中乳蛋白水平并改善牛奶中乳蛋白致敏性一直是学术界的焦点，也是生产者、消费者和乳品加工企业关注的热点。因此，通过分子遗传育种等现代手段改善牛奶中乳蛋白的品质是奶牛育种的重要目标之一。

乳蛋白主要是由酪蛋白和乳清蛋白组成，是由乳腺上皮细胞通过细胞膜上氨基酸转运载体吸收转运血液中游离氨基酸而合成。乳蛋白的合成主要在粗面内质网的核糖体上开始，然后由信号肽引导进入内质网中，并在内质网和高尔基体内进行磷酸化和糖基化等化学修饰，再由分泌泡转运至上皮细胞顶膜，通过胞吐方式释放到细胞外的腺泡腔。乳蛋白的合成受到多个信号通路及多种激素相互影响，如：氨基酸转运、蛋白酪氨酸激酶2-信号转导子和转录激活因子5(janus tyrosine kinase 2-signal transducer and activator of transcription 5, Jak2-Stat5)和雷帕霉素靶蛋白(mammalian target of rapamycin, mTOR)等信号通路，胰岛素、催乳素和生长激素等。由此可见，各通路中的信号传递和激素进入细胞后的信号传递过程至关重要。

28(vacuolar protein sorting 28)是一种液泡蛋白分选基因，是内体蛋白分选转运装置Ⅰ(endosomal sorting complex required for transport-Ⅰ，ESCRT-Ⅰ)的亚单位，与ESCRTs的稳定性有关。泛素化是细胞中以单泛素蛋白和多泛素蛋白链为信号分子，将泛素蛋白结合到底物蛋白特定位点的一种翻译后修饰过程。ESCRTs不仅能识别泛素化膜蛋白，将其转运至溶酶体进行降解，还参与蛋白酶体对胞质中多泛素化蛋白的降解。在课题组前期研究中，发现28基因与奶牛乳蛋白率和乳蛋白量均存在极显著关联(<0.01)，并发现该基因在奶牛的乳腺组织中特异性高表达，但28基因对乳蛋白的调控尚无相关研究。因此，本研究以28基因为奶牛乳蛋白性状的重要候选基因，采用RNAi干扰技术和iTRAQ比较蛋白质组学，分析并试图阐明28基因调控乳蛋白合成的分子机制，为奶牛产奶性状的分子育种研究提供理论基础。

1 材料与方法

1.1 材料

奶牛原代乳腺上皮细胞(BMECs)为原实验室冻存；DMEM/F12培养基、DMEM培养基、Opti-MEM培养基、胎牛血清、双抗(penicillin-streptomycin solution)购自美国Gibco公司；TRIzol购自美国Invitrogen公司；转染试剂X-treme GENE siRNA Transfection Reagent、SYBR Green Mix购自美国Roche公司；溶酶体抑制剂氯喹CQ(chloroquine，C668)购自美国Sigma-Aldrich公司；蛋白酶体抑制剂环氧霉素Epox(epoxomicin，BU-4061T)购自美国MCE公司；泛素抗体(ubiquition，sc53509)购自美国Santa Cruz公司；β-actin抗体(66009-Ig)购自美国Proteintech公司；HRS标记二抗(7076S)购自美国Cell Signaling Technology公司；蛋白酶体活性检测试剂盒(Proteasome-GloChymo-trypsin-Like, Trypsin-Like and Caspase-Like Cell-Based Assays，G1180)购自美国Promega公司。

1.2 细胞培养

将BMECs培养于DMEM/F12全培养基中(含有10%胎牛血清，100 U·mL青霉素-链霉素)，并将其置于温度37 ℃，CO浓度5%的细胞培养箱中。当细胞汇合度达到100%时，用0.25%胰酶进行消化并将BMECs以2.5×10个·孔传代接种至6孔板中，以备后续试验使用。

1.3 BMECs中VPS28基因的敲降

培养BMECs 24 h，当细胞汇合度达80%时使用siRNA串联片段对28基因进行敲降，首先，用150 μL Opti-MEM I Reduced Serum Medium分别稀释2 μg siRNA和20 μL X-treme GENE siRNA Transfection Reagent；然后将二者轻柔混匀并静置20 min；最后，缓慢滴加入细胞培养液中。6 h后更换为完全培养基，培养72 h后收集细胞。

1.4 VPS28基因及乳蛋白合成相关基因的表达量检测

结合KEGG网站(KEGG-Table of Contents, http://www.genome.jp/kegg/kegg2.html)，本研究选择乳蛋白合成、泛素化-溶酶体和泛素化-蛋白酶体通路中11个基因作为候选基因，采用qRT-PCR方法检测这11个相关基因mRNA表达量的变化(具体基因和引物信息见表1)。以为内参基因，采用比较Ct值法即2(△△Ct =目的基因△Ct值-内参基因△Ct值)表示目的基因的相对表达量。采用美国ABI公司SYBR Green I进行定量表达检测。反应体系及反应程序参照SYBR Green I说明书，在ABI 7500荧光定量PCR仪上进行扩增。

表1 目的基因与引物信息

1.5 泛素蛋白表达水平的检测

使用siRNA串联片段对BMECs中28基因进行敲降后，收集细胞提取总蛋白，利用Western blot技术检测敲降前后BMECs中泛素蛋白的表达水平。配制10%聚丙烯酰胺凝胶，按40 μg·孔上样，70 V电泳0.5 h至分离胶处，调整电压为120 V 继续电泳1 h。采用湿转法转膜，200 mA 1 h。采用5%脱脂奶粉封闭1 h后，分别加入泛素抗体(1∶200) 和β-actin抗体(1∶1 000)4 ℃孵育过夜，1×TBST洗涤3次，5 min·次。二抗(1∶2 000)摇晃孵育1 h，1×TBST洗涤3次，5 min·次，用ECL发光液显影，用Image J软件进行灰度值分析。

1.6 蛋白酶体活性的检测

利用蛋白酶体活性检测试剂盒检测28敲降组和未敲降对照组中胰蛋白酶样、糜蛋白酶样和半胱天冬酶样蛋白酶活性。配制荧光素酶检测液，室温孵育30 min。将准备好的细胞悬液加入白壁96孔板中，并加入等体积的检测液，室温孵育10 min 后检测相对荧光信号值。

1.7 BMECs中溶酶体和蛋白酶体的抑制

取10 mg溶酶体抑制剂氯喹CQ，加入0.387 7 mL超纯水，充分混匀配置为50 mmol·LCQ处理液，置于-20 ℃保存；取100 μg蛋白酶体抑制剂Epox，溶于180.27 μL DMSO中，充分混匀配置为1 mmol·L的Epox处理液，置于-20 ℃保存。培养BMECs 24 h，当6孔板中细胞汇合度达80%时，在不同孔中分别加入2 μL CQ或20 μL Epox孵育细胞，24 h收集细胞。

1.8 比较蛋白质组学的分析

使用siRNA串联片段对BMECs28基因进行敲降后，收集细胞提取总蛋白，利用iTRAQ比较蛋白质组学分析敲降前后BMECs的差异表达蛋白。利用在线软件DAVID(https://david.ncifcrf.gov/)对分析得到的差异蛋白(≤0.8或≥ 1.25倍)进行富集分析，并利用在线软件STRING(https://www.string-db.org/)对差异蛋白进行互作分析。

1.9 数据分析

所有试验数据均采用3次独立重复试验结果，以“均数±标准差”表示，两组数据差异使用SPSS 21.0软件进行独立样本检验统计分析，<0.05 时认为具有统计学意义。

2 结 果

2.1 敲降VPS28基因后乳蛋白相关基因、泛素蛋白的表达水平及蛋白酶体活性

采用串联片段的siRNA对BMECs中28基因进行敲降，在成功敲降72%后(图1A、1B)，检测了11个候选基因的表达趋势，结果如图1C所示，酪蛋白基因1S1、2、3、核糖体蛋白8、泛素蛋白基因及蛋白酶体相关基因3、5表达水平均显著上调(<0.05)，12表达水平显著下调(<0.05)，13、2L差异不显著；泛素蛋白表达水平及蛋白酶体活性如图1D、1E所示，在敲降28基因后，泛素蛋白显著上调(<0.05)，蛋白酶体活性下调。结果提示，28基因可能通过泛素-蛋白酶体及泛素-溶酶体信号通路改变细胞中核糖体合成乳蛋白的水平。

A. 空白细胞(10×)；B. BMECs中转染siRNAs敲降VPS28基因(10×)；C. 敲降VPS28基因后BMECs中乳蛋白相关基因的相对表达水平；D. 敲降VPS28基因后BMECs中泛素蛋白的相对表达水平；E. 敲降VPS28基因后BMECs中蛋白酶体相对活性A. Blank cells (10×); B. Knockdown of VPS28 gene by transfection of siRNAs in BMECs (10×); C. Relative expression levels of milk protein-related genes in BMECs after knockdown of the VPS28 gene; D. Relative expression levels of ubiquitin proteins in BMECs after knockdown of the VPS28 gene; E. Relative proteasome activity in BMECs after knockdown of the VPS28 gene图1 敲降BMECs中VPS28基因及检测乳蛋白相关基因、泛素蛋白的表达水平及蛋白酶体活性Fig.1 The knockdown of the VPS28 gene with siRNAs in BMECs and detection of the mRNA expression levels of milk protein related genes, ubiquitin protein, and proteasome activity

2.2 抑制蛋白酶体和溶酶体活性后乳蛋白相关基因表达水平

为验证28基因可以通过泛素-溶酶体及泛素-蛋白酶体信号通路改变细胞中乳蛋白的合成，本研究通过蛋白酶体抑制剂Epox及溶酶体抑制剂CQ抑制细胞中蛋白酶体及溶酶体活性，并检测乳蛋白及核糖体相关基因的表达水平。结果如图2所示，抑制细胞中蛋白酶体活性后，细胞中酪蛋白基因11、2、3的相对表达量显著上调(<0.05)，核糖体蛋白13的相对表达量显著下调(<0.05)；抑制细胞中溶酶体活性后酪蛋白相关基因表达不显著，核糖体相关基因8显著上调(<0.05)。结果表明，泛素-溶酶体及泛素-蛋白酶体信号通路确实可以改变核糖体活性及乳蛋白的合成。

图2 抑制BMECs中蛋白酶体及溶酶体活性后乳蛋白相关基因的相对表达水平Fig.2 Relative mRNA expression levels of milk protein-related genes after inhibition of proteasome and lysosome activities in BMECs

2.3 敲降BMECs中VPS28基因后的蛋白谱分析

为进一步确定28基因通过泛素化途径调控BMECs中乳蛋白合成的调控水平及通路，本研究利用iTRAQ技术对敲降28基因后的BMECs进行了蛋白谱分析。通过与未敲降的对照组相比，共发现129个差异表达蛋白，其中75个显著下调(下调倍数≤0.8倍)，54个显著上调(上调倍数≥1.25倍)。对差异蛋白进行GO分析，如图3A、B所示，在生物学过程中，下调蛋白主要富集在细胞质翻译(GO:0002181)、高密度脂蛋白粒子组装(GO:0043691)、胆固醇反向运输(GO:0043691)及翻译(GO:0006412)，上调蛋白主要富集在基因表达正向调控(GO:0010628)及过时氧化还原过程(GO:0055114)；在细胞组分中，下调蛋白主要富集在乳糜微粒(GO:0042627)、胞质大核糖体亚基(GO:0022625)、粘着斑(GO:0005925)、细胞膜(GO:0016020)及细胞核(GO:0005730)，上调蛋白主要富集在细胞外空间(GO:0005615)、浆膜的重要组成部分(GO:0005887)、膜结构组分(GO:0016021)、胞外外泌体(GO:0070062)及内质网膜组分(GO:0005789)；在分子功能中，下调蛋白主要富集在核糖体结构组成(GO:0003735)、rRNA结合(GO:0019843)、磷脂酰胆碱固醇-甾醇O-酰基转移酶活化剂(GO:0060228)、胆固醇转运体活性(GO:0017127)、胆固醇酯化的正向调节(GO:0010873)及胆固醇结合(GO:0015485)，上调蛋白主要富集在核苷酸结合(GO:0000166)、ATP结合(GO:0005524)、poly-A-RNA结合(GO:0044822)及金属离子结合(GO:0005506)。对差异蛋白进行KEGG通路分析，如图3C所示，下调蛋白主要富集在核糖体(bta03010)、溶酶体通路(bta04142)、剪切体(bta03040)等信号通路，上调蛋白主要富集在内质网的蛋白质加工(bta04141)、加压素调控的水重吸收(bta04962)、RNA转运(bta04962)和阿尔兹海默症(bta05010)等通路中。对差异蛋白进行STRING互作分析，如图3D所示，差异蛋白主要集中在核糖体相关蛋白的互作上。结果表明，在BMECs中敲降28基因可以显著降低核糖体中蛋白质的翻译及溶酶体的功能，同时可以提高细胞中基因的表达转录水平、内质网的蛋白质加工合成及泛素化介导的蛋白质降解过程。

A. 上调差异表达蛋白的GO功能注释；B. 下调差异表达蛋白的GO功能注释；C. 差异表达蛋白的KEGG通路分析；D. 差异表达蛋白互作网络分析A. GO functional annotation of up-regulated differentially expressed proteins; B. GO functional annotation of down-regulated differentially expressed proteins; C. Analysis of the KEGG pathway of differentially expressed proteins; D. Analysis of interaction networks of differentially expressed proteins图3 VPS28干扰前后BMECs中差异表达蛋白的富集分析Fig.3 Enrichment analysis of differentially expressed proteins in BMECs before and after VPS28 interference

3 讨 论

本研究以28基因为影响奶牛乳蛋白性状的重要功能候选基因，在细胞水平对该基因进行功能研究，试图解释28基因调控乳蛋白合成的重要作用，为该基因的功能研究及奶牛乳蛋白性状相关分子遗传标记的筛选提供重要理论依据。

VPS28是真核细胞中内体蛋白分选转运装置ESCRTs的亚单位，属于ESCRT-Ⅰ的顶端帽子结构。VPS28可以通过GLUE结构域直接识别膜蛋白上的泛素蛋白并与之相结合，然后召集ESCRT-0、ESCRT-Ⅱ和ESCRT-Ⅲ组合形成超级复合体ESCRTs，再通过多囊体将泛素化的膜蛋白运输至溶酶体中进行降解，由此构成泛素-溶酶体通路，完成细胞中膜蛋白的泛素化降解过程。同时，研究发现，泛素-溶酶体通路可以通过泛素化影响泛素-蛋白酶体通路的活性，对细胞中泛素-蛋白酶体通路介导的胞质泛素化蛋白降解起到一定调控作用，尤其是当细胞中蛋白酶体活性受限后，泛素-溶酶体信号通路可被激活以补偿、维持细胞中蛋白质的稳态。泛素化是由泛素激活酶(E1)、泛素结合酶(E2)和泛素连接酶(E3)介导的，将泛素蛋白结合到底物蛋白特定位点的一种翻译后修饰过程，可以控制细胞中蛋白质的半衰期和表达水平。泛素化不仅可以作为信号传递过程中的重要调节蛋白，调控包括转录调控、细胞应答及新陈代谢等多种细胞活动，还可以直接识别并标记内质网中错误折叠的蛋白质，使其被有效降解。Mercier和Gaye在1982年研究发现,乳蛋白的合成主要在粗面内质网的核糖体开始，然后由信号肽引导进入内质网中，并在内质网和高尔基体内进行磷酸化和糖基化等化学修饰，再由分泌泡转运至上皮细胞顶膜，最后通过胞吐方式释放到腺泡腔。28基因可能通过泛素-溶酶体和泛素-蛋白酶体信号通路影响核糖体的合成过程和内质网、高尔基体的加工过程，进而影响奶牛乳腺上皮细胞中乳蛋白的合成。Liu和Zhang等指出，血液中游离的氨基酸、脂肪酸、甘油等物质通过主动或被动运输进入乳腺上皮细胞，为保证进入细胞内的物质正确，这些物质需要通过ESCRTs的分选。综合以上研究推测，泛素化信号通路在调控奶牛乳腺上皮细胞合成乳蛋白和乳脂的作用机制相同。

为验证28基因可能通过泛素-溶酶体和泛素-蛋白酶体信号通路影响核糖体中乳蛋白的合成，本研究首先利用RNAi干扰技术对BMECs中的28基因进行有效敲降，然后检测酪蛋白、核糖体和蛋白酶体等相关基因的表达水平以及泛素蛋白的表达水平和蛋白酶体活性。结果发现，敲降28基因后可显著改变细胞中酪蛋白和核糖体相关基因的表达水平，并可显著提高泛素蛋白的表达水平和降低蛋白酶体的活性，与前期研究结果相同，这表明28基因被敲降后可以通过下调蛋白酶体活性而提高泛素蛋白在细胞内的蓄积，这也提示了28基因可能是通过泛素化影响了蛋白酶体活性，进而调控核糖体中乳蛋白的合成。为验证这一结果，本研究利用溶酶体抑制剂和蛋白酶体抑制剂抑制BMECs中溶酶体和蛋白酶体的活性，再次检测酪蛋白和核糖体相关基因的表达，结果发现了相同的趋势，说明28基因确实可以通过泛素-溶酶体和泛素-蛋白酶体信号通路影响BMECs中乳蛋白的合成。

为进一步确定28基因通过泛素化作用对BMECs中乳蛋白合成的调控水平及调控机制，本研究利用iTRAQ技术对敲降前后的BMECs进行了比较蛋白质组学分析。GO分析表明，下调蛋白主要集中在核糖体翻译、细胞质翻译和转运等过程中，而上调蛋白主要分布在细胞膜、细胞质以及细胞外，具有多种结合活性。KEGG分析表明，差异表达蛋白富集在核糖体、溶酶体、剪接体、内质网的蛋白质加工、加压素调控的水重吸收和RNA转运等相关通路中，其中，下调蛋白主要富集在核糖体和溶酶体。核糖体不仅是细胞中翻译蛋白的主要场所，还积极参与了蛋白质的折叠过程，对于各种功能性蛋白质的合成至关重要。溶酶体是细胞中酸性的单层膜细胞器，不仅含有多种水解酶，用于降解细胞中的膜蛋白质、脂类、核酸和多糖等物质，还参与了细胞的胞吞、胞吐和外泌体分泌等过程。当细胞中溶酶体的功能受到抑制可导致线粒体应激，由此引发细胞中蛋白质失衡，同时显著降低细胞的代谢活性，此外，有研究表明，溶酶体可以代偿性发挥蛋白酶体的泛素化胞质蛋白的降解功能，维持细胞中的蛋白质平衡和内环境的稳态，因此，溶酶体可能在游离氨基酸向乳腺上皮细胞转运、乳蛋白合成和顶浆分泌过程中均起到重要调节和决定性作用。RNA剪接是真核生物基因表达调控的重要环节之一，是“中心法则”的关键步骤之一，而负责执行这一过程的是细胞核内一个巨大的且高度动态变化的分子机器—剪接体(spliceosome)。剪接体在真核生物进化中极为保守，对于真核生物维持正常的生命活动至关重要。一个基因转录出的前体mRNA可以通过RNA剪接形成若干种mRNA，极大地丰富真核生物蛋白质组的多样性，其在乳蛋白的功能性蛋白质合成中起到重要作用。因此，本研究表明，敲降28基因可能通过下调核糖体、溶酶体和剪接体相关蛋白的表达，减少乳蛋白的合成和分泌。

加压素主要调控细胞中水的重吸收，可以调控细胞器及细胞膜的通透性，参与乳蛋白的分泌及胞内运输等过程。RNA转运指RNA的核浆穿梭过程，RNA不同的亚细胞定位决定了其不同的生物学功能。研究表明，RNA转运参与了mRNA翻译、蛋白合成和基因调控等多种生物学过程，因此该通路与乳蛋白的合成调控密切相关。内质网是由单层膜构成的囊状、泡状和管状系统。内质网的蛋白质加工是影响内质网中蛋白质折叠的途径，正确折叠的蛋白质被包装成运输囊泡，转运至高尔基体，错误折叠的蛋白质则被泛素分子标记，进一步转运到溶酶体或蛋白酶体进行降解。因此，本研究在敲降28基因后，发现BMECs中加压素介导水分重吸收、内质网中蛋白质加工、RNA转运及泛素化介导蛋白质降解等通路的相关蛋白表达上调，与前人研究结果一致。但28基因对细胞中线粒体活性及细胞中所有代谢产物的具体影响还需要后续深入研究。

4 结 论

本研究以28基因为奶牛乳蛋白性状的重要候选基因，采用RNAi干扰技术和iTRAQ比较蛋白质组学进行分析，敲降28基因后发现酪蛋白、核糖体和蛋白酶体等相关基因及泛素蛋白表达水平显著上调，蛋白酶体活性下调；抑制蛋白酶体后酪蛋白相关基因发生同样变化趋势；进一步采用iTRAQ比较蛋白质组学分析发现，差异表达蛋白主要富集在核糖体、溶酶体、剪接体等通路中，表明敲降28基因可以通过泛素化信号通路影响BMECs中乳蛋白的合成，本研究结果也为奶牛产奶性状的分子育种研究提供了理论基础。