Coronavirus disease 2019 (COVID-19) has become a pandemic with increasing numbers of cases worldwide. SARS-CoV-2, the causative virus of COVID-19, is mainly transmitted through respiratory droplets or through direct and indirect contact with an infected person. The possibility of potential faecal-oral transmission was investigated in this study. We collected 258 faecal specimens from nine provinces in China and detected the nucleic acid of SARS-CoV-2 using real-time RT–PCR. Vero cells were used to isolate the virus from SARS-CoV-2 nucleic acid positive samples, after which sequencing of Spike gene in eight samples was performed. In all, 93 of 258 (36%) stool samples were positive for SARS-CoV-2 RNA. The positive rates of critical, severe, moderate, and mild patients were 54.4%, 56.1%, 30.8%, and 33.3%, respectively. The content of nucleic acid increased within 2 weeks after the onset of the disease. From the perspective of clinical typing, the nucleic acid can be detected in the faeces of critical patients within two weeks and until four to five weeks in the faeces of severe and mild patients. SARS-CoV-2 was isolated from stool specimens of two severe patients. Four non-synonymous mutations in Spike gene were newly detected in three stool samples. A small number of patients had strong faecal detoxification ability. The live virus in faeces could be an important source of contamination, which may lead to infection and further spread in areas with poor sanitary conditions. The findings of this study have public health significance and they should be considered when formulating disease control strategies.

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.

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⁻³ substitutions site^-1 year⁻¹ 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.

Coxsackievirus A6 (CVA6), a member of species A enterovirus, is associated with outbreaks of hand-foot-and-mouth disease and causes a large nationwide burden of disease. However, the molecular pathogenesis of CVA6 remains unclear. In the present study, we established a suckling Institute of Cancer Research (ICR) mouse infection model to explore the neural pathogenicity of CVA6. Five-day-old mice infected with CVA6 strain F219 showed lethargy and paralysis, and died 5 or 6 days after infection via IM injection. Cerebral edema and neuronal cell swelling were observed in the infected brain tissue, and we found that the CVA6 VP1 antigen could co-localize with GFAP-positive astrocytes in infected mouse brain using an immunofluorescence assay. CVA6 strain F219 can also infect human glioma (U251) cells. Transcriptome analysis of brain tissues from infected mice and infected U251 cells showed that significantly differentially expressed genes were enriched in antiviral and immune response and neurological system processes. These results indicate that CVA6 could cause neural pathogenesis and provide basic data for exploring the mechanism of how host–cell interactions affect viral replication and pathogenesis. Importance: Coxsackievirus A6 (CVA6) surpasses the two main pathogens, enterovirus 71 (EV-A71) and coxsackievirus A16 (CVA16), which are the leading pathogens causing HFMD in many provinces of China. In our study, CVA6 infection caused neurogenic pathogenesis in a neonatal murine model, manifesting as cerebral edema and neuronal cell swelling, CVA6 VP1 antigen could co-localize with GFAP-positive astrocytes in the infected mouse brain. Based on CVA6-infected brain tissue and U251 cell transcriptome analysis, we found upregulated antiviral and immune response-related genes such as Zbp1, Usp18, Oas2, Irf7, Ddx60, Ifit3, Ddx58, and Isg15, while the neurological system process-related genes were downregulated, including Fcrls, Ebnrb, Cdk1, and Anxa5.

As the coronavirus disease 2019 (COVID‐19) pandemic is still ongoing and severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) variants are circulating worldwide, an increasing number of breakthrough infections are being detected despite the good efficacy of COVID‐19 vaccines. Data on 88 COVID‐19 breakthrough cases (breakthrough infections group) and 41 unvaccinated cases (unvaccinated group) from June 1 to August 22, 2021, were extracted from a cloud database established at Beijing Ditan Hospital to evaluate the clinical, immunological, and genomic characteristics of COVID‐19 breakthrough infections. Among these 129 COVID‐19 cases, 33 whole genomes were successfully sequenced, of which 23 were Delta variants, including 15 from the breakthrough infections group. Asymptomatic and mild cases predominated in both groups, but two patients developed severe disease in the unvaccinated group. The median time of viral shedding in the breakthrough infections group was significantly lower than that in the unvaccinated group (p = 0.003). In the breakthrough infections group, the IgG titers showed a significantly increasing trend (p = 0.007), and the CD4 + T lymphocyte count was significantly elevated (p = 0.018). For people infected with the Delta variant in the two groups, no significant difference was observed in either the quantitative reverse‐transcription polymerase chain reaction results or viral shedding time. In conclusion, among vaccinated patients, the cases of COVID‐19 vaccine breakthrough infections were mainly asymptomatic and mild, IgG titers were significantly increased and rose rapidly, and the viral shedding time was shorter.

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.