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

我國野牡丹屬植物的表型多樣性研究

2014-04-29 17:10:08戴小紅等
熱帶作物學報 2014年10期

戴小紅等

摘 要 對國產野牡丹屬6種16個居群的22個表型性狀進行觀測,通過方差分析、聚類分析、因子分析對其表型變異、種間親緣關系等進行探討,以揭示野牡丹屬植物的表型多樣性,并為其提供分類依據。結果表明:野牡丹屬種間與種內居群間在數量性狀上均存在豐富的表型變異,15個數量性狀的變異系數為6.77%~34.50%,其中花序花朵數、小苞片寬、葉柄長、花梗長的變異較大,果徑、花徑、花萼長、花瓣長和葉長寬比的變異較小;質量性狀的多樣性主要存在于種間,種內居群間在葉柄顏色、花色上存在差異。表型聚類結果顯示當平均距離為34時,16個居群可分為3類,展毛野牡丹、多花野牡丹、細葉野牡丹和野牡丹為1類,毛菍、地菍各為1類,同種下各居群沒有嚴格依地理距離聚類。因子分析結果顯示,花部形態變異是野牡丹屬植物表型分化的主要因素,基出脈、花梗長、莖被毛的變異對其表型分化的作用較小,在進行表型分類時可避免次要性狀的干擾。

關鍵詞 野牡丹屬;表型多樣性;聚類分析;因子分析

中圖分類號 Q944 文獻標識碼 A

Abstract In order to reveal the phenotypic diversity of Melastoma L. and provide classification evidence, twenty-two phenotypic traits of six species and sixteen populations of Melastoma L. in China were investigated. Phenotypic variations between and within species as well as phylogenetic relationship between species were discussed by using variance analysis, cluster analysis and factor analysis. The results showed that there were rich phenotypic variations between and within species of Melastoma L. in quantitative traits. The variation coefficients of 15 quantitative traits ranged from 6.77% to 34.50%. Among these traits, flower number of inflorescence, bractlet width, petiole length and stalk length varied more while fruit diameter, flower diameter, calyx length, petal length and length-width ratio of leaf varied less. The diversity of qualitative traits mainly existed between species, and colors of petiole and flower also varied in populations. According to cluster analysis by twenty-two phenotypic traits, 16 populations of Melstoma L. could be roughly divided into three groups when the average distance was about 34. Besides M. normale, M. affine, M. intermedium and M. candidum gathering in one group, M. sanguineum and M. dodecandrum was divided into two groups respectively. Populations within species did not clustered strictly by geographical distances. Factor analysis indicated that floral morphological variations might be the main causes of phenotypic differentiation in Melastoma L., on which the variations of basal vein, stalk length and pubescence on stem effected lightly. Therefore, some secondary phenotypic traits of Melastoma L. could be paid little attention in morphological classification.

Key words Melastoma; Phenotypic diversity; Cluster analysis; Factor analysis

doi 10.3969/j.issn.1000-2561.2014.10.026

表型是基因表達與所處環境交互作用的結果,是生物遺傳受制于生態環境條件的遺傳表征,表型特征的不同在很大程度上能夠反映同種植物之間的相似性以及不同種植物之間的本質差異,表型多樣性是基因組DNA編碼序列與其調控序列多樣性的具體表現,因此是遺傳多樣性研究最基本的組成部分[1-2]。與分子標記的遺傳多樣性研究相比,表型是對完整基因系統或基因型的全面度量,利用表型多樣性差異研究群體的遺傳多樣性具有簡便易行、快速經濟的優點[3-4]。作為種質資源研究的重要組成部分,表型多樣性研究對種質資源的鑒定與保存、遺傳改良及開發利用均具有重要意義。

野牡丹屬(Melastoma L.)是一個典型的熱帶亞洲至熱帶澳洲分布式樣的屬,全屬約100種,分布中心在東南亞,中國有9種1變種,多見于長江流域以南各省(區),以廣東、廣西種數最多[5-6]。該屬植物富含鞣質、酚類和黃酮類成分,具有自由基清除、抗氧化、抗癌等多種生物活性[7-10],藥用價值較高。野牡丹屬也是野牡丹科(Melastomataceae)中最具觀賞價值的一個屬[11],屬下植物類型多樣、花色艷麗、花期長,是一類頗具華南鄉土特色的園林植物,部分種的果實還可供提取色素[12]或食用,具備良好的開發利用前景。

近年來,不少學者針對野牡丹屬植物開展了資源調查[11,13-15]、栽培繁殖[16-19]、傳粉生物學[20-21]、生理特性[22-23]等一系列研究。長期以來,學術界對于該屬植物的分類存在一定爭議,部分學者從種子、花粉、葉片的表面特征與顯微形態[24-26]等方面對野牡丹屬植物進行了考察,為該屬植物的分類提供了一定的科學依據。但迄今為止還未見有關該屬植物表型多樣性的報道。本研究以6種16個居群的野牡丹屬植物為對象開展表型多樣性調查與分析,了解其表型變異特征,旨在為該屬植物的分類、資源的保存與利用等提供依據。

1 材料與方法

1.1 供試材料

供試材料為6種野牡丹屬植物,共16個居群,材料及來源見表1。

1.2 方法

1.2.1 數量性狀觀測 調查各居群成熟期植株的葉長、葉寬、葉長寬比、葉柄長、花序花朵數、花梗長、小苞片長、小苞片寬、花萼長、萼裂片長、萼裂片寬、花徑、花瓣長、花瓣寬、果徑等15個數量性狀。根據群體大小,每個居群隨機挑選3~10個個體進行測定,每個個體選取5片成熟葉,用游標卡尺測量葉長、葉寬、葉柄長,計算葉長寬比;每個居群取5~10個花序統計花朵數,并隨機選擇10~20朵花測量花梗長、小苞片長、小苞片寬、花萼長、萼裂片長、萼裂片寬、花徑、花瓣長和花瓣寬;每個居群選取5~30個果實進行果徑測量。

1.2.2 質量性狀觀測 選取7個質量性狀進行觀測,為便于數量化和統計分析,對各質量性狀進行賦值(表2)。

1.3 統計分析

以上述野牡丹屬植物表型性狀的觀測結果為原始數據,在SAS9.3軟件中進行方差分析、相關分析、聚類分析和因子分析。其中方差分析采用GLM程序,以變異系數(CV)表示性狀的離散程度;計算數量性狀間的Pearson相關系數,考察其相關關系;采用類平均法(Average method)進行聚類分析,建立聚類樹狀圖;因子分析采用的提取方法為主成分分析法(Principal component analysis),旋轉方法為最大方差法(VARIMAX),以確定造成表型差異的主要性狀因子。

2 結果與分析

2.1 數量性狀多樣性

對野牡丹屬植物的15個數量性狀進行方差分析的結果顯示(表3),野牡丹屬植物種間和部分種內居群間在葉長、葉寬、葉長寬比、葉柄長、花序花朵數、小苞片長、小苞片寬、花梗長、花萼長、萼裂片長、萼裂片寬、花徑、花瓣長、花瓣寬、果徑這15個主要數量性狀上存在極顯著差異,表現出豐富的表型多樣性。15個數量性狀中,極差較大的性狀包括葉長、花徑、葉寬;極差較小的性狀為葉長寬比、萼裂片寬。極差在一定程度上反映了數量性狀的離散程度,但易受極端值的影響。從變異系數來看,15個數量性狀的變異系數幅度為6.77%~34.50%,各數量性狀按其變異系數由大到小為:花序花朵數>小苞片寬>葉柄長>花梗長>小苞片長>葉寬>萼裂片寬>葉長>花瓣寬>萼裂片長>葉長寬比>花瓣長>花萼長>花徑>果徑,可見各性狀的離散程度不同,變異系數越大,數據的離散程度越高,性狀越不穩定。由此可知,野牡丹屬植物的花序花朵數、小苞片寬、葉柄長、花梗長的變異幅度較大,多樣性較豐富,果徑、花徑、花萼長、花瓣長、葉長寬比的變異幅度較小,性狀較穩定。

