When genetic distance matters: measuring genetic differentiation at microsatellite loci in whole-genome scans of recent and incipient mosquito species

Complete Anopheles funestus mitogenomes reveal an ancient history of mitochondrial lineages and their distribution in southern and central Africa

Anopheles funestus s.s. is a primary vector of malaria in sub-Saharan Africa. Despite its important role in human Plasmodium transmission, evolutionary history, genetic diversity, and population structure of An. funestus in southern and central Africa remains understudied. We deep sequenced, assembled, and annotated the complete mitochondrial genome of An. funestus s.s. for the first time, providing a foundation for further genetic research of this important malaria vector species. We further analyzed the complete mitochondrial genomes of 43 An. funestus s.s. from three sites in Zambia, Democratic Republic of the Congo, and Tanzania. From these 43 mitogenomes we identified 41 unique haplotypes that comprised 567 polymorphic sites. Bayesian phylogenetic reconstruction confirmed the co-existence of two highly divergent An. funestus maternal lineages, herein defined as lineages I and II, in Zambia and Tanzania. The estimated coalescence time of these two mitochondrial lineages is ~500,000 years ago (95% HPD 426,000–594,000 years ago) with subsequent independent diversification. Haplotype network and phylogenetic analysis revealed two major clusters within lineage I, and genetic relatedness of samples with deep branching in lineage II. At this time, data suggest that the lineages are partially sympatric. This study illustrates that accurate retrieval of full mitogenomes of Anopheles vectors enables fine-resolution studies of intraspecies genetic relationships, population differentiation, and demographic history. Further investigations on whether An. funestus mitochondrial lineages represent biologically meaningful populations and their potential implications for malaria vector control are warranted.

Geographic strain differentiation of Schistosoma japonicum in the Philippines using microsatellite markers

BACKGROUND

Microsatellites have been found to be useful in determining genetic diversities of various medically-important parasites which can be used as basis for an effective disease management and control program. In Asia and Africa, the identification of different geographical strains of Schistosoma japonicum, S. haematobium and S. mansoni as determined through microsatellites could pave the way for a better understanding of the transmission epidemiology of the parasite. Thus, the present study aims to apply microsatellite markers in analyzing the populations of S. japonicum from different endemic areas in the Philippines for possible strain differentiation.

METHODOLOGY/ PRINCIPAL FINDINGS

Experimental mice were infected using the cercariae of S. japonicum collected from infected Oncomelania hupensis quadrasi snails in seven endemic municipalities. Adult worms were harvested from infected mice after 45 days of infection and their DNA analyzed against ten previously characterized microsatellite loci. High genetic diversity was observed in areas with high endemicity. The degree of genetic differentiation of the parasite population between endemic areas varies. Geographical separation was considered as one of the factors accounting for the observed difference between populations. Two subgroups have been observed in one of the study sites, suggesting that co-infection with several genotypes of the parasite might be present in the population. Clustering analysis showed no particular spatial structuring between parasite populations from different endemic areas. This result could possibly suggest varying degrees of effects of the ongoing control programs and the existing gene flow in the populations, which might be attributed to migration and active movement of infected hosts from one endemic area to another.

CONCLUSIONS/ SIGNIFICANCE

Based on the results of the study, it is reasonable to conclude that genetic diversity could be one possible criterion to assess the infection status in highly endemic areas. Genetic surveillance using microsatellites is therefore important to predict the ongoing gene flow and degree of genetic diversity, which indirectly reflects the success of the control program in schistosomiasis-endemic areas.

Population dynamics of Anopheles nuneztovari in Colombia

Anopheles nuneztovari is an important Colombian malaria vector widespread on both sides of the Andean Mountains, presenting morphological, behavioral and genetic heterogeneity throughout the country. The aim of this study was to evaluate whether the population structure and distribution of An. nuneztovari in Colombia are associated with ecological and physical barriers present in a heterogeneous landscape. Further, differences in behavior were addressed. A total of 5392 specimens of An. nuneztovari were collected. Mitochondrial and nuclear marker analyses detected subdivision among the northwest-west, northeast and east populations. For both markers, isolation by distance (~53%) and isolation by resistance (>30%) were determinants of population genetic differentiation. This suggests that physical barriers, geographical distance and ecological differences on both sides of the Andean Mountains promoted the genetic differentiation and population subdivision of An. nuneztovari in Colombia. This species showed the highest biting activity after 20:00h; indoor and outdoor preferences were found in all localities. These results indicated that the most effective interventions for controlling vector populations on both sides of the Andes need to be region-specific.

Behavior and population structure of Anopheles darlingi in Colombia

Anopheles darlingi is a widely distributed and important malaria vector in Colombia. Biogeographical and ecological heterogeneity across the Colombian distribution led to the hypothesis of behavioral and genetic differentiation among A. darlingi populations. A total of 2017 A. darlingi specimens were collected during 222 h of sampling. This vector was the most abundant anopheline species in most of the localities sampled. Subdivision between samples collected west and east of the Andes was indicated by 1) mitochondrial COI and nuclear CAD sequences from NW-W and CE-S populations (COI ΦST=0.48761-0.81974, CAD FST=0.11319-0.21321), 2) a COI haplotype network, and 3) SAMOVA. Endo- and exophagy were detected in populations west of the Andes, whereas exophagy was evident in PTG, a locality east of the Andes. Isolation by resistance was significant for COI and explained 26% of the genetic differentiation. We suggest that at a macrogeographic scale, the Andes influence the differentiation of A. darlingi in Colombia and may drive divergence, and, at a microgeographic scale, ecological differences have a significant impact on structure. These data could constitute a baseline for the design of effective vector interventions, locality-specific for the east and similar for panmictic populations west of the Andes.

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

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

A new multiplex SNP genotyping assay for detecting hybridization and introgression between the M and S molecular forms of Anopheles gambiae

The M and S forms of Anopheles gambiae have been the subject of intense study, but are morphologically indistinguishable and can only be identified using molecular techniques. PCR-based assays to distinguish the two forms have been designed and applied widely. However, the application of these assays towards identifying hybrids between the two forms, and backcrossed hybrids in particular, has been problematic as the currently available diagnostic assays are based on single locus and/or are located within a multicopy gene. Here, we present an alternative genotyping method for detecting hybridization and introgression between M and S molecular forms based on a multilocus panel of single-nucleotide polymorphisms (SNPs) fixed between the M and S forms. The panel of SNPs employed is located in so-called islands of divergence leading us to describe this method as the 'Divergence Island SNP' (DIS) assay. We show this multilocus SNP genotyping approach can robustly and accurately detect F1 hybrids as well as backcrossed individuals.

