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

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

Pathological Characteristics of Echovirus 30 Infection in a Mouse Model.

Echovirus 30 (E30), a member of species B enterovirus, is associated with outbreaks of aseptic meningitis and has become a global health emergency. However, the pathogenesis of E30 remains poorly understood due to the lack of appropriate animal models. In this study, we established a mouse infection model to explore the pathogenicity of E30. The 2-day-old IFNAR-/- mice infected with E30 strain WZ16 showed lethargy and paralysis, and some died. Obvious pathological changes were observed in the skeletal muscle, brain tissue, and other tissues, with the highest viral load in the skeletal muscles. Transcriptome analysis of brain and skeletal muscle tissues from infected mice showed that significant differentially expressed genes were enriched in complement response and neuropathy-related pathways. Using immunofluorescence assay, we found that the viral double-stranded RNA (dsRNA) was detected in the mouse brain region and could infect human glioma (U251) cells. These results indicated that E30 affects the nervous system, and they provide a theoretical basis for understanding its pathogenesis. IMPORTANCE Echovirus 30 (E30) infection causes a wide spectrum of diseases with mild symptoms, such as hand, foot, and mouth disease (HFMD), acute flaccid paralysis, and aseptic meningitis and other diseases, especially one of the most common pathogens causing aseptic meningitis outbreaks. We established a novel mouse model of E30 infection by inoculating neonatal mice with clinical isolates of E30 and observed the pathological changes induced by E30. Using the E30 infection model, we found complement responses and neuropathy-related genes in the mice tissues at the transcriptome level. Moreover, we found that the viral dsRNA localized in the mouse brain and could replicate in human glioma cell line U251 rather than in the neuroblastoma cell line, SK-N-SH.

Whole-genome analysis of coxsackievirus B3 reflects its genetic diversity in China and worldwide.

BACKGROUND: Coxsackievirus B3 (CVB3) has emerged as an active pathogen in myocarditis, aseptic meningitis, hand, foot, and mouth disease (HFMD), and pancreatitis, and is a heavy burden on public health. However, CVB3 has not been systematically analyzed with regard to whole-genome diversity and recombination. Therefore, this study was undertaken to systematically examine the genetic characteristics of CVB3 based on its whole genome. METHODS: We combined CVB3 isolates from our national HFMD surveillance and global sequences retrieved from GenBank. Phylogenetic analysis was performed to examine the whole genome variety and recombination forms of CVB3 in China and worldwide. RESULTS: Phylogenetic analysis showed that CVB3 strains isolated worldwide could be classified into clusters A-E based on the sequence of the entire VP1 region. The predominant CVB3 strains in China belonged to cluster D, whereas cluster E CVB3 might be circulated globally compared to other clusters. The average nucleotide substitution rate in the P1 region of CVB3 was 4.82 × 10-3 substitutions/site/year. Myocarditis was more common with cluster A. Clusters C and D presented more cases of acute flaccid paralysis, and cluster D may be more likely to cause HFMD. Multiple recombination events were detected among CVB3 variants, and there were twenty-three recombinant lineages of CVB3 circulating worldwide. CONCLUSIONS: Overall, this study provides full-length genomic sequences of CVB3 isolates with a wide geographic distribution over a long-term time scale in China, which will be helpful for understanding the evolution of this pathogen. Simultaneously, continuous surveillance of CVB3 is indispensable to determine its genetic diversity in China as well as worldwide.

Temporal phylogeny and molecular characterization of echovirus 30 associated with aseptic meningitis outbreaks in China.

BACKGROUND: An outbreak of aseptic meningitis occurred from June to August 2016, in Inner Mongolia Autonomous Region, China. METHODS: To determine its epidemiological characteristics, etiologic agent, and possible origin, specimens were collected for virus isolation and identification, followed by molecular epidemiological analysis. RESULTS: A total of 363 patients were clinically diagnosed from June 1st to August 31st 2016, and most cases (63.1%, n = 229) were identified between June 22nd and July 17th, with children aged 6 to 12 years constituting the highest percentage (68.9%, n = 250). All viral isolates from this study belonged to genotype C of echovirus 30 (E30), which dominated transmission in China. To date, two E30 transmission lineages have been identified in China, of which Lineage 2 was predominant. We observed fluctuant progress of E30 genetic diversity, with Lineage 2 contributing to increased genetic diversity after 2002, whereas Lineage 1 was significant for the genetic diversity of E30 before 2002. CONCLUSIONS: We identified the epidemiological and etiological causes of an aseptic meningitis outbreak in Inner Mongolia in 2016, and found that Lineage 2 played an important role in recent outbreaks. Moreover, we found that Gansu province could play an important role in E30 spread and might be a possible origin site. Furthermore, Fujian, Shandong, Taiwan, and Zhejiang provinces also demonstrated significant involvement in E30 evolution and persistence over time in China.

