Metagenomic analysis of individual mosquito viromes reveals the geographical patterns and drivers of viral diversity

Mosquito transmitted viruses are responsible for an increasing burden of human disease. Despite this, little is known about the diversity and ecology of viruses within individual mosquito hosts. Here, using a meta-transcriptomic approach, we determined the viromes of 2,438 individual mosquitoes (81 species), spanning ~4,000 km along latitudes and longitudes in China. From these data we identified 393 viral species associated with mosquitoes, including 7 (putative) species of arthropod-borne viruses (that is, arboviruses). We identified potential mosquito species and geographic hotspots of viral diversity and arbovirus occurrence, and demonstrated that the composition of individual mosquito viromes was strongly associated with host phylogeny. Our data revealed a large number of viruses shared among mosquito species or genera, enhancing our understanding of the host specificity of insect-associated viruses. We also detected multiple virus species that were widespread throughout the country, perhaps reflecting long-distance mosquito dispersal. Together, these results greatly expand the known mosquito virome, linked viral diversity at the scale of individual insects to that at a country-wide scale, and offered unique insights into the biogeography and diversity of viruses in insect vectors.

Metagenomic analysis of individual mosquitos reveals the ecology of insect viruses

Mosquito transmitted viruses are responsible for an increasing burden of human disease. Despite this, little is known about the diversity and ecology of viruses within individual mosquito hosts. Using a meta-transcriptomic approach, we analysed the virome of 2,438 individual mosquitos (79 species), spanning ~4000 km along latitudes and longitudes in China. From these data we identified 393 core viral species associated with mosquitos, including seven (putative) arbovirus species. We identified potential species and geographic hotspots of viral richness and arbovirus occurrence, and demonstrated that host phylogeny had a strong impact on the composition of individual mosquito viromes. Our data revealed a large number of viruses shared among mosquito species or genera, expanding our knowledge of host specificity of insect-associated viruses. We also detected multiple virus species that were widespread throughout the country, possibly facilitated by long-distance mosquito migrations. Together, our results greatly expand the known mosquito virome, linked the viral diversity at the scale of individual insects to that at a country-wide scale, and offered unique insights into the ecology of viruses of insect vectors.

RAPIDprep: A Simple, Fast Protocol for RNA Metagenomic Sequencing of Clinical Samples

Emerging infectious disease threats require rapid response tools to inform diagnostics, treatment, and outbreak control. RNA-based metagenomics offers this; however, most approaches are time-consuming and laborious. Here, we present a simple and fast protocol, the RAPIDprep assay, with the aim of providing a cause-agnostic laboratory diagnosis of infection within 24 h of sample collection by sequencing ribosomal RNA-depleted total RNA. The method is based on the synthesis and amplification of double-stranded cDNA followed by short-read sequencing, with minimal handling and clean-up steps to improve processing time. The approach was optimized and applied to a range of clinical respiratory samples to demonstrate diagnostic and quantitative performance. Our results showed robust depletion of both human and microbial rRNA, and library amplification across different sample types, qualities, and extraction kits using a single workflow without input nucleic-acid quantification or quality assessment. Furthermore, we demonstrated the genomic yield of both known and undiagnosed pathogens with complete genomes recovered in most cases to inform molecular epidemiological investigations and vaccine design. The RAPIDprep assay is a simple and effective tool, and representative of an important shift toward the integration of modern genomic techniques with infectious disease investigations.

New Detection of Locally Acquired Japanese Encephalitis Virus Using Clinical Metagenomics, New South Wales, Australia

In the context of an emerging Japanese encephalitis outbreak within Australia, we describe a novel locally acquired case in New South Wales. A man in his 70s had rapidly progressive, fatal meningoencephalitis, diagnosed as caused by Japanese encephalitis virus by RNA-based metagenomic next-generation sequencing performed on postmortem brain tissue.

Potent monoclonal antibody-mediated neutralization of a divergent Hendra virus variant

Hendra virus (HeV) and Nipah virus (NiV) are deadly zoonotic Henipaviruses (HNVs) responsible for recurrent outbreaks in humans and domestic species of highly fatal (50 to 95%) disease. A HeV variant (HeV-g2) of unprecedented genetic divergence has been identified in two fatally diseased horses, and in two flying fox species in regions of Australia not previously considered at risk for HeV spillover. Given the HeV-g2 divergence from HeV while retaining equivalent pathogenicity and spillover potential, understanding receptor usage and antigenic properties is urgently required to guide One Health biosecurity. Here, we show that the HeV-g2 G glycoprotein shares a conserved receptor tropism with prototypic HeV and that a panel of monoclonal antibodies recognizing the G and F glycoproteins potently neutralizes HeV-g2– and HeV G/F–mediated entry into cells. We determined a crystal structure of the Fab fragment of the hAH1.3 antibody bound to the HeV G head domain, revealing an antigenic site associated with potent cross-neutralization of both HeV-g2 and HeV. Structure-guided formulation of a tetravalent monoclonal antibody (mAb) mixture, targeting four distinct G head antigenic sites, results in potent neutralization of HeV and HeV-g2 and delineates a path forward for implementing multivalent mAb combinations for postexposure treatment of HNV infections.

Novel Hendra Virus Variant Circulating in Black Flying Foxes and Grey-Headed Flying Foxes, Australia

A novel Hendra virus variant, genotype 2, was recently discovered in a horse that died after acute illness and in Pteropus flying fox tissues in Australia. We detected the variant in flying fox urine, the pathway relevant for spillover, supporting an expanded geographic range of Hendra virus risk to horses and humans.

Novel Hendra Virus Variant Detected by Sentinel Surveillance of Horses in Australia

We identified and isolated a novel Hendra virus (HeV) variant not detected by routine testing from a horse in Queensland, Australia, that died from acute illness with signs consistent with HeV infection. Using whole-genome sequencing and phylogenetic analysis, we determined the variant had ≈83% nt identity with prototypic HeV. In silico and in vitro comparisons of the receptor-binding protein with prototypic HeV support that the human monoclonal antibody m102.4 used for postexposure prophylaxis and current equine vaccine will be effective against this variant. An updated quantitative PCR developed for routine surveillance resulted in subsequent case detection. Genetic sequence consistency with virus detected in grey-headed flying foxes suggests the variant circulates at least among this species. Studies are needed to determine infection kinetics, pathogenicity, reservoir-species associations, viral-host coevolution, and spillover dynamics for this virus. Surveillance and biosecurity practices should be updated to acknowledge HeV spillover risk across all regions frequented by flying foxes.

Metagenomic Identification of Viral Sequences in Laboratory Reagents

Metagenomic next-generation sequencing has transformed the discovery and diagnosis of infectious disease, with the power to characterise the complete 'infectome' (bacteria, viruses, fungi, parasites) of an individual host organism. However, the identification of novel pathogens has been complicated by widespread microbial contamination in commonly used laboratory reagents. Using total RNA sequencing ("metatranscriptomics") we documented the presence of contaminant viral sequences in multiple 'blank' negative control sequencing libraries that comprise a sterile water and reagent mix. Accordingly, we identified 14 viral sequences in 7 negative control sequencing libraries. As in previous studies, several circular replication-associated protein encoding (CRESS) DNA virus-like sequences were recovered in the blank control libraries, as well as contaminating sequences from the Totiviridae, Tombusviridae and Lentiviridae families of RNA virus. These data suggest that viral contamination of common laboratory reagents is likely commonplace and can comprise a wide variety of viruses.

The spatial-temporal dynamics of respiratory syncytial virus infections across the east-west coasts of Australia during 2016-17

Respiratory syncytial virus (RSV) is an important human respiratory pathogen. In temperate regions, a distinct seasonality is observed, where peaks of infections typically occur in early winter, often preceding the annual influenza season. Infections are associated with high rates of morbidity and mortality and in some populations exceed that of influenza. Two subtypes, RSV-A and RSV-B, have been described, and molecular epidemiological studies have shown that both viruses mostly co-circulate. This trend also appears to be the case for Australia; however, previous genomic studies have been limited to cases from one Eastern state—New South Wales. As such, the broader spatial patterns and viral traffic networks across the continent are not known. Here, we conducted a whole-genome study of RSV comparing strains across eastern and Western Australia during the period January 2016 to June 2017. In total, 96 new RSV genomes were sequenced, compiled with previously generated data, and examined using a phylodynamic approach. This analysis revealed that both RSV-A and RSV-B strains were circulating, and each subtype was dominated by a single genotype, RSV-A ON1-like and RSV-B BA10-like viruses. Some geographical clustering was evident in strains from both states with multiple distinct sub-lineages observed and relatively low mixing across jurisdictions, suggesting that endemic transmission was likely seeded from imported, unsampled locations. Overall, the RSV phylogenies reflected a complex pattern of interactions across multiple epidemiological scales from fluid virus traffic across global and regional networks to fine-scale local transmission events.

Cellular activation, differentiation and proliferation influence the dynamics of genetically-intact proviruses over time

Human immunodeficiency virus (HIV) persists in cells despite antiretroviral therapy; however, the influence of cellular mechanisms such as activation, differentiation, and proliferation upon the distribution of proviruses over time is unclear. To address this, we used full-length sequencing to examine proviruses within memory CD4⁺ T-cell subsets longitudinally in 8 participants. Over time, the odds of identifying a provirus increased in effector and decreased in transitional memory cells. In all subsets, more activated (HLA-DR–expressing) cells contained a higher frequency of intact provirus, as did more differentiated cells such as transitional and effector memory subsets. The proportion of genetically identical proviruses increased over time, indicating that cellular proliferation was maintaining the persistent reservoir; however, the number of genetically identical proviral clusters in each subset was stable. As such, key biological processes of activation, differentiation, and proliferation influence the dynamics of the HIV reservoir and must be considered during the development of any immune intervention.

Genomic Insights Into the Pathogenicity of a Novel Biofilm-Forming Enterococcus sp. Bacteria (Enterococcus lacertideformus) Identified in Reptiles

Whole genome analysis of a novel species of enterococci, Enterococcus lacertideformus, causing multi-systemic and invariably fatal disease in critically endangered Christmas Island reptiles was undertaken to determine the genetic elements and potential mechanisms conferring its pathogenic nature, biofilm-forming capabilities, immune recognition avoidance, and inability to grow in vitro. Comparative genomic analyses with related and clinically significant enterococci were further undertaken to infer the evolutionary history of the bacterium and identify genes both novel and absent. The genome had a G + C content of 35.1%, consisted of a circular chromosome, no plasmids, and was 2,419,934 bp in length (2,321 genes, 47 tRNAs, and 13 rRNAs). Multi-locus sequence typing (MLST), and single nucleotide polymorphism (SNP) analysis of multiple E. lacertideformus samples revealed they were effectively indistinguishable from one another and highly clonal. E. lacertideformus was found to be located within the Enterococcus faecium species clade and was closely related to Enterococcus villorum F1129D based on 16S rDNA and MLST house-keeping gene analysis. Antimicrobial resistance (DfreE, EfrB, tetM, bcrRABD, and sat4) and virulence genes (Fss3 and ClpP), and genes conferring tolerance to metals and biocides (n = 9) were identified. The detection of relatively few genes encoding antimicrobial resistance and virulence indicates that this bacterium may have had no exposure to recently developed and clinically significant antibiotics. Genes potentially imparting beneficial functional properties were identified, including prophages, insertion elements, integrative conjugative elements, and genomic islands. Functional CRISPR-Cas arrays, and a defective prophage region were identified in the genome. The study also revealed many genomic loci unique to E. lacertideformus which contained genes enriched in cell wall/membrane/envelop biogenesis, and carbohydrate metabolism and transport functionality. This finding and the detection of putative enterococcal biofilm determinants (EfaAfs, srtC, and scm) may underpin the novel biofilm phenotype observed for this bacterium. Comparative analysis of E. lacertideformus with phylogenetically related and clinically significant enterococci (E. villorum F1129D, Enterococcus hirae R17, E. faecium AUS0085, and Enterococcus faecalis OG1RF) revealed an absence of genes (n = 54) in E. lacertideformus, that encode metabolic functionality, which potentially hinders nutrient acquisition and/or utilization by the bacterium and precludes growth in vitro. These data provide genetic insights into the previously determined phenotype and pathogenic nature of the bacterium.

Identification of SARS-CoV-2 Nucleocapsid and Spike T-Cell Epitopes for Assessing T-Cell Immunity

Developing optimal T-cell response assays to severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) is critical for measuring the duration of immunity to this disease and assessing the efficacy of vaccine candidates. These assays need to target conserved regions of SARS-CoV-2 global variants and avoid cross-reactivity to seasonal human coronaviruses. To contribute to this effort, we employed an in silico immunoinformatics analysis pipeline to identify immunogenic peptides resulting from conserved and highly networked regions with topological importance from the SARS-CoV-2 nucleocapsid and spike proteins. A total of 57 highly networked T-cell epitopes that are conserved across geographic viral variants were identified from these viral proteins, with a binding potential to diverse HLA alleles and 80 to 100% global population coverage. Importantly, 18 of these T-cell epitope derived peptides had limited homology to seasonal human coronaviruses making them promising candidates for SARS-CoV-2-specific T-cell immunity assays. Moreover, two of the NC-derived peptides elicited effector/polyfunctional responses of CD8+ T cells derived from SARS-CoV-2 convalescent patients.IMPORTANCE The development of specific and validated immunologic tools is critical for understanding the level and duration of the cellular response induced by SARS-CoV-2 infection and/or vaccines against this novel coronavirus disease. To contribute to this effort, we employed an immunoinformatics analysis pipeline to define 57 SARS-CoV-2 immunogenic peptides within topologically important regions of the nucleocapsid (NC) and spike (S) proteins that will be effective for detecting cellular immune responses in 80 to 100% of the global population. Our immunoinformatics analysis revealed that 18 of these peptides had limited homology to circulating seasonal human coronaviruses and therefore are promising candidates for distinguishing SARS-CoV-2-specific immune responses from pre-existing coronavirus immunity. Importantly, CD8+ T cells derived from SARS-CoV-2 survivors exhibited polyfunctional effector responses to two novel NC-derived peptides identified as HLA-binders. These studies provide a proof of concept that our immunoinformatics analysis pipeline identifies novel immunogens which can elicit polyfunctional SARS-CoV-2-specific T-cell responses.

