Next Article in Journal
Plant Defense and Viral Counter-Defense during Plant–Geminivirus Interactions
Next Article in Special Issue
Genomic Epidemiology and Transmission Dynamics of Global Coxsackievirus B4
Previous Article in Journal
Efficacy Validation of SARS-CoV-2-Inactivation and Viral Genome Stability in Saliva by a Guanidine Hydrochloride and Surfactant-Based Virus Lysis/Transport Buffer
Previous Article in Special Issue
Evolutionary and Genetic Recombination Analyses of Coxsackievirus A6 Variants Associated with Hand, Foot, and Mouth Disease Outbreaks in Thailand between 2019 and 2022
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Coxsackievirus A6 Infection Causes Neurogenic Pathogenesis in a Neonatal Murine Model

1
WHO WPRO Regional Polio Reference Laboratory, National Health Commission Key Laboratory for Biosafety, National Health Commission Key Laboratory for Medical Virology, National Institute for Viral Disease Control and Prevention, Chinese Center for Disease Control and Prevention, 155 Changbai Road, Beijing 102206, China
2
Center for Biosafety Mega-Science, Chinese Academy of Sciences, Wuhan 430071, China
*
Author to whom correspondence should be addressed.
These authors have contributed equally to this work.
Viruses 2023, 15(2), 511; https://doi.org/10.3390/v15020511
Submission received: 4 January 2023 / Revised: 3 February 2023 / Accepted: 9 February 2023 / Published: 12 February 2023
(This article belongs to the Special Issue Coxsackieviruses and Associated Diseases)

Abstract

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

1. Introduction

Enteroviruses (EVs) are found commonly worldwide and cause various diseases such as hand, foot, and mouth disease (HFMD), viral encephalitis, aseptic meningitis, and acute flaccid paralysis [1]. It is worth noting that out of the annual reported HFMD cases in mainland China since 2013, the proportion caused by coxsackievirus A6 (CVA6) has increased every year, and at present, has surpassed enterovirus 71 (EV-A71) and coxsackievirus A16 (CVA16), the leading pathogens in many provinces of China [2,3]. CVA6-related HFMD presents several clinical features, such as atypical herpes and nail removal, which are sometimes complicated by severe clinical symptoms such as meningitis and encephalitis [4,5,6,7]. However, no vaccine or effective drugs are present [8], and the pathogenesis of CVA6-associated diseases such as HFMD, meningitis, and encephalitis is unclear [9].
Several previous studies have reported that CVA6 murine models using one- to ten-day-old BALB/c or Institute of Cancer Research (ICR) mouse strains through IP, IM, and IC injection exhibited neurotropism and triggered systemic manifestations [8,9,10,11]. However, changes in biological processes at the tissue level and in their regulation at the molecular level during CVA6 infection are rare. Several studies have analyzed the transcriptomic expression profiles of EV infections, mainly EV-A71 and CVA16, which provide clues regarding the pathogenic mechanisms.
In the present study, we established a neonatal mouse model of CVA6 infection for CVA6-related neurological disease research. In addition, we used this model to describe the neurogenic pathogenesis of CVA6 and describe the dysregulation at the molecular level in mouse brain tissue and human glioma (U251) cell level via RNA sequencing (RNA-Seq). Our study established a 5-day-old CVA6-infected ICR mouse model, and provided basic data for exploring the mechanism of host–cell interaction affecting viral neurogenic pathogenesis.

2. Results

2.1. Establishment of a Neonatal ICR Mouse Model of CVA6 Infection

In this study, to establish a mouse model for studying CVA6 pathogenesis, we selected 5-day-old ICR mice through IC, IM, or IP injection, to inoculate CVA6 with serially diluted 10-fold (103 to 106 50% tissue culture infective dose [TCID50] per mouse) (n = 7 per group) and one negative group (n = 7 per group). Their clinical symptoms were observed, and the differences among the three routes of infection were compared.
The average body weights of the mice in the experimental groups infected with 105–106 TCID50 of CVA6 were significantly different from those of the mice in the 103–104 TCID50 of CVA6 and the negative control groups via all three injection routes (Figure 1A–C). All mice infected with 105 to 106 TCID50 via the IM route died 5 to 6 days post infection (dpi); however, under these two doses, some infected mice died after IC and IP routes. At the 104 TCID50 infection dose level, some infected mice died after IM and IC routes, and the mortality of the former was higher than that of the latter; nevertheless, all infected mice survived via the IP route. At the 103 TCID50 infection dose level, comparing all the infected mice that survived via IC and IP route, only small partially infected mice died via IM route (Figure 1D–F). Moreover, the average clinical scores tended to be consistent with average body weight and survival rates. All mice infected with 105–106 TCID50 via the IM route developed significant hindlimb paralysis, gradually lost weight, and died within 5–6 dpi, with a clinical score of 5 between 5–6 dpi (Figure 1G–I). These results indicated that mice infected with 105 to 106 TCID50 via the IM route were suitable for the infection mouse model with a stable time of disease onset, mortality rate, and good reproducibility. Therefore, in this study, we selected 105 TCID50 via the IM route for further investigation.
However, under these two doses, some infected mice infected by IC and IP routes died.

2.2. Effective Replication and Pathological Changes of CVA6 in Infected Mice Brain and Muscular Tissues

