Next Article in Journal
Avian Influenza Virus Subtype H9N2 Affects Intestinal Microbiota, Barrier Structure Injury, and Inflammatory Intestinal Disease in the Chicken Ileum
Next Article in Special Issue
The First Isolation and Whole Genome Sequencing of Murray Valley Encephalitis Virus from Cerebrospinal Fluid of a Patient with Encephalitis
Previous Article in Journal
Paradoxical Effect of Chloroquine Treatment in Enhancing Chikungunya Virus Infection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

A Novel Hepadnavirus Identified in an Immunocompromised Domestic Cat in Australia

1
Sydney School of Veterinary Science, Faculty of Science, University of Sydney, Sydney, NSW 2006, Australia
2
School of Life and Environmental Sciences and Sydney Medical School, Charles Perkins Centre, University of Sydney, Sydney, NSW 2006, Australia
3
Marie Bashir Institute for Infectious Diseases and Biosecurity, The University of Sydney, Sydney, NSW 2006, Australia
4
School of Animal and Veterinary Sciences, Faculty of Science, The University of Adelaide, Adelaide, Roseworthy, SA 5371, Australia
5
Vets@Acacia Gardens, Quakers Hill, NSW 2006, Australia
*
Author to whom correspondence should be addressed.
The authors contributed equally to this article.
Submission received: 27 April 2018 / Revised: 10 May 2018 / Accepted: 14 May 2018 / Published: 17 May 2018
(This article belongs to the Special Issue Emerging Viruses)

Abstract

:
High-throughput transcriptome sequencing allows for the unbiased detection of viruses in host tissues. The application of this technique to immunosuppressed animals facilitates the detection of viruses that might otherwise be excluded or contained in immunocompetent individuals. To identify potential viral pathogens infecting domestic cats we performed high-throughput transcriptome sequencing of tissues from cats infected with feline immunodeficiency virus (FIV). A novel member of the Hepadnaviridae, tentatively named domestic cat hepadnavirus, was discovered in a lymphoma sample and its complete 3187 bp genome characterized. Phylogenetic analysis placed the domestic cat hepadnavirus as a divergent member of mammalian orthohepadnaviruses that exhibits no close relationship to any other virus. DNA extracted from whole blood from pet cats was positive for the novel hepadnavirus by PCR in 6 of 60 (10%) FIV-infected cats and 2 of 63 (3.2%) FIV-uninfected cats. The higher prevalence of hepadnavirus viraemia detected in FIV-infected cats mirrors that seen in human immunodeficiency virus-infected humans coinfected with hepatitis B virus. In summary, we report the first hepadnavirus infection in a carnivore and the first in a companion animal. The natural history, epidemiology and pathogenic potential of domestic cat hepadnavirus merits additional investigation.

