Species-specific RT-PCR amplification of human enteroviruses: a tool for rapid species identification of uncharacterized enteroviruses

Evaluation of VP4-VP2 sequencing for molecular typing of human enteroviruses

Enteroviruses and rhinoviruses are highly diverse, with over 300 identified types. Reverse transcription-polymerase chain reaction (RT-PCR) assays targeting their VP1, VP4, and partial VP2 (VP4-pVP2) genomic regions are used for detection and identification. The VP4-pVP2 region is particularly sensitive to RT-PCR detection, making it efficient for clinical specimen analysis. However, a standard type identification method using this region is lacking. This study aimed to establish such a method by examining the divergence of VP4-pVP2 amino acid sequences between enterovirus and rhinovirus prototypes. Pairwise analysis of 249 types indicated a 95% threshold for enterovirus intra-species identification but not for rhinovirus prototypes. Protein BLAST search analyses of representative enterovirus prototypes, including EV-A71, EV-D68, CVA6, CVA10, CVA16, and polioviruses (PVs), validated the 95% threshold for typing, with a few exceptions such as PV1-PV2 and CVA6-CVA10, as well as some EV-C types. This study proposes a criterion for typing based on VP4-pVP2 amino acids, which can aid in rapid enterovirus diagnosis during routine clinical or environmental surveillance and emergency outbreaks. Our research confirms the reliability of the suggested VP4-pVP2-based threshold for typing and its potential application in laboratory settings.

Genetic variation and evolutionary characteristics of Echovirus 11: new variant within genotype D5 associated with neonatal death found in China

Echovirus 11 (E11) has gained attention owing to its association with severe neonatal infections. From 2018 to 2023, a surge in severe neonatal cases and fatalities linked to a novel variant of genotype D5 was documented in China, France, and Italy. However, the prevention and control of E11 variants have been hampered by limited background data on the virus circulation and genetic variance. Therefore, the present study investigated the circulating dynamics of E11 and the genetic variation and molecular evolution of genotype D5 through the collection of strains from the national acute flaccid paralysis (AFP) and hand, foot, and mouth disease (HFMD) surveillance system in China during 2000-2022 and genetic sequences published in the GenBank database. The results of this study revealed a prevalent dynamic of E11 circulation, with D5 being the predominant genotype worldwide. Further phylogenetic analysis of genotype D5 indicated that it could be subdivided into three important geographic clusters (D5-CHN1: 2014-2019, D5-CHN2: 2016-2022, and D5-EUR: 2022-2023). Additionally, variant-specific (144) amino acid mutation sites and positive-selection pressure sites (132, 262) were identified in the VP1 region. Cluster-specific recombination patterns were also identified, with CVB5, E6, and CVB4 as the major recombinant viruses. These findings provide a preliminary landscape of E11 circulation worldwide and basic scientific data for further study of the pathogenicity of E11 variants.

Epidemiology of Hand, Foot, and Mouth Disease and Genetic Characterization of Coxsackievirus A16 in Shenyang, Liaoning Province, China, 2013-2023

Hand, foot, and mouth disease (HFMD), a common childhood infection caused by enterovirus, poses a serious public health concern in China. We collected and analyzed epidemiological data on 62,133 HFMD cases in Shenyang City, Liaoning Province, from 2013 to 2023. The average annual incidence was 76.12 per 100,000 person-years; 99.45% of cases were mild, while 0.55% were severe. Only one patient died. HFMD infections peaked annually in July. Children in kindergartens and scattered children accounted for 44.6% and 42.2% of cases, respectively. Real-time RT-PCR detection of enteroviruses in 5534 patient samples revealed 3780 positives, of which 25.1% were CVA16-positive. Positives were randomly sampled, yielding 240 VP1 sequences of CVA16. Phylogenetic tree results showed that all VP1 sequences belonged to the B1 sub-genogroup. However, the sub-genogroup prevalence varied over time: from 2013 to 2014 and 2019 to 2021, the predominant sub-genogroup was B1a, while it was B1b from 2015 to 2018. Further phylogenetic analyses showed substantial divergence between B1a branches in CVA16, suggesting possible turnover of the B1a sub-genogroup in CVA16 due to evolution. This study provides epidemiological data on HFMD in Shenyang, and provides a phylogenetic analysis of CVA16, offering a theoretical basis for preventing and controlling HFMD in Shenyang City.

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.

Identification and genetic characterization of a recently identified enterovirus C116 in China

Enterovirus C116 (EV-C116) is a new member of the enterovirus C group which is closely associated with several infectious diseases. Although sporadic studies have detected EV-C116 in clinical samples worldwide, there is currently limited information available. In this study, two EV-C-positive fecal specimens were detected in apparently healthy children, which harbored low abundance, through meta-transcriptome sequencing. Based on the prototypes of several EV-Cs, two lineages were observed. Lineage 1 included many types that could not cause EV-like cytopathic effect in cell culture. Three genogroups of EV-C116 were divided in the maximum likelihood tree, and the two strains in this study (XZ2 and XZ113) formed two different lineages, suggesting that EV-C116 still diffuses worldwide. Obvious inter-type recombination events were observed in the XZ2 strain, with CVA22 identified as a minor donor. However, another strain (XZ113) underwent different recombination situations, highlighting the importance of recombination in the formation of EV-Cs biodiversity. The EV-C116 strains could propagate in rhabdomyosarcoma cell cultures at low titer; however, EV-like cytopathic effects were not observed. HEp-2, L20B, VERO, and 293T cell lines did not provide an appropriate environment for EV-C116 growth. These results challenge the traditional recognition of the uncultured nature of EV-C116 strains and explain the difficulty of clinical detection.

Genetic characterization and molecular epidemiological analysis of enterovirus C99 in China

Enterovirus C99 (EV-C99) is a newly identified EV serotype within the species Enterovirus C. Few studies on EV-C99 have been conducted globally. More information and research on EV-C99 are needed to assess its genetic characteristics, phylogenetic relationships, and associations with enteroviral diseases. Here, the phylogenetic characteristics of 11 Chinese EV-C99 strains have been reported. The full-length genomic sequences of these 11 strains show 79.4-80.5% nucleotide identity and 91.7-94.3% amino acid (aa) identity with the prototype EV-C99. A maximum likelihood phylogenetic tree constructed based on the entire VP1 coding region identified 13 genotypes (A-M), revealing a high degree of variation among the EV-C99 strains. Phylogeographic analysis showed that the Xinjiang Uygur Autonomous Region is an important source of EV-C99 epidemics in various regions of China. Recombination analysis revealed inter-serotype recombination events of 16 Chinese EV-C99 strains in 5' untranslated regions and 3D regions, resulting in the formation of a single recombination form. Additionally, the Chinese strain of genotype J showed rich aa diversity in the P1 region, indicating that the genotype J of EV-C99 is still going through variable dynamic changes. This study contributes to the global understanding of the EV-C99 genome sequence and holds substantial implications for the surveillance of EV-C99.

Characterizing enterovirus C96 genome and phylodynamics analysis

Enterovirus C96 (EV-C96) is a recently discovered serotype belonging to enterovirus C species. It had been isolated from patients with acute flaccid paralysis, hand, foot, and mouth disease, diarrhea, healthy people, or environmental specimens. Despite increasing reports of the virus, the small number of full-length genomes available for EV-C96 has limited molecular epidemiological studies. In this study, newly collected rare EV-C96 strains in China from 1997 to 2020 were combined with sequences available in GenBank for comprehensive analyses. Sequence analysis revealed that the nucleotide sequence similarity of EV-C96 and the prototype strain (BAN00-10488) was 75%-81.8% and the amino acid sequence similarity was 85%-94.9%. EV-C96 had a high degree of genetic variation and could be divided into 15 genogroups. The mean evolutionary rate was 5.16 × 10-3 substitution/site/year, and the most recent common ancestor was dated to 1925. A recombination analysis revealed that EV-C96 may be a recombinant derived from other serotypes in the EV-C group in the nonstructural protein coding region. This comprehensive and integrated analysis of the whole genome sequence of EV-C96 provides valuable data for further studies on the molecular epidemiology of EV-C96 worldwide.

Molecular Epidemiology and Evolution of Coxsackievirus A14

