Genetic characterization and molecular epidemiological analysis of enterovirus C99 in China.

Enterovirus C99 (EV-C99) is a newly identified EV serotype within the species Enterovirus C. Few studies on EV-C99 have been conducted globally. More information and research on EV-C99 are needed to assess its genetic characteristics, phylogenetic relationships, and associations with enteroviral diseases. Here, the phylogenetic characteristics of 11 Chinese EV-C99 strains have been reported. The full-length genomic sequences of these 11 strains show 79.4-80.5% nucleotide identity and 91.7-94.3% amino acid (aa) identity with the prototype EV-C99. A maximum likelihood phylogenetic tree constructed based on the entire VP1 coding region identified 13 genotypes (A-M), revealing a high degree of variation among the EV-C99 strains. Phylogeographic analysis showed that the Xinjiang Uygur Autonomous Region is an important source of EV-C99 epidemics in various regions of China. Recombination analysis revealed inter-serotype recombination events of 16 Chinese EV-C99 strains in 5' untranslated regions and 3D regions, resulting in the formation of a single recombination form. Additionally, the Chinese strain of genotype J showed rich aa diversity in the P1 region, indicating that the genotype J of EV-C99 is still going through variable dynamic changes. This study contributes to the global understanding of the EV-C99 genome sequence and holds substantial implications for the surveillance of EV-C99.

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.

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.

Characterizing enterovirus C96 genome and phylodynamics analysis.

Enterovirus C96 (EV-C96) is a recently discovered serotype belonging to enterovirus C species. It had been isolated from patients with acute flaccid paralysis, hand, foot, and mouth disease, diarrhea, healthy people, or environmental specimens. Despite increasing reports of the virus, the small number of full-length genomes available for EV-C96 has limited molecular epidemiological studies. In this study, newly collected rare EV-C96 strains in China from 1997 to 2020 were combined with sequences available in GenBank for comprehensive analyses. Sequence analysis revealed that the nucleotide sequence similarity of EV-C96 and the prototype strain (BAN00-10488) was 75%-81.8% and the amino acid sequence similarity was 85%-94.9%. EV-C96 had a high degree of genetic variation and could be divided into 15 genogroups. The mean evolutionary rate was 5.16 × 10-3 substitution/site/year, and the most recent common ancestor was dated to 1925. A recombination analysis revealed that EV-C96 may be a recombinant derived from other serotypes in the EV-C group in the nonstructural protein coding region. This comprehensive and integrated analysis of the whole genome sequence of EV-C96 provides valuable data for further studies on the molecular epidemiology of EV-C96 worldwide.

Determination of the median effective dose of sufentanil for inhibiting the laryngeal mask insertion response in geriatric patients: a prospective, double-blinded, dose-response trial.

BACKGROUND: Laryngeal mask airway(LMA) have been widely used in clinical practice. Irritation to the patient during the insertion of a laryngeal mask can cause hemodynamic fluctuations, which is particularly unsafe for geriatric patients. We used probit regression analysis to determine the median effective dose of sufentanil to inhibit the response to LMA insertion in geriatric patients. METHODS: A total of 90 patients were selected for the study using the following inclusion criteria: age ≥ 65 years old, ASA grade I-III, and scheduled to undergo intravenous general anesthesia with LMA insertion. Each patient received a dose of sufentanil for anesthesia induction in one of six levels: 0.05, 0.1, 0.15, 0.2, 0.25, or 0.3 µg kg-1. LMA insertion was scored with a 3-point, 6-category scale, with scores ≥ 16 indicating effective LMA insertion, and < 16 indicating ineffective LMA insertion. Mean arterial blood pressure (MAP), heart rate (HR), and bispectral index (BIS) were recorded 1 min before induction (T1), 1 min after induction (T2), 1 min after LMA insertion (T3), and 5 min after LMA insertion (T4) in each group. In addition, the plasma norepinephrine (NE) levels and adverse reactions were measured at T2 and T3 in each dosage group. RESULTS: Probit regression analysis showed that the ED50 of sufentanil inhibiting the response to LMA insertion in geriatric patients was 0.18 µg kg-1 (95% CI: 0.16-0.21 µg kg-1), and the ED95 was 0.31 µg kg-1 (95% CI: 0.27-0.38 µg kg-1), and the probit(p) = -2.34 + 12.90 × ln(Dose)([Formula: see text] = 0.725, p = 0.948). Among all the patients, the number of effective LMA insertions was 57 (group A), and the number of ineffective LMA insertions was 33 (group B). The MAP, HR, and NE in group B were significantly higher than in group A at T3. CONCLUSIONS: Sufentanil can effectively inhibit the patient's response to LMA insertion, with stable hemodynamics and small stress response. The ED50 and ED95 were 0.18 µg kg-1 (95% CI: 0.16-0.21 µg kg-1) and 0.31 µg kg-1(95% CI: 0.27-0.38 µg kg-1), respectively. TRIAL REGISTRATION: This study was registered in the Chinese Clinical Trial Registry (Registration number: ChiCTR2100051827 ) on October 6, 2021.

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.

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.

