Skip to main content
Download PDF

Genetic diversity and cross-species transmissibility of bat-associated picornaviruses from Spain

Virology Journal volume 21, Article number: 193 (2024) Cite this article

Abstract

Background

Emerging zoonotic diseases arise from cross-species transmission events between wild or domesticated animals and humans, with bats being one of the major reservoirs of zoonotic viruses. Viral metagenomics has led to the discovery of many viruses, but efforts have mainly been focused on some areas of the world and on certain viral families.

Methods

We set out to describe full-length genomes of new picorna-like viruses by collecting feces from hundreds of bats captured in different regions of Spain. Viral sequences were obtained by high-throughput Illumina sequencing and analyzed phylogenetically to classify them in the context of known viruses. Linear discriminant analysis (LDA) was performed to infer likely hosts based on genome composition.

Results

We found five complete or nearly complete genomes belonging to the family Picornaviridae, including a new species of the subfamily Ensavirinae. LDA suggested that these were true vertebrate viruses, rather than viruses from the bat diet. Some of these viruses were related to picornaviruses previously found in other bat species from distant geographical regions. We also found a calhevirus genome that most likely belongs to a proposed new family within the order Picornavirales, and for which genome composition analysis suggested a plant host.

Conclusions

Our findings describe new picorna-like viral species and variants circulating in the Iberian Peninsula, illustrate the wide geographical distribution and interspecies transmissibility of picornaviruses, and suggest new hosts for calheviruses.

Background

Emerging zoonotic viruses originate from cross-species transmission events involving wild or domesticated animals, particularly mammals, bats being one of the main natural reservoirs of viruses. Several human-infective coronaviruses, filoviruses, paramyxoviruses, lyssaviruses, bunyaviruses, and hantaviruses have a recent bat origin [1,2,3,4,5], and bats might be the source of older cross-species transmissions that led to the evolution of measles, mumps, and hepatitis C [6, 7]. Moreover, bats have transmitted viruses of high concern to domestic animals, such as swine acute diarrhea syndrome coronavirus [8]. This prominent role may simply reflect the fact that bats are one of the most abundant and diverse family of mammals, with more than 1200 species described [9, 10]. However, it has also been suggested that bats tend to carry more viruses than other mammals due to their highly gregarious and migratory lifestyles, or particularities in their immune system favoring the establishment of persistent infections [11].

The characterization of viruses circulating in wildlife is a daunting task, and although we only know a very small fraction of the virosphere, significant progresses have been made recently [12,13,14]. Virus discovery has been massively enhanced by metagenomics, which has become the standard approach to describe viral diversity in nature, as well as a powerful surveillance tool [15, 16]. The role of bats as likely reservoirs for emerging human coronavirus diseases, particularly COVID-19, has stimulated research focusing on this family of viruses. However, bats are also frequent carriers of other types of viruses that have been less intensively studied, including retroviruses, picornaviruses, herpesviruses and astroviruses [17]. Reflecting this research bias, coronaviruses account for about 50% of the total entries in the Bat-associated virus database [18]. Another important bias in the field is due to differences in research effort around the world, as almost half of the known bat viruses have been reported in China and other parts of Asia.

The Iberian Peninsula is home to more than 30 species of native bats [19], which have been reported to carry lyssaviruses [20], coronaviruses [21], and filoviruses [22], among other types of viruses. Although the Iberian Peninsula and the rest of Europe do not stand out as a global zoonotic hotspot [23], certain risk factors, such as frequent contact between humans and infected animals in disturbed habitats, are likely to occur in this region. In particular, extensive land-use changes and habitat alterations in the Iberian Peninsula may cause disturbances in the distribution of bat colonies and bat-associated viromes, potentially leading to increased spillover risks [24]. In this study, we applied viral metagenomics to bat feces sampled in different regions of Spain and found a relatively high frequency of picornavirus-like sequences. The order Picornavirales comprises nine families, including the family Picornaviridae, which consists of at least 68 genera containing 158 species, although new groups have been proposed and are currently awaiting approval by the International Committee for the Taxonomy of Viruses (ICTV). Several picornaviruses are pathogenic to humans and other mammals and vertebrates, causing diseases including meningitis, encephalitis, hepatitis, gastroenteritis, myocarditis, paralysis and respiratory infections [25]. Here, we describe five complete or nearly complete genomes belonging to bat-associated picornaviruses from different regions of Spain, as well as a calhevirus belonging to a proposed new family of picorna-like viruses [26].

Methods

Sample collection

Fecal samples were collected from 202 individuals belonging to 23 different bat species from May to October 2022 in different Spanish regions (Cantabria, Castelló, Lugo, Murcia, Salamanca, Teruel, and València; Fig. 1). Bats were captured using nylon mist nets and a harp trap (Austbat), and individuals were briefly placed in cotton bags to recover fresh fecal samples. Then, samples were collected separately for each individual, placed in tubes containing 500 µL of 1X phosphate-buffered saline (PBS), kept cold initially, and then at − 20 °C until they were transported to the laboratory and stored at -80 °C.

