Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y-chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least-divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi-absence of Y-chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non-neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.

The role of Wolbachia bacteria in reproductive incompatibilities and hybrid zones of Diabrotica beetles and Gryllus crickets

A rickettsial bacterium in the genus Wolbachia is the cause of a unidirectional reproductive incompatibility observed between two major beetle pests of maize, the western corn rootworm, Diabrotica virgifera virgifera, and the Mexican corn rootworm, D. v. zeae. These subspecies are allopatric except for two known regions of sympatry in Texas and Mexico. We demonstrate that populations of D. v. virgifera, with the exception of two populations in southern Arizona, are infected with a strain of Wolbachia. Populations of D. v. zeae are not infected. Treatment of D. v. virgifera with tetracycline eliminated the Wolbachia and removed the reproductive incompatibility. Similar patterns of reproductive incompatibility exist among taxa of the cricket genus Gryllus. Gryllus assimilis, G. integer, G. ovisopis, G. pennsylvanicus, and G. rubens are infected with Wolbachia whereas G. firmus is usually not. Populations of G. rubens and G. ovisopis carry the same Wolbachia strain, which is distinct from that of G. integer. G. pennsylvanicus is infected with two Wolbachia strains, that found in G. rubens and one unique to G. pennsylvanicus. Moreover, a proportion of G. pennsylvanicus individuals harbors both strains. Wolbachia may have influenced speciation in some members of the genus Gryllus by affecting the degree of hybridization between species. Given that Wolbachia infections are relatively common in insects, it is likely that other insect hybrid zones may be influenced by infections with Wolbachia.

Strong selection against hybrids at a hybrid zone in the Ensatina ring species complex and its evolutionary implications

The analysis of interactions between lineages at varying levels of genetic divergence can provide insights into the process of speciation through the accumulation of incompatible mutations. Ring species, and especially the Ensatina eschscholtzii system exemplify this approach. The plethodontid salamanders E. eschscholtzii xanthoptica and E. eschscholtzii platensis hybridize in the central Sierran foothills of California. We compared the genetic structure across two transects (southern and northern Calaveras Co.), one of which was resampled over 20 years, and examined diagnostic molecular markers (eight allozyme loci and mitochondrial DNA) and a diagnostic quantitative trait (color pattern). Key results across all studies were: (1) cline centers for all markers were coincident and the zones were narrow, with width estimates of 730 m to 2000 m; (2) cline centers at the northern Calaveras transect were coincident between 1981 and 2001, demonstrating repeatability over five generations; (3) there were very few if any putative F1s, but a relatively high number of backcrossed individuals in the central portion of transects: and (4) we found substantial linkage disequilibrium in all three studies and strong heterozygote deficit both in northern Calaveras, in 2001, and southern Calaveras. Both linkage disequilibrium and heterozygote deficit showed maximum values near the center of the zones. Using estimates of cline width and dispersal, we infer strong selection against hybrids. This is sufficient to promote accumulation of differences at loci that are neutral or under divergent selection, but would still allow for introgression of adaptive alleles. The evidence for strong but incomplete isolation across this centrally located contact is consistent with theory suggesting a gradual increase in postzygotic incompatibility between allopatric populations subject to divergent selection and reinforces the value of Ensatina as a system for the study of divergence and speciation at multiple stages. © 2005 The Society for the Study of Evolution. All rights reserved.

Genomic architecture of adaptive color pattern divergence and convergence in Heliconius butterflies

Identifying the genetic changes driving adaptive variation in natural populations is key to understanding the origins of biodiversity. The mosaic of mimetic wing patterns in Heliconius butterflies makes an excellent system for exploring adaptive variation using next-generation sequencing. In this study, we use a combination of techniques to annotate the genomic interval modulating red color pattern variation, identify a narrow region responsible for adaptive divergence and convergence in Heliconius wing color patterns, and explore the evolutionary history of these adaptive alleles. We use whole genome resequencing from four hybrid zones between divergent color pattern races of Heliconius erato and two hybrid zones of the co-mimic Heliconius melpomene to examine genetic variation across 2.2 Mb of a partial reference sequence. In the intergenic region near optix, the gene previously shown to be responsible for the complex red pattern variation in Heliconius, population genetic analyses identify a shared 65-kb region of divergence that includes several sites perfectly associated with phenotype within each species. This region likely contains multiple cis-regulatory elements that control discrete expression domains of optix. The parallel signatures of genetic differentiation in H. erato and H. melpomene support a shared genetic architecture between the two distantly related co-mimics; however, phylogenetic analysis suggests mimetic patterns in each species evolved independently. Using a combination of next-generation sequencing analyses, we have refined our understanding of the genetic architecture of wing pattern variation in Heliconius and gained important insights into the evolution of novel adaptive phenotypes in natural populations.

