王春花 黃威 周帥鋒 于盼輝 段素霞 張益 高立冬 馬學(xué)軍

102206 北京，中國(guó)疾病預(yù)防控制中心病毒病預(yù)防控制所(王春花、張益、馬學(xué)軍)；410005 長(zhǎng)沙，湖南省疾病預(yù)防控制中心(黃威、周帥鋒、高立冬)；050017 石家莊，河北醫(yī)科大學(xué)(于盼輝、段素霞)

·論著·

1株與手足口病病原混合感染的札如病毒全基因組序列測(cè)定與分析

王春花 黃威 周帥鋒 于盼輝 段素霞 張益 高立冬 馬學(xué)軍

102206 北京，中國(guó)疾病預(yù)防控制中心病毒病預(yù)防控制所(王春花、張益、馬學(xué)軍)；410005 長(zhǎng)沙，湖南省疾病預(yù)防控制中心(黃威、周帥鋒、高立冬)；050017 石家莊，河北醫(yī)科大學(xué)(于盼輝、段素霞)

王春花、黃威為共同第一作者

目的 了解手足口病樣本中除腸道病毒外的其他病毒病原，并分析其分子特征與基因類型。方法 采用宏基因組學(xué)結(jié)合生物信息學(xué)的方法對(duì)1例手足口病患兒樣本進(jìn)行深度測(cè)序分析。結(jié)果 該病例為腸道病毒71型與札如病毒混合感染。反轉(zhuǎn)錄聚合酶鏈反應(yīng)(RT-PCR)驗(yàn)證樣本中札如病毒核酸為陽(yáng)性，經(jīng)組裝拼接獲得全基因組序列為7 429 bp，含有3個(gè)開(kāi)放閱讀框。序列進(jìn)化分析表明，該株病毒與浙江株Hu/G1/Zhejiang1/China/2014的同源性最高為99.4%，屬于札如病毒GI基因型，是當(dāng)前我國(guó)主要的感染型。結(jié)論 利用深度測(cè)序方法鑒定出手足口病樣本混合感染病原為札如病毒，序列分析表明其為我國(guó)流行的優(yōu)勢(shì)毒株。該研究為今后我國(guó)手足口病和腹瀉相關(guān)疾病的病原檢測(cè)提供了參考依據(jù),為了解我國(guó)札如病毒感染的分子特征提供了序列參考。

Fund programs: National Key Research and Development Plan (2016TFC1202700, 2016YFC120090)

札如病毒(Sapovirus, SaV)屬于杯狀病毒科的成員，與諾如病毒(Norovirus, NoV)一樣，是非細(xì)菌性急性胃腸炎的主要病原體之一，可引起不同年齡段的人群發(fā)生腹瀉，特別是5歲以下的兒童[1]。與諾如病毒相比，札如病毒引發(fā)的感染沒(méi)有諾如病毒嚴(yán)重，暴發(fā)也不如諾如病毒常見(jiàn)，陽(yáng)性率也比諾如病毒低，故以往有關(guān)札如病毒的報(bào)道和相關(guān)研究較少。另外，札如病毒可與其他病毒發(fā)生混合感染，如輪狀病毒、星狀病毒等腸道病毒[2-5]。

札如病毒基因組為單股正鏈RNA，基因組核酸序列全長(zhǎng)約7.3-7.5 Kb，含有3個(gè)開(kāi)放閱讀框(Open Reading Frames, ORF)：ORF1、ORF2與ORF3，其中ORF1編碼RNA依賴的RNA聚合酶(RNA-dependent RNA polymerase, RdRp)和衣殼蛋白VP1。根據(jù)這兩個(gè)蛋白區(qū)的基因序列多樣性可將札如病毒分為GⅠ、GⅡ、GⅢ、GⅣ和GV 5個(gè)遺傳組[6]。其中GⅠ和GⅤ組又可進(jìn)一步分為8個(gè)和5個(gè)基因型[7]。近年來(lái)由札如病毒感染而引起的疫情在不斷增加，研究發(fā)現(xiàn)該病毒感染在全球范圍內(nèi)分布較廣，并且涵蓋多種基因型[8-12]。因此，對(duì)札如病毒進(jìn)行序列分析和基因型別鑒定，可為了解我國(guó)人群中札如病毒的流行趨勢(shì)與病毒感染情況和監(jiān)測(cè)病毒變異具有重要意義。本研究對(duì)1株與手足口病病原混合感染的札如病毒進(jìn)行了全基因組序列的測(cè)定與分析，并對(duì)其病毒核酸進(jìn)行了分子驗(yàn)證。

