Search

Abstract Over the last decade metagenomic sequencing has facilitated an increasing number of virome-scale studies, leading to an exponential expansion in understanding of virus diversity. This is partially driven by the decreasing costs of metagenomic sequencing, improvements in computational tools for revealing novel viruses, and an increased understanding of the key role that viruses play in human and animal health. A central concern associated with this remarkable increase in the number of virome-scale studies is the lack of broadly accepted “gold standards” for reporting the data and results generated. This is of particular importance for animal virome studies as there are a multitude of nuanced approaches for both data presentation and analysis, all of which impact the resulting outcomes. As such, the results of published studies can be difficult to contextualise and may be of reduced utility due to reporting deficiencies. Herein, we aim to address these reporting issues by outlining recommendations for the presentation of virome data, encouraging a transparent communication of findings that can be interpreted in evolutionary and ecological contexts.

Abstract Despite a rapid expansion of Porcine reproductive and respiratory syndrome virus (PRRSV) sublineage 8.7 over recent years, very little is known about the patterns of virus evolution, dispersal, and the factors influencing this dispersal. Relying on a national PRRSV surveillance project established over 20 years ago, we expand the available genomic data of sublineage 8.7 from China. We perform independent interlineage and intralineage recombination analyses for the entire study period, which showed a heterogeneous recombination pattern. A series of Bayesian phylogeographic analyses uncover the role of Guangdong as an important infection hub within Asia. The spatial spread of PRRSV is highly linked with a composite of human activities and the heterogeneous provincial distribution of the swine industry, largely propelled by the smaller-scale Chinese rural farming systems in the past years. We sequence all four available modified live vaccines (MLVs) and perform genomic analyses with publicly available data, of which our results suggest a key “leaky” period spanning 2011–2017 with two concurrent amino acid mutations in ORF1a 957 and ORF2 250. Overall, our study provides an in-depth overview of the evolution, transmission dynamics, and potential leaky status of HP-PRRS MLVs, providing critical insights into new MLV development.

ABSTRACT Companion animals such as cats and dogs harbor diverse microbial communities that can potentially impact human health due to close and frequent contact. To better characterize their total infectomes and assess zoonotic risks, we characterized the overall infectomes of companion animals (cats and dogs) and evaluated their potential zoonotic risks. Meta-transcriptomic analyses were performed on 239 samples from cats and dogs collected across China, identifying 24 viral species, 270 bacterial genera, and two fungal genera. Differences in the overall microbiome and infectome composition were compared across different animal species (cats or dogs), sampling sites (rectal or oropharyngeal), and health status (healthy or diseased). Diversity analyses revealed that viral abundance was generally higher in diseased animals compared to healthy ones, while differences in microbial composition were mainly driven by sampling site, followed by animal species and health status. Disease association analyses validated the pathogenicity of known pathogens and suggested potential pathogenic roles of previously undescribed bacteria and newly discovered viruses. Cross-species transmission analyses identified seven pathogens shared between cats and dogs, such as alphacoronavirus 1, which was detected in both oropharyngeal and rectal swabs albeit with differential pathogenicity. Further analyses showed that some viruses, like alphacoronavirus 1, harbored multiple lineages exhibiting distinct pathogenicity, tissue, or host preferences. Ultimately, a systematic evolutionary screening identified 27 potential zoonotic pathogens in this sample set, with far more bacterial than viral species, implying potential health threats to humans. Overall, our meta-transcriptomic analysis reveals a landscape of actively transcribing microorganisms in major companion animals, highlighting key pathogens, those with the potential for cross-species transmission, and possible zoonotic threats.IMPORTANCEThis study provides a comprehensive characterization of the entire community of infectious microbes (viruses, bacteria, and fungi) in companion animals like cats and dogs, termed the “infectome.” By analyzing hundreds of samples from across China, the researchers identified numerous known and novel pathogens, including 27 potential zoonotic agents that could pose health risks to both animals and humans. Notably, some of these zoonotic pathogens were detected even in apparently healthy pets, highlighting the importance of surveillance. The study also revealed key microbial factors associated with respiratory and gastrointestinal diseases in pets, as well as potential cross-species transmission events between cats and dogs. Overall, this work sheds light on the complex microbial landscapes of companion animals and their potential impacts on animal and human health, underscoring the need for monitoring and management of these infectious agents.

c-FLIP functions as a dual regulator of apoptosis and inflammation, yet its implications in Zika virus (ZIKV) infection remain partially understood, especially in the context of ZIKV-induced congenital Zika syndrome (CZS) where both apoptosis and inflammation play pivotal roles. Our findings demonstrate that c-FLIP promotes ZIKV infection in placental cells and myeloid-derived macrophages, involving inflammation and caspase-8/3-mediated apoptosis. Moreover, our observations reveal that c-FLIP augments ZIKV infection in multiple tissues, including blood cell, spleen, uterus, testis, and the brain of mice. Notably, the partial deficiency of c-FLIP provides protection to embryos against ZIKV-induced CZS, accompanied by a reduction in caspase-3-mediated apoptosis. Additionally, we have found a distinctive parental effect of c-FLIP influencing ZIKV replication in fetal heads. In summary, our study reveals the critical role of c-FLIP as a positive regulator in caspase-8/3-mediated apoptosis during ZIKV infection, significantly contributing to the development of CZS.

Abstract Bats are reservoir hosts for many zoonotic viruses. Despite this, relatively little is known about the diversity and abundance of viruses within individual bats, and hence the frequency of virus co-infection and spillover among them. We characterize the mammal-associated viruses in 149 individual bats sampled from Yunnan province, China, using an unbiased meta-transcriptomics approach. This reveals a high frequency of virus co-infection (simultaneous infection of bat individuals by multiple viral species) and spillover among the animals studied, which may in turn facilitate virus recombination and reassortment. Of note, we identify five viral species that are likely to be pathogenic to humans or livestock, based on phylogenetic relatedness to known pathogens or in vitro receptor binding assays. This includes a novel recombinant SARS-like coronavirus that is closely related to both SARS-CoV and SARS-CoV-2. In vitro assays indicate that this recombinant virus can utilize the human ACE2 receptor such that it is likely to be of increased emergence risk. Our study highlights the common occurrence of co-infection and spillover of bat viruses and their implications for virus emergence.

Emerging zoonoses of wildlife origin caused by previously unknown agents are one of the most important challenges for human health. The Qinghai-Tibet Plateau represents a unique ecological niche with diverse wildlife that harbours several human pathogens and numerous previously uncharacterized pathogens. In this study, we identified and characterized a novel arenavirus (namely, plateau pika virus, PPV) from plateau pikas (Ochotona curzoniae) on the Qinghai-Tibet Plateau by virome analysis. Isolated PPV strains could replicate in several mammalian cells. We further investigated PPV pathogenesis using animal models. PPV administered via an intraventricular route caused trembling and sudden death in IFNαβR-/- mice, and pathological inflammatory lesions in brain tissue were observed. According to a retrospective serological survey in the geographical region where PPV was isolated, PPV-specific IgG antibodies were detected in 8 (2.4%) of 335 outpatients with available sera. Phylogenetic analyses revealed that this virus was clearly separated from previously reported New and Old World mammarenaviruses. Under the co-speciation framework, the estimated divergence time of PPV was 77–88 million years ago (MYA), earlier than that of OW and NW mammarenaviruses (26-34 MYA).

