李國(guó)忠，陳創(chuàng)夫

（新疆石河子大學(xué)動(dòng)物科技學(xué)院，新疆石河子 832000）

高分辨率溶解曲線分析結(jié)合COLD-PCR快速檢測(cè)綿羊低豐度FecB基因

李國(guó)忠，陳創(chuàng)夫*

（新疆石河子大學(xué)動(dòng)物科技學(xué)院，新疆石河子 832000）

本文利用高分辨率熔解曲線分析結(jié)合COLD-PCR建立綿羊骨形態(tài)發(fā)生蛋白受體IB（BMPR-IB）基因低豐度A746G點(diǎn)突變的檢測(cè)方法，為FecB基因大面積高通量篩查提供方便快捷和低成本的檢測(cè)技術(shù)。實(shí)驗(yàn)分別提取綿羊BMPR-IB基因A746G點(diǎn)突變陽(yáng)性純合子和陰性純合子的基因組，分別PCR擴(kuò)增223 bp的目的片段，測(cè)定濃度和純度后，分別稀釋至0.05 ng/μL，將陽(yáng)性純合子用陰性純合子進(jìn)行系列稀釋作為模板標(biāo)準(zhǔn)品。根據(jù)模板標(biāo)準(zhǔn)品序列設(shè)計(jì)合成3對(duì)引物，以標(biāo)準(zhǔn)品為模板進(jìn)行COLD-PCR，其產(chǎn)物進(jìn)行高分辨率熔解曲線分析，生成的歸一化溫度轉(zhuǎn)移差異圖（Normlized and Temp-Shifted Differines Plot）即為標(biāo)準(zhǔn)曲線圖，作為待測(cè)樣品FecB檢測(cè)的依據(jù)，對(duì)其可靠重復(fù)性做了進(jìn)一步驗(yàn)證。結(jié)果表明：2對(duì)引物COLD-PCR HRM標(biāo)準(zhǔn)曲線圖曲線分離清晰可靠，檢測(cè)下限均為0.5%，與傳統(tǒng)qPCR HRM相比，靈敏度提高到4倍。本實(shí)驗(yàn)建立的綿羊低豐度FecB基因檢測(cè)方法可以用于綿羊FecB基因大面積高通量快速篩查工作。

綿羊；FecB基因；COLD-PCR；高分辨率熔解曲線分析

綿羊產(chǎn)羔率屬于閾性狀，由多個(gè)基因控制。其中，綿羊骨形態(tài)發(fā)生蛋白受體IB（Bone Morphogenetic Protein Receptor IB, BMPR-IB）基因發(fā)生第六外顯子A746G單堿基突變后被稱作FecB基因，是影響綿羊產(chǎn)羔率的最重要的主效基因，增羔效應(yīng)最明顯，應(yīng)用最廣泛。Davis等[1]報(bào)道，一個(gè)拷貝的FecB基因能增加 1.3～1.6個(gè)卵，雙拷貝能增加 2.7～3.0個(gè)卵；對(duì)于產(chǎn)羔數(shù)而言，第1個(gè)拷貝能增加窩產(chǎn)羔數(shù)0.7～0.9只，第2個(gè)拷貝增加0.4～0.5只。

FecB基因不僅是控制我國(guó)小尾寒羊、湖羊和策勒黑羊等多胎綿羊品種的主效基因，而且有的單胎品種也可能有攜帶FecB的少數(shù)個(gè)體，孫紅霞[2]在寧夏灘羊中檢測(cè)到低頻率的FecB基因。檢測(cè)出這些攜帶有FecB基因的羊只留作種用，對(duì)提高產(chǎn)羔率有非常重要的意義。目前報(bào)道的檢測(cè)方法有AS-PCR、RFLP-PCR、測(cè)序和 HRM等[3-5]。這些方法檢測(cè)靈敏度較低，一般只用于單只羊的檢測(cè)，如果用于羊群大面積篩查，則費(fèi)時(shí)費(fèi)工，成本太高。本研究先通過(guò)低變性溫度共擴(kuò)增PCR（Coamplication at Lower Denaturation Temperature PCR，COLD-PCR）對(duì)FecB進(jìn)行有效富集，然后利用高分辨率熔解曲線分析（High Resolution DNA-melting Analysis，HRM）進(jìn)行快速檢測(cè)，建立了一種低豐度FecB檢測(cè)方法，適用于綿羊FecB基因的高通量大面積篩查工作，并快速鎖定找到攜帶有FecB的羊只，以便選種留種之參考。該法方便快捷，準(zhǔn)確高效，成本低廉。

