Skip to main content

High prevalence of IgG antibodies to Ebola virus in the Efé pygmy population in the Watsa region, Democratic Republic of the Congo

Abstract

Background

Factors related to the natural transmission of Ebola virus (EBOV) to humans are still not well defined. Results of previous sero-prevalence studies suggest that circulation of EBOV in human population is common in sub-Saharan Africa. The Efé pygmies living in Democratic Republic of the Congo are known to be exposed to potential risk factors of EBOV infection such as bush meat hunting, entry into caves, and contact with bats. We studied the pygmy population of Watsa region to determine seroprevalence to EBOV infection and possible risks factors.

Method

Volunteer participants (N = 300) aged 10 years or above were interviewed about behavior that may constitute risk factors for transmission of EBOV, including exposures to rats, bats, monkeys and entry into caves. Samples of venous blood were collected and tested for IgG antibody against EBOV by enzyme-linked immunosorbent assay (ELISA). The χ2-test and Fisher’s exact test were used for the comparison of proportions and the Student’s t-test to compare means. The association between age group and anti-EBOV IgG prevalence was analysed by a nonparametric test for trend.

Results

The prevalence of anti-EBOV IgG was 18.7 % overall and increased significantly with age (p = 0.023). No association was observed with exposure to risk factors (contacts with rats, bats, monkeys, or entry into caves).

Conclusions

The seroprevalence of IgG antibody to EBOV in pygmies in Watsa region is among the highest ever reported, but it remains unclear which exposures might lead to this high infection rate calling for further ecological and behavioural studies.

Peer Review reports

Background

The Ebolavirus genus, which is in the Filoviridae family, comprises five distinct virus species: Zaire ebolavirus (EBOV), Sudan ebolavirus (SUDV), Tai Forest ebolavirus (TAFV), Reston ebolavirus (RESTV), and Bundibugyo ebolavirus (BDBV) [1]. EBOV, SUDV, and BDBV have caused outbreaks of Ebola virus disease (EVD) in humans with case fatality proportions ranging from 25 to 90 % depending on species [2].

There is currently no specific treatment or vaccine against EBOV approved for human use. Human-to-human transmission of the virus is controlled by the implementation of strict public health procedures, including isolation of probable or confirmed cases. However, interventions that aim to prevent the primary introduction of EBOV into human populations are difficult to implement since all the determinants of virus spill over from its natural reservoirs into the human population are still unknown.

Until recently, EVD outbreak was considered as an emerging zoonotic viral disease that occurred mainly in rural areas of Central Africa. But in December 2013, an Ebola outbreak has emerged in West Africa which was unprecedented in magnitude and spread [3]. Historically, outbreaks of EVD have occurred in sub-Saharan Africa including in Sudan, Uganda, Gabon, Côte d’Ivoire, Democratic Republic of the Congo (DRC), and the Republic of the Congo. The majority of human outbreaks have been associated with the manipulation of carcasses of infected great apes or other wildlife in Gabon and Republic of the Congo [4]. In other settings where non-human primates (NHP) or other wildlife were not identified as the cause of outbreaks (e.g., in DRC, Uganda, and Sudan), epidemiological investigations based on information available from index cases have suggested that bats might be a source of EVD outbreak [5, 6]. This hypothesis has been supported by the discovery of EBOV viral gene sequences from fruit bat species, and others bats have been found to be seropositive for EBOV antigens [7]. Fruit bats were also suggested as possible original source of infection based on the ecological investigation of the current EVD outbreak in Western African countries [8]. However, it remains unclear how bats or other carriers actually transmit the virus to humans, NHPs, or other non-identified hosts, and the parameters that lead to the occurrence of outbreaks remain to be determined.

Several epidemiological studies have been carried out to determine the prevalence of past EBOV infections and associated risk factors in human populations in countries where EVD outbreaks have occurred (e.g., DRC and Gabon) [913] and in countries without a history of epidemics (e.g., Cameroon, Central Africa and Liberia) [1419]. Together, these studies have demonstrated that, although EVD outbreaks are sporadic, exposure to EBOV or ebola-like virus is not a rare event in human populations in sub-Saharan Africa and that seroprelavence is higher in the rainforest ecosystem and in hunter populations.

