Search

Malayan pangolin SARS-CoV-2-related coronavirus (SARSr-CoV-2) is closely related to SARS-CoV-2. However, little is known about its pathogenicity in pangolins. Using CT scans we show that SARSr-CoV-2 positive Malayan pangolins are characterized by bilateral ground-glass opacities in lungs in a similar manner to COVID-19 patients. Histological examination and blood gas tests are indicative of dyspnea. SARSr-CoV-2 infected multiple organs in pangolins, with the lungs the major target, and histological expression data revealed that ACE2 and TMPRSS2 were co-expressed with viral RNA. Transcriptome analysis indicated that virus-positive pangolins were likely to have inadequate interferon responses, with relative greater cytokine and chemokine activity in the lung and spleen. Notably, both viral RNA and viral proteins were detected in three pangolin fetuses, providing initial evidence for vertical virus transmission. In sum, our study outlines the biological framework of SARSr-CoV-2 in pangolins, revealing striking similarities to COVID-19 in humans.

A group convened and led by the Virus Evolution Working Group of the World Health Organization reports on its deliberations and announces a naming scheme that will enable clear communication about SARS-CoV-2 variants of interest and concern.

The ongoing pandemic spread of a new human coronavirus, SARS-CoV-2, which is associated with severe pneumonia/disease (COVID-19), has resulted in the generation of tens of thousands of virus genome sequences. The rate of genome generation is unprecedented, yet there is currently no coherent nor accepted scheme for naming the expanding phylogenetic diversity of SARS-CoV-2. Here, we present a rational and dynamic virus nomenclature that uses a phylogenetic framework to identify those lineages that contribute most to active spread. Our system is made tractable by constraining the number and depth of hierarchical lineage labels and by flagging and delabelling virus lineages that become unobserved and hence are probably inactive. By focusing on active virus lineages and those spreading to new locations, this nomenclature will assist in tracking and understanding the patterns and determinants of the global spread of SARS-CoV-2.

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.

The high numbers of COVID-19 cases and deaths in Brazil have made Latin America an epicentre of the pandemic. SARS-CoV-2 established sustained transmission in Brazil early in the pandemic, but important gaps remain in our understanding of virus transmission dynamics at a national scale. We use 17,135 near-complete genomes sampled from 27 Brazilian states and bordering country Paraguay. From March to November 2020, we detected co-circulation of multiple viral lineages that were linked to multiple importations (predominantly from Europe). After November 2020, we detected large, local transmission clusters within the country. In the absence of effective restriction measures, the epidemic progressed, and in January 2021 there was emergence and onward spread, both within and abroad, of variants of concern and variants under monitoring, including Gamma (P.1) and Zeta (P.2). We also characterized a genomic overview of the epidemic in Paraguay and detected evidence of importation of SARS-CoV-2 ancestor lineages and variants of concern from Brazil. Our findings show that genomic surveillance in Brazil enabled assessment of the real-time spread of emerging SARS-CoV-2 variants.

Most of our understanding of the ecology and evolution of avian influenza A virus (AIV) in wild birds is derived from studies conducted in the northern hemisphere on waterfowl, with a substantial bias towards dabbling ducks. However, relevant environmental conditions and patterns of avian migration and reproduction are substantially different in the southern hemisphere. Through the sequencing and analysis of 333 unique AIV genomes collected from wild birds collected over 15 years we show that Australia is a global sink for AIV diversity and not integrally linked with the Eurasian gene pool. Rather, AIV are infrequently introduced to Australia, followed by decades of isolated circulation and eventual extinction. The number of co-circulating viral lineages varies per subtype. AIV haemagglutinin (HA) subtypes that are rarely identified at duck-centric study sites (H8-12) had more detected introductions and contemporary co-circulating lineages in Australia. Combined with a lack of duck migration beyond the Australian-Papuan region, these findings suggest introductions by long-distance migratory shorebirds. In addition, we found no evidence of directional or consistent patterns in virus movement across the Australian continent. This feature corresponds to patterns of bird movement, whereby waterfowl have nomadic and erratic rainfall-dependant distributions rather than consistent intra-continental migratory routes. Finally, we detected high levels of virus gene segment reassortment, with a high diversity of AIV genome constellations across years and locations. These data, in addition to those from other studies in Africa and South America, clearly show that patterns of AIV dynamics in the Southern Hemisphere are distinct from those in the temperate north.Author SummaryA result of the ever-growing poultry industry is a dramatic global increase in the incidence of high pathogenicity avian influenza virus outbreaks. In contrast, wild birds are believed to be the main reservoir for low pathogenic avian influenza A virus. Due to intensive research and surveillance of AIV in waterfowl in the Northern Hemisphere, we have a better understanding of AIV ecology and evolution in that region compared to the Southern Hemisphere, which are characterised by different patterns of avian migration and ecological conditions. We analysed 333 unique AIV genomes collected from wild birds in Australia to understand how Australia fits into global AIV dynamics and how viruses are maintained and dispersed within the continent of Australia. We show that the Southern Hemisphere experiences differing evolutionary dynamics to those seen in Northern Hemisphere with Australia representing a global sink for AIV.

The diversity of lagoviruses (Caliciviridae) in Australia has increased considerably in recent years. By the end of 2017, five variants from three viral genotypes were present in populations of Australian rabbits, while prior to 2014 only two variants were known. To understand the evolutionary interactions among these lagovirus variants, we monitored their geographical distribution and relative incidence over time in a continental-scale competition study. Within 3 years of the incursion of rabbit haemorrhagic disease virus 2 (RHDV2, denoted genotype GI.1bP-GI.2 [polymerase genotype]P-[capsid genotype]) into Australia, two novel recombinant lagovirus variants emerged: RHDV2-4e (genotype GI.4eP-GI.2) in New South Wales and RHDV2-4c (genotype GI.4cP-GI.2) in Victoria. Although both novel recombinants contain non-structural genes related to those from benign, rabbit-specific, enterotropic viruses, these variants were recovered from the livers of both rabbits and hares that had died acutely. This suggests that the determinants of host and tissue tropism for lagoviruses are associated with the structural genes, and that tropism is intricately connected with pathogenicity. Phylogenetic analyses demonstrated that the RHDV2-4c recombinant emerged independently on multiple occasions, with five distinct lineages observed. Both the new RHDV2-4e and -4c recombinant variants replaced the previous dominant parental RHDV2 (genotype GI.1bP-GI.2) in their respective geographical areas, despite sharing an identical or near-identical (i.e. single amino acid change) VP60 major capsid protein with the parental virus. This suggests that the observed replacement by these recombinants was not driven by antigenic variation in VP60, implicating the non-structural genes as key drivers of epidemiological fitness. Molecular clock estimates place the RHDV2-4e recombination event in early to mid-2015, while the five RHDV2-4c recombination events occurred from late 2015 through to early 2017. The emergence of at least six viable recombinant variants within a 2-year period highlights the high frequency of these events, detectable only through intensive surveillance, and demonstrates the importance of recombination in lagovirus evolution.

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.