Viruses of the genus Orthohepadnavirus, family Hepadnaviridae, are partially double-stranded DNA viruses that infect a variety of mammals including primates, bats and rodents. The type species, hepatitis B virus (HBV), is a major public health problem; chronic HBV infection increases the risk of liver diseases including cirrhosis and hepatocellular carcinoma, and an estimated 257 million people are currently living with HBV [1].
We identified a new hepadnavirus in a domestic cat from Australia. In 2016, a seven year-old male-neutered domestic shorthair cat was presented for vomiting and weight loss. On physical examination a large mid-abdominal mass was palpated. Fine needle aspirate cytology of the mass was consistent with large cell lymphoma. The owners requested euthanasia and necropsy and consented to tissue collection approved by the University of Sydney animal ethics committee, (2014/626). Multicentric, large cell, high-grade, B-cell lymphoma was diagnosed on immunohistopathology of formalin-fixed paraffin embedded tissue. Infection with the immunosuppressive lentivirus, feline immunodeficiency virus (FIV), was diagnosed by serology for anti-FIV antibody (Witness FIV/FeLV, Zoetis) and PCR of tumour DNA confirmed by sequencing as described previously [2]. Lymphoma tissue from this case was utilized in a virus discovery project.
Total RNA extracted from frozen lymphoma was enriched for non-ribosomal RNA and Illumina TruSeq paired-end libraries were prepared. Briefly, 30 mg of lymphoma was homogenised using 1.4 mm ceramic beads (OMNI Homogeniser, Kennesaw, GA, USA), total RNA was extracted using the Qiagen RNeasy mini kit (Qiagen, Hilden, Germany) followed by on-column DNA removal to remove genomic DNA. RNA quality, examined using an Agilent Bioanalyzer 2100 (Agilent Technologies, Mulgrave, VIC, Australia), demonstrated RNA integrity numbers greater than eight. Library preparation was performed using a TruSeq RNA library preparation kit (Illumina, San Diego, CA, USA) and ribosomal RNA was removed using a Ribo-Zero gold epidemiology kit (Illumina, USA). RNA sequencing of 100 bp pair-end libraries on an Illumina HiSeq 2500 platform (San Diego, CA, USA) yielded 321,297,408 total reads. Sequencing reads were uploaded onto the NCBI Sequence Read Archive (SRA) database under the BioProject submission SUB3718346.
To identify potential viral transcripts, sequencing reads were assembled de novo using the Trinity software (version 2.1.1) [3]. The resulting contigs were then blasted against the non-redundant nucleotide database (nt) using blastn program with an E-value threshold of 1 × 10−10 and against the non-redundant protein database (nr) using Diamond version 0.8.28.90 with an E-value threshold of 1 × 10−4 [4]. To discover viruses at low abundance, we performed Diamond blastx analyses on sequence reads against the nr database. The virus-associated contigs and reads were recovered based on the taxonomic annotation of each blast hit. Each positive hit was subsequently confirmed by manually inspecting the nucleotide or amino acid sequence alignment with their closest relatives. To exclude the potential presence of endogenous virus (EVEs) in the blast results, we also blasted the reads/contigs against the cat genome sequences available at the reference genome database and the whole genome shotgun database. In total, eight HBV-like reads were identified that shared 73–94% amino acid identity with the core protein, surface protein, and polymerase of known hepadnaviruses. This suggested the presence of a novel virus that was divergent from currently known hepadnaviruses.
To confirm the presence of the hepadnavirus, DNA was extracted from frozen lymphoma and examined with a panhepadnavirus PCR assay targeting the polymerase ORF [5]. A 250 bp product was obtained and sequenced, thereby validating the RNA sequencing. To amplify the entire genome of the novel hepadnavirus from lymphoma tissue, a five primer-set PCR was designed based on recovered reads (Table 1). PCR reactions were performed using Bioline MyTaq HotStart polymerase (Bioline, Australia) based on the manufacturer’s recommended protocol, using 200–500 ng of template DNA, 1× MyTaq buffer, a final primer concentration of 0.4 µM in a 50 µL reaction and an annealing temperature ranging from 55 to 58 °C. For primer sets Hgap-F/R and Cir5-F/Cir4-R a final primer concentration of 0.3 µM was used. The amplicons were sequenced (Macrogen, Seoul, South Korea) and analysed using Geneious version 11.2.0 and the entire 3187 bp circular genome was assembled (Figure 1). The genome contained the polymerase (P), surface (S), core (C) and X ORFs typical of orthohepadnaviruses. Notably, the complete P protein exhibited between 63.2% to 68.7% amino acid similarity to known orthohepadnaviruses. Such a large genetic distance merits assignment of a new species within the genus Orthohepadnavirus, which we tentatively named domestic cat hepadnavirus. The complete genome sequence of domestic cat hepadnavirus has been deposited in GenBank, accession number MH307930.
To determine the evolutionary relationships of domestic cat hepadnavirus, we compared it to reference sequences of P proteins from the entire Hepadnaviridae family downloaded from GenBank. To show the phylogenetic diversity of vertebrate hepadnavirusues in the P protein tree, homologous sequences representative of the viral diversity in other vertebrates [6,7] were added to the alignment. In the case of the more divergent C and S ORFs, data sets of orthohepadnaviruses only were constructed, utilizing the hepadnavirus from the bluegill (GenBank accession KX058433) as an outgroup [6]. However, because of the small size of the X ORF gene alignment, the lack of an appropriate outgroup for this ORF, and because it is largely contained within the other ORFs, the X ORF was excluded from the analysis. All amino acid sequences were aligned using MAFFT (version 7) [8], and all ambiguously aligned regions were removed using TrimAl [9]. This process resulted in final data sets of the following sizes: P protein = 30 taxa, 419 amino acids; C protein = 16 taxa, 182 amino acids; S ORF = 16 taxa, 255 amino acids.
Phylogenetic trees of all three data sets were inferred using the maximum likelihood approach available in PhyML 3.0, employing the LG + Γ model of amino acid substitution and a SPR branch-swapping algorithm [6]. In the case of the P ORF phylogeny (Figure 2A) it is striking that although domestic cat hepadnavirus clearly groups with the mammalian orthohepadnaviruses, it was the most basal member of this group, including those viruses recently identified in bats [10], and exhibited no close phylogenetic relationship to any other hepadnavirus [11]. This divergent position was also well supported by SH-like branch supports (value of 0.84) separating domestic cat hepadnavirus from the other mammalian hepadnaviruses. Less phylogenetic resolution (i.e., lower SH-like branch supports at key nodes) and differing topologies (such as the position of the tent-making bat) are apparent in the shorter C (Figure 2B) and S (Figure 2C) ORF alignments, although domestic cat hepadnavirus was consistently distinct from all other orthohepadnaviruses. Although this analysis clearly shows that domestic cat hepadnavirus is evolutionarily distinct, its exact phylogenetic position is uncertain.
Hepadnaviruses are typically hepatotropic, although viraemia is common in immunosuppressed hosts [12]. To investigate whether this immunosuppressed cat was viraemic for the novel hepadnavirus, PCR was performed on DNA extracted from autologous frozen whole blood using primer set Hgap-F/R (Table 1). A positive result on PCR was confirmed by sequencing. Next, a molecular survey for hepadnavirus viraemia was performed using DNA extracted from stored whole blood from adult pet cats in Australia (University of Sydney Animal Ethics Committee approvals N00/7-2013/3/6029 and 2016/1002). The novel hepadnavirus was detected in whole blood from 6 of 60 (10%) FIV-infected cats and 2 of 63 (3.2%) uninfected cats. Confirmation of this finding in larger studies is awaited. Nonetheless, it is noteworthy that a similar pattern is seen among HIV infected humans where HBV is a common co-pathogen [13]. For example, in USA and Europe the prevalence of persistent HBV antigenaemia is 1% in the general population compared with 7–10% among HIV-infected homosexual men who are more likely to be exposed to and less likely to clear HBV [14].
In conclusion, we report the first hepadnavirus infection in a member of the mammalian order Carnivora and the first report of a hepadnavirus in a companion species, with hepadnavirus viraemia commonly detected among FIV-infected cats. Future investigation of the natural history and epidemiology of domestic cat hepadnavirus and its impact on feline health is clearly warranted.

