2. Abstract Metagenomic next-generation sequencing has transformed the discovery and diagnosis of infectious disease, with the power to characterize 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 libraries of ‘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 libraries, as well as contaminating sequences from the RNA virus families Totiviridae, Tombusviridae and Lentiviridae . These data suggest that the contamination of common laboratory reagents is likely widespread and can comprise a wide variety of viruses. 3. Data summary The authors confirm all supporting data, code and protocols have been provided within the article or through supplementary data files. 1.5 Repositories The viral genome sequence data generated in this study has been deposited in the NCBI database under accession numbers MZ824225 - MZ824237 . Sequence reads are available at the public Sequence Read Archive (SRA) database with accession SRX6803604 and under the BioProject accession PRJNA735051 reference numbers SRR14737466-71 and BioSample numbers SAMN20355437-40.

Abstract Viral pathogens are being increasingly described in association with mass morbidity and mortality events in reptiles. However, our knowledge of reptile viruses and their role in population health remains limited. Herein, we describe a 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 of bearded dragon tissue 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 contrast, the newly discovered BDchPV showed approximately 35-40% identity in the non-structural (NS) protein to parvoviruses obtained from tilapia fish and crocodiles in China. Subsequent specific PCR assays detected BDCV exclusively and comprehensively within animals with proliferative pulmonary lesions and respiratory disease. This study expands our understanding of viral diversity in the context of diseased reptiles in captivity.

ABSTRACT 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, we performed 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. 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.

Abstract 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 viral 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 based on amino acid sequence data. Using this combination of meta-transcriptomics and protein structure prediction we identified the RdRp (PB1) gene segment of a divergent negative-sense RNA virus in a native Australian gecko ( Geyra lauta ) that was confirmed by PCR and Sanger sequencing. Phylogenetic analysis identified the Gecko articulavirus (GECV) as 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 .

Summary The HIV-1 reservoir is composed of cells harboring latent proviruses that are capable of contributing to viremia upon antiretroviral treatment ( ART) interruption. Although this reservoir is known to be maintained by clonal expansion, the contribution of large, infected cell clones to residual viremia and viral rebound remains underexplored. Here, we conducted an extensive analysis on four ART-treated individuals who underwent an analytical treatment interruption ( ATI ). We performed subgenomic (V1-V3 env ), near full-length proviral and integration site sequencing, and used multiple displacement amplification to sequence both the integration site and provirus from single HIV-infected cells. We found eight proviruses that could phylogenetically be linked to plasma virus obtained before or during the ATI. This study highlights a role for HIV-infected cell clones in the maintenance of the replication-competent reservoir and suggests that infected cell clones can directly contribute to rebound viremia upon ATI.

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

Abstract The RNA virus family Flaviviridae harbours several 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 of which have a diverse host range. Members of the genus Hepacivirus are associated with a diverse array of animal species, including humans and 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 including humans. Using a combination of meta-transcriptomic and whole genome sequencing we identified four novel hepaciviruses and one novel variant of a known virus, in five species of native Australian wildlife, expanding our knowledge of the diversity in this important group of RNA viruses. The infected hosts comprised native Australian marsupials and birds, as well as a native gecko ( Gehyra lauta ). The addition of these novel viruses led to the identification of a distinct marsupial clade within the hepacivirus phylogeny that also included an engorged Ixodes holocyclus tick collected while feeding on Australian long-nosed bandicoots ( Perameles nasuta ). Gecko and avian associated hepacivirus lineages were also identified. In addition, by mining the short-read archive (SRA) database we identified another five novel members of Flaviviridae , comprising three new hepaciviruses from avian and primate hosts, as well as two primate-associated pegiviruses. The large-scale phylogenetic analysis of these novel hepacivirus and pegivirus genomes provides additional support for virus-host co-divergence over evolutionary time-scales.

Hepatitis delta virus (HDV) is the smallest known RNA virus and encodes 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. Notably, the novel viruses identified here share HDV-like genomic features such as a small genome size of ~1.7kb in length, circular genomes, 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 HDAgs. Notably, 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 in non-human species, but shed light on its origin and evolutionary history.

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. We used bulk RNA-sequencing (meta-transcriptomics) to assess the diversity and expression levels of functionally active resistance genes in the microbiome of birds with aquatic behavior. We sampled birds across a range of habitats, from penguins in Antarctica to ducks in a wastewater treatment plant in Australia. This revealed 81 antibiotic resistance genes in birds from all localities, including β-lactam, tetracycline and chloramphenicol resistance in Antarctica, and genes typically associated with multidrug resistance plasmids in areas with high human impact. Notably, birds feeding at a wastewater treatment plant carried the greatest resistance gene burden, suggesting that human waste, even if it undergoes treatment, contributes to the spread of antibiotic resistance genes to the wild. Differences in resistance gene burden also reflected the birds' ecology, taxonomic group and microbial functioning. Ducks, which feed by dabbling, carried a higher abundance and diversity of resistance genes than turnstones, avocets and penguins, that usually prey on more pristine waters. In sum, this study helps to reveal the complex factors explaining the distribution of resistance genes and their exchange routes between humans and wildlife.

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 which 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 subviral agents.