用于羊草基因分型的SNP分子標(biāo)記技術(shù)研究

侯莉娟，齊　曉，齊冬梅，陳雙燕*，劉公社*（.中國科學(xué)院植物研究所北方資源植物重點實驗室，北京00093；2.中國農(nóng)業(yè)科學(xué)院草原研究所，內(nèi)蒙古呼和浩特0000；3.全國畜牧總站，北京0025）

用于羊草基因分型的SNP分子標(biāo)記技術(shù)研究

侯莉娟1**，齊曉2，3**，齊冬梅1，陳雙燕1*，劉公社1*
（1.中國科學(xué)院植物研究所北方資源植物重點實驗室，北京100093；2.中國農(nóng)業(yè)科學(xué)院草原研究所，內(nèi)蒙古呼和浩特010010；3.全國畜牧總站，北京100125）

以我國廣泛分布的多年生優(yōu)勢禾本科牧草羊草為研究材料，針對羊草目前尚無SNP標(biāo)記技術(shù)的報道，基于實驗室已有的轉(zhuǎn)錄組數(shù)據(jù)，挖掘其中的SNP分子標(biāo)記，建立一種利用高分辨率熔解曲線（high-resolution melting，HRM）進(jìn)行羊草SNP基因分型的方法。結(jié)果表明，從羊草轉(zhuǎn)錄組數(shù)據(jù)中獲得41843個SNPs，轉(zhuǎn)換的數(shù)量是顛換數(shù)量的2倍；建立了基于SNP的HRM基因分型方法，設(shè)計了112對SNP引物，其中有98對（87.5%）能擴(kuò)增出明亮單一條帶，擴(kuò)增條帶大小符合預(yù)期的有72對（64.29%），適宜進(jìn)行HRM基因分型的引物有46對，占41.07%；用其中的7對引物可將羊草48份種質(zhì)區(qū)分開，并繪制出指紋圖譜。通過對部分HRM基因分型結(jié)果進(jìn)行測序驗證，表明測序分析的基因型結(jié)果與HRM基因分型結(jié)果一致，證明HRM基因分型結(jié)果是可靠的。本研究首次建立的基于SNP的HRM基因分型方法，可用于羊草種質(zhì)資源的指紋圖譜繪制、育種親本的選配、基因關(guān)聯(lián)分析、分子標(biāo)記輔助育種、功能分子標(biāo)記開發(fā)、新品種的保護(hù)等育種工作。

羊草；轉(zhuǎn)錄組；SNP標(biāo)記；HRM基因分型；指紋圖譜

http：//cyxb.lzu.edu.cn

侯莉娟，齊曉，齊冬梅，陳雙燕，劉公社.用于羊草基因分型的SNP分子標(biāo)記技術(shù)研究.草業(yè)學(xué)報，2016，25（2）：105-113.

HOU Li-Juan，QI Xiao，QI Dong-Mei，CHEN Shuang-Yan，LIU Gong-She.A study of SNP molecular marker technology for Leymuschinensis genotyping.Acta Prataculturae Sinica，2016，25（2）：105-113.

單核苷酸多態(tài)性（single nucleotide polymorphisms，SNPs）指在個體基因組上任何一個核苷酸的變異形成的DNA序列多態(tài)性［1］。與其他分子標(biāo)記（如SSR）相比［2］，SNP分子標(biāo)記具有明顯的優(yōu)勢。SNPs數(shù)量多，分布廣泛，多態(tài)性豐富，密度高。在植物基因組內(nèi)，平均每21 bp就出現(xiàn)一個SNP。一般只有兩種堿基組成，是二等位基因，只需+/-的分析，適于快速、規(guī)?；Y查基因組。某些位于基因內(nèi)部的SNP有可能直接影響蛋白質(zhì)結(jié)構(gòu)或表達(dá)水平。SNP等位基因頻率容易估計，易于基因分型，遺傳穩(wěn)定性比SSR高，能夠進(jìn)行高通量自動化基因分型檢測［3-5］。因此，SNP被認(rèn)為是應(yīng)用前景最好的遺傳標(biāo)記之一。