Author Contributions

M.A., M.S., V.R.B., E.C.H. and J.A.B. designed and conducted the study, and wrote the manuscript; S.L., B.J. and B.H. contributed essential data and reviewed the manuscript.

Funding

This research was funded by the Morris Animal Foundation (D15FE-001), a private donation to the University of Sydney and the Cat Protection Society of NSW. ECH is funded by an ARC Australian Laureate Fellowship (FL170100022).

Acknowledgments

We thank Sydney Informatics Core Research Facility and the University of Sydney HPC service at The University of Sydney for providing bioinformatics resources that contributed to the results reported in this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organisation Global Hepatitis Report. 2017. Available online: http://www.who.int/hepatitis/publications/global-hepatitis-report2017/en/ (accessed on 27 April 2018).
  2. McLuckie, A.; Barrs, V.; Lindsay, S.; Aghazadeh, M.; Sangster, C.; Beatty, J. Molecular diagnosis of felis catus gammaherpesvirus 1 (FcaGHV1) infection in cats of known retrovirus status with and without lymphoma. Viruses 2018, 10, 128. [Google Scholar] [CrossRef] [PubMed]
  3. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from rna-seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644. [Google Scholar] [CrossRef] [PubMed]
  4. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using diamond. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef] [PubMed]
  5. Wang, B.; Yang, X.L.; Li, W.; Zhu, Y.; Ge, X.Y.; Zhang, L.B.; Zhang, Y.Z.; Bock, C.T.; Shi, Z.L. Detection and genome characterization of four novel bat hepadnaviruses and a hepevirus in china. Virol. J. 2017, 14, 40. [Google Scholar] [CrossRef] [PubMed]
  6. Dill, J.A.; Camus, A.C.; Leary, J.H.; Di Giallonardo, F.; Holmes, E.C.; Ng, T.F.F. Distinct viral lineages from fish and amphibians reveal the complex evolutionary history of hepadnaviruses. J. Virol. 2016, 90, 7920–7933. [Google Scholar] [CrossRef] [PubMed]
  7. Lauber, C.; Seitz, S.; Mattei, S.; Suh, A.; Beck, J.; Herstein, J.; Borold, J.; Salzburger, W.; Kaderali, L.; Briggs, J.A.G.; et al. Deciphering the origin and evolution of hepatitis B viruses by means of a family of non-enveloped fish viruses. Cell Host Microbe 2017, 22, 387–399. [Google Scholar] [CrossRef] [PubMed]
  8. 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] [PubMed]
  9. Capella-Gutiérrez, S.; Silla-Martínez, J.M.; Gabaldón, T. Trimal: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 2009, 25, 1972–1973. [Google Scholar] [CrossRef] [PubMed]
  10. Nie, F.Y.; Lin, X.D.; Hao, Z.Y.; Chen, X.N.; Wang, Z.X.; Wang, M.R.; Wu, J.; Wang, H.W.; Zhao, G.; Ma, R.Z.; et al. Extensive diversity and evolution of hepadnaviruses in bats in China. Virology 2018, 514, 88–97. [Google Scholar] [CrossRef] [PubMed]
  11. Drexler, J.F.; Geipel, A.; Konig, A.; Corman, V.M.; van Riel, D.; Leijten, L.M.; Bremer, C.M.; Rasche, A.; Cottontail, V.M.; Maganga, G.D.; et al. Bats carry pathogenic hepadnaviruses antigenically related to hepatitis B virus and capable of infecting human hepatocytes. Proc. Natl. Acad. Sci. USA 2013, 110, 16151–16156. [Google Scholar] [CrossRef] [PubMed]
  12. ICTV. Family—Hepadnaviridae. In Virus Taxonomy; Elsevier: San Diego, CA, USA, 2012; pp. 445–455. [Google Scholar]
  13. Lisco, A.; Vanpouille, C.; Margolis, L. Coinfecting viruses as determinants of hiv disease. Curr. HIV/AIDS Rep. 2009, 6, 5–12. [Google Scholar] [CrossRef] [PubMed]
  14. McGovern, B.H. The epidemiology, natural history and prevention of hepatitis B: Implications of hiv coinfection. Antivir. Ther. 2007, 12, H3–H13. [Google Scholar] [PubMed]
