Genomic islands of speciation in Anopheles gambiae

Single-feature polymorphism discovery by computing probe affinity shape powers

Background

Single-feature polymorphism (SFP) discovery is a rapid and cost-effective approach to identify DNA polymorphisms. However, high false positive rates and/or low sensitivity are prevalent in previously described SFP detection methods. This work presents a new computing method for SFP discovery.

Results

The probe affinity differences and affinity shape powers formed by the neighboring probes in each probe set were computed into SFP weight scores. This method was validated by known sequence information and was comprehensively compared with previously-reported methods using the same datasets. A web application using this algorithm has been implemented for SFP detection. Using this method, we identified 364 SFPs in a barley near-isogenic line pair carrying either the wild type or the mutant uniculm2 (cul2) allele. Most of the SFP polymorphisms were identified on chromosome 6H in the vicinity of the Cul2 locus.

Conclusion

This SFP discovery method exhibits better performance in specificity and sensitivity over previously-reported methods. It can be used for other organisms for which GeneChip technology is available. The web-based tool will facilitate SFP discovery. The 364 SFPs discovered in a barley near-isogenic line pair provide a set of genetic markers for fine mapping and future map-based cloning of the Cul2 locus.

Comparative genomic and population genetic analyses indicate highly porous genomes and high levels of gene flow between divergent helianthus species

While speciation can be found in the presence of gene flow, it is not clear what impact this gene flow has on genome- and range-wide patterns of differentiation. Here we examine gene flow across the entire range of the common sunflower, H. annuus, its historically allopatric sister species H. argophyllus and a more distantly related, sympatric relative H. petiolaris. Analysis of genotypes at 26 microsatellite loci in 1015 individuals from across the range of the three species showed substantial introgression between geographically proximal populations of H. annuus and H. petiolaris, limited introgression between H. annuus and H. argophyllus, and essentially no gene flow between the allopatric pair, H. argophyllus and H. petiolaris. Analysis of sequence divergence levels among the three species in 1420 orthologs identified from EST databases identified a subset of loci showing extremely low divergence between H. annuus and H. petiolaris and extremely high divergence between the sister species H. annuus and H. argophyllus, consistent with introgression between H. annuus and H. petiolaris at these loci. Thus, at many loci, the allopatric sister species are more genetically divergent than the more distantly related sympatric species, which have exchanged genes across much of the genome while remaining morphologically and ecologically distinct.

Chromosomal inversions and species differences: when are genes affecting adaptive divergence and reproductive isolation expected to reside within inversions?

Many factors can promote speciation, and one which has received much attention is chromosomal inversions. A number of models propose that the recombination suppressing effects of inversions facilitate the maintenance of differences between interbreeding populations in genes affecting adaptive divergence and reproductive isolation. These models predict that such genes will disproportionately reside within inversions, rather than in collinear regions. This hypothesis has received some support, but exceptions exist. Additionally, the effects of known low levels of recombination within inversions on these models are uninvestigated. Here, simulations are used to compare the maintenance of genetic differences between populations following secondary contact and hybridization in different inversion models. We compare regions with no recombination within them to regions with low recombination and to collinear regions with free recombination. Our most general finding is that the low levels of recombination within an inversion often result in the loss of accentuated divergence in inverted regions compared to collinear ones. We conclude that inversions can facilitate the maintenance of species differences under some conditions, but that large or qualitative differences between inverted and collinear regions need not occur. We also find that strong selection facilitates maintenance of divergence in a manner analogous to inversions.

Rapid fixation of non-native alleles revealed by genome-wide SNP analysis of hybrid tiger salamanders

Background

Hybrid zones represent valuable opportunities to observe evolution in systems that are unusually dynamic and where the potential for the origin of novelty and rapid adaptation co-occur with the potential for dysfunction. Recently initiated hybrid zones are particularly exciting evolutionary experiments because ongoing natural selection on novel genetic combinations can be studied in ecological time. Moreover, when hybrid zones involve native and introduced species, complex genetic patterns present important challenges for conservation policy. To assess variation of admixture dynamics, we scored a large panel of markers in five wild hybrid populations formed when Barred Tiger Salamanders were introduced into the range of California Tiger Salamanders.

Results

At three of 64 markers, introduced alleles have largely displaced native alleles within the hybrid populations. Another marker (GNAT1) showed consistent heterozygote deficits in the wild, and this marker was associated with embryonic mortality in laboratory F2's. Other deviations from equilibrium expectations were idiosyncratic among breeding ponds, consistent with highly stochastic demographic effects.

Conclusion

While most markers retain native and introduced alleles in expected proportions, strong selection appears to be eliminating native alleles at a smaller set of loci. Such rapid fixation of alleles is detectable only in recently formed hybrid zones, though it might be representative of dynamics that frequently occur in nature. These results underscore the variable and mosaic nature of hybrid genomes and illustrate the potency of recombination and selection in promoting variable, and often unpredictable genetic outcomes. Introgression of a few, strongly selected introduced alleles should not necessarily affect the conservation status of California Tiger Salamanders, but suggests that genetically pure populations of this endangered species will be difficult to maintain.

High, clustered, nucleotide diversity in the genome of Anopheles gambiae revealed through pooled-template sequencing: implications for high-throughput genotyping protocols

Background

Association mapping approaches are dependent upon discovery and validation of single nucleotide polymorphisms (SNPs). To further association studies in Anopheles gambiae we conducted a major resequencing programme, primarily targeting regions within or close to candidate genes for insecticide resistance.

Results

