Origin and evolution analysis and genetic characteristics of echovirus 9 in China

Background: Echovirus 9 (E9) is associated with a wide variety of diseases and medical conditions, and the clinical symptoms of sporadic cases caused by E9 often are severe. With a high global prevalence, E9 has caused multiple outbreaks worldwide. However, little is known about the genetic and geographic population dynamics of E9.

Method: A total of 131 VP1 gene sequences, including15 generated in this study and 116 obtained from GenBank, were used to coestimate time-resolved phylogenies to infer viral evolution and transmission in worldwide. Overlapping fragments representing whole genomes were amplified by reverse transcription polymerase chain reaction (RT-PCR) using specific primers. Then, we reported the genetic characteristics of fifteen E9 strains in the Chinese Mainland. Similarity plots and bootscanning analysis were used to determine recombination patterns of E9.

Results: The estimated mean evolutionary rate of global E9 VP1 gene was 4.278 × 10^-3 substitutions per site per year (95% confidence interval [CI], 3.822 × 10^-3/site/year to 4.710 × 10^-3/site/year), and the common ancestor of E9 likely emerged around 1868 (95% CI, 1840 to 1892). The full-length genomic sequences of the fifteen E9 strains showed 76.9-79.6% nucleotide identity and 95.3-95.9% amino acid identity with E9 Barty strain. 11 of 15 E9 whole genome sequence present four recombination patterns, and E9 recombinants have extensive genetic exchanges in the 2C and P3 regions with other Enterovirus B (EV-B) circulated in China. Four of six E9 strains were temperature sensitive, and two were temperature resistant, and a comparative genomics analysis suggested that 411, 865 and 867 amino acid substitution in the P1 region was related to temperature sensitivity.

Conclusion: This study highlights a persistent transmission network of E9 in worldwide, provides valuable information regarding the molecular epidemiology of E9.