Genetic Characterization of Enterovirus A71 Circulating in Africa

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.

Phylogeographic dynamics and molecular characteristics of Enterovirus 71 in China

Enterovirus 71 (EV71) and coxsackievirus (CV-A16) are the major etiological agents of hand, foot and mouth disease (HFMD). This report reviewed the full-length genomic sequences of EV71 identified in different provinces of China between 1998 and 2019 (a total of 312) in addition to eight worldwide reference genomes to address the genomic evolution and genetic events. The main prevalent EV71 strians in China are C4 genotypes, co-circulating with a few A, B5, C1, and C2 subgenotypes. A new emerging subgenotype in China was identified and classified as B6 genotype. Phylogeographic analysis revealed multiple branches, where a Jiangsu strain 2006-52-9 (GenBank ID: KP266579.1) was linked to different subgenotypes through multiple long mutant branches, including the CV-A16 viruses through the A genotype. Furthermore, identification of 28 natural recombination events suggests that the emergence of new genotypes are associated with intratypic recombination involving EV71 strains and intertypic recombination between EV71 and CV-A16 strains. Compared with the structural proteins, the non-structural proteins of EV71 seem to be highly variable with the highest variable regions of peptidase C3 (3C protein), P2A, and the N-terminus of RNA-dependent RNA polymerase. This study updates the phylogenetic and phylogeographic information of EV71 and provides clues to the emergence of new genotypes of EV71 based on genetics.

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.

Molecular Epidemiology of Enterovirus A71 in Surveillance of Acute Flaccid Paralysis Cases in Senegal, 2013-2020

Enterovirus A71 (EV-A71) is a non-polio enterovirus that currently represents a major public health concern worldwide. In Africa, only sporadic cases have been reported. Acute flaccid paralysis and environmental surveillance programs have been widely used as strategies for documenting the circulation of polio and non-polio enteroviruses. To date, little is known about the molecular epidemiology of enterovirus A71 in Africa where resources and diagnostic capacities are limited. To fill this gap in Senegal, a total of 521 non-polio enterovirus isolates collected from both acute flaccid paralysis (AFP) and environmental surveillance (ES) programs between 2013 and 2020 were screened for enterovirus A71 using real-time RT-PCR. Positive isolates were sequenced, and genomic data were analyzed using phylogeny. An overall rate of 1.72% (9/521) of the analyzed isolates tested positive for enterovirus A71. All positive isolates originated from the acute flaccid paralysis cases, and 44.4% (4/9) of them were isolated in 2016. The nine newly characterized sequences obtained in our study included eight complete polyprotein sequences and one partial sequence of the VP1 gene, all belonging to the C genogroup. Seven out of the eight complete polyprotein sequences belonged to the C2 subgenotype, while one of them grouped with previous sequences from the C1 subgenotype. The partial VP1 sequence belonged to the C1 subgenotype. Our data provide not only new insights into the recent molecular epidemiology of enterovirus A71 in Senegal but also point to the crucial need to set up specific surveillance programs targeting non-polio enteroviruses at country or regional levels in Africa for rapid identification emerging or re-emerging enteroviruses and better characterization of public health concerns causing acute flaccid paralysis in children such as enterovirus A71. To estimate the real distribution of EV-A71 in Africa, more sero-epidemiological studies should be promoted, particularly in countries where the virus has already been reported.

Insight into the Enterovirus A71: A review

Enterovirus A71 is a major causative pathogen of hand, foot and mouth disease. It has become a global public health threat, and is especially important for infants and young children in the Asian-Pacific countries. The enterovirus A71 is a non-enveloped virus of the Picornaviridae family having a single-stranded positive-sense RNA genome of about 7.4 kb which encodes the structural and nonstructural proteins. Currently there are no US FDA-approved vaccines or antiviral therapy available against enterovirus A71 infection. Although enterovirus A71 vaccines have been licenced in China, clinically approved vaccines for widespread vaccination programs are lacking. Substantial progress has recently been achieved on understanding the structure and function of enterovirus A71 proteins together with information on the viral genetic diversity and geographic distribution. The present review is intended to provide an overview on our current understanding of the molecular biology and epidemiology of enterovirus A71 which will aid the development of vaccines, therapeutics and other control strategies so as to bolster the preparedness for future enterovirus A71 outbreaks.

