The last bastion? X chromosome genotyping of Anopheles gambiae species pair males from a hybrid zone reveals complex recombination within the major candidate ‘genomic island of speciation’

Speciation with gene flow may be aided by reduced recombination helping to build linkage between genes involved in the early stages of reproductive isolation. Reduced recombination on chromosome X has been implicated in speciation within the Anopheles gambiae complex, species of which represent the major Afrotropical malaria vectors. The most recently diverged, morphologically indistinguishable, species pair, A. gambiae and Anopheles coluzzii, ubiquitously displays a ‘genomic island of divergence’ spanning over 4 Mb from chromosome X centromere, which represents a particularly promising candidate region for reproductive isolation genes, in addition to containing the diagnostic markers used to distinguish the species. Very low recombination makes the island intractable for experimental recombination studies, but an extreme hybrid zone in Guinea Bissau offers the opportunity for natural investigation of X‐island recombination. SNP analysis of chromosome X hemizygous males revealed: (i) strong divergence in the X‐island despite a lack of autosomal divergence; (ii) individuals with multiple‐recombinant genotypes, including likely double crossovers and localized gene conversion; (iii) recombination‐driven discontinuity both within and between the molecular species markers, suggesting that the utility of the diagnostics is undermined under high hybridization. The largely, but incompletely protected nature of the X centromeric genomic island is consistent with a primary candidate area for accumulation of adaptive variants driving speciation with gene flow, while permitting some selective shuffling and removal of genetic variation.     

Insecticide susceptibility of Anopheles coluzzii and Anopheles gambiae mosquitoes in Ibadan,Southwest Nigeria

The emergence of insecticide resistance in Anopheles (Diptera: Culicidae) mosquitoes has great implications for malaria control in Nigeria. This study aimed to determine the dynamics of insecticide susceptibility levels and the frequency of knock‐down resistance (kdr) mutations (L1014F) in wild Anopheles coluzzii Coetzee & Wilkerson sp. n. and Anopheles gambiae Giles from the Ojoo and Bodija areas of Ibadan, in southwest Nigeria. Insecticide susceptibility to pyrethroids, organophosphates, carbamates and organochlorines was assessed using World Health Organization (WHO) bioassays. A subset of the mosquitoes exposed to pyrethroids and DDT was used for species and molecular form identification; kdr genotyping was determined using the TaqMan real‐time polymerase chain reaction assay. The mosquitoes were resistant to pyrethroids and DDT but completely susceptible to organophosphates and carbamates. Bodija samples (n = 186) consisted of An. gambiae (91.4%) and An. coluzzii (8.1%) and included one An. coluzzii/An. gambiae hybrid specimen. All mosquitoes screened in Ojoo (n = 26) were An. gambiae. The 1014F kdr mutation was detected at frequencies of 24.5 and 5.8% in Bodija and Ojoo, respectively. No correlation was observed between kdr genotypes and resistance phenotypes. The results indicate that metabolic resistance probably plays an important role in the development of resistance and highlight the need to implement insecticide resistance management strategies.     

Chromosomal inversions and ecotypic differentiation in Anopheles gambiae: the perspective from whole‐genome sequencing

The molecular mechanisms and genetic architecture that facilitate adaptive radiation of lineages remain elusive. Polymorphic chromosomal inversions, due to their recombination‐reducing effect, are proposed instruments of ecotypic differentiation. Here, we study an ecologically diversifying lineage of Anopheles gambiae, known as the Bamako chromosomal form based on its unique complement of three chromosomal inversions, to explore the impact of these inversions on ecotypic differentiation. We used pooled and individual genome sequencing of Bamako, typical (non‐Bamako) An. gambiae and the sister species Anopheles coluzzii to investigate evolutionary relationships and genomewide patterns of nucleotide diversity and differentiation among lineages. Despite extensive shared polymorphism and limited differentiation from the other taxa, Bamako clusters apart from the other taxa, and forms a maximally supported clade in neighbour‐joining trees based on whole‐genome data (including inversions) or solely on collinear regions. Nevertheless, F_ST outlier analysis reveals that the majority of differentiated regions between Bamako and typical An. gambiae are located inside chromosomal inversions, consistent with their role in the ecological isolation of Bamako. Exceptionally differentiated genomic regions were enriched for genes implicated in nervous system development and signalling. Candidate genes associated with a selective sweep unique to Bamako contain substitutions not observed in sympatric samples of the other taxa, and several insecticide resistance gene alleles shared between Bamako and other taxa segregate at sharply different frequencies in these samples. Bamako represents a useful window into the initial stages of ecological and genomic differentiation from sympatric populations in this important group of malaria vectors.     

