Genetic characterization and molecular epidemiological analysis of novel enterovirus EV-B80 in China

Genomic characteristics of human respiratory syncytial virus from children in China during 2017-2020

Acute lower respiratory tract infections (ALRTIs) are a leading cause of mortality in young children worldwide due to human respiratory syncytial virus (RSV). The aim of this study was to monitor genetic variations in RSV and provide genomic data support for RSV prevention and control. A total of 105 complete RSV genome sequences were determined during 2017-2020. Phylogenetic analysis showed that all of the RSVA sequences were of genotype ON1, and all of the RSVB sequences were of genotype BA9. Notably, a phylogenetic tree based on the whole genome had more branches than a tree based on the G gene. In comparison to the RSV prototype sequences, 71.43% (50/70) of the ON1 sequences had five amino acid substitutions (T113I, V131N, N178G, H258Q, and H266L) that occurred simultaneously, and 68.57% (24/35) of the BA9 genotype sequences had 12 amino acid substitutions, four of which (A131T, T137I, T288I, and T310I) occurred simultaneously. In the F gene, there were 19 amino acid substitutions, which were mainly located in the antigenic sites Ø, II, V, and VII. Other amino acid substitutions were found in the NS1, NS2, P, SH, and L proteins. No significant evidence of recombination was found in any of the sequences. These findings provide important data that will be useful for prevention, control, and vaccine development against RSV.

Changing Epidemiology of Hand, Foot, and Mouth Disease Causative Agents and Contributing Factors

Hand, foot, and mouth disease (HFMD) is a common viral infection primarily affecting children. It causes vesicles on the skin and inside the mouth. Although most cases get better on their own, severe cases can lead to complications such as brain stem encephalitis, meningoencephalitis, acute flaccid paralysis, and pulmonary edema. Hand, foot, and mouth disease is caused by various enteroviruses, with enterovirus A71 (EV-A71) and coxsackievirus A16 being the most common. However, recent studies have shown a shift in the molecular epidemiology of HFMD-causing pathogens, with coxsackievirus A6 and coxsackievirus A10 causing more infections. In addition, extensive recombination events have been identified among enterovirus strains, which may have a role in faster evolution and extinction of dominant enterovirus serotypes. Other strains of enterovirus can also cause severe complications, and there has been an increase in mortality associated with brain stem encephalitis in children under 3 years of age and teenagers. Currently, there are no effective antiviral therapies available to treat enterovirus infections. Vaccines against EV-A71 have been approved and are now used in mainland China. Studying the changing epidemiology of HFMD pathogens and the evolution patterns of its causative agents is crucial in developing effective prevention and control strategies. Increased interest in the molecular epidemiology of HFMD causative agents has led to a better understanding of the critical drivers of HFMD outbreaks, which can inform efforts to prevent and control the disease.

Molecular Genotyping of Circulating Enterovirus in the Lazio Region from 2012 to 2023

Enteroviruses (EVs) are ubiquitous viruses that circulate worldwide, causing sporadic or epidemic infections, typically during the summer and fall. They cause a broad spectrum of illnesses, ranging from an unspecified febrile clinical presentation to a severe illness. EVs are recognized to be the most frequent etiological agents of aseptic meningitis in children. However, as the infection is usually mild and self-limiting, it remains underestimated, and the epidemiology of EVs is poorly understood. To date, no vaccine or effective therapy for all types of enteroviruses is available, and EVs constitute a public health concern. Here, we investigated the molecular epidemiology of EV strains circulating in the Lazio region over a 10-year time span (2012-2023) by using a sequence-typing approach and phylogenetic analysis. The epidemiological trend of EV infection has undergone changes during the SARS-CoV-2 pandemic (2020-2021), which resulted in a modification in terms of the number of diagnosed cases and seasonality. From 2022, the circulation of EVs showed a behavior typical of the pre-pandemic period, although changes in predominantly circulating strains have been noted. Both epidemic and sporadic circulation events have been characterized in the Lazio region. Further analyses are needed to better characterize any strain with higher potential pathogenic power and to identify possible recombinant strains.

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.

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.

