Genomic evidence of spatially structured gene flow and divergent insecticide resistance backgrounds of the malaria vector Anopheles funestus in Tanzania.

Population genetic analysis of mosquitoes is important for understanding the distribution of insecticide resistance alleles, devising sustainable control approaches, and understanding how vector populations are structured in space. Anopheles funestus is the dominant malaria vector in most parts of East and Southern Africa. To better understand its population genetic structure in Tanzania, we sequenced the genomes of 334 individual An. funestus mosquitoes from 11 regions across the country. Signs of reduced migration between western and eastern cohorts across the semi-arid central region containing the Rift Valley suggest a partial barrier to gene flow between these populations. This was evidenced by population structure between the eastern and western cohorts, as well as asynchronous selective sweeps and copy number variant profiles at the Cyp9k1 gene and Cyp6p gene cluster. Eastern cohorts, despite having less diversity and greater inbreeding, also share genetic histories characterized by low genome-wide Fst values with those in the west. This suggests that the barrier to gene flow is porous and likely represents continuous spatial structure rather than a complete barrier to migration. The observed population disconnectedness should be considered for insecticide deployment, resistance management, and the rollout of novel genetic-based vector control approaches. These findings provide the most detailed study of Tanzanian An. funestus population structure and resistance genetics to date. Future research should examine the epidemiological relevance of this partial discontinuity in gene flow and whether these populations have different malaria transmission abilities.