Animal species are able to acquire new genetic material via hybridization and subsequent introgression. However, little is known about how foreign genomic material is incorporated into a population over time and what genes are susceptible to introgression. Here, we follow the closely related mosquito sister species Anopheles coluzzii and Anopheles gambiae in a sympatric natural population in Mali at multiple time points spanning a period of 25 years. During this period, we observed the temporary breakdown of mating barriers, which allowed us to explore the fate of alleles that crossed the species boundary in a natural population. Whole genome sequencing of 74 individuals revealed introgression within only 34 genes (0.26% of total genes) from A. gambiae to A. coluzzii, the majority contained within a 4 Mb region on the 2L chromosome which includes the insecticide resistance gene (AGAP004707). We designed a genotyping assay to follow 25 of the 34 introgressed alleles over time and found that all A. gambiae alleles, except four, reached a frequency of 50% in the A. coluzzii population within 4 years (~50 generations) and increased to ~80% within 6 years (~75 generations). However, the frequency of all introgressed alleles, except three, decreased to ~60% in 2016. This suggests an ongoing process of purifying selection in the population against DNA of foreign ancestry, except for alleles that are under positive selection, resulting in a complex genomic landscape. This study shows that stable introgression is limited to only specific genes even within closely related species.

Knockdown resistance ( kdr) and CYP9K1 genotypes were detected by a MOLDI-TOF based SNP genotyping assay (Sequenom iPLEX) in samples of Anopheles gambiae collected at 13 sites throughout the Union of the Comoros and Dar es Salaam, Tanzania during February and March 2011. All A. gambiae specimens collected in the Comoros were homozygous for the susceptible kdr alleles (+/+) while 96% of A. gambiae from Dar es Salaam were homozygous for the East African kdr resistant genotype (E/E). In contrast, all specimens from Dar es Salaam and the Comoros were homozygous for the cyp3 allele (c3/c3) at the CYP9K1 locus; the locus has been implicated in metabolic resistance against pyrethroid insecticides in West Africa. All specimens had typical A. gambiae genotypes for SNPs within the divergence Islands on all three chromosomes. Although further spatial and temporal studies are needed, the distribution of kdr genotypes between the Comoros and Tanzania further supports isolation of the Comoros populations from A. gambiae populations on mainland Africa .

No abstract