中國花生品種更替引發(fā)的SSR位點遺傳多樣性變化趨勢分析

王 瑾，李玉榮，程增書，陳四龍，宋亞輝，張朋娟

（河北省農(nóng)林科學院糧油作物研究所河北省作物遺傳育種實驗室，石家莊 050035）

中國花生品種更替引發(fā)的SSR位點遺傳多樣性變化趨勢分析

王 瑾，李玉榮*，程增書，陳四龍，宋亞輝，張朋娟

（河北省農(nóng)林科學院糧油作物研究所河北省作物遺傳育種實驗室，石家莊 050035）

【目的】旨在了解自20世紀50年代以來我國花生品種更替所引發(fā)的SSR位點遺傳多樣性變化，以期為花生育種提供參考。【方法】選用154對SSR引物對68個大面積推廣種植品種進行檢測。【結果】共獲得173個位點，檢測到872個等位變異，平均為5.04個等位變異/位點，遺傳多樣性指數(shù)的變化范圍為0.014～0.881，平均0.477。20條染色體中，b07遺傳多樣性最高，染色體a04最低。品種更替過程中，20世紀80年代品種更替對SSR位點遺傳多樣性影響最為顯著，與20世紀70年代及以前品種相比，表現(xiàn)為等位基因遺傳豐富度增加、多樣性指數(shù)增加、品種間遺傳距離略有降低。20世紀80年代以后，SSR位點的等位基因豐富度增加、多樣性指數(shù)和品種間遺傳距離均無明顯變化。檢測到5個等位變異數(shù)隨年代增加減少的SSR位點。聚類分析結果顯示類群分布與系譜及地理來源相關，與品種更替年代無關。【結論】本研究結果表明，在花生品種更替過程中，主栽品種等位基因豐富度增加，而等位基因分布均勻度尚未產(chǎn)生顯著性改變。

花生；品種更替；主栽品種；SSR；遺傳多樣性

花生（Arachishypogaea L.）脂肪含量38%～60%，蛋白質(zhì)含量24%～36%，糖分含量20%左右，含有豐富的維生素 B2、PP、A、D、E，鈣和鐵等，是深受中國人喜愛的油料作物。目前，中國花生年均總產(chǎn)1342萬t，占世界花生總產(chǎn)的39.47%，居世界首位。種植面積為505.6萬公頃，占世界花生面積的19.12%，僅次于印度居世界第二位[1]?；ㄉN植效益高，是促進農(nóng)民增收的重要作物，在農(nóng)業(yè)供給側結構改革中占有重要地位。

自20世紀50年代以來，我國花生經(jīng)歷了5次品種更替，每次更替均顯著提高了單位面積的產(chǎn)量，使花生生產(chǎn)得到飛速發(fā)展。20世紀50年代，伏花生和獅頭企等品種的鑒定和大面積推廣，實現(xiàn)了我國花生品種的第一次更替；20世紀70年代末，隨著花17、粵油551、天府3號和開農(nóng)8號等品種的育成推廣，實現(xiàn)了第二次品種更替；20世紀80年代，隨著海華1號、白沙1016、天府7號和冀油4號等品種的育成推廣，實現(xiàn)了第三次品種更替；20世紀90年代，隨著粵油223、8130、冀花2號和中花4號等品種的育成推廣，實現(xiàn)了第四次品種更替；進入21世紀以來，高產(chǎn)、優(yōu)質(zhì)和抗病等專用型花生新品種逐步替代非專用型花生品種，實現(xiàn)了我國花生品種的第五次更替[2]。然而，近年來花生品種改良遺傳進度放緩，2011年花生單產(chǎn)為3502.50 kg/hm2，2015年為3561.75 kg/hm2，每年增幅僅為0.34%（http://www.zzys.moa.gov.cn/）。遺傳基礎狹窄是造成遺傳改良進度放緩的可能原因，育種過程中品種的更替通常會引起遺傳多樣性降低[3-4]。

