Host traits shape virome composition and virus transmission in wild small mammals.

Bats, rodents, and shrews are the most important animal sources of human infectious diseases. However, the evolution and transmission of viruses among them remain largely unexplored. Through the meta-transcriptomic sequencing of internal organ and fecal samples from 2,443 wild bats, rodents, and shrews sampled from four Chinese habitats, we identified 669 viruses, including 534 novel viruses, thereby greatly expanding the mammalian virome. Our analysis revealed high levels of phylogenetic diversity, identified cross-species virus transmission events, elucidated virus origins, and identified cases of invertebrate viruses in mammalian hosts. Host order and sample size were the most important factors impacting virome composition and patterns of virus spillover. Shrews harbored a high richness of viruses, including many invertebrate-associated viruses with multi-organ distributions, whereas rodents carried viruses with a greater capacity for host jumping. These data highlight the remarkable diversity of mammalian viruses in local habitats and their ability to emerge in new hosts.

The evolutionary history of vertebrate RNA viruses.

Our understanding of the diversity and evolution of vertebrate RNA viruses is largely limited to those found in mammalian and avian hosts and associated with overt disease. Here, using a large-scale meta-transcriptomic approach, we discover 214 vertebrate-associated viruses in reptiles, amphibians, lungfish, ray-finned fish, cartilaginous fish and jawless fish. The newly discovered viruses appear in every family or genus of RNA virus associated with vertebrate infection, including those containing human pathogens such as influenza virus, the Arenaviridae and Filoviridae families, and have branching orders that broadly reflected the phylogenetic history of their hosts. We establish a long evolutionary history for most groups of vertebrate RNA virus, and support this by evaluating evolutionary timescales using dated orthologous endogenous virus elements. We also identify new vertebrate-specific RNA viruses and genome architectures, and re-evaluate the evolution of vector-borne RNA viruses. In summary, this study reveals diverse virus-host associations across the entire evolutionary history of the vertebrates.

Author Correction: The evolutionary history of vertebrate RNA viruses.

Change history: In this Article, author Li Liu should be associated with affiliation number 5 (College of Marine Sciences, South China Agricultural University, Guangzhou, Guangdong, China), rather than affiliation number 4 (Wenzhou Center for Disease Control and Prevention, Wenzhou, Zhejiang, China). This has been corrected online.

Redefining the invertebrate RNA virosphere.

Current knowledge of RNA virus biodiversity is both biased and fragmentary, reflecting a focus on culturable or disease-causing agents. Here we profile the transcriptomes of over 220 invertebrate species sampled across nine animal phyla and report the discovery of 1,445 RNA viruses, including some that are sufficiently divergent to comprise new families. The identified viruses fill major gaps in the RNA virus phylogeny and reveal an evolutionary history that is characterized by both host switching and co-divergence. The invertebrate virome also reveals remarkable genomic flexibility that includes frequent recombination, lateral gene transfer among viruses and hosts, gene gain and loss, and complex genomic rearrangements. Together, these data present a view of the RNA virosphere that is more phylogenetically and genomically diverse than that depicted in current classification schemes and provide a more solid foundation for studies in virus ecology and evolution.

Outbreak of COVID-19 in a family, Wenzhou, China.

Since December 2019, China has experienced a widespread outbreak of COVID-19. However, at the early stage of outbreak, investigations revealed a variety of patterns resulting in the transmission of COVID-19. Thus, it is essential to understand the transmission types and the potential for sustained human-to-human transmission. Moreover, the information regarding the characteristics of transmission helps in coordinating the current screening programme, and controlling and containing measures, and also, helps in deciding the appropriate quarantine duration. Thus, this investigation reports an outbreak of COVID-19 in a family residing in Wenzhou, Zhejiang, China during the month of January-February 2020.

Discovery of a Highly Divergent Coronavirus in the Asian House Shrew from China Illuminates the Origin of the Alphacoronaviruses.

