Next Article in Journal
Epidemiologic Investigation of Two Welder’s Anthrax Cases Caused by Bacillus cereus Group Bacteria: Occupational Link Established by Environmental Detection
Previous Article in Journal
Prescription of Anti-Spike Monoclonal Antibodies in COVID-19 Patients with Resistant SARS-CoV-2 Variants in Italy
Previous Article in Special Issue
High CD169 Monocyte/Lymphocyte Ratio Reflects Immunophenotype Disruption and Oxygen Need in COVID-19 Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bendiocarb and Malathion Resistance in Two Major Malaria Vector Populations in Cameroon Is Associated with High Frequency of the G119S Mutation (Ace-1) and Overexpression of Detoxification Genes

by
Idriss Nasser Ngangue-Siewe
1,2,
Paulette Ndjeunia-Mbiakop
2,3,
Nelly Armanda Kala-Chouakeu
2,4,
Roland Bamou
2,4,
Abdou Talipouo
2,3,
Landre Djamouko-Djonkam
2,4,
John Vontas
5,6,
Konstantinos Mavridis
5,
Jeannette Tombi
3,
Timoléon Tchuinkam
4,
Jean Arthur Mbida-Mbida
1 and
Christophe Antonio-Nkondjio
2,7,*
1
Laboratory of Animal Biology and Physiology, University of Douala, Douala P.O. Box 24157, Cameroon
2
Laboratoire de Recherche sur le Paludisme, Organisation de Coordination pour la lutte contre les Endémies en Afrique Centrale (OCEAC), Yaoundé P.O. Box 288, Cameroon
3
Faculty of Sciences, University of Yaoundé I, Yaoundé P.O. Box 812, Cameroon
4
Vector-Borne Diseases Laboratory of the Applied Biology and Ecology Research Unit (VBID-URBEA), Department of Animal Biology, Faculty of Science of the University of Dschang, Dschang P.O. Box 067, Cameroon
5
Institute of Molecular Biology and Biotechnology, Foundation for Research and Technology-Hellas, 70013 Heraklion, Greece
6
Pesticide Science Laboratory, Department of Crop Science, Agricultural University of Athens, 11855 Athens, Greece
7
Vector Biology, Liverpool School of Tropical Medicine, Pembroke Place, Liverpool L3 5QA, UK
*
Author to whom correspondence should be addressed.
Submission received: 23 June 2022 / Revised: 19 July 2022 / Accepted: 21 July 2022 / Published: 23 July 2022
(This article belongs to the Collection Immunological Responses and Immune Defense Mechanism)

Abstract

:
The spread of pyrethroid resistance in malaria vectors is a major threat affecting the performance of current control measures. However, there is still not enough information on the resistance profile of mosquitoes to carbamates and organophosphates which could be used as alternatives. The present study assessed the resistance profile of Anopheles gambiae s.l. to bendiocarb and malathion, at the phenotypic and molecular levels, in different eco-epidemiological settings in Cameroon. Anopheles gambiae s.l. mosquitoes were collected from four eco-epidemiological settings across the country and their susceptibility level to bendiocarb and malathion was determined using WHO tubes bioassays. The ace-1 target site G119S mutation was screened by PCR. Reverse Transcription quantitative PCR 3-plex TaqMan assays were used to quantify the level of expression of eight genes associated with metabolic resistance. Resistance to malathion and/or bendiocarb was recorded in all study sites except in mosquitoes collected in Kaélé and Njombé. The Ace-1 (G119S) mutation was detected in high frequencies (>40%) in Kékem and Santchou. Both An. gambiae and An. coluzzii were detected carrying this mutation. The cytochrome P450s gene Cyp6p3 associated with carbamate resistance and the glutathione S-transferase gene Gste2 associated with organophosphate resistance were found to be overexpressed. Genes associated with pyrethroid (Cyp6m2, Cyp9k1, Cyp6p3) and organochlorine (Gste2, Cyp6z1, Cyp6m2) and cuticle resistance (Cyp4g16) were also overexpressed. The rapid spread of resistance to organophosphates and carbamates could seriously compromise future control strategies based on IRS. It is therefore becoming important to assess the magnitude of bendiocarb and malathion resistance countrywide.

1. Background

Since their discovery, synthetic insecticides have been largely used for the prevention of vector-borne diseases, due to their high efficacy and ease of use [1]. The fight against malaria vectors in Africa mainly relies on the use of insecticide-based interventions, such as indoor residual spraying (IRS) and insecticide-treated mosquito nets (ITNs).
Pyrethroids, due to their high efficacy against insects and low toxicity to mammals, including humans, are the only insecticide recommended for bed net impregnation [2]. Over 3 billion impregnated bed nets have been distributed across the world since 2010 [3]. However, overreliance on pyrethroids in both public health and agriculture has exerted intensive selective pressure on mosquitoes and has led to the emergence and expansion of pyrethroid resistance in mosquito populations [4,5,6,7]. The increasing resistance of vector species to pyrethroids is a major challenge for malaria control [8]. In several countries in Africa, kdr target site mutations, associated with pyrethroid resistance, are close to fixation and cytochrome P450 genes, encoding for pyrethroid metabolizers, are overexpressed in mosquito populations [9]. The rapid spread of pyrethroid resistance could reverse progress achieved so far in malaria control [10,11].
Carbamates and organophosphates are considered as possible alternatives for controlling pyrethroid-resistant mosquitoes [12,13]. In order to manage insecticide resistance and maintain the effectiveness of current control tools, the World Health Organization (WHO) recommends the implementation of insecticide resistance management strategies, including rotation strategies with the use of two or more insecticides with different mode of action at different periods or the combination of different control interventions [14]. Since 2006, several African countries have started combining long-lasting insecticidal nets (LLINs) and indoor residual spray (IRS) with carbamate or organophosphate for effective control of malaria vector populations [8,15]. Although promising results had been initially recorded, with improved control of vectors and malaria transmission [8], there is a growing concern since carbamate and to some extent organophosphate resistance has begun to be widely distributed in West Africa [16,17,18] and is currently also spreading in Central Africa [19,20]. LLINs are the main tools used for malaria control in Cameroon, but the country is considering complementing the use of LLINs by the addition of indoor residual spraying with organophosphate in places experiencing a high malaria burden or high seasonal transmission of the disease [8]. Carbamates and organophosphates are largely used in agriculture but have never been used in public health in Cameroon and there is, so far, not enough data on the resistance profile of mosquitoes to these insecticides. Resistance to carbamates and organophosphates is associated with the presence of the G119S mutation and the overexpression of detoxification genes [18,19,21,22]. Alarmingly, recent studies in Cameroon have indicated the emergence of resistance to carbamate and organophosphates, associated with increased frequency of the G119S mutation in An. gambiae s.l. [23,24]. Still, there is not enough information on its geographical scale and intensity.
The present study aims to characterize the resistance status of Anopheles gambiae s.l. to bendiocarb (carbamate) and malathion (organophosphate) in four epidemiological settings in Cameroon and to explore the molecular basis of resistance.

