999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于SSR標記分析貴州50個玉米品種的雜合位點

2017-12-01 09:54:55楊文鵬王明春
西南農業學報 2017年11期

王 偉,關 琦,宋 欣,楊文鵬*,王明春

(1.貴州省農業科學院 旱糧研究所,貴州 貴陽 550006;2. 貴州大學 農學院,貴州 貴陽 550025; 3.沈陽市翔宇中學,遼寧 沈陽 110000;4.威海火炬高技術產業開發區農業經濟發展局,山東 威海 264209)

基于SSR標記分析貴州50個玉米品種的雜合位點

王 偉1,關 琦2,3,宋 欣2,4,楊文鵬1*,王明春1

(1.貴州省農業科學院 旱糧研究所,貴州 貴陽 550006;2. 貴州大學 農學院,貴州 貴陽 550025; 3.沈陽市翔宇中學,遼寧 沈陽 110000;4.威海火炬高技術產業開發區農業經濟發展局,山東 威海 264209)

【目的】為推進貴州玉米種質資源的有效利用、遺傳改良以及新品種的選育。【方法】用多態性豐富的85個SSR標記,選取50個來自貴州省內已審定的玉米品種進行分析。【結果】①共檢測出471個SSR等位基因,SSR位點的等位基因數平均為5.54個,PIC值平均為0.58,表明所用SSR標記具有豐富的遺傳多態性。②umc2048、bnlg1885、bnlg1083、bnlg2144、phi092、umc2094、umc1726、nc005、umc1491、umc1619和bnlg2228等11個SSR位點的He和Ho差值較高(gt;0.2),受人為選擇和近交等的影響較大。③雜合位點的品種數占總品種數的比值gt;0.5的SSR標記有47個,其中比值gt;0.8的SSR標記位點有8個。④第1染色體的20~55 cM區段,第4染色體的380~590 cM區段,第5染色體的430~495 cM區段,第6染色體的25~40和40~90 cM 2個區段,第7染色體的395~505 cM區段,第8染色體的75~130 cM區段和第9染色體的60~185 cM區段,是雜合位點的重要富集區域。【結論】這些雜合位點富集區域對貴州50個玉米雜交種的雜種優勢貢獻較大,富集區域含有的優良基因,有利于雜種優勢的基因組分子標記輔助選擇。

SSR;雜合位點;雜交種;玉米

【研究意義】玉米是貴州省重要的飼料和糧食作物,其常年播種面積約占全省糧食總面積的1/4,其產量約占全省糧食總產量的26 %。因此,玉米生產發展的快慢在貴州省糧食生產中具有舉足輕重的作用[1]。實踐證明,提高玉米產量最有效的途徑是新品種的選育[2-3]。遺傳多樣性分析是進行新品種選育的基礎,也是進行品種資源評價和合理利用的重要內容[4]。遺傳多樣性研究可以了解種質資源的遺傳背景、遺傳結構及親緣關系,為作物種質資源的利用提供參考和依據[5]。【前人研究進展】隨著現代分子生物學的發展,分子標記技術在玉米種質研究中越來越受到重視,主要有RFLP[6-7]、RAPD[8-9]、AFLP[10-11]、SSR[12-14]和SNP[15-16]等分子標記技術。其中,SSR標記由于其多態性高、穩定性好、共顯性、操作簡單和花費相對較低等優點,已廣泛應用于玉米自交系遺傳多樣性等研究。【本研究切入點】筆者等曾對貴州省2000年以來審定的部分玉米品種的親本自交系進行遺傳多樣性研究[17],【擬解決的關鍵問題】筆者擬進一步對審定的部分雜交種進行雜合位點分析,旨在探索貴州雜交種的基因組雜合位點及其富集區域,為推進貴州玉米種質資源的有效利用、遺傳改良以及新品種的選育提供理論依據。

1 材料與方法

1.1 材料

