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

皂莢種質資源SCoT遺傳多樣性分析及指紋圖譜的構建

2017-05-30 10:48:04張安世駱揚范定臣張中海
廣西植物 2017年11期

張安世 駱揚 范定臣 張中海

摘要: 采用SCoT標記分析了18個皂莢種質的遺傳多樣性,并采用UPGMA法對18個皂莢種質進行了聚類分析。在此基礎上,通過篩選出的多態性條帶構建了18個皂莢種質的SCoT指紋圖譜。擴增結果表明:從51個SCoT引物中篩選了15個引物進行PCR擴增,共擴增出226條帶,其中多態性條帶216條,多態性比率為96.61%。各引物多態性信息含量(PIC)、觀測等位基因數(Na)、有效等位基因數(Ne)、Neis基因多樣性指數(H)和Shannons信息指數(I)的平均值分別為0.875 9、1.964 9、1.440 1、0.272 6、0.426 1。18個皂莢種質的遺傳相似系數在0.491 4~0.938 1之間,表明供試材料之間具有較豐富的遺傳多樣性。聚類分析結果表明: 在遺傳相似系數為0.60處可將18個皂莢種質分為3組,其中野皂莢單獨為一組,山皂莢和皂莢T聚為一組,其它皂莢材料聚為一組。利用3個引物擴增的8個多態性位點構建了18份皂莢種質資源的DNA指紋圖譜,可以將其區分并精準鑒定。該研究結果為皂莢種質的鑒定和新品種選育提供了一定的理論依據。

關鍵詞: 皂莢, SCoT, DNA指紋圖譜, 遺傳多樣性

中圖分類號: Q949.9, S718.46

文獻標識碼: A

文章編號: 10003142(2017)11137808

Abstract: Genetic diversity of eighteen Gleditsia sinensis germplasms was analyzed based on SCoT markers, the cluster analysis on eighteen G. sinensis germplasms was carried out by UPGMA method, and the SCoT fingerprints of eighteen G. sinensis germplasms were constructed by the polymorphic bands. The results of amplification showed that, 226 SCoT bands were obtained from fifteen primers selected from 51 SCoT primers, including 216 polymorphic bands, with a polymorphism rate of 96.61%. The average of the polymorphism information content(PIC), observed number of alleles(Na), effective number of alleles(Ne), Neis gene diversity (H) and Shannons information index (I) were 0.875 9, 1.964 9, 1.440 1, 0.272 6 and 0.426 1, respectively. The genetic similarity coefficients of these samples were between 0.491 4 and 0.938 1. Therefore, it indicates that there are rich genetic diversities among these materials. The results of cluster analysis showed that, eighteen G. sinensis germplasms could be divided into three groups when the genetic similarity coefficient was 0.60. G. heterophylla formed the first group, G. melanacantha and ZaojiaT formed the second group, and the others formed the third group. DNA fingerprints of eighteen samples were constructed from eight polymorphic loci amplified by three primers, and these materials could be distinguished and identified accurately. All these results provide the important information for the identification and breeding for new cultivars of G. sinensis.

Key words: Gleditsia sinensis, SCoT, DNA fingerprint, genetic diversity

