Molecular Characterization of Mycolactone Producing Mycobacteria from Aquatic...

[Full title: Molecular Characterization of Mycolactone Producing Mycobacteria from Aquatic Environments in Buruli Ulcer Non-Endemic Areas in Côte d’Ivoire]

International Journal of Environmental Research and Public Health Article Molecular Characterization of Mycolactone Producing Mycobacteria from Aquatic Environments in Buruli Ulcer Non-Endemic Areas in Côte d’Ivoire Marcellin B. Tano 1,2 , Christelle Dassi 2,3 , Lydia Mosi 2,4,5, *, Marina Koussémon 1 and Bassirou Bonfoh 2 1 Department of Food Sciences and Technology, Université Nangui Abrogoua, Abidjan, 02 BP 801 Abidjan 01, Côte d’Ivoire; tanobroumarcellin@yahoo.fr (M.B.T.); marinakoussemon@yahoo.fr (M.K.) 2 Department of Environment and Health, Centre Suisse de Recherches Scientifiques en Côte d’Ivoire (CSRS), Adiopodoumé, 01 BP 1303, Abidjan 01, Côte d’Ivoire; christelledassi@gmail.com (C.D.); bassirou.bonfoh@csrs.ci (B.B.) 3 UFR Biosciences, Université Félix Houphouët Boigny, Abidjan, 01 BP V 34 Abidjan 01, Côte d’Ivoire 4 Department of Biochemistry, Cell and Molecular Biology, University of Ghana, Legon, P. O. Box LG 54, Legon, Ghana 5 West African Centre for Cell Biology of Infectious Pathogens, University of Ghana, Legon, P. O. Box LG 54, Legon, Ghana * Correspondence: lmosi@hotmail.com or lmosi@ug.edu.gh; Tel.: +233-54-0890352 Academic Editor: Miklas Scholz Received: 3 October 2016; Accepted: 23 January 2017; Published: 11 February 2017 Abstract: Non-tuberculous mycobacteria (NTM), particularly mycolactone producing mycobacteria (MPM), are bacteria found in aquatic environments causing skin diseases in humans like Buruli ulcer (BU). Although the causative agent for BU, Mycobacterium ulcerans has been identified and associated with slow-moving water bodies, the real transmission route is still unknown. This study aimed to characterize MPMs from environmental aquatic samples collected in a BU non-endemic community, Adiopodoumé, in Côte d’Ivoire. Sixty samples were collected in four types of matrices (plant biofilms, water filtrate residues, plant detritus and soils) from three water bodies frequently used by the population. Using conventional polymerase chain reaction (PCR), MPMs were screened for the 16 S ribosomal RNA (rRNA) mycobacterial gene, the IS 2404 insertion sequence, and MPM enoyl reductase (ER) gene. Variable Number Tandem Repeat (VNTR) typing with loci 6, 19, mycobacterial interspersed repetitive unit 1 (MIRU 1) and sequence type 1(ST 1) was performed to discriminate between different MPMs. Our findings showed 66.7%, 57.5% and 43.5% of positivity respectively for 16 S rRNA, IS 2404 and ER. MPM discrimination using VNTR typing did not show any positivity and therefore did not allow precise MPM distinction. Nevertheless, the observed contamination of some water bodies in a BU non-endemic community by MPMs suggests the possibility of pathogen dissemination and transmission to humans. These aquatic environments could also serve as reservoirs that should be considered during control and prevention strategies Keywords: non-tuberculous mycobacteria; mycolactone producing mycobacteria; environment; Buruli ulcer 1. Introduction Non-tuberculous mycobacteria (NTM) are microorganisms which persist in the environment [ 1 ] and are notorious for causing serious opportunistic diseases, particularly in immunocompromised patients [ 2 ]. Mycolactone producing mycobacteria (MPM) are a sub-group of these NTMs which Int. J. Environ. Res. Public Health 2017 , 14 , 178; doi:10.3390/ijerph 14020178 www.mdpi.com/journal/ijerph