Phylogenetic characteristics and recombination analysis of echovirus 5 associated with severe acute respiratory infection in China

IMPORTANCE

This study is the first report of echovirus 5 (E5) associated with severe acute respiratory infection and obtained the first E5 whole-genome sequence in China. Combined with the sequences available in the GenBank database, the first genotyping, phylogenetic characteristics, recombination, and genetic evolutionary analysis of E5 was performed in this study. Our findings providing valuable information on global E5 molecular epidemiology.

Molecular Epidemiology and Evolution of Coxsackievirus A14

As the proportion of non-enterovirus 71 and non-coxsackievirus A16 which proportion of composition in the hand, foot, and mouth pathogenic spectrum gradually increases worldwide, the attention paid to other enteroviruses has increased. As a member of the species enterovirus A, coxsackievirus A14 (CVA14) has been epidemic around the world until now since it has been isolated. However, studies on CVA14 are poor and the effective population size, evolutionary dynamics, and recombination patterns of CVA14 are not well understood. In this study, 15 CVA14 strains were isolated from HFMD patients in mainland China from 2009 to 2019, and the complete sequences of CVA14 in GenBank as research objects were analyzed. CVA14 was divided into seven genotypes A-G based on an average nucleotide difference of the full-length VP1 coding region of more than 15%. Compared with the CVA14 prototype strain, the 15 CVA14 strains showed 84.0-84.7% nucleotide identity in the complete genome and 96.9-97.6% amino acid identity in the encoding region. Phylodynamic analysis based on 15 CVA14 strains and 22 full-length VP1 sequences in GenBank showed a mean substitution rate of 5.35 × 10-3 substitutions/site/year (95% HPD: 4.03-6.89 × 10-3) and the most recent common ancestor (tMRCA) of CVA14 dates back to 1942 (95% HPD: 1930-1950). The Bayesian skyline showed that the effective population size had experienced a decrease-increase-decrease fluctuation since 2004. The phylogeographic analysis indicated two and three possible migration paths in the world and mainland China, respectively. Four recombination patterns with others of species enterovirus A were observed in 15 CVA14 strains, among which coxsackievirus A2 (CVA2), coxsackievirus A4 (CVA4), coxsackievirus A6 (CVA6), coxsackievirus A8 (CVA8), and coxsackievirus A12 (CVA12) may act as recombinant donors in multiple regions. This study has filled the gap in the molecular epidemiological characteristics of CVA14, enriched the global CVA14 sequence database, and laid the epidemiological foundation for the future study of CVA14 worldwide.

Molecular characterization of a novel clade echovirus 3 isolated from patients with hand-foot-and-mouth disease in southwest China

Echovirus 3 (E3) belongs to the species Enterovirus B. Currently, three nearly whole-genome sequences of E3 are available in GenBank in China. In this study, we determined the whole genomic sequences of six E3 strains isolated from the stools of patients with hand-foot-and-mouth disease in Southwest China in 2022. Their nucleotide and amino acid sequences shared 82.1%-86.4% and 96.6%-97.2% identity with the prototype Morrisey strain, respectively, and showed 87.1% and 97.2% mutual identity. The six E3 strains are not clustered with other Chinese strains and formed a novel subgenotype (C6) with the recent American and British strains. Recombination analyses revealed that intertype recombination had occurred in the 2 C and 3D regions of the six E3 strains with coxsackieviruses B5 and B4, respectively. This study augments the nearly whole-genome sequences of E3 in the GenBank database and extends the molecular characterization of this virus in China.

Molecular Identification of Enteric Viruses in Domestic Animals in Northeastern Gabon, Central Africa

