• <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)系的探討
    黄色配什么色好看| 日韩不卡一区二区三区视频在线| 香蕉精品网在线| 久久久精品免费免费高清| 亚洲精品日韩在线中文字幕| 日日啪夜夜爽| 国产成人精品婷婷| freevideosex欧美| 亚洲在久久综合| 热re99久久精品国产66热6| 日本爱情动作片www.在线观看| 日日摸夜夜添夜夜添av毛片| 王馨瑶露胸无遮挡在线观看| 久久精品国产自在天天线| 欧美老熟妇乱子伦牲交| 天堂俺去俺来也www色官网| 日韩强制内射视频| 成人毛片a级毛片在线播放| 午夜影院在线不卡| 熟女人妻精品中文字幕| 亚洲av男天堂| av免费观看日本| 欧美精品国产亚洲| 夜夜看夜夜爽夜夜摸| 亚洲国产精品999| 久久久精品免费免费高清| 又黄又爽又刺激的免费视频.| 丝袜在线中文字幕| 日本黄色片子视频| 在线 av 中文字幕| 国语对白做爰xxxⅹ性视频网站| 免费观看无遮挡的男女| 亚洲情色 制服丝袜| 天美传媒精品一区二区| 99九九线精品视频在线观看视频| 国产成人91sexporn| 亚洲国产av新网站| 日韩制服骚丝袜av| 成人国产麻豆网| 久久久久精品久久久久真实原创| 久久99一区二区三区| 免费av不卡在线播放| av有码第一页| 女性生殖器流出的白浆| 日韩,欧美,国产一区二区三区| 九草在线视频观看| 欧美区成人在线视频| 午夜福利在线观看免费完整高清在| 肉色欧美久久久久久久蜜桃| 国产av码专区亚洲av| 国产熟女欧美一区二区| 插逼视频在线观看| 国产成人精品久久久久久| 国产极品天堂在线| 日本与韩国留学比较| 免费人妻精品一区二区三区视频| 亚洲真实伦在线观看| 精华霜和精华液先用哪个| 国产极品粉嫩免费观看在线 | 色婷婷久久久亚洲欧美| 成年人午夜在线观看视频| 久久久久精品性色| 国产黄色视频一区二区在线观看| 在线免费观看不下载黄p国产| 国内精品宾馆在线| 国产片特级美女逼逼视频| 亚洲欧美一区二区三区国产| 最近2019中文字幕mv第一页| 欧美丝袜亚洲另类| 亚洲精品国产av蜜桃| 亚洲人成网站在线播| 啦啦啦啦在线视频资源| 久热久热在线精品观看| 韩国av在线不卡| 欧美丝袜亚洲另类| 欧美成人午夜免费资源| 免费黄频网站在线观看国产| 亚洲第一av免费看| 国产又色又爽无遮挡免| 看免费成人av毛片| 成人亚洲欧美一区二区av| 亚洲电影在线观看av| 久久热精品热| 深夜a级毛片| 欧美日韩精品成人综合77777| 69精品国产乱码久久久| 永久网站在线| 亚洲av中文av极速乱| 国产精品偷伦视频观看了| 黄片无遮挡物在线观看| 久久国产乱子免费精品| 亚洲精品日韩av片在线观看| 少妇人妻久久综合中文| 多毛熟女@视频| 看免费成人av毛片| 亚洲精品一区蜜桃| 丰满饥渴人妻一区二区三| 99久久综合免费| 一二三四中文在线观看免费高清| 亚洲欧洲日产国产| 精品国产乱码久久久久久小说| 丰满乱子伦码专区| 日韩伦理黄色片| 精品久久久久久久久亚洲| 日韩中字成人| 免费人妻精品一区二区三区视频| 日本av免费视频播放| 亚洲精品自拍成人| 99视频精品全部免费 在线| 久久久久久久久久人人人人人人| 插逼视频在线观看| 美女大奶头黄色视频| 日韩av在线免费看完整版不卡| 午夜日本视频在线| 精品国产一区二区久久| 大码成人一级视频| 777米奇影视久久| 国产国拍精品亚洲av在线观看| 色哟哟·www| 国产黄片视频在线免费观看| 91久久精品国产一区二区成人| 丰满少妇做爰视频| 麻豆成人午夜福利视频| 亚洲人与动物交配视频| 国产亚洲一区二区精品| av女优亚洲男人天堂| 国产精品一区www在线观看| 老司机亚洲免费影院| 亚洲人与动物交配视频| 国产中年淑女户外野战色| 亚洲精品国产av蜜桃| 丝瓜视频免费看黄片| 久久精品国产a三级三级三级| 又爽又黄a免费视频| 91成人精品电影| 夜夜爽夜夜爽视频| 亚洲三级黄色毛片| 人妻 亚洲 视频| 成年人免费黄色播放视频 | 日韩人妻高清精品专区| 亚洲激情五月婷婷啪啪| 亚洲精品视频女| 国产午夜精品久久久久久一区二区三区| 亚洲美女视频黄频| 国产欧美另类精品又又久久亚洲欧美| 国产极品天堂在线| 精华霜和精华液先用哪个| 精品久久久久久电影网| 国产亚洲欧美精品永久| 在线观看国产h片| videossex国产| 国产精品久久久久久精品电影小说| 丰满饥渴人妻一区二区三| 在线观看人妻少妇| 久久毛片免费看一区二区三区| 国产在线免费精品| 国产永久视频网站| 免费观看a级毛片全部| 一本色道久久久久久精品综合| 久久人人爽av亚洲精品天堂| 免费观看无遮挡的男女| 下体分泌物呈黄色| 亚洲国产成人一精品久久久| 国产色爽女视频免费观看| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品日韩在线中文字幕| av免费观看日本| 亚洲精品色激情综合| 日本猛色少妇xxxxx猛交久久| 嫩草影院新地址| 美女cb高潮喷水在线观看| 韩国高清视频一区二区三区| 在线观看一区二区三区激情| 久久久精品免费免费高清| 全区人妻精品视频| 国产精品成人在线| 日本wwww免费看| 国产高清国产精品国产三级| 国模一区二区三区四区视频| 三级国产精品片| 又爽又黄a免费视频| 在线精品无人区一区二区三| 精品国产一区二区三区久久久樱花| 国产永久视频网站| 午夜日本视频在线| 一区在线观看完整版| 欧美日韩精品成人综合77777| 一级,二级,三级黄色视频| 一区二区av电影网| 成人漫画全彩无遮挡| 国内精品宾馆在线| 亚洲国产成人一精品久久久| 国产亚洲av片在线观看秒播厂| a级一级毛片免费在线观看| 国产av码专区亚洲av| 欧美日韩精品成人综合77777| 国产一区二区在线观看日韩| 久久国产精品大桥未久av | 女人精品久久久久毛片| 亚洲美女视频黄频| 亚洲av电影在线观看一区二区三区| 中国国产av一级| 丝袜在线中文字幕| 日本91视频免费播放| 97超碰精品成人国产| 亚洲av二区三区四区| 日韩人妻高清精品专区| 欧美精品一区二区免费开放| 又粗又硬又长又爽又黄的视频| 18禁在线无遮挡免费观看视频| 色视频www国产| 美女中出高潮动态图| 久久久国产精品麻豆| 婷婷色麻豆天堂久久| 国产成人91sexporn| 丰满迷人的少妇在线观看| 一级毛片我不卡| 国产亚洲一区二区精品| 精品久久久噜噜| 久久女婷五月综合色啪小说| av在线观看视频网站免费| xxx大片免费视频| 国产成人a∨麻豆精品| 热re99久久精品国产66热6| 日韩电影二区| 九九在线视频观看精品| 国产精品国产三级专区第一集| 99久久精品热视频| 国产日韩一区二区三区精品不卡 | 在线观看免费日韩欧美大片 | 深夜a级毛片| 国产男人的电影天堂91| av国产久精品久网站免费入址| 99热6这里只有精品| 高清在线视频一区二区三区| 亚洲精华国产精华液的使用体验| 国产免费一级a男人的天堂| 久久精品久久久久久噜噜老黄| 在线观看三级黄色| 校园人妻丝袜中文字幕| 五月天丁香电影| 久久国产精品男人的天堂亚洲 | 国产欧美日韩综合在线一区二区 | 欧美 亚洲 国产 日韩一| 高清毛片免费看| 黄色一级大片看看| 最新的欧美精品一区二区| 午夜久久久在线观看| 一级毛片我不卡| 亚洲av不卡在线观看| 少妇精品久久久久久久| 国产成人精品福利久久| 日产精品乱码卡一卡2卡三| 多毛熟女@视频| 亚洲人与动物交配视频| 成人美女网站在线观看视频| 黑丝袜美女国产一区| 欧美bdsm另类| 日本av手机在线免费观看| 国产日韩欧美视频二区| 日韩在线高清观看一区二区三区| 亚洲精品自拍成人| 免费av中文字幕在线| 亚洲精品日本国产第一区| 老熟女久久久| 久久鲁丝午夜福利片| 国产亚洲5aaaaa淫片| 男女边吃奶边做爰视频| 国产精品嫩草影院av在线观看| 亚洲人成网站在线播| 高清视频免费观看一区二区| 婷婷色麻豆天堂久久| 伦理电影免费视频| 乱系列少妇在线播放| 在线观看人妻少妇| 人体艺术视频欧美日本| 七月丁香在线播放| 大片免费播放器 马上看| 久热久热在线精品观看| 九草在线视频观看| av福利片在线| 在线观看美女被高潮喷水网站| 亚洲精品成人av观看孕妇| h视频一区二区三区| 精品人妻偷拍中文字幕| 18禁裸乳无遮挡动漫免费视频| 精品少妇内射三级| 蜜桃久久精品国产亚洲av| 精品人妻熟女毛片av久久网站| 日本av免费视频播放| 国产精品女同一区二区软件| 国产乱来视频区| 精品一品国产午夜福利视频| 亚洲精品一二三| 男的添女的下面高潮视频| 看非洲黑人一级黄片| 日韩中文字幕视频在线看片| 一本大道久久a久久精品| 精品亚洲成a人片在线观看| 精品酒店卫生间| 少妇被粗大猛烈的视频| 亚洲色图综合在线观看| 中文字幕制服av| 成人免费观看视频高清| 日韩制服骚丝袜av| 少妇猛男粗大的猛烈进出视频| 97在线人人人人妻| 色网站视频免费| 丰满人妻一区二区三区视频av| 五月天丁香电影| 99久久精品一区二区三区| 国产精品一区二区性色av| 日韩三级伦理在线观看| 91精品国产国语对白视频| 又爽又黄a免费视频| 亚洲欧美精品自产自拍| 亚洲国产欧美在线一区| 简卡轻食公司| 亚洲人成网站在线观看播放| av在线观看视频网站免费| 91aial.com中文字幕在线观看| 欧美少妇被猛烈插入视频| 国产精品一二三区在线看| 人妻制服诱惑在线中文字幕| 国产成人91sexporn| 亚洲va在线va天堂va国产| 欧美成人午夜免费资源| 国产熟女欧美一区二区| 精品久久久久久久久av| 亚洲精品成人av观看孕妇| 不卡视频在线观看欧美| 精品久久久久久电影网| 精品久久久精品久久久| 国产色爽女视频免费观看| 91午夜精品亚洲一区二区三区| 国国产精品蜜臀av免费| 三级国产精品片| 久久精品国产亚洲av涩爱| av免费观看日本| 国语对白做爰xxxⅹ性视频网站| 青春草视频在线免费观看| 国产亚洲5aaaaa淫片| 国精品久久久久久国模美| 高清在线视频一区二区三区| 国产女主播在线喷水免费视频网站| 国产男女内射视频| 国产中年淑女户外野战色| 国产在线视频一区二区| 18禁裸乳无遮挡动漫免费视频| 亚洲欧美一区二区三区国产| 99re6热这里在线精品视频| 少妇丰满av| 午夜福利影视在线免费观看| 在线观看免费视频网站a站| 久久99精品国语久久久| 丰满迷人的少妇在线观看| 欧美成人午夜免费资源| 精品熟女少妇av免费看| 国产精品三级大全| 久热久热在线精品观看| 国产av码专区亚洲av| 一级毛片电影观看| 自线自在国产av| a级毛片在线看网站| 国产高清有码在线观看视频| 久久精品国产亚洲av天美| 在线看a的网站| 人人妻人人添人人爽欧美一区卜| 亚洲av在线观看美女高潮| 男女免费视频国产| 免费大片18禁| 精品人妻一区二区三区麻豆| 天堂中文最新版在线下载| 国产精品一区二区在线不卡| 亚洲国产成人一精品久久久| 插逼视频在线观看| 搡老乐熟女国产| 狠狠精品人妻久久久久久综合| 欧美区成人在线视频| 多毛熟女@视频| 久久久久久人妻| 成年美女黄网站色视频大全免费 | 精品久久国产蜜桃| 男人狂女人下面高潮的视频| 色网站视频免费| 欧美性感艳星| 一级爰片在线观看| 亚洲精品日韩av片在线观看| 久久国内精品自在自线图片| 国产在线男女| 最新中文字幕久久久久| 欧美 亚洲 国产 日韩一| 人人妻人人爽人人添夜夜欢视频 | 精品久久久精品久久久| 亚洲无线观看免费| www.色视频.com| 韩国av在线不卡| 精品国产露脸久久av麻豆| 激情五月婷婷亚洲| 欧美成人午夜免费资源| 日韩av不卡免费在线播放| 国产精品久久久久久av不卡| 亚洲精华国产精华液的使用体验| 亚洲av成人精品一二三区| av女优亚洲男人天堂| 国产精品一区二区三区四区免费观看| 少妇的逼好多水| 亚洲国产成人一精品久久久| 男的添女的下面高潮视频| 美女福利国产在线| 秋霞伦理黄片| 菩萨蛮人人尽说江南好唐韦庄| 国产中年淑女户外野战色| .国产精品久久| 国产精品欧美亚洲77777| 99热这里只有精品一区| 国产免费一级a男人的天堂| 高清不卡的av网站| 久久婷婷青草| 又大又黄又爽视频免费| 97精品久久久久久久久久精品| 看十八女毛片水多多多| 在线观看国产h片| 欧美一级a爱片免费观看看| 亚洲欧美成人综合另类久久久| 在线免费观看不下载黄p国产| 久久久久久久久久人人人人人人| 青春草亚洲视频在线观看| 亚洲av在线观看美女高潮| 久久免费观看电影| 成人无遮挡网站| 亚洲国产精品999| 自拍欧美九色日韩亚洲蝌蚪91 | 深夜a级毛片| 热99国产精品久久久久久7| 看免费成人av毛片| √禁漫天堂资源中文www| 十八禁高潮呻吟视频 | 99热6这里只有精品| 成人综合一区亚洲| 日产精品乱码卡一卡2卡三| 女性被躁到高潮视频| 免费人成在线观看视频色| 久久久久精品久久久久真实原创| 熟女人妻精品中文字幕| 国产高清不卡午夜福利| 国国产精品蜜臀av免费| 亚洲精品一区蜜桃| 69精品国产乱码久久久| 五月天丁香电影| av天堂中文字幕网| 欧美日韩视频高清一区二区三区二| 日本欧美视频一区| 日本色播在线视频| 国产老妇伦熟女老妇高清| 国产精品一区二区在线观看99| 9色porny在线观看| 青春草视频在线免费观看| 女人精品久久久久毛片| 国产91av在线免费观看| 日韩一本色道免费dvd| 精品久久久精品久久久| 99国产精品免费福利视频| 女人精品久久久久毛片| 综合色丁香网| 成人午夜精彩视频在线观看| 亚洲一级一片aⅴ在线观看| 18禁动态无遮挡网站| 亚洲精品视频女| 2021少妇久久久久久久久久久| 久久久久久久久久久免费av| 我要看日韩黄色一级片| 秋霞伦理黄片| av又黄又爽大尺度在线免费看| 欧美成人午夜免费资源| 黄片无遮挡物在线观看| 麻豆成人午夜福利视频| 肉色欧美久久久久久久蜜桃| 亚洲,一卡二卡三卡| 91精品一卡2卡3卡4卡| 大片电影免费在线观看免费| 亚洲精品一区蜜桃| 亚洲精品日韩av片在线观看| 99热全是精品| 亚洲国产精品一区二区三区在线| 18禁在线无遮挡免费观看视频| 亚洲久久久国产精品| 一级毛片电影观看| 男人添女人高潮全过程视频| 成人国产av品久久久| 老熟女久久久| 日韩中字成人| 另类精品久久| 亚洲精品,欧美精品| 国产片特级美女逼逼视频| 你懂的网址亚洲精品在线观看| 亚洲精品久久午夜乱码| 赤兔流量卡办理| 2018国产大陆天天弄谢| 日本91视频免费播放| 国内揄拍国产精品人妻在线| 热re99久久精品国产66热6| 日韩av不卡免费在线播放| 亚洲欧美中文字幕日韩二区| 欧美日韩一区二区视频在线观看视频在线| av在线播放精品| 内射极品少妇av片p| 免费av中文字幕在线| av专区在线播放| 男女边吃奶边做爰视频| 精品一区二区三卡| 国产精品国产三级国产av玫瑰| 国产深夜福利视频在线观看| 99热国产这里只有精品6| 性色av一级| 久久99蜜桃精品久久| 国产中年淑女户外野战色| 中文天堂在线官网| 男人爽女人下面视频在线观看| 精品国产一区二区三区久久久樱花| 日日摸夜夜添夜夜爱| 亚洲国产日韩一区二区| 免费观看性生交大片5| 中文资源天堂在线| 久久久精品94久久精品| 久久婷婷青草| 日韩av免费高清视频| 日韩欧美一区视频在线观看 | 精品人妻一区二区三区麻豆| 少妇人妻久久综合中文| 久久久久久人妻| 看非洲黑人一级黄片| 特大巨黑吊av在线直播| 69精品国产乱码久久久| 久热久热在线精品观看| 日本色播在线视频| 丝袜在线中文字幕| 日本av免费视频播放| 国产精品一二三区在线看| 日本与韩国留学比较| 妹子高潮喷水视频| 国产精品成人在线| 国产69精品久久久久777片| 国产精品无大码| 妹子高潮喷水视频| 婷婷色综合大香蕉| 超碰97精品在线观看| 日产精品乱码卡一卡2卡三| 国产精品国产三级国产专区5o| freevideosex欧美| 秋霞在线观看毛片| 自拍偷自拍亚洲精品老妇| 午夜激情福利司机影院| 久久久午夜欧美精品| 久久精品久久久久久久性| 国产一区二区三区av在线| 国产亚洲一区二区精品| 国产精品成人在线| 晚上一个人看的免费电影| 中文欧美无线码| 国产精品蜜桃在线观看| 久久精品久久久久久噜噜老黄| 晚上一个人看的免费电影| 最黄视频免费看| 丝瓜视频免费看黄片| 桃花免费在线播放| 丁香六月天网| xxx大片免费视频| 99热这里只有是精品50| av又黄又爽大尺度在线免费看| av专区在线播放| 青春草亚洲视频在线观看| 日韩在线高清观看一区二区三区| 久久午夜综合久久蜜桃| 十八禁网站网址无遮挡 | 国产在线免费精品| 大片免费播放器 马上看| 美女脱内裤让男人舔精品视频| 老司机影院成人| a级毛片在线看网站| 99热这里只有是精品在线观看| 精华霜和精华液先用哪个| 亚洲精品乱码久久久v下载方式| 亚洲伊人久久精品综合| 又爽又黄a免费视频| 国产精品国产三级国产av玫瑰| 黄色欧美视频在线观看| 有码 亚洲区| 日韩一区二区视频免费看| av天堂久久9| 久久鲁丝午夜福利片| 成人无遮挡网站| 制服丝袜香蕉在线| 国产黄色视频一区二区在线观看| 22中文网久久字幕| 99热这里只有是精品在线观看| 蜜臀久久99精品久久宅男| 国产精品蜜桃在线观看| 日本vs欧美在线观看视频 | 97超视频在线观看视频| 国产色婷婷99| 久久午夜综合久久蜜桃| 黄色一级大片看看| 一区二区三区乱码不卡18| 男女免费视频国产| 亚洲精品第二区| 婷婷色综合大香蕉| 91成人精品电影| 国产亚洲午夜精品一区二区久久| 蜜桃在线观看..| 天美传媒精品一区二区| 岛国毛片在线播放| 国产精品久久久久成人av| 亚洲精品视频女| 久久久久久久大尺度免费视频|