Evolution of novel mimicry rings facilitated by adaptive introgression in tropical butterflies

Understanding the genetic basis of phenotypic variation and the mechanisms involved in the evolution of adaptive novelty, especially in adaptive radiations, is a major goal in evolutionary biology. Here, we used whole‐genome sequence data to investigate the origin of the yellow hindwing bar in the Heliconius cydno radiation. We found modular variation associated with hindwing phenotype in two narrow noncoding regions upstream and downstream of the cortex gene, which was recently identified as a pigmentation pattern controller in multiple species of Heliconius. Genetic variation at each of these modules suggests an independent control of the dorsal and ventral hindwing patterning, with the upstream module associated with the ventral phenotype and the downstream module with the dorsal one. Furthermore, we detected introgression between H. cydno and its closely related species Heliconius melpomene in these modules, likely allowing both species to participate in novel mimicry rings. In sum, our findings support the role of regulatory modularity coupled with adaptive introgression as an elegant mechanism by which novel phenotypic combinations can evolve and fuel an adaptive radiation.     

A longitudinal study on Anopheles mosquito larval abundance in distinct geographical and environmental settings in western Kenya

Background

Malaria vector control in Africa depends upon effective insecticides in bed nets and indoor residual sprays. This study investigated the extent of insecticide resistance in Anopheles gambiae s.l., Anopheles gambiae s.s. and Anopheles arabiensis in western Kenya where ownership of insecticide-treated bed nets has risen steadily from the late 1990s to 2010. Temporal and spatial variation in the frequency of a knock down resistance (kdr) allele in A. gambiae s.s. was quantified, as was variation in phenotypic resistance among geographic populations of A. gambiae s.l.

Methods

To investigate temporal variation in kdr frequency, individual specimens of A. gambiae s.s. from two sentinel sites were genotyped using RT-PCR from 1996-2010. Spatial variation in kdr frequency, species composition, and resistance status were investigated in additional populations of A. gambiae s.l. sampled in western Kenya in 2009 and 2010. Specimens were genotyped for kdr as above and identified to species via conventional PCR. Field-collected larvae were reared to adulthood and tested for insecticide resistance using WHO bioassays.

Results

Anopheles gambiae s.s. showed a dramatic increase in kdr frequency from 1996 - 2010, coincident with the scale up of insecticide-treated nets. By 2009-2010, the kdr L1014S allele was nearly fixed in the A. gambiae s.s. population, but was absent in A. arabiensis. Near Lake Victoria, A. arabiensis was dominant in samples, while at sites north of the lake A. gambiae s.s was more common but declined relative to A. arabiensis from 2009 to 2010. Bioassays demonstrated that A. gambiae s.s. had moderate phenotypic levels of resistance to DDT, permethrin and deltamethrin while A. arabiensis was susceptible to all insecticides tested.

Conclusions

The kdr L1014S allele has approached fixation in A. gambiae s.s. populations of western Kenya, and these same populations exhibit varying degrees of phenotypic resistance to DDT and pyrethroid insecticides. The near absence of A. gambiae s.s. from populations along the lakeshore and the apparent decline in other populations suggest that insecticide-treated nets remain effective against this mosquito despite the increase in kdr allele frequency. The persistence of A. arabiensis, despite little or no detectable insecticide resistance, is likely due to behavioural traits such as outdoor feeding and/or feeding on non-human hosts by which this species avoids interaction with insecticide-treated nets.     

Limited gene exchange between two sister species of leaf beetles within a hybrid zone in the Alps