A novel subgenotype C6 Enterovirus A71 originating from the recombination between subgenotypes C4 and C2 strains in mainland China

Recombination plays important roles in the genetic diversity and evolution of Enterovirus A71 (EV-A71). The phylogenetics of EV-A71 in mainland China found that one strain DL71 formed a new subgenotype C6 with unknown origin. This study investigated the detailed genetic characteristics of the new variant. DL71 formed a distinct cluster within genotype C based on the genome and individual genes (5′UTR, VP4, VP1, 2A, 2B, 2C, 3D, and 3′UTR). The average genetic distances of the genome and individual genes (VP3, 2A, 2B, 2C, 3A, 3C, and 3D) between DL71 and reference strains were greater than 0.1. Nine recombination events involving smaller fragments along DL71 genome were detected. The strains Fuyang-0805a (C4) and Tainan/5746/98 (C2) were identified as the parental strains of DL71. In the non-recombination regions, DL71 had higher identities with Fuyang-0805a than Tainan/5746/98, and located in the cluster with C4 strains. However, in the recombination regions, DL71 had higher identities with Tainan/5746/98 than Fuyang-0805a, and located in the cluster with C2 strains. Thus, DL71 was a novel multiple inter-subgenotype recombinant derived from the dominant subgenotype C4 and the sporadic subgenotype C2 strains. Monitoring the emergence of new variants by the whole-genome sequencing remains essential for preventing disease outbreaks and developing new vaccines.

Molecular characterization of non-polio enteroviruses isolated from acute flaccid paralysis patients in Uganda

Enteroviruses (EVs) are RNA viruses that can cause many clinical syndromes including acute flaccid paralysis (AFP). Within the global polio laboratory network, EVs are categorized either as polioviruses or non-polio enteroviruses (NPEVs). Specific NPEVs have been described in polio-like residual paralytic events in AFP patients. Retrospective analysis of 112 NPEV isolates from AFP patients was performed and thirty one NPEV types were identified of which 91% were Enterovirus B and 9% were Enterovirus A species. The NPEVs were distributed across the country with most patients in the eastern region (41/89; 46.1%). The highest proportion of patients were children less than 5 years (77/89; 86.5%) and male patients were more common (54/89; 60.7%). Echovirus 11 (11/89; 12.4%) was frequently observed and phylogenetic analysis of these sequences revealed high diversity. Coxsackievirus B5 (CV-B5), CV-B6, E21, and EV-B69 were only seen in patients with residual paralysis. Analyses of the EV-A71 sequence indicated a unique genogroup.

High Diversity of Human Non-Polio Enterovirus Serotypes Identified in Contaminated Water in Nigeria

Human enteroviruses (EVs) are highly prevalent in sewage and have been associated with human diseases with complications leading to severe neurological syndromes. We have used a recently developed molecular method to investigate the presence of EVs in eight samples collected in 2017–2018 from water streams contaminated by drainage channels in three different locations in Nigeria. A total of 93 human EV strains belonging to 45 different serotypes were identified, far exceeding the number of strains and serotypes found in similar samples in previous studies. Next generation sequencing analysis retrieved whole-capsid genomic nucleotide sequences of EV strains belonging to all four A, B, C, and D species. Our results further demonstrate the value of environmental surveillance for the detection of EV transmission of both serotypes commonly associated with clinical syndromes, such as EV-A71, and those that appear to circulate silently but could eventually cause outbreaks and disease. Several uncommon serotypes, rarely reported elsewhere, were detected such as EV-A119, EV-B87, EV-C116, and EV-D111. Ten EV serotypes were detected in Nigeria for the first time and two of them, CV-A12 and EV-B86, firstly described in Africa. This method can be expanded to generate whole-genome EV sequences as we show here for one EV-D111 strain. Our data revealed phylogenetic relationships of Nigerian sewage strains with EV strains reported elsewhere, mostly from African origin, and provided new insights into the whole-genome structure of emerging serotype EV-D111 and recombination events among EV-D serotypes.

