Next Article in Journal
The Transcriptional Program of Staphylococcus aureus Phage K Is Affected by a Host rpoC Mutation That Confers Phage K Resistance
Previous Article in Journal
How Do Time Since Diagnosis and Sociodemographic Factors Influence Attitudes Towards HIV Status Disclosure in People Living with HIV in Poland? Data from Go Holistic Go Beyond Project
Previous Article in Special Issue
Phylodynamic and Epistatic Analysis of Coxsackievirus A24 and Its Variant
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 16
Issue 11
10.3390/v16111772
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Recent Molecular Epidemiology of Echovirus 11 Throughout North and West Africa Resulted in the First Identification of a Recombinant Strain from an Acute Flaccid Paralysis Case in West Africa

by
Ndack Ndiaye
1,*,†,
Fatou Diène Thiaw
1,
Adamou Lagare
2,
Thérèse Sinare
3,
Mohamed Lemine Diakité
4,
Serigne Fallou Mbacké Ngom
1,
Ousmane Kébé
1,
Issifi Kollo Abdoulkader
2,
Gassim Cissé
5,
Mohamed Dia
1,
Hermann Nodji Djimadoum
6,
Christelle Ouedraogo Neya
3,
Rakia Boubakar
2,
Issaka Ouedraogo
3,
Landoh Dadja Essoya
7,
Ndongo Dia
1,
Amadou Alpha Sall
1,
Ousmane Faye
1 and
Martin Faye
1,*,†
1
Virology Department, Institut Pasteur de Dakar, 36 Avenue Pasteur, Dakar 220, Senegal
2
Centre de Recherche Médicale et Sanitaire (CERMES), 634 Bd de la Nation, Niamey YN034, Niger
3
Ministry of Health and Public Hygiene of Burkina Faso, Ouagadougou 7009, Kadiogo, Burkina Faso
4
Ministry of Health of Mauritania, Avenue Gamel Abdel, Nouakchott 115, Mauritania
5
Ministry of Health and Public Hygiene of Guinea, G77P+56P Boulevard de Commerce, Conakry 585, Guinea
6
World Health Organisation Country Office in Mauritania, ILOT K 140-141 Tevragh-Zeina, Route de la Corniche Ouest, Nouakchott 320, Mauritania
7
World Health Organisation Country Office in Guinea, G8Q8+JC6, Corniche N, Conakry 817, Guinea
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Viruses 2024, 16(11), 1772; https://doi.org/10.3390/v16111772
Submission received: 18 September 2024 / Revised: 8 November 2024 / Accepted: 9 November 2024 / Published: 13 November 2024
(This article belongs to the Special Issue An Update on Enterovirus Research)
Browse Figures
Review Reports Versions Notes

Abstract

:
Echovirus 11 has emerged as a major public health concern, causing sepsis in neonates in many European countries in recent years. In Africa, especially West Africa, where resources and diagnostic capacities are limited, only sporadic cases have been reported. To better understand the recent molecular epidemiology of E11 in West Africa, we characterized twenty-three echovirus 11 strains isolated through the acute flaccid paralysis and environmental surveillance systems for polio from 2013 to 2023, using high-throughput sequencing. Our data are noteworthy due to identifying for the first time a recombinant strain from an acute flaccid paralysis case and represent the first focus to date on molecular characterization of echovirus 11 in West Africa. Moreover, our data show that echovirus 11 diverged from 1970 (95% HPD range, 1961–1979) and evolved into four distinct clades, with the virus spread from West Africa to Europe, exhibiting two introductions in France around 2017, from Senegal and Guinea. Furthermore, the in silico analysis reveals four non-conservative amino acid substitutions in the VP1 sequences of the European strains associated with neonatal sepsis in newborns and a conserved amino acid motif in the VP1 protein toward enterovirus genotypes. Our data provide new insights into the epidemiology of echovirus 11 and point to the crucial need to implement specific surveillance programs targeting non-polio enteroviruses for the rapid identification of emerging or re-emerging enterovirus species, particularly in Africa.

1. Introduction

Enteroviruses (EVs) are the most common circulating viruses worldwide and their transmission generally occurs through the fecal–oral route but may also occur via respiratory droplets [1]. Other modes of transmission, such as respiratory and materno-fetal routes, have also been reported [2]. They are associated with a wide spectrum of illnesses, ranging from mild non-specific symptoms, febrile illness, and hand, foot and mouth disease (HFMD) to severe neurological CNS infections such as meningitis, encephalitis and acute flaccid paralysis (AFP) [2]. EVs are a large group of positive-sense single-stranded RNA viruses belonging to the Enterovirus genus from the family Picornaviridae [3]. Small in size, at 15–30 nm, the EV genome has approximately 7500 nucleotides (nt) and is composed of a large open reading frame (ORF) flanked by 5′ and 3′ untranslated regions (UTRs) [3,4]. The ORF can be translated into a 2189-amino-acid-long polyprotein and then cleaved into three polyprotein precursors, P1, P2, and P3. P1 encodes the structural proteins VP4, VP2, VP3, and VP1, while P2 and P3 encode the non-structural proteins 2A, 2B, and 2C, and 3A, 3B, 3C, and 3D, respectively [4]. Sequencing of the VP1 gene is mainly used for EV genotyping as it provides a perfect correlation with the genotype of the strains and an excellent discrimination of all the taxonomic ranks (from species to variants) [5].
Based on the VP1 protein, which is the most diversified region, EVs include fifteen species, of which viruses in four EV species (EV-A to D) and three rhinovirus species (RV-A to C) infect humans [6]. Certain viruses, such as polioviruses, Coxsackievirus A (CVA), Coxsackievirus B (CVB) and echoviruses, have been classified or assigned to EV species [7].
EVs of the B species (EV-B), particularly echoviruses, represent the most common species associated with infections of the central nervous system (CNS) in Asia, Europe, and the USA during the last decade [8,9,10]. Previously reported outbreaks of aseptic meningitis (AM) mainly involved the CVB genotypes and the echovirus types E4, 6, 9, 11, 13, 18 and 30 [11].
Recently, the World Health Organization (WHO) reported an increasing number of severe neonatal infections associated with E11 in many European countries, notably in France, Italy, Croatia, Spain, Sweden, and the United Kingdom [12,13]. These infections were associated with severe neonatal sepsis, hepatic impairment and multi-organ failure in France [14] and severe hepatitis in Italy [15]. E11 is one of the echovirus types most frequently responsible for neonatal infections [16,17,18], and outbreaks due to E11 in newborns are generally associated with high morbidity and mortality [19,20].
To date, there is a scarcity of data on the epidemiology of E11 conducted in Africa, where knowledge and resources are limited. However, previous data from Senegal provided evidence of their circulation in children with AFP [21]. To fill this gap, we investigated the recent molecular epidemiology of E11 in West Africa using genomic sequencing and phylodynamic analyses to better understand the evolutionary mechanisms that have driven E11’s genetic evolution.

2. Materials and Methods

2.1. Ethical Statement

This study is a component of the Global Polio Eradication Initiative. This study did not involve human participants but rather the use of cell culture isolates of viruses recovered from stool specimens of AFP cases collected through routine poliomyelitis surveillance activities at the instigation of the Word Health Organization (WHO) for public health purposes. All the technical and ethical aspects were approved by the WHO and the Ministries of Health of the countries concerned. The protocol and oral consent were determined as routine surveillance activity by the steering committee of the WHO in compliance with all the applicable national regulations governing the protection of human subjects.

2.2. Data Collection

