CRISPR-Cas9 and Cas12a target site richness reflects genomic diversity in natural populations of Anopheles gambiae and Aedes aegypti mosquitoes.

Due to limitations in conventional disease vector control strategies including the rise of insecticide resistance in natural populations of mosquitoes, genetic control strategies using CRISPR gene drive systems have been under serious consideration. The identification of CRISPR target sites in mosquito populations is a key aspect for developing efficient genetic vector control strategies. While genome-wide Cas9 target sites have been explored in mosquitoes, a precise evaluation of target sites focused on coding sequence (CDS) is lacking. Additionally, target site polymorphisms have not been characterized for other nucleases such as Cas12a, which require a different DNA recognition site (PAM) and would expand the accessibility of mosquito genomes for genetic engineering. We undertook a comprehensive analysis of potential target sites for both Cas9 and Cas12a nucleases within the genomes of natural populations of Anopheles gambiae and Aedes aegypti from multiple continents. We demonstrate that using two nucleases increases the number of targets per gene. Also, we identified differences in nucleotide diversity between North American and African Aedes populations, impacting the abundance of good target sites with a minimal degree of polymorphisms that can affect the binding of gRNA. Lastly, we screened for gRNAs targeting sex-determination genes that could be widely applicable for developing field genetic control strategies. Overall, this work highlights the utility of employing both Cas9 and Cas12a nucleases and underscores the importance of designing universal genetic strategies adaptable to diverse mosquito populations.

Genomic signatures of local adaptation in recent invasive Aedes aegypti populations in California.

BACKGROUND: Rapid adaptation to new environments can facilitate species invasions and range expansions. Understanding the mechanisms of adaptation used by invasive disease vectors in new regions has key implications for mitigating the prevalence and spread of vector-borne disease, although they remain relatively unexplored. RESULTS: Here, we integrate whole-genome sequencing data from 96 Aedes aegypti mosquitoes collected from various sites in southern and central California with 25 annual topo-climate variables to investigate genome-wide signals of local adaptation among populations. Patterns of population structure, as inferred using principal components and admixture analysis, were consistent with three genetic clusters. Using various landscape genomics approaches, which all remove the confounding effects of shared ancestry on correlations between genetic and environmental variation, we identified 112 genes showing strong signals of local environmental adaptation associated with one or more topo-climate factors. Some of them have known effects in climate adaptation, such as heat-shock proteins, which shows selective sweep and recent positive selection acting on these genomic regions. CONCLUSIONS: Our results provide a genome wide perspective on the distribution of adaptive loci and lay the foundation for future work to understand how environmental adaptation in Ae. aegypti impacts the arboviral disease landscape and how such adaptation could help or hinder efforts at population control.

Next-generation gene drive for population modification of the malaria vector mosquito, Anopheles gambiae.

A Cas9/guide RNA-based gene drive strain, AgNosCd-1, was developed to deliver antiparasite effector molecules to the malaria vector mosquito, Anopheles gambiae The drive system targets the cardinal gene ortholog producing a red-eye phenotype. Drive can achieve 98 to 100% in both sexes and full introduction was observed in small cage trials within 6 to 10 generations following a single release of gene-drive males. No genetic load resulting from the integrated transgenes impaired drive performance in the trials. Potential drive-resistant target-site alleles arise at a frequency <0.1, and five of the most prevalent polymorphisms in the guide RNA target site in collections of colonized and wild-derived African mosquitoes do not prevent cleavage in vitro by the Cas9/guide RNA complex. Only one predicted off-target site is cleavable in vitro, with negligible deletions observed in vivo. AgNosCd-1 meets key performance criteria of a target product profile and can be a valuable component of a field-ready strain for mosquito population modification to control malaria transmission.

Genome-wide divergence among invasive populations of Aedes aegypti in California.