Although shrews are one of the largest groups of mammals, little is known about their role in the evolution and transmission of viral pathogens, including coronaviruses (CoVs). We captured 266 Asian house shrews (Suncus murinus) in Jiangxi and Zhejiang Provinces, China, during 2013 to 2015. CoV RNA was detected in 24 Asian house shrews, with an overall prevalence of 9.02%. Complete viral genome sequences were successfully recovered from the RNA-positive samples. The newly discovered shrew CoV fell into four lineages reflecting their geographic origins, indicative of largely allopatric evolution. Notably, these viruses were most closely related to alphacoronaviruses but sufficiently divergent that they should be considered a novel member of the genus Alphacoronavirus, which we denote Wénchéng shrew virus (WESV). Phylogenetic analysis revealed that WESV was a highly divergent member of the alphacoronaviruses and, more dramatically, that the S gene of WESV fell in a cluster that was genetically distinct from that of known coronaviruses. The divergent position of WESV suggests that coronaviruses have a long association with Asian house shrews. In addition, the genome of WESV contains a distinct NS7 gene that exhibits no sequence similarity to genes of any known viruses. Together, these data suggest that shrews are natural reservoirs for coronaviruses and may have played an important and long-term role in CoV evolution.IMPORTANCE The subfamily Coronavirinae contains several notorious human and animal pathogens, including severe acute respiratory syndrome coronavirus, Middle East respiratory syndrome coronavirus, and porcine epidemic diarrhea virus. Because of their genetic diversity and phylogenetic relationships, it has been proposed that the alphacoronaviruses likely have their ultimate ancestry in the viruses residing in bats. Here, we describe a novel alphacoronavirus (Wénchéng shrew virus [WESV]) that was sampled from Asian house shrews in China. Notably, WESV is a highly divergent member of the alphacoronaviruses and possesses an S gene that is genetically distinct from those of all known coronaviruses. In addition, the genome of WESV contains a distinct NS7 gene that exhibits no sequence similarity to those of any known viruses. Together, these data suggest that shrews are important and longstanding hosts for coronaviruses that merit additional research and surveillance.

Divergent Viruses Discovered in Arthropods and Vertebrates Revise the Evolutionary History of the Flaviviridae and Related Viruses.

UNLABELLED: Viruses of the family Flaviviridae are important pathogens of humans and other animals and are currently classified into four genera. To better understand their diversity, evolutionary history, and genomic flexibility, we used transcriptome sequencing (RNA-seq) to search for the viruses related to the Flaviviridae in a range of potential invertebrate and vertebrate hosts. Accordingly, we recovered the full genomes of five segmented jingmenviruses and 12 distant relatives of the known Flaviviridae ("flavi-like" viruses) from a range of arthropod species. Although these viruses are highly divergent, they share a similar genomic plan and common ancestry with the Flaviviridae in the NS3 and NS5 regions. Remarkably, although these viruses fill in major gaps in the phylogenetic diversity of the Flaviviridae, genomic comparisons reveal important changes in genome structure, genome size, and replication/gene regulation strategy during evolutionary history. In addition, the wide diversity of flavi-like viruses found in invertebrates, as well as their deep phylogenetic positions, suggests that they may represent the ancestral forms from which the vertebrate-infecting viruses evolved. For the vertebrate viruses, we expanded the previously mammal-only pegivirus-hepacivirus group to include a virus from the graceful catshark (Proscyllium habereri), which in turn implies that these viruses possess a larger host range than is currently known. In sum, our data show that the Flaviviridae infect a far wider range of hosts and exhibit greater diversity in genome structure than previously anticipated. IMPORTANCE: The family Flaviviridae of RNA viruses contains several notorious human pathogens, including dengue virus, West Nile virus, and hepatitis C virus. To date, however, our understanding of the biodiversity and evolution of the Flaviviridae has largely been directed toward vertebrate hosts and their blood-feeding arthropod vectors. Therefore, we investigated an expanded group of potential arthropod and vertebrate host species that have generally been ignored by surveillance programs. Remarkably, these species contained diverse flaviviruses and related viruses that are characterized by major changes in genome size and genome structure, such that these traits are more flexible than previously thought. More generally, these data suggest that arthropods may be the ultimate reservoir of the Flaviviridae and related viruses, harboring considerable genetic and phenotypic diversity. In sum, this study revises the traditional view on the evolutionary history, host range, and genomic structures of a major group of RNA viruses.