2. Results

2.1. Bioassays Results

A total of 960 Anopheles specimens from laboratory strains (Kisumu and Ngousso) and 2490 field specimens from Kékem (n = 160), Njombé (n = 200), Belabo (n = 200), Kaélé (n = 400), Tibati (n = 400), Bertoua (n = 570) and Santchou (n = 560) were exposed to determine their susceptibility profile against bendiocarb and malathion insecticide. The two laboratory strains were all susceptible to bendiocarb and malathion with mortality rate of 100%, demonstrating that the impregnated paper was good. The susceptibility profile was year- and insecticide-dependent in each locality (Figure 1). During the year 2021, in Santchou and Kekem, all tested mosquitoes were resistant to both insecticides (Figure 1). In Belabo, Bertoua and Tibati mosquitoes were found to be resistant to bendiocarb (mortality rate: 80%; 77.77%; 76% respectively) and fully susceptible to malathion (mortality rate: 98%; 99%; 100% respectively). Mosquitoes from Kaélé and Njombé were fully susceptible to both insecticides (mortality rate > 98%) (Figure 1).
In addition, the mortality rate of mosquitoes varied significantly with year, except in Santchou, where no significant difference was found regarding Bendiocab (χ2 = 3.558, p = 0.1687). Samples tested in 2020 in Bertoua and Tibati were found to be significantly more susceptible to bendiocarb (mortality rate 95% to 100%) compared to those collected in 2019 and 2021 (mortality rate varying from 31% to 78%).

2.2. Molecular Identification of Species of the Anopheles gambiae Complex

The genotyping of 637 specimens for species identification showed the presence of three members of the An. gambiae complex in the study sites, i.e., An. coluzzii, An. gambiae and An. arabiensis. In Tibati only An. coluzzii was recorded, whereas in Kaélé, An. arabiensis, An. gambiae and An. coluzzii were present with, however, a predominance of An. coluzzii. In Santchou and Kékem An. gambiae was the only species recorded, while in Bertoua, Njombé and Belabo both An. gambiae and An. coluzzii were recorded. Anopheles gambiae was the predominant species in Bertoua and Belabo, whereas An. coluzzii was slightly more abundant than An. gambiae in Njombé (Table 1).

2.3. Screening of Target Site Mutations (Ace-1 G119S)

The (G119S) mutation was detected in all sites. In Njombé the G119S mutation was detected mainly in bendiocarb and malathion survivors (75.5–100.0%) and at a very low frequency (4.5%) in the unexposed population. In Kékem bendiocarb survivors showed a frequency of 84.6% and unexposed mosquitoes a lower frequency (41.1%). The Bélabo population presented a frequency of 54.1% and 72.3% in bendiocarb and malathion survivors, respectively; it was also detected in unexposed mosquitoes in a frequency of 10.25%. In Tibati, a frequency of 81.3% was recorded in Bendiocarb survivors, whereas it was not detected in unexposed mosquitoes. In Bertoua, bendiocarb and malathion survivors showed frequencies of 71.0% and 77.1%, respectively; unexposed mosquitoes exhibited a frequency of 36.1%. In the Santchou population, the mutation was found in frequencies > 85.0% in bendiocarb and malathion survivors, as well as unexposed mosquitoes (Table 2).

2.4. Expression Analysis of Detoxification Genes

The expression profile of different detoxification genes known to confer resistance was analyzed. Cyp6p3, Cyp6m2, Cyp9k1, Cyp6p4, Cyp6z1, Gste2 and Cyp4g16 were among genes that were significantly overexpressed in at least one study population. More precisely, in the Njombé population (Figure 2A), Cyp9k1 and Cyp6z1 were significantly overexpressed (>8.0 and >1.6 folds, respectively). In Kékem the Cyp6m2 and Cyp9k1 genes were both overexpressed by >3.0 folds (Figure 2B). The Kaélé population (Figure 2C) showed a marked overexpression of Cyp6p4 (>10 folds) and Gste2 (>5.0 folds). In Bélabo, Cyp6p3 and Cyp6p4 were overexpressed by >1.90 folds (Figure 2D). The Tibati population (Figure 2E) showed an overexpression of only Cyp4g16 (>1.60 folds). In Bertoua, Cyp9k1 and Cyp6p4 were significantly overexpressed (>3.0 and >2.0 folds, respectively). The Santchou population showed a >1.6-fold overexpression of Cyp6z1 (Figure 2G).

3. Discussion

