Contrasting dynamics of two incursions of low-pathogenicity avian influenza virus into Australia.

The current panzootic of high pathogenicity avian influenza virus H5N1 demonstrates how viral incursions can have major ramifications for wildlife and domestic animals. Herein, we describe the recent incursion into Australia of two low pathogenicity avian influenza virus subtypes, H4 and H10, that exhibited contrasting evolutionary dynamics. Viruses detected from national surveillance and disease investigations between 2020 and 2022 revealed 27 genomes, 24 of which have at least one segment more closely related to Eurasian or North American avian influenza lineages than those already circulating in Australia. Phylogenetic analysis revealed that H4 viruses circulating in shorebirds represent a recent incursion from Asia that is distinct from those circulating concurrently in Australian waterfowl. Analysis of the internal segments further demonstrates exclusive, persistent circulation in shorebirds. This contrasts with H10, where a novel lineage has emerged in wild waterfowl, poultry, and captive birds across Australia and has likely replaced previously circulating H10 lineages through competitive exclusion. Elucidating different dynamics for avian influenza incursions supports effective disease risk identification and communication that better informs disease preparedness and response.

Using artificial intelligence to document the hidden RNA virosphere.

Current metagenomic tools can fail to identify highly divergent RNA viruses. We developed a deep learning algorithm, termed LucaProt, to discover highly divergent RNA-dependent RNA polymerase (RdRP) sequences in 10,487 metatranscriptomes generated from diverse global ecosystems. LucaProt integrates both sequence and predicted structural information, enabling the accurate detection of RdRP sequences. Using this approach, we identified 161,979 potential RNA virus species and 180 RNA virus supergroups, including many previously poorly studied groups, as well as RNA virus genomes of exceptional length (up to 47,250 nucleotides) and genomic complexity. A subset of these novel RNA viruses was confirmed by RT-PCR and RNA/DNA sequencing. Newly discovered RNA viruses were present in diverse environments, including air, hot springs, and hydrothermal vents, with virus diversity and abundance varying substantially among ecosystems. This study advances virus discovery, highlights the scale of the virosphere, and provides computational tools to better document the global RNA virome.

A virus associated with the zoonotic pathogen Plasmodium knowlesi causing human malaria is a member of a diverse and unclassified viral taxon.

Apicomplexa are single-celled eukaryotes that can infect humans and include the mosquito-borne parasite Plasmodium, the cause of malaria. Increasing rates of drug resistance in human-only Plasmodium species are reducing the efficacy of control efforts and antimalarial treatments. There are also rising cases of P. knowlesi, the only zoonotic Plasmodium species that causes severe disease and death in humans. Thus, there is a need to develop additional innovative strategies to combat malaria. Viruses that infect non-Plasmodium spp. Disease-causing protozoa have been shown to affect pathogen life cycle and disease outcomes. However, only one virus (Matryoshka RNA virus 1) has been identified in Plasmodium, and none have been identified in zoonotic Plasmodium species. The rapid expansion of the known RNA virosphere using structure- and artificial intelligence-based methods suggests that this dearth is due to the divergent nature of RNA viruses that infect protozoa. We leveraged these newly uncovered data sets to explore the virome of human-infecting Plasmodium species collected in Sabah, east (Borneo) Malaysia. We identified a highly divergent RNA virus in two human-infecting P. knowlesi isolates that is related to the unclassified group 'ormycoviruses'. By characterising fifteen additional ormycoviruses identified in the transcriptomes of arthropods we show that this group of viruses exhibits a complex ecology at the arthropod-mammal interface. Through the application of artificial intelligence methods, we then demonstrate that the ormycoviruses are part of a diverse and unclassified viral taxon. This is the first observation of an RNA virus in a zoonotic Plasmodium species. By linking small-scale experimental data to large-scale virus discovery advances, we characterise the diversity and genomic architecture of an unclassified viral taxon. This approach should be used to further explore the virome of disease-causing Apicomplexa and better understand how protozoa-infecting viruses may affect parasite fitness, pathobiology, and treatment outcomes.

Virology-The next fifty years.

Virology has made enormous advances in the last 50 years but has never faced such scrutiny as it does today. Herein, we outline some of the major advances made in virology during this period, particularly in light of the COVID-19 pandemic, and suggest some areas that may be of research importance in the next 50 years. We focus on several linked themes: cataloging the genomic and phenotypic diversity of the virosphere; understanding disease emergence; future directions in viral disease therapies, vaccines, and interventions; host-virus interactions; the role of viruses in chronic diseases; and viruses as tools for cell biology. We highlight the challenges that virology will face moving forward-not just the scientific and technical but also the social and political. Although there are inherent limitations in trying to outline the virology of the future, we hope this article will help inspire the next generation of virologists.

Farmed fur animals harbour viruses with zoonotic spillover potential.

Animals such as raccoon dogs, mink and muskrats are farmed for fur and are sometimes used as food or medicinal products1,2, yet they are also potential reservoirs of emerging pathogens3. Here we performed single-sample metatranscriptomic sequencing of internal tissues from 461 individual fur animals that were found dead due to disease. We characterized 125 virus species, including 36 that were novel and 39 at potentially high risk of cross-species transmission, including zoonotic spillover. Notably, we identified seven species of coronaviruses, expanding their known host range, and documented the cross-species transmission of a novel canine respiratory coronavirus to raccoon dogs and of bat HKU5-like coronaviruses to mink, present at a high abundance in lung tissues. Three subtypes of influenza A virus-H1N2, H5N6 and H6N2-were detected in the lungs of guinea pig, mink and muskrat, respectively. Multiple known zoonotic viruses, such as Japanese encephalitis virus and mammalian orthoreovirus4,5, were detected in guinea pigs. Raccoon dogs and mink carried the highest number of potentially high-risk viruses, while viruses from the Coronaviridae, Paramyxoviridae and Sedoreoviridae families commonly infected multiple hosts. These data also reveal potential virus transmission between farmed animals and wild animals, and from humans to farmed animals, indicating that fur farming represents an important transmission hub for viral zoonoses.

A new lineage nomenclature to aid genomic surveillance of dengue virus.

Dengue virus (DENV) is currently causing epidemics of unprecedented scope in endemic settings and expanding to new geographical areas. It is therefore critical to track this virus using genomic surveillance. However, the complex patterns of viral genomic diversity make it challenging to use the existing genotype classification system. Here, we propose adding 2 sub-genotypic levels of virus classification, named major and minor lineages. These lineages have high thresholds for phylogenetic distance and clade size, rendering them stable between phylogenetic studies. We present assignment tools to show that the proposed lineages are useful for regional, national, and subnational discussions of relevant DENV diversity. Moreover, the proposed lineages are robust to classification using partial genome sequences. We provide a standardized neutral descriptor of DENV diversity with which we can identify and track lineages of potential epidemiological and/or clinical importance. Information about our lineage system, including methods to assign lineages to sequence data and propose new lineages, can be found at: dengue-lineages.org.

Identification of Corynebacterium ulcerans and Erysipelothrix sp. in Malayan pangolins-a potential threat to public health?

