广东省一株新型肠道病毒EV-B83型的全基因组序列进化分析及重组分析

本研究对2001年广东省急性弛缓性麻痹(Acute flaccid paralysis,AFP)病例中分离到的一株肠道病毒AFP341_GD-CHN_2001株进行了高通量测序后,获取该全基因组序列,经鉴定该毒株为EV-B83。为了解该广东省首株EV-B83分离株的全基因组特征及其重组特点,本研究对AFP341_GD-CHN_2001毒株通过最大似然法构建系统进化树,利用bootscanning分析该毒株与其它型别肠道病毒的重组位点,最后根据该重组位点分区段构建系统进化树,进一步验证重组分析结果。全长VP1序列构建的分子进化树表明,AFP341_GD-CHN_2001毒株与EV-B83原型株及柬埔寨SEP025_KH_2012分离株形成一簇;初步构建P1、P2、P3编码区序列系统进化树,提示该分离株在其分子进化中可能存在基因重组,而bootscanning分析表明非结构蛋白编码区是其发生基因重组的区域。根据bootscanning分析的断裂点对该毒株的全基因组序列分成3段构建系统进化树,该系统进化分析结果进一步明确该分离株于2B、2C、3D非结构编码区存在型间重组,与AFP341_GD-CHN_2001发生重组的毒株为CV-A9分离株NSW-V20_AU_2008、CV-B2分离株NSW-V53_AU_2010以及CV-B3分离株SSM_JL-CHN_2006。当前GenBank数据库中仅有来自美国以及中国云南的两株EV-B83全基因组序列,缺乏对该基因型的分子进化特征描述,本研究中的AFP341_GD-CHN_2001毒株为广东省首株EV-B83分离株,扩展了EV-B83数据信息,为今后EV-B83的研究提供了基础性资料。     

2009年云南省2株肠道病毒71型分离株全基因组序列遗传特性分析

对2009年云南省肠道病毒71型分离株KMM09和KM186-09进行全基因组序列测序,并与我国及其它国家流行的EV71基因型进行比较和进化分析。KMM09和KM186-09基因组长为7 409bp,编码2 193个氨基酸,VP1系统进化分析显示2009年云南分离株属于C4基因型的C4a亚型。在结构区,与其它基因型相比较,C基因型之间的核苷酸和氨基酸的同源性高于其它基因型;而在非结构区,C4与B基因型和CA16原型株G10同源性高于其它C基因亚型。通过RDP3重组软件和blast比对分析,发现EV71C4基因型与B3基因型,与CA16原型株G10的基因组在非结构区存在重组。EV71全基因组序列的比较和分析,对了解引起我国手足口病暴发或流行C4基因亚型EV71毒株的遗传特性具有重要意义。     

人肠道病毒A71型分离株全基因组序列及毒力相关位点分析

为了解人肠道病毒A71型(Enterovirus A71,EV-A71)的基因重组和遗传进化特征,探索其毒力相关位点。本研究通过对郴州市3株EV-A71分离株进行全基因组序列测定和分析,应用RDP软件(版本4.1.6)将其与231条EV-A71全基因组序列进行比较分析,检测可能的重组事件及重组位点;采用SimPlot软件(版本号3.5.1)对重组序列、2条亲本序列和其他A组肠道病毒(EV-A)代表株序列进行相似性分析;采用MEGA软件(版本5.2)构建系统发生树;结果发现郴州市3株EV-A71毒株均检测到重组信号,其共同的亲本主序列为广东省2009年毒株,亲本亚序列为中国台湾地区2008年毒株;我国C4b进化分支EV-A71代表株(分离于1998年的SHZH98株)和C4a进化分支EV-A71代表株(分离于2008年的安徽阜阳株)均与CVA4、CVA14及CVA16原型株在3D区域检测到型间重组信号。采用SPSS软件(版本19.0)对不同病例类型的变异位点作卡方检验以筛选出可能的毒力相关位点。发现死亡病例组与重症病例组间未找到具有统计学意义的突变位点;而死亡与重症病例合并组与一般病例组进行比较,共发现18个突变位点具有统计学意义,这些突变位点仅分布于P2区和P3区,而P1区未见。我国大陆的EV-A71流行株与EV-A其他毒株的型间重组早在1998年前即已经发生,并且在2008~2012年的手足口病暴发疫情中,EV-A71的型内重组频繁发生。此外,EV-A71毒力的改变可能与非衣壳蛋白区的变异相关,因此应该关注非衣壳蛋白基因改变对病毒毒力的影响。     

