Search

Hand, foot, and mouth disease (HFMD), a common childhood infection caused by enterovirus, poses a serious public health concern in China. We collected and analyzed epidemiological data on 62,133 HFMD cases in Shenyang City, Liaoning Province, from 2013 to 2023. The average annual incidence was 76.12 per 100,000 person-years; 99.45% of cases were mild, while 0.55% were severe. Only one patient died. HFMD infections peaked annually in July. Children in kindergartens and scattered children accounted for 44.6% and 42.2% of cases, respectively. Real-time RT-PCR detection of enteroviruses in 5534 patient samples revealed 3780 positives, of which 25.1% were CVA16-positive. Positives were randomly sampled, yielding 240 VP1 sequences of CVA16. Phylogenetic tree results showed that all VP1 sequences belonged to the B1 sub-genogroup. However, the sub-genogroup prevalence varied over time: from 2013 to 2014 and 2019 to 2021, the predominant sub-genogroup was B1a, while it was B1b from 2015 to 2018. Further phylogenetic analyses showed substantial divergence between B1a branches in CVA16, suggesting possible turnover of the B1a sub-genogroup in CVA16 due to evolution. This study provides epidemiological data on HFMD in Shenyang, and provides a phylogenetic analysis of CVA16, offering a theoretical basis for preventing and controlling HFMD in Shenyang City. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Coxsackievirus B5 (CVB5) is an important enterovirus B species (EV-Bs) type. We used the full-length genomic sequences of 53 viral sequences from the national hand, foot, and mouth disease surveillance network in the Chinese mainland (2001–2021). Among them, 69 entire VP1 coding region nucleotide sequences were used for CVB5 genotyping and genetic evolution analysis. Phylogenetic analysis based on a data set of 448 complete VP1 sequences showed that CVB5 could be divided into four genotypes (A-D) worldwide. Sequences from this study belonged to genotypes B and D, which dominated transmission in the Chinese mainland. Two transmission lineages of CVB5 have been discovered in the Chinese mainland, lineage 2 was predominant. Markov chain Monte Carlo analysis indicated that the tMRCA of CVB5 in the Chinese mainland could be traced to 1955, while the global trend could be traced to 1862, 93 years earlier than China. The evolution rate of CVB5 was higher in the Chinese mainland than worldwide. The spatiotemporal dynamics analysis of CVB5 assessed that virus transportation events were relatively active in Central, Northeast, North and Northwest China. Recombination analysis revealed frequent intertypic recombination in the non-structural region of CVB5 genotypes B and D with the other EV-Bs, revealing eight recombination lineages. Our study showed the molecular evolution and phylogeography of CVB5 that could provide valuable information for disease prevention. [ABSTRACT FROM AUTHOR]

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)( χ 2 = 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. [ABSTRACT FROM AUTHOR]

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)( χ 2 = 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Background China implemented the globally synchronized switch from trivalent oral poliovirus vaccine (tOPV) to bivalent OPV (bOPV) on May 1, 2016. During April 2018 to May 2019, the first outbreak caused by type 2 circulating vaccine-derived poliovirus (cVDPV2) after the switch occurred in Xinjiang and Sichuan, China. Methods. We performed sequence analysis of VP1 and the whole genome to determine the genomic characteristics of type 2 cVDPVs, and carried out coverage surveys to assess the risk of viral propagation. Surveillance for environment and acute flaccid paralysis was intensified to enhance case ascertainment. Results. Comparison of the complete genomes between early (Xinjiang strain) and late strains (Sichuan strains) revealed that recombination pattern and reverse mutation of attenuation sites had been fixed early, but the mutations of the neutralizing antigenic sites were introduced over the circulation. The Markov Chain Monte Carlo tree showed that the cVDPV2 initial infection was April 2016, earlier than the switch. So, we speculated that the cVDPV2 was originated from tOPV recipients and spread among children with a low level of immunity against the type 2. Conclusions The detection of this outbreak combined acute flaccid paralysis (AFP) surveillance with environmental surveillance (ES) indicates that ES should be expanded geographically to further complement AFP surveillance. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Hand, foot, and mouth disease (HFMD) is a common global epidemic. From 2008 onwards, many HFMD outbreaks caused by coxsackievirus A6 (CV-A6) have been reported worldwide. Since 2013, with a dramatically increasing number of CV-A6-related HFMD cases, CV-A6 has become the predominant HFMD pathogen in mainland China. Phylogenetic analysis based on the VP1 capsid gene revealed that subtype D3 dominated the CV-A6 outbreaks. Here, we performed a large-scale (near) full-length genetic analysis of global and Chinese CV-A6 variants, including 158 newly sequenced samples collected extensively in mainland China between 2010 and 2018. During the global transmission of subtype D3 of CV-A6, the noncapsid gene continued recombining, giving rise to a series of viable recombinant hybrids designated evolutionary lineages, and each lineage displayed internal consistency in both genetic and epidemiological features. The emergence of lineage –A since 2005 has triggered CV-A6 outbreaks worldwide, with a rate of evolution estimated at 4.17 × 10−3 substitutions site-1 year−1 based on a large number of monophyletic open reading frame sequences, and created a series of lineages chronologically through varied noncapsid recombination events. In mainland China, lineage –A has generated another two novel widespread lineages (–J and –L) through recombination within the enterovirus A gene pool, with robust estimates of occurrence time. Lineage –A, –J, and –L infections presented dissimilar clinical manifestations, indicating that the conservation of the CV-A6 capsid gene resulted in high transmissibility, but the lineage-specific noncapsid gene might influence pathogenicity. Potentially important amino acid substitutions were further predicted among CV-A6 variants. The evolutionary phenomenon of noncapsid polymorphism within the same subtype observed in CV-A6 was uncommon in other leading HFMD pathogens; such frequent recombination happened in fast-spreading CV-A6, indicating that the recovery of deleterious genomes may still be ongoing within CV-A6 quasispecies. CV-A6-related HFMD outbreaks have caused a significant public health burden and pose a great threat to children's health; therefore, further surveillance is greatly needed to understand the full genetic diversity of CV-A6 in mainland China. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Coxsackievirus A8 (CV-A8) is one of the pathogens associated with hand, foot and mouth disease (HFMD) and herpangina (HA), occasionally leading to severe neurological disorders such as acute flaccid paralysis (AFP). Only one study aimed at CV-A8 has been published to date, and only 12 whole-genome sequences are publicly available. In this study, complete genome sequences from 11 CV-A8 strains isolated from HFMD patients in extensive regions from China between 2013 and 2018 were determined, and all sequences from GenBank were retrieved. A phylogenetic analysis based on a total of 34 complete VP1 sequences of CV-A8 revealed five genotypes: A, B, C, D and E. The newly emerging genotype E presented a highly phylogenetic divergence compared with the other genotypes and was composed of the majority of the strains sequenced in this study. Markov chain Monte Carlo (MCMC) analysis revealed that genotype E has been evolving for nearly a century and somehow arose in approximately 2010. The Bayesian skyline plot showed that the population size of CV-A8 has experienced three dynamic fluctuations since 2001. Amino acid residues of VP1100N, 103Y, 240T and 241V, which were embedded in the potential capsid loops of genotype E, might enhance genotype E adaption to the human hosts. The CV-A8 whole genomes displayed significant intra-genotypic genetic diversity in the non-capsid region, and a total of six recombinant lineages were detected. The Chinese viruses from genotype E might have emerged recently from recombining with European CV-A6 strains. CV-A8 is a less important HFMD pathogen, and the capsid gene diversity and non-capsid recombination variety observed in CV-A8 strains indicated that the constant generation of deleterious genomes and a constant selection pressure against these deleterious mutations is still ongoing within CV-A8 quasispecies. It is possible that CV-A8 could become an important pathogen in the HFMD spectrum in the future. Further surveillance of CV-A8 is greatly needed. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

When facing massive statistical data, the k-means algorithm is very difficult to satisfy the need of data processing as it lacks an effective parallel mechanism. This paper proposes an improved k-means algorithm (IMR-KCA) to conduct clustering analysis based on medical data employing MapReduce computing framework. Through analyzing the defects of vast redundancy in the traditional k-means algorithms, a selection model is firstly proposed to simplify the computations with multiple clustering centers. Based on several proposed theorems, we prove the correctness of this selection model. Second, this paper provides a method to calculate the distances from extreme points to central points, and the original Euclidean distance is replaced with Manhattan distance. For this simplification, a group of theorems are proposed to prove the correctness. Next, we provide a group of implementation algorithms to complete the parallelism of the clustering computation employing the MapReduce framework. Finally, the experimental results illustrate that IMR-KCA is more reliable and efficient than the direct parallelization of the traditional clustering algorithms based on MapReduce. Copyright © 2017 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]