A tick-borne segmented RNA virus contains genome segments derived from unsegmented viral ancestors.

Although segmented and unsegmented RNA viruses are commonplace, the evolutionary links between these two very different forms of genome organization are unclear. We report the discovery and characterization of a tick-borne virus--Jingmen tick virus (JMTV)--that reveals an unexpected connection between segmented and unsegmented RNA viruses. The JMTV genome comprises four segments, two of which are related to the nonstructural protein genes of the genus Flavivirus (family Flaviviridae), whereas the remaining segments are unique to this virus, have no known homologs, and contain a number of features indicative of structural protein genes. Remarkably, homology searching revealed that sequences related to JMTV were present in the cDNA library from Toxocara canis (dog roundworm; Nematoda), and that shared strong sequence and structural resemblances. Epidemiological studies showed that JMTV is distributed in tick populations across China, especially Rhipicephalus and Haemaphysalis spp., and experiences frequent host-switching and genomic reassortment. To our knowledge, JMTV is the first example of a segmented RNA virus with a genome derived in part from unsegmented viral ancestors.

Diversity and evolution of avian influenza viruses in live poultry markets, free-range poultry and wild wetland birds in China.

The wide circulation of novel avian influenza viruses (AIVs) highlights the risk of pandemic influenza emergence in China. To investigate the prevalence and genetic diversity of AIVs in different ecological contexts, we surveyed AIVs in live poultry markets (LPMs), free-range poultry and the wetland habitats of wild birds in Zhejiang and Hubei provinces. Notably, LPMs contained the highest frequency of AIV infection, and the greatest number of subtypes (n = 9) and subtype co-infections (n = 14), as well as frequent reassortment, suggesting that they play an active role in fuelling AIV transmission. AIV-positive samples were also identified in wild birds in both provinces and free-range poultry in one sampling site close to a wetland region in Hubei. H9N2, H7N9 and H5N1 were the most commonly sampled subtypes in the LPMs from Zhejiang, whilst H5N6 and H9N2 were the dominant subtypes in the LPMs from Hubei. Phylogenetic analyses of the whole-genome sequences of 43 AIVs revealed that three reassortant H5 subtypes were circulating in LMPs in both geographical regions. Notably, the viruses sampled from the wetland regions and free-range poultry contained complex reassortants, for which the origins of some segments were unclear. Overall, our study highlights the extent of AIV genetic diversity in two highly populated parts of central and south-eastern China, particularly in LPMs, and emphasizes the need for continual surveillance.

Taxonomy of the order Mononegavirales: update 2016.

In 2016, the order Mononegavirales was emended through the addition of two new families (Mymonaviridae and Sunviridae), the elevation of the paramyxoviral subfamily Pneumovirinae to family status (Pneumoviridae), the addition of five free-floating genera (Anphevirus, Arlivirus, Chengtivirus, Crustavirus, and Wastrivirus), and several other changes at the genus and species levels. This article presents the updated taxonomy of the order Mononegavirales as now accepted by the International Committee on Taxonomy of Viruses (ICTV).

Reemergence and Autochthonous Transmission of Dengue Virus, Eastern China, 2014.

In 2014, 20 dengue cases were reported in the cities of Wenzhou (5 cases) and Wuhan (15 cases), China, where dengue has rarely been reported. Dengue virus 1 was detected in 4 patients. Although most of these cases were likely imported, epidemiologic analysis provided evidence for autochthonous transmission.

Phylogeny and origins of hantaviruses harbored by bats, insectivores, and rodents.

