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

    中間偃麥草全基因組的選擇

    2022-10-17 08:19:52鵬,
    江蘇農業(yè)科學 2022年18期
    關鍵詞:麥草準確度表型

    郭 鵬, 曹 晟

    (天津農學院計算機與信息工程學院,天津 300384)

    中間偃麥草[(Host) Nevski]原產于東歐,是一種根系發(fā)達、抗寒性強、耐旱、耐鹽、適應性廣、再生性好的優(yōu)良牧草。作為多年生牧草,中間偃麥草產量高,已經在畜牧業(yè)發(fā)達的國家得到了廣泛應用。中間偃麥草的果實可食用,其谷物制成的食品在美國某些市場上已有銷售。此外,中間偃麥草對小麥銹病、白粉病等病害免疫,對黑穗病、葉枯病和根腐病等病害具有高抗性,能與小麥雜交,已經成為小麥遺傳改良中具有重要利用價值的野生親本之一。李振聲院士獲得了偃麥草與小麥的雜交種子,并將其廣泛用作與小麥遠緣雜交的野生親本材料,通過雜交將其優(yōu)良基因轉移進小麥,已經成為小麥遺傳改良的有效途徑。

    黑龍江省分子細胞遺傳與遺傳育種重點實驗室研究團隊對中間偃麥草種質資源在直接馴化改良、遠緣雜交改良等方面進行了研究,并對利用基因組學育種技術進行中間偃麥草遺傳育種研究的重點和選育策略進行了展望。全基因組選擇使用單核苷酸多態(tài)性(single nucleotide polymorphism,SNP)標記信息和表型特征值估計候選個體育種值,根據育種值的大小對候選個體進行篩選。作為育種領域的新興技術,與傳統的僅使用個體表型的估計育種值相比,全基因組選擇具有準確度更高、世代間隔短、育種效率高等優(yōu)點,目前已經被越來越多地應用在農作物的育種領域。

    在國外,Zhang等利用基因組選擇技術構建了基因標記數為3 883個、大小為1 126個的育種群體,并用不同的基因組選擇方法進行了研究,證明了使用基因組選擇技術進行中間偃麥草遺傳育種的可行性。Bajgain等使用全基因組關聯分析尋找中間偃麥草表型發(fā)育相關的基因座,利用顯著位點作為固定效應預測育種值,使用rrBLUP在SNP標記數為8 899個的數據集中進行育種值的估計,結果表明,該方法可使預測能力提高14%??紤]到環(huán)境因素對表型發(fā)育的影響,Bajgain等此后使用GBLUP、貝葉斯方法研究加性、顯性效應以及環(huán)境效應對中間偃麥草性狀發(fā)育的影響,在標記數為 8 899 個、群體大小為451個的基因數據集中進行了穗質量、穗長、每個花序小穗數、產量等9種性狀的全基因組選擇研究。Crain等研究了中間偃麥草 3 658 個個體的46個性狀,對18 357個SNP標記的基因組數據使用GBLUP方法進行了研究,發(fā)現預測能力隨著訓練群體規(guī)模、標記數量的增加而提高。

    本研究對標記數為23 495個、研究群體大小為5 521個的中間偃麥草數據采用隨機方式抽取個體構建訓練群體和校驗群體,使用交叉驗證的方式計算候選育種個體全基因組估計育種值和育種準確度,對比GBLUP、BayesA、BayesB、BayesCπ法的基因組育種值估計結果,以期找到適合中間偃麥草性狀育種分析的全基因組選擇方法來提高育種效率,為中間偃麥草全基因組育種提供參考。

    1 材料與方法

    本研究所用數據是從美國國立生物技術信息中心(NCBI)網站(https://www.ncbi.nlm.nih.gov/bioproject/)上下載的美國堪薩斯州薩利納市土地研究所的中間偃麥草的公共數據。在試驗中選擇中間偃麥草的自由脫粒率、穗產量、株高、種子質量、每個花序小穗數、穗長、落粒性等7種性狀進行研究,研究群體包含5 521個個體,基因組數據包括23 495個SNP標記位點。

    生物育種中的遺傳力也稱遺傳率,表示遺傳因素對表型發(fā)育所起作用的大小。由于表型數據中存在表型值缺失或無效的情況,因此需要對數據進行篩選,經過篩選后具有有效表型數據的7種性狀數據值的平均值、標準差、最大值和最小值的統計結果及對應的遺傳力如表1所示,其中自由脫粒、穗產量、株高、種子質量、每個花序小穗數、穗長、落粒性具有有效表型數據和SNP標記數據的群體中的個體數分別為5 019、5 230、5 488、4 703、5 013、4 512、4 694個。

    表1 表型數據及遺傳力

    1.1 統計模型

    本研究選用的統計模型如下所示,

    式中:為個體的性狀表型值;為標記的數量;為平均表型值;為第個分子標記的效應值;為個體第個標記的基因型編碼(編碼符號為0、1或2);為個體的隨機殘差效應值。

    1.2 BayesA

    1.3 BayesB

    BayesB假設只有部分位點效應能夠影響到遺傳效果,在所有位點的遺傳方差分布中,只有少數位點具有遺傳方差。因此,在BayesB中使用參數控制位點是否具有標記效應。

    1.4 BayesCπ

    BayesCπ假設所有位點具有1個公共方差,而不是每個位點都有自己的方差,并且假設控制標記效應為0的參數()未知。公共方差符合自由度=42、縮放因子為的尺度逆卡方分布,其中的推導過程與BayesB中縮放因子()的推導過程相同。概率()未知,并且假設符合(0,1)間的均勻分布,符合概率為(1-)的SNP效應由多元分布的混合模型計算得到,標記效應值也是使用正態(tài)分布取樣的方式進行估計的。

    春秋末期,宗法制逐步被破壞,導致天子與各諸侯國的關系也發(fā)生了重大變化。天子作為天下共主的地位和權威進一步喪失?!按呵飼r猶嚴祭祀,重聘享,而七國則無其事矣;春秋時,猶論宗姓氏族,而七國則無一言及之矣;春秋時,猶宴會賦詩,而七國則不聞矣;春秋時猶赴告策書,而七國則無有矣。邦無定交,土無定主?!盵23]P715宗法、會盟、祭祀、等諸多方面都發(fā)生了重大變化。

    1.5 GBLUP

    在GBLUP中,基于系譜的親緣關系矩陣由遺傳關系矩陣(G陣)代替,然后使用最佳線性無偏估計計算育種值。G矩陣的公式如下:

    式中:為標記的數量;為標記的等位基因頻率;為SNP標記效應的中心關聯矩陣。

    1.6 GEBV計算

    全基因組估計育種值(genomic estimated breeding value,Gebv)使用如下公式計算:

    式中:為個體的基因組估計育種值。

    1.7 交叉驗證

    本研究選用隨機掩蔽交叉驗證方法將中間偃麥草的表型數據、基因型數據分成校驗集合和訓練集合,使用5倍交叉驗證方法對群體進行分析。每次隨機抽取約1/5個體構成校驗集合,剩余的個體構成訓練集合,具體分組見表2。對于每種性狀,重復基因組育種值估計10次,以10次估計結果的平均準確度作為全基因組估計的準確度。

    表2 訓練群體和校驗群體分組情況

    數據分組后,使用訓練群體中的表型數據、基因型數據估計位點效應值,以BayesA為例,交叉驗證全基因組估計實現流程如圖1所示。驗證階段使用估計的效應值,校驗群體中的基因型編碼計算校驗群體中個體的育種值。

    1.8 估計準確性標準

    將全基因組估計所得育種值與表型數據求相關系數,作為評判標準,計算公式:

    2 結果與分析

    2.1 試驗環(huán)境與試驗相關參數

    本研究的環(huán)境:AMD RyZen 5 1600六核(12線程)中央處理器(CPU),32 G DDR4 2 666 MHz內存,希捷2 TB(1萬轉)硬盤;中文版64位Windows 7操作系統。

    在試驗中,BayesA、BayesB、BayesCπ中的burnin和總取樣運行次數分別設置成400 0、20 000。=4012、=0.002 0作為BayesA中的先驗分布參數;=4.2作為BayesCπ中的先驗分布參數,=4234、=0.042 9作為BayesB中的先驗分布參數,在BayesB、BayesCπ算法中,有效基因位占比()=0.01。

    2.2 估計準確度

    本試驗采用5倍交叉驗證,經過10次運行后,計算每次運行產生的估計準確度。由于訓練群體和校驗群體的構建采用隨機抽取的方式選擇個體,不同個體的SNP標記數據和表型數據存在差異,使得每次最終估計的育種值準確度存在差異,因此選用育種值估計準確度的“均值±標準偏差”作為最終估計結果。由7種中間偃麥草全基因組育種值估計的準確度(表3)可以看出,GBLUP預測的準確度最低,標準偏差也最低,說明GBLUP的性能最穩(wěn)定。在用BayesB方法得出的結果中,4個最優(yōu),1個次優(yōu);在用BayesCπ方法得出的結果中,2個最優(yōu),2個次優(yōu);在用BayesA方法得出的結果中,1個最優(yōu),3個次優(yōu);在用GBLUP方法得出的結果中,1個次優(yōu)??傮w而言,BayesB方法的估計準確度表現最優(yōu)異。

    表3 全基因組育種值估計的準確度

    2.3 運行時間

    在10次交叉驗證運行后,產生的10次運行時間如表4所示??梢钥闯?,同一性狀的3種貝葉斯方法的運行時間均高于GBLUP法的運行時間,而且3種貝葉斯方法的運行時間也有差異,排序為BayesB>BayesA>BayesCπ。

    表4 不同方法的運行時間

    在分析單次貝葉斯全基因組估計的運行時間時,選擇7種性狀的第1次運行時間進行比較。如圖2所示,BayesA法的運行時間大約是BayesB法的75%,BayesCπ法的運行時間大約是BayesB法的60%。

    2.4 中間偃麥草性狀全基因組估計育種值

    由于不同全基因組選擇方法估計的中間偃麥草育種值不同,因此選擇10個不同個體,使用對各性狀估計精度最高的方法計算全基因組估計育種值,結果見表5。在育種實踐中,不同性狀的表型值間存在較大差異,使得估計所得育種值在數量級上的差異較大。因此在育種過程中可以針對某一種性狀進行個體的選擇,也可以針對所有性狀對候選個體進行取舍。在對所有研究性狀的估計育種值進行選擇時,考慮到不同育種值間數量級的差異,通常情況下對不同性狀賦予不同的權重,采用加權求和的方式或者根據綜合公式計算綜合育種值,再根據綜合育種值篩選候選個體。

    表5 中間偃麥草不同性狀全基因組的估計育種值

    3 討論

    育種過程的關鍵環(huán)節(jié)是估計育種值,作為新興的育種技術,全基因組選擇在育種值估計準確度、縮短育種世代間隔、遺傳力低的性狀的選育中表現出了優(yōu)異性能。中間偃麥草作為一種優(yōu)良的牧草,對小麥病害具有高抗性,甚至具有免疫性,是雜交小麥重要的野生親本,進行中間偃麥草的全基因組育種研究既能促進優(yōu)良牧草的改良,又能為雜交小麥育種提供重要的遺傳親本,具有重要的實用價值。

    本試驗對中間偃麥草自由脫粒率、穗產量、株高、種子質量、每個花序小穗數、穗長、落粒性7種性狀的全基因組進行了選擇性分析。結果顯示,從育種值估計準確度上看,BayesB法對自由脫粒率、穗產量、株高、落粒性4種性狀的估計準確度最高;BayesCπ法對種子質量、每個花序小穗數的估計準確度最高;BayesA法對穗長性狀的估計準確度最高。究其原因,全基因組選擇育種值估計的準確度受到群體大小、全基因組選擇方法、遺傳力、種群數據等因素的影響,很少有某種方法能夠保證在所有性狀育種值的估計中都能獲得最優(yōu)的準確度。

    GBLUP使用遺傳關系矩陣代替親緣關系矩陣,利用最佳線性無偏估計方法估計育種值,相校于傳統的基于表型的BLUP方法,更能反映個體發(fā)育過程中性狀的遺傳性能,在保持BLUP較短運行時間的情況下,準確性更高,是當前動植物育種領域全基因組選擇的重要工具。GBLUP假設所有位點都具有遺傳效應,并未考慮單個SNP標記效應,在估計過程中不會出現由于隨機取樣產生的偶然性因素,因此性能穩(wěn)定,中間偃麥草7種性狀的全基因組估計結果的偏差最小。

    貝葉斯方法在基因組選擇過程中求解每個SNP位點的效應值,更能反映基因位點對遺傳影響的真實情況。BayesA假設所有位點對性狀發(fā)育產生遺傳效應,BayesB、BayesCπ假設只有少量基因位點可以影響性狀的發(fā)育,其余的位點效應可以忽略不計,從目前全基因組關聯分析的結果可以看出,少量顯著位點效應影響性狀的發(fā)育更能反映遺傳效果。BayesB假設每個位點都具有自己的方差,BayesCπ假設所有位點都有1個公共的方差,BayesB更能反映SNP效應的真實情況,因此在中間偃麥草7種性狀全基因組育種值的估計精度中,BayesB法得出的最優(yōu)結果最多。

    就運行時間而言,使用SNP標記數據計算遺傳關系矩陣占據了GBLUP絕大部分運行時間,在計算機運行環(huán)境相同的情況下,全基因組選擇群體大小、個體SNP標記數量決定了基因型數據,進而主導著GBLUP的運行時間。BayesA、BayesB、BayesCπ全基因組選擇方法假設SNP標記效應值符合正態(tài)分布,該正態(tài)分布的參數符合尺度逆卡方分布,通過隨機取樣的方式從尺度逆卡方分布中產生正態(tài)分布函數的參數,然后使用馬爾科夫鏈蒙特卡羅(Markov chain Monte Carlo,MCMC)取樣技術(BayesA使用Gibbs取樣,BayesB、BayesCπ使用Metropolis Hastings取樣),利用正態(tài)分布取樣函數產生標記效應值。MCMC的運行過程包括burn-in訓練階段和sampling取樣階段,為了消除取樣偶然性的影響,burn-in階段、sampling階段的迭代次數通常設成較大的值,貝葉斯全基因組選擇方法屬于高度密集型計算任務,往往需要較長的運行時間。BayesCπ假設所有位點具有1個公共的方差,BayesB、BayesA假設每個位點具有自己的方差,需要為每個位點估計正態(tài)分布方差,而且解決相同的問題,Metropolis Hastings取樣運行時間比Gibbs 取樣運行時間更長,因此在中間偃麥草群基因組選擇中,GBLUP的運行時間最短。在3種貝葉斯方法中,BayesB的運行時間最長,BayesCπ的運行時間最短。貝葉斯方法的全基因組選擇準確度高,但漫長的運行時間是育種實踐應用過程中的不利因素,使用更高主頻的高性能處理器結合并行計算技術縮短運行時間,將有助于提高貝葉斯方法全基因組選擇技術的實用性。

    4 結論

    本研究使用GBLUP、BayesA、BayesB和BayesCπ進行中間燕麥草的自由脫粒率、穗產量、株高、種子質量、每個花序小穗數、穗長、落粒性7種性狀的全基因組選擇,經過5倍的交叉驗證發(fā)現,在4種方法育種值的估計結果中,BayesB的準確度最優(yōu),運行時間最長。進行中間偃麥草全基因組的選擇需要針對不同性狀選擇最優(yōu)的全基因組選擇方法。在育種實踐中,綜合考慮不同方法育種性狀估計準確度和算法運行時間才能在中間偃麥草育種中高效地利用全基因組選擇技術,提高育種效率。

    猜你喜歡
    麥草準確度表型
    小麥-中間偃麥草2A/6St代換系014-459的分子細胞遺傳學鑒定
    作物學報(2022年2期)2022-11-06 12:08:56
    一根麥草
    快樂語文(2020年14期)2020-07-04 00:02:00
    幕墻用掛件安裝準確度控制技術
    建筑科技(2018年6期)2018-08-30 03:40:54
    建蘭、寒蘭花表型分析
    現代園藝(2017年21期)2018-01-03 06:41:32
    動態(tài)汽車衡準確度等級的現實意義
    閃爍
    西部(2016年6期)2016-05-14 19:09:59
    GABABR2基因遺傳變異與肥胖及代謝相關表型的關系
    慢性乙型肝炎患者HBV基因表型與血清學測定的臨床意義
    72例老年急性白血病免疫表型分析
    高爐重量布料準確度的提高
    天津冶金(2014年4期)2014-02-28 16:52:58
    精品一区二区三区av网在线观看| 黄色女人牲交| 国产精品精品国产色婷婷| a在线观看视频网站| 成年免费大片在线观看| 亚洲午夜理论影院| 久久久久亚洲av毛片大全| 国产精品,欧美在线| 性色av乱码一区二区三区2| 国产麻豆成人av免费视频| 伦理电影免费视频| 精品高清国产在线一区| 国产高清视频在线播放一区| 免费在线观看完整版高清| 精品免费久久久久久久清纯| 午夜激情福利司机影院| 99国产精品99久久久久| 久久欧美精品欧美久久欧美| 亚洲av电影在线进入| 久9热在线精品视频| 正在播放国产对白刺激| 日韩中文字幕欧美一区二区| 久久99热这里只有精品18| 久久人妻福利社区极品人妻图片| 国产成人欧美在线观看| 亚洲最大成人中文| 男女下面进入的视频免费午夜 | 久久热在线av| 国产精品1区2区在线观看.| 亚洲自拍偷在线| 露出奶头的视频| 免费在线观看黄色视频的| 精品一区二区三区av网在线观看| 999精品在线视频| 免费在线观看完整版高清| 国产成人影院久久av| 亚洲第一电影网av| 动漫黄色视频在线观看| 老汉色av国产亚洲站长工具| 精品熟女少妇八av免费久了| 热re99久久国产66热| 成人av一区二区三区在线看| 亚洲av成人av| 久久午夜亚洲精品久久| 成人欧美大片| 欧美又色又爽又黄视频| 99国产综合亚洲精品| 99热6这里只有精品| 啦啦啦韩国在线观看视频| 精品不卡国产一区二区三区| 免费电影在线观看免费观看| 欧美日韩一级在线毛片| 午夜日韩欧美国产| 精品久久久久久,| 精品国产乱码久久久久久男人| 在线看三级毛片| 看片在线看免费视频| 悠悠久久av| 久久久水蜜桃国产精品网| 午夜精品在线福利| 中文字幕高清在线视频| 国产亚洲精品第一综合不卡| 香蕉丝袜av| 国产区一区二久久| 激情在线观看视频在线高清| 男女做爰动态图高潮gif福利片| 亚洲黑人精品在线| 亚洲天堂国产精品一区在线| 999精品在线视频| www.999成人在线观看| 一本久久中文字幕| 亚洲欧美日韩高清在线视频| 久久精品人妻少妇| 少妇裸体淫交视频免费看高清 | 色尼玛亚洲综合影院| 国产精品国产高清国产av| 巨乳人妻的诱惑在线观看| 黄频高清免费视频| 国产精品 欧美亚洲| 亚洲一区中文字幕在线| 91麻豆精品激情在线观看国产| www.www免费av| 国产又黄又爽又无遮挡在线| 国产一级毛片七仙女欲春2 | 欧美 亚洲 国产 日韩一| 久久久国产欧美日韩av| 欧美日韩乱码在线| 精品欧美一区二区三区在线| 日本成人三级电影网站| 人妻丰满熟妇av一区二区三区| 美女国产高潮福利片在线看| 久久狼人影院| 精品免费久久久久久久清纯| 久久午夜综合久久蜜桃| 丝袜在线中文字幕| 嫩草影院精品99| 国产日本99.免费观看| 欧美日本视频| 男女下面进入的视频免费午夜 | 久久久久久久久中文| 少妇粗大呻吟视频| 亚洲狠狠婷婷综合久久图片| 成人国语在线视频| 免费无遮挡裸体视频| 国产精品久久久人人做人人爽| 久久久久国内视频| 国产精品一区二区免费欧美| 伊人久久大香线蕉亚洲五| 91成年电影在线观看| 制服诱惑二区| 久久精品国产综合久久久| 亚洲真实伦在线观看| 午夜亚洲福利在线播放| 女人被狂操c到高潮| 亚洲国产精品合色在线| 久热这里只有精品99| 色精品久久人妻99蜜桃| 亚洲人成网站高清观看| 成人国产一区最新在线观看| 精品免费久久久久久久清纯| 精品第一国产精品| 精品人妻1区二区| 免费女性裸体啪啪无遮挡网站| 成年女人毛片免费观看观看9| 亚洲va日本ⅴa欧美va伊人久久| 一级a爱片免费观看的视频| 欧美zozozo另类| 成人18禁在线播放| 悠悠久久av| 精品久久久久久久毛片微露脸| 成人国语在线视频| 黄网站色视频无遮挡免费观看| 波多野结衣高清作品| 久99久视频精品免费| 精品国产国语对白av| 日韩成人在线观看一区二区三区| 男女视频在线观看网站免费 | 国产精品1区2区在线观看.| 18禁国产床啪视频网站| x7x7x7水蜜桃| 亚洲欧美精品综合一区二区三区| 国产私拍福利视频在线观看| 两人在一起打扑克的视频| 久久久久久久久中文| 亚洲第一av免费看| 一本久久中文字幕| 久久 成人 亚洲| 欧美日本亚洲视频在线播放| 欧美成人一区二区免费高清观看 | 精品日产1卡2卡| 欧美日韩精品网址| 亚洲,欧美精品.| 国产真实乱freesex| 中文字幕久久专区| 欧美又色又爽又黄视频| 免费高清视频大片| 高清毛片免费观看视频网站| 亚洲欧洲精品一区二区精品久久久| 国产精品1区2区在线观看.| 亚洲国产精品999在线| 制服丝袜大香蕉在线| aaaaa片日本免费| 人人妻人人看人人澡| 男人舔女人的私密视频| 色综合亚洲欧美另类图片| 一夜夜www| 久9热在线精品视频| 欧美久久黑人一区二区| 岛国在线观看网站| 此物有八面人人有两片| 一进一出好大好爽视频| 免费女性裸体啪啪无遮挡网站| 麻豆国产av国片精品| 欧美另类亚洲清纯唯美| 精品久久久久久久人妻蜜臀av| 少妇 在线观看| 精华霜和精华液先用哪个| 亚洲狠狠婷婷综合久久图片| av福利片在线| 国产精品电影一区二区三区| 丝袜在线中文字幕| 国产一区二区激情短视频| 91国产中文字幕| 丁香六月欧美| 国产一区在线观看成人免费| 99精品久久久久人妻精品| 两性午夜刺激爽爽歪歪视频在线观看 | 无限看片的www在线观看| 国产蜜桃级精品一区二区三区| 久久久国产成人免费| 超碰成人久久| 天堂√8在线中文| 首页视频小说图片口味搜索| 久久婷婷成人综合色麻豆| 日本免费a在线| 国产又黄又爽又无遮挡在线| 巨乳人妻的诱惑在线观看| 两性夫妻黄色片| 国产伦在线观看视频一区| av欧美777| 日韩有码中文字幕| 日本成人三级电影网站| 51午夜福利影视在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲真实伦在线观看| 在线观看午夜福利视频| 老司机深夜福利视频在线观看| 亚洲精品国产区一区二| 日韩精品免费视频一区二区三区| 久久国产精品男人的天堂亚洲| 窝窝影院91人妻| 免费在线观看影片大全网站| 国产精品二区激情视频| 欧美日韩福利视频一区二区| 好看av亚洲va欧美ⅴa在| 久久精品夜夜夜夜夜久久蜜豆 | 成人三级黄色视频| 免费看美女性在线毛片视频| 国产成人欧美在线观看| 一二三四在线观看免费中文在| 国产91精品成人一区二区三区| 色精品久久人妻99蜜桃| 丁香欧美五月| 一个人免费在线观看的高清视频| 久久久国产成人免费| 一级毛片高清免费大全| 午夜福利18| 最近最新中文字幕大全免费视频| 法律面前人人平等表现在哪些方面| 久久国产精品影院| 长腿黑丝高跟| 亚洲人成伊人成综合网2020| 欧美日韩瑟瑟在线播放| 欧美激情 高清一区二区三区| 国产精品美女特级片免费视频播放器 | 国产野战对白在线观看| 午夜亚洲福利在线播放| 日本在线视频免费播放| 亚洲三区欧美一区| 亚洲激情在线av| 俺也久久电影网| 久久国产精品人妻蜜桃| 黄色 视频免费看| 国产午夜精品久久久久久| 国产97色在线日韩免费| 777久久人妻少妇嫩草av网站| 美女高潮喷水抽搐中文字幕| 一级作爱视频免费观看| 两个人看的免费小视频| 国产乱人伦免费视频| 97人妻精品一区二区三区麻豆 | 三级毛片av免费| 一级毛片高清免费大全| 伊人久久大香线蕉亚洲五| 女生性感内裤真人,穿戴方法视频| 禁无遮挡网站| 50天的宝宝边吃奶边哭怎么回事| 12—13女人毛片做爰片一| 亚洲熟妇熟女久久| 久久久精品国产亚洲av高清涩受| 日日摸夜夜添夜夜添小说| www.自偷自拍.com| 免费观看人在逋| 亚洲av电影不卡..在线观看| 国产男靠女视频免费网站| 久久九九热精品免费| 成人国产综合亚洲| 日本a在线网址| 免费一级毛片在线播放高清视频| 99精品欧美一区二区三区四区| 精品日产1卡2卡| 亚洲av成人不卡在线观看播放网| 99riav亚洲国产免费| 色婷婷久久久亚洲欧美| 成人亚洲精品一区在线观看| 午夜免费鲁丝| 亚洲av美国av| 欧美大码av| 日韩 欧美 亚洲 中文字幕| x7x7x7水蜜桃| 欧美+亚洲+日韩+国产| 大香蕉久久成人网| 欧美不卡视频在线免费观看 | 午夜老司机福利片| 丝袜在线中文字幕| 又大又爽又粗| 在线观看66精品国产| 日本免费a在线| 黄色毛片三级朝国网站| 天堂√8在线中文| 级片在线观看| 国产欧美日韩一区二区三| 亚洲第一欧美日韩一区二区三区| ponron亚洲| 麻豆av在线久日| 日本熟妇午夜| 精品免费久久久久久久清纯| 国产免费男女视频| or卡值多少钱| 国产成人精品久久二区二区免费| 国产三级黄色录像| 亚洲 欧美一区二区三区| 欧美激情 高清一区二区三区| 精品久久久久久成人av| 99re在线观看精品视频| 久久精品国产99精品国产亚洲性色| 日本免费一区二区三区高清不卡| 国产av又大| 亚洲欧美精品综合一区二区三区| 欧美黑人巨大hd| 亚洲专区字幕在线| 国内毛片毛片毛片毛片毛片| 黑人操中国人逼视频| 变态另类丝袜制服| 窝窝影院91人妻| 国产精品二区激情视频| 亚洲欧美日韩高清在线视频| 亚洲熟妇中文字幕五十中出| 精品国产乱子伦一区二区三区| 国产私拍福利视频在线观看| 久久午夜综合久久蜜桃| 国产91精品成人一区二区三区| 亚洲男人的天堂狠狠| 亚洲国产毛片av蜜桃av| 性欧美人与动物交配| 韩国精品一区二区三区| 一进一出抽搐动态| 亚洲 国产 在线| 99在线视频只有这里精品首页| 一区二区三区激情视频| 午夜福利视频1000在线观看| 久久精品91蜜桃| 日韩高清综合在线| 一边摸一边做爽爽视频免费| 久久久久国内视频| 一区二区三区激情视频| 两性夫妻黄色片| 国产精品一区二区精品视频观看| 亚洲黑人精品在线| 亚洲专区国产一区二区| 这个男人来自地球电影免费观看| 黄色成人免费大全| 国产成人影院久久av| 婷婷精品国产亚洲av| 大型av网站在线播放| 天天一区二区日本电影三级| 99国产精品一区二区蜜桃av| 国产精品一区二区精品视频观看| 免费看十八禁软件| 欧美三级亚洲精品| 欧美成狂野欧美在线观看| 中文亚洲av片在线观看爽| 视频在线观看一区二区三区| 久久国产精品影院| 大香蕉久久成人网| 麻豆久久精品国产亚洲av| 国产精华一区二区三区| 亚洲精品国产区一区二| 久99久视频精品免费| 成人av一区二区三区在线看| 一进一出抽搐动态| 亚洲精品一区av在线观看| aaaaa片日本免费| 99国产极品粉嫩在线观看| 午夜激情福利司机影院| 免费看a级黄色片| 很黄的视频免费| 人人妻,人人澡人人爽秒播| 精品一区二区三区四区五区乱码| 好男人电影高清在线观看| bbb黄色大片| 亚洲全国av大片| 国产极品粉嫩免费观看在线| 精品人妻1区二区| 国产视频内射| 亚洲第一av免费看| 久久伊人香网站| 午夜a级毛片| 午夜成年电影在线免费观看| 日韩欧美一区二区三区在线观看| 国产人伦9x9x在线观看| 宅男免费午夜| 国产99久久九九免费精品| 啦啦啦 在线观看视频| 丁香六月欧美| 极品教师在线免费播放| 色婷婷久久久亚洲欧美| 国产精品久久视频播放| 动漫黄色视频在线观看| 欧美日韩一级在线毛片| 女人高潮潮喷娇喘18禁视频| 久久草成人影院| 一级a爱片免费观看的视频| 在线天堂中文资源库| 亚洲人成网站高清观看| 美女高潮到喷水免费观看| 丁香六月欧美| 叶爱在线成人免费视频播放| 久久久久国产精品人妻aⅴ院| 亚洲一码二码三码区别大吗| 精品少妇一区二区三区视频日本电影| 99riav亚洲国产免费| 欧美三级亚洲精品| 18禁黄网站禁片午夜丰满| 国产精品影院久久| 老司机午夜福利在线观看视频| 国产又爽黄色视频| 国产久久久一区二区三区| 午夜福利一区二区在线看| 91大片在线观看| 亚洲成av片中文字幕在线观看| 啦啦啦观看免费观看视频高清| 亚洲精品美女久久av网站| 99精品在免费线老司机午夜| 天天添夜夜摸| 欧美国产精品va在线观看不卡| 国产精品久久久久久人妻精品电影| 亚洲av五月六月丁香网| 在线观看免费日韩欧美大片| 日韩免费av在线播放| 久久久精品国产亚洲av高清涩受| 97人妻精品一区二区三区麻豆 | 村上凉子中文字幕在线| 午夜久久久在线观看| 精品国产美女av久久久久小说| 久久国产精品人妻蜜桃| 亚洲中文av在线| 91九色精品人成在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国语自产精品视频在线第100页| 俺也久久电影网| 99久久久亚洲精品蜜臀av| 人人妻人人澡人人看| 国产三级在线视频| 日韩欧美在线二视频| 日韩欧美三级三区| 亚洲av中文字字幕乱码综合 | 国产精品1区2区在线观看.| 午夜免费激情av| 老汉色∧v一级毛片| 人人妻人人澡人人看| 在线视频色国产色| √禁漫天堂资源中文www| 在线免费观看的www视频| 免费观看精品视频网站| 黄色视频,在线免费观看| 日韩精品中文字幕看吧| 亚洲五月婷婷丁香| 亚洲av成人不卡在线观看播放网| 精品日产1卡2卡| 午夜福利18| 成年免费大片在线观看| 亚洲中文字幕日韩| 亚洲五月婷婷丁香| 国产成人一区二区三区免费视频网站| 亚洲一码二码三码区别大吗| 国产免费男女视频| 女生性感内裤真人,穿戴方法视频| 久久午夜综合久久蜜桃| 亚洲,欧美精品.| 在线播放国产精品三级| 999精品在线视频| АⅤ资源中文在线天堂| 18禁美女被吸乳视频| 热re99久久国产66热| 亚洲精品美女久久久久99蜜臀| 精品国产乱子伦一区二区三区| 久久天堂一区二区三区四区| 欧美日韩亚洲国产一区二区在线观看| 日韩中文字幕欧美一区二区| 男人的好看免费观看在线视频 | 午夜激情av网站| 亚洲一码二码三码区别大吗| 超碰成人久久| 免费观看人在逋| 国产熟女xx| 久久久久久久午夜电影| 窝窝影院91人妻| 变态另类丝袜制服| 精品国产乱子伦一区二区三区| 深夜精品福利| 黑人欧美特级aaaaaa片| 他把我摸到了高潮在线观看| √禁漫天堂资源中文www| avwww免费| 欧美日韩亚洲国产一区二区在线观看| 97人妻精品一区二区三区麻豆 | 国产精品99久久99久久久不卡| 99国产精品一区二区蜜桃av| 成人国产综合亚洲| 久久久水蜜桃国产精品网| 看黄色毛片网站| 在线天堂中文资源库| 精品国内亚洲2022精品成人| 国产欧美日韩精品亚洲av| 啦啦啦观看免费观看视频高清| 中文字幕人妻丝袜一区二区| 国产精品一区二区免费欧美| 香蕉av资源在线| 日韩欧美国产在线观看| 欧美黑人精品巨大| svipshipincom国产片| 一区福利在线观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲欧美日韩无卡精品| 日本撒尿小便嘘嘘汇集6| 天天躁夜夜躁狠狠躁躁| 99在线视频只有这里精品首页| 亚洲中文字幕日韩| 日韩高清综合在线| 国产精品免费视频内射| 日本 欧美在线| 国产精品免费视频内射| 国产精品久久视频播放| 少妇熟女aⅴ在线视频| 日本a在线网址| 在线十欧美十亚洲十日本专区| 日韩欧美免费精品| av电影中文网址| 国产亚洲精品久久久久5区| 国产伦在线观看视频一区| 免费看美女性在线毛片视频| 777久久人妻少妇嫩草av网站| 亚洲人成电影免费在线| 久久狼人影院| 免费一级毛片在线播放高清视频| 免费搜索国产男女视频| 大型黄色视频在线免费观看| 日本三级黄在线观看| 成人免费观看视频高清| 亚洲欧美精品综合一区二区三区| 三级毛片av免费| 婷婷六月久久综合丁香| 国产成人精品久久二区二区免费| 精品久久久久久久末码| 国产精品久久久久久亚洲av鲁大| 精品久久久久久久久久久久久 | 黄片播放在线免费| 少妇裸体淫交视频免费看高清 | 人人妻人人澡欧美一区二区| 日韩欧美国产在线观看| 欧美日韩亚洲综合一区二区三区_| 亚洲中文字幕日韩| www国产在线视频色| 美女高潮喷水抽搐中文字幕| 国产免费男女视频| 欧美色视频一区免费| 巨乳人妻的诱惑在线观看| 怎么达到女性高潮| а√天堂www在线а√下载| 亚洲精品国产区一区二| 色综合欧美亚洲国产小说| 日韩三级视频一区二区三区| 欧美日韩精品网址| x7x7x7水蜜桃| 免费电影在线观看免费观看| 国内精品久久久久精免费| 久久久水蜜桃国产精品网| 真人做人爱边吃奶动态| 色综合站精品国产| netflix在线观看网站| 国产黄色小视频在线观看| 午夜日韩欧美国产| 成人午夜高清在线视频 | www.www免费av| 老司机靠b影院| 免费在线观看成人毛片| 美女大奶头视频| 欧美国产日韩亚洲一区| 在线永久观看黄色视频| 亚洲天堂国产精品一区在线| 亚洲第一青青草原| 欧美一级a爱片免费观看看 | 香蕉丝袜av| 激情在线观看视频在线高清| 男女视频在线观看网站免费 | 亚洲成人久久性| 日韩一卡2卡3卡4卡2021年| 18禁黄网站禁片免费观看直播| 日本在线视频免费播放| 中文字幕人妻丝袜一区二区| 波多野结衣巨乳人妻| 国产一区二区三区视频了| 精品熟女少妇八av免费久了| 亚洲精品一区av在线观看| 中出人妻视频一区二区| 国产免费av片在线观看野外av| 国产在线精品亚洲第一网站| www.999成人在线观看| 成人欧美大片| 搡老妇女老女人老熟妇| 麻豆久久精品国产亚洲av| 欧洲精品卡2卡3卡4卡5卡区| 99国产综合亚洲精品| 欧美亚洲日本最大视频资源| 国产激情久久老熟女| 黄色视频,在线免费观看| 亚洲精品国产区一区二| 欧洲精品卡2卡3卡4卡5卡区| av在线天堂中文字幕| 国产成人欧美| 日本三级黄在线观看| 久久人妻福利社区极品人妻图片| 每晚都被弄得嗷嗷叫到高潮| 久久伊人香网站| 精品国产乱码久久久久久男人| 一区二区三区高清视频在线| 午夜成年电影在线免费观看| 国内毛片毛片毛片毛片毛片| 午夜视频精品福利| 制服人妻中文乱码| 男女之事视频高清在线观看| 好男人电影高清在线观看| 校园春色视频在线观看| 深夜精品福利|