Comparative Analysis of RNA Virome Composition in Rabbits and Associated Ectoparasites.

Ectoparasites play an important role in virus transmission among vertebrates. Little, however, is known about the nature of those viruses that pass between invertebrates and vertebrates. In Australia, flies and fleas support the mechanical transmission of two viral biological controls against wild rabbits-rabbit hemorrhagic disease virus (RHDV) and myxoma virus. We compared virome compositions in rabbits and these ectoparasites, sequencing total RNA from multiple tissues and gut contents of wild rabbits, fleas collected from these rabbits, and flies trapped sympatrically. Meta-transcriptomic analyses identified 50 novel viruses from multiple RNA virus families. Rabbits and their ectoparasites were characterized by markedly different viromes, with virus abundance greatest in flies. Although viral contigs from six virus families/groups were found in both rabbits and ectoparasites, they clustered in distinct host-dependent lineages. A novel calicivirus and a picornavirus detected in rabbit cecal content were vertebrate specific; the newly detected calicivirus was distinct from known rabbit caliciviruses, while the picornavirus clustered with sapeloviruses. Several picobirnaviruses were also identified that fell in diverse phylogenetic positions, compatible with the idea that they are associated with bacteria. Further comparative analysis revealed that the remaining viruses found in rabbits, and all those from ectoparasites, were likely associated with invertebrates, plants, and coinfecting endosymbionts. While no full genomes of vertebrate-associated viruses were detected in ectoparasites, small numbers of reads from rabbit astrovirus, RHDV, and other lagoviruses were present in flies. This supports a role for flies in the mechanical transmission of RHDV, while their involvement in astrovirus transmission merits additional exploration.IMPORTANCE Ectoparasites play an important role in the transmission of many vertebrate-infecting viruses, including Zika and dengue viruses. Although it is becoming increasingly clear that invertebrate species harbor substantial virus diversity, it is unclear how many of the viruses carried by invertebrates have the potential to infect vertebrate species. We used the European rabbit (Oryctolagus cuniculus) as a model species to compare virome compositions in a vertebrate host and known associated ectoparasite mechanical vectors, in this case, fleas and blowflies. In particular, we aimed to infer the extent of viral transfer between these distinct types of host. Our analysis revealed that despite extensive viral diversity in both rabbits and associated ectoparasites, and the close interaction of these vertebrate and invertebrate species, biological viral transmission from ectoparasites to vertebrate species is rare. We did, however, find evidence to support the idea of a role of blowflies in transmitting viruses without active replication in the insect.

Rabbit Hemorrhagic Disease Virus 2 (RHDV2; GI.2) Is Replacing Endemic Strains of RHDV in the Australian Landscape within 18 Months of Its Arrival.

Rabbit hemorrhagic disease virus 2 (RHDV2; Lagovirus GI.2) is a pathogenic calicivirus that affects European rabbits (Oryctolagus cuniculus) and various hare (Lepus) species. GI.2 was first detected in France in 2010 and subsequently caused epidemics in wild and domestic lagomorph populations throughout Europe. In May 2015, GI.2 was detected in Australia. Within 18 months of its initial detection, GI.2 had spread to all Australian states and territories and rapidly became the dominant circulating strain, replacing Rabbit hemorrhagic disease virus (RHDV/GI.1) in mainland Australia. Reconstruction of the evolutionary history of 127 Australian GI.2 isolates revealed that the virus arrived in Australia at least several months before its initial description and likely circulated unnoticed in wild rabbit populations in the east of the continent prior to its detection. GI.2 sequences isolated from five hares clustered with sequences from sympatric rabbit populations sampled contemporaneously, indicating multiple spillover events into hares rather than an adaptation of the Australian GI.2 to a new host. Since the presence of GI.2 in Australia may have wide-ranging consequences for rabbit biocontrol, particularly with the release of the novel biocontrol agent GI.1a/RHDVa-K5 in March 2017, ongoing surveillance is critical to understanding the interactions of the various lagoviruses in Australia and their impact on host populations.IMPORTANCE This study describes the spread and distribution of Rabbit hemorrhagic disease virus 2 (GI.2) in Australia since its first detection in May 2015. Within the first 18 months following its detection, RHDV2 spread from east to west across the continent and became the dominant strain in all mainland states of Australia. This has important implications for pest animal management and for owners of pet and farmed rabbits, as there currently is no effective vaccine available in Australia for GI.2. The closely related RHDV (GI.1) is used to control overabundant wild rabbits, a serious environmental and agricultural pest in this country, and it is currently unclear how the widespread circulation of GI.2 will impact ongoing targeted wild rabbit management operations.