1 材料與方法

1.1 樣本采集與處理 2015年11月份在湖南省某醫(yī)院收集一名臨床診斷為輕癥手足口病的住院患兒糞便約10 g，立即冷凍至-20℃，然后低溫運(yùn)送至實(shí)驗(yàn)室待檢測(cè)。在生物安全柜中取出黃豆粒大小糞便溶于無(wú)菌磷酸鹽緩沖液(PBS)，經(jīng)超聲和離心后制備成20%的糞便懸液，分裝成200 μl/管，保存于-80℃。

1.2 核酸酶消化與病毒核酸提取 取200 μl便懸液進(jìn)行超速離心(32 000 r/min，3 h)，富集糞便樣本中的病毒。用100 μl PBS重懸沉淀，加入核酸酶cocktail，37℃孵育90 min。核酸酶cocktail包含10U TURBOTMDNase (Ambion), 25U Benzonase? nuclease (Sigma), 5U DNase I (TaKaRa), 2 μg RNase A (TaKaRa), 5U RNase T1 (Thermo Scientic), 4 U RNase One (Promega), 4 U Micrococcal nuclease (New England Biolabs)和15 U Mung Bean nuclease (Promega)。然后根據(jù)Qiaamp Viral RNA Mini Kit(Qiagen)提供的實(shí)驗(yàn)操作步驟提取樣本總RNA，分裝后保存于-80℃。1.3 cDNA文庫(kù)構(gòu)建 將以上獲得的總RNA利用SuperScript? Ⅲ Reverse Transcriptase kit (Invitrogen)聯(lián)合半隨機(jī)引物GGCAGGACCTCTGATACAGGNNNNNN，合成第一鏈cDNA。然后按文獻(xiàn)[12]中的方法依次進(jìn)行第二鏈cDNA合成和PCR擴(kuò)增。

1.4 深度測(cè)序 純化后的PCR產(chǎn)物以100 ng input按Ion Plus Fragment Library kit (Life Technologies)說(shuō)明書(shū)構(gòu)建測(cè)序文庫(kù)，主要步驟為末端修復(fù)、接頭連接、E-Gel片段篩選和qPCR定量。最后以26 pM文庫(kù)進(jìn)行油包水?dāng)U增、ES富集和Ion Torrent PGM測(cè)序儀測(cè)序。

1.5 深度測(cè)序數(shù)據(jù)分析 采用本實(shí)驗(yàn)室自主研發(fā)的病毒鑒定流程(Virus Identification Pipeline, VIP)[13]分析深度測(cè)序數(shù)據(jù)。使用CLC Genomic Workbench 9.0進(jìn)行病毒序列拼接。

1.6 札如病毒的驗(yàn)證 按One-Step RT-PCR kit(Qiagen)說(shuō)明書(shū)，用札如病毒特異性引物(SaV-F, 5′-CTCGCCACCTACRAWGCBTGGTT-3′,SaV-R, 5′-CGGRCYTCAAAVSTACCBCCCCA-3′)[14]檢測(cè)樣本核酸。

1.7 系統(tǒng)進(jìn)化分析 從GenBank下載目前已有的札如病毒全基因組序列或參考株的基因序列，進(jìn)行多序列比對(duì)分析。采用MEGA5.0軟件[15]和鄰接法(Neighbor-Joining)進(jìn)行系統(tǒng)進(jìn)化分析，并用Bootstrap重復(fù)抽樣1 000次對(duì)進(jìn)化樹(shù)進(jìn)行驗(yàn)證。

2 結(jié)果

2.1 深度測(cè)序數(shù)據(jù)概述 經(jīng)去除多克隆和低質(zhì)量數(shù)據(jù)，本研究共獲得66 266條reads,平均測(cè)序讀長(zhǎng)為151 bp。經(jīng)VIP分析，結(jié)果顯示病毒相關(guān)reads所占比例為29.5%。部分序列匹配到了腸道病毒71型和札如病毒且平均覆蓋度為99.9%和100%(圖1)，提示樣本中存在腸道病毒71型和札如病毒核酸，該手足口病患兒可能為腸道病毒71型和札如病毒混合感染。

圖1 腸道病毒71型與札如病毒的基因組覆蓋度和測(cè)序深度分析Fig.1 Viral genome coverage and sequencing depth analysis for enterovirus A71 and sapovirus