Using two pools of mosquito template DNA we sequenced over 300 kbp across 660 distinct amplicons of the An. gambiae genome. Comparison of SNPs identified from pooled templates with those from individual sequences revealed a very low false positive rate. False negative rates were much higher and mostly resulted from SNPs with a low minor allele frequency. Pooled-template sequencing also provided good estimates of SNP allele frequencies. Allele frequency estimation success, along with false positive and negative call rates, improved significantly when using a qualitative measure of SNP call quality. We identified a total of 7062 polymorphic features comprising 6995 SNPs and 67 indels, with, on average, a SNP every 34 bp; a high rate of polymorphism that is comparable to other studies of mosquitoes. SNPs were significantly more frequent in members of the cytochrome p450 mono-oxygenases and carboxy/cholinesterase gene-families than in glutathione-S-transferases, other detoxification genes, and control genomic regions. Polymorphic sites showed a significantly clustered distribution, but the degree of SNP clustering (independent of SNP frequency) did not vary among gene families, suggesting that clustering of polymorphisms is a general property of the An. gambiae genome.

Conclusion

The high frequency and clustering of SNPs has important ramifications for the design of high-throughput genotyping assays based on allele specific primer extension or probe hybridisation. We illustrate these issues in the context of the design of Illumina GoldenGate assays.

The population genomics of trans-specific inversion polymorphisms in Anopheles gambiae

In the malaria mosquito Anopheles gambiae polymorphic chromosomal inversions may play an important role in adaptation to environmental variation. Recently, we used microarray-based divergence mapping combined with targeted resequencing to map nucleotide differentiation between alternative arrangements of the 2La inversion. Here, we applied the same technique to four different polymorphic inversions on the 2R chromosome of An. gambiae. Surprisingly, divergence was much lower between alternative arrangements for all 2R inversions when compared to the 2La inversion. For one of the rearrangements, 2Ru, we successfully mapped a very small region (approximately 100 kb) of elevated divergence. For the other three rearrangements, we did not identify any regions of significantly high divergence, despite ample independent evidence from natural populations of geographic clines and seasonal cycling, and stable heterotic polymorphisms in laboratory populations. If these inversions are the targets of selection as hypothesized, we suggest that divergence between rearrangements may have escaped detection due to retained ancestral polymorphism in the case of the youngest 2R rearrangements and to extensive gene flux in the older 2R inversion systems that segregate in both An. gambiae and its sibling species An. arabiensis.

Adaptive radiations: from field to genomic studies

Adaptive radiations were central to Darwin's formation of his theory of natural selection, and today they are still the centerpiece for many studies of adaptation and speciation. Here, we review the advantages of adaptive radiations, especially recent ones, for detecting evolutionary trends and the genetic dissection of adaptive traits. We focus on Aquilegia as a primary example of these advantages and highlight progress in understanding the genetic basis of flower color. Phylogenetic analysis of Aquilegia indicates that flower color transitions proceed by changes in the types of anthocyanin pigments produced or their complete loss. Biochemical, crossing, and gene expression studies have provided a wealth of information about the genetic basis of these transitions in Aquilegia. To obtain both enzymatic and regulatory candidate genes for the entire flavonoid pathway, which produces anthocyanins, we used a combination of sequence searches of the Aquilegia Gene Index, phylogenetic analyses, and the isolation of novel sequences by using degenerate PCR and RACE. In total we identified 34 genes that are likely involved in the flavonoid pathway. A number of these genes appear to be single copy in Aquilegia and thus variation in their expression may have been key for floral color evolution. Future studies will be able to use these sequences along with next-generation sequencing technologies to follow expression and sequence variation at the population level. The genetic dissection of other adaptive traits in Aquilegia should also be possible soon as genomic resources such as whole-genome sequencing become available.

Inferring selection in the Anopheles gambiae species complex: an example from immune-related serine protease inhibitors

BACKGROUND

Mosquitoes of the Anopheles gambiae species complex are the primary vectors of human malaria in sub-Saharan Africa. Many host genes have been shown to affect Plasmodium development in the mosquito, and so are expected to engage in an evolutionary arms race with the pathogen. However, there is little conclusive evidence that any of these mosquito genes evolve rapidly, or show other signatures of adaptive evolution.

METHODS

Three serine protease inhibitors have previously been identified as candidate immune system genes mediating mosquito-Plasmodium interaction, and serine protease inhibitors have been identified as hot-spots of adaptive evolution in other taxa. Population-genetic tests for selection, including a recent multi-gene extension of the McDonald-Kreitman test, were applied to 16 serine protease inhibitors and 16 other genes sampled from the An. gambiae species complex in both East and West Africa.

RESULTS

Serine protease inhibitors were found to show a marginally significant trend towards higher levels of amino acid diversity than other genes, and display extensive genetic structuring associated with the 2La chromosomal inversion. However, although serpins are candidate targets for strong parasite-mediated selection, no evidence was found for rapid adaptive evolution in these genes.

CONCLUSION

It is well known that phylogenetic and population history in the An. gambiae complex can present special problems for the application of standard population-genetic tests for selection, and this may explain the failure of this study to detect selection acting on serine protease inhibitors. The pitfalls of uncritically applying these tests in this species complex are highlighted, and the future prospects for detecting selection acting on the An. gambiae genome are discussed.

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.

SNEP: Simultaneous detection of nucleotide and expression polymorphisms using Affymetrix GeneChip

BACKGROUND

High-density short oligonucleotide microarrays are useful tools for studying biodiversity, because they can be used to investigate both nucleotide and expression polymorphisms. However, when different strains (or species) produce different signal intensities after mRNA hybridization, it is not easy to determine whether the signal intensities were affected by nucleotide or expression polymorphisms. To overcome this difficulty, nucleotide and expression polymorphisms are currently examined separately.

RESULTS

We have developed SNEP, a new method that allows simultaneous detection of both nucleotide and expression polymorphisms. SNEP involves a robust statistical procedure based on the idea that a nucleotide polymorphism observed at the probe level can be regarded as an outlier, because the nucleotide polymorphism can reduce the hybridization signal intensity. To investigate the performance of SNEP, we used three species: barley, rice and mice. In addition to the publicly available barley data, we obtained new rice and mouse data from the strains with available genome sequences. The sensitivity and false positive rate of nucleotide polymorphism detection were estimated based on the sequence information. The robustness of expression polymorphism detection against nucleotide polymorphisms was also investigated.

