Search

In human‐dominated landscapes, connectivity is crucial for maintaining demographically stable mammalian populations. Here, we provide a comprehensive noninvasive genetic study for the brown bear population in the Hellenic Peninsula. We analyze its population structuring and connectivity, estimate its population size throughout its distribution, and describe its phylogeography in detail for the first time. Our results, based on 150 multilocus genotypes and on 244‐bp sequences of the mtDNA control region, show the population is comprised by three highly differentiated genetic clusters, consistent with geographical populations of Pindos, Peristeri, and Rhodope. By detecting two male bears with Rhodopean ancestry in the western demes, we provide strong evidence for the ongoing genetic connectivity of the geographically fragmented eastern and western distributions, which suggests connectivity of the larger East Balkan and Pindos‐Dinara populations. Total effective population size (Ne) was estimated to be 199 individuals, and total combined population size (NC) was 499, with each cluster showing a relatively high level of genetic variability, suggesting that migration has been sufficient to counteract genetic erosion. The mtNDA results were congruent with the microsatellite data, and the three genetic clusters were matched predominantly with an equal number of mtDNA haplotypes that belong to the brown bear Western mitochondrial lineage (Clade 1), with two haplotypes being globally new and endemic. The detection of a fourth haplotype that belongs to the Eastern lineage (Clade 3a1) in three bears from the western distribution places the southernmost secondary contact zone between the Eastern and Western lineages in Greece and generates new hypotheses about postglacial maxima migration routes. This work indicates that the genetic composition and diversity of Europe's low‐latitude fringe population are the outcome of ancient and historical events and highlight its importance for the connectivity and long‐term persistence of the species in the Balkans.

Hybridization is a common evolutionary process with multiple possible outcomes. In vertebrates, interspecific hybridization has repeatedly generated parthenogenetic hybrid species. However, it is unknown whether the generation of parthenogenetic hybrids is a rare outcome of frequent hybridization between sexual species within a genus or the typical outcome of rare hybridization events. Darevskia is a genus of rock lizards with both hybrid parthenogenetic and sexual species. Using capture sequencing, we estimate phylogenetic relationships and gene flow among the sexual species, to determine how introgressive hybridization relates to the origins of parthenogenetic hybrids. We find evidence for widespread hybridization with gene flow, both between recently diverged species and deep branches. Surprisingly, we find no signal of gene flow between parental species of the parthenogenetic hybrids, suggesting that the parental pairs were either reproductively or geographically isolated early in their divergence. The generation of parthenogenetic hybrids in Darevskia is, then, a rare outcome of the total occurrence of hybridization within the genus, but the typical outcome when specific species pairs hybridize. Our results question the conventional view that parthenogenetic lineages are generated by hybridization in a window of divergence. Instead, they suggest that some lineages possess specific properties that underpin successful parthenogenetic reproduction.

Inversions often underlie complex adaptive traits, but the genic targets inside them are largely unknown. Gene expression profiling provides a powerful way to link inversions with their phenotypic consequences. We examined the effects of the Cf-Inv(1) inversion in the seaweed fly Coelopa frigida on gene expression variation across sexes and life stages. Our analyses revealed that Cf-Inv(1) shapes global expression patterns but the extent of this effect is variable with much stronger effects in adults than larvae. Furthermore, within adults, both common as well as sex specific patterns were found. The vast majority of these differentially expressed genes mapped to Cf-Inv(1). However, genes that were differentially expressed in a single context (i.e. in males, females or larvae) were more likely to be located outside of Cf-Inv(1). By combining our findings with genomic scans for environmentally associated SNPs, we were able to pinpoint candidate variants in the inversion that may underlie mechanistic pathways that determine phenotypes. Together the results in this study, combined with previous findings, support the notion that the polymorphic Cf-Inv(1) inversion in this species is a major factor shaping both coding and regulatory variation resulting in highly complex adaptive effects.

Genetic incompatibilities contribute to reproductive isolation between many diverging populations, but it is still unclear to what extent they play a role if divergence happens with gene flow. In contact zones between the "Crab" and "Wave" ecotypes of the snail Littorina saxatilis, divergent selection forms strong barriers to gene flow, while the role of post‐zygotic barriers due to selection against hybrids remains unclear. High embryo abortion rates in this species could indicate the presence of such barriers. Post‐zygotic barriers might include genetic incompatibilities (e.g. Dobzhansky–Muller incompatibilities) but also maladaptation, both expected to be most pronounced in contact zones. In addition, embryo abortion might reflect physiological stress on females and embryos independent of any genetic stress. We examined all embryos of >500 females sampled outside and inside contact zones of three populations in Sweden. Females' clutch size ranged from 0 to 1,011 embryos (mean 130 ± 123), and abortion rates varied between 0% and 100% (mean 12%). We described female genotypes by using a hybrid index based on hundreds of SNPs differentiated between ecotypes with which we characterized female genotypes. We also calculated female SNP heterozygosity and inversion karyotype. Clutch size did not vary with female hybrid index, and abortion rates were only weakly related to hybrid index in two sites but not at all in a third site. No additional variation in abortion rate was explained by female SNP heterozygosity, but increased female inversion heterozygosity added slightly to increased abortion. Our results show only weak and probably biologically insignificant post‐zygotic barriers contributing to ecotype divergence, and the high and variable abortion rates were marginally, if at all, explained by hybrid index of females., High abortion rates (12% in average) among females of the snail Littorina saxatilis might suggest postzygotic barriers to gene flow between ecotypes living in adjacent microhabitats. However, sampling several microhabitat transects did not more than marginally support increased rates of abortion in hybrid zones.

