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

    基于表型性狀構(gòu)建中國花生地方品種骨干種質(zhì)

    2020-03-23 06:30:12閆彩霞李春娟宋秀霞孫全喜苑翠玲趙小波單世華
    作物學(xué)報(bào) 2020年4期
    關(guān)鍵詞:骨干表型種質(zhì)

    閆彩霞 王 娟 張 浩 李春娟 宋秀霞 孫全喜 苑翠玲 趙小波 單世華,*

    基于表型性狀構(gòu)建中國花生地方品種骨干種質(zhì)

    閆彩霞1王 娟1張 浩1李春娟1宋秀霞2孫全喜1苑翠玲1趙小波1單世華1,*

    1山東省花生研究所, 山東青島 266100;2菏澤市牡丹區(qū)農(nóng)業(yè)農(nóng)村局, 山東菏澤 274000

    中國花生地方品種遺傳多樣性豐富, 是花生新品種選育的重要親本來源。本研究以種質(zhì)庫保存的2741份地方種質(zhì)為材料, 基于種植區(qū)劃和植物學(xué)類型分組, 平方根法確定取樣量, 組內(nèi)按13個表型數(shù)據(jù)進(jìn)行UPGMA聚類分析, 類內(nèi)隨機(jī)取樣, 構(gòu)建骨干種質(zhì)。利用檢驗(yàn)、測驗(yàn)、卡方測驗(yàn)、極差、表型保留比例、表型相關(guān)性等對骨干種質(zhì)代表性進(jìn)行檢驗(yàn)和評價(jià); 并利用主成分分析和直方圖對骨干種質(zhì)進(jìn)行確認(rèn)。結(jié)果表明, 構(gòu)建了包含259份種質(zhì)的中國花生地方品種骨干種質(zhì), 占全部種質(zhì)的9.4%, 包括多粒型14份、珍珠豆型85份、龍生型42份、普通型103份、中間型15份。在< 0.05概率條件下, 骨干種質(zhì)13個性狀的均值、方差、變異系數(shù)、香農(nóng)指數(shù)與全部種質(zhì)無顯著差異, 且保留了全部種質(zhì)的分布范圍、表型保留比例和表型相關(guān)性; 二者的植物學(xué)類型組成和生態(tài)分布是一致的, 具有相似的遺傳結(jié)構(gòu)和分布頻率。建立的骨干種質(zhì)很好地代表了全部種質(zhì)的遺傳變異和群體結(jié)構(gòu), 可為花生種質(zhì)創(chuàng)新和優(yōu)異等位基因發(fā)掘奠定良好的基礎(chǔ)。

    花生; 地方品種; 表型性狀; 骨干種質(zhì); 代表性評價(jià)

    花生栽培種(L.)起源于南美洲, 隸屬于豆科花生屬(), 一年生, 異源四倍體(AABB), 推測是由二倍體野生種經(jīng)過自然雜交和染色體加倍而形成的[1]。自多途徑引進(jìn)后, 在我國多變的氣候及復(fù)雜的地理環(huán)境下形成了遺傳多樣性較為豐富的種質(zhì)資源。目前國家中期種質(zhì)庫中的保有量已達(dá)7000多份, 居世界第3位[2], 其中地方品種共4638份, 因其變異廣泛、適應(yīng)性好、配合力高及對病害、蟲害及逆境等的抗性, 是花生品種改良的重要親本來源。追溯建國以來育成花生品種的系譜, 有25個地方品種直接或間接參與育成了85%的品種, 其中40%的親緣來源于伏花生、徐州68-4、獅頭企、粵油551、徐州402和粵選58。地方品種成為花生育種的“骨干親本”[3], 但是, 在育種實(shí)踐中用作親本的不足百份, 特別是多粒型種質(zhì)和龍生型種質(zhì)有效利用甚少[4]。同時, 花生育種存在著嚴(yán)重的“近親繁殖”, 導(dǎo)致品種遺傳基礎(chǔ)狹窄, 適應(yīng)性、抗病(逆)性減退, 產(chǎn)量上很難取得重要突破, 因此, 拓寬現(xiàn)有花生品種的遺傳基礎(chǔ)已迫在眉睫。

    地方品種具有豐富的多樣性, 這為品種選育和遺傳研究提供了廣闊的遺傳基礎(chǔ)。然而, 其龐大的資源數(shù)量反而阻礙了作物的品種改良和新品種選育。構(gòu)建核心種質(zhì)(core collection)或骨干種質(zhì)(key germplasm), 即以最小的資源份數(shù)最大限度地代表整個資源的遺傳多樣性[5], 無疑是解決這一矛盾的有效策略, 這為地方種質(zhì)的深入評價(jià)、創(chuàng)新利用和基因挖掘開辟了新的途徑。迄今為止, 基于表型性狀和分子標(biāo)記數(shù)據(jù), 已構(gòu)建了水稻、小麥、玉米、大豆、棉花、燕麥、谷子、芝麻、糜子、黍稷等[6]多種作物的核心(骨干)種質(zhì)。核心種質(zhì)的取樣比例、取樣策略以及有效性評價(jià)等理論研究也取得了進(jìn)展, 為核心(骨干)種質(zhì)的構(gòu)建及代表性評價(jià)提供了堅(jiān)實(shí)的理論依據(jù)。

    作物的表型多樣性是遺傳多樣性的外在表現(xiàn), 國內(nèi)外研究機(jī)構(gòu)均利用表型數(shù)據(jù)構(gòu)建了花生核心種質(zhì), 如美國的831份核心種質(zhì)[7]和112份微核心種質(zhì)[8], 國際半干旱地區(qū)熱帶作物研究所(ICRISAT)的1704份核心種質(zhì)[9]和184份微核心種質(zhì)[10], 姜慧芳等[2]構(gòu)建的298份小核心種質(zhì), 并且開展了大量的性狀評價(jià)與抗病(蟲或逆)性鑒定[11-17]。本研究以建國以來收集保存的2741份代表性地方品種為材料, 以13個表型性狀數(shù)據(jù)為基礎(chǔ), 構(gòu)建中國花生地方品種骨干種質(zhì), 對其代表性進(jìn)行評價(jià), 為地方品種的保護(hù)與利用、花生種質(zhì)的創(chuàng)新利用及品種的遺傳改良提供理論依據(jù)。

    1 材料與方法

    1.1 材料

    2741份有明確地理來源的中國栽培花生地方品種, 其有關(guān)數(shù)據(jù)來自山東省花生研究所編寫的《中國花生品種資源目錄》及中國農(nóng)業(yè)科學(xué)院油料作物研究所編寫的《中國花生品種資源目錄(續(xù)編一、二)》。涉及生育期、百果重、百仁重、出仁率、株高、株型、開花習(xí)性、分枝型、植物學(xué)類型9個農(nóng)藝性狀和粗蛋白含量、粗脂肪含量、油酸含量、亞油酸含量4個品質(zhì)性狀。質(zhì)量性狀根據(jù)其表型進(jìn)行賦值(表1), 數(shù)量性狀數(shù)據(jù)采用0.5個標(biāo)準(zhǔn)差為間距進(jìn)行標(biāo)準(zhǔn)化, 分為10級(1級≤-2, 10級 >+2, 中間每級間差0.5,為性狀平均值,為標(biāo)準(zhǔn)差)。

    表1 花生質(zhì)量性狀賦值

    1.2 分組及取樣比例的確定

    根據(jù)中國栽培花生的七大地理種植區(qū)劃(黃河流域花生區(qū)、長江流域花生區(qū)、東南沿?;ㄉ鷧^(qū)、云貴高原花生區(qū)、黃土高原花生區(qū)、東北花生區(qū)和西北花生區(qū))和五大植物學(xué)類型(多粒型、珍珠豆型、龍生型、普通型、中間型)分組。采取常用的平方根法(square root strategy, S法)系統(tǒng)取樣, 即分組取樣量由整個組內(nèi)資源份數(shù)的平方根值占各組平方根之和的比例來決定。

    其中,N為第組的取樣數(shù),n為第組的品種數(shù),n為第組的品種數(shù),為總分組數(shù),為總品種數(shù)。

    1.3 聚類分析

    采用SPSS 22.0統(tǒng)計(jì)分析軟件, 以組為單位采用類平均數(shù)法(UPGMA)進(jìn)行各性狀的系統(tǒng)聚類分析, 類內(nèi)利用Microsoft Excel提供的Randbetween函數(shù)隨機(jī)取樣。

    1.4 骨干種質(zhì)的代表性檢測

    利用SPSS 22.0中的測驗(yàn)、測驗(yàn)分別判斷全部種質(zhì)和骨干種質(zhì)13個性狀的均值、方差、變異系數(shù)及Shannon-Weaver多樣性指數(shù)(′)是否有差異, 卡平方(c2)測驗(yàn)用于檢驗(yàn)兩者的植物學(xué)類型組成、生態(tài)分布及13個性狀的表型分布頻率是否一致。變異系數(shù)(CV) = (標(biāo)準(zhǔn)偏差SD/平均值Mean)×100%, 用于比較全部種質(zhì)和骨干種質(zhì)的離散程度。

    其中,P表示某性狀第級的分布頻率,為總分級數(shù), 用于評價(jià)2個群體的遺傳多樣性水平。極差分析和表型相關(guān)性分析分別用于確定全部種質(zhì)的變異范圍和性狀相關(guān)性是否在骨干種質(zhì)中得到了相應(yīng)保持。表型保留比例(the ratio of phenotypic retention, RPR)用于檢測骨干種質(zhì)中是否保留了全部種質(zhì)足夠的變異。

    其中,M為全部種質(zhì)中某性狀的表型變異數(shù),M為骨干種質(zhì)中該性狀的表型變異數(shù),為性狀總數(shù))。

    1.5 骨干種質(zhì)的確認(rèn)

    利用主成分分析和直方圖比較全部種質(zhì)與骨干種質(zhì)基于主成分的樣品分布圖及13個性狀的分布頻率, 對構(gòu)建的骨干種質(zhì)有效性進(jìn)行確認(rèn)。

    2 結(jié)果與分析

    2.1 骨干種質(zhì)的構(gòu)建、植物學(xué)組成及生態(tài)分布

    全部花生種質(zhì)被分成了26組, 每組種質(zhì)從1到623份不等。每組的取樣量為1~47份。UPGMA法聚類分析后隨機(jī)取樣, 得到242份種質(zhì), 初步評價(jià)后補(bǔ)充極值材料和特殊性狀種質(zhì)17份, 提取9.4%的中國花生地方品種構(gòu)建骨干種質(zhì)。該骨干種質(zhì)包含全部種質(zhì)的植物學(xué)類型, 其中多粒型14份(5.4%)、珍珠豆型85份(32.8%)、龍生型42份(16.2%)、普通型103份(39.8%)、中間型15份(5.8%)。卡方測驗(yàn)表明, 骨干種質(zhì)代表了全部種質(zhì)的植物學(xué)類型組成(c2=1.600,=0.809)(表2)。另外, 骨干種質(zhì)包含了全部種質(zhì)的生態(tài)分布, 其中黃河流域花生區(qū)75份(29.0%)、長江流域花生區(qū)75份(29.0%)、東南沿?;ㄉ鷧^(qū)63份(24.3%)、云貴高原花生區(qū)16份(6.2%)、黃土高原花生區(qū)13份(5.0%)、東北花生區(qū)15份(5.8%)和西北花生區(qū)2份(0.8%), 骨干種質(zhì)代表了全部種質(zhì)的生態(tài)分布(c2=5.232,=0.514), 其中, 西北花生區(qū)由于種質(zhì)較少, 僅有的2份種質(zhì)被全部取樣(表2)。

    2.2 全部種質(zhì)與骨干種質(zhì)的平均值、方差和極差比較

    測驗(yàn)表明, 4個形態(tài)性狀和9個數(shù)量性狀的平均值在全部種質(zhì)與骨干種質(zhì)中無顯著差異, 且骨干種質(zhì)中植物學(xué)類型、百果重和亞油酸含量這3個性狀的平均值都大于全部種質(zhì)的平均值(表3)。測驗(yàn)表明, 除出仁率的方差在全部種質(zhì)與骨干種質(zhì)中差異顯著外, 其余12個性狀的方差均為齊性(表3)。此外, 骨干種質(zhì)大部分性狀的方差高于全部種質(zhì), 表明骨干種質(zhì)遺傳冗余度明顯減小, 變異率更高。極差分析表明, 全部種質(zhì)4個形態(tài)性狀的變異范圍100%保留在骨干種質(zhì)中; 除百仁重(保留范圍64.3%)、出仁率(保留范圍77.2%)外, 其余7個性狀變異范圍的86%~100%保留在骨干種質(zhì)中(表4)。由此可見, 骨干種質(zhì)對全部種質(zhì)的性狀的變異幅度具有良好的代表性。

    2.3 變異系數(shù)和Shannon-Weaver多樣性指數(shù)的比較

    Shannon-Weaver多樣性指數(shù)和變異系數(shù)常用來比較不同樣品的表型特征、等位基因的豐富度和均勻度。從表4可以看出, 全部種質(zhì)和骨干種質(zhì)13個性狀的′和變異系數(shù)是非常相似的。全部種質(zhì)的平均′為1.631±0.153, 變異系數(shù)為29.009±5.964; 骨干種質(zhì)的這2個指標(biāo)分別為1.672±0.154和30.125±5.765, 骨干種質(zhì)略高于全部種質(zhì)。成對雙樣本測驗(yàn)表明, 骨干種質(zhì)′極顯著高于全部種質(zhì)(=0.002), 變異系數(shù)在2個群體中差異不顯著(=0.123), 說明骨干種質(zhì)的樣本足夠大, 且有效去除了總資源中的冗余, 保留了全部種質(zhì)的遺傳多樣性, 變異均勻度顯著提高。

    表2 骨干種質(zhì)和全部種質(zhì)的植物學(xué)類型組成、生態(tài)分布及其卡方測驗(yàn)

    > 0.05表明差異不顯著。> 0.05 means the difference is insignificant.

    表3 全部種質(zhì)與骨干種質(zhì)13個性狀平均值和方差的比較

    > 0.05表示差異不顯著,< 0.05表示差異顯著, NS表示差異不顯著。

    > 0.05 means the difference is insignificant;< 0.05 means the difference is significant; NS: not significant.

    表4 全部種質(zhì)與骨干種質(zhì)極差、變異系數(shù)和遺傳多樣性指數(shù)的比較

    > 0.05表明差異不顯著,< 0.01表明差異極顯著。

    > 0.05 means the difference is insignificant;< 0.01 means the difference is positively significant.

    2.4 分布頻率和表型保留比例的比較

    從表5可以看出, 全部種質(zhì)13個性狀的102個表型分級均包含在骨干種質(zhì)中。對2個群體13個性狀的分布頻率進(jìn)行Chi平方測驗(yàn), 差異均不顯著, 表明2個樣本的性狀分布是一致的, 骨干種質(zhì)可代表全部種質(zhì)的變異。植物學(xué)類型、百仁重和粗脂肪含量的表型保留比例較大, 是補(bǔ)充了一些特殊種質(zhì)和極值材料所致; 其余10個性狀的表型保留比例均比較合適, 表明骨干種質(zhì)保留了全部種質(zhì)豐富的變異, 且豐度更高。

    表5 全部種質(zhì)與骨干種質(zhì)13個性狀的分布頻率和表型保留比例的比較

    > 0.05表明差異不顯著,< 0.01表明差異極顯著。RPR: 表型保留比例。

    > 0.05 means the difference is insignificant;< 0.01 means the difference is positively significant. RPR: the ratio of phenotypic retention.

    2.5 表型相關(guān)分析

    一個具有代表性的骨干種質(zhì)除應(yīng)具有較小的遺傳冗余外, 還應(yīng)保留原群體固有的性狀間的遺傳關(guān)聯(lián)。對13個性狀的表型相關(guān)性分析表明, 64對性狀在全部種質(zhì)中呈極顯著正相關(guān), 其中有59對與骨干種質(zhì)中的相關(guān)性是一致的, 39對在骨干種質(zhì)保持了極顯著正相關(guān)(表6), 因此, 骨干種質(zhì)較好地保持了全部種質(zhì)的表型相關(guān)性。不同性狀間的極顯著相關(guān)性意味著在育種實(shí)踐中, 可以不直接篩選難測量的性狀, 而優(yōu)先篩選與其極顯著相關(guān)的易測量性狀, 從而較快地實(shí)現(xiàn)育種目的。比如, 由于分枝型和油酸含量極顯著負(fù)相關(guān)(全部種質(zhì)中=-0.789**, 骨干種質(zhì)=-0.685**), 那么在雜交后代中要獲得高油酸的單株, 就可以優(yōu)先篩選分枝較多的單株, 直到形成穩(wěn)定的品系后, 再測定其油酸含量。

    表6 在全部種質(zhì)和骨干種質(zhì)中均顯著相關(guān)的性狀

    (續(xù)表6)

    **表示相關(guān)達(dá)到0.01極顯著性水平,*表示相關(guān)達(dá)到0.05顯著水平。

    **indicates significant correlation at≤ 0.01;*indicates significant correlation at≤ 0.05.

    2.6 骨干種質(zhì)的確認(rèn)

    利用主成分分析對所構(gòu)建的骨干種質(zhì)進(jìn)行確認(rèn), 從骨干種質(zhì)和全部種質(zhì)基于第1、第2主成分的樣品分布圖可知(圖1), 全部種質(zhì)中大量的樣品集中在散點(diǎn)圖的左右方并存在較重的相互重疊, 表明這些樣品存在較高的遺傳相似性, 群體的遺傳冗余程度高。9.4%的平方根法取樣后, 骨干種質(zhì)樣品分布的重疊程度得到了顯著降低, 但仍保留了全部種質(zhì)的幾何形狀和特征, 且較多外圍的個體入選到骨干種質(zhì)中, 表明建立的骨干種質(zhì)既去除了全部種質(zhì)的大部分遺傳冗余, 又確保了骨干種質(zhì)的代表性。此外, 繪制了4個植物學(xué)性狀與9個數(shù)量性狀表型分布的直方圖, 可以看出, 植物學(xué)類型、開花習(xí)性、株型和分枝型4個植物學(xué)性狀在全部種質(zhì)和骨干種質(zhì)中的分布非常吻合; 9個數(shù)量性狀在全部種質(zhì)和骨干種質(zhì)中都表現(xiàn)出廣泛的變異, 且大部分變異的分布頻率較為一致, 均基本符合正態(tài)分布(圖2)。因此, 骨干種質(zhì)很好地保留了全部種質(zhì)的遺傳多樣性和群體結(jié)構(gòu), 確保了骨干種質(zhì)的有效性。

    圖1 全部種質(zhì)與9.4%取樣比例骨干種質(zhì)的樣品主成分分布圖

    A: 全部種質(zhì)的樣品分布圖; B: 骨干種質(zhì)的樣品分布圖。

    A: scatter diagram for entire collection; B: scatter diagram for the key germplasm.

    3 討論

    3.1 取樣策略對骨干種質(zhì)構(gòu)建的影響

    要建立質(zhì)量高、代表性強(qiáng)的骨干種質(zhì), 采用的取樣策略很關(guān)鍵。為使構(gòu)建的骨干種質(zhì)能更有效地代表全部種質(zhì)最大的遺傳多樣性, 一般采用2種取樣策略, 即逐步聚類法與分層取樣+系統(tǒng)聚類法。前者多見于基礎(chǔ)數(shù)據(jù)為分子標(biāo)記數(shù)據(jù), 少見于農(nóng)藝數(shù)據(jù)。由于具有最大相似系數(shù)或最近遺傳距離的種質(zhì)能夠聚在一起, 從而篩除相同或相近的種質(zhì)。如劉娟等[18]基于ISSR分子數(shù)據(jù), 采用UPGMA多次聚類抽樣法+位點(diǎn)優(yōu)先取樣策略構(gòu)建了包含31份種質(zhì)的新疆野杏核心種質(zhì)。張春雨等[19]基于SSR分子數(shù)據(jù), 采用UPGMA多次聚類法+位點(diǎn)優(yōu)先取樣策略構(gòu)建了包含25份種質(zhì)的新疆野蘋果核心種質(zhì), 并利用基因多樣度、香農(nóng)指數(shù)、主坐標(biāo)及SRAP數(shù)據(jù)分析確認(rèn)了其代表性。徐益等[20]則是將農(nóng)藝性狀聚類分析和SSR分子標(biāo)記聚類分析相結(jié)合, 遴選出84份黃麻核心種質(zhì), 可以最大限度地代表300份黃麻種質(zhì)資源的遺傳多樣性。劉遵春等[21]基于15個數(shù)量性狀數(shù)據(jù), 采用最短距離法逐步聚類+優(yōu)先取樣法構(gòu)建了包含60份材料的新疆野蘋果核心種質(zhì)。后者常見于以農(nóng)藝數(shù)據(jù)為基礎(chǔ)數(shù)據(jù)的, 以分子數(shù)據(jù)的較少??紤]到作物遺傳多樣性的發(fā)生、發(fā)展層次和分布極不均勻, 一般在建立核心種質(zhì)之前, 將資源歸類, 實(shí)行系統(tǒng)分組取樣, 從而更有效地達(dá)到以最小的重復(fù)代表最大的遺傳多樣性[22]。常見的分組方法基于植物學(xué)分類、系譜起源、地理分布、農(nóng)業(yè)生態(tài)區(qū)、生長習(xí)性等。在系統(tǒng)分組的基礎(chǔ)上, 基于表型或分子標(biāo)記數(shù)據(jù)計(jì)算的遺傳距離進(jìn)行聚類分析, 將材料劃分成不同的類群, 類內(nèi)一般采用隨機(jī)取樣法或位點(diǎn)優(yōu)先取樣法。潘英華等[23]先按照地理分布將普通野生水稻分組, 再利用SSR數(shù)據(jù)進(jìn)行多次UPGMA聚類, 構(gòu)建了包含351份種質(zhì)的廣西普通野生稻核心種質(zhì)。常利芳等[24]先將144份材料分成超甜玉米和普甜玉米2組, 再根據(jù)SSR數(shù)據(jù)進(jìn)行UPGMA聚類分析+位點(diǎn)優(yōu)先取樣, 構(gòu)建了包含33份種質(zhì)的甜玉米核心種質(zhì)。任麗平等[25]先按品種特點(diǎn)、地理來源、生長習(xí)性分組, 再結(jié)合EST-STS標(biāo)記和SSR標(biāo)記分析, 在500余份甘藍(lán)型油菜種質(zhì)中遴選出87份核心種質(zhì)。劉艷陽等[26]根據(jù)地理來源分組, 組內(nèi)比例法聚類抽樣, 結(jié)合SSR標(biāo)記位點(diǎn)優(yōu)先取樣策略進(jìn)行UPGMA逐步聚類, 在5020份芝麻種質(zhì)資源中遴選出501份核心種質(zhì)。本研究按中國花生地方品種的七大地理種植區(qū)劃和五大植物學(xué)類型, 將全部種質(zhì)劃分為26個組, 再根據(jù)13個性狀對組內(nèi)材料進(jìn)行UPGMA聚類分析, 類內(nèi)隨機(jī)取樣, 構(gòu)建的骨干種質(zhì)的均值、方差、極差、變異系數(shù)、香農(nóng)指數(shù)、分布頻率、表型保留比例等與全部種質(zhì)基本無顯著差異, 且很好地保留了全部種質(zhì)的表型相關(guān)性。

    3.2 組內(nèi)取樣量與取樣比例的確定

    作物遺傳結(jié)構(gòu)的差異, 導(dǎo)致在分組聚類的基礎(chǔ)上組內(nèi)取樣量的確定方法也各不相同。常用的有平方根法、比例法、對數(shù)法及多樣性法。劉三才等[27]針對普通小麥核心種質(zhì)進(jìn)行抽樣方法比較表明, 平方根法能提高優(yōu)良類別的頻率, 有利于實(shí)現(xiàn)小麥核心種質(zhì)在育種上的利用。此外, 合理的取樣比例也是構(gòu)建骨干種質(zhì)的重要環(huán)節(jié), 其大小與原始群體的數(shù)量規(guī)模、評價(jià)數(shù)據(jù)類型及物種遺傳結(jié)構(gòu)有關(guān)。目前, 國內(nèi)外所構(gòu)建各類作物核心種質(zhì)的取樣比例基本為5%~30%, 一般在10%左右。Brown[28]提出樣品數(shù)不少于3000時, 以5%~10%的取樣比例就可以代表原始群體70%的變異。Reddy等[29]構(gòu)建鷹嘴豆核心種質(zhì)過程中, 原始種質(zhì)為16,991份, 核心種質(zhì)為1956份, 取樣比例約為10%。Diwan等[30]對美國一年生苜蓿資源的研究表明, 7%是適宜的篩選比例。魏興華等[31]研究450份浙江秈稻地方種質(zhì)的變異, 建立了12.5%的核心種質(zhì)。Zewdie等[32]構(gòu)建的高粱核心種質(zhì)也采用了10%的取樣比例。本研究從2741份花生地方品種中選取了242份骨干種質(zhì), 初步評價(jià)后補(bǔ)充極值材料和特殊性狀種質(zhì)17份, 取樣比例為9.4%, 基本符合核心種質(zhì)的規(guī)模。

    3.3 骨干種質(zhì)的評價(jià)與確認(rèn)

    骨干種質(zhì)的評價(jià)就是檢驗(yàn)其代表性和有效性。李自超等[33]認(rèn)為遺傳多樣性指數(shù)、表型方差、表型分布頻率、變異系數(shù)、表型保留比率等是衡量骨干種質(zhì)的重要參數(shù)。本研究通過檢驗(yàn)、測驗(yàn)、卡方測驗(yàn)對骨干種質(zhì)與全部種質(zhì)比較表明, 多數(shù)性狀的均值、方差、變異系數(shù)、香農(nóng)指數(shù)、分布頻率等無顯著差異; 除百仁重和出仁率外, 大部分性狀的變異范圍均保留在骨干種質(zhì)中。另外, 二者在植物學(xué)組成(c2=1.600,=0.809)、生態(tài)分布(c2=5.232,=0.514)上也有很好的一致性, 表型保留比例和表型相關(guān)性都得到了穩(wěn)定的保持, 可見本研究構(gòu)建的骨干種質(zhì)代表了全部種質(zhì)的遺傳多樣性。出仁率的方差在2個群體間較大, 可能與出仁率的數(shù)據(jù)缺失較多有關(guān)。骨干種質(zhì)保留了百仁重64.3%的變異范圍, 可適當(dāng)補(bǔ)充部分超大仁材料; 出仁率的保留范圍為77.2%, 是由一份出仁率95.20%的種質(zhì)比同級其他種質(zhì)高出約10%所致, 差異在可接受范圍內(nèi)。利用主成分分析和直方圖對骨干種質(zhì)進(jìn)一步確認(rèn), 259份種質(zhì)的遺傳冗余較小, 離散程度和分布特點(diǎn)與全部種質(zhì)一致; 骨干種質(zhì)的性狀變異廣泛, 分布頻率與全部種質(zhì)較符合, 證實(shí)了骨干種質(zhì)的代表性和有效性。

    4 結(jié)論

    基于種植區(qū)劃和植物學(xué)類型分組, 平方根法確定取樣量, 組內(nèi)按表型數(shù)據(jù)進(jìn)行UPGMA聚類分析, 類內(nèi)隨機(jī)取樣, 構(gòu)建了中國花生地方品種的259份骨干種質(zhì), 占全部種質(zhì)的9.4%, 其均值、方差、極差、變異系數(shù)、香農(nóng)指數(shù)、表型保留比例等與全部種質(zhì)無顯著差異, 植物學(xué)類型組成和生態(tài)分布一致, 且保持了全部種質(zhì)的表型分布頻率和表型相關(guān)性。主成分分析和直方圖進(jìn)一步確認(rèn)了骨干種質(zhì)的遺傳多樣性和群體結(jié)構(gòu)。本研究建立的骨干種質(zhì)具有很好的代表性。

    [1] Kochert G, Halward T, Branch W D, Simpson C E. RFLP variability in peanut () cultivars and wild species., 1991, 81: 565–570.

    [2] 姜慧芳, 任小平, 廖伯壽, 黃家權(quán), 陳本銀. 中國花生核心種質(zhì)的建立. 武漢植物學(xué)研究, 2007, 25: 289–293. Jiang H F, Ren X P, Liao B S, Huang J Q, Chen B Y. Establishment of peanut core collection in China., 2007, 25: 289–293 (in Chinese with English abstract).

    [3] 沈一, 鄂志國, 劉永惠, 陳志德. 中國花生品種及其系譜數(shù)據(jù)庫的構(gòu)建. 中國油料作物學(xué)報(bào), 2015, 37: 571–575. Shen Y, E Z G, Liu Y H, Chen Z D. Database construction of Chinese peanut varieties and their genealogy., 2015, 37: 571–575 (in Chinese with English abstract).

    [4] 姜慧芳, 段乃雄. 花生種質(zhì)資源在育種中的利用. 中國種業(yè), 1998, (2): 24–25. Jiang H F, Duan N X. The application of germplasm resources in peanut breeding., 1998, (2): 24–25 (in Chinese with English abstract).

    [5] Frankel O H, Brown A H D. Plant genetic resources today: a critical appraisal. In: Holden J H W, Williams J T, eds. Crop Genetic Resources: Conservation and Evaluation. London: George Allen and Unwin, 1984. pp 249–257.

    [6] 賈繼增, 高麗鋒, 趙光耀, 周文斌, 張衛(wèi)健. 作物基因組學(xué)與作物科學(xué)革命. 中國農(nóng)業(yè)科學(xué), 2015, 48: 3316–3332. Jia J Z, Gao L F, Zhao G Y, Zhou W B, Zhang W J. Crop genomics and crop science revolutions., 2015, 48: 3316–3332 (in Chinese with English abstract).

    [7] Holbrook C C, Anderson W F, Pittman R N. Selection of a core collection from the U.S. germplasm collection of peanut., 1993, 33: 859–861.

    [8] Holbrook C C, Dong W B. Development and evaluation of a mini core collection for the U.S. peanut germplasm collection., 2005, 45: 1540–1544.

    [9] Upadhyaya H D, Bramel P J, Ortiz R, Singh S. Developing a mini core of peanut for utilization of genetic resources., 2002, 42: 2150–2156.

    [10] Upadhyaya H D, Ortiz R, Bramel P J, Sube S. Development of a groundnut core collection using taxonomical, geographical and morphological descriptors., 2003, 50: 139–148.

    [11] Chamberlin K D C, Hassan A M, Mark E P. Evaluation of the U.S. peanut mini core collection using a molecular marker for resistance toJagger., 2010, 172: 109–115.

    [12] Mukri G, Hajisaheb L, Nadaf R S B, Gowda M V C, Upadhyaya H D, Sujay V. Phenotypic and molecular dissection of ICRISAT mini core collection of peanut (L.) for high oleic acid., 2012, 131: 418–422.

    [13] Upadhyaya H D, Dwivedi S L, Vadez V, Hamidou F, Singh S, Varshney R K, Liao B. Multiple resistant and nutritionally dense germplasm identified from mini core collection in peanut., 2013, 54: 679–693.

    [14] Sudini H, Upadhyaya H, Reddy S V, Mangala U N, Kumar K V. Resistance to late leaf spot and rust diseases in ICRISAT’s mini core collection of peanut (L.)., 2015, 44: 557–566.

    [15] Hovav R, Badani H, Ginzberg I, Hovav R, Badani H, Ginzberg I, Chedvat I, Brand Y, Galili S. Evaluation of a peanut collection for shell-colour traits in two diverse soil types., 2012, 131: 148–154.

    [16] 黃莉, 任小平, 張曉杰, 陳玉寧, 姜慧芳. ICRISAT花生微核心種質(zhì)農(nóng)藝性狀和黃曲霉抗性關(guān)聯(lián)分析. 作物學(xué)報(bào), 2012, 38: 935–946. Huang L, Ren X P, Zhang X J, Chen Y N, Jiang H F. Association analysis of agronomic traits and resistance toin the ICRISAT peanut mini-core collection., 2012, 38: 935–946 (in Chinese with English abstract).

    [17] 任小平, 廖伯壽, 張曉杰, 雷永, 黃家權(quán), 晏立英, 陳玉中, 姜慧芳. 中國花生核心種質(zhì)中高油酸材料的分布和遺傳多樣性. 植物遺傳資源學(xué)報(bào), 2011, 12: 513–518. Ren X P, Liao B S, Zhang X J, Lei Y, Huang J Q, Yan L Y, Chen Y Z, Jiang H F. Distributing and genetic diversity of high oleic acid germplasm in peanut (L.) core collection of China., 2011, 12: 513–518 (in Chinese with English abstract).

    [18] 劉娟, 廖康, 趙世榮, 曹倩, 孫琪, 劉歡. 利用ISSR分子標(biāo)記構(gòu)建新疆野杏核心種質(zhì)資源. 中國農(nóng)業(yè)科學(xué), 2015, 48: 2017–2028. Liu J, Liao K, Zhao S R, Cao Q, Sun Q, Liu H. The core collection construction of Xinjiang wild apricot based on ISSR molecular markers., 2015, 48: 2017–2028 (in Chinese with English abstract).

    [19] 張春雨, 陳學(xué)森, 張艷敏, 苑兆和, 劉遵春, 王延齡, 林群. 采用分子標(biāo)記構(gòu)建新疆野蘋果核心種質(zhì)的方法. 中國農(nóng)業(yè)科學(xué), 2009, 42: 597–604. Zhang C Y, Chen X S, Zhang Y M, Yuan Z H, Liu Z C, Wang Y L, Lin Q. A method for constructing core collection of Malussieve-rsii using molecular markers., 2009, 42: 597–604 (in Chinese with English abstract).

    [20] 徐益, 張列梅, 郭艷春, 祁建民, 張力嵐, 方平平, 張立武. 黃麻核心種質(zhì)的遴選. 作物學(xué)報(bào), 2019, 45: 1672–1681. Xu Y, Zhang L M, Guo Y C, Qi J M, Zhang L L, Fang P P, Zhang L W. Core collection screening of a germplasm population in jute (spp.)., 2019, 45: 1672–1681 (in Chinese with English abstract).

    [21] 劉遵春, 張春雨, 張艷敏, 張小燕, 吳傳金, 王海波, 石俊, 陳學(xué)森. 利用數(shù)量性狀構(gòu)建新疆野蘋果核心種質(zhì)的方法. 中國農(nóng)業(yè)科學(xué), 2010, 43: 358–370. Liu Z C, Zhang C Y, Zhang Y M, Zhang X Y, Wu C J, Wang H B, Shi J, Chen X S. Study on method of constructing core collection of Malussieversii based on quantitative traits., 2010, 43: 358–370 (in Chinese with English abstract).

    [22] Zeuli P L S, Qualset C O. Evaluation of five strategies for obtaining a core subset from a large genetic resource collection of durum wheat., 1993, 87: 295–304.

    [23] 潘英華, 徐志健, 梁云濤. 廣西普通野生稻群體結(jié)構(gòu)解析與核心種質(zhì)構(gòu)建. 植物遺傳資源學(xué)報(bào), 2018, 19: 498–509. Pan Y H, Xu Z J, Liang Y T. Genetic structure and core collection of common wild rice (Griff.) in Guangxi., 2018, 19: 498–509 (in Chinese with English abstract).

    [24] 常利芳, 白建榮, 李銳, 張叢卓, 張效梅, 楊瑞娟. 基于SSR標(biāo)記構(gòu)建甜玉米群體的核心種質(zhì). 玉米科學(xué), 2018, 26(3): 40–49.Chang L F, Bai J R, Li R, Zhang C Z, Zhang X M, Yang R J. Construction of a core collection of sweet corn populations based on SSR markers., 2018, 26(3): 40–49 (in Chinese with English abstract).

    [25] 任麗平, 倪西源, 黃吉祥, 雷偉俠, 曹明富, 趙堅(jiān)義. 甘藍(lán)型油菜一個代表性核心種質(zhì)的遴選. 中國農(nóng)業(yè)科學(xué), 2008, 41: 3521–3531. Ren L P, Ni X Y, Huang J X, Lei W X, Cao M F, Zhao J Y. Core collection of a representative germplasm population in., 2008, 41: 3521–3531 (in Chinese with English abstract).

    [26] 劉艷陽, 梅鴻獻(xiàn), 杜振偉, 武軻, 鄭永戰(zhàn), 崔向華, 鄭磊. 基于表型和SSR分子標(biāo)記構(gòu)建芝麻核心種質(zhì). 中國農(nóng)業(yè)科學(xué), 2017, 50: 2433–2441. Liu Y Y, Mei H X, Du Z W, Wu K, Zheng Y Z, Cui X H, Zheng L. Construction of core collection of sesame based on phenotype and molecular markers., 2017, 50: 2433–2441 (in Chinese with English abstract).

    [27] 劉三才, 曹永生, 鄭殿升, 劉春華, 陳夢英. 普通小麥核心種質(zhì)抽樣方法的比較. 麥類作物學(xué)報(bào), 2001, 21(2): 42–45. Liu S C, Cao Y S, Zheng D S, Liu C H, Chen M Y. Comparing of strategies for developing core collection from common wheat., 2001, 21(2): 42–45 (in Chinese with English abstract).

    [28] Brown A H D. Core collections: a practical approach to genetic resources management., 1989, 31: 818–824.

    [29] Reddy L J, Upadhyaya H D, Gowda C L L, Sube S. Development of core collection in pigeonpea [(L.) Millspaugh] using geographic and qualitative morphological descriptors., 2005, 52: 1049–1056.

    [30] Diwan N, McIntosh M S, Bauchan G R. Methods of developing a core collection of annual Medicago species., 1995, 90: 755–761.

    [31] 魏興華, 顏啟傳, 應(yīng)存山, 張麗華, 章林平. 建立浙江地方秈型稻種資源的核心樣品的研究. 中國水稻科學(xué), 1999, 13(2): 81–85.Wei X H, Yan Q C, Ying C S, Zhang L H, Zhang L P. A core collection of Zhejiang traditionalrice germplasm., 1999, 13(2): 81–85 (in Chinese with English abstract).

    [32] Zewdie Y, Tong N, Bosland P. Establishing a core collection of Capsicum, using a cluster analysis with enlightened selection of accessions., 2004, 51: 147–151.

    [33] 李自超, 張洪亮, 曹永生, 裘宗恩, 魏興華, 湯圣祥, 余萍, 王象坤. 中國地方稻種資源初級核心種質(zhì)取樣策略研究. 作物學(xué)報(bào), 2003, 29: 20–24.Li Z C, Zhang H L, Cao Y S, Qiu Z E, Wei X H, Tang S X, Yu P, Wang X K. Studies on the sampling strategy for initial core collection of Chinese ingenious rice., 2003, 29: 20–24 (in Chinese with English abstract).

    Developing the key germplasm of Chinese peanut landraces based on phenotypic traits

    YAN Cai-Xia1, WANG Juan1, ZHANG Hao1, LI Chun-Juan1, SONG Xiu-Xia2, SUN Quan-Xi1, YUAN Cui-Ling1, ZHAO Xiao-Bo1, and SHAN Shi-Hua1,*

    1Shandong Peanut Research Institute, Qingdao 266100, Shandong, China;2Agricultural and Rural Bureau of Mudan District, Heze 274000, Shandong, China

    Chinese peanut landraces are important parent resources in peanut breeding due to their abundant genetic diversity. In this study, a total of 2741 original accessions from peanut seed bank were divided into 26 groups based on their botanical variety and ecological distribution. The key accessions were established based on the analysis of 13 phenotypic traits by the square root strategy, UPGMA clustering within groups and random sampling in individual clusters, and evaluated by-test,-test, Chi-squared test, ranging, the ratio of phenotypic retention, and phenotypic correlation analyses. Finally, the principal components analysis (PCA) and the histogram analysis were used to re-confirm the key germplasm. The total of 259 as a key germplasm was selected, accounting for 9.4% of total accessions, which included 14 of var., 85 of var., 42 of var., 103 of var.and 15 of irregular type. There were no significant differences (< 0.05) in means, variance, coefficient of variation, and Shannon-weaver diversity index for 13 phenotypic traits between key germplasm and entire collection. The key germplasm preserved the distribution range, the ratio of phenotypic retention and the phenotypic correlation of primary collection, with similar composition of botanical variety and ecological distribution. PCA and the histogram confirmed the homogeneity of genetic structure and distribution frequency between two collections. Thus, this key germplasm can represent the genetic variability and population structure of entire collection, and enhance innovation of peanut genetic resources and exploitation of elite alleles.

    peanut; landrace; phenotypic traits; key germplasm; representative evaluation

    2019-07-17;

    2019-12-26;

    2020-01-17.

    10.3724/SP.J.1006.2020.94101

    單世華, E-mail: shansh1971@163.com, Tel: 0532-87629307

    E-mail: cxyan335@sina.com, Tel: 0532-87626756

    本研究由泰山學(xué)者特聘專家項(xiàng)目(ts201712080), 中央引導(dǎo)地方科技發(fā)展專項(xiàng)資金, 農(nóng)業(yè)科研杰出人才及其創(chuàng)新團(tuán)隊(duì)培養(yǎng)(13190194), 青島市民生科技計(jì)劃項(xiàng)目(17-3-3-49-nsh), 山東省農(nóng)業(yè)良種工程項(xiàng)目(2017LZN033, 2017LZGC003), 山東省農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系項(xiàng)目(SDAIT-04-02)和山東省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)科技創(chuàng)新工程項(xiàng)目(CXGC2016A01)資助。

    This study was supported by the Taishan Scholars Project (ts201712080), the Central Guidance for Local Science and Technology, the Outstanding Talents and Innovation Team in Agricultural Research (13190194), the Qingdao Science and Technology Plan for the Public Benefit (17-3-3-49-nsh), the Fine Breeding Project of Shandong Province (2017LZN033, 2017LZGC003), the Shandong Agriculture Research System (SDAIT-04-02), and the Agricultural Science and Technological Innovation Project of Shandong Academy of Agricultural Science (CXGC2016A01).

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

    猜你喜歡
    骨干表型種質(zhì)
    華南地區(qū)最大農(nóng)作物種質(zhì)資源保護(hù)庫建成
    核心研發(fā)骨干均16年以上!創(chuàng)美克在產(chǎn)品研發(fā)上再發(fā)力
    亞麻抗白粉病種質(zhì)資源的鑒定與篩選
    建蘭、寒蘭花表型分析
    骨干風(fēng)采展示
    貴州玉米種質(zhì)資源遺傳多樣性及核心種質(zhì)庫構(gòu)建
    紅錐種質(zhì)早期生長表現(xiàn)
    GABABR2基因遺傳變異與肥胖及代謝相關(guān)表型的關(guān)系
    慢性乙型肝炎患者HBV基因表型與血清學(xué)測定的臨床意義
    72例老年急性白血病免疫表型分析
    一级,二级,三级黄色视频| 欧美精品高潮呻吟av久久| 午夜91福利影院| 免费观看av网站的网址| 日本91视频免费播放| 性少妇av在线| 日韩中文字幕视频在线看片| 亚洲天堂av无毛| 亚洲综合色网址| 啦啦啦中文免费视频观看日本| 制服人妻中文乱码| 春色校园在线视频观看| 国产免费视频播放在线视频| 青青草视频在线视频观看| 岛国毛片在线播放| 国产在线一区二区三区精| 黄片无遮挡物在线观看| 国产黄色视频一区二区在线观看| 免费黄网站久久成人精品| 免费少妇av软件| 国产黄色视频一区二区在线观看| 男男h啪啪无遮挡| 久久精品亚洲av国产电影网| 香蕉丝袜av| 伦理电影免费视频| 亚洲精华国产精华液的使用体验| 日韩一卡2卡3卡4卡2021年| 午夜福利在线免费观看网站| 国产精品 欧美亚洲| 视频在线观看一区二区三区| 岛国毛片在线播放| 欧美激情极品国产一区二区三区| 高清黄色对白视频在线免费看| 午夜福利在线观看免费完整高清在| 肉色欧美久久久久久久蜜桃| 永久网站在线| 亚洲人成电影观看| 国产熟女午夜一区二区三区| 亚洲视频免费观看视频| 在线 av 中文字幕| 侵犯人妻中文字幕一二三四区| 美国免费a级毛片| 男人操女人黄网站| 在现免费观看毛片| 日本欧美视频一区| 你懂的网址亚洲精品在线观看| 青青草视频在线视频观看| 久久久a久久爽久久v久久| av福利片在线| 国产免费福利视频在线观看| 水蜜桃什么品种好| 考比视频在线观看| 国产熟女午夜一区二区三区| 国产亚洲av片在线观看秒播厂| 中文天堂在线官网| 中文字幕av电影在线播放| 男女下面插进去视频免费观看| 深夜精品福利| 久久国内精品自在自线图片| 久久国内精品自在自线图片| 日韩一区二区三区影片| 女人被躁到高潮嗷嗷叫费观| 日本欧美国产在线视频| 99re6热这里在线精品视频| 色哟哟·www| 天美传媒精品一区二区| av有码第一页| 国产成人欧美| 亚洲精品久久成人aⅴ小说| 欧美黄色片欧美黄色片| 伦理电影大哥的女人| 日韩一区二区视频免费看| 亚洲三区欧美一区| 97人妻天天添夜夜摸| 九色亚洲精品在线播放| 久久久精品94久久精品| 亚洲一区二区三区欧美精品| 只有这里有精品99| www.精华液| 一本—道久久a久久精品蜜桃钙片| 国产成人一区二区在线| 国产视频首页在线观看| 欧美日韩视频精品一区| 26uuu在线亚洲综合色| 国产黄频视频在线观看| 亚洲欧美一区二区三区国产| 亚洲国产av影院在线观看| 啦啦啦在线免费观看视频4| 中文欧美无线码| 美女国产高潮福利片在线看| 又粗又硬又长又爽又黄的视频| 国产一区二区三区av在线| 观看av在线不卡| 久久久久精品久久久久真实原创| 香蕉国产在线看| 久久精品国产a三级三级三级| 日韩精品免费视频一区二区三区| 一边摸一边做爽爽视频免费| 久热久热在线精品观看| 免费黄频网站在线观看国产| 欧美变态另类bdsm刘玥| 伦精品一区二区三区| 日韩欧美一区视频在线观看| 日韩一区二区视频免费看| 九草在线视频观看| 人体艺术视频欧美日本| 久久精品人人爽人人爽视色| 欧美精品一区二区大全| 中文天堂在线官网| 99热国产这里只有精品6| 一边亲一边摸免费视频| 久久久久人妻精品一区果冻| 日日爽夜夜爽网站| 欧美精品亚洲一区二区| 下体分泌物呈黄色| 成人二区视频| xxx大片免费视频| 国产熟女午夜一区二区三区| 美女福利国产在线| 另类精品久久| 亚洲欧美日韩另类电影网站| 久久人人爽av亚洲精品天堂| 日韩在线高清观看一区二区三区| 18禁动态无遮挡网站| 在线免费观看不下载黄p国产| 嫩草影院入口| 少妇人妻久久综合中文| 国产黄频视频在线观看| 美女脱内裤让男人舔精品视频| 久久久久久久亚洲中文字幕| 国产在线免费精品| 波多野结衣一区麻豆| 久久久精品94久久精品| 亚洲成人手机| 免费播放大片免费观看视频在线观看| 色哟哟·www| av.在线天堂| 电影成人av| 日韩,欧美,国产一区二区三区| 午夜日韩欧美国产| 伦理电影大哥的女人| 日韩制服骚丝袜av| 日韩熟女老妇一区二区性免费视频| 精品酒店卫生间| 免费高清在线观看日韩| av网站免费在线观看视频| 黄色一级大片看看| 你懂的网址亚洲精品在线观看| 国产在线免费精品| 美女高潮到喷水免费观看| 在线天堂最新版资源| 免费在线观看黄色视频的| 免费av中文字幕在线| av不卡在线播放| 香蕉精品网在线| 久久久国产精品麻豆| 国产男人的电影天堂91| 国产成人欧美| 亚洲一码二码三码区别大吗| 男女国产视频网站| 精品99又大又爽又粗少妇毛片| 天堂中文最新版在线下载| 国产精品三级大全| 日日啪夜夜爽| 在线观看免费日韩欧美大片| 国产精品嫩草影院av在线观看| 亚洲精华国产精华液的使用体验| 桃花免费在线播放| 国产日韩欧美视频二区| 日韩熟女老妇一区二区性免费视频| 777米奇影视久久| 亚洲精品成人av观看孕妇| 日韩,欧美,国产一区二区三区| 中国三级夫妇交换| 另类精品久久| 国产女主播在线喷水免费视频网站| 欧美精品av麻豆av| 亚洲国产成人一精品久久久| 91午夜精品亚洲一区二区三区| 黄色怎么调成土黄色| 97在线人人人人妻| 中文欧美无线码| 亚洲欧美清纯卡通| 免费在线观看完整版高清| 日本wwww免费看| 水蜜桃什么品种好| 国产av码专区亚洲av| 自拍欧美九色日韩亚洲蝌蚪91| 久久久精品区二区三区| 国产成人精品无人区| 黑人欧美特级aaaaaa片| 亚洲欧美色中文字幕在线| 最新的欧美精品一区二区| 精品少妇内射三级| 在线 av 中文字幕| 国产成人精品在线电影| 夫妻午夜视频| 高清在线视频一区二区三区| 中文欧美无线码| 婷婷色综合www| 少妇的逼水好多| 日韩伦理黄色片| 亚洲国产精品成人久久小说| 午夜免费男女啪啪视频观看| 可以免费在线观看a视频的电影网站 | 在线观看美女被高潮喷水网站| 久久精品国产a三级三级三级| a级毛片在线看网站| 午夜免费鲁丝| 日韩精品免费视频一区二区三区| 26uuu在线亚洲综合色| 丝袜脚勾引网站| 久久久久国产网址| 亚洲精品国产色婷婷电影| 国产欧美日韩一区二区三区在线| 丝袜人妻中文字幕| 超碰成人久久| 成年人免费黄色播放视频| 日韩av免费高清视频| 999久久久国产精品视频| 天天操日日干夜夜撸| 啦啦啦在线免费观看视频4| 秋霞伦理黄片| 久久久久精品久久久久真实原创| 丰满饥渴人妻一区二区三| 久久国产精品大桥未久av| 青春草亚洲视频在线观看| www.自偷自拍.com| 日本-黄色视频高清免费观看| 可以免费在线观看a视频的电影网站 | 亚洲欧美精品自产自拍| 男人操女人黄网站| 深夜精品福利| 青春草国产在线视频| 男女午夜视频在线观看| 免费大片黄手机在线观看| 欧美精品一区二区免费开放| 少妇人妻久久综合中文| 亚洲成人手机| 丝瓜视频免费看黄片| 高清视频免费观看一区二区| 波野结衣二区三区在线| 哪个播放器可以免费观看大片| 亚洲综合色网址| 国产精品久久久久成人av| 丝袜喷水一区| 免费女性裸体啪啪无遮挡网站| 18禁裸乳无遮挡动漫免费视频| av国产久精品久网站免费入址| 国产精品一二三区在线看| 午夜日韩欧美国产| 久久久久久久亚洲中文字幕| 伦理电影大哥的女人| 午夜日本视频在线| 精品久久久久久电影网| 嫩草影院入口| 精品99又大又爽又粗少妇毛片| 成人免费观看视频高清| 只有这里有精品99| 久久国产亚洲av麻豆专区| 国产无遮挡羞羞视频在线观看| 高清视频免费观看一区二区| 婷婷色综合大香蕉| 啦啦啦在线免费观看视频4| 亚洲av在线观看美女高潮| 人体艺术视频欧美日本| 国产成人精品婷婷| av女优亚洲男人天堂| 国产日韩欧美亚洲二区| 国产熟女欧美一区二区| 国产精品国产三级专区第一集| 亚洲欧洲日产国产| 97人妻天天添夜夜摸| 老女人水多毛片| 波野结衣二区三区在线| 久久这里只有精品19| 久久国内精品自在自线图片| 国产精品一二三区在线看| 黄色配什么色好看| 免费观看在线日韩| 91aial.com中文字幕在线观看| 色视频在线一区二区三区| 超碰97精品在线观看| 不卡视频在线观看欧美| 免费高清在线观看视频在线观看| 综合色丁香网| 国产精品国产三级国产专区5o| 免费看不卡的av| 欧美日韩成人在线一区二区| 国产片内射在线| 最近的中文字幕免费完整| 考比视频在线观看| 婷婷色综合大香蕉| 我要看黄色一级片免费的| 日韩电影二区| 久久女婷五月综合色啪小说| 中文精品一卡2卡3卡4更新| 成年人午夜在线观看视频| 亚洲精品久久成人aⅴ小说| 日本色播在线视频| 9191精品国产免费久久| 国产精品二区激情视频| 各种免费的搞黄视频| 亚洲精品中文字幕在线视频| 亚洲av日韩在线播放| 99久久综合免费| 欧美老熟妇乱子伦牲交| av网站在线播放免费| 99久久人妻综合| 如日韩欧美国产精品一区二区三区| 日韩 亚洲 欧美在线| 久久精品国产亚洲av高清一级| 亚洲婷婷狠狠爱综合网| 九草在线视频观看| 啦啦啦视频在线资源免费观看| 国产成人一区二区在线| 青青草视频在线视频观看| 午夜福利一区二区在线看| 国产成人精品婷婷| 国产淫语在线视频| 十八禁网站网址无遮挡| 在线观看国产h片| 成人午夜精彩视频在线观看| 高清av免费在线| 久久精品国产亚洲av天美| 男人操女人黄网站| 成年人午夜在线观看视频| 一级毛片黄色毛片免费观看视频| 高清欧美精品videossex| 欧美激情极品国产一区二区三区| 高清不卡的av网站| 国产av一区二区精品久久| 国产97色在线日韩免费| tube8黄色片| 桃花免费在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 观看av在线不卡| 精品国产超薄肉色丝袜足j| 最近最新中文字幕大全免费视频 | 国产精品一国产av| 女的被弄到高潮叫床怎么办| 中文字幕人妻熟女乱码| 女的被弄到高潮叫床怎么办| 日韩精品有码人妻一区| 亚洲精品一区蜜桃| 亚洲av成人精品一二三区| 18禁观看日本| 日韩制服骚丝袜av| 精品亚洲成国产av| 大码成人一级视频| 国产又爽黄色视频| 欧美中文综合在线视频| 亚洲视频免费观看视频| 人人澡人人妻人| 人人妻人人爽人人添夜夜欢视频| 久久影院123| 午夜免费男女啪啪视频观看| 久久热在线av| 亚洲久久久国产精品| 男女高潮啪啪啪动态图| 国产伦理片在线播放av一区| 亚洲精品久久午夜乱码| 免费看不卡的av| 久久ye,这里只有精品| 久久毛片免费看一区二区三区| 丰满乱子伦码专区| 日韩在线高清观看一区二区三区| 在线观看美女被高潮喷水网站| 深夜精品福利| 黄色毛片三级朝国网站| 国产免费现黄频在线看| 精品卡一卡二卡四卡免费| 五月开心婷婷网| 人人妻人人添人人爽欧美一区卜| 大香蕉久久网| 午夜日韩欧美国产| 免费观看在线日韩| 人妻少妇偷人精品九色| 婷婷色麻豆天堂久久| 午夜老司机福利剧场| 人人妻人人爽人人添夜夜欢视频| 国产激情久久老熟女| av在线播放精品| 国产不卡av网站在线观看| 国产成人91sexporn| 狂野欧美激情性bbbbbb| 国产亚洲av片在线观看秒播厂| 日韩一区二区三区影片| av电影中文网址| 亚洲,欧美,日韩| 最近最新中文字幕免费大全7| 五月天丁香电影| 亚洲,欧美,日韩| 欧美日韩亚洲国产一区二区在线观看 | 婷婷成人精品国产| av.在线天堂| 激情五月婷婷亚洲| 丝袜在线中文字幕| 男的添女的下面高潮视频| 一本大道久久a久久精品| 成年人免费黄色播放视频| 精品一区在线观看国产| 欧美亚洲日本最大视频资源| 欧美最新免费一区二区三区| 日韩成人av中文字幕在线观看| 日韩视频在线欧美| 美女主播在线视频| av在线播放精品| 高清视频免费观看一区二区| 多毛熟女@视频| av在线老鸭窝| 婷婷色综合大香蕉| 丝袜脚勾引网站| 亚洲成人av在线免费| 91在线精品国自产拍蜜月| 久久狼人影院| xxx大片免费视频| 性色av一级| 人体艺术视频欧美日本| 在线观看国产h片| 日韩一区二区视频免费看| av免费观看日本| 亚洲欧美清纯卡通| 黄色配什么色好看| 黑人猛操日本美女一级片| 国产乱人偷精品视频| freevideosex欧美| 一级毛片我不卡| 亚洲在久久综合| 久久ye,这里只有精品| 欧美在线黄色| 成人国产麻豆网| 一边摸一边做爽爽视频免费| 国产女主播在线喷水免费视频网站| videossex国产| 欧美激情极品国产一区二区三区| 两个人看的免费小视频| 午夜福利乱码中文字幕| 日韩大片免费观看网站| 中文字幕制服av| 超色免费av| 久久精品夜色国产| 狠狠婷婷综合久久久久久88av| 久久午夜福利片| 国产成人一区二区在线| 多毛熟女@视频| 欧美97在线视频| 26uuu在线亚洲综合色| 亚洲成人一二三区av| 成人毛片a级毛片在线播放| 色播在线永久视频| 久久精品人人爽人人爽视色| 91精品国产国语对白视频| 国产一区二区三区综合在线观看| 在线观看国产h片| 亚洲欧美成人精品一区二区| 成人漫画全彩无遮挡| 久久av网站| 欧美国产精品一级二级三级| 日韩电影二区| 日韩成人av中文字幕在线观看| 夫妻性生交免费视频一级片| 欧美日韩精品网址| 观看av在线不卡| 久久久国产欧美日韩av| 国产成人精品婷婷| 嫩草影院入口| 国产一区二区三区综合在线观看| 国产日韩一区二区三区精品不卡| 80岁老熟妇乱子伦牲交| 国产亚洲欧美精品永久| 国产日韩欧美在线精品| 哪个播放器可以免费观看大片| 黄网站色视频无遮挡免费观看| 亚洲精品视频女| 777久久人妻少妇嫩草av网站| 最近最新中文字幕大全免费视频 | 国产黄频视频在线观看| 五月天丁香电影| 满18在线观看网站| 色播在线永久视频| videosex国产| 久久韩国三级中文字幕| 美女高潮到喷水免费观看| 国产成人免费无遮挡视频| 黄色一级大片看看| 啦啦啦视频在线资源免费观看| 日韩,欧美,国产一区二区三区| 校园人妻丝袜中文字幕| 高清视频免费观看一区二区| 天堂8中文在线网| 亚洲视频免费观看视频| 亚洲av电影在线观看一区二区三区| 日韩在线高清观看一区二区三区| 制服人妻中文乱码| 18禁国产床啪视频网站| 久久综合国产亚洲精品| 国产日韩欧美在线精品| 国产男人的电影天堂91| 男女无遮挡免费网站观看| 国产亚洲午夜精品一区二区久久| 赤兔流量卡办理| av一本久久久久| 老司机影院成人| 高清欧美精品videossex| 国产精品不卡视频一区二区| 9色porny在线观看| 考比视频在线观看| 亚洲国产精品一区三区| 深夜精品福利| 国产高清国产精品国产三级| 国产免费又黄又爽又色| 大香蕉久久成人网| 久久久久久人人人人人| 日韩制服骚丝袜av| 一个人免费看片子| 少妇人妻久久综合中文| 精品亚洲成a人片在线观看| xxx大片免费视频| 在线观看人妻少妇| 久久久久视频综合| 久久国产精品男人的天堂亚洲| 国产黄色免费在线视频| 超碰成人久久| 免费观看在线日韩| 2021少妇久久久久久久久久久| 国产福利在线免费观看视频| 侵犯人妻中文字幕一二三四区| 欧美日韩精品成人综合77777| 久久久久久人人人人人| 桃花免费在线播放| 精品亚洲成国产av| 色吧在线观看| 人妻系列 视频| 久久精品久久久久久噜噜老黄| 国产黄色视频一区二区在线观看| 国产免费福利视频在线观看| 水蜜桃什么品种好| 视频在线观看一区二区三区| 毛片一级片免费看久久久久| 90打野战视频偷拍视频| 一级毛片我不卡| 午夜福利在线观看免费完整高清在| 亚洲国产色片| 免费人妻精品一区二区三区视频| 久久国内精品自在自线图片| 亚洲精品aⅴ在线观看| 99久国产av精品国产电影| 亚洲av日韩在线播放| 午夜激情av网站| 视频在线观看一区二区三区| 日韩中文字幕欧美一区二区 | 欧美日韩视频高清一区二区三区二| 考比视频在线观看| 中国国产av一级| 十八禁网站网址无遮挡| 美女福利国产在线| 欧美日韩一区二区视频在线观看视频在线| 成人18禁高潮啪啪吃奶动态图| 国产欧美亚洲国产| 国产免费视频播放在线视频| 国产av码专区亚洲av| 亚洲国产精品成人久久小说| freevideosex欧美| 飞空精品影院首页| 日韩一卡2卡3卡4卡2021年| 国产精品麻豆人妻色哟哟久久| 免费久久久久久久精品成人欧美视频| 亚洲视频免费观看视频| 一区二区三区四区激情视频| 九草在线视频观看| 一本大道久久a久久精品| av网站免费在线观看视频| 日韩中文字幕欧美一区二区 | 亚洲视频免费观看视频| 日本爱情动作片www.在线观看| 夫妻性生交免费视频一级片| 国产精品嫩草影院av在线观看| 亚洲av国产av综合av卡| 秋霞在线观看毛片| 99香蕉大伊视频| 久热久热在线精品观看| 亚洲综合色惰| 一级爰片在线观看| 亚洲av电影在线观看一区二区三区| 黄色一级大片看看| 免费在线观看视频国产中文字幕亚洲 | 久久精品国产亚洲av高清一级| 久久久精品区二区三区| 日日撸夜夜添| 蜜桃在线观看..| 中国国产av一级| videosex国产| 五月天丁香电影| 男人操女人黄网站| 91成人精品电影| 久久久精品94久久精品| 欧美人与善性xxx| av卡一久久| 在线看a的网站| 一二三四在线观看免费中文在| 亚洲色图 男人天堂 中文字幕| videos熟女内射| 国产人伦9x9x在线观看 | 搡女人真爽免费视频火全软件| 九九爱精品视频在线观看| 亚洲av.av天堂| 精品一区二区三卡| 久久ye,这里只有精品| 国产在线视频一区二区| 国产精品久久久久久av不卡| 99精国产麻豆久久婷婷| 成年人午夜在线观看视频| 国产精品偷伦视频观看了| 亚洲伊人色综图|