Development of a real-time RT-PCR assay for the detection of pan-human parechoviruses.

BACKGROUND: Parechoviruses (PeV-As), which constitute a new genus within the family Picornaviridae, have been associated with numerous localized outbreaks of serious diseases, such as coryza, pneumonia, maculopapular exanthem, and conjunctivitis. However, to the best of our knowledge, only a few laboratories worldwide conduct tests for the identification of this group of viruses. Therefore, in this study, we aimed to develop and validate a real-time RT-PCR assay for the identification of PeV-As. METHODS: To design and validate a real-time PCR primer-probe targeting the 5'-UTR region of PeV-As, the 5'-UTR sequences of PeV-As available in GenBank were aligned using the MUSCLE algorithm in MEGA v7.0. Thereafter, the highly conserved 5'-UTR region was selected, and its primer-probe sequence was designed using Primer Premier v5.0. This primer-probe sequence was then evaluated for specificity, sensitivity, and repeatability, and for its validation, it was tested using fecal samples from 728 healthy children living in Beijing (China). RESULTS: The PeV-A real-time RT-PCR assay detected only the RNA-positive standards of PeV-A genotypes (1-8, 14, 17, and 18), whereas 72 serotypes of non-PeV-A EV viruses were undetected. In addition, the VP1 region of these 11 PeV-A genotypes that tested positive were amplified using the primers designed in this study. Typing results indicated that eight, one, and two strains of the 11 were PeV-A1, PeV-A4, and PeV-A6, respectively. We also determined and presented the genetic characterization and phylogenetic analyses results corresponding to these 11 VP1 region sequences. Furthermore, real-time RT-PCR assay showed good sensitivity with LOD of 102 copies/µL. Positive results in eight parallel experiments at each concentration gradient from 107 copies/µL to 102 copies/µL, indicating good repeatability. CONCLUSION: Our findings suggested that the real-time RT-PCR assay developed in this study can be applied for routine PeV-A identification. We detected PeV-A1, 4 and 6 genotypes in the 728 faecal samples using this method. Additionally, we believe that our results will serve as a foundation for further studies on PeV-As and facilitate the expansion of the gene sequence information available in GenBank.

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

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

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

BACKGROUND: Coxsackievirus B3 (CV-B3) is usually associated with aseptic meningitis and myocarditis; however, the association between CV-B3 and hand, foot, and mouth disease (HFMD) has not been clearly demonstrated, and the phylogenetic dynamics and transmission history of CV-B3 have not been well summarized. METHOD: Two HFMD outbreaks caused by CV-B3 were described in Hebei Province in 2012 and in Shandong Province in 2016 in China. To analyze the epidemiological features of two CV-B3 outbreaks, a retrospective analysis was conducted. All clinical specimens from CV-B3 outbreaks were collected and disposed according to the standard procedures supported by the WHO Global Poliovirus Specialized Laboratory. EV genotyping and phylogenetic analysis were performed to illustrate the genetic characteristics of CV-B3 in China and worldwide. RESULTS: Two transmissible lineages (lineage 2 and 3) were observed in Northern China, which acted as an important "reservoir" for the transmission of CV-B3. Sporadic exporting and importing of cases were observed in almost all regions. In addition, the global sequences of CV-B3 showed a tendency of geographic-specific clustering, indicating that geographic-driven adaptation plays a major role in the diversification and evolution of CV-B3. CONCLUSIONS: Overall, our study indicated that CV-B3 is a causative agent of HFMD outbreak and revealed the phylogenetic dynamics of CV-B3 worldwide, as well as provided an insight on CV-B3 outbreaks for effective intervention and countermeasures.

Synergetic association between coxsackievirus A16 genotype evolution and recombinant form shifts.

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

Diverse Head-to-Tail Sequences in the Circular Genome of Human Bocavirus Genotype 1 among Children with Acute Respiratory Infections Implied the Switch of Template Chain in the Rolling-Circle Replication Model.