Population genetic models of evolution along linear environmental gradients cannot explain why adaptation stops at ecological margins. This is because, unless models impose reductions in carrying capacity at species’ edges, the dominant effect of gene flow is to increase genetic variance and adaptive potential rather than swamping local adaptation. This allows the population to match even very steep changes in trait optima. We extend our previous simulations to explore two nonlinear models of ecological gradients: (a) a sigmoid (steepening) gradient and (b) a linear gradient with a flat centre of variable width. We compare the parameter conditions that allow local adaptation and range expansion from the centre, with those that permit the persistence of a perfectly adapted population distributed across the entire range. Along nonlinear gradients, colonization is easier, and extinction rarer, than along a linear gradient. This is because the shallow environmental gradient near the range centre does not cause gene flow to increase genetic variation, and so does not result in reduced population density. However, as gradient steepness increases, gene flow inflates genetic variance and reduces local population density sufficiently that genetic drift overcomes local selection, creating a finite range margin. When a flat centre is superimposed on a linear gradient, gene flow increases genetic variation dramatically at its edges, leading to an abrupt reduction in density that prevents niche expansion. Remarkably local interruptions in a linear ecological gradient (of a width much less than the mean dispersal distance) can prevent local adaptation beyond this flat centre. In contrast to other situations, this effect is stronger and more consistent where carrying capacity is high. Practically speaking, this means that habitat improvement at patch margins will make evolutionary rescue more likely. By contrast, even small improvements in habitat at patch centres may confine populations to limited areas of ecological space.

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.

Both classical and recent studies suggest that chromosomal inversion polymorphisms are important in adaptation and speciation. However, biases in discovery and reporting of inversions make it difficult to assess their prevalence and biological importance. Here, we use an approach based on linkage disequilibrium among markers genotyped for samples collected across a transect between contrasting habitats to detect chromosomal rearrangements de novo. We report 17 polymorphic rearrangements in a single locality for the coastal marine snail, Littorina saxatilis. Patterns of diversity in the field and of recombination in controlled crosses provide strong evidence that at least the majority of these rearrangements are inversions. Most show clinal changes in frequency between habitats, suggestive of divergent selection, but only one appears to be fixed for different arrangements in the two habitats. Consistent with widespread evidence for balancing selection on inversion polymorphisms, we argue that a combination of heterosis and divergent selection can explain the observed patterns and should be considered in other systems spanning environmental gradients. This article is protected by copyright. All rights reserved.

Zones of secondary contact between closely related taxa are a common legacy of the Quaternary ice ages. Despite their abundance, the factors that keep species apart and prevent hybridization are often unknown. Here, we study a very narrow contact zone between three closely related butterfly species of the Erebia tyndarus species complex. Using genomic data, we first determined whether gene flow occurs and then assessed whether it might be hampered by differences in chromosome number between some species. We found interspecific gene flow between sibling species that differ in karyotype by one chromosome. Conversely, only F1 hybrids occurred between two species that have the same karyotype, forming a steep genomic cline. In a second step, we fitted clines to phenotypic, ecological and parasitic data to identify the factors associated with the genetic cline. We found clines for phenotypic data and the prevalence of the endosymbiont parasite Wolbachia to overlap with the genetic cline, suggesting that they might be drivers for separating the two species. Overall, our results highlight that some gene flow is possible between closely related species despite different chromosome numbers, but that other barriers restrict such gene flow.

A fundamental issue in speciation research is to evaluate phenotypic variation and the genomics driving the evolution of reproductive isolation between sister taxa. Above all, hybrid zones are excellent study systems for researchers to examine the association of genetic differentiation, phenotypic variation and the strength of selection. We investigated two contact zones in the marine gastropod Littorina saxatilis and utilized landmark-based geometric morphometric analysis together with amplified fragment length polymorphism (AFLP) markers to assess phenotypic and genomic divergence between ecotypes under divergent selection. From genetic markers, we calculated the cline width, linkage disequilibrium and the average effective selection on a locus. Additionally, we conducted an association analysis linking the outlier loci and phenotypic variation between ecotypes and show that a proportion of outlier loci are associated with key adaptive phenotypic traits.

The role of natural selection in promoting reproductive isolation has received substantial renewed interest within the last two decades. As a consequence, the study of ecological speciation has become an extremely productive research area in modern evolutionary biology. Recent innovations in sequencing technologies offer an unprecedented opportunity to study the mechanisms involved in ecological speciation. Genome scans provide significant insights but have some important limitations; efforts are needed to integrate them with other approaches to make full use of the sequencing data deluge. An international conference 'Advances in Ecological Speciation' organized by the University of Porto (Portugal) aimed to review current progress in ecological speciation. Using some of the examples presented at the conference, we highlight the benefits of integrating ecological and genomic data and discuss different mechanisms of parallel evolution. Finally, future avenues of research are suggested to advance our knowledge concerning the role of natural selection in the establishment of reproductive isolation during ecological speciation. © 2013 John Wiley & Sons Ltd., The AES conference was organized within the framework of a research project (FCOMP-01-0124-FEDER-014272, PTDC/BIA-EVF/113805/2009) financed by FEDER through COMPETE and by Portuguese national funds through FCT (Fundação para a Ciência e a Tecnologia). This meeting was part of the European Science Foundation's Research Networking Programme ‘Frontiers of Speciation Research’ (FroSpects, www.iiasa.ac.at/Research/EEP/FroSpects) and was cosponsored by the REGPOT-CT-2011-286431-CIBIO-NEW-GEN (PI: Nuno Ferrand) within the framework of the program SP4 – Capacities, Coordination and Support Action from the European Commission (EU FP7) and by the University of Porto. RF, CP, JMF and MM are financed by FCT under the Programa Operacional Potencial Humano – Quadro de Referência Estratégico Nacional from the European Social Fund and the Portuguese Ministério da Educação e Ciência through postdoctoral fellowships (SFRH/BPD/26384/2006 and SFRH/BPD/89313/2012, SFRH/BPD/28869/2006, SFRH/BPD/43264/2008 and SFRH/BPD/46407/2008, respectively). JG is currently supported by an ‘Isidro Parga Pondal’ fellowship (Xunta de Galicia). RKB is supported by NERC and by the Swedish Research Council.