Genome sequencing and its use in public health responses to COVID-19

Human history has been shaped by the heavy burden of infectious disease pandemics. Yet, despite the bitter lessons learned from history, even those in living memory such as the 1918 influenza pandemic and HIV/AIDS epidemic, COVID-19 stands unique in the sudden, immense health and economic impacts to the global human population. While the costs have been great, and the long-term consequences are still being revealed, the urgent need for action has also brought forward rapid developments and innovations to combat COVID-19 and better prepare us for future infectious disease outbreaks. One such area has been the widespread adoption of whole genome sequencing to inform public health responses. Genome sequencing during the COVID-19 pandemic has become key to tracking the spread of SARS-CoV-2 at all scales, to such a degree that terms such as genomics, mutations, variants and clusters are now common vernacular to politicians, health officials and the general public. This article provides a commentary on the genesis and evolution of SARS-CoV-2 genome sequencing, and its critical on-going role in the public health response to the COVID-19 pandemic.

Analytical validity of nanopore sequencing for rapid SARS-CoV-2 genome analysis

Viral whole-genome sequencing (WGS) provides critical insight into the transmission and evolution of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Long-read sequencing devices from Oxford Nanopore Technologies (ONT) promise significant improvements in turnaround time, portability and cost, compared to established short-read sequencing platforms for viral WGS (e.g., Illumina). However, adoption of ONT sequencing for SARS-CoV-2 surveillance has been limited due to common concerns around sequencing accuracy. To address this, here we perform viral WGS with ONT and Illumina platforms on 157 matched SARS-CoV-2-positive patient specimens and synthetic RNA controls, enabling rigorous evaluation of analytical performance. We report that, despite the elevated error rates observed in ONT sequencing reads, highly accurate consensus-level sequence determination was achieved, with single nucleotide variants (SNVs) detected at >99% sensitivity and >99% precision above a minimum ~60-fold coverage depth, thereby ensuring suitability for SARS-CoV-2 genome analysis. ONT sequencing also identified a surprising diversity of structural variation within SARS-CoV-2 specimens that were supported by evidence from short-read sequencing on matched samples. However, ONT sequencing failed to accurately detect short indels and variants at low read-count frequencies. This systematic evaluation of analytical performance for SARS-CoV-2 WGS will facilitate widespread adoption of ONT sequencing within local, national and international COVID-19 public health initiatives.

Genetic diversity of SARS-CoV-2 and clinical, epidemiological characteristics of COVID-19 patients in Hanoi, Vietnam

A second cluster of COVID-19 cases imported from Europe occured in Vietnam from early March 2020. We describe 44 SARS-CoV-2 RT-PCR positive patients (cycle threshold value <30) admitted to the National Hospital for Tropical Diseases in Hanoi between March 6 and April 15 2020. Whole SARS-CoV-2 genomes from these patients were sequenced using Illumina Miseq and analysed for common genetic variants and relationships to local and globally circulating strains. Results showed that 32 cases were Vietnamese with a median age of 37 years (range 15-74 years), and 23 were male. Most cases were acquired outside Vietnam, mainly from the UK (n = 15), other European countries (n = 14), Russia (n = 6) and countries in Asia (n = 3). No cases had travelled from China. Forty-one cases had symptoms at admission, typically dry cough (n = 36), fever (n = 20), sore throat (n = 14) and diarrhoea (n = 12). Hospitalisation was long with a median of 25 days, most commonly from 20-29 days. All SARS-CoV-2 genomes were similar (92-100% sequence homology) to the reference sequence Wuhan_1 (NC_045512), and 32 strains belonged to the B.1.1 lineage. The three most common variants were linked, and included C3037T, C14408T (nsp12: P323L) and A23403G (S: D614G) mutations. This group of mutations often accompanied variant C241T (39/44 genomes) or GGG 28881..28883 AAC (33/44 genomes). The prevalence of the former reflected probable European origin of viruses, and the transition D614G was dominant in Vietnam. New variants were identified; however, none could be associated with disease severity.

Meta-Transcriptomic Discovery of a Divergent Circovirus and a Chaphamaparvovirus in Captive Reptiles with Proliferative Respiratory Syndrome

Viral pathogens are being increasingly described in association with mass morbidity and mortality events in reptiles. However, our knowledge of reptile viruses remains limited. Herein, we describe the meta-transcriptomic investigation of a mass morbidity and mortality event in a colony of central bearded dragons (Pogona vitticeps) in 2014. Severe, extensive proliferation of the respiratory epithelium was consistently found in affected dragons. Similar proliferative lung lesions were identified in bearded dragons from the same colony in 2020 in association with increased intermittent mortality. Total RNA sequencing identified two divergent DNA viruses: a reptile-infecting circovirus, denoted bearded dragon circovirus (BDCV), and the first exogeneous reptilian chaphamaparvovirus—bearded dragon chaphamaparvovirus (BDchPV). Phylogenetic analysis revealed that BDCV was most closely related to bat-associated circoviruses, exhibiting 70% amino acid sequence identity in the Replicase (Rep) protein. In contrast, in the nonstructural (NS) protein, the newly discovered BDchPV showed approximately 31%–35% identity to parvoviruses obtained from tilapia fish and crocodiles in China. Subsequent specific PCR assays revealed BDCV and BDchPV in both diseased and apparently normal captive reptiles, although only BDCV was found in those animals with proliferative pulmonary lesions and respiratory disease. This study expands our understanding of viral diversity in captive reptiles.

Meta-transcriptomic identification of Trypanosoma spp. in native wildlife species from Australia

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.

Novel hepaci- and pegi-like viruses in native Australian wildlife and non-human primates

