Non-polio enteroviruses from acute flaccid paralysis surveillance in Shandong Province, China, 1988-2013.

Enteroviruses (EVs) are important human pathogens associated with various clinical syndromes. This study represents an overview of non-polio enteroviruses (NPEVs) isolated from acute flaccid paralysis (AFP) surveillance in Shandong Province, China from 1988 to 2013. Altogether 792 and 170 NPEV isolates were isolated from stool specimens of 9263 AFP cases and 1059 contacts, respectively. Complete VP1 sequencing and typing on all 962 isolates revealed 53 NPEV types in which echovirus (E) 6 (7.6%), E14 (7.6%), E11 (7.4%), coxsackievirus (CV) B3 (7.4%), E25 (5.6%), CVB5 (4.9%), E7 (4.5%) and EV-A71 (4.4%) were the eight most commonly reported serotypes. Distinct summer-fall seasonality was observed, with June-October accounting for 79.3% of isolation from AFP cases with known month of specimen collection. Increase of isolation of EV-A71 and CVA--the predominant pathogens for the hand, foot, and mouth disease--was observed in recent years. Sequence analysis on VP1 coding region of EV-A71 and E6 suggested Shandong strains had great genetic divergence with isolates from other countries. The results described in this study provide valuable information on the circulation and emergence of different EV types in the context of limited EV surveillance in China.

Enteroviruses (EVs) (genus Enterovirus, family Picornaviridae) are common human pathogens. EV infection is usually asymptomatic, but sometimes it is associated with diverse clinical syndromes ranging from minor febrile illness to severe, potentially fatal diseases such as aseptic meningitis, encephalitis, paralysis, myocarditis and neonatal enteroviral sepsis12. EV infections maintain to be an important public health problem. It is estimated that 10–15 million symptomatic enterovirus infections occurred in the United States each year3.

EVs are small, nonenveloped, single-stranded RNA viruses. EV RNA is approximately 7.5 kb long and encodes a polyprotein that is processed to yield the mature structural (VP1 to VP4) and nonstructural proteins1. EV serotypes have traditionally been classified into echoviruses (E), coxsackieviruses (CV) group A and B, and polioviruses (PV) by neutralization test, and molecular typing method based on VP1 sequences has led to the discovery of a large number of new EV types456. So far, human EVs comprises more than 100 types which are classified into 4 species, EV-A to EV-D7.

PV infection is known to be associated with acute paralytic poliomyelitis. The standard approach recommended by World Health Organization (WHO) for polio surveillance is the detection and investigation of acute flaccid paralysis (AFP) cases, which includes standardized virological analysis of faecal samples of the patient, and/or sometimes those from contacts. Systematic AFP surveillance was initiated in China in 1994, and earlier experience on a small scale in Shandong Province could be dated back to 1988. Since 1996, annual reported AFP cases remained stable around 400. Besides PVs, a considerable proportion of non-polio enteroviruses (NPEVs) can be obtained from AFP surveillance as well. Since there is no specific EV surveillance system in mainland China yet, the AFP surveillance is the only data source for understanding EV circulations.

In China, large-scale outbreaks of hand, foot, and mouth disease (HFMD) occurred repeatedly since 2007, and subgenogroup C4 was revealed to be the sole viral genetic lineage circulating in mainland China since then8. However, existence of C2 subgenogroup of EV-A71 in Shandong in 1996 was demonstrated in AFP surveillance9, suggesting retrospective investigation on EVs from AFP surveillance can provide valuable information on EV circulation. A study in Yunnan Province has identified 46 NPEV serotypes from AFP surveillance in 1997–2000 and 200410. However, to the best of our knowledge, no similar study has been reported in other provinces of China yet.

Shandong is a coastal province with a large population (~96 million) and major ports that could potentially serve as portals for importation of exogenous viruses. Here we describe the molecular typing and transmission patterns of all NPEV isolates detected during 26 years surveillance along with the phylogenetic analysis of two pre-dominant types EV-A71 and E6.

Results

