Molecular diversity of genes related to biological rhythms (period and timeless) and insecticide resistance (Na _V and ace-1) in Anopheles darlingi

Abstract

BACKGROUND

Malaria is a public health concern in the Amazonian Region, where Anopheles darlingi is the main vector of Plasmodium spp. Several studies hypothesised the existence of cryptic species in An. darlingi, considering variations in behaviour, morphological and genetic aspects. Determining their overall genetic background for vector competence, insecticide resistance, and other elements is essential to better guide strategies for malaria control.

OBJECTIVES

This study aimed to evaluate the molecular diversity in genes related to behaviour and insecticide resistance, estimating genetic differentiation in An. darlingi populations from Amazonian localities in Brazil and Pacific Colombian region.

METHODS

We amplified, cloned and sequenced fragments of genes related to behaviour: timeless (tim) and period (per), and to insecticide resistance: voltage-gated sodium channel (Na _V ) and acetylcholinesterase (ace-1) from 516 An. darlingi DNA samples from Manaus, Unini River, Jaú River and Porto Velho - Brazil, and Chocó - Colombia. We discriminated single nucleotide polymorphisms (SNPs), determined haplotypes and evaluate the phylogenetic relationship among the populations.

FINDINGS

The genes per, tim and ace-1 were more polymorphic than Na _V . The classical kdr and ace-1 ^R mutations were not observed. Phylogenetic analyses suggested a significant differentiation between An. darlingi populations from Brazil and Colombia, except for the Na _V gene. There was a geographic differentiation within Brazilian populations considering per and ace-1.

CONCLUSIONS

Our results add genetic data to the discussion about polymorphisms at population levels in An. darlingi. The search for insecticide resistance-related mechanisms should be extended to more populations, especially from localities with a vector control failure scenario.

Key words:
malaria vector; behavioural genes; insecticide resistance; population genetics

Malaria is one of the deadliest tropical diseases. In 2021, 247 million malaria cases occurred worldwide, and in the same year, Brazil recorded 163,585 cases.¹ In this country, malaria prevails in the Amazonian Region, including the states of Acre, Amapá, Amazonas, Maranhão, Mato Grosso, Pará, Rondônia, Roraima and Tocantins.² In Colombia, in 2020, 81,363 malaria cases were reported.³ Regarding the parasite, in Latin America, malaria is caused mainly by Plasmodium vivax and Plasmodium falciparum parasites, transmitted through infected Anopheles spp. females.¹Anopheles (Nyssorhynchus) darlingi Root 1926 is an important malaria vector in South America with a wide distribution from the south of Mexico to the north of Argentina and from the eastern side of the Andes Mountains to the Atlantic Coast.^4,5,6,7 Specially in Colombia, An. darlingi is also found in the western side of the Andes.⁸

The primary tool for reducing the density of Anopheles spp. mosquito populations is based on neurotoxic insecticides, mainly indoor residual spraying (IRS) and long-lasting insecticide-treated nets (LLITN). Pyrethroids are the primarily used compounds globally.¹ Brazilian and Colombian governmental campaigns deploy pyrethroids as the active ingredient in IRS-based applications and LLITN bednets in malaria-endemic regions. In addition, the organophosphate malathion also probably reaches Anopheles mosquitoes in urban centres because this compound has been used against arboviruses vectors in these regions.⁹ The intense and continuous use of chemicals selects insecticide-resistant mosquitoes, posing a severe threat to malaria control.^10,11 There is a vast literature concerning insecticide resistance and the underlying molecular mechanisms in anophelines from African and Asian countries.^12-19 On the other hand, the status of susceptibility or resistance to insecticides is scarcely known in An. darlingi populations, despite its importance in the cycle of malaria, especially in the Amazonian Region.

Pyrethroids and DDT target the voltage-gated sodium channel (Na _V ) of arthropods, causing the knockdown effect: “paralysis and consecutive death”.^20,21 Single nucleotide substitutions in the Na _V gene related to resistance to the knockdown effect are referred to as kdr mutations.^7,22,23 The substitutions Leu to Phe or Ser in the codon 1014 are the most found kdr mutations in diverse Anopheles species,²⁴ including the neotropical An. albimanus and An. albitarsis s.s.^18,25 Organophosphates (OPs) and carbamates target the acetylcholinesterase enzyme, whereas mutations in the ace-1 coding gene are associated with resistance to OPs insecticide class. The substitution G119S in the ace-1 gene was found in resistant Culex pipiens^26,27 and species of the Anopheles gambiae complex cryptic species.^26,28 Analyses of the nucleotide diversity spanning these mutations are essential to tracking the evolutionary dynamics of arising and dispersal of insecticide resistance mechanisms.^{13,16,29,30,31}

The broad geographic occurrence of An. darlingi, associated with differentiation in behavioural, morphological, and genetic traits,^4,32,33 raised the hypothesis that this species would comprise a group of cryptic species. For instance, there are reports describing variations in the egg (exochorion) morphology,³⁴ wings morphometry,^35,36,37 hematophagy activity behaviour, and lifespan.^38,39 Besides, polymorphism in the banding patterns of polytene chromosomes also suggested cytogenetic variations.^40,41,42 At a molecular level, several markers showed diversity among populations along with geographic variation, such as ITS2 in rDNA,⁴³cytochrome oxidase I gene (COI) in mtDNA,^44,45 microsatellites^46-51 and single nucleotide polymorphisms (SNPs).^4,52,53 Despite these multiple behavioural, morphological, and genetic distinctions, some studies claim that these differences are divergences among populations, not enough to evidence An. darlingi as a complex of cryptic species.^{4,35,36,38,45,47,50,51,54-56}

Genes that influence biological rhythm activities, like mating and locomotion,^57,58 are helpful molecular markers to evidence insect speciation. In this context, the genes timeless (tim) and period (per) are good examples in population genetic studies with Drosophila and also insects of sanitary importance, including Anopheles.^59,60 Analysis of nucleotide variation in tim supported the division of Anopheles (kertesia) cruzii into two cryptic species, where one group occurs in Bahia State (northeast of Brazil) and the other in the Southern and Southeastern Brazilian regions.⁶¹ In the Nyssorhynchus subgenus, An. triannulatus s.l. was divided into two distinct clusters: one comprising An. halophylus and An. triannulatus species C, and the second represented by An. triannulatus s.s.⁶² Variations in the per gene, summed with a series of biochemical and behavioural observations, contributed to separating the species Lutzomyia longipalpis, the primary sandfly vector of Leishmania infantum in Brazil, into several groups related to courtship song.^59,63 Populational studies with other sandflies, Lutzomyia umbratilis, Lutzomyia intermedia, and Lutzomya whitmani corroborated the per gene as an excellent molecular marker for the analysis of speciation processes.^64,65,66

Discrimination among cryptic species of insect vectors has great epidemiological importance once different species may have distinct vector capacity, therefore requiring specific vector control strategies.⁶⁷ This study explored the molecular diversity of genes related to behaviour and insecticide resistance, estimating genetic differentiation in An. darlingi populations from Brazilian and Colombian localities.

MATERIALS AND METHODS