1 材料與方法

1.1 實(shí)驗(yàn)材料 羊只血液采集于新疆石河子市屠宰場(chǎng)。用事先含有1 mL ACD 抗凝劑的無(wú)菌EP管分別收集15只宰殺湖羊和2只宰殺哈薩克羊的靜脈血液4 mL，-20℃凍存。

1.2 儀器 Light Cycler?480定量PCR儀，BioDoc -It?Imager System凝膠成像系統(tǒng)，DYY－4C電泳儀，BIO-rad C1000 PCR儀，Eppendorf 5415D超速離心機(jī)，Nanodrop 2000分光光度計(jì)。

1.3 試劑 通用型柱式基因組提取試劑盒、快速瓊脂糖凝膠DNA回收試劑盒和Goadstar DNA Polymerase購(gòu)自于北京康為世紀(jì)公司，PrimeSTAR?HS DNA Polymerase with GC Buffer購(gòu)自寶生物工程(大連)有限公司，MgCl2購(gòu)自北京天根生化科技有限公司，LC Green Plus+購(gòu)自美國(guó)Idaho公司。

1.4 基因組的提取 血液基因組的提取按照試劑盒說(shuō)明進(jìn)行，用分光光度計(jì)測(cè)定濃度和純度，滿足實(shí)驗(yàn)要求。

1.5 FecB基因檢測(cè) 首先根據(jù)FecB基因序列（GenBank登錄號(hào)：AF357007.1）設(shè)計(jì)1對(duì)引物，靶向第六外顯子，PCR產(chǎn)物長(zhǎng)度為179 bp。引物序列：F：5'- GGAAACCCTGAACATCGCTAATAC-3'；R：5'-GTCGCTATGGGGAAGTTTGGATG-3'。分別以提取的各綿羊基因組為模板，用高保真酶進(jìn)行普通PCR。PCR反應(yīng)體系 （50 μL）：2×PrimeSTAR GC Buffer（Mg2+plus）25 μL； dNTP Mixture （2.5 mmol/L each） 4 μL；上下游引物（25 umol/L）各0.4 μL；基因組模板1 μL； PrimeSTAR HS DNA Polymerase （2.5 U/μL）0.5 μL；雙蒸水補(bǔ)足50 μL。PCR反應(yīng)程序：95℃ 4 min；98℃ 10 s，54℃ 10 s，72℃ 20 s，共35個(gè)循環(huán)；72℃ 5 min。

最后按照試劑盒操作說(shuō)明，回收PCR產(chǎn)物，用分光光度計(jì)測(cè)定純度和濃度，送生物工程（上海）股份有限公司測(cè)序，選取質(zhì)量良好的FecB陽(yáng)性純合子和陰性純合子基因組，-20℃凍存?zhèn)溆谩?/p>

1.6 定量PCR標(biāo)準(zhǔn)品模板的制備 首先根據(jù)綿羊FecB基因序列，靶向其第六外顯子設(shè)計(jì)合成1對(duì)引物223-L/223-R，PCR產(chǎn)物長(zhǎng)223 bp，引物序列：223-L：5'-TCAGATGGTGAAACAGATTGGAAA-3'；223-R:5'-AGAAAGAGAGGAAAGCTAGGAAAC C-3'。然后分別以經(jīng)過(guò)測(cè)序檢測(cè)為FecB陽(yáng)性純合子和陰性純合子的備用基因組為模板，用高保真酶PCR擴(kuò)增目的序列。PCR反應(yīng)體系 ( 50 μL )：2×PrimeSTAR GC Buffer(Mg2+plus) 25 μL ；dNTP Mixture (2.5 mmol/L each) 5 μL ；上下游引物（25 umol/L）各0.5μ L ；基因組模板1 μL； PrimeSTAR HS DNA Polymerase (2.5 U/μL) 0.6μL ；雙蒸水補(bǔ)足50 μL。PCR反應(yīng)程序：95℃ 4 min；98℃ 10 s，58℃ 10 s，72℃ 20 s，共30個(gè)循環(huán)；72℃ 5 min。

PCR產(chǎn)物進(jìn)行凝膠電泳，按照快速瓊脂糖凝膠DNA回收試劑盒操作說(shuō)明，回收223 bp的目的序列，用Nanodrop 2000分光光度計(jì)測(cè)定純度和濃度，優(yōu)先選取純度最高的回收產(chǎn)物用于后續(xù)實(shí)驗(yàn)。

