基于磁珠富集法開發的15個三角魴多態性微衛星標記的特性分析

劉凱，謝楠，馮曉宇，馬恒甲

(杭州市農業科學研究院 水產研究所，浙江 杭州 310024)

三角魴(Megalobramaterminalis)，隸屬鯉形目(Cypriniformes)、鯉科(Cyprinidae)、鲌亞科(Culterinae)、魴屬(Megalobrama)，分布于中國嶺南以北各大水系[1]。由于過度捕撈及生態環境變化，三角魴野生資源量嚴重下降。錢塘江流域由于地方對三角魴種質資源保護、增殖放流等工作的開展，擁有一定的資源量，且建有全國唯一的國家級三角魴原種場[2]。

有關三角魴的遺傳數據較少，因此開展遺傳研究對保護三角魴至關重要。微衛星是高多態性的共顯性分子標記，是魚類遺傳多樣性研究的理想分子標記[3-5]。有學者將從團頭魴、鳊魚等開發的微衛星引物跨種用于三角魴進行遺傳多樣性研究[5-7]，但未見有關三角魴微衛星分子標記開發的報道。

本研究利用探針(CA)10，在三角魴基因組PCR文庫中，通過磁珠富集[8]的方法開發了15個具有多態性的微衛星分子標記，該項工作的開展將有利于三角魴遺傳多樣性研究，為三角魴的品種選育打下基礎。

1 材料與方法

1.1 三角魴基因組PCR文庫構建

三角魴取自杭州市農業科學研究院水產研究所(國家級錢塘江三角魴原種場)，采集胸鰭后用無水乙醇保存備用。基因組DNA采用苯酚-氯仿提取法提取[9]，提取后的基因組DNA用限制性內切酶RsaⅠ 37 ℃消化2 h，利用T4-DNA連接酶將接頭21-mer和25-mer[10]與酶切后的DNA片段在25 ℃條件下連接4 h，連接所用接頭由上海桑尼生物科技有限公司合成。用連有接頭的DNA片段作為模板，接頭21-mer作為引物，進行PCR擴增，創建基因組PCR文庫。反應程序為：94 ℃預變性3 min；94 ℃ 30 s，56 ℃ 30 s，72 ℃ 1 min，30個循環；最后72 ℃延伸5 min。

1.2 篩選富集三角魴微衛星文庫

將三角魴基因組PCR文庫10 μL、10 μmol·L-1生物素標記的(CA)10探針1 μL、20×SSC 19.5 μL和雙蒸水34.5 μL混合均勻。98 ℃下變性5 min后，55 ℃下雜交20 min。將雜交完畢的雜交液加入已平衡好的磁珠(每100 μL磁珠，分別用1 mL 1×TE和1 mL 6×SSC洗滌2次，每次5 min，最后用35 μL 6×SSC洗滌后保存備用)，25 ℃溫育30 min，并輕輕搖動，使生物素和鏈霉親和素結合。溫育結束后，將離心管放置到磁力架上，去除溶液。依次用洗液Ⅰ(2×SSC，0.1%SDS)、洗液Ⅱ(1×SSC)洗滌磁珠，去除不含微衛星的序列，30 μL雙蒸水洗脫后室溫靜置10 min，釋放出含有微衛星序列的單鏈DNA。用所獲得的單鏈DNA作為模板，以21-mer作為引物進行PCR擴增。PCR程序為：94 ℃預變性3 min；94 ℃ 30 s，56 ℃ 30 s，72 ℃ 1 min，30個循環；最后72 ℃延伸5 min。將PCR產物同T-載體(pGEM-T vector，購于Promega公司)連接，用大腸埃希菌DH5a感受態細胞進行轉化，得到三角魴微衛星文庫，經過藍白斑篩選，挑取陽性克隆進行測序分析。

1.3 三角魴微衛星的分析

以陽性克隆測序獲得的序列作為參考序列，使用SSR篩選軟件MISA進行SSR篩選。篩選標準為1個堿基重復10次及10次以上，2個堿基重復6次及6次以上，3～6個堿基重復5次及5次以上，2個微衛星之間的距離小于100 bp的時候，2個微衛星組成1個復合微衛星。用SSR出現頻率和SSR平均分布距離來描述SSR，出現頻率前2位的重復單元定義為優勢重復單元。計算公式分別為：SSR出現頻率=搜索到的SSR數量/EST序列數量；SSR平均分布頻率=EST序列總堿基數/搜索到的SSR數量。

1.4 三角魴微衛星的篩選

