Intercity spread of echovirus 6 in shandong province, china: application of environmental surveillance in tracing circulating enteroviruses

Environmental surveillance is an effective approach in investigating circulating enteroviruses and had been conducted in the cities of Jinan and Linyi since February 2008 and April 2010, respectively. This study analyzed 46 sewage samples collected in the two cities in 2011 and found that echovirus 6 (E6) was the predominant serotype, with 134 isolates (65 in Jinan and 69 in Linyi) from 23 (50%) samples. This differs from the 2010 data that found 29 E6 isolates in Jinan and only 3 in Linyi. Phylogenetic analysis of the VP1 coding region showed that all environmental E6 samples from 2008 to 2011 (n = 167) segregated into two lineages and revealed an increase in VP1 gene diversity in 2011, suggesting that the increased number of E6 detections reflects a real epidemic in the two cities. Most Linyi isolates (n = 61, or 88%) in 2011 segregated into sublineage 1a, together with 18 Jinan isolates in 2011. Interestingly, the ancestral VP1 sequence of sublineage 1a inferred using the maximum-likelihood method had 100% identity with the sequence of one environmental isolate from Jinan in August 2010, suggesting an intercity spread from Jinan to Linyi. By Bayesian phylodynamic methods, the most recent common ancestor of Linyi isolates in sublineage 1a dated back to 24 December 2010, revealing that this sublineage was likely imported into Linyi from August to December in 2010. This study demonstrates that environmental surveillance is a sensitive method in tracing transmission pathways of circulating enteroviruses among different regions and reveals that E6-associated aseptic meningitis is an emerging concern in China.     

Environmental Surveillance of Human Enteroviruses in Shandong Province,China, 2008 to 2012: Serotypes,Temporal Fluctuation,and Molecular Epidemiology

Environmental surveillance is an effective approach in investigating the circulation of polioviruses (PVs) and other human enteroviruses (EVs) in the population. The present report describes the results of environmental surveillance conducted in Shandong Province, China, from 2008 to 2012. A total of 129 sewage samples were collected, and 168 PVs and 1,007 nonpolio enteroviruses (NPEVs) were isolated. VP1 sequencing and typing were performed on all isolates. All PV strains were Sabin-like, with the numbers of VP1 substitutions ranging from 0 to 7. The NPEVs belonged to 19 serotypes, and echovirus 6 (E6), E11, coxsackievirus B3 (CVB3), E3, E12, and E7 were the six main serotypes, which accounted for 18.3%, 14.8%, 14.5%, 12.9%, 9.0%, and 5.7% of NPEVs isolated, respectively. Typical summer-fall peaks of NPEV were observed in the monthly distribution of isolation, and an epidemic pattern of annual circulation was revealed for the common serotypes. Phylogenetic analysis was performed on environmental CVB3 and E3 strains with global reference strains and local strains from aseptic meningitis patients. Shandong strains formed distinct clusters, and a close relationship was observed between local environmental and clinical strains. As an EV-specific case surveillance system is absent in China and many other countries, continuous environmental surveillance should be encouraged to investigate the temporal circulation and phylogeny of EVs in the population.     

2010年7～12月肠道病毒B组河南株的血清学分布

本文对河南省2010年7~12月的HFMD监测标本进行了肠道病毒B组的血清型分布研究。来自HFMD病例的阳性病毒分离物进行分子分型方法鉴定并进行VP1完整编码区核苷酸序列测定和分析。所获VP1全序与其他B组各基因型代表株和中国大陆株进行比较并构建系统发生树。共获得14株HEV-B河南株,分为E1、E6 、E11、 E13、 E25、 E30共6个血清型。VP1系统进化分析显示2010年河南HEV-B分离株与原型株亲缘关系较远,E25、E11和E6均与山东株亲缘关系最近,E1和E13均与云南株亲缘关系最近,E30与2008年河南株亲缘关系最近。E6病毒存在两种基因型的共循环现象。     

Elucidation of Echovirus 30's Origin and Transmission during the 2012 Aseptic Meningitis Outbreak in Guangdong,China, through Continuing Environmental Surveillance

