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

    基于形態(tài)學(xué)性狀和SSR標(biāo)記的花生品種遺傳多樣性分析和特異性鑒定

    2019-12-25 10:58:36徐振江饒得花李少雄劉海燕陳小平梁炫強(qiáng)洪彥彬
    作物學(xué)報 2019年1期
    關(guān)鍵詞:類群形態(tài)學(xué)花生

    劉 洪 徐振江 饒得花 魯 清 李少雄 劉海燕 陳小平 梁炫強(qiáng) 洪彥彬,*

    ?

    基于形態(tài)學(xué)性狀和SSR標(biāo)記的花生品種遺傳多樣性分析和特異性鑒定

    劉 洪1,**徐振江1,**饒得花1魯 清2,3李少雄2,3劉海燕2,3陳小平2,3梁炫強(qiáng)2,3洪彥彬2,3,*

    1華南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院, 廣東廣州 510642;2廣東省農(nóng)業(yè)科學(xué)院作物研究所, 廣東廣州 510640;3廣東省農(nóng)作物遺傳改良重點實驗室, 廣東廣州 510640

    以101份南方花生區(qū)試品種為材料, 利用形態(tài)學(xué)性狀和SSR標(biāo)記進(jìn)行品種遺傳多樣性分析和特異性鑒定。結(jié)果表明, 29個形態(tài)學(xué)性狀中有7個無多樣性, 其余22個的多樣性指數(shù)為0.23~0.77, 平均為0.43。在相似系數(shù)為0.76處, 將供試品種劃分為七大類群, 同一育種單位的品種傾向于聚在一起。用40個SSR標(biāo)記共檢測出167個等位基因, 單個標(biāo)記檢測的等位基因數(shù)2~6個, 平均為4.18個。標(biāo)記的多態(tài)性信息量(PIC)差異較大, 最大為0.79, 最小為0.26, 平均為0.55。在相似系數(shù)為0.70處, 供試品種可被劃分為六大類群, 同一省份育成的品種多聚為一類。Mantel檢驗發(fā)現(xiàn)品種間的形態(tài)學(xué)性狀和SSR標(biāo)記的相似系數(shù)矩陣相關(guān)性弱(= 0.36), SSR標(biāo)記無法取代形態(tài)學(xué)性狀單獨用于花生品種特異性鑒定, 但兩者相結(jié)合能有效提高花生品種特異性鑒定的準(zhǔn)確性。

    花生; 形態(tài)學(xué); SSR; 遺傳多樣性; DUS

    花生栽培種(L.)屬落花生屬()花生區(qū)組(section), 異源四倍體(2= 4= 40, AABB), 是世界范圍內(nèi)廣泛栽培的重要油料和經(jīng)濟(jì)作物。因其適應(yīng)性廣、豐產(chǎn)性好、營養(yǎng)豐富、經(jīng)濟(jì)效益高、適于加工, 在促進(jìn)農(nóng)民增收和維護(hù)國家糧油安全中發(fā)揮重要作用。我國不但是花生種植、消費和出口大國, 同時也是花生育種強(qiáng)國。截止2017年7月, 共有409個花生品種申請植物新品種權(quán), 其中授權(quán)品種116個[1]。

    中國是國際植物新品種保護(hù)聯(lián)盟(UPOV)的成員國, 依據(jù)植物新品種保護(hù)條例[2]的規(guī)定, 植物新品種必需具備特異性(distinctness)、一致性(uniformity)和穩(wěn)定性(stability) (簡稱DUS), 才能授予品種權(quán)。特異性測試是DUS測試的核心, 而準(zhǔn)確篩選近似品種是特異性測試的基本要求。目前花生近似品種篩選仍以形態(tài)學(xué)鑒定為基礎(chǔ), 雖然該方法簡單、直接, 但容易受環(huán)境條件和人為因素的影響。此外, 由于過分強(qiáng)調(diào)花生品種的豐產(chǎn)性, 育種家通常選用高配合力親本進(jìn)行雜交育種, 導(dǎo)致育成品種親緣關(guān)系近, 后代遺傳變異停留在近交水平[3], 從而加大了表型篩選近似品種的難度。尤其近幾年品種侵權(quán)事件不斷增加, 同種多名和同名多種等不良現(xiàn)象的涌現(xiàn), 給花生新品種保護(hù)工作和品種管理帶來很大困難。

    與形態(tài)學(xué)標(biāo)記相比, 分子標(biāo)記能對各發(fā)育時期的個體、組織、器官甚至細(xì)胞做檢測, 既不受環(huán)境的影響, 也不受基因表達(dá)與否的限制, 數(shù)量豐富, 遺傳穩(wěn)定, 對生物體的影響表現(xiàn)為“中性”[4]。分子標(biāo)記的所有這些特性, 奠定其廣泛應(yīng)用前景, 并成為植物新品種測試審查和品種鑒別的重要技術(shù)手段[5-6]。早前Tommasini等[7]利用SSR標(biāo)記對油菜品種特異性的測試表明, SSR標(biāo)記可作為油菜DUS測試體系的有效補(bǔ)充。隨后, Arens等[8]在番茄上開發(fā)8個與抗性高度相關(guān)的分子標(biāo)記, 并用于番茄品種特異性鑒定。結(jié)果表明, 該套標(biāo)記具有滿足品種DUS測試要求的潛力, 可作為當(dāng)前DUS測試中抗性鑒定的補(bǔ)充方法。近期, Jones和Mackay[9]通過高密度SNP標(biāo)記分析大麥品種的基因型并預(yù)測其DUS性狀表明, 基因型預(yù)測在大麥特異性鑒定上具有一定的可行性。

    盡管分子標(biāo)記在品種鑒別上發(fā)揮重要作用, 但目前分子標(biāo)記仍未能取代形態(tài)學(xué)性狀, 單獨用于DUS測試中申請品種的特異性鑒定。這是由于分子標(biāo)記與形態(tài)學(xué)性狀之間尚不具備足夠的相關(guān)性。品種之間的DNA指紋即使存在差異, 其形態(tài)表現(xiàn)也可能相同; 而在某些情況下, 品種之間形態(tài)學(xué)上盡管存在很大差異, 也可能具有相同的DNA指紋, 如突變體。此外, 根據(jù)UPOV公約的要求, 新品種申請保護(hù)除了要求特異性之外, 還必須要求一致性和穩(wěn)定性, 而關(guān)于一致性的分子驗證問題, 由于DNA檢測技術(shù)的靈敏性較高, 難于發(fā)現(xiàn)在所有的標(biāo)記位點都穩(wěn)定的品種, 因此, 利用分子標(biāo)記鑒定品種一致性的標(biāo)準(zhǔn)可能較田間表型一致性的標(biāo)準(zhǔn)要低。但是合適的標(biāo)準(zhǔn)和接受概率如何確定, 對于雜交種來說目前仍是一個難題。此外, 對于一致性的確定應(yīng)該選用多少個體, 是否和DUS 測試需要的田間檢測樣本相同, 如何判別差異個體的容許誤差等, 都仍需在實踐中探索[10]。

    由于分子標(biāo)記單獨應(yīng)用于DUS測試仍存在諸多技術(shù)難題, 形態(tài)學(xué)性狀和分子標(biāo)記的結(jié)合將是近期DUS測試研究的主要方向。將分子標(biāo)記和形態(tài)指標(biāo)綜合利用起來描述PVP (Plant Variety Protection)當(dāng)中的品種性狀, 既有利于認(rèn)識品種之間的遺傳關(guān)系, 也可以滿足植物新品種保護(hù)中明確分類的需要。

    鑒于此, 本研究擬同時利用形態(tài)學(xué)性狀和分子標(biāo)記分析南方花生區(qū)試品種的遺傳多樣性, 并通過聚類分析鑒定花生品種的特異性, 為花生DUS測試中近似品種的快速篩選以及品種質(zhì)量管理提供技術(shù)支撐。

    1 材料與方法

    1.1 材料

    選用1999—2015年南方花生區(qū)試品種(表1), 共101個, 由廣東省農(nóng)業(yè)科學(xué)院作物研究所提供。

    1.2 試驗方法

    1.2.1 形態(tài)學(xué)性狀調(diào)查 2016年3月至2017年8月在華南農(nóng)業(yè)大學(xué)五山試驗基地按隨機(jī)區(qū)組設(shè)計和種植試材, 每個小區(qū)畦長3 m, 畦面寬1.4 m, 行距25 cm, 株距20 cm, 重復(fù)2次。每個品種以40株為樣本, 按農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn)NY/T 2237-2012[11]調(diào)查29個形態(tài)學(xué)性狀, 其中質(zhì)量性狀9個、假質(zhì)量性狀3個、數(shù)量性狀17個。通過目測調(diào)查質(zhì)量性狀和假質(zhì)量性狀, 以1/2、1/2/3或者1/9級記錄; 用直尺或游標(biāo)卡尺測量長度, 電子秤測量質(zhì)量, 對照標(biāo)準(zhǔn)品種, 按1~9級記錄。

    表1 供試品種數(shù)量及來源

    1.2.2 DNA提取方法 取花生幼葉, 采用新型快速植物基因組DNA提取試劑盒(DP3111, BioTeke, 中國)提取花生基因組DNA。以1%瓊脂糖凝膠電泳檢測DNA質(zhì)量, 用超微量分光光度計Nanodrop 2000 (Thermo Scientific, USA)檢測樣本濃度。將所有樣本DNA用無菌水統(tǒng)一稀釋至25 ng μL–1,-20℃保存?zhèn)溆谩?/p>

    1.2.3 SSR標(biāo)記選擇及PCR檢測 從Shirasawa等[12]構(gòu)建的包含3693個標(biāo)記的花生栽培種綜合遺傳圖譜中篩選出40個多態(tài)性高、擴(kuò)增穩(wěn)定、帶型清晰的SSR標(biāo)記(每個染色體2個), SSR引物由北京六合 華大基因科技股份有限公司合成。

    PCR反應(yīng)體系為10 μL, 包含DNA模板1 μL、2×PowerPCR MasterMix (PR1702, BioTeke, 中國) 5 μL、上下游引物各0.4 μL、無菌水3.2 μL。PCR擴(kuò)增程序為95℃預(yù)變性5 min; 95℃變性30 s, 55℃退火30 s, 72℃延伸30 s, 35個循環(huán); 72℃延伸10 min, 擴(kuò)增產(chǎn)物經(jīng)Fragment Analyzer全自動毛細(xì)管電泳系統(tǒng)(FSV2-CE, AATI, USA)檢測后將數(shù)據(jù)結(jié)果自動存于ProSize 2.0系統(tǒng)軟件。

    1.3 數(shù)據(jù)統(tǒng)計與分析

    用NTSYS-pc2.11軟件分析數(shù)據(jù), 在Similarity模塊中選擇Qualitative data分別對形態(tài)學(xué)性狀和SSR標(biāo)記的原始矩陣求遺傳相似系數(shù)(genetic similarity, GS), 并獲得相應(yīng)的GS矩陣。在Clustering模塊中選擇SAHN采用UPGMA進(jìn)行聚類分析, 并通過Tree plot模塊生成聚類圖。在Graphics模塊中選擇Matrix comparison plot對形態(tài)學(xué)性狀和SSR標(biāo)記的相似系數(shù)矩陣之間的相關(guān)性進(jìn)行Mantel檢驗。

    2 結(jié)果與分析

    2.1 形態(tài)學(xué)性狀多樣性分析

    對101個花生品種29個表型性狀進(jìn)行描述性統(tǒng)計(表2)表明, 有7個形態(tài)學(xué)性狀(植株開花習(xí)性、主莖開花習(xí)性、小葉性狀、花冠顏色、植株生長習(xí)性、莢果籽仁數(shù)、種皮顏色數(shù)量)在供試品種間無多樣性, 其余22個性狀Simpson多樣性指數(shù)為0.23~0.77, 平均為0.43, 多樣性指數(shù)最高的性狀為百仁重。

    表2 供試品種的形態(tài)學(xué)性狀特征

    盡管總體上花生品種的形態(tài)學(xué)性狀存在一定的多樣性, 但同一育種機(jī)構(gòu)育成的花生品種表型較相似, 遺傳多樣性低。例如“粵油”系列品種多數(shù)具有植株偏矮、葉色深綠、生育期中等、莢果大、果紋深、籽仁率偏低、莢果縊縮程度弱等特點。

    利用22個具有多樣性的形態(tài)學(xué)性狀對供試品種進(jìn)行聚類分析, 從圖1可知, 在相似系數(shù)為0.76處, 可將供試品種劃分為七大類群, 類群I包含10個品種, 其中“湛油”系列品種9個、“粵油”1個。類群II包含28個品種, 其中“粵油”15個、“仲愷花”4個、“賀油”3個、“閩花”和“航花”各2個、“汕油”和“?;ā备?個。類群III包含30個品種, 其中“泉花”8個、“龍花”7個、“莆花”和“閩花”各5個、“?;ā焙汀敖鸹ā备?個、“賀油”1個。類群IV包含17個品種, 其中“汕油”9個、“仲愷花”4個、“賀油”3個、“閩花”1個。類群V包含12個品種, 其中“桂花”11個、“賀油”1個。類群VI和類群VII各包含2個品種, 分別為“湘花120”、“湘花2008”和“云花生4號”、“云花生12號”。同一育種機(jī)構(gòu)育成的花生品種傾向于聚在一起, 例如湖南農(nóng)業(yè)大學(xué)和云南省農(nóng)業(yè)科學(xué)院分別育成的“湘花”和“云花生”系列品種均單獨聚為一類, 廣西農(nóng)業(yè)科學(xué)院經(jīng)濟(jì)作物研究所育成的11個“桂花”品種均聚在一起, 湛江市農(nóng)業(yè)科學(xué)研究院育成的9個“湛油”品種也聚為一類。此外, 同一省份的花生品種也傾向于聚在一起。例如, 類群III包含30個品種, 其中29個來自福建省。類群II包含29個品種, 其中22個來自廣東省。

    2.2 SSR多態(tài)性分析

    利用40個多態(tài)性高、帶型容易識別的SSR標(biāo)記對101個供試品種基因組DNA進(jìn)行擴(kuò)增, 共檢測出167個等位基因, 單個標(biāo)記檢測的等位基因數(shù)為2~6個, 平均為4.18個。標(biāo)記的多態(tài)性信息量(PIC)差異較大, 最大為0.79, 最小為0.26, 平均為0.55 (表3)。

    基于40對SSR引物擴(kuò)增的167個等位基因, 對101個花生品種進(jìn)行聚類分析。由圖3可知, 在相似系數(shù)為0.70處, 供試材料可被劃分為六大類群。類群Ⅰ包含30個花生品種, 其中“泉花”系列品種8個、“龍花”7個、“莆花”和“閩花”各5個、“?;ā焙汀敖鸹ā备?個、“賀油”1個。類群II包含45個品種, 其中“粵油”15個、“汕油”品種9個、“仲愷花”8個、“賀油”4個、“湛油”和“閩花”各3個、“航花”2個、“?;ā?個。類群Ⅲ包含8個品種, 其中“湛油”6個、“粵油”和“汕油”各1個。類群IV包含14個品種, 其中“桂花”11個、“賀油”3個。類群V包含2個品種, 分別為“湘花120”和“湘花2008”。類群VI也只包含2個品種, 分別為“云花生4號”和“云花生12號”。同一省份育成的品種傾向于聚為一類。例如類群I有27個品種來自福建省, 類群II有37個品種來自廣東省, 而類型III、IV、V、VI中的品種則分別來自廣東省、廣西壯族自治區(qū)、湖南省和云南省。

    2.3 花生品種特異性鑒定

    比較形態(tài)學(xué)性狀和SSR標(biāo)記的聚類結(jié)果發(fā)現(xiàn), 兩者在總體上傾向一致, 同一地方育種的品種多聚為一類。但具體到品種與品種間的遺傳相似性, 兩者卻存在較大差異。供試品種之間的形態(tài)學(xué)性狀相似系數(shù)為0.62~1.00, 平均為0.75, SSR標(biāo)記相似系數(shù)為0.55~1.00, 平均為0.70。經(jīng)Mantel檢驗, 供試品種間的形態(tài)學(xué)性狀相似系數(shù)與SSR標(biāo)記相似系數(shù)的相關(guān)性弱(= 0.36)(圖3), 說明通過形態(tài)學(xué)和SSR鑒定品種的特異性, 結(jié)果將不完全一致。由圖1和圖2可知, 通過形態(tài)學(xué)和SSR可區(qū)分的品種均為97個, 無法區(qū)分的品種為4個。其中“賀油11”與“賀油9號”均無法通過形態(tài)學(xué)和SSR區(qū)分(品種間的相似系數(shù)=1.00), 而“粵油7號”和“粵油410”能通過SSR區(qū)分, 但不能通過形態(tài)學(xué)區(qū)分, 相反, “汕油誘1號”與“汕油輻1號”能通過形態(tài)學(xué)區(qū)分, 但不能通過SSR標(biāo)記區(qū)分。

    由于“賀油11”與“賀油9號”在形態(tài)學(xué)和SSR標(biāo)記上均沒有顯著性差異, 初步推斷為同一品種。盡管“粵油7號”和“粵油410”在形態(tài)學(xué)上相同(“粵油7號”為“粵油410”母本), 但SSR標(biāo)記存在差異, 共檢測出9個差異位點, 品種間的相似系數(shù)為0.83, 因此可判斷“粵油7號”和“粵油410”為不同花生品種。此外, “汕油誘1號”和“汕油輻1號”是“汕油212”分別通過EMS (甲基磺酸乙酯)化學(xué)誘變和Co(鈷)輻射誘變選育出來的品種, SSR檢測發(fā)現(xiàn)2個品種的DNA指紋一致。但由于兩者在形態(tài)學(xué)上存在較大差異, 相似性系數(shù)僅為0.73, 通過形態(tài)學(xué)可區(qū)分開來, 因此, 綜合SSR和形態(tài)學(xué)信息, 可判斷“汕油誘1號”和“汕油輻1號”是來自“汕油212”的突變體。

    3 討論

    3.1 花生品種的遺傳多樣性

    本研究調(diào)查供試品種的29個形態(tài)學(xué)性狀, 結(jié)果有7個性狀無遺傳多樣性, 其余性狀的Simpson多樣性指數(shù)平均僅為0.43。說明當(dāng)前培育的品種在表型上存在同質(zhì)化趨勢。相比之下, 40個SSR標(biāo)記均能在區(qū)試品種中檢測出多態(tài)性, 且多態(tài)性信息量平均為0.55, 遠(yuǎn)高于形態(tài)學(xué)性狀的Simpson多樣性指數(shù)(兩者計算公式一致)。本研究中SSR標(biāo)記比形態(tài)學(xué)性狀的遺傳多樣性更豐富, 其原因之一可能是形態(tài)學(xué)性狀在育種過程中受到強(qiáng)烈的人工選擇, 而SSR標(biāo)記無法通過肉眼觀察, 在育種過程中受到人工選擇的影響較小。

    圖1 形態(tài)學(xué)性狀聚類分析圖

    圖3 形態(tài)學(xué)性狀和SSR標(biāo)記遺傳相似性系數(shù)矩陣比較

    形態(tài)學(xué)聚類顯示, 同一單位培育的品種傾向于聚在一起。造成這種局面的原因有2個: (1) 品種存在親緣關(guān)系。分析歷年南方花生區(qū)試品種的系譜發(fā)現(xiàn), 多數(shù)區(qū)試品種的親本之一為參試單位早期培育的品種, 同一單位培育的品種普遍存在親緣關(guān)系。(2) 不同育種單位的選育方法和目標(biāo)不同。不同育種單位在長期的育種過程中形成了自己的育種經(jīng)驗, 對目標(biāo)性狀的選擇存在偏好, 容易培育出自家特色的花生品種。例如汕頭市科學(xué)研究所傾向于選擇株型直立、莖稈纖細(xì)、成熟期底部落葉明顯、中果型、果多、果紋中等、結(jié)莢集中的品種。

    3.2 花生品種特異性鑒定

    分子標(biāo)記與形態(tài)學(xué)性狀之間不具備足夠的相關(guān)性是當(dāng)前分子標(biāo)記單獨應(yīng)用于植物新品種特異性鑒定的瓶頸。由于形態(tài)學(xué)性狀多為數(shù)量性狀, 受微效多基因控制, 而SSR標(biāo)記通常位于基因組中的非轉(zhuǎn)錄區(qū), 被認(rèn)為是“中性“標(biāo)記, 不具備明顯的生物學(xué)功能。因此, 理論上尋找與形態(tài)學(xué)性狀高度相關(guān)的SSR標(biāo)記存在很大難度。本研究通過Mantel檢驗發(fā)現(xiàn)品種間的形態(tài)學(xué)性狀相似系數(shù)與SSR標(biāo)記相似系數(shù)的相關(guān)性很弱(= 0.36), 說明供試的SSR標(biāo)記暫時無法取代形態(tài)學(xué)性狀, 單獨用于花生品種的特異性鑒定。

    近年來有多項研究表明, 分子標(biāo)記和形態(tài)學(xué)性狀可在植物分類[13]、品種鑒別[14]和遺傳多樣性分析[15-16]中起到互補(bǔ)作用, 尤其對種子真實性鑒定、同種異名和異種同名鑒別[17-18]。理論上, 常規(guī)雜交育種過程中, 親本染色體的重組、交換使雜交后代在遺傳組成上幾乎有無限多的可能性。但由于品種的選育多集中在產(chǎn)量、品質(zhì)、抗性、株型等性狀上, 導(dǎo)致與上述性狀連鎖緊密的DNA標(biāo)記位點的多態(tài)性降低。然而在花生基因組中仍存在著大量與上述性狀不在同一染色體上或連鎖不緊密的DNA標(biāo)記位點, 其多態(tài)性可以將通過常規(guī)雜交育種手段選育出的品種區(qū)分開。例如本研究中用形態(tài)學(xué)性狀沒法將“粵油7號”和“粵油410”區(qū)分開, 但用SSR標(biāo)記則能清晰地鑒定出來。但DNA標(biāo)記在品種的特異性鑒定上并非萬能, 對于一些突變體, 特別是點突變體, 有時就無法利用DNA標(biāo)記區(qū)分開[17]。相反, 這時根據(jù)表型突變性狀, 能輕易與其他品種區(qū)分開。例如本研究中SSR標(biāo)記無法區(qū)分突變體“汕油誘1號”和“汕油輻1號”, 通過形態(tài)學(xué)性狀則可輕易將其區(qū)分開。因此, 綜合形態(tài)學(xué)信息和DNA標(biāo)記信息能有效提高品種特異性鑒定的準(zhǔn)確性。

    4 結(jié)論

    南方花生區(qū)試品種形態(tài)學(xué)性狀的遺傳多樣性比SSR標(biāo)記低, 同一育種單位或者同一省份的品種傾向于聚在一起?;ㄉ鶶SR標(biāo)記與形態(tài)學(xué)性狀不具備足夠的相關(guān)性, 暫時無法取代形態(tài)學(xué)性狀單獨用于品種特異性鑒定, 但兩者相結(jié)合能有效提高花生品種特異性鑒定的準(zhǔn)確性。

    [1] 中華人民共和國農(nóng)業(yè)農(nóng)村部種子管理局. 中國種業(yè)大數(shù)據(jù)平臺. 北京:中華人民共和國農(nóng)業(yè)農(nóng)村部, 2017 [2018-04-01]. http://202.127.42.47:6009/Home/BigDataIndex. Bureau of Seed Management of Ministry of Agriculture and Rural Affairs of the People’s Republic of China. China Seed Industry Big Data Platform. Beijing: Ministry of Agriculture and Rural Affairs of the People’s Republic of China, 2017 [2018-04-01]. http://202.127.42.47:6009/Home/BigDataIndex.

    [2] 中華人民共和國國務(wù)院. 中華人民共和國植物新品種保護(hù)條例. 北京: 中國農(nóng)業(yè)出版社, 1997. pp 4–5.The State Council of the People’s Republic of China. Regulations of the People’s Republic of China Plant Variety Protection. Beijing: China Agriculture Press, 1997. pp 4–5 (in Chinese).

    [3] 魯清, 李少雄, 陳小平, 周桂元, 洪彥彬, 李海芬, 梁炫強(qiáng). 我國南方產(chǎn)區(qū)花生育種現(xiàn)狀、存在問題及育種建議. 中國油料作物學(xué)報, 2017, 39: 556–566. Lu Q, Li S X, Chen X P, Zhou G Y, Hong Y B, Li H F, Liang X Q. Current situation, problems and suggestions of peanut breeding in southern China., 2017, 39: 556–566 (in Chinese with English abstract).

    [4] Hayward A C, Tollenaere R, Dalton-Morgan J, Batley J. Molecular marker applications in plants., 2015, 1245: 13–27.

    [5] Cockram J, Jones H, Norris C, O'Sullivan D M. Evaluation of diagnostic molecular markers for DUS phenotypic assessment in the cereal crop, barley (ssp.L.)., 2012, 125: 1735–1749

    [6] Jones H, Norris C, Smith D, Cockram J, Lee D, O’Sullivan D M, Mackay I. Evaluation of the use of high-density SNP genotyping to implement UPOV Model 2 for DUS testing in barley., 2013, 126: 901–911.

    [7] Tommasini L, Batley J, Arnold G M, Cooke R J, Donini P, Lee D, Law J R, Lowe C, Moule C, Trick M, Edwards K J. The deve-lopment of multiplex simple sequence repeat (SSR) markers to complement distinctness, uniformity and stability testing of rape (L.) varieties., 2003, 106: 1091–1101.

    [8] Arens P, Mansilla C, Deinum D, Cavellini L, Moretti A, Rolland S, van der Schoot H, Calvache D, Ponz F, Collonnier C, Mathis R, Smilde D, Caranta C, Vosman B. Development and evaluation of robust molecular markers linked to disease resistance in tomato for distinctness, uniformity and stability testing., 2010, 120: 655–664.

    [9] Jones H, Mackay I. Implications of using genomic prediction within a high-density SNP dataset to predict DUS traits in barley., 2015, 128: 2461–2470.

    [10] 滕海濤, 呂波, 趙久然, 徐巖, 王鳳格, 堵苑苑, 楊坤, 唐浩, 李祥羽. 利用DNA指紋圖譜輔助植物新品種保護(hù)的可能性. 生物技術(shù)通報, 2009, (1): 1–6. Teng H T, Lyu B, Zhao J R, Xu Y, Wang F G, Du Y Y, Yang K, Tang H, Li X Y. DNA fingerprint profile involved in plant variety protection practice., 2009, (1): 1–6 (in Chinese with English abstract).

    [11] 中華人民共和國農(nóng)業(yè)農(nóng)村部. 植物新品種特異性、一致性和穩(wěn)定性測試指南: 花生. 北京: 中國農(nóng)業(yè)出版社, 2012. pp 4–7. Ministry of Agriculture and Rural Affairs of the People’s Republic of China. Guidelines for the Conduct of Tests for Distinctness, Uniformity and Stability: Peanut (L.). Beijing: China Agriculture Press, 2012. pp 4–7 (in Chinese).

    [12] Shirasawa K, Bertioli D J, Varshney R K, Moretzsohn M C, Leal-Bertioli S C M, Thudi M, Pandey M K, Rami J F, Foncéka D, Gowda M V C, Qin H, Guo B, Hong Y, Liang X, Hirakawa H, Tabata S, Isobe S. Integrated consensus map of cultivated peanut and wild relatives reveals structures of the A and B genomes ofand divergence of the legume genomes., 2013, 20: 173–184.

    [13] Heilmann-Clausen J, Christensen M, Fr?slev T G, Kj?ller R. Taxonomy of Tricholoma in northern Europe based on ITS sequence data and morphological characters., 2017, 38: 38–57.

    [14] Smykal P, Horacek J, Dostalova R, Hybl M. Variety discrimination in pea (L.) by molecular, biochemical and morphological markers., 2008, 49: 155–166.

    [15] Ebrahimi A, Zarei A, Zamani F M, Lawson S. Evaluation of genetic variability among “Early Mature”using microsatellite markers and morphological traits., 2017, 5: e3834.

    [16] Rebaa F, Abid G, Aouida M, Abdelkarim S, Aroua I, Muhovski Y, Baudoin J P, M’hamdi M, Sassi K, Jebara M. Genetic variability in Tunisian populations of faba bean (L. var) assessed by morphological and SSR markers., 2017, 23: 397–409.

    [17] 李瑞峰, 高鵬, 朱子成, 欒非時. 基于形態(tài)學(xué)標(biāo)記及SSR標(biāo)記的甜瓜主栽品種分類鑒定研究. 中國蔬菜, 2014, (6): 20–27. Li R F, Gao P, Zhu Z C, Luan F S. Studies on classification and identification based on morphological markers and SSR markers for elite varieties ofL., 2014, (6): 20–27 (in Chinese with English abstract).

    [18] 李鶴, 郭世榮, 束勝,徐揚(yáng), 孫錦. 砧用南瓜種質(zhì)資源形態(tài)學(xué)性狀與SSR標(biāo)記分析. 園藝學(xué)報, 2014, 41: 1379–1390. Li H, Guo S R, Shu S, Xu Y, Sun J. Germplasm resources analysis of rootstock-used pumpkins by phenotype and SSR., 2014, 41: 1379–1390 (in Chinese with English abstract).

    Genetic diversity analysis and distinctness identification of peanut cultivars based on morphological traits and SSR markers

    LIU Hong1,**, XU Zhen-Jiang1,**, RAO De-Hua1, LU Qing2,3, LI Shao-Xiong2,3, LIU Hai-Yan2,3, CHENXiao-Ping2,3, LIANG Xuan-Qiang2,3, and HONG Yan-Bin2,3,*

    1College of Agriculture, South China Agricultural University, Guangzhou 510642, Guangdong, China;2Crops Research Institute, Guangdong Academy of Agricultural Sciences, Guangzhou 510640, Guangdong, China;3Guangdong Provincial Key Laboratory of Crops Genetics and Improvement, Guangzhou 510640, Guangdong, China

    The genetic diversity and distinctness of 101 peanut varieties participated in the South China peanut field trial were analyzed based on their morphological traits and SSR markers. Among 29 morphological traits seven were no difference while 22 demonstrated diversity indexes ranging from 0.23 to 0.77, with an average of 0.43. The varieties were divided into seven groups at the similarity coefficient of 0.76, and those released by the same breeding institution tended to converge into one cluster. Molecular characterization with 40 highly informative SSRs generated a total of 167 alleles ranging from 2 to 6 (averaged 4.18) alleles per marker. The polymorphism information content (PIC) of these markers varied from 0.79 to 0.26 with an average of 0.55/marker. The varieties were divided into six groups based on the similarity coefficient of 0.70, and those released by the same province tended to converge into one cluster. Mantel testing revealed that the correlations of the similarity coefficient matrixes between the morphological traits and SSR markers were weak (= 0.36), implying that SSR markers are unable to replace morphological traits to be solely adopted to identify the distinctness of peanut varieties, but the combination of morphological traits and SSR markers will effectively increase the accuracy of distinctiveness identification.

    peanut; morphology; SSR; genetic diversity; DUS

    2018-04-27;

    2018-08-20;

    2018-09-19.

    10.3724/SP.J.1006.2019.84060

    通信作者(Corresponding author): 洪彥彬, E-mail: hongyanbin@gdaas.cn

    **同等貢獻(xiàn)(Contributed equally to this work)

    劉洪, E-mail: laoliuhongscau@163.com; 徐振江, E-mail: zhenjiangxu521@scau.edu.cn

    本研究由國家自然科學(xué)基金項目(31771841), 2015年品種資源保護(hù)項目(2015-18), 廣東省科技計劃項目(2013B020301014, 2013B050800021, 2017A030311007, 2016B020201003, 2015A030313565)和廣東省現(xiàn)代農(nóng)業(yè)創(chuàng)新聯(lián)盟建設(shè)項目(2016LM3161)資助。

    This study was supported by the National Natural Science Foundation of China (31771841), 2015 Variety Resource Protection Project (2015-18), Science and Technology Plan Project of Guangdong Province (2013B020301014, 2013B050800021, 2017A030311007, 2016B020201003), and Modern Agricultural Science and Technology Innovation Alliance Construction Project of Guangdong Province (2016LM3161).

    URL:http://kns.cnki.net/kcms/detail/11.1809.S.20180918.1354.008.html

    猜你喜歡
    類群形態(tài)學(xué)花生
    掏花生
    薏苡種質(zhì)資源ISSR分子標(biāo)記篩選及親緣關(guān)系分析
    黑石頂自然保護(hù)區(qū)土壤動物功能類群對季節(jié)動態(tài)的響應(yīng)
    到底埋在哪棵樹下
    花生去哪兒了
    醫(yī)學(xué)微觀形態(tài)學(xué)在教學(xué)改革中的應(yīng)用分析
    雨后的露營
    2011年春夏季黃海和東海微型浮游動物類群組成及其攝食的研究
    數(shù)學(xué)形態(tài)學(xué)濾波器在轉(zhuǎn)子失衡識別中的應(yīng)用
    龍灣自然保護(hù)區(qū)森林土壤動物群落多樣性及功能類群
    欧美色视频一区免费| 日本黄大片高清| av片东京热男人的天堂| 啪啪无遮挡十八禁网站| 国产精品av久久久久免费| 听说在线观看完整版免费高清| 亚洲国产欧美一区二区综合| 在线观看免费午夜福利视频| 欧美另类亚洲清纯唯美| 两个人看的免费小视频| 免费一级毛片在线播放高清视频| 婷婷亚洲欧美| 男人舔女人的私密视频| 国产高清激情床上av| 久久久国产精品麻豆| 国产精品久久久久久人妻精品电影| 99热6这里只有精品| 国产人伦9x9x在线观看| 亚洲最大成人中文| 熟女电影av网| 日韩欧美国产一区二区入口| 一区福利在线观看| 亚洲精品美女久久av网站| 亚洲中文av在线| 国产伦人伦偷精品视频| 精品一区二区三区视频在线 | 欧美zozozo另类| 精品国产乱码久久久久久男人| 国内揄拍国产精品人妻在线| 亚洲真实伦在线观看| 亚洲欧美日韩无卡精品| 免费av毛片视频| 亚洲中文av在线| cao死你这个sao货| 久久午夜亚洲精品久久| 一进一出抽搐动态| 精品久久久久久成人av| 亚洲av成人精品一区久久| 久久精品亚洲精品国产色婷小说| 又紧又爽又黄一区二区| 伊人久久大香线蕉亚洲五| 一a级毛片在线观看| 深夜精品福利| 叶爱在线成人免费视频播放| 久久精品国产99精品国产亚洲性色| 午夜日韩欧美国产| 国产三级中文精品| 国产精品女同一区二区软件 | 日韩欧美精品v在线| 舔av片在线| 热99re8久久精品国产| 美女被艹到高潮喷水动态| 99久久久亚洲精品蜜臀av| 长腿黑丝高跟| 99热6这里只有精品| 日韩欧美在线乱码| 男女午夜视频在线观看| 免费搜索国产男女视频| 真实男女啪啪啪动态图| ponron亚洲| 麻豆国产97在线/欧美| 麻豆国产av国片精品| 国产综合懂色| 三级国产精品欧美在线观看 | 久久久久免费精品人妻一区二区| 亚洲精品一区av在线观看| 久久中文字幕人妻熟女| 欧美另类亚洲清纯唯美| 日本黄色视频三级网站网址| 国产精品一及| 精品熟女少妇八av免费久了| 亚洲一区二区三区不卡视频| 最新中文字幕久久久久 | 9191精品国产免费久久| 丁香欧美五月| 精品久久久久久久末码| 夜夜躁狠狠躁天天躁| 国产99白浆流出| 日本精品一区二区三区蜜桃| 免费在线观看亚洲国产| 又黄又爽又免费观看的视频| 丰满的人妻完整版| 熟女人妻精品中文字幕| 亚洲无线观看免费| 久久精品夜夜夜夜夜久久蜜豆| 免费人成视频x8x8入口观看| 亚洲av成人精品一区久久| 身体一侧抽搐| 久久久水蜜桃国产精品网| 欧美色视频一区免费| 国产av麻豆久久久久久久| 97碰自拍视频| 岛国在线免费视频观看| 国产在线精品亚洲第一网站| 99国产极品粉嫩在线观看| 国产精品自产拍在线观看55亚洲| 国产亚洲欧美98| 中文资源天堂在线| 一区二区三区激情视频| e午夜精品久久久久久久| 久久草成人影院| 中文资源天堂在线| 伊人久久大香线蕉亚洲五| 丰满的人妻完整版| 在线观看美女被高潮喷水网站 | 国产精品亚洲美女久久久| 成人av一区二区三区在线看| 亚洲国产欧洲综合997久久,| 国产亚洲av嫩草精品影院| 国产激情久久老熟女| 欧美高清成人免费视频www| 身体一侧抽搐| av天堂在线播放| 丰满的人妻完整版| 亚洲精品美女久久av网站| 日本免费a在线| 欧美成狂野欧美在线观看| 国产单亲对白刺激| 啦啦啦观看免费观看视频高清| 少妇熟女aⅴ在线视频| 黑人巨大精品欧美一区二区mp4| 国产精品久久久久久亚洲av鲁大| 国产精品久久久久久精品电影| 99国产综合亚洲精品| 久久午夜亚洲精品久久| 成熟少妇高潮喷水视频| 中文字幕av在线有码专区| 精品国产超薄肉色丝袜足j| 亚洲美女视频黄频| 欧美成人一区二区免费高清观看 | 久久精品人妻少妇| 精品免费久久久久久久清纯| svipshipincom国产片| 脱女人内裤的视频| 淫秽高清视频在线观看| 日本熟妇午夜| 动漫黄色视频在线观看| 国产乱人视频| 搡老岳熟女国产| 欧美一区二区国产精品久久精品| 久久久色成人| 精品国内亚洲2022精品成人| 国内少妇人妻偷人精品xxx网站 | xxxwww97欧美| 看免费av毛片| 成人国产综合亚洲| 日韩中文字幕欧美一区二区| 精品人妻1区二区| 99热只有精品国产| 欧美激情在线99| 国产成人av教育| 欧美大码av| 欧美不卡视频在线免费观看| 97碰自拍视频| АⅤ资源中文在线天堂| 日本 欧美在线| 久久精品国产99精品国产亚洲性色| 精品久久久久久久久久免费视频| 久久久久国产精品人妻aⅴ院| 婷婷精品国产亚洲av在线| 日本 av在线| 日本一二三区视频观看| 成人永久免费在线观看视频| 亚洲精品色激情综合| 在线观看免费午夜福利视频| 婷婷亚洲欧美| 久久精品综合一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 午夜福利在线观看免费完整高清在 | 小蜜桃在线观看免费完整版高清| 亚洲美女视频黄频| av片东京热男人的天堂| 夜夜看夜夜爽夜夜摸| 国产精品爽爽va在线观看网站| 老司机福利观看| 精品一区二区三区视频在线观看免费| 淫秽高清视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 一进一出抽搐gif免费好疼| 黄频高清免费视频| 一级毛片精品| 久久精品国产清高在天天线| 欧美黄色淫秽网站| 99久久久亚洲精品蜜臀av| 国产黄色小视频在线观看| 级片在线观看| 18禁黄网站禁片免费观看直播| www.www免费av| 中文字幕av在线有码专区| 一区福利在线观看| av女优亚洲男人天堂 | 性色avwww在线观看| 真人做人爱边吃奶动态| 国产91精品成人一区二区三区| 女警被强在线播放| 亚洲专区国产一区二区| 国产精品日韩av在线免费观看| 99久久综合精品五月天人人| 精品一区二区三区视频在线观看免费| 国内毛片毛片毛片毛片毛片| 可以在线观看毛片的网站| 五月伊人婷婷丁香| 国产美女午夜福利| 夜夜看夜夜爽夜夜摸| 国产精品女同一区二区软件 | 99国产综合亚洲精品| av黄色大香蕉| 啦啦啦免费观看视频1| 精品欧美国产一区二区三| 特级一级黄色大片| 日韩三级视频一区二区三区| 亚洲第一电影网av| 少妇人妻一区二区三区视频| 18禁黄网站禁片午夜丰满| 国语自产精品视频在线第100页| 国产一区二区激情短视频| 精品午夜福利视频在线观看一区| 999久久久国产精品视频| 色播亚洲综合网| 成人一区二区视频在线观看| 在线观看一区二区三区| 免费在线观看成人毛片| 成人av一区二区三区在线看| 国产成人影院久久av| 国产伦一二天堂av在线观看| 哪里可以看免费的av片| 亚洲国产精品久久男人天堂| 欧美av亚洲av综合av国产av| 亚洲色图 男人天堂 中文字幕| 天天一区二区日本电影三级| 亚洲精品粉嫩美女一区| 男女那种视频在线观看| 美女免费视频网站| 国内精品久久久久久久电影| 国产亚洲精品久久久久久毛片| 国内毛片毛片毛片毛片毛片| 成人精品一区二区免费| 深夜精品福利| 亚洲avbb在线观看| 欧美日本亚洲视频在线播放| 国产精品美女特级片免费视频播放器 | 欧美性猛交黑人性爽| svipshipincom国产片| 午夜免费成人在线视频| 久9热在线精品视频| 老司机在亚洲福利影院| 亚洲av电影在线进入| 最新美女视频免费是黄的| 亚洲乱码一区二区免费版| 国产亚洲精品久久久com| 日韩国内少妇激情av| 久久香蕉国产精品| 日本一本二区三区精品| 丰满的人妻完整版| www.精华液| 日韩欧美国产一区二区入口| 青草久久国产| 在线a可以看的网站| 亚洲性夜色夜夜综合| 亚洲av成人精品一区久久| 在线观看一区二区三区| 亚洲第一欧美日韩一区二区三区| 国产av在哪里看| 禁无遮挡网站| 小说图片视频综合网站| 99久久无色码亚洲精品果冻| 亚洲一区二区三区不卡视频| 国产亚洲欧美在线一区二区| 精品一区二区三区av网在线观看| 午夜免费观看网址| 人妻久久中文字幕网| 伊人久久大香线蕉亚洲五| 亚洲五月婷婷丁香| 麻豆成人av在线观看| 观看美女的网站| 亚洲 欧美一区二区三区| 丰满的人妻完整版| 亚洲精品456在线播放app | 欧美成狂野欧美在线观看| 久久性视频一级片| 一a级毛片在线观看| 成人特级av手机在线观看| 黄色片一级片一级黄色片| 免费观看精品视频网站| 男人舔女人的私密视频| 国产精品精品国产色婷婷| 香蕉国产在线看| 成人精品一区二区免费| 色哟哟哟哟哟哟| 亚洲熟妇熟女久久| 在线看三级毛片| 色老头精品视频在线观看| 免费电影在线观看免费观看| 午夜福利视频1000在线观看| 脱女人内裤的视频| 欧美+亚洲+日韩+国产| 在线观看美女被高潮喷水网站 | 一区二区三区高清视频在线| 18禁国产床啪视频网站| 亚洲成av人片免费观看| 人妻久久中文字幕网| 日韩av在线大香蕉| av片东京热男人的天堂| 中文字幕久久专区| 三级国产精品欧美在线观看 | 人妻久久中文字幕网| 九九在线视频观看精品| 长腿黑丝高跟| 久久久久性生活片| 香蕉久久夜色| 久久精品aⅴ一区二区三区四区| 一边摸一边抽搐一进一小说| 国产探花在线观看一区二区| 久久久久亚洲av毛片大全| 日本免费a在线| 日韩欧美在线二视频| 两个人的视频大全免费| 午夜亚洲福利在线播放| 国产精品久久电影中文字幕| 不卡一级毛片| 亚洲狠狠婷婷综合久久图片| 亚洲男人的天堂狠狠| 99国产精品一区二区三区| 丁香六月欧美| 国产精品久久久av美女十八| 麻豆国产av国片精品| 国内精品久久久久精免费| 五月伊人婷婷丁香| 男人和女人高潮做爰伦理| 欧美日韩国产亚洲二区| 国产精品久久久av美女十八| 久久久久九九精品影院| 99久久99久久久精品蜜桃| 国产真人三级小视频在线观看| 免费观看的影片在线观看| 在线十欧美十亚洲十日本专区| 色综合站精品国产| 久久午夜综合久久蜜桃| 99久久久亚洲精品蜜臀av| 亚洲av熟女| 国产av麻豆久久久久久久| 12—13女人毛片做爰片一| av在线天堂中文字幕| 日本撒尿小便嘘嘘汇集6| 国产伦在线观看视频一区| 亚洲片人在线观看| tocl精华| svipshipincom国产片| 少妇熟女aⅴ在线视频| 黄片小视频在线播放| 国产成人一区二区三区免费视频网站| 国产毛片a区久久久久| 岛国视频午夜一区免费看| x7x7x7水蜜桃| 国产精品98久久久久久宅男小说| 日本精品一区二区三区蜜桃| 久9热在线精品视频| 老司机在亚洲福利影院| 国产欧美日韩精品亚洲av| 成人三级做爰电影| 国产高潮美女av| 亚洲黑人精品在线| 在线a可以看的网站| 麻豆av在线久日| 三级男女做爰猛烈吃奶摸视频| а√天堂www在线а√下载| 中文字幕久久专区| 午夜福利高清视频| 黄色成人免费大全| 午夜影院日韩av| 国产成人精品无人区| 国产av一区在线观看免费| 又粗又爽又猛毛片免费看| 日本a在线网址| 高潮久久久久久久久久久不卡| 一级a爱片免费观看的视频| 可以在线观看的亚洲视频| 久久国产精品人妻蜜桃| 首页视频小说图片口味搜索| 一本综合久久免费| 亚洲成人中文字幕在线播放| 嫁个100分男人电影在线观看| 成人精品一区二区免费| 国产主播在线观看一区二区| 国产美女午夜福利| 国语自产精品视频在线第100页| 国产欧美日韩精品亚洲av| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美 国产精品| 偷拍熟女少妇极品色| 欧美大码av| 亚洲自拍偷在线| 一本一本综合久久| 国产精品国产高清国产av| 男女床上黄色一级片免费看| 中出人妻视频一区二区| 变态另类成人亚洲欧美熟女| 久久久精品欧美日韩精品| 国产久久久一区二区三区| 国产99白浆流出| 欧美成人性av电影在线观看| 国产av麻豆久久久久久久| 国产伦精品一区二区三区视频9 | www日本黄色视频网| 久久午夜亚洲精品久久| 日本a在线网址| 午夜福利成人在线免费观看| 欧美一区二区精品小视频在线| 男人的好看免费观看在线视频| 又黄又爽又免费观看的视频| 99re在线观看精品视频| 欧美成狂野欧美在线观看| 亚洲国产精品合色在线| 亚洲av成人av| 法律面前人人平等表现在哪些方面| 亚洲专区中文字幕在线| 欧美黑人巨大hd| 中文资源天堂在线| 1024香蕉在线观看| 精品99又大又爽又粗少妇毛片 | 婷婷六月久久综合丁香| 欧美色欧美亚洲另类二区| 亚洲一区二区三区色噜噜| 久久久久久久精品吃奶| 免费一级毛片在线播放高清视频| 久久精品aⅴ一区二区三区四区| 亚洲欧美激情综合另类| 亚洲精品在线美女| 亚洲精华国产精华精| www日本在线高清视频| 国产高清视频在线观看网站| 精华霜和精华液先用哪个| 精品欧美国产一区二区三| 亚洲av五月六月丁香网| 亚洲性夜色夜夜综合| 99热这里只有是精品50| 九九久久精品国产亚洲av麻豆 | 怎么达到女性高潮| 精品国产美女av久久久久小说| 欧美大码av| АⅤ资源中文在线天堂| 国产熟女xx| 日本a在线网址| 三级男女做爰猛烈吃奶摸视频| 成年人黄色毛片网站| 久久久久久九九精品二区国产| xxxwww97欧美| 男女做爰动态图高潮gif福利片| 日本一本二区三区精品| 亚洲狠狠婷婷综合久久图片| 曰老女人黄片| 精品欧美国产一区二区三| 神马国产精品三级电影在线观看| 一卡2卡三卡四卡精品乱码亚洲| 丰满的人妻完整版| 国内久久婷婷六月综合欲色啪| 男女下面进入的视频免费午夜| 午夜激情欧美在线| 国产99白浆流出| 十八禁人妻一区二区| 国内精品一区二区在线观看| 成人特级av手机在线观看| 日本一本二区三区精品| 啦啦啦免费观看视频1| aaaaa片日本免费| 国产高清视频在线播放一区| www.www免费av| 白带黄色成豆腐渣| 国内毛片毛片毛片毛片毛片| 看片在线看免费视频| 日本精品一区二区三区蜜桃| 国产精品 欧美亚洲| 搡老岳熟女国产| 亚洲成a人片在线一区二区| 日本黄大片高清| 桃红色精品国产亚洲av| 免费无遮挡裸体视频| 黄频高清免费视频| 91九色精品人成在线观看| 久久久久久国产a免费观看| 看免费av毛片| 97超视频在线观看视频| 国产高清videossex| 久久精品夜夜夜夜夜久久蜜豆| 午夜成年电影在线免费观看| 亚洲av成人av| 超碰成人久久| 禁无遮挡网站| 日本与韩国留学比较| 亚洲国产精品久久男人天堂| 女生性感内裤真人,穿戴方法视频| 两性午夜刺激爽爽歪歪视频在线观看| 禁无遮挡网站| 国产av在哪里看| 日韩欧美在线二视频| 99在线人妻在线中文字幕| 搞女人的毛片| 给我免费播放毛片高清在线观看| 最好的美女福利视频网| 国产人伦9x9x在线观看| 中文字幕熟女人妻在线| 国产精品1区2区在线观看.| 亚洲人成电影免费在线| 国产真人三级小视频在线观看| 日本一本二区三区精品| 精品国产乱码久久久久久男人| 精品不卡国产一区二区三区| 久久久久性生活片| 精品国产美女av久久久久小说| 国产激情偷乱视频一区二区| 99国产精品一区二区蜜桃av| 我的老师免费观看完整版| 首页视频小说图片口味搜索| 88av欧美| 欧美绝顶高潮抽搐喷水| 成人午夜高清在线视频| 免费无遮挡裸体视频| 母亲3免费完整高清在线观看| 可以在线观看毛片的网站| 老熟妇乱子伦视频在线观看| 亚洲精品久久国产高清桃花| 99视频精品全部免费 在线 | 老司机深夜福利视频在线观看| 精品免费久久久久久久清纯| 亚洲欧美日韩东京热| 欧美在线一区亚洲| 国产毛片a区久久久久| 久久中文字幕一级| 成人永久免费在线观看视频| 制服丝袜大香蕉在线| 中文资源天堂在线| 综合色av麻豆| 9191精品国产免费久久| a级毛片a级免费在线| 精品国产乱子伦一区二区三区| 国产91精品成人一区二区三区| 国产av一区在线观看免费| 一级a爱片免费观看的视频| 亚洲avbb在线观看| 午夜福利欧美成人| 久久草成人影院| 午夜福利18| 熟女少妇亚洲综合色aaa.| 久久99热这里只有精品18| 国产精品永久免费网站| 深夜精品福利| 精品国产乱子伦一区二区三区| 久久久久久久久久黄片| 亚洲人成网站高清观看| 国产私拍福利视频在线观看| 十八禁人妻一区二区| 我的老师免费观看完整版| av女优亚洲男人天堂 | 操出白浆在线播放| 久久久久久人人人人人| bbb黄色大片| 婷婷精品国产亚洲av在线| 国产97色在线日韩免费| 国产三级在线视频| 欧美极品一区二区三区四区| 亚洲自拍偷在线| 亚洲欧洲精品一区二区精品久久久| 变态另类丝袜制服| 久久久久免费精品人妻一区二区| 欧美激情在线99| 亚洲人成网站在线播放欧美日韩| 国产亚洲av嫩草精品影院| 悠悠久久av| 午夜激情欧美在线| 啦啦啦免费观看视频1| 99热这里只有精品一区 | 国产成人aa在线观看| 人人妻人人澡欧美一区二区| 国产精品 欧美亚洲| 丰满人妻一区二区三区视频av | 成人亚洲精品av一区二区| 午夜福利18| 最近视频中文字幕2019在线8| 亚洲专区字幕在线| 99热精品在线国产| 女生性感内裤真人,穿戴方法视频| 亚洲av片天天在线观看| 久久久久久国产a免费观看| 无限看片的www在线观看| 久久久久久久精品吃奶| 一区二区三区激情视频| 亚洲欧美日韩高清专用| 99国产极品粉嫩在线观看| 国产精品 欧美亚洲| 国产成人福利小说| 一个人看视频在线观看www免费 | a在线观看视频网站| 成人国产一区最新在线观看| 午夜福利成人在线免费观看| 国产成年人精品一区二区| 韩国av一区二区三区四区| 欧美黑人巨大hd| 老熟妇仑乱视频hdxx| 老司机午夜福利在线观看视频| 久久中文看片网| 中文资源天堂在线| 久久天躁狠狠躁夜夜2o2o| 精品电影一区二区在线| 校园春色视频在线观看| 亚洲熟妇中文字幕五十中出| 中文资源天堂在线| 久久久久国内视频| 午夜福利高清视频| 久久久国产成人精品二区| 亚洲片人在线观看| 亚洲av成人精品一区久久| 露出奶头的视频| 国产毛片a区久久久久| 两性夫妻黄色片| 神马国产精品三级电影在线观看|