CONCLUSION

SNEP performed well regardless of the genome size and showed a better performance for nucleotide polymorphism detection, when compared with other previously proposed methods. The R-software 'SNEP' is available at http://www.ism.ac.jp/~fujisawa/SNEP/.

Distinctiveness in the face of gene flow: hybridization between specialist and generalist gartersnakes

Patterns of divergence and polymorphism across hybrid zones can provide important clues as to their origin and maintenance. Unimodal hybrid zones or hybrid swarms are composed predominantly of recombinant individuals whose genomes are patchworks of alleles derived from each parental lineage. In contrast, bimodal hybrid zones contain few identifiable hybrids; most individuals fall within distinct genetic clusters. Distinguishing between hybrid swarms and bimodal hybrid zones can be important for taxonomic and conservation decisions regarding the status and value of hybrid populations. In addition, the causes of bimodality are important in understanding the generation and maintenance of biological diversity. For example, are distinct clusters mostly reproductively isolated and co-adapted gene complexes, or can distinctiveness be maintained by a few 'genomic islands' despite rampant gene flow across much of the genome? Here we focus on three patterns of distinctiveness in the face of gene flow between gartersnake taxa in the Great Lakes region of North America. Bimodality, the persistence of distinct clusters of genotypes, requires strong barriers to gene flow and supports recognition of distinct specialist (Thamnophis butleri) and generalist (Thamnophis radix) taxa. Concordance of DNA-based clusters with morphometrics supports the hypothesis that trophic morphology is a key component of divergence. Finally, disparity in the level of differentiation across molecular markers (amplified fragment length polymorphisms) indicates that distinctiveness is maintained by strong selection on a few traits despite high gene flow currently or in the recent past.

Array-based genotyping in S.cerevisiae using semi-supervised clustering

MOTIVATION

Microarrays provide an accurate and cost-effective method for genotyping large numbers of individuals at high resolution. The resulting data permit the identification of loci at which genetic variation is associated with quantitative traits, or fine mapping of meiotic recombination, which is a key determinant of genetic diversity among individuals. Several issues inherent to short oligonucleotide arrays -- cross-hybridization, or variability in probe response to target -- have the potential to produce genotyping errors. There is a need for improved statistical methods for array-based genotyping.

RESULTS

We developed ssGenotyping (ssG), a multivariate, semi-supervised approach for using microarrays to genotype haploid individuals at thousands of polymorphic sites. Using a meiotic recombination dataset, we show that ssG is more accurate than existing supervised classification methods, and that it produces denser marker coverage. The ssG algorithm is able to fit probe-specific affinity differences and to detect and filter spurious signal, permitting high-confidence genotyping at nucleotide resolution. We also demonstrate that oligonucleotide probe response depends significantly on genomic background, even when the probe's specific target sequence is unchanged. As a result, supervised classifiers trained on reference strains may not generalize well to diverged strains; ssG's semi-supervised approach, on the other hand, adapts automatically.

Learning about modes of speciation by computational approaches

How often do the early stages of speciation occur in the presence of gene flow? To address this enduring question, a number of recent papers have used computational approaches, estimating parameters of simple divergence models from multilocus polymorphism data collected in closely related species. Applications to a variety of species have yielded extensive evidence for migration, with the results interpreted as supporting the widespread occurrence of parapatric speciation. Here, we conduct a simulation study to assess the reliability of such inferences, using a program that we recently developed MCMC estimation of the isolation-migration model allowing for recombination (MIMAR) as well as the program isolation-migration (IM) of Hey and Nielsen (2004). We find that when one of many assumptions of the isolation-migration model is violated, the methods tend to yield biased estimates of the parameters, potentially lending spurious support for allopatric or parapatric divergence. More generally, our results highlight the difficulty in drawing inferences about modes of speciation from the existing computational approaches alone.

Recombination and speciation: loci near centromeres are more differentiated than loci near telomeres between subspecies of the European rabbit (Oryctolagus cuniculus)

Recent empirical and theoretical studies suggest that regions of restricted recombination play an important role in the formation of new species. To test this idea, we studied nucleotide variation in two parapatric subspecies of the European rabbit (Oryctolagus cuniculus). We surveyed five loci near centromeres, where recombination is expected to be suppressed, and five loci near telomeres, where recombination is expected to be higher. We analyzed this multilocus data set using a divergence-with-gene flow framework and we report three main findings. First, we estimated that these subspecies diverged approximately 1.8 MYA and maintained large effective population sizes (O. c. algirus N(e) approximately 1,600,000 and O. c. cuniculus N(e) approximately 780,000). Second, we rejected a strict allopatric model of divergence without gene flow; instead, high rates of gene flow were inferred in both directions. Third, we found different patterns between loci near centromeres and loci near telomeres. Loci near centromeres exhibited higher levels of linkage disequilibrium than loci near telomeres. In addition, while all loci near telomeres showed little differentiation between subspecies, three of five loci near centromeres showed strong differentiation. These results support a view of speciation in which regions of low recombination can facilitate species divergence in the presence of gene flow.

Y not introgress? Insights into the genetics of speciation in European rabbits

Hybridization between genetically divergent populations is common in nature. By comparing the rate of gene flow throughout the genome, loci that impede genetic mixing, and therefore underlie reproductive isolation, can be identified, providing key insights into the process of speciation. In a previous issue of Molecular Ecology, Geraldes et al. (2008) report the geographical distribution of Y chromosome lineages in wild populations of rabbits from across the Iberian Peninsula and southern France. Y chromosomes showed a high level of differentiation between rabbit subspecies on either side of a hybrid zone, despite clear evidence for gene flow at other loci. This provocative result suggests a role for the Y chromosome in reproductive isolation, and adds to a growing list of nascent species with mosaic patterns of differentiation across their genomes.

Divergence in the face of gene flow: the case of two newts (amphibia: salamandridae)

