Filovirus-reactive antibodies in humans and bats in Northeast India imply zoonotic spillover

Abstract

Bats are reservoirs for several zoonotic pathogens, including filoviruses. Recent work highlights the diversity of bat borne filoviruses in Asia. High risk activities at the bat-human interface pose the threat of zoonotic virus transmission. We present evidence for prior exposure of bat harvesters and two resident fruit bat species to filovirus surface glycoproteins by screening sera in a multiplexed serological assay. Antibodies reactive to two antigenically distinct filoviruses were detected in human sera and to three individual filoviruses in bats in remote Northeast India. Sera obtained from Eonycteris spelaea bats showed similar patterns of cross-reactivity as human samples, suggesting them as the species responsible for the spillover. In contrast, sera from Rousettus leschenaultii bats reacted to two different virus glycoproteins. Our results indicate circulation of several filoviruses in bats and the possibility for filovirus transmission from bats to humans.

About the Authors:

Pilot Dovih

Roles Conceptualization, Data curation, Investigation, Writing – original draft, Writing – review & editing

Affiliations National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore, India, Sastra University, School of Chemistry and Biotechnology, Thanjavur, Tamil Nadu, India

Eric D. Laing

Roles Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

Affiliation: Uniformed Services University of the Health Sciences, Department of Microbiology and Immunology, Bethesda, Maryland, United States of America

Yihui Chen

Roles Data curation, Formal analysis, Methodology, Visualization, Writing – review & editing

Affiliation: Duke-National University of Singapore Medical School, Programme in Emerging Infectious Diseases, Singapore

Dolyce H. W. Low

Roles Investigation, Project administration, Writing – review & editing

Affiliations Duke-National University of Singapore Medical School, Programme in Emerging Infectious Diseases, Singapore, National University of Singapore, Graduate School for Integrative Sciences and Engineering, Singapore

B. R. Ansil

Roles Investigation, Visualization, Writing – review & editing

Affiliations National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore, India, Manipal Academy of Higher Education, Manipal, Karnataka, India

Xinglou Yang

Roles Resources, Writing – review & editing

Affiliation: Wuhan Institute of Virology, Department of Emerging Infectious Diseases, Wuhan, China

Zhengli Shi

Roles Resources, Writing – review & editing

Affiliation: Wuhan Institute of Virology, Department of Emerging Infectious Diseases, Wuhan, China

ORCID logo http://orcid.org/0000-0001-8089-163X

Christopher C. Broder

Roles Funding acquisition, Investigation, Resources, Writing – review & editing

Affiliation: Uniformed Services University of the Health Sciences, Department of Microbiology and Immunology, Bethesda, Maryland, United States of America

Gavin J. D. Smith

Roles Supervision, Writing – review & editing

Affiliation: Duke-National University of Singapore Medical School, Programme in Emerging Infectious Diseases, Singapore

Martin Linster

Roles Formal analysis, Investigation, Writing – original draft, Writing – review & editing

Affiliation: Duke-National University of Singapore Medical School, Programme in Emerging Infectious Diseases, Singapore

Uma Ramakrishnan

Roles Funding acquisition, Supervision, Writing – review & editing

Affiliation: National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bangalore, India

Ian H. Mendenhall

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Duke-National University of Singapore Medical School, Programme in Emerging Infectious Diseases, Singapore

ORCID logo http://orcid.org/0000-0003-4250-6459

Introduction

Filoviruses are causative agents of viral haemorrhagic disease in humans and non-human primates although virus spillover is rare [1]. There are ten distinct filoviruses classified into four genera, Ebolavirus: Ebola virus, Bundibugyo virus, Taï forest virus, Sudan virus, Reston virus and Bombali virus; Marburgvirus: Marburg virus and Ravn virus; Cuevavirus: Lloviu virus; and Dianlovirus: Měnglà virus [2–4].

Bats are the proposed natural reservoir of filoviruses, involved in enzootic virus maintenance and zoonotic virus transmission to susceptible hosts [4]. The majority of described filoviruses are endemic in the African continent, although filovirus-specific antibodies were detected in bats from Bangladesh [5], the Philippines [6], and Singapore [7]. The genome of a novel filovirus, Měnglà virus, was detected in bats from China [8] and is the second Asiatic filovirus described after Reston virus [9]. Lloviu virus was discovered in Spain in 2011 and detected in Hungary in 2016 [10, 11]. Bats are hunted by humans across Africa and Asia, and at least 167 bat species are consumed [12]. High-risk activities, such as bat hunting and mining in bat-dwelling caves, pose a threat of cross-species filovirus transmission [13].

