Ecological niche contributes to the persistence of the Western x Glaucous-winged Gull hybrid zone

Hybrid zones occur in nature when populations with limited reproductive barriers overlap in space. Many hybrid zones persist over time, and different models have been proposed to explain how selection can maintain hybrid zone stability. More empirical studies are needed to elucidate the role of ecological adaptation in maintaining stable hybrid zones. Here, we investigated the role of exogenous factors in maintaining a hybrid zone between western gulls (Larus occidentalis) and glaucous-winged gulls (L. glaucescens). We used ecological niche models (ENMs) and niche similarity tests to quantify and examine the ecological niches of western gulls, glaucous-winged gulls, and their hybrids. We found evidence of niche divergence between all three groups. Our results best support the bounded superiority model, providing further evidence that exogenous selection favoring hybrids may be an important factor in maintaining this stable hybrid zone.

Variance effective population size is affected by census size in sub-structured populations

Measurement of allele frequency shifts between temporally spaced samples has long been used for assessment of effective population size (N_e ), and this 'temporal method' provides estimates of N_e referred to as variance effective size (N_eV ). We show that N_eV of a local population that belongs to a sub-structured population (a metapopulation) is determined not only by genetic drift and migration rate (m), but also by the census size (N_c ). The realized N_eV of a local population can either increase or decrease with increasing m, depending on the relationship between N_e and N_c in isolation. This is shown by explicit mathematical expressions for the factors affecting N_eV derived for an island model of migration. We verify analytical results using high-resolution computer simulations, and show that the phenomenon is not restricted to the island model migration pattern. The effect of N_c on the realized N_eV of a local subpopulation is most pronounced at high migration rates. We show that N_c only affects local N_eV , whereas N_eV for the metapopulation as a whole, inbreeding (N_eI ), and linkage disequilibrium (N_eLD ) effective size are all independent of N_c . Our results provide a possible explanation to the large variation of N_e /N_c ratios reported in the literature, where N_e is frequently estimated by N_eV . They are also important for the interpretation of empirical N_e estimates in genetic management where local N_eV is often used as a substitute for inbreeding effective size, and we suggest an increased focus on metapopulation N_eV as a proxy for N_eI .

Interspecific introgression of MHC genes in Triturus newts: Evidence from multiple contact zones

The major histocompatibility complex (MHC) genes are central to the adaptive immune response in vertebrates. Selection generally maintains high MHC variation because the spectrum of recognized pathogens depends on MHC polymorphism. Novel alleles favoured by selection originate by interallelic recombination or de novo mutations but may also be acquired by introgression from related species. However, the extent and prevalence of MHC introgression remain an open question. In this study, we tested for MHC introgression in six hybrid zones formed by six Triturus newt species. We sequenced and genotyped the polymorphic second exons of the MHC class I and II genes and compared their interspecific similarity at various distances from the centre of the hybrid zone. We found evidence for introgression of both MHC classes in the majority of examined hybrid zones, with support for a more substantial class I introgression. Furthermore, the overall MHC allele sharing outside of hybrid zones was elevated between pairs of Triturus species with abutting ranges, regardless of the phylogenetic distance between them. No effect of past hybrid zone movement on MHC allele sharing was found. Finally, using previously published genome-wide data, we demonstrated that MHC introgression was more extensive than genome-wide introgression, supporting its adaptive potential. Our study thus provides evidence for the prevalence of MHC introgression across multiple Triturus hybrid zones, indicating that MHC introgression between divergent hybridizing species may be widespread and adaptive.

Absence of heterosis in hybrid crested newts