Geoghegan and Holmes describe the history of evolutionary ideas in the study of viruses, showing that two different approaches to studying virus evolution—the comparative and the experimental—were both established in seminal papers published in the late... RNA viruses are diverse, abundant, and rapidly evolving. Genetic data have been generated from virus populations since the late 1970s and used to understand their evolution, emergence, and spread, culminating in the generation and analysis of many thousands of viral genome sequences. Despite this wealth of data, evolutionary genetics has played a surprisingly small role in our understanding of virus evolution. Instead, studies of RNA virus evolution have been dominated by two very different perspectives, the experimental and the comparative, that have largely been conducted independently and sometimes antagonistically. Here, we review the insights that these two approaches have provided over the last 40 years. We show that experimental approaches using in vitro and in vivo laboratory models are largely focused on short-term intrahost evolutionary mechanisms, and may not always be relevant to natural systems. In contrast, the comparative approach relies on the phylogenetic analysis of natural virus populations, usually considering data collected over multiple cycles of virus–host transmission, but is divorced from the causative evolutionary processes. To truly understand RNA virus evolution it is necessary to meld experimental and comparative approaches within a single evolutionary genetic framework, and to link viral evolution at the intrahost scale with that which occurs over both epidemiological and geological timescales. We suggest that the impetus for this new synthesis may come from methodological advances in next-generation sequencing and metagenomics.

The endosymbiont bacterium Wolbachia is associated with multiple mutualistic effects on insect biology, including nutritional and antiviral properties. Wolbachia naturally occurs in Drosophila fly species, providing an operational model host to study how virome composition may be impacted by its presence. Drosophila simulans populations can carry a variety of Wolbachia strains. In particular, the wAu strain of Wolbachia has been associated with strong antiviral protection under experimental conditions. We used D. simulans sampled from the Perth Hills, Western Australia, to investigate the potential virus protective effect of the wAu strain on individual wild-caught flies. Our data revealed no appreciable variation in virus composition and abundance between Wolbachia infected/uninfected individuals associated with the presence/absence of wAu. However, it remains unclear whether wAu might impact viral infection and host survival by increasing tolerance rather than inducing complete resistance. These data also provide new insights into the natural virome diversity of D. simulans. Despite the small number of individuals sampled, we identified a repertoire of RNA viruses, including Nora virus, Galbut virus, Chaq virus, Thika virus and La Jolla virus, that have been identified in other Drosophila species. In addition, we identified five novel viruses from the families Reoviridae, Tombusviridae, Mitoviridae and Bunyaviridae. Overall, this study highlights the complex interaction between Wolbachia and RNA virus infections and provides a baseline description of the natural virome of D. simulans.

The ongoing pandemic spread of a new human coronavirus, SARS-CoV-2, which is associated with severe pneumonia/disease (COVID-19), has resulted in the generation of tens of thousands of virus genome sequences. The rate of genome generation is unprecedented, yet there is currently no coherent nor accepted scheme for naming the expanding phylogenetic diversity of SARS-CoV-2. Here, we present a rational and dynamic virus nomenclature that uses a phylogenetic framework to identify those lineages that contribute most to active spread. Our system is made tractable by constraining the number and depth of hierarchical lineage labels and by flagging and delabelling virus lineages that become unobserved and hence are probably inactive. By focusing on active virus lineages and those spreading to new locations, this nomenclature will assist in tracking and understanding the patterns and determinants of the global spread of SARS-CoV-2.

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.

Summary The unprecedented pandemic of pneumonia caused by a novel coronavirus, SARS-CoV-2, in China and beyond has had major public health impacts on a global scale [1,2]. Although bats are regarded as the most likely natural hosts for SARS-CoV-2 [3], the origins of the virus remain unclear. Here, we report a novel bat-derived coronavirus, denoted RmYN02, identified from a metagenomics analysis of samples from 227 bats collected from Yunnan Province in China between May and October, 2019. Notably, RmYN02 shares 93.3% nucleotide identity with SARS-CoV-2 at the scale of the complete virus genome and 97.2% identity in the 1ab gene, in which it is the closest relative of SARS-CoV-2 reported to date. In contrast, RmYN02 showed low sequence identity (61.3%) to SARS-CoV-2 in the receptor binding domain (RBD) and might not bind to angiotensin-converting enzyme 2 (ACE2). Critically, and in a similar manner to SARS-CoV-2, RmYN02 was characterized by the insertion of multiple amino acids at the junction site of the S1 and S2 subunits of the spike (S) protein. This provides strong evidence that such insertion events can occur naturally in animal betacoronaviruses.

Our knowledge of the diversity and evolution of the virosphere will likely increase dramatically with the study of microbial eukaryotes, including the microalgae within which few RNA viruses have been documented. By combining total RNA sequencing with sequence and structural-based homology detection, we identified 18 novel RNA viruses in cultured samples from two major groups of microbial algae: the chlorophytes and the chlorarachniophytes. Most of the RNA viruses identified in the green algae class Ulvophyceae were related to the Tombusviridae and Amalgaviridae viral families commonly associated with land plants. This suggests that the evolutionary history of these viruses extends to divergence events between algae and land plants. Seven Ostreobium sp-associated viruses exhibited sequence similarity to the mitoviruses most commonly found in fungi, compatible with horizontal virus transfer between algae and fungi. We also document, for the first time, RNA viruses associated with chlorarachniophytes, including the first negative-sense (bunya-like) RNA virus in microalgae, as well as a distant homolog of the plant virus Virgaviridae, potentially signifying viral inheritance from the secondary chloroplast endosymbiosis that marked the origin of the chlorarachniophytes. More broadly, these data suggest that the scarcity of RNA viruses in algae results from limited investigation rather than their absence.

Due to the scope and impact of the COVID-19 pandemic there exists a strong desire to understand where the SARS-CoV-2 virus came from and how it jumped species boundaries to humans. Molecular evolutionary analyses can trace viral origins by establishing relatedness and divergence times of viruses and identifying past selective pressures. However, we must uphold rigorous standards of inference and interpretation on this topic because of the ramifications of being wrong. Here, we dispute the conclusions of Xia (2020. Extreme genomic CpG deficiency in SARS-CoV-2 and evasion of host antiviral defense. Mol Biol Evol. doi:10.1093/molbev/masa095) that dogs are a likely intermediate host of a SARS-CoV-2 ancestor. We highlight major flaws in Xia’s inference process and his analysis of CpG deficiencies, and conclude that there is no direct evidence for the role of dogs as intermediate hosts. Bats and pangolins currently have the greatest support as ancestral hosts of SARS-CoV-2, with the strong caveat that sampling of wildlife species for coronaviruses has been limited.

In January 2020, a novel betacoronavirus (family Coronaviridae), named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was identified as the etiological agent of a cluster of pneumonia cases occurring in Wuhan City, Hubei Province, China1,2. The disease arising from SARS-CoV-2 infection, coronavirus disease 2019 (COVID-19), subsequently spread rapidly causing a worldwide pandemic. Here we examine the added value of near real-time genome sequencing of SARS-CoV-2 in a subpopulation of infected patients during the first 10 weeks of COVID-19 containment in Australia and compare findings from genomic surveillance with predictions of a computational agent-based model (ABM). Using the Australian census data, the ABM generates over 24 million software agents representing the population of Australia, each with demographic attributes of an anonymous individual. It then simulates transmission of the disease over time, spreading from specific infection sources, using contact rates of individuals within different social contexts. We report that the prospective sequencing of SARS-CoV-2 clarified the probable source of infection in cases where epidemiological links could not be determined, significantly decreased the proportion of COVID-19 cases with contentious links, documented genomically similar cases associated with concurrent transmission in several institutions and identified previously unsuspected links. Only a quarter of sequenced cases appeared to be locally acquired and were concordant with predictions from the ABM. These high-resolution genomic data are crucial to track cases with locally acquired COVID-19 and for timely recognition of independent importations once border restrictions are lifted and trade and travel resume.

