Inversion 2La is associated with enhanced desiccation resistance in Anopheles gambiae

Background

Anopheles gambiae, the principal vector of malignant malaria in Africa, occupies a wide range of habitats. Environmental flexibility may be conferred by a number of chromosomal inversions non-randomly associated with aridity, including 2La. The purpose of this study was to determine the physiological mechanisms associated with the 2La inversion that may result in the preferential survival of its carriers in hygrically-stressful environments.

Methods

Two homokaryotypic populations of A. gambiae (inverted 2La and standard 2L+^a) were created from a parental laboratory colony polymorphic for 2La and standard for all other known inversions. Desiccation resistance, water, energy and dry mass of adult females of both populations were compared at several ages and following acclimation to a more arid environment.

Results

Females carrying 2La were significantly more resistant to desiccation than 2L+^a females at emergence and four days post-emergence, for different reasons. Teneral 2La females had lower rates of water loss than their 2L+^a counterparts, while at four days, 2La females had higher initial water content. No differences in desiccation resistance were found at eight days, with or without acclimation. However, acclimation resulted in both populations significantly reducing their rates of water loss and increasing their desiccation resistance. Acclimation had contrasting effects on the body characteristics of the two populations: 2La females boosted their glycogen stores and decreased lipids, whereas 2La females did the contrary.

Conclusion

Variation in rates of water loss and response to acclimation are associated with alternative arrangements of the 2La inversion. Understanding the mechanisms underlying these traits will help explain how inversion polymorphisms permit exploitation of a heterogeneous environment by this disease vector.

2La chromosomal inversion enhances thermal tolerance of Anopheles gambiae larvae

BACKGROUND

The mosquito Anopheles gambiae is broadly distributed throughout sub-Saharan Africa and this contributes to making it the most efficient vector of malaria on the continent. The pervasiveness of this species is hypothesized to originate in local adaptations facilitated by inversion polymorphisms. One inversion, named 2La, is strongly associated with aridity clines in West and Central Africa: while 2La is fixed in arid savannas, the 2L+a arrangement is predominantly found in the rainforest. Ability to survive high temperature exposure is an essential component of aridity tolerance, particularly in immature stages that are restricted to shallow puddles. Toward deciphering the role of the 2La inversion in local adaptation, the present investigation focused on variation in larval and pupal thermo-tolerance in two populations dissimilar solely in 2La arrangement.

METHODS

A laboratory colony of A. gambiae that is polymorphic for 2La but standard for all other known inversions was used to create 2 homokaryotypic populations (2L+a and 2La). The survival of 4th instar larvae and pupae from both populations was then tested following exposure to thermal stress with and without prior heat hardening.

RESULTS

Larvae responded identically to a 40 degrees C heat stress, with about 50% of larvae dying after 1.5-2 h and few larvae surviving a 3 h stress. When heat hardened prior to the thermal stress, thermo-tolerance of both larval populations increased, with 2La 24 h survival significantly exceeding that of 2L+a. Pupae were generally more thermo-tolerant than larvae, although 2La pupae were less so than 2L+a. Heat hardening had no positive effect on pupal thermo-tolerance.

CONCLUSION

The increased thermo-tolerance observed in 2La larvae following heat hardening suggests higher responsiveness (i.e., thermal sensitivity) of the inverted karyotype. By responding more drastically to the heat shock, 2La larvae are better equipped to resist the potentially lethal temperatures that occur in arid habitats. The lower survival of 2La pupae compared with 2L+a may reflect the cost of this sensitivity, whereby the thermal resistance mechanisms prevent successful completion of metamorphosis. The costs and benefits of thermal resistance are discussed in light of the climates characterizing either end of the 2La frequency cline.

Living at the edge: biogeographic patterns of habitat segregation conform to speciation by niche expansion in Anopheles gambiae

BACKGROUND

Ongoing lineage splitting within the African malaria mosquito Anopheles gambiae is compatible with ecological speciation, the evolution of reproductive isolation by divergent natural selection acting on two populations exploiting alternative resources. Divergence between two molecular forms (M and S) identified by fixed differences in rDNA, and characterized by marked, although incomplete, reproductive isolation is occurring in West and Central Africa. To elucidate the role that ecology and geography play in speciation, we carried out a countrywide analysis of An. gambiae M and S habitat requirements, and that of their chromosomal variants, across Burkina Faso.

RESULTS

Maps of relative abundance by geostatistical interpolators produced a distinct pattern of distribution: the M-form dominated in the northernmost arid zones, the S-form in the more humid southern regions. Maps of habitat suitability, quantified by Ecological Niche Factor Analysis based on 15 eco-geographical variables revealed less contrast among forms. M was peculiar as it occurred proportionally more in habitat of marginal quality. Measures of ecological niche breadth and overlap confirmed the mismatch between the fundamental and realized patterns of habitat occupation: forms segregated more than expected from the extent of divergence of their environmental envelope--a signature of niche expansion. Classification of chromosomal arm 2R karyotypes by multilocus genetic clustering identified two clusters loosely corresponding to molecular forms, with 'mismatches' representing admixed individuals due to shared ancestral polymorphism and/or residual hybridization. In multivariate ordination space, these karyotypes plotted in habitat of more marginal quality compared to non-admixed, 'typical', karyotypes. The distribution of 'typical' karyotypes along the main eco-climatic gradient followed a consistent pattern within and between forms, indicating an adaptive role of inversions at this geographical scale.

CONCLUSION

Ecological segregation between M and S is consistent with niche expansion into marginal habitats by chromosomal inversion variants during early lineage divergence; presumably, this process is promoted by inter-karyotype competition in the higher-quality core habitat. We propose that the appearance of favourable allelic combinations in other regions of suppressed recombination (e.g. pericentromeric portions defining speciation islands in An. gambiae) fosters development of reproductive isolation to protect linkage between separate chromosomal regions.

Ecological niche partitioning between Anopheles gambiae molecular forms in Cameroon: the ecological side of speciation

BACKGROUND

Speciation among members of the Anopheles gambiae complex is thought to be promoted by disruptive selection and ecological divergence acting on sets of adaptation genes protected from recombination by polymorphic paracentric chromosomal inversions. However, shared chromosomal polymorphisms between the M and S molecular forms of An. gambiae and insufficient information about their relationship with ecological divergence challenge this view. We used Geographic Information Systems, Ecological Niche Factor Analysis, and Bayesian multilocus genetic clustering to explore the nature and extent of ecological and chromosomal differentiation of M and S across all the biogeographic domains of Cameroon in Central Africa, in order to understand the role of chromosomal arrangements in ecological specialisation within and among molecular forms.

RESULTS