Pygmies are hunter-gatherer populations living in tropical rainforests in Africa and are known for their nomadic life style that potentially exposes them to EBOV [20]. Ebola seroprevalence studies have been conducted in pygmy populations living in Central African countries including Gabon, Cameroon, and Central African Republic [10, 14, 15]. In DRC, where at least six documented outbreaks of EVD have occurred in the last two decades, several pygmy groups live in the equatorial rainforests, but there are no published data regarding their exposures to EBOV. In 2002, a Marburg haemorrhagic fever (MHF) seroprevalence survey was conducted in the Efé pygmy population living in the Watsa region in north-eastern DRC, where an MHF outbreak occurred between 1998 and 2001 [21]. The main victims of this MHF outbreak were gold miners from the Gorumbwa mine in Durba village and their family members [22]. There was zero prevalence of IgG anti-Marburg antibody in 300 pygmies [21], even though they were reported to be more exposed to risk factors (e.g., hunting and butchering bush meat, entry into caves, and contacts with bats) associated with filovirus infection [7, 2325] compared to the general, non-pygmy population.

To explore whether the Efé pygmy lifestyle in the Watsa region is related to exposure to EBOV, we investigated previously collected serum samples for evidence of past EBOV infection and tested associations with potential risk factors.

Methods

The area, study population, interviews, and blood collections are described elsewhere [21]. In brief, the Watsa region in north-eastern DRC had approximately 180,000 inhabitants including 4,000 pygmies living mainly in the southern part of the region. The study participants came from 39 different settlements and voluntarily reported to study sites for interview. During a three days survey in August 2002, after obtaining informed verbal consent, participants were asked about their subsistence activities and contact with wild animals as potential risk factors of primary filovirus transmission. Exposures to persons presumably sick with haemorrhagic fever (define as severe illness with high fever and bleeding from the nose, mouth, and/or anus) in the hospital and at home, receiving invasive medical or traditional treatment were asked as potential risk factors of secondary filovirus transmission.

We collected a blood sample from each volunteer and stored sera were tested for IgG antibodies to EBOV antigen by enzyme-linked immunosorbent assay (ELISA) at the National Institute for Communicable Diseases, Johannesburg as previously described [26]. Briefly, ELISA plates were coated with EBOV antigens diluted 1:1000 in PBS, overnight at +4 °C. Uninfected Vero cell culture antigens were coated in the same conditions on ELISA plates as controls. Diluted sera 1:400 in 5 % non-fat milk in PBS-Tween 20 (0.1 %) were incubated in the wells overnight at +4 °C. Peroxidase-labelled antibody to human IgG was added to wells and the TMB detector system was used to detect binding. We computed the corrected optical density (OD) as the optical density of the EBOV-coated well subtracted by the OD of the corresponding control well. Corrected ELISA optical density values were expressed as percent positivity (PP) of a human serum sample confirmed positive for EBOV and used as an internal control. Cut-off value for recording positive results were deliberately selected to be stringent at 2 × (mean +3 standard deviations (SD)) PP values determined for 60 serum samples from South African controls subjects who were almost certain to be seronegative.

The interview data were recorded in a database using EpiInfo 6.04 (Centers of Disease Control and Prevention, Atlanta, Georgia, USA), and analysed with STATA 12 (College Station, Texas, USA). When comparing characteristics between seropositive and seronegative individuals the χ2-test and Fisher’s exact test were used to compare proportions and the Student t-test was used to compare means. We computed SD for means and confidence intervals (CI) of proportions using the exact method. Odds Ratios (OR) and exact 95 % confidence intervals (CI) were used to access the association between risk factors and EBOV IgG seroprevalence. The association between age group and anti-EBOV IgG prevalence was analysed by a nonparametric Cuzick’s test for trend [27].

The representative of the Ministry of Health in the Watsa region and the ethics committee of the Institute of Tropical Medicine in Antwerp approved the study protocol.

Results

As previously published, a total of 300 pygmies were enrolled into the study and provided blood samples [21]. The proportion of male and female individuals was comparable. EBOV IgG antibodies were found in 56 of the 300 participants, resulting in an overall prevalence of anti-EBOV IgG in the pygmy population of 18.7 % (95 % CI 14.4–23.5 %): 20 % (95 % CI 13.9–27.3 %) in males and 17 % (95 % CI 11.7–24.4 %) in females. The difference between men and women was not statistically significant in any age group (χ2-test p = 0.35).

The mean age of the study participants was 32 years (±SD 14.6). Antibody-positive participants were significantly older than antibody-negative participants (35.9 ± 14.7 vs. 31.3 ± 14.5; Student t-test p = 0.03). EBOV antibody-positive participants were seen in all age groups, the youngest IgG EBOV-positive participant being 12 and the oldest 65 years old. A significant (Cuzick’s test p = 0.023) linear increase in EBOV IgG prevalence was noted with increasing age group (Fig. 1).

Fig. 1
figure 1

Age distribution of individuals with IgG EBOV antibodies in the Watsa area, DRC (p = 0.023)

There was no statistically significant difference in EBOV seroprevalence between groups exposed or unexposed to potential risk factors of primary filovirus transmission (hunting, entering caves, and contact with wild animals) or of secondary transmission in households or hospitals (Table 1).

Table 1 Antibodies to EBOV and potential risk factors for Ebola virus disease fever in pygmies residing in the Watsa Region, DRC

The IgG EBOV prevalence was higher in study participants who had experienced a febrile haemorrhagic syndrome at least once in their life than those who had not, but the difference was not statistically significant (22 % vs. 16 %; χ2-test p = 0.18).

Discussion

Here, we examined the seroprevalence of EBOV in a pygmy population in the Watsa area of DRC to establish whether their previously reported increased exposure to risk factors than the general population was associated with filovirus infection.

There was a very high EBOV antibody prevalence (18.7 %) in this study population, consistent with previous studies of pygmy populations in Cameroon (14.5 %) [15], Central African Republic (7.0 %), [14] and Gabon (11.1 %) [10]. Ebola ELISA was used as the serological test in the studies in Central African Republic and Gabon, whereas an immunofluorescence antibody (IFA) test was used in the Cameroon study. The ELISA test we used had Zaire ebolavirus species as source of antigen. This assay is more specific and sensitive for the detection of past infection than the standard IFA test [28, 29] and has the potential to detect cross-species ebolavirus antibodies [30]. Therefore, the pygmy seroprevalence reflects a relatively high level of exposure to EBOV or to other cross-reactive ebolavirus species in the study area. Since the EBOV species is known to be very pathogenic with a high case fatality proportion, this seroprevalence might underestimate the level of exposure to EBOV by assessing only the survivors. Alternately, our study population may have been exposed to other known or unknown ebolavirus species of lower pathogenicity, which may cause cross-reaction with the assay and contribute to the high seroprevalence [30, 31].

This study was a retrospective survey using stored pygmy sera, and comparison sera from the general sedentary population in the Watsa area were unavailable. However, two EBOV seroprevalence studies using the same ELISA assay conducted in the sedentary population in Kikwit city (the epicenter of a major EVD outbreak in 1995) and neighbouring villages, DRC, reported a prevalence of 2.2 % and 9.3 %, respectively [9]. The fact that the Efé pygmy population in Watsa tended to have a higher EBOV prevalence than non-pygmy inhabitants of Kikwit city and surrounding villages is in line with studies in Cameroon, Gabon, and Centre African Republic [10, 14, 15].

We found that EBOV IgG seroprevalence in Watsa pygmies increases significantly with age and is maximal in the 60 years and older age group. Other EBOV seroprevalence studies have shown that seroprevalence was highest in individuals aged between 20 and 40 years [9, 10, 14, 15]. The different distributions of seroprevalence by age may be due to differences in lifestyle between study populations. Non-pygmies may be more likely to enter deep into the rain forest for economic survival activities and to be exposed to a variety of putative EBOV animal hosts when they are in the economically most active 20–40 year age group. The Efé pygmies, however, who live in the rain forest all their life, may accumulate exposure to putative EBOV animal hosts throughout life, which is reflected by a linear seroprevalence increase with age. The fact that EBOV seroprevalence in pygmies in RCA [14] and Gabon [10] was highest in the 20-40 year old age group may reflect the trend that Aka and Baka pygmies in those areas have started to abandon their traditional semi-nomadic lifestyle in favour of the sedentary farmer’s life [32].

Here, there was no association between risk factors for primary or secondary transmission of filoviruses and EBOV antibody status. Some of the non-significant differences observed may be real, but our study may have been insufficiently powered to detect them. For instance, categorisation of the level of exposure for participants in contact with someone suffering from haemorrhagic fever might reveal significant differences. Furthermore, our study was limited by the fact that study participants were not selected randomly but were self-recruited, although we have no reason to believe that this resulted in selection bias.

Conclusion

Pygmy populations are known to be particularly at risk of filovirus exposure compared to sedentary populations. Here, pygmies in the Watsa region of DRC were at high risk and had one of the highest reported seroprevalence of EBOV IgG antibodies using a well-established ELISA assay. Since we could not identify risk factors that could explain this high prevalence, further ecological and behavioural studies might be necessary to better understand determinants of EBOV exposure in this population.

Abbreviations

BDBV: bundibugyo ebolavirus, CI: confidence intervals, DRC: Democratic Republic of the Congo, EBOV: zaire ebolavirus, ELISA: enzyme-linked immunosorbent assay, EVD: ebola virus disease, IFA: immunofluorescence antibody, MHF: marburg Hemorrhagic Fever, NHP: non-human primates, OD: optical density, OR: odds ratios, PP: percent positivity, RESTV: reston ebolavirus, SD: standard deviations, SUDV: sudan ebolavirus, TAFV: tai forest ebolavirus

