Genome-wide association studies unveil major genetic loci driving insecticide resistance in Anopheles funestus in four eco-geographical settings across Cameroon.

Abstract

Background: Insecticide resistance is jeopardising malaria control efforts in Africa. Deciphering the evolutionary dynamics of mosquito populations country-wide is essential for designing effective and sustainable national and subnational tailored strategies to accelerate malaria elimination efforts. Here, we employed genome-wide association studies through pooled template sequencing to compare four eco-geographically different populations of the major vector, Anopheles funestus, across a South North transect in Cameroon, aiming to identify genomic signatures of adaptive responses to insecticides.

Results: Our analysis revealed limited population structure within Northern and Central regions (F_ST<0.02), suggesting extensive gene flow, while populations from the Littoral/Coastal region exhibited more distinct genetic patterns (F_ST>0.049). Greater genetic differentiation was observed at known resistance-associated loci, resistance-to-pyrethroids 1 (rp1) (2R chromosome) and CYP9 (X chromosome), with varying signatures of positive selection across populations. Allelic variation between variants underscores the pervasive impact of selection pressures, with rp1 variants more prevalent in Central and Northern populations (F_ST>0.3), and the CYP9 associated variants more pronounced in the Littoral/Coastal region (F_ST =0.29). Evidence of selective sweeps was supported by negative Tajima's D and reduced genetic diversity in all populations, particularly in Central (Elende) and Northern (Tibati) regions. Genomic variant analysis identified novel missense mutations and signatures of complex genomic alterations such as duplications, deletions, transposable element (TE) insertions, and chromosomal inversions, all associated with selective sweeps. A 4.3 kb TE insertion was fixed in all populations with Njombe Littoral/Coastal population, showing higher frequency of CYP9K1 (G454A), a known resistance allele and TE upstream compared to elsewhere.

Conclusion: Our study uncovered regional variations in insecticide resistance candidate variants, emphasizing the need for a streamlined DNA-based diagnostic assay for genomic surveillance across Africa. These findings will contribute to the development of tailored resistance management strategies crucial for addressing the dynamic challenges of malaria control in Cameroon.