台湾蠛蠓rDNA-ITS序列测定与分析
陈敏1, 肖颖2, 林妙端3, 黄恩炯1
1.福建国际旅行卫生保健中心,福州 350001
2.宁德出入境检验检疫局,宁德 352100
3.三明出入境检验检疫局综合技术服务中心,三明 365000
通讯作者:黄恩炯,Email:enjiong@163.com
摘要

目的 探讨rDNA-ITS基因序列在吸血蠓分子鉴定中的应用前景。方法 通过PCR扩增获取台湾蠛蠓rDNA-ITS基因,测定分析ITS1和ITS2基因序列,构建系统发育树。结果 台湾蠛蠓ITS1和ITS2片段长度分别为310 bp与360 bp,与汤斯维尔铗蠓的同源性最大,分别为71%和92%。系统发育分析显示,基于ITS1、ITS2基因序列的分子鉴定结果与传统形态学鉴定结果一致。结论 rDNA-ITS基因序列可作为蠓科昆虫分子鉴定的分子标记。

关键词: 台湾蠛蠓; rDNA; ITS1; ITS2
中图分类号:R384.5 文献标志码:A 文章编号:1002-2694(2018)08-0728-05
Determination and analysis of rDNA-ITS genes from Lasiohelea taiwana
CHEN Min1, XIAO Ying2, LIN Miao-duan3, HUANG En-jiong1
1.Fujian International Travel Healthcare Center, Fuzhou 350001, China
2.Ningde Entry-Exit Inspection and Quarantine Bureau, Ningde 352100, China
3.Comprehensive Technology Service Center of Sanming Entry-Exit Inspection and Quarantine Bureau, Sanming 365000, China
Corresponding author: Huang En-jiong, Email: enjiong@163.com
Abstract

The ITS1 and ITS2 genes sequences of Lasiohelea taiwana were amplified and determined to investigate the feasibility of rDNA-ITS in the molecular identification of blood-sucking midges. PCR was employed in the present study and phylogenetic tree was analyzed according to the ITS1 and ITS2 genes sequences. It was showed that the framents of ITS1 and ITS2 are 310 bp and 360 bp, respectively, sharing the greatest homology of 71% and 92% with Forcipomyia townsvillensis. It was demonstrated in the phylogenetic analysis that the molecular identification based on ITS1 and ITS2 gene sequences is consistent with the morphological identification. Therefore, the sequences of rDNA-ITS can be the molecular marker in the molecular identification of biting midges.

Keyword: Lasiohelea taiwana; rDNA; ITS1; ITS2

台湾蠛蠓(Lasiohelea taiwana)又名台湾铗蠓(Forcipomyia taiwana)、南方蠛蠓, 是闽台等地重要的吸血蠓。该蠓在南方诸省均有分布, 1957年首次从福建乙脑流行区捕获的台湾蠛蠓中分离出乙型脑炎病毒[1], 而后中山医学院又从广州采获的台湾蠛蠓中又分离出乙型脑炎病毒[2]。闽、粤两地的调查均认为该蠓的数量消长曲线与乙型脑炎病毒的流行季节吻合, 以及本种蠛蠓嗜吸人血, 又刺叮乙脑病毒的宿主动物-猪等原因, 因而认为该蠛蠓可能是当地乙型脑炎的媒介昆虫之一, 但至今未得到证实。