ABSTRACT Porcine reproductive and respiratory syndrome virus (PRRSV) is an important swine pathogen affecting the global swine industry. The first open reading frame 5 (ORF5)-based genetic lineage classification system describing global PRRSV-2 genetic diversity was introduced over a decade ago. Although refinements have been proposed for the predominant lineage in the USA (lineage 1), PRRSV-2 phylogenetic classification system at international levels has not been thoroughly evaluated and updated since 2010. In this study, based on analysis of 82,237 global ORF5 sequences reported during 1989–2021, we classified PRRSV-2 into 11 genetic lineages (L1‒L11) and 21 sublineages (L1A‒L1F, L1H‒L1J, L5A‒L5B, L8A‒L8E, and L9A‒L9E). The proposed classification system is flexible for growth if additional lineages, sublineages, or more granular classifications are needed. For example, for more fine-scale epidemiological investigation, L1C was further divided into five groups (L1C.1‒L1C.5), with L1C.5 corresponding to the recently emerged L1C variant. Comparison between restriction fragment length polymorphism (RFLP) typing and phylogenetic classification revealed the inaccuracy of using RFLP to determine PRRSV-2 genetic relatedness in most scenarios. Genetic homology of six commercial PRRSV-2 vaccines to each lineage/sublineage and detection frequency of vaccine-like viruses were determined. Global geographic distribution of each lineage/sublineage was presented. Temporal dynamic changes of PRRSV-2 in the USA during 1989–2021 were investigated by analyzing 73,092 ORF5 sequences. In summary, this study refined PRRSV-2 ORF5-based phylogenetic classification and investigated the geographic distribution and temporal changes of PRRSV-2 at the lineage/sublineage levels. The refined classification system and molecular epidemiology data in this study will be invaluable for future characterization of PRRSV-2. IMPORTANCE In this study, comprehensive analysis of 82,237 global porcine reproductive and respiratory syndrome virus type 2 (PRRSV-2) open reading frame 5 sequences spanning from 1989 to 2021 refined PRRSV-2 genetic classification system, which defines 11 lineages and 21 sublineages and provides flexibility for growth if additional lineages, sublineages, or more granular classifications are needed in the future. Geographic distribution and temporal changes of PRRSV-2 were investigated in detail. This is a thorough study describing the molecular epidemiology of global PRRSV-2. In addition, the reference sequences based on the refined genetic classification system are made available to the public for future epidemiological and diagnostic applications worldwide. The data from this study will facilitate global standardization and application of PRRSV-2 genetic classification.

Ticks, an arthropod known for transmitting various pathogens such as viruses, bacteria, and fungi, pose a perpetual public health concern. A total of 2,570 ticks collected from Nujiang Prefecture in Yunnan Province between 2017 and 2022 were included in the study. Through the meta-transcriptomic sequencing of four locally distributed tick species, we identified 13 RNA viruses belonging to eight viral families, namely, Phenuiviridae, Nairoviridae, Peribunyaviridae, Flaviviridae, Chuviridae, Rhabdoviridae, Orthomyxoviridae, and Totiviridae. The most prevalent viruses were members of the order Bunyavirales, including three of Phenuiviridae, two were classified as Peribunyaviridae, and one was associated with Nairoviridae. However, whether they pose a threat to human health still remains unclear. Indeed, this study revealed the genetic diversity of tick species and tick-borne viruses in Nujiang Prefecture based on COI gene and tick-borne virus research. These data clarified the genetic evolution of some RNA viruses and furthered our understanding of the distribution pattern of tick-borne pathogens, highlighting the importance and necessity of monitoring tick-borne pathogens.

Abstract Background Viruses play critical roles in the marine environment because of their interactions with an extremely broad range of potential hosts. Many studies of viruses in seawater have been published, but viruses that inhabit marine animals have been largely neglected. Oysters are keystone species in coastal ecosystems, yet as filter-feeding bivalves with very large roosting numbers and species co-habitation, it is not clear what role they play in marine virus transmission and coastal microbiome regulation. Results Here, we report a Dataset of Oyster Virome (DOV) that contains 728,784 nonredundant viral operational taxonomic unit contigs (≥ 800 bp) and 3473 high-quality viral genomes, enabling the first comprehensive overview of both DNA and RNA viral communities in the oyster Crassostrea hongkongensis. We discovered tremendous diversity among novel viruses that inhabit this oyster using multiple approaches, including reads recruitment, viral operational taxonomic units, and high-quality virus genomes. Our results show that these viruses are very different from viruses in the oceans or other habitats. In particular, the high diversity of novel circoviruses that we found in the oysters indicates that oysters may be potential hotspots for circoviruses. Notably, the viruses that were enriched in oysters are not random but are well-organized communities that can respond to changes in the health state of the host and the external environment at both compositional and functional levels. Conclusions In this study, we generated a first “knowledge landscape” of the oyster virome, which has increased the number of known oyster-related viruses by tens of thousands. Our results suggest that oysters provide a unique habitat that is different from that of seawater, and highlight the importance of filter-feeding bivalves for marine virus exploration as well as their essential but still invisible roles in regulating marine ecosystems. Video Abstract

Dynamic changes of the paired heavy and light chain B cell receptor (BCR) repertoire provide an essential insight into understanding the humoral immune response post-SARS-CoV-2 infection and vaccination. However, differences between the endogenous paired BCR repertoire kinetics in SARS-CoV-2 infection and previously recovered/naïve subjects treated with the inactivated vaccine remain largely unknown. We performed single-cell V(D)J sequencing of B cells from six healthy donors with three shots of inactivated SARS-CoV-2 vaccine (BBIBP-CorV), five people who received the BBIBP-CorV vaccine after having recovered from COVID-19, five unvaccinated COVID-19 recovered patients and then integrated with public data of B cells from four SARS-CoV-2-infected subjects. We discovered that BCR variable (V) genes were more prominently used in the SARS-CoV-2 exposed groups (both in the group with active infection and in the group that had recovered) than in the vaccinated groups. The VH gene that expanded the most after SARS-CoV-2 infection was IGHV3-33, while IGHV3-23 in the vaccinated groups. SARS-CoV-2-infected group enhanced more BCR clonal expansion and somatic hypermutation than the vaccinated healthy group. A small proportion of public clonotypes were shared between the SARS-CoV-2 infected, vaccinated healthy, and recovered groups. Moreover, several public antibodies had been identified against SARS-CoV-2 spike protein. We comprehensively characterize the paired heavy and light chain BCR repertoire from SARS-CoV-2 infection to vaccination, providing further guidance for the development of the next-generation precision vaccine.

Phlebotomus chinensis sandfly is a neglected insect vector in China that is well-known for carrying Leishmania. Recent studies have expanded its pathogen repertoire with two novel arthropod-borne phleboviruses capable of infecting humans and animals. Despite these discoveries, our knowledge of the general pathogen diversity and overall microbiome composition of this vector species is still very limited. Here we carried out a meta-transcriptomics analysis that revealed the actively replicating/transcribing RNA viruses, DNA viruses, bacteria, and eukaryotic microbes, namely, the “total microbiome”, of several sandfly populations in China. Strikingly, “microbiome” made up 1.8% of total non-ribosomal RNA and comprised more than 87 species, among which 70 were novel, including divergent members of the genera Flavivirus and of the family Trypanosomatidae. Importantly, among these microbes we were able to reveal four distinguished types of human and/or mammalian pathogens, including two phleboviruses (hedi and wuxiang viruses), one novel Spotted fever group rickettsia, as well as a member of Leishmania donovani complex, among which hedi virus and Leishmania each had > 50% pool prevalence rate and relatively high abundance levels. Our study also showed the ubiquitous presence of an endosymbiont, namely Wolbachia, although no anti-viral or anti-pathogen effects were detected based on our data. In summary, our results uncovered the much un-explored diversity of microbes harboured by sandflies in China and demonstrated that high pathogen diversity and abundance are currently present in multiple populations, implying disease potential for exposed local human population or domestic animals.