Hybridization and introgression across different ploidy levels in the Neotropical orchids Epidendrum fulgens and E-puniceoluteum (Orchidaceae)

The hypothesis of gene flow between species with large differences in chromosome numbers has rarely been tested in the wild, mainly because species of different ploidy are commonly assumed to be reproductively isolated from each other because of instantaneous and strong postzygotic barriers. In this study, a broad-scale survey of molecular variation was carried out between two orchid species with different ploidy levels: Epidendrum fulgens (2n = 2x = 24 chromosomes) and Epidendrum puniceoluteum (2n = 4x = 52 chromosomes). To test the strength of their reproductive barriers, we investigated the distribution of genetic variation in sympatric and allopatric populations of these two species and conducted crossing experiments. Nuclear and plastid microsatellite loci were used to genotype 463 individuals from eight populations across the geographical range of both species along the Brazilian coastal plain. All six sympatric populations analysed presented hybrid zones, indicating that hybridization between E. fulgens and E. puniceoluteum is a common phenomenon. Bayesian assignment analysis detected the presence of F(1) and F(2) individuals and also signs of introgression, demonstrating a high potential for interspecific gene flow. Introgression occurs preferentially from E. fulgens to E. puniceoluteum. Pure parental individuals of both species display strong genotype-habitat associations, indicating that environment-dependent selection could be acting in all hybrid zones. This study suggests that hybridization and introgression are evolutionary processes playing a role in the diversification of Epidendrum and indicates the importance of investigations of hybrid zones in understanding reproductive barriers and speciation processes in Neotropical orchid species.

Phylogenetic analysis of beak and feather disease virus across a host ring-species complex

Pathogens have been hypothesized to play a major role in host diversity and speciation. Susceptibility of hybrid hosts to pathogens is thought to be a common phenomenon that could promote host population divergence and subsequently speciation. However, few studies have tested for pathogen infection across animal hybrid zones while testing for codivergence of the pathogens in the hybridizing host complex. Over 8 y, we studied natural infection by a rapidly evolving single-strand DNA virus, beak and feather diseases virus (BFDV), which infects parrots, exploiting a host-ring species complex (Platycercus elegans) in Australia. We found that host subspecies and their hybrids varied strikingly in both BFDV prevalence and load: both hybrid and phenotypically intermediate subspecies had lower prevalence and load compared with parental subspecies, while controlling for host age, sex, longitude and latitude, as well as temporal effects. We sequenced viral isolates throughout the range, which revealed patterns of genomic variation analogous to Mayr's ring-species hypothesis, to our knowledge for the first time in any host-pathogen system. Viral phylogeny, geographic location, intraspecific host density, and parrot community diversity and composition did not explain the differences in BFDV prevalence or load between subpopulations. Overall, our analyses suggest that functional host responses to infection, or force of infection, differ between subspecies and hybrids. Our findings highlight the role of host hybridization and clines in altering host-pathogen interactions, dynamics that can have important implications for models of speciation with gene flow, and offer insights into how pathogens may adapt to diverging host populations.

Avifauna do estado do Acre: composição, distribuição geográfica e conservação

