
==== Front
Int J Microbiol
Int J Microbiol
IJMICRO
International Journal of Microbiology
1687-918X
1687-9198
Wiley

10.1155/2024/7796913
Research Article
Genetic Diversity of Human Enterovirus in Kazakhstan, during 2022
Kamalova Dinara 1 2
Akhmetova Assel 1
Amirgazin Asylulan 1
Sytnik Igor 1
Rudenko Viktoriya 3
Yessimkhanova Gulzhan 3
Bayesheva Dinagul 3
Yegorov Sergey 1 4
https://orcid.org/0000-0002-0307-1053
Shevtsov Alexander ncbshevtsov@gmail.com
1
1 National Center for Biotechnology, Astana, Kazakhstan
2 L. N. Gumilyov Eurasian National University, Astana, Kazakhstan
3 Astana Medical University, Astana, Kazakhstan
4 Department of Biology School of Sciences and Humanities Nazarbayev University, Astana, Kazakhstan
Academic Editor: Todd R. Callaway

2024
26 8 2024
2024 779691322 12 2023
12 6 2024
15 7 2024
Copyright © 2024 Dinara Kamalova et al.
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Enteroviral infection is a common cause of aseptic meningitis, herpangina, and hand, foot, and mouth disease in children. Limited data are available on the enteroviral subtypes associated with hospitalization for these conditions in Kazakhstan. We collected cerebrospinal fluid (CSF) and nasopharyngeal swabs (NSW) from children (N = 152, median age = 8 years) hospitalized with symptoms of aseptic meningitis (AM, N = 139) or herpangina (HA, N = 13) disease. We then genotyped enteroviral subtypes associated with AM (n = 50) and HA (n = 9) using next-generation sequencing (NGS) on the viral protein 1 (VP1), followed up by whole-genome sequencing of the isolated viral species. All identified EVs were species B EV, consisting of five echoviruses (E6, E9, E11, E21, and E25) and three coxsackieviruses (CVA9, CVB3, and CVB5) serotypes within the cohort. The most abundant EVs were CVA9 (38.5%), CVB5 (21.5%), and E6 (13.8%). Most HA samples (6/9) were genotyped with coxsackievirus CVA9, while AM was associated with a variety of both echovirus and coxsackievirus serotypes. The results suggest that coxsackievirus CVA9 may be the dominant serotype circulating in the HA population, while AM is more diverse in terms of circulating echovirus and coxsackievirus serotypes. Further studies are needed to determine the clinical implications of these findings and to investigate potential differences in disease severity or outcomes associated with different EV serotypes.

Ministry of Education and Science of the Republic of KazakhstanAP08052252
==== Body
pmc1. Introduction

Human enteroviruses (HEVs) are small single stranded RNA viruses of the family Picornaviridae (https://www.picornaviridae.com) that cause different types of infectious illnesses in humans affecting millions of people every year in the world. Symptoms can vary from mild to severe and can result in syndromes like respiratory diseases, herpangina, hand-foot-and-mouth disease, aseptic meningitis, acute hemorrhagic conjunctivitis, and myocarditis. [1, 2]. Herpangina and meningitis in children are significant illnesses caused by enteroviruses because of the potential severity and impact on children's health. Herpangina is highly contagious and spread very quickly, and it mostly affects younger children and is characterized by fever, sore throat, and small ulcers in the mouth and at the back of the throat [3, 4]. The ulcers can be painful and cause discomfort for children when drink and eat, which increases the risk of dehydration [3]. Symptoms can be confused with other illnesses; therefore, accurate diagnosis is significant for proper management and treatment [1, 4]. Herpangina complications are rare; however, in some cases, they can cause brainstem encephalitis and acute flaccid paralysis.

Based on serological tests, enteroviruses (EVs) were originally divided into coxsackieviruses A and B, echoviruses, and polioviruses subgroups [5]. Species of the genus Enterovirus are divided based on their pathogenicity, to human (Enterovirus A–D and Rhinovirus A–C) and animal species (Enterovirus E–J) [6]. However, with the development of molecular biology methods, classification is based on genetic characteristics and phylogenetic relationships between strains [7], where each EV species includes different serotype groups. Molecular typing of HEV serotypes is important for many reasons, such as in case of severe clinical disorders (some serotypes can cause fatal neurological disease), for the discrimination of polio and nonpolio enteroviruses, as well as for epidemiological and phylogenetic studies (discovering newly occurring serotypes) [8].

The EVs structure comprises nonenveloped positive stranded RNA molecule (approximately, 7400 nt long), which is surrounded by the structural proteins VP1 to VP4 (each in sixty copies) that form icosahedral capsid structure of the virus [9]. The genome of enterovirus is a single RNA molecule with small virus-encoded protein VPg at the 5′ end and untranslated region (UTR) from both 5′ and 3′ ends [10]. Many research studies have shown that VP1 capsid protein is a good marker for enterovirus serotypes' differentiation and is widely used for PCR-based typing [10–12]. Recently known detection and classification methods are based on reverse-transcription polymerase chain reaction (RT-PCR) [13].

The HEV-B species consists of the vast majority of EV serotypes, and according to King's [14] virus taxonomy, it includes coxsackievirus A9 (CVA9), coxsackieviruses B (CVB) 1–6, echoviruses (EV) 1–7, 9, 11–21, 24–27, 29–33, and enterovirus 69. Moreover, only coxsackievirus A9 from all coxsackieviruses belongs to HEV-B, while others belong to HEV-A and HEV-C [14].

Enteroviruses can be transmitted among humans in different ways: direct contact (fecal-oral way) and indirect contact (contaminated equipment) or drinking water from infected reservoirs. Infections predominate in the summer and autumn seasons, however, can occur sporadically during the year. According to previous research, younger children are known to be affected more often comparing to adults [15, 16]. Although, in most cases, infections caused by EVs occur with mild symptoms as respiratory or gastrointestinal illnesses, some of them can cause severe complications for patients, from central nervous system invasion to meningitis, encephalitis, and paralysis [17].

Understanding the molecular signatures of HEVs is essential for effective diagnosis and treatment protocols. The situation with HEV infection molecular typing in Kazakhstan is not yet fully studied, and there is no whole genome data available yet. This research study describes clinical diseases and associated HEV serotypes, as well as whole genome sequencing of 65 EV isolates collected from children in central infectious diseases hospital, located in Astana in 2022. The aim of the study is to investigate the molecular epidemiology of occurring HEV infection and the prevalence of circulating serotypes in Kazakhstan. In this study, we analyzed 65 sequences of HEV and compared with the publicly available NCBI database to understand important evolutionary dynamics and infection spread.

2. Materials and Methods

2.1. Patients and Selection Criteria

This study provides a retrospective analysis of 152 cases of hospitalization of children in the 3rd Children's Infectious Diseases Hospital in Astana from May 14 to October 12, 2022, with a preliminary diagnosis of enterovirus infection. Children's Infectious Diseases Hospital Astana is the only specialized institution providing inpatient treatment for childhood infections in Astana (population 1.35 million people).

Enterovirus infections were classified by clinical phenotypes, using the following definitions: herpangina (HA) and aseptic meningitis (AM) caused by enterovirus.

Specimens for the current study were collected, stored, and transported according to the CDC guidelines (https://www.cdc.gov/non-polio-enterovirus/lab-testing/specimen-collection.html). Patient specimens for “rapid” molecular diagnostics were processed in the rapid laboratory, while samples for “delayed” molecular diagnostics and “routine” diagnostics were stored at −20°C.Nasopharyngeal and nasal swabs were collected within the first three days of illness. Swabs were collected using a dry, sterile probe with a cotton swab, which, after collection, was placed in a disposable tube containing 500 μl of transport medium (0.9% sodium chloride solution or phosphate buffer).

Cerebrospinal fluid (CSF) was collected under aseptic conditions during the first days of the disease when clinical indications for lumbar puncture were present. A volume of 1 ml was collected in a sterile disposable plastic tube.

Samples from patients selected for “rapid” molecular diagnostic group were processed in the rapid diagnostic laboratory as soon as they were received. Specimens for “delayed” molecular diagnostic testing from subjects randomized to the “conventional” diagnostic group were stored at −20°C. Delivery to the laboratory was made within 12 hours in thermal containers with refrigeration at a temperature not exceeding +4°C.

Herpangina and hand-foot-mouth diseases (HFMD) cause similar symptoms with ulcers in or around the mouth, with HFMD causing vesicular lesions on the arms, legs, knees, or buttocks. A nasopharyngeal swab (NSW) was taken from patients using a sterile cotton swab and placed in a transport medium (DMEM with 50 μg/ml of gentamicin) and stored at minus 70°C until RNA extraction.

Aseptic meningitis was accompanied by an increase in body temperature up to 38−40°C, the development of neck rigidity, headaches, and photophobia was observed, and some patients experienced vomiting, loss of appetite, diarrhea, rash, pharyngitis, and myalgia. For patients with signs of meningitis, according to the standard protocol, a cerebrospinal fluid (CSF) puncture was prescribed, which was divided into 3 equal aliquots by 200 μl. According to a standard diagnostic protocol, two CSF sample aliquots were used for enterovirus diagnosis. One sample was analyzed using bacteriology and cytology methods at the hospital laboratory. The second aliquot was sent for laboratory analysis using RT-PCR kit AmpliSens® Human Enterovirus-FL (AmpliSens, Russia) in the RSE “Center for Sanitary and Epidemiological Expertise” Medical Centre Hospital of President's affairs Administration of the Republic of Kazakhstan. The third part of the CSF was stored at minus 80°C for additional studies and genotyping. In case of presence of vesicular lesions of the oral cavity, a nasopharyngeal swab was also collected from patients.

In total, for the specified period, 139 children were hospitalized with aseptic meningitis diagnosis confirmed and the presence of enterovirus RNA in cerebrospinal fluid samples detected. Herpangina was diagnosed in 13 children patients, and enterovirus RNA detected from NSW samples.

2.2. RNA Extraction

Viral RNA was extracted from 200 μl CSF or NSW samples in transport medium using GeneJET kit Viral DNA and RNA Purification Kit (ThermoScientific, Lithuania) according to the manufacturer's instructions.

2.3. Whole Genome Sequencing and Assembly of Genomes

Almost full-length genomes of enteroviruses were amplified as described previously by Isaacs et al. [18], except that the readily made reaction mix from BioMaster LR HS-PCR 2x (Biolabmix, Russia) was used for amplification step. Second-stage PCR products were purified with AMPure XP magnetic beads (Beckman Coulter, USA) in a 1 : 1 ratio, with elution in 20 μl of DNase/RNase-free distilled water. Libraries were prepared using the Illumina® DNA Prep kit, (M) Tagmentation (96 Samples) (Illumina, Catalog #20018705) with double barcoding. Libraries were sequenced on Illumina MiSeq using MiSeq Reagent Kit v3 (600-cycle) (Catalog #LMS -102-3003). Reagents were used according to the manufacturers' instructions.

The quality control of the obtained reads was carried out using the FastQC program [19]. Before denovo assembly, reads were trimmed using SeqTK (with options: -b 20-e3) and sickle (with options: -t sanger-q 30-l 200-g), followed by deduplication in the FastP program (with options: -- dedup -- dup _ calc _ accuracy 6 -- disable _ quality _ filtering). To assemble the genomes, Megahit v1.2.9 was used with varying the length of κ-mers individually for each sample, selecting the assembly with the longest lengths for subsequent contigs analysis. The resulting contigs were then used as corresponding references to obtain consensus sequences using BWA. Variants were determined using FreeBayes, and consensuses were generated by BCFtools consensus. The ends of assemblies with less than 30x coverage were removed.

2.4. Determination of Serotypes

Serotypes were determined by constructing a phylogenetic tree from the database generated with VP1 sequences of enteroviruses from GenBank (https://www.ncbi.nlm.nih.gov/genbank/). The database was selected using keywords such as “enterovirus,” “VP1,” and “human enterovirus B.” In total, approximately 26,000 genome sequence data with VP1 regions were downloaded. Using the generated database and our sequences, a phylogenetic tree was reconstructed, which made it possible to filter the 81 isolates most closely related to our sequences from the database. Further analysis was performed with 146 enterovirus isolates, including 65 obtained within the current study, and Enterovirus 69 was selected as an outgroup (GenBank: AY302560.1) [20]. The statistical selection of the nucleotide substitution model for phylogenetic analysis was performed using ModelTest -NG (v0.1.7) [21], and the best fitting model was selected for the reconstruction of the phylogenies. Phylogenetic tree was built using the maximum likelihood (ML) method with the bootstrap random sampling method, n = 1000. Estimated bootstrap (bs) values greater than 50% were shown for the tree nodes. The Interactive Tree of Life (iTOL) [22] online tool was used for interpretation, visualization, and annotation of the phylogenetic tree.

2.5. Genotyping by VP1 Region

Representative datasets were downloaded separately for each serotype group to describe the molecular epidemiology of HEV-B serotypes in Kazakhstan. Phylogenetic relatedness of enterovirus isolates collected within the current study with the most recent genotyping and subgenotyping data was compared. Phylogenetic dendrograms were reconstructed based on entire VP1 sequences of the international strains downloaded from GenBank.

Genotyping of coxsackievirus A9 (CVA9) was based on Zhao et al. [23] algorithm, that included 110 international sequences and 25 obtained within the current study. Coxsackievirus B5 (CVB5) genotyping of 134 isolates, including 14 Kazakhstani, was performed using the same algorithm as described in He et al. [24]. Coxsackievirus B3 (CVB3) strains were analyzed based on Yang et al. [25], including 250 worldwide strains and 1 isolate from Kazakhstan. Genetic diversity of echoviruses group worldwide datasets and the relatedness of Kazakhstani isolates were reconstructed according to Cheng et al. [26] for Echovirus 6 (E6) (n = 73); Zhang et al. [27] for Echovirus 9 (E9) (n = 57); Li et al. [28] and Grapin et al. [29] for Echovirus 11 (E11) (n = 140). However, there are no currently available genotyping schemes for the Echovirus 21 (E21) and 25 (E25) groups. Therefore, entire VP1 coding sequences of available E21 and E25 serotypes were downloaded from GenBank (for August 15th, 2023) and phylogenetic dendrograms reconstructed, including sequences obtained within the current study. Genotypes were distinguished based on estimated divergence between/within groups, calculating pairwise genetic distances in MEGA X [30]. Genotypes determined with >15% difference between and <15% within groups.

The pairwise genetic distances for E21 and E25 serotypes were estimated within each serotype group using the Kimura 2-parameter model, with the gamma distribution model (shape parameter = 5) in MEGA X [30]. Enterovirus sequences were classified into so-called genotypes within each serotype (shown in Figures S1–S5) if they shared >85% sequence identity based within the VP1 capsid gene.

Sequences were aligned in Mafft (v7.520) [31], and the best-fitting substitution model was selected (for each serotype group) using the ModelTest-NG (v0.1.7) [21] program. Maximum likelihood tree reconstructed for VP1 sequences with bootstrapping of replicates 1000 times using RAxML-NG (v1.2.0). The ML tree was reconstructed to demonstrate the genetic diversity of the worldwide strains and the allocation of sequenced Kazakhstani enterovirus strains to particular genotype and subgenotype group. The Interactive Tree of Life (iTOL) (v1.0) [22] was used for illustration and annotation of the ML tree.

2.6. Analysis of Whole Genome Sequencing Data

Phylogenetic analysis of 65 whole genome sequences from 58 infected patients in Kazakhstan was performed. The nucleotide substitution model for evolutionary analysis was performed in ModelTest-NG (v0.1.7) [21], and the best-fitting model was selected as “GTR + I + G.” A phylogenetic tree was reconstructed using the maximum likelihood (ML) method with the bootstrap random sampling method in RAxML-NG (v1.2.0), with n = 1000 bootstraps. Statistically significant bootstrap (bs) values greater than 50% were shown for the tree nodes. To detect if a sufficient “temporal signal” was present in the data, TempEst (v1.5.3) [32] software was used. Phylogenetic tree was visualized and annotated in the Interactive Tree Of Life (iTOL) [22] online tool. Pairwise genetic distances were estimated using the ape package in R [33].

3. Results

3.1. Samples' Overview

During the period from May 14 to October 12, 2022, 152 patients with a preliminary diagnosis of enterovirus infection were admitted to the Children's Infectious Diseases Hospital No. 3 in Astana. Schematic representation of the samples collected is shown in Figure 1.

The largest number of patient visits occurred during the summer period, as shown in Figure 2; 67% of the studied patients were hospitalized in June-July. There was a decrease in the number of requests during August, followed by a peak increase in September (23.7%). The number of children under 10 years of age predominates and equals 74.3% (n = 113) of the total number of admitted patients, compared with the number of children aged over 10 years of 25.7% (n = 39) (average age of patients: 8.1 years). Moreover, among patients under the age of 10 years, more than 59% of patients were male (n = 67). Around 91.5% (139/152) of children were hospitalized with a diagnosis of aseptic meningitis and with a diagnosis of herpangina 8.5% (13/152).

Laboratory cytology results were available for 131 of 152 patients. The median value of blood leukocytes was 10.36 cells/L, with 25th and 75th percentile values of 7.99 cells/L and 13.27 cells/L, respectively. The median percentage of lymphocytes and neutrophils in CSF samples were 53% and 42%, respectively, with the values of the 25th and 75th percentiles being 20% and 77% for lymphocytes, and 22.5% and 79% for neutrophils, respectively.

4. Results of Whole-Genome Sequencing and Identification of Serotypes

Nearly complete whole-genome sequences of 65 enterovirus isolates were amplified and assembled from RNA samples collected from 58 patients. From all patients, 50 were diagnosed with aseptic meningitis and 8 with enterovirus infection herpangina. From 50 meningitis patients, 32 CSF samples and 24 NSW samples were selected for whole genome sequencing analysis. Nine sequences were obtained from NSW samples from 8 patients with enterovirus herpangina (Table 1).

Consensus phylogenetic tree was built using the ML algorithm for the VP1 sequence region (Figures 3(a) and 3(b)). Phylogenetic clades were mostly separated by enterovirus serotypes, which is consistent with previous enterovirus research studies [34].

These strains were labeled using the following format: “GenBank accession number”/“Serotype”/“Country of origin”/“Year of isolation” (strains with partially available metadata were also included). The prototype Enterovirus B69 was used as an outgroup. Bootstrap values >50% are shown as blue circles.

Based on phylogenetic analysis of the VP1 gene, 32 sequences obtained from CSF samples of 29 patients with meningitis were clustered into 6 serotypes: CVA9, CVB5, E6, E9, E11, and E21 (Figure 4). 24 sequences obtained from NSW samples of 21 patients with meningitis were clustered into 7 serotypes: CVA9, CVB5, CVB3, E6, E9, E21, and E25. Two CSF samples were found to be infected with both E21/CVB5 and E11/CVB5, respectively. In one patient, CVA9 was identified in a collected CSF sample and two serotypes were identified in the NSW sample (E9/CVA9) (Table 1).

Of the 9 sequences from 8 patients with enterovirus infection, 6 were identified as CVA9, 2 as CVB5, and 1 as E6, with one sample found to be coinfected with CVB5/E6 (Table 1, Figure 4).

4.1. Genotyping of Strains by the VP1 Gene

Subsequently, we analyzed genotypes based on the VP1 gene among certain serotypes. The obtained sequences were grouped with Asian and European genotype lineages (Figures S1–S8). The 25 coxsackievirus A9 samples clustered as a separate clade into genotype I and are genetically closest to the 2013 Russian sequence. 14 Coxsackievirus B5 sequences formed a cluster in genotype D and are genetically close to isolates circulating in China in 2015–2018. Nine echovirus 6 sequences cluster into genotype E and are genetically close to European and American lineages. Seven Echovirus 9 sequences cluster into a separate clade within genotype F, together with European lineages. Three Echovirus 11 sequences were included in the D5 subgenotype. This subgenotype also includes lineage 1 strains circulating in France and Spain in 2022 and 2023 and associated with neonatal infection with liver failure [29, 35, 36]. This lineage was clustered within a separate clade in subgenotype D5. A single coxsackievirus B3 sequence was clustered into the most prevalent genotype E.

For the analysis of Echovirus 21 genotypes, 31 complete VP1 sequences were downloaded from GenBank (sequences deposited until August 15, 2023). Current analysis included 3 sequences from Kazakhstan obtained as part of this study collected from CSF and NSW samples from patients with aseptic meningitis. Based on the phylogenetic dendrogram and percentage of identity between strains, E21 serotype can be divided into 10 genotypes A–J, the maximum difference in nucleotide sequences within one genotype was 12%, and the minimum average difference in nucleotide sequences between genotypes was 17% (Tables S1 and S2). Genotype A is represented by one echovirus 21 prototype strain Farina. Genotypes B, C, and F are represented by single sequences from Madagascar, China, and France. For genotypes D, E, and I, geographic relationships with strains from China, Central African Republic, and India are monitored. Genotype J combines strains from Asia and Europe. Three sequences from Kazakhstan were included in the H genotype, together with strains from China and the United Kingdom.

For the analysis of Echovirus 25 genotypes, 91 complete VP1 sequences were obtained from GenBank (sequences deposited until August 15, 2023). The analysis included 3 sequences from Kazakhstan obtained as part of this study. E25 can be divided into 8 genotypes A–H, and the identity of nucleotide sequences in the genotypes exceeded 85%. Kazakhstani strains are included in genotype A and are in a separate clade, which may subsequently become a subgenotype. Genotypes B and C include strains isolated from China and the USA, respectively. Genotypes D and H are the largest and combine strains circulating on the Eurasian and American continents. Genotype E combines two strains from the USA and Nigeria. Genotypes F and G include strains isolated in Europe.

4.2. Phylogenetic Analysis by Complete Coding Sequence

Based on phylogenetic analysis of the complete coding sequence, 65 sequences were clustered according to specific serotypes (Figure 5).

The 25 CVA9 sequences formed three separate clusters, and the third cluster was represented by a single sequence. We did not find any features in the distribution of sequences in the clades in comparison with the sampling time. For example, sequences collected in June are evenly represented across all contributions, including a third clade represented by a single sequence from NSW from a patient with meningitis that is most genetically distant from all of them by over 3.7% (Table S3).

The CVA9 sequences of EV-H1-vir2 and EV-38 isolated from the same patient with CSF and NSW were 100% identical. The 14 CVB5 sequences formed two clusters, representing sequences collected from patients that were hospitalized in different periods. All 9 E9 sequences were collected in June and formed two clusters, differing from each other by 8% (Table S4).

The E6 sequences are the most genetically homogeneous, with less than 2% variability among themselves (Table S5).

5. Discussion

In Kazakhstan, cases of enterovirus infection are recorded annually, with peak values in the summer-autumn period [37]. In our study, more than 90% of hospitalized cases of enterovirus infection were due to aseptic meningitis, and only 8.5% were due to enterovirus herpangina. These data do not reflect the true picture of the spread of enterovirus HA in Kazakhstan since the classic mild course of herpangina may remain unreported due to home or local hospital treatment. The high number of cases of aseptic meningitis in our study correlates with previous study data, in which enteroviruses were identified as etiological agents in 73% of cases of encephalitis and meningitis in Kazakhstan, with an average incidence of 14 per 100,000 population, and for children under 15 years of age, the incidence is 35.9 per 100,000 children [38]. Despite the high incidence rate, information on circulating serotypes is limited. This is the first study to describe the genetic features of enteroviruses circulating in Kazakhstan based on nearly complete whole-genome data.

The highest incidence was recorded in the summer months, and more than 67% of patients were hospitalized in June-July months, followed by a decrease in hospitalization cases in August and a rise in incidence in September. Peaks of incidence in the summer-autumn period are typical for the northern hemisphere [39]. Of the 152 patients, an almost complete genome sequence of enteroviruses was obtained in 38% (58 patients), while real-time PCR confirmed the presence of enterovirus infection in all patients. The lower efficiency of genome-wide amplification in comparison with the original article [18] may be due to the lack of information about the viral load in the studied samples, RNA degradation, and differences in the amplification reagents used. The assembly revealed coinfection in 6 of the 58 patients. In patients with aseptic meningitis, the presence of two serotypes was established in 2 samples from CSF (E21/CVB5 and E11/CVB5), 2 samples from NSW (E6/CVA9), and in one patient, 3 serotypes were obtained (CSF - Coxsackievirus A9 and Echovirus E9/Coxsackievirus A9 in a sample from NSW), while the complete coding sequences of Coxsackievirus A9 from NSW and CSF were completely identical. The patient with HA was found to be coinfected with NSW CVB5/E6. The introduction of molecular genetic tests and sequencing schemes for genotyping enteroviruses has expanded our understanding of the epidemiology of enteroviruses and led to the discovery of coinfection. Coinfection is most often detected in patients with HFMD syndrome and can reach 34%, while simultaneous persistence of different serotypes can aggravate the severity of the infection and lead to atypical cases involving the respiratory and cardiovascular systems [40–42]. In the current study, we did not record any specific clinical course in patients with coinfections, which may be due to the limited sample size.

The most prevalent enterovirus genotype circulating in Astana during the seasonal rise of enterovirus infection in 2022 is CVA9, which was detected in 18 of 50 patients with aseptic meningitis and 6 of 8 patients with herpangina. Coxsackievirus A9, like many enteroviruses, exhibits a wide range of clinical manifestations, including aseptic meningitis (AM), hand, foot, and mouth disease (HFMD), acute flaccid paralysis (AFP), and persistent diarrhoea. Generalized febrile exanthema can cause pathologies in the placenta in pregnant women and severe generalized infections in newborns [43–46]. In some outbreaks of aseptic meningitis, CVA9 was the main serovariant in China [47], Canada [48], and South Africa [49]. CVA9 belongs to the virus serovars that were most often isolated in serous meningitis in Russia, while in 2022, there was an increase in the incidence of meningitis associated with CVA9 [50, 51]. CVA9 isolated in our study belongs to genotype I, which is mainly represented by Russian strains [23]. Kazakhstan and Russia have the longest land border in the world and close economic and cultural ties, which may facilitate the circulation of identical genotypes.

Coxsackievirus B5 in our study was the second most abundant enterovirus identified in 12 people with aseptic meningitis and in 2 patients with enterovirus herpangina. In global distribution, CVB5 is not a common serovariant in severe neurological complications, including acute flaccid paralysis, encephalitis, and aseptic meningitis [52]. However, in Poland [53] and China [24, 54, 55], CVB5 has been the predominant serovariant in several outbreaks in recent years. As shown in Farshadpour and Taherkhani [56], enterovirus was identified as a main cause of aseptic meningitis among studied population in Iran. Echovirus 30 was emerged as a predominant strain that caused aseptic meningitis in children in their first year of life [56]. Genotyping clustered the sequences of the current study into genotype D, which is predominantly represented by sequences from China. Given the migration pattern of CVB5 from Northeast China to Northwest China and East China to Northwest China [24], it is tempting to hypothesize a shift in the CVB5 circulation boundaries from China to Kazakhstan.

Echoviruses that were predominant in our study were E6 and E9. Echovirus 6 is one of the five enteroviruses that were diagnosed as the etiological agent of neurological complications, including acute flaccid paralysis (AFP), encephalitis, and aseptic meningitis [52, 57, 58]. In our case, E6 was detected in 8 patients with aseptic meningitis, of which only one in CSF and one case with herpangina, while coinfection with other enteroviruses was identified in 4 patients. High levels of E6 coinfection may promote recombination changes, requiring monitoring of the serovariant at the genetic level and tracking changes in virulence [59].

Nevertheless, official statistical data indicate a seasonal increase in enterovirus infections in 2022. Regular genotype studies in enterovirus infections will improve the control and differentiation of outbreaks from seasonal increases and assess the prevalence of dominant genotypes and the association of the clinical course with enterovirus genotypes in Kazakhstan. Analysis of global data shows that echovirus 30, echovirus 6, and echovirus 13 are most commonly detected in enterovirus meningitis [52]. Echoviruses also tend to predominate in meningitis cases recorded during endemic periods. A four-year follow-up in Poland shows the predominance of echovirus 6 and echovirus 30 in enterovirus meningitis cases [60]. In China, E6, E18, E11, E9, E30, and CVB5 were identified as the most common genotypes in aseptic meningitis between 2009 and 2018 [61–63]. In Latin American countries, aseptic meningitis is also usually dominated by echoviruses E30 and E6 [64]. In our case, the predominance of the CVA9 genotype is observed in aseptic meningitis and herpangina, which is usually characteristic of outbreaks [47–49]. Nevertheless, official statistical data indicate a seasonal increase in enterovirus infections in 2022. Regular genotype studies in enterovirus infections will improve the control and differentiation of outbreaks from seasonal increases and assess the prevalence of dominant genotypes and the association of the clinical course with enterovirus genotypes in Kazakhstan.

Current study had several limitations, such as the strategy of laboratory diagnosis protocols in hospital laboratories requires sending the samples for molecular detection analysis using RT-PCR, which in some cases is very time consuming, and doctors might have to wait for the results before the treatment. Reliable rapid diagnosis of enterovirus infections is essential in children; therefore, the use of sensitive screening tests for EV infection detection is important in clinical settings and hospitals. However, in most cases, symptoms are mild in patients' with enteroviral infections and require only supportive medications aimed at relieving the symptoms because it is challenging to distinguish between viral infections with similar symptoms based only on clinical representation. Given the prevalence of asymptomatic EV infections, merely detecting EV in a specimen does not necessarily indicate the cause of disease.

Identifying enterovirus serotypes is challenging and requires sequencing of the strains or sequencing of the VP1 gene, which is not available in some regions of the country. Genotyping can be more informative compared to serotyping, and utilizing the nucleotide sequence of VP1 can serve as a proxy for antigenic typing, aiding in distinguishing between enterovirus serotypes. This study aimed to explore the genetic diversity of enterovirus circulating in the area and to give a broader understanding of the processes underlying this transmission of enteroviral infection among children in Kazakhstan.

6. Conclusions

The reason for the seasonal rise in enterovirus infection in Astana in 2022 is eight serovars of group B enteroviruses; regular monitoring with genotyping is essential to understand any changes in their population structure. The central Eurasian location of Kazakhstan facilitates the circulation of both European and Asian genotypes of enteroviruses. The detection of coinfection in CSF and NSW, as well as the identification of different serotypes in the patients' CSF and NSW samples, suggests the possibility of misidentification of the etiological agent of aseptic meningitis during virological examination of NSW and stool samples. The introduction of NGS genotyping will expand knowledge about the possible coinfection of different enteroviruses and, consequently, recombination patterns in them.

Acknowledgments

This work was supported by the Ministry of Education and Science of the Republic of Kazakhstan to frame project AP08052252 for 2021–2023 to Igor Sytnik.

Data Availability

Sequences obtained within the current study were submitted to GenBank under OR678399-OR678463 accession numbers.

Ethical Approval

The clinical research described in the article was approved by the local bioethical committee of Astana Medical University under approval number 6, dated January 31, 2022. The study was conducted in accordance with the principles outlined in the Helsinki Declaration of 1975 and its amendments from 2005, which ensure ethical guidelines for medical research involving human subjects.

Consent

Written informed consent was obtained from the parents or guardians of the children participated in the study. Laboratory tests for hospitalized patients were conducted within the framework of standard protocol of diagnostics and treatment approved by the Joint Committee on the quality of medical services of the Ministry of Health of the Republic of Kazakhstan from “16” October 2020, Protocol no. 117.

Disclosure

Our study is also aligned with these goals. In our comprehensive study, we concluded all available high-quality clinical studies, and the presented research is the first study of the genetic variation and genotypic diversity of the main pathogen of viral meningitis in Central Eurasia. Our results will influence future diagnostics and supposedly request more specific clinical testing for coxsackieviruses, which appeared to be the predominant causes of aseptic meningitis in children in Kazakhstan. We confirm that this work is original and has not been published elsewhere, nor is it currently under consideration for publication elsewhere. A preprint has previously been published [https://www.researchsquare.com/article/rs-3580377/v1] [DOI: 10.21203/rs.3.rs-3580377/v1] [65].

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors' Contributions

Dinara Kamalova conducted formal analysis and validation of the data, writing, reviewing, editing, and contribution to the final version of the manuscript. Akhmetova Assel performed conceptualization, data curation, writing of the original draft, writing, reviewing, editing, and contribution to the final version of the manuscript. Amirgazin Asylulan conducted conceptualization, formal analysis, validation of the data, data curation, writing, reviewing, editing, and contribution to the final version of the manuscript. Igor Sytnik conducted conceptualization, data curation, writing of the original draft, and contribution to the final version of the manuscript. Rudenko Viktoriya conducted data collection and processing of patient data and interpretation of results. Yessimkhanova Gulzhan conducted data collection and processing of patient data and interpretation of results. Bayesheva Dinagul conducted data collection and processing of patient data and interpretation of results. Yegorov Sergey conducted writing, reviewing, editing, and contribution to the final version of the manuscript. Alexander Shustov conducted writing of the original draft, formal analysis, validation of the data, data curation, writing, reviewing, editing, and contribution to the final version of the manuscript. Dinara Kamalova and Akhmetova Assel are the co-first authors.

Supplementary Materials

Supplementary Materials Figures S1–S8 represent maximum likelihood (ML) phylogenetic trees reconstructed for specific human enterovirus B serotypes CVA9, CVB5, CVB3, E6, E9, E11, E21, and E25, respectively, using the whole VP1 gene sequences of collected from available worldwide sequences (of the serogroup) and from infected patients in Kazakhstan. Within each serotype, the genotypes' classification is shown based on within (<15%) and between clades' diversity (>15%); bootstrap values >50% shown as blue circles. Tables S1-S2 demonstrate the estimates of evolutionary divergence between the sequences of serotype E21 and E25, respectively. The number of base substitutions per site from between sequences is shown. A lower triangle shows between clades' (genotype groups) divergence in nucleotides. Diagonals show within clades' (genotype groups) divergence in nucleotides within each serotype. Tables S3–S5 represent estimates of pairwise evolutionary divergence between sequences within serotypes CVA9, CVB5, and E6, respectively. The number of base substitutions per site between sequences is shown.

Figure 1 Flowchart showing the sampling scheme of the EV strains collected from infected patients in May-October 2022.

Figure 2 (a, b). Distribution of samples collected from patients for EV diagnosis: (a) by sampling time; (b) by age groups. X-axis represents number of samples collected from patients from total number of samples (n = 152). Y-axis shows (a) months; (b) age groups. Numbers on histogram bars show exact number of samples collected in female (orange)/male (blue) patients.

Figure 3 (a) Phylogenetic dendrograms were reconstructed by using the whole VP1 gene sequence of 65 human enterovirus B sequences collected from infected patients in Kazakhstan and 81 sequences downloaded from available worldwide human enterovirus B strains of serotypes CVA9, CVB3, and CVB5. (b) Phylogenetic dendrograms reconstructed by using the whole VP1 gene sequence of 65 human enterovirus B sequences collected from infected patients in Kazakhstan and 81 sequences downloaded from available worldwide human enterovirus B strains of serotypes E6, E9, E11, E21, and E25.

Figure 4 EV serotypes identified in patients with diagnosed aseptic meningitis disease and herpangina by sample types.

Figure 5 Phylogenetic analysis of the 65 complete coding sequences obtained in the current study.

Table 1 Enterovirus results obtained in current study with metadata from patients.

#	Diagnosis new	Age	Sex	Date_hosp	Date_sample	Headache	Fever	Vomit	Weakness	Rush	Severity	Sample type	Sample ID	Serotype	
1	Aseptic meningitis	5.0	F	9/11/2022	11.09.2022	Y	38	Y	Y	N	N	CSF	EV-105	CVA9	
2	Aseptic meningitis	13.0	M	6/3/2022	03.06.2022	Y	38.5	Y	Y	N	N	CSF	EV-12	CVB5	
3	Aseptic meningitis	5.0	M	6/1/2022	02.06.2022	Y	39	Y	Y	N	N	CSF	EV-14	E11	
4	Aseptic meningitis	9.0	M	6/14/2022	14.06.2022	Y	37.5	Y	Y	N	N	CSF	EV-16	Echo21	
5	Aseptic meningitis	7.0	M	6/14/2022	14.06.2022	Y	37.2	Y	Y	N	N	CSF1	EV-17 (vir1-vir2)	E11	
CSF2	Echo21	
6	Aseptic meningitis	13.0	F	6/10/2022	12.06.2022	Y	38	Y	Y	N	N	CSF	EV-18	CVB5	
7	Aseptic meningitis	5.0	F	6/12/2022	01.06.2022	Y	 	N	Y	N	N	CSF	EV-2	E9	
8	Aseptic meningitis	10.0	M	6/9/2022	09.06.2022	Y	38	Y	Y	N	N	CSF1	EV-22 (vir1-vir2)	E11	
CSF2	CVB5	
9	Aseptic meningitis	6.0	M	6/8/2022	08.06.2022	Y	38.2	Y	Y	N	N	CSF	EV-25	CVA9	
10	Aseptic meningitis	8.0	M	6/15/2022	15.06.2022	Y	37.6	Y	Y	N	N	CSF	EV-28	E9	
11	Aseptic meningitis	8.0	M	6/16/2022	16.06.2022	Y	39.9	Y	Y	N	N	CSF	EV-31	CVA9	
12	Aseptic meningitis	8.0	M	6/15/2022	15.06.2022	Y	37	Y	Y	N	N	CSF	EV-32	E9	
13	Aseptic meningitis	13.0	M	6/23/2022	23.06.2022	Y	37.5	Y	Y	N	N	CSF	EV-38	CVA9	
14	Aseptic meningitis	4.0	F	6/16/2022	17.06.2022	Y	38.5	Y	Y	N	N	CSF	EV-39	E9	
15	Aseptic meningitis	8.0	F	6/20/2022	20.06.2022	Y	37.8	Y	Y	N	N	CSF	EV-41	CVA9	
16	Aseptic meningitis	6.0	M	6/20/2022	20.06.2022	Y	38.9	Y	Y	Y	N	CSF	EV-42	CVA9	
17	Aseptic meningitis	9.0	F	6/27/2022	27.06.2022	Y	39	Y	Y	N	N	CSF	EV-45	CVA9	
18	Aseptic meningitis	5.0	M	6/27/2022	27.06.2022	Y	38	Y	Y	N	N	CSF	EV-49	CVA9	
19	Aseptic meningitis	12.0	M	6/29/2022	29.06.2022	Y	39	Y	Y	N	N	CSF	EV-50	CVB5	
20	Aseptic meningitis	13.0	F	6/29/2022	29.06.2022	Y	38.4	Y	Y	N	N	CSF	EV-53	CVA9	
21	Aseptic meningitis	8.0	M	5/29/2022	29.05.2022	Y	37.8	Y	Y	N	N	CSF	EV-6	CVA9	
22	Aseptic meningitis	9.0	M	7/2/2022	02.07.2022	Y	39	Y	Y	N	N	CSF	EV-64	CVB5	
23	Aseptic meningitis	7.0	F	7/10/2022	10.07.2022	Y	38.6	Y	Y	N	N	CSF	EV-69	CVB5	
24	Aseptic meningitis	17.0	M	5/31/2022	31.05.2022	Y	39.1	Y	Y	N	N	CSF	EV-7	CVA9	
25	Aseptic meningitis	12π	M	9/6/2022	06.09.2022	Y	38.3	Y	Y	N	N	CSF	EV-A1-2	CVB5	
26	Aseptic meningitis	9π	F	6/24/2022	24.06.2022	Y	37.8	Y	Y	N	N	SWAB	EV-A2	E9	
27	Aseptic meningitis	12π	M	9/18/2022	18.09.2022	Y	38	Y	Y	N	N	SWAB	EV-A2-2	CVB5	
28	Aseptic meningitis	9.0	M	6/29/2022	28.06.2022	Y	37.2	Y	Y	N	N	SWAB	EV-A3	E25	
29	Herpangina	6.0	F	7/5/2022	06.07.2022	N	38.5	Y	Y	Y	N	SWAB	EV-A4	CVA9	
30	Aseptic meningitis	3.5	M	6/21/2022	23.06.2022	Y	38	Y	Y	N	N	SWAB	EV-B1	CVA9	
31	Aseptic meningitis	4.0	F	9/10/20222	11.09.2022	Y	39	Y	Y	N	N	CSF	EV-B-2	CVB5	
32	Aseptic meningitis	8.0	F	9/8/20222	08.09.2022	Y	38.5	Y	Y	N	N	SWAB	EV-B2-2	CVB3	
33	Aseptic meningitis	9.0	M	6/28/20222	28.06.2022	Y	38.5	Y	Y	N	N	SWAB	EV-B3 (vir1-vir2)	E6	
SWAB	CVA9	
34	Aseptic meningitis	10.0	M	7/9/20222	09.07.2022	Y	37.5	Y	Y	N	N	SWAB	EV-B4	E25	
35	Aseptic meningitis	4.0	M	6/20/20222	20.06.2022	Y	37.8	Y	Y	N	N	SWAB	EV-C1	CVB5	
36	Aseptic meningitis	7.0	M	7/14/20222	17/14/20222	Y	39	Y	Y	N	N	SWAB	EV-C1-2	E6	
37	Herpangina	1.3	F	9/12/20222	13.09.2022	N	38.3	N	Y	Y	N	SWAB	EV-C2-2	CVA9	
38	Aseptic meningitis	8.0	M	6/27/20222	27.06.2022	Y	38	Y	Y	N	N	SWAB	EV-C3	CVA9	
39	Aseptic meningitis	5.0	F	7/17/20222	11.07.2022	Y	38.8	Y	Y	N	N	SWAB	EV-C4	E6	
40	Herpangina	1.9	M	6/21/20222	22.06.2022	N	39	N	Y	N	N	SWAB	EV-D1	CVA9	
41	Aseptic meningitis	5.0	F	7/15/20222	15.07.2022	Y	 	N	Y	N	N	SWAB	EV-D1-2	Echo21	
42	Aseptic meningitis	4.0	M	6/23/20222	23.06.2022	Y	38.8	Y	Y	N	N	SWAB	EV-D2	E25	
43	Herpangina	3.6	F	6/28/20222	29.06.2022	N	39	N	Y	N	N	SWAB	EV-D3	CVA9	
44	Herpangina	6.0	M	7/12/20222	12.07.2022	Y	39	Y	Y	N	N	SWAB	EV-D4	CVA9	
45	Aseptic meningitis	14.0	M	6/23/20222	23.06.2022	Y	37.5	Y	Y	N	N	SWAB	EV-E1	E9	
46	Aseptic meningitis	9.0	M	9/23/20222	24.09.2022	Y	37.3	Y	Y	N	N	CSF	EV-E1-2	CVB5	
47	Herpangina	2.3	M	6/23/20222	29.06.2022	N	39	N	Y	N	N	SWAB1	EV-E2 (vir1-vir2)	CVB5	
 	Herpangina	2.3	M	6/23/20222	29.06.2022	N	39	N	Y	N	N	SWAB2	 	E6	
48	Aseptic meningitis	6.0	M	6/27/20222	27.06.2022	Y	38	Y	Y	N	N	SWAB	EV-E3	E6	
49	Aseptic meningitis	13.0	M	7/13/20222	13.07.2022	Y	38	Y	Y	N	N	SWAB	EV-E4	E6	
50	Aseptic meningitis	7.0	M	6/22/20222	23.06.2022	Y	38	Y	Y	N	N	SWAB	EV-F1	E6	
51	Aseptic meningitis	5.0	M	9/9/20222	09.09.2022	Y	37.8	Y	Y	N	N	CSF	EV-F1-2	CVA9	
52	Aseptic meningitis	2.4	M	6/23/20222	24.06.2022	N	39.6	N	Y	Y	N	SWAB	EV-F2	E9	
53	Herpangina	15.0	M	6/28/20222	29.06.2022	N	38.5	Y	Y	N	N	SWAB	EV-F3	CVA9	
54	Herpangina	1.2	F	6/22/20222	23.06.2022	Y	39.2	Y	Y	N	N	SWAB	EV-G1	CVB5	
55	Aseptic meningitis	4.0	F	9/9/20222	09.09.2022	Y	39	Y	Y	N	N	CSF	EV-G1-2	E6	
56	Aseptic meningitis	5.0	F	6/29/20222	29.06.2022	Y	40	Y	Y	N	N	SWAB1	EV-G3 (vir1-vir2)	E6	
SWAB2	CVA9	
57	Aseptic meningitis	13.0	M	6/23/20222	23.06.2022	Y	37.5	Y	Y	N	N	SWAB1	EV-H1 (vir1-vir2)/E	E9	
SWAB2	CVA9	
CSF	CVA9	
58	Aseptic meningitis	12.0	F	7/4/20222	05.07.2022	Y	39	Y	Y	N	N	SWAB	EV-H3	CVA9
==== Refs
1 Park K. Lee B. Baek K. Enteroviruses isolated from herpangina and hand-foot-and-mouth disease in Korean children Virology Journal 2012 9 p. 205
2 Sawyer M. H. Enterovirus infections: diagnosis and treatment Seminars in Pediatric Infectious Diseases 2002 13 1 40 47 12118843
3 Corsino C. B. Ali R. Linklater D. R. Herpangina StatPearls 2023 Treasure Island, FL, USA StatPearls Publishing
4 Romero J. R. Hand, foot, and mouth disease and Herpangina 2022 https://www.uptodate.com/contents/hand-foot-and-mouth-disease-and-herpangina
5 Hyypia T. Hovi T. Knowles N. J. Stanway G. Classification of enteroviruses based on molecular and biological properties Journal of General Virology 1997 78 Pt 1 1 11 9010279
6 de Crom S. C. Rossen J. W. van Furth A. M. Obihara C. C. Enterovirus and parechovirus infection in children: a brief overview European Journal of Pediatrics 2016 175 8 1023 1029 27156106
7 Pöyry T. Kinnunen L. Hyypiä T. Genetic and phylogenetic clustering of enteroviruses Journal of General Virology 1996 77 Pt 8 1699 1717 8760417
8 Nasri D. Bouslama L. Omar S. Typing of human enterovirus by partial sequencing of VP2 Journal of Clinical Microbiology 2007 45 8 2370 2379 17537940
9 Laitinen O. H. Svedin E. Kapell S. Nurminen A. Hytönen V. P. Flodström-Tullberg M. Enteroviral proteases: structure, host interactions and pathogenicity Reviews in Medical Virology 2016 26 4 251 267 27145174
10 Joffret M. L. Polston P. M. Razafindratsimandresy R. Bessaud M. Heraud J. M. Delpeyroux F. Whole genome sequencing of enteroviruses species A to D by high-throughput sequencing: application for viral mixtures Frontiers in Microbiology 2018 9 p. 2339
11 Caro V. Guillot S. Delpeyroux F. Crainic R. Molecular strategy for “serotyping” of human enteroviruses Journal of General Virology 2001 82 Pt 1 79 91 11125161
12 Oberste M. S. Maher K. Kilpatrick D. R. Flemister M. R. Brown B. A. Pallansch M. A. Typing of human enteroviruses by partial sequencing of VP1 Journal of Clinical Microbiology 1999 37 5 1288 1293 10203472
13 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 Journal of Virology 1999 73 3 1941 1948 9971773
14 King A. M. Q. Family picornaviridae Virus Taxonomy, Seventh Report of the International Committee on Taxonomy of viruses 2000
15 Bubba L. Martinelli M. Pellegrinelli L. A 4-year study on epidemiologic and molecular characteristics of human parechoviruses and enteroviruses circulating in children younger than 5 Years in northern Italy The Pediatric Infectious Disease Journal 2017 36 1 13 19 10.1097/inf.0000000000001344 2-s2.0-84991494489 27749661
16 Stellrecht K. A. Harding I. Woron A. M. Lepow M. L. Venezia R. A. The impact of an enteroviral RT-PCR assay on the diagnosis of aseptic meningitis and patient management Journal of Clinical Virology 2002 25 Suppl 1 S19 S26 12091078
17 Majer A. McGreevy A. Booth T. F. Molecular pathogenicity of enteroviruses causing neurological disease Frontiers in Microbiology 2020 11 p. 540
18 Isaacs S. R. Kim K. W. Cheng J. X. Amplification and next generation sequencing of near full-length human enteroviruses for identification and characterisation from clinical samples Scientific Reports 2018 8 1 11889
19 Andrews S. FastQC: a quality control tool for high throughput sequence data 2010 http://www.bioinformatics.babraham.ac.uk/projects/fastqc/
20 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 Journal of Virology 2004 78 2 855 867 14694117
21 Darriba D. Posada D. Kozlov A. M. Stamatakis A. Morel B. Flouri T. ModelTest-NG: a new and scalable tool for the selection of DNA and protein evolutionary models Molecular Biology and Evolution 2020 37 1 291 294 31432070
22 Letunic I. Bork P. Interactive Tree of Life (iTOL): an online tool for phylogenetic tree display and annotation Bioinformatics 2007 23 1 127 128 17050570
23 Zhao H. Wang J. Chen J. Molecular epidemiology and evolution of coxsackievirus A9 Viruses 2022 14 4
24 He Y. Wei H. Wei L. Molecular epidemiology reveals the Co-circulation of two genotypes of coxsackievirus B5 in China Viruses 2022 14 12
25 Yang Q. Yan D. Song Y. Whole-genome analysis of coxsackievirus B3 reflects its genetic diversity in China and worldwide Virology Journal 2022 19 1 p. 69
26 Cheng W. Ji T. Zhou S. Molecular epidemiological characteristics of echovirus 6 in mainland China: extensive circulation of genotype F from 2007 to 2018 Archives of Virology 2021 166 5 1305 1312 33638089
27 Zhang M. Guo W. Xu D. Molecular characterization of echovirus 9 strains isolated from hand-foot-and-mouth disease in Kunming, Yunnan Province, China Scientific Reports 2022 12 1 p. 2293
28 Li J. Yan D. Chen L. Multiple genotypes of Echovirus 11 circulated in mainland China between 1994 and 2017 Scientific Reports 2019 9 1 10583
29 Grapin M. Mirand A. Pinquier D. Severe and fatal neonatal infections linked to a new variant of echovirus 11, France, July 2022 to April 2023 Euro Surveillance 2023 28 22
30 Kumar S. Stecher G. Li M. Knyaz C. Tamura K. Mega X: molecular evolutionary genetics analysis across computing platforms Molecular Biology and Evolution 2018 35 6 1547 1549 29722887
31 Katoh K. Standley D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability Molecular Biology and Evolution 2013 30 4 772 780 23329690
32 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 Evolution 2016 2 1 p. vew007
33 R Core Team R: A Language and Environment for Statistical Computing 2022 Vienna, Austria R Foundation for Statistical Computing https://www.R-project.org
34 Oberste M. S. Maher K. Kennett M. L. Molecular epidemiology and genetic diversity of echovirus type 30 (E30): genotypes correlate with temporal dynamics of E30 isolation Journal of Clinical Microbiology 1999 37 12 3928 3933 10565909
35 Piralla A. Borghesi A. Di Comite A. Fulminant echovirus 11 hepatitis in male non-identical twins in northern Italy, April 2023 Euro Surveillance 2023 28 24
36 World Health Organization Enterovirus Infection-France 2023 Geneva, Switzerland World Health Organization https://www.who.int/emergencies/disease-outbreak-news/item/2023-DON469
37 Begaidarova R. K. Devdariani K. G. Baygutanova G. Z. Starikov Y. G. Dzhshibekova Z. Z. Zolotareva O. A. Clinical and laboratory features of enteroviral meningitis caused by echo-30 in children of different age groups Medicine and Ecology 2015 2 75 47 54
38 Bumburidi Y. Utepbergenova G. Yerezhepov B. Etiology of acute meningitis and encephalitis from hospital-based surveillance in South Kazakhstan oblast, February 2017-January 2018 PLoS One 2021 16 5 e0251494
39 Brouwer L. Moreni G. Wolthers K. C. Pajkrt D. World-wide prevalence and genotype distribution of enteroviruses Viruses 2021 13 3 10.3390/v13030434
40 Guo W. Xu D. Cong S. Co-infection and enterovirus B: post EV-A71 mass vaccination scenario in China BMC Infectious Diseases 2022 22 1 p. 671
41 Ooi M. H. Wong S. C. Podin Y. Human enterovirus 71 disease in Sarawak, Malaysia: a prospective clinical, virological, and molecular epidemiological study Clinical Infectious Diseases 2007 44 5 646 656 17278054
42 Yang F. Du J. Hu Y. Enterovirus coinfection during an outbreak of hand, foot, and mouth disease in Shandong, China Clinical Infectious Diseases 2011 53 4 400 401 21785005
43 Batcup G. Holt P. Hambling M. H. Gerlis L. M. Glass M. R. Placental and fetal pathology in Coxsackie virus A9 infection: a case report Histopathology 1985 9 11 1227 1235 10.1111/j.1365-2559.1985.tb02802.x 2-s2.0-0022378509 4085986
44 Huang Y. C. Chu Y. H. Yen T. Y. Clinical features and phylogenetic analysis of coxsackievirus A9 in northern Taiwan in 2011 BMC Infectious Diseases 2013 13 p. 33
45 Johnson J. E. Vogt P. K. Cell entry by non-enveloped viruses Current Topics in Microbiology and Immunology 2010 343
46 Shabani A. Makvandi M. Samarbafzadeh A. Echovirus 30 and coxsackievirus A9 infection among young neonates with sepsis in Iran Iranian Journal of Microbiology 2018 10 4 258 265 30483379
47 Cui A. Yu D. Zhu Z. An outbreak of aseptic meningitis caused by coxsackievirus A9 in Gansu, the People’s Republic of China Virology Journal 2010 7 1 p. 72
48 Pabbaraju K. Wong S. Chan E. N. Tellier R. Genetic characterization of a Coxsackie A9 virus associated with aseptic meningitis in Alberta, Canada in 2010 Virology Journal 2013 10 p. 93
49 Smuts H. Cronje S. Thomas J. Brink D. Korsman S. Hardie D. Molecular characterization of an outbreak of enterovirus-associated meningitis in mossel Bay, South Africa, December 2015–January 2016 BMC Infectious Diseases 2018 18 1 p. 709
50 Itani T. M. Chalapa V. I. Slautin V. N. Non-polio enterovirus surveillance in the ural federal district and western Siberia, 2022: is there a need for a vaccine? Vaccines 2023 11 10 p. 1588
51 Novikova N. A. Golitsina L. N. Fomina S. G. Efimov E. I. Molecular monitoring of non-polio enteroviruses in European territory of Russia in 2008–2011 Zhurnal Mikrobiologii, Epidemiologii i Immunobiologii 2013 1 75 78 23805657
52 Suresh S. Rawlinson W. D. Andrews P. I. Stelzer-Braid S. Global epidemiology of nonpolio enteroviruses causing severe neurological complications: a systematic review and meta-analysis Reviews in Medical Virology 2020 30 1 e2082
53 Toczylowski K. Wieczorek M. Bojkiewicz E. Wietlicka-Piszcz M. Gad B. Sulik A. Pediatric enteroviral central nervous system infections in Bialystok, Poland: epidemiology, viral types, and drivers of seasonal variation Viruses 2020 12 8
54 Liu H. Zhang H. Zhang M. Epidemiological and etiological characteristics of viral meningitis for hospitalized pediatric patients in Yunnan, China Medicine (Baltimore) 2022 101 26 e29772
55 Wang H. Y. Li Y. Xu A. Q. Identification and phylogenic analysis of Coxsackie-virus B5 during an outbreak of aseptic meningitis in Shandong Zhonghua Liuxingbingxue Zazhi 2010 31 1 64 68 20302702
56 Farshadpour F. Taherkhani R. Molecular epidemiology of enteroviruses and predominance of echovirus 30 in an Iranian population with aseptic meningitis Journal of NeuroVirology 2021 27 3 444 451 33788142
57 Centers for Disease Control and Prevention Nonpolio enterovirus and human parechovirus surveillance—United States, 2006–2008 Morbidity & Mortality Weekly Report 2010 59 48 1577 1580 21150865
58 Khetsuriani N. Kutateladze T. Zangaladze E. High degree of genetic diversity of non-polio enteroviruses identified in Georgia by environmental and clinical surveillance, 2002–2005 Journal of Medical Microbiology 2010 59 Pt 11 1340 1347 20671086
59 Cabrerizo M. Trallero G. Simmonds P. Recombination and evolutionary dynamics of human echovirus 6 Journal of Medical Virology 2014 86 5 857 864 24114692
60 Wieczorek M. Figas A. Krzysztoszek A. Enteroviruses associated with aseptic meningitis in Poland, 2011–2014 Polish Journal of Microbiology 2016 65 2 231 235
61 Tao Z. Wang H. Li Y. Liu G. Xu A. Lin X. Molecular epidemiology of human enterovirus associated with aseptic meningitis in Shandong Province, China, 2006–2012 PLoS One 2014 9 2
62 Wang J. Meng M. Xu H. Analysis of enterovirus genotypes in the cerebrospinal fluid of children associated with aseptic meningitis in Liaocheng, China, from 2018 to 2019 BMC Infectious Diseases 2021 21 1 p. 405 33933008
63 Zhu Y. Zhou X. Liu J. Molecular identification of human enteroviruses associated with aseptic meningitis in Yunnan province, Southwest China SpringerPlus 2016 5 1 p. 1515
64 Ramalho E. Sousa I. Jr. Burlandy F. Identification and phylogenetic characterization of human enteroviruses isolated from cases of aseptic meningitis in Brazil, 2013–2017 Viruses 2019 11 8 p. 690
65 Kamalova D. Akhmetova A. Amirgazin A. Sytnik I. Shevtsov A. Genetic diversity of enteroviral meningitis, herpangina and hand, foot, and mouth disease in Kazakhstan Preprint (Version 1) Available at Research Square 2023
