The complete genome of Enterovirus 71 China strain

Five overlapping clones covering the full genome of Enterovirus 71 China strain SHZH98 were obtained and then the sequences were determined by the chain termination method. It showed that the full length of EV71 SHZH98 genome (not including Poly A tail) is 7408 bp. There are some diversities on the lengths and sequences of 5′ UTR and 3′ UTR between SHZH98 and the other EV71 strains. In P1 capsid region, which is closely associated with viral immunogenicity, EV71 strain SHZH98 shares the highest homology with Taiwan strains; but in P2 and P3 non-structural gene regions there are higher identities with Coxsakievirus A16 and EV71 strains MS, BrCr than with Taiwan strains. Phylogenetic tree constructed by structural gene region indicates that China strain SHZH98 has a closer relationship with Taiwan strains, however, in the non-coding region it has a closer relationship with Coxsakievirus A16, EV71 strains MS and BrCr. EV71 China strain was analyzed at the molecular level. The results will contribute to the basic study on enteroviruses and the EV71 prevention in China.     

我国分离的肠道病毒71型(SHZH03病毒株)全基因组核苷酸序列分析

对肠道病毒71型(enterovirus 71,EV71)中国(深圳)分离株SHZH03进行了全基因组(未包括多聚腺苷尾)7406个碱基的核苷酸序列测定.结果表明,SHZH03株与其它肠道病毒71型毒株相比,在编码区没有核苷酸的缺失和插入,其5′UTR和3′UTR区的长度和序列有一定的差异.核苷酸同源性比较结果表明,在P1区SHZH03株与SHZH98株、中国台湾流行株(TW2086、TW2272)的同源性较高(分别为92.5%,90.1%和87.9%),与新加坡流行株SIN5666、SIN5865及标准株MS、BrCr的同源性则在81%左右,而与Coxsackievirus A16(Cox.A16)的同源性最低(63.6%).氨基酸同源性比较结果表明,在P1区SHZH03株与Cox. A16的同源性最低,但在P2和P3区SHZH03株与Cox.A16的同源性最高.P1区的遗传进化分析表明,SHZH03株和中国台湾1998年流行的EV71毒株的亲缘关系较近,属于同一型(genogroup),而与标准株BrCr和MS的亲缘关系较远.上述结果有助于肠道病毒71型的基础研究和中国对于EV71所致疾病的预防.     

柯萨奇病毒A组16型中国分离株(Cox.A16 SHZH00-1)全基因组序列测定及分析

对分离自2000年中国深圳地区手足口病患儿粪便标本的柯萨奇病毒A组16型(Coxsackievirus A16,Cox．A16)病毒SHZH00—1株进行全基因组序列测定及分析后发现,Cox．A16SHZH00—1的基因组(未包括多聚腺苷酸尾)长度为7410bp,其中5’端非编码区(5’UTR)长为745bp,病毒基因组编码区全长6582个核苷酸,编码一个含2193个氨基酸残基的多聚蛋白,3’端非编码区长为83bp。Cox．A16SHZH00—1基因组的结构与Cox．A16亚洲地方株Tainan-5079-98(AFl77911)十分相近,整个基因组的核苷酸同源性为97．5％,氨基酸同源性为98．8％;而与Cox．A16国际标准株G10(U05876)差别较大,两者核苷酸同源性仅为79．1％,氨基酸同源性为94．5％。Cox．A16SHZH00—1与两株肠道病毒71型SHZH03(AY465356)和BrCr(U22521)比较,核苷酸同源性均不超过80％。     

肠道病毒71型SHZH03株VP1基因的进化分析

构建了肠道病毒71型（EV71）中国（深圳）分离株SHZH03全基因组的8个相互重叠的克隆,对其全基因组7406bp的核苷酸进行序列测定,利用DNA—Star软件分析外壳蛋白基因VP1的遗传进化。结果表明,SHZH03和SHZH98与亚洲流行株中的台湾1998年流行株、日本1999年流行株的遗传距离较近,而与新加坡2000年和2001年流行株的遗传距离较远;SHZH03株与一些欧洲流行株有较大的差异。以上结果说明我国深圳地区流行的肠道病毒71型有可能来源于台湾1998年的EV71大规模流行时的毒株。     

