Diverse RNA viruses in the venom-related microenvironment of different animal phyla.

Venom is known as the source of natural antimicrobial products. Previous studies have largely focused on the expression of venom-related genes and the biochemical components of venom. With the advent of metagenomic sequencing, many more microorganisms, especially viruses, have been identified in highly diverse environments. Herein, we investigated the RNA virome in the venom-related microenvironment through analysis of a large volume of venom-related RNA-sequencing data mined from public databases. From this, we identified viral sequences belonging to thirty-six different viruses, of which twenty-two were classified as 'novel' as they exhibited less than 90 per cent amino acid identity to known viruses in the RNA-dependent RNA polymerase. Most of these novel viruses possessed genome structures similar to their closest relatives, with specific alterations in some cases. Phylogenetic analyses revealed that these viruses belonged to at least twenty-two viral families or unclassified groups, some of which were highly divergent from known taxa. Although further analysis failed to find venom-specific viruses, some viruses seemingly had much higher abundance in the venom-related microenvironment than in other tissues. In sum, our study provides insights into the RNA virome of the venom-related microenvironment from diverse animal phyla.

Diversity and connectedness of brine shrimp viruses in global hypersaline ecosystems.

Brine shrimp (Artemia) has existed on Earth for 400 million years and has major ecological importance in hypersaline ecosystems. As a crucial live food in aquaculture, brine shrimp cysts have become one of the most important aquatic products traded worldwide. However, our understanding of the biodiversity, prevalence and global connectedness of viruses in brine shrimp is still very limited. A total of 143 batches of brine shrimp (belonging to seven species) cysts were collected from six continents including 21 countries and more than 100 geographic locations worldwide during 1977-2019. In total, 55 novel RNA viruses were identified, which could be assigned to 18 different viral families and related clades. Eleven viruses were dsRNA viruses, 16 were +ssRNA viruses, and 28 were-ssRNA viruses. Phylogenetic analyses of the RNA-directed RNA polymerase (RdRp) showed that brine shrimp viruses were often grouped with viruses isolated from other invertebrates and fungi. Remarkably, most brine shrimp viruses were related to those from different hosts that might feed on brine shrimp or share the same ecological niche. A notable case was the novel brine shrimp noda-like virus 3, which shared 79.25% (RdRp) and 63.88% (capsid proteins) amino acid identity with covert mortality nodavirus (CMNV) that may cause losses in aquaculture. In addition, both virome composition and phylogenetic analyses revealed global connectedness in certain brine shrimp viruses, particularly among Asia and Northern America. This highlights the incredible species diversity of viruses in these ancient species and provides essential data for the prevalence of RNA viruses in the global aquaculture industry. More broadly, these findings provide novel insights into the previously unrecognized RNA virosphere in hypersaline ecosystems worldwide and demonstrate that human activity might have driven the global connectedness of brine shrimp viruses.

Metatranscriptomic Sequencing Reveals Host Species as an Important Factor Shaping the Mosquito Virome.

Mosquitoes are important vector hosts for numerous viral pathogens and harbor a large number of mosquito-specific viruses as well as human-infecting viruses. Previous studies have mainly focused on the discovery of mosquito viruses, and our understanding of major ecological factors associated with virome structure in mosquitoes remains limited. We utilized metatranscriptomic sequencing to characterize the viromes of five mosquito species sampled across eight locations in Yunnan Province, China. This revealed the presence of 52 viral species, of which 19 were novel, belonging to 15 viral families/clades. Of particular note was Culex hepacivirus 1, clustering within the avian clade of hepaciviruses. Notably, both the viromic diversity and abundance of Aedes genus mosquitoes were significantly higher than those of the Culex genus, while Aedes albopictus mosquitoes harbored a higher diversity than Aedes aegypti mosquitoes. Our findings thus point to discernible differences in viromic structure between mosquito genera and even between mosquito species within the same genus. Importantly, such differences were not attributable to differences in sampling between geographical location. Our study also revealed the ubiquitous presence of the endosymbiont bacterium Wolbachia, with the genetic diversity and abundance also varying between mosquito species. In conclusion, our results suggested that the mosquito host species play an important role in shaping the virome's structure. IMPORTANCE This study revealed the huge capability of mosquitoes in harboring a rich diversity of RNA viruses, although relevant studies have characterized the intensively unparalleled diversity of RNA viruses previously. Furthermore, our findings showed discernible differences not only in viromic structure between mosquito genera and even between mosquito species within the same genus but also in the genetic diversity and abundance of Wolbachia between different mosquito populations. These findings emphasize the importance of host genetic background in shaping the virome composition of mosquitoes.

