劉靖 魏俞虹 肖榕

摘 要：擬環(huán)紋豹蛛Pardosa pseudoannulata是稻田害蟲重要的自然天敵，在綠色防控中有重要作用。本研究基于擬環(huán)紋豹蛛轉(zhuǎn)錄組數(shù)據(jù)，使用微衛(wèi)星鑒定工具軟件（MISA） 挖掘和分析其簡單序列重復(fù)（SSR）位點(diǎn)。結(jié)果表明：擬環(huán)紋豹蛛轉(zhuǎn)錄組Unigene序列中共檢測(cè)到6474個(gè)SSR位點(diǎn)，平均長度為12.36 bp，分布于4411條Unigene上。單核苷酸重復(fù)類型為擬環(huán)紋豹蛛轉(zhuǎn)錄組SSR位點(diǎn)的主導(dǎo)基序類型，共5025個(gè)，占總SSR位點(diǎn)的77.62%;在擬環(huán)紋豹蛛轉(zhuǎn)錄組序列識(shí)別的6474個(gè)SSR中共發(fā)現(xiàn)了53種重復(fù)基元，從單核苷酸到五核苷酸重復(fù)基元分別有2種、6種、27種、15種和3種，其中（A/T）基元類型占明顯優(yōu)勢(shì)，占總SSR的74.44%。擬環(huán)紋豹蛛轉(zhuǎn)錄組SSR位點(diǎn)豐富，有較高的多態(tài)性潛能，可作為開發(fā)SSR標(biāo)記的有效來源，為擬環(huán)紋豹蛛的遺傳多樣性分析、功能基因的開發(fā)利用等方面的研究奠定基礎(chǔ)。

關(guān)鍵詞：擬環(huán)紋豹蛛;轉(zhuǎn)錄組;SSR;重復(fù)基元

中圖分類號(hào)：Q346+.5;Q959.1

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1008-0457（2021）04-0031-07

國際DOI編碼：10.15958/j.cnki.sdnyswxb.2021.04.005

Abstract：The wolf spider Pardosa pseudoannulata is an important biological control agent of rice pests and plays an important role in green prevention and control.In this study，based on the transcriptome data of P.pseudoannulata，we used MISA （microsatellite identification tool） to analyze its SSR （simple sequence repeat） loci.The results showed that a total of 6474 SSR loci were detected in the Unigene sequence of the transcriptome data，with an average length of 12.36 bp，distributed on 4411 Unigenes.The Mono-nucleotide repeat type were the dominanttype，its number was 5025，accounting for 77.62% of the total SSR loci.The 6474 SSR lociwere identified 53 repetitive motifs，there were 2，6，27，15 and 3 repetitive motifs from Mono-nucleotide to Hexa-nucleotide.（A/T）was dominant repetitive motif，accounting for 74.44% of the total SSR.The SSR loci of transcriptome in P.pseudoannulata were rich and high polymorphism potential. Our results can be used as an effective source for the development of SSR marker and provide foundation for research of genetic diversity and development functional genes of P.pseudoannulata.

Keywords：Pardosa pseudoannulata;transcriptome;SSR;repeat motif