Emerging infectious diseases, such as severe acute respiratory syndrome (SARS) and Zika virus disease, present a major threat to public health1–3. Despite intense research efforts, how, when and where new diseases appear are still a source of considerable uncertainty. A severe respiratory disease was recently reported in Wuhan, Hubei province, China. As of 25 January 2020, at least 1,975 cases had been reported since the first patient was hospitalized on 12 December 2019. Epidemiological investigations have suggested that the outbreak was associated with a seafood market in Wuhan. Here we study a single patient who was a worker at the market and who was admitted to the Central Hospital of Wuhan on 26 December 2019 while experiencing a severe respiratory syndrome that included fever, dizziness and a cough. Metagenomic RNA sequencing4 of a sample of bronchoalveolar lavage fluid from the patient identified a new RNA virus strain from the family Coronaviridae, which is designated here ‘WH-Human 1’ coronavirus (and has also been referred to as ‘2019-nCoV’). Phylogenetic analysis of the complete viral genome (29,903 nucleotides) revealed that the virus was most closely related (89.1% nucleotide similarity) to a group of SARS-like coronaviruses (genus Betacoronavirus, subgenus Sarbecovirus) that had previously been found in bats in China5. This outbreak highlights the ongoing ability of viral spill-over from animals to cause severe disease in humans., Phylogenetic and metagenomic analyses of the complete viral genome of a new coronavirus from the family Coronaviridae reveal that the virus is closely related to a group of SARS-like coronaviruses found in bats in China.

Here, we report the detection of a novel alphavirus in Australian mosquitoes, provisionally named Yada Yada virus (YYV). Phylogenetic analysis indicated that YYV belongs to the mosquito-specific alphavirus complex. The assembled genome is 11,612 nucleotides in length and encodes two open reading frames.

A recent article in BMC Genomics describes a new bioinformatics tool, HumanMycobiomeScan, to classify fungal taxa in metagenomic samples. This tool was used to characterize the gut mycobiome of hunter-gatherers and Western populations, resulting in the identification of a range of fungal species in the vast majority of samples. In the HumanMycobiomeScan pipeline, sequence reads are mapped against a reference database containing fungal genome sequences only. We argue that using reference databases comprised of a single taxonomic group leads to an unacceptably high number of false-positives due to: (i) mapping to conserved genetic regions in reference genomes, and (ii) sequence contamination in the assembled reference genomes. To demonstrate this, we replaced the HumanMycobiomeScan’s fungal reference database with one containing genome sequences of amphibians and reptiles and re-analysed their case study. The classification pipeline recovered all species present in the reference database, revealing turtles (Geoemydidae), bull frogs (Pyxicephalidae) and snakes (Colubridae) as the most abundant herpetological taxa in the human gut. We also re-analysed their case study using a kingdom-agnostic pipeline. This revealed that while the gut of hunter-gatherers and Western subjects may be colonized by a range of microbial eukaryotes, only three fungal families were retrieved. These results highlight the pitfalls of using taxon-specific reference databases for metagenome classification, even when they are comprised of curated whole genome data. We propose that databases containing all domains of life provide the most suitable option for metagenomic species profiling, especially when targeting microbial eukaryotes.

Summary Background In late December, 2019, patients presenting with viral pneumonia due to an unidentified microbial agent were reported in Wuhan, China. A novel coronavirus was subsequently identified as the causative pathogen, provisionally named 2019 novel coronavirus (2019-nCoV). As of Jan 26, 2020, more than 2000 cases of 2019-nCoV infection have been confirmed, most of which involved people living in or visiting Wuhan, and human-to-human transmission has been confirmed. Methods We did next-generation sequencing of samples from bronchoalveolar lavage fluid and cultured isolates from nine inpatients, eight of whom had visited the Huanan seafood market in Wuhan. Complete and partial 2019-nCoV genome sequences were obtained from these individuals. Viral contigs were connected using Sanger sequencing to obtain the full-length genomes, with the terminal regions determined by rapid amplification of cDNA ends. Phylogenetic analysis of these 2019-nCoV genomes and those of other coronaviruses was used to determine the evolutionary history of the virus and help infer its likely origin. Homology modelling was done to explore the likely receptor-binding properties of the virus. Findings The ten genome sequences of 2019-nCoV obtained from the nine patients were extremely similar, exhibiting more than 99·98% sequence identity. Notably, 2019-nCoV was closely related (with 88% identity) to two bat-derived severe acute respiratory syndrome (SARS)-like coronaviruses, bat-SL-CoVZC45 and bat-SL-CoVZXC21, collected in 2018 in Zhoushan, eastern China, but were more distant from SARS-CoV (about 79%) and MERS-CoV (about 50%). Phylogenetic analysis revealed that 2019-nCoV fell within the subgenus Sarbecovirus of the genus Betacoronavirus, with a relatively long branch length to its closest relatives bat-SL-CoVZC45 and bat-SL-CoVZXC21, and was genetically distinct from SARS-CoV. Notably, homology modelling revealed that 2019-nCoV had a similar receptor-binding domain structure to that of SARS-CoV, despite amino acid variation at some key residues. Interpretation 2019-nCoV is sufficiently divergent from SARS-CoV to be considered a new human-infecting betacoronavirus. Although our phylogenetic analysis suggests that bats might be the original host of this virus, an animal sold at the seafood market in Wuhan might represent an intermediate host facilitating the emergence of the virus in humans. Importantly, structural analysis suggests that 2019-nCoV might be able to bind to the angiotensin-converting enzyme 2 receptor in humans. The future evolution, adaptation, and spread of this virus warrant urgent investigation. Funding National Key Research and Development Program of China, National Major Project for Control and Prevention of Infectious Disease in China, Chinese Academy of Sciences, Shandong First Medical University.

Ticks (order: Ixodida) are a highly diverse and ecologically important group of ectoparasitic blood-feeding organisms. One such species, the seabird tick (Ixodes uriae), is widely distributed around the circumpolar regions of the northern and southern hemispheres. It has been suggested that Ix. uriae spread from the southern to the northern circumpolar region millions of years ago and has remained isolated in these regions ever since. Such a profound biographic subdivision provides a unique opportunity to determine whether viruses associated with ticks exhibit the same evolutionary patterns as their hosts. To test this, we collected Ix. uriae specimens near a Gentoo penguin (Pygoscelis papua) colony at Neko harbour, Antarctica, and from migratory birds—the Razorbill (Alca torda) and the Common murre (Uria aalge)—on Bonden island, northern Sweden. Through meta-transcriptomic next-generation sequencing we identified 16 RNA viruses, seven of which were novel. Notably, we detected the same species, Ronne virus, and two closely related species, Bonden virus and Piguzov virus, in both hemispheres indicating that there have been at least two cross-circumpolar dispersal events. Similarly, we identified viruses discovered previously in other locations several decades ago, including Gadgets Gully virus, Taggert virus and Okhotskiy virus. By identifying the same or closely related viruses in geographically disjunct sampling locations we provide evidence for virus dispersal within and between the circumpolar regions. In marked contrast, our phylogenetic analysis revealed no movement of the Ix. uriae tick hosts between the same locations. Combined, these data suggest that migratory birds are responsible for the movement of viruses at both local and global scales., Author summary As host populations diverge, so may those microorganisms, including viruses, that are dependent on those hosts. To examine this key issue in host-microbe evolution we compared the co-phylogenies of the seabird tick, Ixodes uriae, and their RNA viruses sampled from the far northern and southern hemispheres. Despite the huge geographic distance between them, phylogeographic analysis reveals that the same and closely related viruses were found both within and between the northern and southern circumpolar regions, most likely reflecting transfer by virus-infected migratory birds. In contrast, genomic data suggested that the Ix. uriae populations were phylogenetically distinct between the northern and southern hemispheres. This work emphasises the importance of migratory birds and ticks as vectors and sources of virus dispersal and introduction at both the local and global scales.