最后用滅菌雙蒸水將回收的223 bpPCR產(chǎn)物稀釋至0.05 ng/μL；然后按照一定的比例混合，使得FecB陽(yáng)性率分別為100%、50%、25%、20%、15%、10%、5%、2.5%、2%、1%、0.5%、0.25%、0.1%和0%，作為后續(xù)定量PCR的模板標(biāo)準(zhǔn)品。

1.7 引物設(shè)計(jì)合成 根據(jù)標(biāo)準(zhǔn)品模板的序列設(shè)計(jì)合成3對(duì)引物，使FecB基因的A746G突變單堿基位于擴(kuò)增產(chǎn)物的大概中間位置。引物HPLC純化。普通PCR檢測(cè)引物的擴(kuò)增效率和特異性。引物序列和擴(kuò)增片段長(zhǎng)度見表1。

1.8 qPCR HRM基因分型 用合成的3對(duì)引物和100%、50%和0%的模板標(biāo)準(zhǔn)品，進(jìn)行普通定量PCR結(jié)合HRM分型實(shí)驗(yàn)。此操作在Light Cycler?480上閉管完成。

1.8.1 qPCR qPCR的體系為：5×Goadstar Buffer 4 μL，超純dNTP Mixture（2.5 mmol/L each）0.2 μL，上下游引物（25 μmol/L ）各0.2 μL，標(biāo)準(zhǔn)品模板0.3 ng，Goadstar Taq Polymerase （5 U/μL）0.2 μL，10xLC Green Plus+ 1 μL，雙蒸水補(bǔ)足20 μL。qPCR 反應(yīng)程序：95℃ 10 min； 95℃ 20 s，退火（表1）25 s（Fluorescence Reading），72℃ 15 s，共35個(gè)循環(huán)；72℃ 3 min。

表1 普通定量PCR（qPCR）和COLD-PCR引物

1.8.2 HRM 分型 HRM分析程序：95℃ 30 s，40℃ 4 min，65℃ 1 s，90℃（25次Fluorescence Reading/℃），40℃ 1 min。Light Cycler?480 自帶軟件進(jìn)行掃描分析，生成Normlized and Temp-Shifted Differienc Plot圖，依據(jù)該圖分型。

1.8.3 qPCR HRM基因分型的靈敏度 100%、50%、25%、20%、10%、5%、2%、1%和0%標(biāo)準(zhǔn)品為模板，用引物Primer1和Primer3進(jìn)行qPCR HRM基因分型，以檢測(cè)qPCR HRM的靈敏度。PCR反應(yīng)體系條件和HRM分析程序同上。

1.9 COLD-PCR HRM檢測(cè)低豐度FecB基因

1.9.1 COLD-PCR條件的優(yōu)化 COLD-PCR HRM檢測(cè)效果受Tc值、Mg2+濃度、雜交時(shí)間、熒光染料、模板和引物質(zhì)量等多方面的影響[6]。本實(shí)驗(yàn)Tc值根據(jù)前人經(jīng)驗(yàn)[7]，以Tm-1為基礎(chǔ)，在其上下游設(shè)定Tc梯度，通過(guò)COLD-PCR擴(kuò)增效率和 HRM分析的效果確定最佳Tc值。Mg2+濃度優(yōu)化：以qPCR體系含Mg2+濃度分別為1.5、2.0、3.0、3.5 mmol/L，優(yōu)化最佳濃度。雜交時(shí)間設(shè)為40 s，2 、4 min 和7 min，進(jìn)行優(yōu)化；LC Green Plus+量均為0.5x。

1.9.2 COLD-PCR COLD-PCR體系：和qPCR相同。PCR反應(yīng)程序：預(yù)變性95℃ 10 min； 先運(yùn)行5個(gè)普通PCR循環(huán)：95℃ 20 s，退火溫度（表1）25 s（Fluorescence Reading），72℃ 15 s；然后35個(gè)COLD-PCR循環(huán)：95℃ 20 s，68℃ 40 s，Tc溫度（表2）20 s ，退火溫度（表1）25 s（Fluorescence Reading），72℃ 15 s；72℃ 5 min。

1.9.3 HRM 分析 HRM 分析程序與前面相同。Light Cycler?480自帶軟件進(jìn)行Gene Scaning分析，生成FecB基因檢測(cè)標(biāo)準(zhǔn)曲線圖。得到清晰的標(biāo)準(zhǔn)曲線后，重復(fù)一次實(shí)驗(yàn)，檢測(cè)其穩(wěn)定重復(fù)性。

