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

    貴州地方白羽雞種MSTN基因SNP位點快速篩查及其蛋白功能預(yù)測

    2020-11-02 02:41龍廣麗牟騰慧劉洋饒永超楊德鳳譚光輝林家棟傅筑蔭張福平
    南方農(nóng)業(yè)學(xué)報 2020年8期
    關(guān)鍵詞:多態(tài)性生物信息學(xué)

    龍廣麗 牟騰慧 劉洋 饒永超 楊德鳳 譚光輝 林家棟 傅筑蔭 張福平

    摘要:【目的】明確貴州地方白羽雞種肌肉生長抑制素(MSTN)基因的多態(tài)性,篩選出改良白羽雞肉質(zhì)品質(zhì)的分子遺傳標(biāo)記,為選擇培育優(yōu)良肉質(zhì)雞種提供參考依據(jù)?!痉椒ā恳耘d義矮腳雞白羽品系、羅曼蛋雞及興義矮腳雞白羽品系與羅曼蛋雞雜交的自交后代(F2)為研究對象,提取3個雞種的血液基因組DNA構(gòu)建DNA混合池,PCR擴增MSTN基因全部外顯子序列,直接測序篩選SNP位點并進行生物信息學(xué)分析,探究SNP位點突變對MSTN基因mRNA二級結(jié)構(gòu)及編碼蛋白結(jié)構(gòu)和功能的影響?!窘Y(jié)果】在雞MSTN基因第1外顯子(Exon-1)篩查到6個SNPs位點,分別為G51A、A60G、G195C、A234G、C297T和C324T,且6個SNPs位點均為同義突變。在A234G位點處檢測到興義矮腳雞白羽系和F2代存在多態(tài)性,但羅曼蛋雞未發(fā)生變異;在C297T位點處檢測到羅曼蛋雞存在多態(tài)性,以C為優(yōu)勢等位基因,而興義矮腳雞白羽品系和F2代均未發(fā)生變異;G51A、G60A、G159C和C324T等4個SNPs位點在3個雞種中均存在多態(tài)性。A234G位點突變對雞MSTN基因mRNA二級結(jié)構(gòu)穩(wěn)定性無影響,G51A位點突變致使MSTN基因mRNA二級結(jié)構(gòu)最小自由能降低但未改變其構(gòu)象,A60G、G195C、C297T和C324T等4個SNPs位點突變均引起MSTN基因mRNA二級結(jié)構(gòu)和最小自由能的改變。雞MSTN基因編碼蛋白為不穩(wěn)定的脂溶性蛋白,存在信號肽,但不存在跨膜區(qū)域,其二級結(jié)構(gòu)由無規(guī)則卷曲、β-折疊和α-螺旋組成,且以無規(guī)則卷曲占比最高?!窘Y(jié)論】MSTN基因在貴州地方白羽雞種中的多態(tài)性較豐富,其中A60G、G195C、C297T和C324T等4個SNPs位點突變會導(dǎo)致基因mRNA二級結(jié)構(gòu)的構(gòu)象發(fā)生改變,可作為改良白羽雞肉質(zhì)品質(zhì)的分子遺傳標(biāo)記進一步研究。

    關(guān)鍵詞: 貴州地方白羽雞種;MSTN基因;肉質(zhì)性狀;SNP位點;多態(tài)性;生物信息學(xué)

    中圖分類號: S831.89? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻標(biāo)志碼: A 文章編號:2095-1191(2020)08-1840-09

    Rapid screening and protein function prediction of MSTN gene SNP sites in Guizhou local chicken breeds

    LONG Guang-li1,2,3, MOU Teng-hui1,2,3, LIU Yang1,2,3, RAO Yong-chao1,2,3,

    YANG De-feng1,2,3, TAN Guang-hui1,2, LIN Jia-dong1,2, FU Zhu-yin2,3,

    ZHANG Fu-ping1,2,3*

    (1College of Animal Science, Guizhou University/Key Laboratory of Animal Genetics, Breeding and Reproduction in the Plateau Mountainous Region, Ministry of Education/Key Laboratory of Animal Genetics, Breeding and Reproduction in Guizhou, Guiyang? 550025, China; 2Institute of Poultry, Guizhou University, Guiyang? 550025, China;

    3Scientific Research Chicken Farm of Guizhou University, Guiyang? 550025, China)

    Abstract:【Objective】To identify the polymorphism of the myostatin (MSTN) gene of Guizhou local white feather chicken breeds, screen out ideal molecular genetic markers which could improve the meat quality of white feather chi-ckens, and provide a reference for selecting and breeding excellent broiler chicken breeds. 【Method】Taking Xingyi bantam white feather line, Roman laying chicken and self-progeny(F2) of cross between Xingyi bantam white feather line and Roman laying chicken as the research object, the blood genomic DNA of three chicken breeds was extracted to construct a DNA mixing pool. The PCR amplification of all exon sequences of MSTN gene was directly sequenced and screened. Then the SNP sites were bioinformatically analyzed to explore the impact of SNP site mutations on the se-condary structure of MSTN gene mRNA and the structure and function of the encoded protein. 【Result】Six SNPs sites were screened in the first exon(Exon-1) of the chicken MSTN gene, which were G51A, A60G, G195C, A234G, C297T and C324T, and the six SNPs sites were all synonymous mutations. At the A234G locus, polymorphisms were detected in the Xingyi bantam white feather line and F2 generation, but no mutations occurred in the Roman laying chickens; polymorphisms were detected in the Roman laying chickens at the C297T locus, with C as the dominant position. There were no mutations in Xingyi bantam white feather line and F2 generation. At four SNPs, including G51A, G60A, G159C and C324T, polymorphisms were all detected in three chicken breeds. A234G site mutation had no effect on the stability of the secondary structure of chicken MSTN gene mRNA. G51A site mutation reduced the minimum free energy of mRNA secon-dary structure but did not change its conformation. Four SNPs site mutations including A60G, G195C, C297T and C324T caused changes in the secondary structure and minimum free energy of MSTN gene mRNA. The chicken MSTN gene encoding protein was an unstable fat-soluble protein with a signal peptide but no transmembrane regions. Its secondary structure consisted of random coils, β-fold and α-helixes, with random coils accounting for the highest proportion. 【Conclusion】The MSTN gene has abundant polymorphisms in Guizhou local white feather chicken breeds. Among them, four SNPs mutations such as A60G, G195C, C297T and C324T cause changes in the conformation of the gene mRNA secondary structure, which can be used to do further studies as molecular genetic markers improving the quality of white feather chicken.

    Key words: Guizhou local white feather chicken breeds; MSTN gene; meat quality traits; SNP site; polymorphism; bioinformatics

    Foundation item: Guizhou Science and Technology Planning Project(QKHZC〔2016〕2507)

    0 引言

    【研究意義】興義矮腳雞白羽系主產(chǎn)于貴州黔西南州興義市,其脛短身矮、體軀勻稱、全身羽毛白色,為隱性白羽雞,在生產(chǎn)性能和肉質(zhì)方面具有許多優(yōu)良特性,但產(chǎn)蛋性能不高(朱麗莉等,2011;李林洪,2019)。羅曼蛋雞的產(chǎn)蛋率高,蛋品質(zhì)優(yōu)良(劉雅正,2014),但至今針對其肉質(zhì)方面的相關(guān)研究甚少。優(yōu)勢雜交可顯著提高畜禽的生產(chǎn)性能和經(jīng)濟效益(王芬和王鵬,2009),目前國內(nèi)隱性白羽雞在雜交育種中主要用作母本(李俊英,2005)。選用興義矮腳雞白羽系作為配套母本與羅曼蛋雞雜交,其后代具有較高的生產(chǎn)性能,但在生產(chǎn)性能改良過程中興義矮腳雞白羽系雜交后代的體型發(fā)生明顯變化,因此,亟待篩選出白羽雞肉質(zhì)品質(zhì)的理想分子遺傳標(biāo)記,在不改變地方雞種優(yōu)質(zhì)性能的前提下有效提高白羽雞的產(chǎn)肉性能。【前人研究進展】肌肉生長抑制素(Myostatin,MSTN)又稱為生長分化因子8(Growth differentiation factor-8,GDF-8),可抑制成肌細胞的活化、增殖和分化(McFarland et al.,2007)。MSTN基因是轉(zhuǎn)化生長因子β(TGF-β)超家族成員之一,其特異性表達于畜禽肌肉組織中(McPherron et al.,1997)。MSTN基因作為胞外信號分子能與成肌細胞膜上的受體結(jié)合而引起受體自身磷酸化,啟動細胞內(nèi)一系列信號傳導(dǎo)途徑,作用于生肌決定因子靶基因的調(diào)控區(qū),調(diào)節(jié)肌肉組成蛋白基因表達,進而對骨骼肌的生長發(fā)育起負調(diào)控作用(Lee,2004;Grade et al.,2009;趙義龍等,2017)。MSTN基因突變或缺失可使家畜表現(xiàn)出肌肉超級發(fā)達的雙肌性狀(鄭云,2007;楊永生等,2012;趙振華等,2015),在負向調(diào)控肌細胞的過程中若MSTN活性喪失或減弱,則無法正常調(diào)節(jié)肌纖維沉積,導(dǎo)致動物肌肉過度發(fā)育(Grobet et al.,1997;Bellinge et al.,2005;張賢嫻等,2018)。近年來的相關(guān)研究顯示,MSTN基因作為畜禽生長發(fā)育及肉用性狀相關(guān)的候選基因,其在胚胎發(fā)育過程中的免疫中和作用可促進肌肉生長(Kim et al.,2006),同時參與脂肪代謝過程,抑制脂肪形成(Guo et al.,2008),進而調(diào)節(jié)動物肌肉與脂肪的比例,在改善畜禽肉質(zhì)方面具有極高的研究價值(曹婷等,2017)。Fontanesi等(2011)研究發(fā)現(xiàn)家兔MSTN基因SNP位點突變與其肝臟和胴體量相關(guān);Xu等(2013)研究表明北京鴨MSTN基因多態(tài)性與其胸肌肌肉厚度顯著相關(guān);Zhang等(2015)研究發(fā)現(xiàn),MSTN基因c.234G>A突變位點可作為輔助育種分子標(biāo)記,以提高雞的繁殖效率;Mitrofanova等(2017)在分析雞MSTN基因SNP多態(tài)性時發(fā)現(xiàn),AA基因型個體的活重明顯高于其他基因型個體。可見,MSTN基因是繁育優(yōu)良肉質(zhì)雞種的一個重要候選基因?!颈狙芯壳腥朦c】雞MSTN基因全長5493 bp,包含3個外顯子(劉雪穎,2008),但至今有關(guān)白羽雞MSTN基因多態(tài)性及其蛋白結(jié)構(gòu)的研究鮮見報道?!緮M解決的關(guān)鍵問題】以興義矮腳雞白羽品系、羅曼蛋雞及興義矮腳雞白羽品系與羅曼蛋雞雜交的自交后代(F2)為研究對象,構(gòu)建DNA混合池,PCR擴增MSTN基因全部外顯子序列,直接測序篩選SNP位點并進行生物信息學(xué)分析,明確SNP位點突變對MSTN基因mRNA二級結(jié)構(gòu)及編碼蛋白結(jié)構(gòu)和功能的影響,旨在篩選出改良白羽雞肉質(zhì)性狀的分子遺傳標(biāo)記,為選擇培育優(yōu)良肉質(zhì)雞種提供參考依據(jù)。

    1 材料與方法

    1. 1 試驗材料

    興義矮腳雞白羽品系(B)、羅曼蛋雞(L)及F2代各60份血液樣品,均采自貴州大學(xué)科研雞場,所有血液樣品均置于肝素鈉抗凝真空采血管中,-20 ℃保存?zhèn)溆?。血液基因組DNA提取試劑盒、GoldView核酸染料、DL2000 DNA Marker和2×Taq Master Mix均購自生工生物工程(上海)股份有限公司;TBE緩沖液、雙蒸水和瓊脂糖等試劑購自上海鼎國生物技術(shù)有限公司。

    1. 2 基因組DNA提取

    利用血液基因組DNA提取試劑盒提取血液樣品基因組DNA,以1.0%瓊脂糖凝膠電泳檢測提取效果,用紫外分光光度計測量各DNA樣品濃度,測量3次取平均值,并將所有DNA樣品濃度調(diào)至100 ng/μL,各取10.0 μL構(gòu)建不同雞種的DNA混合池。

    1. 3 引物設(shè)計及PCR擴增

    根據(jù)GenBank已公布的雞MSTN基因序列(登錄號NC_006094.5),運用Primer Premier 5.0設(shè)計引物(表1),并委托生工生物工程(上海)股份有限公司合成。PCR反應(yīng)體系25.0 μL:2×Taq Master Mix 12.5 μL,上、下游引物各1.0 μL,DNA模板1.0 μL,雙蒸水9.5 μL。擴增程序:94 ℃預(yù)變性3 min;94 ℃ 30 s,60 ℃ 30 s,72 ℃ 30 s,進行35個循環(huán);72 ℃延伸10 min,4 ℃保存。PCR擴增產(chǎn)物用1.0%瓊脂糖凝膠電泳進行檢測。

    1. 4 序列校對及等位基因頻率估算

    利用DNASTAR中的MegAlign和EditSeq程序結(jié)合測序峰圖,拼接得到白羽雞種MSTN基因編碼區(qū)(CDS)序列,然后與GenBank已公布的雞MSTN基因CDS序列進行比對,鑒定篩選出SNP位點。運用SeqMan分析測序峰圖,采用MWSnap v3.00.74中的標(biāo)尺測量各SNP位點等位基因的相應(yīng)峰高,并按照崔建勛等(2005)、祖盤玉等(2019)的方法估算等位基因頻率。

    fi=hi/(h1+h2

    式中,fi表示SNP位點某等位基因頻率,h1和h2分別表示測序峰圖某等位基因1峰和2峰的高度。

    1. 5 生物信息學(xué)分析

    參照翁吉梅等(2018)的研究方法對3個雞種的MSTN基因進行生物信息學(xué)分析。采用RNAfold和RNA structure預(yù)測不同SNP位點及其突變前后對MSTN基因mRNA二級結(jié)構(gòu)的影響,并以PredictProtein預(yù)測SNP位點突變前后的MSTN蛋白二級結(jié)構(gòu)。

    2 結(jié)果與分析

    2. 1 PCR擴增結(jié)果

    PCR擴增得到雞MSTN基因的全部外顯子序列(Exon-1、Exon-2和Exon-3),經(jīng)1.0%瓊脂糖凝膠電泳檢測顯示,電泳條帶單一、明亮清晰,無拖帶等現(xiàn)象,說明擴增效果良好,且片段大小與預(yù)期結(jié)果相符。圖1為雞MSTN基因Exon-1的PCR擴增電泳檢測結(jié)果。

    2. 2 測序分析結(jié)果

    以雞MSTN基因(登錄號NC_006094.5)為參考序列,采用DNA混合池(每個品種60份血液樣品,各取1.0 μL)結(jié)合PCR直接測序法對雞MSTN基因外顯子序列上的SNP位點進行鑒定,再運用DNASTAR等在線軟件對所得序列進行比對分析,結(jié)果在雞MSTN基因Exon-1上篩查到6個SNPs位點(圖2),分別是G51A、A60G、G195C、A234G、C297T和C324T,且6個SNPs位點均為同義突變。

    2. 3 等位基因頻率估算結(jié)果

    利用MWSnap v3.00.74中的標(biāo)尺測量不同雞種MSTN基因6個SNPs位點的等位基因峰值,并估算等位基因頻率,結(jié)果(表2)顯示,在A234G位點處檢測到興義矮腳雞白羽系和F2代存在多態(tài)性,但羅曼蛋雞未發(fā)生變異;在G159C位點處檢測到羅曼蛋雞的G和C基因頻率差異不明顯,而在興義矮腳雞白羽系和F2代均以G為優(yōu)勢等位基因;在C297T位點處檢測到羅曼蛋雞存在多態(tài)性,以C為優(yōu)勢等位基因,而興義矮腳雞白羽品系和F2代均未發(fā)生變異。G51A、A60G、G159C和C324T等4個SNPs位點在3個雞種中均存在多態(tài)性。

    2. 4 生物信息學(xué)分析結(jié)果

    2. 4. 1 mRNA二級結(jié)構(gòu)預(yù)測分析結(jié)果 采用RNAfold預(yù)測不同雞種MSTN基因SNP位點突變前后各基因型的mRNA二級結(jié)構(gòu),結(jié)果如圖3所示。A234G位點突變未引起MSTN基因mRNA二級結(jié)構(gòu)構(gòu)象及最小自由能的改變,即對其穩(wěn)定性無影響;A60G、G195C和C297T位點突變致使MSTN基因mRNA二級結(jié)構(gòu)的最小自由能分別增加至-282.90、 -281.90和-281.20 kcal/mol,結(jié)構(gòu)穩(wěn)定性降低;而G51A和C324T位點突變致使MSTN基因mRNA二級結(jié)構(gòu)的最小自由能分別降低為-284.30和-286.10 kcal/mol,結(jié)構(gòu)穩(wěn)定性增強。此外,A60G、G195C、C297T和C324T位點突變會導(dǎo)致MSTN基因mRNA二級結(jié)構(gòu)的構(gòu)象發(fā)生改變,但G51A和A234G位點突變未改變MSTN基因mRNA二級結(jié)構(gòu)的構(gòu)象。

    2. 4. 2 MSTN蛋白理化性質(zhì)預(yù)測分析結(jié)果 ProtParam預(yù)測分析結(jié)果表明,雞MSTN蛋白由375個氨基酸殘基構(gòu)成,其相對分子量為42707.04 Da,總原子數(shù)為5983,分子式為C1903H2986N510O562S22,理論等電點(pI)為6.93,正、負電荷氨基酸殘基均有46個。在含有的常見20種氨基酸中,所占比例最高的是亮氨酸(Leu),占總氨基酸的8.0%;所占比例最低的是組氨酸(His),占1.3%;無吡咯賴氨酸(Pyl)和硒半胱氨酸(Sec)。雞MSTN蛋白不穩(wěn)定系數(shù)為43.68,屬于不穩(wěn)定性蛋白;其平均親水性為(-0.392),屬于脂溶性蛋白。

    2. 4. 3 信號肽預(yù)測分析結(jié)果 采用SignalP 4.0對雞MSTN蛋白信號肽進行預(yù)測,結(jié)果(圖4)表明,C分值為0.124,Y分值為0.260,S分值為0.560(若S分值平均值高于0.500,則為分泌蛋白,且存在信號肽),說明雞MSTN蛋白存在信號肽,屬于分泌蛋白。

    2. 4. 4 MSTN蛋白跨膜結(jié)構(gòu)預(yù)測結(jié)果 以TMHMM 2.0預(yù)測分析雞MSTN蛋白跨膜結(jié)構(gòu),結(jié)果發(fā)現(xiàn)組成MSTN蛋白的375個氨基酸殘基均在膜外(圖5),不存在跨膜區(qū)域。

    2. 4. 5 MSTN蛋白亞細胞定位分析結(jié)果 PSORT II Prediction的蛋白亞細胞定位結(jié)果表明,雞MSTN蛋白在細胞質(zhì)(47.8%)、細胞核(21.7%)、細胞膜(13.0%)和細胞骨架(4.3%)中發(fā)揮生物學(xué)作用。

    2. 4. 6 MSTN蛋白二級結(jié)構(gòu)分析結(jié)果 SOPMA預(yù)測結(jié)果表明,在雞MSTN蛋白二級結(jié)構(gòu)中所占比例最高的是無規(guī)則卷曲(圖6),屬于混合型蛋白。其中,α-螺旋主要形成于第2~18、40~42、46~63、74~79、86~96、111~112、209~217、290~293和327~328位氨基酸殘基;β-折疊在第115~120、136~141、150~153、155~160、169~179、187~195、203~207、225~231、238~240、252~257、280~289、296~299、304~305、307~311、325~326、343~345、349~353、359~364和368~371位氨基酸殘基處分布較多;無規(guī)則卷曲則在第1~2、19~39、43~45、64~73、80~85、97~110、113~114、121~135、142~149、154、161~168、180~186、196~202、208、218~224、232~237、241~251、258~279、294~295、300~303、306、312~324、329~342、346~348、354~358、365~367和372~375位氨基酸殘基處分布較多。

    3 討論

    在雞的肌肉發(fā)生、生長和發(fā)育過程中,MSTN基因發(fā)揮著重要作用,其CDS序列具有豐富多態(tài)性,對雞的生長性狀有顯著影響(Dushyanth et al.,2016)。張冉和林浴霜(2012)研究發(fā)現(xiàn),在雞胚發(fā)育過程中不同雞品種及不同發(fā)育階段的MSTN基因表達量存在明顯差異,說明MSTN對不同雞品種肌肉發(fā)育的調(diào)節(jié)能力存在差異,而導(dǎo)致肌纖維數(shù)量及肌肉含量不同。左斌等(2015)研究表明,MSTN基因主要影響雞的生長性狀,可作為地方雞品種選育的候選基因。王彥欽和邵勇鋼(2017)研究發(fā)現(xiàn),MSTN基因多態(tài)性與拜城油雞體重、屠體重、半凈膛重、全凈膛重、胸肌重和腿肌重等生產(chǎn)性能指標(biāo)密切相關(guān),故推測該基因可作為拜城油雞選育的分子遺傳標(biāo)記。此外,有研究發(fā)現(xiàn)MSTN基因多態(tài)性對12周齡雞屠宰及肉質(zhì)性能的各項指標(biāo)均有影響,說明MSTN基因不僅在骨骼肌的生長發(fā)育過程中發(fā)揮作用,還可能參與脂質(zhì)代謝和脂肪沉積(顧志良等,2003)。MSTN基因是提高畜禽肌肉產(chǎn)量和改善肉品質(zhì)研究的首選基因(王巧巧等,2018),作為肌肉發(fā)育的負向調(diào)控因子,可通過影響MyoD、Myogenin和Myf5等基因控制肌細胞從G1期到G2期及G2期到M期的轉(zhuǎn)變,進而影響肌肉的發(fā)育過程(Feldman et al.,2006)。李維等(2017)研究發(fā)現(xiàn),MSTN基因?qū)Τ嗨疄豕请u肉質(zhì)指標(biāo)中的pH和剪切力具有明顯影響,推測該基因參與調(diào)控赤水烏骨雞肉質(zhì)性狀??梢?,MSTN基因是影響雞肉質(zhì)性狀的主效基因之一,探究MSTN基因SNP位點突變對基因表達的影響,對提高白羽雞的養(yǎng)殖效益具有重要意義。

    本研究在雞MSTN基因Exon-1上發(fā)現(xiàn)G51A、A60G、G195C、A234G、C297T和C324T等6個SNPs位點,均為同義突變,未改變編碼氨基酸序列,由于密碼子的簡并性,同義密碼子的使用頻率存在差異,但對MSTN基因表達是否產(chǎn)生影響還有待進一步探究。此外,在G159C位點檢測到羅曼蛋雞的G和C基因頻率差異不明顯,而在興義矮腳雞白羽系和F2代中均以G為優(yōu)勢等位基因,說明羅曼蛋雞與其他2個雞種在該位點處存在明顯差異。在MSTN基因Exon-2和Exon-3上未發(fā)現(xiàn)SNP位點,與毛亮等(2009)對藏雞MSTN基因、溫彥濤等(2012)對優(yōu)質(zhì)雞S3系MSTN基因的多態(tài)性檢測結(jié)果相似,但與張麗等(2018)對靜寧雞MSTN基因進行SNP位點篩查的結(jié)果相比,G51A、C297T和C324T是新發(fā)現(xiàn)的SNP位點,可能是相對于Exon-1,Exon-2和Exon-3更加保守所致。雞MSTN蛋白理化性質(zhì)預(yù)測結(jié)果顯示,亮氨酸的占比最高,組氨酸的占比最低,且未發(fā)現(xiàn)吡咯賴氨酸和硒半胱氨酸;其不穩(wěn)定系數(shù)為43.68,即MSTN蛋白屬于不穩(wěn)定性蛋白,與任衛(wèi)合等(2018)的研究結(jié)果相似。

    CDS序列上的核苷酸突變可能引起基因mRNA二級結(jié)構(gòu)及其編碼蛋白結(jié)構(gòu)改變,進而影響其生物學(xué)功能(李倩,2013)。采用RNAfold對雞MSTN基因mRNA二級結(jié)構(gòu)的預(yù)測結(jié)果顯示,A234G位點突變未引起MSTN基因mRNA二級結(jié)構(gòu)改變;A60G、G195C和C297T位點突變導(dǎo)致MSTN基因mRNA二級結(jié)構(gòu)的最小自由能增加,mRNA二級結(jié)構(gòu)穩(wěn)定性降低;G51A和C324T位點突變則致使MSTN基因mRNA二級結(jié)構(gòu)的最小自由能降低,mRNA二級結(jié)構(gòu)穩(wěn)定性增強。此外,雞MSTN蛋白二級結(jié)構(gòu)中以無規(guī)則卷曲結(jié)構(gòu)數(shù)量最多,屬于混合型蛋白,與任衛(wèi)合等(2018)、張麗等(2018)的研究結(jié)果相似,即均以無規(guī)則卷曲占比最高。無規(guī)則卷曲的增加易改變蛋白質(zhì)構(gòu)象,進而影響其生物學(xué)功能(鄭海軍等,2013),但由于使用的預(yù)測軟件不同,致使其他結(jié)構(gòu)區(qū)域的占比略有差異,因此雞MSTN蛋白的結(jié)構(gòu)與功能還需進一步探究。

    4 結(jié)論

    MSTN基因在貴州地方白羽雞種中的多態(tài)性較豐富,其中A60G、G195C、C297T和C324T等4個SNPs位點突變會導(dǎo)致基因mRNA二級結(jié)構(gòu)的構(gòu)象發(fā)生改變,可作為改善白羽雞肉質(zhì)品質(zhì)的分子遺傳標(biāo)記進一步研究。

    參考文獻:

    曹婷,周漢林,荀文娟,施力光,周雄,侯冠彧. 2017. MSTN基因?qū)ωi骨骼肌發(fā)育調(diào)控的作用及其研究進展[J]. 基因組學(xué)與應(yīng)用生物學(xué),36(4):1511-1517. [Cao T,Zhou H L,Xun W J,Shi L G,Zhou X,Hou G Y. 2017. The effect of MSTN gene in the regulation of skeletal muscle development of pig and its research progress[J]. Genomics and Applied Biology,36(4):1511-1517.]

    崔建勛,杜紅麗,張細權(quán). 2005. 利用DNA池和測序技術(shù)快速篩查SNPs及估算基因頻率[J]. 遺傳學(xué)報,32(4):372-377. [Cui J X,Du H L,Zhang X Q. 2005. Rapidly scree-ning SNPs and estmiating allelic frequencies by DNA pooling and sequencing[J]. Acta Genetica Sinica,32(4):372-377.]

    顧志良,朱大海,李寧,李輝,鄧學(xué)梅,吳常信. 2003. 雞Myostatin基因單核苷酸多態(tài)性與骨骼肌和脂肪生長的關(guān)系[J]. 中國科學(xué)(C輯),33(3):273-280. [Gu Z L,Zhu D H,Li N,Li H,Deng X M,Wu C X. 2003. Relationship between single nucleotide polymorphism of chicken Myostatin gene and skeletal muscle and fat growth[J]. Science in China (Series C),33(3):273-280.]

    李俊英. 2005. 不同優(yōu)質(zhì)雞品種與隱性白羽肉雞雜交組合的屠體及肉品質(zhì)性狀分析[D]. 北京:中國農(nóng)業(yè)大學(xué). [Li J Y. 2005. Analysis of carcass and meat quality traits of crosses of different quality chicken and recessive white plymouth rock[D]. Beijing:China Agricultural University.]

    李洪林. 2019. 貴州地方雞種隱性白羽與慢羽分子標(biāo)記研究與應(yīng)用[D]. 貴陽:貴州大學(xué). [Li H L. 2019. Study and application of recessive white feather and slow feather molecular markers in Guizhou local chicken breeds[D]. Guiyang:Guizhou University.]

    李維,祖盤玉,李洪林,林家棟,張福平. 2017. 赤水烏骨雞肌肉生長抑制素基因多態(tài)性及其與肉質(zhì)性狀的關(guān)聯(lián)分析[J]. 中國畜牧獸醫(yī),44(4):1115-1121. [Li W,Zu P Y,Li H L,Lin J D,Zhang F P. 2017. Polymorphisms of MSTN gene and its association with meat quality traits in Chishui black-bone chicke[J]. China Animal Husbandry & Veterinary Medicine,44(4):1115-1121.]

    李倩. 2013. mRNA二級結(jié)構(gòu)與基因功能和基因必要性的關(guān)系[D]. 楊凌:西北農(nóng)林科技大學(xué). [Li Q. 2013. Correlations among mRNA secondary structrue,gene functions and essentiality[D]. Yangling:Northwest A & F University.]

    劉雪穎. 2008. 藏雞三個基因的SNPs檢測及其與生長性狀的關(guān)聯(lián)分析和金水烏雞資源群體的構(gòu)建[D]. 武漢:華中農(nóng)業(yè)大學(xué). [Liu X Y. 2008. Tibetan chicken SNPs detection of three chicken genes and their association with growth traits and construction of resource population for Jinshui Silky[D]. Wuhan:Huazhong Agricultural University.]

    劉雅正. 2014. 羅曼蛋雞營養(yǎng)需要[J]. 國外畜牧學(xué):豬與禽,(5):47. [Liu Y Z. 2014. Nutritional needs of Luoman Laying hens[J]. Foreign Animal Husbandry: Pigs and Avians,(5):47.]

    毛亮,徐亞歐,麻麗霞. 2009. 藏雞肌肉生長抑制素基因克隆及序列分析[J]. 西南民族大學(xué)學(xué)報(自然科學(xué)版),35(4):772-775. [Mao L,Xu Y O,Ma L X. 2009. Study on cloning of the MSTN gene in the Tibetan gallus[J]. Journal of Southwest University for Nationalities(Natural Science Edition),35(4):772-775.]

    任衛(wèi)合,任穩(wěn)穩(wěn),馬博妍,高成宏,李宏旭,劉麗霞,曹忻,張麗. 2018. 大耳白兔MSTN基因SNPs位點篩查及生物信息學(xué)分析[J]. 浙江農(nóng)業(yè)學(xué)報,30(12):2018-2023. [Ren W H,Ren W W,Ma B Y,Gao C H,Li H X,Liu L X,Cao X,Zhang L. 2018. SNPs screening and bioinforma-tics analysis of MSTN gene in white rabbits[J]. Acta Agriculturae Zhejiangensis,30(12):2018-2023.]

    王芬,王鵬. 2009. 朗德鵝與昌圖豁鵝雜交后代性能試驗[J]. 中國畜禽種業(yè),(10):144-145. [Wang F,Wang P. 2009. Performance test of hybrid offspring of Landes goose and Changtu goose[J]. The Chinese Livestock and Poultry Breeding,(10):144-145.]

    王巧巧,王亞南,唐輝,王金忠,高翔,李顯耀. 2018. 雞肌纖維發(fā)育相關(guān)基因研究進展[J]. 家禽科學(xué),(5):43-48. [Wang Q Q,Wang Y N,Tang H,Wang J Z,Gao X,Li X Y. 2018. Advances in research on genes related to chi-cken muscle fiber development[J]. Poultry Science,(5):43-48.]

    王彥欽,邵勇鋼. 2017. 拜城油雞肌肉生長抑制素基因多態(tài)性與生產(chǎn)性能的相關(guān)性研究[J]. 中國畜牧獸醫(yī),44(4):1054-1060. [Wang Y Q,Shao Y G. 2017. Study on the association relationship between the polymorphism of MSTN gene and production performance in Baicheng fatty chicken[J]. China Animal Husbandry & Veterinary Medicine,44(4):1054-1060.]

    溫彥濤,武子寅,趙振華,李春苗,黎壽豐,王杏龍. 2012. 雞MSTN基因多態(tài)性及其與屠體性狀的關(guān)聯(lián)分析[J]. 中國家禽,34(16):29-32. [Wen Y T,Wu Z Y,Zhao Z H,Li C M,Li S F,Wang X L. 2012. Polymorphisms in exon 1 of myostatin gene and its relationship with carcass traits in chicken[J]. China Poultry,34(16):29-32.]

    翁吉梅,羅錚,徐偉,封竣琪,蔡惠芬,羅衛(wèi)星,陳說. 2018. 黔北麻羊FTO基因SNPs篩選及其生物信息學(xué)分析[J]. 中國畜牧雜志,54(9):39-45. [Weng J M,Luo Z,Xu W,F(xiàn)eng J Q,Cai H F,Luo W X,Chen S. 2018. Screen and bioinformatics analysis of FTO gene SNPs in Qianbei Ma goat[J]. Chinese Journal of Animal Science,54(9):39-45.]

    楊永生,賀建華,鄧惠中,謝紅兵,江碧波,羅佳捷. 2012. 肌肉生長抑制素對動物肌肉、脂肪和骨骼的影響[J]. 動物營養(yǎng)學(xué)報,24(2):220-225. [Yang Y S,He J H,Deng H Z,Xie H B,Jiang B B,Luo J J. 2012. Myostain:Effects on muscle,fat and bone of animals[J]. Chinese Journal of Animal Nutrition,24(2):220-225.]

    張麗,劉麗霞,戴洪偉,陳紅,王瑞,岳炳輝. 2018. 靜寧雞肌肉生長抑制素基因單核苷酸多態(tài)性篩查及生物信息學(xué)分析[J]. 浙江大學(xué)學(xué)報(農(nóng)業(yè)與生命科學(xué)版),44(5):629-637. [Zhang L,Liu L X,Dai H W,Chen H,Wang R,Yue B H. 2018. Single nucleotide polymorphism scree-ning and bioinformatics analysis of myostatin gene in Jingning chicken[J]. Journal of Zhejiang University(Agri-culture & Life Sciences),44 (5):629-637.]

    張冉,林浴霜. 2012. 不同品種雞MSTN基因表達差異的研究[J]. 中國家禽,34(6):61-63. [Zhang R,Lin Y S. 2012. Study on the difference of MSTN gene expression in different breeds of chickens[J]. China Poultry,34(6):61-63.]

    張賢嫻,冉金山,任鵬,李菁菁,尹玲倩,許應(yīng)鋒,陳異,黎志強,劉益平. 2018. 雞肌肉生長抑制素基因單核苷酸多態(tài)性檢測及其與生長性狀的關(guān)聯(lián)性分析[C]//中國畜牧獸醫(yī)學(xué)會. 中國畜牧獸醫(yī)學(xué)會2018年學(xué)術(shù)年會禽病學(xué)分會第十九次學(xué)術(shù)研討會論文集:416. [Zhang X X,Ran J S,Ren P,Li J J,Yin L Q,Xu Y F,Chen Y,Li Z Q,Liu Y P. 2018. Detection of single nucleotide polymorphism of chicken Myostatin gene and its association analysis with growth traits[C]//Chinese Society of Animal Husbandry and Veterinary Society. The 19th Academic Seminar of 2018 Annual Conference Avian Diseases.]

    趙義龍,黃金鳳,張琪智,買爾外提·波拉提,趙金香,矯繼峰. 2017. 肌肉生長抑制素在動物生產(chǎn)中的應(yīng)用[J]. 四川畜牧獸醫(yī),(8):29-31. [Zhao Y L,Huang J F,Zhang Q Z,Mauriti·Bolati,Zhao J X,Jiao J F. 2017. Myostatin application in animal production[J]. Sichuan Animal & Veterinary Medicine,(8):29-31.]

    趙振華,黎壽豐,黃華云,李春苗,王錢保,薛龍崗. 2015. MSTN基因突變對優(yōu)質(zhì)雞肌肉生長的影響[J]. 安徽農(nóng)業(yè)大學(xué)學(xué)報,42(5):733-737. [Zhao Z H,Li S F,Huang H Y,Li C M,Wang Q B,Xue L G. 2015. The myostatin gene(MSTN) and its relationship with muscle fiber traits in chickens[J]. Journal of Anhui Agricultural University,42(5):733-737.]

    鄭海軍,朱榮,葛春蕾,蔡鮮,徐望,陳偉,陸源. 2013. 人白細胞介素-29的生物信息學(xué)分析[J]. 中國生物制品學(xué)雜志,26(2):209-212. [Zheng H J,Zhu R,Ge C L,Cai X,Xu W,Chen W,Lu Y. 2013. Bioinformatics of human interleukin-29[J]. Chinese Journal of Biologicals,26(2):209-212.]

    鄭云. 2007. 鵝Myostatin基因單核苷酸多態(tài)性與生產(chǎn)性能相關(guān)的研究[D]. 揚州:揚州大學(xué). [Zheng Y. 2007. Study on SNP of the Myostatin gene and its association with production performance in geese[D]. Yangzhou:Yangzhou University.]

    朱麗莉,王文濤,傅筑蔭,陳彬. 2011. 興義矮腳雞白羽系(群)的生產(chǎn)性能觀測[J]. 中國家禽,33(8):53-54. [Zhu L L,Wang W T,F(xiàn)u Z Y,Chen B. 2011. Observation on the production performance of the white feather strain of Xingyi bantam chicken(group)[J]. China Poultry,33(8):53-54.]

    祖盤玉,李維,林家棟,李洪林,劉洋,牟騰慧,龍廣麗,簡華峰,張福平. 2019. 赤水烏骨雞TYR基因多態(tài)性及生物信息學(xué)分析[J]. 南方農(nóng)業(yè)學(xué)報,50(12):2806-2811. [Zu P Y,Li W,Lin J D,Li H L,Liu Y,Mou T H,Long G L,Jian H F,Zhang F P. 2019. Polymorphism and bioinformatics analysis of TYR gene in Chishui black-bone chicken[J]. Journal of Southern Agriculture,50(12):2806-2811.]

    左斌,任冰冰,蔣小松,熊應(yīng)龍,王康環(huán),蔣利,劉光偉,王海,徐亞歐. 2015. 瀘寧雞和米易雞MSTN基因多態(tài)性及其與生長性狀的關(guān)聯(lián)性[J]. 貴州農(nóng)業(yè)科學(xué),43(7):26-30. [Zuo B,Ren B B,Jiang X S,Xiong Y L,Wang K H,Jiang L,Liu G W,Wang H,Xu Y O. 2015. Polymorphism of myostatin(MSTN) gene and correlation between MSTN polymorphism and growth traitss in Luning and Miyi chicken[J]. Guizhou Agricultural Sciences,43(7):26-30.]

    Bellinge R H S,Liberles D A,Iaschi S P A,Obrien P A,Tay G K. 2005. Myostatin and its implications on animal breeding:A review[J]. Animal Genetics,36(1):1-6.

    Dushyanth K,Bhattacharya T K,Shukla R,Chatterjee R N,Sitaramamma T,Paswan C,Guru Vishnu P. 2016. Gene expression and polymorphism of myostatin gene and its association with growth traits in chicken[J]. Animal Biotechnology,27(4):269-277.

    Fontanesi L,Scotti E,F(xiàn)rabetti A,F(xiàn)ornasini D,Picconi A,Russo V. 2011. Identification of polymorphisms in the rabbit (Oryctolagus cuniculus) myostatin(MSTN) gene and association analysis with finishing weight in a commercial rabbit population[J]. Animal Genetics,42(3):339. doi:10.1111/j.1365-2052.2010.02163.x.

    Feldman B J,Streeper R S,F(xiàn)arese R V Jr,Yamamoto K R. 2006. Myostatin modulates adipogenesis to generate adipocytes with favorable metabolic effects[J]. Proceedings of the National Academy of Sciences of the United States of America,103(42):15675-15680.

    Grade C V C,Salerno M S,Schubert F R,Dietrich S,Alvares L E. 2009. An evolutionarily conserved Myostatin proximal promoter/enhancer confers basal levels of transcription and spatial specificity in vivo[J]. Development Genes and Evolution,219(9-10):497-508.

    Grobet L,Martin L J,Poncelet D,Pirottin D,Brouwers B,Riquet J,Schoeberlein A,Dunner S,Ménissier F,Massabanda J,F(xiàn)ries R,Hanset R,Georges M. 1997. A deletion in the bovine myostatin gene causes the double-muscled phenotype in cattle[J]. Nature Genetics,17(1):71-74.

    Guo W,F(xiàn)lanagan J,Jasuja J,Kirkland J,Jiang L,Bhasin S. 2008. The effects of myostatin on adipogenic differentiation of human bone marrow-derived mesenchymal stem cells are mediated through cross-communication between Smad3 and Wnt/beta-catenin signaling pathways[J]. The Journal of Biological Chemistry,283(14):9136-9145.

    Kim Y S,Bobbili N K,Lee Y K,Jin H J,Dunn M A. 2006. Production of a polyclonal anti-myostatin antibody and the effects of in ovo administration of the antibody on posthatch broiler growth and muscle mass[J]. Poultry Science,86(6):1196-1205.

    Lee S J. 2004. Regulation of muscle mass by myostatin[J]. Annual Review of Cell and Developmental Biology,20:61-86. doi:10.1146/annurev.cellbio.20.012103.135836.

    McFarland D C,Velleman S G,Pesall J E,Liu C N. 2007. The role of myostatin in chicken(Gallus domesticus) myogenic satellite cell proliferation and differentiation[J]. General and Comparative Endocrinology,151(3):351-357.

    McPherron A C,Lawler A M,Lee S J. 1997. Regulation of skeletal muscle mass in mice by a new TGF-beta superfamily member[J]. Nature,387(6628):83-90.

    Mitrofanova O V,Dementeva N V,Krutikova A A,Yurchenko O P,Vakhrameev A B,Terletskiy V P. 2017. Association of polymorphic variants in MSTN,PRL,and DRD2 genes with intensity of young animal growth in Pushkin breed chickens[J]. Cytology and Genetics,51(3):179-184.

    Xu T S,Gu L H,Zhang X H,Ye B G,Liu X L,Hou S S. 2013. Characterization of myostatin gene(MSTN) of Pekin duck and the association of its polymorphism with breast muscle traits[J]. Genetics and Molecular Research,12(3):3166-3177.

    Zhang G X,Zhang T,Wei Y,Ding F X,Zhang L,Wang J Y. 2015. Functional identification of an exon 1 substitution in the myostatin gene and its expression in breast and leg muscle of the Bian chicken[J]. British Poultry Science,56(6):639-644.

    (責(zé)任編輯 蘭宗寶)

    收稿日期:2019-12-29

    基金項目:貴州省科技計劃項目(黔科合支撐〔2016〕2507號)

    作者簡介:*為通訊作者,張福平(1978-),博士,副教授,主要從事家禽育種研究工作,E-mail:zfu-1010@126.com。龍廣麗(1994-),研究方向為畜牧養(yǎng)殖,E-mail:632816729@qq.com

    猜你喜歡
    多態(tài)性生物信息學(xué)
    APOE基因多態(tài)性與老年動脈粥樣硬化性腦梗死嚴重程度及預(yù)后相關(guān)性分析
    基因多態(tài)性與老年高血壓的研究進展
    淺談醫(yī)學(xué)院校生物信息學(xué)專業(yè)青年教師規(guī)范培訓(xùn)模式的建立
    “PBL+E—learning”教學(xué)模式探索
    移動教學(xué)在生物信息學(xué)課程改革中的應(yīng)用
    中醫(yī)大數(shù)據(jù)下生物信息學(xué)的發(fā)展及教育模式淺析
    生物信息學(xué)課堂危機及對策研究
    西藏砂生槐EST—SSR引物開發(fā)及多態(tài)性檢測
    番鴨A—FABP基因內(nèi)含子多態(tài)性與生產(chǎn)性能的關(guān)聯(lián)分析
    TGF—β1基因多態(tài)性與糖尿病腎病患者易感性關(guān)系的探討
    h日本视频在线播放| 亚洲,欧美,日韩| av女优亚洲男人天堂| 亚洲婷婷狠狠爱综合网| 免费av毛片视频| 亚洲国产精品sss在线观看| 国产av码专区亚洲av| 国产又色又爽无遮挡免| 成人亚洲精品av一区二区| 欧美成人一区二区免费高清观看| 中国国产av一级| 久久婷婷人人爽人人干人人爱| 欧美精品一区二区大全| h日本视频在线播放| 看十八女毛片水多多多| 欧美成人精品欧美一级黄| 一区二区三区免费毛片| 青春草亚洲视频在线观看| a级毛片免费高清观看在线播放| 国产成年人精品一区二区| 亚洲国产高清在线一区二区三| 亚洲18禁久久av| 美女黄网站色视频| 波多野结衣高清无吗| 久久久久免费精品人妻一区二区| a级毛色黄片| 乱人视频在线观看| 激情 狠狠 欧美| 亚洲成人中文字幕在线播放| 免费黄色在线免费观看| 一本一本综合久久| 国产精品嫩草影院av在线观看| 18禁在线播放成人免费| 免费观看人在逋| 亚洲国产精品成人综合色| 少妇的逼好多水| 亚洲乱码一区二区免费版| 欧美一区二区精品小视频在线| av福利片在线观看| www.av在线官网国产| 精品久久久久久久人妻蜜臀av| 特级一级黄色大片| 国产精品久久久久久久久免| 三级毛片av免费| 国产精华一区二区三区| 亚洲成人av在线免费| 好男人在线观看高清免费视频| 国产精品爽爽va在线观看网站| 一边亲一边摸免费视频| 综合色丁香网| 国产高清有码在线观看视频| 婷婷色综合大香蕉| 少妇人妻一区二区三区视频| 午夜福利在线在线| 三级经典国产精品| 天堂中文最新版在线下载 | 国产精品人妻久久久久久| 亚洲国产成人一精品久久久| 国产色爽女视频免费观看| kizo精华| 亚洲在久久综合| 免费av观看视频| 国产精品国产三级国产专区5o | 精品无人区乱码1区二区| 久久精品国产亚洲网站| 中国国产av一级| 国国产精品蜜臀av免费| 我要看日韩黄色一级片| 亚洲最大成人av| 九色成人免费人妻av| 国产精品av视频在线免费观看| 丝袜喷水一区| 亚洲人成网站在线播| 日本黄色片子视频| 能在线免费看毛片的网站| 最近视频中文字幕2019在线8| av在线观看视频网站免费| 欧美精品一区二区大全| 在线观看av片永久免费下载| 成年女人看的毛片在线观看| 在线a可以看的网站| 大香蕉久久网| 精品久久久久久久人妻蜜臀av| 久久这里有精品视频免费| 国产欧美日韩精品一区二区| 欧美变态另类bdsm刘玥| 欧美+日韩+精品| 色综合色国产| 人体艺术视频欧美日本| 日本黄色视频三级网站网址| 我的老师免费观看完整版| 日本-黄色视频高清免费观看| 精品国产一区二区三区久久久樱花 | 亚洲最大成人中文| 日韩一本色道免费dvd| 久久久午夜欧美精品| 国产av码专区亚洲av| 免费看a级黄色片| 亚洲熟妇中文字幕五十中出| 久久久久久大精品| 一个人观看的视频www高清免费观看| 成人无遮挡网站| av在线播放精品| 日本av手机在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 毛片一级片免费看久久久久| 日韩高清综合在线| 国语对白做爰xxxⅹ性视频网站| 欧美成人a在线观看| 男插女下体视频免费在线播放| 国产亚洲精品久久久com| 精品无人区乱码1区二区| 免费在线观看成人毛片| 欧美成人免费av一区二区三区| 成人亚洲精品av一区二区| 久久精品国产亚洲av天美| av在线老鸭窝| 尾随美女入室| 日韩视频在线欧美| 午夜视频国产福利| 亚洲在线自拍视频| 看片在线看免费视频| 亚洲性久久影院| av天堂中文字幕网| 麻豆久久精品国产亚洲av| 1000部很黄的大片| 久久精品人妻少妇| 波多野结衣高清无吗| 少妇人妻一区二区三区视频| 亚洲美女视频黄频| 国产亚洲91精品色在线| 国产视频首页在线观看| 精品免费久久久久久久清纯| 99久久无色码亚洲精品果冻| 人妻少妇偷人精品九色| 99久国产av精品| 精品99又大又爽又粗少妇毛片| 国产精品麻豆人妻色哟哟久久 | 精品国产露脸久久av麻豆 | 女人被狂操c到高潮| 99热这里只有精品一区| 久久久国产成人精品二区| 少妇人妻精品综合一区二区| a级毛色黄片| 欧美极品一区二区三区四区| 日韩欧美国产在线观看| 国产色婷婷99| 免费一级毛片在线播放高清视频| 亚洲在久久综合| 日韩,欧美,国产一区二区三区 | 日本与韩国留学比较| 国产一区二区三区av在线| 国产亚洲5aaaaa淫片| av播播在线观看一区| 久久久久久九九精品二区国产| 国产爱豆传媒在线观看| 91午夜精品亚洲一区二区三区| 久久精品久久久久久噜噜老黄 | 在现免费观看毛片| 成人鲁丝片一二三区免费| 亚洲国产精品久久男人天堂| 又黄又爽又刺激的免费视频.| 可以在线观看毛片的网站| 啦啦啦观看免费观看视频高清| 天堂av国产一区二区熟女人妻| 成人综合一区亚洲| 亚洲欧美成人精品一区二区| 成人高潮视频无遮挡免费网站| 欧美激情国产日韩精品一区| 岛国毛片在线播放| 高清午夜精品一区二区三区| 国产高潮美女av| 久久精品人妻少妇| 免费观看a级毛片全部| 国产精品人妻久久久影院| 91久久精品国产一区二区三区| 一区二区三区四区激情视频| 午夜福利视频1000在线观看| 国产精品嫩草影院av在线观看| 日韩 亚洲 欧美在线| 国产精品伦人一区二区| 亚洲不卡免费看| 成人综合一区亚洲| 久久精品熟女亚洲av麻豆精品 | 男人和女人高潮做爰伦理| av女优亚洲男人天堂| 免费看日本二区| 久久久a久久爽久久v久久| 中文资源天堂在线| 亚洲国产精品合色在线| 午夜福利网站1000一区二区三区| av国产免费在线观看| 爱豆传媒免费全集在线观看| 嫩草影院精品99| 六月丁香七月| 精品人妻熟女av久视频| 少妇熟女欧美另类| 中文字幕亚洲精品专区| 十八禁国产超污无遮挡网站| 日韩,欧美,国产一区二区三区 | 乱人视频在线观看| 亚洲va在线va天堂va国产| 午夜福利视频1000在线观看| 国产精品一二三区在线看| 亚洲国产成人一精品久久久| 精品久久久久久久人妻蜜臀av| 亚洲国产精品合色在线| 22中文网久久字幕| 久久婷婷人人爽人人干人人爱| 噜噜噜噜噜久久久久久91| 亚洲欧美清纯卡通| 国产亚洲一区二区精品| 久久精品久久久久久久性| 天堂网av新在线| 在线a可以看的网站| 99视频精品全部免费 在线| 成人亚洲欧美一区二区av| 菩萨蛮人人尽说江南好唐韦庄 | 久久久亚洲精品成人影院| 国产黄片美女视频| 人人妻人人看人人澡| 亚洲第一区二区三区不卡| 亚洲成人精品中文字幕电影| 亚洲一级一片aⅴ在线观看| 男女国产视频网站| 一级毛片我不卡| 久99久视频精品免费| 99热这里只有精品一区| 国产午夜精品一二区理论片| 国产白丝娇喘喷水9色精品| 欧美日韩一区二区视频在线观看视频在线 | 日韩在线高清观看一区二区三区| 最近中文字幕2019免费版| 日日啪夜夜撸| 亚洲熟妇中文字幕五十中出| 午夜福利网站1000一区二区三区| 亚洲性久久影院| 九九在线视频观看精品| 五月伊人婷婷丁香| 综合色av麻豆| 国产中年淑女户外野战色| 午夜久久久久精精品| 亚洲最大成人中文| 亚洲精品日韩在线中文字幕| 亚洲av.av天堂| 两个人视频免费观看高清| 欧美xxxx性猛交bbbb| 久久人妻av系列| 日本熟妇午夜| 只有这里有精品99| 五月玫瑰六月丁香| 一个人观看的视频www高清免费观看| 欧美xxxx性猛交bbbb| av在线老鸭窝| 成人亚洲欧美一区二区av| 人人妻人人澡欧美一区二区| 一边亲一边摸免费视频| 少妇猛男粗大的猛烈进出视频 | 岛国毛片在线播放| 国产乱人视频| 日韩三级伦理在线观看| 国产成人a∨麻豆精品| 欧美成人免费av一区二区三区| 国产午夜精品一二区理论片| 亚洲中文字幕一区二区三区有码在线看| 久久精品夜夜夜夜夜久久蜜豆| 18禁在线播放成人免费| 好男人视频免费观看在线| 韩国av在线不卡| 精品久久久噜噜| 亚洲成色77777| 97人妻精品一区二区三区麻豆| 国产精品日韩av在线免费观看| 一级毛片我不卡| 老司机福利观看| 少妇的逼水好多| 精品熟女少妇av免费看| 91av网一区二区| 国产视频首页在线观看| 好男人在线观看高清免费视频| 久久精品国产鲁丝片午夜精品| 不卡视频在线观看欧美| 自拍偷自拍亚洲精品老妇| 成年女人永久免费观看视频| 99久久精品国产国产毛片| 亚洲美女搞黄在线观看| 欧美不卡视频在线免费观看| 91aial.com中文字幕在线观看| 国产亚洲午夜精品一区二区久久 | 熟妇人妻久久中文字幕3abv| 国产三级中文精品| av又黄又爽大尺度在线免费看 | 国产精品久久久久久久久免| 国产精品人妻久久久久久| 国内精品美女久久久久久| 久久久久久伊人网av| 91久久精品国产一区二区成人| 少妇熟女欧美另类| 国产免费又黄又爽又色| 精品一区二区三区人妻视频| 26uuu在线亚洲综合色| 免费观看人在逋| 麻豆精品久久久久久蜜桃| 看黄色毛片网站| 国产精品一区二区性色av| 国产成人免费观看mmmm| 美女内射精品一级片tv| 午夜福利成人在线免费观看| 欧美激情国产日韩精品一区| ponron亚洲| a级毛片免费高清观看在线播放| 麻豆成人午夜福利视频| 观看免费一级毛片| 人体艺术视频欧美日本| 国产午夜精品久久久久久一区二区三区| 亚洲国产欧美在线一区| 国模一区二区三区四区视频| 日本三级黄在线观看| 麻豆成人av视频| 久久热精品热| 国产高清视频在线观看网站| 日本欧美国产在线视频| 一级爰片在线观看| 亚洲欧美精品自产自拍| 美女国产视频在线观看| 一区二区三区免费毛片| 国产乱人偷精品视频| 国产黄片美女视频| 久久精品夜夜夜夜夜久久蜜豆| 最近中文字幕高清免费大全6| 三级毛片av免费| 亚洲在久久综合| 在线观看一区二区三区| 亚洲va在线va天堂va国产| 久久久国产成人免费| 日本一本二区三区精品| 全区人妻精品视频| 黑人高潮一二区| 在线免费十八禁| 中文字幕久久专区| 身体一侧抽搐| 九色成人免费人妻av| 国产人妻一区二区三区在| 老师上课跳d突然被开到最大视频| 日本-黄色视频高清免费观看| 久久人人爽人人爽人人片va| 乱系列少妇在线播放| 国产又色又爽无遮挡免| 深夜a级毛片| 久久国内精品自在自线图片| 亚洲,欧美,日韩| 亚洲自偷自拍三级| 亚洲一区高清亚洲精品| 91狼人影院| 日韩中字成人| 六月丁香七月| 又黄又爽又刺激的免费视频.| 国产极品精品免费视频能看的| 身体一侧抽搐| 日韩欧美三级三区| 国内精品宾馆在线| 日韩国内少妇激情av| 男女下面进入的视频免费午夜| 亚洲婷婷狠狠爱综合网| 欧美97在线视频| av播播在线观看一区| 99在线视频只有这里精品首页| 亚洲性久久影院| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 99热这里只有是精品在线观看| 成人亚洲欧美一区二区av| 在线观看美女被高潮喷水网站| 一卡2卡三卡四卡精品乱码亚洲| 国产精品国产三级专区第一集| 国产成人a∨麻豆精品| 亚洲丝袜综合中文字幕| 国产三级在线视频| 天天躁日日操中文字幕| 久久久精品大字幕| 亚洲伊人久久精品综合 | 天天一区二区日本电影三级| 爱豆传媒免费全集在线观看| 蜜桃亚洲精品一区二区三区| 插逼视频在线观看| 国产av不卡久久| 欧美3d第一页| 免费av不卡在线播放| 97在线视频观看| 久久草成人影院| 岛国毛片在线播放| 国产成人精品婷婷| 色视频www国产| 国产精品.久久久| 高清日韩中文字幕在线| 乱系列少妇在线播放| 黄色一级大片看看| 国产精品一及| 国产老妇伦熟女老妇高清| 国产精品不卡视频一区二区| 亚洲精品成人久久久久久| 亚洲精品国产av成人精品| 成人欧美大片| 亚洲精品日韩在线中文字幕| 国产黄色小视频在线观看| 男女那种视频在线观看| 男人舔奶头视频| 欧美激情在线99| 高清在线视频一区二区三区 | 国产一区二区三区av在线| 国产精品日韩av在线免费观看| 男人狂女人下面高潮的视频| 亚洲人成网站在线观看播放| 2022亚洲国产成人精品| 亚洲最大成人av| 中文乱码字字幕精品一区二区三区 | 免费观看性生交大片5| 直男gayav资源| 中文亚洲av片在线观看爽| a级毛片免费高清观看在线播放| 成人亚洲精品av一区二区| 国产午夜福利久久久久久| 欧美一区二区亚洲| 午夜久久久久精精品| 日本午夜av视频| 国产成人aa在线观看| 日韩国内少妇激情av| 爱豆传媒免费全集在线观看| 国产单亲对白刺激| 亚洲av日韩在线播放| 欧美成人一区二区免费高清观看| 日韩在线高清观看一区二区三区| 97超视频在线观看视频| 女的被弄到高潮叫床怎么办| 午夜福利在线在线| 一本一本综合久久| 九草在线视频观看| 最近视频中文字幕2019在线8| 夜夜爽夜夜爽视频| 国产成人免费观看mmmm| 免费看av在线观看网站| 寂寞人妻少妇视频99o| av视频在线观看入口| 欧美人与善性xxx| 国产精品精品国产色婷婷| 日本免费一区二区三区高清不卡| 偷拍熟女少妇极品色| 九草在线视频观看| 亚洲va在线va天堂va国产| 纵有疾风起免费观看全集完整版 | 久久久久网色| 日本与韩国留学比较| 欧美极品一区二区三区四区| 免费不卡的大黄色大毛片视频在线观看 | 国产色婷婷99| 亚洲av电影在线观看一区二区三区 | 欧美丝袜亚洲另类| 人人妻人人澡欧美一区二区| 99热全是精品| 日韩av不卡免费在线播放| 国产淫片久久久久久久久| 熟女电影av网| 校园人妻丝袜中文字幕| 天堂av国产一区二区熟女人妻| av免费在线看不卡| 一级av片app| 国产三级在线视频| 国产中年淑女户外野战色| 边亲边吃奶的免费视频| 好男人视频免费观看在线| 国内精品一区二区在线观看| 我要搜黄色片| 亚洲国产日韩欧美精品在线观看| 国产精品一区二区三区四区久久| 亚洲天堂国产精品一区在线| 日本wwww免费看| 欧美成人精品欧美一级黄| 国产一区二区三区av在线| 建设人人有责人人尽责人人享有的 | 91精品一卡2卡3卡4卡| 长腿黑丝高跟| 欧美97在线视频| 亚洲18禁久久av| 欧美一级a爱片免费观看看| 亚洲av免费在线观看| 男女那种视频在线观看| 97超视频在线观看视频| 真实男女啪啪啪动态图| 最近最新中文字幕免费大全7| 日日干狠狠操夜夜爽| 欧美日韩国产亚洲二区| 一区二区三区免费毛片| 亚洲激情五月婷婷啪啪| 免费黄网站久久成人精品| 亚洲熟妇中文字幕五十中出| 国产精品99久久久久久久久| 男人的好看免费观看在线视频| ponron亚洲| 国产精品蜜桃在线观看| 亚洲熟妇中文字幕五十中出| 亚洲在久久综合| 久久久久九九精品影院| 国产亚洲5aaaaa淫片| av线在线观看网站| 色综合站精品国产| 一个人看的www免费观看视频| 国产真实乱freesex| 免费观看a级毛片全部| 丝袜喷水一区| 欧美一区二区亚洲| 高清午夜精品一区二区三区| 建设人人有责人人尽责人人享有的 | 卡戴珊不雅视频在线播放| 成年女人看的毛片在线观看| 干丝袜人妻中文字幕| 免费电影在线观看免费观看| 久久久久免费精品人妻一区二区| 日韩高清综合在线| 国产av在哪里看| 精品久久久久久久久久久久久| 欧美性猛交黑人性爽| 大香蕉久久网| 精品人妻熟女av久视频| 亚洲欧洲国产日韩| 午夜视频国产福利| 亚洲五月天丁香| 亚洲激情五月婷婷啪啪| av线在线观看网站| 日韩强制内射视频| 小说图片视频综合网站| 国产精品国产三级国产专区5o | 少妇的逼水好多| 嫩草影院入口| 日韩欧美 国产精品| 亚洲av不卡在线观看| 婷婷六月久久综合丁香| 亚洲国产精品国产精品| 99热网站在线观看| 亚洲人成网站高清观看| 国产精品,欧美在线| 国产日韩欧美在线精品| 亚洲精品色激情综合| 天堂影院成人在线观看| 亚洲18禁久久av| 成人午夜精彩视频在线观看| 舔av片在线| 亚洲欧美日韩东京热| 一区二区三区四区激情视频| .国产精品久久| 禁无遮挡网站| 国产精品一区二区性色av| 国产人妻一区二区三区在| 大话2 男鬼变身卡| 美女cb高潮喷水在线观看| 国产av码专区亚洲av| 一本久久精品| 欧美日韩综合久久久久久| 性色avwww在线观看| 插逼视频在线观看| 最后的刺客免费高清国语| 99视频精品全部免费 在线| 国内精品宾馆在线| 亚洲精品aⅴ在线观看| 国产精品嫩草影院av在线观看| 日本-黄色视频高清免费观看| videossex国产| 国产精品.久久久| 欧美成人一区二区免费高清观看| 夫妻性生交免费视频一级片| 国产高清有码在线观看视频| 搞女人的毛片| 久久这里只有精品中国| 欧美成人午夜免费资源| 黄片wwwwww| 99热这里只有精品一区| 久久99蜜桃精品久久| 又爽又黄无遮挡网站| 亚洲精品日韩av片在线观看| 啦啦啦啦在线视频资源| 91aial.com中文字幕在线观看| 免费看日本二区| 男人舔奶头视频| 在线免费观看不下载黄p国产| 精品人妻偷拍中文字幕| av女优亚洲男人天堂| 日韩av不卡免费在线播放| 精品无人区乱码1区二区| 精品人妻偷拍中文字幕| 熟女人妻精品中文字幕| 亚洲中文字幕一区二区三区有码在线看| 大香蕉97超碰在线| 中文亚洲av片在线观看爽| 国产亚洲精品久久久com| 久久精品91蜜桃| 午夜精品在线福利| ponron亚洲| 一级黄片播放器| 亚洲成人精品中文字幕电影| 男人狂女人下面高潮的视频| 亚洲自拍偷在线| 成人二区视频| 国产精品无大码| 如何舔出高潮| 久久精品国产自在天天线| 精品人妻一区二区三区麻豆| 日韩大片免费观看网站 | 一个人看的www免费观看视频| 18禁裸乳无遮挡免费网站照片| 大话2 男鬼变身卡| 三级毛片av免费| eeuss影院久久| 国产高清三级在线| 国产亚洲最大av| 成人高潮视频无遮挡免费网站| 美女内射精品一级片tv| 永久免费av网站大全|