擬環(huán)紋豹蛛Pardosa pseudoannulata屬于蛛形綱Arachnida蜘蛛目Araneae狼蛛科Lycosidae豹蛛屬Pardosa，在我國分布較廣，是稻田生態(tài)系統(tǒng)中優(yōu)勢(shì)種蜘蛛，尤其在我國南方的廣東[1]、湖南[2]、湖北[3]、云南[4]、貴州[5]等地區(qū)常年可見。擬環(huán)紋豹蛛屬于中型蜘蛛，成蛛體長可達(dá)到13 mm，雌雄蛛均可多次交配，雌蛛每次平均產(chǎn)卵156粒，孵化的幼蛛群集在雌蛛背面3～7 d后才離開，幼蛛期47～158 d，成蛛期121～236 d，雌蛛較雄蛛長21～62 d[2，6-7]。擬環(huán)紋豹蛛是稻田害蟲的重要捕食性天敵，其空間生態(tài)位寬[6，8]，能夠在水面、稻株上、地面上等多種環(huán)境捕食褐飛虱Nilaparvata lugens、白背飛虱Sogatella furcifera、稻縱卷葉螟Cnaphalocrocis medinalis、黑尾葉蟬Nephotettix bipunctatus等害蟲。在稻田中捕食飛虱和葉蟬，90%是狼蛛科蜘蛛的功勞[9-13]。擬環(huán)紋豹蛛的性情兇猛，專捕活蟲，即使飽食，也會(huì)把獵物咬死后才休止，食量也大，2齡幼蛛每頭每天平均捕食黑腹果蠅Drosophila melanogaster 4.36頭，成蛛每頭每天捕食水蠅Hydrellia sp.可達(dá)30頭[8，14]，雌成蛛對(duì)亞洲玉米螟Ostrinia furnacalis 2齡幼蟲的捕食上限為116頭[15]，王智等[16]研究表明低劑量農(nóng)藥對(duì)擬環(huán)紋豹蛛捕食褐飛虱的能力有增強(qiáng)作用。本課題組前期對(duì)擬環(huán)紋豹蛛的求偶交配行為[17]、選擇行為[18]和溫度脅迫[19-20]和SNP位點(diǎn)[21]等方面進(jìn)行了研究。目前，擬環(huán)紋豹蛛轉(zhuǎn)錄組SSR位點(diǎn)分析未見報(bào)道。

微衛(wèi)星（microsatellite），又稱為簡單序列重復(fù)（simple sequence repeat，SSR），是指以1～6個(gè)核苷酸為重復(fù)單元串聯(lián)組成的重復(fù)序列，串聯(lián)重復(fù)10～60次組成的簡單重復(fù)序，微衛(wèi)星具有高多態(tài)性、共顯性、高可重復(fù)性、數(shù)量豐富等優(yōu)點(diǎn)[22-23]。因此，微衛(wèi)星在遺傳分析中可以作為一種理想的分子標(biāo)記，已被廣泛應(yīng)用于生物遺傳研究中，在苦蕎Fagopyrum tartaricum[24]、密花香薷Elsholtzia densa[25]、密葉紅豆杉Taxus fuana[26]、金佛山方竹Chimonobambusa utilis[27]等植物;麥紅吸漿蟲Sitodiplosis mosellana[28]、蚜蟲[29]、小胸鱉甲Microdera punctipennis[30]、桔小實(shí)蠅Bactrocera dorsalis[31]、溝渠豹蛛Pardosa laura[32]等動(dòng)物中進(jìn)行了報(bào)道。本文基于擬環(huán)紋豹蛛轉(zhuǎn)錄組測(cè)序數(shù)據(jù)，對(duì)其SSR特征進(jìn)行分析，旨在為SSR標(biāo)記在該物種的遺傳多樣性分析、功能基因的開發(fā)利用等方面提供參考。

1?材料與方法

1.1?供試蟲源

擬環(huán)紋豹蛛的亞成蛛采集于廣州市華南農(nóng)業(yè)大學(xué)試驗(yàn)田（113.7°E，23.17°N）。亞成蛛采回后單頭放入玻璃試管（ 長100 mm，直徑25 mm） 中飼養(yǎng)，試管底部用一塊蘸水的海綿保濕，管口塞有海綿塞，防止蜘蛛逃逸，置于人工氣候箱內(nèi)（25±0.5）℃，相對(duì)濕度RH為60%±10%，光周期為14 L：10 D飼養(yǎng)。每周飼喂2次，每次飼喂4～6頭活家蠅Musca domestica成蟲，同時(shí)給試管底部的海綿添加足量的水以保持濕潤。每隔1 d觀察記錄擬環(huán)紋豹蛛蛻皮及成熟情況，成熟兩周后的蜘蛛用液氮速凍后置-80 ℃冰箱備用。

1.2?SSR位點(diǎn)的檢測(cè)