Relationships between phylogenetic relatedness, hybrid zone spatial structure, the amount of interspecific gene flow and population demography were investigated, with the newt genus Triturus as a model system. In earlier work, a bimodal hybrid zone of two distantly related species combined low interspecific gene flow with hybrid sterility and heterosis was documented. Apart from that, a suite of unimodal hybrid zones in closely related Triturus showed more or less extensive introgressive hybridization with no evidence for heterosis. We here report on population demography and interspecific gene flow in two Triturus species (T. macedonicus and T. ivanbureschi in Serbia). These are two that are moderately related, engage in a heterogeneous uni-/bimodal hybrid zone and hence represent an intermediate situation. This study used 13 diagnostic nuclear genetic markers in a population at the species contact zone. This showed that all individuals were hybrids, with no parentals detected. Age, size and longevity and the estimated growth curves are not exceeding that of the parental species, so that we conclude the absence of heterosis in T. macedonicus-T. ivanbureschi. Observations across the genus support the hypothesis that fertile hybrids allocate resources to reproduction and infertile hybrids allocate resources to growth. Several Triturus species hybrid zones not yet studied allow the testing of this hypothesis.

THE DYNAMICS OF TWO HYBRID ZONES IN APPALACHIAN SALAMANDERS OF THE GENUS PLETHODON

Two zones of intergradation between populations of Plethodon have been studied for 18 and 20 years, respectively. The data consist of systematic scores of colors, made at least twice annually. Near Heintooga Overlook in the Balsam Mountains (Great Smoky Mountains National Park), the salamanders' cheeks are gray. Proceeding north toward the Smokies, there is increasing frequency and intensity of red color at two, four, and six miles. There has been no change in the scores at any location. The width of the zone and our failure to detect any change can be explained by assuming neutrality of the character and random diffusion during the probable time since contact between the two intergrading forms, which most likely took place after the Hypsithermal Interval, 8,000-5,000 BP. At Coweeta Hydrologic Laboratory in the Nantahala Mountains, Plethodon jordani and P. glutinosus hybridize at intermediate elevations. The lateral white spots of glutinosus decrease and the red on the legs of jordani increases with elevation from 685 m to 1,052 m. At the higher elevation, the proportion of animals scored as "pure" jordani declined significantly from 1974 to 1990, an indication that the hybrid zone is spreading upward. The rate of spread is too great to be explained by random diffusion, so selection for glutinosus characters is the best explanation. The rate of spread of the hybrid zone indicates that hybridization began 60-65 years ago, at the end of the time of intense timbering. Such human disturbances have caused hybridization in other organisms.

Stabilization of a salamander moving hybrid zone

When related species meet upon postglacial range expansion, hybrid zones are frequently formed. Theory predicts that such zones may move over the landscape until equilibrium conditions are reached. One hybrid zone observed to be moving in historical times (1950–1979) is that of the pond‐breeding salamanders Triturus cristatus and Triturus marmoratus in western France. We identified the ecological correlates of the species hybrid zone as elevation, forestation, and hedgerows favoring the more terrestrial T. marmoratus and pond density favoring the more aquatic T. cristatus. The past movement of the zone of ca. 30 km over three decades has probably been driven by the drastic postwar reduction of the “bocage” hedgerow landscape, favoring T. cristatus over T. marmoratus. No further hybrid zone movement was observed from 1979 to the present. To explain the changing dynamics of the hybrid zone, we propose that it stalled, either because an equilibrium was found at an altitude of ca. 140 m a.s.l. or due to pond loss and decreased population densities. While we cannot rule out the former explanation, we found support for the latter. Under agricultural intensification, ponds in the study area are lost at an unprecedented rate of 5.5% per year, so that remaining Triturus populations are increasingly isolated, hampering dispersal and further hybrid zone movement.

Concordant morphological and molecular clines in a contact zone of the Common and Spined toad (Bufo bufo and B. spinosus) in the northwest of France

BACKGROUND

Hybrid zones are regions where individuals of two species meet and produce hybrid progeny, and are often regarded as natural laboratories to understand the process of species formation. Two microevolutionary processes can take place in hybrid zones, with opposing effects on population differentiation. Hybridization tends to produce genetic homogenization, reducing species differences, whereas the presence of mechanisms of reproductive isolation result in barriers to gene flow, maintaining or increasing differences between taxa.

RESULTS