Eleven microsatellite markers were used to determine the genetic population structure and spread of Aedes aegypti (Stegomyia aegypti) (Diptera: Culicidae) in Pakistan using mosquitoes collected from 13 different cities. There is a single genetic cluster of Ae. aegypti in Pakistan with a pattern of isolation by distance within the population. The low level of isolation by distance suggests the long-range passive dispersal of this mosquito, which may be facilitated by the tyre trade in Pakistan. A decrease in genetic diversity from south to north suggests a recent spread of this mosquito from Karachi. A strong negative correlation between genetic distance and the quality of road connections shows that populations in cities connected by better road networks are less differentiated, which suggests the human-aided passive dispersal of Ae. aegypti in Pakistan. Dispersal on a large spatial scale may facilitate the strategy of introducing transgenic Ae. aegypti or intracellular bacteria such as Wolbachia to control the spread of dengue disease in Pakistan, but it also emphasizes the need for simple measures to control container breeding sites.

A pattern of greater divergence in mating traits between sister-species pairs with overlapping ranges than between allopatric species pairs is expected if reinforcement commonly contributes to speciation. Few large-scale comparative analyses have addressed this prediction, especially for genital form. Here, we show that penial morphology follows the predicted pattern in 40 robustly identified sister-species pairs in the marine gastropod subfamily Littorininae. Further work is needed to exclude other processes that may contribute to genital divergence between sympatric species, but the clear pattern we observe strongly suggests a role for genital form in reproductive isolation in this large clade.

Understanding the drivers of speciation is critical to interpreting patterns of biodiversity. The identification of the genetic changes underlying adaptation and reproductive isolation is necessary to link barriers to gene flow to the causal origins of divergence. Here, we present a novel approach to the genetics of speciation, which should complement the commonly used approaches of quantitative trait locus mapping and genome-wide scans for selection. We present a large-scale candidate gene approach by means of sequence capture, applied to identifying the genetic changes underlying reproductive isolation in the pea aphid, a model system for the study of ecological speciation. Targeted resequencing enabled us to scale up the candidate gene approach, specifically testing for the role of chemosensory gene families in host plant specialization. Screening for the signature of divergence under selection at 172 candidate and noncandidate loci, we revealed a handful of loci that show high levels of differentiation among host races, which almost all correspond to odorant and gustatory receptor genes. This study offers the first indication that some chemoreceptor genes, often tightly linked together in the genome, could play a key role in local adaptation and reproductive isolation in the pea aphid and potentially other phytophagous insects. Our approach opens a new route toward the functional genomics of ecological speciation.

Sexual selection has the potential to contribute to population divergence and speciation. Most studies of sexual selection in Drosophila have concentrated on a single signaling modality, usually either courtship song or cuticular hydrocarbons (CHCs), which can act as contact pheromones. We have examined the relationship between both signal types and reproductive success using F1–3 offspring of wild-collected flies, raised in the lab. We used two populations of the Holarctic species Drosophila montana that represent different phylogeographic clades that have been separate for ca. 0.5 million years (MY), and differ to some extent in both traits. Here, we characterize the nature and identify the targets of sexual selection on song, CHCs, and both traits combined within the populations. Three measures of courtship outcome were used as fitness proxies. They were the probability of mating, mating latency, and the production of rejection song by females, and showed patterns of association with different traits that included both linear and quadratic selection. Courtship song predicted courtship outcome better than CHCs and the signal modalities acted in an additive rather than synergistic manner. Selection was generally consistent in direction and strength between the two populations and favored males that sang more vigorously. Sexual selection differed in the extent, strength, and nature on some of the traits between populations. However, the differences in the directionality of selection detected were not a good predictor of population differences. In addition, a character previously shown to be important for species recognition, interpulse interval, was found to be under sexual selection. Our results highlight the complexity of understanding the relationship between within-population sexual selection and population differences. Sexual selection alone cannot predict differences between populations.

Poleward range expansions are widespread responses to recent climate change and are crucial for the future persistence of many species. However, evolutionary change in traits such as colonization history and habitat preference may also be necessary to track environmental change across a fragmented landscape. Understanding the likelihood and speed of such adaptive change is important in determining the rate of species extinction with ongoing climate change. We conducted an amplified fragment length polymorphism (AFLP)-based genome scan across the recently expanded UK range of the Brown Argus butterfly, Aricia agestis, and used outlier-based (DFDIST and BayeScan) and association-based (Isolation-By-Adaptation) statistical approaches to identify signatures of evolutionary change associated with range expansion and habitat use. We present evidence for (i) limited effects of range expansion on population genetic structure and (ii) strong signatures of selection at approximately 5% AFLP loci associated with both the poleward range expansion of A. agestis and differences in habitat use across longestablished and recently colonized sites. Patterns of allele frequency variation at these candidate loci suggest that adaptation to new habitats at the range margin has involved selection on genetic variation in habitat use found across the long-established part of the range. Our results suggest that evolutionary change is likely to affect species’ responses to climate change and that genetic variation in ecological traits across species’ distributions should be maximized to facilitate range shifts across a fragmented landscape, particularly in species that show strong associations with particular habitats.

How common is speciation-with-gene-flow? How much does gene flow impact on speciation? To answer questions like these requires understanding of the common obstacles to evolving reproductive isolation in the face of gene flow and the factors that favour this crucial step. We provide a common framework for the ways in which gene flow opposes speciation and the potential conditions that may ease divergence. This framework is centred on the challenge shared by most scenarios of speciation-with-gene-flow, i.e. the need for coupling among different components of reproductive isolation. Using this structure, we review and compare the factors favouring speciation with the intention of providing a more integrated picture of speciation-with-gene-flow.

