南彥斌, 許嘉誠, 何啟玥, 潘學(xué)能, 周淵濤

(青海大學(xué)農(nóng)牧學(xué)院, 青海 西寧 810016)

微衛(wèi)星(Microsatellite)又名簡單重復(fù)序列(Single sequence reperts,SSR),它由1～6個(gè)核苷酸堿基串聯(lián)重復(fù)而形成,在真核生物基因組中分布比較廣泛[1]。微衛(wèi)星具有幾個(gè)非常顯著的特征,如可重復(fù)性高,共顯性,多態(tài)性高,數(shù)量豐富,對(duì)基因組有較好的覆蓋性等[2]。因此,其被廣泛應(yīng)用于群體遺傳結(jié)構(gòu)[3]、生物遺傳多樣性[4]、系統(tǒng)發(fā)生和遺傳圖譜繪制[5]等研究領(lǐng)域。雖然微衛(wèi)星標(biāo)記優(yōu)點(diǎn)突出但是無法掩蓋傳統(tǒng)的SSR引物開發(fā)方法成本極高且成功率低的缺點(diǎn)[6]。目前,高通量測(cè)序技術(shù)和生物信息學(xué)的的快速發(fā)展,使得SSR標(biāo)記開發(fā)更加高效便捷[7]。基于轉(zhuǎn)錄組數(shù)據(jù)已成功挖掘出多種昆蟲的SSR位點(diǎn)。譬如李敏等[8]從中華大仰蝽(Notonectachinensis) 37 801條unigenes中搜索到3124個(gè)SSR位點(diǎn),利用半定量實(shí)驗(yàn)篩選出16個(gè)SSR位點(diǎn),他們認(rèn)為SSR新標(biāo)記的開發(fā)在中華大仰蝽的遺傳多樣性分析基因圖譜構(gòu)建中具有重要作用。王定鋒等[9]基于轉(zhuǎn)錄組數(shù)據(jù)挖掘灰茶尺蠖(Ectropisgrisescens)微衛(wèi)星標(biāo)記,發(fā)現(xiàn)灰茶尺蠖轉(zhuǎn)錄組中SSR位點(diǎn)數(shù)量大、出現(xiàn)頻率高、基元類型十分豐富,這些信息將為今后更好地利用SSR標(biāo)記研究灰茶尺蠖種群間的遺傳分化和制定灰茶尺蠖綜合防治措施在分子生物學(xué)方面提供借鑒。韓小紅[10]等通過對(duì)星天牛(Anoplophorachinensis)轉(zhuǎn)錄組SSR位點(diǎn)特征分析,發(fā)現(xiàn)其轉(zhuǎn)錄組中SSR的主要重復(fù)類型為單堿基重復(fù),其次是三堿基重復(fù)。SSR位點(diǎn)信息和位點(diǎn)特征是昆蟲分子標(biāo)記研究的重要內(nèi)容。然而截至目前,對(duì)于青海草原毛蟲微衛(wèi)星位點(diǎn)的研究還處于空白狀態(tài)。

草原毛蟲(Gynaephora)俗稱草原毒蛾,隸屬于鱗翅目毒蛾科。在世界范圍內(nèi),目前有記載的草原毛蟲屬昆蟲共有15種,而我國獨(dú)占8種且大部分分布在青藏高原上,對(duì)青藏高原高寒甸草產(chǎn)生了嚴(yán)重的危害[11]。該蟲通常一年發(fā)生一代,其1齡幼蟲在土塊、牛糞和草根下越冬,與國外的草原毛蟲生活史周期相比國內(nèi)的生活史周期較短[12-13]。由于青海草原毛蟲對(duì)生境的要求比較特殊,故而在國內(nèi)對(duì)它的研究較為有限,國外的相關(guān)研究也僅限于毒蛾科的其他蟲種[14-15]。近十年從分子生物學(xué)的角度出發(fā)關(guān)于青海草原毛蟲的研究,主要是從腸道菌、線粒體基因、序列變化和轉(zhuǎn)錄組基因表達(dá)等方面解析青海草原毛蟲適應(yīng)高海拔環(huán)境機(jī)制[15-18]。近年來有研究者從生物防治的角度開始探究青海草原毛蟲的防治。如楊忠岐等研究金小蜂對(duì)青海草原毛蟲的寄生作用[19],陳青紅等研究4種生物藥劑對(duì)青海草原毛蟲的防控效果[20],白海濤等利用昆蟲病原線蟲防治青海草原毛蟲21等。目前,隨著國家對(duì)生態(tài)環(huán)境的重視程度的增加,未來借助于分子生物學(xué)研究,推動(dòng)青海草原毛蟲的生物防治是大勢(shì)所趨。本研究通過高通量測(cè)序技術(shù)獲得青海草原毛蟲的轉(zhuǎn)錄組數(shù)據(jù),結(jié)合生物信息學(xué)方法,對(duì)青海草原毛蟲的SSR位點(diǎn)進(jìn)行挖掘并進(jìn)行序列特征的統(tǒng)計(jì)與分析,設(shè)計(jì)和篩選能夠穩(wěn)定擴(kuò)增的SSR引物,為青海草原毛蟲的分子鑒定與遺傳多樣性分析提供參考。