Rapid discrimination between Anopheles gambiae s.s. and Anopheles arabiensis by High-Resolution Melt (HRM) analysis

There is a need for more cost-effective options to more accurately discriminate among members of the Anopheles gambiae complex, particularly An. gambiae and Anopheles arabiensis. These species are morphologically indistinguishable in the adult stage, have overlapping distributions, but are behaviorally and ecologically different, yet both are efficient vectors of malaria in equatorial Africa. The method described here, High-Resolution Melt (HRM) analysis, takes advantage of minute differences in DNA melting characteristics, depending on the number of incongruent single nucleotide polymorphisms in an intragenic spacer region of the X-chromosome-based ribosomal DNA. The two species in question differ by an average of 13 single-nucleotide polymorphisms giving widely divergent melting curves. A real-time PCR system, Bio-Rad CFX96, was used in combination with a dsDNA-specific dye, EvaGreen, to detect and measure the melting properties of the amplicon generated from leg-extracted DNA of selected mosquitoes. Results with seven individuals from pure colonies of known species, as well as 10 field-captured individuals unambiguously identified by DNA sequencing, demonstrated that the method provided a high level of accuracy. The method was used to identify 86 field mosquitoes through the assignment of each to the two common clusters with a high degree of certainty. Each cluster was defined by individuals from pure colonies. HRM analysis is simpler to use than most other methods and provides comparable or more accurate discrimination between the two sibling species but requires a specialized melt-analysis instrument and software.

Chromosome inversions, genomic differentiation and speciation in the African malaria mosquito Anopheles gambiae

The African malaria vector, Anopheles gambiae, is characterized by multiple polymorphic chromosomal inversions and has become widely studied as a system for exploring models of speciation. Near complete reproductive isolation between different inversion types, known as chromosomal forms, has led to the suggestion that A. gambiae is in early stages of speciation, with divergence evolving in the face of considerable gene flow. We compared the standard chromosomal arrangement (Savanna form) with genomes homozygous for j, b, c, and u inversions (Bamako form) in order to identify regions of genomic divergence with respect to inversion polymorphism. We found levels of divergence between the two sub-taxa within some of these inversions (2Rj and 2Rb), but at a level lower than expected and confined near the inversion breakpoints, consistent with a gene flux model. Unexpectedly, we found that the majority of diverged regions were located on the X chromosome, which contained half of all significantly diverged regions, with much of this divergence located within exons. This is surprising given that the Bamako and Savanna chromosomal forms are both within the S molecular form that is defined by a locus near centromere of X chromosome. Two X-linked genes (a heat shock protein and P450 encoding genes) involved in reproductive isolation between the M and S molecular forms of A. gambiae were also significantly diverged between the two chromosomal forms. These results suggest that genes mediating reproductive isolation are likely located on the X chromosome, as is thought to be the case for the M and S molecular forms. We conclude that genes located on the sex chromosome may be the major force driving speciation between these chromosomal forms of A. gambiae.

Evidence for X-linked introgression between molecular forms of Anopheles gambiae from Angola

The M and S molecular forms of the African malaria vector Anopheles gambiae (Diptera: Culicidae) are morphologically identical incipient species in which reproductive isolation is incomplete, enabling low-level gene flow between forms. In an attempt to find differences between the M and S forms, sequence variation was studied at loci along the X chromosome in adult female An. gambiae from Angola. A high proportion of M form specimens from Angola (79% of the 456 X chromosomes sampled) were found to contain a 16-bp insertion in intron 4 of the X-linked GPRCCK1 locus, relative to the AgamP3 release of the An. gambiae PEST genome sequence. The insertion was in Hardy-Weinberg equilibrium in Angolan M form populations. The same insertion was found in all S form specimens examined, regardless of where in Africa they were sampled, but was absent from a sample of M form specimens collected in Ghana, Bioko and Mali. In M form specimens from Angola, there was an association between alleles at the GPRCCK1 locus and those at a microsatellite locus, AGXH678, close to the centromere of the X chromosome, with significant linkage disequilibrium between loci separated by 0.472 Mbp (P < 0.033). We show that the insertion results from introgression from the S form into the M form, rather than from the retention of an ancestral character. Gene flow from the S to M form could allow genes of adaptive value to be transferred, including those conferring insecticide resistance and others influencing ecology and behaviour, and thus malaria transmission and control. We discuss factors that may have led to this introgression event.

No evidence for biased co-transmission of speciation islands in Anopheles gambiae

Genome-scale scans have revealed highly heterogeneous levels of divergence between closely related taxa in many systems. Generally, a small number of regions show high differentiation, with the rest of the genome showing no or only low levels of divergence. These patterns have been interpreted as evidence for ongoing speciation-with-gene-flow, with introgression homogenizing the whole genome except loci involved in reproductive isolation. However, as the number of selected loci increases, the probability of introgression at unselected loci decreases unless there is a transmission ratio distortion causing an over-representation of specific combinations of alleles. Here we examine the transmission of three 'speciation islands' that contain fixed differences between the M and S forms of the mosquito, Anopheles gambiae. We made reciprocal crosses between M and S parents and genotyped over 2000 F(2) individuals, developing a hierarchical likelihood model to identify specific genotypes that are under- or over-represented among the recombinant offspring. Though our overall results did not match the expected number of F(2) genotypes, we found no biased co-transmission among M or S alleles in the three islands. Our likelihood model did identify transmission ratio distortion at two of the three islands, but this distortion was small (approx. 3%) and in opposite directions for the two islands. We discuss how our results impinge on hypotheses of current gene flow between M and S and ongoing speciation-with-gene-flow in this system.

Comparison of wing geometry data and genetic data for assessing the population structure of Aedes aegypti

Aedes aegypti is the most important vector of dengue viruses in tropical and subtropical regions. Because vaccines are still under development, dengue prevention depends primarily on vector control. Population genetics is a common approach in research involving Ae. aegypti. In the context of medical entomology, wing morphometric analysis has been proposed as a strong and low-cost complementary tool for investigating population structure. Therefore, we comparatively evaluated the genetic and phenotypic variability of population samples of Ae. aegypti from four sampling sites in the metropolitan area of São Paulo city, Brazil. The distances between the sites ranged from 7.1 to 50 km. This area, where knowledge on the population genetics of this mosquito is incipient, was chosen due to the thousands of dengue cases registered yearly. The analysed loci were polymorphic, and they revealed population structure (global F(ST)=0.062; p<0.05) and low levels of gene flow (Nm=0.47) between the four locations. Principal component and discriminant analyses of wing shape variables (18 landmarks) demonstrated that wing polymorphisms were only slightly more common between populations than within populations. Whereas microsatellites allowed for geographic differentiation, wing geometry failed to distinguish the samples. These data suggest that microevolution in this species may affect genetic and morphological characters to different degrees. In this case, wing shape was not validated as a marker for assessing population structure. According to the interpretation of a previous report, the wing shape of Ae. aegypti does not vary significantly because it is stabilised by selective pressure.

