Short‐term hybridisation activates Tnt1 and Tto1 Copia retrotransposons in wild tuber‐bearing Solanum species

Interspecific hybridisation in tuber‐bearing species of Solanum is a common phenomenon and represents an important source of variability, crucial for adaptation and speciation of potato species. In this regard, the effects of interspecific hybridisation on retrotransposon families present in the genomes, and their consequent effects on generation of genetic variability in wild tuber‐bearing Solanum species, are poorly characterised. The aim of this study was to analyse the activity of retrotransposons in inter‐ and intraspecific hybrids between S. kurtzianum and S. microdontum, obtained by controlled crosses, and the effects on morphological, genetic and epigenetic variability. For genetic and epigenetic analysis, S‐SAP (sequence‐specific amplification polymorphism) and TMD (transposon methylation display) techniques were used, respectively, with specific primers for Tnt1 and Tto1 retrotransposon families (Order LTR, Superfamily Copia). The results indicate that at morphological level, interspecific hybrid genotypes differ from their parental species, whereas derived intraspecific hybrids do not. In both cases, we observed significant reductions in pollen grain viability, and a negative correlation with Tnt1 mobility. Both retrotransposons, Tto1 and Tnt1, were mobilised in the genotypes analysed, with mobility ranging from 0 to 7.8%. Furthermore, at the epigenetic level, demethylation was detected in the vicinity of Tnt1 and Tto1 in the hybrids compared with the parental genotypes. These patterns were positively correlated with the activity of the retrotransposons. The results suggest a possible mechanism through which hybridisation events generate genetic variability in tuber‐bearing species of Solanum through retrotranposon activation.     

Genetic and phenotypic variation across a hybrid zone between ecologically divergent tree squirrels (Tamiasciurus)

A hybrid zone along an environmental gradient should contain a clinal pattern of genetic and phenotypic variation. This occurs because divergent selection in the two parental habitats is typically strong enough to overcome the homogenizing effects of gene flow across the environmental transition. We studied hybridization between two parapatric tree squirrels (Tamiasciurus spp.) across a forest gradient over which the two species vary in coloration, cranial morphology and body size. We sampled 397 individuals at 29 locations across a 600‐km transect to seek genetic evidence for hybridization; upon confirming hybridization, we examined levels of genetic admixture in relation to maintenance of phenotypic divergence despite potentially homogenizing gene flow. Applying population assignment analyses to microsatellite data, we found that Tamiasciurus douglasii and T. hudsonicus form two distinct genetic clusters but also hybridize, mostly within transitional forest habitat. Overall, based on this nuclear analysis, 48% of the specimens were characterized as T. douglasii, 9% as hybrids and 43% as T. hudsonicus. Hybrids appeared to be reproductively viable, as evidenced by the presence of later‐generation hybrid genotypes. Observed clines in ecologically important phenotypic traits—fur coloration and cranial morphology—were sharper than the cline of putatively neutral mtDNA, which suggests that divergent selection may maintain phenotypic distinctiveness. The relatively recent divergence of these two species (probably late Pleistocene), apparent lack of prezygotic isolating mechanisms and geographic coincidence of cline centres for both genetic and phenotypic variation suggest that environmental factors play a large role in maintaining the distinctiveness of these two species across the hybrid zone.     

Parallel plumage colour evolution and introgressive hybridization in wheatears

Genetic and phenotypic mosaics, in which various phenotypes and different genomic regions show discordant patterns of species or population divergence, offer unique opportunities to study the role of ancestral and introgressed genetic variation in phenotypic evolution. Here, we investigated the evolution of discordant phenotypic and genetic divergence in a monophyletic clade of four songbird taxa—pied wheatear (O. pleschanka), Cyprus wheatear (Oenanthe cypriaca), and western and eastern subspecies of black‐eared wheatear (O. h. hispanica and O. h. melanoleuca). Phenotypically, black back and neck sides distinguish pied and Cyprus wheatears from the white‐backed/necked black‐eared wheatears. Meanwhile, mitochondrial variation only distinguishes western black‐eared wheatear. In the absence of nuclear genetic data, and given frequent hybridization among eastern black‐eared and pied wheatear, it remains unclear whether introgression is responsible for discordance between mitochondrial divergence patterns and phenotypic similarities, or whether plumage coloration evolved in parallel. Multispecies coalescent analyses of about 20,000 SNPs obtained from RAD data mapped to a draft genome assembly resolve the species tree, provide evidence for the parallel evolution of colour phenotypes and establish western and eastern black‐eared wheatears as independent taxa that should be recognized as full species. The presence of the entire admixture spectrum in the Iranian hybrid zone and the detection of footprints of introgression from pied into eastern black‐eared wheatear beyond the hybrid zone despite strong geographic structure of ancestry proportions furthermore suggest a potential role for introgression in parallel plumage colour evolution. Our results support the importance of standing heterospecific and/or ancestral variation in phenotypic evolution.     

Experimental alteration of DNA methylation affects the phenotypic plasticity of ecologically relevant traits in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Heritable phenotypic variation in plants can be caused not only by underlying genetic differences, but also by variation in epigenetic modifications such as DNA methylation. However, we still know very little about how relevant such epigenetic variation is to the ecology and evolution of natural populations. We conducted a greenhouse experiment in which we treated a set of natural genotypes of Arabidopsis thaliana with the demethylating agent 5-azacytidine and examined the consequences of this treatment for plant traits and their phenotypic plasticity. Experimental demethylation strongly reduced the growth and fitness of plants and delayed their flowering, but the degree of this response varied significantly among genotypes. Differences in genotypes’ responses to demethylation were only weakly related to their genetic relatedness, which is consistent with the idea that natural epigenetic variation is independent of genetic variation. Demethylation also altered patterns of phenotypic plasticity, as well as the amount of phenotypic variation observed among plant individuals and genotype means. We have demonstrated that epigenetic variation can have a dramatic impact on ecologically important plant traits and their variability, as well as on the fitness of plants and their ecological interactions. Epigenetic variation may thus be an overlooked factor in the evolutionary ecology of plant populations.     

Genetic admixture despite ecological segregation in a North African sparrow hybrid zone (Aves,Passeriformes, Passer domesticus × Passer hispaniolensis)

Under different environmental conditions, hybridization between the same species might result in different patterns of genetic admixture. Particularly, species pairs with large distribution ranges and long evolutionary history may have experienced several independent hybridization events over time in different zones of overlap. In birds, the diverse hybrid populations of the house sparrow (Passer domesticus) and the Spanish sparrow (Passer hispaniolensis) provide a striking example. Throughout their range of sympatry, these two species do not regularly interbreed; however, a stabilized hybrid form (Passer italiae) exists on the Italian Peninsula and on several Mediterranean islands. The spatial distribution pattern on the Eurasian continent strongly contrasts the situation in North Africa, where house sparrows and Spanish sparrows occur in close vicinity of phenotypically intermediate populations across a broad mosaic hybrid zone. In this study, we investigate patterns of divergence and admixture among the two parental species, stabilized and nonstabilized hybrid populations in Italy and Algeria based on a mitochondrial marker, a sex chromosomal marker, and 12 microsatellite loci. In Algeria, despite strong spatial and temporal separation of urban early‐breeding house sparrows and hybrids and rural late‐breeding Spanish sparrows, we found strong genetic admixture of mitochondrial and nuclear markers across all study populations and phenotypes. That pattern of admixture in the North African hybrid zone is strikingly different from i) the Iberian area of sympatry where we observed only weak asymmetrical introgression of Spanish sparrow nuclear alleles into local house sparrow populations and ii) the very homogenous Italian sparrow population where the mitogenome of one parent (P. domesticus) and the Z‐chromosomal marker of the other parent (P. hispaniolensis) are fixed. The North African sparrow hybrids provide a further example of enhanced hybridization along with recent urbanization and anthropogenic land‐use changes in a mosaic landscape.     

Analysis of hybridization between tetraploid <Emphasis Type="Italic">Lotus corniculatus</Emphasis> and diploid <Emphasis Type="Italic">Lotus stepposus</Emphasis> (Fabaceae-Loteae): morphological and molecular aspects

We analyzed variability of morphological characters and genetic polymorphism of inter-simple sequence repeat (ISSR) markers in nine natural populations of three Lotus species from Eastern Europe, aiming to provide insights into the nature of the species L. ucrainicus. Nuclear ribosomal internal transcribed spacer (nrITS) was used as an additional molecular marker for a small subset of accessions. Analysis of variance, and principal coordinate and principal component analyses were applied for morphological data study. Cluster analysis [unweighted pair-group method with arithmetic mean (UPGMA)], principal coordinate analysis, and analysis of population genetic structure were used for ISSR pattern study. Morphological and genetic (ISSR, nrITS) evidence suggested hybrid origin of L. × ucrainicus as a result of hybridization between tetraploid species L. corniculatus and diploid species L. stepposus. We conclude that L. × ucrainicus may represent a case of hybrid speciation in statu nascendi, occurring before our very eyes.     

