Emerging viral pathogens are a significant concern, with potential consequences for human, animal and environmental health. Over the past several decades, many novel viruses have been found in animals, including birds, and often pose a significant threat to vulnerable species. However, despite enormous interest in virus research, little is known about virus communities (viromes) in Australian Neophema birds. Therefore, this study was designed to characterise the viromes of Neophema birds and track the evolutionary relationships of recently emerging psittacine siadenovirus F (PsSiAdV-F) circulating in the critically endangered, orange-bellied parrot (OBP, Neophema chrysogaster), using a viral metagenomic approach. This study identified 16 viruses belonging to the families Adenoviridae, Circoviridae, Endornaviridae, Picobirnaviridae and Picornaviridae. In addition, this study demonstrated a potential evolutionary relationship of a PsSiAdV-F sequenced previously from the critically endangered OBP. Strikingly, five adenoviral contigs identified in this study show the highest identities with human adenovirus 2 and human mastadenovirus C. This highlights an important and unexpected aspects of the avian virome and warrants further studies dedicated to this subject. Finally, the findings of this study emphasise the importance of testing birds used for trade or in experimental settings for potential pathogens to prevent the spread of infections.
Over the past several decades, a significant effort has been made to better understand viromes in human and animals, but our understanding of the diversity of viruses in birds is not well understood. The vast majority of research aimed at describing viruses in birds has been on zoonotic viruses (e.g., avian influenza A virus)1, viruses that cause an economic loss in commercial poultry (e.g., Newcastle disease virus, infectious bronchitis virus)2,3, or viruses that result in notable mortality in wild birds (e.g., Wellfleet bay virus, and beak and feather disease virus)4,5,6,7,8. However, little is known about the viromes of birds that are neither economically important nor harbour any known zoonotic pathogens. There is a very limited understanding on the viromes of Australian birds chosen for trade or in vivo experiments, and order Psittaciformes (parrots) is the most traded group among avian orders9. The elegant parrot (Neophema elegans) and the scarlet-chested parrot (Neophema splendida) used in this study were purchased from a commercial trader without knowing any history of diseases.
Metagenomics, or metatranscriptomics, is a relatively new technique that enables the detection and characterisation of entire viromes in animals10,11,12, rather than on single species of virus in isolation. Before metagenomics, there was limited understanding of the viromes present in animal and human hosts. With the advent of metagenomics, information on eukaryotic and prokaryotic viruses, and even on viruses that infect other viruses, has increased8,13,14,15,16,17. However, little is known about the viruses that infect birds of the genus Neophema, an Australian genus with six parrot species, including one of the most critically endangered: the orange-bellied parrot (Neophema chrysogaster). Almost nothing is known about the viruses of the elegant parrot (Neophema elegans) and the scarlet-chested parrot (Neophema splendida) of the genus Neophema, a number of which were housed in La Trobe Animal Research and Teaching Facility for the ‘Parrot Genome Sequencing Project’. The handling of these birds may facilitate the transmission of infectious diseases and may serve as a source of exposure to other birds and to humans involved in their care18,19. Consequently, birds used for trade or in experiments may potentially cause the emergence of pathogens in captive birds, their dispersal into disease-free regions, and the transmission of zoonotic diseases.
In addition, a psittacine siadenovirus F (PsSiAdV-F) has recently been characterised from the critically endangered, orange-bellied parrot (Neophema chrysogaster)20. Subsequent phylogenetic analyses demonstrated that novel PsSiAdV-F infecting orange-bellied parrot evolved prior to all known members of siadenoviruses except the frog siadenovirus A, and did not show any obvious close evolutionary relationship20. Therefore, the aims of this study were not only to identify the known and unknown viruses that may exist in the two species of Neophema birds but also to track the evolutionary relationships of recently emerging PsSiAdV-F circulating in the critically endangered, orange-bellied parrot.
Members of the family Adenoviridae are non-enveloped double-stranded DNA (dsDNA) viruses with a linear genome that ranges from about 26 kb to 46 kb21,22,23. They have been divided into six genera (Aviadenovirus, Atadenovirus, Siadenovirus, Mastadenovirus, Ichtadenovirus and Testadenovirus)24,25,26,27,28,29. They are recognised to cause a wide spectrum of illnesses, ranging from asymptomatic infections to severe illnesses and death in many animals and humans, and are integrated into surveillance programs given their importance in public health30,31. Multiple adenovirus sequences of the genus Siadenovirus and Mastadenovirus were detected in this study.
The complete genome of a psittacine siadenovirus F, PsSiAdV-F, detected in this study was a linear dsDNA molecule of 25,616 bp in length (coverage 127.58×), with a balanced G + C content (55.6%). Based on the complete genome identities, the known AdV genomes that were most closely related to the PsSiAdV-F strain S10/AU were PsSiAdV-F strain OBP2209 (PsSiAdV-F; 99.8%), followed by PsSiAdV-F strain WVL19065-01E (PsSiAdV-F; 99.7%), skua siadenovirus A (SuAdV-A; 57.0%) and raptor siadenovirus A (RAdV-A; 55.6%)28 (Supplementary Table S1). The PsSiAdV-F strain S10/AU genome had 30 predicted methionine-initiated open reading frames (ORFs) encoding proteins that were annotated as putative genes and were numbered from left to right (Table 1). Comparative analysis of the protein sequences encoded by the predicted ORFs, using BLASTX and BLASTP, identified homologs with significant protein sequence similarity for 25 ORFs (Table 1), while five ORFs (ORFs A-E) were found to be unique according to the BLAST database. Among the predicted protein coding ORFs of the PsSiAdV-F strain S10/AU genome, 25 were homologs to a recently sequenced PsSiAdV-F strain OBP2209 gene products (Table 1). The amino acid sequence similarity of predicted ORFs was very high compared to other siadenoviruses, ranging from 98.7 to 100% (Table 1).
Phylogenetic analyses based on two non-structural (polymerase and pTP) and two structural (penton and hexon) protein sequences clearly supported the inclusion of the newly sequenced PsSiAdV-F strain S10/AU in the genus Siadenovirus. In the resulting ML tree, based on concatenated amino acid sequences of four selected AdVs genes, the sequenced PsSiAdV-F strain S10/AU was positioned in a distinct subclade with other two PsSiAdV-F strains (e.g., OBP2209 and WVL19065-01E) (100% bootstrap support) (Fig. 1). Considering the genome-wide identities (> 99%) and phylogenetic position among three strains of PsSiAdV-F, it is reasonable to postulate that the PsSiAdV-F originates from the same ancestor.
Strikingly, this study detected five contigs that were shown to be homologs with human adenoviruses. The largest contig of 958 nt (coverage 8.25×, GenBank accession no. MZ364301) spanned part of the human adenovirus 2 penton gene. From further analysis, the human adenovirus 2 penton gene sequence had a > 99% (930/933) nucleotide match to EU128937.1 (Supplementary Table S2). At the amino acid level, this segment had a 296/298 (99%) amino acid match to the penton gene of the human adenovirus 2 (GenBank accession no. CAC67478.1). Phylogenetically the detected penton gene of the human adenovirus 2 occupied the same clade dominated by other human mastadenoviruses (Fig. S1). Another contig of 322 nt (coverage 4.60×, GenBank accession no. MZ364305) was also part of the adenovirus penton gene and showed the highest nucleotide identity (98.8%) with another human mastadenovirus C (GenBank accession no. AHN92668.1). There were three further contigs of 435, 430 and 405 nt (coverage 6.41×, 8.19× and 5.27×, respectively; GenBank accession no MZ364303, MZ364302 and MZ364304, respectively) spanning part of the human mastadenovirus C capsid protein precursor pIIIa gene, encapsidation protein 52 K gene and capsid protein precursor pVI gene, respectively. The similarities of assembled sequences were found to be a 430/435 (98.9%) nucleotide match to the human mastadenovirus C (GenBank accession no. MH121111.1) (Table S3 and Fig. S2), a 386/390 (99.0%) nucleotide match to the human mastadenovirus C (GenBank accession no. MN088492.1) (Table S4), and a 403/403 (100%) nucleotide match to the human mastadenovirus C (GenBank accession no. MN088492.1), respectively.
Beak and feather disease virus (BFDV) is a member of the Circoviridae family and has a relatively simple but compact circular, ambisense single-stranded DNA (ssDNA) genome of approximately 2.0 kb encoding a replicase (Rep) and a single capsid protein (Cap) which facilitates whole-genome viral epidemiological analysis6,7,32,33. BFDV is highly genetically diverse and infects a large number of psittacine and non-psittacine bird species5,33,34,35,36,37,38,39. A complete genome sequence of BFDV (2250 nt in length, coverage 171.10×) was detected in this study (Fig. 2a). As shown in Fig. 2a, the genome of BFDV (GenBank accession no. MZ364298) contained two bidirectional ORFs encoding the putative Rep and Cap proteins. A BLASTP search in GenBank based on the protein sequence of Rep showed the highest identity of 98.27% (query coverage, 69%; E-value, 0.0) to a BFDV isolated from an Australian psittacine bird, the long-billed corella (Cacatua tenuirostris), in 2010 (GenBank accession no. KF385420.1)7, whereas Cap showed the highest identity of 96.76% (query coverage, 100%; E-value, 7.0 × 10−177) to a BFDV isolated from an Australian psittacine bird, the little corella (Cacatua sanguinea), in 2014 (GenBank accession no. KY189059.1). Phylogenetic analysis based on selected completed genome sequences of BFDV showed that the sequenced BFDV strain in this study clustered in a distinct subclade with other BFDVs isolated from Australian cockatoos (Fig. S3), which suggests that the BFDV sequenced in this study is likely to have originated from the same ancestor.
Two complete CRESS-DNA genomes of 1217 nt (coverage 27.40×, Fig. 2b) and 1094 nt (coverage 4.80×, Fig. 2c) in length were detected in this study. As shown in Fig. 2 (b and c), both the CRESS-DNA genomes encoded a putative Rep gene, while CRESS-DNA isolate001 also contained an additional hypothetical gene. A BLASTP search in GenBank, based on the protein sequence of Rep of CRESS-DNA isolate001, showed the highest identity of 80.94% (query coverage, 99%; E-value, 3.0 × 10−171) to a CRESS-DNA virus isolated from seabass tissue in the United States of America in 2017 (accession no. AXQ65592.1), whereas Rep of CRESS-DNA isolate002 showed the highest identity of 80.92% (query coverage, 100%; E-value, 9.0 × 10−175) to a CRESS-DNA virus isolated from seabass tissue in the United States of America in 2017 (accession no. AXQ65592.1). At the genomic level, both the CRESS-DNA genomes sequenced in this study showed highest nucleotide identity (84.27%) with a CRESS-DNA virus isolated from seabass tissue in the United States of America in 2017 (accession no. MH648902.1). Phylogenetically both the CRESS-DNA virus detected in this study positioned in a subclade with other CRESS-DNA viruses sequenced from seabass and dragonfly (Procordulia grayi)40 in the United States of America and New Zealand, respectively (Fig. S4).
The family Endornaviridae includes viruses with linear, single-stranded, positive-sense RNA genomes that range from 9.7 to 17.6 kb and have been reported as infecting plants, fungi and oomycetes41. In this study, two fragments of psittacine alphaendornavirus (contigs of 697 nt, and 613 nt in length, coverage 4.83× and 4.93×, GenBank accession no. MZ561459 and MZ561460, respectively) were detected. Both the fragments of psittacine alphaendornavirus showed the highest nucleotide identities (> 99.0%, Figs. S5 and S6) with the polyprotein of Helianthus annuus alphaendornavirus (GenBank accession no. NC_040799.1), which was isolated from sunflowers (Helianthus annuus) in China in 2016. Phylogenetically, two fragments of a psittacine alphaendornavirus detected in the Neophema birds were evolutionarily linked with the alphaendornaviruses isolated from sunflowers (Helianthus annuus) in China (Fig. S7).
Viruses of the family Picobirnaviridae are non-enveloped, bi-segmented, double-stranded RNA (dsRNA) viruses of approximately 2.4 to 2.7 kb in length, which contain two functional genes, a structural capsid gene and an RNA-dependent RNA polymerase (RdRp) gene42. Evidence of picobirnavirus (PBV) sequences have been reported in humans, invertebrates, environmental water samples and a wide range of animals, including birds worldwide12,42,43,44,45,46,47. In this study, two contigs of 444 nt and 345 nt in length (coverage 4.86× and 11.15×, respectively), were detected and were found to be psittacine picobirnaviruses. The contig of 444 nt spans part of the RdRp gene of the human picobirnavirus sp. (GenBank accession no. KU892528.1); however, it showed only a 7% query sequence coverage with identities of 100%. Another contig of 345 nt had a 100% nucleotide match (query coverage of 8%) to the capsid protein gene of the human picobirnavirus sp. (GenBank accession no. KU892525.1). However, this study was unable to detect any protein sequence similarity of the detected contigs using a BLAST search. In the resultant phylogenetic tree, the psittacine picobirnaviruses positioned in a subclade dominated by picobirnaviruses sequenced from Australian Shelduck45, chickens48 and humans (Fig. S8). However, bootstrap support for most of nodes in this phylogeny was poor and further sequencing followed by clarification of the evolutionary relationship among these viruses is required.
Picornaviruses are non-enveloped, positive-sense, single-stranded RNA viruses, which infect many vertebrates, including mammals and birds. Recent studies have detected several picorna-like viruses from Australian wild birds12. In this study, a contig of 1135 nt in length (coverage 10.17×, GenBank accession no. MZ645220) was detected, which was part of the hypothetical protein 1 gene of Hubei picorna-like virus (GenBank accession no. KX883953.1) and showed the highest nucleotide identity (94.03%) between them. This contig had a 94.6% match to the hypothetical protein 1 gene of Hubei picorna-like virus at the amino acid level, followed by a 69.7% match to the structural polyprotein of soybean thrips bicistronic virus 1 (GenBank accession no. QQN90113.1) (Fig. 3a). A phylogenetic tree, using the hypothetical protein 1 gene of Hubei picorna-like virus, showed that the picorna-like virus detected in this study occupied the same clade as Hubei picorna-like virus 51 (Fig. 3b), sequenced from flying insects13.
This study used a viral metagenomic approach to help determine viral diversity in apparently healthy captive psittacine birds housed for the ‘Parrot Genome Sequencing Project’. Many captive birds are sourced from the wild and are used for various in vivo experiments, but their virome structure is largely unknown. Accordingly, this study identified many known and unknown viruses in the families Adenoviridae, Circoviridae, Endornaviridae, Picobirnaviridae and Picornaviridae, as well as unclassified CRESS-DNA viruses, in apparently healthy Australian birds of the genus Neophema.
The main aims of this study were not only to understand the virome structure of apparently healthy Australian Neophema birds but also to track the evolutionary history of a recently sequenced PsSiAdV-F in the critically endangered Australian parrot, the orange-bellied parrot (Neophema chrysogaster)20. Notably, the PsSiAdV-F sequenced from this study was very similar (a nucleotide identity of > 99% at genome level) to the recently sequenced two PsSiAdV-F strains of OBP2209 and WVL19065-01E from the orange-bellied parrot and African grey parrot (Psittacus erithacus), respectively, and phylogenetically all the three strains of PsSiAdV-F were positioned in a monophyletic subclade (100% bootstrap support) (Fig. 1) that is basal to all known avian siadenoviruses. This suggests that all the avian siadenoviruses evolved from the ancestral bird that gave rise to PsSiAdV-F, which is consistent with a recent study that reported that PsSiAdV-F was basal to all avian siadenoviruses20. Evidence that this ancestral species may have been a primitive Australian or African parrot includes the observation that PsSiAdV-F has only been detected in psittacine species where it causes predominantly subclinical infections. The PsSiAdV-F is also widespread in captive populations of at least three species of Australian parrot: the orange-bellied parrot (Neophema chrysogaster), the scarlet-chested parrot (Neophema splendida) and the Bourke’s parrot (Neopsephotus bourkii)20,49,50. The finding that PsSiAdV-F originated in Australian psittacine birds needs to be confirmed by it being identified in other Australian psittacine species, and, to date, surveys of psittacine birds in Australia have not detected it.
Except for a recent study in which the authors reported a mastadenovirus sequence (a variant of murine adenovirus 2) in Australian wild birds species26, evidence of mammalian adenoviruses (such as human adenovirus 2 and human mastadenoviruses) in avian host species is scarce. For the first time, this study reported five different contigs, which were found to be highly similar (> 99% identity) to the human adenovirus 2 and human mastadenoviruses. It is impossible to elucidate the host–pathogen dynamics of mammalian adenoviruses from this case alone, but it is evident that captive birds have close interaction with humans. This, therefore, suggests that these birds are potentially exposed to viruses from humans and other animal species, which requires further investigation.
Circoviruses are small, non-enveloped icosahedral viruses with circular ssDNA genomes of only 2 kb, making them the smallest autonomously replicating pathogens. Numerous studies have shown that all Psittaciformes, including threatened Australian psittacine bird species, can be infected by a circovirus (e.g., beak and feather disease virus, BFDV)5,6,33,37,51,52. Therefore, the finding of the BFDV genome in this study supports the presence of BFDV as an endemic virus in Australian birds. In addition, the CRESS-DNA virus has a small circular rep-encoding ssDNA (CRESS-DNA) genome that encodes a replication-associated protein (Rep)53,54. However, there is lack of adequate information on the existence of the CRESS-DNA virus, as well as the pathology they may develop in psittacine birds. This study found evidence of two divergent CRESS-DNA virus genomes, which were evolutionarily related with the CRESS-DNA virus isolated from seabass tissue in the United States of America in 2017 (nucleotide identity ranging from 80 to 84%), and this study cannot rule out that they may also represent food contaminants and environmental pollution. In the absence of any CRESS-DNA virus sequence from Australian psittacine birds, this study suggests that more structured sampling that includes additional hosts could help better understand the ecology and epizootiology of potentially pathogenic CRESS-DNA viruses.
Birds may serve as reservoirs, intermediate hosts or transporters of parasites and pathogens11,12,55. The similarity of some of the detected sequences of the avian and non-avian host-associated RNA viruses to the viruses that were characterised in China (e.g., Endornaviridae)13, and Belgium (e.g., Picobirnaviridae), highlights that birds may harbour a diverse range of viruses that need to be investigated further. It is also noteworthy to mention that during the BLASTN and BLASTX search, none of the top-matched sequences were similar to or related to those viruses from Australia.
Using viral metagenomic or meta-transcriptomic approaches, picornaviruses were previously reported in wild birds without understanding any pathology they may cause11,12,56,57. However, information on picornaviruses harboured by apparently healthy captive birds is scarce. The picornavirus detected in this study was shown to be closely phylogenetically related to a Hubei picorna-like virus 51 sequenced from flying insects13 (Fig. 3b). While further work is required to understand the ecology and pathogenicity of the viruses detected in this study, numerous studies have suggested that the presence of immunosuppressive diseases such as BFDV may increase the likelihood that subclinical infections by other viruses progress to clinical disease in birds infected with multiple agents58,59. In addition, BFDV infection is often associated with clinical evidence of acquired immunodeficiency, leading to a variety of secondary or opportunistic infections18,60,61.
The findings of non-avian host-associated viruses detected in this study were endornavirus and picornavirus, which agree with the previous reports by others, where a large number of non-avian associated viruses were detected in Australian birds11,12. Some aspects of this case are difficult to explain fully without conducting ethically debatable experimental virus-transmission experiments. A possible scenario for transmission could be ingestion of insect vector contaminated food. Alternatively, they may have been eaten infected insects before entering into the pet- trade. Moreover, this may also indicate that many diverse viruses may be harboured by the closely related avian host species, not detected yet.
Overall, this study demonstrates a significant viral diversity in Neophema birds. Notably, it shows the evidence of genetically and phylogenetically similar psittacine siadenovirus (PsSiAdV-F) circulating between critically endangered orange-bellied parrots and other Neophema bird species. This study demonstrates that apparently healthy captive birds can sustain a significant viral diversity. Furthermore, the presence of viral pathogens in these birds presents a concerning example of the ease with which such infectious agents may spread to other captive birds, and the humans in their proximity. Finally, the findings of this study emphasise the importance of testing captive birds for possible pathogens to prevent the spread of potentially deadly diseases.
Sampling and ethical approval
In 2020, four fresh faecal samples were collected from two different species, the elegant parrot (Neophema elegans) and the scarlet-chested parrot (Neophema splendida) that were housed in La Trobe Animal Research and Teaching Facility for the ‘Parrot Genome Sequencing Project’. These samples were stored at − 80 °C within 1 h of their collection and kept in these conditions until their processing. Bird sampling was obtained following approved guidelines set by the Australian Code of Practice for the Care and Use of Animals for Scientific Purposes and approved by the La Trobe University Animal Ethics Committee (research permit number AEC19035) and Department of Environment, Land, Water and Planning (permit number 10009300).
Virus enrichment and virus nucleic acid extraction
Elimination of potential impurities, such as host cells, bacteria, food particles and free nucleic acids, from the faecal samples, followed by enrichment of virus particles was performed as per stated methods12, with miror variations. Briefly, the faecal materials were aseptically resuspended and homogenised vigorously in sterile phosphate-buffered saline (PBS) (1:10) and centrifuged at 2500 × g for 90 min at 4 °C. The supernatant was filtered using a 0.80 μm syringe filter and the filtrate processed downstream. The samples were then ultracentrifuged at 178,000 g for 1 h (30 psi for 1 h) at 4 °C using the Hitachi Ultracentrifuge CP100NX. The supernatant was discarded, and the pellet was suspended in 130 µL of sterile PBS. The filtrates were then nuclease-treated using 2 µL of benzonase nuclease (25–29 U/µL, purity > 90%, Millipore) and 1 µL of micrococcal nuclease (2,000,000 gel units/mL, New England Biolabs) and incubated at 37 °C for 2 h. The nuclease reaction was stopped by adding 3 µL of 500 mM ethylenediaminetetraacetic acid (EDTA). Viral nucleic acids were extracted using the QIAamp Viral RNA Mini Kit (Qiagen, Valencia, CA, USA), without adding any carrier RNA, which allowed the extraction of both viral DNA and RNA simultaneously. The quantity and quality of the isolated nucleic acids were determined using Nanodrop and an Agilent Tape Station (Agilent Technologies, Mulgrave, VIC, Australia) by the Genomic Platform, La Trobe University.
Before library construction, extracted nucleic acids were subjected to cDNA synthesis, and amplification was carried out using the Whole Transcriptome Amplification Kit (WTA2, Sigma-Aldrich, Darmstadt, Germany) as per manufacturer instructions. Amplified PCR products were then purified using the Wizard® SV Gel and PCR Clean-Up kit (Promega, Madison, WI, USA). The quantity and quality of the purified product were checked using a Qubit dsDNA high sensitivity assay kit with Qubit Fluorometer v3.0 (Thermo Fisher Scientific, Waltham, MA, USA).
The library construction was performed as a pool that contained four samples using the Illumina DNA Prep (Illumina, San Diego, CA, USA) as per kit instructions, starting with 250 ng of DNA as measured by Qubit (Invitrogen). The quality and quantity of the prepared library was assessed by the Australian Genome Research Facility (AGRF), Melbourne, Australia. The prepared library was normalised and pooled in equimolar quantities. The quality and quantity of the final pooled library were further assessed as described above before sequencing by the facility. According to the manufacturer’s instructions, cluster generation and sequencing of the pooled library was performed with read lengths of 150-bp paired-end on Illumina® NovaSeq chemistry.
The resulting 52.1 million raw sequence reads were used for a quality control check using CLC Genomics Workbench (version 9.5.4). Preliminary quality evaluation for all raw reads was generated, pre-processed to remove ambiguous base calls and poor-quality reads and trimmed to remove the Illumina adapter sequences. Trimmed sequence reads were mapped against the chicken genome (Gallus, GenBank accession number NC_006088) to remove likely host DNA contamination. In addition, reads were further mapped to Escherichia coli bacterial genomic sequence (GenBank accession no. U00096) to remove possible bacterial contamination. Unmapped reads were used as input data for de novo assembly using a SPAdes assembler (version 3.10.1)62 under the ‘careful’ parameter in the LIMS-HPC system (a High-Performance Computer specialised for genomics research in La Trobe University). The resulting contigs were compared against the nonredundant nucleotide and protein databases on GenBank using BLASTN and BLASTX63, respectively, with an E-value threshold of 1 × 10−5 to remove potential false positives. Contigs that were significant BLAST hits with bacteria, eukaryotes or fungi were filtered out to remove non-viral reads. Virus contigs of interest greater than 300 nucleotides (nt) were imported in Geneious software (Biomatters Ltd., New Zealand, version 10.2.2) for further functional analysis. Average coverage of the viral contigs were calculated using the clean raw reads in CLC Genomics Workbench (version 9.5.4).
The detected complete genome of adenovirus was annotated as per stated protocol23 using Geneious software (version 10.2.2, Biomatters, New Zealand). Briefly, open reading frames (ORFs) longer than 30 amino acids, with a methionine start codon (ATG) and minimal overlap with other ORFs (not exceeding 50% of one of the genes), were selected and annotated. All the detected circovirus genomes were annotated using Geneious software (version 10.2.2, Biomatters, New Zealand), where representative circovirus genomes were used as reference guidelines6,7,38,64. Similarity BLAST searches were performed on the predicted ORFs and were annotated as potential genes if predicted ORFs showed significant sequence similarity to known viral or cellular genes (E-value threshold of 1 × 10−5)63. Detected partial genomes or genes content were annotated according to the BLASTN and BLASTX63 search results.
Comparative genomics and phylogenetic analyses
Genomic features of the newly sequenced viral genomes were visualised using Geneious (version 10.2.2). Sequence similarity percentages between representative viruses were determined using tools available in Geneious (version 10.2.2).
For phylogenetic analyses, representative viral genome or gene sequences were downloaded from GenBank, and virus-specific trees were constructed using CLC Genomics Workbench (version 9.5.4) and Geneious software (version 10.2.2, Biomatters, New Zealand). Amino acid sequences of protein-coding genes and nucleotide sequences of the selected partial genes were aligned using the MAFTT L-INS-I algorithm implemented in Geneious (version 7.388)65. To determine the best-fit model to construct phylogenetic analyses, a model test was performed using CLC Genomics Workbench (version 9.5.4) using default parameters, favouring a general-time-reversible model with gamma distribution rate variation and a proportion of invariable sites (GTR + G + I). Phylogenetic analyses for nucleotide and protein sequences were performed using the GTR and WAG substitution model, respectively, with 1000 bootstrap support in CLC Genomics Workbench (version 9.5.4).
All sequences analysed have been deposited in NCBI GenBank under the accession numbers MZ364296, MZ364298-MZ364305, MZ561459-MZ561460, MZ645220. Raw sequencing data from this study has been deposited in the NCBI Sequence Read Achieve (SRA) under the BioProject ID: PRJNA750905 (BioSample accessions: SAMN20500864) (http://www.ncbi.nlm.nih.gov/sra/).
Lupiani, B. & Reddy, S. M. The history of avian influenza. Comp. Immunol. Microbiol. Infect. Dis. 32, 311–323. https://doi.org/10.1016/j.cimid.2008.01.004 (2009).
Alexander, D. J. Newcastle disease and other avian paramyxoviruses. OIE Rev. Sci. Tech. 19, 443–462. https://doi.org/10.20506/rst.19.2.1231 (2000).
Fabricant, J. The early history of infectious bronchitis. Avian Dis. 42, 648–650. https://doi.org/10.2307/1592697 (1998).
Allison, A. B. et al. Cyclic avian mass mortality in the northeastern United States is associated with a novel orthomyxovirus. J. Virol. 89, 1389–1403. https://doi.org/10.1128/jvi.02019-14 (2015).
Raidal, S. R., Sarker, S. & Peters, A. Review of psittacine beak and feather disease and its effect on Australian endangered species. Aust. Vet. J. 93, 466–470. https://doi.org/10.1111/avj.12388 (2015).
Sarker, S. et al. Mutability dynamics of an emergent single stranded DNA virus in a naïve host. PLoS ONE 9, e85370. https://doi.org/10.1371/journal.pone.0085370 (2014).
Sarker, S. et al. Phylogeny of beak and feather disease virus in cockatoos demonstrates host generalism and multiple-variant infections within Psittaciformes. Virology 460–461, 72–82. https://doi.org/10.1016/j.virol.2014.04.021 (2014).
Chang, W. S. et al. Metatranscriptomic analysis of virus diversity in urban wild birds with paretic disease. J. Virol. https://doi.org/10.1128/jvi.00606-20 (2020).
Chan, D. T. C., Poon, E. S. K., Wong, A. T. C. & Sin, S. Y. W. Global trade in parrots—Influential factors of trade and implications for conservation. Glob. Ecol. Conserv. 30, e01784. https://doi.org/10.1016/j.gecco.2021.e01784 (2021).
Shi, M., Zhang, Y. Z. & Holmes, E. C. Meta-transcriptomics and the evolutionary biology of RNA viruses. Virus Res. 243, 83–90. https://doi.org/10.1016/j.virusres.2017.10.016 (2018).
Vibin, J., Chamings, A., Klaassen, M. & Alexandersen, S. Metagenomic characterisation of additional and novel avian viruses from Australian wild ducks. Sci. Rep. 10, 22284. https://doi.org/10.1038/s41598-020-79413-9 (2020).
Vibin, J. et al. Metagenomics detection and characterisation of viruses in faecal samples from Australian wild birds. Sci. Rep. 8, 8686. https://doi.org/10.1038/s41598-018-26851-1 (2018).
Shi, M. et al. Redefining the invertebrate RNA virosphere. Nature 540, 539–543. https://doi.org/10.1038/nature20167 (2016).
Temmam, S. et al. Characterization of viral communities of biting midges and identification of novel Thogotovirus species and Rhabdovirus genus. Viruses 8, 77. https://doi.org/10.3390/v8030077 (2016).
Lim, E. S. et al. Early life dynamics of the human gut virome and bacterial microbiome in infants. Nat. Med. 21, 1228–1234. https://doi.org/10.1038/nm.3950 (2015).
Zablocki, O. et al. High-level diversity of tailed phages, eukaryote-associated viruses, and virophage-like elements in the metaviromes of antarctic soils. Appl. Environ. Microbiol. 80, 6888–6897. https://doi.org/10.1128/aem.01525-14 (2014).
Roux, S. et al. Ecogenomics of virophages and their giant virus hosts assessed through time series metagenomics. Nat. Commun. 8, 858. https://doi.org/10.1038/s41467-017-01086-2 (2017).
Sutherland, M. et al. Disease surveillance in wild Victorian cacatuids reveals co-infection with multiple agents and detection of novel avian viruses. Vet. Microbiol. 235, 257–264. https://doi.org/10.1016/j.vetmic.2019.07.012 (2019).
Amery-Gale, J. et al. Avian viral surveillance in Victoria, Australia, and detection of two novel avian herpesviruses. PLoS ONE 13, e0194457. https://doi.org/10.1371/journal.pone.0194457 (2018).
Athukorala, A. et al. Genomic characterisation of a highly divergent siadenovirus (Psittacine siadenovirus F) from the critically endangered orange-bellied parrot (Neophema chrysogaster). Viruses https://doi.org/10.3390/v13091714 (2021).
Ahi, Y. S. & Mittal, S. K. Components of adenovirus genome packaging. Front. Microbiol. 7, 1503–1503. https://doi.org/10.3389/fmicb.2016.01503 (2016).
Tan, B. et al. Novel bat adenoviruses with low G+C content shed new light on the evolution of adenoviruses. J. Gen. Virol. 98, 739–748. https://doi.org/10.1099/jgv.0.000739 (2017).
Athukorala, A., Forwood, J. K., Phalen, D. N. & Sarker, S. Molecular characterisation of a novel and highly divergent Passerine adenovirus 1. Viruses https://doi.org/10.3390/v12091036 (2020).
Doszpoly, A., Harrach, B., LaPatra, S. & Benkő, M. Unconventional gene arrangement and content revealed by full genome analysis of the white sturgeon adenovirus, the single member of the genus Ichtadenovirus. Infect. Genet. Evol. 75, 103976. https://doi.org/10.1016/j.meegid.2019.103976 (2019).
Harrach, B., Tarján, Z. L. & Benkő, M. Adenoviruses across the animal kingdom: A walk in the zoo. FEBS Lett. 593, 3660–3673. https://doi.org/10.1002/1873-3468.13687 (2019).
Vaz, F. F. et al. Opportunistic sampling of wild native and invasive birds reveals a rich diversity of adenoviruses in Australia. Virus Evol. https://doi.org/10.1093/ve/veaa024 (2020).
Lange, C. E. et al. First evidence of a new simian adenovirus clustering with human mastadenovirus F viruses. Virol. J. 16, 147. https://doi.org/10.1186/s12985-019-1248-z (2019).
Kovács, E. R. & Benkő, M. Complete sequence of raptor adenovirus 1 confirms the characteristic genome organization of siadenoviruses. Infect. Genet. Evol. 11, 1058–1065. https://doi.org/10.1016/j.meegid.2011.03.021 (2011).
Kaján, G. L., Doszpoly, A., Tarján, Z. L., Vidovszky, M. Z. & Papp, T. Virus–host coevolution with a focus on animal and human DNA Viruses. J. Mol. Evol. 88, 41–56. https://doi.org/10.1007/s00239-019-09913-4 (2020).
Binder, A. M. et al. Human adenovirus surveillance—United States, 2003–2016. MMWR Morb. Mortal. Wkly. Rep. 66, 1039–1042. https://doi.org/10.15585/mmwr.mm6639a2external (2017).
Li, Y. et al. Host range, prevalence, and genetic diversity of adenoviruses in bats. J. Virol. 84, 3889–3897. https://doi.org/10.1128/jvi.02497-09 (2010).
Sarker, S. et al. Structural insights into the assembly and regulation of distinct viral capsid complexes. Nat. Commun. 7, 13014. https://doi.org/10.1038/ncomms13014 (2016).
Sarker, S., Forwood, J. K., Ghorashi, S. A., Peters, A. & Raidal, S. R. Beak and feather disease virus genotypes in Australian parrots reveal flexible host-switching. Aust. Vet. J. 93, 471–475. https://doi.org/10.1111/avj.12389 (2015).
Sarker, S., Forwood, J. K. & Raidal, S. R. Beak and feather disease virus: biology and resultant disease. Wiki J. Sci. 3, 7. https://doi.org/10.15347/wjs/2020.007 (2020).
Sarker, S., Lloyd, C., Forwood, J. & Raidal, S. R. Forensic genetic evidence of beak and feather disease virus infection in a Powerful Owl, Ninox strenua. Emu 116, 71–74. https://doi.org/10.1071/MU15063 (2016).
Fogell, D. J., Martin, R. O. & Groombridge, J. J. Beak and feather disease virus in wild and captive parrots: An analysis of geographic and taxonomic distribution and methodological trends. Arch. Virol. https://doi.org/10.1007/s00705-016-2871-2 (2016).
Das, S. et al. Evolution of circoviruses in lorikeets lags behind its hosts. Mol. Phylogenet. Evol. 100, 281–291. https://doi.org/10.1016/j.ympev.2016.04.024 (2016).
Sarker, S. et al. Evidence of a deep viral host switch event with beak and feather disease virus infection in rainbow bee-eaters (Merops ornatus). Sci. Rep. 5, 14511. https://doi.org/10.1038/srep14511 (2015).
Sarker, S., Das, S., Ghorashi, S. A., Forwood, J. K. & Raidal, S. R. Molecular characterization of genome sequences of beak and feather disease virus from the Australian twenty-eight parrot (Barnardius zonarius semitorquatus). Genome Announc. 2, e01214–e01255. https://doi.org/10.1128/genomeA.01255-14 (2014).
Dayaram, A. et al. Diverse circular replication-associated protein encoding viruses circulating in invertebrates within a lake ecosystem. Infect. Genet. Evol. 39, 304–316. https://doi.org/10.1016/j.meegid.2016.02.011 (2016).
Valverde, R. A. et al. ICTV virus taxonomy profile: Endornaviridae. J. Gen. Virol. 100, 1204–1205. https://doi.org/10.1099/jgv.0.001277 (2019).
Delmas, B. et al. ICTV virus taxonomy profile: Picobirnaviridae. J. Gen. Virol. 100, 133–134. https://doi.org/10.1099/jgv.0.001186 (2019).
Malik, Y. S. et al. Epidemiology, phylogeny, and evolution of emerging enteric Picobirnaviruses of animal origin and their relationship to human strains. Biomed. Res. Int. 2014, 780752. https://doi.org/10.1155/2014/780752 (2014).
Ganesh, B., Masachessi, G. & Mladenova, Z. Animal picobirnavirus. Virusdisease 25, 223–238. https://doi.org/10.1007/s13337-014-0207-y (2014).
Wille, M. et al. Virus-virus interactions and host ecology are associated with RNA virome structure in wild birds. Mol. Ecol. 27, 5263–5278. https://doi.org/10.1111/mec.14918 (2018).
Wille, M., Shi, M., Klaassen, M., Hurt, A. C. & Holmes, E. C. Virome heterogeneity and connectivity in waterfowl and shorebird communities. ISME J. 13, 2603–2616. https://doi.org/10.1038/s41396-019-0458-0 (2019).
Huaman, J. L. et al. Molecular epidemiology and characterization of picobirnavirus in wild deer and cattle from Australia: Evidence of genogroup I and II in the upper respiratory Tract. Viruses 13, 1492 (2021).
Silva, R. R. et al. Genogroup I avian picobirnavirus detected in Brazilian broiler chickens: A molecular epidemiology study. J. Gen. Virol. 95, 117–122. https://doi.org/10.1099/vir.0.054783-0 (2014).
Yang, N. et al. Psittacid Adenovirus-2 infection in the critically endangered orange-bellied parrot (Neophema chrysogastor): A key threatening process or an example of a host-adapted virus?. PLoS ONE 14, e0208674. https://doi.org/10.1371/journal.pone.0208674 (2019).
Phalen, D. N. et al. A survey of a mixed species aviary provides new insights into the pathogenicity, diversity, evolution, host range, and distribution of psittacine and passerine adenoviruses. Avian Pathol. 48, 437–443. https://doi.org/10.1080/03079457.2019.1617835 (2019).
Raidal, S. R. & Peters, A. Psittacine beak and feather disease: Ecology and implications for conservation. Emu 118, 80–93. https://doi.org/10.1080/01584197.2017.1387029 (2018).
Eastwood, J. R. et al. Phylogenetic analysis of beak and feather disease virus across a host ring-species complex. Proc. Natl. Acad. Sci. U. S. A. 111, 14153–14158. https://doi.org/10.1073/pnas.1403255111 (2014).
Breitbart, M. et al. ICTV virus taxonomy profile: Circoviridae. J. Gen. Virol. 98, 1997–1998. https://doi.org/10.1099/jgv.0.000871 (2017).
Kaszab, E. et al. Occurrence and genetic diversity of CRESS DNA viruses in wild birds: A Hungarian study. Sci. Rep. 10, 7036. https://doi.org/10.1038/s41598-020-63795-x (2020).
Hubálek, Z. An annotated checklist of pathogenic microorganisms associated with migratory birds. J. Wildl. Dis. 40, 639–659. https://doi.org/10.7589/0090-3558-40.4.639 (2004).
Wille, M. et al. Virus–virus interactions and host ecology are associated with RNA virome structure in wild birds. Mol. Ecol. 27, 5263–5278. https://doi.org/10.1111/mec.14918 (2018).
Wille, M., Shi, M., Hurt, A. C., Klaassen, M. & Holmes, E. C. RNA virome abundance and diversity is associated with host age in a bird species. Virology 561, 98–106. https://doi.org/10.1016/j.virol.2021.06.007 (2021).
Ramis, A. et al. Diagnosis of psittacine beak and feather disease (PBFD) viral infection, avian polyomavirus infection, adenovirus infection and herpesvirus infection in psittacine tissues using DNA in situ hybridization. Avian Pathol. 23, 643–657. https://doi.org/10.1080/03079459408419034 (1994).
Raidal, S. R., Cross, G. M., Tomaszewski, E., Graham, D. L. & Phalen, D. N. A serologic survey for avian polyomavirus and Pacheco’s disease virus in Australian cockatoos. Avian Pathol. 27, 263–268. https://doi.org/10.1080/03079459808419334 (1998).
Jacobson, E. R. et al. Feather and beak dystrophy and necrosis in cockatoos: Clinicopathologic evaluations. J. Am. Vet. Med. Assoc. 189, 999–1005 (1986).
Latimer, K. S., Niagro, F. D., Steffens, W. L., Ritchie, B. W. & Campagnoli, R. P. Polyomavirus encephalopathy in a Ducorps’ cockatoo (Cacatua ducorpsii) with psittacine beak and feather disease. J. Vet. Diagn. Invest. 8, 291–295. https://doi.org/10.1177/104063879600800303 (1996).
Bankevich, A. et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. https://doi.org/10.1089/cmb.2012.0021 (2012).
Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403–410. https://doi.org/10.1016/S0022-2836(05)80360-2 (1990).
Liu, Q. et al. Viral metagenomics revealed diverse CRESS-DNA virus genomes in faeces of forest musk deer. Virol. J. 17, 61. https://doi.org/10.1186/s12985-020-01332-y (2020).
Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. https://doi.org/10.1093/molbev/mst010 (2013).
Subir Sarker is the recipient of an Australian Research Council Discovery Early Career Researcher Award (Grant Number DE200100367) funded by the Australian Government. The author like to acknowledge Dr Anthony Chamings, Geelong Centre for Emerging Infectious Diseases, Geelong, Victoria for his valuable suggestion on methodology.
The author declares no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Sarker, S. Metagenomic detection and characterisation of multiple viruses in apparently healthy Australian Neophema birds. Sci Rep 11, 20915 (2021). https://doi.org/10.1038/s41598-021-00440-1