1 材料與方法

1.1 供試?yán)ハx

本試驗(yàn)中的青海草原毛蟲樣品于2020年7月采自青海省海北州(36°59′ N,100°52′ E),設(shè)3個(gè)生物學(xué)重復(fù),每組3頭4齡幼蟲及雌雄成蟲,使用液氮快速冷凍后于—80℃超低溫冰箱中保存。取5只單頭4齡幼蟲用于SSR引物驗(yàn)證,使用TIANamp Genomic DNA Kit基因組DNA試劑盒提取總DNA,經(jīng)過2%的瓊脂糖凝膠電泳和微量核酸分光光度計(jì)(北京凱奧科技發(fā)展有限公司)檢測(cè),結(jié)果滿足進(jìn)一步實(shí)驗(yàn)要求,置于-20℃保存?zhèn)溆谩?/p>

1.2 總RNA提取、建庫及數(shù)據(jù)評(píng)估

委托北京諾禾致源科技股份有限公司通過Trizol法提取青海草原毛蟲樣品RNA,經(jīng)純度和濃度檢測(cè)合格后等物質(zhì)的量混合,采用Illumina HiSeqTM2000高通量測(cè)序平臺(tái)對(duì)青海草原毛蟲雌雄成蟲和4齡幼蟲進(jìn)行全長轉(zhuǎn)錄組測(cè)序、文庫構(gòu)建和數(shù)據(jù)評(píng)估。在得到原始數(shù)據(jù)后通過SMRTlink軟件獲得未組裝的Subreads,進(jìn)行校正,之后產(chǎn)生環(huán)形一致性序列(Circular-consensus sequence,CCS)。根據(jù)3′和5′引物及PloyA特征在CCS是否存在,從而確定全長序列和非全長序列,并分類找出全長非嵌合序列,通過同種類型聚類法(Iterative isoform-clustering,ICE)可將相似的全長非嵌合序列(Full length non-chimeric,FLNC)聚成一簇,得到Cluster,每個(gè)Cluster獲得一條一致性序列,結(jié)合非全長序列對(duì)其進(jìn)行校對(duì)(Polishing),篩選得到高質(zhì)量序列用于后續(xù)分析[22]。

1.3 青海草原毛蟲轉(zhuǎn)錄組unigenes的功能注釋

在Trinity v2.0.2拼接的基礎(chǔ)上,使用Corset v4.6軟件進(jìn)行轉(zhuǎn)錄本聚類,去除冗余轉(zhuǎn)錄本,得到非冗余unigenes。對(duì)聚類得到的unigenes使用blast2 go v2.5軟件同GO,KOG和KEGG數(shù)據(jù)庫比對(duì),進(jìn)行基因功能注釋。

1.4 SSR的篩選

利用MISA軟件(Micro satellite identification tool,www.pgrc.ipk-gatresleben.de/misa)對(duì)青海草原毛蟲轉(zhuǎn)錄組中1 Kb以上的unigene進(jìn)行SSR位點(diǎn)分析。對(duì)應(yīng)的各個(gè)unitsize的最少重復(fù)次數(shù)分別為10,6,5,5,5和5,重復(fù)單元的長度大小分別是1,2,3,4,5和6 bp。兩個(gè)相鄰的微衛(wèi)星重復(fù)單元間的間隔長度至少為100 bp。

