Evidence of multiple mechanisms providing carbamate and organophosphate resistance in field An. gambiae population from Atacora in Benin
© Aïkpon et al.; licensee BioMed Central Ltd. 2014
Received: 28 September 2014
Accepted: 24 November 2014
Published: 2 December 2014
Insecticide resistance in Anopheles gambiae s.l is a major concern to malaria vector control programmes. In West Africa, resistance is mainly due to target–site insensitivity arising from a single point mutation. Metabolic-based resistance mechanisms have also been implicated and are currently being investigated in west Africa. The aim of this study is to better understand the origins of carbamate and organophosphate resistance in An. gambiae population from Atacora, Benin in West Africa.
Anopheles mosquitoes were reared from larvae collected in two districts (Kouandé and Tanguiéta) of the Atacora department. Mosquitoes were then exposed to WHO impregnated papers. Four impregnated papers were used: carbamates (0.1% bendiocarb, 0.1% propoxur) and organophosphates (0.25% pirimiphos methyl, 1% fenitrothion). PCR assays were run to determine the members of the An. gambiae complex, as well as phenotypes for insensitive acetylcholinesterase (AChE1). Biochemical assays were also carried out to detect any increase in the activity of enzyme typically involved in insecticide metabolism (oxidase, esterase and glutathion-S-transferase).
769 female of An. gambiae mosquitoes from Kouandé and Tanguiéta were exposed to bendiocarb, propoxur, pirimiphos methyl and fenitrothion. Bioassays showed resistance with low mortality to bendiocarb (78.57% to 80.17%), propoxur (77.21% to 89.77%), and fenitrothion (89.74% to 92.02%). On the other hand, the same populations of An. gambiae from Kouandé and Tanguiéta showed high susceptibility to pirimiphos methyl with recorded mortality of 99.02% and 100% respectively. The low rate of ace-1R allele frequency (3.75% among survivors and 0.48% among dead) added to the high proportion of homozygous susceptible specimens which survived the WHO bioassays (8/28), suggest that the ace-1 mutation could not entirely explain Anopheles gambiae resistance to carbamate and organophosphate. Biochemical assays suggest that resistance in this population is mediated by metabolic resistance with elevated level of GST, MFO and NSE compared to a susceptible strain An. gambiae Kisumu.
Anopheles gambiae populations resistance from Atacora is multifactorial and includes target-site mutation and metabolic mechanism. The co-implication of both resistance mechanisms in An. gambiae s.l may be a serious obstacle for the future success of malaria control operations based on LLINs and IRS.
Malaria is a main cause of mortality and morbidity and Anopheles gambiae is one of the major vectors of this desease in sub-Saharan Africa . The current effective vector control tools include the use of Long Lasting Insecticide Nets (LLIN) and Indoor Residual Spraying (IRS) . The success of such interventions requires a good knowledge of vector populations particulary their susceptibility status to the main insecticides used for such control programmes in order to detect and monitor resistance to these insecticides.
Pyrethroids are the only group of insecticides currently recommended for net treatment and the other classes of insecticides (organochlorine, carbamate and organophosphate) are applied for IRS ,. The main problem with LLINs and IRS is the development of insecticide resistance, particularly pyrethroid-resistance by several populations of An. gambiae-. Recently, the emergence of resistance in populations of An. gambiae to common classes of insecticides used in public health has been reported in many countries in Africa, including Côte d’Ivoire ,, Kenya , Benin ,, Niger , Burkina Faso , Mali , Nigeria ,, South Africa  and Cameroon . More recently in Benin, several studies pointed out An. gambiae resistance to bendiocarb , and showed that the mean Ace-1 mutation frequency had increased significantly from 2010 to 2012 after two years of IRS campaign using bendiocarb in Atacora region in the northern part of the country ,. However, these studies cleary showed the evidence of bendiocarb resistance in An. gambiae populations but mentioned the ace-1R mutation could not entirely explain the resistance to bendiocarb observed.
The aim of this study is to better understand the origins of carbamate and organophosphate resistance in An. gambiae population from Atacora. We conducted bioassays, molecular and biochemical analysis to update resistance status in An. gambiae and detect other potential mechanisms of resistance that might have evolved.
Anopheles gambiae s.l. larvae were collected in 2 districts and in each district, two villages were selected randomly. At each locality chosen, Anopheline larvae were collected from various natural breeding sites including ground pools, gutters, puddles and abandoned potholes, at the begining of the rainy season during June 2014. Water was scooped using a plastic scoop and poured into small transparent plastic bowls. A strainer was used to sieve and pool together the third and fourth instar larvae in order to have sufficient adult emergence of the same physiological age. The mosquito larvae collected were transported in well labelled plastic bottles to the laboratory of the Centre de Recherche Entomologique de Cotonou, Benin (CREC) where they were maintained at 28 ± 2 C and 72 ± 5% relative humidity. A laboratory susceptible strain of An. gambiae Kisumu was used as a reference strain to compare the susceptibility levels of the field populations.
Insecticide susceptibility tests
Mosquitoes collected were assayed using WHO discriminating dosages with four insecticides of technical grade quality: two carbamates (0.1% propoxur, 0.1% bendiocarb) and two organophosphates (1% fenitrothion, 0.25% pirimiphos méthyl). Four batches of 25 unfed females, aged 2–5 days, were exposed to the diagnostic doses of insecticide treated papers for 60 min at 27 ± 1°C and 80% relative humidity. The twenty-five females of An. gambiae were introduced into each tube and monitored at different time intervals (10, 15, 20, 30, 45, 60 minutes), the number “knocked-down” were recorded. After one hour exposure, mosquitoes were transferred into holding tubes and provided with cotton wool soaked with a 10% honey solution. Batches exposed to untreated papers were used as control. Mortalities were recorded after 24 hours and the susceptibility status of the population was graded according to the WHO protocol . Dead and surviving mosquitoes from this bioassay were kept separately in eppendorf tubes containing silica gel and stored at −20°C for further molecular analysis.
Species identification and PCR detection of Ace-1R mutation
Live and dead specimens of An. gambiae from the bioassay tests were subjected to the An. gambiae species specific PCR assays for species identification ,. The PCR-RFLP diagnostic test was used to detect the presence of G119S mutation (ace.1R gene). Mosquito genomic DNA was amplified using the primers Ex3AGdir 5’ GATCGTGGACACCGTGTTCG3’ and Ex3AGrev 5‘AGGATGGCCCGCTGGAACAG3’ according to . One microlitre of total DNA extracted from a single mosquito was used as a template in a 25 ml PCR reaction containing Taq DNA polymerase buffer, 0.2 mM dNTP and 10 pmol of each primer. The PCR conditions were 94°C for 5 min and then 35 cycles of (94°C for 30 s, 54°C for 30 s and 72°C for 30 s) with a final 5 min extension at 72°C. Fifteen microlitres of PCR product were digested with 5U of AluI restriction enzyme (Promega) in a final volume of 25 ml. The PCR fragments were fractionated on a 2% agarose gel stained with ethidium bromide and visualized under UV light.
Biochemical assays were performed to compare the levels of activity of mixed function oxidases (MFO), non-specific esterases (NSE) using a-naphtyl acetate as a substrate and glutathione S-transferases (GST)  in the An. gambiae s.s. susceptible Kisumu and the field population from Tanguiéta and Kouandé districts. Mosquitoes used for the biochemical analysis have not been exposed to any insecticides prior to the assay.
Detoxifying enzyme activities were measured on single mosquitoes (N = 50) from each test locality, which were stored at 280uC within 24 h from emergence (above). Each mosquito was ground on ice in 200 μl of distilled water and the homogenate was centrifuged at 1,4000 rpm for 2 mins. Two 10 ml replicates of supernatant were transferred into two adjacent wells of a microtitre plate for non-specific esterases (NSEs), gluthatione Stransferases (GSTs). Monooxygenases assays were performed with two 20 ml replicates of supernatant.
Non-specific esterases (NSEs)
Non-specific esterase activity was measured using a-naphtol acetate (aNa) and b-naphtol acetate (bNa). In each replicate well, 90 ml of phosphate buffer (PBS, pH = 6.5) and 100 ml of 0.6 M aNa (or bNa) were added to 10 ml of centrifuged mosquito homogenate. After 30 min incubation, 100 ml of Fast Garnett BC solution (8 g Fast Garnett Salt +10 ml distilled water) was added to stop the reaction. The concentration of the final product was determined at 550 nm as an endpoint calculated from standard curves of a- and b-Naphtol, respectively.
To measure glutathione- S-transferase (GST) activity in mosquitoes, 200 ml of GSH/CDNB working solution (100 ml of an extemporaneous solution of 0.6% weight/volume reduced glutathione in 0.1 M sodium phosphate buffer pH = 6.5 + 0.013 g of 1-chloro-2, 4-dinitrobenzene diluted in 1 ml of 70% methanol) were added to each replicate of mosquito homogenate. The reaction was read at 340 nm immediately as a kinetic assay for 5 minutes. An extinction coefficient of 5.76 mM−1 (corrected for a path length of 0.6 cm) was used to convert absorbance values to moles of product. GST specific activity was reported as the rate of formation of GSH produced in mmol.min−1.mg−1 protein.
Oxydases (Cytochrome P450)
Cytochrome P450 activity was determined using the heme-peroxidase assay according to Brogdon et al. . The assay detects the elevation in the amount of heme, which is then converted into equivalent units of cytochrome P450. Eighty ml of 0.625 M potassium phosphate buffer (pH = 7.2) were added to 20 ml of mosquito homogenate together with 200 ml Tetramethyl Benzidine solution (0.011 g 3,3’,5,5’ Tetramethyl Benzidine in 5 ml of 70% methanol +15 ml sodium acetate buffer 0.25 M pH = 5.0); 25 ml of 3% hydrogen peroxide were then added and the mixture was incubated for 30 min at room temperature. Absorbance was read at 630 nm and values calculated from a standard curve of cytochrome C.
Mortality rate is > 98%: the population was considered fully susceptible
Mortality rates ranged between 90 - 98%: resistance suspected in the population
Mortality rates < 90%, the population was considered resistant to the tested insecticides
Allelic frequencies of G119S mutation were analysed using the version 1.2 of Genepop .
To assess if the mutation frequencies were identical across populations, the test of genotypic differentiation was performed . Biochemical assay data (enzymatic activity per mg of protein) of wild specimens of An. gambiae s.l were compared to the Kisumu susceptible strain using Mann–Whitney non-parametric U- test (Statistica software).
Bendiocarb and organophosphate susceptibility status of Anopheles gambiae s.l
PCR- species identification
One hundred mosquitoes were identidied to Anopheles gambiae species. All PCR analysis identifying An. gambiae s.l species in this study confirmed that all mosquitoes belonging to An. gambiae s.l were An. gambiae.
Involvement of ace-1R mutation in providing carbamate and organophosphate resistance
This study highlighted the diversity of carbamate and organophosphate resistance in Anopheles gambiae s.l. populations from Atacora in Benin. The results presented here show that carbamate and organohosphate resistance in sampled Anopheles gambiae populations is multifactorial and includes target-site mutation and metabolic mechanism in this region.
An. gambiae displayed large variations in resistance levels to carbamates and organophosphates. Although the wild populations were all resistant to bendiocarb, resistance was less marked to propoxur and fenitrothion, at WHO diagnostic concentrations. However, all these populations were very susceptible to pirimiphos-methyl. This resistance of the mosquito population to carbamate and fenitrothion would be due to the strong selective pressure that represents the use of insecticides in households for public health purposes, notably IRS using bendiocarb  and massive quantities of carbamates and organophosphates in agricultural settings in the department of Atacora. Indeed, in the cotton growing areas in Atacora, farmers use huge amounts of insecticides to avoid substantial yield reduction of their crops. Several studies showed that agricultural practices seem to have contributed to the emergence of insecticide resistance in Anopheles populations ,,.
Cross-resistance to organophosphates and carbamates suggests the involvement of their common target-site: AChE-1 . The high proportion of homoygous susceptible specimens, which survived the WHO bioassays, added to the low rate of ace-1R allele frequency suggest that the ace-1 mutation could not entirely explain Anopheles gambiae resistance to carbamate and organophosphate.
In the present study, the low ace-1 mutation frequency in sampled An. gambiae populations associated with the resistance to carbamate and organophosphate strongly supports the involvement of metabolic resistance based on the high activities of typical enzymes NSE, GST and MFO. Similar findings have been reported in Culex quinquefasciatus and Anopheles gambiae in four other sites in Benin , where greater oxidase and esterase activities were observed in C.x quinquefasciatus and An. gambiae, but where ace-1 was absent . The same findings were also reported in the field experimental station of Pito (Cameroun) ,.
Considering the mutiple insecticide resistance mechanisms in An. gambiae s.l population from Atacora, resistance management becomes a particulary important issue. Insecticide rotation or synergistic use of different classes of insecticides to enhance the killing effects are being considered as alternative methods of vector control to slow down the spread of resistance genes.
The present study revealed the simultaneous presence of multiple resistance mechanisms in malaria vector An. gambiae s.l population from Atacora. The co-implication of both metabolic and ace-1 resistance mechanisms in An. gambiae s.l may be a serious obstacle for the future success of malaria control operations based on LLINs and IRS. Altenative innovative vector control tools are therefore urgently needed to manage these complex mechanisms and to complement insecticide-based strategies.
This work was financially supported by PMI (President’s Malaria Initiative) through USAID. We thank the team of CREC for their technical assistance during field work. We also thank the people of Atacora for their collaboration.
- Gillies MT, Coetzee M: A supplement to the anophelinae Africa South of the Sahara (Afrotropical region), Johannesburg, South Africa. S Afr Inst Med Res. 1987, 55: 1-143.Google Scholar
- Beier JC, Keating J, Githure JI, Macdonald MB, Impoinvil DE, Novak RJ: Integrated vector management for malaria control. Malaria J 2008, 7(supp1):54.,Google Scholar
- Guillet P, Chandre F, Mouchet J: L’utilisation des insecticides en santé publique: état et perspectives. Med Mal Infect. 1997, 27: 552-557. 10.1016/S0399-077X(97)80117-5.View ArticleGoogle Scholar
- Devine GJ, Ogusuku E: Adaptability is key when monitoring insecticide resistance.Bull World Health Organ 2009, 87:887.,PubMed CentralView ArticlePubMedGoogle Scholar
- Elissa N, Mouchet J, Riviere F, Meunier JY, Yao K: Resistance of Anopheles gambiae s.s. to pyrethroids in Côte-d’Ivoire. Ann Soc Belge Med Trop. 1993, 73: 291-294.Google Scholar
- Akogbéto M, Yakoubou S: Résistance des vecteurs du paludisme vis-à-vis des pyréthrinoïdes utilisés pour l’imprégnation des moustiquaires au Bénin, Afrique de l’Ouest. Bull Soc Pathol Exot. 1999, 92: 123-130.PubMedGoogle Scholar
- Chandre F, Darriet F, Duchon S, Finot L, Manguin S, Carnevale P, Guillet P: Modifications of pyrethroid effects associated with kdr mutation in Anopheles gambiae. Med Vet Entomol. 2000, 14: 81-88. 10.1046/j.1365-2915.2000.00212.x.View ArticlePubMedGoogle Scholar
- Akogbeto M, Djouaka RF, Kinde-Gazard DA: Screening of pesticide residues in soil water samples from agricultural settings. Malaria J 2006, 5:22.,View ArticleGoogle Scholar
- Ahoua-Alou LP, Koffi AA, Adja MA, Tla E, Kouassi PK, Kone M, Chandre F: Distribution of ace-1R and resistance to carbamates and organophosphates in Anopheles gambiae s.s. populations from Côte d’Ivoire.Malaria J 2010, 9:167.,View ArticleGoogle Scholar
- Vulule JM, Beach RF, Atieli FK, McAllister JC, Brogdon WG, Roberts JM, Mwangi RW, Hawley WA: Elevated oxidase and esterase levels associated with permethrin tolerance in Anopheles gambiae from Kenyan villages using permethrin impregnated nets. Med Vet Entomol. 1999, 13: 239-244. 10.1046/j.1365-2915.1999.00177.x.View ArticlePubMedGoogle Scholar
- Aïzoun N, Aïkpon R, Padonou GG, Oussou O, Oké-Agbo F, Gnanguenon V, Ossè R, Akogbéto M: Mixed-function oxidases and esterases associated with permethrin, deltamethrin and bendiocarb resistance in Anopheles gambiae s.l in the south-north transect Benin, West Africa.Parasit Vectors 2013, 6:223.,PubMed CentralView ArticlePubMedGoogle Scholar
- Yadouleton AW, Asidi A, Djouaka RF, Braïma J, Agossou CD, Akogbeto MC: Development of vegetable farming: a cause of the emergence of insecticide resistance in populations of Anopheles gambiae in urban areas of Benin.Malaria J 2009, 8:103.,View ArticleGoogle Scholar
- Czeher C, Labbo R, Arzika I, Duchemin JB: Evidence of increasing Leu-Phe knockdown resistance mutation in Anopheles gambiae from Niger following a nationwide longlasting insecticide-treated nets implementation.Malaria J 2008, 7:189.,View ArticleGoogle Scholar
- Diabate A, Baldet T, Chandre F, Guiguemde RT, Brengues C, Guillet P, Hemingway J, Hougard JM: First report of the kdr mutation in Anopheles gambiae M form from Burkina Faso, West Africa. Parassitologia. 2002, 44: 157-158.PubMedGoogle Scholar
- Fanello C, Petrarca V, Della Torre A, Santolamazza F, Dolo G, Coulibaly M, Alloueche A, Curtis CG, Toure YT, Coluzzi M: The pyrethroid knock-down resistance gene in the Anopheles gambiae complex in Mali and further indication of incipient speciation within An. gambiae s.s. Insect Mol Biol. 2003, 12: 241-245. 10.1046/j.1365-2583.2003.00407.x.View ArticlePubMedGoogle Scholar
- Awolola TS, Brooke BD, Koekemoer LL, Coetzee M: Resistance of the malaria vector Anopheles gambiae s.s. to pyrethroid insecticides, in south-western Nigeria. Ann Trop Med Parasitol. 2002, 96: 849-852. 10.1179/000349802125002581.View ArticlePubMedGoogle Scholar
- Oduola AO, Idowu ET, Oyebola MK, Adeogun AO, Olojede JB, Otubanjo OA, Awolola TS: Evidence of carbamate resistance in urban populations of Anopheles gambiae s.s. mosquitoes resistant to DDT and deltamethrin insecticides in Lagos, South-Western Nigeria.Parasit Vectors 2012, 5:116.,PubMed CentralView ArticlePubMedGoogle Scholar
- Hargreaves K, Koerkemoer LL, Brooke B, Hunt RH, Mthembu J, Coetzee M: Anopheles funestus resistant to pyrethroid insecticides in South Africa. Med Vet Entomol. 2000, 14: 181-189. 10.1046/j.1365-2915.2000.00234.x.View ArticlePubMedGoogle Scholar
- Etang J, Manga L, Chandre F, Guillet P, Fondjo E, Mimpfoundi R, Toto JC, Fontenille D: Insecticide susceptibility status of Anopheles gambiae s.l. (Diptera: Culicidae) in the Republic of Cameroon. J Med Entomol. 2003, 40: 491-497. 10.1603/0022-2585-40.4.491.View ArticlePubMedGoogle Scholar
- Aïkpon R, Agossa R, Ossè R, Oussou O, Aïzoun N, Oké-Agbo F, Akogbéto M: Bendiocarb resistance in Anopheles gambiae s.l. populations from Atacora department in Benin, West Africa: a threat for malaria vector control.Parasit Vectors 2013, 6:192.,PubMed CentralView ArticlePubMedGoogle Scholar
- Aïzoun N, Aïkpon R, Gnaguenon V, Oussou O, Agossa F, Padonou GG, Akogbéto M: Status of organophosphate and carbamate resistance in Anopheles gambiae sensu lato from the south and north Benin, West Africa.Parasit Vectors 2013, 6:274.,PubMed CentralView ArticlePubMedGoogle Scholar
- Aïkpon R, Aïzoun N, Sovi A, Oussou O, Govoetchan R, Gnanguenon V, Oké-Agbo R, Ossè R, Akogbéto M: Increase of Ace-1 resistance allele in the field population of Anopheles gambiae following a large scale indoor residual spraying (IRS) implementation using bendiocarb in Atacora region in Benin, West Africa. J Cell Anim Biol. 2014, 8 (1): 15-22. 10.5897/JCAB2013.0399.View ArticleGoogle Scholar
- Aïkpon R, Aïzoun N, Ossè R, Oké-Agbo R, Oussou O, Govoetchan R, Sovi A, Akogbéto M: Seasonal variation of Ace-1R mutation in Anopheles gambiae s. l. populations from Atacora region in Benin, West Africa. J Entomol Nematol. 2014, 6 (1): 14-18. 10.5897/JEN2013.0085.View ArticleGoogle Scholar
- Guidelines for Testing Mosquito Adulticides Intended for Indoor Residual Spraying (IRS) and Insecticide Treated Nets (ITNs). 2006, WHO/CDS/NTD/ WHOPES/GCDDP/2006.3, GenevaGoogle Scholar
- Scott J, Brogdon W, Collins F: Identification of single specimens of the Anopheles gambiae complex by PCR. Am J Trop Med Hyg. 1993, 49: 520-529.PubMedGoogle Scholar
- Favia G, Lanfrancotti A, Spanos L, Siden-Kiamos I, Louis C: Molecular characterization of ribosomal DNA polymorphisms discriminating among chromosomal forms of Anopheles gambiae s.s. Insect Mol Biol. 2001, 10: 19-23. 10.1046/j.1365-2583.2001.00236.x.View ArticlePubMedGoogle Scholar
- Weill M, Lutfalla G, Mogensen K, Chandre F, Berthomieu A, Berticat C, Pasteur N, Philips A, Fort P, Raymond M: Insecticide resistance in mosquito vectors. Nature. 2003, 423: 136-137. 10.1038/423136b.View ArticlePubMedGoogle Scholar
- Hemingway J: Techniques to Detect Insecticide Resistance Mechanisms (Field and Laboratory Manual). Document WHO/CDS/CPC/MAL/98.6. 1998, World Health Organ, Geneva, SwitzerlandGoogle Scholar
- Brogdon WG, McAllister JC: Insecticide resistance and vector control. Emerg Infect Dis. 1998, 4: 605-613. 10.3201/eid0404.980410.PubMed CentralView ArticlePubMedGoogle Scholar
- Raymond M, Rousset F: Genepop (version 1.2), population genetics software for exact tests and eucumenicism. J Heredity. 1995, 86: 248-249.Google Scholar
- Goudet J, Raymond M, De Meeüs T, Rousset F: Testing differentiation in diploid populations. Genetics. 1996, 144: 1933-1940.PubMed CentralPubMedGoogle Scholar
- Diabaté A, Baldet T, Chandre F, Akogbeto M, Guiguemde TR, Darriet F, Brengues C, Guillet P, Hemingway J, Small GJ, Hougard JM: The role of agricultural use of insecticides in resistance to pyrethroids in Anopheles gambiae s.l. in Burkina Faso. Am J Trop Med Hyg. 2002, 67: 617-622.PubMedGoogle Scholar
- Weill M, Malcolm C, Chandre F, Mogensen K, Berthomieu A, Marquine M, Raymond M: The unique mutation in Ace-1 giving high insecticide resistance is easily detectable in mosquito vectors. Insect Mol Biol. 2004, 13: 1-7. 10.1111/j.1365-2583.2004.00452.x.View ArticlePubMedGoogle Scholar
- Corbel V, N’Guessan R, Brengues C, Chandre F, Djogbénou L, Martin T, Akogbéto M, Hougard JM, Rowland M: Multiple insecticide resistance mechanisms in Anopheles gambiae and Culex quinquefasciatus from Benin, West Africa. Acta Trop. 2007, 101: 207-216. 10.1016/j.actatropica.2007.01.005.View ArticlePubMedGoogle Scholar
- Etang J, Manga L, Toto JC, Guillet P, Fondjo E, Chandre F: Spectrum of metabolic-based resistance to DDT and pyrethroids in Anopheles gambiae s.l. populations from Cameroon. J Vect Ecol. 2007, 32: 123-133. 10.3376/1081-1710(2007)32[123:SOMRTD]2.0.CO;2.View ArticleGoogle Scholar
- Chouaibou M, Simard F, Chandre F, Etang J, Darriet F, Hougard JM: Efficacy of bifenthrin-impregnated bednets against Anopheles funestus and pyrethroid-resistant Anopheles gambiae in North Cameroon.Malar J 2006, 5:77.,PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.