貴州省2000年以來審定的50個雜交玉米品種(表1)。其中,包括普通玉米品種41個(黃粒品種36個、白粒品種4個、黃白粒品種1個),糯玉米品種8個,甜玉米品種1個。

1.2 DNA提取與檢測

采用楊文鵬[18]改進的CTAB法提取幼苗葉片總DNA,用分光光度計檢測DNA的濃度與質量,將DNA樣品濃度稀釋至10 ng/μl,-20 ℃保存。

1.3 SSR標記分析

從www.maizegdb.org網站選取均勻分布在10條染色體上的251個SSR標記,并查找出其PCR引物序列,由上海捷瑞生物工程公司合成引物。采用王偉等[19]優化的方法進行PCR擴增和電泳檢測。

1.4 統計分析

將電泳圖譜上同一位置清晰出現的條帶記為1,無條帶記為0,缺失記為9。利用Systat10.2軟件以歐式方法聚類。利用CERVUS2.0軟件計算多態性信息量(Polymorphic Imformation Content,簡稱PIC)、每位點平均等位基因數(A)、預期雜合度(He)和實測雜合度(Ho)。將SSR標記帶型以A、H和U重新建立數據庫。與親本(父本或母本)一致的帶型記為A,雜合帶型記為H,缺失記為U,用GGT32軟件參照IBM2 2008 Neighbors圖譜繪制圖示基因型,統計雜合位點所占比例及其所在染色體區段,分析雜合位點富集區域。

表1 貴州50個玉米雜交種的來源

2 結果與分析

2.1 SSR位點的遺傳多態性

從251對SSR引物中篩選出PCR擴增效果好、多態性豐富的引物85對(表2)。85個SSR標記在供試的50個品種中共檢測出471個等位基因,每個位點的等位基因數平均為5.54個,變幅為2~17個。高于平均等位基因數的SSR位點有33個,主要集中在1.01~1.04,1.07~1.08,2.01,2.06,2.08,3.05,4.02,4.05,5.06~5.07,6.00~6.02,6.05,6.07,7.02~7.04,8.01,8.03,8.05,8.09,9.07和10.05~10.07區段,說明具有較高等位基因數SSR位點的分布具有趨同性,某些染色體區段具有較高的等位變異。SSR標記位點的平均多態信息量(PIC)為0.58,變幅為0.17~0.90。PICgt;0.5的SSR標記有43個,說明這些SSR標記具有豐富的遺傳多態性。

由圖1可見,在供試雜交種中,SSR標記位點的He和Ho差值大部分在0.2以下,說明這些SSR位點受人為選擇和近交等的影響較小。只有umc2048、bnlg1885、bnlg1083、bnlg2144、phi092、umc2094、umc1726、nc005、umc1491、umc1619、bnlg2228共11個SSR位點的He和Ho差值在0.2以上,說明這些SSR位點受人為選擇和近交等的影響較大。

表2 85個SSR標記在50個玉米雜交種中的遺傳多態性 實測雜合度

續表2 Continued table 2

編號No.標記Marker位置Bin等位基因數Alleles多態性信息量PIC實測雜合度Actualheteroz-ygosity預期雜合度Expectedheteroz-ygosity編號No.標記Marker位置Bin等位基因數Alleles多態性信息量PIC實測雜合度Actualheteroz-ygosity預期雜合度Expectedheteroz-ygosity24umc19733.0570.700.620.7467umc11308.0540.490.660.5925umc22663.0630.460.540.5368umc17288.0630.580.700.6626umc20813.0830.320.360.3869phi2333768.0960.580.620.6227phi0473.0930.400.490.4870umc12799.0040.520.330.5728umc20483.1050.740.210.7871umc23369.0240.670.840.7329umc10174.0140.570.640.6472umc12719.0350.320.320.3530umc12944.02110.800.780.8273umc15709.0440.320.360.3631nc0054.0580.750.510.7974umc10789.0540.580.180.6532phi0924.0830.440.940.5575umc21349.0540.480.460.5733umc15324.1030.370.270.4376umc23469.0630.530.060.6234umc14915.0030.560.370.6477bnlg1289.0790.840.810.8735nc0075.0140.650.530.7178phi04110.0040.440.410.5536phi1135.0330.450.450.5679umc138010.0030.500.440.5737umc16245.0420.370.550.4980umc203410.0230.380.460.4438umc20265.0540.380.280.4281phi05010.0340.380.380.4239phi0855.0630.390.480.5082umc193010.0440.490.340.5440umc10195.06120.840.880.8683umc147710.0580.420.500.4941bnlg18855.0760.700.260.7484bnlg219010.0670.760.430.8042bnlg11185.07100.730.670.7785bnlg145010.0770.800.360.8343bnlg1616.00130.880.880.9平均5.540.580.530.63