肠道病毒71型基因组中小干扰RNA的靶位点预测

肠道病毒71型（EV71）已经在世界范围内有过十多次大的爆发与流行,近年来EV71病毒的流行在亚洲逐渐呈上升的趋势,但是目前尚无有效的治疗措施,因此迫切需要一种治疗EV71的有效药物。本文采用生物信息学的方法,对人类EV71病毒三个不同株型（SHZH03,SHZH98和MS）RNA序列的局域二级结构进行了预测,并从这些病毒株的基因组中分别得到了长度在21～25nt的小干扰RNA靶序列碱基片段。这一结果将有助于治疗EV71药物的开发研究,对预防和控制EV71的爆发和流行也会有重要意义。     

3株肠道病毒71型毒株的全基因组序列分析

目的：对获得的3株肠道病毒71（EV71）型毒株进行全基因组序列测定,并对其进化特点及分型进行初步分析。方法：提取病毒RNA,反转录得到eDNA,PCR分段扩增覆盖病毒全长序列的6个重叠片段（不包括多聚腺苷酸尾）;用软件将3株EV71的备片段序列进行拼接、编辑和校正,随后进行氨基酸翻译及序列比较;用MEGA4.1软件构建系统进化树。结果：获得了3株EV71的全长序列：GDV103株基因组全长7404 nt,包括741bp的5’端非编码区（UTR）、6582bp的病毒基因组编码区（ORF）及81bp的3’UTR;安徽株（Anhui2007）基因组全长7405nt,包括742bp的5＇UTR、6582bp的ORF及81bp的3＇UTR;VR1432株基因组全长7408nt,包括743bp的5’UTR、6582bp的ORF及83bp的3’UTR长。经同源性比对和进化树分析,证实GDV103和安徽株EV71属于C基因型的C4基因亚型。而VR1432株则属于C基因型的C2基因亚型。结论：获得了3株EV71的全长基因组序列,并进一步探讨了其型别,为下一步的干扰素保护宴,哈重定了基础．     

Characterization of full-length enterovirus 71 strains from severe and mild disease patients in northeastern China

Human enterovirus 71 (EV71)-associated hand, foot, and mouth disease (HFMD) has been a leading cause of childhood infection in China since 2008. Epidemic and molecular characteristics of HFMD have been examined in many areas of China, including the central and southern regions. However, clinical and genetic characterization of EV71 in the northeastern region of China is scarce. In this study, a series of analyses were performed on seven full-length EV71 sequences from HFMD patients who had either severe or mild disease. We have determined that these seven circulating EV71 viruses from Changchun, China are actually complex recombinant viruses involving multiple type A human enterovirus (HEV). Classified as EV71 subtype C4 (EV71 C4), these Changchun EV71 viruses contain genetic recombination events between the CA4, CA5, EV71B4 and EV71C1 strains. Most of the structural protein region (P1) of these viruses resembled that of the prototype EV71 C1 strains. The non-structural protein domains (P2 and P3) showed a high degree of similarity with CA4, CA5 and EV71 B4 in different regions. The 5'UTR had unclassified recombination,while partial 3D region of these viruses showed a high degree of similarity to CA16. Phylogenetic analysis of full-length or partial sequences of isolates from severe or mild disease patients in Changchun always formed a single cluster in various phylogenetic analyses of different genomic regions, suggesting that all seven strains originated from one single common ancestor. There was no correlation between viral genomic sequence and virulence. Thus, we found that circulating recombinant forms of EV71 are prevalent among HFMD patients in Northeastern China. The existence of a unique cluster of EV71 related viruses in Northeast China has important implications for vaccine development that would address the increasing prevalence of HFMD.     

Molecular analysis of virulent determinants of enterovirus 71