References

  1. Kuhn JH, Bao Y, Bavari S, Becker S, Bradfute S, et al. Virus nomenclature below the species level: a standardized nomenclature for natural variants of viruses assigned to the family Filoviridae. Arch Virol. 2013;158:301–11.

    Article  CAS  PubMed  Google Scholar 

  2. Feldmann H, Geisbert TW. Ebola haemorrhagic fever. Lancet. 2011;377:849–62.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Baize S, Pannetier D, Oestereich L, Rieger T, Koivogui L, et al. Emergence of Zaire Ebola virus disease in Guinea. N Engl J Med. 2014;371:1418–25.

    Article  CAS  PubMed  Google Scholar 

  4. Leroy EM, Rouquet P, Formenty P, Souquiere S, Kilbourne A, et al. Multiple Ebola virus transmission events and rapid decline of central African wildlife. Science. 2004;303:387–90.

    Article  CAS  PubMed  Google Scholar 

  5. Baron RC, McCormick JB, Zubeir OA. Ebola virus disease in southern Sudan: hospital dissemination and intrafamilial spread. Bull World Health Organ. 1983;61:997–1003.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Leroy EM, Epelboin A, Mondonge V, Pourrut X, Gonzalez JP, et al. Human Ebola outbreak resulting from direct exposure to fruit bats in Luebo, Democratic Republic of Congo, 2007. Vector Borne Zoonotic Dis. 2009;9:723–8.

    Article  PubMed  Google Scholar 

  7. Leroy EM, Kumulungui B, Pourrut X, Rouquet P, Hassanin A, et al. Fruit bats as reservoirs of Ebola virus. Nature. 2005;438:575–6.

    Article  CAS  PubMed  Google Scholar 

  8. Mari Saez A, Weiss S, Nowak K, Lapeyre V, Zimmermann F, et al. Investigating the zoonotic origin of the West African Ebola epidemic. EMBO Mol Med. 2014;7:17–23.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Busico KM, Marshall KL, Ksiazek TG, Roels TH, Fleerackers Y, et al. Prevalence of IgG antibodies to Ebola virus in individuals during an Ebola outbreak, Democratic Republic of the Congo, 1995. J Infect Dis. 1999;179 Suppl 1:S102–107.

    Article  PubMed  Google Scholar 

  10. Heffernan RT, Pambo B, Hatchett RJ, Leman PA, Swanepoel R, et al. Low seroprevalence of IgG antibodies to Ebola virus in an epidemic zone: Ogooue-Ivindo region, Northeastern Gabon, 1997. J Infect Dis. 2005;191:964–8.

    Article  CAS  PubMed  Google Scholar 

  11. Nkoghe D, Padilla C, Becquart P, Wauquier N, Moussavou G, et al. Risk Factors for Zaire ebolavirus--specific IgG in rural Gabonese populations. J Infect Dis. 2011;204 Suppl 3:S768–775.

    Article  CAS  PubMed  Google Scholar 

  12. Ivanoff B, Duquesnoy P, Languillat G, Saluzzo JF, Georges A, et al. Haemorrhagic fever in Gabon. I. Incidence of Lassa, Ebola and Marburg viruses in Haut-Ogooue. Trans R Soc Trop Med Hyg. 1982;76:719–20.

    Article  CAS  PubMed  Google Scholar 

  13. Van der Groen G, Pattyn SR. Measurement of antibodies to Ebola virus in human sera from N. W.-Zaire. Ann Soc Belg Med Trop. 1979;59:87–92.

    PubMed  Google Scholar 

  14. Gonzalez JP, Nakoune E, Slenczka W, Vidal P, Morvan JM. Ebola and Marburg virus antibody prevalence in selected populations of the Central African Republic. Microbes Infect. 2000;2:39–44.

    Article  CAS  PubMed  Google Scholar 

  15. Bouree P, Bergmann JF. Ebola virus infection in man: a serological and epidemiological survey in the Cameroons. Am J Trop Med Hyg. 1983;32:1465–6.

    CAS  PubMed  Google Scholar 

  16. Gonzalez JP, Josse R, Johnson ED, Merlin M, Georges AJ, et al. Antibody prevalence against haemorrhagic fever viruses in randomized representative Central African populations. Res Virol. 1989;140:319–31.

    Article  CAS  PubMed  Google Scholar 

  17. Meunier DM, Johnson ED, Gonzalez JP, Georges-Courbot MC, Madelon MC, et al. Current serologic data on viral hemorrhagic fevers in the Central African Republic. Bull Soc Pathol Exot Filiales. 1987;80:51–61.

    CAS  PubMed  Google Scholar 

  18. Saluzzo JF, Gonzalez JP, Herve JP, Georges AJ, Johnson KM. Preliminary note on the presence of antibodies to Ebola virus in the human population in the eastern part of the Central African Republic. Bull Soc Pathol Exot Filiales. 1980;73:238–41.

    CAS  PubMed  Google Scholar 

  19. Knobloch J, Albiez EJ, Schmitz H. A Serological Survey on Viral Hemorrhagic Fevers in Liberia. Annales De Virologie. 1982;133:125–8.

    Google Scholar 

  20. Johnson ED, Gonzalez JP, Georges A. Filovirus activity among selected ethnic groups inhabiting the tropical forest of equatorial Africa. Trans R Soc Trop Med Hyg. 1993;87:536–8.

    Article  CAS  PubMed  Google Scholar 

  21. Borchert M, Mulangu S, Swanepoel R, Tshomba A, Afounde A, et al. Pygmy populations seronegative for Marburg virus. Emerg Infect Dis. 2005;11:174–7.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Bausch DG, Borchert M, Grein T, Roth C, Swanepoel R, et al. Risk factors for Marburg hemorrhagic fever, Democratic Republic of the Congo. Emerg Infect Dis. 2003;9:1531–7.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Bausch DG, Nichol ST, Muyembe-Tamfum JJ, Borchert M, Rollin PE, et al. Marburg hemorrhagic fever associated with multiple genetic lineages of virus. N Engl J Med. 2006;355:909–19.

    Article  CAS  PubMed  Google Scholar 

  24. Towner JS, Amman BR, Sealy TK, Carroll SA, Comer JA, et al. Isolation of genetically diverse Marburg viruses from Egyptian fruit bats. PLoS Pathog. 2009;5:e1000536.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Morvan JM, Deubel V, Gounon P, Nakoune E, Barriere P, et al. Identification of Ebola virus sequences present as RNA or DNA in organs of terrestrial small mammals of the Central African Republic. Microbes Infect. 1999;1:1193–201.

    Article  CAS  PubMed  Google Scholar 

  26. Ksiazek TG, West CP, Rollin PE, Jahrling PB, Peters CJ. ELISA for the detection of antibodies to Ebola viruses. Journal of Infectious Diseases. 1999;179:S192–8.

    Article  CAS  PubMed  Google Scholar 

  27. Cuzick J. A Wilcoxon-type test for trend. Stat Med. 1985;4:87–90.

    Article  CAS  PubMed  Google Scholar 

  28. Ksiazek TG, West CP, Rollin PE, Jahrling PB, Peters CJ. ELISA for the detection of antibodies to Ebola viruses. J Infect Dis. 1999;179 Suppl 1:S192–198.

    Article  CAS  PubMed  Google Scholar 

  29. Van der Waals FW, Pomeroy KL, Goudsmit J, Asher DM, Gajdusek DC. Hemorrhagic fever virus infections in an isolated rainforest area of central Liberia. Limitations of the indirect immunofluorescence slide test for antibody screening in Africa. Trop Geogr Med. 1986;38:209–14.

    PubMed  Google Scholar 

  30. Macneil A, Reed Z, Rollin PE. Serologic cross-reactivity of human IgM and IgG antibodies to five species of Ebola virus. PLoS Negl Trop Dis. 2011;5:e1175.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Monath TP. Ecology of Marburg and Ebola viruses: speculations and directions for future research. J Infect Dis. 1999;179 Suppl 1:S127–138.

    Article  PubMed  Google Scholar 

  32. Hewlett B, Vandekoppel JMH, Cavallisforza LL. Exploration Ranges of Aka Pygmies of the Central-African-Republic. Man. 1982;17:418–30.

    Article  Google Scholar 