Int. J. Environ. Res. Public Health 2017 , 14 , 178 2 of 10 have been identified as pathogens responsible for infections both in humans and animals. This group comprises Mycobacterium ulcerans , rare select strains of M. marinum , M. liflandii , and M. pseudoshotsii [ 3 ]. Buruli ulcer (BU) is a neglected tropical disease endemic in Australia [ 4 ], South East Asia [ 5 ] and West Africa especially in Benin [ 6 ], Ghana [ 7 ] and Côte d’Ivoire [ 8 ], which is caused by M. ulcerans The disease is very focal in its distribution in all reported countries, with defined endemic and non-endemic areas according to the reports from the various national diseases control centres [ 4 , 9 , 10 ] In 2015, 549 new cases were reported in Côte d’Ivoire, which remains one of most affected countries [ 11 , 12 ] M. ulcerans has been detected in slow moving water bodies in endemic regions [ 13 , 14 ]; however, its precise mode of transmission remains unclear. It is therefore important to assess the distribution of MPMs, particularly M. ulcerans , in both endemic and non-endemic communities in order to decipher the transmission mode to humans and identify the environmental reservoirs of these pathogens Adiopodoumé in Côte d’Ivoire is a non-endemic community where some cases of BU have been reported. The hypothesis that infection results from proximity to contaminated water bodies and socio-cultural practices of individuals along slow moving water bodies [ 13 , 15 ] has been proposed The molecular detection of MPMs by polymerase chain reaction (PCR) amplification from environmental samples has been shown as essential for studying their ecology and transmission [ 10 , 16 , 17 ]. This is because several attempts to directly culture the mycobacterium from the environmental have been largely unsuccessful. Other studies have shown that aquatic plant biofilms, water filtrate residues, plant detritus and soil contain MPM DNA [ 10 , 18 , 19 ]. It has also been demonstrated that vegetables such as lettuce consumed by people could be contaminated by environmental mycobacteria [ 18 ]. The 16 S ribosomal RNA (rRNA) mycobacterial gene and the insertion sequence (IS 2404 ) present in M. ulcerans and other MPM species [ 10 , 20 , 21 ] have been amplified in diverse environmental samples collected from endemic water bodies. Currently, PCR identification of the IS 2404 insertion sequence is one of the World Health Organization (WHO) approved methods for the clinical diagnosis of BU in suspected patients [ 22 ]. Furthermore, a sequence encoding the enoyl reductase (ER) domain of MPM has been included in strain differentiation [ 10 ]. Additionally, variable number tandem repeat (VNTR) typing using the loci 6, 19, sequence type 1 (ST 1) and mycobacterial interspersed repetitive unit 1 (MIRU 1) has been also described as an extremely relevant method for distinguishing between MPMs [ 23 ]. In the present study, a systematic collection of 60 environmental samples from Adiopodoumé, a BU non-endemic area in Côte d’Ivoire, was undertaken as part of a global interand transdisciplinary research study which aimed to identify environmental sources of MPM ( M. ulcerans ) infection and risks associated with BU disease. Using PCR-based detection methods, MPMs were detected in three environmental water bodies suggesting their potential role as a source of contamination to the population in Adiopodoumé 2. Materials and Methods 2.1. Study Site Adiopodoumé is a village located in the southeast of Côte d’Ivoire. This village is surrounded by forest vegetation, several ponds, and the Ebrié lagoon. Its climate is subequatorial with an average ambient temperature of 26 ◦ C and an annual precipitation of about 2000 mm. Adiopodoumé represents a non-endemic community where some suspected cases of BU (13 in 2009 [ 24 ]; 3 in 2014 (data from Adiopodoumé Health Centre) have been reported. The community environment has been modified with soil mining and the building of bridges. The population also uses the water bodies (lagoon, ponds) within this environment for their socio-economic activities like washing dishes, laundry, and gardening.

Int. J. Environ. Res. Public Health 2017 , 14 , 178 3 of 10 2.2. Aquatic Environmental Samples Collection Sample collection was performed according to methods previously described by Narh et al. [ 21 ] with slight modifications. Three water bodies were selected (Figure 1 ) based on human use and associated socio-economic activities Int. J. Environ. Res. Public Health 2017 , 14 , 178 3 of 10 2.2. Aquatic Environmental Samples Collection Sample collection was performed according to methods previously described by Narh et al. [21] with slight modifications. Three water bodies were selected (Figure 1) based on human use and associated socio-economic activities. Figure 1. Study sampling site. Three water bodies were selected based on populations’ use and associated socio-cultural activities. L: Ebrié lagoon; E 1: pond 1; E 2: pond 2. The proposed ecosystem of MPMs within Adiopodoumé community was subsequently described (Figure 2). Sixty environmental samples were collected in four types of matrices (plant biofilms, water filtrate, plant detritus and soil) (Figure S 1). Vegetable leaves ( Celosia argentea , Lactuca sativa , Xanthosoma esculenta , Basella alba , Ipomoea batatas , Hibiscus sabdariffa , Corchorus tridens and Abelmoschus esculentus ) as well as aquatic dominant plants were collected as plant samples. All samples were kept cool and transported to the laboratory for preservation at 4 °C until processing. Legend Area of prospection Adiopodoumé Yopougon Lagoon Ebrié Main road Water bodies Figure 1. Study sampling site. Three water bodies were selected based on populations’ use and associated socio-cultural activities. L: Ebrié lagoon; E 1: pond 1; E 2: pond 2 The proposed ecosystem of MPMs within Adiopodoumé community was subsequently described (Figure 2 ). Sixty environmental samples were collected in four types of matrices (plant biofilms, water filtrate, plant detritus and soil) (Figure S 1). Vegetable leaves ( Celosia argentea , Lactuca sativa , Xanthosoma esculenta , Basella alba , Ipomoea batatas , Hibiscus sabdariffa , Corchorus tridens and Abelmoschus esculentus ) as well as aquatic dominant plants were collected as plant samples. All samples were kept cool and transported to the laboratory for preservation at 4 ◦ C until processing.

Int. J. Environ. Res. Public Health 2017 , 14 , 178 4 of 10 Int. J. Environ. Res. Public Health 2017 , 14 , 178 4 of 10 Figure 2. A proposed ecosystem of mycolactone producing mycobacteria in Adiopodoumé. Three water bodies (pond 1, pond 2 and Ebrié lagoon), frequently used by the community of Adiopodoumé for washing dishes, laundry, and gardening activities, were selected for the study. Pond 1 flowed into pond 2 and they both ran to Ebrié lagoon through the Adiopodoumé Km 17 Bridge. These environmental water bodies represent potential sources of infection by mycolactone producing mycobacteria (MPMs). The transmission of MPMs to humans could occur by direct contact with contaminated environmental samples through skin microtraumatism. 2.3. Molecular Characterization of Mycolactone Producing Mycobacteria DNA extraction from environmental samples, as well as gel-based PCR reactions, were performed following the methods previously described by Williamson et al. and Narh et al. [10,21] with slight modifications. Negative controls (sterile water) were included for the DNA extraction. All gel-based PCR reactions were performed in an A 200 gradient thermal cycler (LonGene, Hangzhou, China). All primers used are described in Table 1. Bovine Serum Albumin (40 ng, Inqaba Biotec, Pretoria, South Africa) was also added to environmental samples to relieve PCR inhibition in the amplification of all target loci. Negative (sterile water) and positive ( M. marinum DL 240490 strain) controls were included for each run. DNA extracted from environmental samples was initially screened for Mycobacterium spp. using mycobacteria-specific primers within the 16 S rRNA gene. Briefly, the reaction was performed in a 25 μL reaction containing 1× PCR buffer (Thermo Scientific, Waltham, MA, USA), 2 mM MgCl 2 , 400 μM each of deoxyribonucleotide (Thermo Scientific), 160 nM each of forward and reverse primers (Inqaba Biotec), 1 U DreamTaq DNA polymerase (Thermo Scientific) and 5 μL of genomic DNA. The reaction was cycled at 95 °C for 3 min followed by 35 cycles each of denaturation at 95 °C for 45 s, annealing at 56 °C for 45 s and extension at 72 °C for 45 s. Final extension was at 72 °C for 10 min and reaction held at 4 °C. Seven microliters (7 uL) of PCR products were run on a 1.5% agarose gel (Thermo Scientific), stained with ethidium bromide (Thermo Scientific) and band sizes were estimated with 100 bp DNA ladder (Thermo Scientific). Pond 1: dish Pond 2: gardening The pond 1 flows to the pond 2 Bridge of Km 17 from Adiopodoumé Lagoon Ebrié: laundry Figure 2. A proposed ecosystem of mycolactone producing mycobacteria in Adiopodoumé Three water bodies (pond 1, pond 2 and Ebrié lagoon), frequently used by the community of Adiopodoumé for washing dishes, laundry, and gardening activities, were selected for the study Pond 1 flowed into pond 2 and they both ran to Ebrié lagoon through the Adiopodoumé Km 17 Bridge These environmental water bodies represent potential sources of infection by mycolactone producing mycobacteria (MPMs). The transmission of MPMs to humans could occur by direct contact with contaminated environmental samples through skin microtraumatism 2.3. Molecular Characterization of Mycolactone Producing Mycobacteria DNA extraction from environmental samples, as well as gel-based PCR reactions, were performed following the methods previously described by Williamson et al. and Narh et al. [ 10 , 21 ] with slight modifications. Negative controls (sterile water) were included for the DNA extraction. All gel-based PCR reactions were performed in an A 200 gradient thermal cycler (LonGene, Hangzhou, China) All primers used are described in Table 1 . Bovine Serum Albumin (40 ng, Inqaba Biotec, Pretoria, South Africa) was also added to environmental samples to relieve PCR inhibition in the amplification of all target loci. Negative (sterile water) and positive ( M. marinum DL 240490 strain) controls were included for each run DNA extracted from environmental samples was initially screened for Mycobacterium spp. using mycobacteria-specific primers within the 16 S rRNA gene. Briefly, the reaction was performed in a 25 µ L reaction containing 1 × PCR buffer (Thermo Scientific, Waltham, MA, USA), 2 mM MgCl 2 , 400 µ M each of deoxyribonucleotide (Thermo Scientific), 160 nM each of forward and reverse primers (Inqaba Biotec), 1 U DreamTaq DNA polymerase (Thermo Scientific) and 5 µ L of genomic DNA The reaction was cycled at 95 ◦ C for 3 min followed by 35 cycles each of denaturation at 95 ◦ C for 45 s, annealing at 56 ◦ C for 45 s and extension at 72 ◦ C for 45 s. Final extension was at 72 ◦ C for 10 min and reaction held at 4 ◦ C. Seven microliters (7 uL) of PCR products were run on a 1.5% agarose gel (Thermo Scientific), stained with ethidium bromide (Thermo Scientific) and band sizes were estimated with 100 bp DNA ladder (Thermo Scientific).

Int. J. Environ. Res. Public Health 2017 , 14 , 178 5 of 10 Positive samples were then screened for the detection of non-tuberculous mycobacteria harboring the IS 2404 insertion sequence using an IS 2404 nested-PCR. Briefly, the first reaction was performed in a 25 µ L reaction containing 1 × PCR buffer (Thermo Scientific), 2 mM MgCl 2 , 300 µ M each of deoxyribonucleotide (Thermo Scientific), 700 nM each of forward and reverse primers (Thermo Scientific), 1 U DreamTaq DNA polymerase (Thermo Scientific) and 5 µ L of genomic DNA. The first reaction was cycled at 95 ◦ C for 2 min followed by 40 cycles each of denaturation at 94 ◦ C for 30 s, annealing at 66 ◦ C for 45 s and extension at 72 ◦ C for 1 min. Final extension was at 72 ◦ C for 10 min and reaction held at 4 ◦ C. The second reaction was performed with the same reagent concentrations as above, but 500 nM each of forward and reverse primers (Thermo Scientific) and 1 µ L of PCR product from the first reaction. The second reaction was cycled at 95 ◦ C for 2 min followed by 35 cycles each of, denaturation at 94 ◦ C for 30 s, annealing at 67.7 ◦ C for 45 s and extension at 72 ◦ C for 45 s Final extension was at 72 ◦ C for 10 min and reaction held at 4 ◦ C Samples which were IS 2404 positive were screened for the presence of the ER gene characterizing MPMs [ 10 , 21 , 23 ]. The ER-PCR reaction was also performed with the same reagent concentrations as described for the 16 S rRNA PCR method and primers described in Table 1 . The reaction was cycled at 95 ◦ C for 2 min followed by 40 cycles each of denaturation at 94 ◦ C for 1 min, annealing at 62.5 ◦ C for 1 min and extension at 72 ◦ C for 1 min. Final extension was at 72 ◦ C for 10 min and reaction held at 4 ◦ C MPM positive samples were processed for VNTR typing with four loci; locus 6, locus 19, MIRU 1 and ST 1 as described in previous studies [ 10 , 21 , 23 ] and with the same reagent concentrations as shown above. The reaction was cycled at 95 ◦ C for 2 min followed by 40 cycles each of denaturation at 94 ◦ C for 1 min, annealing at 58.5 ◦ C for 1 min and extension at 72 ◦ C for 1 min. Final extension was at 72 ◦ C for 10 min and reaction held at 4 ◦ C. In the reaction for ST 1, the annealing was set at 63.1 ◦ C for 1 min Length polymorphism was estimated using PCR product sizes following separation on the agarose gel and recommendations by Williamson et al. and Narh et al. [ 10 , 21 ]. Table 1. Specific primers used for identification of non-tuberculous mycobacteria Primer Name Forward and Reverse Primer Sequences Size of Amplicon (bp) References 16 S rRNA PA: 5 0 -AGAGTTTGATCCTGGCTCAG-3 0 MSHA: 5 0 -AAAAAGCGACAAACCTACGAG-3 0 620 [ 25 ] IS 2404 Nested 1 pGp 1: 5 0 -AGGGCAGCGCGGTGATACGG-3 0 pGp 2: 5 0 -CAGTGGATTGGTGCCGATCGAG-3 0 400 [ 5 ] IS 2404 Nested 2 pGp 3: 5 0 -GGCGCAGATCAACTTCGCGGT-3 0 pGp 4: 5 0 -CTGCGTGGTGCTTTACGCGC-3 0 210 [ 5 ] ER LMF: 5 0 -GAGATCGGTCCCGACGTCTAC-3 0 LMR: 5 0 -GGCTTGACTCATGTCACGTAAG-3 0 420 [ 10 ] Locus 6 R-5 0 -GACATCGAAGAGGTGTGCCGTCT-3 0 F-5 0 -GACCGTCATGTCGTTCGATCCTAGT-3 0 variable [ 10 ] Locus 19 R-5 0 -TGGCGACGATCGAGTCTC-3 0 F-5 0 -CCGACGGATGAATCTGTAGGT-3 0 variable [ 10 ] MIRU 1 R-5 0 -GCCCTCGGGAATGTGGTT-3 0 F-5 0 -GCTGGTTCATGCGTGGAAG-3 0 variable [ 10 ] ST 1 R-5 0 -CGCCACCCGCGGACACAGTCG-3 0 F-5 0 -CTGAGGGGATTTCACGACCAG-3 0 variable [ 10 ] ER: enoyl reductase; MIRU 1: mycobacterial interspersed repetitive unit 1; rRNA: ribosomal RNA; ST 1: sequence type 1 2.4. Statistical Analysis Molecular data were recorded and stored in Microsoft Excel (Microsoft, Redmond, WA, USA) Statistical analysis of average proportions of matrices with positivity for mycobacteria markers was performed using R Version 3.1.2 (R Core Team, Vienna, Austria) and analysis of variance (ANOVA) with one factor and three repetitions. The comparison within matrix positivity was performed using the Fischer test and a p value < 0.05 was considered statistically significant.

Int. J. Environ. Res. Public Health 2017 , 14 , 178 6 of 10 3. Results 3.1. Distribution of Mycobacteria in Aquatic Environment from Adiopodoumé The first screening for mycobacteria among the 60 environmental samples showed 40 (66.7%) positive samples for the gene encoding mycobacterial 16 S rRNA with a 620 bp PCR product (Figure S 2) Of the four different matrices, 12/18 plant biofilms (66.7%) were positive while all the 12 water filtrate residues were positive for this gene. Seven positive plant biofilms were from vegetable leaves ( Celosia argentea , Lactuca sativa , Xanthosoma esculenta , Basella alba , Ipomoea batatas , Hibiscus sabdariffa and Corchorus tridens ). Only 8/15 plant detritus and 8/15 soil samples (53.3%) were positive for the gene coding mycobacterial 16 S rRNA (Table 2 ). Of the eight soil samples that tested positive, five were from the water edge and three were 5 m from the edge of the water body Table 2. Distribution of mycobacteria in aquatic environmental matrices from Adiopodoumé Environmental Samples 16 S rRNA Positives (%) Plant biofilms 12/18 (66.7 a ) Water filtrates 12/12 (100 a ) Plant detritus 8/15 (53.3 a ) Soils 8/15 (53.3 a ) a There was no significant difference between the positive average proportions of the four types of environmental matrices collected and the distribution of mycobacteria in environmental samples could be homogeneous; p value (0.41) > 0.05 for 16 S rRNA 3.2. Distribution of Mycolactone Producing Mycobacteria in Aquatic Environment from Adiopodoumé Using IS 2404 Typing Samples that tested positive for mycobacterial DNA were then screened to identify mycolactone producing mycobacteria by the amplification of the IS 2404 insertion sequence. Twenty-three (57.5%) out of the previous 40 positive environmental samples were positive for IS 2404 insertion sequence with a PCR product size of about 210 bp (Figure S 3). Among these positive samples, 8/12 (66.7%) were plant biofilms with four from vegetables leaves ( Celosia argentea, Lactuca sativa , Hibiscus sabdariffa and Corchorus tridens ), 6/12 (50%) water filtrate residues, 5/8 (62.5%) plant detritus and 4/8 (50%) soil samples (Table 3 ). Table 3. Distribution of IS 2404 positive mycolactone producing mycobacteria in aquatic environmental matrices from Adiopodoumé Environmental Samples IS 2404 Positives (%) Plant biofilms 8/12 (66.7 a ) Water filtrates 6/12 (50 a ) Plant detritus 5/8 (62.5 a ) Soils 4/8 (50 a ) a There was no significant difference between the positive average proportions of the four types of environmental matrices collected and the distribution of mycolactone producing mycobacteria in environmental samples could be homogeneous; p -value (0.72) > 0.05 for IS 2404 3.3. Distribution of Mycolactone Producing Mycobacteria in Aquatic Environment from Adiopodoumé Using ER Typing Amplification of a 420 bp product on the macrolide-lincosamide-streptogramin A (MlsA) domain which encodes the ER domain of the pMUM 001 plasmid of mycolactone producing mycobacteria resulted in only 10 positives out of the 23 IS 2404 positive environmental samples (Figure S 4). Two of these ER-positive samples were aquatic plant biofilms and two water filtrate residues from Ebrié Lagoon. In addition, all the five detritus samples that tested positive for IS 2404 were also ER-positive Only one soil sample from the water body edge was ER-positive (Table 4 ).

Int. J. Environ. Res. Public Health 2017 , 14 , 178 7 of 10 Table 4. Distribution of ER-positive mycolactone producing mycobacteria in aquatic environmental matrices from Adiopodoumé Environmental Samples ER-Positives (%) Plant biofilms 2/8 (25 a ) Water filtrates 2/6 (33.3 a ) Plant detritus 5/5 (100 a ) Soils 1/4 (25 a ) a There was no significant difference between the positive average proportions of the four types of environmental matrices collected and the distribution of mycolactone producing mycobacteria in environmental samples could be homogeneous; p -value (0.14) > 0.05 for ER 3.4. Discrimination between Mycolactone Producing Mycobacteria Species from the Adiopodoumé Aquatic Environment The discrimination of ER positive samples using VNTR typing with locus 6, locus 19, MIRU 1 and ST 1 did not reveal any positive sample. Although VNTR typing was not successful in this, the results above confirm the presence of MPMs in environmental samples collected from Adiopodoumé water bodies 4. Discussion Buruli ulcer is a focal disease and in all affected countries, there are well-defined endemic and non-endemic communities. Cases found in non-endemic communities have been linked with a migration of affected persons from endemic areas [ 4 ]. Although the mode of transmission is still unknown, the causative agent, M. ulcerans has been detected in water bodies of endemic communities The present study focused on the detection of MPM, particularly M. ulcerans from environmental samples collected from aquatic environments in a BU non-endemic community The results showed the presence of mycobacteria in 40 (66.7%) samples after PCR amplification of the gene encoding 16 S rRNA. Stinear et al. [ 17 ] showed the usefulness of 16 S rRNA by asserting that this gene would distinguish between M. ulcerans to M. marinum . This study showed that 100% of water samples tested were positive for the gene coding the 16 S rRNA gene of mycobacteria. This implies that water is the major source of these organisms as reported by Williamson et al. and Narh et al. [ 10 , 21 ] in Ghana with plant biofilms, water, detritus and soils contaminated by mycobacteria. Other authors like Brou et al. [ 14 ] in Côte d’Ivoire and Aiga et al. [ 13 ] in Ghana, have shown that risk factors for BU were linked with water This result also suggests that mycobacteria are found in the samples from the environment The presence of these mycobacteria in Adiopodoumé environment could be explained by the fact that in this community, there are many stagnant water bodies which receive inflows from other slow moving water possibly transporting the pathogen. Some authors as Raghunathan et al. [ 15 ] reported that in Ghana slow moving streams were the main sources of mycobacteria. Adiopodoumé is also located in an urbanized area. In a report by Radomski [ 26 ], urban areas were found to be affected by mycobacterial species due to the proximity of surface water We also included another PCR method based on the amplification of the IS 2404 insertion sequence This sequence has been used for the detection of M. ulcerans and other MPMs in environmental samples because of the high copy number of the insertion sequence (IS) element (213 copies) within their plasmid genome [ 17 ]. The amplification of the IS 2404 insertion sequence in the present study revealed that 23 samples were positive for this marker. This result highlights the presence of MPMs in environmental samples collected from Adiopodoumé water bodies and it is supported by previous studies in Ghana and South Australia [ 10 , 17 , 21 , 27 ]. Furthermore, in Côte d’Ivoire, N’Gazoa-Kacou and her team also assessed M. ulcerans DNA presence by IS 2404 real-time PCR in clinical and environmental samples team [ 28 ]. Stinear et al. conferred that IS 2404 -PCR is a preliminary positive test for the identification of M. ulcerans in environmental samples [ 29 ].

Int. J. Environ. Res. Public Health 2017 , 14 , 178 8 of 10 The present study also highlights the presence of MPMs using IS 2 404 insertion sequence detection in vegetable samples Lactuca sativa , Celosia argentea , Hibiscus sabdariffa and Corchorus tridens These results are supported by those of Marsollier et al. [ 30 ] who showed biofilms formed by M. ulcerans on aquatics plants and an IS 2404 positive plant (Scrophulariaceae family) in Côte d’Ivoire We went on to further detect the presence of the ER gene encoding a polyketide synthase domain of the mycolactone producing plasmid pMUM 001 [ 10 ]. The samples tested positive for the IS 2404 insertion sequence were processed for ER-PCR amplification as a confirmatory test for the identification of M. ulcerans and other MPMs ( M. marinum DL , M. liflandii , and M pseudoshottsii ) [ 21 , 29 ]. Our data showed that from the 23 IS 2404 positive samples, 10 (43.5%) tested positive after ER-PCR, confirming that aquatic environments in Adiopodoumé were contaminated by MPMs. These results are supported by those of Williamson et al. and Narh et al. [ 10 , 21 ] in Ghana, where MPMs were detected in environmental samples from BU non-endemic and endemic areas by ER-PCR Our attempts to further discriminate MPMs by amplification of VNTR loci were not successful in any of the ER-positive environmental samples tested. That may be explained by a low quantity of DNA, from tested environmental samples, not sufficient for amplifying these markers and then typing MPMs. An inverse result was found by Williamson et al. and Narh et al. [ 10 , 21 ] in their studies in Ghana where the same VNTR loci were successfully amplified in environmental samples. Thus, this experiment needs to be repeated under different conditions, using a high quantity of DNA for PCR with confirmation by DNA sequencing and/or more sensitive amplification methods like quantitative PCR, in order to identify more precisely the MPMs involved in aquatic environment contamination in Adiopodoumé community 5. Conclusions The screening of mycolactone producing mycobacteria, using gel-based PCR in environmental samples has shown 66.7% of positivity for 16 S rRNA mycobacterial gene in tested samples. From these positive samples, 57.5% were also positive for the IS 2404 insertion sequence and 43.5% showed positivity for the enoyl reductase gene. Thus, the present study underscores the distribution of MPMs in samples collected within an aquatic environment from Adiopodoumé, a historically non-endemic community for Buruli ulcer disease but where some cases were reported these last years This distribution was homogeneous among the different types of samples collected in this environment Furthermore, VNTR loci typing of these samples did not reveal any positivity that allows the discrimination between M. ulcerans and other MPMs. Nevertheless, these preliminary results showed that more attention should be placed on aquatic environments from BU non-endemic area and environmental samples from there which could be sources of mycobacteria infection in humans using these water bodies and living in these areas Supplementary Materials: The following are available online at http://www.mdpi.com/1660-4601/14/2/178/s 1, Figure S 1: Aquatic environmental sample collection, Figure S 2: Polymerase chain reaction profile obtained after amplification of the 16 S ribosomal RNA (rRNA) gene of mycobacteria; Figure S 3: Polymerase chain reaction profile obtained after amplification of the IS 2404 insertion sequence in non-tuberculous mycobacteria, Figure S 4: Polymerase chain reaction profile obtained after amplification of the enoyl reductase gene of mycolactone producing mycobacteria Acknowledgments: We thank local authorities in the study community as well as populations and field work assistant. We also thank the molecular biology laboratory staff of Centre Suisse de Recherches Scientifiques (CSRS) for assistance with molecular work. The authors acknowledge support from the consortium Afrique One “Ecosystem and Population Health: Expanding Frontiers in Health” which is funded by the Wellcome Trust program (WT 087535 MA). We thank the PASRES program (Programme d’Appui Stratégique de la Recherche Scientifique) for awarding Master student’s fellowships during the study Author Contributions: This work was carried out in collaboration between all authors. C.D., L.M., and B.B conceived and designed the study; C.D. and M.B.T. wrote the protocol; M.B.T. managed the literature searches and performed the experiments; M.B.T., C.D., L.M., M.K., and B.B. analyzed the data and wrote the paper; C.D., L.M., and B.B. contributed with reagents/materials/analysis tools.

Int. J. Environ. Res. Public Health 2017 , 14 , 178 9 of 10 Conflicts of Interest: The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results References 1 Kazda, J. The Chronology of Mycobacteria and the Development of Mycobacterial Ecology. In The Ecology of Mycobacteria , 1 st ed.; Springer Netherlands: Dordrecht, The Netherlands, 2009; pp. 1–11 2 N’Guessan, K.; Kouassi, Y.; Bouzid, S.; Ehuie, P.; Koffi, K.; Oniangue, C.; Aka, N.; Dosso, M. Value and limits of microscopy of exudates in Mycobacterium ulcerans cutaneous infection in Côte d’Ivoire Bull. Soc Pathol. Exot 2001 , 94 , 9–10. [ PubMed ] 3 Ranger, B.S.; Mahrous, E.A.; Mosi, L.; Adusumilli, S.; Lee, R.E.; Colorni, A.; Rhodes, M.; Small, P.L. Globally distributed mycobacterial fish pathogens produce a novel plasmid-encoded toxic macrolide, mycolactone F Infect. Immun 2006 , 74 , 6037–6045. [ CrossRef ] [ PubMed ] 4 Merritt, R.W.; Walker, E.D.; Small, P.L.; Wallace, J.R.; Johnson, P.D.; Benbow, M.E.; Boakye, D.A. Ecology and transmission of Buruli ulcer disease: A systematic review PLoS Negl. Trop. Dis 2010 , 4 , e 911. [ CrossRef ] [ PubMed ] 5 Ablordey, A.; Amissah, D.A.; Aboagye, I.F.; Hatano, B.; Yamazaki, T.; Sata, T.; Ishikawa, K.; Katano, H. Detection of Mycobacterium ulcerans by the loop-mediated isothermal amplification method PLoS Negl Trop. Dis 2012 , 6 , e 1590. [ CrossRef ] [ PubMed ] 6 Debacker, M.; Aguiar, J.; Steunou, C.; Zinsou, C.; Meyers, W.M.; Guédénon, A.; Scott, J.T.; Dramaix, M.; Portaels, F Mycobacterium ulcerans Disease (Buruli Ulcer) in Rural Hospital, Southern Benin, 1997–2001 Emerg. Infect. Dis 2004 , 10 , 1391–1398. [ CrossRef ] [ PubMed ] 7 Amofah, G.; Bonsu, F.; Tetteh, C.; Okrah, J.; Asamoa, K.; Asiedu, K.; Addy, J. Buruli ulcer in Ghana: Results of a national case search Emerg. Infect. Dis 2002 , 8 , 167–170. [ CrossRef ] [ PubMed ] 8 Ahoua, L.; Aka, N.; Ekaza, E.; Bouzid, S.; N’Guessan, R.; Dosso, M. Risk factors for Buruli ulcer in Côte d’Ivoire: Results of a case-control study Afr. J. Biotechnol 2009 , 8 , 536–546 9 Wagner, T.; Benbow, M.E.; Brenden, T.O.; Qi, J.; Johnson, R.C. Buruli ulcer disease prevalence in Benin, West Africa: Associations with land use/cover and the identification of disease clusters Int. J. Health Geogr 2008 , 7 , 25. [ CrossRef ] [ PubMed ] 10 Williamson, H.R.; Benbow, M.E.; Nguyen, K.D.; Beachboard, D.C.; Kimbirauskas, R.K.; McIntosh, M.D.; Quaye, C.; Ampadu, E.O.; Boakye, D.; Merritt, R.W.; et al. Distribution of Mycobacterium ulcerans in Buruli ulcer endemic and non-endemic aquatic sites in Ghana PLoS Negl. Trop. Dis 2008 , 2 , e 205. [ CrossRef ] [ PubMed ] 11 Kanga, J.M.; Kacou, D.E. Aspects épidémiologiques de l’ulcère de Buruli en Côte d’Ivoire: Résultats d’une enquête nationale Bull. Soc. Pathol. Exot 2001 , 94 , 46–51. [ PubMed ] 12 Résumé du bilan 2013 du Programme National de Lutte contre l'ulcère de Buruli. Aspect clinique, épidémiologique et thérapeutique ; PNLUB: Abidjan, Côte d’Ivoire, 2014 13 Aiga, H.; Amano, T.; Cairncross, S.; Adomako, J.; Nanas, O.K.; Coleman, S. Assessing water-related risk factors for Buruli ulcer: A case-control study in Ghana Am. J. Trop. Med. Hyg 2004 , 71 , 387–392. [ PubMed ] 14 Brou, T.; Broutin, H.; Elguero, E.; Asse, H.; Guegan, J.F. Landscape diversity related to Buruli ulcer disease in Côte d’Ivoire PLoS Negl. Trop. Dis 2008 , 2 , e 271. [ CrossRef ] [ PubMed ] 15 Raghunathan, P.L.; Whitney, E.A.; Asamoa, K.; Stienstra, Y.; Taylor, T.H., Jr.; Amofah, G.K.; Ofori-Adjei, D.; Dobos, K.; Guarner, J.; Martin, S.; et al. Risk factors for Buruli ulcer disease ( Mycobacterium ulcerans Infection): Results from a case-control study in Ghana Clin. Infect. Dis 2005 , 40 , 1445–1453. [ CrossRef ] [ PubMed ] 16 Ekaza, E.; Kacou-N’douba, A.; Oniangue, N.C.; Ehuie, P.; N’Guessan, K.R.; Aka, N.; Bouzid, S.A.; Faye-Kette, H.; Dosso, M. Apport de l’amplification génique dans la détection de Mycobacterium ulcerans dans les exsudats et les biopsies cutanées en Côte d’Ivoire Bull. Soc. Pathol. Exot 2004 , 97 , 95–96. [ PubMed ] 17 Stinear, T.P.; Seemann, T.; Pidot, S.; Frigui, W.; Reysset, G.; Garnier, T.; Meurice, G.; Simon, D.; Bouchier, C.; Ma, L.; et al. Reductive evolution and niche adaptation inferred from the genome of Mycobacterium ulcerans , the causative agent of Buruli ulcer Genome Res 2007 , 17 , 192–200. [ CrossRef ] [ PubMed ] 18 Yoder, S.; Argueta, C.; Holtzman, A.; Aronson, T.; Berlin, O.G.W.; Tomasek, P.; Glover, N.; Froman, S.; Stelma, G.J. PCR comparison of Mycobacterium avium isolates obtained from patients and foods Appl. Environ. Microbiol 1999 , 65 , 2650–2653. [ PubMed ]

Int. J. Environ. Res. Public Health 2017 , 14 , 178 10 of 10 19 Lavender, C.J.; Stinear, T.P.; Johnson, P.D.; Azuolas, J.; Benbow, M.E.; Wallace, J.R.; Fyfe, J.A. Evaluation of VNTR typing for the identification of Mycobacterium ulcerans in environmental samples from Victoria, Australia FEMS Microbiol. Lett 2008 , 287 , 250–255. [ CrossRef ] [ PubMed ] 20 Fyfe, J.A.; Lavender, C.J.; Johnson, P.D.R.; Globan, M.; Sievers, A.; Azuolas, J.; Stinear, T.P. Development and application of two multiplex real-time PCR assays for the detection of Mycobacterium ulcerans in clinical and environmental samples Appl. Environ. Microbiol 2007 , 73 , 4733–4740. [ CrossRef ] [ PubMed ] 21 Narh, C.A.; Mosi, L.; Quaye, C.; Dassi, C.; Konan, D.O.; Tay, S.C.K.; de Souza, D.K.; Boakye, D.A.; Bonfoh, B. Source Tracking Mycobacterium ulcerans Infections in the Ashanti Region, Ghana PLoS Negl. Trop. Dis 2015 , 9 , e 0003437. [ CrossRef ] [ PubMed ] 22 Portaels, F.; World Health Organization Laboratory Diagnosis of Buruli Ulcer: A Manual for Health Care Providers ; World Health Organization: Geneva, Switzerland, 2014. Available online: http://www.who.int/ iris/handle/10665/111738 (accessed on 17 January 2017) 23 Williamson, H.; Phillips, R.; Sarfo, S.; Wansbrough-Jones, M.; Small, P. Genetic Diversity of PCR-Positive, Culture-Negative and Culture-Positive Mycobacterium ulcerans Isolated from Buruli Ulcer Patients in Ghana PLoS ONE 2014 , 9 , e 88007 24 Aka, N.; Ekaza, E.; Coulibaly-Ngolo, M.D.; Kouadio, K.; Koffi, L.; Coulibaly, B.; Kodia, M.; Ngazoa-Kakou, S.; N’Guessan, K.R.; Yapo-Crézoit, A.; et al. Buruli ulcer in Côte d’Ivoire: a new disease focus identified in Adiopodoumé, Yopougon commune (Abidjan, Km 17). In Abstracts of the annual meeting of the WHO Global Ulcer Initiative, 22–24 March 2010 ; WHO Headquarters: Geneva, Switzerland, 2010; p. 161 25 Hughes, M.S.; Skuce, R.A.; Beck, L.A.; Neill, S.D. Identification of mycobacteria from animals by restriction enzyme analysis and direct DNA cycle sequencing of polymerase chain reaction-amplified 16 S rRNA gene sequences J. Clin. Microbiol 1993 , 31 , 3216–3222 26 Radomski, N. Sources des mycobactéries non tuberculeuses dans les bassins versants. Ph.D. Thesis, Université Paris-Est, Paris, France, 2011 27 Kirschner, P.; Bottger, E.C. Species identification of mycobacteria using rDNA sequencing. In Mycobacteria Protocols (Methods in Molecular Biology) ; Parish, T.S.N., Ed.; Humana Press Inc.: Totowa, NJ, USA, 1998; Volume 101, pp. 349–361 28 Ngazoa-Kakou, E.S.; Ekaza, E.; Aka, N.; Coulibaly-N’Golo, D.; Coulibaly, B.; Dosso, M. Evaluation of real-time PCR for Mycobacterium ulcerans in endemic region in Côte d’Ivoire Afr. J. Microbiol. Res 2011 , 5 , 2211–2216 29 Stinear, T.P.; Pryor, M.J.; Porter, J.L.; Cole, S.T. Functional analysis and annotation of the virulence plasmid pMUM 001 from Mycobacterium ulcerans Microbiology 2005 , 151 , 683–692. [ CrossRef ] [ PubMed ] 30 Marsollier, L.; Stinear, T.; Aubry, J.; Saint André, J.P.; Robert, R.; Legras, P.; Manceau, A.-L.; Audrain, C.; Bourdon, S.; Kouakou, H.; et al. Aquatic Plants Stimulate the Growth of and Biofilm Formation by Mycobacterium ulcerans in Axenic Culture and Harbor These Bacteria in the Environment Appl. Environ. Microbiol 2004 , 70 , 1097–1103. [ CrossRef ] [ PubMed ] © 2017 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).