key: cord-0017510-4d5nyxe5
authors: Chen, Mengyuan; He, Shuizhen; Yan, Qiang; Zhang, Jianmei; Li, Caiyu; Su, Xiaosong; Zhang, Shiyin; Li, Tingdong; Ge, Shengxiang; Chen, Min; Zhang, Jun; Xia, Ningshao
title: Sporadic hand, foot, and mouth disease cases associated with non-C4 enterovirus 71 strains in Xiamen, China, from 2009 to 2018
date: 2021-05-18
journal: Arch Virol
DOI: 10.1007/s00705-021-05097-9
sha: 752c8f8473e4c0ba8a099b29106078f99b9325a2
doc_id: 17510
cord_uid: 4d5nyxe5

Enterovirus 71 (EV71) has caused large hand, foot, and mouth disease (HFMD) epidemics among young children, and EV71 infection is the leading cause of severe HFMD cases and deaths. In mainland China, the prevalence and risk factors of non-C4 EV71 strains are still unclear. In this study, we monitored non-C4 strains over a 10-year HFMD epidemiological surveillance period in Xiamen. The 5’UTR and VP1 coding region of EV71 strains were amplified by RT-nested PCR and sequenced. Thirty-two non-C4 EV71 strains were identified during 2009-2018. This study provides important information about the prevalence of EV71 in China that will be applicable for development of vaccines and diagnostic reagents as well as establishment of policies for HFMD prevention and control. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s00705-021-05097-9.

Hand, foot and mouth disease (HFMD), caused by human enteroviruses (HEVs), is a common pediatric infectious disease. In recent years, mainland China has experienced a large number of HFMD outbreaks [1, 2] . From January 2008 to June 2019, 22.54 million HFMD cases, resulting in 3,704 deaths, were reported to the national surveillance system for HFMD in China established by the Chinese Center for Disease Control and Prevention (Chinese CDC) (http:// www. nhc. gov. cn/ jkj/ s2907/ new_ list. shtml? tdsou rcetag= s_ pcqq_ aiomsg).

The main pathogens causing HFMD in mainland China are enterovirus 71 (EV71), coxsackievirus A16 (CA16), coxsackievirus A6 (CA6), and coxsackievirus A10 (CA10). Because most of the severe cases involving nervous system syndromes and cardiopulmonary failure are caused by EV71, this virus has received the most attention by clinicians and researchers in many countries. EV71 can be clustered into three genotypes: A, B (subgenotypes B1-B5), and C (subgenotypes C1-C5). In other parts of Asia, several subgenotypes of EV71 have been found to circulate or co-circulate in a same area [3] [4] [5] [6] [7] [8] [9] , whereas in mainland China, almost all EV71 strains belong to subgenotype C4 [2, 10] , and only individual sporadic cases caused by non-C4 strains have been reported [11] [12] [13] . The prevalence and risk factors of non-C4 strains are not known. This might lead to bias in vaccine design, development of diagnostic reagents, and disease control and prevention strategies. In this study, we investigated the prevalence of non-C4 strains in a 10-year epidemic study in Xiamen, a transport hub between Southeast Asia and mainland China that may stand in the forefront of epidemic trends.

From January 2008 to September 2019, 64,309 cases, including 124 severe cases, were reported to the Xiamen CDC. During this period, there was a rising trend in HFMD cases in Xiamen ( Supplementary Fig. S1 Samples were examined using reverse transcription nested-PCR (RT-nested PCR) as reported previously [14] . Samples from 5,315 of these cases were positive for HEVs (73.8%). CA6 (33.63%), EV71 (28.45%), CA16 (23.62%), and CA10 (4.98%) were the primary pathogens detected (Fig. 1 ). Each EV71 strain was subgenotyped. The vast majority belonged to subgenotype C4. However, as VP1 analysis could only distinguish between subgenotypes C4 and C5, a small portion of strains failed to be subgenotyped.

In this study, 55 EV71 samples that could not be subgenotyped previously were analyzed. Nucleic acids were extracted using a GenMagBio Viral DNA/RNA Extraction Kit (GenMagBio, China), and 3-step RT-nested-PCR was performed using non-EV71 VP1 primers (Supplementary Table S1 ). The RT mixture in each tube consisted of 4 μl of M-MLV buffer, 4 μl of a dNTP mixture (10 mM each), 0.6 μl of specific primer (10 μM), 0.2 μl of M-MLV Reverse Transcriptase (TaKaRa), 0.2 μl of Recombinant RNasin Ribonuclease Inhibitor (Promega), 5 μl of viral RNA, and 6 μl of nuclease-free water. The RT reaction was performed at 50 ℃ for 30 min. The first-round PCR mixture in each tube consisted of 2 μl of 10× PCR buffer, 2.4 μl of dNTP mixture (2.5 mM each), 0.4 μl of forward primer (10 μM), 0.4 μl of reverse primer (10 μM), 0.2 μl of rTaq (TaKaRa), 5 μl of cDNA, and 9.6 μl of nuclease-free water. The cycling conditions were as follows: initial denaturation at 95 °C for 10 min, followed by 35 cycles of 94 °C for 40 s, 53 °C for 40 s, and 72 °C for 1 min, and final extension at 72 °C for 10 min. A volume of 2 μl of the first PCR product was used for the second-round PCR with 5 μl of 10× PCR buffer, 6 μl of dNTP mixture (2.5 mM each), 1 μl of forward primer (10 μM), 1 μl of reverse primer (10 μM), 0.5 μl of rTaq (TaKaRa), and 34.5 μl of nuclease-free water and using the same PCR conditions. The amplicons were purified and sequenced. All of the sequences were submitted to the Gen-Bank sequence database (NCBI, USA) with the accession numbers MT861057-MT861109. Phylogenetic trees based on partial 5' UTR and VP1 sequences were constructed using the Jukes-Cantor algorithm and the neighbor-joining method in MEGA 7.0 software (www. megas oftwa re. net). As the partial VP1 sequences were from at different regions of the gene, three separate phylogenetic trees based on partial VP1 sequences were constructed.

The results indicated that these sequences belonged to non-C4 strains that are particularly uncommon in mainland China. Thirty-two non-C4 strains were identified in Xiamen (Supplementary Table S2 ). These cases, all of which were mild, showed clinical manifestations similar to those in patients infected with other HEVs. Based on their VP1 sequences, three B3 strains, 10 B5 strains, one C2 strain and one C5 strain were identified. Among the 15 strains, the 5' UTR and VP1 sequences of four strains clustered into different subgenotypes. This suggested that genotype recombination events had occurred in the region between the 5' UTR and the VP1 gene. Based on the VP1 sequences, most of the B5 strains from Xiamen were closely related to strains detected in Taiwan, while two other Xiamen B5 strains had a close relationship to strains from Southeast Asia. All Xiamen B5 strains isolated from 2010 to 2012 were closely grouped into a sublineage. The Xiamen C2 strain was closely related to strains from the Philippines, while the Xiamen C5 strain was related to strains from Vietnam. The Xiamen B3 strains were segregated into two sublineages.

The VP1 sequences of 17 strains failed to be amplified. Phylogenetic analysis based on 5' UTR sequences showed that one strain was genotype A, two were subgenotype B3, and three were subgenotype B5, while 11 showed similarity to orphan strains from India (Supplementary Table S2 and Fig. 2a-d) . Two C4 strains were excluded based on VP1 phylogenetic analysis. In Asia, HFMD cases are frequently associated with EV71. In Southeast Asia, including Taiwan, several EV71 subgenotypes cocirculated in the same region. However, the epidemic patterns in mainland China are entirely different. Except for isolated cases, all of the EV71 strains belong to subgenotype C4. Rare isolations of non-C4 EV71 strains have been reported in different years and areas [11] [12] [13] . We reported previously that the first imported B5 strains in Xiamen (isolated in 2009) displayed a close genetic relationship to strains present in Taiwan and Singapore in 2006-2008 [15] . Five years later, the second imported B5 strain, which shared the highest nucleotide sequence identity with strains found in Vietnam in 2011-2013, was reported in Chongqing [11] . A C2 strain was isolated in Shandong in 1998, and another was isolated in Beijing in 2015 [12, 13] .

For high sensitivity, reagents used in mainland China are usually developed specifically for detection of EV71 C4 strains. The prevalence and risk factors of non-C4 strains are not clear. Furthermore, it is not known whether EV71 vaccines designed based on C4 strains provide cross-protection against non-C4 EV71 strains. As very few researchers have systematically studied the prevalence of non-C4 EV71 strains in mainland China, the deep surveillance of non-C4 EV71 strains described here will provide important information.

In this study, we first systematically investigated non-C4 EV71 cases in the same region of mainland China.

Thirty-two non-C4 EV71 strains, representing subgenotypes B3, B5, C2, and C3 as well as some strains resembling Indian orphan strains were identified. The non-C4 strains accounted for 2.12% of EV71 strains (32/1512, far greater than in other districts). The non-C4 strains had a close relationship to Southeast Asian strains. There might be two reasons for this. On the one hand, a deep, comprehensive, and long-term surveillance system for monitoring the pathogen spectrum of HFMD has been operating in Xiamen for 10 years. As all subgenotypes of HEVs could be identified using our system, non-C4 EV71 strains were not overlooked. On the other hand, as Xiamen is an important tourist city and a key transportation hub between mainland China and Southeast Asia, Xiamen may stand at the forefront of EV71 epidemic changes. Epidemic patterns of HEVs in Xiamen might be affected by imported HEV strains from Southeast Asia. A study in Taiwan showed that large epidemics occurred shortly after minor strains were identified a few years ago. This suggests that large epidemics caused by non-C4 EV71 strains can occur at any time, and more attention should be given to these strains.

This study provides new understanding of EV71 epidemic patterns in mainland China and also provides important information for formulating EV71 prevention and control strategies and developing vaccines and diagnostic reagents. 

Hand, foot, and mouth disease in China, 2008-12: an epidemiological study

Comparative epidemiology and virology of fatal and nonfatal cases of hand, foot and mouth disease in mainland China from

Epidemiology of enterovirus 71 infections in Taiwan

Molecular epidemiology of enterovirus 71 infection in the central region of Taiwan from

Reemergence of enterovirus 71 in 2008 in Taiwan: dynamics of genetic and antigenic evolution from 1998 to

The largest outbreak of hand; foot and mouth disease in Singapore in 2008: the role of enterovirus 71 and coxsackievirus A strains

Genetic characterization of enterovirus 71 strains circulating in Vietnam in 2012

Hand foot and mouth disease due to enterovirus 71 in Malaysia

Phylogenetic analysis of enterovirus 71 strains isolated during linked epidemics in Malaysia

Co-circulation and genomic recombination of coxsackievirus A16 and enterovirus 71 during a large outbreak of hand, foot, and mouth disease in Central China

Isolation of an imported subgenotype B5 strain of human enterovirus A71 in Chongqing City, China

Identification of a C2 subgenogroup strain of enterovirus 71 in a retrospective study in Shandong Province, China, from 1990 to 2010

Analysis of an imported subgenotype C2 strain of human enterovirus 71 in Beijing

Specific primer amplification of the VP1 region directed by 5' UTR sequence analysis: enterovirus testing and identification in clinical samples from hand-foot-and-mouth disease patients

Isolation and identification of an EV71-B5 strain