South-East Asia is one of the world's richest regions in terms of biodiversity. An understanding of the distribution of diversity and the factors shaping it is lacking, yet essential for identifying conservation priorities for the region's highly threatened biodiversity. Here, we take a large-scale comparative approach, combining data from nine forest-associated Anopheles mosquito species and using statistical phylogeographical methods to disentangle the effects of environmental history, species-specific ecology and random coalescent effects. Spatially explicit modelling of Pleistocene demographic history supports a common influence of environmental events in shaping the genetic diversity of all species examined, despite differences in species' mtDNA gene trees. Populations were periodically restricted to allopatric northeastern and northwestern refugia, most likely due to Pleistocene forest fragmentation. Subsequent southwards post-glacial recolonization is supported by a north-south gradient of decreasing genetic diversity. Repeated allopatric fragmentation and recolonization have led to the formation of deeply divergent geographical lineages within four species and a suture zone where these intraspecific lineages meet along the Thai-Myanmar border. A common environmental influence for this divergence was further indicated by strong support for simultaneous divergence within the same four species, dating to approximately 900 thousand years ago (kya). Differences in the geographical structuring of genetic diversity between species are probably the result of varying species' biology. The findings have important implications for conservation planning; if the refugial regions and suture zone identified here are shared by other forest taxa, the unique and high levels of genetic diversity they house will make these areas conservation priorities.

A comprehensive knowledge on the genetic basis of coloration is crucial to understand how new colour phenotypes arise and how they contribute to the emergence of new species. Variation in melanocortin-1 receptor (Mc1r), a gene that has been reported as a target for repeated evolution in a wide range of vertebrate taxa, was assessed in European ocellated lizards (Lacerta lepida) to search for associations with melanin-based colour phenotypes. Lacerta lepida subspecies’ distribution is associated with the three major bio-climatic regions in the Iberian Peninsula. A nonconserved and derived substitution (T162I) was associated with the L. l. nevadensis phenotype (prevalence of brown scales). Another substitution (S172C) was associated with the presence of black scales in both L. l. lepida and L. l. iberica, but no mutations were found to be associated with the higher proportion of black in L. l. iberica. Extensive genotyping of Mc1r along the contact zone between L. l. nevadensis and L. l. lepida revealed low gene flow (only two hybrids detected). The implications of these findings are discussed in the context of previous knowledge about the evolutionary history of ocellated lizards.

We present an interactive, searchable expressed sequence tag database for the periwinkle snail Littorina saxatilis, an upcoming model species in evolutionary biology. The database is the result of a hybrid assembly between Sanger and 454 sequences, 1290 and 147 491 sequences respectively. Normalized and non-normalized cDNA was obtained from different ecotypes of L. saxatilis collected in the UK and Sweden. The Littorina sequence database (LSD) contains 26 537 different contigs, of which 2453 showed similarity with annotated proteins in UniProt. Querying the LSD permits the selection of the taxonomic origin of blast hits for each contig, and the search can be restricted to particular taxonomic groups. The database allows access to UniProt annotations, blast output, protein family domains (PFAM) and Gene Ontology. The database will allow users to search for genetic markers and identifying candidate genes or genes for expression analyses. It is open for additional deposition of sequence information for L. saxatilis and other species of the genus Littorina. The LSD is available at http://mbio-serv2.mbioekol.lu.se/Littorina/. (Less)

Identification of loci with adaptive importance is a key step to understand the speciation process in natural populations, because those loci are responsible for phenotypic variation that affects fitness in different environments. We conducted an AFLP genome scan in populations of ocellated lizards (Lacerta lepida) to search for candidate loci influenced by selection along an environmental gradient in the Iberian Peninsula. This gradient is strongly influenced by climatic variables, and two subspecies can be recognized at the opposite extremes: L. lepida iberica in the northwest and L. lepida nevadensis in the southeast. Both subspecies show substantial morphological differences that may be involved in their local adaptation to the climatic extremes. To investigate how the use of a particular outlier detection method can influence the results, a frequentist method, DFDIST, and a Bayesian method, BayeScan, were used to search for outliers influenced by selection. Additionally, the spatial analysis method was used to test for associations of AFLP marker band frequencies with 54 climatic variables by logistic regression. Results obtained with each method highlight differences in their sensitivity. DFDIST and BayeScan detected a similar proportion of outliers (3–4%), but only a few loci were simultaneously detected by both methods. Several loci detected as outliers were also associated with temperature, insolation or precipitation according to spatial analysis method. These results are in accordance with reported data in the literature about morphological and life-history variation of L. lepida subspecies along the environmental gradient.

The application of DNA barcoding to dietary studies allows prey taxa to be identified in the absence of morphological evidence and permits a greater resolution of prey identity than is possible through direct examination of faecal material. For insectivorous bats, which typically eat a great diversity of prey and which chew and digest their prey thoroughly, DNA-based approaches to diet analysis may provide the only means of assessing the range and diversity of prey within faeces. Here, we investigated the effectiveness of DNA barcoding in determining the diets of bat species that specialize in eating different taxa of arthropod prey. We designed and tested a novel taxon-specific primer set and examined the performance of short barcode sequences in resolving prey species. We recovered prey DNA from all faecal samples and subsequent cloning and sequencing of PCR products, followed by a comparison of sequences to a reference database, provided species-level identifications for 149/207 (72%) clones. We detected a phylogenetically broad range of prey while completely avoiding detection of nontarget groups. In total, 37 unique prey taxa were identified from 15 faecal samples. A comparison of DNA data with parallel morphological analyses revealed a close correlation between the two methods. However, the sensitivity and taxonomic resolution of the DNA method were far superior. The methodology developed here provides new opportunities for the study of bat diets and will be of great benefit to the conservation of these ecologically important predators.

Genome scans have been used in the studies of ecological speciation to find genomic regions (‘outlier loci’) showing reduced gene flow between divergent populations ⁄ species. High-throughput sequencing (‘454’) offers new opportunities in this field via transcriptome sequencing. Divergent ecotypes of the marine gastropod Littorina saxatilis represent a good example of incipient ecological speciation. We performed a 454-based genome scan between H and M ecotypes of L. saxatilis from the British Isles using cDNA of pooled individuals. Allele frequencies were calculated for 2454 single nucleotide polymorphisms (SNPs), within 572 contigs, and 7% of loci were detected as outliers. Functional annotation of the contigs containing outlier SNPs showed that they included shell matrix and muscle proteins (lithostathine, mucin, titin), proteins involved in energetic metabolism (arginine kinase, NADH dehydrogenase) and reverse transcriptases. Follow-up investigations into these proteins and unannotated outliers will be a promising route in the study of ecological speciation in L. saxatilis.