The study’s objective was to assess the distribution of carbamate and organophosphate resistance in An. gambiae s.l. populations from different epidemiological settings across Cameroon. High resistance to both carbamates and organophosphates was detected and was consistent with previous reports [23,24,25]. Resistance to bendiocarb was largely spread and was detected in five out of the seven sites surveyed, whereas resistance to malathion was detected in the sites of Kékem and Santchou situated in western Cameroon. These sites are semi-rural settings where vast lands are used for the cultivation of vegetables, maize, cocoa, plantains and tomatoes all year round with the intensive application of pesticides [26]. It is possible that high resistance to bendiocarb and malathion detected in these sites is driven by the frequent use of pesticides in agriculture rather than insecticides used in public health. In Tibati, Bertoua and Belabo, mosquitoes were found resistant to bendiocarb, while mosquitoes from Njombé and Kaélé were fully susceptible to bendiocarb and malathion. The heterogeneous pattern of carbamates and organophosphates resistance in An. gambiae s.l. populations across the country could derive from different selective pressure in each site and habitat segregation by the species of the An. gambiae complex [27,28]. Similar geographical patterns of carbamate or organophosphate resistance in members of the An. gambiae complex have been reported in other African countries [22]. Anopheles arabiensis, An. gambiae and An. coluzzii species were recorded in the study sites and their presence followed an eco-epidemiological/site-dependent pattern. More precisely, An. coluzzii was the predominant species in Kaélé, Tibati and Njombé, whereas An. gambiae was predominant in Santchou, Bertoua, Belabo and Kékem. An. arabiensis was detected in low frequency in Kaelé and Tibati. The high frequency of An. coluzzii in the Sahelian and Sudano-Sahelian region of the country (Kaélé and Tibati) is in line with previous studies which indicated that this species is outpacing An. arabiensis and is gradually becoming the predominant species in the Sudano-Sahelian savannah zone [29,30].
In an effort to molecularly profile the evident carbamate and organophosphate phenotypic resistance we measured the frequencies of the Ace-1 G119S mutation, as well as the expression levels of genes associated with metabolic and cuticle resistance.
The Ace-1 G119S mutation was detected at a very high frequency in almost all sites. Studies conducted so far in Africa suggested that the presence of this mutation is mainly associated with resistance to carbamates [17,18,22,31]. The Ace-1 G119S mutation was detected in both An. gambiae and An. coluzzii with, however, a higher frequency in An. gambiae. A similar distribution pattern has been reported in many countries across West Africa [17,22]. The present findings are different from previous studies [23,25], which reported the mutation only in An. gambiae. It is possible that the Ace-1 G119S mutation in An. coluzzii may have emerged during recent years through introgression or independently. Different frequencies of the G119S allele have been reported between An. gambiae and An. coluzzii in areas of sympatry. It has been suggested that the difference in this mutation frequency could result from discrepancies in the selection pressure mosquitoes experienced in aquatic habitats [22]. Despite the fact that the study did not characterize aquatic habitats, studies assessing breeding habitats influence in the cities of Yaoundé and Douala in Cameroon suggested that mosquitoes deriving from agricultural cultivated areas were more tolerant to different insecticides families, including pyrethroids and carbamates, compared to those originating from unpolluted sites [19,32].
In addition to the Ace-1 mutation the expression profiles of genes implicated in metabolic and cuticle resistance were assessed. Regarding carbamate resistance, Cyp6p3 was significantly overexpressed in the Bélabo population, where phenotypic resistance to bendiocarb was also detected. This particular P450 gene has been previously shown to metabolize bendiocarb [33,34]. Regarding organophosphate resistance, in Santchou, a site with low mortality rates for malathion, the overexpression of Cyp6p4 was detected. Cyp6p4 has been previously shown to metabolize malathion, as well as other organophosphates (pyrimiphos-mehtyl, fenitrothion); however, a link to resistance cannot be safely concluded, and organophosphate metabolism by Cyp6p4 could lead either to the activation of more toxic metabolites or detoxification, and the relative activation versus detoxification contribution is not clear. This is also the case for Cyp6p3 and Cyp6m2 with data for the latter pointing towards activation rather than detoxification. The glutathione S-transferase gene gste2 was found to be markedly overexpressed in the Kaélé population. Since, gste2 has been strongly linked with resistance to the organophosphate fenitrothion [34] this constitutes an alarming factor. Our expression analysis revealed also the significant overexpression of genes that are known pyrethroid (Cyp6m2, Cyp9k1, Cyp6p4) and organochlorine (gste2, Cyp6z1, Cyp6m2) metabolizers [9]. Resistance to these insecticides is well documented in Cameroon [35]. All the above show that cross-resistance and multi-resistance is in place in the study’s populations. Regarding cuticle resistance, Cyp4g16, a functional oxidative decarboxylase gene known to catalyze epicuticular lipid biosynthesis and contribute to insecticide resistance via the enrichment of the CHC content, thus reducing insecticide uptake was also found to be upregulated in one population (Tibati) [36].

4. Methods

4.1. Study Site

Mosquitoes were collected from seven localities in Cameroon: Kékem, Njombé, Belabo, Kaélé, Tibati, Bertoua and Santchou (Figure 3). A detailed characteristic of each site is provided in Table 3. Physical characteristics of sites with superscript (a) is published in Kala-Chouakeu et al. [37].

4.2. Collection of Mosquito Larvae, Rearing, and Conservation

The field sampling of anopheline larvae and processing was conducted during the long rainy season in all the sites. It was carried out once yearly in August 2019, August 2020 and November 2021 in Bertoua and Santchou, whereas in the remaining sites the sampling and processing of mosquitoes was performed once or twice (Kékem, Njombé, Belabo, Kaélé, Tibati) because of the unavailability of funds (Figure 3). Larval collections were undertaken in different habitats including temporary water collections, puddles, and semi-permanent sites. In each study site, collected samples from different breeding habitats were pooled per site and reared into adults in the insectary created in the field. Larvae were fed with TetraMin® fish food until pupae. Pupae were collected in a cardboard cup and placed in netting cages for adult emergence. After emergence, adults were offered sugar solution until processing. A subset of 30–40 unexposed, non-blood fed, 3–5-day-old female An. gambiae s.l. from different sites and susceptible strains (the An. gambiae Kisumu strain and the An. coluzzii Ngousso strain) were preserved in RNAlater™ for the characterization of the molecular mechanisms of insecticide resistance. The same sample was also used for species identification. Another subset of about 120–150 mosquitoes per population was used for insecticide bioassays; survivors after exposure to insecticide were preserved in 70% alcohol and used for the confirmation of molecular species identification and detection of G119S.

4.3. Insecticide Bioassay

