RESEARCH
A Bayesian network approach to the study of historical epidemiological databases: modelling meningitis outbreaks in the Niger
Une approche de réseau bayésien de l'étude des bases de données épidémiologiques historiques: la modélisation des poussées de méningite au Niger
Un enfoque de red bayesiana para el estudio de bases históricas de datos epidemiológicos: modelo de brotes de meningitis en Níger
A BeresniakI, *; E BertheratII; W PereaII; G SogaIII; R SouleyIV; D DupontV; S HugonnetII
IData Mining International, Route de l'Aeroport 29-31, CP 221, 1215 Geneva 15, Switzerland
IIEpidemic and Pandemic Alert and Response, World Health Organization, Geneva, Switzerland
IIIWorld Health Organization, Niamey, Niger
IVDepartment of Statistics, Surveillance and Response to Epidemics, Ministry of Health, Niamey, Niger
VData Mining America, Montreal, Canada
ABSTRACT
OBJECTIVE: To develop a tool for evaluating the risk that an outbreak of meningitis will occur in a particular district of the Niger after outbreaks have been reported in other, specified districts of the country.
METHODS: A Bayesian network was represented by a graph composed of 38 nodes (one for each district in the Niger) connected by arrows. In the graph, each node directly influenced each of the "child" nodes that lay at the ends of the arrows arising from that node, according to conditional probabilities. The probabilities between "influencing" and "influenced" districts were estimated by analysis of databases that held weekly records of meningitis outbreaks in the Niger between 1986 and 2005. For each week of interest, each district was given a Boolean-variable score of 1 (if meningitis incidence in the district reached an epidemic threshold in that week) or 0.
FINDINGS: The Bayesian network approach provided important and original information, allowing the identification of the districts that influence meningitis risk in other districts (and the districts that are influenced by any particular district) and the evaluation of the level of influence between each pair of districts.
CONCLUSION: Bayesian networks offer a promising approach to understanding the dynamics of epidemics, estimating the risk of outbreaks in particular areas and allowing control interventions to be targeted at high-risk areas.
RÉSUMÉ
OBJECTIF: Développer un outil permettant d'évaluer le risque d'apparition d'une poussée de méningite dans un district particulier du Niger après le signalement d'autres poussées dans d'autres districts spécifiés du pays.
MÉTHODES: Un réseau bayésien a été représenté par un graphique composé de 38 noeuds (un pour chaque district du Niger) reliés par des flèches. Dans le graphique, chaque noeud a directement influencé chacun des noeuds «enfants» aux extrémités des flèches résultant de ces noeuds, selon les probabilités conditionnelles. Les probabilités entre les districts «influençant» et «influencés» ont été estimées par l'analyse des bases de données qui contenaient les enregistrements hebdomadaires des poussées de méningite au Niger entre 1986 et 2005. Pour chaque semaine d'intérêt, on a attribué à chaque district un score booléen variable de 1 (si l'incidence de méningite dans le district atteignait un seuil épidémique au cours de cette semaine) ou de 0.
RÉSULTATS: L'approche de réseau bayésien a fourni des informations originales et importantes, permettant d'identifier les districts qui influencent le risque de méningite dans d'autres districts (et les districts qui sont influencés par un district particulier) et d'évaluer le niveau d'influence entre chaque couple de districts.
CONCLUSION: Les réseaux bayésiens offrent une approche prometteuse pour comprendre la dynamique des épidémies, évaluant le risque de poussées dans des zones particulières et permettant de cibler les interventions de lutte contre la maladie dans les zones à risque élevé.
RESUMEN
OBJETIVO: Desarrollar una herramienta para evaluar el riesgo de la aparición de un brote de meningitis en un distrito determinado del Níger después de que se haya informado acerca de otros brotes en otros distritos específicos del país.
MÉTODOS: Se representó una red bayesiana con un gráfico compuesto por 38 nodos (uno por cada distrito en el Níger) conectados mediante flechas. En el gráfico, cada nodo influía directamente en cada uno de los nodos de "niños" que se encuentran en los extremos de las flechas que surgen de dicho nodo, de acuerdo con las probabilidades condicionales. Se calcularon las probabilidades entre distritos "influyentes" e "influidos" mediante el análisis de bases de datos que recogían registros semanales de los brotes de meningitis en Níger entre 1986 y 2005. Por cada semana de interés, se adjudicaba a cada uno de los distritos un valor de variable booleana de 1 (si la incidencia de meningitis en el distrito alcanzaba un umbral epidémico en esa semana) o de 0.
RESULTADOS: El enfoque de red bayesiana proporciona información importante y original, lo que permite identificar cuáles son los distritos que influyen en el riesgo de meningitis de otros distritos (y cuáles están bajo la influencia de cualquier distrito determinado) y evaluar el nivel de influencia entre cada par de distritos.
CONCLUSIÓN: Las redes bayesianas ofrecen un enfoque prometedor para entender las dinámicas de las epidemias, permiten calcular el riesgo de brotes en áreas determinadas y concentrar las intervenciones de control objetivo en las áreas de alto riesgo.
Introduction
Throughout the African "meningitis belt", epidemics of meningococcal disease have been reported since the disease was first described, in 1912.1 Every year, western African countries within the Sahelo-Sudanian band experience major outbreaks of meningococcal meningitis, each of which can affect up to 200 000 people, most of them young children.2 Burkina Faso, Mali and the Niger, for example, are regularly hit by meningitis epidemics.3-7 Meningitis is characterized by high levels of seasonal endemicity, with large epidemics of meningococcal meningitis occurring cyclically.8 Each such epidemic typically starts in January and ends in late May.2 The three main pathogens causing bacterial meningitis in Africa are Neisseria meningitidis, Haemophilus influenzae type b and Streptococcus pneumoniae.4, 9, 10
As a major cause of morbidity and mortality in sub-Saharan Africa, meningitis merits specific control measures.11-13 Burkina Faso, Mali and the Niger have each already established a broad surveillance network for collecting data on cases of the disease.4, 6, 11, 12 Since 1986, peripheral health centres in each of these countries have routinely collected data on suspected cases of meningitis and reported each such case to health districts for subsequent analysis at the national health ministry. The information collected in this surveillance now forms important epidemiological databases that could improve our understanding of the nature of meningitis epidemics and allow high-risk areas to be identified. Current strategies for controlling meningitis epidemics in sub-Saharan Africa are based on epidemiological, immunological and logistical considerations. The strategy currently advocated by the World Health Organization (WHO) consists of the early detection of epidemics, the treatment of cases with antibiotics, and mass vaccination to halt the outbreak (if possible, within 4-6 weeks of the epidemic threshold being reached).8, 14, 15 The early detection of meningitis epidemics is based mainly on the observation of the trends in weekly incidence.8, 12, 14, 16-18
The effectiveness of mass vaccination as a strategy for the control of meningitis epidemics has been questioned.19 Given the sporadic nature of the outbreaks, the optimal use of vaccines to control both short-term epidemic and endemic meningococcal disease has been the subject of much debate.20 In particular, the results of several studies in Africa have shown that vaccination during outbreak situations is suboptimal, mainly because populations in resource-poor areas cannot be immunized rapidly enough.20 In addition, although rapid laboratory diagnosis is an essential component in the surveillance of meningococcal epidemics, as it allows decision-makers to select the most appropriate vaccine for mass vaccination,21-23 the resource-poor countries most affected by such epidemics struggle to achieve such diagnosis.21 A new conjugate vaccine that protects against serogroup A meningitis (the form involved in most outbreaks of meningococcal disease in Africa) was approved by the United States Food and Drug Administration in February 2010. Conjugate vaccines prevent both carriage and the transmission of bacteria from person to person.
Extensive descriptive analyses that have been performed on meningitis surveillance databases have provided important information on epidemic cycles, seasonality and the correlations between morbidity and certain co-factors, such as climatic parameters.2, 24-26 While the seasonal and spatial patterns of the disease appear to be linked to climate, the mechanisms responsible for these patterns have still to be elucidated.2, 25
In the present study, surveillance data and a modelling method based on Bayesian networks were used to explore how meningitis incidence in a district of the Niger was influenced by, or influenced, the incidence in any other district. The aim was to develop a method for the optimization of epidemic alerts and the spatial and temporal targeting of immunizations and other interventions for the management of meningitis in the Niger, with the ultimate goal of preventing meningitis epidemics in the country.
Methods
A Bayesian network consists of a graphical model showing the probabilistic relationships between one or more variables of interest. Bayesian methods are particularly valuable whenever there is a need to extract information from data that are uncertain or subject to any kind of error or noise.27, 28 When applied to the forecasting of epidemics, a Bayesian network can allow potential "dependence relationships", such as the probability that an outbreak will occur in one district after it has occurred in certain other districts, to be explored.
A Bayesian network is represented by a graph composed of nodes connected by arrows. Each arrow begins at a "parent" node and ends at a "child" node, with the "parent" directly influencing the connected "child" in some way. The degree of influence between each "parent" and "child" is a conditional probability that can usually be computed from the data. Each node of a Bayesian network has an associated conditional probability table. In an epidemiological analysis, the nodes might represent the districts of a country and the arrows and conditional probability table might show how a disease outbreak in one district is linked to the probability of a disease outbreak in another district. In the example shown (Fig. 1), if district A experienced an outbreak, the probabilities that districts B and C experienced outbreaks would be 85% and 15%, respectively.
Since January 1986, WHO has supervised the weekly collection of data on meningitis incidence in each district of the Niger and these data now form a substantial historical database. The data recorded include district name and population, week number and the reported numbers of suspected cases of meningitis and of deaths attributed to meningitis in that district and week. For the present study, the data from the database for the 14 calendar years from January 1986 to December 1999 were used to construct a Bayesian network, with a node for each of the 38 districts in the Niger. The weekly incidence thresholds set by WHO for epidemic meningitis (i.e. at least 10 cases per 100 000 inhabitants in a district with a population of at least 30 000 or at least five cases per 100 000 in a district with a smaller population) were used to give each district a weekly score of 1 (if an epidemic had occurred) or 0 (if no outbreak had occurred). The corresponding weekly incidence thresholds for "epidemic alert", which could be used for the early detection of a meningitis epidemic, are lower (at least five cases per 100 000 inhabitants and at least two cases per 100 000, respectively).
Two types of analysis were then performed, one (the "first-level" analysis) taking no account of the time taken for an epidemic in one district to influence the development of an epidemic elsewhere (ignoring the timing of the epidemics) and the other (the "second-level" analysis) assuming that an epidemic in one district would influence the development of an epidemic elsewhere 1, 2 or 3 years later. A maximum lag of 3 years was explored because this was both the median period in the cycles of meningitis outbreaks detected in the historical database and the estimated duration of protection resulting from a vaccination programme.
The results of the "first-level" analysis, performed with version 1.0 of the Discoverer software package (Bayesware, Milton Keynes, United Kingdom of Great Britain and Northern Ireland), were used to create a graph of the Bayesian network and included the relevant conditional probability table, showing the probabilities of meningitis outbreaks. Because of the limited capability of the Discoverer software, the analysis was restricted to the first 14 years of the meningitis database for the Niger. The "minimum description length" learning algorithm was used, with network parameters estimated via the maximum-likelihood technique. The "second-level" analysis was performed using the data collected on meningitis in the Niger between January 1986 and December 2005, a β version of the Bayesia package (Bayesia SAS, Laval, France) and the Bayes Net and BNLAT toolboxes from the Matlab package (Mathworks Inc., Natick, United States of America).
Results
Between January 1986 and December 2005, 182 244 suspected cases of meningitis and 14 859 meningitis-related deaths were recorded in the Niger. The cumulative attack rate over the 20 annual cycles of meningitis that occurred over this period was 1919.6 (suspected) cases per 100 000 population. The corresponding mean annual attack rate was 96.0 cases per 100 000 population. The Niger experienced six major national-level epidemics of meningitis between January 1986 and December 2005, with, at times, more than 50 cases recorded per 100 000 population in a single week.
First-level analysis
The results of the first-level analysis, which ignored any time lag between epidemics, are summarized in Fig. 2, with arrows running from the influencing districts to the influenced districts. From this graph, the number of districts influenced by a given district and the number of districts influencing a specific district ( Table 1), can be determined at a glance.
Conditional probabilities were calculated for all districts. As an example, Table 2 presents the conditional probabilities of an outbreak and of no outbreak in the Boboye district, according to the outbreak status of the two influencing districts: Agadez and Birni. These results indicate that, when epidemics occur in both Agadez and Birni, Boboye should be included in any vaccination strategy because the probability of an outbreak in this district is high (96.2%).
Second-level analysis
Fig. 3 (available at: http://www.who.int/bulletin/volumes/90/6/11-086009) illustrates the results of the second-level analysis. The large number of links between the nodes in the uppermost plot in this figure, which shows the links when no lag ("Year 0 to Year 0") or a lag of 1, 2 or 3 years is used, make this plot difficult to interpret. The links become clearer when they are plotted separately for lags of 1, 2 and 3 years, as in the other plots in Fig. 3. Unfortunately, this analysis did not allow the accurate calculation of conditional probabilities, mainly because observations with "no outbreak" status were too few.
Fig. 4 summarizes the results of the analysis of longitudinal surveillance data for those districts in which meningitis incidence reached the epidemic threshold in 2003 and 2004. Arrows again indicate the influencing relationships.
Discussion
The control of epidemic meningitis remains an unresolved problem in Africa, partly because the location of major outbreaks, which might be affected by relatively rare events, is so difficult to predict. The data analysed in the present study were collected over such a long period (up to 20 years) that they were probably affected by unpredicted rare events, such as unusual migratory patterns, armed conflicts and catastrophic dry seasons. In general, valid predictive statistics on such rare events cannot be derived from mathematical models, but data mining methods, when applied to clinical and epidemiological data, may allow previously unpredictable or unknown trends to be made apparent.29 In the Niger, such an approach may well help to elucidate the factors that influence or trigger meningitis outbreaks. Unfortunately, historical epidemiological databases like the one used in the present study are very rare in the context of epidemic diseases, especially in developing countries. While descriptive analyses constitute a standard approach, they do not allow optimal responses to future major public health issues to be planned.
In the present study, although the first-level analysis allowed the conditional probabilities between the influencing and influenced districts to be estimated, it made no allowance for the time needed for the inter-district influences to reveal themselves as changes in meningitis incidence. Unfortunately, when lags of 0, 1, 2, or 3 years were considered in the second-level analysis, too few "no-outbreak" observations were available for an accurate estimation of conditional probabilities. Nevertheless, both methodological approaches appear promising in the planning and timely deployment of vaccination programmes, especially given the limited duration of the protection offered by current vaccines.
A validation test of the Bayesian network model was recently carried out using a 6-year time horizon, with observed surveillance data from 2006, 2007 and 2008 compared with expected model predictions from 2004, 2005 and 2006 (data not shown). Districts reaching epidemic alert or epidemic thresholds in 2004, 2005 and 2006 were considered as potentially influencing other districts in 2006, 2007 and 2008. The proportion of the outbreaks observed in 2006-2008 that were successfully predicted using the data collected in 2004-2006 was 63% when the epidemic threshold was used for the influencing districts and 91% when the alert threshold was used. This promising performance justifies the continued interest in such modelling for guiding immunization planning, including the extension of vaccination campaigns to "influenced" districts when an "influencing" district has reached the epidemic threshold.
More sophisticated approaches should be explored. For example, it would be interesting to analyse the collected data year by year and to build a Bayesian network for each year. Re-sampling simulation techniques could be used to simulate and replicate observations, therefore allowing the construction of a Bayesian network for each simulation, together with a conditional probability table.
In conclusion, a Bayesian network approach offers an innovative and promising technique for extracting meaningful and clinically useful information from a large historical epidemiological database. As suggested by the present study, this method can improve our understanding of the dynamics of epidemic outbreaks, help make health interventions more effective and optimize resource allocations. Unfortunately, the method relies on large volumes of longitudinal records that are rarely available in the context of epidemic diseases. In this respect, the historical database on meningitis in the Niger provided a rare data set that made the testing of this promising approach possible. Unfortunately, no operational conditional probabilities regarding meningitis outbreaks could be determined when a time lag was included. However, further research allowing Bayesian networks to be combined with re-sampling techniques, would allow this issue to be addressed, to the potential benefit of those countries where relevant large epidemiological databases are unavailable. 
Acknowledgements
We thank Michel Lamure, from Claude Bernard Lyon 1 University (France), for his important contributions.
References
1. Woods CW, Armstrong G, Sackey SO, Tetteh C, Bugri S, Perkins BA et al. Emergency vaccination against epidemic meningitis in Ghana: implications for the control of meningococcal disease in West Africa. Lancet 2000;355:30-3. doi: 10.1016/S0140-6736(99)03366-8 PMID: 10615888
2. Sultan B, Labadi K, Guegan J, Janicot S. Climate drives the meningitis epidemics onset in West Africa. PLoS Med 2005;2:e6. doi: 10.1371/journal. pmed.0020006 PMID: 15696216
3. Roberts L. Infectious disease: hitting early, epidemic meningitis ravages Nigeria and Niger. Science 2009;324:20-1. doi: 10.1126/science.324.5923.20
4. Campagne G, Schuchat A, Djibo S, Ousséini A, Cissé L, Chippaux J. Epidemiology of bacterial meningitis in Niamey, Niger, 1981-96. Bull World Health Organ 1999;77:499-508. PMID: 10427935
5. de Chabalier F, Djingarey M, Hassane A, Chippaux J. Meningitis seasonal pattern in Africa and detection of epidemics: a retrospective study in Niger, 1990-98. Trans R Soc Trop Med Hyg 2000;94:664-8. doi: 10.1016/S0035- 9203(00)90224-4 PMID: 11198652
6. Boisier P, Djibo S, Sidikou F, Mindadou H, Kairo K, Djibo A et al. Epidemiological patterns of meningococcal meningitis in Niger in 2003 and 2004: under the threat of N. meningitidis serogroup W135. Trop Med Int Health 2005;10:435-43. doi: 10.1111/j.1365-3156.2005.01394.x
7. Varaine F, Caugant D, Riou J, Konde M, Soga G, Nshimirimana D et al. Meningitis outbreaks and vaccination strategy. Trans R Soc Trop Med Hyg 1997;91:3-7. doi: 10.1016/S0035-9203(97)90371-0 PMID: 9093614
8. Alonso JM, Bertherat E, Perea W, Borrow R, Chanteau S, Cohet C et al. From genomics to surveillance, prevention and control: new challenges for the African meningitis belt. Bull Soc Pathol Exot 2006;99:404-8. PMID: 17253061
9. Chanteau S, Sidikou F, Dijibo S, Moussa A, Mindadou H, Boisier P. Scaling up of PCR-based surveillance of bacterial meningitis in the African meningitis belt: indisputable benefits of multiplex PCR assay in Niger. Trans R Soc Trop Med Hyg 2006;100:677-80. doi: 10.1016/j.trstmh.2005.09.006 PMID: 16359713
10. Parent du Châtelet I, Traore Y, Gessner B, Antignac A, Naccro B, Njanpop-Lafourcade B et al. Bacterial meningitis in Burkina Faso: surveillance using field-based polymerase chain reaction testing. Clin Infect Dis 2005;40:17-25. doi: 10.1086/426436 PMID: 15614687
11. Boisier P, Mainassara H, Sidikou F, Djibo S, Kairo K, Chanteau S. Case-fatality ratio of bacterial meningitis in the African meningitis belt: we can do better. Vaccine 2007;25(Suppl 1):A24-9. doi: 10.1016/j.vaccine.2007.04.036 PMID: 17521784
12. Leake JAD, Kone ML, Yada AA, Barry LF, Traore G, Ware A et al. Early detection and response to meningococcal disease epidemics in sub-Saharan Africa: appraisal of the WHO strategy. Bull World Health Organ 2002;80:342-9. PMID: 12077608
13. Campagne G, Garba A, Schuchat A, Boulanger D, Plikaytis B, Ousseini M et al. Response to conjugate Haemophilus influenza B vaccine among infants in Niamey, Niger. Am J Trop Med Hyg 1998;59:837-42. PMID: 9840608
14. da Silva A, Parent du Châtelet I, Beckr Gaye A, Dompnier J, Seck I. Microeconomic evaluation of a mass preventive immunisation campaign against meningococcal meningitis and yellow fever in Senegal in 1997. Sante 2003;13:215-23. PMID: 15047438
15. Chippaux JP, Soula G, Campagne G, Rey M. Optimizing the response to epidemics of meningococcal meningitis: report of a workshop of experts at CERMES (Niamey, Niger, 12th to 14th January 2998). Sante 1998;8:245-8.
16. de Chabalier F, Hassane A, Chippaux J. Evaluation of surveillance thresholds for prediction of meningitis epidemics using ongoing surveillance data at the district level, in Niger. Trans R Soc Trop Med Hyg 2000;94:251-2. doi: 10.1016/S0035-9203(00)90308-0 PMID: 10974990
17. Chippaux JP, Debois H, Saliou P. Critical review of control strategies for meningococcal meningitis epidemics in sub-Saharan Africa. Bull Soc Pathol Exot 2002;95:37-44. PMID: 12012963
18. Kaninda AV, Belanger F, Lewis R, Batchassi E, Aplogan A, Yakoua Y et al. Effectiveness of incidence thresholds for detection and control of meningococcal meningitis epidemics in northern Togo. Int J Epidemiol 2000;29:933-40. doi: 10.1093/ije/29.5.933 PMID: 11034981
19. Miller MA, Wenger J, Rosenstein N, Perkins B. Evaluation of meningococcal meningitis vaccination strategies for the meningitis belt in Africa. Pediatr Infect Dis J 1999;18:1051-9. doi: 10.1097/00006454-199912000-00005 PMID: 10608623
20. Miller MA, Shahab CK. Review of the cost effectiveness of immunisation strategies for the control of epidemic meningococcal meningitis. Pharmacoeconomics 2005;23:333-43. doi: 10.2165/00019053-200523040- 00004 PMID: 15853434
21. Chanteau S, Rose A, Djibo S, Nato F, Boisier P. Biological diagnosis of meningococcal meningitis in the African meningitis belt: current epidemic strategy and new perspectives. Vaccine 2007;25(Suppl 1):A30-6. doi: 10.1016/j.vaccine.2007.04.037 PMID: 17517452
22. Rose AM, Gerstl S, Mahamane AE-H, Sidikou F, Djibo S, Bonte L et al. Field evaluation of two rapid diagnostic tests for Neisseria meningitidis serogroup A during the 2006 outbreak in Niger. PLoS ONE 2009;4:e7326. doi: 10.1371/journal.pone.0007326 PMID: 19802392
23. Boisier P, Mahamane A, Hamidou A, Sidikou F, Djibo S, Nato F et al. Field evaluation of rapid diagnostic tests for meningococcal meningitis in Niger. Trop Med Int Health 2009;14:111-7. doi: 10.1111/j.1365-3156.2008.02192.x PMID: 19017310
24. Thomson MC, Molesworth AM, Djingarey MH, Yameogo KR, Belanger F, Cuevas LE. Potential of environmental models to predict meningitis epidemics in Africa. Trop Med Int Health 2006;11:781-8. doi: 10.1111/j.1365- 3156.2006.01630.x PMID: 16771998
25. Yaka P, Sultan B, Broutin H, Janicot S, Philippon S, Fourquet N. Relationships between climate and year-to-year variability in meningitis outbreaks: a case study in Burkina Faso and Niger. Int J Health Geogr 2008;7:34. doi: 10.1186/1476-072X-7-34 PMID: 18597686
26. Jackou-Boulama M, Michel R, Ollivier L, Meynard J, Nicolas P, Boutin J. Correlation between rainfall and meningococcal meningitis in Niger. Med Trop (Mars) 2005;65:329-33. French. PMID: 16548483
27. Wilkinson DJ. Bayesian methods in bioinformatics and computational systems biology. Brief Bioinform 2007;8:109-16. doi: 10.1093/bib/bbm007 PMID: 17430978
28. Briggs AH. A Bayesian approach to stochastic cost-effectiveness analysis. Health Econ 1999;8:257-61. doi: 10.1002/(SICI)1099- 1050(199905)8:3<257::AID-HEC427>3.0.CO;2-E PMID: 10348420
29. Lavindrasana J, Cohen G, Depeursinge A, Müller H, Meyer R, Geissbuhler A. Clinical data mining: a review. In: IMIA yearbook of medical informatics 2009. Geneva: International Medical Informatics Association; 2009. pp. 121-33.
Submitted: 10 January 2011
Revised version received: 5 December 2011
Accepted: 5 December 2011
Published online: 20 January 2012
Competing interests: None declared.
* Correspondence to A Beresniak (e-mail: aberesniak@datamining-international.com).