Here we study a contact zone between two hybridizing toad species, Bufo bufo and B. spinosus, through a combination of molecular (12 polymorphic microsatellites, four nuclear and two mitochondrial SNP markers) and morphological data in a transect in the northwest of France. The results show largely concordant clines across markers, defining a narrow hybrid zone of ca. 30 km wide. Most hybrids in the centre of the contact zone are classified as F₂ or backcrossed individuals, with no individuals assigned to the F₁ hybrid class.

CONCLUSIONS

We discuss the implications of these results for our understanding of the evolutionary history of these species. We anticipate that the toad contact zone here described will become an important asset in the study of hybrid zone dynamics and evolutionary biology because of its easy access and the abundance of the species involved.

Homeotic transformations and number changes in the vertebral column of Triturus newts

We explored intraspecific variation in vertebral formulae, more specifically the variation in the number of thoracic vertebrae and frequencies of transitional sacral vertebrae in Triturus newts (Caudata: Salamandridae). Within salamandrid salamanders this monophyletic group shows the highest disparity in the number of thoracic vertebrae and considerable intraspecific variation in the number of thoracic vertebrae. Triturus species also differ in their ecological preferences, from predominantly terrestrial to largely aquatic. Following Geoffroy St. Hilaire's and Darwin's rule which states that structures with a large number of serially homologous repetitive elements are more variable than structures with smaller numbers, we hypothesized that the variation in vertebral formulae increases in more elongated species with a larger number of thoracic vertebrae. We furthermore hypothesized that the frequency of transitional vertebrae will be correlated with the variation in the number of thoracic vertebrae within the species. We also investigated potential effects of species hybridization on the vertebral formula. The proportion of individuals with a number of thoracic vertebrae different from the modal number and the range of variation in number of vertebrae significantly increased in species with a larger number of thoracic vertebrae. Contrary to our expectation, the frequencies of transitional vertebrae were not correlated with frequencies of change in the complete vertebrae number. The frequency of transitional sacral vertebra in hybrids did not significantly differ from that of the parental species. Such a pattern could be a result of selection pressure against transitional vertebrae and/or a bias towards the development of full vertebrae numbers. Although our data indicate relaxed selection for vertebral count changes in more elongated, aquatic species, more data on different selective pressures in species with different numbers of vertebrae in the two contrasting, terrestrial and aquatic environments are needed to test for causality.

Hybrid zones: windows on climate change

Defining the impacts of anthropogenic climate change on biodiversity and species distributions is currently a high priority. Niche models focus primarily on predicted changes in abiotic factors; however, species interactions and adaptive evolution will impact the ability of species to persist in the face of changing climate. Our review focuses on the use of hybrid zones to monitor responses of species to contemporary climate change. Monitoring hybrid zones provides insight into how range boundaries shift in response to climate change by illuminating the combined effects of species interactions and physiological sensitivity. At the same time, the semipermeable nature of species boundaries allows us to document adaptive introgression of alleles associated with response to climate change.

Beam Brook revisited: a molecular study of a historically introduced non-native amphibian (Triturus carnifex) and its limited introgression into native UK Triturus cristatus populations

Human mediated introductions of non-native species can pose a major threat to global biodiversity on several accounts i.e. through competition, the introduction of novel pathogens, and genetic pollution. Where hybridization occurs between two closely related species the F1 offspring are usually phenotypically discernible whereas F2 hybrid, backcrossed or admixed individuals become more difficult to separate. At this point the utilization of molecular methods is required in conservation efforts to differentiate and manage populations. This study demonstrated how a possible threat of hybridization from an introduced non-native (T. carnifex) with a protected native newt species (T. cristatus) could be investigated with molecular tools, and examined the current extent of its genetic introgression over an 80 years period. The results confirmed that hybridization had taken place at the site of introduction (and continues to do so), and that historically limited local dispersal of both non-natives and/or hybrids had occurred sometime in the past. However, the data suggests that although dispersal of hybrids into a local satellite site may still be occuring, hybridization with native species appears limited.

Stable Heliconius butterfly hybrid zones are correlated with a local rainfall peak at the edge of the Amazon basin

