中国部分地区蝙蝠携带病毒的宏基因组学分析

蝙蝠携带有60多种病毒,其中许多对人有高度致病性.为了解中国蝙蝠携带病毒的自然本底、蝙蝠病毒的多样性和挖掘潜在的病毒病原,通过基于Solexa高通量测序的病毒宏基因组学技术对从吉林、云南、湖南采集的蝙蝠组织进行病毒组学研究,获得了11 644 232条读长(Reads),并拼接出44 872条重叠序列(Contig).通过核酸序列注释发现,其中8.2％(4 002/44 872)的重叠序列与病毒相关,能进一步注释到36个病毒科,包括19种脊椎动物病毒、6种植物病毒、4种昆虫病毒和4种噬菌体.通过对重叠序列的遗传进化分析、多序列比对显示,被注释为细小病毒、腺联病毒、博卡病毒、腺病毒、小双节RNA病毒等的重叠序列与已知病毒相似,部分序列却又呈现出明显的序列差异.通过对腺病毒和博卡病毒进一步的PCR扩增证实了此研究方法可靠.旨在了解我国蝙蝠携带病毒组的构成,对建立高效的野生动物源人兽共患病的监测方法提供参考.     

A Preliminary Study of Viral Metagenomics of French Bat Species in Contact with Humans: Identification of New Mammalian Viruses

The prediction of viral zoonosis epidemics has become a major public health issue. A profound understanding of the viral population in key animal species acting as reservoirs represents an important step towards this goal. Bats harbor diverse viruses, some of which are of particular interest because they cause severe human diseases. However, little is known about the diversity of the global population of viruses found in bats (virome). We determined the viral diversity of five different French insectivorous bat species (nine specimens in total) in close contact with humans. Sequence-independent amplification, high-throughput sequencing with Illumina technology and a dedicated bioinformatics analysis pipeline were used on pooled tissues (brain, liver and lungs). Comparisons of the sequences of contigs and unassembled reads provided a global taxonomic distribution of virus-related sequences for each sample, highlighting differences both within and between bat species. Many viral families were present in these viromes, including viruses known to infect bacteria, plants/fungi, insects or vertebrates, the most relevant being those infecting mammals (Retroviridae, Herpesviridae, Bunyaviridae, Poxviridae, Flaviviridae, Reoviridae, Bornaviridae, Picobirnaviridae). In particular, we detected several new mammalian viruses, including rotaviruses, gammaretroviruses, bornaviruses and bunyaviruses with the identification of the first bat nairovirus. These observations demonstrate that bats naturally harbor viruses from many different families, most of which infect mammals. They may therefore constitute a major reservoir of viral diversity that should be analyzed carefully, to determine the role played by bats in the spread of zoonotic viral infections.     

Metagenomic analysis of viruses from bat fecal samples reveals many novel viruses in insectivorous bats in China

Increasing data indicate that bats harbor diverse viruses, some of which cause severe human diseases. In this study, sequence-independent amplification and high-throughput sequencing (Solexa) were applied to the metagenomic analysis of viruses in bat fecal samples collected from 6 locations in China. A total of 8,746,417 reads with a length of 306,124,595 bp were obtained. Among these reads, 13,541 (0.15%) had similarity to phage sequences and 9,170 (0.1%) had similarity to eukaryotic virus sequences. A total of 129 assembled contigs (>100 nucleotides) were constructed and compared with GenBank: 32 contigs were related to phages, and 97 were related to eukaryotic viruses. The most frequent reads and contigs related to eukaryotic viruses were homologous to densoviruses, dicistroviruses, coronaviruses, parvoviruses, and tobamoviruses, a range that includes viruses from invertebrates, vertebrates, and plants. Most of the contigs had low identities to known viral genomic or protein sequences, suggesting that a large number of novel and genetically diverse insect viruses as well as putative mammalian viruses are transmitted by bats in China. This study provides the first preliminary understanding of the virome of some bat populations in China, which may guide the discovery and isolation of novel viruses in the future.     

Metagenomic analysis of the viromes of three North American bat species: viral diversity among different bat species that share a common habitat

Effective prediction of future viral zoonoses requires an in-depth understanding of the heterologous viral population in key animal species that will likely serve as reservoir hosts or intermediates during the next viral epidemic. The importance of bats as natural hosts for several important viral zoonoses, including Ebola, Marburg, Nipah, Hendra, and rabies viruses and severe acute respiratory syndrome-coronavirus (SARS-CoV), has been established; however, the large viral population diversity (virome) of bats has been partially determined for only a few of the ～1,200 bat species. To assess the virome of North American bats, we collected fecal, oral, urine, and tissue samples from individual bats captured at an abandoned railroad tunnel in Maryland that is cohabitated by 7 to 10 different bat species. Here, we present preliminary characterization of the virome of three common North American bat species, including big brown bats (Eptesicus fuscus), tricolored bats (Perimyotis subflavus), and little brown myotis (Myotis lucifugus). In samples derived from these bats, we identified viral sequences that were similar to at least three novel group 1 CoVs, large numbers of insect and plant virus sequences, and nearly full-length genomic sequences of two novel bacteriophages. These observations suggest that bats encounter and disseminate a large assortment of viruses capable of infecting many different animals, insects, and plants in nature.     

Deciphering the bat virome catalog to better understand the ecological diversity of bat viruses and the bat origin of emerging infectious diseases

Studies have demonstrated that ~60%–80% of emerging infectious diseases (EIDs) in humans originated from wild life. Bats are natural reservoirs of a large variety of viruses, including many important zoonotic viruses that cause severe diseases in humans and domestic animals. However, the understanding of the viral population and the ecological diversity residing in bat populations is unclear, which complicates the determination of the origins of certain EIDs. Here, using bats as a typical wildlife reservoir model, virome analysis was conducted based on pharyngeal and anal swab samples of 4440 bat individuals of 40 major bat species throughout China. The purpose of this study was to survey the ecological and biological diversities of viruses residing in these bat species, to investigate the presence of potential bat-borne zoonotic viruses and to evaluate the impacts of these viruses on public health. The data obtained in this study revealed an overview of the viral community present in these bat samples. Many novel bat viruses were reported for the first time and some bat viruses closely related to known human or animal pathogens were identified. This genetic evidence provides new clues in the search for the origin or evolution pattern of certain viruses, such as coronaviruses and noroviruses. These data offer meaningful ecological information for predicting and tracing wildlife-originated EIDs.     

Duration of Maternal Antibodies against Canine Distemper Virus and Hendra Virus in Pteropid Bats

Old World frugivorous bats have been identified as natural hosts for emerging zoonotic viruses of significant public health concern, including henipaviruses (Nipah and Hendra virus), Ebola virus, and Marburg virus. Epidemiological studies of these viruses in bats often utilize serology to describe viral dynamics, with particular attention paid to juveniles, whose birth increases the overall susceptibility of the population to a viral outbreak once maternal immunity wanes. However, little is understood about bat immunology, including the duration of maternal antibodies in neonates. Understanding duration of maternally derived immunity is critical for characterizing viral dynamics in bat populations, which may help assess the risk of spillover to humans. We conducted two separate studies of pregnant Pteropus bat species and their offspring to measure the half-life and duration of antibodies to 1) canine distemper virus antigen in vaccinated captive Pteropus hypomelanus; and 2) Hendra virus in wild-caught, naturally infected Pteropus alecto. Both of these pteropid bat species are known reservoirs for henipaviruses. We found that in both species, antibodies were transferred from dam to pup. In P. hypomelanus pups, titers against CDV waned over a mean period of 228.6 days (95% CI: 185.4–271.8) and had a mean terminal phase half-life of 96.0 days (CI 95%: 30.7–299.7). In P. alecto pups, antibodies waned over 255.13 days (95% CI: 221.0–289.3) and had a mean terminal phase half-life of 52.24 days (CI 95%: 33.76–80.83). Each species showed a duration of transferred maternal immunity of between 7.5 and 8.5 months, which was longer than has been previously estimated. These data will allow for more accurate interpretation of age-related Henipavirus serological data collected from wild pteropid bats.     

Metagenomic sequencing reveals viral abundance and diversity in mosquitoes from the Shaanxi-Gansu-Ningxia region,China