Samples description - We obtained female An. darlingi samples from five Amazonian localities: Estrada do Brasileirinho, Manaus (79 samples), Amazonas State (3º01’16’’S, 59º52’55’’W); Unini River (125 samples) (01º45’46.0’’S, 62º13’39.6’’W) and Jaú River (132 samples) (01º53’2.0’’S, 61º44’31,6’’W), Barcelos, Amazonas State; and Porto Velho (125 samples), Rondônia State (7º18’32’’S, 67º05’42’’W) in Brazil; and from Tagachi-Quibdó (55 samples), Chocó Department (6º7’53.04’’N, 76º26’0.6’’W) in Colombia (Fig. 1). These samples were collected with a Castro aspirator by gently capturing females seeking blooding feeding in protected human landing catches (HLC) during the first hours of the evening.⁶⁸ Mosquitoes were identified based on morphological characters.⁶⁹ Those confirmed as An. darlingi were preserved in silica or ethanol and shipped to the Laboratory for further genomic analyses.

Fig. 1:
Anopheles darlingi sample sites. In red: Tagachi-Quibdó, Choco Department - Colombia. In green, blue, and purple: Unini River, Jaú River, and Manaus, respectively, in Amazonas State. In yellow: Porto Velho, Rondônia State. The Colombian Andean Region is indicated.

We calculated geographic distances among the localities using the free software Qgis, version 2.18.24 (available from: http://www.qgis.org/). We highlight that the collections were not designed specifically for this study, and that the samples were kindly provided for the molecular assays and analyses herein presented. In general, captures occurred indoors, around houses, and in adjacent forest areas.

DNA extraction and sample pooling - Genomic DNA was individually extracted from adult females following⁷⁰ with slight modifications. Mosquitoes were individually macerated in 200 µL TNES buffer [250 mM Tris pH 7.5, 2 M NaCl, 100 mM EDTA, and 2.5% sodium dodecyl sulphate (SDS)] with the addition of 2 µL of 20 mg/mL proteinase K and left for incubation in a 56ºC water bath overnight. After 1 minute centrifugation, we added 100 µL of 5 M NaCl to precipitate the protein content and centrifuged again for 6 min at 15,000 g. The supernatant was transferred to new tubes with a similar volume of 100% isopropanol and centrifuged for 5 min at 15,000 g. The supernatant was discarded, and the DNA pellet was washed with 70% ethanol by 6 min centrifugation at 15,000 g and supernatant discarding. To remove any alcoholic trace, we heated the open tubes for 10 min at 60ºC. Finally, DNA was eluted in TE 0.1X (30 µL) and quantified in a NanoDrop One (ThermoFisher). We made a DNA pool for each population (Manaus, Unini River, Jaú River, Porto Velho and Colombia): 1 µL (around 0,07 µg) of each DNA sample of this respective population: Manaus (N = 79), Unini River (N = 125), Jaú River (N = 132), Porto Velho (N = 125), and Colombia (N = 55). Samples originated in Tagachi - Chocó are referred as Colombia in this paper.

Polymerase chain reaction (PCR) amplification - The primers for Na _V , ace-1, and per fragments were designed for this study, whereas the tim fragment primers were previously available [Supplementary data (Table I)]. PCR amplifications were carried out with one of the three hi-fidelity polymerase kits options: 1) Phusion High-Fidelity PCR Master Mix with GC Buffer (New England, BioLabs); 2) USB Fidelitaq^TM DNA Polymerase (Affymetrix); or 3) Pfu DNA Polymerase (Thermo Scientific), following manufacture instructions and annealing temperature as described in Supplementary data (Table I). All reactions were run in a Veriti Thermocycler (Applied Biosystems), and the PCR products were purified using the magnetic beads kit Agencourt AMPure XP (Beckman Coulter), according to the manufacture instructions. We cloned the purified products with the Kit pJet (Fermentas) into Escherichia coli DH5-α competent cells. The DNA preparations followed the alkaline lysis procedure,⁷¹ and sequencing reactions were performed using the kit BigDye Terminator v3.1 Cycle Sequencing (Applied Biosystems) following standard procedures and submitted to an ABI Prism 3730 (Applied Biosystems) in the DNA Sequencing Facility Platform at Fiocruz.

Sequence analyses - We used the software Geneious 9.1.8⁷² to edit, identify, align, annotate and translate the obtained sequences. The term haplotype is here used to refer to the distinct sequences of the same gene fragment. The polymorphism analyses identified the number of haplotypes (h), polymorphic sites (S), nucleotide diversity (Π), neutral parameter (Ө). We carried out three tests of selective neutrality: Tajima’s D (Tajima 1989), Fu’s Fs (Fu 1997) and Ramos-Onsins and Rozas’ R 2 (2002) with the software DnaSP 5.0.⁷³ For the genetic differentiation analyses, we used the software ProSeq 2.9.1⁷⁴ for calculating fixation index (Fst), number of polymorphic sites (Ss), fixed sites (Sf), and exclusive polymorphic site (Sx, Sy). A Mantel test on geographical [estimated as Ln (km)] and genetic distance [estimated as Fst/(1-Fst)] was performed in Arlequin 3.1.⁷⁵ The correlation coefficient (r) was estimated using 10,000 permutations. To construct phylogenetic relationships and choose the best nucleotide substitution model, we used the JModel Test 2.0.⁷⁶ The maximum likelihood (ML) trees were obtained with MEGA 7.0,⁷⁷ because of its better resolution, and the haplotype networks were achieved using TCS analysis with the software PopART.⁷⁸ Bayesian inferences were carried out using MrBayes 3.2.4.⁷⁹ With this software, we performed for 10-million generation with two parallel searches using nine heated and one cold Markov chain (MCMC). The visualisation and analysis of the MCMC trace files generated through Bayesian phylogenetic inference were performed using Tracer v.1.7.2.⁸⁰ Convergence of the two runs (average standard deviation of split frequencies < 0.01) and likelihood stationarity were checked. Trees were exported to FigTree v1.4.3⁸¹ for visualisation.

RESULTS

We evaluated DNA sequences of fragments of the genes per, tim, Na _V and ace-1, obtained from a total of 516 specimens of An. darlingi from Amazonian localities in Brazil and Colombia. The haplotypes unrevealed of each gene are accessible in the Genbank [sequences and respective accession numbers in the Supplementary data (Table II) (per), Table IV (tim), Table V (Na _V ) and Table VI (ace-1)].

Polymorphism analyses

Biological rhythm genes - We obtained 101 sequences with a 535 bp fragment of the gene per. The fragment included two exons (exon 2 = 22 bp and exon 3 = 437 bp) and an intron (76 bp). We observed 39 variable sites (7%), resulting in 65 haplotypes [Supplementary data (Table II)] with 26 synonymous substitutions and 13 non-synonymous: two and 11, respectively, in exons 2 and 3. Amino acid changes [Supplementary data (Fig. 1)] occurred at the alignment positions: 3 (Leu/Ser), 4 (Asp/Tyr), 15 (Thr/Lys), 18 (Thr/Met), 39 (Gly/Glu), 57 (Gly/Asp), 89 (Gly/Glu), 95 (Arg/Cys), 97 (Val/Ala), 107 (Ser/Leu), 108 (Ser/Asn), 141 (His/Gln), 144 (His/Asn). Supplementary data (Table III) summarises the nucleotide variation for each population, evaluations of the number of haplotypes (H, from 11 to 20), number of polymorphic sites (S, 18 - 23), nucleotide diversity (Π, 0.007 up to 0.011), and neutral parameter (Ө, 0.009 - 0.012). The diversity index was similar in the five populations; nonetheless, the diversity of haplotypes (0.852 up to 0.995), considering the number of sequences at each locality, was higher in Porto Velho (95.2%) and Manaus (89.5%) and lower in Unini River (70.6%), Colombia (56.5%) and Jaú River (52.4%). Most of the haplotypes (89%) were exclusive to one population: Unini River (n = 7), Jaú River (n = 8), Manaus (n = 13), Porto Velho (n = 17) and Colombia (n = 13). Considering all sequences, the per-4 haplotype [Supplementary data (Table II)] was the most frequent, represented by nine sequences: eight in Jaú River and one in Unini River, followed by per-1 with eight sequences: six in Unini River, one in Manaus, and one in Porto Velho.