Adult female An. gambiae s.l. reared from larval collections in different sites were tested alongside the susceptible laboratory strains Kisumu and Ngousso against two insecticides (bendiocarb 0.01% and malathion 5%) following WHO guidelines [38]. An. gambiae s.l. females aged 3–5 days reared from larvae collected on the field were placed in batches of 20 to 25 mosquitoes per tube and left for observation for one hour. After this period, mosquitoes were transferred to four different tubes with insecticide-impregnated papers and exposed for 1 h. The susceptible laboratory strains (An. gambiae (Kisumu strain) and An. coluzzii (Ngousso strain)) were used as control to assess the quality of the impregnated papers. The number of mosquitoes knocked down by the insecticide was recorded after 1h of exposure; then, mosquitoes were fed with a 10% sugar solution and the number of dead mosquitoes was recorded 24 h post-exposure. Mosquitoes were considered resistant when the mortality rate was <90% and susceptible when the mortality rate was ≥98%, and resistance status needed further scrutiny when the mortality rate was <98% and >90% [39].

4.4. Mosquito Processing

4.4.1. Total RNA and DNA Extraction from Mosquito Pools

Using a magnetic beads-based approach with the MagSi kit (MagnaMedics Diagnostics B.V., Geleen, Netherlands), total genomic material (RNA and DNA) was extracted from 687 mosquitoes (N = 90 from Kékem, N = 64 from Njombé, N = 94 from Belabo, N = 70 from Kaélé, N = 90 from Tibati, N = 130 from Santchou and N = 89 from Bertoua, N = 30 from Kisumu and N = 30 from Ngousso, susceptible strains). Mosquitoes were pooled for extractions (up to N = 10 mosquitoes per pool). The quantity and purity of DNA and total RNA were assessed spectrophotometrically via Nanodrop measurements. The average NA concentration was 58.96 ± 5.4 ng/μL and the average A260/A280 ratio, indicative of sample purity, was 2.3 ± 0.05. The quality of RNA was assessed by 1.0% w/v agarose gel electrophoresis. Representative samples are shown in Supplementary Figure S1.

4.4.2. Species Identification

Mosquitoes used for species identification included those exposed to insecticide (survivors, dead) and some not exposed to insecticides. Mosquitoes from different years were analyzed. Species identification at the molecular level was performed using the TaqMan assays described in the Innovative Vector Control Consortium (IVCC) Vector Population Monitoring Tool (VPMT) [39] with modifications described in Wipf et al. [40] for the detection of An. coluzzii and An. gambiae species (Primers: S200_X6.1_F: TCGCCTTAGACCTTGCGTTA, S200_X6.1_R: CGCTTCAAGAATTCGAGATAC, Probes:Pcoluzzii: HEX-ACCGCGCCGCCATACGTAGGA-BHQ1 and Pgambiae: FAM-ATGTCTAATAGTCTCAATAGT -MGB).

4.5. Genotyping of G119S Mutation

The G119S mutation was analyzed using the TaqMan assay (Primers: F: GGCCGTCATGCTGTGGAT and R: GCGGTGCCGGAGTAGA; Probes: Pwt: HEX-TTCGGCGGCGGCT-MGB and Pmut: FAM-TTCGGCGGCAGCT-MGB), as described in VPMT [39]. The percentage of allele frequency for the previously mentioned traits in mosquito pools was calculated with regression models using the protocol described in Mavridis et al., 2018 [41]. Briefly, a standard curve for G119S using plasmid sequences with known % frequencies was constructed. It was then used for the calculation of % G119S frequencies for the unknown populations by transforming the dCt (CtFAM-CtHEX values) to allelic frequencies via the equation of the standard curve’s regression model (Supplementary Figure S2).

4.6. Resistance Gene Expression Analyses

The reverse transcription quantitative PCR (qRT-PCR) 3-plex TaqMan® assays described in Mavridis et al., 2019 [42] were used for the quantification of seven detoxification genes (Cyp6p3, Cyp6m2, Cyp9k1, Cyp6p4, Cyp6z1, Gste2, Cyp6p1) that have been strongly associated with metabolic resistance and one oxidative decarboxylase (Cyp4g16) that is implicated in cuticle resistance [42]. The list of primers and probes that were used are listed in Supplementary Table S1. Reactions were performed in the Viia7 Real-Time PCR system (Applied Biosystems) using the following thermal cycle parameters: 50 °C for 15 min, 95 °C for 3 min and 40 cycles of 95 °C for 3 s and 60 °C for 30 s. The QuantStudioTM Real-Time PCR system v1.3 (Applied Biosystems) software was used for the calculation of Ct values for each reaction.

4.7. Statistical Analysis

Using the method developed by Pfaff [43], fold-changes, 95% CIs and statistical significance were calculated, while graphs were constructed with the SigmaPlot software (v12.0). The mortality rate was expressed as the ratio between the number of mosquitoes that were found dead or not capable of standing on their legs and exposed ones. Confidence intervals were computed using Medcalc. The comparison of mortality rate and fold change was performed using chi square and Student’s t-test, respectively.

5. Conclusions

This study indicates the heterogeneous distribution of carbamate and organophosphate resistance in An. gambiae s.l. populations across Cameroon. Both Ace-1 G119S mutation and the overexpression of genes associated with carbamate/organophosphate resistance were detected. The current expansion of carbamates and organophosphate resistance could seriously compromise future malaria control measures based on the use of organophosphates and carbamates in IRS. Addressing gaps affecting malaria control by improving surveillance activities on the vector is becoming crucial to ensure the success of malaria elimination efforts in Cameroon.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/pathogens11080824/s1; Figure S1: Agarose gel (1.0% w/v) electrophoresis for total RNA of representative samples (N = 6). The presence of distinct ribosomal bands (28S, 18S, 5S from top to bottom) and the absence of degradation products shows that total RNA was intact, suitable for downstream analyses.; Figure S2: Standard curve for G119S using plasmid sequences with known % frequencies (A) and calculation of % G119S for an unknown population of the study as an example (B); Table S1: List of primers and probes used for the multiplex resistance gene expression analyses.

Author Contributions

