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

154份西瓜育種材料果實性狀的主成分分析及聚類分析

2019-12-16 01:41:53于玉紅馬運糧趙小龍李超張琳耿少援
中國瓜菜 2019年11期

于玉紅 馬運糧 趙小龍 李超 張琳 耿少援

摘? ? 要: 為研究西瓜育種材料果實性狀品種間相似度和遺傳多樣性,對154份西瓜育種材料的14個果實性狀指標進行變異分析、主成分分析以及聚類分析。結果表明,14個果實性狀的Shannon多樣性指數變化范圍1.12~2.04,平均值為1.71,變異系數范圍為8.42%~62.51%,其平均值為28.38%;對西瓜果實性狀進行主成分分析,可將8個數量性狀指標歸為4個主成分,其累積貢獻率為83.78%。聚類分析結果表明,當取歐氏距離閾值為32.89時,供試西瓜材料可分為4大類群。上述研究結果表明,西瓜育種材料果實性狀具有較為豐富的多樣性表現,這對于西瓜品種鑒定、種質創新和雜交育種親本選配具有重要的應用價值。

關鍵詞: 西瓜; 育種材料; 果實性狀; 多樣性指數; 聚類分析

Abstract: In order to study the similarity and genetic diversity of fruit traits of watermelon breeding materials, the variation analysis, principal component analysis and clustering analysis of 14 fruit related traits were investigated for 154 watermelon breeding materials. The results show that the Shannon diversity index of 14 fruit traits varied from 1.12 to 2.04, with an average of 1.71, and the coefficient of variation ranged from 8.42% to 62.51%, with an average of 28.38%. Principal component analysis suggested that 12 fruit related traits were summarized into three principal components, and the cumulative contribution rate was 71.99%. By analyzing the principal component of watermelon fruit, 8 quantitative character indexes can be classified into 4 principal components, and their cumulative contribution rate is 83.78%. Based on the cluster analysis,the 154 watermelon breeding materials were divided into four groups in Euclidean distance 32.89 based on 14 biological characteristics. This study indicate that the fruit related traits in watermelon germplasms have rich diversity, it has important application value for watermelon variety identification, germplasm innovation and cross breeding parental selection.

Key words:Watermelon; Breeding materials; Fruit traits; Diversity index; Clustering analysis

隨著人們生活水平的提高,大家對西瓜產品的需求也發生了改變,在追求其高品質的基礎上,更加追求“新、奇、特”等多樣性,這就需要育種工作者掌握豐富的西瓜種質資源,以適應現代農業“優質、多樣、安全、個性”的發展方向。利用親緣關系較遠、較強配合力、利于生產的育種理論,結合國外資源材料的引進,更深入、細致地整理挖掘現有資源是西瓜育種的關鍵。西瓜品質的好壞由果實性狀的差異決定,因此,開展西瓜果實性狀的遺傳分析,有助于西瓜品種的遺傳改良和育種。目前在其他園藝作物資源如杧果、南瓜、甜瓜等的植物學性狀和分子水平都做了較多的遺傳多樣性研究[1-3],近幾年我國越來越注重西瓜資源的保存與創新工作,針對西瓜種質資源遺傳親緣關系和多樣性的研究相繼也有報道,其中馬雙武等[4]分析了我國保存的西瓜資源現狀并提出了建議;為了研究和利用從美國資源庫引進的1 373份西瓜種質資源,范敏等[5]利用聚類分析和主成分分析方法對觀測的12項農藝性狀數據結果進行了分析;為了更準確地了解植物遺傳多樣性,李艷梅等、張愛萍等[6-7]分別利用分子標記技術研究了部分西瓜資源基因組DNA多態性;尚建立、紀海波、陳萍等[8-10]對西瓜的植物學性狀遺傳多樣性進行了研究,結果表明,西瓜不同表型性狀指標的變異系數和多樣性指數較大,表明其具有豐富的遺傳多樣性。

經過多年搜集和積累,筆者所在課題組目前已有400多份西瓜種質資源。對其性狀的了解是進一步利用和創新的基礎?,F對中心可溶性固形物含量較高的154份西瓜種質資源果實性狀的變異性進行了研究,并進行了主成分分析和聚類分析,以期為這些種質資源的利用提供參考和依據。

1 材料與方法