Discovery of common urinary biomarkers for hepatotoxicity induced by carbon tetrachloride, acetaminophen and methotrexate by mass spectrometry-based metabolomics

Liver toxicity represents an important healthcare issue because it causes significant morbidity and mortality and can be difficult to predict before symptoms appear owing to drug therapy or exposure to toxicants. Using metabolomic techniques, we discovered common biomarkers for the prediction of hepatotoxicity in rat urine using mass spectrometry. For this purpose, liver toxicity was induced by 5 days of oral administration of carbon tetrachloride (1 ml kg(-1) per day), acetaminophen (1000 mg kg(-1) per day) and methotrexate (50 mg kg(-1) per day). Serum levels of alkaline phosphatase aspartate aminotransferase, alanine aminotransferase and histopathology in liver tissue were then checked to demonstrate liver toxicity. Global metabolic profiling with UPLC-TOF-MS (ultraperformance liquid chromatography-mass spectrometry), multivariate analysis (partial least square-discriminant analysis, hierarchical analysis) and database searching were performed to discover common biomarkers for liver toxicity induced by these three compounds. Urinary concentrations of the newly discovered biomarkers were then quantified to confirm them as biomarkers of hepatotoxicity with targeted metabolic profiling using GC (gas chromatography)-MS and CE (capillary electrophoresis)-MS. In the results, steroids, amino acids and bile acids were metabolically changed between the control and drug-treated groups in global metabolic profiling; 11β-hydroxyandrosterone, epiandrosterone, estrone, 11-dehydrocorticosterone, glycine, alanine, valine, leucine, dl-ornithine, 3-methylhistidine, cholic acid and lithocholic acid were selected as liver toxicity biomarkers after performing targeted metabolic profiling. In conclusion, we discovered metabolite biomarkers belonging to three different metabolic pathways to check for liver toxicity with mass spectrometry from a metabolomics study that could be used to evaluate hepatotoxicity induced by drugs or other toxic compounds.

Malaria vectors in the Republic of Benin: distribution of species and molecular forms of the Anopheles gambiae complex

Members of the Anopheles gambiae complex are among the best malaria vectors in the world, but their vectorial capacities vary between species and populations. A large-scale sampling of An. gambiae sensu lato was carried out in 2006 and 2007 in various bioclimatic areas of Benin (West Africa). The objective of this study was to collate data on the relative frequencies of species and forms within the An. gambiae complex and to produce a map of their spatial distribution. Sampling took place at 30 sites and 2122 females were analyzed. Two species were identified through molecular methods. The overall collection showed a preponderance of An. gambiae s.s., but unexpectedly, An. arabiensis was reported in the coastal-Guinean bioclimatic area characterized by a mean annual rainfall of >1500 mm where only An. gambiae s.s. was reported previously. Our study of Benin indicates that An. arabiensis would be adapted not only to the urban areas but also to the rural humid regions. Among 1717 An. gambiae s.s., 26.5% were of the M form and 73.3% were S form. Few hybrid specimens between the M and S forms were observed (0.2%). Only the spatial distribution of the M form appears to be mainly a function of bioclimatic area. Factors that influence the distribution of these malaria vectors are discussed. This study underlines the need of further investigations of biological, ecological, and behavioral traits of these species and forms to better appreciate their vectorial capacities. Acquisition of entomological field data appears essential to better estimate the stratification of malaria risk and help improve malaria vector control interventions.

High levels of hybridization between molecular forms of Anopheles gambiae from Guinea Bissau

In the malaria vector Anopheles gambiae Giles sensu stricto, two molecular forms denoted M and S are considered units of incipient speciation within this species. Very low hybrid frequencies and significant genetic differentiation have been found in sympatric M- and S-form populations. We studied the molecular form composition and the degree of genetic differentiation at 15 microsatellites in two samples of An. gambiae collected in two consecutive years from Bissau, Guinea Bissau. High frequencies of M/S hybrids (19-24%) were found in this area. Coincidently, very low levels of genetic differentiation were detected between forms when analysis involved microsatellites mapped at chromosome-3 (mean Fst, 0.000-0.002). The single exception was the X-linked AGXH678, for which high differentiation was measured (Fst, 0.158-0.301). This locus maps near the centromere of chromosome X, a low recombination region in which selection is likely to promote divergence between M and S forms. These results strongly suggest that the degree of isolation between M and S forms, considered the units of incipient speciation within An. gambiae, is not homogenous throughout the species distribution range.

Insertion polymorphisms of SINE200 retrotransposons within speciation islands of Anopheles gambiae molecular forms

Background

SINEs (Short INterspersed Elements) are homoplasy-free and co-dominant genetic markers which are considered to represent useful tools for population genetic studies, and could help clarifying the speciation processes ongoing within the major malaria vector in Africa, Anopheles gambiae s.s. Here, we report the results of the analysis of the insertion polymorphism of a nearly 200 bp-long SINE (SINE200) within genome areas of high differentiation (i.e. "speciation islands") of M and S A. gambiae molecular forms.

Methods

A SINE-PCR approach was carried out on thirteen SINE200 insertions in M and S females collected along the whole range of distribution of A. gambiae s.s. in sub-Saharan Africa. Ten specimens each for Anopheles arabiensis, Anopheles melas, Anopheles quadriannulatus A and 15 M/S hybrids from laboratory crosses were also analysed.

Results

Eight loci were successfully amplified and were found to be specific for A. gambiae s.s.: 5 on 2L chromosome and one on X chromosome resulted monomorphic, while two loci positioned respectively on 2R (i.e. S200 2R12D) and X (i.e. S200 X6.1) chromosomes were found to be polymorphic. S200 2R12D was homozygote for the insertion in most S-form samples, while intermediate levels of polymorphism were shown in M-form, resulting in an overall high degree of genetic differentiation between molecular forms (Fst = 0.46 p < 0.001) and within M-form (Fst = 0.46 p < 0.001). The insertion of S200 X6.1 was found to be fixed in all M- and absent in all S-specimens. This led to develop a novel easy-to-use PCR approach to straightforwardly identify A. gambiae molecular forms. This novel approach allows to overcome the constraints associated with markers on the rDNA region commonly used for M and S identification. In fact, it is based on a single copy and irreversible SINE200 insertion and, thus, is not subjected to peculiar evolutionary patterns affecting rDNA markers, e.g. incomplete homogenization of the arrays through concerted evolution and/or mixtures of M and S IGS-sequences among the arrays of single chromatids.

Conclusion

The approach utilized allowed to develop new easy-to-use co-dominant markers for the analysis of genetic differentiation between M and S-forms and opens new perspectives in the study of the speciation process ongoing within A. gambiae.

The molecular forms of Anopheles gambiae: a phenotypic perspective

The African malaria mosquito Anopheles gambiae is undergoing speciation, being split into the M and S molecular forms. Speciation is the main process promoting biological diversity, thus, new vector species might complicate disease transmission. Genetic differentiation between the molecular forms has been extensively studied, but phenotypic differences between them, the evolutionary forces that generated divergence, and the mechanisms that maintain their genetic isolation have only recently been addressed. Here, we review recent studies suggesting that selection mediated by larval predation and competition promoted divergence between temporary and permanent freshwater habitats. These differences explain the sharp discontinuity in distribution of the molecular forms between rice fields and surrounding savanna, but they can also explain the concurrent cline between humid and arid environments due to the dependence on permanent habitats in the latter. Although less pronounced, differences in adult body size, reproductive output, and longevity also suggest that the molecular forms have adapted to distinct niches. Reproductive isolation between the molecular forms is achieved by spatial swarm segregation, although within-swarm mate recognition appears to play a role in certain locations. The implications of these results to disease transmission and control are discussed and many of the gaps in our understanding are highlighted.

High genetic differentiation between the M and S molecular forms of Anopheles gambiae in Africa

Background

Anopheles gambiae, a major vector of malaria, is widely distributed throughout sub-Saharan Africa. In an attempt to eliminate infective mosquitoes, researchers are trying to develop transgenic strains that are refractory to the Plasmodium parasite. Before any release of transgenic mosquitoes can be envisaged, we need an accurate picture of the differentiation between the two molecular forms of An. gambiae, termed M and S, which are of uncertain taxonomic status.

Methodology/Principal Findings

Insertion patterns of three transposable elements (TEs) were determined in populations from Benin, Burkina Faso, Cameroon, Ghana, Ivory Coast, Madagascar, Mali, Mozambique, Niger, and Tanzania, using Transposon Display, a TE-anchored strategy based on Amplified Fragment Length Polymorphism. The results reveal a clear differentiation between the M and S forms, whatever their geographical origin, suggesting an incipient speciation process.

Conclusions/Significance

Any attempt to control the transmission of malaria by An. gambiae using either conventional or novel technologies must take the M/S genetic differentiation into account. In addition, we localized three TE insertion sites that were present either in every individual or at a high frequency in the M molecular form. These sites were found to be located outside the chromosomal regions that are suspected of involvement in the speciation event between the two forms. This suggests that these chromosomal regions are either larger than previously thought, or there are additional differentiated genomic regions interspersed with undifferentiated regions.

Microsatellite data suggest significant population structure and differentiation within the malaria vector Anopheles darlingi in Central and South America

Background

Anopheles darlingi is the most important malaria vector in the Neotropics. An understanding of A. darlingi's population structure and contemporary gene flow patterns is necessary if vector populations are to be successfully controlled. We assessed population genetic structure and levels of differentiation based on 1,376 samples from 31 localities throughout the Peruvian and Brazilian Amazon and Central America using 5–8 microsatellite loci.

Results

We found high levels of polymorphism for all of the Amazonian populations (mean R_S = 7.62, mean H_O = 0.742), and low levels for the Belize and Guatemalan populations (mean R_S = 4.3, mean H_O = 0.457). The Bayesian clustering analysis revealed five population clusters: northeastern Amazonian Brazil, southeastern and central Amazonian Brazil, western and central Amazonian Brazil, Peruvian Amazon, and the Central American populations. Within Central America there was low non-significant differentiation, except for between the populations separated by the Maya Mountains. Within Amazonia there was a moderate level of significant differentiation attributed to isolation by distance. Within Peru there was no significant population structure and low differentiation, and some evidence of a population expansion. The pairwise estimates of genetic differentiation between Central America and Amazonian populations were all very high and highly significant (F_ST = 0.1859 – 0.3901, P < 0.05). Both the D_A and F_ST distance-based trees illustrated the main division to be between Central America and Amazonia.

Conclusion

We detected a large amount of population structure in Amazonia, with three population clusters within Brazil and one including the Peru populations. The considerable differences in N_e among the populations may have contributed to the observed genetic differentiation. All of the data suggest that the primary division within A. darlingi corresponds to two white gene genotypes between Amazonia (genotype 1) and Central America, parts of Colombia and Venezuela (genotype 2), and are in agreement with previously published mitochondrial COI gene sequences interpreted as incipient species. Overall, it appears that two main factors have contributed to the genetic differentiation between the population clusters: physical distance between the populations and the differences in effective population sizes among the subpopulations.

Mosaic genome architecture of the Anopheles gambiae species complex

Background

Attempts over the last three decades to reconstruct the phylogenetic history of the Anopheles gambiae species complex have been important for developing better strategies to control malaria transmission.

Methodology

We used fingerprint genotyping data from 414 field-collected female mosquitoes at 42 microsatellite loci to infer the evolutionary relationships of four species in the A. gambiae complex, the two major malaria vectors A. gambiae sensu stricto (A. gambiae s.s.) and A. arabiensis, as well as two minor vectors, A. merus and A. melas.

Principal Findings

We identify six taxonomic units, including a clear separation of West and East Africa A. gambiae s.s. S molecular forms. We show that the phylogenetic relationships vary widely between different genomic regions, thus demonstrating the mosaic nature of the genome of these species. The two major malaria vectors are closely related and closer to A. merus than to A. melas at the genome-wide level, which is also true if only autosomes are considered. However, within the Xag inversion region of the X chromosome, the M and two S molecular forms are most similar to A. merus. Near the X centromere, outside the Xag region, the two S forms are highly dissimilar to the other taxa. Furthermore, our data suggest that the centromeric region of chromosome 3 is a strong discriminator between the major and minor malaria vectors.

Conclusions

Although further studies are needed to elucidate the basis of the phylogenetic variation among the different regions of the genome, the preponderance of sympatric admixtures among taxa strongly favor introgression of different genomic regions between species, rather than lineage sorting of ancestral polymorphism, as a possible mechanism.