2.2 雜交種基因組的雜合性

由圖2可見,50個雜交種中10條染色體上的85個SSR標記位點,每個標記位點檢測到具有雜合位點的品種數占總品種數的比值不相同,以第1染色體1.01節段的umc1619最高,為0.98;以第2染色體2.10節段的umc1696和第6染色體6.05節段的umc1474最低,為0.04。

由表3可知,具有雜合位點品種數占總品種數比值gt;0.5的SSR標記位點,第1、2、3、4、5、6、7、8、9和10染色體上分別有8個、4個、3個、4個、3個、8個、6個、8個、2個和1個。經統計,其中比值gt;0.8的SSR標記位點為第1染色體的umc1619、第4染色體的phi092,第5染色體的umc1019,第6染色體的bnlg161和umc1018,第7染色體的umc1782,第8染色體bnlg1194,第9染色體的umc2336。

圖1 SSR標記在供試材料中檢測的Ho及He差值Fig.1 Difference between Ho and He of SSR markers in 50 local hybrid maize varieties

圖2 85個SSR標記的具有雜合位點品種數占總品種數的比值Fig.2 Ratios of hybrid maize varieties with heterozygous locus of 85 SSR markers to total hybrid maize varieties

從圖3可見,第1、2、3、5、8和10染色體的雜合位點比較分散,而第4、6、7和9染色體的雜合位點則相對集中。第4染色體的雜合位點主要集中在70~180和385~590 cM 2個區域。第6染色體的雜合位點主要集中在25~155和255~315 cM 2個區域。第7染色體的雜合位點主要集中在70~135,195~355和395~505 cM 3個區域。第9染色體的雜合位點主要集中在60~185和520~585 cM 2個區域。

由表4可見,第1染色體的20~55 cM區段,第4染色體的380~590 cM區段,第5染色體的430~495 cM區段,第6染色體的25~40和40~90 cM 2個區段,第7染色體的395~505 cM區段,第8染色體的75~130 cM區段和第9染色體的60~185 cM區段,雜合片段比例gt;0.8是非常重要的雜合位點富集區。

3 討 論

D.Botstein等[20]認為,PIC值gt;0.5的共顯性遺傳標記具有高度多態性,0.25lt;PIC值lt;0.5則具有中度多態性,只有lt;0.25為低度多態性。本研究中,85個SSR標記的平均PIC值為0.58,PIC值gt;0.5的SSR標記有53個,0.25lt;PIC值lt;0.5的有31個,lt;0.25的只有1個,說明本試驗選擇的85個SSR標記具有豐富的多態性。

表3 具有雜合位點品種數占總品種數的比值gt;0.5的SSR標記

注:表中括號內數字分別表示染色體上位置和遺傳距離。

Note: The figures in parentheses indicate position on the chromosome and genetic distance respectively.

表4 玉米50個雜交種中重要的雜合位點富集區域

注:與親本(父本或母本)一致帶型的區域為淺色,雜合帶型區域為深色。