人肠道病毒A71型分离株全基因组序列及毒力相关位点分析北大核心CSCD

为了解人肠道病毒A71型（Enterovirus A71,EV-A71）的基因重组和遗传进化特征,探索其毒力相关位点。本研究通过对郴州市3株EV-A71分离株进行全基因组序列测定和分析,应用RDP软件（版本4.1.6）将其与231条EV-A71全基因组序列进行比较分析,检测可能的重组事件及重组位点;采用SimPlot软件（版本号3.5.1）对重组序列、2条亲本序列和其他A组肠道病毒（EV-A）代表株序列进行相似性分析;采用MEGA软件（版本5.2）构建系统发生树;结果发现郴州市3株EV-A71毒株均检测到重组信号,其共同的亲本主序列为广东省2009年毒株,亲本亚序列为中国台湾地区2008年毒株;我国C4b进化分支EV-A71代表株（分离于1998年的SHZH98株）和C4a进化分支EV-A71代表株（分离于2008年的安徽阜阳株）均与CVA4、CVA14及CVA16原型株在3D区域检测到型间重组信号。采用SPSS软件（版本19.0）对不同病例类型的变异位点作卡方检验以筛选出可能的毒力相关位点。发现死亡病例组与重症病例组间未找到具有统计学意义的突变位点;而死亡与重症病例合并组与一般病例组进行比较,共发现18个突变位点具有统计学意义,这些突变位点仅分布于P2区和P3区,而P1区未见。我国大陆的EV-A71流行株与EV-A其他毒株的型间重组早在1998年前即已经发生,并且在2008~2012年的手足口病暴发疫情中,EV-A71的型内重组频繁发生。此外,EV-A71毒力的改变可能与非衣壳蛋白区的变异相关,因此应该关注非衣壳蛋白基因改变对病毒毒力的影响。     

2010年云南省1株ECHO-9病毒全基因序列的遗传特性分析

为了解云南省手足口病患者标本中分离到的一株ECHO-9病毒基因组特征,对2010年云南省ECHO-9病毒分离株MSH-KM812-2010全基因组序列测序,并与GenBank中其它ECHO-9病毒株基因组序列比对和分析。MSH-KM812-2010基因组长为7 424bp,编码2 203个氨基酸,其结构基因区与其它ECHO-9病毒核苷酸和氨基酸的同源性高于其它型别肠道病毒;而在非结构区,与其它肠道病毒血清型的同源性高于ECHO-9病毒株。VP1系统进化分析显示ECHO-9病毒株可形成A、B和C三个进化分支,MSH-KM812-2010株以及其它中国分离株属于C簇。三个分支之间核苷酸序列的差异大于15.0%可将ECHO-9病毒分为三个基因型。通过重组检测软件3(RDP3)和基本局部比对搜索工具(BLAST)比对分析,发现在非结构区可能存在重组。本文首次对我国分离的ECHO-9病毒全基因组序列的测定和分析,对了解ECHO-9病毒遗传特性具有重要意义。     

2014年安徽省柯萨奇病毒A组16型分离株VP1基因特征研究

研究2014年安徽省手足口病(Hand,foot and mouth disease,HFMD)患儿中分离的柯萨奇病毒A组16型(Coxsachivirus A 16,CVA16)毒株VP1区基因特征。收集安徽省2014年1月至11月期间413份HFMD患儿咽拭子标本接种敏感细胞分离肠道病毒,用荧光定量RT-PCR鉴定细胞培养物,对CVA16阳性培养物进行毒株VP1区RT-PCR扩增及核苷酸序列测定和基因特征分析,并与国内外参考序列构建基因亲缘性关系树。共分离鉴定出肠道病毒阳性培养物97份,其中CVA16病毒分离株17株,人肠道病毒71型(Human enterovirus 71,HEV71)分离株76株,其它肠道病毒分离株4株,总体病毒分离率为23.49%(97/413)。CVA16基因亲缘性关系分析显示:安徽省2014分离的17株CVA16毒株都属于B1基因型B1b一个分支,它们之间在核苷酸和氨基酸水平上的同源性分布为分别为95.30%~100%和98.70%~100%,但在B1b分支内形成几个传播链。17株CVA16毒株VP1区核苷酸序列与国内云南、湖南、广东、西藏和江苏地区病毒分离株同源性较高,亲缘性关系近,其中与湖南2013年和广东深圳2014年CVA16病毒分离株同源性最高,为96.40%~99.70%。2014年安徽省分离的CVA16病毒分离株属于B1基因型B1b亚型,为优势流行株;在B1b分支内形成多个小的病毒传播链共同流行。     