The Flaviviridae family of positive-sense RNA viruses contains important pathogens of humans and other animals, including Zika virus, dengue virus, and hepatitis C virus. The Flaviviridae are currently divided into four genera-Hepacivirus, Pegivirus, Pestivirus, and Flavivirus-each with a diverse host range. Members of the genus Hepacivirus are associated with an array of animal species, including humans, non-human primates, other mammalian species, as well as birds and fish, while the closely related pegiviruses have been identified in a variety of mammalian taxa, also including humans. Using a combination of total RNA and whole-genome sequencing we identified four novel hepaci-like viruses and one novel variant of a known hepacivirus in five species of Australian wildlife. The hosts infected comprised native Australian marsupials and birds, as well as a native gecko (Gehyra lauta). From these data we identified a distinct marsupial clade of hepaci-like viruses that also included an engorged Ixodes holocyclus tick collected while feeding on Australian long-nosed bandicoots (Perameles nasuta). Distinct lineages of hepaci-like viruses associated with geckos and birds were also identified. By mining the SRA database we similarly identified three new hepaci-like viruses from avian and primate hosts, as well as two novel pegi-like viruses associated with primates. The phylogenetic history of the hepaci- and pegi-like viruses as a whole, combined with co-phylogenetic analysis, provided support for virus-host co-divergence over the course of vertebrate evolution, although with frequent cross-species virus transmission. Overall, our work highlights the diversity of the Hepacivirus and Pegivirus genera as well as the uncertain phylogenetic distinction between.

Revealing COVID-19 transmission in Australia by SARS-CoV-2 genome sequencing and agent-based modeling

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, China^1,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.

A Divergent Articulavirus in an Australian Gecko Identified Using Meta-Transcriptomics and Protein Structure Comparisons

The discovery of highly divergent RNA viruses is compromised by their limited sequence similarity to known viruses. Evolutionary information obtained from protein structural modelling offers a powerful approach to detect distantly related viruses based on the conservation of tertiary structures in key proteins such as the RNA-dependent RNA polymerase (RdRp). We utilised a template-based approach for protein structure prediction from amino acid sequences to identify distant evolutionary relationships among viruses detected in meta-transcriptomic sequencing data from Australian wildlife. The best predicted protein structural model was compared with the results of similarity searches against protein databases. Using this combination of meta-transcriptomics and protein structure prediction we identified the RdRp (PB1) gene segment of a divergent negative-sense RNA virus, denoted Lauta virus (LTAV), in a native Australian gecko (Gehyra lauta). The presence of this virus was confirmed by PCR and Sanger sequencing. Phylogenetic analysis revealed that Lauta virus likely represents a newly described genus within the family Amnoonviridae, order Articulavirales, that is most closely related to the fish virus Tilapia tilapinevirus (TiLV). These findings provide important insights into the evolution of negative-sense RNA viruses and structural conservation of the viral replicase among members of the order Articulavirales.

Opportunistic sampling of wild native and invasive birds reveals a rich diversity of adenoviruses in Australia

Little is known about the diversity of adenoviruses in wild birds and how they have evolved and are maintained in complex ecosystems. In this study, 409 samples were collected from woodland birds caught for banding (droppings), birds submitted to a wildlife hospital (droppings and tissues), silver gulls (droppings or tissues), and feral pigeons (Columbia livia; oral, cloacal swabs, or tissues) from the Greater Sydney area in NSW, Australia. Additional samples were from native pigeons and doves (swabs) presented to the Healesville Sanctuary, VIC, Australia. Samples were screened for adenovirus DNA using degenerate primers and polymerase chain reaction. Adenovirus sequences were detected in eighty-three samples representing thirty-five novel amino acid sequences. Fourteen novel sequences were atadenoviruses, seven were aviadenoviruses, twelve were siadenoviruses, and one was a mastadenovirus. Sequences from passerine birds were predominately found to form a single lineage within the atadenoviruses, a second lineage in the siadenoviruses, and a third smaller aviadenovirus lineage. These viruses appeared to have co-evolved with a diverse group of woodland birds that share similar habitat. Evidence for host/virus co-evolution in some viruses and a wide host range in others was observed. A high prevalence of adenovirus infection was found in rainbow lorikeets (Trichoglossus haematodus), galahs (Eolophus roseicapilla), and sulphur-crested cockatoos (Cacatua galerita). Sequences were either identical to or mapped to already established lineages in the Aviadenovirus, Siadenovirus, and Atadenovirus genera, suggesting a possible origin of the psittacine adenoviruses in ancestral Australian psittacine birds. The sequences of passerine and psittacine origin provided insight into diversity and structure of the Atadenovirus genus and demonstrated for the first-time viruses of passerine origin in the Aviadenovirus genus. Four unrelated adenovirus sequences were found in silver gull samples (Chroicocephalus novaehollandiae), including one of pigeon origin, suggesting environmental virus exposure. Three pigeon adenovirus types were detected in feral pigeons and infection prevalence was high. Evidence for host switching between invasive species and native species and native species and invasive species was documented. A variant of a murine adenovirus was detected in kidney tissue from two bird species suggesting mouse to bird transmission.

An emergent clade of SARS-CoV-2 linked to returned travellers from Iran

The SARS-CoV-2 epidemic has rapidly spread outside China with major outbreaks occurring in Italy, South Korea, and Iran. Phylogenetic analyses of whole-genome sequencing data identified a distinct SARS-CoV-2 clade linked to travellers returning from Iran to Australia and New Zealand. This study highlights potential viral diversity driving the epidemic in Iran, and underscores the power of rapid genome sequencing and public data sharing to improve the detection and management of emerging infectious diseases.

Novel RNA viruses associated with Plasmodium vivax in human malaria and Leucocytozoon parasites in avian disease

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.

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.

New insights into Sauropsid Papillomaviridae evolution and epizootiology: discovery of two novel papillomaviruses in native and invasive Island geckos

Papillomaviruses cause persistent infections in skin and mucosal membranes and, in at least one species, are also be able to infect a tissue of mesenchymal origin. Infections may either be subclinical or induce proliferative lesions. Of the known papillomaviruses, the majority that have been characterized are from humans and other mammals. Currently, only fifteen complete bird and reptile papillomavirus genomes have been described, and they have been found in birds (n = 11), turtles (n = 2), and snakes (n = 2). Using next-generation sequencing technologies and virus-specific PCR, we have identified two novel papillomavirus genomes, Hemidactylus frenatus Papillomavirus 1 and 2 (HfrePV1, HfrePV2), in the widely distributed and highly invasive Asian house gecko (H.frenatus) and mute gecko (Gehyra mutilata) on Christmas Island and Cocos (Keeling) Islands. HfrePV1 was also detected in critically endangered Lister’s geckos (Lepidodactylus listeri) in their captive breeding colony on Christmas Island. Tissue-containing virus included epidermis, oral mucosa, and liver (HfrePV1) and epidermis, liver, and colon (HfrePV2). Concurrent infections were found in both H.frenatus and G.mutilata. Invasive mourning geckos (Lepidodactylus lugubris) (n = 4), Sri Lankan house geckos (Hemidactylus parvimaculatus) (n = 3), flat-tailed house geckos (Hemidactylus platyurus) (n = 4) from the Cocos Islands, and blue-tailed skinks (Cryptoblepharus egeriae) (n = 10) from Christmas Island were also screened but were not found to be infected. The novel HfrePV1 and HfrePV2 genomes were 7,378 bp and 7,380 bp in length, respectively, and each contained the early (E1, E2, and E7), and late (L1 and L2) open-reading frames. Phylogenetic analysis of the concatenated E1, E2, and L1 proteins from both papillomaviruses revealed that they clustered with, but were basal to, the Sauropsida clade containing bird and reptile viruses. This study sheds light on the evolution of papillomaviruses and the distribution of pathogens in a highly invasive species impacting endangered populations of geckos.

