Identification of a novel bat papillomavirus by metagenomics

Full-genome sequencing of dozens of new DNA viruses found in Spanish bat feces

Bats are natural hosts of multiple viruses, many of which have clear zoonotic potential. The search for emerging viruses has been aided by the implementation of metagenomic tools, which have also enabled the detection of unprecedented viral diversity. Currently, this search is mainly focused on RNA viruses, which are largely over-represented in databases. To compensate for this research bias, we analyzed fecal samples from 189 Spanish bats belonging to 22 different species using viral metagenomics. This allowed us to identify 52 complete or near-complete viral genomes belonging to the families Adenoviridae, Circoviridae, Genomoviridae, Papillomaviridae, Parvoviridae, Polyomaviridae and Smacoviridae. Of these, 30 could constitute new species, doubling the number of viruses currently described in Europe. These findings open the door to a more thorough analysis of bat DNA viruses and their zoonotic potential.

IMPORTANCE

Metagenomics has become a fundamental tool to characterize the global virosphere, allowing us not only to understand the existing viral diversity and its ecological implications but also to identify new and emerging viruses. RNA viruses have a higher zoonotic potential, but this risk is also present for some DNA virus families. In our study, we analyzed the DNA fraction of fecal samples from 22 Spanish bat species, identifying 52 complete or near-complete genomes of different viral families with zoonotic potential. This doubles the number of genomes currently described in Europe. Metagenomic data often produce partial genomes that can be difficult to analyze. Our work, however, has characterized a large number of complete genomes, thus facilitating their taxonomic classification and enabling different analyses to be carried out to evaluate their zoonotic potential. For example, recombination studies are relevant since this phenomenon could play a major role in cross-species transmission.

Novel Virus Identification through Metagenomics: A Systematic Review

Metagenomic Next Generation Sequencing (mNGS) allows the evaluation of complex microbial communities, avoiding isolation and cultivation of each microbial species, and does not require prior knowledge of the microbial sequences present in the sample. Applications of mNGS include virome characterization, new virus discovery and full-length viral genome reconstruction, either from virus preparations enriched in culture or directly from clinical and environmental specimens. Here, we systematically reviewed studies that describe novel virus identification through mNGS from samples of different origin (plant, animal and environment). Without imposing time limits to the search, 379 publications were identified that met the search parameters. Sample types, geographical origin, enrichment and nucleic acid extraction methods, sequencing platforms, bioinformatic analytical steps and identified viral families were described. The review highlights mNGS as a feasible method for novel virus discovery from samples of different origins, describes which kind of heterogeneous experimental and analytical protocols are currently used and provides useful information such as the different commercial kits used for the purification of nucleic acids and bioinformatics analytical pipelines.

Vasomotor effects of noradrenaline, 5-hydroxytryptamine, angiotensin II, bradykinin, histamine, and acetylcholine on the bat (Rhinolophus ferrumequinum) basilar artery

The responsiveness of the basilar artery to intrinsic vasoactive substances is species-specific and can be a unique characteristic. We investigated the responsiveness of the bat (Rhinolophus ferrumequinum) basilar artery to noradrenaline (NA), 5-hydroxytryptamine (5-HT), angiotensin (Ang) II, bradykinin (BK), histamine (His), and acetylcholine (ACh). NA, 5-HT, Ang II, and BK induced contraction, whereas His and ACh induced relaxation, in a concentration-dependent manner. The NA cumulative concentration-response curve was shifted to the right in parallel with phentolamine (an α-antagonist). However, propranolol, a β-antagonist, had no significant effect. The 5-HT curve was shifted to the right in parallel by ketanserin (a 5-HT₂ antagonist) and methiothepin (a 5-HT₁ and 5-HT₂ antagonist). Losartan (an AT₁ antagonist) shifted the Ang II curve to the right, whereas PD123319 (an AT₂ antagonist) had no significant effect. L-NA, indomethacin, and des-Arg⁹-[Leu⁸]-BK (a B₁ antagonist) did not significantly affect BK-induced contractions. HOE140 (a B₂ antagonist) shifted the BK concentration-response curve to the right. The His curve was shifted to the right weakly by diphenhydramine (an H₁ antagonist) and strongly by cimetidine (a H₂ antagonist). ACh-induced relaxation was significantly inhibited by L-NA, atropine, and pFHHSiD (a muscarinic M₃ antagonist), whereas pirenzepine and methoctramine (muscarinic M₁ and M₂ antagonists, respectively) showed no significant effects. At a resting vascular tone, L-NA-induced contraction and indomethacin induced relaxation. These results suggest that α-adrenergic, 5-HT₁, 5-HT₂, AT₁, and B₂ receptors might be important in arterial contraction, whereas M₃ and H₂ (>H₁) receptors might modify these contractions, inducing relaxation.

Bat virome research: the past, the present and the future

Bats have been increasingly recognised as an exceptional reservoir for emerging zoonotic viruses for the past few decades. Recent studies indicate that the unique bat immune system may be partially responsible for their ability to co-exist with viruses with minimal or no clinical diseases. In this review, we discuss the history and importance of bat virome studies and contrast the vast difference between such studies before and after the introduction of next generation sequencing (NGS) in this area of research. We also discuss the role of discovery serology and high-throughput single cell RNA-seq in future bat virome research.

The virome of German bats: comparing virus discovery approaches

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.

The heart microbiome of insectivorous bats from Central and South Eastern Europe

Host associated microbiome not only may affect the individual health-status or provide insights into the species- or group specific bacterial communities but may act as early warning signs in the assessment of zoonotic reservoirs, offering clues to predict, prevent and control possible episodes of emerging zoonoses. Bats may be carriers and reservoirs of multiple pathogens such as viruses, bacteria and parasites, showing in the same time robust immunity against many of them. The microbiota plays a fundamental role on the induction, training and function of the host immune system and the immune system has largely evolved in order to maintain the symbiotic relationship of the host with these diverse microbes. Thus, expanding our knowledge on bat-associated microbiome it can be usefully in understanding bats' outstanding immune capacities. The aim of this study was to investigate the presence of different bacterial communities in heart tissue of insectivorous bats, Nyctalus noctula, Pipistrellus pipistrellus and Rhinoplophus hipposideros, from Central and Eastern Europe using high-throughput sequencing of variable regions of the 16S rRNA. In addition, species-specific PCRs were used to validate the presence of the vector-borne pathogens Bartonella spp. and Rickettsia spp. In this study we identified a wide variety of bacterial groups, with the most abundant phyla being Proteobacteria and Firmicutes. The results showed that at individual level, the year or location had no effect on the diversity and composition of the microbiome, however host species determined both structure and abundance of the bacterial community. We report the presence of vector-borne bacteria Bartonella spp. in samples of N. noctula and indications of Rickettsia spp. in R. hipposideros. Our results provide a first insight into the bacterial community found in heart tissue of bats from Central and South Eastern Europe.

