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

基于表型性狀的新疆大果沙棗多樣性分析及種下等級(jí)的劃分

2015-05-06 16:35:06翟申修黃俊華鄭堯宋丹華
天津農(nóng)業(yè)科學(xué) 2014年9期

翟申修 黃俊華 鄭堯 宋丹華

摘 要:為劃分出新疆地區(qū)大果沙棗的栽培變種,選擇了具代表性的5個(gè)地區(qū)的38棵大果沙棗為研究對(duì)象,運(yùn)用統(tǒng)計(jì)分析方法對(duì)大果沙棗的15個(gè)表型性狀進(jìn)行了分析與比較。結(jié)果表明:大果沙棗6個(gè)質(zhì)量性狀Shannon指數(shù)的平均值是0.966,Simpson 指數(shù)的平均值是0.410,大果沙棗的質(zhì)量性狀豐富度較高;9個(gè)數(shù)量性狀變異系數(shù)平均值為15.426%。通過(guò)主成分分析,前4個(gè)主成分的累計(jì)貢獻(xiàn)率為71.149%,果實(shí)、果核大小、果形指數(shù)及千粒質(zhì)量等9個(gè)性狀可作為劃分種下等級(jí)的主要依據(jù)。聚類(lèi)分析后可將38棵大果沙棗植株劃分為4個(gè)栽培變種。

關(guān)鍵詞:大果沙棗;表型性狀;栽培變種

中圖分類(lèi)號(hào):S793.9 文獻(xiàn)標(biāo)識(shí)碼:A DOI 編碼:10.3969/j.issn.1006-6500.2014.09.025

Abstract: For purpose of classifying infraspecific taxon of E.moorcroftii Wall.Ex Schlecht. in Xinjiang,38 accessions of germplasm resources from five representative were as of Xinjiang were choosed as subjects. With the method of statistical analysis 15 indicators of fruit and kernel of E.moorcroftii Wall.Ex Schlecht were compared and analyzed.The results showed that the mean value of 6 quality traits of the Shannon index in E.moorcroftii Wall.Ex Schlecht.was 0.966,the mean value of its Simpson index was 0.410,the quality traits of E.moorcroftii Wall.Ex Schlecht.were richness;the average of the coefficient of variation of 9 quantitative traits was 15.426%. Through principal component analysis, the cumulative contribution of the first 4 principal components was 71.149%,9 characters could be used as the main basis for dividing a infraspecific taxon.The cluster analysis of 38 individual plants could be divided into 4 groups.

Key words: E.moorcroftii Wall.Ex Schlecht.; phenotypic traits; cultivar

胡頹子科胡頹子屬(Elaeagnus L.)植物約有80多種,廣布于亞洲東部和東南部的亞熱帶和溫帶地區(qū),少數(shù)種分布于亞洲其他地區(qū)及歐洲溫帶地區(qū),北美也有見(jiàn)分布。在我國(guó)約有55種,廣布于全國(guó)各地,但長(zhǎng)江流域及以南地區(qū)更為常見(jiàn)。我國(guó)西北新疆地區(qū)分布的胡頹子屬植物全部屬于落葉組,其中,大果沙棗(E.moorcroftii Wall.Ex Schlecht.)在《中國(guó)沙漠植物志》中被記載[1]。在我國(guó)西北地區(qū),胡頹子屬落葉組植物常被俗稱(chēng)為沙棗。而胡頹子科胡頹子屬落葉組植物多具有抗旱、抗風(fēng)沙、耐鹽堿、耐貧瘠等特點(diǎn),并且是良好的藥用資源,在干旱及鹽堿地區(qū)擁有很好的發(fā)展前景[2-17]。其中,大果沙棗實(shí)品質(zhì)較高,是良好的蜜源植物,也是南疆防風(fēng)固沙林、護(hù)田林以及四旁綠化的主要樹(shù)種[3-5],具有較高的利用價(jià)值,應(yīng)用前景較好。