Genetic population structure of Anopheles gambiae in Equatorial Guinea

BACKGROUND

Patterns of genetic structure among mosquito vector populations in islands have received particular attention as these are considered potentially suitable sites for experimental trials on transgenic-based malaria control strategies. In this study, levels of genetic differentiation have been estimated between populations of Anopheles gambiae s.s. from the islands of Bioko and Annobón, and from continental Equatorial Guinea (EG) and Gabon.

METHODS

Genotyping of 11 microsatellite loci located in chromosome 3 was performed in three island samples (two in Bioko and one in Annobón) and three mainland samples (two in EG and one in Gabon). Four samples belonged to the M molecular form and two to the S-form. Microsatellite data was used to estimate genetic diversity parameters, perform demographic equilibrium tests and analyse population differentiation.

RESULTS

High levels of genetic differentiation were found between the more geographically remote island of Annobón and the continent, contrasting with the shallow differentiation between Bioko island, closest to mainland, and continental localities. In Bioko, differentiation between M and S forms was higher than that observed between island and mainland samples of the same molecular form.

CONCLUSION

The observed patterns of population structure seem to be governed by the presence of both physical (the ocean) and biological (the M-S form discontinuity) barriers to gene flow. The significant degree of genetic isolation between M and S forms detected by microsatellite loci located outside the "genomic islands" of speciation identified in A. gambiae s.s. further supports the hypothesis of on-going incipient speciation within this species. The implications of these findings regarding vector control strategies are discussed.

Locus- and Population-Specific Selection and Differentiation between Incipient Species of Anopheles gambiae

Anopheles gambiae, the primary mosquito vector of malaria in sub-Saharan Africa, is divided into 2 sympatric incipient species known as M form and S form. Recent genomic analysis of each form revealed that differentiation between forms is clustered into 3 unlinked regions of the genome. Here, we expand the investigation of these "genomic islands of speciation" to multiple populations, including all of the genes across one of the islands. Differentiation between the M and S forms in 2 of the islands is complete across all individuals in all populations, confirming that the M and S forms are reproductively isolated taxa. Differentiation at the third island (on chromosome 2R) is limited to Cameroon populations. There is reduced variation in the M form in Cameroon at this location and increased divergence to the outgroup Anopheles arabiensis, supporting an association of adaptation with reproductive isolation.

Insertion polymorphism of transposable elements and population structure of Anopheles gambiae M and S molecular forms in Cameroon

The insertion polymorphism of five transposable element (TE) families was studied by Southern blots in several populations of the M and S molecular forms of the mosquito Anopheles gambiae sensu stricto from southern Cameroon. We showed that the mean TE insertion site number and the within-population insertion site polymorphism globally differed between the M and S molecular forms. The comparison of the TE insertion profiles of the populations revealed a significant differentiation between these two molecular forms (0.163 < Phi(ST) < 0.371). We cloned several insertions of a non-LTR retrotransposon (Aara8) that were fixed in one form and absent in the other one. The only insertion that could be clearly located on a chromosome arm mapped to cytological division 6 of chromosome X, confirming the importance of this region in the ongoing speciation between the M and S molecular forms.

Higher genetic variation estimated by microsatellites compared to isoenzyme markers in Aedes aegypti from Rio de Janeiro

Aedes aegypti populations from five districts in Rio de Janeiro were analyzed using five microsatellites and six isoenzyme markers, to assess the amount of variation and patterns of gene flow at local levels. Microsatellite loci were polymorphic enough to detect genetic differentiation of populations collected at small geographic scales (e.g. within a city). Ae. aegypti populations were highly differentiated as well in the city center as in the outskirt. Thus, dengue virus propagation by mosquitoes could be as efficient in the urban area as in the outskirt of Rio de Janeiro, the main entry point of dengue in Brazil.

Reduced recombination rate and genetic differentiation between the M and S forms of Anopheles gambiae s.s

Genetic differentiation between the largely sympatric molecular forms M and S of Anopheles gambiae appears mostly limited to division 6 and part of division 5 of the X chromosome. This region is adjacent to the centromere and includes the rDNA that was used to define these forms. This localized differentiation between populations that experience gene flow strongly suggests that this region contains genes responsible for reproductive isolation. Regions adjacent to centromeres are known to experience less recombination in several species and it has recently been suggested that low recombination rates can facilitate the accumulation and maintenance of isolation genes in partially isolated populations. Therefore, we measured the recombination rate in division 5D/6 directly and estimate that it is at least 16-fold reduced across this region compared to the remainder of the X chromosome. Additionally, sequence data from four loci from field-collected mosquitoes from several West African countries show very strong differentiation between the molecular forms in division 5D/6, whereas none was observed in two loci elsewhere on the X chromosome. Furthermore, genetic variation was substantially lower in division 5D/6 compared to the two reference loci, and the inferred genealogies of the division 5D/6 genes show patterns consistent with selective sweeps. This suggests that the reduced recombination rate has increased the effect of selection on this region and that our data are consistent with the hypothesis that reduced recombination rates can play a role in the accumulation of isolation genes in the face of gene flow.

Ecological zones rather than molecular forms predict genetic differentiation in the malaria vector Anopheles gambiae s.s. in Ghana

The malaria mosquito Anopheles gambiae s.s. is rapidly becoming a model for studies on the evolution of reproductive isolation. Debate has centered on the taxonomic status of two forms (denoted M and S) within the nominal taxon identified by point mutations in the X-linked rDNA region. Evidence is accumulating that there are significant barriers to gene flow between these forms, but that the barriers are not complete throughout the entire range of their distribution. We sampled populations from across Ghana and southern Burkina Faso, West Africa, from areas where the molecular forms occurred in both sympatry and allopatry. Neither Bayesian clustering methods nor F(ST)-based analysis of microsatellite data found differentiation between the M and S molecular forms, but revealed strong differentiation among different ecological zones, irrespective of M/S status and with no detectable effect of geographical distance. Although no M/S hybrids were found in the samples, admixture analysis detected evidence of contemporary interform gene flow, arguably most pronounced in southern Ghana where forms occur sympatrically. Thus, in the sampled area of West Africa, lack of differentiation between M and S forms likely reflects substantial introgression, and ecological barriers appear to be of greater importance in restricting gene flow.

Species and populations of the Anopheles gambiae complex in Cameroon with special emphasis on chromosomal and molecular forms of Anopheles gambiae s.s