An aseptic meningitis outbreak occurred in Luoding City of Guangdong, China, in 2012, and echovirus type 30 (ECHO30) was identified as the major causative pathogen. Environmental surveillance indicated that ECHO30 was detected in the sewage of a neighboring city, Guangzhou, from 2010 to 2012 and also in Luoding City sewage samples (6/43, 14%) collected after the outbreak. In order to track the potential origin of the outbreak viral strains, we sequenced the VP1 genes of 29 viral strains from clinical patients and environmental samples. Sequence alignments and phylogenetic analyses based on VP1 gene sequences revealed that virus strains isolated from the sewage of Guangzhou and Luoding cities matched well the clinical strains from the outbreak, with high nucleotide sequence similarity (98.5% to 100%) and similar cluster distribution. Five ECHO30 clinical strains were clustered with the Guangdong environmental strains but diverged from strains from other regions, suggesting that this subcluster of viruses most likely originated from the circulating virus in Guangdong rather than having been more recently imported from other regions. These findings underscore the importance of long-term, continuous environmental surveillance and genetic analysis to monitor circulating enteroviruses.     

Epidemics and Frequent Recombination within Species in Outbreaks of Human Enterovirus B-Associated Hand,Foot and Mouth Disease in Shandong China in 2010 and 2011

The epidemiology and molecular characteristics of human enterovirus B (HEV-B) associated with hand, foot and mouth disease (HFMD) outbreaks in China are not well known. In the present study, we tested 201 HEV isolates from 233 clinical specimens from patients with severe HFMD during 2010–2011 in Linyi, Shandong, China. Of the 201 isolates, 189 were fully typed and 18 corresponded to HEV-B species (six serotypes CVA9, CVB1, CVB4, Echo 6, Echo 25 and Echo 30) using sensitive semi-nested polymerase chain reaction analysis of VP1 gene sequences. Phylogenetic analysis based on the VP1 region showed that eight E30SD belonged to a novel sub-genogroup D2; E25SD belonged to a novel sub-genogroup D6; E6SD belonged to sub-lineage C6 and five CVB1SD belonged to subgroup 4C; and B4SD belonged sub-lineage D2. The full viral genomes of the CVB1SD, E6SD, E25SD and E30SD isolates were sequenced. Analysis of phylogenetic and similarity plots indicated that E25SD recombined with E25-HN-2, E30FDJS03 and E4AUS250 at noncontiguous P2A–P3D regions, while E30SD, E30FDJ03, E25-HN-2 and E9 DM had shared sequences in discrete regions of P2 and P3. Both E6SD and B1SD shared sequences with E1-HN, B4/GX/10, B5-HN, and A9-Alberta in contiguous regions of most of P2 and P3. Genetic algorithm recombination detection analysis further confirmed the existence of multiple potential recombination points. In conclusion, analysis of the complete genomes of E25SD, E30SD, CVB1SD and E6SD isolated from HFMD patients revealed that they formed novel subgenogroup. Given the prevalence and recombination of these viruses in outbreaks of HFMD, persistent surveillance of HFMD-associated HEV-B pathogens is required to predict potential emerging viruses and related disease outbreaks.     

2010年山东省环境污水中脊髓灰质炎病毒的监测及其基因特征分析

本研究在山东省开展了脊髓灰质炎病毒(Poliovirus,PV)的外环境监测,从济南、临沂两地采集污水标本,浓缩处理后进行病毒分离,对分离到的PV采用中和试验进行血清定型,并对其VP1及3D区进行序列测定,分析其基因突变和重组情况。2010年,共采集污水标本32份,PV阳性10份,阳性率31.3%;分离到18株PV(PV1型3株,PV2型9株,PV3型6株),均为疫苗相关株,VP1完整编码区核苷酸变异数在0~4个之间,在3株PV2型病毒和4株PV3型病毒的基因组中发现重组;对VP1区影响神经毒力的减毒位点分析发现,PV1型病毒中有1株在nt 2 749发生突变(A→G),PV2型病毒中有1株在nt2 908发生A→G突变,3株在nt2 909发生U→C突变,6株PV3型病毒全部在nt2 493发生C→U突变。环境污水中可以分离到PV,其基因重组率和主要减毒位点的回复突变率较高,未发现脊灰野毒株和疫苗衍生株脊灰病毒(Vaccine-derived poliovirus,VDPV)。     

Isolation and Characterization of a Type 2 Vaccine-Derived Poliovirus from Environmental Surveillance in China, 2012

Environmental surveillance of poliovirus on sewage has been conducted in Shandong Province, China since 2008. A type 2 vaccine-derived poliovirus (VDPV) with 7 mutations in VP1 coding region was isolated from the sewage collected in the city of Jinan in December 2012. The complete genome sequencing analysis of this isolate revealed 25 nucleotide substitutions, 7 of which resulted in amino acid alteration. No evidence of recombination with other poliovirus serotypes was observed. The virus did not lose temperature sensitive phenotype at 40°C. An estimation based on the evolution rate of the P1 coding region suggested that evolution time of this strain might be 160–176 days. VP1 sequence analysis revealed that this VDPV strain is of no close relationship with other local type 2 polioviruses (n = 66) from sewage collected between May 2012 and June 2013, suggesting the lack of its circulation in the local population. The person who excreted the virus was not known and no closely related virus was isolated in local population via acute flaccid paralysis surveillance. By far this is the first report of VDPV isolated from sewage in China, and these results underscore the value of environmental surveillance in the polio surveillance system even in countries with high rates of OPV coverage.     

Molecular epidemiology of echoviruses 11 and 13, based on an environmental surveillance conducted in Toyama Prefecture, 2002-2003

Nineteen echovirus 11 (E11) and 12 E13 isolates were isolated from three rivers in Toyama Prefecture, Japan, during an environmental surveillance conducted from April 2002 to March 2003. The nucleotide sequences of E13 isolates were closely related to those from patients with aseptic meningitis, with less than 1.3% divergence in the VP1 region of the viral capsid gene, and belonged to the same clade responsible for a worldwide outbreak that started in 2000. In contrast, E11 isolates were clustered into three genomic groups and were not closely related to echovirus strains isolated from patients. These results suggest that the combination of both virus isolation from environmental sources and phylogenetic analysis could be complementary assessment approaches to trace prevalent and minor circulating enteroviruses in the human population.     

Molecular Epidemiology of Echoviruses 11 and 13, Based on an Environmental Surveillance Conducted in Toyama Prefecture, 2002-2003

Nineteen echovirus 11 (E11) and 12 E13 isolates were isolated from three rivers in Toyama Prefecture, Japan, during an environmental surveillance conducted from April 2002 to March 2003. The nucleotide sequences of E13 isolates were closely related to those from patients with aseptic meningitis, with less than 1.3% divergence in the VP1 region of the viral capsid gene, and belonged to the same clade responsible for a worldwide outbreak that started in 2000. In contrast, E11 isolates were clustered into three genomic groups and were not closely related to echovirus strains isolated from patients. These results suggest that the combination of both virus isolation from environmental sources and phylogenetic analysis could be complementary assessment approaches to trace prevalent and minor circulating enteroviruses in the human population.     

Environmental surveillance system to track wild poliovirus transmission

Eradication of poliomyelitis from large metropolis cities in India has been difficult due to high population density and the presence of large urban slums. Three paralytic poliomyelitis cases were reported in Mumbai, India, in 1999 and 2000 in spite of high immunization coverage and good-quality supplementary immunization activities. We therefore established a systematic environmental surveillance study by weekly screening of sewage samples from three high-risk slum areas to detect the silent transmission of wild poliovirus. In 2001, from among the 137 sewage samples tested, wild poliovirus type 1 was isolated from 35 and wild poliovirus type 3 was isolated from 1. Acute flaccid paralysis (AFP) surveillance indicated one case of paralytic poliomyelitis from the city. Phylogenetic analysis with complete VP1 sequences revealed that the isolates from environmental samples belonged to four lineages of wild polioviruses recently isolated from poliomyelitis cases in Uttar Pradesh and not to those previously isolated from AFP cases in Mumbai. Wild poliovirus thus introduced caused one case of paralytic poliomyelitis. The virus was detected in environmental samples 3 months before. It was found that wild polioviruses introduced several times during the year circulated in Mumbai for a limited period before being eliminated. Environmental surveillance was found to be sensitive for the detection of wild poliovirus silent transmission. Nucleotide sequence analysis helped identify wild poliovirus reservoir areas.     