As the proportion of non-enterovirus 71 and non-coxsackievirus A16 which proportion of composition in the hand, foot, and mouth pathogenic spectrum gradually increases worldwide, the attention paid to other enteroviruses has increased. As a member of the species enterovirus A, coxsackievirus A14 (CVA14) has been epidemic around the world until now since it has been isolated. However, studies on CVA14 are poor and the effective population size, evolutionary dynamics, and recombination patterns of CVA14 are not well understood. In this study, 15 CVA14 strains were isolated from HFMD patients in mainland China from 2009 to 2019, and the complete sequences of CVA14 in GenBank as research objects were analyzed. CVA14 was divided into seven genotypes A-G based on an average nucleotide difference of the full-length VP1 coding region of more than 15%. Compared with the CVA14 prototype strain, the 15 CVA14 strains showed 84.0-84.7% nucleotide identity in the complete genome and 96.9-97.6% amino acid identity in the encoding region. Phylodynamic analysis based on 15 CVA14 strains and 22 full-length VP1 sequences in GenBank showed a mean substitution rate of 5.35 × 10-3 substitutions/site/year (95% HPD: 4.03-6.89 × 10-3) and the most recent common ancestor (tMRCA) of CVA14 dates back to 1942 (95% HPD: 1930-1950). The Bayesian skyline showed that the effective population size had experienced a decrease-increase-decrease fluctuation since 2004. The phylogeographic analysis indicated two and three possible migration paths in the world and mainland China, respectively. Four recombination patterns with others of species enterovirus A were observed in 15 CVA14 strains, among which coxsackievirus A2 (CVA2), coxsackievirus A4 (CVA4), coxsackievirus A6 (CVA6), coxsackievirus A8 (CVA8), and coxsackievirus A12 (CVA12) may act as recombinant donors in multiple regions. This study has filled the gap in the molecular epidemiological characteristics of CVA14, enriched the global CVA14 sequence database, and laid the epidemiological foundation for the future study of CVA14 worldwide.

Enterovirus 71-Associated Infection in South Vietnam: Vaccination Is a Real Solution

Hand-foot-and-mouth disease (HFMD) is the most common enteroviral infection in South-East Asia. When evaluating the role of enterovirus 71 (EVA71) as an etiological agent of infectious disease in South Vietnam, we revealed a high proportion of EVA71 among identified species A enteroviruses found in 3542 samples from HFMD cases; 125 samples from cases of enteroviral meningitis; and 130 samples from acute flaccid paralysis (AFP) cases. These represent 50%, 54.8%, and 51.5%, respectively. According to molecular analysis, 90% of EVA71 were attributed to genotype C4 and 10% were attributed to genotype B5. The predominance of EVA71 circulation among the population proves the need to strengthen surveillance (with monitoring of enterovirus circulation for facilitation of HFMD outbreak prediction) and to increase the effectiveness of preventative measures by the implementation of vaccination against EVA71-associated infections. A phase III trial of a Taiwanese vaccine (EV71vac) in Taiwan and South Vietnam showed its safety, tolerability, and efficacy in children aged 2-71 months. This B4 genotype-based vaccine, which features cross-protection against B5 and C4 genotypes, and other existing EV71 vaccines can serve as a good approach to solving the HFMD problem, which is so important for Vietnam.

Genomic Epidemiology and Transmission Dynamics of Global Coxsackievirus B4

The aim of this study was to determine the global genetic diversity and transmission dynamics of coxsackievirus B4 (CVB4) and to propose future directions for disease surveillance. Next-generation sequencing was performed to obtain the complete genome sequence of CVB4, and the genetic diversity and transmission dynamics of CVB4 worldwide were analyzed using bioinformatics methods such as phylogenetic analysis, evolutionary dynamics, and phylogeographic analysis. Forty complete genomes of CVB4 were identified from asymptomatic infected individuals and hand, foot, and mouth disease (HFMD) patients. Frequent recombination between CVB4 and EV-B multiple serotypes in the 3D^pol region was found and formed 12 recombinant patterns (A-L). Among these, the CVB4 isolated from asymptomatic infected persons and HFMD patients belonged to lineages H and I, respectively. Transmission dynamics analysis based on the VP1 region revealed that CVB4 epidemics in countries outside China were dominated by the D genotype, whereas the E genotype was dominant in China, and both genotypes evolved at a rate of > 6.50 × 10^-3 substitutions/site/year. CVB4 spreads through the population unseen, with the risk of disease outbreaks persisting as susceptible individuals accumulate. Our findings add to publicly available CVB4 genomic sequence data and deepen our understanding of CVB4 molecular epidemiology.

Genetic characteristics of Coxsackievirus A6 from children with hand, foot and mouth disease in Beijing, China, 2017-2019

OBJECT

To investigate the evolution and genetic characteristics of Coxsackievirus A6 (CVA6) which acted as the predominant pathogen of hand, foot and mouth disease (HFMD) in children in Beijing, China, 2017-2019.

METHODS

Throat swab specimens were collected for general Enterovirus (EV), enterovirus A71 (EV-A71) and CVA16 detection by Real-time PCR. These general EV-positive samples were identified by semi-nested RT-PCR method and sequencing. The CVA6 VP1 gene and genome sequences were amplified and sequenced. The phylogenetic, variation and recombination analyses were performed.

RESULTS

A total of 1721 HFMD patients were enrolled in this study, with the male to female ratio of 1.62:1. The majority of cases were less than five years, which accounted for 73.50%. The overall detection rate of EV was 88.32% (1520/1721). A total of 8 EV types were identified, including CVA6 (55.86%), CVA16 (26.32%), EV-A71 (2.24%), CVA10 (2.04%), CVA4 (1.05%), CVA5 (0.59%), CVA2 (0.33%), and CVA8 (0.07%), while 175 (11.51%) EV were untyped. The main pathogen of HFMD was CVA6 from 2017 to 2018, while CVA6 and CVA16 were the main causative pathogens in 2019. The nucleotide and amino acid sequence identities of the 120 CVA6 complete VP1 gene sequences in this study were 91.2%-100.0% and 97.7%-100.0%, respectively. Compared with the prototype strain (Gdula) of CVA6, the nucleotide and amino acid sequence identities were 81.7%-84% and 94.7%-96.3%, respectively. The phylogenetic tree indicated that all 120 CVA6 sequences belonged to sub-genotype D3, while 119 CVA6 sequences belonged to evolutionary branch D3a, except one from 2017 belonged to D3b. Recombination analysis based on the complete genome sequences showed that potential multiple recombination may have occurred in 2B and 3D protein coding regions with EV-A114.

CONCLUSIONS

The main pathogens of HFMD were CVA6 and CVA16 in Beijing, China, 2017-2019; while these CVA6, as recombination strains, belonged to the D3a evolutionary branch.

Molecular Epidemiology Reveals the Co-Circulation of Two Genotypes of Coxsackievirus B5 in China

Coxsackievirus B5 (CVB5) is an important enterovirus B species (EV-Bs) type. We used the full-length genomic sequences of 53 viral sequences from the national hand, foot, and mouth disease surveillance network in the Chinese mainland (2001-2021). Among them, 69 entire VP1 coding region nucleotide sequences were used for CVB5 genotyping and genetic evolution analysis. Phylogenetic analysis based on a data set of 448 complete VP1 sequences showed that CVB5 could be divided into four genotypes (A-D) worldwide. Sequences from this study belonged to genotypes B and D, which dominated transmission in the Chinese mainland. Two transmission lineages of CVB5 have been discovered in the Chinese mainland, lineage 2 was predominant. Markov chain Monte Carlo analysis indicated that the tMRCA of CVB5 in the Chinese mainland could be traced to 1955, while the global trend could be traced to 1862, 93 years earlier than China. The evolution rate of CVB5 was higher in the Chinese mainland than worldwide. The spatiotemporal dynamics analysis of CVB5 assessed that virus transportation events were relatively active in Central, Northeast, North and Northwest China. Recombination analysis revealed frequent intertypic recombination in the non-structural region of CVB5 genotypes B and D with the other EV-Bs, revealing eight recombination lineages. Our study showed the molecular evolution and phylogeography of CVB5 that could provide valuable information for disease prevention.

Genomic Epidemiology and Phylodynamic Analysis of Enterovirus A71 Reveal Its Transmission Dynamics in Asia

