Detection and Genetic Characterization of Viruses Present in Free-Ranging Snow Leopards Using Next-Generation Sequencing

Snow leopards inhabit the cold, arid environments of the high mountains of South and Central Asia. These living conditions likely affect the abundance and composition of microbes with the capacity to infect these animals. It is important to investigate the microbes that snow leopards are exposed to detect infectious disease threats and define a baseline for future changes that may impact the health of this endangered felid. In this work, next-generation sequencing is used to investigate the fecal (and in a few cases serum) virome of seven snow leopards from the Tost Mountains of Mongolia. The viral species to which the greatest number of sequences reads showed high similarity was rotavirus. Excluding one animal with overall very few sequence reads, four of six animals (67%) displayed evidence of rotavirus infection. A serum sample of a male and a rectal swab of a female snow leopard produced sequence reads identical or closely similar to felid herpesvirus 1, providing the first evidence that this virus infects snow leopards. In addition, the rectal swab from the same female also displayed sequence reads most similar to feline papillomavirus 2, which is the first evidence for this virus infecting snow leopards. The rectal swabs from all animals also showed evidence for the presence of small circular DNA viruses, predominantly Circular Rep-Encoding Single-Stranded (CRESS) DNA viruses and in one case feline anellovirus. Several of the viruses implicated in the present study could affect the health of snow leopards. In animals which are under environmental stress, for example, young dispersing individuals and lactating females, health issues may be exacerbated by latent virus infections.

Introduction

Infectious diseases can affect the abundance and distribution of animals by reducing survival and reproduction (1, 2). Even small changes in these parameters can substantially increase the extinction risk, for example, in species with slow reproduction or where populations are small (3). This applies to many large carnivore populations as they exhibit a low reproductive output, occur at low densities naturally, and populations are often further reduced and isolated as a consequence of habitat destruction, overexploitation of prey species, and human persecution (4). Reduced survival can also disrupt the social system, where the replacement of dominant males can result in infanticide of young, thus further reducing population numbers (4–6). The detrimental effects that a disease outbreak can cause a felid population are well-exemplified by the canine distemper virus outbreak in Serengeti lions (Panthera leo) in 1994. The outbreak, probably originating from domestic dogs (Canis familiaris), claimed about 30% of the lion population and is one of the most cited examples of the potential impact that disease can have on felid populations (7). In contrast to the Serengeti lions, which are frequently observed, such an outbreak could occur unnoticed in snow leopards and other less studied felids.

The snow leopard (Panthera uncia) is a large felid inhabiting the high mountains of South and Central Asia. The species distribution range covers 1.2–1.6 million km², spanning over 12 countries (8). Snow leopards appear to utilize relatively large territories (9, 10) and occur at low densities (0.9–1.8 adults/100 km²) (11). The cold, arid environment inhabited by snow leopards likely has lower microbial abundance than in more temperate and mesic habitats (12). Consequently, snow leopards should encounter disease agents less frequently than many other carnivores and may therefore exhibit lower intrinsic levels of immunity, rendering it vulnerable to disease outbreaks (13, 14). Snow leopards are likely susceptible to most infectious diseases that are known to affect the domestic cat (Felis catus), in addition, spill over from other felids and prey animals contribute to the spectrum of infectious diseases that could affect the health of snow leopards (13, 14). Direct routes of transmission are likely most common, these include both intraspecies contact (mating, fighting, socializing) and interspecies contact with wild and domestic prey, other carnivores, and scavengers. To a degree, indirect routes of transmission such as water holes, carcasses, marking sites, and human settlements could also play a role in disease transmission.

Except for a few incidental reports and a recent study focusing on tick-borne bacteria from ticks collected from snow leopards and protozoan infections evidenced by serology (15), there are no published data on infectious diseases in free-ranging snow leopards, owing to their remote and inaccessible habitat, combined with the species' secretive nature. Accordingly, there is very limited information regarding the prevalence and thus potential threat of infectious diseases to snow leopards or which microorganisms are most commonly found in the species (14). Therefore, to begin mapping the occurrence of infectious agents in free-ranging snow leopards, we collected samples from eight animals captured in a telemetry study (16). In the present study, we report viral sequences that were obtained by next-generation sequencing (NGS) of rectal swab samples and, for a few animals, also serum samples from snow leopards in the Tost Mountains of Mongolia.

Materials and Methods

Study Area

This study was conducted in the Tost Mountains (43°N, 100°E), a relatively isolated range of mountain massifs (1,600–2,500 m.a.s.l.) in the Gobi Desert in southern Mongolia. Temperatures range from 38°C in the summer to −35°C in the winter, and the annual precipitation is <130 mm of which most falls as rain from June to August. An estimated 10–14 adult snow leopards inhabit the 1,700 km² large area (17) where they prey mainly on ibex (Capra sibirica), domestic goats (Capra aegagrus hircus), and argali sheep (Ovis ammon) (18). The human population consists of ~90 semi-nomadic herder families who move seasonally with their livestock, comprising of ~32,000 goats and sheep (Ovis aries), ~1,100 camels (Camelus bactrianus), and ~120 horses (Equus ferus caballus). In addition, most families have at least one dog that either follow the livestock or roam freely, usually spending the night close to the family's camp. Sympatric predators include gray wolf (Canis lupus), Eurasian lynx (Lynx lynx), red fox (Vulpes vulpes), and marten (Martes spp.).