The 737 bp tim sequenced fragment included two exons (177 and 142 bp, respectively, exons 3 and 4) and two introns (82 and 336 bp, respectively, introns 3 and 4). In the total of 86 sequences, we observed 28 polymorphic sites (~ 3,8%), resulting in 34 haplotypes [Supplementary data (Table IV)], with three non-synonymous substitutions, all in exon 3. Amino acid changes [Supplementary data (Fig. 2)] occurred at three different alignment sites: 10 (Glu/Arg), 41 (Gly/Asp), and 55 (Thr/Ala), all in Colombia. Supplementary data (Table III) outlines nucleotide variation for each population (H: 8-12, S: 8-18, Π: 0.002-0.006 and Ө: 0.003 - 0.007). As observed in per, Brazilian populations showed a similar number of polymorphic sites and haplotype amount in the tim sequences, less diverse than in the Colombian population. The diversity indexes were similar in all populations, and the diversity of haplotypes was higher in Colombia (66.7%), followed by Unini River and Porto Velho (62.5%), Manaus (60%), and Jaú River (38.1%). Most of the haplotypes (76.5%) were exclusively found in one population: Colombia (12), Manaus (4), Unini River (4), Jaú River (3), and Porto Velho (3). The haplotype tim-1 was the most frequent (23 sequences): ten in Jaú River, seven in Manaus, four in Unini River, and two in Porto Velho, followed by Tim-6 (7 sequences): three in Jaú River and Porto Velho and one in Unini River [Supplementary data (Table IV)].

Insecticide resistance genes - The Na _V fragment included two exons (88 and 177 bp, respectively, in exons 20 and 21) and one intron (74 bp). In the 143 obtained sequences, we observed 16 polymorphic sites (~ 13%), resulting in 17 haplotypes [Supplementary data (Table V)]. Amino acid changes [Supplementary data (Fig. 3)] occurred at ten different alignment sites: 5 (Trp/Arg), 10 (Val/Ala), 16 (Ile/Thr), 17 (Pro/Leu), 40 (Ser/Pro), 53 (Asp/Asn), 66 (Ile/Thr), 68 (Arg/Cys), 69 (Phe/Leu), 79 (Asn/Asp). When numbered according to the Musca domestica sodium channel protein, these substitutions are usually W991R, V996A, I1002T, P1003L, S1026P, D1039N, I1052T, R1054C, F1055L, N1065D, respectively, none of which ever related to some known kdr mutation. Supplementary data (Table III) summarises nucleotide variation for each population (S: 14-66, H: 2-7, Π: 0.000-0.001, Ө: 0.000 - 0.003). In general, nucleotide variation was low and similar amongst populations. Unini River presented the highest diversity of haplotypes (22.7%), followed by Manaus (21.1%), Jaú River (18.2%), Porto Velho (14.3%), and Colombia (10.6%). Most haplotypes were exclusively found in one population (88.2%): Colombia (6), Unini River (3), Jaú River (3), Manaus (2), and Porto Velho (1). The haplotype Na _V -1 was the most frequent, with 122 sequences: 60 in Colombia (49.2%), 19 in Jaú River (15.5%), 17 in Unini River (13.9%), and 13 in Manaus and Porto Velho (10.7%). The second most frequent was Na _V -9, with five sequences: three in Manaus (60%) and two in Jaú River (40%) [Supplementary data (Table V)].

The 344 bp sequenced fragment of ace-1 included one exon (exon 5), represented by 99 sequences with 25 polymorphic sites (~ 7%), resulting in 36 haplotypes [Supplementary data (Table VI)]. Amino acid changes occurred at four different alignment sites [Supplementary data (Fig. 4)]: 2 (Val/Ile), 28 (Met/Thr), 55 (Ala/Ser), and 107 (Glu/Val). According to the amino acid numbering in the acetylcholinesterase protein of M. domestica, these substitutions are V57I M83T, A110S, and E162V, respectively, none so far reported related to organophosphates and carbamates resistance. Supplementary data (Table III) outlines nucleotide variation for each population (S: 8-11, H: 8-13, Π: 0.006-0.007, Ө, 0.006-0.009). Porto Velho presented the highest haplotype diversity (70.6%), followed by Manaus (55.6%), Jaú River (54.2%), Unini River, and Colombia (50%). Most of the haplotypes (63.9%) were exclusively found in one population: Colombia (6), Unini River (6), Porto Velho (5), Jaú River (4), and Manaus (3). The haplotype Ace-5 was the most frequent [Supplementary data (Table VI)], represented by 28 sequences: eight in Unini River (28.5%), seven in Jaú River and Manaus (25%), five in Porto Velho (18%), and one in Colombia (3,5%). Ace-2 was the second most frequent haplotype, with eight sequences: two in both Jaú River and Porto Velho and four in Colombia. None of these haplotypes carried the ace-1 mutant allele (119S), providing no evidence of ace-1 duplication.

Selective neutrality test - The neutrality tests indicated that in all cases, values were not significant after Bonferroni’s correction, indicating no obvious departures from neutrality [Supplementary data (Table VII)].

Genetic differentiation

Biological rhythms genes - Analyses with the gene per detected significant F _st values (p < 0.001), higher than 0.15 in all comparisons between Colombia or Jaú River with all other populations [Supplementary data (Table VIII)]. Interestingly, populations from Manaus and Porto Velho, which are around 940 km apart, showed a low F _st (0.027; p > 0.05), while Unini River and Jaú River populations, distant in around 5 km, presented a significant differentiation (F _st = 0.190, p < 0.001). In general, Colombian was the most differentiated population (F _st = 0.208 - 0.475; p < 0.001), followed by Jaú River (F_st = 0.178 - 0.190; p < 0.001) [Supplementary data (Table VIII)]. There were no fixed substitution sites in all comparisons [Supplementary data (Table IX)]. We observed between two and 13 exclusive variable polymorphisms, mostly in the Colombian population. In short, the per fragment differed Colombian from Brazilian populations and separated Jaú River from the other Brazilian ones.

Regarding the tim gene, all comparisons involving An. darlingi populations from Brazil showed low genetic differentiation. However, all populations compared with the Colombian presented significant F _st values (p < 0.001), ranging from 0.347 to 0.389 [Supplementary data (Table VIII)]. No fixed substitution sites were found. Like the per gene, more exclusive polymorphic sites were observed in the comparisons between the Colombian and Brazilian populations. Hence, tim gene analyses also differ An. darlingi populations from Colombia and Brazilian.