SummaryThe unprecedented epidemic of pneumonia caused by a novel coronavirus, HCoV-19, in China and beyond has caused public health concern at a global scale. Although bats are regarded as the most likely natural hosts for HCoV-191,2, the origins of the virus remain unclear. Here, we report a novel bat-derived coronavirus, denoted RmYN02, identified from a metagenomics analysis of samples from 227 bats collected from Yunnan Province in China between May and October, 2019. RmYN02 shared 93.3% nucleotide identity with HCoV-19 at the scale of the complete virus genome and 97.2% identity in the 1ab gene in which it was the closest relative of HCoV-19. In contrast, RmYN02 showed low sequence identity (61.3%) to HCoV-19 in the receptor binding domain (RBD) and might not bind to angiotensin-converting enzyme 2 (ACE2). Critically, however, and in a similar manner to HCoV-19, RmYN02 was characterized by the insertion of multiple amino acids at the junction site of the S1 and S2 subunits of the Spike (S) protein. This provides strong evidence that such insertion events can occur in nature. Together, these data suggest that HCoV-19 originated from multiple naturally occurring recombination events among those viruses present in bats and other wildlife species.

Background Wildlife species carry a remarkable diversity of trypanosomes. The detection of trypanosome infection in native Australian fauna is central to understanding their diversity and host-parasite associations. The implementation of total RNA sequencing (meta-transcriptomics) in trypanosome surveillance and diagnosis provides a powerful methodological approach to better understand the host species distribution of this important group of parasites. Methods We implemented a meta-transcriptomic approach to detect trypanosomes in a variety of tissues (brain, liver, lung, skin, gonads) sampled from native Australian wildlife, comprising four marsupials (koala, Phascolarctos cinereus; southern brown bandicoot, Isoodon obesulus; swamp wallaby, Wallabia bicolor; bare-nosed wombat, Vombatus ursinus), one bird (regent honeyeater, Anthochaera phrygia) and one amphibian (eastern dwarf tree frog, Litoria fallax). Samples corresponded to both clinically healthy and diseased individuals. Sequencing reads were de novo assembled into contigs and annotated. The evolutionary relationships among the trypanosomatid sequences identified were determined through phylogenetic analysis of 18S rRNA sequences. Results We detected trypanosome sequences in all six species of vertebrates sampled, with positive samples in multiple organs and tissues confirmed by PCR. Phylogenetic analysis indicated that the trypanosomes infecting marsupials were related to those previously detected in placental and marsupial mammals, while the trypanosome in the regent honeyeater grouped with avian trypanosomes. In contrast, we provide the first evidence for a trypanosome in the eastern dwarf tree frog that was phylogenetically distinct from those described in other amphibians. Conclusions To our knowledge, this is the first meta-transcriptomic analysis of trypanosomes in native Australian wildlife, expanding the known genetic diversity of these important parasites. We demonstrated that RNA sequencing is sufficiently sensitive to detect low numbers of Trypanosoma transcripts and from diverse hosts and tissues types, thereby representing an effective means to detect trypanosomes that are divergent in genome sequence.

Overlapping genes are commonplace in viruses and play an important role in their function and evolution. However, aside from studies on specific groups of viruses, relatively little is known about the extent and nature of gene overlap and its determinants in viruses as a whole. Here, we present an extensive characterisation of gene overlap in viruses through an analysis of reference genomes present in the NCBI virus genome database. We find that over half the instances of gene overlap are very small, covering

Little is known about the factors that shape the ecology of RNA viruses in nature. Wild birds are an important case in point, as other than influenza A virus, avian samples are rarely tested for viruses, especially in the absence of overt disease. Using bulk RNA-sequencing ("meta-transcriptomics"), we revealed the viral diversity present in Australian wild birds through the lens of the ecological factors that may determine virome structure and abundance. A meta-transcriptomic analysis of four Anseriformes (waterfowl) and Charadriiformes (shorebird) species sampled in temperate and arid Australia revealed the presence of 27 RNA virus genomes, 18 of which represent newly described species. The viruses identified included a previously described gammacoronavirus and influenza A viruses. Additionally, we identified novel virus species from the families Astroviridae, Caliciviridae, Reoviridae, Rhabdoviridae, Picobirnaviridae and Picornaviridae. We noted differences in virome structure that reflected underlying differences in location and influenza A infection status. Red-necked Avocets (Recurvirostra novaehollandiae) from Australia's arid interior possessed the greatest viral diversity and abundance, markedly higher than individuals sampled in temperate Australia. In Ruddy Turnstones (Arenaria interpres) and dabbling ducks (Anas spp.), viral abundance and diversity were higher and more similar in hosts that were positive for influenza A infection compared to those that were negative for this virus, despite samples being collected on the same day and from the same location. This study highlights the extent and diversity of RNA viruses in wild birds and lays the foundation for understanding the factors that determine virome structure in wild populations.

Overlapping genes in viruses maximize the coding capacity of their genomes and allow the generation of new genes without major increases in genome size. Despite their importance, the evolution and function of overlapping genes are often not well understood, in part due to difficulties in their detection. In addition, most bioinformatic approaches for the detection of overlapping genes require the comparison of multiple genome sequences that may not be available in metagenomic surveys of virus biodiversity. We introduce a simple new method for identifying candidate functional overlapping genes using single virus genome sequences. Our method uses randomization tests to estimate the expected length of open reading frames and then identifies overlapping open reading frames that significantly exceed this length and are thus predicted to be functional. We applied this method to 2548 reference RNA virus genomes and find that it has both high sensitivity and low false discovery for genes that overlap by at least 50 nucleotides. Notably, this analysis provided evidence for 29 previously undiscovered functional overlapping genes, some of which are coded in the antisense direction suggesting there are limitations in our current understanding of RNA virus replication.

Mutation. The word naturally conjures fears of unexpected and freakish changes. Ill-informed discussions of mutations thrive during virus outbreaks, including the ongoing spread of SARS-CoV-2. In reality, mutations are a natural part of the virus life cycle and rarely impact outbreaks dramatically.