In the Northeast Indian state of Nagaland, local ethnic groups have conducted bat harvests for at least seven generations as a source of food and traditional medicine. These bat hunters are exposed to saliva, blood, and excreta from the bat species Rousettus leschenaultii and Eonycteris spelaea. We conducted a serological survey of both hunted bat species and human hunters to study if humans have been exposed to filoviruses potentially originating from bats.

Materials and methods

Ethics statement

All study participants provided written informed consent by signing a form in their native language Nagamese. All human samples and surveys were collected under National Centre of Biological Sciences (NCBS) IEC permit 7/001 and National University of Singapore (NUS) IRB permit N-17-034E. Negative control sera were collected under the NUS IRB permit number H-18-029. All bats were sampled under the NCBS, Tata Institute of Fundamental Research IACUC permit #UR-6/2014, which adheres to provisions of the Prevention of Cruelty to Animals Act (1960) and the Breeding of and Experiments on Animals Rules (1998) and the NUS IACUC permit B16-0159 under the National Advisory Committee for Laboratory Animal Research (NACLAR) guidelines in Singapore.

Sample and data collection

In 2017, 85 individuals participating in an annual bat harvest in Mimi village (Fig 1) were provided a paper-based survey to record their gender, age, occupation and number of times involved in the bat harvest. Blood of consenting volunteers was collected in a serum separation tube (Vacutainer, Becton Dickinson, New Jersey, USA). Bat blood from E. spelaea (n = 16) and R. leschenaultii (n = 30) was collected by cardiac puncture after being sacrificed by the harvesters.

[Figure omitted. See PDF.]

Fig 1. Geographical map of the border region between India and Myanmar.