The viral loads in 5-day-old mice inoculated with 105 TCID50 via the IM route were measured at 1, 3, and 5 dpi. The results showed that the viral loads in the muscular tissues ranged from 103.5 TCID50 at 1 dpi to 106.25 TCID50 at 5 dpi, which was consistent with many other studies indicating the muscular tissues, the strongest tissue tropism of EVs (Figure 2A). However, the viral loads in the brain tissues ranged from 101.25 TCID50 at 1 dpi to 103.75 TCID50 at 5 dpi, which indicated effective replication of CVA6 in suckling mouse brain (Figure 2B).
Subsequently, these two infected mouse tissues at 5 dpi underwent hematoxylin and eosin (H&E) staining for analysis of histopathological changes. Compared with the control group, in infected mice, hind limb skeletal muscle fibers exhibited severe necrosis and rupture, as well as inflammatory edema (Figure 2C–E). The encephalitis symptoms present in infected mice with cerebral edema and neuronal cell swelling were similar to the neurological symptoms in humans (Figure 2F–H).

2.3. CVA6 Antigen Could Co-Localize with GFAP-Positive Astrocytes in Infected Mice Brain

The pathological manifestations of fatal cases caused by CVA6 are characterized by inflammatory damage to the central nervous system (CNS) [12]. Many studies have shown that astrocytes are emerging as pivotal regulators of CNS inflammatory responses [13,14]. GFAP (glial fibrillary acidic protein) is a class-III intermediate filament, a cell-specific marker, which differentiates astrocytes from other glial cells during the development of the central nervous system. Thus, in this study, we conducted immunofluorescence analysis of the potential targeted cells in the infected mouse brain [15]. Immunofluorescence staining studies showed that compared to NeuN-positive neuronal cells, the CVA6 VP1 antigen co-localized with GFAP-positive astrocytes (Figure 3A,B). This result indicates that the potential preferential infection of astrocytes by CVA6 may play a key role in viral neurological pathogenesis. Therefore, the exact molecular mechanisms should be further explored.

2.4. CVA6 Strain F219 Replicates Well in Human Glioma (U251) Cells

Since we found that the CVA6 VP1 antigen could co-localize with GFAP-positive astrocytes in the infected mouse brain, we further determined whether CVA6 could replicate in human neuronal cell lines. We measured viral replication dynamics in the human glioma cell line U251 and measured the viral titer at five time points after infection with a multiplicity of infection (MOI) of 5. CVA6 could cause cytopathic effects (CPE) after 24 h of infection and was replicated well in U251 cells via growth curves (Figure 4A,B). Consistent with the growth curve, immunofluorescence staining revealed that the CVA6 VP2 capsid protein could be detected in U251 cells and that it co-localized with GFAP (Figure 4C).

2.5. CVA6 Infection Increases the Defense Response to Virus Genes and Decreases the Number of Neurological System Genes in the Brain and U251 Cells at the Transcriptome Level

The above results demonstrated the dynamic virus titers and pathological changes at both the infected mouse brain tissue and human glioma (U251) cell levels. However, the molecular mechanisms underlying CVA6 pathogenesis are not yet fully understood. It is well known that viruses act as strictly parasitic organisms, and investigation of the mechanism of how host–cell interactions affect viral replication and pathogenesis is essential.
In this study, we described CVA6 infection, especially changes in differential host factor expression after infection, in human glioma (U251) cells and infected mouse brain using transcriptome sequencing analysis (Figure 5A). We used Hiplot Enhanced MA v0.1.0 analysis, which showed that the significantly differentially expressed genes (DEGs) were interferon-stimulated genes (ISGs), such as Ifit1, GBP4, etc., and inflammatory factors, such as CXCL10, Stat1, etc., which indicated that CVA6 maintained effective infections in human glioma (U251) cells and mouse brain (Figure 5B,C). The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEGs in the RNA-Seq of infected U251 cells showed that antiviral or inflammation-related signaling pathways, such as Type I interferon signaling pathway, toll-like receptor signaling pathway, cytokine–cytokine receptor interaction, and so on, and the neurological system, process pathway enrichments (Figure 5D,E). RNA-Seq of infected mouse brain showed antiviral or inflammation-related signaling pathways, such as the Type I interferon signaling pathway, RIG-I-like or Toll-like receptor signaling pathway, NF-kappa B signaling, and other neurological complications, such as neuroactive ligand–receptor interaction, neuron projection, glioma, etc. (Figure 5F,G). Notably through further visualization analysis of mutual DEGs, we found that some genes associated with defense response to viruses pathway (Oas2, Dtx3l, Gbp7, Ddx58, Usp18, Rtp4, Ifit3, Nlrc5, Rsad2, Oasl1, Irf7, Bst2, Trim34a, Ddx60, Isg15, Ifi204, Slfn8, Mx1, Stat1, Ifitm3, Parp9, Ifit1, Zbp1, and Itgam) were significantly upregulated, while most genes associated with neurological system process pathway (Rasgrp3, Abca9, Gpr37l1, Pla2g7, Fut9, Aqp4, Ednrb, Cdk1, Cenpf, Top2a, Hmgb2, S100b, Anxa5, Hebp5, and Sptssb) were significantly downregulated compared to the other three genes associated with neuro-inflammation pathways (Mertk, C1qb, and Fcrls) which were significantly upregulated (Figure 5H).
Taken together, we conclude that neurotropic CVA6 preferentially infects and replicates in astrocytes of the neonatal murine model. CVA6 could cause pathogenesis in the ICR mouse brain tissues and the human glioma (U251) cells to the detailed dysregulated genes with notable features showing upregulation of antiviral and immune response-related genes such as Zbp1, Usp18, Oas2, Irf7, Ddx60, Ifit3, Ddx58, and Isg15, while downregulating neurological system process-related genes, such as Fcrls, Ebnrb, Cdk1, and Anxa5.

3. Discussion

