李 慧， 郎坤玲， 沈長朋， 李 潔， 陶云荔， 褚 棟

青島農(nóng)業(yè)大學(xué)農(nóng)學(xué)與植物保護學(xué)院,山東省植物病蟲害綜合防控重點實驗室,山東 青島 266109

微衛(wèi)星DNA(microsatellite DNA)一般指基因組中由短的重復(fù)單元(一般為1～6個堿基)組成的DNA串聯(lián)重復(fù)序列，被稱為短串聯(lián)重復(fù)序列(short tandem repeats，STRs)或簡單重復(fù)序列(simple sequence repeats，SSRs)。微衛(wèi)星DNA廣泛分布于真核生物的基因組中，在原核生物的基因組中也有少量分布(羅文永等，2003)。微衛(wèi)星DNA符合孟德爾遺傳模式，共顯性遺傳，通常具有豐富的多態(tài)性，易于檢測，是一類很好的分子標(biāo)記(代金霞，2005；何平，1998；劉佳妮等，2008)。微衛(wèi)星DNA分子標(biāo)記有廣泛的應(yīng)用價值，現(xiàn)已被應(yīng)用于入侵害蟲的起源、入侵途徑及模式和種群擴散等研究中(褚棟等，2007、2012；高長生等，2011；Chuetal.，2011)。

微衛(wèi)星DNA的獲得方式主要有幾種途徑：(1)篩選基因文庫法，用含有微衛(wèi)星的探針在基因文庫中篩選含有微衛(wèi)星序列的陽性克隆并測序，再根據(jù)微衛(wèi)星DNA兩端的側(cè)翼序列設(shè)計引物，最后用該引物對應(yīng)的微衛(wèi)星DNA位點做定性分析(Rassmannetal.，1991)。此方法工作量大，效率低。(2)微衛(wèi)星富集法，用含微衛(wèi)星序列的探針進行雜交富集，構(gòu)建其富集文庫(Kandpaletal.，1994; Karagyozovetal.，1993)。這種方法操作較復(fù)雜，大量擴增目的片段的同時，產(chǎn)生了大量的冗余序列(張增翠和侯喜林，2004)。(3)搜索GenBank、EMBL 和 DDBJ 等公共數(shù)據(jù)庫以獲得含SSR序列的方法，開發(fā)的微衛(wèi)星位點的多態(tài)性一般不是很高(段惠生等，2012)。近年來，由于錨定PCR技術(shù)簡便、高效、多態(tài)性好，在微衛(wèi)星篩選中得到了廣泛應(yīng)用(盛良明等，2007；張俊鵬等，2012)。該方法是由Fisheretal.(1996)提出，用錨定簡并微衛(wèi)星引物對基因組DNA進行擴增。此方法是一種分離微衛(wèi)星的簡單、快捷的方法，避免了SSR富集的隨機性，具有明確的目的性和較強的富集能力，提高了微衛(wèi)星的分離效率，且通常每個序列都含有2個完整的微衛(wèi)星位點，開發(fā)效率極高(張增翠和侯喜林，2004)。

目前，桃小食心蟲CarposinasasakiiMatsumura、桃蛀螟Conogethespunctiferalis(Guenée)、玉米螟Ostrinianubilalis(Hübner)、二點委夜蛾P(guān)roxenuslepigone(Moschler)、花薊馬Frankliniellaintonsa(Trybom)、黃胸薊馬Thripshamaiiensis(Morgan)、棕櫚薊馬ThripspalmiKarny、斑翅果蠅Drosophilasuzukii(Matsurmura)、稻水象甲LissorhoptrusoryzophilusKuschel、扶桑綿粉蚧PhenacoccussolenopsisTinsley等這些重要農(nóng)業(yè)害蟲的微衛(wèi)星位點尚未見報道。本試驗用錨定PCR技術(shù)對上述10種重要農(nóng)業(yè)害蟲進行微衛(wèi)星DNA位點篩選，并對微衛(wèi)星DNA序列進行分析比較，旨在探索該方法在農(nóng)業(yè)害蟲微衛(wèi)星DNA篩選中的有效性，為進一步挖掘農(nóng)業(yè)害蟲的微衛(wèi)星位點及研究其種群遷移擴散等奠定基礎(chǔ)。

