RIDGE, a tool tailored to detect gene flow barriers across species pairs

Characterizing the processes underlying reproductive isolation between diverging lineages is central to understanding speciation. Here, we present RIDGE-Reproductive Isolation Detection using Genomic polymorphisms-a tool tailored for quantifying gene flow barrier proportion and identifying the relevant genomic regions. RIDGE relies on an Approximate Bayesian Computation with a model-averaging approach to accommodate diverse scenarios of lineage divergence. It captures heterogeneity in effective migration rate along the genome while accounting for variation in linked selection and recombination. The barrier detection test relies on numerous summary statistics to compute a Bayes factor, offering a robust statistical framework that facilitates cross-species comparisons. Simulations revealed RIDGE's efficiency in capturing signals of ongoing migration. Model averaging proved particularly valuable in scenarios of high model uncertainty where no migration or migration homogeneity can be wrongly assumed, typically for recent divergence times <0.1 2Ne generations. Applying RIDGE to four published crow data sets, we first validated our tool by identifying a well-known large genomic region associated with mate choice patterns. Second, while we identified a significant overlap of outlier loci using RIDGE and traditional genomic scans, our results suggest that a substantial portion of previously identified outliers are likely false positives. Outlier detection relies on allele differentiation, relative measures of divergence and the count of shared polymorphisms and fixed differences. Our analyses also highlight the value of incorporating multiple summary statistics including our newly developed outlier ones that can be useful in challenging detection conditions.

Spatiotemporal dynamics in butterfly hybrid zones

Evaluating whether hybrid zones are stable or mobile can provide novel insights for evolution and conservation biology. Butterflies exhibit high sensitivity to environmental changes and represent an important model system for the study of hybrid zone origins and maintenance. Here, we review the literature exploring butterfly hybrid zones, with a special focus on their spatiotemporal dynamics and the potential mechanisms that could lead to their movement or stability. We then compare different lines of evidence used to investigate hybrid zone dynamics and discuss the strengths and weaknesses of each approach. Our goal with this review is to reveal general conditions associated with the stability or mobility of butterfly hybrid zones by synthesizing evidence obtained using different types of data sampled across multiple regions and spatial scales. Finally, we discuss spatiotemporal dynamics in the context of a speciation/divergence continuum, the relevance of hybrid zones for conservation biology, and recommend key topics for future investigation.

Acknowledging more biodiversity without more species

Delimiting and naming biodiversity is a vital step toward wildlife conservation and research. However, species delimitation must be consistent across biota so that the limited resources available for nature protection can be spent effectively and objectively. To date, newly discovered lineages typically are either left undescribed and thus remain unprotected or are being erroneously proposed as new species despite mixed evidence for completed speciation, in turn contributing to the emerging problem of taxonomic inflation. Inspired by recent conceptual and methodological progress, we propose a standardized workflow for species delimitation that combines phylogenetic and hybrid zone analyses of genomic datasets ("genomic taxonomy"), in which phylogeographic lineages that do not freely admix are ranked as species, while those that have remained fully genetically compatible are ranked as subspecies. In both cases, we encourage their formal taxonomic naming, diagnosis, and description to promote social awareness toward biodiversity. The use of loci throughout the genome overcomes the unreliability of widely used barcoding genes when phylogeographic patterns are complex, while the evaluation of divergence and reproductive isolation unifies the long-opposed concepts of lineage species and biological species. We suggest that a shift in conservation assessments from a single level (species) toward a two-level hierarchy (species and subspecies) will lead to a more balanced perception of biodiversity in which both intraspecific and interspecific diversity are valued and more adequately protected.

Large X-effects are absent in torrent frogs with nascent sex chromosomes

Sex chromosomes are popularized as a special role in driving speciation. However, the empirical evidence from natural population processes has been limited to organisms with degenerated sex chromosomes, where hemizygosity is mainly considered to act as the driver of reproductive isolation. Here, we examined several hybrid zones of torrent frog Amolops mantzorum species complex, using an approach by mapping species-diagnostic loci onto the reference genome to compare sex-linked versus autosomal patterns of introgression. We find little support in sex-linked incompatibilities for large X-effects for these populations in hybrid zones with homomorphic sex chromosomes, due to the absence of the hemizygous effects. As expected, the large X-effects were not found in those with heteromorphic but newly evolved sex chromosomes, owing to the absence of strong genetic differences between X and Y chromosomes. The available data so far on amphibians suggest little role for sex-linked genes in speciation. The large X-effects in those with nascent sex chromosomes may not be as ubiquitous as presumed across the animal kingdom.

An allozyme polymorphism is associated with a large chromosomal inversion in the marine snail Littorina fabalis

Understanding the genetic targets of natural selection is one of the most challenging goals of population genetics. Some of the earliest candidate genes were identified from associations between allozyme allele frequencies and environmental variation. One such example is the clinal polymorphism in the arginine kinase (Ak) gene in the marine snail Littorina fabalis. While other enzyme loci do not show differences in allozyme frequencies among populations, the Ak alleles are near differential fixation across repeated wave exposure gradients in Europe. Here, we use this case to illustrate how a new sequencing toolbox can be employed to characterize the genomic architecture associated with historical candidate genes. We found that the Ak alleles differ by nine nonsynonymous substitutions, which perfectly explain the different migration patterns of the allozymes during electrophoresis. Moreover, by exploring the genomic context of the Ak gene, we found that the three main Ak alleles are located on different arrangements of a putative chromosomal inversion that reaches near fixation at the opposing ends of two transects covering a wave exposure gradient. This shows Ak is part of a large (3/4 of the chromosome) genomic block of differentiation, in which Ak is unlikely to be the only target of divergent selection. Nevertheless, the nonsynonymous substitutions among Ak alleles and the complete association of one allele with one inversion arrangement suggest that the Ak gene is a strong candidate to contribute to the adaptive significance of the inversion.

Revisiting a classic hybrid zone: Movement of the northern flicker hybrid zone in contemporary times

Natural hybrid zones have provided important insights into the evolutionary process, and their geographic dynamics over time can help to disentangle the underlying biological processes that maintain them. Here, we leverage replicated sampling of an identical transect across the hybrid zone between yellow-shafted and red-shafted flickers in the Great Plains to assess its stability over ∼60 years (1955-1957 to 2016-2018). We identify a ∼73-km westward shift in the hybrid zone center toward the range of the red-shafted flicker, but find no associated changes in width over our sampling period. In fact, the hybrid zone remains remarkably narrow, suggesting some kind of selective pressure maintains the zone. By comparing to previous work in the same geographic region, it appears likely that the movement in the hybrid zone has occurred in the years since the early 1980s. This recent movement may be related to changes in climate or land management practices that have allowed westward movement of yellow-shafted flickers into the Great Plains.

Morphometry of two cryptic tree frog species at their hybrid zone reveals neither intermediate nor transgressive morphotypes

Under incomplete reproductive isolation, secondary contact of diverged allopatric lineages may lead to the formation of hybrid zones that allow to study recombinants over several generations as excellent systems of genomic interactions resulting from the evolutionary forces acting on certain genes and phenotypes. Hybrid phenotypes may either exhibit intermediacy or, alternatively, transgressive traits, which exceed the extremes of their parents due to epistasis and segregation of complementary alleles. While transgressive morphotypes have been examined in fish, reptiles, birds, and mammals, studies in amphibians are rare. Here, we associate microsatellite-based genotypes with morphometrics-based morphotypes of two tree frog species of the Hyla arborea group, sampled across a hybrid zone in Poland, to understand whether the genetically differentiated parental species also differ in morphology between each other and their hybrids and whether secondary contact leads to the evolution of intermediate or transgressive morphotypes. Using univariate approaches, explorative multivariate methods (principal component analyses) as well as techniques with prior grouping (discriminant function analyses), we find that morphotypes of both parental species and hybrids differ from each other. Importantly, hybrid morphotypes are neither intermediate nor transgressive but found to be more similar to H. orientalis than to H. arborea.

Consequences of Hybridization in Mammals: A Systematic Review

Hybridization, defined as breeding between two distinct taxonomic units, can have an important effect on the evolutionary patterns in cross-breeding taxa. Although interspecific hybridization has frequently been considered as a maladaptive process, which threatens species genetic integrity and survival via genetic swamping and outbreeding depression, in some cases hybridization can introduce novel adaptive variation and increase fitness. Most studies to date focused on documenting hybridization events and analyzing their causes, while relatively little is known about the consequences of hybridization and its impact on the parental species. To address this knowledge gap, we conducted a systematic review of studies on hybridization in mammals published in 2010-2021, and identified 115 relevant studies. Of 13 categories of hybridization consequences described in these studies, the most common negative consequence (21% of studies) was genetic swamping and the most common positive consequence (8%) was the gain of novel adaptive variation. The total frequency of negative consequences (49%) was higher than positive (13%) and neutral (38%) consequences. These frequencies are biased by the detection possibilities of microsatellite loci, the most common genetic markers used in the papers assessed. As negative outcomes are typically easier to demonstrate than positive ones (e.g., extinction vs hybrid speciation), they may be over-represented in publications. Transition towards genomic studies involving both neutral and adaptive variation will provide a better insight into the real impacts of hybridization.

Disentangling the components of triploid block and its fitness consequences in natural diploid-tetraploid contact zones of Arabidopsis arenosa

Hybrid seed inviability (HSI) is an important mechanism of reproductive isolation and speciation. HSI varies in strength among populations of diploid species but it remains to be tested whether similar processes affect natural variation in HSI within ploidy-variable species (triploid block). Here we used extensive endosperm, seed and F₁ -hybrid phenotyping to explore HSI variation within a diploid-autotetraploid species. By leveraging 12 population pairs from three ploidy contact zones, we tested for the effect of interploidy crossing direction (parent of origin), ploidy divergence and spatial arrangement in shaping reproductive barriers in a naturally relevant context. We detected strong parent-of-origin effects on endosperm development, F₁ germination and survival, which was also reflected in the rates of triploid formation in the field. Endosperm cellularization failure was least severe and F₁ -hybrid performance was slightly better in the primary contact zone, with genetically closest diploid and tetraploid lineages. We demonstrated overall strong parent-of-origin effects on HSI in a ploidy variable species, which translate to fitness effects and contribute to interploidy reproductive isolation in a natural context. Subtle intraspecific variation in these traits suggests the fitness consequences of HSI are predominantly a constitutive property of the species regardless of the evolutionary background of its populations.

The genic view of hybridization in the Anthropocene

Human impact is noticeable around the globe, indicating that a new era might have begun: the Anthropocene. Continuing human activities, including land-use changes, introduction of non-native species and rapid climate change, are altering the distributions of countless species, often giving rise to human-mediated hybridization events. While the interbreeding of different populations or species can have detrimental effects, such as genetic extinction, it can be beneficial in terms of adaptive introgression or an increase in genetic diversity. In this paper, I first review the different mechanisms and outcomes of anthropogenic hybridization based on literature from the last five years (2016-2020). The most common mechanisms leading to the interbreeding of previously isolated taxa include habitat change (51% of the studies) and introduction of non-native species (34% intentional and 19% unintentional). These human-induced hybridization events most often result in introgression (80%). The high incidence of genetic exchange between the hybridizing taxa indicates that the application of a genic view of speciation (and introgression) can provide crucial insights on how to address hybridization events in the Anthropocene. This perspective considers the genome as a dynamic collection of genetic loci with distinct evolutionary histories, giving rise to a heterogenous genomic landscape in terms of genetic differentiation and introgression. First, understanding this genomic landscape can lead to a better selection of diagnostic genetic markers to characterize hybrid populations. Second, describing how introgression patterns vary across the genome can help to predict the likelihood of negative processes, such as demographic and genetic swamping, as well as positive outcomes, such as adaptive introgression. It is especially important to not only quantify how much genetic material introgressed, but also what has been exchanged. Third, comparing introgression patterns in pre-Anthropocene hybridization events with current human-induced cases might provide novel insights into the likelihood of genetic swamping or species collapse during an anthropogenic hybridization event. However, this comparative approach remains to be tested before it can be applied in practice. Finally, the genic view of introgression can be combined with conservation genomic studies to determine the legal status of hybrids and take appropriate measures to manage anthropogenic hybridization events. The interplay between evolutionary and conservation genomics will result in the constant exchange of ideas between these fields which will not only improve our knowledge on the origin of species, but also how to conserve and protect them.

A holocenic and dynamic hybrid zone between two cactophilic Drosophila species in a coastal lowland plain of the Brazilian Atlantic Forest

Hybridization and introgression are processes that contribute to shaping biological diversity. The factors promoting the formation of these processes are multiples but poorly explored in a biogeographical and ecological context. In the southeast coastal plain of the Brazilian Atlantic Forest, a hybrid zone was described between two closely related cactophilic species, Drosophila antonietae and D. serido. Here, we revisited and analysed specimens from this hybrid zone to evaluate its temporal and spatial dynamic. We examined allopatric and sympatric populations of the flies using independent sources of data such as mitochondrial and nuclear sequences, microsatellite loci, morphometrics of wings and male genitalia, and climatic niche models. We also verified the emergence of the flies from necrotic tissues of collected cacti to verify the role of host association for the population dynamics. Our results support the existence of a hybrid zone due to secondary contact and limited to the localities where the two species are currently in contact. Furthermore, we detected asymmetric bidirectional introgression and the maintenance of the species integrity, ecological association and morphological characters, suggesting selection and limited introgression. Considering our paleomodels, probably this hybrid zone is recent and the contact occurred during the Holocene to the present day, favoured by range expansion of their populations due to expansion of open and dry areas in eastern South America during palaeoclimatic and geomorphological events.

Demographic history shapes genomic ancestry in hybrid zones

Demographic factors such as migration rate and population size can impede or facilitate speciation. In hybrid zones, reproductive boundaries between species are tested and demography mediates the opportunity for admixture between lineages that are partially isolated. Genomic ancestry is a powerful tool for revealing the history of admixed populations, but models and methods based on local ancestry are rarely applied to structured hybrid zones. To understand the effects of demography on ancestry in hybrids zones, we performed individual-based simulations under a stepping-stone model, treating migration rate, deme size, and hybrid zone age as parameters. We find that the number of ancestry junctions (the transition points between genomic regions with different ancestries) and heterogenicity (the genomic proportion heterozygous for ancestry) are often closely connected to demographic history. Reducing deme size reduces junction number and heterogenicity. Elevating migration rate increases heterogenicity, but migration affects junction number in more complex ways. We highlight the junction frequency spectrum as a novel and informative summary of ancestry that responds to demographic history. A substantial proportion of junctions are expected to fix when migration is limited or deme size is small, changing the shape of the spectrum. Our findings suggest that genomic patterns of ancestry could be used to infer demographic history in hybrid zones.

