中国媒介生物学及控制杂志  2015, Vol. 26 Issue (6): 550-554

文章信息

刘飞, 郭秀霞, 程鹏, 王海防, 刘丽娟, 杨培培, 赵玉强, 王怀位, 公茂庆
LIU Fei, GUO Xiu-xia, CHENG Peng, WANG Hai-fang, LIU Li-juan, YANG Pei-pei, ZHAO Yu-qiang, WANG Huai-wei, GONG Mao-qing
山东省淡色库蚊核糖体DNA第二内转录间隔区基因多态性分析
Genetic polymorphism analysis on rNDA-ITS2 from different geographical strains of Culex pipiens pallens in Shandong province
中国媒介生物学及控制杂志, 2015, 26(6): 550-554
Chin J Vector Biol & Control, 2015, 26(6): 550-554
10.11853/j.issn.1003.4692.2015.06.003

文章历史

收稿日期: 2015-08-21
网络出版时间: 2015-10-14 14:39
山东省淡色库蚊核糖体DNA第二内转录间隔区基因多态性分析
刘飞1,2,3, 郭秀霞1, 程鹏1, 王海防1, 刘丽娟1, 杨培培1, 赵玉强1, 王怀位1, 公茂庆1    
1 山东省医学科学院, 山东省寄生虫病防治研究所, 山东省医药卫生媒介生物学重点实验室, 山东 济宁 272033;
2 济南大学山东省医学科学院医学与生命科学学院, 济南 250012;
3 济宁市第一人民医院, 山东 济宁 272033
摘要: 目的 探讨山东省不同地区淡色库蚊核糖体DNA第二内转录间隔区(rDNA-ITS2)基因序列特征,分析其遗传变异。方法 在山东省的济南、济宁、青岛、聊城、临沂、淄博市以及江苏省南京市采集淡色库蚊成蚊,并取实验室选育的抗氯氰菊酯与抗敌敌畏品系成蚊,提取蚊虫DNA,PCR扩增淡色库蚊的rDNA-ITS2,PCR产物纯化后连接至pMD19-T载体进行基因测序,生物信息学分析rDNA-ITS2 基因多态性。结果 PCR特异扩增山东省6个不同地理株及南京市和实验室2个不同品系的淡色库蚊rDNA-ITS2区域,获得460 bp左右的扩增片段。测序结果显示淡色库蚊rDNA-ITS2出现36个变异位点,表现出基因的多态性,其rDNA-ITS2基因单核苷酸存在转换、颠换和缺失。结论 山东省不同地理区域和不同抗性的淡色库蚊在rDNA-ITS2区出现遗传多态性,可能与蚊虫生态环境有一定的关系。
关键词: 淡色库蚊     核糖体DNA第二内转录间隔区     基因多态性    
Genetic polymorphism analysis on rNDA-ITS2 from different geographical strains of Culex pipiens pallens in Shandong province
LIU Fei1,2,3, GUO Xiu-xia1, CHENG Peng1, WANG Hai-fang1, LIU Li-juan1, YANG Pei-pei1, ZHAO Yu-qiang1, WANG Huai-wei1, GONG Mao-qing1     
1 Shandong Institute of Parasitic Diseases, Shandong Academy of Medical Sciences, Department of Medical Entomology, Shandong Province Key Laboratory of Health and Vector Biology, Jining 272033, Shandong Province, China;
2 School of Medicine and Life Sciences, University of Jinan, Shandong Academy of Medical Sciences, Jinan 250012, Shandong Province, China;
3 Shandong Jining No. 1 People's Hospital, Jining 272033, Shandong Province, China
Abstract: Objective To explore gene sequence characteristics of rNDA-ITS2 in different populations of Culex pipiens pallens from Shandong province and analyze the genetic variation of these mosquitoes. Methods Adult-female mosquitoes were collected from different cities: Jinan, Jining, Qingdao, Liaocheng, Linyi, Zibo from Shandong and a sample from Nanjing. The pyrethroid-resistant and dichlorvos-resistant strains were reared in the laboratory. Single mosquito genomic DNA was extracted from Cx. pipiens pallens of these different populations and rNDA-ITS2 of Cx. pipiens pallens were specifically amplified by PCR. The PCR products were cloned into the pMD-19T simple vectors and sequenced. Genetic polymorphisms of them were analyzed by biological information systems. Results The obtained rDNA-ITS2 sequences were approximate 460 bp in length, with 36 mutation sites. Culex pipiens pallens from different geographical environment showed gene polymorphisms in sequence analysis of rNDA-ITS2. There were transition,transversion and deletion in the mono nucleotide of rNDA-ITS2 gene. Conclusion Genetic polymorphisms presented on rNDA-ITS2 of Cx. pipiens pallens from different geographical regions in Shandong province and different resistance strains in the lab were probably associated with the ecological environment.
Key words: Culex pipiens pallens     rDNA-ITS2     Genetic polymorphism    