Fig. 1
figure 1

Bat captures for virome analysis. Names of the localities where sampling was carried out and the number of fecal samples collected are indicated

Full size image

Sample processing and RNA extraction

A fraction of the collected samples was combined into a total of 26 pools, each including between one and 15 samples from the same bat species (Additional file 1: Table S1). Before starting sample processing, each pool was spiked with 105 plaque-forming units (PFU) of vesicular stomatitis virus (VSV) to assess final viral recovery efficiency. Then, each pool was homogenized in a Precellys Evolution tissue homogenizer (Bertin) using 2 mL tubes and 1.4 mm ceramic beads. The homogenization protocol consisted of three cycles of 30 s each at 6500 rpm, with a 10 s pause between cycles. After homogenization, samples were centrifuged twice at 20,000 g for 3 min at 4 °C. Supernatants were then filtered using Minisart cellulose acetate syringe filters with a 1.2 µm pore size (Sartorius) and transferred to ultra-clean 2 mL tubes. For RNA extraction, 250 µL of the filtrate were initially cleaned with Trizol LS reagent (Invitrogen) and then 280 µL from this first step were used for final extraction using the QIAamp Viral RNA minikit (Qiagen). RNA was eluted in a final volume of 40 µL and stored at − 80 °C.

Sequencing and classification of viral sequences

Extracted RNA was used for library preparation using the stranded mRNA preparation kit (Illumina), skipping mRNA enrichment steps, and starting at the fragmentation step. Samples were subjected to paired-end sequencing using a NextSeq 550 device with read length of 150 bp at each end. Raw reads were de-duplicated and quality filtered using fastp v0.23.2 [27]. A quality trimming threshold of 20 was used, and reads below 70 nt in length were removed from the data set using fastp v0.23.2. De novo sequence assembly was performed using SPAdes v3.15.4 [28] with the meta option, and MEGAHIT v1.2.9 [29] with default parameters. Then, contigs assembled with both methods were clustered to remove replicates or small replicates of larger contigs using CD-HIT v4.8.1 [30], and those shorter than 1000 nt were removed. The obtained clustered sequences were taxonomically classified using Kaiju v1.9.0 [31] with the subset of NCBI nr protein database downloaded on June 6, 2023. Clustered sequences were also analyzed using Virsorter2 v2.2.4 [32] to detect viral contigs. Then, viral contigs were analyzed with CheckV v1.0.1 [33] using the CheckV database v1.5 to further assess their quality. Finally, contigs assigned to the Picornaviridae family or associated with picornaviruses were selected.

Genome annotation and phylogenetic analysis

Sequences similar to each contig of interest were searched using DIAMOND v2.0.15.153 [34] with the blastx option and the NCBI nr database downloaded on June 7, 2023. Then, the 100 sequences closest to each contig were extracted and analyzed for association with vertebrate-infecting viruses, while protein domains were annotated using Interproscan v5.63-95.0 [35] using the Pfam database v35.0. Open reading frames (ORFs) were predicted using ORFfinder. Nucleotide-based multiple alignments were performed using MAFFT v7.490 [36] and, where indicated, pairwise identities were calculated using the Sequence Demarcation Toolkit (SDT) v1.2 [37]. Coverage statistics of selected viral contigs were calculated by remapping the trimmed and filtered reads to each contig using Bowtie2 v2.2.5 [38]. In addition, selected viral contigs were compared to NCBI databases using BLAST to obtain identity values and refine annotations. For phylogenetic analysis of P1 and 3Dpol regions, nucleotide sequences were aligned with Clustal Omega v1.4.1 [39] and trimmed with ClipKIT [40]. The best substitution model was inferred using the ModelFinder feature in IQTREE [41], which was the GTR + I + G4 model for the P1 and 3Dpol Picornaviridae family phylogenies, GTR + I + G4 + F for the whole polyprotein Picornaviridae family phylogeny, and Rtrev + G4 for the Picornavirales phylogeny. Phylogenetic tree inference was performed using Mr.Bayes v3.2.7 [42] or IQ-TREE (v2.0.3) [43]. Branch support was estimated with 10,000 ultra-fast bootstrapping replicates and 1,000 bootstrap replicates for the SH-like approximate likelihood ratio test.

PCR confirmation

The origin of the new viruses detected using metagenomics was determined by individual PCR analysis of the different samples included in each pool. For this purpose, RNA isolation was independently performed for each sample with the QIAamp Viral RNA mini kit (Qiagen) following manufacturer's instructions. RNA was eluted in a volume of 30 µL, and 4 µL were used to obtain first strand cDNA with Superscript IV reverse transcriptase (Invitrogen) using specific reverse primers for each virus. Then, amplification was performed using NZYTaq II Green Master Mix (NZYTech) and specific primers (Additional file 3: Table S3). PCR products were visualized by electrophoresis on a 1% agarose gel using Green Safe Premium (NZYTech).