Cases and NPEV isolation

The AFP surveillance was conducted in all 138 counties of Shandong Province and more than 600 sentinel hospitals were included. Hospitalized patients < 15 years of age were reported, and male to female ratio was 1.9:1 for all cases from 1990 to 2013. The data for the cases in AFP surveillance in the years 1988 and 1989 were unavailable. Nevertheless, 4 NPEV strains were isolated during that period. During 1990–2013, a total of 9263 AFP cases were reported. EV isolation was performed on all stool specimens and 788 NPEV strains were isolated from stool specimens of these cases. Also, specimens from 1059 contacts of AFP cases were collected during this period and 170 NPEVs were isolated. The annual numbers of reported AFP cases and isolates were illustrated in Fig. 1. AFP surveillance in China increased substantially throughout the 1990s, but the surveillance and report before 1995 was not so active, resulting in the low level of NPEV isolation. It should also be notified that the absence of information of contacts and related isolates before 2007 (Fig. 1) are due to the unavailable information, not indicating that no specimens from contacts were collected at that time.

Figure 1

Annual numbers of AFP cases and contacts and isolation in Shandong Province, China, 1988–2013.

The monthly distributions of AFP cases and NPEV isolation were illustrated in Fig. 2. AFP reports increased in summer and autumn months with a peak in July (1100 [11.9%] of 9255 reports with known month during 1990–2013). However, a more prominent seasonality of NPEV detection from AFP cases was observed with June–October accounting for 79.3% (629 of 792) of isolation with known month of specimen collection.

Figure 2

Numbers of AFP cases and NPEV isolates from AFP cases by month of specimen collection — Shandong, 1990–2013.

Serotype

Molecular typing was performed successfully on all 962 isolates from AFP cases and contacts during 1988–2013, and a total of 53 NPEV types were identified (Table 1, 2). Of these, the 16 most commonly reported serotypes (E6, E14, E11, CVB3, E25, CVB5, E7, EV-A71, E3, E30, CVB1, CVA9, CVB2, CVA24, E13 and CVA4) accounted for 74.6% of all isolation. The eight most common serotypes accounted for approximately half of all reports (49.3%). A total of 19 isolates belonging to 9 new EV types was identified including EV-A76, EV-A90, EV-B73, EV-B74, EV-B75, EV-B80, EV-B87, EV-B97 and EV-C96 (Table 1, 2). Co-infection of different serotypes was observed in 15 patients and 1 contact.

Table 1

Annual isolation of different EV types from AFP cases in Shandong Province, China, 1988–2013

		Year
Species	Type	88	89	90	91	92	93	94	95	96	97	98	99	00	01	02	03	04	05	06	07	08	09	10	11	12	13	Total
A	CVA2	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	1	0	1	0	5	8
	CVA4	0	0	0	0	0	0	0	0	1	0	2	0	0	0	0	1	0	0	1	0	1	3	3	5	0	2	19
	CVA5	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	1
	CVA6	0	0	0	0	1	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	1	0	2	0	0	1	6
	CVA8	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	1
	CVA10	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	1	0	1	0	1	1	3	0	8
	CVA12	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	1
	CVA14	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	1
	CVA16	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	1	0	2
	EV-A71	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	1	1	2	2	2	1	3	1	6	2	10	32
	EV-A76	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	1	0	0	0	0	0	0	0	0	2
	EV-A90	0	0	0	0	0	0	0	0	0	0	0	0	0	2	0	1	0	0	0	0	0	0	0	0	0	0	3
	Total	0	0	0	0	1	0	0	0	4	0	2	0	0	2	0	3	3	3	5	3	4	7	7	16	6	18	84
