Genomic population structure and insecticide resistance mechanisms in the malaria vector An. coluzzii across contrasting bioclimatic zones in West Africa

Genomic population structure and insecticide resistance mechanisms in the malaria vector An. coluzzii across contrasting bioclimatic zones in West Africa | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Genomic population structure and insecticide resistance mechanisms in the malaria vector An. coluzzii across contrasting bioclimatic zones in West Africa Enock K. Amoako, Kelly L. Bennett, Anastasia Hernandez-Koutoucheva, and 10 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-7878288/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 07 Jan, 2026 Read the published version in BMC Genomics → Version 1 posted 11 You are reading this latest preprint version Abstract Environmental barriers influencing the movement of insect vectors can govern adaptive gene flow, including the dispersal of insecticide resistance mechanisms that compromise population control. We sought to understand population connectivity of the major malaria vector, An. coluzzii , across the different bioclimatic zones of West Africa using SNPs from whole genomes and inversion karyotypes previously associated with environmental adaptation. We identified restricted gene flow between populations from northern savannah and southern forested regions. Using Ghana as a case study, we found marked differences in insecticide resistance profiles across the different bioclimatic zones suggesting that population connectivity impacts on adaptive allele sharing. Greater evidence for target site pyrethroid and metabolic cross-resistance in the North reflects differences in insecticide use across the country. We also observed distinct resistance mechanisms in the coastal region of Greater Accra which may result from intense urban agricultural activity. Overall, findings suggest that environmental conditions restrict An. coluzzii gene flow to impact the geographical distribution of molecular insecticide resistance. Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Introduction Insect vectors with wide geographical distributions are subject to varying selection pressures driven by diverse ecological conditions. Both geographical barriers and ecological selection on insect vectors can promote variation in physiological and behavioral traits that impact human disease transmission. For example, transmission can be impacted by differences in geographical species distribution, blood-feeding behavior, dispersal and habitat use[ 1 ]. Ecological divergence often impacts on gene flow to form genetically differentiated ecotypes or species with different abilities to act as disease vectors e.g., tsetse flies that are vectors for sleeping sickness[ 2 ] or mosquitoes that transmit malaria[ 3 ]. Differences in population connectivity underpinned by vector ecology also impacts the sharing of adaptive alleles across geographical space. For example, these could be adaptations which promote behavioral avoidance or resistance mechanisms, such as insecticide target site or metabolic resistance, which compromise the population control strategies public health relies on[ 4 ]. There is therefore a need to understand how geography, ecology and anthropogenetic selection pressures, such as the application of insecticides, interact to influence the effectiveness of disease reduction strategies. Ghana, with its rich tapestry of ecological landscapes ranging from coastal Savannah, through dense tropical rainforests, to Sahelian Savannah is representative of the north to south environmental gradient that spans the countries of central West Africa. Its diverse ecological conditions impact the distribution of disease vectors, including major malaria vectors in the Anopheles mosquito genus[ 5 ]. Malaria burden and seasonal peaks in incidence vary markedly across Ghana resulting in ~ 500,000 infections a year[ 6 ]. The key malaria vectors in central West Africa include members of the Anopheles gambiae s.l. complex and Anopheles funestus s.l species group. Within the former is the major malaria vector An. coluzzii . Due to its ecological plasticity, it has a widespread distribution across Ghana including in the northern Savannah where the malaria burden is highest[ 7 , 8 ]. It also exhibits considerable variation in physiological and behavioral traits across the different bioclimatic zones of the West African subregion. These include differences in mating behaviour[ 9 ], dispersal strategies[ 10 ] and phenotypes promoting thermal tolerance including entering into a state of dormancy during dry conditions, i.e., aestivation[ 11 ]. An. coluzzii also showcases genomic variability across contrasting environmental conditions. For example, a striking cline in 2La and 2Rb inversion frequencies across the bioclimatic zones of West and Central Africa have been observed in Nigeria, Burkina Faso and Cameroon[ 12 – 15 ]. While populations in arid, savanna regions exhibit higher frequencies of the inverted 2La and 2Rb karyotypes, those in humid, forested areas show lower frequencies[ 16 , 17 ]. Experimental work has also associated these inversion karyotypes with feeding behavior[ 18 ] and thermal tolerance[ 19 ]. Furthermore, these karyotypes can impact resistance to Plasmodium parasite infection[ 20 ] and to insecticides used for population control[ 21 , 22 ], highlighting that environmental adaptation can impact directly on factors influencing malaria transmission. A major hurdle to malaria control is the rapid evolution of insecticide resistance in vector populations[ 23 ] which undermines the effectiveness of control tools including insecticide-treated nets (ITNs) and indoor residual spraying (IRS)[ 24 , 25 ]. Resistance mechanisms include target site resistance and metabolic resistance (e.g., heightened activity of detoxification enzymes). Target site resistance is conferred by non-synonymous substitutions in the voltage-gated sodium channel ( Vgsc ) gene[ 26 , 27 ], coding for the protein targeted by widely-used pyrethroid insecticides, the acetylcholinesterase ( Ace-1 ) gene[ 28 , 29 ] which codes for the protein target of carbamates and organophosphates, and the GABA receptor subunit Rdl gene, coding for the target of dieldrin[ 30 , 31 ]. Metabolic resistance is provided by increased activity of enzyme products determined by key gene classes such as cytochrome P450 monooxygenases (P450s)[ 32 , 33 ], esterases including carboxylesterases (Coe)[ 34 – 36 ], and glutathione S-transferases (Gst) [ 36 ]and uridine diphosphate (UDP)-glycosyltransferases[ 37 ]. Research conducted in Burkina Faso and Ghana has revealed that An. coluzzii populations in urban areas often exhibit higher levels of insecticide resistance phenotypes compared to rural populations[ 38 , 39 ]. The urban/rural disparity has been attributed to differential exposure to insecticides used in agriculture and public health in peri-urban and urban settings. The distribution of insecticide resistance mechanisms is therefore expected to vary across geographical space and with environmental characteristics. Moreover, resistance loci like target-site substitutions and copy number variants at metabolic resistance loci vary geographically, often reflecting localized selection pressures[ 40 – 42 ]. Furthermore, scans for signals of recent selection have highlighted potential differences in insecticide resistance loci including Vgsc and cytochrome P450’s between forest and non-forest ecosystems in southern Ghana, suggesting differential selection pressures in the contrasting bioclimates[ 43 ]. Understanding the geographic distribution of these resistance mechanisms and how environmental factors influence the sharing of adaptive alleles is crucial for designing targeted and effective control strategies. Advances in genomics have revolutionized our ability to investigate population structure and molecular insecticide resistance mechanisms in malaria vectors at high resolution[ 44 – 46 ]. Here we use whole-genome sequencing of samples collected in Ghana, alongside data available in the MalariaGEN Vector Observatory ( https://www.malariagen.net/vobs/ ), to elucidate population structure across the different bioclimatic zones of West Africa and investigate if ecological variation promotes population structure in An. coluzzii . Using SNP, haplotype and CNV variation data, we also explore how molecular insecticide resistance varies across geographical space using Ghana as a case study to determine if connected populations are more likely to share adaptive alleles resulting from the selection pressures imposed by front-line public health. Methods Mosquito sampling and identification Adult wild caught mosquitoes were sampled across Ghana from different ecological zones in cross-sectional studies from 2016 to 2018. Mosquito sampling was done in Navrongo, Upper East Region of the northern Savannah, Adansi in the Ashanti Region and Koforidua in the Eastern Region of the middle transitional forest, Madina in Greater Accra and Takoradi in Western Region of the coastal Savannah of Ghana. Anopheline mosquitoes were caught using Human Landing Catches (HLCs) at night, CDC light traps hanged overnight, and pyrethroid spray collections of resting mosquitoes in the early mornings. Adult female An. gambiae s.l were identified using morphological keys[ 47 ] and stored singly in 70% ethanol in 96 well PCR plates. These were then shipped to the Wellcome Sanger Institute Genomic Surveillance Unit. DNA was extracted using the Qiagen DNeasy blood and Tissue Kit (Qiagen Science, MD, USA) according to the manufacturer’s instructions before sequencing and genotyping. Genomic data from other West African countries were openly available through the Vector Observatory ( https://www.malariagen.net/vobs/ ), a collaborative project to obtain genomics data of Anopheles mosquitoes for malaria control. Sequencing and variant calling Using previously described protocols described by the Anopheles gambiae 1000 Genomes (Ag1000G) phase 3 project (The Anopheles gambiae 1000 Genomes Consortium 2018, 2020), Illumina paired-end sequencing was performed with HiSeq 2000 and HiSeq X technologies. Briefly, 100-150bp sequencing reads were aligned to the AgamP4 PEST reference genome using BWA version 0.7.15[ 48 ] and single nucleotide polymorphisms (SNP) called using GATK version 3.7.0[ 49 ]. The resulting data was quality controlled to only include individuals with ≥ 10X overall median coverage and with data across greater than 50% of the reference genome. Samples identified as cross-contaminated by a percentage of ≥ 4.5% were excluded as defined by Ag1000G phase 3 project protocols. Sites where SNP calling and genotyping were expected to be unreliable based on previous analyses of Mendelian inheritance in laboratory crosses were also excluded from analysis. Haplotypes were phased using a WDL implementation of read-backed phasing using WhatsHap v1.0[ 50 ] and statistical phasing using SHAPEIT v4.2.1[ 51 ] ( https://github.com/malariagen/pipelines ). Copy number variants (CNVs) for each individual were called based on the copy number state inferred across windows of the genome. The copy number state was normalised using a Gaussian hidden Markov model (HMM) implemented by hmmlearn ( https://github.com/hmmlearn/hmmlearn ) as described previously[ 52 ]. A CNV was called across regions where at least five contiguous genomic windows had a predicted copy number > 2 or > 1 for the X chromosome in males. To increase reliability, only CNV calls with a high likelihood > 1000 and low coverage variance < 0.35 based on the HMM were retained. All preceding analyses were performed with the malariagen_data python package. Taxonomic status Taxa were provisionally assigned to the samples using Ancestry Informative Markers (AIMS). These are a set of SNPs previously described by the Anopheles gambiae 1000 Genomes project as exclusive to each taxonomic group based on data generated from the Anopheles 16 genomes project[ 53 ]. A set of 2612 markers were used to differentiate sister species An. gambiae / An. coluzzii from the more divergent An. arabiensis , identified when the fraction of arabiensis-like alleles was > 0.6. A set of 700 AIMs were used to differentiate An. gambiae from An. coluzzii with samples scored as An. gambiae when the fraction of coluzzii-like calls was 0.9. Individuals in-between these fractions represent other taxa. The taxonomic status of individuals was then confirmed with both PCA and an unrooted Neighbour-joining tree using 100,000 biallelic SNPs evenly spread across chromosome three because this region is unaffected by structural rearrangements such as inversions. Chosen SNPs had a minor allele frequency greater than 0.2% and no missing data. The Neighbour-Joining tree was constructed using city block distance. Population structure To investigate population structure, Principal Component Analysis (PCA) dimensionality reduction and Neighbour-joining trees were constructed for An. coluzzii using the 3L chromosome as applied for taxonomic analysis. We also performed PCA focused on the allele counts at two inverted regions of the genome previously shown to segregate among An. coluzzii from different environments[ 16 , 54 ]. These regions are located on the 2L (2L:20,528,089 − 42,165,182) and 2R (2R:19,444,433 − 26,313,071) chromosomes. Inversion karyotype frequencies within each population cohort were assessed by typing individuals for inversion status using correlated tag SNPs as previously described[ 55 ]. We also assessed genomic differentiation among population cohorts using Hudson’s pairwise F ST using the 3L chromosome arm[ 56 ]. To explore the genetic diversity and demography of populations, informative summary statistics were then calculated including Nucleotide diversity (θπ), Watterson's theta ( θ W ) and Tajima’s D using cohorts with the malariagen_data python package. Statistics were only calculated for population cohorts with a minimum of ten individuals. Insecticide resistance To investigate the presence of known target site mutations associated with insecticide resistance we calculated amino acid substitution frequencies at genomic sites of interest for each population cohort. These were based on the occurrence of non-synonymous SNPs at an appreciable frequency present at greater than 5%. Regions included the voltage-gated sodium channel (Vgsc; AGAP004707) as the target gene of pyrethroids, the glutathione S-transferase gene which confers resistance to DDT ( Gste2 ; AGAP009194), the Resistance to dieldrin gene ( Rdl ; AGAP006028) and the organophosphate target gene, acetylcholinesterase ( Ace1 ; AGAP001356). We also calculated the frequencies of copy number variants (CNVs) for genes associated with insecticide resistance and present at greater than 5% frequency. These included the Cytochrome P450 gene’s (AGAP002862-AGAP002870, AGAP000818, AGAP008212-AGAP008219), carboxylesterases (AGAP006228, AGAP006723-AGAP006728), Ace1 (AGAP001356) and Gste2 (AGAP009194). Selection scans To identify novel regions of the genome under selection we utilised the H12 homozygosity statistic across windows of the genome[ 57 ] with statistical peaks representing either a hard or soft sweep. The statistic was calculated using an optimal size of 1500 windows. This was identified by plotting the distribution of H12 values across different window sizes and identifying when values fell below 0.1 for the 95th percentile. Diplotype clustering To investigate which variants are associated with clusters of diplotypes under selection (i.e., regions of diploid genotypes), hierarchical clustering was performed using city block genetic distance and complete linkage using the malariagen_python package[ 34 ]. Observed amino acid substitutions and CNV variants were plotted onto the resulting dendrogram. Variants uniquely associated with clusters of diplotypes with low heterozygosity are candidates under selection. Results Population Sampling Anopheles gambiae sensu lato samples were collected from 2016 to 2018 across the arid savannah region of the Upper East in northern Ghana. Of 1663 samples sent for whole genome sequencing, 1473 individuals passed quality controls as defined by the Anopheles gambiae 1000 Genomes (AG1000G) project[ 44 , 58 ]. From these samples, 1324 mosquitoes were identified as An. coluzzii through PCA and AIMs analysis generated during data curation[ 44 , 58 ] (Table 1 ). In addition, we included 486 An. coluzzii genomes from a previous study in Ghana[ 40 ] (Table 1 ). The sample sets included mosquitoes collected from 2012 to 2018 across five administrative districts located in the wet humid and deciduous forest region of South Ghana (Fig. 1 ,, Fig. S1 , Supplementary Table 1). For comparison across the West African region, we used data made publicly available in the AG1000G Phase 3 resource ( https://malariagen.github.io/vector-data/ag3/ag3.0.html ). The sample sets included An. coluzzii collected from 2004 to 2014, originating from the dry savannah regions of Burkina Faso (n = 135) and Mali (n = 90) or from the wet humid region of southern Cote d’Ivoire (n = 80) (Fig. 1 , Fig. S1 ). Sample sets were classified for analysis as a population cohort based on the administrative district and year of collection. Samples sequenced during this study had a median coverage of 35X which generated 162,714,957 SNPs on alignment to the AgamP3 genome. Of these SNPs, 52,946,551 were biallelic and segregated within the samples. Table 1 The number of An. coluzzii collected from each sampling location in Ghana as part of the present study and details of other previously published datasets from West Africa included in analysis[ 40 , 44 , 58 ]. Country Location Latitude Longitude Year No. samples Study ID Ghana Upper East Region 10.85 -1.11 2018 563 This study 2017 557 This study 2016 204 This study Ashanti Region 6.18 -1.62 2018 274 Lucas et al. 2023 Central Region 6.01 -1.87 2018 148 Lucas et al. 2023 5.61 -1.55 2012 14 AG1000G Phase 3 Eastern Region 6.09 -0.26 2012 1 AG1000G Phase 3 Greater Accra Region 5.67 -0.22 2012 14 AG1000G Phase 3 Western Region 4.91 -1.77 2012 24 AG1000G Phase 3 Burkina Faso Hauts-Bassins 11.151 -4.235 2012 82 AG1000G Phase 3 2014 53 AG1000G Phase 3 Mali Koulikouro 12.68 -7.84 2004 11 AG1000G Phase 3 2014 27 AG1000G Phase 3 Segou 13.2 -6.13 2004 25 AG1000G Phase 3 Sikasso 10.83 -7.81 2012 27 AG1000G Phase 3 Cote d'Ivoire Lagunes 5.898 -4.823 2012 80 AG1000G Phase 3 Geographical population Structure of An. coluzzii in Ghana Population connectivity and consequent gene flow can influence the spread of adaptive alleles associated with insecticide resistance or vector competence for disease transmission. To first assess population connectivity of An. coluzzii across the different bioclimatic zones of Ghana, we implemented both principal components analysis (PCA) and neighbor-joining trees (NJT) using 100,000 single nucleotide polymorphisms (SNPs) spanning the 3L chromosome region (3L:15,000,000–41,000,000), which is free from structural rearrangements that could bias inference. Both the principal component analysis (PCA) and neighbor-joining tree (NJT) presented two major clusters of An. coluzzii (Fig. 2 A and Fig. S2 ). One cluster included samples from the northern arid savannah region of the Upper East admin region only. The second included all populations from southern Ghana including the wet and humid deciduous forest regions of Greater Accra, Ashanti, Central, Eastern and Western Ghana. Although some southern comparisons were from a different collection time to the northern cohorts, two southern cohorts from Ashanti and the Central region were sampled in the same year. Furthermore, since an older 2012 cohort from the Central region clusters with a later cohort from the same region in 2019, we can observe that the pattern of population structure between the north and south of Ghana has held across time. Genetic differentiation was also lower within northern and southern Ghana cohorts (F ST 0.000–0.005) but higher between these regions (F ST 0.008–0.013), suggesting restricted gene flow between the north and south of the country (Fig. 3 ). Additionally, we observed that An. coluzzii from Greater Accra from 2012 diverged from the main cluster of individuals from southern Ghana on the PCA plot and appeared as a cluster of particularly closely related individuals in the Neighbor-Joining tree. This comparison included southern cohorts collected from the same time point. F ST was also higher when comparing Greater Accra with the other southern population cohorts (F ST 0.005). Although nucleotide diversity values did not differ substantially between any population cohort (Fig. S3 ), the Greater Accra cohort had a slightly higher Tajima’s D indicating it may have experienced greater genetic drift. Overall findings suggest high connectivity among southern populations but restricted gene flow between the different bioclimatic regions of northern and southern Ghana. To investigate further whether population structure in An. coluzzii could be associated with different bioclimatic zones more widely across West Africa, we extended our analysis to include the West African populations of Burkina Faso, Cote d’Ivoire and Mali available through the Vector Observatory ( https://www.malariagen.net/vobs/ ). Once again, we observed two major PCA and NJT clusters (Fig. 2 B and Fig. S4 ). One included northern Ghana and the arid savannah regions of Mali and Burkina Faso. Another included southern Ghana along with southern Cote d’Ivoire, which both experience similar wet humid climate conditions. In support of these findings, F ST was lower between the northern cohorts from Ghana, Burkina Faso and Mali (F ST 0.000-0.001) and between southern Ghana and Cote d’Ivoire (F ST 0.005–0.006) (Fig. 3 ). In contrast, F ST was higher between northern and southern comparisons (F ST 0.008 − 0.0015). These results support the notion that gene flow between An. coluzzii is restricted between populations found in the arid savannah and wet humid forest environments of West Africa. To further investigate gene flow among An. coluzzii in west and central Africa, we investigated two inverted regions of chromosome two previously associated with different climate conditions[ 54 ]. We used tagging SNPs, correlated with inversion status, to assess the frequency of the different inversion karyotypes in each population cohort[ 19 ]. One major 2La inversion karyotype was shared between An. coluzzii from the northern arid regions of Ghana, Burkina Faso and Mali, suggesting unrestricted gene flow among these population cohorts (Fig. 4 A). The inverted 2La inversion prevalent at 96–100% frequency in the northern populations has been associated with arid environments and thermotolerance[ 59 , 60 ] (Table 2 ). In contrast, inverted 2La was between 0–16% frequency in the southern cohorts while the standard 2La + karyotype associated with mesic environments[ 54 ] was dominant at 33–100% frequency. Analysis of the 2Rb inversion (2R:19,444,433 − 26,313,071) presented similar findings in that the southern and northern populations had different dominant inversion karyotypes (Fig. 4 B). However, findings differed in that all three 2Rb karyotypes were present in the northern regions while one form prevailed in the southern populations. Frequency analysis based on inversion tagging SNPs revealed that the dominant southern karyotype was the standard 2Rb + chromosomal form present at 93–100% (Table 2 ). Although all three karyotypes were present in the northern populations, the inverted 2Rb karyotype and heterozygote form reached higher frequencies (11–63% and 32–73%, respectively) than standard 2Rb+ (4–46%) and were more prevalent on average in the north (37% and 45% for 2La and the heterozygote, respectively) than the south (0% and 4%, respectively). Findings are consistent with the observation that inverted 2Rb generally appears at higher frequencies in arid environments while standard 2Rb + is associated with mesic environments[ 54 ] and supports restricted gene flow between the northern Sahelian and southern forest regions of West Africa characterised by different bioclimatic zones. Table 2 Percentage of 2La and 2Rb inversion karyotypes in population cohorts across the bioclimatic zones of West Africa using tagging SNPs for inversion karyotypes[ 19 ]. Country Location Year Inversion karyotype % 2La+ 2La+/2La 2La 2Rb+ 2Rb+/2Rb 2Rb Cote d'Ivoire Lagunes 2012 41 45 14 93 8 0 Ghana Western Region 2012 33 58 8 96 1 0 Central Region 2012 48 36 16 100 3 0 2018 34 50 16 98 3 0 Greater Accra Region 2012 64 29 7 100 0 0 Eastern Region 2012 100 0 0 100 0 0 Ashanti Region 2018 46 45 9 96 10 0 Upper East 2016 0 1 99 12 50 39 2017 0 1 99 8 40 52 2018 0 0 100 7 34 60 Burkina Faso Hauts-Bassins 2012 0 2 98 46 41 12 2014 0 4 96 42 40 19 Mali Koulikoro 2004 0 0 100 9 73 18 2014 0 4 96 4 33 63 Segou 2004 0 4 96 8 32 60 Sikasso 2012 0 4 96 30 59 11 Resistance to Insecticides To investigate geographical differences in insecticide resistance, we focused our analysis on the An. coluzzii data from Ghana in West Africa, which included country-wide sampling and the most recent timepoints for comparison. First, substitutions associated with target-site were investigated. Amino acid allele frequencies of three genes that encode for insecticide binding targets were computed: Vgsc (AGAP004707), Rdl (AGAP006028), and Ace1 (AGAP001356). The kdr allele Vgsc-L995F established as conferring insecticide resistance to pyrethroids[ 27 , 58 ] was identified in all populations but differed between the northern and southern populations with frequencies ranging from 57–64% and 86–92%, respectively (Fig. 5 ). Similar frequencies of Vgsc- L995F within the latter range were observed in the southern cohorts across the different years of collection ranging from 2012 to 2018. Additionally, we observed a double substitution, Vgsc-V402L and I1527T , at 36–43% frequency in northern Ghana and 8–25% in the south, including directly comparable cohorts sampled from the different locations in 2018. The double substitution previously observed in An. coluzzii from Burkina Faso and Kenya[ 42 , 61 ] provides pyrethroid and DDT resistance at a reduced fitness cost compared to L995F [ 62 ]. As a result, the substitution pair is expected to replace L995F as the dominant Vgsc insecticide resistance mechanism. Similar to a study from Burkina Faso[ 61 ], we observed the V402L and I1527T substitution pair increasing in frequency. Frequencies increased 7% from 2016 to 2018 in northern Ghana while the frequency of L995F reduced over the same period, advocating for a shift in the main Vgsc resistance allele in West Africa. We also observed both the allele pairs A296G / T345M and A296S / T345S at the Rdl locus, which are known to confer resistance to organochlorines such as dieldrin[ 63 , 64 ] (Fig. 5 ). However, we found a difference in their geographical distribution. We observed A296G and T345M in all cohorts from southern Ghana only, while A296S and T345S were only found in northern Ghana. This provides support for our finding that gene flow is restricted across Ghana, including at the 2La inversion region on which Rdl is located. To date, the Rdl-A296S / T345S allele pair has only been reported in Burkina Faso[ 61 ] which shares a border with northern Ghana and supports our finding of population connectivity across the northern arid regions of West Africa. To investigate whether the different Rdl substitution pairs were associated with the chromosomal inversion karyotype on which the gene is located, we calculated the frequencies of Rdl substitutions for individuals with either the 2La or 2La + inversion. As expected, all individuals from the Upper East had the inverted 2La inversion karyotype and also the A296S and T345S substitution pair only found in this region (Fig. S5 ). Individuals from southern Ghana with the standard 2La inversion had the A296G and T345M substitution pair associated with the region, with substitution frequencies for 2La population cohorts ranging from 16–58%, but individuals with the inverted 2La inversion did not have either substitution pair. Interestingly, this was the case for individuals collected from the same cohort from Ashanti and the Central Region, suggesting an association of the A296G and T345M substitution with the 2La inversion karyotype. Finally, we observed the G280S mutation and linked CNVs at the Ace1 gene which are associated with resistance to organophosphates and carbamates[ 65 , 66 ] (Fig. 5 ). These were present at > 5% in the population cohort from Greater Accra from 2012 only which formed a somewhat genetically differentiated group upon comparison with the other southern population cohorts on PCA analysis, including other southern populations sampled at the same time (Fig. 2 A). Metabolic Resistance To further investigate metabolic insecticide resistance in An. coluzzii across Ghana, we calculated the frequency of individuals within each population cohort with a copy number greater than two for genes associated with insecticide resistance, including cytochrome P450’s (AGAP002862-AGAP002870, AGAP000818, AGAP008212-AGAP008219), carboxylesterases (AGAP006228, AGAP006723-AGAP006728) and Gste2 (AGAP009194). We also calculated amino acid substitution frequencies for the latter since the Gste2 - I114T mutation is known to increase the activity of Gste2 [ 36 ]. We found CNVs at the cytochrome P450 cluster Cyp 6 on chromosome two (Fig. 6 ). Amplifications at the Cyp6aa1 (AGAP002862) and 2 (AGAP013128) regions were present across Ghana at 60–84% frequency in the more recent population cohorts from 2018. This contrasts with the earlier sampling points from 2012, which presented frequencies between 0–16%, suggesting that metabolic resistance has risen in the country. In addition, we observed duplications at the carboxylesterase cluster Coeae2g-7g , which have been associated with resistance to pirimiphos methyl[ 35 ] at a higher frequency in northern Ghana (12–28%) than southern Ghana (0–6%), including cohorts collected in the same year (Fig. 6 ). Gste2 CNV amplifications and the I114 substitution were also present at higher frequencies in northern Ghana (17–23% and 49–53%, respectively) compared to southern Ghana (4–8% to 35–38%, respectively) (Fig. 6 ). The latter had similar I114 frequencies across the different collection years. The exception was An. coluzzii from Greater Accra in 2012, which presented similar frequencies of both CNVs (21%) and the Gste2 - I114T substitution (46%) to the more recently sampled northern populations. Moreover, we observed CNVs at the cytochrome Cyp9k1 at 21% frequency in Greater Accra, while frequencies were ≤ 5% in all other population cohorts (Fig. 6 ), including southern cohorts collected at the same time. Overall, our findings indicate that metabolic resistance varies across Ghana, with the northern populations in general more impacted than the southern populations, although the region of Greater Accra had particularly high CNV frequencies at insecticide resistance loci. Signals of selection To identify signals of recent selection that may indicate novel insecticide resistance mechanisms, we calculated the H12 homozygosity statistic across windows of the genome and identified statistical peaks[ 67 ]. We observed a number of novel selection signals at loci which have not yet been functionally validated to confer insecticide resistance. We identified a selection peak over the Keap1 locus (AGAP003645: 2R:40,926,195 − 40,945,169) in all cohorts from southern Ghana (Fig. 7 and Fig. S7 ). Keap1 regulates the formation of the transcription factor Maf-S, known to trigger the expression of multiple metabolic resistance genes including cytochrome p450’s and glutathione S-transferases[ 68 ]. A selection signal has recently been observed at this gene in An. arabiensis from Kenya[ 69 ] and associated with deltamethrin treated survivors from Tanzania[ 35 ]. Keap1 therefore provides a good candidate for functional validation. In addition, we observed selection signals unique to the cohort from Greater Accra (Fig. 7 ). We observed a selection signal close to the cytochrome Cyp12f (AGAP008019 3R:4,324,183-4,326,568) which is differentially expressed in permethrin and DDT resistant strains of An. gambiae [ 70 ]. Furthermore, a selection signal was present near a UDP glucuronosyltransferase (AGAP028055 3R:2,836,386-2,838,097) which is a class of uridine diphosphate (UDP)-glycosyltransferase (UGT) detoxification enzymes involved in xenobiotic metabolism[ 71 ]. Findings of selection signals over novel genes with a possible link to insecticide resistance in only Greater Accra in addition to our observation of high frequencies of known target-site and metabolic resistance mechanisms in this cohort, suggests that this region is under particularly high insecticide resistance pressure. The selection signals we observed at known insecticide resistance loci agreed with the CNV and substitution frequencies generated for each population cohort. For example, we observed a selection peak at the Cyp6 gene cluster (2R), Vgsc (2L) and Gste2 (3R) for all population cohorts for which we also observed high frequencies of CNVs and resistance associated substitutions (Fig. S6 -8). A signal at Rdl (2L) was apparent for the populations from the Upper East and Greater Accra, which both had relatively high frequencies of Rdl resistance associated substitutions. Additionally, Greater Accra had a signal at Ace1 (2R) which was the only population cohort with CNVs and the G280S substitution at this locus. Interestingly, there was also a selection peak at Coeae2f and Coeae2g-7g on 2L for Greater Accra only although we only observed CNVs at appreciable frequencies in An. coluzzii from northern Ghana. It could be that the selection signal in Greater Accra is driven by another mechanism i.e., a substitution, although none have been identified as conferring resistance to date. Although functional studies have not yet associated substitutions at Keap1 with insecticide resistance, unique SNPs have been observed in haplotypes under selection for An. arabiensis from Kenya[ 69 ]. These included the substitution D780N and the stop gain mutation E762, which could result in loss of function and prevent the repression of detox gene expression in the absence of stress[ 46 ]. Therefore, we used diplotype clustering[ 34 ] to assess whether CNV duplications or SNPs were associated with haplotypes under selection in An. coluzzii from Ghana. We observed two haplotypes with low heterozygosity in southern Ghana for which we also observed a selection signal at the Keap1 locus (Fig. S9 ). The two haplotypes were associated with a particular set of SNPs, including V631I , V816F and V1001L and G788R and A943V , respectively. In particular, the latter three SNPs are at higher frequencies in southern Ghana (43–52%) than northern Ghana (11–19%) and present at similar frequencies, suggesting they are either linked or have a synergistic effect (Fig. S10 ). These SNPs are different from those previously observed in East African An. arabiensis , but different molecular responses to insecticides are commonly observed across different taxa[ 61 , 67 ]. The role of the variants we observed in insecticide resistance could be determined through validation studies and also provide a focus for longitudinal studies tracking the frequency and geographical spread of candidate variants. We also used diplotype clustering to further investigate genomic variation at the Carboxylesterase cluster Coeae2g-7g since we observed a selection signal for An. coluzzii from Greater Accra despite the absence of CNVs. However, the dendrogram did not reveal any low heterozygosity clusters expected to result from selection on the region (Fig. S11 ). Increased sampling of the Greater Accra Region may improve resolution to detect selected variants. Discussion Using whole genome sequence variation data from West African An. coluzzii , we observed population structure and divergent insecticide resistance mechanisms driven by gene flow, environmental variation, and differences in selection pressure from population control. We found two genetically different groups of An. coluzzii based on population structure analysis of a large number of SNPs and two chromosomal inversion regions, both with similar levels of nucleotide diversity. These two populations may represent An. coluzzii subject to distinct ecological selection pressures in the arid savannah region of northern West Africa and the wet deciduous forest region of the south. Isolation by distance could explain the pattern of population structure we have observed. However, since our sites in the south were ~ 500 km from those in the north and we did not have samples from intermediate sites, we were unable to determine whether this pattern reflects mosquito dispersal. However, the pattern of population structure between An. coluzzii from the different bioclimatic zones was present across West Africa where An. coluzzii are also expected to undergo long-range windborne dispersal to find suitable oviposition sites during the dry season[ 72 ]. Furthermore, there are no major geographical barriers expected to prevent mosquito dispersal across the investigated region. Therefore, another explanation for the population structure we observed is that reduced habitat suitability or local environmental adaptation restricts gene flow between northern (i.e Sahelian) and southern (ie. forest) West Africa. The former is supported by the geographical distribution of An. coluzzii which is bimodal because they predominate along the savannah and coastal regions[ 73 ], including in Ghana where aridity is increased in the north and southernmost coastal regions compared to the central zone[ 74 ]. However, it could also be that environmental adaptation leads to population specific mating or a fitness cost that restricts gene flow. Genomic evidence for environmental adaptation in An. coluzzii has been previously described based on the segregation of chromosome two inversions across West and Central Africa[ 59 ]. Similarly, we found that the inversion karyotypes of 2La and 2Rb segregated among An. coluzzii from northern and southern West Africa. Our findings support previous observations from Nigeria, Burkina Faso and Cameroon[ 12 , 13 , 15 , 17 ] that the inverted 2La and 2Rb predominates in arid regions while standard karyotypes are common in humid forest environments. Such inversions have the potential to preserve locally selected alleles by reducing recombination between heterozygotes and are thought to be major mechanisms driving speciation and the formation of Anopheles ecotypes[ 3 , 75 ]. Ecological studies comparing An. coluzzii across West Africa are lacking and therefore further work is required to determine whether distinct phenotypic and/or behavioral differences are observed between cohorts from the different bioclimatic zones. However, a preference for mating with the local population has been observed for An. coluzzii sourced from different ecozones in Burkina Faso under experimental conditions, which are analogous to those in other West African countries[ 9 ]. Distinct phenotypic differences have also been observed in summer dormancy, i.e., aestivation, from contrasting arid and mesic environments[ 11 , 76 ]. Yet the question remains over whether aestivation in An. coluzzii has a genomic basis and can be linked to habitat use or geographical origin or whether it is a plastic response to dry season conditions, i.e., it is induced by the environment similar to the diapause response in other mosquitoes[ 76 , 77 ]. It is also unknown to what extent aestivating An. coluzzii can increase population structure between cohorts. Overall, further experimental work linking ecologically driven phenotypes to genomic variation is required to untangle whether An. coluzzi i from contrasting environments in West Africa represent ecotypes adapted to different environmental conditions, but our findings raise this possibility. The presence of ecologically distinct forms of An. coluzzii would have consequences for population control and malaria transmission since they will exhibit different life history traits and different responses to insecticide use. Indeed, we found that restricted gene flow between An. coluzzii across the different bioclimatic zones of northern and southern Ghana has led to different insecticide resistance profiles. Findings support previous work which found F ST outlier regions at target site and metabolic resistance genes on comparison of forest and non-forest populations in Ghana, although comparisons did not include sites from northern Ghana[ 43 ]. Overall, we found that An. coluzzii from northern Ghana exhibited greater evidence for both target site and metabolic resistance than the southern populations, excepting the Greater Accra region. Along with higher frequencies of the known Vgsc -L995F substitution associated with pyrethroid resistance, northern Ghana also had the double substitution Vgsc -V402L and I1527T previously reported from nearby connected Burkina Faso[ 61 ]. Since this substitution pair provides pyrethroid resistance at a lower fitness cost than L995F[ 62 ], we expect it to rapidly spread across northern West Africa where we observed unrestricted gene flow. In addition, we observed higher frequencies of carboxylesterase CNVs and GSTE substitutions and CNVs in northern Ghana compared to the south, associated with multiple classes of insecticides including pirimiphos-methyl[ 40 ]. Although bed net use is widespread across the country, northern Ghana has a high malaria burden and has been additionally targeted with high coverage of long-lasting insecticide treated nets and IRS campaigns with propoxur between 2012 and 2016 [ 78 , 79 ]. Despite a lower prevalence of molecular insecticide resistance, An. coluzzii from southern Ghana are still impacted mainly through target site resistance at Vgsc -L955F. In addition, we found that the populations in southern Ghana were uniquely impacted by a selection signal at Keap1 . This gene is also under selection in East Africa and has the potential to impact the expression of detoxification genes to metabolise insecticides[ 46 , 69 ]. For example, a reduction in the regulatory action of Keap1 decreases mortality to organophosphates but increases mortality to pyrethroids and DDT[ 68 ]. Although it is unknown whether substitutions at the gene impact on insecticide resistance, we have identified several SNPs including V631I, V816F and V1001L and G788R and A943V which can be targeted for experimental work, including testing whether these variants have a synergistic effect. Although we also observed differences in Rdl substitutions between different bioclimatic zones and inversion karyotypes, the source of selection on this locus is unclear since dieldrin is no longer used in vector control. Possibly dieldrin use is ongoing, i.e., in agriculture, or the locus provides cross-resistance to another used insecticide. We found that An. coluzzii from Greater Accra on the southern coast of Ghana were somewhat genetically divergent from other southern populations based on PCA and F ST population differentiation statistics. In addition, they had unique molecular insecticide resistance mechanisms and novel selection signals on chromosome three. For example, we found appreciable frequencies of Cyp9k1 CNVs, the Ace1 -G280S substitution and Ace1 CNVs at this locus in Greater Accra only, suggesting resistance to multiple classes of insecticides. In addition, we found novel signals of selection at the cytochrome Cyp12f linked to permethrin and DDT resistance in An. gambiae [ 70 ] and a UDP glucuronosyltransferase analogous to a UGT gene that confers resistance to the pyrethroid lambda-cyhalthrin and the neonicotinoid imidacloprid in fruitflys[ 71 ]. Findings suggest that An. coluzzii are under intense selection pressure from insecticides in Greater Accra. This notion is supported by bioassay data, which suggested particularly high resistance to pyrethroids and carbamates in An. gambiae from the region[ 80 ]. However, the findings of phenotypic resistance are not directly comparable with the genomic data since our sampling was from 2012 while the bioassay experiment was conducted in 2017. The findings are somewhat surprising given that malaria incidence and bed net use in Greater Accra is comparatively low to elsewhere in Ghana[ 81 , 82 ]. However, Greater Accra is a populated urban area where there could be substantial use of household insecticide sprays for personal protection[ 83 ]. Furthermore, mosquito populations encounter increased pollution in urban environments, including habitats contaminated with insecticides from urban agriculture[ 73 , 80 ]. For example, pollution from urban agriculture can result in similar levels of resistance when compared to populations from rural and cultivated areas[ 84 – 86 ]. Further study is required to identify the drivers of selection pressure within the city of Greater Accra. Even so, it is concerning that we have observed multiple and unique resistance mechanisms likely to reduce the efficacy of malaria control. Since we have found unrestricted gene flow across southern Ghana, any novel insecticide resistance mechanism has the potential to spread. However, this will also depend on selection pressure across the country, since resistance is often accompanied by a fitness cost[ 87 , 88 ]. Experimental work through bioassays and/or association studies are required to confirm that the novel regions under selection confer a resistance mechanism and to pinpoint the genomic variation underlying the phenotype. Our results demonstrate that gene flow among vector populations is important in influencing the distribution of insecticide resistance mechanisms. A full understanding of mosquito population structure, both regional and large-scale, is required to predict the success of gene drive technologies and how novel insecticide resistance mechanisms will spread when they arise in response to population control measures. This includes a greater understanding of how environmental conditions influence mosquito dispersal and connectivity, and how this may alter with climate change. A strong understanding of population connectivity is particularly important given the introduction of new technologies introduced to combat the rise in metabolic resistance, i.e., dual active ingredient (AI) nets. New technologies impose novel selection pressures and may be quickly challenged given Anopheles propensity for a rapid evolutionary response[ 61 , 89 , 90 ]. Any novel insecticide resistance mechanism will need to be quickly managed to maintain efficacy on their introduction. In tandem, large-scale routine monitoring of the temporal and geographical distribution of molecular insecticide resistance mechanisms across West Africa will be essential for a targeted defense. Declarations Funding The MalariaGEN Vector Observatory is supported by multiple institutes and funders. The Wellcome Sanger Institute’s participation was supported by funding from Wellcome (220540/Z/20/A, 'Wellcome Sanger Institute Quinquennial Review 2021-2026') and the Bill & Melinda Gates Foundation (INV-001927 and INV-068808). The Liverpool School of Tropical Medicine's participation was supported by the National Institute of Allergy and Infectious Diseases ([NIAID] R01-AI116811), with additional support from the Medical Research Council (MR/P02520X/1). The latter grant is a UK-funded award and is part of the EDCTP2 programme supported by the European Union. Martin Donnelly is supported by a Royal Society Wolfson Fellowship (RSWF\FT\180003). The Pan-African Mosquito Control Association’s participation was funded by the Bill and Melinda Gates Foundation (INV-031595). Lucas N. Amenga-Etego is supported by the Bill & Melinda Gates Foundation (INV-050873) and the National Institute of Health and Care Research, UK (grant number: NIHR134717). Author contributions AM, CC, GA and LAE conceptualised and designed the study, interpreted the data and assisted in drafting the manuscript. EKA, ID, CMM, SB, VAA, CA, KLM, designed the study and conducted sample collection, processing, and data collection. EKA, AHK and KLB conducted the data analysis and interpretation and assisted in drafting the manuscript. All the authors read, reviewed and approved this manuscript. Competing interests The authors declare that they have no competing interests. No chatbots or artificial intelligence tools were used in any of these studies. Materials and correspondence Enock Kofi Amoako ; Lucas Amenga Etego; Data availability The sequences of the samples identified in this study were submitted to the European Nucleotide Archive (ENA) (Project: PRJEB2141, accessions ERR2656751-ERR9796298). Acknowledgements This study was supported by the MalariaGEN Vector Observatory which is an international collaboration working to build capacity for malaria vector genomic research and surveillance, and involves contributions by the following institutions and teams. Wellcome Sanger Institute: Lee Hart, Kelly L. Bennett, Anastasia Hernandez-Koutoucheva, Jon Brenas, Menelaos Ioannidis, Chris Clarkson, Alistair Miles, Julia Jeans, Paballo Chauke, Victoria Simpson, Eleanor Drury, Osama Mayet, Sónia Gonçalves, Katherine Figueroa, Tom Maddison, Kevin Howe, Mara Lawniczak; Liverpool School of Tropical Medicine: Eric Lucas, Sanjay Nagi, Martin Donnelly; Broad Institute of Harvard and MIT: Jessica Way, George Grant; Pan-African Mosquito Control Association: Jane Mwangi, Edward Lukyamuzi, Sonia Barasa, Ibra Lujumba, Elijah Juma. The authors would like to thank the staff of the Wellcome Sanger Genomic Surveillance unit and the Wellcome Sanger Institute Sample Logistics, Sequencing and Informatics facilities for their contributions. The MalariaGEN Vector Observatory is supported by multiple institutes and funders. The Wellcome Sanger Institute’s participation was supported by funding from Wellcome (220540/Z/20/A, 'Wellcome Sanger Institute Quinquennial Review 2021-2026') and the Bill & Melinda Gates Foundation (INV-001927 and INV-068808). The Liverpool School of Tropical Medicine's participation was supported by the National Institute of Allergy and Infectious Diseases ([NIAID] R01-AI116811), with additional support from the Medical Research Council (MR/P02520X/1). The latter grant is a UK-funded award and is part of the EDCTP2 programme supported by the European Union. Martin Donnelly is supported by a Royal Society Wolfson Fellowship (RSWF\FT\180003). The Pan-African Mosquito Control Association’s participation was funded by the Bill and Melinda Gates Foundation (INV-031595). References Mwima R, Hui T-YJ, Nanteza A, Burt A, Kayondo JK. Potential persistence mechanisms of the major Anopheles gambiae species complex malaria vectors in sub-Saharan Africa: a narrative review. Malar J. 2023;22:336. De Meeûs T, Bouyer J, Ravel S, Solano P. Ecotype Evolution in Glossina palpalis Subspecies, Major Vectors of Sleeping Sickness. PLoS Negl Trop Dis. 2015;9:e0003497-. Small ST, Costantini C, Sagnon N, Guelbeogo MW, Emrich SJ, Kern AD, et al. Standing genetic variation and chromosome differences drove rapid ecotype formation in a major malaria mosquito. Proceedings of the National Academy of Sciences. 2023;120:e2219835120. 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:207–11. de Souza D, Kelly-Hope L, Lawson B, Wilson M, Boakye D. Environmental Factors Associated with the Distribution of Anopheles gambiae s.s in Ghana; an Important Vector of Lymphatic Filariasis and Malaria. PLoS One. 2010;5:e9927. Donkor E, Kelly M, Eliason C, Amotoh C, Gray DJ, Clements ACA, et al. A Bayesian spatio-temporal analysis of malaria in the Greater Accra region of Ghana from 2015 to 2019. Int J Environ Res Public Health. 2021;18:6080. Hinne IA, Attah SK, Mensah BA, Forson AO, Afrane YA. Larval habitat diversity and Anopheles mosquito species distribution in different ecological zones in Ghana. Parasit Vectors. 2021;14:193. Ghana Statistical Service. Ghana Multiple Indicator Cluster Survey with an Enhanced Malaria Module and Biomarker. 2011. Nignan C, Poda BS, Sawadogo SP, Maïga H, Dabiré KR, Gnankine O, et al. Local adaptation and colonization are potential factors affecting sexual competitiveness and mating choice in Anopheles coluzzii populations. Sci Rep. 2022;12:636. Lehmann T, Weetman D, Huestis DL, Yaro AS, Kassogue Y, Diallo M, et al. Tracing the origin of the early wet-season Anopheles coluzzii in the Sahel. Evol Appl. 2017;10:704–17. Yaro AS, Traoré AI, Huestis DL, Adamou A, Timbiné S, Kassogué Y, et al. Dry season reproductive depression of Anopheles gambiae in the Sahel. J Insect Physiol. 2012;58:1050–9. Coluzzi M, Sabatini A, Petrarca V, Di Deco MA. Chromosomal differentiation and adaptation to human environments in the Anopheles gambiae complex. Trans R Soc Trop Med Hyg. 1979;73:483–97. Costantini C, Ayala D, Guelbeogo WM, Pombi M, Some CY, Bassole IHN, et al. Living at the edge: biogeographic patterns of habitat segregation conform to speciation by niche expansion in Anopheles gambiae. BMC Ecol. 2009;9:16. Simard F, Ayala D, Kamdem GC, Pombi M, Etouna J, Ose K, et al. Ecological niche partitioning between Anopheles gambiae molecular forms in Cameroon: the ecological side of speciation. BMC Ecol. 2009;9:17. Adeogun AO, Popoola KOK, Brooke BD, Olakiigbe AK, Awolola ST. Polymorphic inversion 2La frequencies associated with ecotypes in populations of Anopheles coluzzii from Southwest Nigeria. Sci Afr. 2021;12:e00746. Ayala D, Acevedo P, Pombi M, Dia I, Boccolini D, Costantini C, et al. Chromosome inversions and ecological plasticity in the main African malaria mosquitoes. Evolution (N Y). 2017;71:686–701. Simard F, Ayala D, Kamdem GC, Pombi M, Etouna J, Ose K, et al. Ecological niche partitioning between Anopheles gambiae molecular forms in Cameroon: the ecological side of speciation. BMC Ecol. 2009;9:1–24. Coluzzi M, Sabatini A, Petrarca V, Di Deco MA. Chromosomal differentiation and adaptation to human environments in the Anopheles gambiae complex. Trans R Soc Trop Med Hyg. 1979;73:483–97. Love RR, Pombi M, Guelbeogo MW, Campbell NR, Stephens MT, Dabire RK, et al. Inversion Genotyping in the Anopheles gambiae Complex Using High-Throughput Array and Sequencing Platforms. G3 Genes|Genomes|Genetics. 2020;10:3299–307. Petrarca V, Beier JC. Intraspecific chromosomal polymorphism in the Anopheles gambiae complex as a factor affecting malaria transmission in the Kisumu area of Kenya. Am J Trop Med Hyg. 1992;46:229–37. Adeogun AO, Brooke BD, Olayanju DR, Adegbehingbe K, Oyeniyi TA, Olakiigbe AK, et al. Test for association between dieldrin resistance and 2La inversion polymorphism in Anopheles coluzzii from Lagos, Nigeria. 2019. Tropical biomedicine, 36;3:587–593. Brooke BD, Hunt RH, Coetzee M. Resistance to dieldrin + fipronil assorts with chromosome inversion 2La in the malaria vector Anopheles gambiae. Med Vet Entomol. 2000;14:190–4. Clarkson CS, Weetman D, Essandoh J, Yawson AE, Maslen G, Manske M, et al. Adaptive introgression between Anopheles sibling species eliminates a major genomic island but not reproductive isolation. Nat Commun. 2014;5:4248. 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–96. Hemingway J, Ranson H, Magill A, Kolaczinski J, Fornadel C, Gimnig J, et al. Averting a malaria disaster: will insecticide resistance derail malaria control? The Lancet. 2016;387:1785–8. Davies TGE, Field LM, Usherwood PNR, Williamson MS. A comparative study of voltage-gated sodium channels in the Insecta: implications for pyrethroid resistance in Anopheline and other Neopteran species. Insect Mol Biol. 2007;16:361–75. Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Bergé JB, Devonshire AL, et al. Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998;7:179–84. Essandoh J, Yawson AE, Weetman D. Acetylcholinesterase (Ace-1) target site mutation 119S is strongly diagnostic of carbamate and organophosphate resistance in Anopheles gambiae ss and Anopheles coluzzii across southern Ghana. Malar J. 2013;12:404. Weill M, Malcolm C, Chandre F, Mogensen K, Berthomieu A, Marquine M, et al. The unique mutation in ace‐1 giving high insecticide resistance is easily detectable in mosquito vectors. Insect Mol Biol. 2004;13:1–7. Hancock PA, Ochomo E, Messenger LA. Genetic surveillance of insecticide resistance in African Anopheles populations to inform malaria vector control. Trends Parasitol. 2024;40:604–18. Du W, Awolola TS, Howell P, Koekemoer LL, Brooke BD, Benedict MQ, et al. Independent mutations in the Rdl locus confer dieldrin resistance to Anopheles gambiae and An. arabiensis. Insect Mol Biol. 2005;14:179–83. Edi C V, Djogbénou L, Jenkins AM, Regna K, Muskavitch MAT, Poupardin R, et al. CYP6 P450 Enzymes and ACE-1 Duplication Produce Extreme and Multiple Insecticide Resistance in the Malaria Mosquito Anopheles gambiae. PLoS Genet. 2014;10:e1004236. Ibrahim SS, Riveron JM, Stott R, Irving H, Wondji CS. The cytochrome P450 CYP6P4 is responsible for the high pyrethroid resistance in knockdown resistance-free Anopheles arabiensis. Insect Biochem Mol Biol. 2016;68:23–32. Nagi SC, Lucas ER, Egyir-Yawson A, Essandoh J, Dadzie S, Chabi J, et al. Parallel Evolution in Mosquito Vectors—A Duplicated Esterase Locus is Associated With Resistance to Pirimiphos-methyl in Anopheles gambiae. Mol Biol Evol. 2024;41:msae140. Lucas ER, Nagi SC, Kabula B, Batengana B, Kisinza W, Egyir-Yawson A, et al. Copy number variants underlie major selective sweeps in insecticide resistance genes in Anopheles arabiensis. PLoS Biol. 2024;22:e3002898. Riveron JM, Yunta C, Ibrahim SS, Djouaka R, Irving H, Menze BD, et al. A single mutation in the GSTe2 gene allows tracking of metabolically based insecticide resistance in a major malaria vector. Genome Biol. 2014;15:R27. Pu J, Chung H. New and emerging mechanisms of insecticide resistance. Curr Opin Insect Sci. 2024;63:101184. Dabiré RK, Namountougou M, Sawadogo SP, Yaro LB, Toé HK, Ouari A, et al. Population dynamics of Anopheles gambiae s.l. in Bobo-Dioulasso city: bionomics, infection rate and susceptibility to insecticides. Parasit Vectors. 2012;5:127. Perugini E, Pichler V, Guelbeogo WM, Micocci M, Poggi C, Manzi S, et al. Longitudinal survey of insecticide resistance in a village of central region of Burkina Faso reveals co-occurrence of 1014F, 1014S and 402L mutations in Anopheles coluzzii and Anopheles arabiensis. Malar J. 2024;23:250. Lucas ER, Nagi SC, Egyir-Yawson A, Essandoh J, Dadzie S, Chabi J, et al. Genome-wide association studies reveal novel loci associated with pyrethroid and organophosphate resistance in Anopheles gambiae and Anopheles coluzzii. Nat Commun. 2023;14:4946. Ibrahim SS, Muhammad A, Hearn J, Weedall GD, Nagi SC, Mukhtar MM, et al. Molecular drivers of insecticide resistance in the Sahelo-Sudanian populations of a major malaria vector Anopheles coluzzii. BMC Biol. 2023;21:125. Kamau L, Bennett KL, Ochomo E, Herren J, Agumba S, Otieno S, et al. The Anopheles coluzzii range extends into Kenya: detection, insecticide resistance profiles and population genetic structure in relation to conspecific populations in West and Central Africa. Malar J. 2024;23:122. Dennis TPW, Essandoh J, Mable BK, Viana MS, Yawson AE, Weetman David. Signatures of adaptation at key insecticide resistance loci in Anopheles gambiae in Southern Ghana revealed by reduced-coverage WGS. Sci Rep. 2024;14:8650. Consortium A gambiae 1000G. Genetic diversity of the African malaria vector Anopheles gambiae. Nature. 2017;552:96. Clarkson CS, Miles A, Harding NJ, Lucas ER, Battey CJ, Amaya-Romero JE, et al. Genome variation and population structure among 1142 mosquitoes of the African malaria vector species Anopheles gambiae and Anopheles coluzzii. Genome Res. 2020;30:1533–46. Ingham VA, Tennessen JA, Lucas ER, Elg S, Yates HC, Carson J, et al. Integration of whole genome sequencing and transcriptomics reveals a complex picture of the reestablishment of insecticide resistance in the major malaria vector Anopheles coluzzii. PLoS Genet. 2021;17:e1009970. Coetzee M. Key to the females of Afrotropical Anopheles mosquitoes (Diptera: Culicidae). Malar J. 2020;19:70. Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25:1754–60. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297–303. Martin M, Ebert P, Marschall T. Read-Based Phasing and Analysis of Phased Variants with WhatsHap. In: Peters BA, Drmanac R, editors. Haplotyping: Methods and Protocols. New York, NY: Springer US; 2023. p. 127–38. Delaneau O, Marchini J, McVean GA, Donnelly P, Lunter G, Marchini JL, et al. Integrating sequence and array data to create an improved 1000 Genomes Project haplotype reference panel. Nat Commun. 2014;5:3934. Lucas ER, Miles A, Harding NJ, Clarkson CS, Lawniczak MKN, Kwiatkowski DP, et al. Whole-genome sequencing reveals high complexity of copy number variation at insecticide resistance loci in malaria mosquitoes. Genome Res. 2019;29:1250–61. Neafsey DE, Waterhouse RM, Abai MR, Aganezov SS, Alekseyev MA, Allen JE, et al. Highly evolvable malaria vectors: the genomes of 16 Anopheles mosquitoes. Science (1979). 2015;347:1258522. Ayala D, Ullastres A, González J. Adaptation through chromosomal inversions in Anopheles. Front Genet. 2014;5:129. Love RR, Pombi M, Guelbeogo MW, Campbell NR, Stephens MT, Dabire RK, et al. Inversion Genotyping in the Anopheles gambiae Complex Using High-Throughput Array and Sequencing Platforms. G3 Genes|Genomes|Genetics. 2020;10:3299–307. Hudson RR, Slatkin M, Maddison WP. Estimation of levels of gene flow from DNA sequence data. Genetics. 1992;132:583–9. Garud NR, Messer PW, Buzbas EO, Petrov DA. Recent Selective Sweeps in North American Drosophila melanogaster Show Signatures of Soft Sweeps. PLoS Genet. 2015;11:e1005004. Clarkson CS, Miles A, Harding NJ, O’Reilly AO, Weetman D, Kwiatkowski D, et al. The genetic architecture of target‐site resistance to pyrethroid insecticides in the African malaria vectors Anopheles gambiae and Anopheles coluzzii. Mol Ecol. 2021;30:5303–17. Ayala D, Zhang S, Chateau M, Fouet C, Morlais I, Costantini C, et al. Association mapping desiccation resistance within chromosomal inversions in the African malaria vector Anopheles gambiae. Mol Ecol. 2019;28:1333–42. Ibrahim SS, Mukhtar MM, Muhammad A, Wondji CS. 2La paracentric chromosomal inversion and overexpressed metabolic genes enhance thermotolerance and pyrethroid resistance in the major malaria vector Anopheles gambiae. Biology. 2021;10:518. Kientega M, Clarkson CS, Traoré N, Hui T-YJ, O’Loughlin S, Millogo A, et al. Whole-genome sequencing of major malaria vectors reveals the evolution of new insecticide resistance variants in a longitudinal study in Burkina Faso. Malar J. 2023;:280. Williams J, Cowlishaw R, Sanou A, Ranson H, Grigoraki L. In vivo functional validation of the V402L voltage gated sodium channel mutation in the malaria vector An. gambiae. Pest Manag Sci. 2022;78:1155–63. Du W, Awolola TS, Howell P, Koekemoer LL, Brooke BD, Benedict MQ, et al. Independent mutations in the Rdl locus confer dieldrin resistance to Anopheles gambiae and An. arabiensis. Insect Mol Biol. 2005;14:179–83. Wondji CS, Dabire RK, Tukur Z, Irving H, Djouaka R, Morgan JC. Identification and distribution of a GABA receptor mutation conferring dieldrin resistance in the malaria vector Anopheles funestus in Africa. Insect Biochem Mol Biol. 2011;41:484–91. Lucas ER, Rockett KA, Lynd A, Essandoh J, Grisales N, Kemei B, et al. A high throughput multi-locus insecticide resistance marker panel for tracking resistance emergence and spread in Anopheles gambiae. Sci Rep. 2019;9:13335. Mitchell SN, Rigden DJ, Dowd AJ, Lu F, Wilding CS, Weetman D, et al. Metabolic and target-site mechanisms combine to confer strong DDT resistance in Anopheles gambiae. PLoS One. 2014;9:e92662. Mwinyi SH, Bennett KL, Nagi SC, Kabula B, Matowo J, Weetman D, et al. Genomic Analysis Reveals a New Cryptic Taxon Within the Anopheles gambiae Complex With a Distinct Insecticide Resistance Profile in the Coast of East Africa. Mol Ecol. 2025;:e17762. Ingham VA, Pignatelli P, Moore JD, Wagstaff S, Ranson H. The transcription factor Maf-S regulates metabolic resistance to insecticides in the malaria vector Anopheles gambiae. BMC Genomics. 2017;18:1–11. Polo B, Bennett KL, Barasa S, Brenas J, Agumba S, Mwangangi J, et al. Genomic surveillance reveals geographical heterogeneity and differences in known and novel insecticide resistance mechanisms in Anopheles arabiensis across Kenya. 2024. David J-P, Strode C, Vontas J, Nikou D, Vaughan A, Pignatelli PM, et al. The Anopheles gambiae detoxification chip: a highly specific microarray to study metabolic-based insecticide resistance in malaria vectors. Proceedings of the National Academy of Sciences. 2005;102:4080–4. Logan RAE, Mäurer JB, Wapler C, Ingham VA. Uridine diphosphate (UDP)-glycosyltransferases (UGTs) are associated with insecticide resistance in the major malaria vectors Anopheles gambiae s.l. and Anopheles funestus. Sci Rep. 2024;14:19821. Huestis DL, Dao A, Diallo M, Sanogo ZL, Samake D, Yaro AS, et al. Windborne long-distance migration of malaria mosquitoes in the Sahel. Nature. 2019;574:404–8. Tene‐Fossog B, Fotso‐Toguem YG, Amvongo‐Adjia N, Ranson H, Wondji CS. Temporal variation of high‐level pyrethroid resistance in the major malaria vector Anopheles gambiae sl in Yaoundé, Cameroon, is mediated by target‐site and metabolic resistance. Med Vet Entomol. 2022;36:247–59. Kudom AA. Larval ecology of Anopheles coluzzii in Cape Coast, Ghana: water quality, nature of habitat and implication for larval control. Malar J. 2015;14:447. Ayala FJ, Coluzzi M. Chromosome speciation: Humans, Drosophila, and mosquitoes. Proceedings of the National Academy of Sciences. 2005;102 suppl_1:6535–42. Hidalgo K, Siaussat D, Braman V, Dabiré KR, Simard F, Mouline K, et al. Comparative physiological plasticity to desiccation in distinct populations of the malarial mosquito Anopheles coluzzii. Parasit Vectors. 2016;9:1–13. Armbruster PA. Photoperiodic diapause and the establishment of Aedes albopictus (Diptera: Culicidae) in North America. J Med Entomol. 2016;53:1013–23. Gogue C, Wagman J, Tynuv K, Saibu A, Yihdego Y, Malm K, et al. An observational analysis of the impact of indoor residual spraying in Northern, Upper East, and Upper West Regions of Ghana: 2014 through 2017. Malar J. 2020;19:1–13. Tiedje KE, Oduro AR, Bangre O, Amenga-Etego L, Dadzie SK, Appawu MA, et al. Indoor residual spraying with a non-pyrethroid insecticide reduces the reservoir of Plasmodium falciparum in a high-transmission area in northern Ghana. PLOS Global Public Health. 2022;2:e0000285. Pwalia R, Joannides J, Iddrisu A, Addae C, Acquah-Baidoo D, Obuobi D, et al. High insecticide resistance intensity of Anopheles gambiae (sl) and low efficacy of pyrethroid LLINs in Accra, Ghana. Parasit Vectors. 2019;12:1–9. Kawaguchi K, Donkor E, Lal A, Kelly M, Wangdi K. Distribution and risk factors of malaria in the Greater Accra Region in Ghana. Int J Environ Res Public Health. 2022;19:12006. Aheto JMK, Menezes LJ, Takramah W, Cui L. Modelling spatiotemporal variation in under-five malaria risk in Ghana in 2016–2021. Malar J. 2024;23:102. Silva Martins WF, Reid E, Tomlinson S, Evans G, Gibson J, Guy A, et al. Improving the efficiency of aerosolized insecticide testing against mosquitoes. Sci Rep. 2023;13:6281. Chabi J, Eziefule MC, Pwalia R, Joannides J, Obuobi D, Amlalo G, et al. Impact of urban agriculture on the species distribution and insecticide resistance profile of Anopheles gambiae ss and Anopheles coluzzii in Accra Metropolis, Ghana. Advances in Entomology. 2018;6:198. Tchigossou G, Dossou C, Tepa-Yotto G, Koto M, Atoyebi SM, Tossou E, et al. Resistance to neonicotinoids is associated with metabolic detoxification mechanisms in Anopheles coluzzii from agricultural and urban sites in southern Benin. Frontiers in Tropical Diseases. 2024;5:1339811. Antonio-Nkondjio C, Fossog BT, Ndo C, Djantio BM, Togouet SZ, Awono-Ambene P, et al. 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. Nkahe DL, Kopya E, Djiappi-Tchamen B, Toussile W, Sonhafouo-Chiana N, Kekeunou S, et al. Fitness cost of insecticide resistance on the life-traits of a Anopheles coluzzii population from the city of Yaoundé, Cameroon. Wellcome Open Res. 2020;5:171. Gul H, Gadratagi BG, Güncan A, Tyagi S, Ullah F, Desneux N, et al. Fitness costs of resistance to insecticides in insects. Front Physiol. 2023;14:1238111. Toé KH, N’Falé S, Dabiré RK, Ranson H, Jones CM. The recent escalation in strength of pyrethroid resistance in Anopheles coluzzi in West Africa is linked to increased expression of multiple gene families. BMC Genomics. 2015;16:146. Njoroge H, van’t Hof A, Oruni A, Pipini D, Nagi SC, Lynd A, et al. Identification of a rapidly-spreading triple mutant for high-level metabolic insecticide resistance in Anopheles gambiae provides a real-time molecular diagnostic for antimalarial intervention deployment. Mol Ecol. 2022;31:4307–18. Additional Declarations No competing interests reported. Supplementary Files Table1samplecollection.csv Table22La2Rbinversionkaryotypes.csv TableS1samplingtable.csv FigS5SupplementaryfigureRdlinversion.tif FigS9keap1diplotypeclustering.tif FigS10keapaa.tif FigS11Supplementaryfigurecoeaediplotypeclustering.tif FigS2Ghananjt.tif FigS3diversity.tif FigS7SignalsChr3.tif FigS1WestAfricanCountriesMap.tif FigS4Westafricanjt.tif FigS6SignalsChr2.tif FigS8signalsChrX.tif Cite Share Download PDF Status: Published Journal Publication published 07 Jan, 2026 Read the published version in BMC Genomics → Version 1 posted Editorial decision: Revision requested 10 Nov, 2025 Reviews received at journal 03 Nov, 2025 Reviews received at journal 30 Oct, 2025 Reviewers agreed at journal 28 Oct, 2025 Reviewers agreed at journal 24 Oct, 2025 Reviewers agreed at journal 24 Oct, 2025 Reviewers invited by journal 23 Oct, 2025 Editor invited by journal 22 Oct, 2025 Editor assigned by journal 20 Oct, 2025 Submission checks completed at journal 20 Oct, 2025 First submitted to journal 16 Oct, 2025 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-7878288","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":539396089,"identity":"5ab154c7-e5c3-42ac-bac2-1d8a6ce3797f","order_by":0,"name":"Enock K. Amoako","email":"","orcid":"","institution":"West African Centre for Cell Biology of Infectious Pathogens (WACCBIP), University of Ghana","correspondingAuthor":false,"prefix":"","firstName":"Enock","middleName":"K.","lastName":"Amoako","suffix":""},{"id":539396090,"identity":"5dccf8c0-22d4-4bba-93a5-ca08425ff0ba","order_by":1,"name":"Kelly L. Bennett","email":"","orcid":"","institution":"Genomic Surveillance Unit, Wellcome Sanger Institute, Hinxton","correspondingAuthor":false,"prefix":"","firstName":"Kelly","middleName":"L.","lastName":"Bennett","suffix":""},{"id":539396091,"identity":"ba4336bc-b9e6-4596-8f49-853342b94c5d","order_by":2,"name":"Anastasia Hernandez-Koutoucheva","email":"","orcid":"","institution":"Genomic Surveillance Unit, Wellcome Sanger Institute, Hinxton","correspondingAuthor":false,"prefix":"","firstName":"Anastasia","middleName":"","lastName":"Hernandez-Koutoucheva","suffix":""},{"id":539396092,"identity":"de58dd1d-1076-43e1-b811-e8a3e0201170","order_by":3,"name":"Isaiah Debrah","email":"","orcid":"","institution":"West African Centre for Cell Biology of Infectious Pathogens (WACCBIP), University of Ghana","correspondingAuthor":false,"prefix":"","firstName":"Isaiah","middleName":"","lastName":"Debrah","suffix":""},{"id":539396093,"identity":"715d0238-963c-4a80-8281-fc8ba636ffc9","order_by":4,"name":"Collins M. Morang’a","email":"","orcid":"","institution":"West African Centre for Cell Biology of Infectious Pathogens (WACCBIP), University of Ghana","correspondingAuthor":false,"prefix":"","firstName":"Collins","middleName":"M.","lastName":"Morang’a","suffix":""},{"id":539396094,"identity":"850a544e-e366-42f9-9bbe-cf77f73730d1","order_by":5,"name":"Stephen Binaansim","email":"","orcid":"","institution":"West African Centre for Cell Biology of Infectious Pathogens (WACCBIP), University of Ghana","correspondingAuthor":false,"prefix":"","firstName":"Stephen","middleName":"","lastName":"Binaansim","suffix":""},{"id":539396095,"identity":"7fe40eae-ad90-496a-8faf-31868dd41325","order_by":6,"name":"Victor A. Asoala","email":"","orcid":"","institution":"Navrongo Health Research Centre, Ghana Health Service","correspondingAuthor":false,"prefix":"","firstName":"Victor","middleName":"A.","lastName":"Asoala","suffix":""},{"id":539396098,"identity":"51ad70d8-cbb0-4015-8c5a-4b0d9f587810","order_by":7,"name":"Cristina Ariana","email":"","orcid":"","institution":"Genomic Surveillance Unit, Wellcome Sanger Institute, Hinxton","correspondingAuthor":false,"prefix":"","firstName":"Cristina","middleName":"","lastName":"Ariana","suffix":""},{"id":539396100,"identity":"be501d64-89ee-42ac-9a87-4def2a9d64af","order_by":8,"name":"Keziah L. Malm","email":"","orcid":"","institution":"National Malaria Elimination Programme, Ministry of Health","correspondingAuthor":false,"prefix":"","firstName":"Keziah","middleName":"L.","lastName":"Malm","suffix":""},{"id":539396104,"identity":"2e3f0dea-d8f7-41a3-9ba0-7c7b9a0831fb","order_by":9,"name":"Gordon Awandare","email":"","orcid":"","institution":"West African Centre for Cell Biology of Infectious Pathogens (WACCBIP), University of Ghana","correspondingAuthor":false,"prefix":"","firstName":"Gordon","middleName":"","lastName":"Awandare","suffix":""},{"id":539396105,"identity":"326b346a-62a8-45ee-8867-4f85a7d5a4ee","order_by":10,"name":"Alistair Miles","email":"","orcid":"","institution":"Genomic Surveillance Unit, Wellcome Sanger Institute, Hinxton","correspondingAuthor":false,"prefix":"","firstName":"Alistair","middleName":"","lastName":"Miles","suffix":""},{"id":539396106,"identity":"2ef79718-d480-463b-a8c5-f5b150133b77","order_by":11,"name":"Chris S. Clarkson","email":"","orcid":"","institution":"Genomic Surveillance Unit, Wellcome Sanger Institute, Hinxton","correspondingAuthor":false,"prefix":"","firstName":"Chris","middleName":"S.","lastName":"Clarkson","suffix":""},{"id":539396107,"identity":"1b90b59d-91b5-4116-bf6b-bdd294a5e18c","order_by":12,"name":"Lucas N. Amenga-Etego","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABAElEQVRIiWNgGAWjYBACPhTeBzCZwMDA2IBbCxuMwQPEjDNI1sLMQ5QW9tOJjwt+2cnZs7dffGzbVsvAz55jwPh1Bx4tPLmbjWf2JRvz8JwpNs5tO84g2fPGgFn2DD6H5W6T5u1hTuyRyEmTzm07xmBwI8eAWbINjxb+t9t/8/bU1/fIv0mTtgRqsSeoRSJ3GzPPj8MJPBLsx6QZ22oYDCSAfvmIV8vbzdK8DccNe87kMBv2nDvAI3HmWcFhRjxa+PlzN37m+VMtz95+/OGDH2V1cvztyRsf/sSjBQwgZvIYAInD4KiBkHjBHxDB/gBI1EHM+EFQyygYBaNgFIwgAADuCk8ewQwldQAAAABJRU5ErkJggg==","orcid":"","institution":"West African Centre for Cell Biology of Infectious Pathogens (WACCBIP), University of Ghana","correspondingAuthor":true,"prefix":"","firstName":"Lucas","middleName":"N.","lastName":"Amenga-Etego","suffix":""}],"badges":[],"createdAt":"2025-10-16 13:38:12","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-7878288/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-7878288/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1186/s12864-025-12508-7","type":"published","date":"2026-01-07T15:57:36+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":95172887,"identity":"b9a36c49-6a18-4106-a10b-86870fdb24ef","added_by":"auto","created_at":"2025-11-05 06:40:29","extension":"docx","order_by":0,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":9155632,"visible":true,"origin":"","legend":"","description":"","filename":"ghanaformatkelly11025.docx","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/a6af52d27380565cf0d4b4cf.docx"},{"id":95172873,"identity":"5f031344-9d41-4ea2-8156-1b81d3269807","added_by":"auto","created_at":"2025-11-05 06:40:29","extension":"tif","order_by":1,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":292018,"visible":true,"origin":"","legend":"","description":"","filename":"Fig1Map.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/f955af52c0938d50e6443761.tif"},{"id":95172881,"identity":"6401661a-6db2-4197-a69f-5fd0eebacdf5","added_by":"auto","created_at":"2025-11-05 06:40:29","extension":"tif","order_by":2,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":971692,"visible":true,"origin":"","legend":"","description":"","filename":"Fig2PCA.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/fea0abc255ca0e71fe7cdcae.tif"},{"id":95228162,"identity":"60f50ea0-c0cb-40d9-9714-1262c2d56b33","added_by":"auto","created_at":"2025-11-05 16:33:27","extension":"csv","order_by":3,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":803,"visible":true,"origin":"","legend":"","description":"","filename":"Table1samplecollection.csv","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/51ea60a875c513ee3e99254b.csv"},{"id":95172900,"identity":"d8002747-5fe1-4419-8b16-8cb5df221bf8","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"tif","order_by":4,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":820052,"visible":true,"origin":"","legend":"","description":"","filename":"Fig3fst.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/267299128a6bbedad6982575.tif"},{"id":95172882,"identity":"bae4c463-1ee0-483d-9b6c-a05e960b4345","added_by":"auto","created_at":"2025-11-05 06:40:29","extension":"csv","order_by":5,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":630,"visible":true,"origin":"","legend":"","description":"","filename":"Table22La2Rbinversionkaryotypes.csv","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/386dfde7f966b256f4dabdb3.csv"},{"id":95172910,"identity":"34d12cee-8b35-4a78-9f58-e61a18beeb31","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"tif","order_by":6,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":2324850,"visible":true,"origin":"","legend":"","description":"","filename":"Fig4Inversionfreq.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/43c9cd4de22fd34e0dfe897e.tif"},{"id":95172891,"identity":"f0ac40b2-0710-4061-85f2-1063562658ab","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"tif","order_by":7,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1061474,"visible":true,"origin":"","legend":"","description":"","filename":"Fig5targetfreqs.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/bdaf9f95b994384a87bf95b5.tif"},{"id":95226708,"identity":"c259bfd8-3514-4ac0-bf5f-375c93d1cbe7","added_by":"auto","created_at":"2025-11-05 16:31:39","extension":"tif","order_by":8,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1170256,"visible":true,"origin":"","legend":"","description":"","filename":"Fig6metabolicfreq.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/515e601c1275749500faea85.tif"},{"id":95172894,"identity":"382e1f03-9d25-49d6-b46b-c7cdfc2b24a4","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"tif","order_by":9,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1381428,"visible":true,"origin":"","legend":"","description":"","filename":"Fig7gaselection.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/bab628c96d6d2d77325ebd77.tif"},{"id":95226696,"identity":"f8a4a2d2-fafe-4659-b697-d78ff2d0153f","added_by":"auto","created_at":"2025-11-05 16:31:39","extension":"json","order_by":10,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":12720,"visible":true,"origin":"","legend":"","description":"","filename":"a68d8e493f7c43caa2cabe2e5e7e3f95.json","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/60b42fecfb7e0dc8460bd2ff.json"},{"id":95172904,"identity":"eb8819c1-a94a-4d87-a055-c75176b9307a","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"tif","order_by":11,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":373376,"visible":true,"origin":"","legend":"","description":"","filename":"FigS10keapaa.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/9b9408dc450c4df7f9ab6977.tif"},{"id":95172921,"identity":"744f1b0e-6d8d-4d12-8fa0-0d429eca6bec","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"tif","order_by":12,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":407892,"visible":true,"origin":"","legend":"","description":"","filename":"FigS11Supplementaryfigurecoeaediplotypeclustering.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/d9877ca8eca78e50533bcabb.tif"},{"id":95228699,"identity":"c870c3d4-61c4-4d0e-a6d0-91809a2608c0","added_by":"auto","created_at":"2025-11-05 16:34:04","extension":"tif","order_by":13,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":350976,"visible":true,"origin":"","legend":"","description":"","filename":"FigS1WestAfricanCountriesMap.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/1439ad1d3f2d7305f2ce6bbe.tif"},{"id":95228326,"identity":"0d85aad0-4bd3-420b-983c-268ada317534","added_by":"auto","created_at":"2025-11-05 16:33:37","extension":"tif","order_by":14,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":351428,"visible":true,"origin":"","legend":"","description":"","filename":"FigS2Ghananjt.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/87bfc370fc27c59c60f590eb.tif"},{"id":95227622,"identity":"363f3c53-d4da-45fc-a44e-d109107feb64","added_by":"auto","created_at":"2025-11-05 16:32:40","extension":"tif","order_by":15,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1285918,"visible":true,"origin":"","legend":"","description":"","filename":"FigS3diversity.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/f49d03e8116cb14f5f5f6a30.tif"},{"id":95227521,"identity":"4d962eb3-113e-4d3f-ba69-5458478d1833","added_by":"auto","created_at":"2025-11-05 16:32:35","extension":"tif","order_by":16,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":390482,"visible":true,"origin":"","legend":"","description":"","filename":"FigS4Westafricanjt.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/660f75901271b36785469765.tif"},{"id":95172953,"identity":"8d734c14-afc9-44ad-b88e-a7f8b019e4f9","added_by":"auto","created_at":"2025-11-05 06:40:32","extension":"tif","order_by":17,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":122832,"visible":true,"origin":"","legend":"","description":"","filename":"FigS5SupplementaryfigureRdlinversion.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/39846bc809849fa7c7662e5f.tif"},{"id":95227524,"identity":"2f25324c-8b4f-47da-ad7b-062fcbda4169","added_by":"auto","created_at":"2025-11-05 16:32:35","extension":"tif","order_by":18,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1343040,"visible":true,"origin":"","legend":"","description":"","filename":"FigS6SignalsChr2.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/708f414abf14740918197a3c.tif"},{"id":95226920,"identity":"8dbfc71e-827d-4540-9620-19c4fcb3f142","added_by":"auto","created_at":"2025-11-05 16:31:53","extension":"tif","order_by":19,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":996192,"visible":true,"origin":"","legend":"","description":"","filename":"FigS7SignalsChr3.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/05e59dd74a8af08a0b39e28e.tif"},{"id":95226669,"identity":"5b6b86b5-9e5a-4f39-b3ca-8ef3598a2c9f","added_by":"auto","created_at":"2025-11-05 16:31:36","extension":"tif","order_by":20,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1542520,"visible":true,"origin":"","legend":"","description":"","filename":"FigS8signalsChrX.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/29806d29f144c331c5942025.tif"},{"id":95228205,"identity":"ce0338b1-d816-423e-a161-a651ebba937f","added_by":"auto","created_at":"2025-11-05 16:33:28","extension":"tif","order_by":21,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":346410,"visible":true,"origin":"","legend":"","description":"","filename":"FigS9keap1diplotypeclustering.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/7535abe24757364abc2568f4.tif"},{"id":95227770,"identity":"a9b1f661-363e-4371-932e-6df1b0351231","added_by":"auto","created_at":"2025-11-05 16:32:56","extension":"csv","order_by":22,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":125394,"visible":true,"origin":"","legend":"","description":"","filename":"TableS1samplingtable.csv","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/ae262ba156e39507219928cb.csv"},{"id":95312470,"identity":"81d349ea-574f-4a3c-adf3-8f53848ec229","added_by":"auto","created_at":"2025-11-06 15:49:29","extension":"xml","order_by":23,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":237818,"visible":true,"origin":"","legend":"","description":"","filename":"a68d8e493f7c43caa2cabe2e5e7e3f951enriched.xml","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/0ed405b5b6be37d560726056.xml"},{"id":95172937,"identity":"49877bb4-b26f-4434-923a-0ff1fa0453a1","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"tif","order_by":24,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":292018,"visible":true,"origin":"","legend":"","description":"","filename":"Fig1Map.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/c11fdd111660754d8f4c92ee.tif"},{"id":95172903,"identity":"0500645f-2b5e-44ad-a57a-1314651ce92f","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"tif","order_by":25,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":971692,"visible":true,"origin":"","legend":"","description":"","filename":"Fig2PCA.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/f6a72230f9f0c35ca2d1e072.tif"},{"id":95226800,"identity":"72ced38f-2f84-494e-956b-ca6d9697ca23","added_by":"auto","created_at":"2025-11-05 16:31:44","extension":"tif","order_by":26,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":820052,"visible":true,"origin":"","legend":"","description":"","filename":"Fig3fst.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/352db57c090db1964fcf64e6.tif"},{"id":95228917,"identity":"ab0aa618-9bdb-4034-a8b7-8182e9082146","added_by":"auto","created_at":"2025-11-05 16:34:16","extension":"tif","order_by":27,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":2324850,"visible":true,"origin":"","legend":"","description":"","filename":"Fig4Inversionfreq.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/4f4ededc8ca41fbcd4f92dee.tif"},{"id":95226773,"identity":"185cfd04-da7e-4761-ba5a-759898135d8f","added_by":"auto","created_at":"2025-11-05 16:31:43","extension":"tif","order_by":28,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1061474,"visible":true,"origin":"","legend":"","description":"","filename":"Fig5targetfreqs.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/cfbbc34de083c14fb510ccd1.tif"},{"id":95172933,"identity":"cba24a65-1e01-458c-b20c-5880df5c0f54","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"tif","order_by":29,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1170256,"visible":true,"origin":"","legend":"","description":"","filename":"Fig6metabolicfreq.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/a6e65c26e86bb9a9a32b571a.tif"},{"id":95227659,"identity":"10ee0a1d-a9fd-4ceb-b246-4ddc6a81968b","added_by":"auto","created_at":"2025-11-05 16:32:42","extension":"tif","order_by":30,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1381428,"visible":true,"origin":"","legend":"","description":"","filename":"Fig7gaselection.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/0f9cad3d553caafd4fec3b0f.tif"},{"id":95172935,"identity":"324afdfc-c9aa-41dc-87a3-3472ec428d7e","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"jpeg","order_by":31,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":371148,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage1.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/9148abdd58e8d6a3ae3fb1d7.jpeg"},{"id":95172923,"identity":"1b6fc9c0-80c4-4435-9c98-cbfe2e4d902c","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"jpeg","order_by":32,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":587405,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage10.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/a05077c3fcc4421052344207.jpeg"},{"id":95172959,"identity":"6cbc4ede-37aa-4846-8cd9-a2f396aa4490","added_by":"auto","created_at":"2025-11-05 06:40:32","extension":"png","order_by":33,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":377357,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage11.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/16f45583b76a2fe45e464421.png"},{"id":95172898,"identity":"260a9d5d-598c-41c9-8866-658c9f9e59eb","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"png","order_by":34,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":76534,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage12.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/0eb67984964b2579b1c7f05f.png"},{"id":95172925,"identity":"e362e705-e9cd-4243-b1cc-c6e1b4cc7ce3","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"jpeg","order_by":35,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":799025,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage13.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/b7fa35293528f6089a087819.jpeg"},{"id":95172929,"identity":"43bf18c2-57bf-43d7-b733-467f4f7e2802","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"jpeg","order_by":36,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":726620,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage14.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/a58cd638ac806d571d41ab75.jpeg"},{"id":95227882,"identity":"1b0ea78f-009b-4423-a549-ec3bb4626157","added_by":"auto","created_at":"2025-11-05 16:33:03","extension":"jpeg","order_by":37,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":826105,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage15.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/7745a39e6839f5fac1829b59.jpeg"},{"id":95226654,"identity":"3a6113a3-15b5-49bb-95d5-35f6375b94c8","added_by":"auto","created_at":"2025-11-05 16:31:36","extension":"png","order_by":38,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":190632,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage16.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/e16b7e447d5527f796a4f9cd.png"},{"id":95226540,"identity":"563d8579-b531-444e-bac6-d80aabdf4788","added_by":"auto","created_at":"2025-11-05 16:31:21","extension":"png","order_by":39,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":240703,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage17.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/7444f099f24117c25e0d44dc.png"},{"id":95228655,"identity":"d3416ed4-b205-4c7b-99e0-89dceb421a97","added_by":"auto","created_at":"2025-11-05 16:34:01","extension":"png","order_by":40,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":232021,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage18.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/6a0d74597bdde2e5e87e2bcb.png"},{"id":95172939,"identity":"6e9fcc51-0b83-4bdf-940b-9f9a9e4c0441","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"jpeg","order_by":41,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":375724,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage2.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/ff6d9c014518e87610088b80.jpeg"},{"id":95228295,"identity":"9ec58997-35ed-415d-9607-1c668eff089d","added_by":"auto","created_at":"2025-11-05 16:33:35","extension":"jpeg","order_by":42,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":609799,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage3.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/d9d0cb25379f0c162de90892.jpeg"},{"id":95172927,"identity":"4e5c49ba-1ad9-4970-aba0-a2d22786b4b5","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"png","order_by":43,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":978057,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage4.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/af0dd8352ba1e6b3c5c8ba5a.png"},{"id":95226618,"identity":"e7d04676-797c-4d7f-bb5e-7a31344fedb5","added_by":"auto","created_at":"2025-11-05 16:31:28","extension":"jpeg","order_by":44,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":253855,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage5.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/c420f7edc97fd2e5ba9e86a6.jpeg"},{"id":95228734,"identity":"b3da9da2-217c-48f8-9800-507c5b1ac4a4","added_by":"auto","created_at":"2025-11-05 16:34:05","extension":"jpeg","order_by":45,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":263959,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage6.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/f70fedbe27fd94948326b23b.jpeg"},{"id":95226810,"identity":"65765f0a-dfc5-48a7-8a3f-bd6e69cf41e4","added_by":"auto","created_at":"2025-11-05 16:31:44","extension":"png","order_by":46,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":875499,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage7.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/3b155112a84525d8df739420.png"},{"id":95172905,"identity":"120c22c5-8b0c-41b1-9451-ab1de0a3d343","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"png","order_by":47,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":273048,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage8.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/2a645524eac2da419471f4b9.png"},{"id":95226828,"identity":"01d44576-e1ec-49c9-b3c2-eca8e4beb7d4","added_by":"auto","created_at":"2025-11-05 16:31:45","extension":"png","order_by":48,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":264694,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage9.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/ae49dea8fafb980cfb4b5e7d.png"},{"id":95226422,"identity":"a423cbe1-4de7-44fd-b93f-e5a58bda3a70","added_by":"auto","created_at":"2025-11-05 16:31:05","extension":"png","order_by":49,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":88185,"visible":true,"origin":"","legend":"","description":"","filename":"OnlineFig1Map.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/6a0179fec7f39293210e1280.png"},{"id":95226484,"identity":"8e8ac0cd-7d34-4718-9ab0-1252097d2cde","added_by":"auto","created_at":"2025-11-05 16:31:13","extension":"png","order_by":50,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":110026,"visible":true,"origin":"","legend":"","description":"","filename":"OnlineFig2PCA.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/90a90a0186e0b3cfdc589a0b.png"},{"id":95172916,"identity":"b2cfb5be-142e-4106-b4ae-b3946d09dc9d","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"png","order_by":51,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":94499,"visible":true,"origin":"","legend":"","description":"","filename":"OnlineFig3fst.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/953a05c1add310cb1d771b4f.png"},{"id":95228154,"identity":"c439f8ce-a200-47dc-acb4-09fbaf1edc98","added_by":"auto","created_at":"2025-11-05 16:33:27","extension":"png","order_by":52,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":124770,"visible":true,"origin":"","legend":"","description":"","filename":"OnlineFig4Inversionfreq.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/946f30cf65670ad244a180c6.png"},{"id":95172909,"identity":"ad890f7d-8ecd-44cd-ab4c-0001fdcf13a8","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"png","order_by":53,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":143102,"visible":true,"origin":"","legend":"","description":"","filename":"OnlineFig5targetfreqs.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/5a12b59c19519f9ade1a54e5.png"},{"id":95226727,"identity":"50d6fcaf-4595-4cbf-b7a3-c8d3e21bb0e8","added_by":"auto","created_at":"2025-11-05 16:31:40","extension":"png","order_by":54,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":160074,"visible":true,"origin":"","legend":"","description":"","filename":"OnlineFig6metabolicfreq.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/b5d4b19d57707c1227c0fa33.png"},{"id":95226758,"identity":"0fede5b1-f11c-4740-941b-a64e7e16eb5d","added_by":"auto","created_at":"2025-11-05 16:31:43","extension":"png","order_by":55,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":103171,"visible":true,"origin":"","legend":"","description":"","filename":"OnlineFig7gaselection.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/10fd54a7f6f1a1e0c6538170.png"},{"id":95226811,"identity":"e820bbc2-b296-4a95-848c-1f1b041af44f","added_by":"auto","created_at":"2025-11-05 16:31:44","extension":"png","order_by":56,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":114556,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage1.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/978a51148fff38be16ef07d0.png"},{"id":95172918,"identity":"b95a95bb-279e-4bd1-97db-25ff57198cd1","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"png","order_by":57,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":117729,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage10.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/086f099bb70a83d721ae9722.png"},{"id":95226516,"identity":"b14f6f13-a7d3-413c-ba81-7bc1df26a59e","added_by":"auto","created_at":"2025-11-05 16:31:18","extension":"png","order_by":58,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":50598,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage11.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/028c673920ace87348421dec.png"},{"id":95172931,"identity":"1c3319b4-825f-44bf-a144-82fc2f336e32","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"png","order_by":59,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":23114,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage12.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/4191709779e25fdfe798c9fd.png"},{"id":95172942,"identity":"1e0b65c6-0934-4546-a7c2-c92dfa9e42d1","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"png","order_by":60,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":149843,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage13.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/145176108dc9815735436793.png"},{"id":95172955,"identity":"d9780cea-e24a-48be-a700-41ade83bd418","added_by":"auto","created_at":"2025-11-05 06:40:32","extension":"png","order_by":61,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":120290,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage14.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/4e55210450044be475c56688.png"},{"id":95172965,"identity":"651cd585-cb75-4a3a-98a1-cec4f10edb48","added_by":"auto","created_at":"2025-11-05 06:40:32","extension":"png","order_by":62,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":161359,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage15.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/13d4fb2f6894ae091fdd0e7a.png"},{"id":95227783,"identity":"bf2d5787-3e74-4549-9433-9a1efdd01523","added_by":"auto","created_at":"2025-11-05 16:32:56","extension":"png","order_by":63,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":43953,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage16.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/94566c07a3573ba614546e7a.png"},{"id":95226672,"identity":"3ebb3a14-6411-4971-9da7-13c298d32876","added_by":"auto","created_at":"2025-11-05 16:31:36","extension":"png","order_by":64,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":44914,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage17.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/1b3fdef7de58d2daccd652c7.png"},{"id":95226877,"identity":"daea9c47-d4a2-436d-a40c-437bad82695e","added_by":"auto","created_at":"2025-11-05 16:31:49","extension":"png","order_by":65,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":46968,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage18.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/5d050d20d45be6a6cf192c6d.png"},{"id":95172936,"identity":"6b262cd5-51d3-4c87-a168-7eecb3a9e75f","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"png","order_by":66,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":85971,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage2.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/4b955c48092b432d99667c0b.png"},{"id":95172961,"identity":"eb7f4849-7358-41a6-bfd5-1ed7f13e664b","added_by":"auto","created_at":"2025-11-05 06:40:32","extension":"png","order_by":67,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":137628,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage3.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/07b3c2286ef0d43971fd0a4b.png"},{"id":95228298,"identity":"ad7c0672-407b-4068-b706-f414b6a0853c","added_by":"auto","created_at":"2025-11-05 16:33:35","extension":"png","order_by":68,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":123486,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage4.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/06f33ee6729df47fd4162308.png"},{"id":95172951,"identity":"b44ea445-c072-41d8-8465-7b7c93103c57","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"png","order_by":69,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":186236,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage5.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/7df0a9c8278bcdf68ba4686f.png"},{"id":95172971,"identity":"c7ae3274-cc6a-4ca9-b8e6-ad03b5a370a8","added_by":"auto","created_at":"2025-11-05 06:40:32","extension":"png","order_by":70,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":208782,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage6.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/fdb51a74c48e6d4401663f21.png"},{"id":95172945,"identity":"bd54a5bc-e038-431a-9c5c-f7790783860b","added_by":"auto","created_at":"2025-11-05 06:40:31","extension":"png","order_by":71,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":102888,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage7.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/8c1ee7f0a011e35f5e9f9909.png"},{"id":95228193,"identity":"074e625c-cd50-4a5c-8516-0776090fd7dd","added_by":"auto","created_at":"2025-11-05 16:33:28","extension":"png","order_by":72,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":63889,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage8.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/b41474f96f4b076d00efa2d8.png"},{"id":95228158,"identity":"bb298461-29d4-4877-adc0-f9650b86a172","added_by":"auto","created_at":"2025-11-05 16:33:27","extension":"png","order_by":73,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":70252,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage9.png","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/2dee7aa3949ed1ed37f5997b.png"},{"id":95172957,"identity":"ee461552-1378-4577-a406-17abdc5f01b4","added_by":"auto","created_at":"2025-11-05 06:40:32","extension":"xml","order_by":74,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":239137,"visible":true,"origin":"","legend":"","description":"","filename":"a68d8e493f7c43caa2cabe2e5e7e3f951structuring.xml","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/54b3a48b11fab04148db7bf3.xml"},{"id":95172967,"identity":"98183a8b-d7d2-4ffd-8bfc-2f0ad1b41d7d","added_by":"auto","created_at":"2025-11-05 06:40:32","extension":"html","order_by":75,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":257212,"visible":true,"origin":"","legend":"","description":"","filename":"earlyproof.html","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/0ddfadc3801567b4efa5e6f6.html"},{"id":95172872,"identity":"04753b43-4a23-4d7c-9ddc-336d8fc31f7f","added_by":"auto","created_at":"2025-11-05 06:40:29","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":553957,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eA\u003c/strong\u003e: Map of West Africa highlighting countries included in the analysis including\u003cstrong\u003e \u003c/strong\u003eGhana (green), Burkina Faso (light blue), Côte d'Ivoire (gray), and Mali (dark blue). \u003cstrong\u003eB:\u003c/strong\u003eMap of Ghana and its bioclimatic zones including the Northern Savannah (white), Middle Belt (light gray), and Southern Forest (dark gray).\u003c/p\u003e","description":"","filename":"Fig1Map.tif.jpg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/0149e6d92aa6ab0393b6550f.jpg"},{"id":95227043,"identity":"fcf6e613-7f9d-451e-a108-7b13ed80577e","added_by":"auto","created_at":"2025-11-05 16:32:01","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":723605,"visible":true,"origin":"","legend":"\u003cp\u003ePlotting from Principal Components Analysis to investigate population structure in \u003cem\u003eAn. coluzzii\u003c/em\u003e collections from \u003cstrong\u003eA\u003c/strong\u003e. Ghana and \u003cstrong\u003eB\u003c/strong\u003e. from Ghana and the surrounding West African countries of Burkina Faso, Cote d’Ivoire and Mali using 10,000 SNPs on chromosome 3L. \u0026nbsp;Grey circles represent 95% confidence centroid ellipses of two K-means clusters.\u003c/p\u003e","description":"","filename":"Fig2PCA.tif.jpg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/c89e2d81c82c5719085e9d77.jpg"},{"id":95172874,"identity":"53601ff5-d927-4717-bd44-660e127655b1","added_by":"auto","created_at":"2025-11-05 06:40:29","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":496654,"visible":true,"origin":"","legend":"\u003cp\u003eHudson’s pairwise F\u003csub\u003eST\u003c/sub\u003e between population cohorts from Ghana and nearby West and Central African countries including Burkina Faso, Mali and Cote d’Ivoire. \u0026nbsp;Lighter shades of green indicate higher values of F\u003csub\u003eST\u003c/sub\u003e.\u003c/p\u003e","description":"","filename":"Fig3fst.tif.jpg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/1488aefe97682e81f5a52b35.jpg"},{"id":95172879,"identity":"53f0ccf7-2438-47fe-abb3-2dff3c32dea0","added_by":"auto","created_at":"2025-11-05 06:40:29","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":993147,"visible":true,"origin":"","legend":"\u003cp\u003eThe frequency of karyotypes of the A. 2La and B. 2Rb inversion region in \u003cem\u003eAn. coluzzii\u003c/em\u003e from Ghana (GH) and the surrounding West African countries of Burkina Faso (BF), Mali (ML), and Cote d’Ivoire (CI). Southern forest locations include the Western, Central, Greater Accra, Eastern and Ashanti regions in Ghana and Lagunes in Cote d’Ivoire. Northern Savannah locations include the Upper East region in Ghana and locations in Burkina Faso, Mali, and Cote d’Ivoire.\u003c/p\u003e","description":"","filename":"Fig4Inversionfreq.tif.jpg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/0eaec33002a6242d8d964d71.jpg"},{"id":95227069,"identity":"f5252208-8e37-40bb-af08-682b95376edd","added_by":"auto","created_at":"2025-11-05 16:32:02","extension":"jpg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":732560,"visible":true,"origin":"","legend":"\u003cp\u003eAmino acid and CNV frequencies at target site insecticide resistance loci including \u003cem\u003eVgsc\u003c/em\u003e, \u003cem\u003eRdl \u003c/em\u003eand \u003cem\u003eAce1 \u003c/em\u003ein population cohorts from Ghana. \u0026nbsp;Darker shades of red indicate higher frequencies.\u003c/p\u003e","description":"","filename":"Fig5targetfreqs.tif.jpg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/d792f36a25b4203d8dd7efe9.jpg"},{"id":95226677,"identity":"43ce3153-c322-49b0-873e-b1a4956a4be3","added_by":"auto","created_at":"2025-11-05 16:31:37","extension":"jpg","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":750218,"visible":true,"origin":"","legend":"\u003cp\u003eCNV and amino acid frequencies at metabolic insecticide resistance loci including the \u003cem\u003eCyp6 \u003c/em\u003egene cluster, \u003cem\u003eCyp6m2, Cyp9k1, \u003c/em\u003ethe carboxylesterase gene cluster \u003cem\u003eCoeae2-7G \u003c/em\u003eand \u003cem\u003eGste\u003c/em\u003e. \u0026nbsp;Darker shades of red indicate higher frequencies.\u003c/p\u003e","description":"","filename":"Fig6metabolicfreq.tif.jpg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/2e2c4411629ac15212fe6341.jpg"},{"id":95228871,"identity":"d3c74106-c609-4297-bdb1-f92782fdd072","added_by":"auto","created_at":"2025-11-05 16:34:13","extension":"jpg","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":797382,"visible":true,"origin":"","legend":"\u003cp\u003eH12 selection scans for chromosomes 2RL, 3RL and X of the \u003cem\u003eAn. coluzzii \u003c/em\u003epopulation from the southern region of Greater Accra region in 2018. High values of H12 indicate a signal of positive selection. Signals at both known and novel loci putatively involved in insecticide resistance are annotated on the plot.\u003c/p\u003e","description":"","filename":"Fig7gaselection.tif.jpg","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/66584e0677c67d5384281fb4.jpg"},{"id":100070201,"identity":"05f35ae6-1b02-49b0-b0af-952fac7ff67b","added_by":"auto","created_at":"2026-01-12 16:16:57","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":6258121,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/62289fae-b807-4272-94d1-d8711a255b1c.pdf"},{"id":95172871,"identity":"1780c646-f2b6-43c6-b4de-bc7c51e4b3f4","added_by":"auto","created_at":"2025-11-05 06:40:29","extension":"csv","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":803,"visible":true,"origin":"","legend":"","description":"","filename":"Table1samplecollection.csv","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/b183881d1b76e7c862f4d035.csv"},{"id":95172877,"identity":"5c496eb7-e422-4112-a366-ce3a7ea4566f","added_by":"auto","created_at":"2025-11-05 06:40:29","extension":"csv","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":630,"visible":true,"origin":"","legend":"","description":"","filename":"Table22La2Rbinversionkaryotypes.csv","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/fa05c13af74afece136dcf39.csv"},{"id":95228851,"identity":"f710a2b9-3cc3-4f89-a553-51498e38fd90","added_by":"auto","created_at":"2025-11-05 16:34:11","extension":"csv","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":125394,"visible":true,"origin":"","legend":"","description":"","filename":"TableS1samplingtable.csv","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/a83d482111e0fdc1b2b0e442.csv"},{"id":95172885,"identity":"b88eb363-c893-4425-b967-30287f309ff4","added_by":"auto","created_at":"2025-11-05 06:40:29","extension":"tif","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":122832,"visible":true,"origin":"","legend":"","description":"","filename":"FigS5SupplementaryfigureRdlinversion.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/331e85515a7ab21d9249522c.tif"},{"id":95228168,"identity":"74d0df40-3ad5-4bc7-813a-d1dc14996271","added_by":"auto","created_at":"2025-11-05 16:33:27","extension":"tif","order_by":4,"title":"","display":"","copyAsset":false,"role":"supplement","size":346410,"visible":true,"origin":"","legend":"","description":"","filename":"FigS9keap1diplotypeclustering.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/e5735cbc736acde59bd91745.tif"},{"id":95172890,"identity":"2a287f0f-f113-46fb-8f9c-e09242069d1f","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"tif","order_by":5,"title":"","display":"","copyAsset":false,"role":"supplement","size":373376,"visible":true,"origin":"","legend":"","description":"","filename":"FigS10keapaa.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/385a87e7fa529f57b75a0fd2.tif"},{"id":95172888,"identity":"b89bda94-e687-424b-9a48-562318300817","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"tif","order_by":6,"title":"","display":"","copyAsset":false,"role":"supplement","size":407892,"visible":true,"origin":"","legend":"","description":"","filename":"FigS11Supplementaryfigurecoeaediplotypeclustering.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/63284ae2e851422fda6f7a44.tif"},{"id":95172915,"identity":"38c0ca4d-75ec-46f8-943e-7f0b626c9dd1","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"tif","order_by":7,"title":"","display":"","copyAsset":false,"role":"supplement","size":351428,"visible":true,"origin":"","legend":"","description":"","filename":"FigS2Ghananjt.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/5dbb5b975140c373bd0ff9b1.tif"},{"id":95172897,"identity":"727b13e4-5412-4178-a881-341da7ecfa91","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"tif","order_by":8,"title":"","display":"","copyAsset":false,"role":"supplement","size":1285918,"visible":true,"origin":"","legend":"","description":"","filename":"FigS3diversity.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/b6b91d7780f5b49c9ac47a75.tif"},{"id":95226923,"identity":"09d6cd3a-2ebd-4869-a01e-b0d6958f3235","added_by":"auto","created_at":"2025-11-05 16:31:53","extension":"tif","order_by":9,"title":"","display":"","copyAsset":false,"role":"supplement","size":996192,"visible":true,"origin":"","legend":"","description":"","filename":"FigS7SignalsChr3.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/14620a2cbdb427e32046cc58.tif"},{"id":95172906,"identity":"cc6e2870-447b-456c-8070-c00d33b0f5f8","added_by":"auto","created_at":"2025-11-05 06:40:30","extension":"tif","order_by":10,"title":"","display":"","copyAsset":false,"role":"supplement","size":350976,"visible":true,"origin":"","legend":"","description":"","filename":"FigS1WestAfricanCountriesMap.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/693599ea8d50361b2876dd5e.tif"},{"id":95228671,"identity":"c90a923c-18ba-4a41-8ceb-228d634bc507","added_by":"auto","created_at":"2025-11-05 16:34:02","extension":"tif","order_by":11,"title":"","display":"","copyAsset":false,"role":"supplement","size":390482,"visible":true,"origin":"","legend":"","description":"","filename":"FigS4Westafricanjt.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/17d4e53cffee3b15f802a4ed.tif"},{"id":95228872,"identity":"e2621f95-bf52-47a1-bc6b-ef8e0f6a1e72","added_by":"auto","created_at":"2025-11-05 16:34:13","extension":"tif","order_by":12,"title":"","display":"","copyAsset":false,"role":"supplement","size":1343040,"visible":true,"origin":"","legend":"","description":"","filename":"FigS6SignalsChr2.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/f9f3c0ccae148e316eea3f70.tif"},{"id":95226719,"identity":"b0901b9f-6f22-44d5-b36d-ff8b5bc49441","added_by":"auto","created_at":"2025-11-05 16:31:40","extension":"tif","order_by":13,"title":"","display":"","copyAsset":false,"role":"supplement","size":1542520,"visible":true,"origin":"","legend":"","description":"","filename":"FigS8signalsChrX.tif","url":"https://assets-eu.researchsquare.com/files/rs-7878288/v1/4ae8a04fc32eae17e00822b9.tif"}],"financialInterests":"No competing interests reported.","formattedTitle":"Genomic population structure and insecticide resistance mechanisms in the malaria vector An. coluzzii across contrasting bioclimatic zones in West Africa","fulltext":[{"header":"Introduction","content":"\u003cp\u003eInsect vectors with wide geographical distributions are subject to varying selection pressures driven by diverse ecological conditions. Both geographical barriers and ecological selection on insect vectors can promote variation in physiological and behavioral traits that impact human disease transmission. For example, transmission can be impacted by differences in geographical species distribution, blood-feeding behavior, dispersal and habitat use[\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e]. Ecological divergence often impacts on gene flow to form genetically differentiated ecotypes or species with different abilities to act as disease vectors e.g., tsetse flies that are vectors for sleeping sickness[\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e] or mosquitoes that transmit malaria[\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]. Differences in population connectivity underpinned by vector ecology also impacts the sharing of adaptive alleles across geographical space. For example, these could be adaptations which promote behavioral avoidance or resistance mechanisms, such as insecticide target site or metabolic resistance, which compromise the population control strategies public health relies on[\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e]. There is therefore a need to understand how geography, ecology and anthropogenetic selection pressures, such as the application of insecticides, interact to influence the effectiveness of disease reduction strategies.\u003c/p\u003e\u003cp\u003eGhana, with its rich tapestry of ecological landscapes ranging from coastal Savannah, through dense tropical rainforests, to Sahelian Savannah is representative of the north to south environmental gradient that spans the countries of central West Africa. Its diverse ecological conditions impact the distribution of disease vectors, including major malaria vectors in the \u003cem\u003eAnopheles\u003c/em\u003e mosquito genus[\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e]. Malaria burden and seasonal peaks in incidence vary markedly across Ghana resulting in ~\u0026thinsp;500,000 infections a year[\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. The key malaria vectors in central West Africa include members of the \u003cem\u003eAnopheles gambiae s.l.\u003c/em\u003e complex and \u003cem\u003eAnopheles funestus s.l\u003c/em\u003e species group. Within the former is the major malaria vector \u003cem\u003eAn. coluzzii\u003c/em\u003e. Due to its ecological plasticity, it has a widespread distribution across Ghana including in the northern Savannah where the malaria burden is highest[\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e, \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e]. It also exhibits considerable variation in physiological and behavioral traits across the different bioclimatic zones of the West African subregion. These include differences in mating behaviour[\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e], dispersal strategies[\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e] and phenotypes promoting thermal tolerance including entering into a state of dormancy during dry conditions, i.e., aestivation[\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]. \u003cem\u003eAn. coluzzii\u003c/em\u003e also showcases genomic variability across contrasting environmental conditions. For example, a striking cline in 2La and 2Rb inversion frequencies across the bioclimatic zones of West and Central Africa have been observed in Nigeria, Burkina Faso and Cameroon[\u003cspan additionalcitationids=\"CR13 CR14\" citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e]. While populations in arid, savanna regions exhibit higher frequencies of the inverted 2La and 2Rb karyotypes, those in humid, forested areas show lower frequencies[\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e, \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e]. Experimental work has also associated these inversion karyotypes with feeding behavior[\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e] and thermal tolerance[\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e]. Furthermore, these karyotypes can impact resistance to \u003cem\u003ePlasmodium\u003c/em\u003e parasite infection[\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e] and to insecticides used for population control[\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e, \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e], highlighting that environmental adaptation can impact directly on factors influencing malaria transmission.\u003c/p\u003e\u003cp\u003eA major hurdle to malaria control is the rapid evolution of insecticide resistance in vector populations[\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e] which undermines the effectiveness of control tools including insecticide-treated nets (ITNs) and indoor residual spraying (IRS)[\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e, \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e]. Resistance mechanisms include target site resistance and metabolic resistance (e.g., heightened activity of detoxification enzymes). Target site resistance is conferred by non-synonymous substitutions in the voltage-gated sodium channel (\u003cem\u003eVgsc\u003c/em\u003e) gene[\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e, \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e], coding for the protein targeted by widely-used pyrethroid insecticides, the acetylcholinesterase (\u003cem\u003eAce-1\u003c/em\u003e) gene[\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e, \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e] which codes for the protein target of carbamates and organophosphates, and the GABA receptor subunit \u003cem\u003eRdl\u003c/em\u003e gene, coding for the target of dieldrin[\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e, \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e]. Metabolic resistance is provided by increased activity of enzyme products determined by key gene classes such as cytochrome P450 monooxygenases (P450s)[\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e, \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e], esterases including carboxylesterases (Coe)[\u003cspan additionalcitationids=\"CR35\" citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e], and glutathione S-transferases (Gst) [\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e]and uridine diphosphate (UDP)-glycosyltransferases[\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e]. Research conducted in Burkina Faso and Ghana has revealed that \u003cem\u003eAn. coluzzii\u003c/em\u003e populations in urban areas often exhibit higher levels of insecticide resistance phenotypes compared to rural populations[\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e, \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e]. The urban/rural disparity has been attributed to differential exposure to insecticides used in agriculture and public health in peri-urban and urban settings. The distribution of insecticide resistance mechanisms is therefore expected to vary across geographical space and with environmental characteristics. Moreover, resistance loci like target-site substitutions and copy number variants at metabolic resistance loci vary geographically, often reflecting localized selection pressures[\u003cspan additionalcitationids=\"CR41\" citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e]. Furthermore, scans for signals of recent selection have highlighted potential differences in insecticide resistance loci including \u003cem\u003eVgsc\u003c/em\u003e and cytochrome P450\u0026rsquo;s between forest and non-forest ecosystems in southern Ghana, suggesting differential selection pressures in the contrasting bioclimates[\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e]. Understanding the geographic distribution of these resistance mechanisms and how environmental factors influence the sharing of adaptive alleles is crucial for designing targeted and effective control strategies.\u003c/p\u003e\u003cp\u003eAdvances in genomics have revolutionized our ability to investigate population structure and molecular insecticide resistance mechanisms in malaria vectors at high resolution[\u003cspan additionalcitationids=\"CR45\" citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e]. Here we use whole-genome sequencing of samples collected in Ghana, alongside data available in the MalariaGEN Vector Observatory (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.malariagen.net/vobs/\u003c/span\u003e\u003cspan address=\"https://www.malariagen.net/vobs/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), to elucidate population structure across the different bioclimatic zones of West Africa and investigate if ecological variation promotes population structure in \u003cem\u003eAn. coluzzii\u003c/em\u003e. Using SNP, haplotype and CNV variation data, we also explore how molecular insecticide resistance varies across geographical space using Ghana as a case study to determine if connected populations are more likely to share adaptive alleles resulting from the selection pressures imposed by front-line public health.\u003c/p\u003e"},{"header":"Methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e\u003ch2\u003eMosquito sampling and identification\u003c/h2\u003e\u003cp\u003eAdult wild caught mosquitoes were sampled across Ghana from different ecological zones in cross-sectional studies from 2016 to 2018. Mosquito sampling was done in Navrongo, Upper East Region of the northern Savannah, Adansi in the Ashanti Region and Koforidua in the Eastern Region of the middle transitional forest, Madina in Greater Accra and Takoradi in Western Region of the coastal Savannah of Ghana. \u003cem\u003eAnopheline\u003c/em\u003e mosquitoes were caught using Human Landing Catches (HLCs) at night, CDC light traps hanged overnight, and pyrethroid spray collections of resting mosquitoes in the early mornings. Adult female \u003cem\u003eAn. gambiae s.l\u003c/em\u003e were identified using morphological keys[\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e] and stored singly in 70% ethanol in 96 well PCR plates. These were then shipped to the Wellcome Sanger Institute Genomic Surveillance Unit. DNA was extracted using the Qiagen DNeasy blood and Tissue Kit (Qiagen Science, MD, USA) according to the manufacturer\u0026rsquo;s instructions before sequencing and genotyping. Genomic data from other West African countries were openly available through the Vector Observatory (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.malariagen.net/vobs/\u003c/span\u003e\u003cspan address=\"https://www.malariagen.net/vobs/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), a collaborative project to obtain genomics data of \u003cem\u003eAnopheles\u003c/em\u003e mosquitoes for malaria control.\u003c/p\u003e\u003c/div\u003e\n\u003ch3\u003eSequencing and variant calling\u003c/h3\u003e\n\u003cp\u003eUsing previously described protocols described by the \u003cem\u003eAnopheles gambiae\u003c/em\u003e 1000 Genomes (Ag1000G) phase 3 project (The Anopheles gambiae 1000 Genomes Consortium 2018, 2020), Illumina paired-end sequencing was performed with HiSeq 2000 and HiSeq X technologies. Briefly, 100-150bp sequencing reads were aligned to the AgamP4 PEST reference genome using BWA version 0.7.15[\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e] and single nucleotide polymorphisms (SNP) called using GATK version 3.7.0[\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e]. The resulting data was quality controlled to only include individuals with \u0026ge;\u0026thinsp;10X overall median coverage and with data across greater than 50% of the reference genome. Samples identified as cross-contaminated by a percentage of \u0026ge;\u0026thinsp;4.5% were excluded as defined by Ag1000G phase 3 project protocols. Sites where SNP calling and genotyping were expected to be unreliable based on previous analyses of Mendelian inheritance in laboratory crosses were also excluded from analysis.\u003c/p\u003e\u003cp\u003eHaplotypes were phased using a WDL implementation of read-backed phasing using WhatsHap v1.0[\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e] and statistical phasing using SHAPEIT v4.2.1[\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e] (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://github.com/malariagen/pipelines\u003c/span\u003e\u003cspan address=\"https://github.com/malariagen/pipelines\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003cspan type=\"Underline\" class=\"Underline\" name=\"Emphasis\"\u003e).\u003c/span\u003e Copy number variants (CNVs) for each individual were called based on the copy number state inferred across windows of the genome. The copy number state was normalised using a Gaussian hidden Markov model (HMM) implemented by hmmlearn (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://github.com/hmmlearn/hmmlearn\u003c/span\u003e\u003cspan address=\"https://github.com/hmmlearn/hmmlearn\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) as described previously[\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e]. A CNV was called across regions where at least five contiguous genomic windows had a predicted copy number\u0026thinsp;\u0026gt;\u0026thinsp;2 or \u0026gt;\u0026thinsp;1 for the X chromosome in males. To increase reliability, only CNV calls with a high likelihood\u0026thinsp;\u0026gt;\u0026thinsp;1000 and low coverage variance\u0026thinsp;\u0026lt;\u0026thinsp;0.35 based on the HMM were retained. All preceding analyses were performed with the malariagen_data python package.\u003c/p\u003e\n\u003ch3\u003eTaxonomic status\u003c/h3\u003e\n\u003cp\u003eTaxa were provisionally assigned to the samples using Ancestry Informative Markers (AIMS). These are a set of SNPs previously described by the \u003cem\u003eAnopheles gambiae\u003c/em\u003e 1000 Genomes project as exclusive to each taxonomic group based on data generated from the \u003cem\u003eAnopheles\u003c/em\u003e 16 genomes project[\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e]. A set of 2612 markers were used to differentiate sister species \u003cem\u003eAn. gambiae\u003c/em\u003e/\u003cem\u003eAn. coluzzii\u003c/em\u003e from the more divergent \u003cem\u003eAn. arabiensis\u003c/em\u003e, identified when the fraction of arabiensis-like alleles was \u0026gt;\u0026thinsp;0.6. A set of 700 AIMs were used to differentiate \u003cem\u003eAn. gambiae\u003c/em\u003e from \u003cem\u003eAn. coluzzii\u003c/em\u003e with samples scored as \u003cem\u003eAn. gambiae\u003c/em\u003e when the fraction of coluzzii-like calls was \u0026lt;\u0026thinsp;0.12 and \u003cem\u003eAn. coluzzii\u003c/em\u003e where this fraction was \u0026gt;\u0026thinsp;0.9. Individuals in-between these fractions represent other taxa. The taxonomic status of individuals was then confirmed with both PCA and an unrooted Neighbour-joining tree using 100,000 biallelic SNPs evenly spread across chromosome three because this region is unaffected by structural rearrangements such as inversions. Chosen SNPs had a minor allele frequency greater than 0.2% and no missing data. The Neighbour-Joining tree was constructed using city block distance.\u003c/p\u003e\n\u003ch3\u003ePopulation structure\u003c/h3\u003e\n\u003cp\u003eTo investigate population structure, Principal Component Analysis (PCA) dimensionality reduction and Neighbour-joining trees were constructed for \u003cem\u003eAn. coluzzii\u003c/em\u003e using the 3L chromosome as applied for taxonomic analysis. We also performed PCA focused on the allele counts at two inverted regions of the genome previously shown to segregate among \u003cem\u003eAn. coluzzii\u003c/em\u003e from different environments[\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e, \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e]. These regions are located on the 2L (2L:20,528,089\u0026thinsp;\u0026minus;\u0026thinsp;42,165,182) and 2R (2R:19,444,433\u0026thinsp;\u0026minus;\u0026thinsp;26,313,071) chromosomes. Inversion karyotype frequencies within each population cohort were assessed by typing individuals for inversion status using correlated tag SNPs as previously described[\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e]. We also assessed genomic differentiation among population cohorts using Hudson\u0026rsquo;s pairwise F\u003csub\u003eST\u003c/sub\u003e using the 3L chromosome arm[\u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e]. To explore the genetic diversity and demography of populations, informative summary statistics were then calculated including Nucleotide diversity (θπ), Watterson's theta (\u003cem\u003eθ\u003c/em\u003e\u003csub\u003e\u003cem\u003eW\u003c/em\u003e\u003c/sub\u003e) and Tajima\u0026rsquo;s D using cohorts with the malariagen_data python package. Statistics were only calculated for population cohorts with a minimum of ten individuals.\u003c/p\u003e\n\u003ch3\u003eInsecticide resistance\u003c/h3\u003e\n\u003cp\u003eTo investigate the presence of known target site mutations associated with insecticide resistance we calculated amino acid substitution frequencies at genomic sites of interest for each population cohort. These were based on the occurrence of non-synonymous SNPs at an appreciable frequency present at greater than 5%. Regions included the voltage-gated sodium channel (Vgsc; AGAP004707) as the target gene of pyrethroids, the glutathione S-transferase gene which confers resistance to DDT (\u003cem\u003eGste2\u003c/em\u003e; AGAP009194), the Resistance to dieldrin gene (\u003cem\u003eRdl\u003c/em\u003e; AGAP006028) and the organophosphate target gene, acetylcholinesterase (\u003cem\u003eAce1\u003c/em\u003e; AGAP001356). We also calculated the frequencies of copy number variants (CNVs) for genes associated with insecticide resistance and present at greater than 5% frequency. These included the Cytochrome P450 gene\u0026rsquo;s (AGAP002862-AGAP002870, AGAP000818, AGAP008212-AGAP008219), carboxylesterases (AGAP006228, AGAP006723-AGAP006728), \u003cem\u003eAce1\u003c/em\u003e (AGAP001356) and \u003cem\u003eGste2\u003c/em\u003e (AGAP009194).\u003c/p\u003e\u003cdiv id=\"Sec8\" class=\"Section2\"\u003e\u003ch2\u003eSelection scans\u003c/h2\u003e\u003cp\u003eTo identify novel regions of the genome under selection we utilised the H12 homozygosity statistic across windows of the genome[\u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e57\u003c/span\u003e] with statistical peaks representing either a hard or soft sweep. The statistic was calculated using an optimal size of 1500 windows. This was identified by plotting the distribution of H12 values across different window sizes and identifying when values fell below 0.1 for the 95th percentile.\u003c/p\u003e\u003c/div\u003e\n\u003ch3\u003eDiplotype clustering\u003c/h3\u003e\n\u003cp\u003eTo investigate which variants are associated with clusters of diplotypes under selection (i.e., regions of diploid genotypes), hierarchical clustering was performed using city block genetic distance and complete linkage using the malariagen_python package[\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e]. Observed amino acid substitutions and CNV variants were plotted onto the resulting dendrogram. Variants uniquely associated with clusters of diplotypes with low heterozygosity are candidates under selection.\u003c/p\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec11\" class=\"Section2\"\u003e\u003ch2\u003ePopulation Sampling\u003c/h2\u003e\u003cp\u003e\u003cem\u003eAnopheles\u003c/em\u003e gambiae \u003cem\u003esensu lato\u003c/em\u003e samples were collected from 2016 to 2018 across the arid savannah region of the Upper East in northern Ghana. Of 1663 samples sent for whole genome sequencing, 1473 individuals passed quality controls as defined by the \u003cem\u003eAnopheles gambiae\u003c/em\u003e 1000 Genomes (AG1000G) project[\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e, \u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e58\u003c/span\u003e]. From these samples, 1324 mosquitoes were identified as \u003cem\u003eAn. coluzzii\u003c/em\u003e through PCA and AIMs analysis generated during data curation[\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e, \u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e58\u003c/span\u003e] (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). In addition, we included 486 \u003cem\u003eAn. coluzzii\u003c/em\u003e genomes from a previous study in Ghana[\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e] (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). The sample sets included mosquitoes collected from 2012 to 2018 across five administrative districts located in the wet humid and deciduous forest region of South Ghana (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e,, Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e, Supplementary Table\u0026nbsp;1). For comparison across the West African region, we used data made publicly available in the AG1000G Phase 3 resource (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://malariagen.github.io/vector-data/ag3/ag3.0.html\u003c/span\u003e\u003cspan address=\"https://malariagen.github.io/vector-data/ag3/ag3.0.html\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003cspan type=\"Underline\" class=\"Underline\" name=\"Emphasis\"\u003e).\u003c/span\u003e The sample sets included \u003cem\u003eAn. coluzzii\u003c/em\u003e collected from 2004 to 2014, originating from the dry savannah regions of Burkina Faso (n\u0026thinsp;=\u0026thinsp;135) and Mali (n\u0026thinsp;=\u0026thinsp;90) or from the wet humid region of southern Cote d\u0026rsquo;Ivoire (n\u0026thinsp;=\u0026thinsp;80) (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e). Sample sets were classified for analysis as a population cohort based on the administrative district and year of collection. Samples sequenced during this study had a median coverage of 35X which generated 162,714,957 SNPs on alignment to the AgamP3 genome. Of these SNPs, 52,946,551 were biallelic and segregated within the samples.\u003c/p\u003e\u003cp\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e\u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e\u003cdiv class=\"CaptionContent\"\u003e\u003cp\u003eThe number of \u003cem\u003eAn. coluzzii\u003c/em\u003e collected from each sampling location in Ghana as part of the present study and details of other previously published datasets from West Africa included in analysis[\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e, \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e, \u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e58\u003c/span\u003e].\u003c/p\u003e\u003c/div\u003e\u003c/caption\u003e\u003ccolgroup cols=\"7\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e\u003cp\u003eCountry\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c2\"\u003e\u003cp\u003eLocation\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c3\"\u003e\u003cp\u003eLatitude\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c4\"\u003e\u003cp\u003eLongitude\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c5\"\u003e\u003cp\u003eYear\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c6\"\u003e\u003cp\u003eNo. samples\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c7\"\u003e\u003cp\u003eStudy ID\u003c/p\u003e\u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003eGhana\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eUpper East Region\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e10.85\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-1.11\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2018\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e563\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eThis study\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2017\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e557\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eThis study\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2016\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e204\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eThis study\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eAshanti Region\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e6.18\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-1.62\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2018\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e274\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eLucas \u003cem\u003eet al.\u003c/em\u003e 2023\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eCentral Region\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e6.01\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-1.87\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2018\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e148\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eLucas \u003cem\u003eet al.\u003c/em\u003e 2023\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e5.61\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-1.55\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2012\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e14\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eAG1000G Phase 3\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eEastern Region\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e6.09\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-0.26\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2012\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e1\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eAG1000G Phase 3\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eGreater Accra Region\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e5.67\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-0.22\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2012\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e14\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eAG1000G Phase 3\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eWestern Region\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e4.91\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-1.77\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2012\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e24\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eAG1000G Phase 3\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003eBurkina Faso\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eHauts-Bassins\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e11.151\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-4.235\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2012\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e82\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eAG1000G Phase 3\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2014\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e53\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eAG1000G Phase 3\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003eMali\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eKoulikouro\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e12.68\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-7.84\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2004\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e11\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eAG1000G Phase 3\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2014\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e27\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eAG1000G Phase 3\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSegou\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e13.2\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-6.13\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2004\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e25\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eAG1000G Phase 3\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSikasso\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e10.83\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-7.81\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2012\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e27\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eAG1000G Phase 3\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003eCote d'Ivoire\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eLagunes\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e5.898\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-4.823\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2012\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e80\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eAG1000G Phase 3\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/colgroup\u003e\u003c/table\u003e\u003c/div\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003cb\u003eGeographical population Structure of\u003c/b\u003e \u003cb\u003eAn. coluzzii\u003c/b\u003e \u003cb\u003ein Ghana\u003c/b\u003e\u003c/p\u003e\u003cp\u003ePopulation connectivity and consequent gene flow can influence the spread of adaptive alleles associated with insecticide resistance or vector competence for disease transmission. To first assess population connectivity of \u003cem\u003eAn. coluzzii\u003c/em\u003e across the different bioclimatic zones of Ghana, we implemented both principal components analysis (PCA) and neighbor-joining trees (NJT) using 100,000 single nucleotide polymorphisms (SNPs) spanning the 3L chromosome region (3L:15,000,000\u0026ndash;41,000,000), which is free from structural rearrangements that could bias inference. Both the principal component analysis (PCA) and neighbor-joining tree (NJT) presented two major clusters of \u003cem\u003eAn. coluzzii\u003c/em\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e2\u003c/span\u003eA and Fig. \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e). One cluster included samples from the northern arid savannah region of the Upper East admin region only. The second included all populations from southern Ghana including the wet and humid deciduous forest regions of Greater Accra, Ashanti, Central, Eastern and Western Ghana. Although some southern comparisons were from a different collection time to the northern cohorts, two southern cohorts from Ashanti and the Central region were sampled in the same year. Furthermore, since an older 2012 cohort from the Central region clusters with a later cohort from the same region in 2019, we can observe that the pattern of population structure between the north and south of Ghana has held across time. Genetic differentiation was also lower within northern and southern Ghana cohorts (F\u003csub\u003eST\u003c/sub\u003e 0.000\u0026ndash;0.005) but higher between these regions (F\u003csub\u003eST\u003c/sub\u003e 0.008\u0026ndash;0.013), suggesting restricted gene flow between the north and south of the country (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e3\u003c/span\u003e). Additionally, we observed that \u003cem\u003eAn. coluzzii\u003c/em\u003e from Greater Accra from 2012 diverged from the main cluster of individuals from southern Ghana on the PCA plot and appeared as a cluster of particularly closely related individuals in the Neighbor-Joining tree. This comparison included southern cohorts collected from the same time point. F\u003csub\u003eST\u003c/sub\u003e was also higher when comparing Greater Accra with the other southern population cohorts (F\u003csub\u003eST\u003c/sub\u003e 0.005). Although nucleotide diversity values did not differ substantially between any population cohort (Fig. \u003cspan refid=\"MOESM3\" class=\"InternalRef\"\u003eS3\u003c/span\u003e), the Greater Accra cohort had a slightly higher Tajima\u0026rsquo;s D indicating it may have experienced greater genetic drift. Overall findings suggest high connectivity among southern populations but restricted gene flow between the different bioclimatic regions of northern and southern Ghana.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eTo investigate further whether population structure in \u003cem\u003eAn. coluzzii\u003c/em\u003e could be associated with different bioclimatic zones more widely across West Africa, we extended our analysis to include the West African populations of Burkina Faso, Cote d\u0026rsquo;Ivoire and Mali available through the Vector Observatory (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.malariagen.net/vobs/\u003c/span\u003e\u003cspan address=\"https://www.malariagen.net/vobs/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). Once again, we observed two major PCA and NJT clusters (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e2\u003c/span\u003eB and Fig. \u003cspan refid=\"MOESM4\" class=\"InternalRef\"\u003eS4\u003c/span\u003e). One included northern Ghana and the arid savannah regions of Mali and Burkina Faso. Another included southern Ghana along with southern Cote d\u0026rsquo;Ivoire, which both experience similar wet humid climate conditions. In support of these findings, F\u003csub\u003eST\u003c/sub\u003e was lower between the northern cohorts from Ghana, Burkina Faso and Mali (F\u003csub\u003eST\u003c/sub\u003e 0.000-0.001) and between southern Ghana and Cote d\u0026rsquo;Ivoire (F\u003csub\u003eST\u003c/sub\u003e 0.005\u0026ndash;0.006) (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e3\u003c/span\u003e). In contrast, F\u003csub\u003eST\u003c/sub\u003e was higher between northern and southern comparisons (F\u003csub\u003eST\u003c/sub\u003e 0.008\u0026thinsp;\u0026minus;\u0026thinsp;0.0015). These results support the notion that gene flow between \u003cem\u003eAn. coluzzii\u003c/em\u003e is restricted between populations found in the arid savannah and wet humid forest environments of West Africa.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eTo further investigate gene flow among \u003cem\u003eAn. coluzzii\u003c/em\u003e in west and central Africa, we investigated two inverted regions of chromosome two previously associated with different climate conditions[\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e]. We used tagging SNPs, correlated with inversion status, to assess the frequency of the different inversion karyotypes in each population cohort[\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e]. One major 2La inversion karyotype was shared between \u003cem\u003eAn. coluzzii\u003c/em\u003e from the northern arid regions of Ghana, Burkina Faso and Mali, suggesting unrestricted gene flow among these population cohorts (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e4\u003c/span\u003eA). The inverted 2La inversion prevalent at 96\u0026ndash;100% frequency in the northern populations has been associated with arid environments and thermotolerance[\u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e59\u003c/span\u003e, \u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e60\u003c/span\u003e] (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). In contrast, inverted 2La was between 0\u0026ndash;16% frequency in the southern cohorts while the standard 2La\u0026thinsp;+\u0026thinsp;karyotype associated with mesic environments[\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e] was dominant at 33\u0026ndash;100% frequency. Analysis of the 2Rb inversion (2R:19,444,433\u0026thinsp;\u0026minus;\u0026thinsp;26,313,071) presented similar findings in that the southern and northern populations had different dominant inversion karyotypes (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e4\u003c/span\u003eB). However, findings differed in that all three 2Rb karyotypes were present in the northern regions while one form prevailed in the southern populations. Frequency analysis based on inversion tagging SNPs revealed that the dominant southern karyotype was the standard 2Rb\u0026thinsp;+\u0026thinsp;chromosomal form present at 93\u0026ndash;100% (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). Although all three karyotypes were present in the northern populations, the inverted 2Rb karyotype and heterozygote form reached higher frequencies (11\u0026ndash;63% and 32\u0026ndash;73%, respectively) than standard 2Rb+ (4\u0026ndash;46%) and were more prevalent on average in the north (37% and 45% for 2La and the heterozygote, respectively) than the south (0% and 4%, respectively). Findings are consistent with the observation that inverted 2Rb generally appears at higher frequencies in arid environments while standard 2Rb\u0026thinsp;+\u0026thinsp;is associated with mesic environments[\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e] and supports restricted gene flow between the northern Sahelian and southern forest regions of West Africa characterised by different bioclimatic zones.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e\u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e\u003cdiv class=\"CaptionContent\"\u003e\u003cp\u003ePercentage of 2La and 2Rb inversion karyotypes in population cohorts across the bioclimatic zones of West Africa using tagging SNPs for inversion karyotypes[\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e].\u003c/p\u003e\u003c/div\u003e\u003c/caption\u003e\u003ccolgroup cols=\"9\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c8\" colnum=\"8\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c9\" colnum=\"9\"\u003e\u003c/div\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e\u003cp\u003eCountry\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c2\"\u003e\u003cp\u003eLocation\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c3\"\u003e\u003cp\u003eYear\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colspan=\"2\" nameend=\"c5\" namest=\"c4\"\u003e\u003cp\u003eInversion karyotype %\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c6\"\u003e\u0026nbsp;\u003c/th\u003e\u003cth align=\"left\" colname=\"c7\"\u003e\u0026nbsp;\u003c/th\u003e\u003cth align=\"left\" colname=\"c8\"\u003e\u0026nbsp;\u003c/th\u003e\u003cth align=\"left\" colname=\"c9\"\u003e\u0026nbsp;\u003c/th\u003e\u003c/tr\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/th\u003e\u003cth align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/th\u003e\u003cth align=\"left\" colname=\"c3\"\u003e\u0026nbsp;\u003c/th\u003e\u003cth align=\"left\" colname=\"c4\"\u003e\u003cp\u003e2La+\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c5\"\u003e\u003cp\u003e2La+/2La\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c6\"\u003e\u003cp\u003e2La\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c7\"\u003e\u003cp\u003e2Rb+\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c8\"\u003e\u003cp\u003e2Rb+/2Rb\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c9\"\u003e\u003cp\u003e2Rb\u003c/p\u003e\u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e\u003cb\u003eCote d'Ivoire\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e\u003cb\u003eLagunes\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2012\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e41\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e45\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e14\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e93\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e8\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e\u003cb\u003eGhana\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e\u003cb\u003eWestern Region\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2012\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e33\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e58\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e8\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e96\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e1\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e\u003cb\u003eCentral Region\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2012\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e48\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e36\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e16\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e100\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e3\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2018\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e34\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e50\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e16\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e98\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e3\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e\u003cb\u003eGreater Accra Region\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2012\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e64\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e29\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e7\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e100\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e\u003cb\u003eEastern Region\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2012\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e100\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e100\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e\u003cb\u003eAshanti Region\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2018\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e46\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e45\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e9\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e96\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e10\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e\u003cb\u003eUpper East\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2016\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e1\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e99\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e12\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e50\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e39\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2017\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e1\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e99\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e8\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e40\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e52\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2018\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e100\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e7\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e34\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e60\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e\u003cb\u003eBurkina Faso\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e\u003cb\u003eHauts-Bassins\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2012\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e2\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e98\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e46\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e41\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e12\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2014\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e4\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e96\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e42\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e40\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e19\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e\u003cb\u003eMali\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e\u003cb\u003eKoulikoro\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2004\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e100\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e9\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e73\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e18\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2014\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e4\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e96\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e4\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e33\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e63\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e\u003cb\u003eSegou\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2004\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e4\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e96\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e8\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e32\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e60\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e\u003cb\u003eSikasso\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e\u003cb\u003e2012\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e4\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e96\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e\u003cp\u003e30\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c8\"\u003e\u003cp\u003e59\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c9\"\u003e\u003cp\u003e11\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/colgroup\u003e\u003c/table\u003e\u003c/div\u003e\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec12\" class=\"Section2\"\u003e\u003ch2\u003eResistance to Insecticides\u003c/h2\u003e\u003cp\u003eTo investigate geographical differences in insecticide resistance, we focused our analysis on the \u003cem\u003eAn. coluzzii\u003c/em\u003e data from Ghana in West Africa, which included country-wide sampling and the most recent timepoints for comparison. First, substitutions associated with target-site were investigated. Amino acid allele frequencies of three genes that encode for insecticide binding targets were computed: \u003cem\u003eVgsc\u003c/em\u003e (AGAP004707), \u003cem\u003eRdl\u003c/em\u003e (AGAP006028), and \u003cem\u003eAce1\u003c/em\u003e (AGAP001356). The \u003cem\u003ekdr\u003c/em\u003e allele \u003cem\u003eVgsc-L995F\u003c/em\u003e established as conferring insecticide resistance to pyrethroids[\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e, \u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e58\u003c/span\u003e] was identified in all populations but differed between the northern and southern populations with frequencies ranging from 57\u0026ndash;64% and 86\u0026ndash;92%, respectively (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e5\u003c/span\u003e). Similar frequencies of Vgsc-\u003cem\u003eL995F\u003c/em\u003e within the latter range were observed in the southern cohorts across the different years of collection ranging from 2012 to 2018. Additionally, we observed a double substitution, \u003cem\u003eVgsc-V402L\u003c/em\u003e and \u003cem\u003eI1527T\u003c/em\u003e, at 36\u0026ndash;43% frequency in northern Ghana and 8\u0026ndash;25% in the south, including directly comparable cohorts sampled from the different locations in 2018. The double substitution previously observed in \u003cem\u003eAn. coluzzii\u003c/em\u003e from Burkina Faso and Kenya[\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e, \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e] provides pyrethroid and DDT resistance at a reduced fitness cost compared to \u003cem\u003eL995F\u003c/em\u003e[\u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e62\u003c/span\u003e]. As a result, the substitution pair is expected to replace \u003cem\u003eL995F\u003c/em\u003e as the dominant \u003cem\u003eVgsc\u003c/em\u003e insecticide resistance mechanism. Similar to a study from Burkina Faso[\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e], we observed the \u003cem\u003eV402L\u003c/em\u003e and \u003cem\u003eI1527T\u003c/em\u003e substitution pair increasing in frequency. Frequencies increased 7% from 2016 to 2018 in northern Ghana while the frequency of \u003cem\u003eL995F\u003c/em\u003e reduced over the same period, advocating for a shift in the main \u003cem\u003eVgsc\u003c/em\u003e resistance allele in West Africa.\u003c/p\u003e\u003cp\u003eWe also observed both the allele pairs \u003cem\u003eA296G\u003c/em\u003e/\u003cem\u003eT345M and A296S\u003c/em\u003e/\u003cem\u003eT345S\u003c/em\u003e at the \u003cem\u003eRdl\u003c/em\u003e locus, which are known to confer resistance to organochlorines such as dieldrin[\u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e63\u003c/span\u003e, \u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e64\u003c/span\u003e] (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e5\u003c/span\u003e). However, we found a difference in their geographical distribution. We observed \u003cem\u003eA296G\u003c/em\u003e and \u003cem\u003eT345M\u003c/em\u003e in all cohorts from southern Ghana only, while \u003cem\u003eA296S\u003c/em\u003e and \u003cem\u003eT345S\u003c/em\u003e were only found in northern Ghana. This provides support for our finding that gene flow is restricted across Ghana, including at the 2La inversion region on which \u003cem\u003eRdl\u003c/em\u003e is located. To date, the \u003cem\u003eRdl-A296S\u003c/em\u003e/\u003cem\u003eT345S\u003c/em\u003e allele pair has only been reported in Burkina Faso[\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e] which shares a border with northern Ghana and supports our finding of population connectivity across the northern arid regions of West Africa.\u003c/p\u003e\u003cp\u003eTo investigate whether the different \u003cem\u003eRdl\u003c/em\u003e substitution pairs were associated with the chromosomal inversion karyotype on which the gene is located, we calculated the frequencies of \u003cem\u003eRdl\u003c/em\u003e substitutions for individuals with either the 2La or 2La\u0026thinsp;+\u0026thinsp;inversion. As expected, all individuals from the Upper East had the inverted 2La inversion karyotype and also the \u003cem\u003eA296S\u003c/em\u003e and \u003cem\u003eT345S\u003c/em\u003e substitution pair only found in this region (Fig. \u003cspan refid=\"MOESM5\" class=\"InternalRef\"\u003eS5\u003c/span\u003e). Individuals from southern Ghana with the standard 2La inversion had the \u003cem\u003eA296G\u003c/em\u003e and \u003cem\u003eT345M\u003c/em\u003e substitution pair associated with the region, with substitution frequencies for 2La population cohorts ranging from 16\u0026ndash;58%, but individuals with the inverted 2La inversion did not have either substitution pair. Interestingly, this was the case for individuals collected from the same cohort from Ashanti and the Central Region, suggesting an association of the \u003cem\u003eA296G\u003c/em\u003e and \u003cem\u003eT345M\u003c/em\u003e substitution with the 2La inversion karyotype.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eFinally, we observed the \u003cem\u003eG280S\u003c/em\u003e mutation and linked CNVs at the \u003cem\u003eAce1\u003c/em\u003e gene which are associated with resistance to organophosphates and carbamates[\u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e65\u003c/span\u003e, \u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e66\u003c/span\u003e] (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e5\u003c/span\u003e). These were present at \u0026gt;\u0026thinsp;5% in the population cohort from Greater Accra from 2012 only which formed a somewhat genetically differentiated group upon comparison with the other southern population cohorts on PCA analysis, including other southern populations sampled at the same time (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e2\u003c/span\u003eA).\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec13\" class=\"Section2\"\u003e\u003ch2\u003eMetabolic Resistance\u003c/h2\u003e\u003cp\u003eTo further investigate metabolic insecticide resistance in \u003cem\u003eAn. coluzzii\u003c/em\u003e across Ghana, we calculated the frequency of individuals within each population cohort with a copy number greater than two for genes associated with insecticide resistance, including cytochrome P450\u0026rsquo;s (AGAP002862-AGAP002870, AGAP000818, AGAP008212-AGAP008219), carboxylesterases (AGAP006228, AGAP006723-AGAP006728) and \u003cem\u003eGste2\u003c/em\u003e (AGAP009194). We also calculated amino acid substitution frequencies for the latter since the \u003cem\u003eGste2\u003c/em\u003e-\u003cem\u003eI114T\u003c/em\u003e mutation is known to increase the activity of \u003cem\u003eGste2\u003c/em\u003e[\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e]. We found CNVs at the cytochrome P450 cluster \u003cem\u003eCyp\u003c/em\u003e6 on chromosome two (Fig.\u0026nbsp;\u003cspan refid=\"Fig11\" class=\"InternalRef\"\u003e6\u003c/span\u003e). Amplifications at the \u003cem\u003eCyp6aa1\u003c/em\u003e (AGAP002862) and 2 (AGAP013128) regions were present across Ghana at 60\u0026ndash;84% frequency in the more recent population cohorts from 2018. This contrasts with the earlier sampling points from 2012, which presented frequencies between 0\u0026ndash;16%, suggesting that metabolic resistance has risen in the country. In addition, we observed duplications at the carboxylesterase cluster \u003cem\u003eCoeae2g-7g\u003c/em\u003e, which have been associated with resistance to pirimiphos methyl[\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e] at a higher frequency in northern Ghana (12\u0026ndash;28%) than southern Ghana (0\u0026ndash;6%), including cohorts collected in the same year (Fig.\u0026nbsp;\u003cspan refid=\"Fig11\" class=\"InternalRef\"\u003e6\u003c/span\u003e). \u003cem\u003eGste2\u003c/em\u003e CNV amplifications and the \u003cem\u003eI114\u003c/em\u003e substitution were also present at higher frequencies in northern Ghana (17\u0026ndash;23% and 49\u0026ndash;53%, respectively) compared to southern Ghana (4\u0026ndash;8% to 35\u0026ndash;38%, respectively) (Fig.\u0026nbsp;\u003cspan refid=\"Fig11\" class=\"InternalRef\"\u003e6\u003c/span\u003e). The latter had similar \u003cem\u003eI114\u003c/em\u003e frequencies across the different collection years. The exception was \u003cem\u003eAn. coluzzii\u003c/em\u003e from Greater Accra in 2012, which presented similar frequencies of both CNVs (21%) and the \u003cem\u003eGste2\u003c/em\u003e-\u003cem\u003eI114T\u003c/em\u003e substitution (46%) to the more recently sampled northern populations. Moreover, we observed CNVs at the cytochrome \u003cem\u003eCyp9k1\u003c/em\u003e at 21% frequency in Greater Accra, while frequencies were \u0026le;\u0026thinsp;5% in all other population cohorts (Fig.\u0026nbsp;\u003cspan refid=\"Fig11\" class=\"InternalRef\"\u003e6\u003c/span\u003e), including southern cohorts collected at the same time. Overall, our findings indicate that metabolic resistance varies across Ghana, with the northern populations in general more impacted than the southern populations, although the region of Greater Accra had particularly high CNV frequencies at insecticide resistance loci.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec14\" class=\"Section2\"\u003e\u003ch2\u003eSignals of selection\u003c/h2\u003e\u003cp\u003eTo identify signals of recent selection that may indicate novel insecticide resistance mechanisms, we calculated the H12 homozygosity statistic across windows of the genome and identified statistical peaks[\u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e67\u003c/span\u003e]. We observed a number of novel selection signals at loci which have not yet been functionally validated to confer insecticide resistance. We identified a selection peak over the \u003cem\u003eKeap1\u003c/em\u003e locus (AGAP003645: 2R:40,926,195\u0026thinsp;\u0026minus;\u0026thinsp;40,945,169) in all cohorts from southern Ghana (Fig.\u0026nbsp;\u003cspan refid=\"Fig13\" class=\"InternalRef\"\u003e7\u003c/span\u003e and Fig. \u003cspan refid=\"MOESM7\" class=\"InternalRef\"\u003eS7\u003c/span\u003e). \u003cem\u003eKeap1\u003c/em\u003e regulates the formation of the transcription factor Maf-S, known to trigger the expression of multiple metabolic resistance genes including cytochrome p450\u0026rsquo;s and glutathione S-transferases[\u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e68\u003c/span\u003e]. A selection signal has recently been observed at this gene in \u003cem\u003eAn. arabiensis\u003c/em\u003e from Kenya[\u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e69\u003c/span\u003e] and associated with deltamethrin treated survivors from Tanzania[\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e]. \u003cem\u003eKeap1\u003c/em\u003e therefore provides a good candidate for functional validation. In addition, we observed selection signals unique to the cohort from Greater Accra (Fig.\u0026nbsp;\u003cspan refid=\"Fig13\" class=\"InternalRef\"\u003e7\u003c/span\u003e). We observed a selection signal close to the cytochrome \u003cem\u003eCyp12f\u003c/em\u003e (AGAP008019 3R:4,324,183-4,326,568) which is differentially expressed in permethrin and DDT resistant strains of \u003cem\u003eAn. gambiae\u003c/em\u003e[\u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e70\u003c/span\u003e]. Furthermore, a selection signal was present near a UDP glucuronosyltransferase (AGAP028055 3R:2,836,386-2,838,097) which is a class of uridine diphosphate (UDP)-glycosyltransferase (UGT) detoxification enzymes involved in xenobiotic metabolism[\u003cspan citationid=\"CR71\" class=\"CitationRef\"\u003e71\u003c/span\u003e]. Findings of selection signals over novel genes with a possible link to insecticide resistance in only Greater Accra in addition to our observation of high frequencies of known target-site and metabolic resistance mechanisms in this cohort, suggests that this region is under particularly high insecticide resistance pressure.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eThe selection signals we observed at known insecticide resistance loci agreed with the CNV and substitution frequencies generated for each population cohort. For example, we observed a selection peak at the \u003cem\u003eCyp6\u003c/em\u003e gene cluster (2R), \u003cem\u003eVgsc\u003c/em\u003e (2L) and \u003cem\u003eGste2\u003c/em\u003e (3R) for all population cohorts for which we also observed high frequencies of CNVs and resistance associated substitutions (Fig. \u003cspan refid=\"MOESM6\" class=\"InternalRef\"\u003eS6\u003c/span\u003e-8). A signal at \u003cem\u003eRdl\u003c/em\u003e (2L) was apparent for the populations from the Upper East and Greater Accra, which both had relatively high frequencies of \u003cem\u003eRdl\u003c/em\u003e resistance associated substitutions. Additionally, Greater Accra had a signal at \u003cem\u003eAce1\u003c/em\u003e (2R) which was the only population cohort with CNVs and the G280S substitution at this locus. Interestingly, there was also a selection peak at \u003cem\u003eCoeae2f\u003c/em\u003e and \u003cem\u003eCoeae2g-7g\u003c/em\u003e on 2L for Greater Accra only although we only observed CNVs at appreciable frequencies in \u003cem\u003eAn. coluzzii\u003c/em\u003e from northern Ghana. It could be that the selection signal in Greater Accra is driven by another mechanism i.e., a substitution, although none have been identified as conferring resistance to date.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eAlthough functional studies have not yet associated substitutions at \u003cem\u003eKeap1\u003c/em\u003e with insecticide resistance, unique SNPs have been observed in haplotypes under selection for \u003cem\u003eAn. arabiensis\u003c/em\u003e from Kenya[\u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e69\u003c/span\u003e]. These included the substitution D780N and the stop gain mutation E762, which could result in loss of function and prevent the repression of detox gene expression in the absence of stress[\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e]. Therefore, we used diplotype clustering[\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e] to assess whether CNV duplications or SNPs were associated with haplotypes under selection in \u003cem\u003eAn. coluzzii\u003c/em\u003e from Ghana. We observed two haplotypes with low heterozygosity in southern Ghana for which we also observed a selection signal at the \u003cem\u003eKeap1\u003c/em\u003e locus (Fig. \u003cspan refid=\"MOESM9\" class=\"InternalRef\"\u003eS9\u003c/span\u003e). The two haplotypes were associated with a particular set of SNPs, including \u003cem\u003eV631I\u003c/em\u003e, \u003cem\u003eV816F\u003c/em\u003e and \u003cem\u003eV1001L\u003c/em\u003e and \u003cem\u003eG788R\u003c/em\u003e and \u003cem\u003eA943V\u003c/em\u003e, respectively. In particular, the latter three SNPs are at higher frequencies in southern Ghana (43\u0026ndash;52%) than northern Ghana (11\u0026ndash;19%) and present at similar frequencies, suggesting they are either linked or have a synergistic effect (Fig. \u003cspan refid=\"MOESM10\" class=\"InternalRef\"\u003eS10\u003c/span\u003e). These SNPs are different from those previously observed in East African \u003cem\u003eAn. arabiensis\u003c/em\u003e, but different molecular responses to insecticides are commonly observed across different taxa[\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e, \u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e67\u003c/span\u003e]. The role of the variants we observed in insecticide resistance could be determined through validation studies and also provide a focus for longitudinal studies tracking the frequency and geographical spread of candidate variants.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eWe also used diplotype clustering to further investigate genomic variation at the Carboxylesterase cluster \u003cem\u003eCoeae2g-7g\u003c/em\u003e since we observed a selection signal for \u003cem\u003eAn. coluzzii\u003c/em\u003e from Greater Accra despite the absence of CNVs. However, the dendrogram did not reveal any low heterozygosity clusters expected to result from selection on the region (Fig. \u003cspan refid=\"MOESM11\" class=\"InternalRef\"\u003eS11\u003c/span\u003e). Increased sampling of the Greater Accra Region may improve resolution to detect selected variants.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eUsing whole genome sequence variation data from West African \u003cem\u003eAn. coluzzii\u003c/em\u003e, we observed population structure and divergent insecticide resistance mechanisms driven by gene flow, environmental variation, and differences in selection pressure from population control. We found two genetically different groups of \u003cem\u003eAn. coluzzii\u003c/em\u003e based on population structure analysis of a large number of SNPs and two chromosomal inversion regions, both with similar levels of nucleotide diversity. These two populations may represent \u003cem\u003eAn. coluzzii\u003c/em\u003e subject to distinct ecological selection pressures in the arid savannah region of northern West Africa and the wet deciduous forest region of the south. Isolation by distance could explain the pattern of population structure we have observed. However, since our sites in the south were ~\u0026thinsp;500 km from those in the north and we did not have samples from intermediate sites, we were unable to determine whether this pattern reflects mosquito dispersal. However, the pattern of population structure between \u003cem\u003eAn. coluzzii\u003c/em\u003e from the different bioclimatic zones was present across West Africa where \u003cem\u003eAn. coluzzii\u003c/em\u003e are also expected to undergo long-range windborne dispersal to find suitable oviposition sites during the dry season[\u003cspan citationid=\"CR72\" class=\"CitationRef\"\u003e72\u003c/span\u003e]. Furthermore, there are no major geographical barriers expected to prevent mosquito dispersal across the investigated region. Therefore, another explanation for the population structure we observed is that reduced habitat suitability or local environmental adaptation restricts gene flow between northern (i.e Sahelian) and southern (ie. forest) West Africa. The former is supported by the geographical distribution of \u003cem\u003eAn. coluzzii\u003c/em\u003e which is bimodal because they predominate along the savannah and coastal regions[\u003cspan citationid=\"CR73\" class=\"CitationRef\"\u003e73\u003c/span\u003e], including in Ghana where aridity is increased in the north and southernmost coastal regions compared to the central zone[\u003cspan citationid=\"CR74\" class=\"CitationRef\"\u003e74\u003c/span\u003e]. However, it could also be that environmental adaptation leads to population specific mating or a fitness cost that restricts gene flow. Genomic evidence for environmental adaptation in \u003cem\u003eAn. coluzzii\u003c/em\u003e has been previously described based on the segregation of chromosome two inversions across West and Central Africa[\u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e59\u003c/span\u003e]. Similarly, we found that the inversion karyotypes of 2La and 2Rb segregated among \u003cem\u003eAn. coluzzii\u003c/em\u003e from northern and southern West Africa. Our findings support previous observations from Nigeria, Burkina Faso and Cameroon[\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e, \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e, \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e, \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e] that the inverted 2La and 2Rb predominates in arid regions while standard karyotypes are common in humid forest environments. Such inversions have the potential to preserve locally selected alleles by reducing recombination between heterozygotes and are thought to be major mechanisms driving speciation and the formation of \u003cem\u003eAnopheles\u003c/em\u003e ecotypes[\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e, \u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e75\u003c/span\u003e]. Ecological studies comparing \u003cem\u003eAn. coluzzii\u003c/em\u003e across West Africa are lacking and therefore further work is required to determine whether distinct phenotypic and/or behavioral differences are observed between cohorts from the different bioclimatic zones. However, a preference for mating with the local population has been observed for \u003cem\u003eAn. coluzzii\u003c/em\u003e sourced from different ecozones in Burkina Faso under experimental conditions, which are analogous to those in other West African countries[\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e]. Distinct phenotypic differences have also been observed in summer dormancy, i.e., aestivation, from contrasting arid and mesic environments[\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e, \u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e76\u003c/span\u003e]. Yet the question remains over whether aestivation in \u003cem\u003eAn. coluzzii\u003c/em\u003e has a genomic basis and can be linked to habitat use or geographical origin or whether it is a plastic response to dry season conditions, i.e., it is induced by the environment similar to the diapause response in other mosquitoes[\u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e76\u003c/span\u003e, \u003cspan citationid=\"CR77\" class=\"CitationRef\"\u003e77\u003c/span\u003e]. It is also unknown to what extent aestivating \u003cem\u003eAn. coluzzii\u003c/em\u003e can increase population structure between cohorts. Overall, further experimental work linking ecologically driven phenotypes to genomic variation is required to untangle whether \u003cem\u003eAn. coluzzi\u003c/em\u003ei from contrasting environments in West Africa represent ecotypes adapted to different environmental conditions, but our findings raise this possibility.\u003c/p\u003e\u003cp\u003eThe presence of ecologically distinct forms of \u003cem\u003eAn. coluzzii\u003c/em\u003e would have consequences for population control and malaria transmission since they will exhibit different life history traits and different responses to insecticide use. Indeed, we found that restricted gene flow between \u003cem\u003eAn. coluzzii\u003c/em\u003e across the different bioclimatic zones of northern and southern Ghana has led to different insecticide resistance profiles. Findings support previous work which found F\u003csub\u003eST\u003c/sub\u003e outlier regions at target site and metabolic resistance genes on comparison of forest and non-forest populations in Ghana, although comparisons did not include sites from northern Ghana[\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e]. Overall, we found that \u003cem\u003eAn. coluzzii\u003c/em\u003e from northern Ghana exhibited greater evidence for both target site and metabolic resistance than the southern populations, excepting the Greater Accra region. Along with higher frequencies of the known \u003cem\u003eVgsc\u003c/em\u003e-L995F substitution associated with pyrethroid resistance, northern Ghana also had the double substitution \u003cem\u003eVgsc\u003c/em\u003e-V402L and I1527T previously reported from nearby connected Burkina Faso[\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e]. Since this substitution pair provides pyrethroid resistance at a lower fitness cost than L995F[\u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e62\u003c/span\u003e], we expect it to rapidly spread across northern West Africa where we observed unrestricted gene flow. In addition, we observed higher frequencies of carboxylesterase CNVs and GSTE substitutions and CNVs in northern Ghana compared to the south, associated with multiple classes of insecticides including pirimiphos-methyl[\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e]. Although bed net use is widespread across the country, northern Ghana has a high malaria burden and has been additionally targeted with high coverage of long-lasting insecticide treated nets and IRS campaigns with propoxur between 2012 and 2016 [\u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e78\u003c/span\u003e, \u003cspan citationid=\"CR79\" class=\"CitationRef\"\u003e79\u003c/span\u003e]. Despite a lower prevalence of molecular insecticide resistance, \u003cem\u003eAn. coluzzii\u003c/em\u003e from southern Ghana are still impacted mainly through target site resistance at \u003cem\u003eVgsc\u003c/em\u003e-L955F. In addition, we found that the populations in southern Ghana were uniquely impacted by a selection signal at \u003cem\u003eKeap1\u003c/em\u003e. This gene is also under selection in East Africa and has the potential to impact the expression of detoxification genes to metabolise insecticides[\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e, \u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e69\u003c/span\u003e]. For example, a reduction in the regulatory action of \u003cem\u003eKeap1\u003c/em\u003e decreases mortality to organophosphates but increases mortality to pyrethroids and DDT[\u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e68\u003c/span\u003e]. Although it is unknown whether substitutions at the gene impact on insecticide resistance, we have identified several SNPs including V631I, V816F and V1001L and G788R and A943V which can be targeted for experimental work, including testing whether these variants have a synergistic effect. Although we also observed differences in \u003cem\u003eRdl\u003c/em\u003e substitutions between different bioclimatic zones and inversion karyotypes, the source of selection on this locus is unclear since dieldrin is no longer used in vector control. Possibly dieldrin use is ongoing, i.e., in agriculture, or the locus provides cross-resistance to another used insecticide.\u003c/p\u003e\u003cp\u003eWe found that \u003cem\u003eAn. coluzzii\u003c/em\u003e from Greater Accra on the southern coast of Ghana were somewhat genetically divergent from other southern populations based on PCA and F\u003csub\u003eST\u003c/sub\u003e population differentiation statistics. In addition, they had unique molecular insecticide resistance mechanisms and novel selection signals on chromosome three. For example, we found appreciable frequencies of \u003cem\u003eCyp9k1\u003c/em\u003e CNVs, the \u003cem\u003eAce1\u003c/em\u003e-G280S substitution and \u003cem\u003eAce1\u003c/em\u003e CNVs at this locus in Greater Accra only, suggesting resistance to multiple classes of insecticides. In addition, we found novel signals of selection at the cytochrome \u003cem\u003eCyp12f\u003c/em\u003e linked to permethrin and DDT resistance in \u003cem\u003eAn. gambiae\u003c/em\u003e[\u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e70\u003c/span\u003e] and a UDP glucuronosyltransferase analogous to a UGT gene that confers resistance to the pyrethroid lambda-cyhalthrin and the neonicotinoid imidacloprid in fruitflys[\u003cspan citationid=\"CR71\" class=\"CitationRef\"\u003e71\u003c/span\u003e]. Findings suggest that \u003cem\u003eAn. coluzzii\u003c/em\u003e are under intense selection pressure from insecticides in Greater Accra. This notion is supported by bioassay data, which suggested particularly high resistance to pyrethroids and carbamates in \u003cem\u003eAn. gambiae\u003c/em\u003e from the region[\u003cspan citationid=\"CR80\" class=\"CitationRef\"\u003e80\u003c/span\u003e]. However, the findings of phenotypic resistance are not directly comparable with the genomic data since our sampling was from 2012 while the bioassay experiment was conducted in 2017. The findings are somewhat surprising given that malaria incidence and bed net use in Greater Accra is comparatively low to elsewhere in Ghana[\u003cspan citationid=\"CR81\" class=\"CitationRef\"\u003e81\u003c/span\u003e, \u003cspan citationid=\"CR82\" class=\"CitationRef\"\u003e82\u003c/span\u003e]. However, Greater Accra is a populated urban area where there could be substantial use of household insecticide sprays for personal protection[\u003cspan citationid=\"CR83\" class=\"CitationRef\"\u003e83\u003c/span\u003e]. Furthermore, mosquito populations encounter increased pollution in urban environments, including habitats contaminated with insecticides from urban agriculture[\u003cspan citationid=\"CR73\" class=\"CitationRef\"\u003e73\u003c/span\u003e, \u003cspan citationid=\"CR80\" class=\"CitationRef\"\u003e80\u003c/span\u003e]. For example, pollution from urban agriculture can result in similar levels of resistance when compared to populations from rural and cultivated areas[\u003cspan additionalcitationids=\"CR85\" citationid=\"CR84\" class=\"CitationRef\"\u003e84\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR86\" class=\"CitationRef\"\u003e86\u003c/span\u003e]. Further study is required to identify the drivers of selection pressure within the city of Greater Accra. Even so, it is concerning that we have observed multiple and unique resistance mechanisms likely to reduce the efficacy of malaria control. Since we have found unrestricted gene flow across southern Ghana, any novel insecticide resistance mechanism has the potential to spread. However, this will also depend on selection pressure across the country, since resistance is often accompanied by a fitness cost[\u003cspan citationid=\"CR87\" class=\"CitationRef\"\u003e87\u003c/span\u003e, \u003cspan citationid=\"CR88\" class=\"CitationRef\"\u003e88\u003c/span\u003e]. Experimental work through bioassays and/or association studies are required to confirm that the novel regions under selection confer a resistance mechanism and to pinpoint the genomic variation underlying the phenotype.\u003c/p\u003e\u003cp\u003eOur results demonstrate that gene flow among vector populations is important in influencing the distribution of insecticide resistance mechanisms. A full understanding of mosquito population structure, both regional and large-scale, is required to predict the success of gene drive technologies and how novel insecticide resistance mechanisms will spread when they arise in response to population control measures. This includes a greater understanding of how environmental conditions influence mosquito dispersal and connectivity, and how this may alter with climate change. A strong understanding of population connectivity is particularly important given the introduction of new technologies introduced to combat the rise in metabolic resistance, i.e., dual active ingredient (AI) nets. New technologies impose novel selection pressures and may be quickly challenged given \u003cem\u003eAnopheles\u003c/em\u003e propensity for a rapid evolutionary response[\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e, \u003cspan citationid=\"CR89\" class=\"CitationRef\"\u003e89\u003c/span\u003e, \u003cspan citationid=\"CR90\" class=\"CitationRef\"\u003e90\u003c/span\u003e]. Any novel insecticide resistance mechanism will need to be quickly managed to maintain efficacy on their introduction. In tandem, large-scale routine monitoring of the temporal and geographical distribution of molecular insecticide resistance mechanisms across West Africa will be essential for a targeted defense.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe MalariaGEN Vector Observatory is supported by multiple institutes and funders. The Wellcome Sanger Institute’s participation was supported by funding from Wellcome \u0026nbsp;(220540/Z/20/A, 'Wellcome Sanger Institute Quinquennial Review 2021-2026') and the Bill \u0026amp; Melinda Gates Foundation (INV-001927 and INV-068808). The Liverpool School of Tropical Medicine's participation was supported by the National Institute of Allergy and Infectious Diseases ([NIAID] R01-AI116811), with additional support from the Medical Research Council (MR/P02520X/1). The latter grant is a UK-funded award and is part of the EDCTP2 programme supported by the European Union. Martin Donnelly is supported by a Royal Society Wolfson Fellowship (RSWF\\FT\\180003). The Pan-African Mosquito Control Association’s participation was funded by the Bill and Melinda Gates Foundation (INV-031595). Lucas N.\u0026nbsp;Amenga-Etego is supported by the Bill \u0026amp; Melinda Gates Foundation (INV-050873) and the National Institute of Health and Care Research, UK (grant number: NIHR134717).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAM, CC, GA and LAE conceptualised and designed the study, interpreted the data and assisted in drafting the manuscript. EKA, ID, CMM, SB, VAA, CA, KLM, \u0026nbsp; designed the study and conducted sample collection, processing, and data collection. EKA, AHK and KLB \u0026nbsp;conducted the data analysis and interpretation and assisted in drafting the manuscript. All the authors read, reviewed and approved this manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no competing interests.\u003c/p\u003e\n\u003cp\u003eNo chatbots or artificial intelligence tools were used in any of these studies.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMaterials and correspondence\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eEnock Kofi Amoako\u003cstrong\u003e;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eLucas Amenga Etego;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData availability\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe sequences of the samples identified in this study were submitted to the European Nucleotide Archive (ENA) (Project: PRJEB2141, accessions ERR2656751-ERR9796298).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis study was supported by the MalariaGEN Vector Observatory which is an international collaboration working to build capacity for malaria vector genomic research and surveillance, and involves contributions by the following institutions and teams. Wellcome Sanger Institute: Lee Hart, Kelly L. Bennett, Anastasia Hernandez-Koutoucheva, Jon Brenas, Menelaos Ioannidis, Chris Clarkson, Alistair Miles, Julia Jeans, Paballo Chauke, Victoria Simpson, Eleanor Drury, Osama Mayet, Sónia Gonçalves, Katherine Figueroa, Tom Maddison, Kevin Howe, Mara Lawniczak; Liverpool School of Tropical Medicine: Eric Lucas, Sanjay Nagi, Martin Donnelly; Broad Institute of Harvard and MIT: Jessica Way, George Grant; Pan-African Mosquito Control Association: Jane Mwangi, Edward Lukyamuzi, Sonia Barasa, Ibra Lujumba, Elijah Juma. The authors would like to thank the staff of the Wellcome Sanger Genomic Surveillance unit and the Wellcome Sanger Institute Sample Logistics, Sequencing and Informatics facilities for their contributions.\u003c/p\u003e\n\u003cp\u003eThe MalariaGEN Vector Observatory is supported by multiple institutes and funders. The Wellcome Sanger Institute’s participation was supported by funding from Wellcome \u0026nbsp;(220540/Z/20/A, 'Wellcome Sanger Institute Quinquennial Review 2021-2026') and the Bill \u0026amp; Melinda Gates Foundation (INV-001927 and INV-068808). The Liverpool School of Tropical Medicine's participation was supported by the National Institute of Allergy and Infectious Diseases ([NIAID] R01-AI116811), with additional support from the Medical Research Council (MR/P02520X/1). The latter grant is a UK-funded award and is part of the EDCTP2 programme supported by the European Union. Martin Donnelly is supported by a Royal Society Wolfson Fellowship (RSWF\\FT\\180003). The Pan-African Mosquito Control Association’s participation was funded by the Bill and Melinda Gates Foundation (INV-031595).\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eMwima R, Hui T-YJ, Nanteza A, Burt A, Kayondo JK. Potential persistence mechanisms of the major Anopheles gambiae species complex malaria vectors in sub-Saharan Africa: a narrative review. Malar J. 2023;22:336.\u003c/li\u003e\n\u003cli\u003eDe Mee\u0026ucirc;s T, Bouyer J, Ravel S, Solano P. Ecotype Evolution in Glossina palpalis Subspecies, Major Vectors of Sleeping Sickness. PLoS Negl Trop Dis. 2015;9:e0003497-.\u003c/li\u003e\n\u003cli\u003eSmall ST, Costantini C, Sagnon N, Guelbeogo MW, Emrich SJ, Kern AD, et al. Standing genetic variation and chromosome differences drove rapid ecotype formation in a major malaria mosquito. Proceedings of the National Academy of Sciences. 2023;120:e2219835120.\u003c/li\u003e\n\u003cli\u003eBhatt 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:207\u0026ndash;11. \u003c/li\u003e\n\u003cli\u003ede Souza D, Kelly-Hope L, Lawson B, Wilson M, Boakye D. Environmental Factors Associated with the Distribution of Anopheles gambiae s.s in Ghana; an Important Vector of Lymphatic Filariasis and Malaria. PLoS One. 2010;5:e9927.\u003c/li\u003e\n\u003cli\u003eDonkor E, Kelly M, Eliason C, Amotoh C, Gray DJ, Clements ACA, et al. A Bayesian spatio-temporal analysis of malaria in the Greater Accra region of Ghana from 2015 to 2019. Int J Environ Res Public Health. 2021;18:6080.\u003c/li\u003e\n\u003cli\u003eHinne IA, Attah SK, Mensah BA, Forson AO, Afrane YA. Larval habitat diversity and Anopheles mosquito species distribution in different ecological zones in Ghana. Parasit Vectors. 2021;14:193. \u003c/li\u003e\n\u003cli\u003eGhana Statistical Service. Ghana Multiple Indicator Cluster Survey with an Enhanced Malaria Module and Biomarker. 2011.\u003c/li\u003e\n\u003cli\u003eNignan C, Poda BS, Sawadogo SP, Ma\u0026iuml;ga H, Dabir\u0026eacute; KR, Gnankine O, et al. Local adaptation and colonization are potential factors affecting sexual competitiveness and mating choice in Anopheles coluzzii populations. Sci Rep. 2022;12:636. \u003c/li\u003e\n\u003cli\u003eLehmann T, Weetman D, Huestis DL, Yaro AS, Kassogue Y, Diallo M, et al. Tracing the origin of the early wet-season Anopheles coluzzii in the Sahel. Evol Appl. 2017;10:704\u0026ndash;17. \u003c/li\u003e\n\u003cli\u003eYaro AS, Traor\u0026eacute; AI, Huestis DL, Adamou A, Timbin\u0026eacute; S, Kassogu\u0026eacute; Y, et al. Dry season reproductive depression of Anopheles gambiae in the Sahel. J Insect Physiol. 2012;58:1050\u0026ndash;9. \u003c/li\u003e\n\u003cli\u003eColuzzi M, Sabatini A, Petrarca V, Di Deco MA. Chromosomal differentiation and adaptation to human environments in the Anopheles gambiae complex. Trans R Soc Trop Med Hyg. 1979;73:483\u0026ndash;97. \u003c/li\u003e\n\u003cli\u003eCostantini C, Ayala D, Guelbeogo WM, Pombi M, Some CY, Bassole IHN, et al. Living at the edge: biogeographic patterns of habitat segregation conform to speciation by niche expansion in Anopheles gambiae. BMC Ecol. 2009;9:16. \u003c/li\u003e\n\u003cli\u003eSimard F, Ayala D, Kamdem GC, Pombi M, Etouna J, Ose K, et al. Ecological niche partitioning between Anopheles gambiae molecular forms in Cameroon: the ecological side of speciation. BMC Ecol. 2009;9:17. \u003c/li\u003e\n\u003cli\u003eAdeogun AO, Popoola KOK, Brooke BD, Olakiigbe AK, Awolola ST. Polymorphic inversion 2La frequencies associated with ecotypes in populations of Anopheles coluzzii from Southwest Nigeria. Sci Afr. 2021;12:e00746. \u003c/li\u003e\n\u003cli\u003eAyala D, Acevedo P, Pombi M, Dia I, Boccolini D, Costantini C, et al. Chromosome inversions and ecological plasticity in the main African malaria mosquitoes. Evolution (N Y). 2017;71:686\u0026ndash;701. \u003c/li\u003e\n\u003cli\u003eSimard F, Ayala D, Kamdem GC, Pombi M, Etouna J, Ose K, et al. Ecological niche partitioning between Anopheles gambiae molecular forms in Cameroon: the ecological side of speciation. BMC Ecol. 2009;9:1\u0026ndash;24. \u003c/li\u003e\n\u003cli\u003eColuzzi M, Sabatini A, Petrarca V, Di Deco MA. Chromosomal differentiation and adaptation to human environments in the Anopheles gambiae complex. Trans R Soc Trop Med Hyg. 1979;73:483\u0026ndash;97. \u003c/li\u003e\n\u003cli\u003eLove RR, Pombi M, Guelbeogo MW, Campbell NR, Stephens MT, Dabire RK, et al. Inversion Genotyping in the Anopheles gambiae Complex Using High-Throughput Array and Sequencing Platforms. G3 Genes|Genomes|Genetics. 2020;10:3299\u0026ndash;307. \u003c/li\u003e\n\u003cli\u003ePetrarca V, Beier JC. Intraspecific chromosomal polymorphism in the Anopheles gambiae complex as a factor affecting malaria transmission in the Kisumu area of Kenya. Am J Trop Med Hyg. 1992;46:229\u0026ndash;37.\u003c/li\u003e\n\u003cli\u003eAdeogun AO, Brooke BD, Olayanju DR, Adegbehingbe K, Oyeniyi TA, Olakiigbe AK, et al. Test for association between dieldrin resistance and 2La inversion polymorphism in Anopheles coluzzii from Lagos, Nigeria. 2019. Tropical biomedicine, 36;3:587\u0026ndash;593.\u003c/li\u003e\n\u003cli\u003eBrooke BD, Hunt RH, Coetzee M. Resistance to dieldrin\u0026emsp;+\u0026emsp;fipronil assorts with chromosome inversion 2La in the malaria vector Anopheles gambiae. Med Vet Entomol. 2000;14:190\u0026ndash;4. \u003c/li\u003e\n\u003cli\u003eClarkson CS, Weetman D, Essandoh J, Yawson AE, Maslen G, Manske M, et al. Adaptive introgression between Anopheles sibling species eliminates a major genomic island but not reproductive isolation. Nat Commun. 2014;5:4248.\u003c/li\u003e\n\u003cli\u003eRanson 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\u0026ndash;96.\u003c/li\u003e\n\u003cli\u003eHemingway J, Ranson H, Magill A, Kolaczinski J, Fornadel C, Gimnig J, et al. Averting a malaria disaster: will insecticide resistance derail malaria control? The Lancet. 2016;387:1785\u0026ndash;8.\u003c/li\u003e\n\u003cli\u003eDavies TGE, Field LM, Usherwood PNR, Williamson MS. A comparative study of voltage-gated sodium channels in the Insecta: implications for pyrethroid resistance in Anopheline and other Neopteran species. Insect Mol Biol. 2007;16:361\u0026ndash;75. \u003c/li\u003e\n\u003cli\u003eMartinez-Torres D, Chandre F, Williamson MS, Darriet F, Berg\u0026eacute; JB, Devonshire AL, et al. Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998;7:179\u0026ndash;84. \u003c/li\u003e\n\u003cli\u003eEssandoh J, Yawson AE, Weetman D. Acetylcholinesterase (Ace-1) target site mutation 119S is strongly diagnostic of carbamate and organophosphate resistance in Anopheles gambiae ss and Anopheles coluzzii across southern Ghana. Malar J. 2013;12:404.\u003c/li\u003e\n\u003cli\u003eWeill M, Malcolm C, Chandre F, Mogensen K, Berthomieu A, Marquine M, et al. The unique mutation in ace‐1 giving high insecticide resistance is easily detectable in mosquito vectors. Insect Mol Biol. 2004;13:1\u0026ndash;7.\u003c/li\u003e\n\u003cli\u003eHancock PA, Ochomo E, Messenger LA. Genetic surveillance of insecticide resistance in African Anopheles populations to inform malaria vector control. Trends Parasitol. 2024;40:604\u0026ndash;18. \u003c/li\u003e\n\u003cli\u003eDu W, Awolola TS, Howell P, Koekemoer LL, Brooke BD, Benedict MQ, et al. Independent mutations in the Rdl locus confer dieldrin resistance to Anopheles gambiae and An. arabiensis. Insect Mol Biol. 2005;14:179\u0026ndash;83.\u003c/li\u003e\n\u003cli\u003eEdi C V, Djogb\u0026eacute;nou L, Jenkins AM, Regna K, Muskavitch MAT, Poupardin R, et al. CYP6 P450 Enzymes and ACE-1 Duplication Produce Extreme and Multiple Insecticide Resistance in the Malaria Mosquito Anopheles gambiae. PLoS Genet. 2014;10:e1004236.\u003c/li\u003e\n\u003cli\u003eIbrahim SS, Riveron JM, Stott R, Irving H, Wondji CS. The cytochrome P450 CYP6P4 is responsible for the high pyrethroid resistance in knockdown resistance-free Anopheles arabiensis. Insect Biochem Mol Biol. 2016;68:23\u0026ndash;32. \u003c/li\u003e\n\u003cli\u003eNagi SC, Lucas ER, Egyir-Yawson A, Essandoh J, Dadzie S, Chabi J, et al. Parallel Evolution in Mosquito Vectors\u0026mdash;A Duplicated Esterase Locus is Associated With Resistance to Pirimiphos-methyl in Anopheles gambiae. Mol Biol Evol. 2024;41:msae140.\u003c/li\u003e\n\u003cli\u003eLucas ER, Nagi SC, Kabula B, Batengana B, Kisinza W, Egyir-Yawson A, et al. Copy number variants underlie major selective sweeps in insecticide resistance genes in Anopheles arabiensis. PLoS Biol. 2024;22:e3002898.\u003c/li\u003e\n\u003cli\u003eRiveron JM, Yunta C, Ibrahim SS, Djouaka R, Irving H, Menze BD, et al. A single mutation in the GSTe2 gene allows tracking of metabolically based insecticide resistance in a major malaria vector. Genome Biol. 2014;15:R27. \u003c/li\u003e\n\u003cli\u003ePu J, Chung H. New and emerging mechanisms of insecticide resistance. Curr Opin Insect Sci. 2024;63:101184. \u003c/li\u003e\n\u003cli\u003eDabir\u0026eacute; RK, Namountougou M, Sawadogo SP, Yaro LB, To\u0026eacute; HK, Ouari A, et al. Population dynamics of Anopheles gambiae s.l. in Bobo-Dioulasso city: bionomics, infection rate and susceptibility to insecticides. Parasit Vectors. 2012;5:127.\u003c/li\u003e\n\u003cli\u003ePerugini E, Pichler V, Guelbeogo WM, Micocci M, Poggi C, Manzi S, et al. Longitudinal survey of insecticide resistance in a village of central region of Burkina Faso reveals co-occurrence of 1014F, 1014S and 402L mutations in Anopheles coluzzii and Anopheles arabiensis. Malar J. 2024;23:250. \u003c/li\u003e\n\u003cli\u003eLucas ER, Nagi SC, Egyir-Yawson A, Essandoh J, Dadzie S, Chabi J, et al. Genome-wide association studies reveal novel loci associated with pyrethroid and organophosphate resistance in Anopheles gambiae and Anopheles coluzzii. Nat Commun. 2023;14:4946. \u003c/li\u003e\n\u003cli\u003eIbrahim SS, Muhammad A, Hearn J, Weedall GD, Nagi SC, Mukhtar MM, et al. Molecular drivers of insecticide resistance in the Sahelo-Sudanian populations of a major malaria vector Anopheles coluzzii. BMC Biol. 2023;21:125. \u003c/li\u003e\n\u003cli\u003eKamau L, Bennett KL, Ochomo E, Herren J, Agumba S, Otieno S, et al. The Anopheles coluzzii range extends into Kenya: detection, insecticide resistance profiles and population genetic structure in relation to conspecific populations in West and Central Africa. Malar J. 2024;23:122. \u003c/li\u003e\n\u003cli\u003eDennis TPW, Essandoh J, Mable BK, Viana MS, Yawson AE, Weetman David. Signatures of adaptation at key insecticide resistance loci in Anopheles gambiae in Southern Ghana revealed by reduced-coverage WGS. Sci Rep. 2024;14:8650. \u003c/li\u003e\n\u003cli\u003eConsortium A gambiae 1000G. Genetic diversity of the African malaria vector Anopheles gambiae. Nature. 2017;552:96.\u003c/li\u003e\n\u003cli\u003eClarkson CS, Miles A, Harding NJ, Lucas ER, Battey CJ, Amaya-Romero JE, et al. Genome variation and population structure among 1142 mosquitoes of the African malaria vector species Anopheles gambiae and Anopheles coluzzii. Genome Res. 2020;30:1533\u0026ndash;46.\u003c/li\u003e\n\u003cli\u003eIngham VA, Tennessen JA, Lucas ER, Elg S, Yates HC, Carson J, et al. Integration of whole genome sequencing and transcriptomics reveals a complex picture of the reestablishment of insecticide resistance in the major malaria vector Anopheles coluzzii. PLoS Genet. 2021;17:e1009970.\u003c/li\u003e\n\u003cli\u003eCoetzee M. Key to the females of Afrotropical Anopheles mosquitoes (Diptera: Culicidae). Malar J. 2020;19:70.\u003c/li\u003e\n\u003cli\u003eLi H, Durbin R. Fast and accurate short read alignment with Burrows\u0026ndash;Wheeler transform. Bioinformatics. 2009;25:1754\u0026ndash;60. \u003c/li\u003e\n\u003cli\u003eMcKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297\u0026ndash;303.\u003c/li\u003e\n\u003cli\u003eMartin M, Ebert P, Marschall T. Read-Based Phasing and Analysis of Phased Variants with WhatsHap. In: Peters BA, Drmanac R, editors. Haplotyping: Methods and Protocols. New York, NY: Springer US; 2023. p. 127\u0026ndash;38.\u003c/li\u003e\n\u003cli\u003eDelaneau O, Marchini J, McVean GA, Donnelly P, Lunter G, Marchini JL, et al. Integrating sequence and array data to create an improved 1000 Genomes Project haplotype reference panel. Nat Commun. 2014;5:3934.\u003c/li\u003e\n\u003cli\u003eLucas ER, Miles A, Harding NJ, Clarkson CS, Lawniczak MKN, Kwiatkowski DP, et al. Whole-genome sequencing reveals high complexity of copy number variation at insecticide resistance loci in malaria mosquitoes. Genome Res. 2019;29:1250\u0026ndash;61.\u003c/li\u003e\n\u003cli\u003eNeafsey DE, Waterhouse RM, Abai MR, Aganezov SS, Alekseyev MA, Allen JE, et al. Highly evolvable malaria vectors: the genomes of 16 Anopheles mosquitoes. Science (1979). 2015;347:1258522.\u003c/li\u003e\n\u003cli\u003eAyala D, Ullastres A, Gonz\u0026aacute;lez J. Adaptation through chromosomal inversions in Anopheles. Front Genet. 2014;5:129.\u003c/li\u003e\n\u003cli\u003eLove RR, Pombi M, Guelbeogo MW, Campbell NR, Stephens MT, Dabire RK, et al. Inversion Genotyping in the Anopheles gambiae Complex Using High-Throughput Array and Sequencing Platforms. G3 Genes|Genomes|Genetics. 2020;10:3299\u0026ndash;307. \u003c/li\u003e\n\u003cli\u003eHudson RR, Slatkin M, Maddison WP. Estimation of levels of gene flow from DNA sequence data. Genetics. 1992;132:583\u0026ndash;9. \u003c/li\u003e\n\u003cli\u003eGarud NR, Messer PW, Buzbas EO, Petrov DA. Recent Selective Sweeps in North American Drosophila melanogaster Show Signatures of Soft Sweeps. PLoS Genet. 2015;11:e1005004.\u003c/li\u003e\n\u003cli\u003eClarkson CS, Miles A, Harding NJ, O\u0026rsquo;Reilly AO, Weetman D, Kwiatkowski D, et al. The genetic architecture of target‐site resistance to pyrethroid insecticides in the African malaria vectors Anopheles gambiae and Anopheles coluzzii. Mol Ecol. 2021;30:5303\u0026ndash;17.\u003c/li\u003e\n\u003cli\u003eAyala D, Zhang S, Chateau M, Fouet C, Morlais I, Costantini C, et al. Association mapping desiccation resistance within chromosomal inversions in the African malaria vector Anopheles gambiae. Mol Ecol. 2019;28:1333\u0026ndash;42.\u003c/li\u003e\n\u003cli\u003eIbrahim SS, Mukhtar MM, Muhammad A, Wondji CS. 2La paracentric chromosomal inversion and overexpressed metabolic genes enhance thermotolerance and pyrethroid resistance in the major malaria vector Anopheles gambiae. Biology. 2021;10:518.\u003c/li\u003e\n\u003cli\u003eKientega M, Clarkson CS, Traor\u0026eacute; N, Hui T-YJ, O\u0026rsquo;Loughlin S, Millogo A, et al. Whole-genome sequencing of major malaria vectors reveals the evolution of new insecticide resistance variants in a longitudinal study in Burkina Faso. Malar J. 2023;:280. \u003c/li\u003e\n\u003cli\u003eWilliams J, Cowlishaw R, Sanou A, Ranson H, Grigoraki L. In vivo functional validation of the V402L voltage gated sodium channel mutation in the malaria vector An. gambiae. Pest Manag Sci. 2022;78:1155\u0026ndash;63.\u003c/li\u003e\n\u003cli\u003eDu W, Awolola TS, Howell P, Koekemoer LL, Brooke BD, Benedict MQ, et al. Independent mutations in the Rdl locus confer dieldrin resistance to Anopheles gambiae and An. arabiensis. Insect Mol Biol. 2005;14:179\u0026ndash;83. \u003c/li\u003e\n\u003cli\u003eWondji CS, Dabire RK, Tukur Z, Irving H, Djouaka R, Morgan JC. Identification and distribution of a GABA receptor mutation conferring dieldrin resistance in the malaria vector Anopheles funestus in Africa. Insect Biochem Mol Biol. 2011;41:484\u0026ndash;91. \u003c/li\u003e\n\u003cli\u003eLucas ER, Rockett KA, Lynd A, Essandoh J, Grisales N, Kemei B, et al. A high throughput multi-locus insecticide resistance marker panel for tracking resistance emergence and spread in Anopheles gambiae. Sci Rep. 2019;9:13335.\u003c/li\u003e\n\u003cli\u003eMitchell SN, Rigden DJ, Dowd AJ, Lu F, Wilding CS, Weetman D, et al. Metabolic and target-site mechanisms combine to confer strong DDT resistance in Anopheles gambiae. PLoS One. 2014;9:e92662.\u003c/li\u003e\n\u003cli\u003eMwinyi SH, Bennett KL, Nagi SC, Kabula B, Matowo J, Weetman D, et al. Genomic Analysis Reveals a New Cryptic Taxon Within the Anopheles gambiae Complex With a Distinct Insecticide Resistance Profile in the Coast of East Africa. Mol Ecol. 2025;:e17762.\u003c/li\u003e\n\u003cli\u003eIngham VA, Pignatelli P, Moore JD, Wagstaff S, Ranson H. The transcription factor Maf-S regulates metabolic resistance to insecticides in the malaria vector Anopheles gambiae. BMC Genomics. 2017;18:1\u0026ndash;11.\u003c/li\u003e\n\u003cli\u003ePolo B, Bennett KL, Barasa S, Brenas J, Agumba S, Mwangangi J, et al. Genomic surveillance reveals geographical heterogeneity and differences in known and novel insecticide resistance mechanisms in Anopheles arabiensis across Kenya. 2024. \u003c/li\u003e\n\u003cli\u003eDavid J-P, Strode C, Vontas J, Nikou D, Vaughan A, Pignatelli PM, et al. The Anopheles gambiae detoxification chip: a highly specific microarray to study metabolic-based insecticide resistance in malaria vectors. Proceedings of the National Academy of Sciences. 2005;102:4080\u0026ndash;4.\u003c/li\u003e\n\u003cli\u003eLogan RAE, M\u0026auml;urer JB, Wapler C, Ingham VA. Uridine diphosphate (UDP)-glycosyltransferases (UGTs) are associated with insecticide resistance in the major malaria vectors Anopheles gambiae s.l. and Anopheles funestus. Sci Rep. 2024;14:19821. \u003c/li\u003e\n\u003cli\u003eHuestis DL, Dao A, Diallo M, Sanogo ZL, Samake D, Yaro AS, et al. Windborne long-distance migration of malaria mosquitoes in the Sahel. Nature. 2019;574:404\u0026ndash;8. \u003c/li\u003e\n\u003cli\u003eTene‐Fossog B, Fotso‐Toguem YG, Amvongo‐Adjia N, Ranson H, Wondji CS. Temporal variation of high‐level pyrethroid resistance in the major malaria vector Anopheles gambiae sl in Yaound\u0026eacute;, Cameroon, is mediated by target‐site and metabolic resistance. Med Vet Entomol. 2022;36:247\u0026ndash;59.\u003c/li\u003e\n\u003cli\u003eKudom AA. Larval ecology of Anopheles coluzzii in Cape Coast, Ghana: water quality, nature of habitat and implication for larval control. Malar J. 2015;14:447. \u003c/li\u003e\n\u003cli\u003eAyala FJ, Coluzzi M. Chromosome speciation: Humans, Drosophila, and mosquitoes. Proceedings of the National Academy of Sciences. 2005;102 suppl_1:6535\u0026ndash;42. \u003c/li\u003e\n\u003cli\u003eHidalgo K, Siaussat D, Braman V, Dabir\u0026eacute; KR, Simard F, Mouline K, et al. Comparative physiological plasticity to desiccation in distinct populations of the malarial mosquito Anopheles coluzzii. Parasit Vectors. 2016;9:1\u0026ndash;13.\u003c/li\u003e\n\u003cli\u003eArmbruster PA. Photoperiodic diapause and the establishment of Aedes albopictus (Diptera: Culicidae) in North America. J Med Entomol. 2016;53:1013\u0026ndash;23.\u003c/li\u003e\n\u003cli\u003eGogue C, Wagman J, Tynuv K, Saibu A, Yihdego Y, Malm K, et al. An observational analysis of the impact of indoor residual spraying in Northern, Upper East, and Upper West Regions of Ghana: 2014 through 2017. Malar J. 2020;19:1\u0026ndash;13.\u003c/li\u003e\n\u003cli\u003eTiedje KE, Oduro AR, Bangre O, Amenga-Etego L, Dadzie SK, Appawu MA, et al. Indoor residual spraying with a non-pyrethroid insecticide reduces the reservoir of Plasmodium falciparum in a high-transmission area in northern Ghana. PLOS Global Public Health. 2022;2:e0000285.\u003c/li\u003e\n\u003cli\u003ePwalia R, Joannides J, Iddrisu A, Addae C, Acquah-Baidoo D, Obuobi D, et al. High insecticide resistance intensity of Anopheles gambiae (sl) and low efficacy of pyrethroid LLINs in Accra, Ghana. Parasit Vectors. 2019;12:1\u0026ndash;9.\u003c/li\u003e\n\u003cli\u003eKawaguchi K, Donkor E, Lal A, Kelly M, Wangdi K. Distribution and risk factors of malaria in the Greater Accra Region in Ghana. Int J Environ Res Public Health. 2022;19:12006.\u003c/li\u003e\n\u003cli\u003eAheto JMK, Menezes LJ, Takramah W, Cui L. Modelling spatiotemporal variation in under-five malaria risk in Ghana in 2016\u0026ndash;2021. Malar J. 2024;23:102.\u003c/li\u003e\n\u003cli\u003eSilva Martins WF, Reid E, Tomlinson S, Evans G, Gibson J, Guy A, et al. Improving the efficiency of aerosolized insecticide testing against mosquitoes. Sci Rep. 2023;13:6281. \u003c/li\u003e\n\u003cli\u003eChabi J, Eziefule MC, Pwalia R, Joannides J, Obuobi D, Amlalo G, et al. Impact of urban agriculture on the species distribution and insecticide resistance profile of Anopheles gambiae ss and Anopheles coluzzii in Accra Metropolis, Ghana. Advances in Entomology. 2018;6:198.\u003c/li\u003e\n\u003cli\u003eTchigossou G, Dossou C, Tepa-Yotto G, Koto M, Atoyebi SM, Tossou E, et al. Resistance to neonicotinoids is associated with metabolic detoxification mechanisms in Anopheles coluzzii from agricultural and urban sites in southern Benin. Frontiers in Tropical Diseases. 2024;5:1339811.\u003c/li\u003e\n\u003cli\u003eAntonio-Nkondjio C, Fossog BT, Ndo C, Djantio BM, Togouet SZ, Awono-Ambene P, et al. Anopheles gambiae distribution and insecticide resistance in the cities of Douala and Yaound\u0026eacute; (Cameroon): influence of urban agriculture and pollution. Malar J. 2011;10:154. \u003c/li\u003e\n\u003cli\u003eNkahe DL, Kopya E, Djiappi-Tchamen B, Toussile W, Sonhafouo-Chiana N, Kekeunou S, et al. Fitness cost of insecticide resistance on the life-traits of a Anopheles coluzzii population from the city of Yaound\u0026eacute;, Cameroon. Wellcome Open Res. 2020;5:171.\u003c/li\u003e\n\u003cli\u003eGul H, Gadratagi BG, G\u0026uuml;ncan A, Tyagi S, Ullah F, Desneux N, et al. Fitness costs of resistance to insecticides in insects. Front Physiol. 2023;14:1238111.\u003c/li\u003e\n\u003cli\u003eTo\u0026eacute; KH, N\u0026rsquo;Fal\u0026eacute; S, Dabir\u0026eacute; RK, Ranson H, Jones CM. The recent escalation in strength of pyrethroid resistance in Anopheles coluzzi in West Africa is linked to increased expression of multiple gene families. BMC Genomics. 2015;16:146. \u003c/li\u003e\n\u003cli\u003eNjoroge H, van\u0026rsquo;t Hof A, Oruni A, Pipini D, Nagi SC, Lynd A, et al. Identification of a rapidly-spreading triple mutant for high-level metabolic insecticide resistance in Anopheles gambiae provides a real-time molecular diagnostic for antimalarial intervention deployment. Mol Ecol. 2022;31:4307\u0026ndash;18. \u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"bmc-genomics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"gics","sideBox":"Learn more about [BMC Genomics](http://bmcgenomics.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/gics","title":"BMC Genomics","twitterHandle":"#BMCGenomics","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"","lastPublishedDoi":"10.21203/rs.3.rs-7878288/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-7878288/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eEnvironmental barriers influencing the movement of insect vectors can govern adaptive gene flow, including the dispersal of insecticide resistance mechanisms that compromise population control. We sought to understand population connectivity of the major malaria vector, \u003cem\u003eAn. coluzzii\u003c/em\u003e, across the different bioclimatic zones of West Africa using SNPs from whole genomes and inversion karyotypes previously associated with environmental adaptation. We identified restricted gene flow between populations from northern savannah and southern forested regions. Using Ghana as a case study, we found marked differences in insecticide resistance profiles across the different bioclimatic zones suggesting that population connectivity impacts on adaptive allele sharing. Greater evidence for target site pyrethroid and metabolic cross-resistance in the North reflects differences in insecticide use across the country. We also observed distinct resistance mechanisms in the coastal region of Greater Accra which may result from intense urban agricultural activity. Overall, findings suggest that environmental conditions restrict \u003cem\u003eAn. coluzzii\u003c/em\u003e gene flow to impact the geographical distribution of molecular insecticide resistance.\u003c/p\u003e","manuscriptTitle":"Genomic population structure and insecticide resistance mechanisms in the malaria vector An. coluzzii across contrasting bioclimatic zones in West Africa","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-11-05 06:40:24","doi":"10.21203/rs.3.rs-7878288/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2025-11-10T19:57:33+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-11-03T19:04:00+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-10-30T08:34:42+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"74470171928012631846308803146202059283","date":"2025-10-28T09:12:43+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"159638331570617695158159474755234569539","date":"2025-10-24T07:45:03+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"249636478759369392606943382952406608776","date":"2025-10-24T06:38:35+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2025-10-24T00:49:10+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2025-10-22T18:19:42+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2025-10-20T05:37:26+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2025-10-20T05:37:07+00:00","index":"","fulltext":""},{"type":"submitted","content":"BMC Genomics","date":"2025-10-16T13:23:48+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-genomics","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"gics","sideBox":"Learn more about [BMC Genomics](http://bmcgenomics.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/gics","title":"BMC Genomics","twitterHandle":"#BMCGenomics","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"1ed78d50-5545-4102-8ca3-c9f296179de2","owner":[],"postedDate":"November 5th, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2026-01-12T16:11:23+00:00","versionOfRecord":{"articleIdentity":"rs-7878288","link":"https://doi.org/10.1186/s12864-025-12508-7","journal":{"identity":"bmc-genomics","isVorOnly":false,"title":"BMC Genomics"},"publishedOn":"2026-01-07 15:57:36","publishedOnDateReadable":"January 7th, 2026"},"versionCreatedAt":"2025-11-05 06:40:24","video":"","vorDoi":"10.1186/s12864-025-12508-7","vorDoiUrl":"https://doi.org/10.1186/s12864-025-12508-7","workflowStages":[]},"version":"v1","identity":"rs-7878288","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-7878288","identity":"rs-7878288","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}