Enterovirus A71 (EV-A71) is one of the main pathogens causing hand, foot, and mouth disease (HFMD) outbreaks in Asian children under 5 years of age. In severe cases, it can cause neurological complications and be life-threatening. In this study, 200 newly sequenced EV-A71 whole-genome sequences were combined with 772 EV-A71 sequences from GenBank for large-scale analysis to investigate global EV-A71 epidemiology, phylogeny, and Bayesian phylodynamic characteristics. Based on the phylogenetic analysis of the EV-A71 3D^pol region, six new evolutionary lineages (lineages B, J, K, O, P, and Q) were found in this study, and the number of evolutionary lineages was expanded from 11 to 17. Temporal dynamics and recombination breakpoint analyses based on genotype C revealed that recombination of nonstructural protein-coding regions, including 3D^pol, is an important reason for the emergence of new lineages. The EV-A71 epidemic in the Asia-Pacific region is complex, and phylogeographic analysis found that Vietnam played a key role in the spread of subgenotypes B5 and C4. The origin of EV-A71 subgenotype C4 in China is East China, which is closely related to the prevalence of subgenotype C4 in the south and throughout China. Selection pressure analysis revealed that, in addition to VP1 amino acid residues VP1-98 and VP1-145, which are associated with EV-A71 pathogenicity, amino acid residues VP1-184 and VP1-249 were also positively selected, and their functions still need to be determined by biology and immunology. This study aimed to provide a solid theoretical basis for EV-A71-related disease surveillance and prevention, antiviral research, and vaccine development through a comprehensive analysis. IMPORTANCE EV-A71 is one of the most important pathogens causing HFMD outbreaks; however, large-scale studies of EV-A71 genomic epidemiology are currently lacking. In this study, 200 new EV-A71 whole-genome sequences were determined. Combining these with 772 EV-A71 whole-genome sequences in the GenBank database, the evolutionary and transmission characteristics of global and Asian EV-A71 were analyzed. Six new evolutionary lineages were identified in this study. We also found that recombination in nonstructural protein-coding regions, including 3D^pol, is an important cause for the emergence of new lineages. The results provided a solid theoretical basis for EV-A71-related disease surveillance and prevention, antiviral research, and vaccine development.

Enteroviral Infections in Infants

Enteroviruses (EVs) are major pathogens in young infants. These viruses were traditionally classified into the following four subgenera: polio, coxsackie A and B, and echoviruses. Now that poliomyelitis seems to be controlled in most parts of the world, coxsackie and echoviruses are gaining more attention because (i) the structural and pathophysiological similarities and (ii) the consequent possibilities in translational medicine. Enteroviruses are transmitted mainly by oral and fecal-oral routes; the clinical manifestations include a viral prodrome including fever, feeding intolerance, and lethargy, which may be followed by exanthema; aseptic meningitis and encephalitis; pleurodynia; myopericarditis; and multi-system organ failure. Laboratory diagnosis is largely based on reverse transcriptase-polymerase chain reaction, cell culture, and serology. Prevention and treatment can be achieved using vaccination, and administration of immunoglobulins and antiviral drugs. In this article, we have reviewed the properties of these viruses, their clinical manifestations, and currently available methods of detection, treatment, and prognosis.

Description of non-polio enteroviruses identified in two national surveillance programmes in South Africa

Background

Human enteroviruses (EV) consist of 106 serotypes and four species: EV-A, EV-B, EV-C and EV-D. Enteroviruses cause clinical symptoms varying from severe to mild. Knowledge of EV burden in South Africa is limited, and as non-polio EV are important causes of acute flaccid paralysis (AFP) and meningitis, information on the circulating serotypes is vital.

Methods

Between 2010 and 2012, a total of 832 stool and viral isolate specimens were obtained from two national surveillance programmes at the National Institute for Communicable Diseases: the Rotavirus Sentinel Surveillance Programme (RSSP) and the AFP surveillance programme. Real-time polymerase chain reaction and Sanger sequencing were performed to detect and serotype EV.

Results

Non-polio EV were detected in 446 specimens, of which 308 were sequenced. Stool specimens yielded a greater variety of serotypes than viral cultures. EV-B viruses were predominant (58.44%), whilst EV-C viruses were detected in 31% of the specimens tested. South African prevalence for these viruses was higher than other countries, such as France with less than 2%, and Spain and the United States with less than 10%. The most common serotype detected was Enterovirus 99 (EV-C, 8.63%), which has not been reported in other regions.

Conclusion

Direct sequencing from stool specimens yields a broader, more comprehensive description of EV infections compared to sequencing from viral cultures. Disease-associated serotypes were detected, but only in small numbers. This study provides a baseline for EV strain circulation; however, surveillance needs to be expanded to improve EV knowledge in South Africa.

Coxsackievirus B4: an underestimated pathogen associated with a hand, foot, and mouth disease outbreak

In order to discover the causes of a coxsackievirus B4 (CV-B4)-associated hand, foot, and mouth disease (HFMD) outbreak and to study the evolutionary characteristics of the virus, we sequenced isolates obtained during an outbreak for comparative analysis with previously sequenced strains. Phylogenetic and evolutionary dynamics analysis was performed to examine the genetic characteristics of CV-B4 in China and worldwide. Phylogenetic analysis showed that CV-B4 originated from a common ancestor in Shandong. CV-B4 strains isolated worldwide could be classified into genotypes A-E based on the sequence of the VP1 region. All CV-B4 strains in China belonged to genotype E. The global population diversity of CV-B4 fluctuated substantially over time, and CV-B4 isolated in China accounted for a significant increase in the diversity of CV-B4. The average nucleotide substitution rate in VP1 of Chinese CV-B4 (5.20 × 10^-3 substitutions/site/year) was slightly higher than that of global CV-B4 (4.82 × 10^-3 substitutions/site/year). This study is the first to investigate the evolutionary dynamics of CV-B4 and its association with an HFMD outbreak. These findings explain both the 2011 outbreak and the global increase in CV-B4 diversity. In addition to improving our understanding of a major outbreak, these findings provide a basis for the development of surveillance strategies.

Epidemiological and molecular characteristics of circulating CVA16, CVA6 strains and genotype distribution in hand, foot and mouth disease cases in 2017 to 2018 from Western India

Hand, Foot, and Mouth disease (HFMD) is a mild exanthematous and febrile disease occurs in children aged ≤10 years old. The present study highlights clinical, epidemiological characteristics, distribution of enterovirus (EV) types, and sub genotypes in HFMD cases reported during 2017 to 2018 in Western India. A total of 93 clinical samples collected from 68 HFMD cases were included. The presence of EV-RNA was determined by 5'UTR based nested reverse transcription polymerase chain reaction followed by molecular typing, sub genotyping by VP1/2A junction or VP1, full VP1 gene amplification, and phylogenetic analysis. The study reports 80.64% (75/93) EV positivity and 94.66% (71/75) typing rate, with a predominant circulation of CVA16 and CVA6 strains. Sequence analysis revealed the presence of coxsackievirus (CV)A16 (57.7%), CVA6 (40.8%), and Echo1 (1.4%) strains. EV infections were predominantly observed in children aged 1 to 3 years old (43.9%). Although cases were reported throughout the year, peaked in July (15.8%) and August (24.6%) months and persisted till September (19.3%). All the CVA16 and CVA6 positive strains were genotyped using full VP1 gene amplification. All CVA16 Indian strains (n = 41) were clustered with rarely reported B1c sub genotype and CVA6 strains (n = 29) with E2 sub-lineage. The study highlights the genetic characteristics of circulating CVA16, CVA6, and Echo1 strains in HFMD cases from Western India. The emergence of CVA16 B1c genotype and sub-lineage E2 of CVA6 strains and their constant circulation further demands systemic surveillance studies on HFMD from different parts of India to facilitate the rapid diagnosis of CVA16 and CVA6 strains using the molecular and serological based approach and for intervention strategies.

Direct Diagnosis of Echovirus 12 Meningitis Using Metagenomic Next Generation Sequencing

The current point-of-care diagnosis of enterovirus meningitis does not identify the viral genotype, which is prognostic. In this case report, more than 81% of an Echovirus 12 genome were detected and identified by metagenomic next-generation sequencing, directly from the cerebrospinal fluid collected in a 6-month-old child with meningeal syndrome and meningitis: introducing Echovirus 12 as an etiological agent of acute meningitis in the pediatric population.

A novel interspecies recombinant enterovirus (Enterovirus A120) isolated from a case of acute flaccid paralysis in China

EV-A120 is a recently identified serotype of the enterovirus A species. Only one full-length genomic sequence is currently available in GenBank, and very few studies have been conducted on EV-A120 globally. Thus, additional information and research on EV-A120 are needed to explore its genetic characteristics, phylogeny, and relationship with enteroviral disease. In this study, we report the phylogenetic characteristics of a EV-A120 strain (Q0082/XZ/CHN/2000) from Tibet, China. The amino acid sequence similarity and nucleotide sequence similarity of the full-length genomic sequence of this EV-A120 strain and the EV-A120 prototype strain were 96.3% and 79.9%, respectively, showing an evolutionary trend. Recombination analysis found intraspecies recombination in the 5′ -UTR, 2B, 2C, and 3D regions. Serum neutralization testing of the EV-A120 (Q0082) strain was also carried out. Low serum-positive rates and geometric mean titres (GMTs) indicated that the extent of EV-A120 transmission and exposure in the population was very limited compared with that in the outbreaks of EV-A71 and CV-A16 in China since 2008. The EV-A120 strain (Q0082) is non-temperature sensitive, indicating its potential to spread in the population. In summary, this study reports the full-length genomic sequence of EV-A120 and provides important information for its global molecular epidemiology.