皂莢屬(Gleditsia Linn.) 系豆科蘇木亞科, 全世界約12種, 主要分布于亞洲、美洲和熱帶非洲(顧萬春等,2003)。中國產8種, 引進1種, 其中中國皂莢是分布在我國的特有種, 廣泛分布于我國東北、華北、華東、華南等地, 多生長于平原、山谷及丘陵地區。皂莢能抗旱、耐鹽、耐高溫,具有多種多樣的生態屬性,遺傳多樣性程度高,是經濟林、用材林、防護林及園林綠化的理想樹種(顧萬春等,2003;楊海東,2003)。皂莢具有很高的藥用價值,皂莢、皂刺、皂根、皂葉均可入藥,特別是棘刺是中醫治療乳腺癌、肺癌等多種癌癥的常用配伍藥之一, 被列為抗癌中草藥(蔣建新等,2003)。同時,皂莢還是優良的城市抗污染劑和日用化工原料(王薊花等,2008), 殺蟲、殺鼠效果也非常好(張宏利等,2013)。目前國內外對皂莢的研究主要集中在皂莢果實與皂刺中活性成分(張宏利等,2013; 余俐佳等, 2016)、引種栽培(范定臣等,2015)及生態屬性(楊海東,2003)等方面,并開展了相關的遺傳學研究(蘭彥平和顧萬春,2006; 李偉等,2013)。最近,分子標記技術也已開始在皂莢研究中得到應用。邢俊連等(2017)成功地進行了皂莢ESTSSR引物的開發與篩選,且證實這些引物在其近緣種間具有很高的通用性,為以后開展皂莢種質資源的遺傳多樣性研究提供有效的分析方法和手段。李偉等(2017)利用AFLP分子標記技術對10個南方皂莢群體進行了遺傳多樣性分析,結果表明群體內的變異是皂莢遺傳變異的主體,皂莢群體的遺傳多樣性與遺傳結構的形成不僅與其分布廣泛、種子特性及生活史有關,而且與人為的砍伐、引種、生境片段化等因素有重要關系,并由此提出了皂莢的保護策略。就目前而言,分子標記技術在皂莢的系統進化、種質資源鑒定、指紋圖譜構建等方面還缺乏有效利用。

隨著分子生物學的發展,大量的分子標記技術被開發利用,每種分子標記技術都各有其特點。目標起始密碼子多態性(start codon targeted polymorphism, SCoT)是由Collard & Mackill(2009)開發的新目的基因分子標記技術,并首先在水稻上得以應用,是一種基于單引物擴增反應的分子標記技術。SCoT具有操作簡單、通用性良好、多態性豐富、成本低廉等優點,同時又能對性狀進行跟蹤,有利于分子輔助育種,已在多種植物上得到成功應用(楊祥燕等,2013; 王健勝等,2015; 王發明等,2017)。因此,本試驗采用SCoT標記對河南省18份皂莢種質資源進行遺傳多樣性分析,并構建DNA指紋圖譜,以期為皂莢種質資源的鑒定和優良品種的選育等研究提供科學依據。

1材料與方法

1.1 材料和試劑

材料包括3個種共18份,具體見表1。除野皂莢、山皂莢和皂莢T為實生苗外,其余的15份材料是以單株產量高、單刺(或單果)平均重量重的皂莢植株為接穗、野皂莢為砧木(濟科除外)嫁接而獲得的無性系, 其中部分品種已經通過河南省

林木品種審定委員會審定。另外皂莢T是從土耳其引進種子的3年生實生苗,曾被當做皂莢品種,但其所表現出的形態特征與皂莢有較大差異,而與山皂莢極為相似。所有材料均隨機選取3株,每株采集2片健康、幼嫩葉片,置于-80 ℃超低溫冰箱保存備用。

試劑為2×Taq MasterMix(含有Taq DNA Polymerase, 2×Taq PCR Buffer, 3 mmol·L1 MgCl2和400 μmol·L1 dNTP mix)購自北京康為世紀生物科技有限公司,SCoT引物根據蘇亞春等(2012)公布的40條引物序列由金唯智生物科技(北京)有限公司合成。

1.2 皂莢基因組DNA的提取

每個皂莢品種取3片不同植株的幼嫩葉片等量混合,采用改良CTAB法(張安世等,2009)提取皂莢基因組DNA,并將模板DNA濃度稀釋至20 ng·μL1,保存于-20 ℃備用。

1.3 SCoTPCR分析

選用40個SCoT單引物(SC1SC40)和11個由SCoT單引物搭配而成的引物組合共計51個引物對供試材料進行擴增。反應體積 10 μL,其中DNA 1.0 μL,引物 1.4 μL,2×Taq MasterMix 5.0 μL,RNaseFree water 2.6 μL。SCoTPCR擴增程序:首先94 ℃預變性5 min;然后進行40個循環,包括94 ℃變性1 min,50 ℃復性1 min,72 ℃延伸1.5 min。最后72 ℃延伸8 min,4 ℃保存備用。擴增產物用1.5%瓊脂糖凝膠分離。

1.4 數據統計分析

根據電泳結果進行條帶的統計分析。在電泳圖相同位置上若出現DNA條帶記為“1”, 無記為“0”,形成1,0矩陣,將此矩陣導入POPGENE1.32軟件進行多態性百分率(PPL)、觀測等位基因數(Na)、

有效等位基因數(Ne)、Neis基因多樣性(H)和Shannons信息指數(I) 等遺傳多樣性參數分析,參照黃秀等(2014)方法計算引物多態性信息含量(PIC),通過NTSYSpc 2.0軟件依據UPGMA法進行聚類分析,同時參照云天海等(2013)的方法構建18個皂莢種質DNA指紋圖譜。

2結果與分析

2.1 SCoT引物篩選與多態性分析

選用51個SCoT引物對18份皂莢種質材料進行了PCR擴增,最終篩選出了條帶清晰、多態性高、重復性好的15個引物,其中,包括10個SCoT單引物和5個引物組合。所用引物序列及多態性統計結果見表2。15個引物共擴增出226個條帶,其中多態性條帶216個,多態位點百分率(PPL)為96.61%。多態性信息含量(PIC)的變化范圍為0.783 0~0.925 6,平均值為0.875 9。說明所選引物在18個皂莢種質間具有很高的SCoT多態性。觀測等位基因數(Na)的變化范圍為1.812 5~2.000,平均值為1.964 6;有效等位基因數(Ne)的變化范圍為1.319 1~1.591 0,平均值為1.440 1;Neis基因多樣性指數的變化范圍為0.215 4~0.345 0,平均值為0.272 6;Shannons信息指數的變化范圍為0.349 0~0.516 8,平均值為0.426 1。上述結果表明,18個皂莢種質間存在較豐富的遺傳多樣性。同時,利用SPSS17.0軟件對4個主要遺傳多樣性參數多態性信息含量、有效等位基因數、Neis 基因多樣性指數和Shannons信息指數進行非參數Kruskal Wallis Test獨立樣本檢驗,結果顯示18個皂莢種質之間遺傳多樣性水平存在顯著性差異(X2=49.795,P=0.000﹤0.001),表明上述分析得到的各參數值能很好地說明皂莢種質間遺傳多樣性水平。

2.2 遺傳相似性和聚類分析

利用NTSYSpc軟件計算種質間的遺傳相似系數(Gs),結果表明18份皂莢種質兩兩間的Gs在0.491 4~0.938 1之間,平均值為0.711 4,變幅為0.446 7,說明供試材料間存在較大的遺傳差異。其中焦科4和焦科5的Gs最大(0.938 1),親緣關系最近,山皂莢和豫皂2的Gs最小(0.491 4),親緣關系最遠。

利用NTSYSpc軟件對18份皂莢種質進行聚類分析。結果表明(圖2),在Gs為0.60處可將18分皂莢材料分為3組。第1組為野皂莢;第2組為山皂莢和皂莢T;第3組共有15個皂莢種質材料:密刺、碩刺、皂莢H、豫皂1、懷皂王2、博科、焦科1、焦科2、焦科3、濟科、嵩刺1、太行2、豫皂2、焦科4和焦科5。

2.3 DNA指紋圖譜的構建

通過篩選的15個引物對18份皂莢種質的擴增結果分析,選取其中SC1、SC13和SC23 3個引物擴增的8個多態性位點構建了18份皂莢種質的DNA指紋圖譜(圖3)。每份材料都有唯一的指紋圖譜,可以將18份皂莢種質區分并準確鑒定。

3討論與結論

3.1 皂莢種質資源的遺傳多樣性分析

利用分子標記技術分析皂莢種質資源的遺傳多樣性,了解不同皂莢種質之間的遺傳差異,對于皂莢的新品種選育和種質資源鑒定具有重要意義。在遺傳多樣性研究中,多態性是評價遺傳多樣性的重要依據之一。在本研究中,通過SCoT標記對18份皂莢種質資源的遺傳變異分析表明,多態性位點百分率 (PPL)高達96.61%,多態性信息量(PIC)為0.875 9,處于高水平(高:PIC >

0.5; 適中: 0.5 > PIC > 0.25; 低: PIC < 0.25)(徐玉仙等,2015),說明本研究篩選的引物在供試材料間具有很高的多態性,可以有效的用于皂莢的遺傳多樣性分析。同樣,有效等位基因數、Neis 基因多樣性和 Shannons 信息指數分別為1.440 1、0.272 6和0.426 1,也處于較高水平。另外,18個皂莢種質間的遺傳相似系數(Gs)在0.491 4~0.938 1之間,變幅達0.446 7,說明供試材料間存在較豐富的多樣性。通過對多態性信息量、有效等位基因數、Neis 基因多樣性指數和Shannons信息指數,4個主要的遺傳多樣性參數進行的非參數Kruskal Wallis Test獨立樣本檢驗結果也證實了這一結論。

3.2 皂莢種質資源的聚類分析

本研究涉及野皂莢、山皂莢和皂莢共3個種18份材料。聚類結果表明,在Gs為0.60處可將18份皂莢種質分為3組。其中,野皂莢單獨為一組,山皂莢和皂莢T聚為一組,其余15份皂莢種質材料聚為一組。在上述第3組的15份皂莢種質中,焦科4和焦科5是以皂莢莢果為目的選育的品種,其余均以皂刺為目的選育。焦科4和焦科5的接穗均為莢果較大的不同皂莢品種,在所有供試材料中親緣關系最近,在其15份皂莢種質一組聚為一個分支,其它13份材料聚為另一分支,與預期相符。濟科是以嵩刺1為接穗、山皂莢為砧木的嫁接品種,兩者首先聚在一起,也與實際情況相吻合。從聚類圖中還能看出部分種質材料具有一定的地理來源一致性特征。如采集于河南修武的5份材料除焦科4和焦科5聚在一起外,焦科1、焦科2和焦科3也聚在一起;采集于河南博愛的5份材料也分別聚于2個分支,豫皂1、懷皂王2和博科為一個分支,太行2和豫皂2為另一分支;采集于河南嵩縣的3份材料除嵩刺1外,其余2份材料密刺和碩刺聚在一起。另外,皂莢T是從土耳其引進種子的三年生播種苗,起初曾被作為皂莢品種,但其三年生播種苗無論從葉的形態、皮孔大小及形狀以及嫩枝的顏色都與山皂莢極為相似,因此認為該種子在引進時可能鑒定有誤。本研究結果顯示山皂莢和皂莢T聚為一組,兩者具有很近的親緣關系,因此,應將皂莢T歸為山皂莢而非皂莢品種,同時還要結合其它標記技術及其以后的生長過程中所表現出的其它形態學特征進行進一步的驗證。

3.3 皂莢種質資源鑒定與DNA指紋圖譜構建

分子標記是以DNA多態性為基礎的遺傳標記,能直接反映植物間的遺傳差異,具有高效、快捷、準確度高、信息量大等優點,已經成為植物種質鑒定的有效方法,也是數字化身份證構建的有效工具,已在大量的種質資源研究中得到應用。SCoT作為一種新型的DNA分子標記技術,在植物種質鑒定和DNA指紋圖譜構建方面已經得到了初步應用。陳紅等(2014)利用SCoT分子標記技術對71份貴州桃種質進行了分析,結果表明,來自貴陽和黔南州荔波縣的兩份青桃資源(20號和39號)親緣關系較遠,說明它們為同名異物。同時也證實自同一縣份的血桃(31號和32號)和白花桃(16號和19號)不完全是同一樣品,從而在分子水平理清了供試材料的親緣關系,并對部分桃資源名稱混亂的現象進行了糾正,為貴州桃種質資源的科學保存與利用提供了依據。林清等(2013)通過對5個SCoT引物在46份供試芥菜種質中擴增所得的16個多態性DNA譜帶進行統計,初步構建了芥菜種質的指紋圖譜,能夠準確地鑒定這46份芥菜種質。本研究利用SCoT技術,選用3個引物擴增的8個多態性位點構建了18份皂莢種質的DNA指紋圖譜,為皂莢種質間的鑒別提供了重要的科學依據。

綜上所述,通過SCoT標記對18份皂莢種質資源的遺傳變異分析表明,18份皂莢種質材料間具有較豐富遺傳多樣性;聚類分析與皂莢的傳統分類基本吻合,并對可疑材料皂莢T進行了合理歸類;利用3個SCoT引物擴增的8個多態性位點構建的DNA指紋圖譜具有唯一性,可以利用此圖譜對供試材料進行鑒定。

參考文獻:

CHEN H, YANG X, AN HM, 2014. Genetic diversity of peach accessions in Guizhou analysed by SCoT markers [J]. Acta Bot BorealOccident Sin, 34(8):1559-1564. [陳紅,楊鑫,安華明,2014.貴州桃種質資源遺傳多樣性的SCoT分析 [J].西北植物學報,34(8):1559-1564.]

COLLARD BCY,MACKILL DJ, 2009. Start codon targeted (SCoT) polymorphism: a simple,novel DNA marker technique for generating genetargeted markers in plants [J]. Plant Mol Biol Rep, 27:86-93

FAN DC, MA QZ, YANG W, et al,2015. Cultivation techniques of Gleditsia sinensis in the central plains [M].Zhengzhou: The Yellow River Water Conservancy Press:72-83 [范定臣,馬群智,楊偉,等, 2015. 中原地區皂莢栽培技術 [M].鄭州:黃河水利出版社:72-83]

GU WC, LAN YP, SUN CL, 2003. Reseach adnances and utilization development of Dleditsia sinensis in world [J]. Sci Silv Sin,39(4):127-133. [顧萬春,蘭彥平,孫翠玲, 2003. 世界皂莢(屬) 的研究與開發利用 [J].林業科學, 39(4):127-133.]

HUANG X, ZENG J, NIE G, et al,2014. Construction of fingerprinting and genetic diversity analysis of Hemarthia cultivars by ESTSSR [J]. J Trop & Subtrop Bot,22(2):165-171. [黃秀,曾捷,聶剛,等,2014.牛鞭草品ESTSS指紋圖譜構建及遺傳多樣性分析 [J].熱帶亞熱帶植物學報, 22(2):165 -171.]

JIANG JX, ZHANG WM, ZHU LW, et al,2003. Chemical utilization of Chinese Gleditsia sinensis resources [J]. Chin Wild Plant Resourc, 22(6):9-11,26. [蔣建新,張衛明,朱莉偉,等, 2003. 我國皂莢資源的化學利用 [J]. 中國野生植物資源, 22(6): 9-11, 26.]

LAN YP, GU WC, 2006. Geographical variation of morphologic characteristics of Gleditsia sinensis seeds and legumes in the North region [J]. Sci Silv Sin, 42(7):47-51. [蘭彥平,顧萬春, 2006. 北方地區皂莢種子及莢果形態特征的地理變異 [J]. 林業科學, 42(7):47-51.]

LI W, LIN FR, ZHENG YQ, et al, 2013. Phenotypic diversity of pods and seeds in natural populations of Gleditsia sinensis in southern China [J]. Chin J Plant Ecol, 37 (1): 61-69. [李偉,林富榮,鄭勇奇,等, 2013. 皂莢南方天然群體種實表型多樣性 [J]. 植物生態學報, 37(1): 61-69.]

LI W, LIN FR, ZHENG YQ, et al, 2017. Genetic diversity of ten Gleditsia sinersis populations from Southern China [J]. For Res, 30(1):46-52. [李偉,林富榮,鄭勇奇,等, 2017. 10個南方皂莢群體遺傳多樣性的AFLP分析 [J]. 林業科學研究,30(1):46-52.]

