Skip to main content
Advertisement
  • Loading metrics

Mycobacterium ulcerans Ecological Dynamics and Its Association with Freshwater Ecosystems and Aquatic Communities: Results from a 12-Month Environmental Survey in Cameroon

  • Andrés Garchitorena ,

    andres.garchitorena@ird.fr

    Affiliations UMR MIVEGEC 5290 CNRS - IRD - Université de Montpellier I - Université de Montpellier II, Montpellier, France, Ecole des Hautes Etudes en Santé Publique, Rennes, France, Service de Mycobactériologie, Centre Pasteur du Cameroun, Réseau International des Instituts Pasteur, Yaoundé, Cameroun

  • Benjamin Roche,

    Affiliations Ecole des Hautes Etudes en Santé Publique, Rennes, France, UMMISCO, UMI IRD-UPMC 209, Bondy, France

  • Roger Kamgang,

    Affiliation Service de Mycobactériologie, Centre Pasteur du Cameroun, Réseau International des Instituts Pasteur, Yaoundé, Cameroun

  • Joachim Ossomba,

    Affiliation Service de Mycobactériologie, Centre Pasteur du Cameroun, Réseau International des Instituts Pasteur, Yaoundé, Cameroun

  • Jérémie Babonneau,

    Affiliation ATOMycA, Inserm Avenir Team, CRCNA, Inserm U892, 6299 CNRS and LUNAM, CHU and Université d'Angers, Angers, France

  • Jordi Landier,

    Affiliations Service de Mycobactériologie, Centre Pasteur du Cameroun, Réseau International des Instituts Pasteur, Yaoundé, Cameroun, Institut Pasteur, Unité d'Epidemiologie de Maladies Emergentes, Paris, France

  • Arnaud Fontanet,

    Affiliations Institut Pasteur, Unité d'Epidemiologie de Maladies Emergentes, Paris, France, Conservatoire National des Arts et Métiers, Paris, France

  • Antoine Flahault,

    Affiliations Centre Virchow-Villermé, Descartes School of Medicine, Université Sorbonne Paris-Cité, Paris, France, Global Health Institute, School of Medicine, University of Geneva, Geneva, Switzerland

  • Sara Eyangoh ,

    Contributed equally to this work with: Sara Eyangoh, Jean-François Guégan, Laurent Marsollier

    Affiliation Service de Mycobactériologie, Centre Pasteur du Cameroun, Réseau International des Instituts Pasteur, Yaoundé, Cameroun

  • Jean-François Guégan ,

    Contributed equally to this work with: Sara Eyangoh, Jean-François Guégan, Laurent Marsollier

    Affiliations UMR MIVEGEC 5290 CNRS - IRD - Université de Montpellier I - Université de Montpellier II, Montpellier, France, Ecole des Hautes Etudes en Santé Publique, Rennes, France

  • Laurent Marsollier

    Contributed equally to this work with: Sara Eyangoh, Jean-François Guégan, Laurent Marsollier

    Affiliation ATOMycA, Inserm Avenir Team, CRCNA, Inserm U892, 6299 CNRS and LUNAM, CHU and Université d'Angers, Angers, France

Abstract

Background

Mycobacterium ulcerans (MU) is the agent responsible for Buruli Ulcer (BU), an emerging skin disease with dramatic socioeconomic and health outcomes, especially in rural settings. BU emergence and distribution is linked to aquatic ecosystems in tropical and subtropical countries, especially to swampy and flooded areas. Aquatic animal organisms are likely to play a role either as host reservoirs or vectors of the bacilli. However, information on MU ecological dynamics, both in space and time, is dramatically lacking. As a result, the ecology of the disease agent, and consequently its mode of transmission, remains largely unknown, which jeopardizes public health attempts for its control. The objective of this study was to gain insight on MU environmental distribution and colonization of aquatic organisms through time.

Methodology/Principal Findings

Longitudinal sampling of 32 communities of aquatic macro-invertebrates and vertebrates was conducted from different environments in two BU endemic regions in Cameroon during 12 months. As a result, 238,496 individuals were classified and MU presence was assessed by qPCR in 3,084 sample-pools containing these aquatic organisms. Our study showed a broad distribution of MU in all ecosystems and taxonomic groups, with important regional differences in its occurrence. Colonization dynamics fluctuated along the year, with the highest peaks in August and October. The large variations observed in the colonization dynamics of different taxonomic groups and aquatic ecosystems suggest that the trends shown here are the result of complex ecological processes that need further investigation.

Conclusion/Perspectives

This is the largest field study on MU ecology to date, providing the first detailed description of its spatio-temporal dynamics in different aquatic ecosystems within BU endemic regions. We argue that coupling this data with fine-scale epidemiological data through statistical and mathematical models will provide a major step forward in the understanding of MU ecology and mode of transmission.

Author Summary

Buruli ulcer, caused by the pathogen Mycobacterium ulcerans (MU), is a severe skin disease occurring in tropical and subtropical countries. Strongly associated to freshwater ecosystems and especially swampy and flooded areas, transmission of this bacterium within ecosystems and across multiple aquatic organisms is still an enigma. Here, we studied in depth the temporal and spatial variations of MU presence in freshwater ecosystems and aquatic organisms in two areas of Cameroon where Buruli ulcer is endemic. We found MU widely present across ecosystems and taxonomic groups along the year and we described a general trend for MU persistence in the environment. Moreover, the colonization dynamics of aquatic ecosystems suggest that each kind of ecosystem may have distinct favourable times of the year for MU presence. In addition to setting the scene for a preventive approach for humans based on ecosystem characteristics, this study suggests that MU transmission is the result of complex ecological processes between biotic and environmental factors. Such results call for an integrative approach in order to disentangle the respective contributions of aquatic organisms and environmental conditions on MU presence and persistence in the environment.

Introduction

Mycobacterium ulcerans (MU) is the agent responsible of Buruli ulcer (BU), an emerging human skin disease affecting human populations in tropical and subtropical regions [1]. While effective treatment is available through a combination of rifampicin-streptomycin for small lesions, with additional surgery required in some cases, early access to treatment is often an issue, especially in poor rural areas where most of the disease burden accumulates [2][4]. Absence or delay of treatment may cause irreversible deformity, long-term disabilities, extensive skin lesions, and even severe secondary infections [5]. Public health efforts for disease control require early detection of cases, but MU ecology and the conditions triggering human infection are poorly understood, which undermines our capacity to detect areas at risk.

Buruli ulcer has been present in Cameroon since the first reported cases in 1969 from the Centre Province, in the districts of Akonolinga and Ayos [6]. The highest BU prevalence in this region dominated by tropical rainforest is distributed along the Nyong River basin, where swampy and flooded areas prevail [7]. A second endemic site appeared in 2004 in Bankim (Adamaoua Province), a region at the border with Nigeria in a transition zone between forest and savannah. Within this region, the construction of a dam in 1989 resulted in a large area of flooded land in the district, and BU cases are mostly concentrated between this dam and the Mbam River [8].