1.1 材料

試驗材料所選用的154份西瓜種質資源(見表1)均為開封市蔬菜科學研究所通過多代選擇獲得的自交純系栽培種,并來自于不同地域,其中日本11份、韓國5份、荷蘭3份、美國1份,其余來自于國內,臺灣25份、湖南6份、北京11份、上海3份、河南75份、新疆6份、山東4份、湖北1份、河北3份。

1.2 方法

試驗分別于2014—2017年在河南開封花生莊進行,試驗地土壤肥力中等,采用高壟雙覆膜的栽培模式,肥水管理與生產相同,每年每份材料種植20株,隨機排列,株行距0.3 m×1.8 m,雙蔓整枝,每株留1果,每份材料調查5株,數據資料為3~4年的平均值,觀察記載6個質量性狀和8個質量性狀統計分析,包括:果實形狀、果皮底色、果皮覆紋形狀、果皮覆紋顏色、果肉色澤、果肉質地、果皮厚度、果實中心可溶性固形物含量、果實邊部可溶性固形物含量、果實縱徑、果實橫徑、單瓜質量、果形指數、中邊可溶性固形物含量差14個性狀,其中數量性狀為5個單瓜的平均值,將質量性狀按照《西瓜種質資源描述規范和數據標準》[11]進行分級和賦值,具體見表2,與數量性狀一起建立原始矩陣。

1.3 數據處理與分析方法

將果實縱徑、果實橫徑、果皮厚度、中心可溶性固形物含量、邊部可溶性固形物含量、單瓜質量、果形指數、中邊可溶性固形物含量差8個性狀的測量值用Excel計算出平均值、標準差、變異系數。根據平均值、標準差將各性狀測量數據分為10級,從第1級{Xi<(x-2δ)}到第10級{Xi≥(x +2δ)},其中i表示級數,x和δ分別表示調查性狀的平均值和標準差,Xi為性狀i級的取值,每0.5δ為1級,統計各性狀在10個級別中的分布頻率。每級的相對頻率(Pi)用于計算多樣性指數。利用Shannon-weaver遺傳多樣性指數來衡量性狀遺傳多樣性大小。計算公式為:H=-ΣPilnPi。其中Pi為某性狀第i級內材料份數占總份數的百分比,ln為自然對數。

用DPS統計分析軟件進行數量形狀的主成分分析,再使用DPS統計分析軟件根據歐氏距離法計算各樣品間的遺傳距離,同時作數據標準化處理按照離差平方和法建立聚類分析樹狀圖。

2 結果與分析

2.1 西瓜果實性狀的變異性和多樣性分析

由表3可知,154份西瓜材料涵蓋了西瓜的大部分類型,體現了材料的多樣性(表3)。果實形狀以圓形果實為主,有92份,占總數的59.7%,橢圓形的和圓柱形的分別有22份和24份,占14.3%和15.6%,橄欖形的有16份,占10.4%。果皮顏色以綠皮類型居多,占40.9%,其次是淺綠和墨綠皮類型,分別占23.4%和16.2%,綠白皮和深綠皮各有11份,僅有2份黃皮;從覆紋顏色來看,以墨綠色最多,占49.4%,其次是綠色,占24.0%,有8份西瓜材料沒有覆紋;從覆紋形狀來看,網條和齒條、條帶基本一樣多。從果肉顏色來看,以紅瓤類型為主,占37.7%,大紅瓤類型占20.9%,黃瓤和桃紅色類型同樣多,占7.1%,粉紅瓤類型占9.7%,還有1份乳白瓤類型;從果肉質地來看,51.3%的西瓜材料果肉質地為脆,質地為軟和酥脆的同樣多。

由表4可知,154份西瓜材料14個表型性狀中,變異系數平均值為28.38%,其中,果實形狀變異范圍最大,變異系數為62.51%;其次是中邊可溶性固形物含量差、單瓜質量、果肉質地、果皮覆紋形狀、果皮厚度、果肉顏色、果皮覆紋顏色變異范圍較大,變異系數分別為37.53%、35.44%、34.01%、32.32%、30.87%、27.15%、25.07%,這幾個性狀的變異系數均在25%以上,表明這些性狀個體間的差異較大,說明在育種改良上還有很大空間。中心可溶性固形物含量變異系數最小,變異系數為8.42%,在現有品種的基礎上被改良的空間較小,這與目前我們以高糖西瓜作為主要育種目標有很大關系。