ABSTRACT Nelson Bay reovirus (NBV) is an emerging zoonotic virus that can cause acute respiratory disease in humans. These viruses are mainly discovered in Oceania, Africa, and Asia, and bats have been identified as their main animal reservoir. However, despite recent expansion of diversity for NBVs, the transmission dynamics and evolutionary history of NBVs are still unclear. This study successfully isolated two NBV strains (MLBC1302 and MLBC1313) from blood-sucking bat fly specimens (Eucampsipoda sundaica) and one (WDBP1716) from the spleen specimen of a fruit bat (Rousettus leschenaultii), which were collected at the China-Myanmar border area of Yunnan Province. Syncytia cytopathic effects (CPE) were observed in BHK-21 and Vero E6 cells infected with the three strains at 48 h postinfection. Electron micrographs of ultrathin sections showed numerous spherical virions with a diameter of approximately 70 nm in the cytoplasm of infected cells. The complete genome nucleotide sequence of the viruses was determined by metatranscriptomic sequencing of infected cells. Phylogenetic analysis demonstrated that the novel strains were closely related to Cangyuan orthoreovirus, Melaka orthoreovirus, and human-infecting Pteropine orthoreovirus HK23629/07. Simplot analysis revealed the strains originated from complex genomic reassortment among different NBVs, suggesting the viruses experienced a high reassortment rate. In addition, strains successfully isolated from bat flies also implied that blood-sucking arthropods might serve as potential transmission vectors. IMPORTANCE Bats are the reservoir of many viral pathogens with strong pathogenicity, including NBVs. Nevertheless, it is unclear whether arthropod vectors are involved in transmitting NBVs. In this study, we successfully isolated two NBV strains from bat flies collected from the body surface of bats, which implies that they may be vectors for virus transmission between bats. While the potential threat to humans remains to be determined, evolutionary analyses involving different segments revealed that the novel strains had complex reassortment histories, with S1, S2, and M1 segments highly similar to human pathogens. Further experiments are required to determine whether more NBVs are vectored by bat flies, their potential threat to humans, and transmission dynamics.

Yunnan province in China shares its borders with three neighboring countries: Myanmar, Vietnam, and Laos. The region is characterized by a diverse climate and is known to be a suitable habitat for various arthropods, including midges which are notorious for transmitting diseases which pose significant health burdens affecting both human and animal health. A total of 431,100 midges were collected from 15 different locations in the border region of Yunnan province from 2015 to 2020. These midges were divided into 37 groups according to the collection year and sampling site. These 37 groups of midges were then homogenized to extract nucleic acid. Metatranscriptomics were used to analyze their viromes. Based on the obtained cytochrome C oxidase I gene (COI) sequences, three genera were identified, including one species of Forcipomyia, one species of Dasyhelea, and twenty-five species of Culicoides. We identified a total of 3199 viruses in five orders and 12 families, including 1305 single-stranded positive-stranded RNA viruses (+ssRNA) in two orders and seven families, 175 single-stranded negative-stranded RNA viruses (−ssRNA) in two orders and one family, and 1719 double-stranded RNA viruses in five families. Six arboviruses of economic importance were identified, namely Banna virus (BAV), Japanese encephalitis virus (JEV), Akabane virus (AKV), Bluetongue virus (BTV), Tibetan circovirus (TIBOV), and Epizootic hemorrhagic disease virus (EHDV), all of which are capable, to varying extents, of causing disease in humans and/or animals. The survey sites in this study basically covered the current distribution area of midges in Yunnan province, which helps to predict the geographic expansion of midge species. The complexity and diversity of the viral spectrum carried by midges identified in the study calls for more in-depth research, which can be utilized to monitor arthropod vectors and to predict the emergence and spread of zoonoses and animal epidemics, which is of great significance for the control of vector-borne diseases.

Between 7 December 2022 and 28 February 2023, China experienced a new wave of COVID-19 that swept across the entire country and resulted in an increasing amount of respiratory infections and hospitalizations. The purpose of this study is to reveal the intensity and composition of coinfecting microbial agents. In total, 196 inpatients were recruited from The Third People’s Hospital of Shenzhen, and 169 respiratory and 73 blood samples were collected for metagenomic next-generation sequencing. The total “Infectome” was characterized and compared across different groups defined by the SARS-CoV-2 detection status, age groups, and severity of disease. Our results revealed a total of 22 species of pathogenic microbes (4 viruses, 13 bacteria, and 5 fungi), and more were discovered in the respiratory tract than in blood. The diversity of the total infectome was highly distinguished between respiratory and blood samples, and it was generally higher in patients that were SARS-CoV-2-positive, older in age, and with more severe disease. At the individual pathogen level, HSV-1 seemed to be the major contributor to these differences observed in the overall comparisons. Collectively, this study reveals the highly complex respiratory infectome and high-intensity coinfection in patients admitted to the hospital during the period of the 2023 COVID-19 pandemic in China.

Abstract Background Ticks (order Ixodida) are ectoparasites, vectors and reservoirs of many infectious agents affecting humans and domestic animals. However, the lack of information on tick genomic diversity leaves significant gaps in the understanding of the evolution of ticks and associated bacteria. Results We collected > 20,000 contemporary and historical (up to 60 years of preservation) tick samples representing a wide range of tick biodiversity across diverse geographic regions in China. Metagenomic sequencing was performed on individual ticks to obtain the complete or near-complete mitochondrial (mt) genome sequences from 46 tick species, among which mitochondrial genomes of 23 species were recovered for the first time. These new mt genomes data greatly expanded the diversity of many tick groups and revealed five cryptic species. Utilizing the same metagenomic sequence data we identified divergent and abundant bacteria in Haemaphysalis, Ixodes, Dermacentor and Carios ticks, including nine species of pathogenetic bacteria and potentially new species within the genus Borrelia. We also used these data to explore the evolutionary relationship between ticks and their associated bacteria, revealing a pattern of long-term co-divergence relationship between ticks and Rickettsia and Coxiella bacteria. Conclusions In sum, our study provides important new information on the genetic diversity of ticks based on an analysis of mitochondrial DNA as well as on the prevalence of tick-borne pathogens in China. It also sheds new light on the long-term evolutionary and ecological relationships between ticks and their associated bacteria. Graphical Abstract

SARS-CoV Spike (S) protein shares considerable homology with SARS-CoV-2 S, especially in the conserved S2 subunit (S2). S protein mediates coronavirus receptor binding and membrane fusion, and the latter activity can greatly influence coronavirus infection. We observed that SARS-CoV S is less effective in inducing membrane fusion compared with SARS-CoV-2 S. We identify that S813T mutation is sufficient in S2 interfering with the cleavage of SARS-CoV-2 S by TMPRSS2, reducing spike fusogenicity and pseudoparticle entry. Conversely, the mutation of T813S in SARS-CoV S increased fusion ability and viral replication. Our data suggested that residue 813 in the S was critical for the proteolytic activation, and the change from threonine to serine at 813 position might be an evolutionary feature adopted by SARS-2-related viruses. This finding deepened the understanding of Spike fusogenicity and could provide a new perspective for exploring Sarbecovirus' evolution.

Abstract Background The global pork industry is continuously affected by infectious diseases that can result in large-scale mortality, trade restrictions, and major reductions in production. Nevertheless, the cause of many infectious diseases in pigs remains unclear, largely because commonly used diagnostic tools fail to capture the full diversity of potential pathogens and because pathogen co-infection is common. Results We used a meta-transcriptomic approach to systematically characterize the pathogens in 136 clinical cases representing different disease syndromes in pigs, as well as in 12 non-diseased controls. This enabled us to simultaneously determine the diversity, abundance, genomic information, and detailed epidemiological history of a wide range of potential pathogens. We identified 34 species of RNA viruses, nine species of DNA viruses, seven species of bacteria, and three species of fungi, including two novel divergent members of the genus Pneumocystis. While most of these pathogens were only apparent in diseased animals or were at higher abundance in diseased animals than in healthy animals, others were present in healthy controls, suggesting opportunistic infections. Importantly, most of the cases examined here were characterized by co-infection with more than two species of viral, bacterial, or fungal pathogens, some with highly correlated occurrence and abundance levels. Examination of clinical signs and necropsy results in the context of relevant pathogens revealed that a multiple-pathogen model was better associated with the data than a single-pathogen model was. Conclusions Our data demonstrate that most of the pig diseases examined were better explained by the presence of multiple rather than single pathogens and that infection with one pathogen can facilitate infection or increase the prevalence/abundance of another. Consequently, it is generally preferable to consider the cause of a disease based on a panel of co-infecting pathogens rather than on individual infectious agents. Video abstract

China is one of the countries with the richest wildlife population. The large variety of widely distributed species act as natural or susceptible hosts for numerous infectious diseases. It is estimated that there are more than 1.2 million unknown virus species in China, and there might be 10,000–30,000 unknown bacteria in wild mammals on the Qinghai-Tibet Plateau alone. There are no less than 600,000 species of animal-borne parasites and approximately 2 million species of fungi worldwide. With rapid economic growth and globalization, humans and wildlife interact more frequently, which enhances the probability of wildlife-borne pathogens infecting humans. The occurrence of animal-borne infectious diseases will become the “new normal” we have to face in the future. Therefore, research should be carried out on wildlife-borne microorganisms and the prevention and control of emerging infectious diseases to establish an analytical framework and an evaluation technology system for risk assessment and early warning of potential animal-borne emerging infectious diseases. This will not only improve our understanding of wildlife-borne microbial communities but also enable in-depth analysis, discovery, early warning, and even prediction of major animal-borne emerging infectious diseases that might occur in the future. Furthermore, this research will reduce response times, minimize the social and economic impact and losses, enable interventions related to the emergence or spread of the disease as early as possible, and comprehensively improve our management of infectious disease outbreaks.

Abstract RNA viruses (realm: Riboviria), including RNA phages and eukaryote‐infecting RNA viruses, are essential components of marine ecosystems. A large number of marine RNA viruses have been discovered in the last two decades because of the rapid development of next‐generation sequencing (NGS) technology. Indeed, the combination of NGS and state‐of‐the‐art meta‐omics methods (viromics, the study of all viruses in a specific environment) has led to a fundamental understanding of the taxonomy and genetic diversity of RNA viruses in the sea, suggesting the complex ecological roles played by RNA viruses in this complex ecosystem. Furthermore, comparisons of viromes in the context of highly variable marine niches reveal the biogeographic patterns and ecological impact of marine RNA viruses, whose role in global ecology is becoming increasingly clearer. In this review, we summarize the characteristics of the global marine RNA virosphere and outline the taxonomic hierarchy of RNA viruses with a specific focus on their ancient evolutionary history. We also review the development of methodology and the major progress resulting from its applications in RNA viromics. The aim of this review is not only to provide an in‐depth understanding of multifaceted aspects of marine RNA viruses, but to offer future perspectives on developing a better methodology for discovery, and exploring the evolutionary origin and major ecological significance of marine RNA virosphere.

SARS-CoV-2, the causative agent for COVID-19, infect human mainly via respiratory tract, which is heavily inhabited by local microbiota. However, the interaction between SARS-CoV-2 and nasopharyngeal microbiota, and the association with metabolome has not been well characterized. Here, metabolomic analysis of blood, urine, and nasopharyngeal swabs from a group of COVID-19 and non-COVID-19 patients, and metagenomic analysis of pharyngeal samples were used to identify the key features of COVID-19. Results showed lactic acid, l-proline, and chlorogenic acid methyl ester (CME) were significantly reduced in the sera of COVID-19 patients compared with non-COVID-19 ones. Nasopharyngeal commensal bacteria including Gemella morbillorum, Gemella haemolysans and Leptotrichia hofstadii were notably depleted in the pharynges of COVID-19 patients, while Prevotella histicola, Streptococcus sanguinis, and Veillonella dispar were relatively increased. The abundance of G. haemolysans and L. hofstadii were significantly positively associated with serum CME, which might be an anti-SARS-CoV-2 bacterial metabolite. This study provides important information to explore the linkage between nasopharyngeal microbiota and disease susceptibility. The findings were based on a very limited number of patients enrolled in this study; a larger size of cohort will be appreciated for further investigation.

ABSTRACT Porcine reproductive and respiratory syndrome (PRRS) is one of the most economically important diseases affecting the pig-raising industry. The PRRS virus (PRRSV) has high genetic diversity, partly owing to viral recombination. Some individual recombinant type 2 PRRSV (PRRSV-2) strains have been detected; however, the sequence composition characteristics of recombination hot spots and potential driving forces for recombinant PRRSV-2 are still unreported. Therefore, all available genomic sequences of PRRSV-2 (n = 949, including 29 genomes sequenced in this study) from 11 countries from 1991 to 2021 were collected and analyzed. The results revealed that the dominant major recombinant parent has been converted from lineage 3 (L3) to L1 since 2012. The recombination hot spots were located at nucleotides (nt) 7900 to 8200 (in NSP9, encoding viral RNA-dependent RNA polymerase) and nt 12500 to nt 13300 (in ORF2-ORF4, mean ORF2 to ORF4); no AU-rich characteristics were found in the recombination hot spots. Based on infectious clones of L1 and L8 PRRSV-2, recombinant PRRSVs were generated by switching complete or partial NSP9 (harboring the recombination hot spot). The results showed that recombinant PRRSVs based on the L1 backbone, but not the L8 backbone, acquired a higher replication capacity in pig primary alveolar macrophages. These findings will help to understand the reason behind the dominance of L1-based recombination in PRRSV-2 strains and provide new clues for an in-depth study of the recombination mechanism of PRRSV-2. IMPORTANCE Recombination is an important driver of the genetic shifts that are tightly linked to the evolution of RNA viruses. Viral recombination contributes substantially to the emergence of new variants, alterations in virulence, and pathogenesis. PRRSV is genetically diverse, partly because of extensive recombination. In this study, we analyzed interlineage recombination based on available genomic sequences of PRRSV-2 from 1991 to 2021. The study revealed the temporal and geographical distribution of recombinant PRRSVs and the recombination hot spot’s location and showed that artificially constructed recombinant PRRSVs (harboring a high-frequency region) had more viral genomic copies than their parental virus, indicating that dominant recombination was shaped by a tendency to benefit viral replication. This finding will enrich our understanding of PRRSV recombination and provide new clues for an in-depth study of the recombination mechanism.