BACKGROUND: In the summer of 2013, Aedes aegypti Linnaeus was first detected in three cities in central California (Clovis, Madera and Menlo Park). It has now been detected in multiple locations in central and southern CA as far south as San Diego and Imperial Counties. A number of published reports suggest that CA populations have been established from multiple independent introductions. RESULTS: Here we report the first population genomics analyses of Ae. aegypti based on individual, field collected whole genome sequences. We analyzed 46 Ae. aegypti genomes to establish genetic relationships among populations from sites in California, Florida and South Africa. Based on 4.65 million high quality biallelic SNPs, we identified 3 major genetic clusters within California; one that includes all sample sites in the southern part of the state (South of Tehachapi mountain range) plus the town of Exeter in central California and two additional clusters in central California. CONCLUSIONS: A lack of concordance between mitochondrial and nuclear genealogies suggests that the three founding populations were polymorphic for two main mitochondrial haplotypes prior to being introduced to California. One of these has been lost in the Clovis populations, possibly by a founder effect. Genome-wide comparisons indicate extensive differentiation between genetic clusters. Our observations support recent introductions of Ae. aegypti into California from multiple, genetically diverged source populations. Our data reveal signs of hybridization among diverged populations within CA. Genetic markers identified in this study will be of great value in pursuing classical population genetic studies which require larger sample sizes.

Introgression between Anopheles gambiae and Anopheles coluzzii in Burkina Faso and its associations with kdr resistance and Plasmodium infection.

BACKGROUND: Insecticide resistance in Anopheles coluzzii mosquitoes has become widespread throughout West Africa including in Burkina Faso. The insecticide resistance allele (kdr or L1014F) is a prime indicator that is highly correlated with phenotypic resistance in West Africa. Studies from Benin, Ghana and Mali have suggested that the source of the L1014F is introgression of the 2L divergence island via interspecific hybridization with Anopheles gambiae. The goal of this study was to characterize local mosquito populations in the Nouna Department, Burkina Faso with respect to: (i) the extent of introgression between An. coluzzii and An. gambiae, (ii) the frequency of the L1014F mutation and (iii) Plasmodium infection rates. METHODS: A total of 95 mosquitoes were collected from ten sites surrounding Nouna town in Kossi Province, Burkina Faso in 2012. The species composition, the extent of introgression in An. coluzzii mosquitoes and their Plasmodium infection rates were identified with a modified version of the "Divergence Island SNP" (DIS) genotyping assay. RESULTS: The mosquito collection contained 70.5% An. coluzzii, 89.3% of which carried a 3 Mb genomic region on the 2L chromosome with L1014F insecticide resistance mutation that was introgressed from An. gambiae. In addition, 22.4% in the introgressed An. coluzzii specimens were infected with Plasmodium falciparum, whereas none of the non-introgressed ("pure") An. coluzzii were infected. CONCLUSION: This paper is the first report providing divergence island SNP genotypes for natural population of Burkina Faso and corresponding Plasmodium infection rates. These observations warrant further study and could have a major impact on future malaria control strategies in Burkina Faso.

The Genetic Basis of Host Preference and Resting Behavior in the Major African Malaria Vector, Anopheles arabiensis.

Malaria transmission is dependent on the propensity of Anopheles mosquitoes to bite humans (anthropophily) instead of other dead end hosts. Recent increases in the usage of Long Lasting Insecticide Treated Nets (LLINs) in Africa have been associated with reductions in highly anthropophilic and endophilic vectors such as Anopheles gambiae s.s., leaving species with a broader host range, such as Anopheles arabiensis, as the most prominent remaining source of transmission in many settings. An. arabiensis appears to be more of a generalist in terms of its host choice and resting behavior, which may be due to phenotypic plasticity and/or segregating allelic variation. To investigate the genetic basis of host choice and resting behavior in An. arabiensis we sequenced the genomes of 23 human-fed and 25 cattle-fed mosquitoes collected both in-doors and out-doors in the Kilombero Valley, Tanzania. We identified a total of 4,820,851 SNPs, which were used to conduct the first genome-wide estimates of "SNP heritability" for host choice and resting behavior in this species. A genetic component was detected for host choice (human vs cow fed; permuted P = 0.002), but there was no evidence of a genetic component for resting behavior (indoors versus outside; permuted P = 0.465). A principal component analysis (PCA) segregated individuals based on genomic variation into three groups which were characterized by differences at the 2Rb and/or 3Ra paracentromeric chromosome inversions. There was a non-random distribution of cattle-fed mosquitoes between the PCA clusters, suggesting that alleles linked to the 2Rb and/or 3Ra inversions may influence host choice. Using a novel inversion genotyping assay, we detected a significant enrichment of the standard arrangement (non-inverted) of 3Ra among cattle-fed mosquitoes (N = 129) versus all non-cattle-fed individuals (N = 234; χ2, p = 0.007). Thus, tracking the frequency of the 3Ra in An. arabiensis populations may be of use to infer selection on host choice behavior within these vector populations; possibly in response to vector control. Controlled host-choice assays are needed to discern whether the observed genetic component has a direct relationship with innate host preference. A better understanding of the genetic basis for host feeding behavior in An. arabiensis may also open avenues for novel vector control strategies based on driving genes for zoophily into wild mosquito populations.

The fate of genes that cross species boundaries after a major hybridization event in a natural mosquito population.

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.

Adaptive introgression in an African malaria mosquito coincident with the increased usage of insecticide-treated bed nets.

Animal species adapt to changes in their environment, including man-made changes such as the introduction of insecticides, through selection for advantageous genes already present in populations or newly arisen through mutation. A possible alternative mechanism is the acquisition of adaptive genes from related species via a process known as adaptive introgression. Differing levels of insecticide resistance between two African malaria vectors, Anopheles coluzzii and Anopheles gambiae, have been attributed to assortative mating between the two species. In a previous study, we reported two bouts of hybridization observed in the town of Selinkenyi, Mali in 2002 and 2006. These hybridization events did not appear to be directly associated with insecticide-resistance genes. We demonstrate that during a brief breakdown in assortative mating in 2006, A. coluzzii inherited the entire A. gambiae-associated 2L divergence island, which includes a suite of insecticide-resistance alleles. In this case, introgression was coincident with the start of a major insecticide-treated bed net distribution campaign in Mali. This suggests that insecticide exposure altered the fitness landscape, favoring the survival of A. coluzzii/A. gambiae hybrids, and provided selection pressure that swept the 2L divergence island through A. coluzzii populations in Mali. We propose that the work described herein presents a unique description of the temporal dynamics of adaptive introgression in an animal species and represents a mechanism for the rapid evolution of insecticide resistance in this important vector of human malaria in Africa.

Spatiotemporal dynamics of gene flow and hybrid fitness between the M and S forms of the malaria mosquito, Anopheles gambiae.

The M and S forms of Anopheles gambiae have been the focus of intense study by malaria researchers and evolutionary biologists interested in ecological speciation. Divergence occurs at three discrete islands in genomes that are otherwise nearly identical. An "islands of speciation" model proposes that diverged regions contain genes that are maintained by selection in the face of gene flow. An alternative "incidental island" model maintains that gene flow between M and S is effectively zero and that divergence islands are unrelated to speciation. A "divergence island SNP" assay was used to explore the spatial and temporal distributions of hybrid genotypes. Results revealed that hybrid individuals occur at frequencies ranging between 5% and 97% in every population examined. A temporal analysis revealed that assortative mating is unstable and periodically breaks down, resulting in extensive hybridization. Results suggest that hybrids suffer a fitness disadvantage, but at least some hybrid genotypes are viable. Stable introgression of the 2L speciation island occurred at one site following a hybridization event.

Anopheles darlingi polytene chromosomes: revised maps including newly described inversions and evidence for population structure in Manaus.

Salivary gland polytene chromosomes of 4th instar Anopheles darlingi Root were examined from multiple locations in the Brazilian Amazon. Minor modifications were made to existing polytene photomaps. These included changes to the breakpoint positions of several previously described paracentric inversions and descriptions of four new paracentric inversions, two on the right arm of chromosome 3 and two on the left arm of chromosome 3 that were found in multiple locations. A total of 18 inversions on the X (n = 1) chromosome, chromosome 2 (n = 7) and 3 (n = 11) were scored for 83 individuals from Manaus, Macapá and Porto Velho municipalities. The frequency of 2Ra inversion karyotypes in Manaus shows significant deficiency of heterozygotes (p < 0.0009). No significant linkage disequilibrium was found between inversions on chromosome 2 and 3. We hypothesize that at least two sympatric subpopulations exist within the An. darlingi population at Manaus based on inversion frequencies.

Complex genome evolution in Anopheles coluzzii associated with increased insecticide usage in Mali.

In certain cases, a species may have access to important genetic variation present in a related species via adaptive introgression. These novel alleles may interact with their new genetic background, resulting in unexpected phenotypes. In this study, we describe a selective sweep on standing variation on the X chromosome in the mosquito Anopheles coluzzii, a principal malaria vector in West Africa. This event may have been influenced by the recent adaptive introgression of the insecticide resistance gene known as kdr from the sister species Anopheles gambiae. Individuals carrying both kdr and a nearly fixed X-linked haplotype, encompassing at least four genes including the P450 gene CYP9K1 and the cuticular protein CPR125, have rapidly increased in relative frequency. In parallel, a reproductively isolated insecticide-susceptible A. gambiae population (Bamako form) has been driven to local extinction, likely due to strong selection from increased insecticide-treated bed net usage.

Colonization of malaria vectors under semi-field conditions as a strategy for maintaining genetic and phenotypic similarity with wild populations.

BACKGROUND: Malaria still accounts for an estimated 207 million cases and 627,000 deaths worldwide each year. One proposed approach to complement existing malaria control methods is the release of genetically-modified (GM) and/or sterile male mosquitoes. As opposed to laboratory colonization, this requires realistic semi field systems to produce males that can compete for females in nature. This study investigated whether the establishment of a colony of the vector Anopheles arabiensis under more natural semi-field conditions can maintain higher levels of genetic diversity than achieved by laboratory colonization using traditional methods. METHODS: Wild females of the African malaria vector An. arabiensis were collected from a village in southern Tanzania and used to establish new colonies under different conditions at the Ifakara Health Institute. Levels of genetic diversity and inbreeding were monitored in colonies of An. arabiensis that were simultaneously established in small cage colonies in the SFS and in a large semi-field (SFS) cage and compared with that observed in the original founder population. Phenotypic traits that determine their fitness (body size and energetic reserves) were measured at 10(th) generation and compared to founder wild population. RESULTS: In contrast to small cage colonies, the SFS population of An. arabiensis exhibited a higher degree of similarity to the founding field population through time in several ways: (i) the SFS colony maintained a significantly higher level of genetic variation than small cage colonies, (ii) the SFS colony had a lower degree of inbreeding than small cage colonies, and (iii) the mean and range of mosquito body size in the SFS colony was closer to that of the founding wild population than that of small cage colonies. Small cage colonies had significantly lower lipids and higher glycogen abundances than SFS and wild population. CONCLUSIONS: Colonization of An. arabiensis under semi-field conditions was associated with the retention of a higher degree of genetic diversity, reduced inbreeding and greater phenotypic similarity to the founding wild population than observed in small cage colonies. Thus, mosquitoes from such semi-field populations are expected to provide more realistic representation of mosquito ecology and physiology than those from small cage colonies.

Defining Genetic, Taxonomic, and Geographic Boundaries Among Species of the Psorophora confinnis (Diptera: Culicidae) Complex in North and South America.

The Psorophora confinnis complex is currently composed of three species--Psorophora confinnis sensu stricto (Lynch Arribalzaga) in South America, Psorophora columbiae (Dyar and Knab) in North America, and Psorophora jamaicensis (Theobald) in the Caribbean. Members of the complex are of considerable importance as vectors of arboviruses, for example, Venezuelan equine encephalitis virus, and are significant biting pests throughout their range. The biological and geographic boundaries of Ps. confinnis and Ps. columbiae are unclear. In fact, the name Ps. columbiae is presently designated as "provisional." In this article, we aim to clarify the taxonomy and geographic distributions of species within the Ps. confinnis complex. A population genetics approach was employed using gene and genotypic frequency data at 26 isozyme loci. The results suggest that the Ps. confinnis complex in North and South America is composed of four species. Ps. confinnis s.s. and Ps. columbiae are distinct species in South and North America, respectively. Populations in Colombia, South America, formally designated as Ps. funiculus (Dyar) and populations in the southwestern United States and western Mexico, formally designated Ps. toltecum (Dyar and Knab), are distinct species. Psorophora toltecum and Psorophora funiculus species names should be resurrected from synonymy. In addition we identified a Ps. columbiae and Ps. toltecum hybrid zone in central Texas in a region described as being one of 13 North American suture zones, being geographical areas in which closely related species occur in sympatry and frequently hybridize.

Resolving confusion in the native distribution of Aedes (Hulecoeteomyia) koreicus (Diptera: Culicidae).

Recent studies on invasive Aedes koreicus (Edwards 1917) have been conducted to elucidate the source population introduced to Europe. However, current information about the native range of Ae. koreicus is not consistent. The purpose of this study is to resolve confusion in the native distribution of Ae. koreicus by reviewing available literature from the first description of the species in its native range in 1917 to the first invasion in Europe in 2008. Aedes koreicus have been recorded in China, Japan, Korea, and eastern Russia. The 2 existing records of Ae. koreicus from Hokkaido, Japan, however, is likely due to the misidentification of 2 different morphologically similar species, Ae. koreicoides (Sasa, Kano & Hayashi 1950) and Aedes japonicus (Theobald 1901). Upon re-examination of published records, we conclude that the native distribution of Ae. koreicus is confined to continental eastern Asian regions, specifically China, Korea, and eastern Russia.

A preliminary investigation of the relationship between water quality and Anopheles gambiae larval habitats in Western Cameroon.

BACKGROUND: Water quality and anopheline habitat have received increasing attention due to the possibility that challenges during larval life may translate into adult susceptibility to malaria parasite infection and/or insecticide resistance. METHODS: A preliminary study of Anopheles gambiae s.s. larval habitats in the north-west and south-west regions of Cameroon was conducted in order to detect associations between An. gambiae s.s. molecular form and 2La inversion distributions with basic water quality parameters. Water quality was measured by temperature, pH, conductivity, total dissolved solids (TDS) at seven sites in Cameroon and one site in Selinkenyi, Mali. RESULTS: Principal components and correlation analyses indicated a complex relationship between 2La polymorphism, temperature, conductivity and TDS. Cooler water sites at more inland locations yielded more S form larvae with higher 2La inversion polymorphism while warmer water sites yielded more M form larvae with rare observations of the 2La inversion. DISCUSSION: More detailed studies that take into account the population genetics but also multiple life stages, environmental data relative to these life stages and interactions with both humans and the malaria parasite may help us to understand more about how and why this successful mosquito is able to adapt and diverge, and how it can be successfully managed.

The potential invasion into North America and Europe by non-native mosquito, Aedes koreicus (Diptera: Culicidae).

Aedes koreicus (Edward, 1917) (Diptera: Culicidae), a mosquito species native to East Asia, has spread to parts of Europe and Central Asia since 2008. The species shares ecological characteristics with Aedes japonicus (Theobald, 1901) (Diptera: Culicidae), which has already successfully invaded and established in North America and Europe. Given these similarities, it is plausible that Ae. koreicus may also invade North America in the future. However, the invasion of Ae. koreicus may be masked or have delayed detection due to their similar morphologies with Ae. japonicus. This study highlights the potential risks of invasion of Ae. koreicus into North America, especially in the northeastern United States, and for further expansion in Europe. We used the maximum entropy model to identify areas with a high likelihood of presence in North America and Europe using comprehensive occurrence records from East Asia, Central Asia, and Europe. We have identified 15 additional countries in Europe and 7 states in the United States that will likely have suitable environments for Ae. koreicus. Additionally, we reviewed the morphological characteristics of Ae. koreicus and Ae. japonicus and provided morphological keys to distinguish the 2 species. Morphological results contradicting previous studies suggested that finding the origin by morphological comparison between Ae. koreicus populations may need re-evaluation. The information presented here will be useful for researchers and public health professionals in high-risk areas to be informed about morphological characteristics to distinguish Ae. koreicus from similar-looking Ae. japonicus. These tools will allow more careful monitoring of the potential introduction of this highly invasive species.

The first genome sequence of Anopheles squamous from Macha, Zambia.

Despite efforts to minimize the impacts of malaria and reduce the number of primary vectors, malaria has yet to be eliminated in Zambia. Understudied vector species may perpetuate malaria transmission in pre-elimination settings. Anopheles squamosus is one of the most abundantly caught mosquito species in southern Zambia and has previously been found with Plasmodium falciparum sporozoites, a causal agent of human malaria. This species may be a critical vector of malaria transmission, however, there is a lack of genetic information available for An. squamosus. We report the first genome data and the first complete mitogenome (Mt) sequence of An. squamosus. The sequence was extracted from one individual mosquito from the Chidakwa area in Macha, Zambia. The raw reads were obtained using Illumina Novaseq 6000 and assembled through NOVOplasty alignment with related species. The length of the An. squamosus Mt was 15,351 bp, with 77.9 % AT content. The closest match to the whole mitochondrial genome in the phylogenetic tree is the African malaria mosquito, Anopheles gambiae. Its genome data is available through National Center for Biotechnology Information (NCBI) Sequencing Reads Archive (SRA) with accession number SRR22114392. The mitochondrial genome was deposited in NCBI GenBank with the accession number OP776919. The ITS2 containing contig sequence was deposited in GenBank with the accession number OQ241725. Mitogenome annotation and a phylogenetic tree with related Anopheles mosquito species are provided.

Field Evaluation of In2Care Mosquito Traps to Control Aedes aegypti and Aedes albopictus (Diptera: Culicidae) in Hawai'i Island.

Aedes aegypti Linnaeus and Aedes albopictus Skuse are vectors of dengue virus and responsible for multiple autochthonous dengue outbreaks in Big Island, Hawai'i. Control of Ae. aegypti and Ae. albopictus has been achieved in In2Care trap trials, which motivated us to investigate this potential control approach in the Big Island. Our In2Care trial was performed in the coastal settlement of Miloli'i in the southwest of Big Island where both Ae. aegypti and Ae. albopictus are found. This trial starting in the second week of July and ending in the last week of October 2019 fell within the traditional wet season in Miloli'i. No significant reduction in egg or adult counts in our treatment areas following 12 wk of two In2Care trap placements per participating household were observed. In fact, an increase in numbers of adults during the trial reached levels that required the local mosquito abatement program to stop the In2Care trap trial and institute a thorough source reduction and treatment campaign. The source reduction campaign revealed a large variety and quantity of water sources competed with the oviposition cups we had placed, which likely lowered the chances of our oviposition cups being visited by pyriproxyfen-contaminated Aedes adults exiting the In2Care traps.

Complete mitogenome sequence of Aedes (Hulecoeteomyia) japonicus japonicus from Hawai'i Island.

We report the first complete mitogenome (Mt) sequence of Aedes japonicus japonicus (Diptera: Culicidae). The sequence was extracted from one adult from the Big Island of Hawai'i Island. The length of the Ae. japonicus japonicus Mt was 16,528bp with 78.1% AT content. Its sequence is most similar to the Mt sequence of Aedes koreicus with 90.81% sequence identity. This is the first full Mt sequence available for this species and provides important genetic resource for studying population genetics and dynamics of this important invasive mosquito species.

Anopheles gambiae on remote islands in the Indian Ocean: origins and prospects for malaria elimination by genetic modification of extant populations.

The mosquito Anopheles gambiae s.s. is a primary malaria vector throughout sub-Saharan Africa including the islands of the Comoros archipelago (Anjouan, Grande Comore, Mayotte and Mohéli). These islands are located at the northern end of the Mozambique Channel in eastern Africa. Previous studies have shown a relatively high degree of genetic isolation between the Comoros islands and mainland populations of A. gambiae, but the origin of the island populations remains unclear. Here, we analyzed phylogenetic relationships among island and mainland populations using complete mitochondrial genome sequences of individual A. gambiae specimens. This work augments earlier studies based on analysis of the nuclear genome. We investigated the source population of A. gambiae for each island, estimated the number of introductions, when they occurred and explored evidence for contemporary gene flow between island and mainland populations. These studies are relevant to understanding historical patterns in the dispersal of this important malaria vector and provide information critical to assessing their potential for the exploration of genetic-based vector control methods to eliminate this disease. Phylogenetic analysis and haplotype networks were constructed from mitogenome sequences of 258 A. gambiae from the four islands. In addition, 112 individuals from seven countries across sub-Saharan Africa and Madagascar were included to identify potential source populations. Our results suggest that introduction events of A. gambiae into the Comoros archipelago were rare and recent events and support earlier claims that gene flow between the mainland and these islands is limited. This study is concordant with earlier work suggesting the suitability of these oceanic islands as appropriate sites for conducting field trial releases of genetically engineered mosquitoes (GEMs).