High Heterogeneity of Echoviruses in Brazilian Children with Acute Gastroenteritis

Echoviruses (E) are a diverse group of viruses responsible for various pathological conditions in humans including aseptic meningitis, myocarditis, and acute flaccid paralysis. The detection and identification of echovirus genotypes in clinical samples is challenging due to its high genetic diversity. Here, we report the complete genome sequences of nine echoviruses, obtained by next-generation sequencing of 238 fecal samples from individuals with gastroenteritis in regions of Brazil. Detected viruses were classified into six genotypes: Three E1 sequences (BRA/TO-028, BRA/TO-069 and BRA/TO-236), one E3 (BRA/TO-018), one E11 (BRA/TO-086), one E20 (BRA/TO-016), two E29 (BRA/TO-030 and BRA/TO-193), and one E30 sequence (BRA/TO-032). Phylogenetic analysis indicated that the echoviruses E1 and E29 circulating in Brazil are divergent from strains circulating worldwide. The genotype diversity identified in our study may under-represent the total echovirus diversity in Brazil because of the small sample size and the restricted geographical distribution covered by the survey.

One-Step Reverse-Transcription Recombinase Polymerase Amplification Using Lateral Flow Strips for the Detection of Coxsackievirus A6

Hand, foot, and mouth disease (HFMD) is a common infectious disease affecting mainly children under 5 years of age. Coxsackievirus A6 (CVA-6), a major causative pathogen of HFMD, has caused outbreaks in recent years. Currently, no effective vaccine or antiviral treatments are available. In this study, one-step reverse-transcription recombinase polymerase amplification (RT-RPA), combined with a disposable lateral flow strip (LFS) assay, was developed to detect CVA-6. This assay can be performed in less than 35 min at 37°C without expensive instruments, and the result can be observed directly with the naked eye. The sensitivity of the RT-RPA-LFS was 10 copies per reaction, which was comparable to that of the conventional real-time quantitative polymerase chain reaction (qPCR) assays. Moreover, the assay specificity was 100%. The clinical performance of the RT-RPA-LFS assay was evaluated using 142 clinical samples, and the coincidence rate between RT-RPA-LFS and qPCR was 100%. Therefore, our RT-RPA-LFS assay provides a simple and rapid approach for point-of-care CVA-6 diagnosis.

Molecular Epidemiological, Serological, and Pathogenic Analysis of EV-B75 Associated With Acute Flaccid Paralysis Cases in Tibet, China

Enterovirus B75 (EV-B75) is a newly identified serotype of the enterovirus B species. To date, only 112 cases related to EV-B75 have been reported worldwide, and research on EV-B75 is still limited with only two full-length genome sequences available in GenBank. The present study reported seven EV-B75 sequences from a child with acute flaccid paralysis and six asymptomatic close contacts in Shigatse, Tibet. Phylogenetic analysis revealed that the Tibetan strain was possibly imported from neighboring India. Seroepidemiological analyses indicated that EV-B75 has not yet caused a large-scale epidemic in Tibet. Similarity plots and boot scanning analyses revealed frequent intertypic recombination in the non-structural region of all seven Tibet EV-B75 strains. All seven Tibetan strains were temperature-sensitive, suggesting their poor transmissibility in the environment. Overall, though the seven Tibetan strains did not cause large-scale infection, prevention and control of the novel enterovirus cannot be underestimated.

Prevalence and Bayesian Phylogenetics of Enteroviruses Derived From Environmental Surveillance Around Polio Vaccine Switch Period in Shandong Province, China

We present the results of environmental surveillance for poliovirus (PV) and non-poliovirus (NPEV) around the switch from trivalent to bivalent oral polio-vaccine (OPV) which occurred in China in May 2016. Sewage samples were collected in Jinan and Linyi city from 2015 to 2017. Enterovirus (EV) isolation, VP1 amplification, Sanger sequencing, and phylogenetic analyses were performed. Among105 sewage samples (36 in Jinan and 69 in Linyi), 101 were positive for EV, with 74.3% (78/105) PV-positive samples and 90.5% (95/105) NPEV-positive samples. A total of 893 EV isolates were obtained, including 326 (36.5%) PVs and 567 (63.5%) NPEVs. Echovirus (E) -11 was the most common serotype out of 18 detected NPEV types (120/567), followed by E-3 (75/567) and E-6 (74/567). PV2 vanished and PV3 came to be the ascendant PV type in sewage after May 2016. Eight PV isolates were judged as pre-vaccine-derived poliovirus (pre-VDPV) and no VDPV or wild PV isolates were monitored. Bayesian phylogenetics demonstrated global E-11 originated in 1876 and evolved with the estimated rate of 4.63 × 10^-3 nucleotide substitutions per site per year (s/s/y). Multiple circulating clusters that originated at different times were coexisting in Shandong province. The most recently common ancestor of global coxsackievirus B5 could date back to 1867, at the evolutionary rate of 3.95 × 10^-3 s/s/y. In conclusion, our study described the changes of PVs and NPEVs around the polio vaccine switch period and provided meaningful global molecular epidemiological data for further studies of EV-related diseases among the population.

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.

Cell culture demonstrates superior sensitivity over one step real time RT PCR and nested VP1 amplification for Enteroviruses

This study evaluated and compared the sensitivity profile of routine cell culture, nested VP1 amplification and one step real time RT PCR for Enteroviruses. Serially diluted spiked samples of four model viruses (EV71, CVA16, CVB5 and PV1) and 32 true positive samples including Poliovirus (PV1 & PV3), Coxsackie virus (CVB5, CVB3, CVB1 & CVA4, 10, 16), Echovirus (Echo 6, 7, 11, 13, 18, 25 & 30) and Enterovirus 71 (E71), and 32 true negative stool samples were subjected to cell culture, nested RT PCR and one step real time RT PCR. The result of sensitivity test indicated superior sensitivity with one step real time RT PCR (75 %, 24/32) against cell culture (71.9 %, 23/32) and nested RT PCR (65.6 %, 21/32). The most specific test was cell culture (100 %, 32/32), followed by nested RT PCR (96.9 %, 31/32). Positive predictive values were 100 %: 23/23, 95.5 %; 21/22 and 88.9 %; 24/27, for cell culture, nested RT PCR and one step real time RT PCR, respectively, and one step real time RT PCR had the highest negative predictive value (78.4 %, 29/37). Overall result indicate relatively high analytical sensitivity with all the tests, suggesting superior performance by cell culture. Therefore, cell culture is the gold standard. However, considering intensive nature of cell cultures and prolong window for results, it is wise to consider one step real time RT PCR in routine diagnosis for its added advantages. Meanwhile, selecting a combination of tests can maximize detection, depending on the laboratory strength.

Molecular typing and characterization of a novel genotype of EV-B93 isolated from Tibet, China

EV-B93 is a novel serotype within the Enterovirus B species and is uncommon worldwide. Currently, only one full-length genomic sequence (the prototype strain) has been deposited in the GenBank database. In this study, three EV-B93 were identified, including one from an acute flaccid paralysis (AFP) patient (named 99052/XZ/CHN/1999, hereafter XZ99052) and two from healthy children (named 99096/XZ/CHN/1999 and 99167/XZ/CHN/1999, hereafter XZ99096 and XZ99167, respectively) from Tibet in 1999 during the polio eradication program. The identity between the nucleotide and amino acid sequences of the Tibet EV-B93 strain and the EV-B93 prototype strain is 83.2%–83.4% and 96.8%–96.9%, respectively. The Tibet EV-B93 strain was found to have greater nucleotide sequence identity in the P3 region to another enterovirus EV-B107 as per a phylogenetic tree analysis, which revealed that recombination occurred. Seroepidemiology data showed that EV-B93 has not produced an epidemic in Tibet and there may be susceptible individuals. The three Tibet EV-B93 strains are temperature-resistant with prognosticative virulence, suggesting the possibility of a potential large-scale outbreak of EV-B93. The analyzed EV-B93 strains enrich our knowledge about this serotype and provide valuable information on global EV-B93 molecular epidemiology. What is more, they permit the appraisal of the serotype's potential public health impact and aid in understanding the role of recombination events in the evolution of enteroviruses.

Neonatal hand, foot, and mouth disease due to coxsackievirus A6 in Shanghai

BACKGROUND