African swine fever (ASF) outbreak have caused tremendous economic loss to the pig industry in China since its emergence in August 2018. Previous studies revealed that many published sequences are not suitable for detailed analyses due to the lack of data regarding quality parameters and methodology, and outdated annotations. Thus, high-quality genomes of highly pathogenic strains that can be used as references for early Chinese ASF outbreaks are still lacking, and little is known about the features of intra-host variants of ASF virus (ASFV). In this study, a full genome sequencing of clinical samples from the first ASF outbreak in Guangdong in 2018 was performed using MGI (MGI Tech Co., Ltd., Shenzhen, China) and Nanopore sequencing platforms, followed by Sanger sequencing to verify the variations. With 22 sequencing corrections, we obtained a high-quality genome of one of the earliest virulent isolates, GZ201801_2. After proofreading, we improved (add or modify) the annotations of this isolate using the whole genome alignment with Georgia 2007/1. Based on the complete genome sequence, we constructed the methylation profiles of early ASFV strains in China and predicted the potential 5mC and 6mA methylation sites, which are likely involved in metabolism, transcription, and replication. Additionally, the intra-host single nucleotide variant distribution and mutant allele frequency in the clinical samples of early strain were determined for the first time and found a strong preference for A and T substitution mutation, non-synonymous mutations, and mutations that resulted in amino acid substitutions into Lysine. In conclusion, this study provides a high-quality genome sequence, updated genome annotation, methylation profile, and mutation spectrum of early ASFV strains in China, thereby providing a reference basis for further studies on the evolution, transmission, and virulence of ASFV.

Abstract B cell response plays a critical role against SARS-CoV-2 infection. However, little is known about the diversity and frequency of the paired SARS-CoV-2 antigen-specific BCR repertoire after SARS-CoV-2 infection. Here, we performed single-cell RNA sequencing and VDJ sequencing using the memory and plasma B cells isolated from five convalescent COVID-19 patients, and analyzed the spectrum and transcriptional heterogeneity of antibody immune responses. Via linking BCR to antigen specificity through sequencing (LIBRA-seq), we identified a distinct activated memory B cell subgroup (CD11c high CD95 high ) had a higher proportion of SARS-CoV-2 antigen-labeled cells compared with memory B cells. Our results revealed the diversity of paired BCR repertoire and the non-stochastic pairing of SARS-CoV-2 antigen-specific immunoglobulin heavy and light chains after SARS-CoV-2 infection. The public antibody clonotypes were shared by distinct convalescent individuals. Moreover, several antibodies isolated by LIBRA-seq showed high binding affinity against SARS-CoV-2 receptor-binding domain (RBD) or nucleoprotein (NP) via ELISA assay. Two RBD-reactive antibodies C14646P3S and C2767P3S isolated by LIBRA-seq exhibited high neutralizing activities against both pseudotyped and authentic SARS-CoV-2 viruses in vitro. Our study provides fundamental insights into B cell response following SARS-CoV-2 infection at the single-cell level.

Cats harbor many important viral pathogens, and the knowledge of their diversity has been greatly expanded thanks to increasingly popular molecular sequencing techniques. While the diversity is mostly described in numerous regionally defined studies, there lacks a global overview of the diversity for the majority of cat viruses, and therefore our understanding of the evolution and epidemiology of these viruses was generally inadequate. In this study, we analyzed 12,377 genetic sequences from 25 cat virus species and conducted comprehensive phylodynamic analyses. It revealed, for the first time, the global diversity for all cat viruses known to date, taking into account highly virulent strains and vaccine strains. From there, we further characterized and compared the geographic expansion patterns, temporal dynamics and recombination frequencies of these viruses. While respiratory pathogens such as feline calicivirus showed some degree of geographical panmixes, the other viral species are more geographically defined. Furthermore, recombination rates were much higher in feline parvovirus, feline coronavirus, feline calicivirus and feline foamy virus than the other feline virus species. Collectively, our findings deepen the understanding of the evolutionary and epidemiological features of cat viruses, which in turn provide important insight into the prevention and control of cat pathogens.

Abstract Background Prior studies reported that 5 ~ 32% COVID-19 patients were critically ill, a situation that poses great challenge for the management of the patients and ICU resources. We aim to identify independent risk factors to serve as prediction markers for critical illness of SARS-CoV-2 infection. Methods Fifty-two critical and 200 non-critical SARS-CoV-2 nucleic acid positive patients hospitalized in 15 hospitals outside Wuhan from January 19 to March 6, 2020 were enrolled in this study. Multivariable logistic regression and LASSO logistic regression were performed to identify independent risk factors for critical illness. Results Age older than 60 years, dyspnea, respiratory rate > 24 breaths per min, leukocytosis > 9.5 × 109/L, neutrophilia > 6.3 × 109/L, lymphopenia 3.53, fibrinogen > 4 g/L, d-dimer > 0.55 μg/mL, blood urea nitrogen > 7.1 mM, elevated aspartate transaminase, elevated alanine aminotransferase, total bilirubin > 21 μM, and Sequential Organ Failure Assessment (SOFA) score ≥ 2 were identified as risk factors for critical illness. LASSO logistic regression identified the best combination of risk factors as SOFA score, age, dyspnea, and leukocytosis. The Area Under the Receiver-Operator Curve values for the risk factors in predicting critical illness were 0.921 for SOFA score, 0.776 for age, 0.764 for dyspnea, 0.658 for leukocytosis, and 0.960 for the combination of the four risk factors. Conclusions Our findings advocate the use of risk factors SOFA score ≥ 2, age > 60, dyspnea and leukocytosis > 9.5 × 109/L on admission, alone or in combination, to determine the optimal management of the patients and health care resources.

We isolated 17 viral strains capable of causing cytopathic effects in mammalian cells and death in neonatal mice from sand flies in China. Phylogenetic analysis showed that these strains belonged to the genus Phlebovirus. These findings highlight the need to control this potentially emerging virus to help safeguard public health.