The discovery of severe acute respiratory syndrome-coronavirus-2-like and Middle East respiratory syndrome-coronavirus-like viruses in Malayan pangolins has raised concerns about their potential role in the spread of zoonotic diseases. Herein, we describe the isolation and whole-genome sequencing of potentially zoonotic two bacterial pathogens from diseased Malaysian pangolins (Manis javanica)-Corynebacterium ulcerans and Erysipelothrix sp. The newly identified species were designated as C. ulcerans P69 and Erysipelothrix sp. P66. C. ulcerans P69 exhibited 99.2% whole-genome nucleotide identity to human bacterial isolate 4940, suggesting that it might have zoonotic potential. Notably, C. ulcerans P69 lacked the diphtheria toxin (tox) gene that is widely used in vaccines to protect humans from corynebacterial infection, which suggests that the current vaccine may be of limited efficacy against this pangolin strain. C. ulcerans P69 also contains other known virulence-associated genes such as pld and exhibits resistance to several antibiotics (erythromycin, clindamycin, penicillin G, gentamicin, tetracycline), which may affect its effective control. Erysipelothrix sp. P66 was closely related to Erysipelothrix sp. strain 2-related strains, exhibiting 98.8% whole-genome nucleotide identity. This bacterium is lethal in mice, and two commercial vaccines failed to protect its challenge, such that it could potentially pose a threat to the swine industry. Overall, this study highlights that, in addition to viruses, pangolins harbor bacteria that may pose a potential threat to humans and domestic animals, and which merit attention. IMPORTANCE: This study firstly reports the presence of two potentially zoonotic bacteria, Corynebacterium ulcerans and Erysipelothrix sp., in diseased Malaysian pangolins collected in 2019. The pangolin C. ulcerans is lethal in mice and resists many antibiotics. It clustered with a lethal human strain but lacked the diphtheria toxin gene. Diphtheria toxin is widely used as a vaccine around the world to protect humans from the infection of corynebacteria. The lack of the tox gene suggests that the current vaccine may be of limited efficacy against this pangolin strain. The pangolin Erysipelothrix sp. is the sister clade of Erysipelothrix rhusiopathiae. It is lethal in mice, and two commercial vaccines failed to protect the mice against challenge with the pangolin Erysipelothrix sp., such that this strain could potentially pose a threat to the swine industry. These findings emphasize the potential threat of pangolin bacteria.

Genetic tracing of market wildlife and viruses at the epicenter of the COVID-19 pandemic.

Zoonotic spillovers of viruses have occurred through the animal trade worldwide. The start of the COVID-19 pandemic was traced epidemiologically to the Huanan Seafood Wholesale Market. Here, we analyze environmental qPCR and sequencing data collected in the Huanan market in early 2020. We demonstrate that market-linked severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genetic diversity is consistent with market emergence and find increased SARS-CoV-2 positivity near and within a wildlife stall. We identify wildlife DNA in all SARS-CoV-2-positive samples from this stall, including species such as civets, bamboo rats, and raccoon dogs, previously identified as possible intermediate hosts. We also detect animal viruses that infect raccoon dogs, civets, and bamboo rats. Combining metagenomic and phylogenetic approaches, we recover genotypes of market animals and compare them with those from farms and other markets. This analysis provides the genetic basis for a shortlist of potential intermediate hosts of SARS-CoV-2 to prioritize for serological and viral sampling.

Mapping glycoprotein structure reveals Flaviviridae evolutionary history.

Viral glycoproteins drive membrane fusion in enveloped viruses and determine host range, tissue tropism and pathogenesis1. Despite their importance, there is a fragmentary understanding of glycoproteins within the Flaviviridae2, a large virus family that include pathogens such as hepatitis C, dengue and Zika viruses, and numerous other human, animal and emergent viruses. For many flaviviruses the glycoproteins have not yet been identified, for others, such as the hepaciviruses, the molecular mechanisms of membrane fusion remain uncharacterized3. Here we combine phylogenetic analyses with protein structure prediction to survey glycoproteins across the entire Flaviviridae. We find class II fusion systems, homologous to the Orthoflavivirus E glycoprotein in most species, including highly divergent jingmenviruses and large genome flaviviruses. However, the E1E2 glycoproteins of the hepaciviruses, pegiviruses and pestiviruses are structurally distinct, may represent a novel class of fusion mechanism, and are strictly associated with infection of vertebrate hosts. By mapping glycoprotein distribution onto the underlying phylogeny, we reveal a complex evolutionary history marked by the capture of bacterial genes and potentially inter-genus recombination. These insights, made possible through protein structure prediction, refine our understanding of viral fusion mechanisms and reveal the events that have shaped the diverse virology and ecology of the Flaviviridae.

VirID: Beyond Virus Discovery-An Integrated Platform for Comprehensive RNA Virus Characterization.

