魏中艷 李慧慧 李 駿 Yasir A. Gamar 馬巖松 邱麗娟,*

?

應(yīng)用SNP精準(zhǔn)鑒定大豆種質(zhì)及構(gòu)建可掃描身份證

魏中艷1李慧慧1李 駿2Yasir A. Gamar1馬巖松3邱麗娟1,*

1國家農(nóng)作物基因資源與遺傳改良重大科學(xué)工程/ 農(nóng)業(yè)部種質(zhì)資源利用重點實驗室/ 中國農(nóng)業(yè)科學(xué)院作物科學(xué)研究所, 北京100081;2廣西大學(xué)農(nóng)學(xué)院, 廣西南寧 530004;3黑龍江省農(nóng)業(yè)科學(xué)院大豆研究所, 黑龍江哈爾濱 150086

為加強(qiáng)大豆種質(zhì)資源管理及品種保護(hù), 本研究構(gòu)建了一套基于單核苷酸多態(tài)性(SNP)標(biāo)記的快速鑒定體系。選用分布于13個基因的23個SNP標(biāo)記對599份大豆表型精準(zhǔn)鑒定種質(zhì)進(jìn)行基因型分析。結(jié)果表明, SNP標(biāo)記的遺傳多樣性指數(shù)范圍0~0.722, 其中3個SNP在供試材料中不存在堿基差異, 2個SNP為特異等位變異。從具有多態(tài)性的20個SNP中選出14個(GlySNP14)用于種質(zhì)鑒定, 其中12個為高多態(tài)性SNP, 2個為特異性SNP。模擬結(jié)果表明, GlySNP14的種質(zhì)鑒別能力(750份種質(zhì))高于任意相同數(shù)目標(biāo)記構(gòu)成的SNP隨機(jī)組合(最多361份種質(zhì))。GlySNP14在599份種質(zhì)中共形成了176個單倍型, 其中100份種質(zhì)具有獨特的單倍型。結(jié)合這些種質(zhì)的其他屬性, 構(gòu)建了100份種質(zhì)由38個數(shù)字組成的身份證, 前10位數(shù)字為種質(zhì)屬性信息碼, 包括品種類別、來源地等; 11~38位數(shù)字為品種分子指紋碼, 代表品種的特異分子信息, 并最終以一維碼和二維碼的形式表示, 為種質(zhì)資源簡易管理與保護(hù)利用提供了有效途徑。

大豆; SNP; 遺傳多樣性; DNA指紋圖譜; 品種ID

大豆(L.)是豆科(Fabaceae)大豆屬(Glycine)一年生草本作物。大豆起源于中國已得到國際學(xué)者普遍認(rèn)可[1-2], 我國大豆種質(zhì)資源豐富, 保存的數(shù)量居全世界之首[3]。大豆也是我國重要的經(jīng)濟(jì)作物之一, 在全國都有廣泛種植[4], 到目前為止, 已育成和推廣大豆品種超過1800個[5]。隨著大豆種質(zhì)資源和品種的交換, 種子市場中“同物異名”和“套牌”、“冒牌”現(xiàn)象日趨嚴(yán)重, 這不僅給市場管理和消費者帶來困擾, 更會給品種的選育、經(jīng)營及種植等相關(guān)人員和單位帶來經(jīng)濟(jì)損失。因此, 種質(zhì)資源和品種身份證的構(gòu)建具有重要意義。

作物品種資源種子身份證的構(gòu)建已從形態(tài)標(biāo)記向高通量分子鑒定技術(shù)發(fā)展。傳統(tǒng)的鑒定方法主要是形態(tài)標(biāo)記、生理標(biāo)記和生化指標(biāo), 隨著品種資源數(shù)目的增加, 僅僅基于這些特性已經(jīng)難以對現(xiàn)有品種準(zhǔn)確鑒定[6]。因此急需開發(fā)一套準(zhǔn)確可靠、簡單易行的快速鑒定體系, 為作物新品種審定及種質(zhì)材料評價提供科學(xué)依據(jù)。20世紀(jì)以來, DNA分析技術(shù)開始被大量應(yīng)用于植物學(xué)研究, 開發(fā)出RFLP、RAPD、AFLP、SSR以及SNP、SRAP等一系列標(biāo)記[7]。分子標(biāo)記不但能夠節(jié)省常規(guī)田間調(diào)查和收集整理數(shù)據(jù)的時間, 而且具有不受環(huán)境影響、鑒別品種準(zhǔn)確和變異極豐富等優(yōu)點[8], 利用分子標(biāo)記構(gòu)建DNA指紋圖譜進(jìn)行品種真?zhèn)涡砸殉蔀榘l(fā)展趨勢。鑒于方法的穩(wěn)定和有效性, 國際植物品種權(quán)保護(hù)聯(lián)盟(UPOV)在BMT測試指南草案中已將構(gòu)建DNA指紋數(shù)據(jù)庫的標(biāo)記方法確定為SSR和SNP[9]。較SSR標(biāo)記相比, SNP具有針對性強(qiáng)、變異來源豐富、潛在數(shù)量巨大等特點。Jung等[10]首次在辣椒中分析438對COSII引物, 選出40個可以鑒別79個熱帶商業(yè)品種和17個甜椒品種的SNP標(biāo)記。Shirasawa 等[11]研究發(fā)現(xiàn), 利用8個SNP標(biāo)記可以區(qū)分43個水稻品種。王大莉[12]用72個SNP標(biāo)記構(gòu)建了22個香菇主栽品種的DNA指紋圖譜。利用SNP分子標(biāo)記鑒定作物的指紋圖譜, 對作物品種資源管理與品種保護(hù)利用及評價無疑將起到積極的作用。