1 材料與方法

1. 1 試驗材料及DNA提取

本研究中所用的物種均放于95%乙醇中，-20 ℃下保存?；蚪MDNA的提取方法參照Chuetal.(2005)，供試蟲體樣本數(shù)量均為1頭，提取的DNA放于-20 ℃下保存。

1. 2 錨定PCR擴增

以提取的基因組DNA為PCR反應(yīng)的模板，使用根據(jù)需要設(shè)計的5′錨定簡并引物NNNNNNNNKKVRVRV(CA)6進行擴增。建立50 μL反應(yīng)體系，體系中含有DNA模板3 μL、buffer 5 μL、dNTPs 1 μL、簡并引物1 μL(10 μmol·L-1)、ddH2O 39. 5 μL、Taq酶0. 5 μL。PCR反應(yīng)程序：94 ℃預(yù)變性4 min，35個循環(huán)(94 ℃變性30 s，55 ℃復(fù)性45 s，72 ℃延伸2 min)，72 ℃再延伸7 min。

1. 3 克隆與測序

將擴增產(chǎn)物經(jīng)2%瓊脂糖凝膠電泳檢測，使用試劑盒(北京全式金生物技術(shù)有限公司)切膠回收長度500～750 bp的DNA，與PMD18-T連接并轉(zhuǎn)入大腸桿菌Trans1-T1感受態(tài)細胞，涂抹于含有氨芐西林(0. 1 mg·mL-1)的LB瓊脂平板上。挑取單菌落，放于含有氨芐西林(0. 1 mg·mL-1)的LB液體培養(yǎng)基中振蕩培養(yǎng)。取2 μL菌液，用通用引物M13對插入片段進行PCR擴增。13 μL反應(yīng)體系中含有菌液2 μL、buffer 1. 3 μL、dNTPs 0. 26 μL、M13通用引物0. 26 μL(10 μmol·L-1)、ddH2O 9. 3 μL、Taq酶0. 13 μL。PCR程序：94 ℃預(yù)變性4 min，35個循環(huán)(94 ℃變性30 s，55 ℃復(fù)性45 s，72 ℃延伸2 min)，72 ℃再延伸7 min。PCR產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測，挑選陽性克隆。

1. 4 序列測定與分析

隨機選取片段大小為500～750 bp的50個陽性克隆進行測序。利用Sequencher 5. 0 Demo軟件(http:∥www.genecodes.com/)去除冗余序列，再用SSR Hunter1. 3軟件對DNA序列進行SSR位點的搜索，設(shè)置重復(fù)數(shù)至少為5，構(gòu)成重復(fù)元件的核苷酸數(shù)最多為6個。

2 結(jié)果與分析

2. 1 簡并引物的擴增及PCR產(chǎn)物的克隆與測序

用含有CA重復(fù)單元的5′錨定簡并引物擴增10個物種的DNA，擴增片段較彌散，多集中于200～2000 bp。經(jīng)PCR篩選，各個物種的克隆片段在500～750 bp之間。陽性克隆率在10個物種中有較大的差異，最高為棕櫚薊馬(95. 8%)，最低為黃胸薊馬(20. 9%)，平均陽性克隆率為57. 8%。10個物種的微衛(wèi)星比率在70. 8%～100%之間，最高為扶桑綿粉蚧(100%)，最低的是二點委夜蛾(70. 8%),微衛(wèi)星平均比率為91. 3%。克隆效率差異較大，最高為稻水象甲(88. 2%)，最低為黃胸薊馬(19. 9%)，克隆平均效率為52. 8%(表1)。