O estado do Acre faz fronteiras internacionais com o Peru e a Bolívia e nacionais com os estados do Amazonas e Rondônia. O Acre está localizado nas terras baixas da Amazônia sul‐ocidental, próximo ao sopé dos Andes, dentro de uma região considerada megadiversa da Amazônia brasileira. Apesar disso, a região ainda é pouco conhecida e considerada prioritária para a realização de novos levantamentos biológicos. Com o intuito de contribuir para o conhecimento da avifauna do sudoeste amazônico, este estudo teve como principais objetivos responder as seguintes questões: (a) Quantas e quais são as espécies de aves do estado do Acre? (b) Como as espécies estão distribuídas dentro do estado do Acre? e (c) Qual o estado de conservação das espécies residentes no estado do Acre? A metodologia para responder a estas questões contemplou: (a) uma ampla revisão bibliográfica; (b) dois anos de levantamento em campo, incluindo registros e a coletas de espécimes testemunhos; (c) a confecção do mapa de distribuição de cada táxon (incluindo espécies e subespécies); (d) a distribuição dos táxons pelas três grandes regiões interfluviais do Estado (leste, central e oeste); (e) a identificação de zonas de contato e hibridização, baseada na distribuição dos táxons parapátridos dentro do Estado; (f) o cálculo da distribuição potencial dos táxons dentro do Acre, baseado na extrapolação da área ocupada por cada unidade ecológica (fitofisionomia) onde eles foram registrados; (g) o cálculo da meta de conservação de cada táxon residente no Estado e (h) uma análise de lacunas, baseada na sobreposição dos mapas de distribuição potencial de cada táxon com o das Áreas Protegidas do Estado. A análise de lacuna foi realizada tendo três diferentes cenários como referência: (a) primeiro cenário ‐ levou em consideração todas as Áreas Protegidas; (b) o segundo cenário – levou em consideração apenas as Unidades de Conservação ‐ UCs de Proteção Integral e (c) terceiro cenário – levou em consideração apenas as UCs de Proteção Integral + as de Uso Sustentável (exclusas as Terras Indígenas). A revisão bibliográfica e os levantamentos (históricos e de campo) tiveram início em agosto de 2005 e se estenderam até dezembro de 2007. Após o término da revisão bibliográfica e das expedições em campo, foram compilados 7.141 registros de aves para o todo o estado do Acre. Destes, 4. 623 são de espécimes coletados, dos quais, 2.295 (49,6%) são oriundos de coletas feitas durante a realização deste estudo. Confirmou‐se para o Acre a presença de 655 espécies biológicas, distribuídas em 73 Famílias e 23 Ordens. Como consequência direta deste estudo, cinco novas espécies foram acrescentadas à lista de aves brasileiras. Registrou‐se também, 59 espécies migratórias, das quais, 30 (50,8%) são migrantes neárticas, 11 (18,6%) foram consideradas como migrantes intratropicais e 18 (30,5%) como migrantes austrais. De todas as espécies registradas no Estado, 44 são endêmicas do centro de endemismo Inambari. Dos 556 táxons de aves florestais residentes no Acre, 72,8% (405) distribui‐se nas três sub‐regiões do Estado; 10,0% (56) foi registrado apenas na sub‐região oeste; 5,3% (30) apenas na sub‐região leste e 0,5% (03) apenas na sub‐região central. Ao menos seis pares de táxons irmãos apresentaram padrão de distribuição alopátrida e 15 conjuntos de táxons apresentaram distribuição parapátrida dentro do Estado. Foram identificadas duas zonas de contato secundário (leste/oeste) e duas possíveis zonas de hibridização (leste/oeste) dentro do Estado. As análises de lacunas mostraram que no primeiro cenário, 87,1% dos táxons atingiram 100% da meta de conservação; 12% ficaram em classes de conservação intermediárias, ou seja, em lacuna parcial de proteção e apenas 0,8% ficaram em lacuna total de proteção. No segundo cenário, apenas 0,6% dos táxons atingiram 100% da meta de conservação; 97,6% ficaram em lacuna parcial e 1,8% ficaram em lacuna total de proteção. No terceiro cenário, 73,5% dos táxons atingiram 100% da meta de conservação; 25,5% ficaram em lacuna parcial e apenas 0,8% ficaram em lacuna total de proteção. As principais conclusões obtidas a partir deste estudo foram: (a) que a riqueza avifaunística do estado do Acre é bastante expressiva, porém, o número de espécies detectadas deverá aumentar à medida que novos levantamentos forem realizados; (b) que os rios Purus e Juruá não são as barreiras físicas que determinam o padrão de distribuição da maioria das aves residentes no estado do Acre; (c) que a presença de zonas de contato secundário, não coincidentes com o curso dos dois principais rios do Estado, dá suporte a ideia de que fatores não ligados a uma barreira física devem estar atuando na manutenção do padrão de distribuição atual de alguns táxons de aves residentes no Acre; (d) que o número de espécies “desprotegidas” ou em “lacuna parcial de proteção” entre a avifauna do Acre é muito baixo quando todo o sistema de Áreas Protegidas é levado em consideração, porém este número aumenta com a exclusão das Unidades de Conservação de Uso Sustentável e das Terras Indígenas; (e) que as aves restritas às campinas e campinaranas do oeste do Acre são as únicas que se encontram em lacuna do sistema de áreas protegidas do Acre, indicando a necessidade de se criar uma ou mais Unidades de Conservação para proteger este habitat específico.

Variação geográfica, zonas de intergradação e especiação no complexo Icterus Cayanensis-Chrysocephalus (Aves : Icteridae)