當(dāng)蜘蛛成熟14 d后，使用TRIzol試劑盒提取蜘蛛的總RNA，經(jīng)質(zhì)量檢測(cè)合格的樣品，送北京百邁客生物科技有限公司利用Illumina HiseqTM 2500平臺(tái)進(jìn)行轉(zhuǎn)錄組測(cè)序，測(cè)序結(jié)果用Trinity軟件進(jìn)行數(shù)據(jù)組裝獲得Unigene。利用MISA（MIcroSAtellite identification tool，http：//pgrc.ipk-gatersleben.de/misa/misa.html）軟件對(duì)擬環(huán)紋豹蛛轉(zhuǎn)錄組中的Unigene進(jìn)行SSR位點(diǎn)檢測(cè)。檢測(cè)的重復(fù)序列從單堿基重復(fù)至六堿基重復(fù)，篩選參數(shù)：單核苷酸序列重復(fù)次數(shù)不低于10次，二核苷酸序列重復(fù)次數(shù)不低于6次，三核苷酸序列、四核苷酸序列、五核苷酸序列和六核苷酸序列重復(fù)次數(shù)則不低于5次。若2個(gè)SSR間距小于100 bp，則組成一個(gè)復(fù)合SSR位點(diǎn)。

2?結(jié)果與分析

2.1?擬環(huán)紋豹蛛轉(zhuǎn)錄組分布概況

組裝后的擬環(huán)紋豹蛛轉(zhuǎn)錄組數(shù)據(jù)，共獲得67 725條Unigene，篩選長度1 kb以上的13 762條Unigene，利用MISA軟件分析發(fā)現(xiàn)，SSR重復(fù)序列總共有6474個(gè)，分布在4411條Unigene上，擬環(huán)紋豹蛛轉(zhuǎn)錄組Unigene序列平均每4581 bp就出現(xiàn)一個(gè)SSR位點(diǎn)。包含1個(gè)以上SSR位點(diǎn)數(shù)的Unigene數(shù)目1347條（表1）?？梢?，擬環(huán)紋豹蛛轉(zhuǎn)錄組所包含SSR位點(diǎn)較豐富，且分布較廣泛。

2.2?擬環(huán)紋豹蛛轉(zhuǎn)錄組SSR位點(diǎn)不同重復(fù)類型分析

對(duì)擬環(huán)紋豹蛛轉(zhuǎn)錄組搜索的SSR位點(diǎn)進(jìn)行了核苷酸重復(fù)類型分類，共有5種重復(fù)類型（表2）。不同重復(fù)類型形成的SSR位點(diǎn)數(shù)目相差較大，其中單核苷酸重復(fù)位點(diǎn)最多，為5025個(gè)，占總SSR位點(diǎn)的77.62%，總長度為62 117 bp，平均每477 bp出現(xiàn)1個(gè)單核苷酸重復(fù)序列，單核苷酸平均長度為12.36 bp;二核苷酸重復(fù)SSR位點(diǎn)，為985個(gè)，占總SSR位點(diǎn)的15.21%，總長度為13 868 bp，平均每2138 bp出現(xiàn)1個(gè)二核苷酸重復(fù)序列，二核苷酸平均長度為14.08 bp;三核苷酸重復(fù)SSR位點(diǎn)，為421個(gè)，占總SSR位點(diǎn)的6.50%，總長度為6918 bp，平均每4287 bp出現(xiàn)1個(gè)三核苷酸重復(fù)序列，三核苷酸平均長度為16.43 bp;四核苷酸重復(fù)SSR位點(diǎn)，為40個(gè)，占總SSR位點(diǎn)的0.62%，總長度為824 bp，平均每35 990 bp出現(xiàn)1個(gè)四核苷酸重復(fù)序列，四核苷酸平均長度為20.60 bp;五核苷酸重復(fù)SSR位點(diǎn)最少，為3個(gè)，占總SSR位點(diǎn)的0.05%，總長度為90 bp，平均每329 506 bp出現(xiàn)1個(gè)五核苷酸重復(fù)序列，五核苷酸平均長度為30.00 bp。不同重復(fù)類型形成SSR頻率差異較大，其中單核苷酸重復(fù)類型SSR發(fā)生頻率最高，每100條Unigene有36.51個(gè)SSR位點(diǎn)出現(xiàn);五核苷酸重復(fù)類型SSR發(fā)生頻率最低，100條Unigene僅有0.02個(gè)SSR位點(diǎn)出現(xiàn)（表2）。

2.3?擬環(huán)紋豹蛛轉(zhuǎn)錄組SSR基元類型分析