We studied the geographical distribution of species, chromosomal, and molecular forms of the Anopheles gambiae Giles (Diptera: Culicidae) complex in 23 sites in Cameroon, Central Africa. Almost all the specimens collected in the four northern-most arid sites were Anopheles arabiensis. Anopheles melas was found in a rural locality surrounded by mangrove swamps, on the Atlantic Coast. In total, 1,525 An. gambiae s.s. females were identified down to their molecular form, and inversion polymorphisms on polytene chromosomes were scored from 186 half-gravid females. The Forest chromosomal form, with standard arrangements almost fixed on both arms of chromosome-2, was the only one observed in the southern, more humid localities. Karyotypes typical of Savanna and Mopti were recorded northwards, in the humid savannas of the Adamawa Province. The molecular forms M and S were widespread throughout Cameroon, and assort independently from the chromosomal forms. S-form populations were characterized by karyotypes typical of Forest and Savanna chromosomal forms, and M-form populations were characterized by karyotypes typical of Forest, Savanna, and Mopti. No M/S hybrid patterns were detected, although M and S mosquitoes were sympatric in 15 sites, providing further evidence for positive assortative mating within molecular forms. The observed ecogeographical distribution of M and S was peculiar: the ecological parameters involved in this distribution still need to be clarified as well as the possible role of competitive exclusion between chromosomally homosequential molecular forms. No difference was observed in host preference or in Plasmodium falciparum infection rates between sympatric M and S populations.

Genomic islands of speciation in Anopheles gambiae

The African malaria mosquito, Anopheles gambiae sensu stricto (A. gambiae), provides a unique opportunity to study the evolution of reproductive isolation because it is divided into two sympatric, partially isolated subtaxa known as M form and S form. With the annotated genome of this species now available, high-throughput techniques can be applied to locate and characterize the genomic regions contributing to reproductive isolation. In order to quantify patterns of differentiation within A. gambiae, we hybridized population samples of genomic DNA from each form to Affymetrix GeneChip microarrays. We found that three regions, together encompassing less than 2.8 Mb, are the only locations where the M and S forms are significantly differentiated. Two of these regions are adjacent to centromeres, on Chromosomes 2L and X, and contain 50 and 12 predicted genes, respectively. Sequenced loci in these regions contain fixed differences between forms and no shared polymorphisms, while no fixed differences were found at nearby control loci. The third region, on Chromosome 2R, contains only five predicted genes; fixed differences in this region were also verified by direct sequencing. These “speciation islands” remain differentiated despite considerable gene flow, and are therefore expected to contain the genes responsible for reproductive isolation. Much effort has recently been applied to locating the genes and genetic changes responsible for reproductive isolation between species. Though much can be inferred about speciation by studying taxa that have diverged for millions of years, studying differentiation between taxa that are in the early stages of isolation will lead to a clearer view of the number and size of regions involved in the genetics of speciation. Despite appreciable levels of gene flow between the M and S forms of A. gambiae, we were able to isolate three small regions of differentiation where genes responsible for ecological and behavioral isolation are likely to be located. We expect reproductive isolation to be due to changes at a small number of loci, as these regions together contain only 67 predicted genes. Concentrating future mapping experiments on these regions should reveal the genes responsible for reproductive isolation between forms.

Using DNA microarrays, the authors identify 3 small regions of the genome that differ between two forms of hybridizing mosquitoes; regions that are likely to contain the genes responsible for reproductive isolation.

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.

On the distribution and genetic differentiation of Anopheles gambiae s.s. molecular forms

This paper summarises published and unpublished data on the spatial and temporal distribution, and on the genetic characterisation of molecular forms M and S of Anopheles gambiae s.s. The two forms are characterised by a high level of gene-flow restriction, by a largely overlapping geographical and temporal distribution, and by a low degree of genetic differentiation. Floating paracentric inversions on chromosome-2 are shown to be shared by the two forms, although with very different frequencies of alternative arrangements, confirming that these inversions are most probably involved in ecotypic adaptation, rather than in the building of reproductive barriers. Further studies and tools are needed to throw light on the genetic and biological differentiation of M and S to improve the knowledge of the real composition of the vector system, of its demography, population genetics and dynamics, also in view of the possible consequences on the transmission of human pathogens in sub-Saharan Africa. Preliminary results and perspectives of the use of transposable element insertion sites as markers of genetic differentiation and tools for population genetic studies are discussed.

Effect of seminal fluids in mating between M and S forms of Anopheles gambiae

Previous studies have shown that sympatric populations of M and S molecular forms of Anopheles gambiae sensu stricto exhibit strong assortative mating. In the few documented cases of cross-mating between M and S forms, females that mated with amale of the alternative form were often also mated with a male of their own form. A potential explanation for the association between cross-mating and double mating could be that male accessory gland or sperm proteins that are responsible for inducing refractoriness to further mating by females have diverged between the M and S forms. This mechanism of postmating reproductive isolation would have important implications for our understanding of the speciation processes in the An. gambiae complex. We tested for this mechanism, by comparing the likelihood of mating, feeding, and laying eggs, as well as the fertility of females presented with males of their own form or the alternate form in the laboratory. We also compared the likelihood of remating in cross-mated and assortatively-mated females, and we analyzed their progeny to unravel patterns of sperm precedence. We found that cross-mated females differed from assortatively-mated females only in terms of egg-hatching rate and larval survival but that these effects could be attributed to hybrid vigor rather than differential response to seminal products. Cross-mating between forms was not associated with remating behavior. These results indicate that the sex proteins responsible for inhibiting further insemination and triggering the gonotrophic cycle in females have not diverged between these M and S populations. We discuss alternative explanations for the patterns of cross-mating and multiple mating observed in the field.

Gene flow between chromosomal forms of the malaria vector Anopheles funestus in Cameroon, Central Africa, and its relevance in malaria fighting

Knowledge of population structure in a major vector species is fundamental to an understanding of malaria epidemiology and becomes crucial in the context of genetic control strategies that are being developed. Despite its epidemiological importance, the major African malaria vector Anopheles funestus has received far less attention than members of the Anopheles gambiae complex. Previous chromosomal data have shown a high degree of structuring within populations from West Africa and have led to the characterization of two chromosomal forms, "Kiribina" and "Folonzo." In Central Africa, few data were available. We thus undertook assessment of genetic structure of An. funestus populations from Cameroon using chromosomal inversions and microsatellite markers. Microsatellite markers revealed no particular departure from panmixia within each local population and a genetic structure consistent with isolation by distance. However, cytogenetic studies demonstrated high levels of chromosomal heterogeneity, both within and between populations. Distribution of chromosomal inversions was not random and a cline of frequency was observed, according to ecotypic conditions. Strong deficiency of heterokaryotypes was found in certain localities in the transition area, indicating a subdivision of An. funestus in chromosomal forms. An. funestus microsatellite genetic markers located within the breakpoints of inversions are not differentiated in populations, whereas in An. gambiae inversions can affect gene flow at marker loci. These results are relevant to strategies for control of malaria by introduction of transgenes into populations of vectors.

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.