Human pegivirus in brain tissue of a patient with encephalitis

We describe a case of meningoencephalitis in which meta-transcriptomic (RNA) sequencing detected human pegivirus (HPgV) in brain tissue, cerebrospinal fluid, and serum in the absence of other pathogens. This is the first detection of HPgV in antemortem brain tissue, although it is uncertain whether HPgV is responsible for the observed encephalitis.

Metatranscriptomics as a tool to identify fungal species and subspecies in mixed communities - a proof of concept under laboratory conditions

High-throughput sequencing (HTS) enables the generation of large amounts of genome sequence data at a reasonable cost. Organisms in mixed microbial communities can now be sequenced and identified in a culture-independent way, usually using amplicon sequencing of a DNA barcode. Bulk RNA-seq (metatranscriptomics) has several advantages over DNA-based amplicon sequencing: it is less susceptible to amplification biases, it captures only living organisms, and it enables a larger set of genes to be used for taxonomic identification. Using a model mock community comprising 17 fungal isolates, we evaluated whether metatranscriptomics can accurately identify fungal species and subspecies in mixed communities. Overall, 72.9% of the RNA transcripts were classified, from which the vast majority (99.5%) were correctly identified at the species level. Of the 15 species sequenced, 13 were retrieved and identified correctly. We also detected strain-level variation within the Cryptococcus species complexes: 99.3% of transcripts assigned to Cryptococcus were classified as one of the four strains used in the mock community. Laboratory contaminants and/or misclassifications were diverse, but represented only 0.44% of the transcripts. Hence, these results show that it is possible to obtain accurate species- and strain-level fungal identification from metatranscriptome data as long as taxa identified at low abundance are discarded to avoid false-positives derived from contamination or misclassifications. This study highlights both the advantages and current challenges in the application of metatranscriptomics in clinical mycology and ecological studies.

A survey of a mixed species aviary provides new insights into the pathogenicity, diversity, evolution, host range, and distribution of psittacine and passerine adenoviruses

A Bourke's parrot (Neopsephotus bourkii) originating from an aviary in Australia, containing two species of parrots, five species of finch and a species of dove, was presented for necropsy. The Bourke's parrot died from gastritis caused by Macrorhabdus ornithogaster, but also had an interstitial nephritis and ureteritis with adenovirus-like inclusion bodies within collecting duct epithelial cells. The adenovirus causing the lesions was shown to be Psittacine adenovirus-2 (PsAdV-2) using a PCR assay specific for adenoviruses and sequencing of amplicons. A survey of droppings from other birds in the aviary using the same PCR assay with amplicon sequencing found a high prevalence of infection of PsAdV-2 in Bourke's and scarlet-chested parrots (Neophema splendida). PsAdV-2 was also present in droppings from a Namaqua dove (Oena capensis). Gouldian finches (Erythrura gouldiae), red-billed firefinches (Lagonosticta senegala), and red-throated parrot finches (Erythrura psittacea) were shedding Gouldian finch adenovirus-1 (GFAdV-1). Two novel adenoviruses, an atadenovirus and a siadenovirus, were detected in the droppings from long-tailed finches (Poephila acuticauda). Kidney tissue from three of four scarlet-chested parrots submitted for necropsy from a second aviary were also positive for PsAdv-2. These findings and previously reported findings of widespread PsAdv-2 infection in captive orange-bellied parrots (Neophemia chrysogaster) raise the possibility that PsAdV-2 is enzootic in Australian aviculture. This represents the first report of GFAdV-1 in Australia and first identification of infection in finch species other than the Gouldian finch. Identification of two novel adenoviruses in long-tailed finches suggests that other novel adenoviruses are circulating in other finch species. RESEARCH HIGHLIGHTS Psittacine adenovirus-2 was present in high prevalence in two Australian aviaries. Gouldian finch adenovirus-1 (GFAdV-1) was detected in Australia for the first time. The host range of GFAdV-1 host range was expanded to other finch species. Novel atadenovirus and siadenovirus were detected in Estrildid finches.

Novel hepatitis D-like agents in vertebrates and invertebrates

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.

Identification of diverse arthropod associated viruses in native Australian fleas

Fleas are important vectors of zoonotic disease. However, little is known about the natural diversity and abundance of flea viruses, particularly in the absence of disease associations, nor the evolutionary relationships among those viruses found in different parasitic vector species. Herein, we present the first virome scale study of fleas, based on the meta-transcriptomic analysis of 52 fleas collected along the eastern coast of Australia. Our analysis revealed 18 novel RNA viruses belonging to nine viral families with diverse genome organizations, although the majority (72%) possessed single-stranded positive-sense genomes. Notably, a number of the viruses identified belonged to the same phylogenetic groups as those observed in ticks sampled at the same locations, although none were likely associated with mammalian infection. Overall, we identified high levels of genomic diversity and abundance of viruses in the flea species studied, and established that fleas harbor viruses similar to those seen to other vectors.

Meta-transcriptomics reveals a diverse antibiotic resistance gene pool in avian microbiomes

BACKGROUND

Antibiotic resistance is rendering common bacterial infections untreatable. Wildlife can incorporate and disperse antibiotic-resistant bacteria in the environment, such as water systems, which in turn serve as reservoirs of resistance genes for human pathogens. Anthropogenic activity may contribute to the spread of bacterial resistance cycling through natural environments, including through the release of human waste, as sewage treatment only partially removes antibiotic-resistant bacteria. However, empirical data supporting these effects are currently limited. Here we used bulk RNA-sequencing (meta-transcriptomics) to assess the diversity and expression levels of functionally viable resistance genes in the gut microbiome of birds with aquatic habits in diverse locations.

RESULTS

We found antibiotic resistance genes in birds from all localities, from penguins in Antarctica to ducks in a wastewater treatment plant in Australia. Comparative analysis revealed that birds feeding at the wastewater treatment plant carried the greatest resistance gene burden, including genes typically associated with multidrug resistance plasmids as the aac(6)-Ib-cr gene. Differences in resistance gene burden also reflected aspects of bird ecology, taxonomy, and microbial function. Notably, ducks, which feed by dabbling, carried a higher abundance and diversity of resistance genes than turnstones, avocets, and penguins, which usually prey on more pristine waters.

CONCLUSIONS