2.2 數量性狀間的相關性分析

對15個數量性狀進行相關分析的結果見表4。在葉部性狀中,葉長與葉寬、葉長與葉柄長、葉寬與葉柄長均為極顯著正相關,葉寬與葉長寬比為極顯著負相關。葉長、葉寬、葉柄長與除花朵數、花梗長外的其余花部性狀以及果徑間均具有極顯著正相關關系,可見葉片性狀對花、果性狀具有促進作用。葉長寬比與花朵數為顯著正相關,與小苞片長、寬則為顯著負相關。在花部性狀中,小苞片長、小苞片寬、花萼長、萼裂片長、萼裂片寬、花徑、花瓣長、花瓣寬相互間均為極顯著正相關。花梗長與葉長、萼裂片長、花徑、花瓣長寬等為極顯著或顯著正相關;花梗長還與花朵數具有顯著正相關關系,長花梗不僅有利于較多花朵的聚生,也便于花朵為適應傳粉而進行方向調整。除花朵數、花梗長以外,其余花部性狀與果徑間也為極顯著正相關。

2.3 質量性狀多樣性

野牡丹屬植物質量性狀的多樣性主要存在于種間。以株型來說,除地菍和細葉野牡丹為低矮亞灌木外,其余種均為株高≥50 cm的直立灌木。從被毛情況來看,地菍的莖常無毛,細葉野牡丹、多花野牡丹、展毛野牡丹、野牡丹的莖上多被糙伏毛或短柔毛,毛菍的莖上則具有開展的長粗毛。從葉片特征來看,地菍葉片的主脈兩側各有1~2條側脈,葉基廣楔形,頂端急尖或廣急尖;野牡丹葉片的主脈兩側各有2~3條側脈,葉基圓形至心形,頂端急尖;其余種的葉片主脈兩側各有2條側脈,葉基多近圓形,葉尖漸尖或長漸尖。種內居群間在株型、莖被毛、基出脈、葉尖和葉基形態上的差異不大。葉柄顏色和花色性狀在野牡丹屬植物種間和部分種內居群間均有一定差異,如地菍的葉柄顏色為紅褐色,花色為粉紅色;細葉野牡丹的葉柄顏色有深淺之分,花色淡紫或粉紅;展毛野牡丹和多花野牡丹的葉柄為綠色、淺褐色或深褐色,前者花色多為紫紅色,后者花色為粉紅色或淡紫色;野牡丹的葉柄多為淺褐色,花色淡紫至紫紅;毛菍的葉柄顏色也為淺褐色,花色淡紫。種內不同居群間存在的差異表明野牡丹屬植物的葉柄顏色和花色性狀可能受居群環境因子的影響較大。

2.4 表型聚類分析

采用類平均法得到野牡丹屬不同居群的表型聚類結果(圖1),當聚類之間的平均距離為34時,可將16個居群大致分為3類,第1類為地菍,由P9、P10這2個居群聚為一類,為貼地匍匐生長的低矮亞灌木,葉、花、果明顯較小,與其它種的親緣關系可能較遠;第2類為毛菍,由P16居群自成一類,植株高大,被開展的長粗毛,葉、花、果均較大,與其它種有明顯區別;第3類包含了展毛野牡丹、多花野牡丹、細葉野牡丹和野牡丹,其植株、葉、花、果的大小介于第1、2類之間。在第3類中,展毛野牡丹與多花野牡丹最先聚在一起,其次是細葉野牡丹,野牡丹最后聚入,表明展毛野牡丹和多花野牡丹的表型相似度最高,二者與野牡丹的表型差異較大。如圖1所示,同種下的各居群均能較好地聚為一類。在展毛野牡丹種內,地理位置較近的居群較先聚類,地理位置最遠的居群表型變異最明顯,其表型分化可能在一定程度上與生境斷裂有關;而多花野牡丹、細葉野牡丹、野牡丹種內各居群并沒有嚴格依地理距離聚類,其最先聚在一起的并不是地理位置較近的居群,這可能是因為其居群間相對孤立,表型性狀的變異不連續。