Hantaviruses are among the most important zoonotic pathogens of humans and the subject of heightened global attention. Despite the importance of hantaviruses for public health, there is no consensus on their evolutionary history and especially the frequency of virus-host co-divergence versus cross-species virus transmission. Documenting the extent of hantavirus biodiversity, and particularly their range of mammalian hosts, is critical to resolving this issue. Here, we describe four novel hantaviruses (Huangpi virus, Lianghe virus, Longquan virus, and Yakeshi virus) sampled from bats and shrews in China, and which are distinct from other known hantaviruses. Huangpi virus was found in Pipistrellus abramus, Lianghe virus in Anourosorex squamipes, Longquan virus in Rhinolophus affinis, Rhinolophus sinicus, and Rhinolophus monoceros, and Yakeshi virus in Sorex isodon, respectively. A phylogenetic analysis of the available diversity of hantaviruses reveals the existence of four phylogroups that infect a range of mammalian hosts, as well as the occurrence of ancient reassortment events between the phylogroups. Notably, the phylogenetic histories of the viruses are not always congruent with those of their hosts, suggesting that cross-species transmission has played a major role during hantavirus evolution and at all taxonomic levels, although we also noted some evidence for virus-host co-divergence. Our phylogenetic analysis also suggests that hantaviruses might have first appeared in Chiroptera (bats) or Soricomorpha (moles and shrews), before emerging in rodent species. Overall, these data indicate that bats are likely to be important natural reservoir hosts of hantaviruses.

Cross-species transmission in the speciation of the currently known murinae-associated hantaviruses.

To gain more insight into the phylogeny of Dabieshan virus (DBSV), carried by Niviventer confucianus and other Murinae-associated hantaviruses, genome sequences of novel variants of DBSV were recovered from Niviventer rats trapped in the mountainous areas of Wenzhou, China. Genetic analyses show that all known genetic variants of DBSV, including the ones identified in this study, are distinct from other Murinae-associated hantaviruses. DBSV variants show geographic clustering and high intraspecies diversity. The data suggest that DBSV is a distinct species in the genus Hantavirus. Interestingly, DBSV shows the highest sequence identity to Hantaan virus (HTNV), with a >7% difference in the sequences of the N, GPC, and L proteins, while N. confucianus is more closely related to Rattus norvegicus (the host of Seoul virus [SEOV]) than to Apodemus agrarius (the host of HTNV and Saaremaa virus [SAAV]). Further genetic analyses of all known Murinae-associated hantaviruses (both established and tentative species) show that many of them, including DBSV, may have originated from host switching. The estimation of evolutionary rates and divergence time supports the role of cross-species transmission in the evolution of Murinae-associated hantaviruses. The detection of positive selection suggests that genetic drift may contribute to the speciation of Murinae-associated hantaviruses and that adaptation has a role as well.

Migration of Norway rats resulted in the worldwide distribution of Seoul hantavirus today.

Despite the worldwide distribution, most of the known Seoul viruses (SEOV) are closely related to each other. In this study, the M and the S segment sequences of SEOV were recovered from 130 lung tissue samples (mostly of Norway rats) and from six patient serum samples by reverse transcription-PCR. Genetic analysis revealed that all sequences belong to SEOV and represent 136 novel strains. Phylogenetic analysis of all available M and S segment sequences of SEOV, including 136 novel Chinese strains, revealed four distinct groups. All non-Chinese SEOV strains and most of the Chinese variants fell into the phylogroup A, while the Chinese strains originating from mountainous areas clustered into three other distinct groups (B, C, and D). We estimated that phylogroup A viruses may have arisen only within the last several centuries. All non-Chinese variants appeared to be directly originated from China. Thus, phylogroup A viruses distributed worldwide may share a recent ancestor, whereas SEOV seems to be as diversified genetically as other hantaviruses. In addition, all available mitochondrial DNA (mtDNA) sequences of Norway rats, including our 44 newly recovered mtDNA sequences, were divided into two phylogenetic groups. The first group, which is associated with the group A SEOV variants, included most of rats from China and also all non-Chinese rats, while the second group consisted of a few rats originating only from mountain areas in China. We hypothesize that an ancestor of phylogroup A SEOV variants was first exported from China to Europe and then spread through the New World following the migration of Norway rats.

Coronaviruses in wild animals sampled in and around Wuhan at the beginning of COVID-19 emergence.