BackgroundMosquitoes host and transmit numerous arthropod-borne viruses (arboviruses) that cause disease in both humans and animals. Effective surveillance of virome profiles in mosquitoes is vital to the prevention and control of mosquito-borne diseases in northwestern China, where epidemics occur frequently.MethodsMosquitoes were collected in the Shaanxi-Gansu-Ningxia region (Shaanxi Province, Gansu Province, and Ningxia Hui Autonomous Region) of China from June to August 2019. Morphological methods were used for taxonomic identification of mosquito species. High-throughput sequencing and metagenomic analysis were used to characterize mosquito viromes.ResultsA total of 22,959 mosquitoes were collected, including Culex pipiens (45.7%), Culex tritaeniorhynchus (40.6%), Anopheles sinensis (8.4%), Aedes (5.2%), and Armigeres subalbatus (0.1%). In total, 3,014,183 (0.95% of clean reads) viral sequences were identified and assigned to 116 viral species (including pathogens such as Japanese encephalitis virus and Getah virus) in 31 viral families, including Flaviviridae, Togaviridae, Phasmaviridae, Phenuiviridae, and some unclassified viruses. Mosquitoes collected in July (86 species in 26 families) showed greater viral diversity than those from June and August. Culex pipiens (69 species in 25 families) and Culex tritaeniorhynchus (73 species in 24 families) carried more viral species than Anopheles sinensis (50 species in 19 families) or Aedes (38 species in 20 families) mosquitoes.ConclusionViral diversity and abundance were affected by mosquito species and collection time. The present study elucidates the virome compositions of various mosquito species in northwestern China, improving the understanding of virus transmission dynamics for comparison with those of disease outbreaks.     

Metagenomic analysis of viral community in the Yangtze River expands known eukaryotic and prokaryotic virus diversity in freshwater

Viruses in aquatic ecosystems are characterized by extraordinary abundance and diversity. Thus far, there have been limited studies focused on viral communities in river water systems. Here, we investigated the virome of the Yangtze River Delta using viral metagenomic analysis. The compositions of viral communities from six sampling sites were analyzed and compared. By using library construction and next generation sequencing, contigs and singlet reads similar to viral sequences were classified into 17 viral families, including nine dsDNA viral families, four ssDNA viral families and four RNA viral families. Statistical analysis using Friedman test suggested that there was no significant difference among the six sampling sites (P > 0.05). The viromes in this study were all dominated by the order Caudovirales, and a group of Freshwater phage uvFW species were particularly prevalent among all the samples. The virome from Nanjing presented a unique pattern of viral community composition with a relatively high abundance of family Parvoviridae. Phylogenetic analyses based on virus hallmark genes showed that the Caudovirales order and CRESS-DNA viruses presented high genetic diversity, while viruses in the Microviridae and Parvoviridae families and the Riboviria realm were relatively conservative. Our study provides the first insight into viral community composition in large river ecosystem, revealing the diversity and stability of river water virome, contributing to the proper utilization of freshwater resource.     

Metagenomic Profile of the Viral Communities in Rhipicephalus spp. Ticks from Yunnan,China

Besides mosquitoes, ticks are regarded as the primary source of vector-borne infectious diseases. Indeed, a wide variety of severe infectious human diseases, including those involving viruses, are transmitted by ticks in many parts of the world. To date, there are no published reports on the use of next-generation sequencing for studying viral diversity in ticks or discovering new viruses in these arthropods from China. Here, Ion-torrent sequencing was used to investigate the presence of viruses in three Rhipicephalus spp. tick pools (NY-11, NY-13, and MM-13) collected from the Menglian district of Yunnan, China. The sequencing run resulted in 3,641,088, 3,106,733, and 3,871,851 reads in each tick pool after trimming. Reads and assembled contiguous sequences (contigs) were subject to basic local alignment search tool analysis against the GenBank database. Large numbers of reads and contigs related to known viral sequences corresponding to a broad range of viral families were identified. Some of the sequences originated from viruses that have not been described previously in ticks. Our findings will facilitate better understanding of the tick virome, and add to our current knowledge of disease-causing viruses in ticks living under natural conditions.     

Viral metagenomics analysis of planktonic viruses in East Lake, Wuhan, China