表1 微衛(wèi)星文庫的篩選結(jié)果Table 1 The screening result of microsatellite library

*陽性克隆率指陽性克隆占所有克隆的百分比；**微衛(wèi)星比率指含有微衛(wèi)星的克隆占所有陽性克隆的百分比；***克隆效率指含有微衛(wèi)星的克隆占所有克隆的百分比。*Positive clone rate refers to the percentage of positive clones in all clones；**Microsatellite rate refers to the percentage of microsatellites identified in all positive clones；***Cloning efficiency refers to the percentage of microsatellites identified in all clones.

2. 2 微衛(wèi)星類型

去除重復(fù)序列后進行微衛(wèi)星搜索和分類統(tǒng)計，結(jié)果表明，冗余率(重復(fù)序列的數(shù)量/總序列的數(shù)量)在不同物種間差異較大，最低為棕櫚薊馬(2. 3%)，最高為扶桑綿粉蚧(71. 7%)(圖1)。各個物種平均每條序列含微衛(wèi)星位點數(shù)為1. 9～2. 4，桃小食心蟲、扶桑綿粉蚧和稻水象甲的平均位點數(shù)小于2，其他物種均大于2。

從微衛(wèi)星位點(核苷酸堿基組成)分布類型的角度分析，除簡并引物中所含有的CA/TG重復(fù)序列外，還發(fā)現(xiàn)其他類型的重復(fù)序列，如：二核苷酸AT、GC、CT/AG等3種重復(fù)單元；三核苷酸AAC、GGT、GCA、GCC、TCC等10種重復(fù)單元；四核苷酸TCGC、GACA、GTGC、ATCT、TGTC等5種重復(fù)單元。不同物種以二核苷酸CA/TG重復(fù)序列數(shù)目最多，在89. 2%～100%之間，二核苷酸(CT/AG、AT、GC)、三核苷酸和四核苷酸所占比例較低。其中，花薊馬、扶桑綿粉蚧和稻水象甲二核苷酸的比例為100%，扶桑綿粉蚧的CA/TG重復(fù)為100%(表2)。從微衛(wèi)星位點(重復(fù)次數(shù))分布類型的角度分析，10個物種中，具有最高重復(fù)次數(shù)的是玉米螟，達39次，其最低重復(fù)次數(shù)為5(SSR hunter的最低重復(fù)次數(shù)設(shè)置為5次)。大部分物種的SSR序列較短，如桃小食心蟲的SSR片段長度為10～36 bp，玉米螟的SSR片段長度為10～78 bp。各個物種的平均重復(fù)次數(shù)為6. 7～8. 9次(圖2)。

根據(jù)Weber(1990)的微衛(wèi)星分類標(biāo)準(zhǔn)，本研究中涉及的10個物種的微衛(wèi)星類型90%以上都是完全型，不完全型和復(fù)合型的SSR所占比例較少。不完全型的微衛(wèi)星比率最高的物種為玉米螟(5. 7%)，復(fù)合型的微衛(wèi)星比率最高的物種為桃小食心蟲(5. 0%)。其中，扶桑綿粉蚧和桃蛀螟的微衛(wèi)星類型均為完全型(圖3)。

圖1 序列冗余率分布

表2 微衛(wèi)星位點類型分布(基于核苷酸堿基數(shù)目)Table 2 Distribution of microsatellite types (based on nucleotide numbers)

*包括重復(fù)單位CTG、GCA、GCC、TCC、GAA、AGA、TTA、TCC、AAC、CCA;**包括重復(fù)單位TGTC、GTGC、ATCT、TCTG、GACA。*including the repeat unit as follows: CTG, GCA, GCC, TCC, GAA, AGA, TTA, TCC, AAC, CCA;**including the repeat unit as follows: TGTC, GTGC, ATCT, TCTG, GACA.

3 討論