Analysing genomic variation within and between sister species is a first step towards understanding species boundaries. We focused on two sister species of cold‐resistant leaf beetles, Gonioctena quinquepunctata and G. intermedia, whose ranges overlap in the Alps. A previous study of DNA sequence variation had revealed multiple instances of mitochondrial genome introgression in this region, suggesting recent hybridization between the two species. To evaluate the extent of gene exchange resulting from these hybridization events, we sampled individuals of both species inside and outside the hybrid zone and analysed genomic variation among them using RAD‐seq markers. Individual levels of introgression in the nuclear genome were estimated first by defining species‐specific SNPs (displaying a fixed difference between species) a priori and second by using model‐based methods. Both types of analyses indicated little gene exchange, if any, between species at the level of the nuclear genome. Whereas the first method suggested slightly more gene flow, we argue that it has likely overestimated introgression in the phylogeographic context of this study. We conclude that strong intrinsic barriers prevent genetic exchange at the level of the nuclear genome between the two species. The apparent discrepancy observed between introgression occurring in the nuclear and mitochondrial genomes could be explained by selection acting in favour of the latter. Also, these results have consequences for the phylogeographic study of each species, since we can assume that genetic diversity in the overlapping portion of their ranges is not the product of introgression.     

Hybridization between two parapatric ranid frog species in the northern Sierra Nevada,California, USA

Contact zones between species provide a unique opportunity to test whether taxa can hybridize or not. Cross‐breeding or hybridization between closely related taxa can promote gene flow (introgression) between species, adaptation, or even speciation. Though hybridization events may be short‐lived and difficult to detect in the field, genetic data can provide information about the level of introgression between closely related taxa. Hybridization can promote introgression between species, which may be an important evolutionary mechanism for either homogenization (reversing initial divergence between species) or reproductive isolation (potentially leading to speciation). Here, we used thousands of genetic markers from nuclear DNA to detect hybridization between two parapatric frog species (Rana boylii and Rana sierrae) in the Sierra Nevada of California. Based on principal components analysis, admixture, and analysis of heterozygosity at species diagnostic SNPs, we detected two F1 hybrid individuals in the Feather River basin, as well as a weak signal of introgression and gene flow between the frog species compared with frog populations from two other adjacent watersheds. This study provides the first documentation of hybridization and introgression between these two species, which are of conservation concern.     

Indirect evidence that agricultural pesticides select for insecticide resistance in the malaria vector Anopheles gambiae

We investigated the possible relationship between the agricultural use of insecticides and the emergence of insecticide resistance. Bioassays were conducted using simulated mosquito larval habitats and well known Anopheles gambiae strains. Soil samples were collected from vegetable production areas in Benin, including one site with insecticide use, one site where insecticides had not been used for two months, and a third where insecticides had not been used. Pupation and emergence rates were very low in pyrethroid‐susceptible strains when exposed to soil that had been recently exposed to insecticides. Pupation and emergence rates in strains with the kdr mutation alone or both the kdr and Ace‐1 mutations were much higher. Overall, strains with the kdr mutation survived at higher rates compared to that without kdr mutation. Although this study is observational, we provide indirect evidence indicating that soils from agricultural areas contain insecticide residues that can play a role in the emergence of insecticide resistance in Anopheles. This aspect should be taken into account to better utilize the insecticide in the context of integrated pest management programs.     

Nuclear introgression without mitochondrial introgression in two turtle species exhibiting sex‐specific trophic differentiation

Despite the presence of reproductive barriers between species, interspecific gene introgression has been documented in a range of natural systems. Comparing patterns of genetic introgression in biparental versus matrilineal markers can potentially reveal sex‐specific barriers to interspecific gene flow. Hybridization has been documented in the freshwater turtles Graptemys geographica and G. pseudogeographica, whose ranges are largely sympatric. Morphological differentiation between the species is restricted to females, with female G. geographica possessing large heads and jaws compared to the narrow heads of G. pseudogeographica females. If hybrid females are morphologically intermediate, they may be less successful at exploiting parental feeding niches, thereby limiting the introgression of maternally inherited, but not biparental, molecular markers. We paired sequence data with stable isotope analysis and examined sex‐specific genetic introgression and trophic differentiation in sympatric populations of G. geographica and G. pseudogeographica. We observed introgression from G. pseudogeographica into G. geographica at three nuclear loci, but not at the mitochondrial locus. Analysis of ?¹⁵N and ?¹³C was consistent with species differences in trophic positioning in females, but not males. These results suggest that ecological divergence in females may reduce the opportunity for gene flow in this system.     

Hybridization and mitochondrial genome introgression between Rana chensinensis and R. kukunoris