Potential role of viral metagenomics as a surveillance tool for the early detection of emerging novel pathogens

Since the early times, human beings have always been faced with deadly microbial infections, both bacterial and viral. Pathogens such as viruses are always evolving owing to the processes of antigenic shift and drift. Such viral evolution results in the emergence of new types and serovars that prove deadly for humans-like influenza pandemics, severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS). The pandemic of novel coronavirus SARS-CoV-2 is the recent example. It has resulted in a great loss of human lives and a serious burden on economy across the globe. To counter such situations, a system should exist for the early detection of emerging viral pathogens. This will help prevent possible outbreaks and save human lives. Most of such deadly novel viruses and viral outbreaks are known to be originated from animal hosts. Regular monitoring of potential hot spots of such emerging microbes, such as zoos and animal markets, through metagenomics could help assess the presence of new viruses and pathogens. In this review, we focus on the potential of viral metagenomics and propose a surveillance system based on it for the early detection and hence prevention of such emerging viral infections.

A Preliminary Study of the Virome of the South American Free-Tailed Bats (Tadarida brasiliensis) and Identification of Two Novel Mammalian Viruses

Bats provide important ecosystem services as pollinators, seed dispersers, and/or insect controllers, but they have also been found harboring different viruses with zoonotic potential. Virome studies in bats distributed in Asia, Africa, Europe, and North America have increased dramatically over the past decade, whereas information on viruses infecting South American species is scarce. We explored the virome of Tadarida brasiliensis, an insectivorous New World bat species inhabiting a maternity colony in Rosario (Argentina), by a metagenomic approach. The analysis of five pooled oral/anal swab samples indicated the presence of 43 different taxonomic viral families infecting a wide range of hosts. By conventional nucleic acid detection techniques and/or bioinformatics approaches, the genomes of two novel viruses were completely covered clustering into the Papillomaviridae (Tadarida brasiliensis papillomavirus type 1, TbraPV1) and Genomoviridae (Tadarida brasiliensis gemykibivirus 1, TbGkyV1) families. TbraPV1 is the first papillomavirus type identified in this host and the prototype of a novel genus. TbGkyV1 is the first genomovirus reported in New World bats and constitutes a new species within the genus Gemykibivirus. Our findings extend the knowledge about oral/anal viromes of a South American bat species and contribute to understand the evolution and genetic diversity of the novel characterized viruses.

A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster

BACKGROUND

An ongoing outbreak of pneumonia associated with a novel coronavirus was reported in Wuhan city, Hubei province, China. Affected patients were geographically linked with a local wet market as a potential source. No data on person-to-person or nosocomial transmission have been published to date.

METHODS

In this study, we report the epidemiological, clinical, laboratory, radiological, and microbiological findings of five patients in a family cluster who presented with unexplained pneumonia after returning to Shenzhen, Guangdong province, China, after a visit to Wuhan, and an additional family member who did not travel to Wuhan. Phylogenetic analysis of genetic sequences from these patients were done.

FINDINGS

From Jan 10, 2020, we enrolled a family of six patients who travelled to Wuhan from Shenzhen between Dec 29, 2019 and Jan 4, 2020. Of six family members who travelled to Wuhan, five were identified as infected with the novel coronavirus. Additionally, one family member, who did not travel to Wuhan, became infected with the virus after several days of contact with four of the family members. None of the family members had contacts with Wuhan markets or animals, although two had visited a Wuhan hospital. Five family members (aged 36-66 years) presented with fever, upper or lower respiratory tract symptoms, or diarrhoea, or a combination of these 3-6 days after exposure. They presented to our hospital (The University of Hong Kong-Shenzhen Hospital, Shenzhen) 6-10 days after symptom onset. They and one asymptomatic child (aged 10 years) had radiological ground-glass lung opacities. Older patients (aged >60 years) had more systemic symptoms, extensive radiological ground-glass lung changes, lymphopenia, thrombocytopenia, and increased C-reactive protein and lactate dehydrogenase levels. The nasopharyngeal or throat swabs of these six patients were negative for known respiratory microbes by point-of-care multiplex RT-PCR, but five patients (four adults and the child) were RT-PCR positive for genes encoding the internal RNA-dependent RNA polymerase and surface Spike protein of this novel coronavirus, which were confirmed by Sanger sequencing. Phylogenetic analysis of these five patients' RT-PCR amplicons and two full genomes by next-generation sequencing showed that this is a novel coronavirus, which is closest to the bat severe acute respiatory syndrome (SARS)-related coronaviruses found in Chinese horseshoe bats.

INTERPRETATION

Our findings are consistent with person-to-person transmission of this novel coronavirus in hospital and family settings, and the reports of infected travellers in other geographical regions.

FUNDING

The Shaw Foundation Hong Kong, Michael Seak-Kan Tong, Respiratory Viral Research Foundation Limited, Hui Ming, Hui Hoy and Chow Sin Lan Charity Fund Limited, Marina Man-Wai Lee, the Hong Kong Hainan Commercial Association South China Microbiology Research Fund, Sanming Project of Medicine (Shenzhen), and High Level-Hospital Program (Guangdong Health Commission).

A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster

BACKGROUND

An ongoing outbreak of pneumonia associated with a novel coronavirus was reported in Wuhan city, Hubei province, China. Affected patients were geographically linked with a local wet market as a potential source. No data on person-to-person or nosocomial transmission have been published to date.

METHODS

In this study, we report the epidemiological, clinical, laboratory, radiological, and microbiological findings of five patients in a family cluster who presented with unexplained pneumonia after returning to Shenzhen, Guangdong province, China, after a visit to Wuhan, and an additional family member who did not travel to Wuhan. Phylogenetic analysis of genetic sequences from these patients were done.

FINDINGS

From Jan 10, 2020, we enrolled a family of six patients who travelled to Wuhan from Shenzhen between Dec 29, 2019 and Jan 4, 2020. Of six family members who travelled to Wuhan, five were identified as infected with the novel coronavirus. Additionally, one family member, who did not travel to Wuhan, became infected with the virus after several days of contact with four of the family members. None of the family members had contacts with Wuhan markets or animals, although two had visited a Wuhan hospital. Five family members (aged 36-66 years) presented with fever, upper or lower respiratory tract symptoms, or diarrhoea, or a combination of these 3-6 days after exposure. They presented to our hospital (The University of Hong Kong-Shenzhen Hospital, Shenzhen) 6-10 days after symptom onset. They and one asymptomatic child (aged 10 years) had radiological ground-glass lung opacities. Older patients (aged >60 years) had more systemic symptoms, extensive radiological ground-glass lung changes, lymphopenia, thrombocytopenia, and increased C-reactive protein and lactate dehydrogenase levels. The nasopharyngeal or throat swabs of these six patients were negative for known respiratory microbes by point-of-care multiplex RT-PCR, but five patients (four adults and the child) were RT-PCR positive for genes encoding the internal RNA-dependent RNA polymerase and surface Spike protein of this novel coronavirus, which were confirmed by Sanger sequencing. Phylogenetic analysis of these five patients' RT-PCR amplicons and two full genomes by next-generation sequencing showed that this is a novel coronavirus, which is closest to the bat severe acute respiatory syndrome (SARS)-related coronaviruses found in Chinese horseshoe bats.

INTERPRETATION

Our findings are consistent with person-to-person transmission of this novel coronavirus in hospital and family settings, and the reports of infected travellers in other geographical regions.

FUNDING

The Shaw Foundation Hong Kong, Michael Seak-Kan Tong, Respiratory Viral Research Foundation Limited, Hui Ming, Hui Hoy and Chow Sin Lan Charity Fund Limited, Marina Man-Wai Lee, the Hong Kong Hainan Commercial Association South China Microbiology Research Fund, Sanming Project of Medicine (Shenzhen), and High Level-Hospital Program (Guangdong Health Commission).

A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster

BACKGROUND

An ongoing outbreak of pneumonia associated with a novel coronavirus was reported in Wuhan city, Hubei province, China. Affected patients were geographically linked with a local wet market as a potential source. No data on person-to-person or nosocomial transmission have been published to date.

METHODS

In this study, we report the epidemiological, clinical, laboratory, radiological, and microbiological findings of five patients in a family cluster who presented with unexplained pneumonia after returning to Shenzhen, Guangdong province, China, after a visit to Wuhan, and an additional family member who did not travel to Wuhan. Phylogenetic analysis of genetic sequences from these patients were done.

FINDINGS

From Jan 10, 2020, we enrolled a family of six patients who travelled to Wuhan from Shenzhen between Dec 29, 2019 and Jan 4, 2020. Of six family members who travelled to Wuhan, five were identified as infected with the novel coronavirus. Additionally, one family member, who did not travel to Wuhan, became infected with the virus after several days of contact with four of the family members. None of the family members had contacts with Wuhan markets or animals, although two had visited a Wuhan hospital. Five family members (aged 36-66 years) presented with fever, upper or lower respiratory tract symptoms, or diarrhoea, or a combination of these 3-6 days after exposure. They presented to our hospital (The University of Hong Kong-Shenzhen Hospital, Shenzhen) 6-10 days after symptom onset. They and one asymptomatic child (aged 10 years) had radiological ground-glass lung opacities. Older patients (aged >60 years) had more systemic symptoms, extensive radiological ground-glass lung changes, lymphopenia, thrombocytopenia, and increased C-reactive protein and lactate dehydrogenase levels. The nasopharyngeal or throat swabs of these six patients were negative for known respiratory microbes by point-of-care multiplex RT-PCR, but five patients (four adults and the child) were RT-PCR positive for genes encoding the internal RNA-dependent RNA polymerase and surface Spike protein of this novel coronavirus, which were confirmed by Sanger sequencing. Phylogenetic analysis of these five patients' RT-PCR amplicons and two full genomes by next-generation sequencing showed that this is a novel coronavirus, which is closest to the bat severe acute respiatory syndrome (SARS)-related coronaviruses found in Chinese horseshoe bats.

INTERPRETATION

Our findings are consistent with person-to-person transmission of this novel coronavirus in hospital and family settings, and the reports of infected travellers in other geographical regions.

FUNDING

The Shaw Foundation Hong Kong, Michael Seak-Kan Tong, Respiratory Viral Research Foundation Limited, Hui Ming, Hui Hoy and Chow Sin Lan Charity Fund Limited, Marina Man-Wai Lee, the Hong Kong Hainan Commercial Association South China Microbiology Research Fund, Sanming Project of Medicine (Shenzhen), and High Level-Hospital Program (Guangdong Health Commission).

Emerging horizon for bat borne viral zoonoses

Bats are the only flying placental mammals that constitute the second largest order of mammals and present all around the world except in Arctic, Antarctica and a few oceanic islands. Sixty percent of emerging infectious diseases originating from animals are zoonotic and more than two-thirds of them originate in wildlife. Bats were evolved as a super-mammal for harboring many of the newly identified deadly diseases without any signs and lesions. Their unique ability to fly, particular diet, roosting behavior, long life span, ability to echolocate and critical susceptibility to pathogens make them suitable host to harbor numerous zoonotic pathogens like virus, bacteria and parasite. Many factors are responsible for the emergence of bat borne zoonoses but the most precipitating factor is human intrusions. Deforestation declined the natural habitat and forced the bats and other wild life to move out of their niche. These stressed bats, having lost foraging and behavioral pattern invade in proximity of human habitation. Either directly or indirectly they transmit the viruses to humans and animals. Development of fast detection modern techniques for viruses from the diseased and environmental samples and the lessons learned in the past helped in preventing the severity during the latest outbreaks.

Diversity and Evolution of Viral Pathogen Community in Cave Nectar Bats (Eonycteris spelaea)