14個性狀指標的多樣性指數平均值為1.71,邊部可溶性固形物含量多樣性指數最大,為2.04,果實形狀多樣性指數最小,為1.12??v徑、橫徑、單瓜質量、果實中心可溶性固形物含量、果皮厚度、果肉顏色、中邊可溶性固形物含量差與邊部可溶性固形物含量的多樣性指數高于平均值,表明供試西瓜材料在以上8個性狀中每組(根據平均值與標準差將各個性狀觀測值分為10組)內的分布較為均勻。果實形狀的多樣性指數最小,表明供試西瓜材料在這個性狀上表現型相對較少,大部分為圓形瓜,占比達92%。

2.2 西瓜果實數量性狀的主成分分析

由于西瓜果實性狀之間存在相關性,因此可用主成分分析的方法確定各個性狀指標對于果實綜合性狀的影響。利用主成分分析評價作物品種各綜合指標的主成分,可為合理選配親本及快速篩選評價指標提供理論依據。通過對西瓜果實8個數量性狀進行主成分分析,以特征向量值大于1為標準,可將8個性狀歸為4個主成分,代表了果實性狀的主要遺傳信息。由表5可知,第一主成分的貢獻率為32.3%,除中心可溶性固形物含量和邊部可溶性固形物含量的載荷符號為負向外,其他性狀的的載荷符號都是正向,說明在其他性狀都適宜的情況下,中心可溶性固形物含量和邊部可溶性固形物含量是影響西瓜果實整體品質的主要性狀,這和我們平時在選育新品種中過于追求高糖度西瓜新品種有關。第二主成分的貢獻率為23.0%,單瓜質量和橫徑2個指標為負向載荷符號,說明單瓜質量的降低是影響果實產量的主要性狀,而橫徑又是影響西瓜外觀形象的主要性狀特征,這與我們近幾年注重選育圓形瓜有關。第三主成分的貢獻率為14.9%,第四主成分的貢獻率為13.6%。4個主成分特征值累計貢獻率達到83.8%。因此可用4個主成分代替原來8個指標對西瓜果實性狀進行評價。以上結果表明,果實邊部可溶性固形物含量、單瓜質量、縱徑、果形指數、中心可溶性固形物含量、中邊可溶性固形物含量差、橫徑是影響所選西瓜品種果實性狀的主要因子,可作為劃分品種類群的依據。

2.3 西瓜果實性狀的聚類分析

利用聚類分析對品種資源進行分析,可以真實反映品種的綜合性狀,可為育種取材提供客觀依據。由數量遺傳學的基本原理可知,親緣關系越近,越早聚合為一類。反之,親緣關系越遠的品種,越晚聚合為一類。

采用系統聚類法,利用歐氏距離,離差平方和法,根據14個果實表型性狀指標對154份西瓜種質材料進行聚類分析,圖1給出了西瓜種質資源類型間的距離及遺傳差異。從圖1可看出,當歐氏距離為32.89時,可分為4類。

第1類為大果類型,包括37份材料,21份來自于河南,其余引自于外地。果實形狀為圓柱形,果皮顏色以綠色和淺綠色為主,5份綠白皮,1份黃綠,2份墨綠和2份深綠色;果皮覆紋顏色以綠色和深綠色為主,覆紋形狀有19份為網條性狀,13份為齒條;瓤色以紅色為主,2份黃色,5份大紅色,果肉質地以脆為主,7份沙瓤,7份軟瓤;縱徑平均為22.3 cm,橫徑平均為13.4 cm,果形指數1.5,果皮厚度平均為1.1 cm,中心可溶性固形物含量平均為12.2%,邊部可溶性固形物含量平均為9.0%,中邊可溶性固形物含量差平均值為2.9,單瓜質量4.7 kg。