獲得的陽性克隆在上海桑尼生物科技有限公司進行測序。利用Primer 3 interface modules對獲得的陽性克隆序列進行預處理后，利用Primer 3進行SSR引物的批量設計，引物設計的參數是Tm為60 ℃，引物長度為20 bp，所有成功設計的引物由上海桑尼生物科技有限公司合成。將三角魴的基因組DNA樣品作為模板，對成功設計的引物進行TouchDown PCR擴增來初步篩選能擴增的引物。采用Hotstart PCR Master Mix試劑盒(購自上海近岸科技有限公司)進行TouchDown PCR擴增。擴增條件如下：94 ℃預變性3 min；94 ℃ 30 s，退火30 s，72 ℃ 20 s，起始退火溫度65 ℃，每個循環降低1 ℃，共10個循環；94 ℃30 s，55 ℃退火30 s，72 ℃20 s，共20個循環；72 ℃延伸5 min，4 ℃保存。

1.5 三角魴微衛星的分型

微衛星引物為已初步篩選能擴增的引物，由上海桑尼生物科技有限公司合成，上游引物5’端加上FAM熒光標記(表1)。利用AmpliTaq Gold Fast PCR Master Mix試劑盒(購自Applied Biosystems公司)進行TouchDown PCR擴增。擴增條件如下：94 ℃預變性3 min；94 ℃ 30 s，退火30 s，72 ℃ 20 s，起始退火溫度65 ℃，每個循環降低11 ℃，共10個循環；94 ℃ 30 s，55 ℃退火30 s，72 ℃ 20 s，共20個循環；72 ℃延伸30 min，4 ℃保溫，-20 ℃保存備用。PCR擴增產物送上海桑尼生物科技有限公司，通過ABI Prism3730 xl測序儀進行毛細管電泳。以LIZ-500為分子量內標，通過GeneMaker ver. 2.2軟件讀取微衛星擴增產物的分子量數據。

1.6 數據處理

利用Cervus ver. 3.0.7[11]分析三角魴微衛星的等位基因數、觀測雜合度、期望雜合度，多態信息含量及進行Hardy-Weinberg平衡檢驗。其中，Hardy-Weinberg平衡檢驗采用Yates修正的方差檢驗[12]，同時采用Bonferroni修正進行顯著性判別[13]。利用Micro-Checker ver. 2.2.3[14]進行空等位基因的判別并計算空等位基因頻率。

2 結果與分析

本研究中，隨機篩選陽性克隆600個，測序得到陽性克隆367個，其中255個克隆含有微衛星序列。255個微衛星序列中，含有微衛星座位355個，其中以AC/GT為單位的重復序列327個，占微衛星座位的92.1%，重復次數在20次以上的微衛星146個，占總數的41.1%，10～20次的202個，占總數的56.9%。

在255個微衛星序列中，除了一些微衛星序列因本身結構或兩端太短不能設計引物外，其他微衛星利用引物設計軟件Primer Premier 3.0設計引物15對，均為有效引物，具有多態性(表1)。

15個具有多態性的微衛星分子標記具體信息見表2。等位基因數量范圍為3～16，觀察雜合度為0.000 0～0.818 2，期望雜合度為0.081 8～0.922 6，多態信息含量為0.078 9～0.899 9。其中，具有高度多態性(PIC>0.5)的有6個標記，4個標記具有中度多態性(0.25

表2 三角魴微衛星多態性位點特征分析

注：Hardy-Weinberg平衡概率值小于經Bonferroni修正的P值(P<0.003 6)則判定偏離平衡，計算時最小期望頻率值取2，ND表示頻率值偏小而未進行Hardy-Weinberg檢驗。

3 討論

本研究中有Mt257、Mt284、Mt301、Mt309、Mt355微衛星位點未能進行Hardy-Weinberg平衡檢驗，這可能是由于該位點雜合度過低。Mt153位點未能進行Hardy-Weinberg平衡檢驗的原因可能是由于空等位基因頻率過高。此外，3個微衛星位點偏離了Hardy-Weinberg平衡，這種情況的發生可能是由于研究中樣本數量過少，該位點上存在基因流動，即Wahlund效應[16]。

利用Micro-Checker 軟件分析表明，樣品中存在空等位基因位點。低頻率空等位基因常常在許多動物的微衛星研究中發現[17-18]，頻率低于0.2是合理范圍[19]。本研究中，有Mt19、Mt153、Mt246、Mt317、Mt355微衛星位點空等位基因頻率高于0.2，可能是由于stutter峰存在或者等位基因丟失而造成的。

本研究報道了利用磁珠富集法開發的15個具有多態性的微衛星分子標記，這些新分離的分子標記將有助于三角魴種群遺傳變異、遺傳結構、種群資源保護及其他方面的研究。