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

Genotype F of Echovirus 25 with multiple recombination pattern have been persistently and extensively circulating in Chinese mainland

Echovirus 25 (E25), a member of the Enterovirus B (EV-B) species, can cause aseptic meningitis (AM), viral meningitis (VM), and acute flaccid paralysis (AFP). However, systematic studies on the molecular epidemiology of E25, especially those concerning its evolution and recombination, are lacking. In this study, 18 strains of E25, isolated from seven provinces of China between 2009 and 2018, were collected based on the Chinese hand, foot, and mouth disease (HFMD) surveillance network, and 95 sequences downloaded from GenBank were also screened. Based on the phylogenetic analysis of 113 full-length VP1 sequences worldwide, globally occurring E25 strains were classified into 9 genotypes (A-I), and genotype F was the dominant genotype in the Chinese mainland. The average nucleotide substitution rate of E25 was 6.08 × 10-3 substitutions/site/year, and six important transmission routes were identified worldwide. Seventeen recombination patterns were determined, of which genotype F can be divided into 9 recombination patterns. A positive selector site was found in the capsid protein region of genotype F. Recombination analysis and pressure selection analysis for genotype F showed multiple recombination patterns and evolution characteristics, which may be responsible for it being the dominant genotype in the Chinese mainland. This study provides a theoretical basis for the subsequent prevention and control of E25.

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⁻² substitution/site/year (95% highest posterior density (HPD): 8.17 × 10⁻³ to 1.4 × 10⁻² 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.

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

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

Global spatiotemporal transmission patterns of human enterovirus 71 from 1963 to 2019

Enterovirus 71 (EV71) can cause large outbreaks of hand, foot, and mouth disease (HFMD) and severe neurological diseases, which is regarded as a major threat to public health, especially in Asia-Pacific regions. However, the global spatiotemporal spread of this virus has not been identified. In this study, we used large sequence datasets and a Bayesian phylogenetic approach to compare the molecular epidemiology and geographical spread patterns of different EV71 subgroups globally. The study found that subgroups of HFMD presented global spatiotemporal variation, subgroups B0, B1, and B2 have caused early infections in Europe and America, and then subgroups C1, C2, C3, and C4 replaced B0-B2 as the predominant genotypes, especially in Asia-Pacific countries. The dispersal patterns of genotype B and subgroup C4 showed the complicated routes in Asia and the source might in some Asian countries, while subgroups C1 and C2 displayed more strongly supported pathways globally, especially in Europe. This study found the predominant subgroup of EV71 and its global spatiotemporal transmission patterns, which may be beneficial to reveal the long-term global spatiotemporal transmission patterns of human EV71 and carry out the HFMD vaccine development.

Molecular epidemiology and recombination of Enterovirus A71 in mainland China from 1987 to 2017

Enterovirus A71 (EV-A71) is an important pathogen of severe hand, foot, and mouth disease (HFMD) in young children. This study aimed to retrospectively analyze the molecular epidemiology and recombination of EV-A71 in mainland China during 1987–2017. Phylogenetic tree showed that besides the previously reported subgenotypes A, B5, C0, C2, C3, and C4, a new subgenotype C6 emerged in mainland China. Recombination analysis indicated that C4 EV-A71 was derived from a common ancestor as a “double-recombinant” virus by intertypic recombination between C EV-A71 and CVA4, CVA5, CVA14, and CVA16 strains in P3 region and intratypic recombination between C and B EV-A71 strains in P2 region. The B5 EV-A71 shared high similarity with C EV-A71 in P1 region while it contained an unidentified sequence in P2 and P3 regions with two possible recombination patterns: one occurred between C4 EV-A71 and CVA3, CVA5, CVA6, CVA10, and CVA12 stains with one breakpoint in 3C, and the other occurred between C1, C2, C3, and C5 EV-A71 and CVA4, CVA5, CVA14, and CVA16 strains with two breakpoints in the 2A/2B junction and 3C. The C2 EV-A71 was probably a recombinant virus between C4 EV-A71 and CVA8 strains with two breakpoints located in the 5′UTR and 2A/2B junction. Moreover, an incredible recombination of C6 EV-A71 occurred between C4 and C2 EV-A71 with multiple breakpoints. Thus, continuous studies on EV-A71 genome characteristics are still useful and essential for monitoring emergence of new viruses and preventing HFMD outbreaks.

Genetic recombination in fast-spreading coxsackievirus A6 variants: a potential role in evolution and pathogenicity

Hand, foot, and mouth disease (HFMD) is a common global epidemic. From 2008 onwards, many HFMD outbreaks caused by coxsackievirus A6 (CV-A6) have been reported worldwide. Since 2013, with a dramatically increasing number of CV-A6-related HFMD cases, CV-A6 has become the predominant HFMD pathogen in mainland China. Phylogenetic analysis based on the VP1 capsid gene revealed that subtype D3 dominated the CV-A6 outbreaks. Here, we performed a large-scale (near) full-length genetic analysis of global and Chinese CV-A6 variants, including 158 newly sequenced samples collected extensively in mainland China between 2010 and 2018. During the global transmission of subtype D3 of CV-A6, the noncapsid gene continued recombining, giving rise to a series of viable recombinant hybrids designated evolutionary lineages, and each lineage displayed internal consistency in both genetic and epidemiological features. The emergence of lineage –A since 2005 has triggered CV-A6 outbreaks worldwide, with a rate of evolution estimated at 4.17 × 10⁻³ substitutions site^-1 year⁻¹ based on a large number of monophyletic open reading frame sequences, and created a series of lineages chronologically through varied noncapsid recombination events. In mainland China, lineage –A has generated another two novel widespread lineages (–J and –L) through recombination within the enterovirus A gene pool, with robust estimates of occurrence time. Lineage –A, –J, and –L infections presented dissimilar clinical manifestations, indicating that the conservation of the CV-A6 capsid gene resulted in high transmissibility, but the lineage-specific noncapsid gene might influence pathogenicity. Potentially important amino acid substitutions were further predicted among CV-A6 variants. The evolutionary phenomenon of noncapsid polymorphism within the same subtype observed in CV-A6 was uncommon in other leading HFMD pathogens; such frequent recombination happened in fast-spreading CV-A6, indicating that the recovery of deleterious genomes may still be ongoing within CV-A6 quasispecies. CV-A6-related HFMD outbreaks have caused a significant public health burden and pose a great threat to children’s health; therefore, further surveillance is greatly needed to understand the full genetic diversity of CV-A6 in mainland China.

The Effect of Sample Bias and Experimental Artefacts on the Statistical Phylogenetic Analysis of Picornaviruses

Statistical phylogenetic methods are a powerful tool for inferring the evolutionary history of viruses through time and space. The selection of mathematical models and analysis parameters has a major impact on the outcome, and has been relatively well-described in the literature. The preparation of a sequence dataset is less formalized, but its impact can be even more profound. This article used simulated datasets of enterovirus sequences to evaluate the effect of sample bias on picornavirus phylogenetic studies. Possible approaches to the reduction of large datasets and their potential for introducing additional artefacts were demonstrated. The most consistent results were obtained using “smart sampling”, which reduced sequence subsets from large studies more than those from smaller ones in order to preserve the rare sequences in a dataset. The effect of sequences with technical or annotation errors in the Bayesian framework was also analyzed. Sequences with about 0.5% sequencing errors or incorrect isolation dates altered by just 5 years could be detected by various approaches, but the efficiency of identification depended upon sequence position in a phylogenetic tree. Even a single erroneous sequence could profoundly destabilize the whole analysis by increasing the variance of the inferred evolutionary parameters.