All species are restricted in their distribution. Currently, ecological models can only explain such limits if patches vary in quality, leading to asymmetrical dispersal, or if genetic variation is too low at the margins for adaptation. However, population genetic models suggest that the increase in genetic variance resulting from dispersal should allow adaptation to almost any ecological gradient. Clearly therefore, these models miss something that prevents evolution in natural populations. We developed an individual-based simulation to explore stochastic effects in these models. At high carrying capacities, our simulations largely agree with deterministic predictions. However, when carrying capacity is low, the population fails to establish for a wide range of parameter values where adaptation was expected from previous models. Stochastic or transient effects appear critical around the boundaries in parameter space between simulation behaviours. Dispersal, gradient steepness, and population density emerge as key factors determining adaptation on an ecological gradient.

Transitions from sexual to asexual reproduction are often coupled with elevations in ploidy. As a consequence, the importance of ploidy per se for the maintenance and spread of asexual populations is unclear. To examine the effects of ploidy and asexual reproduction as independent determinants of the success of asexual lineages, we sampled diploid sexual, diploid asexual, and triploid asexual Eucypris virens ostracods across a European wide range. Applying nuclear and mitochondrial markers, we found that E. virens consists of genetically highly differentiated diploid sexual populations, to the extent that these sexual clades could be considered as cryptic species. All sexual populations were found in southern Europe and North Africa and we found that both diploid asexual and triploid asexual lineages have originated multiple times from several sexual lineages. Therefore, the asexual lineages show a wide variety of genetic backgrounds and very strong population genetic structure across the wide geographic range. Finally, we found that triploid, but not diploid, asexual clones dominate habitats in northern Europe. The limited distribution of diploid asexual lineages, despite their shared ancestry with triploid asexual lineages, strongly suggests that the wider geographic distribution of triploids is due to elevated ploidy rather than to asexuality.

The two species of the palm genus Howea (Arecaceae) from Lord Howe Island, a minute volcanic island in the Tasman Sea, are now regarded as one of the most compelling examples of sympatric speciation, although this view is still disputed by some authors. Population genetic and ecological data are necessary to provide a more coherent and comprehensive understanding of this emerging model system. Here, we analyse data on abundance, juvenile recruitment, pollination mode and genetic variation and structure in both species. We find that Howea forsteriana is less abundant than Howea belmoreana. The genetic data based on amplified fragment length polymorphisms markers indicate similar levels of variation in the two species, despite the estimated census population size of H. belmoreana being three times larger than that of H. forsteriana. Genetic structure within species is low although some weak isolation by distance is detectable. Gene flow between species appears to be extremely limited and restricted to early-generation hybrids - only three admixed individuals, classified as F2s or first generation backcrosses to a parental species, were found among sampled palms. We conclude that speciation in Howea was indeed sympatric, although under certain strict definitions it may be called parapatric.

International audience; Understanding how speciation can take place in the presence of homogenizing gene flow remains a major challenge in evolutionary biology. In the early stages of ecological speciation, reproductive isolation between populations occupying different habitats is expected to be concentrated around genes for local adaptation. These genomic regions will show high divergence while gene exchange in other regions of the genome should continue relatively unimpaired, resulting in low levels of differentiation. The problem is to explain how speciation progresses from this point towards complete reproductive isolation, allowing genome-wide divergence. A new study by Via and West (2008) on speciation between host races of the pea aphid, Acyrthosiphon pisum, introduces the mechanism of 'divergence hitchhiking' which can generate large 'islands of differentiation' and facilitate the build-up of linkage disequilibrium, favouring increased reproductive isolation. This idea potentially removes a major stumbling block to speciation under continuous gene flow.

Amplified fragment length polymorphism (AFLP) fingerprint data are now commonly collected using DNA sequencers. AFLP genotypes are still often scored by eye from such data - a time-consuming, error-prone and subjective process. We present a semi-automated method of genotyping sequencer-collected AFLPs at predefined fragment locations (loci) within the fingerprint. Our method uses thresholds of AFLP-polymerase chain reaction-product fluorescence intensity (peak height) in order to: (i) exclude AFLP loci that are likely to contribute high rates of error to data sets, and (ii) determine the AFLP phenotype (fragment presence or absence) at the retained loci. Error rate analysis is an integral part of this process and is used to determine optimal thresholds that minimize genotyping error, while maximizing the numbers of retained loci. We show that application of this method to a large AFLP data set allows genotype calls that are rapid, objective and repeatable, facilitating the extraction of reliable genotype data for molecular ecological studies.

Chemoperception plays a key role in adaptation and speciation in animals, and the senses of olfaction and gustation are mediated by gene families which show large variation in repertoire size among species. In Drosophila, there are around 60 loci of each type and it is thought that ecological specialization influences repertoire size, with increased pseudogenization of loci. Here, we analyse the size of the gustatory and olfactory repertoires among the genomes of 12 species of Drosophila. We find that repertoire size varies substantially and the loci are evolving by duplication and pseudogenization, with striking examples of lineage-specific duplication. Selection analyses imply that the majority of loci are subject to purifying selection, but this is less strong in gustatory loci and in loci prone to duplication. In contrast to some other studies, we find that few loci show statistically significant evidence of positive selection. Overall genome size is strongly correlated with the proportion of duplicated chemoreceptor loci, but genome size, specialization and endemism may be interrelated in their influence on repertoire size.