目前, 台湾蠛蠓的鉴定主要依靠传统的形态学鉴定, 需要进行玻片标本制作, 步骤多、周期长、耗时费力, 因此, 亟需寻求一种准确而新的分类方法。rDNA广泛分布于各种生物细胞中, 在功能上具有高度的保守性和良好的时钟性, 其测定序列常用于物种分类和生物间系统发育的研究, 尤其是物种的种间分类研究[3, 4, 5]。目前rDNA作为遗传标记在寄生虫鉴定中的应用日趋广泛, 前人对此作了详尽的描述[6, 7, 8]。ITS是位于rDNA上18S和28S基因之间的区域片段, 包括ITS1和ITS2两段序列, 由于该区域不加入成熟核糖体, 所以受到的选择压力较小, 进化速度较其他区域快, 具有种内变异小而种间变异大的特性, 加上协同进化作用保证了该片段在基因组不同单元间的一致性, 因而适合进行各种分子操作, 是较理想的种间鉴定和系统发育分析的遗传标志[9, 10, 11, 12]。近几年来国内外对ITS序列标记在蠓科昆虫中的研究和应用越来越多, rDNA的ITS1和ITS2序列在蠓类昆虫的鉴定中发挥了重要的作用[13, 14, 15, 16, 17, 18, 19, 20, 21, 22]。本研究以台湾蠛蠓为对象, 通过PCR扩增其ITS基因, 并对ITS基因序列进行分析, 以期为建立吸血蠓的分子鉴定技术奠定基础。

1 材料与方法
1.1 试验对象

台湾蠛蠓采自福州周边地区, 试虫均为雌虫, 标本均用95%乙醇浸泡保存。取出保存的标本, 将头、生殖器部分和翅切下, 用于常规形态鉴定, 剩余部分移入Eppendorf管中-20 ℃保存。每份标本组织单管保存并编号, 编号与用于形态学鉴定的标本编号一致, 确保标本组织的唯一性。形态学鉴定为台湾蠛蠓的标本组织用于后续DNA提取。同个地点采集的标本每次20~30只。

1.2 基因组DNA提取

参照昆虫基因组DNA提取试剂盒说明书提取台湾蠛蠓DNA, 略有改动。

1.3 ITS基因PCR扩增

1.3.1 ITS1基因 扩增引物参考文献[15]:Pan CulF 5'-GTAGGTGAACCTGCGGAAGG-3', Pan CulR 5'-TGCGGTCTTCATCGAACCCAT-3'[15]。反应体系和反应条件见表1表2

表1 ITS1基因PCR反应体系 Tab.1 PCR reaction system of ITS1 genes
表2 ITS1基因PCR反应条件 Tab.2 PCR reaction conditions of ITS1 genes

1.3.2 ITS2基因 扩增引物参考文献[20]:5.8S 5'-GATGAAGACCGCAGCAAACT-3', 28S 5'-ATTTGGGGGTAGTCACACAT-3'。 PCR反应体系同ITS1, PCR反应条件见表3

表3 ITS1基因PCR反应条件 Tab.3 PCR reaction conditions of ITS1 genes
1.4 电泳检测

取5 μ L PCR产物, 在浓度为1.2%琼脂糖凝胶中进行电泳, 电泳缓冲液为0.5× TBE。电泳结束后于凝胶成像仪中观察、拍照并记录结果。

1.5 序列测定与分析

PCR产物的纯化与测序委托上海生工生物工程公司完成。用DNAMAN、BioEdit、clustalX、MEGA等软件进行序列分析并构建系统进化树。

2 结 果
2.1 rDNA-ITS基因的扩增

台湾蠛蠓rDNA-ITS1、ITS2基因PCR扩增产物经琼脂糖凝胶电泳后, 与DL-1000 Marker 进行比较, 发现ITS1和 ITS2分别在约310 bp、约360 bp左右处有一条特异性的明亮条带(图1)。

图1 台湾蠛蠓rDNA-ITS基因扩增结果
M: DL-1000 Marker; 泳道1-4: ITS1基因; 泳道5-8: ITS2基因
Fig.1 PCR amplification of rDNA-ITS gene of Lasiohelea taiwana

2.2 台湾蠛蠓ITS1、ITS2序列分析

将所测得的序列, 通过NCBI BLAST在线分析系统与GenBank收录的全部序列进行比较, 结果表明:台湾蠛蠓ITS1株间同源性96%~98%, 与汤斯维尔铗蠓Forcipomyia townsvillensis的同源性最大, 有172个碱基相同, 同源性为71%; 其ITS2片段株间同源性99%~100%, 同样是与汤斯维尔铗蠓的同源性最大, 有144个相同碱基, 同源性为92%。