1.5 SSR引物設(shè)計(jì)與驗(yàn)證

在1.4節(jié)篩選到的SSR位點(diǎn)的基礎(chǔ)上,運(yùn)用Primer3軟件[23]對(duì)引物進(jìn)行批量設(shè)計(jì)。目標(biāo)擴(kuò)增片段必須設(shè)置為包含SSR起始-3 bp,終止為+6 bp,擴(kuò)增片段的長度為700～300 bp。引物的長度設(shè)為18～25 bp,最適長度為22 bp,引物最大能允許的不能識(shí)別的堿基數(shù)為1個(gè),引物的最適退火溫度為60℃,退火溫度設(shè)置為55～65℃,上下游引物間的退火溫度差異不能超過3℃,引物末端穩(wěn)定性不能超過250[24]。

為了驗(yàn)證引物的穩(wěn)定性,從批量設(shè)計(jì)的引物中隨機(jī)選出25對(duì)引物,交由生工生物工程(上海)股份有限公司。PCR反應(yīng)體系為25 μL:青海草原毛蟲幼蟲DNA模板1 μL,上下游引物(10 μmol·L-1)各加1 μL,2×EasyTaq?PCR SuperMix 12.5 μL,無核酸酶水9.5 μL。反應(yīng)程序:94℃預(yù)變性5 min;94℃變性30 s,57℃退火30 s,72℃延伸40 s,總共35個(gè)循環(huán);72℃延伸7 min,保存在4℃下。擴(kuò)增后的產(chǎn)物用2%瓊脂糖凝膠電泳檢測(cè)(100V,30 min)。

2 結(jié)果與分析

2.1 青海草原毛蟲轉(zhuǎn)錄組測(cè)序分析

采用Illumina HiSeqTM2000高通量測(cè)序平臺(tái)對(duì)青海草原毛蟲雌雄成蟲、4齡幼蟲轉(zhuǎn)錄組數(shù)據(jù)測(cè)序,獲得71 038 419條原始reads;測(cè)序獲得unigenes序列總長度為69 070 017 bp,堿基錯(cuò)誤率為0.02%;Q20(質(zhì)量值≥20堿基所占百分比)含量平均約為98.18%,Q30(質(zhì)量值≥30堿基所占百分比)含量平均約為94.60%,測(cè)序質(zhì)量較好。采用Trinity軟件對(duì)clean reads數(shù)據(jù)進(jìn)行拼接組裝,組裝后獲得63 335條all unigenes。其中,N90(拼接轉(zhuǎn)錄本按照長度從長到短排序,累加轉(zhuǎn)錄本的長度,到不小于總長90%的拼接轉(zhuǎn)錄本的長度)和N50(拼接轉(zhuǎn)錄本按照長度從長到短排序,累加轉(zhuǎn)錄本的長度,到不小于總長50%的拼接轉(zhuǎn)錄本的長度)的長度分別為440 bp和1 714 bp,平均長度為1 091 bp。所有63 335條unigenes序列用于后續(xù)的SSR搜索。

2.2 青海草原毛蟲轉(zhuǎn)錄組的GO注釋