Multilocus clines between Müllerian mimetic races of Heliconius butterflies provide a classic example of the maintenance of hybrid zones and their importance in speciation. Concordant hybrid zones in the mimics Heliconius erato and H. melpomene in northern Peru were carefully documented in the 1980s, and this prior work now permits a historical analysis of the movement or stasis of the zones. Previous work predicted that these zones might be moving toward the Andes due to selective asymmetry. Extensive deforestation and climate change might also be expected to affect the positions and widths of the hybrid zones. We show that the positions and shapes of these hybrid zones have instead remained remarkably stable between 1985 and 2012. The stability of this interaction strongly implicates continued selection, rather than neutral mixing following secondary contact. The stability of cline widths and strong linkage disequilibria (gametic correlation coefficients Rmax = 0.35-0.56 among unlinked loci) over 25 years suggest that mimetic selection pressures on each color pattern locus have remained approximately constant (s ≈ 0.13-0.40 per locus in both species). Exceptionally high levels of precipitation at the edge of the easternmost Andes may act as a population density trough for butterflies, trapping the hybrid zones at the foot of the mountains, and preventing movement. As such, our results falsify one prediction of the Pleistocene Refugium theory: That the ranges of divergent species or subspecies should be centered on regions characterized by maxima of rainfall, with hybrid zones falling in more arid regions between them.

Data concatenation, Bayesian concordance and coalescent-based analyses of the species tree for the rapid radiation of Triturus newts

The phylogenetic relationships for rapid species radiations are difficult to disentangle. Here we study one such case, namely the genus Triturus, which is composed of the marbled and crested newts. We analyze data for 38 genetic markers, positioned in 3-prime untranslated regions of protein-coding genes, obtained with 454 sequencing. Our dataset includes twenty Triturus newts and represents all nine species. Bayesian analysis of population structure allocates all individuals to their respective species. The branching patterns obtained by data concatenation, Bayesian concordance analysis and coalescent-based estimations of the species tree differ from one another. The data concatenation based species tree shows high branch support but branching order is considerably affected by allele choice in the case of heterozygotes in the concatenation process. Bayesian concordance analysis expresses the conflict between individual gene trees for part of the Triturus species tree as low concordance factors. The coalescent-based species tree is relatively similar to a previously published species tree based upon morphology and full mtDNA and any conflicting internal branches are not highly supported. Our findings reflect high gene tree discordance due to incomplete lineage sorting (possibly aggravated by hybridization) in combination with low information content of the markers employed (as can be expected for relatively recent species radiations). This case study highlights the complexity of resolving rapid radiations and we acknowledge that to convincingly resolve the Triturus species tree even more genes will have to be consulted.

The distribution of the crested and marbled newt species (Amphibia: Salamandridae: Triturus) – an addition to the New Atlas of Amphibians and Reptiles of Europe

In the recently published New Atlas of Amphibians and Reptiles of Europe (Sillero et al., 2014a), the distribution of the newt genusTrituruswas not resolved at the level of the species. The main reason for this was the lack of high quality distribution data from in and around the parapatric contact zones between species, where interspecific hybridization occurs. We are working extensively onTriturusand the (particularly genetic) data we have accumulated allow us to map the individualTriturusspecies at the appropriate scale. We here provide a database composed of distribution data for the individual species, at generally high resolution, particularly from in and around contact zones. Based on this database we produce maps at the 50 × 50 km UTM grid resolution as used in the new atlas and highlight those grid cells in which more than oneTriturusspecies occurs.

Multilocus phylogeography of the common lizard Zootoca vivipara at the Ibero-Pyrenean suture zone reveals lowland barriers and high-elevation introgression

BACKGROUND

The geographic distribution of evolutionary lineages and the patterns of gene flow upon secondary contact provide insight into the process of divergence and speciation. We explore the evolutionary history of the common lizard Zootoca vivipara (= Lacerta vivipara) in the Iberian Peninsula and test the role of the Pyrenees and the Cantabrian Mountains in restricting gene flow and driving lineage isolation and divergence. We also assess patterns of introgression among lineages upon secondary contact, and test for the role of high-elevation trans-mountain colonisations in explaining spatial patterns of genetic diversity. We use mtDNA sequence data and genome-wide AFLP loci to reconstruct phylogenetic relationships among lineages, and measure genetic structure.