對(duì)大果沙棗表型性狀的多樣性研究,不僅有助于了解大果沙棗各變異特征對(duì)區(qū)域環(huán)境變化的適應(yīng)能力,而且還可以為群落分類(lèi)、群落結(jié)構(gòu)特征的揭示等提供重要的依據(jù)[6]。根據(jù)初步調(diào)查,大果沙棗個(gè)體間果實(shí)千粒質(zhì)量已存在較大差異。為了選育品質(zhì)優(yōu)良的大果沙棗,篩選了部分表型性狀作為其種下等級(jí)的劃分依據(jù),探究其栽培變種的分類(lèi)方式,并最終劃分出確切的種下等級(jí)。目前,國(guó)內(nèi)文獻(xiàn)中對(duì)胡頹子屬落葉組種的劃分依據(jù)以相對(duì)較穩(wěn)定的表型性狀為主 [2,11-20],因此筆者也選擇表型性狀作為大果沙棗種下等級(jí)的劃分依據(jù)。種下等級(jí)的劃分有助于篩選大果沙棗中品質(zhì)優(yōu)良的栽培變種,對(duì)種質(zhì)資源的開(kāi)發(fā)利用具有重要意義[7-10],可為今后的進(jìn)一步研究奠定基礎(chǔ)。

1 材料和方法

1.1 試驗(yàn)材料采集

本研究采集了全疆范圍內(nèi)6個(gè)地區(qū)的胡頹子屬落葉組植物,篩選了大果沙棗分布較廣的新疆南北部以及東部地區(qū),總體來(lái)說(shuō),大果沙棗的采樣點(diǎn)涵蓋了新疆的5個(gè)地區(qū)(圖1)。對(duì)比《中國(guó)沙漠植物志》、《中國(guó)植物志》[1,11]中對(duì)沙棗和尖果沙棗的描述,以及文獻(xiàn)[12-14]中對(duì)東方沙棗的描述,大果沙棗與其他種之間的主要差別在于花盤(pán)有毛,果核長(zhǎng)1.5~1.8 cm以及枝少具刺或無(wú)刺等3個(gè)性狀,可作為區(qū)分大果沙棗的主要依據(jù),此外因取樣的廣泛性需要,將果核長(zhǎng)大于1.8 cm且其他性狀均符合的大果沙棗植株也選為研究對(duì)象。新疆阿勒泰、阿克蘇、喀什、吐魯番及哈密5個(gè)地區(qū)的共計(jì)38株大果沙棗種實(shí)樣本均采自防護(hù)林及住宅周邊地區(qū),屬于栽培品種;目前新疆地區(qū)所見(jiàn)野生種均為尖果沙棗。采集取樣分為2個(gè)時(shí)期,每年10月份采集果期標(biāo)本,次年采集同株花期標(biāo)本。

1.2 性狀測(cè)定方法

抽樣測(cè)定表型性狀為枝葉、果實(shí)、果核及花的15個(gè)形態(tài)性狀:枝有無(wú)刺、果實(shí)縱徑(cm)、果實(shí)橫徑(cm)、果實(shí)果形指數(shù)、果柄長(zhǎng)(cm)、果面顏色、果實(shí)形狀、覆蓋鱗片程度、果實(shí)千粒質(zhì)量(g)、果核縱徑(cm)、果核橫徑(cm)、果核果形指數(shù)、果核千粒質(zhì)量(g)、柱頭彎曲形狀、柱頭與雄蕊對(duì)比。葉面通過(guò)解剖鏡觀察,數(shù)量化性狀果柄長(zhǎng)、果實(shí)縱徑、果實(shí)橫徑、果核橫徑及果核縱徑等用直尺、游標(biāo)卡尺測(cè)定;千粒質(zhì)量用電子天平測(cè)定。根據(jù)果核長(zhǎng)和果核寬計(jì)算出果核果形指數(shù);根據(jù)果實(shí)長(zhǎng)和果實(shí)寬計(jì)算出果實(shí)果形指數(shù)。質(zhì)量性狀制定統(tǒng)一標(biāo)準(zhǔn)后,通過(guò)1,2,3…等數(shù)字表示質(zhì)量性狀的不同水平,例如:枝有刺標(biāo)記為1,枝無(wú)刺標(biāo)記為2。另外,在測(cè)定柱頭彎曲形狀過(guò)程中發(fā)現(xiàn),大果沙棗同株上所表現(xiàn)的性狀并不穩(wěn)定,往往呈現(xiàn)“J”狀與環(huán)狀同時(shí)存在的情況,故研究中采用其同株中出現(xiàn)次數(shù)較多的柱頭形狀作為試驗(yàn)數(shù)據(jù)。每單株測(cè)定花及果實(shí)樣本各10個(gè),每個(gè)樣本重復(fù)測(cè)定3次,取平均值[15]。