The development of high-throughput genome sequencing enables accurate measurements of levels of sub-consensus intra-host virus genetic diversity and analysis of the role played by natural selection during cross-species transmission. We analysed the natural and experimental evolution of rabies virus (RABV), an important example of a virus that is able to make multiple host jumps. In particular, we (i) analyzed RABV evolution during experimental host switching with the goal of identifying possible genetic markers of host adaptation, (ii) compared the mutational changes observed during passage with those observed in natura, and (iii) determined whether the colonization of new hosts or tissues requires adaptive evolution in the virus. To address these aims, animal infection models (dog and fox) and primary cell culture models (embryo brain cells of dog and fox) were developed and viral variation was studied in detail through deep genome sequencing. Our analysis revealed a strong unidirectional host evolutionary effect, as dog-adapted rabies virus was able to replicate in fox and fox cells relatively easily, while dogs or neuronal dog cells were not easily susceptible to fox adapted-RABV. This suggests that dog RABV may be able to adapt to some hosts more easily than other host variants, or that when RABV switched from dogs to red foxes it lost its ability to adapt easily to other species. Although no difference in patterns of mutation variation between different host organs was observed, mutations were common following both in vitro and in vivo passage. However, only a small number of these mutations also appeared in natura, suggesting that adaptation during successful cross-species virus transmission is a complex, multifactorial evolutionary process., Author summary Understanding the mechanisms that underpin the cross-species transmission and host adaptation of rabies virus (RABV) remains an important part of the ongoing goal to reduce and eliminate rabies. We utilized next-generation sequencing to perform a deep comparative analysis of the genomic evolution of RABV subpopulations during host adaptation in culture and in animals, with the aim of determining the molecular mechanisms involved in the host-species or tissue adaptation of rabies virus. In particular, we aimed to determine whether experimental evolution can recapitulate evolution in nature. Our results suggest that a limited number of mutations that appeared following both in vitro and in vivo passage were observed in natura. This study also suggests that dog RABV may be able to adapt to some hosts more easily than other host variants.

Although RNA viruses are likely to exist in every species of cellular life, our knowledge of their biodiversity and evolution has been focused on those that can cause disease in humans and in economically important animals and plants. Recently published studies, however, have shown that every healthy organism can carry viruses. In this study, we present the unbiased discovery of highly divergent RNA viruses in European rodents and rabbits using Next Generation Sequencing (NGS). Tissue samples from different organs were collected from UK rodents and French rabbits and were initially screened for viruses with PCR. Following up on preliminary data, the positive samples were sent for NGS to acquire full genomes and perform unbiased virus discovery. Our findings encompass a number of novel viruses including astrovirus, rotavirus A, hantavirus, picornavirus, coronavirus, paramyxovirus etc. Among the novel viruses, the picornavirus interestingly showed 50 % and 34 % similarity to the closest relative for the non-structural and the structural protein respectively at amino acid level, suggesting a potentially novel genus within the family. This study shows the presence of highly divergent RNA viruses in European rodents and rabbits enriching thus our current knowledge of their origin and evolution. Finally, this study shows that these animal species can be the reservoirs of RNA viruses notorious for cross-species transmission, increasing the risk of a spill-over to humans or livestock.

Hepatitis delta virus (HDV) is the smallest known RNA virus, encoding a single protein. Until recently, HDV had only been identified in humans, where it is strongly associated with co-infection with hepatitis B virus (HBV). However, the recent discovery of HDV-like viruses in metagenomic samples from birds and snakes suggests that this virus has a far longer evolutionary history. Herein, using additional meta-transcriptomic data, we show that highly divergent HDV-like viruses are also present in fish, amphibians, and invertebrates, with PCR and Sanger sequencing confirming the presence of the invertebrate HDV-like viruses. Notably, the novel viruses identified here share genomic features characteristic of HDV, such as a circular genome of only approximately 1.7 kb in length, and self-complementary, unbranched rod-like structures. Coiled-coil domains, leucine zippers, conserved residues with essential biological functions, and isoelectronic points similar to those in the human hepatitis delta virus antigens (HDAgs) were also identified in the putative non-human viruses. Importantly, none of these novel HDV-like viruses were associated with hepadnavirus infection, supporting the idea that the HDV–HBV association may be specific to humans. Collectively, these data not only broaden our understanding of the diversity and host range of HDV, but also shed light on its origin and evolutionary history.

Emerging viruses have the potential to impose substantial mortality, morbidity and economic burdens on human populations. Tracking the spread of infectious diseases to assist in their control has traditionally relied on the analysis of case data gathered as the outbreak proceeds. Here, we describe how many of the key questions in infectious disease epidemiology, from the initial detection and characterization of outbreak viruses, to transmission chain tracking and outbreak mapping, can now be much more accurately addressed using recent advances in virus sequencing and phylogenetics. We highlight the utility of this approach with the hypothetical outbreak of an unknown pathogen, 'Disease X', suggested by the World Health Organization to be a potential cause of a future major epidemic. We also outline the requirements and challenges, including the need for flexible platforms that generate sequence data in real-time, and for these data to be shared as widely and openly as possible.

Two closely related caliciviruses cocirculate in Australia: rabbit hemorrhagic disease virus (RHDV) and rabbit calicivirus Australia 1 (RCV-A1). RCV-A1 causes benign enteric infections in the European rabbit ( Oryctolagus cuniculus ) in Australia and New Zealand, while its close relative RHDV causes a highly pathogenic infection of the liver in the same host. The comparison of these viruses provides important information on the nature and trajectory of virulence evolution, particularly as highly virulent strains of RHDV may have evolved from nonpathogenic ancestors such as RCV-A1. To determine the evolution of RCV-A1 we sequenced the full-length genomes of 44 RCV-A1 samples isolated from healthy rabbits and compared key evolutionary parameters to those of its virulent relative, RHDV. Despite their marked differences in pathogenicity and tissue tropism, RCV-A1 and RHDV have evolved in a very similar manner. Both viruses have evolved at broadly similar rates, suggesting that their dynamics are largely shaped by high background mutation rates, and both exhibit occasional recombination and an evolutionary environment dominated by purifying selection. In addition, our comparative analysis revealed that there have been multiple changes in both virulence and tissue tropism in the evolutionary history of these and related viruses. Finally, these new genomic data suggest that either RCV-A1 was introduced into Australia after the introduction of myxoma virus as a biocontrol agent in 1950 or there was drastic reduction of the rabbit population, and hence of RCV-A1 genetic diversity, perhaps coincident with the emergence of myxoma virus. IMPORTANCE The comparison of closely related viruses that differ profoundly in propensity to cause disease in their hosts offers a powerful opportunity to reveal the causes of changes in virulence and to study how such changes alter the evolutionary dynamics of these pathogens. Here we describe such a novel comparison involving two closely related RNA viruses that cocirculate in Australia, the highly virulent rabbit hemorrhagic disease virus (RHDV) and the nonpathogenic rabbit calicivirus Australia 1 (RCV-A1). Both viruses infect the European rabbit, but they differ in virulence, tissue tropism, and mechanisms of transmission. Surprisingly, and despite these fundamental differences, RCV-A1 and RHDV have evolved at very similar (high) rates and with strong purifying selection. Furthermore, candidate key mutations were identified that may play a role in virulence and/or tissue tropism and therefore warrant further investigation.