Download references

Acknowledgements

We thank the members of the Watsa Rural Health Zone Central Office, particularly Mr Mwimba Arajebo, for making demographic data of the zone available, and the interviewers who made this investigation possible despite difficult working conditions.

We thank Dr Julius Lutwama at Uganda Virus Research Institute, Entebbe, Uganda for granting temporary storage of the samples and facilitating their shipment to Johannesburg.

We are indebted to the study participants for their trust and availability. We gratefully acknowledge editorial assistance as part of the Nextgenediting Global Initiative (www.nextgenediting.com).

Funding

This work was supported by ‘Fonds voor Wetenschappelijk Onderzoek—Vlaanderen’ [1.5.188.01], Antwerp; the Framework Agreement between the Belgian Directorate for Development Co-operation and the Institute of Tropical Medicine, Antwerp and the Wellcome Trust Grant [WT087546MA] to the Southern African Centre for Infectious Diseases and Surveillance (SACIDS). These funding institutions have played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

The data supporting our findings can be found in the Additional file 1.

Authors’ contributions

Conceived and designed the study: SM, MB. Analysed the data: SM, MB. Performed the laboratory experiments: JP, RS. Contributed data collection and analysis: SM, AT, AA AK, JJM and PVS. Wrote the paper: SM, MB. All authors read and approved the final manuscript.

Competing interest

The authors have declare that no competing interests exist.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Ethical clearance for the study protocol was obtained by the ethics committee of the Institute of Tropical Medicine in Antwerp and the representative of the Ministry of Health in Watsa. All participants provided verbal informed consent for participation in the study. Where participants were children, a parent or guardian provided the consent.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Sabue Mulangu.

Additional file

Additional file 1:

Efe pygmy, Watsa, DRC: EBOV exposures data. (XLSX 49 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://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. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mulangu, S., Borchert, M., Paweska, J. et al. High prevalence of IgG antibodies to Ebola virus in the Efé pygmy population in the Watsa region, Democratic Republic of the Congo. BMC Infect Dis 16, 263 (2016). https://doi.org/10.1186/s12879-016-1607-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-016-1607-y

Keywords