Nucleotide composition analysis

Nucleotide composition analysis was performed as previously described [26, 44]. Briefly, picorna-like genome sequences were downloaded from Virus-Host DB [45] on December, 2023 (www.genome.jp/virushostdb). Viruses defined as infecting more than one host type (i.e. mammals, insects, or plants) were removed and those with segmented genomes were concatenated, resulting in a total of 1065 genomes. Mononucleotide frequencies, dinucleotide frequencies and dinucleotide bias were calculated using the Biostrings package v.2.70.1 in R v4.3.2. These three sets of metrics were used to predict viral host through linear discriminant analysis using the R MASS package. Finally, the predictive power of the model was tested by leave-one-out cross-validation. Accession numbers, reported host, mono- and dinucleotide frequencies, dinucleotide bias, and resulting linear discriminants are reported in Additional file 4: Table S4.

Results and discussion

Overview

We obtained droppings from 202 animals belonging to 23 different bat species and captured in seven different regions of Spain between May and October 2022 (Fig. 1; Additional file 1: Table S1). Fecal material from individuals of the same species was pooled before sample processing. After RNA extraction and library preparation, Illumina sequencing generated between 4.5 and 48 million raw reads per pool (Additional file 2: Table S2).

Quality-filtered reads were de novo assembled and the resulting contigs were analyzed to detect viral sequences. In total, 430 high-quality complete or nearly complete viral genomes were obtained, the vast majority of which corresponded to bacteriophages. Of these, six corresponded to picorna-like sequences of at least 7kb, representing complete or quasi-complete viral genomes. These included five picornaviruses genomes, while the sixth was a new calhevirus. Proposed names for these viruses are provided in Table 1. Since the pools analysed contained samples from different regions, RNA extraction and RT-PCR were performed for each sample (Additional file 3: Table S3) in order to assign the geographical origin of described picorna-like genomes.

Table 1 Summary of picorna-like viral genomes found in bat feces from Spain
Full size table

We then annotated and characterized the genome structure of each picornavirus (Fig. 2), analysed characteristic sequence motifs (Table 2), and performed a Bayesian phylogenetic analysis to define their relationship to those of ICTV-approved viral species, using both the P1 and 3Dpol sequences for the Picornaviridae members (Fig. 3) and the RdRp sequence for the calhevirus (Fig. 4). However, since most viral genomes have not yet been officially classified, we also looked for the most similar sequences in Genbank to ascertain whether our newly described viruses constituted novel species according to ICTV picornavirus species demarcation criteria [46]. Finally, in an attempt to rule out a diet or environmental origin, we used linear discriminant analysis (LDA) to infer the most likely host of each virus based on nucleotide composition and dinucleotide biases (Fig. 5, Additional file 4: Table S4). Using 1065 picorna-like viruses (794 mammalian, 34 fish, 59 avian, 99 plant, and 79 insect viruses), leave-one-out cross-validation indicated that vertebrate, plant, and insect viruses can be confidently distinguished, whereas the separation between mammalian, avian and fish viruses is less clear (Fig. 5A). This model was then applied to identify the most likely host of the newly described bat-associated viruses. For the sake of clarity, the following sections will detail the above analyses individually for each new virus described.

Fig. 2
figure 2

Genome organization of viruses identified in this study belonging to Picornaviridae family. Open reading frames for each mature protein (L, VP1-4, 2A-C, 3A-D) are shown as colored regions and mapped to P1, P2, or P3 polyproteins accordingly

Full size image
Table 2 Picornavirus conserved sequence motifs in newly reported viruses
Full size table
Fig. 3
figure 3

Phylogenetic trees of the Picornaviridae family for the P1 and 3Dpol regions. The newly described viruses are placed in the context of ICTV-approved viruses. Ampivirus A1 (accession KP770140) was selected as an outgroup for both trees. Phylogenetic inference was performed with the Bayesian MCMC algorithm using a GTR + I + G4 substitution model, and 12 and 8 million generations for P1 and 3Dpol phylogenies, respectively. Posterior probabilities for each node are indicated (values of 1 are omitted)

Full size image
Fig. 4
figure 4

Phylogenetic analysis of the RdRp region for Caliciviridae, Dicistroviridae, Iflaviridae, Marnaviridae, Noraviridae, Picornaviridae, and Secoviridae families, previously known calheviruses, and Castelló calhevirus. Phylogenetic inference was performed with the Bayesian MCMC algorithm using 10 million generations and a Rtrev + G4 substitution model. Posterior probabilities for each node are indicated (values of 1 are omitted)

Full size image
Fig. 5
figure 5

Linear discriminant analysis for the prediction of picornaviruses hosts. A. Confusion matrix obtained with the training dataset of 1065 viral genomes. B. Posterior probabilities for the classification of novel picorna-like genomes into each class (values below 0.01 are omitted). C. Biplots of linear discriminant analysis. Novel picorna-like sequences described in this study are shown with black dots

Full size image

New member of the subfamily Ensavirinae in Barbastella bats