蚊虫作为一类重要的医学昆虫,能传播疟疾、流行性乙型脑炎(乙脑)、登革热和丝虫病等多种虫媒传染病[1],媒介研究是蚊媒病流行病学研究和预防控制的首要工作之一。淡色库蚊(Culex pipiens pallens)是我国重要的吸血骚扰性蚊虫,在北方地区分布广泛,是乙脑和丝虫病的重要传播媒介。不同的蚊媒种群对于病原体的传播具有种属特异性,进而影响其流行病学特征,对蚊媒种属进行遗传变异研究和分子鉴别对于蚊媒疾病的防治具有重要意义。多拷贝的核糖体DNA是在蚊虫近缘种的鉴别和系统关系分析中应用较为成功和理想的分类指标。研究表明,蚊虫rDNA单位有18S、5.8S、28S RNA编码区、基因间隔区(ICS)、第1和第2转录间隔区(ITSl、ITS2 )、外转录间隔区(EST)组成,编码区的序列高度保守,间隔区的序列存在种间差异,其中ITS区不加入成熟核糖体,受到选择压力小,进化的速率较快,因此可以从其序列中得到较多的信息[2]。本实验根据核糖体DNA第二内转录间隔区(rDNA-ITS2)基因特点设计特异性引物,PCR扩增获得山东省不同地区不同地理环境下淡色库蚊的ITS2基因序列特征,探讨不同群体ITS2的基因分化及遗传差异,为山东省蚊媒监测及疾病防治工作提供依据。

1 材料与方法 1.1 材料 1.1.1 淡色库蚊样本

分别在济南(AJN)、济宁(JN)、青岛(QD)、聊城(LC)、临沂(LY)、淄博(ZB)以及南京(NJ)现场采集淡色库蚊雌蚊,实验室选育的抗敌敌畏(DDV)和抗氯氰菊酯品系(Pyr)淡色库蚊常规饲养。

1.1.2 试剂与主要仪器

dNTP、LA Taq聚合酶、pMD19-T Vector试剂盒、大肠埃希菌感受态细胞Escherichia coli JM109、DNA分子质量标准(DL-2000、1 kp Ladder)、T4 DNA连接酶均购自TaKaRa公司;质粒小提试剂盒和凝胶回收试剂盒购于Axygen公司;其余试剂为国产分析纯。T Gradient Thermoblock PCR扩增仪(Biometra),Mini-Protean Tetra Cell电泳仪(Bio-Rad),Smart Spec3000紫外分光光度计(Bio-Rad)。

1.2 方法 1.2.1 蚊虫基因组DNA提取