Multilevel analyses of genetic differentiation in Anopheles gambiae s.s. reveal patterns of gene flow important for malaria-fighting mosquito projects

Malaria control projects based on the introduction and spread of transgenes into mosquito populations depend on the extent of isolation between those populations. On the basis of the distribution of paracentric inversions, Anopheles gambiae has been subdivided into five subspecific chromosomal forms. Estimating gene flow between and within these forms of An. gambiae presents a number of challenges. We compared patterns of genetic divergence (F(ST)) between sympatric populations of the Bamako and Mopti forms at five sites. We used microsatellite loci within the j inversion on chromosome 2, which is fixed in the Bamako form but absent in the Mopti form, and microsatellites on chromosome 3, a region void of inversions. Estimates of genetic diversity and F(ST)'s suggest genetic exchanges between forms for the third chromosome but little for the j inversion. These results suggest a role for the inversion in speciation. Extensive gene flow within forms among sites resulted in populations clustering according to form despite substantial gene flow between forms. These patterns underscore the low levels of current gene flow between chromosomal forms in this area of sympatry. Introducing refractoriness genes in areas of the genome void of inversions may facilitate their spread within forms but their passage between forms may prove more difficult than previously thought.

Natural hybridization and the evolution of domesticated, pest and disease organisms

The role of natural hybridization in the evolutionary history of numerous species is well recognized. The impact of introgressive hybridization and hybrid speciation has been documented especially in plant and animal assemblages. However, there remain certain areas of investigation for which natural hybridization and its consequences remain under-studied and under-appreciated. One such area involves the evolution of organisms that positively or negatively affect human populations. In this review, I highlight exemplars of how natural hybridization has contributed to the evolution of (i) domesticated plants and animals; (ii) pests; (iii) human disease vectors; and (iv) human pathogens. I focus on the effects from genetic exchange that may lead to the acquisition of novel phenotypes and thus increase the beneficial or detrimental (to human populations) aspects of the various taxa.

Quantitative trait loci in Anopheles gambiae controlling the encapsulation response against Plasmodium cynomolgi Ceylon

Background

Anopheles gambiae females are the world's most successful vectors of human malaria. However, a fraction of these mosquitoes is refractory to Plasmodium development. L3-5, a laboratory selected refractory strain, encapsulates transforming ookinetes/early oocysts of a wide variety of Plasmodium species. Previous studies on these mosquitoes showed that one major (Pen1) and two minor (Pen2, Pen3) autosomal dominant quantitative trait loci (QTLs) control the melanotic encapsulation response against P. cynomolgi B, a simian malaria originating in Malaysia.

Results

We have investigated the response of L3-5 to infection with P. cynomolgi Ceylon, a different but related parasite species, in crosses with the susceptible strain 4Arr. Refractoriness to this parasite is incompletely recessive. Infection and genotyping of F2 intercross females at genome-spanning microsatellite loci revealed that 3 autosomal QTLs control encapsulation of this species. Two loci map to the regions containing Pen2 and Pen3. The novel QTL maps to chromosome 3R, probably to polytene division 32 or 33. Thus the relative contribution of any QTL to oocyst encapsulation varies with the species of parasite. Further, different QTLs were most readily identified in different F2 families. This, like the F1 data, suggests that L3-5 is not genetically homogeneous and that somewhat different pathways may be used to achieve an encapsulation response.

Conclusion

We have shown here that different QTLs are involved in responses against different Plasmodium parasites.

MOLECULAR AND QUANTITATIVE TRAIT VARIATION WITHIN AND AMONG POPULATIONS OF THE INTERTIDAL COPEPOD TIGRIOPUS CALIFORNICUS

While molecular and quantitative trait variation may be theoretically correlated, empirical studies using both approaches frequently reveal discordant patterns, and these discrepancies can contribute to our understanding of evolutionary processes. Here, we assessed genetic variation in six populations of the copepod Tigriopus californicus. Molecular variation was estimated using five polymorphic microsatellite loci, and quantitative variation was measured using 22-life history and morphometric characters. Within populations, no correlation was found between the levels of molecular variation (heterozygosity) and quantitative variation (heritability). Between populations, quantitative subdivision (Q(ST)) was correlated with molecular subdivision when measured as F(ST) but not when measured as R(ST). Unlike most taxa studied to date, the overall level of molecular subdivision exceeded the level of quantitative subdivision (F(ST) = 0.80, R(ST) = 0.89, Q(ST) = 0.30). Factors that could contribute to this pattern include stabilizing or fluctuating selection on quantitative traits or accelerated rates of molecular evolution.

Genetic delineation of sibling species of the pest fruit fly Bactocera (Diptera: Tephritidae) using microsatellites

Using a large set of microsatellites, the genetic relationships between three closely related Australian fruit fly species, Bactrocera tryoni (Froggatt), B. neohumeralis (Hardy) and B. aquilonis(May) were investigated. Bactrocera tryoni and B. neohumeralis are sympatric, while B. aquilonisis allopatric to both. The sympatric species, B. tryoni and B. neohumeralis, were found to be genetically distinct. It is likely that despite differences in mating time between these two species, some gene flow still occurs. In contrast, the sibling species B. tryoni and B. aquilonis were found to be closely related, despite allopatry. The level of genetic divergence was similar to that found within eastern Australian populations of B. tryoni. Consideration of all available genetic data suggests that this similarity is not due to recent (i.e. within the last 30 years) displacement of B. aquilonis by B. tryoni from the B. aquilonis region (north-western Australia). Instead the data suggests that, at least in the areas sampled, asymmetrical hybridization may have occurred over a longer timescale.

Locus-specific genetic differentiation at Rw among warfarin-resistant rat (Rattus norvegicus) populations