Over the last several decades, no emerging virus has had a profound impact on the world as the SARS-CoV-2 that emerged at the end of 2019 has done. To know where severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) originated from and how it jumped into human population, we immediately started a surveillance investigation in wild mammals in and around Wuhan when we determined the agent. Herein, coronaviruses were screened in the lung, liver, and intestinal tissue samples from fifteen raccoon dogs, seven Siberian weasels, three hog badgers, and three Reeves's muntjacs collected in Wuhan and 334 bats collected around Wuhan. Consequently, eight alphacoronaviruses were identified in raccoon dogs, while nine betacoronaviruses were found in bats. Notably, the newly discovered alphacoronaviruses shared a high whole-genome sequence similarity (97.9 per cent) with the canine coronavirus (CCoV) strain 2020/7 sampled from domestic dog in the UK. Some betacoronaviruses identified here were closely related to previously known bat SARS-CoV-related viruses sampled from Hubei province and its neighbors, while the remaining betacoronaviruses exhibited a close evolutionary relationship with SARS-CoV-related bat viruses in the RdRp gene tree and clustered together with SARS-CoV-2-related bat coronaviruses in the M, N and S gene trees, but with relatively low similarity. Additionally, these newly discovered betacoronaviruses seem unlikely to bind angiotensin-converting enzyme 2 because of the deletions in the two key regions of their receptor-binding motifs. Finally, we did not find SARS-CoV-2 or its progenitor virus in these animal samples. Due to the high circulation of CCoVs in raccoon dogs in Wuhan, more scientific efforts are warranted to better understand their diversity and evolution in China and the possibility of a potential human agent.

RNA viromes from terrestrial sites across China expand environmental viral diversity.

Environmental RNA viruses are ubiquitous and diverse, and probably have important ecological and biogeochemical impacts. Understanding the global diversity of RNA viruses is limited by sampling biases, dependence on cell culture and PCR for virus discovery, and a focus on viruses pathogenic to humans or economically important animals and plants. To address this knowledge gap, we generated metatranscriptomic sequence data from 32 diverse environments in 16 provinces and regions of China. We identified 6,624 putatively novel virus operational taxonomic units from soil, sediment and faecal samples, greatly expanding known diversity of the RNA virosphere. These newly identified viruses included positive-sense, negative-sense and double-strand RNA viruses from at least 62 families. Sediments and animal faeces were rich sources of viruses. Virome compositions were affected by local environmental factors, including organic content and eukaryote species abundance. Notably, environmental factors had a greater impact on the abundance and diversity of plant, fungal and bacterial viruses than of animal viromes. Our data confirm that RNA viruses are an integral part of both terrestrial and aquatic ecosystems.

Hantaviruses in small mammals and humans in the coastal region of Zhejiang Province, China.

To gain more insights into the epidemiology of hantaviruses in the coastal region of Zhejiang Province, China, the morbidity and mortality of hemorrhagic fever with renal syndrome (HFRS) were analyzed in two coastal areas: Cixi (hilly terrain) and Wenzhou (mountainous terrain). More HFRS cases have been reported in Cixi than in Wenzhou. Annual incidence rate of HFRS in Cixi had been on the level of approximately 1.5/100,000 from 1968 (when the first HFRS case was reported) to 2007, with the highest incidence rate of 8.54/100,000 in 1999. The annual incidence rate in Wenzhou has been relatively low, less than 0.5/100,000 since the first HFRS case was reported in 1981. A total of 461 rodents and 199 shrews were captured in these two areas. Hantavirus antibodies were detected in 16 of 241 (6.64%) Rattus norvegicus and 13 of 122 (10.66%) R. flavipectus. Interestingly, hantavirus antigens were identified in 6 of 196 (3.06%) Suncus murinus. Genetic analysis showed that partial M and S segment sequences recovered from rats in the two regions belong to Seoul virus (SEOV) and can be assigned into two genetic lineages. SEOV variants that belong to these two lineages of viruses are distributed widely in China and have been found outside China. As most trapped rodents were rats and SEOV was the only hantavirus detected, these results suggested that SEOV plays an important role in human hantavirus infections. They also reinforce the need for vigilance in preventing HFRS caused by hantaviruses in the coastal region.

Extensive genetic diversity and host range of rodent-borne coronaviruses.