Mitochondrial genome (mito‐genome) introgression among metazoans is commonplace, and several biological processes may promote such introgression. We examined two proposed processes for the mito‐genome introgression between Rana chensinensis and R. kukunoris: natural hybridization and sex‐biased dispersal. We sampled 477 individuals from 28 sites in the potential hybrid zone in the western Tsinling Mountains. Mitochondrial gene (cyt‐b) trees were used to examine the introgression events. Microsatellite DNA loci, cyt‐b and morphological data were used to identify hybrids and to examine the extent of natural hybridization. We detected rampant bidirectional introgressions, both ancient and recent, between the two species. Furthermore, we found a wide hybrid zone, and frequent and asymmetric hybridization. The hybrid zone cline analysis revealed a clear mitochondrial–nuclear discordance; while most nuclear markers displayed similar and steep clines, cyt‐b had a displaced cline centre and a more gradual and wider cline. We also detected strong and asymmetric historical maternal gene flow across the hybrid zone. This widespread hybridization and detected low mito‐nuclear conflicts may, at least partially, explain the high frequency of introgression. Lastly, microsatellite data and population genetic methods were used to assess sex‐biased dispersal. A weak pattern of female‐biased dispersal was detected in both species, suggesting it may not play an important role in the observed introgression. Our data are consistent with the hybridization hypothesis, but support for the sex‐biased dispersal hypothesis is weak. We further suggest that selective advantages of the R. kukunoris‐type mito‐genome in thermal adaptation may also contribute to the introgression between the two species.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.     

Behavioural response to combined insecticide and temperature stress in natural populations of Drosophila melanogaster

Insecticide resistance evolves extremely rapidly, providing an illuminating model for the study of adaptation. With climate change reshaping species distribution, pest and disease vector control needs rethinking to include the effects of environmental variation and insect stress physiology. Here, we assessed how both long‐term adaptation of populations to temperature and immediate temperature variation affect the genetic architecture of DDT insecticide response in Drosophila melanogaster. Mortality assays and behavioural assays based on continuous activity monitoring were used to assess the interaction between DDT and temperature on three field‐derived populations from climate extremes (Raleigh for warm temperate, Tasmania for cold oceanic and Queensland for hot tropical). The Raleigh population showed the highest mortality to DDT, whereas the Queensland population, epicentre for derived alleles of the resistance gene Cyp6g1, showed the lowest. Interaction between insecticide and temperature strongly affected mortality, particularly for the Tasmanian population. Activity profiles analysed using self‐organizing maps show that the insecticide promoted an early response, whereas elevated temperature promoted a later response. These distinctive early or later activity phases revealed similar responses to temperature and DDT dose alone but with more or less genetic variance depending on the population. This change in genetic variance among populations suggests that selection particularly depleted genetic variance for DDT response in the Queensland population. Finally, despite similar (co)variation between traits in benign conditions, the genetic responses across population differed under stressful conditions. This showed how stress‐responsive genetic variation only reveals itself in specific conditions and thereby escapes potential trade‐offs in benign environments.     

Malaria vector composition and insecticide susceptibility status in Guinea Conakry,West Africa

This study provides data on malaria vector species composition and insecticide susceptibility status from three localities in Guinea Conakry. A total of 497 mosquitoes were collected resting indoors and morphologically identified as belonging to the Anopheles gambiae complex. The majority of these were An. gambiae s.s. (99.6%), but a small percentage (0.4%) were identified as Anopheles arabiensis. Thirty‐four Anopheles funestus s.s. were also collected. The molecular S form of An. gambiae s.s. was predominant over the M form in Siguiri (95%) and Boffa (97.4%), whereas at Mt Nimba the M form was more abundant (61.4%) than the S form (38.1%). One hybrid M/S specimen was recorded from Mt Nimba. Siguiri populations showed high levels of resistance to DDT, dieldrin and bendiocarb. Anopheles gambiae from Boffa were largely susceptible to the insecticides tested. At Mt Nimba, resistance to DDT and bendicocarb was detected. Biochemical enzyme analysis showed that an altered acetylcholinesterase is operating in the field at low levels. The frequency of the 1014F kdr allele in the An. gambiae S form was 0.24 at Siguiri and 0.14 at Mt Nimba. A single RR specimen was found in the M form. The heterogeneity in species composition and resistance profiles between sites requires vector control interventions to be tailored to each site based on the data collected from ongoing monitoring and surveillance.     