中国6株狂犬病病毒街毒株全基因组测序与分析

实验研究中对分离于中国的6株狂犬病病毒街毒株进行了全基因组测序,对基因组的5个结构基因(N、P、M、G和L)的核苷酸和推断的氨基酸序列以及非编码区序列进行了分析与比较,并与来自GenBank的40株毒株从全基因组水平进行了分子进化分析。所测6株中国狂犬病病毒街毒株的全基因组核苷酸序列长度介于11 907 nt(CQ92)和11 924 nt(SH06和gg4)之间,基因组结构相同,用全基因组和不同的结构基因构建的进化树拓扑结构相似,基因组3′和5′末端高度保守而且末端11个核苷酸互补配对,5个结构基因的保守性依次是NLMGP,核苷酸同源性的最小值依次分别是81.9%、81.7%、80.7%、78.3%和76.7%。     

山羊关节炎脑炎病毒甘肃株全基因组克隆和序列分析

根据GenBank中发表的山羊关节炎脑炎病毒CAEV-CO株（caprine arthritis encephalitis virus-CO,CAEV-CO）的全基因组序列（序列号：NC-001463）,设计合成了7对引物,对CAEV甘肃株的全基因组进行了PCR扩增,并对扩增产物进行了克隆和测序。结果表明,CAEV-甘肃株基因组全长为9186nt。与国际标准毒株CAEV—CO株相比,在编码区有12个碱基的缺失,二者的核苷酸同源性为91．0％;在非编码区有27个核苷酸的插人,核苷酸同源性为97．0％。gag、pol、蛋白Q、tat、env基因编码的氨基酸同源性分别为94．6％、94．7％、87．9％、94．7％和91．5％。     

福建一例HAstV-5型星状病毒的全基因组测序与重组分析

为了探究一例福建省检出的HAstV-5型星状病毒2013/Fuzhou/85毒株基因组分子结构特点,本研究采用PCR分段扩增、测序、拼接的方法,获得2013/Fuzhou/85毒株基因组序列全长6 803bp:5’端和3’端均有85bp非编码区;中间3个开放阅读框:ORF1a长2 802bp(86～2 887nt),编码非结构蛋白丝氨酸蛋白酶;ORF1b长1 548bp(2 827～4 374nt),编码非结构蛋白RNA聚合酶;ORF2长2 352bp(4 367～6 718nt),编码结构蛋白衣壳蛋白前体。目前,GenBank中仅有两株HAstV-5型星状病毒全基因组序列:中国辽宁毒株(JQ403108)和巴西哥亚尼亚毒株(DQ028633),2013/Fuzhou/85毒株和中国辽宁毒株核苷酸相似度最高,达94.4%。对该HAstV-5型星状病毒3个开放阅读框分别构建系统进化树,发现ORF1a与HAstV-1(JF327666)相似度最高,ORF1b和ORF2与HAstV-5(JQ403108)相似度最高,提示其有可能存在重组,用Simplot软件进行重组分析,重组位点位于2 741bp,在ORF1a和ORF1b重叠区的上游。本研究中对2013/Fuzhou/85毒株的全基因组测序和重组分析,可以为星状病毒的重组和遗传进化规律研究提供参考。     

2012年江苏省柯萨奇病毒A组14型分离株全基因序列分析