B	CVA9	0	0	0	1	0	0	0	0	0	4	4	0	2	4	0	0	4	1	4	0	0	0	0	0	1	1	26
	CVB1	0	0	1	0	0	0	1	1	0	0	0	0	3	0	0	0	0	1	6	0	0	1	0	4	1	0	19
	CVB3	0	0	1	0	0	1	1	1	7	0	0	0	9	6	13	1	2	7	0	0	9	1	0	1	4	0	64
	CVB2	0	0	0	0	0	0	1	0	0	2	0	0	2	1	0	2	2	0	0	1	1	3	0	3	0	1	19
	CVB4	0	0	0	0	0	0	0	0	0	1	0	0	3	0	1	3	2	0	0	1	0	0	3	0	0	1	15
	CVB5	0	1	0	0	4	9	0	3	2	0	2	0	0	0	6	1	0	2	1	0	0	8	2	0	1	1	43
	E1	0	0	1	3	0	0	0	0	0	2	0	0	3	1	0	2	2	0	1	0	0	0	1	1	1	0	18
	E2	0	0	0	0	0	0	1	0	1	0	0	0	0	0	0	0	0	0	0	1	0	0	0	1	0	0	4
	E3	0	0	1	0	0	0	0	1	0	1	0	0	1	0	0	5	4	0	2	1	2	4	1	1	5	0	29
	E4	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1
	E6	1	0	0	1	0	0	0	3	0	4	5	1	8	4	3	0	0	1	1	6	3	1	8	8	0	0	58
	E7	0	0	0	0	0	0	1	0	0	3	2	1	1	0	3	9	1	3	4	2	0	1	0	0	1	1	33
	E9	0	0	0	0	0	0	0	0	0	0	1	1	1	1	0	3	1	2	0	0	0	0	5	0	1	2	18
	E11	0	1	0	4	2	1	5	1	5	5	0	2	7	2	0	4	0	0	3	2	10	3	2	0	3	1	63
	E12	0	0	0	3	0	1	0	2	0	1	8	0	0	0	1	0	1	0	1	1	0	0	0	0	0	0	19
	E13	0	0	0	3	0	0	0	0	1	1	0	0	2	0	0	4	0	0	7	2	0	0	3	0	0	0	23
	E14	0	0	0	0	0	2	0	0	0	2	2	1	1	4	2	6	9	5	5	2	3	1	2	0	2	3	52
	E15	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	1
	E16	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1
	E17	0	0	0	0	0	0	0	1	1	0	0	0	0	1	1	0	0	0	0	0	0	0	0	0	0	0	4
	E18	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	2	0	3
	E19	0	0	0	0	1	0	1	0	0	1	0	0	0	1	3	5	0	1	0	0	0	3	0	0	0	1	17
	E20	0	0	0	2	0	0	1	0	0	0	0	2	0	5	0	1	0	0	0	0	0	0	0	0	0	0	11
	E21	0	0	0	2	0	0	0	1	0	1	0	0	0	1	1	2	0	0	2	0	0	0	0	0	4	0	14
	E24	0	0	0	0	0	0	0	0	0	0	1	1	0	0	0	2	0	0	2	0	0	0	0	0	0	0	6
	E25	0	0	0	0	0	0	2	1	0	2	3	0	0	3	1	7	1	4	6	3	2	2	0	4	0	0	41
	E29	0	0	0	3	0	0	0	0	1	5	6	0	1	0	0	0	1	0	0	0	4	0	0	0	0	0	21
	E30	0	0	0	1	0	1	0	0	0	0	1	1	0	1	0	6	3	2	4	0	4	3	0	1	3	0	31
	E31	0	0	0	0	0	0	0	0	0	2	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	2
	E33	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1
	EV-B73	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	1
	EV-B74	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	1
	EV-B75	0	0	0	0	0	0	0	0	0	2	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	2
	EV-B80	0	0	0	0	0	0	0	0	0	1	2	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	4
	EV-B87	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	1
	EV-B97	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1
	Total	1	2	4	23	7	15	14	17	19	40	37	11	45	35	35	63	36	31	49	22	38	31	27	24	29	12	667