微衛(wèi)星標記 SSR（simple sequence repeats）是從分子水平上探討遺傳多樣性的實用工具。已被廣泛應用于花生的遺傳多樣性分析、遺傳圖譜的構建、QTL定位和雜交種鑒定等方面[1，5-7]。由于SSR標記具有多態(tài)性豐度高的特點，致使傳統(tǒng)銀染法檢測容易產(chǎn)生實驗誤差，同一試驗中不易準確區(qū)分差異小的等位變異；不同試驗之間更難判斷等位變異是否一致。毛細管電泳檢測實現(xiàn)了數(shù)據(jù)收集和處理的自動化，克服了銀染法的不足，有效的降低了人為試驗誤差[8]，提高了試驗精度，并使不同試驗結果有效整合成為可能。

近年來，國內(nèi)外學者利用SSR標記對野生花生[9-10]、栽培花生的不同類型間[10-13]、同一栽培花生類型的不同品種間[14-15]和花生核心種質(zhì)[16-19]進行了遺傳多樣性分析。然而，由于年代久遠，花生品種更替過程中涉及的一些品種分布分散、不易收集，針對歷史上5次品種更替前后，主栽品種遺傳多樣性變化的專門研究尚未見報道。

本研究對中華人民共和國成立以來育成并大面積推廣的不同時期花生主栽品種進行系統(tǒng)搜集整理，并應用SSR標記毛細管電泳檢測，對我國花生主栽品種遺傳多樣性變化進行分析，揭示品種更替對花生主栽品種SSR位點遺傳多樣性演變的影響，為花生遺傳育種和品種推廣等提供依據(jù)和參考。

1 材料和方法

1.1 試驗材料

本研究選取我國主栽花生品種共計68份為材料（見表1），品種主要來源于本單位種質(zhì)資源圃和國家油料作物中期庫。其中包含北方花生區(qū)（河北、河南、江蘇、山東、安徽）38個花生品種，長江流域花生區(qū)（四川、湖北、湖南）12個花生品種和南方花生區(qū)（福建、廣東、廣西）18個花生品種。通過參考《中國花生品種及其系譜》[2]和《全國農(nóng)作物主要品種推廣情況統(tǒng)計表》及在地方農(nóng)科院相關人員的幫助下，將材料分為4個時期進行分析（由于材料數(shù)目的限制，本研究中將70年代及以前的材料歸為一類進行分析），其中包含20世紀70年代及以前花生品種7個，80年代花生品種15個，90年代和21世紀花生品種均為23個。

1.2 試驗方法

剪取來源于15～20個單株的0.1 g新鮮花生葉片用液氮研成粉末，置于2mL離心管中，然后采用生工生物工程（上海）有限公司的基因組DNA快速抽提試劑盒（植物）進行DNA樣品制備。DNA提取完成后用分光光度計測量其濃度和純度，并將部分DNA稀釋成20 ng/μL，置于4℃?zhèn)溆谩?/p>

本研究中共選取多態(tài)性高且覆蓋花生20對染色體的154對SSR引物對材料進行掃描，引物的序列及所在染色體參考文獻[5，20-22]。采用普通引物加M13通用引物的方法檢測，即引物合成時前導鏈5′端添加M13序列[23]。采用FAM、VIC、NED和PET共4種不同顏色的熒光染料直接標記M13互補鏈的每個堿基。PCR 反應體系（15μL） 為：7.5μL 2×power Taq PCR Master Mix（美國 Microanalysis Inc.），0.2μL 2μmol/L Primers 1（攜帶 M13），1μL 2μmol/L Primers 2，0.8μL 2μmol/L M13，2μL DNA，3.5μL ddH2O。反應程序為：94℃預變性5min；94℃變性45 s；50℃退火 45 s；72℃延伸 45 s，共 34個循環(huán)；94℃變性 45 s；53℃退火 45 s；72℃延伸 45 s，共 8個循環(huán)；72℃延伸10min。PCR擴增產(chǎn)物通過ABI3730 DNA分析儀進行毛細管電泳檢測，檢測結果利用Genemapper 3.2軟件（Applied Biosystems）進行數(shù)據(jù)采集。

表1 主栽品種及其育成年代Table 1 The origin and the years released of peanut elite cultivars

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