1.3 統(tǒng)計(jì)分析方法

利用Shannon指數(shù)和Simpson指數(shù)作為大果沙棗質(zhì)量性狀多樣性水平的評(píng)價(jià)指標(biāo)。

用變異系數(shù)CV表示性狀離散程度[18]。

本文所用統(tǒng)計(jì)分析軟件為SPSS-statistics-v19。

2 結(jié)果與分析

2.1 質(zhì)量性狀多樣性分析

由表1數(shù)據(jù)可知,6個(gè)質(zhì)量性狀Shannon指數(shù)的平均值是0.966,最高的是覆蓋鱗片程度,達(dá)到 1.609,最低的是枝有無(wú)刺,僅為0.562;從不同來(lái)源性狀來(lái)分析,果實(shí)>花部>枝,Shannon指數(shù)均值分別為1.238,0.761和0.562;Shannon指數(shù)大于1.000的性狀有果面顏色和鱗片覆蓋程度,表明這兩個(gè)性狀的多樣性程度相對(duì)較高[13]。而Simpson 指數(shù)的平均值是0.410,最高的也是覆蓋鱗片程度,為0.611,最低的也是枝有無(wú)刺,為0.229; 從不同來(lái)源性狀來(lái)分析,果實(shí)>花部>枝,Simpson指數(shù)均值分別為 0.535,0.313 和 0.229;Simpson指數(shù)大于0.500的性狀也是果面顏色和鱗片覆蓋程度。結(jié)果表明,大果沙棗質(zhì)量性狀變異中果實(shí)的多樣性程度較高,花與枝的質(zhì)量性狀相對(duì)穩(wěn)定。

2.2 數(shù)量性狀多樣性分析

2.2.1 變異系數(shù) 由表2可知,大果沙棗的數(shù)量性狀變異系數(shù)平均值為15.426%,其中,果實(shí)及果核的千粒質(zhì)量變異系數(shù)最高,分別達(dá)到了25.730%和22.974%;果實(shí)千粒質(zhì)量最小值為525.000 g,最大值為1 562.000 g,變化幅度較大;果柄長(zhǎng)的變異系數(shù)為19.999%,相對(duì)來(lái)說(shuō)也較高[16];而果實(shí)縱橫徑分別為12.687%和13.549%,果核縱橫徑分別為9.092%和13.119%,說(shuō)明果實(shí)、果核差異不大。其余數(shù)量性狀的變異系數(shù)則分布在9.000%~14.000%之間,總體來(lái)說(shuō)大果沙棗種質(zhì)資源中數(shù)量性狀存在一定的變異,因此,在所選的大果沙棗品種之中,可能已經(jīng)產(chǎn)生了栽培變種。

2.2.2 均值、標(biāo)準(zhǔn)差與多重比較 從表型性狀的多重比較來(lái)看,不同地區(qū)的大果沙棗果實(shí)與果核縱橫徑、果形指數(shù)、果核千粒質(zhì)量等性狀呈現(xiàn)差異較大的狀況,哈密地區(qū)的果實(shí)千粒質(zhì)量平均值最高,達(dá)到1 354.67 g,阿勒泰最低,為634.50 g;喀什地區(qū)的果實(shí)及果核果形指數(shù)最低,而阿克蘇地區(qū)果形指數(shù)則較高。綜上可見(jiàn),大果沙棗的數(shù)量性狀已存在差異,可對(duì)其種下等級(jí)的劃分提供依據(jù)。