Detection and Circulation of a Novel Rabbit Hemorrhagic Disease Virus in Australia.

The highly virulent rabbit hemorrhagic disease virus (RHDV) has been widely used in Australia and New Zealand since the mid-1990s to control wild rabbits, an invasive vertebrate pest in these countries. In January 2014, an exotic RHDV was detected in Australia, and 8 additional outbreaks were reported in both domestic and wild rabbits in the 15 months following its detection. Full-length genomic analysis revealed that this virus is a recombinant containing an RHDVa capsid gene and nonstructural genes most closely related to nonpathogenic rabbit caliciviruses. Nationwide monitoring efforts need to be expanded to assess if the increasing number of different RHDV variants circulating in the Australian environment will affect biological control of rabbits. At the same time, updated vaccines and vaccination protocols are urgently needed to protect pet and farmed rabbits from these novel rabbit caliciviruses.

Benign Rabbit Caliciviruses Exhibit Evolutionary Dynamics Similar to Those of Their Virulent Relatives.

UNLABELLED: 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.

Benign Rabbit Calicivirus in New Zealand.

The Czech v351 strain of rabbit hemorrhagic disease virus (RHDV1) is used in Australia and New Zealand as a biological control agent for rabbits, which are important and damaging introduced vertebrate pests in these countries. However, nonpathogenic rabbit caliciviruses (RCVs) can provide partial immunological cross-protection against lethal RHDV infection and thus interfere with effective rabbit biocontrol. Antibodies that cross-reacted against RHDV antigens were found in wild rabbits before the release of RHDV1 in New Zealand in 1997, suggesting that nonpathogenic RCVs were already present in New Zealand. The aim of this study was to confirm the presence of nonpathogenic RCV in New Zealand and describe its geographical distribution. RCV and RHDV antibody assays were used to screen serum samples from 350 wild rabbits from 14 locations in New Zealand. The serological survey indicated that both RCV and RHDV are widespread in New Zealand wild rabbits, with antibodies detected in 10 out of 14 and 12 out of 14 populations, respectively. Two closely related RCV strains were identified in the duodenal tissue from a New Zealand wild rabbit (RCV Gore-425A and RCV Gore-425B). Both variants are most closely related to Australian RCV strains, but with 88% nucleotide identity, they are genetically distinct. Phylogenetic analysis revealed that the New Zealand RCV strains fall within the genetic diversity of the Australian RCV isolates, indicating a relatively recent movement of RCVs between Australia and New Zealand.IMPORTANCE Wild rabbits are important and damaging introduced vertebrate pests in Australia and New Zealand. Although RHDV1 is used as a biological control agent, some nonpathogenic RCVs can provide partial immunological cross-protection against lethal RHDV infection and thus interfere with its effectiveness for rabbit control. The presence of nonpathogenic RCVs in New Zealand wild rabbits has been long hypothesized, but earlier attempts to isolate a New Zealand RCV strain have been unsuccessful. Therefore, it is important to determine if such nonpathogenic viruses exist in New Zealand rabbits, especially considering the proposed introduction of new RHDV strains into New Zealand as biocontrols.

Comparative Phylodynamics of Rabbit Hemorrhagic Disease Virus in Australia and New Zealand.

UNLABELLED: The introduction of rabbit hemorrhagic disease virus (RHDV) into Australia and New Zealand during the 1990s as a means of controlling feral rabbits is an important case study in viral emergence. Both epidemics are exceptional in that the founder viruses share an origin and the timing of their release is known, providing a unique opportunity to compare the evolution of a single virus in distinct naive populations. We examined the evolution and spread of RHDV in Australia and New Zealand through a genome-wide evolutionary analysis, including data from 28 newly sequenced RHDV field isolates. Following the release of the Australian inoculum strain into New Zealand, no subsequent mixing of the populations occurred, with viruses from both countries forming distinct groups. Strikingly, the rate of evolution in the capsid gene was higher in the Australian viruses than in those from New Zealand, most likely due to the presence of transient deleterious mutations in the former. However, estimates of both substitution rates and population dynamics were strongly sample dependent, such that small changes in sample composition had an important impact on evolutionary parameters. Phylogeographic analysis revealed a clear spatial structure in the Australian RHDV strains, with a major division between those viruses from western and eastern states. Importantly, RHDV sequences from the state where the virus was first released, South Australia, had the greatest diversity and were diffuse throughout both geographic lineages, such that this region was likely a source population for the subsequent spread of the virus across the country. IMPORTANCE: Most studies of viral emergence lack detailed knowledge about which strains were founders for the outbreak or when these events occurred. Hence, the human-mediated introduction of rabbit hemorrhagic disease virus (RHDV) into Australia and New Zealand from known starting stocks provides a unique opportunity to understand viral evolution and emergence. Within Australia, we revealed a major phylogenetic division between viruses sampled from the east and west of the country, with both regions likely seeded by viruses from South Australia. Despite their common origins, marked differences in evolutionary rates were observed between the Australian and New Zealand RHDV, which led to conflicting conclusions about population growth rates. An analysis of mutational patterns suggested that evolutionary rates have been elevated in the Australian viruses, at least in part due to the presence of low-fitness (deleterious) variants that have yet to be selectively purged.

Identification and characterization of antibody-binding epitopes on the norovirus GII.3 capsid.

Genotype II.3 (GII.3) noroviruses are a major cause of sporadic gastroenteritis, particularly in children. The greater incidence of GII.3 noroviruses in the pediatric population compared to the adult demographic suggests development of herd immunity to this genotype, possibly as a consequence of limited evolution of immune epitopes. This study aimed to identify and characterize immune epitopes on the GII.3 capsid protein and to determine the level of immune cross-reactivity within the genotype. A panel of seven GII.3 virus-like particles (VLPs), representing norovirus strains isolated during 1975 to 2008, was tested by enzyme-linked immunosorbent assay (ELISA) for reactivity with human sera and a rabbit anti-GII.3 strain-specific polyclonal serum generated against the 2008 GII.3 VLP. Immunoprecipitation of protease-digested GII.3 VLPs and sequencing of bound peptides via mass spectrometry were used to locate epitopes on the capsid. Two epitopes were investigated further using Mimotopes technology. Serum binding studies demonstrated complete intragenotype GII.3 cross-reactivity using both human and rabbit serum. Six immunoreactive regions containing epitopes were located on the GII.3 capsid protein, two within each capsid domain. Epitopes in the S and P1 domains were highly conserved within GII.3 noroviruses. P2 domain epitopes were variable and contained evolutionarily important residues and histo-blood group antigen (HBGA) binding residues. In conclusion, anti-GII.3 antibody-binding epitopes are highly cross-reactive and mostly conserved within GII.3 strains. This may account for the limited GII.3 prevalence in adults and suggests that a GII.3 strain may be a valuable inclusion in a multivalent pediatric targeted VLP vaccine. Exploration of norovirus immune epitopes is vital for effective vaccine design. IMPORTANCE This study represents an important contribution to the understanding of norovirus immunology in a pediatric genotype. The high cross-reactivity and conservation of GII.3 epitopes suggest development of herd immunity against GII.3 and indicate that a GII.3 strain would be a valuable inclusion in a pediatric targeted multivalent vaccine. Immunological understanding of pediatric norovirus strains is important since norovirus vaccines will likely target high-risk groups such as the pediatric population.

The importance of intergenic recombination in norovirus GII.3 evolution.

Norovirus genotype II.3 (GII.3) strains are a major cause of sporadic gastroenteritis. Intergenic recombination between the capsid and RNA-dependent RNA polymerase (RdRp) genes is common and results in the acquisition of an alternative RdRp genotype. This study aimed to explore the evolution of the GII.3 capsid gene, focusing on the influence of intergenic recombination. The capsid genes from six GII.3 norovirus strains, collected from Australian children between 2001 and 2010, were sequenced and aligned with 66 GII.3 capsid sequences from GenBank, spanning 1975 to 2010. The GII.3 capsid gene evolved at a rate of 4.16 × 10(-3) to 6.97 × 10(-3) nucleotide substitutions/site/year from 1975 to 2010 and clustered into five temporally sequential lineages. Clustering of the GII.3 capsid gene sequences was associated with intergenic recombination and switches between RdRp genotypes GII.3, GII.a, GII.b, GII.12, and an undefined ancestral RdRp. Comparison of the substitution rate of the GII.3 and GII.b RdRps suggested that RdRp switching allows a higher evolutionary rate, leading to increased genetic diversity and adaptability. Alignment of GII.3 capsid sequences revealed 36 lineage-specific conserved amino acid substitutions, four of which were under positive selection. Many conserved substitutions were within predicted antibody binding regions and close to host attachment factor binding sites. In conclusion, evolution of GII.3 noroviruses was primarily driven by intergenic recombination. The acquisition of new RdRps may lead to a faster mutation rate and increased genetic diversity, improving overall GII.3 fitness.

Genomic and phylogenetic features of the Picobirnaviridae suggest microbial rather than animal hosts.

The RNA virus family Picobirnaviridae has traditionally been associated with the gastrointestinal systems of terrestrial mammals and birds, with the majority of viruses detected in animal stool samples. Metatranscriptomic studies of vertebrates, invertebrates, microbial communities, and environmental samples have resulted in an enormous expansion of the genomic and phylogenetic diversity of this family. Yet picobirnaviruses remain poorly classified, with only one genus and three species formally ratified by the International Committee of Virus Taxonomy. Additionally, an inability to culture picobirnaviruses in a laboratory setting or isolate them in animal tissue samples, combined with the presence of bacterial genetic motifs in their genomes, suggests that these viruses may represent RNA bacteriophage rather than being associated with animal infection. Utilising a data set of 2,286 picobirnaviruses sourced from mammals, birds, reptiles, fish, invertebrates, microbial communities, and environmental samples, we identified seven consistent phylogenetic clusters likely representing Picobirnavirus genera that we tentatively name 'Alpha-', 'Beta-', 'Gamma-', 'Delta-', 'Epsilon-', 'Zeta-', and 'Etapicobirnavirus'. A statistical analysis of topological congruence between virus-host phylogenies revealed more frequent cross-species transmission than any other RNA virus family. In addition, bacterial ribosomal binding site motifs were more enriched in Picobirnavirus genomes than in the two groups of established RNA bacteriophage-the Leviviricetes and Cystoviridae. Overall, our findings support the hypothesis that the Picobirnaviridae have bacterial hosts and provide a lower-level taxonomic classification for this highly diverse and ubiquitous family of RNA viruses.

The diverse liver viromes of Australian geckos and skinks are dominated by hepaciviruses and picornaviruses and reflect host taxonomy and habitat.

