梁梓強，梁士可，劉婷婷，李廣宏，王方海

(中山大學有害生物控制與資源利用國家重點實驗室//昆蟲學研究所，廣東 廣州510275)

白背飛虱18S核糖體基因克隆及系統(tǒng)發(fā)育*

梁梓強，梁士可，劉婷婷，李廣宏，王方海

(中山大學有害生物控制與資源利用國家重點實驗室//昆蟲學研究所，廣東 廣州510275)

以白背飛虱Sogatellafurcifera長翅雌蟲為材料，提取其基因組DNA，以已知的半翅目昆蟲18S rDNA全序列和白背飛虱18S rDNA部分序列為參考，通過PCR擴增出白背飛虱18S rDNA的未知部分序列，最后通過拼接，首次得到了白背飛虱18S rDNA的全長序列，長度為1 891 bp，堿基組成比例均衡，分別為：24.1% A；23.7% C；27.7% G；24.5% T。將此全長序列與分屬6個目(半翅目、雙翅目、膜翅目、直翅目、鞘翅目、等翅目)昆蟲的18S rDNA序列進行比較，一共發(fā)現(xiàn)有4個保守區(qū)域，其中第二個區(qū)域發(fā)生插入或缺失的情況最少。以18S rDNA第二保守區(qū)域進行系統(tǒng)發(fā)育分析，發(fā)現(xiàn)原同翅目昆蟲多數(shù)均與半翅目昆蟲菜蝽親緣關系較近，但原同翅目飛虱科昆蟲，如白背飛虱、褐飛虱則與半翅目昆蟲菜蝽的親緣關系較遠，表明近年來將原同翅目昆蟲全部劃歸為半翅目還有待商榷之處。

核糖體；18S rDNA；白背飛虱

白背飛虱Sogatellafurcifera(Horvath)屬半翅目Hemiptera飛虱科Delphacidae,是糧食作物水稻的重要害蟲。白背飛虱成蟲常會出現(xiàn)長、短兩種不同的翅型個體，長翅型個體與短翅型個體相比有較強的活動能力，可以進行遠距離遷飛擴散，但是繁殖能力較弱，白背飛虱不同翅型個體的比率是預測其種群動態(tài)變化的關鍵指標之一[1]。

昆蟲基因組存在DNA甲基化現(xiàn)象[2]，同時DNA的甲基化與昆蟲的多型現(xiàn)象密切相關[3]。白背飛虱中，DNA甲基化的研究主要圍繞其對性別分化[4-5]和翅型分化[6]的作用展開了探索。利用甲基化敏感性代表性差異分析法(MS-RDA)分析白背飛虱長、短翅型雌成蟲基因組，發(fā)現(xiàn)兩個特異的DNA甲基化片段分別與18S和28S rDNA高度同源[6]。因此獲得核糖體基因全長序列將有助于進一步分析其甲基化的具體式樣和水平，找出不同翅型個體間存在的差異，從而有可能弄明白核糖體DNA甲基化對白背飛虱翅型分化的具體作用和機制。

真核生物核糖體是細胞內合成蛋白質的場所，由大、小兩個亞基組成，主要成分為核糖體RNA和蛋白質。核糖體RNA具有自我復制和剪切功能。20世紀80年代以來，伴隨著序列分析軟件和計算機的發(fā)展，核糖體基因被廣泛應用于昆蟲的分子系統(tǒng)學研究[7]。

昆蟲分類經(jīng)典的方法主要以形態(tài)學為基礎，相對簡單且快捷，但限制形態(tài)學鑒定的因素較多，如昆蟲生活環(huán)境、分類經(jīng)驗不足等，所以只靠形態(tài)特征進行分類，難以保證絕對的準確性。核糖體基因中18S rDNA具有序列高度保守性，并且許多生物的18S rDNA已被發(fā)表[8-9]，所以比較適合用于各類生物系統(tǒng)發(fā)生等進化關系的研究。

本研究參考了大量半翅目昆蟲的18S rDNA全序列信息，克隆和分析了白背飛虱18S rDNA全長序列，為日后研究白背飛虱核糖體基因全序列甲基化水平和翅型分化的關系奠定基礎。同時，我們將白背飛虱和六個目昆蟲的18S rDNA序列進行比較分析，并根據(jù)18S rDNA的同源保守區(qū)進行了進化樹的描繪，能為分子系統(tǒng)學研究提供參考。

1 材料和方法

1.1 供試昆蟲

白背飛虱成蟲采自廣州市華南農業(yè)大學水稻試驗田。將0.01 g長翅雌成蟲置于液氮中研磨成粉，其基因組DNA的提取采用Universal Genomic DNA Extraction Kit(TaKaRa)，具體方法和步驟參照該試劑盒的說明書。

1.2 引物設計

從GenBank上搜索已經(jīng)發(fā)表的白背飛虱18S rDNA序列(登錄號JF773150，1 733 bp；登錄號JX556779，1 281 bp；登陸號HM017250，1 200 bp)和白背飛虱18S-ITS1-5.8S-ITS2-28S核糖體序列(AB625607，AB625606，AB625605)，拼接成一條長度為1 856 bp的18S rDNA序列。

根據(jù)已經(jīng)公布的半翅目昆蟲的18S rDNA全序列(表1)，利用Megalign軟件分析表1中的8個序列和上述已拼接好的18S rDNA序列，可以得知拼接的白背飛虱18S rDNA不是完整序列，缺失5′端部分。

表1 半翅目昆蟲18S rDNA全序列的登錄號
Table 1 Complete sequence of the Hemiptera 18S rDNA

SpeciesGenBankaccessionNo.PhilaenusspumariusU06480SpissistilusfestinusU06477HenestarisoschaniniAY324853LygushesperusU06476PachygronthaantennataAY324852HemiowoodwardiawilsoniAF131198OkanaganautahensisU06478HackeriellaveitchiAF004766

為了確定白背飛虱18S rDNA的5′端序列，我們分別對柑橘木虱(GCA_000475195.1)和褐飛虱(GCA_000757685.1)的基因組進行分析，找到其18S rDNA及其上游1 000 bp序列，根據(jù)18S rDNA 和上游序列的保守區(qū)域設計了一對引物用于擴增白背飛虱18S rDNA缺失的5′端序列(圖1)，引物具體如下：

18S-F：5′ACGAATCATGAGCTGCTGAC 3′

18S-R：5′ATGGCTTAATCTTTGAGACAAG 3′

18S-F和18S-R預期擴增的序列大小為515 bp。

圖1 PCR擴增白背飛虱18S rDNA 5′序列Fig.1 PCR strategy to amplify 5′sequence of S.furcifera 18S rDNA 灰色部分為白背飛虱18S rDNA上游序列；黑色部分為白背飛虱18S rDNA序列； 矩形部分為序列保守區(qū)域；18S-F和18S-R為擴增白背飛虱18S rDNA 5′端引物

1.3 PCR及檢測

PCR的反應體系為：模板DNA 2 μL，上、下游引物各4 μL，Premix Taq(TaKaRa)25 μL，ddH2O 15 μL。PCR循環(huán)條件為：94 ℃ 30 s； 94 ℃ 20 s，51 ℃ 30 s，72 ℃ 24 s，30個循環(huán)；72 ℃ 8 min，4 ℃ ∞。反應結束后，取5 μL PCR產(chǎn)物進行w=1%凝膠電泳檢測，然后將目的條帶純化回收，克隆并測序。

1.4 18S rDNA序列比對及系統(tǒng)發(fā)育分析

從GenBank中下載了昆蟲綱6個目共9個18S rDNA的序列(表2)，結合白背飛虱18S rDNA全序列，利用DNAstar和MEGA5軟件進行分析。

表2 昆蟲18S rDNA序列來源
Table 2 The source of insect 18S rDNA sequence

OrderSpeciesGenBankaccessionNo.鞘翅目黃粉蟲TenebriomolitorX07801雙翅目果蠅DrosophilamelanogasterNR133559膜翅目蚜繭蜂LipolexisscutellarisAY216699直翅目中華稻蝗OxyachinensisAY037173等翅目澳白蟻MastotermesdarwiniensisAY491152半翅目褐飛虱NilaparvatalugensJF773148半翅目沫蟬PhilaenusspumariusU06480半翅目褐色橘蚜ToxopteracitricidaAY216697半翅目菜蝽EurydemamaracandicaJX997807

2 結果與分析

2.1 白背飛虱18S rDNA基因克隆及序列分析

白背飛虱18S rDNA 5′端序列的PCR擴增結果見圖2。將PCR產(chǎn)物克隆后進行測序，測序完成后，與步驟1.2中已拼接出的白背飛虱18S rDNA序列進行再次拼接。

根據(jù)表1中8個半翅目昆蟲的18S rDNA 基因序列，比較分析18S rDNA基因的起始位點(圖3)，最終確定了白背飛虱18S rDNA全序列，共1 891 bp，其GenBank的登錄號為KU518352。白背飛虱18S rDNA全序列堿基組成分別為：24.1% A；23.7% C；27.7% G；24.5% T，堿基組成比例均衡。

2.2 昆蟲18S rDNA同源性及系統(tǒng)發(fā)育關系分析

通過使用DNAstar中的MegAlign程序對9個

18S rDNA序列進行排序處理，比對發(fā)現(xiàn)4個具有較高同源性的區(qū)域，分別相當于白背飛虱18S rDNA的1～202，335～629, 779～1 431和1 510～1 760 bp的序列。4個同源區(qū)域中，發(fā)生插入或缺失的幾率低，主要為轉換，而在非同源區(qū)域中，發(fā)生插入或缺失的幾率大大增加，如雙翅目或膜翅目昆蟲在此區(qū)域就存在較多的插入序列。

圖2 白背飛虱18S rDNA基因5′端序列擴增結果Fig.2 The PCR products of S.furcifera 18S rDNA 5′sequenceM: DL-2000 DNA Marker ; 1:擴增產(chǎn)物

基于以上的分析，如果根據(jù)18S rDNA全序列進行系統(tǒng)發(fā)育分析，結果就會不太可靠。由于在第二個同源區(qū)域中，出現(xiàn)插入或缺失的現(xiàn)象最少，所以本研究根據(jù)第二個同源區(qū)域進行系統(tǒng)發(fā)育樹的構建，圖4為所構建的系統(tǒng)進化樹。

根據(jù)圖4所示，我們發(fā)現(xiàn)白背飛虱和褐飛虱位于同一分支上，與雙翅目的果蠅親緣關系最接近，與鞘翅目的黃粉蟲最遠。而同屬半翅目的沫蟬、褐色橘蚜、菜蝽3種昆蟲親緣關系接近，但與白背飛虱和褐飛虱的親緣關系較遠。

圖3 半翅目18S rDNA基因起始位點Fig.3 The initiation site of Hemiptera 18S rDNA黑色箭頭指向的核苷酸為半翅目18S rDNA基因起始位點

圖4 基于昆蟲18S rDNA保守區(qū)域的系統(tǒng)發(fā)生樹Fig.4 Phylogenetic analysis of the 18S rDNA conserved domain

3 討 論

核糖體基因一直以來都是人們研究的對象，18S rDNA基因同樣也是白背飛虱研究的熱點之一[6, 10]。過去我們雖然能得到很多有關白背飛虱18S rDNA基因的信息，但是這些基因信息均為部分片段，使得所做的研究具有局限性。通過對白背飛虱雌雄個體和長、短翅個體的MSAP圖譜的分析，發(fā)現(xiàn)DNA甲基化有可能參與了白背飛虱性別和翅型分化的調控[11]，其中不同翅型白背飛虱基因組中存在與18S rDNA同源的差異DNA甲基化片段。本研究首次得到了白背飛虱18S rDNA全序列，為18S rDNA基因的甲基化分析提供可靠的序列參考，并且為今后研究18S rDNA基因甲基化式樣和水平對白背飛虱性別和翅型分化的調控打下基礎。

堿基組合比例對系統(tǒng)發(fā)育的構建有較大的影響[12]，白背飛虱及其它用于構建進化樹的昆蟲的18S rDNA的堿基組合比例均無偏異，故從18S rDNA基因角度可較為真實反映各目昆蟲的進化情況。原屬同翅目的白背飛虱、褐飛虱、沫蟬、褐色橘蚜與半翅目的菜蝽的親緣關系各有不同，沫蟬、褐色橘蚜與菜蝽較接近，而白背飛虱和褐飛虱與菜蝽則較遠。如果從18S rDNA基因進化角度來看，近年來把原同翅目昆蟲全部劃歸為半翅目，似乎還值得商榷，因為飛虱科昆蟲與半翅目昆蟲間親緣關系并不接近。

[1] 周崇高，馮婭琳，陸潮峰，等. 稻飛虱翅型分化的研究進展 [J]. 湖北農業(yè)科學，2014，53(17)：3985-3990. ZHOU C G, FENG Y L, LU C F, et al. Advances in Wing Dimorphism of Rice Planthoppers [J]. Hubei Agricultural Sciences, 2014, 53(17): 3985-3990.

[2] DREWELL R A, LO N, OXLEY P R, et al. Kin conflict in insect societies: a new epigenetic perspective [J]. Trends in Ecology and Evolution, 2012, 27(7): 367-373.

[3] PATALANO S, HORE T A, REIK W. Shifting behavior: epigenetic reprogramming in eusocial insects [J]. Current Opinion in Cell Biology, 2012, 24: 367-373.

[4] ZHANG M, CHEN J L, LIANG S K, et al. Differentially methylated genomic fragments related with sexual dimorophism of rice pests,Sogatellafurcifera[J]. Insect Science, 2015, 22：731-738.

[5] ZHANG M, CHEN J L, ZHOU X S, et al. Different genomic DNA methylation patterns between male and female adults of white-backed planthoppersSogatellafurcifera[J]. Journal of Asia-Pacific Entomology, 2014,17: 917-921.

[6] LIANG S K, LIANG Z Q, ZHOU X S, et al. CpG methylated ribosomal RNA genes in relation to wing polymorphism in the rice pestSogatellafurcifera[J]. Journal of Asia-Pacific Entomology, 2015, 18: 471-475.

[7] 劉殿鋒，蔣國芳. 核基因序列在昆蟲分子系統(tǒng)學上的應用 [J].動物分類學報，2005，30(17)：484-492. LIU D F , JIANG G F. The application of nuclear genes sequences in insect molecular systematic [J]. Acta Zootaxonomica Sinica, 2005, 30(17): 484-492.

[8] NEEB J, van de PEER Y, HENDRICKS L, et al. Compliation of small ribosomal subunit RNA sequences [J]. Nucleic Acids Res, 1990, 18(Suppl): 2237-2317.

[9] KJER K M. Aligned 18S and insect phylogeny [J]. Syst Biol, 2004, 53(3): 506-514.

[10] ZHOU X S, CHEN J L, ZHANG M, et al. Differential DNA Methylation between two wing phenotypes adults of Sogatella furcifera [J]. Genesis, 2013, 51(12): 819-826.

[11] 張梅，陳家林，周曉穗, 等. MSAP在水稻害蟲白背飛虱中的應用研究 [J]. 中山大學學報(自然科學版)，2015, 54(1)：98-102. ZHANG M, CHEN J L, ZHOU X S, et al. The Application of MSAP inSogatellafurcifera(Horvath) [J]. Acta Scientiarum Naturalium Universitatis SunYatseni, 2015, 54(1): 98-102.

[12] LOCKHART P J, HOWE C J, BRYANT D A S, et al. Substitutional bias confounds inference of cyanelle origins from sequence data [J]. J Mol Evol, 1992, 34(2): 153-162.

Clonging and phylogenetic analysis of 18S rDNA gene fromSogatellafurcifera(Horvath)

LIANGZiqiang,LIANGShike,LIUTingting,LIGuanghong,WANGFanghai

(State Key Laboratory for Biocontrol and Institute of Entomology, Sun Yat-sen University, Guangzhou 510275, China)

The long wing female adults ofSogatellafurciferawere selected to extract their genomic DNA. According to the complete 18S rDNA gene sequence from Hemiptera, the unknown 18S rDNA sequence ofS.furciferawas amplified by PCR. The complete 18S rDNA gene sequence ofS.furciferawas obtained for the first time by splicing the unknown and known 18S rDNA sequences, whole sequence length was 1 891 bp. Comparing 18S rDNA sequence fromS.furciferaand those from other insects, four conserved regions were found and the second region had the least insertion or deletion. Based on the second region, phylogenetic analysis was made and the results showed that most Homoptera insects used in this study had a close relatives withEurydemamaracandica(Hemiptera) except Delphacida.

ribosome; 18S rDNA;Sogatellafurcifera

2016-05-16 基金項目：廣東省自然科學基金(S2022020002353,031628)；國家自然科學基金(311771844)

梁梓強(1991年生)，男；研究方向：昆蟲分子生物學；E-mail：616661857@qq.com

王方海(1965年生)，男；研究方向：昆蟲生理生化 ；E-mail：lsswfh@mail.sysu.edu.cn

10.13471/j.cnki.acta.snus.2017.01.019

Q965

A

0529-6579(2017)01-0121-04