Understanding the process of divergence requires the quantitative characterization of patterns of gene flow between diverging taxa. New and powerful coalescent-based methods give insight into these processes in unprecedented details by enabling the reconstruction of the temporal distribution of past gene flow. Here, we use sequence variation at eight nuclear markers and mitochondrial DNA (mtDNA) in multiple populations to study diversity, divergence, and gene flow between two subspecies of a salamander, the smooth newt (Lissotriton vulgaris kosswigi and Lissotriton vulgaris vulgaris) in Turkey. The ranges of both subspecies encompass mainly the areas of this important glacial refugial area. Populations in refugia where species have been present for a long time and differentiated in situ should better preserve the record of past gene flow than young populations in postglacial expansion areas. Sequence diversity in both subspecies was substantial (nuclear pi(sil) = 0.69% and 1.31%). We detected long-term demographic stability in these refugial populations with large effective population sizes (N(e)) of the order of 1.5-3 x 10(5) individuals. Gene trees and the isolation with migration (IM) analysis complemented by tests of nested IM models showed that despite deep, pre-Pleistocene divergence of the studied newts, asymmetric introgression from vulgaris to kosswigi has occurred, with signatures of recent gene flow in mtDNA and an anonymous nuclear marker, and evidence for more ancient introgression in nuclear introns. The distribution of migration times raises the intriguing possibility that even the initial divergence may have occurred in the face of gene flow.

Hybridization, ecological races and the nature of species: empirical evidence for the ease of speciation

Species are generally viewed by evolutionists as 'real' distinct entities in nature, making speciation appear difficult. Charles Darwin had originally promoted a very different uniformitarian view that biological species were continuous with 'varieties' below the level of species and became distinguishable from them only when divergent natural selection led to gaps in the distribution of morphology. This Darwinian view on species came under immediate attack, and the consensus among evolutionary biologists today appears to side more with the ideas of Ernst Mayr and Theodosius Dobzhansky, who argued 70 years ago that Darwin was wrong about species. Here, I show how recent genetic studies of supposedly well-behaved animals, such as insects and vertebrates, including our own species, have supported the existence of the Darwinian continuum between varieties and species. Below the level of species, there are well-defined ecological races, while above the level of species, hybridization still occurs, and may often lead to introgression and, sometimes, hybrid speciation. This continuum is evident, not only across vast geographical regions, but also locally in sympatry. The existence of this continuum provides good evidence for gradual evolution of species from ecological races and biotypes, to hybridizing species and, ultimately, to species that no longer cross. Continuity between varieties and species not only provides an excellent argument against creationism, but also gives insight into the process of speciation. The lack of a hiatus between species and ecological races suggests that speciation may occur, perhaps frequently, in sympatry, and the abundant intermediate stages suggest that it is happening all around us. Speciation is easy!

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.

Exploring the origin and degree of genetic isolation of Anopheles gambiae from the islands of São Tomé and Príncipe, potential sites for testing transgenic-based vector control

The evolutionary processes at play between island and mainland populations of the malaria mosquito vector Anopheles gambiae sensu stricto are of great interest as islands may be suitable sites for preliminary application of transgenic-based vector control strategies. São Tomé and Príncipe, located off the West African coast, have received such attention in recent years. This study investigates the degree of isolation of An. gambiae s.s. populations between these islands and the mainland based on mitochondrial and ribosomal DNA molecular data. We identify possible continental localities from which these island populations derived. For these purposes, we used F_ST values, haplotype networks, and nested clade analysis to estimate migration rates and patterns. Haplotypes from both markers are geographically widespread across the African continent. Results indicate that the populations from São Tomé and Príncipe are relatively isolated from continental African populations, suggesting they are promising sites for test releases of transgenic individuals. These island populations are possibly derived from two separate continental migrations. This result is discussed in the context of the history of the African slave trade with respect to São Tomé and Príncipe.

The evolution of TEP1, an exceptionally polymorphic immunity gene in Anopheles gambiae

Background

Host-parasite coevolution can result in balancing selection, which maintains genetic variation in the susceptibility of hosts to parasites. It has been suggested that variation in a thioester-containing protein called TEP1 (AGAP010815) may alter the ability of Anopheles mosquitoes to transmit Plasmodium parasites, and high divergence between alleles of this gene suggests the possible action of long-term balancing selection. We studied whether TEP1 is a case of an ancient balanced polymorphism in an animal immune system.

Results

We found evidence that the high divergence between TEP1 alleles is the product of genetic exchange between TEP1 and other TEP loci, i.e. gene conversion. Additionally, some TEP1 alleles showed unexpectedly low variability.

Conclusion

The TEP1 gene appears to be a chimera produced from at least two other TEP loci, and the divergence between TEP1 alleles is probably not caused by long-term balancing selection, but is instead due to two independent gene conversion events from one of these other genes. Nevertheless, TEP1 still shows evidence of natural selection, in particular there appears to have been recent changes in the frequency of alleles that has diminished polymorphism within each allelic class. Although the selective force driving this dynamic was not identified, given that susceptibility to Plasmodium parasites is known to be associated with allelic variation in TEP1, these changes in allele frequencies could alter the vectoring capacity of populations.

Review. The strength and genetic basis of reproductive isolating barriers in flowering plants

Speciation is characterized by the evolution of reproductive isolation between two groups of organisms. Understanding the process of speciation requires the quantification of barriers to reproductive isolation, dissection of the genetic mechanisms that contribute to those barriers and determination of the forces driving the evolution of those barriers. Through a comprehensive analysis involving 19 pairs of plant taxa, we assessed the strength and patterns of asymmetry of multiple prezygotic and postzygotic reproductive isolating barriers. We then reviewed contemporary knowledge of the genetic architecture of reproductive isolation and the relative role of chromosomal and genic factors in intrinsic postzygotic isolation. On average, we found that prezygotic isolation is approximately twice as strong as postzygotic isolation, and that postmating barriers are approximately three times more asymmetrical in their action than premating barriers. Barriers involve a variable number of loci, and chromosomal rearrangements may have a limited direct role in reproductive isolation in plants. Future research should aim to understand the relationship between particular genetic loci and the magnitude of their effect on reproductive isolation in nature, the geographical scale at which plant speciation occurs, and the role of different evolutionary forces in the speciation process.

