This is the first country-wide surveillance of bat-borne viruses in Kenya spanning from 2012–2015 covering sites perceived to have medium to high level bat-human interaction. The objective of this surveillance study was to apply a non-invasive approach using fresh feces to detect viruses circulating within the diverse species of Kenyan bats. We screened for both DNA and RNA viruses; specifically, astroviruses (AstVs), adenoviruses (ADVs), caliciviruses (CalVs), coronaviruses (CoVs), flaviviruses, filoviruses, paramyxoviruses (PMVs), polyomaviruses (PYVs) and rotaviruses. We used family-specific primers, amplicon sequencing and further characterization by phylogenetic analysis. Except for filoviruses, eight virus families were detected with varying distributions and positive rates across the five regions (former provinces) studied. AstVs (12.83%), CoVs (3.97%), PMV (2.4%), ADV (2.26%), PYV (1.65%), CalVs (0.29%), rotavirus (0.19%) and flavivirus (0.19%). Novel CalVs were detected in Rousettus aegyptiacus and Mops condylurus while novel Rotavirus-A-related viruses were detected in Taphozous bats and R. aegyptiacus. The two Rotavirus A (RVA) strains detected were highly related to human strains with VP6 genotypes I2 and I16. Genotype I16 has previously been assigned to human RVA-strain B10 from Kenya only, which raises public health concern, particularly considering increased human-bat interaction. Additionally, 229E-like bat CoVs were detected in samples originating from Hipposideros bats roosting in sites with high human activity. Our findings confirm the presence of diverse viruses in Kenyan bats while providing extended knowledge on bat virus distribution. The detection of viruses highly related to human strains and hence of public health concern, underscores the importance of continuous surveillance.
Citation: Cecilia Waruhiu, Sheila Ommeh, Vincent Obanda, Bernard Agwanda, Francis Gakuya, Xing-Yi Ge, Xing-Lou Yang, Li-Jun Wu, Ali Zohaib, Ben Hu, Zheng-Li Shi. Molecular detection of viruses in Kenyan bats and discovery of novel astroviruses, caliciviruses and rotaviruses[J]. VIROLOGICA SINICA, 2017, 32 (2): 101-114 https://doi.org/10.1007/s12250-016-3930-2
Received: 9 December, 2016; Accepted: 15 February 2017; Published: 6 April 2017
Copyright: © The Author(s) 2017. This article is published with open access at Springerlink.com
Data Availability: All relevant data are within the paper and its Supporting Information files.
Corresponding author: Sheila Ommeh, Phone: +254-7-07007028, Email: email@example.com ORCID: 0000-0002-0780-0399; Zheng-Li Shi, Phone: +86-27-87197240, Email: firstname.lastname@example.org ORCID: 0000-0001-8089-163X.
Bats (Order Chiroptera) are a group of highly diverse mammalian species that occurs throughout the world and accounts for more than 25% of mammals (Olival, 2016). Bats occur in large colonies and have evolved unique adaptations such as flight (O’Shea et al., 2014) which allow them to travel over long distances. Over time, bats have extended their roost sites from forested areas, caves and other isolated sites to habitats with substantial human activity. This can be attributed to anthropogenic activities such as urbanization and deforestation that have resulted in destruction of natural bat habitats (Plowright et al., 2015).
In the last decade, bats have been recognized as an important natural reservoir for zoonotic viruses affecting both humans and livestock (Plowright et al., 2015). Both RNA and DNA virus families have been detected in both frugivorous and insectivorous bats including RNA viruses: coronaviruses, filoviruses, lyssaviruses, paramyxoviruses, flaviviruses, astroviruses, caliciviruses; DNA viruses: adenoviruses, herpesviruses, poxviruses, polyomaviruses and papillomavirus among others (Shi, 2013). Over 200 virus species have been reported in bats (Moratelli et al., 2015) and additional novel viruses continue to be detected, yet the global diversity of viruses in bats is still thought to be largely unknown (Anthony et al., 2013). In many instances, some of the viruses harbored in bats such as; Nipah virus, severe acute respiratory syndrome coronavirus (SARS-CoV), multiple strains of lyssavirus, Marburg virus and Ebola virus have spilled over to humans and livestock causing severe disease in different parts of the world (Chua et al., 2000; Leroy et al., 2005; Wang et al., 2005).
Africa is one of the major hotspots of emerging infectious diseases (Jones et al., 2008) hence signifying need for enhance surveillance efforts as well as conservation efforts aimed at mitigating the risks of zoonotic disease emergence. Studies conducted in African countries such as Ghana, Gabon, Cameroon among others have reported circulation of several virus families in bats, such as;Filoviridae, Coronaviridae, Reoviridae, Paramyxoviridae, Adenoviridae and Polyomaviridae (Baker et al., 2013; Corman et al., 2015; Drexler et al., 2009; Yinda et al., 2016). These findings, as well as occurrence of outbreaks linked to bats, such as Ebola virus disease, affirm the recommendation of Jones et al. (2008) , on the necessity for targeted research in Africa as one of the high-risk zones for emergence of zoonotic pathogens (Gatherer, 2014; Marí Saéz et al., 2014). Information of the geographical distribution and genetic diversity of viruses harbored by bats in Kenya is quite limited, yet bat is one of the mammalian orders associated with a growing list of deadly emerging pathogens. Although there have been some studies on bat viruses carried out between 2007 to 2011, they are insufficient to display the epidemiological status of bat viruses in Kenya (Kading et al., 2013; Kuzmin et al., 2010a; Kuzmin et al., 2010b; Tao et al., 2013; Tong et al., 2009). As early as 1982, Marburg virus (MARV) had been linked to the death of a febrile ill patient in Kenya (Smith et al., 1982) and a similar fatal case of MARV infection occurred a few years later in 1987 (Johnson et al., 1996) was linked to a visit to a bat cave. In 2007, MARV was detected in tissue of Egyptian fruit bat Rousettus aegyptiacus (Kuzmin, et al., 2010). Notably, the MARV strain detected by Kuzmin (2010) differed genetically from ones reported in 1980 and 1987 suggesting circulation of several strains in the same cave and/or its environments. In a more significant case, a close encounter with a bat was reported to have resulted in the death of a tourist due to infection with Duvenhage virus, a bat rabies virus (Van Thiel et al., 2009), presenting evidence of direct virus transmission from bats to humans.
An overarching quest in virus discovery is to detect viruses and also determine transmission pathways especially to humans. Bats are likely to shed some viruses through feces and urine, materials that can be highly infectious to humans. Therefore, this study focused on collecting freshly voided bat feces for molecular detection of multiple virus families. Detection of the virus in these samples, will further strengthen the recommended use of non-invasive methods for virus discovery thereby furthering conservation efforts. In addition, presence of viruses in bat feces reiterates the risk of virus transmission to humans, especially through direct hand contact with freshly voided feces/urine. These factors are important especially in understanding the epidemiology of bat-borne viruses as humans increase their encroachment to wildlife habitats in most African countries including Kenya. In this study, we reveal the virus diversity from bat fecal samples collected in different habitats and in different geographical locations of Kenya.
Study area and sample collection
Fresh bat feces were collected from 48 sites within 17 counties in 5-formerprovincial regions of Kenya (Figure 1) within the period of November 2012 to November 2015. The sampled bat colonies included roosts in caves, trees, inhabited and abandoned buildings (offices and homes) all with relatively high level of human-bat interaction. Fresh fecal samples (One fecal pellet per tube) were collected and placed in RNALater preservative solution (QIAGEN, Germany) as described previously (Ge et al., 2011). Briefly, to collect the fresh fecal pellets, clean polythene sheets measuring 2.0×2.0 m were spread under known bat roosting sites at about 18:00 hours and left to remain until the next morning at about 5:30–6:00 hours the next morning when the pellets were picked into vials that were placed in cooled iceboxes. Samples were transported to the laboratory and stored –80℃ until they were analyzed. Bat species composition in particular roost sites were identified by an experienced chiro- pterologist and further classified by mitochondrion cytochrome b amplification and sequencing.
Sample processing and Virus detection
The samples were completely homogenized by vortexing for 15 seconds and then centrifuged at 12000g for 1 minute to obtain a clear supernatant. Viral nucleic acid was extracted from 200μL of the supernatant using the high pure Viral RNA kit (Roche, Germany) following manufacturer’s instructions. Virus screening was performed by RT-PCR using the Invitrogen OneStep-RT PCR kit using gene specific primers (Supplementary Table S1). Virus families screened included Coronaviridae (CoV), Filoviridae, Reoviridae, Astroviridae (AstV), Adenoviridae (ADV), Polyomaviridae (PYV), Paramyxoviridae (PMV) and Caliciviridae (CalV). The RT-PCR was performed in a 25 μL reaction mix containing 12.5 μL PCR 2 × reaction mix Buffer, 10pmol of each primer, 2.5 mol/L MgSO4, 20 U RNase Inhibitor, 1 μL SuperScript III/ Platinum Taq Enzyme Mix, and 5 μL RNA. A nested/hemi-nested PCR was done using 1uL of the first-round PCR product as a template in a 25 μL reaction mix contained 2.5 μL PCR reaction buffer, 5 pmol of each primer, 50 mmol/L MgCl2, 0.5 mmol/L dNTP and 0.1 μL Platinum Taq Enzyme (Invitrogen). The PCR products were visualized on 1% agarose gel after electrophoresis. The amplicons were gel purified and sequenced directly using the ABI prism sequencer or cloned using pGEM-T Easy Vector System before sequencing.
Sequence and phylogenetic analysis
Sequences were edited using Bio-edit program and preliminary identification done by GenBank Blast searches on the non-redundant database. Sequences were aligned using ClustalW and phylogenetic trees constructed using the maximum likelihood method with an appropriate model and 1000 bootstrap replications. Phylogenetic analysis using MEGA v5.0 (Tamura et al., 2011) software was carried out with representative sequences of the related viral families available in public databases.
A total of 3334 bat fecal samples were collected from various sites in 17 counties in Kenya between 2012 and 2015 (Figure 1) and preserved in RNALater solution. We identified 8 virus families including;Astroviridae, Adenoviridae, Caliciviridae, Coronaviridae, Flaviviridae, Paramyxoviridae, Polyomaviridae and Reoviridae (Table 1). We did not detect filoviruses in any of the samples. The Coastal region of Kenya had the richest diversity since all the 8 virus families were detected. Virus families CoVs, AstVs, ADVs and PMVs were most spatially distributed as they occurred in all the 5-former provincial regions. For the first time in Kenya, AstVs and CalVs were detected in bats, while new genotypes of RVA are reported in R. aegyptiacus and T. mauritianus.
We detected two rotaviruses closely related to group A rotaviruses in two different bat species from two populations that are geographically far from each other, specifically 2980/BatRVA-R. aegyptiacus in Meru and 322/BatRVA-Taphozous mauritianus bat in Kwale region. Phylogenetic analysis based on the partial sequence of VP6 gene revealed that 2980/BatRVA was more closely related to a human RVA strain (G3P) isolated from Kenya in 1997 (Ghosh et al., 2011), while 322/BatRVA was closely related to a human RVA strain (G6P) isolated in Belgium in 1997 (Matthijnssens et al., 2008) (Figure 2). Based on RotaC analysis, the genotypes for the VP6 of 322/BatRVA and 2890/BatRVA were I2 and I16 respectively. On sequence similarity comparison, the partial sequences of 2980/BatRVA and 322/BatRVA were seen to have nucleotide sequence similarities ranging from 69.67 to 87.6% and 69.07 to 94.68% with selected representative sequences of human and animal rotavirus A. The two sequences shared a 79.59% nucleotide similarity to each other and 79.86% and 79.25% nt similarity to the previously reported Eidolon helvum-derived RVA-KE4852/07 from Kenya (Supplementary Table S2).
A maximum likelihood tree of rotaviruses group A based on partial (776bp) VP6 gene.Porcine Rotavirus C (Accession no. KC164679.1) was used as an outgroup. Detected bat rotavirus strains are depicted in boldface and identified with diamond node markers (◆) while other bat-derived rotaviruses are depicted by circular node marker (○). Alignments were performed using ClustalW and bootstrap values were determined using 1000 replicates on Mega5. Detected taxa are named in the following pattern: Virus strain/Identification code/Bat Host species/Country/Specific collection location/Year of collection.
Coronavirus RNA was detected in 3.97% of samples (80 out of 2,014) tested and were present in all the sampled provinces. To determine virus diversity, we carried out a phylogenetic analysis based on a 440 bp fragment of the RNA-dependent RNA-polymerase (RdRp) gene region.Figure 3 shows a maximum likelihood phylogeny of our CoV sequences reconstructed with other previously identified CoVs from animals and humans from the GenBank. The BtCoVs clustered within the alphacoronavirus and betacoronavirus groups. In the alphaco ronavirus group, the viruses detected in this study were mapped into 6 clusters five of which had host restriction and the other one had mixed host species. The cluster with mixed host species included Rhinolophus and Hipposideros bats, two genera that previously were classified within the same family (Teeling et al., 2005). In the betacoronavirus group, two of the clusters were HKU-9 related CoVs while the third formed a branch related to the SARS-L CoVs denoting it as a distant relative of severe acute respiratory syndrome (SARS) CoV.
Novel strains of caliciviruses were detected in two sites; Kilifi and Marsabit that are geographically distant and ecologically different, in two bat genera;Mops species and R. aegyptiacus. The sequences were identified as sapoviruses on the basis of Blastx analysis of the protein database. Two sequences derived from M. condylurus bats (2228/BtCalV/Mops and 2255/BtCalV/Mops) shared 100% amino acid sequence similarity and a 57.77% amino acid sequence similarity with the 22/BtCalV/R. aegyptiacus (Supplementary Table S3). Viruses (2228/ BtCalV/Mops and 2255/BtCalV/Mops were identical sequences originating from the same colony site. These two sequences from Mops spp had 21.11%–64.44% nucleotide sequence identity to other members of the sapovirus group available in the GenBank while sample 22/R. aegyptiacus/BtCalV had 26.66%–55.55% sequence similarity. Upon phylogenetic analysis, the three BtCalVs detected in this study clustered with other reported bat sapoviruses; however, 22/R. aegyp tiacus/BtCalV formed an independent branch, indicating higher variability (Figure 4).
A maximum likelihood tree of caliciviruses based on a partial amino acid sequence. of the RdRp gene (91aa). The Novel bat caliciviruses detected in this study are shown in boldface and marked with a diamond (◆) on the nodes. Detected taxa are named in the following pattern: Virus strain/Identification code/Bat Host species/Country/Specific collection location/Year of collection.
A total of 132 samples tested positive for astroviruses (12.83% positive rate) with two samples from E. helvum bat having co-infection with two variant strains. The sequences revealed a noteworthy viral diversity in the astroviruses present in bats in Kenya. All bat taxa screened in this study were found to harbor astroviruses except Epomophorus species, an observation that may be attributed to low sample numbers for this group. To the best of our knowledge this is the first reporting of bat astroviruses in Kenya as well as first detection of astroviruses in E. helvum, R. aegyptiacus, Cardioderma cor, Coleura afra, Chaerophon pumilis, Triaenops affer and M. condylurus. A phylogenetic analysis of the partial RdRp gene (350 nt) showed that the viral sequences detected in the Kenyan bats clustered together with the other reported bat astroviruses (Supplementary Figure S1) from Asia and Europe (Chu et al., 2008; Drexler et al., 2011; Hu et al., 2014; Kemenesi et al., 2014). The BtAstVs detected in Kenya were found to be very diverse with sequences clustering within several clades of the phylogenetic tree (Supplementary Figure S1). Further sequence comparison showed that the nucleotide identities among the Kenyan-BtAstVs ranged between 25.77 to 83.16%.
Analysis of bat derived adenovirus sequences from this study alongside those available in the GenBank clustered the 26 identified Kenyan-BtADVs within the bat Mastadenovirus genus (Supplementary Figure S2). The sequences occurred in 5 different clusters with cluster 1 comprising entirely of Kenyan bat-ADVs from R. aegyptiacus, E. helvum, C. afra and Otomops martiensenii. Cluster 2 comprised of two sequences from E. helvum that formed a clade with human ADV3/E and chim panzee/Simian ADV. These two sequences were identical (100% nt identities) and had 67.14% to 71.49% nt identities with members of the phylogenetically related taxa. Clusters 3, 4 and 5 were closely related to other reported bat-derived taxa and were obtained from Rhinolophus species, C. cor, Miniopterus spp and C. afra.
We detected 57 paramyxovirus sequences from 5 different bat species, a positive rate of 2.4%. These viruses were found in three of the five regions sampled (Table 1), with the coastal region displaying the highest positive rates. On phylogenetic reconstruction of our sequences with other detected animal and human paramyxoviruses, our viruses clustered into three distinct groupings (Supplementary Figure S3). Cluster 1 comprised of viruses derived from different bat host species and was grouped within the Jeilong-related paramyxoviruses group. Clusters 2 and 3 comprised of viruses detected in E. helvum and fell within the Rubula-related and Henipa-related virus clusters. BtPMV3851, which clustered within the Henipa-related bat-PMV group, shared 97.85% nt identity with a BatPMV detected in the urine of E. helvum bush-meat isolated in DRC in 2009 (Accession number HE647836.1).
Polyomavirus and Flavivirus
Polyomaviruses sequences were detected in four bat species tested;O. martiensenni, E. helvum, Myotis tricolor and Miniopterus minor. The positive rate was 1.65% much lower in comparison to the 11.8% positive rate reported by Tao et al. (2013) in a study of PYVs in Kenyan bats. Result of Blastn and phylogenetic analysis showed that most of the detected viruses were closely related to other reported bat polyomaviruses (Supplementary Figure S4). The virus sequences formed six separate clusters with 5 being closely related to other bat polyomaviruses and one cluster showing a distant relationship with non-human primate derived sequences.
Two flavivirus sequences were detected in E. helvum bats in the coastal region of Kenya. The nucleotide sequences were most closely related to Jugra and Bagaza viruses (88%–89% nt similarities) upon GenBank searches as well as phylogenetic analysis.
Over 100 bat species (Patterson et al., 2012) exist in Kenya and occur in over a wide range of ecological regions. Due to rapid growth of human population and various anthropogenic activities, several incidences of bat colonies roosting in buildings, farms and other human inhabited areas are common. Such sites, as well as other natural bat habitats including as caves also with medium to high human activity, motivated our sampling design. Over the past few decades, about 75% of newly emerging and re-emerging pathogens largely linked to wildlife have been described (Woolhouse, 2002). Numerous viruses which in several cases occur asymptomatically in bats have been found to be of public health and veterinary health concern (Calisher et al., 2006). In several cases such as SARS-CoV, Ebola virus and Nipah virus, spillover events have resulted in disease cases in both livestock and humans (Chua et al., 2000; Liu et al., 2015; Zhong et al., 2003). Spillover from bats to humans has been shown to be direct (Van Thiel et al., 2009) or in some cases linked to availability of intermediate hosts (Wang et al., 2005).
Here we report the detection of eight viral families from bat fecal samples collected across Kenya. We detected viral sequences related to Adenoviridae, Astroviridae, Caliciviridae, Coronaviridae, Flaviviridae, Paramyxoviridae, Polyomaviridae and Reoviridae (specifically Rotavirus genus). The detection rates for the different viral families varied across the country with counties in the coastal region displaying the highest virus diversity (Table 1). The high virus richness in the coastal region may be attributed to a higher diversity in bat species screened. Among the detected sequences there were novel rotaviruses, caliciviruses and astroviruses. Similar to other previous studies carried out to screen these viruses in bats, we report high detection rates for AstVs compared to RVA and CalV (Hu et al., 2014; Kemenesi et al., 2014; Tse et al., 2012).
Rotavirus A sequences have previously been detected in various bat species in a number of countries including Kenya (E. helvum), China (Rhinolophus hipposideros and Aselliscus stoliczkanus), France (Myotis mystacinus) and recently in Cameroon (E. helvum) and Brazil (multiple species) (Asano et al., 2016; Dacheux et al., 2014; Esona et al., 2011; He et al., 2013; Xia et al., 2014; Yinda et al., 2016). Group H rotaviruses which have previously been reported in humans and pigs have also been recently detected in bats in South Korea (Kim et al., 2016). Previous studies, reported that the detected bat RVAs were potentially reassortant strains carrying genes from the human RVA hence indicating a high possibility of interactions between bat and human strains of the virus. Partial sequences of BatRVA VP6 gene amplified in this study were found to be quite closely related to human Rotavirus A. Viruses 322/BatRVA and 2980/ BatRVA detected from T. mauritianus and R. aegyptiacus, respectively, were most closely related to RVA-VP6 from Belgium (EF554119.1) and Kenya (HM 627557.1) which have the genotype constellations G6-P-I2-R2-C2-M2-A3-N2-T6-E2-H3 and G3-P-I16-R8-C5-M5-A5-N5-T5-E13-H5 respectively. The 322/BatRVA was assigned genotype I2 which has been described in several species including but not limited to cow, simian, rabbit as well as human but not in bats. VP6 genotypes described in bats are I15 in Kenya and Cameroon, I8 in China and I22 in Cameroon (Esona et al., 2010; He et al., 2013; Yinda et al., 2016). The 2980/ BatRVA was assigned I16 genotype with 87.4% nt similarity to the Kenyan human B10 isolate, which is notably close to the 85% cut-off point. The I16 genotype has only been previously assigned to the Kenyan human isolate, hence our virus is the second sequence to enter this group. This alludes to reassortment of human and bat virus strains supporting previous research findings. These two viruses were genotypically distinct from each other, as they shared a nucleotide similarity of 79.9% which is lower than the set cut-off values of 85% for VP6 gene (Matthijnssens et al., 2008). Further sequence characterization of the virus strains detected in our study is essential to provide additional information on specific genotype. Our study provides the first report of batRVA in T. mauritianus and R. aegyptiacus. The presence of several strains of rotaviruses in multiple species of bats shows that bats are natural reservoirs of rotavirus. The alleged reassortment with human and other animal rotaviruses could indicate a strong potential of emergence and re-emergence of novel and known infectious rotaviruses.
Diverse strains of coronavirus have been described in bats. Of these strains, SARS-CoV and MERS CoV have caused substantial disease in humans. These two viruses cause acute respiratory distress syndrome (ARDS) with high mortality rates (Hu et al., 2015). SARS-CoV was first described in China in 2002 (Zhong et al., 2003) while MERS-CoV first emerged in Saudi Arabia in 2012 (Zaki et al., 2012). Since then several SARS-CoV and MERS-CoV related viruses have been discovered in other areas of the world. Our results display a high diversity of bat CoVs in the Kenyan population including the detection of 229E-related CoV and a distant relative of the SARS-CoV. The 229E-related CoVs sequences were detected in bats of Hipposideros species from two different sites: Taita Taveta and Meru counties that are far apart. The colonies roost in buildings within wildlife conservation areas, with a high level of human activity. On phylogenetic analysis as well as nucleotide and amino acid sequence similarity comparisons, our 229ECoV-like sequences were more closely related to Alpaca-CoV which was first described in the USA from camelid species (Crossley et al., 2012). The sequences also show a distant relationship to other bat-derived 229E-like CoVs from Ghana (Corman et al., 2015). These finding support the hypothesis of ancestral origin of human 229E virus from bats with a camelid intermediate host.
Bat caliciviruses have previously been reported in Hipposideros pomona bats from China (Tse et al., 2012), Myotis daubentonii, Myotis alcuthoe and Eptesicus serotonus from Hungary (Kemenesi et al., 2014). Here, we found Sapovirus sequences in two different species (Mops condylurus and Rousettus aegyptiacus) which have not been previously reported to harbor caliciviruses. The viruses were detected within the same province (Former Eastern province) but from two sites (over 300KM apart) that vary in ecology, one being a cave and the other a household attic. Sequence comparison and phylogenetic analyses show the detected sequences to be divergent from previously reported BtCalVs. The viruses from Mops species segregated with other reported bat CalVs but were genetic variants with 21.11–64.44% nucleotide sequence similarities with other reported bat sapoviruses. The virus sequence derived from R. aegyptiacus species which was more genetically distant formed an independent branch related to the bat sapovirus cluster. To our knowledge, our results are the first documented report of BatCalVs in Kenya and Africa. Given that we obtained the viral sequences from two species in separate sites, there is a high likelihood that more strains of CalVs are circulating in bat populations.
The prevalence and diversity of astroviruses in the Kenyan bat population was high, recording the highest detection rate and most widely distributed virus species (by location and host species) compared to all others de tected in this study. The diversity was reflected in both within host species and site-origin of the viruses. Phylogenetic analysis based on partial RdRp gene revealed a significantly high diversity among the detected strains of BtAstVs (Supplementary Figure S1). Unlike previous reports (Hu et al., 2014), the strains of viruses present in the Kenyan bats exhibit only partial host restrictions with clades having sequences from various bat species. Sequences of BtAstVs from E. helvum and R. aegyptiacus irrespective of location of origin, display a close relatedness between strains in comparison to sequences derived from other bat species. Further phylogenetic comparison of our BtAstV sequences with other non-bat derived astroviruses reveal a distant relationship of R. aegyptiacus derived sequences to porcine AstVs and E. helvum derived sequences to feline astroviruses, while the other species remain clustered with bat-derived AstVs. Nucleotide identities of R. aegyptiacus to porcine AstVs ranges from 60.13–64.26% and E. helvum to feline AstVs ranges from 51.20–70.44%, values which are higher than similar comparisons to BatAstVs (Data not shown).
A previous study reported BtADVs in Chaerophon and Otomops bat species in Kenya (Conrardy et al., 2014). We detected BtADVs in 3 additional species;M. minor.M. tricolor, E. helvum as well as in O. martiansenii. The BtADVs mapped to 5 clusters, all four of which were within clades containing other BtADVs (Supplementary Figure S2). One of the clusters containing two BtADVs from E. helvum species, infers a distant relationship to human ADVs and simian ADV. This is consistent with results obtained upon nucleotide sequence comparison on NCBI (BAT062/BtADV 89% nt identity to human ADV-E on Blastn and 69% identities to Simian ADV on Blastx; while BAT085/BtADV had 85%nt identity to Simian ADV3 on Blastn and 69% to simian ADV on Blastx). These two E. helvum derived ADV sequences share 95.5% nucleotide identity.
Paramyxoviruses have previously been detected in five bat species in Kenya including;C. cor, Chaerophon pumulis, Miniopterus species, O. martiensseni and R. aegyptiacus (Conrardy et al., 2014). Our results indicated a wider host range with detection of viral sequences in six additional bat species including;E. helvum, T. mauritianus, Hipposideros species, C. afra, M. condylurus and Rhinolophus species. Consistent with the previous study, O. martiensenii displayed the highest detection rate of PMVs. In both studies, Henipa-related viruses were detected in fruit bats; in E. helvum in our study and R. aegyptiacus in the study by Conrardy et al. (2014) . These two separate studies report presence of Henipa-related viruses in fecal swabs of fruit bats in Kenya. The E. helvum colony found positive for this virus roosted within a highly populated market place and may pose a public health risk. Henipaviruses have been reported to cause serious disease in humans and have a high mortality rate of 70% (Aljofan, 2015). Detection of Henipa-related virus in E. helvum and R. aegyptiacus
Polyomavirus and Flavivirus
The PYVs detected in this study were seen to be Orthopolyomaviruses based on Blastn and phylogenetic analysis. Similar to a previous study, further analysis of the partial VP1 gene of the detected PYVs showed relatedness to bat PYVs available in the GenBank, as well as to non-human primate PYVs (Tao et al., 2013). These results infer a historical common ancestor for the polyomaviruses bat and non-human primates.
Interestingly, flavivirus sequences were also detected in the E. helvum colony roosting near a market. These two flavivirus sequences were nearly identical, bearing 93.58% nucleotide sequence identity. The sequences had 83.33% and 87.17% nt sequence identity with their closest characterized flavivirus relative; Jugra virus (DQ859066.1). Jugra virus which has previously been detected in Cynopterus brachiotis (lesser short-nosed fruit bat) is a mosquito transmissible virus and a distant relative of Dengue virus (Calisher et al., 2006). Dengue virus disease is endemic in East Africa, including Kenya (Sutherland et al., 2011). Further characterization of these detected virus strains is important to understand their genome as well as their relationship to other disease-causing flaviviruses such as Dengue virus.
Detection of viruses in bat fecal material is of high significance in terms of public health as it is one of the most direct ways that bat pathogens can be transmitted to humans. For instance, touching bat feces and surfaces contaminated with bat feces is a risk for human infection especially if skin on hands are broken. Fecal droppings from fruit bats may contaminate fruits on trees and vegetables, hence consumption of raw or unwashed fruits and vegetables may be a health risk. Specifically, rotavirus, calicivirus and astroviruses which are transmitted via the oral-fecal route and cause significant health burden in the third world nations. The high similarity of the bat RVA to the human strain of the virus raises genuine concerns of possible transmission considering the cohabitation of humans and many bat species. Our results provide new and extended knowledge with regards to the virus diversity in Kenyan bat species. We recommend sustained surveillance to provide status of most current circulating pathogens and hence inform preparedness in cases of bat-related emerging and re-emerging disease outbreaks. The role of bats in maintenance and circulation of the detected novel viruses (astro, calici- and rota-viruses) is still largely unknown and further studies ought to be carried out to this end.
We are grateful for the Government of Kenya who permitted this study, specifically the Directors of Kenya Wildlife Service and National Museums of Kenya. We also appreciate the field and laboratory teams in Kenya and China who assisted in many different ways towards the success of this work. This work was funded by Sino-Africa Joint Research Center (SAJC201313 and SAJC201605).
COMPLIANCE WITH ETHICS AND GUIDELINES
The authors declare they have no conflicts of interest. The study was approved by the Kenya Wildlife Service, research permit, KWS/BRM/5001. All institutional and National guidelines for care and handling use of animals were followed.
Z-LS, SO and CW designed the study. Z-LS, SO, BA and VO coordinated the study. Z-LS, SO, VO, BA, FG, XYG, XLY, LJW and CW participated and supported in field sampling. CW and AZ carried out molecular studies under the supervision of HB. CW analyzed the data. CW and Z-LS drafted the manuscript. All authors read and approved the final manuscript
This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
- . Aljofan M. 2015. Henipaviruses are not yet history. Future Virol, 10: 507-515.
- . Anthony SJ, Epstein JH, Murray KA, Navarrete-Macias I, Zambrana-Torrelio CM., Solovyov A.Ojeda-Flores R, Arrigo NC, Islam A, Ali-Khan S, Hosseini P, Bogich TL, Olival KJ, Sanchez-Leon MD, Karesh WB, Goldstein T, Luby SP, Morse SS, Mazet JAK, Daszak P, Lipkin WI. 2013. A strategy to estimate unknown viral diversity in mammals. mBio, 4: e00598.
- . Asano KM, Gregori F, Hora AS, Scheffer KC, Fahl WO, Iamamoto K, Mori E, Silva FDF, Taniwaki SA, Brandão PE. 2016. Group A rotavirus in Brazilian bats: description of novel T15 and H15 genotypes. Arch of Virol, 161: 3225.
- . Baker KS, Leggett RM, Bexfield NH, Alston M, Daly G, Todd S, Tachedjian M, Holmes CEG, Crameri S, Wang L-F, Heeney JL, Suu-Ire R, Kellam P, Cunningham AA, Wood JLN, Caccamo M, Murcia PR. 2013. Metagenomic study of the viruses of African straw-coloured fruit bats: Detection of a chiropteran poxvirus and isolation of a novel adenovirus. Virol J, 441: 95-106.
- . Calisher CH, Childs JE, Field HE, Holmes KV, Schountz T. 2006. Bats: Important reservoir hosts of emerging viruses. Clin Microbiol Rev, 19: 531-545.
- . Chu DKW, Poon LLM, Guan Y, Peiris JSM. 2008. Novel astroviruses in insectivorous bats. J Virol, 82: 9107-9114.
- . Chua KB, Bellini WJ, Rota PA, Harcourt BH, Tamin A, Lam SK, Ksiazek TG, Rollin PE, Zaki SR, Shieh WJ, Goldsmith CS, Gubler DJ, Roehrig JT, Eaton B, Gould R, Olson J, Field H, Daniels P, Ling AE, Peters CJ, Anderson LJ, Mahy BW. 2000. Nipah virus: a recently emergent deadly paramyxovirus. Science, 288: 1432-1435.
- . Conrardy C, Tao Y, Kuzmin IV, Niezgoda M, Agwanda B, Breiman RF, Anderson LJ, Rupprecht CE, Tong S. 2014. Molecular detection of adenoviruses, rhabdoviruses, and para myxoviruses in bats from Kenya. Am J Trop Med. Hyg, 91: 258-266.
- . Corman VM, Baldwin HJ, Tateno AF, Zerbinati RM, Annan A, Owusu M, Nkrumah EE, Maganga GD, Oppong S, Adu-Sarkodie Y, Vallo P, Ribeiro LV, Filho F, Leroy EM, Thiel V, Hoek L, Poon LLM, Tschapka M, Drexler JF. 2015. Evidence for an ancestral association of human coronavirus 229E with bats. J Virol, 89: 11858-11870.
- . Crossley BM, Mock RE, Callison SA, & Hietala SK. 2012. Identification and characterization of a novel alpaca respiratory coronavirus most closely related to the human coronavirus 229E. Viruses, 4: 3689-3700.
- . Dacheux L, Cervantes-Gonzalez M, Guigon G, Thiberge JM, Vandenbogaert M, Maufrais C, Caro V, Bourhy H. 2014. A preliminary study of viral metagenomics of French bat species in contact with humans: Identification of new mammalian viruses. PLoS ONE, 9: e87194.
- . Drexler JF, Corman VM, Gloza-Rausch F, Seebens A, Annan A, Ipsen A, Krupa T, Muller MA, Kalko EKV, Adu-Sarkodie Y, Oppong S, Drosten C. 2009. Henipavirus RNA in African bats. PloS One, 4: e6367.
- . Drexler JF, Corman VM, Wegner T, Tateno AF, Zerbinati RM, Gloza-Rausch F, Seebens A, Muller MA, Drosten C. 2011. Amplification of emerging viruses in a bat colony. Emerg Infect Dis, 17: 449-456.
- . Esona MD, Banyai K, Foytich K, Freeman M, Mijatovic-Rustempasic S, Hull J, Kerin T, Steele AD, Armahd GE, Geyer A, Page N, Agbaya VA, Forbi JC, Aminui M, Gautama R, Seheri LM, Nyangao J, Glass R, Bowena MD, Gentsch JR. 2011. Genomic characterization of human rotavirus G10 strains from the African Rotavirus Network: Relationship to animal rotaviruses. Infect Genet Evol, 11: 237-241.
- . Esona MD, Mijatovic-Rustempasic S, Conrardy C, Tong S, Kuzmin IV, Agwanda B, Breiman RF, Manyai K, Niezgoda M, Rupprecht E, Gentsch JR, Bowen MD. 2010. Reassortant group a rotavirus from straw-colored fruit bat (Eidolon helvum) . Emerg Infect Dis, 16: 1844-1852.
- . Gatherer D. 2014. The 2014 Ebola virus disease outbreak in West Africa. J Gen Virol, 95: 1619-1624.
- . Ge X, Li J, Peng C, Wu L, Yang X, Wu Y, Zhang Y, Shi Z. 2011. Genetic diversity of novel circular ssDNA viruses in bats in China. J Gen Virol, 92: 2646-2653.
- . Ghosh S, Gatheru Z, Nyangao J, Adachi N, Urushibara N, Kobayashi N. 2011. Full genomic analysis of a simian SA11-like G3P. Infect Genet Evol, 11: 57-63.
- . He B, Yang F, Yang W, Zhang Y, Feng Y, Zhou J, Xie J, Feng Y, Bao X, Guo H, Li Y, Kia L, Li N, Matthijnssens J, Zhang H, Tu C. 2013. Characterization of a Novel G3P. J Virol, 87: 12357-12366.
- . Hu B, Chmura AA, Li J, Zhu G, Desmond JS, Zhang Y, Zhang W, Epstein JH, Daszak P, Shi Z. 2014. Detection of diverse novel astroviruses from small mammals in China. J Gen Virol, 95: 2442-2449.
- . Hu B, Ge X, Wang LF, Shi Z. 2015. Bat origin of human co ronaviruses. Virol J, 12: 221.
- . Johnson E, Johnson B, Silverstein D, Tukei P, Geisbert T, Sanchez A, Jahrling P. 1996. Characterization of a new Marburg virus isolated from a 1987 fatal case in Kenya. Arch Virol, 11: 101-114.
- . Jones KE, Patel NG, Levy M. A, Storeygard A, Balk D, Gittleman JL, Daszak P. 2008. Global trends in emerging infectious diseases. Nature, 451: 990-993.
- . Kading RC, Gilbert AT, Mossel EC, Crabtree MB, Kuzmin IV, Niezgoda M, Agwanda B, Markotter W, Weil MR, Montgomery JM, Rupprecht CE, Miller BR. 2013. Isolation and molecular characterization of Fikirini rhabdovirus, a novel virus from a Kenyan bat. J Gen Virol, 94: 2393-2398.
- . Kemenesi G, Dallos B, Gorfol T, Boldogh S, Estok P, Kurucz K, Kutas A, Foldes F, Oldal M, Nemeth V, Martella V, Banyai K, Jakab F. 2014. Molecular survey of RNA viruses in Hungarian bats: discovering novel astroviruses, coronaviruses, and ca liciviruses. Vector-Borne Zoonot Dis, 14: 846-855.
- . Kim HK, Yoon SW, Kim DJ, Koo BS, Noh JY, Kim JH, Choi YG, Na W, Chang KT, SOng D, Jeong DG. 2016. Detection of Severe Acute Respi ratory Syndrome-Like, Middle East Respiratory Syndrome-Like bat coronaviruses and group H Rotavirus in faeces of Korean bats. Transbound Emerg Dis, 63: 365-372.
- . Kuzmin IV, Mayer AE, Niezgoda M, Markotter W, Agwanda B, Breiman RF, Rupprecht CE. 2010a. Shimoni bat virus, a new representative of the Lyssavirus genus. Virus Res, 149:197-210.
- . Kuzmin IV, Niezgoda M, Franka R, Agwanda B, Markotter W, Breiman RF, Shieh W-J, Zaki SR, Rupprecht CE. 2010b. Marburg virus in fruit bat, Kenya. Emerg Infect Dis, 16: 352-354.
- . Leroy EM, Kumulungui B, Pourrut X, Rouquet P, Hassanin A, Yaba P, Delicat A, Paweska JT, Gonzalez JP, Swanepoel R. 2005. Fruit bats as reservoirs of Ebola virus. Nature, 438: 575-576.
- . Liu WB, Li ZX, Du Y, Wen G. 2015. Ebola virus disease: from epidemiology to prophylaxis. Mil Med Res, 2: 7.
- . Matthijnssens J, Ciarlet M, Heiman E, Arijs I, Delbeke T, McDonald SM, Palombo EA, Iturriza-Gomara M, Maes P, Patton JT, Rahman M, Van Ranst M. 2008a. Full genome-based classification of rotaviruses reveals a common origin between human Wa-Like and porcine rotavirus strains and human DS-1-like and bovine rotavirus strains. J Virol, 82: 3204-3219.
- . Matthijnssens J, Ciarlet M, Rahman M, Attoui H, Bányai K, Estes MK, Gentsch JR, Itturiza-Gomara M, Kirkwood C, Martella V, Merens PPC, Nakagomi O, Patton JT, Ruggeri FM, Saif LJ, Santos N, Steyer A, Taniguchi K, Desselberger U, Van Ranst M. 2008b. Recommendations for the classification of group A rotaviruses using all 11 genomic RNA segments. Arch Virol, 153: 1621-1629.
- . Moratelli R, Calisher CH. 2015. Bats and zoonotic viruses: Can we confidently link bats with emerging deadly viruses? Mem. I. Oswaldo-Cruz, 110: 1-22.
- . O’Shea TJ, Cryan PM, Cunningham AA, Fooks AR, Hayman DTS, Luis AD, Peel AJ, Plowright RK, Wood JLN. 2014. Bat flight and zoonotic viruses. Emerg Infect Dis, 20: 741-745.
- . Olival . 2016. To cull, or not to cull, bat is the question. EcoHealth, 13: 6-8.
- . Patterson BD, Webala PW. 2012. Keys to the bats (Mammalia: Chiroptera) of East Africa. Fieldiana Life and Earth Sciences, 12: 6.
- . Plowright RK, Eby P, Hudson PJ, Smith IL, Westcott D, Bryden WL, Martin G, Tabor GM, Skerratt LF, Anderson DL, Crameri G, Quammen D, Jordan D, Freeman P, Wand LF, Epstein JH, Marsh GA, Kung NY, McCallum H. 2015. Ecological dynamics of emerging bat virus spillover. Proc Biol Sci, 282: 2014 -2124.
- . Shi Z. 2013. Emerging infectious diseases associated with bat viruses. Sci China Ser C, 56: 678-682.
- . Smith DH, Isaacson M, Johnson KM, Bagshawe A, Johnson BK, Swanapoel R, Johnson KM, Killey M, Bagshawe A, Siongok T, Keruga WK. 1982. Marburg-Virus Disease in Kenya. Lancet, 319: 816-820.
- . Sutherland LJ, Cash AA, Huang YJS, Sang RC, Malhotra I, Moormann AM, King CL, Weaver SC, King CH, LaBeaud AD. 2011. Serologic evidence of arboviral infections among humans in Kenya. Am J Trop Med Hyg, 85: 158-161.
- . Tamura K, Peterson D, Peterson N, Stecher G, Nei M & Kumar S. 2011. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol, 28: 2731-2739.
- . Tao Y, Shi M, Conrardy C, Kuzmin IV, Recuenco S, Agwanda B, Alvares DA, Ellison JA, Gilbert AT, Moran D, Niezgoda M, Lindblade KA, Holmes EC, Brieman RF, Rupprecht CE, Tong S. 2013. Discovery of diverse polyomaviruses in bats and the evolutionary history of the Polyomaviridae. J Gen Virol, 94: 738-748.
- . Teeling EC, Springer MS, Madsen O, Bates P, O’Brien SJ, Murphy WJ.. 2005. A molecular phylogeny for bats illuminates bioge ography and the fossil record. Science, 307: 580-584.
- . Tong S, Conrardy C, Ruone S, Kuzmin IV, Guo X, Tao Y, Niezgoda M, Haynes L, Agwanda B, Breiman RF, Anderson LJ, Rupprecht CE. 2009. Detection of novel SARS-like and other coronaviruses in bats from Kenya. Emerg Infect Dis, 15: 482-485.
- . Tse H, Chan WM, Li KSM, Lau SKP, Woo PCY & Yuen KY. 2012. Discovery and genomic characterization of a novel bat sapovirus with unusual genomic features and phylogenetic position. PLoS ONE, 7: e34987.
- . Van Thiel PPAM, De Bie RMA, Eftimov F, Tepaske R, Zaaijer HL, Van Doornum GJJ, Schutten M, Osterhaus ADM, Majoie CBL, Aronica E, Fehlner-Gardiner C, Wandeler AI, Kager PA. 2009. Fatal human rabies due to duvenhage virus from a bat in Kenya: Failure of treatment with coma-induction, ketamine, and antiviral drugs. PLoS Neglect Trop. D, 3: 1-8.
- . Wang M, Yan M, Xu H, Liang W, Kan B, Zheng B, Chen H, Zheng H, Xu Y, Zhang E, Wang H, Ye J, Li G, Li M, Cui Z, Liu YF, Guo RT, Liu XN, Zhan LH, Zhou DH, Zhao A, Hai R, Yu D, Guan Y, Xu J. 2005. SARS-CoV infection in a restaurant from palm civet. Emerg Infect Dis, 11: 1860-1865.
- . Woolhouse MEJ. 2002. Population biology of emerging and re-emerging pathogens. Trends Microbiol, 10: 3-7.
- . Xia L, Fan Q, He B, Xu L, Zhang F, Hu T, Wang Y, Li N, Qiu W, Zheng Y, Matthijnssens J, Tu C. 2014. The complete genome sequence of a G3P Chinese bat rotavirus suggests multiple bat rotavirus inter-host species transmission events. Infect Genet Evol, 28: 1-4.
- . Yinda CK, Zeller M, Maes P, Deboutte W, Beller L, Heylen E, Ghogomu SM, Ranst MV, Matthijnssens J. 2016. Evidence for reassortment of highly divergent novel rotaviruses from bats in Cameroon, without evidence for human interspecies transmissions. Sci Rep, 6: 34209.
- . Zaki AM, Van Boheemen S, Bestebroer TM, Osterhaus AD, Fouchier RA. 2012. Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia. New Eng J Med, 367: 1814-1820.
- . Zhong NS, Zheng BJ, Li YM, Poon LM, Xie ZH, Chan KH, Li PH, Tan SY, Chang Q, Xie JP, Liu XQ, XU J, Li DX, Yuen KY, Peiris JSM, Guan Y. 2003. Epidemiology and cause of severe acute respiratory syndrome (SARS) in Guangdong, People’s Republic of China, in February, 2003. Lancet, 362: 1353-1358.