In recent years CVA6 has caused HFMD epidemics in mainland China have become a serious public health concern with a large number of severe cases [12]. Nevertheless, there is no effective vaccine or drug for CVA6 infection. Although many studies have reported the clinical medicine associated with central nervous system disease and epidemiological studies of CVA6 characterized by highly frequent recombination events correlated with pathogenicity, the molecular pathogenesis is deemed insufficient [2,6,13,14]. Animal models can be used to explore the pathogenic mechanisms of viral infection; since suckling murine and primate animals are suitable hosts for human enteroviruses, previous studies have shown that animal models can be used to explore possible mechanisms of enterovirus infection [15,16,17,18]. Several previous studies have reported CVA6 murine models of one- to ten-day-old BALB/c or ICR mouse strains. In this study, we established a 5-day-old CVA6-infected ICR mouse model and described the neurogenic pathogenesis. We conclude that neurotropic CVA6 preferentially infects and replicates in astrocytes of the neonatal mouse model, and also replicates well in human glioma (U251) cells. Finally, CVA6 infection increases the defense response to viral genes and reduced the number of neurological system genes in the brain and U251 cells at the transcriptome level.
CVA6 infection can cause neurological and systemic complications such as HFMD, aseptic meningitis, and encephalitis [9,10,11]. In this study, the brain edema, neuronal cell swelling, and high viral titers in the infected mice also led to significant astrocyte proliferation, microgliosis, and immune cell aggregation [19]. Previous studies have also shown that enterovirus 71 (EV-A71) preferentially infects astrocytes, and that echovirus 11 or 30 can infect human glioma (U251) cells, but not the neuroblastoma cell line SK-N-SH [20,21,22]. We also found CVA6 antigen co-localization with GFAP-positive astrocytes, which may be consistent with the enterovirus-infected brain leading to gliosis [23]. Further studies are needed to explore the consequences of neuronal death and dysfunction.
Previous studies have reported that EV-A71, CVA16, and echovirus 30 can result in the differential expression of genes in the nervous system development pathways [21,24,25,26]. In this study, CVA6 infected 5-day-old ICR mice by intramuscular injection and human glioma (U251) cells, and differentially expressed genes (DEGs) were screened using RNA-seq analysis. The results of this study indicate that after CVA6 infection, DEGs that are enriched in neurological system processes, such as Fcrls, Ebnrb, Cdk1, and Anxa5, which play important roles in glial cell growth and migration, were significantly downregulated. Astrocytes secrete key proteins and metabolites used to fuel neurons as well as in the biogenesis of neurotransmitters [27,28]. EV infection systematically causes neural inflammation and injury in the brain via dysfunctional astrocytes that disrupt brain homeostasis [29].
However, DEGs were enriched in antiviral and immune responses that were upregulated in response to CVA6 infection. Our analysis showed that Zbp1, Usp18, Oas2, Irf7, Ddx60, Ifit3, Ddx58, and Isg15 were significantly upregulated after CVA6 infection. CVA6 infection leads to the increased expression of ISGs, including Oas2, Irf7, Ddx60, Ifit3, Ddx58, and Isg15, which exhibit immunopathogenic potential and are implicated in neural inflammation [30]. Z-DNA-binding protein 1 (Zbp1) can be induced by type I interferons (IFNs) to upregulate IFN expression after viral infection, thus activating the expression of interferon regulatory factor 1 (IRF1) [31]. In addition, Zbp1 was activated by the receptor-interacting serine/threonine-protein kinase (RIPK) signaling pathway to alter neuronal cell metabolism and inhibit viral genome replication [32]. Previous studies have shown that ubiquitin specific peptidase 18 (Usp18) is involved in the negative regulation of IFN-induced inflammation and is enriched in the mitochondria of virus-infected cells to promote the polyubiquitination and aggregation of mitochondrial antiviral signaling protein (MAVS) [33]. This also provides insight into the specific mechanisms of CVA6 infection. Since the innate immune response is the first line of viral defense, our transcriptome sequencing data revealed dysregulation of neuropathy-related pathways [34]. Therefore, we used a five-day-old mouse model to clarify the pathological characteristics of CVA6 and found that antiviral and immune responses and neurological system process pathways could be involved in the pathogenesis of CVA6 at the molecular level. These findings provide a theoretical basis for understanding the pathogenic mechanisms.

4. Materials and Methods

4.1. Ethics Statement

ICR mice were purchased from SPF Biotechnology Co., Ltd. (Beijing, China). All animal experiments were approved by the Ethics Review Committee of the National Institute for Viral Disease Prevention and Control, Chinese Center for Disease Control and Prevention (approval No. 20201022059).

4.2. Virus, Cells, and Mice

One CVA6 clinical isolate, F219, was isolated from a severe case during an outbreak of HFMD in Yunnan in 2016. In addition to the typical clinical manifestations of HFMD, such as fever, rash on hands and feet, this patient also presented some nervous system symptoms, such as headache, poor spirit, lethargy, etc.
Human rhabdomyosarcoma (RD; ATCC) cells were cultured in Dulbecco’s modified Eagle medium (DMEM) (Gibco, CA, USA) containing 10% fetal bovine serum (FBS) (Gibco, USA) and 1% penicillin-streptomycin (HyClone, Logan, UT, USA) at 37 °C with 5% CO2. In addition, human glioma cells (U251; ATCC) were cultured in a 1:1 mixture of DMEM and F-12 (Gibco, California, USA), supplemented with 10% FBS. After the degree of cell fusion reached 80–90%, the medium was replaced with medium containing 2% FBS (maintenance media; MM) and incubated with the virus to obtain CPE. After three freeze-thaw cycles, all viruses in the supernatant were filtered through a 0.2-mm filter and stored at −80 °C.
All ICR mice were raised in independently ventilated cages and were provided with sufficient food and water.

