Genomic analyses revealed low genetic variation in the intron-exon boundary of the doublesex gene within the natural populations of An. gambiae s.l. in Burkina Faso

BackgroundThe recent success of a population control gene drive targeting the doublesex gene in Anopheles gambiae paved the way for developing self-sustaining and self-limiting genetic control strategies targeting the sex determination pathway to reduce and/or distort the reproductive capacity of insect vectors. However, targeting these genes for genetic control requires a better understanding of their genetic variation in natural populations to ensure effective gene drive spread. Using whole genome sequencing (WGS) data from the Ag1000G project (Ag3.0, 3.4 and 3.8), and Illumina pooled amplicon sequencing, we investigated the genetic polymorphism of the intron-4-exon-5 boundary of the doublesex gene in the natural populations of An. gambiae sensu lato (s.l.). ResultsThe analyses showed a very low variant density at the gRNA target sequence of the Ag(QFS)1 gene drive (previously called dsxFCRISPRh) within the populations of West and East Africa. However, populations from the forest area in Central Africa exhibited four SNP at frequencies ranging from 0.011 to 0.26. The SNP (2R:48714641[C > T]) at high frequencies, i.e. 0.26 is identified within the An. coluzzii population from Angola. The analyses also identified 90 low frequency (1 - 5%) SNPs in the genomic region around the gRNA target sequence (intron-4-exon-5 boundary). Three of these SNPs (2R:48714472 A > T; 2R:48714486 C > A; 2R:48714516 C > T) were observed at frequencies higher than 5% in the UTR region of the doublesex gene. The results also showed a very low variant density and constant nucleotide diversity over a five-year survey in natural An. gambiae s.l. populations of Burkina Faso. ConclusionThese findings will guide the implementation of doublesex-targeted gene drives to support the current control tools in malaria elimination efforts. Our methods can be applied to efficiently monitor the evolution of any sequence of interest in a natural population via pooled amplicon sequencing, surpassing the need for WGS.