Salinity, not genetic incompatibilities, limits the establishment of the invasive hybrid cattail Typha × glauca in coastal wetlands

Hybrids of a single pair of parent species can be much more common in some geographical regions than in others. The reasons for this are not well understood, but could help explain processes such as species diversification or the range expansion of invasive hybrids. The widespread cattails Typha latifolia and T. angustifolia seldom hybridize in some parts of their range, but in other areas produce the dominant hybrid T. × glauca. We used a combination of field and greenhouse experiments to investigate why T. × glauca has invaded wetlands in the Laurentian Great Lakes region of southern Ontario, Canada, but is much less common in the coastal wetlands of Nova Scotia (NS) in eastern Canada. One potentially important environmental difference between these two regions is salinity. We therefore tested three hypotheses: (1) T. latifolia and T. angustifolia in NS are genetically incompatible; (2) the germination or growth of T. × glauca is reduced by salinity; and (3) T. latifolia, a main competitor of T. × glauca, is locally adapted to saline conditions in NS. Our experiments showed that NS T. latifolia and T. angustifolia are genetically compatible, and that saline conditions do not impede growth of hybrid plants. However, we also found that under conditions of high salinity, germination rates of hybrid seeds were substantially lower than those of NS T. latifolia. In addition, germination rates of NS T. latifolia were higher than those of Ontario T. latifolia, suggesting local adaptation to salinity in coastal wetlands. This study adds to the growing body of literature which identifies the important roles that local habitat and adaptation can play in the distributions and characteristics of hybrid zones.

Salinity, not genetic incompatibilities, limits the establishment of the invasive hybrid cattail Typha × glauca in coastal wetlands

Hybrids of a single pair of parent species can be much more common in some geographical regions than in others. The reasons for this are not well understood, but could help explain processes such as species diversification or the range expansion of invasive hybrids. The widespread cattails Typha latifolia and T. angustifolia seldom hybridize in some parts of their range, but in other areas produce the dominant hybrid T. × glauca. We used a combination of field and greenhouse experiments to investigate why T. × glauca has invaded wetlands in the Laurentian Great Lakes region of southern Ontario, Canada, but is much less common in the coastal wetlands of Nova Scotia (NS) in eastern Canada. One potentially important environmental difference between these two regions is salinity. We therefore tested three hypotheses: (1) T. latifolia and T. angustifolia in NS are genetically incompatible; (2) the germination or growth of T. × glauca is reduced by salinity; and (3) T. latifolia, a main competitor of T. × glauca, is locally adapted to saline conditions in NS. Our experiments showed that NS T. latifolia and T. angustifolia are genetically compatible, and that saline conditions do not impede growth of hybrid plants. However, we also found that under conditions of high salinity, germination rates of hybrid seeds were substantially lower than those of NS T. latifolia. In addition, germination rates of NS T. latifolia were higher than those of Ontario T. latifolia, suggesting local adaptation to salinity in coastal wetlands. This study adds to the growing body of literature which identifies the important roles that local habitat and adaptation can play in the distributions and characteristics of hybrid zones.

How do species barriers decay? Concordance and local introgression in mosaic hybrid zones of mussels