隨著新一代測序技術(shù)的發(fā)展，轉(zhuǎn)錄組測序已經(jīng)成為鑒定非模式植物SNP多態(tài)性分子標(biāo)記的重要數(shù)據(jù)源?；谵D(zhuǎn)錄組挖掘SNP分子標(biāo)記具有顯著的優(yōu)點，因為這些SNP變異位于功能基因上，可能直接與植物表型相關(guān)。目前基于轉(zhuǎn)錄組測序開發(fā)SNP標(biāo)記已在水稻（Oryza sativa）、玉米（Zea mays）、高粱（Sorghum bicolor）、大麥（Hordeum vulgare）、小麥（Triticum aestivum）等作物以及苜蓿（Medicago sativa）、高羊茅（Festuca arundinacea）和白三葉（Trifolium repens）等牧草中廣泛開展，開發(fā)的SNP被廣泛應(yīng)用于遺傳圖譜構(gòu)建和遺傳多樣性分析中［6-13］。高分辨率熔解曲線（high-resolution melting，HRM）是在PCR基礎(chǔ)上通過測定DNA雙鏈熔解曲線變化來檢測SNP變異的方法。HRM操作簡便，在PCR結(jié)束后添加LC Green染料，通過監(jiān)測升溫過程中熒光染料與PCR擴(kuò)增產(chǎn)物的結(jié)合情況來實現(xiàn)的，在升溫過程中，雙鏈解開，熒光染料從解鏈的分子上釋放，同時熒光強度降低，所以通過熒光強度和曲線的變化上就可以進(jìn)行SNP基因分型。該方法高度靈敏、簡單有效［14］，已經(jīng)在多倍體苜蓿、多年生黑麥草（Lolium perenne）、土豆（Solanum tuberosum）等物種中應(yīng)用，進(jìn)行基因分型和定位性狀［5，15-16］。

羊草（Leymus chinensis）是我國廣泛分布的具有優(yōu)勢的多年生禾本科牧草，目前在種質(zhì)資源的精細(xì)評價和創(chuàng)新利用方面還比較薄弱，亟須采用新一代SNP分子標(biāo)記技術(shù)對這些豐富的資源進(jìn)行精細(xì)評價，以便發(fā)掘其中有價值的資源，加速牧草育種工作進(jìn)程。針對羊草目前尚無單核苷酸多態(tài)性SNP標(biāo)記技術(shù)的報道，本研究基于實驗室建立的羊草轉(zhuǎn)錄組數(shù)據(jù)，挖掘其中的SNP分子標(biāo)記，首次建立了一種基于SNP標(biāo)記的HRM分析方法，用于羊草種質(zhì)資源的基因分型和精細(xì)評價。

1　材料與方法

1.1材料

本研究所用到的羊草材料取自中國科學(xué)院植物研究所羊草種質(zhì)資源圃，于2014-2015年進(jìn)行實驗。

1.2藥品和試劑

廣譜植物基因組DNA快速提取試劑盒購自北京博邁德基因技術(shù)有限公司；PCR擴(kuò)增kit購自北京全式金公司；LC Green飽和熒光染料（美國Idaho公司，貨號BCH M-ASY-0005），購自思博全科技有限公司；Light Scanner專用礦物油為Sigma公司產(chǎn)品；p MD19-T simple載體購自Ta KaRa公司。SanPrep柱式DNA膠回收試劑盒購自生工生物工程上海（股份）有限公司；引物合成和克隆測序由北京諾賽基因組研究中心有限公司完成；其他所用普通化學(xué)試劑均為國產(chǎn)分析純。

1.3方法

1.3.1發(fā)掘羊草SNP根據(jù)實驗室前期獲得的羊草轉(zhuǎn)錄組數(shù)據(jù)，包含了豐富的組織器官和各種脅迫處理的材料［17-20］。采用生物信息學(xué)技術(shù)發(fā)掘羊草SNP，首先使用BWA軟件將所有Clean reads定位到87214條Uni-gene上，然后利用samtools去call SNP，BWA和samtools均使用軟件推薦的參數(shù)［21］。