East Lake (Lake Donghu), located in Wuhan, China, is a typical city freshwater lake that has been experiencing eutrophic conditions and algal blooming during recent years. Marine and fresh water are considered to contain a large number of viruses. However, little is known about their genetic diversity because of the limited techniques for culturing viruses. In this study, we conducted a viral metagenomic analysis using a high-throughput sequencing technique with samples collected from East Lake in Spring, Summer, Autumn, and Winter. The libraries from four samples each generated 234,669, 71,837, 12,820, and 34,236 contigs (> 90 bp each), respectively. The genetic structure of the viral community revealed a high genetic diversity covering 23 viral families, with the majority of contigs homologous to DNA viruses, including members of Myoviridae, Podoviridae, Siphoviridae, Phycodnaviridae, and Microviridae, which infect bacteria or algae, and members of Circoviridae, which infect invertebrates and vertebrates. The highest viral genetic diversity occurred in samples collected in August, then December and June, and the least diversity in March. Most contigs have low-sequence identities with known viruses. PCR detection targeting the conserved sequences of genes (g20, psbA, psbD, and DNApol) of cyanophages further confirmed that there are novel cyanophages in the East Lake. Our viral metagenomic data provide the first preliminary understanding of the virome in one freshwater lake in China and would be helpful for novel virus discovery and the control of algal blooming in the future.     

Bat Rabies in Guatemala

Rabies in bats is considered enzootic throughout the New World, but few comparative data are available for most countries in the region. As part of a larger pathogen detection program, enhanced bat rabies surveillance was conducted in Guatemala, between 2009 and 2011. A total of 672 bats of 31 species were sampled and tested for rabies. The prevalence of rabies virus (RABV) detection among all collected bats was low (0.3%). Viral antigens were detected and infectious virus was isolated from the brains of two common vampire bats (Desmodus rotundus). RABV was also isolated from oral swabs, lungs and kidneys of both bats, whereas viral RNA was detected in all of the tissues examined by hemi-nested RT-PCR except for the liver of one bat. Sequencing of the nucleoprotein gene showed that both viruses were 100% identical, whereas sequencing of the glycoprotein gene revealed one non-synonymous substitution (302_T,S). The two vampire bat RABV isolates in this study were phylogenetically related to viruses associated with vampire bats in the eastern states of Mexico and El Salvador. Additionally, 7% of sera collected from 398 bats demonstrated RABV neutralizing antibody. The proportion of seropositive bats varied significantly across trophic guilds, suggestive of complex intraspecific compartmentalization of RABV perpetuation.     

Diversity of Bartonella and Rickettsia spp. in Bats and Their Blood-Feeding Ectoparasites from South Africa and Swaziland

In addition to several emerging viruses, bats have been reported to host multiple bacteria but their zoonotic threats remain poorly understood, especially in Africa where the diversity of bats is important. Here, we investigated the presence and diversity of Bartonella and Rickettsia spp. in bats and their ectoparasites (Diptera and Siphonaptera) collected across South Africa and Swaziland. We collected 384 blood samples and 14 ectoparasites across 29 different bat species and found positive samples in four insectivorous and two frugivorous bat species, as well as their Nycteribiidae flies. Phylogenetic analyses revealed diverse Bartonella genotypes and one main group of Rickettsia, distinct from those previously reported in bats and their ectoparasites, and for some closely related to human pathogens. Our results suggest a differential pattern of host specificity depending on bat species. Bartonella spp. identified in bat flies and blood were identical supporting that bat flies may serve as vectors. Our results represent the first report of bat-borne Bartonella and Rickettsia spp. in these countries and highlight the potential role of bats as reservoirs of human bacterial pathogens.     

Viral Communities Among Sympatric Vampire Bats and Cattle

Vampire bats are the only mammals known to feed exclusively on blood from other animals, often from domestic cattle. We tested the hypothesis that the adaptation of vampire bats to hematophagy would have resulted in shared viral communities among vampire bats and cattle, as a direct result of historic spillover events occurring due to hematophagy. We analyzed the presence of different viruses in sample populations of sympatric bat and prey populations and searched for shared viruses between taxa. A limited number of DNA viral groups were detected within each species. However, there was no evidence for a shared viral community among the vampire bat and cattle populations tested.     

Coronavirus Infection and Diversity in Bats in the Australasian Region

Following the SARS outbreak, extensive surveillance was undertaken globally to detect and identify coronavirus diversity in bats. This study sought to identify the diversity and prevalence of coronaviruses in bats in the Australasian region. We identified four different genotypes of coronavirus, three of which (an alphacoronavirus and two betacoronaviruses) are potentially new species, having less than 90% nucleotide sequence identity with the most closely related described viruses. We did not detect any SARS-like betacoronaviruses, despite targeting rhinolophid bats, the putative natural host taxa. Our findings support the virus-host co-evolution hypothesis, with the detection of Miniopterus bat coronavirus HKU8 (previously reported in Miniopterus species in China, Hong Kong and Bulgaria) in Australian Miniopterus species. Similarly, we detected a novel betacoronavirus genotype from Pteropus alecto which is most closely related to Bat coronavirus HKU9 identified in other pteropodid bats in China, Kenya and the Philippines. We also detected possible cross-species transmission of bat coronaviruses, and the apparent enteric tropism of these viruses. Thus, our findings are consistent with a scenario wherein the current diversity and host specificity of coronaviruses reflects co-evolution with the occasional host shift.     