迄今為止, 受制于SNP標(biāo)記自動化篩選及檢測成本高等問題, SNP標(biāo)記在大豆種質(zhì)資源的鑒定方面鮮有報道, 利用與表型性狀相關(guān)的功能性SNP鑒定資源, 且最終以條形碼的形式表現(xiàn)的研究更未見報道。本研究利用與主要農(nóng)藝性狀密切相關(guān)的功能基因SNP標(biāo)記, 精準(zhǔn)鑒定不同生態(tài)區(qū)的大豆種質(zhì), 分析其DNA指紋圖譜, 同時結(jié)合大豆品種屬性, 構(gòu)建其分子身份證, 以期為大豆遺傳研究、分子輔助育種和品種鑒定提供依據(jù)。

1 材料與方法

1.1 供試材料

北方種質(zhì)208份, 于2008年至2010年連續(xù)3年種植于北方地區(qū)3個不同的地點(黑龍江、吉林、內(nèi)蒙古); 黃淮種質(zhì)245份, 于2008年至2010年連續(xù)3年種植于黃淮地區(qū)3個不同的地點(山東、河南、河北); 南方種質(zhì)146份, 于2008年至2010年連續(xù)3年種植于南方地區(qū)3個不同的地點(湖北、江西、廣西)(圖1)。并對其茸毛色、結(jié)莢習(xí)性、百粒重、熟期、胞囊線蟲抗性等表型進(jìn)行多年多點鑒定。

1.2 DNA提取和純化

萌發(fā)后取每份大豆種質(zhì)5個單株的幼葉混樣, 使用Fermentas試劑盒提取基因組DNA。其后加2 μL稀釋20倍的RNase A, 放到37°C培養(yǎng)箱中15 min, 去除DNA中的RNA, 于–80°C保存待送樣。

圖1 大豆精準(zhǔn)鑒定種質(zhì)的種植分布圖

1.3 SNP標(biāo)記鑒定

選用的23個SNP位點中, 3個來自大豆結(jié)莢習(xí)性相關(guān)基因; 2個來自大豆茸毛色相關(guān)基因; 7個來自大豆熟期相關(guān)基因,,和; 4個來自大豆SCN抗性相關(guān)基因和; 6個來自大豆SMV抗性相關(guān)基因、、和; 1個來自大豆百粒重相關(guān)基因(表1)。采用Infmium芯片技術(shù)(Bioyong Technologies Inc.), 鑒定全基因組或特定SNP位點。

1.4 SNP多樣性計算

利用PowerMaker 3.25軟件[13]計算SNP標(biāo)記在精準(zhǔn)鑒定種質(zhì)中的分布頻率及遺傳相似系數(shù); 參照Poole[14]介紹的Shannon-Weaver公式, 多樣性指數(shù)= –∑PlnP, 式中P為每個位點的等位變異頻率。

1.5 SNP鑒別能力模擬

利用14個SNP在599份大豆種質(zhì)中的等位變異頻率, 通過R程序模擬分析(http://www.r-project. org/)。共設(shè)置8個模擬群體, 群體大小分別為50、100、500、1000、1500、2000、2200和2500。在20個SNP中(除去3個沒有堿基差異的標(biāo)記: TT14, PRO16, SER46)隨機(jī)選取14個SNP對上述8個群體進(jìn)行鑒別分析, 該過程一共模擬1000次。最后利用每個群體1000次模擬結(jié)果的平均值比較分析。

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

利用Flapjack軟件[15]分析單倍型; 利用MEGA5.1軟件對供試大豆品種進(jìn)行聚類分析, 并確定材料的遺傳結(jié)構(gòu); 利用條碼在線生成器(http:// barcode.tec-it.com/barcode-generator.aspx)對品種身份證進(jìn)行條碼轉(zhuǎn)換。

2 結(jié)果與分析

2.1 SNP標(biāo)記多樣性分析

本研究23個SNP標(biāo)記共來自于13個基因, 分別與結(jié)莢習(xí)性、茸毛色、熟期、百粒重、大豆胞囊線蟲病與花葉病毒病抗性相關(guān)(表1), 用這些標(biāo)記分析599份種質(zhì), 每個SNP標(biāo)記其遺傳多樣性指數(shù)變異范圍為0~0.722, 平均值為0.312。在23個SNP標(biāo)記中, 9個發(fā)生顛換, 變異范圍為0.50%~36.93%; 5個發(fā)生缺失, 變異范圍為0.67%~99.16%; 6個為轉(zhuǎn)換, 變異范圍為0.17%~43.79%。絕大多數(shù)位點都有2~3種基因型, 只有3個位點僅有一種基因型, 分別是TT14、PRO16和SER46。

利用599份種質(zhì)的23個SNP鑒定數(shù)據(jù)分析表明, 精準(zhǔn)鑒定種質(zhì)間的平均遺傳相似性系數(shù)為0.7923, 以“中豆24”與“茶色豆”遺傳差異最大, 為0.0435, 而“Magnolid”與“Nova”遺傳差異最小, 為0.9783。根據(jù)中國大豆3個生態(tài)區(qū)劃分(表2), 北方生態(tài)區(qū)(吉林、黑龍江和內(nèi)蒙古3個省、區(qū)) 208個品種間的平均遺傳相似系數(shù)為0.8295; 黃淮海生態(tài)區(qū)(包括山東、河南與河北3個省) 245個品種間的平均遺傳相似系數(shù)為0.7941; 南方生態(tài)區(qū)(包括湖北、江西和廣西3個省、區(qū)) 146個品種間的平均遺傳相似系數(shù)為0.7982。表明品種間遺傳差異以北方生態(tài)區(qū)最小, 南方生態(tài)區(qū)次之, 黃淮海生態(tài)區(qū)最大。

2.2 GlySNP14的構(gòu)建

通過對23個SNP位點的等位變異分析, 發(fā)現(xiàn)兩類特殊的SNP, 一類是2個具有特異等位變異的SNP位點(TT12和AT22), 僅用該位點即可區(qū)分“半野生”和“黃豆<2>”, 也可將這2份種質(zhì)與其他種質(zhì)區(qū)分開來。另一類是3個沒有堿基差異的SNP位點, 包括TT14、PRO16和SER46, 這3個位點不能區(qū)分任何一種參試的種質(zhì), 不適宜用于身份證構(gòu)建(表1)。具有多態(tài)性的20個SNP對599份種質(zhì)進(jìn)行單倍型分析, 共形成225個單倍型(圖2-A), 其中142個為特異單倍型(圖2B), 即20個SNP可以將這142份種質(zhì)完全區(qū)分。

在去除上述兩類特殊等位變異的基礎(chǔ)上, 進(jìn)一步對18個多態(tài)性SNP進(jìn)行組合分析, 采用隨機(jī)組合和按SNP多樣性指數(shù)由高到低逐一組合兩種方法, 選出鑒別品種能力最佳標(biāo)記組合GlySNP14。如圖3所示, 選用多樣性指數(shù)高的SNP組合的鑒別種質(zhì)效率明顯優(yōu)于隨機(jī)組合, 選用多樣性指數(shù)最高的2個SNP標(biāo)記PRO19和SER43進(jìn)行組合, 可以鑒別出1份種質(zhì), 平均每個標(biāo)記鑒別0.5份種質(zhì), 而隨機(jī)選用AR2和AR5進(jìn)行組合, 不能鑒別出任何種質(zhì); 在上述基礎(chǔ)上, 增加TT8與PRO19和SER43進(jìn)行組合, 可以鑒別出3份種質(zhì), 平均每個SNP鑒別效率為0.5,而利用隨機(jī)選取的AR7與AR2和AR5進(jìn)行組合, 僅能鑒別出1份種質(zhì), 平均每個SNP鑒別效率為0.33。當(dāng)多樣性指數(shù)高的前12個SNP標(biāo)記(AR2、AR5、AR7、TT8、TT11、PRO19、AT21、AT23和SER32、36、43、64)進(jìn)行組合時, 可以鑒別98份大豆種質(zhì), 每個SNP的平均鑒別效率最高, 達(dá)到8.17, 即該12個SNP標(biāo)記為最佳組合。利用上述12個SNP標(biāo)記與具有特異等位變異的2個SNP標(biāo)記(TT12和AT22)相結(jié)合組成GlySNP14, 在599份大豆種質(zhì)中共形成176個單倍型, 其中100個為特異單倍型。即利用GlySNP14可以將100份大豆種質(zhì)完全區(qū)分, 同時也可以將該100份大豆種質(zhì)從其他種質(zhì)種中區(qū)分出來。

表1 SNP標(biāo)記多樣性分析

1TV: 顛換;2TS: 轉(zhuǎn)換;3D: 缺失;4ND: 數(shù)據(jù)缺失。

1TV: transversion;2TS: transition;3D: deletion;4ND: no data recorded.

表2 大豆3個主要生態(tài)區(qū)種質(zhì)資源的遺傳相似性

圖2 基于SNP標(biāo)記的單倍型分析

A: 599份大豆品種的單倍型分析; B: 特異單倍型分析。

A: Analysis of haplotype for 599 soybean varieties; B: Analysis of specific haplotype.

圖3 SNP標(biāo)記組合分析

2.3 GlySNP14的聚類分析

GlySNP14可以將100份大豆種質(zhì)完全區(qū)分開(圖4), 共分為8組, A與F組分別以地方和選育品種為主, 各占56.25%和76.02%; B組的3個品種分別是1個美國春大豆和2個地方品種; C、G和H組均由國外和選育品種組成; E組10個品種中3個為美國春大豆, 7個為黃淮夏大豆; 最后D組為混合組, 20個品種中40%為東北春大豆, 國外品種和南方春大豆各占30%。

2.4 GlySNP14鑒別能力的模擬分析

SNP大多只有兩種堿基形式, 被視為二等位標(biāo)記[12]。雖然14個SNP理論上有214種組合方式, 但在現(xiàn)實種質(zhì)或群體材料中形成的組合方式要比理論上低很多。利用每個SNP在599份種質(zhì)中形成的等位變異頻率, 對選出的14個SNP的品種鑒別能力的模擬分析表明, 這14個SNP最多可鑒別750個不同種質(zhì)。為了驗證GlySNP14的鑒別效率, 從20個SNP中(除去3個沒有堿基差異的標(biāo)記TT14、PRO16和SER46)隨機(jī)選出14個SNP進(jìn)行鑒別能力模擬分析, 該過程模擬1000次。結(jié)果如圖5所示, 隨機(jī)選擇的任意14個SNP組合, 最多只能鑒別361個品種。在8個不同大小的模擬群體中, GlySNP14的鑒別能力均強(qiáng)于隨機(jī)選擇的SNP標(biāo)記。

2.5 SNP指紋圖譜構(gòu)建

為構(gòu)建身份證, 對大豆種質(zhì)的SNP數(shù)據(jù)進(jìn)行數(shù)字化編碼。在14個SNP標(biāo)記中共有10種基因型, 分別為AA、GG、CC、TT、DELDEL、TA、CT、CG、GA和C.DEL, 用數(shù)字1~5對堿基A、G、C、T和DEL編碼。例如, 地方品種“茶色豆”的14個SNP標(biāo)記基因型分別為GG、AA、AA、CC、TT、GG、TT、CC、TT、GG、CC、GG、TT和AA, 轉(zhuǎn)換成28位的指紋編碼為2211113344224433442233224411, 其中第1~4位的“2211”表示前2個SNP位點AR2和AR5基因型為GG和AA, 第5~8位的“1133”表示位列第3、第4位的SNP位點基因型為AA與CC, 其余20位的編碼以此類推。

圖4 基于SNP標(biāo)記的100份大豆品種聚類分析圖

圖5 隨機(jī)選擇的14個SNP與GlySNP鑒別能力的模擬分析

2.6 大豆種質(zhì)商品碼構(gòu)建

除SNP指紋編碼外, 進(jìn)一步對大豆品種的商品信息編碼。基本商品信息分為3個部分: (1)作物及品種屬性代碼。利用7位數(shù)字標(biāo)識作物及品種屬性, 其中1~6位表示大豆種屬(大豆屬于大田作物中的豆科作物); 第7位表示大豆栽培、野生類型, 即“1”表示栽培大豆, “2”表示野生大豆。(2)區(qū)域代碼。用于表示大豆品種的來源地區(qū), 以各省市的行政區(qū)劃代碼組成, 即河北為13、山東為37、黑龍江為23等, 國外品種以00表示。(3)品種類別代碼。用“1、2”表示大豆的地方與選育兩種類型。例如, 大豆品種“青仁黑豆”的商品碼為“0102011221”, 其中前7位“0102011”為作物及品種屬性代碼, 其中第1、第2位的“01”表示大田作物; 第3、第4位的“02”表示豆科; 第5、第6位的“01”表示大豆; 第7位的“1”表示栽培豆。第8、第9位的“22”是行政區(qū)劃代碼, 表示該品種來源地為吉林省。第10位“1”表示該品種為地方品種。

2.7 大豆種質(zhì)身份證構(gòu)建

大豆種質(zhì)身份證由上述的商品碼與指紋碼構(gòu)成, 總數(shù)為38位。以“北豐14”為例, 其身份證號為01020112322211115541223444111133224141, 其作物屬性為大田作物-豆科-大豆-栽培種(01-02-01-1); 品種來源為黑龍江選育品種(23-2); 品種DNA指紋為2211115541223444111133224141, 表示其14個SNP分子標(biāo)記基因型分別為GG、AA、AA、DELDEL、TA、GG、CT、TT、AA、AA、CC、GG、TA和TA。依據(jù)此方法完成了其他大豆種質(zhì)身份證的構(gòu)建。利用條碼在線生成軟件分別生成大豆品種身份證的一維碼與二維碼。

圖6 “北豐14”品種身份證及條形碼示意圖

A: “北豐14”品種身份證構(gòu)成; B: “北豐14”品種身份證條形碼。

A: Composition of “Beifeng 14” variety ID; B: Bar code of “Beifeng 14” variety ID.

利用條碼在線生成軟件分別生成100份大豆品種身份證的一維碼與二維碼?！氨必S14”的身份證條碼如圖6-B, 其余99份品種的部分身份證條碼見表3。將本研究中大豆品種身份證以一維碼或二維碼的形式來標(biāo)示相應(yīng)的種子, 就可以對種子的相關(guān)信息快速識別和規(guī)范管理, 也對相應(yīng)品種的知識產(chǎn)權(quán)提供了保護(hù)和科學(xué)依據(jù)。

3 討論

3.1 品種DNA指紋圖譜的不同表現(xiàn)形式

利用分子標(biāo)記構(gòu)建DNA指紋圖譜, 已逐漸向快速高效監(jiān)測方向發(fā)展, 并趨向于數(shù)字化。趙洪錕等[16]構(gòu)建吉林大豆骨干親本及主推品種RAPD標(biāo)記指紋圖譜, 采用的是分子標(biāo)記在瓊脂糖凝膠電泳中的遷移率大小及條帶的有無。隨后, 從定性描述轉(zhuǎn)為二位碼, 將分子標(biāo)記條帶的有無轉(zhuǎn)化為“0”和“1”編碼成數(shù)字指紋[17]。但上述二位碼大多屬于單純的遺傳變異編碼, 并未與品種信息相結(jié)合。陸徐忠等[18]利用SSR標(biāo)記鑒定水稻品種, 并將SSR標(biāo)記信息與商品信息相結(jié)合構(gòu)建了供試水稻品種身份證, 并最終以條形碼的形式表現(xiàn)。該方法構(gòu)建的品種身份證表示形式簡便, 易于監(jiān)測。但前人所使用的SSR標(biāo)記多為隨機(jī)標(biāo)記, 無法與表型性狀相聯(lián)系。相比之下, SNP具有針對性強(qiáng), 變異來源豐富, 且潛在數(shù)量巨大等特點。然而, 利用SNP標(biāo)記, 尤其是像本研究利用功能標(biāo)記構(gòu)建作物身份證, 并最終以條形碼的形式表現(xiàn)的研究尚未見報道。

表3 部分大豆品種的身份證條碼信息

本研究構(gòu)建的大豆種質(zhì)資源身份證, 不同于遺傳上單純的分子指紋圖譜。該身份證由38個數(shù)字組成, 包括品種屬性碼和分子指紋碼, 其中品種屬性碼由作物屬性、品種來源、品種類別3個部分組成, 最終將品種身份證編碼轉(zhuǎn)換為條碼形式(一維碼和二維碼)。利用該方法構(gòu)建的大豆品種身份證可以應(yīng)用在多種途徑上, 例如市場上商品種子的管理, 將品種身份證條形碼或二維碼標(biāo)識于種子包裝上, 可以方便種子信息的獲得與防偽。同時本試驗構(gòu)建身份證用的分子標(biāo)記是主要農(nóng)藝性狀功能基因SNP標(biāo)記, 包括結(jié)莢習(xí)性、茸毛色、熟期、百粒重、胞囊線蟲病與花葉病毒病抗性。上述農(nóng)藝性狀均與大豆品質(zhì)和農(nóng)業(yè)生產(chǎn)及農(nóng)民利益息息相關(guān), 有望用于進(jìn)一步豐富品種分子身份證信息。

3.2 SNP標(biāo)記與SSR標(biāo)記指紋圖譜的比較

指紋圖譜構(gòu)建所用標(biāo)記的選擇向快速、成本低和高通量檢測方向發(fā)展。國際植物品種權(quán)保護(hù)聯(lián)盟(UPOV)確定構(gòu)建DNA指紋數(shù)據(jù)庫的標(biāo)記包括SSR和SNP[9]。SSR具有多態(tài)性高的優(yōu)點, 但相比之下, SNP標(biāo)記在基因組中更為豐富。Jones等[19]同時利用88個SSR標(biāo)記和187個SNP標(biāo)記對58個玉米自交系和2個雜交種進(jìn)行基因型鑒定分析表明, SSR標(biāo)記的平均等位變異數(shù)目是SNP標(biāo)記的2.5倍, 但平均雜合率比SNP高0.5倍; 相比之下, SNP的可重復(fù)性(98.7%)顯著高于SSR標(biāo)記(91.7%)。張昆鵬等[20]比較9個油菜品種的SSR指紋圖譜和SNP指紋圖譜發(fā)現(xiàn), 基于SNP標(biāo)記計算的遺傳系數(shù)更精確, 結(jié)果更可靠。這表明SNP標(biāo)記的高重復(fù)性和準(zhǔn)確性能夠保證不同條件下檢測結(jié)果的正確性和可比性, 這對于種質(zhì)和品種的鑒定、品種權(quán)的保護(hù)與管理都是非常重要的技術(shù)支撐。

大豆身份證構(gòu)建的相關(guān)研究中, SSR標(biāo)記是鑒別品種的常用方法。但相比之下, SNP在大豆中分布更為廣泛。到目前為止, SNP在栽培大豆中的總數(shù)達(dá)到4~5×106個[21], 野生大豆PAN-GENOME中存在3.63~4.72×106個[23], 幾乎每1 kbp就可發(fā)現(xiàn)4個SNP。雖然與SSR相比, SNP等位變異類型有限, 但這一不足可以通過提高SNP數(shù)量來克服; 由于SNP大多只有兩種堿基形式, 被視為二等位標(biāo)記, 具有簡化基因分型方法及后續(xù)數(shù)字編碼的優(yōu)勢。與SSR標(biāo)記比, SNP是核苷酸自身變異, 不需要讀取分子量大小, 采用測序方法直接獲得等位變異基因型, 遺傳穩(wěn)定性高, 更適合于數(shù)字化建庫。

最少SNP標(biāo)記的選擇是降低構(gòu)建指紋圖譜成本的核心。與SSR的等位基因數(shù)多相比, SNP需要的數(shù)量比較多。Yoon等[6]比較了23個SNP核心標(biāo)記與Song等[22]開發(fā)的13個SSR標(biāo)記鑒定大豆種質(zhì)的效率, 結(jié)果表明, 平均每個SSR位點的等位變異數(shù)為7.80, 平均多樣性為0.73, 可以將101個品種完全鑒別; 每個SNP位點的平均等位變異數(shù)為2, 平均多樣性為0.45, 可鑒別132份供試品種; 而對23個SNP位點鑒別能力模擬分析結(jié)果表明, 最多能鑒別2200份種質(zhì)。相比之下, 本研究選用的14個標(biāo)記模擬分析中最多可以鑒別750份種質(zhì), 但GlySNP14標(biāo)記組合在實際的599份大豆種質(zhì)中鑒別出100份種質(zhì)。這可能是由于選用的標(biāo)記為重要性狀相關(guān)標(biāo)記, 有些標(biāo)記為功能基因標(biāo)記, 多態(tài)性相對較低有關(guān), 但這些標(biāo)記在鑒別品種的同時, 具有提供基因型數(shù)據(jù)的優(yōu)點。這為今后作物種質(zhì)SNP 指紋圖譜庫的構(gòu)建提供了有益的啟示。

4 結(jié)論

從23個功能基因相關(guān)SNP中挑選出14個, 精準(zhǔn)鑒定了大豆種質(zhì)品種。GlySNP14在599份大豆品種中鑒定出100份種質(zhì)并對其構(gòu)建了品種身份證, 該身份證由品種的屬性信息與分子指紋信息兩部分組成, 并最終以條碼形式(一維碼和二維碼)表示, 為大豆種質(zhì)資源管理與保護(hù)利用提供了便利。

[1] Fukuda Y. Cytogenetical studies on the wild and cultivated Manchurian soybeans (L.)., 1933, 6: 489–506

[2] Hymowitz T, Newell C A. Taxonomy of the genus Glycine, domestication and uses of soybeans., 1981, 35(3): 272–288

[3] 宋喜娥, 李英慧, 常汝鎮(zhèn), 郭平毅, 邱麗娟. 中國栽培大豆((L.) Merr.)微核心種質(zhì)的群體結(jié)構(gòu)與遺傳多樣性. 中國農(nóng)業(yè)科學(xué), 2010, 43: 2209–2219 Song X E, Li Y H, Chang R Z, Guo P Y, Qiu L J. Population structure and genetic diversity of mini core collection of cultivated soybean ((L.) Merr.) in China., 2010, 43: 2209–2219 (in Chinese with English abstract)

[4] 趙團(tuán)結(jié), 蓋鈞鎰. 栽培大豆起源與演化研究進(jìn)展. 中國農(nóng)業(yè)科學(xué), 2004, 37: 954–962 Zhao T J, Gai J Y. The origin and evolution of cultivated soybean [(L.) Merr.]., 2003, 37: 954–962 (in Chinese with English abstract)

[5] 王彩潔, 孫石, 金素娟, 李偉, 吳存祥, 侯文勝, 韓天富. 中國大豆主產(chǎn)區(qū)不同年代大面積種植品種的遺傳多樣性分析. 作物學(xué)報, 2013, 39: 1917–1926 Wang C J, Sun S, Jin S J, Li W, Wu T X, Hou W S, Han T F. Genetic Diversity analysis of widely-planted soybean varieties from different decades and major production regions in China., 2013, 39: 1917–1926 (in Chinese with English abstract)

[6] Yoon M S, Song Q J, Choi I Y, Specht J E, Hyten D, Cregan P. BARCSoySNP23: a panel of 23 selected SNPs for soybean cultivar identification., 2007, 114: 885–899

[7] 宋婉, 續(xù)九如. 果樹種質(zhì)資源鑒定及 DNA 指紋圖譜應(yīng)用研究進(jìn)展. 北京林業(yè)大學(xué)學(xué)報, 2000, 22(1): 76–80 Song W, Xu J R. Progress in application of germplasm identification and fingerprinting of fruit trees., 2000, 22(1): 76–80 (in Chinese with English abstract)

[8] 邱福林, 莊杰云, 華澤田, 王彥榮, 程式華. 北方雜交粳稻骨干親本遺傳差異的SSR標(biāo)記檢測. 中國水稻科學(xué), 2005, 19: 101–104 Qiu F L, Zhuang J Y, Hua Z T, Wang Y R, Cheng S H. Inspect of genetic differentiation of main parents of japonica hybrid rice in the Northern China by Simple Sequence Repeats (SSR)., 2005, 19: 101–104 (in Chinese with English abstract)

[9] 滕海濤, 呂波, 趙久然, 徐巖, 王鳳格, 堵苑苑, 楊坤, 唐浩, 李祥羽. 利用 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)

[10] Jung J, Park S, Liu W Y, Kang B. Discovery of single nucleotide polymorphism in Capsicum and SNP markers for cultivar identification., 2010, 175: 91–107

[11] Shirasawa K, Monna L, Kishitani S, Nishio T. Single nucleotide polymorphisms in randomly selected genes among japonica rice (L.) varieties identified by PCR-RF-SSCP., 2004, 11: 275–283

[12] 王大莉. 香菇栽培品種SNP指紋圖譜庫的構(gòu)建. 華中農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 湖北武漢, 2012 Wang D L. Construction of a SNP Fingerprint Databases forCultivars. MS Thesis of Huazhong Agricultural University, Wuhan, China, 2012 (in Chinese with English abstract)

[13] Liu, K, Muse S V. PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics, 2005, 21: 2128–2129

[14] Poole R W. Introduction to quantitative ecology. McGraw-Hill, New York, 1974

[15] Milne I, Shaw P, Stephen G, Bayer M, Cardle L, Thomas W T, Flavell A J, Flavell M D, Marshall D. Flapjack-graphical genotype visualization., 2010, 26: 3133–3134

[16] 趙洪錕, 李啟云. 吉林省大豆骨干親本及主推品種 DNA 指紋圖譜的構(gòu)建. 中國油料作物學(xué)報, 2000, 22(4): 12–14 Zhao H K, Li Q Y. Jilin province soybean backbone parents and construction of main varieties of DNA fingerprint., 2000, 22(4): 12–14 (in Chinese with English abstract)

[17] 匡猛, 楊偉華, 許紅霞, 王延琴, 周大云, 馮新愛. 中國棉花主栽品種DNA指紋圖譜構(gòu)建及SSR標(biāo)記遺傳多樣性分析. 中國農(nóng)業(yè)科學(xué), 2011, 44: 20–27 Kuang M, Yang W H, Xu H X, Wang Y Q, Zhou D Y, Feng X A. Construction of DNA fingerprinting and analysis of genetic diversity with SSR markers for cotton major cultivars in China.,2011, 44: 20–27 (in Chinese with English abstract)

[18] 陸徐忠, 倪金龍, 李莉, 汪秀峰, 馬卉, 張小娟, 楊劍波. 利用SSR分子指紋和商品信息構(gòu)建水稻品種身份證. 作物學(xué)報, 2014, 40: 823–829 Lu X Z, Ni J L, Li L, Wang X F, Ma H, Zhang X J, Yang J B. Construction of rice variety indentity using SSR fingerprint and commodity information., 2014, 40: 823–829 (in Chinese with English abstract)

[19] Jones E S, Sullivan H, Bhattramakki D, Smith J. A comparison of simple sequence repeat and single nucleotide polymorphism marker technologies for the genotypic analysis of maize (L.)., 2007, 115: 361–371

[20] 張昆鵬. 利用SNP標(biāo)記構(gòu)建油菜品種指紋圖譜及定位下卷葉性狀基因的研究. 南京農(nóng)業(yè)大學(xué)博士學(xué)位論文, 江蘇南京, 2013 Zhang K P. Studies on Rapeseed Variety Fingerprints and Mapping of Gene of the Down-curly Leaf by Use of SNP Markers inL. PhD Dissertation of Nanjing Agricultural University, Nanjing, China, 2013 (in Chinese with English abstract)

[21] Song Q J, Quigley C V, Nelson R L, Carter T, Boerma H, Strachan J, Cregan P. A selected set of trinucleotide simple sequence repeat markers for soybean cultivar identification., 1999, 12: 207–220

[22] Zhu Y L, Song Q J, Hyten D L, Van Tassell C, Matukumalli L, Grimm D, Hyatt S, Fickus E, Young N, Cregan P. Single- nucleotide polymorphisms in soybean., 2003, 163: 1123–1134

[23] Li Y, Zhou G, Ma J, Jiang W, Jin L, Zhang Z, Guo Y, Zhang J, Sui Y, Zheng L, Zhang S, Zuo Q, Shi X, Li Y, Zhang W, Hu Y, Kong G, Hong H, Tan B, Song J, Liu Z, Wang Y, Ruan H, Yeung C K L, Liu J, Wang H, Zhang L, Guan R, Wang K, Li W, Chen S, Chang R, Jiang Z, Jackson S A, Li R, Qiu L. De novo assembly of soybean wild relatives for pan-genome analysis of diversity and agronomic traits., 2014, 32: 1045–1052

Accurate Identification of Varieties by Nucleotide Polymorphisms and Establishment of Scannable Variety IDs for Soybean Germplasm

WEI Zhong-Yan1, LI Hui-Hui1, LI Jun2, Yasir A. Gamar1, MA Yan-Song3, and QIU Li-Juan1, *

1National Key Facility for Crop Gene Resources and Genetic Improvement (NFCRI) / Key Laboratory of Germplasm Utilization, Ministry of Agriculture / Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China;2College of Agriculture, Guangxi University, Nanjing 530004, Guangxi China;3Soybean Research Institute, Heilongjiang Academy of Agricultural Sciences, Harbin 150086, Heilongjiang, China

In order to strengthen the management of soybean germplasm and variety protection, SNP markers were developed to establish the identity of soybean varieties. A set of 23 SNP markers distributed in 13 genes were used to discriminate genotypes of 599 soybean varieties, grown in most of the soybean producing areas in China. Fourteen SNPs with high polymorphism selected from the 23 SNPs (GlySNP14) showed the improved variety identification capability, compared with any combination of 14 random SNPs. A simulated experiment confirmed that GlySNP14 could effectively distinguish of 750 soybean varieties, while the combination of random by selected from the 14 SNPs could distinguish only 361 varieties. The established ID of soybean varieties in this study contained a 38-digit serial code, which can be used for the accurate identification of soybean varieties and meet the requirements of genetic resources protection.

soybean; SNP; genetic diversity; DNA fingerprint; variety ID

2017-06-11;

2017-11-21;

2017-12-11.

10.3724/SP.J.1006.2018.000315

本研究由農(nóng)產(chǎn)品質(zhì)量安全監(jiān)管(種子管理)項目(2130109)和國家農(nóng)作物資源共享平臺(NICGR2016大豆)項目資助。

This study was supported by the project of Quality and Safety Supervision of agriculture food (Seed Management) (2130109) and National Infrastructure for Crop Germplasm Resources (NICGR2016 Soybean).

Corresponding author邱麗娟, E-mail: qiulijuan@caas.cn, Tel: 86-10-82105843, Fax: 86-10-82105843

E-mail: w_zhongyan@163.com

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