I.N.N.-S., P.N.-M. and N.A.K.-C.: participated in data collection; performed laboratory analysis, analyzed and interpreted the data; wrote the manuscript; revised the manuscript. R.B., A.T. and L.D.-D.: participated in data collection, laboratory analysis and revised the manuscript. J.V. and K.M.: participated in laboratory analysis and revised the manuscript. J.T., T.T. and J.A.M.-M.: revised the manuscript. C.A.-N.: Conceived and designed the study protocol; critically revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work received financial support from the Bill & Melinda Gates Foundation, Grant ID: OPP1210340 to CAN, the European Union’s Horizon 2020 Research and Innovation Programme under grant agreement No. 731060 (INFRAVEC2) and the Panafrican Mosquito Association (PAMCA). The funding body did not have any role in the experimental design, collection of data, analysis, and interpretation of data or in the writing of the manuscript.

Institutional Review Board Statement

The study was conducted under the ethical clearance No. 2020/04/1209/CE/CNERSH/SP delivered by the Cameroon National Ethics (CNE) Committee for Research on Human Health.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data have been made available in the manuscript and supplementary files.

Acknowledgments

We are grateful to the administrative and traditional authorities, mosquito collectors and the population of Kaélé, Tibati, Santchou, Bertoua, Njombé, Kékem and Belabo for their participation in the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. van den Berg, H.; da Silva Bezerra, H.S.; Al-Eryani, S.; Chanda, E.; Nagpal, B.N.; Knox, T.B.; Velayudhan, R.; Yadav, R.S. Recent Trends in Global Insecticide Use for Disease Vector Control and Potential Implications for Resistance Management. Sci. Rep. 2021, 11, 23867. [Google Scholar] [CrossRef] [PubMed]
  2. Zaim, M.; Aitio, A.; Nakashima, N. Safety of Pyrethroid-Treated Mosquito Nets. Med. Vet. Entomol. 2000, 14, 1–5. [Google Scholar] [CrossRef]
  3. World Malaria Report 2020. Available online: https://www.who.int/publications-detail-redirect/9789240015791 (accessed on 1 June 2022).
  4. Antonio-Nkondjio, C.; Sonhafouo-Chiana, N.; Ngadjeu, C.S.; Doumbe-Belisse, P.; Talipouo, A.; Djamouko-Djonkam, L.; Kopya, E.; Bamou, R.; Awono-Ambene, P.; Wondji, C.S. Review of the Evolution of Insecticide Resistance in Main Malaria Vectors in Cameroon from 1990 to 2017. Parasit. Vectors 2017, 10, 472. [Google Scholar] [CrossRef] [PubMed]
  5. Djiappi-Tchamen, B.; Nana-Ndjangwo, M.S.; Mavridis, K.; Talipouo, A.; Nchoutpouen, E.; Makoudjou, I.; Bamou, R.; Mayi, A.M.P.; Awono-Ambene, P.; Tchuinkam, T.; et al. Analyses of Insecticide Resistance Genes in Aedes Aegypti and Aedes Albopictus Mosquito Populations from Cameroon. Genes 2021, 12, 828. [Google Scholar] [CrossRef] [PubMed]
  6. Ranson, H.; Lissenden, N. Insecticide Resistance in African Anopheles Mosquitoes: A Worsening Situation That Needs Urgent Action to Maintain Malaria Control. Trends Parasitol. 2016, 32, 187–196. [Google Scholar] [CrossRef]
  7. Talipouo, A.; Mavridis, K.; Nchoutpouen, E.; Djiappi-Tchamen, B.; Fotakis, E.A.; Kopya, E.; Bamou, R.; Kekeunou, S.; Awono-Ambene, P.; Balabanidou, V.; et al. High Insecticide Resistance Mediated by Different Mechanisms in Culex Quinquefasciatus Populations from the City of Yaoundé, Cameroon. Sci. Rep. 2021, 11, 7322. [Google Scholar] [CrossRef]
  8. Hemingway, J.; Vontas, J.; Poupardin, R.; Raman, J.; Lines, J.; Schwabe, C.; Matias, A.; Kleinschmidt, I. Country-Level Operational Implementation of the Global Plan for Insecticide Resistance Management. Proc. Natl. Acad. Sci. USA 2013, 110, 9397–9402. [Google Scholar] [CrossRef] [Green Version]
  9. Vontas, J.; Katsavou, E.; Mavridis, K. Cytochrome P450-Based Metabolic Insecticide Resistance in Anopheles and Aedes Mosquito Vectors: Muddying the Waters. Pestic. Biochem. Physiol. 2020, 170, 104666. [Google Scholar] [CrossRef] [PubMed]
  10. Bhatt, S.; Weiss, D.J.; Cameron, E.; Bisanzio, D.; Mappin, B.; Dalrymple, U.; Battle, K.; Moyes, C.L.; Henry, A.; Eckhoff, P.A.; et al. The Effect of Malaria Control on Plasmodium Falciparum in Africa between 2000 and 2015. Nature 2015, 526, 207–211. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Churcher, T.S.; Lissenden, N.; Griffin, J.T.; Worrall, E.; Ranson, H. The Impact of Pyrethroid Resistance on the Efficacy and Effectiveness of Bednets for Malaria Control in Africa. eLife 2016, 5, e16090. [Google Scholar] [CrossRef]
  12. Agossa, F.R.; Aïkpon, R.; Azondékon, R.; Govoetchan, R.; Padonou, G.G.; Oussou, O.; Oké-Agbo, F.; Akogbéto, M.C. Efficacy of Various Insecticides Recommended for Indoor Residual Spraying: Pirimiphos Methyl, Potential Alternative to Bendiocarb for Pyrethroid Resistance Management in Benin, West Africa. Trans. R. Soc. Trop. Med. Hyg. 2014, 108, 84–91. [Google Scholar] [CrossRef]
  13. Blaschke, T.F.; Lumpkin, M.; Hartman, D. The World Health Organization Prequalification Program and Clinical Pharmacology in 2030. Clin. Pharmacol. Ther. 2020, 107, 68–71. [Google Scholar] [CrossRef]
  14. World Health Organization. World Health Statistics 2012; World Health Organization: Geneva, Switzerland, 2012. [Google Scholar]
  15. Dengela, D.; Seyoum, A.; Lucas, B.; Johns, B.; George, K.; Belemvire, A.; Caranci, A.; Norris, L.C.; Fornadel, C.M. Multi-Country Assessment of Residual Bio-Efficacy of Insecticides Used for Indoor Residual Spraying in Malaria Control on Different Surface Types: Results from Program Monitoring in 17 PMI/USAID-Supported IRS Countries. Parasit. Vectors 2018, 11, 71. [Google Scholar] [CrossRef]
  16. Aïkpon, R.; Agossa, F.; 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Djogbénou, L.; Weill, M.; Hougard, J.-M.; Raymond, M.; Akogbéto, M.; Chandre, F. Characterization of Insensitive Acetylcholinesterase (Ace-1R) in Anopheles Gambiae (Diptera: Culicidae): Resistance Levels and Dominance. J. Med. Entomol. 2007, 44, 805–810. [Google Scholar] [CrossRef] [PubMed]
  18. Edi, C.V.A.; Koudou, B.G.; Jones, C.M.; Weetman, D.; Ranson, H. Multiple-Insecticide Resistance in Anopheles Gambiae Mosquitoes, Southern Côte d’Ivoire. Emerg. Infect. Dis. 2012, 18, 1508–1511. [Google Scholar] [CrossRef] [PubMed]
  19. Antonio-Nkondjio, C.; Poupardin, R.; Tene, B.F.; Kopya, E.; Costantini, C.; Awono-Ambene, P.; Wondji, C.S. Investigation of Mechanisms of Bendiocarb Resistance in Anopheles Gambiae Populations from the City of Yaoundé, Cameroon. Malar. J. 2016, 15, 424. [Google Scholar] [CrossRef] [Green Version]
  20. Menze, B.D.; Wondji, M.J.; Tchapga, W.; Tchoupo, M.; Riveron, J.M.; Wondji, C.S. Bionomics and Insecticides Resistance Profiling of Malaria Vectors at a Selected Site for Experimental Hut Trials in Central Cameroon. Malar. J. 2018, 17, 317. [Google Scholar] [CrossRef]
  21. Aïkpon, R.; Sèzonlin, M.; Ossè, R.; Akogbéto, M. Evidence of Multiple Mechanisms Providing Carbamate and Organophosphate Resistance in Field An. Gambiae Population from Atacora in Benin. Parasit. Vectors 2014, 7, 568. [Google Scholar] [CrossRef]
  22. Essandoh, J.; Yawson, A.E.; Weetman, D. Acetylcholinesterase (Ace-1) Target Site Mutation 119S Is Strongly Diagnostic of Carbamate and Organophosphate Resistance in Anopheles Gambiae s.s. and Anopheles Coluzzii across Southern Ghana. Malar. J. 2013, 12, 404. [Google Scholar] [CrossRef] [Green Version]
  23. Bamou, R.; Sonhafouo-Chiana, N.; Mavridis, K.; Tchuinkam, T.; Wondji, C.S.; Vontas, J.; Antonio-Nkondjio, C. Status of Insecticide Resistance and Its Mechanisms in Anopheles Gambiae and Anopheles Coluzzii Populations from Forest Settings in South Cameroon. Genes 2019, 10, 741. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Binyang, A.J.; Elanga-Ndille, E.; Tene-Fossog, B.; Ndo, C.; Nouage, L.; Assatse, T.; Fotso-Toguem, Y.; Tabue, R.; Zeukeng, F.; Nguiffo, D.N.; et al. Distribution of Acetylcholinesterase (Ace-1R) Target-Site G119S Mutation and Resistance to Carbamates and Organophosphates in Anopheles Gambiae Sensu Lato Populations from Cameroon. Parasit. Vectors 2022, 15, 53. [Google Scholar] [CrossRef]
  25. Elanga-Ndille, E.; Nouage, L.; Ndo, C.; Binyang, A.; Assatse, T.; Nguiffo-Nguete, D.; Djonabaye, D.; Irwing, H.; Tene-Fossog, B.; Wondji, C.S. The G119S Acetylcholinesterase (Ace-1) Target Site Mutation Confers Carbamate Resistance in the Major Malaria Vector Anopheles Gambiae from Cameroon: A Challenge for the Coming IRS Implementation. Genes 2019, 10, 790. [Google Scholar] [CrossRef] [Green Version]
  26. Pouokam, G.B.; Lemnyuy Album, W.; Ndikontar, A.S.; Sidatt, M.E.H. A Pilot Study in Cameroon to Understand Safe Uses of Pesticides in Agriculture, Risk Factors for Farmers’ Exposure and Management of Accidental Cases. Toxics 2017, 5, 30. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Kamdem, C.; Fossog, B.T.; Simard, F.; Etouna, J.; Ndo, C.; Kengne, P.; Boussès, P.; Etoa, F.-X.; Awono-Ambene, P.; Fontenille, D.; et al. Anthropogenic Habitat Disturbance and Ecological Divergence between Incipient Species of the Malaria Mosquito Anopheles Gambiae. PLoS ONE 2012, 7, e39453. [Google Scholar] [CrossRef] [Green Version]
  28. Tene Fossog, B.; Ayala, D.; Acevedo, P.; Kengne, P.; Ngomo Abeso Mebuy, I.; Makanga, B.; Magnus, J.; Awono-Ambene, P.; Njiokou, F.; Pombi, M.; et al. Habitat Segregation and Ecological Character Displacement in Cryptic African Malaria Mosquitoes. Evol. Appl. 2015, 8, 326–345. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Ibrahim, S.S.; Fadel, A.N.; Tchouakui, M.; Terence, E.; Wondji, M.J.; Tchoupo, M.; Kérah-Hinzoumbé, C.; Wanji, S.; Wondji, C.S. High Insecticide Resistance in the Major Malaria Vector Anopheles Coluzzii in Chad Republic. Infect. Dis. Poverty 2019, 8, 100. [Google Scholar] [CrossRef]
  30. Fadel, A.N.; Ibrahim, S.S.; Tchouakui, M.; Terence, E.; Wondji, M.J.; Tchoupo, M.; Wanji, S.; Wondji, C.S. A Combination of Metabolic Resistance and High Frequency of the 1014F Kdr Mutation Is Driving Pyrethroid Resistance in Anopheles Coluzzii Population from Guinea Savanna of Cameroon. Parasit. Vectors 2019, 12, 263. [Google Scholar] [CrossRef]
  31. Asidi, A.N.; N’ Guessan, R.; Koffi, A.A.; Curtis, C.F.; Hougard, J.-M.; Chandre, F.; Corbel, V.; Darriet, F.; Zaim, M.; Rowland, M.W. Experimental Hut Evaluation of Bednets Treated with an Organophosphate (Chlorpyrifos-Methyl) or a Pyrethroid (Lambdacyhalothrin) Alone and in Combination against Insecticide-Resistant Anopheles Gambiae and Culex Quinquefasciatus Mosquitoes. Malar. J. 2005, 4, 25. [Google Scholar] [CrossRef] [Green Version]
  32. Antonio-Nkondjio, C.; Fossog, B.T.; Ndo, C.; Djantio, B.M.; Togouet, S.Z.; Awono-Ambene, P.; Costantini, C.; Wondji, C.S.; 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. [Google Scholar] [CrossRef] [Green Version]
  33. Yunta, C.; Hemmings, K.; Stevenson, B.; Koekemoer, L.L.; Matambo, T.; Pignatelli, P.; Voice, M.; Nász, S.; Paine, M.J.I. Cross-Resistance Profiles of Malaria Mosquito P450s Associated with Pyrethroid Resistance against WHO Insecticides. Pestic. Biochem. Physiol. 2019, 161, 61–67. [Google Scholar] [CrossRef] [PubMed]
  34. Adolfi, A.; Poulton, B.; Anthousi, A.; Macilwee, S.; Ranson, H.; Lycett, G.J. Functional Genetic Validation of Key Genes Conferring Insecticide Resistance in the Major African Malaria Vector, Anopheles Gambiae. Proc. Natl. Acad. Sci. USA 2019, 116, 25764–25772. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Nwane, P.; Etang, J.; Chouaibou, M.; Toto, J.C.; Kerah-Hinzoumbé, C.; Mimpfoundi, R.; Awono-Ambene, H.P.; Simard, F. Trends in DDT and Pyrethroid Resistance in Anopheles Gambiaes.s. Populations from Urban and Agro-Industrial Settings in Southern Cameroon. BMC Infect. Dis. 2009, 9, 163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Balabanidou, V.; Kampouraki, A.; MacLean, M.; Blomquist, G.J.; Tittiger, C.; Juárez, M.P.; Mijailovsky, S.J.; Chalepakis, G.; Anthousi, A.; Lynd, A.; et al. Cytochrome P450 Associated with Insecticide Resistance Catalyzes Cuticular Hydrocarbon Production in Anopheles Gambiae. Proc. Natl. Acad. Sci. USA 2016, 113, 9268–9273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kala Chouakeu, N.A.; Ngingahi, L.G.; Bamou, R.; Talipouo, A.; Ngadjeu, C.S.; Mayi, M.P.A.; Kopya, E.; Awono-Ambene, P.; Tchuinkam, T.; Antonio Nkondjio, C. Knowledge, Attitude, and Practices (KAP) of Human Populations towards Malaria Control in Four Ecoepidemiological Settings in Cameroon. J. Trop. Med. 2021, 2021, e9925135. [Google Scholar] [CrossRef] [PubMed]
  38. World Health Organization. Test Procedures for Insecticide Resistance Monitoring in Malaria Vector Mosquitoes, 2nd ed.; World Health Organization: Geneva, Switzerland, 2018; ISSN 978-92-4-151157-5. [Google Scholar]
  39. Bass, C.; Nikou, D.; Vontas, J.; Williamson, M.S.; Field, L.M. Development of High-Throughput Real-Time PCR Assays for the Identification of Insensitive Acetylcholinesterase (Ace-1R) in Anopheles Gambiae. Pestic. Biochem. Physiol. 2010, 96, 80–85. [Google Scholar] [CrossRef]
  40. Wipf, N.C.; Duchemin, W.; Kouadio, F.-P.A.; Fodjo, B.K.; Sadia, C.G.; Mouhamadou, C.S.; Vavassori, L.; Mäser, P.; Mavridis, K.; Vontas, J.; et al. Multi-Insecticide Resistant Malaria Vectors in the Field Remain Susceptible to Malathion, despite the Presence of Ace1 Point Mutations. PLoS Genet. 2022, 18, e1009963. [Google Scholar] [CrossRef]
  41. Mavridis, K.; Wipf, N.; Müller, P.; Traoré, M.M.; Muller, G.; Vontas, J. Detection and Monitoring of Insecticide Resistance Mutations in Anopheles Gambiae: Individual vs. Pooled Specimens. Genes 2018, 9, 479. [Google Scholar] [CrossRef] [Green Version]
  42. Mavridis, K.; Wipf, N.; Medves, S.; Erquiaga, I.; Müller, P.; Vontas, J. Rapid Multiplex Gene Expression Assays for Monitoring Metabolic Resistance in the Major Malaria Vector Anopheles Gambiae. Parasit. Vectors 2019, 12, 9. [Google Scholar] [CrossRef] [Green Version]
  43. Pfaffl, M.W. A New Mathematical Model for Relative Quantification in Real-Time RT–PCR. Nucleic Acids Res. 2001, 29, e45. [Google Scholar] [CrossRef]