RESULTS

The main genetic split in mtDNA corresponds generally to the French and Spanish sides of the Pyrenees as previously reported, in contrast to genome-wide AFLP data, which show a major division between NW Spain and the rest. Both types of markers support the existence of four distinct and geographically congruent genetic groups, which are consistent with major topographic barriers. Both datasets reveal the presence of three independent contact zones between lineages in the Pyrenean region, one in the Basque lowlands, one in the low-elevation mountains of the western Pyrenees, and one in the French side of the central Pyrenees. The latter shows genetic evidence of a recent, high-altitude trans-Pyrenean incursion from Spain into France.

CONCLUSIONS

The distribution and age of major lineages is consistent with a Pleistocene origin and a role for both the Pyrenees and the Cantabrian Mountains in driving isolation and differentiation of Z. vivipara lineages at large geographic scales. However, mountain ranges are not always effective barriers to dispersal, and have not prevented a recent high-elevation trans-Pyrenean incursion that has led to asymmetrical introgression among divergent lineages. Cytonuclear discordance in patterns of genetic structure and introgression at contact zones suggests selection may be involved at various scales. Suture zones are important areas for the study of lineage formation and speciation, and our results show that biogeographic barriers can yield markedly different phylogeographic patterns in different vertebrate and invertebrate taxa.

Postglacial species displacement in Triturus newts deduced from asymmetrically introgressed mitochondrial DNA and ecological niche models

BACKGROUND

If the geographical displacement of one species by another is accompanied by hybridization, mitochondrial DNA can introgress asymmetrically, from the outcompeted species into the invading species, over a large area. We explore this phenomenon using the two parapatric crested newt species, Triturus macedonicus and T. karelinii, distributed on the Balkan Peninsula in south-eastern Europe, as a model.

RESULTS

We first delimit a ca. 54,000 km(2) area in which T. macedonicus contains T. karelinii mitochondrial DNA. This introgression zone bisects the range of T. karelinii, cutting off a T. karelinii enclave. The high similarity of introgressed mitochondrial DNA haplotypes with those found in T. karelinii suggests a recent transfer across the species boundary. We then use ecological niche modeling to explore habitat suitability of the location of the present day introgression zone under current, mid-Holocene and Last Glacial Maximum conditions. This area was inhospitable during the Last Glacial Maximum for both species, but would have been habitable at the mid-Holocene. Since the mid-Holocene, habitat suitability generally increased for T. macedonicus, whereas it decreased for T. karelinii.

CONCLUSION

The presence of a T. karelinii enclave suggests that T. karelinii was the first to colonize the area where the present day introgression zone is positioned after the Last Glacial Maximum. Subsequently, we propose T. karelinii was outcompeted by T. macedonicus, which captured T. karelinii mitochondrial DNA via introgressive hybridization in the process. Ecological niche modeling suggests that this replacement was likely facilitated by a shift in climate since the mid-Holocene. We suggest that the northwestern part of the current introgression zone was probably never inhabited by T. karelinii itself, and that T. karelinii mitochondrial DNA spread there through T. macedonicus exclusively. Considering the spatial distribution of the introgressed mitochondrial DNA and the signal derived from ecological niche modeling, we do not favor the hypothesis that foreign mitochondrial DNA was pulled into the T. macedonicus range by natural selection.

Habitat association and seasonality in a mosaic and bimodal hybrid zone between Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae)

