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.

主站蜘蛛池模板: 四虎成人精品| 国产不卡一级毛片视频| 九九热精品在线视频| 成人综合在线观看| 欧美日韩亚洲国产主播第一区| 国产原创第一页在线观看| 伊人久久大香线蕉综合影视| 国产成人AV男人的天堂| 丝袜美女被出水视频一区| 欧美国产成人在线| 青青草久久伊人| www.99在线观看| 2021国产精品自产拍在线观看 | 在线观看国产精品第一区免费| 久青草免费视频| 国产精品成人一区二区不卡 | 国产欧美日韩精品综合在线| 国模在线视频一区二区三区| 一本久道久综合久久鬼色| 中文国产成人久久精品小说| 欧美亚洲国产视频| 蜜桃臀无码内射一区二区三区| 国产在线视频导航| 日韩精品高清自在线| 久久精品丝袜| 97在线碰| 国产精品视频公开费视频| 成人在线不卡| 国产一区免费在线观看| 亚洲无线一二三四区男男| 国产成人一二三| 国产成人区在线观看视频| 71pao成人国产永久免费视频 | 国产91全国探花系列在线播放| 又黄又湿又爽的视频| AV色爱天堂网| 国产日韩欧美一区二区三区在线| 国产成人精品18| 亚洲男人的天堂久久香蕉网| 国内精品自在自线视频香蕉| 亚洲国产精品人久久电影| 日韩中文无码av超清| 91色在线视频| 欧美日韩v| 欧美天天干| 午夜福利在线观看成人| 综合人妻久久一区二区精品| 欧美日本中文| 亚洲视频免| 久久亚洲综合伊人| 午夜啪啪福利| 久久天天躁夜夜躁狠狠| 国产97视频在线观看| 精品国产网站| 三级国产在线观看| 国产99视频免费精品是看6| 免费精品一区二区h| 亚洲欧美精品一中文字幕| 色综合久久无码网| 亚洲黄网视频| 看看一级毛片| 91在线激情在线观看| 国产国模一区二区三区四区| 在线免费看片a| 亚洲无码熟妇人妻AV在线| 91精品国产丝袜| 色网站在线视频| 免费又黄又爽又猛大片午夜| 国产又粗又猛又爽视频| 天天色天天综合网| 精久久久久无码区中文字幕| swag国产精品| 久久性妇女精品免费| 精品無碼一區在線觀看 | 黄色网址免费在线| 91精品啪在线观看国产| 久久人人妻人人爽人人卡片av| 午夜限制老子影院888| 久久亚洲美女精品国产精品| 国产福利在线免费观看| 国产精品网址你懂的| 久久不卡国产精品无码|