在擬環(huán)紋豹蛛轉(zhuǎn)錄組序列識(shí)別的6474個(gè)SSR中共發(fā)現(xiàn)了53種重復(fù)基元，其中單核苷酸重復(fù)基元有2種、二核苷酸重復(fù)基元有6種、三核苷酸重復(fù)基元有27種、四核苷酸重復(fù)基元有15種、五核苷酸重復(fù)基元有3種。不同基元類型形成SSR位點(diǎn)數(shù)目存在差異，從總體來看，隨著SSR基元堿基數(shù)目的增加，SSR位點(diǎn)數(shù)呈下降趨勢(shì)。其中單核苷酸基元（A/T）發(fā)生頻率最高為35.02%，其次是（AT/AT）的2.81%、（TA/TA）的2.03%、（C/G）的1.50%、（AC/GT）的0.80%等（表3）。

在單核苷酸重復(fù)基元中，（A/T）基元類型占明顯優(yōu)勢(shì)，共形成4819個(gè)SSR位點(diǎn)，占單核苷酸SSR位點(diǎn)的95.90%，總SSR位點(diǎn)的74.44%，屬于單堿基的優(yōu)勢(shì)重復(fù)基元（圖1，表3）。在二核苷酸重復(fù)基元中（AT/AT）和（TA/TA）基元類型占明顯優(yōu)勢(shì)。（AT/AT）基元類型共形成387個(gè)SSR位點(diǎn)，占二核苷酸SSR位點(diǎn)的39.29%，占總SSR位點(diǎn)的5.98%;其次是（TA/TA）基元類型，共形成280個(gè)SSR位點(diǎn)，占二核苷酸SSR位點(diǎn)的28.43%，占總SSR位點(diǎn)的4.32%（圖1，表3）。在三核苷酸重復(fù)基元中（AAT/ATT）、（TCA/TGA）、（TAA/TTA）基元類型占明顯優(yōu)勢(shì)。（AAT/ATT）基元類型形成46個(gè)SSR位點(diǎn)，占三核苷酸SSR位點(diǎn)的10.93%，占總SSR位點(diǎn)的0.71%;（TCA/TGA）基元類型形成39個(gè)SSR位點(diǎn)，占三核苷酸SSR位點(diǎn)的9.26%，占總SSR位點(diǎn)的0.60%;（TAA/TTA）基元類型形成37個(gè)SSR位點(diǎn)，占三核苷酸SSR位點(diǎn)的8.79%，占總SSR位點(diǎn)的0.57%（圖1，表3）。在四核苷酸重復(fù)基元中（TAAA/TTTA）和（AAAT/ATTT）基元類型占明顯優(yōu)勢(shì)。（TAAA/TTTA）基元類型共形成8個(gè)SSR位點(diǎn)，占四核苷酸SSR位點(diǎn)的20.00%，占總SSR位點(diǎn)的0.12%;其次是（AAAT/ATTT）基元類型，共形成7個(gè)SSR位點(diǎn)，占四核苷酸SSR位點(diǎn)的17.50%，占總SSR位點(diǎn)的0.11%（圖1，表3）。在五核苷酸重復(fù)基元中（AACCT/AGGTT）、（TTTTA/TAAAA）、（ATATA/TATAT）基元類型，每一種類型都只形成1個(gè)SSR位點(diǎn)，占五核苷酸SSR位點(diǎn)的33.33%，占總SSR位點(diǎn)的0.02%（圖1，表3）。

2.4?擬環(huán)紋豹蛛轉(zhuǎn)錄組SSR重復(fù)次數(shù)分析

