Evolutionary Diversity of Coxsackievirus A6 Causing Severe Hand, Foot, and Mouth Disease - China, 2012-2023.

Introduction: Coxsackievirus A6 (CVA6) has emerged as a significant pathogen responsible for severe cases of hand, foot, and mouth disease (HFMD). This study aims to delineate the demographic characteristics and analyze the viral evolution of severe HFMD associated with CVA6, thereby assisting in its surveillance and management. Methods: In this investigation, 74 strains of CVA6 were isolated from samples collected from severe HFMD cases between 2012 and 2023. The VP1 gene sequences of CVA6 were amplified and analyzed to assess population historical dynamics and evolutionary characteristics using BEAST, DnaSP6, and PopART. Results: A significant portion (94.4%) of severe CVA6-associated HFMD cases (51 out of 54, with 20 lacking age information) were children under 5 years old. Among the 74 CVA6 strains analyzed, 72 belonged to the D3a sub-genotype, while only two strains were D2 sub-genotype. The average genetic distance between VP1 sequences prior to 2015 was 0.027, which increased to 0.051 when compared to sequences post-2015. Historical population dynamics analysis indicated three significant population expansions of severe CVA6-associated HFMD during 2012-2013, 2013-2014, and 2019-2020, resulting in the formation of 65 distinct haplotypes. Consistent with the MCC tree findings, transitioning between regional haplotypes required multiple base substitutions, showcasing an increase in population diversity during the evolutionary process (from 14 haplotypes in 2013 to 55 haplotypes over the subsequent decade). Conclusions: CVA6, associated with severe HFMD, is evolving and presents a risk of outbreak occurrence. Thus, enhanced surveillance of severe HFMD is imperative.

Genomic epidemiology of CVA10 in Guangdong, China, 2013-2021.

Hand, Foot and Mouth Disease (HFMD) is a highly contagious viral illness primarily affecting children globally. A significant epidemiological transition has been noted in mainland China, characterized by a substantial increase in HFMD cases caused by non-Enterovirus A71 (EV-A71) and non-Coxsackievirus A16 (CVA16) enteroviruses (EVs). Our study conducts a retrospective examination of 36,461 EV-positive specimens collected from Guangdong, China, from 2013 to 2021. Epidemiological trends suggest that, following 2013, Coxsackievirus A6 (CVA6) and Coxsackievirus A10 (CVA10) have emerged as the primary etiological agents for HFMD. In stark contrast, the incidence of EV-A71 has sharply declined, nearing extinction after 2018. Notably, cases of CVA10 infection were considerably younger, with a median age of 1.8 years, compared to 2.3 years for those with EV-A71 infections, possibly indicating accumulated EV-A71-specific herd immunity among young children. Through extensive genomic sequencing and analysis, we identified the N136D mutation in the 2 A protein, contributing to a predominant subcluster within genogroup C of CVA10 circulating in Guangdong since 2017. Additionally, a high frequency of recombination events was observed in genogroup F of CVA10, suggesting that the prevalence of this lineage might be underrecognized. The dynamic landscape of EV genotypes, along with their potential to cause outbreaks, underscores the need to broaden surveillance efforts to include a more diverse spectrum of EV genotypes. Moreover, given the shifting dominance of EV genotypes, it may be prudent to re-evaluate and optimize existing vaccination strategies, which are currently focused primarily target EV-A71.

Evolution, recombination and geographic spreading of global Coxsackievirus A6.

BACKGROUND: The increasing incidence of hand, foot, and mouth disease (HFMD) associated with Coxsackievirus A6 (CVA6) has become a very significant public health problem. The aim of this study is to investigate the recombination, geographic transmission, and evolutionary characteristics of the global CVA6. METHODS: From 2019 to 2022, 73 full-length CVA6 sequences were obtained from HFMD patients in China and analyzed in combination with 1032 published whole genome sequences. Based on this dataset, the phylogenetic features, recombinant diversity, Bayesian phylodynamic characteristics, and key amino acid variations in CVA6 were analyzed. RESULTS: The four genotypes of CVA6, A, D, E, and F, are divided into 24 recombinant forms (RFs, RF-A - RF-X) based on differences in the P3 coding region. The eastern China region plays a key role in the dissemination of CVA6 in China. VP1-137 and VP1-138 are located in the DE loop on the surface of the CVA6 VP1 protein, with the former being a highly variable site and the latter having more non-synonymous substitutions. CONCLUSIONS: Based on whole genome sequences, this study contributes to the CVA6 monitoring, early warning, and the pathogenic mechanism by studying recombination diversity, geographical transmission characteristics, and the variation of important amino acid sites.