We obtained a 7.0 kb picornavirus sequence from the feces of Barbastella barbastellus animals captured in Rábade and Begonte, two localities in the province of Lugo. The genome organization of this viral genome was typical of the subfamily Ensavirinae, a classification that was confirmed by phylogenetic analysis of the P1, 3Dpol, and full-length genome sequences (Additional file 5: Figure S1). The new sequence exhibited 63.5%, 66.4%, and 74.5% amino acid identity to Bat sapelovirus strain Henan-Mr221 (GenBank accession MF352431) in polyprotein, P1, and 2C + 3CD sequences, respectively. This situates the newly described genome close to the boundary for defining a new species according to ICTV criteria for sapeloviruses (less than 70%, 64%, and 70% identity in the polyprotein, P1 and 2C + 3CD amino acid sequences). While the P1 region was more closely related to that of sapelovirus and diresapivirus, the 3Dpol and full genome phylogenies placed it closer to parabovirus and felipivirus. We hence tentatively called this virus Lugo bat picornavirus. The cross-validated LDA model assigned a 93% probability for the virus to infect mammals, hence providing confidence that this was a true bat picornavirus and not a virus present in the diet or the environment.

Variants of the Scotophilus kuhlii picornavirus in different Myotis bats

Two similar sequences of 7.7 kb sharing > 94% amino acid identity in the P1 and 2C + 3CD regions were obtained from a Myotis myotis individual captured in Navarredonda de la Rinconada (province of Salamanca) and a M. blythii animal dwelling in a cave called Cueva de la Almagra (Murcia) 460 km away from the other sample. Their genome structure was again typical of the subfamily Ensavirinae, as confirmed by phylogenetic analysis of the P1 and 3Dpol sequences. Phylogenetic reconstruction using full-length genomes also placed both sequences into the Ensavirinae subfamily (Additional file 5: Figure S1). Both sequences exhibited canonical motifs in 2C and 3D proteins, but showed an alanine (GxCGx10AxHxxG) residue instead of the more widely conserved glycine residue in the 3C motif (GxCGx10GxHxxG). The alanine in the active site of the 3C proteinase is present in other members of the Ensavirinae subfamily, including all ICTV accepted species of the Parabovirus, Felipivirus and Diresapivirus genera. According to our phylogenetic analysis, the P1 sequences of both viruses were more closely related to those of felipivirus, boosepivirus, and parabovirus, while 3Dpol sequences fell into the genus Boosepivirus. We thus conservatively called these viruses Navarredonda bat picornavirus and Almagra bat picornavirus. The LDA model assigned a 100% probability of Navarredonda bat picornavirus to infect mammals, whereas this probability was 82% for Almagra bat coronavirus. The closest Genbank match in a blastp analysis, with > 90% amino acid identity in the P1 and 2C + 3CD, corresponded to Scotophilus kuhlii picornavirus, which was discovered in bat feces in 2014 in Vietnam (Genbank accession MZ544244, unpublished study). Therefore, according to ICTV guidelines, these viruses should be defined as variants of the same species. However, we notice that they were found in bats from different genera (Scotophilus and two Myotis species) and in very distant regions, suggesting a global distribution and broad host range for this virus.

A variant of mischivirus B1 found in the feces of a Miniopterus bat from Murcia

A picornavirus sequence of 8.5 kb was detected in two Miniopterus scheibersii bats from an arid location called Cabezo Gordo, in the province of Murcia. This sequence displayed typical conserved motifs in 2C, 3C and 3D proteins, and was found to possess a NPG|P motif in the 2A polypeptide, which is present in the Mischivirus genus (Additional file 5: Figure S1) among other genera. Phylogenetic analysis of the P1 and 3Dpol sequences confirmed that this virus belongs to the Mischivirus genus, with higher posterior probabilities in the 3Dpol phylogeny. Phylogenetic analysis using full-length genomes also placed this sequence as a member of the Mischivirus genus (Additional file 5: Figure S1). A blastp analysis of entire polyprotein revealed a 96.2% identity to mischivirus B1, a virus found in the feces of M. schreibersii in Hungary [47], indicating that our virus constitutes a variant of mischivirus B1. Genome composition analysis identified birds as the most likely host of the newly described virus with an 80% probability, versus only a 18% for mammals. Thus, an environmental origin from bird material cannot be ruled out, although it should be noted that the LDA model showed poor ability to discern between different types of vertebrate hosts. We hence called this virus Cabezo Gordo bat-associated mischivirus. In any case, it is highly unlikely that this virus came from the diet since M. scheibersii is an insectivorous bat.

A divergent variant of Wenzhou bat picornavirus 7 found in Pipistrellus bats from Valencia