2.3 性狀的主成分分析

通過(guò)主成分分析(表4)得出,前4個(gè)主成分的貢獻(xiàn)率分別為26.080%,18.962%,15.944%,10.162%,主成分累計(jì)貢獻(xiàn)率為71.149%,前4個(gè)主成分可以代表原始因子的大部分信息[17]。根據(jù)入選主成分得分系數(shù)得出,第1個(gè)主成分主要受果實(shí)及果核大小與千粒質(zhì)量的影響;而第2個(gè)主成分主要受枝的特征以及果實(shí)形狀影響;第3主成分受果實(shí)與果核果形指數(shù)的影響;第4主成分主要受柱頭彎曲形狀的影響。總體看來(lái),果實(shí)與果核大小及質(zhì)量、花的特征等性狀是影響力較大的性狀。

2.4 聚類(lèi)分析

在SPSS軟件系統(tǒng)中,采用最長(zhǎng)距離法(furthest neighbor method)進(jìn)行聚類(lèi)分析。

根據(jù)需要,依據(jù)分析結(jié)果將所采大果沙棗擬分為4個(gè)栽培變種(圖2)。栽培變種1包括HM5,HM30,HM34,……,KSH42,KSH43等共計(jì)13株,該栽培變種主要特點(diǎn)是:果實(shí)較大,果實(shí)千粒質(zhì)量在1 404~1 562 g之間,平均值為1 493 g,柱頭多長(zhǎng)于雄蕊,果實(shí)形狀多為卵形。栽培變種2包括KSH2,KSH13,KSH23,……,HM33,HM35等13株,栽培變種2果實(shí)形狀橢圓形居多,果實(shí)千粒質(zhì)量在1 110~1 392 g之間,平均值為1 247 g,柱頭全部長(zhǎng)于雄蕊,柱頭多數(shù)環(huán)狀。栽培變種3包括KSH35,KSH44,KSH49,……,TLF7等7株,栽培變種3主要特征為:果實(shí)千粒質(zhì)量在525~767 g之間,平均值為677 g,柱頭全部長(zhǎng)于雄蕊,柱頭彎曲形狀鉤狀與環(huán)狀均有,果面顏色棕紅色居多。栽培變種4包括KSH50,KSH51,HM36,……,TLF3等5株,栽培變種4主要特征為:柱頭彎曲形狀以“J”狀居多,果實(shí)千粒質(zhì)量在930~1 037 g之間,平均值為980 g,果實(shí)與果核果形指數(shù)較高,果實(shí)果形指數(shù)平均1.48,果核果形指數(shù)平均3.74。

現(xiàn)根據(jù)聚類(lèi)結(jié)果制定其植物檢索表,如下。

1.果實(shí)千粒質(zhì)量1 100~1 570 g,鱗片覆蓋程度低。

2.果實(shí)多呈卵形,柱頭長(zhǎng)度部分短于雄蕊,果核橢圓形 栽培變種1

2.果實(shí)多呈橢圓形,柱頭全部長(zhǎng)于雄蕊,果核窄橢圓形 栽培變種2

1.果實(shí)千粒質(zhì)量500~1 100 g,鱗片覆蓋程度中等

3.枝少具刺或無(wú)刺,果核果形指數(shù)2.8~3.5 栽培變種3

3.枝無(wú)刺,果核果形指數(shù)3.5~4.2

栽培變種4

3 結(jié)論與討論

3.1 對(duì)大果沙棗的界定條件需拓展

本研究以《中國(guó)沙漠植物志》中對(duì)大果沙棗的描述作為參考,篩選了大果沙棗植株作為研究材料,采樣過(guò)程中發(fā)現(xiàn),目前新疆地區(qū)的大果沙棗栽培品種中已出現(xiàn)部分個(gè)體的果核縱徑不能滿足1.5~1.8 cm的條件,通過(guò)對(duì)比于瑋瑋等[19]在《大果沙棗和尖果沙棗植物學(xué)特征比較研究》一文中對(duì)大果沙棗的描述,本研究將果核縱徑大于1.8 cm但其他性狀均相符的植株也列入選擇范圍,根據(jù)其文中對(duì)大果沙棗植物學(xué)特征分析,可確定本研究所選植物材料均為大果沙棗,故《中國(guó)沙漠植物志》中所述大果沙棗的果核縱徑應(yīng)包含大于1.8 cm的個(gè)體。

3.2 大果沙棗的表型性狀豐富

Shannon指數(shù)和Simpson指數(shù)越高,代表性狀的豐富程度和均勻程度越大。根據(jù)吳根松等[20]對(duì)梅花表型多樣性研究所述,Shannon指數(shù)和Simpson 指數(shù)分別達(dá)到1.000和0.500即為多樣性程度高。大果沙棗6個(gè)質(zhì)量性狀Shannon指數(shù)的平均值是0.966,最高的是覆蓋鱗片程度,達(dá)到 1.609,說(shuō)明它的質(zhì)量性狀已經(jīng)具有較高的豐富度;而性狀Simpson 指數(shù)的平均值是0.410,最高的也是覆蓋鱗片程度,為0.611,代表其質(zhì)量性狀的多樣度也較大。從不同來(lái)源性狀的Simpson指數(shù)均值來(lái)分析,果實(shí)>花部>枝,表明果實(shí)部分?jǐn)?shù)量性狀的多樣性程度最高。大果沙棗數(shù)量性狀的變異系數(shù)平均值為15.426%,果實(shí)及果核的千粒質(zhì)量變異系數(shù)最高,分別達(dá)到了25.730%和22.974%,說(shuō)明大果沙棗數(shù)量性狀也出現(xiàn)了一定變異,其中果實(shí)及果核千粒質(zhì)量變異程度相對(duì)較大。綜上分析,大果沙棗質(zhì)量性狀與數(shù)量性狀均產(chǎn)生了豐富的變異。

3.3 所采大果沙棗可劃分為4個(gè)栽培變種

通過(guò)主成分分析,抽取大果沙棗的9種性狀作為分類(lèi)依據(jù),通過(guò)變異較大的果實(shí)千粒質(zhì)量、鱗片覆蓋程度、果實(shí)形狀、柱頭與雄蕊比、枝有無(wú)刺、果核果形指數(shù)等性狀,最終可將大果沙棗劃分為4個(gè)栽培變種:栽培變種1和2的果實(shí)千粒質(zhì)量為1 100~1 570 g,鱗片覆蓋程度低,而前者果實(shí)多呈卵形,柱頭長(zhǎng)度部分短于雄蕊,果核橢圓形,后者果實(shí)多呈橢圓形,柱頭全部長(zhǎng)于雄蕊,果核窄橢圓形;栽培變種3和4的果實(shí)千粒質(zhì)量為500~1 100 g,鱗片覆蓋程度中等,其中,栽培變種3枝少具刺或無(wú)刺,果核果形指數(shù)為2.8~3.5,栽培變種4則枝無(wú)刺,果核果形指數(shù)介于3.5~4.2之間。其中,大果沙棗的果實(shí)縱橫徑分別為12.687%和13.549%,果核縱橫徑分別為9.092%和13.119%,低于平均水平,差異并不明顯,不能作為區(qū)分栽培變種的依據(jù)。聚類(lèi)分析所劃分的結(jié)果表明,在大果沙棗的果實(shí)千粒質(zhì)量甚至果面顏色等性狀中已有明顯差異,存在于栽培變種間的性狀變異對(duì)本研究栽培變種的劃分具有極其重要的意義,因?yàn)榉植荚陬?lèi)型間的變異不僅能夠反映地理、生殖隔離上的變異,而且它還是種內(nèi)多樣性的重要組成部分[21]。因此,以大果沙棗中變異較大的性狀為依據(jù),可劃分出新疆地區(qū)大果沙棗的栽培變種。

3.4 果實(shí)千粒質(zhì)量栽培變種間差異大

大果沙棗可被分為4個(gè)栽培變種,而對(duì)其栽培變種間表型性狀進(jìn)行比對(duì)后發(fā)現(xiàn),其果實(shí)千粒質(zhì)量差異較大,變異幅度在500~1 570 g之間。栽培變種1果實(shí)較大,果實(shí)千粒質(zhì)量在1 404~1 562 g之間,平均值為1 493 g;栽培變種2果實(shí)千粒質(zhì)量在1 110~1 392 g之間,平均值為1 247 g;栽培變種3果實(shí)千粒質(zhì)量在525~767 g之間,平均值為677 g;栽培變種4果實(shí)千粒質(zhì)量在930~1 037 g之間,平均值為980 g。結(jié)果表明,果實(shí)千粒質(zhì)量在大果沙棗栽培變種間存在較大變異。

參考文獻(xiàn):

[1] 中國(guó)科學(xué)院蘭州沙漠研究所.中國(guó)沙漠植物志(第2卷)[M].北京:科學(xué)出版社,1987:392-395.

[2] 黃俊華,買(mǎi)買(mǎi)提江.新疆胡頹子屬植物(Elaeagnus)分類(lèi)探討[J].植物研究,2005,25(3):268-271.

[3] 江發(fā)壽,謝建新,劉金榮,等.沙棗的營(yíng)養(yǎng)成分分析及沙棗油的理化常數(shù)測(cè)定[J]. 石河子大學(xué)學(xué)報(bào):自然科學(xué)版,2002,6(1):20-22.

[4] 李磊,賈志清,寧虎森,等.水分脅迫下2種沙棗的抗旱性比較[J].林業(yè)科學(xué)研究,2009,22(3):335-342.

[5] 齊曼·尤努斯,李秀霞,李陽(yáng),等.鹽脅迫對(duì)大果沙棗膜脂過(guò)氧化和保護(hù)酶活性的影響[J].干旱區(qū)研究,2005,22(4):87-91.

[6] 趙紅洋,李玉霖,王新源,等.科爾沁沙地52種植物葉片性狀變異特征研究[J].中國(guó)沙漠,2010(6):1 292-1 298.

[7] 石翠芳,孫智達(dá),謝筆鈞.沙棗果肉原花青素的提取、純化及抗氧化性能的研究[J].農(nóng)業(yè)工程學(xué)報(bào),2005,22(3):158-161.

[8] 孫景寬,張文輝,劉新成.干旱脅迫對(duì)沙棗和孩兒拳頭的生理特性的影響[J].西北植物學(xué)報(bào),2008,28(9):1 868-1 874.

[9] 辛艷偉.沙棗的開(kāi)發(fā)和利用[J].安徽農(nóng)業(yè)科學(xué),2007(2):399-400,402.

[10] 聶小蘭.沙棗的研究現(xiàn)狀及展望[J].北方園藝,2007(4):67-69.

[11] 中國(guó)科學(xué)院中國(guó)植物志編輯委員會(huì).中國(guó)植物志[M].北京:科學(xué)出版社,1983:41-42.

[12] 蘇維埃社會(huì)主義共和國(guó)科馬洛夫植物研究所.蘇聯(lián)植物原材料(第二冊(cè))[M].圣彼得堡:蘇聯(lián)科學(xué)院,1957:97-98.

[13] 蘇維埃社會(huì)主義共和國(guó)科馬洛夫植物研究所.蘇聯(lián)植物志(第十五冊(cè))[M].圣彼得堡:蘇聯(lián)科學(xué)院,1949:518-525.

[14]塔吉克蘇維埃社會(huì)主義共和國(guó)科學(xué)院植物研究所.塔吉克蘇維埃社會(huì)主義共和國(guó)植物區(qū)系(第六冊(cè))[M].列寧格勒:科學(xué)出版社,1981:9-671.