4.3. Animal Infection Experiments

Five-day-old ICR mice were infected with CVA6 strain F219 at 103, 104, 105, or 106 50% TCID50 through IC, IM, and IP routes (n = 7 per group); control mice were inoculated with MM and reared separately from the infected mice. The weights, survival rates, and clinical scores of the mice were recorded daily. The grade of clinical disease was scored as follows: 0, no disease; 1, ruffled fur; 2, weight loss; 3, single limb paralysis; 4, hind limb paralysis; and 5, moribund, or dead [35].

4.4. Virus Isolation and Titers in Mouse Tissues

Brain and hind limb muscles were harvested from 5-day-old ICR mice via IM injection of F219 at 1, 3, and 5 dpi (n = 5 per time point). All samples were added to phosphate-buffered saline containing 1% penicillin-streptomycin and crushed in a tissue grinder (Scientz, Ningbo, China) for three freeze-thaw cycles before centrifugation to obtain the supernatant. The supernatant collected at each time point was serially diluted 10-fold and inoculated into RD cells that were cultured in 96-well plates. After 5 days of continuous observation, virus titers were detected using the TCID50 assay.

4.5. Growth Curves

U251 cells were infected at an MOI of 5 with the CVA6 for multicycle replication under the conditions of 5% CO2 and 37 °C. After 2 h of incubation, the cells were washed three times with PBS and the infection medium was added. At 2, 6, 12, 24, and 48 h post infection, supernatants were collected, and viral titers were determined in RD cells.

4.6. Histopathological and Immunofluorescence Assay

After IM injection of F219 in 5-day-old ICR mice, the brains and hind limb muscles of the experimental and control groups were harvested at 5 dpi. The tissues were fixed in formalin buffer for 24 h, dehydrated, embedded in paraffin, and sliced into 5-μm-thick sections. H&E staining was performed after deparaffinization of the tissue sections with xylene. Finally, all images were observed using an AxioCam MRc5 (Carl Zeiss, Berlin, Germany) at magnifications of 200 or 400×.
For immunofluorescence analysis, tissues were fixed in 10% paraformaldehyde solution for 24 h, embedded in optimal cutting temperature compound (OCT), and sectioned. OCT-embedded sections were stained to evaluate CVA6 VP1 co-localization with NeuN-positive neurons or GFAP-positive neuroglial cells. Rabbit polyclonal anti-VP1 CVA6 antibody (1:100 dilution; GeneTex, San Antonio, TX, USA, Cat No. GTX132346, USA), mouse monoclonal anti-GFAP antibody (1:100 dilution; Abcam, Cambridge, UK, Cat No. ab4648, United Kingdom), and mouse monoclonal anti-NeuN antibody (1:100 dilution; Abcam, Cat No. ab104224, United Kingdom) were incubated at 4 °C overnight.
The CVA6-infected U251 cells were detected using rabbit polyclonal anti-VP1 of the CVA6 antibody (1:100 dilution; GeneTex, Cat No. GTX132346, USA) and mouse monoclonal anti-GFAP antibody (1:100 dilution; Abcam, Cat No. ab4648, United Kingdom). Fluorescence Dylight 680, Goat Anti-Rabbit IgG and Dylight 488, and Goat Anti-Mouse IgG antibodies (1:500 dilution) were added and incubated at room temperature (20 °C) for 1 h. Finally, all images were observed using a Leica laser scanning confocal microscope (Leica, Wetzlar, Germany) at magnifications of 200× or 400×.

4.7. Library Preparation for Transcriptome Sequencing

CVA6 was infected via IM inoculation with 105 TCID50 in 5-day-old ICR mice and in the control group (n = 3 per group). TRIzol reagent was used to extract total RNA from the mouse brain tissues. Qubit 2.0 (Invitrogen, Waltham, MA, USA) was used to accurately quantify RNA concentrations and Bioanalyzer 2100 (Agilent, Santa Clara, CA, USA) was used to detect RNA integrity number (RIN ≥ 8). Briefly, the first strand of cDNA was synthesized using random primer reverse transcription and the second-strand cDNAs were amplified using RNase H and DNA polymerase. Finally, the enriched and purified cDNAs were sequenced on an Illumina NovaSeq 6000 platform. Gene expression levels were quantified using a feature count (version 1.5.0).

4.8. RNA-Seq Data Bioinformatics Analysis

The Deseq2 software was used for differential gene expression analysis, and a significance cut-off of p < 0.05 and Log2 fold change ≥ ±1.5 was used as the standard for screening differentially expressed genes. To clearly understand the biological significance of differentially expressed genes, GO and KEGG pathway analyses were performed. GO analysis was carried out for three aspects: cellular component (CC), biological process (BP), and molecular function (MF). KEGG pathway analyses were used to correlate the catalog of differentially expressed genes with system functions at higher levels of cells, species, and ecosystems to better understand the molecular response networks of protein-coding genes. These were then enriched and analyzed using the online tool Database for Annotation, Visualization, and Integrated Discovery (DAVID) Bioinformatics Resources 6.8.

4.9. Statistical Analysis

All statistical analyses were performed using GraphPad Prism 8.0 (GraphPad Software). The survival rates of different groups of mice were assessed using the log-rank (Mantel–Cox) test. The Wilcoxon test was used to analyze the clinical score curves. The tissue viral titers in different mice were analyzed by two-tailed Student’s t-tests, and p values < 0.05 were considered significant.

Author Contributions

Conceptualization, Y.Z. (Yong Zhang) and Q.S.; methodology, Q.S. and J.L.; software, Q.S.; validation, Q.S., J.L., R.W., T.S., Y.Z. (Yanjun Zong), C.W., Y.L., X.L., Y.S. and Y.Z. (Yong Zhang); formal analysis, Q.S.; investigation, Q.S.; resources, Y.Z. (Yong Zhang); data curation, Q.S. and J.L.; writing—original draft preparation, Q.S.; writing—review and editing, Y.Z. (Yong Zhang) and Q.S.; visualization, Q.S.; supervision, Y.Z. (Yong Zhang); project administration, Y.Z. (Yong Zhang); funding acquisition, Y.Z. (Yong Zhang). All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by National Key Research and Development Program of China (Project No. 2021YFC2302003).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Ethics Committee of National Institute for Viral Disease Control and Prevention, Chinese Center for Disease Control and Prevention (protocol code No.20201022059).

Informed Consent Statement

Not applicable.

Data Availability Statement

All sequencing data in this study are available in the Sequence Read Archive (SRA) database under the accession numbers PRJNA929099 and PRJNA929115.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lefkowitz, E.J.; Dempsey, D.M.; Hendrickson, R.C.; Orton, R.J.; Siddell, S.G.; Smith, D.B. Virus taxonomy: The database of the International Committee on Taxonomy of Viruses (ICTV). Nucleic Acids Res. 2018, 46, D708–D717. [Google Scholar] [CrossRef]
  2. Song, Y.; Zhang, Y.; Han, Z.; Xu, W.; Xiao, J.; Wang, X.; Wang, J.; Yang, J.; Yu, Q.; Yu, D.; et al. Genetic recombination in fast-spreading coxsackievirus A6 variants: A potential role in evolution and pathogenicity. Virus Evol. 2020, 6, veaa048. [Google Scholar] [CrossRef]
  3. Ji, T.; Han, T.; Tan, X.; Zhu, S.; Yan, D.; Yang, Q.; Song, Y.; Cui, A.; Zhang, Y.; Mao, N.; et al. Surveillance, epidemiology, and pathogen spectrum of hand, foot, and mouth disease in mainland of China from 2008 to 2017. Biosaf. Health 2019, 1, 32–40. [Google Scholar] [CrossRef]
  4. Hussain, K.; Chawla, S.; Muir, D.; Roberts, N. Clinical manifestations of coxsackievirus A6 what a dermatologist needs to know. Int. J. Dermatol. 2022. [Google Scholar] [CrossRef]
  5. Li, Y.; Xiong, T.; Meng, Y.; Zhao, R.; Yu, G.; Liu, J.; Xiao, J.; Wang, K.; Wu, T.; Wei, S.; et al. Risk factors for severe hand, foot, and mouth disease infected with Coxsackievirus A6: A hospital-based case-control study. J. Med. Virol. 2020, 92, 3144–3150. [Google Scholar] [CrossRef] [PubMed]
  6. Yang, X.; Li, Y.; Zhang, C.; Zhan, W.; Xie, J.; Hu, S.; Chai, H.; Liu, P.; Zhao, H.; Tang, B.; et al. Clinical features and phylogenetic analysis of severe hand-foot-and-mouth disease caused by Coxsackievirus A6. Infect. Genet. Evol. 2020, 77, 104054. [Google Scholar] [CrossRef]
  7. Xu, S.; Li, H.; Qiao, P.; Xu, G.; Zhao, D.; Lin, X.; Qin, Y.; Yu, H.; Zhang, X.; Zhang, W.; et al. Neonatal hand, foot, and mouth disease due to coxsackievirus A6 in Shanghai. BMC Pediatr. 2020, 20, 364. [Google Scholar] [CrossRef]
  8. Qian, S.S.; Wei, Z.N.; Jin, W.P.; Wu, J.; Zhou, Y.P.; Meng, S.L.; Guo, J.; Wang, Z.J.; Shen, S. Efficacy of a coxsackievirus A6 vaccine candidate in an actively immunized mouse model. Emerg. Microbes Infect. 2021, 10, 763–773. [Google Scholar] [CrossRef] [PubMed]
  9. Jiang, Z.; Zhang, Y.; Lin, H.; Cheng, Q.; Lu, X.; Liu, W.; Zhou, R.; Zhong, B.; Tian, X. A 10-Day-Old Murine Model of Coxsackievirus A6 Infection for the Evaluation of Vaccines and Antiviral Drugs. Front. Immunol. 2021, 12, 665197. [Google Scholar] [CrossRef] [PubMed]
  10. Zhang, Z.; Dong, Z.; Wei, Q.; Carr, M.J.; Li, J.; Ding, S.; Tong, Y.; Li, D.; Shi, W. A Neonatal Murine Model of Coxsackievirus A6 Infection for Evaluation of Antiviral and Vaccine Efficacy. J. Virol. 2017, 91, e02450–e02467. [Google Scholar] [CrossRef] [Green Version]
  11. Li, D.; Sun, T.; Tao, L.; Ji, W.; Zhu, P.; Liang, R.; Zhang, Y.; Chen, S.; Yang, H.; Jin, Y.; et al. A mouse-adapted CVA6 strain exhibits neurotropism and triggers systemic manifestations in a novel murine model. Emerg. Microbes Infect. 2022, 11, 2248–2263. [Google Scholar] [CrossRef]
  12. Zhao, T.S.; Du, J.; Sun, D.P.; Zhu, Q.R.; Chen, L.Y.; Ye, C.; Wang, S.; Liu, Y.Q.; Cui, F.; Lu, Q.B. A review and meta-analysis of the epidemiology and clinical presentation of coxsackievirus A6 causing hand-foot-mouth disease in China and global implications. Rev. Med. Virol. 2020, 30, e2087. [Google Scholar] [CrossRef]
  13. Wang, S.H.; Wang, A.; Liu, P.P.; Zhang, W.Y.; Du, J.; Xu, S.; Liu, G.C.; Zheng, B.S.; Huan, C.; Zhao, K.; et al. Divergent Pathogenic Properties of Circulating Coxsackievirus A6 Associated with Emerging Hand, Foot, and Mouth Disease. J. Virol. 2018, 92, e00303–e00318. [Google Scholar] [CrossRef] [PubMed]
  14. Chen, M.; Zuo, X.; Tan, Y.; Ju, Y.; Bi, F.; Wang, H.; Chen, M. Six amino acids of VP1 switch along with pandemic of CV-A6-associated HFMD in Guangxi, southern China, 2010–2017. J. Infect. 2019, 78, 323–337. [Google Scholar] [CrossRef]
  15. Liu, L.; Zhao, H.; Zhang, Y.; Wang, J.; Che, Y.; Dong, C.; Zhang, X.; Na, R.; Shi, H.; Jiang, L.; et al. Neonatal rhesus monkey is a potential animal model for studying pathogenesis of EV71 infection. Virology 2011, 412, 91–100. [Google Scholar] [CrossRef]
  16. Nagata, N.; Shimizu, H.; Ami, Y.; Tano, Y.; Harashima, A.; Suzaki, Y.; Sato, Y.; Miyamura, T.; Sata, T.; Iwasaki, T. Pyramidal and extrapyramidal involvement in experimental infection of cynomolgus monkeys with enterovirus 71. J. Med. Virol. 2002, 67, 207–216. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, S.Y.; Lin, T.L.; Chen, H.Y.; Lin, T.S. Early and rapid detection of enterovirus 71 infection by an IgM-capture ELISA. J. Virol. Methods 2004, 119, 37–43. [Google Scholar] [CrossRef]
  18. Nagata, N.; Iwasaki, T.; Ami, Y.; Tano, Y.; Harashima, A.; Suzaki, Y.; Sato, Y.; Hasegawa, H.; Sata, T.; Miyamura, T.; et al. Differential localization of neurons susceptible to enterovirus 71 and poliovirus type 1 in the central nervous system of cynomolgus monkeys after intravenous inoculation. J. Gen. Virol. 2004, 85, 2981–2989. [Google Scholar] [CrossRef] [PubMed]
  19. Matschke, J.; Lutgehetmann, M.; Hagel, C.; Sperhake, J.P.; Schroder, A.S.; Edler, C.; Mushumba, H.; Fitzek, A.; Allweiss, L.; Dandri, M.; et al. Neuropathology of patients with COVID-19 in Germany: A post-mortem case series. Lancet Neurol. 2020, 19, 919–929. [Google Scholar] [CrossRef]
  20. Zhang, G.; Li, J.; Sun, Q.; Zhang, K.; Xu, W.; Zhang, Y.; Wu, G. Pathological Features of Echovirus-11-Associated Brain Damage in Mice Based on RNA-Seq Analysis. Viruses 2021, 13, 2477. [Google Scholar] [CrossRef]
  21. Li, J.; Zhang, G.; Sun, Q.; Zhang, K.; Lu, H.; Xiao, J.; Han, Z.; Zhao, H.; Xu, W.; Zhang, Y.; et al. Pathological Characteristics of Echovirus 30 Infection in a Mouse Model. J. Virol. 2022, 96, e0012922. [Google Scholar] [CrossRef] [PubMed]
  22. Feng, M.; Guo, S.; Fan, S.; Zeng, X.; Zhang, Y.; Liao, Y.; Wang, J.; Zhao, T.; Wang, L.; Che, Y.; et al. The Preferential Infection of Astrocytes by Enterovirus 71 Plays a Key Role in the Viral Neurogenic Pathogenesis. Front. Cell Infect. Microbiol. 2016, 6, 192. [Google Scholar] [CrossRef] [PubMed]
  23. Lai, R.H.; Chow, Y.H.; Chung, N.H.; Chen, T.C.; Shie, F.S.; Juang, J.L. Neurotropic EV71 causes encephalitis by engaging intracellular TLR9 to elicit neurotoxic IL12-p40-iNOS signaling. Cell Death Dis. 2022, 13, 328. [Google Scholar] [CrossRef] [PubMed]
  24. Sun, Q.; Li, J.; Zhang, B.; Wang, R.; Wang, C.; Li, X.; Liu, Y.; Zhang, Y. Transcriptome sequencing analysis of echovirus 30 infection reveals its potential pathogenesis. Front. Microbiol. 2022, 13, 958385. [Google Scholar] [CrossRef] [PubMed]
  25. Hu, Y.; Yang, Z.; Wang, S.; Sun, D.; Zhong, M.; Wen, M.; Song, J.; Zhang, Y. Comprehensive analysis of differential expression profiles via transcriptome sequencing in SH-SY5Y cells infected with CV-A16. PLoS ONE 2020, 15, e0241174. [Google Scholar] [CrossRef]
  26. Hu, Y.; Xu, Y.; Huang, Z.; Deng, Z.; Fan, J.; Yang, R.; Ma, H.; Song, J.; Zhang, Y. Transcriptome sequencing analysis of SH-SY5Y cells infected with EV71 reveals the potential neuropathic mechanisms. Virus Res. 2020, 282, 197945. [Google Scholar] [CrossRef]
  27. Colombo, E.; Farina, C. Astrocytes: Key Regulators of Neuroinflammation. Trends Immunol. 2016, 37, 608–620. [Google Scholar] [CrossRef] [PubMed]
  28. Molofsky, A.V.; Kelley, K.W.; Tsai, H.H.; Redmond, S.A.; Chang, S.M.; Madireddy, L.; Chan, J.R.; Baranzini, S.E.; Ullian, E.M.; Rowitch, D.H. Astrocyte-encoded positional cues maintain sensorimotor circuit integrity. Nature 2014, 509, 189–194. [Google Scholar] [CrossRef]
  29. Luo, Z.; Su, R.; Wang, W.; Liang, Y.; Zeng, X.; Shereen, M.A.; Bashir, N.; Zhang, Q.; Zhao, L.; Wu, K.; et al. EV71 infection induces neurodegeneration via activating TLR7 signaling and IL-6 production. PLoS Pathog. 2019, 15, e1008142. [Google Scholar] [CrossRef]
  30. Zhou, Z.; Ren, L.; Zhang, L.; Zhong, J.; Xiao, Y.; Jia, Z.; Guo, L.; Yang, J.; Wang, C.; Jiang, S.; et al. Heightened Innate Immune Responses in the Respiratory Tract of COVID-19 Patients. Cell Host Microbe 2020, 27, 883–890.e2. [Google Scholar] [CrossRef]
  31. Kuriakose, T.; Zheng, M.; Neale, G.; Kanneganti, T.D. IRF1 Is a Transcriptional Regulator of ZBP1 Promoting NLRP3 Inflammasome Activation and Cell Death during Influenza Virus Infection. J. Immunol. 2018, 200, 1489–1495. [Google Scholar] [CrossRef] [PubMed]
  32. Daniels, B.P.; Kofman, S.B.; Smith, J.R.; Norris, G.T.; Snyder, A.G.; Kolb, J.P.; Gao, X.; Locasale, J.W.; Martinez, J.; Gale, M., Jr.; et al. The Nucleotide Sensor ZBP1 and Kinase RIPK3 Induce the Enzyme IRG1 to Promote an Antiviral Metabolic State in Neurons. Immunity 2019, 50, 64–76.e4. [Google Scholar] [CrossRef] [PubMed]
  33. Hou, J.; Han, L.; Zhao, Z.; Liu, H.; Zhang, L.; Ma, C.; Yi, F.; Liu, B.; Zheng, Y.; Gao, C. USP18 positively regulates innate antiviral immunity by promoting K63-linked polyubiquitination of MAVS. Nat. Commun. 2012, 12, 2970. [Google Scholar] [CrossRef] [PubMed]
  34. Singh, H.; Koury, J.; Kaul, M. Innate Immune Sensing of Viruses and Its Consequences for the Central Nervous System. Viruses 2021, 13, 170. [Google Scholar] [CrossRef] [PubMed]
  35. VanBlargan, L.A.; Himansu, S.; Foreman, B.M.; Ebel, G.D.; Pierson, T.C.; Diamond, M.S. An mRNA Vaccine Protects Mice against Multiple Tick-Transmitted Flavivirus Infections. Cell Rep. 2018, 25, 3382–3392.e3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Body weight, survival rate, and clinical score of suckling mice in inoculation route and F-virus challenge dose experiments. Five-day-old ICR mice were infected with strain F219 at 103, 104, 105, and 106 TCID50 by intracranial injection (A,D,G), intraperitoneal injection (B,E,H), and intramuscular injection (C,F,I), respectively (n = 7 per group); control mice were inoculated with cell culture maintenance medium (MM). Body weight (AC), survival rate (DF) and clinical scores (GI) of 5-day-old mice inoculated via IC, IP, and IM were monitored daily. Clinical scores of CVA6-infected and control mice. 0, no disease; 1, ruffled fur; 2, weight loss; 3, single limb paralysis; 4, hindlimb paralysis; 5, moribund or dead.
