Next Article in Journal
Binding Mechanism Elucidation of the Acute Respiratory Disease Causing Agent Adenovirus of Serotype 7 to Desmoglein-2
Previous Article in Journal
Regulation of Polyomavirus Transcription by Viral and Cellular Factors
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 10
10.3390/v12101073
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Wei-Shan Chang
1,
Ci-Xiu Li
1,
Jane Hall
2,
John-Sebastian Eden
1,3,
Timothy H. Hyndman
4,
Edward C. Holmes
1,* and
Karrie Rose
2,*
1
Marie Bashir Institute for Infectious Diseases and Biosecurity, School of Life and Environmental Sciences and School of Medical Sciences, The University of Sydney, Sydney, NSW 2006, Australia
2
Australian Registry of Wildlife Health, Taronga Conservation Society Australia, Mosman, NSW 2088, Australia
3
Westmead Institute for Medical Research, Centre for Virus Research, Westmead, NSW 2145, Australia
4
School of Veterinary Medicine, Murdoch University, Murdoch, WA 6150, Australia
*
Authors to whom correspondence should be addressed.
Viruses 2020, 12(10), 1073; https://0-doi-org.brum.beds.ac.uk/10.3390/v12101073
Submission received: 26 August 2020 / Revised: 19 September 2020 / Accepted: 22 September 2020 / Published: 25 September 2020
(This article belongs to the Section Animal Viruses)
Browse Figures
Versions Notes

Abstract

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

1. Introduction

Squamates (scaled reptiles) are one of the most diverse groups of vertebrate fauna in Australia [1]. Bearded dragons (Pogona spp.), including central bearded dragons (Pogona vitticeps), are native reptiles that inhabit wide geographic areas throughout the arid and semiarid regions of Australia. Bearded dragons are increasingly popular exotic pets, and are now one of the most common companion lizard species in Australia, Europe, Asia, and North America [2]. Bearded dragons are also ideal for laboratory studies: their well-established history in captivity makes them suitable for investigations of reptile physiology, and they have unusual patterns of chromosomal, genomic, and temperature-based sex-determination, including sex reversal with specific behavioral consequences [3,4]. Despite their significance, our knowledge of reptile pathogens, particularly viruses, remains limited. Indeed, to date, the only viruses known to infect bearded dragons are double-strand DNA viruses from the families Iridoviridae [5,6,7,8] and Adenoviridae [9,10], single-stranded DNA viruses from the Parvoviridae [11,12,13] and RNA viruses from the Paramyxoviridae [14,15].
The diagnosis of microbial infections in reptiles is complex and often associated with multiple synergistic and predisposing factors. It is likely that undetected viruses contribute to a number of infectious disease outbreaks in captive reptiles, including bearded dragons. To this end, we investigated two unusual mortality events of bearded dragons in a research colony in Australia in 2014 and 2020. In 2014, approximately 40 central bearded dragons died after emergence from brumation in a colony of over 400 animals located in Canberra, Australia. Affected animals were either found dead or were found listless and died within 24 h. In 2020, the same colony experienced low-grade post-brumation mortality followed several months later by star gazing and poor body condition in juveniles and acute death in an adult bearded dragon.
Brumation is an extreme hypometabolic state used by some reptiles to cope with low or unpredictable food availability and unfavorable seasonal conditions for the duration of winter. A number of specific protective measures have been identified during this hibernation in bearded dragons, including increased neuroprotection in the brain, maintenance of heart function through hypertrophy, and upgrading antioxidant capacity and mitochondrial maintenance by skeletal muscle atrophy [16]. Emergence from brumation may also represent a period of increased disease susceptibility in reptiles, as shown by downregulated transcription of several genes responsible for microbial pathogen defense, cellular and oxidative stress, and cell differentiation and growth [17].
Herein, we employed PCR testing and meta-transcriptomic approaches combined with gross and microscopic pathology to investigate the possible involvement of viral pathogens associated with unknown mortality events in a captive colony of bearded dragons.

2. Materials and Methods

2.1. Sample Collection and Processing