目的探讨2012年江苏省手足口病(HFMD)病原体CVA14病毒的全基因序列特点。方法采集江苏省邳州市2012年HFMD患儿的咽拭子或肛拭子样本,分离培养得到2株CVA14病毒,利用RT-PCR扩增CVA14全长基因。利用DNAStar软件包的Meg Align进行核苷酸序列分析,构建系统进化树。结果扩增得到覆盖CVA14全基因组的3个片段,拼接后为CVA14全基因组,分别命名为PZ05Y/JS/2012和PZ15G/JS/2012。同源性分析结果显示,这2株分离株之间同源性为99.4%,原型株(G-14)的核苷酸同源性为84.4%。进化树结果显示,所有的CVA14形成A、B和C共3个分支。结论获得了江苏省CVA14的全基因序列,分析了该毒株的进化特点,补充了中国CVA14病毒的基因信息。     

一起无菌性脑膜炎疫情中埃可病毒30型的全基因组分析

对引起一起无菌性脑膜炎的埃可病毒30型（（Echovirus 30,E30）毒株进行全基因组测序,并分析其遗传变异和分子进化特征。提取病毒RNA,荧光RT-PCR确定为肠道病毒,再扩增其VP1区,测序确定为E30,设计针对E30全基因组的引物,RT-PCR扩增和序列测定获得全基因序列。利用DNAMAN 9.0、MEGA X、RDP 5和SimPlot 3.5.1软件分析全基因序列。E30无锡株基因组全长7 425 bp个核苷酸（nt）,5′端和3′端分别为743 nt和97 nt的UTR,二个UTR之间为一个6 585 nt长的开放阅读框,编码一个含2 195个氨基酸（aa）的多聚蛋白。与GenBank中基因组序列比对,最同源的是毒株USA/2017/CA-RGDS-1048 （基因登录号：MN153801）,核苷酸同源性为97.5%,氨基酸同源性为99.1%。VP1基因进化树分析显示,E30无锡株属于h型,并可归类于h型下分的基因簇GroupⅣ。种系进化分析和同源性分析提示E30无锡株可上溯到中国江苏毒株FDJS03毒株。分析还发现E30无锡株在非结构区存在重组,重组序列可能来自E3...     

急性弛缓性麻痹病例中Ⅱ型脊髓灰质炎疫苗病毒变异株基因特征分析

研究Ⅱ型脊髓灰质炎(脊灰)疫苗变异株的基因特征,为我国使用口服脊灰减毒活疫苗/脊灰灭活疫苗使用策略,维持无脊灰状态和全球最终消灭脊灰提供科学依据。根据型内鉴定的检测结果,从2000~2001年AFP病例分离到的Ⅱ型脊灰疫苗变异株中选取有聚集性的5株病毒进行全基因组序列测定(贵州省3株、山东省2株),并进行核苷酸、氨基酸同源性分析。贵州省3株病毒全基因组序列完全一致,但与SabinⅢ型病毒发生重组,重组区域在3A区(nt5343~5353);与疫苗株相比,Ⅱ型区域变异10个碱基,其中VP1区变异4个,与SabinⅡ型株核苷酸同源性为99·56%,氨基酸同源性99·34%;Ⅲ型区域变异9个碱基。山东省2株病毒全基因序列共享16个突变位点,没有发生重组,与SabinⅡ型株相比,VP1区分别变异7个和4个碱基,核苷酸同源性分别为99·22%和99·56%,氨基酸同源性分别为99·0%和99·67%。上述5株病毒在重要的减毒位点nt481、nt2909均发生突变。此研究中5株病毒分属于两个不同的传播链,但是共享nt481、nt2909、nt2992三个突变位点,这3个突变位点不在重组区域内,他们的共同作用可能是影响病毒传播力的重要因素,但目前尚无证据证明脊灰疫苗病毒型间重组会增加病毒的毒力及传播力。     

浙江地区鼬獾和犬源狂犬病病毒分离株全基因组测序与分析

