• <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
    国产精品久久久久久精品电影 | 老鸭窝网址在线观看| 在线观看午夜福利视频| 日本撒尿小便嘘嘘汇集6| 九色国产91popny在线| 中文字幕精品免费在线观看视频| 一区二区三区高清视频在线| 精品人妻在线不人妻| 中文字幕另类日韩欧美亚洲嫩草| 久久久国产成人免费| 黄色丝袜av网址大全| 亚洲专区字幕在线| 大型黄色视频在线免费观看| 天堂动漫精品| 久久久久国内视频| 国产激情久久老熟女| 国产精品亚洲美女久久久| 纯流量卡能插随身wifi吗| 欧美激情 高清一区二区三区| 免费人成视频x8x8入口观看| www国产在线视频色| 丁香六月欧美| 国产一级毛片七仙女欲春2 | a级毛片在线看网站| 亚洲免费av在线视频| 脱女人内裤的视频| 国产亚洲精品久久久久久毛片| 国产不卡一卡二| 国产一区二区三区在线臀色熟女| 男人操女人黄网站| 国产aⅴ精品一区二区三区波| 少妇熟女aⅴ在线视频| 乱人伦中国视频| 日韩三级视频一区二区三区| 成年人黄色毛片网站| 色尼玛亚洲综合影院| 欧美 亚洲 国产 日韩一| 人人妻人人澡欧美一区二区 | 亚洲一区中文字幕在线| 日本欧美视频一区| 免费高清在线观看日韩| 热99re8久久精品国产| 一级毛片精品| 日本免费一区二区三区高清不卡 | 咕卡用的链子| 婷婷精品国产亚洲av在线| 老汉色av国产亚洲站长工具| 亚洲色图 男人天堂 中文字幕| 亚洲第一电影网av| 777久久人妻少妇嫩草av网站| 亚洲国产精品999在线| 90打野战视频偷拍视频| 在线观看免费午夜福利视频| 亚洲情色 制服丝袜| 又紧又爽又黄一区二区| 黄色女人牲交| 国产精品久久视频播放| 久久精品国产99精品国产亚洲性色 | 午夜福利影视在线免费观看| 国产精品影院久久| 90打野战视频偷拍视频| 久久这里只有精品19| 国产成人免费无遮挡视频| 悠悠久久av| 国产精品,欧美在线| 一级,二级,三级黄色视频| www国产在线视频色| 99久久99久久久精品蜜桃| 女性被躁到高潮视频| 亚洲情色 制服丝袜| 一级毛片高清免费大全| 好看av亚洲va欧美ⅴa在| 久久久久久免费高清国产稀缺| 亚洲欧美日韩无卡精品| 视频区欧美日本亚洲| 久久精品91无色码中文字幕| 老司机午夜福利在线观看视频| 人妻丰满熟妇av一区二区三区| 老熟妇乱子伦视频在线观看| 国产av在哪里看| 亚洲精品美女久久av网站| 69av精品久久久久久| 久久久精品国产亚洲av高清涩受| 免费在线观看日本一区| 亚洲精品一卡2卡三卡4卡5卡| 日韩欧美国产一区二区入口| 欧美黄色片欧美黄色片| 亚洲狠狠婷婷综合久久图片| 免费看十八禁软件| 精品无人区乱码1区二区| 91精品国产国语对白视频| 久久香蕉精品热| 男女下面进入的视频免费午夜 | 亚洲国产日韩欧美精品在线观看 | 亚洲最大成人中文| 国产精品影院久久| 精品乱码久久久久久99久播| 久久久久久久久久久久大奶| x7x7x7水蜜桃| 久久精品国产99精品国产亚洲性色 | 欧美激情久久久久久爽电影 | 国产欧美日韩综合在线一区二区| 窝窝影院91人妻| 午夜福利高清视频| 亚洲片人在线观看| 一区在线观看完整版| 老熟妇仑乱视频hdxx| 亚洲国产精品久久男人天堂| 丰满的人妻完整版| 欧美黑人精品巨大| 99re在线观看精品视频| 国产男靠女视频免费网站| 精品一品国产午夜福利视频| 免费在线观看完整版高清| 成人国产一区最新在线观看| 99国产极品粉嫩在线观看| 一级a爱片免费观看的视频| 国产精品一区二区三区四区久久 | 久久久国产成人精品二区| 中出人妻视频一区二区| 99香蕉大伊视频| www国产在线视频色| 欧美大码av| 国产1区2区3区精品| 中文字幕高清在线视频| 欧美日韩亚洲国产一区二区在线观看| 国产亚洲av嫩草精品影院| 国产黄a三级三级三级人| 成人特级黄色片久久久久久久| 级片在线观看| www国产在线视频色| 一个人观看的视频www高清免费观看 | 久久精品国产99精品国产亚洲性色 | 日韩欧美一区视频在线观看| 亚洲精品国产色婷婷电影| 99在线人妻在线中文字幕| 1024香蕉在线观看| 淫秽高清视频在线观看| 在线天堂中文资源库| 老熟妇仑乱视频hdxx| 成人国产综合亚洲| 最好的美女福利视频网| 成人永久免费在线观看视频| 淫妇啪啪啪对白视频| 我的亚洲天堂| 成人欧美大片| 精品久久久久久成人av| 99国产精品99久久久久| 12—13女人毛片做爰片一| 午夜福利一区二区在线看| 免费人成视频x8x8入口观看| 欧美一级毛片孕妇| 少妇熟女aⅴ在线视频| 久久人人97超碰香蕉20202| 亚洲精品美女久久久久99蜜臀| www.www免费av| 美女扒开内裤让男人捅视频| 色婷婷久久久亚洲欧美| 亚洲国产精品sss在线观看| or卡值多少钱| 精品国产一区二区久久| 色综合婷婷激情| 91字幕亚洲| 老汉色av国产亚洲站长工具| 久久久国产成人免费| 淫妇啪啪啪对白视频| 亚洲电影在线观看av| 欧美午夜高清在线| 久久久久久久精品吃奶| 久久久久国产一级毛片高清牌| 91在线观看av| 乱人伦中国视频| or卡值多少钱| 日本vs欧美在线观看视频| 久久国产亚洲av麻豆专区| 亚洲精品中文字幕在线视频| 久久精品国产清高在天天线| 久久久久国内视频| 黑人欧美特级aaaaaa片| 日本黄色视频三级网站网址| 88av欧美| 久久香蕉激情| 久久 成人 亚洲| 午夜两性在线视频| 精品国产美女av久久久久小说| av网站免费在线观看视频| 精品午夜福利视频在线观看一区| 亚洲专区字幕在线| 亚洲精品一区av在线观看| 国产精品爽爽va在线观看网站 | 欧美日韩乱码在线| 首页视频小说图片口味搜索| 亚洲午夜理论影院| 午夜福利高清视频| 露出奶头的视频| 久久人妻熟女aⅴ| 亚洲免费av在线视频| 久久精品国产综合久久久| 亚洲精品中文字幕在线视频| 99精品欧美一区二区三区四区| 18禁裸乳无遮挡免费网站照片 | 如日韩欧美国产精品一区二区三区| 高清在线国产一区| 亚洲一区二区三区色噜噜| 深夜精品福利| 成年人黄色毛片网站| 成人欧美大片| 久久人人精品亚洲av| 搡老妇女老女人老熟妇| 波多野结衣一区麻豆| 校园春色视频在线观看| 黑人巨大精品欧美一区二区mp4| 亚洲精品久久成人aⅴ小说| 99国产精品一区二区蜜桃av| 欧美大码av| 婷婷精品国产亚洲av在线| 免费女性裸体啪啪无遮挡网站| 欧美另类亚洲清纯唯美| 色哟哟哟哟哟哟| 99在线人妻在线中文字幕| 免费在线观看黄色视频的| 日韩大尺度精品在线看网址 | 久久午夜亚洲精品久久| 亚洲人成伊人成综合网2020| 国产免费男女视频| 国产又爽黄色视频| 日本 欧美在线| 色尼玛亚洲综合影院| 午夜两性在线视频| 啦啦啦 在线观看视频| 久久久国产成人免费| 精品国产乱子伦一区二区三区| 国产成人欧美| 大陆偷拍与自拍| 亚洲免费av在线视频| 亚洲自偷自拍图片 自拍| 国产av在哪里看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲avbb在线观看| 亚洲欧美日韩另类电影网站| 在线免费观看的www视频| tocl精华| 韩国av一区二区三区四区| 美女大奶头视频| 欧美成人一区二区免费高清观看 | 亚洲人成电影免费在线| 国产97色在线日韩免费| 成人国产综合亚洲| 亚洲男人天堂网一区| 99在线视频只有这里精品首页| 一级毛片精品| 91成人精品电影| 精品免费久久久久久久清纯| 国产人伦9x9x在线观看| 欧美亚洲日本最大视频资源| 淫秽高清视频在线观看| 狂野欧美激情性xxxx| 一本综合久久免费| 国产又爽黄色视频| 久久久久九九精品影院| 中文字幕最新亚洲高清| 国产亚洲精品综合一区在线观看 | 看免费av毛片| 久久久久久人人人人人| 亚洲人成77777在线视频| 国产野战对白在线观看| 男女之事视频高清在线观看| 国产精品一区二区免费欧美| 婷婷精品国产亚洲av在线| 免费人成视频x8x8入口观看| 国产在线精品亚洲第一网站| 日韩精品青青久久久久久| 中文字幕另类日韩欧美亚洲嫩草| 法律面前人人平等表现在哪些方面| 亚洲中文字幕日韩| 视频区欧美日本亚洲| 中文字幕色久视频| 国产精品98久久久久久宅男小说| 日韩三级视频一区二区三区| e午夜精品久久久久久久| 国产成人免费无遮挡视频| 国产又爽黄色视频| 国产熟女午夜一区二区三区| av电影中文网址| 日本 av在线| 真人一进一出gif抽搐免费| 黑人操中国人逼视频| 国产xxxxx性猛交| 亚洲五月天丁香| 精品一区二区三区av网在线观看| 少妇粗大呻吟视频| 搡老妇女老女人老熟妇| 久久人人爽av亚洲精品天堂| 精品人妻1区二区| 午夜亚洲福利在线播放| 国产精品永久免费网站| 999久久久精品免费观看国产| 日韩三级视频一区二区三区| 亚洲色图 男人天堂 中文字幕| 啦啦啦韩国在线观看视频| 中文亚洲av片在线观看爽| 日本 av在线| 别揉我奶头~嗯~啊~动态视频| 午夜福利影视在线免费观看| 18禁国产床啪视频网站| 亚洲成a人片在线一区二区| 国产精品一区二区精品视频观看| 最好的美女福利视频网| 91成年电影在线观看| √禁漫天堂资源中文www| 精品国产亚洲在线| 精品国产乱码久久久久久男人| 欧洲精品卡2卡3卡4卡5卡区| 国产高清视频在线播放一区| 日韩av在线大香蕉| 免费一级毛片在线播放高清视频 | 久久亚洲真实| 国产精品亚洲av一区麻豆| 色哟哟哟哟哟哟| 日韩成人在线观看一区二区三区| 91成年电影在线观看| 少妇的丰满在线观看| 免费无遮挡裸体视频| 欧美av亚洲av综合av国产av| 成年人黄色毛片网站| 一级毛片精品| 免费av毛片视频| 一区在线观看完整版| 亚洲自偷自拍图片 自拍| 岛国视频午夜一区免费看| 久9热在线精品视频| av电影中文网址| 男女下面进入的视频免费午夜 | 老汉色av国产亚洲站长工具| 欧美成人免费av一区二区三区| 久热爱精品视频在线9| 长腿黑丝高跟| 黄色片一级片一级黄色片| 国产在线精品亚洲第一网站| av中文乱码字幕在线| 在线观看日韩欧美| 黑人欧美特级aaaaaa片| 日韩成人在线观看一区二区三区| 亚洲人成电影观看| 韩国av一区二区三区四区| 宅男免费午夜| 女人被狂操c到高潮| 欧美亚洲日本最大视频资源| 日日干狠狠操夜夜爽| 久久青草综合色| 美女免费视频网站| 久久香蕉激情| 侵犯人妻中文字幕一二三四区| 99国产综合亚洲精品| 午夜成年电影在线免费观看| 亚洲少妇的诱惑av| 免费看十八禁软件| 国产aⅴ精品一区二区三区波| 一进一出抽搐gif免费好疼| av在线天堂中文字幕| 婷婷丁香在线五月| 国产欧美日韩一区二区三区在线| a级毛片在线看网站| 免费一级毛片在线播放高清视频 | 国产三级黄色录像| 日本三级黄在线观看| 99久久久亚洲精品蜜臀av| 亚洲人成77777在线视频| 国产高清videossex| 国产三级在线视频| 欧美成狂野欧美在线观看| 在线十欧美十亚洲十日本专区| 亚洲中文av在线| 精品国产一区二区久久| av免费在线观看网站| 99国产精品免费福利视频| 亚洲第一欧美日韩一区二区三区| 激情视频va一区二区三区| 香蕉丝袜av| 一级毛片女人18水好多| 日本 av在线| av福利片在线| 精品不卡国产一区二区三区| 精品国产乱码久久久久久男人| 久久人人爽av亚洲精品天堂| 色在线成人网| 1024视频免费在线观看| 久久人妻av系列| 一a级毛片在线观看| 18禁黄网站禁片午夜丰满| 色婷婷久久久亚洲欧美| 美女高潮到喷水免费观看| 欧美成人免费av一区二区三区| 精品电影一区二区在线| 色精品久久人妻99蜜桃| 久久久久久久久中文| 亚洲视频免费观看视频| 国产亚洲欧美98| 国产午夜福利久久久久久| 日本a在线网址| 窝窝影院91人妻| 欧美午夜高清在线| 久久伊人香网站| 老汉色∧v一级毛片| av在线天堂中文字幕| 两个人看的免费小视频| 一级a爱片免费观看的视频| 自拍欧美九色日韩亚洲蝌蚪91| 久久欧美精品欧美久久欧美| 久久久久九九精品影院| 激情在线观看视频在线高清| 欧美日韩亚洲综合一区二区三区_| 欧美最黄视频在线播放免费| 亚洲视频免费观看视频| 国产精品永久免费网站| av在线播放免费不卡| 亚洲第一电影网av| 国产亚洲精品第一综合不卡| 老司机午夜十八禁免费视频| 精品不卡国产一区二区三区| 免费观看人在逋| 精品卡一卡二卡四卡免费| 两个人免费观看高清视频| 久久久久精品国产欧美久久久| 国产精品久久久久久亚洲av鲁大| 69精品国产乱码久久久| 在线观看66精品国产| 中文字幕另类日韩欧美亚洲嫩草| 午夜日韩欧美国产| 天堂√8在线中文| 色尼玛亚洲综合影院| 久久久国产欧美日韩av| 18禁黄网站禁片午夜丰满| 久久久久久免费高清国产稀缺| 天堂动漫精品| 制服人妻中文乱码| 男女下面进入的视频免费午夜 | 国产精品 欧美亚洲| 久久人人精品亚洲av| x7x7x7水蜜桃| 国产午夜福利久久久久久| 18美女黄网站色大片免费观看| 色在线成人网| 欧美成狂野欧美在线观看| 一二三四在线观看免费中文在| 美女高潮到喷水免费观看| 很黄的视频免费| 一本综合久久免费| 国产又色又爽无遮挡免费看| 天堂影院成人在线观看| 国产精品久久久av美女十八| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲成人免费电影在线观看| 中亚洲国语对白在线视频| netflix在线观看网站| 亚洲熟妇熟女久久| 亚洲精品av麻豆狂野| 18禁国产床啪视频网站| 十分钟在线观看高清视频www| 久久人人精品亚洲av| 久久香蕉激情| 亚洲一区二区三区色噜噜| 一本久久中文字幕| 黑人巨大精品欧美一区二区mp4| 69av精品久久久久久| av电影中文网址| 精品久久久久久久毛片微露脸| 50天的宝宝边吃奶边哭怎么回事| 久久精品国产亚洲av香蕉五月| 成人国语在线视频| 亚洲精品在线观看二区| 精品一区二区三区四区五区乱码| 午夜视频精品福利| 国内毛片毛片毛片毛片毛片| 亚洲第一电影网av| 日本 av在线| 亚洲av电影不卡..在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 欧美色欧美亚洲另类二区 | 又大又爽又粗| 亚洲专区字幕在线| 中文字幕另类日韩欧美亚洲嫩草| 久久精品亚洲精品国产色婷小说| 在线天堂中文资源库| 欧美成人免费av一区二区三区| 最好的美女福利视频网| 免费不卡黄色视频| www.自偷自拍.com| 久久久久国内视频| 亚洲专区国产一区二区| 午夜免费激情av| 亚洲电影在线观看av| av超薄肉色丝袜交足视频| 国产熟女xx| cao死你这个sao货| 一二三四在线观看免费中文在| 中亚洲国语对白在线视频| 国产精品1区2区在线观看.| 欧美在线一区亚洲| 国产精品综合久久久久久久免费 | 一二三四在线观看免费中文在| 别揉我奶头~嗯~啊~动态视频| 成人三级黄色视频| 黄色a级毛片大全视频| 亚洲狠狠婷婷综合久久图片| 一级a爱视频在线免费观看| av福利片在线| 一a级毛片在线观看| 国产成人av激情在线播放| 国产又色又爽无遮挡免费看| 老司机福利观看| 久热这里只有精品99| 午夜福利成人在线免费观看| 亚洲精品美女久久久久99蜜臀| 99精品在免费线老司机午夜| 亚洲男人天堂网一区| 日本三级黄在线观看| 亚洲中文av在线| 精品高清国产在线一区| 久久久精品欧美日韩精品| 一本综合久久免费| 嫩草影院精品99| 国语自产精品视频在线第100页| 一级,二级,三级黄色视频| 色婷婷久久久亚洲欧美| 免费搜索国产男女视频| 波多野结衣巨乳人妻| 国产精品精品国产色婷婷| 国产伦一二天堂av在线观看| 亚洲av第一区精品v没综合| 亚洲一区二区三区不卡视频| 久久人妻福利社区极品人妻图片| 欧美性长视频在线观看| 亚洲欧美一区二区三区黑人| 亚洲成人久久性| 可以在线观看的亚洲视频| 神马国产精品三级电影在线观看 | 日韩精品免费视频一区二区三区| 九色国产91popny在线| 国产精品av久久久久免费| 欧美国产精品va在线观看不卡| 亚洲性夜色夜夜综合| 91麻豆av在线| 国产精品久久久久久精品电影 | 国产日韩一区二区三区精品不卡| 深夜精品福利| 日日夜夜操网爽| 黄色视频,在线免费观看| 一本大道久久a久久精品| 亚洲第一青青草原| 国产在线精品亚洲第一网站| а√天堂www在线а√下载| 国产精品乱码一区二三区的特点 | 日本在线视频免费播放| 老司机靠b影院| x7x7x7水蜜桃| 在线观看www视频免费| 麻豆av在线久日| 12—13女人毛片做爰片一| 99riav亚洲国产免费| 亚洲成人久久性| 欧美日韩福利视频一区二区| 久久精品91蜜桃| 18禁裸乳无遮挡免费网站照片 | 久久亚洲精品不卡| 美女扒开内裤让男人捅视频| 又紧又爽又黄一区二区| 免费在线观看黄色视频的| 欧美黑人精品巨大| 在线观看日韩欧美| av天堂在线播放| 亚洲av电影不卡..在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 欧洲精品卡2卡3卡4卡5卡区| 日韩视频一区二区在线观看| 精品久久久久久久久久免费视频| 欧美老熟妇乱子伦牲交| 久久精品国产99精品国产亚洲性色 | 黄色女人牲交| 午夜亚洲福利在线播放| 国产成人欧美| 色尼玛亚洲综合影院| 俄罗斯特黄特色一大片| 在线观看www视频免费| 老司机福利观看| 亚洲中文av在线| 久久精品亚洲精品国产色婷小说| 大香蕉久久成人网| 色在线成人网| www.熟女人妻精品国产| 亚洲熟妇中文字幕五十中出| 午夜激情av网站| 日韩大码丰满熟妇| 侵犯人妻中文字幕一二三四区| 亚洲五月天丁香| 天天一区二区日本电影三级 | 精品久久久久久成人av| 亚洲精品在线观看二区| 自拍欧美九色日韩亚洲蝌蚪91| 精品一区二区三区视频在线观看免费| 久久久精品国产亚洲av高清涩受| 美女午夜性视频免费| 法律面前人人平等表现在哪些方面| 99在线人妻在线中文字幕| 精品高清国产在线一区| 亚洲av熟女| 韩国精品一区二区三区| 法律面前人人平等表现在哪些方面| 999久久久精品免费观看国产| 亚洲人成伊人成综合网2020| 神马国产精品三级电影在线观看 | 99国产精品一区二区蜜桃av| 国产亚洲精品久久久久久毛片| 香蕉国产在线看| 久久婷婷成人综合色麻豆|