Whole‐genome sequencing reveals absence of recent gene flow and separate demographic histories for Anopheles punctulatus mosquitoes in Papua New Guinea

Anopheles mosquitoes are the vectors of several human diseases including malaria. In many malaria endemic areas, several species of Anopheles coexist, sometimes in the form of related sibling species that are morphologically indistinguishable. Determining the size and organization of Anopheles populations, and possible ongoing gene flow among them is important for malaria control and, in particular, for monitoring the spread of insecticide resistance alleles. However, these parameters have been difficult to evaluate in most Anopheles species due to the paucity of genetic data available. Here, we assess the extent of contemporary gene flow and historical variations in population size by sequencing and de novo assembling the genomes of wild‐caught mosquitoes from four species of the Anopheles punctulatus group of Papua New Guinea. Our analysis of more than 50 Mb of orthologous DNA sequences revealed no evidence of contemporary gene flow among these mosquitoes. In addition, investigation of the demography of two of the An. punctulatus species revealed distinct population histories. Overall, our analyses suggest that, despite their similarities in morphology, behaviour and ecology, contemporary sympatric populations of An. punctulatus are evolving independently.     

Subspecies delineation amid phenotypic,geographic and genetic discordance in a songbird

Understanding the processes that drive divergence within and among species is a long‐standing goal in evolutionary biology. Traditional approaches to assessing differentiation rely on phenotypes to identify intra‐ and interspecific variation, but many species express subtle morphological gradients in which boundaries among forms are unclear. This intraspecific variation may be driven by differential adaptation to local conditions and may thereby reflect the evolutionary potential within a species. Here, we combine genetic and morphological data to evaluate intraspecific variation within the Nelson's (Ammodramus nelsoni) and salt marsh (Ammodramus caudacutus) sparrow complex, a group with populations that span considerable geographic distributions and a habitat gradient. We evaluated genetic structure among and within five putative subspecies of A. nelsoni and A. caudacutus using a reduced‐representation sequencing approach to generate a panel of 1929 SNPs among 69 individuals. Although we detected morphological differences among some groups, individuals sorted along a continuous phenotypic gradient. In contrast, the genetic data identified three distinct clusters corresponding to populations that inhabit coastal salt marsh, interior freshwater marsh and coastal brackish–water marsh habitats. These patterns support the current species‐level recognition but do not match the subspecies‐level taxonomy within each species—a finding which may have important conservation implications. We identified loci exhibiting patterns of elevated divergence among and within these species, indicating a role for local selective pressures in driving patterns of differentiation across the complex. We conclude that this evidence for adaptive variation among subspecies warrants the consideration of evolutionary potential and genetic novelty when identifying conservation units for this group.     

Transcriptomic differences between euryhaline and stenohaline malaria vector sibling species in response to salinity stress

Evolution of osmoregulatory systems is a key factor in the transition of species between fresh‐ and saltwater habitats. Anopheles coluzzii and Anopheles merus are stenohaline and euryhaline malaria vector mosquitoes belonging to a larger group of sibling species, the Anopheles gambiae complex, which radiated in Africa within the last 2 million years. Comparative ecological genomics of these vector species can provide insight into the mechanisms that permitted the rapid radiation of this species complex into habitats of contrasting salinity. Here, we use RNA‐Seq to investigate gene expression differences between An. coluzzii and An. merus after briefly exposing both young and old larval instars of each species to either saltwater (SW) or freshwater (FW). Our study aims to identify candidate genes and pathways responsible for the greater SW tolerance of An. merus. Our results are congruent with the ability of gene induction to mediate salinity tolerance, with both species showing increasing amounts of differential gene expression between SW and FW as salt concentrations increase. Besides ion transporters such as AgAE2 that may serve as effectors for osmoregulation, we also find mitogen‐activated protein kinases that may serve in a phosphorylation signalling pathway responding to salinity, and report potential cross‐talk between the mosquito immune response and osmoregulation. This study provides a key step towards applying the growing molecular knowledge of these malaria vectors to improve understanding of their ecological tolerances and habitat occupancy.     

Diversification and asymmetrical gene flow across time and space: lineage sorting and hybridization in polytypic barking frogs