Note: The tint and dark area presents the identical banding pattern and heterozygous banding pattern compared with parents (male or female parent) respectively.

雜合度可反映群體在多個位點上的遺傳變異[21],分為預期雜合度和實測雜合度,兩者越接近,表明該位點受外來選擇及近交等因素的影響較小。本研究中,85個SSR標記在50個玉米雜交種中檢測到的預期雜合度和實測雜合度差值不盡相同,多數差值在0.2以下,說明多數SSR位點受人為選擇和近交等因素的影響較小,可能這些SSR位點在基因組上比較保守,在玉米遺傳改良過程中被保留。而SSR標記umc2048、bnlg1885、bnlg1083、bnlg2144、phi092、umc2094、umc1726、nc005、umc1491、umc1619、bnlg2228共11個位點的差值在0.2以上,說明這些SSR位點受人為選擇和近交等因素的影響較大,可能是玉米遺傳改良過程中的熱點區域。

本研究50個玉米雜交種中,85個SSR標記位點檢測到雜合位點的個體數占總樣品數的比值不相同,比值gt;0.8的主要集中在少數SSR標記位點,如umc1619(1.01)、phi092(4.08)、umc1019(5.06)、bnlg161(6.00)、umc1018(6.00)、umc1782(7.04)、bnlg1194(8.01)和umc2336(9.02)。依據雜種優勢理論,雜合體優于純合體,那么這些雜合位點可能對供試材料的雜種優勢貢獻較大,或者含有更多的有利基因。另外,雜合片段gt;0.8的圖示基因型區段,主要富集在第1染色體的20~55 cM區段、第4染色體的380~590 cM區段、第5染色體的430~495 cM區段、第6染色體的25~40和40~90 cM 2個區段、第7染色體的395~505 cM區段、第8染色體的75~130 cM區段和第9染色體的60~185 cM區段。上述雜合位點在染色體上的富集區域,可能對貴州50個玉米雜交種的雜種優勢貢獻較大,因此在進行玉米自交系選育和雜交組配時,上述位點的遺傳特性值得考慮。由于相鄰位點間的連鎖關系,這些富集區域可能有利于雜種優勢的基因組分子標記輔助選擇。

郭戰勇等[22]利用單片段代換系的測交群體,在2種環境中,定位了11個與玉米籽粒性狀相關的雜種優勢位點,其中粒長1個,粒厚4個,百粒重6個。在這11個雜種優勢位點所在區段中,有9個區段全部或部分落在本研究雜合位點比值gt;0.5的圖示基因型區段之內,占81.82 %,其中有2個區段全部或部分落在本研究雜合位點比值gt;0.8的圖示基因型區段之內,占22.22 %。彭倩等[23]同樣利用單片段代換系的測交群體,在兩種環境中,定位了23個與玉米產量和果穗性狀相關的雜種優勢位點,其中穗長4個,穗粗4個,穗行數4個,行粒數7個,產量4個。在這23個雜種優勢位點所在區段中,有17個區段全部或部分落在本研究雜合位點比值gt;0.5的圖示基因型區段之內,占73.91 %,其中有1個區段全部落在本研究雜合位點比值大于0.8的圖示基因型區段之內,占5.88 %。另外,郭戰勇等[22-23]在單個環境和兩個環境定位的玉米籽粒、果穗性狀和產量性狀等的雜種優勢位點連鎖染色體片段有139個。本研究,雜合比值位點gt;0.5的圖示基因型區段有47個,其中有36個區段與他們定位的片段完全或部分重疊,占76.60 %。其中,雜合位點比值gt;0.8的圖示基因型區段有8個,有4個與其定位的片段完全或部分重疊,占50 %。因此,本項目研究結果對貴州玉米的遺傳改良和雜種優勢利用的分子標記輔助育種具有重要意義。