第2類為中大果類型,包括56份材料,38份來自河南,其余引自外地。果實圓形,果皮顏色多樣化,以綠色為主,其他有淺綠、深綠、墨綠、綠白色,還有1份黃色果皮;果皮覆紋顏色以墨綠色為主,其次為深綠色和綠色,1份深黃色;覆紋形狀有條帶、網條、齒條3種;瓤色主要有紅色和大紅色;果肉質地以脆為主,其次為硬;縱徑平均值16.2 cm,橫徑平均值15.9 cm,果皮厚度平均值為1.0 cm,果皮較厚,中心可溶性固形物含量平均為11.7%,糖度高,邊部可溶性固形物含量平均為8.9%,單瓜質量4.5 kg。

第3類為中果類型,包括8份材料,5份來自于河南,2份來自于湖南,1份來自于臺灣。果實形狀為圓形,只有1份為橢圓形;果皮墨綠色;無果皮覆紋;瓤色多樣化,有紅色、大紅色、粉紅、橘紅、淺黃、深黃色;果肉質地脆;縱徑平均為17.6 cm,橫徑平均為16.4 cm,果皮厚度平均為1.1 cm,果皮厚,中心可溶性固形物含量平均為12.4%,邊部可溶性固形物含量平均為9.1%,單瓜質量4.2 kg。

第4類為小果類型,包括53份材料,12份來自于河南,其他均為外地引進,大部分引自臺灣。果實形狀為28份圓形,其余為橢圓形和橄欖形;果皮顏色以淺綠色和綠色為主,其余有深綠、墨綠,淺黃色1份;果皮覆紋顏色以墨綠色為主,有27份,其次為綠色,有16份,其余為深綠色,有12份,黃色1份;覆紋形狀有條帶、齒條、網條;瓤色多樣化,有乳白色、淺黃、黃色、深黃色、桃紅色、粉紅色、紅色、大紅色;果肉質地以脆為主;縱徑平均為17.0 cm,橫徑平均為13.6 cm,果皮厚度平均為0.8 cm,果皮較薄,中心可溶性固形物含量平均為12.4%,糖度高,邊部可溶性固形物含量平均為10.3%,中邊可溶性固形物含量差很小,單瓜質量2.7 kg。

3 討 論

從上述結果可以看出,利用數量分類法能直觀反映品種間親緣關系,雜交育種時,可參考聚類分析鑒定結果,明確區分不同類群資源特點,直接選擇所需親本或親緣關系較遠的西瓜材料雜交,有利于擴大遺傳背景,提高育種效率。

變異系數是反映生物體形態遺傳多樣性的重要遺傳參數,潘存祥等[12]對國內外783份西瓜種質資源24個表型性狀進行了遺傳研究,發現24個表型性狀多樣性指數平均值為1.68,平均變異系數為31.19%,其中種子覆紋顏色變異系數最大(70.90%),第一雌花節位變異系數最小(0.48%)。紀海波等[9]對768份西瓜種質資源的13個數量性狀和11個質量性狀進行研究,其平均多樣性指數為1.6,變異系數均值為32.41%,其中種子斑紋形狀的變異系數最大(73.32%),開放雌花到莖生長點距離的變異系數最?。?.32%)。郭祿芹等[13]對167份西瓜種質材料的的遺傳多樣性進行分析,質量性狀的多樣性指數平均值為1.16,數量性狀多樣性指數平均值為1.71,其18個數量性狀的變異系數變化范圍為16.43%~90.73%,平均值為34.12%。而本研究中多樣性指數平均值為1.71,與前者的研究結果基本相同,變異系數平均值為28.38%,變異系數的變幅范圍為8.4%~62.51%,介于前者的研究中間,體現了試驗所用西瓜材料的多樣性,中心可溶性固形物含量變異系數最低,為8.4%,這可能與目前我們以高糖西瓜作為主要育種目標有很大關系。

通過對西瓜果實8個數量性狀進行主成分分析,可將8個性狀歸為4個主成分。4個主成分特征值累計貢獻率達到83.8%,4個主成分可代替原來8個指標對西瓜果實性狀進行評價。趙衛星等[14]將7個西瓜果實性狀歸納為4個主成分,累計貢獻率為93.94%,并認為果實質量、果皮厚、中心可溶性固形物含量、中邊可溶性固形物含量差是影響所選西瓜品種數量性狀的主要因子,果實邊部可溶性固形物含量和縱徑是次要因子。王志強等[15]將6個數量性狀綜合成為3個主成分,且這3個主成分構成的信息可解釋總信息量的84.30%,認為對果實數量性狀影響最大的是果實質量、中心可溶性固形物含量和果皮硬度3個性狀指標。朱子成等[16]的研究表明,前3個主成分特征值的累計貢獻率達87.84%,認為對西瓜果實影響最大的是單果質量、可溶性固形物含量和果形指數。而筆者研究認為,果實邊部可溶性固形物含量、單瓜質量、縱徑、果形指數、中心可溶性固形物含量、中邊可溶性固形物含量差、橫徑是影響所選西瓜品種果實性狀的主要因子,可作為劃分品種類群的依據,因此在育種選擇過程中,高品質是我們以后育種的主要目標。