比较ITS1序列时, 除了本试验所测台湾蠛蠓序列外, 还选用了GenBank所查得的几种蠓相应的ITS1序列。8种蠓的ITS1序列长度范围为 271 ~469 bp, G+C含量为26.79%~34.53%。台湾蠛蠓的ITS1序列长度与两铗蠓的序列长度相当, 均在270 ~300 bp之间, 而库蠓的ITS1序列长度在不同的蠓种间相差很大。台湾蠛蠓的ITS1序列长度与G+C含量均介于库蠓与铗蠓之间(表4)。

表4 8种蠓ITS1基因序列组成 Tab.4 The sequences composition of ITS1 genes from 8 biting midges

比较ITS2序列时, 将台湾蠛蠓ITS2序列与GenBank所查得的3种库蠓、2种铗蠓以及双翅目中的白纹伊蚊(Aedes albopictus)、有羽摇蚊(Chironomus plumosus)的已知序列进行比较。8种双翅目昆虫的ITS2序列长度范围为 205 ~435 bp, G+C含量为20.87%~56.35%。台湾蠛蠓的ITS2序列长度与两铗蠓的序列长度相当, 均在300 bp左右。库蠓的ITS2序列长度较短, 范围在200 ~250 bp, 而白纹伊蚊、有羽摇蚊的ITS2序列长度明显大于蠓类, 分别为424 bp与435 bp(表5)。

表5 8种双翅目昆虫ITS2的登录号及序列组成 Tab.5 The sequences composition of ITS2 genes from 8 dipterons
2.3 台湾蠛蠓ITS1、ITS2的系统发育分析

采用最大似然法(Maximum Likelihood)构建台湾蠛蠓的ITS1、ITS2序列系统进化树(图2, 图3)。结果显示, 基于ITS1、ITS2基因分类结果与形态学传统分类结果一致, 在系统进化树上聚类成库蠓属、蠛蠓属和铗蠓属3个分支。在系统进化树上各蠓种均先与同属蠓种聚类, 再与不同属其他蠓种聚类。

图2 台湾蠛蠓rDNA-ITS1基因核苷酸序列系统进化树Fig.2 Phylogenetic tree of rDNA-ITS1 gene sequences of Lasiohelea taiwana

图3 台湾蠛蠓rDNA-ITS2基因核苷酸序列系统进化树Fig.3 Phylogenetic tree of rDNA-ITS2 gene sequences of Lasiohelea taiwana

3 讨 论

目前, 对台湾蠛蠓的识别还是采用传统的形态学鉴定方法, 但是在病媒生物监测中, 很有可能得到的是不完整的标本。在这种情况下, 依靠传统的形态学鉴定的办法难以迅速得到结果, 而分子鉴定技术能很好地解决这一难题。

rDNA是生物体内最保守的基因之一, 其内转录间隔区ITS具有极大的种间甚至种内变异性, 被认为是发展敏感、特异检测系统的理想靶物。Mathieu等[16]利用rDNA-ITS1作为靶标基因建立了不显库蠓复合体(obsoletus complex)分子鉴定方法; Perrin等[15]对残肢库蠓(C. imicola)进行ITS1序列测定和系统发育分析, 种内同源性高达99.8%; 对蠓ITS2的研究也都表明, ITS2作为靶标基因不仅有利于鉴定蠓新种, 还能有效辅助澄清库蠓亚属里的同种异名[20, 21], 分子鉴定与形态学鉴定相佐证能够解决蠓形态学鉴定中的难题。

本研究中, 从ITS1和ITS2的片段长度上看, 台湾蠛蠓的片段长度均较接近铗蠓, 与库蠓及其他双翅目昆虫如白纹伊蚊、有羽摇蚊的ITS1和ITS2长度区分明显。蠛蠓属原是铗蠓属内的一个亚属, 虞以新[23]根据其特有的形态将其独立为一个属。本研究所得台湾蠛蠓的ITS1、ITS2进行系统发育树显示台湾蠛蠓和铗蠓属先聚类于同一分支, 从分子生物学方面验证了形态分类学结果。

而库蠓ITS序列长度在不同的蠓种间相差很大, 在ITS2系统发育树上, 变翅库蠓(C. variipennis)与其他库蠓不聚类于同一分支, 这与国外学者利用COⅡ 基因研究蠓科分子分类时, 发现C. saundersi不与其他库蠓聚类于同一分支的现象一致[24], 可能是由于库蠓种类众多, 种间基因变异大的原因所致。这也验证了真核生物ITS序列具有高度变异性的特点, 在物种进化和遗传关系研究中具有重要的价值。

总之, rDNA-ITS可有效应用于蠓类的分子鉴定, 而且能够提供足够的遗传信息, 但在进行分类和种属鉴定时, 不能仅考虑片段的长度, 应结合序列组成和同源性进行综合评价。

The authors have declared that no competing interests exist.

参考文献
[1] 吴皎如, 吴树吟. 由糠蚊科蠛蠓属 Lasiohelea台湾蠛蠓分离出乙型脑炎病毒[J]. 微生物学报, 1957, 5(1) : 22-26. 10.13343/j.cnki.wsxb.1957.01.010 [本文引用:1]
[2] 容瓘, 柯小麟. 从广州市天然界蠛蠓分离流行性乙型脑炎病毒的试验[J]. 微生物学报, 1962, 8(3) : 280-284. [本文引用:1]
[3] Dover G. Molecular drive: a cohesive mode of species evolution[J]. Nature, 1982, 299(5879) : 111-117. 10.1017/S1464793104006657 [本文引用:1]
[4] Hillis DM, Davis SK. Ribosomal DNA: intraspecific polymorphism, concerted evolution and phylogeny reconstruction[J]. Syst Zool, 1988, 37(1) : 63-66. DOI: 10.2307/2413191 [本文引用:1]
[5] Hillis DM, Dixon MT. Ribosomal DNA: molecular evolution and phylogenetic inference[J]. Q Rev Biol, 1991, 66(4) : 411-453. DOI: 10.1086/417338 [本文引用:1]
[6] Gajahar AA, Marquardt WC, Hall R, et al. Ribosomal RNA sequence of Sarcocystis muris, Theileria annulata and Crypthecodinium cohnii reveal evolutionary relationship among apicomlexan, dinoflagellates and ciliates[J]. Mol Biochem Parasitol, 1991, 45(1) : 147-154 10.1016/0166-6851(91)90036-6 [本文引用:1]
[7] Ellis J, Hefford C, Baverstock PR, et al. Ribosomal DNA sequence comparison of Babesia and Theileria[J]. Mol Biochem Parasitol, 1992, 54(1) : 97-96. 10.1016/0166-6851(92)90097-4 [本文引用:1]
[8] Allsopp MT, Cavalier-Smith T, De Waal DT, et al. Phylogeny and evolution of the piroplasms[J]. Parasitology, 1994, 108: 147-152. DOI: 10.1017/S0031182000068232 [本文引用:1]
[9] Paskewitz SM, Wesson DM, Collins FH. The internal transcribed spacers of ribosomal DNA in five members of the Anopheles gambiae species complex[J]. Insect Mol Biol, 1993, 2(4) : 247-257. DOI: 10.1111/j.1365-2583.1994.tb00144.x [本文引用:1]
[10] Kuhls K, Mauricio IL, Pratlong F, et al. Analysis of ribosomal DNA internal transcribed spacer sequences of the Leishmania donovani complex[J]. Microbes Infect, 2005, 7(11-12) : 1224-1234. DOI: 10.1016/j.micinf.2005.04.009 [本文引用:1]
[11] 牛庆丽, 罗建勋, 殷宏. 转录间隔区(ITS)在寄生虫分子生物学分类中的应用及其进展[J]. 中国兽医寄生虫病, 2008, 16(4): 41-47. DOI: 10.3969/j.issn.1674-6422.2008.04.009 [本文引用:1]
[12] 李晓娟, 杨毅梅. 核糖体DNA内转录间隔区序列标记在寄生虫学中的应用[J]. 大理学院学报, 2008, 7(4) : 79-81. DOI: 10.3969/j.issn.1672-2345.2008.04.029 [本文引用:1]
[13] Cetre-Sossah C, Baldet T, Delecolle JC, et al. Molecular detection of Culicoides spp. and Culicoides imicola, the principal vector of bluetongue (BT) and African horse sickness (AHS) in Africa and Europe[J]. Vet Res, 2004, 35(3) : 325-337. DOI: 10.1051/vetres:2004015 [本文引用:1]
[14] Cetre-Sossah C, Mathieu B, Setier-Rio ML, et al. Development and evaluation of a real-time quantitative PCR assay for Culicoides imicola, one of the main vectors of bluetongue (BT) and African horse sickness (AHS) in Africa and Europe[J]. Res Vet Sci, 2008, 85(2) : 372-382. DOI: 10.1016/j.rvsc.2007.12.001 [本文引用:1]
[15] Perrin A, Cetre-Sossah C, Mathieu B, et al. Phylogenetic analysis of Culicoides species from France based on nuclear ITS1-rDNA sequences[J]. Med Vet Entomol, 2006, 20(2) : 219-228. DOI: 10.1111/j.1365-2915.2006.00616.x [本文引用:3]
[16] Mathieu B, Perrin A, Baldet T, et al. Molecular identification of western European species of Obsoletus complex (Diptera: Ceratopogonidae) by an internal transcribed spacer-1 rDNA multiplex polymerase chain reaction assay[J]. J Med Entomol, 2007, 44(6) : 1019-1025. DOI: 10.1093/jmedent/44.6.1019 [本文引用:2]
[17] Stephan A, Clausen PH, Bauer B, et al. PCR identification of Culicoides dewulfi midges (Diptera: Ceratopogonidae), potential vectors of bluetongue in Germany[J]. Parasitol Res, 2009, 105(2) : 367-371. DOI: 10.1007/s00436-009-1407-z [本文引用:1]
[18] Matsumoto Y, Yanase T, Tsuda T, et al. Charceterization of internal transcribed spacer (ITS1)-ITS2 region of ribosomal RNA gene from 25 species of Culicoides biting midges (Diptera: Ceratopogonidae) in Japan[J]. J Med Entomol, 2009, 46(5) : 1099-1108. DOI: 10.1603/033.046.0517 [本文引用:1]
[19] Li GQ, Hu YL, Kanu S, et al. PCR amplification and sequencing of ITS1 rDNA of Culicoides arakawae[J]. Vet Parasitol, 2003, 112(1/2) : 101-108. DOI: 10.1016/S0304-4017(02)00409-0 [本文引用:1]
[20] Gomulski LM, Meiswinkel R, Delecolle JC, et al. Phylogeny of the subgenus Culicoides and related species in Italy, inferred from internal transcribed spacer 2 ribosomal DNA sequences[J]. Med Vet Entomo, 2006, 20(2) : 229-238. DOI: 10.1111/j.1365-2915.2006.00620.x [本文引用:3]
[21] Monaco F, Benedetto L, Marcello VD, et al. Development and preliminary evaluation of a real-time polymerase chain reaction for the identification of Cuilcoides obsoletus sensu strictu, C. scoticus and C. montanus in the Obsoletus Complex in Italy[J]. Vet Ital, 2010, 46(2) : 215-220. [本文引用:2]
[22] 胡友兰, 李国清. 荒川库蠓rDNA ITS2的序列测定与分析[J]. 动物学报, 2003, 49(2) : 277-280. DOI: 10.3969/j.issn.1674-5507.2003.02.019 [本文引用:1]
[23] 虞以新, 刘康南. 中国蠛蠓的研究[M]. 北京: 科学出版社, 1982: 1-84. [本文引用:1]
[24] Bechenbach AT, Borkent A. Molecular analysis of the biting midges (Diptera: Ceratopogonidae), based on mitochondrial cytochrome oxidase subunit 2[J]. Mol Phylogenet Evol, 2003, 27: 21-35. DOI: 10.1016/S1055-7903(02)00395-0 [本文引用:1]