2.2 RT-PCR驗(yàn)證 用札如病毒特異性引物對(duì)樣本核酸進(jìn)行RT-PCR檢測(cè)，瓊脂糖電泳擴(kuò)增產(chǎn)物，結(jié)果顯示樣本核酸泳道有434 bp的目標(biāo)條帶。將該產(chǎn)物送公司進(jìn)行Sanger法測(cè)序，序列分析表明該樣本為札如病毒陽(yáng)性。

圖2 札如病毒VP1基因(A)和全基因組(B)核苷酸序列的遺傳進(jìn)化分析Fig.2 Phylogenetic analysis of the nucleotide sequence in VP 1 gene (A) and complete genome (B) of sapovirus.

2.3 札如病毒核酸序列分析 將深度測(cè)序數(shù)據(jù)進(jìn)行組裝拼接，獲得札如病毒全基因組序列，命名為SapovirusHu/GI/Hunan/China/2016，并將其全長(zhǎng)序列采用Sequin完成序列注釋和生成測(cè)序提交文件后上傳GenBank，獲得序列登錄號(hào)為KX980412。全基因組序列分析表明：該株病毒基因組全長(zhǎng)7 429 bp，包含3個(gè)ORF，其中ORF1編碼VP1蛋白。GC含量為50.81%，各堿基比例為：A(25.6%)、C(25.7%)、G(25.1%)、T(23.6%)。經(jīng)與GenBank中現(xiàn)有的札如病毒全基因組序列進(jìn)行多序列比對(duì)分析，結(jié)果顯示Sapovirus Hu/GI/Hunan/China/2016與GenBank中唯一的中國(guó)人源札如病毒Hu/G1/Zhejiang1/China/2014(KT327081)的同源性最高，核苷酸序列一致性為99.4%，氨基酸序列一致性為99.5%。

2.4 系統(tǒng)進(jìn)化樹(shù)分析 從GenBank中下載札如病毒VP1基因的代表性序列和札如病毒全基因組的代表性序列，與本研究獲得的病毒株相關(guān)序列進(jìn)行多序列比對(duì)分析，構(gòu)建系統(tǒng)進(jìn)化樹(shù)。如圖2所示，基于VP1基因(1 686 bp)和基于全基因組序列構(gòu)建的進(jìn)化樹(shù)，Sapovirus Hu/GI/Hunan/China/2016與其他代表株聚集為GI-GV 5個(gè)分支，并且均與Hu/G1/Zhejiang1/China/2014(KT327081)處于同一個(gè)遺傳進(jìn)化分支，表明該株病毒與Hu/G1/Zhejiang1/China/2014(KT327081)距離最近，同屬于GI基因型。

3 討論

本研究采用宏基因組學(xué)方法對(duì)一例手足口病輕癥患者的糞便樣本進(jìn)行了深度測(cè)序分析，結(jié)果表明該患者為腸道病毒71型與札如病毒混合感染。自1976年首次在日本扎幌市的一次暴發(fā)性腹瀉中發(fā)現(xiàn)該病毒以來(lái)，札如病毒與其他病毒混合感染的情況在世界范圍內(nèi)較為常見(jiàn)。我國(guó)在札如病毒陽(yáng)性樣本中檢測(cè)到了輪狀病毒和星狀病毒，混合感染率為20%。王顏歌等[16]在札如病毒陽(yáng)性樣本中分別檢測(cè)到了輪狀病毒、星狀病毒和諾如病毒，構(gòu)成比為18.75%。Wang等[17]在日本急性腸胃炎嬰兒和兒童的糞便樣本中也檢測(cè)到了札如病毒與輪狀病毒的混合感染，混合感染率為11.76%。近期有研究報(bào)道，札如病毒與諾如病毒的不同型別存在雙重感染或三重感染，并且此類混合感染主要為海鮮類食源性引發(fā)[18-20]。據(jù)我們所知，目前已報(bào)道的混合感染均來(lái)自于急性腸胃炎病例，并且均為札如病毒與腹瀉相關(guān)病毒共感染的現(xiàn)象。本研究首次在手足口病病例中發(fā)現(xiàn)了手足口病病原與札如病毒共感染的情況，為以后手足口病和腹瀉相關(guān)疾病的病原檢測(cè)提供了參考依據(jù)。

本研究還對(duì)混合感染的札如病毒序列進(jìn)行了組裝和拼接，獲得了7 429 bp全長(zhǎng)基因組序列。并將該基因組序列與GenBank中現(xiàn)存的參考株全序進(jìn)行了多序列比對(duì)和系統(tǒng)發(fā)育分析。結(jié)果顯示，無(wú)論核苷酸序列還是氨基酸序列，Hu/GI/Hunan/China/2016與札如病毒浙江株[21]的同源性最高；無(wú)論是基于VP1基因序列還是基于基因組全長(zhǎng)序列構(gòu)建的進(jìn)化樹(shù)，Hu/GI/Hunan/China/2016與札如病毒浙江株均位于相同的進(jìn)化分支。序列分析與進(jìn)化樹(shù)均表明Hu/GI/Hunan/China/2016為札如病毒GI型，與國(guó)內(nèi)多個(gè)地區(qū)流行的型別相同。我國(guó)深圳、海南、廣西、吉林、上海、陜西、浙江等地，札如病毒GI型是流行的優(yōu)勢(shì)毒株[4,16,20]。除GI型外，深圳等地區(qū)還報(bào)道了札如病毒GII型、GIV型。可見(jiàn)，札如病毒的流行具有多樣性。因此，本研究對(duì)札如病毒的序列進(jìn)行分析和型別鑒定對(duì)該病毒的診斷、監(jiān)測(cè)和防控具有重要意義。

[1] Svraka S, Vennema H, van der VB, et al. Epidemiology and genotype analysis of emerging sapovirus-associated infections across Europe[J]. J Clin Microbiol, 2010, 48(6): 2191-2198.doi:10.1128/JCM.02427-09.

[2] Okada M, Shinozaki K, Ogawa T, et al. Molecular epidemiology and phylogenetic analysis of Sapporo-like viruses[J]. Arch Virol, 2002, 147(7): 1445-1451. doi:10.1007/s00705-002-0821-7.

[3] Phan TG, Tuan AN, Nishimura S, et al. Etiologic agents of acute gastroenteritis among Japanese infants and children: virus diversity and genetic analysis of sapovirus[J]. Arch Virol, 2005, 150(7):1415-1424.doi:10.1007/s00705-005-0514-0.

[4] 常昭瑞,靳淼,劉娜,等.我國(guó)九省區(qū)2006年札如病毒流行狀況及基因序列分析[J].病毒學(xué)報(bào), 2009, 25(2): 113-115.doi:CNKI:SUN:BDXB.0.2009-02-006.

[5] Yan H, Nguyen TA, Phan TG, et al. Development of RT-multiplex PCR assay for detection of adenovirus and group A and C rotaviruses in diarrheal fecal specimens from children in China[J]. Kansenshogaku Zasshi, 2004, 78(8): 699-709.

[6] Dufkova L, Scigalkova I, Moutelikova R, et al. Genetic diversity of porcine sapoviruses, kobuviruses, and astroviruses in asymptomatic pigs: an emerging new sapovirus GIII genotype[J]. Arch Virol, 2013, 158(3):549-558.doi:10.1007/s00705-012-1528-z.

[7] Akihara S, Phan TG, Nguyen TA, et al. Identification of sapovirus infection among Japanese infants in a day care center[J]. J Med Virol, 2005, 77(4): 595-601.doi:10.1002/jmv.20497.

[8] Sdiri-Loulizi K, Gharbi-Khélifi H, de Rougemont A, et al. Acute infantile gastroenteritis associated with human enteric viruses in Tunisia[J]. J Clin Microbiol, 2008, 46(4):1349-1355.doi:10.1128/JCM.02438-07.

[9] Soares CC, Maciel de Albuquerque MC, Maranh?o AG, et al. Astrovirus detection in sporadic cases of diarrhea among hospitalized and non-hospitalized children in Rio De Janeiro, Brazil, from 1998 to 2004[J]. J Med Virol, 2008, 80(1):113-117.doi:10.1002/jmv.21053.

[10] Vernacchio L, Vezina RM, Mitchell AA, et al. Diarrhea in American infants and young children in the community setting: incidence, clinical presentation and microbiology[J]. Pediatr Infect Dis J, 2006, 25(1): 2-7.

[11] Lasure N, Gopalkrishna V. Epidemiological profile and genetic diversity of sapoviruses (SaVs) identified in children suffering from acute gastroenteritis in Pune, Maharashtra, Western India, 2007-2011[J]. Epidemiol Infect, 2017, 145(1): 106-114.doi:10.1017/S0950268816001953.

[12] Neill JD, Bayles DO, Ridpath JF. Simultaneous rapid sequencing of multiple RNA virus genomes[J]. J Virol Methods, 2014, 201: 68-72.doi:10.1016/j.jviromet.2014.02.016.

[13] Y Li, Wang H, Nie K, et al. VIP: an integrated pipeline for metagenomics of virus identification and discovery[J]. Sci Rep, 2016, 6: 23774. doi:10.1038/srep23774.