The Eurasian badger (Meles meles) has been implicated in the transmission of bovine tuberculosis (TB, caused by Mycobacterium bovis) to cattle. However, evidence suggests that attempts to reduce the spread of TB among cattle in Britain by culling badgers have mixed effects. A large-scale field experiment (the randomized badger culling trial, RBCT) showed that widespread proactive badger culling reduced the incidence of TB in cattle within culled areas but that TB incidence increased in adjoining areas. Additionally, localized reactive badger culling increased the incidence of TB in cattle. It has been suggested that culling-induced perturbation of badger social structure may increase individual movements and elevate the risk of disease transmission between badgers and cattle. Field studies support this hypothesis, by demonstrating increases in badger group ranges and the prevalence of TB infection in badgers following culling. However, more evidence on the effect of culling on badger movements is needed in order to predict the epidemiological consequences of this control strategy. Here, analysis of the genetic signatures of badger populations in the RBCT revealed increased dispersal following culling. While standard tests provided evidence for greater dispersal after culling, a novel method indicated that this was due to medium- and long-distance dispersal, in addition to previously reported increases in home-range size. Our results also indicated that, on average, badgers infected with M. bovis moved significantly farther than did uninfected badgers. A disease control strategy that included culling would need to take account of the potentially negative epidemiological consequences of increased badger dispersal.

By contrast with a multitude of laboratory studies on the social organization of fish, relatively little is know about the size, composition and dynamics of free-ranging fish shoals. We give an overview of the available information on fish shoals and assess to what degree the predictions made from laboratory studies are consistent with field data. The section on shoal choice behaviour in the laboratory is structured so that the evidence for different shoaling preferences is discussed in the context of their mechanisms and functions. Predictions based on experiments in captivity regarding preferences for conspecifics, individuals of similar body length and unparasitized fish were highly consistent with field observations on free-ranging shoals whereas preferences for familiar conspecifics and kin remain to be conclusively demonstrated in the field. In general, there is a shortage of studies in which shoaling preferences have been investigated both in the laboratory and the field, and field studies have so far been largely descriptive revealing little about the underlying mechanisms of observed patterns. Given the great importance of fish shoals both in fundamental and applied research, an advancement of our knowledge of their social organization should significantly contribute to a better understanding of a whole range of topics including reciprocal altruism, group-living and self-organization.

The Drosophila virilis species group offers valuable opportunities for studying the roles of chromosomal re-arrangements and mating signals in speciation. The 13 species are divided into two subgroups, the montana and virilis 'phylads'. There is greater differentiation among species within the montana phylad in both karyotype and acoustic signals than exists among members of the virilis phylad. Drosophila montana is a divergent species which is included in the montana phylad. Here, we analyse the phylogeography of D. montana to provide a framework for understanding divergence of acoustic signals among populations. We analysed mitochondrial sequences corresponding to the cytochrome oxidase I and cytochrome oxidase II genes, as well as 16 microsatellite loci, from 108 lines of D. montana covering most of the species' range. The species shows a clear genetic differentiation between North American and Scandinavian populations. Microsatellite allele frequencies and mitochondrial DNA haplotypes gave significant FST values between populations from Canada, USA and Finland. A Bayesian analysis of population structure based on the microsatellite frequencies showed four genetically distinct groups, corresponding to these three populations plus a small sample from Japan. A network based on mitochondrial haplotypes showed two Finnish clades of very different shape and variability, and another clade with all sequences from North America and Japan. All D. montana populations showed evidence of demographic expansion but the patterns inferred by coalescent analysis differed between populations. The divergence times between Scandinavian and North American clades were estimated to range from 450,000 to 900,000 years with populations in Canada and the USA possibly representing descendants of different refugial populations. Long-term separation of D. montana populations could have provided the opportunity for differentiation observed in male signal traits, especially carrier frequency of the song, but relaxation of sexual selection during population expansion may have been necessary.

Hybridization between roach Rutilus rutilus, bream Abramis brama and rudd Scardinius erythrophthaimus was investigated using morphological and genetic analysis employing both nuclear (ITS1) and mitochondrial (cytochrome b) markers. Allele-specific amplification (ASA) reactions for both markers were developed and ITS1 sequence data for all three species are presented. Sequencing detected two ITS1 haplotypes within both roach and bream which most likely evolved in isolation and were subsequently brought together as a result of restocking by anglers. Analysis of cloned hybrid nuclear ribosomal DNA revealed evidence of recombination between parental ITS1 sequences. ASA proved to be an effective method for identifying hybrids and detected multiple ITS1 copies in fishes identified as purebred by morphological analysis. In addition this suggests post-Fl hybridization and introgression may be occurring between roach and bream, and rudd and bream, although some barriers appear to be suppressing backcrosses within the hybrid population. Analysis of the hybrid population demonstrated that hybridization has occurred in both directions.

During autumn ‘swarming’, large numbers of temperate bats chase each other in and around underground sites. Swarming has been proposed to be a mating event, allowing interbreeding between bats from otherwise isolated summer colonies. We studied the population structure of the Natterer's bat (Myotis nattereri), a swarming species in northern England, by sampling bats at seven sites in two swarming areas and at 11 summer colonies. Analysis of molecular variance (amova) and genetic assignment analyses showed that the swarming areas (60 km apart) support significantly different populations. A negative correlation was found between the distance of a summer colony from a swarming area and the assignment of bats to that area. High gene diversity was found at all sites (HE = 0.79) suggesting high gene flow. This was supported by a low FST (0.017) among summer colonies and the absence of isolation by distance or substructure among colonies which visit one swarming area. The FST, although low, was significantly different from zero, which could be explained by a combination of female philopatry and male-mediated gene flow through mating at swarming sites with bats from other colonies. Modelling suggested that if effective size of the summer colonies (Ne) was low to moderate (10–30), all mating must occur at the swarming sites to account for the observed FST. If the Ne was higher (50), in addition to random mating during swarming, there may be nonrandom mating at swarming sites or some within-colony mating. Conservation of swarming sites that support potentially large populations is discussed.