本研究表明，10個物種中陽性克隆率最高為95. 8%，平均陽性克隆率為57. 8%，高于通過探針雜交法獲得的陽性克隆率(安建東等，2011；郝卓然等，2012;吉亞杰和張德興，2004)。本研究中，微衛(wèi)星比率最高為100%，微衛(wèi)星平均比率為91. 3%，高于通過探針雜交法獲得的微衛(wèi)星比率(程月琴等，2013;張國彥和翟保平，2008)。盛良明等(2007)利用錨定PCR技術(shù)測定結(jié)果表明，陽性克隆率和微衛(wèi)星比率均達到了100%；張俊鵬等(2012)也證明利用該技術(shù)篩選的微衛(wèi)星比率極高，達97. 7%。徐波等(2012)對使用磁珠富集法與錨定PCR法開發(fā)背瘤麗蚌Lamprotulaleai(Gray)微衛(wèi)星標(biāo)記進行比較，發(fā)現(xiàn)錨定PCR技術(shù)獲得的陽性克隆比率高。此外，本試驗也發(fā)現(xiàn)10個物種的冗余率(2. 3%～71. 7%)、假陽性率(0～29. 2%)之間差異較大。

圖2 微衛(wèi)星重復(fù)次數(shù)分布圖Fig.2 Distribution of repeat number of microsatellite in the species analysed CN=Carposina nipponensis, CP=Conogethes punctiferalis, ON=Ostrinianubilalis, PL=Proxenus lepigone, FI=Frankliniella intonsa, TH=Thrips hawaiiensis, TP=Thrips palmi, DS=Drosophila suzukii, LO=Lissorhoptrus oryzophilus, PS=Phenacoccus solenopsis.

圖3 微衛(wèi)星類型分布[基于Weber(1990)標(biāo)準(zhǔn)]

在所有動物基因組中，(CA)n微衛(wèi)星的含量最為豐富(Brenneretal.,1993)。因此，本研究使用(CA)n簡并引物來研究10個物種，以獲得更多的微衛(wèi)星位點。本試驗中，二核苷酸、三核苷酸、四核苷酸重復(fù)單元在鱗翅目、纓翅目、雙翅目、同翅目和鞘翅目中差異不明顯，二核苷酸在不同目間的比例為95. 1%～100%，三核苷酸的比例為0～4. 9%，四核甘酸為0～1. 73%，較高的二核苷酸比例可能與錨定引物篩選方法有關(guān)。

本試驗所設(shè)計的錨定簡并引物中含有6個微衛(wèi)星重復(fù)，測序結(jié)果表明，獲得的克隆兩端大部分至少含有6個二核苷酸重復(fù)，并且每條序列大部分都含有2個以上的微衛(wèi)星位點，說明錨定簡并引物擴增可以有效獲得微衛(wèi)星DNA位點。本試驗中分離的微衛(wèi)星的平均重復(fù)數(shù)為6. 35～7. 50次，最高重復(fù)次數(shù)為34次。馬麗華(2008)通過該方法獲得的微衛(wèi)星的平均重復(fù)數(shù)為7. 19次。而通過磁珠富集法獲得的微衛(wèi)星的重復(fù)次數(shù)大部分都在10次以上(安建東等，2011；連灝等，2012)。這些結(jié)果表明,錨定PCR技術(shù)獲得的微衛(wèi)星重復(fù)數(shù)低于磁珠富集法(馬麗華，2008；徐波等，2012)。

Weber(1990)按照微衛(wèi)星核心序列的不同將微衛(wèi)星分為3種類型：完全型(perfect)、不完全型(imperfect)和復(fù)合型(compound)。本試驗中,10個物種的微衛(wèi)星類型90%以上都是完全型，不完全型和復(fù)合型的微衛(wèi)星DNA所占比例較小。完全型微衛(wèi)星DNA是動物基因組中比例最大的結(jié)構(gòu)類型，而不完全型和復(fù)合型微衛(wèi)星比例的高低在不同物種中存在差異(連灝等，2012；馬雅軍等，2008；王蕾等，2009；張麗等，2009)。