1.3.2設(shè)計SNP引物根據(jù)生物信息學(xué)分析發(fā)掘的SNP位點，利用DNAMAN軟件設(shè)計SNP引物。具體方法如下：首先確定SNP位點所在序列的方向，如果序列是反向的，需要進(jìn)行反向互補，然后與水稻基因組（http：//rice.plantbiology.msu.edu/）進(jìn)行BLASTN比較，確認(rèn)序列的結(jié)構(gòu)，避免擴(kuò)增片段跨越兩個外顯子。在SNP位點的上下游分別設(shè)計正向引物和反向引物，引物設(shè)計的參數(shù)為引物長度18～24 bp；Tm 62～66℃，盡可能保證上下游引物的Tm值一致，一般不超過2℃；GC含量40%～60%，45%～55%最佳；產(chǎn)物大小為50～200 bp。

1.3.3建立高分辨率熔解曲線方法該方法的建立包括以下條件的摸索：篩選有效擴(kuò)增的SNP引物對，確定HRM-PCR反應(yīng)體系與程序，確定適宜的DNA模板，分析LC Green染料是否對PCR擴(kuò)增有影響，優(yōu)化HRM熔解曲線分析。具體步驟為：首先在PCR擴(kuò)增儀（Bio-Rad T100）上進(jìn)行PCR反應(yīng)，之后在反應(yīng)體系內(nèi)加入LC Green飽和熒光染料，并在PCR儀內(nèi)進(jìn)行1 min短暫的錯配反應(yīng)，然后將PCR產(chǎn)物（96/孔板）插入LightScanner儀器上進(jìn)行HRM分析。

1.3.4HRM用于羊草種質(zhì)基因分型利用建立的HRM方法對48份羊草種質(zhì)材料進(jìn)行基因分型。具體步驟如下：以羊草基因組DNA為模板，采用待檢測SNP引物進(jìn)行PCR擴(kuò)增。PCR反應(yīng)體系10μL：40 ng模板DNA、1.5 mmol/L Mg2+、400μmol/L d NTPs、0.5 U TaqDNA聚合酶、1μmol/L引物及10×PCR buffer。PCR反應(yīng)程序：94℃預(yù)變性2 min；然后94℃變性30 s，72℃復(fù)性30 s，72℃延伸30 s（復(fù)性溫度從72℃到58℃每個循環(huán)降落1℃），然后進(jìn)行30個循環(huán)94℃變性30 s，58℃復(fù)性30 s，72℃延伸30 s，最后72℃延伸5 min。之后在PCR擴(kuò)增產(chǎn)物中加入1μL LC Green，然后進(jìn)行95℃30 s，25℃30 s。此程序結(jié)束后，將反應(yīng)產(chǎn)物放在Light-Scanner儀器上進(jìn)行HRM分析，對羊草不同種質(zhì)進(jìn)行基因分型。

1.3.5測序驗證HRM分型結(jié)果選擇HRM熔解峰和熔解曲線較理想的部分SNP引物對擴(kuò)增產(chǎn)物進(jìn)行測序，以驗證HRM分型結(jié)果是否可靠。具體步驟如下：將待測序SNP引物的HRM不同基因分型的PCR產(chǎn)物連接到p MD19-T載體上，轉(zhuǎn)化大腸桿菌DH5α，利用ABI 3730測序儀進(jìn)行測序，根據(jù)測序峰圖和等位基因出現(xiàn)的頻率確定基因型，并與HRM基因分型結(jié)果進(jìn)行比較，以驗證引物基因分型的準(zhǔn)確性。

表1　羊草SNP類型及引物設(shè)計情況Table1　Various types of SNP and designed primers in L.chinensis

2　結(jié)果與分析

2.1從羊草轉(zhuǎn)錄組數(shù)據(jù)中發(fā)掘SNP