Knowledge of the genetic basis of divergence in mating signal characters that contribute to reproductive isolation is critical to understanding speciation. Here, we describe a semi-automated system for characterizing grasshopper acoustic signals. We used this system to study the genetic basis of divergence in three male calling song components [echeme (EL), syllable (SL) and phrase (PL) lengths] between Chorthippus brunneus and C. jacobsi, two species of grasshoppers that hybridize in northern Spain. We also studied the number of pegs in the stridulatory file. For all characters, additive effects accounted for most of the genetic differentiation between species. However, the three song components also showed small but significant epistatic effects. No sex linkage was detected. Wright–Castle–Lande estimates of the minimum numbers of genetic factors underlying song and peg number divergence were low: peg number (ne ¼ 5.87 ± 5.84), SL (ne ¼ 2.37 ± 4.79) and PL (ne ¼ 0.87 ± 0.86). On the other hand, EL appeared to be controlled by many genes. These results suggest that divergence in SL and PL might be driven by sexual selection whereas EL might not be under selection. This is consistent with experimental results on female song preference in related species. However, the fact that few factors appear to underlie the differences in peg number is surprising. Peg number is not closely related to song characteristics. It often varies between closely related grasshopper species and it has been assumed to be a neutral character. The biometrical approaches used here tend to underestimate the number of factors influencing a trait but provide valuable background for subsequent quantitative trait loci analyses.

Genetic differentiation among populations of Anopheles jeyporiensis was examined using 76 mtDNA COII sequences from 16 sites throughout southern China and northern Vietnam. The COII sequences are AT-rich (74.58%) and reveal high levels of diversity with 39 of 685 sites polymorphic and 50 different haplotypes present. Genetic variation is high within populations and significant geographical structure was detected at both population and regional levels. In the larger samples, the distributions of haplotypes suggest recent population expansion.

Many hybrid zones contain a deficit of hybrid genotypes relative to expectations from tension zone models. This is often associated with separation of parental genotypes into distinct habitats (mosaicism), but sometimes parentals can be found co-occuring in the same local population (bimodality). In both cases, prezygotic isolation may play an important role in determining the genotypic composition of the zone. Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae) meet and form a complex hybrid zone in northern Spain. Analysis of stridulatory peg numbers reveals partial spatial and seasonal isolation in a 25 km 2 area of the zone: C. jacobsi phenotypes predominate in June and July and are present in both valley and mountain habitats; C. brunneus phenotypes predominate in August and are restricted to valley habitats, always in sympatry with C. jacobsi. Strong assortative mating was observed in laboratory mating experiments. Spatial, seasonal and behavioural isolation combine to produce strong premating isolation in the study area. These results suggest a role for premating isolation in maintaining both the mosaic structure and bimodality of this hybrid zone.

1. The grasshopper species Chorthippus brunneus and C. jacobsi form a complex mosaic hybrid zone in northern Spain. Two mark–release–recapture studies were carried out near the centre of the zone in order to make direct estimates of lifetime dispersal. 2. A model framework based on a simple random walk in homogeneous habitat was extended to include the estimation of philopatry and flying propensity. Each model was compared with the real data, correcting for spatial and temporal biases in the data sets. 3. All four data sets (males and females at each site) deviated significantly from a random walk. Three of the data sets showed strong philopatry and three had a long dispersal tail, indicating a low propensity to move further than predicted by the random walk model. 4. Neighbourhood size estimates were 76 and 227 for the two sites. These estimates may underestimate effective population size, which could be increased by the long tail to the dispersal function. The random walk model overestimates lifetime dispersal and hence the minimum spatial scale of adaptation. 5. Best estimates of lifetime dispersal distance of 7–33 m per generation were considerably lower than a previous indirect estimate of 1344 m per generation. This discrepancy could be influenced by prezygotic isolation, an inherent by-product of mosaic hybrid zone structure.

The phylogenetic relationships of fifteen Oriental and two Afrotropical taxa of the Myzomyia Series of Anopheles subgenus Cellia and two outgroup species, An. maculatus (Neocellia Series) and An. dirus A (Neomyzomyia Series), were inferred from nucleotide sequences of the entire 685 bp of the mitochondrial cytochrome oxidase subunit II locus (COII) and 372 bp of the third domain (D3) of the 28S rDNA locus, both separately and together. Alignment of the D3 sequences was achieved with the aid of secondary structure comparisons, and the pattern of nucleotide substitution was best explained by the GTR + I + G model for either separate or combined datasets. Maximum likelihood and maximum parsimony analyses robustly identified five monophylies: An. fluviatilis U and T; An. fluviatilis U and T + An. minimus A, C, E and #157 + An. leesoni; An. filipinae + An. mangyanus; An. filipinae + An. mangyanus + An. aconitus; and An. culicifacies A and B. The results confirm the specific status of An. flavirostris, the close relationship of An. leesoni with the Minimus Complex, and the exclusion of An. jeyporiensis, An. culicifacies s.l and An. funestus from the Minimus Group. All of the species classified as members of the Minimus Group on morphological grounds formed a single clade, which comprised two subgroups: the Minimus Subgroup, including An. minimus s.l., An. fluviatilis s.l., An. leesoni and An. flavirostris, and the Aconitus Subgroup, including An. filipinae, An. mangyanus, An. aconitus, An. pampanai and An. varuna. However, these clades are only weakly supported by the present dataset.

In some species, populations with few founding individuals can be resilient to extreme inbreeding. Inbreeding seems to be the norm in the common bed bug, Cimex lectularius, a flightless insect that, nevertheless, can reach large deme sizes and persist successfully. However, bed bugs can also be dispersed passively by humans, exposing inbred populations to gene flow from genetically distant populations. The introduction of genetic variation through this outbreeding could lead to increased fitness (heterosis) or be costly by causing a loss of local adaptation or exposing genetic incompatibility between populations (outbreeding depression). Here, we addressed how inbreeding within demes and outbreeding between distant populations impact fitness over two generations in this re-emerging public health pest. We compared fitness traits of families that were inbred (mimicking reproduction following a founder event) or outbred (mimicking reproduction following a gene flow event). We found that outbreeding led to increased starvation resistance compared to inbred families, but this benefit was lost after two generations of outbreeding. No other fitness benefits of outbreeding were observed in either generation, including no differences in fecundity between the two treatments. Resilience to inbreeding is likely to result from the history of small founder events in the bed bug. Outbreeding benefits may only be detectable under stress and when heterozygosity is maximized without disruption of coadaptation. We discuss the consequences of these results both in terms of inbreeding and outbreeding in populations with genetic and spatial structuring, as well as for the recent resurgence of bed bug populations.

