Phylogeography of Coxsackievirus A16 Reveals Global Transmission Pathways and Recent Emergence and Spread of a Recombinant Genogroup

Synergetic association between coxsackievirus A16 genotype evolution and recombinant form shifts

Coxsackievirus A16 (CVA16) is a major pathogen that causes hand, foot, and mouth disease (HFMD). The recombination form (RF) shifts and global transmission dynamics of CVA16 remain unknown. In this retrospective study, global sequences of CVA16 were retrieved from the GenBank database and analyzed using comprehensive phylogenetic inference, RF surveys, and population structure. A total of 1,663 sequences were collected, forming a 442-sequences dataset for VP1 coding region analysis and a 345-sequences dataset for RF identification. Based on the VP1 coding region used for serotyping, three genotypes (A, B, and D), two subgenotypes of genotype B (B1 and B2), and three clusters of subgenotype B1 (B1a, B1b, and B1c) were identified. Cluster B1b has dominated the global epidemics, B2 disappeared in 2000, and D is an emerging genotype dating back to August 2002. Globally, four oscillation phases of CVA16 evolution, with a peak in 2013, and three migration pathways were identified. Europe, China, and Japan have served as the seeds for the global transmission of CVA16. Based on the 3D coding region of the RFs, five clusters of RFs (RF-A to -E) were identified. The shift in RFs from RF-B and RF-C to RF-D was accompanied by a change in genotype from B2 to B1a and B1c and then to B1b. In conclusion, the evolution and population dynamics of CVA16, especially the coevolution of 3D and VP1 genes, revealed that genotype evolution and RF replacement were synergistic rather than stochastic.

Molecular evolutionary dynamics of enterovirus A71, coxsackievirus A16 and coxsackievirus A6 causing hand, foot and mouth disease in Thailand, 2000-2022

Hand, foot and mouth disease (HFMD) is a public health threat worldwide, particularly in the Asia-Pacific region. Enterovirus A71 (EV-A71), coxsackievirus A16 (CVA16), and CVA6 are the major pathogens causing HFMD outbreaks in several countries, including Thailand. We retrieved 385 VP1 nucleotide sequences, comprising 228 EV-A71, 33 CVA16, and 124 CVA6, deposited in the databases between 2000 and 2022 for molecular evolutionary characterization using Bayesian phylogeny. All EV-A71 identified belonged to genotype B, subgenotypes B4, and B5, and to genotype C, subgenotypes C1, C2, C4a, C4b, and C5. The analyzes demonstrated these viruses' co-circulation and subgenotypic changes throughout the past two decades. The CVA16 was grouped in genotype B1, predominantly subgenotype B1a, and the CVA6 was grouped in subgenotype D3, clades 1-4. The tMRCA of EV-A71 genotypes B and C, CVA16 B1, and CVA6 D3 dated 1993.79, 1982.62, 1995.86, and 2007.31, respectively, suggesting that the viruses were likely introduced and cryptically circulated in Thailand before the HFMD cases were recognized. We demonstrated these viruses' fluctuation and cyclical pattern throughout the two decades of observation. This study provided insight into evolutionary dynamics concerning molecular epidemiology and supported the selection of current genotype-matched vaccines, vaccine development, and implementation.

Phylogenetics and phylogeographic characteristics of coxsackievirus A16 in hand foot and mouth disease and herpangina cases collected in Beijing, China from 2019 to 2021

