Skip to main content

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

Abstract

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.

Background

There is a growing body of research into the ecology of infectious diseases of wildlife. Despite this, some host taxa, geographic areas, and pathogens are routinely neglected in wildlife, particularly those that do not threaten human or livestock health [1]. Similarly, there has been insufficient investigation of “biodiversity diseases” where there are no perceived direct agriculture or human health links. These include such high profile examples as chytrid fungus infection in amphibians and white nose syndrome in bats [2]. Bias over the nature of emerging disease in wildlife, particularly before a diagnosis has been established, can rob us of the opportunity to conduct science-based risk assessments into potential collective health threats and impacts. Wildlife diseases may have implications for the loss of natural biomass and important ecosystem services such as clean water and pollination. Wildlife diseases are also central to a One Health perspective because they can threaten native species with extinction [3], can be used as biological controls against invasive pests, and may act as conduits for viruses to move to domestic species and ultimately humans [4].

Although common and widely distributed across Australia, the brushtail possum (Trichosurus vulpecula) is protected under conservation legislation. In contrast, it is a highly successful invasive vertebrate pest in New Zealand, where it was introduced from Tasmania in the 1830s to build a local fur industry [5]. The species then rapidly adapted to the natural environment of New Zealand, causing devastating destruction of native forests and wildlife. Possums also function as an important reservoir of the bacteria Mycobacterium bovis [6] and Leptospira sp. [7, 8]. Brushtail possums in New Zealand are estimated to cause $NZ35M/year in agricultural losses and the shared costs of attempted control exceed $100 M/year (https://predatorfreenz.org/resources/introduced-predator-facts/possum-facts/).

Wobbly possum disease (WPD) is a severe neurological disease of Australian brushtail possums in New Zealand, caused by a virus of family Arteriviridae (genus Dipartevirus, order Nidovirales) termed wobbly possum disease virus (WPDV) (https://predatorfreenz.org/resources/introduced-predator-facts/possum-facts/). The genome of WPDV is similar to those of other arteriviruses, with the typical gene organisation of ORF1a, ORF1b, glycoproteins 2,3,4,5, membrane protein and nucleocapsid protein. WPDV was first recognized in New Zealand in diseased possums present at a research facility in 1995, and has subsequently been found in a free-living population in the same country [9, 10]. In New Zealand, WPD presents as severe, multifocal neurological signs, including incoordination, loss of balance, head tilt, circling, difficulty in climbing, abnormal gait, aimless wandering and can include central blindness. Most animals are found thin or emaciated and are either anorexic or exhibit abnormal behaviours such as daytime feeding [10,11,12]. There are no known treatments for WPD and affected animals presented to wildlife carers most often die or are euthanized. Histologically, cases of WPD in New Zealand are similar to those from experimentally induced infection with WPDV and include non-suppurative meningoencephalitis, perivascular inflammation in multiple organs, most consistently liver kidney, and salivary gland [13]. While considerable data is available on WPD in New Zealand [9, 10, 13, 14], WPD in brushtail possums in Australia has been sparsely described and investigated [15, 16].

Serological surveys of archival cases in New Zealand suggested an estimated WPD prevalence of 30% in the free-living possum population [17]. Sporadic cases and rare outbreaks of WPD have been reported in mainland Australia and Tasmania based on similar clinical presentations, histopathologic changes and exclusion of differential diagnoses such as trauma, bacterial infection, haemorrhage from rodenticide ingestion, Angiostrongylus cantonensis or Toxoplasma gondii infection. Previous investigations diagnosed WPD in 21 of 31(68%) brushtail possums with neurological dysfunction in the Sydney basin, New South Wales, between October 1998 and June 2010 [16]. Although a viral aetiology for mainland Australian WPD has been suspected, a causative agent has not been established and there is considerable uncertainty surrounding the microbial causes, potential transmission routes, carrier states, and other elements of host-virus ecology.

Arteriviruses tend to be host species-specific and cause a number of distinct disease syndromes with unique pathogenesis and epidemiology. Notable pathogens in this family include equine arteritis virus, porcine reproductive and respiratory syndrome virus, and simian haemorrhagic fever virus [18]. These viruses are phylogenetically distinct from the recently emerging Nidovirales from the family Coronaviridae that have been identified as pathogens of Australian reptiles, causing respiratory disease in shingleback lizards [19] and a mass mortality in a species of freshwater turtle associated with Bellinger River virus [20]. Phylogenetically, WPDV falls into what appears to be a divergent branch lineage of the Arteriviridae [21]. Current hypotheses for the origin of WPDV in New Zealand are that the virus evolved and emerged in possums through cross-species transmission after translocation, or the virus is possum-specific and was carried by animals during their introduction from Australia [5]. The identification of WPDV in native animals on the Australian mainland would be important evidence supporting the latter.

Hepaciviruses are a genus of positive-sense, single-stranded RNA viruses from the family Flaviviridae, the best characterised of which is Hepatitis C virus (HCV) due to its association with hepatitis and hepatocellular carcinoma in humans worldwide [22]. However, other hepaciviruses have been identified in dogs [23], horses [24], rodents [25], bats [26] and ducks [27], and new hepaciviruses are regularly described in a variety of animal species [22]. For example, a novel hepacivirus was recently identified in bald eagles associated with mass mortality events in the U.S. [28]. To date, however, hepaciviruses have not been reported from any marsupial species. Indeed, despite their ubiquity and proximity to human populations, little is known about the possum virome, with previous descriptions only covering enteroviruses [29], an adenovirus [30], and a papillomavirus [31].

Herein, we review archived Australian cases of WPD to build a more robust syndrome description and, where frozen tissues were available, applied total RNA sequencing (“meta-transcriptomics”) to determine the infectious aetiology and potential origins of WPD in mainland Australian possums, particularly in comparison with syndromes and agents previously identified in animals from New Zealand. Accordingly, this study provides new insights into the diversity and origins of a novel group of viral pathogens.

Methods

Sample collection

Samples were collected between May 1999 and June 2010 from nine brushtail possums found within the greater Sydney basin. All animals were handled under a series of NSW Office of Environment and Heritage Licences to Rehabilitate Injured, Sick or Orphaned Protected Wildlife (#MWL000100542). Samples were collected from affected possums post mortem, immediately after euthanasia via intravenous barbiturates delivered while the animals were under gaseous anaesthetic (2% isofluorane in 1 L/min Oxygen - Isofluorane 100%, Zoetis, Australia) for veterinary examination. Fresh brain, liver and kidney were collected aseptically and frozen at -80 °C. A range of tissues were fixed in 10% neutral buffered formalin, processed in ethanol, embedded with paraffin blocks, sectioned, stained with haematoxylin and eosin and permanently mounted with a cover slip. Samples were collected under the Opportunistic Sample Collection Program of the Taronga Animal Ethics Committee, and under scientific licences #SL10469 and SL100104 issued by the NSW Office of Environment and Heritage.

Historical case review

All brushtail possum cases within the Australian Registry of Wildlife Health since systematic record keeping began in 1981 were evaluated for signalment, clinical signs and histological lesions. Animals with a combination of visual impairment, severe depression or central nervous system disturbance in conjunction with non-suppurative lesions in the central nervous system, optic tracts or eyes were considered to fit the WPD syndrome description. The severity of non-suppurative lesions in the meninges, cerebral cortex, brainstem, cerebellum, liver, kidney and eyes were graded on a scale of 0–4 where 0 represented no discernible lesions and 4 represented severe and extensive lymphoplasmacytic inflammation. Necrosis was graded on a scale of 0–4 where 1 represented mild, multifocal single cell necrosis and 4 characterised extensive malacia. Retinal atrophy was similarly assessed on a scale of 0–4 with grade 4 signifying a nearly aceullular ganglion cell layer.

Pathogen discovery using meta-transcriptomics

For RNA extraction, archival tissues including brain, liver and kidney samples of animals were processed using the RNeasy Plus Mini Kit (Qiagen, Germany). RNA concentration and integrity were measured using a NanoDrop spectrophotometer (ThermoFisher Scientific, USA) and TapeStation (Agilent). Samples were then pooled in equal proportions based on tissue type and/or individual cases for different purposes (Table 1). Illumina TruSeq stranded RNA libraries were prepared on the pooled samples following rRNA depletion using the RiboZero Gold kit (Epidemiology). Paired-end (100 bp) sequencing of the rRNA-depleted RNA libraries was then performed on an Illumina HiSeq 2500 system at the Australian Genome Research Facility (AGRF), Melbourne.

Table 1 Sumsmary of the meta-transcriptomic sequence data generated here

The meta-transcriptomic analytical pipeline was constructed based on the methods previously used in our group [32,33,34]. Accordingly, RNA sequencing reads were trimmed of low quality bases and adapter sequences and de novo assembled using Trinity 2.1.1 [35]. Assembled sequence contigs were annotated using both nucleotide and protein BLAST searches against the NCBI non-redundant sequence database. To identify low abundance organisms, the sequence reads were also annotated directly using a Diamond Blastx search against the NCBI virus RefSeq viral protein database (with an e-value cutoff of < 10− 5). Open reading frames were predicted from the viral contigs in Geneious v11.1.2 [36] with gene annotation and functional predictions made against the Conserved domain databases (CDD) [37] and Geneious. Virus read abundance was assessed using a read mapping approach available in the BBmap program [38].

PCR assays and genome sequencing of WPDV and possum hepacivirus

Our metagenomic analysis suggested the presence of both WPDV and a novel hepacivirus (see Results). To confirm their presence in each archival tissue from both diseased or healthy possums, the initial PCR primers were designed based on aligned sequence reads generated from RNA-seq libraries (Table 1). Subsequently, SuperScript IV VILO cDNA synthesis system (Invitrogen) was used to reverse transcribe the RNA from individual cases. The cDNA generated from the sampled tissues was used for viral specific PCRs targeting regions identified by RNA-Seq. All PCR assays were performed using Platinum SuperFi DNA polymerase (Invitrogen) with a final concentration of 0.2 μM for both forward and reverse primers.

Long, overlapping RT-PCR assays were also developed to complete the genome of two representative Australian WPDVs (Additional file 1: Table S1). Similarly, additional sets of PCR primers for hepacivirus genome sequencing were designed from RNA-sequencing data of individual affected cases (Table 1). All PCRs were performed from VILO cDNA prepared for initial confirmatory assays using Platinum SuperFi (Invitrogen) as per the manufacturer’s instructions. Amplicons were visualized using agarose gel electrophoresis and confirmed by both Sanger and high throughput sequencing (Nextera XT library preparation on the Illumina MiSeq platform – v2 300 cycle kit), with all procedures performed following the manufacturer’s instructions.

Wobbly possum disease virus qRT-PCR

A quantitative RT-PCR assay targeting the conserved RdRp domain was developed to identify viruses of all three different lineages of WPDV (New Zealand, Australia 1 and 2). The qRT-PCR primers and probe sequences (Integrated DNA Technologies, USA) were designed as described in SI Table 1. PrimeTime Gene Expression Mastermix (Integrated DNA Technologies, USA) was used to perform the assays. For absolute quantification, the pooled genome amplicons were prepared as standards across aerial 10-fold dilution series of known quantities and used to determine virus copy number.

Phylogenetic analysis

To determine the evolutionary relationships of WDPV and the novel hepacivirus identified in this study, we performed phylogenetic analyses of amino acid sequences of the NSP2 protein (the RNA-dependent RNA polymerase (RdRp) of arteriviruses) and the complete polyprotein of hepacivirus, respectively. All sequences were aligned using MAFFT version 7 [39] with the L-INS-i algorithm, with all ambiguously aligned regions removed using TrimAL (v.1.4.1) [40]. Maximum likelihood trees of both data sets were then estimated using IQ-TREE 1.6.7 [41] utilising the best-fit model of amino acid substitution (LG + F + I + Γ4). Statistical support for individual nodes was assessed using a bootstrap approach with 1000 replicates, and all trees were midpoint rooted. Finally, phylogenetic trees were visualized using FigTree v1.4.3 (http://tree.bio.ed.ac.uk/software/figtree, latest download date (March. 2018)).

Nucleotide sequence accession numbers

The RNA sequencing data generated in this study have been deposited in the GenBank Sequence Read Archive under accession numbers PRJNA587034. All consensus genome sequences of identified viruses have been uploaded in GenBank under accession numbers MN635447 to MN635449.

Results

Clinical and histological description of Australian wobbly possum disease virus cases

Signalment, gross and histopathology data were collated from 474 brushtail possums, 49 of which fit the syndrome description for WPD (see Additional file 3: Video). A table chronologically summarising the clinical and histological findings of WPD in brushtail possums in Australia is included as Additional file 2: Table S2. The index case of WPD on mainland Australia was a juvenile male possum found in Mosman, NSW in December 1983, while the index case in Tasmania was November 1985 based on Registry records. The majority of WPD cases occurred in adult female possums (35 female, 9 male, 6 unknown sex, 42 adult, 3 juvenile, 4 unknown age).

Blindness was the most common presenting sign (n = 35) and was characterised by dilated and unresponsive pupils and lack of a menace reflex. Depression or docility were noted in six animals. A further 4 animals were found moribund. Clinical signs attributable to the nervous system included ataxia (n = 9), circling (n = 4), strabismus (n = 1), paralysis (n = 2), nystagmus (n = 1), and knuckling (n = 1), and these animals were most often concurrently blind. Some WPD affected animals were maintained in care for many weeks, with stable or slowly progressive signs. Ophthalmological examination of blind possums often revealed an optic disc that was pale and lacked the normal fundus vascular tuft. The body condition of affected possums ranged from excellent to emaciated. The most common cause of death was euthanasia.

Common microscopic lesions in WPD affected possums are illustrated in Fig. 1 and included mild to severe non-suppurative leptomeningitis, multifocally extending into the perivascular spaces of the brainstem, cerebrum and cerebellum. Cellular infiltrates ranged from 1 to 12 cells deep, and were composed of lymphocytes, plasma cells and smaller numbers of macrophages. Small foci of gliosis were multifocally evident in the neuropil of severely affected animals. Neutrophils were seen among the inflammatory infiltrates in only two animals. Necrosis was a less consistent finding and ranged between mild, multifocal single cell necrosis of nerve cell bodies, particularly affecting cerebellar Purkinje cells, to foci of malacia or cerebellar folial atrophy. Wallerian degeneration and spongiotic change were observed in the optic tracts of the brain, and optic nerve, often in association with non-suppurative inflammation in the optic nerve, perineurium, retrobulbar fat, and sclera. Ocular lesions included mononuclear cell infiltrates, as above, evident variably within the choroid, retina, ciliary body and iris. Retinal atrophy was evident as loss of cellularity within the ganglion cell and nuclear layers, and occasionally a loss of rods and cones. Lymphoplasmacytic infiltrates in renal and hepatic tissue were only evident within animals originating from Tasmania, except in rare cases where concurrent disease could account for lesions.

Fig. 1
figure 1

Histopathologic findings in Australian wobbly possum disease cases (hematoxylin and eosin). a Focal gliosis (black asterix), non-suppurative cellular infiltrates distending the leptomeninges (white asterix), and surrounding blood vessels (black arrowhead) within the brainstem. b Lymphocytes, plasma cells and scattered macrophages around a neural blood vessel. c Non-suppurative infiltrates in retrobulbar adipose tissue (black arrow), the perineurium of the optic nerve (black arrow head), and the scleral perivascular space (white arrow head). Retinal atrophy with an acellular ganglion cell layer (black asterix and inset). d Distention of the choroid layer with non-suppurative inflammation (arrow). e Optic nerve with Wallerian degeneration, illustrated by a macrophage within an axonal chamber (black arrowhead), and a non-suppurative perivascular infiltrate (white arrowhead)

Meta-transcriptomic discovery of novel WPDV and hepacivirus in Australia brushtail possums

In total, eight rRNA-depleted RNA sequencing libraries were constructed, generating from 23,480,309 to 26,509,764 reads per pool (Table 1). The first (Pool 1) contained the brain RNA from three suspected WPDV-infected possums. Our initial analysis did not identify any assembled viral contigs including any related to WPDV or other arteriviruses. As low abundance pathogens may not have sufficient coverage for assembly, the trimmed sequence reads were then further annotated directly (i.e. without assembly) against the NCBI NR database using Diamond BlastX. This analysis revealed six sequence reads (three sequence pairs) with homology to WPDV from New Zealand (GenBank accession: NC_026811). In addition to the WPDV reads, eight sequence reads from a novel hepacivirus were present in the same library. Further RNA sequencing libraries were then prepared from the three individual cases found to be positive for WPDV (of which all were co-infected with the hepacivirus) that included libraries from multiple tissues where available (Table 1: cases 3619, 7613, 2545). To determine if any other organisms were present in the remaining cases, we also sequenced a hepacivirus monoinfection (Table 1: case 2345) as well as a pool of RNAs from those cases negative for both WPDV and hepacivirus (Table 1: pool 2). Overall, the read abundance of WPDV and hepaciviruses in each library was low, ranging from 0.000007 to 0.000343% (Table 1).

Identification of divergent lineages of WPDV

Individual PCR assays were designed based on each of the three arterivirus-like sequence pairs identified from the initial RNA library and used to confirm the presence of WPD-like virus in specific cases. However, these PCR assays selectively amplified the viruses from the WPD-affected cases suggesting extensive genetic variation and the presence of different lineages of WPDVs in Australia. Accordingly, a more sensitive qRT-PCR assay was developed based on the conserved RdRp domain and specifically designed to target all WPDVs in this study including the NZ prototype strain (GenBank: JN116253.3). WPDV sequences were confirmed by qRT-PCR for the three WPD-affected possums including individual brain, liver and kidney tissues. Importantly, this pan-WPDV qRT-PCR assay did not identify any new cases and the tissues tested from the control possums (including case 1–4) were all negative. From this, we were also able to determine viral loads, which for WPD-positive tissues ranged from 21 to 4777 copies per pg RNA and varied according to different tissue types. The highest level of WPDV expression was observed in liver tissues from case 2525 (4777 copies/pg RNA) and from case 3619 (4125 copies/pg RNA), followed by kidney (570–1830 copies/pg RNA) and brain from each diseased possum (17–71 copies/pg RNA).

Comparative genomics and phylogenetics of WPDVs

PCR and sequencing of the WPDV replicase gene showed that the three individual possums carried distinctive (i.e. case-specific) nucleotide variation in WPDV, comprising two distinct lineages (Fig. 2). In addition, the WPDVs from different organs of each individual animal formed discrete clusters, strongly suggesting that these results were not due to cross-sample contamination. Strikingly, pairwise alignment of the partial replicase polyprotein 1ab nucleotide genome (291 bp) revealed that the two lineages of Australian WPDV strains exhibited almost the same level of sequence divergence between them as they did to the New Zealand strain (~ 70% nucleotide sequence similarity and 77% amino acid sequence similarity), indicative of substantial intra-specific virus diversity.

Fig. 2
figure 2

Identification of distinct lineages of WPDV in Australia. a Phylogenetic analysis of nucleotide sequences representing a partial (291 bp) replicase polyprotein 1ab region of WPDV, showing the NZ reference strains compared to the three Australian cases identified here. Two major lineages of the Australian viruses were identified - denoted Lineages 1 and 2 - and coloured in blue and red, respectively. Nodes show bootstrap values from 1000 replicates. The scale bar shows the number of nucleotide substitutions per site. b Pairwise comparisons of nucleotide similarity (%) across the viral genome

Long overlapping PCRs were used to determine the near complete genomes of two representative WPDV strains, here termed WPDV-AU1 (case 2545) and AU2 (case 3619). In comparison to the WPDV-NZ strain, the two WPDV-AU genomes retained the arterivirus-like features, encoding a replicase polyprotein ORF1ab, glycoprotein (GP) 2, 3, 4, 5, and membrane protein (MP) and nucleocapsid protein (NP) (Fig. 3a). The data generated here enabled us to determine the phylogenetic position of WPDV (both the Australian and New Zealand strain) within the Arteriviridae as a whole through a phylogenetic analysis of 1787 amino acids of the RdRp protein (NSP2). Notably, WPDV, derived from a marsupial, fell as a sister-group to those arteriviruses previously documented in placental mammals (Fig. 3b). In addition, the overall phylogeny of the Arteriviridae generally mirrors that of the vertebrate hosts from which they were sampled, with those viruses sampled from fish falling in the basal position, indicative of virus-host co-divergence for the entire evolutionary history of vertebrates.

Fig. 3
figure 3

Genome organisation and evolutionary relationships of WPDV to other arteriviruses. a Genome features and comparison of three lineages of WPDV - WPDV-AU1, WPDV-AU2 and WPDV-NZ (GenBank accession: NC_026811). b Phylogenetic analysis of the arterivirus RdRp protein. Nodes show bootstrap values obtained from 1000 replicates. Scale bar shows the number of amino acid substitutions per site. The tree was midpoint rooted for clarity only. c Pairwise nucleotide similarities (%) among the three WPDVs

A novel possum hepacivirus

Interestingly, our meta-transcriptomic analysis also identified a novel hepacivirus that we have termed possum hepacivirus. This virus retains the classical hepacivirus features of a single-stranded, positive-sense RNA genome of 9296 bp in length that encodes a single polyprotein. By further hepacivirus-specific PCR confirmation, we revealed the presence of the virus in liver, brain and kidney from three WPDV-affected possums and another two non-WPDV detected possums in our study, thereby revealing virus co-infection in animals with WPD (Table 2). A phylogenetic analysis of the complete polyprotein of hepacivirus (3260 amino acids) revealed that possum hepacivirus was most closely related to Norway rat hepacivirus 2 (accession number: YP_009109558.1) identified from brown rats (Rattus norvegicus) in New York City, US, with which it shared 41.2% sequence similarity in the polyprotein (Fig. 4).

Table 2 Summary of case details, diagnosis and PCR confirmation of WPDV and Brushtail possum hepacivirus infection status.
Fig. 4
figure 4

Evolutionary relationships of the novel hepacivirus identified from brushtailed possums. a Genome depiction of possum hepacivirus compared to its closest relative - Norway rat hepacivirus 2 (GenBank accession: NC_025673). b Phylogenetic analysis of the polyprotein gene (RdRp protein) of hepaciviruses. The novel possum hepacivirus is shown in red. Bootstrap values > 70% are shown for key nodes (1000 replicates). The tree was midpoint rooted for clarity only. The scale bar shows the number of amino acid substitutions per site

Discussion

Our retrospective investigation and meta-transcriptomic analysis with subsequent PCR confirmation revealed that brushtailed possums in mainland Australia that presented with a disease syndrome similar to that of WPD seen in New Zealand were infected with an arterivirus related to WPDV. In addition, three of the five animals presenting with WPD had evidence of co-infection with a novel hepacivirus, although whether this contributed to the symptoms and lesions observed is unclear and clearly requires additional investigation.

Previously, strong evidence for the association between WPD and a viral aetiology was obtained using primary cell culture system for virus isolation [14], WPD-specific PCR assays and in situ hybridisation [42]. By reproducing WPD in healthy possums with purified WPDV isolate inoculation, the histological changes highly met the natural infection status, fulfilling Koch’s postulates of causation. In addition, the establishment of the RT-qPCR and indirect ELISA assays for molecular and serological survey assists understanding disease prevalence in the wild population in New Zealand [42].

While the clinical and histological presentation of WPD as described here and by others varies geographically in wild possum populations in Australia, no active monitoring scheme is in place to explore the prevalence of the syndrome and diversity of associated viruses. Histologically, WPD in historical cases from Tasmania consistently had non-suppurative hepatitis and nephritis, similar to that described in the New Zealand form of WPD, and distinct from mainland Australian cases, which lack that feature. Importantly, we have designed a sensitive qRT-PCR assay that is able to detect divergent lineages of WPDV, including, WPDV-AU1, AU2, which should merit broader application in epidemiological studies of viral and syndrome prevalence.

Previous work [13] revealed that the viral loads of WPDV were significantly higher from the liver of experimental infected possums than other organs including kidney, brain, salivary glands and bladder, which is consistent with our qRT-PCR. Of note, these results also show consistently that the lowest viral loads were in brain tissues. This was also apparent in the WPDV affected case 7614, which contained no detectable WPDV reads in the brain tissue by RNA-seq, but was clearly WPDV positive by PCR and Sanger assays. Although the qRT-PCR assays in our study successfully detected various strains, the viral loads among different tissues were generally lower than reference results from experimentally infected possums in the New Zealand studies [42]. This finding may reflect the often chronic nature of infection in diseased Australian possums, which were frequently examined several weeks after the onset of disease, but is also likely explained by the archival nature of our samples and hence relatively low preservation might affect the RNA integrity due to degradation. However, since our testing number is comparatively small, speculation on the overall distribution of WPDVs in Australian possum populations remains challenging. The collection of additional samples from different geographical locations, closely related, and sympatric species, larger serological and molecular surveys, in conjunction with population and ecological information will clearly assist understanding virus variation and transmission dynamics in the wild.

Our phylogenetic analysis reveals that WPDV formed a distinct lineage and that there has likely been long-term co-divergence between arteriviruses and their vertebrate hosts over many millions of years. As it is also clear that our sampling of the Arteriviridae is sparse, it is inevitable that more will be discovered through expanded metagenomic studies of diverse vertebrates. More challenging will be determining exactly how WPDV became established in brushtailed possums in both Australia and New Zealand. Reports indicated that WPDV, or a virus causing similar neurological symptoms, has been present in free-living possum populations in New Zealand as far back as 1999 [10]. Our investigations track WPD cases back to 1983 and 1985 on mainland Australia and Tasmania respectively. As no evolutionary rate is available for WPDV it is currently impossible to perform a direct estimation of divergence times. However, the substantial sequence divergence between the Australian mainland and New Zealand strains, perhaps sufficient to warrant their classification as separate virus species, suggests that their separation is relatively old and may have occurred close to when possums were introduced from Tasmania into New Zealand in the 1830s. Additional testing of Tasmanian brushtail possums will be an integral to addressing the question of whether WPDV was translocated to New Zealand with possums, or emerged and evolved afterwards.

A particularly notable aspect of this study was the presence of a novel possum hepacivirus also identified in our WPD-affected and non WPDV-detected possums. Most animal hepaciviruses are associated with chronic infection and strong hepatotropism, leading to hepatitis, cirrhosis, and severe hepatopathy [22]. More recently, a divergent Wenling shark virus (WLSV) hepacivirus was discovered in the liver of graceful catsharks (Proscyllium habereri) using the same bulk meta-transcriptomic approach as utilized here [43]. As the first documented hepacivirus identified from a marsupial, the novel Brushtail possum hepacivirus revealed here expands the host range of this important group of viruses and again highlights the potential missing diversity of genus Hepacivirus. In addition, that the Brushtail possum hepacivirus falls in a clade of rodent hepaciviruses indicates that cross-species transmission has played a key role in shaping phylogenetic patterns. Clearly, determining the natural infection routes and pathogenesis of this virus in possums merits additional work. For example, the establishment of in situ hybridisation based techniques will help reveal the nature of host-pathogen interactions. Whether the novel possum hepacivirus identified here contributes to clinical disease or reduced fitness, alone or in conjunction with WPDV, and the status of both agents in wild populations remains largely unexplored, but is clearly a key area for future study including the presence in New Zealand animals or any of the inoculums used from previous challenge studies. More generally, an enhanced understanding of the roles of these and closely related organisms will shed important light on virome of non-eutherian mammals.

Conclusions

The meta-transcriptomic approach described has elucidated the pathogen likely responsible for a disease syndrome first detected in Australian mainland wildlife 36 years ago. Factors that may have contributed to this diagnostic delay and difficulty include the presentation in a common species, the sporadic rather than outbreak-based emergence of the syndrome, and the potential pre-conception that the disease posed no threat to human or livestock health. Importantly, meta-transcriptomics is becoming more rapid and cost-effective, and has demonstrated its capacity to remove barriers that delay or obstruct traditional wildlife disease investigation. The data and knowledge generated from meta-transcriptomics studies such as this will inform risk assessments addressing the potential threats and impacts of emergent pathogens within a One Health paradigm.

Availability of data and materials

The RNA sequencing data generated in this study have been deposited in the GenBank Sequence Read Archive under accession numbers PRJNA587034. All consensus genome sequences of identified viruses have been deposited in GenBank under accession numbers MN635447 to MN635449.

Abbreviations

NCBI:

National Center for Biotechnology Information

NZ:

New Zealand

RT-PCR:

Reverse transcription polymerase chain reaction

WPD:

Wobbly possum disease

WPDV:

Wobbly possum disease virus

References

  1. Tompkins DM, Carver S, Jones ME, Krkosek M, Skerratt LF. Emerging infectious diseases of wildlife: a critical perspective. Trends Parasitol. 2015;31:149–59.

    Article  PubMed  Google Scholar 

  2. Grogan LF, Berger L, Rose K, Grillo V, Cashins SD, Skerratt LF. Surveillance for emerging biodiversity diseases of wildlife. PLoS Pathog. 2014;10:e1004015.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. Gurevitch J, Padilla DK. Are invasive species a major cause of extinctions? Trends Ecol Evol. 2004;19:470–4.

    Article  PubMed  Google Scholar 

  4. Di Giallonardo F, Holmes EC. Exploring host-pathogen interactions through biological control. PLoS Pathog. 2015;11:e1004865.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. King CM. The handbook of New Zealand mammals. The mammal society, New Zealand branch. Wellington: Oxford University Press; 1990.

    Google Scholar 

  6. Coleman JD, Coleman MC, Warburton B. Trends in the incidence of tuberculosis in possums and livestock, associated with differing control intensities applied to possum populations. New Zealand Vet J. 2006;54:52–60.

    Article  CAS  Google Scholar 

  7. Sanhueza JM, Heuer C, Wilson PR, Benschop J, Collins-Emerson JM. Seroprevalence and risk factors for leptospira seropositivity in beef cattle, sheep and deer farmers in New Zealand. Zoonoses Public Health. 2017;64:370–80.

    Article  CAS  PubMed  Google Scholar 

  8. Eymann J, Smythe LD, Symonds ML, Dohnt MF, Barnett LJ, Cooper DW, et al. Leptospirosis serology in the common brushtail possum (Trichosurus vulpecula) from urban Sydney. Australia J Wildl Dis. 2007;43:492–7.

    Article  PubMed  Google Scholar 

  9. Dunowska M, Biggs PJ, Zheng T, Perrott MR. Identification of a novel nidovirus associated with a neurological disease of the Australian brushtail possum (Trichosurus vulpecula). Vet Microbiol. 2012;156:418–24.

    Article  CAS  PubMed  Google Scholar 

  10. Mackintosh CG, Crawford JL, Thompson EG, McLeod BJ, Gill JM, O'Keefe JS. A newly discovered disease of the brushtail possum: wobbly possum syndrome. New Zealand Vet J. 1995;43:126.

    Article  CAS  Google Scholar 

  11. Perrott MR, Meers J, Cooke MM, Wilks CR. A neurological syndrome in a free-living population of possums (Trichosurus vulpecula). New Zealand Vet J. 2000;48:9–15.

    Article  CAS  Google Scholar 

  12. Giles J, Perrott M, Roe W, Dunowska M. The aetiology of wobbly possum disease: reproduction of the disease with purified nidovirus. Virology. 2016;491:20–6.

    Article  CAS  PubMed  Google Scholar 

  13. Perrott MR, Wilks CR, Meers J. Routes of transmission of wobbly possum disease. New Zealand Vet J. 2000;48:3–8.

    Article  CAS  Google Scholar 

  14. Giles J, Perrott M, Roe W, Shrestha K, Aberdein D, Morel P, et al. Viral RNA load and histological changes in tissues following experimental infection with an arterivirus of possums (wobbly possum disease virus). Virology. 2018;522:73–80.

    Article  CAS  PubMed  Google Scholar 

  15. Giles JC, Perrott MR, Dunowska M. Primary possum macrophage cultures support the growth of a nidovirus associated with wobbly possum disease. J Virol Meth. 2015;222:66–71.

    Article  CAS  Google Scholar 

  16. Hartley WJ. Central nervous system disorders in the brush-tail possum in eastern Australia. New Zealand Vet J. 1993;41:44–5.

    Article  Google Scholar 

  17. Ma G, Dennis M, Rose K, Spratt D, Spielman D. Tawny frogmouths and brushtail possums as sentinels for Angiostrongylus cantonensis, the rat lungworm. Vet Parasitol. 2013;192:158–65.

    Article  PubMed  Google Scholar 

  18. Giles JC, Johnson W, Jones G, Heuer C, Dunowska M. Development of an indirect ELISA for detection of antibody to wobbly possum disease virus in archival sera of Australian brushtail possums (Trichosurus vulpecula) in New Zealand. New Zealand Vet J. 2018;66:186–93.

    Article  CAS  Google Scholar 

  19. Kuhn JH, Lauck M, Bailey AL, Shchetinin AM, Vishnevskaya TV, Bao Y, et al. Reorganization and expansion of the nidoviral family Arteriviridae. Arch Virol. 2016;161:755–68.

    Article  CAS  PubMed  Google Scholar 

  20. O'Dea MA, Jackson B, Jackson C, Xavier P, Warren K. Discovery and partial genomic characterisation of a novel nidovirus associated with respiratory disease in wild shingleback lizards (Tiliqua rugosa). PLoS One. 2016;11:e0165209.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Zhang J, Finlaison DS, Frost MJ, Gestier S, Gu X, Hall J, et al. Identification of a novel nidovirus as a potential cause of large scale mortalities in the endangered Bellinger River snapping turtle (Myuchelys georgesi). PLoS One. 2018;13:e0205209.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Gulyaeva A, Dunowska M, Hoogendoorn E, Giles J, Samborskiy D, Gorbalenya AE. Domain organization and evolution of the highly divergent 5′ coding region of genomes of arteriviruses, including the novel possum nidovirus. J Virol. 2017;91:e02096–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Hartlage AS, Cullen JM, Kapoor A. The strange, expanding world of animal hepaciviruses. Annu Rev Virol. 2016;3:353–75.

    Article  CAS  Google Scholar 

  24. Kapoor A, Simmonds P, Gerold G, Qaisar N, Jain K, Henriquez JA, et al. Characterization of a canine homolog of hepatitis C virus. Proc Natl Acad Sci U S A. 2011;108:11608–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Burbelo PD, Dubovi EJ, Simmonds P, Medina JL, Henriquez JA, Mishra N, et al. Serology-enabled discovery of genetically diverse hepaciviruses in a new host. J Virol. 2012;86:6171–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Firth C, Bhat M, Firth MA, Williams SH, Frye MJ, Simmonds P, et al. Detection of zoonotic pathogens and characterization of novel viruses carried by commensal Rattus norvegicus in New York City. mBio. 2014;5:e01933–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Quan PL, Firth C, Conte JM, Williams SH, Zambrana-Torrelio CM, Anthony SJ, et al. Bats are a major natural reservoir for hepaciviruses and pegiviruses. Proc Natl Acad Sci U S A. 2013;110:8194–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Chu L, Jin M, Feng C, Wang X, Zhang D. A highly divergent hepacivirus-like flavivirus in domestic ducks. J Gen Virol. 2019;100:1234–40.

    Article  PubMed  Google Scholar 

  29. Zheng T. Characterisation of two enteroviruses isolated from Australian brushtail possums (Trichosurus vulpecula) in New Zealand. Arch Virol. 2007;152:191–8.

    Article  CAS  PubMed  Google Scholar 

  30. Thomson D, Meers J, Harrach B. Molecular confirmation of an adenovirus in brushtail possums (Trichosurus vulpecula). Virus Res. 2002;83:189–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Perrott MR, Meers J, Greening GE, Farmer SE, Lugton IW, Wilks CR. A new papillomavirus of possums (Trichosurus vulpecula) associated with typical wart-like papillomas. Arch Virol. 2000;145:1247–55.

    Article  CAS  PubMed  Google Scholar 

  32. Eden JS, Rose K, Ng J, Shi M, Wang Q, Sintchenko V, et al. Francisella tularensis ssp. holarctica in ringtail possums, Australia. Emerg Infect Dis. 2017;23:1198–201.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Shi M, Lin XD, Chen X, Tian JH, Chen LJ, Li K, et al. The evolutionary history of vertebrate RNA viruses. Nature. 2018;556:197–202.

    Article  CAS  PubMed  Google Scholar 

  34. Shi M, Lin XD, Tian JH, Chen LJ, Chen X, Li CX, et al. Redefining the invertebrate RNA virosphere. Nature. 2016;540:539–43.

    Article  CAS  PubMed  Google Scholar 

  35. Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat Protoc. 2013;8:1494–512.

    Article  CAS  PubMed  Google Scholar 

  36. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Marchler-Bauer A, Derbyshire MK, Gonzales NR, Lu S, Chitsaz F, Geer LY, et al. CDD: NCBI's conserved domain database. Nucleic Acids Res. 2015;43:D222–6.

    Article  CAS  PubMed  Google Scholar 

  38. BBMap short-read aligner, and other bioinformatics tools. 2016. Available at: http://sourceforge.net/projects/bbmap/.

  39. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25:1972–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Nguyen LT, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32:268–74.

    Article  CAS  PubMed  Google Scholar 

  42. Dunowska M, Gopakumar G, Perrott MR. Development of a real-time reverse transcription PCR assay for detection of a novel nidovirus associated with a neurological disease of the Australian brushtail possum (Trichosurus vulpecula). New Zealand Vet J. 2013;61:286–91.

    Article  CAS  Google Scholar 

  43. Shi M, Lin XD, Vasilakis N, Tian JH, Li CX, Chen LJ, et al. Divergent viruses discovered in arthropods and vertebrates revise the evolutionary history of the Flaviviridae and related viruses. J Virol. 2016;90:659–69.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This research was supported by Taronga Conservation Society Australia, Taronga Conservation Science Initiative, and New South Wales National Parks and Wildlife Service. Our thanks are extended to Jane Hall for diagnostic support and data management, Dr. Hannah Bender for graphic design of the photomicrographic plate, and to Drs. Cathy Shilton, Cheryl Sangster, Lydia Tong, Richard Montali, Rod Reece, Shannon Donahoe, Dr. Andrew Davis Ferguson and Dr. David Obendorf for their diagnostic contributions on individual cases. We also acknowledge The University of Sydney HPC service (Artemis) for providing high-performance computing resources that have contributed to the research reported in this paper.

Funding

ECH is supported by an ARC Australian Laureate Fellowship (FL170100022).

Author information

Authors and Affiliations

Author notes

  1. William J. Hartley is deceased. This paper is dedicated to his memory.

    • William J. Hartley
Authors

Contributions

Conceived the project - JSE, KR, ECH; collected samples - WJH, KR; performed laboratory work - WSC, JSE, KR; analysed the data - WSC, JSE, MS, KR; wrote the paper - WSC, JSE, MS, KR, ECH. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Edward C. Holmes.

Ethics declarations

Ethics approval and consent to participate

Samples were collected under the Opportunistic Sample Collection Program of the Taronga Animal Ethics Committee, and under scientific licences #SL10469 and SL100104 issued by the NSW Office of Environment and Heritage.

Consent for publication

With the exception of WJH (deceased), all authors have agreed to submission of the final version of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

PCR primers and qPCR primers for WPDV used in this study.

Additional file 2: Table S2.

Presentation and pathology of wobbly possum disease in brushtail possums in Australia. This table chronologically summarises the clinical signs and histological findings in brushtail possums considered to fit the syndrome description for wobbly possum disease. Most possums originate from New South Wales on mainland Australia, except for those denoted with a T that originate from Tasmania.

Additional file 3: Video S1.

Video showing the clinical manifestations of wobbly possum disease syndrome in an Australian brushtail possum.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chang, WS., Eden, JS., Hartley, W.J. et al. Metagenomic discovery and co-infection of diverse wobbly possum disease viruses and a novel hepacivirus in Australian brushtail possums. One Health Outlook 1, 5 (2019). https://doi.org/10.1186/s42522-019-0006-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s42522-019-0006-x

Keywords