Distribution of human cases around the world seems to be closely related to freshwater ecosystems, especially to areas of slow flowing or stagnant waters [9][14]. Furthermore, emergence of cases in many parts of the world has been associated to the creation of swamps and flooded areas either naturally after heavy rains [13] or under the pressure of human action, i.e. construction of dams or irrigation [8], [15]. Micro-aerobic conditions may promote MU growth [16] and genomic analyses suggest that MU has adapted to a restricted environmental niche, possibly an arthropod [17][19]. Favorable conditions in these types of environment are likely to drive MU growth and persistence and may ultimately affect the transmission to human populations. A direct transmission could take place from the environment where MU is present through existing wounds or passive inoculation [20][23]. However, a direct link between the type of ecosystem and MU abundance in the environment has never been shown.

The role of aquatic communities of macro-invertebrates and vertebrates as a fundamental part of the aquatic ecosystem on MU ecology and transmission is also unclear. Following detection of MU in the environment from abiotic, i.e. water, soil [24], [25] and biotic samples, i.e. plants, fishes [26], tadpoles [27], insect larvae [28], snails [29], and water bugs [30], it has been suggested that bacteria present in the aquatic environment (water, plant biofilm, mud, and detritus) could be concentrated by filtering and grazing invertebrates and then be transmitted through predation up to higher levels of the aquatic trophic web [31]. In addition, some specific taxonomic groups could act as keystone species in the transmission of MU within the aquatic ecosystem [32]. Finally, water bugs of the families Belostomatidae and Naucoridae (Order Hemiptera), which are voracious predators of aquatic organisms may get colonized through this trophic web and transmit the bacteria to humans through biting [30], [33][35].

In order to better understand such a complex disease system, it is essential to address its changes over time and space. Freshwater ecosystems are highly dynamic with seasonal variations in abiotic and biotic parameters impacting on aquatic community assemblages and structures [36], [37]. However, comprehensive field studies performed in Africa to date have addressed temporal dynamics but in only one taxonomic order, Hemiptera water bugs [38], or have focused on aquatic communities but neglecting their temporal dimension [39], [40]. As a result, detailed information on temporal dynamics of MU persistence and spread in the whole aquatic community is dramatically lacking.

Here, we address this issue by performing a large-scale sampling of multiple aquatic communities over space and time in two BU endemic areas of Cameroon, Akonolinga and Bankim. This study aims to improve knowledge on MU environmental distribution and colonization of aquatic organisms throughout the year, with two specific objectives:

  1. To compare MU spatio-temporal distribution in various aquatic ecosystems including swamps, flooded areas, rivers and streams, from two regions with distinct environmental characteristics.
  2. To characterize MU colonization of aquatic communities of macro-invertebrates and vertebrates and its temporal dynamics along the year.

Materials and Methods

Study area and sample sites

Between June 2012 and May 2013, periodic sampling of aquatic communities was performed in Akonolinga and Bankim, two regions in Cameroon where BU is endemic [7], [8]. In order to track colonization dynamics, monthly samples were collected in Akonolinga. In addition, sampling was performed every three months in Bankim, allowing a description of a wider range of environmental characteristics (savannah and tropical rainforest). Within each region, selection of survey sites was done in a two-step procedure. Initially, we classified the villages in each region based on (i) BU human prevalence and (ii) surrounding environmental conditions, according to national health data and land cover data respectively. We pre-selected a number of villages that represented a gradient in both of these parameters within each region. In order to evaluate the relevance of these sites for the study, this pre-selection was followed by on-site visits of all water bodies surrounding the villages and discussions with the local population and health authorities (accessibility, land-use change, human use, persistence throughout the year, etc.).

In all, 32 water sites were selected (16 in each region), including a large variety of streams, rivers, swamps and flooded areas. Streams were defined as bodies of water with a current and were clearly confined within a bed of up to 30 m wide. They included both rainforest streams in Akonolinga and rainforest and savannah streams in Bankim. Rivers were larger than streams, and their margin was highly variable depending on the season, being up to several hundred meters wide in periods of intensive rainfall. They included the Nyong and Mfoumou rivers in Akonolinga, but not the Mbam river in Bankim due to very strong currents that prevented appropriate sampling. We considered as swamps all permanent wetlands with stagnant or very slow flowing waters, many of which were created as the result of roads blocking the natural course of a stream. Finally, flooded areas were temporary bodies of stagnant water formed either naturally after heavy rains in flat areas of forest or savannah, or artificially as in the case of the Mapé Dam in Bankim.

Aquatic sampling

Sampling in each region was performed between 8am and 4pm during 5 consecutive days. In order to ensure comparability of the results, identical methods were carried out by the same persons for all sites throughout the study. In each water body, 4 locations were chosen in areas of slow water flow and among the dominant aquatic vegetation. The sample was limited to those places accessible by a person with waders (depth max. 1.50 m). At each location, 5 sweeps were done with a metallic dip net (32×32 cm, 1 mm mesh size) within a surface of 1 m2 and at different depth levels (down to a depth of 1 m). All the material collected was placed into a bucket with water and passed through a 3-layer filter (32×32 cm grid; 20, 5 and 1 mm mesh sizes, respectively) with abundant water. The material in the first two layers was placed in white rectangular basins, and visible aquatic organisms were identified on site, classified and stored separately into tubes with 70% ethanol. The material contained in the last layer, a mixture of plant debris and small invertebrates, was put into 150 ml flasks with 95% ethanol and brought to the laboratory, where identification of all other individuals in the community (larger than 1 mm) was done with the use of a binocular microscope.

Entomological classification and PCR pool design

Aquatic macro-invertebrates were classified down to the family level whenever possible, using taxonomic keys provided in the Guide to the Freshwater Invertebrates of Southern Africa series [41][47] and other relevant literature [48][51]. In order to avoid cross-contamination between samples, all the equipment used in the classification (forceps, basins, gloves, Petri dishes, etc.) was discarded or decontaminated with NaOH 1 M at the end of each sample classification.

Individuals from the same sample were pooled for PCR analysis by groups of aquatic organisms belonging to the same taxonomic group. Two pooling strategies were used to fulfill the purposes of our study. First, for all sites, we tested a total of 6 sample-pools for each month and each site in order to better describe spatio-temporal dynamics of MU presence. For this, we chose the 5 most abundant taxonomic groups in all sites (to allow for comparability of results) plus a sixth group that was different in each site (to gain representation of all groups), and we pooled all individuals of the same group. Second, we chose 10 sites, 5 sites in each region, for which we applied a more in-depth molecular analysis every 3 months in order to have a better characterization of MU presence in taxonomic groups. Within these sites, all individuals of each taxonomic group were distributed in 4 sample-pools, and all taxonomic groups were tested. The same 10 sites were used along the year and this subgroup presented a similar geographical and environmental variability as the larger group of 32 sites. A maximum of 2 g of pool weight was established in order to avoid excessive inhibition during the qPCR analysis. For each sample-pool, composition, number of individuals and weight of the pool were recorded.