Anopheles gambiae complex along The Gambia river, with particular reference to the molecular forms of An. gambiae s.s

BACKGROUND

The geographic and temporal distribution of M and S molecular forms of the major Afrotropical malaria vector species Anopheles gambiae s.s. at the western extreme of their range of distribution has never been investigated in detail.

MATERIALS AND METHODS

Collections of indoor-resting An. gambiae s.l. females were carried out along a ca. 400 km west to east transect following the River Gambia from the western coastal region of The Gambia to south-eastern Senegal during 2005 end of rainy season/early dry season and the 2006 rainy season. Specimens were identified to species and molecular forms by PCR-RFLP and the origin of blood-meal of fed females was determined by ELISA test.

RESULTS

Over 4,000 An. gambiae s.l. adult females were collected and identified, 1,041 and 3,038 in 2005 and 2006, respectively. M-form was mainly found in sympatry with Anopheles melas and S-form in the western part of the transect, and with Anopheles arabiensis in the central part. S-form was found to prevail in rural Sudan-Guinean savannah areas of Eastern Senegal, in sympatry with An. arabiensis. Anopheles melas and An. arabiensis relative frequencies were generally lower in the rainy season samples, when An. gambiae s.s. was prevailing. No large seasonal fluctuations were observed for M and S-forms. In areas where both M and S were recorded, the frequency of hybrids between them ranged from to 0.6% to 7%.

DISCUSSION

The observed pattern of taxa distribution supports the hypothesis of a better adaptation of M-form to areas characterized by water-retaining alluvial deposits along the Gambia River, characterized by marshy vegetation, mangrove woods and rice cultivations. In contrast, the S-form seems to be better adapted to free-draining soil, covered with open woodland savannah or farmland, rich in temporary larval breeding sites characterizing mainly the eastern part of the transect, where the environmental impact of the Gambia River is much less profound and agricultural activities are mainly rain-dependent. Very interestingly, the observed frequency of hybridization between the molecular forms along the whole transect was much higher than has been reported so far for other areas.

CONCLUSION

The results support a bionomic divergence between the M and S-forms, and suggest that the western extreme of An. gambiae s.s. geographical distribution may represent an area of higher-than-expected hybridization between the two molecular forms.

Genomic analysis of differentiation between soil types reveals candidate genes for local adaptation in Arabidopsis lyrata

Serpentine soil, which is naturally high in heavy metal content and has low calcium to magnesium ratios, comprises a difficult environment for most plants. An impressive number of species are endemic to serpentine, and a wide range of non-endemic plant taxa have been shown to be locally adapted to these soils. Locating genomic polymorphisms which are differentiated between serpentine and non-serpentine populations would provide candidate loci for serpentine adaptation. We have used the Arabidopsis thaliana tiling array, which has 2.85 million probes throughout the genome, to measure genetic differentiation between populations of Arabidopsis lyrata growing on granitic soils and those growing on serpentinic soils. The significant overrepresentation of genes involved in ion transport and other functions provides a starting point for investigating the molecular basis of adaptation to soil ion content, water retention, and other ecologically and economically important variables. One gene in particular, calcium-exchanger 7, appears to be an excellent candidate gene for adaptation to low Ca∶Mg ratio in A. lyrata.

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.

Ecological reproductive isolation of coast and inland races of Mimulus guttatus

Adaptive divergence due to habitat differences is thought to play a major role in formation of new species. However it is rarely clear the extent to which individual reproductive isolating barriers related to habitat differentiation contribute to total isolation. Furthermore, it is often difficult to determine the specific environmental variables that drive the evolution of those ecological barriers, and the geographic scale at which habitat-mediated speciation occurs. Here, we address these questions through an analysis of the population structure and reproductive isolation between coastal perennial and inland annual forms of the yellow monkeyflower, Mimulus guttatus. We found substantial morphological and molecular genetic divergence among populations derived from coast and inland habitats. Reciprocal transplant experiments revealed nearly complete reproductive isolation between coast and inland populations mediated by selection against immigrants and flowering time differences, but not postzygotic isolation. Our results suggest that selection against immigrants is a function of adaptations to seasonal drought in inland habitat and to year round soil moisture and salt spray in coastal habitat. We conclude that the coast and inland populations collectively comprise distinct ecological races. Overall, this study suggests that adaptations to widespread habitats can lead to the formation of reproductively isolated species.

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 Yin and Yang of linkage disequilibrium: mapping of genes and nucleotides conferring insecticide resistance in insect disease vectors

Genetic technologies developed in the last 20 years have lead to novel and exciting methods to identify genes and specific nucleotides within genes that control phenotypes in field collected organisms. In this review we define and explain two of these methods: linkage disequilibrium (LD) mapping and quantitative trait nucleotide (QTN) mapping. The power to detect valid genotype-phenotype associations with LD or QTN mapping depends critically on the extent to which segregating sites in a genome assort independently. LD mapping depends on markers being in disequilibrium with the genes that condition expression of the phenotype. In contrast, QTN mapping depends critically upon most proximal loci being at equilibrium. We show that both patterns actually exist in the genome of Anapheles gambiae, the most important malaria vector in sub-Saharan Africa while segregating sites appear to be largely in equilibrium throughout the genome of Aedes aegypti, the vector of Dengue and Yellow fever flaviviruses. We discuss additional approaches that will be needed to identify genes and nucleotides that control phenotypes in field collected organisms, focusing specifically on ongoing studies of genes conferring resistance to insecticides.

Assessment of alternative mating strategies in Anopheles gambiae: Does mating occur indoors?