Enterovirus 71 (EV71) is the most important causative agent of hand, foot and mouth disease (HFMD) in children. In most cases, it is a self-limiting illness. However some EV71 infectious cases can develop severe clinical outcomes, such as encephalitis, meningitis, poliomyelitis like paralysis, and even death. To identify the determinants of virulence, the deduced amino acid sequence of polyprotein and nucleotide sequence of 5'-NTR and 3'-NTR in 25 SC-EV71 strains (strains from severe cases) and 31 MC-EV71 strains (strains from mild cases) were analyzed. Results showed four amino acids on two positions (Gly(P710)/Gln(P710)/Arg(P710) and Glu(P729)) on the DE and EF loop of VP1, one (Lys(P930)) on the surface of protease 2A and four nucleotides on three positions (G(P272), U(P488) and A(P700)/U(P700)) in the 5'-NTR region are associated with EV71 virulent phenotype. Predicted secondary structure of RNA using the consensus sequence of 5'-NTR by RNAStructure showed the mutation of nucleotide at position 488 in strain BJ08-Z004-3 (position 491 in prototype strain BrCr) can result in the discrepancy of an additional pair of nucleotides and thus change the stability of the second structure of IRES. Fragment base content analysis showed that in the region 696 to 714 bp at the 5'-NTR, where the A(P700)/U(P700) was located, the nucleotide constitution ratios differed significantly between SC-EV71 and MC-EV71 strains. In conclusion, comparative genomic analysis showed that virulence of EV71 strains are mainly determined by the amino acids on two positions of VP1, one position of protease 2A and the nucleotides on three positions in 5'-NTR.     

Mutations in the non-structural protein region contribute to intra-genotypic evolution of enterovirus 71

Background

Clinical manifestations of enterovirus 71 (EV71) range from herpangina, hand-foot-and-mouth disease (HFMD), to severe neurological complications. Unlike the situation of switching genotypes seen in EV71 outbreaks during 1998–2008 in Taiwan, genotype B5 was responsible for two large outbreaks in 2008 and 2012, respectively. In China, by contrast, EV71 often persists as a single genotype in the population and causes frequent outbreaks. To investigate genetic changes in viral evolution, complete EV71 genome sequences were used to analyze the intra-genotypic evolution pattern in Taiwan, China, and the Netherlands.

Results

Genotype B5 was predominant in Taiwan’s 2008 outbreak and was re-emergent in 2012. EV71 strains from both outbreaks were phylogenetically segregated into two lineages containing fourteen non-synonymous substitutions predominantly in the non-structural protein coding region. In China, genotype C4 was first seen in 1998 and caused the latest large outbreak in 2008. Unlike shifting genotypes in Taiwan, genotype C4 persisted with progressive drift through time. A majority of non-synonymous mutations occurred in residues located in the non-structural coding region, showing annual increases. Interestingly, genotype B1/B2 in the Netherlands showed another stepwise evolution with dramatic EV71 activity increase in 1986. Phylogeny of the VP1 coding region in 1971–1986 exhibited similar lineage turnover with genotype C4 in China; however, phylogeny of the 3D-encoding region indicated separate lineage appearing after 1983, suggesting that the 3D-encoding region of genotype B2 was derived from an unidentified ancestor that contributed to intra-genotypic evolution in the Netherlands.

Conclusions

Unlike VP1 coding sequences long used for phylogenetic study of enteroviruses due to expected host immune escape, our study emphasizes a dominant role of non-synonymous mutations in non-structural protein regions that contribute to (re-)emergent genotypes in continuous stepwise evolution. Dozens of amino acid substitutions, especially in non-structural proteins, were identified via genetic changes driven through intra-genotypic evolution worldwide. These identified substitutions appeared to increase viral fitness in the population, affording valuable insights not only for viral evolution but also for prevention, control, and vaccine against EV71 infection.     

Co-Circulation and Genomic Recombination of Coxsackievirus A16 and Enterovirus 71 during a Large Outbreak of Hand,Foot, and Mouth Disease in Central China

A total of 1844 patients with hand, foot, and mouth disease (HFMD), most of them were children of age 1–3-year-old, in Central China were hospitalized from 2011 to 2012. Among them, 422 were infected with coxsackievirus A16 (CVA16), 334 were infected with enterovirus 71 (EV71), 38 were co-infected with EV71 and CVA16, and 35 were infected with other enteroviruses. Molecular epidemiology analysis revealed that EV71 and CVA16 were detected year-round, but EV71 circulated mainly in July and CVA16 circulated predominantly in November, and incidence of HFMD was reduced in January and February and increased in March. Clinical data showed that hyperglycemia and neurologic complications were significantly higher in EV71-infected patients, while upper respiratory tract infection and C-reactive protein were significantly higher in CVA16-associated patients. 124 EV71 and 80 CVA16 strains were isolated, among them 56 and 68 EV71 strains were C4a and C4b, while 25 and 55 CVA16 strains were B1a and B1b, respectively. Similarity plots and bootscan analyses based on entire genomic sequences revealed that the three C4a sub-genotype EV71 strains were recombinant with C4b sub-genotype EV71 in 2B–2C region, and the three CVA16 strains were recombinant with EV71 in 2A–2B region. Thus, CVA16 and EV71 were the major causative agents in a large HFMD outbreak in Central China. HFMD incidence was high for children among household contact and was detected year-round, but outbreak was seasonal dependent. CVA16 B1b and EV71 C4b reemerged and caused a large epidemic in China after a quiet period of many years. Moreover, EV71 and CVA16 were co-circulated during the outbreak, which may have contributed to the genomic recombination between the pathogens. It should gain more attention as there may be an upward trend in co-circulation of the two pathogens globally and the new role recombination plays in the emergence of new enterovirus variants.     

Interactome analysis of the EV71 5′ untranslated region in differentiated neuronal cells SH‐SY5Y and regulatory role of FBP3 in viral replication

Enterovirus 71 (EV71), a single‐stranded RNA virus, is one of the most serious neurotropic pathogens in the Asia‐Pacific region. Through interactions with host proteins, the 5′ untranslated region (5′UTR) of EV71 is important for viral replication. To gain a protein profile that interact with the EV71 5′UTR in neuronal cells, we performed a biotinylated RNA‐protein pull‐down assay in conjunction with LC–MS/MS analysis. A total of 109 proteins were detected and subjected to Database for Annotation, Visualization and Integrated Discovery (DAVID) analyses. These proteins were found to be highly correlated with biological processes including RNA processing/splicing, epidermal cell differentiation, and protein folding. A protein–protein interaction network was constructed using the STRING online database to illustrate the interactions of those proteins that are mainly involved in RNA processing/splicing or protein folding. Moreover, we confirmed that the far‐upstream element binding protein 3 (FBP3) was able to bind to the EV71 5′UTR. The redistribution of FBP3 in subcellular compartments was observed after EV71 infection, and the decreased expression of FBP3 in host neuronal cells markedly inhibited viral replication. Our results reveal various host proteins that potentially interact with the EV71 5′UTR in neuronal cells, and we found that FBP3 could serve as a positive regulator in host cells.     

Complete nucleotide sequences of all three poliovirus serotype genomes. Implication for genetic relationship, gene function and antigenic determinants

The complete nucleotide sequences of the genomes of the type 2 ( P712 , Ch, 2ab ) and type 3 (Leon 12a1b ) poliovirus vaccine strains were determined. Comparison of the sequences with the previously established genome sequence of type 1 (LS-c, 2ab ) poliovirus vaccine strain revealed that 71% of the nucleotides in the genome RNAs were common, that the 5' and 3' termini of the genomes were highly homologous, and that more than 80% of the nucleotide differences in the coding region occurred in the third letter position of in-phase codons, resulting in a low frequency of amino acid difference. These results strongly suggested that the serotypes of poliovirus derived from a common prototype. A comparison of the amino acid sequences predicted from the genome sequences showed highest variation in the capsid protein region, whereas non-structural proteins are highly conserved. Initiation of polyprotein synthesis occurs in all three strains more than 740 nucleotides downstream from the 5' end. An analysis of the non-coding region suggests that small peptides that could potentially originate from this region are conserved. The amino acid sequences immediately surrounding the cleavage signals, however, show a higher than average degree of variation. The analysis of the amino acid sequences of the capsid protein VP1 of all serotypes has led to the prediction of potential antigenic sites on the virion involved in neutralization.     