DNA extraction and purification from pools of aquatic organisms

Pooled individuals were all ground together and homogenized in 50 mM NaOH solution using Tissue Lyser II (QIAGEN). Tissue homogenates were heated at 95°C for 20 min. DNA from homogenized insect tissues was purified using QIAquick 96 PCR Purification Kit (QIAGEN), according to manufacturer's recommendations. 10% negative controls were included for extraction and purification.

Detection of MU DNA by quantitative PCR

Oligonucleotide primer and TaqMan probe sequences were selected from the GenBank IS2404 sequence [52] and the ketoreductase B (KR) domain of the mycolactone polyketide synthase (mls) gene (Table 1) from the plasmid pMUM001 [52], [53]. QPCR mixtures contained 5 µl of template DNA, 0.3 µM concentration of each primer, 0.25 µM concentration of the probe, and Brilliant II QPCR master Mix Low Rox (Agilent Technologies) in a total volume of 25 µl. Amplification and detection were performed with Thermocycler (Chromo 4, Bio-Rad) using the following program: 1 cycle of 50°C for 2 min, 1 cycle of 95°C for 15 min, 40 cycles of 95°C for 15 s and 60°C for 1 min. DNA extracts were tested at least in duplicates and the 10% negative controls were included in each assay. Quantitative readout assays were set up, based on external standard curve with MU (strain 1G897) DNA serially diluted over 5 logs (from 106 to 102 U/ml). Samples were considered positive only if both the gene sequence encoding the ketoreductase B domain (KR) of the mycolactone polyketide synthase and IS2404 sequence were detected, with threshold cycle (Ct) values strictly <35 cycles.

thumbnail
Table 1. Primers and probes used to detect M. ulcerans DNA sequences by Taq Man real-time PCR.

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

Data analysis

All statistical analyses were conducted using R statistical software, version 2.14.0 [54]. Maps were created using ArcGIS 10.0 and information displayed in them was obtained from the USGS Shuttle Radar Topography Mission (elevation data) [55], IFORA project (hydrographic network) and Institut National de Cartographie du Cameroun (roads). Data on rainfall was obtained from the NASA Tropical Rainfall Measuring Mission [56]. Pearson Chi Square tests were used to compare proportions of positive sample-pools coming from different types of ecosystems and p-values were computed by Monte-Carlo simulation. One-sample proportions tests with continuity correction were used to calculate the confidence intervals of the proportions. Associations of MU colonization dynamics of taxonomic groups or MU colonization of different ecosystems with rainfall patterns were investigated by calculating the cross-correlation of the time series two by two.

Results

Global distribution of M. ulcerans

Distribution of M. ulcerans in aquatic ecosystems.

MU was broadly distributed within both regions, and was found at least once in more than 80% of sites sampled during the year, with different distribution patterns for each region (Figure 1). In Akonolinga, MU was detected in all sites at least once during the year regardless of the geographical location or the type of ecosystem sampled. MU distribution in Bankim was more restricted, with 4 out of 16 sites found negative all year long, notably from streams in the northern part of the region. Overall, the proportion of positive sample-pools (hereafter defined as “pool positivity” or “pool prevalence”) ranged from 0 to 25% in the different sites, with the highest rates distributed along the road in the southern part of Bankim between the Mapé Dam and the Mbam River, and close to the basin of the Nyong river in Akonolinga (in swamps and streams nearby).

thumbnail
Figure 1. M. ulcerans spatial distribution in water bodies sampled in Cameroon from June 2012 to Mai 2013.

Maps show regional distribution of M. ulcerans in water bodies sampled in Bankim (Left) and Akonolinga (Right). Each circle is a site and colors represent the type of ecosystem sampled. The size of the circles varies according to the percentage of pools that were qPCR positive to both KR and IS2404. Inlet figures illustrate a map of Cameroon with the location of Yaoundé, the capital city (dark triangle) and locations of Bankim and Akonolinga (dark squares).

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

Aquatic ecosystems with stagnant waters appeared to be associated with higher MU presence (Figure 2). We found MU in aquatic organisms from all four types of aquatic ecosystems sampled, with an average of 7.7% of positive sample-pools across ecosystems. Overall, positivity rate was 4.9% in rivers, 4.6% in flooded areas, 10.0% in swamps and 6.2% in streams. We found that swamps had significantly higher positivity than all other ecosystems in Bankim, with positivity in swamps 3 and 5 times higher than in streams and flooded areas respectively (χ2 test, p-value <0.0001 for both). However, no significant differences in MU presence were found for any given environment in Akonolinga, although positivity in flooded areas and swamps was slightly higher.

thumbnail
Figure 2. Overall distribution of M.ulcerans in aquatic ecosystems.

Bars represent total proportion of M. ulcerans DNA positive sample-pools from each type of ecosystem in Akonolinga (blue) and Bankim (red). Whiskers indicate 95% Confidence intervals for the proportions. Asterisks represent significant differences in positivity between ecosystems within each region (χ2 test, p-value<0.0001).

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

Distribution of M. ulcerans in the aquatic community.

A total number of 238,496 individuals were collected and classified over the course of the study, 200,918 in Akonolinga and 37,578 in Bankim. According to the pooling strategy described above, 145,255 of those (61%) were distributed in 3,084 sample-pools and analyzed by qPCR. 65 distinctive taxa were identified (Table S1). 85% of the whole aquatic community overall was made up of only 5 taxonomic orders: Coleoptera, Diptera, Ephemeroptera, Odonata and Hemiptera (Table 1). Among these, the most abundant families were Baetidae (18%), Noteridae (12%), Chironomidae (11%) and Hydrophilidae (7%). Aquatic vertebrates (fishes, tadpoles) and semi-aquatic or terrestrial orders (Araneae, Lepidoptera larvae, Collembola) represented 4% and 2%, respectively. Among aquatic ecosystems, water bodies with standing and slow flowing waters had less biodiversity in terms of number of orders, and were dominated by the 5 most abundant orders mentioned above (Table S2). Conversely, streams had higher biodiversity, with a larger proportion of other groups such as Decapoda and Trichoptera.

MU was present in nearly all taxonomic groups of the aquatic community and it was approximately evenly distributed among the whole aquatic community (Table 2). Pool prevalence for most of the groups was between 5–15%. Larvae of the order Lepidoptera had the highest pool prevalence overall (13.6%), followed by Annelida (12.3%) and Hemiptera (11.4%). However, regional differences in MU distribution should be noted: most of positive Lepidoptera and Annelida came from Bankim, where positivity was nearly 3 times higher for both groups than in Akonolinga (20.8 and 17.7% in Bankim compared to 5.0 and 6.7% in Akonolinga, respectively), although these differences were not significant. The lowest pool prevalence among positive groups was found in Acari (2.8%), Mollusca (3.3%) and Araneae (5.6%). Finally, we failed to detect MU only in two taxonomic groups: Trichoptera (89 pools tested, 1,434 individuals) and Collembola (28 pools tested, 79 individuals).

