In silico karyotyping of chromosomally polymorphic malaria mosquitoes in the Anopheles gambiae complex

Chromosomal inversion polymorphisms play an important role in adaptation to environmental heterogeneities. For mosquito species in the Anopheles gambiae complex that are significant vectors of human malaria, paracentric inversion polymorphisms are abundant and are associated with ecologically and epidemiologically important phenotypes. Improved understanding of these traits relies on determining mosquito karyotype, which currently depends upon laborious cytogenetic methods whose application is limited both by the requirement for specialized expertise and for properly preserved adult females at specific gonotrophic stages. To overcome this limitation, we developed sets of tag single nucleotide polymorphisms (SNPs) inside inversions whose biallelic genotype is strongly correlated with inversion genotype. We leveraged 1,347 fully sequenced An. gambiae and Anopheles coluzzii genomes in the Ag1000G database of natural variation. Beginning with principal components analysis (PCA) of population samples, applied to windows of the genome containing individual chromosomal rearrangements, we classified samples into three inversion genotypes, distinguishing homozygous inverted and homozygous uninverted groups by inclusion of the small subset of specimens in Ag1000G that are associated with cytogenetic metadata. We then assessed the correlation between candidate tag SNP genotypes and PCA-based inversion genotypes in our training sets, selecting those candidates with >80% agreement. Our initial tests both in held-back validation samples from Ag1000G and in data independent of Ag1000G suggest that when used for in silico inversion genotyping of sequenced mosquitoes, these tags perform better than traditional cytogenetics, even for specimens where only a small subset of the tag SNPs can be successfully ascertained

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 4(th) 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

Habitat segregation and ecological character displacement in cryptic African malaria mosquitoes

This is an open access article under the terms of the Creative Commons Attribution License.-- et al.Understanding how divergent selection generates adaptive phenotypic and population diversification provides a mechanistic explanation of speciation in recently separated species pairs. Towards this goal, we sought ecological gradients of divergence between the cryptic malaria vectors Anopheles coluzzii and An. gambiae and then looked for a physiological trait that may underlie such divergence. Using a large set of occurrence records and eco-geographic information, we built a distribution model to predict the predominance of the two species across their range of sympatry. Our model predicts two novel gradients along which the species segregate: distance from the coastline and altitude. Anopheles coluzzii showed a 'bimodal' distribution, predominating in xeric West African savannas and along the western coastal fringe of Africa. To test whether differences in salinity tolerance underlie this habitat segregation, we assessed the acute dose-mortality response to salinity of thirty-two larval populations from Central Africa. In agreement with its coastal predominance, Anopheles coluzzii was overall more tolerant than An. gambiae. Salinity tolerance of both species, however, converged in urban localities, presumably reflecting an adaptive response to osmotic stress from anthropogenic pollutants. When comparing degree of tolerance in conjunction with levels of syntopy, we found evidence of character displacement in this trait.Funding was provided by the Institut de Recherche pour le Dèveloppement, the National Institutes of Health (Projects R01 AI63508, and R01 AI76584 to NJB) and the Wellcome Trust (Intermediate Fellowship in Public Health and Tropical Medicine N° WT086423MA to CAN). PA was funded by the Spanish Ministerio de Economía y Competitividad (MINECO) and Universidad de Castilla-La Mancha (UCLM) through a ‘Ramón y Cajal’ contract(RYC-2012-11970). INAM was supported by a fellowship from the Universidad Nacional de Guinea Ecuatorial.Peer Reviewe

Highly specific PCR-RFLP assays for karyotyping the widespread 2Rb inversion in malaria vectors of the Anopheles gambiae complex

Background: Chromosomal inversion polymorphisms play a role in adaptation to heterogeneous environments. Inversion polymorphisms are implicated in the very high ecological flexibility of the three main malaria vector species of the Afrotropical Anopheles gambiae complex, facilitating the exploitation of anthropogenic environmental modifications and promoting a strong association with humans. In addition to extending the species' spatial and temporal distribution, inversions are associated with epidemiologically relevant mosquito behavior and physiology, underscoring their medical importance. We here present novel PCR-RFLP based assays strongly predictive of genotype for the cosmopolitan 2Rb inversion in An. coluzzii and An. gambiae, a development which overcomes the numerous constraints inherent to traditional cytological karyotyping. Methods: We designed PCR-RFLP genotyping assays based on tag SNPs previously computationally identified as strongly predictive (> 95%) of 2Rb genotype. We targeted those tags whose alternative allelic states destroyed or created the recognition site of a commercially available restriction enzyme, and designed assays with distinctive cleavage profiles for each inversion genotype. The assays were validated on 251 An. coluzzii and 451 An. gambiae cytologically karyotyped specimens from nine countries across Africa and one An. coluzzii laboratory colony. Results: For three tag SNPs, PCR-RFLP assays (denoted DraIII, MspAI, and TatI) reliably produced robust amplicons and clearly distinguishable electrophoretic profiles for all three inversion genotypes. Results obtained with the DraIII assay are ≥ 95% concordant with cytogenetic assignments in both species, while MspAI and TatI assays produce patterns highly concordant with cytogenetic assignments only in An. coluzzii or An. gambiae, respectively. Joint application of species-appropriate pairs of assays increased the concordance levels to > 99% in An. coluzzii and 98% in An. gambiae. Potential sources of discordance (e.g. imperfect association between tag and inversion, allelic dropout, additional polymorphisms in the restriction target site, incomplete or failed restriction digestion) are discussed. Conclusions: The availability of highly specific, cost effective and accessible molecular assays for genotyping 2Rb in An. gambiae and An. coluzzii allows karyotyping of both sexes and all developmental stages. These novel tools will accelerate deeper investigations into the role of this ecologically and epidemiologically important chromosomal inversion in vector biology.[Figure not available: see fulltext.

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