Evidence of hand, foot, and mouth disease (HFMD) in neonates is limited. The aim of this study was to evaluate the clinical symptoms, pathogens, possible transmission routes, and prognosis of neonatal HFMD in Shanghai.

METHODS

This was a case-control study based on the HFMD registry surveillance system. All neonates and infected family members were enrolled between 2016 and 2017 in Shanghai. Neonates with HFMD were followed for at least half a year. Detailed questionnaires, medical history, and physical examination were recorded. Routine blood examination, liver and renal function, immunophenotypes of peripheral blood lymphocytes (CD3, CD4, and CD8 T-cells; NK cells), immunoglobulin (Ig) M, IgG, and IgA, and cytokine interleukin (IL-1β, IL-2R, IL-6, IL-8, IL-10, and TNF-α) levels were measured. All rectal swab specimens were collected and genotyped for enterovirus, and phylogenetic analysis based on the VP1 sequences of coxsackievirus A6 (CV-A6) was performed to investigate molecular and evolutionary characteristics. T-test or nonparametric test was used to evaluate the differences. Logistic analysis was applied to calculate the risk of clinical manifestations in the group of HFMD neonates and their paired siblings.

RESULTS

There were 16 neonates among the 12,608 diagnosed patients with HFMD, accounting for 0.13%. All neonatal infections were transmitted by other members of the family, mainly the elder siblings, and were caused by CV-A6. CV-A6 was the emerging and predominant causative agent of HFMD in Shanghai. None of the neonates with HFMD experienced fever, onychomadesis, or severe complications. However, two elder sibling patients showed lethargy, and one developed hypoperfusion. In the elder siblings with HFMD, the proportion of white blood cells was generally higher than in neonates with HFMD. The immunologic function of the neonates with HFMD was basically normal. The levels of inflammatory markers were higher in both neonates and elder siblings with HFMD compared to age-matched controls. The clinical symptoms receded about 1 week after onset. None of the neonates had sequelae.

CONCLUSIONS

In our study, CV-A6 infection in neonates was benign, but had the character of family clustering. Due to the two-child policy in China, elder siblings may be the main route of HFMD transmission.

Detection of enteroviruses in urban sewage by next generation sequencing and its application in environmental surveillance

Environmental surveillance has been used successfully in monitoring enterovirus (EV) circulation; however cell culture method may introduce a selective bias in those EV strains that are recovered from the environment. In this study, urban sewage samples were collected monthly in Jinan, China in 2018 and concentrated via membrane adsorption/elution method. A P1 seminested RT-PCR (RT-snPCR) and NGS method was developed, by which amplicons of 4000 nucleotide in length covering the entire P1 region of EVs were obtained from sewage concentrates and were further analyzed by next generation sequencing (NGS). In addition, for each sewage concentrate, two other assays - cell culture and NGS based partial VP1 amplicon sequencing - were conducted in parallel and compared. The results showed that the P1 RT-snPCR and NGS method generated the most data, with 32 serotypes identified belonging to species EV-A (n = 11), EV-B (n = 14), and EV-C (n = 7). These serotypes covered all those detected from the methods of cell culture (n = 10) and partial VP1 amplicon sequencing (n = 16). EV serotypes from acute flaccid paralysis surveillance correlated with those from sewage. Phylogenetic analysis on coxsackievirus B5, a common pathogen of meningitis, revealed close genetic relationship between environmental and clinical sequences. These results demonstrate sewage contains different EVs related to a variety of diseases. Traditional cell culture method underestimates the existence of some serotypes. NGS based environmental surveillance provides data which are consistent with those from clinical diseases, greatly improves our understanding on the actual circulation in the population, and should be encouraged for public health surveillance.

Whole-Genome Sequencing of Human Enteroviruses from Clinical Samples by Nanopore Direct RNA Sequencing

Enteroviruses are small RNA viruses that affect millions of people each year by causing an important burden of disease with a broad spectrum of symptoms. In routine diagnostic laboratories, enteroviruses are identified by PCR-based methods, often combined with partial sequencing for genotyping. In this proof-of-principle study, we assessed direct RNA sequencing (DRS) using nanopore sequencing technology for fast whole-genome sequencing of viruses directly from clinical samples. The approach was complemented by sequencing the corresponding viral cDNA via Illumina MiSeq sequencing. DRS of total RNA extracted from three different enterovirus-positive stool samples produced long RNA fragments, covering between 59% and 99.6% of the most similar reference genome sequences. The identification of the enterovirus sequences in the samples was confirmed by short-read cDNA sequencing. Sequence identity between DRS and Illumina MiSeq enterovirus consensus sequences ranged between 94% and 97%. Here, we show that nanopore DRS can be used to correctly identify enterovirus genotypes from patient stool samples with high viral load and that the approach also provides rich metatranscriptomic information on sample composition for all life domains.

Environmental Surveillance Complements Case-Based Surveillance of Acute Flaccid Paralysis in Polio Endgame Strategy 2019-2023

The Polio Endgame Strategy 2019-2023 has been developed. However, more effective and efficient surveillance activities should be conducted with the preparedness of emergence for vaccine-derived poliovirus (VDPV) or wild poliovirus (WPV). We reviewed the impact of the case-based acute flaccid paralysis (AFP) surveillance (1991 to 2018) and environmental surveillance (2011 to 2018) in polio eradication in Shandong province of China. Clinical characteristics of AFP cases and enterovirus (EV) investigation of research samples were assessed. During the period, 10,224 AFP cases were investigated, and 352 sewage samples were collected. The nonpolio AFP rate sustained at over 2.0/100,000 since 1997. Of 10,224 cases, males and young children experienced a higher risk of severe diseases, and 68.5% suffered lower limb paralysis. We collected 1,707 EVs from AFP cases, including 763 polioviruses and 944 nonpolio enteroviruses (NPEVs). No WPV was isolated since 1992. The AFP surveillance showed high sensitivity in detecting 143 vaccine-associated paralytic poliomyelitis (VAPP) cases and 6 VDPVs. For environmental surveillance, 217 (61.6%) samples were positive for poliovirus, and altogether, 838 polioviruses and 2,988 NPEVs were isolated. No WPV was isolated in environmental surveillance, although one VDPV2 was identified. Phylogenetic analysis revealed environmental surveillance had the capacity to detect a large scope of NPEVs. The case-based AFP surveillance will be indispensable for detecting VAPP cases and VDPV circulation in countries using oral polio vaccine. Environmental surveillance is advantageous in identifying EV circulation and responding to ongoing circulating VDPV outbreaks and should be expanded to complement the AFP surveillance.IMPORTANCE Interrupting wild poliovirus transmission and stopping circulating vaccine-derived poliovirus (cVDPV) outbreaks have been proposed as two global goals by the World Health Organization in the Global Polio Eradication Initiative (GPEI). This analysis, based on the 28-year acute flaccid paralysis (AFP) surveillance and 8-year environmental surveillance, provides continued high-quality surveillance performance in achieving the GPEI and detecting the circulation of enterovirus. Given the ongoing cVDPV outbreaks in the world, we present the surveillance capacity of environmental surveillance in capturing enterovirus circulation. The final poliovirus (especially VDPV) elimination has become increasingly complex, and the case-based AFP surveillance alone will lead to difficulties in early detecting dynamics of poliovirus transmission and monitoring the extent of environmental circulation. This study goes beyond previous work to provide a detailed comprehensive evaluation of the enterovirus surveillance and can be used to formulate a set of implementation plan and performance indicators for environmental surveillance.

A Large-Scale Outbreak of Echovirus 30 in Gansu Province of China in 2015 and Its Phylodynamic Characterization

Background

Echovirus 30 (E-30) has been investigated and reported worldwide and is closely associated with several infectious diseases, including encephalitis; myocarditis; and hand, foot, and mouth disease. Although many E-30 outbreaks associated with encephalitis have been reported around the world, it was not reported in northwest China until 2015.

Methods

The clinical samples, including the feces, serum, throat swabs, and cerebrospinal fluid, were collected for this study and were analyzed for diagnosis. E-30 was isolated and processed according to the standard procedures. The epidemiological and phylogenetic analysis were performed to indicate the characteristics of E-30 outbreaks and phylodynamics of E-30 in China.

Results

The E-30 outbreaks affected nine towns of Gansu Province in 2015, starting at a school of Nancha town and spreading to other towns within 1 month. The epidemiological features showed that children aged 6–15 years were more susceptible to E-30 infection. The genotypes B and C cocirculated in the world, whereas the latter dominated the circulation of E-30 in China. The genome sequences of this outbreak present 99.3–100% similarity among these strains, indicating a genetic-linked aggregate outbreak of E-30 in this study. Two larger genetic diversity expansions and three small fluctuations of E-30 were observed from 1987 to 2016 in China, which revealed the oscillating patterns of E-30 in China. In addition, the coastal provinces of China, such as Zhejiang, Fujian, and Shandong, were initially infected, followed by other parts of the country. The E-30 strains isolated from mainland of China may have originated from Taiwan of China in the last century.

