• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    宰后成熟期間能量物質(zhì)、pH值和肌原纖維小片化指數(shù)對(duì)秦川牛肉嫩度的影響及其機(jī)理

    2022-07-02 03:49:30李亞蕾羅瑞明張杏亞馬思麗
    食品科學(xué) 2022年11期
    關(guān)鍵詞:秦川牛嫩度肌原纖維

    羅 輝,楊 波,李亞蕾*,羅瑞明,張杏亞,馬思麗,姬 琛

    (寧夏大學(xué)食品與葡萄酒學(xué)院,寧夏 銀川 750021)

    秦川牛是陜甘寧地區(qū)優(yōu)勢(shì)畜種,近5 年存欄量在300萬(wàn) 頭左右,年均出欄量有100萬(wàn) 頭左右,年產(chǎn)肉量在50萬(wàn) t以上,具有一定地域特色。秦川牛肉質(zhì)細(xì)致、肉味濃郁、營(yíng)養(yǎng)全價(jià),肉用性能好。其背最長(zhǎng)肌具有瘦肉比例高、大理石花紋豐富等特點(diǎn),可用于制作高檔牛排,提升附加值,因此本研究選用秦川牛背最長(zhǎng)肌為研究對(duì)象。

    目前肉類(lèi)行業(yè)無(wú)論是在活體動(dòng)物中還是在屠體中缺乏一種快速、客觀的技術(shù)來(lái)預(yù)測(cè)評(píng)估牛肉貯藏過(guò)程中的品質(zhì)指標(biāo)變化。解決這一技術(shù)難題能提高牛肉行業(yè)的效率,并提高工業(yè)生產(chǎn)更高質(zhì)量肉制品的能力。為了實(shí)現(xiàn)這一目標(biāo),必須確定影響肉質(zhì)性狀變化的潛在機(jī)制,為此,研究肉質(zhì)的科研工作者普遍認(rèn)為,在試圖開(kāi)發(fā)預(yù)測(cè)肉品質(zhì)的技術(shù)之前,首先要揭示負(fù)責(zé)調(diào)控表觀變化的基因和蛋白質(zhì)網(wǎng)絡(luò)。這樣便于準(zhǔn)確地管理肉類(lèi)質(zhì)量性狀,并為生產(chǎn)者開(kāi)發(fā)基因選擇方案,以優(yōu)化種牛品質(zhì)。在過(guò)去的幾十年里,人們?yōu)閷?shí)現(xiàn)這一目標(biāo)已經(jīng)做出了許多嘗試,特別是利用蛋白質(zhì)組學(xué)技術(shù)來(lái)研究肉類(lèi)品質(zhì)變化,以更好地揭示牛肉嫩度的潛在機(jī)制。Sawdy等研究表明牛背最長(zhǎng)肌貯藏7 d后嫩度與肌球蛋白重鏈斷裂有關(guān),在肌肉嫩化過(guò)程中,肌球蛋白輕鏈1、2和重鏈1的表達(dá)量明顯降低;李升升利用非標(biāo)記定量(label-free quantitative,LFQ)蛋白質(zhì)組學(xué)技術(shù)研究牦牛平滑肌貯藏過(guò)程嫩度的變化機(jī)理,結(jié)果表明SLC25A4、Calpastatin、TPM3、UQCRC1、LMOD1、THBS4和ZYX 7種蛋白可能是牦牛平滑肌貯藏期嫩度變化的指示蛋白;Jia等也利用LFQ蛋白質(zhì)組學(xué)方法對(duì)牛胸最長(zhǎng)肌肌漿蛋白進(jìn)行研究發(fā)現(xiàn),在貯藏的前期peroxiredoxin-6表達(dá)量高的組隨著貯藏時(shí)間的延長(zhǎng)嫩度也會(huì)顯著增高,由此可以通過(guò)檢測(cè)活體組織或屠宰后不久牛胸最長(zhǎng)肌中peroxiredoxin-6含量來(lái)預(yù)測(cè)嫩度變化,因此peroxiredoxin-6可作為肉類(lèi)嫩度的潛在蛋白標(biāo)記物;Poleti等利用高清晰度質(zhì)譜儀(high definition mass spectrometry,HDMS)檢測(cè)Nellore牛背最長(zhǎng)肌蛋白表達(dá)量,發(fā)現(xiàn)118種蛋白質(zhì)的豐度發(fā)生了顯著變化,經(jīng)過(guò)生物學(xué)分析發(fā)現(xiàn)這些蛋白主要富集在丙酮酸代謝、分解代謝、肌動(dòng)蛋白結(jié)合、細(xì)胞骨架蛋白結(jié)合和鈣離子結(jié)合等途徑,由此推斷能量代謝和細(xì)胞骨架結(jié)構(gòu)的變化可能影響到嫩度。在幾十年的研究中,人們利用蛋白質(zhì)組學(xué)技術(shù)不斷揭示了調(diào)控嫩度變化的指示蛋白,調(diào)控嫩度的蛋白質(zhì)圖譜也在逐漸完善,嫩度變化預(yù)測(cè)技術(shù)也有了部分理論的支撐。但不同牛品種之間表現(xiàn)也會(huì)有較大差異,秦川牛背最長(zhǎng)肌嫩度變化機(jī)理報(bào)道鮮有發(fā)現(xiàn),本研究希望能進(jìn)一步完善這一理論基礎(chǔ)。

    本研究中利用4維-LFQ(4 dimension-LFQ,4D-LFQ)蛋白質(zhì)組學(xué)技術(shù)篩選秦川牛差異蛋白,4D-LFQ蛋白質(zhì)組學(xué)技術(shù)相對(duì)于3維蛋白質(zhì)組學(xué)在質(zhì)荷比(/)、離子強(qiáng)度和保留時(shí)間3個(gè)維度添加了離子淌度(mobility)。離子淌度為牛肉肌肉組織這種復(fù)雜樣品體系的蛋白組檢測(cè)提供了更多可能,可根據(jù)離子的截面和形狀進(jìn)一步對(duì)蛋白進(jìn)行分離鑒定,使蛋白組檢測(cè)數(shù)據(jù)更為準(zhǔn)確、全面。

    1 材料與方法

    1.1 材料與試劑

    25 月齡秦川公牛 甘肅省平?jīng)鍪袥艽h旭康食品有限責(zé)任公司。

    蛋白酶抑制劑 德國(guó)Calbiochem 公司;胰酶美國(guó)Promega公司;乙腈 美國(guó)Fisher Chemical公司;甲酸瑞士Fluka公司;三氟乙酸、尿素、碘代乙酰胺、二硫蘇糖醇、三乙基碳酸氫銨、三氟乙酸 美國(guó)Sigma公司;BCA試劑盒 上海碧云天生物技術(shù)有限公司;pH校準(zhǔn)液上海儀電科學(xué)儀器股份有限公司。

    1.2 儀器與設(shè)備

    超聲波細(xì)胞破碎機(jī) 那艾精密儀器(上海)有限公司;205便捷式pH計(jì) 德國(guó)德圖集團(tuán);MiniVac Alpha冷凍離心機(jī) 美國(guó)Scan speed公司;timsTOF Pro質(zhì)譜儀、NanoElute超高效液相色譜儀 德國(guó)Bruker公司;Forma 994超低溫冰箱 美國(guó)Thermo公司;TA-XT plus質(zhì)構(gòu)儀英國(guó)Stable Micro Systems公司。

    1.3 方法

    1.3.1 樣品采集

    選用3 頭生長(zhǎng)發(fā)育正常、體壯無(wú)病、25 月齡左右的去勢(shì)秦川公牛,屠宰后取牛兩側(cè)背最長(zhǎng)肌并分割為質(zhì)量約(100f10)g的肉樣,隨機(jī)分成5 組并置于透氣性為23.5 g/(m?d)的聚乙烯保鮮袋中。立即取其中一組樣品檢測(cè)各項(xiàng)指標(biāo)(視為宰后第0天樣品),將蛋白質(zhì)組學(xué)、能量物質(zhì)檢測(cè)和剪切力檢測(cè)實(shí)驗(yàn)所用樣品置于液氮中保存。其他4 組貯藏于0~4 ℃的冰箱中,分別在第2、4、6、8天時(shí)取出,測(cè)定相應(yīng)指標(biāo)。

    1.3.2 能量物質(zhì)含量的測(cè)定

    稱(chēng)取50 mg樣品置于離心管中,加入4 ℃預(yù)冷的1 mL 混合液(含0.6 mol/L高氯酸和1 mmol/L乙二胺四乙酸),勻漿2 min后4 ℃、10 000h條件下離心10 min,取700 μL上清液加入0.2 mL 1 mol/L NaOH。搖勻后反應(yīng)5 min,取0.5 mL混合溶液用聚碳酸酯濾膜過(guò)濾后于-80 ℃冰箱保存?zhèn)溆?。采用高相液相色譜法檢測(cè)ATP、ADP、AMP、NADH的含量。流動(dòng)相為220 mmol/L硫酸鉀緩沖液(含10%甲醇并使用四丁基氫氧化銨調(diào)節(jié)至pH 6.5),等度洗脫,流速為0.8 mL/min,檢測(cè)波長(zhǎng)為254 nm。

    1.3.3 剪切力測(cè)定

    參考文獻(xiàn)[16]測(cè)定樣品剪切力,以表征樣品的嫩度。取出貯藏在0~4 ℃環(huán)境下的肉樣,插入熱電偶溫度計(jì)至肉樣中央并置于蒸煮袋中,在85 ℃水浴鍋中加熱,待肉樣中心溫度達(dá)80 ℃時(shí)取出,流動(dòng)自來(lái)水沖洗蒸煮袋將肉樣冷卻至室溫,切除表面,避開(kāi)筋腱等部分,平行于肌纖維方向用口徑為1.27 cm的取樣器切割制作肉柱,用質(zhì)構(gòu)儀“V”型剪切刀檢測(cè)肉柱,剪切刀垂直于肌纖維方向,剪切速率與剪切距離分別為1.5 mm/s、40 mm,每個(gè)肉樣檢測(cè)3 次,結(jié)果取平均值。

    1.3.4 pH值測(cè)定

    參照文獻(xiàn)[17]的方法,pH計(jì)使用前用pH校準(zhǔn)液校正,取4 cmh4 cmh4 cm樣品(60 g左右)背最長(zhǎng)肌,探頭沿肌纖維方向插入肉樣內(nèi)部,每個(gè)樣品取5個(gè)點(diǎn)測(cè)定,結(jié)果取平均值。

    1.3.5 4D-LFQ蛋白質(zhì)組學(xué)技術(shù)檢測(cè)蛋白質(zhì)量濃度

    參照文獻(xiàn)[18]的方法提取蛋白質(zhì)。4 ℃預(yù)冷后的研缽放入適量樣品和液氮研磨成粉末狀,加樣品4 倍體積的裂解緩沖液(等體積的8 mol/L尿素和1%蛋白質(zhì)酶抑制劑)。超聲裂解后移入離心管,4 ℃、12 000h條件下離心10 min,取上清液利用BCA試劑盒測(cè)定蛋白質(zhì)量濃度。

    參照文獻(xiàn)[19]的方法對(duì)提取的蛋白進(jìn)行胰酶酶解。取上述蛋白溶液加入二硫蘇糖醇溶液,調(diào)節(jié)蛋白質(zhì)濃度為5 mmol/L。56 ℃避光環(huán)境下還原30 min再加入100 mmol/L碘代乙酰胺溶液,將蛋白質(zhì)濃度調(diào)節(jié)至1 mmol/L,25 ℃避光環(huán)境下靜止12 min。用四乙基溴化銨溶液將樣品中的尿素稀釋到2 mol/L以下,按照蛋白質(zhì)量的2%添加胰酶,37 ℃酶解12 h,再按照蛋白質(zhì)量的1%添加胰酶,37 ℃酶解4 h。

    液相色譜-質(zhì)譜聯(lián)用分析:使用NanoElute超高效液相系統(tǒng)進(jìn)行肽段分離,流動(dòng)相A為0.1%(體積分?jǐn)?shù))甲酸水溶液,用來(lái)溶解蛋白片段;流動(dòng)相B為0.1%(體積分?jǐn)?shù))甲酸-乙腈溶液;流速為300 nL/min。梯度洗脫程序?yàn)椋?~3 min,4%~6% B;3~70 min,6%~25% B;70~82 min,25%~32% B;82~87 min,32%~80% B;87~90 min,80%~100% B。肽段分離后注入到capillary離子源中進(jìn)行電離,電壓設(shè)置為1.4 kV,使用飛行時(shí)間質(zhì)譜對(duì)肽段母離子及其二級(jí)碎片進(jìn)行檢測(cè)和分析,設(shè)置掃描范圍100~1 700/。通過(guò)平行累積串行碎裂模式(parallel accumulation serial fragmentation,PASEF)進(jìn)行數(shù)據(jù)采集。將采集后一級(jí)質(zhì)譜用10 次PASEF模式采集母離子電荷數(shù)在0~5范圍獲得二級(jí)譜圖,為了避免母離子的重復(fù)掃描串聯(lián)質(zhì)譜掃描的動(dòng)態(tài)排除時(shí)間設(shè)置為24 s。

    1.3.6 數(shù)據(jù)庫(kù)檢索

    二級(jí)質(zhì)譜的數(shù)據(jù)檢索采用Maxquant(v1.6.6.0)軟件進(jìn)行。各項(xiàng)參數(shù)設(shè)置如下:選擇Bos_taurus_9913_PR_20190915.fasta(37 948 條序列)為數(shù)據(jù)庫(kù),加入常見(jiàn)的污染庫(kù)排除污染蛋白質(zhì)的影響,加入反庫(kù)特征信息匹配排除計(jì)算隨機(jī)匹配造成的假陽(yáng)性率(false discovery rate,F(xiàn)DR)影響;酶切方式為T(mén)rypsin/P,Main search與First search的二級(jí)碎片離子的質(zhì)量誤差容忍度為0.04 Da,一級(jí)母離子質(zhì)量誤差容忍度設(shè)為40h10Da;肽段最大修飾數(shù)、漏切位點(diǎn)數(shù)和肽段最小長(zhǎng)度分別設(shè)置為5、2、7。固定修飾為半胱氨酸烷基化,可變修飾為蛋白N端的乙?;图琢虬彼岬难趸k钠ヅ鋱D譜鑒定的FDR和蛋白鑒定均設(shè)置為1%。

    1.3.7 生物信息學(xué)分析

    篩選出的差異蛋白質(zhì)利用軟件ClusterProfiler(3.4.4)進(jìn)行生物信息學(xué)分析,比對(duì)后獲得差異基因的基因本體(gene ontology,GO)(包括生物學(xué)過(guò)程(biological process,BP)、細(xì)胞組分(cellular component,CC))功能信息。利用String 10.0數(shù)據(jù)庫(kù)(http://string-db.org/)和Cytoscape 3.6.1軟件繪制差異蛋白相互作用網(wǎng)絡(luò)圖。

    1.3.8 肌原纖維小片化指數(shù)的測(cè)定

    稱(chēng)取肉樣1 g使用冷凍研磨儀攪碎后置于離心管中,加入10 mL肌原纖維小片化指數(shù)(myofibril fragmentation index,MFI)緩沖液(含100 mmol/L氯化鉀、11.2 mmol/L磷酸氫二鉀、8.8 mmol/L磷酸二氫鉀、1 mmol/L乙二醇雙四乙酸、1 mmol/L氯化鎂、1 mmol/L疊氮化鈉),4 ℃、1 200 r/min勻漿3 min,然后4 ℃、1 000h離心15 min收集上清液;加入10 mL MFI緩沖液重懸沉淀并按上述步驟再次勻漿后離心,收集上清液;再加入5 mL預(yù)冷MFI緩沖溶液使沉淀充分懸浮后用200 目篩過(guò)濾,再用5 mL MFI緩沖液清洗離心管,然后用200 目篩過(guò)濾,合并上清液和濾液,并采用考馬斯亮藍(lán)法測(cè)定混合溶液的蛋白質(zhì)量濃度。用MFI緩沖液調(diào)整蛋白質(zhì)量濃度為0.5 mg/mL左右,在595 nm波長(zhǎng)處測(cè)定吸光度,每個(gè)樣品重復(fù)3 次,取平均值帶入下式計(jì)算MFI。

    1.4 數(shù)據(jù)處理及分析

    所有數(shù)據(jù)均以平均值±標(biāo)準(zhǔn)差表示,采用SPSS 25軟件進(jìn)行單因素分析,采用鄧肯檢驗(yàn)進(jìn)行顯著性分析,<0.05表示差異顯著,<0.01表示差異極顯著。采用Origin 2018軟件作圖,采用微生信在線平臺(tái)(http://www.bioinformatics.com.cn/)進(jìn)行秦川牛肉背最長(zhǎng)肌差異蛋白數(shù)據(jù)分析并作圖。采用SPSS 25軟件分別對(duì)差異蛋白表達(dá)量與剪切力、能量物質(zhì)含量進(jìn)行Pearson相關(guān)性分析,篩選出與嫩度相關(guān)的蛋白質(zhì)。

    2 結(jié)果與分析

    2.1 秦川牛背最長(zhǎng)肌貯藏過(guò)程中ATP、ADP、AMP、NADH含量的變化

    采用高效液相色譜法測(cè)定了秦川牛背最長(zhǎng)肌不同貯藏期能量物質(zhì)含量。如表1所示,隨著貯藏時(shí)間的延長(zhǎng),ATP的含量呈顯著下降趨勢(shì)(<0.05),由0 d時(shí)的(10.15f0.60)μmol/g下降至8 d時(shí)的(0.50f0.20)μmol/g,其中在0~2 d下降速率較快,在2~8 d下降速率較為緩慢,第8天時(shí)基本檢測(cè)不出ATP,ADP、AMP的含量在整個(gè)貯藏期(0~8 d)均呈現(xiàn)下降趨勢(shì)。秦川牛宰后肌肉細(xì)胞中ATP進(jìn)行著不可逆的降解過(guò)程,由于ATP含量的減少,ADP、AMP含量隨之減少。貯藏過(guò)程中NADH含量變化可反映組織細(xì)胞中氧化還原和能量變化狀態(tài)。實(shí)驗(yàn)結(jié)果中NADH含量逐漸降低,說(shuō)明宰后肌肉細(xì)胞發(fā)生變化,使細(xì)胞呼吸功能失調(diào),氧化磷酸化進(jìn)程受阻,引起NADH含量出現(xiàn)顯著下降的現(xiàn)象。

    表1 秦川牛背最長(zhǎng)肌不同貯藏時(shí)間下ATP、ADP、AMP、NADH含量變化Table 1 Changes in ATP, ADP, AMP and NADH contents in longissimus dorsi muscle of Qinchuan cattle with storage time

    2.2 秦川牛背最長(zhǎng)肌貯藏過(guò)程剪切力的變化

    嫩度是影響秦川牛肉口感的重要品質(zhì)屬性,也是影響消費(fèi)者購(gòu)買(mǎi)欲的重要指標(biāo)。咀嚼破碎肉制品時(shí)口腔所受的抵抗力大小可稱(chēng)為肉制品的嫩度,實(shí)驗(yàn)結(jié)果以剪切力表征嫩度。如圖1所示,秦川牛肉在貯藏過(guò)程中剪切力總體呈先上升后下降的趨勢(shì)(<0.01),其下降幅度大于上升幅度。第0天時(shí)剪切力為(120.25f5.15)N,在第4天時(shí)達(dá)到最大值,為(157.94f2.53)N,第8天時(shí)又降至(93.41f5.08)N,通過(guò)以上剪切力的變化可以看出秦川牛背最長(zhǎng)肌在貯藏中的嫩度變化情況,在0~4 d時(shí)逐漸降低,第4天時(shí)到達(dá)最低點(diǎn),之后逐漸升高,直至第8天。

    圖1 秦川牛背最長(zhǎng)肌貯藏期間剪切力變化Fig.1 Changes in shear force of longissimus dorsi muscle with storage time

    2.3 秦川牛背最長(zhǎng)肌貯藏過(guò)程中pH值變化

    秦川牛背最長(zhǎng)肌在貯藏過(guò)程中pH值變化情況如圖2所示,宰后0 d時(shí)pH值為6.78f0.02,在0~2 d貯藏期內(nèi)下降速率較快,2~4 d貯藏期下降速度放緩后到達(dá)最低點(diǎn)5.37f0.03,4~8 d時(shí)呈緩慢上升趨勢(shì),8 d時(shí)上升至5.50f0.04。整體下降幅度大于上升幅度。無(wú)氧條件下,機(jī)體通過(guò)無(wú)氧糖酵解產(chǎn)能的同時(shí)產(chǎn)生乳酸;此外,宰后ATP水解產(chǎn)生磷酸,這些酸性物質(zhì)長(zhǎng)期積累并無(wú)法分解與轉(zhuǎn)運(yùn),使得肌肉組織的pH值迅速下降。

    圖2 秦川牛背最長(zhǎng)肌貯藏期間pH值的變化Fig.2 Changes in pH in longissimus dorsi muscle of Qinchuan cattle with storage time

    2.4 蛋白質(zhì)組學(xué)數(shù)據(jù)分析及關(guān)鍵蛋白質(zhì)篩選結(jié)果

    實(shí)驗(yàn)中采用4D-LFQ對(duì)貯藏0、4、8 d的樣品蛋白質(zhì)豐富度變化進(jìn)行定量分析,利用Maxquant(v1.6.6.0)軟件對(duì)所得的蛋白質(zhì)二級(jí)結(jié)構(gòu)進(jìn)行質(zhì)譜檢索,共檢索確定蛋白1 149個(gè),設(shè)置1.200 倍為差異顯著倍數(shù),共篩選出120個(gè)差異顯著的蛋白質(zhì),為進(jìn)一步確定數(shù)據(jù)的可信度,對(duì)差異蛋白進(jìn)行了聚類(lèi)分析與主成分分析。

    秦川牛背最長(zhǎng)肌差異蛋白豐富度層次聚類(lèi)分析如圖3A所示,各組的3個(gè)平行樣品均可聚為一類(lèi),說(shuō)明樣品重復(fù)性較好,且不同貯藏時(shí)間樣品差異明顯;同時(shí),4、8 d的樣品距離更近,與0 d樣品相區(qū)分。差異蛋白的主成分分析結(jié)果如圖3B所示,PC1、PC2的累計(jì)貢獻(xiàn)率達(dá)到72.23%,因此這兩個(gè)主成分可較好地代表差異蛋白表達(dá)的信息。PC1能將0 d與其他樣品區(qū)分開(kāi),PC2能將4 d與8 d樣品區(qū)分開(kāi)。綜上可判斷,秦川牛背最長(zhǎng)肌差異蛋白0、4、8 d樣品重復(fù)性較好且有顯著差異,具有一定代表性和可信度。

    圖3 秦川牛肉背最長(zhǎng)肌差異蛋白整體數(shù)據(jù)分析結(jié)果Fig.3 Analysis of differential proteins between longissimus dorsi muscle of Qinchuan cattle stored for different periods

    采用SPSS 25軟件對(duì)差異蛋白表達(dá)量與剪切力進(jìn)行Pearson相關(guān)性分析,篩選出與嫩度變化相關(guān)的11種差異蛋白,如表2所示。

    表2 秦川牛背最長(zhǎng)肌與嫩度相關(guān)的差異蛋白Table 2 Differential proteins related to beef tenderness in longissimus dorsi muscle of Qinchuan cattle

    2.5 差異蛋白GO注釋結(jié)果

    對(duì)秦川牛肉背最長(zhǎng)肌篩選出的差異蛋白進(jìn)行GO分析后,所涉及到的生物過(guò)程與細(xì)胞組分前10的條目情況如圖4所示,生物過(guò)程主要集中在氧化還原過(guò)程(ATP5F1D、NDUFA6、NDUFB5、SUCLA2)、線粒體組織(ATP5F1D、NDUFA6、NDUFB5)、NADH脫氫酶符合體組裝(NDUFA6、NDUFB5)、細(xì)胞呼吸(ATP5F1D、SUCLA2)等過(guò)程;細(xì)胞組分主要集中在催化復(fù)合物(NDUFA6、PPP3R1、NDUFB5、HNRNPK、PSMD13)、線粒體蛋白質(zhì)復(fù)合物(ATP5F1D、NDUFA6、NDUFB5、ATP5F1C)、細(xì)胞器內(nèi)膜(ATP5F1D、NDUFA6、NDUFB5、ATP5F1C)、肌纖維膜(PPP3R1、CAMK2D)等。

    圖4 秦川牛肉背最長(zhǎng)肌差異蛋白GO分析結(jié)果Fig.4 Gene ontology (GO) analysis of differential proteins in longissimus dorsi muscle of Qinchuan cattle

    2.6 差異蛋白相互作用網(wǎng)絡(luò)圖

    圖5是利用Cytoscape 3.6.1軟件和String 10.0數(shù)據(jù)庫(kù)制作的差異蛋白相互作用網(wǎng)絡(luò)圖。此網(wǎng)絡(luò)圖聚類(lèi)系數(shù)為0.691,蛋白質(zhì)相互作用富集值為1.39h10。如圖5所示,ATP5D、ATP5C1、SUCLA2、SUCLG1、NDUFB5、NDUFA6的相互作用較為強(qiáng)烈,通過(guò)GO分析發(fā)現(xiàn)在氧化還原過(guò)程,線粒體組織條目有較密切互作關(guān)系,分析發(fā)現(xiàn)以上6種蛋白質(zhì)組成的集團(tuán)均參與調(diào)解宰后能量代謝,最終通過(guò)ATP含量和pH值變化影響嫩度變化;CAMK2D與PPP3R1的相關(guān)作用較為強(qiáng)烈,分析發(fā)現(xiàn)兩者具有協(xié)同關(guān)系,共同參與細(xì)胞凋亡,促進(jìn)MFI的增大,肌原纖維受到破壞,最終使得嫩度降低;CTSD、HNRNPK、PSMD13暫未與其他蛋白發(fā)現(xiàn)相互作用。

    圖5 差異蛋白相互作用網(wǎng)絡(luò)圖Fig.5 Interaction network between differential proteins

    2.7 pH值、能量物質(zhì)含量及相關(guān)差異蛋白對(duì)嫩度的影響

    利用SPSS 25軟件對(duì)表2中11種蛋白質(zhì)的表達(dá)量與能量物質(zhì)含量進(jìn)行Pearson相關(guān)性分析,篩選出6種與能量變化相關(guān)的差異蛋白質(zhì),結(jié)果如表3所示,它們分別為ATP5F1D、ATP5F1C、NDUFA6、NDUFA6、SUCLA2、SUCG1。

    表3 差異蛋白質(zhì)表達(dá)量與能量物質(zhì)含量的相關(guān)性Table 3 Correlation between expression of differential proteins and energy substance contents

    表2中ATP5F1D是ATP合酶的亞基δ,ATP5F1C是ATP合酶亞基γ,這兩種蛋白質(zhì)參與ATP合成酶的組裝,是ATP合酶5個(gè)亞基中的2個(gè),對(duì)ATP合酶催化活性有重要作用,承擔(dān)著離子傳遞的作用。ATP含量合酶參與ATP的合成與消耗,所以ATP合酶直接影響機(jī)體細(xì)胞的供能和耗能,其表達(dá)量也會(huì)直接反映細(xì)胞中能量需求情況和三羧酸循環(huán)產(chǎn)能的情況。由表3可知,ATP5F1D、ATP5F1C表達(dá)量與能量基礎(chǔ)物質(zhì)ATP含量呈極顯著負(fù)相關(guān)(<0.01)。表2中NDUFA6為NADH脫氫酶的附件亞基。NADH脫氫酶參與線粒體內(nèi)膜上的呼吸作用,對(duì)電子從NADH傳遞給輔酶Q有催化作用,NDUFA6表達(dá)量與能量基礎(chǔ)物質(zhì)ATP含量呈顯著負(fù)相關(guān)(<0.05),與NADH含量呈極顯著正相關(guān)(<0.01)。ATP5F1D和NDUFA6在0~4 d的貯藏過(guò)程中均有高水平的表達(dá),ATP5F1D上調(diào)1.819 倍,NDUFA6上調(diào)1.467 倍。兩者均是糖酵解過(guò)程中關(guān)鍵酶的組件,其高水平的表達(dá)說(shuō)明機(jī)體細(xì)胞在切斷氧氣供應(yīng)時(shí),能提高參與產(chǎn)能的酶含量,利用無(wú)氧糖酵解過(guò)程產(chǎn)能,在此過(guò)程中大量產(chǎn)生乳酸,使得肌肉pH值下降。

    琥珀酰輔酶A合成酶在檸檬酸循環(huán)中發(fā)揮作用,在琥珀酸硫激酶的催化下,琥珀酰輔酶A高能鍵水解生成琥珀酸,同時(shí)使GDP磷酸化生成GTP,GTP再與ADP作用生成ATP。表2中SUCLA2作為琥珀酸-CoA連接酶[ADP形成]亞基β,提供酶的核苷酸特異性,并與底物琥珀酸結(jié)合,而輔酶A和磷酸的結(jié)合位點(diǎn)在α亞基。而篩選出SUCG1就是琥珀酸-CoA連接酶[ADP形成]亞基α。通過(guò)表3也可以看出,SUCLA2的表達(dá)量與ATP、NADH含量呈極顯著正相關(guān)(<0.01),與ADP、AMP顯著正相關(guān)(<0.05);SUCG1表達(dá)量與ATP含量呈含量顯著正相關(guān)(<0.05)。表2中NDUFB5與NDUFA6同為NADH脫氫酶的附件亞基,在磷酸氧化過(guò)程中有著同等作用,從表3中也可以看出NDUFB5的表達(dá)量與ATP、NADH含量均呈顯著負(fù)相關(guān)(<0.05)。從表2可以看出,SUCLA2的表達(dá)量在0~8 d調(diào)整0.762 倍,NDUFB5的表達(dá)量在4~8 d調(diào)整0.497 倍,SUCG1的表達(dá)量在4~8 d調(diào)整0.742 倍。SUCLA2、NDUFB5、NDUFA6是在三羧酸循環(huán)、呼吸電子鏈傳遞、通過(guò)化學(xué)滲透耦合進(jìn)行ATP合成、解偶聯(lián)蛋白質(zhì)中起關(guān)鍵作用的蛋白,在4~8 d的貯藏期內(nèi)均有下調(diào)趨勢(shì),說(shuō)明無(wú)氧糖酵解過(guò)程可能會(huì)受到抑制,能量和乳酸的生成量相應(yīng)減少。

    秦川牛屠宰后血液流動(dòng)停止,肌肉細(xì)胞通過(guò)調(diào)節(jié)以上蛋白表達(dá)量,啟動(dòng)糖酵解過(guò)程產(chǎn)能,本實(shí)驗(yàn)數(shù)據(jù)顯示,0~4 d時(shí)秦川牛肉的ATP的含量減少至(1.20f0.45)μmol/g,為初始狀態(tài)的11.82%,在能量物質(zhì)ATP下降到宰前初始含量的20%以下時(shí),肌纖維發(fā)生交聯(lián),肌肉會(huì)出現(xiàn)僵直現(xiàn)象,由圖1可知,在0~4 d貯藏期秦川牛肉剪切力也隨之上升,4 d后能量物質(zhì)含量下降速率有所減慢,對(duì)肌肉收縮的影響可能會(huì)減小。無(wú)氧糖酵解產(chǎn)生能量的同時(shí)產(chǎn)生乳酸;此外宰后ATP水解產(chǎn)生磷酸,這些酸性物質(zhì),使得肌肉組織的pH值迅速下降。4 d時(shí)pH值為5.37f0.03,低pH值改變了質(zhì)子的穩(wěn)態(tài),降低了肌原纖維蛋白之間的靜電斥力,導(dǎo)致肌肉纖維不可逆?zhèn)认蚴湛s(肌節(jié)延伸),使得肉嫩度變差;低pH值使纖維蛋白暴露出更多水解位點(diǎn),有助于水解酶的水解作用;低pH值會(huì)降低溶酶體膜的穩(wěn)定性,使得肌肉細(xì)胞釋放內(nèi)源水解酶,同時(shí)pH值的降低會(huì)使肌肉細(xì)胞處于酸性環(huán)境,打破細(xì)胞中的離子平衡,促進(jìn)鈣離子的釋放,Ca濃度的提高會(huì)激活鈣蛋白酶,降解纖維蛋白,促使嫩度升高,但pH值的降低也會(huì)使得-鈣激活酶發(fā)生自溶現(xiàn)象,降低酶活性,降低肌原纖維蛋白的水解速率;也有研究表明pH值小于5.8時(shí),組織酶D對(duì)水解鴕鳥(niǎo)纖維有促進(jìn)作用,由此可判斷pH值會(huì)影響到內(nèi)源酶的催化活性,對(duì)肌原纖維蛋白的降解速率產(chǎn)生一定影響。通過(guò)圖6也可以看出,剪切力和ATP含量、pH值均呈顯著負(fù)相關(guān)(<0.05)。綜上可推測(cè),0~4 d能量物質(zhì)含量的減少和pH值的降低會(huì)使得肌肉的僵直,肌肉嫩度降低,4~8 d能量物質(zhì)含量變化較小,對(duì)肌肉的影響也會(huì)相應(yīng)降低,長(zhǎng)期低pH值下的各種途徑綜合作用對(duì)嫩度提高的速率和嫩度提高的程度均有影響。

    圖6 秦川牛背最長(zhǎng)肌能量物質(zhì)含量、pH值和剪切力之間的相關(guān)系數(shù)Fig.6 Correlation coefficients between energy substance contents, pH and shear force

    2.8 與MFI相關(guān)的差異蛋白對(duì)宰后秦川牛肉嫩度變化的影響

    秦川牛肉在貯藏過(guò)程中MFI變化情況如圖7所示,在整個(gè)貯藏過(guò)程中MFI呈顯著上升趨勢(shì)(<0.05),由0 d時(shí)28.50f5.1持續(xù)增長(zhǎng)至8 d時(shí)的99.98f1.20,變化了250.81%。此結(jié)果與李升升對(duì)牦牛平滑肌嫩度形成機(jī)理研究中的MFI增加趨勢(shì)一致。

    圖7 秦川牛背最長(zhǎng)肌不同貯藏時(shí)間MFI變化情況Fig.7 Changes in MFI of longissimus muscle of Qinchuan cattle with storage time

    利用SPSS 25軟件將表2中11種與嫩度相關(guān)的蛋白質(zhì)表達(dá)量與與MFI進(jìn)行相關(guān)性分析,篩選出5種(CTSD、PSMD13、PPP3R1、CAMK2D、HNRNPK)與MFI有關(guān)的關(guān)鍵蛋白質(zhì),相關(guān)性分析結(jié)果如表4所示。

    表4 秦川牛背最長(zhǎng)肌貯藏期間剪切力、MFI與差異蛋白含量相關(guān)系數(shù)Table 4 Correlation coefficients between shear force, MFI and differential protein contents

    CTSD為組織蛋白酶D,位于溶酶體,pH值降低等因素會(huì)影響溶酶體膜穩(wěn)定性,使得CTSD被釋放,后經(jīng)過(guò)ADAM30D切割后可激活活性,在蛋白質(zhì)的降解與細(xì)胞凋亡中起作用,CTSD在0~8 d中的表達(dá)量上調(diào)1.303 倍(表2),CTSD表達(dá)量與剪切力呈顯著負(fù)相關(guān)(<0.05),與MFI呈顯著正相關(guān)(<0.05),CTSD的高表達(dá)可能會(huì)出現(xiàn)肌肉細(xì)胞更高水平的自噬。PSMD13是26S蛋白酶體非ATP酶調(diào)節(jié)亞基13,26S蛋白酶體非ATP酶可以去除可能損害細(xì)胞功能的錯(cuò)誤折疊或受損的蛋白質(zhì),還可以通過(guò)去除不再需要的功能蛋白質(zhì),維持蛋白質(zhì)動(dòng)態(tài)平衡,同時(shí)參與細(xì)胞凋亡的進(jìn)程,PDMD13在0~8 d中上調(diào)1.336 倍(表2),PSMD13表達(dá)量與剪切力呈極顯著負(fù)相關(guān)(<0.01),與MFI呈顯著正相關(guān)(<0.05),PSMD13表達(dá)量發(fā)生變化會(huì)影響到細(xì)胞凋亡和受損蛋白質(zhì)去除。HNRNPK是與mRNA結(jié)合的主要蛋白質(zhì),可能在核內(nèi)不均一RNA(heterogeneous nuclear RNA,hnRNA)的核代謝中通過(guò)轉(zhuǎn)錄激活和抑制發(fā)揮作用,在糖基化中起促進(jìn)作用,作為抑癌基因/的轉(zhuǎn)錄輔激活因子,在/和/誘導(dǎo)中起作用;在轉(zhuǎn)錄抑制方面,它通過(guò)與/誘導(dǎo)的非編碼RNA(-)相互作用發(fā)揮作用,這種相互作用是誘導(dǎo)細(xì)胞凋亡所發(fā)生的,HNRNPK在0~8 d調(diào)整0.722 倍(表2),HNRNPK表達(dá)量與剪切力呈顯著正相關(guān)(<0.05),與MFI呈顯著負(fù)相關(guān)(<0.05),HNRNPK在貯藏期內(nèi)的異常表達(dá),會(huì)直接影響到肌細(xì)胞的凋亡進(jìn)程。

    研究中發(fā)現(xiàn)了兩種關(guān)鍵蛋白質(zhì)在細(xì)胞凋亡進(jìn)程中具有協(xié)同作用。PPP3R1是鈣調(diào)磷酸酶B亞基1型,是一種鈣調(diào)節(jié)蛋白刺激的蛋白磷酸酶,提供細(xì)胞對(duì)Ca的敏感性,在4~8 d上調(diào)1.309 倍,提高了細(xì)胞對(duì)鈣離子的敏感性。CAMK2D是鈣/鈣調(diào)素依賴(lài)性蛋白激酶II型亞基δ,CAMK2D通過(guò)靶向鈣離子通道、轉(zhuǎn)運(yùn)蛋白參與細(xì)胞鈣離子內(nèi)流、肌漿網(wǎng)鈣釋放、肌漿網(wǎng)鈣攝取,影響細(xì)胞的鈣穩(wěn)態(tài)和興奮收縮偶聯(lián)的調(diào)節(jié),CAMK2D表達(dá)量在4~8 d上調(diào)1.506 倍。由表4可以看出,PPP3R1表達(dá)量與剪切力呈極顯著負(fù)相關(guān)(<0.01),與MFI呈極顯著正相關(guān)(<0.01);CAMK2D表達(dá)量與剪切力呈顯著負(fù)相關(guān)(<0.05),與MFI呈顯著正相關(guān)(<0.05)。初步判斷PPP3R1和CAMK2D在細(xì)胞凋亡進(jìn)程中具有協(xié)同作用,CAMK2D促進(jìn)細(xì)胞對(duì)Ca的吸收,提高細(xì)胞中Ca的濃度;PPP3R1可提高細(xì)胞對(duì)鈣離子的敏感性,Ca超負(fù)荷的積累是觸發(fā)細(xì)胞凋亡的因素,這兩種蛋白同時(shí)在4~8 d時(shí)高表達(dá),共同作用促使細(xì)胞凋亡進(jìn)程的發(fā)生,降低肌原纖維結(jié)構(gòu)強(qiáng)度,使得嫩度提高。

    MFI可直接反映肌纖維的降解程度,增大原因是由于Z線蛋白的降解和肌原纖維細(xì)胞凋亡。宰后無(wú)氧糖酵解導(dǎo)致pH值降低,使得細(xì)胞離子失衡,鈣離子濃度升高;降低溶酶體膜與細(xì)胞膜的穩(wěn)定性,釋放并激活內(nèi)源酶,CTSD、PSMD13、PPP3R1、CAMK2D、HNRNPK 5種參與細(xì)胞凋亡、蛋白質(zhì)水解的蛋白豐富度發(fā)生了改變,加快了肌原纖維蛋白的水解與細(xì)胞凋亡的發(fā)生,促使MFI升高了250.81%,降低了肌原纖維骨架強(qiáng)度,使得秦川牛背最長(zhǎng)肌嫩度升高。

    2.9 能量物質(zhì)、pH值、MFI對(duì)秦川牛肉背最長(zhǎng)肌嫩度變化不同時(shí)間段產(chǎn)生影響的分析

    秦川牛宰后能量和pH值變化較為迅速,對(duì)肌纖維收縮的影響較早,0~4 d時(shí)肌肉僵直對(duì)嫩度變化影響較大,4~8 d時(shí)隨著能量物質(zhì)趨于穩(wěn)定,對(duì)肌肉收縮的影響減小,對(duì)肌肉嫩度的影響也會(huì)隨之降低。低pH值與5種參與細(xì)胞凋亡有關(guān)的內(nèi)源酶相互作用使得MFI升高,是長(zhǎng)期積累的過(guò)程,當(dāng)MFI升高到一定程度,肌原纖維骨架強(qiáng)度顯著降低時(shí),才會(huì)使得肌肉嫩度升高,由圖1剪切力數(shù)據(jù)可以看出MFI作用主要集中在4~8 d時(shí),所以MFI作用具有滯后性。

    3 結(jié) 論

    秦川牛肉在貯藏過(guò)程中6種蛋白(ATP5F1D、ATP5F1C、NDUFB5、NDUFA6、SUCLA2、SUCG1)通過(guò)呼吸電子鏈傳遞、通過(guò)化學(xué)滲透耦合進(jìn)行ATP合成、解偶聯(lián)蛋白質(zhì)等途徑影響宰后糖酵解過(guò)程,ATP含量降至初始狀態(tài)(0 d)的11.82%,引起肌纖維發(fā)生交聯(lián),肌肉會(huì)出現(xiàn)僵直,造成肌肉剪切力增大,嫩度下降。能量物質(zhì)含量宰后變化較快且后期穩(wěn)定,所以對(duì)肌肉嫩度影響主要集中在0~4 d時(shí)。宰后酸性物質(zhì)的積累造成肌肉pH值下降,導(dǎo)致肌原纖維蛋白之間的靜電斥力減小,影響了肌肉收縮;參與細(xì)胞鈣離子濃度、促使肌原纖維受體暴露、調(diào)節(jié)內(nèi)源酶活性等過(guò)程使得嫩度升高。5種蛋白(PPP3R1、CAMK2D、HNRNPK、PSMD13、CTSD)通過(guò)肌纖維膜、膜蛋白復(fù)合物、細(xì)胞器內(nèi)膜、鈣信號(hào)通路調(diào)控肌纖維細(xì)胞的凋亡和蛋白的水解,造成肌原纖維結(jié)構(gòu)被破壞,使得MFI逐漸升高,肌原纖維骨架強(qiáng)度降低,促使肌肉嫩度升高。長(zhǎng)期低pH值和5種與細(xì)胞凋亡和蛋白水解有關(guān)的內(nèi)源酶相互作用下的MFI升高是逐漸積累的過(guò)程,作用時(shí)間滯后,從剪切力數(shù)據(jù)可以看出,MFI途徑對(duì)嫩度降低的影響主要集中在4 d后。

    猜你喜歡
    秦川牛嫩度肌原纖維
    秦川牛遺傳資源保護(hù)和開(kāi)發(fā)利用工作的思考
    秦川牛改良效果研究和建議
    飼糧中添加不同抗氧化劑對(duì)熟化期羊肉嫩度的影響及動(dòng)力學(xué)分析
    不同嫩度羊肉中鈣蛋白酶的差異
    不同嫩度茶鮮葉原料稀土元素含量的變化
    秦川牛種質(zhì)資源保護(hù)與肉用選育改良
    肌原纖維蛋白與大豆分離蛋白復(fù)合體系乳化性的研究
    TG酶協(xié)同超高壓處理對(duì)雞胸肉中肌原纖維蛋白凝膠品質(zhì)的影響
    NaCl濃度對(duì)肌原纖維蛋白-食用膠混合物功能特性的影響
    氧化程度對(duì)肌原纖維蛋白熱凝膠及理化特性的影響
    麻豆国产av国片精品| 亚洲 国产 在线| 中文资源天堂在线| 日日夜夜操网爽| 成熟少妇高潮喷水视频| 亚洲国产欧洲综合997久久,| 国产精品一区二区三区四区久久| 午夜a级毛片| 我的老师免费观看完整版| 国产精品一区二区三区四区免费观看 | 熟女电影av网| 亚洲精品一卡2卡三卡4卡5卡| 精品午夜福利视频在线观看一区| 国产不卡一卡二| 国内精品一区二区在线观看| 在线播放无遮挡| 欧美3d第一页| 国产男人的电影天堂91| 国产精品久久久久久精品电影| 此物有八面人人有两片| h日本视频在线播放| 人人妻人人看人人澡| 亚洲中文日韩欧美视频| 韩国av一区二区三区四区| 亚洲三级黄色毛片| 亚洲,欧美,日韩| 国产精品野战在线观看| 色5月婷婷丁香| 欧美日本视频| 99热精品在线国产| 亚洲在线观看片| 国产精品不卡视频一区二区| 免费av不卡在线播放| 婷婷精品国产亚洲av| 国产成人av教育| 免费在线观看成人毛片| 亚洲图色成人| 日韩在线高清观看一区二区三区 | eeuss影院久久| 99在线人妻在线中文字幕| 两个人视频免费观看高清| 亚洲av免费高清在线观看| 老司机福利观看| 中国美白少妇内射xxxbb| 一个人看的www免费观看视频| 国产午夜精品久久久久久一区二区三区 | 一进一出抽搐动态| 99久久成人亚洲精品观看| 91在线观看av| 精品不卡国产一区二区三区| 久久人人精品亚洲av| 美女高潮喷水抽搐中文字幕| 欧美黑人巨大hd| 午夜激情欧美在线| 欧美高清性xxxxhd video| 中国美白少妇内射xxxbb| 色哟哟哟哟哟哟| 午夜福利18| 日韩欧美三级三区| 成人午夜高清在线视频| 久久午夜亚洲精品久久| 亚洲欧美激情综合另类| 99热这里只有精品一区| 日本在线视频免费播放| 精品久久久噜噜| 国产欧美日韩精品一区二区| 尤物成人国产欧美一区二区三区| 亚洲成av人片在线播放无| 老熟妇乱子伦视频在线观看| 三级国产精品欧美在线观看| 精品人妻视频免费看| 免费不卡的大黄色大毛片视频在线观看 | 成人无遮挡网站| 极品教师在线免费播放| 一个人观看的视频www高清免费观看| 国模一区二区三区四区视频| 亚洲欧美清纯卡通| 麻豆成人午夜福利视频| 99久国产av精品| 中文在线观看免费www的网站| 国产又黄又爽又无遮挡在线| 亚洲三级黄色毛片| 老熟妇乱子伦视频在线观看| 亚洲色图av天堂| 午夜福利在线观看免费完整高清在 | 日本免费a在线| 91久久精品国产一区二区成人| 欧美成人一区二区免费高清观看| 简卡轻食公司| 桃色一区二区三区在线观看| 精品人妻一区二区三区麻豆 | 一级黄片播放器| 日韩欧美在线二视频| 最近最新中文字幕大全电影3| 免费人成在线观看视频色| 亚洲午夜理论影院| 一区二区三区激情视频| 亚洲 国产 在线| 老熟妇仑乱视频hdxx| 亚洲成人中文字幕在线播放| 91麻豆精品激情在线观看国产| 久久人人爽人人爽人人片va| 琪琪午夜伦伦电影理论片6080| 天美传媒精品一区二区| 中文字幕熟女人妻在线| 日本爱情动作片www.在线观看 | 亚洲美女搞黄在线观看 | 国产一区二区三区视频了| 欧美bdsm另类| 丰满乱子伦码专区| 亚洲黑人精品在线| 18禁裸乳无遮挡免费网站照片| 国产69精品久久久久777片| 国产伦精品一区二区三区四那| 国产精品亚洲一级av第二区| 午夜福利18| 色综合站精品国产| 日本一二三区视频观看| 婷婷丁香在线五月| 久久久久精品国产欧美久久久| 国产精品不卡视频一区二区| 永久网站在线| 婷婷六月久久综合丁香| 丰满乱子伦码专区| 久久久久久久精品吃奶| 免费搜索国产男女视频| 亚洲欧美清纯卡通| 久久草成人影院| 亚洲三级黄色毛片| 亚洲av五月六月丁香网| 国产三级中文精品| 国产精品一及| 日本a在线网址| bbb黄色大片| 淫妇啪啪啪对白视频| 免费无遮挡裸体视频| 国产男靠女视频免费网站| 久久精品影院6| 亚洲无线在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲经典国产精华液单| 精品人妻1区二区| 淫秽高清视频在线观看| 给我免费播放毛片高清在线观看| 中文字幕免费在线视频6| 色5月婷婷丁香| 日本与韩国留学比较| 丰满的人妻完整版| 国产国拍精品亚洲av在线观看| 欧美色视频一区免费| 精品久久久噜噜| 日本黄大片高清| 欧美性猛交黑人性爽| 麻豆精品久久久久久蜜桃| 国产精品不卡视频一区二区| 日本免费一区二区三区高清不卡| 国产欧美日韩精品亚洲av| 亚洲va在线va天堂va国产| 一卡2卡三卡四卡精品乱码亚洲| 小蜜桃在线观看免费完整版高清| 国产乱人伦免费视频| 大型黄色视频在线免费观看| 久久久成人免费电影| 一夜夜www| 精品午夜福利视频在线观看一区| 午夜免费男女啪啪视频观看 | 日本免费一区二区三区高清不卡| 欧美性感艳星| 日本撒尿小便嘘嘘汇集6| 18+在线观看网站| 亚洲精品在线观看二区| 狂野欧美激情性xxxx在线观看| 精品人妻熟女av久视频| 99热这里只有是精品50| 国产高潮美女av| 亚洲熟妇中文字幕五十中出| 好男人在线观看高清免费视频| 淫妇啪啪啪对白视频| 久久久精品大字幕| 免费观看在线日韩| 婷婷六月久久综合丁香| 中国美女看黄片| 国产av一区在线观看免费| 日韩一区二区视频免费看| 91久久精品电影网| 免费高清视频大片| 亚洲欧美清纯卡通| 乱码一卡2卡4卡精品| 99精品在免费线老司机午夜| 色综合色国产| 最新中文字幕久久久久| 久久久久久大精品| 日韩欧美精品免费久久| 色av中文字幕| 大又大粗又爽又黄少妇毛片口| 亚洲专区国产一区二区| av国产免费在线观看| 俄罗斯特黄特色一大片| 在线免费十八禁| 真实男女啪啪啪动态图| 亚洲专区国产一区二区| 日本与韩国留学比较| 亚洲精品国产成人久久av| 能在线免费观看的黄片| 最近最新中文字幕大全电影3| 白带黄色成豆腐渣| 搡老岳熟女国产| 日本免费a在线| 99热6这里只有精品| 成人综合一区亚洲| 人妻夜夜爽99麻豆av| 免费大片18禁| 国产精品一区二区免费欧美| 永久网站在线| АⅤ资源中文在线天堂| 美女高潮的动态| 亚洲精华国产精华液的使用体验 | 在线天堂最新版资源| 俄罗斯特黄特色一大片| 成人美女网站在线观看视频| 日本色播在线视频| 99在线视频只有这里精品首页| 搞女人的毛片| 亚洲精品粉嫩美女一区| 亚洲第一电影网av| 久久人人精品亚洲av| 久久久久九九精品影院| 欧美绝顶高潮抽搐喷水| 有码 亚洲区| 蜜桃久久精品国产亚洲av| 久久久久久九九精品二区国产| 欧美zozozo另类| 成人综合一区亚洲| 草草在线视频免费看| 亚洲成人久久爱视频| 免费av不卡在线播放| 九九热线精品视视频播放| 亚洲美女搞黄在线观看 | 一级a爱片免费观看的视频| 91久久精品国产一区二区三区| av福利片在线观看| 无遮挡黄片免费观看| 乱人视频在线观看| 成熟少妇高潮喷水视频| 欧美三级亚洲精品| 露出奶头的视频| 91精品国产九色| 黄色一级大片看看| 日日啪夜夜撸| 欧美性猛交╳xxx乱大交人| 一卡2卡三卡四卡精品乱码亚洲| 99久久九九国产精品国产免费| 97碰自拍视频| 成年女人永久免费观看视频| 国内久久婷婷六月综合欲色啪| 99国产极品粉嫩在线观看| 午夜福利视频1000在线观看| 欧美激情国产日韩精品一区| 国产成年人精品一区二区| 色av中文字幕| 99国产极品粉嫩在线观看| 国产亚洲精品久久久com| 黄色一级大片看看| av在线天堂中文字幕| 亚洲在线观看片| 成人国产麻豆网| 亚洲成a人片在线一区二区| 夜夜看夜夜爽夜夜摸| 两个人视频免费观看高清| 精品久久久噜噜| 国产高清激情床上av| 美女高潮的动态| 一进一出抽搐gif免费好疼| 一边摸一边抽搐一进一小说| 天天一区二区日本电影三级| 黄色视频,在线免费观看| 亚洲三级黄色毛片| 免费电影在线观看免费观看| 女人被狂操c到高潮| 1000部很黄的大片| 性插视频无遮挡在线免费观看| 国产伦在线观看视频一区| 床上黄色一级片| 午夜久久久久精精品| eeuss影院久久| 在线看三级毛片| 亚洲中文字幕一区二区三区有码在线看| 啪啪无遮挡十八禁网站| 日韩国内少妇激情av| 久久精品综合一区二区三区| 中文字幕久久专区| 我的女老师完整版在线观看| 欧美高清性xxxxhd video| 99国产极品粉嫩在线观看| 亚洲最大成人手机在线| 亚洲人成网站在线播放欧美日韩| 99国产精品一区二区蜜桃av| 搞女人的毛片| 中文字幕熟女人妻在线| 亚洲av中文字字幕乱码综合| 88av欧美| 九色国产91popny在线| 国产高清激情床上av| 免费看日本二区| 人人妻,人人澡人人爽秒播| 亚洲欧美日韩高清专用| 亚洲电影在线观看av| 嫩草影院入口| 成人国产综合亚洲| 精品免费久久久久久久清纯| 韩国av一区二区三区四区| 国产一区二区在线av高清观看| 观看免费一级毛片| 国产精品久久久久久精品电影| 亚洲av不卡在线观看| 亚洲乱码一区二区免费版| 51国产日韩欧美| 999久久久精品免费观看国产| 中亚洲国语对白在线视频| 成人高潮视频无遮挡免费网站| 国产精品爽爽va在线观看网站| 99热精品在线国产| 国产男人的电影天堂91| 亚洲人成网站高清观看| 色吧在线观看| 国产综合懂色| 亚洲一级一片aⅴ在线观看| 亚洲图色成人| 精品久久久久久久久久免费视频| 日日撸夜夜添| 在线播放国产精品三级| 国产伦精品一区二区三区视频9| 成人二区视频| 国产一区二区三区av在线 | 成年女人永久免费观看视频| 成人欧美大片| 天堂√8在线中文| 真实男女啪啪啪动态图| 在线免费观看的www视频| av黄色大香蕉| 国产精品一及| 国产 一区 欧美 日韩| 久久精品夜夜夜夜夜久久蜜豆| 免费看av在线观看网站| 村上凉子中文字幕在线| 97超级碰碰碰精品色视频在线观看| 久久热精品热| 国产精品一区www在线观看 | 亚洲欧美精品综合久久99| 在线观看av片永久免费下载| 色播亚洲综合网| 少妇的逼水好多| 成熟少妇高潮喷水视频| 精品久久久久久久人妻蜜臀av| 男人舔女人下体高潮全视频| 在线观看午夜福利视频| 天堂av国产一区二区熟女人妻| 热99在线观看视频| 午夜爱爱视频在线播放| 一个人看的www免费观看视频| 久久精品国产清高在天天线| 久久久色成人| 久久久久久久久久成人| 亚洲av免费在线观看| 99热这里只有是精品在线观看| 黄色日韩在线| 国产一区二区三区在线臀色熟女| 直男gayav资源| 欧美又色又爽又黄视频| 成人无遮挡网站| 搡老熟女国产l中国老女人| 丰满人妻一区二区三区视频av| 欧美色欧美亚洲另类二区| 国产老妇女一区| 免费看美女性在线毛片视频| 99热这里只有是精品50| 亚洲精品成人久久久久久| 午夜爱爱视频在线播放| 三级男女做爰猛烈吃奶摸视频| 国产精品精品国产色婷婷| 九色国产91popny在线| 精品国产三级普通话版| 亚洲不卡免费看| 久久婷婷人人爽人人干人人爱| 亚洲不卡免费看| 久久婷婷人人爽人人干人人爱| 免费黄网站久久成人精品| 精品久久久久久久久亚洲 | 精品人妻1区二区| 97超级碰碰碰精品色视频在线观看| av在线观看视频网站免费| 两人在一起打扑克的视频| 亚洲真实伦在线观看| 很黄的视频免费| 亚洲真实伦在线观看| 欧美丝袜亚洲另类 | 亚洲不卡免费看| 国产免费av片在线观看野外av| 久久久久久久久久久丰满 | 精品日产1卡2卡| 美女免费视频网站| 直男gayav资源| 我的女老师完整版在线观看| 免费在线观看日本一区| 露出奶头的视频| 中文资源天堂在线| 岛国在线免费视频观看| 亚洲专区国产一区二区| 18禁裸乳无遮挡免费网站照片| av.在线天堂| 男人和女人高潮做爰伦理| 国产精品一区二区性色av| 乱系列少妇在线播放| 成人国产一区最新在线观看| 久久久国产成人精品二区| 日本色播在线视频| 久久久久精品国产欧美久久久| 精品乱码久久久久久99久播| 国产精品久久电影中文字幕| 12—13女人毛片做爰片一| 亚洲五月天丁香| 成人高潮视频无遮挡免费网站| 欧美在线一区亚洲| 国产视频一区二区在线看| 一本一本综合久久| 人人妻,人人澡人人爽秒播| 国产单亲对白刺激| 国产伦精品一区二区三区视频9| 国产aⅴ精品一区二区三区波| 亚洲熟妇中文字幕五十中出| 91在线观看av| 91久久精品国产一区二区成人| 国产老妇女一区| 最近在线观看免费完整版| 中文字幕av成人在线电影| 日本黄色片子视频| 婷婷六月久久综合丁香| 黄色欧美视频在线观看| 欧美性猛交╳xxx乱大交人| 人妻制服诱惑在线中文字幕| 久久午夜福利片| 色哟哟哟哟哟哟| 很黄的视频免费| www日本黄色视频网| 高清在线国产一区| 变态另类丝袜制服| 国产精品99久久久久久久久| 欧美成人一区二区免费高清观看| 成人av一区二区三区在线看| 成人鲁丝片一二三区免费| 国产高清激情床上av| 国产在线精品亚洲第一网站| 亚洲狠狠婷婷综合久久图片| 欧美精品啪啪一区二区三区| 日韩高清综合在线| 嫩草影视91久久| 国产v大片淫在线免费观看| 久久亚洲精品不卡| 大又大粗又爽又黄少妇毛片口| 久久久精品欧美日韩精品| 国产精品无大码| 国产 一区精品| 啦啦啦啦在线视频资源| 日本与韩国留学比较| 中文字幕精品亚洲无线码一区| 22中文网久久字幕| 精品一区二区三区视频在线| 国产一级毛片七仙女欲春2| 久久久久九九精品影院| 性欧美人与动物交配| 成人毛片a级毛片在线播放| 国产伦人伦偷精品视频| 老司机深夜福利视频在线观看| 色哟哟·www| 欧美国产日韩亚洲一区| 黄色女人牲交| 日本黄色片子视频| 不卡一级毛片| 99久久精品国产国产毛片| 欧美最新免费一区二区三区| 免费观看人在逋| 国产伦精品一区二区三区四那| 五月玫瑰六月丁香| 国产精品一区二区免费欧美| 日日撸夜夜添| 国产在视频线在精品| 91av网一区二区| 老司机午夜福利在线观看视频| 日韩高清综合在线| 欧美不卡视频在线免费观看| 欧美成人性av电影在线观看| 国产精品99久久久久久久久| 偷拍熟女少妇极品色| 欧美高清成人免费视频www| 成年女人永久免费观看视频| 看片在线看免费视频| 麻豆国产97在线/欧美| 国产探花在线观看一区二区| 亚洲av成人av| 精品久久久噜噜| 国产国拍精品亚洲av在线观看| 国产高清三级在线| 亚洲18禁久久av| 国产 一区 欧美 日韩| 精品一区二区免费观看| 少妇裸体淫交视频免费看高清| 乱系列少妇在线播放| 亚洲人成网站在线播放欧美日韩| 中文字幕熟女人妻在线| 国产精品不卡视频一区二区| 超碰av人人做人人爽久久| 欧美精品啪啪一区二区三区| 少妇裸体淫交视频免费看高清| 国产成人a区在线观看| 香蕉av资源在线| 国产精品一区二区免费欧美| 琪琪午夜伦伦电影理论片6080| 亚洲午夜理论影院| 亚洲性夜色夜夜综合| 日日摸夜夜添夜夜添av毛片 | 亚洲在线观看片| 琪琪午夜伦伦电影理论片6080| 日韩在线高清观看一区二区三区 | 久久久国产成人精品二区| 2021天堂中文幕一二区在线观| 欧美xxxx性猛交bbbb| 国产高清视频在线播放一区| 日韩 亚洲 欧美在线| 一区二区三区高清视频在线| 中亚洲国语对白在线视频| 亚洲av二区三区四区| 国产精品久久视频播放| 在线观看66精品国产| 高清毛片免费观看视频网站| 日韩中文字幕欧美一区二区| 亚洲欧美清纯卡通| 国产午夜福利久久久久久| 亚洲av免费高清在线观看| 12—13女人毛片做爰片一| 久99久视频精品免费| 此物有八面人人有两片| 中文字幕精品亚洲无线码一区| 精品人妻1区二区| 免费看av在线观看网站| www日本黄色视频网| 99热这里只有精品一区| 久久99热这里只有精品18| 午夜久久久久精精品| 免费电影在线观看免费观看| 精品久久久久久久久久免费视频| 国产欧美日韩精品亚洲av| 麻豆一二三区av精品| 狂野欧美激情性xxxx在线观看| 色5月婷婷丁香| 一个人免费在线观看电影| 国产毛片a区久久久久| 久久久久国内视频| 男人舔奶头视频| 又爽又黄a免费视频| 最近视频中文字幕2019在线8| 国产高清不卡午夜福利| 国内精品久久久久久久电影| 丰满的人妻完整版| 亚洲性久久影院| 亚洲av免费高清在线观看| eeuss影院久久| 国产精品野战在线观看| 欧美成人一区二区免费高清观看| 2021天堂中文幕一二区在线观| 最近最新中文字幕大全电影3| 制服丝袜大香蕉在线| 国产免费一级a男人的天堂| 午夜福利成人在线免费观看| 亚洲av日韩精品久久久久久密| av福利片在线观看| 国产69精品久久久久777片| 亚洲熟妇中文字幕五十中出| 黄色一级大片看看| 国产中年淑女户外野战色| 99热网站在线观看| 伦精品一区二区三区| 亚洲最大成人手机在线| 男人舔女人下体高潮全视频| 中文字幕av在线有码专区| 久久亚洲真实| 婷婷精品国产亚洲av在线| 老女人水多毛片| 日本a在线网址| 大型黄色视频在线免费观看| 尤物成人国产欧美一区二区三区| 18禁在线播放成人免费| 免费在线观看影片大全网站| 九九爱精品视频在线观看| 国产精品精品国产色婷婷| 身体一侧抽搐| 亚洲美女搞黄在线观看 | 日本黄色片子视频| 长腿黑丝高跟| h日本视频在线播放| 久9热在线精品视频| 国产高清视频在线播放一区| 欧美3d第一页| 久久久久久伊人网av| 18+在线观看网站| 日本欧美国产在线视频| 男女下面进入的视频免费午夜| 亚洲午夜理论影院| 亚洲欧美日韩高清专用| 久久这里只有精品中国| 欧美+日韩+精品| 国产精品美女特级片免费视频播放器| 免费av毛片视频| 一个人看视频在线观看www免费| 麻豆国产av国片精品| 国产蜜桃级精品一区二区三区|