與親本(父本或母本)一致帶型的區域為淺色,雜合帶型區域為深色圖3 50個玉米雜交種10條染色體的圖示基因型Fig.3 Graphical genotype of 10 chromosomes in 50 hybrid maize varieties

4 結 論

本研究用85個多態性豐富的SSR標記,分析了50個貴州省審定的玉米品種,共檢測出471個SSR等位基因,每位點等位基因數平均為5.54個,PIC值平均為0.58,具有豐富的遺傳多態性。其中,umc2048、bnlg1885、bnlg1083、bnlg2144、phi092、umc2094、umc1726、nc005、umc1491、umc1619和bnlg2228共11個SSR位點的He和Ho差值gt;0.2,受人為選擇和近交的影響較大。雜合位點的品種數占總品種數的比值gt;0.5的SSR標記有47個,其中比值大于0.8的SSR標記位點有8個;雜合位點比值gt;0.8的圖示基因型區段有第1染色體的20~55 cM區段,第4染色體的380~590 cM區段,第5染色體的430~495 cM區段,第6染色體的25~40 cM和40~90 cM 2個區段,第7染色體的395~505 cM區段,第8染色體的75~130 cM區段和第9染色體的60~185 cM區段,是重要的雜合位點富集區域。這些雜合位點富集區域對貴州50個玉米雜交種的雜種優勢貢獻較大,富集區域含有的優良基因,有利于雜種優勢的基因組分子標記輔助選擇。

[1]陳澤輝. 貴州玉米種質改良及雜交種的選育研究進展[J].貴州農業科學,2005,33(B09):16-19.

[2]張世煌. 玉米種質創新和商業育種策略[J].玉米科學,2006,14(4):1-3,6.

[3]張仁和,薛吉全,趙 靜,等.陜西省玉米品種雜種優勢群和雜優模式分析[J].玉米科學,2007,15(S1):173-176.

[4]徐愛遐,馬朝芝,肖恩時,等.中國西部芥菜型油菜遺傳多樣性研究[J].作物學報,2008,34(5):754-763.

[5]孫 健,王敬國,劉化龍,等.黑龍江省主栽水稻品種的遺傳多樣性分析[J].作物雜志,2011(1):63-67.

[6]Livini C, Ajmone-Marsan P, Melchinger A E, et al. Genetic diversity of maize inbred lines within and among heterotic groups revealed by RFLPs[J]. Theoretical and Applied Genetics,1992,84(1/2):17-25.

[7]Messmer M M, Melchinger A E, Lee M, et al. Genetic diversity among progenitors and elite lines from the Iowa Stiff Stalk Synthetic (BSSS) maize population: comparison of allozyme and RFLP data[J]. Theoretical and Applied Genetics,1991,83(1):97-107.

[8]Carvalho V P, Ruas C F, Ferreira J M, et al. Genetic diversity among maize (ZeamaysL.) landraces assessed by RAPD markers[J]. Genetics and Molecular Biology,2004,27(2):228-236.

[9]Shieh G J, Thseng F S. Genetic diversity of Tainan-white maize inbred lines and prediction of single cross hybrid performance using RAPD markers[J]. Euphytica, 2002, 124(3): 307-313.

[10]Li Y, Shi Y S, Song Y C, et al. Analysis of genetic diversity in maize inbred lines based on AFLP markers[J].Maydica,2004,49:89-96.

[11]Vuylsteke M, Mank R, Brugmans B, et al. Further characterization of AFLP?data as a tool in genetic diversity assessments among maize (ZeamaysL.) inbred lines[J].Molecular Breeding,2000,6(3):265-276.

[12]Dev J, Rathour R. Genetic diversity among maize inbreds adapted to hill regions of India revealed by functional SSR markers [J].Maize Genetics Conference Abstracts,2008,50:P186.

[13]Warburton M L, Reif J C, Frisch M, et al. Genetic diversity in CIMMYT nontemperate maize germplasm: landraces, open pollinated varieties, and inbred lines[J]. Crop Science, 2008, 48(2):617-624.

