Abstract
The straw-coloured fruit bat, Eidolon helvum, is Africa’s most widely distributed and commonly hunted fruit bat, often living in close proximity to human populations. This species has been identified as a reservoir of potentially zoonotic viruses, but uncertainties remain regarding viral transmission dynamics and mechanisms of persistence. Here we combine genetic and serological analyses of populations across Africa, to determine the extent of epidemiological connectivity among E. helvum populations. Multiple markers reveal panmixia across the continental range, at a greater geographical scale than previously recorded for any other mammal, whereas populations on remote islands were genetically distinct. Multiple serological assays reveal antibodies to henipaviruses and Lagos bat virus in all locations, including small isolated island populations, indicating that factors other than population size and connectivity may be responsible for viral persistence. Our findings have potentially important public health implications, and highlight a need to avoid disturbances that may precipitate viral spillover.
Similar content being viewed by others
Introduction
Recent studies have demonstrated the potential of bats to act as reservoirs of zoonotic pathogens (as reviewed in Hayman et al.1). One example is the common and conspicuous straw-coloured fruit bat (Eidolon helvum), which has been identified as a reservoir host for Lagos bat virus (LBV, family Rhabdoviridae, genus Lyssavirus)2 and henipaviruses (family Paramyxoviridae)3 in mainland Africa. E. helvum is a gregarious, predominantly tree-roosting species and large roosts (sometimes numbering more than one million bats) frequently exist in close proximity to large human settlements, including Accra (Ghana), Abidjan (Côte d'Ivoire), Dar es Salaam (Tanzania), Lagos (Nigeria) and Kampala (Uganda)4.
Much of the serological evidence for zoonotic pathogens in bats comes from single cross-sectional studies, with few conducted longitudinally or across a representative proportion of the entire species range. However, longitudinal surveys of E. helvum colonies in Ghana have demonstrated relatively high roost-level seroprevalences to LBV over multiple years, which increase with bat age5. These findings indicate endemic circulation with horizontal transmission, making E. helvum a true reservoir host of LBV in that country. Moreover, neutralizing antibodies to LBV have also been detected in cross-sectional serological surveys in Kenya6 and Nigeria7 and LBV has been isolated from a small number of sick or dead wild E. helvum bats in Nigeria, Senegal and Kenya (as reviewed in Banyard et al.2).
Old World fruit bats (Pteropodidae) are the principal reservoir hosts of henipaviruses8, with flying fox populations (Pteropus spp.) found to harbour Nipah virus (NiV) in Southeast Asia, and both Hendra virus (HeV) and Cedar virus (CedPV) in Australia. NiV and HeV are highly pathogenic in humans and other mammals, yet the recently discovered CedPV differs in its apparent apathogenicity in laboratory animal species9. Cross-neutralizing antibodies to HeV and NiV have been detected in sympatric Pteropus spp. and Madagascan fruit bats (E. dupreanum)10, and Hayman et al.3 first documented antibodies to henipaviruses in bats outside the range of Pteropus spp, with a 40% seroprevalence being found in E. helvum in Ghana. These serological findings were recently supported by the detection of henipavirus-like RNA in E. helvum in Ghana and central Africa11,12,13; yet, although a full genome sequence for one of these African henipavirus-like viruses was obtained13, live viruses have not yet been isolated.
These findings collectively highlight the potential for zoonotic pathogen spillover from E. helvum to humans, with routes of infection being via urine12, faeces13 or the hunting and preparation of bat meat for food14. However, no such spillovers have been reported for LBV or African henipaviruses. This might be because spillover has not yet occurred, or it might reflect poor medical surveillance capabilities in much of Africa, and the lack of availability of specific diagnostic assays15.
Much is yet to be understood regarding the host response to natural lyssavirus and henipavirus infections in bats; experimental inoculations have yielded inconsistent results across individuals and studies. Bats infected with lyssaviruses may or may not develop clinical signs corresponding to those seen in other mammals (as reviewed in Banyard et al.2), whereas no clinical illness has been observed in bats infected with henipaviruses8. Acute antibody responses have been observed for both viruses after experimental infection, with boosted titres upon reinfection8,16. An assumption could follow that these infections are immunizing in bats, however, seroconversion is not universally observed, and therefore this remains open to challenge. Typically, pathogens causing acute immunizing infections require large host population sizes and a ‘critical community size’ (CCS) for persistence is expected unless birth rates are very high.
Many uncertainties also remain regarding the specific viral transmission dynamics in E. helvum. Key aspects of this species’ ecology might further increase potential for viral persistence within populations. In particular, it is a migratory species that comprises both permanent and seasonal colonies across much of sub-Saharan Africa4 and a small number of offshore islands, including those in the Gulf of Guinea17 (Fig. 1). However, the widespread and continuous distribution represented in Fig. 1 over-simplifies a more intricate distribution pattern, comprising aggregated populations across a connected, rather than continuous, landscape18. Annual seasonal migrations result in abrupt fluctuations in the size of permanent colonies, and also in the formation of solely seasonal colonies. For example, the largest known E. helvum colony in Kasanka National Park in Central Zambia is populated rapidly each year to reach an estimated 1.5 million individuals19, and persists for just 2 ½ months. Satellite telemetry studies indicate that these bats are capable of migrating vast distances (for example, up to 370 km in one night and ~2,500 km over 5 months)20. It has been suggested that migration occurs along a ‘north–south’ axis, with seasonal movements following latitudinal shifts of the Intertropical Convergence Zone weather system20,21; however, the routes and drivers of migrations are not fully understood. Such large-scale movements are expected to lead to widespread gene flow, and it has been argued that extensive genetic mixing among wildlife populations may increase the potential for viral epidemics22. Therefore, to characterize viral infection dynamics in wildlife populations, information on host population structure and connectivity is needed. Indeed, Plowright et al.23 suggested that a large, weakly coupled asynchronous metapopulation structure might be necessary for population-level persistence of HeV with either acute ‘explosive’, or slow ‘smouldering’ epidemics resulting from spatial heterogeneity in population herd immunity. We recently demonstrated evidence of exposure to henipaviruses in the small, isolated population of E. helvum on the Gulf of Guinea island of Annobón, indicating that a metapopulation model may not be required for persistence of all henipaviruses24. The persistence of lyssaviruses in some temperate insectivorous bat species has been shown to depend on certain life history traits, including hibernation and birth pulses25, but persistence mechanisms in non-hibernating species, such as E. helvum, are unknown.
To determine the extent of genetic and epidemiological connectivity among E. helvum populations, and thus gain better understanding of viral transmission dynamics and zoonotic risk, here we combine genetic and serological analyses of populations across Africa. We use mitochondrial (mtDNA) and nuclear DNA analyses to characterize the range-wide metapopulation structure of E. helvum, and hypothesize that this would inform our understanding of viral dynamics across the population. Together with serological analyses, we assess the epidemiological consequences of this structure for the species’ ability to act as a reservoir host of the potentially zoonotic viruses, LBV and henipaviruses.
Results
Sampling
Samples (including wing membrane biopsies, blood and urine) were obtained from 2,084 individual E. helvum bats across continental Africa and the Gulf of Guinea islands. In addition, pooled urine samples were collected from beneath some colonies. Details of sampling locations (Fig. 1 and Supplementary Data 1) and sample sizes for genetic, serological and urine analyses (Table 1) are provided.
Microsatellite and mtDNA genetic analyses
Overall, results from multiple analyses presented below showed that E. helvum forms a panmictic population across its continental range, with no evidence of isolation by distance (IBD) or structuring according to migratory routes. Bats on the offshore island of Bioko were found to be part of this panmictic population; however, the more isolated island populations in the Gulf of Guinea were genetically distinct from one another and from the continental population.
Of 114 unique cytochrome b (cytb) haplotypes identified from 544 individuals, 75% were singletons (only found in a single individual across all populations, Table 2). Haplotype diversity, molecular diversity, allelic richness and observed heterozygosity were all higher within continental with Bioko (CB) populations than in isolated island (iIS) populations. Nucleotide diversity was low across all populations, but particularly so in Príncipe and Annobón.
Structure among populations assessed by pairwise FST (using microsatellite data) and φST (using mtDNA data) values gave similar results, with near-zero, non-significant values among CB populations, contrasting with larger, significant values between iIS and CB populations (Supplementary Table S1). Each island population was also significantly differentiated from one another. These results were supported by analysis of molecular variance (AMOVA), where maximal structure among groups (high FCT and φCT values) and minimal structure among populations within groups (low FSC and φSC), were observed when populations were separated into three (CB, São Tomé with Príncipe (STP), Annobón) or four (CB, São Tomé, Príncipe, Annobón) groups (Table 3, analyses 7 and 8). IBD analyses detected no positive correlation between genetic distance (Slatkin’s linearized φST and FST) and log geographical distance in any mtDNA or microsatellite analyses (Fig. 2). This finding was consistent when latitude was ignored and longitudinal distances were used in the analyses, accounting for presumed north–south migration routes of E. helvum21.
A Bayesian phylogeny (Supplementary Fig. S1) and median joining haplotype network (Fig. 3) both recovered three main E. helvum clades. The star-like network was characterized by a few common haplotypes, surrounded by many haplotypes present in only 1–5 individuals. Thorough spatial mixing was evident, with the central haplotype (Hap2) being shared by 85 bats representing all CB populations plus a single bat from Annobón. Most bats from the isolated island (iIS) populations (253/272; 93%) were divided between two haplotypes at opposite ends of the network (Hap8, predominantly Annobón, and Hap111, predominantly STP; Supplementary Fig. S2).
Consistent with these results, Bayesian clustering of individual genotypes revealed three clusters (K=3) based on mean likelihood (log P (X|K) values (Fig. 4)), corresponding to populations from CB, STP and Annobón. With increasing values of K, the STP and Annobón clusters remained unchanged, and the CB cluster became increasingly subdivided into multiple clusters of approximately equal proportion (Fig. 3), again indicative of a strong signature of a single panmictic CB population. Analyses run with CB or iIS samples as separate data sets did not reveal additional clusters. Using these three clusters as prior population information to identify potential migrants among clusters, STRUCTURE assignment tests (admixture analyses based on nuclear data), indicated that 19/502 individuals were ‘admixed’ (that is, had an assignment probability (p) to any one main cluster of 0.8>p>0.2). No bats were classified as recent (first generation) migrants (Supplementary Table S2).
Isolation-with-migration models and approximate Bayesian computation were unsuccessful in obtaining reliable estimates of gene flow between these islands, as a result of lack of convergence or unrealistically large estimates of effective population size, respectively.
LBV serological analyses
Using modified fluorescent antibody virus neutralization assays, neutralizing antibodies to LBV were detected in all continental and island locations (Table 1), yet seroprevalences showed significant variation by geographical location. A strikingly low LBV seroprevalence relative to other locations was observed in the Annobón population (χ2=66.5, P<0.001), but seroprevalences in Bioko, São Tomé and Príncipe were not significantly different from mainland populations. Excluding Annobón and populations with sample sizes that were insufficient to allow a reliable seroprevalence to be calculated (Malawi, Zambia and Uganda; n=12, 9 and 4, respectively), the mean LBV seroprevalence was 34% (95% CI: 32−37%) and the range of adult seroprevalences was 24–51% (Supplementary Data 1). In the Annobón population, neutralizing antibodies to LBV were detected in 1 of 72 (1.4%, 0.0–7.5%) bats sampled in 2010 (ref. 24), and in 6 of 49 (14%, 7–27%) bats sampled in 2011.
Henipavirus serological analyses
Antibodies binding to NiV soluble G glycoproteins were detected using Luminex microsphere binding assays in all populations sampled (Table 1). In contrast to the LBV results, henipavirus seroprevalences in all Gulf of Guinea islands (including Annobón) were similar to those in continental populations. Excluding populations with very small sample sizes, as above, the mean henipavirus seroprevalence was 42% (39−44%), with adult seroprevalences ranging from 29−60% (Supplementary Data 1). Using virus neutralization tests (VNTs), a NiV seroprevalence of 5% (11/222, 3−9%) was detected in bats sampled from Tanzania and 1.7% (2/118, 0.5−6%) in bats from Annobón. For bats from Bioko, São Tomé and Príncipe, NiV VNTs were performed on a subset of the samples (those with binding assay median fluorescence intensities >750 (n=49, 20 and 39, respectively)), of which 32, 50 and 51% were neutralizing, respectively.
For both LBV and henipaviruses, no significant differences in seroprevalence were detected between males and females.
Urine analyses
Polymerase chain reactions (PCRs) performed on E. helvum urine samples from Tanzania, Uganda, Malawi, Zambia and Annobón detected paramyxovirus polymerase gene sequences in 3/23 extraction pools (from Ugandan and Tanzanian sampling sites, Table 1). These showed close relationships with sequences detected previously in E. helvum in Ghana12 (Fig. 5). One PCR-positive pooled sample from Tanzania comprises urine expressed directly from the bladders of six individual E. helvum, all of which were seronegative for henipaviruses using microsphere binding assays and VNTs.
Discussion
In this study, using data from both mtDNA and microsatellite markers, we demonstrate that the population of E. helvum is panmictic across its continental African range. An absence of IBD indicated that gene flow was no more likely to occur among neighbouring populations than distant populations of >4,500 km, making E. helvum the largest reported panmictic unit of any mammal, and one of the largest of any vertebrate, exceeded only by the bigeye tuna (Thunnus obesus; >8,000 km)26,27 and the Kentish plover (Chadadrius alexandrines; >10,000 km)28. Even present day human populations retain genetic structure over such large distances29. In fact, the range of E. helvum extends further north and west of the sampling sites in this study, so additional sampling is required to assess whether panmixia extends across this range; a distance of >6,500 km.
The hypothesis that greater genetic differentiation might exist across migratory pathways (on an east–west axis) than along migratory pathways (on a north–south axis) was not supported by our results, probably either because gene flow between distinct migratory populations homogenises allele frequencies, or because E. helvum migration is opportunistic and tracks changes in available food resources rather than following defined migratory routes.
Included in the panmictic E. helvum population are bats on the near-shore island of Bioko in the Gulf of Guinea (which separated from the African continent ~7,000 years ago). Our results indicate that the 32 km stretch of ocean that separates Bioko from the continent is not a significant barrier to dispersal, as might be expected given that individuals are capable of covering such distances during foraging bouts20.
In contrast to the panmictic continental and Bioko (CB) population, populations on the three more isolated Gulf of Guinea islands (São Tomé, Príncipe and Annobón) showed evidence of genetic isolation. This accords with the results from studies of other Gulf of Guinea island taxa, including other species of bat30, bird (for example, Melo et al.31), and reptile (for example, Jesus et al.32). Although E. helvum is a long-range migrant and has been observed as a vagrant on islands 570 km from the African coastline33, the strong genetic structure detected among the island and CB population clusters and the absence of genetic evidence (using assignment tests) of recent migrants between these clusters, indicate that dispersal between clusters (with successful mating) is likely to be rare: no first generation migrants were detected, although some individuals may have been second or third generation migrants. Additional support for population genetic results in this and other fruit bat species comes from genetic studies of external parasites and their pathogens34,35, including detection of congruence between population genetic structure of external parasites and their hosts.
The Gulf of Guinea ocean channels are likely to have provided a barrier to initial colonization and inter-island dispersal. Also, we found that the smallest discrete population of E. helvum (on the island of Annobón) showed genetic divergence and is truly isolated. Our mtDNA and microsatellite results are consistent with those of a previous study that found that E. helvum on Annobón showed differences in morphological traits and allozyme frequencies compared with other islands17. However, while Juste et al.17 also concluded that a lack of phenetic differentiation on Bioko, São Tomé and Príncipe suggested gene flow among the islands, our use of multiple nuclear and mtDNA markers, provides further insight. For example, while São Tomé shows greater connectivity with Príncipe (<150 km apart) than with either the CB or the Annobón populations (which lie >220 km away), the distance separating these two populations is still a substantial barrier to inter-island gene flow, as shown by significant pairwise φST and FST values. As São Tomé and Príncipe are within the same cluster, it is not possible to identify migrants between these two islands using assignment tests. Other genetic methods to estimate gene flow and demographic history between multiple populations, including isolation-with-migration models and approximate Bayesian computation, were unsuccessful in obtaining reliable and credible estimates of gene flow between these islands, suggesting that even for our substantial data sets, modelling of low rates of gene flow using current techniques and assumptions is not robust.
Although genetic analyses cannot replace direct studies on individual bat movements and demographic connectivity, they can contribute to a broader perspective upon which epidemiological studies on transmission and maintenance of viruses among and within populations can be based36. The strong genetic clustering observed here makes it likely that the separation of E. helvum into three distinct genetic population clusters (CB, STP and AN) is echoed as at least three epidemiologically distinct populations. A freely mixing, panmictic continental population would likely facilitate viral transmission among E. helvum colonies across this range. Our serology results are consistent with this, with henipavirus and LBV antibodies being detected across all continental sampling sites at seroprevalences similar to those previously observed for henipaviruses in Ghana3 and for LBV in Ghana, Kenya and Nigeria5,6,7,37.
Further support for the conclusion that distant continental populations may belong to a single epidemiological unit was provided by high nucleotide sequence identities between paramyxoviral sequences detected in E. helvum urine samples from Uganda and Tanzania and those already reported from Ghana12. In that study and others13,38, a diverse range of paramyxovirus sequences, including henipavirus-like sequences, were detected within single E. helvum populations. Further sampling efforts to enable exploration of viral sequence diversity across all the sites studied here would help determine whether different virus variants are maintained by each of these distinct epidemiological units and whether viral diversity may have a role in within-population viral persistence. Additional data are required to fully understand how virus variants are maintained within E. helvum populations.
Although genetic differentiation and isolation of E. helvum in the STP cluster was expected to be reflected epidemiologically, perhaps with an absence of antibodies on these islands due to restricted population sizes, we found that seroprevalences to both viruses were comparable with those on the mainland. These data suggest that: population sizes on each island are sufficient to maintain LBV and henipaviruses and are above the CCS required for persistence (although this concept requires further theoretical exploration for animal populations where birth rates, and hence population sizes, are highly seasonal); sufficient movement may occur between the two islands to maintain a larger epidemiologically connected population; alternative hosts may be involved; or our original assumptions on transmission and persistence may need re-examination (see below). The use of satellite telemetry has been enlightening in other fruit bat species39,40 and would be required to definitively assess movement patterns of bats on these two islands. However, dispersal between São, Tomé and Príncipe was suggested by our observation of a single asynchronous birth on Príncipe in the absence of other pregnant or lactating females, but which was contemporaneous with the presence of neonates on São, Tomé. The asynchronous Príncipe birth is highly unusual for a species, which employs delayed implantation to facilitate a highly synchronized birth pulse. If the two populations are connected via dispersal, the asynchrony in reproductive seasons between São Tomé and Príncipe could facilitate viral persistence by staggered introduction of susceptible individuals, via birth, into the population. Finally, LBV has been detected in bats of several species in Africa, with ranges overlapping that of E. helvum2, but the role that interspecies transmission has in the maintenance of LBV in its host populations remains a gap in our knowledge. Of all these species, only Rousettus aegyptiacus, the Egyptian fruit bat, is present on São, Tomé and Príncipe. This is a cave-roosting species, and mixed colonies with E. helvum are unlikely, although these two species might mix at feeding sites. LBV has been isolated from R. aegyptiacus on two occasions (as reviewed in Kuzmin et al.41), and seroprevalence levels comparable with those reported in E. helvum were detected in Kenya6. On São, Tomé and Príncipe, R. aegyptiacus, or indeed other species, may facilitate the persistence of LBV in E. helvum.
Although findings from the CB and STP populations could be consistent with a metapopulation model of persistence, as proposed for HeV in Australia23 and NiV in Malaysia42, our results from Annobón indicate that this appears unnecessary for the persistence of henipaviruses or LBV in E. helvum. On Annobón, E. helvum is the only bat species confirmed to be currently present and has a population size of only ~2,500 (ref. 24). Surprisingly, and in contrast to findings in other, less-isolated island systems42, the henipavirus seroprevalence in the Annobón E. helvum population was within the range of that observed in both the CB and the STP populations. Conversely, the Annobón LBV seroprevalence was much lower than in other populations. Although evidence of infection with lyssaviruses has been reported in other island bat species (for example, refs 43, 44), the bat populations in those studies were either much larger, within flight distance of continental bat populations, and/or hosted multiple sympatric bat species. In Annobón, all LBV seropositive individuals were adult, and further longitudinal studies are required to determine whether LBV is persistently maintained on this island (that is, the population size is greater than the CCS), or whether these findings represent a single epidemic wave subsequent to introduction of the virus from another population. Unfortunately, deriving a quantitative estimate for the CCS is problematic, particularly for virus-host systems where little information is available regarding host demographics, virus transmission mechanisms and within-host immune responses45. For both LBV and henipaviruses, important areas of future study include viral diversity and phylogeography, within-host persistence and immunity, incubation periods and frequency- versus density-dependent transmission.
Multiple henipavirus-like sequences have been previously reported in E. helvum11,12,13. In the absence of isolation or full genomic characterization, it cannot be definitely confirmed whether these sequences represent true henipaviruses. However, a phylogenetic analysis undertaken here, incorporating the most recently isolated henipavirus (CedPV in Australia) and sequence fragments from bat paramyxoviruses worldwide (Fig. 6) demonstrate that two virus sequences from E. helvum in Gabon13 fall within the clade of currently identified henipaviruses. These sequences therefore likely represent true African henipaviruses.
This study took a multidisciplinary approach, combining ecological, genetic and serological studies, to explore the ways in which the structure, dynamics and connectivity of E. helvum populations across Africa affects the viral transmission dynamics within them. These critical population-level processes are expected to be important in determining viral persistence within populations, and yet, although the three genetically distinct populations identified here are also highly likely to be separated epidemiologically, each of these population clusters is capable of maintaining henipaviruses and LBV, apparently without the need for a metapopulation model of persistence via migration and reintroduction.
The findings presented here have potentially important implications for public health. The large population sizes of E. helvum, its tendency to roost and feed in close proximity to human populations, its extensive distribution across Africa and its frequent harvesting for bushmeat, present numerous opportunities for the exposure of people to excreta, tissues and body fluids from these bats. The widespread presence of potentially zoonotic viruses in this species across Africa might therefore be of significant public health concern. Despite the possibility for undiagnosed spillover, the lack of detection makes it unlikely that pathogenic henipaviruses from E. helvum are regularly crossing the species barrier and undergoing significant sustained transmission in humans at this point in time. Spillover of NiV into pig populations in Malaysia may have occurred at least once prior to the detection of a major outbreak46, and therefore, detection of henipavirus antibodies in pigs in Ghana47 warrants further study. Although no human cases of LBV infection have been reported, this virus causes clinical rabies in other mammalian hosts2, and may not be detected as a cause of human rabies unless specific molecular-based LBV assays are performed.
Changes in bat–human interactions and bat–domestic animal interactions are hypothesized to be a catalyst for the zoonotic spillover of novel viruses from wildlife. Stressors, such as habitat loss, land-use change and increasing bat–human interactions may precipitate viral spillover from bats to other species23. Understanding viral persistence and the potential for spillover in African bat populations in the face of extensive hunting, logging and human population growth is of central importance for both public health and conservation, especially as these processes can be expected to increase over time.
Methods
Sampling
All fieldwork was undertaken under permits granted by national and local authorities, with ethical approval from the Zoological Society of London Ethics Committee (project reference WLE/0489). Personal protective equipment (long clothing, facemasks, eye protection and gloves) was worn during sample collection. Sampling was conducted in geographically widespread E. helvum populations along longitudinal and latitudinal axes across the species’ range (Fig. 1, Supplementary Data 1). In São Tomé, bats were obtained in collaboration with local hunters, who hunted at roost sites during the day or feeding sites at night. Elsewhere, bats were captured at the roost with mist nets (6–18 m; 38 mm) as they departed the roost site at dusk, or returned at dawn.
Female reproductive status was assigned as non-reproductive, pregnant or lactating, assessed visually or via abdominal palpation. Age was assessed by morphological characteristics and all individuals could be allocated into one of four age classes: neonate (<2 months), juvenile (J; 2—<6 months), sexually immature (SI; 6—<24 months) or adult (A; ≥24 months). For a subset of samples, the timing of sampling allowed further classification of SI individuals into 6-month age groups SI.1, SI.2 and SI.3 (6<12, 12—<18, 18—<24 months, respectively).
Genetic and blood samples were collected under manual restraint. Wing membrane biopsies (4 mm) were placed into 70% alcohol. Up to 1 ml blood was collected from the propatagial vein using a citrated 1-ml syringe and placed into a plain 1.5 ml Eppendorf tube. Pooled urine samples (up to 500 μl) were collected by pipette from plastic sheeting placed under E. helvum colonies in Tanzania and Uganda at dawn12, or directly from individual bats (in Tanzania, Malawi, Zambia and Annobón), and frozen at −80 °C without preservative. ‘Populations’ were initially defined arbitrarily based on national borders related to roost location.
Molecular methods
Genomic DNA was extracted from E. helvum tissues using DNeasy Blood and Tissue Kits (Qiagen, Crawley, West Sussex, UK) and was supplied for one E. dupreanum bat from Madagascar by the Institut Pasteur de Madagascar. Multiplexed genotyping was performed using 18 loci in six multiplexed reactions (TSY, FWB, MNQX, AgPK, AcAfAi, AdAh) using a Type-it Multiplex PCR Master Mix (Qiagen). From 20 E. helvum loci developed in a previous study48, Loci E and Ae were discarded due to difficulty in scoring or high error rates and data pertaining to locus Ag were re-binned and re-scored, correcting earlier issues with allelic dropout. Positive and negative controls were included on each plate. Amplification of mtDNA cytb gene fragments from continental samples used generic primers L14722 (5′-CGA AGC TTG ATA TGA AAA ACC ATC GTT G-3′)49 and H15149 (5′- AAA CTG CAG CCC CTC AGA ATG ATA TTT GTC CTC A-3′)50 in 20 μl reactions, containing 0.1–1 ng template DNA, 0.2 μM of each primer, 0.25 mM of each dNTP, 1.5 mM MgCl2, 0.25 μl of Taq polymerase (Invitrogen), and 0.2 μl 10 × reaction buffer and with the following conditions: 5 min at 94 °C; 40 cycles of 1 min at 93 °C, 1 min at 54 °C and 2 min at 72 °C; then 7 min at 72 °C. Although these generic primers were adequate with continental samples (8% PCR failure), amplification from isolated Gulf of Guinea island samples was less successful (48% PCR failure). Shortened primers (EhM2814 (5′-GCT TGA TAT GAA AAA CCA TCG TTG-3′) and EhM2815 (5′-CAG CCC CTC AGA ATG ATA TTT GT-3′) resulted in successful amplification when using Microzone MegaMix-Gold reagent (Microzone UK). PCRs were performed in 20 μl reactions, containing 2 ng template DNA, 0.25 μM of each primer and 10 μl MegaMix-Gold, using the following conditions: 5 min at 95 °C; 33 cycles of 30 s at 95 °C, 30 s at 53 °C and 45 s at 72 °C. PCR products were sequenced in both directions, aligned, manually checked and trimmed to 397 bp. No sequence differences were detected in 38 samples sequenced using both primer pairs, so data were combined.
RNA was extracted from urine samples using the MagMAX viral RNA isolation kit (Life Technologies, Paisley, UK), and the presence of paramyxovirus polymerase gene RNA was tested for using two heminested RT–PCRs (PAR-F2: 5′-GTT GCT TCA ATG GTT CAR GGN GAY AA-3′, PAR-R: 5′-GCT GAA GTT ACI GGI TCI CCD ATR TTN C-3′)12,51.
Genetic data analyses
After removing non-independent samples (known or suspected offspring of other individuals within the data set), cytb analyses and microsatellite analyses (at 17 loci) were performed on data from 544 and 502 individuals, respectively (Table 1). Abbreviations for population groupings used in analyses are CT (all continental populations), CB (all continental populations plus Bioko), IS (all four island populations), iIS (three isolated island populations (São Tomé, Príncipe and Annobón)) and STP (São Tomé and Príncipe) (Supplementary Fig. S3).
The statistical power of the microsatellite and mtDNA data sets to reject a null hypothesis of genetic homogeneity was assessed using the software POWSIM52. Values from the empirical data sets (number of populations, population sample sizes, number of loci and allele frequencies) were used to simulate 1,000 random sets of 12 subpopulations with expected FST values of 0.001−0.01. As mtDNA is haploid, the sample size was halved for mtDNA analyses. Power calculations indicated that the inability to detect population structure among CB populations was not because of insufficient power within the data set, and the estimated probability of falsely detecting significant differentiation was in line with the typically accepted 0.05 cutoff. Analyses were therefore continued as described below.
For microsatellite data, departures from Hardy–Weinberg equilibrium (HWE) and the presence of linkage disequilibrium among loci were assessed using FSTAT v2.9 (ref. 53) and GENEPOP v4.0.10 (ref. 54), respectively. Significance levels were adjusted for multiple testing using the false discovery rate method. Genetic diversity for each population and region was assessed by calculating observed heterozygosity (HO), expected heterozygosity (HE), and average allelic richness (RS) in FSTAT. Population structure was assessed by calculating pairwise FST values between populations and by AMOVA, as implemented in the software ARLEQUIN v3.5 (ref. 55). Significance levels were obtained with 10,000 permutations. Data were tested for the presence of IBD by regressing natural logarithm-transformed geographical distances between sampling sites (in km) against Slatkin’s linearized FST (FST/(1-FST)). Statistical significance was assessed using a Mantel test with 10,000 permutations in ARLEQUIN.
Bayesian clustering analyses were performed 20 times for each value of K (K=1–13, representing the number of populations) for 1.5 × 106 iterations with 500,000 burn-in steps using the admixture model with correlated allele frequencies in STRUCTURE56. Analyses were repeated for separate continental and island data sets. Symmetric similarity coefficients were used to assess consistency among replicate runs for each value of K using the Greedy algorithm of CLUMPP v1.1 (ref. 57), and only runs with symmetric similarity coefficients >0.8 were included in further analyses. Individual membership coefficients from replicate runs were visualized graphically using the software DISTRUCT v1.1 (ref. 58). To ensure that some loci not in HWE in the Bioko populations (see results) were not affecting clustering from this population, analyses were repeated separately with data from loci in or out of HWE in Bioko. No difference was seen in the results, and therefore remaining analyses were run with 16 loci. Assignment tests were performed in STRUCTURE and admixture was assessed using the USEPOPINFO option, using the clustering partition with the optimal mean log likelihood value as prior population information. Based on their assignment probability, P, individuals were considered non-migrant (P>0.8), admixed (0.2>P>0.8) or a recent migrant (P<0.2)59. STRUCTURE and CLUMPP analyses were performed using the CamGrid distributed computing resource. Comparable analyses were performed using spatially explicit methods; however, the results were consistent and are not presented here.
For mtDNA, in addition to AMOVA and IBD analyses, descriptive parameters of genetic diversity were calculated in the software DnaSP v5.10 (ref. 60). Rarefaction down to the minimum sample size was used to calculate haplotypic richness (a measure of diversity standardized across population sample sizes) using the software RAREFAC61. Pairwise φST values were calculated in ARLEQUIN and significance values were adjusted for multiple comparisons using the false discovery rate method. Median joining networks (MJNs) were constructed in the software NETWORK v4.6 (ref. 62). For comparison, statistical parsimony networks were constructed using TCS63, with a 95% parsimony connection limit; however, the results were consistent and are not presented here. A phylogeny of unique cytb haplotypes was reconstructed by Bayesian inference in MRBAYES v3.1.2 (ref. 64), using the E. dupreanum cytb sequence as an outgroup (which was found to be 91% (360/397 bp) identical to the consensus E. helvum cytb sequence). The most appropriate substitution model (GTR+I) was selected using PAUP* v4.0b10 (ref. 65) and MODELTEST v3.7 (ref. 66). MRBAYES was run with four simultaneous chains, sampled every 100 generations and the first 25% of trees were discarded as burn-in. Generations were added until the s.d. of split frequencies was below 0.015 (10 × 106 generations).
The relative contributions of isolation and gene flow (migration) on observed levels of population divergence were estimated using an isolation-with-migration model in IMa2 (ref. 67). Once priors had been optimized, analyses were run until stationarity was reached, which took ~2–3 months and 1.7–46 million steps, depending on sample size, before genealogy sampling commenced. Genealogy information was saved every 100 steps, and sampling was continued until ~100,000 genealogies were available for each pairwise comparison (~1 month, depending on sample size). Eight competing colonization scenarios were explored by analysing microsatellite and mtDNA data using approximate Bayesian computation methods in the software DIYABC v 1.0 (ref. 68). Eight different colonization scenarios were considered.
To construct a phylogenetic analysis of known henipaviruses and henipavirus-like viruses globally and other known Paramyxovirinae, sequences of a 559 bp segment of the polymerase gene were obtained from GenBank (Supplementary Table S3). Phylogenetic trees from these sequences and of viral sequences from urine samples analysed in this study were constructed using MRBAYES under the GTR+I+G model.
Serological analyses
The number of samples analysed using various serological assays for HeV, NiV and LBV is shown in Table 1. Antibodies against LBV (LBV.NIG56-RV1) were detected using a modified fluorescent antibody virus neutralization assay37, using the LBVNig56 isolate. Samples were tested in duplicate using threefold serial dilutions and titres corresponding to 100% neutralization of virus input are reported as IC100 endpoint reciprocal dilutions and were considered positive at >1:9.
Antibodies against henipaviruses (HeV and NiV) were detected using Luminex multiplexed microsphere binding assays and VNTs using purified recombinant expressed henipavirus soluble G glycoproteins69, which were conjugated to internally coloured and distinguishable microspheres, allowing multiplexing. Antibody binding to each microsphere was detected after conjugation of bound antibodies with biotinylated protein A and fluorescent streptavidin-R-phycoerythrin. Binding results are given as MFI values of at least 100 microspheres for each virus type, and an MFI >500 was considered positive70. Alternative, lower cutoffs were also considered based on results from mixture model analyses70. These resulted in higher seroprevalences, but no overall change in patterns to the higher, more conservative, cutoff presented here. In VNTs, samples exhibiting virus neutralization at dilutions of ≥1:10 were considered positive. Stronger results were consistently observed in NiV binding assays and VNTs24, so only NiV results are reported here. χ2 tests were used to detect significant (P<0.05) variations in seroprevalences.
Additional information
Accession numbers: All sequences reported in this study have been deposited in the GenBank nucleotide database under accession numbers KC164869 to KC164982 (Eidolon helvum cytochrome b mtDNA) and JX870901 to JX870903 (paramyxovirus polymerase gene sequences from E. helvum urine).
How to cite this article: Peel, A. J. et al. Continent-wide panmixia of an African fruit bat facilitates transmission of potentially zoonotic viruses. Nat. Commun. 4:2770 doi: 10.1038/ncomms3770 (2013).
References
Hayman, D. T. S. et al. Ecology of zoonotic infectious diseases in bats: current knowledge and future directions. Zoonoses Public Health 60, 2–21 (2012).
Banyard, A. C., Hayman, D. T. S., Johnson, N., McElhinney, L. M. & Fooks, A. R. Bats and lyssaviruses. Adv. Virus Res. 79, 239–289 (2011).
Hayman, D. T. S. et al. Evidence of Henipavirus Infection in West African Fruit Bats. PLoS One 3, e2739 (2008).
Mickleburgh, S., Hutson, A., Bergmans, W., Fahr, J. & Racey, P. A. Eidolon helvum. In: IUCN 2008. 2008 IUCN Red List of Threatened Species. Available from: www.iucnredlist.org. Downloaded on 6 February 2011 (2008).
Hayman, D. T. S. et al. Endemic Lagos bat virus infection in Eidolon helvum. Epidemiol. Infect. 140, 2163–2171 (2012).
Kuzmin, I. V. et al. Lagos bat virus in Kenya. J. Clin. Microbiol. 46, 1451–1461 (2008).
Dzikwi, A. et al. Evidence of Lagos bat virus circulation among Nigerian fruit bats. J. Wild. Dis. 46, 267 (2010).
Halpin, K. et al. Pteropid bats are confirmed as the reservoir hosts of henipaviruses: a comprehensive experimental study of virus transmission. Am J. Trop. Med. Hyg. 85, 946–951 (2011).
Marsh, G. A. et al. Cedar virus: a novel henipavirus isolated from Australian bats. PLoS Pathog. 8, e1002836 (2012).
Iehlé, C. et al. Henipavirus and Tioman virus antibodies in Pteropodid bats, Madagascar. Emerg. Infect. Dis. 13, 159–161 (2007).
Drexler, J. F. et al. Henipavirus RNA in African bats. PLoS One 4, e6367 (2009).
Baker, K. S. et al. Co-circulation of diverse paramyxoviruses in an urban African fruit bat population. J. Gen. Virol. 93, 850–856 (2012).
Drexler, J. F., Corman, V. & Müller, M. Bats host major mammalian paramyxoviruses. Nat. Commun. 3, 1–12 (2012).
Kamins, A. O. et al. Uncovering the fruit bat bushmeat commodity chain and the true extent of fruit bat hunting in Ghana, West Africa. Biol. Conserv. 144, 3000–3008 (2011).
Mallewa, M. et al. Rabies encephalitis in malaria-endemic area, Malawi, Africa. Emerg. Infect. Dis. 13, 136–139 (2007).
Turmelle, A. S., Jackson, F. R., Green, D., McCracken, G. F. & Rupprecht, C. E. Host immunity to repeated rabies virus infection in big brown bats. J. Gen. Virol. 91, 2360–2366 (2010).
Juste, J., Ibanez, C. & Machordom, A. Morphological and allozyme variation of Eidolon helvum (Mammalia: Megachiroptera) in the islands of the Gulf of Guinea. Biol. J. Linn. Soc. 71, 359–378 (2000).
Bergmans, W. Taxonomy and biogeography of African fruit bats (Mammalia, Megachiroptera). 3. The genera Scotonycteris Matshcie, 1894, Casinycteris Thomas, 1910, Pteropus Brisson, 1762, and the fruit bat Eidolon Rafinesque, 1815. Beaufortia 40, 111–177 (1990).
Sørensen, U. & Halberg, K. Mammoth roost of nonbreeding straw-coloured fruit bat Eidolon helvum (Kerr, 1792) in Zambia. Afr. J. Ecol. 39, 213–215 (2001).
Richter, H. & Cumming, G. First application of satellite telemetry to track African straw-coloured fruit bat migration. J. Zool. 275, 172–176 (2008).
Thomas, D. The annual migrations of three species of West African fruit bats (Chiroptera: Pteropodidae). Can J. Zool. 61, 2266–2272 (1983).
Hess, G. Disease in metapopulation models: implications for conservation. Ecology 77, 1617–1632 (1996).
Plowright, R. K. et al. Urban habituation, ecological connectivity and epidemic dampening: the emergence of Hendra virus from flying foxes (Pteropus spp.). Proc. R. Soc. B. 278, 3703–3712 (2011).
Peel, A. J. et al. Henipavirus neutralising antibodies in an isolated island population of African fruit bats. PLoS One 7, e30346 (2012).
George, D. B. et al. Host and viral ecology determine bat rabies seasonality and maintenance. Proc. Natl Acad. Sci. USA 108, 10208–10213 (2011).
Gonzalez, E. G., Beerli, P. & Zardoya, R. Genetic structuring and migration patterns of Atlantic bigeye tuna, Thunnus obesus (Lowe, 1839). BMC Evol. Biol. 8, 252 (2008).
Appleyard, S. A., Ward, R. D. & Grewe, P. M. Genetic stock structure of bigeye tuna in the Indian Ocean using mitochondrial DNA and microsatellites. J Fish Biol. 60, 767–770 (2002).
Küpper, C. et al. High gene flow on a continental scale in the polyandrous Kentish plover Charadrius alexandrinus. Mol. Ecol. 21, 5864–5879 (2012).
Ramachandran, S. et al. Support from the relationship of genetic and geographic distance in human populations for a serial founder effect originating in Africa. Proc. Natl Acad. Sci. USA 102, 15942–15947 (2005).
Juste, J. Allozyme variation of the Egyptian Rousette (Rousettus egyptiacus; chiroptera pteropodidae) in the Gulf of Guinea (West-Central Africa). Biochem System Ecol. 24, 499–508 (1996).
Melo, M., Warren, B. H. & Jones, P. J. Rapid parallel evolution of aberrant traits in the diversification of the Gulf of Guinea white‐eyes (Aves, Zosteropidae). Mol. Ecol. 20, 4953–4967 (2011).
Jesus, J. et al. Phylogenetic relationships of African green snakes (genera Philothamnus and Hapsidophrys) from Sao Tome, Principe and Annobon islands based on mtDNA sequences, and comments on their colonization and taxonomy. Herpetol J 19, 41–48 (2009).
Jiménez, S. & Hazevoet, C. J. First record of straw-coloured fruit bat Eidolon helvum (Kerr, 1792) for the Cape Verde Islands. Zool Caboverd 1, 116–118 (2010).
Olival, K. J. et al. Lack of population genetic structure and host specificity in the bat fly, Cyclopodia horsfieldi, across species of Pteropus bats in Southeast Asia. Parasit. Vectors 6, 1–18 (2013).
Billeter, S. A. et al. Bartonella species in bat flies (Diptera: Nycteribiidae) from western Africa. Parasitology 139, 324–329 (2012).
Biek, R. & Real, L. A. The landscape genetics of infectious disease emergence and spread. Mol. Ecol. 19, 3515–3531 (2010).
Hayman, D. T. S. et al. Antibodies against Lagos bat virus in megachiroptera from West Africa. Emerg. Infect. Dis. 14, 926–928 (2008).
Weiss, S. et al. Henipavirus-related sequences in fruit bat bushmeat, Republic of Congo [letter]. Emerg. Infect. Dis. 18,, (2012).
Breed, A. C., Field, H. E., Smith, C. S., Edmonston, J. & Meers, J. Bats without borders: long-distance movements and implications for disease risk management. Ecohealth 7, 204–212 (2010).
Epstein, J. H. et al. Pteropus vampyrus, a hunted migratory species with a multinational home-range and a need for regional management. J Appl Ecol 46, 991–1002 (2009).
Kuzmin, I. V. et al. Bats, emerging infectious diseases, and the rabies paradigm revisited. Emerg. Health Threats J doi:10.3402/ehtj.v4i0 7159 (2011).
Rahman, S. A. et al. Risk factors for Nipah virus infection among pteropid bats, Peninsular Malaysia. Emerg. Infect. Dis. 19, 51–60 (2013).
Wright, A., Rampersad, J., Ryan, J. & Ammons, D. Molecular characterization of rabies virus isolates from Trinidad. Vet Microbiol 87, 95–102 (2002).
Arguin, P. M. et al. Serologic evidence of Lyssavirus infections among bats, the Philippines. Emerg. Infect. Dis. 8, 258–262 (2002).
Lloyd-Smith, J. O. et al. Should we expect population thresholds for wildlife disease? Trends Ecol. Evol. 20, 511–519 (2005).
Pulliam, J. R. C., Field, H. E. & Olival, K. J. & Henipavirus Ecology Research Group. Nipah virus strain variation. Emerg. Infect. Dis. 11, 1978–1979 author reply 1979 (2005).
Hayman, D. T. S. et al. Antibodies to Henipavirus or Henipa-Like viruses in domestic pigs in Ghana, West Africa. PLoS One 6, e25256 (2011).
Peel, A. J., Rossiter, S. J., Wood, J. L. N., Cunningham, A. A. & Sargan, D. R. Characterization of microsatellite loci in the straw-colored fruit bat, Eidolon helvum (Pteropodidae). Conserv. Genet. Resour. 2, 279–282 (2010).
Juste, J. et al. Phylogeography of African Fruitbats (Megachiroptera). Mol Phylogenet Evol 13, 596–604 (1999).
Kocher, T. et al. Dynamics of mitochondrial DNA evolution in animals: amplification and sequencing with conserved primers. Proc. Natl Acad. Sci. USA 86, 6196–6200 (1989).
Tong, S., Chern, S.-W. W., Li, Y., Pallansch, M. A. & Anderson, L. J. Sensitive and broadly reactive reverse transcription-PCR assays to detect novel paramyxoviruses. J. Clin. Microbiol. 46, 2652–2658 (2008).
Ryman, N. & Palm, S. POWSIM: a computer program for assessing statistical power when testing for genetic differentiation. Mol. Ecol. Notes 6, 600–602 (2006).
Goudet, J. FSTAT. (Version 1.2): A computer program to calculate F-statistics. J. Hered. 86, 485–486 (1995).
Rousset, F. genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol. Ecol. Resour. 8, 103–106 (2008).
Excoffier, L. & Lischer, H. E. L. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567 (2010).
Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 155, 945–959 (2000).
Jakobsson, M. & Rosenberg, N. A. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics 23, 1801–1806 (2007).
Rosenberg, N. DISTRUCT: a program for the graphical display of population structure. Mol. Ecol. Notes 4, 137–138 (2004).
Vonholdt, B. M. et al. A novel assessment of population structure and gene flow in grey wolf populations of the Northern Rocky Mountains of the United States. Mol. Ec.ol 19, 4412–4427 (2010).
Librado, P. & Rozas, J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 1451–1452 (2009).
Petit, R. J., Mousadik, E. l.,A. & Pons, O. Identifying Populations for Conservation on the Basis of Genetic Markers. Conserv. Biol. 12, 844–855 (2008).
Bandelt, H., Forster, P. & Rohl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16, 37–48 (1999).
Clement, M., Posada, D. & Crandall, K. TCS: a computer program to estimate gene genealogies. Mol. Ecol. 9, 1657–1659 (2000).
Ronquist, F. & Huelsenbeck, J. P. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19, 1572–1574 (2003).
Swofford, D. PAUP*. Phylogenetic Analysis Using Parsimony (and Other Methods). Version 4.06b (2001).
Posada, D. & Crandall, K. MODELTEST: testing the model of DNA substitution. Bioinformatics 14, 817–818 (1998).
Hey, J. Isolation with migration models for more than two populations. Mol. Biol. Evol. 27, 905–920 (2010).
Cornuet, J. et al. Inferring population history with DIY ABC: a user-friendly approach to approximate Bayesian computation. Bioinformatics 24, 2713–2719 (2008).
Bossart, K. N. et al. Neutralization assays for differential henipavirus serology using Bio-Plex protein array systems. J. Virol. Methods 142, 29–40 (2007).
Peel, A. J. et al. Use of cross-reactive serological assays for detecting novel pathogens in wildlife: assessing an appropriate cutoff for henipavirus assays in African bats. J. Virol. Methods 295–303 (2013).
Acknowledgements
The authors thank the governments of Ghana, Tanzania, Malawi, Zambia, Uganda, Equatorial Guinea and São Tomé and Príncipe for facilitating this research. For their invaluable support in planning and implementing the fieldwork, we also thank the Zambian Wildlife Authority; Kasanka Trust; Malawi Ministry of Tourism, Wildlife and Culture, Malawi Ministry of Agriculture, Irrigation and Food Security; Tanzania Ministry of Livestock Development and Fisheries; Tanzanian Wildlife Research Institute; Sokoine University of Agriculture, Tanzania; Makerere University, Uganda; Equatorial Guinea Ministerio de Agricultura y Bosques; Equatorial Guinea Instituto del Desarrollo Forestal y Gestión de las Áreas Protegidas; Universidad Nacional de Guinea Ecuatorial; São Tomé and Príncipe Ministério de Agricultura, Desenvolivimento Rural e Pesca; Ecosystèmes Forestiers d'Afrique Centrale; Associacão Monte Pico, São Tomé and Príncipe; Alex Torrance; Lucrecia Bilé Osa Ahara;Iñaki Rodriguez Prieto; Andrés Fernandez Loras; Heidi Ruffler; Ricardo Castro Cesar de Sa; Ricardo Faustino de Lima; and Mariana Carvalho. Additional genetic samples were provided by Javier Juste, Ivan Kuzmin, Guy-Crispin Gembu Tungaluna and Jean-Michel Heraud. We also thank Daniel Horton and David Selden for constructive discussions and technical assistance, respectively. The work was supported by funding from the University of Cambridge, Institute of Zoology, CSIRO Australian Animal Health Laboratory and the UK Department for Environment, Food and Rural Affairs. The Charles Slater Trust, Zebra Foundation for Veterinary Zoological Education and Isaac Newton Trust provided grants to A.J.P. D.T.S.H. and K.S.B. were funded by the Wellcome Trust. D.T.S.H., A.R.F. and J.L.N.W. are funded by the Research and Policy for Infectious Disease Dynamics (RAPIDD) programme of the Science and Technology Directorate, Department of Homeland Security, and Fogarty International Center, USA. D.T.S.H. has a David H. Smith fellowship. C.C.B. is partially funded by National Institutes of Health, USA, grant AI054715. T.L. was funded by Lincoln Park Zoo, Chicago, US. A.R.F. was supported by the Department for Environment, Food and Rural Affairs (Defra) project SEV3500. A.A.C. is supported by a Royal Society Wolfson Research Merit award.
Author information
Authors and Affiliations
Contributions
Conceived the research: A.J.P., D.R.S., K.S.B., D.T.S.H., .R..S, A.R.F., S.J.R., J.L.N.W., A.A.C. Collected samples: A.J.P., K.S.B., D.T.S.H., R.S., T.L. Provided reagents: C.C.B., L.W., A.R.F. Performed experiments: A.J.P., K.S.B., D.T.S.H., J.A.B., G.C. Analysed data: A.J.P., D.R.S., K.S.B., D.T.S.H., J.A.B., G.C., A.R.F., S.J.R., J.L.N.W., A.A.C. Wrote or revised the manuscript: A.J.P., K.S.B., D.T.S.H., J.A.B., G.C., C.C.B., T.L., L.W., A.R.F., S.J.R., J.L.N.W., A.A.C. J.L.N.W. and A.A.C. contributed equally to this work.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Figures and Tables
Supplementary Figures S1-S3 and Supplementary Tables S1-S3. (PDF 830 kb)
Supplementary Data 1
Details of Eidolon helvum sampling locations and sample types used in this study (XLSX 17 kb)
Rights and permissions
About this article
Cite this article
Peel, A., Sargan, D., Baker, K. et al. Continent-wide panmixia of an African fruit bat facilitates transmission of potentially zoonotic viruses. Nat Commun 4, 2770 (2013). https://doi.org/10.1038/ncomms3770
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/ncomms3770
This article is cited by
Climate change increases cross-species viral transmission risk
Nature (2022)
Henipaviruses: an expanding global public health concern?
GeroScience (2022)
Fruit bats in flight: a look into the movements of the ecologically important Eidolon helvum in Tanzania
One Health Outlook (2020)
Rescue and characterization of recombinant cedar virus, a non-pathogenic Henipavirus species
Virology Journal (2018)
Support for viral persistence in bats from age-specific serology and models of maternal immunity
Scientific Reports (2018)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.