key: cord-0689196-rr4rcc1x
authors: Kohl, Claudia; Brinkmann, Annika; Radonić, Aleksandar; Dabrowski, Piotr Wojtek; Mühldorfer, Kristin; Nitsche, Andreas; Wibbelt, Gudrun; Kurth, Andreas
title: The virome of German bats: comparing virus discovery approaches
date: 2021-04-01
journal: Sci Rep
DOI: 10.1038/s41598-021-86435-4
sha: 197c43e048522323c5c4265cbae87dce8ed7bf91
doc_id: 689196
cord_uid: rr4rcc1x

Bats are known to be reservoirs of several highly pathogenic viruses. Hence, the interest in bat virus discovery has been increasing rapidly over the last decade. So far, most studies have focused on a single type of virus detection method, either PCR, virus isolation or virome sequencing. Here we present a comprehensive approach in virus discovery, using all three discovery methods on samples from the same bats. By family-specific PCR screening we found sequences of paramyxoviruses, adenoviruses, herpesviruses and one coronavirus. By cell culture we isolated a novel bat adenovirus and bat orthoreovirus. Virome sequencing revealed viral sequences of ten different virus families and orders: three bat nairoviruses, three phenuiviruses, one orbivirus, one rotavirus, one orthoreovirus, one mononegavirus, five parvoviruses, seven picornaviruses, three retroviruses, one totivirus and two thymoviruses were discovered. Of all viruses identified by family-specific PCR in the original samples, none was found by metagenomic sequencing. Vice versa, none of the viruses found by the metagenomic virome approach was detected by family-specific PCRs targeting the same family. The discrepancy of detected viruses by different detection approaches suggests that a combined approach using different detection methods is necessary for virus discovery studies.

. Comparison of results from virus discovery studies conducted. The viral sequences obtained by family-specific PCR are marked with "A" and viral sequences obtained by virome sequencing are marked with "B". Adenovirus and reoviruses were also isolated. The "C" indicates the finding of the respective viruses found by other studies. E the results from this study, NZ New Zealand. www.nature.com/scientificreports/ ity criteria. These reads could not be remapped to reference sequences, revealed non-plausible BLASTx/BLASTn results (i.e. genomic DNA of bats) or shared 100 percent identity to sequences obtained before (possible crosscontamination). Supplementary Table S2 summarizes the number of filtered reads (length/quality) obtained per pool with the number of viral reads (w/o phages) finally checked for quality and allocated by Diamond 38 , MEGAN6 and further analysis. In Fig. 2 the overall number of viral reads in all pools per virus family is summarized.

Viruses of nine different families and orders (Parvoviridae, Picornaviridae, Totiviridae, Mononegavirales, Reoviridae, Bunyavirales, Tymovirus, Retroviridae and phages) were confirmed within the bat samples by virome sequencing and further analysis. For each family, selected viruses were retested by specific PCRs in the individual bats' organs. For the viruses of highest interest, phylogenetic reconstructions were calculated as described in the methods section.

All viruses were confirmed by sequence assembly and comparison to reference strains and to each other as described in the methods section. Most of the viruses were tested back with designed specific primers, PCR and Sanger sequencing. The individual length of obtained contigs, the accession number of the reference sequence and the pairwise identity on nt and aa level are summarized in Table 3 . The accession numbers of all obtained viral sequences are available in Table 4 . The full description of the results is given in the supplemental section.

Parvoviridae. Parvoviruses (n = 5) were found in pools 1, 2, 3 and 9 belonging to the subspecies Densovirinae (146,313 reads) and Parvovirinae (20,296 reads). The virus sequences were named after the order of appearance and in relation to reference strains: Blatella Germanica densovirus-like virus 1, Blatella Germanica densovirus- Table 2 . Results PCR screening. Bats were screened for the presence of virus nucleic acids by different PCR assays. pan generic family-specific assay, qPCR quantitative real-time PCR, cPCR conventional specific PCR. *#/novel/size: #, number of positive samples; novel, number of novel viruses obtained; size, bat sample size. H Screened in collaboration with our Hungarian colleagues (see Ref. 33 ). Fig. 7 . Laurel Lake-like viruses 1 and 2 with other members of phenuiviruses is displayed in the supplemental section ( Supplementary Fig. S3 ).