These transcriptome data suggest that human waste, even if it undergoes treatment, might contribute to the spread of antibiotic resistance genes to the wild. Differences in microbiome functioning across different bird lineages may also play a role in the antibiotic resistance burden carried by wild birds. In summary, we reveal the complex factors explaining the distribution of resistance genes and their exchange routes between humans and wildlife, and show that meta-transcriptomics is a valuable tool to access functional resistance genes in whole microbial communities.

Metagenomic discovery and co-infection of diverse wobbly possum disease viruses and a novel hepacivirus in Australian brushtail possums

BACKGROUND

Australian brushtail possums (Trichosurus vulpecula) are an introduced pest species in New Zealand, but native to Australia where they are protected for biodiversity conservation. Wobbly possum disease (WPD) is a fatal neurological disease of Australian brushtail possums described in New Zealand populations that has been associated with infection by the arterivirus (Arteriviridae) wobbly possum disease virus (WPDV-NZ). Clinically, WPD-infected possums present with chronic meningoencephalitis, choroiditis and multifocal neurological symptoms including ataxia, incoordination, and abnormal gait.

METHODS

We conducted a retrospective investigation to characterise WPD in native Australian brushtail possums, and used a bulk meta-transcriptomic approach (i.e. total RNA-sequencing) to investigate its potential viral aetiology. PCR assays were developed for case diagnosis and full genome recovery in the face of extensive genetic variation.

RESULTS

We identified genetically distinct lineages of arteriviruses from archival tissues of WPD-infected possums in Australia, termed wobbly possum disease virus AU1 and AU2. Phylogenetically, WPDV-AU1 and WPDV-AU2 shared only ~ 70% nucleotide similarity to each other and the WPDV-NZ strain, suggestive of a relatively ancient divergence. Notably, we also identified a novel and divergent hepacivirus (Flaviviridae) - the first in a marsupial - in both WPD-infected and uninfected possums, indicative of virus co-infection.

CONCLUSIONS

We have identified marsupial-specific lineages of arteriviruses in mainland Australia that are genetically distinct from that in New Zealand, in some cases co-infecting animals with a novel hepacivirus. Our study provides new insight into the hidden genetic diversity of arteriviruses, the capacity for virus co-infection, and highlights the utility of meta-transcriptomics for disease investigation in a One Health context.

A Divergent Hepatitis D-Like Agent in Birds

Hepatitis delta virus (HDV) is currently only found in humans and is a satellite virus that depends on hepatitis B virus (HBV) envelope proteins for assembly, release, and entry. Using meta-transcriptomics, we identified the genome of a novel HDV-like agent in ducks. Sequence analysis revealed secondary structures that were shared with HDV, including self-complementarity and ribozyme features. The predicted viral protein shares 32% amino acid similarity to the small delta antigen of HDV and comprises a divergent phylogenetic lineage. The discovery of an avian HDV-like agent has important implications for the understanding of the origins of HDV and sub-viral agents.

Virus-virus interactions and host ecology are associated with RNA virome structure in wild birds

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.

Amplification of Near Full-length HIV-1 Proviruses for Next-Generation Sequencing

The Full-Length Individual Proviral Sequencing (FLIPS) assay is an efficient and high-throughput method designed to amplify and sequence single, near full-length (intact and defective), HIV-1 proviruses. FLIPS allows determination of the genetic composition of integrated HIV-1 within a cell population. Through identifying defects within HIV-1 proviral sequences that arise during reverse transcription, such as large internal deletions, deleterious stop codons/hypermutation, frameshift mutations, and mutations/deletions in cis acting elements required for virion maturation, FLIPS can identify integrated proviruses incapable of replication. The FLIPS assay can be utilized to identify HIV-1 proviruses that lack these defects and are therefore potentially replication-competent. The FLIPS protocol involves: lysis of HIV-1 infected cells, nested PCR of near full-length HIV-1 proviruses (using primers targeted to the HIV-1 5' and 3' LTR), DNA purification and quantification, library preparation for Next-generation Sequencing (NGS), NGS, de novo assembly of proviral contigs, and a simple process of elimination for identifying replication-competent proviruses. FLIPS provides advantages over traditional methods designed to sequence integrated HIV-1 proviruses, such as single-proviral sequencing. FLIPS amplifies and sequences near full-length proviruses enabling replication competency to be determined, and also uses fewer amplification primers, preventing the consequences of primer mismatches. FLIPS is a useful tool for understanding the genetic landscape of integrated HIV-1 proviruses, especially within the latent reservoir, however, its utilization can extend to any application in which the genetic composition of integrated HIV-1 is required.

Viral Discovery in the Invasive Australian Cane Toad (Rhinella marina) Using Metatranscriptomic and Genomic Approaches

Cane toads are poisonous amphibians that were introduced to Australia in 1935 for insect control. Since then, their population has increased dramatically, and they now threaten many native Australian species. One potential method to control the population is to release a cane toad virus with high mortality rates, yet few cane toad viruses have been characterized. This study samples cane toads from different Australian locations and uses an RNA sequencing and computational approach to find new viruses. We report novel complete picornavirus and retrovirus sequences that were genetically similar to viruses infecting frogs, reptiles, and fish. Using data generated in other studies, we show that these viral sequences are present in cane toads from distinct Australian locations. Three sequences related to circoviruses were also found in the toad genome. The identification of new viral sequences will aid future studies that investigate their prevalence and potential as agents for biocontrol.

Cane toads are a notorious invasive species, inhabiting over 1.2 million km² of Australia and threatening native biodiversity. The release of pathogenic cane toad viruses is one possible biocontrol strategy yet is currently hindered by the poorly described cane toad virome. Metatranscriptomic analysis of 16 cane toad livers revealed the presence of a novel and full-length picornavirus, Rhimavirus A (RhiV-A), a member of a reptile- and amphibian-specific cluster of the Picornaviridae basal to the Kobuvirus-like group. In the combined liver transcriptome, we also identified a complete genome sequence of a distinct epsilonretrovirus, Rhinella marina endogenous retrovirus (RMERV). The recently sequenced cane toad genome contains 8 complete RMERV proviruses as well as 21 additional truncated insertions. The oldest full-length RMERV provirus was estimated to have inserted 1.9 million years ago (MYA). To screen for these viral sequences in additional toads, we analyzed publicly available transcriptomes from six diverse Australian locations. RhiV-A transcripts were identified in toads sampled from three locations across 1,000 km of Australia, stretching to the current Western Australia (WA) invasion front, while RMERV transcripts were observed at all six sites. Finally, we scanned the cane toad genome for nonretroviral endogenous viral elements, finding three sequences related to small DNA viruses in the family Circoviridae. This shows ancestral circoviral infection with subsequent genomic integration. The identification of these current and past viral infections enriches our knowledge of the cane toad virome, an understanding of which will facilitate future work on infection and disease in this important invasive species.