2.5 表型的因子分析

表型性狀往往包含多個變量,采用因子分析對其進行降維,可以簡化分析過程,更好地描述變異構成特征。對16個居群的22個表型性狀進行因子分析,前4個因子的特征值、方差貢獻率、累積方差貢獻率如表5。由表5可知,這4個因子對方差的累積貢獻率達到了89.339%,已足夠描述野牡丹屬植物的表型性狀;其中因子1的方差貢獻率高達50.856%,可以將其看作是影響野牡丹屬植物表型分化的主要因子。對輸出的方差最大正交旋轉載荷矩陣(表6)進行分析得出,因子1在花萼長、萼裂片寬、花瓣寬、花徑、葉長、花瓣長、萼裂片長、小苞片寬、小苞片長上具有較大載荷,除葉長外,其余均為描述花部性狀的指標,可見花部形態的變異是引起野牡丹屬植物表型分化的主要因素,可將因子1看作花部形態因子;因子1在葉長、葉柄長、果徑、葉片寬、株型上也具有較大載荷,表明這些性狀與花部形態性狀緊密相關。因子1已基本涵蓋了野牡丹屬植物的主要形態特征,可作為該屬植物分類的重要依據。因子2中葉尖形態和葉長寬比是主要的分化指標,突出了葉形,也有利于分類識別,可將因子2定義為葉形因子;基出脈在因子2中具有較大的負載荷,其與葉形具有一定的負相關,一般來說,葉長寬比較大、葉片細長者基出脈數偏少,反之則多。因子3在花色、花朵數上具有較大載荷,可看作花部質量因子;葉基形態在因子3中同樣具有較大載荷,其與花色、花朵數之間可能存在一定相關性。因子4可定義為葉柄顏色因子。葉基形態和葉柄顏色均為葉部質量性狀,是對葉片特征的補充描述。

3 討論與結論

表型多樣性是遺傳多樣性與環境多樣性的綜合體現[4]。盡管同工酶標記和分子標記已廣泛應用于植物遺傳多樣性的研究,但表型性狀的研究仍然是遺傳多樣性研究最基本的方法和途徑。近年來表型多樣性研究一直比較活躍,在披堿草屬[2]、紫丁香[4]、寬刺薔薇[27]、胭脂花[28]、白花樹[29]、睡蓮屬[30]、岷江百合[31]、陸地棉[32]等植物多樣性研究中行之有效,較好地揭示了群體的遺傳規律和變異,為其遺傳多樣性的客觀評價提供了科學依據。

表型性狀變異往往具有適應和進化上的意義[2],豐富的表型變異蓄積提高了群體適應多種環境的可能性,也為種的形成提供了基礎。野牡丹屬植物分布較廣,生境多樣,其種間和種內居群間均表現出豐富的表型多樣性。本研究結果顯示,野牡丹屬植物在果橫徑、花直徑、花萼筒長、花瓣長、葉長寬比上的變異系數較小,表明該屬的花、果大小和葉形性狀較穩定,種間相似性較高;但在其余表型性狀上種間差異明顯,尤以花序花朵數、小苞片長寬、葉柄長、花梗長等性狀的變異幅度較大。野牡丹屬植物種內居群間也存在明顯的表型分化,但大多集中在易受環境影響以及由此產生的個體生長差異所引起的性狀變異,如葉片大小、葉柄長、花梗長、葉柄顏色、花色等的差異。

