- Open Access
Status of pyrethroid resistance in Anopheles gambiae s. s. M form prior to the scaling up of Long Lasting Insecticidal Nets (LLINs) in Adzopé, Eastern Côte d’Ivoire
© Ahoua Alou et al.; licensee BioMed Central Ltd. 2012
- Received: 19 April 2012
- Accepted: 4 December 2012
- Published: 11 December 2012
The growing development of pyrethroid resistance constitutes a serious threat to malaria control programmes and if measures are not taken in time, resistance may compromise control efforts in the foreseeable future. Prior to Long Lasting Insecticidal Nets (LLINs) distribution in Eastern Cote d’Ivoire, we conducted bioassays to inform the National Malaria Control Programme of the resistance status of the main malaria vector, Anopheles gambiae s. s. and the need for close surveillance of resistance.
Larvae of An. gambiae s. s. were collected in two areas of Adzopé (Port-Bouët and Tsassodji) and reared to adults. WHO susceptibility tests with impregnated filter papers were carried out to detect resistance to three pyrethroids commonly used to develop LLINs: permethrin 1%, deltamethrin 0.05% and lambda-cyhalothrin 0.05%. Molecular assays were conducted to detect M and S forms and the L1014F kdr allele in individual mosquitoes.
Resistance, at various degrees was detected in both areas of Adzopé. Overall, populations of An. gambiae at both sites surveyed showed equivalent frequency of the L1014F kdr allele (0.67) but for all tested pyrethroids, there were significantly higher survival rates for mosquitoes from Tsassodji (32–58%) than those from Port-Bouët (3–32%) (p < 0.001), indicating the implication of resistance mechanisms other than kdr alone. During the survey period (May–June) in this forested area of Côte d’Ivoire, An. gambiae s. s. found were exclusively of the M form and were apparently selected for pyrethroid resistance through agricultural and household usage of insecticides.
Prior to LLINs scaling up in Eastern Côte d’Ivoire, resistance was largely present at various levels in An. gambiae. Underlying mechanisms included the high frequency of the L1014F kdr mutation and other unidentified components, probably metabolic detoxifiers. Their impact on the efficacy of the planned strategy (LLINs) in the area should be investigated alongside careful monitoring of the trend in that resistance over time. The need for alternative insecticides to supplement or replace pyrethroids on nets must be stressed.
- Indoor Residual Spray
- Malaria Vector
- National Malaria Control Programme
- Pyrethroid Resistance
Malaria vector control strategies rely heavily upon the use of insecticide treated nets (ITNs) and indoor residual spraying (IRS). Pyrethroids are the most commonly used insecticides for IRS and also are the only compounds currently approved by the World Health Organization Pesticide Evaluation Scheme (WHOPES) for ITNs. Unfortunately, the growing development of pyrethroid resistance threatens to undermine malaria control programmes. Although the epidemiological significance of pyrethroid resistance has yet to be demonstrated, there are indications, at least in experimental huts and rural households, that ITNs are losing their protective power[4, 5].
A well characterized mechanism of pyrethroid resistance in the malaria vector Anopheles gambiae is pyrethroids target site insensitivity in the voltage-gated sodium channel that induces knockdown resistance (kdr)[6, 7]. In An. gambiae s. s. two alternative amino acid substitutions at the same position (L1014F and L1014S) confer resistance. The first mutation, involving a leucine-to-phenylalanine substitution originally found in West Africa is commonly termed L1014F kdr (kdr-w), whereas the latter mutation found in East Africa and characterized by a serine substitution at the same position is termed L1014S kdr (kdr-e). First detected among An. gambiae field populations from Côte d’Ivoire and Burkina-Faso, the L1014Fkdr is now widespread across West Africa and in some parts of Central and Eastern Africa[3, 8, 9]. The L1014F kdr mutation was initially detected in the S molecular form of An. gambiae s. s. but has now been reported in both S and M forms from West and Central Africa[8, 11, 12].
In Côte d’Ivoire, the National Malaria Control Programme (NMCP) strategies are based on effective case management and high coverage of populations with Long Lasting Insecticidal nets (LLINs), particularly for children <5 yrs old and pregnant women. Aided by the Global Fund Initiatives (Round 6), the NMCP implemented in 2010 a free mass distribution of LLINs aimed at reducing malaria morbidity and mortality. However, several studies have demonstrated an increase in the frequency of the knockdown resistance gene mutations in An. gambiae s. s. following a nationwide LLINs implementation[13–15]. Increasing resistance in malaria vectors may have important implications for vector control programmes, mainly when these are based on the scaling up of LLINs or IRS.
Adzopé, a district hyper-endemic for malaria is one of 19 health districts in Côte d’Ivoire targeted for the distribution of LLINs. Malaria vector susceptibility/resistance status to commonly used insecticides is well documented across the country (North-south-West) where LLINs are being scaled up, except the Eastern district of Adzopé.
This paper therefore reports, prior to LLINs distribution, the pyrethroid susceptibility/resistance status and the prevalence of the L1014F kdr allele in the main malaria vector, An. gambiae s. s. found in the area. The study would inform the NMCP of the resistance status of the main malaria vector, An. gambiae s. s. in the area and the need for its close surveillance.
The study was conducted at Adzopé (6°10 N, 3°85 W) a forested district in Eastern Côte d’Ivoire.
Tsassodji, located in the heart of Adzopé city, with no farming practices causing water retention, and;
Port-Bouët, on the outskirts of the city, with irrigated vegetable production and rice field.
Larvae from Tsassodji were mainly collected from temporary breeding sites such as puddles, shallow wells, gutters, footprints, tyre tracks maintained by rainfall during the collection period. By contrast, at Port-Bouët, the larvae were sampled in rice and vegetable areas.
Pyrethroid susceptibility tests
During the raining season, between May and June 2007, An. gambiae larvae were collected from the breeding sites at each location and reared to the adult stage. Four batches of 20–25 unfed female mosquitoes, 2–3 day-old, were exposed to diagnostic concentrations of permethrin 25/75 (1%), deltamethrin (0.05%) and lambda-cyhalothrin (0.05%) filter papers using WHO standard cylinder kits. Mosquitoes exposed to untreated filter papers served as controls. The number of mosquitoes knocked down at regular intervals during the exposure period was scored and time to knock down 50% and 95% of the exposed mosquitoes (KDT50 and KDT95) were determined. At the end of the exposure period, mosquitoes were transferred to holding tubes and given access to sugar food (10% solution). Mortality was scored 24 h after the holding period. A laboratory susceptible An. gambiae strain Kisumu was used as a reference. Abbott’s formula was used to correct the observed mortality when the rates were between 5% and 20%. All specimens were individually kept on silica gel in Eppendorf tubes for molecular analysis.
Genomic DNA was extracted from individual mosquitoes according to Collins et al. and used for polymerase chain reaction (PCR) analysis. Anopheles gambiae s. l. were identified to species according to Scott et al. then An. gambiae s. s. to molecular forms according to Favia et al.. Both live and dead specimens from the bioassay tests were assayed using allele-specific PCR for the detection of the kdr L1014F mutation as per Martinez-Torres et al.. For each site, the PCR was run on an equal number of mosquito samples, randomly selected among the live and dead individuals after the exposure to the pyrethroid-treated papers.
The WHO criteria for discriminating individuals for susceptibility/resistance status were applied: 98–100% mortality indicating susceptibility; 80–97% mortality requiring confirmation of resistance; and <80% mortality suggesting resistance. Knockdown data were analyzed using the PoloPlus 1.0 software (LeOra Software). KDT50 and KDT95 were generated by means of a log-time probit model. Genotype frequencies at the L1014F kdr locus were compared to Hardy-Weinberg expectations by using the exact probability test developed in Genepop 4.0 software. The genotypic differentiation of the L1014F kdr locus in An. gambiae s. s. populations from the two sites was also tested using the Fischer exact test with the same software.
Knockdown times (KDTs) and mortality rates of Anopheles gambiae M form after exposure to diagnostic concentrations of pyrethroids on filter papers
KDT50s for all field populations (Tsassodji and Port-Bouët) increased significantly compared to the baseline susceptible strain Kisumu. With both mosquito populations, KDT50s were significantly higher for permethrin (>60 min) than for deltamethrin (33.8–50.4 min) and lambda-cyhalothrin (44.7–55.9 min). The Resistance Ratios (RR50s) as assessed by the ratios of the knock down times were moderate with permethrin (6.5–7.9-fold) to low with deltamethrin and lambda-cyhalothrin (1.6–2.4-fold). Clearly, deltamethrin and lambda-cyhalothrin showed greater toxicity than permethrin to both An. gambiae populations. The highest levels of resistance to pyrethroids were observed among populations from the Tsassodji area.
Distribution of the kdr-w allele in Anopheles gambiae M form from Port-Bouët and Tsassodji
L1014F kdr mutation
Genotypic differentiation test analysis showed that there was no significant difference between the distribution of the L1014F kdr allele in suburban area with irrigated agricultural practices (Port-Bouët) and central area without agriculture (Tsassodji) (χ2 = 0.00; df = 2; p = 1.00) despite greater tendency of mosquitoes to survive at Tsassodji than Port Bouët in cylinder bioassays. The L1014F kdr gene frequencies were found to be in Hardy-Weinberg equilibrium in both natural populations of An. gambiae s. s. from Port-Bouët (p = 0.13) and Tsassodji (p = 1.00).
The growing development of insecticide resistance constitutes a serious threat to malaria control programmes and if measures are not taken in time, resistance may compromise control effort in the foreseeable future[4, 22]. Monitoring the development of vector resistance in the field prior to the implementation of any malaria vector control initiative is of paramount importance.
This study documented the susceptibility/resistance status of An. gambiae to three important pyrethroids commonly used to treat mosquito nets currently distributed across sub Saharan Africa, including Côte d’Ivoire (permethrin, deltamethrin and lambda-cyhalothrin). The focus sites of interest were two distinct areas of Adzopé, a city in Eastern Côte d’Ivoire receiving LLINs donated by the Global Fund initiative. The study was carried out prior to the nets distribution to inform the National Malaria Control Programme of the resistance status of vectors and the need for close surveillance of the resistance phenomenon.
The results presented here show that An. gambiae s. s. in the forested areas of Port-Bouët and Tsassodji in Adzopé was exclusively of the M-form. This agrees with a previous study conducted in forest areas of southern Côte d’Ivoire. However, the study was conducted only during May-June and we cannot exclude the occurrence of the S form of An. gambiae s. s. during the rest of the year.
Both vector populations have developed various levels of resistance to the three pyrethroids tested. In the bioassays, An. gambiae s. s. populations sampled in suburban Port-Bouët area of Adzopé city where more controlled agricultural practices with irrigation system exist, showed higher vulnerability to pyrethroid deposits compared to samples that were collected in the heart of the city (Tsassodji) with no such agricultural practice.
This bioassay observation contrasted with the molecular results: the L1014F kdr mutation was detected at both sites but the frequencies of the allele were equivalent (0.67). We do not overlook the implication of the kdr mutation in pyrethroid resistance observed in mosquitoes at both sites as this is also supported by the increased in knockdown time (KDT50) relative to the normal Kisumu strain. Higher KDT50 values in field populations of mosquitoes have been suggested to provide an early indication of the involvement of kdr gene in phenotypic resistance[17, 24]. However, with same kdr rate (67%) but phenotypic difference in expression of resistance by bioassay mortality between both An. gambiae s. s. populations suggests the co-existence of both kdr and other mechanisms, probably enzyme detoxifiers such as esterases, monooxygenases or GSTs. Over-expression of cytochrome P450 genes associated with pyrethroid resistance is most common in An. gambiae s. s., sometimes in association with the L1014F kdr allele[25–28]. No further investigation was conducted to detect additional mechanisms to kdr conferring resistance, but one must not preclude any metabolically mediated mechanisms in the pyrethroid resistance observed at Adzopé. Integrated investigations, which allow detection of target sites mutations and metabolic detoxification agents, should be stressed in order to provide a more comprehensive insight into the genetic basis and the mechanisms responsible for the resistance phenotype in these mosquito populations.
The study confirmed the spread of pyrethroid resistance in An. gambiae s. s. first detected in Côte d’Ivoire, and now in all western African countries investigated[8, 30]. This is the first instance of pyrethroid resistance recorded in An. gambiae s. s. from the eastern part of Côte d’Ivoire, particularly in Adzopé. This resistance may be explained by the selection pressure from both agricultural and domestic usage of insecticides. Farmers in Port-Bouët admitted to use pyrethroids for crop protection. With agricultural practices, the amount of insecticides being applied to the environment is greatly increasing, and may have a pronounced effect on the mosquito ecology and resistance[31–34]. The protection measures against mosquito bites at both sites in Adzopé are mainly domestic aerosols and mosquito coils. Port-Bouët, the suburban area is densely populated with lower educated classes compared to the higher living standard of people at Tsassodji with modern housing and wealth. The household use of products may explain the resistance level observed at Tsassodji as previously reported in rural Côte d’Ivoire.
Prior to LLINs selective distribution to pregnant women and children <5 yrs in Adzopé, Côte d’Ivoire, resistance to a range of pyrethroids (deltamethrin, permethrin, lambda-cyhalothrin) commonly used to treat these nets was detected at various levels in An. gambiae. Mechanisms underlying resistance included high frequency of L1014F kdr and other unidentified components, probably metabolic detoxifiers. Their impact on the efficacy of the planned strategy (LLINs) in the area should be investigated alongside careful monitoring of the trend in that resistance over time. The need for alternative insecticides to supplement or replace pyrethroids on nets must be stressed.
Authors thank Aboubacar Koné and JB Assamoi for their technical assistance during mosquito collections and laboratory analysis. They are also grateful to Bamoro Coulibaly for generating the map of the study sites.
- WHO: The World Malaria Report 2011. 2011, Geneva: World Health OrganizationGoogle Scholar
- WHO: The work of the African Network on Vector Resistance to insecticides, 2000–2004. 2005, Brazzaville: World Health Organization Regional Office for AfricaGoogle Scholar
- Ranson H, N’Guessan R, Lines J, Moiroux N, Nkuni Z, Corbel V: Pyrethroid resistance in African anopheline mosquitoes: what are the implications for malaria control?. Trends Parasitol. 2011, 27: 91-98. 10.1016/j.pt.2010.08.004.View ArticlePubMedGoogle Scholar
- N’Guessan R, Corbel V, Akogbeto M, Rowland M: Reduced efficacy of insecticide-treated nets and indoor residual spraying for malaria control in pyrethroid resistance area, Benin. Emerg Infect Dis. 2007, 13: 199-206. 10.3201/eid1302.060631.PubMed CentralView ArticlePubMedGoogle Scholar
- Sharp BL, Ridl FC, Govender D, Kuklinski J, Kleinschmidt I: Malaria vector control by indoor residual insecticide spraying on the tropical island of Bioko, Equatorial Guinea. Malar J. 2007, 6: 52-10.1186/1475-2875-6-52.PubMed CentralView ArticlePubMedGoogle Scholar
- Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Berge JB, Devonshire AL, Guillet P, Pasteur N, Pauron D: Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s. s. Insect Mol Biol. 1998, 7: 179-184. 10.1046/j.1365-2583.1998.72062.x.View ArticlePubMedGoogle Scholar
- Ranson H, Jensen B, Vulule JM, Wang X, Hemingway J, Collins FH: Identification of a point mutation in the voltage-gated sodium channel gene of Kenyan Anopheles gambiae associated with resistance to DDT and pyrethroids. Insect Mol Biol. 2000, 9: 491-497. 10.1046/j.1365-2583.2000.00209.x.View ArticlePubMedGoogle Scholar
- Santolamazza F, Calzetta M, Etang J, Barrese E, Dia I, Caccone A, Donnelly MJ, Petrarca V, Simard F, Pinto J, della Torre A: Distribution of knock-down resistance mutations in Anopheles gambiae molecular forms in West and West-Central Africa. Malar J. 2008, 7: 74-10.1186/1475-2875-7-74.PubMed CentralView ArticlePubMedGoogle Scholar
- Vezenegho SB, Brooke BD, Hunt RH, Coetzee M, Koekemoer LL: Malaria vector composition and insecticide susceptibility status in Guinea Conakry, West Africa. Med Vet Entomol. 2009, 23: 326-334. 10.1111/j.1365-2915.2009.00840.x.View ArticlePubMedGoogle Scholar
- Chandre F, Brengues C, DossouYovo J, Manguin S, Darriet F, Diabaté A, Carnevale P, Guillet P: Current distribution of pyrethroid resistance gene (kdr) in Anopheles gambiae complex from West Africa and further evidence from reproductive isolation of Mopti form. Parassitogia. 1999, 41: 319-322.Google Scholar
- Weill M, Chandre F, Brengues C, Manguin S, Akogbeto M, Pasteur N, Guillet P, Raymond M: The kdr mutation occur in Mopti form of Anopheles gambiae s. s. through introgression. Insect Mol Biol. 2000, 9: 451-455. 10.1046/j.1365-2583.2000.00206.x.View ArticlePubMedGoogle Scholar
- Etang J, Fondjo E, Chandre F, Morlais I, Brengues C, Nwane P, Chouaibou M, Ndjemai H, Simard F: First report of knockdown mutations in the malaria vector Anopheles gambiae from Cameroon. Am J Trop Med Hyg. 2006, 74: 795-797.PubMedGoogle Scholar
- Stump AD, Atieli FK, Vulule JM, Besansky NJ: Dynamics of the pyrethroid knockdown resistance allele in Western Kenyan populations of Anopheles gambiae in response to insecticide-treated bed net trials. Am J Trop Med Hyg. 2004, 70: 591-596.PubMedGoogle 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 long-lasting insecticide-treated nets implementation. Malar J. 2008, 7: 189-10.1186/1475-2875-7-189.PubMed CentralView ArticlePubMedGoogle Scholar
- Trape JF, Tall A, Diagne N, Ndiath O, Ly AB: Malaria morbidity and pyrethroid resistance after the introduction of insecticide-treated bednets and artemisinin-based combination therapies: a longitudinal study. Lancet Infect Dis. 2011, 11: 925-932. 10.1016/S1473-3099(11)70194-3.View ArticlePubMedGoogle Scholar
- Assi SB, Ouatara L, Adja AM, Bla KB, Tchiekoi NB, Lingué N, Ahoua Alou PL, Koffi AA, Yavo W, Koné M, Djaman AJ: Epidemiology of malaria in urban area (Adzopé city, Côte d’Ivoire). Am Eur J Sci Res. 2010, 5: 94-100.Google Scholar
- WHO: Test procedures for insecticide resistance monitoring in malaria vectors, bio-efficacy, and persistence of insecticides on treated surfaces. 1998, Geneva, Switzerland: WHO/CDS/CPC/MAL/98.12Google Scholar
- Collins FH, Finnerty V, Petrarca V: Ribosomal DNA probes differentiate five cryptic species in the Anopheles gambiae s. s. complex. Parasitology. 1988, 30: 231-240.Google Scholar
- Scott JA, Brogdon WG, Collins FH: Identification of single specimens of Anopheles gambiae complex by the polymerase chain reaction. 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
- Raymond M, Rousset F: Genepop (version 1.2), population genetics software for exact tests and eucumenicism. J Heredity. 1995, 86: 248-249.Google Scholar
- Asidi A, N’Guessan R, Akogbeto M, Curtis C, Rowland M: Loss of household protection from use of insecticide-treated nets against pyrethroid-resistant mosquitoes. Benin. Emerg Infect Dis. 2012, 18: 1101-1106.View ArticlePubMedGoogle Scholar
- Ahoua Alou PL, Koffi AA, Adja AM, Tia E, Kouassi KP, Koné M, Chandre F: Distribution of ace-1R and resistance to carbamates and organophosphates in Anopheles gambiae s. s. populations from Côte d'Ivoire. Malar J. 2010, 9: 167-10.1186/1475-2875-9-167.PubMed CentralView ArticlePubMedGoogle Scholar
- Chandre F, Darriet F, Manga L, Akogbeto M, Faye O, Mouchet J, Guillet P: Status of pyrethroid resistance in Anopheles gambiae s. l. Bull World Health Organ. 1999, 77: 230-235.PubMed CentralPubMedGoogle 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 mechanism 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
- Müller P, Warr E, Stevenson BJ, Pignatelli PM, Morgan JC, Steven A, Yawson AE, Mitchell SN, Ranson H, Hemingway J, Paine MJI, Donnelly MJ: Field-caught permethrin-resistant Anopheles gambiae overexpress CYP6P3, a P450 that metabolises pyrethroids. PLoS Genet. 2008, 4: e1000286-10.1371/journal.pgen.1000286.PubMed CentralView ArticlePubMedGoogle Scholar
- Djouaka RF, Bakare AA, Coulibaly ON, Akogbéto MC, Ranson H, Hemingway J, Strode C: Expression of the cytochrome P450s, CYP6P3 and CYP6M2 are significantly elevated in multiple pyrethroid resistant populations of Anopheles gambiae s. s. from Southern Benin and Nigeria. BMC Genomics. 2008, 9: 538-10.1186/1471-2164-9-538.PubMed CentralView ArticlePubMedGoogle Scholar
- Mitchell SN, Stevenson BJ, Müller P, Wilding CS, Egyir-Yawson A, Field SG, Hemingway J, Paine MJI, Ranson H, Donnelly MJ: Identification and validation of a gene causing cross-resistance between insecticide classes in Anopheles gambiae from Ghana. PNAS. 2012, 109: 6147-6152. 10.1073/pnas.1203452109.PubMed CentralView ArticlePubMedGoogle Scholar
- Elissa N, Mouchet J, Rivière F, Meunier JY, Yao K: Resistance of Anopheles gambiaes.s. to pyrethroids in Côte d’Ivoire. Ann Soc Belge Méd Trop. 1993, 73: 291-294.Google Scholar
- Chandre F, Darriet F, Manguin S, Brengues C, Carnevale P, Guillet P: Pyrethroid cross resistance spectrum among population of Anopheles gambiae from Côte d’Ivoire. J Am Mosq Control Assoc. 1999, 15: 53-59.PubMedGoogle Scholar
- Diabaté A, Baldet T, Chandre F, Akogbéto M, Guiguemdé 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
- Chouaïbou M, Etang J, Brévault T, Nwane P, Hinzoumbé CK, Mimpfoundi R, Simard F: Dynamics of insecticide resistance in the malaria vector Anopheles gambiae s. l. from an area of extensive cotton cultivation in Northern Cameroon. Trop Med Int Health. 2008, 13: 476-486. 10.1111/j.1365-3156.2008.02025.x.View ArticlePubMedGoogle Scholar
- Yadouleton AWM, 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. Malar J. 2009, 8: 103-10.1186/1475-2875-8-103.PubMed CentralView ArticlePubMedGoogle Scholar
- Antonio-Nkondjio C, Fossog BT, Ndo C, Djantio BM, Togouet SZ, Awono-Ambene P, Costantini C, Wondji CS, Ranson H: Anopheles gambiae distribution and insecticide resistance in the cities of Douala and Yaoundé (Cameroon): influence of urban agriculture and pollution. Malar J. 2011, 10: 154-10.1186/1475-2875-10-154.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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.