[14]Enoki H, Sato H, Koinuma K. SSR analysis of genetic diversity among maize inbred lines adapted to cold regions of Japan[J]. Theoretical and Applied Genetics, 2002, 104(8):1270-1277.

[15]Yan J, Shah T, Warburton M L, et al. Genetic characterization and linkage disequilibrium estimation of a global maize collection using SNP markers[J]. PloS one, 2009, 4(12):e8451.

[16]Van Inghelandt D, Melchinger A E, Lebreton C, et al. Population structure and genetic diversity in a commercial maize breeding program assessed with SSR and SNP markers[J]. Theoretical and Applied Genetics, 2010, 120(7):1289-1299.

[17]楊文鵬, 關 琦, 楊留啟, 等. 貴州70份玉米自交系的SSR標記遺傳多樣性及其雜種優勢群分析[J]. 植物遺傳資源學報, 2011, 12(2): 241-248.

[18]楊文鵬.玉米一個高賴氨酸突變基因的定位及o2胚乳異雄核受精和o2位點等位變異的分析[J].武漢:華中農業大學,2005.

[19]王 偉,楊文鵬,關 琦,等.玉米SSR分子標記技術操作規程的優化[J].安徽農業科學,2008,36(11):4459-4465.

[20]D Botstein, R L White, M Skolnick,R W Davis. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. American Journal of Human Genetics, 1980, 32(3):314.

[21]孔繁玲.植物數量遺傳學[M].北京:中國農業大學出版社,2006:407.

[22]郭戰勇,呂盼晴,張向歌, 等.利用單片段代換系的測交群體定位玉米籽粒性狀雜種優勢位點[J].中國農業科學, 2016, 49(4): 621-631.

[23]彭 倩,薛亞東,張向歌,等.利用單片段代換系測交群體定位玉米產量相關性狀的雜種優勢位點[J].作物學報, 2016, 42(4): 482-491.

(責任編輯 劉忠麗)

HeterozygousLocusof50LocalHybridMaizeVarietiesBasedonSSRMarkerAnalysisinGuizhou

WANG Wei1, GUAN Qi2,3, SONG Xin2,4, YANG Wen-peng1*, WANG Ming-chun1

(1.Guizhou Institute of Upland Crops, Guizhou Academy of Agricultural Sciences, Guiyang Guizhou 550006, China; 2.College of Agriculture, Guizhou University, Guiyang Guizhou 550025, China; 3.Xiangyu Middle School, Liaoning Shenyang 110000, China; 4. Bureau of Agricultural Economic Development, Weihai Torch High-tech Industrial Development Zone, Weihai Shandong 264209,China)

【Objective】This study aimed to improve effective utilization, genetic improvement and variety breeding of maize germplasm resources in Guizhou.【Method】50 approved local maize varieties in Guizhou are analyzed by 85 SSR markers with abundance polymorphism. 【Result】(i) 471 SSR alleles are detected and the average value and PIC of alleles on SSR loci 5.54/locus and 0.58/locus, which indicates that the used SSR markers have the abundant genetic polymorphism. (ii) The difference in He and Ho among 11 SSR loci of umc2048, bnlg1885, bnlg1083, bnlg2144, phi092, umc2094, umc1726, nc005, umc1491, umc1619 and bnlg2228 is high (gt;0.2), which indicates the He and Ho are influenced by artificial selection and inbreeding greatly. (iii) There are 47 SSR marker loci with gt; 0.5 ratio of varieties with heterozygous locus in total varieties. Of which, the ratio of 8 SSR marker loci is gt; 0.8. (iv) The segments of 20-55 cM on Chromosome 1,380-590 cM on Chromosome 4,430-495 cM on Chromosome 5, 25-40 cM and 40-90 cM on Chromosome 6, 395-505 cM on Chromosome 7, 75-130 cM on Chromosome 8 and 60-185 cM on Chromosome 9 are the important enrichment areas of heterozygous loci. 【Conclusion】The enrichment areas of heterozygous loci have a great contribution to heterosis of 50 hybrid maize varieties in Guizhou. The enrichment enrichment areas with good genes are beneficial to genome marker-assisted selection (MAS) in heterosis of maize.