Bats are unique mammals, exhibit distinctive life history traits and have unique immunological approaches to suppression of viral diseases upon infection. High-throughput next-generation sequencing has been used in characterizing the virome of different bat species. The cave nectar bat, Eonycteris spelaea, has a broad geographical range across Southeast Asia, India and southern China, however, little is known about their involvement in virus transmission. Here we investigate the diversity and abundance of viral communities from a colony of Eonycteris spelaea residing in Singapore. Our results detected 47 and 22 different virus families from bat fecal and urine samples, respectively. Among these, we identify a large number of virus families including Adenoviridae, Flaviviridae, Reoviridae, Papillomaviridae, Paramyxoviridae, Parvoviridae, Picornaviridae, and Polyomaviridae. In most cases, viral sequences from Eonycteris spelaea are genetically related to a group of bat viruses from other bat genera (e.g., Eidolon, Miniopterus, Rhinolophus and Rousettus). The results of this study improve our knowledge of the host range, spread and evolution of several important viral pathogens. More significantly, our findings provide a baseline to study the temporal patterns of virus shedding and how they correlate with bat phenological trends.

Detection of adenovirus, papillomavirus and parvovirus in Brazilian bats of the species Artibeus lituratus and Sturnira lilium

Bats play a significant role in maintaining their ecosystems through pollination, dispersal of seeds, and control of insect populations, but they are also known to host many microorganisms and have been described as natural reservoirs for viruses with zoonotic potential. The diversity of viruses in these animals remains largely unknown, however, because studies are limited by species, location, virus target, or sample type. Therefore, the aim of this study was to detect fragments of viral genomes in bat samples. We performed high-throughput sequencing analysis and specific PCR and RT-PCR on pools of anal and oropharyngeal swabs from Artibeus lituratus and Sturnira lilium collected in southern Brazil. As a result, a member of the family Adenoviridae related to human adenovirus C was detected in anal swabs from S. lilium. In addition, we detected a papillomavirus in an anal swab from A. lituratus. Our analyses also allowed the detection of adenoviruses and parvoviruses in oropharyngeal swabs collected from A. lituratus. These results increase our knowledge about viral diversity and illustrate the importance of conducting virus surveillance in bats.

Molecular identification of Betacoronavirus in bats from Sardinia (Italy): first detection and phylogeny

Bats may be natural reservoirs for a large variety of emerging viruses, including mammalian coronaviruses (CoV). The recent emergence of severe acute respiratory syndrome-associated coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) in humans, with evidence that these viruses may have their ancestry in bats, highlights the importance of virus surveillance in bat populations. Here, we report the identification and molecular characterization of a bat β-Coronavirus, detected during a viral survey carried out on different bat species in the island of Sardinia (Italy). Cutaneous, oral swabs, and faecal samples were collected from 46 bats, belonging to 15 different species, and tested for viral presence. Coronavirus RNA was detected in faecal samples from three different species: the greater horseshoe bat (Rhinolophus ferrumequinum), the brown long-eared bat (Plecotus auritus), and the European free-tailed bat (Tadarida teniotis). Phylogenetic analyses based on RNA-dependent RNA polymerase (RdRp) sequences assigned the detected CoV to clade 2b within betacoronaviruses, clustering with SARS-like bat CoVs previously reported. These findings point to the need for continued surveillance of bat CoV circulating in Sardinian bats, and extend the current knowledge on CoV ecology with novel sequences detected in bat species not previously described as β-Coronavirus hosts.

Using noninvasive metagenomics to characterize viral communities from wildlife

Microbial communities play an important role in organismal and ecosystem health. While high-throughput metabarcoding has revolutionized the study of bacterial communities, generating comparable viral communities has proven elusive, particularly in wildlife samples where the diversity of viruses and limited quantities of viral nucleic acid present distinctive challenges. Metagenomic sequencing is a promising solution for studying viral communities, but the lack of standardized methods currently precludes comparisons across host taxa or localities. Here, we developed an untargeted shotgun metagenomic sequencing protocol to generate comparable viral communities from noninvasively collected faecal and oropharyngeal swabs. Using samples from common vampire bats (Desmodus rotundus), a key species for virus transmission to humans and domestic animals, we tested how different storage media, nucleic acid extraction procedures and enrichment steps affect viral community detection. Based on finding viral contamination in foetal bovine serum, we recommend storing swabs in RNAlater or another nonbiological medium. We recommend extracting nucleic acid directly from swabs rather than from supernatant or pelleted material, which had undetectable levels of viral RNA. Results from a low-input RNA library preparation protocol suggest that ribosomal RNA depletion and light DNase treatment reduce host and bacterial nucleic acid, and improve virus detection. Finally, applying our approach to twelve pooled samples from seven localities in Peru, we showed that detected viral communities saturated at the attained sequencing depth, allowing unbiased comparisons of viral community composition. Future studies using the methods outlined here will elucidate the determinants of viral communities across host species, environments and time.

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

Environmental disturbances in the Neotropics (e.g., deforestation, agriculture intensification, urbanization) contribute to an increasing risk of cross-species transmission of microorganisms and to disease outbreaks due to changing ecosystems of reservoir hosts. Although Amazonia encompasses the greatest diversity of reservoir species, the outsized viral population diversity (virome) has yet to be investigated. Here, through a metagenomic approach, we identified 10,991 viral sequences in the saliva and feces of two bat species, Desmodus rotundus (hematophagous), trapped in two different caves surrounded by primary lowland forest, and Molossus molossus (insectivorous), trapped in forest and urban habitats. These sequences are related to 51 viral families known to infect a wide range of hosts (i.e., bacteria, plants, insects and vertebrates). Most viruses detected reflected the diet of bat species, with a high proportion of plant and insect-related viral families for M. molossus and a high proportion of vertebrate-related viral families for D. rotundus, highlighting its influence in shaping the viral diversity of bats. Lastly, we reconstructed the phylogenetic relationships for five vertebrate-related viral families (Nairoviridae, Circoviridae, Retroviridae, Herpesviridae, Papillomaviridae). The results showed highly supported clustering with other viral sequences of the same viral family hosted by other bat species, highlighting the potential association of viral diversity with the host’s diet. These findings provide significant insight into viral bat diversity in French Guiana belonging to the Amazonian biome and emphasize that habitats and the host’s dietary ecology may drive the viral diversity in the bat communities investigated.

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

Accumulating studies have shown that bats could harbor various important pathogenic viruses that could be transmitted to humans and other animals. Extensive metagenomic studies of different organs/tissues from bats have revealed a large number of novel or divergent viruses. To elucidate viral diversity and epidemiological and phylogenetic characteristics, six pooled fecal samples from bats were generated (based on bat species and geographic regions characteristic for virome analysis). These contained 500 fecal samples from six bat species, collected in four geographic regions. Metagenomic analysis revealed a plethora of divergent viruses originally found in bats. Multiple contigs from influenza A virus and coronaviruses in bats shared high identity with those from humans, suggesting possible cross-species transmission, whereas a number of contigs, whose sequences were taxonomically classifiable within Alphapapillomavirus, Betaretrovirus, Alpharetrovirus, Varicellovirus, Cyprinivirus, Chlorovirus and Cucumovirus had low identity to viruses in existing databases, which indicated possible evolution of novel viral species. None of the established caliciviruses and picornaviruses were found in the 500 fecal specimens. Papillomaviruses with high amino acid identity were found in Scotophilus kuhlii and Rhinolophus blythi, challenging the hypotheses regarding the strict host specificity and co-evolution of papillomaviruses. Phylogenetic analysis showed that four bat rotavirus A strains might be tentative G3 strains, according to the Rotavirus Classification Working Group classification.

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

A number of PVs have been described in bats but to the best of our knowledge not from feces. Using a previously described NetoVIR protocol, Eidolon helvum pooled fecal samples (Eh) were treated and sequenced by Illumina next generation sequencing technology. Two complete genomes of novel PVs (EhPV2 and EhPV3) and 3 partial sequences (BATPV61, BATPV890a and BATPV890b) were obtained and analysis showed that the EhPV2 and EhPV3 major capsid proteins cluster with and share 60-64% nucleotide identity with that of Rousettus aegyptiacus PV1, thus representing new species of PVs within the genus Psipapillomavirus. The other PVs clustered in different branches of our phylogenetic tree and may potentially represent novel species and/or genera. This points to the vast diversity of PVs in bats and in Eidolon helvum bats in particular, therefore adding support to the current concept that PV evolution is more complex than merely strict PV-host co-evolution.

Molecular methods for identification and characterization of novel papillomaviruses

Papillomaviruses (PV) are a remarkably heterogeneous family of small DNA viruses that infect a wide variety of vertebrate species and are aetiologically linked with the development of various neoplastic changes of the skin and mucosal epithelia. Based on nucleotide similarity, PVs are hierarchically classified into genera, species and types. Novel human PV (HPV) types are given a unique number only after the whole genome has been cloned and deposited with the International HPV Reference Center. As of 9 March 2015, 200 different HPV types, belonging to 49 species, had been recognized by the International HPV Reference Center. In addition, 131 animal PV types identified from 66 different animal species exist. Recent advances in molecular techniques have resulted in an explosive increase in the identification of novel HPV types and novel subgenomic HPV sequences in the last few years. Among PV genera, the γ-PV genus has been growing most rapidly in recent years with 80 completely sequenced HPV types, followed by α-PV and β-PV genera that have 65 and 51 recognized HPV types, respectively. We reviewed in detail the contemporary molecular methods most often used for identification and characterization of novel PV types, including PCR, rolling circle amplification and next-generation sequencing. Furthermore, we present a short overview of 12 and 10 novel HPV types recently identified in Sweden and Slovenia, respectively. Finally, an update on the International Human Papillomavirus Reference Center is provided.

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

Bats harbour a diverse array of viruses, including significant human pathogens. Extensive metagenomic studies of material from bats, in particular guano, have revealed a large number of novel or divergent viral taxa that were previously unknown. New Zealand has only two extant indigenous terrestrial mammals, which are both bats, Mystacina tuberculata (the lesser short-tailed bat) and Chalinolobus tuberculatus (the long-tailed bat). Until the human introduction of exotic mammals, these species had been isolated from all other terrestrial mammals for over 1 million years (potentially over 16 million years for M. tuberculata). Four bat guano samples were collected from M. tuberculata roosts on the isolated offshore island of Whenua hou (Codfish Island) in New Zealand. Metagenomic analysis revealed that this species still hosts a plethora of divergent viruses. Whilst the majority of viruses detected were likely to be of dietary origin, some putative vertebrate virus sequences were identified. Papillomavirus, polyomavirus, calicivirus and hepevirus were found in the metagenomic data and subsequently confirmed using independent PCR assays and sequencing. The new hepevirus and calicivirus sequences may represent new genera within these viral families. Our findings may provide an insight into the origins of viral families, given their detection in an isolated host species.

A systematic approach to novel virus discovery in emerging infectious disease outbreaks

The discovery of novel viruses is of great importance to human health-both in the setting of emerging infectious disease outbreaks and in disease syndromes of unknown etiology. Despite the recent proliferation of many efficient virus discovery methods, careful selection of a combination of methods is important to demonstrate a novel virus, its clinical associations, and its relevance in a timely manner. The identification of a patient or an outbreak with distinctive clinical features and negative routine microbiological workup is often the starting point for virus hunting. This review appraises the roles of culture, electron microscopy, and nucleic acid detection-based methods in optimizing virus discovery. Cell culture is generally slow but may yield viable virus. Although the choice of cell line often involves trial and error, it may be guided by the clinical syndrome. Electron microscopy is insensitive but fast, and may provide morphological clues to choice of cell line or consensus primers for nucleic acid detection. Consensus primer PCR can be used to detect viruses that are closely related to known virus families. Random primer amplification and high-throughput sequencing can catch any virus genome but cannot yield an infectious virion for testing Koch postulates. A systematic approach that incorporates carefully chosen combinations of virus detection techniques is required for successful virus discovery.

Immune status, strain background, and anatomic site of inoculation affect mouse papillomavirus (MmuPV1) induction of exophytic papillomas or endophytic trichoblastomas