Environmental Surveillance System To Track Wild Poliovirus Transmission

Eradication of poliomyelitis from large metropolis cities in India has been difficult due to high population density and the presence of large urban slums. Three paralytic poliomyelitis cases were reported in Mumbai, India, in 1999 and 2000 in spite of high immunization coverage and good-quality supplementary immunization activities. We therefore established a systematic environmental surveillance study by weekly screening of sewage samples from three high-risk slum areas to detect the silent transmission of wild poliovirus. In 2001, from among the 137 sewage samples tested, wild poliovirus type 1 was isolated from 35 and wild poliovirus type 3 was isolated from 1. Acute flaccid paralysis (AFP) surveillance indicated one case of paralytic poliomyelitis from the city. Phylogenetic analysis with complete VP1 sequences revealed that the isolates from environmental samples belonged to four lineages of wild polioviruses recently isolated from poliomyelitis cases in Uttar Pradesh and not to those previously isolated from AFP cases in Mumbai. Wild poliovirus thus introduced caused one case of paralytic poliomyelitis. The virus was detected in environmental samples 3 months before. It was found that wild polioviruses introduced several times during the year circulated in Mumbai for a limited period before being eliminated. Environmental surveillance was found to be sensitive for the detection of wild poliovirus silent transmission. Nucleotide sequence analysis helped identify wild poliovirus reservoir areas.     

Environmental Surveillance of Poliovirus in Sewage Water around the Introduction Period for Inactivated Polio Vaccine in Japan

Environmental virus surveillance was conducted at two independent sewage plants from urban and rural areas in the northern prefecture of the Kyushu district, Japan, to trace polioviruses (PVs) within communities. Consequently, 83 PVs were isolated over a 34-month period from April 2010 to January 2013. The frequency of PV isolation at the urban plant was 1.5 times higher than that at the rural plant. Molecular sequence analysis of the viral VP1 gene identified all three serotypes among the PV isolates, with the most prevalent serotype being type 2 (46%). Nearly all poliovirus isolates exhibited more than one nucleotide mutation from the Sabin vaccine strains. During this study, inactivated poliovirus vaccine (IPV) was introduced for routine immunization on 1 September 2012, replacing the live oral poliovirus vaccine (OPV). Interestingly, the frequency of PV isolation from sewage waters declined before OPV cessation at both sites. Our study highlights the importance of environmental surveillance for the detection of the excretion of PVs from an OPV-immunized population in a highly sensitive manner, during the OPV-to-IPV transition period.     

柯萨奇病毒A组2、6、8、12型山东地方株型别鉴定及其基因特征分析

在前期研究中,对山东省1989～2011年急性弛缓性麻痹(Acute flaccid paralysis,AFP)和手足口病(Hand,foot and mouth disease,HFMD)患者标本中分离到的部分人类肠道病毒(Human enteroviruses,HEVs)进行分子定型,发现了8个HEV-A血清型。为进一步探究HEV-A组病毒在山东省的基因型分布及分子进化特征,本研究继续对剩余分离株进行核酸提取,以HEV-A组特异性引物对VP1完整编码区进行序列扩增、测序和分子定型,又鉴定出4个血清型,分别为柯萨奇病毒A组(Coxsackievirus A,CVA)2、6、8、12型,共7株病毒。同源性分析显示,7株HEV-A山东分离株与其原型株核苷酸同源性介于80.8%～85.0%之间。系统进化树表明,CVA8、12型山东分离株与国内分离株存在一定亲缘关系,而CVA2、6型山东株同国内外分离株亲缘关系均较远,并且CVA2、6型在山东省内可能存在多个传播链。本研究将我省HEV-A组的血清型别增加至12个,进一步丰富了HEV-A山东地方株的遗传进化特征,并且发现CVA2、6、8、12血清型在山东省甚至全国范围内存在不同基因特征的传播链的共循环。     

Nucleotide sequence analysis of VP7 gene of Indian isolates of bluetongue virus vis-à-vis other serotypes from different parts of the world.