Conclusion

The highly similar E-30 genomes in this outbreak showed an aggregate outbreak of E-30, with nine towns affected. Our results suggested that, although the genetic diversity of E-30 oscillates, the dominant lineages of E-30 in China has complex genetic transmission. The coastal provinces played an important role in E-30 spread, which implied further development of effective countermeasures. This study provides a further insight into the E-30 outbreak and transmission and illustrates the importance of valuable surveillance in the future.

Phylogenetic characteristics and molecular epidemiological analysis of novel enterovirus EV-B83 isolated from Tibet, China

Enterovirus B83 (EV-B83) is a new member of the enterovirus B group. Currently, there are only two full-length genomic sequences of EV-B83 in the GenBank database and few VP1 region sequences. The aetiology and epidemiology of EV-B83 is unclear. 24 stool specimens were collected from twelve AFP patients and 298 stool specimens were collected from 298 healthy children in support of polio eradication activities in Tibet in 1999. Two polioviruses (isolated by L20B cell) and one non-polio enterovirus (isolated by RD cell) were isolated from AFP patients and nine polioviruses (isolated by L20B cell) and 90 non-polio enteroviruses (isolated by RD cell) were isolated from health children. Through molecular typing, we confirmed that the six of non-polio enteroviruses belong to EV-B83. The sequence similarity between the VP1 region of the Tibet isolates and that of the EV-B83 prototype strain was 80%. The maximum-likelihood phylogenetic tree of the partial VP1 region in EV-B83 demonstrated that EV-B83 formed four genotypes globally during the evolution process. The six Tibet EV-B83 strains formed the D genotype alone. Recombination analysis of Tibet EV-B83 showed that CV-B4, CV-A9, EV-B80, and EV-B106 may act as recombinant donors in multiple regions. The serum neutralization test showed that the antibody-positive rate was 58.8% and GMT was 1:19.70, which was higher than the previously reported results of EV-B106 and EV-B80. Temperature sensitivity test results showed that the six Tibet EV-B83 strains were temperature-insensitive with stronger virulence and potential infectivity, which was consistent with the results of the serum neutralization test. This study enriched the genome-wide sequence, epidemiological characteristics, and provided basic data for the follow-up study of EV-B83.

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.

Pathogenic characteristics of hand, foot and mouth disease in Shaanxi Province, China, 2010-2016

Hand, foot, and mouth disease (HFMD) is a common childhood illness caused by enteroviruses. We analyzed the pathogenic characteristics of HFMD in Shaanxi province, China, during 2010-2016. Clinical samples were collected from HFMD cases. Real-time PCR and RT-PCR were used to identify the enterovirus(EVs) serotypes. Viral RNA sequences were amplified using RT-PCR and compared by phylogenetic analysis. Descriptive epidemiological methods were used to analyze. A total of 16,832 HFMD positive cases were confirmed in the laboratory. EV-A71 and CV-A16 were the main pathogens in 2010. EV-A71 was the dominant pathogen in the periods of 2011 to 2012 and 2014, 2016. In 2013 and 2015, other EVs increased greatly, in which CV-A6 was the predominant pathogen. EV-A71 was more frequently detected in deaths and severe cases. Phylogenetic analysis revealed that EV-A71 belonged to the C4a evolution branch of C4 sub-genotype and CV-A16 belonged to the B1a or B1b evolution branch of B1 sub-genotype, whereas CV-A6 strains were assigned to D2 or D3 sub-genotype. The pathogen spectrum of HFMD has changed in 7 years, and the major serotypes EV-A71, CV- A16 and CV- A6 alternated or co-circulated. Long-term surveillance and research of EVs should be strengthened for the prevention and control of HFMD.

Environmental Surveillance Can Dynamically Track Ecological Changes in Enteroviruses

Environmental surveillance can be used to trace enteroviruses shed from human stool using a sewer network that is independent of symptomatic or asymptomatic infection. In this study, the local transmission of enteroviruses was analyzed using two wastewater treatment plants, which were relatively close to each other (15 km), designated as sentinels. Influent was collected at both sentinels once a month from 2013 to 2016, and viruses were isolated. Using neutralizing tests with type-specific polyclonal antisera and molecular typing, 933 isolates were identified as enteroviruses. Our results showed that the frequency of virus isolation varied for each serotype at the two sentinels in a time-dependent manner. Because echovirus 11 (Echo11) and coxsackievirus B5 isolates showed a high frequency and were difficult to distinguish, they were further grouped into various lineages based on the VP1 amino acid sequences. The prevalence of each lineage was visualized using multidimensional scaling. The results showed that Echo11 isolates of the same lineage were isolated continuously, similar to coxsackievirus B5 isolates of three lineages. Conversely, Echo1, Echo13, Echo18, Echo19, Echo20, Echo29, and Echo33 were isolated only once each. Our findings suggested that if an enterovirus is imported into the population, it may result in small-scale transmission, whereas if there are initially many infected individuals, it may be possible for the virus to spread to a wide area, beyond the local community, over time. In addition, our findings could provide insights into risk assessment of transmission for importation of poliovirus in polio-free countries and regions.IMPORTANCE In this study, we showed that environmental enterovirus surveillance can be used to monitor the propagation of nonpolio enteroviruses in addition to poliovirus detection. Since epidemiological studies of virus transmission based on the past were performed using specimens from humans, there were limitations to research design, such as specimen collection for implementation on a large-scale target population. However, environmental monitoring can dynamically track the ecological changes in enteroviruses in the region by monitoring viruses in chronological order and targeting the population within the area by monitoring viruses over time. We observed differences in the transmission of echovirus 11 and coxsackievirus B5 in the region according to lineage in a time-dependent manner and with a multidimensional scaling pattern.

Rapid and Cost-Efficient Enterovirus Genotyping from Clinical Samples Using Flongle Flow Cells

Enteroviruses affect millions of people worldwide and are of significant clinical importance. The standard method for enterovirus identification and genotyping still relies on Sanger sequencing of short diagnostic amplicons. In this study, we assessed the feasibility of nanopore sequencing using the new flow cell “Flongle” for fast, cost-effective, and accurate genotyping of human enteroviruses from clinical samples. PCR amplification of partial VP1 gene was performed from multiple patient samples, which were multiplexed together after barcoding PCR and sequenced multiple times on Flongle flow cells. The nanopore consensus sequences obtained from mapping reads to a reference database were compared to their Sanger sequence counterparts. Using clinical specimens sampled over different years, we were able to correctly identify enterovirus species and genotypes for all tested samples, even when doubling the number of barcoded samples on one flow cell. Average sequence identity across sequencing runs was >99.7%. Phylogenetic analysis showed that the consensus sequences achieved with Flongle delivered accurate genotyping. We conclude that the new Flongle-based assay with its fast turnover time, low cost investment, and low cost per sample represents an accurate, reproducible, and cost-effective platform for enterovirus identification and genotyping.

The rare enterovirus c99 and echovirus 29 strains in Brazil: potential risks associated to silent circulation

Human enteroviruses (EVs) are associated with a wide spectrum of human diseases. Here we report the complete genome sequences of one EV-C99 strain and one E29 strain obtained from children suffering from acute gastroenteritis, without symptoms of enteroviral syndromes. This is the first report of EV-C99 in South America, and the second E29 genome described worldwide. Continuous surveillance on EVs is vital to provide further understanding of the circulation of new or rare EV serotypes in the country. The present study also highlights the capacity of EVs to remain in silent circulation in populations.

Multiple genotypes of Echovirus 11 circulated in mainland China between 1994 and 2017

Echovirus 11 (E-11) is one of the most frequently isolated enteroviruses causing meningitis and other diseases such as hand, foot, and mouth disease (HFMD) and acute flaccid paralysis (AFP). Fifty-nine newly determined E-11 VP1 sequences from the China AFP and HFMD surveillance network and 500 E-11 VP1 sequences obtained from the GenBank database, which were associated with 12 categories of diseases, were screened for phylogenetic analysis. Based on the standard method of genotype classification, E-11 strains circulated worldwide were reclassified into six genotypes as A, B, C, D, E, and F, in which genotype F is newly divided, and genotypes A and C are further divided into A1–5 and C1–4 by this research, whereas genotype D was still divided into D1–5 as in a previous study of Oberste et al. Sub-genotype A1 was the predominant sub-genotype in mainland China between 2008–2017, whereas sub-genotype D5 was the predominant sub-genotype circulated outside China from 1998–2014. However, genotype and sub-genotype spectra showed statistical significance among AFP and HFMD cases (χ² = 60.86, P < 0.001), suggesting that different genotypes might have a tendency to cause different diseases. Strengthening the surveillance of E-11 might provide further information about pathogenic evolution or specific nucleotide mutation associated with different clinical diseases.