Lizards have diverse ecologies and evolutionary histories, and represent a promising group to explore how hosts shape virome structure and virus evolution. Yet, little is known about the viromes of these animals. In Australia, squamates (lizards and snakes) comprise the most diverse order of vertebrates, and Australia hosts the highest diversity of lizards globally, with the greatest breadth of habitat use. We used meta-transcriptomic sequencing to determine the virome of nine co-distributed, tropical lizard species from three taxonomic families in Australia and analyzed these data to identify host traits associated with viral abundance and diversity. We show that lizards carry a large diversity of viruses, identifying more than thirty novel, highly divergent vertebrate-associated viruses. These viruses were from nine viral families, including several that contain well known pathogens, such as the Flaviviridae, Picornaviridae, Bornaviridae, Iridoviridae, and Rhabdoviridae. Members of the Flaviviridae were particularly abundant across species sampled here, largely belonging to the genus Hepacivirus: fourteen novel hepaciviruses were identified, broadening the known diversity of this group and better defining its evolution by uncovering new reptilian clades. The evolutionary histories of the viruses studied here frequently aligned with the biogeographic and phylogenetic histories of the hosts, indicating that exogenous viruses may help infer host evolutionary history if sampling is strategic and sampling density high enough. Notably, analysis of alpha and beta diversity revealed that virome composition and richness in the animals sampled here was shaped by host taxonomy and habitat. In sum, we identified a diverse range of reptile viruses that broadly contributes to our understanding of virus-host ecology and evolution.

Australian terrestrial environments harbour extensive RNA virus diversity.

Australia is home to a diverse range of unique native fauna and flora. To address whether Australian ecosystems also harbour unique viruses, we performed meta-transcriptomic sequencing of 16 farmland and sediment samples taken from the east and west coasts of Australia. We identified 2460 putatively novel RNA viruses across 18 orders, the vast majority of which belonged to the microbe-associated phylum Lenarviricota. In many orders, such as the Nodamuvirales and Ghabrivirales, the novel viruses identified here comprised entirely new clades. Novel viruses also fell between established genera or families, such as in the Cystoviridae and Picornavirales, while highly divergent lineages were identified in the Sobelivirales and Ghabrivirales. Viral read abundance and alpha diversity were influenced by sampling site, soil type and land use, but not by depth from the surface. In sum, Australian soils and sediments are home to remarkable viral diversity, reflecting the biodiversity of local fauna and flora.

Comment on Shah et al. Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses, 1988-2021. Viruses 2023, 15, 815.

Shah and colleagues [...].

Divergent hepaciviruses, delta-like viruses, and a chu-like virus in Australian marsupial carnivores (dasyurids).

Although Australian marsupials are characterised by unique biology and geographic isolation, little is known about the viruses present in these iconic wildlife species. The Dasyuromorphia are an order of marsupial carnivores found only in Australia that include both the extinct Tasmanian tiger (thylacine) and the highly threatened Tasmanian devil. Several other members of the order are similarly under threat of extinction due to habitat loss, hunting, disease, and competition and predation by introduced species such as feral cats. We utilised publicly available RNA-seq data from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database to document the viral diversity within four Dasyuromorph species. Accordingly, we identified fifteen novel virus sequences from five DNA virus families (Adenoviridae, Anelloviridae, Gammaherpesvirinae, Papillomaviridae, and Polyomaviridae) and three RNA virus taxa: the order Jingchuvirales, the genus Hepacivirus, and the delta-like virus group. Of particular note was the identification of a marsupial-specific clade of delta-like viruses that may indicate an association of deltaviruses with marsupial species. In addition, we identified a highly divergent hepacivirus in a numbat liver transcriptome that falls outside of the larger mammalian clade. We also detect what may be the first Jingchuvirales virus in a mammalian host-a chu-like virus in Tasmanian devils-thereby expanding the host range beyond invertebrates and ectothermic vertebrates. As many of these Dasyuromorphia species are currently being used in translocation efforts to reseed populations across Australia, understanding their virome is of key importance to prevent the spread of viruses to naive populations.