The Indian state of Nagaland and Mimi village are indicated. The map was created using QGIS v2.18.7 software (https://qgis.org/en/site/). The India shapefile was downloaded from the India Remote Sensing and GIS website (http://www.indianremotesensing.com/) and the Bangladesh and Myanmar shapefiles were downloaded from DIVA-GIS (http://www.diva-gis.org/gdata). All layers were in the geographic coordinate system WGS 84 and all software and map layers used are open access.

https://doi.org/10.1371/journal.pntd.0007733.g001

Blood samples were centrifuged at 1,500 x g for 10 min and sera were stored at 4°C until transport to NCBS, where they were held at -80°C. All sera were gamma-irradiated at 20KGy with the Blood Irradiator 2000 (BRIT, Mumbai, India) and heat-inactivated at 56°C for 60 min prior to screening. Pooled kidney, lung and spleen samples obtained from E. spelaea (n = 34) and R. leschenaultii (n = 69) were collected and stored individually in RNALater (Sigma-Aldrich). Personnel handling potentially infectious material in the field wore N95 particulate respirators, surgical gowns, face-shields and were double-gloved. Surfaces were disinfected with 3% Virkon solution. Needles and scalpel blades were single use and disposed in sharps containers. Sealed containers were autoclaved at the Healthcare and Research Center of the Naga Hospital Authority Kohima, Nagaland, India. Small amounts of tissue (lung, spleen, kidney) was excised, combined by individual bat and homogenized in AVL Buffer (Qiagen) at NCBS. An aliquot of homogenate was used for downstream PCR analysis, while another aliquot of homogenate was used for NGS analysis.

PCR, next generation sequencing and serology

RNA was extracted with PureLink RNA Mini Kit (Invitrogen), and cDNA was synthesized using SuperScript III Reverse Transcriptase (Invitrogen). Samples were tested with a nested pan-filovirus PCR assay targeting the L gene as reported previously [14]. The positive control for the PCR was synthesized by Integrated DNA Technologies based on a region of the L gene of Bundibugyo virus (Genbank Accession: KU182911). The plasmid was amplified and extracted after transformation into competent cells using a NucleoBond Xtra Maxi Plus EF kit (Macherey-Nagel, Düren, Germany) and following the manufacturer’s instruction. Purified plasmid was quantified, 10-fold serially diluted and validated with primers prior to use in the assay. Aliquots of homogenate for NGS analysis were further pooled by species (E. spelaea n = 34 and R. leschenaultii n = 34). NGS libraries were made (S1 Appendix) and validated by bioanalyzer, then sequenced on a HiSeqX Illumina machine with 2 x 150 bp reads by Medgenome Labs Ltd. (Bangalore, India) [15,16].

Human and bat sera samples were screened in a filovirus multiplex microsphere immunoassay as previously described [17]. Recombinant ectodomains of envelope attachment glycoproteins (GPe) from Ebola virus (EBOV), Bundibugyo virus (BDBV), Taï forest virus (TAFV), Sudan virus (SUDV), Reston virus (RESTV), Marburg virus (MARV), Ravn virus (RAVV), Lloviu virus (LLOV) and Měnglà virus (MLAV) (Table 1) were expressed in a mammalian cell-culture system [18, 19]. In 2018, purified, oligomeric GPe antigens (minus MLAV) were coupled to MagPlex microspheres (Luminex, Austin, TX, USA) and bat and human samples were diluted at 1:100 in PBS and run on a Bio-Plex 200 system (Bio-Rad, Hercules, California, USA) in duplicate. After sera incubation with GPe-coupled microspheres, samples were washed, incubated with biotinylated-Protein A and biotinylated Protein G (1:1 ratio) (Thermo Fisher Scientific, Waltham, MA, USA), washed and then finally incubated with streptavidin-phycoerythrin (PE) (Bio-Rad). After the discovery of Měnglà virus, the serum samples were re-run in 2019 with all GPe antigens and individual serum samples were diluted 1:100 for human and 1:250 for bat sera in PBS. Median fluorescence intensities (MFI) were measured using a MAGPIX machine (Bio-Rad) (S1, S2 and S3 Tables). Cell culture supernatant from a GPe untransfected cell line, was prepared and included in the multiplex immunoassay as a mock antigen sample to normalize non-specific antisera reactivity. Due to the absence of negative sera from the study site, we obtained seven presumptively negative human sera samples from a sample bank at Duke-NUS Medical School, Singapore. These were tested using a MAGPIX machine following the technical details described above in eight technical replicates to determine the variation of individual samples in repeat measurements to individual GPe (S4 Table).

[Figure omitted. See PDF.]

Table 1. Virus name, host and location of isolation, and accession numbers for recombinant filovirus attachment glycoproteins (GPe) used in multiplex serological binding assays.

https://doi.org/10.1371/journal.pntd.0007733.t001

Phylogenetic and statistical analysis

Next generation sequencing data was analysed as previously described [20]. Briefly, FASTQ files were trimmed for quality at a PHRED score of 20 and were then analysed in DIAMOND using the NCBI nr reference database [21]. DIAMOND outputs were analysed in MEGAN to determine the sequence similarly [22]. In the absence of negative serum controls from Nagaland, two independent methods were employed to define positive and negative cut-offs. The MFI values of the mock antigen were subtracted from each GPe MFI and the values were transformed to be positive, with the lowest number being 1. A log-normal model was fitted to the MFI data and the parameters of the fit were estimated in order to calculate the 95th percentile of the log-normal distribution (S1 Fig). To control for ebolavirus cross-reactivity [23], a separate cutoff was implemented at the three-fold change above the arithmetic mean of the mock-adjusted scaled MFI. Positive samples were defined as exceeding both thresholds. The 2018 Bio-Plex data was analysed with the same two statistical methods, but not mock-adjusted (S2 Fig). The spread of MFI values for duplicate measurements of each individual sample was plotted (S3 and S4 Figs). Mean negative human sera MFI values were plotted and the standard deviation for each sample was presented (S5 Fig). All analyses and visualizations were implemented in R 3.5.1 [24]; the R code can be retrieved from the authors upon request.

Results

The majority of bat hunters were between 18 and 50 years of age, male, and participated at least eleven times in the harvest (Table 2). All bat tissues tested were PCR-negative for filovirus-specific nucleic acid. There were a total of 13,993,300 reads from the R. leschenaultii NGS dataset and 7,975,905 reads from the E. spelaea NGS dataset and no filovirus sequences were identified. In our 2019 serum screen (that included MLAV), we detected filovirus-reactive sera in 5.9% (5/85) of human samples, 6.2% (1/16) of E. spelaea samples, and 13.3% (4/30) of R. leschenaultii samples. The highest MFI values, corresponding to levels of filovirus-specific serum IgG, were detected for EBOV-GPe in human and E. spelaea sera, and for MLAV- and RAVV-GPe in R. leschenaultii sera (Fig 2). Our results suggest human exposure to two antigenically distinct filoviruses, the first group of sera (H10, H27, H30, H40) being reactive to EBOV-, BDBV-, and SUDV-GPe, and one individual serum (H45) singly reactive to MARV-GPe. An individual E. spelaea serum sample (E34) reacted to the EBOV-, SUDV- and TAFV-GPe, displaying a similar cross-reactivity pattern as seen for the first group of human sera and was previously reported for filovirus-positive E. spelaea samples from Singapore [7]. The 2019 screening results corroborated positive samples that were detected when screened in 2018 using a Bio-Plex 200 system (S4 Fig; H27, H30, H40, E34, R39), strengthening interpretation of positive sera samples screened in two different years with two different Luminex xMAP-based machines. There was minimal intra- and inter-individual variation in the negative human sera samples and these were all well below the MFI cutoff values for each GPe (S5 Fig).

[Figure omitted. See PDF.]

Fig 2.

MFI values for sera obtained from humans (A), Eonycteris spelaea (B), Rousettus leschenaultii (C). Antibodies reactive to filovirus GPe from Ebola virus (EBOV), Bundibugyo virus (BDBV), Taï Forest virus (TAFV), Sudan ebolavirus (SUDV), Reston virus (RESTV), Lloviu virus (LLOV), Měnglà virus (MLAV), Marburg virus (MARV), and Ravn virus (RAVV) are quantified in a bead-based fluorescence assay. Grey dots represent individual samples. A boxplot is overlaid to indicate median, quartiles and extremes of the sample distribution. A black dashed line indicates the cutoff determined from a single lognormal curve-fit and a black solid black line the three-fold increase over the mean. A cladogram in panel A indicates the phylogenetic relationships of individual filovirus GPe based on their amino acid sequence.

https://doi.org/10.1371/journal.pntd.0007733.g002

[Figure omitted. See PDF.]

Table 2. Basic demographic information on human study population.

https://doi.org/10.1371/journal.pntd.0007733.t002

Discussion

Despite the growing evidence that filoviruses are present in South and Southeast Asia, there has been a historical absence of outbreaks of filovirus haemorrhagic fever in this region. Reasons why clusters of human filovirus infections have not yet been detected in Asia, include (a) ecological barriers prevent zoonotic transmission, (b) viruses are unable to sustain transmission between humans or (c) an uncharacterized diversity of non-pathogenic, antigenically-related filoviruses exist and cause asymptomatic infection in humans. Human populations with wildlife contact and no history of Ebola virus disease in Uganda [25] and the Democratic Republic of Congo [26] were reportedly ebolavirus seropositive. Similarly, here we report the presence of filovirus (e.g. ebolavirus, marburgvirus and dianlovirus) reactive antibodies in both human (e.g. bat hunters) and bat populations in Northeast India, a region with no historical record of Ebola virus disease.

Cross-reactivity between EBOV, BDBV and SUDV of the tested samples is in agreement with a previous report [23]. The mammalian cell-culture expression system to produce oligomeric, native-like GPe capture antigens in this multiplex assay provides several benefits compared to peptide-based antigen assays, including the retention of post-translational modifications (i.e. glycosylation) and native quaternary structures allowing capture of conformational-dependent antibodies. The use of GPe from all presently described filoviruses—with the exception of the recently described Bombali virus—allows for simultaneous detection and antigenic differentiation of virus species-specific IgGs and the identification of cross-reactive IgG responses. Most ebolavirus serology surveillance studies are unable to address the complex biology of known and unknown filoviruses in terms of cross-reactivity of specific antibodies. This filovirus serological assay addresses many limitations of previously employed assays by using oligomeric, native-like virus antigens in a multiplex manner and represents an improved biosurveillance tool. Establishing thresholds for low sample number serum sets that lack both positive and negative controls is challenging, but our efforts to employ two independent statistical methods and two machine platforms yielded congruent results and is in agreement with prior approaches to estimate seropositivity [27, 28].

Furthermore, this study describes serum reactivity to MLAV, incidentally in the same bat species (R. leschenaultii) mentioned in the initial report [2]. Reactivity to RESTV, which circulates endemically among bats, pigs and monkeys in the Philippines [29], and causes subclinical infection in animal care takers and slaughterhouse workers in the Philippines [9], was not detected in our study. Interestingly, MLAV and RAVV positive R. leschenaultii sera suggest circulation of two distinct filoviruses within the same species, which is serologically distinct from reactivity with E. spelaea and human samples. These results are concordant with previous findings of exposure to filoviruses antigenically closely related to EBOV, BDBV, and SUDV in E. spelaea in Singapore [17].

Two proposed mechanisms for sustained virus infection in the studied bat species are frequent co-roosting with other bats and the introduction of large numbers of susceptible juveniles into the population [30]. The two bat species sampled in this study, R. leschenaultii and E. spelaea, roost in large colonies in caves with rolling parturition patterns [31]. Though we have serological evidence of filovirus exposure, there was no genomic data detected. There are several reasons why this may be, including; small sample size, low virus copy numbers, uncertain epidemiological shedding periodicity, and high filovirus genetic diversity that is not captured by the primers employed here. Our results reinforce the need to select sentinel sites for virus surveillance at the human-animal interface and highlights some of the gaps in our understanding of filovirus transmission and ecology.

Supporting information

S1 Appendix. Technical description of next generation sequencing sample preparation, library preparation, sequencing, and bioinformatic pipeline.

https://doi.org/10.1371/journal.pntd.0007733.s001

(DOCX)

S1 Fig.

Histograms for each GPe (Ebola virus (EBOV), Bundibugyo virus (BDBV), Taï Forest virus (TAFV), Sudan virus (SUDV), Reston virus (RESTV), Lloviu virus (LLOV), Marburg virus (MARV), and Ravn virus (RAVV), Měnglà virus (MLAV) and mock antigen (MOCK)) and sera from humans (A), Eonycteris spelaea (B), and Rousettus leschenaultii (C). Lognormal distribution representing the best fit of all samples are indicated by solid black lines. A solid black vertical lines indicates 3-fold over mean and a dotted black line denotes cutoff established by lognormal curve fitting.

https://doi.org/10.1371/journal.pntd.0007733.s002

(TIF)

S2 Fig.

MFI values for sera obtained from humans (A), Eonycteris spelaea (B), Rousettus leschenaultii (C) screened in 2018 on a Bio-Plex machine. Antibodies reactive to filovirus GPe from Ebola virus (EBOV), Bundibugyo virus (BDBV), Taï Forest virus (TAFV), Sudan ebolavirus (SUDV), Reston virus (RESTV), Lloviu virus (LLOV), Marburg virus (MARV), and Ravn virus (RAVV) are quantified in a bead-based fluorescence assay. Grey dots represent individual samples. A boxplot is overlaid to indicate median, quartiles and extremes of the sample distribution. A black dashed line indicates the cutoff determined from a single lognormal curve-fit and a black solid black line the three-fold increase over the mean.

https://doi.org/10.1371/journal.pntd.0007733.s003

(TIF)

S3 Fig.

For the 2018 Bio-Plex dataset, mean values (horizontal lines) and spread of the two individual measurements are shown (vertical lines) for sera from human (A), Eonycteris spelaea (B) and Rousettus leschenaultii (C).

https://doi.org/10.1371/journal.pntd.0007733.s004

(TIF)

S4 Fig.

For the 2019 MAGPIX dataset, mean values (horizontal lines) and spread of the two individual measurements are shown (vertical lines) for sera from Eonycteris spelaea (A) and Rousettus leschenaultii (B).

https://doi.org/10.1371/journal.pntd.0007733.s005

(TIF)

S5 Fig. Mean values and standard deviation of seven normal human serum samples from healthy volunteers (varying colours).

Samples were tested in our assay with the indicated antigens in eight technical replicates.

https://doi.org/10.1371/journal.pntd.0007733.s006

(TIF)

S1 Table. Raw MFI values for human samples diluted at 1:100.

https://doi.org/10.1371/journal.pntd.0007733.s007

(XLSX)

S2 Table. Raw MFI values for Eonycteris spelaea samples diluted at 1:250.

https://doi.org/10.1371/journal.pntd.0007733.s008

(XLSX)

S3 Table. Raw MFI values for Rousettus leschenaultii samples diluted at 1:250.

https://doi.org/10.1371/journal.pntd.0007733.s009

(XLSX)

S4 Table. Negative human sera control values for seven individuals with eight technical replicates.

https://doi.org/10.1371/journal.pntd.0007733.s010

(XLSX)

Acknowledgments

We would like to thank Shoba M Bomrr for his assistance in coordinating field work in Mimi as well as the porters assisting in transporting equipment and supplies to the field site. We thank Alan Hitch for his assistance with statistical analysis and Vanessa Paynter who assisted with laboratory work. We express our gratitude to Lianying Yan and Spencer Sterling for excellent technical assistance. The opinions or assertions contained herein are the private ones of the authors and are not to be construed as official or reflecting the views of the Department of Defense, the Department of the Navy, or the Uniformed Services University of the Health Sciences, and no official endorsement should be inferred. The project or effort depicted was or is sponsored by the Department of the Defense, Defense Threat Reduction Agency. The content of the information does not necessarily reflect the position or the policy of the federal government, and no official endorsement should be inferred.

Citation: Dovih P, Laing ED, Chen Y, Low DHW, Ansil BR, Yang X, et al. (2019) Filovirus-reactive antibodies in humans and bats in Northeast India imply zoonotic spillover. PLoS Negl Trop Dis 13(10): e0007733. https://doi.org/10.1371/journal.pntd.0007733

References

1. Olival KJ, Hayman DT. Filoviruses in bats: current knowledge and future directions. Viruses. 2014;6(4):1759–88. pmid:24747773

2. Yang X-L, Tan CW, Anderson DE, Jiang R-D, Li B, Zhang W, et al. Characterization of a filovirus (Měnglà virus) from Rousettus bats in China. Nature Microbiology. 2019;4(3):390–5. pmid:30617348

3. Goldstein T, Anthony SJ, Gbakima A, Bird BH, Bangura J, Tremeau-Bravard A, et al. The discovery of Bombali virus adds further support for bats as hosts of ebolaviruses. Nature microbiology. 2018:1.

4. Hayman DT. Bats as viral reservoirs. Annual review of virology. 2016;3:77–99. pmid:27578437

5. Olival KJ, Islam A, Yu M, Anthony SJ, Epstein JH, Khan SA, et al. Ebola virus antibodies in fruit bats, bangladesh. Emerg Infect Dis. 2013;19(2):270–3. pmid:23343532

6. Taniguchi S, Watanabe S, Masangkay JS, Omatsu T, Ikegami T, Alviola P, et al. Reston Ebolavirus antibodies in bats, the Philippines. Emerg Infect Dis. 2011;17(8):1559–60. pmid:21801651

7. Laing ED, Mendenhall IH, Linster M, Low DHW, Chen Y, Yan L, et al. Serologic Evidence of Fruit Bat Exposure to Filoviruses, Singapore, 2011–2016. Emerging Infectious Diseases. 2018;24(1):122–6.

8. Yang XL, Tan CW, Anderson DE, Jiang RD, Li B, Zhang W, et al. Characterization of a filovirus (Mengla virus) from Rousettus bats in China. Nature microbiology. 2019;4(3):390–5. pmid:30617348

9. Miranda MEG, Miranda NLJ. Reston ebolavirus in humans and animals in the Philippines: a review. The Journal of infectious diseases. 2011;204(suppl_3):S757–S60.

10. Kemenesi G, Kurucz K, Dallos B, Zana B, Földes F, Boldogh S, et al. Re-emergence of Lloviu virus in Miniopterus schreibersii bats, Hungary, 2016. Emerging microbes & infections. 2018;7(1):1–4.

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

12. Mildenstein T, Tanshi I, Racey PA. Exploitation of bats for bushmeat and medicine. Bats in the Anthropocene: conservation of bats in a changing world: Springer; 2016. p. 325–75.

13. Kurpiers LA, Schulte-Herbrüggen B, Ejotre I, Reeder DM. Bushmeat and emerging infectious diseases: lessons from Africa. Problematic Wildlife: Springer; 2016. p. 507–51.

14. He B, Feng Y, Zhang H, Xu L, Yang W, Zhang Y, et al. Filovirus RNA in Fruit Bats, China. Emerg Infect Dis. 2015;21(9):1675–7. pmid:26291173

15. Hansen TA, Mollerup S, Nguyen N-p, White NE, Coghlan M, Alquezar-Planas DE, et al. High diversity of picornaviruses in rats from different continents revealed by deep sequencing. Emerging microbes & infections. 2016;5(1):1–8.

16. Lewandowska DW, Zagordi O, Geissberger F-D, Kufner V, Schmutz S, Böni J, et al. Optimization and validation of sample preparation for metagenomic sequencing of viruses in clinical samples. Microbiome. 2017;5(1):94. pmid:28789678

17. Laing ED, Mendenhall IH, Linster M, Low DHW, Chen Y, Yan L, et al. Serologic Evidence of Fruit Bat Exposure to Filoviruses, Singapore, 2011–2016. Emerg Infect Dis. 2018;24(1):114–7.

18. Bossart KN, Crameri G, Dimitrov AS, Mungall BA, Feng YR, Patch JR, et al. Receptor binding, fusion inhibition, and induction of cross-reactive neutralizing antibodies by a soluble G glycoprotein of Hendra virus. J Virol. 2005;79(11):6690–702. pmid:15890907

19. Chan YP, Yan L, Feng YR, Broder CC. Preparation of recombinant viral glycoproteins for novel and therapeutic antibody discovery. Methods Mol Biol. 2009;525:31–58, xiii. pmid:19252850

20. Mendenhall IH, Wen DLH, Jayakumar J, Gunalan V, Wang L, Mauer-Stroh S, et al. Diversity and Evolution of Viral Pathogen Community in Cave Nectar Bats (Eonycteris spelaea). Viruses. 2019;11(3):250.

21. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60. pmid:25402007

22. Huson DH, Beier S, Flade I, Gorska A, El-Hadidi M, Mitra S, et al. MEGAN Community Edition—Interactive Exploration and Analysis of Large-Scale Microbiome Sequencing Data. PLoS Comput Biol. 2016;12(6):e1004957. pmid:27327495

23. MacNeil A, Reed Z, Rollin PE. Serologic cross-reactivity of human IgM and IgG antibodies to five species of Ebola virus. PLoS neglected tropical diseases. 2011;5(6):e1175. pmid:21666792

24. Team RC. R: A language and environment for statistical computing. 2013.

25. Smiley Evans T, Tutaryebwa L, Gilardi KV, Barry PA, Marzi A, Eberhart M, et al. Suspected Exposure to Filoviruses Among People Contacting Wildlife in Southwestern Uganda. J Infect Dis. 2018.

26. Mulangu S, Alfonso VH, Hoff NA, Doshi RH, Mulembakani P, Kisalu NK, et al. Serologic Evidence of Ebolavirus Infection in a Population With No History of Outbreaks in the Democratic Republic of the Congo. The Journal of Infectious Diseases. 2018;217(4):529–37. pmid:29329455

27. Peel AJ, McKinley TJ, Baker KS, Barr JA, Crameri G, Hayman DT, et al. Use of cross-reactive serological assays for detecting novel pathogens in wildlife: assessing an appropriate cutoff for henipavirus assays in African bats. Journal of virological methods. 2013;193(2):295–303. pmid:23835034

28. De Nys HM, Kingebeni PM, Keita AK, Butel C, Thaurignac G, Villabona-Arenas C-J, et al. Survey of Ebola Viruses in Frugivorous and Insectivorous Bats in Guinea, Cameroon, and the Democratic Republic of the Congo, 2015–2017. Emerging infectious diseases. 2018;24(12):2228. pmid:30307845

29. Cantoni D, Hamlet A, Michaelis M, Wass MN, Rossman JS. Risks Posed by Reston, the Forgotten Ebolavirus. mSphere. 2016;1(6).

30. Hayman DT. Biannual birth pulses allow filoviruses to persist in bat populations. Proceedings Biological sciences / The Royal Society. 2015;282(1803):20142591.

31. Mendenhall IH, Borthwick S, Neves ES, Low D, Linster M, Liang B, et al. Identification of a Lineage D Betacoronavirus in Cave Nectar Bats (Eonycteris spelaea) in Singapore and an Overview of Lineage D Reservoir Ecology in SE Asian Bats. Transbound Emerg Dis. 2016.