O complexo Icterus cayanensis-chtysocephalus apresenta um intrincado padrão de variação em plumagem e tamanho corpóreo. São reconhecidos, tradicionalmente, para o grupo seis táxons: Icterus chtysocephalus, I. cayanensis cayanensís. I. cayanensis tibialis, I. cayanensis tibialis, I. cayanensis valenciobuenoi, I. cayanensis periporphyrus e I. cayanensis pyrrhopterus, que se substituem geograficamente ao longo de grande parte da América do Sul. Neste estudo foi feita a descrição dos padrões de variação geográfica. Foram diagnosticadas quatro espécies, à luz do conceito filogenético de espécie: Icterus cayanensis (Amazônia Meridional), Icterus chrysocephalus (Amazônia Setentrional), Icterus tibialis (Caatinga) e Icterus pyrrhopterus (Chaco); os táxons I. cayanensis valenciobuenoi e I. cayanensis periporphyrus foram sinonimizados. Entre as formas amazonicas (chrysocephalus e cayanensis) foi detectada a presença de uma zona híbrida mais extensa do que aquela reportada na literatura. No Brasil Central foi diagnosticada a maior zona de intergradação conhecida para aves, com aproximadamente 2.300 km de extensão, produto do intercruzamento entre Icterus tibialis e Icterus pyrrhopterus, formas distribuídas pela Caatinga e Chaco, respectivamente. Postula-se que as zonas de intergradação diagnosticadas neste estudo são produto do intercruzamento de populações previamente diferenciadas em isolamento geográfico.

Strong postzygotic isolation prevents introgression between two hybridizing Neotropical orchids, Epidendrum denticulatum and E-fulgens

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Muster und Prozesse in einer Microcebus-Hybridzone in Südostmadagaskar

Im Rahmen dieser Doktorarbeit wurde in zwei Schwerpunktanalysen mit eine Teil- und Gesamtdatensatz die Untersuchung der Hybridisierung zwischen den beiden Microcebus-Arten M. murinus und M. griseorufus im Ökoton Südostmadagaskars umfangreich und vertieft untersucht. Für die genetischen Analysen wurden die maternal vererbte mitochondriale Hypervariable Region I (HVR 1) und neun nukleäre biparental vererbte Mikrosatellitenmarker eingesetzt. Als weiterer Datensatz wurden morphometrische Daten verwendet. Für die erste Schwerpunktanalyse wurde ein bereits vorhandener Teildatensatz (Hapke 2005 & Gligor 2006) mit Daten von insgesamt 162 Individuen aus neun Populationen der Dornbuschzone, der Übergangswaldzone und des Küstenwaldgebietes eingesetzt. In der zweiten Schwerpunktanalyse wurde eine umfangreiche Untersuchung der Microcebus griseoruus-M. murinus- Hybridzone vorgenommen. Für diese detaillierte Charakterisierung der Hybridzone wurde eine ausgedehnte und fein auflösende Probennahme in einem als Kernzone definierten Bereich, der die gesamte Übergangswaldzone und die dazu benachbarten Dornbuschgebiete umfasste, durchgeführt. Die morphometrischen und genetischen Daten der neu beprobten Individuen dieser Kernzone wurden mit den Daten des Teildatensatzes und weiteren Daten aus Küstenwaldpopulationen (Hapke 2005) zu einem Gesamtdatensatz zusammengefasst. Die Integration des Teildatensatzes in den Gesamtdatensatz erforderte umfassende und zeitintensive Labor- und Analysearbeiten, die im Rahmen dieser Doktorarbeit durchgeführt wurden. Der Gesamtdatensatz umfasste insgesamt 569 Individuen der Gattung Microcebus aus 29 Untersuchungsstandorten. Die mit beiden Datensätzen durchgeführte Analyse morphometrischer Daten zeigte deutlich, dass die Mehrzahl der Individuen aus der Übergangswaldzone einen intermediären Morphotyp aufweist. Durch die mit den Daten des Teildatensatzes durchgeführten Bayes’schen Clusteranalysen und Assignment-Tests, das vornehmlich in den Populationen der Übergangszone beobachtete signifikante Kopplungsungleichgewicht und Heterozygotendefizit, die festgestellte Verteilung der mitochondrialen Haplotypen und das kontrastierende Muster zwischen nukleären Mikrosatellitengenotypen und mitochondrialen Haplotypen in den Übergangswaldpopulationen konnte erstmals das Vorkommen einer Hybridzone zwischen Microcebus-Arten wissenschaftlich fundiert festgestellt werden. Die Ergebnisse dieser Schwerpunktanalyse wurden in der Fachzeitschrift Molecular Ecology publiziert (Gligor et al. 2009). Die in der ersten Schwerpunktanalyse festgestellte Hybridzone konnte durch die zweite Schwerpunktanalyse mit den genetischen und morphometrischen Daten des Gesamtdatensatzes nicht nur bestätigt werden, sondern auch auf die gesamte Übergangswaldzone erweitert werden. Ferner wurden starke Hinweise auf eine Hybridisierung beider Microcebus-Arten an einigen Dornbuschstandorten der Kernzone gefunden. Durch die große Datenmenge des Gesamtdatensatzes, vor allem aus der Kernzone des Untersuchungsgebietes, war es möglich eine fundierte Charakterisierung der Microcebus griseoruus-M. murinus- Hybridzone durchzuführen. Die Übereinstimmung der Hybridzone mit dem beobachteten Vegetationsmosaik zusammen mit den Ergebnissen der PCA, der PCoA und der Bayes’schen Clusteranalyse sprechen für das Modell der „Mosaik Hybridzone“, während die Einzelbetrachtung der mosaikartig verteilten intermediären Übergangswälder eine hohe Abundanz der Hybride aufzeigte und somit eher das „Bounded Hybrid Superiority model“ unterstützt. Der gewählte geographische Beprobungsmaßstab könnte somit einen Einfluss auf die beobachtete Struktur einer Hybridzone haben. Eines der markantesten Muster in der Hybridzone ist das stark kontrastierende cyto-nukleäre Muster. Der seit ca. 3000 Jahren fortschreitende Klimawandel in Südmadagaskar und die damit verbundene Expansion des Verbreitungsgebietes der Art Microcebus griseorufus nach Osten, das in dieser Arbeit festgestellte „male-biased dispersal“ bei M. griseorufus und der Einfluss exogener Selektion sprechen stark für eine massive asymmetrische nukleäre Genintrogression von M. griseorufus-Allelen in M. murinus-Populationen, verbunden mit einer potentiellen Verdrängung der Art M. murinus aus der Übergangswaldzone. In den jeweiligen Kerngebieten Dornbusch und Küstenwald bleibt jedoch die Diskretheit beider Arten gewahrt.