Tissues were collected from two disease outbreaks in a research colony in Canberra, Australia. Five dead bearded dragons were submitted for postmortem examination during a mortality event in 2014. In 2020, 12 live, sick animals (including 8 central bearded dragons (Pogona vitticeps), 3 jacky lizards (Amphibolurus muricatus), and 1 grassland earless dragon (Tympanocryptis pinguicolla)) were examined, sedated with 20 mg/kg of im alfaxalone (Alfaxan® CD-RTU 10 mg/mL, Jurox, Australia), delivered into the biceps to rapidly induce a surgical plane of anesthesia, followed by 1 ml/kg of iv pentobarbitone (Lethobarb, 325 mg/mL, Virbac, Australia) to effect euthanasia. Gross postmortem examinations were conducted immediately, and fresh portions of cervical spinal cord, liver, heart, spleen, 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, sectioned, stained with hematoxylin and eosin, and mounted with a cover slip prior to examination by light microscopy. Giemsa, Ziehl–Neelsen, and periodic acid-Schiff stains were applied to a subset of embedded tissue samples to identify bacteria, and to exclude the presence of mycobacteria and fungi. Aerobic, anaerobic, and fungal cultures were undertaken on lung and liver samples aseptically collected from the five dragons that died in 2014.
Proliferation of respiratory epithelium was graded on a scale of 0–4, where 0 indicated no discernable lesion and 4 represented severe and extensive proliferation of infundibular and mesobronchial epithelium. Samples were collected under the Opportunistic Sample Collection Program of the Taronga Conservation Society Australia’s Animal Ethics Committee (Approval #3b/12/18, 20 December 2018) and under scientific license #SL100104 issued by the NSW Office of Environment and Heritage and valid 2 May 2011 to 30 April 2021. Samples were submitted to Murdoch University for PCR testing of pneumotropic viruses, and to the University of Sydney for meta-transcriptomic analysis and further PCR testing.
For PCR testing for pneumotropic viruses, total nucleic acid was extracted from tissues using the MELTTM Total Nucleic Acid Isolation System (Ambion, TX, USA) according to the manufacturer’s instructions. Total nucleic acid was eluted into 30 µL elution buffer. For metatranscriptomics and PCR testing for parvoviruses and circoviruses, the RNA extracted from the spinal cord, liver, lung, and kidney of diseased animals was processed using the RNeasy Plus Mini Kit (Qiagen, Hilden, Germany). Initially, frozen tissue was partially thawed and submerged in RLT lysis buffer containing 1% β-mercaptoethanol and 0.5% reagent DX before being homogenized using a hand-operated TissueRupture (Qiagen). To pool samples in equal proportions, RNA concentrations and integrity were validated using a NanoDrop spectrophotometer (ThermoFisher Scientific, MA, USA) and a TapeStation (Agilent, CA, USA). Illumina TruSeq stranded RNA libraries were prepared on the pooled samples following rRNA depletion using a RiboZero Gold rRNA removal kit (Epidemiology). Finally, 150 bp paired-end sequencing of the rRNA-depleted RNA library was generated on an Illumina NovaSeq 6000 service system at Australian Genome Research Facility (AGRF), Melbourne.

2.2. PCR Testing for Pneumotropic Viruses

Nucleic acid from the kidney, liver, and lung was pooled and tested for sunshineviruses, orthoreoviruses, and ferlaviruses using one-step reverse transcription (RT)-PCR. For each virus genus, 1 µL extracted nucleic acid was added to 0.8 µL SuperScript® III RT/Platinum® Taq Mix (Invitrogen, Victoria, Australia), 10 µL 2× Reaction Mix, and 1 µM (final concentration) of each primer, and was then made up to a final volume of 20 µL. For nested PCRs, 1 µL PCR product was used as template for the second round of amplification using Platinum® PCR Supermix (Invitrogen) in a final volume of 20 µL. For sunshineviruses, the primer pair SunshineS2–SunshineAS2 was used [18]. For orthoreoviruses, the primer pairs 1607F–2608R and 2090F–2334R were used for the first and second rounds, respectively [19]. For ferlaviruses, the primer pairs L5–L6 and L7–L8 were used, respectively [20]. Sunshine Coast virus [21], a reptile orthoreovirus (kindly provided by Dr Rachel Marschang) and a ferlavirus (American Type Culture Collection VR-1408) were used as positive controls. PCR products were visualised with agarose gel electrophoresis.

2.3. Pathogen Discovery Using Meta-transcriptomics

Our meta-transcriptomic approach to pathogen discovery was based on those previously employed by our group [22,23]. RNA sequencing reads were trimmed for quality using Trimmomatic [24] before de novo assembly with Trinity, version 2.5.1 [25]. Assembled sequence contigs were annotated using both nucleotide BLAST searches against the NCBI nt database and Diamond Blastx against the NCBI nr database [26], with e-value cut-offs of <10−6 and <10−5, respectively. Open reading frames were predicted from the potential viral contigs using Geneious v11.1.5 [27], with gene annotation and functional predictions made against the Conserved Domain Database (CDD) [28]. To evaluate the virus abundance and coverage, reads were mapped back to the genome with BBmap [29].

2.4. PCR Testing for Chaphamaparvoviruses and Circoviruses

Virus-specific PCR primers were designed based on identified transcripts from the RNA-seq data (Table 1). Accordingly, the extracted RNA (5 µL) from each tissue of individual cases was reverse transcribed using the SuperScript IV VILO cDNA synthesis system (Invitrogen). The cDNA generated from each sampled tissue (2 µL) was used for viral specific PCRs targeting regions identified by RNA-Seq. For BDCV, the primer pairs BDCV_F2 and BDCV_R1 were used (cycling conditions = 98 °C × 1 m, 35 × (98 °C × 10 s, 62 °C × 10 s, 72 °C × 45 s)). For BDchPV, the primers BDChPV_F4 and BDChPV_R3 were used (cycling conditions 98 °C × 1 m, 35 × (98 °C × 10s, 62.2 °C × 10 s, 72 °C × 60 s)). All PCRs were performed using Platinum SuperFi DNA polymerase (Invitrogen) with a final concentration of 0.2 μM for both forward and reverse primers. PCR products were visualized with agarose gel electrophoresis and Sanger sequenced at the AGRF.

2.5. Phylogenetic Analysis

Phylogenetic trees were estimated based on MAFFT [30] alignments of the conserved nonstructural region of the viruses identified here, with all analyses utilizing representative members of each virus family taken from NCBI/GenBank. Maximum likelihood (ML) phylogenies were inferred using IQ-TREE (version 1.4.3) [31], employing the best-fit LG+F+Γ4 model of amino acid substitution on the nonstructural (NS) protein (717 amino acids) and the LG+F+I+Γ4 model on the replicase (rep) protein (274 amino acids). Statistical support for individual nodes was estimated via bootstrap resampling (1000 replicates). Data were visualized using Figtree 1.4.2 (http://tree.bio.ed.ac.uk/software/figtree/). The genome sequence of BDCV and BDchPV are available on NCBI/GenBank (accession numbers: MT732118 and MT73219). Raw sequencing reads are available at the Sequence Read Archive (SRA) under accession PRJNA644669.

2.6. Mining the Sequence Read Archive (SRA)

To help determine if BDCV and BDchPV were present in other bearded dragon or other reptile species, we screened publicly accessible high-throughput sequencing data available on the NCBI SRA database (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra). Accordingly, a large RNA-seq data set was obtained using the NCBI SRA toolkit (version 2.9.2) from the bearded dragon genus Pogona (NCBI taxid:52201). Retrieved FASTQ reads were then subjected to a blast analysis using Diamond v.0.9.25 [32] against customized databases containing the core genes from reference circoviruses and parvoviruses, employing an e-value of cut-off of 1 × e−5.

3. Results

3.1. Clinical and Histopathological Findings of Diseased Captive Reptiles

The mass mortality and morbidity of up to 40 adult central bearded dragons in 2014, and intermittent ill health and star gazing in hatchlings from the same research colony in 2020, were both accompanied by unusual and severe proliferation of pulmonary epithelium. In total, 17 cases were investigated in this study. The severe proliferation of pulmonary epithelium was found in 7 cases, including 5 bearded dragons from the 2014 outbreak, and 1 bearded dragon and 1 Jacky lizard from 2020 mortality events. The infundibular mucosa, which is normally squamous, exhibited metaplasia ranging from columnar to pseudostratified columnar and ciliate. The mesobronchial epithelium was similarly thickened, pseudostratified, and dysplastic with multifocal loss of cilia. Amorphous basophilic inclusions were evident within enlarged nuclei of infundibular epithelial cells in two out of seven affected animals. Eosinophilic intracytoplasmic inclusion bodies were evident within the respiratory epithelium of three out of seven affected animals. Rarely binucleate epithelial cells and syncytia were identified within mesobronchial epithelium.
Similar inclusions were evident within hepatocytes of three out of five dragons and renal tubular epithelium of two out of five dragons from the initial 2014 mortality event. Pasteurellaceae bacteria of unknown species were isolated from the lung and liver of each of the five dragons examined in 2014, and Aeromonas hydrophila was isolated in the lung of two of these animals; however, these were not considered contributory to pulmonary proliferation nor inclusion body formation. A summary of the signalment, gross, and histological findings in each animal in the initial 2014 outbreak is summarized in Table S1, and the histological lesions are shown in Figure 1.

3.2. PCR Testing for Pneumotropic Viruses

All samples were PCR-negative for sunshineviruses, ferlaviruses, and orthoreoviruses.

3.3. Metatranscriptomic Pathogen Discovery

Sequencing of a pooled rRNA-depleted RNA-seq library resulted in 81,188,754 raw reads. After filtering, 54,804,486 paired-end trimmed reads were generated, which were then de novo assembled into 266,118 contigs. Analyses of these read data revealed the presence of two DNA viruses, denoted here as bearded dragon circovirus (BDCV) and bearded dragon chaphamaparvovirus (BDchPV), with read abundances of 0.0006% and 0.003%, respectively. Virus-targeted PCR assays (Table 1) and Sanger sequencing were used to recover the full genome sequences of these two viruses. Virus-specific PCR was also used on spinal cord, lung, liver, and kidney to provide insight into the distribution of these viruses. No other reptilian associated circoviruses and parvoviruses were detected using SRA screening.

3.4. Genome Characterization of a Novel Bearded Dragon Circovirus

We identified an abundant circoviruslike contig (1693 bp, 304 transcripts per million (TPM)) from our pooled lung, kidney, and liver samples. The typical features of the genus Circovirus are nonenveloped, icosahedral, single-stranded circular DNA (ssDNA) genomes of approximately 1.5–2 kb in size. The full genome of BDCV was recovered through PCR assays and Sanger sequencing. The BDCV genome comprises 1761 bp of circular DNA with a 43.8% GC content, encoding rep and cap genes of 308 and 218 amino acid residues, respectively.
The replicase protein of circoviruses introduces an endonucleolytic nick within the typical stem-loop structure with nonamer sequence 5′-AGTATTAC-3′ in the intergenic region of the genome, thereby initiating the rolling circle replication (RCR). Comparison with other members of the Circoviridae revealed the presence of the classical conserved motifs, including RCR motifs, as well as SF3 helicase motifs (Table 2 and Figure 2).
To determine the evolutionary history of BDCV with respect to other circoviruses, we performed a phylogenetic analysis of the Rep protein. Notably, BDCV clustered robustly (99% bootstrap support; 77.85% amino acid pairwise identity in Rep genes) with a bat-associated virus lineage that includes Tadarida brasiliensis circovirus 1 (accession number YP_009170674.1) and Rhinolophus ferrumequinum circovirus 1 (YP_009506275.1), recently identified in metagenomic analyses of bat gut samples from Brazil and China, respectively. At the nucleotide level, the BDCV genome showed 69.33% and 71.53% sequence similarity with Tadarida brasiliensis circovirus 1 and Rhinolophus ferrumequinum circovirus 1, respectively. Since the species demarcation threshold for circoviruses is <80% genomewide pairwise identity, we suggest that BDCV likely represents a new species within the Circoviridae.

3.5. Identification of a Highly Divergent Bearded Dragon Chaphamaparvovirus

We identified two abundant chaphamaparvoviruslike transcripts in the RNA-seq library. The Parvoviridae are a family of small, nonenveloped, ssDNA animal viruses with a linear genome of 4–6 kb in length. Chaphamaparvoviruses are a recently defined genus in which a variety of component virus species have now been identified [33,34,35,36]. Four sets of bridging PCR assays were designed to recover a near complete virus genome. This revealed that the near complete virus genome comprised 4181 nt with two distinct ORFs that encode the nonstructural protein (NS, 633 aa) and the structural protein (VP, 600 aa). We further utilized specific primers to amplify the targeted NS (BDchPV_F2 and R1) and VP region (BDchPV_F4 and R3) for virus screening in different organs (Table 1).
Two conserved motifs were identified in the NS protein, including the putative endonuclease metal coordination motif “HLH” and the helicase motif “GPASTGKS” at amino acid positions 110–112 and 310–317, respectively (Figure 3). Additionally, the potential accessory ORF1 (positions 74–358, 95 aa) and ORF3 (p15-like) (positions 136–543, 136 aa), previously identified in most other amniote-associated chaphamaparvoviruses, were both detected in BDchPV. The conserved PLA2 and G-rich motifs widely present in other members of Parvoviridae were absent from the VP protein, as is the case for other chaphamaparvoviruses [37,38].
Phylogenetic analysis based on the predicted amino acid sequence of the complete NS protein revealed that BDchPV fell within the Chaphamaparvovirus lineage, clustering with viruses associated with fish and crocodiles. Tilapia parvovirus is a recently identified Chaphamaparvovirus isolated from the feces of farmed tilapia and crocodiles in China [36]. However, the amino acid sequences of the novel BDchPV identified here shared only 32.7% and 31.6% pairwise identity in NS gene and VP gene with this virus (accession no: MN162688.1), respectively (Figure 4).

3.6. Prevalence of BDCV and BDchPV through PCR Screening

The nucleic acid extracted from archived bearded dragon cases in the 2014 (n = 5) and captive reptile cases in 2020 (n = 12) outbreaks were screened using BDCV and BDchPV specific primers, respectively (Table 3). The circovirus was detected from kidney and lung samples from one bearded dragon (case: 10043.3) in the original respiratory disease outbreak in 2014, and of three out of twelve liver samples from the second outbreak (2020) were BDCV positive. In total, BDchPV was detected in for liver samples, one brain sample, and one lung sample, including two original respiratory cases in 2014 and three additional cases in 2020 (Table 3). Collectively, the case prevalence of both viruses from the two outbreaks was 29.4% (five out of seventeen) for BDCV and 29.4% (five out of seventeen) for BDchPV. Sanger sequencing results from PCR products showed 95%–99% nucleotide identity to the index cases.

4. Discussion

Our knowledge of the viruses of herpetofauna is expanding rapidly. Herein, we describe the meta-transcriptomic discovery and characterization of two complete genomes of DNA viruses, as well as their prevalence from mass mortality events of captive reptiles.
Circoviruses are small DNA viruses that infect a wide range of vertebrates, with most pathogenic consequences reported in association with immunocompromised or immunosuppressed hosts, including canine [39], psittacine [40], porcine [41], and mink [42] species. In 2019, a circovirus was detected in the liver and gut of black-headed pythons (Aspidites melanocephalus) with spinal osteopathy [43]. The virus was then identified in tissues of a Boelen’s python (Morelia boeleni) and two annulated tree boas (Corallus annulatus). That investigation claimed to be the first to describe a reptilian circovirus, although other circo-like viruses had previously been described in testudines and squamates [44,45]. In all cases, the association between the circovirus and disease was either weak or absent. Endogenous circoviruses have also been detected in snake genomes [46]. It is notable that the bearded dragon circovirus newly identified here clustered with a bat-associated circoviruses lineage. However, these bat viruses were identified through metagenomic sequencing of environmental related samples (i.e., fecal viromes) such that their true hosts remain unclear.
Parvoviruses are associated with disease in a variety of host species, ranging from canines [47] to livestock [48,49,50,51], rodents [36], and humans [52]. Reptile parvoviruses were identified in several snake species (snake adeno-associated virus) [11,53,54,55,56,57] and some lizard species including bearded dragons [11,12]. To date, most of the exogenous reptilian parvoviruses identified belong to the genus Dependoparvovirus, and it has been thought that these viruses require helper viruses (usually adenoviruses) for replication [13], although more recent work is challenging this notion [58]. Chaphamaparvoviruses (ChPVs) are a newly identified genus of parvoviruses [59] identified in rodents, birds, pigs, bats, Tasmanian devils, dogs, cats, primates, and even invertebrates [33,34,35,36,37,60]. Although a ChPV was previously identified in crocodile faeces [36], it is possible that this originated in the tilapia fish fed to these crocodiles. Our documentation of a novel and exogenous reptile chaphamaparvovirus in central bearded dragons extends the host range of these parvoviruses.
Interestingly, coinfection with circoviruses and parvoviruses is increasingly reported in different animal species, often in association with immunosuppression and exacerbated disease severity. For example, previous studies have indicated that Porcine circovirus 2 (PCV-2) associated diseases are augmented by concurrent viral infections such as Porcine parvovirus, and these viruses may serve as important cofactors in the pathogenesis of Porcine multisystemic wasting syndrome (PMWS) [61]. Similarly, Canine circovirus 1 (CaCV-1) and Canine parvovirus 2 (CPV-2) infection are linked to recurrent outbreaks of bloody diarrhea and sudden death in puppies [62], and a novel lorikeet chaphamaparvovirus coinfected with psittacine circovirus (beak and feather disease virus) in wild birds has been described [63]. Recently, a similar coexistence of Tasmanian devil-associated circovirus and a chaphamaparvovirus was identified in a Tasmania devil metagenomic virome study, but with no disease association [64].
Despite our description of these pathogens, their roles in the disease syndromes described here remain uncertain and require additional investigation. Indeed, the examination of feeding and husbandry practices surrounding brumation is warranted, particularly since the animals that did not undergo brumation appeared to be unaffected during the mortality event. Although the association between the two viruses and proliferative lung disease was not clearly established, the role of these viruses as predisposing factors for disease or reduced fitness merits further investigation.
Reptiles comprise a considerable proportion of the vertebrate biomass of many ecosystems and play an important role in the food-web, and the delivery of ecosystem services through pollination, seed dispersal, and pest control. Australia has the highest diversity of lizards of any country, yet our understanding of the biological threats to these animals is scarce. It has been difficult to obtain an accurate understanding of the health status of wildlife that are free-living, maintained for research, education, or species recovery, using traditional diagnostic techniques. This study illustrates the value of metagenomic approaches to disease investigation as a complement to traditional histology and pathogen culture. The identification of bearded dragon circovirus and the bearded dragon chaphamaparvovirus provides insight into reptilian viral diversity and reptile health, and undoubtedly merits broader surveillance.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/12/10/1073/s1. Table S1: Presentation, pathology in central bearded dragon morbidity and mortality event, 2014.

Author Contributions

Conceptualization, E.C.H., K.R. and W.-S.C.; methodology, W.-S.C., C.-X.L., T.H.H., J.H. and K.R.; software, W.-S.C.; laboratory analysis, W.-S.C., C.-X.L., J.-S.E. and K.R.; formal analysis, W.-S.C.; investigation, W.-S.C.; resources, K.R. and E.C.H.; data curation, W.-S.C.; writing—original draft preparation, W.-S.C. and K.R.; writing—review and editing, W.-S.C., C.-X.L, J.H., J.-S.E., T.H.H., E.C.H. and K.R.; visualization, W.-S.C. and K.R.; supervision, E.C.H., K.R. and J.-S.E.; funding acquisition, E.C.H., K.R. All authors have read and agreed to the published version of the manuscript.

Funding

E.C.H. is supported by an ARC Australian Laureate Fellowship (FL170100022). The Taronga Conservation Science Initiative contributed to this Pathogen Discovery Program.

Acknowledgments

This research was supported by Taronga Conservation Society Australia and the Taronga Foundation. We thank Wendy Ruscoe, Jacqui Richardson, and Arthur Georges (University of Canberra) for their contribution of specimens, history, and context. Paul Thompson (Taronga Wildlife Hospital) is acknowledged for his role in microbial culture. We acknowledge the University of Sydney HPC service (Artemis) for providing high-performance computing resources that contributed to the research reported in this paper. We acknowledge Chris Jolly (University of New South Wales) for his advice for the manuscript and thank Wen-Lin Wang (National Taiwan University) for her advice and support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tingley, R.; Macdonald, S.L.; Mitchell, N.J.; Woinarski, J.C.Z.; Meiri, S.; Bowles, P.; Cox, N.A.; Shea, G.M.; Böhm, M.; Chanson, J.; et al. Geographic and taxonomic patterns of extinction risk in Australian squamates. Biol. Conserv. 2019, 238, 108203. [Google Scholar] [CrossRef]
  2. Cogger, H.G. Reptiles and Amphibians of Australia, 7th ed.; CSIRO Publishing: Collingwood, VIC, Australia, 2014; 1033p. [Google Scholar]
  3. Deveson, I.W.; Holleley, C.E.; Blackburn, J.; Marshall Graves, J.A.; Mattick, J.S.; Waters, P.D.; Georges, A. Differential intron retention in Jumonji chromatin modifier genes is implicated in reptile temperature-dependent sex determination. Sci. Adv. 2017, 3, e1700731. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Holleley, C.E.; O’Meally, D.; Sarre, S.D.; Marshall Graves, J.A.; Ezaz, T.; Matsubara, K.; Azad, B.; Zhang, X.; Georges, A. Sex reversal triggers the rapid transition from genetic to temperature-dependent sex. Nature 2015, 523, 79–82. [Google Scholar] [CrossRef] [PubMed]
  5. Stöhr, A.C.; Blahak, S.; Heckers, K.O.; Wiechert, J.; Behncke, H.; Mathes, K.A.; Günther, P.; Zwart, P.; Ball, I.; Rüschoff, B.; et al. Ranavirus infections associated with skin lesions in lizards. Vet. Res. 2013, 44, 84. [Google Scholar] [CrossRef] [Green Version]
  6. Tamukai, K.; Tokiwa, T.; Kobayashi, H.; Une, Y. Ranavirus in an outbreak of dermatophilosis in captive inland bearded dragons (Pogona vitticeps). Vet. Dermatol. 2016, 27, 99-e28. [Google Scholar] [CrossRef]
  7. Grosset, C.; Wellehan, J.F., Jr.; Owens, S.D.; McGraw, S.; Gaffney, P.M.; Foley, J.; Childress, A.L.; Yun, S.; Malm, K.; Groff, J.M.; et al. Intraerythrocytic iridovirus in central bearded dragons (Pogona vitticeps). J. Vet. Diagn. Investig. 2014, 26, 354–364. [Google Scholar] [CrossRef] [Green Version]
  8. Papp, T.; Marschang, R.E. Detection and characterization of invertebrate iridoviruses found in reptiles and prey insects in Europe over the past two decades. Viruses 2019, 11, 600. [Google Scholar] [CrossRef] [Green Version]
  9. Benge, S.L.; Hyndman, T.H.; Funk, R.S.; Marschang, R.E.; Schneider, R.; Childress, A.L.; Wellehan, J.F. Identification of helodermatid adenovirus 2 in a captive central bearded dragon (Pogona Vitticeps), wild gila monsters (Heloderma Suspectum), and a death adder (Acanthophis Antarcticus). J. Zoo Wildl. Med. 2019, 50, 238–242. [Google Scholar]
  10. Infectious Diseases and Pathology of Reptiles: Color Atlas and Text, 2nd ed.; Jacobson, E.R.; Garner Michael, M. (Eds.) CRC Press/Taylor & Francis: Boca Raton, FL, USA, 2020; 1030p. [Google Scholar]
  11. Pénzes, J.J.; Pham, H.T.; Benkő, M.; Tijssen, P.; Benk, M. Novel parvoviruses in reptiles and genome sequence of a lizard parvovirus shed light on Dependoparvovirus genus evolution. J. Gen. Virol. 2015, 96, 2769–2779. [Google Scholar] [CrossRef]
  12. Kim, D.Y.; Mitchell, M.A.; Bauer, R.W.; Poston, R.; Cho, D.Y. An outbreak of adenoviral infection in inland bearded dragons (Pogona vitticeps) coinfected with dependovirus and coccidial protozoa (Isospora sp.). J. Vet. Diagn. Investig. 2002, 14, 332–334. [Google Scholar] [CrossRef] [Green Version]
  13. Jacobson, E.R.; Kopit, W.; Kennedy, F.A.; Funk, R.S. Coinfection of a bearded dragon, Pogona vitticeps, with adenovirus- and dependovirus-like Viruses. Vet. Pathol. 1996, 33, 343–346. [Google Scholar] [CrossRef] [PubMed]
  14. Crossland, N.A.; DiGeronimo, P.M.; Sokolova, Y.; Childress, A.L.; Wellehan, J.F.X.; Nevarez, J.; Paulsen, D. Pneumonia in a captive central bearded dragon with concurrent detection of helodermatid adenovirus 2 and a novel mycoplasma species. Vet. Pathol. 2018, 55, 900–904. [Google Scholar] [CrossRef] [Green Version]
  15. Abbas, M.D.; Ball, I.; Ruckova, Z.; Öfner, S.; Stöhr, A.C.; Marschang, R.E. Virological screening of bearded dragons (Pogona vitticeps) and the first detection of paramyxoviruses in this species. J. Herpetol. Med. Surg. 2012, 22, 86–90. [Google Scholar] [CrossRef]
  16. Capraro, A.; O’Meally, D.; Waters, S.A.; Patel, H.R.; Georges, A.; Waters, P.D. Waking the sleeping dragon: Gene expression profiling reveals adaptive strategies of the hibernating reptile Pogona vitticeps. BMC Genom. 2019, 20, 460. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Drake, K.K.; Aiello, C.M.; Bowen, L.; Lewison, R.L.; Esque, T.C.; Nussear, K.E.; Waters, S.C.; Hudson, P.J. Complex immune responses and molecular reactions to pathogens and disease in a desert reptile (Gopherus agassizii). Ecol. Evol. 2019, 9, 2516–2534. [Google Scholar] [CrossRef] [Green Version]
  18. Hyndman, T.H.; Johnson, R.S. Evidence for the vertical transmission of Sunshine virus. Vet. Microbiol. 2015, 175, 179–184. [Google Scholar] [CrossRef]
  19. Wellehan, J.F., Jr.; Childress, A.L.; Marschang, R.E.; Johnson, A.J.; Lamirande, E.W.; Roberts, J.F.; Vickers, M.L.; Gaskin, J.M.; Jacobson, E.R. Consensus nested PCR amplification and sequencing of diverse reptilian, avian, and mammalian orthoreoviruses. Vet. Microbiol. 2009, 133, 34–42. [Google Scholar] [CrossRef]
  20. Ahne, W.; Batts, W.N.; Kurath, G.; Winton, J.R. Comparative sequence analyses of sixteen reptilian paramyxoviruses. Virus Res. 1999, 63, 65–74. [Google Scholar] [CrossRef]
  21. Hyndman, T.H.; Marschang, R.E.; Wellehan, J.F.X., Jr.; Nicholls, P.K. Isolation and molecular identification of Sunshine virus, a novel paramyxovirus found in Australian snakes. Infect. Genet. Evol. 2012, 12, 1436–1446. [Google Scholar] [CrossRef] [Green Version]
  22. Shi, M.; Lin, X.-D.; Chen, X.; Tian, J.-H.; Chen, L.-J.; Li, K.; Wang, W.; Eden, J.-S.; Shen, J.-J.; Liu, L.; et al. The evolutionary history of vertebrate RNA viruses. Nature 2018, 556, 197–202. [Google Scholar] [CrossRef]
  23. Shi, M.; Lin, X.-D.; Vasilakis, N.; Tian, J.-H.; Li, C.-X.; Chen, L.-J.; Eastwood, G.; Diao, X.-N.; Chen, M.-H.; Chen, X.; et al. Divergent viruses discovered in arthropods and vertebrates revise the evolutionary history of the Flaviviridae and related viruses. J. Virol. 2015, 90, 659–669. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M.; et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protocol. 2013, 8, 1494–1512. [Google Scholar] [CrossRef] [PubMed]
  26. Mount, D.W. Using the Basic Local Alignment Search Tool (BLAST). CSH Protoc. 2007. [Google Scholar] [CrossRef]
  27. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef]
  28. Marchler-Bauer, A.; Derbyshire, M.K.; Gonzales, N.R.; Lu, S.; Chitsaz, F.; Geer, L.Y.; Geer, R.C.; He, J.; Gwadz, M.; Hurwitz, D.I.; et al. CDD: NCBI’s conserved domain database. Nucleic Acids Res. 2015, 43, 222–226. [Google Scholar] [CrossRef] [Green Version]
  29. Bushnell, B. BBMap Short-Read Aligner, And Other Bioinformatics Tools. 2016. Available online: https://sourceforge.net/projects/bbmap/ (accessed on 28 May 2020).
  30. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  31. Nguyen, L.T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef]
  32. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef]
  33. Fahsbender, E.; Altan, E.; Seguin, M.A.; Young, P.; Estrada, M.; Leutenegger, C.M.; Delwart, E. Chapparvovirus DNA found in 4% of dogs with diarrhea. Viruses 2019, 11, 398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Li, Y.; Gordon, E.; Idle, A.; Altan, E.; Seguin, M.A.; Estrada, M.; Deng, X.; Delwart, E. Virome of a feline outbreak of diarrhea and vomiting includes bocaviruses and a novel chapparvovirus. Viruses 2020, 12, 506. [Google Scholar] [CrossRef] [PubMed]
  35. Yang, S.; Liu, Z.; Wang, Y.; Li, W.; Fu, X.; Lin, Y.; Shen, Q.; Wang, X.; Wang, H.; Zhang, W. A novel rodent Chapparvovirus in feces of wild rats. Virol. J. 2016, 13, 133. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Du, J.; Wang, W.; Chan, J.F.; Wang, G.; Huang, Y.; Yi, Y.; Zhu, Z.; Peng, R.; Hu, X.; Wu, Y.; et al. Identification of a novel ichthyic parvovirus in marine species in Hainan Island, China. Front. Microbiol. 2019, 10, 2815. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Lee, Q.; Padula, M.P.; Pinello, N.; Williams, S.H.; O’Rourke, M.B.; Fumagalli, M.J.; Orkin, J.D.; Song, R.; Shaban, B.; Brenner, O.; et al. Murine and related chapparvoviruses are nephro-tropic and produce novel accessory proteins in infected kidneys. PLoS Pathog. 2020, 16, e1008262. [Google Scholar] [CrossRef] [Green Version]
  38. Souza, W.M.; Romeiro, M.F.; Fumagalli, M.J.; Modha, S.; de Araujo, J.; Queiroz, L.H.; Durigon, E.L.; Figueiredo, L.T.M.; Murcia, P.R.; Gifford, R.J. Chapparvoviruses occur in at least three vertebrate classes and have a broad biogeographic distribution. J. Gen. Virol. 2017, 98, 225–229. [Google Scholar] [CrossRef]
  39. Kapoor, A.; Dubovi, E.J.; Henriquez-Rivera, J.A.; Lipkin, W.I. Complete genome sequence of the first canine circovirus. J. Virol. 2012, 86, 7018. [Google Scholar] [CrossRef] [Green Version]
  40. Raidal, S.R.; Sarker, S.; Peters, A. Review of psittacine beak and feather disease and its effect on Australian endangered species. Aust. Vet. J. 2015, 93, 466–470. [Google Scholar] [CrossRef]
  41. Allan, G.; Krakowka, S.; Ellis, J.; Charreyre, C. Discovery and evolving history of two genetically related but phenotypically different viruses, porcine circoviruses 1 and 2. Virus Res. 2012, 164, 4–9. [Google Scholar] [CrossRef]
  42. Lian, H.; Liu, Y.; Li, N.; Wang, Y.; Zhang, S.; Hu, R. Novel circovirus from mink, China. Emerg. Infect. Dis. 2014, 20, 1548–1550. [Google Scholar] [CrossRef]
  43. Altan, E.; Kubiski, S.V.; Burchell, J.; Bicknese, E.; Deng, X.; Delwart, E. The first reptilian circovirus identified infects gut and liver tissues of black-headed pythons. Vet. Res. 2019, 50, 35. [Google Scholar] [CrossRef] [Green Version]
  44. Tarján, Z.L.; Pénzes, J.J.; Tóth, R.P.; Benkő, M. First detection of circovirus-like sequences in amphibians and novel putative circoviruses in fishes. Acta Vet. Hung. 2014, 62, 134–144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Marton, S.; Ihász, K.; Lengyel, G.; Farkas, S.L.; Dán, Á.; Paulus, P.; Banyai, K.; Fehér, E. Ubiquiter circovirus sequences raise challenges in laboratory diagnosis: The case of honey bee and bee mite, reptiles, and free living amoebae. Acta Microbiol. Immunol. Hung. 2015, 62, 57–73. [Google Scholar] [PubMed] [Green Version]
  46. Gilbert, C.; Meik, J.M.; Dashevsky, D.; Card, D.C.; Castoe, T.A.; Schaack, S. Endogenous hepadnaviruses, bornaviruses and circoviruses in snakes. Proc. Biol. Sci. 2014, 281, 20141122. [Google Scholar] [CrossRef] [PubMed]
  47. Dunham, S.; Daly, J.M. Canine parvovirus study. Vet Rec 2016, 179, 51. [Google Scholar] [CrossRef] [PubMed]
  48. Feng, B.; Xie, Z.; Deng, X.; Xie, L.; Xie, Z.; Huang, L.; Fan, Q.; Luo, S.; Huang, J.; Zhang, Y.; et al. Genetic and phylogenetic analysis of a novel parvovirus isolated from chickens in Guangxi, China. Arch. Virol. 2016, 161, 3285–3289. [Google Scholar] [CrossRef]
  49. Zhang, J.; Liu, P.; Wu, Y.; Wang, M.; Jia, R.; Zhu, D.; Liu, M.; Yang, Q.; Wu, Y.; Zhao, X.-X.; et al. Growth characteristics of the novel goose parvovirus SD15 strain in vitro. BMC Vet. Res. 2019, 15, 63. [Google Scholar] [CrossRef]
  50. Palinski, R.M.; Mitra, N.; Hause, B.M. Discovery of a novel Parvovirinae virus, porcine parvovirus 7, by metagenomic sequencing of porcine rectal swabs. Virus Genes 2016, 52, 564–567. [Google Scholar] [CrossRef]
  51. Miłek, D.; Woźniak, A.; Stadejek, T. The detection and genetic diversity of novel porcine parvovirus 7 (PPV7) on Polish pig farms. Res. Vet. Sci. 2018, 120, 28–32. [Google Scholar] [CrossRef]
  52. Fahsbender, E.; Charlys da-Costa, A.; Elise Gill, D.; Augusto de Padua Milagres, F.; Brustulin, R.; Julio Costa Monteiro, F.; Octavio da Silva Rego, M.; Soares D’Athaide Ribeiro, E.; Cerdeira Sabino, E.; Delwart, E. Plasma virome of 781 Brazilians with unexplained symptoms of arbovirus infection include a novel parvovirus and densovirus. PLoS ONE 2020, 15, e0229993. [Google Scholar] [CrossRef] [Green Version]
  53. Farkas, S.L.; Zádori, Z.; Benkő, M.; Essbauer, S.; Harrach, B.; Tijssen, P. A parvovirus isolated from royal python (Python regius) is a member of the genus Dependovirus. J. Gen. Virol. 2004, 85, 555–561. [Google Scholar] [CrossRef]
  54. Wozniak, E.J.; DeNardo, D.F.; Brewer, A.; Wong, V.; Tarara, R.P. Identification of adenovirus- and dependovirus-like agents in an outbreak of fatal gastroenteritis in captive born California mountain kingsnakes, Lampropeltis zonata multicincta. J. Herpetol. Med. Surg. 2000, 10, 4–7. [Google Scholar] [CrossRef]
  55. Ahne, W.; Scheinert, P. Reptilian viruses: Isolation of parvovirus-like particles from corn snake Elapha guttata (Colubridae). J. Vet. Med. 1989, 36, 409–412. [Google Scholar] [CrossRef] [PubMed]
  56. Jacobson, E.R.; Gaskin, J.M.; Gardiner, C.H. Adenovirus-like infection in a boa constrictor. J. Am. Vet. Med. Assoc. 1985, 187, 1226–1227. [Google Scholar] [PubMed]
  57. Farkas, S.G.J. First Hungarian report of inclusion body hepatitis associated with adenoviruses and secondary parvovirus infection in an Indonesian pit-viper [Parias (Trimeresurus) hageni]. Magy. Állatorvosok Lapja 2008, 130, 755–761. [Google Scholar]
  58. Pénzes, J.J.; Benkő, M. Novel parvovirus from the worm lizard Trogonophis wiegmanni—First virus ever detected in amphisbaenian hosts. Acta Vet. Hung. 2014, 62, 284–292. [Google Scholar] [CrossRef] [Green Version]
  59. Pénzes, J.J.; Söderlund-Venermo, M.; Canuti, M.; Eis-Hübinger, A.M.; Hughes, J.; Cotmore, S.F.; Harrach, B. Reorganizing the family Parvoviridae: A revised taxonomy independent of the canonical approach based on host association. Arch. Virol. 2020, 165, 2133–2146. [Google Scholar] [CrossRef]
  60. Pénzes, J.J.; de Souza, W.M.; Agbandje-McKenna, M.; Gifford, R.J. An ancient lineage of highly divergent parvoviruses infects both vertebrate and invertebrate hosts. Viruses 2019, 11, 525. [Google Scholar] [CrossRef] [Green Version]
  61. Ellis, J.A.; Bratanich, A.; Clark, E.G.; Allan, G.; Meehan, B.; Haines, D.M.; Harding, J.; West, K.H.; Krakowka, S.; Konoby, C.; et al. Coinfection by porcine circoviruses and porcine parvovirus in pigs with naturally acquired postweaning multisystemic wasting syndrome. J. Vet. Diagn. Invest. 2000, 12, 21–27. [Google Scholar] [CrossRef]
  62. Thaiwong, T.; Wise, A.G.; Maes, R.K.; Mullaney, T.; Kiupel, M. Canine circovirus 1 (CaCV-1) and canine parvovirus 2 (CPV-2): Recurrent dual infections in a papillon breeding colony. Vet. Pathol. 2016, 53, 1204–1209. [Google Scholar] [CrossRef] [Green Version]
  63. Chang, W.-S.; Eden, J.-S.; Hall, J.; Shi, M.; Rose, K.; Holmes, E.C. Meta-transcriptomic analysis of virus diversity in urban wild birds with paretic disease. J. Virol. 2020, 94, e00606-20. [Google Scholar] [CrossRef]
  64. Chong, R.; Shi, M.; Grueber, C.E.; Holmes, E.C.; Hogg, C.J.; Belov, K.; Barrs, V.R. Fecal viral diversity of captive and wild Tasmanian devils characterized using virion-enriched metagenomics and metatranscriptomics. J. Virol. 2019, 93, e00205-19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Viruses 12 01073 g001 550