The emergence and reemergence of rapidly evolving RNA viruses—particularly those responsible for respiratory diseases, such as influenza viruses and coronaviruses—pose a significant threat to global health, including the potential of major pandemics. Importantly, recent advances in high-throughput genome sequencing enable researchers to reveal the genomic diversity of these viral pathogens at much lower cost and with much greater precision than they could before. In particular, the genome sequence data generated allow inferences to be made on the molecular basis of viral emergence, evolution, and spread in human populations in real time. In this review, we introduce recent computational methods that analyze viral genomic data, particularly in combination with metadata such as sampling time, geographic location, and virulence. We then outline the insights these analyses have provided into the fundamental patterns and processes of evolution and emergence in human respiratory RNA viruses, as well as the major challenges in such genomic analyses.

Influenza A virus is characterized by high genetic diversity. However, most of what is known about influenza evolution has come from consensus sequences sampled at the epidemiological scale that only represent the dominant virus lineage within each infected host. Less is known about the extent of within-host virus diversity and what proportion of this diversity is transmitted between individuals. To characterize virus variants that achieve sustainable transmission in new hosts, we examined within-host virus genetic diversity in household donor-recipient pairs from the first wave of the 2009 H1N1 pandemic when seasonal H3N2 was co-circulating. Although the same variants were found in multiple members of the community, the relative frequencies of variants fluctuated, with patterns of genetic variation more similar within than between households. We estimated the effective population size of influenza A virus across donor-recipient pairs to be approximately 100-200 contributing members, which enabled the transmission of multiple lineages, including antigenic variants.

Eukaryotes of the genus Plasmodium cause malaria, a parasitic disease responsible for substantial morbidity and mortality in humans. Yet, the nature and abundance of any viruses carried by these divergent eukaryotic parasites is unknown. We investigated the Plasmodium virome by performing a meta-transcriptomic analysis of blood samples taken from patients suffering from malaria and infected with P. vivax, P. falciparum or P. knowlesi. This resulted in the identification of a narnavirus-like sequence, encoding an RNA polymerase and restricted to P. vivax samples, as well as an associated viral segment of unknown function. These data, confirmed by PCR, are indicative of a novel RNA virus that we term Matryoshka RNA virus 1 (MaRNAV-1) to reflect its analogy to a "Russian doll": a virus, infecting a parasite, infecting an animal. Additional screening revealed that MaRNAV-1 was abundant in geographically diverse P. vivax derived from humans and mosquitoes, strongly supporting its association with this parasite, and not in any of the other Plasmodium samples analyzed here nor Anopheles mosquitoes in the absence of Plasmodium. Notably, related bi-segmented narnavirus-like sequences (MaRNAV-2) were retrieved from Australian birds infected with a Leucocytozoon—a genus of eukaryotic parasites that group with Plasmodium in the Apicomplexa subclass hematozoa. Together, these data support the establishment of two new phylogenetically divergent and genomically distinct viral species associated with protists, including the first virus likely infecting Plasmodium parasites. As well as broadening our understanding of the diversity and evolutionary history of the eukaryotic virosphere, the restriction to P. vivax may be of importance in understanding P. vivax-specific biology in humans and mosquitoes, and how viral co-infection might alter host responses at each stage of the P. vivax life-cycle., Author summary While parasites are a major cause of human disease, they can themselves be infected by viruses. We asked whether three of the major malaria-causing parasites in humans—Plasmodium vivax, P. falciparum and P. knowlesi—were also infected by viruses. To this end we performed total RNA-Sequencing (“meta-transcriptomics”) on human blood samples infected with these Plasmodium species. This resulted in the discovery of an abundant bi-segmented virus—Matryoshka RNA virus 1 (MaRNAV-1)—in all P. vivax samples tested (but no other Plasmodium species) that contains a replicase segment related to those of narnaviruses, arguably the simplest type of RNA viruses discovered to date. By screening for MaRNAV-1 in a larger set of Plasmodium species we revealed a strong specificity between this virus and P. vivax, as well as the presence of a related virus—MaRNAV-2—in avian Leucocytozoon hematozoa parasites. This is the first discovery of a Plasmodium-associated virus and will assist in revealing the deep evolutionary history of RNA viruses and our understanding of Plasmodium biology and disease processes.

Dengue virus (DENV) transmission by mosquitoes is a time-dependent process that begins with the consumption of an infectious blood-meal. DENV infection then proceeds stepwise through the mosquito from the midgut to the carcass, and ultimately to the salivary glands, where it is secreted into saliva and then transmitted anew on a subsequent bite. We examined viral kinetics in tissues of the Aedes aegypti mosquito over a finely graded time course, and as per previous studies, found that initial viral dose and serotype strain diversity control infectivity. We also found that a threshold level of virus is required to establish body-wide infections and that replication kinetics in the early and intermediate tissues do not predict those of the salivary glands. Our findings have implications for mosquito GMO design, modeling the contribution of transmission to vector competence and the role of mosquito kinetics in the overall DENV epidemiological landscape., Author summary DENV infection in the mosquito is a complex and dynamic process. Following ingestion of an infected blood meal, DENV enters the mosquito midgut epithelial cells, where it replicates. Subsequently, the virus disseminates and infects other tissues, including hemocytes, fat body and reproductive organs, ultimately reaching the salivary glands. The kinetics of infection are influenced by genetic variation in the virus. Comparisons between strains within single serotypes, have revealed variation in infection rates in mosquitoes. To explore the role of infectious dose, serotype and tissue in viral infection kinetics we sampled DENV loads in populations of infected mosquitoes over numerous, sequential time-points. We reveal that the kinetics of DENV infection in the midgut, carcass and salivary glands of the mosquito Aedes aegypti are strikingly different among the strains selected for this study, and that these differences are also driven by the initial infectious dose.

Lymphoma is one of the most common malignancies in domestic cats. The lymphomagenicpotential of Felis catus gammaherpesvirus 1 (FcaGHV1), a common infection in domestic cats,is unknown. In other species, including humans, cellular transformation by gammaherpesvirusesis typically mediated by viral genes expressed during latency. We analysed tumour RNA, fromdiffuse large B-cell lymphomas (DLBCL) appearing in cats coinfected with FcaGHV1 and felineimmunodeficiency virus (FIV) (n = 10). Analysis was done by high throughput transcriptomesequencing and reverse transcription PCR. A limited repertoire of FcaGHV transcripts was identifiedin five tumors, including homologs of oncogenic latency-associated transcripts, latency-associatednuclear antigen (LANA, ORF73) and vFLIP (F7), lytic genes (ORF50, ORF6, ORF59, F10), and an ORFunique to FcaGHV1, F20. In situ hybridization of FIV-associated DLBCLs (n = 9), post-transplantlymphomas (n = 6) and high-grade B and T-cell intestinal lymphomas (n = 8) identified a single casein which FcaGHV1 nucleic acid was detectable. These results demonstrate that FcaGHV1 transcriptscan be detected in some FIV-associated lymphomas, but with a low copy number, precludingassessment of a potential role for FcaGHV1 in lymphomagenesis. Future investigation of the FcaGHV1transcriptome in clinical samples might employ viral enrichment and greater sequencing depth toenhance the retrieval of viral reads. Our results suggest prioritization of a subset of intestinal T-celltumors and large granular lymphocyte lymphoma, for study.