Sex-specific clines support incipient speciation in a common European mammal

Hybrid zones provide excellent opportunities to study processes and mechanisms underlying reproductive isolation and speciation. Here we investigated sex-specific clines of molecular markers in hybrid zones of morphologically cryptic yet genetically highly-diverged evolutionary lineages of the European common vole (Microtus arvalis). We analyzed the position and width of four secondary contact zones along three independent transects in the region of the Alps using maternally (mitochondrial DNA) and paternally (Y-chromosome) inherited genetic markers. Given male-biased dispersal in the common vole, a selectively neutral secondary contact would show broader paternal marker clines than maternal ones. In a selective case, for example, involving a form of Haldane’s rule, Y-chromosomal clines would not be expected to be broader than maternal markers because they are transmitted by the heterogametic sex and thus gene flow would be restricted. Consistent with the selective case, paternal clines were significantly narrower or at most equal in width to maternal clines in all contact zones. In addition, analyses using maximum likelihood cline-fitting detected a shift of paternal relative to maternal clines in three of four contact zones. These patterns suggest that processes at the contact zones in the common vole are not selectively neutral, and that partial reproductive isolation is already established between these evolutionary lineages. We conclude that hybrid zone movement, sexual selection and/or genetic incompatibilities are likely associated with an unusual unidirectional manifestation of Haldane’s rule in this common European mammal.

Hybridization and adaptive radiation