Another picornavirus sequence of 8.5 kb was found in the feces of four Pipistrellus kuhlii individuals from different pools obtained in Fontanars, in the province of València. The sequence exhibited a genome layout matching the Kobuvirus genus, and protein 2A displayed an H-box (HWAV) upstream from a NC-box that is also typical of this genus [48]. Phylogenetic analysis with P1, 3Dpol, and full-length genome sequences (Additional file 5: Figure S1) confirmed that this was a member of the Kobuvirus genus. Genome composition analysis identified birds as the most likely hosts with an 88% probability, versus only a 12% for mammals. Thus, we tentatively called this virus Fontanars bat-associated kobuvirus. As above, it is unlikely that this virus came from the diet, since P. kuhlii is also an insectivorous bat. An environmental origin from bird material cannot be discarded but appears to be less likely than an incorrect classification given the poor ability of the LDA model to discern between avian and mammalian hosts. The closest ICTV-accepted species in the phylogeny was aichivirus F1, but with a high divergence. The best blastp hit corresponded to Wenzhou bat picornavirus 7 (Genbank accession OQ716048), a picornavirus recently discovered in P. abramus from China [49]. These two sequences shared 72.8% sequence identity across the whole polyprotein, while percentage identities for P1 and 2C + 3CD were 77.6% and 74.7%, respectively. Hence, according to ICTV criteria, our sequence and Wenzhou bat picornavirus 7 should be considered members of the same species. However, they are relatively divergent and were found in different Pipistrellus species from very distant regions, suggesting that this virus is present globally, displays high genetic diversity, and shows interspecies transmissibility.

A plant-like calhevirus found in a Tadarida bat from Castelló

A 10 kb picorna-like sequence was obtained from the feces of Tadarida teniotis bats captured in Castelló. This sequence exhibited quite dissimilar motifs in the 2C and 3C proteins compared to picornaviruses, showing a GxxRxGKT motif in the 2C protein instead of the widely conserved GxxGxGKS, and a serine residue in the catalytic site of the 3C protease (GxSGx10GxYxxA) instead of the typical cysteine. The best hit in a blastn analysis using the entire sequence (68.7% nucleotide identity along a 28% coverage span) corresponded to Jinngmen bat picorna-like virus 1, a virus recently discovered in Myotis chinensis [49]. All other top hits belonged to calhevirus sequences [26, 50], with similar nucleotide identities and slightly lower coverages. The non-structural ORF1 of our sequence shared a 51.0% amino acid identity to that of Jingmen bat picorna-like virus 1, whereas the capsid polyproteins exhibited 53.6%, indicating that this should be a novel calhevirus species, which we called Castelló calhevirus. Genome composition analysis identified plants as the most likely hosts with 98% probability. T. teniotis is yet another insectivorous bat, but ingestion of plant material associated with insects cannot be discarded.

Conclusions

Investigation of the virome of 23 different bat species across Spain’s geography led to the identification of six complete or quasi-complete new viral genomes belonging to the Picornavirales order. Our results provide information about the host, geographical range, and genetic diversity of bat-associated picornaviruses. Two of the described viral genomes may constitute new species according to ICTV criteria. One of these viruses, Lugo bat picornavirus, is likely to infect bats, whereas the other is a divergent calhevirus that may infect plants, according to genome composition analysis. Previously annotated calheviruses were isolated from human, shrew, and rat feces, and were inferred to infect arthropod hosts using genome composition analysis [26, 50, 51]. The newly described Castelló calhevirus might constitute the first member of this family infecting plants.

We also found a new kobuvirus in multiple Pipistrellus kuhlii individuals, suggesting that this was a prevalent virus amongst the local bat population of València, although this abundance might be related to a certain outbreak in these bats. Importantly, members of the Kobuvirus genus can cause gastroenteritis and diarrhea in humans and other mammals such as pigs [52, 53]. The newly detected kobovirus was highly divergent from its closest ICTV-accepted relative (Aichivirus F2), and was closer to a virus recently discovered in a different species of the same bat genus in China, the similarity between these two viruses being close to the species demarcation threshold (ca. 70% amino acid identity).

Three additional viral genomes belonged more clearly to previously described species. Two of these genomes, obtained from different species of Myotis in different regions of Spain, showed > 90% amino acid identity with a picornavirus described in 2014 in lesser Asiatic yellow bats (Scotophilus kuhlii) from Vietnam, showing that these sequences pertain to a geographically widespread and genetically diverse picornavirus capable of infecting bats from different genera. Another picornavirus genome found in the feces of M. schreibersii was shown to be a variant of mischivirus B1 with 96.2% amino acid identity to another viral sequence obtained from the same bat species in Hungary in 2013 [47]. Mischiviruses have also been reported in canine and shrew-derived samples, suggesting that this genus can infect a broad range of mammals.

Both the Fontanars bat-associated kobuvirus and the Cabezo Gordo bat-associated mischivirus most probable hosts were predicted to be birds, rather than mammals, despite being isolated from bat fecal samples. Although dietary or environmental origins are unlikely considering their insectivorous nature, such origins cannot be fully disregarded. For this purpose, further studies into the virome of both Spanish bat and bird communities could assist in unraveling the tropism range of these newly identified viruses. In that direction, further research effort is warranted to confirm the predicted bird tropism, which could increase their zoonotic potential, also considering the expanded cross-species infectivity and transmissibility between bats species reported here for the Fontanars-bat associated kobuvirus, and Navarredonda and Almagra bat picornaviruses (Additional file 5).

Overall, our results show that picornaviruses are distributed widely in the Spanish bat population. For some of these picornaviruses, we found evidence that variants of the same virus can be found in distant parts of the world and can infect divergent bat species and, potentially, other vertebrates. In light of these findings, special protection of the natural habitats of Spanish bats is warranted to minimize zoonosis risk. Further research efforts should be implemented towards characterizing the abundance of picornaviruses and other viral families in different types of hosts and different regions of the world in order to obtain a more complete picture of their diversity, cross-species transmissibility, and associated human spillover risks.

Availability of data and materials

Raw sequence reads were deposited in the GenBank Sequence Read Archive under accession number SRR27912329, SRR27912332, SRR27912333, SRR27912339, SRR27912350, SRR27912351, and SRR27962820-SRR27962839. The viral contigs corresponding to complete or nearly complete genomes were deposited in Genbank under the accession numbers PP654851-PP654856.

References

  1. Geng R, Zhou P. Severe acute respiratory syndrome (SARS) related coronavirus in bats. Anim Dis. 2021;1:4.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Kohl C, Kurth A. European bats as carriers of viruses with zoonotic potential. Viruses. 2014;6:3110–28.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Shipley R, Wright E, Selden D, Wu G, Aegerter J, Fooks AR, et al. Bats and viruses: emergence of novel lyssaviruses and association of bats with viral zoonoses in the EU. Trop Med Infect Dis. 2019;4:31.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Simons RRL, Gale P, Horigan V, Snary EL, Breed AC. Potential for introduction of bat-borne zoonotic viruses into the EU: a review. Viruses. 2014;6:2084–121.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Wang L-F, Anderson DE. Viruses in bats and potential spillover to animals and humans. Curr Opin Virol. 2019;34:79–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Drexler JF, Corman VM, Müller MA, Maganga GD, Vallo P, Binger T, et al. Bats host major mammalian paramyxoviruses. Nat Commun. 2012;3:796.

    Article  PubMed  Google Scholar 

  7. Quan P-L, Firth C, Conte JM, Williams SH, Zambrana-Torrelio CM, Anthony SJ, et al. Bats are a major natural reservoir for hepaciviruses and pegiviruses. Proc Natl Acad Sci U S A. 2013;110:8194–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Zhou P, Fan H, Lan T, Yang X-L, Shi W-F, Zhang W, et al. Fatal swine acute diarrhoea syndrome caused by an HKU2-related coronavirus of bat origin. Nature. 2018;556:255–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Mollentze N, Streicker DG. Viral zoonotic risk is homogenous among taxonomic orders of mammalian and avian reservoir hosts. Proc Natl Acad Sci USA. 2020;117:9423–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Tsagkogeorga G, Parker J, Stupka E, Cotton JA, Rossiter SJ. Phylogenomic analyses elucidate the evolutionary relationships of bats. Curr Biol. 2013;23:2262–7.

    Article  CAS  PubMed  Google Scholar 

  11. Brook CE, Dobson AP. Bats as “special” reservoirs for emerging zoonotic pathogens. Trends Microbiol. 2015;23:172–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Holmes EC. The ecology of viral emergence. Annu Rev Virol. 2022;9:173–92.

    Article  CAS  PubMed  Google Scholar 

  13. Carlson CJ, Zipfel CM, Garnier R, Bansal S. Global estimates of mammalian viral diversity accounting for host sharing. Nat Ecol Evol. 2019;3:1070–5.

    Article  PubMed  Google Scholar 

  14. Carroll D, Daszak P, Wolfe ND, Gao GF, Morel CM, Morzaria S, et al. The global virome project. Science. 2018;359:872–4.

    Article  CAS  PubMed  Google Scholar 

  15. Bassi C, Guerriero P, Pierantoni M, Callegari E, Sabbioni S. Novel virus identification through metagenomics: a systematic review. Life (Basel). 2022;12:2048.

    CAS  PubMed  Google Scholar 

  16. Roev GV, Borisova NI, Chistyakova NV, Agletdinov MR, Akimkin VG, Khafizov K. Unlocking the viral universe: metagenomic analysis of bat samples using next-generation sequencing. Microorganisms. 2023;11:2532.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Cui X, Fan K, Liang X, Gong W, Chen W, He B, et al. Virus diversity, wildlife-domestic animal circulation and potential zoonotic viruses of small mammals, pangolins and zoo animals. Nat Commun. 2023;14:2488.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Database of Bat-Associated Viruses (DBatVir), accessed on July 23, 2020. [Internet]. 2020. Available from: www.mgc.ac.cn/cgi-bin/DBatVir/main.cgi

  19. Moraga-Fernández A, Sánchez-Sánchez M, Queirós J, Lopes AM, Vicente J, Pardavila X, et al. A study of viral pathogens in bat species in the Iberian Peninsula: identification of new coronavirus genetic variants. Int J Vet Sci Med. 2022;10:100–10.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Aréchiga Ceballos N, Vázquez Morón S, Berciano JM, Nicolás O, Aznar López C, Juste J, et al. Novel lyssavirus in bat, Spain. Emerg Infect Dis. 2013;19:793–5.

    Article  PubMed  Google Scholar 

  21. Falcón A, Vázquez-Morón S, Casas I, Aznar C, Ruiz G, Pozo F, et al. Detection of alpha and betacoronaviruses in multiple Iberian bat species. Arch Virol. 2011;156:1883–90.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Negredo A, Palacios G, Vázquez-Morón S, González F, Dopazo H, Molero F, et al. Discovery of an ebolavirus-like filovirus in europe. PLoS Pathog. 2011;7:e1002304.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Allen T, Murray KA, Zambrana-Torrelio C, Morse SS, Rondinini C, Di Marco M, et al. Global hotspots and correlates of emerging zoonotic diseases. Nat Commun. 2017;8:1124.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Serra-Cobo J, López-Roig M. Bats and emerging infections: an ecological and virological puzzle. Adv Exp Med Biol. 2017;972:35–48.

    Article  PubMed  Google Scholar 

  25. Andino R, Kirkegaard K, Macadam A, Racaniello VR, Rosenfeld AB. The picornaviridae family: knowledge gaps, animal models, countermeasures, and prototype pathogens. J Infect Dis. 2023;228:S427–45.

    Article  PubMed  Google Scholar 

  26. Kapoor A, Simmonds P, Lipkin WI, Zaidi S, Delwart E. Use of nucleotide composition analysis to infer hosts for three novel picorna-like viruses. J Virol. 2010;84:10322–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34:i884–90.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 2017;27:824–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Li D, Liu C-M, Luo R, Sadakane K, Lam T-W. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31:1674–6.

    Article  CAS  PubMed  Google Scholar 

  30. Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22:1658–9.

    Article  CAS  PubMed  Google Scholar 

  31. Menzel P, Ng KL, Krogh A. Fast and sensitive taxonomic classification for metagenomics with Kaiju. Nat Commun. 2016;7:11257.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Guo J, Bolduc B, Zayed AA, Varsani A, Dominguez-Huerta G, Delmont TO, et al. VirSorter2: a multi-classifier, expert-guided approach to detect diverse DNA and RNA viruses. Microbiome. 2021;9:37.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Nayfach S, Camargo AP, Schulz F, Eloe-Fadrosh E, Roux S, Kyrpides NC. CheckV assesses the quality and completeness of metagenome-assembled viral genomes. Nat Biotechnol. 2021;39:578–85.

    Article  CAS  PubMed  Google Scholar 

  34. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12:59–60.

    Article  CAS  PubMed  Google Scholar 

  35. Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Muhire BM, Varsani A, Martin DP. SDT: a virus classification tool based on pairwise sequence alignment and identity calculation. PLoS One. 2014;9:e108277.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7:539.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Steenwyk JL, Buida TJ, Li Y, Shen X-X, Rokas A. ClipKIT: a multiple sequence alignment trimming software for accurate phylogenomic inference. PLoS Biol. 2020;18:e3001007.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Nguyen L-T, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32:268–74.

    Article  CAS  PubMed  Google Scholar 

  42. Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.

    Article  CAS  PubMed  Google Scholar 

  43. Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Teeling E, editor. Mol Biol Evol. 2020;37:1530–4.

  44. Yinda CK, Zell R, Deboutte W, Zeller M, Conceição-Neto N, Heylen E, et al. Highly diverse population of Picornaviridae and other members of the Picornavirales, in Cameroonian fruit bats. BMC Genomics. 2017;18:249.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Virus-Host Database [Internet]. [cited 2024 Apr 4]. Available from: https://www.genome.jp/virushostdb/

  46. Zell R, Delwart E, Gorbalenya AE, Hovi T, King AMQ, Knowles NJ, et al. ICTV virus taxonomy profile: picornaviridae. J Gen Virol. 2017;98:2421–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Kemenesi G, Zhang D, Marton S, Dallos B, Görföl T, Estók P, et al. Genetic characterization of a novel picornavirus detected in Miniopterus schreibersii bats. J Gen Virol. 2015;96:815–21.

    Article  CAS  PubMed  Google Scholar 

  48. Yang X, Cheng A, Wang M, Jia R, Sun K, Pan K, et al. Structures and corresponding functions of five types of picornaviral 2A proteins. Front Microbiol. 2017;8:1373.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Chen Y-M, Hu S-J, Lin X-D, Tian J-H, Lv J-X, Wang M-R, et al. Host traits shape virome composition and virus transmission in wild small mammals. Cell. 2023;186:4662-4675.e12.

    Article  CAS  PubMed  Google Scholar 

  50. Sasaki M, Orba Y, Ueno K, Ishii A, Moonga L, Hang’ombe BM, et al. Metagenomic analysis of the shrew enteric virome reveals novel viruses related to human stool-associated viruses. J Gen Virol. 2015;96:440–52.

    Article  CAS  PubMed  Google Scholar 

  51. Siqueira JD, Dominguez-Bello MG, Contreras M, Lander O, Caballero-Arias H, Xutao D, et al. Complex virome in feces from Amerindian children in isolated Amazonian villages. Nat Commun. 2018;9:4270.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Zhang Q, Hu R, Tang X, Wu C, He Q, Zhao Z, et al. Occurrence and investigation of enteric viral infections in pigs with diarrhea in China. Arch Virol. 2013;158:1631–6.

    Article  CAS  PubMed  Google Scholar 

  53. Bergallo M, Galliano I, Montanari P, Rassu M, Daprà V. Aichivirus in children with diarrhea in Northern Italy. Intervirology. 2017;60:196–200.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank members of the Virus Evolution lab for helpful comments and suggestions.