Coxsackievirus A16 (CV-A16) is a significant pathogen responsible for causing hand foot and mouth disease (HFMD) and herpangina (HA). This study aimed to investigate the recent evolution and spread of CV-A16 by monitoring HFMD and HA cases in 29 hospitals across 16 districts in Beijing from 2019 to 2021. The first five cases of HFMD and the first five cases of HA each month in each hospital were included in the study. Real-time reverse transcription polymerase chain reaction was used to identify CV-A16, CV-A6, and EV-A71. From each district, two to four CV-A16 positive samples with a relatively long sampling time interval every month were selected for sequencing. A total of 3344 HFMD cases and 2704 HA cases were enrolled in this study, with 76.0% (2541/3344) of HFMD and 45.4% (1227/2704) of HA cases confirmed to be infected by enterovirus. Among the EV-positive samples, CV-A16 virus was detected in 33.61% (854/2541) of HFMD cases and 13.4% (165/1227) of HA cases, with the predominant cluster being B1a. Both B1a and B1b had a co-circulation of local and imported strains, with different origin time (1993 vs. 1995), different global distribution (14 countries vs. 10 countries), and different transmission centers but mainly distributed in the southern and eastern regions of Beijing. Strengthening surveillance of HFMD in southern and eastern regions will improve the prevention and control efficiency of enterovirus infections.

Current status of hand-foot-and-mouth disease

Hand-foot-and-mouth disease (HFMD) is a viral illness commonly seen in young children under 5 years of age, characterized by typical manifestations such as oral herpes and rashes on the hands and feet. These symptoms typically resolve spontaneously within a few days without complications. Over the past two decades, our understanding of HFMD has greatly improved and it has received significant attention. A variety of research studies, including epidemiological, animal, and in vitro studies, suggest that the disease may be associated with potentially fatal neurological complications. These findings reveal clinical, epidemiological, pathological, and etiological characteristics that are quite different from initial understandings of the illness. It is important to note that HFMD has been linked to severe cardiopulmonary complications, as well as severe neurological sequelae that can be observed during follow-up. At present, there is no specific pharmaceutical intervention for HFMD. An inactivated Enterovirus A71 (EV-A71) vaccine that has been approved by the China Food and Drug Administration (CFDA) has been shown to provide a high level of protection against EV-A71-related HFMD. However, the simultaneous circulation of multiple pathogens and the evolution of the molecular epidemiology of infectious agents make interventions based solely on a single agent comparatively inadequate. Enteroviruses are highly contagious and have a predilection for the nervous system, particularly in child populations, which contributes to the ongoing outbreak. Given the substantial impact of HFMD around the world, this Review synthesizes the current knowledge of the virology, epidemiology, pathogenesis, therapy, sequelae, and vaccine development of HFMD to improve clinical practices and public health efforts.

Epidemiology of hand, foot, and mouth disease and the genetic characteristics of Coxsackievirus A16 in Taiyuan, Shanxi, China from 2010 to 2021

Hand, foot, and mouth disease (HFMD) is a common childhood infectious disease caused by human enteroviruses (EV). This study aimed to describe the epidemiological features of HFMD and the genetic characteristics of Coxsackievirus A16 (CVA16) in Taiyuan, Shanxi, China, from 2010 to 2021. Descriptive epidemiological methods were used to analyze the time and population distribution of HFMD and the genetic characteristics of CVA16. Except being affected by the COVID-19 epidemic in 2020, HFMD epidemics were sporadic from January to March each year, and began to increase in April, with a major epidemic peak from May to August, which declined in September, followed by a secondary peak from October to December. The prevalence of EV infection was the highest in children aged one to five years (84.42%), whereas its incidence was very low in children under one year of age (5.48%). Enterovirus nucleic acid was detected by real-time reverse transcription polymerase chain reaction in 6641 clinical specimens collected from patients with HFMD from 2010 to 2021, and 4236 EV-positive specimens were detected, including 988 enterovirus A71 (EV-A71), 1488 CVA16, and 1760 other enteroviruses. CVA16 remains prevalent and has co-circulated with other EVs in Taiyuan from 2010 to 2021. A phylogenetic tree constructed based on the VP1 region showed that all CVA16 strains belonged to two different clades of the B1 genotype, B1a and B1b. They showed a nucleotide similarity of 86.5-100%, and an amino acid similarity of 96.9-100%. Overall, these findings add to the global genetic resources of CVA16, demonstrate the epidemiological characteristics of HFMD as well as the genetic features of CVA16 in Taiyuan City during 2010-2021, and provide supporting evidence for the prevention and control of HFMD.