Tymovirus. Tymoviruses (n = 2) were found in pools 3 and 5 (423 reads). The found virus sequences share the highest identity with the strains Bombyx mori latent virus and Grapevine Red Globe virus. The virus sequences were named after the order of appearance and in relation to reference strains: Bombyx mori latent-like virus and Grapevine Red Globe-like virus. Obtained sequences are available under the following accession numbers: MN851298 and MN851299. 

Virome sequencing of German bat carcasses of 16 bat species resulted in a high variety of confirmed viral sequences of parvoviruses, picornaviruses, totiviruses, reoviruses, nairoviruses, phenuiviruses, tymoviruses, retroviruses and several phages. In addition to the confirmed viral sequences, several sequencing reads were identified that shared a high homology and identity to other viruses. These high-identity reads were excluded as they were very likely false-positive results (i.e. poxviruses) (Fig. 1) . The occurrence of false-positive results in virome sequencing is well known and has been described before 39 . Table 1 compares the results obtained from different bat virome studies with our findings. Virome profiles found in our study are generally comparable to those in other studies. Several Picornaviruses, Parvoviruses, Retroviruses, Tymoviruses and Totiviruses were identified in the nine distinct virome profiles of European bats (Supplementary Tables S3, S4 , S5, S6, S7, S8, S9, S10 and S11). However, some of the viruses detected in our study are particularly interesting as they are phylogenetically closely related to viruses that can cause diseases in humans, or they are the first description of these viruses in www.nature.com/scientificreports/ certain bat species or within the European geographical range. The following discussion gives a summary of the viruses of highest interest. The full discussion on these viruses can be found in the supplemental section. Six Bunyaviruses were identified in this study (nairoviruses and phenuiviruses). Nairoviridae are a family within the order Bunyavirales and named after the type species Nairobi sheep disease virus 40 . The majority of nairoviruses is transmitted by ticks and several are capable of causing severe diseases in humans and animals 40 . Bat nairoviruses have been reported before: Ahun nairovirus (KF170224) and Gossas virus (KR534878) 21, 28, 41 . Other nairoviruses of bats have been described and form two monophyletic genogroups within the nairoviruses, Keterah and Kasokero 41 .

In this study three novel nairoviral sequences were detected and confirmed in tissues from German bats; these are related to Issyk-Kul virus (Id 95% nt; 99% aa), Sapphire II virus (Id 85% nt; 54% aa) and Avalon Bres virus (Id 71% nt; 50% aa). Nine out of twelve Eptesicus nilssonii bats in pool 3 were infected with a yet undescribed Issyk-Kul virus strain (Supplementary Table S5 ). The phylogenetic reconstruction clearly allocated this novel strain to already described Issyk-Kul viruses within the Keterah genogroup (Fig. 6) 41 . Issyk-Kul virus had first been isolated in 1970 from a Nyctalus noctula bat in Kyrgyzstan and later on in Tajikistan and Kazakhstan 42, 43 . Issyk-Kul virus was described to cause sporadic febrile outbreaks in humans with headache, myalgia and nausea 42, 44 . It is assumed that Issyk-Kul virus can be transmitted by tick bites and exposure to bat feces and urine, eventually 42 . The Issyk-Kul-like virus described here was found predominantly in liver, spleen and lung tissues of the respective bats, indicating systemic infection of bats instead of solely passaging intestinal tick content. We named the novel Issyk-Kul virus strain "PbGER" after its origin in Prackenbach, Germany ( Table 4 ). The novel Issyk-Kul virus Table S7 ). Phylogenetic reconstruction indicated that Sapphire II-like virus is related to Sapphire II virus and clusters with the Dera Ghazi Khan genogroup usually associated with birds 41 (Fig. 5 ). Sapphire II virus was isolated from swallowed ticks in 1972 and was not reported to cause any diseases in humans 46 . This is the first description of this genotype in bats. We named the Sapphire II-like virus Berlin bat nairovirus (BbnV) after the origin of the corresponding bats.