IMPORTANCE Cane toads are poisonous amphibians that were introduced to Australia in 1935 for insect control. Since then, their population has increased dramatically, and they now threaten many native Australian species. One potential method to control the population is to release a cane toad virus with high mortality rates, yet few cane toad viruses have been characterized. This study samples cane toads from different Australian locations and uses an RNA sequencing and computational approach to find new viruses. We report novel complete picornavirus and retrovirus sequences that were genetically similar to viruses infecting frogs, reptiles, and fish. Using data generated in other studies, we show that these viral sequences are present in cane toads from distinct Australian locations. Three sequences related to circoviruses were also found in the toad genome. The identification of new viral sequences will aid future studies that investigate their prevalence and potential as agents for biocontrol.

No detectable effect of Wolbachia wMel on the prevalence and abundance of the RNA virome of Drosophila melanogaster

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 wMel strain of Wolbachia can improve host survival following viral infection. However, it is unclear whether wMel 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 wMel-positive and wMel-negative D. melanogaster populations sampled from six locations in Australia. In addition, we examined the impact of wMel on individual flies by obtaining transcriptome data from 20 wMel-positive and 20 wMel-negative D. melanogaster from the location (Melbourne) with highest density of wMel. 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 wMel. Although the viral loads of picorna-like viruses are reduced by wMel 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.

A complex mosaic of enteroviruses shapes community-acquired hand, foot and mouth disease transmission and evolution within a single hospital

Human enteroviruses (EV) pose a major risk to public health. This is especially so in the Asia-Pacific region where increasing numbers of hand, foot and mouth disease (HFMD) cases and large outbreaks of severe neurological disease associated with EV-A71 have occurred. Despite their importance, key aspects of the emergence, epidemiology and evolution of EVs remain unclear, and most studies of EV evolution have focused on a limited number of genes. Here, we describe the genomic-scale evolution of EV-A viruses sampled from pediatric patients with mild disease attending a single hospital in western Sydney, Australia, over an 18-month period. This analysis revealed the presence of eight viral serotypes-Coxsackievirus (CV) A2, A4, A5, A6, A8, A10, A16 and EV-A71-with up to four different serotypes circulating in any 1 month. Despite an absence of large-scale outbreaks, high levels of geographical and temporal mixing of serotypes were identified. Phylogenetic analysis revealed that multiple strains of the same serotype were present in the community, and that this diversity was shaped by multiple introductions into the Sydney population, with only a single lineage of CV-A6 exhibiting in situ transmission over the entire study period. Genomic-scale analyses also revealed the presence of novel and historical EV recombinants. Notably, our analysis revealed no association between viral phylogeny, including serotype, and patient age, sex, nor disease severity (for uncomplicated disease). This study emphasizes the contribution of EV-A viruses other than EV-A71 to mild EV disease including HFMD in Australia and highlights the need for greater surveillance of these viruses to improve strategies for outbreak preparedness and vaccine design.

Identification of Genetically Intact HIV-1 Proviruses in Specific CD4+ T Cells from Effectively Treated Participants

Latent replication-competent HIV-1 persists in individuals on long-term antiretroviral therapy (ART). We developed the Full-Length Individual Proviral Sequencing (FLIPS) assay to determine the distribution of latent replication-competent HIV-1 within memory CD4⁺ T cell subsets in six individuals on long-term ART. FLIPS is an efficient, high-throughput assay that amplifies and sequences near full-length (∼9 kb) HIV-1 proviral genomes and determines potential replication competency through genetic characterization. FLIPS provides a genome-scale perspective that addresses the limitations of other methods that also genetically characterize the latent reservoir. Using FLIPS, we identified 5% of proviruses as intact and potentially replication competent. Intact proviruses were unequally distributed between T cell subsets, with effector memory cells containing the largest proportion of genetically intact HIV-1 proviruses. We identified multiple identical intact proviruses, suggesting a role for cellular proliferation in the maintenance of the latent HIV-1 reservoir.

Virological Sampling of Inaccessible Wildlife with Drones

There is growing interest in characterizing the viromes of diverse mammalian species, particularly in the context of disease emergence. However, little is known about virome diversity in aquatic mammals, in part due to difficulties in sampling. We characterized the virome of the exhaled breath (or blow) of the Eastern Australian humpback whale (Megaptera novaeangliae). To achieve an unbiased survey of virome diversity, a meta-transcriptomic analysis was performed on 19 pooled whale blow samples collected via a purpose-built Unmanned Aerial Vehicle (UAV, or drone) approximately 3 km off the coast of Sydney, Australia during the 2017 winter annual northward migration from Antarctica to northern Australia. To our knowledge, this is the first time that UAVs have been used to sample viruses. Despite the relatively small number of animals surveyed in this initial study, we identified six novel virus species from five viral families. This work demonstrates the potential of UAVs in studies of virus disease, diversity, and evolution.

The evolutionary history of vertebrate RNA viruses

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

The transmission dynamics and diversity of human metapneumovirus in Peru

Background

The transmission dynamics of human metapneumovirus (HMPV) in tropical countries remain unclear. Further understanding of the genetic diversity of the virus could aid in HMPV vaccine design and improve our understanding of respiratory virus transmission dynamics in low‐ and middle‐income countries.

Materials & Methods

We examined the evolution of HMPV in Peru through phylogenetic analysis of 61 full genome HMPV sequences collected in three ecologically diverse regions of Peru (Lima, Piura, and Iquitos) during 2008‐2012, comprising the largest data set of HMPV whole genomes sequenced from any tropical country to date.

Results

We revealed extensive genetic diversity generated by frequent viral introductions, with little evidence of local persistence. While considerable viral traffic between non‐Peruvian countries and Peru was observed, HMPV epidemics in Peruvian locales were more frequently epidemiologically linked with other sites within Peru. We showed that Iquitos experienced greater HMPV traffic than the similar sized city of Piura by both Bayesian and maximum likelihood methods.

Conclusions

There is extensive HMPV genetic diversity even within smaller and relatively less connected cities of Peru and this virus is spatially fluid. Greater diversity of HMPV in Iquitos compared to Piura may relate to higher volumes of human movement, including air traffic to this location.

Emerging recombinant noroviruses identified by clinical and waste water screening

Norovirus is estimated to cause 677 million annual cases of gastroenteritis worldwide, resulting in 210,000 deaths. As viral gastroenteritis is generally self-limiting, clinical samples for epidemiological studies only partially represent circulating noroviruses in the population and is biased towards severe symptomatic cases. As infected individuals from both symptomatic and asymptomatic cases shed viruses into the sewerage system at a high concentration, waste water samples are useful for the molecular epidemiological analysis of norovirus genotypes at a population level. Using Illumina MiSeq and Sanger sequencing, we surveyed circulating norovirus within Australia and New Zealand, from July 2014 to December 2016. Importantly, norovirus genomic diversity during 2016 was compared between clinical and waste water samples to identify potential pandemic variants, novel recombinant viruses and the timing of their emergence. Although the GII.4 Sydney 2012 variant was prominent in 2014 and 2015, its prevalence significantly decreased in both clinical and waste water samples over 2016. This was concomitant with the emergence of multiple norovirus strains, including two GII.4 Sydney 2012 recombinant viruses, GII.P4 New Orleans 2009/GII.4 Sydney 2012 and GII.P16/GII.4 Sydney 2012, along with three other emerging strains GII.17, GII.P12/GII.3 and GII.P16/GII.2. This is unusual, as a single GII.4 pandemic variant is generally responsible for 65–80% of all human norovirus infections at any one time and predominates until it is replaced by a new pandemic variant. In sumary, this study demonstrates the combined use of clinical and wastewater samples provides a more complete picture of norovirus circulating within the population.