At the end of 2019 Wuhan witnessed an outbreak of “atypical pneumonia” that later developed into a global pandemic. Metagenomic sequencing rapidly revealed the causative agent of this outbreak to be a novel coronavirus denoted SARS-CoV-2. To provide a snapshot of the pathogens in pneumonia-associated respiratory samples from Wuhan prior to the emergence of SARS-CoV-2, we collected bronchoalveolar lavage fluid samples from 408 patients presenting with pneumonia and acute respiratory infections at the Central Hospital of Wuhan between 2016 and 2017. Unbiased total RNA sequencing was performed to reveal their “total infectome”, including viruses, bacteria and fungi. We identified 35 pathogen species, comprising 13 RNA viruses, 3 DNA viruses, 16 bacteria and 3 fungi, often at high abundance and including multiple co-infections (13.5%). SARS-CoV-2 was not present. These data depict a stable core infectome comprising common respiratory pathogens such as rhinoviruses and influenza viruses, an atypical respiratory virus (EV-D68), and a single case of a sporadic zoonotic pathogen–Chlamydia psittaci. Samples from patients experiencing respiratory disease on average had higher pathogen abundance than healthy controls. Phylogenetic analyses of individual pathogens revealed multiple origins and global transmission histories, highlighting the connectedness of the Wuhan population. This study provides a comprehensive overview of the pathogens associated with acute respiratory infections and pneumonia, which were more diverse and complex than obtained using targeted PCR or qPCR approaches. These data also suggest that SARS-CoV-2 or closely related viruses were absent from Wuhan in 2016–2017. Author summary Acute respiratory infections and pneumonia are a significant public health concern on a global scale. This was brought sharply into focus by the emergence of COVID-19 at the end of 2019 in Wuhan, China. Far less is known, however, about the total “infectomes” associated with respiratory infection and pneumonia: that is, all the viruses, bacteria and fungi that might be associated with these important human diseases. Herein, we describe total infectomes in bronchoalveolar lavage fluid samples taken from 408 patients presenting with pneumonia and acute respiratory infection in a single hospital in Wuhan between 2016 and 2017. From these samples we identified 35 pathogen species, comprising 13 RNA viruses, 3 DNA viruses, 16 bacteria and 3 fungi, often at high abundance and including multiple co-infections. We found no evidence for SARS-CoV-2. These data reveal a stable core infectome comprising common respiratory pathogens as well as their history of global transmission.

Rhesus macaques (Macaca mulatta) are the most widely distributed species of Old World monkey and are frequently used as animal models to study human health and disease. Their gastrointestinal microbial community likely plays a major role in their physiology, ecology and evolution. Herein, we compared the fecal microbiome and antibiotic resistance genes in 15 free-ranging and 81 zoo-captive rhesus macaques sampled from two zoos in China, using both 16S amplicon sequencing and whole genome shotgun DNA sequencing approaches. Our data revealed similar levels of microbial diversity/richness among the three groups, although the composition of each group differed significantly and were particularly marked between the two zoo-captive and one wild groups. Zoo-captive animals also demonstrated a greater abundance and diversity of antibiotic genes. Through whole genome shotgun sequencing we also identified a mammalian (simian) associated adenovirus. Overall, this study provides a comprehensive analysis of resistomes and microbiomes in zoo-captive and free-ranging monkeys, revealing that semi-captive wildlife might harbor a higher diversity of antimicrobial resistant genes.

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 Wēnzhōu 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.

ABSTRACTFrom December 2019, an outbreak of unusual pneumonia was reported in Wuhan with many cases linked to Huanan Seafood Market that sells seafood as well as live exotic animals. We investigated two patients who developed acute respiratory syndromes after independent contact history with this market. The two patients shared common clinical features including fever, cough, and multiple ground-glass opacities in the bilateral lung field with patchy infiltration. Here, we highlight the use of a low-input metagenomic next-generation sequencing (mNGS) approach on RNA extracted from bronchoalveolar lavage fluid (BALF). It rapidly identified a novel coronavirus (named 2019-nCoV according to World Health Organization announcement) which was the sole pathogens in the sample with very high abundance level (1.5% and 0.62% of total RNA sequenced). The entire viral genome is 29,881 nt in length (GenBank MN988668 and MN988669, Sequence Read Archive database Bioproject accession PRJNA601736) and is classified into β-coronavirus genus. Phylogenetic analysis indicates that 2019-nCoV is close to coronaviruses (CoVs) circulating in Rhinolophus (Horseshoe bats), such as 98.7% nucleotide identity to partial RdRp gene of bat coronavirus strain BtCoV/4991 (GenBank KP876546, 370 nt sequence of RdRp and lack of other genome sequence) and 87.9% nucleotide identity to bat coronavirus strain bat-SL-CoVZC45 and bat-SL-CoVZXC21. Evolutionary analysis based on ORF1a/1b, S, and N genes also suggests 2019-nCoV is more likely a novel CoV independently introduced from animals to humans.

ABSTRACTCirculating in China and 158 other countries and areas, the ongoing COVID-19 outbreak has caused devastating mortality and posed a great threat to public health. However, efforts to identify effectively supportive therapeutic drugs and treatments has been hampered by our limited understanding of host immune response for this fatal disease. To characterize the transcriptional signatures of host inflammatory response to SARS-CoV-2 (HCoV-19) infection, we carried out transcriptome sequencing of the RNAs isolated from the bronchoalveolar lavage fluid (BALF) and peripheral blood mononuclear cells (PBMC) specimens of COVID-19 patients. Our results reveal distinct host inflammatory cytokine profiles to SARS-CoV-2 infection in patients, and highlight the association between COVID-19 pathogenesis and excessive cytokine release such as CCL2/MCP-1, CXCL10/IP-10, CCL3/MIP-1A, and CCL4/MIP1B. Furthermore, SARS-CoV-2 induced activation of apoptosis and P53 signalling pathway in lymphocytes may be the cause of patients’ lymphopenia. The transcriptome dataset of COVID-19 patients would be a valuable resource for clinical guidance on anti-inflammatory medication and understanding the molecular mechansims of host response.

Abstract Background Australian brushtail possums (Trichosurus vulpecula) are an introduced pest species in New Zealand, but native to Australia where they are protected for biodiversity conservation. Wobbly possum disease (WPD) is a fatal neurological disease of Australian brushtail possums described in New Zealand populations that has been associated with infection by the arterivirus (Arteriviridae) wobbly possum disease virus (WPDV-NZ). Clinically, WPD-infected possums present with chronic meningoencephalitis, choroiditis and multifocal neurological symptoms including ataxia, incoordination, and abnormal gait. Methods We conducted a retrospective investigation to characterise WPD in native Australian brushtail possums, and used a bulk meta-transcriptomic approach (i.e. total RNA-sequencing) to investigate its potential viral aetiology. PCR assays were developed for case diagnosis and full genome recovery in the face of extensive genetic variation. Results We identified genetically distinct lineages of arteriviruses from archival tissues of WPD-infected possums in Australia, termed wobbly possum disease virus AU1 and AU2. Phylogenetically, WPDV-AU1 and WPDV-AU2 shared only ~ 70% nucleotide similarity to each other and the WPDV-NZ strain, suggestive of a relatively ancient divergence. Notably, we also identified a novel and divergent hepacivirus (Flaviviridae) - the first in a marsupial - in both WPD-infected and uninfected possums, indicative of virus co-infection. Conclusions We have identified marsupial-specific lineages of arteriviruses in mainland Australia that are genetically distinct from that in New Zealand, in some cases co-infecting animals with a novel hepacivirus. Our study provides new insight into the hidden genetic diversity of arteriviruses, the capacity for virus co-infection, and highlights the utility of meta-transcriptomics for disease investigation in a One Health context.