Analysis of Cathepsin and Furin Proteolytic Enzymes Involved in Viral Fusion Protein Activation in Cells of the Bat Reservoir Host

Bats of different species play a major role in the emergence and transmission of highly pathogenic viruses including Ebola virus, SARS-like coronavirus and the henipaviruses. These viruses require proteolytic activation of surface envelope glycoproteins needed for entry, and cellular cathepsins have been shown to be involved in proteolysis of glycoproteins from these distinct virus families. Very little is currently known about the available proteases in bats. To determine whether the utilization of cathepsins by bat-borne viruses is related to the nature of proteases in their natural hosts, we examined proteolytic processing of several viral fusion proteins in cells derived from two fruit bat species, Pteropus alecto and Rousettus aegyptiacus. Our work shows that fruit bat cells have homologs of cathepsin and furin proteases capable of cleaving and activating both the cathepsin-dependent Hendra virus F and the furin-dependent parainfluenza virus 5 F proteins. Sequence analysis comparing Pteropus alecto furin and cathepsin L to proteases from other mammalian species showed a high degree of conservation; however significant amino acid variation occurs at the C-terminus of Pteropus alecto furin. Further analysis of furin-like proteases from fruit bats revealed that these proteases are catalytically active and resemble other mammalian furins in their response to a potent furin inhibitor. However, kinetic analysis suggests that differences may exist in the cellular localization of furin between different species. Collectively, these results indicate that the unusual role of cathepsin proteases in the life cycle of bat-borne viruses is not due to the lack of active furin-like proteases in these natural reservoir species; however, differences may exist between furin proteases present in fruit bats compared to furins in other mammalian species, and these differences may impact protease usage for viral glycoprotein processing.     

Genetic Characterization of Betacoronavirus Lineage C Viruses in Bats Reveals Marked Sequence Divergence in the Spike Protein of Pipistrellus Bat Coronavirus HKU5 in Japanese Pipistrelle: Implications for the Origin of the Novel Middle East Respiratory Syndrome Coronavirus

While the novel Middle East respiratory syndrome coronavirus (MERS-CoV) is closely related to Tylonycteris bat CoV HKU4 (Ty-BatCoV HKU4) and Pipistrellus bat CoV HKU5 (Pi-BatCoV HKU5) in bats from Hong Kong, and other potential lineage C betacoronaviruses in bats from Africa, Europe, and America, its animal origin remains obscure. To better understand the role of bats in its origin, we examined the molecular epidemiology and evolution of lineage C betacoronaviruses among bats. Ty-BatCoV HKU4 and Pi-BatCoV HKU5 were detected in 29% and 25% of alimentary samples from lesser bamboo bat (Tylonycteris pachypus) and Japanese pipistrelle (Pipistrellus abramus), respectively. Sequencing of their RNA polymerase (RdRp), spike (S), and nucleocapsid (N) genes revealed that MERS-CoV is more closely related to Pi-BatCoV HKU5 in RdRp (92.1% to 92.3% amino acid [aa] identity) but is more closely related to Ty-BatCoV HKU4 in S (66.8% to 67.4% aa identity) and N (71.9% to 72.3% aa identity). Although both viruses were under purifying selection, the S of Pi-BatCoV HKU5 displayed marked sequence polymorphisms and more positively selected sites than that of Ty-BatCoV HKU4, suggesting that Pi-BatCoV HKU5 may generate variants to occupy new ecological niches along with its host in diverse habitats. Molecular clock analysis showed that they diverged from a common ancestor with MERS-CoV at least several centuries ago. Although MERS-CoV may have diverged from potential lineage C betacoronaviruses in European bats more recently, these bat viruses were unlikely to be the direct ancestor of MERS-CoV. Intensive surveillance for lineage C betaCoVs in Pipistrellus and related bats with diverse habitats and other animals in the Middle East may fill the evolutionary gap.     