Next Generation Sequencing and Bioinformatics Methodologies for Infectious Disease Research and Public Health: Approaches, Applications, and Considerations for Development of Laboratory Capacity

Next generation sequencing (NGS) combined with bioinformatics has successfully been used in a vast array of analyses for infectious disease research of public health relevance. For instance, NGS and bioinformatics approaches have been used to identify outbreak origins, track transmissions, investigate epidemic dynamics, determine etiological agents of a disease, and discover novel human pathogens. However, implementation of high-quality NGS and bioinformatics in research and public health laboratories can be challenging. These challenges mainly include the choice of the sequencing platform and the sequencing approach, the choice of bioinformatics methodologies, access to the appropriate computation and information technology infrastructure, and recruiting and retaining personnel with the specialized skills and experience in this field. In this review, we summarize the most common NGS and bioinformatics workflows in the context of infectious disease genomic surveillance and pathogen discovery, and highlight the main challenges and considerations for setting up an NGS and bioinformatics-focused infectious disease research public health laboratory. We describe the most commonly used sequencing platforms and review their strengths and weaknesses. We review sequencing approaches that have been used for various pathogens and study questions, as well as the most common difficulties associated with these approaches that should be considered when implementing in a public health or research setting. In addition, we provide a review of some common bioinformatics tools and procedures used for pathogen discovery and genome assembly, along with the most common challenges and solutions. Finally, we summarize the bioinformatics of advanced viral, bacterial, and parasite pathogen characterization, including types of study questions that can be answered when utilizing NGS and bioinformatics.

Development of a New Internally Controlled One-Step Real-Time RT-PCR for the Molecular Detection of Enterovirus A71 in Africa and Madagascar

Enterovirus A71 (EV-A71) is a leading cause of hand-foot-and-mouth disease (HFMD) and can be associated with severe neurological complications. EV-A71 strains can be classified into seven genogroups, A-H, on the basis of the VP1 capsid protein gene sequence. Genogroup A includes the prototype strain; genogroups B and C are responsible of major outbreaks worldwide, but little is known about the others, particularly genogroups E and F, which have been recently identified in Africa and Madagascar, respectively. The circulation of EV-A71 in the African region is poorly known and probably underestimated. A rapid and specific assay for detecting all genogroups of EV-A71 is required. In this study, we developed a real-time RT-PCR assay with a competitive internal control (IC). The primers and TaqMan probe specifically target the genomic region encoding the VP1 capsid protein. Diverse EV-A71 RNAs were successfully amplified from the genogroups A, B, C, D, E, and F, with similar sensitivity and robust reproducibility. Neither cross reaction with other EVs nor major interference with the competitive IC was detected. Experimentally spiked stool and plasma specimens provided consistent and reproducible results, and validated the usefulness of the IC for demonstrating the presence of PCR inhibitors in samples. The analysis in an African laboratories network of 1889 untyped enterovirus isolates detected 15 EV-A71 of different genogroups. This specific real-time RT-PCR assay provides a robust and sensitive method for the detection of EV-A71 in biological specimens and for the epidemiological monitoring of EV-A71 including its recently discovered genogroups.

Identification and molecular characterization of the first complete genome sequence of Human Parechovirus type 15

Using a metagenomics approach, we have determined the first full-length genome sequence of a human parechovirus type 15 (HPeV15) strain, isolated from a child with acute flaccid paralysis and co-infected with EV-A71. HPeV15 is a rarely reported type. To date, no full-length genome sequence of HPeV15 is available in the GenBank database, where only limited VP1 sequences of this virus are available. Pairwise comparisons of the complete VP1 nucleotide and deduced amino acid sequences revealed that the study strain belongs to type 15 as it displayed 79.6% nucleotide and 93.4% amino acid identity with the HPeV15 prototype strain. Comparative analysis of available genomic regions and phylogenetic analysis using the P2 and P3 coding regions revealed low nucleotide identity to HPeV reference genomes. Phylogenetic and similarity plot analyses showed that genomic recombination events might have occurred in the UTRs and nonstructural region during HPeV15 evolution. The study strain has high similarity features with different variants of HPeV3 suggesting intertypic recombination. Our data contributes to the scarce data available on HPeVs in Africa and provides valuable information for future studies that aim to understand the evolutionary history, molecular epidemiology or biological and pathogenic properties of HPeV15.