分子遺傳標(biāo)記是近來現(xiàn)代遺傳標(biāo)記學(xué)發(fā)展較快的領(lǐng)域之一，微衛(wèi)星則被認(rèn)為是各類分子遺傳標(biāo)記中最有價值的一種。隨著生物統(tǒng)計學(xué)研究的不斷深入，微衛(wèi)星在度量品種遺傳多樣性、估測品種間遺傳距離及構(gòu)建系統(tǒng)發(fā)生樹等研究中顯示出巨大優(yōu)勢，其應(yīng)用前景非常廣闊(徐興莉和楊虎，2011)。本研究采用5′錨定PCR技術(shù)大大降低了研究成本,提高了篩選效率，是從物種基因組中分離微衛(wèi)星位點的一個較好策略。

安建東, 黃家興, 董捷, 周冰峰. 2011. 火紅熊蜂微衛(wèi)星標(biāo)記的篩選及種特異性分析. 昆蟲學(xué)報, 54(12): 1423-1432.

程月琴, 焦振彬, 張佩, 葉永忠, 王紅衛(wèi). 2013. 地黃微衛(wèi)星富集文庫構(gòu)建及特性分析. 種子, 32(5): 12-16.

褚棟, 李顯春, 張友軍. 2012. 基于微衛(wèi)星標(biāo)記對中國Q型煙粉虱早期入侵種群與B型煙粉虱種群的遺傳結(jié)構(gòu)分析. 昆蟲學(xué)報, 55(12): 1376-1385.

褚棟, 張友軍, 萬方浩. 2007. 分子標(biāo)記技術(shù)在入侵生態(tài)學(xué)研究中的應(yīng)用. 應(yīng)用生態(tài)學(xué)報, 18(6): 1383-1387.

代金霞. 2005. 微衛(wèi)星DNA標(biāo)記技術(shù)及其應(yīng)用. 農(nóng)業(yè)科學(xué)研究, 26(1): 67-70, 79.

段惠生, 張安盛, 趙傳志, 于毅, 褚棟. 2012. 西花薊馬EST-SSR信息分析標(biāo)記篩選及其與Genomic-SSR的多態(tài)性比較. 昆蟲學(xué)報, 55(6): 634-640.

高長生, 國棟, 劉國霞, 陶云荔, 張友軍, 褚棟. 2011. Q型煙粉虱東地中海種群遺傳多樣性的mtCOI與SSR分析. 昆蟲學(xué)報, 54(12): 1416-1422.

郝卓然, 梁利群, 常玉梅, 任波. 2012. 扁吻魚微衛(wèi)星分子標(biāo)記的篩選及特征分析. 華北農(nóng)學(xué)報, 27(增刊): 40-45.

何平. 1998. 真核生物中的微衛(wèi)星及其應(yīng)用遺傳. 遺傳, 20(4): 42-47.

吉亞杰, 張德興. 2004. 鱗翅目昆蟲基因組中微衛(wèi)星DNA的特征以及對其分離的影響. 動物學(xué)報, 50(4): 608-614.

連灝, 石米娟, 杜富寬, 江遙, 黃容, 廖蘭杰, 汪亞平, 朱作言. 2012. 草魚(AG)微衛(wèi)星標(biāo)記克隆及特征分析. 水生生物學(xué), 36(1): 29-34.

劉佳妮, 桂富榮, 李正躍. 2008. SSR分子標(biāo)記技術(shù)在入侵昆蟲學(xué)研究中的運用. 植物保護, 34(3): 7-11.

羅文永, 胡駿, 李曉方. 2003. 微衛(wèi)星序列及其應(yīng)用. 遺傳, 25(5): 615-619.

馬麗華. 2008. 應(yīng)用5′錨定PCR開發(fā)谷子微衛(wèi)星標(biāo)記. 石家莊: 河北師范大學(xué).