由于野牡丹屬植物分布較廣,形態特征變異較大,長期以來學術界對該屬植物的分類并未形成統一意見。根據Meyer 2001年所做的修訂,英文版《Flora of China》(2007年)認可該屬植物有22種,在中國境內分布的僅為5種,其中將多花野牡丹、展毛野牡丹和野牡丹(M. candidum)合并為同1新種野牡丹(M. malabathricum)[33],這一點目前在國內仍存在爭議。楊利平等[18]對國產5種野牡丹屬植物種子特征的初步研究表明,多花野牡丹、展毛野牡丹、野牡丹在種子千粒重、發芽率、種子表面微形態等方面存在較大差異。劉雪凝等[19]對5種野牡丹屬植物花粉形態的研究結果表明,野牡丹、展毛野牡丹、多花野牡丹的花粉極面觀、孔溝延伸程度和花粉外壁紋飾差異顯著,推斷野牡丹是最為原始的類型,展毛野牡丹和多花野牡丹均較其進化。陳進燎等[20]對野牡丹屬6種植物葉片表面形態特征進行掃描電鏡觀察,結果顯示多花野牡丹、野牡丹存在較大差異,經聚類分析發現多花野牡丹、展毛野牡丹與細葉野牡丹的親緣關系較近,與野牡丹的親緣關系較遠,結合花期、花型、果型的比較差異,推斷多花野牡丹和野牡丹應該為2種植物。本研究通過對野牡丹屬植物表型性狀的聚類分析結果表明,展毛野牡丹和多花野牡丹的表型相似度最高,二者與野牡丹的表型相似度較之細葉野牡丹要低,據此推測展毛野牡丹和多花野牡丹可能具有最近的親緣關系,二者與細葉野牡丹的親緣關系較近,與野牡丹的親緣關系則可能較遠,這與陳進燎等的研究結果基本一致。由此筆者認為將多花野牡丹與展毛野牡丹合并具有一定的合理性,但其與野牡丹是否能歸入同1個種值得商榷。表型性狀兼具穩定性和變異性,既是遺傳變異的表征,也易受環境的影響,因此本研究中基于表型相似度的野牡丹屬部分種間親緣關系的推斷可能不完全是其遺傳關系的真實反映,還需結合細胞學、孢粉學、分子生物學等方法進一步論證。

本研究中因子分析的結果顯示,花部形態的變異是野牡丹屬植物表型分化的主要因素,其次是葉、果和株型的變異,部分性狀如基出脈、花梗長、莖被毛的變異對野牡丹屬植物表型分化的作用較小,這表明用于野牡丹屬植物分類鑒定的表型性狀并非越多越好,應著重注意葉、花、果的形態、株型、花序花朵數、花色等較重要的性狀,避免過多次要性狀的干擾。

參考文獻

[1] 葛 頌, 洪德元. 遺傳多樣性及其檢測方法[M]. //中國科學院生物多樣性委員會. 生物多樣性研究的原理與方法. 北京: 中國科學技術出版社出版, 1994: 123-140.

[2] 祁 娟, 曹文俠, 閆偉紅. 披堿草屬野生居群表型多樣性及其與環境關系研究[J]. 西北植物學報, 2013, 33(5): 1 027-1 033.

[3] 顧萬春. 統計遺傳學[M]. 北京: 科學出版社, 2004: 165.

[4] 明 軍, 顧萬春. 紫丁香表型多樣性研究[J]. 林業科學研究, 2006, 19(2): 199-204.

[5] 陳 介. 中國野牡丹科野牡丹屬植物的研究[J]. 華南農學院學報, 1983, 4(1): 31-36.

[6] 陳 介. 野牡丹科[M]. //中國科學院中國植物志編輯委員會. 中國植物志(第53卷). 北京: 科學出版社, 1984: 152-161.