测定浙江地区狂犬病病毒分离株(鼬獾和犬)全基因组序列,从分子水平对病毒进行遗传变异特征分析,了解狂犬病病毒在浙江的流行和变异情况以及目前浙江流行株的遗传学背景,以丰富中国狂犬病病毒街毒流行株的全基因组信息。脑内接种1~2日乳鼠分离狂犬病病毒,RT-PCR反应测定浙江地区狂犬病病毒分离株全基因组核苷酸序列,并进行编码蛋白和序列相似性比较及种系发生分析。测序获得狂犬病病毒浙江淳安鼬獾分离株F02、F04和松阳犬分离株D01、D02全基因组核苷酸序列信息:基因组全长11 923和11 925 nts,前导序列Leader长58nts,5个ORF为:NP(1 353 nts);PP(894 nts);MP(609 nts);GP(1 575 nts);LP(6 386 nts),N-P-M-G间隔序列长2、5、5 nts;G-L基因间的伪基因Ψ长423 nts;Trailer尾长70 nts。核酸BLAST及多重序列比对分析显示浙江地区4个狂犬病病毒分离株的全基因组序列的组成和结构符合弹状病毒科狂犬病病毒属的特征;中国毒株之间特别是浙江同种动物狂犬病病毒之间各个基因区域核苷酸与氨基酸序列相似性最高,浙江病毒全基因组序列编码蛋白氨基酸序列相似性高于核苷酸序列相似性,说明蛋白质编码基因的核苷酸变异大多属于同义突变;浙江病毒负链RNA基因组5个基因编码氨基酸的长度没有变异,5个编码蛋白仅表现较少的序列变化;浙江病毒与本研究选择的代表性引用街毒株或者来自街毒的减毒株的变异位点和变异类型相似,多重序列相似性的比较和种系发生分析显示所分离的狂犬病病毒浙江街毒株均属于基因1型,具有较独特的中国地域性特点,故本研究中的浙江地区分离株极有可能是自然界中固有的街毒株。     

一株辛德毕斯病毒的基因组序列测定及分析

目的：对引进的一株辛德毕斯病毒的基因组序列进行测定,阐明其与已报道毒株序列的关系。方法：对辛德毕斯病毒基因组编码区进行分段RT-PCR扩增,对非编码区采用RACE法进行扩增,将扩增产物直接进行测序,应用DNAStar软件将测序结果拼接得到基因组序列,采用MEGA3.1软件对9株辛德毕斯病毒基因组序列进行系统进化发生树的构建。结果与结论：此株辛德毕斯病毒基因组共11663nt,编码3745个氨基酸残基,其中5＇端的2/3基因组编码4种非结构蛋白NSp1、NSp2、NSp3和NSp4,3＇端的1/3基因组编码5种结构蛋白E1、E2、E3、6K和C;结构基因和非结构基因之间有48nt的连接区为非翻译区;病毒基因组5＇末端和3＇末端分别有59、318nt的非编码区;序列同源性分析结果表明,此株病毒与S.A.AR86株的同源性最高,两者核苷酸序列的同源性为99.7%,氨基酸序列的同源性为99.6%,而与本室保存的另一辛德毕斯病毒MEI株的遗传进化关系稍远,系统进化发生树处于不同分支上。     

三株新城疫广西分离株全基因组序列的测定与分析

根据GenBank上所公布的新城疫病毒(NDV)的全基因组序列,设计了8对引物,运用RT-PCR方法获取了3株广西地方强毒株GX7/02、GX9/03和GX11/03的全基因组序列,并对其进行了比较分析。此3株病毒的全基因组序列均由15192个碱基组成,与GenBank公布的ZJ1、U.S/Largo/71、Italy/2736/00等7个毒株的全基因组序列长度相同,比LaSota、Clone-30和B1的全基因组序列多出6个核苷酸,此6个核苷酸位于np基因的非编码区内,相对与NDV毒株LaSota、Clone-30和B1序列的1647~1648nt位。通过序列的比较分析,发现GX7/02、GX9/03和GX11/03与ZJ1毒株的同源性较高,而与LaSota、Clone-30和B1等毒株的同源性较低。     

三株新城疫广西分离株全基因组序列的测定与分析

根据GenBank上所公布的新城疫病毒(NDV)的全基因组序列,设计了8对引物,运用RT-PCR 方法获取了3株广西地方强毒株GX7/02、GX9/03和GX11/03的全基因组序列,并对其进行了比较分析.此3株病毒的全基因组序列均由15192个碱基组成,与GenBank公布的ZJ1、U.S/Largo/71、Italy/2736/00等7个毒株的全基因组序列长度相同,比LaSota、Clone-30和B1的全基因组序列多出6个核苷酸,此6个核苷酸位于np基因的非编码区内,相对与NDV毒株LaSota、Clone-30和B1序列的1647～1648nt 位.通过序列的比较分析,发现GX7/02、GX9/03和GX11/03与ZJ1毒株的同源性较高,而与LaSota、Clone-30和B1等毒株的同源性较低.     