C	CVA17	0	0	2	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	3
	CVA20	0	0	0	2	0	0	1	2	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	6
	CVA21	0	0	0	0	0	0	0	0	0	6	0	2	0	1	0	1	0	0	0	0	0	0	0	0	0	0	10
	CVA24	1	0	0	0	0	1	3	1	0	1	0	1	0	0	0	1	0	0	0	1	2	1	4	1	1	1	20
	EV-C96	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0	0	2
	Total	1	0	2	3	0	1	4	3	0	8	0	3	0	1	0	2	0	1	0	1	2	1	4	1	1	2	41
Total		2	2	6	26	8	16	18	20	23	48	39	14	45	38	35	68	39	35	54	26	44	39	38	41	36	32	792

Table 2

Annual isolation of different EV types from contacts of AFP cases in Shandong Province, China, 1988–2013

		Year
Species	Type	93	01	05	07	08	09	10	11	12	13	Total
A	CVA2	0	0	0	0	0	0	0	1	0	0	1
	CVA4	0	0	0	0	1	0	1	2	0	0	4
	CVA10	0	0	0	0	0	0	0	2	0	0	2
	EV-A71	0	0	0	0	0	1	2	7	0	0	10
	Total	0	0	0	0	1	1	3	12	0	0	17
B	CVA9	0	0	0	0	0	0	2	1	1	0	4
	CVB1	0	0	0	1	0	0	2	13	0	0	16
	CVB3	0	0	0	0	4	0	0	3	0	0	7
	CVB2	0	0	0	0	0	1	3	2	0	2	8
	CVB4	0	0	0	2	0	0	3	0	0	0	5
	CVB5	2	0	0	0	0	2	0	0	0	0	4
	E1	0	0	0	0	0	0	0	1	0	0	1
	E3	0	0	0	6	0	0	1	1	1	4	13
	E6	1	1	0	2	0	0	5	6	0	0	15
	E7	2	0	0	0	5	0	1	0	0	2	10
	E9	0	0	0	0	0	0	4	0	0	0	4
	E11	0	0	1	3	0	1	2	0	0	1	8
	E12	2	0	0	0	0	0	0	1	0	0	3
	E14	1	0	0	1	1	4	6	4	4	0	21
	E15	0	0	0	1	0	0	0	1	0	0	2
	E16	0	0	0	0	0	0	1	0	0	0	1
	E20	0	1	0	0	0	0	0	0	0	0	1
	E21	0	0	0	0	0	0	0	0	1	0	1
	E24	0	0	0	1	0	0	0	0	0	0	1
	E25	0	0	0	1	3	3	2	4	0	0	13
	E30	0	0	0	0	2	1	3	1	0	0	7
	Total	8	2	1	18	15	12	35	38	7	9	145
C	CVA24	0	0	0	1	2	1	1	1	0	0	6
	EV-C96	0	0	0	0	0	2	0	0	0	0	2
	Total	0	0	0	1	2	3	1	1	0	0	8
Total		8	2	1	19	18	16	39	51	7	9	170

During 1988–2013, NPEVs belonging to species A, B and C accounted for 10.5% (101 of 962), 84.4% (812 of 962) and 5.1% (49 of 962) of total isolation, respectively. No EV-D virus was identified. EV-B is the major species in most years of the study period. However, dramatic increase in EV-A constitution in AFP cases was observed in recent years (Table 1). Especially in 2011 and 2013, EV-A accounted for 39.0% (16 of 41) and 56.3% (18 of 32) of annual isolation, respectively. Common HFMD pathogens, such as EV-A71, CVA2, CVA4, CVA6, CVA10 and CVA16 were all present in the 26-y surveillance. Interestingly, several EV-A viruses had been observed in AFP surveillance before the onset of national HFMD epidemic in China in 2007. The first isolation of EV-A71, CVA2, CVA4, CVA6 and CVA10 can be dated back to the year 1996, 1996, 1996, 1992 and 2004, respectively.