SSR; Heterozygous; Locus; Hybrid; Maize

S513

A

1001-4829(2017)11-2417-08

10.16213/j.cnki.scjas.2017.11.005

2017-07-11

貴州省科技計劃項目[黔科合重大專項字(2013)6022號,黔科合NY(2014)3020];貴州省農業科學院院專項[黔農科院院專項(2014)008];貴州省科研機構服務企業行動計劃項目[黔科合服企(2014)4002];貴州省產業技術體系項目(GZCYTX2013-0701);貴州省農業科學院自主創新專項[黔農科院自主創新科研專項字(2014)006]

王 偉(1982-),男,安徽臨泉人,副研究員,碩士研究生,從事玉米分子遺傳學與育種研究,E-mail:wwmaize@126.com,*為通訊作者:楊文鵬,E-mail:ywpmaize@126.com。

主站蜘蛛池模板: 亚洲国内精品自在自线官| 欧洲高清无码在线| 丝袜美女被出水视频一区| 久久成人18免费| 日韩小视频网站hq| 黄色网页在线播放| 国产aaaaa一级毛片| 黑色丝袜高跟国产在线91| 强奷白丝美女在线观看| 久久久久久久蜜桃| 人妻21p大胆| 国产成人凹凸视频在线| 国产成人久久综合777777麻豆| 97超爽成人免费视频在线播放 | 久青草网站| 在线观看国产小视频| AV不卡在线永久免费观看| 日韩成人午夜| 午夜爽爽视频| 亚洲免费福利视频| 日韩视频精品在线| 伊人久久大线影院首页| 老司机午夜精品视频你懂的| 中文字幕欧美日韩高清| 欧美影院久久| 人人妻人人澡人人爽欧美一区| 日韩免费成人| 免费av一区二区三区在线| 91久久青青草原精品国产| 日本a级免费| 91福利在线观看视频| 91精品情国产情侣高潮对白蜜| 免费无遮挡AV| 99热这里只有精品5| 国产无码网站在线观看| 久久semm亚洲国产| 久久婷婷色综合老司机| 亚洲二三区| 中国成人在线视频| 精品少妇人妻无码久久| 啦啦啦网站在线观看a毛片| 亚洲精品成人7777在线观看| 久久动漫精品| 精品久久高清| 一级成人a毛片免费播放| 亚洲va在线观看| 精品国产乱码久久久久久一区二区| 超碰aⅴ人人做人人爽欧美 | 波多野结衣一区二区三区AV| 国产又黄又硬又粗| 久久国产高清视频| 久青草国产高清在线视频| 国产精女同一区二区三区久| 国产成人精品免费视频大全五级| av在线无码浏览| 婷婷亚洲视频| 国产在线视频自拍| aaa国产一级毛片| 国产免费怡红院视频| 黄色网站在线观看无码| 久久精品午夜视频| 国产一级特黄aa级特黄裸毛片| 99成人在线观看| 色精品视频| 国产极品美女在线播放| 久久综合色播五月男人的天堂| 欧美一区二区精品久久久| 国产97公开成人免费视频| 亚洲无限乱码一二三四区| 天堂岛国av无码免费无禁网站| 国产精品美女自慰喷水| 91视频首页| 香蕉伊思人视频| 尤物特级无码毛片免费| 99久久免费精品特色大片| 天天综合网色中文字幕| 亚洲一区二区三区国产精品| 乱人伦中文视频在线观看免费| 国产凹凸视频在线观看| 中文天堂在线视频| 亚洲人成在线精品| 欧日韩在线不卡视频|