通過對羊草轉(zhuǎn)錄組數(shù)據(jù)進(jìn)行生物信息學(xué)分析，如表1所示，共獲得41781個SNPs，其中堿基顛換（transversion，即嘌呤與嘧啶之間的替代）有8種形式（AC/AT/CA/CG/GC/GT/TA/TG），共計13053個；轉(zhuǎn)換（transition，即嘌呤與嘌呤，或嘧啶與嘧啶之間的替換）有4種形式（AG/GA/CT/TC），共計28728個；轉(zhuǎn)換的數(shù)量約是顛換的2倍，符合一般規(guī)律。在8種形式顛換和4種形式轉(zhuǎn)換SNP變異中各選取一些SNP位點設(shè)計引物，共設(shè)計了112對引物，包含112 個SNP位點。

2.2建立HRM分析方法

2.2.1PCR體系和擴(kuò)增程序HRM基因分型對模板DNA比較敏感，為了保證模板DNA的均質(zhì)化，本實驗所有DNA樣品采用統(tǒng)一試劑盒提取和純化，然后將所有DNA的濃度調(diào)至統(tǒng)一的標(biāo)準(zhǔn)濃度（本研究以濃度40 ng/μL為標(biāo)準(zhǔn)）。結(jié)果發(fā)現(xiàn)DNA濃度在40～120 ng/μL的范圍內(nèi)都能擴(kuò)增出清晰明亮的條帶。為了節(jié)省DNA量，達(dá)到多次使用的目的，將PCR反應(yīng)體系中DNA模板濃度設(shè)定為40 ng/μL。PCR反應(yīng)體系中引物濃度為1μmol/L，與常規(guī)PCR要求的一致，但是對于個別引物，擴(kuò)增產(chǎn)物中如果有引物二聚體或非特異性擴(kuò)增條帶，則將引物濃度降低至0.25μmol/L，可降低引物二聚體。

由于各引物對的Tm值和適宜退火溫度不一致，摸索每對引物的適宜退火溫度比較繁瑣，同時為了提高實驗效率，需要確定一個廣譜穩(wěn)定的PCR擴(kuò)增程序，適合所有引物對的擴(kuò)增。本研究經(jīng)過摸索，設(shè)定了15個循環(huán)的降落PCR，退火溫度每個循環(huán)降落1℃，從72℃降落到58℃，之后在58℃的退火溫度下再進(jìn)行30個循環(huán)的擴(kuò)增比較合適（圖1），擴(kuò)增產(chǎn)物清晰明亮，提高了引物的擴(kuò)增效率。

2.2.2SNP引物對的篩選為了驗證所合成的SNP引物對能否擴(kuò)增出目的條帶，對其進(jìn)行了PCR擴(kuò)增及電泳檢測，結(jié)果如表2所示。112對SNP引物中，能夠擴(kuò)增出清晰明亮單一條帶的引物有98對，占87.5%；其中，擴(kuò)增條帶大小符合預(yù)期的引物有72對，占64.29%。然后對72對引物進(jìn)行HRM分析，以便篩選出能夠用于HRM基因分型的SNP引物對。通過觀察HRM的熔解曲線和熔解峰，發(fā)現(xiàn)HRM熔解曲線有明顯拐角，熔解峰高聳、單一的引物可以用于HRM分析（圖2），共挑選出46對引物可用于HRM基因分型（表2）。

圖1　循環(huán)數(shù)對PCR擴(kuò)增產(chǎn)物的影響Fig.1　Effect of cycle number to PCR products

圖2　可用于HRM分析的SNP引物的熔解曲線和熔解峰示例Fig.2　Example of melting curves and melting peaks of SNP primer pairs for HRM analysis

表2　SNP引物對的篩選Table2　Screening of SNP primer pairs