Head-to-tail sequences have been reported in human bocavirus (HBoV) 1-4. To reveal their features and functions, HBoV DNA was screened among respiratory specimens from pediatric patients with an acute respiratory infection (ARI) between April 2020 and December 2022, followed by HBoV genotyping. Head-to-tail sequences were detected using nested PCR, TA cloning, and Sanger sequencing, and these findings were confirmed by mNGS and amplicon sequencing. The secondary structure was predicted using the Mfold web server. The results indicated that head-to-tail sequences were detected in 42 specimens through TA cloning from 351 specimens positive for HBoV1 DNA, yielding 92 sequences into 32 types and 2 categories. Additionally, head-to-tail sequences were detected in 16 specimens by amplicon sequencing, yielding 60 sequences categorized into 23 types. The 374nt type, detected in 13 specimens, contains variants 374a and 374b, which differ in the unpaired loop regions of the palindrome or complementary reverse sequences, implying a switch of template chains during the replication process. The mNGS results in three specimens confirmed the presence of circular genome in copies below 1%. In conclusion, head-to-tail sequences of HBoV1 were common in children with ARI and were highly diverse in length and sequences. The variants may be generated by the switch of the template chain in the rolling-circle replication model.

The Cumulative Variations of Respiratory Syncytial Virus Fusion Protein (F) in Ten Consecutive Years in China.

BACKGROUND: Variations in the fusion (F) protein of respiratory syncytial virus (RSV) with main antigenic sites I-V and Ø may affect the development of RSV vaccines and therapies. METHODS: In the study, 30 respiratory specimens positive for RSV were randomly selected from children with acute lower respiratory infections (ALRI) in Beijing every year from 2012 to 2021 for F gene sequencing. Then, 300 F gene sequences and 508 uploaded to GenBank from China were subjected to phylogenetic analysis. RESULTS: The results indicated the nucleotide identities were 95.4-100% among 446 sequences of RSV A, and 96.3-100% among 362 of RSV B. The most common variant loci were N80K (100.00%) and R213S (97.76%) for site Ø, and V384I/T (98.43%) for site I among sequences of RSV A, and M152I (100.00%), I185V (100.00%), and L172Q/H (94.48%) for site V, and R202Q (99.45%) for site Ø among sequences of RSV B. N276S appears in 95.29% sequences of RSV A, while S276N and N262 I/S appear in 1.38% and 0.55% sequences of RSV B, respectively. No variation was found in all sequences at the binding sites of 14N4 and motavizumab. CONCLUSIONS: There were cumulative variations of the RSV F gene, especially at some binding sites of antigenic sites.

High-Frequency Recombination of Human Adenovirus in Children with Acute Respiratory Tract Infections in Beijing, China.

Recombination events in human adenovirus (HAdV) have led to some new highly pathogenic or infectious types. It is vital to monitor recombinant HAdVs, especially in children with acute respiratory tract infections (ARIs). In the retrospective study, HAdV positive specimens were collected from pediatric patients with ARIs during 2015 to 2021, then typed by sequence analysis of the penton base, hexon and fiber gene sequence. For those with inconsistent typing results, a modified method with species-specific primer sets of a fiber gene sequence was developed to distinguish co-infections of different types from recombinant HAdV infections. Then, plaque assays combined with meta-genomic next-generation sequencing (mNGS) were used to reveal the HAdV genomic characteristics. There were 466 cases positive for HAdV DNA (2.89%, 466/16,097) and 350 (75.11%, 350/466) successfully typed with the most prevalent types HAdV-B3 (56.57%, 198/350) and HAdV-B7 (32.00%, 112/350), followed by HAdV-C1 (6.00%, 21/350). Among 35 cases (7.51%, 35/466) with inconsistent typing results, nine cases were confirmed as co-infections by different types of HAdVs, and 26 cases as recombinant HAdVs in six genetic patterns primarily clustered to species C (25 cases) in pattern 1-5, or species D (1 case) in pattern 6. The novel recombinant HAdV of species D was identified with multiple recombinant events among HAdV-D53, HAdV-D64, and HAdV-D8, and officially named as HAdV-D115. High-frequency recombination of HAdVs in six genetic recombination patterns were identified among children with ARIs in Beijing. Specifically, there is a novel Adenovirus D human/CHN/S8130/2023/115[P22H8F8] designed as HAdV D115.

Evolution, recombination and geographic spreading of global Coxsackievirus A6.