Astroviruses (AstVs), enteroviruses (EVs), and caliciviruses (CaVs) infect several vertebrate taxa. Transmitted through the fecal-oral route, these enteric viruses are highly resistant and can survive in the environment, thereby increasing their zoonotic potential. Here, we screened for AstVs, EVs, and CaVs to investigate the role of domestic animals in the emergence of zoonoses, because they are situated at the human/wildlife interface, particularly in rural forested areas in Central Africa. Rectal swabs were obtained from 123 goats, 41 sheep, and 76 dogs in 10 villages located in northeastern Gabon. Extracted RNA reverse-transcribed into cDNA was used to detect AstVs, EVs, and CaVs by amplification of the RNA-dependent RNA polymerase (RdRp), or capsid protein (VP1) gene using PCR. A total of 23 samples tested positive, including 17 goats for AstVs, 2 goats, 2 sheep, 1 dog for EVs, and 1 dog for CaVs. Phylogenetic analyses revealed that AstV RdRp sequences clustered with sheep-, goat-, or bovine-related AstVs. In addition, one goat and two sheep VP1 sequences clustered with caprine/ovine-related Evs within the Enterovirus G species, and the CaV was a canine vesivirus. However, human-pathogenic Evs, EV-B80 and EV-C99, were detected in goats and dogs, raising questions on the maintenance of viruses able to infect humans.

Molecular evolution and antigenic drift of type 3 iVDPVs excreted from a patient with immunodeficiency in Ningxia, China