<p>Abstract</p> <p>Background</p> <p><it>Anopheles gambiae</it>, 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.</p> <p>Methods</p> <p>Two homokaryotypic populations of <it>A. gambiae </it>(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.</p> <p>Results</p> <p>Females carrying 2La were significantly more resistant to desiccation than 2L+^afemales at emergence and four days post-emergence, for different reasons. Teneral 2La females had lower rates of water loss than their 2L+^acounterparts, 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.</p> <p>Conclusion</p> <p>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.</p

Effective population size of Anopheles funestus chromosomal forms in Burkina Faso

BACKGROUND: As Anopheles funestus is one of the principal Afro-tropical malaria vectors, a more complete understanding of its population structure is desirable. In West and Central Africa, An. funestus population structure is complicated by the coexistence of two assortatively mating chromosomal forms. Effective population size (N(e)) is a key parameter in understanding patterns and levels of intraspecific variation, as it reflects the role of genetic drift. Here, N(e )was estimated from both chromosomal forms, Kiribina and Folonzo, in Burkina Faso. METHODS: Short-term N(e )was estimated by evaluating variation at 16 microsatellite loci across temporal samples collected annually from 2000–2002. Estimates were based on standardized variance in allele frequencies or a maximum likelihood method. Long-term N(e )was estimated from genetic diversity estimates using mtDNA sequences and microsatellites. RESULTS: For both forms, short-term and long-term N(e )estimates were on the order of 10(3 )and 10(5), respectively. Long-term N(e )estimates were larger when based on loci from chromosome 3R (both inside and outside of inversions) than loci outside of this arm. CONCLUSION: N(e )values indicate that An. funestus is not subject to seasonal bottlenecks. Though not statistically different because of large and overlapping confidence intervals, short-term N(e )estimates were consistently smaller for Kiribina than Folonzo, possibly due to exploitation of different breeding sites: permanent for Folonzo and intermittent for Kiribina. The higher long-term N(e )estimates on 3R, the arm carrying the two inversions mainly responsible for defining the chromosomal forms, give natural selection broader scope and merit further study

Phylogenetic analysis and temporal diversification of mosquitoes (Diptera: Culicidae) based on nuclear genes and morphology

<p>Abstract</p> <p>Background</p> <p>Phylogenetic analyses provide a framework for examining the evolution of morphological and molecular diversity, interpreting patterns in biogeography, and achieving a stable classification. The generic and suprageneric relationships within mosquitoes (Diptera: Culicidae) are poorly resolved, making these subjects difficult to address.</p> <p>Results</p> <p>We carried out maximum parsimony and maximum likelihood, including Bayesian, analyses on a data set consisting of six nuclear genes and 80 morphological characters to assess their ability to resolve relationships among 25 genera. We also estimated divergence times based on sequence data and fossil calibration points, using Bayesian relaxed clock methods. Strong support was recovered for the basal position and monophyly of the subfamily Anophelinae and the tribes Aedini and Sabethini of subfamily Culicinae. Divergence times for major culicid lineages date to the early Cretaceous.</p> <p>Conclusions</p> <p>Deeper relationships within the family remain poorly resolved, suggesting the need for additional taxonomic sampling. Our results support the notion of rapid radiations early in the diversification of mosquitoes.</p

Seasonal distribution of Anopheles funestus chromosomal forms from Burkina Faso

Abstract Background Previous studies of Anopheles funestus chromosomal inversion polymorphisms in Burkina Faso showed large departures from Hardy-Weinberg equilibrium and linkage disequilibrium among inversions located on different chromosomes, implying the existence of two taxonomic units ("chromosomal forms") with limited genetic flow. One chromosomal form, named Folonzo, is highly polymorphic for alternative rearrangements of 3R a, 3R b, 2R a, and 3L a; the other, Kiribina, is predominantly characterized by the standard arrangement of these inversions. To investigate the temporal distribution of these chromosomal forms, further collections were carried out in two villages near Ouagadougou where they are found in sympatry. Methods Chromosomal karyotypes were determined from indoor-resting, half-gravid females sampled within and across six breeding seasons, from December 1998 to April 2007. Results As expected, the pattern of chromosomal polymorphism in An. funestus was consistent with assortatively mating Folonzo and Kiribina forms. When samples were assigned to each chromosomal form, their relative abundance varied within successive breeding seasons in a repeating pattern of temporal variability. Relative abundance of the Folonzo form was correlated with climatic variables related to temperature and rainfall. Conclusion The relative abundance of Folonzo and Kiribina forms of An. funestus likely reflects different larval ecologies that are linked to varying climatic conditions. Further analysis of the bionomics of these vectors is recommended in light of its relevance to vector control.http://deepblue.lib.umich.edu/bitstream/2027.42/112459/1/12936_2009_Article_987.pd

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

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 nonrandomly 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 BMC Evolutionary Biology 2008, 8:309 http://www.biomedcentral.com/1471-2148/8/309 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