利用PowerMarkerV3.25軟件[24]統(tǒng)計每個SSR位點的等位變異豐富度（A）、多態(tài)性信息指數(shù)（PIC）、品種間的遺傳距離（D）等。同時利用此軟件進行群體間的遺傳距離的估算并以除權配對法（UPGMA）進行聚類分析。用FigtreeV1.4.3（http://tree.bio.ed.ac.uk/software/figtree/）軟件對輸出圖文件進行編輯并繪制聚類圖。

多態(tài)性信息指數(shù)：

其中，pij表示第i位點第j個等位變異在全部供試材料中出現(xiàn)的頻率。

全基因組的遺傳豐富度：

其中，Aij表示基因組第i個基因組第j個位點的等位變異數(shù)目；n表示檢測位點總數(shù)。

全基因組的多樣性指數(shù)：

其中，PICi表示第i個位點的多態(tài)性信息指數(shù)；n表示檢測到的位點總數(shù)。

年代內(nèi)品種間遺傳距離：

其中，Dij表示基因組第i個位點與第j個位點的遺傳距離；n表示檢測位點總數(shù)。

2 結果與分析

2.1 品種更替前后主栽品種遺傳多樣性變化

2.1.1 參試品種多樣性參數(shù)分析

選用154對SSR引物對主栽品種進行檢測，共獲得173個位點，其中Ah1TC3B05等19對引物擴增得到2個位點。173個SSR位點共檢測到872個等位變異，位點的平均等位變異數(shù)為5.04個，等位變異最小為2個，分別為seg12E03、GA79和TC3G05等32個位點。等位變異數(shù)目最多為15個，為標記seq16F10。等位基因中有200個稀有等位基因（出現(xiàn)頻率小于0.05），占總檢出等位基因數(shù)的22.94%。共檢測到1108個雜合位點，占所有多態(tài)性位點（11150個）的9.94%。主栽品種遺傳多態(tài)性指數(shù)的變化范圍為0.014～0.881，平均0.477，多樣性最高位點為TC1A08，最低位點為Ah-408。

對20條染色體的等位基因遺傳豐富度和多樣性指數(shù)比較，不同染色體上的SSR位點等位基因遺傳豐富度差別較大，如圖1所示，染色體b07最高為7.20個，其次為a10、a01和b04染色體，分別達到了7.00個、6.42個和6.00個，其余染色體都低于6個，a04、a09和b08最低均為4.00個。多樣性指數(shù)分析結果表明，染色體b07最高為0.65，其次為a07為0.56，染色體b06最低0.28，染色體a04較低為0.37。結合二者進行分析，20條染色體中，b07遺傳多樣性最高，染色體a04較低。

圖1 20條染色體的等位基因豐富度和多樣性指數(shù)比較Figure 1 The comparisons of genetic richness of allele and gene diversity index in twenty chromosomes

2.1.2 4次品種更替引起的多樣性變化

為探討我國花生主栽品種遺傳多樣性的演變趨勢，在全基因組遺傳豐富度、多樣性指數(shù)和年代內(nèi)品種間遺傳距離等方面對4個群體進行評價分析，結果如表2所示，我國花生全基因組遺傳豐富度是逐漸增大的趨勢，由20世紀70年代及之前的2.80增加到21世紀初期的4.06。同時其多樣性指數(shù)也由20世紀70年代及以前的0.42增大到20世紀80年代的0.45，但20世紀90年代和21世紀初期均為0.46與80年代相比，差異不大。以上結果表明，近年來隨著育種的發(fā)展，我國花生基因組遺傳豐富度有所增加，品種遺傳基礎得到拓寬。但品種的遺傳多樣性自20世紀80年代以來，并未得到較大提高。

對4個年代主栽品種的年代內(nèi)品種間遺傳距離進行統(tǒng)計計算（見表2），不同年代間遺傳距離變異范圍為0.5223～0.5407。其中20世紀70年代及以前的7個主栽品種遺傳距離最大為0.5407，到20世紀80年代減小為0.5284，但20世紀80年代以后，主栽品種的遺傳距離接近。

表2 不同階段花生主栽品種遺傳多樣性分析Table 2 The genetic diversity of peanut elite cultivars at different period

2.1.3 不同SSR位點對品種更替響應程度分析