Insecticide resistance genes - All Na _V gene fragment comparisons did not show any significant genetic differentiation between Brazilian pairwise populations (F _st from 0 to 0.06; p > 0.05). The Colombian population differed with Unini River (F _st = 0.035, p < 0.01) and with Manaus (F _st = 0.057, p < 0.05) populations [Supplementary data (Table VIII)]. There was one shared polymorphic site between Unini River and Manaus (Supp. Material Table IX), no one fixed substitution site and low exclusive variation in all populations.

In general, the ace-1 gene evidenced low and non-significant Fst values, suggesting low genetic differentiation among Brazilian populations. The exception was the pair Unini River and Porto Velho populations, that although they are nearly 1,130 km apart, showed a significant differentiation (F _st = 0.164, p < 0.001) [Supplementary data (Table VIII)]. All comparisons between the Colombia and a Brazilian populations resulted in high and significant differentiation index values (F _st = 0.054 - 0.330, p < 0.05). There were elevated shared polymorphic sites, no fixed substitution sites, and low exclusive variable sites among all populations [Supplementary data (Table IX)]. Like with the biological rhythms genes, ace-1 analyses differentiated the An. darlingi from Colombia from all Brazilian populations.

Geographic x genetic differentiation - The Mantel test [Supplementary data (Tables X-XI)] did not evidence any statistically significant correlation between genetic and geographical distances among Brazilian populations, therefore not supporting the hypothesys of Isolation by Distance model. Nevertheless, we observed positive correlations for tim (p = 0.024) and with a p near to significance value in ace (p = 0.057) when Colombia was included.

Phylogenetic analysis - The TCS haplotype networks in most cases separated the Colombian from the other populations, except for ace-1 [Supplementary data (Figs 5-8)] . We carried out ML analysis and Bayesian inference of phylogenetic trees. The ML (Fig. 2) and Bayesian phylogenetic trees [Supplementary data (Figs 9-12) considered the JModel Test, with the best-fit models of nucleotide substitution: General Time Reversible model for per fragment, Kimura 2 parameters for tim, Hasegawa-Kishino-Yano model for Na _V gene and Tamura 3 parameters for ace-1. The Bayesian phylogenetic trees of per and tim genes showed a clustering of most Colombian sequences with posterior probabilities of 0.8083 and 0.9967, respectively [Supplementary data (Figs 9, 10)]. No sub-clustering among Brazilian populations was observed. For the insecticide resistance-related genes, most Na _V sequences were grouped homogeneously, showing a low variation level as expected, once it is known that the Na _V is highly conserved, even among distinct taxa. However, ace-1 showed heterogeneous distribution with no supported clusters.

Fig. 2:
phylogenetic trees to per (A), tim (B), Na _V (C) and ace-1 (D). Numbers on nodes represent the percentage bootstrap values on 1000 replications. In green: Unini River (Uni), blue: Jaú River (Jau), yellow: Porto Velho (Pve), purple: Manaus (Man); red: Colombia (Col).

DISCUSSION

The hypothesis that An. darlingi is a complex of cryptic species still needs more robust support, despite several studies evidencing behavioural, morphological, biochemical and genetic differences.^{4,6,35,36,38,44,47,50-52,55,56} In this study, we presented evidence of genetic differentiation among An. darlingi Amazonian populations from Brazil and a population of Tagachí in the department of Chocó - Colombia, with genes commonly used in molecular behaviour studies and markers for insecticide resistance mechanisms.

As a whole, the genes per, tim and ace-1 showed similar polymorphic diversity, with low nucleotide diversity indexes (0.00257-0.01171). The gene per is also related to sexual selection⁸² and is expected to be highly conserved among distinct populations. Indeed, populational studies with other insect disease vectors also showed low nucleotide diversity in the per gene, as in An. gambiae s.l. populations from Burkina Faso, Senegal and Cameroon (Π 0.0031)⁸² and the sandfly Lutzomya umbralitis (Π 0.0031) from Manacapuru, in Amazonas State.⁶⁶ The di-aminoacid repetition Thr-Gly in the Period protein of D. melanogaster is related to thermostability, and the frequency of haplotypes with a specific number of repetitions varies according to geographic origin in the Northern hemisphere (Thr-Gly₁₇ and Thr-Gly₂₀). The Thr-Gly₂₀ repetition block was suggested to be an adaptation to cold weather and high latitudes, while Thr-Gly₁₇ is adapted to hot weather, as they are more frequent in the Mediterranean Region.^83,84,85 Our study detected a block of repetitions with the CAT codon in the An. darlingi per gene, coding for 12 histidines (His₁₂) in all per haplotypes. However, there was a substitution inside this poly-His region in the haplotypes per34 from Colombia (His to Glu) and per63 from Manaus (His to Asn). It will be interesting to further elucidate the phenotypical meaning of this polymorphism.

Polymorphisms found in the tim gene contributed to identifying sibling cryptic species in sympatric samples of the Triannulatus complex, divided into An. triannulatus s.s., An. halophylus and An. triannulatus, in the central part of Brazil.⁶² Herein we showed similar diversity indexes for tim in all evaluated populations, however with high and significative F _st pairwise values differentiating Brazilian and Colombian populations. The tim polymorphism analyses in the main simian and human Plasmodium vectors found in the Atlantic Rain Forest in the Brazilian littoral evidenced the division of An. (Kerteszia) cruzii into one distinct cluster in each Northeastern and Southeastern Regions.⁶¹ To our knowledge, the only previous analyses of this gene in An. darlingi explored genetic differences in the host seek biting behaviour of endophily and exophily in Portuchuelo and Macapá, in the northeastern Amazon, with low but relevant differences between these two groups.⁸⁶ Further studies focusing on the variation of tim in more An. darlingi populations, considering geologic and climate aspects, will help to elucidate those findings.

The Na _V fragment revealed the lowest nucleotide diversity (Π, 0.00042 - 0.00173) compared to the other genes. A previous study on this gene showed low nucleotide variability in Lu. longipalpis and Lutzomyia cruzi, with specific variations in each species.⁸⁷ It is expected because few mutations are permissive in the Na _V protein, given this highly conserved functional role in the physiology of the nervous system.⁸⁸ Unlike per and tim, Na _V and ace-1 presented some similar haplotypes between Brazilian and Colombian mosquitoes suggesting gene flow among them and ancestral polymorphism. Interestingly, the Porto Velho population revealed two highly polymorphic Na _V haplotypes, with all variation observed in the intron. Future studies in that locality, considering collections under distinct circumstances (inside/outside the houses, different moments along the night, etc.), and including samples from other sites will provide better clues whether this polymorphism is linked to different behavioural aspects.

Neurotoxic insecticide usage is the primary strategy in vector control, mainly because of the use of LLITN, which reduced significantly malaria cases in Africa.⁸⁹ Resistance management studies associated with vector population genetics are essential to understand better the origins and dispersion of resistance mechanisms.⁶⁷ The low nucleotide diversity profile in the ace-1 gene that we observed in An. darlingi populations match similar profiles observed among populations of An. gambiae s.l. (Π 0.00634)⁹⁰ and Culex pipiens (Π 0.024).⁹¹ Gene duplication events on ace-1 were described in Anopheles and Culex species.^26,90,91,92 This phenomenon may enhance genetic diversity in this gene, as the different copies may evolve independently. The ace-1 ^R allele (G119S), associated with organophosphate resistance, also results in a high fitness cost to insect development and reproduction.⁹⁰ Continued pressure has selected duplications containing the wild-type allele and an “organophosphate-resistant” allele, thus maintaining insecticide resistance and decreasing deleterious fitness effects.⁹² This same organophosphate resistance has not been described in Amazonian populations of An. darlingi,⁹³ however the ace-1 ^R allele was found in An. albimanus from Peru.⁹³ Although this mutant allele was not found among our samples, the high diversity observed and the presence of at least non-synonymous substitutions suggest that we should monitor variations in the insecticide susceptibility in these populations, especially in those closer to urban centres where organophosphates are used against Aedes mosquitoes.⁹

The classic L1014F kdr mutation is reported in at least 14 Anopheles species, including the Neotropical An. albitarsis s.s. in southern Brazil²⁵ and on An. albimanus in Central and North America,^18,94 selected due to the intense pressure exerted by vector control interventions with pyrethroids in IRS and LLITN²⁴ or in proximity to agricultural areas.²⁵ We did not find any non-synonymous variation in the Na _V gene fragment of An. darlingi. Further studies must consider evaluating the sequence beyond the fragment containing the 1014 site in the Na _V genes, as the substitution Asn to Tyr found in the 1575 Na _V site of An. gambiae, for example.^16,24 Due to the kdr recessive trait, the mutant allele tends to remain under low frequencies for a long time, until it reaches the tipping-point, from which the kdr frequency will increase under exponential speed.⁹⁵ The surveillance of known kdr mutations and screenings in search for non-described variations in the Na _V gene are essential to be investigated in An. darlingi populations before kdr alleles reach levels able to affect the chemical control efficacy. We must also consider that in several regions, An. darlingi populations are exophilic and exophagic,^6,96,97,98 meaning that the mosquito contact with insecticide in IRS and impregnated nets can be minimal. Therefore, future studies should consider testing samples collected indoors and outside the houses. Additionally, recent studies have shown differences between An. darlingi populations from urban and rural areas due to human interventions, deforestation and seasonal climate variations influencing on mosquito distribution. These changes may impact directly on vector control and surveillance strategies.^54,97

A previous study with An. darlingi populations throughout Central and South American localities using sequences of COI fragment evidenced a well distinguished differentiation of the pool of this mitochondrial DNA haplotypes from Latin American and South American populations.⁴⁴ Here, analyses with per, tim and ace-1 fragments grouped separated Brazilian and Colombian An. darlingi populations. The Mantel Test evidenced significant isolation by distance, at least for the tim and ace-1 fragments among Brazilian and the Colombian populations. Our samples from Colombia are from a locality in the western Andean Region (Fig. 1), and a previous study demonstrated that the Andes mountain is a potent barrier that prevents the geneflow between An. darlingi populations from eastern and western sides of these montains.⁵⁶ The lack of gene flow and distinct selection processes through time is an example of a pathway toward speciation, as discussed in previous studies.^{4,37,44,48,50,51} For instance, wing morphological variations observed in An. darlingi between indoor and outdoor mosquitoes in Colombia’s western and eastern regions were attributed to incipient species differentiation.³⁷ The Amazonian hydrographic basin is also an important barrier to the geneflow, as demonstrated by microsatellites with An. darlingi populations from northern and southern Amazonian riversides.⁴⁶ Here, we showed a significant F _st differentiation based on the per gene between the Unini and Jaú populations. Although this was the less distant pair of localities evaluated (50 km), they are communities in distinct riversides.

We added An. darlingi molecular data important for population genetic analyses and polymorphisms in behavioural and insecticide resistance-related genes. These data can contribute to the discussion if An. darlingi is or not a complex of cryptic species, in addition to previous and future studies with a broader sampling, considering both geographic and behavioural parameters. Understanding species complex dynamics has a great epidemiologic significance once they may present distinct vector capacity and different responses to control strategies.⁶⁷

In conclusion - Our results add genetic data about differentiation among An. darlingi populations, considering per, tim, and ace-1 genes. These results alone are not enough to support the hypothesis about the existence of cryptic species in An. darlingi, however highlights genetic differences between Brazilian and Colombian populations. New non-synonymous variations were identified in the ace-1 gene and the classical kdr mutation was not observed in the Na _V gene.

ACKNOWLEDGEMENTS

To the “Rede de Plataformas Tecnológicas Fiocruz (RPT/Fiocruz/RJ) - DNA Sequencing Platform” for their invaluable support in this research. We also thank the dedicated Entomology field professionals who contributed to the collection process. Furthermore, we would like to acknowledge the exceptional work done by Heloisa Diniz (Núcleo de Atividades de Extensão, IOC/Fiocruz) in enhancing the figures presented in this paper.

REFERENCES

• ¹
  WHO - World Health Organization. World Malaria Report 2022. Gêneva: World Health Organization; 2022. 372 pp.
• ²
  PAHO/WHO - Pan American Health Organization/World Health Organization. Epidemiological update: malaria in the Americas. Washington, DC: PAHO/WHO: 2019.
• ³
  INS - Instituto Nacional de Salud. Protocolo de vigilancia de malaria. 2022. Available from: https://www.ins.gov.co/buscador-eventos/Lineamientos/Pro_Malaria%202022.pdf
  » https://www.ins.gov.co/buscador-eventos/Lineamientos/Pro_Malaria%202022.pdf
• ⁴
  Emerson KJ, Conn JE, Bergo ES, Randel MA, Sallum MAM. Brazilian Anopheles darlingi Root (Diptera: Culicidae) clusters by major biogeographical region. PLoS One. 2015; 10(7): 1-15.
• ⁵
  Araújo MS, Andrade AO, dos Santos NAC, Pereira DB, Costa GS, Medeiros de Paulo PF, et al. Brazil's first free-mating laboratory colony of Nyssorhynchus darlingi. Rev Soc Bras Med Trop. 2019; 52: e20190159.
• ⁶
  Campos M, Alonso DP, Conn JE, Vinetz JM, Emerson KJ. Genetic diversity of Nyssorhynchus (Anopheles) darlingi related to biting behaviour in western Amazon. Parasit Vectors. 2019; 12(242): 1-9.
• ⁷
  Lol JC, Castañeda D, Mackenzie-Impoinvil L, Romero CG, Lenhart A, Padilla N. Development of molecular assays to detect target-site mechanisms associated with insecticide resistance in malaria vectors from Latin America. Malar J. 2019; 18: 1-9.
• ⁸
  Montoya-Lerma J, Solarte YA, Giraldo-Calderón GI, Quiñones ML, Ruiz-López F, Wilkerson RC, et al. Malaria vector species in Colombia - A review. Mem Inst Oswaldo Cruz. 2011; 106(Suppl. 1): 223-38.
• ⁹
  Campos KB, Martins AJ, Rodovalho CM, Bellinato DF, Dias LS, Macoris MLG, et al. Assessment of the susceptibility status of Aedes aegypti (Diptera: Culicidae) populations to pyriproxyfen and malathion in a nation-wide monitoring of insecticide resistance performed in Brazil from 2017 to 2018. Parasit Vectors. 2020; 13(1): 531.
• ¹⁰
  Zaim M, Guillet P. Alternative insecticides: an urgent need. Trends Parasitol. 2002; 18(4): 161-3.
• ¹¹
  Coleman M, Hemingway J, Gleave KA, Wiebe A, Gething PW, Moyes CL. Developing global maps of insecticide resistance risk to improve vector control. Malar J. 2017; 16(86): 9.
• ¹²
  Pinto J, Lynd A, Elissa N, Donnelly MJ, Costa C, Gentile G, et al. Co-occurrence of East and West African kdr mutations suggests high levels of resistance to pyrethroid insecticides in Anopheles gambiae from Libreville, Gabon. Med Vet Entomol. 2006; 20(1): 27-32.
• ¹³
  Santolamazza F, Mancini E, Simard F, Qi Y, Tu Z, della Torre A. Insertion polymorphisms of SINE200 retrotransposons within speciation islands of Anopheles gambiae molecular forms. Malar J. 2008; 7(163): 1-10.
• ¹⁴
  Singh OP, Dykes CL, Das MK, Pradhan S, Bhatt RM, Agrawal OP, et al. Presence of two alternative kdr-like mutations, L1014F and L1014S, and a novel mutation, V1010L, in the voltage gated Na+ channel of Anopheles culicifacies from Orissa, India. Malar J. 2010; 9(146): 1-6.
• ¹⁵
  Badolo A, Traore A, Jones CM, Sanou A, Flood L, Guelbeogo WM, et al. Three years of insecticide resistance monitoring in Anopheles gambiae in Burkina Faso: resistance on the rise? Malar J. 2012; 11(232): 1-11.
• ¹⁶
  Jones CM, Liyanapathirana M, Agossa FR, Weetman D, Ranson H, Donnelly MJ, et al. Footprints of positive selection associated with a mutation (N1575Y) in the voltage-gated sodium channel of Anopheles gambiae. PNAS. 2012; 109(17): 6614-9.
• ¹⁷
  Kang S, Jung J, Lee S, Hwang H, Kim W. The polymorphism and the geographical distribution of the knockdown resistance (kdr) of Anopheles sinensis in the Republic of Korea. Malar J. 2012; 11(151): 1-8.
• ¹⁸
  Lol JC, Castellanos ME, Liebman KA, Lenhart A, Pennington PM, Padilla NR. Molecular evidence for historical presence of knock-down resistance in Anopheles albimanus, a key malaria vector in Latin America. Parasit Vectors. 2013; 6(268): 7.
• ¹⁹
  Zouré AA, Badolo A, Francis F. Resistance to insecticides in Anopheles gambiae complex in West Africa: a review of the current situation and the perspectives for malaria control. Int J Trop Insect Sci. 2021; 41: 1-13.
• ²⁰
  Busvine JR. Mechanism of resistance to insecticides in housefly. Nature. 1951; 168: 193-5.
• ²¹
  Milani R. Comportamento mendeliano della resistenza alla azione del DDT e correlazione tra abbattimento e mortalitá in Musca domestica L. Rev Parasitol. 1954; 15: 513-42.
• ²²
  Verhaeghen K, Bortel WV, Trung HD, Sochantha T, Keokenchanh K, Coosemans M. Knockdown resistance in Anopheles vagus, An. sinensis, An. paraliae and An. peditaeniatus populations of the Mekong region. Parasit Vectors. 2010; 3(59): 1-12.
• ²³
  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? Cell Press. 2011; 27(2): 91-8.
• ²⁴
  Silva APB, Santos JMM, Martins AJ. Mutations in the voltage-gated sodium channel gene of anophelines and their association with resistance to pyrethroids - A review. Parasit Vectors. 2014; 7(450): 1-14.
• ²⁵
  Braga TA, Loureiro AC, Lima JBP, Martins AJR. Insecticide resistance in Anopheles Albitarsis s.s from a rice production field, with first record of Kdr mutation in this species. Research Square. 2021. Available from: https://assets.researchsquare.com/files/rs-139005/v1/4a9d7910-4788-412e-b4e1-15eb8bcef6d0.pdf?c=1649852483
  » https://assets.researchsquare.com/files/rs-139005/v1/4a9d7910-4788-412e-b4e1-15eb8bcef6d0.pdf?c=1649852483
• ²⁶
  Weill M, Malcolm C, Chandre F, Mogensen K, Berthomieu A, Marquine M, et al. The unique mutation in ace-1giving high insecticide resistance is easily detectable in mosquito vectors. Insect Mol Biol. 2004; 13(1): 1-7.
• ²⁷
  Tmimi FZ, Faraj C, Bkhache M, Mounaji K, Failloux A. Insecticide resistance and target site mutations (G119S ace-1 and L1014F kdr) of Culex pipiens in Morocco. Parasit Vectors. 2018; 11(51): 1-9.
• ²⁸
  Alout H, Djogbénou L, Berticat C, Chandre F, Weill M. Comparison of Anopheles gambiae and Culex pipiens acetycholinesterase 1 biochemical properties. Comp Biochem Physiol. 2008; 150: 271-7.
• ²⁹
  Pinto J, Lynd A, Vicente JL, Santolamazza F, Randle NP, Gentile G, et al. Multiple origins of knockdown resistance mutations in the Afrotropical mosquito vector Anopheles gambiae. PLoS One. 2007; 2: 1-11.
• ³⁰
  Mitri C, Markianos K, Guelbeogo WM, Bischoff E, Gneme A, Eiglmeier K, et al. The kdr-bearing haplotype and susceptibility to Plasmodium falciparum in Anopheles gambiae: genetic correlation and functional testing. Malar J. 2015; 14(391): 1-11.
• ³¹
  Cosme LV, Gloria-Soria A, Caccone A, Powell JR, Martins AJ. Evolution of kdr haplotypes in worldwide populations of Aedes aegypti: independent origins of the F1534C kdr mutation. PLoS Negl Trop Dis. 2020; 14(4):e0008219.
• ³²
  Consoli RAG, Lourenço-de-Oliveira R. Principais mosquitos de importância sanitária no Brasil. Rio de Janeiro: Fiocruz; 1994. 228 pp.
• ³³
  Forattini OP. Culicidologia médica. Vol. 2. São Paulo: EDUSP; 2002.
• ³⁴
  Almeida F, Suesdek L, Motoki MT, Bergo ES, Sallum MAM. Morphometric comparisons of the scanning electron micrographs of the eggs of Anopheles (Nyssorhynchus) darlingi Root (Diptera: Culicidae). Acta Trop. 2014; 139: 115-22.
• ³⁵
  Manguin S, Wilkerson RC, Conn JE, Rubio-Palis Y, Danoff-Burg JA, Roberts DR. Populations structure of the primary malaria vector in South America, Anopheles darlingi using isozyme, random amplified polymorphic DNA, internal transcribed spacer 2, and morphologic markers. Am J Trop Med. 1999; 60(3): 364-76.
• ³⁶
  Motoki MT, Suesdek L, Bergo ES, Sallum MAM. Wing geometry of Anopheles darlingi Root (Diptera: Culicidae) in five major Brazilian ecoregions. Infect Genet Evol. 2012; 12(6): 1246-52.
• ³⁷
  Pacheco MA, González R, Brochero HL. Anopheles darlingi Root 1926 (Diptera:Culicidae): varioaciones mormétricas em alas y patas de poblaciones de Colombia. Biomédica. 2017; 37(2): 124-34.
• ³⁸
  Rosa-Freitas MG, Broomfield G, Priestman A, Milligan PJ, Momen H, Molyneux DH. Cuticular hydrocarbons, isoenzymes and behavior of three populations of Anopheles darlingi from Brazil. J Am Mosq Control Assoc. 1992; 8(4): 357-66.
• ³⁹
  Chu VM, Sallum MAM, Moore TE, Lainhart W, Schlichting CD, Conn JE. Regional variation in life historytraits and plastic responses to temperature of the major malaria vector Nyssorhynchus darlingi in Brazil. Sci Rep. 2019; 9(5356): 1-11.
• ⁴⁰
  Schreiber G, Guedes AS. Cytological aspects of the taxonomy of Anophelines (subgenus Nyssorhynchus). Bull World Health Organ. 1961; 24(4-5): 657-8.
• ⁴¹
  Kreutzer R, Kitzmiller J, Ferreira E. Inversion polymorphismin the salivary gland chromossomes of Anopheles darlingi Root. Mosq News. 1972; 32(4): 555-65.
• ⁴²
  Tadei WP, Dos Santos JMM, Rabbani MG. Biologia de anofelinos amazônicos. V. Polimorfismo cromossômico de Anopheles darlingi Root (Diptera, Culicidae). Acta Amaz. 1982; 12: 353-69.
• ⁴³
  Malafronte RS, Marrelli MT, Marinotti O. Analysis of ITS2 DNA sequences from Brazilian Anopheles darlingi (Diptera: Culicidae). J Med Entomol. 1999; 36(5): 631-4.
• ⁴⁴
  Mirabello L, Conn JE. Molecular population genetics of the malaria vector Anopheles darlingi in Central and South America. Heredity (Edinb). 2006; 96: 311-21.
• ⁴⁵
  Pedro PM, Sallum MAM. Spatial expansion and population structure of the neotropical malaria vector, Anopheles darlingi (Diptera: Culicidae). Phylogeography a Neotrop Malar Vector. 2009; 97: 854-66.
• ⁴⁶
  Conn JE, Vineis JH, Bollback JP, Onyabe DY, Wilkerson RC, Póvoa MM. Population structure of the malaria vector Anopheles darlingi in a malaria-endemic region of eastern amazonian Brazil. Am J Trop Hyg. 2006; 74(5): 798-806.
• ⁴⁷
  Scarpassa VM, Conn JE. Population genetic structure of the major malaria vector Anopheles darlingi (Diptera: Culicidae) from the Brazilian Amazon, using microsatellite markers. Mem Inst Oswaldo Cruz. 2007; 102(3): 319-27.
• ⁴⁸
  Mirabello L, Vineis JH, Yanoviak SP, Scarpassa VM, Póvoa MM, Padilla N, et al. Microsatellite data suggest significant population structure and differentiation within the malaria vector Anopheles darlingi in Central and South America. BMC Ecol. 2008; 8(3): 1-15.
• ⁴⁹
  Gutiérrez LA, Gómez GF, González JJ, Castro MI, Luckhart S, Conn JE, et al. Microgeographic genetic variation of the malaria vector Anopheles darlingi Root (Diptera: Culicidae) from Córdoba and Antioquia, Colombia. Am J Trop Med Hyg. 2010; 83(1): 38-47.
• ⁵⁰
  Angêlla AF, Salgueiro P, Gil LHS, Vicente JL, Pinto J, Ribolla PEM. Seasonal genetic partitioning in the neotropical malaria vector, Anopheles darlingi. Malar J. 2014; 13(203): 1-10.
• ⁵¹
  Rosero CY, Jaramillo GI, Gonzalez R, Cardenas H. Genetic differentiation of Colombian populations of Anopheles darlingi Root (Diptera: Culicidae). Neotrop Entomol. 2017; 46(5): 487-98.
• ⁵²
  Prussing C, Emerson KJ, Bickersmith AS, Sallum MAM, Conn JE. Minimal genetic differentiation of the malaria vector Nyssorhynchus darlingi associated with forest cover level in Amazonian Brazil. PLoS One. 2019; 14(11): 1-16.
• ⁵³
  Chu VM, Sallum MAM, Moore TE, Emerson KJ, Schlichting CD, Conn JE. Evidence for family-level variation of phenotypic traits in response to temperature of Brazilian Nyssorhynchus darlingi. Parasit Vectors. 2020; 13(55): 1-15.
• ⁵⁴
  Campos M, Conn JE, Alomso DP, Vinetz JM, Emerson KJ, Ribolla PEM. Microgeographical structure in the major Neotropical malaria vector Anopheles darlingi using microsatellites and SNP markers. Parasit Vectors. 2017; 10(76): 1-8.
• ⁵⁵
  Lainhart W, Bickersmith SA, Nadler KJ, Moreno M, Saavedra MP, Chu VM, et al. Evidence for temporal population replacement and the signature of ecological adaptation in a major Neotropical malaria vector in Amazonian Peru. Malar J. 2015; 14(375): 17.
• ⁵⁶
  González R, Wilkerson R, Suárez MF, García F, Gallego G, Cárdenas H, et al. A population genetics study of Anopheles darlingi (Diptera: Culicidae) from Colombia based on random amplified polymorphic DNA-polymerase chain reaction and amplified fragment lenght polymorphism markers. Mem Inst Oswaldo Cruz. 2007; 102(3): 255-62.
• ⁵⁷
  Sakai T, Ishida N. Circadian rhythms of female mating activity governed by clock genes in Drosophila. Proc Natl Acad Sci. 2001; 98(16): 9221-5.
• ⁵⁸
  Tauber E, Roe H, Costa R, Hennessy M, Kyriacou CP. Temporal mating isolation driven by a behavioural gene in Drosophila. Curr Biol. 2003; 13: 140-5.
• ⁵⁹
  Araki AS, Vogoder FM, Bauzer LGSR, Ferreira GEM, Souza NA, Araújo IB, et al. Molecular and behavioural differentiation among Brazilian populations of Lutzomya longipalpis (Diptera:Psychodidae:Phlebotominae). PLoS Negl Trop Dis. 2009; 3: 1-12.
• ⁶⁰
  Bauzer LG, Souza NA, Ward RD, Kyriacou CP, Peixoto AA. The period gene and genetic differentiation between three Brazilian populations of Lutzomyia longipalpis. Insect Mol Biol. 2002; 11: 315-23.
• ⁶¹
  Rona LD, Carvalho-Pinto CJ, Gentile C, Grisard EC, Peixoto AA. Assessing the molecular divergence between Anopheles (Kerteszia) cruzii populations from Brazil using the timeless gene: further evidence of a species complex. Malar J. 2009; 8(60): 1-10.
• ⁶²
  Silva-do-Nascimento T, Pitaluga LDR, Peixoto AA, Lourenço-de-Oliveira R. Molecular divergence in the timeless and cpr genes among three sympatric cryptic species of the Anopheles triannulatus complex. Mem Inst Oswaldo Cruz. 2011; 106(Suppl. 1): 218-22.
• ⁶³
  Araki AS, Ferreira GEM, Mazzoni CJ, Souza NA, Machado RC, Bruno RV, et al. Multilocus analysis of divergence and introgression in sympatric and allopatric sibling species of then Lutzomyia longipalpis complex in Brazil. PLoS Negl Trop Dis. 2013; 7(10): 1-13.
• ⁶⁴
  Sawyer LA, Sandrelli F, Pasetto C, Peixoto AA, Rosato E, Costa R, et al. The period gene Thr-Gly polymorphism in Australian and African Drosophila melanogaster populations: implications for selection. Genetics. 2006; 174(1): 465-80.
• ⁶⁵
  Mazzoni CJ, Souza NA, Andrade-Coelho C, Kyriacou CP, Peixoto AA. Molecular polymorphism, differentiation and introgression in the period gene between Lutzomyia intermedia and Lutzomyia whitmani. BMC Evol Biol. 2006; 6(85): 1-11.
• ⁶⁶
  Freitas MTS, Ríos-Velasquez CM, Silva LG, Costa CRL, Marcelino A, Leal-Balbino TC, et al. Phenotypic and genotypic variations among three allopatric populations of Lutzomyia umbratilis, main vector of Leishmania guyanensis. Parasit Vectors. 2015; 8(448): 1-10.
• ⁶⁷
  Bickford D, Lohman DJ, Sodhi NS, Ng PKL, Meier R, Winker K, et al. Cryptic species as a window on diversity and conservation. Trends Ecol Evol. 2007; 22(3): 148-55.
• ⁶⁸
  WHO - World Health Organization. Malaria entomology and vector control - Learner's Guide. Geneva: World Health Organization; 2002. 107 pp.
• ⁶⁹
  Faran ME, Linthicum KJ. A handbook of the Amazonian species of Anopheles (Nyssorhynchus) (Diptera:Culicidae). Mosq Syst. 1981; 13(1): 1-81.
• ⁷⁰
  Martins J, Solomon SE, Mikheyev AS, Mueller UG, Ortiz A, Bacci M. Nuclear mitochondrial-like sequences in ants: evidence from Atta cephalotes (Formicidae: Attini). Insect Mol Biol. 2007; 16(6): 777-84.
• ⁷¹
  Sambrook J, Russel DW. Molecular cloning. A laboratory manual. New York: Cold Spring Harbor; 2001.
• ⁷²
  Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012; 28(12): 1647-9.
• ⁷³
  Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009; 25(11): 1451-2.
• ⁷⁴
  Filatov DA. Processing and population genetic analysis of multigenic datasets with ProSeq3 software. Bioinformatics. 2009; 25(23): 3189-90.
• ⁷⁵
  Excoffier L, Laval G, Schneider S. Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005; 1: 47-50.
• ⁷⁶
  Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008; 25(7): 1253-6.
• ⁷⁷
  Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016; 3(7): 1870-74.
• ⁷⁸
  Leigh JW, Bryant D. PopART: full-feature software for haplotype network construction. Methods Ecol Evol. 2015; 6(9): 1110-6.
• ⁷⁹
  Ronquist F, Teslenko M, Van der Mark P, Ayres D, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012; 61: 539-42.
• ⁸⁰
  Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. Posterior summarisation in Bayesian phylogenetics using Tracer 1.7. Syst Biol. 2018; 67(5): 901-4.
• ⁸¹
  Rambaut A. FigTree v1.4.4 2006-2018: tree figure drawing tool. 2018. Available from: http://tree.bio.ed.ac.uk/software/figtree/
  » http://tree.bio.ed.ac.uk/software/figtree/
• ⁸²
  Morlais I, Poncon N, Simard F, Cohuet A, Fontenille D. Intraspecific nucleotide variant in Anopheles gambiae: new insights into biology of malaria vectors. Am J Trop Med Hyg. 2004; 71(6): 795-802.
• ⁸³
  Costa R, Peixoto AA, Barbujani G, Kyriacou CP. A latitudinal cline in a Drosophila clock gene. Proc R Soc Biol Sci. 1992; 250(1327): 43-9.
• ⁸⁴
  Rosato E, Peixoto A, Costa R, Kyriacou CP. Linkage disequilibrium, mutational analysis and natural selection in the repetitive region of the clock gene, period, in Drosophila melanogaster. Genet Res Camb. 1997; 69: 89-99.
• ⁸⁵
  Kyriacou CP, Peixoto AA, Sandrelli F, Costa R, Tauber E. Clines in clock genes: fine-tuning circadian rhythms to the environment. Trends Genet. 2007; 24(3): 124-32.
• ⁸⁶
  Bottino RC. Genética de populações de Anopheles darlingi e Anopheles marajoara (Diptera: Culicidae), utilizando o gene timeless como marcador molecular [Dissertação de Mestrado em Biologia Parasitária]. Rio de Janeiro: Instituto Oswaldo Cruz; 2007.
• ⁸⁷
  Lins RMMA, Souza NA, Brazil RP, Maingon RDC, Peixoto AA. Fixed differrences in the paralytic gene define two lineages within the Lutzomyia longipalpis complex producing different types of courtship songs. PLoS One. 2012; 7(9): 1-8.
• ⁸⁸
  Ffrench-Constant RH, Pittendrigh B, Vaughan A, Anthony N. Why are there so few resistance-associated mutations in insecticide target genes? Philos Trans R Soc Lond B Biol Sci. 1998; 353(1376): 1685-93.
• ⁸⁹
  Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015; 526(7572): 207-11.
• ⁹⁰
  Djogbénou L, Chandre F, Berthomieu A, Dabiré R, Koffi A, Alout H, et al. Evidence of introgression of the ace-1R mutation and of the ace-1 duplication in West African Anopheles gambiae s. s. PLoS One. 2008; 5: 1-7.
• ⁹¹
  Labbé P, Berthomieu A, Berticat C, Alout H, Raymond M, Lenormand T, et al. Independent duplications of the acetylcholinesterase gene conferring insecticide resistance in the Mosquito Culex pipiens. Mol Biol Evol. 2007; 24(4): 1056-67.
• ⁹²
  Djogbénou L, Labbé P, Chandre F, Pasteur N, Weill M. Ace-1 duplication in Anopheles gambiae: a challenge for malaria control. Malar J. 2009; 8(70): 1-6.
• ⁹³
  Liebman KA, Pinto J, Valle J, Palomino M, Vizcaino L, Brogdon W, et al. Novel mutations on the ace-1 gene of the malaria vector Anopheles albimanus provide evidence for balancing selection in an area of high insecticide resistance in Peru. Malar J. 2015; 14(74): 1-10.
• ⁹⁴
  Orjuela LI, Alvarez-Diaz DA, Morales JÁ, Grisales N, Ahumada ML, Venegas J, et al. Absence of knockdown mutations in pyrethroid and DDT resistant populations of the main malaria vectors in Colombia. Malar J. 2019; 18(384): 1-9.
• ⁹⁵
  WHO - World Health Organization. Global plan for insecticide resistance management in malaria vectors. Geneva: World Health Organization; 2012. 32 pp.
• ⁹⁶
  Hiwat H, Bretas G. Ecology of Anopheles darlingi Root with respect to vector importance: a review. Parasit Vectors. 2011; 4(177): 1-13.
• ⁹⁷
  Moreno M, Saavedra MP, Bickersmith SA, Lainhart W, Tong C, Alava F, et al. Implications for changes in Anopheles darlingi biting behaviour in three communities in the peri-Iquitos region of Amazonian Peru. Malar J. 2015; 14(290): 1-11.
• ⁹⁸
  Vezenegho SB, Adde A, de Santi VP, Issaly J, Carinci R, Gaborit P, et al. High malaria transmission in a forested malaria focus in French Guiana: How can exophagic Anopheles darlingi thwart vector control and prevention measures? Mem Inst Oswaldo Cruz. 2016; 111(9): 561-9.

Publication Dates

• Publication in this collection
  10 July 2023
• Date of issue
  2023

History

• Received
  06 July 2022
• Accepted
  22 May 2023