Populations may diverge at fitness-related genes as a result of adaptation to local conditions. The ability to detect this divergence by marker-based genomic scans depends on the relative magnitudes of selection, recombination, and migration. We survey rat (Rattus norvegicus) populations to assess the effect that local selection with anticoagulant rodenticides has had on microsatellite marker variation and differentiation at the warfarin resistance gene (Rw) relative to the effect on the genomic background. Initially, using a small sample of 16 rats, we demonstrate tight linkage of microsatellite D1Rat219 to Rw by association mapping of genotypes expressing an anticoagulant-rodenticide-insensitive vitamin K 2,3-epoxide reductase (VKOR). Then, using allele frequencies at D1Rat219, we show that predicted and observed resistance levels in 27 populations correspond, suggesting intense and recent selection for resistance. A contrast of F(ST) values between D1Rat219 and the genomic background revealed that rodenticide selection has overwhelmed drift-mediated population structure only at Rw. A case-controlled design distinguished these locus-specific effects of selection at Rw from background levels of differentiation more effectively than a population-controlled approach. Our results support the notion that an analysis of locus-specific population genetic structure may assist the discovery and mapping of novel candidate loci that are the object of selection or may provide supporting evidence for previously identified loci.

Genetic markers for study of the anopheline vectors of human malaria

Human malaria is truly a disease of global proportions and is one of the most broadly distributed vector-borne infections. Anopheline mosquitoes are the exclusive vectors of human malaria. A handful of species predominate as the most notorious malaria vectors, but the species and forms involved in the transmission of human malaria world-wide are incredibly diverse. Many of the anophelines that vector malaria exist as members of species complexes that often contain vector and non-vector species. Additionally, single anopheline species often exhibit significant heterogeneity across the species' range. This phenotypic and genotypic plasticity exacerbates the difficulties in identification of vector populations and implementation of effective surveillance and control strategies. Polytene chromosome investigations were among the first to provide researchers with tangible genetic markers that could be used to differentiate between what are now recognised as species and chromosomal forms of anopheline mosquitoes. The advent of the polymerase chain reaction gave access to the molecular genetics of genomes and the techniques that followed have facilitated investigation of the genetics of individual specimens or population size samples. The variety and number of genetic markers available for the study of malaria vectors has literally exploded in the last 10 years. Markers have expanded from the 'traditional tools' to include a vast array of molecular markers. Contemporary markers range from what are now referred to as 'classical genetic markers' to methods used to detect and identify single nucleotide polymorphisms and finally to highly polymorphic markers. One of the greatest advantages of this wide variety of genetic markers is that researchers may choose to utilise any combination of markers or techniques to address multifaceted questions relating to malaria transmission. These molecular markers have proven useful in a wide variety of applications including molecular taxonomy, evolutionary systematics, population genetics, genetic mapping, and investigation of defined phenotypes.

The genome sequence of the malaria mosquito Anopheles gambiae

Anopheles gambiae is the principal vector of malaria, a disease that afflicts more than 500 million people and causes more than 1 million deaths each year. Tenfold shotgun sequence coverage was obtained from the PEST strain of A. gambiae and assembled into scaffolds that span 278 million base pairs. A total of 91% of the genome was organized in 303 scaffolds; the largest scaffold was 23.1 million base pairs. There was substantial genetic variation within this strain, and the apparent existence of two haplotypes of approximately equal frequency ("dual haplotypes") in a substantial fraction of the genome likely reflects the outbred nature of the PEST strain. The sequence produced a conservative inference of more than 400,000 single-nucleotide polymorphisms that showed a markedly bimodal density distribution. Analysis of the genome sequence revealed strong evidence for about 14,000 protein-encoding transcripts. Prominent expansions in specific families of proteins likely involved in cell adhesion and immunity were noted. An expressed sequence tag analysis of genes regulated by blood feeding provided insights into the physiological adaptations of a hematophagous insect.

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.

Genetic differentiation in the African malaria vector, Anopheles gambiae s.s., and the problem of taxonomic status

Of the seven recognized species of the Anopheles gambiae complex, A. gambiae s.s. is the most widespread and most important vector of malaria. It is becoming clear that, in parts of West Africa, this nominal species is not a single panmictic unit. We found that the internal transcribed spacer (ITS) of the X-linked rDNA has two distinct sequences with three fixed nucleotide differences; we detected no heterozygotes at these three sites, even in areas of sympatry of the two ITS types. The intergenic spacer (IGS) of this region also displays two distinct sequences that are in almost complete linkage disequilibrium with the distinct ITS alleles. We have designated these two types as S/type I and M/type II. These rDNA types correspond at least partly to the previously recognized chromosomal forms. Here we expand the geographic range of sampling to 251 individuals from 38 populations. Outside of West Africa, a single rDNA type, S/type I, corresponds to the Savanna chromosomal form. In West Africa, both types are often found in a single local sample. To understand if these findings might be due to unusual behavior of the rDNA region, we sequenced the same region for 46 A. arabiensis, a sympatric sibling species. No such distinct discontinuity was observed for this species. Autosomal inversions in one chromosome arm (2R), an insecticide resistance gene on 2L, and this single X-linked region indicate at least two genetically differentiated subpopulations of A. gambiae. Yet, rather extensive studies of other regions of the genome have failed to reveal genetic discontinuity. Evidently, incomplete genetic isolation exists within this single nominal species.

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.

Comparative genome analysis and pathway reconstruction

Pathway reconstruction builds on genome and biochemical data with the aim of reconstructing higher level interactions between identified enzymes in a specific genome, in particular the different enzyme pathways (species or individual/patient). Metabolite flow in a pathway is analyzed by different tools, such as elementary mode analysis. This reveals key enzymes and pharmacological targets in the enzyme network. An overview of bioinformatic tools and algorithms for these tasks, application examples and recent results from these techniques are presented. Target selection, drug development and optimization can all be sped up using these approaches.

A comparative genomic analysis of two distant diptera, the fruit fly, Drosophila melanogaster, and the malaria mosquito, Anopheles gambiae

Genome evolution entails changes in the DNA sequence of genes and intergenic regions, changes in gene numbers, and also changes in gene order along the chromosomes. Genes are reshuffled by chromosomal rearrangements such as deletions/insertions, inversions, translocations, and transpositions. Here we report a comparative study of genome organization in the main African malaria vector, Anopheles gambiae, relative to the recently determined sequence of the Drosophila melanogaster genome. The ancestral lines of these two dipteran insects are thought to have separated approximately 250 Myr, a long period that makes this genome comparison especially interesting. Sequence comparisons have identified 113 pairs of putative orthologs of the two species. Chromosomal mapping of orthologous genes reveals that each polytene chromosome arm has a homolog in the other species. Between 41% and 73% of the known orthologous genes remain linked in the respective homologous chromosomal arms, with the remainder translocated to various nonhomologous arms. Within homologous arms, gene order is extensively reshuffled, but a limited degree of conserved local synteny (microsynteny) can be recognized.