Papillomaviruses (PVs) induce papillomas, premalignant lesions, and carcinomas in a wide variety of species. PVs are classified first based on their host and tissue tropism and then their genomic diversities. A laboratory mouse papillomavirus, MmuPV1 (formerly MusPV), was horizontally transmitted within an inbred colony of NMRI-Foxn1(nu)/Foxn1nu (nude; T cell deficient) mice of an unknown period of time. A ground-up, filtered papilloma inoculum was not capable of infecting C57BL/6J wild-type mice; however, immunocompetent, alopecic, S/RV/Cri-ba/ba (bare) mice developed small papillomas at injection sites that regressed. NMRI-Foxn1(nu) and B6.Cg-Foxn1(nu), but not NU/J-Foxn1(nu), mice were susceptible to MmuPV1 infection. B6 congenic strains, but not other congenic strains carrying the same allelic mutations, lacking B- and T-cells, but not B-cells alone, were susceptible to infection, indicating that mouse strain and T-cell deficiency are critical to tumor formation. Lesions initially observed were exophytic papillomas around the muzzle, exophytic papillomas on the tail, and condylomas of the vaginal lining which could be induced by separate scarification or simultaneous scarification of MmuPV1 at all four sites. On the dorsal skin, locally invasive, poorly differentiated tumors developed with features similar to human trichoblastomas. Transcriptome analysis revealed significant differences between the normal skin in these anatomic sites and in papillomas versus trichoblastomas. The primarily dysregulated genes involved molecular pathways associated with cancer, cellular development, cellular growth and proliferation, cell morphology, and connective tissue development and function. Although trichoepitheliomas are benign, aggressive tumors, few of the genes commonly associated with basal cell carcinoma or squamous cells carcinoma were highly dysregulated.

A nairovirus isolated from African bats causes haemorrhagic gastroenteritis and severe hepatic disease in mice

Bats can carry important zoonotic pathogens. Here we use a combination of next-generation sequencing and classical virus isolation methods to identify novel nairoviruses from bats captured from a cave in Zambia. This nairovirus infection is highly prevalent among giant leaf-nosed bats, Hipposideros gigas (detected in samples from 16 individuals out of 38). Whole-genome analysis of three viral isolates (11SB17, 11SB19 and 11SB23) reveals a typical bunyavirus tri-segmented genome. The strains form a single phylogenetic clade that is divergent from other known nairoviruses, and are hereafter designated as Leopards Hill virus (LPHV). When i.p. injected into mice, the 11SB17 strain causes only slight body weight loss, whereas 11SB23 produces acute and lethal disease closely resembling that observed with Crimean–Congo Haemorrhagic Fever virus in humans. We believe that our LPHV mouse model will be useful for research on the pathogenesis of nairoviral haemorrhagic disease.

Bats carry viruses that can cause disease in other animals and in humans. Here, Ishii et al. identify new nairoviruses from African bats and show that some of them can produce a severe haemorrhagic disease in laboratory mice that is similar to Crimean–Congo haemorrhagic fever in humans.

Novel papillomaviruses in free-ranging Iberian bats: no virus-host co-evolution, no strict host specificity, and hints for recombination

Papillomaviruses (PVs) are widespread pathogens. However, the extent of PV infections in bats remains largely unknown. This work represents the first comprehensive study of PVs in Iberian bats. We identified four novel PVs in the mucosa of free-ranging Eptesicus serotinus (EserPV1, EserPV2, and EserPV3) and Rhinolophus ferrumequinum (RferPV1) individuals and analyzed their phylogenetic relationships within the viral family. We further assessed their prevalence in different populations of E. serotinus and its close relative E. isabellinus. Although it is frequent to read that PVs co-evolve with their host, that PVs are highly species-specific, and that PVs do not usually recombine, our results suggest otherwise. First, strict virus–host co-evolution is rejected by the existence of five, distantly related bat PV lineages and by the lack of congruence between bats and bat PVs phylogenies. Second, the ability of EserPV2 and EserPV3 to infect two different bat species (E. serotinus and E. isabellinus) argues against strict host specificity. Finally, the description of a second noncoding region in the RferPV1 genome reinforces the view of an increased susceptibility to recombination in the E2-L2 genomic region. These findings prompt the question of whether the prevailing paradigms regarding PVs evolution should be reconsidered.

Viral metagenomics analysis demonstrates the diversity of viral flora in piglet diarrhoeic faeces in China

To investigate the diversity of viral flora, we used metagenomics to study the viral communities in a pooled faecal sample of 27 diarrhoeic piglets from intensive commercial farms in China. The 15 distinct mammalian viruses identified in the pooled diarrhoeic sample were, in order of abundance of nucleic acid sequence, Porcine epidemic diarrhea virus (PEDV), sapovirus, porcine bocavirus-4 (PBoV-4), sapelovirus, torovirus, coronavirus, PBoV-2, stool-associated single-stranded DNA virus (poSCV), astrovirus (AstV), kobuvirus, posavirus-1, porcine enterovirus-9 (PEV-9), porcine circovirus-like (po-circo-like) virus, picobirnavirus (PBV) and Torque teno sus virus 2 (TTSuV-2). The prevalence rate of each virus was verified from diarrhoeic and healthy piglets by PCR assay. A mean of 5.5 different viruses were shed in diarrhoeic piglets, and one piglet was in fact co-infected with 11 different viruses. By contrast, healthy piglets shed a mean of 3.2 different viruses. Compared with samples from healthy piglets, the co-infection of PEDV and PBoV had a high prevalence rate in diarrhoea samples, suggesting a correlation with the appearance of diarrhoea in piglets. Furthermore, we report here for the first time the presence of several recently described viruses in China, and the identification of novel genotypes. Therefore, our investigation results provide an unbiased survey of viral communities and prevalence in faecal samples of piglets.

Viral metagenomics on animals as a tool for the detection of zoonoses prior to human infection?

Many human viral infections have a zoonotic, i.e., wild or domestic animal, origin. Several zoonotic viruses are transmitted to humans directly via contact with an animal or indirectly via exposure to the urine or feces of infected animals or the bite of a bloodsucking arthropod. If a virus is able to adapt and replicate in its new human host, human-to-human transmissions may occur, possibly resulting in an epidemic, such as the A/H1N1 flu pandemic in 2009. Thus, predicting emerging zoonotic infections is an important challenge for public health officials in the coming decades. The recent development of viral metagenomics, i.e., the characterization of the complete viral diversity isolated from an organism or an environment using high-throughput sequencing technologies, is promising for the surveillance of such diseases and can be accomplished by analyzing the viromes of selected animals and arthropods that are closely in contact with humans. In this review, we summarize our current knowledge of viral diversity within such animals (in particular blood-feeding arthropods, wildlife and domestic animals) using metagenomics and present its possible future application for the surveillance of zoonotic and arboviral diseases.