Avalon Bres virus-like sequence was detected in Pipistrellus pipistrellus bats of pool 6 (Supplementary Table S8 ). Phylogenetically, Avalon Bres virus clusters monophyletically with the Sakhalin genogroup. Viruses of these genogroups have not been described before to be associated with bats. However, a serological study showed that several wild-caught bats had antibody responses to CCHFV proteins 47 . We named the Avalon Bres virus-like virus Wittenau bat nairovirus (WbnV) after the origin of the corresponding bats.

Phenuivirus is a family within the order Bunyavirales 40 . The majority of the ten phenuivirus genera is mosquito-borne; however, some genera are transmitted by ticks (i.e. Banyangviruses) and are capable of causing severe diseases in humans and animals. Three phenuiviruses (genus phlebovirus) have been reported to be identified from bats: Malsoor virus, Rift Valley virus and Toscana virus 42, 43 . Malsoor virus was isolated from Rousettus leschenaultii in India and is by phylogenetic reconstruction monophyletic with viruses of the genus Reconstruction was performed via the Bayesian MCMC approach using MrBayes with the following settings (burn-in, 30%; frequency, 100; chain length, 1 million to 10 million, depending on when a standard derivation of below 0.025 was reached) 78 . Hantaan virus was used as an outgroup. Reconstructed trees were visualized using FigTree and posterior probabilities were depicted for each node (http:// tree. bio. ed. ac. uk/ softw are/ figtr ee/). (Fig. 7) . This evolutionary distance could indicate a potential zoonotic transmission of both Malsoor virus and the found Malsoor-like virus to humans. We named the Malsoor-like virus Zwiesel bat banyangvirus (ZbbV) after the origin of the corresponding bats. The novel Zwiesel bat phlebovirus is further analyzed by whole genome sequencing as well as throughout analysis; this is the subject of another study published simultaneously as a spin-off to this study 49 .

Laurel Lake-like virus 1 and Laurel Lake-like virus 2 (genus Laulavirus) were identified in bats from pool 8 (Supplementary Table S10 ). Pool 8 comprises bats of two species, Pipistrellus nathusii and Plecotus aureus. Phylogenetic reconstruction of the Laurel Lake-like virus showed that Laurel Lake virus and Laurel Lake-like virus are quite distanced from the other viruses of the Uukuniemi group, although they are clearly clustering ( Supplementary Fig. S3) 50 . We named the Laurel Lake-like virus 1 Bavarian bat lalavirus (BblV) after the origin of the corresponding bats. We named the Laurel Lake-like virus 2 Munich bat lalavirus (MblV) after the origin of the corresponding bats. Reconstruction was performed via the Bayesian MCMC approach using MrBayes with the following settings (burn-in, 30%; frequency, 200; chain length, 1 million to 10 million, depending on when a standard derivation of below 0.025 was reached) 78 . Hantaan virus was used as an outgroup. Reconstructed trees were visualized using FigTree and posterior probabilities were depicted for each node (http:// tree. bio. ed. ac. uk/ softw are/ figtr ee/). www.nature.com/scientificreports/ Three reovirus sequences (genera orbivirus, rotavirus and orthoreovirus) were identified in the virome data of pools 2, 5 and 8 (Supplementary Tables S4, S5 and S10). The orthoreovirus in pool 8 has been known before (T3/Bat/Germany/342/08) and served as a kind of positive control for this study, as we added the bat tissue from bat 342/08, from which the virus was initially isolated, to the pool 1 (Supplementary Table S10 ). The orbivirus and rotavirus are genera within the subfamily Sedoreovirinae of the family Reoviridae. The orbiviruses comprise several species that are inducing severe diseases in humans and animals (i.e. bluetongue disease or epizootic hemorrhagic disease virus) and are capable of replication in several arthropod and vertebrate hosts 51 . In pool 2 an orbivirus related to the yet unpublished Bat orbivirus from China (AccNo. MH144554.1) (Id 81% aa) and Sathuvachari virus was identified in Nyctalus noctula bats (Supplementary Table S4) (Fig. 3) . Sathuvachari virus has first been isolated in India in 1963, and because of its orbivirus character has tentatively been classified as mosquito-borne although the strain was isolated from a bird (starling) 52 . There have been two further orbiviruses identified from bats (Myotis rickettsi and Rhinolphus ferrumequinum) in China (Acc. No. KX343070.1 and KX161703.1). We named the Sathuvachari-like virus Irlbach bat orbivirus (IboV) after the origin of the corresponding bats.

The bat rotavirus identified in this study was detected in pool 5 which comprises tissues from Pipistrellus pipistrellus bats (Supplementary Table S7 ). Phylogenetic reconstruction allocates the bat rotavirus sequence into a distinct but related clade to rotaviruses type A (Id 58% nt; 57% aa) (Fig. 4) . Numerous rotaviruses have been identified in bats but only one other rotavirus has been identified in European bats 28 . The bat rotavirus from France is similarly related to rotaviruses of group A. The zoonotic potential of these bat rotaviruses related to Figure 7 . Phylogenetic reconstruction of Phenuiviruses glycoprotein-ML, 1 million, 2578 nt. Reconstruction was performed via the Bayesian MCMC approach using MrBayes with the following settings (burn-in, 30%; frequency, 100; chain length, 1 million to 10 million, depending on when a standard derivation of below 0.025 was reached) 78 . Reconstructed trees were visualized using FigTree and posterior probabilities were depicted for each node (http:// tree. bio. ed. ac. uk/ softw are/ figtr ee/). Table S4 ). It was described for these bats to have parasites in their intestines 53 ; it is possible that the sequence originated from a tapeworm. On amino acid level the novel sequence shared highest similarity to Wenzhou tapeworm 1 virus (Id 35% aa) and Midway virus 54, 55 . Phylogenetically the strain was grouped distinct from Midway virus and was more closely related to Borna disease virus and rabies virus. The similarities on nt level are quite low and the tree has to be interpreted with caution (SF2). However, these findings are interesting as bats are often discussed as potential reservoir host of many viruses of the order Mononegavirales (i.e. Ebola virus, Marburg virus, Nipah virus, Hendra virus and lyssaviruses). We named the Wenzhou tapeworm 1-like virus Berlin bat mononegavirus (BbmV) after the origin of the corresponding bats.

The screening of 375 bats with numerous family-specific PCRs resulted in several novel sequences ( Table 2 ). The inoculation of the bat tissues into six different cell lines revealed two novel isolates (Supplementary Table S1 ). Virome sequencing of 189 histopathologically suspicious tissues of these bats resulted again in a high variety of novel virus sequences. Taking all of this into account, it is somewhat surprising to find different viruses in identical samples with these complementary approaches. However, the setup of the study design is insufficient to provide reliable statistics. For the RNA virus screening the organ material was homogenized in RNAlater ® buffer and 30 µl per organ were pooled before extraction. For the DNA virus screening 30 µl of each organ natively frozen at − 80 °C was pooled before extraction. Cell culture isolation was also done with the native samples, although on differing cell culture systems. However, in cell culture only the occurrence of a visible CPE gives a clue of an ongoing virus infection and the absence of CPE does not exclude the potential infection. For the virome sequencing only those individual organs (native − 80 °C) were used that showed alterations in histopathology. If family-specific PCRs resulted in clear amplification and revealed a novel virus sequence after Sanger sequencing, one would expect to be able to obtain sequences of the same virus by virome sequencing of the same sample. In fact, many of the organ bat tissues that were prepared for virome sequencing also revealed positive results in family-specific PCRs. For example, Bat E95/09 (Supplementary Table S7 ) was tested positive for paramyxoviruses by PCR 32 . The respective tissues were subjected to virus purification and virome sequencing within pool 4, but no paramyxoviruses were detected in the virome data. A possible explanation could be that the purification method we have used, TUViD-VM, did not sufficiently purify paramyxoviruses from tissues; but this is exactly what this protocol was designed for 10 . Actually, the whole protocol was validated and tested with a range of viruses (i.e. paramyxovirus, poxvirus, influenza virus and reovirus), with the most remarkable results from paramyxovirus (Sendai virus) 10 . In addition, we have tested the reproducibility and consistency of the TUViD-VM protocol numerous times and found it to be very stable 56 . Furthermore, it has to be considered that we have used pooled tissues for our examinations instead of i.e. single organs. It is possible that this may have affected the detectability of the different viruses. However, we have used pooled tissues for all methods compared and therefore results should be comparable. Taking this into account, we propose that nothing is more important in virome studies than the ratio between targeted and untargeted sequences. For instance, if we apply a family-specific PCR to a mixture of cDNA containing 100 genome equivalents of the targeted paramyxovirus, we will most probably be able to obtain a clean and sufficient amplification of the targeted sequences, although the cDNA contains only the 100 genome equivalents of paramyxovirus over millions of other sequences. In comparison, the approach of virome sequencing differs significantly. When millions of other sequences compete (e.g. host genome) the likelihood of detection of an underrepresented number of viral sequences would be very low. To address and circumvent this obvious fact, the TUViD-VM protocol was applied in this study, as it aims to decrease the amount of background (host) sequences and amplify the viral sequences to increase the detection ratio and likelihood. One accompanying effect we observed was that high numbers of viral reads of other viruses (i.e. more than 900,000 reads for picornaviruses in pool 5) impeded the sequencing of low amounts of other viral amplicons. Similarly, the analyses of PCR-negative results revealed a surprising outcome. All bats were screened for the presence of bunyaviral sequences, especially for phenuiviruses and nairoviruses, by family-specific PCR with negative results. Interestingly, virome sequencing revealed the presence of six novel strains of bunyaviruses with identities ranging from 70 to 95 percent identity on nt level to known virus sequences.

This leads to the conclusion that the outcome of virus discovery studies, whether approached by PCR, metagenomics or cell culture, provides only a fraction of the viruses truly present in a sample. The methods should be used complementarily to give a more complete picture of virus diversity rather than to replace each other-as has often been described for the advantages of metagenomics over classic PCR approaches. The combined data of existing individual studies investigated by different virus discovery protocols with their conclusions are further used to extrapolate and predict outbreak risks and hotspots. As long as we use single virus detection protocols in studies only, the very true diversity of viruses will stay hidden, unrecognized by our insufficient methodological approaches. It might be time for common approaches, where scientists doing virus discovery and modeling work join forces and design studies in a more comparable and meaningful manner.

Study. In compliance with the species protection through the European Commission (https:// ec. europa. eu/ envir onment/ nature/ legis lation/ habit atsdi recti ve/) and the Agreement on the Conservation of Populations of European Bats (www. eurob ats. org), investigative research on bats was permitted by local government bodies 38 . We named the manuscript Virome of German Bats as the bats were found in Germany. However, bats are migrating animals, many of them covering long distances throughout Europe. Therefore, our results could also be expected in other parts of Europe. We herewith confirm that all permits to investigate carcasses of deceased bats and all experimental protocols on the dead bats were approved by the respective local governmental authori- 3] ). All bat carcasses were kindly provided by bat researchers and bat rehabilitation centers from the different geographic regions. In case bats died in care or had to be euthanized for medical reasons, the carcasses were handled as described before 38, [57] [58] [59] . Aliquots of the individual organs were divided between two tubes: one tube with RNAlater ® (Qiagen, Hilden, Germany) for RNA extraction and one tube native frozen at − 80 °C and sent to the Robert Koch Institute for virological examination (Fig. 8 ).

All work was performed at Biosafety level-2 conditions with appropriate precautions. The study was divided into three methodological sections. (1) PCR screening, (2) Virus isolation via cell culture and (3) Virome analysis. Every section had its own difficulties and individual steps for sample preparation explained in the following.

PCR screening for RNA viruses. 30 µl of each of the eight individual organ homogenates (in RNAlater ® ) of one respective bat were pooled (lungs, liver, spleen, heart, brain, salivary gland, intestine and kidney). RNA extraction from RNAlater ® -stabilized samples and cDNA synthesis was performed as described before 32 . The pool containing 240 µl of organ homogenate was then extracted using the PureLinkTM Viral RNA/DNA Mini Kit (Invitrogen, Carlsbad, USA). Extracted RNA was transcribed to cDNA by using the TaqMan ® Reverse Transcription Reagents (Applied Biosystems, Darmstadt, Germany). The cDNA pools of 375 bats were screened by differing family-specific PCR systems for the detection of paramyxoviruses 60 , arenaviruses 61 , coronaviruses 35 , filoviruses 62 , flaviviruses 63 , hantaviruses 64 , nairo-and phleboviruses 65 and influenza viruses 66 .

PCR screening for DNA viruses. 30 µl of each of the eight individual organ homogenates (native frozen at − 80 °C) of one respective bat were pooled (lungs, liver, spleen, heart, brain, salivary gland, intestine, kidney). DNA was extracted using the PureLinkTM Viral RNA/DNA Mini Kit (Invitrogen) as described before 31, 33 . The DNA pools of 375 bats were screened by differing family-specific PCR systems for the detection of adenoviruses 67 , herpesviruses 68 and poxviruses 69, 70 . www.nature.com/scientificreports/ Virus isolation in cell culture. Virus isolation in cell culture was performed from pools containing all eight individual organs of the respective bats natively frozen at − 80 °C as described before 30, 31 . For cell culture isolation different cell culture systems were used; these are listed in Supplementary Table S1 .

Virome sequencing. Selected single organs (native frozen at − 80 °C) of 118 of the 375 individual bats were subjected to further investigation. These 118 bats and their corresponding 189 native frozen organ tissues were chosen by analysis of the pathological findings 71 . Several bats were found to have alterations that may be related to viral infections (e.g., interstitial pneumonia, pulmonary BALT hyperplasia, infiltrates of mononuclear cells in different organs, diffuse enlarged villi and mononuclear intestinal cell infiltrates or catarrhal or hemorrhagic enteritis). Due to the large number of individual samples, organs were allocated to nine pools for virus purification with TUViD-VM and sequencing 10 . Pools were established using the following criteria: (a) Pools should consist of similar sample numbers, (b) Pools should have only one bat species and (c) If number of bats is not sufficient to build a single pool, species are mixed. Table 5 and Fig. 9 summarize the composition of these nine individual pools. 200 µl of each individual organ from the 118 bats was pooled into nine differing pools, depending on histopathological results and according to the species. The composition of the individual pools is given in Table 5 . Each of the nine pools was divided for inoculation onto Vero and Paki cell cultures for virus isolation as described before 72, 73 . The original TUViD-VM purification protocol used in this study comprises the following steps: (1) homogenization, (2) clearing centrifugation, (3) ultracentrifugation through sucrose cushion, (4) ultracentrifugation, (5) digestion, (6) extraction, (7) cDNA and double strand synthesis and (8) random amplification and size selection via agarose gel electrophoresis 10 . Instead of using 200 µl of tissue homogenate as prescribed in the original protocol, 3 ml of organ tissue was used and the protocol was upscaled with small modifications 10 : in most pools the homogenate volume exceeded 3 ml; however, in case the amount of homogenate was too low, the volume was adapted with PBS to 3 ml; for ultracentrifugation, the rotors and pace were adapted accordingly by using the SW32 rotor (Beckman Coulter, Krefeld, Germany) also for the first ultracentrifugation step. From step (4) of the original TUViD-VM protocol onward, the protocol was performed without any further modifications.

After size selection by excision of the area between 500 and 1000 nt following agarose gel electrophoresis, samples were prepared for sequencing. Libraries were built using the NexteraXT protocol as described before 10 . The sequencing reaction was performed on the Illumina HiSeq Sequencer in high throughput mode using 150 cycles.

Raw data was processed as follows ( Fig. 8 ): Adapters were trimmed and filtering of read length and quality was performed. Raw reads were then subtractive mapped against customized databases to remove background reads (bacteria and mammalian genomes). Remaining reads were assembled to contigs using Velvet 74 . Contigs were then compared to the viral protein database (nr NCBI) via the DIAMOND tool and results were visualized in MEGAN6 (-sensitive) 37, 38 . Selected results from MEGAN6 were blasted using BLASTx on the NCBI GenBank database. Full genomes (if available) of the most similar sequences were downloaded to serve as a reference. The initial trimmed reads were then mapped to the reference sequence using Bowtie2 to identify more matching reads and check the correct allocation visually for quality 75 . Virus genomes were also assembled using Bowtie2. Sequences were checked for plausibility by analyzing the closest related virus. In the case of i.e. retroviruses Table 5 . Composition of bat pools. Table summarizes species and organs for the nine individual pools used for metagenomics sequencing and cell culture isolation. Organs written in bold have the highest percentage in the individual pool. # number, L lung, I intestine, S spleen, SG salivary gland, H heart, Li liver, K kidney, B brain. www.nature.com/scientificreports/ the closest related sequences were also from bat species, proving that they are possible to occur. Sequences of highest interest were extended by spanning PCRs and primer walking followed by Sanger sequencing. For all sequences of interest specific primers were designed to confirm the sequences in cDNA by conventional PCR of the individual extracted pools and organ samples eventually. Cycling conditions and primers are available on request. Any bands of targeted height were purified, Sanger sequenced and compared to the original sequence obtained from Next Generation Sequencing (NGS). Final sequences obtained from data analysis and Sanger sequencing were used to reconstruct phylogenetic trees. Sequences were aligned to type species (whenever possible) of the respective viral family (ICTV database) and other related viruses of interest via ClustalW. Alignment quality was checked with the online tool T-Coffee 76 . Model of evolution was predicted via jmodeltest and the model with the best AIC score was picked 77 . The actual reconstruction was performed via the Bayesian MCMC approach using MrBayes with the following settings (burn-in, 30%; frequency, 100; chain length, 1 million to 10 million depending on when a standard derivation of below 0.025 was reached) 78 . Reconstructed trees were visualized using FigTree and posterior probabilities were depicted for each node (http:// tree. bio. ed. ac. uk/ softw are/ figtr ee/). www.nature.com/scientificreports/

Bats: Important reservoir hosts of emerging viruses

Isolation of Hendra virus from pteropid bats: A natural reservoir of Hendra virus

Isolation of genetically diverse Marburg viruses from Egyptian fruit bats

Bats are natural reservoirs of SARS-like coronaviruses

Emerging diseases. Malaysian researchers trace Nipah virus outbreak to bats

A molecular phylogeny for bats illuminates biogeography and the fossil record

Insights into the ancestral organisation of the mammalian MHC class II region from the genome of the pteropid bat

Evolution and comparative analysis of the bat MHC-I region

Comparative analysis of bat genomes provides insight into the evolution of flight and immunity

Protocol for metagenomic virus detection in clinical specimens

Laboratory procedures to generate viral metagenomes

Evaluation of rapid and simple techniques for the enrichment of viruses prior to metagenomic virus discovery

Metagenomic study of the viruses of African straw-coloured fruit bats: Detection of a chiropteran poxvirus and isolation of a novel adenovirus

Evidence for retrovirus and paramyxovirus infection of multiple bat species in China

Metagenomic analysis of viruses from bat fecal samples reveals many novel viruses in insectivorous bats in China

Discovery of novel virus sequences in an isolated and threatened bat species, the New Zealand lesser short-tailed bat (Mystacina tuberculata)

Identification of a novel bat papillomavirus by metagenomics

Bat guano virome: Predominance of dietary viruses from insects and plants plus novel mammalian viruses

Deciphering the bat virome catalog to better understand the ecological diversity of bat viruses and the bat origin of emerging infectious diseases

Metagenomic analysis of the viromes of three North American bat species: Viral diversity among different bat species that share a common habitat

Virome analysis of two sympatric bat species (Desmodus rotundus and Molossus molossus) in French Guiana

Virome analysis for identification of novel mammalian viruses in bat species from Chinese provinces

A single bat species in Cameroon harbors multiple highly divergent papillomaviruses in stool identified by metagenomics analysis

Viral metagenomics of six bat species in close contact with humans in southern China

Novel highly divergent sapoviruses detected by metagenomics analysis in straw-colored fruit bats in Cameroon

Virome profiling of bats from Myanmar by metagenomic analysis of tissue samples reveals more novel mammalian viruses

Virome analysis for identification of novel mammalian viruses in bats from Southeast China

A preliminary study of viral metagenomics of French bat species in contact with humans: Identification of new mammalian viruses

The database of bat-associated viruses

Isolation and characterization of three mammalian orthoreoviruses from European bats

Genome analysis of bat adenovirus 2: Indications of interspecies transmission

Novel paramyxoviruses in free-ranging European bats

Random sampling of the Central European bat fauna reveals the existence of numerous hitherto unknown adenoviruses

Discovery of herpesviruses in bats

Generic detection of coronaviruses and differentiation at the prototype strain level by reverse transcription-PCR and nonfluorescent low-density microarray

New adenovirus in bats

MEGAN analysis of metagenomic data

Fast and sensitive protein alignment using DIAMOND

False-positive results in metagenomic virus discovery: A strong case for follow-up diagnosis

Bunyavirus taxonomy: Limitations and misconceptions associated with the current ICTV criteria used for species demarcation

A global genomic characterization of nairoviruses identifies nine discrete genogroups with distinctive structural characteristics and host-vector associations

Issyk-Kul" virus, a new arbovirus isolated from bats and Argas (Carios) vespertilionis (Latr., 1802) in the Kirghiz S

Complete genomic sequence of Issyk-Kul virus

Outbreak of arbovirus infection in the Tadzhik SSR due to the Issyk-Kul virus (Issyk-Kul fever)

First detection of bat-borne Issyk-Kul virus in

The Dictionary of Virology

Evidence for widespread infection of African bats with Crimean-Congo hemorrhagic fever-like viruses

Malsoor virus, a novel bat phlebovirus, is closely related to severe fever with thrombocytopenia syndrome virus and heartland virus

Zwiesel bat banyangvirus, a potentially zoonotic Huaiyangshan banyangvirus (formerly known as SFTS)-like banyangvirus in Northern bats from

Identification of novel viruses in Amblyomma americanum, Dermacentor variabilis, and Ixodes scapularis ticks

The diversity of the orthoreoviruses: molecular taxonomy and phylogentic divides

A novel mosquito-borne Orbivirus species found in South-east Asia

Diseases and causes of death in European bats: Dynamics in disease susceptibility and infection rates

Redefining the invertebrate RNA virosphere

Isolation and characterization of midway virus: A new tick-borne virus related to Nyamanini

Use of RNALater((R)) preservation for virome sequencing in outbreak settings

Genetic diversity of Pasteurella species isolated from European vespertilionid bats

Yersinia species isolated from bats

Proposal of Vespertiliibacter pulmonis gen. nov., sp. Nov. and two genomospecies as new members of the family Pasteurellaceae isolated from European bats

Sensitive and broadly reactive reverse transcription-PCR assays to detect novel paramyxoviruses

RT-PCR assay for detection of Lassa virus and related Old World arenaviruses targeting the L gene

Rapid molecular strategy for filovirus detection and characterization

Generic RT-nested-PCR for detection of flaviviruses using degenerated primers and internal control followed by sequencing for specific identification

Hantavirus in African wood mouse

Consensus amplification and novel multiplex sequencing method for S segment species identification of 47 viruses of the Orthobunyavirus, Phlebovirus, and Nairovirus genera of the family Bunyaviridae

Diagnostic approach for the differentiation of the pandemic influenza A(H1N1) v virus from recent human influenza viruses by real-time PCR

Detection and analysis of six lizard adenoviruses by consensus primer PCR provides further evidence of a reptilian origin for the atadenoviruses

Detection of new DNA polymerase genes of known and potentially novel herpesviruses by PCR with degenerate and deoxyinosine-substituted primers

Multicolour, multiplex real-time PCR assay for the detection of human-pathogenic poxviruses

GC content-based pan-pox universal PCR assays for poxvirus detection

Diseases in free-ranging bats from Germany

Hervey virus: Study on co-circulation with Henipaviruses in Pteropid bats within their distribution range from Australia to Africa

Isolation of multiple novel paramyxoviruses from pteropid bat urine

Velvet: Algorithms for de novo short read assembly using de Bruijn graphs

Ultrafast and memory-efficient alignment of short DNA sequences to the human genome

A novel method for fast and accurate multiple sequence alignment

jModelTest 2: More models, new heuristics and parallel computing

MRBAYES: Bayesian inference of phylogenetic trees

Open Access funding enabled and organized by Projekt DEAL.

The authors declare no competing interests.

The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-021-86435-4.Correspondence and requests for materials should be addressed to C.K.Reprints and permissions information is available at www.nature.com/reprints.Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.