馬雅軍, 樊勇, 吳靜. 2008. 雷氏按蚊多態(tài)微衛(wèi)星DNA位點的篩選和特征. 寄生蟲與醫(yī)學(xué)昆蟲學(xué)報, 15(3): 150-153.

盛良明, 薛華柏, 王化坤, 徐春明, 王三紅, 章鎮(zhèn). 2007. 利用5′錨定PCR技術(shù)分離枇杷微衛(wèi)星標(biāo)記. 江蘇農(nóng)業(yè)學(xué)報, 23(1): 50-53.

王蕾, 劉繼紅, 張立冬, 孫效文. 2009. 牙鲆基因組(CAG)n微衛(wèi)星DNA特征分析. 中國水產(chǎn)科學(xué), 16(6): 807-815.

徐波, 汪桂玲, 李家樂. 2012. 磁珠富集法與5′錨定PCR法開發(fā)背瘤麗蚌微衛(wèi)星標(biāo)記的比較. 生態(tài)學(xué)雜志, 31(4): 923-930.

徐興莉, 楊虎. 2011. 微衛(wèi)星DNA標(biāo)記技術(shù)的特點及其在動物研究中的應(yīng)用. 畜禽業(yè), (12): 34-35.

張國彥, 翟保平. 2008. 東方粘蟲(Pseudaletiaseparate(Walker))微衛(wèi)星富集文庫的構(gòu)建與分析. 生態(tài)學(xué)報, 28(8): 3860-3867.

張俊鵬, 孫典巧, 王日昕, 徐田軍. 2012. 刺參CT/AG微衛(wèi)星標(biāo)記的篩選. 浙江海洋學(xué)院學(xué)報： 自然科學(xué)版, 31(2): 97-102.

張麗, 樊勇, 馬雅軍. 2009. 中華白蛉微衛(wèi)星DNA序列的分離和多態(tài)位點篩選的初步研究. 中國寄生蟲學(xué)與寄生蟲病雜志, 27(6): 503-507.

張增翠, 侯喜林. 2004. SSR分子標(biāo)記開發(fā)策略及評價. 遺傳, 26(5): 763-768.

Brenner S, Elgar G, Sandford R, Macrae A, Venkatesh B A and Paricio S. 1993. Characterization of the pufferfish (Fugu) genome as a compact model vertebrate genome.Nature, 366: 265-268.

Chu D, Gao C S, De Barro P, Wan F H and Zhang Y J. 2011. Investigation of the genetic diversity of an invasive whitefly (Bemisiatabaci) in China using both mitochondrial and nuclear DNA markers.BulletinofEntomologicalResearch, 101: 467-475.

Chu D, Zhang Y J, Cong B, Xu B Y, Wu Q J and Zhu G R. 2005. Sequences analysis of mtDNA COI gene and molecular phylogeny of different geographical populations ofBemisiatabaci(Gennadius).AgriculturalSciencesinChina, 4: 533-541.

Fisher P J, Gardner R C and Richardson T E. 1996. Single locus microsatellite isolated using 5′-anchored PCR.NucleicAcidsResearch, 24: 4369-4371.

Kandpal R P, Kandpal G and Weissman S M. 1994. Construction of libraries enriched for sequence repeats and jumping clones, and hybridization selection for region-specific markers.ProceedingsoftheNationalAcademyoftheSciencesoftheUnitedStatesofAmerica, 91(1): 88-92.

Karagyozov L, Kalcheva I D and Chapman V M. 1993. Construction of random small-insert genomic libraries highly enriched for simple sequence repeats.NucleicAcidsResearch, 21: 3911-3912.

Rassmann K， Schlotterer C and Tautz D. 1991. Isolation of simple sequence loci for use in polymerase chain reaction-based DNA fingerprinting.Electrophoresis, 12(2-3): 113-118.

Weber J L. 1990. Informativeness of human (dC-dA)n. (dG-dT)npolymorphisms.Genomics, 7: 524-554.