[14] 何雅青,卓菲,張海龍,等.一起札如病毒所致成年人急性胃腸炎暴發(fā)的分子流行病學(xué)研究[J].疾病監(jiān)測(cè), 2012, 27(2):101-103.doi:CNKI:SUN:JBJC.0.2012-02-007.

[15] Tamura K, Peterson D, Peterson N, et al. Mega5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods[J]. Mol Biol Evol, 2011, 28(10): 2731-2739.doi:10.1093/molbev/msr121.

[16] 王顏歌,何雅青,向文龍,等.深圳散發(fā)腹瀉患者札如病毒的檢測(cè)及型別分析[J].南方醫(yī)科大學(xué)學(xué)報(bào), 2011, 31(6):1017-1019.doi:CNKI:SUN:DYJD.0.2011-06-025.

[17] Wang Y, Zhang J, Shen Z. The impact of calicivirus mixed infection in an oyster-associated outbreak during a food festival[J]. J Clin Virol, 2015, 73:55-63.doi: 10.1016/j.jcv.2015.10.004.

[18] Pongsuwanna Y, Tacharoenmuang R, Prapanpoj M, et al. Monthly Distribution of Norovirus and Sapovirus in Viral Gastroenteritis in Thailand[J]. Jpn J Infect Dis, 2016,18.doi: 10.7883/yoken.JJID.2015.440.

[19] Gao Z, Li X, Yan H, et al. Human calicivirus occurrence among outpatients with diarrhea in Beijing, China, between April 2011 and March 2013[J]. J Med Virol, 2015, 87(12):2040-2047.doi: 10.1002/jmv.24265.

[20] Lu L, Jia R, Zhong H, et al. Molecular characterization and multiple infections of rotavirus, norovirus, sapovirus, astrovirus and adenovirus in outpatients with sporadic gastroenteritis in Shanghai, China, 2010-2011[J]. Arch Virol, 2015, 160(5): 1229-1238.doi: 10.1007/s00705-015-2387-1.

[21] Zhou X, Sun Y, Shang X, et al. Complete genome sequence of a sapovirus from a child in Zhejiang, China[J]. Virus Genes, 2016, 52(5):706-710.doi: 10.1007/s11262-016-1343-9.

Complete genome analysis of a sapovirus in a mixed infected patient of hand foot and mouth disease

WangChunhua,HuangWei,ZhouShuaifeng,YuPanhui,DuanSuxia,ZhangYi,GaoLidong,MaXuejue

NationalInstituteforViralDiseaseControlandPrevention,ChineseCenterforDiseaseControlandPrevention,Beijing102206,China(WangCH,ZhangY,MaXJ);HunanProvincialCenterforDiseaseControlandPrevention,Changsha410005,China(HuangW,ZhouSF,GaoLD);HebeiMedicalUniversity,Shijiazhuang050017,China(YuPH,DuanSX)WangChunhuaandHuangWeiarethefirstauthorswhocontributedequallytothearticle

MaXuejun,Email:maxj@ivdc.chinacdc.cn;GaoLidong,Email: 810173358@qq.com

Objective To understand the potential viral pathogens other than enteroviruses existing in samples of hand foot and mouth disease (HFMD) patient and study their molecular feature and genotype. Methods The deep sequencing analysis of a fecal specimen collected from HFMD patient was conducted by metagenomics and bioinformatics. Results Enterovirus A71 and sapovirus mixed infection was found in this case. The nucleic acid of sapovirus was confirmed positive by RT-PCR and the 7 429 bp complete genome sequence of sapovirus was obtained by assembling sequencing reads which consisted of 3 open reading frames. Phylogenetic analysis revealed that this strain of virus should belong to the genotype 1 of sapovirus having a homology of 99.4% with sapovirus Hu/G1/Zhejiang1/China/2014 strain, which is a currently predominant genotype circulating in China. Conclusions The sapovirus, which is a predominant strain circulating in China, was a mixed infected causative agent existing in HFMD sample identified by deep sequencing. This study will serve as a reference for pathogen detection of HFMD and diarrheal related diseases, as well as provide a sequence reference for molecular feature study of sapovirus in China in the future.

Sapovirus; Deep sequencing; Complete genome sequence; Phylogenetic analysis

國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016TFC1202700，2016YFC1200900)

馬學(xué)軍，Email: maxj@ivdc.chinacdc.cn；高立冬，Email: 810173358@qq.com

10.3760/cma.j.issn.1003-9279.2017.01.011

札如病毒；深度測(cè)序；全基因組序列；進(jìn)化分析

2016-05-11)