對173個SSR位點在不同年代花生主栽品種間等位變異的變化進行分析，結果表明，29個位點在不同年代花生主栽品種間的等位變異未發(fā)生改變；69個位點隨著年代的增加等位變異逐漸增大；70個位點隨著年代的增加等位變異出現(xiàn)先增大后減少的趨勢；5個位點的等位變異隨著年代的增加而減少。表3中為隨著年代的增加等位變異減少的5個標記。

表3 SSR位點在不同年代花生品種間等位變異數(shù)目Table 3 The number of the alleles of different periods peanut elite cultivars

對173個位點的等位變異數(shù)與多態(tài)性信息指數(shù)進行相關分析，結果表明二者之間存在極顯著正相關（P＜0.01），但二者之間線性關系中等（相關系數(shù)R＜0.8）。在等位變異相同的情況下，多態(tài)性信息指數(shù)存在較大差異，例如PM31和TC3E05檢測到等位變異數(shù)目都為8個，多態(tài)性信息指數(shù)分別為0.65和0.80。圖2即為2位點等位變異及其頻率分布，從圖中可以看出PM31等位變異分布不均勻，其中200 bp的頻率遠高于其他等位變異，推測此位點可能與某些重要農(nóng)藝性狀基因緊密連鎖，由于受到較強的選擇壓保留下來[3，25]。TC3E05檢測結果與PM31相反，等位變異分布相對均勻。

圖2 PM31和TC3E05位點等位變異及其頻率Figure 2 Alleles and allele frequency of the locus PM31 and TC3E05

2.2 不同年代育成品種聚類分析

68份花生材料基于Nei’s遺傳距離的系統(tǒng)聚類分析結果示于圖3，所選材料分為4個類群。第Ⅰ類群有29個材料，主要為來源于河南（7個，占24.14%）、山東（7個，占 24.14%）、河北（6個，占20.69%）。來源于20世紀70年代及以前、80年代、90年代和21世紀初期的品種分別為2、7、8和12個；第Ⅱ類群包含22個品種，主要來源于廣東（8個，占36.36%）和廣西（5個，占22.73%）。來源于20世紀70年代及以前、80年代、90年代和21世紀初期的品種分別為4、5、10和3個；第Ⅲ類群包含11個品種，其中河南和山東材料最多均為4個，其次為四川2個。來源于20世紀70年代及以前、80年代、90年代和21世紀初期的品種分別為1、2、3和5個；第Ⅳ類群有6個品種，其中包含山東和河南各2個。來源于80年代、90年代和21世紀初期的品種分別為1、2和3個。從每個小類群的組成可見，相同省份的材料并未完全聚類于同一小類群，但具有一定的傾向性，如河北省的品種都在第Ⅰ類群中，廣東和廣西品種大都也在第Ⅱ類群中。第Ⅱ類群主要為獅頭企、獅頭企與伏花生雜交育成品種及其衍生后代。但每個類群中均包含不同年代的品種。從結果可以看出，類群分布與親緣關系及地域相關，與品種更替時間無關。

對來源于不同年代花生主栽品種兩兩間進行遺傳距離的計算，表4可以看出20世紀70年代及以前與其他時期間的遺傳距離較大，為0.2076～0.2386，并隨著年代的增加逐漸增大。20世紀80年代主栽品種與20世紀90年代和21世紀主栽品種之間遺傳距離較小在0.1641～0.1651之間。20世紀90年代和21世紀主栽品種遺傳距離為0.1754。

圖3 68個花生品種的SSR標記遺傳聚類圖Figure 3 Dendrograph of 68 peanut varieties based on SSR genetic distance

表4 不同階段花生主栽品種間遺傳距離分析Table 4 The genetic distance of peanut elite cultivars between different periods

3 討論與結論

3.1 不同研究中花生SSR多態(tài)性的比較