LIN Q, LONG ZJ, HAN GH, et al, 2013. Genetic diversity and fingerprint based on SCoT markers in Brassica juncea(L.) Czern. et Coss. [J]. Chin Veget, (12):31-39. [林清,龍治堅,韓國輝,等, 2013. 基于SCoT標記的芥菜種質遺傳多樣性與指紋圖譜 [J]. 中國蔬菜, (12): 31-39.]

SU YC, LING H, WANG HB, et al, 2012.Optimization of SCoTPCR reaction system, and screening and utilization of polymorphic primers in Sugarcane [J]. Chin J Appl Environ Biol,18(5): 810-818. [蘇亞春,凌輝,王恒波,等, 2012. 甘蔗SCoTPCR反應體系優化與多態性引物篩選及應用 [J]. 應用與環境生物學報, 18(5):810-818.]

WANG FM, LI JW, YE KY, et al, 2017. Comparative analysis on the genetic diversity of 41 Vitis germplasm resources by ISSR and SCoT molecular markers [J]. Guihaia, 37(1):1-8. [王發明,李潔維,葉開玉,等,2017. 41份葡萄種質遺傳多樣性的 ISSR和SCoT對比分析 [J].廣西植物,37(1):1-8.]

WANG JH, TANG J, LI D, et al, 2008. Chemical constituents and bioactivity of Gleditsia plants [J]. Chin Wild Plant Resourc, 27(6):1-3. [王薊花,唐靜,李端,等,2008.皂莢化學成分和生物活性的研究進展 [J]. 中國野生植物資源,27(6):1-3.]

WANG JS, HE JH, CHEN HR, et al, 2015. Genetic diversity analysis of pineapple(Ananas comosus(L.)Merr)germplasm with SCoT molecular markers [J]. J Plant Genet Resourc, 16(4): 848-856. [王健勝,賀軍虎,陳華蕊,等, 2015. 菠蘿種質目標起始密碼子(SCoT) 遺傳多樣性分析 [J]. 植物遺傳資源學報, 16(4): 848-856.]

XING JL, MENG YQ, LIN FR, et al ,2017. Development and assessment of ESTSSR for Gleditsia sinersis based on transcriptome sequences [J]. J Plant Genet Resourc, 18(1):149-155. [邢俊連,孟艷瓊,林富榮,等,2017.皂莢ESTSSR分子標記開發與分析評價 [J].植物遺傳資源學報,18(1):149-155.]

XU YX, ZHANG WW, MO HB, et al, 2015. Genetic diversity analysis of Nelumbo accessions based on ESTSSR markers [J]. Plant Divers Resourc, 37(5): 595-604. [徐玉仙,張微微,莫海波,等,2015.基于ESTSSR標記的蓮屬種質資源遺傳多樣性分析 [J].植物分類與資源學報, 37(5): 595-604.]

YANG HD,2003. Multifunction of Gleditsia sinensis and its application in land afforestation [J]. Guizhou Agric Sci, 31(4) : 73-74. [楊海東, 2003.皂莢的多種功效及其綠化應用 [J]. 貴州農業科學, 31(4) : 73-74.]

YANG XY, CAI YB, HUANG QW, et al, 2013. SCoT fingerprints and genetic variations of the Papaya (Carica papaya L.) major cultivars [J]. Acta Bot BorealOccident Sin, 33(9):1756-1761. [楊祥燕,蔡元保,黃秋偉,等, 2013. 番木瓜主栽品種SCoT指紋圖譜構建及遺傳變異分析 [J].西北植物學報,33(9):1756-1761.]

YU LJ, YANG JY, LIU CC, et al, 2016. Extraction and purification of nonionic surfactant from natural plant and the evaluation of their complex synergism [J]. Fine Chem, 33(3):295-301. [余俐佳,楊江勇,劉長城,等, 2016. 從天然植物原料中提取非離子表面活性物質的方法及其復配增效評價 [J].精細化工, 33(3):295-301.]