The Intrapopulation Genetic Diversity of RNA Virus May Influence the Sensitivity of Chlorine Disinfection

RNA virus populations are not clonal; rather, they comprise a mutant swarm in which sequences are slightly different from the master sequence. Genetic diversity within a population (intrapopulation genetic diversity) is critical for RNA viruses to survive under environmental stresses. Disinfection has become an important practice in the control of pathogenic viruses; however, the impact of intrapopulation genetic diversity on the sensitivity of disinfection, defined as -log10 (postdisinfected infectious titer/predisinfected titer), has not been elucidated. In this study, we serially passaged populations of rhesus rotavirus. We demonstrated that populations with reduced chlorine sensitivity emerged at random and independently of chlorine exposure. Sequencing analysis revealed that compared with sensitive populations, less-sensitive ones had higher non-synonymous genetic diversity of the outer capsid protein gene, suggesting that changes in the amino acid sequences of the outer capsid protein were the main factors influencing chlorine sensitivity. No common mutations were found among less-sensitive populations, indicating that rather than specific mutations, the diversity of the outer capsid protein itself was associated with the disinfection sensitivity and that the disinfection sensitivity changed stochastically. Simulation results suggest that the disinfection sensitivity of a genetically diverse population is destabilized if cooperative viral clusters including multiple sequences are formed. These results advocate that any prevention measures leading to low intrapopulation genetic diversity are important to prevent the spread and evolution of pathogenic RNA viruses in society.

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.

•

Systematic profiles of genetic features of CVA10 near complete genome.

•

First report of the appearance of CVA10 genogroup D in Chinese mainland.

•

Genomic comparisons indicate the potential recombinant origin of CVA10 genogroup D.

Enterovirus Surveillance (EVSurv) in Germany

The major aim of the enterovirus surveillance (EVSurv) in Germany is to prove the absence of poliovirus circulation in the framework of the Global Polio Eradication Program (GPEI). Therefore, a free-of-charge enterovirus diagnostic is offered to all hospitals for patients with symptoms compatible with a polio infection. Within the quality proven laboratory network for enterovirus diagnostic (LaNED), stool and cerebrospinal fluid (CSF) samples from patients with suspected aseptic meningitis/encephalitis or acute flaccid paralysis (AFP) are screened for enterovirus (EV), typing is performed in all EV positive sample to exclude poliovirus infections. Since 2006, ≈200 hospitals from all 16 German federal states have participated annually. On average, 2500 samples (70% stool, 28% CSF) were tested every year. Overall, the majority of the patients studied are children <15 years. During the 15-year period, 53 different EV serotypes were detected. While EV-A71 was most frequently detected in infants, E30 dominated in older children and adults. Polioviruses were not detected. The German enterovirus surveillance allows monitoring of the circulation of clinically relevant serotypes resulting in continuous data about non-polio enterovirus epidemiology.

Molecular epidemiology of coxsackievirus A16 circulating in children in Beijing, China from 2010 to 2019

BACKGROUND

Coxsackievirus A16 (CVA16) is one of the major etiological agents of hand, foot and mouth disease (HFMD). This study aimed to investigate the molecular epidemiology and evolutionary characteristics of CVA16.

METHODS

Throat swabs were collected from children with HFMD and suspected HFMD during 2010-2019. Enteroviruses (EVs) were detected and typed by real-time reverse transcription-polymerase chain reaction (RT-PCR) and RT-PCR. The genotype, evolutionary rate, the most recent common ancestor, population dynamics and selection pressure of CVA16 were analyzed based on viral protein gene (VP1) by bioinformatics software.

RESULTS