利用Blast2GO v2.5進(jìn)行GO注釋,將注釋成功的基因按照GO三個(gè)大類(生物學(xué)過程,細(xì)胞組成,分子功能)的下一層級(jí)進(jìn)行分類。結(jié)果表明(圖1),18 588條unigenes共得到78 461條功能注釋,分布在43個(gè)功能區(qū)。其中,獲得注釋序列數(shù)最多的是生物學(xué)過程分類,有39474條,分布在26個(gè)功能區(qū);其次是分子功能分類的注釋序列,有21 443條,分布在12個(gè)功能區(qū);注釋序列最少的是細(xì)胞組分類別,為17 544條,分布在5個(gè)功能區(qū)。在生物學(xué)過程類別中,細(xì)胞進(jìn)程和代謝進(jìn)程的注釋序列數(shù)占主導(dǎo)地位,分別為10 846和9 355條。在分子功能分類中,注釋到結(jié)合和催化活性的序列數(shù)最多,分別為9 932和7 295條。在細(xì)胞組分分類中,注釋到細(xì)胞結(jié)構(gòu)體和胞內(nèi)部分的序列占主導(dǎo)地位,分別為8 267和4 877條。有11個(gè)功能區(qū)的注釋序列超過2 000條,分別是細(xì)胞進(jìn)程(10 846),結(jié)合(9 932),代謝過程(9 355),細(xì)胞結(jié)構(gòu)體(8 267),催化活性(7 295),生物調(diào)節(jié)(4 085),生物調(diào)控進(jìn)程(3 803),含蛋白復(fù)合物(3 490),定位(3 002)和應(yīng)激響應(yīng)(2 607)。

圖1 青海草原毛蟲unigenes GO功能注釋

2.3 青海草原毛蟲Unigenes的KEGG分類

KEGG代謝通路富集分析結(jié)果表明(圖2),共有11 846條unigenes參與到細(xì)胞進(jìn)程、環(huán)境信息處理、遺傳信息處理、新陳代謝和有機(jī)體系統(tǒng)這五大類生化代謝通路中。其中代謝系統(tǒng)途徑中基因最多(3 346條),它們主要參與碳水化合物代謝和脂質(zhì)代謝等過程,分別為566條和436條。在這34組代謝通路子類別中,信號(hào)轉(zhuǎn)導(dǎo)(Signal transduction)通路的基因最多,為1 302條。這些代謝通路相關(guān)的基因分布于已知的284個(gè)代謝通路,其中富集最多的10條通路分別是核糖體,內(nèi)質(zhì)網(wǎng)蛋白加工,嘌呤代謝,RNA轉(zhuǎn)運(yùn),剪接體,氧化磷酸化,碳代謝,嘧啶代謝,cAMP信號(hào)通路,真核生物的核糖體生物發(fā)生,真核生物中的核糖體生物合成(表1)。

表1 KEGG數(shù)據(jù)庫中單基因的十大代謝途徑富集

2.4 青海草原毛蟲Unigenes的KOG分類

對(duì)有注釋的基因進(jìn)行直系同源分類(圖3),9 521條unigenes共得到10 628條注釋,分為26類,其中,一般功能預(yù)測(cè)和翻譯后修飾的基因最多,分別為1 487條(13.99%) 和1 114條(10.48%);其余依次是T類信號(hào)轉(zhuǎn)導(dǎo)機(jī)制,有1 005條(9.46%);J類翻譯、核糖體結(jié)構(gòu)和生物合成,有864條(8.13%);K類轉(zhuǎn)錄,有644條(6.06%);U類細(xì)胞內(nèi)運(yùn)輸、分泌和囊泡轉(zhuǎn)運(yùn),有639條(6.01%);S類功能未知,有630條(5.93%);A類核酸處理與修飾,有598條(5.63%)。其余的都低于500條,X類未命名蛋白質(zhì),有1條(0.01%)。

圖3 青海草原毛蟲unigenes KOG分類注釋

2.5 青海草原毛蟲轉(zhuǎn)錄組中SSR分布特征

對(duì)青海草原毛蟲轉(zhuǎn)錄組中的63 335條unigenes的數(shù)據(jù)利用MISA軟件進(jìn)行簡單重復(fù)序列(SSR)位點(diǎn)搜索?？偣舱业?2 597個(gè)微衛(wèi)星位點(diǎn),分布于9 851條unigenes中,其微衛(wèi)星出現(xiàn)頻率(含有SSR的unigenes數(shù)量與總unigenes數(shù)量之比)是15.55%,其中含有1個(gè)以上SSR位點(diǎn)的獨(dú)立基因序列有1 852條,以復(fù)合型形式存在的SSR有801個(gè),SSR的分布頻率(SSR個(gè)數(shù)與總unigene數(shù)量之比)為19.89%(表2)。這些SSR基序包含1～5 bp的串聯(lián)重復(fù)序列。