Understanding why some hybrid zones are bimodal and others unimodal can aid in identifying barriers to gene exchange following secondary contact. The hybrid zone between the grasshoppers Chorthippus brunneus and C. jacobsi contains a mix of allopatric parental populations and inter-mingled bimodal and unimodal sympatric populations, and provides an ideal system to examine the roles of local selection and gene flow between populations in maintaining bimodality. However, it is first necessary to confirm, over a larger spatial scale, previously identified associations between population composition and season and habitat. Here we use cline-fitting of one morphological and one song trait along two valley transects, and intervening mountains, to confirm previously identified habitat associations (mountain versus valley) and seasonal changes in population composition. As expected from previous findings of studies on a smaller spatial scale, C. jacobsi dominated mountain habitats and mixed populations dominated valleys, and C. brunneus became more prevalent in August. Controlling for habitat and incorporating into the analysis seasonal changes in cline parameters and the standard errors of parental trait values revealed wider clines than previous studies (best estimates of 6.4 to 24.5 km in our study versus 2.8 to 4.7 km in previous studies) and increased percentage of trait variance explained (52.7% and 61.5% for transects 1 and 2 respectively, versus 17.6%). Revealing such strong and consistent patterns within a complex hybrid zone will allow more focused examination of the causes of variation in bimodality in mixed populations, in particular the roles of local selection versus habitat heterogeneity and gene flow between differentiated populations.

Relationships between spatio-temporal environmental and genetic variation reveal an important influence of exogenous selection in a pupfish hybrid zone

The importance of exogenous selection in a natural hybrid zone between the pupfishes Cyprinodon atrorus and Cyprinodon bifasciatus was tested via spatio-temporal analyses of environmental and genetic change over winter, spring and summer for three consecutive years. A critical influence of exogenous selection on hybrid zone regulation was demonstrated by a significant relationship between environmental (salinity and temperature) and genetic (three diagnostic nuDNA loci) variation over space and time (seasons) in the Rio Churince system, Cuatro Ciénegas, Mexico. At sites environmentally more similar to parental habitats, the genetic composition of hybrids was stable and similar to the resident parental species, whereas complex admixtures of parental and hybrid genotypic classes characterized intermediate environments, as did the greatest change in allelic and genotypic frequencies across seasons. Within hybrids across the entire Rio Churince system, seasonal changes in allelic and genotypic frequencies were consistent with results from previous reciprocal transplant experiments, which showed C. bifasciatus to suffer high mortality (75%) when exposed to the habitat of C. atrorus in winter (extreme temperature lows and variability) and summer (abrupt salinity change and extreme temperature highs and variability). Although unconfirmed, the distributional limits of C. atrorus and C. atrorus-like hybrids appear to be governed by similar constraints (predation or competition). The argument favouring evolutionary significance of hybridization in animals is bolstered by the results of this study, which links the importance of exogenous selection in a contemporary hybrid zone between C. atrorus and C. bifasciatus to previous demonstration of the long-term evolutionary significance of environmental variation and introgression on the phenotypic diversification Cuatro Ciénegas Cyprinodon.

Multilocus species tree analyses resolve the radiation of the widespread Bufo bufo species group (Anura, Bufonidae)

New analytical methods are improving our ability to reconstruct robust species trees from multilocus datasets, despite difficulties in phylogenetic reconstruction associated with recent, rapid divergence, incomplete lineage sorting and/or introgression. In this study, we applied these methods to resolve the radiation of toads in the Bufo bufo (Anura, Bufonidae) species group, ranging from the Iberian Peninsula and North Africa to Siberia, based on sequences from two mitochondrial and four nuclear DNA regions (3490 base pairs). We obtained a fully-resolved topology, with the recently described Bufo eichwaldi from the Talysh Mountains in south Azerbaijan and Iran as the sister taxon to a clade including: (1) north African, Iberian, and most French populations, referred herein to Bufo spinosus based on the implied inclusion of populations from its type locality and (2) a second clade, sister to B. spinosus, including two sister subclades: one with all samples of Bufo verrucosissimus from the Caucasus and another one with samples of B. bufo from northern France to Russia, including the Apennine and Balkan peninsulas and most of Anatolia. Coalescent-based estimations of time to most recent common ancestors for each species and selected subclades allowed historical reconstruction of the diversification of the species group in the context of Mediterranean paleogeography and indicated a long evolutionary history in this region. Finally, we used our data to delimit the ranges of the four species, particularly the more widespread and historically confused B. spinosus and B. bufo, and identify potential contact zones, some of which show striking parallels with other co-distributed species.