A total of 4709 throat swabs were screened. EVs were detected in 3180 samples and 814 were CVA16 positive. More than 81% of CVA16-positive children were under 5 years old. The prevalence of CVA16 showed obvious periodic fluctuations with a high level during 2010-2012 followed by an apparent decline during 2013-2017. However, the activities of CVA16 increased gradually during 2018-2019. All the Beijing CVA16 strains belonged to sub-genotype B1, and B1b was the dominant strain. One B1c strain was detected in Beijing for the first time in 2016. The estimated mean evolutionary rate of VP1 gene was 4.49 × 10^-3 substitution/site/year. Methionine gradually fixed at site-23 of VP1 since 2012. Two sites were detected under episodic positive selection, one of which (site-223) located in neutralizing linear epitope PEP71.

CONCLUSIONS

The dominant strains of CVA16 belonged to clade B1b and evolved in a fast evolutionary rate during 2010-2019 in Beijing. To provide more favorable data for HFMD prevention and control, it is necessary to keep attention on molecular epidemiological and evolutionary characteristics of CVA16.

Coxsackieviruses A6 and A16 associated with hand, foot, and mouth disease in Vietnam, 2008-2017: Essential information for rational vaccine design

Development of multivalent hand, foot, and mouth disease (HFMD) vaccines against enterovirus A71 (EV-A71) and several non-EV-A71 enteroviruses is needed for this life-threatening disease with a huge economic burden in Asia-Pacific countries. Comprehensive studies on the molecular epidemiology and genetic and antigenic characterization of major causative enteroviruses will provide information for rational vaccine design. Compared with molecular studies on EV-A71, that for non-EV-A71 enteroviruses remain few and limited in Vietnam. Therefore, we conducted a 10-year study on the circulation and genetic characterization of coxsackievirus A16 (CV-A16) and CV-A6 isolated from patients with HFMD in Northern Vietnam between 2008 and 2017. Enteroviruses were detected in 2228 of 3212 enrolled patients. Of the 42 serotypes assigned, 28.4% and 22.4% accounted for CV-A6 and CV-A16, being the second and the third dominant serotypes after EV-A71 (31.7%), respectively. The circulation of CV-A16 and CV-A6 showed a wide geographic distribution and distinct periodicity. Phylogenetic analyses revealed that the majority of Vietnamese CV-A6 and CV-A16 strains were located within the largest sub-genotypes or sub-genogroups. These comprised strains isolated from patients with HFMD worldwide during the past decade and the Vietnamese strains have been evolving in a manner similar to the strains circulating worldwide. Amino acid sequences of the putative functional loops on VP1 and other VPs among Vietnamese CV-A6 and CV-A16 isolates were highly conserved. Moreover, the functional loop patterns of VP1 were similar to the dominant patterns found worldwide, except for the T164K substitution on the EF loop in Vietnamese CV-A16. The findings suggest that the development of a universal HFMD vaccine, at least in Vietnam, must target CV-A6 and CV-A16 as two of the three major HFMD-causing serotypes. Vietnamese isolates or their genome sequences can be considered for rational vaccine design.

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.

Global Spread of the B5 Subgenotype EV-A71 and the Phylogeographical Analysis of Chinese Migration Events

The subgenotype B5 of EV-A71 is a widely circulating subgenotype that frequently spreads across the globe. Several outbreaks have occurred in nations, such as Malaysia, Thailand, Vietnam, and Japan. Appearing first in Taiwan, China, the subgenotype has been frequently reported in mainland of China even though no outbreaks have been reported so far. The current study reconstructed the migration of the B5 subgenotype of EV-A71 in China via phylogeographical analysis. Furthermore, we investigated its population dynamics in order to draw more credible inferences. Following a dataset cleanup of B5 subgenotype of EV-A71, we detected earlier B5 subgenotypes of EV-A71 sequences that had been circulating in Malaysia and Singapore since the year 2000, which was before the 2003 outbreak that occurred in Sarawak. The Bayesian inference indicated that the most recent common ancestor of B5 subgenotype EV-A71 appeared in September, 1994 (1994.75). With respect to the overall prevalence, geographical reconstruction revealed that the B5 subgenotype EV-A71 originated singly from single-source cluster and subsequently developed several active lineages. Based on a large amount of data that was accumulated, we conclude that the appearance of the B5 subgenotype of EV-A71 in mainland of China was mainly due to multiple migrations from different origins.