辽宁省首株柯萨奇病毒B组5型全基因组序列分析

研究引起辽宁地区手足口病的柯萨奇病毒B组5型（coxsackievirus B5,CV-B5）基因组特征。对2018年从辽宁省688份肠道病毒核酸阳性的标本中分离到的1株CV-B5进行高通量测序，并对其全基因组进行遗传进化分析。结果表明，CV-B5辽宁分离株与国内流行株的全基因组核苷酸序列同源性为78.5%～97%，氨基酸序列同源性为75.3%～96.7%。基于全基因组的进化分析将CV-B5流行株分为A～D四个基因型，辽宁分离株属于D基因型。通过重组分析发现其在P3区的3D区段发生重组。首次在辽宁地区手足口病患儿中分离出CV-B5，辽宁省分离株（LN2018-23-21/CHN/2018）可能为重组株。     

FMDV OH99株基因组全序列的测定及其基因特征研究

采用RT-PCR方法对FMDV OH99株基因组全序列进行了分子克隆与测序。结果表明OH99株基因组全基因组序列长8040nt,其中5’NCR长1026nt,前导蛋白(L)编码区长603nt。该毒株结构蛋白与非结构蛋白编码区的核苷酸序列为6318nt,3’NCR长93nt,其后是poly(A)尾巴,测序结果表明该结构至少含有56个A。应用分子生物学软件,将OH99株与其它参考毒株进行了序列比较,并对其基因特征、推导的氨基酸序列进行了研究分析。结果显示,在分类地位上OH99株归属于O型FMDV,与OTY TW／97具有较高的同源性,而与其他参考毒株的差异性比较大,而且在基因组功能未知区域和3A编码区域具有两处明显的基因片段缺失现象,其中3A编码区缺失30nt,与OTY TW/97株相同,但功能未知区域的缺失状况与OTY TW／97稍有差异。根据VP1基因序列,对OH99株与参考毒株进行了系统发生树分析,分析结果表明OH99株与0TY TW／97株在同一基因型内,其遗传关系最近,而与其毒株遗传关系较远。     

中国脊髓灰质炎病毒疫苗株基因变异的分析

从急性弛缓性麻痹（ＡＦＰ）病例分离的３２株脊髓灰质炎（脊灰）病毒,在ＶＰ１编码区,经ＰＣＲ－ＲＦＬＰ方法鉴定,并经核酸测序进一步证实力疫苗株,其中５株为Ｉ型疫苗,１６株Ⅱ型疫苗株,１１株Ⅲ型疫苗株,以同样的分析方法检测这些毒株基因组的３Ｄ聚合酶编码区,发现２株Ｉ型疫苗株在３Ｄ区的４３１个核苷酸序列为野毒序列,即这２株Ｓａｂｉｎ１基因型（ＶＰ１）与野毒基因型（３Ｄ）重组株;其余３株Ｉ型疫苗株未发现基     

FMDV OH99株基因组全序列的测定及其基因特征研究

采用RT-PCR方法对FMDV OH99株基因组全序列进行了分子克隆与测序.结果表明OH99株基因组全基因组序列长8040nt,其中5'NCR长1026nt,前导蛋白(L)编码区长603nt.该毒株结构蛋白与非结构蛋白编码区的核苷酸序列为6318nt,3'NCR长93nt,其后是poly(A)尾巴,测序结果表明该结构至少含有56个A.应用分子生物学软件,将OH99株与其它参考毒株进行了序列比较,并对其基因特征、推导的氨基酸序列进行了研究分析.结果显示,在分类地位上OH99株归属于O型FMDV,与OTY TW/97具有较高的同源性,而与其他参考毒株的差异性比较大,而且在基因组功能未知区域和3A编码区域具有两处明显的基因片段缺失现象,其中3A编码区缺失30nt,与OTY TW/97株相同,但功能未知区域的缺失状况与OTY TW/97稍有差异.根据VP1基因序列,对OH99株与参考毒株进行了系统发生树分析,分析结果表明OH99株与OTY TW/97株在同一基因型内,其遗传关系最近,而与其毒株遗传关系较远.