Whether interspecific hybridization is important as a mechanism that generates biological diversity is a matter of controversy. Whereas some authors focus on the potential of hybridization as a source of genetic variation, functional novelty and new species, others argue against any important role, because reduced fitness would typically render hybrids an evolutionary dead end. By drawing on recent developments in the genetics and ecology of hybridization and on principles of ecological speciation theory, I develop a concept that reconciles these views and adds a new twist to this debate. Because hybridization is common when populations invade new environments and potentially elevates rates of response to selection, it predisposes colonizing populations to rapid adaptive diversification under disruptive or divergent selection. I discuss predictions and suggest tests of this hybrid swarm theory of adaptive radiation and review published molecular phylogenies of adaptive radiations in light of the theory. Some of the confusion about the role of hybridization in evolutionary diversification stems from the contradiction between a perceived necessity for cessation of gene flow to enable adaptive population differentiation on the one hand [1], and the potential of hybridization for generating adaptive variation, functional novelty and new species 2, 3 and 4 on the other. Much progress in the genetics 5, 6, 7, 8 and 9 and ecology of hybridization 9, 10 and 11, and in our understanding of the role of ecology in speciation (see Glossary) 12, 13 and 14 make a re-evaluation timely. Whereas botanists traditionally stressed the diversity-generating potential of hybridization 2, 3 and 14, zoologists traditionally saw it as a process that limits diversification [1] and refer to it mainly in the contexts of hybrid zones (Box 1) and reinforcement of reproductive isolation [15]. Judging by the wide distribution of allopolyploidy among plants, many plant species might be of direct hybrid origin or descended from a hybrid species in the recent past [16]. The ability to reproduce asexually might explain why allopolyploid hybrid species are more common in plants than in animals. Allopolyploidy arises when meiotic mismatch of parental chromosomes or karyotypes causes hybrid sterility. Mitotic error, duplicating the karyotype, can restore an asexually maintained hybrid line to fertility. Although bisexual allopolyploid hybrid species are not uncommon in fish [17] and frogs [18], the difficulty with which allopolyploid animals reproduce, typically requiring gynogenesis[19], makes establishment and survival of allopolyploid animal species difficult.

Gene flow and the maintenance of species boundaries

Hybrid zones are regions where individuals from genetically differentiated populations meet and mate, resulting in at least some offspring of mixed ancestry. Patterns of gene flow (introgression) in hybrid zones vary across the genome, allowing assessment of the role of individual genes or genome regions in reproductive isolation. Here, we document patterns of introgression between two recently diverged species of field crickets. We sampled at a very fine spatial scale and genotyped crickets for 110 highly differentiated single nucleotide polymorphisms (SNPs) identified through transcriptome scans. Using both genomic and geographic cline analysis, we document remarkably abrupt transitions (<100 m) in allele frequencies for 50 loci, despite high levels of gene flow at other loci. These are among the steepest clines documented for any hybridizing taxa. Furthermore, the cricket hybrid zone provides one of the clearest examples of the semi-permeability of species boundaries. Comparisons between data from the fine-scale transect and data (for the same set of markers) from sampling a much larger area in a different region of the cricket hybrid zone reveal consistent patterns of introgression for individual loci. The consistency in patterns of introgression between these two distant and distinct regions of the hybrid zone suggests that strong selection is acting to maintain abrupt discontinuities within the hybrid zone and that genomic regions with restricted introgression likely include genes that contribute to nonecological prezygotic barriers.

Tracing reinforcement through asymmetrical partner preference in the European common vole Microtus arvalis

Background The mechanistic basis of speciation and in particular the contribution of behaviour to the completion of the speciation process is often contentious. Contact zones between related taxa provide a situation where selection against hybridization might reinforce separation by behavioural mechanisms, which could ultimately fully isolate the taxa. One of the most abundant European mammals, the common vole Microtus arvalis, forms multiple natural hybrid zones where rapidly diverging evolutionary lineages meet in secondary contact. Very narrow zones of hybridization spanning only a few kilometres and sex-specific gene flow patterns indicate reduced fitness of natural hybrids and incipient speciation between some of the evolutionary lineages. In this study, we examined the contribution of behavioural mechanisms to the speciation process in these rodents by fine-mapping allopatric and parapatric populations in the hybrid zone between the Western and Central lineages and experimental testing of the partner preferences of wild, pure-bred and hybrid female common voles. Results Genetic analysis based on microsatellite markers revealed the presence of multiple parapatric and largely non-admixed populations at distances of about 10 km at the edge of the area of natural hybridization between the Western and Central lineages. Wild females from Western parapatric populations and lab-born F1 hybrids preferred males from the Western lineage whereas wild females of Central parapatric origin showed no measurable preference. Furthermore, wild and lab-born females from allopatric populations of the Western or Central lineages showed no detectable preference for males from either lineage. Conclusions The detected partner preferences are consistent with asymmetrical reinforcement of pre-mating reproductive isolation mechanisms in the European common vole and with earlier results suggesting that hybridization is more detrimental to the Western lineage. As a consequence, these differences in behaviour might contribute to a further geographical stabilization of this moving hybrid zone. Such behavioural processes could also provide a mechanistic perspective for frequently-detected asymmetrical introgression patterns in the largely allopatrically diversifying Microtus genus and other rapidly speciating rodents.