2 結(jié)果分析

2.1 FecB檢測(cè) 經(jīng)測(cè)序，15只湖羊中有2只為FecB陽(yáng)性純合子，哈薩克羊均為陰性純合子。選取濃度高純度高的FecB陽(yáng)性純合子和陰性純合子基因組用于后續(xù)實(shí)驗(yàn)。

2.2 qPCR HRM分型結(jié)果 3對(duì)引物的分型圖清晰，效果較好（圖1）。各基因型PCR產(chǎn)物經(jīng)測(cè)序驗(yàn)證正確無(wú)誤。qPCR HRM檢測(cè)下限一般為5%～10%[8]，本實(shí)驗(yàn)用引物Primer1和Primer3做了檢測(cè)靈敏度分析，結(jié)果顯示，檢測(cè)下限均為2%（圖2）。

2.3 COLD-PCR條件的優(yōu)化 各對(duì)引物PCR產(chǎn)物的Tm值和COLD-PCR最佳Tc值見表2。

經(jīng)濃度梯度的實(shí)驗(yàn)檢測(cè)，本實(shí)驗(yàn)Mg2+的最佳濃度為1.5 mmol/L（Bufffer自帶），額外增加Mg2+使擴(kuò)增效率和特異性均下降；早期報(bào)道的COLDPCR異源雙鏈雜交時(shí)間一般為數(shù)分鐘，近年，有研究用PhusionTMHigh-fidelity Polymerase或HotStar Taq DNA Polymerase等酶及其Buffer，最佳雜交時(shí)間降為30 s[9-11]，這可能和所使用的試劑有關(guān)。本實(shí)驗(yàn)最佳雜交時(shí)間為40 s，如果延長(zhǎng)時(shí)間則PCR特異性下降，熔解峰圖凌亂，這可能與Buffer組成等有關(guān)，降低了雜交所需必要時(shí)間。

表2 2對(duì)引物PCR產(chǎn)物的Tm值和COLD-PCR的 Tc值 ℃

2.4 COLD-PCR HRM檢測(cè)結(jié)果 用Primer1和Primer2引物以及優(yōu)化的最佳條件進(jìn)行COLDPCR HRM分析，結(jié)果見圖3（a和b、c和d、），COLD-PCR HRM檢測(cè)靈敏度均為0.5%，與傳統(tǒng)定量PCR HRM相比，靈敏度提高到4倍。用Primer1引物進(jìn)行的獨(dú)立重復(fù)實(shí)驗(yàn)（圖3 e和f）表明，得到的標(biāo)準(zhǔn)曲線具有重復(fù)穩(wěn)定性。

3 討 論

本實(shí)驗(yàn)建立的qPCR HRM方法，適用于檢測(cè)單只或少數(shù)羊只混合樣品的FecB情況；建立的COLD-PCR HRM低豐度檢測(cè)方法，除此之外，還適用于大面積高通量篩查工作。本文2對(duì)引物的COLD-PCR HRM檢測(cè)方法，靈敏度均為0.5%。羊?yàn)槎扼w，所以，如果將每只羊的基因組等量混合制作混合池，以此作為模板，1個(gè)PCR反應(yīng)樣孔1次最多檢測(cè)100只羊。按此計(jì)算，如果用96孔板，一次PCR反應(yīng)即可完成數(shù)千甚至近萬(wàn)只羊FecB總頻率的檢測(cè)，非常方便快捷，而且成本低廉。如果想要在陽(yáng)性混合樣本中，找出到底哪只或哪些羊只攜帶有FecB，可以將樣本分成若干組，每組再檢測(cè)一次，縮小范圍，這樣，可以快速鎖定并找出陽(yáng)性個(gè)體，最后測(cè)序驗(yàn)證。另外，經(jīng)獨(dú)立實(shí)驗(yàn)驗(yàn)證，所建立的標(biāo)準(zhǔn)曲線圖具有重復(fù)穩(wěn)定性，這為FecB基因

選種育種帶來(lái)便利。

圖1 qPCR HRM 分型結(jié)果