Young species complexes that are widespread across ecologically disparate regions offer important insights into the process of speciation because of their relevance to how local adaptation and gene flow influence diversification. We used mitochondrial DNA and up to 28 152 genomewide single nucleotide polymorphisms from polytypic barking frogs (Craugastor augusti complex) to infer phylogenetic relationships and test for the signature of introgressive hybridization among diverging lineages. Our phylogenetic reconstructions suggest (i) a rapid Pliocene–Pleistocene radiation that produced at least nine distinct lineages and (ii) that geographic features of the arid Central Mexican Plateau contributed to two independent northward expansions. Despite clear lineage differentiation (many private alleles and high between‐lineage F_ST scores), D‐statistic tests, which differentiate introgression from ancestral polymorphism, allowed us to identify two putative instances of reticulate gene flow. Partitioned D‐statistics provided evidence that these events occurred in the same direction between clades but at different points in time. After correcting for geographic distance, we found that lineages involved in hybrid gene flow interactions had higher levels of genetic variation than independently evolving lineages. These findings suggest that the nature of hybrid compatibility can be conserved overlong periods of evolutionary time and that hybridization between diverging lineages may contribute to standing levels of genetic variation.     

Islands and Stepping-Stones: Comparative Population Structure of Anopheles gambiae sensu stricto and Anopheles arabiensis in Tanzania and Implications for the Spread of Insecticide Resistance

Population genetic structures of the two major malaria vectors Anopheles gambiae s.s. and An. arabiensis, differ markedly across Sub-Saharan Africa, which could reflect differences in historical demographies or in contemporary gene flow. Elucidation of the degree and cause of population structure is important for predicting the spread of genetic traits such as insecticide resistance genes or artificially engineered genes. Here the population genetics of An. gambiae s.s. and An. arabiensis in the central, eastern and island regions of Tanzania were compared. Microsatellite markers were screened in 33 collections of female An. gambiae s.l., originating from 22 geographical locations, four of which were sampled in two or three years between 2008 and 2010. An. gambiae were sampled from six sites, An. arabiensis from 14 sites, and both species from two sites, with an additional colonised insectary sample of each species. Frequencies of the knock-down resistance (kdr) alleles 1014S and 1014F were also determined. An. gambiae exhibited relatively high genetic differentiation (average pairwise F_ST = 0.131), significant even between nearby samples, but without clear geographical patterning. In contrast, An. arabiensis exhibited limited differentiation (average F_ST = 0.015), but strong isolation-by-distance (Mantel test r = 0.46, p = 0.0008). Most time-series samples of An. arabiensis were homogeneous, suggesting general temporal stability of the genetic structure. An. gambiae populations from Dar es Salaam and Bagamoyo were found to have high frequencies of kdr 1014S (around 70%), with almost 50% homozygote but was at much lower frequency on Unguja Island, with no. An. gambiae population genetic differentiation was consistent with an island model of genetic structuring with highly restricted gene flow, contrary to An. arabiensis which was consistent with a stepping-stone model of extensive, but geographically-restricted gene flow.     

Host‐associated genomic differentiation in congeneric butterflies: now you see it,now you do not

Ecotypic variation among populations may become associated with widespread genomic differentiation, but theory predicts that this should happen only under particular conditions of gene flow, selection and population size. In closely related species, we might expect the strength of host‐associated genomic differentiation (HAD) to be correlated with the degree of phenotypic differentiation in host‐adaptive traits. Using microsatellite and Amplified Fragment Length Polymorphism (AFLP) markers, and controlling for isolation by distance between populations, we sought HAD in two congeneric species of butterflies with different degrees of host plant specialization. Prior work on Euphydryas editha had shown strong interpopulation differentiation in host‐adapted traits, resulting in incipient reproductive isolation among host‐associated ecotypes. We show here that Euphydryas aurinia had much weaker host‐associated phenotypic differentiation. Contrary to our expectations, we detected HAD in Euphydryas aurinia, but not in E. editha. Even within an E. aurinia population that fed on both hosts, we found weak but significant sympatric HAD that persisted in samples taken 9 years apart. The finding of significantly stronger HAD in the system with less phenotypic differentiation may seem paradoxical. Our findings can be explained by multiple factors, ranging from differences in dispersal or effective population size, to spatial variation in genomic or phenotypic traits and to structure induced by past histories of host‐adapted populations. Other infrequently measured factors, such as differences in recombination rates, may also play a role. Our result adds to recent work as a further caution against assumptions of simple relationships between genomic and adaptive phenotypic differentiation.