Serological investigation on the prevalence of poliovirus in Guangdong province: A cross-sectional study.

In 2019, we conducted a cross-sectional study for polio virus seroprevalence in Guangdong province, China. We assessed the positivity rates of poliomyelitis NA and GMT in serum across various demographic groups, and the current findings were compared with pre-switch data from 2014. Using multistage random sampling method, four counties/districts were randomly selected per city, and within each, one general hospital and two township hospitals were chosen. Healthy individuals coming for medical checkups or vaccination were invited. A total of 1318 individual samples were collected and tested. In non-newborn population, age-dependent positivity rates ranged from 77.8% to 100% for PV1 NA and 70.3% to 98.9% for PV3 NA (p < .01). The lowest GMT values for both types (17.03 and 8.46) occurred in the 20 to <30 years age group, while peak GMTs for PV1 and PV3 were observed in 1 to <2 (340.14) and 0 to <1-year (168.90) age groups, respectively. GMTs for PV1 (P = .002) and PV3 (P = .007) in Eastern Guangdong were lower than those in the other three regions. Male participants showed higher GMTs than females (P = .016 and .033, respectively). In newborn population, both males and females showed higher PV1 NA positivity rates and GMTs compared to PV3 (p < .05). Post-switch PV3 NA positivity rates were higher than pre-switch rates (p = .016). GMTs of both PV1 and PV3 were significantly higher post-switch (p < .001). The positivity rates of NAs and GMTs remain high level, which play an important role in resisting poliomyelitis infection. Effect of the converted immunization program was more pronounced than that before.

Dynamic changes in polioviruses identified by environmental surveillance in Guangzhou, 2009-2021.

Polio cases can be missed by acute flaccid paralysis (AFP) case surveillance alone, emphasizing the importance of environmental surveillance (ES). In this study, to investigate the serotype distribution and epidemiological trends of poliovirus (PV), we characterized PV isolated from domestic sewage in Guangzhou City, Guangdong Province, China from 2009 to 2021. A total of 624 sewage samples were collected from the Liede Sewage Treatment Plant, and the positive rates of PV and non-polio enteroviruses were 66.67% (416/624) and 78.37% (489/624), respectively. After sewage sample treatment, each sewage sample was inoculated in six replicate tubes of three cell lines, and 3370 viruses were isolated during the 13-year surveillance period. Among these, 1086 isolates were identified as PV, including type 1 PV (21.36%), type 2 PV (29.19%), and type 3 PV (49.48%). Based on VP1 sequences, 1057 strains were identified as Sabin-like, 21 strains were high-mutant vaccines, and eight strains were vaccine-derived poliovirus (VDPV). The numbers and serotypes of PV isolates in sewage were influenced by the vaccine switch strategy. After type 2 OPV was removed from the trivalent oral PV (OPV) vaccine and a bivalent OPV (bOPV) was adopted in May 2016, the last type 2 PV strain was isolated from sewage, with no detection thereafter. Type 3 PV isolates increased significantly and became the dominant serotype. Before and after the second vaccine switch in January 2020, that is, from the first dose of IPV and second-fourth doses of bOPV to the first two doses of IPV and third-fourth doses of bOPV, there was also a statistical difference in PV positivity rates in sewage samples. Seven type 2 VDPVs and one type 3 VDPV were identified in sewage samples in 2009-2021, and phylogenetic analysis indicated that all VDPVs isolated from ES in Guangdong are newly discovered VDPVs, different from VDPV previously discovered in China, and were classified as ambiguous VDPV. It is noteworthy that no VDPV cases were reported in AFP case surveillance in the same period. In conclusion, continued PV ES in Guangzhou since April 2008 has been a useful supplement to AFP case surveillance, providing an important basis for evaluating the effectiveness of vaccine immunization strategies. ES improves early detection, prevention, and control; accordingly, this strategy can curb the circulation of VDPVs and provide a strong laboratory basis for maintaining a polio-free status.

Prevalence of Non-Polio Enteroviruses in the Sewage of Guangzhou City, China, from 2013 to 2021.