Measures of linkage disequilibrium among neighbouring SNPs indicate asymmetries across the house mouse hybrid zone

Theory predicts that naturally occurring hybrid zones between genetically distinct taxa can move over space and time as a result of selection and/or demographic processes, with certain types of hybrid zones being more or less likely to move. Determining whether a hybrid zone is stationary or moving has important implications for understanding evolutionary processes affecting interactions in hybrid populations. However, direct observations of hybrid zone movement are difficult to make unless the zone is moving rapidly. Here, evidence for movement in the house mouse Mus musculus domesticus × Mus musculus musculus hybrid zone is provided using measures of LD and haplotype structure among neighbouring SNP markers from across the genome. Local populations of mice across two transects in Germany and the Czech Republic were sampled, and a total of 1301 mice were genotyped at 1401 markers from the nuclear genome. Empirical measures of LD provide evidence for extinction and (re)colonization in single populations and, together with simulations, suggest hybrid zone movement because of either geography-dependent asymmetrical dispersal or selection favouring one subspecies over the other.

Pliocene and Pleistocene diversification and multiple refugia in a Eurasian shrew (Crocidura suaveolens group)

We sequenced 998 base pairs (bp) of mitochondrial DNA cytochrome b and 799 bp of nuclear gene BRCA1 in the Lesser white-toothed shrew (Crocidura suaveolens group) over its geographic range from Portugal to Japan. The aims of the study were to identify the main clades within the group and respective refugia resulting from Pleistocene glaciations. Analyses revealed the Asian lesser white-toothed shrew (C. shantungensis) as the basal clade, followed by a major branch of C. suaveolens, subdivided sensu stricto into six clades, which split-up in the Upper Pliocene and Lower Pleistocene (1.9-0.9 Myr). The largest clade, occurring over a huge range from east Europe to Mongolia, shows evidence of population expansion after a bottleneck. West European clades originated from Iberian and Italo-Balkanic refugia. In the Near East, three clades evolved in an apparent hotspot of refugia (west Turkey, south-west and south-east of the Caucasus). Most clades include specimens of different morphotypes and the validity of many taxa in the C. suaveolens group has to be re-evaluated.

BARRIERS TO SYMPATRY BETWEEN AVIAN SIBLING SPECIES (PARIDAE: BAEOLOPHUS) IN LOCAL SECONDARY CONTACT

Range limits and secondary contact zones often occur at ecotones between major associations of habitat and climate. Therefore, understanding processes that limit sympatry between species in such areas provides an important framework for testing biogeographic and evolutionary hypotheses. Theoretical and empirical work has shown that the evolution of species borders is influenced by a complexity of factors, including gene flow from central to peripheral populations and the ability of species to adapt locally to environmental conditions. However, few studies have used bioclimatic models, combined with molecular and morphological data, to predict geographic range limits in the context of gene flow across a secondary contact zone. In this study, I applied these methods to test specific hypotheses about barriers to sympatry between closely related species where they approach and contact each other. Specifically, I examined the importance of historical isolation, local adaptation, and symmetry of gene flow in limiting sympatry and range expansion of ecologically distinct species across environmental gradients. Molecular (mitochondrial DNA, allozymes), morphological, and bioclimatic data were obtained for two avian sibling species (Baeolophus inornatus and B. ridgwayi) that exist in recent, narrow secondary contact in northern California. These species are broadly allopatric and occupy rangewide associations of oak and pinyon-juniper woodlands, respectively, although B. inornatus also inhabits mixed or juniper woodlands locally. Patterns of molecular variation generally were congruent with morphological and bioclimatic data, and support prior evidence for a history of isolation, adaptation, and divergence in distinctive, species-specific vegetation-climate associations. However, molecular and morphological clines fall east of the limit of oaks, and individuals of B. inornatus in this juniper-associated contact zone experience bioclimates that are more similar to B. ridgwayi than to B. inornatus in oak habitat. Thus, B. inornatus is able to adapt and expand locally into the range of its close relative, but not vice versa. These data support the hypothesis that gene flow is asymmetrical where peripheral populations meet at range boundaries. Physiological differences between species may play an important role in influencing these patterns. Empirical studies that highlight the importance of local adaptation and patterns of gene flow in which closely related species contact across ecotones are central to understanding limits on geographic ranges, sympatry, and introgression-a cornerstone of biogeographic and speciation theory.

