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

    中甸牦牛mtDNA D-Loop區(qū)遺傳多樣性及系統(tǒng)進(jìn)化分析

    2016-11-24 06:03:27涂世英柴志欣鐘金城
    中國草食動(dòng)物科學(xué) 2016年3期
    關(guān)鍵詞:中甸野牛核苷酸

    涂世英,汪 琦,柴志欣,鐘金城

    (1.云南省動(dòng)物衛(wèi)生監(jiān)督所,昆明650000;2.西南民族大學(xué),動(dòng)物遺傳育種學(xué)國家民委-教育部重點(diǎn)實(shí)驗(yàn)室,成都610041;3.西南民族大學(xué),青藏高原研究院,成都610041)

    中甸牦牛mtDNA D-Loop區(qū)遺傳多樣性及系統(tǒng)進(jìn)化分析

    涂世英1,汪 琦2,3,柴志欣2,3,鐘金城2,3

    (1.云南省動(dòng)物衛(wèi)生監(jiān)督所,昆明650000;2.西南民族大學(xué),動(dòng)物遺傳育種學(xué)國家民委-教育部重點(diǎn)實(shí)驗(yàn)室,成都610041;3.西南民族大學(xué),青藏高原研究院,成都610041)

    根據(jù)牦牛線粒體DNA序列設(shè)計(jì)引物,擴(kuò)增獲得了中甸牦牛線粒體D-loop區(qū)全序列,并以山羊?yàn)橥鈱賹ε喛撇糠峙7N(野牦牛、九龍牦牛、麥洼牦牛、西藏牦牛、天祝牦牛、青海牦牛、大通牦牛、美洲野牛、歐洲野牛、印度野牛、亞洲水牛和普通牛)的mtDNA D-loop序列進(jìn)行了系統(tǒng)進(jìn)化分析,以期了解中甸牦牛的遺傳多樣性,從而為中甸牦牛遺傳資源的保護(hù)、開發(fā)和利用奠定理論基礎(chǔ)。結(jié)果表明:①中甸牦牛線粒體D-loop區(qū)序列長在890~910 bp之間,T、C、A、G四種堿基的平均含量分別為28.78%、24.41%、32.34%和14.47%,存在一定的堿基組成偏倚性;②中甸牦牛共有13種單倍型,平均單倍型多樣性(Hd)為0.983 3,平均核苷酸多樣性(Pi)為0.005 34,遺傳多樣性較為豐富;③中甸牦牛與已知牛種mtDNA D-loop區(qū)序列雙參數(shù)距離顯示,其與麥洼牦牛距離最?。?.006),與亞洲水牛距離最大(0.179);④系統(tǒng)進(jìn)化分析顯示,中甸牦牛是我國眾多家牦牛類群中的一支,與麥洼牦牛聚為一類,推測其可能與麥洼牦牛因地理位置存在基因交流,由共同祖先進(jìn)化而來。

    中甸牦牛;D-loop區(qū);牛亞科;系統(tǒng)發(fā)育

    doi:10.3969/j.issn.2095-3887.2016.03.002

    牦牛是青藏高原特有的珍稀家畜品種,是高寒草地不可替代的生物資源,已成為青藏高原地區(qū)畜牧業(yè)發(fā)展的獨(dú)立板塊。中甸牦牛主要分布在海拔3 000 m以上的高寒地區(qū),是云南迪慶高原地區(qū)特有的畜種,主產(chǎn)于中甸縣尼西、格咱、大中甸、小中甸等藏族聚居的鄉(xiāng)鎮(zhèn),在大理州高寒地區(qū)也有分布,數(shù)量以迪慶中甸縣最多,是我國優(yōu)良牦牛類型之一[1]。中甸牦牛毛色以黑色居多,其次為黑白花,公母都有角,角為黑色或灰白,體格粗壯結(jié)實(shí),具有耐寒耐粗飼,乳肉品質(zhì)優(yōu)良和抗病力強(qiáng)等優(yōu)良種質(zhì)特性[2-3]。由于受中甸地區(qū)社會經(jīng)濟(jì)、自然氣候條件以及公犏牛不育等因素制約,目前仍處在用本地黃牛與牦牛雜交繁殖犏牛的階段,因此將來對中甸牦牛的研究主要集中在適宜高寒自然環(huán)境的生理特性和改良遺傳性狀等方面。

    線粒體DNA(mitochondrial DNA,mtDNA)因具有母系遺傳、多態(tài)性豐富以及結(jié)構(gòu)簡單等特點(diǎn),被認(rèn)為是研究物種起源、演化和分類最好的分子遺傳標(biāo)記[4],其中D-loop區(qū)進(jìn)化速度快,變異較大,適用于對亞科內(nèi)屬、種間的系統(tǒng)學(xué)研究[5]。張成福等[6-7]通過測定西藏牦牛類群的D-loop區(qū)序列對西藏牦牛的遺傳多樣性、類群間親緣關(guān)系和系統(tǒng)進(jìn)化關(guān)系進(jìn)行了研究;李齊發(fā)等[8]通過測定九龍牦牛的D-loop區(qū)全序列并對其進(jìn)行了系統(tǒng)發(fā)育分析。目前尚未見云南中甸牦牛線粒體基因的相關(guān)報(bào)道。本研究利用DNA測序技術(shù),測定了中甸牦牛線粒體基因D-Loop區(qū)全序列,并以羊亞科(Caprinae)中的山羊(Capra hircus)作為外類群,以及部分牛品種(野牦牛、九龍牦牛、麥洼牦牛、西藏牦牛、天祝牦牛、青海牦牛、大通牦牛、美洲野牛、歐洲野牛、印度野牛、亞洲水牛和普通牛)進(jìn)行系統(tǒng)進(jìn)化分析,探究中甸牦牛的遺傳多樣性及其物種間的親緣關(guān)系,以期從分子水平上探討牛亞科不同物種間的系統(tǒng)進(jìn)化關(guān)系及中甸牦牛的分類地位,為中甸牦牛遺傳資源的進(jìn)一步保護(hù)、定向育種、開發(fā)和利用提供理論基礎(chǔ)。

    1 材料與方法

    1.1 樣品采集

    在云南省迪慶中甸縣,選取成年健康的中甸牦牛(15頭),采集耳組織樣品,75%乙醇保存帶回實(shí)驗(yàn)室,存于-20℃的冰箱中備用。

    1.2 基因組DNA的提取及檢測

    運(yùn)用動(dòng)物組織基因組提取試劑盒(TianGen生物技術(shù)公司)提取基因組DNA,用1%的瓊脂糖凝膠電泳和紫外分光光度計(jì)分別檢測DNA的純度和濃度,-20℃保存?zhèn)溆谩?/p>

    1.3 引物設(shè)計(jì)與序列擴(kuò)增

    根據(jù)GenBank中公布的普通牛線粒體DNA序列(AF492351.1),利用Primer Premier 5.0設(shè)計(jì)1對引物擴(kuò)增三江黃牛的D-loop區(qū),引物序列F:CTGCAGTCTCAC CCATCAACC;R:GGGGTGTAGATGCTTGC,由英濰捷基(上海)生物科技有限公司合成。

    PCR反應(yīng)體系為25 μL。其中,上、下游引物各1 μL(10 pmol/μL),DNA模板1 μL(200~600 ng/μL),2×long Taq DNA預(yù)混酶12.50 μL(0.625 U),ddH2O(純化水)9.50 μL。

    PCR擴(kuò)增條件:94℃預(yù)變性2 min;94℃變性15 s,57.7℃退火15 s,72℃延伸30 s,35個(gè)循環(huán);72℃延伸3 min,4℃保存。

    1.4 目的基因克隆測序

    PCR產(chǎn)物用1%瓊脂糖凝膠電泳檢測目的片段大小,再用DNA純化試劑盒(由TianGen公司生產(chǎn))進(jìn)行回收分離純化,直接連接到pMDTM19-T載體(由TaKaRa公司生產(chǎn))上,并轉(zhuǎn)化到高效感受態(tài)細(xì)胞DH5α中,LB平板培養(yǎng)基(含有2%X-Gal、1%IPTG和1%Amp)上37℃培養(yǎng)12~16 h,篩選陽性克隆在700 μLLB液體培養(yǎng)基(含1%Amp)中37℃振蕩培養(yǎng)6~7 h。重組子經(jīng)菌液PCR鑒定后,送往英濰捷基(上海)生物科技有限公司進(jìn)行正反雙向測序。

    1.5 牛亞科部分牛種和外類群構(gòu)建D-loop區(qū)序列的來源

    選取牛亞科中具有代表性的牛種和中甸牦牛進(jìn)行系統(tǒng)建群分析,包括野牦牛(Poephagus mutusc)、九龍牦牛(Jiulong yak)、麥洼牦牛(Maiwa yak)、西藏牦牛(Xizang yak)、天祝牦牛(Tianzhu yak)、青海牦牛(Qinghai yak)、大通牦牛(Datong yak)、美洲野牛(Bison bison)、歐洲野牛(Bison bonasus)、印度野牛(Bison gaurus)、亞洲水牛(Bubalus bubalis)、普通牛(Bos taurus)與中甸牦牛(Zhongdian yak)比對分析,并選取與牛亞科親緣關(guān)系較近的羊亞科中的山羊(Capra hircus)作為系統(tǒng)發(fā)育分析的外類群,從GenBank上下載這些相關(guān)牛種D-loop區(qū)的全序列(見表1)。

    1.6 數(shù)據(jù)分析和系統(tǒng)進(jìn)化分析

    選取擴(kuò)增所得中甸牦牛線粒體DNA D-loop區(qū)序列以及從GenBank下載的所有序列用DNAMAN進(jìn)行序列編輯整理,再用clustalX軟件進(jìn)行序列同源性比對,利用Dnasp 5.0軟件統(tǒng)計(jì)三江黃牛線粒體DNAD-loop區(qū)的核苷酸單倍型數(shù)目(number ofhaplotypes)、單倍型多樣度(haplotype diversity)、核苷酸多樣性(nucleotide diversity)、平均核苷酸差異(average number ofnucleotide differences)等,并對所有突變位點(diǎn)進(jìn)行Tajima中性檢驗(yàn)。用MEGA5.0軟件統(tǒng)計(jì)三江黃牛核苷酸序列的堿基組成、轉(zhuǎn)換/顛換比(transition/transversion,Ts/Tv),并用Kimura2-parameter模型計(jì)算遺傳距離,運(yùn)用鄰接法(NJ)與其他各牛種構(gòu)建系統(tǒng)發(fā)育樹,自舉分析(bootstrap test)1 000次重復(fù)抽樣檢驗(yàn)獲得置信度(BP)。

    表1 牛亞科代表性牛種和外類群的D-loop區(qū)序列

    2 結(jié)果與分析

    2.1 中甸牦牛mtDNA D-loop區(qū)的PCR擴(kuò)增

    中甸牦牛mtDNAD-loop區(qū)PCR擴(kuò)增產(chǎn)物凝膠電泳(圖1),目的片段大小1 000 bp左右,與預(yù)期結(jié)果一致,表明擴(kuò)增成功。

    2.2 中甸牦牛mtDNA D-loop區(qū)的核苷酸組成及多樣性

    通過測序可知,15頭中甸牦牛mtDNA D-loop區(qū)全長均在909~911 bp之間,存在堿基插入、缺失和替換,共有51個(gè)變異位點(diǎn),其中15個(gè)插入或缺失位點(diǎn),36個(gè)轉(zhuǎn)化或顛換位點(diǎn),轉(zhuǎn)換/顛換(Ts/Tv)比1.77,表明中甸牦牛核苷酸多態(tài)性的變異主要以轉(zhuǎn)換為主,其中A→G/C→T轉(zhuǎn)換次數(shù)最多,顛換次數(shù)較少。在轉(zhuǎn)換變異類型中以A/G和T/C為主,在插入類型以C堿基為主。單態(tài)突變位點(diǎn)35個(gè),占多態(tài)位點(diǎn)總數(shù)的68.63%,簡約信息位點(diǎn)16個(gè),占多態(tài)位點(diǎn)總數(shù)的31.37%。T、C、A、G四種堿基的平均含量為28.78%、24.41%、32.34%和14.47%,G+C含量為38.88%,A+T的含量為61.12%,說明中甸牦牛mtDNA D-loop區(qū)富含A/T堿基,具有一定偏倚性。

    圖1 中甸牦牛線粒體DNAD-loop區(qū)擴(kuò)增產(chǎn)物

    用Dnasp5.0軟件對15頭中甸牦牛個(gè)體的單倍型進(jìn)行多樣性分析(表2),結(jié)果表明,中甸牦牛共有13種單倍型,平均單倍型多樣性(Hd)為0.983 3,平均核苷酸多樣性(Pi)為0.005 34。對所有的突變位點(diǎn)進(jìn)行Tajima中性檢驗(yàn),得到Tajim’D值:1.465 04>0.1,表明差異不顯著,所有的變異均符合中性進(jìn)化。

    表2 中甸牦牛mtDNAD-loop區(qū)的單倍型分布

    2.3 中甸牦牛與各牛種類群間D-loop區(qū)遺傳多樣性差異

    除印度野牛外,所列各牛種中4種堿基分布情況大致相同(表3),T、C、A、G四種核苷酸的平均含量分別為 28.66%(28.36%~29.18%)、25.11%(24.40%~25.51%)、32.36%(31.69%~32.85%)、13.87%(13.34%~14.47%),其中A/T含量(61.02%)明顯高于G/C含量(38.98%),表明各牛種線粒體DNA D-loop區(qū)均富含堿基A和T。

    牛亞科內(nèi)部分牛種類群間D-loop區(qū)序列差異百分比見表4。由表4可知,中甸牦牛與美洲野牛的序列差異最小,為1.96%,與其他幾種家牦牛的序列差異百分比在5.42%~7.99%之間,與印度野牛差異最大,為25.70%。在牛亞科中,九龍牦牛和大通牦牛序列無差異,麥洼牦牛和天祝牦牛序列也無差異;麥洼牦牛、西藏牦牛、天祝牦牛與青海牦牛的序列間差異較小,為0.11%,而與印度野牛、亞洲水牛的序列差異百分比較大(23.29%~33.58%);這與以往的研究結(jié)果基本一致。

    表3 牛亞科部分牛種和外類群的D-loop區(qū)段的堿基組成%

    表4 牛亞科內(nèi)不同物種間線粒體DNAD-loop區(qū)序列差異百分比%

    2.4 中甸牦牛與各牛種類群間的遺傳距離及系統(tǒng)進(jìn)化分析

    根據(jù)中甸牦牛和各牛種類群間線粒體DNA D-loop區(qū)堿基序列,用Mega 5.0軟件計(jì)算14個(gè)牛種間的遺傳距離(表5)。結(jié)果表明,14個(gè)牛種間的遺傳距離范圍在0~0.192之間,其中中甸牦牛和亞洲水牛的距離最大(0.179),與麥洼牦牛遺傳距離最短(0.006),這表明它們的親緣關(guān)系較近,麥洼牦牛是主要分布于四川省的優(yōu)良草地型地方品種,而中甸牦牛常年在云南高海拔草場上放牧生存,兩者均屬于家牦牛,推測因地理位置接壤可能存在基因交流,親緣關(guān)系上可能由共同的祖先進(jìn)化而來。在其他牛種中,歐洲野牛與普通牛的遺傳距離最遠(yuǎn)(0.186),天祝牦牛與青海牦牛的遺傳距離最短(0.001)。

    根據(jù)各牛種間的遺傳距離,以山羊作為外類群,用鄰接法(NJ)構(gòu)建分子系統(tǒng)進(jìn)化樹(圖2)。由圖2可知,牛亞科和外類群(山羊)分別各自聚為一類,牛亞科內(nèi),亞洲水牛與其他物種各自聚為一類;中甸牦牛和麥洼牦牛聚在一起,再與九龍牦牛、大通牦牛、天祝牦牛和青海牦牛相聚,同屬于家牦牛的分支。中甸牦牛分布于云南省,而麥洼牦牛分布于四川省,兩地的距離較近,由此推測,兩者存在地理上的親緣關(guān)系或是基因交流。此外,野牦牛和西藏牦牛聚為一類,歐洲野牛和普通牛聚為一類,聚類結(jié)果與以往的動(dòng)物分類學(xué)結(jié)果基本一致。

    表5 牛亞科各牛種間Kimura雙參數(shù)遺傳距離

    圖2 基于Kimura雙參數(shù)遺傳距離的各牛種類群間NJ系統(tǒng)進(jìn)化關(guān)系

    3 討論

    遺傳多樣性是所有生物多樣性的重要組成部分,它是生物進(jìn)化和物種分化的基礎(chǔ)[9],也對種群維持繁衍和適應(yīng)生境變化具有重要意義。由于Cytb、l2S rRNA等線粒體編碼基因的變異過程較慢,常用來進(jìn)行種以上等級的系統(tǒng)發(fā)育關(guān)系重建[10],而mtDNA D-Loop區(qū)的序列變異則較快,因此被用來進(jìn)行種以下水平的系統(tǒng)關(guān)系研究。通常動(dòng)物mtDNA進(jìn)化主要以核苷酸轉(zhuǎn)換為主,插入和缺失較少,核苷酸的轉(zhuǎn)換是插入和缺失的2~10倍,而在mtDNAD-Loop區(qū)核苷酸的轉(zhuǎn)換概率更是遠(yuǎn)大于核苷酸插入和缺失的概率。家牦牛豐富的遺傳多樣性對我國乃至全世界牦牛產(chǎn)業(yè)的可持續(xù)發(fā)展起著舉足輕重的作用。

    目前對家牦牛的mtDNA D-loop區(qū)研究主要集中在四川牦牛和西藏牦牛上。賴松家等[5]運(yùn)用測序技術(shù)測定了中國5個(gè)家牦牛品種的mtDNA D-loop區(qū)全序列,表明我國牦牛mtDNA D-loop單倍型類型豐富并有2個(gè)母系來源,分別是橫斷高山型和青藏高原型兩大生態(tài)類型;鐘金城等[11]和肖玉萍等[12]采用RAPD和AFLP標(biāo)記對九龍牦牛、麥洼牦牛、大通牦牛和天祝白牦牛進(jìn)行了分析,結(jié)果表明九龍牦牛的遺傳多樣性最為豐富,天祝白牦牛遺傳多樣性最低,4個(gè)牦牛品種聚為兩大類:九龍牦牛為一類,其他3個(gè)牦牛品種聚為一類,這與本研究聚類分析所得結(jié)論基本一致。鐘金城等[13]和柴志欣等[14]用RAPD引物對西藏11個(gè)牦牛群體進(jìn)行研究,表明西藏牦牛類群內(nèi)部具有較多遺傳變異,推測西藏西部可能是牦牛的起源地之一。

    本研究以動(dòng)物mtDNAD-loop區(qū)為分子標(biāo)記,對中甸牦牛mtDNAD-loop區(qū)進(jìn)行序列分析,結(jié)果顯示,15頭中甸牦牛mtDNA D-loop區(qū)序列長度在890~910 bp之間,個(gè)體間差異較小。通過對中甸牦牛mtDNA D-loop區(qū)序列的比對分析,發(fā)現(xiàn)T、C、A、G四種堿基的平均含量為28.78%、24.41%、32.34%和14.47%,堿基組成偏好性明顯。中甸牦牛平均單倍型多樣性(Hd)為0.983 3,平均核苷酸多樣性(Pi)為0.005 34,序列變異存在堿基插入、缺失和替換,共有51個(gè)變異位點(diǎn),其中15個(gè)插入或缺失位點(diǎn),36個(gè)轉(zhuǎn)化或顛換位點(diǎn),轉(zhuǎn)換/顛換(Ts/Tv)比1.77,表明中甸牦牛種內(nèi)遺傳多樣性豐富,這與田應(yīng)華等[15]利用微衛(wèi)星標(biāo)記技術(shù)對中甸牦牛分析的結(jié)果一致,揭示中甸牦牛群體變異在中國牦牛地方群體中處于較高水平。

    本研究對14個(gè)牛亞科部分牛種進(jìn)行遺傳距離和系統(tǒng)進(jìn)化分析,顯示中甸牦牛是我國家牦牛中的一個(gè)品種,與麥洼牦牛遺傳距離最近,為0.006,說明兩者關(guān)系較近,中甸牦牛產(chǎn)自云南省,而麥洼牦牛產(chǎn)自四川省,兩地的距離較近,并且中甸牦牛和麥洼牦牛同屬草地型牦牛,生活在高山草甸草場,因此推測兩者可能存在基因交流。所構(gòu)建的系統(tǒng)發(fā)育樹顯示,牛亞科中水牛屬最早分化出來,然后是印度野牛、歐洲野牛與普通牛、美洲野牛、家牦牛與野牦牛的分支;聚類分析結(jié)果表明,歐洲野牛與普通牛的親緣關(guān)系較近,另外美洲野牛與牦牛之間親緣關(guān)系較近,這與郭松長等[10]所得出的結(jié)果一致;系統(tǒng)進(jìn)化樹還顯示,西藏牦牛和野牦牛是一個(gè)分支,這也和以往得出的的研究成果相同,現(xiàn)在的家養(yǎng)牦牛起源于我國的西藏,野牦牛則是家養(yǎng)牦牛的近緣野生種[16]。以牦?;鳛榭甲C,野牦牛與家牦牛都曾分布在歐亞大陸東北部,是原始牦牛的后代,由于地殼運(yùn)動(dòng)、氣候變化等原因遷移至我國青藏高原地區(qū),并適應(yīng)了高寒、缺氧等生態(tài)環(huán)境從而生存下來,因此今天的家牦牛和野牦牛都由同一祖先進(jìn)化而來,它們之間不存在祖先和后代的關(guān)系[17]。

    目前對中甸牦牛的遺傳多樣性研究尚未見報(bào)道,本研究首次對中甸牦牛的遺傳多樣性和系統(tǒng)進(jìn)化作了系統(tǒng)全面的研究。中甸牦牛是我國優(yōu)良地方家牦牛品種之一,是我國眾多牦牛類群中的一支,其群體遺傳變異程度較高,遺傳多樣性較豐富,該結(jié)果可為中甸牦牛遺傳資源的進(jìn)一步保護(hù)、定向育種、開發(fā)和利用提供理論基礎(chǔ)。

    [1] 和紹禹,田允波,葛長榮,等.中甸牦牛[J].黃牛雜志,1998,24(2):22-24.

    [2] 云南省畜牧局.云南省家畜家禽品種志[M].昆明:云南科技出版社,1987:88-94.

    [3] 郭憲,閻萍,曾玉峰,等.中國牦牛遺傳資源現(xiàn)狀分析[J].中國畜禽種業(yè),2008(1):60-62.

    [4] Bruford MW,BradleyD G,Luikart G.DNA markers reveal the complexity of livestock domestication[J].Nature Reviews Genetics,2003,4(11):900-910.

    [5] 賴松家,王玲,劉益平,等.中國部分牦牛品種線粒體DNA遺傳多態(tài)性研究[J].遺傳學(xué)報(bào),2005,32(5):463-470.

    [6] 張成福,徐麗娟,姬秋梅,等.西藏牦牛mtDNA D-loop區(qū)的遺傳多樣性及遺傳分化[J].生態(tài)學(xué)報(bào),2012,5(32):1387-1395.

    [7] 宋喬喬,鐘金城,張成福,等.西藏牦牛線粒體DNA的遺傳多樣性及系統(tǒng)進(jìn)化分析[J].獸類學(xué)報(bào),2014,34(4):356-365.

    [8] 李齊發(fā),李隱俠,趙興波,等.牦牛線粒體DNAD-loop區(qū)序列測定及其在牛亞科中分類地位的研究[J].畜牧獸醫(yī)學(xué)報(bào),2008,39(1):1-6.

    [9] 馬志杰,鐘金城,韓建林,等.野牦牛mtDNAD-Loop區(qū)的遺傳多樣性[J].生態(tài)學(xué)報(bào),2009,29(9):4798-4802.

    [10]郭松長,劉建全,祁得林,等.牦牛的分類學(xué)地位及起源研究:mt DNAD-Loop序列的分析[J].獸類學(xué)報(bào),2006,26(4):325-330.

    [11]鐘金城,陳智華,趙素君,等.牦牛生態(tài)類型的分類[J].生態(tài)學(xué)報(bào),2006,26(7):2068-2072.

    [12]肖玉萍,鐘金城,金雙.4個(gè)牦牛品種的RAPD遺傳多樣性研究[J].中國牛業(yè)科學(xué),2007,33(6):5-10.

    [13]鐘金城,柴志欣,姬秋梅,等.西藏牦牛的遺傳多樣性及其系統(tǒng)進(jìn)化研究[J].西南民族大學(xué)學(xué)報(bào):自然科學(xué)版,2011,37(3):368-378.

    [14]柴志欣,趙上娟,姬秋梅,等.西藏牦牛的RAPD遺傳多樣性及其分類研究[J].畜牧獸醫(yī)學(xué)報(bào),2011,42(10):1380-1386.

    [15]田應(yīng)華,錢林東,凌軍,等.云南中甸牦牛的遺傳多樣性研究[J].西南農(nóng)業(yè)學(xué)報(bào),2009,22(3):794-797.

    [16]鐘金城,柴志欣,馬志杰,等.野牦牛線粒體基因組序列測定及其系統(tǒng)進(jìn)化[J].生態(tài)學(xué)報(bào),2015,35(5):1564-1572.

    [17]鐘金城.牦牛遺傳與育種[M].成都:四川科學(xué)技術(shù)出版社,1996:24-68.

    Analysis on Genetic Diversity and Phyletic Evolution of Mitochondrial DNA D-loop Region from Zhongdian Yak

    Tu Shiying1,WangQi2,3,ZhongJincheng2,3,et al
    (1.The Animal Health Supervision,Yunnan Province,Kunming650000,China;2.KeyLaboratoryofAnimal Genetics and BreedingofState Ethnic Affair Commission and MinistryofEduction,Chengdu 610041,China;3.Qinghai-Tibet Plateau College,Southwest Universityfor Nationalities,Chengdu 610041,China)

    The Zhongdian yak mitochondrial D-loop area was cloned using a specific primers designed according to the mitochondrial DNA sequence of yak,which was then used for phylogenetic analysis combined with the mtDNA D-loop sequence of bovine typical cattle breeds,including Poephagus mutusc,Jiulong yak,Maiwa yak,Xizang yak,Tianzhu yak,Qinghai yak,Datong yak,Bison bison,Bison bonasus,Bison gaurus,Bubalus bubalis and Bos Taurus,with that of sheep as the exogenous reference.These data may be beneficial for our understanding the genetic diversity of Zhongdian yak,and put forward the theoretical foundation for the protection,development and utilization of genetic resources of Zhongdian yak.The results showed that,①Zhongdian yak mitochondrial D-loop area sequence length between 890~910 bp.T,C,A,G took the average content of 28.78%,24.41%,32.34%and 14.47%,respectively,and there was a certain base composition of bias;②a total of 13 Zhongdian yak haploid type with the average haploid typediversity(Hd)of 0.983 3,the average nucleotide diversity(Pi)of 0.005 34 was identified;③Zhongdian yak shared the minimum sequence distance with that ofMaiwa yak(0.000 6),and hit the maximum sequence distance on Asian buffalo(0.179),which was used the Kimura two-parameter genetic distance with known cattle mtDNA D-loop region sequences;④phylogenetic evolution analysis found that Zhongdian yak is an independent yak breeds,obtained the Maiwa yak into one group,and speculated that yak presence ofgene flowdue totheir location,it was evolved froma common ancestor.

    Zhongdian yak;D-loop region;bovine;phylogeny

    S823.8+5

    A

    2095-3887(2016)03-0005-06

    2016-03-24

    四川省科技計(jì)劃項(xiàng)目(2015JY0248)

    涂世英(1968-),女,高級畜牧師。研究方向:牦牛生產(chǎn)學(xué)。

    鐘金城(1963-),男,教授,博士。研究方向:動(dòng)物遺傳學(xué)。

    猜你喜歡
    中甸野牛核苷酸
    單核苷酸多態(tài)性與中醫(yī)證候相關(guān)性研究進(jìn)展
    徐長風(fēng):核苷酸類似物的副作用
    肝博士(2022年3期)2022-06-30 02:48:28
    “中甸之爭”:清初麗江改土歸流原因探析
    西藏研究(2021年1期)2021-06-09 08:09:34
    美洲野牛當(dāng)寵物
    Acknowledgment to reviewers—November 2018 to September 2019
    中甸牦牛養(yǎng)殖技術(shù)要點(diǎn)
    小學(xué)生學(xué)習(xí)指導(dǎo)(小軍迷聯(lián)盟)(2018年10期)2018-10-12 01:13:42
    中甸牦牛消化道寄生蟲初步調(diào)查及防制
    兩頭野牛
    廣東人群8q24rs1530300單核苷酸多態(tài)性與非綜合征性唇腭裂的相關(guān)性研究
    人人妻人人看人人澡| 色哟哟·www| 国产伦精品一区二区三区视频9| 精品一区二区免费观看| av视频在线观看入口| 国产精品一区www在线观看 | 亚洲精品在线观看二区| 久久久久久大精品| 亚洲va在线va天堂va国产| 中文字幕免费在线视频6| 国产精品不卡视频一区二区| 一进一出抽搐gif免费好疼| 日日撸夜夜添| 一本久久中文字幕| av在线老鸭窝| 小说图片视频综合网站| 香蕉av资源在线| av视频在线观看入口| 麻豆av噜噜一区二区三区| 欧美潮喷喷水| 日本黄大片高清| 日日撸夜夜添| www.www免费av| 精品久久久久久成人av| 国产白丝娇喘喷水9色精品| 日韩欧美精品v在线| 欧美一区二区国产精品久久精品| 国产成人影院久久av| 国产精品久久电影中文字幕| 亚洲最大成人中文| 欧美日韩精品成人综合77777| 亚洲狠狠婷婷综合久久图片| 简卡轻食公司| 国产精品一区二区免费欧美| 日本成人三级电影网站| www.色视频.com| 亚洲av第一区精品v没综合| 国内精品久久久久久久电影| 亚洲中文日韩欧美视频| 国产精品福利在线免费观看| 我要看日韩黄色一级片| 日本黄色片子视频| 久久精品国产亚洲av涩爱 | 99视频精品全部免费 在线| 99国产极品粉嫩在线观看| 久久久久久伊人网av| 日韩av在线大香蕉| 日韩欧美在线二视频| 一个人看视频在线观看www免费| 日日啪夜夜撸| 色综合站精品国产| 午夜激情欧美在线| 欧美高清成人免费视频www| ponron亚洲| 国产高清视频在线观看网站| 毛片女人毛片| 99热精品在线国产| 嫩草影视91久久| 有码 亚洲区| 国产欧美日韩一区二区精品| 亚洲欧美日韩高清在线视频| 国产极品精品免费视频能看的| 久久99热这里只有精品18| 女人十人毛片免费观看3o分钟| 麻豆av噜噜一区二区三区| 在线观看舔阴道视频| 国产乱人视频| 我要看日韩黄色一级片| 91精品国产九色| 在线播放无遮挡| 此物有八面人人有两片| 久久久久久大精品| 一区二区三区免费毛片| 人妻久久中文字幕网| 欧美日韩瑟瑟在线播放| 免费av毛片视频| 色播亚洲综合网| 中文字幕免费在线视频6| 少妇熟女aⅴ在线视频| 欧美日本亚洲视频在线播放| 一级毛片久久久久久久久女| 女人被狂操c到高潮| 亚洲无线观看免费| 黄片wwwwww| 欧美性猛交黑人性爽| 亚洲一级一片aⅴ在线观看| 日日夜夜操网爽| 一本精品99久久精品77| 国产色爽女视频免费观看| 国产精品电影一区二区三区| 亚洲性久久影院| 久久久久国产精品人妻aⅴ院| 国产在线男女| 又黄又爽又免费观看的视频| 三级国产精品欧美在线观看| 久久精品人妻少妇| 精品久久久久久久久亚洲 | 久久国内精品自在自线图片| 国内精品宾馆在线| 日本黄色片子视频| 色尼玛亚洲综合影院| 波多野结衣巨乳人妻| 欧美区成人在线视频| 国产精品福利在线免费观看| 国产精品国产三级国产av玫瑰| 91久久精品国产一区二区三区| 日本爱情动作片www.在线观看 | 免费大片18禁| 国产高清有码在线观看视频| 亚洲va在线va天堂va国产| 日本黄色片子视频| 久久精品国产99精品国产亚洲性色| 欧美极品一区二区三区四区| 国产久久久一区二区三区| 成人美女网站在线观看视频| 久久久久免费精品人妻一区二区| 亚洲经典国产精华液单| 欧美一区二区精品小视频在线| 亚洲成人中文字幕在线播放| 免费高清视频大片| 亚洲aⅴ乱码一区二区在线播放| 日韩 亚洲 欧美在线| 九九爱精品视频在线观看| 亚洲中文日韩欧美视频| 波野结衣二区三区在线| 日韩人妻高清精品专区| 亚洲欧美日韩无卡精品| 亚洲精品粉嫩美女一区| a在线观看视频网站| 麻豆精品久久久久久蜜桃| 最近最新中文字幕大全电影3| 日韩,欧美,国产一区二区三区 | 色尼玛亚洲综合影院| 国产精品久久久久久久电影| 少妇的逼好多水| 真人一进一出gif抽搐免费| 成人一区二区视频在线观看| 2021天堂中文幕一二区在线观| 日本爱情动作片www.在线观看 | 一级黄色大片毛片| 九九热线精品视视频播放| 日本黄色片子视频| 国产精品久久久久久久电影| 色综合站精品国产| 别揉我奶头 嗯啊视频| 国产亚洲精品久久久久久毛片| 深夜精品福利| 精品无人区乱码1区二区| 伦精品一区二区三区| 午夜福利欧美成人| 亚洲国产色片| 成人高潮视频无遮挡免费网站| 一区福利在线观看| 别揉我奶头~嗯~啊~动态视频| 看黄色毛片网站| 深夜a级毛片| 欧美成人a在线观看| 国产色婷婷99| 自拍偷自拍亚洲精品老妇| 成人午夜高清在线视频| av专区在线播放| 婷婷精品国产亚洲av在线| 69av精品久久久久久| 精品久久久久久久久久免费视频| 琪琪午夜伦伦电影理论片6080| 淫妇啪啪啪对白视频| 亚洲午夜理论影院| 午夜免费成人在线视频| a级毛片免费高清观看在线播放| 欧美激情在线99| 国产麻豆成人av免费视频| 亚洲第一电影网av| 中国美白少妇内射xxxbb| 两个人的视频大全免费| 99久久中文字幕三级久久日本| 舔av片在线| 国产精品美女特级片免费视频播放器| 久久精品国产鲁丝片午夜精品 | 久久久久久久久大av| 91在线精品国自产拍蜜月| 夜夜夜夜夜久久久久| 赤兔流量卡办理| 老熟妇仑乱视频hdxx| 国产高清视频在线观看网站| 亚洲va在线va天堂va国产| 色尼玛亚洲综合影院| 亚洲av二区三区四区| 国产一区二区在线观看日韩| 在线观看美女被高潮喷水网站| 国内精品久久久久久久电影| 色综合婷婷激情| 最新中文字幕久久久久| 国产美女午夜福利| 熟女电影av网| 亚洲无线在线观看| 观看免费一级毛片| 欧美bdsm另类| 热99re8久久精品国产| 亚洲国产精品合色在线| АⅤ资源中文在线天堂| 亚洲久久久久久中文字幕| 亚洲第一区二区三区不卡| 美女xxoo啪啪120秒动态图| 一卡2卡三卡四卡精品乱码亚洲| 在线天堂最新版资源| 99九九线精品视频在线观看视频| 国内久久婷婷六月综合欲色啪| 国产午夜福利久久久久久| 国产精品精品国产色婷婷| 欧美国产日韩亚洲一区| 人人妻,人人澡人人爽秒播| 一区二区三区免费毛片| 亚洲国产日韩欧美精品在线观看| 日本 av在线| 亚洲av免费高清在线观看| 亚洲,欧美,日韩| 丰满的人妻完整版| 欧美成人免费av一区二区三区| 全区人妻精品视频| 欧美激情在线99| 亚洲三级黄色毛片| 最新中文字幕久久久久| 麻豆av噜噜一区二区三区| 在线天堂最新版资源| 最近视频中文字幕2019在线8| 欧美国产日韩亚洲一区| 国产精品无大码| 一本久久中文字幕| 国产黄片美女视频| 欧美一区二区亚洲| 又黄又爽又刺激的免费视频.| 亚洲精品一区av在线观看| 日韩av在线大香蕉| 日韩在线高清观看一区二区三区 | 老司机午夜福利在线观看视频| 中文亚洲av片在线观看爽| 亚洲最大成人av| av黄色大香蕉| 精品人妻视频免费看| 成年人黄色毛片网站| 日韩,欧美,国产一区二区三区 | 亚洲,欧美,日韩| 88av欧美| 国产精品久久电影中文字幕| 国产av在哪里看| 一本久久中文字幕| av福利片在线观看| 69人妻影院| 欧美最黄视频在线播放免费| 国产高清激情床上av| 国产亚洲av嫩草精品影院| 日日夜夜操网爽| 亚洲中文字幕日韩| 人妻制服诱惑在线中文字幕| 一a级毛片在线观看| 久久九九热精品免费| 在线播放无遮挡| 亚洲国产精品sss在线观看| 97热精品久久久久久| 中文字幕人妻熟人妻熟丝袜美| 国产精品久久久久久久电影| 国产一区二区亚洲精品在线观看| xxxwww97欧美| 国产白丝娇喘喷水9色精品| 日本精品一区二区三区蜜桃| 桃色一区二区三区在线观看| 欧美高清成人免费视频www| 亚洲精品在线观看二区| 国产黄色小视频在线观看| 亚洲一区二区三区色噜噜| 久久亚洲精品不卡| 亚洲va在线va天堂va国产| 欧美高清性xxxxhd video| 欧美日韩瑟瑟在线播放| 国产主播在线观看一区二区| 亚洲av熟女| 一卡2卡三卡四卡精品乱码亚洲| 国产黄片美女视频| 国产高清激情床上av| 亚洲欧美日韩卡通动漫| 91精品国产九色| 老熟妇乱子伦视频在线观看| av国产免费在线观看| 天美传媒精品一区二区| 中出人妻视频一区二区| 又黄又爽又免费观看的视频| 精品国产三级普通话版| 亚洲精华国产精华精| 天堂影院成人在线观看| 日本黄大片高清| 欧美一区二区国产精品久久精品| 国内精品久久久久久久电影| 国内精品一区二区在线观看| 亚洲美女视频黄频| 免费观看的影片在线观看| 国产激情偷乱视频一区二区| 嫩草影视91久久| 亚洲va在线va天堂va国产| 欧美精品啪啪一区二区三区| 99九九线精品视频在线观看视频| 在线国产一区二区在线| 亚洲内射少妇av| 老师上课跳d突然被开到最大视频| 欧美日韩精品成人综合77777| 在线免费观看不下载黄p国产 | 日本a在线网址| 波多野结衣巨乳人妻| 亚洲熟妇中文字幕五十中出| 级片在线观看| 日韩欧美一区二区三区在线观看| 别揉我奶头 嗯啊视频| 少妇人妻一区二区三区视频| 夜夜夜夜夜久久久久| 亚洲 国产 在线| 亚洲四区av| 日本成人三级电影网站| 亚洲国产精品成人综合色| 在线播放无遮挡| 精品久久久久久久久av| 免费在线观看日本一区| 91麻豆精品激情在线观看国产| 人妻久久中文字幕网| 国产精品无大码| 99热6这里只有精品| 12—13女人毛片做爰片一| 最后的刺客免费高清国语| 看十八女毛片水多多多| 国产精品一区二区三区四区久久| 别揉我奶头~嗯~啊~动态视频| 欧美成人性av电影在线观看| 嫁个100分男人电影在线观看| 啦啦啦啦在线视频资源| 久久久久国产精品人妻aⅴ院| 免费av毛片视频| 别揉我奶头~嗯~啊~动态视频| 精品日产1卡2卡| 国内精品久久久久精免费| 搞女人的毛片| 不卡一级毛片| 一区福利在线观看| 1024手机看黄色片| 欧美日本视频| 日本色播在线视频| 免费人成视频x8x8入口观看| 人妻丰满熟妇av一区二区三区| 女同久久另类99精品国产91| 日本与韩国留学比较| 麻豆国产97在线/欧美| 91久久精品国产一区二区成人| 午夜精品在线福利| 日本在线视频免费播放| 成人国产综合亚洲| 天堂动漫精品| 久9热在线精品视频| 免费在线观看成人毛片| 如何舔出高潮| 尤物成人国产欧美一区二区三区| 真实男女啪啪啪动态图| 国产精品电影一区二区三区| 91狼人影院| 欧美区成人在线视频| 日本免费一区二区三区高清不卡| 神马国产精品三级电影在线观看| 最好的美女福利视频网| 日本色播在线视频| 国产免费av片在线观看野外av| 亚洲一级一片aⅴ在线观看| 变态另类丝袜制服| 狂野欧美激情性xxxx在线观看| 一级毛片久久久久久久久女| 亚洲av第一区精品v没综合| 日本 欧美在线| 国产精品嫩草影院av在线观看 | 久久亚洲真实| 久久国内精品自在自线图片| 国产淫片久久久久久久久| 少妇人妻精品综合一区二区 | 亚洲av不卡在线观看| 欧美国产日韩亚洲一区| 人妻丰满熟妇av一区二区三区| 小蜜桃在线观看免费完整版高清| 色综合婷婷激情| 欧洲精品卡2卡3卡4卡5卡区| 成人亚洲精品av一区二区| bbb黄色大片| 精品午夜福利视频在线观看一区| 国内久久婷婷六月综合欲色啪| 久久久精品大字幕| 午夜视频国产福利| 欧美一区二区精品小视频在线| 日本爱情动作片www.在线观看 | 欧美高清性xxxxhd video| 真人做人爱边吃奶动态| 又黄又爽又免费观看的视频| 国产极品精品免费视频能看的| 男人舔奶头视频| 男女下面进入的视频免费午夜| 国产男人的电影天堂91| 制服丝袜大香蕉在线| 男插女下体视频免费在线播放| 国产中年淑女户外野战色| 亚洲美女搞黄在线观看 | 欧美丝袜亚洲另类 | 99国产极品粉嫩在线观看| 男女下面进入的视频免费午夜| 国产男人的电影天堂91| 国产伦人伦偷精品视频| 99久久精品国产国产毛片| 日本黄色片子视频| 国产成人a区在线观看| 久久精品国产亚洲av涩爱 | 国产一区二区在线av高清观看| 嫩草影院新地址| 亚洲欧美日韩高清专用| 成熟少妇高潮喷水视频| 免费在线观看日本一区| 国产精品国产高清国产av| 午夜福利高清视频| 久久婷婷人人爽人人干人人爱| 女生性感内裤真人,穿戴方法视频| 丝袜美腿在线中文| 少妇人妻精品综合一区二区 | 中文字幕人妻熟人妻熟丝袜美| 91在线观看av| 永久网站在线| 国产不卡一卡二| 久久久久久大精品| 国产成人影院久久av| 又爽又黄无遮挡网站| 精品国内亚洲2022精品成人| 成年版毛片免费区| 日韩欧美国产一区二区入口| 色吧在线观看| 搡老妇女老女人老熟妇| 久久九九热精品免费| 校园春色视频在线观看| 国产精品人妻久久久久久| 亚洲第一区二区三区不卡| 色噜噜av男人的天堂激情| 国产一区二区激情短视频| 99riav亚洲国产免费| 免费看日本二区| 麻豆一二三区av精品| 熟女电影av网| 日韩精品有码人妻一区| 狂野欧美激情性xxxx在线观看| 国产精品日韩av在线免费观看| 两个人视频免费观看高清| 一进一出抽搐动态| 久久国产精品人妻蜜桃| 成年版毛片免费区| 国产精品无大码| 国产精品爽爽va在线观看网站| 中文字幕人妻熟人妻熟丝袜美| 成人三级黄色视频| 春色校园在线视频观看| 天堂av国产一区二区熟女人妻| 日韩精品青青久久久久久| 偷拍熟女少妇极品色| 如何舔出高潮| 日本黄色视频三级网站网址| 亚洲av第一区精品v没综合| 久久精品国产99精品国产亚洲性色| 日韩欧美在线乱码| 又粗又爽又猛毛片免费看| АⅤ资源中文在线天堂| 国内精品宾馆在线| 真实男女啪啪啪动态图| 国产视频一区二区在线看| 亚洲美女视频黄频| 99久久无色码亚洲精品果冻| 99热这里只有精品一区| 成人三级黄色视频| 亚洲国产精品成人综合色| 中出人妻视频一区二区| 国产午夜精品论理片| 人人妻人人澡欧美一区二区| av在线天堂中文字幕| 两人在一起打扑克的视频| 尾随美女入室| 亚洲乱码一区二区免费版| 舔av片在线| 亚洲在线观看片| 久久99热这里只有精品18| 狂野欧美白嫩少妇大欣赏| 婷婷精品国产亚洲av| 欧美一区二区国产精品久久精品| 成年女人毛片免费观看观看9| 在线观看66精品国产| 欧美性猛交黑人性爽| 欧洲精品卡2卡3卡4卡5卡区| av天堂中文字幕网| 91久久精品国产一区二区成人| 人妻久久中文字幕网| 国产黄片美女视频| 97人妻精品一区二区三区麻豆| 国产高清有码在线观看视频| 欧美极品一区二区三区四区| 看十八女毛片水多多多| 国产蜜桃级精品一区二区三区| www.www免费av| 国产成人a区在线观看| 免费在线观看影片大全网站| 亚洲熟妇熟女久久| 亚洲成人久久爱视频| 亚洲成av人片在线播放无| 99国产极品粉嫩在线观看| 一级毛片久久久久久久久女| 亚洲欧美精品综合久久99| 欧美日韩瑟瑟在线播放| 草草在线视频免费看| 男人的好看免费观看在线视频| 午夜精品一区二区三区免费看| 三级国产精品欧美在线观看| 亚洲精品粉嫩美女一区| 一级黄片播放器| 不卡视频在线观看欧美| 美女免费视频网站| 国产精品久久视频播放| 欧美一级a爱片免费观看看| 日韩欧美 国产精品| 久久久久久大精品| 免费看日本二区| 国产乱人伦免费视频| 国产精品国产三级国产av玫瑰| 日韩亚洲欧美综合| 精品久久久久久久人妻蜜臀av| 亚洲精品在线观看二区| av视频在线观看入口| 亚洲图色成人| 精品一区二区三区视频在线观看免费| 99在线人妻在线中文字幕| 一级黄片播放器| 黄色一级大片看看| 亚洲av免费在线观看| 亚洲成人久久爱视频| 老司机深夜福利视频在线观看| 香蕉av资源在线| 18禁黄网站禁片午夜丰满| 欧美绝顶高潮抽搐喷水| 日韩欧美三级三区| 午夜激情福利司机影院| 国产精品国产高清国产av| 深夜精品福利| 又爽又黄无遮挡网站| 中文亚洲av片在线观看爽| 综合色av麻豆| 极品教师在线视频| 色哟哟哟哟哟哟| 美女xxoo啪啪120秒动态图| 99久久成人亚洲精品观看| 国产三级中文精品| 啪啪无遮挡十八禁网站| 国产在线男女| 99久久精品热视频| 99热6这里只有精品| 在线观看免费视频日本深夜| 黄色视频,在线免费观看| 在线播放国产精品三级| 国产男靠女视频免费网站| 嫁个100分男人电影在线观看| 91在线观看av| 国产成年人精品一区二区| 亚洲精品456在线播放app | 最近视频中文字幕2019在线8| 老师上课跳d突然被开到最大视频| 亚洲精品一区av在线观看| 日本一二三区视频观看| 搞女人的毛片| 国产一区二区激情短视频| 久久精品国产亚洲av天美| 成人国产一区最新在线观看| 国模一区二区三区四区视频| 日日摸夜夜添夜夜添av毛片 | 97超级碰碰碰精品色视频在线观看| 一级黄色大片毛片| 制服丝袜大香蕉在线| 日日夜夜操网爽| 丰满的人妻完整版| 一级av片app| 波多野结衣高清无吗| 亚洲精品久久国产高清桃花| 在线播放无遮挡| 五月玫瑰六月丁香| 亚洲图色成人| 亚洲成人中文字幕在线播放| 狂野欧美白嫩少妇大欣赏| 欧美三级亚洲精品| 毛片一级片免费看久久久久 | 亚洲经典国产精华液单| 欧美日韩乱码在线| 美女被艹到高潮喷水动态| 久久中文看片网| 精品久久久久久成人av| 午夜免费成人在线视频| 1024手机看黄色片| 亚洲中文字幕日韩| 少妇猛男粗大的猛烈进出视频 | 可以在线观看的亚洲视频| 全区人妻精品视频| 色播亚洲综合网| 国产精品99久久久久久久久| 久久久久久久午夜电影| 国内精品久久久久精免费| 国产精品av视频在线免费观看| xxxwww97欧美| 亚洲专区中文字幕在线| 日韩亚洲欧美综合| 国产视频一区二区在线看| 啪啪无遮挡十八禁网站| 少妇熟女aⅴ在线视频| 国产精品女同一区二区软件 | 日本 av在线| 精品午夜福利视频在线观看一区| 国产三级在线视频|