Wolbachia is an endosymbiotic bacterium that can block viral infections in arthropods, generating interest in its potential to control the spread of mosquito-borne disease. Drosophila melanogaster is model organism for Wolbachia infection, and the w Mel strain of Wolbachia can improve host survival following viral infection. However, it is unclear whether w Mel induces anti-viral blocking against the broader native virome of D. melanogaster , or whether the major effect of Wolbachia is a reduction in viral abundance rather than viral clearance. We examined the effect of Wolbachia on viral abundance by comparing the total transcriptome of w Mel-positive and w Mel-negative D. melanogaster populations sampled from six locations in Australia. In addition, we examined the impact of w Mel on individual flies by obtaining transcriptome data from 20 w Mel-positive and 20 w Mel-negative D. melanogaster from the location (Melbourne) with highest density of w Mel. These data revealed high viral abundance in both Wolbachia -positive and -negative populations and individuals. Notably, none of the viral species identified, representing RNA viruses from at least nine families/floating genera, showed evidence of protection by w Mel. Although the viral loads of picorna-like viruses are reduced by w Mel under experimental conditions, we observed no such effect here. These data show that D. melanogaster can harbour abundant RNA viruses regardless of its Wolbachia status and imply that the interaction between Wolbachia and viruses in nature is more complex than simple blocking.

Yersinia pestis , the causative agent of plague, evolved from the closely related pathogen Yersinia pseudotuberculosis . During its emergence, Y. pestis is believed to have acquired its unique pathogenic characteristics through numerous gene gains/losses, genomic rearrangements, and single nucleotide polymorphism (SNP) changes. One such SNP creates a single amino acid variation in the DNA binding domain of PhoP, the response regulator in the PhoP/PhoQ two-component system. Y. pseudotuberculosis and the basal human-avirulent strains of Y. pestis harbor glycines at position 215 of PhoP, whereas the modern human-virulent strains (e.g., KIM and CO92) harbor serines at this residue. Since PhoP plays multiple roles in the adaptation of Y. pestis to stressful host conditions, we tested whether this amino acid substitution affects PhoP activity or the ability of Y. pestis to survive in host environments. Compared to the parental KIM6+ strain carrying the modern allele of phoP ( phoP-S215 ), a derivative carrying the basal allele ( phoP-G215 ) exhibited slightly defective growth under a low-Mg 2+ condition and decreased transcription of a PhoP target gene, ugd , as well as an ∼8-fold increase in the susceptibility to the antimicrobial peptide polymyxin B. The phoP-G215 strain showed no apparent defect in flea colonization, although a phoP -null mutant showed decreased flea infectivity in competition experiments. Our results suggest that the amino acid variation at position 215 of PhoP causes subtle changes in the PhoP activity and raise the possibility that the change in this residue have contributed to the evolution of increased virulence in Y. pestis . IMPORTANCE Y. pestis acquired a single nucleotide polymorphism (SNP) in phoP when the highly human-virulent strains diverged from less virulent basal strains, resulting in an amino acid substitution in the DNA binding domain of the PhoP response regulator. We show that Y. pestis carrying the modern phoP allele has an increased ability to induce the PhoP-regulated ugd gene and resist antimicrobial peptides compared to an isogenic strain carrying the basal allele. Given the important roles PhoP plays in host adaptation, the results raise an intriguing possibility that this amino acid substitution contributed to the evolution of increased virulence in Y. pestis . Additionally, we present the first evidence that phoP confers a survival fitness advantage to Y. pestis inside the flea midgut.

The factors that determine the pattern and rate of spread of influenza virus at a continental-scale are uncertain. Although recent work suggests that influenza epidemics in the United States exhibit a strong geographical correlation, the spatiotemporal dynamics of influenza in Australia, a country and continent of approximately similar size and climate complexity but with a far smaller population, are not known. Using a unique combination of large-scale laboratory-confirmed influenza surveillance comprising >450,000 entries and genomic sequence data we determined the local-level spatial diffusion of this important human pathogen nationwide in Australia. We used laboratory-confirmed influenza data to characterize the spread of influenza virus across Australia during 2007–2016. The onset of established epidemics varied across seasons, with highly synchronized epidemics coinciding with the emergence of antigenically distinct viruses, particularly during the 2009 A/H1N1 pandemic. The onset of epidemics was largely synchronized between the most populous cities, even those separated by distances of >3000 km and those that experience vastly diverse climates. In addition, by analyzing global phylogeographic patterns we show that the synchronized dissemination of influenza across Australian cities involved multiple introductions from the global influenza population, coupled with strong domestic connectivity, rather than through the distinct radial patterns of geographic dispersal that are driven by work-flow transmission as observed in the United States. In addition, by comparing the spatial structure of influenza A and B, we found that these viruses tended to occupy different geographic regions, and peak in different seasons, perhaps indicative of moderate cross-protective immunity or viral interference effects. The highly synchronized outbreaks of influenza virus at a continental-scale revealed here highlight the importance of coordinated public health responses in the event of the emergence of a novel, human-to-human transmissible, virus., Author summary We use a unique combination of large-scale laboratory-confirmed influenza surveillance and genomic sequence data to determine the pattern and determinants of influenza spread through time and space at a continental scale, utilizing 450,000 data entries gathered across Australia over a 10-year period (2007–2016). Our results reveal a remarkable epidemiological synchronicity across this large and geographically diverse continent, and in the face of enormous climatic variation, from the tropical north to the temperate south. Such synchronicity was especially marked during years associated with the emergence of antigenically distinct strains. This pattern of continental synchronicity highlights the importance of a highly coordinated responses in the event of the emergence of a novel, human-to-human transmissible, virus, and hence has broad-scale public health implications.

Molecular clock models are commonly used to estimate evolutionary rates and timescales from nucleotide sequences. The goal of these models is to account for rate variation among lineages, such that they are assumed to be adequate descriptions of the processes that generated the data. A common approach for selecting a clock model for a data set of interest is to examine a set of candidates and to select the model that provides the best statistical fit. However, this can lead to unreliable estimates if all the candidate models are actually inadequate. For this reason, a method of evaluating absolute model performance is critical. We describe a method that uses posterior predictive simulations to assess the adequacy of clock models. We test the power of this approach using simulated data and find that the method is sensitive to bias in the estimates of branch lengths, which tends to occur when using underparameterized clock models. We also compare the performance of the multinomial test statistic, originally developed to assess the adequacy of substitution models, but find that it has low power in identifying the adequacy of clock models. We illustrate the performance of our method using empirical data sets from coronaviruses, simian immunodeficiency virus, killer whales, and marine turtles. Our results indicate that methods of investigating model adequacy, including the one proposed here, should be routinely used in combination with traditional model selection in evolutionary studies. This will reveal whether a broader range of clock models to be considered in phylogenetic analysis.

Enterovirus A71 (EV-A71) has emerged as the most important cause of large outbreaks of severe and sometimes fatal hand, foot and mouth disease (HFMD) across the Asia-Pacific region. EV-A71 outbreaks have been associated with (sub)genogroup switches, sometimes accompanied by recombination events. Understanding EV-A71 population dynamics is therefore essential for understanding this emerging infection, and may provide pivotal information for vaccine development. Despite the public health burden of EV-A71, relatively few EV-A71 complete-genome sequences are available for analysis and from limited geographical localities. The availability of an efficient procedure for whole-genome sequencing would stimulate effort to generate more viral sequence data. Herein, we report for the first time the development of a next-generation sequencing based protocol for whole-genome sequencing of EV-A71 directly from clinical specimens. We were able to sequence viruses of subgenogroup C4 and B5, while RNA from culture materials of diverse EV-A71 subgenogroups belonging to both genogroup B and C was successfully amplified. The nature of intra-host genetic diversity was explored in 22 clinical samples, revealing 107 positions carrying minor variants (ranging from 0 to 15 variants per sample). Our analysis of EV-A71 strains sampled in 2013 showed that they all belonged to subgenogroup B5, representing the first report of this subgenogroup in Vietnam. In conclusion, we have successfully developed a high-throughput next-generation sequencing-based assay for whole-genome sequencing of EV-A71 from clinical samples.