利用聚類分析對品種資源進行分析,可以真實反映品種的綜合性狀,陳萍等[10]根據單果質量、果皮厚度和含糖量對119份西瓜種質材料進行聚類分析,當歐氏距離絕對值最長距離等于5.5時,可分成10類。趙衛星等[14]以主成分中4個主要因子為指標進行聚類分析,將56份西瓜品種分為小果型、中小果型、大果型、中大果型4類。郭祿芹等[13]根據29個表型性狀,利用UPGMA法對167份種質材料進行聚類分析,165號材料單獨位于聚類圖的最外層,說明與其他材料的親緣關系較遠。除該材料外,其他西瓜種質材料可分為4類。不同來源的種質相互交織在一起,并沒有按照來源區分開,說明同一類群間種質材料的遺傳背景相似度的高低和親緣關系的遠近與來源地無顯著相關性,也可能是由于表型數據易受環境影響而造成的誤差所致。筆者基于表型數據的聚類結果與陳萍、趙衛星等[10,14]的研究結果類似。而許彥賓等[17]利用SSR標記分析34份西瓜自交系的遺傳關系認為,所有材料的平均相似系數為0.56,大體可分為3類(A、B和C)。其中,B、C 2類包含4份自交系(X36、X48、X84和X89),其果實較小,均從中國新疆、甘肅等地的品種中自交分離所得,具有較高的遺傳多樣性。A類包含30份自交系,除X38外均為大果型材料,平均相似系數高達0.78,西瓜自交系之間的遺傳分化較小,遺傳背景較為狹窄。筆者基于14個果實表型性狀指標對154份西瓜種質材料進行聚類分析,當歐氏距離為32.89時,可分為4類。其中第4類為小果類型,皮薄、瓜個小、糖度高、品質和風味好,瓤色和皮色類型也較豐富,符合當下人們生活需求,在以后選育工作中具有重要利用價值。

筆者通過對西瓜果實性狀的鑒定與分析,明確了西瓜果實性狀的遺傳多樣性,其研究結果對于西瓜資源表型性狀精準鑒定、種質創新及雜交育種親本選配等方面具有積極的作用。

參考文獻

[1] 顏國榮,王威,馬艷明,等.甜瓜種質資源若干果實數量性狀的差異研究及分級探討[J].新疆農業科學,2006,4(34):328-331.

[2] 周俊國,李新崢,朱月林,等.蔓生型南瓜資源部分植物學數量性狀的評價指標探討[J].植物遺傳資源學報,2007,8(1):30-34.

[3] 馬蔚紅,謝江輝,武紅霞,等.杧果種質資源果實數量性狀評價指標探討[J].果樹學報,2006,23(2):218-222.

[4] 馬雙武,王吉明,邱江濤.我國西瓜甜瓜種質資源收集保存現狀及建議[J].中國西瓜甜瓜,2003,16(5):17-19.

[5] 范敏,宮國義,張瑞麟,等.美國資源庫西瓜種質的初步觀察與數量分類[J].中國西瓜甜瓜,2004,17(4):1-3.

[6] 李艷梅,段會軍,馬峙英.西瓜種質資源的遺傳多樣性及親緣關系的AFLP分析[J].華北農學報,2007,22(S1):177-180.

[7] 張愛萍,王曉武,張岳莉,等.西瓜種質資源遺傳多樣性的SRAP分析[J].中國農學通報,2008,24(4):115-120.

[8] 尚建立,王吉明,郭琳琳,等.西瓜種質資源主要植物學性狀的遺傳多樣性及相關性分析[J].植物遺傳資源學報,2012,13(1):11-15.

