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.