Bluetongue virus (BTV), a member of genus Orbivirus, a family Reoviridae, is a non-enveloped with double shelled structure and ten segmented double stranded (ds) RNA genome. The RNA segment S7 encodes an inner capsid serogroup specific viral protein VP7. To amplify coding region of VP7 gene of BTV, new primers, forward primer (18-38 bp) and reverse primer (1156-1136 bp), were designed using VP7 gene sequences available in GenBank. This primer pair successfully amplified cell culture adapted Indian isolates of BTV belonging to two different serotypes 1 and 18. The coding sequences of two Indian isolates of BTV (BTV-1H and BTV-18B) were cloned into pPCR Script-Amp SK (+) plasmid vector and transformed into XL10-Gold Kan ultracompetent E. coli cells. The positive clones selected by blue-white screening and colony touch PCR were sequenced. The sequence analysis revealed that there was 93-97% nucleotide sequence identity in VP7 gene of three different Indian serotypes of BTV. The VP7 gene sequences of Indian isolates have comparatively less sequence homology (< 80%) with American (US), and French isolates compared to South African (SA), Australian (AUS) and Chinese (PRC) isolates. In silico restriction enzyme profile analysis of VP7 gene sequences revealed that Indian isolates of BTV-1 can be differentiated from other BTV-1 isolates reported from SA, AUS and PRC using TaqI. Similarly the Indian isolates of BTV belonging to three different serotypes can be differentiated using EcoRI, Hae III and TaqI restriction enzymes.     

Molecular evolution of a type 1 wild-vaccine poliovirus recombinant during widespread circulation in China

Type 1 wild-vaccine recombinant polioviruses were isolated from poliomyelitis patients in China from 1991 to 1993. We compared the sequences of 34 recombinant isolates over the 1,353-nucleotide (nt) genomic interval (nt 2480 to 3832) encoding the major capsid protein, VP1, and the protease, 2A. All recombinants had a 367-nt block of sequence (nt 3271 to 3637) derived from the Sabin 1 oral poliovirus vaccine strain spanning the 3'-terminal sequences of VP1 (115 nt) and the 5' half of 2A (252 nt). The remaining VP1 sequences were closely (up to 99.5%) related to those of a major genotype of wild type 1 poliovirus endemic to China up to 1994. In contrast, the non-vaccine-derived sequences at the 3' half of 2A were more distantly related (<90% nucleotide sequence match) to those of other contemporary wild polioviruses from China. The vaccine-derived sequences of the earliest (April 1991) isolates completely matched those of Sabin 1. Later isolates diverged from the early isolates primarily by accumulation of synonymous base substitutions (at a rate of approximately 3.7 x 10(-2) substitutions per synonymous site per year) over the entire VP1-2A interval. Distinct evolutionary lineages were found in different Chinese provinces. From the combined epidemiologic and evolutionary analyses, we propose that the recombinant virus arose during mixed infection of a single individual in northern China in early 1991 and that its progeny spread by multiple independent chains of transmission into some of the most populous areas of China within a year of the initiating infection.     

Rubella epidemics and genotypic distribution of the rubella virus in shandong province, china, in 1999-2010

Background

The rubella vaccine was introduced into the immunization program in 1995 in the Shandong province, China. A series of different rubella vaccination strategies were implemented at different stages of measles control in Shandong province.

Methodology/Principal Findings

The average reported incidence rate of rubella cases remained at a low level in Shandong province after 1999. However, rubella epidemics occurred repeatedly in 2001/2002, 2006, and 2008/2009. The age of the onset of rubella cases gradually increased during 1999–2010, which showed that most cases were found among the 10 years old in 1999 and among the 17 years old in 2010. Phylogenetic analysis was performed and a phylogenetic tree was constructed based on the World Health Organization standard sequence window for rubella virus isolates. All rubella viruses isolated in Shandong province were divided into 4 genotypes: 1E, 1F, 2A, and 2B. Genotype 1E viruses accounted for the majority (79%) of all these viruses. The similarity of nucleotide and amino acid sequences among genotype 1E viruses was 98.2–100% and 99.1–100%, respectively. All Shandong genotype 1E strains, differed from international genotype 1E strains, belonged to cluster 1 and interdigitated with the viruses from other provinces in mainland China. The effective number of infections indicated by a Bayesian skyline plot remained constant from 2001 to 2009.

Conclusions/Significance

The gradual shift of disease burden to an older age group occurred after a rubella-containing vaccine was introduced into the childhood immunization schedule in 1995 in Shandong province. Four genotypes, including 1E, 1F, 2A, and 2B, were found in Shandong province during 2000–2009. Genotype 1E, rather than genotype 1F, became the predominant genotype circulating in Shandong province from 2001. All Shandong genotype 1E viruses belong to the genotype 1E/cluster 1; they have constantly circulated, and co-evolved and co-circulated, with those from other provinces.     

Genomic surveillance of coxsackievirus A10 reveals genetic features and recent appearance of genogroup D in Shanghai,China, 2016-2020

Coxsackievirus A10 (CVA10) is one of the major causative agents of hand, foot and mouth disease (HFMD). To investigate the epidemiological characteristics as well as genetic features of CVA10 currently circulating in Shanghai, China, we collected a total of 9,952 sporadic HFMD cases from January 2016 to December 2020. In the past five years, CVA10 was the fourth prevalent causatives associated with HFMD in Shanghai and the overall positive rate was 2.78%. The annual distribution experienced significant fluctuations over the past five years. In addition to entire VP1 sequencing, complete genome sequencing and recombination analysis of CVA10 isolates in Shanghai were further performed. A total of 64 near complete genomes and 11 entire VP1 sequences in this study combined with reference sequences publicly available were integrated into phylogenetic analysis. The CVA10 sequences in this study mainly belonged to genogroup C and presented 91%-100% nucleotide identity with other Chinese isolates based on VP1 region. For the first time, our study reported the appearance of CVA10 genogroup D in Chinese mainland, which had led to large-scale outbreaks in Europe previously. The recombination analysis showed the recombination break point located between 5,100 nt and 6,700 nt, which suggesting intertypic recombination with CVA16 genogroup D. To conclusion, CVA10 genogroup C was the predominant genogroup in Shanghai during 2016-2020. CVA10 recombinant genogroup D was firstly reported in circulating in Chinese mainland. Continuous surveillance is needed to better understand the evolution relationships and transmission pathways of CVA10 to help to guide disease control and prevention.     

Population Dynamics and Genetic Diversity of C4 Strains of Human Enterovirus 71 in Mainland China, 1998–2010

Background

Since 1997, several countries within the Asian Pacific region have been affected by one or more massive outbreaks of Hand Foot and Mouth Disease (HFMD). Virus typing experiments revealed that these outbreaks were caused by strains of human enterovirus 71 (EV71) belonging to several different, recently emerged subgenogroups. In mainland China, a different situation was observed. The first outbreak, localized in Shangdong Province, was reported in 2007, and was followed by a wide-spread outbreak in mainland China in 2008. Since then, numbers of reported HFMD cases have been persistently high.

Methodology/Principal Findings

To gain insight in the epidemiological behavior of EV71 in China, we studied genetic diversity and EV71 population dynamics to address whether the increase in number of reported EV71 infections reflects a real increase in viral spread or is just the result of increased awareness and surveillance. We used systematically collected VP1 gene sequences of 257 EV71 strains collected in Guangdong province from 2008 to 2010 as part of HFMD surveillance activities, and supplemented them with 305 GenBank EV71 reference stains collected in China from 1998 to 2010. All isolates from Guangdong Province belonged to subgenogroup C4. Viral population dynamics indicated that the increased reporting of HFMD in China since 2007 reflects a real increase in viral spread and continued replacement of viral lineages through time. Amino acid sequence comparisons revealed substitution of amino acid in residues 22, 145 and 289 through time regularly with the VP1 gene of EV71 strains isolated in mainland China from 1998 to 2010.

Conclusions

EV71 strains isolated in mainland China mainly belonged to subgenogroup C4. There was exponential growth of the EV71 virus population in 2007 and 2008. There was amino acid substitution through time regularly with the VP1 gene which possibly increased viral spread and/or ability of the virus to circulate persistently among the Chinese population.