The most commonly isolated serotypes changed over time. Long-term trends of circulation of eight predominant NPEV serotypes were illustrated in Fig. 3, and epidemic pattern of circulation was observed for all serotypes including EV-A71, which remained rare before 2010 whereas increased dramatically since 2011. In 2011 and 2013, EV-A71 accounted for 14.6% (6 of 41) and 31.2% (10 of 32) of annual isolation, respectively. Epidemic pattern of circulation of E6, E11 and CVB5 was most distinct with regular sharp increases every 3–6 years, which usually last for 1–3 years (Fig. 3). No geographical aggregation was observed for these NPEV types.

Figure 3

Annual number of isolation of 8 predominant serotypes from AFP cases, Shandong, 1990–2013.

Annual isolation of E6, E11, E14 and EV-A71 is illustrated in Fig. 3A, and that of CVB3, CVB5, E7 and E25 is illustrated in Fig. 3B.

VP1 sequence analysis of EV-A71

Altogether 42 EV-A71 isolates, 32 from AFP patients and the rest 10 from contacts, were identified during the AFP surveillance in 1988–2013. During the period from 1988 to 2002, no EV-A71 strains were isolated except for one C2 subgenogroup strain in 19969. Since 2003, continuous detection of EV-A71 was observed and they all belonged to C4a subgenogroup. Homologous comparison revealed 91.8% to 99.8% VP1 nucleotide sequence similarities among C4 subgenogroup strains from AFP surveillance, and 81.4% to 83.6% with prototype strain BrCr. Phylogenetic analysis was performed on the complete VP1 sequences of the isolates from AFP surveillance with global reference sequences and those previously obtained from patients with HFMD and aseptic meningitis in Shandong Province11 (Fig. 4). Except for the C2 subgenogroup strain in 1996, Shandong AFP strains were closely related to local HFMD and aseptic meningitis isolates from 2007 to 2010 and belonged to C4a subgenogroup.

Figure 4

Phylogenetic tree based on EV71 VP1 sequences of Shandong isolates and reference strains of genogroup C.

Circle indicates the prototype strain, triangles indicate strains from AFP surveillance, squares indicates Shandong HFMD reference strains from 2007 to 2010, and diamonds indicate Shandong strains from aseptic meningitis patients. The years of isolation are shown as the first two digits in the names of the strains for AFP isolates, and as the last two digits in the names of the strains for strains from HFMD and aseptic meningitis.

VP1 sequence analysis of E6

A total of 73 E6 strains were isolated during the AFP surveillance in 1988–2013. Homologous comparison revealed 77.3% to 100% VP1 nucleotide sequence similarities among themselves and the overall mean p-distance value is 0.162. Phylogenetic analysis was conducted on the complete VP1 sequences of the isolates with global reference E6 sequences. Shandong E6 strains segregated into 4 clusters, A–D (Fig. 5). Cluster A consisted of isolates from 2007 to 2012, indicating their circulation in Shandong is recently established. Nine strains from other provinces of China also belonged to this cluster. Cluster B consisted of isolates from 1988 to 2002. Members of this cluster disappeared for 11 years, suggesting it might become extinct in Shandong. Cluster C was composed of isolates from 1993 to 2011, indicating this cluster had circulated in Shandong province for a long time. Great genetic divergence was observed with other E6 strains (mean p-distance between groups, 0.204). Cluster D consisted of only 1 strain from 2007. It has close relationship with a strain from aseptic meningitis patient in Anhui Province, China in 2005. They were located in a branch consisted of strains from Australia, France, Russia, and Korea. However, no other Chinese strain of this branch was observed, suggesting this cluster might be imported into China recently and circulated at a low level.

Figure 5

Phylogenetic tree based on E6 VP1 sequences of Shandong isolates and global strains.

Red triangles indicate Shandong strains from AFP surveillance, and green squares indicate strains from other provinces of China. Circle indicates the prototype strain and VP1 sequence of E29 strain JV-10 serves as an outgroup.

Discussion