Mating in Anopheles gambiae has been observed only in outdoor swarms. Here we evaluate whether mating also occurs indoors. Mark-release-recapture of virgin males and females in natural houses showed that mating occurred over a single day even when mosquitoes can leave the house through exit traps and without adaptation to laboratory conditions. In these experiments, insemination rate in the M molecular form of An. gambiae (and An. arabiensis) was higher than that of the S form (15 versus 6%). Under these conditions, smaller females of the M form mated more frequently than larger females of that form. Sampling mosquitoes throughout the day showed that both sexes enter houses around sunrise and leave around sunset, staying indoors together from dawn to dusk. In an area dominated by the M form, the daily rate of insemination in samples from exit traps was approximately 5% higher than in those from entry traps, implying that mating occurred indoors. Importantly, frequency of cross mating between the molecular forms was as high as that between members of the same form, indicating that, indoors, assortative mating breaks down. Altogether, these results suggest that indoor mating is an alternative mating strategy of the M molecular form of An. gambiae. Because naturally occurring mating couples have not yet been observed indoors, this conclusion awaits validation.

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.

Genetic hitchhiking in a subdivided population of Mytilus edulis

BACKGROUND

Few models of genetic hitchhiking in subdivided populations have been developed and the rarity of empirical examples is even more striking. We here provide evidences of genetic hitchhiking in a subdivided population of the marine mussel Mytilus edulis. In the Bay of Biscay (France), a patch of M. edulis populations happens to be separated from its North Sea conspecifics by a wide region occupied only by the sister species M. galloprovincialis. Although genetic differentiation between the two M. edulis regions is largely non-significant at ten marker loci (average FST~0.007), a strong genetic differentiation is observed at a single locus (FST = 0.25). We validated the outlier status of this locus, and analysed DNA sequence polymorphism in order to identify the nature of the selection responsible for the unusual differentiation.

RESULTS

We first showed that introgression of M. galloprovincialis alleles was very weak in both populations and did not significantly affect their differentiation. Secondly, we observed the genetic signature of a selective sweep within both M. edulis populations in the form of a star-shaped clade of alleles. This clade was nearly fixed in the North Sea and was segregating at a moderate frequency in the Bay of Biscay, explaining their genetic differentiation. Incomplete fixation reveals that selection was not direct on the locus but that the studied sequence recombined with a positively selected allele at a linked locus while it was on its way to fixation. Finally, using a deterministic model we showed that the wave of advance of a favourable allele at a linked locus, when crossing a strong enough barrier to gene flow, generates a step in neutral allele frequencies comparable to the step observed between the two M. edulis populations at the outlier locus. In our case, the position of the barrier is now materialised by a large patch of heterospecific M. galloprovincialis populations.

CONCLUSION

High FST outlier loci are usually interpreted as being the consequence of ongoing divergent local adaptation. Combining models and data we show that among-population differentiation can also dramatically increase following a selective sweep in a structured population. Our study illustrates how a striking geographical pattern of neutral diversity can emerge from past indirect hitchhiking selection in a structured population. NOTE: Nucleotide sequences reported in this paper are available in the GenBanktrade mark database under the accession numbers EU684165 - EU684228.

SNP discovery and molecular evolution in Anopheles gambiae, with special emphasis on innate immune system

Background

Anopheles innate immunity affects Plasmodium development and is a potential target of innovative malaria control strategies. The extent and distribution of nucleotide diversity in immunity genes might provide insights into the evolutionary forces that condition pathogen-vector interactions. The discovery of polymorphisms is an essential step towards association studies of susceptibility to infection.

Results

We sequenced coding fragments of 72 immune related genes in natural populations of Anopheles gambiae and of 37 randomly chosen genes to provide a background measure of genetic diversity across the genome. Mean nucleotide diversity (π) was 0.0092 in the A. gambiae S form, 0.0076 in the M form and 0.0064 in A. arabiensis. Within each species, no statistically significant differences in mean nucleotide diversity were detected between immune related and non immune related genes. Strong purifying selection was detected in genes of both categories, presumably reflecting strong functional constraints.

Conclusion

Our results suggest similar patterns and rates of molecular evolution in immune and non-immune genes in A. gambiae. The 3,214 Single Nucleotide Polymorphisms (SNPs) that we identified are the first large set of Anopheles SNPs from fresh, field-collected material and are relevant markers for future phenotype-association studies.

Evidence of introgression of the ace-1(R) mutation and of the ace-1 duplication in West African Anopheles gambiae s. s

BACKGROUND

The role of inter-specific hybridisation is of particular importance in mosquito disease vectors for predicting the evolution of insecticide resistance. Two molecular forms of Anopheles gambiae s.s., currently recognized as S and M taxa, are considered to be incipient sibling species. Hybrid scarcity in the field was suggested that differentiation of M and S taxa is maintained by limited or absent gene flow. However, recent studies have revealed shared polymorphisms within the M and S forms, and a better understanding of the occurrence of gene flow is needed. One such shared polymorphism is the G119S mutation in the ace-1 gene (which is responsible for insecticide resistance); this mutation has been described in both the M and S forms of A. gambiae s.s.

METHODS AND RESULTS

To establish whether the G119S mutation has arisen independently in each form or by genetic introgression, we analysed coding and non-coding sequences of ace-1 alleles in M and S mosquitoes from representative field populations. Our data revealed many polymorphic sites shared by S and M forms, but no diversity was associated with the G119S mutation. These results indicate that the G119S mutation was a unique event and that genetic introgression explains the observed distribution of the G119S mutation within the two forms. However, it was impossible to determine from our data whether the mutation occurred first in the S form or in the M form. Unexpectedly, sequence analysis of some resistant individuals revealed a duplication of the ace-1 gene that was observed in both A. gambiae s.s. M and S forms. Again, the distribution of this duplication in the two forms most likely occurred through introgression.

CONCLUSIONS

These results highlight the need for more research to understand the forces driving the evolution of insecticide resistance in malaria vectors and to regularly monitor resistance in mosquito populations of Africa.

Genetics of intrinsic postzygotic isolation in a circumpolar plant species, Draba nivalis (Brassicaceae)