Francisella tularensis ssp. holarctica in Ringtail Possums, Australia

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

The paradox of HBV evolution as revealed from a 16th century mummy

Hepatitis B virus (HBV) is a ubiquitous viral pathogen associated with large-scale morbidity and mortality in humans. However, there is considerable uncertainty over the time-scale of its origin and evolution. Initial shotgun data from a mid-16^th century Italian child mummy, that was previously paleopathologically identified as having been infected with Variola virus (VARV, the agent of smallpox), showed no DNA reads for VARV yet did for hepatitis B virus (HBV). Previously, electron microscopy provided evidence for the presence of VARV in this sample, although similar analyses conducted here did not reveal any VARV particles. We attempted to enrich and sequence for both VARV and HBV DNA. Although we did not recover any reads identified as VARV, we were successful in reconstructing an HBV genome at 163.8X coverage. Strikingly, both the HBV sequence and that of the associated host mitochondrial DNA displayed a nearly identical cytosine deamination pattern near the termini of DNA fragments, characteristic of an ancient origin. In contrast, phylogenetic analyses revealed a close relationship between the putative ancient virus and contemporary HBV strains (of genotype D), at first suggesting contamination. In addressing this paradox we demonstrate that HBV evolution is characterized by a marked lack of temporal structure. This confounds attempts to use molecular clock-based methods to date the origin of this virus over the time-frame sampled so far, and means that phylogenetic measures alone cannot yet be used to determine HBV sequence authenticity. If genuine, this phylogenetic pattern indicates that the genotypes of HBV diversified long before the 16^th century, and enables comparison of potential pathogenic similarities between modern and ancient HBV. These results have important implications for our understanding of the emergence and evolution of this common viral pathogen.

Hepatitis B virus (HBV) exerts formidable morbidity and mortality in humans. We used ancient DNA techniques to recover the complete genome sequence of an HBV from the mummified remains of a child discovered in the 16^th century from Naples, Italy. Strikingly, our analysis of this specimen resulted in two contrasting findings: while the damage patterns lend credence to this HBV sequence being authentically 16^th century, phylogenetic analysis revealed a close relationship to recently sampled viruses as expected if the sequence were a modern contaminant. We reconcile these two observations by showing that HBV evolution over the last ~450 years is characterized by a marked lack of temporal structure that hinders attempts to resolve the evolutionary time-scale of this important human pathogen.

Persistent infections in immunocompromised hosts are rarely sources of new pathogen variants

Many viruses, including human norovirus and influenza, cause self-limiting diseases of short duration. However, infection by the same viruses in an immunocompromised host can result in prolonged illness in the absence of effective treatment. Such persistent infections are often characterized by increased genetic diversity with potentially elevated rates of evolution compared to acute infections, leading to suggestions that immunocompromised hosts represent an important reservoir for the emergence of novel viral variants. Here, we develop a mathematical model that combines epidemiological dynamics with within-host evolution to quantify the relative contribution of immunocompromised hosts to the overall rate of pathogen evolution. Using human norovirus as a case study we show that the majority of evolutionary substitutions are expected to occur in acute infections of immunocompetent hosts. Hence, despite their potential to generate a high level of diversity, infections of immunocompromised hosts likely contribute less to the evolution and emergence of new genetic variants at the epidemiological scale because such hosts are rare and tend to be isolated. This result is robust to variation in key parameters, including the proportion of the population immunocompromised, and provides a means to understand the adaptive significance of mutations that arise during chronic infections in immunocompromised hosts.

Emergent multisystemic Enterococcus infection threatens endangered Christmas Island reptile populations

Multisystemic infections with a morphologically unusual bacterium were first observed in captive critically endangered Lister’s geckos (Lepidodactylus listeri) on Christmas Island in October 2014. Since then the infection was identified in another captive critically endangered lizard species, the blue-tailed skink (Cryptoblepharus egeriae) and two species of invasive geckos; the four clawed gecko (Gehyra mutilata) and Asian house gecko (Hemidactylus frenatus), in a wide geographic range across the east side of the island. The Gram and periodic acid-Schiff positive cocci to diplococci have a propensity to form chains surrounded by a matrix, which ultrastructurally appears to be formed by fibrillar capsular projections. The bacterium was associated with severe and extensive replacement of tissues, but minimal host inflammatory response. Attempts to grow the organism in culture and in embryonated eggs were unsuccessful. Molecular characterisation of the organism placed it as a novel member of the genus Enterococcus. Disease Risk Analyses including this organism should now be factored into conservation management actions and island biosecurity.

Genomic and phenotypic characterization of myxoma virus from Great Britain reveals multiple evolutionary pathways distinct from those in Australia

The co-evolution of myxoma virus (MYXV) and the European rabbit occurred independently in Australia and Europe from different progenitor viruses. Although this is the canonical study of the evolution of virulence, whether the genomic and phenotypic outcomes of MYXV evolution in Europe mirror those observed in Australia is unknown. We addressed this question using viruses isolated in the United Kingdom early in the MYXV epizootic (1954-1955) and between 2008-2013. The later UK viruses fell into three distinct lineages indicative of a long period of separation and independent evolution. Although rates of evolutionary change were almost identical to those previously described for MYXV in Australia and strongly clock-like, genome evolution in the UK and Australia showed little convergence. The phenotypes of eight UK viruses from three lineages were characterized in laboratory rabbits and compared to the progenitor (release) Lausanne strain. Inferred virulence ranged from highly virulent (grade 1) to highly attenuated (grade 5). Two broad disease types were seen: cutaneous nodular myxomatosis characterized by multiple raised secondary cutaneous lesions, or an amyxomatous phenotype with few or no secondary lesions. A novel clinical outcome was acute death with pulmonary oedema and haemorrhage, often associated with bacteria in many tissues but an absence of inflammatory cells. Notably, reading frame disruptions in genes defined as essential for virulence in the progenitor Lausanne strain were compatible with the acquisition of high virulence. Combined, these data support a model of ongoing host-pathogen co-evolution in which multiple genetic pathways can produce successful outcomes in the field that involve both different virulence grades and disease phenotypes, with alterations in tissue tropism and disease mechanisms.