RNA viruses exhibit vast phylogenetic diversity and can significantly impact public health and agriculture. However, current bioinformatics tools for viral discovery from metagenomic data frequently generate false positive virus results, overestimate viral diversity, and misclassify virus sequences. Additionally, current tools often fail to determine virus-host associations, which hampers investigation of the potential threat posed by a newly detected virus. To address these issues we developed VirID, a software tool specifically designed for the discovery and characterization of RNA viruses from metagenomic data. The basis of VirID is a comprehensive RNA-dependent RNA polymerase database to enhance a workflow that includes RNA virus discovery, phylogenetic analysis, and phylogeny-based virus characterization. Benchmark tests on a simulated data set demonstrated that VirID had high accuracy in profiling viruses and estimating viral richness. In evaluations with real-world samples, VirID was able to identify RNA viruses of all types, but also provided accurate estimations of viral genetic diversity and virus classification, as well as comprehensive insights into virus associations with humans, animals, and plants. VirID therefore offers a robust tool for virus discovery and serves as a valuable resource in basic virological studies, pathogen surveillance, and early warning systems for infectious disease outbreaks.

Analysis of the genetic diversity in RNA-directed RNA polymerase sequences: implications for an automated RNA virus classification system.

RNA viruses are characterized by a broad host range and high levels of genetic diversity. Despite a recent expansion in the known virosphere following metagenomic sequencing, our knowledge of the species rank genetic diversity of RNA viruses, and how often they are misassigned and misclassified, is limited. We performed a clustering analysis of 7801 RNA-directed RNA polymerase (RdRp) sequences representing 1897 established RNA virus species. From this, we identified substantial genetic divergence within some virus species and inconsistency in RNA virus assignment between the GenBank database and The International Committee on Taxonomy of Viruses (ICTV). In particular, 27.57% virus species comprised multiple virus operational taxonomic units (vOTUs), including Alphainfluenzavirus influenzae, Mammarenavirus lassaense, Apple stem pitting virus, and Rotavirus A, with each having over 100 vOTUs. In addition, the distribution of average amino acid identity between vOTUs within single assigned species showed a relatively low threshold: <90% and sometimes <50%. However, when only exemplar sequences from virus species were analyzed, 1889 of the ICTV-designated RNA virus species (99.58%) were clustered into a single vOTU. Clustering of the RdRp sequences from different virus species also revealed that 17 vOTUs contained two distinct virus species. These potential misassignments were confirmed by phylogenetic analysis. A further analysis of average nucleotide identity (ANI) values ranging from 70% to 97.5% revealed that at an ANI of 82.5%, 1559 (82.18%) of the 1897 virus species could be correctly clustered into one single vOTU. However, at ANI values >82.5%, an increasing number of species were clustered into two or more vOTUs. In sum, we have identified some inconsistency and misassignment of the RNA virus species based on the analysis of RdRp sequences alone, which has important implications for the development of an automated RNA virus classification system.

Five-year longitudinal surveillance reveals the continual circulation of both alpha- and beta-coronaviruses in Plateau and Gansu pikas (Ochotona spp.) at Qinghai Lake, China1.

The discovery of alphacoronaviruses and betacoronaviruses in plateau pikas (Ochotona curzoniae) expanded the host range of mammalian coronavirus (CoV) to a new order - Lagomorpha. However, the diversity and evolutionary relationships of CoVs in these plateau-region-specific animal population remains uncertain. We conducted a five-year longitudinal surveillance of CoVs harboured by pikas around Qinghai Lake, China. CoVs were identified in 33 of 236 plateau pikas and 2 of 6 Gansu pikas (Ochotona cansus), with a total positivity rate of 14.5%, and exhibiting a wide spatiotemporal distribution across seven sampling sites and six time points. Through meta-transcriptomic sequencing and RT-PCR, we recovered 16 near-complete viral genome sequences. Phylogenetic analyses classified the viruses as variants of either pika alphacoronaviruses or betacoronaviruses endemic to plateau pikas from the Qinghai-Tibet Plateau region. Of particular note, the pika-associated betacoronaviruses may represent a novel subgenus within the genus Betacoronavirus. Tissue tropism, evaluated using quantitative real-time PCR, revealed the presence of CoV in the rectal and/or lung tissues, with the highest viral loads at 103.55 or 102.80 RNA copies/µL. Surface plasmon resonance (SPR) assays indicated that the newly identified betacoronavirus did not bind to human or pika Angiotensin-converting enzyme 2 (ACE2) or Dipeptidyl peptidase 4 (DPP4). The findings highlight the ongoing circulation and broadening host spectrum of CoVs among pikas, emphasizing the necessity for further investigation to evaluate their potential public health risks.

A ~40-kb flavi-like virus does not encode a known error-correcting mechanism.

It is commonly held that there is a fundamental relationship between genome size and error rate, manifest as a notional "error threshold" that sets an upper limit on genome sizes. The genome sizes of RNA viruses, which have intrinsically high mutation rates due to a lack of mechanisms for error correction, must therefore be small to avoid accumulating an excessive number of deleterious mutations that will ultimately lead to population extinction. The proposed exceptions to this evolutionary rule are RNA viruses from the order Nidovirales (such as coronaviruses) that encode error-correcting exonucleases, enabling them to reach genome lengths greater than 40 kb. The recent discovery of large-genome flavi-like viruses (Flaviviridae), which comprise genomes up to 27 kb in length yet seemingly do not encode exonuclease domains, has led to the proposal that a proofreading mechanism is required to facilitate the expansion of nonsegmented RNA virus genomes above 30 kb. Herein, we describe a ~40 kb flavi-like virus identified in a Haliclona sponge metatranscriptome that does not encode a known exonuclease. Structural analysis revealed that this virus may have instead captured cellular domains associated with nucleic acid metabolism that have not been previously found in RNA viruses. Phylogenetic inference placed this virus as a divergent pesti-like lineage, such that we have provisionally termed it "Maximus pesti-like virus." This virus represents an instance of a flavi-like virus achieving a genome size comparable to that of the Nidovirales and demonstrates that RNA viruses have evolved multiple solutions to overcome the error threshold.

Meta-transcriptomic analysis of companion animal infectomes reveals their diversity and potential roles in animal and human disease.

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

Viromes of Antarctic fish resemble the diversity found at lower latitudes.

Antarctica harbours some of the most isolated and extreme environments on Earth, concealing a largely unexplored and unique component of the global animal virosphere. To understand the diversity and evolutionary histories of viruses in these polar species, we determined the viromes of gill metatranscriptomes from 11 Antarctic fish species with 248 samples collected from the Ross Sea region spanning the Perciformes, Gadiformes, and Scorpaeniformes orders. The continent's shift southward and cooling temperatures >20 million years ago led to a reduction in biodiversity and subsequent radiation of some marine fauna, such as the notothenioid fishes. Despite decreased host species richness in polar regions, we revealed a surprisingly complex virome diversity in Ross Sea fish, with the types and numbers of viruses per host species and individuals sampled comparable to that of fish in warmer marine environments with higher host community diversity. We also observed a higher number of closely related viruses likely representing instances of recent and historic host-switching events among the Perciformes (all notothenioids) than in the Gadiformes, suggesting that rapid speciation events within this order generated closely related host species with few genetic barriers to cross-species transmission. Additionally, we identified novel genomic variation in an arenavirus with a split nucleoprotein sequence containing a stable helical structure, indicating potential adaptation of viral proteins to extreme temperatures. These findings enhance our understanding of virus evolution and virus-host interactions in response to environmental shifts, especially in less diverse ecosystems that are more vulnerable to the impacts of anthropogenic and climate changes.

From islands to infectomes: host-specific viral diversity among birds across remote islands.

BACKGROUND: Accelerating biodiversity loss necessitates monitoring the potential pathogens of vulnerable species. With a third of New Zealand's avifauna considered at risk of extinction, a greater understanding of the factors that influence microbial transmission in this island ecosystem is needed. We used metatranscriptomics to determine the viruses, as well as other microbial organisms (i.e. the infectomes), of seven bird species, including the once critically endangered black robin (Petroica traversi), on two islands in the remote Chatham Islands archipelago, New Zealand. RESULTS: We identified 19 likely novel avian viruses across nine viral families. Black robins harboured viruses from the Flaviviridae, Herpesviridae, and Picornaviridae, while introduced starlings (Sturnus vulgaris) and migratory seabirds (Procellariiformes) carried viruses from six additional viral families. Potential cross-species virus transmission of a novel passerivirus (family: Picornaviridae) between native (black robins and grey-backed storm petrels) and introduced (starlings) birds was also observed. Additionally, we identified bacterial genera, apicomplexan parasites, as well as a novel megrivirus linked to disease outbreaks in other native New Zealand birds. Notably, island effects were outweighed by host taxonomy as a significant driver of viral composition, even among sedentary birds. CONCLUSIONS: These findings underscore the value of surveillance of avian populations to identify and minimise escalating threats of disease emergence and spread in these island ecosystems. Importantly, they contribute to our understanding of the potential role of introduced and migratory birds in the transmission of microbes and associated diseases, which could impact vulnerable island-endemic species.

Diverse RNA viruses in the venom-related microenvironment of different animal phyla.

Venom is known as the source of natural antimicrobial products. Previous studies have largely focused on the expression of venom-related genes and the biochemical components of venom. With the advent of metagenomic sequencing, many more microorganisms, especially viruses, have been identified in highly diverse environments. Herein, we investigated the RNA virome in the venom-related microenvironment through analysis of a large volume of venom-related RNA-sequencing data mined from public databases. From this, we identified viral sequences belonging to thirty-six different viruses, of which twenty-two were classified as 'novel' as they exhibited less than 90 per cent amino acid identity to known viruses in the RNA-dependent RNA polymerase. Most of these novel viruses possessed genome structures similar to their closest relatives, with specific alterations in some cases. Phylogenetic analyses revealed that these viruses belonged to at least twenty-two viral families or unclassified groups, some of which were highly divergent from known taxa. Although further analysis failed to find venom-specific viruses, some viruses seemingly had much higher abundance in the venom-related microenvironment than in other tissues. In sum, our study provides insights into the RNA virome of the venom-related microenvironment from diverse animal phyla.

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.

The radiation of New Zealand's skinks and geckos is associated with distinct viromes.

BACKGROUND: New Zealand is home to over 120 native endemic species of skinks and geckos that radiated over the last 20-40 million years, likely driven by the exploitation of diverse habitats formed during the Miocene. The recent radiation of animal hosts may facilitate cross-species virus transmission, likely reflecting their close genetic relationships and therefore relatively low barriers for viruses to emerge in new hosts. Conversely, as animal hosts adapt to new niches, even within specific geographic locations, so too could their viruses. Consequently, animals that have niche-specialised following radiations may be expected to harbour genetically distinct viruses. Through a metatranscriptomic analysis of eight of New Zealand's native skink and gecko species, as well as the only introduced lizard species, the rainbow skink (Lampropholis delicata), we aimed to reveal the diversity of viruses in these hosts and determine whether and how the radiation of skinks and geckos in New Zealand has impacted virus diversity and evolution. RESULTS: We identified a total of 15 novel reptilian viruses spanning 11 different viral families, across seven of the nine species sampled. Notably, we detected no viral host-switching among the native animals analysed, even between those sampled from the same geographic location. This is compatible with the idea that host speciation has likely resulted in isolated, niche-constrained viral populations that have prevented cross-species transmission. Using a protein structural similarity-based approach, we further identified a highly divergent bunya-like virus that potentially formed a new family within the Bunyavirales. CONCLUSIONS: This study has broadened our understanding of reptilian viruses within New Zealand and illustrates how niche adaptation may limit viral-host interactions.