Funding

This research was financially supported by grant PID2020-118602RB-I00 from the Spanish Ministerio de Ciencia e Innovación (MICINN) and grant CIAICO/2022/110 from the Conselleria de Educación, Universidades y Empleo (Generalitat Valenciana). MCS was funded by a Ph.D. fellowship from the Spanish Ministerio de Ciencia e Innovación, the Agencia Estatal de Investigación and the European Social Fund Plus (PRE2021‐099824).

Author information

Authors and Affiliations

  1. Institute for Integrative Systems Biology (I2SysBio), Universitat de València – CSIC, València, Spain

    Marc Carrascosa-Sàez, Jaime Buigues, Iván Andreu-Moreno, Raquel Martínez-Recio, Clàudia Soriano-Tordera, José M. Cuevas & Rafael Sanjuán

  2. Institut Cavanilles de Biodiversitat I Biologia Evolutiva, Universitat de València, València, Spain

    Adrià Viñals & Juan S. Monrós

  3. Department of Genetics, Universitat de València, València, Spain

    José M. Cuevas & Rafael Sanjuán

Authors
  1. Marc Carrascosa-Sàez
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Jaime Buigues
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Adrià Viñals
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Iván Andreu-Moreno
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Raquel Martínez-Recio
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Clàudia Soriano-Tordera
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Juan S. Monrós
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. José M. Cuevas
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Rafael Sanjuán
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