[7] 唐鐵鑫, 吳 鴻, 李啟垓. 野牡丹屬藥用植物黃酮類成分比較[J]. 中藥材, 2007, 30(8): 912-913.

[8] 姚亮亮, 劉新亞. 野牡丹屬植物的化學成分、 藥理活性及臨床應用研究概況[J]. 江西中醫學院學報, 2010, 22(6): 52-55.

[9] 劉 惠, 沈毅華, 劉 文. 野牡丹提取物對血小板聚集的影響[J]. 廣東醫學院學報, 2012, 30(5): 482-483, 487.

[10] 梁春玲, 周玖瑤, 吳俊標, 等. 野牡丹抗炎鎮痛作用及其對小鼠免疫功能影響的研究[J]. 中國藥師, 2012, 15(11): 1 547-1 550

[11] 彭東輝, 張啟翔, 黃俊婷. 中國野牡丹科觀賞植物種質資源及其在福建省的分布初步調查[J]. 中國園林, 2007, 11: 83-88.

[12] 李愛華. 地菍果紅色素的初步研究[J]. 食品工業, 1995, 3: 20-21.

[13] 朱 純. 廣東野牡丹屬植物的分布及在城市生態系統中的作用[J]. 中國野生植物資源, 2008, 27(2): 29-32.

[14] 趙 健, 仇 碩, 李秀娟, 等. 廣西野牡丹科觀賞植物種質資源及開發前景[J]. 河北農業科學, 2008, 12(7): 80-82.

[15] 田建平, 胡遠艷. 海南野生野牡丹科藥用觀賞植物資源的調查[J]. 海南醫學院學報, 2008, 14(1): 24-26, 31.

[16] 馬國華, 張靜峰, 劉 念, 等. 從多花野牡丹和野牡丹花柄直接誘導出芽[J]. 植物生理學通訊, 2004, 40(6): 719.

[17] 楊利平, 劉桂芳, 劉雪凝. 細葉野牡丹的組培快繁[J]. 東北林業大學學報, 2012, 40(9): 25-27.

[18] 林秀香. 細葉野牡丹扦插繁殖技術研究[J]. 西南林學院學報, 2009, 29(1): 42-44.

[19] 彭東輝, 邵偉麗, 陳龍菊, 等. 紫毛野牡丹的無菌播種與試管繁殖[J]. 植物生理學通訊, 2009, 45(12): 1 199-1 200.

[20] 楊利平, 屈 平, 李 燕, 等. 多花野牡丹傳粉生物學特性研究[J]. 河北農業大學學報, 2007, 30(6): 42-45.

[21] 彭東輝, 蘭思仁, 吳沙沙. 中國特有種枝毛野牡丹傳粉生物學及繁育系統研究[J]. 林業科學研究, 2014, 27(1): 11-16.

[22] 劉 惠, 趙 平, 蔡錫安, 等. 坡地和沼澤地野牡丹葉片的解剖特征與氣孔氣體交換特性[J]. 熱帶亞熱帶植物學報, 2005, 13(2): 109-116.

[23] Watanabe T & Osaki M. Influence of aluminum and phosphorus on growth and xylem sap composition in Melastoma malabathricum L.[J]. Plant and Soil, 2001, 237(1): 63-70.

[24] 楊利平, 莊 斌, 蘇正華, 等. 野牡丹屬植物種子特征的初步研究[J]. 植物遺傳資源學報, 2008, 9(2): 248-252.

[25] 劉雪凝, 屈 平. 5種野牡丹屬植物花粉形態研究[J]. 河北農業大學學報, 2012, 35(5): 63-66.

[26] 陳進燎, 蘭思仁, 吳沙沙, 等. 6種野牡丹屬植物葉片表面特征及其分類學意義研究[J]. 福建林學院學報, 2013, 33(2): 106-112.