Nationwide genomic surveillance reveals the prevalence and evolution of honeybee viruses in China.

BACKGROUND: The economic and environmental value of honeybees has been severely challenged in recent years by the collapse of their colonies worldwide, often caused by outbreaks of infectious diseases. However, our understanding of the diversity, prevalence, and transmission of honeybee viruses is largely obscure due to a lack of large-scale and longitudinal genomic surveillance on a global scale. RESULTS: We report the meta-transcriptomic sequencing of nearly 2000 samples of the two most important economic and widely maintained honeybee species, as well as an associated ectoparasite mite, collected across China during 2016-2019. We document the natural diversity and evolution of honeybee viruses in China, providing evidence that multiple viruses commonly co-circulate within individual bee colonies. We also expanded the genomic data for 12 important honeybee viruses and revealed novel genetic variants and lineages associated with China. We identified more than 23 novel viruses from the honeybee and mite viromes, with some exhibiting ongoing replication in their respective hosts. Together, these data provide additional support to the idea that mites are an important reservoir and spill-over host for honeybee viruses. CONCLUSIONS: Our data show that honeybee viruses are more widespread, prevalent, and genetically diverse than previously realized. The information provided is important in mitigating viral infectious diseases in honeybees, in turn helping to maintain sustainable productive agriculture on a global scale. Video Abstract.

Diagnosis and analysis of unexplained cases of childhood encephalitis in Australia using metatranscriptomic sequencing.

Encephalitis is most often caused by a variety of infectious agents identified through diagnostic tests utilizing cerebrospinal fluid. We investigated the clinical characteristics and potential aetiological agents of unexplained encephalitis through metagenomic sequencing of residual clinical samples from multiple tissue types and independent clinical review. Forty-three specimens were collected from 18 encephalitis cases with no cause identified by the Australian Childhood Encephalitis study. Samples were subjected to total RNA sequencing ('metatranscriptomics') to determine the presence and abundance of potential pathogens, and to describe the possible aetiologies of unexplained encephalitis. Using this protocol, we identified five RNA and two DNA viruses associated with human infection from both non-sterile and sterile sites, which were confirmed by PCR. These comprised two human rhinoviruses, two human seasonal coronaviruses, two polyomaviruses and one picobirnavirus. Human rhinovirus and seasonal coronaviruses may be responsible for five of the encephalitis cases. Immune-mediated encephalitis was considered likely in six cases and metatranscriptomics did not identify a possible pathogen in these cases. The aetiology remained unknown in nine cases. Our study emphasizes the importance of respiratory viruses in the aetiology of unexplained child encephalitis and suggests that non-central-nervous-system sampling in encephalitis clinical guidelines and protocols could improve the diagnostic yield.

A Novel Bunyavirus Discovered in Oriental Shrimp (Penaeus chinensis).

Herein, we describe a novel bunyavirus, oriental wenrivirus 1 (OWV1), discovered in moribund oriental shrimp (Penaeus chinensis) collected from a farm in China in 2016. Like most bunyaviruses, OWV1 particles were enveloped, spherical- to ovoid-shaped, and 80-115 nm in diameter. However, its genome was found to comprise four segments of (-)ssRNA. These included an L RNA segment (6,317 nt) encoding an RNA-directed RNA polymerase (RdRp) of 2,052 aa, an M RNA segment (2,978 nt) encoding a glycoprotein precursor (GPC) of 922 aa, an S1 RNA segment (1,164 nt) encoding a nucleocapsid (N) protein of 243 aa, and an S2 RNA segment (1,382 nt) encoding a putative non-structural (NSs2) protein of 401 aa. All the four OWV1 RNA segments have complementary terminal decanucleotides (5'-ACACAAAGAC and 3'-UGUGUUUCUG) identical to the genomic RNA segments of uukuviruses and similar to those of phleboviruses and tenuiviruses in the Phenuiviridae. Phylogenetic analyses revealed that the RdRp, GPC, and N proteins of OWV1 were closely related to Wenzhou shrimp virus 1 (WzSV-1) and Mourilyan virus (MoV) that infect black tiger shrimp (P. monodon). Phylogenetic analyses also suggested that OWV1 could be classified into a second, yet to be established, species of the Wenrivirus genus in the Phenuiviridae. These wenriviruses also clustered with Wenling crustacean virus 7 from shrimps and bunya-like brown spot virus from white-clawed crayfish. Of note there were no homologs of the NSs2 of OWV1 and MoV/WzSV-1 in GenBank, and whether other crustacean phenuiviruses also possess a similar S2 RNA segment warrants further investigation. In addition, we established a TaqMan probe-based reverse-transcription quantitative PCR method for detection of OWV1, and it was detected as 1.17 × 102-1.90 × 107 copies/ng-RNA in gills of 23 out of 32 P. chinensis samples without an obvious gross sign. However, the discovery of OWV1 highlights the expanding genomic diversity of bunyaviruses.

Metagenomic Identification of Viral Sequences in Laboratory Reagents.

Metagenomic next-generation sequencing has transformed the discovery and diagnosis of infectious disease, with the power to characterise the complete 'infectome' (bacteria, viruses, fungi, parasites) of an individual host organism. However, the identification of novel pathogens has been complicated by widespread microbial contamination in commonly used laboratory reagents. Using total RNA sequencing ("metatranscriptomics") we documented the presence of contaminant viral sequences in multiple 'blank' negative control sequencing libraries that comprise a sterile water and reagent mix. Accordingly, we identified 14 viral sequences in 7 negative control sequencing libraries. As in previous studies, several circular replication-associated protein encoding (CRESS) DNA virus-like sequences were recovered in the blank control libraries, as well as contaminating sequences from the Totiviridae, Tombusviridae and Lentiviridae families of RNA virus. These data suggest that viral contamination of common laboratory reagents is likely commonplace and can comprise a wide variety of viruses.

High-resolution Metatranscriptomic Characterization of the Pulmonary RNA Virome After Lung Transplantation.

BACKGROUND: Lung transplantation provides a unique opportunity to investigate the constituents and temporal dynamics of the human pulmonary microbiome after lung transplantation. For methodological reasons, prior studies using metagenomics have detected DNA viruses but not demonstrated the presence of RNA viruses, including those that are common community acquired. In this proof-of-concept study, we aimed to further characterize the pulmonary microbiome after lung transplantation by using metagenomic next-generation sequencing (mNGS), with a particular focus on the RNA virome. METHODS: We performed a single-center longitudinal study of lower respiratory tract RNA viruses and bacteria using bronchoalveolar lavage at postoperative day 1 and week 6 analyzed with total RNA sequencing (metatranscriptomics). Five primary and 5 repeat transplant recipients were recruited. RESULTS: mNGS identified 5 RNA viruses (nil in the normal saline control), including 4 species of human rhinovirus not previously reported in Australia: A7 (HRV-A7), C22 (HRV-C22), B52 (HRV-B52), and B72 (HRV-B72). Overall, 12/20 specimens were virus positive in 7/10 cases. Human parainfluenza virus 3 was the most frequent virus in 7/20 specimens in 5/10 cases. In this small study, we did not detect a significant difference in abundance and diversity of RNA viruses and bacteria at postoperative day 1 and 6 wk, nor differences between retransplant recipients and primary lung transplant recipients. CONCLUSIONS: Our study demonstrates how mNGS can also identify RNA viruses within the human pulmonary virome, including novel RNA viruses, and paves the way for a greater understanding of the complex relationships among the constituents of the pulmonary infectome.

A Novel Virus of Flaviviridae Associated with Sexual Precocity in Macrobrachium rosenbergii.

Since 2010, sexual precocity, a typical sign of the iron prawn syndrome (IPS), resulting in the reduced size of farmed giant freshwater prawns Macrobrachium rosenbergii, has caused substantial production losses. However, the cause of IPS was not clear. We ran tests for eight major shrimp pathogens, but none were detected from IPS-affected prawns. We performed the histopathological examination of tissues and identified an eosinophilic inclusion in the perinuclear cytoplasm of cells in various tissues associated with nervous and endocrinal functions in the compound eyes. A subsequent bioassay with viral extracts of IPS-affected samples reproduced the gross signs of IPS. Metatranscriptomic sequencing identified a novel virus of Flaviviridae in all IPS-affected M. rosenbergii prawns, which was not found in samples without IPS. This virus contains a positive-sense, single-stranded RNA genome of 12,630 nucleotides (nt). Phylogenetic analysis of the conserved RdRp and NS3 domains showed that it may belong to a new genus between Jingmenvirus and Flavivirus. Under transmission electron microscopy (TEM), putative virus particles showed as spherical with a diameter of 40 to 60 nm. In situ hybridization found hybridization signals consistent with the histopathology in the compound eyes from IPS-affected M. rosenbergii. We provisionally name this virus infectious precocity virus (IPV) and propose the binominal Latin name Crustaflavivirus infeprecoquis gen. nov., sp. nov. We developed a nested reverse transcription-PCR diagnostic assay and confirmed that all IPS-affected prawns tested IPV positive but normal prawns tested negative. Collectively, our study revealed a novel virus of Flaviviridae associated with sexual precocity in M. rosenbergii. IMPORTANCE The iron prawn syndrome (IPS), also described as sexual precocity, results in the reduced size of farmed prawns at harvest and significant economic losses. IPS has been frequently reported in populations of farmed Macrobrachium rosenbergii since 2010, but the cause was heretofore unknown. Here, we reported a novel virus identified from prawns with IPS using infection experiments, metatranscriptomic sequencing, and transmission electron microscopy and provisionally named it infectious precocity virus (IPV). Phylogenetic analysis showed that IPV represents a new genus, proposed as Crustaflavivirus gen. nov., in the family Flaviviridae. This study provides novel insight that a viral infection may cause pathological change and sexual maturation and subsequently affect crustacean growth. Therefore, we call for quarantine inspection of IPV in transboundary trade of live M. rosenbergii and enhanced surveillance of IPV in aquaculture in the region and globally.

Bioinformatics resources for SARS-CoV-2 discovery and surveillance.

In early January 2020, the novel coronavirus (SARS-CoV-2) responsible for a pneumonia outbreak in Wuhan, China, was identified using next-generation sequencing (NGS) and readily available bioinformatics pipelines. In addition to virus discovery, these NGS technologies and bioinformatics resources are currently being employed for ongoing genomic surveillance of SARS-CoV-2 worldwide, tracking its spread, evolution and patterns of variation on a global scale. In this review, we summarize the bioinformatics resources used for the discovery and surveillance of SARS-CoV-2. We also discuss the advantages and disadvantages of these bioinformatics resources and highlight areas where additional technical developments are urgently needed. Solutions to these problems will be beneficial not only to the prevention and control of the current COVID-19 pandemic but also to infectious disease outbreaks of the future.

Meta-Transcriptomic Discovery of a Divergent Circovirus and a Chaphamaparvovirus in Captive Reptiles with Proliferative Respiratory Syndrome.

Viral pathogens are being increasingly described in association with mass morbidity and mortality events in reptiles. However, our knowledge of reptile viruses remains limited. Herein, we describe the meta-transcriptomic investigation of a mass morbidity and mortality event in a colony of central bearded dragons (Pogona vitticeps) in 2014. Severe, extensive proliferation of the respiratory epithelium was consistently found in affected dragons. Similar proliferative lung lesions were identified in bearded dragons from the same colony in 2020 in association with increased intermittent mortality. Total RNA sequencing identified two divergent DNA viruses: a reptile-infecting circovirus, denoted bearded dragon circovirus (BDCV), and the first exogeneous reptilian chaphamaparvovirus-bearded dragon chaphamaparvovirus (BDchPV). Phylogenetic analysis revealed that BDCV was most closely related to bat-associated circoviruses, exhibiting 70% amino acid sequence identity in the Replicase (Rep) protein. In contrast, in the nonstructural (NS) protein, the newly discovered BDchPV showed approximately 31%-35% identity to parvoviruses obtained from tilapia fish and crocodiles in China. Subsequent specific PCR assays revealed BDCV and BDchPV in both diseased and apparently normal captive reptiles, although only BDCV was found in those animals with proliferative pulmonary lesions and respiratory disease. This study expands our understanding of viral diversity in captive reptiles.

Molecular epidemiology of rotavirus A and adenovirus among children with acute diarrhea in Hangzhou, China.

BACKGROUND: Rotavirus A (RVA) and adenovirus (Adv) are important causes of acute diarrhea in children. RVAs are classified into G and P genotypes based on viral proteins (VP)7 and VP4 gene and Adv contains over 70 genotypes based on hexon and fiber gene. This study aimed to characterize the molecular epidemiology of RVA and Adv in children with acute diarrhea during 2017-2018 in Hangzhou. METHODS: The stool samples were collected and tested for RVA and Adv by reverse transcription-quantitative PCR (RT-qPCR) assay. The RVA positive samples were detected by RT-PCR for VP7(G) and VP4([P]) genotypes, and the Adv positive samples were detected by PCR for genotyping by the target to hexon gene. RESULTS: Among 228 RVA-positive samples, G9 was detected as the most frequent genotype (195/228, 85.5%), followed by G3 (20/228, 8.8%), G2 (7/228, 3.1%) and G1 (6/228, 2.6%). G9 strains were closely related to strains from China and neighboring countries, as well as the USA. On the other hand, P[8] strains were detected in 219 (96.1%) samples with most closely related to one strain from Malawi, and P[4] in 9 (3.9%) samples. G9P[8] (84.6%, 193/228) was the most prevalent rotavirus A strains, followed by G3P[8] (8.8%, 20/228), G2P[4] (3.1%, 7/228), G1P[8] (2.6%, 6/228) and G9P[4] (0.9%, 2/228). Of 167 Adv-positive cases, 2 different genotypes were identified with 152 (91.0%) of Adv-41and 15 (9%) of Adv-40. All Adv strains were closely related to prototype strains of Adv types 40 and 41 in India. CONCLUSIONS: G9P[8] of RVA and Adv-41 were the most common genotypes that caused children's acute diarrhea in Hangzhou, 2017-2018.

High resolution metagenomic characterization of complex infectomes in paediatric acute respiratory infection.

The diversity of pathogens associated with acute respiratory infection (ARI) makes diagnosis challenging. Traditional pathogen screening tests have a limited detection range and provide little additional information. We used total RNA sequencing ("meta-transcriptomics") to reveal the full spectrum of microbes associated with paediatric ARI. Throat swabs were collected from 48 paediatric ARI patients and 7 healthy controls. Samples were subjected to meta-transcriptomics to determine the presence and abundance of viral, bacterial, and eukaryotic pathogens, and to reveal mixed infections, pathogen genotypes/subtypes, evolutionary origins, epidemiological history, and antimicrobial resistance. We identified 11 RNA viruses, 4 DNA viruses, 4 species of bacteria, and 1 fungus. While most are known to cause ARIs, others, such as echovirus 6, are rarely associated with respiratory disease. Co-infection of viruses and bacteria and of multiple viruses were commonplace (9/48), with one patient harboring 5 different pathogens, and genome sequence data revealed large intra-species diversity. Expressed resistance against eight classes of antibiotic was detected, with those for MLS, Bla, Tet, Phe at relatively high abundance. In summary, we used a simple total RNA sequencing approach to reveal the complex polymicrobial infectome in ARI. This provided comprehensive and clinically informative information relevant to understanding respiratory disease.

Identification of A Novel Papillomavirus Associated with Squamous Cell Carcinoma in A Domestic Cat.

Papillomaviruses infect the skin and mucosal surfaces of diverse animal hosts with consequences ranging from asymptomatic colonization to highly malignant epithelial cancers. Increasing evidence suggests a role for papillomaviruses in the most common cutaneous malignancy of domestic cats, squamous cell carcinoma (SCC). Using total DNA sequencing we identified a novel feline papillomavirus in a nasal biopsy taken from a cat presenting with both nasal cavity lymphoma and recurrent squamous cell carcinoma affecting the nasal planum. We designate this novel virus as Felis catus papillomavirus 6 (FcaPV6). The complete FcaPV6 7453 bp genome was similar to those of other feline papillomaviruses and phylogenetic analysis revealed that it was most closely related to FcaPV3, although was distinct enough to represent a new viral type. Classification of FcaPV6 in a new genus alongside FcaPVs 3, 4 and 5 is supported. Archived excisional biopsy of the SCC, taken 20 months prior to presentation, was intensely positive on p16 immunostaining. FcaPV6, amplified using virus-specific, but not consensus, PCR, was the only papillomavirus detected in DNA extracted from the SCC. Conversely, renal lymphoma, sampled at necropsy two months after presentation, tested negative on FcaPV6-specific PCR. In sum, using metagenomics we demonstrate the presence of a novel feline papillomavirus in association with cutaneous squamous cell carcinoma.

Identification of a Novel Equine Papillomavirus in Semen from a Thoroughbred Stallion with a Penile Lesion.

Papillomaviruses (PVs) have been identified in a wide range of animal species and are associated with a variety of disease syndromes including classical papillomatosis, aural plaques, and genital papillomas. In horses, 13 PVs have been described to date, falling into six genera. Using total RNA sequencing (meta-transcriptomics) we identified a novel equine papillomavirus in semen taken from a thoroughbred stallion suffering a genital lesion, which was confirmed by nested RT-PCR. We designate this novel virus Equus caballus papillomavirus 9 (EcPV9). The complete 7656 bp genome of EcPV9 exhibited similar characteristics to those of other horse papillomaviruses. Phylogenetic analysis based on concatenated E1-E2-L2-L1 amino acid sequences revealed that EcPV9 clustered with EcPV2, EcPV4, and EcPV5, although was distinct enough to represent a new viral species within the genus Dyoiotapapillomavirus (69.35%, 59.25%, and 58.00% nucleotide similarity to EcPV2, EcPV4, and EcPV5, respectively). In sum, we demonstrate the presence of a novel equine papillomavirus for which more detailed studies of disease association are merited.

Molecular epidemiology of enterovirus from children with herpangina or hand, foot, and mouth disease in Hangzhou, 2016.

Enteroviruses (EVs) are the major cause of hand, foot, and mouth disease (HFMD) and herpangina in children. In this study, we conducted a molecular investigation of EVs in throat swab samples from children in Hangzhou, China with a diagnosis of HFMD or herpangina. EVs were detected using one-step real-time RT-PCR, and their serotypes were determined based on partial VP1 gene sequences. The molecular typing results revealed the presence of six different EV serotypes in HFMD cases, including coxsackievirus (CV) A16 (20/30, 66.7%), CVA4 (3/30, 10.0%), CVA6 (3/30, 10.0%), EVA71 (2/30, 6.7%), CVB4 (1/30, 3.3%), and CVB5 (1/30, 3.3%). Eleven different EV serotypes were detected in herpangina cases, among which CVA4 was the most frequently detected serotype (105/170, 61.8%), followed by CVA16 (30/170, 17.6%), CVB4 (9/170, 5.3%), CVA6 (6/170, 3.5%), CVB3 (5/170, 2.9%), CVA10 (3/170, 1.8%), EVA71 (4/170, 2.4%), Echo9 (3/170, 1.8%), CVA9 (2/170, 1.2%), CVB1 (3/170, 1.8%) and CVA5 (1/170, 0.6%). The nucleotide sequence identity of EV strains from the same subtype ranged from 80.7% to 100%, and most of the EVs were closely related to virus strains found in Australia and mainland China. In conclusion, CVA 16 and CVA 4 were the main serotypes causing HFMD and herpangina, respectively, in children in Hangzhou in 2016. Most of these EVs were closely related to virus strains from Australia and mainland China.

Establishment and characterization of an oral gerbil model for a non-mouse-adapted enterovirus 71 strain.

Enterovirus 71 (EV71) is one of the major pathogens causing hand, foot, and mouth disease (HFMD) with neurological and systemic complications worldwide, and it is mostly discovered in infants and young children. It is of great significance to establish suitable animal models of EV71 infection on research of distribution and pathogenesis of the virus. In this study, we established a successful infection of a non-mouse-adapted isolate of EV71 via oral route in 7-day-old Mongolian gerbil (Meriones unguiculatus), all of which were paralyzed and died within 10 days post infection. Analysis of virus loads in twelve tissues showed that virus was first detected in intestine, blood, heart, lung, and muscle one day post-infection, and then in the rest of the tissues/organs within the next few days, among which thymus, spleen, spinal cord and muscles had the highest virus titer at 5 days post infection. Pathological examination showed that severe necrosis was observed in skeletal muscle and spinal cord, and edema was observed in both heart and lung. Comparisons of host gene expression of various tissues from infected and non-infected gerbils revealed a general up-regulation of genes related to anti-viral response and a viral receptor gene (sialic acid-linked glycans), as well as a tissue(gut)-specific up-regulation of genes related to neuronal communication. Collectively, our results showed that EV71 could induce severe neurological complications as well as massive tissue damage all over the body, which indicates that oral infection of 7-day gerbils can be a suitable animal model to study the infection of EV71 in human.