Identification of an Unclassified Paramyxovirus in Coleura afra: A Potential Case of Host Specificity

Bats are known to harbor multiple paramyxoviruses. Despite the creation of two new genera, Aquaparamyxovirus and Ferlavirus, to accommodate this increasing diversity, several recently isolated or characterized viruses remain unclassified beyond the subfamily level. In the present study, among 985 bats belonging to 6 species sampled in the Belinga caves of Gabon, RNA of an unclassified paramyxovirus (Belinga bat virus, BelPV) was discovered in 14 African sheath-tailed bats (Coleura afra), one of which exhibited several hemorrhagic lesions at necropsy, and viral sequence was obtained in two animals. Phylogenetically, BelPV is related to J virus and Beilong virus (BeiPV), two other unclassified paramyxoviruses isolated from rodents. In the diseased BelPV-infected C. afra individual, high viral load was detected in the heart, and the lesions were consistent with those reported in wild rodents and mice experimentally infected by J virus. BelPV was not detected in other tested bat species sharing the same roosting sites and living in very close proximity with C. afra in the two caves sampled, suggesting that this virus may be host-specific for C. afra. The mode of transmission of this paramyxovirus in bat populations remains to be discovered.     

Bioecological Drivers of Rabies Virus Circulation in a Neotropical Bat Community

Introduction

In addition to the commonly accepted importance of the vampire bat in the maintenance and transmission of the rabies virus (RABV) in South America, RABV infection of other species is widely evidenced, challenging their role in the viral cycle.

Methodology / Principles findings

To identify the bioecological drivers of RABV circulation in neotropical bat communities, we conducted a molecular and serological survey on almost 1,000 bats from 30 species, and a 4-year longitudinal survey in two colonies of vampire bats in French Guiana. RABV was molecularly detected in a common vampire and in a frugivorous bat. The sequences corresponded to haematophagous bat-related strains and were close to viruses circulating in the Brazilian Amazon region. Species’ seroprevalence ranged from 0 to 20%, and the risk of seropositivity was higher in bats with a haematophagous diet, living in monospecific colonies and in dense forests. The longitudinal survey showed substantial temporal fluctuations, with individual waves of seroconversions and waning immunity. The high prevalences observed in bat communities, in most habitats and in species that do not share the same microhabitats and bioecological patterns, the temporal variations, and a rather short period of detectable antibodies as observed in recaptured vampires suggest (i) frequent exposure of animals, (ii) an ability of the infected host to control and eliminate the virus, (iii) more relaxed modes of exposure between bats than the commonly assumed infection via direct contact with saliva of infected animals, all of which should be further investigated.

Conclusions / significance

We hypothesize that RABV circulation in French Guiana is mainly maintained in the pristine forest habitats that may provide sufficient food resources to allow vampire bats, the main prevalent species, to survive and RABV to be propagated. However, on the forest edge and in disturbed areas, human activities may induce more insidious effects such as defaunation. One of the ecological consequences is the disappearance of resources for tertiary or secondary consumers. Populations of vampires may then shift to alternative resources such as cattle, domestic animals and humans. Therefore, a good forest status, allowing both a dilution effect in highly rich bat communities and the maintenance of large populations of medium-sized and large mammals used as prey by vampires, should prevent their migration to anthropized areas.     

Viral metagenomics reveals diverse viruses in the fecal samples of children with diarrhea

Diarrhea is the third leading cause of death in developing countries in children under the age of five. About half a million children die of diarrhea every year, most of which in developing countries. Viruses are the main pathogen of diarrhea. In China, the fecal virome of children with diarrhea has been rarely studied. Using an unbiased viral metagenomics approach, we analyzed the fecal virome in children with diarrhea. Many DNA or RNA viruses associated with diarrhea identified in those fecal samples were mainly from six families of Adenoviridae, Astroviridae, Caliciviridae, Parvoviridae, Picornaviridae, and Reoviridae. Among them, the family of Caliciviridae accounts for the largest proportion of 78.42%, following with Adenoviridae (8.94%) and Picornaviridae (8.36%). In addition to those diarrhea-related viruses that have already been confirmed to cause human diarrhea, the viruses not associated with diarrhea were also identified including anellovirus and picobirnavirus. This study increased our understanding of diarrheic children fecal virome and provided valuable information for the prevention and treatment of viral diarrhea in this area.