2.2.3LC Green對PCR擴(kuò)增的影響為了檢驗LC Green染料對PCR擴(kuò)增是否有影響，以確定加入染料的時間點，本研究設(shè)置了4個PCR反應(yīng)，即不加模板對照（A）不加染料（B）、PCR反應(yīng)之前加染料（C）和PCR反應(yīng)之后加染料（D），4個PCR反應(yīng)同時進(jìn)行，只是在染料添加的順序上有區(qū)別。結(jié)果表明（圖3），“B”管和“D”管的電泳結(jié)果均可看到清晰明亮的條帶，而PCR反應(yīng)之前加染料“C”管的條帶消失。說明PCR反應(yīng)前加入LC Green染料對PCR擴(kuò)增結(jié)果有影響。因此，LC Green應(yīng)該在PCR反應(yīng)結(jié)束之后加入，然后在PCR儀中進(jìn)行錯配反應(yīng)（95℃30 s；25℃30 s）。

2.2.4LightScanner檢測時熔解溫度范圍的設(shè)定

在運行HRM檢測前，需要選擇熔解溫度范圍，儀器默認(rèn)設(shè)置為65～95℃。但是對于某些引物對來說，如果將95℃設(shè)定為右邊界，熔解峰還沒有繪制完，反應(yīng)就結(jié)束了。對于這樣的引物對，需要增大熔解溫度右邊界的值，保證熔解峰的完整和分析的可靠性。儀器的最大溫度為98℃，因此在HRM分析時，可以將熔解溫度范圍設(shè)定的盡量大一些（設(shè)定為65～98℃），以保證所有引物對的熔解峰都能繪制完全。圖4所示為兩對代表性引物的熔解峰，A代表了絕大多數(shù)引物對的HRM熔解峰情況，在儀器默認(rèn)熔解溫度范圍內(nèi)，而B代表了個別引物對熔解峰出現(xiàn)較晚的情況。

圖3　LC Green對PCR擴(kuò)增產(chǎn)物的影響Fig.3　Effects of LC Green on PCR amplicons

圖4　兩對代表性引物對12個DNA進(jìn)行HRM分析的熔解峰Fig.4　Melting peaks of the two representative primer pairs for HRM analysis to 12 DNAs

2.3基于SNP的HRM對羊草種質(zhì)基因分型

利用上述篩選的46對HRM引物對，利用其中7對引物（代表基因分型少，4組；基因分型中等，5～7組；基因分型較多，8組）對羊草48個種質(zhì)進(jìn)行了基因分型，結(jié)果表明使用這7對引物就可以將羊草48個種質(zhì)區(qū)分開，繪制出指紋圖譜。如圖5所示，不同引物對的分型結(jié)果用不同的顏色區(qū)分；同一引物對的分型結(jié)果用數(shù)據(jù)條長短區(qū)分，數(shù)據(jù)條長度相同的種質(zhì)，基因型一致，數(shù)據(jù)條長度不同的種質(zhì)，表明在基因分型時被分到不同的組別。舉例來說，YC1和YC2兩份種質(zhì)，除了SNP42引物變異位點存在差異之外，其余6個SNP引物變異位點均沒有差異，表明這兩份種質(zhì)的差異較??；而YC22和YC23兩份種質(zhì)，7個SNP引物的變異位點均存在差異，表明它們的差異較大。

圖5　利用7對SNP引物對羊草48個種質(zhì)繪制的指紋圖譜Fig.5　Fingerprints of 48 L.chinensis DNA using 7 SNP primer pairs

2.4測序驗證HRM分型結(jié)果

根據(jù)上述HRM分型結(jié)果，挑選了3對代表性SNP引物，對其擴(kuò)增產(chǎn)物進(jìn)行測序以驗證HRM分型結(jié)果。從每對引物的每個分組（代表一個基因型）中挑選一個樣品進(jìn)行測序，通過對測序克隆峰圖的分析及克隆間的序列比對，統(tǒng)計SNP位點的堿基及其頻率，確定其對應(yīng)的基因型。測序結(jié)果發(fā)現(xiàn)（表3），除了目標(biāo)SNP位點的變異外，還出現(xiàn)了其他位點的SNP變異；對每對引物所擴(kuò)增產(chǎn)物的SNP變異位點基因型進(jìn)行分析，發(fā)現(xiàn)不同羊草種質(zhì)的基因型都不一致，這與HRM分型結(jié)果一致。以SNP2為例，HRM將羊草種質(zhì)YC3，YC5和YC12分型到a，b和c組中，而測序結(jié)果表明，3個種質(zhì)的目標(biāo)位點67基因型一致，但其他變異位點63和69的基因型有差異。因此，測序數(shù)據(jù)證明HRM分型結(jié)果是可靠的。

表3　3對引物擴(kuò)增產(chǎn)物SNP變異位點的基因型Table3　Genotypes of SNP mutation sites f rom the amplified products of three primer pairs

3　討論

近年來隨著高通量測序技術(shù)的飛速發(fā)展，各類生物信息學(xué)軟件的開發(fā)，使得SNP分子標(biāo)記的開發(fā)更便捷。對于多倍體或沒有參考基因組的物種來說，利用新一代高通量轉(zhuǎn)錄組測序技術(shù)可大量鑒定SNP［22］。目前，利用轉(zhuǎn)錄組測序技術(shù)已經(jīng)在水稻、玉米、高粱、大麥和小麥等農(nóng)作物以及苜蓿、高羊茅和白三葉等牧草中鑒定到大量的SNP［6-13］。本研究通過生物信息學(xué)分析從羊草轉(zhuǎn)錄組測序數(shù)據(jù)中挖掘出41843個SNPs，表明轉(zhuǎn)錄組測序是從多倍體非模式植物中獲得大量SNP的有效方法。

HRM分析是一種快速掃描PCR擴(kuò)增產(chǎn)物內(nèi)SNP變異的高度靈敏、簡單有效的方法［14］。該技術(shù)已經(jīng)成功應(yīng)用于同源四倍體植物苜蓿和土豆的SNP基因分型和多態(tài)性檢測，并且在發(fā)現(xiàn)變異、分析候選基因、比較和關(guān)聯(lián)定位的基因分型等育種工作中也有很好地應(yīng)用潛力［5，16］。本研究利用羊草轉(zhuǎn)錄組中發(fā)掘的SNP，在SNP引物的設(shè)計，HRM-PCR反應(yīng)體系與程序，LC Green染料對PCR擴(kuò)增的影響，HRM熔解曲線分析等方面進(jìn)行了摸索與優(yōu)化，建立了基于EST-SNP的羊草HRM分型方法，與前人研究相比［5，16］，本研究最大的特點是確定了一個廣譜穩(wěn)定的PCR擴(kuò)增程序，提高了實驗效率，并且利用HRM分析成功對羊草進(jìn)行了SNP基因分型，通過測序驗證HRM分型是可靠的。

大型機(jī)械壓力機(jī)大多為二到三級傳動結(jié)構(gòu)，第一級為電動機(jī)通過皮帶帶動主軸轉(zhuǎn)動，第二級為主軸上的小齒輪帶動曲軸上的大齒輪轉(zhuǎn)動。從而通過曲柄滑塊機(jī)構(gòu)驅(qū)動滑塊作上下往復(fù)的直線運動。因此大型機(jī)械壓力機(jī)傳動系統(tǒng)噪聲主要是齒輪副的噪聲。

羊草是異源四倍體，假設(shè)基因內(nèi)某一位點的堿基出現(xiàn)變異，例如C突變成T，則這一位點理論上有5種基因型：CCCC、CCCT、CCTT、CTTT、TTTT。因此對于有1個SNP變異位點的引物對來說，其擴(kuò)增產(chǎn)物在SNP位點的基因型理論上是5種。如果一對SNP引物有2個SNP變異位點（每個變異位點理論上有5種基因型），理論上能區(qū)分52=25個羊草基因型。本研究通過測序發(fā)現(xiàn)，SNP2、SNP14和SNP41引物對擴(kuò)增產(chǎn)物中，除了轉(zhuǎn)錄組數(shù)據(jù)中發(fā)掘的目標(biāo)SNP位點外，還有其他SNP變異位點，分別有3，4和7個SNP變異位點，因此，理論上可以區(qū)分53，54和57個羊草基因型。但是，由于受供檢測的羊草基因型的限制，以及HRM軟件最大分組的限制（最高分組為8組，超過8組的基因型被分在未知組中），HRM實際基因分型能力遠(yuǎn)遠(yuǎn)沒有達(dá)到理論水平。隨著HRM軟件的改進(jìn)，該技術(shù)在基因分型能力上可大大提高。

SNP在育種中的應(yīng)用非常廣泛，包括遺傳多樣性分析、等位基因挖掘、高分辨率遺傳連鎖圖譜構(gòu)建、QTL定位、標(biāo)記輔助選擇、品種鑒定等。本研究所獲得的SNP來源于羊草轉(zhuǎn)錄本，對應(yīng)基因，可利用這些SNP與田間表型性狀進(jìn)行關(guān)聯(lián)分析，發(fā)現(xiàn)與重要性狀關(guān)聯(lián)的SNP標(biāo)記。還可以用來進(jìn)行羊草親本選配、發(fā)展功能分子標(biāo)記、用于羊草分子標(biāo)記輔助育種、構(gòu)建指紋圖譜及用于羊草育成新品種的保護(hù)等，從而推進(jìn)羊草分子育種工作的進(jìn)程。

References：

［1］Primmer C R，Borge T，Lindell J，et al.Single-nucleotide polymorphism characterization in species with limited available sequence information：high nucleotide diversity revealed in the avian genome.Molecular Ecology，2002，11：603-612.

［2］Zhang L J，Liu L L，Chang Z J，et al.Development and functional verification of cSSR markers for recessive genetic male sterility in oats.Acta Prataculturae Sinica，2015，24（7）：146-154.

［3］Edward K，Poole R，Barker G.SNP discovery in plants.In：Henry R（ed）.Plant Genotyping II SNP Technology［M］.CABI Oxfordshire，2008：1-29.

［4］Mir R R，Hiremath P J，Lizarazu O R，et al.Evolving molecular marker technologies in plants：from RFLPs to GBS.Diagnostics in Plant Breeding，2013，11：229-247.

［5］Han Y，Khu D M，Monteros M J.High-resolution melting analysis for SNP genotyping and mapping in tetraploid alfalfa （Medicago sativa L.）.Molecular Breeding，2012，29（2）：489-501.

［7］Allen A M，Barker G L，Berry S T，et al.Transcript-specific，single-nucleotide polymorphism discovery and linkage analysis in hexaploid bread wheat（Triticum aestivum L.）.Plant Biotechnology Journal，2011，9（9）：1086-1099.

［8］Nelson J C，Wang S，Wu Y，et al.Single-nucleotide polymorphism discovery by high-throughput sequencing in sorghum. BMC Genomics，2011，12：352.

［9］Barbazuk W B，Schnable P S.SNP discovery by transcriptome pyrosequencing.Methods in Molecular Biology，2011，729：225-246.

［10］Close T J，Bhat P R，Lonardi S，et al.Development and implementation of high-throughput SNP genotyping in barley.BMC Genomics，2009，10：582.

［11］Li X，Acharya A，F(xiàn)armer A D，et al.Prevalence of single nucleotide polymorphism among 27 diverse alfalfa genotypes as assessed by transcriptome sequencing.BMC Genomics，2012，13：568.

［12］Hand M L，Cogan N O，F(xiàn)orster J W.Genome-wide SNP identification in multiple morphotypes of allohexaploid tall fescue （Festuca arundinacea Schreb）.BMC Genomics，2012，13：219.

［13］Nagy I，Barth S，Mehenni-Ciz J，et al.A hybrid next generation transcript sequencing-based approach to identify allelic and homeolog-specific single nucleotide polymorphisms in allotetraploid white clover.BMC Genomics，2013，14：100.

［14］Cho M H，Ciulla D，Klanderman B J，et al.High-resolution melting curve analysis of genomic and whole-genome amplified DNA.Clinical Chemistry，2008，54：2055-2058.

［15］Studer B，Jensen L，F(xiàn)iil A，et al.‘‘Blind''mapping of genic DNA sequence polymorphisms in Lolium perenne L.by high resolution melting curve analysis.Molecular Breeding，2009，24：191-199.

［16］de Koeyer D，Douglass K，Murphy A，et al.Application of high-resolution DNA melting for genotyping and variant scanning of diploid and autotetraploid potato.Molecular Breeding，2010，25：67-90.

［17］Chen S Y，Huang X，Yan X Q，et al.Transcriptome analysis in sheepgrass：a dominant perennial grass of the Eurasian Steppe.PLoS One，2013，8：e67974.

［18］Chen S Y，Cai Y Y，Zhang L X，et al.Transcriptome analysis reveals common and distinct mechanisms for sheepgrass （Leymus chinensis）responses to defoliation compared to mechanical wounding.PLoS One，2014，9：e89495.

［19］Huang X，Peng X J，Zhang L X，et al.Bovine serum albumin in saliva mediates grazing response in Ley mus chinensis revealed by RNA sequencing.BMC Genomics，2014，15：1126.

［20］Zhou Q Y，Jia J T，Huang X，et al.The large-scale investigation of gene expression in Leymus chinensis stigmas provides a valuable resource for understanding the mechanisms of poaceae self-incompatibility.BMC Genomics，2014，15：399.

［21］Li H.A statistical framework for SNP calling，mutation discovery，association mapping and population genetical parameter estimation from sequencing data.Bioinformatics，2011，27（21）：2987-2993.

［22］Ganal M T，Polley A，Graner E，et al.Large SNP arrays for genotyping in crop plants.Journal of Biosciences，2012，37：821-828.

［2］張麗君，劉龍龍，暢志堅，等.草業(yè)學(xué)報，2015，24（7）：146-154.

A study of SNP molecular marker technology for Leymus chinensis genotyping

HOU Li-Juan1**，QI Xiao2，3**，QI Dong-Mei1，CHEN Shuang-Yan1*，LIU Gong-She1*
1.Key Laboratory of Plant Resources，Institute of Botany，The Chinese Academy of Sciences，Beijing 100093，China；2.Institute of Grassland Research of Chinese Academy of Agricultural Sciences，Hohhot 010010，China；3.National Animal Husbandry Service，Beijing 100125，China

Leymus chinensis，is a widely distributed and often ecologically dominant perennial grass in China，frequently used as research material.Because we know of no previous reports of markers in L.chinensis，we performed single nucleotide polymorphism（SNP）mining of 48 genotypes，based on the transcriptome data in the laboratory and established a method using high resolution melting（HRM）curves for L.chinensis SNP genotyping.A total of 41843 SNPs were detected from the transcriptome data by bioinformatics analysis. Transitions were nearly twice as frequent as transversions.Among 112 pairs of designed SNP primers，98 SNP primer pairs（87.5%）were able to amplify bright single bands，72 primer pairs（64.29%）had a suiTableamplified band size，and 46 SNP primer pairs（41.07%）were suiTablefor HRM genotyping.The 48 germplasm genotypes of L.chinensis could be separated by 7 primer pairs and the fingerprint of each identified.The partial HRM genotyping results were verified by sequencing.The results of genotype by sequencing analysis were consistent with the results of HRM genotyping，so confirming that the HRM genotyping was reliable.HRM for SNP genotyping first established in this study can also be used in other breeding work with L.chinensis，including parent selection，gene correlation analysis，molecular marker assisted breeding，development of functional molecular markers and new variety protection.

Leymus chinensis；transcriptome；SNP marker；HRM genotyping；fingerprint

10.11686/cyxb2015466

2015-09-29；改回日期：2015-10-26

國家草品種區(qū)域試驗項目-羊草新品種“三性”測試（農(nóng)財發(fā)［2015］36號）和國家重點基礎(chǔ)研究發(fā)展計劃（973計劃）項目（2014CB138704）資助。

侯莉娟（1989-），女，山西太原人，碩士。E-mail：jennifer＿h(yuǎn)ou762@163.com。齊曉（1982-），男，河北泊頭人，碩士。E-mail：tq07mms @sina.com。**共同第一作者These authors contributed equally to this work.

Corresponding author.E-mail：sychen@ibcas.ac.cn，liugs@ibcas.ac.cn