The surveillance of the epidemiological and serotype characteristics of hand, foot, mouth disease in Neijiang city, China, 2010-2017: A retrospective study

Hand, foot, and mouth disease (HFMD) is well recognized as one of the major threats to children’s health globally. The increasing complexity of the etiology of HFMD still challenges disease control in China. There is little surveillance of the molecular epidemiological characteristics of the enteroviruses (EVs) that cause HFMD in Neijiang city or the Sichuan Basin area in Southwest China. In this study, demographic and epidemiological information for 14,928 probable HFMD cases was extracted and analyzed to describe the epidemic features of HFMD in Neijiang city from Jan 2010 to Dec 2017. The swab samples of select probable HFMD cases from 2012 to 2017 were tested by reverse transcription (RT) real-time PCR to identify the serotype distribution of EVs, and 110 randomly selected RT-real-time PCR positive samples were then amplified and analyzed for the VP1 or VP4 regions of EVs to further analyze the phylogenetic characteristics of the circulating strains in this area. The eight-year average annual incidence was 49.82 per 100,000 in Neijiang. The incidence rates varied between 19.51 and 70.73 per 100,000, demonstrating peaks of incidence in even-number years (2012, 2014 and 2016). The median age of the probable cases was 27 months and the interquartile range (25^th to 75^th percentile) of ages for the probable HFMD cases was between 14 and 42 months. The male-to-female ratio of the probable HFMD cases was 1.47:1, and scattered children were the major population classification (81.7%). Two epidemic peaks were observed: one major peak between April and July and the other lesser peak between October and December. Of 6513 probable cases tested with RT-real-time PCR, 4015 (61.6%) were positive for enterovirus with the serotype distribution as follows: EV71+, 30.1% (n = 1210); CV-A16+, 28.7% (n = 1154) and a sole pan-enterovirus+, 41.1% (n = 1651). A total of 91 cases (82.7%, 91/110) were successfully amplified and underwent phylogenetic analysis: all EV71+ cases were C4a serotype (n = 23/30); all CV-A16+ cases were B2b serotype (n = 24/30); of 42 sole pan-enterovirus+ samples, 20 were CV-A6, 14 were CV-A10 and the rest within this group were CV-A4 (n = 4), CV-A8 (n = 2), CV-A9 (n = 1) and CV-B3 (n = 1). Our findings provide important evidence that aids the improvement of strategies for vaccination against HFMD and comprehensive disease control in China.

Two Coxsackievirus B3 outbreaks associated with hand, foot, and mouth disease in China and the evolutionary history worldwide

BACKGROUND

Coxsackievirus B3 (CV-B3) is usually associated with aseptic meningitis and myocarditis; however, the association between CV-B3 and hand, foot, and mouth disease (HFMD) has not been clearly demonstrated, and the phylogenetic dynamics and transmission history of CV-B3 have not been well summarized.

METHOD

Two HFMD outbreaks caused by CV-B3 were described in Hebei Province in 2012 and in Shandong Province in 2016 in China. To analyze the epidemiological features of two CV-B3 outbreaks, a retrospective analysis was conducted. All clinical specimens from CV-B3 outbreaks were collected and disposed according to the standard procedures supported by the WHO Global Poliovirus Specialized Laboratory. EV genotyping and phylogenetic analysis were performed to illustrate the genetic characteristics of CV-B3 in China and worldwide.

RESULTS

Two transmissible lineages (lineage 2 and 3) were observed in Northern China, which acted as an important "reservoir" for the transmission of CV-B3. Sporadic exporting and importing of cases were observed in almost all regions. In addition, the global sequences of CV-B3 showed a tendency of geographic-specific clustering, indicating that geographic-driven adaptation plays a major role in the diversification and evolution of CV-B3.

CONCLUSIONS

Overall, our study indicated that CV-B3 is a causative agent of HFMD outbreak and revealed the phylogenetic dynamics of CV-B3 worldwide, as well as provided an insight on CV-B3 outbreaks for effective intervention and countermeasures.

Genetic landscape and macro-evolution of co-circulating Coxsackieviruses A and Vaccine-derived Polioviruses in the Democratic Republic of Congo, 2008-2013

Enteroviruses (EVs) are among the most common viruses infecting humans worldwide but only a few Non-Polio Enterovirus (NPEV) isolates have been characterized in the Democratic Republic of Congo (DR Congo). Moreover, circulating vaccine-derived polioviruses (PVs) [cVDPVs] isolated during multiple outbreaks in DR Congo from 2004 to 2018 have been characterized so far only by the sequences of their VP1 capsid coding gene. This study was carried to i) investigate the circulation and genetic diversity of NPEV and polio vaccine isolates recovered from healthy children and Acute Flaccid Paralysis (AFP) patients, ii) evaluate the occurrence of genetic recombination among EVs belonging to the Enterovirus C species (including PVs) and iii) identify the virological factors favoring multiple emergences of cVDPVs in DR Congo. The biological material considered in this study included i) a collection of 91 Sabin-like PVs, 54 cVDPVs and 150 NPEVs isolated from AFP patients between 2008 and 2012 in DR Congo and iii) a collection of 330 stool specimens collected from healthy children in 2013 in the Kasai Oriental and Maniema provinces of DR Congo. Studied virus isolates were sequenced in four distinct sub-genomic regions 5’-UTR, VP1, 2C^ATPase and 3D^pol. Resulting sequences were compared through comparative phylogenetic analyses. Virus isolation showed that 19.1% (63/330) healthy children were infected by EVs including 17.9% (59/330) of NPEVs and 1.2% (4/330) of type 3 Sabin-like PVs. Only one EV-C type, EV-C99 was identified among the NPEV collection from AFP patients whereas 27.5% of the 69 NPEV isolates typed in healthy children belonged to the EV-C species: CV-A13 (13/69), A20 (5/69) and A17 (1/69). Interestingly, 50 of the 54 cVDPVs featured recombinant genomes containing exogenous sequences in at least one of the targeted non-structural regions of their genomes: 5’UTR, 2C^ATPase and 3D^pol. Some of these non-vaccine sequences of the recombinant cVDPVs were strikingly related to homologous sequences from co-circulating CV-A17 and A20 in the 2C^ATPase region as well as to those from co-circulating CV-A13, A17 and A20 in the 3D^pol region. This study provided the first evidence uncovering CV-A20 strains as major recombination partners of PVs. High quality AFP surveillance, sensitive environmental surveillance and efficient vaccination activities remain essential to ensure timely detection and efficient response to recombinant cVDPVs outbreaks in DR Congo. Such needs are valid for any epidemiological setting where high frequency and genetic diversity of Coxsackieviruses A13, A17 and A20 provide a conducive viral ecosystem for the emergence of virulent recombinant cVDPVs.

The strategy of the Global Polio Eradication Initiative is based on the surveillance of patients suffering from Acute Flaccid Paralysis (AFP) and mass vaccination with live-attenuated vaccine strains of polioviruses (PVs) in endemic areas. However, vaccine strains of PVs can circulate and replicate for a long time when the vaccine coverage of the population is low. Such prolonged circulation and replication of vaccine strains of PVs can result to the emergence of circulating vaccine-derived polioviruses [cVDPVs] that are as virulent as wild PVs. In this study, we performed the molecular characterization of a large collection of 377 virus isolates recovered from paralyzed patients between 2008 and 2012 in DR Congo and healthy children in 2013 in the Kasai Oriental and Maniema provinces of DR Congo. We found that the genetic diversity of enteroviruses of the species Enterovirus C is more important than previously reported. Interestingly, 50 of the 54 cVDPVs featured recombinant genomes containing exogenous sequences of the 2C ATPase and/or 3D polymerase coding genes acquired from co-circulating Coxsackieviruses A13, A17 and A20. Coxsackieviruses A20 strains were identified for the first time as major partners of genetic recombination with co-circulating live-attenuated polio vaccine strains.

Our findings highlight the need to reinforce and maintain high quality surveillance of PVs and efficient immunization activities in order to ensure early detection and control of emerging cVDPVs in all settings where high frequency and diversity of Coxsackieviruses A13, A17 and A20 have been documented.

Phylogenetic analysis and phenotypic characterisatics of two Tibet EV-C96 strains

BACKGROUND

Enterovirus C96 (EV-C96) is a newly named type of enterovirus belonging to species C, and the prototype strain (BAN00-10488) was firstly isolated in 2000 from a stool specimen of a patient with acute flaccid paralysis in Bangladesh. In this study, we report the genomic and phenotypic characteristics of two EV-C96 strains isolated from individuals from the Tibet Autonomous Region of China.

METHODS

Human rhabdomyosarcoma (RD), human laryngeal epidermoid carcinoma (HEp-2), and human cervical cancer (Hela) cells were infected with the Tibet EV-C96 strains, and enterovirus RNA in the cell culture was detected with a real time RT-PCR-based enterovirus screening method. The temperature sensitivity of Tibet EV-C96 strains were assayed on a monolayer of RD cells in 24-well plates. Full-length genome sequencing was performed by a 'primer-walking' strategy, and the evolutionary history of EV-C96 was studied by maximum likelihood analysis.

RESULTS

Strain 2005-T49 grew in all three kinds of cells, and it was not temperature sensitive. In contrast, none of the three cells produced CPE for strain 2012-94H. Phylogenetic analysis of the two Tibetan viruses, other EV-C96 strains, and EV-C prototypes showed that EV-C96 strains were grouped into three clusters (Cluster1-3) based on their VP1 sequences, which may represent three genotypes. Phylogenetic trees based on the P2 and P3 regions highlighted the difference between Chinese EV-C96 strains and the EV-C96 prototype strain BAN-10488. All Chinese strains formed a cluster separate from BAN-10488, which clustered with CV-A1/CV-A22/CV-A19.

CONCLUSIONS

There is genetic variability between EV-C96 strains which suggest that at least few genetic lineages co-exist and there has been some degree of circulation in different geographical regions for some time. Some recombination events must have occurred during EV-C96 evolution as EV-C96 isolates cluster with different EV-C prototype strains in phylogenetic trees in different genomic regions. However, recombination does not seem to have occurred frequently as EV-C96 isolates from different years and locations appear to cluster together in all genomic regions analysed. These findings expand the understanding of the characterization of EV-C96 and are meaningful for the surveillance of the virus.

High Permissiveness for Genetic Exchanges between Enteroviruses of Species A, including Enterovirus 71, Favors Evolution through Intertypic Recombination in Madagascar

Human enteroviruses of species A (EV-A) are the leading cause of hand-foot-and-mouth disease (HFMD). EV-A71 is frequently implicated in HFMD outbreaks and can also cause severe neurological manifestations. We investigated the molecular epidemiological processes at work and the contribution of genetic recombination to the evolutionary history of EV-A in Madagascar, focusing on the recently described EV-A71 genogroup F in particular. Twenty-three EV-A isolates, collected mostly in 2011 from healthy children living in various districts of Madagascar, were characterized by whole-genome sequencing. Eight different types were identified, highlighting the local circulation and diversity of EV-A. Comparative genome analysis revealed evidence of frequent recent intra- and intertypic genetic exchanges between the noncapsid sequences of Madagascan EV-A isolates. The three EV-A71 isolates had different evolutionary histories in terms of recombination, with one isolate displaying a mosaic genome resulting from recent genetic exchanges with Madagascan coxsackieviruses A7 and possibly A5 and A10 or common ancestors. The engineering and characterization of recombinants generated from progenitors belonging to different EV-A types or EV-A71 genogroups with distantly related nonstructural sequences indicated a high level of permissiveness for intertypic genetic exchange in EV-A. This permissiveness suggests that the primary viral functions associated with the nonstructural sequences have been highly conserved through the diversification and evolution of the EV-A species. No outbreak of disease due to EV-A has yet been reported in Madagascar, but the diversity, circulation, and evolution of these viruses justify surveillance of EV-A circulation and HFMD cases to prevent possible outbreaks due to emerging strains.IMPORTANCE Human enteroviruses of species A (EV-A), including EV-A71, are the leading cause of hand-foot-and-mouth disease (HFMD) and may also cause severe neurological manifestations. We investigated the circulation and molecular evolution of EV-A in Madagascar, focusing particularly on the recently described EV-A71 genogroup F. Eight different types, collected mostly in 2011, were identified, highlighting the local circulation and diversity of EV-A. Comparative genome analysis revealed evidence of frequent genetic exchanges between the different types of isolates. The three EV-A71 isolates had different evolutionary histories in terms of recombination. The engineering and characterization of recombinants involving progenitors belonging to different EV-A types indicated a high degree of permissiveness for genetic exchange in EV-A. No outbreak of disease due to EV-A has yet been reported in Madagascar, but the diversity, circulation, and evolution of these viruses justify the surveillance of EV-A circulation to prevent possible HFMD outbreaks due to emerging strains.

Reference Echovirus 7 and 19 Genomes from Nigeria

We describe the genomes of two echovirus isolates from Nigeria as reference enterovirus species B genomes for the region. These echovirus 7 and 19 genomes have 7,411 nucleotides (nt) and 7,426 nt and were recovered from sewage-contaminated water (in 2010) and an acute flaccid paralysis case (in 2014), respectively.

We describe the genomes of two echovirus isolates from Nigeria as reference enterovirus species B genomes for the region. These echovirus 7 and 19 genomes have 7,411 nucleotides (nt) and 7,426 nt and were recovered from sewage-contaminated water (in 2010) and an acute flaccid paralysis case (in 2014), respectively.

Genetic characterization and molecular epidemiological analysis of novel enterovirus EV-B80 in China

Enterovirus B80 (EV-B80) is a newly identified serotype belonging to the enterovirus B species. To date, only two full-length genomic sequences of EV-B80 are available in GenBank, and few studies on EV-B80 have been conducted in China or worldwide. More information and research on EV-B80 is needed to assess its genetic characteristics, phylogenetic relationships, and association with enteroviral diseases. In this study, we report the phylogenetic characteristics of three Xinjiang EV-B80 strains and one Tibet EV-B80 strain in China. The full-length genomic sequences of four strains show 78.8-79% nucleotide identity and 94-94.2% amino acid identity with the prototype of EV-B80, indicating a tendency for evolution. Based on a maximum likelihood phylogenetic tree based on the entire VP1 region, three genotypes (A-C) were defined, revealing the possible origin of EV-B80 strains in the mainland of China. Recombination analysis revealed intraspecies recombinations in all four EV-B80 strains in nonstructural regions along with two recombination patterns. Due to the geographic factor, the coevolution of EV-B strains formed two different patterns of circulation. An antibody seroprevalence study against EV-B80 in two Xinjiang prefectures also showed that EV-B80 strains were widely prevalent in Xinjiang, China, compared to other studies on EV-B106 and EV-B89. All four EV-B80 strains are not temperature sensitive, showing a higher transmissibility in the population. In summary, this study reports the full-length genomic sequences of EV-B80 and provides valuable information on global EV-B80 molecular epidemiology.

Identification of previously untypable RD cell line isolates and detection of EV-A71 genotype C1 in a child with AFP in Nigeria

We previously attempted to identify 96 nonpolio enteroviruses (EVs) recovered in RD cell culture from children <15 years with acute flaccid paralysis in Nigeria. We succeeded in identifying 69 of the isolates. Here, we describe an attempt to identify the remaining 27 isolates. Twenty-six (the 27th isolate was exhausted) isolates/samples that could not be typed previously were further analyzed. All were subjected to RNA extraction, cDNA synthesis, enterovirus 5'-UTR-VP2 PCR assay and a modified VP1 snPCR assay. Both the 5'-UTR-VP2 and VP1 amplicons were sequenced, isolates identified and subjected to phylogenetic analysis.Twenty of the 26 samples analyzed were identified. Altogether, 23 (three samples had co-infection) EV strains were recovered. These belong to 11 EV (one EVA, nine EVB and one EVC) types which were EVA71 genotype C1 (1 strain), CVB3 (7 strains), CVB5 (1 strain), E5 (2 strain), E11 (3 strains), E13 (2 strain), E19 (1 strain), E20 (1 strain), E24 (2 strains), EVB75 (1 strain) and EVC99 (2 strains). Of the 11 EV types, the 5'-UTR-VP2 assay identified seven while the VP1 assay identified 10. Both assays simultaneously detected 7 of the 11 EV types identified in this study with 100% congruence. We successfully identified 20 of 26 samples that were previously untypable. We also provided evidence that suggests a clade of EVA71 genotype C1 might have been circulating in sub-Saharan Africa since 2008. Finally, we showed that the 5'-UTR -VP2 assay might be as valuable as the VP1 assay in EV identification.

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.