Continuous surveillance of enteroviruses (EVs) in urban domestic sewage can timely reflect the circulation of EVs in the environment and crowds, and play a predictive and early warning role in EV-related diseases. To better understand the long-term epidemiological trends of circulating EVs and EV-related diseases, we conducted a 9-year (2013 to 2021) surveillance study of non-polio EVs (NPEVs) in urban sewage in Guangzhou city, China. After concentrating and isolating the viruses from the sewage samples, NPEVs were detected and molecular typing was performed. Twenty-one different NPEV serotypes were identified. The most isolated EVs were echovirus 11 (E11), followed by coxsackievirus (CV) B5, E6, and CVB3. EV species B prevailed in sewage samples, but variations in the annual frequency of different serotypes were also observed in different seasons, due to spatial and temporal factors. E11 and E6 were detected continuously before 2017, and the number of isolates was relatively stable during the surveillance period. However, after their explosive growth in 2018 and 2019, their numbers suddenly decreased significantly. CVB3 and CVB5 had alternating trends; CVB5 was most frequently detected in 2013 to 2014 and 2017 to 2018, while CVB3 was most frequently detected in 2015 to 2016 and 2020 to 2021. Phylogenetic analysis showed that at least two different transmission chains of CVB3 and CVB5 were prevalent in Guangzhou City. Our results show that in the absence of a comprehensive and systematic EV-related disease surveillance system in China, environmental surveillance is a powerful and effective tool to strengthen and further investigate the invisible transmission of EVs in the population. IMPORTANCE This study surveilled urban sewage samples from north China for 9 years to monitor enteroviruses. Samples were collected, processed, and viral identification and molecular typing were performed. We detected 21 different non-polio enteroviruses (NPEVs) with yearly variations in prevalence and peak seasons. In addition, this study is very important for understanding the epidemiology of EVs during the COVID-19 pandemic, as the detection frequency and serotypes of EVs in sewage changed considerably around 2020. We believe that our study makes a significant contribution to the literature because our results strongly suggest that environmental surveillance is an exceptionally important tool, which can be employed to detect and monitor organisms of public health concern, which would otherwise be missed and under-reported by case-based surveillance systems alone.

Viral infection and transmission in a large, well-traced outbreak caused by the SARS-CoV-2 Delta variant.

The SARS-CoV-2 Delta variant has spread rapidly worldwide. To provide data on its virological profile, we here report the first local transmission of Delta in mainland China. All 167 infections could be traced back to the first index case. Daily sequential PCR testing of quarantined individuals indicated that the viral loads of Delta infections, when they first become PCR-positive, were on average ~1000 times greater compared to lineage A/B infections during the first epidemic wave in China in early 2020, suggesting potentially faster viral replication and greater infectiousness of Delta during early infection. The estimated transmission bottleneck size of the Delta variant was generally narrow, with 1-3 virions in 29 donor-recipient transmission pairs. However, the transmission of minor iSNVs resulted in at least 3 of the 34 substitutions that were identified in the outbreak, highlighting the contribution of intra-host variants to population-level viral diversity during rapid spread.

Molecular surveillance of coxsackievirus A16 in southern China, 2008-2019.

A national surveillance system on hand, foot, and mouth disease (HFMD) was launched in 2008 in China. Since then, millions of HFMD cases have been reported each year, with enterovirus A71 (EV-A71), coxsackievirus A16 (CV-A16), and coxsackievirus A6 (CV-A6) as the major causative pathogens. Long-term surveillance of viral infection rates and genetic changes is essential for understanding the disease epidemiology pattern. Here, we analyzed molecular surveillance data on CV-A16 covering a period of 12 years (2008-2019) in Guangdong, China, one of the regions reporting the largest number of HFMD cases. Full VP1 sequences of 456 strains were determined to examine the genetic diversity and changes in the distribution of CV-A16 variants. Our study revealed an irregular pattern of CV-A16 infections in Guangdong. Different from the cyclic epidemics observed in some Asia-Pacific regions, there was a continuously high CV-A16 infection rate from 2008 to 2014, and after a period of lower epidemic activity in 2015-2017, an upsurge of CV-A16 infection was observed in 2018-2019. Cocirculation of subgenotypes B1a and B1b was observed, but while subgenotype B1a was predominant from 2008 to 2012, it appears to have been replaced by B1b, which has circulated as the predominant subgenotype since 2013. Phylogenetic analysis showed that most of the circulating CV-A16 strains are endemic, with occasional transmission between neighboring regions. The re-emergence of B1a in 2016-2019 in Guangdong was likely the result of introduction(s) from Southeast Asia. These results highlight the importance of continuous molecular surveillance from different areas, which will improve our understanding of the origin of the epidemic and facilitate the development of strategies for HFMD disease control.

Emergence of a non vaccine-cognate enterovirus A71 genotype C1 in mainland China.

BACKGROUND: EV-A71 is a common causative agent of hand foot and mouth disease. In mainland China, EV-A71 subgenotype C4 has been the sole circulating genotype since 2008, and was used in the production of multiple licensed vaccines. Here, we report the first detection EV-A71 C1 strains in China. METHODS: Full genomic sequence were obtained. The origin of the EV-A71 C1 strains were tracked down by Bayesian inferences. Recombination was analyzed using Simplot program. And the antigenicity were tested using the microneutralization test. RESULTS: The C1-GD2019 shared high identity with the C1-like lineage recently identified in Europe and was introduced into Guangdong in 2018-2019. Close genetic relatedness between the C1-GD2019 and Europe C1-like strains were observed except for the 3D-3'UTR region. The late showed high similarity with CVA genomes. Antigenic variance was found. The C1-GD2019 could not be effectively neutralized by EV-A71 C4a neutralizing antibody positive samples. CONCLUSION: This is the first report of EV-A71 subgenotype C1 isolated in China. It is a recombinant strain originating from C1-like strains recently identified in Europe and CVA strains. The different antigenicity between the C1 strains and C4a vaccine strains highlighted the importance on closely monitoring the EV-A71 C1 strains in China.

Genomic epidemiology of coxsackievirus A16 in mainland of China, 2000-18.

Hand, foot, and mouth disease (HFMD), which is a frequently reported and concerning disease worldwide, is a severe burden on societies globally, especially in the countries of East and Southeast Asia. Coxsackievirus A16 (CV-A16) is one of the most important causes of HFMD and a severe threat to human health, especially in children under 5 years of age. To investigate the epidemiological characteristics, spread dynamics, recombinant forms (RFs), and other features of CV-A16, we leveraged the continuous surveillance data of CV-A16-related HFMD cases collected over an 18-year period. With the advent of the EV-A71 vaccine since 2016, which targeted the EV-A71-related HFMD cases, EV-A71-related HFMD cases decreased dramatically, whereas the CV-A16-related HFMD cases showed an upward trend from 2017 to October 2019. The CV-A16 strains observed in this study were genetically related and widely distributed in the mainland of China. Our results show that three clusters (B1a-B1c) existed in the mainland of China and that the cluster of B1b dominates the diffusion of CV-A16 in China. We found that eastern China played a decisive role in seeding the diffusion of CV-A16 in China, with a more complex and variant transmission trend. Although EV-A71 vaccine was launched in China in 2016, it did not affect the genetic diversity of CV-A16, and its genetic diversity did not decline, which confirmed the epidemiological surveillance trend of CV-A16. Two discontinuous clusters (2000-13 and 2014-18) were observed in the full-length genome and arranged along the time gradient, which revealed the reason why the relative genetic diversity of CV-A16 increased and experienced more complex fluctuation model after 2014. In addition, the switch from RFs B (RF-B) and RF-C co-circulation to RF-D contributes to the prevalence of B1b cluster in China after 2008. The correlation between genotype and RFs partially explained the current prevalence of B1b. This study provides unprecedented full-length genomic sequences of CV-A16 in China, with a wider geographic distribution and a long-term time scale. The study presents valuable information about CV-A16, aimed at developing effective control strategies, as well as a call for a more robust surveillance system, especially in the Asia-Pacific region.

Molecular epidemiology of enteroviruses associated with severe hand, foot and mouth disease in Shenzhen, China, 2014-2018.

In this study, we investigated the epidemiology and molecular characteristics of enteroviruses associated with severe hand, foot and mouth disease (HFMD) in Shenzhen, China, during 2014-2018. A total of 137 fecal specimens from patients with severe HFMD were collected. Enterovirus (EV) types were determined using real-time reverse transcription polymerase chain reaction (RT-PCR), RT nested PCR, and sequencing. Sequences were analyzed using bioinformatics programs. Of 137 specimens tested, 97 (70.8%), 12 (8.8%), and 10 (7.3%) were positive for EV-A71, coxsackievirus A6 (CVA6), and CVA16, respectively. Other pathogens detected included CVA2 (2.9%, 4/137), CVA10 (2.9%, 4/137), CVA5 (0.7%, 1/137), echovirus 6 (E6) (0.7%, 1/137) and E18 (0.7%, 1/137). The most frequent complication in patients with proven EV infections was myoclonic jerk, followed by aseptic encephalitis, tachypnea, and vomiting. The frequencies of vomiting and abnormal eye movements were higher in EV-A71-infected patients than that in CVA6-infected or CVA16-infected patients. Molecular phylogeny based on the complete VP1 gene revealed no association between the subgenotype of the virus and disease severity. Nevertheless, 12 significant mutations that were likely to be associated with virulence or the clinical phenotype were observed in the 5'UTR, 2Apro, 2C, 3A, 3Dpol and 3'UTR of CVA6. Eight significant mutations were observed in the 5'UTR, 2B, 3A, 3Dpol and 3'UTR of CVA16, and 10 significant mutations were observed in the 5'UTR, VP1, 3A and 3Cpro of CVA10. In conclusion, EV-A71 is still the main pathogen causing severe HFMD, although other EV types can also cause severe complications. Potential virulence or phenotype-associated sites were identified in the genomes of CVA6, CVA16, and CVA10.

Tracking echovirus eleven outbreaks in Guangdong, China: a metatranscriptomic, phylogenetic, and epidemiological study.

In April 2019, a suspect cluster of enterovirus cases was reported in a neonatology department in Guangdong, China, resulting in five deaths. We aimed to investigate the pathogen profiles in fatal cases, the circulation and transmission pattern of the viruses by combining metatranscriptomic, phylogenetic, and epidemiological analyses. Metatranscriptomic sequencing was used to characterize the enteroviruses. Clinical and environmental surveillance in the local population was performed to understand the prevalence and genetic diversity of the viruses in the local population. The possible source(s), evolution, transmission, and recombination of the viruses were investigated by incorporating genomes from the current outbreak, from local retrospective surveillance, and from public databases. Metatranscriptomic analysis identified Echovirus 11 (E11) in three fatal cases. Seroprevalence of neutralization antibody to E11 was 35 to 44 per cent in 3-15 age groups of general population, and the viruses were associated with various clinical symptoms. From the viral phylogeny, nosocomial transmissions were identified and all E11 2019 outbreak strains were closely related with E11 strains circulating in local population 2017-19. Frequent recombination occurred among the 2019 Guangdong E11 outbreak strains and various genotypes in enterovirus B species. This study provides an example of combining advanced genetic technology and epidemiological surveillance in pathogen diagnosis, source(s), and transmission tracing during an infectious disease outbreak. The result highlights the hidden E11 circulation and the risk of viral transmission and infection in the young age population in China. Frequent recombination between Guangdong-like strains and other enterovirus genotypes also implies the prevalence of these emerging E11 strains.

Emerging recombination of the C2 sub-genotype of HFMD-associated CV-A4 is persistently and extensively circulating in China.

Sporadic outbreaks caused by coxsackievirus A4 (CV-A4) have been reported worldwide. To further elucidate the detailed genetic characteristics and evolutionary recombination events of CV-A4, virus samples from nationwide hand, foot and mouth disease (HFMD) surveillance, encompassing 27 out of the 31 provinces in China, were investigated. Comprehensive and systematic phylogenetic analyses were performed by using 29 complete genomes, 142 complete CV-A4 VP1 sequences. Four genotypes (A, B, C and D) and five sub-genotypes (C1-C5) were re-identified based on the complete VP1 sequences. C2 is the predominant sub-genotype of CV-A4 associated with HFMD and has evolved into 3 clusters. Cluster 1 is a major cluster that has been persistently and extensively circulating in China since 2006 and has been associated with all severe cases. All the sequences showed high homology with the CV-A4 prototype in the P1 region, while higher identities with CV-A5, CV-14 and CV-16 in the P2 and P3 regions. Recombination analysis revealed that C2 had two specific genetic recombination patterns with other EV-A prototypes in the 5'-UTR and 3D region compared with C5. These recombination patterns might be associated with the increased transmissibility of C2 viruses, which were obtained due to their persistent and extensive circulation in populations.

The emerging sub-genotype C2 of CoxsackievirusA10 Associated with Hand, Foot and Mouth Disease extensively circulating in mainland of China.

Coxsackievirus A10 (CV-A10) associated with Hand, foot, and mouth disease (HFMD) cases emerged increasingly in recent years. In this study, the samples from nation-wide HFMD surveillance, including 27 out of 31 provinces in China were investigated, and the continuous and extensive virological surveillance, covered 13 years, were conducted to provide a comprehensive molecular characterization analysis of CV-A10. 855 CV-A10 viruses (33 severe cases included), were isolated from HFMD children patients during 2009 to 2016 in China. 164 representative sequences from these viruses, together with 117 CV-A10 sequences downloaded from GenBank based on entire VP1 were recruited in this study. Two new genotypes (F and G) and two sub-genotypes (C1 and C2) were identified. Among 264 Chinese sequences, 9 of them were genotype B, 8 of them were C1, and the other (247) were C2, the predominant sub-genotype in China since 2012. Chinese C2 viruses showed obvious temporal characteristics and can be divided into 3 clusters (cluster 1~3). Cluster 3 viruses was circulating extensively during 2014 and 2016 with more severe cases. It is very necessary and important to continuously conduct the extensive virological surveillance for CV-A10, and further evolutionary studies will provide more evidence on its evolution and virulence.

Association of GII.P16-GII.2 Recombinant Norovirus Strain with Increased Norovirus Outbreaks, Guangdong, China, 2016.

An unusual prevalence of recombinant GII.2 noroviruses (GII.P16-GII.2) in Guangdong, China, at the end of 2016 caused a sharp increase in outbreaks of acute gastroenteritis. This event was another non-GII.4 epidemic that emerged after the GII.17 viruses in 2014 and 2015 and warrants global surveillance.

Isolation and Characterization of a Highly Mutated Chinese Isolate of Enterovirus B84 from a Patient with Acute Flaccid Paralysis.

Enterovirus B84 (EV-B84) is a newly identified serotype within the species Enterovirus B (EV-B). To date, only ten nucleotide sequences of EV-B84 are published and only one full-length genome sequence (the prototype strain) is available in the GenBank database. Here, a highly mutated EV-B84 (strain AFP452/GD/CHN/2004) was recovered from a patient with acute flaccid paralysis in the Guangdong province of China in 2004 making this the first report of EV-B84 in China. Sequence comparison and phylogenetic dendrogram analysis revealed high variation from the global EV-B84 strains (African and Indian strains) and frequent intertypic recombination in the non-structural protein region, suggesting high genetic diversity in EV-B84. The Chinese EV-B84 strain, apparently evolving independently of the other ten strains, strongly suggests that the EV-B84 strain has been circulating for many years. However, the extremely low isolation rate suggests that it is not a prevalent EV serotype in China or worldwide. This study provides valuable information about the molecular epidemiology of EV-B84 in China, and will be helpful in future studies to understand the association of EV-B84 with neurological disorders; it also helps expand the number of whole virus genome sequences of EV-B84 in the GenBank database.

Evaluating the transmission routes of hand, foot, and mouth disease in Guangdong, China.

Although it is an enteroviral infectious disease, recent studies suggest that respiratory transmission might play a role in the transmission of hand, foot, and mouth disease (HFMD). We evaluated the transmission modes (respiratory and fecal-oral transmission) of HFMD among children using a case-control study in Guangdong, China. Our analyses suggested that fecal-oral transmission might be the principal transmission mode of HFMD among children in the study area, and handwashing habits of the children and their parents should be emphasized to control this infection.

[The Complete Sequence Analysis of 18 Strains of Coxsackievirus A6 in Guangdong Province of China].

In this study, we examined the complete genome of coxsackievirus A6 (CVA6) from hand, foot, and mouth disease in Guangdong Province from 2013,and explored the genetic similarities and differences in epidemic and non-epidemic stains of CVA6.Eighteen strains of CVA6 were included in complete genome sequencing, and the sequences were subject to phylogenetic analysis,sequence alignment analysis and genetic recombination analysis using the software DNASTAR6.0,MEGA5.2and SimPlot3.5.1.The results showed that the complete genome of 18 Guangdong CVA6strains ranged from 7390bp to 7392bp.No insertions or deletions were detected in the coding region. There were several insertions and deletions in 5'UTR and 3'UTR.Phylogenetic analysis indicated that the nucleotide and amino acid sequence identity between the 18 complete genomes were 90.5%-99.6% and 97.5%-99.9%,respectively.The strains isolated in2013 could be further divided into two clusters, III and IV, while the strains isolated in 2011 were only present in the IV cluster. Genetic recombination analysis revealed that the Guangdong representative strain of CVA6,GD870/2013,had gene recombination in the P2 and P3regions,while the GD839/2013 strain did not show obvious genetic recombination. Genome-wide analysis of CVA6 revealed that there are two possible transmitted chains, III and IV, in epidemic strains from Guangdong Province in 2013.The transmitted chain Ⅲ originated from the strain with genetic recombination in the P2 and P3regions,whichwas completely different from the chain IV. Transmission of chain IV of CVA6 was only observed in the nonepidemic 2011 strain.