Tunable 360° microwave photonic multichannel phase shifter with frequency quadrupling.

In this paper, a photonic scheme to generate a frequency-quadrupled microwave signal with full-range 360° tunable phase shift is proposed and experimentally demonstrated. Pure ±2nd-order sidebands with the carrier suppressed are generated by a dual-parallel Mach-Zehnder modulator. The two sidebands are separated by a fiber Bragg grating with the polarization state of one sideband rotated by 90 deg via a Faraday rotating mirror and then recombined to obtain a pair of orthogonally polarized wavelengths. The two orthogonally polarized optical sidebands are aligned into the same polarization direction by using a polarizer (Pol). Finally, a radio frequency (RF) signal with frequency quadrupling of a local oscillator signal is obtained by beating the ±2nd-order sidebands at a photodetector. The phase θ of the frequency-quadruped RF signal can be independently and arbitrarily adjusted from 0° to 360° through control of the polarization direction. Experiments are carried out to demonstrate the scheme, and frequency-quadrupled microwave signals at 12 GHz and 16 GHz are generated. A continuous phase shift from 0° to 360° of the frequency-quadrupled signal at 12 GHz is also verified. The proposed scheme not only has the merits of large operation bandwidth and simple structure but also can be extended to multi-channel applications.

The Median Effective Dose of Dexmedetomidine for the Inhibition of Emergence Delirium in Preschool Children Undergoing Tonsillectomy and/or Adenoidectomy: A Retrospective Dose-response Trial.

The incidence of emergence delirium (ED) is higher in preschool children undergoing tonsillectomy and/or adenoidectomy. The purpose of this study was to determine the median effective dose (ED50) of dexmedetomidine (DEX) for the inhibition of ED in preschool children by using probit regression analysis. A total of 140 anesthesia records were retrieved and divided into seven groups based on the infusion rate of DEX: .2, .25, .3, .35, .4, .45, and .5 µg·kg-1·h-1. The Pediatric Anesthesia Emergence Delirium Scale (PAEDS) was used to assess ED in preschool children, and ED was defined as a PAEDS score ≥ 10. Probit regression analysis revealed that the ED50 and ED95 of DEX were .31 µg·kg-1·h-1 (95% CI: .29-.35) and .48 µg·kg-1·h-1 (95% CI: .44-.56), respectively. Probit(p) = -2.84 + 9.28 × ln (Dose), (χ2 = 1.925, P = .859). The PAEDS score was significantly increased in the ED group, and the rate of bradycardia was significantly decreased in the ED group compared with the without ED group (27.3% vs 54.1%, P = .02). DEX can effectively inhibit the ED in preschool children undergoing tonsillectomy and/or adenoidectomy, however, bradycardia was the main complication.

Effects of glycine 64 substitutions in RNA-dependent RNA polymerase on ribavirin sensitivity and pathogenicity of coxsackievirus A6.