Abstract Background Antibiotic resistance is rendering common bacterial infections untreatable. Wildlife can incorporate and disperse antibiotic-resistant bacteria in the environment, such as water systems, which in turn serve as reservoirs of resistance genes for human pathogens. Anthropogenic activity may contribute to the spread of bacterial resistance cycling through natural environments, including through the release of human waste, as sewage treatment only partially removes antibiotic-resistant bacteria. However, empirical data supporting these effects are currently limited. Here we used bulk RNA-sequencing (meta-transcriptomics) to assess the diversity and expression levels of functionally viable resistance genes in the gut microbiome of birds with aquatic habits in diverse locations. Results We found antibiotic resistance genes in birds from all localities, from penguins in Antarctica to ducks in a wastewater treatment plant in Australia. Comparative analysis revealed that birds feeding at the wastewater treatment plant carried the greatest resistance gene burden, including genes typically associated with multidrug resistance plasmids as the aac(6)-Ib-cr gene. Differences in resistance gene burden also reflected aspects of bird ecology, taxonomy, and microbial function. Notably, ducks, which feed by dabbling, carried a higher abundance and diversity of resistance genes than turnstones, avocets, and penguins, which usually prey on more pristine waters. Conclusions These transcriptome data suggest that human waste, even if it undergoes treatment, might contribute to the spread of antibiotic resistance genes to the wild. Differences in microbiome functioning across different bird lineages may also play a role in the antibiotic resistance burden carried by wild birds. In summary, we reveal the complex factors explaining the distribution of resistance genes and their exchange routes between humans and wildlife, and show that meta-transcriptomics is a valuable tool to access functional resistance genes in whole microbial communities.

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

Feline panleukopenia (FPL), a highly contagious and frequently fatal disease of cats, is caused by Feline parvovirus (FPV) and Canine parvovirus (CPV). We characterised the diversity of these Carnivore protoparvovirus 1 variants in 18 faecal samples collected from domestic cats with FPL during an outbreak, using targeted parvoviral DNA metagenomics to a mean depth of >10,000 × coverage per site. All samples comprised FPV alone. Compared with the reference FPV genome, isolated in 1967, 44 mutations were detected. Ten of these were nonsynonymous, including 9 in nonstructural genes and one in VP1/VP2 (Val232Ile), which was the only one to exhibit interhost diversity, being present in five sequences. There were five other polymorphic nucleotide positions, all with synonymous mutations. Intrahost diversity at all polymorphic positions was low, with subconsensus variant frequencies (SVF) of

The family Rhabdoviridae contain important human and mammalian pathogens that are vectored by different arthropod species. The ground supernatants of mosquitoes were used to inoculate in BHK-21 and C6/36 cells for virus isolation. Then, the viral complete genome sequence was obtained and used for phylogenetic analysis. In this study, we observed a cytopathic effect (CPE) in mosquito cells (C6/36) and rod-like virion after inoculating a pool of Armigeres subalbatus samples collected in Shanxi Province, China, in 2019 (SX1916). Meta-transcriptomics sequencing revealed the presence of two distinctive rhabdoviruses with similar abundance levels, namely, Shanxi Armigeres subalbatus rhabdovirus (SXARV) and Shanxi Arboretum virus (SXABTV). Despite the fact that the SXARV genome (9590 nt) was much shorter than that of SXABTV (11,480 nt), both belonged to the Almendravirus group within Rhabdoviridae whose genomes encoded five proteins (N, P, M, G, and L) and a small hydrophobin (U1) and the difference in lengths is mainly caused by a substantially shorter N protein encoded by SXARV. On the phylogenetic tree, SXABTV was closely related (90.7% amino acid identity at L protein) with the Arboretum virus isolated from Psorophora albigenu mosquitoes in Peru in 2014, whereas SXARV was distantly related to Rio Chico virus (63.3% amino acid identity), a genetic distance large enough to be defined as a new species within Rhabdoviridae. Collectively, we report a simultaneous isolation of two related rhabdoviruses from Armigeres subalbatus that marked the circulation of almendraviruses in Shanxi, China.

ABSTRACT Increasing data indicate that insects serve as major reservoirs and vectors of viruses, which account for the continuously increasing ecological burden and infectious disease outbreaks. Uncovering the hidden diversity of viruses in insects will further the understanding of the ecological and evolutionary perspectives in the emergence of insect-associated virus diseases. In this study, we queried transcriptome sequencing (RNA-Seq) data from more than 600 species across 32 insect orders dwelling in different ecological habitats and recovered more than 1,213 RNA viruses that were recapitulated in 40 families, 2 unclassified genera, and many unspecified viral groups. These novel viruses included the well-known insect-associated viruses within Flaviviridae, Picornavirales, Bunyavirales, Mononegavirales, Nidovirales, Reoviridae, and Negevirus. More appeared to form novel clusters within previously described taxa or could be resolved as paraphyletic, including the first astrovirus identified in insects, in which many were sufficiently divergent to warrant the establishment of new virus genera or families. Additionally, some viruses were closely related to the recognized plant-, fungus-, and vertebrate-specific species, implying the importance of relationships between insect behavior and virus spread. Comparative genome analyses also revealed high genomic variability with respect to the flexible gene pool and genome architecture of these newly described viruses, including the evidence for genome reshuffling first discovered in Dicistroviridae. The data reflecting the genetically and ecologically diverse viral populations in insects greatly expand our understanding of RNA viruses in nature and highlight that the biodiversity of RNA viruses remains largely unexplored. IMPORTANCE Insects comprise the largest proportion of animals on earth and are frequently implicated in the transmission of vector-borne diseases. However, considerable attention has been paid to the phytophagous and hematophagous insects, with results that provide insufficient and biased information about the viruses in insects. Here, we have delivered compelling evidence for the exceptional abundance and genetic diversity of RNA viruses in a wide range of insects. Novel viruses were found to cover major categories of RNA viruses, and many formed novel clusters divergent from the previously described taxa, dramatically broadening the range of known RNA viruses in insects. These newly characterized RNA viruses exhibited high levels of genomic plasticity in genome size, open reading frame (ORF) number, intergenic structure, and gene rearrangement and segmentation. This work provides comprehensive insight into the origin, spread, and evolution of RNA viruses. Of course, a large-scale virome project involving more organisms would provide more-detailed information about the virus infections in insects.

We isolated Tamdy virus (TAMV; strain XJ01/TAMV/China/2018) from Hyalomma asiaticum ticks infesting Bactrian camels in Xinjiang, China, in 2018. The genome of the strain showed high nucleotide similarity with previously described TAMV strains from Asia. Our study highlights the potential threat of TAMV to public health in China.

Dear Professor Remi N. Charrel and Professor Jerome Depaquit, we thank you for your interest in our research and for your kind suggestions [...]

Mammal-associated coronaviruses have a long evolutionary history across global bat populations, which makes them prone to be the most likely ancestral origins of coronavirus-associated epidemics and pandemics globally. Limited coronavirus research has occurred at the junction of Europe and Asia, thereby investigations in Georgia are critical to complete the coronavirus diversity map in the region. We conducted a cross-sectional coronavirus survey in bat populations at eight locations of Georgia, from July to October of 2014. We tested 188 anal swab samples, remains of previous pathogen discovery studies, for the presence of coronaviruses using end-point pan-coronavirus RT-PCR assays. Samples positive for a 440 bp amplicon were Sanger sequenced to infer coronavirus subgenus or species through phylogenetic reconstructions. Overall, we found a 24.5% positive rate, with 10.1% for Alphacoronavirus and 14.4% for Betacoronavirus. Albeit R. euryale, R. ferrumequinum, M. blythii and M. emarginatus were found infected with both CoV genera, we could not rule out CoV co-infection due to limitation of the sequencing method used and sample availability. Based on phylogenetic inferences and genetic distances at nucleotide and amino acid levels, we found one putative new subgenus and three new species of Alphacoronavirus, and two new species of Betacoronavirus.