Genetic characterization of VP1 of coxsackieviruses A2, A4, and A10 associated with hand, foot, and mouth disease in Vietnam in 2012-2017: endemic circulation and emergence of new HFMD-causing lineages

While conducting sentinel surveillance of hand, foot, and mouth disease (HFMD) in Vietnam, we found a sudden increase in the prevalence of coxsackievirus A10 (CV-A10) in 2016 and CV-A2 and CV-A4 in 2017, the emergence of which has been reported recently to be associated with various clinical manifestations in other countries. However, there have been only a limited number of molecular studies on those serotypes, with none being conducted in Vietnam. Therefore, we sequenced the entire VP1 genes of CV-A10, CV-A4, and CV-A2 strains associated with HFMD in Vietnam between 2012 and 2017. Phylogenetic analysis revealed a trend of endemic circulation of Vietnamese CV-A10, CV-A4, and CV-A2 strains and the emergence of thus-far undescribed HFMD-causing lineages of CV-A4 and CV-A2. The Vietnamese CV-A10 strains belonged to a genotype comprising isolates from patients with HFMD from several other countries; however, most of the Vietnamese strains were grouped into a local lineage. Recently, emerging CV-A4 strains in Vietnam were grouped into a unique lineage within a genotype comprising strains isolated from patients with acute flaccid paralysis from various countries. New substitutions were detected in the putative BC and HI loops in the Vietnamese CV-A4 strains. Except for one strain, Vietnamese CV-A2 isolates were grouped into a unique lineage of a genotype that includes strains from various countries that are associated with other clinical manifestations. Enhanced surveillance is required to monitor their spread and to specify their roles as etiological agents of HFMD or "HFMD-like" diseases, especially for CV-A4 and CV-A2. Further studies including whole-genome sequencing should be conducted to fully understand the evolutionary changes occurring in these newly emerging strains.

Monitoring of enterovirus diversity in wastewater by ultra-deep sequencing: An effective complementary tool for clinical enterovirus surveillance

In a one-year (October 2014-October 2015) pilot study, we assessed wastewater monitoring with sustained sampling for analysis of global enterovirus (EV) infections in an urban community. Wastewater was analysed by ultra-deep sequencing (UDS) after PCR amplification of the partial VP1 capsid protein gene. The nucleotide sequence analysis showed an unprecedented diversity of 48 EV types within the community, which were assigned to the taxonomic species A (n = 13), B (n = 23), and C (n = 12). During the same period, 26 EV types, of which 22 were detected in wastewater, were identified in patients referred to the teaching hospital serving the same urban population. Wastewater surveillance detected a silent circulation of 26 EV types including viruses reported in clinically rare respiratory diseases. Wastewater monitoring as a supplementary procedure can complement clinical surveillance of severe diseases related to non-polio EVs and contribute to the final stages of poliomyelitis eradication.

Multirecombinant Enterovirus A71 Subgenogroup C1 Isolates Associated with Neurologic Disease, France, 2016-2017

In 2016, an upsurge of neurologic disease associated with infection with multirecombinant enterovirus A71 subgenogroup C1 lineage viruses was reported in France. These viruses emerged in the 2000s; 1 recombinant is widespread. This virus lineage has the potential to be associated with a long-term risk for severe disease among children.

Recombination in Enteroviruses, a Multi-Step Modular Evolutionary Process