To better understand the genetic diversity, host associations and evolution of coronaviruses (CoVs) in China we analyzed a total of 696 rodents encompassing 16 different species sampled from Zhejiang and Yunnan provinces. Based on reverse transcriptase PCR-based CoV screening of fecal samples and subsequent sequence analysis of the RNA-dependent RNA polymerase gene, we identified CoVs in diverse rodent species, comprising Apodemus agrarius, Apodemus chevrieri, Apodemus latronum, Bandicota indica, Eothenomys cachinus, Eothenomys miletus, Rattus andamanensis, Rattus norvegicus, and Rattus tanezumi. CoVs were particularly commonplace in A. chevrieri, with a detection rate of 12.44 per cent (24/193). Genetic and phylogenetic analysis revealed the presence of three groups of CoVs carried by a range of rodents that were closely related to the Lucheng Rn rat CoV (LRNV), China Rattus CoV HKU24 (ChRCoV_HKU24), and Longquan Rl rat CoV (LRLV) identified previously. One newly identified A. chevrieri-associated virus closely related to LRNV lacked an NS2 gene. This virus had a similar genetic organization to AcCoV-JC34, recently discovered in the same rodent species in Yunnan, suggesting that it represents a new viral subtype. Notably, additional variants of LRNV were identified that contained putative non-structural (NS)2b genes located downstream of the NS2 gene that were likely derived from the host genome. Recombination events were also identified in the open reading frame (ORF) 1a gene of Lijiang-71. In sum, these data reveal the substantial genetic diversity and genomic complexity of rodent-borne CoVs, and extend our knowledge of these major wildlife virus reservoirs.

Diverse and abundant resistome in terrestrial and aquatic vertebrates revealed by transcriptional analysis.

Despite increasing evidence that antibiotic resistant pathogens are shared among humans and animals, the diversity, abundance and patterns of spread of antibiotic resistance genes (ARGs) in wildlife remains unclear. We identified 194 ARGs associated with phenotypic resistance to 13 types of antibiotic in meta-transcriptomic data generated from a broad range of lower vertebrates residing in both terrestrial and aquatic habitats. These ARGs, confirmed by PCR, included those that shared high sequence similarity to clinical isolates of public health concern. Notably, the lower vertebrate resistome varied by ecological niche of the host sampled. The resistomes in marine fish shared high similarity and were characterized by very high abundance, distinct from that observed in other habitats. An assessment of ARG mobility found that ARGs in marine fish were frequently co-localized with mobile elements, indicating that they were likely spread by horizontal gene transfer. Together, these data reveal the remarkable diversity and transcriptional levels of ARGs in lower vertebrates, and suggest that these wildlife species might play an important role in the global spread of ARGs.

Newly identified viral genomes in pangolins with fatal disease.

Epizootic pathogens pose a major threat to many wildlife species, particularly in the context of rapidly changing environments. Pangolins (order Pholidota) are highly threatened mammals, in large part due to the trade in illegal wildlife. During July to August 2018 four sick wild pangolins (three Manis javanica and one Manis pentadactyla) exhibiting a variety of clinical symptoms were rescued by the Jinhua Wildlife Protection Station in Zhejiang province, China. Although three of these animals died, fortunately one recovered after 2 weeks of symptomatic treatment. Using meta-transcriptomics combined with reverse transcription polymerase chain reaction (RT-PCR), we identified two novel RNA viruses in two of the dead pangolins. Genomic analysis revealed that these viruses were most closely related to pestiviruses and coltiviruses, although still highly genetically distinct, with more than 48 and 25 per cent sequence divergence at the amino acid level, respectively. We named these Dongyang pangolin virus (DYPV) and Lishui pangolin virus (LSPV) based on the sampling site and hosts. Although coltiviruses (LSPV) are known to be transmitted by ticks, we found no evidence of LSPV in ticks sampled close to where the pangolins were collected. In addition, although DYPV was present in nymph ticks (Amblyomma javanense) collected from a diseased pangolin, they were not found in the local tick population. Epidemiological investigation revealed that both novel viruses might have been imported following the illegal international trade of pangolins. Hence, these data indicate that illegal wildlife trafficking not only threatens the status of pangolin populations, but may also spread epizootic pathogens.