本實(shí)驗(yàn)中，與傳統(tǒng)PCR相比，COLD-PCR將高分辨率熔解曲線分析的靈敏度提高到4倍。COLD-PCR對(duì)反應(yīng)孔溫的精確控制有很高要求，而本實(shí)驗(yàn)所用定量PCR儀溫度設(shè)置只能設(shè)整數(shù)值，所以最佳Tc的摸索和設(shè)定受到了限制，有可能限制了突變富集的程度。 Milbury等[10]在COLD-PCR的基礎(chǔ)上，建立了ice-COLD-PCR；Kit等[11]進(jìn)一步發(fā)展了ice-COLD-PCR，提出了E-ice-COLDPCR的方法，降低了對(duì)定量PCR儀控溫性能的要求，并提高了突變富集程度[12-13]。但這些方法下游結(jié)合HRM檢測(cè)的技術(shù)還未見報(bào)道。

圖3 COLD-PCR HRM檢測(cè)結(jié)果

[1] Davis G H, Montgomery G W, Aiaason A J, et al. Segregation of a major gene influencing fecundity in progeny of Booroola sheep[J] . New Zealand J Agric Res, 1982, 67 (4):525-529.

[2] 孫紅霞. 3個(gè)綿羊品種BMPR-IB、BMP15、GDF9和INHA基因多態(tài)性及與產(chǎn)羔數(shù)的相關(guān)性[D]. 楊凌;西北農(nóng)林科技大學(xué), 2009.

[3] 蔣立斌, 張?jiān)粕? 楊華, 等. 非探針?lè)ǜ叻直媛嗜劢馇€( HRM)對(duì) FecB 基因分型方法的研究[J]. 黑龍江畜牧獸醫(yī), 2015, ( 1):11-14.

[4] 陳曉勇, 孫洪新, 向海, 等. 等位基因特異性 PCR 技術(shù)在綿羊基因SNP分型中的建立[J]. 畜牧與獸醫(yī), 2015, 47 ( 1):49-52.

[5] 儲(chǔ)明星, 狄冉, 葉素成, 等. 綿羊多胎主效基因FecB分子檢測(cè)方法的建立與應(yīng)用[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào), 2009, 17( 1):52-58.

[6] Milbury C A, Li J G. COLD-PCR–Enhanced high-resolution melting enables rapid and selective identification of lowlevel unknown mutations[J]. Clin Chem, 2009, 55(12): 2130-2143.

[7] Milbury C A, Li J, Liu P, et al. COLD-PCR: improving the sensitivity of molecular diagnostics assays[J]. Expert Rev Mol Diagn, 2011, 11(2): 159-169.

[8] Milbury C A, Chen C C, Mamon H, et al. Multiplex amplification coupled with COLD-PCR and high resolutionmelting enables identification of low-abundance mutations in cancer samples with low DNA content[J]. J Mol Diagn, 2011, 13(2):220-232.

[9] Castellanos-Rizaldos E, Milbury C A, Guha M , et al. COLD-PCR Enriches Low-Level Variant DNA Sequences and Increases the Sensitivity of Genetic Testing[J]. Met in Mol Biol, 2014, 1102:623.

[10] Milbury C A, Li J, Makrigiorgos G M. Ice-COLD-PCR enables rapid amplification and robust enrichment for lowabundance unknown DNA mutations[J]. Nucleic Acids Res, 2011, 39(1):e2.

[11] Kit A H, Mazaleyrat N, Daunay A, et al. Sensitive detection of KRAS mutations using enhanced-ice-COLD-PCR mutation enrichment and direct sequence identification[J]. Hum Mutat, 2013, 34(11):1568-1580.

[12] How-Kit A, Tost J. Pyrosequencing?-based identification of low-frequency mutations enriched through enhanced-ice-COLD-PCR[J]. Meth Mole Biol, 2015, 1315:83-101.

[13] Schermelleh L, Haemmer A, Spada F, et al. Ultrasensitive detection and identification of BRAF V600 mutations in fresh frozen, FFPE, and plasma samples of melanoma patients by E-ice-COLD-PCR[J]. Anal Bioanal Chem, 2014, 406(22):5513-5520.

S826.2

A

10.19556/j.0258-7033.2017-06-122

2016-12-07；

2017-04-25

李國(guó)忠（1968-），男，新疆石河子人，博士研究生，主要從事動(dòng)物轉(zhuǎn)基因研究，E-mail: lgz8791@163.com

* 通訊作者：陳創(chuàng)夫（1964-），男，教授，博士生導(dǎo)師，主要從事病原微生物分子機(jī)制研究，E-mail: ccf－xb@ 163.com

圖2 qPCR HRM檢測(cè)靈敏度