thumbnail
Table 2. Overall abundance and M. ulcerans presence in pools of aquatic vertebrates and macro-invertebrates.

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

Ecological dynamics along the year

Monthly fluctuations of M. ulcerans presence in aquatic ecosystems.

MU was present in aquatic ecosystems nearly all year long. In Akonolinga, where samples were collected every month, MU was only absent in May, and in Bankim we detected MU in all four time steps sampled (every three months). In this section, only the dynamics for the 12 months in Akonolinga are shown (Figures S1 and S2 show the dynamics in Bankim). MU presence fluctuated through time (Figure 3), with changes from 0 to 15% in total pool positivity. The largest peak in pool positivity was found in August and October, and we found a progressive drop in pool positivity from October to February.

thumbnail
Figure 3. Monthly distribution of M. ulcerans positivity rate in sample-pools from aquatic ecosystems in Akonolinga from June 2012 to May 2013.

Values indicate the proportion of pools of aquatic organisms collected from a specific ecosystem that were positive to M.ulcerans at a given month. The solid line in black represents the total trend (all ecosystems); dashed lines represent trends for pools from each type of ecosystem. Missing information for flooded areas in February and March is due to lack of water in those sites, which prevented sampling.

https://doi.org/10.1371/journal.pntd.0002879.g003

Each ecosystem had distinct temporal variations and a favorable time of the year for MU presence (Figure 3). In rivers and flooded areas, MU was absent for a long period of time (4 and 8 months respectively) and then experienced a sudden increase in pool positivity (in April and July respectively). As a result, more than half of positive sample-pools in these ecosystems were found in a specific season, the low rainy season for rivers and the low dry season for flooded areas (Figure 4). In swamps and streams, the seasonal effect was less pronounced with presence of MU most of the year and fluctuations in pool positivity that ranged from 0 to 15% for swamps and to 30% in streams. Over one third of positive sample-pools in swamps and streams were found during the low rainy season and around one third were found in another season (high dry season for swamps and low dry season for streams). Of all ecosystems, only MU positivity dynamics for rivers were correlated to rainfall dynamics (Figure S5).

thumbnail
Figure 4. Seasonal distribution of M. ulcerans positive sample-pools from each type of ecosystem.

Inset figure on top indicates the rainfall patterns in Akonolinga from June 2012 to May 2013 and the cutting of the period sampled into two dry seasons (LD and HD; Rainfall <100 mm) and two rainy seasons (LR and HR; Rainfall >100 mm). Bars indicate the proportion of M. ulcerans positive sample-pools from a given season and ecosystem out of the total number of positive sample-pools from that ecosystem.

https://doi.org/10.1371/journal.pntd.0002879.g004

Temporal dynamics of M. ulcerans presence in the aquatic community.

MU colonization dynamics for the different taxonomic groups were highly variable (Figure 5). Hemipterans were the only group positive during 11 months of the year, whereas the order Coleoptera was repeatedly negative for more than half a year (from November to May). The highest peaks in pool positivity at any given month were for Hemiptera in June (>30%) and for Diptera in August (25%). Pool positivity in other orders was lower than 20% all year long. Out of the 5 orders systematically tested for all sites and months for MU presence, none of their colonization dynamics were correlated to rainfall (Figure S6).

thumbnail
Figure 5. Monthly distribution of M. ulcerans positivity rate in pools of aquatic organisms in Akonolinga from June 2012 to May 2013.

Values indicate the proportion of pools of aquatic organisms belonging to a specific taxon that were positive to M.ulcerans at a given month. Only the 5 most abundant taxonomic orders were systematically tested for all sites and months. The positivity dynamics for the rest of sample-pools are grouped as “others”. The solid line in black represents the total trend (all taxonomic groups); dashed lines represent trends for each taxonomic group.

https://doi.org/10.1371/journal.pntd.0002879.g005

Discussion

Despite the great health and socio-economic burden borne by people affected with BU, little is known about the ecology and mode of transmission of this disease. MU is embedded in an environment that is inherently dynamic, but information on spatio-temporal dynamics of MU persistence and spread is dramatically lacking. The results shown here represent a step forward in the understanding of MU ecology. They provide the first account of MU spatio-temporal dynamics in aquatic communities from a variety of ecosystems within BU endemic regions. We show first that MU is ubiquitous within these regions and can be found in all types of freshwater ecosystems, but swampy areas seem more favorable to MU presence, as demonstrated in Bankim. Then, we confirm that MU is present in nearly all taxonomic groups of the aquatic community, but we show that groups common in streams are minimally colonized. Finally, we demonstrate that MU has distinctive temporal dynamics in each ecosystem and taxonomic group, suggesting that MU occurrence is probably driven by complex ecological interactions between environmental abiotic and biotic factors.

We found that MU presence in Bankim was more restricted to the south of this area, especially between the Mapé Dam and the Mbam River, where BU cases concentrate [57] and more swamps and flooded areas prevail. The construction of the dam has been previously associated to the emergence of cases in the area [8], [22] and proximity to the Mbam River was found to be a risk factor in a case-control study [22]. However, our results suggest that swamps created along the road, rather than the flooded areas created artificially by the dam or naturally near the Mbam after heavy rains, are more favorable to the presence of MU. Swamps are characterized by stagnant waters with low oxygen and high temperatures, which may create conditions favorable to MU growth and specific fauna in which to develop [16][19]. Furthermore, while water level and conditions in flooded areas are highly variable throughout the year, swamps are more stable environments, which could influence the differences observed in these two stagnant ecosystems [58].

In contrast, MU is present everywhere across the Akonolinga region and all environments presented very similar positivity, although the highest positivity concentrated near the basin of the Nyong river. While climate, land cover or human modifications of the environment could be behind these disparate regional distributions, it could also reflect a spread of the bacteria over time. Indeed, it is possible that MU initially persists in the most favorable environments (swamps), as in the case of Bankim where cases have been reported for less than 10 years [8], spreading over time to other environments where water conditions and aquatic communities are less favorable and/or intermittent along the year, as in the case of Akonolinga where MU is endemic for more than 40 years [6]. Flying insects could be responsible of this dissemination as previously suggested [38], [59]. Out of the two taxonomic orders that are both aquatic in adult stage and capable of flight (Coleoptera and Hemiptera), only Hemiptera was found positive in all types of ecosystem. Indeed, this group was found positive to MU in 65% of the sites, more than any other group of the aquatic community (table S4).