In this study, a dataset including E11’s VP1, capsid and complete genome sequences retrieved from GenBank (https://www.ncbi.nlm.nih.gov/genbank/) (accessed on 21 February 2024) and the Viral Bioinformatics Resource Centre (BV-BRC) (https://www.bv-brc.org/view/Taxonomy/12058#view_tab=genomes) (accessed on 21 February 2024) [21,22], and the newly characterized E11 sequences obtained from NPEV-positive isolates identified from 2013 to 2023 at the Inter-country Reference Laboratory for Polio in Senegal according to the WHO guidelines [23,24], was analyzed.

2.3. RNA Extraction and Molecular Testing

The viral RNA was extracted from 200 µL of virus isolates using the QIAmp Viral RNA Mini Kit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. The extracted RNA was eluted in a final volume of 60 µL of nuclease-free water. Subsequently, the RNA extracts were confirmed by specific real-time reverse transcriptase polymerase chain reaction (RT-qPCR) using a pan-enterovirus (panEV) RT-PCR assay [25] with the qScript™ XLT One-Step RT-PCR (Quanta Bio, Beverly, MA, USA) according to the manufacturer’s instructions. In addition, differential singleplex RT qPCR assays for the specific detection of echovirus [26] were performed using specific primers and probes with RNA extracts. The experiments were performed using the CFX96 Real-Time PCR system (Bio-Rad, Singapore) and the echovirus-positive RNAs were also stored at -80 °C until further testing.

2.4. Complete Genome Sequencing

Echovirus-positive RNAs with threshold cycle (Ct) values lower than 25 were selected and depleted using mammalian ribosomal RNA-specific primers and RNAseH enzyme (NEB). The purified viral RNA was reverse-transcribed into first-stranded cDNA using random hexamer primers with the Invitrogen SuperScript IV Reverse Transcriptase kit (Thermo Fisher, Waltham, MA, USA). Then, after the second strand cDNA synthesis, the fragmented amplicons were tagged using the “indexes” incorporated into the Nextera XT DNA Library Preparation kit (Illumina, San Diego, CA, USA). The purified libraries were quantified using a Qubit 3.0 fluorometer (Invitrogen Inc., Waltham, MA, USA) and normalized to 200 pM before pooling in a single tube. The pooled libraries were loaded on an Illumina iSeq100 instrument using the iSeq100 reagent v2 for 300 cycles.

2.5. Sequencing of the Entire Coding Region of the Capsid Protein

For isolates that provided partial sequences or failed using the Illumina sequencing method, amplicons of the entire coding region of the capsid protein were generated using the OneTaq One-Step RT-PCR kit (New England Biolabs, Ipswich, MA, United States of America) with a specific PCR method, as previously described [27]. Briefly, the PCR was carried out in a final volume of 50 μL, containing 10 μL of the RNA product, 25 μL of the 2× OneTaq One-Step Reaction Mix, 9 μL of nuclease-free water, 2 μL of forward primer (10 μM), 2 μL of reverse primer (10 μM) and 1 μL OneTaq One-Step Enzyme Mix. This reaction’s master mix was amplified using the following conditions: 15 min at 48 °C, 1 min at 94 °C, followed by 40 cycles of 15 s at 94 °C, 30 s at 55 °C and 4 min at 68 °C and a final extension step at 68 °C for 5 min. The amplicon size of the PCR products was verified using electrophoresis on a 1% agarose gel containing SafeViewTM Classic (Applied Biological Materials, 3671 Viking Way Unit 1, Richmond, BC, Canada) and visualized under UV fluorescence using the Syngene InGenius3 device (Fisher Scientific, Hampton, New Hampshire, États-Unis).
Amplicons of the capsid protein were purified using the AMPure XP magnetic beads and quantified using the dsDNA High Sensitivity Kit on a Qubit 3.0 fluorometer (Thermo Fisher). The purified entire coding regions of the capsid protein were barcoded using the Rapid Barcoding Kit (SQK.110.96) (Oxford Nanopore Technologies, Gosling Building, Edmund Halley Road, Oxford Science Park, UK) with the MRT001 expansion module (Oxford Nanopore Technology) and pooled in a single tube. The libraries were then purified and sequenced on a GridION instrument (Oxford Nanopore Technology, Gosling Building, Edmund Halley Road, Oxford Science Park, UK).

2.6. Sequencing Data Analysis

The raw data from both the Illumina and nanopore experiments were processed, merged into a single FastQ file using a bash command-line pipeline, and consensus sequences were generated using the “Genome detective virus tool” software (version 2.40) [28]. The newly characterized sequences were analyzed using the online Basic Local Alignment Search (BLAST) program (https://blast.ncbi.nlm.nih.gov/Blast.cgi (accessed on 29 February 2024) to assess their homology with previously available data. In addition, the genotype was confirmed using the online Enterovirus Genotyping Tool (RIVM), version 0.1. program (enterovirus/typingtool/).

2.7. Recombination Analysis

To assess the occurrence of recombination events, similarity plotting and bootscan analysis were conducted on our dataset by using the SimPlot program (version 3.5.1) [29] with a 400 nt window moving in 20 nt steps and the Kimura two-parameter method with a transition–transversion ratio of 2 with 1000 resampling. Two related reference sequences of parental genotypes previously available via GenBank and one known outgroup (irrelevant) sequence were included in our dataset and simultaneously analyzed to assess the possible breakpoints of intergenotypic recombination. These reference sequences have been chosen as they were considered to be both representative and phylogenetically informative.
In addition, the RDP5 (Recombination Detection Program version 5) program [30] was also used to detect the potential recombinants confirmed by at least 6 of the 7 selected methods using the default settings, including RDP, GENECONV, Maxchi, Bootscan, Siscan, Chimaera, 3Seq and LARD. Full genome sequences were used as queries for BLASTn and sequences with the highest homologies and complete genomes were used in the recombination analysis.

2.8. Assessment of Selection Pressures

The ratio of nonsynonymous to synonymous amino acid substitutions (dN/dS) is a useful indicator of the strength of natural selection acting on protein-coding genes. This is a widely used method to detect positive selection. The dN/dS statistical test allowed us to distinguish diversifying or positive selection (dN/dS > 1) from negative or purifying selection (dN/dS < 1). For this, a total of 2 alignment partitions were performed corresponding to the entire coding region of the capsid and VP1 protein. As the site model, we used single-likelihood ancestor counting (SLAC), which estimated the difference between the nonsynonymous (dN) and synonymous (dS) rates per codon site at the 0.1 significance level. The fast, unconstrained Bayesian approximation (FUBAR) method, which evaluated the episodic positive selection at each site in the alignment at a posterior probability ≥0.9, was also used [31]. Finally, the mixed effects model of evolution (MEME) was also conducted at a 0.1 significance level for estimating the selective pressure changes among codon sites. All three methods were conducted with the HyPhy package implemented in the Datamonkey web server (http://www.datamonkey.org, accessed 14 March 2024) [32,33]. An episode of positive diversifying selection was considered if it was detected by at least two different methods.

2.9. Phylogenetic Analyses

Multiple alignments were performed using our dataset, including the newly characterized sequences and previous sequences of the capsid protein, the VP1 protein, complete E11 genomes and one outgroup retrieved from GenBank (https://www.ncbi.nlm.nih.gov/genbank, accessed on 26 February 2024), with the ClustalW method implemented in the BioEdit software (version 7.2.5) [34]. The maximum likelihood (ML) phylogenetic tree was inferred using IQ-TREE (version 1.6.12), with the GTR+F+G4 model as the best-fitted nucleotide (nt) substitution model for our dataset, and the branch support was calculated using ultra-fast bootstrapping with 1000 replications [35]. The ML tree topology was visualized with FigTree (version 1.4.4) [36] and the nodes were supported by bootstrap values.
The root-to-tip divergence as well as the regression slopes and correlations analyses were performed using the TempEst v1.5.3 [37] software to assess the reliable temporal signal of the E11 VP1 sequences. Evolutionary analysis of E11 was performed using BEAST2 v2.7.6 [38] to estimate the mean substitution rate and the time of the most recent common ancestor (tMRCA). A maximum clade credibility (MCC) tree was constructed by incorporating a 10% burn-in period. The median heights were generated from a 50 million generation with a Bayesian Markov chain Monte Carlo (MCMC) approach for phylogenetic reconstruction sampling every 1000 generations. In addition, a general time-reversible nucleotide substitution model (GTR) with a gamma distribution prior to each relative substitution rate and a strict clock model were used to infer the echovirus E-11 virus evolutionary timescale. The analysis assumed a strict clock model and a coalescent Bayesian Skyline as the tree priors. Convergence was confirmed by achieving an effective sample size (ESS) greater than 200, using Tracer v1.7.2. [39].

2.10. In Silico Analysis Towards the Vaccine Strains

The amino acid polymorphisms on the E11’s VP1 were assessed using a total of 33 sequences, including African and European strains detected in AFP and sepsis, respectively. In addition, the conserved structures were assessed using the web-based ConSurf server (http://consurf.tau.ac.il/ (accessed on 18 May 2024) [40] with the 21_279_SEN_2021 AFP Senegalese E11 strain protein sequence as a query and the multiple sequence alignment as the input. The multiple sequence alignment (MSA) was built using MAFFT. The degree of conservation was calculated using a Bayesian method with the most-fitted evolutionary substitution model for our dataset, considering nine distinct grades, with grade 1 being the least and grade 9 being the most conserved.

3. Results

3.1. Characteristics of West African E11 Strains Analyzed in This Study

Twenty-three E11 sequences from West African isolates were obtained, including eight VP3-VP1 proteins, two VP1 proteins, twelve entire coding regions of the capsid protein and one complete genome sequence. Moreover, 4 out of these 23 sequences were identified from the environmental surveillance (ES), especially in 2023 (Table 1). A rate of 52.2% (n =12) of sequences originated from Senegal. The median age of these E11-positive AFP cases was 2 years (ranging from 16 months to 5 years) and the sex ratio was 3.6 (with 21.7% missing data (n = 5)). The virus has been isolated in West Africa since 2013, with an increase since 2021. Moreover, it is mainly identified between March and November, which corresponds to the summer and rainy season in Senegal, with a higher rate in September (31.5% (n = 6)). In addition, all the E11 from AFP were identified from single infections, while a mixture was identified for ES involving E19 and E12.

3.2. Recombination and Selection Pressures Assessment

One complete E11 genome sequence from Senegal (21-279-SEN-2021) was generated, with a length of 7357 nucleotides (nt), that encoded a polyprotein of 2452 amino acids. The complete genome was closely related (99% nt similarity) to an E11 strain previously identified in non-human primates (NHP) [41] from Nigeria in 2016 (accession number MN657230.1).
The SimPlot graph revealed a difference in the nt similarity in the non-structural genomic region between the 21-279-SEN-2021 strain and a closely related sequence, suggesting a possible recombination event. Furthermore, the occurrence of a breakpoint at nt position 4062 in the P2 region was detected using the bootscanning method. In addition, this recombination event was confirmed using RDP5 analysis, which showed that the E11 strain 21-279-SEN-2021 sequence from Senegal recombined with the E11 strain isolated in NHP from Nigeria in 2016 (accession number MN657230.1) (Figure S1). However, further study is needed since the Nigerian minor parent strain has also been described as recombinant.
Phylogenetic analyses of the complete genome sequence (A), P1 (B), P2 (C) and P3 (D) coding regions of the 21-279-SEN-2021 strain confirmed the RDP5 findings, with the 21-279-SEN-2021 strain clustering with the NHP E11 sequence from Nigeria (accession number MN657230.1) (Figure 1).
The entire coding regions of the capsid and the VP1 protein were analyzed separately for the estimation of sites under positive diversifying selection, applying four methods to ensure the consistency of these events along the structural gene sequences. Using this approach, we found several sites under strong negative selection in both the entire coding region of the capsid and VP1 protein using the SLAC and FUBAR methods. However, sites under episodic positive selection were only observed for the entire coding region of the capsid protein with the MEME method (p < 0.1), and therefore, this event was not considered to be significant (Table 2).

3.3. Phylogenetic and Evolutionary Analysis

All the newly characterized E11 nucleotide sequences aligned, revealing a high identity at the nucleotide level. Moreover, the phylogenetic ML tree inferences based on the entire coding region of the capsid protein (Figure 2A) and VP1 (Figure 2B) and partial VP1 (Figure 2C) sequences confirmed the data from the BLAST and RIVM and revealed that the E11 sequences were closely related to previous strains, mainly from Africa, and had nucleotide homology ranging from 80 to 100%. Interestingly, our data showed a monophyletic sub-regional clustering of E11 genotype in Africa, which is distinct from the strains that were recently associated with sepsis in Europe (Figure 2).
The root-to-tip analysis revealed a positive temporal signal for E11’s VP1 protein sequences (Figure 3A). The maximum clade credibility (MCC) tree confirmed the genotyping results (Figure 3B) and revealed that the evolutionary substitution rate for E11’s VP1 protein sequence was 7.637 × 10−3 substitutions/site/year. Moreover, this analysis showed that the newly characterized E11 sequences were distinct from recent strains identified from the European countries and were not native to their isolated country. In fact, the routes of virus introductions into the West African continent vary. However, North African (Libya, Algeria) and European (France) countries seem to have been the key geographical sources for these newly characterized E11 sequences. Furthermore, we have not observed a temporal and spatial correlation between clusters, suggesting an active virus circulation between countries.
In fact, our data showed that the estimated time to the most recent common ancestor (tMRCA) of the E11 strains diverged from 1970 (95% HPD range, 1961–1979). Subsequently, the virus evolved continuously in four distinct clades: the overall newly characterized E11 sequence grouped into clades I–III and the recent European strains belonging to clade IV. Our data exhibited virus dynamics not only throughout North and West Africa but also its spread from West Africa to Europe. They showed two E11 introductions in France in around 2017, from Senegal and Guinea.
Subsequently, clade I, including the newly characterized Senegalese sequences, has experienced two introductions to Senegal; one from Ghana through Niger in 2013 (95% HPD, 2012–2014) and another from France through Guinea in 2019 (95% HPD 2018–2020). Clade II, grouping Senegalese and Nigerian strains, originated from Libya in 2009 (95% HPD, 2006–2012) while clade III, including Mauritanian strains, spread from Algeria in 2011 (95% HPD, 2009–2013). In addition, the virus spread from Senegal in 2005 (95% HPD, 2003–2007) and was successively introduced to Mauritania and France. Clade IV seems to be restricted to Europe (Figure 3).

3.4. VP1 Amino Acid Sequence Analysis

In silico analysis of the amino acid sequences of the VP1 protein has shown a polymorphism for European strains recently associated with sepsis in Europe, including four non-conservative substitutions located at amino acid positions 42, 89, 157 and 211 (Figure 4A). These data suggest a probable impact of these amino acid transitions on the viral pathogenesis or biological characteristics. Furthermore, data from the ConSurf prediction exhibited a conserved amino acid motif in the VP1 protein, located at aa positions 162–169 and corresponding to the EF loop (Figure 4B).

4. Discussion

Within the human EV-B species, echoviruses are the most commonly isolated NPEV associated with neurological diseases such as AFP [21,22,42], AM and multisystem hemorrhagic disease [42]. Epidemics of CNS infections involving echoviruses have been previously reported in several countries, with high morbidity and mortality [43]. From 1983 to 2003, case fatality rates ranging from 16.7% to 20% were reported in the United States among patients admitted for CNS infections due to echovirus genotypes such as E6, E7, E11, E20, and E30 [16]. In addition, CNS infections as AM outbreaks mainly involved the E4, E6, E9, E11, E13, E18 and E30 genotypes [6,13].
The recent E11 outbreak in multiple European countries between July 2022 and 2023 was associated with severe acute and fulminant hepatitis in newborns, affecting mainly males [13,14,15]. In Africa, E11 was rarely reported, so only a few data derived from the AFP surveillance and ES are currently available [21,22,41]. Therefore, the surveillance system for AFP and ES could largely contribute to a better understanding of the recent molecular epidemiology of E11 in Africa. Herein, we reported on the genetic diversity and phylodynamics of E11 using high-throughput sequencing, revealing the first complete sequence from an AFP in West Africa to date.
Although echovirus is commonly associated with both epidemic and endemic patterns of infection in individuals of all ages, all the E11-positive isolates from our study included children under 6 years (Table 1), indicating the impact of infections in this vulnerable age group. In addition, recent studies have shown E11’s association with severe disease in neonates, particularly during the first two weeks of life and in those born prematurely [13,14,15]. However, severe complications such as fetal infections in utero have also been associated with disease and/or death [44]. Indeed, human neonatal Fc receptor (FcRn), which is the uncoating receptor for major EV-B [45], is highly expressed on syncytiotrophoblasts, the fetal-derived cells that comprise the outermost cellular barrier of the human placenta and that directly contact maternal blood [46]. Moreover, FcRn can mediate severe E11 infection in suckling mice, but only when they are interferon-deficient, suggesting that interferon signaling may be important for age-related susceptibility [47].
Furthermore, the notion that males are more likely to contract NPEV infection, as previously reported [48], was established in our study (Table 1). Nevertheless, there is a crucial need for further study on sex as a risk factor for enterovirus infections.
Even if NPEVs circulate year-round with a peak in summer [49], we have shown that E11 was mainly identified between March and November (Table 1), suggesting that the prevalence of these viruses could be associated with climate factors such as high temperature and high humidity, as reported in a previous study in China [45].
In our study, four strains were identified from environmental surveillance, demonstrating not only its importance to complement the AFP surveillance but also its sensitivity as an early warning system for pathogens of concern such as NPEV [46].
Interestingly, the occurrence of recombination (Figure 1 and Figure S1) of the newly characterized complete genome sequence with a strain isolated from non-human primates in Nigeria [41] at the non-structural region suggests a potential intratypic recombination during E11’s evolution. In fact, recombination has been reported as a major mechanism participating in the evolution of enteroviruses viruses, especially those belonging to the EV-B species [49,50], and could trigger serious public health problems [51]. In addition, previous data on EV species suggested that intertypic or intratypic recombination frequently occurs in the non-structural regions [52,53], suggesting that these genomic regions represent relatively non-stable units.
The assessment of the selection pressures acting on the E11 genome (Table 2) showed episodes of strong purifying selection in the entire coding region of the capsid and VP1 protein, indicating the stability of these regions. In fact, negative selection pressures, defined as the main evolutionary features avoiding the occurrence of new mutations, were mainly described in the VP1 protein of enteroviruses [54]. Although it could be associated with the pressure from the patient’s immune system, genomic surveillance of EV could be further promoted for rapid identification of possible recombinant strains [54].
The phylogenetic analysis (Figure 2) showed that West African E11 strains formed one group and clustered with previous E11 strains from AFP cases in Ghana in 2017 [55] and from non-human primates in Nigeria in 2016. Similar data have been previously reported for a sample collected in the environment in 2012 [41], suggesting a probable zoo-anthroponotic transmission of enteroviruses from humans to non-human primates. Thus, future studies relying on the assessment of risk factors, such as direct or indirect contact with non-human primates, could be further promoted.
The newly characterized strains showed a high nucleotide identity, which indicated that the West African E11 strains emerged a long time before the strains recently identified in the European Union region. Like other RNA viruses, enteroviruses are characterized by evolutionary mechanisms that shape the genomic diversity by the accumulation of mutations, insertions or deletion due to the lack of proofreading activity of the 3D polymerase [56]. This mechanism is crucial for viral adaptability, dissemination, and pathogenesis [57]. Our data underscore the importance of better monitoring E11 by establishing a sustainable genomic surveillance program for enteroviruses in Africa. To date, there is a low number of complete African genome sequences (n = 502) available in GenBank (https://www.ncbi.nlm.nih.gov/genbank, (accessed on 21 February 2024)), and 96% are poliovirus stains (n = 481). In fact, the monitoring of E11 circulation was mainly limited to the AFP surveillance system and based on only partial VP1 sequences. In addition, an assessment of the pathogenicity of the West African E11 strains in comparison to those recently associated with sepsis in Europe is warranted.
Our Bayesian MCMC analyses (Figure 3) allowed us to compare the mean evolutionary rate of E11 (7.637 × 10−3 substitutions/site/year), which was faster than those of other enterovirus genotypes such as EV-A71 [58]. Furthermore, we identified the virus’s introduction and circulation in Senegal several years ago, suggesting a silent E11 circulation until its detection in AFP cases in 2013. Further, based on our dataset, E11 has been introduced into West Africa from North African (Libya, Algeria) and European (France) countries. In fact, E11 has been frequently identified as the causative agent of outbreaks in these countries [59,60]. Therefore, E11’s introductions may be related to its high international human mobility, as previously reported [61].
Furthermore, as also described for NPEV [62], it was found that amino acid changes on E11’s VP1 region may affect its viral virulence. Considering that the VP1 is notably the receptor-binding protein of enteroviruses and the major neutralizing antigen, it is therefore a major viral factor responsible for cell, tissue and host tropism [63]. Thus, mutations occurring in this region might affect the binding and the uncoating process of E11 and increase its transmission ability.
In addition, one amino acid substitution was located in the critical binding regions for neutralizing antibodies that correspond to the BC loop (located at amino acids 57 to 89), which have been associated with virus escape to the immune system [64]. However, their significance and mechanism require further investigations.
The increasing and continuous identification of emerging and/or re-emerging enteroviruses point to the crucial need for vaccine development. The conserved amino acid motifs in the VP1 protein could not only provide more advantages in the development of a broad-spectrum anti-enteroviral drug but also in vaccine design as it plays a pivotal role in the virus’s life cycle.

5. Conclusions

E11 has emerged as a recent public health concern across Europe and better understanding of its molecular epidemiology in West Africa is crucial for preparedness, monitoring and control. Our study is noteworthy due to reporting for the first time on the molecular epidemiology of E11 circulating in West Africa and the first identification of a recombinant strain isolated from an AFP case in West Africa. In addition, our data provide new insights into E11’s emergence and circulation dynamic, and useful information for disease control and prevention. Besides AFP and ES surveillance programs, more specific hospital-based surveillance targeting E11, including genomic sequencing, could also be promoted to assess the burden of E11 in Africa. More experimental investigations focusing on the pathogenicity of the identified E11 isolates in West Africa are warranted.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/v16111772/s1, Figure S1: Recombination analysis of the whole-genome sequence of the 21-279-SEN-2021 strain using RDP5.

Author Contributions

Conceptualization: N.N. and M.F.; funding acquisition: N.D., A.A.S., O.F. and M.F.; field investigations: A.L., T.S., M.L.D., I.K.A., G.C., H.N.D., C.O.N., R.B., I.O. and L.D.E.; formal analysis: N.N., S.F.M.N., F.D.T., O.K., M.D. and M.F.; supervision: M.F.; writing: N.N. and M.F.; review and editing: N.N., F.D.T., A.L., T.S., M.L.D., S.F.M.N., I.K.A., G.C., H.N.D., C.O.N., R.B., I.O., L.D.E., O.K., M.D., N.D., A.A.S., O.F. and M.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable. As part of the GPEI, our study did not directly involve humans or animals but rather included cell culture isolates from the AFP surveillance and ES collected as part of routine surveillance activities for polio in Senegal for public health purposes. A cleared, oral consent was obtained from all the patients or their parents or relatives.

Informed Consent Statement

The primary stools, sewage samples and NPEV isolates were collected as part of routine surveillance activities for polioviruses. Our protocol was approved by the WHO and the National Ethical Committee at the Ministry of Health and Social Actions in Senegal. A cleared, oral consent was obtained from all the patients or their parents or relatives.

Data Availability Statement

All the data are available in the present manuscript.

Acknowledgments

We would like to thank all the staff of the virology department, and the authors for their support, particularly the team from the Inter-country Reference Center for Poliomyelitis in Senegal and colleagues in partner countries.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Itani, T.; Chalapa, V.; Semenov, A.; Sergeev, A. Laboratory diagnosis of nonpolio enteroviruses: A review of the current literature. Biosaf. Health 2023, 5, 112–119. Available online: https://www.sciencedirect.com/science/article/pii/S259005362200177X (accessed on 21 May 2024). [CrossRef]
  2. Singh, S.; Mane, S.S.; Kasniya, G.; Cartaya, S.; Rahman, M.M.; Maheshwari, A.; Motta, M.; Dudeja, P. Enteroviral Infections in Infants. Newborn 2022, 1, 297–305. Available online: https://pmc.ncbi.nlm.nih.gov/articles/PMC9599990/ (accessed on 21 May 2024). [PubMed]
  3. Zell, R.; Delwart, E.; Gorbalenya, A.E.; Hovi, T.; King, A.M.Q.; Knowles, N.J.; Lindberg, A.M.; Pallansch, M.A.; Palmenberg, A.C.; Reuter, G. ICTV virus taxonomy profile: Picornaviridae. J. Gen. Virol. 2017, 98, 2421. [Google Scholar] [CrossRef] [PubMed]
  4. Knowles, N.J.; Hovi, T.; Hyypiä, T.; King, A.M.Q.; Lindberg, A.M.; Pallansch, M.A.; Palmenberg, A.C.; Simmonds, P.; Skern, T.; Stanway, G.; et al. Picornaviridae. In Virus Taxonomy: Classification and Nomenclature of Viruses: Ninth Report of the International Committee on Taxonomy of Viruses; King, A.M.Q., Adams, M.J., Carstens, E.B., Lefkowitz, E.J., Eds.; Elsevier: San Diego, CA, USA, 2012; pp. 855–880. [Google Scholar]
  5. Oberste, M.S.; Maher, K.; Kilpatrick, D.R.; Pallansch, M.A. Molecular evolution of the human enteroviruses: Correlation of serotype with VP1 sequence and application to picornavirus classification. J. Virol. 1999, 73, 1941–1948. [Google Scholar] [CrossRef] [PubMed]
  6. Pallansch, M.A.; Oberste, M.S.; Whitton, J.L. Enteroviruses: Polioviruses, coxsackieviruses, echoviruses, and newer enteroviruses. In Fields Virology, 6th ed.; Knipe, D.M., Howley, P.M., Cohen, J.I., Griffen, D.G., Lamb, R.A., Martin, M.A., Racaniello, V.R., Roizman, B., Eds.; Lippincott Williams & Wilkins: Philadelphia, PA, USA, 2013; Volume 1, pp. 490–530. [Google Scholar]
  7. Walker, P.J.; Siddell, S.G.; Lefkowitz, E.J.; Mushegian, A.R.; Adriaenssens, E.M.; Dempsey, D.M.; Dutilh, B.E.; Harrach, B.; Harrison, R.L.; Hendrickson, R.C. Changes to virus taxonomy and the Statutes ratified by the International Committee on Taxonomy of Viruses (2020). Virol. Division News 2020, 164, 2737–2748. [Google Scholar] [CrossRef]
  8. Abedi, G.R.; Watson, J.T.; Nix, W.A.; Oberste, M.S.; Gerber, S.I. Enterovirus and parechovirus surveillance—United States, 2014–2016. Morb. Mortal. Wkly. Rep. 2018, 67, 515. [Google Scholar] [CrossRef]
  9. Lafolie, J.; Labbé, A.; L’honneur, A.S.; Madhi, F.; Pereira, B.; Decobert, M.; Adam, M.N.; Gouraud, F.; Faibis, F.; Arditty, F. Assessment of blood enterovirus PCR testing in paediatric populations with fever without source, sepsis-like disease, or suspected meningitis: A prospective, multicentre, observational cohort study. Lancet Infect. Dis. 2018, 18, 1385–1396. [Google Scholar] [CrossRef]
  10. Chen, X.; Li, J.; Guo, J.; Xu, W.; Sun, S.; Xie, Z. An outbreak of echovirus 18 encephalitis/meningitis in children in Hebei Province, China, 2015. Emerg. Microbes Infect. 2017, 6, e54. [Google Scholar] [CrossRef]
  11. Lee, B.E.; Davies, H.D. Aseptic meningitis. Curr. Opin. Infect. Dis. 2007, 20, 272–277. Available online: https://journals.lww.com/co-infectiousdiseases/Fulltext/2007/06000/Aseptic_meningitis.7.aspx (accessed on 12 February 2024). [CrossRef]
  12. Mondiale de la Santé, O.; World Health Organization. WHO’s Health Emergencies Programme: Acute emergencies monthly summary–May 2023–Programme OMS de gestion des situations d’urgence sanitaire: Résumé mensuel des situations d’urgence aiguë–mai 2023. Wkly. Epidemiol. Rec. Relevé Épidémiologique Hebd. 2023, 98, 281–286. [Google Scholar]
  13. Mondiale de la Santé, O.; World Health Organization. Weekly Epidemiological Record, 2023, vol. 98, 26 [full issue]. Wkly. Epidemiol. Rec. Relevé Épidémiologique Hebd. 2023, 98, 279–286. [Google Scholar]
  14. Grapin, M.; Mirand, A.; Pinquier, D.; Basset, A.; Bendavid, M.; Bisseux, M.; Jeannoël, M.; Kireche, B.; Kossorotoff, M.; L’Honneur, A.-S. Severe and fatal neonatal infections linked to a new variant of echovirus 11, France, July 2022 to April 2023. Eurosurveillance 2023, 28, 2300253. [Google Scholar] [CrossRef] [PubMed]
  15. Piralla, A.; Borghesi, A.; Di Comite, A.; Giardina, F.; Ferrari, G.; Zanette, S.; Figar, T.A.; Angelini, M.; Pisoni, C.; Pitrolo, A.M.G. Fulminant echovirus 11 hepatitis in male non-identical twins in northern Italy, April 2023. Eurosurveillance 2023, 28, 2300289. [Google Scholar] [CrossRef] [PubMed]
  16. Khetsuriani, N.; LaMonte, A.; Oberste, M.S.; Pallansch, M. Neonatal enterovirus infections reported to the national enterovirus surveillance system in the United States, 1983–2003. Pediatr. Infect. Dis. J. 2006, 25, 889–893. [Google Scholar] [CrossRef] [PubMed]
  17. Bersani, I.; Auriti, C.; Piersigilli, F.; Dotta, A.; Diomedi-Camassei, F.; Di Pede, A.; Buttinelli, G.; Danhaive, O. Neonatal acute liver failure due to enteroviruses: A 14 years single NICU experience. J. Matern. Fetal Neonatal Med. 2020, 33, 2576–2580. [Google Scholar] [CrossRef]
  18. Davies, D.P.; Hughes, C.A.; MacVicar, J.; Hawkes, P.; Mair, H.J. Echovirus-11 infection in a special-care baby unit. Echovirus-11 Infect. Spec.-Care Baby Unit 1979, 1, 96. [Google Scholar] [CrossRef]
  19. Mostoufizadeh, M.; Lack, E.E.; Gang, D.L.; Perez-Atayde, A.R.; Driscoll, S.G. Postmortem manifestations of echovirus 11 sepsis in five newborn infants. Hum. Pathol. 1983, 14, 818–823. [Google Scholar] [CrossRef]
  20. Modlin, J.F.; Polk, B.F.; Horton, P.; Etkind, P.; Crane, E.; Spiliotes, A. Perinatal echovirus infection: Risk of transmission during a community outbreak. New Engl. J. Med. 1981, 305, 368–371. [Google Scholar] [CrossRef]
  21. Ndiaye, N.; Kébé, O.; Diarra, M.; Thiaw, F.D.; Dia, M.; Dia, N.D.; Sall, A.A.; Fall, M.; Faye, O.; Faye, M. Non-polio enteroviruses circulation in acute flaccid paralysis cases and sewage in Senegal from 2013 to 2021. Int. J. Infect. Dis. 2024, 138, 54–62. Available online: https://www.sciencedirect.com/science/article/pii/S1201971223007816 (accessed on 14 March 2024). [CrossRef]
  22. Fernandez-Garcia, M.D.; Kebe, O.; Fall, A.D.; Ndiaye, K. Identification and molecular characterization of non-polio enteroviruses from children with acute flaccid paralysis in West Africa, 2013–2014. Sci. Rep. 2017, 7, 3808. [Google Scholar] [CrossRef]
  23. World Health Organization. Polio Laboratory Manual; World Health Organization: Geneva, Switzerland, 2004. [Google Scholar]
  24. Lewis, G.D.; Metcalf, T.G. Polyethylene glycol precipitation for recovery of pathogenic viruses, including hepatitis A virus and human rotavirus, from oyster, water, and sediment samples. Appl. Environ. Microbiol. 1988, 54, 1983–1988. [Google Scholar] [CrossRef] [PubMed]
  25. Nijhuis, M.; van Maarseveen, N.; Schuurman, R.; Verkuijlen, S.; de Vos, M.; Hendriksen, K.; van Loon, A.M. Rapid and Sensitive Routine Detection of All Members of the Genus Enterovirus in Different Clinical Specimens by Real-Time PCR. J. Clin. Microbiol. 2002, 40, 3666–3670. [Google Scholar] [CrossRef] [PubMed]
  26. Fujimoto, T.; Izumi, H.; Okabe, N.; Enomoto, M.; Konagaya, M.; Chikahira, M.; Munemura, T.; Taniguchi, K. Usefulness of real-time reverse transcription-polymerase chain reaction for the diagnosis of echovirus aseptic meningitis using cerebrospinal fluid. Jpn. J. Infect. Dis. 2009, 62, 455–457. Available online: https://www.jstage.jst.go.jp/article/yoken/62/6/62_JJID.2009.455/_article/-char/ja/ (accessed on 22 January 2024). [CrossRef]
  27. Majumdar, M.; Martin, J. Detection by direct next generation sequencing analysis of emerging enterovirus D68 and C109 strains in an environmental sample from Scotland. Front. Microbiol. 2018, 9, 1956. Available online: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01956/full (accessed on 21 February 2024). [CrossRef]
  28. Vilsker, M.; Moosa, Y.; Nooij, S.; Fonseca, V.; Ghysens, Y.; Dumon, K.; Pauwels, R.; Alcantara, L.C.; Vanden Eynden, E.; Vandamme, A.-M. Genome Detective: An automated system for virus identification from high-throughput sequencing data. Bioinformatics 2019, 35, 871–873. Available online: https://academic.oup.com/bioinformatics/article-abstract/35/5/871/5075035 (accessed on 13 March 2024). [CrossRef]
  29. Lole, K.S.; Bollinger, R.C.; Paranjape, R.S.; Gadkari, D.; Kulkarni, S.S.; Novak, N.G.; Ingersoll, R.; Sheppard, H.W.; Ray, S.C. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J. Virol. 1999, 73, 152–160. [Google Scholar] [CrossRef]
  30. Martin, D.P.; Varsani, A.; Roumagnac, P.; Botha, G.; Maslamoney, S.; Schwab, T.; Kelz, Z.; Kumar, V.; Murell, B. RDP5: A computer program for analyzing recombination in, and removing signals of recombination from, nucleotide sequence datasets. Virus Evol. 2021, 7, veaa087. Available online: https://academic.oup.com/ve/article-abstract/7/1/veaa087/6020281 (accessed on 11 April 2024). [CrossRef]
  31. Murrell, B.; Moola, S.; Mabona, A.; Weighill, T.; Sheward, D.; Pond, S.K.; FUBAR, K.S. A fast, unconstrained bayesian approximation for inferring selection. Mol. Biol. Evol. 2013, 30, 1196–1205. [Google Scholar] [CrossRef]
  32. Murrell, B.; Wertheim, J.O.; Moola, S.; Weighill, T.; Scheffler, K.; Kosakovsky Pond, S.L. Detecting individual sites subject to episodic diversifying selection. PLoS Genet. 2012, 8, e1002764. Available online: https://journals.plos.org/plosgenetics/article?id=10.1371/journal.pgen.1002764 (accessed on 14 March 2024). [CrossRef]
  33. Kosakovsky Pond, S.L.; Frost, S.D. Not so different after all: A comparison of methods for detecting amino acid sites under selection. Mol. Biol. Evol. 2005, 22, 1208–1222. Available online: https://academic.oup.com/mbe/article-abstract/22/5/1208/1066893 (accessed on 14 March 2024). [CrossRef]
  34. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. In Proceedings of the Nucleic Acids Symposium Series; Oxford University Press: Oxford, UK, 1999; Volume 41, pp. 95–98. [Google Scholar]
  35. Nguyen, L.-T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef] [PubMed]
  36. Rambaut, A. FigTree, a graphical viewer of phylogenetic trees (tree. bio. ed. ac. 691 uk/software/figtree). Inst Evol Biol Univ Edinb. 2014, 1, 2. [Google Scholar]
  37. Rambaut, A.; Lam, T.T.; Max Carvalho, L.; Pybus, O.G. Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evol. 2016, 2, vew007. Available online: https://academic.oup.com/ve/article-abstract/2/1/vew007/1753488 (accessed on 16 July 2024). [CrossRef] [PubMed]
  38. Bouckaert, R.; Heled, J.; Kühnert, D.; Vaughan, T.; Wu, C.-H.; Xie, D.; Suchard, M.A.; Rambaut, A.; Drummond, A.J. BEAST 2: A software platform for Bayesian evolutionary analysis. PLoS Comput. Biol. 2014, 10, e1003537. Available online: https://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1003537 (accessed on 16 July 2024). [CrossRef]
  39. Rambaut, A.; Drummond, A.J.; Xie, D.; Baele, G.; Suchard, M.A. Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst. Biol. 2018, 67, 901–904. Available online: https://academic.oup.com/sysbio/article-abstract/67/5/901/4989127 (accessed on 16 July 2024). [CrossRef]
  40. Celniker, G.; Nimrod, G.; Ashkenazy, H.; Glaser, F.; Martz, E.; Mayrose, I.; Pupko, T.; Ben-Tal, N. ConSurf: Using Evolutionary Data to Raise Testable Hypotheses about Protein Function. Isr. J. Chem. 2013, 53, 199–206. [Google Scholar] [CrossRef]
  41. Faleye, T.O.C.; George, U.E.; Klapsa, D.; Majumdar, M.; Oragwa, A.O.; Adewumi, M.O.; Martin, J.; Adeniji, J.A. Isolation and Genomic Characterization of Echovirus 11 from faeces of a Non-Human Primate in Nigeria. EcoHealth 2020, 17, 461–468. [Google Scholar] [CrossRef]
  42. Brouwer, L.; Moreni, G.; Wolthers, K.C.; Pajkrt, D. World-wide prevalence and genotype distribution of enteroviruses. Viruses 2021, 13, 434. [Google Scholar] [CrossRef]
  43. Kim, H.-J.; Kang, B.; Hwang, S.; Hong, J.; Kim, K.; Cheon, D.-S. Epidemics of viral meningitis caused by echovirus 6 and 30 in Korea in 2008. Virol. J. 2012, 9, 38. [Google Scholar] [CrossRef]
  44. Tassin, M.; Martinovic, J.; Mirand, A.; Peigue-Lafeuille, H.; Picone, O.; Benachi, A.; Vauloup-Fellous, C. A case of congenital echovirus 11 infection acquired early in pregnancy. J. Clin. Virol. 2014, 59, 71–73. Available online: https://www.sciencedirect.com/science/article/pii/S1386653213004812?casa_token=ARcpMvcTHoUAAAAA:_699pnzkTUYbPqayL9a1a8ev6R-pqlj6w5wjyNOFwgy9onC3v2ifcWl3ImOGeMdWwl_xoHWu5pYR (accessed on 15 March 2024). [CrossRef]
  45. Liu, Q.; Xu, W.; Lu, S.; Jiang, J.; Zhou, J.; Shao, Z.; Liu, X.; Xu, L.; Xiong, Y.; Zheng, H.; et al. Landscape of emerging and re-emerging infectious diseases in China: Impact of ecology, climate, and behavior. Front. Med. 2018, 12, 3–22. [Google Scholar] [CrossRef] [PubMed]
  46. Fischer, T.K.; Simmonds, P.; Harvala, H. The importance of enterovirus surveillance in a post-polio world. Lancet Infect. Dis. 2022, 22, e35–e40. Available online: https://www.thelancet.com/journals/laninf/article/PIIS1473-3099(20)30852-5/fulltext (accessed on 18 March 2024). [CrossRef] [PubMed]
  47. Wang, C.; Li, J.; Liu, Y.; Sun, Q.; Liu, Z. Pathogenesis of enterovirus infection in central nervous system. Biosaf. Health 2023, 5, 233–239. [Google Scholar] [CrossRef]
  48. Antona, D.; Kossorotoff, M.; Schuffenecker, I.; Mirand, A.; Leruez-Ville, M.; Bassi, C.; Aubart, M.; Moulin, F.; Lévy-Bruhl, D.; Henquell, C. Severe paediatric conditions linked with EV-A71 and EV-D68, France, May to October 2016. Eurosurveillance 2016, 21, 30402. [Google Scholar] [CrossRef]
  49. Pons-Salort, M.; Oberste, M.S.; Pallansch, M.A.; Abedi, G.R.; Takahashi, S.; Grenfell, B.T.; Grassly, N.C. The seasonality of nonpolio enteroviruses in the United States: Patterns and drivers. Proc. Natl. Acad. Sci. USA 2018, 115, 3078–3083. [Google Scholar] [CrossRef]
  50. Kyriakopoulou, Z.; Pliaka, V.; Amoutzias, G.D.; Markoulatos, P. Recombination among human non-polio enteroviruses: Implications for epidemiology and evolution. Virus Genes 2015, 50, 177–188. [Google Scholar] [CrossRef]
  51. Muslin, C.; Mac Kain, A.; Bessaud, M.; Blondel, B.; Delpeyroux, F. Recombination in enteroviruses, a multi-step modular evolutionary process. Viruses 2019, 11, 859. Available online: https://www.mdpi.com/1999-4915/11/9/859 (accessed on 29 July 2024). [CrossRef]
  52. Oberste, M.S.; Maher, K.; Pallansch, M.A. Evidence for frequent recombination within species human enterovirus B based on complete genomic sequences of all thirty-seven serotypes. J. Virol. 2004, 78, 855–867. [Google Scholar] [CrossRef]
  53. Nikolaidis, M.; Mimouli, K.; Kyriakopoulou, Z.; Tsimpidis, M.; Tsakogiannis, D.; Markoulatos, P.; Amoutzias, G.D. Large-scale genomic analysis reveals recurrent patterns of intertypic recombination in human enteroviruses. Virology 2019, 526, 72–80. [Google Scholar] [CrossRef]
  54. Tee, K.K.; Lam, T.T.-Y.; Chan, Y.F.; Bible, J.M.; Kamarulzaman, A.; Tong, C.Y.W.; Takebe, Y.; Pybus, O.G. Evolutionary Genetics of Human Enterovirus 71: Origin, Population Dynamics, Natural Selection, and Seasonal Periodicity of the VP1 Gene. J. Virol. 2010, 84, 3339–3350. [Google Scholar] [CrossRef]
  55. Marine, R.L.; Ntim, N.A.A.; Castro, C.J.; Attiku, K.O.; Pratt, D.; Duker, E.; Agbosu, E.; Ng, T.F.F.; Gatei, W.; Obodai, E. Strengthening laboratory surveillance of viral pathogens: Experiences and lessons learned building next-generation sequencing capacity in Ghana. Int. J. Infect. Dis. 2019, 81, 231–234. Available online: https://www.sciencedirect.com/science/article/pii/S1201971219300724 (accessed on 18 March 2024). [CrossRef] [PubMed]
  56. Chapman, N.M. Persistent enterovirus infection: Little deletions, long infections. Vaccines 2022, 10, 770. Available online: https://www.mdpi.com/2076-393X/10/5/770 (accessed on 19 March 2024). [CrossRef] [PubMed]
  57. Vignuzzi, M.; Stone, J.K.; Arnold, J.J.; Cameron, C.E.; Andino, R. Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population. Nature 2006, 439, 344–348. Available online: https://www.nature.com/articles/nature04388 (accessed on 17 March 2024). [CrossRef] [PubMed]
  58. Lukashev, A.N.; Vakulenko, Y.A. Molecular evolution of types in non-polio enteroviruses. J. Gen. Virol. 2017, 98, 2968–2981. [Google Scholar] [CrossRef]
  59. Patel, J.R.; Daniel, J.; Mathan, V.I. An epidemic of acute diarrhoea in rural southern India associated with echovirus type 11 infection. Epidemiol. Infect. 1985, 95, 483–492. [Google Scholar] [CrossRef]
  60. Fares, W.; Rezig, D.; Seghier, M.; Ben Yahia, A.; Touzi, H.; Triki, H. Phylogenetic analysis of complete VP1 sequences of echoviruses 11 and 6: High genetic diversity and circulation of genotypes with a wide geographical and temporal range. J. Med. Microbiol. 2011, 60, 1017–1025. [Google Scholar] [CrossRef]
  61. Findlater, A.; Bogoch, I.I. Human mobility and the global spread of infectious diseases: A focus on air travel. Trends Parasitol. 2018, 34, 772–783. Available online: https://www.cell.com/trends/parasitology/fulltext/S1471-4922(18)30142-9 (accessed on 25 July 2024). [CrossRef]
  62. Bek, E.J.; McMinn, P.C. Recent Advances in Research on Human Enterovirus 71. Future Virol. 2010, 5, 453–468. [Google Scholar] [CrossRef]
  63. Sun, H.; Gao, M.; Cui, D. Molecular characteristics of the VP1 region of enterovirus 71 strains in China. Gut Pathog. 2020, 12, 38. [Google Scholar] [CrossRef]
  64. Zhang, Y.; Li, J.; Li, Q. Immune evasion of enteroviruses under innate immune monitoring. Front. Microbiol. 2018, 9, 1866. Available online: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01866/full (accessed on 7 August 2024). [CrossRef]
Viruses 16 01772 g001
Figure 1. Phylogenetic tree of the complete genome sequence (A) (~7381 bp), P1 (~2643 bp) (B), P2 (~2000 bp) (C) and partial P3 (~2250 bp) (D), including 21-279-SEN-2021 and the prototype strain (X80059.1) highlighted in blue and red, respectively. Bootstrap values ≥ 90 are shown on the tree. The scale bar indicates the distances of the branches. The sequences obtained from GenBank included human enterovirus A as an outgroup.
Figure 1. Phylogenetic tree of the complete genome sequence (A) (~7381 bp), P1 (~2643 bp) (B), P2 (~2000 bp) (C) and partial P3 (~2250 bp) (D), including 21-279-SEN-2021 and the prototype strain (X80059.1) highlighted in blue and red, respectively. Bootstrap values ≥ 90 are shown on the tree. The scale bar indicates the distances of the branches. The sequences obtained from GenBank included human enterovirus A as an outgroup.
Viruses 16 01772 g001
Viruses 16 01772 g002
Figure 2. Phylogenetic comparison trees based on (A) 41 entire coding region of the capsid protein sequences (~3800 bp), (B) 68 VP1 sequences (~780 bp) and (C) 70 partial VP1 sequences (~400 bp) from the West African E11 sequences and those obtained from GenBank included human enterovirus A as an outgroup. The characterized sequences are highlighted in blue, purple, green, orange and indigo for Senegal, Guinea, Mauritania, Niger and Burkina Faso, respectively (C). The E11 prototype strain (X80059.1) is highlighted in red and bootstrap values ≥ 90 are shown on the tree. The scale bar indicates the distances of the branches.
Figure 2. Phylogenetic comparison trees based on (A) 41 entire coding region of the capsid protein sequences (~3800 bp), (B) 68 VP1 sequences (~780 bp) and (C) 70 partial VP1 sequences (~400 bp) from the West African E11 sequences and those obtained from GenBank included human enterovirus A as an outgroup. The characterized sequences are highlighted in blue, purple, green, orange and indigo for Senegal, Guinea, Mauritania, Niger and Burkina Faso, respectively (C). The E11 prototype strain (X80059.1) is highlighted in red and bootstrap values ≥ 90 are shown on the tree. The scale bar indicates the distances of the branches.
Viruses 16 01772 g002
Viruses 16 01772 g003
Figure 3. Root-to-tip regression analysis; shown is a linear regression plot for the root-to-tip divergence versus the sampling year (A), Bayesian-inferred maximum clade credibility (MCC) phylogenetic tree based on 47 E11 VP1 (B) and map showing the spatiotemporal and geographical dispersal of the E11 strains included in the present study (C): the grey plots indicate the ocean, blue arrow lines indicate dispersal between West and North Africa, red arrow lines indicate dispersal from West Africa to Europe. The characterized sequences are highlighted in red; the tree branches are colored according to the country of origin. The estimated time to the most recent common ancestor (tMRCA) and posterior probability (PP) are shown in the nodes and the timescale in years is indicated on the x-axis.
Figure 3. Root-to-tip regression analysis; shown is a linear regression plot for the root-to-tip divergence versus the sampling year (A), Bayesian-inferred maximum clade credibility (MCC) phylogenetic tree based on 47 E11 VP1 (B) and map showing the spatiotemporal and geographical dispersal of the E11 strains included in the present study (C): the grey plots indicate the ocean, blue arrow lines indicate dispersal between West and North Africa, red arrow lines indicate dispersal from West Africa to Europe. The characterized sequences are highlighted in red; the tree branches are colored according to the country of origin. The estimated time to the most recent common ancestor (tMRCA) and posterior probability (PP) are shown in the nodes and the timescale in years is indicated on the x-axis.
Viruses 16 01772 g003
Viruses 16 01772 g004
Figure 4. (A) In silico analysis of the amino acid sequences of E11’s VP1 protein. The gaps are indicated by a dash (-) and the conserved amino acid residues by a dot (.). Non-conservative amino acid substitutions in strains recently associated with sepsis in Europe are identified by red boxes. (B) Evolutionary conservation of amino acids in E11’s VP1 protein using the ConSurf server. The amino acids highlighted in the red box indicate the EF loop’s position.
Figure 4. (A) In silico analysis of the amino acid sequences of E11’s VP1 protein. The gaps are indicated by a dash (-) and the conserved amino acid residues by a dot (.). Non-conservative amino acid substitutions in strains recently associated with sepsis in Europe are identified by red boxes. (B) Evolutionary conservation of amino acids in E11’s VP1 protein using the ConSurf server. The amino acids highlighted in the red box indicate the EF loop’s position.
Viruses 16 01772 g004
Table 1. Information on the E11 sequences from West Africa analyzed in this study.
Table 1. Information on the E11 sequences from West Africa analyzed in this study.
Strain ID NumberCollection YearCountryCaseAge (Month)SexIsolated MonthSequence RegionSequence Length (nt)Accession NumberReferenceCo-Infections
13_650_MAU_20132013MauritaniaAFP36MSeptemberCapsid protein3804PP816738This studyNone
13_473_SEN_20132013SenegalAFPN/AN/AN/AVP3-VP1 protein688KY433623[23]None
13_700_MAU_20132013MauritaniaAFP16MSeptemberVP3-VP1 protein876KY433708[23]None
13_647_MAU_20132013MauritaniaAFP60MAugustCapsid protein3896PP8167377This studyNone
13_699_MAU_20132013MauritaniaAFP18MSeptemberCapsid protein3902PP816739This studyNone
13_754_SEN_20132013SenegalAFP24MOctoberCapsid protein4000PP816740This studyNone
13_653_SEN_20132013SenegalAFP36MSeptemberVP3-VP1 protein690KY433626[23]None
14_530_NIG_20142014NigerAFP24MSeptemberCapsid protein3900PP816741This studyNone
14_240_GUI_20142014GuineaAFP24FMayVP3-VP1 protein876KY433711[23]None
14_239_GUI_20142014GuineaAFP36MMayVP3-VP1 protein876KY433710[23]None
14_381_NIG_20142014NigerAFP19MJuneVP3-VP1 protein876KY433712[23]None
14_128_NIG_20142014NigerAFP17FMarchVP3-VP1 protein876KY433709[23]None
17_787_SEN_20172017SenegalAFPN/AN/AN/AVP3-VP1 protein711OM827233[23]None
17_718_SEN_20172017SenegalAFPN/AN/AN/AVP1 protein783PP816742This studyNone
21_279_SEN_20212022SenegalAFP48MMarchComplete genome7357PP816743This studyNone
22_438_SEN_20222022SenegalAFPN/AN/AN/ACapsid protein3900PP816746This studyNone
22_081_SEN_20222022SenegalAFP24FNovemberCapsid protein4003PP816744This studyNone
22_102_SEN_20222022SenegalAFP24MMayCapsid protein3800PP816745This studyNone
22_751_SEN_20222022SenegalAFPN/AN/ASeptemberCapsid protein3800PP816747This studyNone
23_371_BFA_20232023Burkina FasoESN/AN/AJulyCapsid protein4000PP816749This studyE19
23_399_SEN_20232023SenegalESN/AN/AJulyCapsid protein4000PP816750This studyE19
23_368_SEN_20232023SenegalESN/AN/AJulyVP1 protein804PP816748This studyE19
23_415_BFA_20232023Burkina FasoESN/AN/AJulyCapsid protein3987PP816751This studyE12
F = female; M = male; N/A = not applicable; nt = nucleotide.
Table 2. Episodes of positive diversifying selection on the echovirus 11 proteins.
Table 2. Episodes of positive diversifying selection on the echovirus 11 proteins.
ProteinsNumber of Sites Detected by MethodEvidence of Positive Selection
SLAC
(p < 0.1)		FUBAR
(Posterior Probability ≥ 0.9)		MEME
(p < 0.1)
VP1Sites under negative selection (dN/dS < 1)1341550NO
Sites under positive selection (dN/dS > 1)000
CapsidSites under negative selection (dN/dS < 1)87711500NO
Sites under positive selection (dN/dS > 1)006
SLAC = Single-likelihood ancestor counting; FUBAR= Fast unconstrained Bayesian approximation. MEME = Mixed effects model of evolution; dN = non-synonymous substitutions; dS = synonymous substitutions; VP1, Capsid = Viral proteins.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ndiaye, N.; Thiaw, F.D.; Lagare, A.; Sinare, T.; Diakité, M.L.; Ngom, S.F.M.; Kébé, O.; Abdoulkader, I.K.; Cissé, G.; Dia, M.; et al. Recent Molecular Epidemiology of Echovirus 11 Throughout North and West Africa Resulted in the First Identification of a Recombinant Strain from an Acute Flaccid Paralysis Case in West Africa. Viruses 2024, 16, 1772. https://doi.org/10.3390/v16111772

AMA Style

Ndiaye N, Thiaw FD, Lagare A, Sinare T, Diakité ML, Ngom SFM, Kébé O, Abdoulkader IK, Cissé G, Dia M, et al. Recent Molecular Epidemiology of Echovirus 11 Throughout North and West Africa Resulted in the First Identification of a Recombinant Strain from an Acute Flaccid Paralysis Case in West Africa. Viruses. 2024; 16(11):1772. https://doi.org/10.3390/v16111772

Chicago/Turabian Style

Ndiaye, Ndack, Fatou Diène Thiaw, Adamou Lagare, Thérèse Sinare, Mohamed Lemine Diakité, Serigne Fallou Mbacké Ngom, Ousmane Kébé, Issifi Kollo Abdoulkader, Gassim Cissé, Mohamed Dia, and et al. 2024. "Recent Molecular Epidemiology of Echovirus 11 Throughout North and West Africa Resulted in the First Identification of a Recombinant Strain from an Acute Flaccid Paralysis Case in West Africa" Viruses 16, no. 11: 1772. https://doi.org/10.3390/v16111772

APA Style

Ndiaye, N., Thiaw, F. D., Lagare, A., Sinare, T., Diakité, M. L., Ngom, S. F. M., Kébé, O., Abdoulkader, I. K., Cissé, G., Dia, M., Djimadoum, H. N., Neya, C. O., Boubakar, R., Ouedraogo, I., Essoya, L. D., Dia, N., Sall, A. A., Faye, O., & Faye, M. (2024). Recent Molecular Epidemiology of Echovirus 11 Throughout North and West Africa Resulted in the First Identification of a Recombinant Strain from an Acute Flaccid Paralysis Case in West Africa. Viruses, 16(11), 1772. https://doi.org/10.3390/v16111772

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop