Next Article in Journal
Special Issue “Recent Advances in Morbillivirus Vaccine Development and Oncolytic Virotherapy”
Previous Article in Journal
Viral Pathogenesis, Recombinant Vaccines, and Oncolytic Virotherapy: Applications of the Canine Distemper Virus Reverse Genetics System
 
 
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 3
10.3390/v12030340
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:
Brief Report

Coinfections of Novel Polyomavirus, Anelloviruses and a Recombinant Strain of Myxoma Virus-MYXV-Tol Identified in Iberian Hares

by
Ana Águeda-Pinto
1,2,3,
Simona Kraberger
3,4,
Michael C. Lund
4,
Christian Gortázar
5,
Grant McFadden
3,
Arvind Varsani
4,6,* and
Pedro J. Esteves
1,2,7,*
1
CIBIO/InBio—Centro de Investigação em Biodiversidade e Recursos Genéticos, Universidade do Porto, Campus Agrário de Vairão, 4485-661 Vairão, Portugal
2
Departamento de Biologia, Faculdade de Ciências, Universidade do Porto, 4169-007 Porto, Portugal
3
Center for Immunotherapy, Vaccines, and Virotherapy (CIVV), The Biodesign Institute, Arizona State University, Tempe, AZ 85287, USA
4
The Biodesign Center for Fundamental and Applied Microbiomics, Center for Evolution and Medicine and School of Life sciences, Arizona State University, Tempe, AZ 85287, USA
5
SaBio Instituto de Investigación en Recursos Cinegéticos IREC-CSIC-UCLM-JCCM, Ronda de Toledo, 28005 Ciudad Real, Spain
6
Structural Biology Research Unit, Department of Clinical Laboratory Sciences, University of Cape Town, Cape Town 7701, South Africa
7
CITS—Centro de Investigação em Tecnologias da Saúde, IPSN, CESPU, 4585-116 Gandra, Portugal
*
Authors to whom correspondence should be addressed.
Viruses 2020, 12(3), 340; https://0-doi-org.brum.beds.ac.uk/10.3390/v12030340
Submission received: 24 February 2020 / Revised: 11 March 2020 / Accepted: 17 March 2020 / Published: 20 March 2020
(This article belongs to the Section Animal Viruses)
Browse Figures
Versions Notes

Abstract

:
Viruses are ubiquitous in nature; however, very few have been identified in the Leporid species. In the fall of 2018, an outbreak of myxomatosis in Iberian hares (Lepus granatensis) was reported in Spain and a novel recombinant myxoma virus strain (MYXV-Tol) was identified. To investigate variability within the recombinant region of the MYXV-Tol and identify any potential viral coinfections, samples (ear, eyelid or vaginal) of Iberian hares were collected from Spain and analyzed. The presence of the recombinant region of the MYXV-Tol was confirmed in six out of eleven samples analyzed. Additionally, a polyomavirus (family Polyomaviridae), representing a putative new species, and anelloviruses (family Anelloviridae) belonging to two putative species were identified, some as coinfection with the recombinant MYXV-Tol. The two polyomavirus genomes were identified in two hares and share >99% genome-wide identity. Based on the analysis of their large T-antigen, the new polyomavirus clusters in a distant clade from other mammals sharing <64% amino acid identity. A total of 14 anelloviruses were identified, which share 63–99% genome-wide identity. Overall, our results show a coinfection of different DNA viruses in the studied samples and raise awareness regarding the extensive unsampled diversity of viruses in hares.

1. Introduction

New molecular tools and sequencing technologies have revolutionized viral detection, enabling a better exploration of viral diversity within various organisms. However, most studies have focused on viruses that are strongly associated with disease. The bias toward the study of viruses associated with disease is clear in the Leporidae family (Lagomorpha order). In fact, most of the research undertaken has been directed toward the highly pathogenic lagoviruses (family Caliciviridae), such as the rabbit hemorrhagic disease virus [1,2,3], and the leporipoxviruses (family Poxviridae) like the myxoma virus (MYXV) [4,5,6,7]. In a recent study of diseased Iberian hares (Lepus granatensis) from the Toledo province in Spain, a novel recombinant MYXV strain was identified (MYXV-Tol; GenBank Accession MK836424) showing, for the first time, that a pathogenic myxoma virus could also infect and cause myxomatosis in species of the Lepus genus [8]. Historically, MYXV has evolved in natural Sylvilagus hosts, like the South American tapeti and the North American brush rabbit. Nonetheless, when MYXV encounters naïve European Oryctolagus rabbits, the virus causes the lethal myxomatosis disease. The genome of the recombinant MYXV-Tol identified in Iberian hares is ~99% similar to MYXV variants/strains previously reported circulating in European rabbits, with the exception of a recombinant region ~2800 bp in length and three disrupted genes (M009L, M036L and M152R). The recombinant region encodes four additional open reading frames (ORFs), which were more closely related, but not identical, to the MYXV proteins encoded by M060R, M061R, M064R and M065R genes [8,9,10].
Besides the rabbit hemorrhagic disease virus and MYXV, there is limited information on viruses associated with Leporids. A recent study has reported the presence of a novel herpesvirus (Herpesviridae family) in Iberian hare samples collected during the myxomatosis outbreak in Spain [11]. Metagenomics studies targeting hare fecal samples found circular replication-associated proteins encoding single-stranded (CRESS) DNA viruses in the families Genomoviridae and Smacoviridae as well as unclassified groups associated with the European hare (Lepus europaeus) and Snowshoe hare (Lepus americanus) [12,13].
To determine whether the Iberian hares found dead in Spain harbor viruses in addition to the MYXV-Tol strain, a metagenomic approach was used to identify circular DNA viruses. Using this approach, we identified circular DNA viruses for which we designed specific abutting primers to screen and recover full genomes from various Iberian hare samples collected in Spain following a myxomatosis outbreak. In these samples, virus genomes belonging to the Anelloviridae and Polyomaviridae families were identified and recovered. Neither of these have been previously reported in Leporids. The Anelloviridae family is a group of highly diverse circular single-stranded DNA viruses [14,15,16,17,18] with genome sizes ranging from ~2–4 kb. First identified in humans in 1997 [15], anelloviruses have been found in a wide range of animals, including pigs, non-human primates, seals, bats and horses [16,18,19,20,21], and are largely host species-specific. Polyomaviruses on the other hand are circular double-stranded DNA viruses of ~4–7 kb that, like anelloviruses, are host specific and have been identified in various animals [22,23,24]. In our analysis of various sample types from 11 hares, we identified 14 viruses belonging to two new species of anelloviruses and two viruses belonging to a new species of polyomavirus. Furthermore, we highlighted that three of the hares were coinfected with the MYXV-Tol variant, and either both or one of the novel anelloviruses and/or polyomaviruses.

2. Materials and Methods

2.1. Sample Collection and Viral Nucleic Acid Extraction

Samples from eleven symptomatic Iberian hares (19 samples in total from either ear, eyelid or vaginal sample) found dead in the fall of 2018 were collected and stored in RNAlater® (Sigma, St. Louis, MO, USA) at −80 °C, as shown in Table 1. These samples were collected from different regions of Spain: Toledo, Cuenca, Ciudad Real and Madrid. Viral nucleic acid was extracted using 5 g of each sample according to a modified phenol–chloroform extraction protocol [25] described in Águeda-Pinto et al. [8].

2.2. Detection of the MYXV-Tol

To detect the presence of the new MYXV-Tol, primers flanking the unique recombinant region (~2800 bp) were designed based on the previously described MYXV-Tol sequence (GenBank Accession MK836424), shown in Table S1, that yielded a ~4.8 kb amplicon. Viral nucleic acid from each individual sample was used as a template for PCR amplification using Kapa HiFi HotStart DNA polymerase (Kapa Biosystems, Wilmington, MA, USA) following manufactures’ recommendations at an annealing temperature of 60 °C. The amplicons were resolved on a 0.7% agarose gel stained with SYBR Safe (ThermoFisher Scientific, Waltham, MA, USA) and ~4.8 kb size fragments were excised, gel purified and cloned into pJET1.2 plasmid vector (ThermoFisher Scientific, Waltham, MA, USA). The resulting recombinant plasmids were Sanger sequenced by primer walking at Macrogen Inc. (Seoul, Korea). The sequence contigs were assembled using Geneious 11.1 [26].

2.3. Detection and Recovery of Anelloviruses and Polyomavirus Genomes

Previously, the MYXV-Tol variant from sample Lag01_EL, shown in Table 1, was cultured in permissive rabbit cells (RK13, ATCC # CCL-37), purified and sequenced using high-throughput sequencing, which led to the characterization of the full genome of a new myxoma strain, MYXV-Tol, GenBank Accession MK836424 [8].
In this study, viral DNA from the sample Lag01_EL was extracted using the Roche High Pure Viral nucleic acid kit (Roche Diagnostics, Indianapolis, IN, USA). Circular molecules were amplified using rolling circle amplification (RCA) with the Illustra TempliPhi 100 Amplification Kit (GE Healthcare, Chicago, IL, USA). The RCA amplicons were used to generate a 2 × 100 bp Illumina sequencing library and sequenced on an Illumina HiSeq4000 (Illumina, San Diego, CA, USA) at Macrogen Inc. (Seoul, Korea). The paired-end reads were de novo assembled using metaSPAdes 3.12.0 [27]. Contigs with terminal sequence redundancy were assumed to represent circular molecules and all contigs >750 nucleotides (nt) were analyzed using BLASTx [28] against a viral GenBank RefSeq protein database.
Based on these contigs, abutting primers were designed, shown in Table S1, and these primers were used to screen and recover full genomes from the 19 individual samples of 11 animals by PCR. The total DNA extracted from each of the 19 tissue samples was subjected to RCA to amplify circular molecules and 0.5 µL of this was used as a template with Kapa HiFi HotStart DNA polymerase (Kapa Biosystems, Wilmington, MA, USA) following manufactures’ recommendations at an annealing temperature of 60 °C with specific primers for the amplification of the anellovirus and polyomavirus genomes. The amplicons were resolved on a 0.7% agarose gel stained with SYBR Safe (ThermoFisher Scientific, Waltham, MA, USA). Amplicons (~2.5 kb for anelloviruses and ~5.4 kb for polyomaviruses) were excised, gel purified and cloned into pJET1.2 plasmid vector (ThermoFisher Scientific, Waltham, MA, USA). The resulting recombinant plasmids were Sanger sequenced by primer walking at Macrogen Inc. (Seoul, Korea). The sequence contigs were assembled using Geneious 11.1 [26].

2.4. Sequence Analysis of MYXV-Tol Regions, Anelloviruses and Polyomaviruses

2.4.1. MYXV-Tol Recombinant Region Sequences

A dataset of MYXV-Tol recombinant region sequences was assembled with sequences from this study (n = 7) and the two available in GenBank (GenBank accession # MK836424 and MK340973). These were aligned using MUSCLE [29] and any polymorphism, insertions or deletions were identified manually. The nucleotide sequence pairwise identities were determined using SDT [30].

2.4.2. Anellovirus Sequences Analyses

A dataset of the most closely related anelloviruses (n = 4) was created. The ORF1 sequences were extracted and an alignment of the ORF1 amino acid sequences was used to infer a maximum likelihood phylogenetic tree using PHYML [31] with the WAG+G substitution model, determined as the best fit model using ProtTest [32]. Branches with <0.8 aLRT support were collapsed using TreeGraph2 [33] and the resulting phylogenetic tree was midpoint rooted.

2.4.3. Polyomavirus Sequence Analyses

A dataset of representative polyomaviruses (n = 125) was downloaded from GenBank. From these, the large T-antigen and VP1 sequences were extracted. The aligned amino acid sequences of the large T-antigen was used to infer a maximum likelihood phylogenetic tree using PHYML [31] with the rtREV + G + I substitution model, determined as the best fit model using ProtTest [32]. Branches with <0.8 aLRT support were collapsed using TreeGraph2 [33] and the resulting phylogenetic tree was rooted with large T-antigen sequences of fish polyomaviruses [22]. The genome-wide nucleotide, the large T-antigen and the VP1 amino acid sequence pairwise identities were determined using SDT v1.2 [30].

3. Results and Discussion

After the outbreak of myxomatosis in Iberian hares, samples from Spain were collected to evaluate the presence of novel DNA viruses. Previously, the MYXV from one of the Iberian hare samples (Lag01; eyelid) was cultured in permissive rabbit cells (RK13, ATCC # CCL-37). Following purification, the genome of the MYXV was sequenced using high-throughput sequencing, which led to the identification of the recombinant region (~2.8 kb) derived from an unknown poxvirus in the MYXV genome and thus the recombinant strain, MYXV-Tol [8]. Following this report, other studies [9,10] also identified the MYXV containing the same novel ~2.8 kb recombinant region.
In this study, we identify the recombinant region of MYXV-Tol in seven of the 18 Iberian hare tissue samples from four hares (not including sample Lag01 eyelid from which the full MYXV was isolated). Sequence analyses of the recombinant region show that the obtained amplicons are >99% similar to the recombinant region in the MYXV-Tol (GenBank accession # MK836424) strain previously reported [8]. Among these sequences, three single nucleotide polymorphisms (SNPs) and two deletions were identified, shown in Figure 1. All the SNPs resulted in non-synonymous changes. One deletion of three nucleotides (ATC) located in the pox virus host range gene (GenBank accession # MT072320) resulted in the deletion of an amino acid (D) near the C terminus. In the sequence of the recombinant region of the MYXV-Tol from Lag04 V sample (GenBank accession # MT072323) there was a four nucleotide (TATA) deletion in the m009L gene resulting in a frame shift and extension of the open reading frame.
In the 19 Iberian hare tissue samples, a diverse range of unique anelloviruses (n = 14) were recovered, as seen in Table 1, suggesting that these viruses are common in Iberian hares. The hare- derived anelloviruses have a genome organization that has at least three ORFs and a conserved untranslated region (UTR). The hare anelloviruses range in size from 2490 to 2529 nt. Based on the last report of the International Committee on Taxonomy of Viruses (ICTV) [34], the classification of viruses in the family Anelloviridae is based on a global alignment pairwise identity of the ORF1 gene nucleotide sequence, with the cut-off values of >35% divergence for species and >56% divergence for genera. To facilitate the assignment of the anelloviruses at species and genera levels, a nucleotide pairwise comparison was generated, shown in Supplementary Data S1, and based on this, the anelloviruses genomes could be tentatively classified into two species which we refer to as Lepus torque teno virus (LepTTV) 1 and 2 from here on. Of the 14 genomes identified, four form a single LepTTV 1 supported clade, and 10 sequences form a second supported clade comprising of LepTTV 2s, shown in Figure 2. The closest relatives of the anelloviruses identified from hares are those identified in rodents (GenBank accession # MF684738 and MF684741) [35] and mosquitoes (GenBank accession # HQ335083 and HQ335082) [36]. The analysis of the pairwise identity for the ORF1 nucleotide, amino acid and full genome revealed that the LepTTVs share >63%, >56% and >70% pairwise identity with each other, respectively, shown in Figure 2 and Supplementary Data S1. When compared to anelloviruses previously identified in other organisms, LepTTVs share 61–67% (nt), 48–53% (aa), and 61–64% (full genome) pairwise identity. LepTTV 1 was identified in the eyelid (n = 2) and vaginal (n = 2) samples from three individual hares, shown in Table 1. LepTTV 2 was identified in the ear, eyelid (n = 5) and vaginal (n = 2) samples from three individual hares, shown in Table 1.
The polyomaviruses (5386 bp) were detected in two samples, the eyelid of Lag01 and Lag10, shown in Table 1. Typically, polyomavirus genomes are divided into three different regions: a non-coding region that regulates the expression of early and late genes; an early region that encodes the large T-antigen and the small T-antigen; and a late gene region that encodes the capsid proteins VP1, VP2, and VP3 [37,38]. In the genome of the polyomavirus recovered from the Iberian hare samples, we identified all the polyomavirus ORFs. The large T-antigen, VP1 and genome-wide pairwise analysis of the Lepus polyomavirus (LepPyV) and those most closely related revealed that these two LepPyV isolates share >99% identity and are most closely related to the Glis glis polyomavirus (GenBank accession # MG701352) [39] sharing 60% (large T-antigen), 67% (VP1) and 64% (full genome) pairwise identity, shown in Figure 3. The large T-antigen of the hare-derived polyomaviruses recovered in this study is nested within a larger mammal clade but forms a distinct sub-clade. The species demarcation for polyomaviruses is based on >15% divergence of the large T-antigen, thus the hare-derived polyomavirus represents a new species and we tentatively refer to this polyomavirus as Lepus polyomavirus 1 (LepPyV 1) and this would be part of the Betapolyomavirus genus [24,40].
Anelloviruses have a high degree of genetic variability between genotypes and strains and can establish persistent infections in their hosts, with no clear associated pathology [41]. However, it is still poorly understood how the host’s immune system is affected by the continuous exposure to anelloviruses [15,42]. Reports suggest that anelloviruses have a negative impact in several clinically relevant infections by modulating the baseline level of inflammation and interference with molecules of different signal pathways (reviewed in Maggi and Bendinelli [43]). For example, anellovirus ORF2 protein suppresses the NF-κB pathway, a key regulatory element that participates in the synthesis of proinflammatory cytokines like IL-6, IL-8, and cyclo-oxygenase-2 [44]. It was also found that anelloviruses encode miRNAs that could be involved in viral immune evasion and regulation of IFN signaling [45,46]. Studies have shown that the polyomavirus large T-antigen is capable of inactivating some host proteins such as p53, which is responsible for the control of cell-cycle inducing apoptosis [47,48]. Polyomaviruses are also able to prevent lysis of infected cells and initiate oncogenic transformation [38,49,50]. However, without further studies it is not possible to determine whether the coinfections of MYXV-Tol with LepTTV or LepPyV, present in Iberian hares, have any effect on the myxomatosis disease manifestation.

4. Concluding Remarks

In this study, we provide the first report of anelloviruses and polyomaviruses associated with hares. In addition, we identify the presence of the MYXV-Tol variant in four hares, and reveal coinfections in three hares, of either both or one of the novel anelloviruses and/or polyomaviruses. Hare populations have declined since the late 1960s which is thought to be in part due to disease [7,51,52]. The recent documentation of an outbreak of myxoma virus in hares caused by a prevalent recombinant MYXV-Tol [8,9,10] highlights another threat to the species. As an important food source for many predator species, such as the endangered Iberian lynx (Lynx pardinus), evaluating circulating pathogen populations and monitoring disease spread as well as spillover of pathogens is crucial for the survival of the hare and other species higher on the trophic cascade. Overall, this study expands our knowledge of DNA viruses associated with hares, highlighting the unsampled diversity of viruses, and the importance of investigating pathogen complexes within a host to further understand disease dynamics.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/12/3/340/s1. Table S1: Summary of primers used to amplify the MYXV-Tol recombinant cassette sequences, and the anellovirus and polyomavirus genomes. Supplementary Data S1: Pairwise comparison of the full genomes of anelloviruses e and the global alignment pairwise identity matrix of the ORF1 nucleotide used for species demarcation.

Author Contributions

Conceptualization, A.Á.-P., A.V., P.J.E.; Methodology, A.Á.-P., S.K., M.C.L., C.G., G.M., A.V., P.J.E.; Formal Analysis, A.Á.-P., S.K., M.C.L., A.V.; Resources, C.G., G.M., A.V., P.J.E.; Data Curation, A.Á.-P., S.K., A.V.; Writing—Original Draft Preparation, A.Á.-P., S.K., A.V.; Writing—Review and Editing, A.Á.-P., S.K., M.C.L., C.G., G.M., A.V., P.J.E.; Data visualization, A.Á.-P., S.K., A.V.; Supervision, G.M., A.V., P.J.E.; Project Administration, G.M., A.V., P.J.E.; Funding Acquisition, G.M., A.V., P.J.E. All authors have read and agreed to the published version of the manuscript.

Funding

All animal sampling took place post-mortem. According to EU and National legislation (2010/63/UE Directive, the Spanish Royal Decree (53/2013), and the University of Castilla–La Mancha guidelines, no permission or consent was required to conduct the research reported herein. FCT—Foundation for Science and Technology supported the doctoral fellowship of A.Á.-P. (ref. SFRH/BD/128752/2017) and the investigator grant of P.J.E. (IF/00376/2015). This article is a result of the project NORTE-01-0145-FEDER-000007, supported by Norte Portugal Regional Operational Programme (NORTE2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF). This work was also supported by the National Institute of Health (NIH) grant R01 A1080607.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Le Gall-Recule, G.; Lavazza, A.; Marchandeau, S.; Bertagnoli, S.; Zwingelstein, F.; Cavadini, P.; Martinelli, N.; Lombardi, G.; Guerin, J.L.; Lemaitre, E.; et al. Emergence of a new lagovirus related to Rabbit Haemorrhagic Disease Virus. Vet. Res. 2013, 44, 81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Le Gall-Recule, G.; Zwingelstein, F.; Fages, M.P.; Bertagnoli, S.; Gelfi, J.; Aubineau, J.; Roobrouck, A.; Botti, G.; Lavazza, A.; Marchandeau, S. Characterisation of a non-pathogenic and non-protective infectious rabbit lagovirus related to RHDV. Virology 2011, 410, 395–402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Strive, T.; Wright, J.D.; Robinson, A.J. Identification and partial characterisation of a new Lagovirus in Australian wild rabbits. Virology 2009, 384, 97–105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Stanford, M.M.; Werden, S.J.; McFadden, G. Myxoma virus in the European rabbit: Interactions between the virus and its susceptible host. Vet. Res. 2007, 38, 299–318. [Google Scholar] [CrossRef] [Green Version]
  5. Nash, P.; Barrett, J.; Cao, J.X.; Hota-Mitchell, S.; Lalani, A.S.; Everett, H.; Xu, X.M.; Robichaud, J.; Hnatiuk, S.; Ainslie, C.; et al. Immunomodulation by viruses: The myxoma virus story. Immunol. Rev. 1999, 168, 103–120. [Google Scholar] [CrossRef]
  6. Kerr, P.J.; Ghedin, E.; DePasse, J.V.; Fitch, A.; Cattadori, I.M.; Hudson, P.J.; Tscharke, D.C.; Read, A.F.; Holmes, E.C. Evolutionary history and attenuation of myxoma virus on two continents. PLoS Pathog. 2012, 8, e1002950. [Google Scholar] [CrossRef] [Green Version]
  7. Tsokana, C.N.; Sokos, C.; Giannakopoulos, A.; Birtsas, P.; Valiakos, G.; Spyrou, V.; Athanasiou, L.V.; Burriel, A.R.; Billinis, C. European Brown hare (Lepus europaeus) as a source of emerging and re-emerging pathogens of Public Health importance: A review. Vet. Med. Sci. 2020. [Google Scholar] [CrossRef]
  8. Agueda-Pinto, A.; Lemos de Matos, A.; Abrantes, M.; Kraberger, S.; Risalde, M.A.; Gortazar, C.; McFadden, G.; Varsani, A.; Esteves, P.J. Genetic Characterization of a Recombinant Myxoma Virus in the Iberian Hare (Lepus granatensis). Viruses 2019, 11, 530. [Google Scholar] [CrossRef] [Green Version]
  9. Dalton, K.P.; Martin, J.M.; Nicieza, I.; Podadera, A.; de Llano, D.; Casais, R.; Gimenez, S.; Badiola, I.; Aguero, M.; Duran, M.; et al. Myxoma virus jumps species to the Iberian hare. Transbound. Emerg. Dis. 2019, 66, 2218–2226. [Google Scholar] [CrossRef]
  10. Garcia-Bocanegra, I.; Camacho-Sillero, L.; Risalde, M.A.; Dalton, K.P.; Caballero-Gomez, J.; Aguero, M.; Zorrilla, I.; Gomez-Guillamon, F. First outbreak of myxomatosis in Iberian hares (Lepus granatensis). Transbound. Emerg. Dis. 2019, 66, 2204–2208. [Google Scholar] [CrossRef]
  11. Abade dos Santos, F.A.; Monteiro, M.; Pinto, A.; Carvalho, C.L.; Peleteiro, M.C.; Carvalho, P.; Mendonça, P.; Carvalho, T.; Duarte, M.D. First description of a herpesvirus infection in genus Lepus. bioRxiv 2020. [Google Scholar] [CrossRef]
  12. Steel, O.; Kraberger, S.; Sikorski, A.; Young, L.M.; Catchpole, R.J.; Stevens, A.J.; Ladley, J.J.; Coray, D.S.; Stainton, D.; Dayaram, A.; et al. Circular replication-associated protein encoding DNA viruses identified in the faecal matter of various animals in New Zealand. Infect. Genet. Evol. 2016, 43, 151–164. [Google Scholar] [CrossRef] [PubMed]
  13. Kraberger, S.; Waits, K.; Ivan, J.; Newkirk, E.; VandeWoude, S.; Varsani, A. Identification of circular single-stranded DNA viruses in faecal samples of Canada lynx (Lynx canadensis), moose (Alces alces) and snowshoe hare (Lepus americanus) inhabiting the Colorado San Juan Mountains. Infect. Genet. Evol. 2018, 64, 1–8. [Google Scholar] [CrossRef]
  14. Nishiyama, S.; Dutia, B.M.; Stewart, J.P.; Meredith, A.L.; Shaw, D.J.; Simmonds, P.; Sharp, C.P. Identification of novel anelloviruses with broad diversity in UK rodents. J. Gen. Virol. 2014, 95, 1544–1553. [Google Scholar] [CrossRef] [PubMed]
  15. Spandole, S.; Cimponeriu, D.; Berca, L.M.; Mihaescu, G. Human anelloviruses: An update of molecular, epidemiological and clinical aspects. Arch. Virol. 2015, 160, 893–908. [Google Scholar] [CrossRef] [PubMed]
  16. Hrazdilova, K.; Slaninkova, E.; Brozova, K.; Modry, D.; Vodicka, R.; Celer, V. New species of Torque Teno miniviruses infecting gorillas and chimpanzees. Virology 2016, 487, 207–214. [Google Scholar] [CrossRef] [Green Version]
  17. Fahsbender, E.; Burns, J.M.; Kim, S.; Kraberger, S.; Frankfurter, G.; Eilers, A.A.; Shero, M.R.; Beltran, R.; Kirkham, A.; McCorkell, R.; et al. Diverse and highly recombinant anelloviruses associated with Weddell seals in Antarctica. Virus Evol. 2017, 3. [Google Scholar] [CrossRef] [Green Version]
  18. Crane, A.; Goebel, M.E.; Kraberger, S.; Stone, A.C.; Varsani, A. Novel anelloviruses identified in buccal swabs of Antarctic fur seals. Virus Genes 2018, 54, 719–723. [Google Scholar] [CrossRef]
  19. Bigarre, L.; Beven, V.; de Boisseson, C.; Grasland, B.; Rose, N.; Biagini, P.; Jestin, A. Pig anelloviruses are highly prevalent in swine herds in France. J. Gen. Virol. 2005, 86, 631–635. [Google Scholar] [CrossRef]
  20. Cibulski, S.P.; Teixeira, T.F.; de Sales Lima, F.E.; do Santos, H.F.; Franco, A.C.; Roehe, P.M. A Novel Anelloviridae Species Detected in Tadarida brasiliensis Bats: First Sequence of a Chiropteran Anellovirus. Genome Announc. 2014, 2. [Google Scholar] [CrossRef] [Green Version]
  21. Li, L.; Giannitti, F.; Low, J.; Keyes, C.; Ullmann, L.S.; Deng, X.; Aleman, M.; Pesavento, P.A.; Pusterla, N.; Delwart, E. Exploring the virome of diseased horses. J. Gen. Virol. 2015, 96, 2721–2733. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Van Doorslaer, K.; Kraberger, S.; Austin, C.; Farkas, K.; Bergeman, M.; Paunil, E.; Davison, W.; Varsani, A. Fish polyomaviruses belong to two distinct evolutionary lineages. J. Gen. Virol. 2018, 99, 567–573. [Google Scholar] [CrossRef] [PubMed]
  23. Buck, C.B.; Van Doorslaer, K.; Peretti, A.; Geoghegan, E.M.; Tisza, M.J.; An, P.; Katz, J.P.; Pipas, J.M.; McBride, A.A.; Camus, A.C.; et al. The Ancient Evolutionary History of Polyomaviruses. PLoS Pathog. 2016, 12, e1005574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Calvignac-Spencer, S.; Feltkamp, M.C.; Daugherty, M.D.; Moens, U.; Ramqvist, T.; Johne, R.; Ehlers, B.; Polyomaviridae Study Group of the International Committee on Taxonomy of Viruses. A taxonomy update for the family Polyomaviridae. Arch. Virol. 2016, 161, 1739–1750. [Google Scholar] [CrossRef] [Green Version]
  25. Di Pietro, F.; Ortenzi, F.; Tilio, M.; Concetti, F.; Napolioni, V. Genomic DNA extraction from whole blood stored from 15- to 30-years at -20 degrees C by rapid phenol-chloroform protocol: A useful tool for genetic epidemiology studies. Mol. Cell. Probes 2011, 25, 44–48. [Google Scholar] [CrossRef]
  26. 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]
  27. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [Green Version]
  28. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  29. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  30. Muhire, B.M.; Varsani, A.; Martin, D.P. SDT: A virus classification tool based on pairwise sequence alignment and identity calculation. PLoS ONE 2014, 9, e108277. [Google Scholar] [CrossRef]
  31. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. ProtTest 3: Fast selection of best-fit models of protein evolution. Bioinformatics 2011, 27, 1164–1165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Stover, B.C.; Muller, K.F. TreeGraph 2: Combining and visualizing evidence from different phylogenetic analyses. BMC Bioinform. 2010, 11, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Biagini, P. Classification of TTV and related viruses (anelloviruses). Curr. Top. Microbiol. Immunol. 2009, 331, 21–33. [Google Scholar] [CrossRef] [PubMed]
  35. Du, J.; Li, Y.; Lu, L.; Zheng, D.; Liu, B.; Yang, L.; Su, H.; Dong, J.; Sun, L.; Zhu, Y.; et al. Biodiversity of rodent anelloviruses in China. Emerg. Microbes Infect. 2018, 7, 38. [Google Scholar] [CrossRef] [Green Version]
  36. Ng, T.F.; Willner, D.L.; Lim, Y.W.; Schmieder, R.; Chau, B.; Nilsson, C.; Anthony, S.; Ruan, Y.; Rohwer, F.; Breitbart, M. Broad surveys of DNA viral diversity obtained through viral metagenomics of mosquitoes. PLoS ONE 2011, 6, e20579. [Google Scholar] [CrossRef] [Green Version]
  37. Hirsch, H.H.; Steiger, J. Polyomavirus BK. Lancet Infect. Dis. 2003, 3, 611–623. [Google Scholar] [CrossRef]
  38. Dalianis, T.; Hirsch, H.H. Human polyomaviruses in disease and cancer. Virology 2013, 437, 63–72. [Google Scholar] [CrossRef] [Green Version]
  39. Ehlers, B.; Anoh, A.E.; Ben Salem, N.; Broll, S.; Couacy-Hymann, E.; Fischer, D.; Gedvilaite, A.; Ingenhutt, N.; Liebmann, S.; Martin, M.; et al. Novel Polyomaviruses in Mammals from Multiple Orders and Reassessment of Polyomavirus Evolution and Taxonomy. Viruses 2019, 11, 930. [Google Scholar] [CrossRef] [Green Version]
  40. Moens, U.; Calvignac-Spencer, S.; Lauber, C.; Ramqvist, T.; Feltkamp, M.C.W.; Daugherty, M.D.; Verschoor, E.J.; Ehlers, B.; Ictv Report Consortium. ICTV Virus Taxonomy Profile: Polyomaviridae. J. Gen. Virol. 2017, 98, 1159–1160. [Google Scholar] [CrossRef]
  41. Freer, G.; Maggi, F.; Pifferi, M.; Di Cicco, M.E.; Peroni, D.G.; Pistello, M. The Virome and Its Major Component, Anellovirus, a Convoluted System Molding Human Immune Defenses and Possibly Affecting the Development of Asthma and Respiratory Diseases in Childhood. Front. Microbiol. 2018, 9, 686. [Google Scholar] [CrossRef] [Green Version]
  42. McElvania TeKippe, E.; Wylie, K.M.; Deych, E.; Sodergren, E.; Weinstock, G.; Storch, G.A. Increased prevalence of anellovirus in pediatric patients with fever. PLoS ONE 2012, 7, e50937. [Google Scholar] [CrossRef] [PubMed]
  43. Maggi, F.; Bendinelli, M. Immunobiology of the Torque Teno Viruses and Other Anelloviruses; TT Viruses; Springer: Berlin/Heidelberg, Germany, 2009; pp. 65–90. [Google Scholar]
  44. Zheng, H.; Ye, L.; Fang, X.; Li, B.; Wang, Y.; Xiang, X.; Kong, L.; Wang, W.; Zeng, Y.; Ye, L.; et al. Torque teno virus (SANBAN isolate) ORF2 protein suppresses NF-kappaB pathways via interaction with IkappaB kinases. J. Virol. 2007, 81, 11917–11924. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kincaid, R.P.; Burke, J.M.; Cox, J.C.; de Villiers, E.M.; Sullivan, C.S. A human torque teno virus encodes a microRNA that inhibits interferon signaling. PLoS Pathog. 2013, 9, e1003818. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Vignolini, T.; Macera, L.; Antonelli, G.; Pistello, M.; Maggi, F.; Giannecchini, S. Investigation on torquetenovirus (TTV) microRNA transcriptome in vivo. Virus Res. 2016, 217, 18–22. [Google Scholar] [CrossRef]
  47. Tevethia, M.J.; Pipas, J.M.; Kierstead, T.; Cole, C. Requirements for immortalization of primary mouse embryo fibroblasts probed with mutants bearing deletions in the 3’ end of SV40 gene A. Virology 1988, 162, 76–89. [Google Scholar] [CrossRef]
  48. Chao, H.H.; Buchmann, A.M.; DeCaprio, J.A. Loss of p19 (ARF) eliminates the requirement for the pRB-binding motif in simian virus 40 large T antigen-mediated transformation. Mol. Cell. Biol. 2000, 20, 7624–7633. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Cheng, J.; DeCaprio, J.A.; Fluck, M.M.; Schaffhausen, B.S. Cellular Transformation by Simian Virus 40 and Murine Polyoma Virus T Antigens. Semin. Cancer Biol. 2009, 19, 218–228. [Google Scholar] [CrossRef] [Green Version]
  50. Neu, U.; Stehle, T.; Atwood, W.J. The Polyomaviridae: Contributions of virus structure to our understanding of virus receptors and infectious entry. Virology 2009, 384, 389–399. [Google Scholar] [CrossRef] [Green Version]
  51. Ferreras, P.; Travaini, A.; Zapata, S.C.; Delibes, M. Short-term responses of mammalian carnivores to a sudden collapse of rabbits in Mediterranean Spain. Basic Appl. Ecol. 2011, 12, 116–124. [Google Scholar] [CrossRef] [Green Version]
  52. Smith, R.K.; Jennings, N.V.; Harris, S. A quantitative analysis of the abundance and demography of European hares Lepus europaeus in relation to habitat type, intensity of agriculture and climate. Mammal Rev. 2005, 35. [Google Scholar] [CrossRef]
Viruses 12 00340 g001 550
Figure 1. Myxoma virus genome schematic with recombinant region highlighted. The variations identified in the recombinant gene cassette of the sequenced MYXV-Tol variants from various samples are summarized in reference to the MYXV-Tol GenBank accession. SNPs are shown with red highlighted nucleotide and the encoded amino acid in brackets. Deletions are shown with Δ symbol and the deleted amino acid in brackets.
Figure 1. Myxoma virus genome schematic with recombinant region highlighted. The variations identified in the recombinant gene cassette of the sequenced MYXV-Tol variants from various samples are summarized in reference to the MYXV-Tol GenBank accession. SNPs are shown with red highlighted nucleotide and the encoded amino acid in brackets. Deletions are shown with Δ symbol and the deleted amino acid in brackets.
Viruses 12 00340 g001
Viruses 12 00340 g002 550
Figure 2. Maximum likelihood phylogenetic tree of the anelloviruses ORF1 protein sequences of those recovered from hares together with those most closely related, available in GenBank. A pairwise comparison color matrix of the ORF1 amino acid and nucleotide sequence is shown to the right of the phylogenetic tree.
Figure 2. Maximum likelihood phylogenetic tree of the anelloviruses ORF1 protein sequences of those recovered from hares together with those most closely related, available in GenBank. A pairwise comparison color matrix of the ORF1 amino acid and nucleotide sequence is shown to the right of the phylogenetic tree.
Viruses 12 00340 g002
Viruses 12 00340 g003 550
Figure 3. Maximum likelihood phylogenetic tree of the polyomavirus large T-antigen amino acid sequences of polyomaviruses identified in hares and representative sequences of polyomavirus species available in GenBank. The enlarged clade shows those most closely related to the hare polyomavirus. A pairwise comparison color matrix of the large-T protein, VP1 protein and the full genome is provided to the right of the phylogenetic tree.
Figure 3. Maximum likelihood phylogenetic tree of the polyomavirus large T-antigen amino acid sequences of polyomaviruses identified in hares and representative sequences of polyomavirus species available in GenBank. The enlarged clade shows those most closely related to the hare polyomavirus. A pairwise comparison color matrix of the large-T protein, VP1 protein and the full genome is provided to the right of the phylogenetic tree.
Viruses 12 00340 g003
Table
Table 1. Sample information of the hare samples analyzed in this study including geographic origin, tissue type and GenBank accession numbers of the viral sequences.
Table 1. Sample information of the hare samples analyzed in this study including geographic origin, tissue type and GenBank accession numbers of the viral sequences.
PoxvirusAnellovirusPolyomavirus
Animal IDLocationTissue TypeMYXV-TolLepTTV1LepTTV2LepPyV1
Lag01Spain: ToledoEyelidMK836424MN994854MN994855MN994868
Lag02Spain: CuencaEar----
Lag03Spain: Ciudad RealEar----
Lag04Spain: Ciudad RealVaginalMT072323---
Lag05Spain: Ciudad RealEar----
Eyelid ---
Lag06Spain: Ciudad RealEarMT072322---
EyelidMT072321-MN994856-
VaginalMT072320---
Lag07Spain: Ciudad RealEarMT072319---
EyelidMT072318---
Lag10SpainEyelid-MN994857MN994858MN994869
Vaginal-MN994859MN994860
MN994861		-
Lag11Spain: ToledoEyelid--MN994862-
Vaginal----
Lag12Spain: Ciudad RealEar----
Eyelid----
Lag14Spain: MadridEyelid--MN994863
MN994864		-
VaginalMT072317MN994865MN994866
MN994867		-

Share and Cite

MDPI and ACS Style

Águeda-Pinto, A.; Kraberger, S.; Lund, M.C.; Gortázar, C.; McFadden, G.; Varsani, A.; Esteves, P.J. Coinfections of Novel Polyomavirus, Anelloviruses and a Recombinant Strain of Myxoma Virus-MYXV-Tol Identified in Iberian Hares. Viruses 2020, 12, 340. https://0-doi-org.brum.beds.ac.uk/10.3390/v12030340

AMA Style

Águeda-Pinto A, Kraberger S, Lund MC, Gortázar C, McFadden G, Varsani A, Esteves PJ. Coinfections of Novel Polyomavirus, Anelloviruses and a Recombinant Strain of Myxoma Virus-MYXV-Tol Identified in Iberian Hares. Viruses. 2020; 12(3):340. https://0-doi-org.brum.beds.ac.uk/10.3390/v12030340

Chicago/Turabian Style

Águeda-Pinto, Ana, Simona Kraberger, Michael C. Lund, Christian Gortázar, Grant McFadden, Arvind Varsani, and Pedro J. Esteves. 2020. "Coinfections of Novel Polyomavirus, Anelloviruses and a Recombinant Strain of Myxoma Virus-MYXV-Tol Identified in Iberian Hares" Viruses 12, no. 3: 340. https://0-doi-org.brum.beds.ac.uk/10.3390/v12030340

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