Phylogeography and introgressive hybridization: chipmunks (genus Tamias) in the northern Rocky Mountains

If phylogeographic studies are to be broadly used for assessing population-level processes relevant to speciation and systematics, the ability to identify and incorporate instances of hybridization into the analytical framework is essential. Here, we examine the evolutionary history of two chipmunk species, Tamias ruficaudus and Tamias amoenus, in the northern Rocky Mountains by integrating multivariate morphometrics of bacular (os penis) variation, phylogenetic estimation, and nested clade analysis with regional biogeography. Our results indicate multiple examples of mitochondrial DNA introgression layered within the evolutionary history of these nonsister species. Three of these events are most consistent with recent and/or ongoing asymmetric introgression of mitochondrial DNA across morphologically defined secondary contact zones. In addition, we find preliminary evidence where a fourth instance of nonconcordant characters may represent complete fixation of introgressed mitochondrial DNA via a more ancient hybridization event, although alternative explanations of convergence or incomplete sorting of ancestral polymorphisms cannot be dismissed with these data. The demonstration of hybridization among chipmunks with strongly differentiated bacular morphology contradicts long-standing assumptions that variation within this character is diagnostic of complete reproductive isolation within Tamias. Our results illustrate the utility of phylogeographic analyses for detecting instances of reticulate evolution and for incorporating this and other information in the inference of the evolutionary history of species.

Bimodal hybrid zones and speciation

Contact zones exemplify a series of stages in speciation. In unimodal hybrid zones intermediates predominate; in bimodal zones hybrids are rare and parental forms predominate; and finally, species might overlap, but never hybridize. Recent studies show bimodality to be associated strongly with assortative mating or fertilization, and only weakly with overall levels of genetic divergence or intrinsic genomic incompatibility. Ecological divergence across most bimodal hybrid zones suggests that ecology contributes more to speciation than genomic incompatibility. This continuum of stable contact zones provides empirical evidence for a route to speciation, which does not require allopatry.

Mother species-father species: unidirectional hybridization in animals with female choice

As mitochondria are inherited in a matrilinear way, an animal hybrid contains the mitochondrial DNA of its 'mother species'. Of 80 studies that analysed the mitochondrial DNA of at least five hybrid individuals, 50 showed that all hybrids contained the mitochondrial DNA of only one of the two parental species, indicating either mating of females of species A with males of species B but not vice versa (unidirectional hybridization) or the disappearance of one of the two parental mtDNA types after reciprocal hybridization. I review and discuss factors promoting unidirectional or reciprocal hybridization and present a sexual selection hypothesis for unidirectional hybridization. The inequality of the sexes in parental investment leads to the sex that invests more being the more discriminating one. In the presence of conspecific males, females reject allospecific males and, consequently, a male in an environment of both allospecific sexes is unlikely to mate, while in the absence of conspecific males, females sometimes accept fertilizations by males of other species. Thus, hybrid matings are usually between the females of a rare species and the males of a common species, but not vice versa. Copyright 1999 The Association for the Study of Animal Behaviour.

The evolutionary ecology of Daphnia

In order to generate genetic markers from both nuclear and mitochondrial DNA, we used three PCR-based techniques (RAPD, mtDNA-RFLPs, and sequencing of an amplified mtDNA fragment) to illuminate various aspects of the population genetics of large-lake Daphnia species. Estimations of genetic diversity at different taxonomic levels integrated with ecological data revealed insights into the genetic components of the evolutionary process of interspecific hybridization in these Daphnia species, which had previously been documented with allozyme markers. Our new molecular data suggest the occurrence of recent hybridization and backcrossing events, and allow the identification of the maternal species of hybrid clones.