Hand, foot, and mouth disease (HFMD) caused by a group of enteroviruses is a global public health problem. In recent years, coxsackievirus A6 (CVA6) has emerged as an important HFMD agent. Previous studies have shown that mutations of glycine 64 in RNA-dependent RNA polymerase (3D polymerase), which is central to viral replication, cause phenotypic changes such as ribavirin resistance, increased replication fidelity, and virulence attenuation in poliovirus and enterovirus A71. In this study, we constructed CVA6 mutants with G64R, G64S, and G64T substitutions by site-directed mutagenesis in full-length cDNA of an infectious CVA6 strain cloned in pcDNA3.1. Viral RNA was obtained by in vitro transcription, and the rescued virus strains were propagated in RD cells. Sequencing after six passages revealed that G64S and G64T mutations were stably inherited, whereas G64R was genetically unstable and reversed to the wild type. Comparison of the biological characteristics of the wild-type and mutant CVA6 strains in an in vivo model (one-day-old ICR mice) revealed that the pathogenicity of CVA6-G64S and CVA6-G64T was significantly reduced compared to wild-type CVA6. In vitro experiments indicated the mutant CVA6-G64S and CVA6-G64T strains had increased resistance to 0.8 mM ribavirin and a decreased replication rate in the presence of 0.8 mM guanidine hydrochloride. Our results show that mutation of residue 64 reduces CVA6 susceptibility to ribavirin and increases CVA6 susceptibility to guanidine hydrochloride, together with increased replication fidelity and attenuated viral pathogenicity, thus laying a foundation for the development of safe and effective live attenuated CVA6 vaccine.

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.

Evolutionary Diversity of Coxsackievirus A6 Causing Severe Hand, Foot, and Mouth Disease - China, 2012-2023.

Introduction: Coxsackievirus A6 (CVA6) has emerged as a significant pathogen responsible for severe cases of hand, foot, and mouth disease (HFMD). This study aims to delineate the demographic characteristics and analyze the viral evolution of severe HFMD associated with CVA6, thereby assisting in its surveillance and management. Methods: In this investigation, 74 strains of CVA6 were isolated from samples collected from severe HFMD cases between 2012 and 2023. The VP1 gene sequences of CVA6 were amplified and analyzed to assess population historical dynamics and evolutionary characteristics using BEAST, DnaSP6, and PopART. Results: A significant portion (94.4%) of severe CVA6-associated HFMD cases (51 out of 54, with 20 lacking age information) were children under 5 years old. Among the 74 CVA6 strains analyzed, 72 belonged to the D3a sub-genotype, while only two strains were D2 sub-genotype. The average genetic distance between VP1 sequences prior to 2015 was 0.027, which increased to 0.051 when compared to sequences post-2015. Historical population dynamics analysis indicated three significant population expansions of severe CVA6-associated HFMD during 2012-2013, 2013-2014, and 2019-2020, resulting in the formation of 65 distinct haplotypes. Consistent with the MCC tree findings, transitioning between regional haplotypes required multiple base substitutions, showcasing an increase in population diversity during the evolutionary process (from 14 haplotypes in 2013 to 55 haplotypes over the subsequent decade). Conclusions: CVA6, associated with severe HFMD, is evolving and presents a risk of outbreak occurrence. Thus, enhanced surveillance of severe HFMD is imperative.

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.

Evolution, recombination and geographic spreading of global Coxsackievirus A6.

BACKGROUND: The increasing incidence of hand, foot, and mouth disease (HFMD) associated with Coxsackievirus A6 (CVA6) has become a very significant public health problem. The aim of this study is to investigate the recombination, geographic transmission, and evolutionary characteristics of the global CVA6. METHODS: From 2019 to 2022, 73 full-length CVA6 sequences were obtained from HFMD patients in China and analyzed in combination with 1032 published whole genome sequences. Based on this dataset, the phylogenetic features, recombinant diversity, Bayesian phylodynamic characteristics, and key amino acid variations in CVA6 were analyzed. RESULTS: The four genotypes of CVA6, A, D, E, and F, are divided into 24 recombinant forms (RFs, RF-A - RF-X) based on differences in the P3 coding region. The eastern China region plays a key role in the dissemination of CVA6 in China. VP1-137 and VP1-138 are located in the DE loop on the surface of the CVA6 VP1 protein, with the former being a highly variable site and the latter having more non-synonymous substitutions. CONCLUSIONS: Based on whole genome sequences, this study contributes to the CVA6 monitoring, early warning, and the pathogenic mechanism by studying recombination diversity, geographical transmission characteristics, and the variation of important amino acid sites.

Transmission of recombinant enterovirus A76 (EV-A76) in Xinjiang Uighur autonomous region of China.

Enterovirus 76 (EV-A76) is a serotype of enterovirus A and has been rarely reported. In this paper, we present the genetic characteristics of 15 EV-A76 isolates reported to circulate in the Xinjiang Uighur autonomous region of China in 2011. Sequence analysis revealed that all Chinese EV-A76 isolates had high similarity (> 98.3%) in the VP1 region, and five Chinese EV-A76 isolates were selected for whole genome sequencing based on VP1 nucleotide divergence. Similarity plots and boot-scanning analyses revealed frequent intertypic recombination in the nonstructural region of the EV-A76 isolate, as found with the EV-A89 donor sequence (also isolated in Xinjiang). The breakpoint of recombination is around nucleotide 3960, and the recombinant fragments covered part 2C and all P3 regions. This study increases publicly available EV-A76 nucleotide sequence and further our understanding EV-A76 molecular epidemiology.

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.

Epidemiology of Hand, Foot, and Mouth Disease and Genetic Evolutionary Characteristics of Coxsackievirus A10 in Taiyuan City, Shanxi Province from 2016 to 2020.

In recent years, the prevalence of hand, foot, and mouth disease (HFMD) caused by enteroviruses other than enterovirus A71 (EV-A71) and coxsackievirus A16 (CVA16) has gradually increased. The throat swab specimens of 2701 HFMD cases were tested, the VP1 regions of CVA10 RNA were amplified using RT-PCR, and phylogenetic analysis of CVA10 was performed. Children aged 1-5 years accounted for the majority (81.65%) and boys were more than girls. The positivity rates of EV-A71, CVA16, and other EVs were 15.22% (219/1439), 28.77% (414/1439), and 56.01% (806/1439), respectively. CVA10 is one of the important viruses of other EVs. A total of 52 CVA10 strains were used for phylogenetic analysis based on the VP1 region, 31 were from this study, and 21 were downloaded from GenBank. All CVA10 sequences could be assigned to seven genotypes (A, B, C, D, E, F, and G), and genotype C was further divided into C1 and C2 subtypes, only one belonged to subtype C1 and the remaining 30 belonged to C2 in this study. This study emphasized the importance of strengthening the surveillance of HFMD to understand the mechanisms of pathogen variation and evolution, and to provide a scientific basis for HFMD prevention, control, and vaccine development.

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.

Genetic characterization and molecular evolution of type 3 vaccine-derived polioviruses from an immunodeficient patient in China.

In 2013, a case of immunodeficiency vaccine-derived poliovirus (iVDPV) was identified in Jiangxi Province, China. In this study, we purified 14 type 3 original viral isolates from this case and characterized the molecular evolution of these iVDPVs for 298 days. Genetic variants were found in most of the original viral isolates, with complex genetic and evolutionary relationships among the variants. A phylogenetic tree constructed based on the P1 region showed that these iVDPVs were classified into lineage A and B. The dominant lineage B represents a major trend in virus evolution. The nucleotide substitution rate at the third codon position (3CP) estimated by the BEAST program was 1.76 × 10-2 substitutions/site/year (95% HPD: 1.23-2.39 × 10-2). The initial OPV dose was given dating back to March 2013, which was close to the time of the last OPV vaccination, suggesting that OPV infection may have originated with the last dose of vaccine. Recombinant analysis showed that these iVDPVs were inter-vaccine recombinants with two recombination patterns, S3/S2/S1 and S3/S2/S3/S2/S1. Whole genome sequence analysis revealed that key nucleotide sites (C472U, C2034U, U2493C) associated with the attenuated phenotype of Sabin 3 have been replaced. Temperature sensitivity test showed that all tested strains were temperature-sensitive, except for the variant Day11-5. Interestingly, we observed that the variant Day11-5 temperature resistance properties may be associated with the Lys to Met substitution at the VP2-162 site. Serological test and whole genome sequence analysis showed that the seropositivity rate remained high, and mutations in the antigenic sites did not significantly alter neutralization ability.

Genomic and evolutionary characteristics of G3P[8] group a rotavirus strains in China, 2016 to 2018.

Rotavirus A (RVA) G3P[8] is sporadically detected in China, although G9P[8] predominates. To evaluate their genetic composition at the whole-genome level, 24 G3P[8] RVA strains isolated from children under five years were sequenced and characterized. The 24 strains were genotyped as G3-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1, indicating the Wa-like genotype constellation. A maximum clade credibility (MCC) tree for VP7 indicated that G3 had an estimated mean evolutionary rate of 7.279 × 10-4 substitutions/site/year; thus, 3-5 years would pass from the generation of an ancestor virus to the epidemic spread of that virus throughout China. Considering the ongoing prevalence as well as rapid evolution, it is important to monitor G3P[8] RVA epidemics; continuous nationwide surveillance is essential.