A 2.5-year-old pediatric patient with acute flaccid paralysis was diagnosed with primary immunodeficiency (PID) in Ningxia Province, China, in 2011. Twelve consecutive stool specimens were collected from the patient over a period of 10 months (18 February 2011 to 20 November 2011), and 12 immunodeficiency vaccine-derived poliovirus (iVDPV) strains (CHN15017-1 to CHN15017-12) were subsequently isolated. Nucleotide sequencing analysis of the plaque-purified iVDPVs revealed 2%-3.5% VP1-region differences from their parental Sabin 3 strain. Full-length genome sequencing showed they were all Sabin 3/Sabin 1 recombinants, sharing a common 2C-region crossover site, and the two key determinants of attenuation (U472C in the 5' untranslated region and T2493C in the VP1 region) had reverted. Temperature-sensitive experiments demonstrated that the first two iVDPV strains partially retained the temperature-sensitive phenotype's nature, while the subsequent ten iVDPV strains distinctly lost it, possibly associated with increased neurovirulence. Nineteen amino-acid substitutions were detected between 12 iVDPVs and the parental Sabin strain, of which only one (K1419R) was found on the subsequent 10 iVDPV isolates, suggesting this site's potential as a temperature-sensitive determination site. A Bayesian Monte Carlo Markov Chain phylogenetic analysis based on the P1 coding region yielded a mean iVDPV evolutionary rate of 1.02 × 10^-2 total substitutions/site/year, and the initial oral-polio-vaccine dose was presumably administered around June 2009. Our findings provide valuable information regarding the genetic structure, high-temperature growth sensitivity, and antigenic properties of iVDPVs following long-term evolution in a single PID patient, thus augmenting the currently limited knowledge regarding the dynamic changes and evolutionary pathway of iVDPV populations with PID during long-term global replication.

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.

Detection of recombinant breakpoint in the genome of human enterovirus E11 strain associated with a fatal nosocomial outbreak

BACKGROUND

The aim of this study was to characterize the genome of a recombinant Enterovirus associated with severe and fatal nosocomial infection; it was typed as Echovirus 11 (E-11) according to the VP1 gene. Enterovirus infection is generally asymptomatic and self-limited, but occasionally it may progress to a more severe clinical manifestation, as in the case described here. Recombination plays a crucial role in the evolution of Enteroviruses (EVs) and has been recognized as the main driving force behind the emergence of epidemic strains associated with severe infection. Therefore, it is of utmost importance to monitor the circulation of recombinant strains for surveillance purposes.

METHODS

Enterovirus-RNA was detected in the serum and liver biopsy of patients involved in the nosocomial cluster by commercial One-Step qRT-PCR method and the Enterovirus strains were isolated in vitro. The EVs typing was determined by analyzing the partial-length of the 5'UTR and VP1 sequences with the web-based open-access Enterovirus Genotyping Tool Version 0.1. The amplicons targeting 5'UTR, VP1 and overlapping fragments of the entire genome were sequenced with the Sanger method. Phylogenetic analysis was performed comparing the VP1 and the full-genome sequences of our strains against an appropriate reference set of Enterovirus prototypes of the Picornaviridae genera and species retrieved from the Enterovirus Genotyping Tool. Recombination analysis was performed using RDP4 software.

RESULTS

The Neighbor-Joining tree of the VP1 gene revealed that the 4 patients were infected with an identical molecular variant of Echovirus 11 (E-11). While the phylogenetic and the RDP4 analysis of the full-genome sequences provided evidence that it was a chimeric strain between an E-11 and a Coxsackievirus B (CV-B).

CONCLUSIONS

The chimeric structure of the E-11 genome might have contributed to the severe infection and epidemic feature of the strain, but further biological characterizations are needed. The evidence reported in this study, highlights the limit of typing techniques based on the VP1 gene, as they fail to identify the emergence of recombinant strains with potentially more pathogenic or epidemic properties, thus providing only partial information on the epidemiology and pathogenesis of Enteroviruses.

Monsavirus in monkey rectal swab and throat swab specimens in China: Proposal for Posaliviridae as a new family in Picornavirales

Posa-like viruses have been detected in the fecal samples of several host species and are considered unclassified members of Picornavirales. Here, we identified genomic fragments of novel posa-like viruses (monsaviruses) in monkey specimens through next generation sequencing and obtained 11 full-length genomes. This monsavirus shared 88.5-89.2% nucleotide similarity with the Tottori-HG1 strain (GenBank accession LC123275). In total, 713 nucleotide polymorphism sites were identified, indicating their persistent evolution during circulation. The genomic organization and phylogenetic relationship of monsavirus were determined. Subsequent phylogenetic analysis of the conserved replication block of Hel-Pro-RdRp and core RNA-dependent RNA polymerase domain-based analysis of posa-like viruses showed significant separation compared with other known families. Further, posa-like virus genomes possessed the classical replication block of picornavirus in the 5' part of genome and picorna-like capsid domains at the structural coding region of 3' part of genome. Based on these results, we proposed the new family Posaliviridae, within Picornavirales. Four genera, which showed 68.6-75.5% amino acid distances but similar genomic organization including the conserved replication block of Hel-Pro-RdRp, the same order of the genomic coding region, and picorna-like capsid domains, were identified. The flexible genomic organization strategy and a large evolutionary scale of Posaliviridae was explicit. This study provides novel information on monsaviruses and important taxonomic data for the family Posaliviridae.

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.

A novel interspecies recombinant enterovirus (Enterovirus A120) isolated from a case of acute flaccid paralysis in China

EV-A120 is a recently identified serotype of the enterovirus A species. Only one full-length genomic sequence is currently available in GenBank, and very few studies have been conducted on EV-A120 globally. Thus, additional information and research on EV-A120 are needed to explore its genetic characteristics, phylogeny, and relationship with enteroviral disease. In this study, we report the phylogenetic characteristics of a EV-A120 strain (Q0082/XZ/CHN/2000) from Tibet, China. The amino acid sequence similarity and nucleotide sequence similarity of the full-length genomic sequence of this EV-A120 strain and the EV-A120 prototype strain were 96.3% and 79.9%, respectively, showing an evolutionary trend. Recombination analysis found intraspecies recombination in the 5′ -UTR, 2B, 2C, and 3D regions. Serum neutralization testing of the EV-A120 (Q0082) strain was also carried out. Low serum-positive rates and geometric mean titres (GMTs) indicated that the extent of EV-A120 transmission and exposure in the population was very limited compared with that in the outbreaks of EV-A71 and CV-A16 in China since 2008. The EV-A120 strain (Q0082) is non-temperature sensitive, indicating its potential to spread in the population. In summary, this study reports the full-length genomic sequence of EV-A120 and provides important information for its global molecular epidemiology.

Molecular Epidemiological, Serological, and Pathogenic Analysis of EV-B75 Associated With Acute Flaccid Paralysis Cases in Tibet, China

Enterovirus B75 (EV-B75) is a newly identified serotype of the enterovirus B species. To date, only 112 cases related to EV-B75 have been reported worldwide, and research on EV-B75 is still limited with only two full-length genome sequences available in GenBank. The present study reported seven EV-B75 sequences from a child with acute flaccid paralysis and six asymptomatic close contacts in Shigatse, Tibet. Phylogenetic analysis revealed that the Tibetan strain was possibly imported from neighboring India. Seroepidemiological analyses indicated that EV-B75 has not yet caused a large-scale epidemic in Tibet. Similarity plots and boot scanning analyses revealed frequent intertypic recombination in the non-structural region of all seven Tibet EV-B75 strains. All seven Tibetan strains were temperature-sensitive, suggesting their poor transmissibility in the environment. Overall, though the seven Tibetan strains did not cause large-scale infection, prevention and control of the novel enterovirus cannot be underestimated.

Genomic epidemiology of coxsackievirus A16 in mainland of China, 2000-18

Hand, foot, and mouth disease (HFMD), which is a frequently reported and concerning disease worldwide, is a severe burden on societies globally, especially in the countries of East and Southeast Asia. Coxsackievirus A16 (CV-A16) is one of the most important causes of HFMD and a severe threat to human health, especially in children under 5 years of age. To investigate the epidemiological characteristics, spread dynamics, recombinant forms (RFs), and other features of CV-A16, we leveraged the continuous surveillance data of CV-A16-related HFMD cases collected over an 18-year period. With the advent of the EV-A71 vaccine since 2016, which targeted the EV-A71-related HFMD cases, EV-A71-related HFMD cases decreased dramatically, whereas the CV-A16-related HFMD cases showed an upward trend from 2017 to October 2019. The CV-A16 strains observed in this study were genetically related and widely distributed in the mainland of China. Our results show that three clusters (B1a-B1c) existed in the mainland of China and that the cluster of B1b dominates the diffusion of CV-A16 in China. We found that eastern China played a decisive role in seeding the diffusion of CV-A16 in China, with a more complex and variant transmission trend. Although EV-A71 vaccine was launched in China in 2016, it did not affect the genetic diversity of CV-A16, and its genetic diversity did not decline, which confirmed the epidemiological surveillance trend of CV-A16. Two discontinuous clusters (2000-13 and 2014-18) were observed in the full-length genome and arranged along the time gradient, which revealed the reason why the relative genetic diversity of CV-A16 increased and experienced more complex fluctuation model after 2014. In addition, the switch from RFs B (RF-B) and RF-C co-circulation to RF-D contributes to the prevalence of B1b cluster in China after 2008. The correlation between genotype and RFs partially explained the current prevalence of B1b. This study provides unprecedented full-length genomic sequences of CV-A16 in China, with a wider geographic distribution and a long-term time scale. The study presents valuable information about CV-A16, aimed at developing effective control strategies, as well as a call for a more robust surveillance system, especially in the Asia-Pacific region.

Global Spread of the B5 Subgenotype EV-A71 and the Phylogeographical Analysis of Chinese Migration Events

The subgenotype B5 of EV-A71 is a widely circulating subgenotype that frequently spreads across the globe. Several outbreaks have occurred in nations, such as Malaysia, Thailand, Vietnam, and Japan. Appearing first in Taiwan, China, the subgenotype has been frequently reported in mainland of China even though no outbreaks have been reported so far. The current study reconstructed the migration of the B5 subgenotype of EV-A71 in China via phylogeographical analysis. Furthermore, we investigated its population dynamics in order to draw more credible inferences. Following a dataset cleanup of B5 subgenotype of EV-A71, we detected earlier B5 subgenotypes of EV-A71 sequences that had been circulating in Malaysia and Singapore since the year 2000, which was before the 2003 outbreak that occurred in Sarawak. The Bayesian inference indicated that the most recent common ancestor of B5 subgenotype EV-A71 appeared in September, 1994 (1994.75). With respect to the overall prevalence, geographical reconstruction revealed that the B5 subgenotype EV-A71 originated singly from single-source cluster and subsequently developed several active lineages. Based on a large amount of data that was accumulated, we conclude that the appearance of the B5 subgenotype of EV-A71 in mainland of China was mainly due to multiple migrations from different origins.

The Husavirus Posa-Like Viruses in China, and a New Group of Picornavirales

Novel posa-like viral genomes were first identified in swine fecal samples using metagenomics and were designated as unclassified viruses in the order Picornavirales. In the present study, nine husavirus strains were identified in China. Their genomes share 94.1–99.9% similarity, and alignment of these nine husavirus strains identified 697 nucleotide polymorphism sites across their full-length genomes. These nine strains were directly clustered with the Husavirus 1 lineage, and their genomic arrangement showed similar characteristics. These posa-like viruses have undergone a complex evolutionary process, and have a wide geographic distribution, complex host spectrum, deep phylogenetic divergence, and diverse genomic organizations. The clade of posa-like viruses forms a single group, which is evolutionarily distinct from other known families and could represent a distinct family within the Picornavirales. The genomic arrangement of Picornavirales and the new posa-like viruses are different, whereas the posa-like viruses have genomic modules similar to the families Dicistroviridae and Marnaviridae. The present study provides valuable genetic evidence of husaviruses in China, and clarifies the phylogenetic dynamics and the evolutionary characteristics of Picornavirales.

Molecular typing and characterization of a novel genotype of EV-B93 isolated from Tibet, China

EV-B93 is a novel serotype within the Enterovirus B species and is uncommon worldwide. Currently, only one full-length genomic sequence (the prototype strain) has been deposited in the GenBank database. In this study, three EV-B93 were identified, including one from an acute flaccid paralysis (AFP) patient (named 99052/XZ/CHN/1999, hereafter XZ99052) and two from healthy children (named 99096/XZ/CHN/1999 and 99167/XZ/CHN/1999, hereafter XZ99096 and XZ99167, respectively) from Tibet in 1999 during the polio eradication program. The identity between the nucleotide and amino acid sequences of the Tibet EV-B93 strain and the EV-B93 prototype strain is 83.2%–83.4% and 96.8%–96.9%, respectively. The Tibet EV-B93 strain was found to have greater nucleotide sequence identity in the P3 region to another enterovirus EV-B107 as per a phylogenetic tree analysis, which revealed that recombination occurred. Seroepidemiology data showed that EV-B93 has not produced an epidemic in Tibet and there may be susceptible individuals. The three Tibet EV-B93 strains are temperature-resistant with prognosticative virulence, suggesting the possibility of a potential large-scale outbreak of EV-B93. The analyzed EV-B93 strains enrich our knowledge about this serotype and provide valuable information on global EV-B93 molecular epidemiology. What is more, they permit the appraisal of the serotype's potential public health impact and aid in understanding the role of recombination events in the evolution of enteroviruses.

A Large-Scale Outbreak of Echovirus 30 in Gansu Province of China in 2015 and Its Phylodynamic Characterization

Background

Echovirus 30 (E-30) has been investigated and reported worldwide and is closely associated with several infectious diseases, including encephalitis; myocarditis; and hand, foot, and mouth disease. Although many E-30 outbreaks associated with encephalitis have been reported around the world, it was not reported in northwest China until 2015.

Methods

The clinical samples, including the feces, serum, throat swabs, and cerebrospinal fluid, were collected for this study and were analyzed for diagnosis. E-30 was isolated and processed according to the standard procedures. The epidemiological and phylogenetic analysis were performed to indicate the characteristics of E-30 outbreaks and phylodynamics of E-30 in China.

Results

The E-30 outbreaks affected nine towns of Gansu Province in 2015, starting at a school of Nancha town and spreading to other towns within 1 month. The epidemiological features showed that children aged 6–15 years were more susceptible to E-30 infection. The genotypes B and C cocirculated in the world, whereas the latter dominated the circulation of E-30 in China. The genome sequences of this outbreak present 99.3–100% similarity among these strains, indicating a genetic-linked aggregate outbreak of E-30 in this study. Two larger genetic diversity expansions and three small fluctuations of E-30 were observed from 1987 to 2016 in China, which revealed the oscillating patterns of E-30 in China. In addition, the coastal provinces of China, such as Zhejiang, Fujian, and Shandong, were initially infected, followed by other parts of the country. The E-30 strains isolated from mainland of China may have originated from Taiwan of China in the last century.

Conclusion

The highly similar E-30 genomes in this outbreak showed an aggregate outbreak of E-30, with nine towns affected. Our results suggested that, although the genetic diversity of E-30 oscillates, the dominant lineages of E-30 in China has complex genetic transmission. The coastal provinces played an important role in E-30 spread, which implied further development of effective countermeasures. This study provides a further insight into the E-30 outbreak and transmission and illustrates the importance of valuable surveillance in the future.

Tracking echovirus eleven outbreaks in Guangdong, China: a metatranscriptomic, phylogenetic, and epidemiological study

In April 2019, a suspect cluster of enterovirus cases was reported in a neonatology department in Guangdong, China, resulting in five deaths. We aimed to investigate the pathogen profiles in fatal cases, the circulation and transmission pattern of the viruses by combining metatranscriptomic, phylogenetic, and epidemiological analyses. Metatranscriptomic sequencing was used to characterize the enteroviruses. Clinical and environmental surveillance in the local population was performed to understand the prevalence and genetic diversity of the viruses in the local population. The possible source(s), evolution, transmission, and recombination of the viruses were investigated by incorporating genomes from the current outbreak, from local retrospective surveillance, and from public databases. Metatranscriptomic analysis identified Echovirus 11 (E11) in three fatal cases. Seroprevalence of neutralization antibody to E11 was 35 to 44 per cent in 3–15 age groups of general population, and the viruses were associated with various clinical symptoms. From the viral phylogeny, nosocomial transmissions were identified and all E11 2019 outbreak strains were closely related with E11 strains circulating in local population 2017–19. Frequent recombination occurred among the 2019 Guangdong E11 outbreak strains and various genotypes in enterovirus B species. This study provides an example of combining advanced genetic technology and epidemiological surveillance in pathogen diagnosis, source(s), and transmission tracing during an infectious disease outbreak. The result highlights the hidden E11 circulation and the risk of viral transmission and infection in the young age population in China. Frequent recombination between Guangdong-like strains and other enterovirus genotypes also implies the prevalence of these emerging E11 strains.

Phylogenetic characteristics and molecular epidemiological analysis of novel enterovirus EV-B83 isolated from Tibet, China

Enterovirus B83 (EV-B83) is a new member of the enterovirus B group. Currently, there are only two full-length genomic sequences of EV-B83 in the GenBank database and few VP1 region sequences. The aetiology and epidemiology of EV-B83 is unclear. 24 stool specimens were collected from twelve AFP patients and 298 stool specimens were collected from 298 healthy children in support of polio eradication activities in Tibet in 1999. Two polioviruses (isolated by L20B cell) and one non-polio enterovirus (isolated by RD cell) were isolated from AFP patients and nine polioviruses (isolated by L20B cell) and 90 non-polio enteroviruses (isolated by RD cell) were isolated from health children. Through molecular typing, we confirmed that the six of non-polio enteroviruses belong to EV-B83. The sequence similarity between the VP1 region of the Tibet isolates and that of the EV-B83 prototype strain was 80%. The maximum-likelihood phylogenetic tree of the partial VP1 region in EV-B83 demonstrated that EV-B83 formed four genotypes globally during the evolution process. The six Tibet EV-B83 strains formed the D genotype alone. Recombination analysis of Tibet EV-B83 showed that CV-B4, CV-A9, EV-B80, and EV-B106 may act as recombinant donors in multiple regions. The serum neutralization test showed that the antibody-positive rate was 58.8% and GMT was 1:19.70, which was higher than the previously reported results of EV-B106 and EV-B80. Temperature sensitivity test results showed that the six Tibet EV-B83 strains were temperature-insensitive with stronger virulence and potential infectivity, which was consistent with the results of the serum neutralization test. This study enriched the genome-wide sequence, epidemiological characteristics, and provided basic data for the follow-up study of EV-B83.