DBatVir: the database of bat-associated viruses

Emerging infectious diseases remain a significant threat to public health. Most emerging infectious disease agents in humans are of zoonotic origin. Bats are important reservoir hosts of many highly lethal zoonotic viruses and have been implicated in numerous emerging infectious disease events in recent years. It is essential to enhance our knowledge and understanding of the genetic diversity of the bat-associated viruses to prevent future outbreaks. To facilitate further research, we constructed the database of bat-associated viruses (DBatVir). Known viral sequences detected in bat samples were manually collected and curated, along with the related metadata, such as the sampling time, location, bat species and specimen type. Additional information concerning the bats, including common names, diet type, geographic distribution and phylogeny were integrated into the database to bridge the gap between virologists and zoologists. The database currently covers >4100 bat-associated animal viruses of 23 viral families detected from 196 bat species in 69 countries worldwide. It provides an overview and snapshot of the current research regarding bat-associated viruses, which is essential now that the field is rapidly expanding. With a user-friendly interface and integrated online bioinformatics tools, DBatVir provides a convenient and powerful platform for virologists and zoologists to analyze the virome diversity of bats, as well as for epidemiologists and public health researchers to monitor and track current and future bat-related infectious diseases.

Database URL:http://www.mgc.ac.cn/DBatVir/

The role of metagenomics in understanding the human microbiome in health and disease

The term microbiome refers to the genetic material of the catalog of microbial taxa associated with humans. As in all ecosystems, the microbiota reaches a dynamic equilibrium in the human body, which can be altered by environmental factors and external stimuli. Metagenomics is a relatively new field of study of microbial genomes within diverse environmental samples, which is of increasing importance in microbiology. The introduction of this ecological perception of microbiology is the key to achieving real knowledge about the influence of the microbiota in human health and disease. The aim of this review is to summarize the link between the human microbiota (focusing on the intestinal, vaginal, skin, and airway body sites) and health from this ecological point of view, highlighting the contribution of metagenomics in the advance of this field.

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

The prediction of viral zoonosis epidemics has become a major public health issue. A profound understanding of the viral population in key animal species acting as reservoirs represents an important step towards this goal. Bats harbor diverse viruses, some of which are of particular interest because they cause severe human diseases. However, little is known about the diversity of the global population of viruses found in bats (virome). We determined the viral diversity of five different French insectivorous bat species (nine specimens in total) in close contact with humans. Sequence-independent amplification, high-throughput sequencing with Illumina technology and a dedicated bioinformatics analysis pipeline were used on pooled tissues (brain, liver and lungs). Comparisons of the sequences of contigs and unassembled reads provided a global taxonomic distribution of virus-related sequences for each sample, highlighting differences both within and between bat species. Many viral families were present in these viromes, including viruses known to infect bacteria, plants/fungi, insects or vertebrates, the most relevant being those infecting mammals (Retroviridae, Herpesviridae, Bunyaviridae, Poxviridae, Flaviviridae, Reoviridae, Bornaviridae, Picobirnaviridae). In particular, we detected several new mammalian viruses, including rotaviruses, gammaretroviruses, bornaviruses and bunyaviruses with the identification of the first bat nairovirus. These observations demonstrate that bats naturally harbor viruses from many different families, most of which infect mammals. They may therefore constitute a major reservoir of viral diversity that should be analyzed carefully, to determine the role played by bats in the spread of zoonotic viral infections.

Characterization of a novel G3P[3] rotavirus isolated from a lesser horseshoe bat: a distant relative of feline/canine rotaviruses

Bats are considered important animal reservoirs for many viruses pathogenic to humans. An approach based on viral metagenomics was used to study gut specimens from 78 insectivorous bats in Yunnan Province, China. Seventy-four reads were found to be related to group A rotavirus (RVA). Further reverse transcription-PCR screening and viral isolation on cell cultures confirmed the presence of a novel RVA strain, named RVA/Bat-tc/MSLH14/2012/G3P[3], in 1 (6%) of 16 lesser horseshoe bats. Full genomic sequencing analyses showed that MSLH14 possessed the genotype constellation G3-P[3]-I8-R3-C3-M3-A9-N3-T3-E3-H6, which is akin to human and animal rotaviruses believed to be of feline/canine origin. Phylogenetic analysis indicated that VP7 was most closely related to bovine RVA strains from India, whereas VP4 was most closely related to an unusual human RVA strain, CMH222, with animal characteristics isolated in Thailand. The remaining gene segments were only distantly related to a range of animal RVA strains, most of which are believed to be related to feline/canine RVAs. Experimental infection showed that bat RVA strain MSLH14 was highly pathogenic to suckling mice, causing 100% mortality when they were inoculated orally with a titer as low as 5 × 10² 50% tissue culture infective doses. As this virus is not closely related to any known RVA strain, it is tempting to speculate that it is a true bat RVA strain rather than a virus transmitted between species. However, further screening of bat populations, preferably juvenile animals, will be crucial in determining whether or not this virus is widely distributed in the bat population.

Romance of the three domains: how cladistics transformed the classification of cellular organisms

Cladistics is a biological philosophy that uses genealogical relationship among species and an inferred sequence of divergence as the basis of classification. This review critically surveys the chronological development of biological classification from Aristotle through our postgenomic era with a central focus on cladistics. In 1957, Julian Huxley coined cladogenesis to denote splitting from subspeciation. In 1960, the English translation of Willi Hennig’s 1950 work, Systematic Phylogenetics, was published, which received strong opposition from pheneticists, such as numerical taxonomists Peter Sneath and Robert Sokal, and evolutionary taxonomist, Ernst Mayr, and sparked acrimonious debates in 1960–1980. In 1977–1990, Carl Woese pioneered in using small subunit rRNA gene sequences to delimitate the three domains of cellular life and established major prokaryotic phyla. Cladistics has since dominated taxonomy. Despite being compatible with modern microbiological observations, i.e. organisms with unusual phenotypes, restricted expression of characteristics and occasionally being uncultivable, increasing recognition of pervasiveness and abundance of horizontal gene transfer has challenged relevance and validity of cladistics. The mosaic nature of eukaryotic and prokaryotic genomes was also gradually discovered. In the mid-2000s, high-throughput and whole-genome sequencing became routine and complex geneologies of organisms have led to the proposal of a reticulated web of life. While genomics only indirectly leads to understanding of functional adaptations to ecological niches, computational modeling of entire organisms is underway and the gap between genomics and phenetics may soon be bridged. Controversies are not expected to settle as taxonomic classifications shall remain subjective to serve the human scientist, not the classified.