Jack of all nectars, master of most: DNA methylation and the epigenetic basis of niche width in a flower-living yeast

In addition to genetic differences between individuals as a result of nucleotide sequence variation, epigenetic changes that occur as a result of DNA methylation may also contribute to population niche width by enhancing phenotypic plasticity, although this intriguing possibility remains essentially untested. Using the nectar‐living yeast Metschnikowia reukaufii as study subject, we examine the hypothesis that changes in genome‐wide DNA methylation patterns underlie the ability of this fugitive species to exploit a broad resource range in its heterogeneous and patchy environment. Data on floral nectar characteristics and their use by M. reukaufii in the wild were combined with laboratory experiments and methylation‐sensitive amplified polymorphism (MSAP) analyses designed to detect epigenetic responses of single genotypes to variations in sugar environment that mimicked those occurring naturally in nectar. M. reukaufii exploited a broad range of resources, occurring in nectar of 48% of species and 52% of families surveyed, and its host plants exhibited broad intra‐ and interspecific variation in sugar‐related nectar features. Under experimental conditions, sugar composition, sugar concentration and their interaction significantly influenced the mean probability of MSAP markers experiencing a transition from unmethylated to methylated state. Alterations in methylation status were not random but predictably associated with certain markers. The methylation inhibitor 5‐azacytidine (5‐AzaC) had strong inhibitory effects on M. reukaufii proliferation in sugar‐containing media, and a direct relationship existed across sugar × concentration experimental levels linking inhibitory effect of 5‐AzaC and mean per‐marker probability of genome‐wide methylation. Environmentally induced DNA methylation polymorphisms allowed genotypes to grow successfully in extreme sugar environments, and the broad population niche width of M. reukaufii was largely made possible by epigenetic changes enabling genotype plasticity in resource use.     

Adaptive admixture in the West African bovine hybrid zone: insight from the Borgou population

Understanding the adaptive response to environmental fluctuations represents a central issue in evolutionary biology. Population admixture between divergent ancestries has often been considered as an efficient short‐term adaptation strategy. Cattle populations from the West African Bos taurus × Bos indicus hybrid zone represent a valuable resource to characterize the effect of such adaptive admixture at the genome level. We here provide a detailed assessment of the global and local genome ancestries of the Borgou breed, one of the most representative cattle of this hybrid zone. We analysed a large data set consisting of 38 100 SNPs genotyped on 203 Borgou and 591 individuals representative of all the different cattle ancestries. At the global genomic level, we show that Borgou is a stabilized admixed breed whose origin (c. 130 years ago) traces back to the great African rinderpest pandemic, several centuries after the last admixture events, the West African zebus originate from (c. 500 years ago). To identify footprints of adaptive admixture, we combined the identification of signatures of selection and the functional annotation of the underlying genes using systems biology tools. The detection of the SILV coat coloration gene likely under artificial selection may be viewed as a validation of our approach. Overall, our results suggest that the long‐term presence of pathogens and the intermediate environmental conditions are the main acting selective pressures. Our analytical framework can be extended to other model or nonmodel species to understand the process that shapes the patterns of genetic variability in hybrid zones.     

Does genetic distance between parental species influence outcomes of hybridization among coral reef butterflyfishes?

Christmas Island is located at the overlap of the Indian and Pacific Ocean marine provinces and is a hot spot for marine hybridization. Here, we evaluate the ecological framework and genetic consequences of hybridization between butterflyfishes Chaetodon guttatissimus and Chaetodon punctatofasciatus. Further, we compare our current findings to those from a previous study of hybridization between Chaetodon trifasciatus and Chaetodon lunulatus. For both species groups, habitat and dietary overlap between parental species facilitate frequent heterospecific encounters. Low abundance of potential mates promotes heterospecific pair formation and the breakdown of assortative mating. Despite similarities in ecological frameworks, the population genetic signatures of hybridization differ between the species groups. Mitochondrial and nuclear data from C. guttatissimus × C. punctatofasciatus (1% divergence at cyt b) show bidirectional maternal contributions and relatively high levels of introgression, both inside and outside the Christmas Island hybrid zone. In contrast, C. trifasciatus × C. lunulatus (5% cyt b divergence) exhibit unidirectional mitochondrial inheritance and almost no introgression. Back‐crossing of hybrid C. guttatissimus × C. punctatofasciatus and parental genotypes may eventually confound species‐specific signals within the hybrid zone. In contrast, hybrids of C. trifasciatus and C. lunulatus may coexist with and remain genetically distinct from the parents. Our results, and comparisons with hybridization studies in other reef fish families, indicate that genetic distance between hybridizing species may be a factor influencing outcomes of hybridization in reef fish, which is consistent with predictions from terrestrially derived hybridization theory.     