Study Animals and Sampling

Snow leopards were captured using modified Aldrich-style foot snares where cats were immobilized with a combination of medetomidine and tiletamine–zolazepam. See reference (16) for a detailed description of capture procedures. We collected rectal swabs from eight snow leopards (three adult males, one subadult male, and four females) from October 2011 to May 2013 by inserting an Amies charcoal cotton swab (Copan Italia S.p.A.) into the rectum and moving it along the wall. The swabs were immediately placed in the charcoal medium. Blood was collected from the cephalic vein and placed in 4 ml blood serum separating tubes (BD vacutainer, Plymouth, UK). Serum tubes stood for 12–15 h to separate cells and serum, which was then decanted into cryovials. We obtained enough serum for NGS from three individuals. In the camp, swabs were stored in a cool and dark box, and sera were kept at −18°C until the samples could be transported to the National Veterinary Institute in Uppsala, Sweden (within 0.5–2.5 months after capture).

Sample Preparation

Rectal swabs from the eight animals (F3, F7, F8, F9, M1, M7, M9, M10) were immersed in 1,200 μl TE-buffer [100 mM Tris, 10 mM ethylenediaminetetraacetic acid (EDTA), pH 8.0] and homogenized by shaking the tubes vigorously for 1 h at room temperature. After centrifugation (200 g for 5 min), the supernatant was collected and filtered through a 0.45 μm filter to remove particles of bacterium-size and larger. The filtrate was treated with 400 U/ml of DNase I (Roche Applied Science) and 10 μg/ml of RNase A (Invitrogen) in 1× DNase buffer (Roche Applied Science) at 37°C for 2 h to degrade unprotected nucleic acids. Nuclease treatment was also applied to serum samples for three animals (F3, M1, M9) in the same way, and subsequently, DNA and RNA were extracted from 200 μl each of the nuclease-treated samples (the M1 serum sample was also prepared without nuclease treatment) using QIAamp DNA Mini Kit (Qiagen) and a combination of TRIzol (Invitrogen) and RNeasy kit (Qiagen), respectively, following the manufactures' instructions. Sequence-independent single-primer amplification (SISPA) was applied to both RNA and DNA as follows: SuperScript III first-strand synthesis kit (Invitrogen) was used to generate cDNA from the RNA preparation. The reaction was primed by a primer FR20RV-6N (19) following the manufacturer's instructions. Double-stranded DNA was obtained by incubation of the cDNA products with Klenow Fragment DNA polymerase (New England Biolabs) at 37°C for 1 h. Subsequently, the Klenow enzyme was inactivated at 75°C for 10 min. DNA templates were also tagged with the FR20RV-6N primer during a Klenow Fragment reaction at the same conditions. Random amplification of the tagged cDNA and DNA was performed using the primer FR20RV (19) under the following conditions: 10 min at 95°C, followed by 40 cycles of 30 s at 95°C, 30 s at 58°C, and 90 s at 72°C. The reaction was ended with an extra elongation step at 72°C for 10 min. The PCR reaction contained 1× PCR buffer, 2.5 mM MgCl₂, 2.5 mM dNTPs, 0.4 mM primer, and 1.25 U AmpliTaq Gold DNA polymerase (Applied Biosystems). The amplified DNA fragments were further treated with EcoRV (New England Biolabs) to remove the amplification primers and purified by QIAquick PCR purification kit (Qiagen). Concentration was measured with a Qubit fluorometer using Qubit dsDNA HS (High Sensitivity) Assay Kit (Invitrogen), and a 0.2 ng/μl aliquot was prepared for each sample. Nextera XT DNA Library Preparation Kit (Illumina, Inc.) was used to fragment the input DNA and tag the DNA from each sample with a pair of unique index primers by a 12-cycle PCR amplification. The libraries were purified with AMPure XP beads (Sigma), and Agilent High-Sensitivity DNA Kit (Agilent) was used to verify the length distribution of the fragments and for quantification of the libraries. Finally, an equimolar amount (2 nM) of each sample library with sufficient quality and concentration was pooled, thus constituting two pools of seven serum sample preparations and 15 rectal sample preparations, respectively. The pools were denatured and further diluted to a final concentration of 10 pM. Sequencing was performed on a MiSeq desktop sequencer using MiSeq 500 cycles reagent kit 500 (v. 2) (Illumina, Inc.). Library preparation and sequencing were performed according to the manufacturer's instructions.

Bioinformatics

The sequence reads were assigned to species by homology searching the NCBI nt database with the BLASTn algorithm as implemented on a Decypher server (TimeLogic^®, Carlsbad, CA). Before blasting, the sequence reads were quality checked and trimmed using HTStream v.1.0.0 (HTStream, RRID:SCR_018354) (20). To reduce the computational burden, the trimmed reads were first blasted against the VRL section of the NCBI nt database (i.e., the viral sequences) with a cutoff expect value (e-value) of 10⁻⁵. The reads with BLASTn hits to the VRL database were collected with an in-house python script and subsequently blasted (BLASTn) against the whole nt database with the same e-value. Finally, the reads with best hits to viral sequences in the nt database were collected with an in-house python script. This procedure reduced the computational burden to extract the reads with closest homology to viral sequences in the nt NCBI database with about 90%. These selected subsets of sequence reads were then subjected to de novo assembly using the CLC genomics workbench. The reads were mapped back to the assembled contigs, and consensus sequences were extracted by using a low coverage threshold of 5 and a noise threshold of 0.3 (with a secondary nucleotide present above this fraction, a degenerate nucleotide was inserted in the consensus sequence).

The consensus sequences were aligned with related reference sequences for GenBank using a gap open and extension penalty of 10 and 1, respectively. End gaps were free since sequences of different lengths were aligned. Before phylogenetic trees were constructed, the sequences of the alignment were trimmed to the same length. The phylogenetic trees were created using the neighborhood-joining method with the Jukes–Cantor distance measure. The topology was verified using 1,000 bootstrap calculations, and only branches with better than 60% support are shown.

Genotyping of rotavirus A (RVA) were carried out by the reimplementation of the RotaC^2.0 (21) at the NIAID Virus Pathogen Database and Analysis Resource [Virus Pathogen Resource (ViPR), RRID:SCR_012983] (22, 23).

Results

Extracted RNA and DNA from eight rectal swabs and three sera were processed for NGS. Only seven rectal samples and three serum samples produced enough sequence reads for further analysis. The samples produced around 2 million sequence reads each. The inference of virus sequences in the snow leopard rectal swabs samples, and sera, was carried out by finding the most closely homologous sequence in GenBank by the blastn algorithm (24). The results are summarized in Table 1. As expected, most of the reads were bacterial despite the experimental procedures to enrich viral sequences in the extracted nucleic acids. The number of viral reads varied in the samples in the range 38–1,844. Generally, mostly DNA viruses were detected with the notable exception of RVA, a double-stranded RNA virus, which was clearly the most abundant virus species.

TABLE 1

Table 1. Summary of blastn hits to the GenBank nt (database release 228) supported by at least four reads each.

Rotavirus A

The most striking feature of the rectal swab samples is the large abundance of RVA reads. Four of seven samples (F1, M1, M10, and F7) contain RVA with the number of sequence reads in the range 88–1,635. From animal M1 and M10, a 1,036- and 1,161-nt-long fragment of the NSP1 and VP3 gene, respectively, could be assembled, which allowed genotyping of these fragments. The average read coverage of these contigs was 76 and 223 reads, respectively. After mapping back the reads to the contigs and applying a conservative coverage lower limit of 5 reads, the lengths of the contigs were reduced to 771 and 1,127 nt. These trimmed contigs were used for phylogenetic analysis. The NSP1 contig was of genotype A3 and the VP3 contig was of genotype M2, as determined by the RVA genotyping tool available at the ViPR resource (https://www.viprbrc.org). Of the 11 NSP1 genes from feline RVA available at the ViPR database, the most common genotype was A3 with five members. There were also 11 VP3 genes from feline RVA strains represented at the ViPR resource. Of these, two were of the M2 genotype and the remaining nine were all of the M3 genotype. Phylogenetic analysis of the VP3 fragment shows that the snow leopard fragment is most similar to a human isolate from Japan and is found in a clade with rotaviruses isolated from many different host species, including isolates from primates, bovids, cervids, cats, pigs, and horse (Figure 1). The NSP1 sequence is also found in clade containing rotavirus isolates from humans and bovids (Figure 2). It should be noted that the feline isolates available in GenBank are all isolated from domestic cats, and we have not been able to find rotavirus sequences for comparison from wild felids. In general, the rotavirus sequences found in snow leopard are not related to rotavirus sequences isolated from cats available in GenBank (Figures 1, 2).

FIGURE 1

Figure 1. Neighbor-joining phylogenetic tree for a 1,127-nt-long fragment of the rotavirus VP3 gene obtained from the rectal swab of animal M10 together with selected sequences from GenBank (NCBI accessions are given within parentheses). These were selected to include the representative from different hosts with the highest similarity, except that all available isolates from felids were included if not identical. The phylogenetic trees were constructed with the CLC genomics workbench using the neighbor-joining method with the Jukes–Cantor distance measure with 1,000 bootstrap calculations. Only branches with at least 60% bootstrap support are shown. The NCBI accession of the novel snow leopard VP3 sequence is shown in red.

FIGURE 2

Figure 2. Neighbor-joining phylogenetic tree for a 771-nt-long fragment of the rotavirus NSP1 gene obtained from the rectal swab of animal M1 together with selected sequences from GenBank (NCBI accessions are given within parentheses). These were selected to include the representative from different hosts with the highest similarity, except that all available isolates from felids were included if not identical. The phylogenetic trees were constructed with the CLC genomics workbench using the neighbor-joining method with the Jukes–Cantor distance measure with 1,000 bootstrap calculations. Only branches with at least 60% bootstrap support are shown. The NCBI accession of the novel snow leopard NSP1 sequence is shown in red.

Herpesviruses

Most samples from the seven animals had sequence reads that were similar to various herpesviruses (Table 1). However, few reads with relatively high e-values and blastn hits to herpesviruses found in distantly related host species make most of these hits of questionable relevance. The exceptions are the rectal swab from F9 (Table 1) and the serum sample from M1 (not in Table 1 since only two reads were found), which after assembly of the reads contained a 228- and 289-nt contig, respectively, with sequence identical or closely similar to segments of felid herpesvirus 1 (e.g., GenBank accession MH070348).

Papillomaviruses

The swab from one of the female snow leopards (F9) produced 108 reads, which could be assembled into a 489-nt-long contig from the E1 gene most similar to feline papillomavirus type 2 (Table 1). After mapping back the reads on the contigs (average coverage 43 reads) and applying a conservative lower limit for coverage of 5 reads, the length of the contig was reduced to 435 nt. Phylogenetic analysis of this E1 fragment with related papillomaviruses shows that the papillomavirus sequence found in this snow leopard female was most similar to the feline papillomavirus type 2 (Figure 3). The nucleotide identity with feline papillomavirus type 2 (NC_038520) is 74%, and the amino acid identity is 91%.

FIGURE 3

Figure 3. Neighbor-joining phylogenetic tree for a 435-nt-long fragment of the papillomavirus E1 gene obtained from the rectal swab of animal F9 together with selected sequences from GenBank (NCBI accessions are given within parentheses). These were selected to include viral sequences obtained from related host species or for showing high similarity. The phylogenetic trees were constructed with the CLC genomics workbench using the neighbor-joining method with the Jukes–Cantor distance measure with 1,000 bootstrap calculations. Only branches with at least 60% bootstrap support are shown. The NCBI accession of the snow leopard papillomavirus sequence is shown in red.

Circular Rep-Encoding Single-Stranded DNA Viruses and Anelloviruses

The rectal swabs of four of seven animals (M1, M7, F7, M10) contained viral sequences that were very similar to virus Chimp162 (Table 1). This is a Rep-protein sequence obtained from chimpanzee feces sample collected in Uganda and likely originates from a Circular Rep-Encoding Single-Stranded (CRESS) DNA virus (25). The prevalence of CRESS DNA viruses is further manifested in animal F7 by blast hits to Dragonfly larvae-associated circular virus-3 (26), Lake Sarah-associated circular virus-36 (27), and related viruses in the serum of animal F3 and several hits to the Genomoviridae family for animal M10 (Table 1). All these blast hits indicate that CRESS DNA viruses are common in snow leopards and can be found both in the digestive tract and in the blood. Indeed, circoviruses also belong to the CRESS DNA virus group, and the fecal swabs of the snow leopard males M1 and M10 both show evidence of the presence of a circovirus. In the case of M1, sequence reads similar to bat circovirus are found both in the swab sample and in the sera. The fecal swab from the male snow leopard M9 did not show any reads related to CRESS DNA viruses; however, the nucleic acid preparation optimized for DNA viruses showed 35 reads with the highest similarity to the Anelloviridae family. The feline anellovirus displayed the closest similarity with an e-value of 2 × 10⁻¹⁷ (Table 1), which is relatively high, indicating that the putative snow leopard anellovirus only is distantly related to the feline anellovirus.

Other Viruses

Samples from animal F7 and M1 have sequence reads with homology to porcine stool-associated virus 1 (posavirus 1) (28). The sequence similarity is low as measured by the e-values (Table 1), and the number of reads is low. There is thus some evidence that a virus related to posavirus-1 of the Picornavirales order may be present in the snow leopards, but this requires further substantiation. This is also true for the reads found with weak homology to astroviruses (animals M1 and M10; Table 1). Astroviruses are frequently found in stool samples from many mammals, and our study suggests that astrovirus may also infect snow leopards, but further investigations are needed to verify astrovirus infections in snow leopards. Six reads are most similar to a volepox virus in animal F9 (Table 1), but the similarity is relatively low, and the finding needs further support. Poxviruses do not normally infect felids, although cowpox is an exception (29). Snow leopards feed on a large array of prey species and can occasionally kill smaller animals such as voles (30), and the volepox virus reads might indicate infected prey animal rather than infection of the snow leopard.

Discussion

NGS provides a powerful tool for metagenomic investigation and has led to identification of a multitude of novel viruses (31–33). In the present work, we have investigated the intestinal virome of snow leopards by rectal swab sampling and also sera for a subset of the animals. Respiratory viruses are consequently not investigated. In addition, cell-associated viruses of the blood will be detected to a less degree in sera compared to whole blood or buffy coat samples.

The cutoff value for the degree of similarity considered a hit used in the BLAST homology search is set relatively low to make sure that as many potential virus reads as possible are found. This, on the other hand, means that false hits may occur and, if a related virus is absent from the database, the probability for false hits increases. In addition, viruses not infecting the snow leopard but other microorganisms that in turn has infected the snow leopard will also end up as virus hits as well as viruses infecting prey consumed. All these types of presumably false hits have been manually excluded and are not listed in Table 1. The complete list of all virus hits is found in the Supplementary Material.

Reads most similar to RVA are by far the most common viral reads in the fecal swab samples. Many RVA reads were found in animals sampled both in 2011 and 2012, and it appears likely that rotavirus infections are common among snow leopards. Rotaviruses are important enteric pathogens both for humans and animals (34, 35). Reports on rotavirus infections in wild felids are lacking, and the scarce reports in domestic cats have not indicated serious impacts of the disease (36), although exceptions exist (37). In a large epidemiological investigation of RVA in cats housed within 25 rescue catteries across the United Kingdom, only 3% of the 1,727 cats sampled were infected (36). Although only eight snow leopards were sampled, the present study shows that the RVA prevalence of ~67% in snow leopards in the Tost Mountains of Mongolia is surprisingly high compared to existing knowledge. The genotypes of the NSP1 and VP3 genes of RVA could be determined and were found to be A3 and M2. These genotypes have also been found for RVA infecting domestic cats. The significance of RVA infections for the health status of snow leopards remains to be determined. It is known that rotaviruses frequently switch host and can establish productive infections leading to new strains in the novel host (38). It cannot be excluded that the rotaviruses observed here originate from prey either via a host switch event or alternatively from feeding on infected prey. The similarity of the snow leopard rotaviral sequences to bovine rotaviruses rather than feline rotaviruses may indicate that.

Herpesvirus is known to occur in wild felids (39) and can in particular affect the health and survival of juveniles (40). Feline herpesvirus type 1 (FeHV-1) is an important cause of acute upper respiratory tract and ocular disease in cats (41). In the present study, a 289-nt-long contig assembled from a serum sample from animal M1 and a 228-nt-long read from the rectal swab from animal F9 were found to be identical or nearly identical to segments from the genome of the FeHV-1 strain KANS_02. This represents to our knowledge the first indication that FeHV-1 circulates among snow leopards.

It was observed some 20 years ago that different patterns were observed with immunohistochemical screening for oral and cutaneous papillomavirus-induced lesions both for snow leopards kept in captivity and domestic cats, indicating the presence of two distinct papillomavirus species (42). The only snow leopard papillomavirus that has been genetically characterized (UuPV-1) is similar to feline papillomavirus type 1 (43). However, the papillomavirus discovered in the present study is more similar to feline papillomavirus type 2 (Figure 3). There is increasing evidence that feline papillomavirus type 2 plays a significant role in the development of skin cancers of domestic cats (44) and could potentially also cause skin cancer among snow leopards.

With application of NGS, a multitude of small circular DNA viruses have been discovered (45–47). For example, many single-stranded DNA (ssDNA) viruses encoding a replication-associated protein have been discovered, and since they by phylogenetic analysis appear to belong to yet uncharacterized but related virus families, they are collectively referred to as CRESS DNA viruses (46, 47). Among the CRESS DNA viruses detected in vertebrates, the Circoviridae and Genomoviridae are well-established and prevalent. Viruses in the Parvoviridae family, which can cause the serious disease feline infectious enteritis (48), have small linear single-stranded DNA genome and thus do not belong to the CRESS DNA virus group. But they have a similar genome organization, and all encode a homologous Rep-protein. Besides the CRESS DNA viruses, also members of the Anelloviridae are circular single-stranded DNA viruses commonly found infecting vertebrates. In contrast to the CRESS DNA viruses, anelloviruses lack a Rep-coding gene (49).

In the present study, we observed blast hits with best matches to the putative parvovirus NIH-CQV. However, these viral sequences have been shown to be contaminants from the Qiagen extraction spin columns (50) and will not be further discussed. No other parvovirus blast hits were found.

Similar to many other studies [e.g., (25)] of fecal microbiome, many reads from small circular DNA viruses were observed in the present study (Table 1). In particular, virus Chimp162 (25) is prevalent and is found in four of seven fecal swab samples. Another significant finding of small circular DNA virus is feline anellovirus found in one animal in this study. The clinical significance of these types of viruses has however not yet been well-established.

Finally, it should be noted that few RNA virus reads are found in the present study. The main exception is rotavirus, which is a double-stranded RNA virus known to be very stable and resistant in the environment (51, 52). These observations could be interpreted as the RNA to a large extent has not been well-preserved on the rectal swab samples. The rectal swabs were maintained at ambient temperature during several months at the site of collection before they were transported to locations where they could be kept in freezers. Further sampling taking precautions for RNA preservation, for instance by using PrimeStore^® molecular transport medium (Longhorn™ Vaccines and Diagnostics, LLC) or RNAlater™, might be necessary to investigate the presence of other potentially important RNA viruses in snow leopards.

Several of the viruses indicated in the present study such as feline herpesvirus 1, feline papillomavirus 2, and RVA could affect the health of snow leopards (37, 41, 44). The impact that these viruses may have on the health of the snow leopards largely depends on the host immunity, which in turn is modulated by other factors such as availability of prey, suitable habitat, and stress. Thus, in sensitive animals that are under environmental stress, for example, young dispersing individuals and lactating females, any acquired health issues may be exacerbated by latent virus infections, which for example papillomavirus, and herpesvirus are known to establish.

A complication of this study was the isolated location where the animals reside, prohibiting optimal sample maintenance. This fact will bias the viruses observed toward those with DNA genomes and those that possess a more stable capsid structure that protects the nucleic acids better. With these limitations in mind, this is nonetheless a first step toward establishing a baseline characterizing the viral pathogens most commonly found in snow leopards. It is important to continue the sampling in our study area and to expand into other parts of the snow leopard distribution range.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov, PRJNA626440.

Ethics Statement

Capture of the snow leopards followed appropriate ethical standards and ethical approval for the study was obtained from Mongolia's Ministry of Environment and Tourism (Government Building II, United Nations Street 5/2, Ulaanbaatar 15160, Mongolia).

Author Contributions

ÖJ planned the study and carried out the sampling together with PL and JM. KU and MW carried out sample preparations and next-generation sequencing. ML performed bioinformatic analysis of sequence data. ML, ÖJ, and KU wrote the manuscript that was revised by JM and MW. All authors read and approved the final manuscript.

Funding

Funding for the study was obtained through Helsinki Zoo and Foundation Segré, David Shepherd Wildlife Foundation, Disney Conservation Foundation, Kolmården Zoo, Nordens Ark, Columbus Zoo and Aquarium, Disney Conservation Foundation, and individual donations. None of the funding bodies had any role in the design of the study, collection, preparation, or analyses of data or publication of the manuscript.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We are thankful to the Ministry for Environment and Tourism, Government of Mongolia, and the Mongolian Academy of Sciences for supporting our snow leopard work. We also want to thank Dr. Gustaf Samelius for useful discussions.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2020.00645/full#supplementary-material

References

1. Deem SL, Karesh WB, Weisman W. Putting theory into practice: wildlife health in conservation. Conserv Biol. (2001) 15:1224–33. doi: 10.1046/j.1523-1739.2001.00336.x

CrossRef Full Text | Google Scholar

2. Woodroffe R. Managing disease threats to wild mammals. Anim Conserv. (1999) 2:185–93. doi: 10.1111/j.1469-1795.1999.tb00064.x

CrossRef Full Text | Google Scholar

3. Chapron G, Miquelle DG, Lambert A, Goodrich JM, Legendre S, Clobert J. The impact on tigers of poaching versus prey depletion. J Appl Ecol. (2008) 45:1667–74. doi: 10.1111/j.1365-2664.2008.01538.x

CrossRef Full Text | Google Scholar

4. Ripple WJ, Estes JA, Beschta RL, Wilmers CC, Ritchie EG, Hebblewhite M, et al. Status and ecological effects of the world's largest carnivores. Science. (2014) 343:1241484. doi: 10.1126/science.1241484

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Packer C, Pusey AE. Adaptions of female lions to infanticide by incoming males. Am Nat. (1983) 121:716–28. doi: 10.1086/284097

CrossRef Full Text | Google Scholar

6. Swenson JE. Implications of sexually selected infanticide for the hunting of large carnivores. In: Festa-Bianchet M, Apollonio M, editors. Animal Behavior and Wildlife Conservation. Washington, DC: Island Press. (2003). p. 171–91.

Google Scholar

7. Roelke-Parker ME, Munson L, Packer C, Kock R, Cleaveland S, Carpenter M, et al. A canine distemper virus epidemic in serengeti lions (Panthera leo). Nature. (1996) 379:441–5. doi: 10.1038/379441a0

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Jackson RM, Mishra C, McCarthy TM, Ale SB. Snow leopards: conflict and conservation. In: Macdonald DW, Loveridge AJ, editors. Biology and Conservation of Wild. Felids: Oxford University Press. (2010). p. 417–30.

Google Scholar

9. Johansson Ö, Koehler G, Rauset GR, Samelius G, Andrén H, Mishra C, et al. Sex-specific seasonal variation in puma and snow leopard home range utilization. Ecosphere. (2018) 9:e02371. doi: 10.1002/ecs2.2371

CrossRef Full Text | Google Scholar

10. Johansson Ö, Simms A, McCarthy T. From VHF to satellite GPS collars – advancements in snow leopard telemetry. In: McCarthy T, Mallon D, editors. Snow Leopards. Biodiversity of the World: Conservation From Genes to Landscapes. 1st ed. Elsevier: Academic Press. (2016). p. 355–65.

Google Scholar

11. McCarthy T, Mallon D, Sanderson EW, Zahler P, Fisher K. What is a snow leopard? Biogeography and status overview. In: McCarthy T, Mallon D, editors. Snow leopards. Biodiversity of the World: Conservation From Genes to Landscapes. 1st ed. Elsevier: Academic Press. (2016). p. 23–42. doi: 10.1016/B978-0-12-802213-9.00003-1

CrossRef Full Text | Google Scholar

12. Blankinship JC, Niklaus PA, Hungate BA. A meta-analysis of responses of soil biota to global change. Oecologia. (2011) 165:553–65. doi: 10.1007/s00442-011-1909-0

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Murray DL, Kapke CA, Evermann JF, Fuller TK. Infectious disease and the conservation of free-ranging large carnivores. Anim Conserv. (1999) 2:241–54. doi: 10.1111/j.1469-1795.1999.tb00070.x

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Ostrowski S, Gilbert M. Diseases of free-ranging snow leopards and primary prey species. In: McCarthy T, Mallon D, editors. Snow Leopards. Biodiversity of the World: Conservation From Genes to Landscapes. 1st ed. Elsevier: Academic Press. (2016). p. 97–112.

Google Scholar

15. Esson C, Skerratt LF, Berger L, Malmsten J, Strand T, Lundkvist Å, et al. Health and zoonotic Infections of snow leopards Panthera unica in the south gobi desert of Mongolia. Infect Ecol Epidemiol. (2019) 9:1604063. doi: 10.1080/20008686.2019.1604063

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Johansson Ö, Malmsten J, Mishra C, Lkhagvajav P, McCarthy T. Reversible immobilization of free-ranging snow leopards (Panthera uncia) with a combination of medetomidine and tiletamine-zolazepam. J Wildl Dis. (2013) 49:338–46. doi: 10.7589/2012-02-049

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Sharma K, Bayrakcismith R, Tumursukh L, Johansson Ö, Sevger P, McCarthy T, et al. Vigorous dynamics underlie a stable population of the endangered snow leopard Panthera uncia in tost mountains, South Gobi, Mongolia. PLoS ONE. (2014) 9:e101319. doi: 10.1371/journal.pone.0101319

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Johansson Ö, McCarthy T, Samelius G, Andrén H, Tumursukh L, Mishra C. Snow leopard predation in a livestock dominated landscape in Mongolia. Biol Conserv. (2015) 184:251–8. doi: 10.1016/j.biocon.2015.02.003

CrossRef Full Text | Google Scholar

19. Allander T, Emerson SU, Engle RE, Purcell RH, Bukh J. A virus discovery method incorporating DNase treatment and its application to the identification of two bovine parvovirus species. Proc Natl Acad Sci USA. (2001) 98:11609–14. doi: 10.1073/pnas.211424698

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Babraham Bioinformatics. Fast, Streaming QA/QC for High Throughput Sequencing Data. (2018). Available online at: https://github.com/s4hts/HTStream

21. Maes P, Matthijnssens J, Rahman M, Van Ranst M. RotaC: a web-based tool for the complete genome classification of group A rotaviruses. BMC Microbiol. (2009) 9:238. doi: 10.1186/1471-2180-9-238

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Virus Pathogen Resource (ViPR). (2020). Available online at: http://www.viprbrc.org (accessed June 9, 2019).

23. Pickett BE, Sadat EL, Zhang Y, Noronha JM, Squires RB, Hunt V, et al. ViPR: an open bioinformatics database and analysis resource for virology research. Nucleic Acids Res. (2012) 40:D593–8. doi: 10.1093/nar/gkr859

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. (1990) 215:403–10. doi: 10.1016/S0022-2836(05)80360-2

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Li L, Kapoor A, Slikas B, Bamidele OS, Wang C, Shaukat S, et al. Multiple diverse circoviruses infect farm animals and are commonly found in human and chimpanzee feces. J Virol. (2010) 84:1674–82. doi: 10.1128/JVI.02109-09

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Dayaram A, Galatowitsch M, Harding JS, Arguello-Astorga GR, Varsani A. Novel circular DNA viruses identified in procordulia grayi and xanthocnemis zealandica larvae using metagenomic approaches. Infect Genet Evol. (2014) 22:134–41. doi: 10.1016/j.meegid.2014.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Dayaram A, Galatowitsch ML, Arguello-Astorga GR, van Bysterveldt K, Kraberger S, Stainton D, et al. Diverse circular replication-associated protein encoding viruses circulating in invertebrates within a lake ecosystem. Infect Genet Evol. (2016) 39:304–16. doi: 10.1016/j.meegid.2016.02.011

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Shan T, Li L, Simmonds P, Wang C, Moeser A, Delwart E. The fecal virome of pigs on a high-density farm. J Virol. (2011) 85:11697–708. doi: 10.1128/JVI.05217-11

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Fenner F. Adventures with poxviruses of vertebrates. FEMS Microbiol Rev. (2000) 24:123–33. doi: 10.1111/j.1574-6976.2000.tb00536.x

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Mallon D, Harris RB, Wegge P. Snow leopard prey and diet. In: McCarthy T, Mallon D, editors. Snow Leopards. Biodiversity of the World: Conservation From Genes to Landscapes. 1st ed. Elsevier: Academic Press. (2016).

Google Scholar

31. Delwart EL. Viral metagenomics. Rev Med Virol. (2007) 17:115–31. doi: 10.1002/rmv.532

CrossRef Full Text | Google Scholar

32. Granberg F, Karlsson OE, Leijon M, Liu L, Belak S. Molecular approaches to recognize relevant and emerging infectious diseases in animals. Methods Mol Biol. (2015) 1247:109–24. doi: 10.1007/978-1-4939-2004-4_7

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Mokili JL, Rohwer F, Dutilh BE. Metagenomics and future perspectives in virus discovery. Curr Opin Virol. (2012) 2:63–77. doi: 10.1016/j.coviro.2011.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Dodet B, Heseltine E, Saliou P. Rotaviruses in human and veterinary medicine. Trends Microbiol. (1997) 5:176–8. doi: 10.1016/S0966-842X(97)01045-7

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Martella V, Banyai K, Matthijnssens J, Buonavoglia C, Ciarlet M. Zoonotic aspects of rotaviruses. Vet Microbiol. (2010) 140:246–55. doi: 10.1016/j.vetmic.2009.08.028

PubMed Abstract | CrossRef Full Text | Google Scholar

36. German AC, Iturriza-Gomara M, Dove W, Sandrasegaram M, Nakagomi T, Nakagomi O, et al. Molecular epidemiology of rotavirus in cats in the United Kingdom. J Clin Microbiol. (2015) 53:455–64. doi: 10.1128/JCM.02266-14

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Phan TG, Leutenegger CM, Chan R, Delwart E. Rotavirus I in feces of a cat with diarrhea. Virus Genes. (2017) 53:487–90. doi: 10.1007/s11262-017-1440-4

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Matthijnssens J, Ciarlet M, Heiman E, Arijs I, Delbeke T, McDonald SM, et al. Full genome-based classification of rotaviruses reveals a common origin between human Wa-Like and porcine rotavirus strains and human DS-1-like and bovine rotavirus strains. J Virol. (2008) 82:3204–19. doi: 10.1128/JVI.02257-07

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Spencer JA. Survey of antibodies to feline viruses in free-ranging lions. South Afr J Wildlife Res. (1991) 21:59–61.

Google Scholar

40. Gaskell R, Willoughby K. Herpesviruses of carnivores. Vet Microbiol. (1999) 69:73–88. doi: 10.1016/S0378-1135(99)00092-9

CrossRef Full Text | Google Scholar

41. Gaskell R, Dawson S, Radford A, Thiry E. Feline herpesvirus. Vet Res. (2007) 38:337–54. doi: 10.1051/vetres:2006063

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Sundberg JP, Van Ranst M, Montali R, Homer BL, Miller WH, Rowland PH, et al. Feline papillomas and papillomaviruses. Vet Pathol. (2000) 37:1–10. doi: 10.1354/vp.37-1-1

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Rector A, Lemey P, Tachezy R, Mostmans S, Ghim SJ, Van Doorslaer K, et al. Ancient papillomavirus-host co-speciation in felidae. Genome Biol. (2007) 8:R57. doi: 10.1186/gb-2007-8-4-r57

PubMed Abstract | CrossRef Full Text

44. Munday JS, Sharp CR, Beatty JA. Novel viruses: update on the significance of papillomavirus infections in cats. J Feline Med Surg. (2019) 21:409–18. doi: 10.1177/1098612X18808105

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Shulman LM, Davidson I. Viruses with circular single-stranded DNA genomes are everywhere! Annu Rev Virol. (2017) 4:159–80. doi: 10.1146/annurev-virology-101416-041953

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Zhao L, Rosario K, Breitbart M, Duffy S. Eukaryotic circular rep-encoding single-stranded DNA (CRESS DNA) viruses: ubiquitous viruses with small genomes and a diverse host range. Adv Virus Res. (2019) 103:71–133. doi: 10.1016/bs.aivir.2018.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Simmonds P, Adams MJ, Benko M, Breitbart M, Brister JR, Carstens EB, et al. Consensus statement: virus taxonomy in the age of metagenomics. Nat Rev Microbiol. (2017) 15:161–8. doi: 10.1038/nrmicro.2016.177

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Stuetzer B, Hartmann K. Feline parvovirus infection and associated diseases. Vet J. (2014) 201:150–5. doi: 10.1016/j.tvjl.2014.05.027

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Zhang W, Wang H, Wang Y, Liu Z, Li J, Guo L, et al. Identification and genomic characterization of a novel species of feline anellovirus. Virol J. (2016) 13:146. doi: 10.1186/s12985-016-0601-8

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Xu B, Zhi N, Hu G, Wan Z, Zheng X, Liu X, et al. Hybrid DNA virus in Chinese patients with seronegative hepatitis discovered by deep sequencing. Proc Natl Acad Sci USA. (2013) 110:10264–9. doi: 10.1073/pnas.1303744110

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Keswick BH, Pickering LK, DuPont HL, Woodward WE. Survival and detection of rotaviruses on environmental surfaces in day care centers. Appl Environ Microbiol. (1983) 46:813–6. doi: 10.1128/AEM.46.4.813-816.1983

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Estes MK, Graham DY, Smith EM, Gerba CP. Rotavirus stability and inactivation. J Gen Virol. (1979) 43:403–9. doi: 10.1099/0022-1317-43-2-403

PubMed Abstract | CrossRef Full Text | Google Scholar