Figure 1. Body weight, survival rate, and clinical score of suckling mice in inoculation route and F-virus challenge dose experiments. Five-day-old ICR mice were infected with strain F219 at 103, 104, 105, and 106 TCID50 by intracranial injection (A,D,G), intraperitoneal injection (B,E,H), and intramuscular injection (C,F,I), respectively (n = 7 per group); control mice were inoculated with cell culture maintenance medium (MM). Body weight (AC), survival rate (DF) and clinical scores (GI) of 5-day-old mice inoculated via IC, IP, and IM were monitored daily. Clinical scores of CVA6-infected and control mice. 0, no disease; 1, ruffled fur; 2, weight loss; 3, single limb paralysis; 4, hindlimb paralysis; 5, moribund or dead.
Viruses 15 00511 g001
Figure 2. Virus titers in tissues and histopathological changes of 5-day-old mice. Intramuscular injection of the CVA6 strain F219 was performed in 5-day-old ICR mice. Mice were then sacrificed after 1, 3, and 5 days, and virus titers were measured in the brain and muscle tissue (A,B). Histopathological analysis showing muscle fiber rupture (D,E), cerebral edema, and neuronal cell swelling (G,H); no changes were observed in control tissues (C,F). The magnification of (C,D,F,G) is × 200, and the magnification of (E,H) is ×400. The scale bar for (C,D,F,G) is 100 μm, and the scale bar for (E,H) is 50 μm.
Figure 2. Virus titers in tissues and histopathological changes of 5-day-old mice. Intramuscular injection of the CVA6 strain F219 was performed in 5-day-old ICR mice. Mice were then sacrificed after 1, 3, and 5 days, and virus titers were measured in the brain and muscle tissue (A,B). Histopathological analysis showing muscle fiber rupture (D,E), cerebral edema, and neuronal cell swelling (G,H); no changes were observed in control tissues (C,F). The magnification of (C,D,F,G) is × 200, and the magnification of (E,H) is ×400. The scale bar for (C,D,F,G) is 100 μm, and the scale bar for (E,H) is 50 μm.
Viruses 15 00511 g002
Figure 3. Co-localization of viral antigens with astrocyte marker GFAP in infected brain tissues. Suckling mice were injected with CVA6 or cell culture maintenance medium. Brain tissue sections of suckling mice were fixed and immunostained with DAPI (Blue), VP1 (Red), and astrocyte marker GFAP (Green) (AH); immunostaining with DAPI (Blue), VP1 (Red), and neuron marker NeuN (Green) (IP); Magnifications: (A,B) ×400. Scale bar = 20 μm.
Figure 3. Co-localization of viral antigens with astrocyte marker GFAP in infected brain tissues. Suckling mice were injected with CVA6 or cell culture maintenance medium. Brain tissue sections of suckling mice were fixed and immunostained with DAPI (Blue), VP1 (Red), and astrocyte marker GFAP (Green) (AH); immunostaining with DAPI (Blue), VP1 (Red), and neuron marker NeuN (Green) (IP); Magnifications: (A,B) ×400. Scale bar = 20 μm.
Viruses 15 00511 g003
Figure 4. CVA6 infects U251 cells. Cytopathic effect of virus on U251 at 24 h post infection (A); cells were inoculated with CVA6 at 37 °C, and harvested at 2, 6, 12, 24, and 48 h post-inoculation. Virus titers were determined on RD cells (B); U251 cells infected with CVA6 and control were fixed and immunostained with DAPI (Blue), astrocyte marker (GFAP, Yellow), and VP2 (Green). The representative images were acquired using fluorescence confocal microscopy. Magnifications: (A,B), ×200 × 600. Scale bar = 100 μm (C).
Figure 4. CVA6 infects U251 cells. Cytopathic effect of virus on U251 at 24 h post infection (A); cells were inoculated with CVA6 at 37 °C, and harvested at 2, 6, 12, 24, and 48 h post-inoculation. Virus titers were determined on RD cells (B); U251 cells infected with CVA6 and control were fixed and immunostained with DAPI (Blue), astrocyte marker (GFAP, Yellow), and VP2 (Green). The representative images were acquired using fluorescence confocal microscopy. Magnifications: (A,B), ×200 × 600. Scale bar = 100 μm (C).
Viruses 15 00511 g004
Figure 5. Transcriptome analysis of the brain tissues and U251 cells infected with CVA6. Total RNA was extracted from U251 cells and brain tissues of suckling mice and analyzed by sequencing (A). Hiplot Enhanced MA v0.1.0 analysis showing the comparison of the RNA-Seq of differential gene expression of infected U251 cells (B) and that of infected suckling mice brain tissue (C). GO and KEGG pathway enrichment analysis of the RNA-Seq of differentially expressed genes in the infected U251 cells (D,E) and that of infected brain tissue (F,G). The dot plot displays the enrichment score values. Visualization of differentially expressed genes in the brain tissue (H).
Figure 5. Transcriptome analysis of the brain tissues and U251 cells infected with CVA6. Total RNA was extracted from U251 cells and brain tissues of suckling mice and analyzed by sequencing (A). Hiplot Enhanced MA v0.1.0 analysis showing the comparison of the RNA-Seq of differential gene expression of infected U251 cells (B) and that of infected suckling mice brain tissue (C). GO and KEGG pathway enrichment analysis of the RNA-Seq of differentially expressed genes in the infected U251 cells (D,E) and that of infected brain tissue (F,G). The dot plot displays the enrichment score values. Visualization of differentially expressed genes in the brain tissue (H).
Viruses 15 00511 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sun, Q.; Li, J.; Wang, R.; Sun, T.; Zong, Y.; Wang, C.; Liu, Y.; Li, X.; Song, Y.; Zhang, Y. Coxsackievirus A6 Infection Causes Neurogenic Pathogenesis in a Neonatal Murine Model. Viruses 2023, 15, 511. https://doi.org/10.3390/v15020511

AMA Style

Sun Q, Li J, Wang R, Sun T, Zong Y, Wang C, Liu Y, Li X, Song Y, Zhang Y. Coxsackievirus A6 Infection Causes Neurogenic Pathogenesis in a Neonatal Murine Model. Viruses. 2023; 15(2):511. https://doi.org/10.3390/v15020511

Chicago/Turabian Style

Sun, Qiang, Jichen Li, Rui Wang, Tiantian Sun, Yanjun Zong, Congcong Wang, Ying Liu, Xiaoliang Li, Yang Song, and Yong Zhang. 2023. "Coxsackievirus A6 Infection Causes Neurogenic Pathogenesis in a Neonatal Murine Model" Viruses 15, no. 2: 511. https://doi.org/10.3390/v15020511

APA Style

Sun, Q., Li, J., Wang, R., Sun, T., Zong, Y., Wang, C., Liu, Y., Li, X., Song, Y., & Zhang, Y. (2023). Coxsackievirus A6 Infection Causes Neurogenic Pathogenesis in a Neonatal Murine Model. Viruses, 15(2), 511. https://doi.org/10.3390/v15020511

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop