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

40份辣椒種質資源表型性狀的遺傳多樣性分析

2025-08-22 00:00:00王童童郭志偉任廣乾孫強劉賀娟張蓓程相杰靳正雅任福森甄俊琦吳濤李儀
中國瓜菜 2025年8期

中圖分類號:S641.3 文獻標志碼:A 文章編號:1673-2871(2025)08-067-10

Genetic diversity analysis of phenotypic traits in 40 Capsicum annuum L. germplasmresources

WANG Tongtong', GUO Zhiwei',REN Guangqian', SUN Qiang',LIU Hejuan', ZHANG Bei', CHENG Xiangjie', JIN Zhengya2,RENFusen', ZHEN Junqi',WU Tao1,LI Yi3

(1.Xinxiangdulcinnan;ianicuualolo ServiceCeegi;eUtythto

Abstract:To enhance theutilizationof thecollected Capsicum annuumL.germplasmresources and improve breeding efficiency,geneticdiversityanalysis,clusteranalsis,correlationanalysis,andprincipalcomponentanalysiswrecaed out n this study with23 phenotypic traits of40 stable inbredlines.Theresultsshowed thatthe diversity indices of the 23 phenotypictraitsvariedfromOto1.997,ofwhichthehighestdiversityindex wasplant width,andthecoeficientofariation of single fruit mass was the largest among the quantitative traits ( 52.92% ,which indicated that thephenotypic traits of the 40C? annumL.germplasmresources had rich genetic diversity.Theresults ofcluster analysis showed that the 40 C.anmum L.germplasmresources could cluster into two taxa at geneticdistance 20,andthe germplasmresources of taxonIwere mostlycompactplants,with elongatedandthin-skinned fruitsandrelativelylarge fruit shapeindex,which could provide the theoretical basis forthe breeding practiceofthe subsequent diferentfruit shapesand plant shapes.The resultsofcorrelationanalysis showed thatthecorelationcoefficientsof74pairs oftraits weresignificantorextremely significant,indicating that mostofthephenotypictraits interacted witheach other.Theresults ofthe principalcompoent analysis showed that the cumulative contribution rate of the first six principal components was as high as 77.493% ,which couldreflectthe basiccharacteristicsofthephenotypic traits,the twelvephenotypic traitsas fruit transverse diameter, fruit shapeindex,andfruitcro-sectionshapecouldbeusedasimportantindexesforevaluatingthegeneticdiversityof C. annuum L.'s phenotypes and genetic improvement of C annuum L. variety. This study provides an important theoretical reference basis for the utilization and innovation of C annuumL.germplasm resources and variety selection. Key words: Capsicum annuumL.; Germplasm resource; Phenotypic trait; Genetic diversity

辣椒(CapsicumannuumL.)原產于中南美洲熱帶及亞熱帶地區,現已遍布世界多個國家,大約在16世紀后期傳入中國,在中國經過了400多年的發展,已成為中國餐桌上重要的蔬菜之一。辣椒對光周期不敏感,種植范圍廣;種類繁多,從植物學上可分為毛辣椒、長柄辣椒、木本辣椒、一年生辣椒等,其中一年生辣椒又可分為長椒、燈籠椒、圓錐椒、簇生椒;加工方式多樣,可以做干制椒、辣椒粉、油辣椒、剁辣椒、辣椒醬、泡辣椒等;消費功能多,除菜用、調味、加工功能外,還具有藥用、飼用、觀賞等功能。作為中國“菜籃子\"的主要蔬菜作物之一,其產業大、效益好、發展廣。據統計,21世紀以來,辣椒已成為中國種植面積最大的蔬菜[1-3]。

辣椒種質資源的收集與利用能夠為辣椒的品種選育提供大量研究素材,為后續辣椒育種工作打下堅實基礎,遺傳多樣性是物種在演化的過程中不斷積累的結果,對所收集到的種質資源遺傳多樣性進行分析能夠很好地展現出群體間的遺傳差異,是種質資源研究的主要內容之一。隨著生物技術的不斷發展,分析遺傳多樣性的方法越來越多,目前有形態標記、生理標記、細胞標記、分子標記等[3]。

通過表型性狀指標對植物遺傳多樣性進行研究快速、直觀,已在馬鈴薯、谷子、棉花、棗、燕麥、番茄、蘋果等多種物種中廣泛應用[410]。截至目前,關于辣椒種質資源表型性狀遺傳多樣性的研究已有較多報道,劉林婭等1對81份辣椒種質資源的20個表型性狀進行了遺傳多樣性分析,為更好地利用這些種質資源提供了重要參考。吳迪等[12對275份辣椒的18個表型性狀進行了遺傳多樣性分析,表明這些材料具有豐富的多樣性,具有較好的遺傳改良基礎。朱珍花等[3對32份干制辣椒進行了遺傳多樣性分析,為干制辣椒種質資源創新和育種提供了良好的理論依據。筆者以收集到的40份穩定自交系種質作為原始種質資源庫,調查21個表型性狀并進行遺傳多樣性分析,以期為加強辣椒種質資源利用、提高育種效率提供理論依據。

1 材料與方法

1.1材料

試驗材料為蔬菜研究所收集的40份辣椒種質,各種質來源于多個地區,其中河南新鄉5份,天津2份,河南鄭州6份,甘肅19份,河南洛陽、山東、內蒙古赤峰各1份,美國3份,云南2份(表1。各種質材料分別于2022年12月在研發基地的溫室中使用50孔穴盤進行育苗,2023年3一7月以完全隨機排列的方式種植在研發基地的溫室大棚中,采用地膜覆蓋加膜下滴灌栽培模式,病蟲害防治、肥水等田間管理基本一致,株行距為 40cm× 60cm ,小區面積 5.4m2 ,每小區22株,3次重復。

表1種質資源編號及來源

Table1Numberandoriginof germplasmresources

1.2 測定項目和方法

依據《辣椒種質資源描述規范和數據標準》14提供的辣椒表型性狀描述及賦值標準對21項表型性狀進行觀察和測量。各性狀賦值標準見表1,其中株高、株幅、果縱徑、果橫徑、肉厚、心室數、單果質量、葉柄長、葉片長、葉片寬等10個數量性狀直接采用鋼尺、游標卡尺以及電子秤實測數值,果形指數 Ψ:=Ψ 果縱徑/果橫徑,葉形指數 葉片長/葉片寬。每份種質材料隨機取樣10株,果實取發育正常的對椒,葉片于四門斗始花期取植株中部完整且生長正常的最大葉片。

1.3 數據分析

采用Excel2010處理分析數據,得出數量性狀的基本統計量(最小值、最大值、平均值、標準差和變異系數)和質量性狀的分布頻率。遺傳多樣性指數H(Shannon一wiener index)計算公式: 1nPj ,式中 Pj 為第 j 等級中種質數占總種質數的比例。數量性狀參照胡建斌等5的方法,根據數據的平均值( )和標準差 (SD) 將數據分為10個等級,第 i 性狀的第1級 ,第10級 2SD),每一級相差 0.5SD ,每一級的相對頻率用于計算多樣性指數。采用SPSS22.0進行聚類分析、相關性分析和主成分分析。

2 結果與分析

2.1辣椒種質資源質量性狀多樣性分析

統計9個質量性狀的頻率分布并計算各個性狀多樣性指數,結果如表2所示。遺傳多樣性指數變化范圍為0.464~1.180,其中果實橫切面形狀的多樣性指數最高,其次是分枝性(1.067),其他質量性狀的多樣性指數均小于1。花藥顏色中藍和紫占比相近,分別為 52.5% 和 45.0% ,占比最小的淺藍色僅

表2質量性狀頻率分布和多樣性指數

Table2 Distribution frequency and Shannon-wiener index ofquality trait

注:括號中數據表示分布頻率 (%) 。Note:Data in theparentheses indicate the frequency of distribution (%)

2.5% ;株型有半直立 (82.5% 和直立( 17.5% 兩種;葉形有長卵圓形 (77.5% 和披針形 (22.5% 兩種;葉面特征為微皺( 82.5% 或皺 (17.5% );分枝性以弱為主L 45.0% );花柱長度多為雄蕊近等長,占 65.0% ;花梗著生狀態以直立為主,占 67.5% ,而下垂和側生均在15.0% 左右;莖大多表現為無茸毛,占 70.0% ;果實橫切面形狀主要為近四邊形 47.5% )。

2.2辣椒種質資源數量性狀多樣性分析

對辣椒種質資源的12個數量性狀進行田間調查和統計,由表3可知,各數量性狀的多樣性指數均 gt;1 ,說明這些性狀的指數是合格的,變幅為1.404~1.997,說明這些性狀的表型多樣性較為豐富,其中株幅的多樣性指數最高,其次為果橫徑(1.973)和肉厚(1.948)。變異系數最大的是單果質量 (52.92% ),變異范圍為 9.53~154.11g ,其次是果形指數 (47.77% )、肉厚 (37.74% )、果橫徑 33.63% )和果縱徑 (30.12% ),其他數量性狀的變異系數介于17.01%~26.70% ,其中株幅的變異系數最小,說明供試材料在此性狀上表現出穩定的遺傳。

2.3辣椒種質資源表型性狀聚類分析

使用SPSS22.0軟件對辣椒種質資源的21個表型性狀進行聚類分析,如圖1所示,在遺傳距離20處可以將辣椒種質分為I和II2個類群,在遺傳距離15處可以將類群ⅡI再分為2個子類群ⅡI-1和II-2,各類群的主要表型性狀如表4所示。

類群I包含8份種質,類群II包含32份種質。2個類群相比,類群I的果橫徑與葉柄長極顯著大于類群ⅡI,且類群I的株幅、肉厚、葉片長、葉片寬均顯著大于類群II,而類群Ⅱ的果形指數顯著大于類群I。總體來看,類群I的主要特征是葉柄較長,葉片較大,植株較高且伸展幅度較大,果實較粗且厚,果形指數較小;類群II果實偏細長且果皮略薄,果形指數較大,葉片較窄小且葉柄較短,植株較矮且緊湊。

在類群ⅡI的2個子類群中,類群II-1包含10份種質,類群II-2包含22份種質。類群II-2株高和類群II-1相差無幾,株幅極顯著大于類群II-1,且果實橫徑、肉厚、心室數、單果質量均極顯著大于類群ⅡI-1,果縱徑較類群II-1大且差異達到顯著水平,除葉片長、葉片寬、果形指數、葉形指數小于類群II-1外,其余性狀均大于類群II-1,但差異不顯著。總體來看,類群ⅡI-1的主要特征為植株較為緊湊,果實偏小,心室數偏少但果形指數較大,葉片較大且葉形指數較大但葉柄長稍短;類群II-2植株伸展幅度較大,果實整體較大且肉厚,但果形指數小,葉片略小但葉柄長。

2.4辣椒種質表型性狀相關性分析

對辣椒種質資源的表型性狀進行相關性分析,結果(表5)顯示,各性狀之間相互聯系,辣椒種質資源各表型性狀間有74對性狀存在顯著或極顯著相關性。其中株型與株高呈極顯著正相關(0.542);分枝性與株高和果形指數呈極顯著正相關,與果實橫切面形狀、果橫徑、肉厚、心室數、單果質量呈極顯著負相關;葉面特征與花梗著生狀態、葉柄長、葉片長、葉片寬呈極顯著正相關,相關系數分別為0.552、0.561、0.560、0.596;花梗著生狀態與果實橫切面形狀、葉片長、葉片寬呈極顯著正相關,與果縱徑呈極顯著負相關;果實橫切面形狀與果橫徑、肉厚、心室數、單果質量呈極顯著正相關,相關系數分別為0.645、0.473、0.768、0.450,與果形指數呈極顯著負相關(-0.545);株幅與果縱徑、單果質量呈極顯著正相關,相關系數分別為0.519、0.468;果縱徑與果形指數呈極顯著正相關;果橫徑與肉厚、心室數、單果質量、葉柄長呈極顯著正相關,相關系數分別為 0.762,0.662,0.791,0.470 ,與果形指數呈極顯著負相關,相關系數為-0.837;肉厚與心室數、單果質量、葉柄長呈極顯著正相關,相關系數分別為0.596、0.766、0.463,與果形指數呈極顯著負相關(-0.638) ;心室數與單果質量呈極顯著正相關(0.517),與果形指數呈極顯著負相關 (-0.561) ;單果質量與葉柄長呈極顯著正相關(0.428),與果形指數呈極顯著負相關(-0.472);葉柄長與葉片長、葉片寬呈極顯著正相關,與果形指數呈顯著負相關;葉片長與葉片寬呈極顯著正相關。綜上所述,辣椒種質資源的大部分表型性狀之間相互影響,通過合理組配調控株高、葉片等性狀可增加單產。

表3數量性狀多樣性統計

Table3Diversitystatisticsof quantitative traits

圖140份辣椒種質資源表型性狀聚類分析

Fig.1Cluster analysis of phenotypic characters of 40 C.annuum L.germplasm resources

2.5辣椒表型性狀主成分分析

對40份辣椒種質資源表型性狀進行主成分分析(表6),提取特征值大于1.0的前6個主成分,累計貢獻率達 77.493% ,前6個主成分代表了辣椒表型性狀的主要信息量,說明這6個主成分可以反映辣椒種質表型性狀的基本特征。其中,第1主成分的特征值為6.425,貢獻率為 30.597% ,特征向量絕對值較高的有果橫徑、果形指數和果實橫切面形狀,主要反映果實形狀;第2主成分的特征值為3.348,貢獻率為 15.941% ,特征向量絕對值較高的有株高、株型和葉片寬,主要反映株型和葉片形狀;第3主成分的特征值和貢獻率分別為2.513和11.966% ,株幅、單果質量、花梗著生狀態的特征向量絕對值較高;第4主成分的特征值為1.646,貢獻率為 7.840% ,主要反映莖茸毛;第5、6主成分的特征值分別為1.263和1.078,貢獻率分別為 6.014% 和 5.134% ,特征向量絕對值較高的分別為花柱長度和葉形指數。從特征值和貢獻率來看,果橫徑、果形指數、果實橫切面形狀、單果質量、株高、株型、株幅、葉片寬、葉形指數、莖茸毛、花梗著生狀態和花柱長度共12個表型性狀是引起辣椒種質表型差異的主要因素,可以為辣椒種質資源的利用、創新與評價提供一定的參考依據。

表4種質各類群表型性狀

Table4 Phenotypic traitsofdifferent germplasm groups

注:*和**分別表示類群I和類群II在0.05和0.01水平差異顯著; Δ?ΔΔ 分別表示類群II-1和類群II-2在0.05和0.01水平差異顯著。 Note: * and ** indicate significant difference between group Iand group II at 0.05 and 0.01 level,respectively; ▲ and indicate significant dif'erencebetween group II-1 and group II-2 at 0.O5 and 0.01 level, respectively.

3 討論與結論

種質資源是作物遺傳改良和相關基礎研究的物質基礎[。根據作物主要農藝性狀中的數量性狀和質量性狀進行遺傳多樣性分析,能夠為后續品種改良提供良好的理論基礎。因此筆者對從不同地區引進的40份辣椒種質資源的21個表型性狀進行多樣性分析,發現辣椒種質資源表型性狀類型豐富,質量性狀中多樣性指數較高的為果實橫切面形狀(1.180)和分枝性(1.067);數量性狀中多樣性指數最高的是株幅(1.997),果橫徑(1.973)和肉厚(1.948)次之,變異系數最大的是單果質量 (52.92% ),其次是果形指數 47.77% )、肉厚 (37.74% )、果橫徑 (33.63% )和果縱徑 (30.12% ,與質量性狀相比,數量性狀多樣性指數更高,變異更大,遺傳資源更豐富,這與李寧等[7]、趙紅等[18]的研究結果一致。

根據聚類分析結果可以將40份辣椒種質資源劃分為2個類群,類群I的主要特征是葉柄較長,葉片較大,植株較高且伸展幅度較大,果實較粗且肉厚,果形指數較小;類群II包含32份材料,占供試種質材料總數的 80% ,該類群果實偏細長且果皮略薄,果形指數較大,葉片較窄小且葉柄較短,植株較矮且緊湊。郝衛等研究表明,葉形、株型、果形是最關鍵的性狀,本試驗結果與其一致。此外,聚類結果顯示,不同地區的辣椒種質可以聚為一類,這與劉林婭等的研究結果一致,說明不同原始生境并不是造成辣椒種質表型差異的唯一因素,還與種質遺傳分化和基因交流有關。

針對表型性狀間的相關性分析能夠揭示各性狀之間的關聯性,通過重要性狀的表觀選擇和遺傳改良間接地同步改良次要性狀,從而加速育種進程[20]。相關性分析結果表明,各表型性狀之間存在錯綜復雜的關系,其中果形指數與果實橫切面形狀、果橫徑、肉厚、心室數、單果質量等果實性狀及葉柄長呈極顯著負相關,與分枝性和果縱徑呈極顯著正相關,與葉片長、葉片寬呈顯著負相關,因此若想改變辣椒果形,要重視對植株分枝性和葉形的觀察與利用,這與劉林婭等的研究結果一致。育種選擇時可利用農藝性狀間的相關性對一些目的性狀進行選擇,果縱徑、果橫徑、肉厚、單果質量這些影響產量的表型性狀均與分枝性顯著相關,且除果橫徑外,果縱徑、肉厚、單果質量均與株幅呈顯著或極顯著正相關,因此若想提高產量,可多關注分枝性和株幅。

通過相關分析了解性狀之間的相關性,再通過主成分分析,能夠利用降維的方式將多個指標轉化為較少的幾個主成分,從而更加清楚地顯示辣椒各性狀的重要程度[2。筆者提取的特征值大于1.0的前6個主成分累計貢獻率高達 77.493% ,從特征值和貢獻率來看,表型性狀中果橫徑、果形指數、果實橫切面形狀、單果質量、株高、株型、株幅、葉片寬、葉形指數、莖茸毛、花梗著生狀態和花柱長度這12個表型性狀能夠作為評價辣椒表型遺傳多樣性和棘椒品種遺傳改良的重要指標。

綜上所述,筆者通過對40份辣椒種質資源的21個表型性狀進行多樣性分析,揭示了辣椒種質資源的豐富遺傳多樣性。結果表明,數量性狀的多樣性指數和變異系數均大于質量性狀,尤其是株幅、果橫徑和肉厚等性狀表現出較高的遺傳變異潛力。聚類分析將種質資源分為兩個主要類群,進一步揭示了果實形態、植株結構和葉片特征等關鍵性狀的遺傳分化規律。相關性分析和主成分分析表明,果形指數、果實橫切面形狀等性狀與產量性狀,單果質量密切相關,研究結果為辣椒品種的遺傳改良提供了重要依據。未來研究可進一步探討這些性狀的遺傳機制,并結合分子標記輔助選擇,加速辣椒育種進程。

參考文獻

[1] 鄒學校,馬艷青,戴雄澤,等.辣椒在中國的傳播與產業發展[J].園藝學報,2020,47(9):1715-1726.

[2] 鄧惠如,張素勤,耿廣東.辣椒逆境響應基因研究現狀及前景[J/OL].分子植物育種,1-12(2023-11-09)[2025-03-03].https://link.cnki.net/urlid/46.1068.S.20231108.1149.002.

[3] 董文堂,鄧昌蓉,侯全剛,等.辣椒種質資源遺傳多樣性分析及核心種質構建研究進展[J].青海農林科技,2023(3):46-49.

[4] 徐曉,楊夢穎,滿全財,等.195份馬鈴薯種質資源表型性狀綜合評價[J].核農學報,2023,37(9):1710-1722.

[5] 解慧芳,牛靜,邢璐,等.117份谷子核心種質資源表型性狀的遺傳多樣性分析[J].江蘇農業科學,2023,51(13):76-81.

[6] 王業舉,張虎,張博,等.235份陸地棉表型性狀遺傳多樣性分析[J].江蘇農業學報,2023,39(3):636-644.

[7] 賈平平,靳娟,阿布都卡尤木·阿依麥提,等.新疆118個棗品種表型性狀比較分析[J].西北農業學報,2023,32(6):887-898.

[8] 耿小麗,張榕,張少平,等.130份燕麥種質表型性狀多樣性分析及評價[J].草業科學,2020,37(10):2022-2034.

[9] 趙云霞,顏秀娟,王學梅,等.246份番茄種質資源表型性狀的遺傳多樣性[J].江蘇農業科學,2021,49(17):134-140

[10]張學超,任海龍,唐式敏,等.伊犁天山160份野蘋果種質資源表型性狀的遺傳多樣性分析[J].植物遺傳資源學報,2021,22(6):1521-1530.

[11]劉林婭,黃亞成,楊那,等.81份辣椒種質資源表型性狀的遺傳多樣性分析[J].熱帶作物學報,2023,44(4):706-715.

[12]吳迪,付文婷,吳康云,等.275份辣椒種質資源表型性狀的遺傳多樣性分析[J].中國瓜菜,2024,37(9):47-53.

[13]朱珍花,濮丹,楊洋,等.干制辣椒種質資源表型性狀遺傳多樣性分析[J].北方園藝,2024(6):1-8.

[14]李錫香,張寶璽.辣椒種質資源描述規范和數據標準[M].北京:中國農業出版社,2006.

[15]胡建斌,馬雙武,王吉明,等.基于表型性狀的甜瓜核心種質構建[J].果樹學報,2013,30(3):404-411.

[16]黎裕,李英慧,楊慶文,等.基于基因組學的作物種質資源研究:現狀與展望[J].中國農業科學,2015,48(17):3333-3353.

[17]李寧,王飛,姚明華,等.國內外辣椒種質資源表型性狀多樣性及相關性分析[J].辣椒雜志,2015,13(1):8-13.

[18]趙紅,曹亞從,張正海,等.我國辣椒核心種質資源園藝性狀多樣性的分析和評價[J].中國蔬菜,2018(1):25-34.

[19]赫衛,張慧,董延龍,等.辣椒種質資源形態學性狀相關性、主成分與聚類分析[J].北方園藝,2018(4):9-17.

[20]梁永書,占小登,高志強,等.超級稻協優9308衍生群體根系與地上部重要農藝性狀的相關性[J].作物學報,2011,37(10):1711-1723.

[21]裴蕓,徐秀紅,陸錦彪,等.151份貴州地方櫻桃番茄資源的遺傳多樣性分析[J].浙江農業學報,2022,34(2):310-316.

主站蜘蛛池模板: 四虎影视国产精品| 美女无遮挡拍拍拍免费视频| 最新国产成人剧情在线播放| 午夜国产在线观看| 永久免费精品视频| 91久久偷偷做嫩草影院免费看| 2018日日摸夜夜添狠狠躁| 亚洲日韩精品综合在线一区二区| 国产精品三区四区| 三上悠亚精品二区在线观看| 日韩精品一区二区三区中文无码| 国产福利小视频在线播放观看| 老司机aⅴ在线精品导航| 亚洲成人一区在线| 一本一本大道香蕉久在线播放| 日本高清成本人视频一区| 亚洲欧美自拍中文| 国产麻豆va精品视频| 在线观看无码av免费不卡网站| 亚洲国产精品日韩欧美一区| 黄色国产在线| 精品国产自在在线在线观看| 男女性色大片免费网站| 色一情一乱一伦一区二区三区小说| 中文精品久久久久国产网址| 不卡网亚洲无码| 全裸无码专区| 亚洲一区第一页| 在线免费不卡视频| 91国内视频在线观看| 精品久久久久久久久久久| 99这里精品| 中文字幕在线一区二区在线| 国产毛片片精品天天看视频| 在线国产欧美| 欧美另类视频一区二区三区| 一本一道波多野结衣一区二区| 久久精品波多野结衣| 亚洲国产一成久久精品国产成人综合| 天堂在线www网亚洲| 制服丝袜亚洲| 无码在线激情片| 国模沟沟一区二区三区| 国产女人18毛片水真多1| 国产精品香蕉在线观看不卡| 91免费在线看| a毛片基地免费大全| 玖玖精品在线| 国产高潮流白浆视频| 无码av免费不卡在线观看| 91丝袜乱伦| 99视频只有精品| 在线不卡免费视频| 国产精品手机在线播放| 亚洲91在线精品| 亚洲国产成人麻豆精品| 久久semm亚洲国产| 欧美综合一区二区三区| 午夜福利亚洲精品| 91精品啪在线观看国产| 久草视频中文| 91美女视频在线| 中文字幕亚洲乱码熟女1区2区| 久久国产精品影院| yy6080理论大片一级久久| 少妇精品久久久一区二区三区| 日韩欧美网址| 久久a毛片| 国产小视频a在线观看| 在线免费不卡视频| a亚洲视频| 色综合狠狠操| 国产一级毛片yw| 亚洲精品图区| 中文字幕无线码一区| 色综合久久无码网| 婷婷亚洲视频| 91精品啪在线观看国产91| 午夜国产在线观看| 欧美日韩国产高清一区二区三区| 色婷婷成人| 国产精品免费福利久久播放|