Genomic and serologic characterization of enterovirus A71 brainstem encephalitis

Objective

In 2016, Catalonia experienced a pediatric brainstem encephalitis outbreak caused by enterovirus A71 (EV-A71). Conventional testing identified EV in the periphery but rarely in CSF. Metagenomic next-generation sequencing (mNGS) and CSF pan-viral serology (VirScan) were deployed to enhance viral detection and characterization.

Methods

RNA was extracted from the CSF (n = 20), plasma (n = 9), stool (n = 15), and nasopharyngeal samples (n = 16) from 10 children with brainstem encephalitis and 10 children with meningitis or encephalitis. Pathogens were identified using mNGS. Available CSF from cases (n = 12) and pediatric other neurologic disease controls (n = 54) were analyzed with VirScan with a subset (n = 9 and n = 50) validated by ELISA.

Results

mNGS detected EV in all samples positive by quantitative reverse transcription polymerase chain reaction (qRT-PCR) (n = 25). In qRT-PCR-negative samples (n = 35), mNGS found virus in 23% (n = 8, 3 CSF samples). Overall, mNGS enhanced EV detection from 42% (25/60) to 57% (33/60) (p-value = 0.013). VirScan and ELISA increased detection to 92% (11/12) compared with 46% (4/12) for CSF mNGS and qRT-PCR (p-value = 0.023). Phylogenetic analysis confirmed the EV-A71 strain clustered with a neurovirulent German EV-A71. A single amino acid substitution (S241P) in the EVA71 VP1 protein was exclusive to the CNS in one subject.

Conclusion

mNGS with VirScan significantly increased the CNS detection of EVs relative to qRT-PCR, and the latter generated an antigenic profile of the acute EV-A71 immune response. Genomic analysis confirmed the close relation of the outbreak EV-A71 and neuroinvasive German EV-A71. A S241P substitution in VP1 was found exclusively in the CSF.

Using Statistical Phylogenetics for Investigation of Enterovirus 71 Genotype A Reintroduction into Circulation

Neurovirulent enterovirus 71 (EV-A71) caused a massive epidemic in China in 2008-2011. While subgenotype C4 was the major causative agent, a few isolates were almost identical to the prototype EV-A71 strain and belonged to genotype A. This variant was allegedly extinct since 1970, and its identification in this epidemic suggests reintroduction of the archive virus. Regression analysis of genetic distances (TempEst software) was of moderate utility due to the low resolution of classical phylogenetic methods. Bayesian phylogenetic analysis (BEAST software) suggested artificial introduction event based on highly aberrant phylogenetic tree branch rates that differed by over three standard deviations from the mean substitution rate for EV71. Manual nucleotide-level analysis was used to further explore the virus spread pattern after introduction into circulation. Upon reintroduction, the virus accumulated up to seven substitutions in VP1, most of them non-synonymous and located within the capsid's canyon or at its rims, compatible with readaptation of a lab strain to natural circulation.

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.

A decade of sustained selection pressure on two surface sites of the VP1 protein of Enterovirus A71 suggests that immune evasion may be an indirect driver for virulence

Enterovirus A71 (EV-A71) is an emerging pathogen in the Enterovirus A species group. EV-A71 causes hand, foot and mouth disease (HFMD), with virulent variants exhibiting polio-like acute flaccid paralysis and other central nervous system manifestations. We analysed all enterovirus A71 complete genomes with collection dates from 2008 to mid-2018. All sub-genotypes exhibit a strong molecular clock with omega (dN/dS) suggesting strong purifying selection. In sub-genotypes B5 and C4, positive selection can be detected at two surface sites on the VP1 protein, also detected in positive selection studies performed prior to 2008. Toggling of a limited repertoire of amino acids at these positively selected residues over the last decade suggests that EV-A71 may be undergoing a sustained frequency-dependent selection process for immune evasion, raising issues for vaccine development. These same sites have also been previously implicated in virus-host binding and strain-associated severity of HFMD, suggesting that immune evasion may be an indirect driver for virulence (154 words).

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.