Although there are several well-established hypotheses for the origins of postmating isolation during allopatric divergence, there have been very few attempts to determine their relative importance in nature. We have developed an approach based on knowledge of the differing evolutionary histories of populations within species that allows systematic comparison of the predictions of these hypotheses. In previous work, we have applied this methodology to mating signal variation and premating reproductive isolation between populations of the meadow grasshopper Chorthippus parallelus. Here we review the principles behind our approach and report a study measuring postmating isolation in the same set of populations. The populations have known and differing evolutionary histories and relationships resulting from the colonization of northern Europe following the last glaciation. We use a maximum-likelihood analysis to compare the observed pattern of postmating isolation with the predictions of the hypotheses that isolation primarily evolves either as a result of gradual accumulation of mutations in allopatry, or through processes associated with colonization, such as founder events. We also quantify the extent to which degree of postmating isolation can be predicted by genetic distance. Our results suggest that although there is only a weak correlation between genetic distance and postmating isolation, long periods of allopatry do lead to postmating isolation. In contrast to the pattern of premating isolation described in our previous study, colonization does not seem to be associated with increased postmating isolation.

Mosquitoes of the Anopheles minimus group (Diptera: Culicidae) from nine Provinces of southern China were identified morphologically and by molecular characterization, using single-strand conformation polymorphisms (SSCPs) and sequence data for the D3 region of the 28S ribosomal DNA and the mitochondrial COII locus. Species A and C (sensu Green etal., 1990) of the An. minimus complex were found to be sympatric in Yunnan Province. Species A occurs eastward from Yunnan through southern Guangxi, Hainan, Guangdong and Taiwan Provinces, whereas species C occurs northward to northern Guangxi, Guizhou and Sichuan Provinces. Morphological and molecular evidence (based on specimens from the field and four isofemale lines) shows that ,,In. minimus forms A and B (sensu Yu & Li, 1984) are morphological variants of species A, which is accepted as An. minimus Theobald sensu stricto (type-locality: Pokfulam, Hong Kong). The so-called subspecies x of An. minimus (sensu Baba, 1950) is reinterpreted as ,,In. aconitus D6nitz. The distribution and vector status of members of the ,,In. minimus group are discussed in relation to the historical and current transmission of malaria and filariasis in China. Both/In. minimus A and C have been implicated as widespread vectors of malaria, whereas only species A has been found in Hainan, where/In. minimus s.l. was a vector of Bancroftian filariasis. The

The grasshoppers Chorthippus brunneus and Chorthippus jacobsi are easily distinguished by male calling song and the number of stridulatory pegs on the hind femur, and form a mosaic hybrid zone in northern Spain. In this paper, we fit a two-dimensional cline to variation in male calling song characters, which are of particular interest as they are likely to be involved in mate choice by females. As with variation in peg number, local habitat makes only a small contribution in explaining deviations in mean song score from clinal expectations. However, the fitted width of the cline for song characters is significantly narrower than for peg number, suggesting that mating signals may be associated with reduced hybrid fitness in the field and that recombination rates are sufficient to allow clines for different characters to diverge in width. Despite this, estimates for the overall elevation in linkage disequilibrium at the zone center, based on covariance between peg and song characters, reveal a substantial overrepresentation of parental genotypes at the cline center relative to the expectations of a tension zone of similar width. Examination of covariance at individual sites reveals that this inflated estimate of linkage disequilibrium is caused by several sites where the distribution of phenotypes is effectively bimodal. This substantial variation in linkage disequilibrium at the cline center could result from local variation in the strength of assortative mating or selection against hybrids, or may reflect the long-distance colonization of empty habitat from outside the hybrid zone, which would continually create new contacts between parental genotypes at the cline center. Hybrid zones like this, in which strong linkage disequilibrium occurs in some situations but not in others, are of particular relevance to speciation research and allow investigation of the spread of combinations of alleles through different genetic and ecological backgrounds.

Speciation requires the acquisition of reproductive isolation, and the circumstances under which this could evolve are of great interest. Are new species formed after the acquisition of generalized incompatibility arising between physically separated populations, or may they arise as a result of the action of disruptive selection beginning with the divergence of a rather restricted set of gene loci? Here we apply the technique of amplified fragment length polymorphism (AFLP) analysis to an intertidal snail whose populations display a cline in shell shape across vertical gradients on rocky shores. We compare the FST values for 306 AFLP loci with the distribution of FST estimated from a simulation model using values of mutation and migration derived from the data. We find that about 5% of these loci show greater differentiation than expected, providing evidence of the effects of selection across the cline, either direct or indirect through linkage. This is consistent with expectations from nonallopatric speciation models that propose an initial divergence of a small part of the genome driven by strong disruptive selection while divergence at other loci is prevented by gene flow. However, the pattern could also be the result of differential introgression after secondary contact.

There is evidence for introgression between the malaria vectors Anopheles gambiae and Anopheles arabiensis, which belong to the Anopheles gambiae complex. Here we consider their relationship to another member of the complex, Anopheles bwambae. Species identifications were made using rDNA-PCR and revealed one An. bwambae/gambiae hybrid among 459 wild-caught mosquitoes. A 237-nucleotide region of the mitochondrial ND5 gene was sequenced in sixty-nine individuals. The results suggest that there is an An. bwambae specific group of haplotypes. However, two An. bwambae individuals carried haplotypes typical of An. gambiae and An. arabiensis. This is the first evidence for introgression between An. bwambae and other species of the An. gambiae complex and suggests that introgression may be more widespread, also occurring between other members of the complex.