Digest: Transgenerational stress memory mechanisms in Arabidopsis thaliana*

Arabidopsis thaliana accessions have shown genetic diversity and type of stressor to be important determinants of transgenerational stress memory. Alvarez et al. found that certain accessions showed reversible phenotypic plasticity, supporting a model of transgenerational stress memory based upon epigenetic changes. The main proposed epigenetic regulators include DNA methylation, histone modifications, and RNA silencing via small noncoding RNA.     

Secondary contact between Lycaeides idas and L. melissa in the Rocky Mountains: extensive admixture and a patchy hybrid zone

Studies of hybridization have increased our understanding of the nature of species boundaries, the process of speciation, and the effects of hybridization on the evolution of populations and species. In the present study we use genetic and morphological data to determine the outcome and consequences of secondary contact and hybridization between the butterfly species Lycaeides idas and L. melissa in the Rocky Mountains. Admixture proportions estimated from structure and geographical cline analysis indicate L. idas and L. melissa have hybridized extensively in the Rocky Mountains and that reproductive isolation was insufficient to prevent introgression for much of the genome. Geographical patterns of admixture suggest that hybridization between L. idas and L. melissa has led to the formation of a hybrid zone. The hybrid zone is relatively wide, given estimates of dispersal for Lycaeides butterflies, and does not show strong evidence of cline concordance among characters. We believe the structure of the Lycaeides hybrid zone might be best explained by the patchy distribution of Lycaeides, local extinction and colonization of habitat patches, environmental variation and weak overall selection against hybrids. We found no evidence that hybridization in the Rocky Mountains has resulted in the formation of independent hybrid species, in contrast to the outcome of hybridization between L. idas and L. melissa in the Sierra Nevada. Finally, our results suggest that differences in male morphology between L. idas and L. melissa might contribute to isolation, or perhaps even that selection has favoured the spread of L. melissa male genitalia alleles.     

Transpositional reactivation of the <Emphasis Type="Italic">Dart</Emphasis> transposon family in rice lines derived from introgressive hybridization with <Emphasis Type="Italic">Zizania latifolia</Emphasis>

Background  

It is widely recognized that interspecific hybridization may induce "genome shock", and lead to genetic and epigenetic instabilities in the resultant hybrids and/or backcrossed introgressants. A prominent component involved in the genome shock is reactivation of cryptic transposable elements (TEs) in the hybrid genome, which is often associated with alteration in the elements' epigenetic modifications like cytosine DNA methylation. We have previously reported that introgressants derived from hybridization between Oryza sativa (rice) and Zizania latifolia manifested substantial methylation re-patterning and rampant mobilization of two TEs, a copia retrotransposon Tos17 and a MITE mPing. It was not known however whether other types of TEs had also been transpositionally reactivated in these introgressants, their relevance to alteration in cytosine methylation, and their impact on expression of adjacent cellular genes.     

Within-genotype epigenetic variation enables broad niche width in a flower living yeast