近年來，SSR標記已廣泛用于花生種質(zhì)資源的遺傳多樣性檢測中。姜慧芳等[16]對我國298份小核心種質(zhì)進行SSR多樣性分析，花生資源平均多態(tài)性條帶數(shù)為4.8，遺傳多樣性指數(shù)為0.22。Wang M.L.等[19]對美國小核心種質(zhì)112份進行SSR多樣性分析，每個位點檢測到的等位基因數(shù)平均為8.1個，PIC為0.53。近年來，國內(nèi)外學者利用SSR標記對栽培種進行多態(tài)性檢測也屢見報道，Wang H.等[26]對中國、印度和美國的79份栽培種和高代品系進行了分析，等位基因豐富度為4.25，PIC為0.43。任小平分析了我國196份花生栽培種，其等位基因豐富度為 3.0，PIC 為 0.38[1]。K.M.Pandey等[27]對 20 個花生親本進行遺傳多樣性分析，其等位基因豐富度為3.2，PIC為0.31。本研究對我國68份花生主栽品種進行遺傳多樣性分析，等位基因豐富度為5.04，PIC為0.48。不同的研究得到的結果不同，主要是由于以下幾個方面的原因造成的：一是供試材料的選擇，本研究中不同年代主栽品種的等位基因豐富度和遺傳多樣性指數(shù)存在很大差異，花生中也已經(jīng)被證明地方種的多樣性比育成品種豐富[12]。因此，可通過在雜交育種過程中篩選鑒定更多優(yōu)異地方種作親本，進而拓寬我國花生品種的遺傳基礎；二是SSR標記的選擇，選擇擴增產(chǎn)物穩(wěn)定、多態(tài)性高的標記可大大提高檢測效率，節(jié)約檢測時間和檢測成本，因此，可優(yōu)先選擇他人文章中發(fā)表的多態(tài)性標記進行擴增。本研究中選用的標記主要來源于文獻[5，20-22]，都是遺傳圖譜中兩親本間分離的標記；三是檢測方法的選擇，常規(guī)的SSR銀染技術在大規(guī)模、多批次的數(shù)據(jù)收集和分析時仍存在相當大的難度，如不同等位變異的準確識別困難及不同批次反應數(shù)據(jù)難以統(tǒng)一。SSR熒光標記采用了熒光標記技術，利用ABI3700 DNA分析儀進行擴增產(chǎn)物的檢測，采用Genemapper3.2軟件進行檢測數(shù)據(jù)的采集，實現(xiàn)了數(shù)據(jù)收集和處理的自動化，明顯降低實驗結果的判讀錯誤，提高數(shù)據(jù)分析的可靠性和準確度，并且經(jīng)分析熒光標記技術較銀染方法在每個位點上多檢測到3個等位變異，具有更高的靈敏性，更適于進行遺傳多樣性分析和研究[8]。因此，本研究中利用熒光標記技術對材料進行檢測得到較好的多態(tài)性結果。

3.2 品種更替對花生主栽品種遺傳多樣性的影響

我國花生品種改良起步較晚，中華人民共和國成立后生產(chǎn)上利用的大都為農(nóng)家品種，20世紀50年代鑒定的農(nóng)家品種伏花生和獅頭企大面積推廣利用，實現(xiàn)了我國花生品種的第一次更替。后來隨著育種工作的蓬勃開展分別于20世紀70年代、80年代、90年代和21世紀又經(jīng)歷了4次花生品種更替[2]。本研究中，我國花生主栽品種的等位基因豐富度隨著年代的增加逐漸增加，這主要是由于我國花生主栽品種雖然大都有伏花生和獅頭企的血緣，但是隨著年代的增加，在雜交過程中開始利用篩選到的優(yōu)良地方種及國外引進優(yōu)良品種作親本，育種方法也開始由地方種的鑒定篩選到單交、復交、誘變等多種途徑，因此，檢測到的等位基因數(shù)目逐漸遞增。但是在20世紀80年代以后主栽品種的多樣性指數(shù)未出現(xiàn)明顯差異，遺傳距離也是隨著年代的增加而降低，雖在2000年以后略有增加，但較20世紀70年代及以前的品種還是降低明顯。這是由于我國主栽品種雖然在雜交過程中利用了地方品種和國外引進品種，但大都還是利用育成種作親本進行修修補補，育成的品種具有的某些特有等位基因多是稀有基因，對群體遺傳結構的影響不大，因而通過SSR分析結果計算年代內(nèi)品種間遺傳距離時得到數(shù)值并不高。因此，基于以上結果可確認我國花生主栽品種遺傳多樣性自20世紀80年代以來雖全基因組遺傳豐富度得到了提高，但等位基因分布均勻度尚未產(chǎn)生顯著性改變。

3.3 SSR多態(tài)性與選擇牽連效應的關系

張學勇等[25]指出，在基因組中，一些承受強選擇作用的基因在群體中的多樣性顯著降低，同時這些基因附近區(qū)域的遺傳多樣性也明顯下降。在遺傳學中將這種對個別基因的選擇導致其側翼區(qū)域遺傳多樣性降低的現(xiàn)象稱之為選擇牽連效應（hitch hiking effect，也稱選擇搭載效應）。自然選擇和人工選擇使被選擇位點兩側的多樣性降低，在育種過程中育種家感興趣的目標性狀如產(chǎn)量性狀、品質(zhì)性狀、抗病性等相關基因的選擇，可能存在較強的選擇牽連作用。對20條染色體進行比較，a04染色體的遺傳多樣性最低，推測a04染色體在育種過程中承受了較大的選擇壓，在其中也檢測到許多于脂肪酸組成相關的QTL位點[28]。5個等位變異減少的標記中，PM13已被證明與花生油酸含量[29]和黑斑病抗性[30]相關；TC3H07[31]與花生株型有關；PM434和 PM377[32]與花生銹病抗性相關。

3.4 聚類結果與品種更替的關系

利用154對SSR引物擴增結果，將68份材料劃分為4個類群，材料類群劃分在與地理來源相關的同時，也可以發(fā)現(xiàn)與親緣關系有關，例如：第Ⅱ類群主要為獅頭企、獅頭企與伏花生雜交育成品種及其衍生后代。除第四類群外，每一類群中也都包含不同年代的主栽品種。因此，聚類結果中類群的劃分與地理來源和親緣關系有關，與育成年代無關。我們要拓寬育成品種的遺傳多樣性，可以通過挑選不同地理來源及親緣關系遠的材料做親本，不用考慮其育成年代。

以上結果表明，我國花生自20世紀80年代以來的主栽品種間全基因組遺傳豐富度逐漸增大，但遺傳多樣性指數(shù)和年代內(nèi)品種間的遺傳距離變化不明顯，可能是影響花生適應性和產(chǎn)量潛力增加的原因。因此，今后花生育種要在進一步引進國外優(yōu)異資源、鑒定篩選優(yōu)異的地方資源和挖掘野生近緣植物中的優(yōu)良基因的同時，結合現(xiàn)在基因克隆、分子標記和轉基因等工作，不斷拓寬和創(chuàng)新花生育成品種的遺傳基礎，積極推動我國花生育種及花生產(chǎn)業(yè)的蓬勃發(fā)展。

［1］REN X P，JIANG H F，YAN Z Y，et al.Genetic diversity and population structure of the major peanut(Arachis hypogaea L.)cultivars grown in China by SSR markers［J］.Plos One，2014，9（2）：88091.

［2］禹山林.中國花生品種及其系譜［M］.上海：上?？茖W技術出版社，2008：10-11.

［3］郝晨陽，王蘭芬，張學勇，等.我國育成小麥品種的遺傳多樣性演變［J］.中國科學，2005，35（5）：408-415.

［4］劉傳光，張桂權.用SSR標記分析1949-2005年華南地區(qū)常規(guī)秈稻主栽品種遺傳多樣性及變化趨勢［J］.作物學報，2010，36（11）：1843-1852.

［5］GAUTAMI B，F(xiàn)ONCEKA D，PANDEY M K，et al.An international reference consensus genetic map with 897 marker loci based on 11 mapping populations for tetraploid groundnut(Arachis hypogaea L.)［J］.Plos One，2012，7（7）：41213.

［6］PANDEY M K，WANG M L，QIAO L X，et al.Identification of QTLs associated with oil content and mapping FAD2 genes and their relative contribution to oil quality in peanut(Arachis hypogaea L.)［J］.BMC Genetics，2014，15（1）：1-14.