Although rodents are important reservoirs for RNA viruses, to date only one species of rodent coronavirus (CoV) has been identified. Herein, we describe a new CoV, denoted Lucheng Rn rat coronavirus (LRNV), and novel variants of two Betacoronavirus species termed Longquan Aa mouse coronavirus (LAMV) and Longquan Rl rat coronavirus (LRLV), that were identified in a survey of 1465 rodents sampled in China during 2011–2013. Phylogenetic analysis revealed that LAMV and LRLV fell into lineage A of the genus Betacoronavirus, which included CoVs discovered in humans and domestic and wild animals. In contrast, LRNV harbored by Rattus norvegicus formed a distinct lineage within the genus Alphacoronavirus in the 3CLpro, RdRp, and Hel gene trees, but formed a more divergent lineage in the N and S gene trees, indicative of a recombinant origin. Additional recombination events were identified in LRLV. Together, these data suggest that rodents may carry additional unrecognized CoVs., Highlights • A novel coronavirus (CoV) was identified in rodents from China. • Three species of CoVs are co-circulating in rodents in Zhejiang province, China. • Recombination events were observed in the newly identified rodent CoVs.

Koolpinyah virus (KOOLV) isolated from healthy Australian cattle and Yata virus (YATV) isolated from a pool of Mansonia uniformis mosquitoes in the Central African Republic have been tentatively identified as rhabdoviruses. KOOLV was shown previously to be related antigenically to kotonkon virus, an ephemerovirus that has caused an ephemeral fever-like illness in cattle in Nigeria, but YATV failed to react antigenically with any other virus tested. Here we report the complete genome sequences of KOOLV (16,133 nt) and YATV (14,479 nt). Each has a complex genome organisation, with multiple genes, including a second non-structural glycoprotein (GNS) gene and a viroporin (α1) gene, between the G and L genes as is characteristic of ephemeroviruses. Based on an analysis of genome organisation, sequence identity and cross-neutralisation, we demonstrate that both KOOLV and YATV should be classified as two new species in the genus Ephemerovirus.

The coevolution of myxoma virus (MYXV) and wild European rabbits in Australia and Europe is a paradigm for the evolution of a pathogen in a new host species. Genomic analyses have identified the mutations that have characterized this evolutionary process, but defining causal mutations in the pathways from virulence to attenuation and back to virulence has not been possible. Using reverse genetics, we examined the roles of six selected mutations found in Australian field isolates of MYXV that fall in known or potential virulence genes. Several of these mutations occurred in genes previously identified as virulence genes in whole-gene knockout studies. Strikingly, no single or double mutation among the mutations tested had an appreciable impact on virulence. This suggests either that virulence evolution was defined by amino acid changes other than those analyzed here or that combinations of multiple mutations, possibly involving epistatic interactions or noncoding sequences, have been critical in the ongoing evolution of MYXV virulence. In sum, our results show that single-gene knockout studies of a progenitor virus can have little power to predict the impact of individual mutations seen in the field. The genetic determinants responsible for this canonical case of virulence evolution remain to be determined. IMPORTANCE The species jump of myxoma virus (MYXV) from the South American tapeti to the European rabbit populations of Australia and Europe is a canonical example of host-pathogen coevolution. Detailed molecular studies have identified multiple genes in MYXV that are critical for virulence, and genome sequencing has revealed the evolutionary history of MYXV in Australia and Europe. However, it has not been possible to categorically identify the key mutations responsible for the attenuation of or reversion to virulence during this evolutionary process. Here we use reverse genetics to examine the role of mutations in viruses isolated early and late in the Australian radiation of MYXV. Surprisingly, none of the candidate mutations that we identified as likely having roles in attenuation proved to be important for virulence. This indicates that considerable caution is warranted when interpreting the possible role of individual mutations during virulence evolution.

Understanding how microbiomes affect host resistance, parasite virulence, and parasite-associated diseases requires a collaborative effort between parasitologists, microbial ecologists, virologists, and immunologists. We hereby propose the Parasite Microbiome Project to bring together researchers with complementary expertise and to study the role of microbes in host-parasite interactions., Understanding how microbiomes affect host resistance, parasite virulence, and parasite-associated diseases requires a collaborative effort between parasitologists, microbial ecologists, virologists, and immunologists. We hereby propose the Parasite Microbiome Project to bring together researchers with complementary expertise and to study the role of microbes in host-parasite interactions. Data from the Parasite Microbiome Project will help identify the mechanisms driving microbiome variation in parasites and infected hosts and how that variation is associated with the ecology and evolution of parasites and their disease outcomes. This is a call to arms to prevent fragmented research endeavors, encourage best practices in experimental approaches, and allow reliable comparative analyses across model systems. It is also an invitation to foundations and national funding agencies to propel the field of parasitology into the microbiome/metagenomic era.

Viruses use the cellular machinery of their hosts for replication. It has therefore been proposed that the nucleotide and dinucleotide compositions of viruses should match those of their host species. If this is upheld, it may then be possible to use dinucleotide composition to predict the true host species of viruses sampled in metagenomic surveys. However, it is also clear that different taxonomic groups of viruses tend to have distinctive patterns of dinucleotide composition that may be independent of host species. To determine the relative strength of the effect of host versus virus family in shaping dinucleotide composition, we performed a comparative analysis of 20 RNA virus families from 15 host groupings, spanning two animal phyla and more than 900 virus species. In particular, we determined the odds ratios for the 16 possible dinucleotides and performed a discriminant analysis to evaluate the capability of virus dinucleotide composition to predict the correct virus family or host taxon from which it was isolated. Notably, while 81% of the data analyzed here were predicted to the correct virus family, only 62% of these data were predicted to their correct subphylum/class host and a mere 32% to their correct mammalian order. Similarly, dinucleotide composition has a weak predictive power for different hosts within individual virus families. We therefore conclude that dinucleotide composition is generally uniform within a virus family but less well reflects that of its host species. This has obvious implications for attempts to accurately predict host species from virus genome sequences alone. IMPORTANCE Determining the processes that shape virus genomes is central to understanding virus evolution and emergence. One question of particular importance is why nucleotide and dinucleotide frequencies differ so markedly between viruses. In particular, it is currently unclear whether host species or virus family has the biggest impact on dinucleotide frequencies and whether dinucleotide composition can be used to accurately predict host species. Using a comparative analysis, we show that dinucleotide composition has a strong phylogenetic association across different RNA virus families, such that dinucleotide composition can predict the family from which a virus sequence has been isolated. Conversely, dinucleotide composition has a poorer predictive power for the different host species within a virus family and across different virus families, indicating that the host has a relatively small impact on the dinucleotide composition of a virus genome.