PolyC-Binding Protein 1 Interacts with 5′-Untranslated Region of Enterovirus 71 RNA in Membrane-Associated Complex to Facilitate Viral Replication

Enterovirus 71 (EV71) is one causative agent of hand, foot, and mouth disease (HFMD), which may lead to severe neurological disorders and mortality in children. EV71 genome is a positive single-stranded RNA containing a single open reading frame (ORF) flanked by 5′-untranslated region (5′UTR) and 3′UTR. The 5′UTR is fundamentally important for virus replication by interacting with cellular proteins. Here, we revealed that poly(C)-binding protein 1 (PCBP1) specifically binds to the 5′UTR of EV71. Detailed studies indicated that the RNA-binding K-homologous 1 (KH1) domain of PCBP1 is responsible for its binding to the stem-loop I and IV of EV71 5′UTR. Interestingly, we revealed that PCBP1 is distributed in the nucleus and cytoplasm of uninfected cells, but mainly localized in the cytoplasm of EV71-infected cells due to interaction and co-localization with the viral RNA. Furthermore, sub-cellular distribution analysis showed that PCBP1 is located in ER-derived membrane, in where virus replication occurred in the cytoplasm of EV71-infected cells, suggesting PCBP1 is recruited in a membrane-associated replication complex. In addition, we found that the binding of PCBP1 to 5′UTR resulted in enhancing EV71 viral protein expression and virus production so as to facilitate viral replication. Thus, we revealed a novel mechanism in which PCBP1 as a positive regulator involved in regulation of EV71 replication in the host specialized membrane-associated replication complex, which provides an insight into cellular factors involved in EV71 replication.     

Mapping of attenuating sequences of an avirulent poliovirus type 2 strain.

A mouse model for poliomyelitis was used to identify genomic sequences that attenuate neurovirulence of poliovirus strain P2/P712. This type 2 strain is avirulent in primates and mice yet grows as well as virulent strains in cell culture. The approach used was to exchange portions of the genome of the mouse-virulent P2/Lansing strain with the corresponding region from P2/P712 to identify sequences that could attenuate Lansing neurovirulence in mice. A full-length infectious cDNA of P2/P712 was assembled and used to construct recombinants between P2/P712 and P2/Lansing. The results of neurovirulence testing of 11 recombinants indicated that strong attenuating determinants are located in the 5' noncoding region of P2/P712 and a region encoding capsid protein VP1 and 2Apro, 2B, and part of 2C. An attenuating determinant was further localized to between nucleotides 456 and 628 of P2/P712. A third sequence from P2/P712, nucleotides 752 to 2268, encoding VP4, VP2, and part of VP3, was weakly attenuating. The sequence from nucleotide 4454, approximately halfway through the 2C-coding region, to the end of the P2/P712 genome did not contain attenuating determinants. Nucleotide sequence analysis revealed that P2/P712 differs from the type 2 Sabin vaccine strain by only 22 nucleotides. Six differences lead to amino acid changes in the coding region, and four differences are in the 5' noncoding region. These studies show that, like the type 1 and type 3 Sabin vaccine strains, the attenuated type 2 strain P712 contains multiple attenuating sequences, including strongly attenuating sequences in the 5' noncoding region of the genome.     

Genetic analysis of the VP1 region of human enterovirus 71 strains isolated in Fuyang,China, during 2008