［7］王瑾，李玉榮，張嘉楠，等.花生雜交種F1真?zhèn)蔚腟SR標記檢測［C］.武漢：中國油料作物學報，2013（35）：171-175.

［8］郝晨陽，王蘭芬，賈繼增，等.SSR熒光標記和銀染技術的比較分析［J］.作物學報，2005，31（2）：144-149.

［9］MORETZSOHN M C，HOPKINS M S，MITCHELL S E，et al.Genetic diversity of peanut(Arachis hypogaea L.)and its wild relatives based on the analysis of hypervariable regions of the genome［J］.BMC Plant Biology，2004，4（1）：11.

［10］BARKLEY N A，DEAN R E，PITTMAN R N，et al.Genetic diversity of cultivated and wild-type peanuts evaluated with M13-tailed SSR markers and sequencing［J］.Genetics Research，2007，4（89）：93-106.

［11］崔順立，劉立峰，陳煥英，等.河北省花生地方品種基于SSR標記的遺傳多樣性［J］.中國農(nóng)業(yè)科學，2009，42（9）：3346-3353.

［12］姜慧芳，任小平，陳玉寧，等.中國花生地方品種與育成品種的遺傳多樣性［J］.西北植物學報，2011，31（8）：1551-1559.

［13］HE G，MENG R，NEWMAN M，et al.Microsatellites as DNA markers in cultivated peanut(Arachis hypogaea L.)［J］.BMC Plant Biology，2003，3（1）：3.

［14］唐榮華，莊偉建，高國慶，等.珍珠豆型花生的簡單序列重復(SSR)多態(tài)性［J］.中國油料作物學報，2004，26（2）：20-26.

［15］張建成，王傳堂，楊新道.SSR和STS標記在花生栽培品種鑒定中的應用研究［J］.植物遺傳資源學報，2006，7（2）：215-219.

［16］姜慧芳，任小平，張曉杰，等.中國花生小核心種質(zhì)SSR遺傳多樣性［J］.中國油料作物學報，2010，32（4）：472-478.

［17］JIANG H，HUANG L，REN X，et al.Diversity characterization and association analysis of agronomic traits in a Chinese peanut(Arachis hypogaea L.)mini-core collection［J］.Journal of Integrative Plant Biology，2014，56（2）：159-169.

［18］任小平，張曉杰，廖伯壽，等.ICRISAT花生微核心種質(zhì)資源SSR 標記遺傳多樣性分析［J］.中國農(nóng)業(yè)科學，2010，43（14）：2848-2858.

［19］WANG M L，SUKUMARAN S，BARKLEY N A，et al.Population structure and marker trait association analysis of the US peanut(Arachis hypogaea L.)mini-core collection［J］.Theoretical Applied Genetics，2011，123（8）：1307-1317.

［20］SHIRASAWA K，BERTIOLI D J，VARSHNEY R K，et al.Integrated consensus map of cultivated peanut and wild relatives reveals structures of the A and B genomes of Arachis and divergence of the legume genomes［J］.DNA Research，2013，20（2）：173.

［21］HONG Y，CHEN X，LIANG X，et al.A SSR-based composite genetic linkage map for the cultivated peanut(Arachis hypogaea L.)genome［J］.BMC Plant Biology，2010（10）：17.

［22］HUANG L，HE H，CHEN W，et al.Quantitative trait locus analysis of agronomic and quality-related traits in cultivated peanut(Arachis hypogaea L.)［J］.Theoretical Applied Genetics，2015（128）：1103-1115.

［23］OETTINH W，LEE H，F(xiàn)LANDERS D，et al.Linkage analysis with multiplexed short tandem repeat polymorphisms using infrared fluorescence and M13 tailed primers［J］.Genomics，1996，30（3）：450-458.

［24］LIU K，MUSE S V.PowerMarker：an integrated analysis environment for genetic marker analysis［J］.Bioinformatics，2005，21（9）：2128-2129.

［25］張學勇，童依平，游光霞，等.選擇牽連效應分析：發(fā)掘重要基因的新思路［J］.中國農(nóng)業(yè)科學，2006，39（8）：1526-1535.