[27] 楊樹華, 郭 寧, 葛維亞, 等. 新疆東天山地區寬刺薔薇居群表型多樣性分析[J]. 植物遺傳資源學報, 2013, 14(3): 455-461.

[28] 穆繪莉, 游曉會, 石 超, 等. 胭脂花野生群體表型多樣性研究[J]. 北京林業大學學報, 2012, 34(3): 97-102.

[29] 柳新紅, 李因剛, 趙 勛, 等. 白花樹天然群體表型多樣性研究[J]. 林業科學研究, 2011, 24(6): 694-700.

[30] 張海平, 房偉民, 陳發棣, 等. 部分睡蓮屬植物形態性狀的多樣性分析[J]. 南京農業大學學報, 2009, 32(4): 47-52.

[31] 張彩霞, 明 軍, 劉 春, 等. 岷江百合天然群體的表型多樣性[J]. 園藝學報, 2008, 35(8): 1 183-1 188.

[32] 陳 光, 杜雄明. 我國陸地棉基礎種質表型性狀的遺傳多樣性分析[J]. 西北植物學報, 2006, 26(8): 1 649-1 656.

[33] Chen J, Renner S S. Melastomataceae[M]//: Flora of China. Beijing: Science Press, 13: 363-365.

主站蜘蛛池模板: 亚洲人人视频| 免费高清毛片| 2020亚洲精品无码| 久久精品最新免费国产成人| 国产美女91视频| 99国产在线视频| 午夜影院a级片| 国产精品蜜芽在线观看| 欧美成人一级| 国产精品欧美日本韩免费一区二区三区不卡 | 亚洲天堂成人在线观看| 欧洲极品无码一区二区三区| 亚洲va欧美va国产综合下载| 欧美日韩国产一级| 国产真实乱子伦视频播放| 国产福利在线免费| 欧美A级V片在线观看| 亚洲精品无码在线播放网站| 国产综合在线观看视频| 欧美视频二区| 亚洲熟妇AV日韩熟妇在线| 久久久久久久久亚洲精品| 97成人在线观看| 亚洲天天更新| 久久综合丝袜长腿丝袜| 91毛片网| 久久中文电影| 试看120秒男女啪啪免费| 亚洲精品手机在线| 狠狠亚洲婷婷综合色香| av一区二区无码在线| a级毛片一区二区免费视频| 1级黄色毛片| 色噜噜在线观看| 2021国产精品自产拍在线观看| 国产在线麻豆波多野结衣| 精品视频一区在线观看| 亚洲综合在线网| 国产日韩欧美视频| 久久五月天国产自| 好久久免费视频高清| 国产人成在线观看| 国产在线自揄拍揄视频网站| 国产97视频在线观看| 99在线观看国产| 亚洲第一黄片大全| 国产一级视频久久| 午夜精品一区二区蜜桃| 91久久偷偷做嫩草影院免费看| 国精品91人妻无码一区二区三区| 亚洲综合色吧| 欧美日韩导航| 波多野结衣一区二区三区四区视频 | 日本久久网站| 91高清在线视频| 亚洲成a人在线观看| 日韩 欧美 国产 精品 综合| 欧美国产日韩在线| 欧美啪啪精品| 欧美精品1区2区| 999精品视频在线| 国产男人天堂| 少妇精品久久久一区二区三区| 青草娱乐极品免费视频| 久久综合伊人77777| 国产精品久久久免费视频| 国产乱论视频| 日韩高清欧美| 一级爆乳无码av| 中文字幕在线一区二区在线| 午夜啪啪福利| 国产资源免费观看| 国产激情在线视频| 欧美三级自拍| 色婷婷亚洲十月十月色天| a亚洲视频| 99在线视频免费| 国产一在线观看| 国产99视频精品免费观看9e| 欧美激情第一区| 亚洲无码精品在线播放 | 欧美区一区|