The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study 'replicated' instances of secondary contact between closely related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry-informative panel of such SNPs. We then compared their frequencies in newly sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bi-stable variants (Dobzhansky-Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact.

Distinct patterns of hybridization across a suture zone in a coral reef fish (Dascyllus trimaculatus)

Hybrid zones are natural laboratories for investigating the dynamics of gene flow, reproductive isolation, and speciation. A predominant marine hybrid (or suture) zone encompasses Christmas Island (CHR) and Cocos (Keeling) Islands (CKE), where 15 different instances of interbreeding between closely related species from Indian and Pacific Oceans have been documented. Here, we report a case of hybridization between genetically differentiated Pacific and Indian Ocean lineages of the three-spot dascyllus, Dascyllus trimaculatus (Rüppell, 1829). Field observations indicate there are subtle color differences between Pacific and Indian Ocean lineages. Most importantly, population densities of color morphs and genetic analyses (mitochondrial DNA and SNPs obtained via RADSeq) suggest that the pattern of hybridization within the suture zone is not homogeneous. At CHR, both color morphs were present, mitochondrial haplotypes of both lineages were observed, and SNP analyses revealed both pure and hybrid genotypes. Meanwhile, in CKE, the Indian Ocean color morphs were prevalent, only Indian Ocean mitochondrial haplotypes were observed, and SNP analysis showed hybrid individuals with a large proportion (~80%) of their genotypes assigning to the Indian Ocean lineage. We conclude that CHR populations are currently receiving an influx of individuals from both ocean basins, with a greater influence from the Pacific Ocean. In contrast, geographically isolated CKE populations appear to be self-recruiting and with more influx of individuals from the Indian Ocean. Our research highlights how patterns of hybridization can be different at scales of hundreds of kilometers, due to geographic isolation and the history of interbreeding between lineages.

Complex patterns of differentiation and gene flow underly the divergence of aposematic phenotypes in Oophaga poison frogs

Hybridization and introgression can have complex consequences for both species evolution and conservation. Here, we investigated the origin and characteristics of a putative hybrid zone between two South American poison dart frog species, Oophaga anchicayensis and the critically endangered Oophaga lehmanni, which are heavily sought after on the illegal pet market. Using a combination of phenotypic (49 traits) and genomic (ddRADseq) data, we found that the putative hybrids are morphologically distinct from their parental species and confirmed genomic signatures of admixture in these populations. Several lines of evidence (hybrid indices, interspecific hybrid heterozygosity, genomic clines, comparisons with simulated hybrids and demographic modelling) support the conclusion that these populations are not comprised of early-generation hybrids and thus, they probably did not arise as a result of illegal translocations associated with wildlife trafficking. Instead, they probably represent an independent lineage which has persisted through isolation and has only relatively recently re-established gene flow with both parental species. Furthermore, we detected signals of differential introgression from parental species into these hybrid populations which suggest relaxed stabilizing selection on these aposematic colour morphs, potentially via context-dependent female choice. These populations thus provide a fascinating window into the role of hybridization, isolation and female choice in the diversification of South American poison dart frogs. In addition, our results underline the importance of landscape conservation measures to protect, not only known localities of nominal species, but also the phenotypic and genomic variation harbored by admixed lineages which represent crucial repositories for the impressive diversity in this system.

Tracing the footprints of a moving hybrid zone under a demographic history of speciation with gene flow

A lack of optimal gene combinations, as well as low levels of genetic diversity, is often associated with the formation of species range margins. Conservation efforts rely on predictive modelling using abiotic variables and assessments of genetic diversity to determine target species and populations for controlled breeding, germplasm conservation and assisted migration. Biotic factors such as interspecific competition and hybridization, however, are largely ignored, despite their prevalence across diverse taxa and their role as key evolutionary forces. Hybridization between species with well-developed barriers to reproductive isolation often results in the production of offspring with lower fitness. Generation of novel allelic combinations through hybridization, however, can also generate positive fitness consequences. Despite this possibility, hybridization-mediated introgression is often considered a threat to biodiversity as it can blur species boundaries. The contribution of hybridization towards increasing genetic diversity of populations at range margins has only recently gathered attention in conservation studies. We assessed the extent to which hybridization contributes towards range dynamics by tracking spatio-temporal changes in the central location of a hybrid zone between two recently diverged species of pines: Pinus strobiformis and P. flexilis. By comparing geographic cline centre estimates for global admixture coefficient with morphological traits associated with reproductive output, we demonstrate a northward shift in the hybrid zone. Using a combination of spatially explicit, individual-based simulations and linkage disequilibrium variance partitioning, we note a significant contribution of adaptive introgression towards this northward movement, despite the potential for differences in regional population size to aid hybrid zone movement. Overall, our study demonstrates that hybridization between recently diverged species can increase genetic diversity and generate novel allelic combinations. These novel combinations may allow range margin populations to track favourable climatic conditions or facilitate adaptive evolution to ongoing and future climate change.

Temporal and spatial dynamics of competitive parapatry in chewing lice

We synthesize observations from 1979 to 2016 of a contact zone involving two subspecies of pocket gophers (Thomomys bottae connectens and T. b. opulentus) and their respective chewing lice (Geomydoecus aurei and G. centralis) along the Río Grande Valley in New Mexico, U.S.A., to test predictions about the dynamics of the zone. Historically, the natural flood cycle of the Rio Grande prevented contact between the two subspecies of pocket gophers. Flood control measures completed in the 1930s permitted contact, thus establishing the hybrid zone between the pocket gophers and the contact zone between their lice (without hybridization). Since that time, the pocket gopher hybrid zone has stabilized, whereas the northern chewing louse species has replaced the southern louse species at a consistent rate of ~150 m/year. The 0.2-0.8 width of the replacement zone has remained constant, reflecting the constant rate of chewing louse species turnover on a single gopher and within a local pocket gopher population. In contrast, the full width of the replacement zone (northernmost G. centralis to southernmost G. aurei) has increased annually. By employing a variety of metrics of the species replacement zone, we are better able to understand the dynamics of interactions between and among the chewing lice and their pocket gopher hosts. This research provides an opportunity to observe active species replacement and resulting distributional shifts in a parasitic organism in its natural setting.

Homoploid hybridization of plants in the Hengduan mountains region

The Hengduan Mountains Region (HMR) is a major global biodiversity hotspot. Complex tectonic and historical climatic conditions created opportunities for natural interspecific hybridization. Likewise, anthropogenic disturbance potentially raises the frequency of hybridization. Among species studies to date, the frequency of homoploid hybridization appears in the HMR. Of nine taxa in which natural hybridization has been detected, three groups are involved in homoploid hybrid speciation, and species pairs from the remaining six genera suggest that continuous gene flow occurs in hybrid zones. Reproductive isolation may greatly affect the dynamic and architecture of hybrid zones in the HMR. Asymmetrical hybridization and introgression can primarily be attributed to both prezygotic and postzygotic barriers. The frequent observation of such asymmetry may imply that reproductive barrier contributes to maintaining species boundaries in the alpine region. Ecological isolations with environmental disturbance may promote breeding barriers between parental species and hybrids. Hybrid zones may be an important phase for homoploid hybrid speciation. Hybrid zones potentially provided abundant genetic resources for the diversification of the HMR flora. The ecological and molecular mechanisms of control and mediation for natural hybridization will help biologists to understand the formation of biodiversity in the HMR. More researches from ecological and molecular aspects were required in future studies.

Phylogeography of a cryptic speciation continuum in Eurasian spadefoot toads (Pelobates)

Cryptic phylogeographic diversifications provide unique models to examine the role of phylogenetic divergence on the evolution of reproductive isolation, without extrinsic factors such as ecological and behavioural differentiation. Yet, to date very few comparative studies have been attempted within such radiations. Here, we characterize a new speciation continuum in a group of widespread Eurasian amphibians, the Pelobates spadefoot toads, by conducting multilocus (restriction site associated DNA sequencing and mitochondrial DNA) phylogenetic, phylogeographic and hybrid zone analyses. Within the P. syriacus complex, we discovered species-level cryptic divergences (>5 million years ago [My]) between populations distributed in the Near-East (hereafter P. syriacus sensu stricto [s.s.]) and southeastern Europe (hereafter P. balcanicus), each featuring deep intraspecific lineages. Altogether, we could scale hybridizability to divergence time along six different stages, spanning from sympatry without gene flow (P. fuscus and P. balcanicus, >10 My), parapatry with highly restricted hybridization (P. balcanicus and P. syriacus s.s., >5 My), narrow hybrid zones (~15 km) consistent with partial reproductive isolation (P. fuscus and P. vespertinus, ~3 My), to extensive admixture between Pleistocene and refugial lineages (≤2 My). This full spectrum empirically supports a gradual build up of reproductive barriers through time, reversible up until a threshold that we estimate at ~3 My. Hence, cryptic phylogeographic lineages may fade away or become reproductively isolated species simply depending on the time they persist in allopatry, and without definite ecomorphological divergence.

Phylogeography and species limits in the red-shouldered hawk (Buteo lineatus): Characterization of the Northern Florida Suture Zone in birds

The North American red-shouldered hawk, Buteo lineatus, is comprised of two widely allopatric eastern and western populations with an additional well-marked subspecies in the Florida peninsula. The two eastern populations meet in northern Florida, the location of a well-known suture zone in many nonavian organisms. We sequenced the complete mitochondrial ND2 gene and two nuclear introns to investigate its genetic population structure and species status. No mitochondrial haplotypes were shared between the eastern and western populations, and genetic variance among 14 populations was 0.42; almost all of this (0.40) was distributed among the three regions. A clade of haplotypes very common in the Florida peninsula decreased in frequency elsewhere and, when modeled as a hybrid zone, had an estimated width of 1,158 km with a center near Ocala, FL. Ecological niche modeling suggests the western, eastern, and Florida peninsula populations were geographically isolated during the last glacial maximum. We consider these to represent three phylogenetic species. A coalescent analysis incorporating incomplete lineage sorting and gene tree uncertainty also suggested the divergence between the western and eastern populations is consistent with species-level divergence. With the addition of this hawk, four avian species are now known to hybridize along the Gulf Coast of the United States in or near the Northern Florida Suture Zone. The widths of these avian zones vary substantially (176-1,158 km) and appear to reflect magnitude of gene flow, rather than extent of genetic differentiation. None of these birds was suggested as possible exemplars in the original description of the suture zone. Of the six species that were so identified, three have been surveyed to date, but none of those was found to be genetically differentiated.

Morphologically cryptic Amazonian bird species pairs exhibit strong postzygotic reproductive isolation

We possess limited understanding of how speciation unfolds in the most species-rich region of the planet-the Amazon basin. Hybrid zones provide valuable information on the evolution of reproductive isolation, but few studies of Amazonian vertebrate hybrid zones have rigorously examined the genome-wide underpinnings of reproductive isolation. We used genome-wide genetic datasets to show that two deeply diverged, but morphologically cryptic sister species of forest understorey birds show little evidence for prezygotic reproductive isolation, but substantial postzygotic isolation. Patterns of heterozygosity and hybrid index revealed that hybrid classes with heavily recombined genomes are rare and closely match simulations with high levels of selection against hybrids. Genomic and geographical clines exhibit a remarkable similarity across loci in cline centres, and have exceptionally narrow cline widths, suggesting that postzygotic isolation is driven by genetic incompatibilities at many loci, rather than a few loci of strong effect. We propose Amazonian understorey forest birds speciate slowly via gradual accumulation of postzygotic genetic incompatibilities, with prezygotic barriers playing a less important role. Our results suggest old, cryptic Amazonian taxa classified as subspecies could have substantial postzygotic isolation deserving species recognition and that species richness is likely to be substantially underestimated in Amazonia.

Chorthippus parallelus and Wolbachia: Overlapping Orthopteroid and Bacterial Hybrid Zones

Wolbachia is a well-known endosymbiotic, strictly cytoplasmic bacterium. It establishes complex cytonuclear relations that are not necessarily deleterious to its host, but that often result in reproductive alterations favoring bacterial transmission. Among these alterations, a common one is the cytoplasmic incompatibility (CI) that reduces the number of descendants in certain crosses between infected and non-infected individuals. This CI induced by Wolbachia appears in the hybrid zone that the grasshoppers Chorthippus parallelus parallelus (Cpp) and C. p. erythropus (Cpe) form in the Pyrenees: a reputed model in evolutionary biology. However, this cytonuclear incompatibility is the result of sophisticated processes of the co-divergence of the genomes of the bacterial strains and the host after generations of selection and coevolution. Here we show how these genome conflicts have resulted in a finely tuned adjustment of the bacterial strain to each pure orthopteroid taxon, and the striking appearance of another, newly identified recombinant Wolbachia strain that only occurs in hybrid grasshoppers. We propose the existence of two superimposed hybrid zones: one organized by the grasshoppers, which overlaps with a second, bacterial hybrid zone. The two hybrid zones counterbalance one another and have evolved together since the origin of the grasshopper's hybrid zone.

Clines on the seashore: The genomic architecture underlying rapid divergence in the face of gene flow

Adaptive divergence and speciation may happen despite opposition by gene flow. Identifying the genomic basis underlying divergence with gene flow is a major task in evolutionary genomics. Most approaches (e.g., outlier scans) focus on genomic regions of high differentiation. However, not all genomic architectures potentially underlying divergence are expected to show extreme differentiation. Here, we develop an approach that combines hybrid zone analysis (i.e., focuses on spatial patterns of allele frequency change) with system-specific simulations to identify loci inconsistent with neutral evolution. We apply this to a genome-wide SNP set from an ideally suited study organism, the intertidal snail Littorina saxatilis, which shows primary divergence between ecotypes associated with different shore habitats. We detect many SNPs with clinal patterns, most of which are consistent with neutrality. Among non-neutral SNPs, most are located within three large putative inversions differentiating ecotypes. Many non-neutral SNPs show relatively low levels of differentiation. We discuss potential reasons for this pattern, including loose linkage to selected variants, polygenic adaptation and a component of balancing selection within populations (which may be expected for inversions). Our work is in line with theory predicting a role for inversions in divergence, and emphasizes that genomic regions contributing to divergence may not always be accessible with methods purely based on allele frequency differences. These conclusions call for approaches that take spatial patterns of allele frequency change into account in other systems.

Co-occurrence among three divergent plant-castrating fungi in the same Silene host species

The competitive exclusion principle postulates that different species can only coexist in sympatry if they occupy distinct ecological niches. The goal of this study was to understand the geographical distribution of three species of Microbotryum anther-smut fungi that are distantly related but infect the same host plants, the sister species Silene vulgaris and S. uniflora, in Western Europe. We used microsatellite markers to investigate pathogen distribution in relation to host specialization and ecological factors. Microbotryum violaceo-irregulare was only found on S. vulgaris at high elevations in the Alps. Microbotryum lagerheimii could be subdivided into two genetically differentiated clusters, one on S. uniflora in the UK and the second on S. vulgaris in the Alps and Pyrenees. The most abundant pathogen species, M. silenes-inflatae, could be subdivided into four genetic clusters, co-occurring in the Alps, the UK and the Pyrenees, and was found on both S. vulgaris and S. uniflora. All three fungal species had high levels of homozygosity, in agreement with the selfing mating system generally observed in anther-smut fungi. The three pathogen species and genetic clusters had large range overlaps, but occurred at sites with different elevations, temperatures and precipitation levels. The three Microbotryum species thus do not appear to be maintained by host specialization or geographic allopatry, but instead may occupy different ecological niches in terms of environmental conditions.

Molecular and ecological signatures of an expanding hybrid zone

Many species are currently changing their distributions and subsequently form sympatric zones with hybridization between formerly allopatric species as one possible consequence. The damselfly Ischnura elegans has recently expanded south into the range of its ecologically and morphologically similar sister species Ischnura graellsii. Molecular work shows ongoing introgression between these species, but the extent to which this species mixing is modulated by ecological niche use is not known. Here, we (1) conduct a detailed population genetic analysis based on molecular markers and (2) model the ecological niche use of both species in allopatric and sympatric regions. Population genetic analyses showed chronic introgression between I. elegans and I. graellsii across a wide part of Spain, and admixture analysis corroborated this, showing that the majority of I. elegans from the sympatric zone could not be assigned to either the I. elegans or I. graellsii species cluster. Niche modeling demonstrated that I. elegans has modified its environmental niche following hybridization and genetic introgression with I. graellsii, making niche space of introgressed I. elegans populations more similar to I. graellsii. Taken together, this corroborates the view that adaptive introgression has moved genes from I. graellsii into I. elegans and that this process is enabling Spanish I. elegans to occupy a novel niche, further facilitating its expansion. Our results add to the growing evidence that hybridization can play an important and creative role in the adaptive evolution of animals.

POPULATION CAGE EXPERIMENTS WITH A VERTEBRATE: THE TEMPORAL DEMOGRAPHY AND CYTONUCLEAR GENETICS OF HYBRIDIZATION IN GAMBUSIA FISHES

The dynamics of mitochondrial and multilocus nuclear genotypic frequencies were monitored for 2 yr in experimental populations established with equal numbers of two poeciliid fishes (Gambusia affinis and Gambusia holbrooki) that hybridize naturally in the southeastern United States. In replicated "small-pool" populations (experiment I), 1018 sampled individuals at six time periods revealed an initial flush of hybridization, followed by a rapid decline in frequencies of G. affinis nuclear and mitochondrial alleles over 64 wk. Decay of gametic and cytonuclear disequilibria differed from expectations under random mating as well as under a model of assortative mating involving empirically estimated mating propensities. In two replicate "large-pond" populations (experiment II), 841 sampled individuals across four reproductive cohorts revealed lower initial frequencies of F₁ hybrids than in experiment I, but again G. holbrooki alleles achieved high frequencies over four generations (72 wk). Thus, evolution within experimental Gambusia hybrid populations can be extremely rapid, resulting in consistent loss of G. affinis nuclear and cytoplasmic alleles. Concordance in results between experiments and across genetic markers suggests strong directional selection favoring G. holbrooki genotypes. Results are interpreted in light of previous reports of genotype-specific differences in life-history traits, reproductive ecology, patterns of recruitment, and size-specific mortality, and in the context of patterns of introgression previously studied indirectly from spatial observations on cytonuclear genotypes in natural Gambusia populations.

HYBRID ZONE DYNAMICS ARE INFLUENCED BY GENOTYPE-SPECIFIC VARIATION IN LIFE-HISTORY TRAITS: EXPERIMENTAL EVIDENCE FROM HYBRIDIZING GAMBUSIA SPECIES

Results from two experiments are presented that contrast differences in life-history traits and population dynamics between two species of live bearing fishes (Gambusia affinis and G. holbrooki) that hybridize across portions of the southeastern United States. Progeny from parental holbrooki and holbrooki-affinis F₁ crosses exhibited larger lengths at birth, at 15 days, and matured earlier, and at larger size than did progeny from parental affinis and affinis-holbrooki F₁ crosses. Comparisons of experimental populations of affinis, holbrooki, and mixed (affinis + holbrooki) species composition followed over two years revealed that affinis populations consistently exhibited smaller population size, lower carrying capacity, lower recruitment, and larger over-winter mortality than did holbrooki or mixed populations. Evidence for density-dependent reductions in fecundity and concomitant increases in juvenile mortality rates were observed in all populations, but were most pronounced for affinis populations. Genotype-specific differences in life-history traits appear to confer differential advantage to offspring of parental holbrooki origin and F, progeny of holbrooki maternal parentage given the resource availability and the age structure and densities experienced during these experiments. Results have direct implications regarding the rate and direction of evolution within hybrid zones formed by these two species.

GENETICS OF PHYSIOLOGICAL DIFFERENTIATION WITHIN THE MARINE MUSSEL GENUS MYTILUS

Two divergent taxa in the marine mussel genus Mytilus are largely isolated geographically and are routinely exposed to distinctly different thermal environments. We tested the hypothesis that the two taxa are physiologically differentiated with respect to temperature and examined the evolved adaptations allowing one of the taxa to exploit habitats where warm-temperate conditions prevail for prolonged periods. We first analyzed the physiological response to high temperature of mussels collected from a hybrid population containing members of both pure taxa, F, hybrids, and a variety of introgressed genotypes. The experimental temperature of 23°C was chosen to be permissive to the taxon that occurs in warm-temperate regions (Mytilus galloprovincialis) and restrictive to the cold-water taxon (Mytilus edulis). The results show that the two taxa are physiologically differentiated. Under the experimental conditions, M. galloprovincialis exhibited a threefold higher feeding rate and a slightly elevated metabolic rate compared with M. edulis. These differences did not result in a significant difference in net energy balance between the two taxa, probably because of an interaction between physiological response and food availability. However, M. galloprovincialis grew significantly faster in the field, indicating that the physiological differences observed in the laboratory also occur in nature. Numerous introgressed genotypes provided the opportunity to test for cosegregation between the physiological differences and four highly differentiated genetic markers. Two of the markers (esterase and octopine dehydrogenase) cosegregate with variation in feeding rate and shell growth and explained most of the physiological differences observed between taxa. A strong concordance existed between these two loci, suggesting that they may be linked and may mark segregation of the same linkage group. The results suggest that the physiological differentiation between these taxa may be controlled by a few genes (perhaps only one) each with large effect.

NEUTRAL GENE FLOW ACROSS SINGLE LOCUS CLINES

We study the strength of the genetic barrier to neutral gene flow in a general one-locus, two-allele model that includes viability selection as well as fertility selection and premating isolation. We have separately considered adult migration and pollen migration. Our theoretical results suggest that selection against hybrid formation in the form of fertility selection or assortative mating is more effective in preventing introgression of neutral alleles than is reduction in hybrid viability. We argue that in experimental studies of introgression testing of F₁ hybrids is as important as testing of parental forms. To illustrate the utility of this approach, we estimate the strength of the barrier to neutral gene flow between Piriqueta caroliniana and P. viridis, between Iris hexagona and I. fulva, and between Chorthippus brunneus and C. biguttulus.

NARROW HYBRID ZONE BETWEEN TWO SUBSPECIES OF BIG SAGEBRUSH (ARTEMISIA TRIDENTATA: ASTERACEAE). IV. RECIPROCAL TRANSPLANT EXPERIMENTS

Does endogenous or exogenous selection stabilize the big sagebrush (Artemisia tridentata) hybrid zone? After two years of study, our reciprocal transplant experiments showed significant genotype by environment interactions for a number of fitness components, including germination, growth, and reproduction. Hybrids were the most fit within the hybrid garden. In the parental gardens, the native parental taxon was more fit than either the alien parental or hybrids. These results are consistent with the bounded hybrid superiority model, which assumes exogenous selection, but are clearly at odds with the dynamic equilibrium model, which assumes endogenous selection and universal hybrid unfitness.

Genetic structure of Micromeria (Lamiaceae) in Tenerife, the imprint of geological history and hybridization on within-island diversification

Geological history of oceanic islands can have a profound effect on the evolutionary history of insular flora, especially in complex islands such as Tenerife in the Canary Islands. Tenerife results from the secondary connection of three paleo‐islands by a central volcano, and other geological events that further shaped it. This geological history has been shown to influence the phylogenetic history of several taxa, including genus Micromeria (Lamiaceae). Screening 15 microsatellite markers in 289 individuals representing the eight species of Micromeria present in Tenerife, this study aims to assess the genetic diversity and structure of these species and its relation with the geological events on the island. In addition, we evaluate the extent of hybridization among species and discuss its influence on the speciation process. We found that the species restricted to the paleo‐islands present lower levels of genetic diversity but the highest levels of genetic differentiation suggesting that their ranges might have contracted over time. The two most widespread species in the island, M. hyssopifolia and M. varia, present the highest genetic diversity levels and a genetic structure that seems correlated with the geological composition of the island. Samples from M. hyssopifolia from the oldest paleo‐island, Adeje, appear as distinct while samples from M. varia segregate into two main clusters corresponding to the paleo‐islands of Anaga and Teno. Evidence of hybridization and intraspecific migration between species was found. We argue that species boundaries would be retained despite hybridization in response to the habitat's specific conditions causing postzygotic isolation and preserving morphological differentiation.

Is cascade reinforcement likely when sympatric and allopatric populations exchange migrants?

When partially reproductively isolated species come back into secondary contact, these taxa may diverge in mating preferences and sexual cues to avoid maladaptive hybridization, a process known as reinforcement. This phenomenon often leads to reproductive character displacement (RCD) between sympatric and allopatric populations of reinforcing species that differ in their exposure to hybridization. Recent discussions have reinvigorated the idea that RCD may give rise to additional speciation between conspecific sympatric and allopatric populations, dubbing the concept "cascade reinforcement." Despite some empirical studies supporting cascade reinforcement, we still know very little about the conditions for its evolution. In the present article, we address this question by developing an individual-based population genetic model that explicitly simulates cascade reinforcement when one of the hybridizing species is split into sympatric and allopatric populations. Our results show that when sympatric and allopatric populations reside in the same environment and only differ in their exposure to maladaptive hybridization, migration between them generally inhibits the evolution of cascade by spreading the reinforcement alleles from sympatry into allopatry and erasing RCD. Under these conditions, cascade reinforcement only evolved when migration rate between sympatric and allopatric populations was very low. This indicates that stabilizing sexual selection in allopatry is generally ineffective in preventing the spread of reinforcement alleles. Only when sympatric and allopatric populations experienced divergent ecological selection did cascade reinforcement evolve in the presence of substantial migration. These predictions clarify the conditions for cascade reinforcement and facilitate our understanding of existing cases in nature.

Evolutionary consequences of climate-induced range shifts in insects

Range shifts can rapidly create new areas of geographic overlap between formerly allopatric taxa and evidence is accumulating that this can affect species persistence. We review the emerging literature on the short- and long-term consequences of these geographic range shifts. Specifically, we focus on the evolutionary consequences of novel species interactions in newly created sympatric areas by describing the potential (i) short-term processes acting on reproductive barriers between species and (ii) long-term consequences of range shifts on the stability of hybrid zones, introgression and ultimately speciation and extinction rates. Subsequently, we (iii) review the empirical literature on insects to evaluate which processes have been studied, and (iv) outline some areas that deserve increased attention in the future, namely the genomics of hybridisation and introgression, our ability to forecast range shifts and the impending threat from insect vectors and pests on biodiversity, human health and crop production. Our review shows that species interactions in de novo sympatric areas can be manifold, sometimes increasing and sometimes decreasing species diversity. A key issue that emerges is that climate-induced hybridisations in insects are much more widespread than anticipated and that rising temperatures and increased anthropogenic disturbances are accelerating the process of species mixing. The existing evidence only shows the tip of the iceberg and we are likely to see many more cases of species mixing following range shifts in the near future.

Pre-zygotic isolation in the macroalgal genus Fucus from four contact zones spanning 100-10 000 years: a tale of reinforcement?

Hybrid zones provide an ideal natural experiment to study the selective forces driving evolution of reproductive barriers and speciation. If hybrid offspring are less fit than the parental species, pre-zygotic isolating barriers can evolve and strengthen in response to selection against the hybrids (reinforcement). Four contact zones between the intertidal macroalgae Fucus serratus (Fs) and Fucus distichus (Fd), characterized by varying times of sympatry and order of species introduction provide an opportunity to investigate reinforcement. We examined patterns of hybridization and reproductive isolation between Fs and Fd in: (i) northern Norway (consisting of two natural sites, 10 000 years old), (ii) the Kattegat near Denmark (Fd introduced, nineteenth century) and (iii) Iceland (Fs introduced, nineteenth century). Using 10 microsatellites and chloroplast DNA, we showed that hybridization and introgression decreased with increasing duration of sympatry. The two younger contact zones revealed 13 and 24% hybrids and several F 1 individuals, in contrast to the older contact zone with 2-3% hybrids and an absence of F 1s. Cross-fertilization experiments revealed that the reduction in hybridization in the oldest zone is consistent with increased gametic incompatibility.

Creating new evolutionary pathways through bioinvasion: the population genetics of brushtail possums in New Zealand

Rapid increases in global trade and human movement have created novel mixtures of organisms bringing with them the potential to rapidly accelerate the evolution of new forms. The common brushtail possum (Trichosurus vulpecula), introduced into New Zealand from Australia in the 19th century, is one such species having been sourced from multiple populations in its native range. Here, we combine microsatellite DNA- and GIS-based spatial data to show that T. vulpecula originating from at least two different Australian locations exhibit a population structure that is commensurate with their introduction history and which cannot be explained by landscape features alone. Most importantly, we identify a hybrid zone between the two subspecies which appears to function as a barrier to dispersal. When combined with previous genetic, morphological and captive studies, our data suggest that assortative mating between the two subspecies may operate at a behavioural or species recognition level rather than through fertilization, genetic incompatibility or developmental inhibition. Nevertheless, hybridization between the two subspecies of possum clearly occurs, creating the opportunity for novel genetic combinations that would not occur in their natural ranges and which is especially likely given that multiple contact zones occur in New Zealand. This discovery has implications for wildlife management in New Zealand because multiple contact zones are likely to influence the dispersal patterns of possums and because differential susceptibility to baiting with sodium fluoroacetate between possums of different origins may promote novel genetic forms.

Wolbachia effects in natural populations of Chorthippus parallelus from the Pyrenean hybrid zone

We evaluate for the first time the effect of Wolbachia infection, involving two different supergroups, on the structure and dynamics of the hybrid zone between two subspecies of Chorthippus parallelus (Orthoptera) in the Pyrenees. Wolbachia infection showed no effects on female fecundity or a slight increment in females infected by F supergroup, although in the last case it has to be well established. Cytoplasmic incompatibility (CI) is confirmed in crosses carried out in the field between individuals from a natural hybrid population. This CI, registered as the relative reduction in embryo production (sh ), was of sh = 0.355 and sh = 0.286 in unidirectional crosses involving B and F supergroups, respectively. CI also occurred in bidirectional crosses (sh = 0.147) but with a weaker intensity. The transmission rates of the two Wolbachia strains (B and F) were estimated by the optimization of a theoretical model to reach the infection frequencies observed in certain population. To fit this scenario, both supergroups should present transmission rates close to 1. Further, we have simulated the infection dynamics, and hence, the capacity of Wolbachia to structure the population of the host insects and to affect to reproduction and genetic introgression in the hybrid zone. This represents a first example of the influence of Wolbachia in an insect natural hybrid zone.

Segregation of male-sterility alleles across a species boundary

Hybrid zones may serve as bridges permitting gene flow between species, including alleles influencing the evolution of breeding systems. Using greenhouse crosses, we assessed the likelihood that a hybrid zone could serve as a conduit for transfer of nuclear male-sterility alleles between a gynodioecious species and a hermaphroditic species with very rare females in some populations. Segregation patterns in progeny of crosses between rare females of hermaphroditic Schiedea menziesii and hermaphroditic plants of gynodioecious Schiedea salicaria heterozygous at the male-sterility locus, and between female S. salicaria and hermaphroditic plants from the hybrid zone, were used to determine whether male-sterility was controlled at the same locus in the parental species and the hybrid zone. Segregations of females and hermaphrodites in approximately equal ratios from many of the crosses indicate that the same nuclear male-sterility allele occurs in the parent species and the hybrid zone. These rare male-sterility alleles in S. menziesii may result from gene flow from S. salicaria through the hybrid zone, presumably facilitated by wind pollination in S. salicaria. Alternatively, rare male-sterility alleles might result from a reversal from gynodioecy to hermaphroditism in S. menziesii, or possibly de novo evolution of male sterility. Phylogenetic analysis indicates that some species of Schiedea have probably evolved separate sexes independently, but not in the lineage containing S. salicaria and S. menziesii. High levels of selfing and expression of strong inbreeding depression in S. menziesii, which together should favour females in populations, argue against a reversal from gynodioecy to hermaphroditism in S. menziesii.

Alternative forms for genomic clines

Understanding factors regulating hybrid fitness and gene exchange is a major research challenge for evolutionary biology. Genomic cline analysis has been used to evaluate alternative patterns of introgression, but only two models have been used widely and the approach has generally lacked a hypothesis testing framework for distinguishing effects of selection and drift. I propose two alternative cline models, implement multivariate outlier detection to identify markers associated with hybrid fitness, and simulate hybrid zone dynamics to evaluate the signatures of different modes of selection. Analysis of simulated data shows that previous approaches are prone to false positives (multinomial regression) or relatively insensitive to outlier loci affected by selection (Barton's concordance). The new, theory-based logit-logistic cline model is generally best at detecting loci affecting hybrid fitness. Although some generalizations can be made about different modes of selection, there is no one-to-one correspondence between pattern and process. These new methods will enhance our ability to extract important information about the genetics of reproductive isolation and hybrid fitness. However, much remains to be done to relate statistical patterns to particular evolutionary processes. The methods described here are implemented in a freely available package “HIest” for the R statistical software (CRAN; http://cran.r-project.org/).

Patterns of hybridization and asymmetrical gene flow in hybrid zones of the rare Eucalyptus aggregata and common E. rubida

The patterns of hybridization and asymmetrical gene flow among species are important for understanding the processes that maintain distinct species. We examined the potential for asymmetrical gene flow in sympatric populations of Eucalyptus aggregata and Eucalyptus rubida, both long-lived trees of southern Australia. A total of 421 adults from three hybrid zones were genotyped with six microsatellite markers. We used genealogical assignments, admixture analysis and analyses of spatial genetic structure and spatial distribution of individuals, to assess patterns of interspecific gene flow within populations. A high number of admixed individuals were detected (13.9-40% of individuals), with hybrid populations consisting of F(1) and F(2) hybrids and backcrosses in both parental directions. Across the three sites, admixture proportions were skewed towards the E. aggregata genetic cluster (x=0.56-0.65), indicating that backcrossing towards E. aggregata is more frequent. Estimates of long-term migration rates also indicate asymmetric gene flow, with higher migration rates from E. aggregata to hybrids compared with E. rubida. Taken together, these results indicate a greater genetic input from E. aggregata into the hybrid populations. This asymmetry probably reflects differences in style lengths (E. rubida: ~7 mm, E. aggregata: ~4 mm), which can prevent pollen tubes of smaller-flowered species from fertilizing larger-flowered species. However, analyses of fine-scale genetic structure suggest that localized seed dispersal (<40 m) and greater clustering between hybrid and E. aggregata individuals may also contribute to directional gene flow. Our study highlights that floral traits and the spatial distributions of individuals can be useful predictors of the directionality of interspecific gene flow in plant populations.

Mycorrhizal associations and reproductive isolation in three closely related Orchis species

BACKGROUND AND AIMS

The maintenance of species boundaries in sympatric populations of closely related species requires some kind of reproductive isolation that limits gene flow among species and/or prevents the production of viable progeny. Because in orchids mycorrhizal fungi are needed for seed germination and subsequent seedling establishment, orchid-mycorrhizal associations may be involved in acting as a post-mating barrier.

METHODS

We investigated the strength of post-mating barriers up to the seed germination stage acting between three closely related Orchis species (Orchis anthropophora, O. militaris and O. purpurea) and studied the role of mycorrhizal fungi in hybridization by burying seed packets of pure and hybrid seeds. After retrieval and assessment of seed germination, the fungi associating with protocorms originating from hybrid and pure seeds were determined and compared with those associating with adult individuals using DNA array technology.

RESULTS

Whereas pre-zygotic post-mating barriers were rather weak in most crosses, post-zygotic post-mating barriers were stronger, particularly when O. purpurea was crossed with O. anthropophora. Germination trials in the field showed that seed germination percentages of hybrid seeds were in most cases lower than those originating from pure crosses. In all species pair combinations, total post-mating reproductive isolation was asymmetric. Protocorms associated with a smaller range of fungal symbionts than adult plants, but there was considerable overlap in mycorrhizal associations between protocorms and their respective parents.

CONCLUSIONS

Our results suggest that mycorrhizal associations contribute little to reproductive isolation. Pre-mating barriers are probably the main factors determining hybridization rates between the investigated species.

SPECIATION IN MAMMALS AND THE GENETIC SPECIES CONCEPT

We define a genetic species as a group of genetically compatible interbreeding natural populations that is genetically isolated from other such groups. This focus on genetic isolation rather than reproductive isolation distinguishes the Genetic Species Concept from the Biological Species Concept. Recognition of species that are genetically isolated (but not reproductively isolated) results in an enhanced understanding of biodiversity and the nature of speciation as well as speciation-based issues and evolution of mammals. We review criteria and methods for recognizing species of mammals and explore a theoretical scenario, the Bateson-Dobzhansky-Muller (BDM) model, for understanding and predicting genetic diversity and speciation in mammals. If the BDM model is operating in mammals, then genetically defined phylogroups would be predicted to occur within species defined by morphology, and phylogroups experiencing stabilizing selection will evolve genetic isolation without concomitant morphological diversification. Such species will be undetectable using classical skin and skull morphology (Morphological Species Concept). Using cytochrome-b data from sister species of mammals recognized by classical morphological studies, we estimated the number of phylogroups that exist within mammalian species and hypothesize that there will be >2,000 currently unrecognized species of mammals. Such an underestimation significantly affects conclusions on the nature of speciation in mammals, barriers associated with evolution of genetic isolation, estimates of biodiversity, design of conservation initiatives, zoonoses, and so on. A paradigm shift relative to this and other speciation-based issues will be needed. Data that will be effective in detecting these "morphologically cryptic genetic species" are genetic, especially DNA-sequence data. Application of the Genetic Species Concept uses genetic data from mitochondrial and nuclear genomes to identify species and species boundaries, the extent to which the integrity of the gene pool is protected, nature of hybridization (if present), and introgression. Genetic data are unique in understanding species because the use of genetic data 1) can quantify genetic divergence from different aspects of the genome (mitochondrial and nuclear genes, protein coding genes, regulatory genes, mobile DNA, microsatellites, chromosomal rearrangements, heterochromatin, etc.); 2) can provide divergence values that increase with time, providing an estimate of time since divergence; 3) can provide a population genetics perspective; 4) is less subject to convergence and parallelism relative to other sets of characters; 5) can identify monophyly, sister taxa, and presence or absence of introgression; and 6) can accurately identify hybrid individuals (kinship and source of hybrid individuals, F(1)s, backcrosses, direction of hybridization, and in concert with other data identify which hybrids are sterile or fertile). The proposed definition of the Genetic Species Concept is more compatible with a description of biodiversity of mammals than is "reproductively isolated species." Genetic profiles of mammalian species will result in a genetic description of species and mammalian diversity, and such studies are being accelerated by technological advances that reduce cost and increase speed and efficiency of generating genetic data. We propose that this genetic revolution remain museum- and voucher specimen-based and that new names are based on a holotype (including associated tissues) deposited in an accredited museum.

