• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    二棱大麥種質(zhì)的綜合評價

    2018-05-11 09:36:14劉亞楠郭寶健許如根
    麥類作物學報 2018年4期
    關鍵詞:穗長粒重穗數(shù)

    劉亞楠,朱 娟,呂 超,郭寶健,許如根

    (1.揚州大學江蘇省作物遺傳生理國家重點實驗室培育點/糧食作物現(xiàn)代產(chǎn)業(yè)技術協(xié)同創(chuàng)新中心/教育部植物功能基因組學重點實驗室/揚州大學大麥研究所,江蘇揚州 225009;2.江蘇天豐種業(yè)有限公司,江蘇淮安 223001)

    大麥(HordeumvulgareL.)是世界上重要的禾谷類作物,依據(jù)穗的棱型,可分為二棱大麥和六棱大麥。相比六棱大麥,二棱大麥籽粒品質(zhì)更適合用作釀造啤酒的原料。我國的二棱大麥主要分布于北方的春大麥區(qū)和長江中下游的冬大麥區(qū),由于生態(tài)環(huán)境的不同,二棱大麥種質(zhì)性狀存在差異[1]。如江蘇不同時期育成的二棱啤酒大麥品(系)先后經(jīng)歷了早熟、抗病、高產(chǎn)、優(yōu)質(zhì)等階段,在產(chǎn)量、抗倒性、抗病性方面均發(fā)生明顯的改善。大麥產(chǎn)量性狀與株型關系密切[2]。二棱大麥的千粒重與產(chǎn)量的正相關性最大,因而其高產(chǎn)應在足穗基礎上,以較高的千粒重為保障[3]。大麥種質(zhì)資源是大麥種質(zhì)創(chuàng)新、生物學研究的基礎[4],對二棱大麥種質(zhì)資源的遺傳多樣性進行合理的多元分析與評價,有助于挖掘二棱大麥特異種質(zhì)資源,充分利用其優(yōu)良基因,防止育種骨干親本的單一性,加快育種進程。

    通過表型性狀[5-7]、分子標記[8-14]及生化標記[15]等不同水平鑒定,可對二棱大麥種質(zhì)的遺傳多樣性進行分析,從而篩選出農(nóng)藝、品質(zhì)、抗耐性等性狀優(yōu)異的種質(zhì)資源[16-17]。表型性狀的評價方法最直觀,易被育種家掌握。大麥不同性狀代表了不同方面的特性,而且性狀之間往往存在不同程度的相關性,因而采用多性狀對種質(zhì)資源進行綜合評價更具合理性。周 偉等[18]利用主成分分析篩選了二棱大麥種質(zhì)的農(nóng)藝性狀重要組分,通過系統(tǒng)聚類的方法獲得特異性種質(zhì),可作為不同類型大麥的育種親本或中間材料。楊樹明等[19]對國內(nèi)外不同地區(qū)啤酒大麥品種進行農(nóng)藝性狀的鑒定與聚類,認為不同地區(qū)啤酒大麥多樣性分布具有明顯生態(tài)地域性。近年來,麻類、大蒜、小麥、棉花等作物研究逐漸將基于主成分的二維排序分析引入到種質(zhì)鑒定與綜合評價[20-24],揭示了作物種質(zhì)的性狀表現(xiàn)及種質(zhì)間的差異。此外,也有學者利用主成分分析與多元線性回歸構(gòu)建核心種質(zhì)的評價方程對作物種質(zhì)多性狀進行綜合評價[25-27]。這些分析方法在大麥中卻未見報道。本研究通過結(jié)合主成分分析、二維排序分析、回歸分析等多元分析與評價方法,明確二棱大麥種質(zhì)遺傳多樣性,發(fā)掘出綜合育種目標優(yōu)良或特異性種質(zhì)資源,以期為二棱大麥育種研究提供科學依據(jù)。

    1 材料與方法

    1.1 試驗材料

    以揚州大學大麥研究所收集篩選的97份來自國內(nèi)外的二棱大麥品種(系)為參試材料。材料名稱及來源見表1。

    1.2 試驗設計

    2014年秋將試驗材料種植于鹽城市大中農(nóng)場,采用隨機區(qū)組設計,每個材料播種2行,行長1.2 m,行距0.2 m,株距3 cm,每行點播40粒,3次重復。2015年秋將試驗材料分別種植于揚州大學新農(nóng)場及鹽城市大中農(nóng)場,同2014年同規(guī)格種植,常規(guī)栽培管理。

    1.3 性狀測定

    大麥成熟后,每個材料取5株代表植株,按照《大麥種質(zhì)資源描述規(guī)范和數(shù)據(jù)標準》測定株高、穗長、穗下節(jié)間長、主穗粒數(shù)、單株穗數(shù)、單株生物量、單株粒重和千粒重。

    1.4 統(tǒng)計分析

    使用Excel 2007進行基礎數(shù)據(jù)的整理,計算各性狀的均值、最大值、最小值、標準差、變異系數(shù)、遺傳多樣性指數(shù),使用SPSS 21.0軟件進行方差分析、相關分析、二維排序分析、回歸分析,使用Matlab軟件進行主成分分析;使用R語言進行基于馬氏距離類平均法聚類。

    通過計算參試材料的平均數(shù)(μ)和標準差(σ),將考察的目標性狀均劃分為10級,從第1級[xi<(μ-2σ)]到第10級[xi>(μ+2σ)],每間隔0.5σ為一級,計算每一級相對頻率Pi,從而得到遺傳多樣性指數(shù)[28](Shannon-Wiener diversity index,H′)。H′=-∑PilnPi。式中Pi為某性狀第i級別材料份數(shù)占總份數(shù)的百分比。

    2 結(jié)果與分析

    2.1 二棱大麥種質(zhì)主要性狀的表現(xiàn)

    方差分析(表2)表明,二棱大麥株高、穗下節(jié)間長、穗長、單株生物量、主穗粒數(shù)、千粒重、單株粒重7個性狀的環(huán)境、基因型及基因型與環(huán)境互作效應均達到極顯著水平;單株穗數(shù)的環(huán)境和基因型效應也達極顯著水平,但環(huán)境與基因型互作效應不顯著。二棱大麥種質(zhì)穗長、單株穗數(shù)、單株粒重、單株生物量在2年3試點變異均較大,變異系數(shù)為14.00%~26.99%,而株高、穗下節(jié)間長、主穗粒數(shù)、千粒重在二棱大麥種質(zhì)中變異范圍較小,變異系數(shù)為9.42%~14.74%(表3)。其中,2016年揚州新農(nóng)場試點的變異最大,而鹽城試點兩年的變異系數(shù)較小且表現(xiàn)較為一致。單株穗數(shù)和生物量受環(huán)境影響最大,在兩年3試點間差異顯著。不同性狀中,穗下節(jié)間長、主穗粒數(shù)的遺多樣性指數(shù)最高,分別為2.06和2.07;單株生物量的遺傳多樣性指數(shù)最低,為1.81。

    表1二棱大麥種質(zhì)資源信息表

    Table1Informationoftwo-rowedbarleygermplasm

    編號Number品種(系)Variety(line)來源Source編號Number品種(系)Variety(line)來源Source1蘇農(nóng)22Sunong22江蘇Jiangsu50如東5218Rudong5218江蘇Jiangsu2浙農(nóng)大7號Zhenong7浙江Zhejiang51連9723Lian9723江蘇Jiangsu3蘇農(nóng)91?712Sunong91?712江蘇Jiangsu52紅07?456Hong07?456黑龍江Heilongjiang4揚農(nóng)啤4號Yangnongpi4江蘇Jiangsu53紅08?718Hong08?718黑龍江Heilongjiang5創(chuàng)大麥3號Chuangmai3湖北Hubei54浙啤33Zhepi33浙江Zhejiang6揚農(nóng)啤6號Yangnongpi6江蘇Jiangsu55鄂大麥507Edamai507湖北Hubei7蘇農(nóng)16Sunong16江蘇Jiangsu56駐2005?7?9Zhu2005?7?9河南Henan8云啤1號Yunpi1云南Yunnan57蘇啤4號Supi4江蘇Jiangsu9鹽引1號Yanyin1日本Japan58西南86?911Xinan86?911四川Sichuan10矮稈?1Aigan?1江蘇Jiangsu59川52209Chuan52209四川Sichuan1132122湖北Hubei60揚引02Yangyin02江蘇Jiangsu12浙大96?6Zheda96?6浙江Zhejiang61鹽引1號×浙農(nóng)大3號Yanyin1×Zhenongda3江蘇Jiangsu13莆大麥5號Pudamai5福建Fujian62QS江蘇Jiangsu14鄂大麥6號Edamai6湖北Hubei63駐3Zhu3河南Henan15駐05?043?2?1Zhu05?043?2?1河南Henan64駐7Zhu7河南Henan16蘇農(nóng)6472Sunong6472江蘇Jiangsu65周選1號Zhouxuan1河南Henan17蘇引麥3號Suyinmai3上海Shanghai66紅09?784Hong09?784黑龍江Heilongjiang182004?日引1號2004?Riyin1日本Japan67華大麥5號Huadamai5湖北Hubei192004?日引3號2004?Riyin3日本Japan68云啤2號Yunpi2云南Yunnan20揚農(nóng)啤5號Yangnongpi5江蘇Jiangsu69云啤3號Yunpi3云南Yunnan21Nasonijo日本Japan70云啤4號Yunpi4云南Yunnan22揚啤1號Yangpi1江蘇Jiangsu71云啤5號Yunpi5云南Yunnan23莆大麥9號Pudamai9福建Fujian72云啤6號Yunpi6云南Yunnan24揚農(nóng)啤8號Yangnongpi8江蘇Jiangsu73單95168×鹽94130Dan95168×Yan94130江蘇Jiangsu25Harrinonijo日本Japan74滬01?2946Hu01?2946上海Shanghai26浙秀12Zhexiu9浙江Zhejiang75云啤7號Yunpi7云南Yunnan27優(yōu)質(zhì)1號Youzhi1江蘇Jiangsu76云啤9號Yunpi9云南Yunnan28駐9505?1?3Zhu9505?1?3河南Henan77法瓦維特Favorit匈牙利Hungary29蘇啤3號Supi3江蘇Jiangsu78駐5Zhu5河南Henan3091單291dan2江蘇Jiangsu79華2328Hua2328湖北Hubei31陜西釀造2條Shaanxiliangzao2陜西Shanxi80美97?1455Mei97?1455美國USA32通麥10號Tongmai10江蘇Jiangsu81駐6Zhu6河南Henan33鹽99175Yan99175江蘇Jiangsu82華大麥6號Huadamai6湖北Hubei34揚飼麥1號×鑒135Yangjian135江蘇Jiangsu83華大麥7號Huadamai7湖北Hubei35浙35?21Zhe35?21浙江Zhejiang84Gairdner澳大利亞Australia36鄂大麥9706Edamai9706湖北Hubei85Frankin澳大利亞Australia37駐06049?1Zhu06049?1河南Henan86T98189澳大利亞Australia38花11Hua11上海Shanghai87Suyin27日本Japan39花22Hua22上海Shanghai88揚農(nóng)啤9號Yangnongpi9江蘇Jiangsu40揚飼麥3號Yangsimai3江蘇Jiangsu89揚農(nóng)啤10號Yangnongpi10江蘇Jiangsu41揚農(nóng)啤7號Yangnongpi7江蘇Jiangsu90揚農(nóng)啤2號Yangnongpi2江蘇Jiangsu42蘇花2號Suhua2江蘇Jiangsu91揚農(nóng)啤11號Yangnongpi11江蘇Jiangsu43Haronanijo日本Japan92揚農(nóng)啤12號Yangnongpi12江蘇Jiangsu44駐96015?5?3Zhu96015?5?3河南Henan93黃長芒Huangchangmang湖北Hubei45余姚向天二棱Yuyaoxiangtianerleng浙江Zhejiang94滬1154Hu1154上海Shanghai46莆大麥8號Pudamai8福建Fujian95揚0187Yang0187江蘇Jiangsu47英7204Ying7204英國UK961430R江蘇Jiangsu48駐97022Zhu97022河南Henan97泰興9425Taixing9425江蘇Jiangsu49如東6109Rudong6109江蘇Jiangsu

    根據(jù)遺傳多樣性,將二棱大麥8個農(nóng)藝性狀分為10級(圖1)。株高、主穗粒數(shù)及千粒重主要集中在第4~第7級:株高集中主要分布在85.15~105.85 cm范圍,主穗粒數(shù)主要集中在26.4~32.80粒范圍,千粒重主要集中在38.69~47.53 g范圍,三個性狀主要分布范圍的材料數(shù)分別占總體的74.23%、69.07%和75.26%;穗下節(jié)間長、單株穗數(shù)、單株粒重及單株生物量主要集中在第3~第7級:穗下節(jié)間長集中分布在30.08~39.04 cm范圍,單株穗數(shù)集中分別在7.60~10.86個范圍,單株粒重集中分布7.36~11.21 g范圍,單株生物量集中分布在19.41~28.35 g范圍,4個性狀主要分布范圍的材料分別占總體的79.38%、80.41%、80.41%和90.72%。穗長主要集中分布在第4、第5級,即6.07~6.62 cm 范圍,材料數(shù)占總體的50.52%。這表明參試二棱大麥品種主要性狀分布較為集中,存在主要性狀的趨同性問題。

    2.2 二棱大麥種質(zhì)的聚類分析

    基于馬氏距離的類平均法,使用R語言對二棱大麥種質(zhì)的農(nóng)藝性狀綜合表現(xiàn)進行聚類分析,結(jié)果(圖2)表明,97份二棱大麥種質(zhì)在馬氏距離30處可聚為8類。

    第Ⅰ類包括38份種質(zhì),其中江蘇16份,河南8份,湖北4份,浙江4份,上海3份,四川2份,福建1份及日本1份,該類材料來源主要集中在長江中下游冬大麥區(qū),以近年來各省選育推廣的品種為主。該類種質(zhì)植株較矮,穗下節(jié)間較長,穗略短,單株穗數(shù)略少,主穗粒數(shù)較多,千粒重和單株粒重中等,單株生物量略低。

    PH:株高;ILBS:穗下節(jié)間長;SL:穗長;GNPS:主穗粒數(shù);SNPP:單株穗數(shù);TKW:千粒重;GWPP:單株粒重;BWPP:單株生物量。下表和圖2同。1~10不同顏色表示10個分組組別。

    PH:Plant height; ILBS:Internode length below spike; SL:Spike length; GNPS:Grain number per spike; SNPP:Spike number per plant; TKW:Thousand-kernel weight; GWPP:Grain weight per plant.The same in other tables and figures.Ten different colors represent ten groups,respectively.

    圖1二棱大麥種質(zhì)資源主要農(nóng)藝性狀分布

    Fig.1Distributionofagronomictraitsintwo-rowedbarleygermplasm

    表2二棱大麥種質(zhì)資源農(nóng)藝性狀的方差分析

    Table2ANOVAofagronomictraitsoftwo-rowedbarleygermplasm

    變異來源Variationsource株高PH穗下節(jié)間長ILBS穗長SL單株生物量GWPP主穗粒數(shù)GNPS單株穗數(shù)SNPP千粒重TKW單株粒重GWPP環(huán)境Environment344.85??29.21??180.19??199.07??18.39??82.83??31.86??28.85??基因型Genotype23.49??26.19??56.16??4.00??10.97??4.62??15.01??3.35??互作Interaction1.45??2.43??2.02??2.36??1.30??2.202.42??2.31??

    *和**分別表示在0.05和0.01水平上顯著。下同。

    * and ** mean significant difference at 0.05 and 0.01 levels,respectively.The same in other tables.

    表3二棱大麥種質(zhì)資源主要農(nóng)藝性狀在不同試點的表現(xiàn)

    Table3Traitsvariationoftwo-rowedbarleygermplasmindifferentlocations

    性狀Trait2015年鹽城點Yanchenglocationin2015均值Mean變化范圍Range變異系數(shù)CV/%2016年鹽城點Yanchenglocationin2016均值Mean變化范圍Range變異系數(shù)CV/%2016年揚州新農(nóng)場點Yangzhoulocationin2016均值Mean變化范圍Range變異系數(shù)CV/%遺傳多樣性指數(shù)H′株高PH/cm101.03a72.40~127.8010.0687.67b64.73~114.1710.8397.81a72.83~150.3013.232.01穗下節(jié)間長ILBS/cm34.92b26.34~41.799.6235.30ab24.27~43.7611.4636.16a21.68~46.0211.792.06穗長SL/cm6.90b5.01~10.4516.757.11b5.28~10.0715.077.48a5.95~11.4815.891.89主穗粒數(shù)GNPS30.36a22.27~37.7310.6028.98b21.07~46.8013.7229.51ab22.87~44.8711.612.07單株穗數(shù)SNPP9.61b6.97~16.4717.598.69c5.87~13.4819.1210.37a6.43~18.9019.052.01千粒重TKW/g43.07ab33.61~53.399.4244.15a31.53~53.619.9641.86b26.61~53.4814.741.98單株粒重GWPP/g9.84a4.67~16.1320.108.96b5.31~14.7519.7510.12a4.11~19.9526.842.03單株生物量BWPP/g23.55b16.32~52.6818.7521.10c16.07~27.9914.0029.67a17.03~55.8225.391.81

    數(shù)據(jù)均值后面的小寫字母表示各性狀在3試點間的差異顯著(P<0.05)。

    The letter following the average indicate significant difference at different locations.

    第Ⅱ類包括13份種質(zhì),其中云南5份,黑龍江1份,國外引進種質(zhì)7份,該類種質(zhì)以春大麥為主,株高、穗長及單株穗數(shù)居中,穗下節(jié)間較短,千粒重較低,主穗粒數(shù)較少,單株粒重與生物量最低,在江蘇生態(tài)條件的綜合表現(xiàn)較差。

    第Ⅲ類包括28份種質(zhì),其中江蘇8份,湖北5份,浙江、云南、黑龍江、河南種質(zhì)各2份,陜西、福建各1份,日本3份,英國1份,以長江中下游冬大麥品種為主。該類種質(zhì)表現(xiàn)為株高適中,穗下節(jié)間較長,主穗粒數(shù)較多,單株穗數(shù)、千粒重、單株粒重及生物量適中,穗較短,綜合表現(xiàn)中等,與第一類比較近。

    第Ⅳ類有10份種質(zhì),包括江蘇的8份材料、上海的滬1154及福建的莆大麥8號。該類種質(zhì)表現(xiàn)為矮稈,穗下節(jié)間長適中,單株穗數(shù)較多,千粒重與單株粒重較高,穗長適中,生物量略低,綜合表現(xiàn)較好。

    第Ⅴ類包括2份種質(zhì),即云啤9號和美國引進美97-1455。其特征為株高適中,穗下節(jié)間較短,穗長較長,單株穗數(shù)中等,穗粒數(shù)最多,但單株粒重較低,千粒重與生物量最低。

    第Ⅵ類包括3份種質(zhì),即1430R、揚引02及周選1號,均為引進改良品種,其植株較高,穗下節(jié)間與穗長較長,著粒密度稀,穗粒數(shù)適中,單株穗數(shù)較多,單株粒重較高,千粒重最高。

    第Ⅶ類僅包括來自湖北的裸大麥地方品種黃長芒,其高稈,穗下節(jié)間與穗長最長,多穗,主穗粒數(shù)與千粒重適中,單株粒重與生物量最高,單株性狀綜合表現(xiàn)最好。

    第Ⅷ類包含江蘇的優(yōu)質(zhì)1號和上海的滬01-2946,二者均高稈,穗長與穗下節(jié)間較長,主穗粒數(shù)最多,千粒重較高,單株穗數(shù)與單株粒重,生物量表現(xiàn)中等。

    2.3 二棱大麥主要性狀的主成分分析

    二棱大麥種質(zhì)各性狀間存在不同程度的相關性(表4),直接依據(jù)這些表型評價種質(zhì)有失偏頗。為消除性狀間相關性的影響,采用主成分分析進行種質(zhì)評價。按照累計貢獻率大于85%的原則,提取前4個主成分來概括全部性狀的大部分生物學信息(表5)。

    第一主成分的貢獻率為33.19%,千粒重和單株粒重的特征向量值最大,且均為正值,其凝聚的生物學信息跟粒重直接相關,可概括為粒重因子,粒重增加有利于高產(chǎn),第一主成分應大為好。

    第二主成分的貢獻率為22.39%,單株穗數(shù)和主穗粒數(shù)的特征向量值最大,單株穗數(shù)值為正,主穗粒數(shù)值為負,將其生物學信息概括為穗部因子,分蘗增多使生物量顯著增加,不利于大穗大粒的形成,單株穗數(shù)與千粒重呈不顯著負相關,該主成分過大不符合二棱啤酒大麥的籽粒要求,應適中為好。

    第三主成分貢獻率為18.51%,穗長的特征向量值最大,且為正值,故該主成分涵蓋的生物學信息為穗長因子,穗子增長,主穗粒數(shù)增加,影響千粒重的提高,但其負相關效應不顯著,第三主成分應適當偏大為好。

    第四主成分貢獻率為11.12%,株高、穗下節(jié)間長及生物量的特征向量值較大,且為正值,將其生物學信息概括為株型因子,穗下節(jié)間增長有利于產(chǎn)量、穗粒數(shù)、單株穗數(shù)、千粒重的增加,但植株過高易造成倒伏,第四主成分應略小為好。

    圖2 二棱大麥種質(zhì)聚類圖

    Table4Correlativecoefficientsbetweenagronomictraitsoftwo-rowedbarleygermplasm

    性狀Trait株高PH穗下節(jié)間長ILBS穗長SL主穗粒數(shù)GNPS單株穗數(shù)SNPP千粒重TKW單株粒重GWPP穗下節(jié)間長ILBS0.320??穗長SL0.391??-0.087主穗粒數(shù)GNPS0.1730.364??0.282??單株穗數(shù)SNPP-0.004-0.231?0.150-0.449??千粒重TKW-0.0400.425??-0.0580.040-0.198單株粒重GWPP-0.0110.308??0.237?0.1460.425??0.517??單株生物量BWPP0.469??0.345??0.393??0.1170.390??0.314??0.704??

    表5二棱大麥種質(zhì)資源主成分分析的特征值及貢獻率

    Table5EigenvalueandcontributionrateofPCAoftwo-rowedbarleygermplasm

    性狀TraitPV1PV2PV3PV4株高PH-0.038-0.0180.1930.988穗下節(jié)間長ILBS0.357-0.291-0.1650.553穗長SL-0.0110.0130.7260.204主穗粒數(shù)GNPS0.099-0.5660.4070.071單株穗數(shù)SNPP0.0680.6720.181-0.003千粒重TKW0.506-0.131-0.159-0.036單株粒重GWPP0.5240.1990.293-0.057單株生物量BWPP0.3730.2350.3420.506特征值E2.6551.7911.4800.890百分率CR/%33.19322.39318.50611.121累計百分率CCR/%33.19355.58674.09385.213

    2.3.1 基于主成分的二維排序分析

    基于二棱大麥種質(zhì)資源主成分分析的結(jié)果,選取前4個主成分得分值為指標,以第一主成分為橫坐標,分別以第二、三、四主成分為縱坐標繪制二維排序圖(圖3)。第一主成分粒重因子與第二主成分穗部因子二維排序圖(圖3a)中應以分布于的第一、四象限橫坐標軸附近的種質(zhì)為好,主要包括周選1號(65)、蘇啤4號(57)、蘇農(nóng)22(1)、優(yōu)質(zhì)1號(27)、駐9505-1-3(28)、鄂大麥6號(14),這些種質(zhì)大粒、豐產(chǎn)性較好,且單株穗數(shù)與穗粒數(shù)協(xié)調(diào)。第一主成分粒重因子與第三主成分穗長因子二維排序圖(圖3b)中應以第一象限的種質(zhì)較好,主要包括黃長芒(93)、周選1號(65)、揚引02(60)、浙啤33(54)、華2328(79)、揚農(nóng)啤9號(88),這類種質(zhì)表現(xiàn)為長穗、穗粒數(shù)多且粒大、豐產(chǎn)性好。第一主成分粒重因子與第四主成分株型因子二維排序圖(圖3c)中應以第四象限種質(zhì)為好,主要包括揚農(nóng)啤9號(88)、揚農(nóng)啤2號(90)、揚農(nóng)啤12號(92)、揚農(nóng)啤7號(41)、揚農(nóng)啤10號(89)、揚農(nóng)啤6號(6)、浙啤33(54)、滬1154(94)、揚農(nóng)啤5號(20),該類種質(zhì)表現(xiàn)為矮稈抗倒伏且粒大、豐產(chǎn)性好。

    2.3.2 基于主成分的M值綜合評價

    對二棱大麥種質(zhì)資源8個主要農(nóng)藝性狀數(shù)值進行標準化處理,獲得前4個主成分得分,將這4個主成分得分作歸一化處理,以貢獻率計算各主成分權重系數(shù)分別為0.389 5、0.262 8、0.217 1和0.130 6,利用前4個主成分的得分值與其對應的權重乘積之和M值作綜合指標,對各種質(zhì)材料進行綜合評價,M值越高,綜合性狀表現(xiàn)越好。

    圖3 二棱大麥種質(zhì)資源主成分二維排序圖

    二棱大麥種質(zhì)材料排在前十位的分別為黃長芒(93)、揚引02(60)、1430R(96)、周選1號(65)、揚農(nóng)啤7號(41)、如東6109(49)、滬1154(94)、揚農(nóng)啤9號(88)、蘇啤4號(57)、揚農(nóng)啤12號(92)。其中黃長芒為裸大麥地方品種,其他均為江蘇、上海地區(qū)選育的啤用大麥品種。這些種質(zhì)綜合表現(xiàn)較好,可直接在生產(chǎn)上加以利用。

    2.3.3 二棱大麥種質(zhì)多性狀綜合評價的最優(yōu)回歸方程

    結(jié)合綜合得分M值和8個主要農(nóng)藝性狀,通過逐步回歸分析,構(gòu)建二棱大麥種質(zhì)多性狀綜合評價的最優(yōu)回歸方程:Y=-7.660-0.002X1+0.017X2+0.072X3-0.016X4+0.169X5+0.041X6+0.210X7+0.073X8。式中Y代表綜合得分M值,X1、X2、X3、X4、X5、X6、X7、X8分別代表株高、穗下節(jié)間長、穗長、主穗粒數(shù)、單株穗數(shù)、千粒重、單株粒重、單株生物量,其中單株穗數(shù)、單株粒重的偏回歸系數(shù)較大。方程的決定系數(shù)R2為0.998,說明擬合很好,可用于二棱大麥種質(zhì)多性狀的綜合評價。

    3 討 論

    3.1 二棱大麥種質(zhì)主要性狀的表現(xiàn)及其遺傳多樣性

    二棱大麥性狀表型是大麥育種家選育高產(chǎn)品種的依據(jù),也是育種家對產(chǎn)量性狀遺傳機理進行合理解析的重要依據(jù),其遺傳多樣性為二棱大麥基因資源發(fā)掘提供必要的信息[29]。本研究二棱大麥的變異系數(shù)在不同表型性狀及試點間表現(xiàn)不同,穗長、單株穗數(shù)、單株粒重、單株生物量4性狀變異系數(shù)較大,易受環(huán)境影響,而株高、穗下節(jié)間長、主穗粒數(shù)、千粒重4性狀變異系數(shù)較小,說明二棱大麥籽粒性狀(穗粒數(shù)、千粒重)及株型性狀(株高、穗下節(jié)間長)相對穩(wěn)定,是育種的基礎,而育種的突破在于穗長和單株穗數(shù),這與張英虎等[3]、沈會權等[1]的研究結(jié)果較一致。2016年揚州大學新農(nóng)場試點的各性狀的變異系數(shù)普遍高于2015年和2016年的鹽城試點,但各性狀變化趨勢基本相同。

    遺傳多樣性指數(shù)已被廣泛應用于各類作物的多樣性評價中,尤其是核心種質(zhì)的多樣性評價[25-27]。二棱大麥種質(zhì)各性狀的遺傳多樣性指數(shù)的依次為主穗粒數(shù)(2.07)>穗下節(jié)間長(2.06)>單株粒重(2.03)>單株穗數(shù)(2.01)>株高(2.01)>千粒重(1.98)>穗長(1.89)>單株生物量(1.81),本研究二棱大麥各表型性狀的多樣性指數(shù)略高于李守明等[30]的研究結(jié)果,這與本研究選用參試大麥種質(zhì)來源廣泛有關。

    3.2 二棱大麥種質(zhì)不同評價方法的比較

    馬氏距離類平均法聚類可消除性狀間相關性對聚類結(jié)果的影響[31]。本研究將97份二棱大麥種質(zhì)資源分為8類,第Ⅰ類和第Ⅲ類的種質(zhì)較多,占種質(zhì)總數(shù)的一半以上,第Ⅰ類種質(zhì)為矮稈多粒型,第Ⅲ類種質(zhì)為多粒型,第Ⅳ類種質(zhì)為矮稈多穗豐產(chǎn)性型,第Ⅴ類種質(zhì)為多粒長穗型,第Ⅵ類為大粒型,第Ⅶ類種質(zhì)長穗、多穗、豐產(chǎn)性型,第Ⅷ類種質(zhì)為多粒大粒型,但Ⅵ、Ⅶ、Ⅷ類種質(zhì)高稈性狀需要改良,第Ⅱ類種質(zhì)為春大麥型,不適應江蘇生態(tài)條件,綜合表現(xiàn)較差,但從豐富基因來看,可用于個別性狀的改良親本。表型性狀雖然受環(huán)境影響,但性狀的表型鑒定比較直觀;分子標記鑒定結(jié)果比較可靠,但同一性狀往往由多個基因控制,基因又有主效和微效之分,很難利用少數(shù)幾個標記準確鑒定某一性狀的表型。李守明等[30]將大麥種質(zhì)性狀表型、醇溶蛋白和SSR標記聚類結(jié)果比較分析,認為三種聚類結(jié)果的對應關系并不明顯。因此,采取多環(huán)境、多重復鑒定大麥性狀表型,可以適當減少環(huán)境對性狀表型的影響程度,基于表型性狀聚類仍具有其實用性,必要時可以分子標記鑒定作為補充。

    主成分二維排序圖是利用各主成分得分繪制的二維散點圖,以第一主成分為橫坐標,其他若干主成分分別作縱坐標,根據(jù)各主成分的特性篩選相應特性種質(zhì)。本研究以二棱大麥種質(zhì)的大粒、高產(chǎn)為基礎,對兼具穗數(shù)與穗粒數(shù)并重及矮稈特性的3類種質(zhì)進行分析,結(jié)果表明,3類特異性的種質(zhì)在3個二維排序圖中沒有發(fā)生重疊,表明參試材料特點突出,二棱大麥育種仍具有較大的上升空間。二維排序分析對區(qū)分不同種質(zhì)農(nóng)藝、產(chǎn)量性狀優(yōu)劣以及評價特異種質(zhì)方面較為清晰直觀[23],其對種質(zhì)評價的綜合性尚不夠。

    種質(zhì)的綜合評價值(M值)可為育種家提供直觀、便捷、量化的參考指標。M值排列前十位的二棱大麥種質(zhì)中,除了一份地方種質(zhì)外,其余均為江浙滬周圍地區(qū)選育的品種品系。

    通過逐步回歸構(gòu)建了二棱大麥種質(zhì)多性狀綜合評價最優(yōu)方程:Y=-7.660-0.002X1+0.017X2+0.072X3-0.016X4+0.169X5+0.041X6+0.210X7+0.073X8,8個主要性狀均顯著影響二棱大麥種質(zhì)的表型性狀綜合值(M值),均應作為二棱大麥種質(zhì)綜合評價性狀指標。

    參考文獻:

    [1] 沈會權,張英虎,欒海業(yè),等.江蘇二棱大麥農(nóng)藝、籽粒和品質(zhì)性狀的特征及其相關性分析[J].麥類作物學報,2016,36(3):379.

    SHEN H Q,ZHANG Y H,LUAN H Y,etal.Characteristics and correlation of the agronomic,grain and quality trait of the two-row barley in Jiangsu [J].JournalofTriticeaeCrops,2016,36(3):379.

    [2] 李洪濤,王 軍,張燦宏,等.大麥株高及其構(gòu)成因素的雜種優(yōu)勢和配合力分析[J].浙江農(nóng)業(yè)學報,2015,27(2):141.

    LI H T,WANG J,ZHANG C H,etal.Heterosis and combining ability of plant height and its components of barley [J].ActaAgriculturaeZhejiangensis,2015,27(2):141.

    [3] 張英虎,陳 和,陳 健,等.2007-2014年江蘇沿海地區(qū)大麥品種(系)產(chǎn)量與產(chǎn)量構(gòu)成分析[J].浙江農(nóng)業(yè)學報,2015,27(9):1510.

    ZHANG Y H,CHEN H,CHEN J,etal.Analysis of yield and yield components of barley varieties(lines) in coastal areas in Jiangsu during 2007 to 2014 [J].ActaAgriculturaeZhejiangensis,2015,27(9):1510.

    [4] 方桂英,杜向紅,張保軍,等.約旦野生大麥種質(zhì)資源形態(tài)與產(chǎn)量性狀分析[J].西北農(nóng)業(yè)學報,2010,19(6):78.

    FANG G Y,DU X H,ZHANG B J,etal.Analysis of agronomic traits ofHordeumspontaneumC.Koch germplasm from Jordan [J].ActaAgriculturaeBoreali-occidentalisSinica,2010,19(6):78.

    [5] 孫立軍.中國大麥遺傳資源和優(yōu)異種質(zhì)[M].北京:中國農(nóng)業(yè)科技出版社,2001.

    SUN L J.Chinese barley genetic resources and excellent germplasm [M].Beijing:China Agricultural Science and Technology Press,2001.

    [6] WARD D J.Some evolutionary aspects of certain morphologic characters in a world collection of barley [J].TechnicalBulletins,1962:112.

    [7] TOLBERT D M,QUALSET C O,JAIN S K,etal.A diversity analysis of a world collection of barley [J].CropScience,1979,19(6):789.

    [8] BACKES G,HATZ B,JAHOOR A,etal.RFLP diversity within and between major groups of barley in Europe [J].PlantBreeding,2003,122(4):291.

    [9] 侯永翠,顏澤洪,魏育明,等.利用RAPD標記分析大麥種質(zhì)資源的遺傳多樣性[J].植物遺傳資源學報,2005,6(2):145.

    HOU Y C,YAN Z H,WEI Y M,etal.Using RAPD Markers to analysis genetic diversity in barley [J].JournalofPlantGeneticResources,2005,6(2):145.

    [10] 張 鏑,丁 毅.基于AFLP標記的中國西藏近緣野生大麥遺傳多樣性分析[J].遺傳,2007,29(6):725.

    ZHANG D,DING Y.Genetic diversity of wild close relatives of barley in Tibet of China revealed by AFLP [J].Hereditas,2007,29(6):725.

    [11] 張赤紅,張 京.大麥品種資源遺傳多樣性的SSR標記評價[J].麥類作物學報,2008,28(2):214.

    ZHANG C H,ZHANG J.Genetic diversity assessment of barley germplasm resources using SSR markers [J].JournalofTriticeaeCrops,2008,28(2):214.

    [12] WANG J,WANG S H,LAI Y,etal.Genetic diversity and population structure analysis by using SSR markers in barley [J].ActaAgriculturaeNucleataeSinica,2014,28(2):177.

    [13] BUNDOCK P C,HENRY R J.Single nucleotide polymorphism,haplotype diversity and recombination in theIsagene of barley [J].Tag.theoretical&AppliedGenetics.theoretischeUndAngewandteGenetik,2004,109(3):543.

    [14] CRONIN J K,NEVO E.Adaptive climatic molecular evolution in wild barley at theIsadefense locus [J].ProceedingsoftheNationalAcademyofSciencesoftheUnitedStatesofAmerica,2007,104(8):2773.

    [15] ZHANG Q,MAROOF M A,KLEINHOFS A.Comparative diversity analysis of RFLPs and isozymes within and among populations ofHordeumvulgaressp.spontaneum [J].Genetics,1993,134(3):909.

    [16] 陳曉靜,陳 和,陳 健,等.300份國外啤酒大麥引種鑒定與分析[J].江西農(nóng)業(yè)學報,2007,19(12):21.

    CHEN X J,CHEN H,CHEN J,etal.Identification and analysis of 300 foreign malting barleys [J].ActaAgricultureJiangxi,2007,19(12):21.

    [17] 張小燕,張躍進,龐有強.日本大麥種質(zhì)資源鑒定評價[J].西北農(nóng)業(yè)學報,2009,18(6):122.

    ZHANG X Y,ZHANG Y J,PANG Y L.Evaluation on barley germplasm from Japan [J].ActaAgriculturaeBoreali-occidentalisSinica,2009,18(6):122.

    [18] 周 偉.引進二棱大麥種質(zhì)資源的主成分分析與聚類分析[J].河北農(nóng)業(yè)科學,2010,14(7):56.

    ZHOU W.Principal Component analysis and cluster analysis of introducing two-rowed barley germplasm [J].JournalofHebeiAgriculturalSciences,2010,14(7):56.

    [19] 楊樹明,普曉英,張 京,等.不同地區(qū)啤酒大麥品種農(nóng)藝性狀鑒定與分類研究[J].植物遺傳資源學報,2011,12(1):37.

    YANG S M,PU X Y,ZHANG J.Agronomic characteristics and classification of malting barley from different regions [J].JournalofPlantGeneticResources,2011,12(1):37.

    [20] 陶愛芬,祁建民,林培青,等.紅麻優(yōu)異種質(zhì)產(chǎn)量和品質(zhì)性狀主成分聚類分析與綜合評價[J].中國農(nóng)業(yè)科學,2008,41(9):2859.

    TAO A F,QI J M,LIN P Q,etal.Cluster analysis and evaluation of elite kanaf germplasm based on principal components [J].ScientiaAgriculturaSinica,2008,41(9):2859.

    [21] 王利民,張建平,米 君,等.國外引進油用亞麻品種資源農(nóng)藝性狀分析與評價[J].中國油料作物學報,2011,33(4):356.

    WANG L M,ZHANG J P,MI J,etal.Agronomic traits analysis and evaluation of introduced flax germplasm based on principal components and cluster analysis [J].ChineseJournalofOilCropSciences,2011,33(4):356.

    [22] 陳書霞,周 靜,申曉青,等.大蒜種質(zhì)產(chǎn)量和品質(zhì)性狀主成分聚類分析與綜合評價[J].植物遺傳資源學報,2012,13(3):429.

    CHEN S X,ZHOU J,SHEN X Q,etal.Cluster analysis and evaluation of garlic(AlliumsativumL.) germplasm based on principal components [J].JournalofPlantGeneticResources,2012,13(3):429.

    [23] 要燕杰,高 翔,吳 丹,等.小麥農(nóng)藝性狀與品質(zhì)特性的多元分析與評價[J].植物遺傳資源學報,2014,15(1):38.

    YAO Y J,GAO X,WU D,etal.Multivariate analysis and evaluation of agronomic and quality traits based on principal components in wheat [J].JournalofPlantGeneticResources,2014,15(1):38.

    [24] 馬 麒,宿俊吉,寧新柱,等.新疆海島棉育種目標性狀多元分析與種質(zhì)資源評價[J].西南農(nóng)業(yè)學報,2016,29(7):1530.

    MA Q,SU J J,NING X Z,etal.Multivariate analysis on breeding targeted traits of sea island cotton(G.barbadense) bred in Xinjiang and evaluation of germplasm resources [J].SouthwestChinaJournalofAgriculturalSciences,2016,29(7):1530.

    [25] 胡標林,萬 勇,李 霞,等.水稻核心種質(zhì)表型性狀遺傳多樣性分析及綜合評價[J].作物學報,2012,38(5):829.

    HU B L,WANG Y,LI X,etal.Analysis on genetic diversity of phenotypic traits in rice(Oryzasativa) core collection and its comprehensive assessment [J].ActaAgronomicaSinica,2012,38(5):829.

    [26] 代攀虹,孫君靈,何守樸,等.陸地棉核心種質(zhì)表型性狀遺傳多樣性分析及綜合評價[J].中國農(nóng)業(yè)科學,2016,49(19):3694.

    DAI P H,SUN L J,HE S P,etal.Comprehensive evaluation and genetic diversity analysis of phenotypic traits of core collection in upland cotton [J].ScientiaAgriculturaSinica,2016,49(19):3694.

    [27] 王海崗,賈冠清,智 慧,等.谷子核心種質(zhì)表型遺傳多樣性分析及綜合評價[J].作物學報,2016,42(1):19.

    WANG H G,JIA G Q,ZHI H,etal.Phenotypic diversity evaluations of foxtail millet core collections [J].ActaAgronomicaSinica,2016,42(1):19.

    [28] 董玉琛,曹永生,張學勇,等.中國普通小麥初選核心種質(zhì)的產(chǎn)生[J].植物遺傳資源學報,2003,4(1):1.

    DONG Y S,CAO Y S,ZHANG X Y,etal.Establishment of candidate core collections in Chinese common wheat germplasm [J].JournalofPlantGeneticResources,2003,4(1):1.

    [29] CAMPBELL D R.Using phenotypic manipulations to study multivariate selection of floral trait associations [J].AnnalsofBotany,2009,103(9):1557.

    [30] 李守明,梁 維,魏凌基,等.新疆大麥種質(zhì)資源農(nóng)藝性狀和醇溶蛋白遺傳多樣性分析[J].新疆農(nóng)業(yè)科學,2009,46(2):269.

    LI S M,LIANG W,WEI L J,etal.Agronomic trait of germplasm resources and genetic diversity of hordein in Xinjiang barley [J].XinjiangAgriculturalSciences,2009,46(2):269.

    [31] 劉遵春,張春雨,張艷敏,等.利用數(shù)量性狀構(gòu)建新疆野蘋果核心種質(zhì)的方法[J].中國農(nóng)業(yè)科學,2010,43(2):358.

    LIU Z C,ZHANG C Y,ZHANG Y M,etal.Study on method of constructing core collection ofMalussieversiibased on quantitative traits [J].ScientiaAgriculturaSinica,2010,43(2):358.

    猜你喜歡
    穗長粒重穗數(shù)
    不同高低畦種植模式對冬小麥干物質(zhì)積累和產(chǎn)量的影響
    播期和密度對揚麥25產(chǎn)量及產(chǎn)量構(gòu)成因子的影響
    基于GBS測序和連鎖分析的藜麥單株粒重QTL定位
    種子(2023年9期)2023-11-22 13:10:56
    四倍體小麥株高和穗長性狀的QTL定位及其遺傳效應分析
    離體穗培養(yǎng)條件下C、N供給對小麥穗粒數(shù)、粒重及蛋白質(zhì)含量的影響
    玉米自交系京92改良后代單穗粒重的雜種優(yōu)勢研究
    雜交晚粳稻通優(yōu)粳1號產(chǎn)量及構(gòu)成因子分析
    高產(chǎn)小麥構(gòu)成要素的解析與掌握
    水稻株高、穗長和每穗穎花數(shù)的遺傳研究
    水稻穗長和有效穗數(shù)的QTL定位分析
    欧美成人午夜精品| 日本色播在线视频| 亚洲精品一区蜜桃| 国产av国产精品国产| 菩萨蛮人人尽说江南好唐韦庄| 久久免费观看电影| 不卡视频在线观看欧美| 极品少妇高潮喷水抽搐| 国产人伦9x9x在线观看| 在线精品无人区一区二区三| 中文字幕色久视频| 欧美精品高潮呻吟av久久| 亚洲av男天堂| 欧美 日韩 精品 国产| 最新在线观看一区二区三区 | 嫩草影视91久久| 精品一区二区免费观看| 欧美97在线视频| 国产精品一区二区在线不卡| 黄色视频不卡| 黄网站色视频无遮挡免费观看| 丝袜美足系列| 国产精品欧美亚洲77777| 精品一区二区三区四区五区乱码 | 国产在线一区二区三区精| 亚洲欧美一区二区三区久久| 777米奇影视久久| 久久久久精品久久久久真实原创| 日韩人妻精品一区2区三区| 人妻一区二区av| 国产亚洲av片在线观看秒播厂| 一二三四在线观看免费中文在| 人妻人人澡人人爽人人| 成人18禁高潮啪啪吃奶动态图| 中文字幕色久视频| 99香蕉大伊视频| 最新的欧美精品一区二区| 黄频高清免费视频| av有码第一页| 国产毛片在线视频| 久热这里只有精品99| 老熟女久久久| 激情视频va一区二区三区| 久久久久网色| 亚洲国产欧美网| 亚洲精品久久午夜乱码| 大码成人一级视频| 国产精品无大码| 精品国产乱码久久久久久男人| 欧美黑人欧美精品刺激| 成人午夜精彩视频在线观看| 女人被躁到高潮嗷嗷叫费观| 水蜜桃什么品种好| 日本欧美国产在线视频| 国产成人a∨麻豆精品| 欧美日韩视频高清一区二区三区二| 亚洲综合色网址| av福利片在线| 国产精品国产三级专区第一集| 美女脱内裤让男人舔精品视频| 欧美人与善性xxx| 午夜影院在线不卡| 久久婷婷青草| 老司机靠b影院| 国产片特级美女逼逼视频| 国产欧美日韩综合在线一区二区| 国产av一区二区精品久久| 日韩一区二区视频免费看| 成人午夜精彩视频在线观看| 国产野战对白在线观看| 九九爱精品视频在线观看| 老司机靠b影院| 777久久人妻少妇嫩草av网站| 99久国产av精品国产电影| 亚洲成人手机| 日日摸夜夜添夜夜爱| 欧美久久黑人一区二区| 我的亚洲天堂| 欧美日韩亚洲国产一区二区在线观看 | 国产日韩欧美视频二区| www日本在线高清视频| 啦啦啦在线免费观看视频4| 精品第一国产精品| 国产男女内射视频| 亚洲成国产人片在线观看| 国产精品秋霞免费鲁丝片| 欧美人与性动交α欧美精品济南到| 亚洲精品aⅴ在线观看| 精品国产一区二区三区久久久樱花| 久久久久久久久久久免费av| 宅男免费午夜| 亚洲国产毛片av蜜桃av| 男女午夜视频在线观看| 嫩草影视91久久| 亚洲精品国产区一区二| 妹子高潮喷水视频| 热99国产精品久久久久久7| 国产99久久九九免费精品| www.熟女人妻精品国产| 国产精品免费大片| 熟女少妇亚洲综合色aaa.| bbb黄色大片| 91aial.com中文字幕在线观看| 国产精品 国内视频| 母亲3免费完整高清在线观看| 肉色欧美久久久久久久蜜桃| 香蕉丝袜av| 久久亚洲国产成人精品v| 欧美在线黄色| 国产精品香港三级国产av潘金莲 | 咕卡用的链子| 丁香六月欧美| 成人影院久久| 欧美另类一区| 2018国产大陆天天弄谢| 一区二区av电影网| 久久精品久久久久久噜噜老黄| 波多野结衣一区麻豆| 如日韩欧美国产精品一区二区三区| 又大又爽又粗| 亚洲精品日韩在线中文字幕| 亚洲一码二码三码区别大吗| 伊人亚洲综合成人网| 国产精品 国内视频| 精品免费久久久久久久清纯 | 一级毛片 在线播放| 99精品久久久久人妻精品| 人人妻人人澡人人看| 中文字幕人妻熟女乱码| 青草久久国产| 日韩av在线免费看完整版不卡| 永久免费av网站大全| 肉色欧美久久久久久久蜜桃| 亚洲av日韩精品久久久久久密 | 可以免费在线观看a视频的电影网站 | 午夜久久久在线观看| av网站在线播放免费| 国产精品人妻久久久影院| 亚洲av福利一区| 大码成人一级视频| 一边亲一边摸免费视频| 国产男人的电影天堂91| 精品一区二区三区四区五区乱码 | 亚洲专区中文字幕在线 | 最新的欧美精品一区二区| 国产精品 国内视频| 亚洲在久久综合| 午夜免费男女啪啪视频观看| 亚洲,欧美精品.| 欧美人与善性xxx| 亚洲男人天堂网一区| 黄色一级大片看看| 成年人免费黄色播放视频| 老汉色av国产亚洲站长工具| 久久久久精品人妻al黑| 十八禁网站网址无遮挡| 国产在线免费精品| 夫妻性生交免费视频一级片| 中文字幕另类日韩欧美亚洲嫩草| 母亲3免费完整高清在线观看| 汤姆久久久久久久影院中文字幕| 波多野结衣一区麻豆| 久久久久精品久久久久真实原创| 成人三级做爰电影| 91精品伊人久久大香线蕉| 最近的中文字幕免费完整| 99热全是精品| 丝瓜视频免费看黄片| 在线观看一区二区三区激情| 国产午夜精品一二区理论片| 大香蕉久久网| 大话2 男鬼变身卡| 午夜日本视频在线| 国产极品粉嫩免费观看在线| 亚洲精品乱久久久久久| 成年美女黄网站色视频大全免费| 亚洲国产精品一区二区三区在线| 男女免费视频国产| 91国产中文字幕| 极品人妻少妇av视频| 大码成人一级视频| 欧美精品av麻豆av| 久久精品久久精品一区二区三区| 亚洲国产最新在线播放| 丝袜脚勾引网站| 超碰97精品在线观看| av在线app专区| 热99久久久久精品小说推荐| 亚洲中文av在线| 大话2 男鬼变身卡| 香蕉国产在线看| 成人国产麻豆网| 国产精品欧美亚洲77777| 精品少妇久久久久久888优播| 日韩欧美精品免费久久| 97人妻天天添夜夜摸| 如何舔出高潮| 国产野战对白在线观看| 19禁男女啪啪无遮挡网站| 亚洲欧美精品自产自拍| 日韩av不卡免费在线播放| 亚洲天堂av无毛| 性少妇av在线| 亚洲国产中文字幕在线视频| 国产精品一区二区在线不卡| 好男人视频免费观看在线| 水蜜桃什么品种好| 如何舔出高潮| 国产免费又黄又爽又色| 99国产精品免费福利视频| 99国产精品免费福利视频| 成人亚洲欧美一区二区av| 永久免费av网站大全| 香蕉国产在线看| 十分钟在线观看高清视频www| 国产人伦9x9x在线观看| 丝袜美足系列| 男人爽女人下面视频在线观看| 国产亚洲av高清不卡| 欧美日韩亚洲高清精品| 巨乳人妻的诱惑在线观看| 精品人妻在线不人妻| 在线天堂中文资源库| 大话2 男鬼变身卡| 婷婷色综合大香蕉| 熟女少妇亚洲综合色aaa.| 成人漫画全彩无遮挡| 一本大道久久a久久精品| 青春草国产在线视频| 黑人巨大精品欧美一区二区蜜桃| 一级爰片在线观看| 国产一卡二卡三卡精品 | 国产精品久久久久久人妻精品电影 | 夫妻性生交免费视频一级片| 男人操女人黄网站| 国产精品三级大全| 欧美日韩av久久| 香蕉国产在线看| 国产精品无大码| 日日撸夜夜添| 国产精品一区二区在线不卡| 男男h啪啪无遮挡| 97人妻天天添夜夜摸| √禁漫天堂资源中文www| 搡老乐熟女国产| 一本久久精品| 亚洲精品国产区一区二| 搡老岳熟女国产| 伦理电影大哥的女人| 1024视频免费在线观看| 桃花免费在线播放| 我的亚洲天堂| a级毛片黄视频| 乱人伦中国视频| 午夜福利,免费看| 国产女主播在线喷水免费视频网站| 哪个播放器可以免费观看大片| 在线观看三级黄色| 丁香六月欧美| 亚洲精品成人av观看孕妇| av线在线观看网站| 国产精品香港三级国产av潘金莲 | 亚洲第一青青草原| 人成视频在线观看免费观看| 国产精品欧美亚洲77777| 曰老女人黄片| www.自偷自拍.com| 亚洲成人国产一区在线观看 | 午夜福利免费观看在线| 久久免费观看电影| 国产精品亚洲av一区麻豆 | 成人午夜精彩视频在线观看| 人体艺术视频欧美日本| 日韩伦理黄色片| 欧美少妇被猛烈插入视频| 夜夜骑夜夜射夜夜干| 在线观看国产h片| 久久ye,这里只有精品| 欧美日韩国产mv在线观看视频| 韩国高清视频一区二区三区| 青春草视频在线免费观看| 女人久久www免费人成看片| 91国产中文字幕| 美女扒开内裤让男人捅视频| 亚洲婷婷狠狠爱综合网| 天堂俺去俺来也www色官网| 天天影视国产精品| 国产男女超爽视频在线观看| 制服诱惑二区| 久久久久精品国产欧美久久久 | 在线观看www视频免费| 80岁老熟妇乱子伦牲交| 欧美xxⅹ黑人| 免费不卡黄色视频| 观看美女的网站| 日韩一区二区三区影片| 午夜福利一区二区在线看| 亚洲欧美精品自产自拍| 视频区图区小说| 亚洲精品av麻豆狂野| 国产亚洲一区二区精品| 国产成人午夜福利电影在线观看| 亚洲一码二码三码区别大吗| 18禁裸乳无遮挡动漫免费视频| 亚洲成人手机| 中文字幕亚洲精品专区| 国产不卡av网站在线观看| 免费女性裸体啪啪无遮挡网站| 在线观看免费午夜福利视频| 成人漫画全彩无遮挡| 午夜激情av网站| 巨乳人妻的诱惑在线观看| 大话2 男鬼变身卡| 日韩大片免费观看网站| 咕卡用的链子| 国产乱人偷精品视频| 捣出白浆h1v1| 国产精品偷伦视频观看了| av一本久久久久| 一二三四在线观看免费中文在| 成年动漫av网址| 一区二区三区激情视频| 欧美97在线视频| 免费女性裸体啪啪无遮挡网站| 亚洲国产av影院在线观看| 久久韩国三级中文字幕| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产毛片av蜜桃av| 中文字幕另类日韩欧美亚洲嫩草| 日日撸夜夜添| 亚洲国产最新在线播放| 男女午夜视频在线观看| 成人亚洲精品一区在线观看| 香蕉国产在线看| 日本爱情动作片www.在线观看| 天堂俺去俺来也www色官网| 久久天堂一区二区三区四区| 久久天躁狠狠躁夜夜2o2o | 女人精品久久久久毛片| 国产在视频线精品| 欧美av亚洲av综合av国产av | 91精品伊人久久大香线蕉| 免费黄网站久久成人精品| 90打野战视频偷拍视频| 久久这里只有精品19| 日日爽夜夜爽网站| 晚上一个人看的免费电影| 亚洲中文av在线| 国产精品人妻久久久影院| av不卡在线播放| 欧美精品高潮呻吟av久久| 久久99一区二区三区| 中文字幕色久视频| 午夜福利视频在线观看免费| 黄色一级大片看看| 亚洲欧美成人综合另类久久久| 啦啦啦视频在线资源免费观看| 亚洲av综合色区一区| 一本一本久久a久久精品综合妖精| 大香蕉久久成人网| 黄色怎么调成土黄色| 飞空精品影院首页| 岛国毛片在线播放| 在线观看免费午夜福利视频| 久久久久久久精品精品| 热re99久久精品国产66热6| 国产精品香港三级国产av潘金莲 | 黑丝袜美女国产一区| 亚洲免费av在线视频| 中文乱码字字幕精品一区二区三区| 国产精品香港三级国产av潘金莲 | 欧美xxⅹ黑人| 一边亲一边摸免费视频| 久久人人爽人人片av| 亚洲中文av在线| 亚洲综合色网址| 激情视频va一区二区三区| 亚洲中文av在线| 亚洲欧洲国产日韩| 我要看黄色一级片免费的| 免费黄频网站在线观看国产| 亚洲第一区二区三区不卡| 亚洲精品乱久久久久久| 国产精品久久久久成人av| 国产淫语在线视频| 日韩一本色道免费dvd| 黑人欧美特级aaaaaa片| www.自偷自拍.com| 大陆偷拍与自拍| 亚洲精品乱久久久久久| 丝袜脚勾引网站| 亚洲精品一二三| 91aial.com中文字幕在线观看| 天堂8中文在线网| 亚洲精品国产区一区二| 人人妻人人添人人爽欧美一区卜| 中文字幕制服av| 久久女婷五月综合色啪小说| 日本vs欧美在线观看视频| 中文字幕精品免费在线观看视频| 无遮挡黄片免费观看| 精品福利永久在线观看| 啦啦啦在线观看免费高清www| 一二三四在线观看免费中文在| 欧美激情 高清一区二区三区| 久久精品亚洲av国产电影网| 啦啦啦在线免费观看视频4| 亚洲欧美色中文字幕在线| 亚洲精品自拍成人| 久久97久久精品| 在线观看人妻少妇| 一级片免费观看大全| 中国国产av一级| 国产爽快片一区二区三区| 国产av精品麻豆| 多毛熟女@视频| 欧美成人精品欧美一级黄| 国产av码专区亚洲av| 国产激情久久老熟女| 亚洲成人一二三区av| 午夜福利在线免费观看网站| 国产亚洲欧美精品永久| 亚洲欧美激情在线| 美女高潮到喷水免费观看| 国产精品人妻久久久影院| 一二三四在线观看免费中文在| 在线精品无人区一区二区三| 国产一区二区在线观看av| 一边摸一边抽搐一进一出视频| 黄片小视频在线播放| 日韩精品有码人妻一区| 成年人午夜在线观看视频| 久久毛片免费看一区二区三区| 99精品久久久久人妻精品| 欧美人与性动交α欧美软件| 久久天堂一区二区三区四区| 亚洲欧洲国产日韩| 亚洲一区中文字幕在线| 国产在线视频一区二区| 亚洲成国产人片在线观看| 97精品久久久久久久久久精品| 久久鲁丝午夜福利片| 成人国产av品久久久| 夫妻午夜视频| 在线亚洲精品国产二区图片欧美| 美女脱内裤让男人舔精品视频| 男女高潮啪啪啪动态图| videos熟女内射| 亚洲精品久久久久久婷婷小说| 毛片一级片免费看久久久久| 考比视频在线观看| 国产国语露脸激情在线看| 9热在线视频观看99| 美女高潮到喷水免费观看| xxx大片免费视频| 精品国产超薄肉色丝袜足j| av.在线天堂| 日本av手机在线免费观看| 国产xxxxx性猛交| 欧美激情高清一区二区三区 | 男人添女人高潮全过程视频| 国产一区二区三区av在线| www日本在线高清视频| 欧美日本中文国产一区发布| 少妇的丰满在线观看| 悠悠久久av| av网站免费在线观看视频| 卡戴珊不雅视频在线播放| 久久毛片免费看一区二区三区| 操美女的视频在线观看| 日韩一区二区三区影片| 你懂的网址亚洲精品在线观看| 国产精品女同一区二区软件| 午夜影院在线不卡| 99热国产这里只有精品6| 黄色 视频免费看| 国产xxxxx性猛交| 欧美日韩福利视频一区二区| 国产精品免费视频内射| 免费高清在线观看日韩| 一区二区三区四区激情视频| 日韩电影二区| 国产精品秋霞免费鲁丝片| 秋霞在线观看毛片| 亚洲七黄色美女视频| 美女福利国产在线| 深夜精品福利| √禁漫天堂资源中文www| 老熟女久久久| av免费观看日本| 丝袜脚勾引网站| 国产av码专区亚洲av| 国精品久久久久久国模美| kizo精华| 一边亲一边摸免费视频| 人人妻,人人澡人人爽秒播 | 国产一卡二卡三卡精品 | 一本久久精品| 久久精品aⅴ一区二区三区四区| 午夜激情久久久久久久| 亚洲精品国产区一区二| 色精品久久人妻99蜜桃| av不卡在线播放| 国产免费现黄频在线看| 精品一区在线观看国产| 欧美av亚洲av综合av国产av | 国产无遮挡羞羞视频在线观看| 永久免费av网站大全| 欧美乱码精品一区二区三区| 在线天堂中文资源库| 亚洲av电影在线观看一区二区三区| 国产精品香港三级国产av潘金莲 | 日韩不卡一区二区三区视频在线| 亚洲精品成人av观看孕妇| 男女免费视频国产| 精品卡一卡二卡四卡免费| 99re6热这里在线精品视频| 中文字幕最新亚洲高清| 美女中出高潮动态图| 日韩一卡2卡3卡4卡2021年| 国产精品久久久久久久久免| 天天躁夜夜躁狠狠久久av| 亚洲欧美清纯卡通| 久久精品人人爽人人爽视色| 亚洲 欧美一区二区三区| 免费在线观看完整版高清| www.精华液| 午夜影院在线不卡| 七月丁香在线播放| 亚洲精品国产av蜜桃| 国产黄频视频在线观看| 在线 av 中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 在线观看免费高清a一片| 亚洲在久久综合| 观看美女的网站| 成年动漫av网址| 最近最新中文字幕免费大全7| 狠狠精品人妻久久久久久综合| 久久精品亚洲av国产电影网| 最近2019中文字幕mv第一页| 国产伦人伦偷精品视频| 一区二区三区乱码不卡18| 久久精品人人爽人人爽视色| 久久99热这里只频精品6学生| 久久久欧美国产精品| 国产激情久久老熟女| 欧美精品高潮呻吟av久久| 日韩大片免费观看网站| 亚洲少妇的诱惑av| 亚洲美女视频黄频| 纯流量卡能插随身wifi吗| 亚洲人成电影观看| 日韩 亚洲 欧美在线| 亚洲国产欧美一区二区综合| www日本在线高清视频| 欧美xxⅹ黑人| 无限看片的www在线观看| 久久精品人人爽人人爽视色| 亚洲,一卡二卡三卡| 不卡视频在线观看欧美| 日本av手机在线免费观看| 欧美黑人欧美精品刺激| 多毛熟女@视频| 少妇被粗大的猛进出69影院| 一级片'在线观看视频| 男人添女人高潮全过程视频| 亚洲精品乱久久久久久| 在线观看人妻少妇| 国产成人免费观看mmmm| 精品亚洲成a人片在线观看| 性少妇av在线| 国产日韩欧美亚洲二区| 99久久精品国产亚洲精品| 涩涩av久久男人的天堂| 十分钟在线观看高清视频www| 国产亚洲av片在线观看秒播厂| 一级毛片 在线播放| 国产日韩欧美亚洲二区| 亚洲免费av在线视频| 丁香六月天网| 中文字幕色久视频| 国产又爽黄色视频| 亚洲国产中文字幕在线视频| 久热这里只有精品99| 熟女少妇亚洲综合色aaa.| 秋霞伦理黄片| 国产亚洲最大av| 免费在线观看完整版高清| 校园人妻丝袜中文字幕| 亚洲国产日韩一区二区| 午夜福利视频精品| 亚洲av日韩精品久久久久久密 | 国产黄频视频在线观看| 亚洲精品国产区一区二| 美女主播在线视频| 18禁观看日本| 国产麻豆69| 一级爰片在线观看| 亚洲国产欧美在线一区| av视频免费观看在线观看| 久久午夜综合久久蜜桃| 人妻一区二区av| 国产成人av激情在线播放| 久久毛片免费看一区二区三区| 各种免费的搞黄视频| 男人舔女人的私密视频| 亚洲激情五月婷婷啪啪| 高清欧美精品videossex| 亚洲精品美女久久久久99蜜臀 | a级毛片在线看网站| 人妻人人澡人人爽人人| 大香蕉久久网| 日日爽夜夜爽网站| 国产熟女欧美一区二区|