Figure 1. Photomicrographs of pulmonary lesions in central bearded dragons from the 2014 and 2020 outbreak. Thickening of the normally squamous infundibular epithelium. (Haematoxylin and eosin (HE) stain 40×) (A) and the segmental nature of this change observed in an affected animal from the 2020 outbreak (HE 20×) (B) Metaplasia of the normally squamous infundibular epithelium to pseudostratified columnar and multifocally ciliate. Normal lung appears on the right and infundibular epithelium becomes proliferative on the left (HE 1000×) (C) and multifocal epithelial cells with large, irregularly shaped, basophilic nuclei (arrowheads, HE 1000×) (D) Bizarre atrial epithelial cells containing eosinophilic cytoplasmic inclusions (arrowheads, (HE 1000×) (E).
Figure 1. Photomicrographs of pulmonary lesions in central bearded dragons from the 2014 and 2020 outbreak. Thickening of the normally squamous infundibular epithelium. (Haematoxylin and eosin (HE) stain 40×) (A) and the segmental nature of this change observed in an affected animal from the 2020 outbreak (HE 20×) (B) Metaplasia of the normally squamous infundibular epithelium to pseudostratified columnar and multifocally ciliate. Normal lung appears on the right and infundibular epithelium becomes proliferative on the left (HE 1000×) (C) and multifocal epithelial cells with large, irregularly shaped, basophilic nuclei (arrowheads, HE 1000×) (D) Bizarre atrial epithelial cells containing eosinophilic cytoplasmic inclusions (arrowheads, (HE 1000×) (E).
Viruses 12 01073 g001
Viruses 12 01073 g002 550
Figure 2. Genome characterization and phylogeny of bearded dragon circovirus. (A) Genome organization of bearded dragon circovirus. The outermost circles of the metadata ring represent the GC content (%, green), nucleotide polymorphism (%, orange), and read mapping coverage (yellow) of the genome. In the innermost circles, the proteins encoded by the replicase (Rep) and capsid protein (Cap) are labelled in pink and cyan, respectively. (B) Phylogenetic analysis of the Rep gene (nonstructural replicase protein) of circoviruses, including members of the Cyclovirus and Circovirus genera. Bootstrap values > 70% are presented for key nodes (1000 replicates). The tree was midpoint rooted for clarity only. Scale bar shows the number of substitutions per site. (C) Comparison of nonamer sequences with closely related viruses in the Circoviridae.
Figure 2. Genome characterization and phylogeny of bearded dragon circovirus. (A) Genome organization of bearded dragon circovirus. The outermost circles of the metadata ring represent the GC content (%, green), nucleotide polymorphism (%, orange), and read mapping coverage (yellow) of the genome. In the innermost circles, the proteins encoded by the replicase (Rep) and capsid protein (Cap) are labelled in pink and cyan, respectively. (B) Phylogenetic analysis of the Rep gene (nonstructural replicase protein) of circoviruses, including members of the Cyclovirus and Circovirus genera. Bootstrap values > 70% are presented for key nodes (1000 replicates). The tree was midpoint rooted for clarity only. Scale bar shows the number of substitutions per site. (C) Comparison of nonamer sequences with closely related viruses in the Circoviridae.
Viruses 12 01073 g002
Viruses 12 01073 g003 550
Figure 3. Genome of bearded dragon chaphamaparvovirus in comparison to related viruses. The amino acid sequence size of the ORFs of each virus is shown. Yellow and orange boxes refer to the NS and VP genes, respectively. The classic motif features HXH (blue asterisk) and GPXNTGKX (green asterisk) on NS genes are labelled.
Figure 3. Genome of bearded dragon chaphamaparvovirus in comparison to related viruses. The amino acid sequence size of the ORFs of each virus is shown. Yellow and orange boxes refer to the NS and VP genes, respectively. The classic motif features HXH (blue asterisk) and GPXNTGKX (green asterisk) on NS genes are labelled.
Viruses 12 01073 g003
Viruses 12 01073 g004 550
Figure 4. Genome characterization and phylogenetic relationships of bearded dragon chaphamaparvovirus. (A) Genome depiction and reads mapping coverage of bearded dragon chaphamaparvovirus. (B) Phylogenetic analysis of the nonstructural protein gene (NS gene) of parvoviruses, including other chaphamaparvoviruses and members of the Parvovirinae, Densovirinae, and Hamaparvovirine. Bootstrap values >70% were presented for key nodes (1000 replicates). The tree was midpoint rooted for clarity only. Scale bar shows the number of substitutions per site.
Figure 4. Genome characterization and phylogenetic relationships of bearded dragon chaphamaparvovirus. (A) Genome depiction and reads mapping coverage of bearded dragon chaphamaparvovirus. (B) Phylogenetic analysis of the nonstructural protein gene (NS gene) of parvoviruses, including other chaphamaparvoviruses and members of the Parvovirinae, Densovirinae, and Hamaparvovirine. Bootstrap values >70% were presented for key nodes (1000 replicates). The tree was midpoint rooted for clarity only. Scale bar shows the number of substitutions per site.
Viruses 12 01073 g004
Table
Table 1. PCR primers used in this study.
Table 1. PCR primers used in this study.
Primer NameTargetSequence 5′-3′Site StartSite EndTm (°C)GC
(%)
BDCV_F1BDCV_capGTATTGTTTCAGAGATTCGCCG658657.845.5
BDCV_R1ATTGTAATGAAGGCGTTGTTCC1144116557.640.9
BDCV_F2BDCV_capGCCGAAGTGTCTGGTAATACTA97499557.343.5
BDCV_R2GATTGGTTAAGAGTCGGGATTTC44046257.640.9
BDCV_F3BDCV_repAATACCTGACCAGTCTGCCA1551157057.350.0
BDCV_R3CGATGGCGTAGGAACAGATG65964059.355.0
BDCV_F4BDCV_repGGTGAAGTCGATTGTTCCTTTG31233358.445.5
BDCV_R4GTGGGTATTGTTTCAGAGATTCG1633161158.943.5
BDchPV_F1BDchPV_NSCTGAACAATATCTCTTCGTTGGAG21724057.441.7
BDchPV_R1TCCCAGACACCCATATAACAATC1490151257.743.5
BDchPV_F2BDchPV_NSGCTCTCCTGAATATTGTGATATGTC95497857.440.0
BDchPV_R2TTCCCTACCCACTTCTTTCTG2255227556.947.6
BDchPV_F3BDchPV_VPCTAGCTCAGACATCAGATCCAG2004202557.550.0
BDchPV_R3GTCGGTGGAGGATTGTCTATAG3346336757.550.0
BDchPV_F4BDchPV_VPAGCTTTTCCCAAATGACACTG2848286857.042.9
BDchPV_R4CAGAAGACACAACAGTAGGTAAAG4087411057.241.7
BDCV: Bearded Dragon CircoVirus, BDchPV: Bearded Dragon ChaphamaParvoVirus, rep: replicase, cap: capsid, NS: Non-Structural protein, VP: Viral capsid Protein.
Table
Table 2. Rep protein (amino acid) motifs detected in the bearded dragon circovirus.
Table 2. Rep protein (amino acid) motifs detected in the bearded dragon circovirus.
Name of Region Sequence in the CircoviridaeSequence in Bearded Dragon Circovirus
RCR motifRCR motif IFT(L/I)NNFTINN
RCR motif IIPHLQGPHIQG
RCR motif IIIYC(S/x)KYCSK
SF3 helicase motifWalker AG(P/x)(P/x) GxGK(S/t)GEPGSGKS
Walker BuuDDFILDDF
Walker CuTSNITSN
Table
Table 3. PCR testing of captive reptiles from the 2014 and 2020 outbreaks.
Table 3. PCR testing of captive reptiles from the 2014 and 2020 outbreaks.
YearCase no.SpeciesClinical SignOrganPCR Results
BDCVBDchPV
201410043.1CBDRespiratoryKidney--
Lung--
10043.2CBDRespiratoryLiver--
10043.3CBDRespiratoryKidney+-
Lung++
10043.4CBDRespiratoryLiver-+
10043.5CBDRespiratoryLiver+-
202013270.1CBDMild neurologicalBrain--
Liver-+
13270.2CBDNormalBrain--
Liver--
13270.3CBDMild nonspecific illnessBrain--
Liver++
13270.4CBDNeurologicalBrain-+
Liver++
13270.5CBDNormalBrain--
Liver--
13270.6CBDNeurologicalBrainn.a.n.a.
Liver--
13270.7CBDNeurologicalBrainn.a.n.a.
Liver--
13270.8CBDNeurologicalBrainn.a.n.a.
Liver--
13270.9JLNormalBrain--
Liver+-
13270.10JLNeurologicalBrain--
Liver--
13270.11JLNeurologicalBrain--
Liver--
13270.12GEDMoribundBrain--
Liver--
Total: BDCV positive cases: 5/17; BDchPV positive cases: 5/17
CBD: central bearded dragon (Pogona vitticeps); JL: jacky lizard (Amphibolurus muricatus); GED: grassland earless dragon (Tympanocryptis pinguicolla); n.a.: nonapplicable; green: negative; red: positive.

Share and Cite

MDPI and ACS Style

Chang, W.-S.; Li, C.-X.; Hall, J.; Eden, J.-S.; Hyndman, T.H.; Holmes, E.C.; Rose, K. Meta-Transcriptomic Discovery of a Divergent Circovirus and a Chaphamaparvovirus in Captive Reptiles with Proliferative Respiratory Syndrome. Viruses 2020, 12, 1073. https://0-doi-org.brum.beds.ac.uk/10.3390/v12101073

AMA Style

Chang W-S, Li C-X, Hall J, Eden J-S, Hyndman TH, Holmes EC, Rose K. Meta-Transcriptomic Discovery of a Divergent Circovirus and a Chaphamaparvovirus in Captive Reptiles with Proliferative Respiratory Syndrome. Viruses. 2020; 12(10):1073. https://0-doi-org.brum.beds.ac.uk/10.3390/v12101073

Chicago/Turabian Style

Chang, Wei-Shan, Ci-Xiu Li, Jane Hall, John-Sebastian Eden, Timothy H. Hyndman, Edward C. Holmes, and Karrie Rose. 2020. "Meta-Transcriptomic Discovery of a Divergent Circovirus and a Chaphamaparvovirus in Captive Reptiles with Proliferative Respiratory Syndrome" Viruses 12, no. 10: 1073. https://0-doi-org.brum.beds.ac.uk/10.3390/v12101073

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop