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.

Population modification strategies for malaria vector control are uniquely resilient to observed levels of gene drive resistance alleles.

Cas9/guide RNA (gRNA)-based gene drive systems are expected to play a transformative role in malaria elimination efforts., whether through population modification, in which the drive system contains parasite-refractory genes, or population suppression, in which the drive system induces a severe fitness load resulting in population decline or extinction. DNA sequence polymorphisms representing alternate alleles at gRNA target sites may confer a drive-resistant phenotype in individuals carrying them. Modeling predicts that, for observed levels of SGV at potential target sites and observed rates of de novo DRA formation, population modification strategies are uniquely resilient to DRAs. We conclude that gene drives can succeed when fitness costs incurred by drive-carrying mosquitoes are low enough to prevent strong positive selection for DRAs produced de novo or as part of the SGV and that population modification strategies are less prone to failure due to drive resistance.

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.

Practical Application of a Relationship-Based Model to Engagement for Gene-Drive Vector Control Programs.

Engagement is an important component in the advancement of gene-drive vector control research programs as developers look to transition the technology from the laboratory to the field. As research advances and engagement surrounding this novel technology is put into practice, knowledge can be gained from practical experiences and applications in the field. A relationship-based model (RBM) provides a framework for end-user development of engagement programs and strategies. The model places end users at the center of the engagement decision-making processes rather than as recipients of predetermined strategies, methods, and definitions. Successful RBM application for healthcare delivery has previously been demonstrated, and the University of California Malaria Initiative (UCMI) has applied this model to its gene-drive program in the Democratic Republic of São Tomé and Príncipe. The model emphasizes the importance of local leadership in the planning, development, and implementation of all phases of project engagement. The primary aim of this paper is to translate the model from paper to practice and provide a transparent description, using practical examples, of the UCMI program implementation of RBM at its field site. End-user development of the UCMI engagement program provides a unique approach to the development of ethical, transparent, and effective engagement strategies for malaria control programs. This paper may also serve as a reference and example for projects looking to establish an engagement program model that integrates end-user groups in the decision-making processes surrounding engagement.

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.

Conceptual risk assessment of mosquito population modification gene-drive systems to control malaria transmission: preliminary hazards list workshops.

The field-testing and eventual adoption of genetically-engineered mosquitoes (GEMs) to control vector-borne pathogen transmission will require them meeting safety criteria specified by regulatory authorities in regions where the technology is being considered for use and other locales that might be impacted. Preliminary risk considerations by researchers and developers may be useful for planning the baseline data collection and field research used to address the anticipated safety concerns. Part of this process is to identify potential hazards (defined as the inherent ability of an entity to cause harm) and their harms, and then chart the pathways to harm and evaluate their probability as part of a risk assessment. The University of California Malaria Initiative (UCMI) participated in a series of workshops held to identify potential hazards specific to mosquito population modification strains carrying gene-drive systems coupled to anti-parasite effector genes and their use in a hypothetical island field trial. The hazards identified were placed within the broader context of previous efforts discussed in the scientific literature. Five risk areas were considered i) pathogens, infections and diseases, and the impacts of GEMs on human and animal health, ii) invasiveness and persistence of GEMs, and interactions of GEMs with target organisms, iii) interactions of GEMs with non-target organisms including horizontal gene transfer, iv) impacts of techniques used for the management of GEMs and v) evolutionary and stability considerations. A preliminary hazards list (PHL) was developed and is made available here. This PHL is useful for internal project risk evaluation and is available to regulators at prospective field sites. UCMI project scientists affirm that the subsequent processes associated with the comprehensive risk assessment for the application of this technology should be driven by the stakeholders at the proposed field site and areas that could be affected by this intervention strategy.

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).

Differential Plasmodium falciparum infection of Anopheles gambiae s.s. molecular and chromosomal forms in Mali.

BACKGROUND: Anopheles gambiae sensu stricto (s.s.) is a primary vector of Plasmodium falciparum in sub-Saharan Africa. Although some physiological differences among molecular and chromosomal forms of this species have been demonstrated, the relative susceptibility to malaria parasite infection among them has not been unequivocally shown. The objective of this study was to investigate P. falciparum circumsporozoite protein infection (CSP) positivity among An. gambiae s.s. chromosomal and molecular forms. METHODS: Wild An. gambiae from two sites Kela (n=464) and Sidarebougou (n=266) in Mali were screened for the presence of P. falciparum CSP using an enzyme-linked immunosorbent assay (ELISA). Samples were then identified to molecular form using multiple PCR diagnostics (n=713) and chromosomal form using chromosomal karyotyping (n=419). RESULTS: Of 730 An. gambiae sensu lato (s.l.) mosquitoes, 89 (12.2%) were CSP ELISA positive. The percentage of positive mosquitoes varied by site: 52 (11.2%) in Kela and 37 (13.9%) in Sidarebougou. Eighty-seven of the positive mosquitoes were identified to molecular form and they consisted of nine Anopheles arabiensis (21.4%), 46 S (10.9%), 31 M (12.8%), and one MS hybrid (14.3%). Sixty of the positive mosquitoes were identified to chromosomal form and they consisted of five An. arabiensis (20.0%), 21 Savanna (15.1%), 21 Mopti (30.4%), 11 Bamako (9.2%), and two hybrids (20.0%). DISCUSSION: In this collection, the prevalence of P. falciparum infection in the M form was equivalent to infection in the S form (no molecular form differential infection). There was a significant differential infection by chromosomal form such that, P. falciparum infection was more prevalent in the Mopti chromosomal forms than in the Bamako or Savanna forms; the Mopti form was also the most underrepresented in the collection. Continued research on the differential P. falciparum infection of An. gambiae s.s. chromosomal and molecular forms may suggest that Plasmodium - An. gambiae interactions play a role in malaria transmission.

Identification and characterization of single nucleotide polymorphisms (SNPs) in Culex theileri (Diptera: Culicidae).

Culex theileri Theobald (Diptera: Culicidae) is one of the most common mosquito species in northeastern Turkey and serves as a vector for various zoonotic diseases including West Nile virus. Although there have been some studies on the ecology of Cx. theileri, very little genetic data has been made available. We successfully sequenced 11 gene fragments from Cx. theileri specimens collected from the northeastern part of Turkey. On average, we found a Single nucleotide polymorphism every 45 bp. Transitions outnumbered transversions, at a ratio of 2:1. This is the first report of genetic polymorphisms in Cx. theileri and Single nucleotide polymorphism discovered from this study can be used to investigate population structure and gene-environmental interactions.

Single-nucleotide polymorphisms for high-throughput genotyping of Anopheles arabiensis in East and southern Africa.

Anopheles arabiensis Patton is one of the principal vectors of malaria in sub-Saharan Africa, occupying a wide variety of ecological zones. This species is increasingly responsible for malaria transmission in Africa and is becoming the dominant vector species in some localities. Despite its growing importance, little is known about genetic polymorphisms in this species. Multiple sequences of various gene fragments from An. arabiensis isolates from Cameroon were obtained from GenBank. In total, 20 gene fragments containing single-nucleotide polymorphisms (SNPs) at moderate density were selected for direct sequencing from field collected specimens from Tanzania and Zambia. We obtained 301 SNPs in total from the 20 gene fragments, 60 of which were suitable for Illumina GoldenGate SNP genotyping. A greater number of SNPs (n = 185) was suitable for analysis using Sequenom iPLEX, an alternative high-throughput genotyping technology using mass spectrometry. An SNP was present every 59 (+/- 44.5) bases on average. Overall, An. arabiensis from Tanzania and Zambia are genetically closer (mean F(ST) = 0.075) than either is to populations in Cameroon (F(ST, TZ-CM) = 0.250, F(ST,ZA-CM) = 0.372). A fixed polymorphism between East/southern and Central Africa was identified on AGAP000574, a gene on the X chromosome. We have identified SNPs in natural populations of An. arabiensis. SNP densities in An. arabiensis were higher than Anopheles gambiae s.s., suggesting a greater challenge in the development of high-throughput SNP analysis for this species. The SNP markers provided in this study are suitable for a high-throughput genotyping analysis and can be used for population genetic studies and association mapping efforts.