RNA recombination is a major driving force in the evolution and genetic architecture shaping of enteroviruses. In particular, intertypic recombination is implicated in the emergence of most pathogenic circulating vaccine-derived polioviruses, which have caused numerous outbreaks of paralytic poliomyelitis worldwide. Recent experimental studies that relied on recombination cellular systems mimicking natural genetic exchanges between enteroviruses provided new insights into the molecular mechanisms of enterovirus recombination and enabled to define a new model of genetic plasticity for enteroviruses. Homologous intertypic recombinant enteroviruses that were observed in nature would be the final products of a multi-step process, during which precursor nonhomologous recombinant genomes are generated through an initial inter-genomic RNA recombination event and can then evolve into a diversity of fitter homologous recombinant genomes over subsequent intra-genomic rearrangements. Moreover, these experimental studies demonstrated that the enterovirus genome could be defined as a combination of genomic modules that can be preferentially exchanged through recombination, and enabled defining the boundaries of these recombination modules. These results provided the first experimental evidence supporting the theoretical model of enterovirus modular evolution previously elaborated from phylogenetic studies of circulating enterovirus strains. This review summarizes our current knowledge regarding the mechanisms of recombination in enteroviruses and presents a new evolutionary process that may apply to other RNA viruses.

Molecular surveillance of coxsackievirus A16 reveals the emergence of a new clade in mainland China

Coxsackievirus A16 (CV-A16) of the genotypes B1a and B1b have co-circulated in mainland China in the past decades. From 2013 to 2017, a total of 3,008 specimens from 3,008 patients with mild hand, foot, and mouth disease were collected in the present study. Viral RNA was tested for CV-A16 by a real-time RT-PCR method, and complete VP1 sequences and full-length genome sequences of CV-A16 strains from this study were determined by RT-PCR and sequencing. Sequences were analyzed using a series of bioinformatics programs. The detection rate for CV-A16 was 4.1%, 25.9%, 10.6%, 28.1% and 12.9% in 2013, 2014, 2015, 2016 and 2017, respectively. Overall, the detection rate for CV-A16 was 16.5% (497/3008) in this 5-year period in Shenzhen, China. One hundred forty-two (142/155, 91.6%) of the 155 genotype B1 strains in the study belonged to subgenotype B1b, and 13 (13/155, 8.4%) strains belonged to subgenotype B1a. Two strains (CVA16/Shenzhen174/CHN/2017 and CVA16/Shenzhen189/CHN/2017) could not be assigned to a known genotype. Phylogenetic analysis of these two strains and other Chinese CV-A16 strains indicated that these two CV-A16 strains clustered independently in a novel clade whose members differed by 8.4%-11.8%, 8.4%-12.1%, and 14.6%-14.8% in their nucleotide sequences from those of Chinese B1a, B1b, and genotype D strains, respectively. Phylogenetic analysis of global CV-A16 strains further indicated that the two novel CV-A16 strains from this study grouped in a previously uncharacterized clade, which was designated as the subgenogroup B3 in present study. Meanwhile, phylogenetic reconstruction revealed two other new genotypes, B1d and B4, which included a Malaysian strain and two American strains, respectively. The complete genome sequences of the two novel CV-A16 strains showed the highest nucleotide sequence identity of 92.3% to the Malaysian strain PM-15765-00 from 2000. Comparative analysis of amino acid sequences of the two novel CV-A16 strains and their relatives suggested that variations in the nonstructural proteins may play an important role in the evolution of modern CV-A16.

Molecular characterization of Coxsackievirus A16 strains isolated from children with severe hand, foot, and mouth disease in Yunnan, Southwest China, during 2009-2015

Coxsackievirus A16 (CV-A16) commonly causes mild symptoms, but severe diseases, such as aseptic meningitis, encephalitis, and even fatal cases, have been reported. Thirteen CV-A16 strains were isolated from patients with severe hand, foot, and mouth disease in Yunnan, Southwest China, from 2009 to 2015. Subgenotype B1a and B1b of CV-A16 were predominantly circulating the region with B1b the predominant strain in recent years. The mean rate of nucleotide substitution based on the VP1 gene sequence was 4.545 × 10 ^-3  substitution per site per year from 2009 to 2015. These results may help in understanding the genetic diversity of CV-A16 and develop a CV-A16 vaccine.