YUN TH, ZHENG DJ,XIE LS, et al, 2013. Analysis of genetic specificity and DNA fingerprinting establishment for the Hainan Island landraces of Cucurbita moschata [J]. J Plant Genet Resourc, 14(4):679-685. [云天海,鄭道君,謝良商,等, 2013. 中國南瓜海南農家品種間的遺傳特異性分析和 DNA 指紋圖譜構建 [J].植物遺傳資源學報,14(4):679-685.]

ZHANG AS, XING ZF, LIU YY, et al, 2009. Comparison and analysis on the different methods of DNA extration methods from Bryophytes [J]. Henan Sci, 27(5):559-562. [張安世,邢智峰,劉永英,等, 2009. 苔蘚植物DNA不同提取方法的比較分析 [J]. 河南科學, 27(5):559-562.]

ZHANG HL, JIAN LR, HAN CX, et al, 2013. Solation,identification and toxicity of Gleditsiosides O [J]. Acta Bot BorealOccident Sin, 33(6):1234-1238. [張宏利,簡利茹,韓崇選,等, 2013. 皂莢果實中皂苷O的分離鑒定及毒性研究 [J]. 西北植物學報, 33(6): 1234-1238.]

主站蜘蛛池模板: 亚洲av综合网| 成人自拍视频在线观看| 在线播放国产99re| 国产亚洲精品97AA片在线播放| 久久综合伊人 六十路| 伊人久热这里只有精品视频99| 666精品国产精品亚洲| 欧美精品1区2区| 国产高清在线观看| 中文字幕乱码二三区免费| 午夜免费小视频| 亚洲第一极品精品无码| 国产色偷丝袜婷婷无码麻豆制服| AV不卡无码免费一区二区三区| 97超爽成人免费视频在线播放| 看av免费毛片手机播放| 任我操在线视频| 国产黄网永久免费| 亚洲日本韩在线观看| 国产精品19p| 直接黄91麻豆网站| 91久久偷偷做嫩草影院| 无码网站免费观看| 亚洲午夜福利在线| 操国产美女| 激情在线网| 亚洲成年网站在线观看| 国产精品林美惠子在线观看| 国产午夜人做人免费视频中文 | 精品欧美日韩国产日漫一区不卡| 中文字幕波多野不卡一区| 日本欧美精品| 91精品啪在线观看国产| 无码精品一区二区久久久| 亚洲色图综合在线| 亚洲一区免费看| 中文无码精品a∨在线观看| 五月婷婷综合网| 丰满人妻久久中文字幕| 国产日韩欧美视频| 亚洲一区精品视频在线| 天天婬欲婬香婬色婬视频播放| 亚洲二区视频| 国产精品视频导航| 国产在线精彩视频二区| 国产内射在线观看| 国产精品冒白浆免费视频| 一区二区三区四区精品视频| 国产精品久久久久久影院| 91青青视频| 国产在线视频导航| 无码视频国产精品一区二区| 91国语视频| 国产综合欧美| 波多野结衣爽到高潮漏水大喷| 国产精品护士| 欧美成人手机在线视频| 中文字幕永久在线看| 久草热视频在线| 香蕉综合在线视频91| 在线看片中文字幕| 国产精品福利社| 亚国产欧美在线人成| 99久久性生片| 国产成人一级| 欧美亚洲国产日韩电影在线| 久久99蜜桃精品久久久久小说| 色综合天天娱乐综合网| 高清免费毛片| 久久精品波多野结衣| 无码在线激情片| 无码综合天天久久综合网| 久久亚洲中文字幕精品一区| 毛片免费高清免费| 在线a视频免费观看| 日韩一级毛一欧美一国产| 亚洲免费福利视频| 精品国产三级在线观看| 中文字幕66页| 亚洲中文字幕久久精品无码一区| 国产区在线观看视频| 国产成人综合网|