MU is present in nearly every group of the aquatic community and no taxonomic group stands among others as the major host carrier of MU. Aquatic vertebrates and invertebrates, as well as semi-aquatic groups, are positive for IS2404 and KR, with similar pool prevalence. This is in line with the idea of multi-host transmission dynamics and more particularly a transmission through ecological webs, where some species can highly contribute to MU transmission without experiencing a significantly larger positivity [32]. Nevertheless, some patterns arise for several specific taxonomic groups. Firstly, the most positive order in terms of pool prevalence are Lepidoptera larvae (caterpillars), an invertebrate with semi-aquatic families mostly living and feeding on riverine aquatic plants [42]. This finding suggests that some aquatic plants might play an important role on MU persistence and development in the aquatic ecosystem or in ecotone areas, and be a source of infection for herbivorous invertebrates. Indeed, some plants could harbor MU in endemic regions [24], [60] and they stimulate its growth under experimental conditions [61]. Secondly, groups of aquatic invertebrates that were found mainly in streams such as Trichoptera and Decapoda are among the groups with the lowest pool prevalence. These findings support the hypothesis that MU might not be well adapted to environmental conditions in this type of aquatic ecosystems.

Regarding the seasonal dynamics, MU is present in freshwater ecosystems and aquatic organisms throughout the year but there are fluctuations both between seasons and within each season, as previously demonstrated for MU colonization of water bugs [38]. The highest peak in positivity appears in August and October (i.e. over the high rainy season), and then decreases progressively throughout the high dry season (November to February). These findings could be consistent with the idea of a run-off of bacteria into the aquatic environment during periods of intensive rainfall, as previously suggested [24], [62]. However, the lack of correlation between rainfall patterns and the dynamics observed for the various ecosystems and taxonomic orders highlights that more complex interactions might take place within the aquatic community. Differences in feeding habits may explain the distinct colonization dynamics of different orders. For instance, while Hemiptera were found positive all year long (except in May), Coleoptera were repeatedly found negative for more than half a year (Figure 5). These two orders share many common features: they have both larval and adult aquatic stages, many are capable of flight, and their abundance dynamics along the year are very similar (Figure S3 and S4). However, while most families of Hemiptera are voracious predators of aquatic organisms (only Corixidae feed on aquatic plants), families of Coleoptera present a large spectrum of feeding habits that include predators, shredders, scrappers, filtering collectors and omnivorous organisms [41], [42]. Laboratory experiments support the idea of a trophic transmission of MU through predation [30], [34], [63], [64] and a mathematical model studying MU prevalence within 27 aquatic communities in Ghana suggested that a transmission through ecological webs is more likely than a purely environmental acquisition from contaminated water [32]. Our results support the hypothesis that biotic interactions may play a role in MU transmission and that MU dynamics could result from a complex interplay between environmental abiotic factors and variations in community assemblages.

We show that important fluctuations in MU positivity take place within each particular ecosystem. For most sites, we checked for the presence of MU in a given month and site by analysing 6 pools of aquatic organisms. This may be insufficient to demonstrate the absence of the bacteria in the ecosystem, since pool positivity overall was lower than 10%. We attempted to increase the chances to detect MU by pooling all individuals of the most abundant taxonomic orders in the aquatic ecosystem, which allowed us to pool and analyze over 60% of the 238,496 individuals sampled without losing comparability of the results. Furthermore, disparate sampling strategies for each region could be behind the differences found between the types of environment for the two study regions. Bankim was only sampled 4 months of the year as opposed to 12 months in Akonolinga. Therefore, we cannot rule out the possibility that sampling in Bankim may have taken place at appropriate times of the year for swamps but not for the other environments in this region. We tried to avoid this by sampling in Bankim at regular intervals (every three months), therefore capturing a maximum of variability along the year.

This study reinforces the idea that MU persists in a wide range of locations [24], [40] and taxonomic groups [28] and the pool positivity rates described here (nearly 10% overall) are consistent with previous studies [8], [28], [38]. This ubiquity of MU and its persistence in the environment throughout the year contrast with the focal distribution and low number of BU human cases in endemic regions. A possible explanation is that while we are likely to be detecting one (or several) of the MU ecovars present in the environment (previously referred to as mycolactone producing mycobacteria), this does not necessarily imply that we are detecting strains of MU with pathogenic potential to cause BU in humans [17], [65][67]. Future studies comparing the strain diversity of environmental and human samples with molecular techniques such as SNP typing [68], [69] could shed some light on this issue. Furthermore, we rely as previous studies on qPCR amplification of KR and IS2404 sequences as an indicator of the presence of MU, which gives no certainty of whether the DNA detected belongs to viable mycobacteria. The lack of an appropriate technique to culture MU from the environment remains a major limitation of fieldwork studies. Nevertheless, qPCR remains the gold standard for environmental studies on MU ecology [8], [27], [38], [40]. An alternative hypothesis is that while the presence of MU in the environment reflects a potential risk for infection, many environmental and socio-economic factors may need to come together to enable MU transmission to humans. Sero-epidemiological studies have shown that a large proportion of the population living in endemic regions have been exposed to MU, but only a small fraction develop the disease [70]. Therefore, MU might only trigger BU disease under certain environmental conditions (a bacterial concentration threshold and/or contact with a competent, infected vector) or in subpopulations in high contact with potential sources of infection and with increased susceptibility to infection (due to immunity, hygiene, etc.).

In conclusion, this study provides for the first time a detailed characterization through space and time of MU presence in two BU endemic regions with distinct environmental conditions. The understanding of MU ecology to date is still limited, especially regarding the conditions that allow this mycobacterium to persist in the environment and be transmitted to humans. Our study attempts to complete previous approaches by sampling multiple aquatic communities over time in order to better understand the influence of aquatic ecosystems on MU presence and its dynamics along the year. The global trend we describe for MU dynamics could be the result of complex ecological processes, with interactions between environmental abiotic and biotic factors that require deeper analysis, something that is beyond the scope of this paper. However, we believe that coupling data produced by such field studies with fine-scale epidemiological data and integrated through statistical and mathematical models could provide a major step forward in the understanding of MU ecology and BU mode of transmission.

Supporting Information

Figure S1.

Monthly distribution of M. ulcerans positivity rate in pools from aquatic ecosystems in Bankim from June 2012 to March 2013. Values indicate the proportion of pools of aquatic organisms collected from a specific ecosystem that were positive to M.ulcerans at a given month. The solid line in black represents the total trend (all ecosystems); Dashed lines represent trends for pools from each type of ecosystem.

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

(TIF)

Figure S2.

Monthly distribution of M. ulcerans positivity rate in pools of aquatic organisms in Bankim from June 2012 to March 2013. Values indicate the proportion of pools of aquatic organisms belonging to a specific taxon that were positive to M.ulcerans at a given month. Only the 5 most abundant taxonomic orders were systematically tested for all sites and months. The positivity dynamics for the rest of pools are grouped as “others”. The solid line in black represents the total trend (all taxonomic groups); Dashed lines represent trends for each taxonomic group.

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

(TIF)

Figure S3.