Environmental Surveillance Reveals Complex Enterovirus Circulation Patterns in Human Populations

Background

Enteroviruses are common human pathogens occasionally associated with severe disease, notoriously paralytic poliomyelitis caused by poliovirus. Other enterovirus serotypes such as enterovirus A71 and D68 have been linked to severe neurological syndromes. New enterovirus serotypes continue to emerge, some believed to be derived from nonhuman primates. However, little is known about the circulation patterns of many enterovirus serotypes and, in particular, the detailed enterovirus composition of sewage samples.

Methods

We used a next-generation sequencing approach analyzing reverse transcriptase polymerase chain reaction products synthesized directly from sewage concentrates.

Results

We determined whole-capsid genome sequences of multiple enterovirus strains from all 4 A to D species present in environmental samples from the United Kingdom, Senegal, and Pakistan.

Conclusions

Our results indicate complex enterovirus circulation patterns in human populations with differences in serotype composition between samples and evidence of sustained and widespread circulation of many enterovirus serotypes. Our analyses revealed known and divergent enterovirus strains, some of public health relevance and genetically linked to clinical isolates. Enteroviruses identified in sewage included vaccine-derived poliovirus and enterovirus D-68 stains, new enterovirus A71 and coxsackievirus A16 genogroups indigenous to Pakistan, and many strains from rarely reported serotypes. We show how this approach can be used for the early detection of emerging pathogens and to improve our understanding of enterovirus circulation in humans.

Antigenic characteristics and genomic analysis of novel EV-A90 enteroviruses isolated in Xinjiang, China

Enterovirus A90 (EV-A90) is a novel serotype of enterovirus A species that is rarely reported. Here, we isolated five enteroviruses from patients with acute flaccid paralysis in Hotan and Kashgar cities in Xinjiang, China that were identified as EV-A90 by molecular typing. The VP1 sequences of these Xinjiang EV-A90 strains showed 88.4–89% nucleotide sequence identity to the prototype EV-A90 strain; however, genome analysis indicated complex recombination events in P2 and P3 regions. Next, the seroprevalence of EV-A90 was examined in 49 serum specimens collected in Hotan and Kashgar, and 37.5% were EV-A90 antibody positive (>1:8), with a geometric mean titre (GMT) of 1:10.47. The low positive rate and GMT suggest a low-level EV-A90 epidemic in Xinjiang. Two of the five Xinjiang EV-A90 strains were temperature sensitive, and three were temperature resistant, and a comparative genomics analysis suggested that an amino acid substitution (H1799Y) in the 3D^pol region was related to temperature sensitivity. Although the epidemic strength is low, some EV-A90 strains were temperature resistant, which is suggestive of strong virulence and transmission capacity. This study expanded the number of EV-A90 in GenBank and provided basic data that may be useful for studying the molecular epidemiology of EV-A90.

Large-scale HFMD epidemics caused by Coxsackievirus A16 in Bangalore, India during 2013 and 2015

Hand foot and mouth disease (HFMD) is a relatively unreported disease in India. This study was undertaken to characterize the enterovirus type/s associated with two unexpectedly-massive epidemics that occurred in Bangalore, India in 2013 and 2015. Stool samples of 229 children with HFMD living in Northern and Southern areas of Bangalore were tested by RT-PCR; 189 (82.5%) were enterovirus positive. The Indian CV-A16 strains exhibited 98-99% sequence identity with those reported in France and China in the 5' untranslated region. BLAST and phylogenetic analyses of complete genomes of representative Indian isolates revealed that the 2015 epidemic was predominated by an inter-species recombinant between CV-A16 and coxsackievirus B5. The 2013 epidemic was primarily caused by nonrecombinant strains. The CV-A16 strains circulated in India since 2007 and phylogeographic analyses indicated imported cases in France and China. In conclusion, CV-A16-associated HFMD epidemics should be recognized as an emerging public health problem in India.