Enterovirus 71 (EV71) is a common cause of Hand, foot, and mouth disease (HFMD) and may also cause severe neurological diseases, such as encephalitis and poliomyelitis-like paralysis. To examine the genetic diversity of EV71, we determined and analyzed the complete VP1 sequences (891 nucleotides) from nine EV71 strains isolated in Fuyang, China. We found that nine EV71 strains isolated were over 98% homologous at the nucleotide level and 93%-100% homologous to members of the C4 subgenogroup. At the amino acid level, these Fuyang strains were 99% -100% homologous to one another, 97%-100% homologous to members of the C4 subgenogroup, and the histidine(H) at amino acid position 22 was conserved among the Fuyang strains. The results indicate that Fuyang isolates belong to genotype C4, and an H at position 22 appears to be a marker for the Fuyang strains.     

An attenuated strain of enterovirus 71 belonging to genotype a showed a broad spectrum of antigenicity with attenuated neurovirulence in cynomolgus monkeys

Enterovirus 71 (EV71) is a causative agent of hand, foot, and mouth disease and is also sometimes associated with serious neurological disorders. In this study, we characterized the antigenicity and tissue specificity of an attenuated strain of EV71 [EV71(S1-3')], which belongs to genotype A, in a monkey infection model. Three cynomolgus monkeys were inoculated with EV71(S1-3'), followed by lethal challenge with the parental virulent strain EV71(BrCr-TR) via an intravenous route on day 45 postinoculation of EV71(S1-3'). Monkeys inoculated with EV71(S1-3') showed a mild neurological symptom (tremor) but survived lethal challenge by virulent EV71(BrCr-TR) without exacerbation of the symptom. The immunized monkey sera showed a broad spectrum of neutralizing activity against different genotypes of EV71, including genotypes A, B1, B4, C2, and C4. For the strains examined, the sera showed the highest neutralization activity against the homotype (genotype A) and the lowest neutralization activity against genotype C2. The order of decreasing neutralization activity of sera was as follows: A > B1 > C4 > B4 > C2. To examine the tissue specificity of EV71(S1-3'), two monkeys were intravenously inoculated with EV71(S1-3'), followed by examination of virus distribution in the central nervous system (CNS) and extraneural tissues. In the CNS, EV71(S1-3') was isolated only from the spinal cord. These results indicate that EV71(S1-3') acts as an effective antigen, although this attenuated strain was still neurotropic when inoculated via the intravenous route.     

A mouse-adapted enterovirus 71 strain causes neurological disease in mice after oral infection

A mouse-adapted enterovirus 71 (EV71) strain with increased virulence in mice, MP4, was generated after four serial passages of the parental EV71 strain 4643 in mice. Strain MP4 exhibited a larger plaque size, grew more rapidly, and was more cytotoxic in vitro than strain 4643. Although strains 4643 and MP4 both induced apoptosis of SK-N-SH human neuroblastoma cells, MP4 was more virulent than 4643 in 1-day-old mice (50% lethal doses, 10(2) and 10(4) PFU/mouse, respectively). Strain MP4 (5 x 10(6) PFU/mouse), but not 4643, could orally infect 7-day-old mice, resulting in rear-limb paralysis followed by death 5 to 9 days after inoculation with the virus. Histopathologically, neuronal loss and apoptosis were evident in the spinal cords as well as the brain stems of the infected mice. The limb muscles displayed massive necrosis. There was early and transient virus replication in the intestines, whereas the spinal cord, brain, and muscle became the sites of viral replication during the late phase of the infection. Virus transmission occurred among infected and noninfected cagemates, as demonstrated by the occurrence of seroconversion and the presence of viable viruses in the stool samples of the latter. Protection against EV71 challenge was demonstrated following administration of hyperimmune serum 1 day after inoculation with the virus. Nucleotide sequence analysis of the genome of EV71 strain MP4 revealed four nucleotide changes on the 5' untranslated region, three on the VP2 region, and eight on the 2C region, resulting in one and four amino acid substitutions in the VP2 and 2C proteins, respectively.