Feline calicivirus (FCV) causes upper respiratory tract disease (URTD) and sporadic outbreaks of virulent systemic disease (FCV-VSD). The basis for the increased pathogenicity of FCV-VSD viruses is incompletely understood, and antivirals for FCV-VSD have yet to be developed. We investigated the clinicoepidemiology and viral features of three FCV-VSD outbreaks in Australia and evaluated the in vitro efficacy of nitazoxanide (NTZ), 2′-C-methylcytidine (2CMC) and NITD-008 against FCV-VSD viruses. Overall mortality among 23 cases of FCV-VSD was 39%. Metagenomic sequencing identified five genetically distinct FCV lineages within the three outbreaks, all seemingly evolving in situ in Australia. Notably, no mutations that clearly distinguished FCV-URTD from FCV-VSD phenotypes were identified. One FCV-URTD strain likely originated from a recombination event. Analysis of seven amino-acid residues from the hypervariable E region of the capsid in the cultured viruses did not support the contention that properties of these residues can reliably differentiate between the two pathotypes. On plaque reduction assays, dose–response inhibition of FCV-VSD was obtained with all antivirals at low micromolar concentrations; NTZ EC50, 0.4–0.6 µM, TI = 21; 2CMC EC50, 2.7–5.3 µM, TI > 18; NITD-008, 0.5 to 0.9 µM, TI > 111. Investigation of these antivirals for the treatment of FCV-VSD is warranted.

The occurrence of Francisella tularensis outside of endemic areas, such as North America and Eurasia, has been enigmatic. We report the metagenomic discovery and isolation of F. tularensis ssp. holarctica biovar japonica from diseased ringtail possums in Sydney, Australia. This finding confirms the presence of F. tularensis in the Southern Hemisphere.

In a previous study, a metatranscriptomics survey of RNA viruses in several important lower vertebrate host groups revealed huge viral diversity, transforming the understanding of the evolution of vertebrate-associated RNA virus groups. However, the diversity of the DNA and retro-transcribing viruses in these host groups was left uncharacterized. Given that RNA sequencing is capable of revealing viruses undergoing active transcription and replication, we collected previously generated datasets associated with lower vertebrate hosts, and searched them for DNA and retro-transcribing viruses. Our results revealed the complete genome, or “core gene sets”, of 18 vertebrate-associated DNA and retro-transcribing viruses in cartilaginous fishes, ray-finned fishes, and amphibians, many of which had high abundance levels, and some of which showed systemic infections in multiple organs, suggesting active transcription or acute infection within the host. Furthermore, these new findings recharacterized the evolutionary history in the families Hepadnaviridae, Papillomaviridae, and Alloherpesviridae, confirming long-term virus–host codivergence relationships for these virus groups. Collectively, our results revealed reliable and sufficient information within metatranscriptomics sequencing to characterize not only RNA viruses, but also DNA and retro-transcribing viruses, and therefore established a key methodology that will help us to understand the composition and evolution of the total “infectome” within a diverse range of vertebrate hosts.

We report the isolation of a newly recognized phlebovirus, Hedi virus (HEDV), from Phlebotomus chinensis sandflies collected in Shanxi Province, China. The virus’ RNA is comprised of three segments. The greatest amino acid sequence similarity of the three gene segments between this virus and previously recognized phleboviruses is 40.85–63.52%, and the RNA-dependent RNA polymerase (RdRp) amino acid sequence has the greatest similarity (63.52%) to the Rift Valley fever virus (RVFV) ZH-548 strain. Phylogenetic analysis of the amino acid sequence of the virus RdRp indicated that HEDV is close to RVFV and distinct from other phleboviruses, forming its own evolutionary branch. We conclude that it is necessary to increase the monitoring of phleboviruses carried by sandflies in China.

Middle East respiratory syndrome coronavirus (MERS-CoV) was identified on the Arabian Peninsula in 2012 and is still causing cases and outbreaks in the Middle East. When MERS-CoV was first identified, the closest related virus was in bats; however, it has since been recognized that dromedary camels serve as a virus reservoir and potential source for human infections. A total of 376 camels were screened for MERS-Cov at a live animal market in the Eastern Region of the Emirate of Abu Dhabi, UAE. In all, 109 MERS-CoV-positive camels were detected in week 1, and a subset of positive camels were sampled again weeks 3 through 6. A total of 126 full and 3 nearly full genomes were obtained from 139 samples. Spike gene sequences were obtained from 5 of the 10 remaining samples. The camel MERS-CoV genomes from this study represent 3 known and 2 potentially new lineages within clade B. Within lineages, diversity of camel and human MERS-CoV sequences are intermixed. We identified sequences from market camels nearly identical to the previously reported 2015 German case who visited the market during his incubation period. We described 10 recombination events in the camel samples. The most frequent recombination breakpoint was the junctions between ORF1b and S. Evidence suggests MERS-CoV infection in humans results from continued introductions of distinct MERS-CoV lineages from camels. This hypothesis is supported by the camel MERS-CoV genomes sequenced in this study. Our study expands the known repertoire of camel MERS-CoVs circulating on the Arabian Peninsula.Emerging Microbes & Infections (2017) 6, e101; doi:10.1038/emi.2017.89; published online 8 November 2017

Tilapia lake virus (TiLV) has caused mass mortalities in farmed and wild tilapia with serious economic and ecological consequences. Until recently, this virus was the sole member of the Amnoonviridae, a family within the order Articulavirales comprising segmented negative-sense RNA viruses. We sought to identify additional viruses within the Amnoonviridae through total RNA sequencing (meta-transcriptomics) and data mining of published transcriptomes. Accordingly, we sampled marine fish species from both Australia and China and discovered several segments of two new viruses within the Amnoonviridae, tentatively called Flavolineata virus and Piscibus virus, respectively. In addition, by mining vertebrate transcriptome data, we identified nine additional virus transcripts matching to multiple genomic segments of TiLV in both marine and freshwater fish. These new viruses retained sequence conservation with the distantly related Orthomyxoviridae in the RdRp subunit PB1, but formed a distinct and diverse phylogenetic group. These data suggest that the Amnoonviridae have a broad host range within fish and that greater animal sampling will identify additional divergent members of the Articulavirales.

Astroviruses, isolated from numerous avian and mammalian species including humans, are commonly associated with enteritis and encephalitis. Two astroviruses have previously been identified in cats, and while definitive evidence is lacking, an association with enteritis is suggested. Using metagenomic next-generation sequencing of viral nucleic acids from faecal samples, we identified two novel feline astroviruses termed Feline astrovirus 3 and 4. These viruses were isolated from healthy shelter-housed kittens (Feline astrovirus 3; 6448 bp) and from a kitten with diarrhoea that was co-infected with Feline parvovirus (Feline astrovirus 4, 6549 bp). Both novel astroviruses shared a genome arrangement of three open reading frames (ORFs) comparable to that of other astroviruses. Phylogenetic analysis of the concatenated ORFs, ORF1a, ORF1b and capsid protein revealed that both viruses were phylogenetically distinct from other feline astroviruses, although their precise evolutionary history could not be accurately determined due to a lack of resolution at key nodes. Large-scale molecular surveillance studies of healthy and diseased cats are needed to determine the pathogenicity of feline astroviruses as single virus infections or in co-infections with other enteric viruses.