为了明确内蒙古锡林郭勒盟西乌旗羊株和锡林浩特市人株细粒棘球蚴之间的基因型及其遗传变异情况。
方法利用分子生物学技术对采自两地区细粒棘球蚴的线粒体CO1基因与ND1基因进行了克隆和测序,并运用DNAStar5.0中的MegAlign工具对已测出的序列进行了分析。
结果西乌旗羊株细粒棘球蚴与锡林浩特市人株细粒棘球蚴CO1基因片段和ND1基因片段大小分别均为936 bp和895 bp。西乌旗羊株细粒棘球蚴CO1基因序列与新疆羊株细粒棘球蚴的CO1基因序列同源性为99.3%;锡林浩特市人株细粒棘球蚴CO1基因序列与新疆人株细粒棘球蚴的CO1基因序列同源性为98.6%;且西乌旗羊株和锡林浩特市人株细粒棘球蚴ND1基因与国内外已报道的G1型ND1基因序列完全相同。
结论表明内蒙古上述两地区细粒棘球蚴的同源性很高,且基因型均为G1型,这为内蒙古某些地区细粒棘球蚴流行虫株的确定提供了重要依据,同时也对当地细粒棘球蚴病的预防和控制具有重要意义。
The purpose of this study was to clarify the strains of
细粒棘球蚴病,又名包虫病,是世界广泛分布的人兽共患性寄生虫病[ 1, 2],对感染动物的肉产量、肉品质、产奶量及后代繁殖均有严重影响[ 3],更重要的是该病对人类的健康也构成了严重的威胁[ 4, 5, 6, 7],因此一直倍受国内外学者的关注。近年来,相关学者已在基因水平上对细粒棘球蚴的虫株的遗传变异进行了研究。
据报道,虫株的遗传变异可能会影响到生理周期模式、宿主的特异性、感染性、抗原性、传播动力学、治疗药物的敏感性、致病性、诊断结果、流行病学和防控措施[ 8]。因此,确定虫株对制定棘球蚴病的预防和控制策略具有十分重要的意义。细粒棘球绦虫虫株的鉴定比较复杂,以往主要是以其成虫和幼虫的形态学特征,结合生物学、生物化学、同工酶分析和流行病学等方面的资料,作为虫株鉴别的依据。虽然形态学和生物学特征可为细粒棘球绦虫虫株鉴定提供很有价值的资料,但是由于受宿主和环境多种因素的影响,并不能反映不同虫株之间基因水平上的差异。近年来分子生物学技术的飞速发展,极大地促进了人们对细粒棘球绦虫虫株基因型鉴定的研究。众所周知,线粒体DNA(mtDNA)具有进化率快、重组率低、独立进化和基因组精简等特点[ 9, 10],可作为重要的分子标记,用于生物种群的遗传变异研究,尤其适于种以下水平的分类研究[ 11]。细粒棘球绦虫的线粒体基因组共编码2个核糖体RNA(12S、16S)、22个tRNA、1个细胞色素b(Cytb)、3个细胞色素氧化酶亚基(CO1、CO2、CO3)、NADH氧化还原酶6个亚基的基因(ND1、ND2、ND3、ND4、ND4L、ND6)和ATP酶一个亚基的基因(ATPase6)[ 12]。其中编码呼吸链的细胞色素c氧化酶和NADH氧化还原酶普遍存在于原核生物和真核生物细胞的线粒体中,是线粒体基因组中进化速率较慢的,可用以比较其在进化上的相互关系[ 13]。
鉴于以上因素,本研究以内蒙古锡林郭勒盟西乌旗地区检获的羊株和锡林浩特市地区检获的人株细粒棘球蚴作为研究对象,对其CO1与ND1基因的部分序列分别进行分析,并将两种序列与GenBank中已报道的细粒棘球蚴CO1和ND1基因相应序列进行比对,目的在于明确这些地区不同动物宿主间细粒棘球蚴种株的基因型和遗传变异情况及不同地理株之间的基因差异,这对了解家畜和人细粒棘球蚴病的感染来源以及不同宿主间的传播链并制定不同地区细粒棘球蚴病的具体防控措施,具有重要参考意义。
1.1.1 试验样品
细粒棘球蚴包囊采自于内蒙古锡林郭勒盟西乌旗绵羊和锡林浩特市某医院细粒棘球蚴病患者。将采集的样本用注射器将囊液吸出,离心1 min使囊液与原头蚴分离,-20 ℃条件下保存。
1.1.2 菌种、载体及主要试剂
大肠杆菌DH5α感受态细胞、DNA提取试剂盒(DNeasy Blood & Tissue Kit)和DNA凝胶回收试剂盒(AxyPrep DNA Gel Extraction Kit),购自北京市天根生化科技有限公司;Premix Taq、pMD19-T载体、DL2000 DNA Marker、氨苄青霉素(Amp)、5-溴-4-氯-3-吲哚-β-D半乳糖苷(X-gal)、异丙基-β-D-硫代半乳糖苷(IPTG)、10×Loading Buffer,均购自大连宝生物工程有限公司;核酸染料,购自北京市赛百盛基因技术有限公司。
1.2.1 总DNA提取
取实验室保存的细粒棘球蚴原头蚴,按照动物基因组DNA提取试剂盒说明书,提取基因组DNA,并置于-20℃冰箱中保存备用。
1.2.2 CO1与ND1基因扩增
以细粒棘球蚴DNA为模板,用上游引物5'-TTGTTAGGTGGTTTGTCTGA-3'和下游引物5'-GGCCATCATCAAATAAACAT-3',扩增CO1基因部分序列,PCR反应体系为50 μL,内含Permix Taq 25 μL,上下游引物各1.5 μL,模板5 μL,加dH2O补足至50 μL。同时,设置不加DNA模板的阴性对照。PCR扩增条件为:94 ℃预变性5 min后,进入循环;即95 ℃变性1 min,45 ℃退火45 s,72 ℃延伸70 s,共进行35个循环;最后72 ℃延伸10 min。
同样以细粒棘球蚴DNA为模板,用上游引物5'-GGTGGTTGTTTTTGGGTTAG-3'和下游引物5'-GTTTGCTATTGTTAATTAA-3',扩增ND1基因部分序列。PCR反应体系为同上,同时设阴性对照。PCR条件为:94 ℃预变性5 min,95 ℃变性1 min,50 ℃退火50 s,72 ℃延伸70 s,共进行35个循环,循环结束后以72 ℃继续延伸10 min。
1.2.3 PCR产物纯化、克隆与测序
利用胶回收试剂盒对PCR扩增产物的目的片段进行胶回收,将回收的CO1基因与ND1基因PCR扩增产物分别与pMD19-T载体连接。然后,将链接产物分别转化至DH5α感受态细胞中,再分别涂布于LB固体培养基(含有Amp、IPTG、X-gal)表面,37 ℃恒温培养12 h~16 h,后置4℃冰箱显色1 h~2 h。挑取数个白色菌落接种于含有Amp的LB液体培养基中,37 ℃、180 r/min振荡培养过夜,菌液经PCR鉴定后,将阳性菌液送至中美泰和公司测序。
细粒棘球蚴CO1和ND1基因分别经1%琼脂糖电泳后,在约936 bp、895 bp处呈现特异性条带,与目的片段大小基本一致,且无非特异性扩增条带产生,扩增结果见图1、图2:
西乌旗羊株细粒棘球蚴CO1基因与锡林浩特市人株细粒棘球蚴CO1基因片段大小均为936 bp,其A、T、G、C含量分别为15.71%、24.89%、48.72%、10.68%和15.60%、25.11%、47.54%、11.75%。
西乌旗羊株细粒棘球蚴ND1基因与锡林浩特市人株细粒棘球蚴ND1基因片段大小均为895 bp其A、T、G、C含量分别为19.44%、25.92%、46.59%、8.04%和19.44%、25.92%、46.82%、7.82%。
比对结果显示:西乌旗羊株细粒棘球蚴CO1基因与锡林浩特市人株细粒棘球蚴CO1基因分别同GenBank中已登录的新疆羊株细粒棘球蚴的CO1基因序列相比,相应区域同源性分别存在6个和13个变异位点,突变类型均为转换,即G与A置换和T与C置换;同源性分别为99.3%、98.6%(图3):
利用DNAStar5.0中的MegAlign工具对ND1基因的序列进行分析,结果显示:西乌旗羊株细粒棘球蚴的ND1基因序列与国内外已报道的G1型ND1基因序列完全相同;锡林浩特人株细粒棘球蚴的ND1基因序列与G1型ND1基因序列存在2个变异位点,全部为T与C转换;西乌旗羊株细粒棘球蚴和锡林浩特人株细粒棘球蚴均为G1基因型。
现有研究表明,动物mtDNA按照严格的母性遗传方式传代, 分子量小, 拷贝数多, 基因的碱基替代率比单拷贝的核DNA 基因组快5~10 倍,进化速率很快,其中CO1基因和ND1基因相当保守,表现为不同程度的种内遗传差异和种间差异。因此,CO1基因和ND1基因可用于动物性寄生虫的鉴定和不同地理株的区分[ 14, 15]。李明伟等(2005)[ 16]对犬弓首蛔虫的CO1基因部分序列进行了测序,结果表明不同地方虫株的CO1基因表现出一定的差异性。刘国华等(2009)用不同地区的泡状带绦虫的CO1基因序列与GenBank中报道的泡状带绦虫的CO1基因序列进行了对比,结果发现二者CO1基因序列无差异[ 17]。本研究对西乌旗羊株细粒棘球蚴和锡林浩特人株细粒棘球蚴CO1基因进行了测序,并将两株细粒棘球蚴CO1基因与新疆株细粒棘球蚴的CO1基因进行了比较。从序列对比结果可看出,西乌旗羊株细粒棘球蚴和锡林浩特人株细粒棘球蚴的CO1基因与新疆株细粒棘球蚴的CO1基因同源性很高,分别为99.3%、98.6%。
细粒棘球绦虫对中间宿主有广泛的适应性,但在长期进化演变的过程中,虫体与宿主之间在相互适应中也发生着变异。以往以形态学特征为主的虫株鉴定依据受到多种因素(如宿主和环境)的影响,并不能完全反映其基因水平的差异;而基于DNA序列的分析可直接鉴定生物体的不同种、亚种及地理株,能避免宿主和环境等因素的影响[ 18],ND基因是常用的分型基因。国内外学者根据基因型的差异,将细粒棘球绦虫分为10个基因型,即G1~G10。不同基因型虫株对人的感染性和致病性有差别,G1、G2、G3、G5、G6、G7、G8、G9均可感染人,其他基因型虫株对人的致病性有待进一步研究。不同基因型虫株间的核苷酸都有不同程度的差异,如G2和G3与G1间的核苷酸差异很小,被称为细粒棘球绦虫的狭义种;G4与其它虫株间的核苷酸差异在10%,G5与其它虫株间的核苷酸差异在6%以上;G6、G7、G8、G9、G10与G1间的核苷酸差异超过10%。杨俊克(2004)[ 19]对我国三省区(青海省、甘肃省和新疆维吾尔族自治区)细粒棘球蚴基因变异的情况进行了分析,由检测结果可知,各地区分离株ND1基因的变异率为0.1%~0.7%,均属于G1基因型。本研究对内蒙古两地区的细粒棘球蚴虫株进行了测序,经DNAStar5.0软件分析可知,两地区细粒棘球蚴均为G1型,以上数据表明,这些地区人的细粒棘区病有可能通过绵羊—犬—人方式进行传播流行,需要加以注意和防控。
综上所述,本次研究所得的西乌旗羊株细粒棘球蚴和锡林浩特市人株细粒棘球蚴的CO1基因与新疆株细粒棘球蚴的CO1基因相比均有变异位点,突变类型为转换,即不同地理株细粒棘球蚴CO1基因的部分序列略有不同;但在ND1基因方面,与新疆株相比,西乌旗羊株细粒棘球蚴ND1基因序列无变异位点,而锡林浩特市人株细粒棘球蚴ND1基因序列存在2个变异位点,全部为T与C转换。上述研究结果,为确定细粒棘球绦虫在内蒙古某些地区的流行虫株基因型及其可能的传播链情况提供了重要依据。
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|