［26］WANG H，KHERA P，HUANG B，et al.Analysis of genetic diversity and population structure of peanut cultivars and breeding lines from China，India and the US using SSR markers［J］.Journal of Integrative Plant Biology，2016，58（5）：452-465.

［27］PANDEY M K，GAUTAMI B，JAYAKUMAR T，et al.Highly informative genic and genomic SSR markers to facilitate molecular breeding in cultivated groundnut(Arachis hypogaea L.)［J］.Plant Breeding，2012，131（1）：139-147.

［28］WANG M L，KHERA P，PANDEY M K，et al.Genetic mapping of QTLs controlling fatty acids provided insights into the genetic control of fatty acid synthesis pathway in Peanut(Arachis hypogaea L.)［J］.Plos One，2015，10（4）：0119454.

［29］SARVAMANGALA C，GOWDA M V C，VARSHNEY R K，et al.Identification of quantitative trait loci for protein content，oil content and oil quality for groundnut(Arachis hypogaea L.)［J］.Field Crops Research，2011，122（1）：49-59.

［30］GAUTAMI B，PANDEY M K，VADEZ V，et al.Quantitative trait locus analysis and construction of consensus genetic map for drought tolerance traits based on three recombinant inbred line populations in cultivated groundnut(Arachis hypogaea L.)［J］.Molecular Breeding，2012：DOI 10.1007/s11032-011-9660-0.

［31］FONCEKA D，TOSSIM H A，RIVALLAN R，et al.Fostered and left behind alleles in peanut：interspecific QTL mapping reveals footprints of domestication and useful natural variation for breeding［J］.BMC Plant Biology，2012，12（1）：26.

［32］KHEDIKAR Y P，GOWDA M V C，SARVAMANGALA C.A QTL study on late leaf spot and rust revealed one major QTL for molecular breeding for rust resistance in groundnut(Arachis hypogaea L.)［J］.Theoretical Applied Genetics，2010，121（5）：971-984.

Analysis the Genetic Diversity Variation Trends of Peanut Elite Cultivars in China using SSR Markers

WANG Jin，LI YU-rong*，CHENG Zeng-shu，CHEN Si-long，SONG Ya-hui，ZHANG Peng-juan
（Institute of Cereal and Oil Crops，Hebei Academy of Agriculture and Forestry Sciences/Laboratory of Crop Genetics and Breeding of Hebei Province，Shijiazhuang 050035，China)

【Objective】The objective of this study was to detect the genetic diversity in China peanut cultivar replacement since 1950s based on SSR marker in order to provide references for peanut breeding.【Method】The genotype of 68 elite varieties was screened using 154 SSR markers.【Results】As results，173 loci and 872 alleles were detected.The mean number of allele per locus was 5.04.The polymorphism information content(PIC)value ranged from 0.014 to 0.881，with average 0.477.Among 20 chromosomes，b07 had highest PIC value while a04 had lowest one.The peanut cultivar evolution in 1980s had the most significant impact on genetic diversity.The varieties that released in this stage had more allele，higher PIC value and similar genetic distance compared with the varieties that released before 1980s.The varieties that released after 1980s has more allele number than others，but which had no significant impact on PIC value and genetic distance.Furthermore，the allele number of 5 SSR loci showed significant decrease with cultivar evolution.The varieties with related pedigree or geographic origin were assigned into same group，while cluster result had no relationship with released year of varieties.【Conclusion】Our results showed that the allele number had increased while the evenness did not change in the past 50 years.

Arachishypogaea；cultivar evolutions；elite cultivar；SSR marker；genetic diversity

S565.202

A

1000-2650（2017）03-0309-08

10.16036/j.issn.1000-2650.2017.03.004

2017-05-18

國家自然科學基金項目（31301354）；河北省科技計劃項目（16226301D）；國家花生產(chǎn)業(yè)技術體系建設項目（CARS-14）。

王瑾，副研究員，博士，主要從事花生遺傳育種與栽培研究，E-mail：wangjinnky@163.com。*責任作者：李玉榮，研究員，學士，主要從事花生遺傳育種與栽培研究，E-mail：liyrl@163.com。

（本文審稿：任小平；責任編輯：劉詩航；英文編輯：劉詩航）