Hybrid zones as a tool for identifying adaptive genetic variation in outbreeding forest trees: lessons from wild annual sunflowers (Helianthus spp.)

The identification and study of adaptively important genes in forest trees represents a formidable challenge because of their long generation spans. In annual or perennial herbs, formal genetic studies can be employed to identify the quantitative trait loci (QTLs) and/or candidate genes that underlie important traits, and the segregating populations can be transplanted into natural populations to measure the strength and direction of selection. However, the application of these methods to forest trees is difficult, because the creation of appropriate genetic material is extremely time-consuming in long-lived, woody plants, and lifetime fitness estimates are difficult or impossible to obtain. Although QTL mapping should in principle be feasible in wild intraspecific populations (as an alternative to artificial crosses), this approach is less likely to be successful in trees because LD (linkage disequilibrium) will decay quickly in large outbreeding plant populations. Within the present paper, we discuss a modified approach based on natural hybrid zones. We describe the use of wild annual sunflowers (Helianthus spp.) as a model for exploring the hybrid zone approach. Transplanted experimental hybrids allowed us to assess the adaptive value of individual chromosomal blocks in nature, and data on natural Helianthus hybrids suggest that similar approaches are possible in natural hybrid zones. Our results allowed us to test the role of hybridization in the origin of ecological divergence in wild sunflowers. In addition, they have practical implications for identifying adaptively important genes or QTLs in trees. This is exemplified by three temperate forest taxa, Populus (poplars, aspens, cottonwoods), Fraxinus (ash), and Quercus (oak). All three are diploid and important genomic tools are under development. Moreover, all three offer extensive hybrid zones whose likely age can be inferred from fossil data. Age data enables estimates of the size and frequency of chromosomal blocks in hybrids, thereby providing guidance in designing marker-based experiments. We predict that natural hybrid zones will be valuable tools for identifying the QTLs and/or candidate genes responsible for adaptive traits in forest trees.