[15] 王金星,潘剛,王滑,等.西藏核桃葉片和堅(jiān)果表型多樣性及其相關(guān)關(guān)系研究[J].林業(yè)科學(xué)研究,2012(2):236-240.

[16] 肖亮,蔣建雄,易自力,等.廣西省芒野生居群表型多樣性研究[J]. 草業(yè)學(xué)報(bào),2013(4):43-50.

[17] 周蓉,蔣芳玲,梁梅,等.基于表型性狀的番茄品種評(píng)價(jià)和遺傳多樣性分析[J].西北農(nóng)業(yè)學(xué)報(bào),2012(9):95-102.

[18] 楊繼.植物種內(nèi)形態(tài)變異的機(jī)制及其研究方法[J].武漢植物學(xué)研究,1991,9(2):185-195.

[19] 于瑋瑋,王小莉,李慧,等.大果沙棗和尖果沙棗植物學(xué)特征比較研究[J].天津農(nóng)學(xué)院學(xué)報(bào),2012(4):36-38,50.

[20] 吳根松,孫麗丹,郝瑞杰,等.梅花種質(zhì)資源表型多樣性研究[J].安徽農(nóng)業(yè)科學(xué),2011,20:12 008-12 009,12 012.

[21] Daniel L H,Andrew G C.Principles of population genetics[M]. Sunderland: Sinauer Associates,Inc.,USA,1989.

主站蜘蛛池模板: 欧美成人综合在线| h视频在线播放| 亚洲嫩模喷白浆| 综合社区亚洲熟妇p| 亚洲男人的天堂在线| 亚洲综合片| a毛片在线| 午夜电影在线观看国产1区| 亚洲天堂区| 中文无码影院| 精品国产欧美精品v| 91在线丝袜| 国产成人在线小视频| 99热6这里只有精品| 久久人搡人人玩人妻精品一| 国产精品久久久久鬼色| 激情视频综合网| 幺女国产一级毛片| 日韩一级二级三级| 免费激情网站| 伊人AV天堂| 国产h视频在线观看视频| 99re这里只有国产中文精品国产精品 | 日韩二区三区无| 色噜噜狠狠狠综合曰曰曰| 国产成人精品一区二区免费看京| 国产成人久久综合777777麻豆| 青青青国产在线播放| 韩国v欧美v亚洲v日本v| 蜜臀AVWWW国产天堂| 精品一区国产精品| 97视频免费在线观看| 99er这里只有精品| 日韩高清中文字幕| 在线人成精品免费视频| 丁香婷婷激情综合激情| 亚洲成人精品久久| 日韩免费毛片| 精品视频一区二区观看| 国产在线98福利播放视频免费| 国产亚洲欧美在线专区| 国产极品美女在线播放| 亚洲三级网站| 国产精品毛片一区| 国产主播福利在线观看| 久久无码高潮喷水| 成人国产精品网站在线看 | 国产视频入口| 精品国产一区二区三区在线观看| 亚洲成人高清无码| 精品无码日韩国产不卡av| 色妞www精品视频一级下载| 欧美午夜性视频| 亚洲第一成年网| 91视频区| YW尤物AV无码国产在线观看| 婷婷六月综合| 婷婷午夜天| 亚洲欧美在线综合图区| 亚洲日韩每日更新| 免费无码AV片在线观看国产| 亚洲黄色视频在线观看一区| 亚洲精品少妇熟女| 色综合天天综合| 国产一级精品毛片基地| 青青草国产在线视频| 新SSS无码手机在线观看| 国产精品lululu在线观看 | 亚洲最猛黑人xxxx黑人猛交| 992Tv视频国产精品| 精品无码人妻一区二区| 国产在线91在线电影| 青青青伊人色综合久久| 国产成人欧美| 亚洲视频免费在线看| 手机永久AV在线播放| 日本精品中文字幕在线不卡| 114级毛片免费观看| 欧美日本一区二区三区免费| 国产亚洲美日韩AV中文字幕无码成人| 亚洲中文字幕在线精品一区| 婷婷六月激情综合一区|