Since the Global Polio Eradication Initiative (GPEI) launched in 1988, the global incidence of poliomyelitis has dropped a lot. No cases due to wild-type poliovirus 2 (WPV2) have been identified since 199912. However, AFP surveillance is still of great importance in monitoring re-emergence of WPV poliomyelitis in many previously polio-free countries and emergence of vaccine-derived polioviruses (VDPVs)13. During AFP surveillance in Shandong, no WPV has been isolated since 1992, and VPDV had been detected in 2007 (P1, 9 nt; P1, 13 nt), 2009 (P2, 11 nt) and 2011 (P1, 10 nt), respectively. They were all ambiguous VPDV (aVDPV) which were clinical isolates from persons with no known immunodeficiency. In USA, the National Enterovirus Surveillance System (NESS) has provided valuable information for investigating temporal patterns of circulation of different serotypes, guiding outbreak investigations, and identifying targets for development of diagnostic assays and antivirals2. Whereas in China, EV surveillance based on human specimens is very limited and AFP surveillance provided the only information on EV circulation.

Since the correlation of serotype with VP1 sequence was demonstrated5, molecular techniques of EV typing are becoming increasingly available and new EV types continue to be identified1415161718. In the 26-y surveillance, a total of 9 new EV types were identified. The numbers of isolates of each new type ranged from 1 to 4, accounting for a small proportion of total isolation (2.0%, 19/962). According to a previous 5-y study in Yunnan Province10, some of these new EV types such as EV-A76, EV-B75, EV-B80 and EV-C96 were also found. There is a long geographic distance (>1800 km) between the two provinces, and the co-detection of these serotypes suggest they had a wide spread in mainland China, although with low levels of activity. However, the eight most common serotypes in the Yunnan study were E13, E14, E12, CVB3, E2, E6, EV-C96 and E1, which are quite different with the results of this study (E6, E14, E11, CVB3, E25, CVB5, E7 and EV-A71). Considering the results of E11, E13, E7, E6, CVA24, E30, E19 and E29 as the eight most common serotypes in Philippines in a 17-y study19, the variation in constitution of most common serotypes suggests EV serotypes have incongruent epidemic magnitude under different spatio-temporal conditions.

In mainland China, nationwide epidemic of HFMD has occurred repeatedly since 2007. EV-A71 and CVA16 are the two major causative agents in different areas, and other coxsackieviruses of species A, such as CVA2, CVA4, CVA6, CVA10, CVA12, etc, have also been demonstrated to be involved202122. Consistently, an increase of isolation of EV-A71 and other EV-A serotypes from AFP cases was also observed (Table 1, Fig. 3). EV-A71 has been demonstrated to be associated with neurological diseases in young children including aseptic meningitis23, poliomyelitis-like paralysis2425, and fatal brain-stem encephalitis26. The increase in isolation from AFP cases suggests since the large-scale epidemic of HFMD in 2007, EV-A71 associated serious neurological diseases have become a major public health concern in mainland China, which is consistent with the observation of emerging EV-A71 detection in aseptic meningitis patients in a previous study in Shandong Province27.

Shandong Province is located in temperate zone. EV infections are peaked in summer-fall seasons in temperate climates, and the results from our study are consistent with the seasonality (Fig. 2). According to results from EV surveillance in USA in 1970–2005, long-term trends of circulation of different EV serotypes have two major patterns: epidemic and endemic. Epidemic pattern is characterized by substantial fluctuations in circulation levels over time, including large peaks when the serotype was among the most prevalent enteroviruses reported for a given year. Endemic pattern has stable and usually low-levels of circulation with few distinct peaks2. In this study, all eight most common serotypes displayed epidemic pattern of long-term trends, although these data are from AFP surveillance, not specific EV surveillance system. EV-A71 had large-scale nationwide epidemic in current China, and an increase of isolation was observed in AFP cases in this study, whereas no similar trend was observed in the surveillance from USA228. Further surveillance will help understand the long-term trend of EV-A71 circulation in China under the circumstance of continuous epidemics of HFMD. In mainland China, C4 was demonstrated to be the predominant subgenogroup responsible for most documented HFMD outbreaks and epidemic since 19988, and Shandong strains from meningitis patients in 2010 also belonged to subgenogroup C4a11. VP1 sequence analysis in this study revealed close genetic relationship among EV-A71 C4 subgenogroup strains from AFP, HFMD and aseptic meningitis patients, suggesting C4 is an important pathogen of multiple clinical manifestations.