Sterility barriers, ranging from incomplete to fully developed, were recently demonstrated within taxonomic species of the genus Draba, suggesting the existence of numerous, cryptic biological species. Because these taxa are predominately selfers and of Pleistocene origin, it was concluded that hybrid sterility evolved quickly and possibly by genetic drift. Here we used genetic mapping and QTL analyses to determine the genetic basis of hybrid sterility between geographically distant populations of one of these taxonomic species, Draba nivalis. Fifty microsatellite loci were mapped, and QTL analyses identified five loci underlying seed fertility and two underlying pollen fertility. Four of five seed fertility QTLs reduced fertility in heterozygotes, an observation most consistent with drift-based fixation of underdominant sterility loci. However, several nuclear-nuclear interactions were also found, including two that acted like reciprocal translocations with lowest fitness in double heterozygotes, and two that had a pattern of fitness consistent with Bateson-Dobzhansky-Muller incompatibilities. In contrast, pollen fertility QTLs exhibited additive inheritance, with lowest fertility associated with the paternal allele, a pattern of inheritance suggestive of cytonuclear incompatibilities. The results imply that multiple genetic mechanisms underlie the rapid evolution of reproductive barriers in Draba.

Gene flow and the genealogical history of Heliconius heurippa

Background

The neotropical butterfly Heliconius heurippa has a hybrid colour pattern, which also contributes to reproductive isolation, making it a likely example of hybrid speciation. Here we used phylogenetic and coalescent-based analyses of multilocus sequence data to investigate the origin of H. heurippa.

Results

We sequenced a mitochondrial region (CoI and CoII), a sex-linked locus (Tpi) and two autosomal loci (w and sd) from H. heurippa and the putative parental species, H. cydno and H. melpomene. These were analysed in combination with data from two previously sequenced autosomal loci, Dll and Inv. H. heurippa was monophyletic at mtDNA and Tpi, but showed a shared distribution of alleles derived from both parental lineages at all four autosomal loci. Estimates of genetic differentiation showed that H. heurippa is closer to H. cydno at mtDNA and three autosomal loci, intermediate at Tpi, and closer to H. melpomene at Dll. Using coalescent simulations with the Isolation-Migration model (IM), we attempted to establish the incidence of gene flow in the origin of H. heurippa. This analysis suggested that ongoing introgression is frequent between all three species and variable in extent between loci.

Conclusion

Introgression, which is a necessary precursor of hybrid speciation, seems to have also blurred the coalescent history of these species. The origin of Heliconius heurippa may have been restricted to introgression of few colour pattern genes from H. melpomene into the H. cydno genome, with little evidence of genomic mosaicism.

Genomic analysis of adaptive differentiation in Drosophila melanogaster

Drosophila melanogaster shows clinal variation along latitudinal transects on multiple continents for several phenotypes, allozyme variants, sequence variants, and chromosome inversions. Previous investigation suggests that many such clines are due to spatially varying selection rather than demographic history, but the genomic extent of such selection is unknown. To map differentiation throughout the genome, we hybridized DNA from temperate and subtropical populations to Affymetrix tiling arrays. The dense genomic sampling of variants and low level of linkage disequilibrium in D. melanogaster enabled identification of many small, differentiated regions. Many regions are differentiated in parallel in the United States and Australia, strongly supporting the idea that they are influenced by spatially varying selection. Genomic differentiation is distributed nonrandomly with respect to gene function, even in regions differentiated on only one continent, providing further evidence for the role of selection. These data provide candidate genes for phenotypes known to vary clinally and implicate interesting new processes in genotype-by-environment interactions, including chorion proteins, proteins regulating meiotic recombination and segregation, gustatory and olfactory receptors, and proteins affecting synaptic function and behavior. This portrait of differentiation provides a genomic perspective on adaptation and the maintenance of variation through spatially varying selection.

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.

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.

Ancient hybridization and mitochondrial capture between two species of chipmunks

Models that posit speciation in the face of gene flow are replacing classical views that hybridization is rare between animal species. We use a multilocus approach to examine the history of hybridization and gene flow between two species of chipmunks (Tamias ruficaudus and T. amoenus). Previous studies have shown that these species occupy different ecological niches and have distinct genital bone morphologies, yet appear to be incompletely isolated reproductively in multiple areas of sympatry. We compared data from four sequenced nuclear loci and from seven microsatellite loci to published cytochrome b sequences. Interspecific gene flow was primarily restricted to introgression of the T. ruficaudus mitochondrial genome into a sympatric subspecies of T. amoenus, T. a. canicaudus, with the four sequenced nuclear loci showing little to no interspecific allele sharing. Microsatellite data were consistent with high levels of differentiation between the species and also showed no current gene flow between broadly sympatric populations of T. a. canicaudus and T. ruficaudus. Coalescent analyses date the mtDNA introgression event from the mid-Pleistocene to late Pliocene. Overall, these data indicate that introgression has had a minimal impact on the nuclear genomes of T. amoenus and T. ruficaudus despite multiple independent hybridization events. Our findings challenge long-standing assumptions on patterns of reproductive isolation in chipmunks and suggest that there may be other examples of hybridization among the 23 species of Tamias that occur in western North America.

Demographic processes shaping genetic variation

Demographic processes modulate genome-wide levels and patterns of genetic variation via impacting effective population size independently of natural selection. Such processes include the perturbation of population distributions from external events shaping habitat landscape and internal factors shaping the probability of contemporaneous alleles in a population (coalescence). Several patterns have recently emerged: spatial and temporal heterogeneity in population structure have different influences on the persistence of new mutations and genetic variation, multi-locus analyses indicate that gene flow continues to occur during speciation and the incorporation of demographic processes into models of molecular evolution and association genetics approaches has improved statistical power to detect deviations from neutral-equilibrium expectations and decreased false positive rates.

Population genetics of speciation in two closely related wild tomatoes (Solanum section Lycopersicon)

We present a multilocus sequencing study to assess patterns of polymorphism and divergence in the closely related wild tomato species, Solanum peruvianum and S. chilense (Solanum section Lycopersicon, Solanaceae). The data set comprises seven mapped nuclear loci (approximately 9.3 kb of analyzed sequence across loci) and four local population samples per species that cover much of the species' range (between 80 and 88 sequenced alleles across both species). We employ the analytical framework of divergence population genetics (DPG) in evaluating the utility of the "isolation" model of speciation to explain observed patterns of polymorphism and divergence. Whereas the isolation model is not rejected by goodness-of-fit criteria established via coalescent simulations, patterns of intragenic linkage disequilibrium provide evidence for postdivergence gene flow at two of the seven loci. These results suggest that speciation occurred under residual gene flow, implying that natural selection is one of the evolutionary forces driving the divergence of these tomato species. This inference is fully consistent with their recent divergence, conservatively estimated to be <or=0.55 million years. We discuss possible biases in the demographic parameter estimates due to the current restriction of DPG algorithms to panmictic species.

Comparing whole genomes using DNA microarrays

The rapid accumulation of complete genomic sequences offers the opportunity to carry out an analysis of inter- and intra-individual genome variation within a species on a routine basis. Sequencing whole genomes requires resources that are currently beyond those of a single laboratory and therefore it is not a practical approach for resequencing hundreds of individual genomes. DNA microarrays present an alternative way to study differences between closely related genomes. Advances in microarray-based approaches have enabled the main forms of genomic variation (amplifications, deletions, insertions, rearrangements and base-pair changes) to be detected using techniques that are readily performed in individual laboratories using simple experimental approaches.

Hybrid zone origins, species boundaries, and the evolution of wing-pattern diversity in a polytypic species complex of North American admiral butterflies (Nymphalidae: Limenitis)

Hybrid zones present opportunities to study the effects of gene flow, selection, and recombination in natural populations and, thus, provide insights into the genetic and phenotypic changes that occur early in speciation. Here we investigate a hybrid zone between mimetic (Limenitis arthemis astyanax) and nonmimetic (Limenitis arthemis arthemis) populations of admiral butterflies using DNA sequence variation from mtDNA and seven nuclear gene loci. We find three distinct mitochondrial clades within this complex, and observe a strong overall concordance between wing-pattern phenotypes and mitochondrial variation. Nuclear gene genealogies, in contrast, revealed no evidence of exclusivity for either wing-pattern phenotype, suggesting incomplete barriers to gene exchange and/or insufficient time for lineage sorting. Coalescent simulations indicate that gene flow between these two subspecies is highly asymmetric, with the majority of migration occurring from mimetic into nonmimetic populations. Selective sweeps of alleles responsible for mimetic phenotypes may have occurred more than once when mimetic and nonmimetic Limenitis occurred together in the presence of the model (Battus philenor).

A test of the chromosomal theory of ecotypic speciation in Anopheles gambiae

The role of chromosomal inversions in speciation has long been of interest to evolutionists. Recent quantitative modeling has stimulated reconsideration of previous conceptual models for chromosomal speciation. Anopheles gambiae, the most important vector of human malaria, carries abundant chromosomal inversion polymorphism nonrandomly associated with ecotypes that mate assortatively. Here, we consider the potential role of paracentric inversions in promoting speciation in A. gambiae via "ecotypification," a term that refers to differentiation arising from local adaptation. In particular, we focus on the Bamako form, an ecotype characterized by low inversion polymorphism and fixation of an inversion, 2Rj, that is very rare or absent in all other forms of A. gambiae. The Bamako form has a restricted distribution by the upper Niger River and its tributaries that is associated with a distinctive type of larval habitat, laterite rock pools, hypothesized to be its optimal breeding site. We first present computer simulations to investigate whether the population dynamics of A. gambiae are consistent with chromosomal speciation by ecotypification. The models are parameterized using field observations on the various forms of A. gambiae that exist in Mali, West Africa. We then report on the distribution of larvae of this species collected from rock pools and more characteristic breeding sites nearby. Both the simulations and field observations support the thesis that speciation by ecotypification is occurring, or has occurred, prompting consideration of Bamako as an independent species.

Evidence for divergent selection between the molecular forms of Anopheles gambiae: role of predation

BACKGROUND

The molecular forms of Anopheles gambiae are undergoing speciation. They are characterized by a strong assortative mating and they display partial habitat segregation. The M form is mostly found in flooded/irrigated areas whereas the S form dominates in the surrounding areas, but the ecological factors that shape this habitat segregation are not known. Resource competition has been demonstrated between species undergoing divergent selection, but resource competition is not the only factor that can lead to divergence.

RESULTS

In a field experiment using transplantation of first instar larvae, we evaluated the role of larval predators in mediating habitat segregation between the forms. We found a significant difference in the ability of the molecular forms to exploit the different larval sites conditioned on the presence of predators. In absence of predation, the molecular forms outcompeted each other in their respective natural habitats however, the developmental success of the M form was significantly higher than that of the S form in both habitats under predator pressure.

CONCLUSION

Our results provide the first empirical evidence for specific adaptive differences between the molecular forms and stress the role of larval predation as one of the mechanisms contributing to their divergence.

A genomic view of introgression and hybrid speciation

Hybridization in plants and animals is more common and has more complex outcomes than previously realized. Genome-wide analyses of introgression in organisms ranging from oaks to sunflowers to fruit flies show that a substantial fraction of their genomes are permeable to alleles from related species. Hybridization can lead to rapid genomic changes, including chromosomal rearrangements, genome expansion, differential gene expression, and gene silencing, some of which are mediated by transposable elements. These genomic changes may lead to beneficial new phenotypes, and selection for fertility and ecological traits may in turn alter genome structure. Dramatic increases in the availability of genomic tools will produce a new understanding of the genetic nature of species and will resolve a century-old debate over the basis of hybrid vigor, while the natural recombinants found in hybrid zones will permit genetic mapping of species differences and reproductive barriers in nonmodel organisms.