Next generation sequencing of viral RNA genomes

Background

With the advent of Next Generation Sequencing (NGS) technologies, the ability to generate large amounts of sequence data has revolutionized the genomics field. Most RNA viruses have relatively small genomes in comparison to other organisms and as such, would appear to be an obvious success story for the use of NGS technologies. However, due to the relatively low abundance of viral RNA in relation to host RNA, RNA viruses have proved relatively difficult to sequence using NGS technologies. Here we detail a simple, robust methodology, without the use of ultra-centrifugation, filtration or viral enrichment protocols, to prepare RNA from diagnostic clinical tissue samples, cell monolayers and tissue culture supernatant, for subsequent sequencing on the Roche 454 platform.

Results

As representative RNA viruses, full genome sequence was successfully obtained from known lyssaviruses belonging to recognized species and a novel lyssavirus species using these protocols and assembling the reads using de novo algorithms. Furthermore, genome sequences were generated from considerably less than 200 ng RNA, indicating that manufacturers’ minimum template guidance is conservative. In addition to obtaining genome consensus sequence, a high proportion of SNPs (Single Nucleotide Polymorphisms) were identified in the majority of samples analyzed.

Conclusions

The approaches reported clearly facilitate successful full genome lyssavirus sequencing and can be universally applied to discovering and obtaining consensus genome sequences of RNA viruses from a variety of sources.

Safety analysis of a Russian phage cocktail: from metagenomic analysis to oral application in healthy human subjects

Phage therapy has a long tradition in Eastern Europe, where preparations are comprised of complex phage cocktails whose compositions have not been described. We investigated the composition of a phage cocktail from the Russian pharmaceutical company Microgen targeting Escherichia coli/Proteus infections. Electron microscopy identified six phage types, with numerically T7-like phages dominating over T4-like phages. A metagenomic approach using taxonomical classification, reference mapping and de novo assembly identified 18 distinct phage types, including 7 genera of Podoviridae, 2 established and 2 proposed genera of Myoviridae, and 2 genera of Siphoviridae. De novo assembly yielded 7 contigs greater than 30 kb, including a 147-kb Myovirus genome and a 42-kb genome of a potentially new phage. Bioinformatic analysis did not reveal undesired genes and a small human volunteer trial did not associate adverse effects with oral phage exposure.

Animal papillomaviruses

We provide an overview of the host range, taxonomic classification and genomic diversity of animal papillomaviruses. The complete genomes of 112 non-human papillomavirus types, recovered from 54 different host species, are currently available in GenBank. The recent characterizations of reptilian papillomaviruses extend the host range of the Papillomaviridae to include all amniotes. Although the genetically diverse papillomaviruses have a highly conserved genomic lay-out, deviations from this prototypic genome organization are observed in several animal papillomaviruses, and only the core ORFs E1, E2, L2 and L1 are present in all characterized papillomavirus genomes. The discovery of papilloma-polyoma hybrids BPCV1 and BPCV2, containing a papillomaviral late region but an early region encoding typical polyomaviral nonstructural proteins, and the detection of recombination breakpoints between the early and late coding regions of cetacean papillomaviruses, could indicate that early and late gene cassettes of papillomaviruses are relatively independent entities that can be interchanged by recombination.

Identification of a novel human papillomavirus by metagenomic analysis of samples from patients with febrile respiratory illness

As part of a virus discovery investigation using a metagenomic approach, a highly divergent novel Human papillomavirus type was identified in pooled convenience nasal/oropharyngeal swab samples collected from patients with febrile respiratory illness. Phylogenetic analysis of the whole genome and the L1 gene reveals that the new HPV identified in this study clusters with previously described gamma papillomaviruses, sharing only 61.1% (whole genome) and 63.1% (L1) sequence identity with its closest relative in the Papillomavirus episteme (PAVE) database. This new virus was named HPV_SD2 pending official classification. The complete genome of HPV-SD2 is 7,299 bp long (36.3% G/C) and contains 7 open reading frames (L2, L1, E6, E7, E1, E2 and E4) and a non-coding long control region (LCR) between L1 and E6. The metagenomic procedures, coupled with the bioinformatic methods described herein are well suited to detect small circular genomes such as those of human papillomaviruses.

Multiple evolutionary origins of bat papillomaviruses

Infection by papillomaviruses (PVs) has been linked to different types of neoplasias, in both human and non-human hosts. Knowledge about PV diversity is essential to reliably infer the evolutionary history of these pathogens and to elucidate the link between infection and disease. We cloned and sequenced the complete genome of a novel PV, EhelPV1, isolated from hair bulbs from a captive straw-colored fruit bat Eidolon helvum (Pteropodidae, Chiroptera). We also retrieved partial sequences of the E1 and L1 genes from hair bulbs from a captive Indian flying fox Pteropus giganteus (Pteropodidae, Chiroptera). The detected virus (PgigPV1) presumably corresponded to a novel type as well. Maximum likelihood phylogenetic analyses were conducted using a representative collection of 132 PVs. EhelPV1 belonged to the Lambda+Mu-PV crown group and was most closely related to another bat PV, MschPV2. Both fragments of PgigPV1 were placed alongside with EhelPV1. The novel PVs were phylogenetically distant from other previously described bat PVs, namely MrPV1, MschPV1 and RaPV1. We have further characterized the sequence patterns of the E2-binding sites occurring in the upstream regulatory region of Lambda+Mu-PVs. Common fingerprints within this region are shared by certain PVs. However, there is not a sharp correspondence between the repertoire of transcription factor binding sites in the viral regulatory region and host range, tissue tropism or viral life style. Our results reinforce the hypothesis that PVs have undergone an initial radiation prior to the divergence of the mammalian hosts, giving rise to the present-day PV crown groups.