Niche theory is one of the central organizing concepts in ecology. Generally, this theory defines a given species niche as all of the factors that effect the persistence of the species as well as the impact of the species in a given location ( Hutchinson 1957 ; Chase 2011 ). Many studies have argued that phenotypic plasticity enhances niche width because plastic responses allow organisms to express advantageous phenotypes in a broader range of environments ( Bradshaw 1965 ; Van Valen 1965 ; Sultan 2001 ). Further, species that exploit habitats with fine‐grained variation, or that form metapopulations, are expected to develop broad niche widths through phenotypic plasticity ( Sultan & Spencer 2002 ; Baythavong 2011 ). Although a long history of laboratory, greenhouse and reciprocal transplant experiments have provided insight into how plasticity contributes to niche width ( Pigliucci 2001 ), recent advances in molecular approaches allow for a mechanistic understanding of plasticity at the molecular level ( Nicotra et al. 2010 ). In particular, variation in epigenetic effects is a potential source of the within‐genotype variation that underlies the phenotypic plasticity associated with broad niche widths. Epigenetic mechanisms can alter gene expression and function without altering DNA sequence ( Richards 2006 ) and may be stably transmitted across generations ( Jablonka & Raz 2009 ; Verhoeven et al. 2010 ). Also, epigenetic mechanisms may be an important component of an individual’s response to the environment ( Verhoeven et al. 2010 ). While these ideas are intriguing, few studies have made a clear connection between genome‐wide DNA methylation patterns and phenotypic plasticity (e.g. Bossdorf et al. 2010 ). In this issue of Molecular Ecology, Herrera et al. (2012) present a study that demonstrates epigenetic changes in genome‐wide DNA methylation are causally active in a species’ ability to exploit resources from a broad range of environments and are particularly important in harsh environments.     

Changes in DNA-methylation during zygotic embryogenesis in interspecific hybrids of beans (<Emphasis Type="Italic">Phaseolus</Emphasis> ssp.)

Hybrid embryos resulting from crosses between Phaseolus species often fail to reach maturity and some combinations frequently abort at early developmental stages. The genetic or molecular basis for these consistent developmental defects is at present not clear. However, an extremely complex genetic system, thought to be caused by major epigenetic changes associated with gene expression changes, has been shown to be active in plant species. We have investigated DNA methylation in two interspecific hybrids, Phaseolus vulgaris × Phaseolus coccineus and its reciprocal crosses, using methylation sensitive amplification polymorphism (MSAP). The potential use of MSAP for detecting methylation variation during embryogenesis in interspecific hybrids is discussed. Significant differences in the DNA methylation patterns were observed in abortive (interspecific hybrids) and non abortive (parental) genotypes. Taken together, our results strongly suggest that generalized alterations in DNA methylation profiles could play a causative role in early interspecific embryo abortion in vivo. A considerable change in the methylation pattern during embryogenesis could be involved in the disruption of the regulation or maintenance of the embryogenesis process of Phaseolus interspecific hybrids. The results also support the earlier hypothesis that DNA methylation is critical for the regulation of plant embryogenesis and gene expression.     

Manipulation of cytosine methylation does not remove latitudinal clines in two invasive goldenrod species in Central Europe

Invasive species frequently differentiate phenotypically in novel environments within a few generations, often even with limited genetic variation. For the invasive plants Solidago canadensis and S. gigantea, we tested whether such differentiation might have occurred through heritable epigenetic changes in cytosine methylation. In a 2‐year common‐garden experiment, we grew plants from seeds collected along a latitudinal gradient in their non‐native Central European range to test for trait differentiation and whether differentiation disappeared when seeds were treated with the demethylation agent zebularine. Microsatellite markers revealed no population structure along the latitudinal gradient in S. canadensis, but three genetic clusters in S. gigantea. Solidago canadensis showed latitudinal clines in flowering phenology and growth. In S. gigantea, the number of clonal offspring decreased with latitude. Although zebularine had a significant effect on early growth, probably through effects on cytosine methylation, latitudinal clines remained (or even got stronger) in plants raised from seeds treated with zebularine. Thus, our experiment provides no evidence that epigenetic mechanisms by selective cytosine methylation contribute to the observed phenotypic differentiation in invasive goldenrods in Central Europe.     

Rapid evolution in response to introduced predators II: the contribution of adaptive plasticity

Background  

Introductions of non-native species can significantly alter the selective environment for populations of native species, which can respond through phenotypic plasticity or genetic adaptation. We examined phenotypic and genetic responses of Daphnia populations to recent introductions of non-native fish to assess the relative roles of phenotypic plasticity versus genetic change in causing the observed patterns. The Daphnia community in alpine lakes throughout the Sierra Nevada of California (USA) is ideally suited for investigation of rapid adaptive evolution because there are multiple lakes with and without introduced fish predators. We conducted common-garden experiments involving presence or absence of chemical cues produced by fish and measured morphological and life-history traits in Daphnia melanica populations collected from lakes with contrasting fish stocking histories. The experiment allowed us to assess the degree of population differentiation due to fish predation and examine the contribution of adaptive plasticity in the response to predator introduction.     

The epigenetic footprint of poleward range‐expanding plants in apomictic dandelions

Epigenetic modifications, such as DNA methylation variation, can generate heritable phenotypic variation independent of the underlying genetic code. However, epigenetic variation in natural plant populations is poorly documented and little understood. Here, we test whether northward range expansion of obligate apomicts of the common dandelion (Taraxacum officinale) is associated with DNA methylation variation. We characterized and compared patterns of genetic and DNA methylation variation in greenhouse‐reared offspring of T. officinale that were collected along a latitudinal transect of northward range expansion in Europe. Genetic AFLP and epigenetic MS‐AFLP markers revealed high levels of local diversity and modest but significant heritable differentiation between sampling locations and between the southern, central and northern regions of the transect. Patterns of genetic and epigenetic variation were significantly correlated, reflecting the genetic control over epigenetic variation and/or the accumulation of lineage‐specific spontaneous epimutations, which may be selectively neutral. In addition, we identified a small component of DNA methylation differentiation along the transect that is independent of genetic variation. This epigenetic differentiation might reflect environment‐specific induction or, in case the DNA methylation variation affects relevant traits and fitness, selection of heritable DNA methylation variants. Such generated epigenetic variants might contribute to the adaptive capacity of individual asexual lineages under changing environments. Our results highlight the potential of heritable DNA methylation variation to contribute to population differentiation along ecological gradients. Further studies are needed using higher resolution methods to understand the functional significance of such natural occurring epigenetic differentiation.     

Genetic and epigenetic variations associated with adaptation to heterogeneous habitat conditions in a deciduous shrub

Environmentally induced phenotypic plasticity is thought to play an important role in the adaption of plant populations to heterogeneous habitat conditions, and yet the importance of epigenetic variation as a mechanism of adaptive plasticity in natural plant populations still merits further research. In this study, we investigated populations of Vitex negundo var. heterophylla (Chinese chastetree) from adjacent habitat types at seven sampling sites. Using several functional traits, we detected a significant differentiation between habitat types. With amplified fragment length polymorphisms (AFLP) and methylation‐sensitive AFLP (MSAP), we found relatively high levels of genetic and epigenetic diversity but very low genetic and epigenetic differences between habitats within sites. Bayesian clustering showed a remarkable habitat‐related differentiation and more genetic loci associated with the habitat type than epigenetic, suggesting that the adaptation to the habitat is genetically based. However, we did not find any significant correlation between genetic or epigenetic variation and habitat using simple and partial Mantel tests. Moreover, we found no correlation between genetic and ecologically relevant phenotypic variation and a significant correlation between epigenetic and phenotypic variation. Although we did not find any direct relationship between epigenetic variation and habitat environment, our findings suggest that epigenetic variation may complement genetic variation as a source of functional phenotypic diversity associated with adaptation to the heterogeneous habitat in natural plant populations.     

Geographic variation in hybridization and ecological differentiation between three syntopic,morphologically similar species of montane lizards

To understand factors shaping species boundaries in closely related taxa, a powerful approach is to compare levels of genetic admixture at multiple points of contact and determine how this relates to intrinsic and extrinsic factors, such as genetic, morphological and ecological differentiation. In the Australian Alps, the threatened alpine bog skink Pseudemoia cryodroma co‐occurs with two morphologically and ecologically similar congeners, P. entrecasteauxii and P. pagenstecheri, and all three species are suspected to hybridize. We predicted that the frequency of hybridization should be negatively correlated with genetic divergence, morphological differentiation and microhabitat separation. We tested this hypothesis using a mitochondrial locus, 13 microsatellite loci, morphological and microhabitat data and compared results across three geographically isolated sites. Despite strong genetic structure between species, we detected hybridization between all species pairs, including evidence of backcrossed individuals at the two sites where all three species are syntopic. Hybridization frequencies were not consistently associated with genetic, morphological or ecological differentiation. Furthermore, P. entrecasteauxii and P. pagenstecheri only hybridized at the two sites where they are syntopic with P. cryodroma, but not at the largest site where P. cryodroma was not recorded, suggesting that P. cryodroma may serve as a bridging species. This study reveals the complex dynamics within a three species hybrid zone and provides a baseline for assessing the impact of climate change and anthropogenic habitat modification on future hybridization frequencies.