表2 青海草原毛蟲轉(zhuǎn)錄組SSR分布及特征

2.6 青海草原毛蟲轉(zhuǎn)錄組序列中SSR位點(diǎn)重復(fù)類型分析

從青海草原毛蟲轉(zhuǎn)錄組數(shù)據(jù)中可以得到,單核苷酸重復(fù)是最主要的重復(fù),占SSR總數(shù)的71.83%。二核苷酸重復(fù)居其次,占SSR總數(shù)的14.27%。三、四核苷酸重復(fù),分別占SSR總數(shù)的10.59%和2.99%。五核苷酸和六核苷酸重復(fù)相較于前幾種核苷酸的含量很少,分別占SSR總數(shù)的0.29%和0.0.3%。出現(xiàn)頻率最高的是10次重復(fù)基元,共5 265個(gè)SSR位點(diǎn),占總SSR位點(diǎn)的41.77%(表3)。

表3 青海草原毛蟲轉(zhuǎn)錄組數(shù)據(jù)的SSR重復(fù)類型分布

對(duì)5～99次SSR基元重復(fù)次數(shù)進(jìn)行分析結(jié)果表明(圖4),SSR位點(diǎn)重復(fù)在5～10次的有8 335個(gè),占總SSR位點(diǎn)的66.17%;重復(fù)次數(shù)在11～20次的有3 896個(gè),占總SSR位點(diǎn)的30.93%;重復(fù)次數(shù)在21～30次的有267個(gè),占總SSR位點(diǎn)的2.12%;重復(fù)次數(shù)在31～40次的有53個(gè),占總SSR位點(diǎn)的0.42%;重復(fù)次數(shù)在41～50次的有20個(gè),占總SSR位點(diǎn)的0.16%;重復(fù)次數(shù)在51～99次的有26個(gè),占總SSR位點(diǎn)的0.21%。對(duì)單核苷酸重復(fù),二核苷酸重復(fù)和三核苷酸重復(fù)的分析發(fā)現(xiàn),SSR的數(shù)量與重復(fù)次數(shù)成反比出現(xiàn),隨著重復(fù)次數(shù)的增加,SSR的數(shù)量反而在減少。在圖中未呈現(xiàn)的四、五、六核苷酸重復(fù)具有與單、二、三核苷酸重復(fù)相同的趨勢(shì)表現(xiàn)。

圖4 青海草原毛蟲轉(zhuǎn)錄組中SSR基元重復(fù)次數(shù)

2.7 青海草原毛蟲轉(zhuǎn)錄組基元類型分析

青海草原毛蟲轉(zhuǎn)錄組SSR基元類型結(jié)果顯示(圖5和表4),有60個(gè)SSR重復(fù)基元類型,共獲得12 597個(gè)SSR位點(diǎn)。其中,2種類型的基元在單核苷酸重復(fù)中呈現(xiàn),(A/T)n是最主要的基元類型,形成8990個(gè)SSR位點(diǎn),占同類基元SSR位點(diǎn)的99.36%,占總SSR位點(diǎn)的71.37%。4種類型的基元在二核苷酸重復(fù)中呈現(xiàn),(AC/GT)n為主要優(yōu)勢(shì)基元,形成763個(gè)SSR位點(diǎn)(42.46%,6.06%),其次為(AT/AT)n和(AG/CT)n基元,分別形成619個(gè)(34.45%,4.91%)和370個(gè)SSR位點(diǎn)(20.59%,2.94%)。10種類型的基元在三核苷酸重復(fù)中呈現(xiàn),(AAT/ATT)n是最主要的優(yōu)勢(shì)基元,有533個(gè)SSR位點(diǎn)(39.96%,4.23%);其次為(ATC/ATG)n,(AAC/GTT)n,(AAG/CTT)n。分別形成236個(gè)(17.69%,1.87%),142個(gè)(10.64%,1.13%),114個(gè)(8.55%,0.91%)SSR位點(diǎn)。16種類型的基元在四核苷酸重復(fù)中呈現(xiàn),(AAAT/ATTT)n是最主要的優(yōu)勢(shì)基元,有170個(gè)SSR位點(diǎn)(45.09%,1.35%);其他優(yōu)勢(shì)基元分別為(ACAT/ATGT)n、(AGAT/ATCT)n,分別形成57個(gè)(15.12%,0.45%)和44個(gè)(11.67%,0.35%)SSR位點(diǎn)。在五核苷酸重復(fù)中有22個(gè)重復(fù)基元類型,最主要的優(yōu)勢(shì)基元是(AAAAT/ATTTT)n,(AACGT/ACGTT)n,在這兩個(gè)基元中都形成了4個(gè)(11.11%,0.03%)SSR位點(diǎn)。六核苷酸重復(fù)有5個(gè)重復(fù)基元類型,全都形成1個(gè)SSR位點(diǎn)(20%,0.01%)。

表4 青海草原毛蟲轉(zhuǎn)錄組不同重復(fù)基元SSR位點(diǎn)數(shù)及頻率

圖5 青海草原毛蟲轉(zhuǎn)錄組不同重復(fù)類型優(yōu)勢(shì)基元SSR位點(diǎn)數(shù)分布

2.8 青海草原毛蟲SSRs引物驗(yàn)證

本研究從被篩選出的微衛(wèi)星引物中隨機(jī)選取25對(duì),對(duì)5個(gè)不同樣品(海北州同一地點(diǎn)的5個(gè)青海草原毛蟲個(gè)體重復(fù))的青海草原毛蟲DNA樣本進(jìn)行PCR擴(kuò)增,結(jié)果11對(duì)引物在每個(gè)樣品中擴(kuò)增取得的片段與預(yù)期的片段幾乎一致。有效擴(kuò)增率為44%(表5,圖6)。

表5 青海草原毛蟲部分EST-SSR篩選引物信息

圖6 青海草原毛蟲5個(gè)樣品25個(gè)SSR位點(diǎn)引物的PCR擴(kuò)增產(chǎn)物電泳結(jié)果

3 討論

伴隨著高通量測(cè)序技術(shù)的快速發(fā)展,對(duì)昆蟲轉(zhuǎn)錄組學(xué)的研究變得非常廣泛[25]。本研究測(cè)序獲得unigenes序列總長度為69 070 017 bp,平均長度為1 091 bp,N50長度為1 714 bp,N90長度為440 bp。由于N50的值較大,因此其組裝完整性較高,保證了轉(zhuǎn)錄組分析的準(zhǔn)確性。

本研究發(fā)現(xiàn),青海草原毛蟲轉(zhuǎn)錄組數(shù)據(jù)庫中存在大量的微衛(wèi)星,且種類豐富,核苷酸的重復(fù)類型出現(xiàn)了6種。單核苷酸重復(fù)中最主要的重復(fù)類型是A/T,這與對(duì)印度谷螟(Plodiainterpunctella)[2]、黏蟲(Mythimnaseparata)[26]、褐飛虱(Nilaparuatalugens)[27]SSR位點(diǎn)特征的研究結(jié)果基本相同。但本次試驗(yàn)結(jié)果與對(duì)黃粉蟲(Tenebriomolitor)[28]、麥紅吸漿蟲(Sitodiplosismosellana)[29]、沙蔥葉甲(Galerucadaurica)[30]等昆蟲的SSR位點(diǎn)特征研究結(jié)果相比較存在差異。這些研究發(fā)現(xiàn)核苷酸的重復(fù)類型以單核苷酸為主,三核苷酸次之。本研究的結(jié)果是單核苷酸重復(fù)是最主要的重復(fù)類型,二核苷酸次之。該結(jié)果與意大利微蝗(Calliptamusitalicus)[31]、梨小食心蟲(Grapholithamolesta)[32]、沙棘木蠹蛾(Holcocerushippophaecolus)的研究結(jié)果類似。通常情況下,在ESTs里,單核苷酸為主要重復(fù)類型,三核苷酸重復(fù)居于其次。原因在于,與編碼區(qū)其他的重復(fù)基元類型相比,三核苷酸核心基元更加穩(wěn)定,基本不會(huì)出現(xiàn)編碼框滑動(dòng)突變現(xiàn)象[33]。本研究中出現(xiàn)差異,可能是不同樣品之間本身存在差異,從而導(dǎo)致結(jié)果不同。

