Characterization of VP1 sequence of Coxsackievirus A16 isolates by Bayesian evolutionary method

Genetic and Cross Neutralization Analyses of Coxsackievirus A16 Circulating in Taiwan from 1998 to 2021 Suggest Dominant Genotype B1 can Serve as Vaccine Candidate

Coxsackievirus A16 (CVA16) is well known for causing hand-foot-and-mouth disease (HFMD) and outbreaks were frequently reported in Taiwan in the past twenty years. The epidemiology and genetic variations of CVA16 in Taiwan from 1998 to 2021 were analyzed in this study. CVA16 infections usually occurred in early summer and early winter, and showed increased incidence in 1998, 2000-2003, 2005, 2007-2008, and 2010 in Taiwan. Little or no CVA16 was detected from 2017 to 2021. CVA16 infection was prevalent in patients between 1 to 3 years old. A total of 69 isolates were sequenced. Phylogenetic analysis based on the VP1 region showed that CVA16 subgenotype B1 was dominantly isolated in Taiwan from 1998 to 2019, and B2 was identified only from isolates collected in 1999 and 2000. There was a high frequency of synonymous mutations in the amino acid sequences of the VP1 region among CVA16 isolates, with the exception of position 145 which showed positive selection. The recombination analysis of the whole genome of CVA16 isolates indicated that the 5'-untranslated region and the non-structural protein region of CVA16 subgenotype B1 were recombined with Coxsackievirus A4 (CVA4) and enterovirus A71 (EVA71) genotype A, respectively. The recombination pattern of subgenotype B2 was similar to B1, however, the 3D region was similar to EVA71 genotype B. Cross-neutralization among CVA16 showed that mouse antisera from various subgenotypes viruses can cross-neutralize different genotype with high neutralizing antibody titers. These results suggest that the dominant CVA16 genotype B1 can serve as a vaccine candidate for CVA16.

Genome-wide population structure inferences of human coxsackievirus-A; insights the genotypes diversity and evolution

Coxsackievirus-A (CV-A) is a causative agent of Hand Foot Mouth Disease (HFMD) worldwide. It belongs to the Human Enterovirus genus of the family Picornaviridae. The genomics data availability of CV-A samples, isolated from human host across different continental regions, provide an excellent opportunity to study its genetic composition, diversity, and evolutionary events. The complete genome sequences of 424 CV-A isolates were analyzed through a model-based population genetic approach implemented in the STRUCTURE program. Twelve genetically distinct sub-populations were identified for CV-A isolates with a marked Fst distinction of 0.76991 (P-value = 0.00000). Besides, genetically admixed strains were characterized in the G-Id, G-IIIb clusters constituted by the CV-A12 and CV-A6 enterovirus serotypes. The serotypes depicted inter/intra-genotype recombination and episodic positive selection signatures in the structural and non-structural protein-coding regions. The observed genetic composition of CV-A samples was also deduced by the phylogenetic tree analyses, where a uniform genetic structure was inferred for most of the CV-A genotypes. However, the CV-A6 serotype samples genetically stratified into three sub-populations that may lead to the emergence of new lineages in future. These informations may implicate in planning the effective strategies to combat the coxsackievirus-A-mediated infection.

Molecular epidemiology of coxsackievirus A16 circulating in children in Beijing, China from 2010 to 2019

BACKGROUND

Coxsackievirus A16 (CVA16) is one of the major etiological agents of hand, foot and mouth disease (HFMD). This study aimed to investigate the molecular epidemiology and evolutionary characteristics of CVA16.

METHODS

Throat swabs were collected from children with HFMD and suspected HFMD during 2010-2019. Enteroviruses (EVs) were detected and typed by real-time reverse transcription-polymerase chain reaction (RT-PCR) and RT-PCR. The genotype, evolutionary rate, the most recent common ancestor, population dynamics and selection pressure of CVA16 were analyzed based on viral protein gene (VP1) by bioinformatics software.

RESULTS

A total of 4709 throat swabs were screened. EVs were detected in 3180 samples and 814 were CVA16 positive. More than 81% of CVA16-positive children were under 5 years old. The prevalence of CVA16 showed obvious periodic fluctuations with a high level during 2010-2012 followed by an apparent decline during 2013-2017. However, the activities of CVA16 increased gradually during 2018-2019. All the Beijing CVA16 strains belonged to sub-genotype B1, and B1b was the dominant strain. One B1c strain was detected in Beijing for the first time in 2016. The estimated mean evolutionary rate of VP1 gene was 4.49 × 10^-3 substitution/site/year. Methionine gradually fixed at site-23 of VP1 since 2012. Two sites were detected under episodic positive selection, one of which (site-223) located in neutralizing linear epitope PEP71.

CONCLUSIONS

The dominant strains of CVA16 belonged to clade B1b and evolved in a fast evolutionary rate during 2010-2019 in Beijing. To provide more favorable data for HFMD prevention and control, it is necessary to keep attention on molecular epidemiological and evolutionary characteristics of CVA16.

Genetic Variation of Multiple Serotypes of Enteroviruses Associated with Hand, Foot and Mouth Disease in Southern China

Enteroviruses (EVs) species A are a major public health issue in the Asia–Pacific region and cause frequent epidemics of hand, foot and mouth disease (HFMD) in China. Mild infections are common in children; however, HFMD can also cause severe illness that affects the central nervous system. To molecularly characterize EVs, a prospective HFMD virological surveillance program was performed in China between 2013 and 2016. Throat swabs, rectal swabs and stool samples were collected from suspected HFMD patients at participating hospitals. EVs were detected using generic real-time and nested reverse transcription-polymerase chain reactions (RT-PCRs). Then, the complete VP1 regions of enterovirus A71 (EV-A71), coxsackievirus A16 (CVA16) and CVA6 were sequenced to analyze amino acid changes and construct a viral molecular phylogeny. Of the 2836 enrolled HFMD patients, 2,517 (89%) were EV positive. The most frequently detected EVs were CVA16 (32.5%, 819), CVA6 (31.2%, 785), and EV-A71 (20.4%, 514). The subgenogroups CVA16_B1b, CVA6_D3a and EV-A71_C4a were predominant in China and recombination was not observed in the VP1 region. Sequence analysis revealed amino acid variations at the 30, 29 and 44 positions in the VP1 region of EV-A71, CVA16 and CVA6 (compared to the respective prototype strains BrCr, G10 and Gdula), respectively. Furthermore, in 21 of 24 (87.5%) identified EV-A71 samples, a known amino acid substitution (D31N) that may enhance neurovirulence was detected. Our study provides insights about the genetic characteristics of common HFMD-associated EVs. However, the emergence and virulence of the described mutations require further investigation.

Genome Analysis of Coxsackievirus A4 Isolates From Hand, Foot, and Mouth Disease Cases in Shandong, China

Coxsackievirus A4 (CVA4) is one of the most prevalent pathogens associated with hand, foot and mouth disease (HFMD), an acute febrile illness in children, and is also associated with acute localized exanthema, myocarditis, hepatitis and pancreatitis. Despite this, limited CVA4 genome sequences are currently available. Herein, complete genome sequences from CVA4 strains (n = 21), isolated from patients with HFMD in Shandong province, China between 2014 and 2016, were determined and phylogenetically characterized. Phylogenetic analysis of the VP1 gene from a larger CVA4 collection (n = 175) showed that CVA4 has evolved into four separable genotypes: A, B, C, and D; and genotype D could be further classified in to two sub-genotypes: D1 and D2. Each of the 21 newly described genomes derived from isolates that segregated with sub-genotype D2. The CVA4 genomes displayed significant intra-genotypic genetic diversity with frequent synonymous substitutions occurring at the third codon positions, particularly within the P2 region. However, VP1 was relatively stable and therefore represents a potential target for molecular diagnostics assays and also for the rational design of vaccine epitopes. The substitution rate of VP1 was estimated to be 5.12 × 10^-3 substitutions/site/year, indicative of ongoing CVA4 evolution. Mutations at amino acid residue 169 in VP1 gene may be responsible for differing virulence of CVA4 strains. Bayesian skyline plot analysis showed that the population size of CVA4 has experienced several dynamic fluctuations since 1948. In summary, we describe the phylogenetic and molecular characterization of 21 complete genomes from CVA4 isolates which greatly enriches the known genomic diversity of CVA4 and underscores the need for further surveillance of CVA4 in China.

An Environmental Surveillance in Uruguay Reveals the Presence of Highly Divergent Types of Human Enterovirus Species C and a High Frequency of Species A and B Types

Information about Human Enterovirus circulation in Uruguay is scarce. The aim of this study was to generate the first description about their circulation in the country through the study of sewage samples collected before and after the switch from Oral Poliovirus Vaccine to Inactivated Poliovirus Vaccine. Viruses were concentrated by an adsorption-elution to a negatively charged membrane, and real-time quantitative PCR and qualitative PCR methods were used to detect, quantify, and characterize enteroviruses. Positive samples were inoculated in RD cells and two passages were performed. Additionally, RD+ samples were subsequently passed onto L20B cells. Human Enteroviruses were detected in 67.6% of the samples, with concentrations between 4.9 and 6.6 Log₁₀ genomic copies per liter. 10% of positive samples replicated in RD cells, of which none in L20B cells. Molecular characterization of Human Enterovirus strains directly detected from sewage sample concentrates allowed the identification of highly divergent members of species C such as Enterovirus C99 and Coxsackievirus A13, as well as the frequent detection of species A and B members (particularly Coxsackievirus A16 and Echovirus 6, respectively). Other detected types were Coxsackievirus A2, A22, B1, B5, Echovirus 5, and 9. The characterization of viruses isolated in cell culture revealed the presence of Echovirus 6 and Coxsackievirus B3. Despite the absence of poliovirus, a wide circulation of different enterovirus types was evidenced in Uruguayan sewage samples, highlighting that the local populations are exposed to different kinds of diseases originated by several human enterovirus.

Molecular epidemiology and phylogenetics of human enteroviruses: Is there a forest behind the trees?

Enteroviruses are among the best studied small non-enveloped enteric RNA viruses. Most enteroviruses are easy to isolate in cell culture, and many non-polio enterovirus strains were archived worldwide as a byproduct of the WHO poliovirus surveillance system. Common outbreaks and epidemics, most prominently the epidemic of hand-foot-and-mouth disease with severe neurological complications in East and South-East Asia, justify practical interest of non-polio enteroviruses. As a result, there are over 50 000 enterovirus nucleotide sequences available in GenBank. Technical possibilities have been also improving, as Bayesian phylogenetic methods with an integrated molecular clock were introduced a decade ago and provided unprecedented opportunities for phylogenetic analysis. As a result, hundreds of papers were published on the molecular epidemiology of enteroviruses. This review covers the modern methodology, structure, and biases of the sequence dataset available in GenBank. The relevance of the subtype classification, findings of co-circulation of multiple genetic variants, previously unappreciated complexity of viral populations, and global evolutionary patterns are addressed. The most relevant conclusions and prospects for further studies on outbreak emergence mechanisms are discussed.

Molecular evolution of hepatitis B vaccine escape variants in China, during 2000-2016

Hepatitis B vaccine escape variants are the main threat to hepatitis B virus (HBV) infection in vaccination era worldwide. With 215 genotype B HBV and 313 genotype C HBV vaccine escape variants isolated from China during 2000-2016, we reported that genotype B HBV vaccine escape strains diverged in ∼1997 (95% HPD; 1987-2005), while genotype C HBV vaccine escape strains diverged in ∼1976 (95% HPD; 1955-2003). Additionally, the p-distance of genotype C HBV vaccine escape strains was 0.0291±0.0169, which was significantly higher than that in the genotype B HBV (t=131.02, p<0.05). However, genotype B HBV vaccine escape strains evolved more rapidly than genotype C HBV (2.103×10^-3 vs 1.083×10^-3 substitutions/site/year). Bayesian skyline plot analysis showed that the populations of genotype C HBV vaccine escape strains fluctuated more than those in genotype B HBV. Four sites (A5T/S, L21S, T/A126S and T/N131I/A) and 13 sites (N3S, T5A, G10Q/R/E, L21S, T47K/A/V, L98V/P, I/S126N/V/T, Q129H/R/L, T131P/I/N/A, G145A/R, L175S/F, L213I/S, V224A/G) were found to be under positive selection in genotype B and C HBV vaccine escape strains, respectively. More importantly, N3S, L21S, T47K, L98V, I/S126T and L213I mutations were detected in 1 (2.5%), 1 (2.5%), 1 (2.5%), 3 (7.5%), 1 (2.5%), 1 (2.5%) genotype C HBV infected Chinese younger with neonatal HBV vaccination, respectively. Therefore, our results should be valuable in further understanding the molecular evolution of HBV and providing new ideas for the elimination of HBV infection.