擬環(huán)紋豹蛛SSR基元重復(fù)次數(shù)在不同基元類型間存在差異，不同重復(fù)次數(shù)所形成的SSR位點(diǎn)差異較大，每種類型構(gòu)成的SSR位點(diǎn)中，最小重復(fù)次數(shù)的SSR位點(diǎn)最多，且隨著基元堿基數(shù)目的增加，重復(fù)次數(shù)類型呈下降趨勢(shì)（圖2）。單核苷酸重復(fù)基元的重復(fù)次數(shù)類型最多，共有15種重復(fù)次數(shù)，不同重復(fù)次數(shù)所形成的SSR位點(diǎn)數(shù)量差異較大，最小重復(fù)次數(shù)為10次，形成2210個(gè)SSR位點(diǎn)，占單核苷酸SSR位點(diǎn)的43.98%，最大重復(fù)次數(shù)為24次，僅形成1個(gè)SSR位點(diǎn)，占單核苷酸SSR位點(diǎn)的0.02%。二核苷酸重復(fù)基元共有8種重復(fù)次數(shù)類型，其中最小重復(fù)次數(shù)為6次（505個(gè)，51.26%），最大重復(fù)次數(shù)為13次（1個(gè)，0.10%）。三核苷酸重復(fù)基元共有6種重復(fù)次數(shù)類型，其中最小重復(fù)次數(shù)為5次（278個(gè)，66.03%），最大重復(fù)次數(shù)為10次（1個(gè)，0.24%）。四核苷酸重復(fù)基元共有2種重復(fù)次數(shù)類型，最小重復(fù)次數(shù)為5次（34個(gè)，85%），最大重復(fù)次數(shù)為6次（位點(diǎn)數(shù)：6個(gè)，15%）。五核苷酸重復(fù)基元的共有3種重復(fù)次數(shù)類型，其中最小重復(fù)次數(shù)為5次（1個(gè)，33.33%），最大重復(fù)次數(shù)為7次（1個(gè)，33.33%）。

2.5?擬環(huán)紋豹蛛轉(zhuǎn)錄組SSR長度分析

影響SSR多態(tài)性的主要因素是SSR序列的長度，SSR序列長度與其多態(tài)性成正比。SSR序列長度<12 bp時(shí)，SSR多態(tài)性很低;12 bp≤SSR序列長度<20 bp時(shí)，SSR多態(tài)性良好;SSR序列長度≥20 bp時(shí)，SSR多態(tài)性較高[27，33]。對(duì)擬環(huán)紋豹蛛轉(zhuǎn)錄組SSR長度進(jìn)行分析，SSR長度在10～35 bp之間，平均長度為12.95 bp（表2，圖3）。其中數(shù)量最多的SSR序列長度為12～19 bp，占比為47.28%;其次是長度小于12 bp SSR序列，占比為46.18%;然后是SSR序列長度為20～35 bp，占比為6.53%，總體來看，SSR多態(tài)性屬于良好（圖3）。

3?結(jié)論與討論

本研究基于擬環(huán)紋豹蛛轉(zhuǎn)錄組數(shù)據(jù)，共獲得67 725條Unigene，篩選長度1 kb以上的13 762條Unigene，利用MISA軟件共檢測(cè)到6474個(gè)SSR位點(diǎn)，分布于4411條Unigene上，SSR出現(xiàn)頻率為47.04%。擬環(huán)紋豹蛛轉(zhuǎn)錄組SSR出現(xiàn)頻率明顯高于荔枝蒂蛀蟲Conopomorpha sinensis（15.25%）[34]、印度谷螟Plodia interpunctella（8.52%）[35]、小胸鱉甲M.punctipennis（7.94%）[30]、麥紅吸漿蟲S.mo-sellana（3.49%）[28]等昆蟲。不同物種間SSR位點(diǎn)出現(xiàn)頻率的差異，可能是與不同物種基因組特點(diǎn)有關(guān)，也可能和測(cè)序深度、篩選SSR的軟件及檢測(cè)時(shí)選擇的參數(shù)等因素有關(guān)。如，本研究與溝渠豹蛛P.laura（80.11%）[32]的研究相比，研究對(duì)象都屬于豹蛛屬蜘蛛，設(shè)定參數(shù)也相同，均是對(duì)大于1 kb的Unigene進(jìn)行SSR搜索，但檢測(cè)到的SSR數(shù)量和SSR出現(xiàn)的頻率存在明顯的差異。