[9] 紀海波,張玉鑫,李玉明,等.西瓜種質資源主要性狀的表型多樣性[J].西北農林科技大學學報(自然科學版),2013,41(8):155-160.

[10] 陳萍,祝佳強,祝雙田,等.西瓜種質材料果實性狀的變異性和聚類分析[J].北方園藝,2010(17):182-184.

[11] 馬雙武,劉君璞.西瓜種質資源描述規范和數據標準[M].北京:中國農業出版社,2005.

[12] 潘存祥,許勇,紀海波,等.西瓜種質資源表型多樣性及聚類分析[J].植物遺傳資源學報,2015,16(1):59-63.

[13] 郭祿芹,趙世豪,朱華玉,等.167份西瓜種質材料的遺傳多樣性分析[J].中國瓜菜,2018,31(1):5-11.

[14] 趙衛星,李曉慧,徐小利,等.西瓜品種果實性狀的變異性及主成分和聚類分析[J].河南農業科學,2017,46(1):105-109.

[15] 王志強,劉聲鋒,郭守金,等.西瓜種質資源種子性狀的遺傳多樣性和相關性分析[J].中國農學通報,2013,29(25):205-208.

[16] 朱子成,張子幸,裴爽,等.西瓜主要數量性狀的相關性和主成分分析[J].北方園藝,2017(12):35-39.

[17] 許彥賓,王艷玲,溫新惠,等.利用SSR標記分析西瓜自交系的遺傳關系[J].中國瓜菜,2017,30(4):10-13.

主站蜘蛛池模板: 国产91高跟丝袜| 一级毛片基地| 免费又爽又刺激高潮网址 | 一级一毛片a级毛片| 婷婷综合在线观看丁香| 91在线国内在线播放老师| 亚洲天堂网视频| 国产福利一区在线| 国产成人亚洲无码淙合青草| 成·人免费午夜无码视频在线观看| 免费Aⅴ片在线观看蜜芽Tⅴ | 性网站在线观看| 伊人久久大线影院首页| 伊人久久久久久久| 亚洲一欧洲中文字幕在线| 欧美国产中文| 日日拍夜夜嗷嗷叫国产| 久久无码av三级| 国产香蕉97碰碰视频VA碰碰看| 欧美成人免费一区在线播放| 国产va视频| 视频二区国产精品职场同事| 黄色网站不卡无码| 欧美成人影院亚洲综合图| 欧美日本视频在线观看| 色婷婷久久| 国产波多野结衣中文在线播放| 国产成人做受免费视频| 少妇精品久久久一区二区三区| 无码一区中文字幕| 日韩精品亚洲人旧成在线| 国产69精品久久久久孕妇大杂乱| 真人免费一级毛片一区二区| 九九视频免费在线观看| 中国国产高清免费AV片| 日韩高清无码免费| 青青青国产在线播放| 韩国v欧美v亚洲v日本v| 国产精品成人不卡在线观看 | 99精品在线看| 久久久久人妻一区精品| 免费中文字幕一级毛片| 亚洲中文字幕手机在线第一页| 91在线国内在线播放老师| 午夜电影在线观看国产1区| 久草视频精品| 日韩无码视频播放| 国产av剧情无码精品色午夜| a级毛片在线免费观看| 波多野结衣中文字幕一区| 欧美亚洲香蕉| 97精品久久久大香线焦| 成人福利在线视频| 无码av免费不卡在线观看| 精品久久久久成人码免费动漫 | 国产精品 欧美激情 在线播放| 欧美日韩精品在线播放| 91探花国产综合在线精品| 国产美女主播一级成人毛片| 91久久国产成人免费观看| 国产亚洲男人的天堂在线观看| 91福利国产成人精品导航| 国产福利影院在线观看| 欧美成人一级| 免费无码一区二区| 久久国产精品麻豆系列| 国产精品开放后亚洲| 99ri国产在线| 亚洲精品桃花岛av在线| 国产成熟女人性满足视频| 全部免费毛片免费播放| 精品国产亚洲人成在线| 国产办公室秘书无码精品| 国产黑丝视频在线观看| 久久香蕉国产线| 日韩大片免费观看视频播放| 无码在线激情片| 亚洲人免费视频| 国产一区二区精品高清在线观看| 精品国产毛片| 国产精品久久久久无码网站| 久久精品中文字幕少妇|