Gammaretroviruses, novel viruses and pathogenic bacteria in Australian bats with neurological signs, pneumonia and skin lesions.

More than 70 bat species are found in mainland Australia. While most studies of bat viromes focus on sampling seemingly healthy individuals, little is known about the viruses and bacteria associated with diseased bats. We performed traditional diagnostic techniques and metatranscriptomic sequencing on tissue samples from 43 Australian bats, comprising three flying fox (Pteropodidae) and two microbat species experiencing a range of disease syndromes, including mass mortality, neurological signs, pneumonia and skin lesions. Of note, we identified the recently discovered Hervey pteropid gammaretrovirus in a bat with lymphoid leukemia, with evidence of replication consistent with an exogenous virus. The possible association of Hervey pteropid gammaretrovirus with lymphoid leukemia clearly merits additional investigation. One novel picornavirus and at least three new astroviruses and bat pegiviruses were also identified in a variety of tissue types, as well as a number of likely bacterial pathogens or opportunistic infections, most notably Pseudomonas aeruginosa.

Utilizing Molecular Epidemiology and Citizen Science for the Surveillance of Lagoviruses in Australia.

Australia has multiple lagoviruses with differing pathogenicity. The circulation of these viruses was traditionally determined through opportunistic sampling events. In the lead up to the nationwide release of RHDVa-K5 (GI.1aP-GI.1a) in 2017, an existing citizen science program, RabbitScan, was augmented to allow members of the public to submit samples collected from dead leporids for lagovirus testing. This study describes the information obtained from the increased number of leporid samples received between 2015 and 2022 and focuses on the recent epidemiological interactions and evolutionary trajectory of circulating lagoviruses in Australia between October 2020 and December 2022. A total of 2771 samples were tested from January 2015 to December 2022, of which 1643 were lagovirus-positive. Notable changes in the distribution of lagovirus variants were observed, predominantly in Western Australia, where RHDV2-4c (GI.4cP-GI.2) was detected again in 2021 after initially being reported to be present in 2018. Interestingly, we found evidence that the deliberately released RHDVa-K5 was able to establish and circulate in wild rabbit populations in WA. Overall, the incorporation of citizen science approaches proved to be a cost-efficient method to increase the sampling area and enable an in-depth analysis of lagovirus distribution, genetic diversity, and interactions. The maintenance of such programs is essential to enable continued investigations of the critical parameters affecting the biocontrol of feral rabbit populations in Australia, as well as to enable the detection of any potential future incursions.

Comparative evolution of GII.3 and GII.4 norovirus over a 31-year period.

Noroviruses are the most common cause of epidemic gastroenteritis. Genotype II.3 is one of the most frequently detected noroviruses associated with sporadic infections. We studied the evolution of the major capsid gene from seven archival GII.3 noroviruses collected during a cross-sectional study at the Children's Hospital in Washington, DC, from 1975 through 1991, together with capsid sequence from 56 strains available in GenBank. Evolutionary analysis concluded that GII.3 viruses evolved at a rate of 4.16 × 10(-3) nucleotide substitutions/site/year (strict clock), which is similar to that described for the more prevalent GII.4 noroviruses. The analysis of the amino acid changes over the 31-year period found that GII.3 viruses evolve at a relatively steady state, maintaining 4% distance, and have a tendency to revert back to previously used residues while preserving the same carbohydrate binding profile. In contrast, GII.4 viruses demonstrate increasing rates of distance over time because of the continued integration of new amino acids and changing HBGA binding patterns. In GII.3 strains, seven sites acting under positive selection were predicted to be surface-exposed residues in the P2 domain, in contrast to GII.4 positively selected sites located primarily in the shell domain. Our study suggests that GII.3 noroviruses caused disease as early as 1975 and that they evolve via a specific pattern, responding to selective pressures induced by the host rather than presenting a nucleotide evolution rate lower than that of GII.4 noroviruses, as previously proposed. Understanding the evolutionary dynamics of prevalent noroviruses is relevant to the development of effective prevention and control strategies.