Figure 1. Genome organization of the domestic cat hepadnavirus. The complete genome consists of 3187 bp. The innermost circles represent the GC (blue) and AT (green) content of the genome. The proteins encoded by the polymerase, surface, core and X ORFs are labelled, as are the positions of primers used in this study (Table 1).
Figure 1. Genome organization of the domestic cat hepadnavirus. The complete genome consists of 3187 bp. The innermost circles represent the GC (blue) and AT (green) content of the genome. The proteins encoded by the polymerase, surface, core and X ORFs are labelled, as are the positions of primers used in this study (Table 1).
Viruses 10 00269 g001
Figure 2. The Phylogenetic position of domestic cat hepadnavirus within the Hepadnaviridae. (A) Maximum likelihood phylogeny based on the P (polymerase) ORF and incorporating a wide range of vertebrate hepadnaviruses, including Roundleaf bat HBV (YP_009045991), Horseshoe bat HBV (YP_009045995), Long-fingered bat HBV (YP_007677999), Tent-making bat HBV (YP_009045999), Bat HBVs (AVW79974, ARM20233, ARM20221, and AQT40960), Ground squirrel hepatitis virus (NP_040994), Arctic ground squirrel HBV (AAB08032), Woodchuck hepatitis virus (NP_671813); Hepatitis B virus (YP_009173866), Woolly monkey HBV (YP_009175034), Capuchin monkey HBV (AVV68831), Bluegill hepadnavirus (YP_009259541), Astatotilapia metahepadnavirus (ref7), Tetra hepadnavirus (ref7), Icefish hepadnavirus (ref7), Heron HBV (NP_040998), Parrot HBVs (YP_004956864, AFY97702), Ross's goose HBV (YP_024968), Duck HBV (NP_039821), Snow goose HBV (YP_031695), Stork HBV (CAC80811), Sheldgoose HBV (YP_024974), Tinamou HBV (YP_009389524), Tibetan frog hepadnavirus (YP_009259545), White sucker HBV (YP_009165599). This tree was mid-point rooted for clarity only. Phylogenies based on the (B) C ORF (C) and S ORF of mammalian orthohepadnaviruses and rooted by the homologous sequences found in the bluegill. In all phylogenies the host groups of each virus species are indicated by symbols to the right of the tree and branch lengths are scaled according to the number of amino acid substitutions per site. SH-like support values >0.80 for nodal support are shown by asterisks in all cases.
Figure 2. The Phylogenetic position of domestic cat hepadnavirus within the Hepadnaviridae. (A) Maximum likelihood phylogeny based on the P (polymerase) ORF and incorporating a wide range of vertebrate hepadnaviruses, including Roundleaf bat HBV (YP_009045991), Horseshoe bat HBV (YP_009045995), Long-fingered bat HBV (YP_007677999), Tent-making bat HBV (YP_009045999), Bat HBVs (AVW79974, ARM20233, ARM20221, and AQT40960), Ground squirrel hepatitis virus (NP_040994), Arctic ground squirrel HBV (AAB08032), Woodchuck hepatitis virus (NP_671813); Hepatitis B virus (YP_009173866), Woolly monkey HBV (YP_009175034), Capuchin monkey HBV (AVV68831), Bluegill hepadnavirus (YP_009259541), Astatotilapia metahepadnavirus (ref7), Tetra hepadnavirus (ref7), Icefish hepadnavirus (ref7), Heron HBV (NP_040998), Parrot HBVs (YP_004956864, AFY97702), Ross's goose HBV (YP_024968), Duck HBV (NP_039821), Snow goose HBV (YP_031695), Stork HBV (CAC80811), Sheldgoose HBV (YP_024974), Tinamou HBV (YP_009389524), Tibetan frog hepadnavirus (YP_009259545), White sucker HBV (YP_009165599). This tree was mid-point rooted for clarity only. Phylogenies based on the (B) C ORF (C) and S ORF of mammalian orthohepadnaviruses and rooted by the homologous sequences found in the bluegill. In all phylogenies the host groups of each virus species are indicated by symbols to the right of the tree and branch lengths are scaled according to the number of amino acid substitutions per site. SH-like support values >0.80 for nodal support are shown by asterisks in all cases.
Viruses 10 00269 g002
Table 1. Oligonucleotides used to amplify the whole genome of domestic cat hepadnavirus.
Table 1. Oligonucleotides used to amplify the whole genome of domestic cat hepadnavirus.
OligonucleotidesSequenceSize (bp)Tm (°C)
Cir5-F5′-TTGGCACCTGGATTCGCA-3′140057
Cir4-R5′-AGATGTTCCACACTCTTAGCC-3′
Cir8-F5′-TTGGCACCTGGATTCGCA-3′90058
Cir7-R5′-CGTAGACGAAGGACACGTC-3′
Cir7-F5′-CCATCGATTTACACACTTCCCA-3′95057
Cir5-R5′-TGCGAATCCAGGTGCCAA-3′
Cir1-R5′-ATAACCGTATGCTCCGGAAG-3′100055
Hgap-F5′-GTGCTCTGATAACCGTATGCTC-3′23055
Hgap-R5′-CTAGAATGGCTACATGGGGTTAG-3′