按文献[3]操作步骤,提取单蚊基因组DNA,具体步骤: 裂解液(NaCl 0.467 g,EDTA 2.233 g,SDS 0.5 g,Tris-HCl 1.121 g,蔗糖5.476 g,加H2O至100 ml),65 ℃预热。单个蚊虫置于1.5 ml离心管中,加入100 μl预热的裂解液,用塑料研磨棒研磨至匀浆状。65 ℃孵育30 min。加5 mol/L(25 μl)乙酸钾至终浓度1 mol/L,混匀。4 ℃置30 min。13 200 r/min(离心半径10 cm)离心10 min。取上清液,可抽取100~110 μl上清液,不抽沉淀。加2~2.5倍体积预冷无水乙醇(-20 ℃)沉淀核酸。4 ℃静置30 min至3 h。13 200 r/min (离心半径10 cm)离心10 min。弃上清液,加75%乙醇200~300 μl(-4 ℃) 洗涤盐分;弃上清液,37 ℃烘干。加50~100 μl双蒸灭菌水溶解,置-20 ℃冰箱保存待用。紫外分光光度计测DNA浓度及纯度,1%琼脂糖凝胶电泳分析。

1.2.2 淡色库蚊rDNA-ITS2基因的克隆

根据蚊虫rDNA-ITS2序列特征,在相对保守区域设计引物,引物由上海Invitrogen贸易有限公司合成。上游引物:5′- GGG GTA GTC ACA CAT TAT TTG AGG-3′,与淡色库蚊5.8S末端编码区保守序列互补;下游引物:5′-GAA CTG CAG GAC ACA TGA ACA CCG-3′,与淡色库蚊28S前端编码区保守序列互补。以提取的DNA为模板,用上述引物进行PCR扩增。采用50 μl反应体系:ddH2O 38.5 μl、10×LA PCR Buffer 5 μl、10 mmol/L dNTP 1 μl、模板DNA 1 μl、LA Taq聚合酶0.5 μl、10 pmol/L上下游引物各2 μl。反应条件:94 ℃预变性3 min; 94 ℃变性30 s,55 ℃退火45 s,72 ℃延伸1min,共30个循环;最后72 ℃延伸5 min。取扩增产物1 μl进行1%琼脂糖凝胶电泳分析。

1.2.3 PCR产物纯化测序及序列分析

用DNA片段纯化/回收试剂盒对目的DNA片段进行纯化,将纯化后的基因片段连接到pMD-19T simple Vector上,然后转化至大肠埃希菌感受态细胞E. coli JM109,并将转化菌涂布到含有氨苄西林的LB固体培养平板上,37 ℃培养过夜,挑取单菌落进行菌体PCR鉴定,阳性克隆摇菌培养并送Invitrogen公司测序。不同个体淡色库蚊rDNA-ITS2序列用DNAman软件进行序列比对,用ClustalX 1.81和Mega 4.0软件分析不同个体间的基因多态性和构建系统进化树。

2 结 果 2.1 淡色库蚊rDNA-ITS2克隆

以提取的6个不同地区淡色库蚊及实验室培养株的DNA为模板,PCR扩增得到rDNA-ITS2基因片段,1%琼脂糖凝胶电泳结果表明,在500 bp左右可观察到与预期分子质量大小相符的DNA条带(图1)。

注:M. DNA 分子质量标准(DL-2000);1~6. 济南、济宁、青岛、 淄博、临沂、聊城市蚊虫rDNA⁃ITS2 PCR产物。 图 1 淡色库蚊rDNA-ITS2基因PCR产物电泳结果 Figure 1 Analysis of PCR amplified rDNA-ITS2 genes of Cx. pipiens pallens
2.2 菌落PCR鉴定

将纯化的PCR产物与pMD-19T simple Vector连接并转化E. coli JM109,37 ℃过夜培养后挑选单菌落进行PCR鉴定阳性克隆,1%琼脂糖凝胶电泳分析(图2),其中在500 bp左右出现DNA条带的为阳性克隆。

注:M. DNA 分子质量标准(DL-2000);1~20. 不同菌落的PCR 鉴定结果。 图 2 菌落PCR的琼脂糖凝胶电泳分析 Figure 2 Analysis of colony PCR
2.3 淡色库蚊rDNA-ITS2基因测序及序列分析

测序结果发现9个克隆有9种序列类型,长度为460 bp左右。通过序列比对可以看出这些不同地区的淡色库蚊rDNA-ITS2基因序列存在着多态性(图3),淡色库蚊rDNA-ITS2出现36个变异位点,变异率为7.83%。rDNA-ITS2 基因单核苷酸存在转换、颠换和缺失,其中7处出现转换,包括4次A-G转换,3次T-C转换;4处出现颠换,包括1次T-G颠换,2次T-A颠换,1次C-G颠换;25处存在单核苷酸缺失(表1)。系统进化分析表明,山东省不同采样点淡色库蚊rDNA-ITS2基因存在不同分支,遗传距离存在差异,与实验室不同抗性品系rDNA-ITS2基因之间的遗传距离也存在差异,济南和济宁市聚成一类,聊城和淄博市聚成一类,青岛、临沂与南京市各为相对独立的一支,而抗DDV抗性和抗氯氰菊酯抗性的淡色库蚊和临沂地区淡色库蚊最为相近,是另外独立的一支(图4)。

图 3 淡色库蚊rDNA-ITS2基因序列比对结果 Figure 3 Sequence alignment of Cx. pipiens pallens rDNA-ITS2 · 553 ·
表 1 淡色库蚊rDNA-ITS2基因序列位点变异 Table 1 The mutation sites of rDNA-ITS2 gene sequences of Cx. pipiens pallens
图 4 山东省不同地区淡色库蚊rDNA-ITS2系统发育进化树 Figure 4 Phylogenetic analysis of Cx. pipiens pallens from different areas in Shandong province
3 讨 论

在蚊虫分类中,长期以来都是以形态特征作为分类依据,然而许多复合体的相继发现和亲缘种的存在,使得单纯的形态分类学无法对其进行鉴别。它们虽然外部形态相似,但其生态习性、分布区域和媒介效能等却可能有很大差异,并能影响其流行病学特征。因此,在难以用传统的形态学鉴别手段区分形态接近的蚊虫情况下,寻找新的有效方法对其进行正确鉴别具有重要的流行病学意义。随着分子生物学技术的迅猛发展,出现了线粒体DNA、核糖体DNA、微卫星DNA、单核苷酸多态性(SNP)和限制性片段长度多态性(RFLP)、随机扩增多态性(RAPD)等多种分子遗传标记。越来越多的研究人员使用分子生物和生物信息学方法鉴别蚊虫种群,并取得较好研究成果[4, 5, 6, 7]。山东省位于中国东部沿海、黄河下游,属暖温带季风气候类型,有多种蚊类孳生繁殖[8]。淡色库蚊是主要的蚊种,是丝虫病和乙脑的重要传播媒介,严重危害人类健康。充分运用分子生物学技术获得遗传信息,探索山东省蚊虫种群内的基因鉴别标志,对于深入了解媒介遗传多样性、媒介对病原体的易感性具有重要意义。

在蚊虫近缘种的鉴别和系统关系分析中应用较为理想的分类指标是多拷贝的核糖体DNA。rDNA-ITS2序列差异在种内高度保守,这一基因家族通过不等交换(unequal crossingover)和基因转换(gene conversion)经历快速的协同进化导致了重复单位在基因组中的一致性。rDNA-ITS2在生物界存在的普遍性、多拷贝性,在个体和群体内具有较好的均一性、少量样本能有效代表群体的变异情况等优点,很多学者将其作为生物系统进化研究的一个非常有用的分子标志[9, 10, 11, 12]。目前,蚊虫复合体成员种分子鉴别研究中应用rDNA作为分类指标较多,大多数以ITS2为依据[13]

本研究以山东省不同地区的淡色库蚊和实验室培养株为材料,设计rDNA-ITS2特异性引物,提取单蚊DNA,采用PCR方法获得rDNA-ITS2基因序列并进行序列分析。该方法虽简便快捷,但需要提取合格的蚊虫DNA,严格控制反应条件并精确测序。结果显示不同地区的淡色库蚊出现单核苷酸多态性,rDNA-ITS2基因变异率为7.83%,存在单核苷酸转换、颠换和缺失,缺失频率最高,转换频率高于颠换频率。变异位点均位于ITS2基因区,两侧的基因高度保守。不同地区淡色库蚊的rDNA-ITS2基因存在差异性,与实验室两个抗性品系之间也有较明显差异。这些差异与淡色库蚊对病原体易感性之间的关系有待于进一步研究。从生物进化树图中可以看出共形成两大分支,济南和济宁市的淡色库蚊遗传关系相近,聊城与淄博市的淡色库蚊遗传关系相近,青岛、临沂与南京市为相对独立的一支,这种基因的多态性可能与地域差异造成的生态环境以及人群活动和城市建设有关。抗DDV抗性和抗氯氰菊酯抗性的淡色库蚊来源于同一个分支,与两种抗性品系来源于同一种敏感品系的选育相符合,杀虫剂的选择压力使蚊虫向两个方向进化。二者和临沂地区淡色库蚊最为相近,可以预测其来源的敏感品系淡色库蚊可能与临沂地区的淡色库蚊遗传关系相近。

参考文献
[1] 李文刚,赵敏.虫媒传染病流行现状[J].传染病信息, 2011, 24(1):8-11.
[2] 江世宏,孟子烨,陈晓琴,等.核糖体DNA序列分析在昆虫系统学研究中的应用[J].昆虫分类学报, 2008, 30(3):225-238.
[3] Collins FH, Mendez MA, Rasmussen MO, et al.A ribosomal RNA gene probe differentiates member species of the Anopheles gambiae complex[J]. Am J Trop Med Hyg, 1987, 37(1):37-41.
[4] Paula RL, Ximo M.Mosquito species (Diptera:Culicidae) in three ecosystems from the Colombian Andes:identification through DNA barcoding and adult morphology[J]. Zookeys, 2015 (513):39-64.
[5] Ivanescu ML, Acatrinei D, Pavel I, et al.PCR identification of five species from the Anopheles maculipennis complex (Diptera: Culicidae) in North-Eastern Romania[J]. Acta Parasitol, 2015, 60(2):283-289.
[6] García-Vásquez A, Razo-Mendivil U, Rubio-Godoy M. Morphological and molecular description of eight new species of Gyrodactylus von Nordmann, 1832(Platyhelminthes:Monogenea) from poeciliid fishes, collected in their natural distribution range in the Gulf of Mexico slope, Mexico[J]. Parasitol Res, 2015, 114(9):3337-3355.
[7] 阳帆,张仁利,刘阳,等.深圳市白纹伊蚊细胞色素C氧化酶亚基Ⅰ基因多态性分析[J].中国媒介生物学及控制杂志, 2015, 26(4):337-340.
[8] Dai YH, Wang HF, Cheng P, et al.Mosquito species (Diptera: Culicidae) reported from Shandong province, China[J]. J Entomol Sci, 2011, 46(3):247-255.
[9] Sum JS, Lee WC, Amir A, et al.Phylogenetic study of six species of Anopheles mosquitoes in Peninsular Malaysia based on inter-transcribed spacer region 2(ITS2) of ribosomal DNA[J]. Parasit Vectors, 2014, 7:309.
[10] 钟敏,王英,付蒙颖,等.中缅边境地区部分传疟按蚊蚊种的形态学及分子鉴定[J].第三军医大学学报, 2014, 36(24): 2432-2436.
[11] Wassim NM.Secondray structure and sequence of ITS2-rDNA of the Egyptian malaria vector Anopheles pharoensis (Theobald) [J]. J Egypt Soc Parasitol, 2014, 44(1):197-204.
[12] 王海防,王怀位,程鹏,等.嗜血习性同域分化的中华按蚊种群rDNA-ITS2序列分析[J].中国血吸虫病防治杂志, 2014, 26 (5):526-530.
[13] Byrd BD, Harrison BA, Zavortink TJ, et al.Sequence, secondary structure, and phylogenetic analyses of the ribosomal internal transcribed spacer 2(ITS2) in members of the North American Signifera Group of Orthopodomyia (Diptera: Culicidae)[J]. J Med Entomol, 2012, 49(6):1189-1197.