為驗(yàn)證青海草原毛蟲SSR的穩(wěn)定性,本研究利用其轉(zhuǎn)錄組數(shù)據(jù)設(shè)計(jì)出2 714對(duì)引物,從中隨機(jī)選出25對(duì)引物進(jìn)行擴(kuò)增試驗(yàn),結(jié)果表明11對(duì)引物在每個(gè)單個(gè)樣品中都能成功擴(kuò)增,并得到預(yù)期條帶。這些位點(diǎn)的多態(tài)性目前還不明確,需要進(jìn)一步試驗(yàn)評(píng)估,但它的數(shù)量足以開展青海草原毛蟲種群遺傳學(xué)相關(guān)研究。導(dǎo)致不同樣品的EST-SSR位點(diǎn)擴(kuò)增失敗的原因有很多,如:由于基因組DNA中的內(nèi)含子在轉(zhuǎn)錄后沒表現(xiàn)出來,對(duì)擴(kuò)增過程造成干擾從而導(dǎo)致預(yù)期的條帶與實(shí)際擴(kuò)增的條帶相比,會(huì)出現(xiàn)擴(kuò)增條帶遠(yuǎn)遠(yuǎn)大于預(yù)期的現(xiàn)象;在基因組中EST-SSR的表達(dá)豐富度太低,導(dǎo)致產(chǎn)物的濃度低,無法檢測(cè)[2]。

微衛(wèi)星的重要作用體現(xiàn)在基因調(diào)控、蛋白功能、種群遺傳結(jié)構(gòu)、種群擴(kuò)散及遷移等方面[34]。引物的開發(fā)是微衛(wèi)星應(yīng)用的前提,傳統(tǒng)的方法開發(fā)成本太高,不合時(shí)宜。通過在轉(zhuǎn)錄組高通量數(shù)據(jù)中直接發(fā)掘微衛(wèi)星并設(shè)計(jì)引物的方法,大大減少了開發(fā)成本,很大程度上對(duì)微衛(wèi)星的開發(fā)進(jìn)程起到推動(dòng)作用[30]。本研究在青海草原毛蟲轉(zhuǎn)錄組數(shù)據(jù)庫中共篩出12 597個(gè)SSR位點(diǎn),并對(duì)其分布情況及可用性進(jìn)行了詳細(xì)敘述,這些SSR的發(fā)掘,對(duì)于之后研究青海草原毛蟲的遺傳多樣性和功能基因組學(xué)奠定了堅(jiān)實(shí)的基礎(chǔ),也更進(jìn)一步表明通過高通量測(cè)序技術(shù)發(fā)掘昆蟲的SSR引物是一種符合人類需求的研究方法。

4 結(jié)論

本研究通過KOG,GO注釋和KEGG通路三大數(shù)據(jù)庫分析,發(fā)現(xiàn)unigenes注釋主要集中于一般功能預(yù)測(cè)、細(xì)胞進(jìn)程和碳水化合物代謝。本研究共篩出12 597個(gè)SSR位點(diǎn),發(fā)生頻率為19.89%,SSR位點(diǎn)數(shù)量豐富,(A/T)n,(AC/GT)n分別是單核苷酸重復(fù)和二核苷酸的優(yōu)勢(shì)基元,SSR數(shù)量隨重復(fù)次數(shù)增加而降低,重復(fù)次數(shù)類型隨基元序列長度增長而減少。從1 556條青海草原毛蟲unigenes中成功設(shè)計(jì)出2 714對(duì)青海草原毛蟲SSR引物,隨機(jī)挑出25對(duì)引物進(jìn)行PCR驗(yàn)證,有11對(duì)引物擴(kuò)增出目的DNA片段。本研究基于轉(zhuǎn)錄組數(shù)據(jù)成功篩選出青海草原毛蟲微衛(wèi)星位點(diǎn),獲得的青海草原毛蟲SSR位點(diǎn)為進(jìn)一步研究其種群遺傳學(xué)和種群擴(kuò)散及遷移奠定基礎(chǔ)。