擬環(huán)紋豹蛛轉(zhuǎn)錄組搜索的SSR位點(diǎn)進(jìn)行了核苷酸重復(fù)類型分類，共有5種重復(fù)類型，單核苷酸重復(fù)位點(diǎn)為5025個(gè)、二核苷酸重復(fù)SSR位點(diǎn)為985個(gè)、三核苷酸重復(fù)SSR位點(diǎn)為421個(gè)、四核苷酸重復(fù)SSR位點(diǎn)為40個(gè)、五核苷酸重復(fù)SSR位點(diǎn)為3個(gè)。擬環(huán)紋豹蛛轉(zhuǎn)錄組SSR位點(diǎn)的主導(dǎo)基序類型是單核苷酸重復(fù)類型，占總SSR位點(diǎn)的77.62%，高于馬鈴薯甲蟲Leptinotarsa decemlineata（63.23%）[36]、梨網(wǎng)蝽Stephanitis nashi（40.41%）[37]、麥紅吸漿蟲S.mosellana（37.07%）[28]、小胸鱉甲M.punctipennis（47.17%）[30]等昆蟲，與溝渠豹蛛P.laura（72.67%）[32]相差不大，低于窄足真蚋Simulium（Eusimulium）angustipes（87.05%）[38]。擬環(huán)紋豹蛛轉(zhuǎn)錄組中共發(fā)現(xiàn)了53種重復(fù)基元，從單核苷酸到五核苷酸重復(fù)基元分別有2種、6種、27種、15種和3種，其中單核苷酸（A/T）基元類型占總SSR的74.44%，占明顯優(yōu)勢(shì)，這與溝渠豹蛛P.laura[32]、馬鈴薯甲蟲L.decemlineata[36]、麥紅吸漿蟲S.mosellana[28]、小胸鱉甲M.punctipennis[30]等轉(zhuǎn)錄組SSR分析研究結(jié)果相似。二核苷酸重復(fù)基元中，（AT/AT）是優(yōu)勢(shì)基元，分別占總SSR的5.98%，這與溝渠豹蛛P.laura[32]分析研究結(jié)果相似。但在擬環(huán)紋豹蛛轉(zhuǎn)錄組SSR中，不含二核苷酸重復(fù)基元（GC/CG），這與桔小實(shí)蠅B.dorsalis[31]及褐飛虱N.lugens[39]的研究結(jié)果相似，（GC/CG）是非常稀有的重復(fù)基元，其數(shù)量幾乎接近零，推測(cè)（GC/CG）重復(fù)基元在不同物種的基因組中都可能是相當(dāng)稀有的SSR。三核苷酸重復(fù)基元中，（AAT/ATT）基元占優(yōu)勢(shì)，與褐飛虱（AAG/CTT）[39]、麥紅吸漿蟲（AAC/GTT）[28]、小胸鱉甲（CCG/GCC）[30]不同，與溝渠豹蛛[32]相同，（AAT/ATT）基元占優(yōu)勢(shì)，可能是由于高滑變率及其作為終止密碼子的延伸區(qū)引起的[40]。

綜上所述，本研究基于擬環(huán)紋豹蛛轉(zhuǎn)錄組數(shù)據(jù)篩選得到其SSR位點(diǎn)信息，對(duì)不同類型位點(diǎn)的數(shù)量、比例、頻率與重復(fù)次數(shù)，以及對(duì)不同的重復(fù)基元類型的數(shù)量和比例等進(jìn)行了分析，以期為研究擬環(huán)紋豹蛛的遺傳多樣性分析、功能基因的開發(fā)利用等方面的研究奠定基礎(chǔ)。

參?考?文?獻(xiàn)：

[1]?劉雨芳，古德祥，張古忍.廣東省稻田生態(tài)系統(tǒng)中蜘蛛群落多樣性研究[J].蛛形學(xué)報(bào)，2003，12（1）：27-31.

[2]?王智.擬環(huán)紋豹蛛的生物生態(tài)學(xué)研究[J].昆蟲學(xué)報(bào)，2007，50（9）：927-932.

[3]?陳博，文樂雷，趙菊鵬，等.七種稻田常用殺蟲劑對(duì)擬環(huán)紋豹蛛的室內(nèi)風(fēng)險(xiǎn)評(píng)價(jià)[J].植物保護(hù)學(xué)報(bào)，2017，44（6）：1059-1060.

[4]?楊庭榜，楊自忠，石愛民.洱海周邊擬環(huán)紋豹蛛和星豹蛛卵袋中寄生蜂種類及其寄主選擇[J].四川動(dòng)物，2012，31（6）：858-861.

[5]?劉靖，孫良玉，付迪，等.貴州省稻田蜘蛛多樣性調(diào)查及優(yōu)勢(shì)種分析[J].應(yīng)用昆蟲學(xué)報(bào)，2021，58（1）：142-157.

[6]?王洪全，顏亨梅，楊海明.中國稻田蜘蛛群落結(jié)構(gòu)研究初報(bào)[J].蛛形學(xué)報(bào)，1999，8（2）：95-105.

[7]?尹長民，彭賢錦，顏亨梅.湖南動(dòng)物志，蜘蛛類（上下冊(cè)）[M].長沙：湖南科學(xué)技術(shù)出版社，2012：852-854.

[8]?王智，宋大祥，朱明生.稻田蜘蛛和害蟲的生態(tài)位研究[J].華南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2005，26（2）：47-51.

[9]?王洪全，顏亨梅，楊海明.中國稻田蜘蛛生態(tài)與利用研究[J].中國農(nóng)業(yè)科學(xué)，1996，29（5）：69-76.

[10]?王智，顏亨梅，呂志躍，等.稻田蜘蛛優(yōu)勢(shì)種對(duì)飛虱與葉蟬控制力的分析[J].生命科學(xué)研究，2001，5（1）：76-79.

[11]?吳小玲，王智.稻田蜘蛛優(yōu)勢(shì)種和主要害蟲的空間生態(tài)位研究[J].科技創(chuàng)新導(dǎo)報(bào)，2011，8（27）：148-149.

[12]?張古忍，張文慶，古德祥.稻田主要節(jié)肢類捕食者群落的結(jié)構(gòu)和動(dòng)態(tài)[J].中山大學(xué)學(xué)報(bào)論叢，1995，15（2）：33-40.

[13]?顏亨梅，王洪全，楊海明，等.湖南稻田蜘蛛群落生態(tài)特點(diǎn)及控蟲作用研究[J].蛛形學(xué)報(bào)，2002，11（1）：48-50.

[14]?王洪全，周家友，劉貴勻.擬環(huán)紋狼蛛的生物學(xué)研究[J].動(dòng)物學(xué)報(bào)，1982，28（1）：69-79.

[15]?田耀加，梁廣文，曾玲，等.擬環(huán)紋豹蛛對(duì)亞洲玉米螟卵及低齡幼蟲的捕食效應(yīng)[J].中國生物防治學(xué)報(bào)，2012，28（4）：478-483.

[16]?王智，宋大祥，朱明生.低劑量農(nóng)藥作用下擬環(huán)紋豹蛛對(duì)褐飛虱的功能反應(yīng)及搜尋行為[J].昆蟲學(xué)報(bào)，2006，49（2）：295-301.

[17]?肖榕，陳占起，彭宇，等.擬環(huán)紋豹蛛交配和同類相食行為影響因子分析[J].華中昆蟲研究，2012，8（8）：86-89.

[18]?肖榕，曹營帥，武韓，等.從蟲害誘導(dǎo)水稻揮發(fā)物中篩選對(duì)擬環(huán)紋豹蛛有引誘活性的化合物[J].植物保護(hù)學(xué)報(bào)，2018，45（5）：1021-1027.

[19]?XIAO R，WANG L，CAO Y S，et al.Transcriptome response to temperature stress in the wolf spider Pardosa pseudoannulata （Araneae：Lycosidae）[J].Ecology and Evolution，2016，6（11）：3540-3553.

[20]?SUN L Y，LIU J，LI Q，et al.Cloning and differential expression of three heat shock protein genes associated with thermal stress from the wolf spider Pardosa pseudoannulata （Araneae：Lycosidae）[J].Journal of Asia-Pacific Entomology，2021，24（1）：158-166

[21]?王祥銀，孫良玉，劉靖，等.擬環(huán)紋豹蛛響應(yīng)溫度脅迫SNP位點(diǎn)挖掘及其功能注釋分析[J].華北農(nóng)學(xué)報(bào)，2019，34（1）：331-339.

[22]?ESTOUP A，JARNE P，CORNUET J M.Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis[J].Molecular Ecology，2002，11（9）：1591-1604

[23]?黎東海，趙萍.基于轉(zhuǎn)錄組數(shù)據(jù)的齒緣刺獵蝽微衛(wèi)星分子標(biāo)記開發(fā)[J].昆蟲學(xué)報(bào)，2019，62（6）：694-702.

[24]?馬名川，劉龍龍，劉璋，等.苦蕎全基因組SSR位點(diǎn)特征分析與分子標(biāo)記開發(fā)[J].作物雜志，2021，47（1）：38-46.

[25]?富貴，劉玉萍，蘇旭.基于轉(zhuǎn)錄組數(shù)據(jù)的密花香薷SSR位點(diǎn)特征分析[J].西北植物學(xué)報(bào)，2021，41（4）：654-663.

[26]?申響保，朱妍潔，徐剛標(biāo).密葉紅豆杉SSR位點(diǎn)分布特征及分子標(biāo)記開發(fā)[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2021，41（4）：139-147.

[27]?朱瀟，代朝霞，劉艷江，等.金佛山方竹筍籜轉(zhuǎn)錄組SNP和SSR位點(diǎn)信息分析[J].分子植物育種，2021，6（3）：1-19.

[28]?段云，吳仁海，羅禮智，等.麥紅吸漿蟲唾腺EST-SSRs的信息分析及分子標(biāo)記篩選[J].昆蟲學(xué)報(bào)，2011，54（10）：1147-1154.

[29]?王永模，沈佐銳，高靈旺.微衛(wèi)星標(biāo)記及其在蚜蟲種群生物學(xué)研究中的應(yīng)用[J].昆蟲學(xué)報(bào)，2007，50（6）：621-627.

[30]?庫爾班·吐松，陳真，劉小寧，等.荒漠昆蟲小胸鱉甲轉(zhuǎn)錄組微衛(wèi)星的生物信息學(xué)分析[J].基因組學(xué)與應(yīng)用生物學(xué)，2017，36（3）：865-873.

[31]?魏丹丹，石俊霞，張夏瑄，等.基于轉(zhuǎn)錄組數(shù)據(jù)的桔小實(shí)蠅微衛(wèi)星位點(diǎn)信息分析[J].應(yīng)用生態(tài)學(xué)報(bào)，2014，25（6）：1799-1805.

[32]?寧青，王立華，戴余軍，等.基于轉(zhuǎn)錄組的溝渠豹蛛SSR信息分析[J].基因組學(xué)與應(yīng)用生物學(xué)，2020，39（12）：5464-5470.

[33]?TEMNYKH S，DECLERCK G，LUKASHOVA A，et al.Computational and experimental analysis of microsatellites in rice（Oryza sativa L.）：frequency，length variation，transposon associations，and genetic marker potential[J].Genome Res，2001，11（8） ：1441-1452.

[34]?孟翔，胡俊杰，李艷華，等.基于轉(zhuǎn)錄組數(shù)據(jù)的荔枝蒂蛀蟲SSR位點(diǎn)信息分析[J].環(huán)境昆蟲學(xué)報(bào)，2017，39（6）：1219-1224.

[35]?唐培安，陶冶心，薛昊，等.基于轉(zhuǎn)錄組數(shù)據(jù)的印度谷螟微衛(wèi)星位點(diǎn)分析[J].植物保護(hù)，2017，43（3）：43-48.

[36]?劉旸，付開赟，吐爾遜，等.馬鈴薯甲蟲基因組SSR位點(diǎn)分析及引物效率的驗(yàn)證[J].環(huán)境昆蟲學(xué)報(bào)，2018，40（3）：633-644.

[37]?謝瑾燕，周健，黃大莊.基于RNA-seq數(shù)據(jù)開發(fā)梨網(wǎng)蝽SSR引物研究[J].林業(yè)與生態(tài)科學(xué)，2019，34（3）：280-285.

[38]?郭歡，王剛，張樹田，等.基于RNA-seq數(shù)據(jù)的窄足真蚋SSR分子標(biāo)記開發(fā)[J].昆蟲學(xué)報(bào)，2018，61（7）：815-824.

[39]?劉玉娣，侯茂林.褐飛虱EST資源的微衛(wèi)星信息分析[J].昆蟲學(xué)報(bào)，2010，53（3）：239-247.

[40]?胡德輝，焦曉真，桑洪玉，等.葉蟬EST資源的微衛(wèi)星信息分析[J].南方農(nóng)業(yè)學(xué)報(bào)，2013，44（4）：558-561.