Species distribution modelling with presence-only data revealed differences in the ecological niche of both molecular forms and the sibling species, An. arabiensis. The fundamental environmental envelope of the two molecular forms, however, overlapped to a large extent in the rainforest, where they occurred in sympatry. The S form had the greatest niche breadth of all three taxa, whereas An. arabiensis and the M form had the smallest niche overlap. Correspondence analysis of M and S karyotypes confirmed that molecular forms shared similar combinations of chromosomal inversion arrangements in response to the eco-climatic gradient defining the main biogeographic domains occurring across Cameroon. Savanna karyotypes of M and S, however, segregated along the smaller-scale environmental gradient defined by the second ordination axis. Population structure analysis identified three chromosomal clusters, each containing a mixture of M and S specimens. In both M and S, alternative karyotypes were segregating in contrasted environments, in agreement with a strong ecological adaptive value of chromosomal inversions.

CONCLUSION

Our data suggest that inversions on the second chromosome of An. gambiae are not causal to the evolution of reproductive isolation between the M and S forms. Rather, they are involved in ecological specialization to a similar extent in both genetic backgrounds, and most probably predated lineage splitting between molecular forms. However, because chromosome-2 inversions promote ecological divergence, resulting in spatial and/or temporal isolation between ecotypes, they might favour mutations in other ecologically significant genes to accumulate in unlinked chromosomal regions. When such mutations occur in portions of the genome where recombination is suppressed, such as the pericentromeric regions known as speciation islands in An. gambiae, they would contribute further to the development of reproductive isolation.

VectorBase: a data resource for invertebrate vector genomics

VectorBase (http://www.vectorbase.org) is an NIAID-funded Bioinformatic Resource Center focused on invertebrate vectors of human pathogens. VectorBase annotates and curates vector genomes providing a web accessible integrated resource for the research community. Currently, VectorBase contains genome information for three mosquito species: Aedes aegypti, Anopheles gambiae and Culex quinquefasciatus, a body louse Pediculus humanus and a tick species Ixodes scapularis. Since our last report VectorBase has initiated a community annotation system, a microarray and gene expression repository and controlled vocabularies for anatomy and insecticide resistance. We have continued to develop both the software infrastructure and tools for interrogating the stored data.

Chromosomal plasticity and evolutionary potential in the malaria vector Anopheles gambiae sensu stricto: insights from three decades of rare paracentric inversions

Background

In the Anopheles gambiae complex, paracentric chromosomal inversions are non-randomly distributed along the complement: 18/31 (58%) of common polymorphic inversions are on chromosome arm 2R, which represents only ~30% of the complement. Moreover, in An. gambiae sensu stricto, 6/7 common polymorphic inversions occur on 2R. Most of these inversions are considered markers of ecological adaptation that increase the fitness of the carriers of alternative karyotypes in contrasting habitats. However, little is known about the evolutionary forces responsible for their origin and subsequent establishment in field populations.

Results

Here, we present data on 82 previously undescribed rare chromosomal inversions (RCIs) recorded during extensive field sampling in 16 African countries over a 30 year period, which may shed light on the dynamics of chromosomal plasticity in An. gambiae. We analyzed breakpoint distribution, length, and geographic distribution of RCIs, and compared these measures to those of the common inversions. We found that RCIs, like common inversions, are disproportionately clustered on 2R, which may indicate that this arm is especially prone to breakages. However, contrasting patterns were observed between the geographic distribution of common inversions and RCIs. RCIs were equally frequent across biomes and on both sides of the Great Rift Valley (GRV), whereas common inversions predominated in arid ecological settings and west of the GRV. Moreover, the distribution of RCI lengths followed a random pattern while common inversions were significantly less frequent at shorter lengths.

Conclusion

Because 17/82 (21%) RCIs were found repeatedly at very low frequencies – at the same sampling location in different years and/or in different sampling locations – we suggest that RCIs are subject mainly to drift under unperturbed ecological conditions. Nevertheless, RCIs may represent an important reservoir of genetic variation for An. gambiae in response to environmental changes, further testifying to the considerable evolutionary potential hidden within this pan-African malaria vector.

Localization of candidate regions maintaining a common polymorphic inversion (2La) in Anopheles gambiae

Chromosomal inversion polymorphisms are thought to play a role in adaptive divergence, but the genes conferring adaptive benefits remain elusive. Here we study 2La, a common polymorphic inversion in the African malaria vector Anopheles gambiae. The frequency of 2La varies clinally and seasonally in a pattern suggesting response to selection for aridity tolerance. By hybridizing genomic DNA from individual mosquitoes to oligonucleotide microarrays, we obtained a complete map of differentiation across the A. gambiae genome. Comparing mosquitoes homozygous for the 2La gene arrangement or its alternative (2L+^a), divergence was highest at loci within the rearranged region. In the 22 Mb included within alternative arrangements, two ∼1.5 Mb regions near but not adjacent to the breakpoints were identified as being significantly diverged, a conclusion validated by targeted sequencing. The persistent association of both regions with the 2La arrangement is highly unlikely given known recombination rates across the inversion in 2La heterozygotes, thus implicating selection on genes underlying these regions as factors responsible for the maintenance of 2La. Polymorphism and divergence data are consistent with a model in which the inversion is maintained by migration-selection balance between multiple alleles inside these regions, but further experiments will be needed to fully distinguish between the epistasis (coadaptation) and local adaptation models for the maintenance of 2La.

A chromosomal inversion occurs when part of the chromosome breaks, rotates 180 degrees, and rejoins the broken chromosome. The result is a chromosome carrying a segment whose gene order is reversed. Whereas the physical rearrangement itself may have no direct consequences on gene function, recombination between alleles in the rearranged and wild type segments is suppressed. If multiple alleles inside the inverted or original orientations are well adapted to contrasting environmental conditions, suppressed recombination provides a mechanism to keep beneficial allelic combinations from being shuffled between different genetic backgrounds. Working with wild populations of flies, Dobzhansky provided the first evidence that selection was key to maintaining inversion polymorphism. Subsequently, examples of inversion polymorphisms under selection in other organisms have been found, notably in the mosquito that transmits most cases of human malaria, Anopheles gambiae. However, the genes or gene regions conferring fitness advantages have yet to be discovered. In this study, the authors used modern genomics tools to map such regions in an inversion at an unprecedented level of detail, and show that these regions are likely to be responsible for the maintenance of the inversion polymorphism in natural populations. This study lays the groundwork for future efforts to identify the genes themselves and their role in adaptation.

Differential gene expression in incipient species of Anopheles gambiae

A speciation process is ongoing in the primary vector of malaria in Africa, Anopheles gambiae. Assortatively mating incipient species known as the M and S forms differentially exploit larval breeding sites associated with different ecological settings. However, some ongoing gene flow between M and S limits significant genomic differentiation mainly to small centromere-proximal regions on chromosomes X and 2L, termed 'speciation islands' with the expectation that they contain the genes responsible for reproductive isolation. As the speciation islands exhibit reduced recombination and low polymorphism, more detailed genetic analysis using fine-scale mapping is impractical. We measured global gene expression differences between M and S using oligonucleotide microarrays, with the goal of identifying candidate genes that could be involved in this ongoing speciation process. Gene expression profiles were examined in two independent colonies of both forms at each of three developmental periods of interest: fourth instar larvae, virgin females, and gravid females. Patterns were validated on a subset of genes using quantitative real-time reverse transcription polymerase chain reaction of RNA samples from laboratory colonies and wild mosquitoes collected from Cameroon and Burkina Faso. Considered across all three developmental periods, differentially expressed genes represented approximately 1-2% of all expressed genes. Although disproportionately represented in the X speciation island, the vast majority of genes were located outside any speciation island. Compared to samples from the other developmental periods, virgin females were characterized by more than twice as many differentially expressed genes, most notably those implicated in olfaction and potentially, mate recognition.

Genetic exchange in 2La inversion heterokaryotypes of Anopheles gambiae

In the malaria vector Anopheles gambiae, alternative arrangements of chromosome 2 (2La and 2L+(a)) vary in relative frequency along clines of aridity, suggesting the action of natural selection on targets within the inversion. Our long term goal of detecting such targets depends in part on the level of genetic exchange between arrangements. Accordingly, we estimated recombination rates on 2L from the backcross progeny of 2La/+(a) heterokaryotypes and as a control, from 2L+(a) homokaryotypes. In homokaryotypes, the recombination rate was uniform at ~2.0 centimorgans per megabase (cM/Mb). In heterokaryotypes, recombination within the rearranged region was reduced to < 0.5 cM/Mb, with slightly higher but nevertheless reduced levels (< 1.0 cM/Mb) flanking the rearrangement. Yet, gene exchange was recorded between nearly all markers, including those very near the distal inversion breakpoint. These results suggest that reduced recombination is a necessary but not sufficient mechanism for genetic isolation between alternative arrangements, and that the targets of natural selection can be identified against the different chromosomal backgrounds.

PCR-based karyotyping of Anopheles gambiae inversion 2Rj identifies the BAMAKO chromosomal form

Background

The malaria vector Anopheles gambiae is polymorphic for chromosomal inversions on the right arm of chromosome 2 that segregate nonrandomly between assortatively mating populations in West Africa. One such inversion, 2Rj, is associated with the BAMAKO chromosomal form endemic to southern Mali and northern Guinea Conakry near the Niger River. Although it exploits a unique ecology and both molecular and chromosomal data suggest reduced gene flow between BAMAKO and other A. gambiae populations, no molecular markers exist to identify this form.

Methods

To facilitate study of the BAMAKO form, a PCR assay for molecular karyotyping of 2Rj was developed based on sequences at the breakpoint junctions. The assay was extensively validated using more than 700 field specimens whose karyotypes were determined in parallel by cytogenetic and molecular methods. As inversion 2Rj also occurs in SAVANNA populations outside the geographic range of BAMAKO, samples were tested from Senegal, Cameroon and western Guinea Conakry as well as from Mali.

Results

In southern Mali, where 2Rj polymorphism in SAVANNA populations was very low and most of the 2Rj homozygotes were found in BAMAKO karyotypes, the molecular and cytogenetic methods were almost perfectly congruent. Elsewhere agreement between the methods was much poorer, as the molecular assay frequently misclassified 2Rj heterozygotes as 2R+^j standard homozygotes.

Conclusion

Molecular karyotyping of 2Rj is robust and accurate on 2R+^j standard and 2Rj inverted homozygotes. Therefore, the proposed approach overcomes the lack of a rapid tool for identifying the BAMAKO form across developmental stages and sexes, and opens new perspectives for the study of BAMAKO ecology and behaviour. On the other hand, the method should not be applied for molecular karyotyping of j-carriers within the SAVANNA chromosomal form.

Segmental duplication implicated in the genesis of inversion 2Rj of Anopheles gambiae

The malaria vector Anopheles gambiae maintains high levels of inversion polymorphism that facilitate its exploitation of diverse ecological settings across tropical Africa. Molecular characterization of inversion breakpoints is a first step toward understanding the processes that generate and maintain inversions. Here we focused on inversion 2Rj because of its association with the assortatively mating Bamako chromosomal form of An. gambiae, whose distinctive breeding sites are rock pools beside the Niger River in Mali and Guinea. Sequence and computational analysis of 2Rj revealed the same 14.6 kb insertion between both breakpoints, which occurred near but not within predicted genes. Each insertion consists of 5.3 kb terminal inverted repeat arms separated by a 4 kb spacer. The insertions lack coding capacity, and are comprised of degraded remnants of repetitive sequences including class I and II transposable elements. Because of their large size and patchwork composition, and as no other instances of these insertions were identified in the An. gambiae genome, they do not appear to be transposable elements. The 14.6 kb modules inserted at both 2Rj breakpoint junctions represent low copy repeats (LCRs, also called segmental duplications) that are strongly implicated in the recent (∼0.4N_e generations) origin of 2Rj. The LCRs contribute to further genome instability, as demonstrated by an imprecise excision event at the proximal breakpoint of 2Rj in field isolates.

Molecular karyotyping of the 2La inversion in Anopheles gambiae

The African malaria vector Anopheles gambiae is polymorphic for alternative arrangements on the left arm of chromosome 2 (2La and 2L+(a)) that are non-randomly distributed with respect to degree of aridity. Detailed studies on the ecological role of inversion 2La have been hindered by the technical demands of traditional karyotype analysis and by sex- and stage-specific limitations on the availability of polytene chromosomes favorable for analysis. Recent molecular characterization of both inversion breakpoints presented the opportunity to develop a polymerase chain reaction (PCR)-based method for karyotype analysis. Here we report the development of this molecular diagnostic assay and the results of extensive field validation. When tested on 765 An. gambiae specimens sampled across Africa, the molecular approach compared favorably with traditional cytologic methods, correctly scoring > 94% of these specimens. By providing ready access to the 2La karyotype, this tool lays groundwork for future studies of the ecological genomics of this medically important species.

Polymorphism at the defensin gene in the Anopheles gambiae complex: testing different selection hypotheses

Genetic variation in defensin, a gene encoding a major effector molecule of insects immune response was analyzed within and between populations of three members of the Anopheles gambiae complex. The species selected included the two anthropophilic species, An. gambiae and An. arabiensis and the most zoophilic species of the complex, An. quadriannulatus. The first species was represented by four populations spanning its extreme genetic and geographical ranges, whereas each of the other two species was represented by a single population. We found (i) reduced overall polymorphism in the mature peptide region and in the total coding region, together with specific reductions in rare and moderately frequent mutations (sites) in the coding region compared with non-coding regions, (ii) markedly reduced rate of non-synonymous diversity compared with synonymous variation in the mature peptide and virtually identical mature peptide across the three species, and (iii) increased divergence between species in the mature peptide together with reduced differentiation between populations of An. gambiae in the same DNA region. These patterns suggest a strong purifying selection on the mature peptide and probably the whole coding region. Because An. quadriannulatus is not exposed to human pathogens, identical mature peptide and similar pattern of polymorphism across species implies that human pathogens played no role as selective agents on this peptide.

VectorBase: a home for invertebrate vectors of human pathogens

VectorBase () is a web-accessible data repository for information about invertebrate vectors of human pathogens. VectorBase annotates and maintains vector genomes providing an integrated resource for the research community. Currently, VectorBase contains genome information for two organisms: Anopheles gambiae, a vector for the Plasmodium protozoan agent causing malaria, and Aedes aegypti, a vector for the flaviviral agents causing Yellow fever and Dengue fever.

Variation in recombination rate across the X chromosome of Anopheles gambiae

The M and S molecular forms of Anopheles gambiae are considered to be incipient species, despite residual gene exchange. Of the three small genome regions that are strongly differentiated between the molecular forms ("speciation islands"), two are located near centromeres, on the left arm of chromosome 2 and the X chromosome. To test the prediction of reduced recombination in these islands, we estimated recombination rates between microsatellite loci on the X chromosome using two M-form strains. Across most of the chromosome, recombination occurred at approximately 1 centimorgan per megabase (cM Mb(-1)), a value closely matching the genome-wide average estimated for A. gambiae and for other eukaryotes. Recombination was much higher at the telomeric end, > 7 cM Mb(-1). In the speciation island at the centromeric end, recombination was sharply reduced to approximately 0.2 cM Mb(-1), consistent with a role for reduced recombination in maintaining differentiation between nascent species despite gene flow.

Gene finding on the Y: fruitful strategy in Drosophila does not deliver in Anopheles

The Anopheles gambiae genome project yielded almost complete sequences for the autosomes and for a large part of the X chromosome, however, no information for the Y chromosome was obtained. Yet, by design, fragmented Y chromosome sequences should be present in the resulting assembly. Here we report the search for Anopheles Y chromosome genes using a strategy successfully applied for identification of Y genes in Drosophila. A complete set of the unmapped scaffolds was targeted in a broad TBLASTN search using both A. gambiae predicted genes and all proteins from nr database as query sequences. After filtering of the BLAST report, we selected 181 scaffolds possibly containing fragments of Y chromosome genes to experimentally test their Y-linkage. Surprisingly, none of the tested sequences appeared to originate from the Y chromosome. Several factors could account for the failure to detect Y genes, including their different organization in A. gambiae compared to Drosophila and the suboptimal quality of the assembly and annotation of the Anopheles genome. Regardless of the cause, our results illuminate problems associated with the genome analysis of outbred organisms.

Rangewide population genetic structure of the African malaria vector Anopheles funestus

Anopheles funestus is a primary vector of malaria in Africa south of the Sahara. We assessed its rangewide population genetic structure based on samples from 11 countries, using 10 physically mapped microsatellite loci, two per autosome arm and the X (N = 548), and 834 bp of the mitochondrial ND5 gene (N = 470). On the basis of microsatellite allele frequencies, we found three subdivisions: eastern (coastal Tanzania, Malawi, Mozambique and Madagascar), western (Burkina Faso, Mali, Nigeria and western Kenya), and central (Gabon, coastal Angola). A. funestus from the southwest of Uganda had affinities to all three subdivisions. Mitochondrial DNA (mtDNA) corroborated this structure, although mtDNA gene trees showed less resolution. The eastern subdivision had significantly lower diversity, similar to the pattern found in the codistributed malaria vector Anopheles gambiae. This suggests that both species have responded to common geographic and/or climatic constraints. The western division showed signatures of population expansion encompassing Kenya west of the Rift Valley through Burkina Faso and Mali. This pattern also bears similarity to A. gambiae, and may reflect a common response to expanding human populations following the development of agriculture. Due to the presumed recent population expansion, the correlation between genetic and geographic distance was weak. Mitochondrial DNA revealed further cryptic subdivision in A. funestus, not detected in the nuclear genome. Mozambique and Madagascar samples contained two mtDNA lineages, designated clade I and clade II, that were separated by two fixed differences and an average of 2% divergence, which implies that they have evolved independently for approximately 1 million years. Clade I was found in all 11 locations, whereas clade II was sampled only on Madagascar and Mozambique. We suggest that the latter clade may represent mtDNA capture by A. funestus, resulting from historical gene flow either among previously isolated and divergent populations or with a related species.

Divergence with gene flow in Anopheles funestus from the Sudan Savanna of Burkina Faso, West Africa

Anopheles funestus is a major vector of malaria across Africa. Understanding its complex and nonequilibrium population genetic structure is an important challenge that must be overcome before vector populations can be successfully perturbed for malaria control. Here we examine the role of chromosomal inversions in structuring genetic variation and facilitating divergence in Burkina Faso, West Africa, where two incipient species (chromosomal forms) of A. funestus, defined principally by rearrangements of chromosome 3R, have been hypothesized. Sampling across an approximately 300-km east-west transect largely contained within the Sudan-Savanna ecoclimatic zone, we analyzed chromosomal inversions, 16 microsatellite loci distributed genomewide, and 834 bp of the mtDNA ND5 gene. Both molecular markers revealed high genetic diversity, nearly all of which was accounted for by within-population differences among individuals, owing to recent population expansion. Across the study area there was no correlation between genetic and geographic distance. Significant genetic differentiation found between chromosomal forms on the basis of microsatellites was not genomewide but could be explained by chromosome 3R alone on the basis of loci inside and near inversions. These data are not compatible with complete reproductive isolation but are consistent with differential introgression and sympatric divergence between the chromosomal forms, facilitated by chromosome 3R inversions.

Breakpoint structure reveals the unique origin of an interspecific chromosomal inversion (2La) in the Anopheles gambiae complex

Paracentric chromosomal inversions are major architects of organismal evolution and have been associated with adaptations relevant to malaria transmission in anopheline mosquitoes. The processes responsible for their origin and maintenance, still poorly understood, can be illuminated by analysis of inversion breakpoint sequences. Here, we report the breakpoint structure of chromosomal inversion 2La from the principal malaria vector Anopheles gambiae and its relatives in the A. gambiae complex. The distal and proximal breakpoints of the standard (2L+a) arrangement contain gene duplications: full-length genes and their truncated copies at opposite ends. Intact genes without pseudogene copies in the alternative arrangement (2La) imply that 2L+a is derived and was viable despite damage to genes, because duplication preserved gene function. A unique origin for the interspecific 2La inversion was challenged previously by indirect genetic evidence, but breakpoint sequences determined from members of the A. gambiae complex strongly suggest their descent from a single event. The derived position of 2L+a, long considered ancestral in this medically important group, has significant implications for the phylogenetic history and the evolution of vectorial capacity in the A. gambiae complex.

Analysis of the complete mitochondrial DNA from Anopheles funestus: an improved dipteran mitochondrial genome annotation and a temporal dimension of mosquito evolution

Virtually no information regarding timing of deep lineage divergences within mosquito family (Culicidae) exists, which poses an important problem in the postgenomic era. To address this issue, the complete 15,354 bp mitochondrial genome of Anopheles funestus was assembled from both mtDNA and cDNA sequences generated from transcripts of the mtDNA-encoded protein and rRNA genes. Analysis of the transcript information allowed an improved genome annotation, revealing that the translation initiation codon for the cox1 gene is TCG, rather than atypical, longer codons proposed in several other insects. The 5'ends of nad1 and nad5 transcripts begin with TTG and GTG triplets, respectively, which apparently serve as the translation initiators for those genes. We used all the A. funestus mtDNA gene sequences and three other publicly available mosquito mtDNA genomes for the estimation of divergence time points within Culicidae. The maximum likelihood date estimates for the splits between Anopheles and Aedes (approximately 145-200 Mya), between Anopheles subgenera Cellia and Anopheles (approximately 90-106 Mya), and between lineages within subgenus Anopheles (approximately 70-85 Mya) inferred from protein-coding genes are roughly twice as high as the dates based on RNA gene sequences. Although existing evidence does not unequivocally favor one of the alternatives, fossil-based predictions of the age of the family Culicidae are in better agreement with dates inferred from protein-coding genes.

Genetic structure of Anopheles gambiae populations on islands in northwestern Lake Victoria, Uganda

Background

Alternative means of malaria control are urgently needed. Evaluating the effectiveness of measures that involve genetic manipulation of vector populations will be facilitated by identifying small, genetically isolated vector populations. The study was designed to use variation in microsatellite markers to look at genetic structure across four Lake Victoria islands and two surrounding mainland populations and for evidence of any restriction to free gene flow.

Methods

Four Islands (from 20–50 km apart) and two surrounding mainland populations (96 km apart) were studied. Samples of indoor resting adult mosquitoes, collected over two consecutive years, were genotyped at microsatellite loci distributed broadly throughout the genome and analysed for genetic structure, effective migration (Nem) and effective population size (Ne).

Results

Ne estimates showed island populations to consist of smaller demes compared to the mainland ones. Most populations were significantly differentiated geographically, and from one year to the other. Average geographic pair-wise FST ranged from 0.014–0.105 and several pairs of populations had Ne m < 3. The loci showed broad heterogeneity at capturing or estimating population differences.

Conclusion

These island populations are significantly genetically differentiated. Differences reoccurred over the study period, between the two mainland populations and between each other. This appears to be the product of their separation by water, dynamics of small populations and local adaptation. With further characterisation these islands could become possible sites for applying measures evaluating effectiveness of control by genetic manipulation.

Chromosomal evidence of incipient speciation in the Afrotropical malaria mosquito Anopheles funestus

The analysis of chromosomal polymorphism of paracentric inversions in anopheline mosquitoes has often been instrumental to the discovery of sibling species complexes and intraspecific genetic heterogeneities associated with incipient speciation processes. To investigate the population structure of Anopheles funestus Giles (Diptera: Culicidae), one of the three most important vectors of human malaria in sub-Saharan Africa, a three-year survey of chromosomal polymorphism was carried out on 4,638 karyotyped females collected indoors and outdoors from two villages of central Burkina Faso. Large and temporally stable departures from Hardy-Weinberg equilibrium due to significant deficits of heterokaryotypes were found irrespective of the place of capture, and of the spatial and temporal units chosen for the analysis. Significant linkage disequilibrium was observed among inversion systems on independently assorting chromosomal arms, indicating the existence of assortative mating phenomena. Results were consistent with the existence of two chromosomal forms characterized by contrasting degrees of inversion polymorphism maintained by limitations to gene flow. This hypothesis was supported by the reestablishment of Hardy-Weinberg and linkage equilibria when individual specimens were assigned to each chromosomal form according to two different algorithms. This pattern of chromosomal variability is suggestive of an incipient speciation process in An. funestus populations from Burkina Faso.

Centromere-proximal differentiation and speciation in Anopheles gambiae

The M and S molecular forms of Anopheles gambiae are undergoing speciation as they adapt to heterogeneities in the environment, spreading malaria in the process. We hypothesized that their divergence despite gene flow is facilitated by reduced recombination at the centromeric (proximal) end of the X chromosome. We sequenced introns from 22 X chromosome genes in M and S from two locations of West Africa where the forms are sympatric. Generally, in both forms nucleotide diversity was high distally, lower proximally, and very low nearest the centromere. Conversely, differentiation between the forms was virtually zero distally and very high proximally. Pairwise comparisons to a close relative, the sibling species Anopheles arabiensis, demonstrated uniformly high divergence regardless of position along the X chromosome, suggesting that this pattern is not purely mechanical. Instead, the pattern observed for M and S suggests the action of divergent natural selection countering gene flow only at the proximal end of the X chromosome, where recombination is reduced. Comparison of sites with fixed differences between M and S to the corresponding sites in A. arabiensis revealed that derived substitutions had been fixed in both forms, further supporting the hypothesis that both have been under selection. These derived substitutions are fixed in the two West African samples and in samples of S from western and coastal Kenya, suggesting that selection occurred before the forms expanded to their current ranges. Our findings are consistent with a role for suppressed genetic recombination in speciation of A. gambiae.

SINE insertion polymorphism on the X chromosome differentiates Anopheles gambiae molecular forms

Polymorphic SINE insertions can be useful markers for assessing population structure and differentiation. Maque is a family of SINE elements which, based on bioinformatic analysis, was suggested to have been active recently in Anopheles gambiae, the major vector of malaria. Here, we report the development of polymorphic Maque insertions as population genetic markers in A. gambiae, and the use of these markers to better characterize divergence on the X chromosome between A. gambiae M and S molecular forms in populations from Burkina Faso and Mali. Our data are consistent with the recent activity of Maque. Phylogenetic analysis suggests that at least two recently active lineages may have a role in mediating genome evolution. We found differences in element insertion frequency and sequence between the M and S populations analysed. Significant differentiation was observed between these two groups across a 6 Mb region at the proximal (centromeric) end of the X chromosome. Locus-specific F(ST) values ranged from 0.14 to 1.00 in this region, yet were not significantly different from zero in more distal locations on the X chromosome; the trend was consistent in populations from both geographical locales suggesting that differentiation is not due to local adaptation. Strong differentiation between M and S at the proximal end of the X chromosome, but not outside this region, suggests the action of selection counteracting limited gene flow between these taxa and supports their characterization as incipient species.

Molecular differentiation between chromosomally defined incipient species of Anopheles funestus

Anopheles funestus Giles is one of the most important vectors of malaria in sub-Saharan Africa. The population structure of this mosquito in Burkina Faso, West Africa based on chromosomal inversion data led to the description of two chromosomal forms, Kiribina and Folonzo. Because both forms co-occur in the same locales yet differ significantly, both in the frequency of inverted arrangements on chromosome arms 3R and 2R and in vectorial capacity, they were hypothesized to be emerging species with at least partial barriers to gene flow. This hypothesis would be strengthened by molecular evidence of differentiation between Kiribina and Folonzo at loci outside chromosomal inversions. We surveyed molecular variation in sympatric populations of the two forms using sequences from the mitochondrial ND5 gene and genotypes at sixteen microsatellite loci distributed across the genome. Both classes of marker revealed slight but significant differentiation between the two forms (mtDNA F(ST) = 0.023, P < 0.001; microsatellite F(ST) = 0.004, P < 0.001; R(st) = 0.009, P = 0.002). Locus-by-locus analysis of the microsatellite data showed that significant differentiation was not genome-wide, but could be attributed to five loci on chromosome 3R (F(ST) = 0.010, P < 0.001; R(st) = 0.016, P = 0.002). Importantly, three of these loci are outside of, and in linkage equilibrium with, chromosomal inversions, suggesting that differentiation between chromosomal forms extends beyond the inversions themselves. The slight overall degree of differentiation indicated by both marker classes is likely an underestimate because of recent population expansion inferred for both Folonzo and Kiribina. The molecular evidence from this study is consistent with the hypothesis of incipient speciation between Kiribina and Folonzo.

The mosquito genome: organization, evolution and manipulation

Apart from the genetic flexibility of the vectors, impediments to the control of vector-borne diseases include the rapid spread of drug resistance throughout parasite populations, the increasing movement of people to and from disease-endemic regions and the limited funds and public health infrastructures of most developing countries. The widely used residual insecticides and antiparasitic drugs have been inadequate solutions to the problem of vector-borne disease control. New approaches are needed. The enormous impact of recent developments in molecular genetics on the understanding of basic biology and human disease has stimulated a re-examination of the prospects for genetic manipulation of vector populations as a means for reducing or eliminating vector-borne diseases, especially malarial. Although control scenarios that exploit this technology may never be realized, Nora Besansky and Frank Collins emphasize that the increase in knowledge of basic mosquito biology on which these ideas depend will inevitably stimulate novel approaches to the control of mosquito-borne diseases.

Sex-linked differentiation between incipient species of Anopheles gambiae

Emerging species within the primary malaria vector Anopheles gambiae show different ecological preferences and significant prezygotic reproductive isolation. They are defined by fixed sequence differences in X-linked rDNA, but most previous studies have failed to detect large and significant differentiation between these taxa elsewhere in the genome, except at two other loci on the X chromosome near the rDNA locus. Hypothesizing that this pericentromeric region of the X chromosome may be accumulating differences faster than other regions of the genome, we explored the pattern and extent of differentiation between A. gambiae incipient species and a sibling species, A. arabiensis, from Burkina Faso, West Africa, at 17 microsatellite loci spanning the X chromosome. Interspecific differentiation was large and significant across the entire X chromosome. Among A. gambiae incipient species, we found some of the highest levels of differentiation recorded in a large region including eight independent loci near the centromere of the X chromosome. Outside of this region, no significant differentiation was detected. This pattern suggests that selection is playing a role in the emergence of A. gambiae incipient species. This process, associated with efficient exploitation of anthropogenic modifications to the environment, has public health implications as it fosters the spread of malaria transmission both spatially and temporally.

Isolation and characterization of Y chromosome sequences from the African malaria mosquito Anopheles gambiae

The karyotype of the African malaria mosquito Anopheles gambiae contains two pairs of autosomes and a pair of sex chromosomes. The Y chromosome, constituting approximately 10% of the genome, remains virtually unexplored, despite the recent completion of the A. gambiae genome project. Here we report the identification and characterization of Y chromosome sequences of total length approaching 150 kb. We developed 11 Y-specific PCR markers that consistently yielded male-specific products in specimens from both laboratory colony and natural populations. The markers are characterized by low sequence polymorphism in samples collected across Africa and by presence in more than one copy on the Y. Screening of the A. gambiae BAC library using these markers allowed detection of 90 Y-linked BAC clones. Analysis of the BAC sequences and other Y-derived fragments showed massive accumulation of a few transposable elements. Nevertheless, more complex sequences are apparently present on the Y; these include portions of an approximately 48-kb-long unmapped AAAB01008227 scaffold from the whole genome shotgun assembly. Anopheles Y appears not to harbor any of the genes identified in Drosophila Y. However, experiments suggest that one of the ORFs from the AAAB01008227 scaffold represents a fragment of a gene with male-specific expression.

Satellite DNA from the Y chromosome of the malaria vector Anopheles gambiae

Satellite DNA is an enigmatic component of genomic DNA with unclear function that has been regarded as "junk." Yet, persistence of these tandem highly repetitive sequences in heterochromatic regions of most eukaryotic chromosomes attests to their importance in the genome. We explored the Anopheles gambiae genome for the presence of satellite repeats and identified 12 novel satellite DNA families. Certain families were found in close juxtaposition within the genome. Six satellites, falling into two evolutionarily linked groups, were investigated in detail. Four of them were experimentally confirmed to be linked to the Y chromosome, whereas their relatives occupy centromeric regions of either the X chromosome or the autosomes. A complex evolutionary pattern was revealed among the AgY477-like satellites, suggesting their rapid turnover in the A. gambiae complex and, potentially, recombination between sex chromosomes. The substitution pattern suggested rolling circle replication as an array expansion mechanism in the Y-linked 53-bp satellite families. Despite residing in different portions of the genome, the 53-bp satellites share the same monomer lengths, apparently maintained by molecular drive or structural constraints. Potential functional centromeric DNA structures, consisting of twofold dyad symmetries flanked by a common sequence motif, have been identified in both satellite groups.

Revisiting the role of introgression vs shared ancestral polymorphisms as key processes shaping genetic diversity in the recently separated sibling species of the Anopheles gambiae complex

The role of interspecific hybridisation in the evolution of pest species is poorly understood. In mosquito disease vectors this is of particular importance due to the evolution of insecticide resistance and the proposed release of transgenic strains that are refractory to the malaria parasite. In this study, we apply population genetic methods in a novel manner to determine whether mitochondrial DNA sequences have introgressed between the closely related African malaria vectors Anopheles gambiae and A. arabiensis. Our results suggest that speciation was geologically recent and ancestral haplotypes at the ND5 locus are retained in both species. In addition, comparing haplotype frequencies in allopatric and sympatric populations, suggest locale specific unidirectional introgression of mitochondria from A. arabiensis into A. gambiae.

Dynamics of the pyrethroid knockdown resistance allele in western Kenyan populations of Anopheles gambiae in response to insecticide-treated bed net trials

Permethrin and DDT resistance in Anopheles gambiae s.s. associated with a leucine-serine knockdown resistance (kdr) mutation in the voltage-gated sodium channel gene was discovered recently in western Kenya where a large scale permethrin-impregnated bed net (ITN) program has been implemented. Collections of An. gambiae s.l. were made from intervention and control villages prior to and after onset of the program. The kdr genotypes were determined using allele-specific polymerase chain reaction diagnostic tests. In An. gambiae s.s., the frequency of the kdr mutation prior to ITN introduction was approximately 3-4% in western Kenya and zero in samples from the coast. After ITN introduction, the kdr mutation increased in ITN and neighboring villages from approximately 4% to approximately 8%, but remained unchanged in villages at least 20 km distant and was not detected in coastal Kenya. The identical leucine-serine mutation was found in a single An. arabiensis individual among 658 tested. The leucine-phenylalanine kdr mutation common in west African An. gambiae populations was not detected in An. gambiae s.l. from Kenya. Implications for the population structure and control of An. gambiae are discussed.

A microsatellite map of the African human malaria vector Anopheles funestus

Microsatellite markers and chromosomal inversion polymorphisms are useful genetic markers for determining population structure in Anopheline mosquitoes. In Anopheles funestus (2N = 6), only chromosome arms 2R, 3R, and 3L are known to carry polymorphic inversions. The physical location of microsatellite markers with respect to polymorphic inversions is potentially important information for interpreting population genetic structure, yet none of the available marker sets have been physically mapped in this species. Accordingly, we mapped 32 polymorphic A. funestus microsatellite markers to the polytene chromosomes using fluorescent in situ hybridization (FISH) and identified 16 markers outside of known polymorphic inversions. Here we provide an integrated polytene chromosome map for A. funestus that includes the breakpoints of all known polymorphic inversions as well as the physical locations of microsatellite loci developed to date. Based on this map, we suggest a standard set of 16 polymorphic microsatellite markers that are distributed evenly across the chromosome complement, occur predominantly outside of inversions, and amplify reliably. Adoption of this set by researchers working in different regions of Africa will facilitate metapopulation analyses of this primary malaria vector.

Semipermeable species boundaries between Anopheles gambiae and Anopheles arabiensis: evidence from multilocus DNA sequence variation

Attempts to reconstruct the phylogenetic history of the Anopheles gambiae cryptic species complex have yielded strongly conflicting results. In particular, An. gambiae, the primary African malaria vector, is variously placed as a sister taxon to either Anopheles arabiensis or Anopheles merus. The recent divergence times for members of this complex complicate phylogenetic analysis, making it difficult to unambiguously implicate interspecific gene flow, versus retained ancestral polymorphism, as the source of conflict. Using sequences at four unlinked loci, which were determined from multiple specimens within each of five species in the complex, we found contrasting patterns of sequence divergence between the X chromosome and the autosomes. The isolation model of speciation assumes a lack of gene flow between species since their separation. This model could not be rejected for An. gambiae and An. arabiensis, although the data fit the model poorly. On the other hand, evidence from gene trees supports genetic introgression of chromosome 2 inversions between An. gambiae and An. arabiensis, and also points to more broad scale genetic exchange of autosomal sequences between this species pair. That such exchange has been relatively recent is suggested not only by the lack of fixed differences at three autosomal loci but also by the sharing of full haplotypes at two of the three loci, which is in contrast to several fixed differences and considerably deeper divergence on the X. The proposed acquisition by An. gambiae of sequences from the more arid-adapted An. arabiensis may have contributed to the spread and ecological dominance of this malaria vector.

Inversions and gene order shuffling in Anopheles gambiae and A. funestus

In tropical Africa, Anopheles funestus is one of the three most important malaria vectors. We physically mapped 157 A. funestus complementary DNAs (cDNAs) to the polytene chromosomes of this species. Sequences of the cDNAs were mapped in silico to the A. gambiae genome as part of a comparative genomic study of synteny, gene order, and sequence conservation between A. funestus and A. gambiae. These species are in the same subgenus and diverged about as recently as humans and chimpanzees. Despite nearly perfect preservation of synteny, we found substantial shuffling of gene order along corresponding chromosome arms. Since the divergence of these species, at least 70 chromosomal inversions have been fixed, the highest rate of rearrangement of any eukaryote studied to date. The high incidence of paracentric inversions and limited colinearity suggests that locating genes in one anopheline species based on gene order in another may be limited to closely related taxa.

Speciation within Anopheles gambiae--the glass is half full

Restrictions to gene flow among molecular forms of the mosquito Anopheles gambiae sensu stricto reveal an ongoing speciation process affecting the epidemiology of malaria in sub-Saharan Africa.

Molecular systematics of Anopheles: from subgenera to subpopulations

The century-old discovery of the role of Anopheles in human malaria transmission precipitated intense study of this genus at the alpha taxonomy level, but until recently little attention was focused on the systematics of this group. The application of molecular approaches to systematic problems ranging from subgeneric relationships to relationships at and below the species level is helping to address questions such as anopheline phylogenetics and biogeography, the nature of species boundaries, and the forces that have structured genetic variation within species. Current knowledge in these areas is reviewed, with an emphasis on the Anopheles gambiae model. The recent publication of the genome of this anopheline mosquito will have a profound impact on inquiries at all taxonomic levels, supplying better tools for estimating phylogeny and population structure in the short term, and ultimately allowing the identification of genes and/or regulatory networks underlying ecological differentiation, speciation, and vectorial capacity.

Structure and evolution of mtanga, a retrotransposon actively expressed on the Y chromosome of the African malaria vector Anopheles gambiae

Here we report the discovery of a novel family of long terminal repeat (LTR)-retrotransposons designated mtanga-Y, specific to the Y chromosome of the African malaria vector, Anopheles gambiae. mtanga-Y elements represent the first Y-linked sequences and the first members of the Ty1-copia superfamily of retrotransposons described from this mosquito. Analysis of a full-length 4,284-bp element revealed the presence of two intact overlapping open reading frames bounded by LTRs of 119 bp. Evidence suggests that the elements are capable of retrotransposition, as transcripts and potential replication intermediates (one-LTR circles) were detected. However, the approximately 12 copies of mtanga-Y appear to be clustered rather than dispersed on the Y chromosome. Absent from the Y chromosome of four sibling species (A. arabiensis, A. quadriannulatus, A. melas, and A. merus), similar, but often defective, mtanga elements are present elsewhere in these genomes, as well as in A. gambiae. These data are consistent with a relatively recent invasion of the A. gambiae Y chromosome by an intact element. The presence of functional mtanga-Y elements suggests that the Y chromosome may be a source, not just a sink, for retrotransposons.

Toward understanding Anophelinae (Diptera, Culicidae) phylogeny: insights from nuclear single-copy genes and the weight of evidence

A phylogeny of the mosquito subfamily Anophelinae was inferred from fragments of two protein-coding nuclear genes, G6pd (462 bp) and white (801 bp), and from a combined data set (2,136 bp) that included a portion of the mitochondrial gene ND5 and the D2 region of the ribosomal 28S gene. Sixteen species from all three anopheline genera and six Anopheles subgenera were sampled, along with six species of other mosquitoes used as an outgroup. Each of four genes analyzed individually recovered the same well-supported clades; topological incongruence was limited to unsupported or poorly supported nodes. As assessed by the incongruence length difference test, most of the conflicting signal was contributed by third codon positions. Strong structural constraints, as observed in white and G6pd, apparently had little impact on phylogenetic inference. Compared with the other genes, white provided a superior source of phylogenetic information. However, white appears to have experienced accelerated rates of evolution in few lineages, the affinities of which are therefore suspect. In combined analyses, most of the inferred relationship were well-supported and in agreement with previous studies: monophyly of Anophelinae, basal position of Chagasia, monophyly of Anopheles subgenera, and subgenera Nyssorhynchus + Kerteszia as sister taxa. The results suggested also monophyletic origin of subgenera Cellia + Anopheles, and the white gene analysis supported genus Bironella as a sister taxon to Anopheles. The present data and other available evidence suggest a South American origin of Anophelinae, probably in the Mesozoic; a rapid diversification of Bironella and basal subgeneric lineages of Anopheles, potentially associated with the breakup of Gondwanaland; and a relatively recent and rapid dispersion of subgenus Anopheles.

Evolution of mitochondrial and ribosomal gene sequences in anophelinae (Diptera: Culicidae): implications for phylogeny reconstruction

In this study, two mitochondrial genes, cyt b and ND5, and the D2 expansion segment of the 28S nuclear ribosomal gene were used to reconstruct a phylogeny of the mosquito subfamily Anophelinae. The ingroup consisted of all three genera of Anophelinae and five of six subgenera of Anopheles. Six genera of Culicinae were used as the outgroup. Extreme conservation at the protein level coupled with rapid saturation of synonymous positions probably accounted for the lack of meaningful phylogenetic signal in the cyt b gene. In contrast, abundant variation at all codon positions of the ND5 gene allowed recovery of the basal and most of the recent relationships. Phylogenetic analysis of D2 produced results consistent with those of ND5. Combined analysis indicated well-supported monophyletic Anophelinae (with Chagasia basal), Anopheles + Bironella, and subgeneric clades within the genus Anopheles. Moreover, subgenera Nyssorhynchus and Kerteszia were supported as a monophyletic lineage. The Kishino-Hasegawa test could not reject the monophyly of Anopheles, whereas the recently proposed hypothesis of close affinity of Bironella to the subgenus Anopheles was rejected by the analyses of ND5 and combined data sets. The lack of resolution of Bironella and Anopheles clades, or basal relationships among subgeneric clades within Anopheles, suggests their rapid diversification. Recovery of relationships consistent with morphology and previous molecular studies provides evidence of substantial phylogenetic signal in D2 and ND5 genes at levels of divergence from closely related species to subfamily in mosquitoes.

Patterns of DNA sequence variation in chromosomally recognized taxa of Anopheles gambiae: evidence from rDNA and single-copy loci

Patterns of DNA sequence variation in the ribosomal DNA (rDNA) second internal transcribed spacer (ITS2) and five unlinked single-copy nuclear loci were examined for evidence of reproductive isolation among four chromosomally recognized taxa of Anopheles gambiae from West Africa: Savanna, Bamako, Mopti and Forest, as well as sibling species An. arabiensis and An. merus. Included among the single-copy loci were three sequence-tagged random amplified polymorphic DNA (RAPD) loci, two of which (R15 and R37) had been reported as discriminating between Mopti and other chromosomal forms. Each of the five single-copy sequences were highly polymorphic in most samples. However, the R15 and R37 loci had no diagnostic value, and therefore are not recommended as tools in recognition of field-collected An. gambiae chromosomal forms. Although pairwise comparisons between species generally revealed significant levels of differentiation at all five loci, variation was not partitioned by chromosomal form within An. gambiae at any single-copy locus examined. The few exceptions to these trends appear related to a location either inside or nearby chromosomal inversions. At the tryptophan oxygenase locus inside inversion 2Rb, variation was structured only by inversion orientation and not by taxonomic designation even between An. gambiae and An. arabiensis, providing the first molecular evidence that the 2Rb inversion was transferred between species by introgressive hybridization. By contrast, the rDNA showed fixed differences between species and a difference diagnostic for Mopti, consistent with effective, if not complete, reproductive isolation. The apparent disagreement between the data from this locus and multiple single-copy loci within An. gambiae may be explained by the much lower effective population size of rDNA, owing to concerted evolution, which confers increased sensitivity at much shorter divergence times. Taken together with the accompanying reports by della Torre et al. (2001), Favia et al. (2001) and Gentile et al. (2001), our data suggest that neutral molecular markers may not have the sensitivity required to detect isolation between these recently established taxa.

Ribosomal DNA internal transcribed spacer (ITS2) sequences differentiate Anopheles funestus and An. rivulorum, and uncover a cryptic taxon

Differentiation among the closely related Afrotropical species comprising the Funestus Group is difficult by traditional taxonomic measures. Anopheles rivulorum is the second most abundant and widespread species in the Funestus Group, and is occasionally collected indoors along with the dominant member and major malaria vector, An. funestus. The prospect of misidentification of An. rivulorum as An. funestus prompted the development of a rapid, polymerase chain reaction (PCR)-based method for identifying these two species. The ribosomal internal transcribed spacer 2 (ITS2) was amplified from thirty-five specimens of An. rivulorum collected from the extremes of its range: Eastern Africa (Kenya), Southern Africa (South Africa) and Western Africa (Burkina Faso). The ITS2 region of An. rivulorum ( approximately 380 bp) is sufficiently different in size from the ITS2 of An. funestus ( approximately 700 bp) that these species can be distinguished by agarose gel electrophoresis of PCR products without further manipulation. Comparison of the An. rivulorum and An. funestus ITS2 nucleotide sequences revealed such extensive divergence that meaningful alignment was impossible, except for a 25 bp island near the 5' end. Intraspecific sequence comparisons revealed no variation among An. rivulorum individuals collected from the same country. However, sequence divergence was 2% between specimens from South Africa and Kenya, and nearly tenfold higher ( approximately 19%) between specimens from Burkina Faso and either South Africa or Kenya, an unprecedented level of intraspecific ITS2 divergence in Anopheles. Taken together, these data suggest that the Burkina Faso sample is not An. rivulorum, but rather a cryptic taxon within the Funestus Group.

The Rift Valley complex as a barrier to gene flow for Anopheles gambiae in Kenya: the mtDNA perspective

Descriptions of A. gambiae population structure based on microsatellite loci and mitochondrial DNA (mtDNA) were incongruent. High differentiation of populations was measured across the Rift Valley by microsatellites, but no differentiation was detected based on mtDNA. To resolve this conflict, we compared the old data to new mtDNA data using the same specimen previously genotyped in microsatellite loci. Analysis of a larger number of mtDNA sequences resulted in high and significant differentiation between populations across the Rift Valley. We developed a method to assess whether differentiation across the Rift Valley was generated by pure drift rather than mutation-drift, based on DNA sequence data. Applying this method to the mtDNA data suggested that pure drift was the primary force generating differentiation between the populations across the Rift, while mutation-drift generated differentiation across the continent. Given adequate sample size, mtDNA provided congruent results with microsatellite loci.