E6 is one of the most frequently isolated serotype in many regions of the world. Its infection includes aseptic meningitis, meningoencephalitis, rashes, gastrointestinal illnesses, neonatal hepatitis and pneumonitis. In China, outbreaks and sporadic aseptic meningitis cases caused by E6 have been reported recently2729. In the VP1 phylogenetic tree, Shandong strains were grouped into four clusters. Except that the only member of cluster D had close relationship with strains from France, Australia, Russia and South Korea, the rest 3 clusters contained Chinese strains exclusively, suggesting the existence of geographic barrier between mainland China and other regions of the world. Frequent travel might increase importation or exportation of different E6 lineages.

Co-infection was detected in 15 patients and 1 contact. It was identified in two kinds of situations. First, in this study, entire VP1 sequence of an EV strain was obtained via combination of 5′ end and 3′ end sequences. So, if molecular typing on 5′ end and 3′ end PCR sequences revealed different types, this indicates a mixture of different EV types in one isolate. Second, RD and HEp-2 cell lines have been demonstrated to possess different sensitivity to different serotypes. RD is sensitive to echovirus and EV-A CVAs, while HEp-2 is sensitive to CVBs and some EV-C types. So, isolation of different EV types from a single specimen was observed from different cell lines. In comprehending the surveillance results of this study, the limitation of virus isolation should be considered. Some serotypes, such as some coxsackie A viruses of EV-C species, cannot produce visible CPE in these cell lines. So, these viruses cannot be recovered in the cell culture method. Hence, the actual EV distribution in the population might not be completely identical with that reflected by this study. Also, this study describes EV typing results from AFP surveillance. Since some NPEVs, such as EV-A71, had been demonstrated to be the pathogens of AFP, the results of this study might be different from usual EV surveillance (infections often milder) and do not necessarily directly reflect to general EV circulation.

China has a vast territory and a large population. The EV circulation in Shandong Province is not necessarily consistent with that of other regions in China. However, EV has the ability of transmission over long distance, and frequent transmissions over cities or provinces in China have been observed previously830. So, the results of serotype spectrum, long-term pattern of circulation and molecular epidemiology can serve as a reflection of the situation in other parts of Eastern China. Further similar studies from other regions will help understand the EV circulation in the country.

Methods

Ethics statement

Ethical approval was given by the Ethics Review Committee of Shandong Center for Disease Control and Prevention, and the methods were carried out in accordance with the principles of the Declaration of Helsinki. Written informed consents for the use of their clinical samples were obtained from all subjects (the legal guardians of the patients and contacts).

Virus isolation

Stool samples from the AFP patients and contacts were collected and processed according to standard procedures recommended by the WHO31. Contacts were collected for AFP cases with (1) <5 years of age, (2) <3 doses of OPV immunization or unknown OPV history, (3) unqualified stool specimen, or (4) clinical diagnosis of poliomyelitis. RD (human rhabdomyosarcoma cell line), HEp-2 (human epidermoid carcinoma cell line) and L20B (mouse cell line expressing the gene for the human cellular receptor for poliovirus) cell lines were used for virus isolation31. All cell lines were gifts from the WHO Global Poliovirus Specialized Laboratory in the USA and were all originally purchased from the American Type Culture Collection. A total of 200 μl of chloroform-treated stool solution was added to each of the cell culture tubes, and the inoculated cells were examined daily. After 7 d, the tubes were frozen, thawed, and re-passaged, and another 7-d examination was performed. Infected cell cultures were harvested and used for further examination until complete cytopathic effect was obtained. To avoid cross contamination, cell tubes of normal L20B, RD and HEp-2 cells served as negative controls. Isolates from RD or HEp-2 cell lines were re-passaged to L20B cell line, and were designated as NPEV if no CPE was observed.

VP1 amplification, sequencing and molecular typing

Total RNA was extracted from 140 μl of the infected cell culture using QIAamp viral RNA mini kit (Qiagen, Valencia, CA, USA) and reverse transcription-PCR (RT-PCR) was performed using Access RT-PCR System (Promega, USA) according to the manufacturers' procedures.

Primer pairs 486/488, 487/489 and 040/011 for EV-A strains, and primer pairs 008/013, 490/492, 187/011, 491/493 and 012/011 for EV-B strains were used for amplification of partial VP1 sequences4532. As for EV-C VP1 amplification, previously designed primer pairs 494/496, 495/497 and 040/011432 are not always efficient. So, we designed primer pairs UF1/UR1, UF2/UR2 and 040/DR1 (Table 3) in this study based on alignment of all available EV-C complete genome sequences from GenBank. These primers turned out to be efficient not only for all EV-C strains described in this study (CVA17, CVA20, CVA21, CVA24 and EV-C96), but also for Sabin PV1–3 and other available CVA13 and EV-C99 strains (from healthy population) of EV-C in our lab. Nevertheless, due to limited EVs in our lab, their performance on other EV-C viruses beyond the listed 10 serotypes had not been testified. Combination of 5′ end and 3′ end sequences yielded entire VP1 coding sequences. As the serotype and species were unknown before typing, EV-B primers were first used for RT-PCR amplification. And EV-A and EV-C primers were used in case of negative results by using EV-B primers. In order to avoid cross contamination, a RT-PCR reaction using the RNA extracted from normal RD cell served as a blank control, and a negative control containing all the components of the reaction except for the template was also included.

Table 3

Primers used for PCR amplification of the VP1 region of EV-C strains

Primer	Sequence (5′ to 3′)	Gene	Position*
UF1	GATGAYTWYACIGMIGGIGG	VP3	2309–2328
UF2	ATGKIYATIYTIGSITTYGT	VP3	2387–2406
UR1	CCRTCITARAARTGISWRTAIGC	VP1	3111–3089
UR2	YTGCCAIGTRTARTCRTCCC	VP1	3007–2988
0405	ATGTAYRTICCIMCIGGIGC	VP1	2951–2970
DR1	GCYWTRTTYTGRTGICCRAA	2A	3408–3389

*Nucleotide sequence coordinates are given relative to the genome sequence of PV1-Sabin (AY184219) for orientation only.

PCR positive products were purified and sequenced bi-directionally with the BigDye Terminator v3.0 Cycle Sequencing kit (Applied Biosystems, Foster City, CA), and sequences were analyzed by ABI 3130 genetic analyzer (Applied Biosystems, Hitachi, Japan). Molecular typing based on VP1 sequences was performed using online Enterovirus Genotyping Tool version 0.133.

Homologous comparison and phylogenetic analysis

Sequence analysis was performed on two predominant serotypes, EV-A71 and E6. The entire list of VP1 sequences of global reference E6 and EV-A71 strains used in sequence analysis can be found as Supplementary Table S1 online. Nucleotide sequence alignments were carried out by BioEdit 7.0.5.3 software34. Phylogenetic trees were constructed by Mega 4.0 using neighbor-joining method after estimation of genetic distance using the Kimura two-parameter method35. A bootstrapping test was performed with 1,000 duplicates.

Author Contributions

Z.T., H.W., Y.L. (Yao Liu), Y.L. (Yan Li), A.X. and L.S. conceived the study and drafted the paper, P.J., X.L., M.L., S.W., Y.Z., Y.F., Q.F., H.Y., J.Y., P.C. and W.L. gathered and analyzed the data, and G.L., F.J, L.F. and P.X. helped to interpret results and contributed to the writing. All authors reviewed the manuscript.

Additional Information

Nucleotide accession numbers. VP1 sequences used in this study were deposited in GenBank under accession numbers GQ253420–GQ253423, GQ329778–GQ329785, JQ326298–JQ326306, HQ829954–HQ829961, HQ399492–HQ399495 and KJ772420–KJ772493.