Abundance dynamics of aquatic organisms in Akonolinga from June 2012 to May 2013. Abundance values are normalized within each group by dividing abundance for a given month by the maximal abundance for that group. The solid line in black represents the total trend (all taxonomic groups); Dashed lines represent trends for each taxonomic group. Only the 5 most abundant orders are represented. The rest of orders are grouped as “others.

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

(TIF)

Figure S4.

Abundance dynamics of aquatic organisms in Bankim from June 2012 to March 2013. Abundance values are normalized within each group by dividing abundance for a given month by the maximal abundance for that group. The solid line in black represents the total trend (all taxonomic groups); Dashed lines represent trends for each taxonomic group. Only the 5 most abundant orders are represented. The rest of orders are grouped as “others.

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

(TIF)

Figure S5.

Temporal cross-correlation of monthly rainfall distribution and M. ulcerans positivity rate in pools from aquatic ecosystems in Akonolinga from June 2012 to May 2013. Vertical bars indicate the strength of the correlation between the two series for a given lag (in months). Horizontal dashed blue lines represent the threshold of statistical significance.

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

(TIF)

Figure S6.

Temporal cross-correlation of monthly rainfall distribution and M. ulcerans positivity rates in pools of aquatic organisms in Akonolinga from June 2012 to May 2013. Vertical bars indicate the strength of the correlation between the two series for a given lag (in months). Horizontal blue dashed lines represent the threshold of statistical significance.

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

(TIF)

Table S1.

Overall abundance of aquatic vertebrates and macro-invertebrates at the lowest classification level achieved. Results are given for Akonolinga (12 months of sampling) and Bankim (4 months of sampling). Abundance indicates total number of individual organisms collected of each taxonomic group.

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

(PDF)

Table S2.

Total and relative abundance of aquatic vertebrates and macro-invertebrates in aquatic ecosystems. Results are given for Akonolinga (12 months of sampling) and Bankim (4 months of sampling). Total abundance indicates total number of individual organisms collected of each taxonomic group in a given ecosystem. Relative abundance (in brackets) indicates the percentage of individuals from each ecosystem belonging to that taxonomic group.

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

(PDF)

Table S3.

Detection of M. ulcerans DNA by qPCR with KR and IS2404 sequences. Results are given for pools of aquatic vertebrates and invertebrates from Akonolinga (12 months of sampling) and Bankim (4 months of sampling). All pools were tested for the KR sequence. KR positive pools where then confirmed by IS2404. Only sample-pools positive to both sequences were considered positive.

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

(PDF)

Table S4.

Distribution of M. ulcerans positive macro-invertebrates and vertebrates over space and time. Values indicate the number of sites and months where a taxonomic group has been found positive to M. ulcerans DNA by both IS2404 and KR out of the total number of sites and months where the group has been tested.

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

(PDF)

Acknowledgments

We are grateful to the staff of Médecins sans Frontières and Dr. Gervais Folefack (Ministry of Health) in Akonolinga, Dr. Philippe Le Gall (IRD) and Dr. Gaëtan Texier (CPC) for their support in the initial stages of the fieldwork. We are thankful to Myriam Benhalima, Afi Leslie Kaiyven, Francis E. Ntsonzock and the staff of the Mycobacteriology service at Centre Pasteur in Cameroon for their invaluable help in different phases of the study. Thanks are due to the IFORA project, and more particularly Dr. Hervé Chevillotte (IRD). We also thank Kevin Carolan and Estelle Marion for critical review of the manuscript and insightful discussions on the subject.

Author Contributions

Conceived and designed the experiments: AG BR SE JFG LM. Performed the experiments: AG RK JO JB. Analyzed the data: AG BR JFG. Contributed reagents/materials/analysis tools: SE JFG LM. Wrote the paper: AG BR RK JO JB JL AFo AFl SE JFG LM.

References

  1. 1. WHO (2008) Buruli ulcer: progress report, 2004–2008. Geneva, Switzerland.
  2. 2. Mulder AA, Boerma RP, Barogui Y, Zinsou C, Johnson RC, et al. (2008) Healthcare seeking behaviour for Buruli ulcer in Benin: a model to capture therapy choice of patients and healthy community members. Trans R Soc Trop Med Hyg 102: 912–920.
  3. 3. Aujoulat I, Johnson C, Zinsou C, Guédénon A, Portaels F (2003) Psychosocial aspects of health seeking behaviours of patients with Buruli ulcer in southern Benin. Trop Med Int Health 8: 750–759.
  4. 4. Muela Ribera J, Peeters Grietens K, Toomer E, Hausmann-Muela S (2009) A word of caution against the stigma trend in neglected tropical disease research and control. PLoS Negl Trop Dis 3: e445.
  5. 5. Debacker M, Aguiar J, Steunou C, Zinsou C, Meyers WM, et al. (2004) Mycobacterium ulcerans disease (Buruli ulcer) in rural hospital, Southern Benin, 1997–2001. Emerg Infect Dis 10: 1391–1398.
  6. 6. Ravisse P, Rocques MC, Le Bourthe F, Tchuembou CJ, Menard JJ (1975) Une affection méconnue au Cameroun, l'ulcère à Mycobactérie. Med Trop 35: 471–474.
  7. 7. Porten K, Sailor K, Comte E, Njikap A, Sobry A, et al. (2009) Prevalence of Buruli ulcer in Akonolinga health district, Cameroon: results of a cross sectional survey. PLoS Negl Trop Dis 3: e466.
  8. 8. Marion E, Landier J, Boisier P, Marsollier L, Fontanet A, et al. (2011) Geographic Expansion of Buruli Ulcer Disease, Cameroon. Emerg Infect Dis 17: 551–553.
  9. 9. Johnson RC, Makoutodé M, Sopoh GE, Elsen P, Gbovi J, et al. (2005) Buruli Ulcer Distribution in Benin. Emerg Infect Dis 11: 500–501.
  10. 10. Brou T, Broutin H, Elguero E, Asse H, Guegan J-F (2008) Landscape diversity related to Buruli ulcer disease in Côte d'Ivoire. PLoS Negl Trop Dis 2: e271.
  11. 11. Duker AA, Carranza EJ, Hale M (2004) Spatial dependency of Buruli ulcer prevalence on arsenic-enriched domains in Amansie West District, Ghana: implications for arsenic mediation in Mycobacterium ulcerans infection. Int J Health Geogr 3: 19.
  12. 12. Wagner T, Benbow ME, Brenden TO, Qi J, Johnson RC (2008) Buruli ulcer disease prevalence in Benin, West Africa: associations with land use/cover and the identification of disease clusters. Int J Health Geogr 7: 25.
  13. 13. Barker DJP (1972) The distribution of Buruli disease in Uganda. Trans R Soc Trop Med Hyg 66: 867–874
  14. 14. Uganda Buruli Group (1971) Epidemiology of Mycobacterium ulcerans infection (buruli ulcer) at Kinyara, Uganda. Trans R Soc Trop Med Hyg 65: 763–775.
  15. 15. Veitch MGK, Johnson PDR, Flood PE, Leslie DE, Street AC, et al. (1997) A large localized outbreak of Mycobacterium ulcerans infection on a temperate southern Australian island. Epidemiol Infect 119: 313–318.
  16. 16. Palomino JC, Obiang a M, Realini L, Meyers WM, Portaels F (1998) Effect of oxygen on growth of Mycobacterium ulcerans in the BACTEC system. J Clin Microbiol 36: 3420–3422.
  17. 17. Doig KD, Holt KE, Fyfe J a M, Lavender CJ, Eddyani M, et al. (2012) On the origin of Mycobacterium ulcerans, the causative agent of Buruli ulcer. BMC Genomics 13: 258.
  18. 18. Käser M, Rondini S, Naegeli M, Stinear T, Portaels F, et al. (2007) Evolution of two distinct phylogenetic lineages of the emerging human pathogen Mycobacterium ulcerans. BMC Evol Biol 7: 177.
  19. 19. Stinear TP, Seemann T, Pidot S, Frigui W, Reysset G, et al. (2007) Reductive evolution and niche adaptation inferred from the genome of Mycobacterium ulcerans, the causative agent of Buruli ulcer. Genome Res 17: 192–200.
  20. 20. Raghunathan PL, Whitney E a S, Asamoa K, Stienstra Y, Taylor TH, et al. (2005) Risk factors for Buruli ulcer disease (Mycobacterium ulcerans Infection): results from a case-control study in Ghana. Clin Infect Dis 40: 1445–1453.
  21. 21. Pouillot R, Matias G, Wondje CM, Portaels F, Valin N, et al. (2007) Risk factors for buruli ulcer: a case control study in Cameroon. PLoS Negl Trop Dis 1: e101.
  22. 22. Landier J, Piam F, Noumen-djeunga B, Sime J, Wantong G, et al. (2011) Adequate Wound Care and Use of Bed Nets as Protective Factors against Buruli Ulcer?: Results from a Case Control Study in Cameroon. PLoS Negl Trop Dis 5: e1392.
  23. 23. Nackers F, Johnson RC, Glynn JR, Zinsou C, Tonglet R, et al. (2007) Environmental and health-related risk factors for Mycobacterium ulcerans disease (Buruli ulcer) in Benin. Am J Trop Med Hyg 77: 834–836.
  24. 24. Williamson HR, Benbow ME, Campbell LP, Johnson CR, Sopoh G, et al. (2012) Detection of Mycobacterium ulcerans in the Environment Predicts Prevalence of Buruli Ulcer in Benin. PLoS Negl Trop Dis 6: e1506.
  25. 25. Roberts B, Hirst R (1997) Immunomagnetic separation and PCR for detection of Mycobacterium ulcerans. J Clin Microbiol 35: 2709–2711.
  26. 26. Eddyani M, Ofori-adjei D, Teugels G, Weirdt D De, Boakye D, et al. (2004) Potential Role for Fish in Transmission of Mycobacterium ulcerans Disease (Buruli Ulcer): an Environmental Study. Appl Environ Microbiol 70: 5679–5681.
  27. 27. Willson SJ, Kaufman MG, Merritt RW, Williamson HR, Malakauskas DM, et al. (2013) Fish and amphibians as potential reservoirs of Mycobacterium ulcerans, the causative agent of Buruli ulcer disease. Infect Ecol Epidemiol 3: 19946.
  28. 28. Williamson HR, Benbow ME, Nguyen KD, Beachboard DC, Kimbirauskas RK, et al. (2008) Distribution of Mycobacterium ulcerans in buruli ulcer endemic and non-endemic aquatic sites in Ghana. PLoS Negl Trop Dis 2: e205.
  29. 29. Marsollier L, Se T, Aubry J, Merritt RW, Andre J Saint, et al. (2004) Aquatic Snails, Passive Hosts of Mycobacterium ulcerans. Society 70: 6296–6298.
  30. 30. Marsollier L, Robert R, Aubry J, Andre J Saint, Kouakou H, et al. (2002) Aquatic Insects as a Vector for Mycobacterium ulcerans. Appl Environ Microbiol 68: 4623–4628.
  31. 31. Merritt RW, Walker ED, Small PLC, Wallace JR, Johnson PDR, et al. (2010) Ecology and transmission of Buruli ulcer disease: a systematic review. PLoS Negl Trop Dis 4: e911.
  32. 32. Roche B, Eric Benbow M, Merritt R, Kimbirauskas R, McIntosh M, et al. (2013) Identifying the Achilles heel of multi-host pathogens: the concept of keystone “host” species illustrated by Mycobacterium ulcerans transmission. Environ Res Lett 8: 045009.
  33. 33. Marsollier L, Aubry J, Coutanceau E, André J-P Saint, Small PL, et al. (2005) Colonization of the salivary glands of Naucoris cimicoides by Mycobacterium ulcerans requires host plasmatocytes and a macrolide toxin, mycolactone. Cell Microbiol 7: 935–943.
  34. 34. Marsollier L, André J-PS, Frigui W, Reysset G, Milon G, et al. (2007) Early trafficking events of Mycobacterium ulcerans within Naucoris cimicoides. Cell Microbiol 9: 347–355.
  35. 35. Marsollier L, Deniaux E, Brodin P, Marot A, Wondje CM, et al. (2007) Protection against Mycobacterium ulcerans lesion development by exposure to aquatic insect saliva. PLoS Med 4: e64.
  36. 36. Florencio M, Gómez-Rodríguez C (2013) Competitive exclusion and habitat segregation in seasonal macroinvertebrate assemblages in temporary ponds. Freshw Sci 32: 650–662.
  37. 37. Angélibert S, Marty P, Céréghino R, Giani N (2004) Seasonal variations in the physical and chemical characteristics of ponds: implications for biodiversity conservation. Aquat Conserv Mar Freshw Ecosyst 14: 439–456.
  38. 38. Marion E, Eyangoh S, Yeramian E, Doannio J, Landier J, et al. (2010) Seasonal and Regional Dynamics of M. ulcerans Transmission in Environmental Context: Deciphering the Role of Water Bugs as Hosts and Vectors. PLoS Negl Trop Dis 4: 10.
  39. 39. Benbow ME, Williamson H, Kimbirauskas R, Mcintosh MD, Kolar R, et al. (2008) Aquatic Invertebrates as Unlikely Vectors of Buruli Ulcer Disease. Emerg Infect Dis 14: 1247–1254.
  40. 40. Benbow ME, Kimbirauskas R, McIntosh MD, Williamson H, Quaye C, et al.. (2013) Aquatic Macroinvertebrate Assemblages of Ghana, West Africa: Understanding the Ecology of a Neglected Tropical Disease. Ecohealth.
  41. 41. Stals R, Moor I de (2007) Guides to the Freshwater Invertebrates of Southern Africa. Volume 10: Coleoptera. Water Research Commision.
  42. 42. De Moor I, Day J, de Moor F (2004) Guide to the Freshwater Invertebrates of Southern Africa. Volume 8: Insecta II: Hemiptera, Megaloptera, Neuroptera, Trichoptera, Lepidoptera and Orthoptera. Water Research Commision.
  43. 43. Day J, Harrison A, de Moor I (2003) Guide to the Freshwater Invertebrates of Southern Africa. Volume 9: Insecta III (Diptera). Water Research Commision.
  44. 44. De Moor I, Day J, de Moor F (2003) Guide to the Freshwater Invertebrates of Southern Africa. Volume 7: Insecta I (Ephemeroptera, Odonata & Plecoptera). Water Research Commision.
  45. 45. Day J, de Moor I (2002) Guide to the Freshwater Invertebrates of Southern Africa. Volume 6: Arachnida and Mollusca (Araneae, Water Mites and Mollusca). Water Research Commision.
  46. 46. Day J, de Moor I (2002) Guide to the Freshwater Invertebrates of Southern Africa. Volume 5: Non-Arthropods. Water Research Commision.
  47. 47. Fay J (2001) Guide to the Freshwater Invertebrates of Southern Africa. Volume 4: Crustacea III: Bathynellacea, Amphipoda, Isopoda, Spelaeogriphacea, Tanaidaceae and Decapoda. Water Research Commision.
  48. 48. Moisan J (2010) Guide d'identification des principaux macroinvertébrés benthiques d'eau douce du Québec, 2010 – Surveillance volontaire des cours d'eau peu profonds. Direction du suivi de l'état de l'environnement, ministère du Développement durable, de l'Environnement et des Parcs.
  49. 49. Scholth CH, Holm E (1985) Insects of Southern Africa. Durban: Butterworths.
  50. 50. Durand JR, Lévêque C (1980) Flore et faune aquatiques de l'Afrique sahélo-soudanienne. Volume 1. ORSTOM.
  51. 51. Durand JR, Lévêque C (1981) Flore et faune aquatiques de l'Afrique sahelo-soudanienne. Volume 2. ORSTOM.
  52. 52. Rondini S, Troll H, Bodmer T, Pluschke G (2003) Development and Application of Real-Time PCR Assay for Quantification of Mycobacterium ulcerans DNA. J Clin Microbiol 41: 4231–4237.
  53. 53. Fyfe J a M, Lavender CJ, Johnson PDR, Globan M, Sievers A, et al. (2007) Development and application of two multiplex real-time PCR assays for the detection of Mycobacterium ulcerans in clinical and environmental samples. Appl Environ Microbiol 73: 4733–4740.
  54. 54. R Development Core Team (2011) R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.
  55. 55. US Geological Survey (n.d.) Shuttle Radar Topography Mission 90m. Available: http://srtm.usgs.gov/index.php. Accessed 20 December 2013.
  56. 56. NASA Goddard Space Flight Center (n.d.) Tropical Rainfall Measuring Mission. Available: http://disc.sci.gsfc.nasa.gov/precipitation/documentation/guide/TRMM_dataset.gd.shtml#datadesaccess. Accessed 1 December 2013.
  57. 57. Bratschi MW, Bolz M, Minyem JC, Grize L, Wantong FG, et al. (2013) Geographic distribution, age pattern and sites of lesions in a cohort of Buruli ulcer patients from the Mapé Basin of Cameroon. PLoS Negl Trop Dis 7: e2252.
  58. 58. Morris A, Gozlan R, Marion E, Marsollier L, Andreou D, et al. (2014) First detection of Mycobacterium ulcerans in environmental samples from South America. PLoS Negl Trop Dis 8: e2660.
  59. 59. Merritt RW, Benbow ME, Small PL (2005) Unraveling an emerging disease associated with disturbed aquatic environments: the case of Buruli ulcer. Front Ecol Environ 3: 323–331.
  60. 60. McIntosh M, Williamson H, Benbow ME, Kimbirauskas R, Quaye C, et al.. (2014) Associations Between Mycobacterium ulcerans and Aquatic Plant Communities of West Africa: Implications for Buruli Ulcer Disease. Ecohealth. doi: 10.1007/s10393-013-0898-3.
  61. 61. Marsollier L, Stinear T, Aubry J, Paul J, Andre S, et al. (2004) Aquatic Plants Stimulate the Growth of and Biofilm Formation by Mycobacterium ulcerans in Axenic Culture and Harbor These Bacteria in the Environment. Society 70: 1097–1103.
  62. 62. Hayman J (1991) Postulated Epidemiology of Mycobacterium ulcerans Infection. Int J Epidemiol 20: 1093–1098.
  63. 63. Mosi L, Williamson H, Wallace JR, Merritt RW, Small PLC (2008) Persistent association of Mycobacterium ulcerans with West African predaceous insects of the family belostomatidae. Appl Environ Microbiol 74: 7036–7042.
  64. 64. Marsollier L, Brodin P, Jackson M, Korduláková J, Tafelmeyer P, et al. (2007) Impact of Mycobacterium ulcerans biofilm on transmissibility to ecological niches and Buruli ulcer pathogenesis. PLoS Pathog 3: e62.
  65. 65. Lavender CJ, Stinear TP, Johnson PDR, Azuolas J, Benbow ME, et al. (2008) Evaluation of VNTR typing for the identification of Mycobacterium ulcerans in environmental samples from Victoria, Australia. FEMS Microbiol Lett 287: 250–255.
  66. 66. Pidot SJ, Asiedu K, Käser M, Fyfe J a M, Stinear TP (2010) Mycobacterium ulcerans and Other Mycolactone-Producing Mycobacteria Should Be Considered a Single Species. PLoS Negl Trop Dis 4: e663.
  67. 67. Tobias NJ, Doig KD, Medema MH, Chen H, Haring V, et al. (2013) Complete genome sequence of the frog pathogen Mycobacterium ulcerans ecovar Liflandii. J Bacteriol 195: 556–564.
  68. 68. Röltgen K, Qi W, Ruf M-T, Mensah-Quainoo E, Pidot SJ, et al. (2010) Single nucleotide polymorphism typing of Mycobacterium ulcerans reveals focal transmission of buruli ulcer in a highly endemic region of Ghana. PLoS Negl Trop Dis 4: e751.
  69. 69. Vandelannoote K, Jordaens K, Bomans P, Leirs H, Durnez L, et al. (2014) Insertion Sequence Element Single Nucleotide Polymorphism Typing Provides Insights into the Population Structure and Evolution of Mycobacterium ulcerans across Africa. Appl Environ Microbiol 80: 1197–1209.
  70. 70. Yeboah-Manu D, Röltgen K, Opare W, Asan-Ampah K, Quenin-Fosu K, et al. (2012) Sero-Epidemiology as a Tool to Screen Populations for Exposure to Mycobacterium ulcerans. PLoS Negl Trop Dis 6: e1460.