Figure 1. Mortality rate (%) of field collected and laboratory strains of Anopheles gambiae s.l. to bendiocarb 0.01% (carbamate) (A) and malathion 5% (organophosphate) (B) over time. At least 20–25 mosquitoes in four replicates were exposed per insecticides per site and during each collection period. Color represents the collection period (blue: August 2019; orange: August 2020 and grey: November 2021). The line (90%) represents the threshold of susceptibility (resistant when the mortality rate was <90% and susceptible when the mortality rate was >90%. The upper script (*) p < 0.05, (**) p < 0.01; (NS) non-significant. The bars represent standard error.
Figure 1. Mortality rate (%) of field collected and laboratory strains of Anopheles gambiae s.l. to bendiocarb 0.01% (carbamate) (A) and malathion 5% (organophosphate) (B) over time. At least 20–25 mosquitoes in four replicates were exposed per insecticides per site and during each collection period. Color represents the collection period (blue: August 2019; orange: August 2020 and grey: November 2021). The line (90%) represents the threshold of susceptibility (resistant when the mortality rate was <90% and susceptible when the mortality rate was >90%. The upper script (*) p < 0.05, (**) p < 0.01; (NS) non-significant. The bars represent standard error.
Pathogens 11 00824 g001
Figure 2. Expression analysis of genes associated with insecticide resistance in the study’s populations (parts AG). Error bars indicate 95% CIs. * denotes genes that showed statistically significant upregulation. ▲ denotes genes that showed consistent upregulation compared to both susceptible strains.
Figure 2. Expression analysis of genes associated with insecticide resistance in the study’s populations (parts AG). Error bars indicate 95% CIs. * denotes genes that showed statistically significant upregulation. ▲ denotes genes that showed consistent upregulation compared to both susceptible strains.
Pathogens 11 00824 g002
Figure 3. Map of Cameroon showing the study sites. Number of collection periods (* = once, ** = twice, *** = thrice).
Figure 3. Map of Cameroon showing the study sites. Number of collection periods (* = once, ** = twice, *** = thrice).
Pathogens 11 00824 g003
Table 1. Distribution of An. gambiae s.l. species at the different study sites.
Table 1. Distribution of An. gambiae s.l. species at the different study sites.
Population/SitesSites CharacteristicsSample SizeSpecies ID% Frequency
Kisumu Susceptible strain40An. gambiae100
Ngousso 40An. coluzzii100
TibatiSahelo-Sudanese (humid savannah)80An. coluzzii100
KaéléSahelian zone An. arabiensis10.0
60An. gambiae4.0
An. coluzzii86.0
KékemHighland grassfields80An. gambiae100
Santchou120An. gambiae100
BertouaForest zone79An. gambiae87.0
An. coluzzii13.0
Njombé54An. gambiae40.9
An. coluzzii59.1
Belabo84An. gambiae82.0
An. coluzzii18.0
Table 2. Ace-1 G119S mutation allelic frequencies in the different samples (mean ± SE).
Table 2. Ace-1 G119S mutation allelic frequencies in the different samples (mean ± SE).
PopulatioNCategorySample Size% G119S (ace-1)
KisumuSusceptible lab strain400.0 ± 0.0
NgoussoSusceptible lab strain400.0 ± 0.0
NjombéBendiocarb survivors375.5
Malathion survivors1100.0
Unexposed204.5 ± 2.5
KékemBendiocarb survivors2084.6 ± 5.1
Unexposed2041.1 ± 18.3
BélaboBendiocarb survivors2054.1 ± 9.6
Malathion survivors472.3
Unexposed2010.25 ± 5.9
TibatiBendiocarb survivors2081.3 ± 18.7
Unexposed200.0 ± 0.0
BertouaBendiocarb survivors1871.0 ± 9.9
Malathion survivors177.1
Unexposed2036.1 ± 9.4
SantchouBendiocarb survivors2087.3 ± 1.75
Malathion survivors4086.93 ± 8.0
Unexposed2087.9 ± 12.2
Table 3. Physical characteristics of the study sites.
Table 3. Physical characteristics of the study sites.
Study sitesKékemNjombéBelaboKaélé aTibati aSantchou aBertoua a
Administrative regionWestLittoralEastFar NorthAdamawaWestEast
Coordinates5°10′ N,
10°02′ E
4°64′ N,
9°67′ E
4°56′ N,
13°18′ E
10°50′ N,
14°56′ E
12°37′ N,
12°37′ E
5°58′ N,
9°58′ E
4°34′ N,
13°41′ E
DomainHighland GrassfieldsForest zoneForest zoneSahelian zoneSahelo soudanese Highland GrassfieldsForest zone
ClimateEquatorialEquatorialSubtropicalSahelianTropical humidEquatorialSubtropical
SeasonsDry season (November to March), rainy season (April to October)Dry season (December to February), rainy season (March to November) Dry season (December to March and July), rainy season (March to June and August to November)Dry season (October to May), rainy season (June to September)Dry season (November to April), rainy season (May to October) Dry season (November to March), rainy season (April to October), Dry season (December to March and July), rainy season (March to June and August to November)
VegetationGrasslandGrasslandEvergreen degraded forest Wooded SavanahGrassy SavanahGrasslandEvergreen degraded forest
a Physical characteristics of sites published in Kala-Chouakeu et al. [37].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ngangue-Siewe, I.N.; Ndjeunia-Mbiakop, P.; Kala-Chouakeu, N.A.; Bamou, R.; Talipouo, A.; Djamouko-Djonkam, L.; Vontas, J.; Mavridis, K.; Tombi, J.; Tchuinkam, T.; et al. Bendiocarb and Malathion Resistance in Two Major Malaria Vector Populations in Cameroon Is Associated with High Frequency of the G119S Mutation (Ace-1) and Overexpression of Detoxification Genes. Pathogens 2022, 11, 824. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens11080824

AMA Style

Ngangue-Siewe IN, Ndjeunia-Mbiakop P, Kala-Chouakeu NA, Bamou R, Talipouo A, Djamouko-Djonkam L, Vontas J, Mavridis K, Tombi J, Tchuinkam T, et al. Bendiocarb and Malathion Resistance in Two Major Malaria Vector Populations in Cameroon Is Associated with High Frequency of the G119S Mutation (Ace-1) and Overexpression of Detoxification Genes. Pathogens. 2022; 11(8):824. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens11080824

Chicago/Turabian Style

Ngangue-Siewe, Idriss Nasser, Paulette Ndjeunia-Mbiakop, Nelly Armanda Kala-Chouakeu, Roland Bamou, Abdou Talipouo, Landre Djamouko-Djonkam, John Vontas, Konstantinos Mavridis, Jeannette Tombi, Timoléon Tchuinkam, and et al. 2022. "Bendiocarb and Malathion Resistance in Two Major Malaria Vector Populations in Cameroon Is Associated with High Frequency of the G119S Mutation (Ace-1) and Overexpression of Detoxification Genes" Pathogens 11, no. 8: 824. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens11080824

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop