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

    綿羊Y染色體特異性引物及SNPs的篩選

    2018-08-17 02:08:00曹學濤裴生偉張晉李發(fā)弟李剛李萬宏樂祥鵬
    中國農業(yè)科學 2018年15期
    關鍵詞:薩福克巖羊綿羊

    曹學濤,裴生偉,張晉,李發(fā)弟,2,李剛,李萬宏,樂祥鵬

    ?

    綿羊Y染色體特異性引物及SNPs的篩選

    曹學濤1,裴生偉1,張晉3,李發(fā)弟1,2,李剛3,李萬宏1,樂祥鵬1

    (1蘭州大學草地農業(yè)生態(tài)系統(tǒng)國家重點實驗室/蘭州大學農業(yè)農村部草牧業(yè)創(chuàng)新重點實驗室/蘭州大學草地農業(yè)科技學院,蘭州 730020;2甘肅省肉羊繁育生物技術工程實驗室,甘肅民勤 733300;3甘肅省家畜繁殖中心,甘肅武威 733000)

    【目的】哺乳動物Y染色體雄性特異區(qū)在減數分裂過程中不與X染色體發(fā)生重組,遵循嚴格的父系遺傳,是研究父系遺傳多樣性的重要遺傳資源;此外絕大多數Y染色體基因主要或者特異性的在睪丸組織中表達,并在精子發(fā)生和雄性繁殖力方面扮演著至關重要的作用。由于Y染色體測序極其困難,造成很多物種Y染色體序列很少。因此本文基于目前現有??苿游颵染色體引物信息,旨在鑒定綿羊Y染色體特異性引物,比較綿羊Y染色體基因片段與巖羊、牛、山羊和牦牛Y染色體的差異,同時篩選不同綿羊品種在Y染色體基因片段內的SNPs,將找到的Y-SNPs和綿羊的睪丸大小進行相關性分析,為鑒定綿羊Y染色體基因片段奠定基礎,并為今后綿羊Y染色體單倍型的構建、綿羊胚胎性別早期鑒定和公綿羊繁殖力相關分子標記的篩選提供科學依據?!痉椒ā扛鶕壳拔墨I公布的29對牛科動物Y染色體特異性引物序列,以母綿羊和水分別作陰性和空白對照,以公綿羊DNA為模板驗證引物的雄性特異性;確定綿羊Y染色體特異引物后,利用DNA混合池測序結合限制性長度片段多態(tài)性等技術在薩??搜颍╪=146)、白薩福克羊(n=91)、東弗里升羊(n=6)、特克塞爾羊(n=72)、南非肉用美利奴羊(n=17)、杜泊羊(n=32)、湖羊(n=55)、藏羊(n=34)、灘羊(n=43)和巖羊(n=14)公羊群體中進行Y-SNPs掃描。用Chromas和DNASTAR等軟件分析混合池測序的結果,利用DNAman軟件將綿羊Y染色體基因片段與牦牛、山羊、黃牛和巖羊進行同源性分析。同時,測量薩??搜颉姿_??搜颉|弗里升羊、特克塞爾羊、南非肉用美利奴羊、杜泊羊周歲的陰囊圍,利用SPSS 19.0軟件分析SNPs位點與公羊陰囊圍的相關性?!窘Y果】在分析的29對引物中,6對引物為綿羊Y染色體特異性引物,分別能擴增出3、6、、、11和6片段。17對引物未能出現擴增條帶,6對引物在母羊DNA中出現擴增條帶。通過比對發(fā)現綿羊6個基因片段與巖羊、牛、山羊、牦牛的同源性在81.51%—98.84%。此外,首次在薩??撕桶姿_福克羊群體的11片段中鑒定得到一個G>A的突變,通過I酶切分析發(fā)現在薩??搜?、白薩??搜蛑杏?種基因型(AA、GG),在其他7個綿羊品種中只有GG基因型。在白薩??搜蛉后w中,GG基因型的頻率為0.747,AA基因型的頻率為0.253;在薩??搜蛉后w中GG基因型的頻率為0.986,AA基因型的頻率為0.014;說明在薩??搜蚝桶姿_??搜蛉后w中,GG基因型為優(yōu)勢基因型。關聯分析顯示在白薩??搜蛉后w中,GG基因型個體周歲的陰囊圍顯著的高于AA基因型(=0.029)?!窘Y論】鑒定得到6個綿羊Y染色體基因片段,它們與牛、山羊、牦牛和巖羊具有較高的同源性,表明Y染色體基因在進化過程中具有一定的保守性。首次在薩福克和白薩??搜蛉后w11片段中找到一個Y-SNP(G>A),其與白薩??搜蛑軞q的睪丸大小緊密相關。

    Y染色體;綿羊;單核苷酸多態(tài)性

    0 引言

    【研究意義】Y染色體基因在研究父系遺傳多樣性、維持雄性特性,精子發(fā)生和雄性繁殖力起著關鍵作用,Y染色體分子標記已經運用于綿羊父系起源進化研究及父系單倍型構建[1-3]。然而Y染色體測序困難,造成綿羊Y染色體基因序列極度貧乏。因此本研究基于前期??苿游颵染色體研究成果,進行綿羊Y染色體基因片段鑒定和Y-SNPs的篩選,以期為優(yōu)秀種公羊的早期選擇、胚胎性別的早期鑒定和綿羊品種的父系遺傳多樣性的保護提供理論依據。【前人研究進展】哺乳動物的性染色體,X和Y染色體,是由一對原始的常染色體進化而來[4]。Y染色體是雄性動物特有的性染色體,由兩部分組成:擬常染色體區(qū)域和Y染色體雄性特異區(qū)(male specific region of Y chromosome, MSY)組成。MSY區(qū)在減數分裂過程中不與X染色體發(fā)生重組交換,遵循嚴格的父系遺傳。Y染色體結構非常復雜,包含有大量的重復序列和回文結構,使得其序列的組裝極為困難[5]。到目前為止,綿羊Y染色體只有、等少數基因的部分序列已公布[6-8]。MEADOWS等利用23對人、牛、綿羊Y染色體特異性引物擴得到綿羊、、、、等Y染色體基因的部分片段[9]。FENG等擴增得到巖羊ZFY基因639 bp的片段[7]。劉帥兵克隆得到ZFY基因整個mRNA全長[10]。熊勇等克隆得到藏系綿羊ZFY基因447 bp的片段[11]。蔣利等利用山羊SRY基因序列設計引物,擴增得到阿勒泰羊、藏綿羊SRY基因整個編碼區(qū)的雄性特異性片段,研究表明SRY基因在不同綿羊品種間具有高度的同源性[12]?!颈狙芯壳腥朦c】目前只有人類、猩猩、獼猴、小鼠和牛等少數物種的Y染色體完成了測序[3,5,13-15]。關于綿羊Y染色體,少數研究者只擴增出基因的部分片段,未見克隆得到完整的Y染色體基因序列?!緮M解決的關鍵問題】本研究基于牛科動物現有Y染色體特異引物鑒定綿羊Y染色體單拷貝基因,掃描Y染色體基因在不同綿羊品種中的SNPs,以期找到Y染色體分子標記,以期為篩選不同羊品種Y染色體SNPs提供基礎,也為綿羊早期胚胎性別的鑒定提供了分子標記,為進一步研究綿羊Y染色體基因序列奠定基礎。

    1 材料與方法

    試驗于2016年10月至2017年5月在蘭州大學草地農業(yè)生態(tài)系統(tǒng)國家重點實驗室完成。

    1.1 試驗材料

    496份公綿羊和14份公巖羊的血液樣本,包括采自北京澳鑫牧業(yè)有限公司的薩福克羊(n=146)、白薩??搜颍╪=91)、東弗里升羊(n=6)、特克塞爾羊(n=72)、南非肉用美利奴羊(n=17)、杜泊羊(n=32),甘肅金昌中天羊業(yè)有限公司的湖羊(n=55),甘肅甘南藏族自治州的藏羊(n=34),寧夏鹽池縣的灘羊(n=43),巖羊采自新疆動物園野生動物保護中心。采集3只湖羊母羊血液(甘肅金昌中天羊業(yè)有限公司)及3只巖羊母羊血液,作為陰性對照。同時在收集血樣時測定薩??搜?、白薩福克羊、東弗里升羊、特克塞爾羊、南非肉用美利奴羊、杜泊羊周歲(周歲±15—20 d)的陰囊圍直徑。

    1.2 DNA的提取和DNA池的構建

    采用常規(guī)的酚-氯仿抽提法提取綿羊血液基因組DNA[16]。用Nano drop 2000(Thermo,美國)檢測提取的DNA濃度和純度,并用1%的瓊脂糖凝膠電泳檢測DNA的質量。將質量檢測合格綿羊和巖羊DNA統(tǒng)一稀釋為50 ng·μL-1,每個個體DNA取10 μL,以同一個品種盡量在一個DNA池的原則,每10個個體混成一個DNA池,共計51個DNA池,混合均勻后-20℃保存?zhèn)溆谩?/p>

    1.3 綿羊Y染色體特異性引物的篩選

    查閱相關文獻可知,目前在牛、牦牛、山羊等牛科動物Y染色體研究有29對引物(表1),分別能擴增出MSY區(qū)的、、、、、等基因。本研究的29對引物均由北京六合華大基因科技有限公司合成。實驗以公羊DNA混合池為模板,以母羊DNA為陰性對照,滅菌水為空白對照,通過PCR擴增確定綿羊Y染色體特異性引物。PCR擴增體系為50 μL:模板2 μL,上、下游引物各1 μL,2×EasyTaq酶25 μL和ddH2O 21μL。擴增反應條件:94℃預變性30 s;94℃變性30 s,退火30 s,72℃延伸1 min/1kb,34個循環(huán);72℃延伸5 min;4℃保存。PCR產物用1%瓊脂糖凝膠電泳檢測。

    表1 引物序信息

    表中僅列出綿羊Y染色體特異性引物的退火溫度The annealing temperatures were listed for those ovine Y-specific primers

    1.4 綿羊Y染色體特異性基因片段SNPs篩選

    將上述鑒定為綿羊Y染色體特異性引物,以DNA混合池為模板的PCR產物送北京六合華大基因科技有限公司進行測序。DNAstar和Chromas軟件用于測序結果分析,DNAMAN用于多序列比對及同源性分析。經過混合池測序,在11片段內鑒定得到一個G>A突變,利用Watcut軟件分析(http://watcut. uwaterloo.ca/template.php?act=snp_new)發(fā)現該突變位點可被限制性內切酶188I識別。酶切反應體系:10×NEB Buffer 3 μL,PCR產物4 μL和5U188I內切酶,加滅菌水至25 μL。酶切反應條件為37℃溫育3 h。取酶切消化液10 μL加樣于2.5%瓊脂糖凝膠上進行電泳。

    1.5 數據處理與統(tǒng)計分析

    利用SPSS 19.0軟件分析基因位點與綿羊睪丸大小的相關性。先對數據進行描述性統(tǒng)計分析,確定是否存在離群值,根據數據特征,利用單因素方差分析(ANOVA)探究基因型效應。

    2 結果

    2.1 綿羊Y染色體特異性引物

    通過PCR擴增,本研究發(fā)現29對引物中6對引物,包括3、6、、、11、6確定為Y染色體特異性引物(圖1),只在公羊DNA上出現擴增條帶,在母羊DNA及水中沒有擴增產物。此外17對引物未能出現擴增條帶或非特異性擴增,6對引物母羊DNA能擴增出條帶(電子版附圖S1)。

    2.2 綿羊Y染色體特異性DNA片段與其他物種的同源性分析

    通過對上述6對特異性引物的PCR產物進行測序,分析發(fā)現3片段大小為614 bp,該片段與山羊第11外顯子和部分CDS區(qū)相似度較高,相似度為98.84%。6片段大小為461 bp,該片段與山羊部分CDS區(qū)相似度較高,相似度為98.84%。片段大小為444 bp,該片段與牛第26外顯子相似度較高,相似度為84.96%。片段大小為256 bp,該片段與牛第19外顯子相似度較高,相似度為81.51%。片段大小為440 bp,該片段與山羊SRY基因序列相似度較高,相似度為90.41%,該基因片段的部分片段(237 bp)與NCBI公布的綿羊的SRY基因CDS區(qū)(AY604735.1)相似度達100%。3片段大小為621 bp,該片段與牛ZFY基因第11外顯子和部分CDS區(qū)相似度較高,相似度為97.67%。

    用DNAMAN生物軟件,將綿羊3、6、、、11、6基因序列與山羊、牛、牦牛和巖羊相應基因序列進行序列比對(電子版附圖S2—S7),計算不同物種基因序列間的同源性大?。ū?)。結果表明,綿羊3基因片段與巖羊、山羊、牛、牦牛各物種間同源性大小分別為98.84%、97.89%、97.08%和95.79%;綿羊SRY6基因片段與巖羊、山羊、牛、牦牛各物種間同源性大小分別為95.25%、93.93%、90.54%和90.32%;綿羊USP9Y基因片段與牛、牦牛各物種間同源性大小分別為84.96%、84.67%;綿羊UTY基因片段與巖羊、黃牛、牦牛各物種間同源性同為96.77%、81.51%和81.51%;綿羊SRY11基因片段與巖羊、山羊、牛、牦牛各物種間同源性大小分別為92.63%、90.41%、84.97%和84.75%;綿羊ZFY6基因片段與巖羊、山羊、牛、牦牛各物種間同源性大小分別為99.06%、89.11%、97.67%和86.76%。由于USP9Y基因是從巖羊母羊DNA中擴增出,表明該片段并不是巖羊Y染色體片段。

    A:ZFY3基因PCR擴增結果;B:SRY6基因PCR擴增結果;C:USP9Y基因PCR擴增結果;D:UTY基因PCR擴增結果;E:SRY11基因PCR擴增結果;F:ZFY6基因PCR擴增結果;1-3公羊DNA,4-6母羊DNA,7-9 ddH2O,M.DNA相對分子質量標準DL2000

    表2 綿羊Y染色體基因與不同物種同源基因的比較

    2.3 綿羊Y染色體SNPs的篩選及基因效應分析

    通過測序,3、6、、、6基因序列均未發(fā)現突變位點。在薩??撕桶姿_??搜蛉后w中,11片段混合池測序出現雙峰,進一步對其中一個混池的每個個體分別進行測序,驗證為G>A的突變(圖2)。對薩福克及白薩??说膫€體進行PCR-RFLP檢測,PCR產物經188I酶切后產生AA和GG 2種基因型(圖3)。AA基因型個體的擴增產物能被188I酶切后將產生4個片段,分別為54、63、151和172 bp;GG基因型個體的擴增產物能被188I酶切后將產生3個片段,分別為63、172和205bp。11片段G>A突變位點在薩??搜颉姿_??搜颉|弗里升羊、特克塞爾羊、南非肉用美利奴羊、杜泊羊、湖羊、藏羊、灘羊公羊群體中的基因型頻率、等位基因頻率見表3。結合混合池測序及酶切結果發(fā)現,在薩??搜?、白薩福克羊中存在AA和GG基因型,其中在白薩福克羊群體中AA基因型頻率為0.253,GG基因型頻率為0.747;在薩??搜蛉后w中AA基因型頻率為0.014,GG基因型頻率為0.986。結果顯示GG基因型在薩??搜蚝桶姿_??搜蛉后w中為優(yōu)勢基因型;而東弗里升羊、特克塞爾羊、南非肉用美利奴羊、杜泊羊、湖羊、藏羊、灘羊公羊群體中只有GG基因型。通過對巖羊個體的混合池測序,在11片段也未發(fā)現突變位點。性狀關聯分析表明,在白薩??搜蛉后w中,GG基因型個體的睪丸大小顯著的高于AA基因型(表4)。

    A. AA基因型; B. GG基因型

    表3 SRY11片段G>A 位點在不同綿羊品種中的基因型和等位基因頻率

    1,3:AA基因型;2,4:GG基因型;M為DNA相對分子質量標準DL2000

    表4 SRY11片段G>A位點與白薩福克羊睪丸大小之間的關聯分析

    3 討論

    在目前研究綿羊起源進化及單倍型分析中,更多的是以線粒體基因遺傳多樣性為基礎的母系遺傳研究,而利用以Y-SNPs為基礎的父系遺傳的研究較少。綿羊的父系及母系遺傳相結合才能更加全面、系統(tǒng)的分析綿羊起源進化、單倍型及分類。本研究篩選得到了6對綿羊Y染色體特異性引物,均為首次擴增得到,與MEADOWS等篩選得到9對綿羊Y染色體特異性引物無重復[9]。(sex-determining region on Y chromosome,)位于哺乳動物Y染色體,是性別決定直接相關的基因。該基因從人的Y染色體上首次分離得到,并被確定為睪丸決定因子[24]。MEADOWS等通過測序在SRY基因5′啟動區(qū)發(fā)現8個SNPs,結合Y染色微衛(wèi)星鑒定綿羊及野生綿羊單倍型[25]。MEADOWS等發(fā)現得到綿羊第一外顯子與牛的同源性為94.7%,并在5′啟動區(qū)鑒定得到了一個A>G的突變[9]。在本次研究擴增得到了兩個SRY基因片段,其與牛的同源性分別為90.54%和84.97%,這與前人的研究相似,說明綿羊和牛SRY基因具有較高的同源性。本研究在11片段鑒定到1個G>A的同義突變,通過188I酶切分析,發(fā)現只有在薩??搜颉姿_??搜蛑杏蠥A和GG 2種基因型,而在其他7個綿羊品種和巖羊中只有GG基因型,暗示該位點可能只存在于薩福克羊和白薩??搜虻难y(tǒng)中,但還需要在更多的綿羊品種群體中進行驗證。此外該突變位點與白薩??搜蛑軞q的睪丸大小緊密相關,即GG基因型個體的睪丸大小顯著大于AA型個體。本研究為首次發(fā)現Y染色體突變與綿羊繁殖性狀相關。該位點在更大的白薩福克羊群體進行驗證后可用于白薩??朔N公羊的選育。

    ZFY基因即Y連鎖的鋅指蛋白因子(zinc finger protein, Y-linked),位于Y染色體短臂(Yp11.3)上,由11個外顯子和1個隨機重復區(qū)域的13個“鋅-指”結構組成,在雄性動物睪丸的生長發(fā)育過程中發(fā)揮作用[26-27]。毛德才研究發(fā)現牦牛、猩猩、人、牛等9個物種間的核苷酸序列有較高同源性[28]。熊勇等研究表明藏系綿羊ZFY基因與豬、摩弗倫羊、狼、人類等其他物種之間具有很高的同源性[11]。LAWSON等通過進化樹分析發(fā)現:綿羊與巖羊、山羊、牛的親緣關系較近[8]。ARKADI克隆得到了整個牛ZFY基因整個編碼區(qū)及馬和豬ZFY基因編碼區(qū)部分序列,通過比對發(fā)現ZFY基因編碼區(qū)在3個物種的同源性較高[29]。本研究發(fā)現綿羊3和6片段與巖羊、山羊、牛、牦牛等物種間同源性大小在89.11%—99.6%,表明該基因在不同物種間也具有較高的同源性,與前人的研究基本一致。AASEN等[30]利用片段的擴增及酶切鑒別人、牛、綿羊和山羊的胚胎性別,該研究的結果也可用于綿羊胚胎性別的早期鑒定。

    SATOU等[31]研究發(fā)現USP9Y基因的突變影響小鼠睪丸重量。有學者發(fā)現在牛USP9Y基因第26內含子發(fā)現了一個81 bp的插入缺失,并用于構建牛Y染色體單倍型及區(qū)分瘤牛和普通牛[21,32-33]。本研究在USP9Y基因片段中沒有檢測到遺傳變異,今后需要在更多的綿羊品種進行Y-SNPs掃描。先前的研究利用BAC的方法在絨山羊上鑒定到USP9Y基因348 bp的片段,其與家貓、人、恒河猴、黑猩猩等的同源性在89%以上[34]。本研究擴增得到的綿羊USP9Y基因與牛、牦牛各物種間同源性大小分別為84.96%、84.67%,表明該基因在進化過程中具有一定的保守性。USP9Y基因在巖羊母羊DNA中擴增出,表明該片段并不是巖羊Y染色體片段,因此本研究的結果可作為區(qū)分巖羊和綿羊的分子標記。

    UTY基因即Y連鎖的泛轉錄三角形四肽重復因子,在人上的研究發(fā)現UTY基因有6個不同的mRNA轉錄本,UTY蛋白被認為是未成熟的組織適合性抗原,可啟動男性特異性細胞毒性T淋巴細胞反應[35-36]。MEADOWS等擴增得到了UTY基因外顯子1,通過比對發(fā)現與牛的同源性為98.1%[9]。本研究擴增得到的綿羊UTY基因與巖羊、牛、牦牛各物種間同源性達81.51%—96.77%,表明該基因在進化過程中相對保守。GINJA等擴增牛UTY基因19外顯子,發(fā)現有一個C>A的突變[23]。然而在本研究中在UTY基因上并未檢測到突變位點。

    本研究首次在發(fā)現一個G>A的突變,在3、6、、、6片段均未發(fā)現突變,其原因可能與樣本量大小及品種差異有關。但是本研究的結果為今后我國綿羊品種父系遺傳多樣性的研究奠定了基礎。當然,為了系統(tǒng)和全面研究綿羊Y染色體,必須依賴綿羊Y染色體完整的序列,完成綿羊Y染色體測序是研究綿羊父系遺傳多樣性和綿羊雄性繁殖力育種亟待解決的難題。

    4 結論

    本研究首次擴增得到綿羊、、和4個基因的6個片段。通過比對發(fā)現綿羊3、6、、、6片段與巖羊、牛、山羊、牦牛的具有較高的同源性。首次在11片段發(fā)現G>A的突變,該突變位點只存在于薩??撕桶姿_??斯蛉后w中,并與白薩??搜虿G丸大小緊密相關。

    A:2 PCR擴增結果;B:4 PCR擴增結果;C:5 PCR擴增結果;D:9 PCR擴增結果;E:1 PCR擴增結果;F:9 PCR擴增結果;1-3公羊DNA,4-6母羊DNA,7-9ddH2O,M.DNA相對分子質量標準DL2000

    A.PCR products of3 fragment; B.PCR products of6 fragment; C.PCR products offragment; D.PCR products offragment; E.PCR products of11 fragment; F.PCR products of3 fragment;1-3 rams DNA;4-6 ewes DNA;7-9 ddH2O;M.DL2000 marker

    圖S1 綿羊Y染色體非特異性PCR擴增結果

    Fig. S1 PCR products of Y nonspecific fragment in sheep

    圖S2 不同物種3基因序列比對結果

    Fig. S2 Comparison of ovine3 gene sequence with other species

    圖S3 不同物種6基因序列比對結果

    Fig. S3 Comparison of ovine6 gene sequence with other species

    圖S4 不同物種基因序列比對結果

    Fig. S4 Comparison of ovinegene sequence with other specie

    圖S5 不同個物種基因序列比對結果

    Fig. S5 Comparison of ovinegene sequence with other specie

    圖S6 不同物種11基因序列比對結果

    Fig. S6 Comparison of ovine11 gene sequence with other specie

    圖S7 不同物種6基因序列比對結果

    Fig. S7 Comparison of ovine6 gene sequence with other specie

    [1] MIN Z, PENG W F, YANG G L, LV F H, LIU M J, LI W R, LIU Y G, LI J Q, FENG W, SHEN Z Q. Y chromosome haplotype diversity of domestic sheep () in northern Eurasia.2014, 45(6): 903-907.

    [2] FERENCAKOVIC M, CURIK I, PéREZ-PARDAL L, ROYO L J, CUBRIC-CURIK V, FERNáNDEZ I, áLVAREZ I, KOSTELIC A, SPREM N, KRAPINEC K. Mitochondrial DNA and Y-chromosome diversity in East Adriatic sheep.2013, 44(2): 184-192.

    [3] CHANG T C, YANG Y, RETZEL E F, LIU W S. Male-specific region of the bovine Y chromosome is gene rich with a high transcriptomic activity in testis development.2013, 110(30): 12373-12378.

    [4] GRAVES J A. Sex chromosome specialization and degeneration in mammals., 2006, 124(5):901-914.

    [5] KIM T Y, SHIN K D, SONG G M. The male-specific region of the human Y chromosome is a mosaic of discrete sequence classes.,2003, 423(6942): 825-837.

    [6] PAYEN E, PAILHOUX E, ABOU M R, GIANQUINTO L, KIRSZENBAUM M, LOCATELLI A, COTINOT C. Characterization of ovine SRY transcript and developmental expression of genes involved in sexual differentiation.,1996, 40(3): 567-575.

    [7] FENG J, LAJIA C, TAYLOR D J, WEBSTER M S. Genetic distinctiveness of endangered dwarf blue sheep (schaeferi): evidence from mitochondrial control region and Y-linkedintron sequences.,2001, 92(1): 9-15.

    [8] LAWSON L J, HEWITT G M. Comparison of substitution rates inandintrons of sheep and goat related species supports the hypothesis of male-biased mutation rates.2002, 54(1), 54-61.

    [9] MEADOWS J R, HAWKEN R J, KIJAS J W. Nucleotide diversity on the ovine Y chromosome.2004, 35(5): 379-385.

    [10] 劉帥兵. 綿羊基因干擾載體篩選及性控效果的驗證[D]. 石河子: 石河子大學, 2014.

    LIU S B. The RNAi vector screening and verification of the effects of sexual control fromgene[D]. Shihezi: Shihezi University, 2014. (in Chinese)

    [11] 熊勇, 徐亞歐. 藏系綿羊/基因片段進化分析. 生物技術2010, 20(6): 1-5.

    XIONG Y, XU Y O. Evolution analysis the Tibetan sheepgene fragment.2010, 20(6): 1-5. (in Chinese)

    [12] 蔣利, 王康環(huán), 王海, 任冰冰, 劉光偉, 徐亞歐. 阿勒泰羊、藏綿羊(賈洛類群、歐拉類群)基因序列分析及親緣關系研究. 西南民族大學學報(自然科學版)2013, 39(4): 487-494.

    JIANG L, WANG K H, WANG H, REN B B, LIU G W, XU Y O. Sequence analysis of the Sry gene and genetic relationship of Altay Sheep, Tibetan sheep(Jialuo population, Oula population).2013, 39(4): 487-494. (in Chinese)

    [13] HUGHES J F, SKALETSKY H, PYNTIKOVA T, GRAVES T A, DAALEN S K M V, MINX P J, FULTON R S, MCGRATH S D, LOCKE D P, FRIEDMAN C. Chimpanzee and human Y chromosomes are remarkably divergent in structure and gene content.2010, 463(7280): 536-539.

    [14] HUGHES J F, ROZEN S. Genomics and genetics of human and primate y chromosomes.2012, 13(1): 83-108.

    [15] SOH Y Q S, ALF?LDI J, PYNTIKOVA T, BROWN L, GRAVES T, MINX P, FULTON R, KREMITZKI C, KOUTSEVA N, MUELLER J. Sequencing the mouse Y chromosome reveals convergent gene acquisition and amplification on both sex chromosomes.2014, 159(4): 800-813.

    [16] SAMBROCK J, RUSSEL D., 3rd ed. New York: Cold Spring Harbor Laboratory Press, 2001: 463-485.

    [17] LI R, WANG S Q, XU S Y, HUANG J P, WANG F Q, MA Z J, DANG R H, LAN X Y, CHEN H, LEI C Z. Novel Y-chromosome polymorphisms in chinese domestic yak.2014, 45(3): 449-452.

    [18] NIJMAN I J, VAN BOXTEL D C J, VAN CANN L M, MARNOCH Y, CUPPEN E, LENSTRA J A. Phylogeny of Y chromosomes from bovine species.2008, 24(5): 723-726.

    [19] CHENG H, SHI H, ZHOU R, GUO Y, LI L, LIU J, YONG J, KUDO T, SUTOU S. Characterization of Bovidae sex-determining gene.2001, 33(6): 687-694.

    [20] VERKAARE L, NIJMAN I J, BEEKE M, HANEKAMP E, LENSTRA J A. Maternal and paternal lineages in cross-breeding bovine species. Has wisent a hybrid origin?2004, 21(7): 1165-1170.

    [21] BONFIGLIO S, DE G A, TESFAYE K, GRUGNI V, SEMINO O, FERRETTI L. A novelpolymorphism allowing a rapid and unambiguous classification ofY chromosomes into haplogroups.2012, 43(5): 611-613.

    [22] G?THERSTR?M A, ANDERUNG C, HELLBOR L G, ELBURG R, SMITH C, DAN G B, ELLEGREN H. Cattle Domestication in the Near East Was Followed by Hybridization with Aurochs Bulls in Europe.2005, 272(1579): 2345-2350.

    [23] GINJA C, TELO D G L, PENEDO M C. Y chromosome haplotype analysis in Portuguese cattle breeds using SNPs and STRs.2015, 100(2): 148.

    [24] ANDREW P B, SINCLAIR H, PALMER M S, HAWKINS J R, GRIFFITHS B L, SMITH M J, FOSTER J W, FRISCHAUF A F, LOVELL-BADGE R, GOOGFELLOW P N. A gene from the human sex-determining region encodes a protein with homology to a conserved DNA-binding motif.1900, 346(6281): 5.

    [25] MEADOWS J R, KIJAS J W. Re-sequencing regions of the ovine Y chromosome in domestic and wild sheep reveals novel paternal haplotypes.2009, 40(1): 119-123.

    [26] AO A, ERICKSON R P, WINSTON R M, HANDYSIDE A H. Transcription of paternal Y-linked genes in the human zygote as early as the pronucleate stage.1994, 2(4): 281-287.

    [27] NAGAMINE C M, CHAN K, KOZAK C A, LAU Y F. Chromosome mapping and expression of a putative testis-determining gene in mouse.1989, 243(4887): 80-83.

    [28] 毛德才, 徐亞歐, 熊勇, 曹特, 蘇暢. 牦牛與其他物種基因片段間的進化關系. 中國生物化學與分子生物學報2007, 23(6): 450-456.

    MAO D C, XU Y O, XIONG Y, CAO T, SU C. Evolution analysis ofgene fragments between yak and other species., 2007, 23(6): 450-456. (in Chinese)

    [29] POLOUMIENKO A. Cloning and comparative analysis of the bovine, porcine, and equine sex chromosome genesand.2004, 47(1): 74-83.

    [30] AASEN E, MEDRANO J F. Amplification of theandgenes for sex identification in humans, cattle, sheep and goats.1990, 8(12): 1279.

    [31] SATOU K, SUTO J I. Effect of the Y chromosome on testis weight in mice.2015, 77(6): 753.

    [32] 馬志杰, 胡雙龍, 李瑞哲, 陳生梅, 雷初朝, 晁生玉. 柴達木黃牛Y染色體單倍型組構成及父系起源——基于基因多態(tài)性的分析. 中國畜牧雜志2017, 53(6): 36-39.

    MA Z J, HU S L, LI R Z , CHEN S M, LEI C Z, CHAO S Y. Y Chromosome haplogroups and paternal origin of qaidam cattle: Based on the polymorphism analysis of USP9Y gene., 2017, 53(6): 36-39. (in Chinese)

    [33] LI R, XIE W M, CHANG Z H, WANG S Q, DANG R H, LAN X Y, CHEN H, LEI C Z. Y chromosome diversity and paternal origin of Chinese cattle.2013, 40(12): 6633-6636.

    [34] 肖紅梅, 劉志紅, 張文廣, 李金泉. 絨山羊基因的BAC篩選與鑒定. 中國畜牧獸醫(yī)2012, 39(4): 31-34.

    XIAO H M, LIU Z H, ZHANG W G, LI J Q. Screening and identification of cashmere goatsgene., 2012, 39(4): 31-34. (in Chinese)

    [35] LAASER I, THEIS F J, DE ANGELIS M H, KOLB H J, ADAMSKI J. Huge splicing frequency in human Y chromosomalgene.2011, 15(3): 141.

    [36] IVANOV R, HOL S, AARTS T, HAGENBEEK A, SLAGER E H, EBELING S.-specific TCR-transfer generates potential graft-versus-leukaemia effector T cells.2005, 129(3): 392.

    (責任編輯 林鑒非)

    Screening of Y Chromosome Specific Primers and Y-SNPs in Sheep

    CAO XueTao1, PEI ShengWei1, ZHANG Jin3, LI FaDi1,2, Li Gang3, LI WanHong1, YUE XiangPeng1

    (1State Key Laboratory of Grassland Agro-ecosystems, Lanzhou University/Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs/ College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730020;2Engineering Laboratory of Mutton Sheep Breeding and Re-production Biotechnology in Gansu Province, Minqin 733300, Gansu;3Animal Breeding Centre of Gansu Province, Wuwei 733000, Gansu)

    【Objective】The male specific region of the mammalian Y chromosome (MSY)does not recombine with X chromosome during meiosis process, which is an important genetic resource for analyzing paternal genetic diversity due to its strict father to son inherited character. In addition, most of genes on MSY exclusively or predominantly express in the testis, indicating they may play essential roles in spermatogenesis and male reproduction. Since it is extremely difficult to sequence of entire Y chromosome, many species have very few Y chromosome sequences. Therefore, the current study was conducted to select ovine Y chromosome specific primers and Y-SNP based on the previous primer information used in bovidae species, and to compare the Y-fragment sequences similarity among sheep, bovine, goat, yak and bharal. Meanwhile, the Y-SNP was associated with sheep scrotal circumference to supply scientific basis for constructing ovine Y-haplotypes, identifying molecular markers for embryo sex and male reproductive traits in the future.【Method】Based on the investigation of available references about bovidae Y chromosome, 29 pairs of Y-primers reported in cattle, yak, goat were selected to amplify rams DNA using ewes DNA and ddH2O as negative controls. Subsequently, the Y-SNPs within ovine Y-specific fragment identified of different sheep breeds were investigated by DNA sequencing of DNA pooling and PCR-RFLP methods, including Suffolk sheep (n=146),White Suffolk sheep (n=91),East Friesian sheep (n=6), Texel sheep (n=72), South African Mutton Merino (n=17), Dorper sheep (n=32), Hu sheep (n=55),Tibetan sheep (n=34), Tan sheep (n=43), and bharal (n=14). Chromas and DNASTAR were used to analyze the results of DNA-pool sequencing, and DNAman was used for homology analysis of yak, goat, cattle and bharal. Meanwhile, the correlation analysis between the11 gene fragment polymorphisms and the testis size was performed by SPSS 19.0.【Result】The results showed that 6 out of 29 pairs of primers analyzed were ovine Y-specific, which could amplify3,6,,,11 and6 fragments, respectively. However, 17 pairs of primers failed to show amplification bands, and 6 pairs of primers showed amplified bands in the DNA of the ewes. The similarity of them among sheep, bharal, cattle, goat and yak ranged from 81.51% to 98.84%. In addition, a Y-SNP (G>A) with in11 fragment was first identified in the Suffolk and white Suffolk sheep.According to RCR-RFLP analysis, two genotypes (AA, GG) were detected in the Suffolk sheep and white Suffolk sheep, while only the GG genotype was found in the other seven sheep breeds. The genotypic frequencies of the GG and AA were 0.747 and 0.253 in White Suffolk sheep, respectively, while they were 0.986 and 0.014 in Suffolk sheep, indicating the dominant genotype was GG genotype in White Suffolk sheep and Suffolk sheep. Association analysis suggested that the testis size of the GG genotype was significantly higher than those of the AA genotype in the white Suffolk sheep population (=0.029).【Conclusion】In this study, six pairs of ovine Y-specific primers were identified, and the Y-linked fragment identified in ovine showed a high similarity with cattle, goats and yaks, indicating certain conservation in the evolutionary process. In addition, a Y-SNP was found to be specific in White Suffolk sheep and Suffolk sheep, which was closely associated with the testis size of white Suffolk sheep.

    Y chromosome; sheep; single nucleotide polymorphism (SNP)

    2017-12-25;

    2018-04-17

    甘肅省農業(yè)生物技術研究與應用開發(fā)項目(GNSW-2015-24,GNSW-2016-22)、甘肅省重點研發(fā)計劃(17YF1NA066)、甘肅省農業(yè)科技創(chuàng)新項目(GNCX-2014-41)、國家肉羊產業(yè)技術體系(CARS-38)和長江學者和創(chuàng)新團隊發(fā)展計劃(IRT13019)

    曹學濤,E-mail:caoxt16@ lzu.edu.cn。信作者樂祥鵬,E-mail:lexp@lzu.edu.cn

    10.3864/j.issn.0578-1752.2018.15.0014

    猜你喜歡
    薩???/a>巖羊綿羊
    巖羊
    戰(zhàn)勝狼的巖羊
    中外文摘(2019年10期)2019-11-12 19:59:26
    數綿羊
    戰(zhàn)勝狼的巖羊
    戰(zhàn)勝狼的巖羊
    中外管理(2019年5期)2019-05-22 13:27:32
    數綿羊
    奔跑的綿羊
    幼兒畫刊(2018年7期)2018-07-24 08:26:10
    薩福克羊種質特性及利用狀況
    薩???、德國美利奴羊與小尾寒羊雜交羔羊育肥性能的研究
    草食家畜(2016年3期)2016-07-25 08:58:40
    巧計得綿羊
    亚洲精品一卡2卡三卡4卡5卡| 日韩中文字幕欧美一区二区| 免费在线观看视频国产中文字幕亚洲| 咕卡用的链子| 不卡av一区二区三区| 在线十欧美十亚洲十日本专区| 国产精品乱码一区二三区的特点 | 少妇被粗大的猛进出69影院| 精品国产美女av久久久久小说| 欧美激情 高清一区二区三区| 在线天堂中文资源库| 亚洲人成77777在线视频| 中文字幕另类日韩欧美亚洲嫩草| 免费看十八禁软件| 国内久久婷婷六月综合欲色啪| 精品高清国产在线一区| 久久欧美精品欧美久久欧美| 国产精品自产拍在线观看55亚洲| www国产在线视频色| 中文字幕av电影在线播放| 亚洲中文字幕日韩| 亚洲精品在线美女| 自线自在国产av| 国产高清激情床上av| 国产精品久久久久久精品电影 | 亚洲熟妇中文字幕五十中出| 美女高潮到喷水免费观看| 少妇裸体淫交视频免费看高清 | 亚洲国产日韩欧美精品在线观看 | 亚洲男人的天堂狠狠| 国产在线精品亚洲第一网站| 亚洲精品中文字幕在线视频| 国产一区二区三区综合在线观看| 真人一进一出gif抽搐免费| 成年女人毛片免费观看观看9| 欧美日本视频| 巨乳人妻的诱惑在线观看| 91av网站免费观看| 国产精品1区2区在线观看.| 纯流量卡能插随身wifi吗| 美女高潮喷水抽搐中文字幕| 两性夫妻黄色片| 正在播放国产对白刺激| 免费在线观看影片大全网站| 亚洲精品美女久久久久99蜜臀| 嫁个100分男人电影在线观看| www.熟女人妻精品国产| 亚洲欧美激情在线| 久热这里只有精品99| 久久久久久久久久久久大奶| 制服人妻中文乱码| 黄色 视频免费看| 99精品欧美一区二区三区四区| a在线观看视频网站| 中出人妻视频一区二区| 成人国产综合亚洲| 国产私拍福利视频在线观看| 好男人在线观看高清免费视频 | 12—13女人毛片做爰片一| 88av欧美| 黄色片一级片一级黄色片| 丁香六月欧美| 精品电影一区二区在线| 叶爱在线成人免费视频播放| www国产在线视频色| 欧美性长视频在线观看| 午夜福利一区二区在线看| 久久天躁狠狠躁夜夜2o2o| 伦理电影免费视频| 久久精品91蜜桃| 亚洲视频免费观看视频| 午夜精品久久久久久毛片777| 脱女人内裤的视频| 成人av一区二区三区在线看| 亚洲狠狠婷婷综合久久图片| 校园春色视频在线观看| 国产视频一区二区在线看| 精品乱码久久久久久99久播| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成+人综合+亚洲专区| 久久久国产欧美日韩av| 99国产精品一区二区蜜桃av| 欧美大码av| 中文亚洲av片在线观看爽| 宅男免费午夜| 国产国语露脸激情在线看| 丁香六月欧美| 一级a爱片免费观看的视频| 国产成人精品在线电影| 精品免费久久久久久久清纯| 国产av又大| 精品日产1卡2卡| 精品一区二区三区av网在线观看| 好看av亚洲va欧美ⅴa在| 欧美一区二区精品小视频在线| 在线免费观看的www视频| 国产成人精品在线电影| 欧美成狂野欧美在线观看| 久久草成人影院| 我的亚洲天堂| 精品不卡国产一区二区三区| 岛国视频午夜一区免费看| 欧美乱妇无乱码| 看免费av毛片| 欧美av亚洲av综合av国产av| 久久久久九九精品影院| 免费一级毛片在线播放高清视频 | 久久国产乱子伦精品免费另类| 色在线成人网| 久久国产乱子伦精品免费另类| 男女午夜视频在线观看| 国产精品久久久久久精品电影 | 51午夜福利影视在线观看| 嫁个100分男人电影在线观看| 99久久综合精品五月天人人| 国产激情欧美一区二区| 国产精品99久久99久久久不卡| 久久人人97超碰香蕉20202| 欧美午夜高清在线| 免费女性裸体啪啪无遮挡网站| 在线观看日韩欧美| 欧美另类亚洲清纯唯美| 欧美大码av| 看免费av毛片| 一本大道久久a久久精品| 黑人巨大精品欧美一区二区mp4| 精品乱码久久久久久99久播| 日韩有码中文字幕| 色婷婷久久久亚洲欧美| 深夜精品福利| 桃色一区二区三区在线观看| 手机成人av网站| 久久精品国产清高在天天线| 午夜免费激情av| 夜夜夜夜夜久久久久| av片东京热男人的天堂| 国产精品 欧美亚洲| 成人18禁在线播放| 午夜日韩欧美国产| 精品一区二区三区视频在线观看免费| 午夜福利欧美成人| 久久久久久大精品| 国产成人影院久久av| 亚洲欧美精品综合久久99| 91精品国产国语对白视频| 激情视频va一区二区三区| 亚洲片人在线观看| 久久久精品欧美日韩精品| 禁无遮挡网站| www.精华液| 黑人操中国人逼视频| 99久久久亚洲精品蜜臀av| 久久国产乱子伦精品免费另类| 美国免费a级毛片| 亚洲专区中文字幕在线| 亚洲av电影在线进入| 三级毛片av免费| 日日爽夜夜爽网站| 亚洲精品国产区一区二| 1024香蕉在线观看| 一边摸一边抽搐一进一小说| 国产97色在线日韩免费| 国产成人精品久久二区二区免费| tocl精华| 黑人欧美特级aaaaaa片| 一边摸一边做爽爽视频免费| 波多野结衣av一区二区av| 国产精品二区激情视频| 精品熟女少妇八av免费久了| 成年版毛片免费区| 国产成人精品在线电影| xxx96com| 国产不卡一卡二| 91国产中文字幕| 久久精品国产综合久久久| 在线观看免费视频网站a站| 国产色视频综合| 热re99久久国产66热| 国产一区在线观看成人免费| 中文字幕人妻熟女乱码| 美国免费a级毛片| 在线永久观看黄色视频| av片东京热男人的天堂| 久久婷婷成人综合色麻豆| 国产伦人伦偷精品视频| 怎么达到女性高潮| 亚洲美女黄片视频| 午夜精品久久久久久毛片777| 又黄又粗又硬又大视频| 天堂√8在线中文| 大码成人一级视频| 国产午夜福利久久久久久| or卡值多少钱| 两个人看的免费小视频| 日本欧美视频一区| 如日韩欧美国产精品一区二区三区| 99久久99久久久精品蜜桃| 久热这里只有精品99| 欧美日韩福利视频一区二区| 亚洲色图av天堂| 99在线人妻在线中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 999精品在线视频| 香蕉丝袜av| 亚洲天堂国产精品一区在线| 波多野结衣av一区二区av| 午夜福利视频1000在线观看 | 19禁男女啪啪无遮挡网站| 中文字幕色久视频| 午夜福利高清视频| 久久亚洲真实| 成人永久免费在线观看视频| 亚洲精品国产精品久久久不卡| 国产亚洲欧美精品永久| 免费看a级黄色片| 18美女黄网站色大片免费观看| 波多野结衣巨乳人妻| 欧美黑人精品巨大| 少妇的丰满在线观看| 一区二区三区激情视频| 免费搜索国产男女视频| 亚洲午夜理论影院| 99在线人妻在线中文字幕| 亚洲avbb在线观看| 午夜福利18| 在线观看免费日韩欧美大片| 手机成人av网站| videosex国产| 91老司机精品| 亚洲激情在线av| 欧美丝袜亚洲另类 | 午夜老司机福利片| 久久国产精品影院| 精品免费久久久久久久清纯| 中文字幕av电影在线播放| 一区二区日韩欧美中文字幕| 在线av久久热| 好男人电影高清在线观看| 免费看十八禁软件| 不卡av一区二区三区| 国产黄a三级三级三级人| 成熟少妇高潮喷水视频| 亚洲自拍偷在线| 真人一进一出gif抽搐免费| 美女国产高潮福利片在线看| 香蕉丝袜av| 免费观看精品视频网站| 日本 av在线| 三级毛片av免费| 午夜两性在线视频| av超薄肉色丝袜交足视频| 搡老熟女国产l中国老女人| 中文字幕人妻丝袜一区二区| 制服人妻中文乱码| av视频免费观看在线观看| 欧美丝袜亚洲另类 | 成人亚洲精品av一区二区| 亚洲专区中文字幕在线| 久9热在线精品视频| 国产精品免费视频内射| 大型av网站在线播放| 99久久99久久久精品蜜桃| 人人妻人人爽人人添夜夜欢视频| 欧美成人免费av一区二区三区| 99久久久亚洲精品蜜臀av| 久久精品91无色码中文字幕| 免费在线观看视频国产中文字幕亚洲| 国产精品自产拍在线观看55亚洲| 99国产精品99久久久久| 日韩三级视频一区二区三区| 黑人欧美特级aaaaaa片| 91麻豆精品激情在线观看国产| 成年女人毛片免费观看观看9| 大型黄色视频在线免费观看| 国产亚洲av嫩草精品影院| 757午夜福利合集在线观看| 国产成人精品在线电影| 我的亚洲天堂| 人人妻人人爽人人添夜夜欢视频| 亚洲天堂国产精品一区在线| 他把我摸到了高潮在线观看| 日本在线视频免费播放| 久久精品aⅴ一区二区三区四区| 国产区一区二久久| 老熟妇仑乱视频hdxx| 午夜免费激情av| 欧美不卡视频在线免费观看 | 国产三级黄色录像| 一级毛片女人18水好多| 十八禁网站免费在线| 亚洲精品国产色婷婷电影| 天天躁狠狠躁夜夜躁狠狠躁| 99久久精品国产亚洲精品| 午夜视频精品福利| 一级作爱视频免费观看| 欧美日韩黄片免| 一级片免费观看大全| 无限看片的www在线观看| 久久久国产精品麻豆| 日韩大尺度精品在线看网址 | 91精品国产国语对白视频| 人人妻人人爽人人添夜夜欢视频| 天堂影院成人在线观看| 两个人免费观看高清视频| 亚洲成国产人片在线观看| 91字幕亚洲| 禁无遮挡网站| 999久久久精品免费观看国产| 极品教师在线免费播放| 日日摸夜夜添夜夜添小说| 国产精品免费一区二区三区在线| 丰满人妻熟妇乱又伦精品不卡| 日韩成人在线观看一区二区三区| 成人国语在线视频| 国产野战对白在线观看| 首页视频小说图片口味搜索| 最好的美女福利视频网| 免费在线观看亚洲国产| 女人高潮潮喷娇喘18禁视频| 久久这里只有精品19| 老汉色∧v一级毛片| 日本黄色视频三级网站网址| 在线观看午夜福利视频| 国产亚洲精品久久久久久毛片| 国产高清有码在线观看视频 | 国产97色在线日韩免费| 色在线成人网| 可以免费在线观看a视频的电影网站| 国产精品影院久久| 国产麻豆69| 久久香蕉精品热| 日本 欧美在线| 少妇熟女aⅴ在线视频| 香蕉丝袜av| 国产熟女xx| 久久久精品欧美日韩精品| 亚洲色图综合在线观看| 亚洲,欧美精品.| 美女扒开内裤让男人捅视频| 久久亚洲精品不卡| 精品久久蜜臀av无| 免费少妇av软件| 女人被狂操c到高潮| 国产三级在线视频| 欧美中文日本在线观看视频| 久久久久久人人人人人| 国产亚洲av嫩草精品影院| 亚洲欧美激情在线| 亚洲色图 男人天堂 中文字幕| 黑人巨大精品欧美一区二区mp4| 亚洲国产欧美一区二区综合| 九色亚洲精品在线播放| 午夜福利在线观看吧| 精品福利观看| 最近最新免费中文字幕在线| 亚洲人成网站在线播放欧美日韩| 高清毛片免费观看视频网站| 成人特级黄色片久久久久久久| 日韩高清综合在线| 中文亚洲av片在线观看爽| 又大又爽又粗| 不卡一级毛片| 免费在线观看黄色视频的| 丝袜在线中文字幕| 国产成人av教育| 欧美一区二区精品小视频在线| 亚洲一卡2卡3卡4卡5卡精品中文| 嫩草影视91久久| 99久久国产精品久久久| 国产亚洲精品第一综合不卡| 精品欧美一区二区三区在线| 久久精品aⅴ一区二区三区四区| 法律面前人人平等表现在哪些方面| 午夜a级毛片| 成熟少妇高潮喷水视频| 国产av一区二区精品久久| 日本三级黄在线观看| 国产av一区二区精品久久| 黄色丝袜av网址大全| 亚洲视频免费观看视频| 国产成人av激情在线播放| 天堂动漫精品| 一级毛片高清免费大全| 亚洲av美国av| 亚洲第一欧美日韩一区二区三区| ponron亚洲| 日韩大码丰满熟妇| 亚洲国产精品合色在线| 久久天躁狠狠躁夜夜2o2o| 女同久久另类99精品国产91| 真人一进一出gif抽搐免费| 99久久久亚洲精品蜜臀av| 免费在线观看黄色视频的| 日韩精品青青久久久久久| 国产麻豆69| 国产三级在线视频| 欧美日韩亚洲国产一区二区在线观看| 久热这里只有精品99| 级片在线观看| 岛国视频午夜一区免费看| 日韩大尺度精品在线看网址 | 色老头精品视频在线观看| av网站免费在线观看视频| 久热爱精品视频在线9| 国产精品久久视频播放| 精品一区二区三区视频在线观看免费| 午夜精品久久久久久毛片777| 欧美日韩黄片免| www.www免费av| 亚洲最大成人中文| 亚洲国产高清在线一区二区三 | 色婷婷久久久亚洲欧美| 国产三级黄色录像| 人妻丰满熟妇av一区二区三区| 国产免费av片在线观看野外av| av片东京热男人的天堂| 好男人电影高清在线观看| 一区在线观看完整版| 国产精品99久久99久久久不卡| 久久久久久国产a免费观看| 伦理电影免费视频| av天堂在线播放| 国产欧美日韩综合在线一区二区| 欧美日韩黄片免| 久久天躁狠狠躁夜夜2o2o| netflix在线观看网站| 一边摸一边抽搐一进一出视频| 欧美日韩亚洲综合一区二区三区_| 夜夜躁狠狠躁天天躁| 少妇被粗大的猛进出69影院| 性欧美人与动物交配| 黄网站色视频无遮挡免费观看| 亚洲aⅴ乱码一区二区在线播放 | 一级,二级,三级黄色视频| 99国产精品一区二区三区| 午夜福利免费观看在线| 50天的宝宝边吃奶边哭怎么回事| av视频免费观看在线观看| 91麻豆av在线| 精品人妻1区二区| 欧美日韩福利视频一区二区| 禁无遮挡网站| 天天一区二区日本电影三级 | 看黄色毛片网站| 国产三级黄色录像| 精品国内亚洲2022精品成人| 亚洲久久久国产精品| 国产成人啪精品午夜网站| 啦啦啦免费观看视频1| 欧美黄色淫秽网站| 99国产极品粉嫩在线观看| 免费在线观看黄色视频的| 老汉色∧v一级毛片| 久久久久久久久久久久大奶| 757午夜福利合集在线观看| av在线天堂中文字幕| 免费一级毛片在线播放高清视频 | 国产精品一区二区在线不卡| 亚洲精品久久成人aⅴ小说| 国产单亲对白刺激| 波多野结衣高清无吗| 一级作爱视频免费观看| 一区二区三区激情视频| 最近最新免费中文字幕在线| 日韩精品中文字幕看吧| 少妇裸体淫交视频免费看高清 | 国产一区二区激情短视频| 亚洲在线自拍视频| 国产91精品成人一区二区三区| 成熟少妇高潮喷水视频| 国产三级在线视频| 免费在线观看黄色视频的| 日韩免费av在线播放| 两性夫妻黄色片| 亚洲中文字幕日韩| 电影成人av| 岛国在线观看网站| 91av网站免费观看| 国产亚洲精品av在线| 亚洲无线在线观看| 亚洲精品av麻豆狂野| 黄色视频不卡| 黄色毛片三级朝国网站| 啦啦啦 在线观看视频| 日韩大尺度精品在线看网址 | 夜夜躁狠狠躁天天躁| 亚洲成av片中文字幕在线观看| 国产aⅴ精品一区二区三区波| 国产亚洲av高清不卡| 18禁美女被吸乳视频| 精品国产国语对白av| 女生性感内裤真人,穿戴方法视频| 最新美女视频免费是黄的| 精品国产美女av久久久久小说| 日本三级黄在线观看| 99re在线观看精品视频| 精品国产一区二区久久| 亚洲专区中文字幕在线| 69精品国产乱码久久久| 久久香蕉精品热| 电影成人av| 激情在线观看视频在线高清| 色综合站精品国产| 国产99久久九九免费精品| 国产一区二区三区综合在线观看| 最新美女视频免费是黄的| 亚洲激情在线av| 国产片内射在线| 91麻豆精品激情在线观看国产| 欧美大码av| 91麻豆精品激情在线观看国产| 国产片内射在线| 国产精品,欧美在线| 久久影院123| 精品国产美女av久久久久小说| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲自拍偷在线| 可以免费在线观看a视频的电影网站| 天天躁狠狠躁夜夜躁狠狠躁| 少妇的丰满在线观看| 啦啦啦免费观看视频1| 美女扒开内裤让男人捅视频| 久热这里只有精品99| 国产片内射在线| 精品乱码久久久久久99久播| 欧美一级a爱片免费观看看 | 欧美日韩瑟瑟在线播放| 亚洲aⅴ乱码一区二区在线播放 | av片东京热男人的天堂| 国产野战对白在线观看| 国产单亲对白刺激| 国产真人三级小视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久国产a免费观看| 老熟妇乱子伦视频在线观看| 午夜精品久久久久久毛片777| 99re在线观看精品视频| 一级毛片精品| 一级,二级,三级黄色视频| 韩国精品一区二区三区| 在线播放国产精品三级| 极品教师在线免费播放| 动漫黄色视频在线观看| 狂野欧美激情性xxxx| 亚洲熟妇中文字幕五十中出| www日本在线高清视频| 一级作爱视频免费观看| 亚洲片人在线观看| 一级毛片高清免费大全| 一边摸一边抽搐一进一小说| 国产精品影院久久| 成人18禁高潮啪啪吃奶动态图| 香蕉丝袜av| 身体一侧抽搐| 亚洲avbb在线观看| 90打野战视频偷拍视频| 久热爱精品视频在线9| 日韩一卡2卡3卡4卡2021年| 欧美激情 高清一区二区三区| 亚洲专区字幕在线| 麻豆国产av国片精品| 午夜a级毛片| 热re99久久国产66热| 这个男人来自地球电影免费观看| 亚洲精品在线观看二区| 欧美久久黑人一区二区| 国产高清有码在线观看视频 | 黑人操中国人逼视频| 大香蕉久久成人网| 啦啦啦免费观看视频1| 久久 成人 亚洲| 1024视频免费在线观看| 村上凉子中文字幕在线| 中文字幕久久专区| 亚洲国产欧美一区二区综合| av中文乱码字幕在线| 淫妇啪啪啪对白视频| 69av精品久久久久久| 欧美+亚洲+日韩+国产| 制服人妻中文乱码| 亚洲七黄色美女视频| 亚洲 国产 在线| 制服人妻中文乱码| 狠狠狠狠99中文字幕| 国产高清视频在线播放一区| 国产精品永久免费网站| 久久人妻福利社区极品人妻图片| 免费在线观看影片大全网站| svipshipincom国产片| 国内毛片毛片毛片毛片毛片| 老司机靠b影院| 国产激情久久老熟女| 黑丝袜美女国产一区| 亚洲精品国产区一区二| 午夜成年电影在线免费观看| 叶爱在线成人免费视频播放| 悠悠久久av| 黄色丝袜av网址大全| 日本免费a在线| 午夜视频精品福利| 级片在线观看| 日本a在线网址| 日韩精品青青久久久久久| 两个人看的免费小视频| 日韩av在线大香蕉| 最好的美女福利视频网| 久久青草综合色| 亚洲精品久久国产高清桃花| 搡老岳熟女国产| 在线国产一区二区在线| 欧美乱妇无乱码| 日本 av在线| 亚洲专区字幕在线| 国产国语露脸激情在线看| 欧美黄色片欧美黄色片|