Faecal virome of the Australian grey-headed flying fox from urban/suburban environments contains novel coronaviruses, retroviruses and sapoviruses.

Bats are important reservoirs for viruses of public health and veterinary concern. Virus studies in Australian bats usually target the families Paramyxoviridae, Coronaviridae and Rhabdoviridae, with little known about their overall virome composition. We used metatranscriptomic sequencing to characterise the faecal virome of grey-headed flying foxes from three colonies in urban/suburban locations from two Australian states. We identified viruses from three mammalian-infecting (Coronaviridae, Caliciviridae, Retroviridae) and one possible mammalian-infecting (Birnaviridae) family. Of particular interest were a novel bat betacoronavirus (subgenus Nobecovirus) and a novel bat sapovirus (Caliciviridae), the first identified in Australian bats, as well as a potentially exogenous retrovirus. The novel betacoronavirus was detected in two sampling locations 1375 km apart and falls in a viral lineage likely with a long association with bats. This study highlights the utility of unbiased sequencing of faecal samples for identifying novel viruses and revealing broad-scale patterns of virus ecology and evolution.

Characterization of norovirus strains in Australian children from 2006 to 2008: prevalence of recombinant strains.

Noroviruses are highly infectious and are the most common cause of gastroenteritis outbreaks. Genotype II.4 strains have been the dominant type identified in adults, however the genotype distribution in children is less clearly defined. This study aimed to detect and genotype norovirus strains infecting children hospitalized with acute gastroenteritis in Melbourne, Australia from 2006 to 2008. Stool samples were collected from 272 children admitted to the Royal Children's Hospital, Melbourne, Australia, with non-rotavirus acute gastroenteritis between April 2006 and December 2008. Using RT-PCR, norovirus was detected in 36% of samples. Strains were genetically characterized via analysis of regions from both the capsid gene and the RNA dependent RNA polymerase (RdRp) gene, to investigate genotype distribution and incidence of recombination. Typing based on the capsid gene (n = 70) detected GII.4 (49%) and GII.3 (46%) as the most predominant genotypes. Strains with a GII.4 capsid were usually assigned a GII.4 RdRp, whereas most strains identified as GII.3 based on capsid typing were assigned a GIIb RdRp (71%). The GII.3/GIIb represent recombinant strains. Sequence analysis of the putative recombination breakpoint was performed for three representative suspected recombinants: GII.3/GIIb (n = 2) and GII.3/GII.12 (n = 1). Recombination analysis confirmed these strains as recombinants and identified putative breakpoints adjacent to the ORF1/ORF2 junction. This study highlights the importance of norovirus infection as a cause of pediatric gastroenteritis. It also reinforces the high circulation of recombinant strains causing disease in children, particularly the GII.3/GIIb strain.

Current challenges to virus discovery by meta-transcriptomics.

Meta-transcriptomic next-generation sequencing has transformed virus discovery, dramatically expanding our knowledge of the known virosphere. Nevertheless, the use of meta-transcriptomics for virus discovery faces important challenges. As this technology becomes more widely adopted, the proportion of viral sequences in public databases with incorrect (e.g. mis-assignment of host) or limited information (e.g. lacking taxonomic classification) is likely to grow, limiting their utility in bioinformatic pipelines for virus discovery. In addition, we currently lack the bioinformatic tools that can accurately identify viruses showing little or no sequence similarity to database viruses or those that represent likely reagent contaminants. Herein, we outline some of the challenges to effective meta-transcriptomic virus discovery as well as their potential solutions.

Frequent intergenotypic recombination between the non-structural and structural genes is a major driver of epidemiological fitness in caliciviruses.

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.