BACKGROUND: The increasing incidence of hand, foot, and mouth disease (HFMD) associated with Coxsackievirus A6 (CVA6) has become a very significant public health problem. The aim of this study is to investigate the recombination, geographic transmission, and evolutionary characteristics of the global CVA6. METHODS: From 2019 to 2022, 73 full-length CVA6 sequences were obtained from HFMD patients in China and analyzed in combination with 1032 published whole genome sequences. Based on this dataset, the phylogenetic features, recombinant diversity, Bayesian phylodynamic characteristics, and key amino acid variations in CVA6 were analyzed. RESULTS: The four genotypes of CVA6, A, D, E, and F, are divided into 24 recombinant forms (RFs, RF-A - RF-X) based on differences in the P3 coding region. The eastern China region plays a key role in the dissemination of CVA6 in China. VP1-137 and VP1-138 are located in the DE loop on the surface of the CVA6 VP1 protein, with the former being a highly variable site and the latter having more non-synonymous substitutions. CONCLUSIONS: Based on whole genome sequences, this study contributes to the CVA6 monitoring, early warning, and the pathogenic mechanism by studying recombination diversity, geographical transmission characteristics, and the variation of important amino acid sites.

Genetic Evolution and Variation of Human Adenovirus Serotype 31 Epidemic Strains in Beijing, China, during 2010-2022.

Human adenovirus serotype 31 (HAdV-31) is closely associated with gastroenteritis in children and can cause fatal systemic disseminated diseases in immunocompromised patients. The lack of genomic data for HAdV-31, especially in China, will greatly limit research on its prevention and control. Sequencing and bioinformatics analyses were performed for HAdV-31 strains from diarrheal children in Beijing, China, during 2010-2022. Three capsid protein genes (hexon, penton, and fiber) were obtained in 37 cases, including one in which the whole genome was sequenced. HAdV-31 strains clustered into three distinct clades (I-III) in a phylogenetic tree constructed based on concatenated genes and the whole genome; the endemic strains only gathered into clade II, and most of the reference strains clustered into clade I. Compared with penton and hexon, fiber had a faster evolutionary rate (1.32 × 10-4 substitutions/site/year), an earlier divergence time (1697), lower homology (98.32-100% at the amino acid level), and greater genetic variation (0.0032). Four out of the six predicted positive selection pressure codons were also in the knob of fiber. These results reveal the molecular evolution characteristics and variations of HAdV-31 in Beijing, and fiber may be one of the main evolution driving forces.

Genomic Epidemiology and Transmission Dynamics of Global Coxsackievirus B4.

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

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

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

Molecular analysis of Coxsackievirus A24 variant isolates from three outbreaks of acute hemorrhagic conjunctivitis in 1988, 1994 and 2007 in Beijing, China.

Coxsackievirus A24 variant (CVA24v) is a major pathogen that causes continued outbreaks and pandemics of acute hemorrhagic conjunctivitis (AHC). In China, the first confirmed outbreak of CVA24v-related AHC occurred in Beijing in 1988, followed by another two significant outbreaks respectively in 1994 and 2007, which coincides with the three-stage dynamic distribution of AHC in the world after 1970s. To illustrate the genetic characteristics of CVA24v in different periods, a total of 23 strains were isolated from those three outbreaks and the whole genome of those isolations were sequenced and analyzed. Compared with the prototype strain, the 23 strains shared four nucleotide deletions in the 5' UTR except the 0744 strain isolated in 2007. And at the 98th site, one nucleotide insertion was found in all the strains collected from 2007. From 1994 to 2007, amino acid polarity in the VP1 region at the 25th and the 32nd site were changed. Both the 3C and VP1 phylogenetic tree indicated that isolates from 1988 and 1994 belonged to Genotype III (GIII), and 2007 strains to Genotype IV (GIV). According to the Bayesian analysis based on complete genome sequence, the most recent common ancestors for the isolates in 1988, 1994 and 2007 were respectively estimated around October 1987, February 1993 and December 2004. The evolutionary rate of the CVA24v was estimated to be 7.45 â€‹× â€‹10-3 substitutions/site/year. Our study indicated that the early epidemic of CVA24v in Chinese mainland was the GIII. Point mutations and amino acid changes in different genotypes of CVA24v may generate intensity differences of the AHC outbreak. CVA24v has been evolving constantly with a relatively rapid rate.

Molecular Characteristics and Genetic Evolution of Echovirus 33 in Mainland of China.

Echovirus, a member of the Enterovirus B (EV-B) family, has led to numerous outbreaks and pandemics, causing a broad spectrum of diseases. Based on the national hand, foot, and mouth disease (HFMD) surveillance system, seven strains of echovirus 33 (E33) were isolated from Mainland of China between 2010 and 2018. The whole genomes of these strains were isolated and sequenced, and phylogenetic trees were constructed based on the gene sequences in different regions of the EV-B prototype strains. It was found that E33 may be recombined in the P2 and P3 regions. Five genotypes (A-E) were defined based on the entire VP1 region of E33, of which the C gene subtype was the dominant gene subtype at present. Recombinant analysis showed that genotype C strains likely recombined with EV-B80, EV-B85, E13, and CVA9 in the P2 and P3 regions, while genotype E had the possibility of recombination with CVB3, E3, E6, and E4. Results of Bayesian analysis indicated that E33 may have appeared around 1955 (95% confidence interval: 1945-1959), with a high evolutionary rate of 1.11 × 10-2 substitution/site/year (95% highest posterior density (HPD): 8.17 × 10-3 to 1.4 × 10-2 substitution/site/year). According to spatial transmission route analysis, two significant transmission routes were identified: from Australia to India and from Oman to Thailand, which the E33 strain in Mainland of China likely introduced from Mexico and India. In conclusion, our study fills the gaps in the evolutionary analysis of E33 and can provide important data for enterovirus surveillance.

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

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

Genetic characterization and molecular epidemiology of Coxsackievirus A12 from mainland China during 2010-2019.

Coxsackievirus A12 (CVA12) is an enterovirus that has been isolated in many countries in recent years. However, studies on CVA12 are limited, and its effective population size, evolutionary dynamics and recombination patterns have not been clarified now. In this study, we described the phylogenetic characteristics of 16 CVA12 strains isolated from pediatric HFMD patients in mainland China from 2010 to 2019. Comparison of the nucleotide sequences and amino acid sequences with the CVA12 prototype strain revealed that the 16 CVA12 strains are identical in 78.8-79% and 94-94.2%, respectively. A phylodynamic analysis based on the 16 full-length VP1 sequences from this study and 21 sequences obtained from GenBank revealed a mean substitution rate of 6.61 × 10-3 substitutions/site/year (95% HPD: 5.16-8.20 × 10-3), dating the time to most recent common ancestor (tMRCA) of CVA12 back to 1946 (95% HPD: 1942-1947). The Bayesian skyline plot showed that the effective population size has experienced twice dynamic fluctuations since 2007. Phylogeographic analysis identified two significant migration pathways, indicating the existence of cross-provincial transmission of CVA12 in mainland China. Recombination analysis revealed two recombination patterns between 16 CVA12 strains and other EV-A, suggesting that there may be extensive genetic exchange between CVA12 and other enteroviruses. In summary, a total of 16 full-length CVA12 strains were reported in this study, providing valuable references for further studies of CVA12 worldwide.

Molecular Epidemiology and Evolution of Coxsackievirus A9.

Nineteen CVA9 isolates were obtained between 2010 and 2019 from six provinces of mainland China, using the HFMD surveillance network established in China. Nucleotide sequencing revealed that the full-length VP1 of 19 CVA9 isolates was 906 bases encoding 302 amino acids. The combination of the thresholds of the phylogenetic tree and nucleotide divergence of different genotypes within the same serotype led to a value of 15-25%, and enabled CVA9 worldwide to be categorized into ten genotypes: A-J. The phylogenetic tree showed that the prototype strain was included in genotype A, and that the B, C, D, E, H, and J genotypes disappeared during virus evolution, whereas the F, I, and G genotypes showed co-circulation. Lineage G was the dominant genotype of CVA9 and included most of the strains from nine countries in Asia, North America, Oceania, and Europe. Most Chinese strains belonged to the G genotype, suggesting that the molecular epidemiology of China is consistent with that observed worldwide. The 165 partial VP1 strains (723 nt) showed a mean substitution rate of 3.27 × 10-3 substitution/site/year (95% HPD range 2.93-3.6 × 10-3), dating the tMRCA of CVA9 back to approximately 1922 (1911-1932). The spatiotemporal dynamics of CVA9 showed the spread of CVA9 obviously increased in recent years. Most CVA9 isolates originated in USA, but the epidemic areas of CVA9 are now concentrated in the Asia-Pacific region, European countries, and North America. Recombination analysis within the enterovirus B specie (59 serotypes) revealed eight recombination patterns in China at present, CVB4, CVB5, E30, CVB2, E11, HEV106, HEV85, and HEV75. E14, and E6 may act as recombinant donors in multiple regions. Comparison of temperature sensitivity revealed that temperature-insensitive strains have more amino acid substitutions in the RGD motif of the VP1 region, and the sites T283S, V284M, and R288K in the VP1 region may be related to the temperature tolerance of CVA9.