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.

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.

Importation and circulation of rubella virus lineages 1E-L2 and 2B-L2c between 2018 and 2021 in China: Virus evolution and spatial-temporal transmission characteristics.

To better understand the importation and circulation patterns of rubella virus lineages 1E-L2 and 2B-L2c circulating in China since 2018, 3,312 viral strains collected from 27 out of 31 provinces in China between 2018 and 2021 were sequenced and analyzed with the representative international strains of lineages 1E-L2 and 2B-L2c based on genotyping region. Time-scale phylogenetic analysis revealed that the global lineages 1E-L2 and 2B-L2c presented distinct evolutionary patterns. Lineage 1E-L2 circulated in relatively limited geographical areas (mainly Asia) and showed geographical and temporal clustering, while lineage 2B-L2c strains circulated widely throughout the world and exhibited a complicated topology with several independently evolved branches. Furthermore, both lineages showed extensive international transmission activities, and phylogeographic inference provided evidence that lineage 1E-L2 strains circulating in China possibly originated from Japan, while the source of lineage 2B-L2c isolated since 2018 is still unclear. After importation into China in 2018, the spread of lineage 1E-L2 presented a three-stage transmission pattern from southern to northern China, whereas lineage 2B-L2c spread from a single point in western China to all the other four regions. These two transmission patterns allowed both imported lineages to spread rapidly across China during the 2018-9 rubella epidemic and eventually established endemic circulations. This study provides critical scientific data for rubella control and elimination in China and worldwide.

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 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.

Highly diverse ribonucleic acid viruses in the viromes of eukaryotic host species in Yunnan province, China.

Background: The diversity in currently documented viruses and their morphological characteristics indicates the need for understanding the evolutionary characteristics of viruses. Notably, further studies are needed to obtain a comprehensive landscape of virome, the virome of host species in Yunnan province, China. Materials and methods: We implemented the metagenomic next-generation sequencing strategy to investigate the viral diversity, which involved in 465 specimens collected from bats, pangolins, monkeys, and other species. The diverse RNA viruses were analyzed, especially focusing on the genome organization, genetic divergence and phylogenetic relationships. Results: In this study, we investigated the viral composition of eight libraries from bats, pangolins, monkeys, and other species, and found several diverse RNA viruses, including the Alphacoronavirus from bat specimens. By characterizing the genome organization, genetic divergence, and phylogenetic relationships, we identified five Alphacoronavirus strains, which shared phylogenetic association with Bat-CoV-HKU8-related strains. The pestivirus-like virus related to recently identified Dongyang pangolin virus (DYPV) strains from dead pangolin specimens, suggesting that these viruses are evolving. Some genomes showed higher divergence from known species (e.g., calicivirus CS9-Cali-YN-CHN-2020), and many showed evidence of recombination events with unknown or known strains (e.g., mamastroviruses BF2-astro-YN-CHN-2020 and EV-A122 AKM5-YN-CHN-2020). The newly identified viruses showed extensive changes and could be assigned as new species, or even genus (e.g., calicivirus CS9-Cali-YN-CHN-2020 and iflavirus Ifla-YN-CHN-2020). Moreover, we identified several highly divergent RNA viruses and estimated their evolutionary characteristics among different hosts, providing data for further examination of their evolutionary dynamics. Conclusion: Overall, our study emphasizes the close association between emerging viruses and infectious diseases, and the need for more comprehensive surveys.

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.

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.

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.

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.

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.

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.

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.

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.

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.