草履蚧不同寄主和地理種群遺傳分化的RAPD分析

劉軍俠，劉全超，姜文虎，宋曉英，高寶嘉，*

1河北農業大學林學院，保定　0710002河北省林木種質資源與森林保護重點實驗室，保定　071000

?

草履蚧不同寄主和地理種群遺傳分化的RAPD分析

劉軍俠1，2，劉全超1，姜文虎1，宋曉英1，高寶嘉1，2，*

1河北農業大學林學院，保定071000
2河北省林木種質資源與森林保護重點實驗室，保定071000

摘要:應用隨機擴增多態性DNA(RAPD)技術對草履蚧保定、石家莊、邯鄲16個不同寄主地理種群遺傳多樣性和種群分化進行研究，結果顯示4個RAPD引物共擴增出41個多態性位點，多態位點比率為100％。遺傳距離指數在0.701—0.4360，平均為0.2395。其中以邯鄲楓楊和邯鄲垂絲海棠為寄主的草履蚧種群遺傳距離最小(0.0701);以石家莊紫葉李和邯鄲木槿為寄主的種群遺傳距離最大(0.4360)。遺傳一致度系數在0.6466—0.9290。說明草履蚧不同種群遺傳多樣性豐富并存在遺傳差異。聚類分析結果表明草履蚧種群遺傳多樣性同時受到寄主和地理因素的雙重影響，且不同寄主草履蚧種群已產生明顯的遺傳分化。

關鍵詞:草履蚧;寄主植物;種群分化; RAPD

劉軍俠，劉全超，姜文虎，宋曉英，高寶嘉.草履蚧不同寄主和地理種群遺傳分化的RAPD分析.生態學報，2016，36(2): 298-305.

Liu J X，Liu Q C，Jiang W H，Song X Y，Gao B J.RAPD analysis of genetic differentiation in the different hosts and geographic populations of Drosicha corpulenta.Acta Ecologica Sinica，2016，36(2): 298-305.

隨著分子生物學發展，應用分子系統分析的方法研究昆蟲種群分化和遺傳多樣性，從分子水平探討寄主、環境因子等對昆蟲適應進化的影響，以明確昆蟲適應進化的分子機制。蚧蟲作為一類營相對固定生活的刺吸類害蟲，呈聚集型空間分布，擴散能力弱，其發生與寄主種群以及地理環境有著密切的關系［1-4］，在長期的系統發育中形成了相對穩定的區系組成以適應不同環境，不同的寄主種群、生態地理種群間差異明顯。高寶嘉等對皺大球蚧(Eulecanium kuwana i)和瘤堅大球蚧(E.gigan tean)地理種群的RAPD遺傳分化的研究發現，已發生種內遺傳變異，但尚未達到亞種水平［5］。Burban等通過對松干蚧的研究指出松干蚧種群變異與寄主植物的遺傳結構有關［6］。本研究選取河北省不同草履蚧Drosicha corpulenta(kawana)種群，對其遺傳分化進行RAPD分析，揭示寄主植物、地理環境等生態條件對蚧蟲遺傳分化的影響，為深入探討植物昆蟲協同進化關系和昆蟲適應潛力提供依據。

1材料和方法

1.1供試昆蟲

2011年4—5月于保定、石家莊、邯鄲不同寄主上采集雌成蟲(表1)。不同寄主相距50m以上，每種寄主東、南、西、北方向的上、中、下部枝條分別取樣。將采集的蚧蟲饑餓24h以上，使其消化道內的物質充分消化，以防外源DNA干擾實驗結果;－20℃冰箱中保存。

表1　供試草履蚧材料來源Table 1 Sources of D.corpulenta materials for trial

1.2主要試劑與儀器

本實驗所用的TaqDNA聚合酶、dNTP、蛋白酶K、RAPD引物、DNA標記(DL2000)均由上海生物工程有限公司生產。

1.3研究方法

1.3.1提取DNA

參照Boyce等的方法［7］，并做適當修改。取1只蚧蟲于滅菌的研缽(每個種群取28個樣本)，加入液氮迅速研磨成粉末，放入1.5mL的離心管中，加入700μL預冷的2×CTAB(0.1mol/L Tris2-HCl，1.4mol/L NaCl，0.2mol/L EDTA，0.2％β巰基乙醇，0.05mol/L CT AB，pH8.3)，加入蛋白酶K至終濃度為100μg/mL，65℃下水浴1h，用12000r/min離心5min，取上清液，加入等體積得酚∶氯仿∶異戊醇(25∶24∶1)輕輕混勻，反復抽提兩次。再經等體積的異丙醇沉淀1h，4℃下10000r/min離心10min，去上清，70％的乙醇清洗1—2次，干燥后加入適量TE溶液溶解DNA，加入RNaesA至終濃度為2μg/L后，37℃下水浴1h，－20℃保存待用。