Share and Cite

MDPI and ACS Style

Aghazadeh, M.; Shi, M.; Barrs, V.R.; McLuckie, A.J.; Lindsay, S.A.; Jameson, B.; Hampson, B.; Holmes, E.C.; Beatty, J.A. A Novel Hepadnavirus Identified in an Immunocompromised Domestic Cat in Australia. Viruses 2018, 10, 269. https://0-doi-org.brum.beds.ac.uk/10.3390/v10050269

AMA Style

Aghazadeh M, Shi M, Barrs VR, McLuckie AJ, Lindsay SA, Jameson B, Hampson B, Holmes EC, Beatty JA. A Novel Hepadnavirus Identified in an Immunocompromised Domestic Cat in Australia. Viruses. 2018; 10(5):269. https://0-doi-org.brum.beds.ac.uk/10.3390/v10050269

Chicago/Turabian Style

Aghazadeh, Mahdis, Mang Shi, Vanessa R. Barrs, Alicia J. McLuckie, Scott A. Lindsay, Barbara Jameson, Bronte Hampson, Edward C. Holmes, and Julia A. Beatty. 2018. "A Novel Hepadnavirus Identified in an Immunocompromised Domestic Cat in Australia" Viruses 10, no. 5: 269. https://0-doi-org.brum.beds.ac.uk/10.3390/v10050269

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