MCS integrated data analyses and drafted the manuscript. JB obtained viral contigs and analyzed data. AV obtained bat feces samples. IAM performed LDA. RMR prepared samples for sequencing. CST contributed to viral detection by PCR. JSM supervised work and obtained bat feces samples. JMC obtained funding, supervised work, and revised the manuscript. RS obtained funding, supervised work, and co-wrote the manuscript.

Corresponding authors

Correspondence to José M. Cuevas or Rafael Sanjuán.

Ethics declarations

Ethics approval and consent to participate

Samples consisted of feces from wild animals captured using nylon mist nets or a harp trap. Bats were kept briefly in cotton bags until fresh fecal samples were obtained. According to the European Directive regulating the protection of animals used for scientific purposes (2010/62/EU, Article 1), subsequently transposed into Spanish legislation (Royal Decree 53/213, 1 February, Article 2), procedures used in this study are not subject to the condition of animal experimentation and therefore do not require a IACUC approval document, but specifically a permit from the competent regional authority. The necessary permits from the Generalitat Valenciana for the sampling of wild bats were granted under Exp. 2022-VS (FAU22_009).

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

. Table S1. Detailed list of specimens and samples examined in this study. Detailed list of sampled specimens, including collection date, species name, age, sex, type of trap used and pooling

Additional file 2

. Table S2. Illumina reads and viral contigs obtained. Sequencing data for each pool, including raw and clean reads, number of contigs, number of viral contigs, and complete or high-quality RNA genomes

Additional file 3

. Table S3. Specific primers used for PCR confirmation. Target viral genome and target genomic region, 5’–3’ sequence, and amplicon length are indicated

Additional file 4

. Table S4. Nucleotide composition-based linear discriminant model. Detailed information for the model used, including accession numbers, reported host, mono- and dinucleotide frequencies, dinucleotide bias, and resulting linear discriminants

Additional file 5

. Figure S1. Phylogenetic tree of the Picornaviridae family for the complete polyprotein. The newly described viruses are placed in the context of ICTV-approved viruses. Ampivirus A1 (accession KP770140) was used as outgroup. Phylogenetic inference was done using maximum likelihood and a GTR+I+G4+F substitution model. Bootstrap values are shown in grey for each node

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Carrascosa-Sàez, M., Buigues, J., Viñals, A. et al. Genetic diversity and cross-species transmissibility of bat-associated picornaviruses from Spain. Virol J 21, 193 (2024). https://doi.org/10.1186/s12985-024-02456-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12985-024-02456-1

Keywords

Virology Journal

ISSN: 1743-422X