1.3.2 RAPD擴增及產物檢測

經預試驗篩選出4條可用引物，其引物名稱及序列見表2。參照Williams的方法［8］，反應體系25μL，其中模板50ng，10×Buffer2.5μL，Mg2+114—210 mmol/L(因座位而異)，dNTP200μmol/L，引物(0.8μmol/L)，Taq酶1.5U反應條件為: 94℃預變性1min; 92℃變性1 min，37℃退火1 min(因座位而異)，72℃延伸2 min，45個循環; 72℃再延伸10min; 4℃終止反應。

擴增產物用1.5％瓊脂糖凝膠(含1％核酸顯色劑)分離，電泳緩沖液為0.5×TBE，上樣量為10μL。電泳結束后，用凝膠成像儀拍照。

1.3.3數據分析

RAPD為顯性標記，同一引物的擴增產物在瓊脂糖凝膠電泳中遷移律相同的條帶被認為是同源性的，屬于同一位點的條帶按清晰可見的強帶或反復出現的弱帶記為1，無帶記為0，形成二元數據。應用Popstream32 和NTSYSpc2.1軟件，計算種群遺傳多樣性參數并根據遺傳相似系數進行UPGMA種群聚類分析。

表2　隨機引物編號及序列Table 2 RAPD primer codes and their sequences

2結果與分析

2.1遺傳多態位點分析

引物擴增結果見圖1，不同種群遺傳多樣性統計見表3。供試引物對16個種群448個個體進行擴增，共產生41條清晰、穩定的片段，平均每條引物擴增10.25條片段，各片斷分子量大小在200—2000bp之間，種群內和種群間并未產生各自特征性帶譜。實驗共檢測出41個多態位點，多態位點比率為100％，說明草履蚧種群遺傳多樣性豐富。但不同種群的擴增產物多態性結果不同，各種群多態位點比率在31.71％—63.41％之間，其中邯鄲紫薇寄主種群多態位點比率最低(31.71％)，邯鄲紫葉李寄主種群(63.41％)和邯鄲桑樹寄主種群(60.98％)多態位點最高，遺傳多樣性較其它種群豐富。

圖1　引物H9對16個草履蚧種群的RAPD-PCR結果Fig.1　RAPD-PCR results of primer H9 for 16 populations of D.corpulenta1—4: HL、HF、HZ、HS; 5—8: HB、HW、HC、HM; 9—12: SC、SX、SZ、SB; 13—16: BY; BB; BZ; BBT; 17: DL2000

表3　不同種群遺傳多樣性統計Table 3 Statistics on the genetic diversity of different populations

2.2遺傳多樣性分析

Nei指數估計的不同種群間遺傳距離和遺傳一致度結果見表4。由表4可知，16個種群Shannon信息指數的變化范圍在0.1708—0.2967之間，平均值為0.2306。邯鄲桑樹種群的Shannon信息指數最高，邯鄲紫薇種群最低。16個種群Nei's基因多樣性指數變化范圍在0.1101—0.1991之間，平均值為0.1511。邯鄲桑樹種群的基因多樣性最高，保定紫葉碧桃種群最低，這與Shannon信息指數估計的遺傳多樣性稍有差異。說明同一物種不同種群遺傳多樣性存在差異。

2.3遺傳相似度與遺傳距離分析

由表4可知，16個草履蚧種群的遺傳相似度在0.6466—0.9290之間，平均值為0.7873;遺傳距離在0.0701—0.4360之間，平均值為0.2395，說明不同種群間遺傳差異較大。其中以邯鄲楓楊和以邯鄲垂絲海棠為寄主的草履蚧種群最相似，遺傳距離最小(0.0701);以石家莊紫葉李和以邯鄲木槿為寄主的種群遺傳距離最大(0.4360)，說明二者的遺傳差異最大。

紫葉李邯鄲種群與石家莊種群遺傳距離為0.3330，與保定種群遺傳距離為0.3426;紫葉碧桃邯鄲種群與石家莊種群遺傳距離為0.3170，與保定種群遺傳距離為0.3816;即同種寄主不同地理種群隨地理距離的增加遺傳距離增加。

2.4種群遺傳相似度聚類分析

用NTSYSpc2.1對草履蚧種群的遺傳相似系數進行UPGMA聚類分析，進一步分析種群間遺傳關系(圖2)。結果顯示:同一地理寄主種群絕大多數個體都能聚為一類，說明同一種群中個體間的遺傳特性相似程度高，遺傳穩定性;草履蚧相同寄主種群首先聚類，遺傳相似系數在0.67水平，說明寄主因素對草履蚧的種群分化產生明顯的影響;保定、石家莊、邯鄲的紫葉碧桃種群各自聚類，且遺傳距離較遠;同樣以上3地的紫葉李種群也各自聚類，說明地理因素對種群的分化產生影響;草履蚧不同寄主的三大地理種群并未產生明顯地理聚類，說明相對于地理因素，寄主因素對草履蚧種群分化的影響更加明顯。同一地區不同寄主種群聚類見圖3、圖4、圖5，可見同一地區同一寄主的個體聚為一類，而同一地區不同寄主間聚類規律不明顯。

3結論與討論

本文通過以上對草履蚧不同寄主和地理種群遺傳分化的RAPD分析，表明: Shannon信息指數和Nei's基因多樣性指數顯示草履蚧不同地理種群以及不同寄主種群遺傳多樣性豐富程度不同，草履蚧不同寄主種群之間存在不同程度的遺傳差異，并已產生種群分化。草履蚧相同寄主不同地理種群之間存在遺傳差異，地理因素對草履蚧種群分化產生影響。聚類分析結果顯示，草履蚧種群分化同時受到寄主因素和地理因素的雙重影響，相對于地理因素，寄主對草履蚧種群分化的影響更加明顯。

草履蚧寄主植物廣泛，主動擴散距離有限，因此其發生與寄主植物關系十分密切。本研究表明草履蚧不同寄主種群之間已經產生了明顯的種群分化，這與草履蚧寄主生態條件、環境脅迫有密切關系。寄主植物自身的物理、化學等因素可影響植食性昆蟲的生物學習性和遺傳特性，促進昆蟲種群分化，而昆蟲對寄主植物具有選擇性，對于不同寄主植物喜好程度也不盡相同［9-10］。昆蟲學家認為植食性昆蟲可因嗜食的寄主不同而造成遺傳隔離，由此形成不同的物種。Takada、Blackman、Weber等研究表明煙蚜在不同寄主植物上存在著種群分化［11-13］。本實驗中寄主植物分類科屬不盡相同，其生物學特性和化學成分各異，我們推測寄主植物種類不同是造成草履蚧種群分化的重要原因。

此外，大氣污染以及昆蟲群落結構［14-15］等環境脅迫均能影響寄主植物與昆蟲的關系，本實驗寄主植物除受到草履蚧危害外，同時受到其他不同種昆蟲取食。不同寄主植物上昆蟲的種類、數量以及植物所處環境的不同可能都對草履蚧種群分化期作用，下一步有待從多因素分析探討草履蚧種群分化原因。

圖2　16個草履蚧種群聚類分析圖Fig.2 The UPGMA dendrogram cluster of 16 Drosicha corpul populations種群代號同表1(See table 1 for explanation of population code);為了方便比較查看，本圖中每個種群采用5個樣本的聚類分析;其聚類分析趨勢與全部個體聚類圖相一致，具有代表性

圖3　石家莊草履蚧種群聚類分析圖Fig.3 The UPGMA dendrogram cluster of Drosicha corpul populations of Shijiazhuang

圖4　保定草履蚧種群聚類分析圖Fig.4 The UPGMA dendrogram cluster of Drosicha corpul populations of Baoding

圖5　邯鄲草履蚧種群聚類分析圖Fig.5 The UPGMA dendrogram cluster of Drosicha corpul populations of Handan

參考文獻(References):

［1］梁隱泉，高寶嘉，甄志先，王江柱，劉志群，牛建忠.棗園昆蟲群落及其與日本龜蠟蚧發生的關系.應用生態學，2006，17(3): 472-476.

［2］薛皎亮，謝映平，范三紅，李艷芳.城市污染促進瘤堅大球蚧種群增長的營養學研究.林業科學.2004，40(1): 111-116.

［3］崔巍，高寶嘉.華北經濟樹種主要蚧蟲及其防治.北京:林業出版社，1995.

［4］Beuning L L，Murphy P，Wu E，Batchelor T A，Morris B A M.Molecular-based approach to the differentiation of mealybug(Hemiptera: Pseudococcidae)species.Journal of Economic Entomology，1999，92(2): 463-472.

［5］高寶嘉，李夢，王春和，高立杰.皺大球蚧(Eulecanium kuwanai)和瘤堅大球蚧(Eulecanium gigantean)地理種群的RAPD遺傳分化.生態學報，2007，27(3): 918-923.

［6］Burban C，Petit R J，Carcreff E，Jactel H.Rangewide variation of the maritime pine bast scale Matsucoccus feytaudi Duc.(Homoptera: Matsucoccidae)in relation to the genetic structure of its host.Molecular Ecology，1999，10(8): 1593-1602.

［7］Boyce T M，Zwick M E，Aquado C F.Mitochondrial DNA in the bark weevils: size，structure and heteroplasmy.Genetics，1989，123(4): 825-836.

［8］Williams J G，Kubelik A R，Livak K J，Rafalski J A，Tinqey S V.DNA polymorphisms amplified by arbitrary primers are useful as genetic markers.Nucleic Acids Research，1990，18(22): 6531-6535.

［9］Kelley S T，Farrell B D，Mitton J B.Effects of specialization on genetic differentiation in sister species of bark beetles.Heredity，2000，84(2): 218-227.

［10］Shivayogeshwara B，Mallikharjunaiah H，Krishnaprasad N K，Shetty M V N.Integrated management of spodoptera litura Fabricius(Noctuidae: Lepidoptera)in FCV tobacco crop.Tobacco Research，1991，17(2): 59-61.

［11］Takada H.Differences in adaptation to tobacco between host forms of Myzus persicae(Sulzer)(Homoptera: Aphididae).Japanese Journal of Applied Entomology and Zoology，1988，32(2): 115-119.

［12］Blackman R L.Morphological discrimination of a tobacco-feeding form from Myzus persicae(Sulzer)(Homoptera: Aphididae)and a key to new world Myzus(nectarosiphon)species.Bulletin of Entomological Research，1987，77: 713-730.

［13］Weber G.Genetic variability in host plant adaptation of the green peach aphid，Myzus persicae.Entomologia Experimentalis et Applicata，1985，38(1): 49-56.

［14］Nelson B W，John M.Lead residues in Eastern Tent caterpillars(Malacosoma americanum)and their host plant(Prunus serotina)close to a major highway.Environment Entomology，1980，9(1): 10-12.

［15］Alstad D N，Edmunds G F Jr，Weinstein L H.Effects of air pollutants on insect populations.Annual Review of Entomology，1982，27: 369-384.

RAPD analysis of genetic differentiation in the different hosts and geographic populations of Drosicha corpulenta

LIU Junxia1，2，LIU Quanchao1，JIANG Wenhu1，SONG Xiaoying1，GAO Baojia1，2，*
1 College of Forestry，Agricultural University of Hebei，Baoding 071000，China
2 Key Laboratory of Germplasm Resources of Forest Trees and Forest Protection College of Forest of Hebei，Baoding 071000，China

Abstract:The genetic diversity and population differentiation among geographic populations of 16 different Drosicha corpulenta host collected from Baoding，Shijiazhuang，and Handan in Hebei were analyzed by RAPD.The results showed that four RAPD primers were amplified out of 41 polymorphism loci，and that the percent of polymorphic loci was 100％.The genetic distance index of the populations ranged from 0.0701 to 0.4360 and the average value was 0.2395.The genetic distance index for D.corpulenta populations on Pterocarya stenoptera and Malus halliana(Voss.)Koehne in Handan was the smallest(0.0701).The genetic distance index for D.corpulenta populations on Prunus cerasifera in Shijiazhuang and Hibiscus syriacus in Handan was the largest of all，and the similarity coefficients ranged from 0.6466 to 0.9290.The genetic diversity of the different D.corpulenta populations was classified as abundant，and genetic differences were found.According to the results of cluster analysis，the genetic diversity of D.corpulenta populations was affected by both host and geographic location，and the different host populations of D.corpulenta showed obvious genetic differentiation.

Key Words:Drosicha corpulenta; host plant; population differentiation; RAPD

*通訊作者

Corresponding author.E-mail: baojiagao@ 163.com

收稿日期:2013-06-13;網絡出版日期: 2015-06-10

基金項目:河北省自然基金(C2009000592)

DOI:10.5846/stxb201306131697