Site-specific temporal variation of population dynamics in subalpine endemic plant species.

Endemic plants in high mountains are projected to be at high risk because of climate change. Temporal demographic variation is a major factor affecting population viability because plants often occur in small, isolated populations. Because isolated populations tend to exhibit genetic differentiation, analyzing temporal demographic variation in multiple populations is required for the management of high mountain endemic species. We examined the population dynamics of an endemic plant species, Primula farinosa subsp. modesta, in four subalpine sites over six years. Stage-based transition matrices were constructed, and temporal variation in the projected population growth rate (λ) was analyzed using life table response experiments (LTREs). The variation in λ was primarily explained by the site × year interaction rather than the main effects of the site and year. The testing sites exhibited inconsistent patterns in the LTRE contributions of the vital rates to the temporal deviation of λ. However, within sites, growth or stasis had significant negative correlations with temporal λ deviation. Negative correlations among the contributions of vital rates were also detected within the two testing sites, and the removal of the correlations alleviated temporal fluctuations in λ. The response of vital rates to yearly environmental fluctuations reduced the temporal variation of λ. Such effects manifested especially at two sites where plants exhibited higher plasticity than plants at other sites. Site-specific temporal variation implies that populations of high mountain species likely exhibit asynchronous temporal changes, and multiple sites need to be evaluated for their conservation.

Fiber-optic refractometer for in vivo sugar concentration measurements of low-nectar-producing flowers.

Sugar concentration in floral nectars is an assessment required in several diverse fields of application. The widely used analysis, consisting of nectar extraction with a microcapillary and sugar concentration measurement with a light refractometer, is not reliable when the nectar is secreted in small quantities, unextractable with a microcapillary. Ancillary methods adopted in such cases are destructive, rather complicated and often provide much less precise and accurate results. The microscopic-size, low cost and biocompatibility of optical fibers were exploited to deliver light directly inside the flower with minimal invasiveness and measure instantaneously the refractometric properties of the nectar without extracting it. After comparing the new and old methods using two known nectariferous species, the new approach was validated on Primula palinuri, whose nectar is unextractable with microcapillaries. The fiber-optic probe was able to measure the nectar refractive index in P. palinuri flowers making it possible to highlight a previously undetected significant trend of the sugar concentration throughout the long anthesis of the single flowers. Changes in nectar concentrations are similar in both longistylous and brevistylous flowers. The fiber-optic refractometer is an advancement of light refractometer analysis. Further customization of the laboratory set-up into portable equipment will boost applications.

Conservation implications of asymmetric introgression and reproductive barriers in a rare primrose species.

BACKGROUND: Primula is a large genus of flowering herbs well known for their heterostyly. Currently few natural hybrids are known and reproductive barriers in this genus in the wild have received little attention. However, there is instance of hybridization between rare and widely-spread species, and conservation implications of such situation is poorly understood. In the present study, we investigated hybridization patterns and reproductive barriers between a wide spread species, Primula poissonii and a rare species P. anisodora, of which only three populations are currently known. RESULTS: Pollinator-mediated reproductive isolation was strong between parental species but not significant between hybrids and parental species. Hand pollination experiments showed significant reduction of both fruit- and seed-set for heterospecific pollination as compared with conspecific pollination for both parental species. Furthermore, hybrids had higher fruit- and seed-set when pollinated with P. anisodora pollen as opposed to P. poissonii pollen. Microsatellites identified backcrosses to P. anisodora in two of the three populations of P. anisodora, and additionally more individuals of P. anisodora showed introgression from P. poissonii than vice versa. CONCLUSIONS: These results provide evidence for potential genetic swamping of the P. anisodora populations, which could pose a serious threat for this locally endemic species.

Nectar robbers influence the trait-fitness relationship of Primula secundiflora.

The trait-fitness relationship influences the strength and direction of floral evolution. To fully understand and predict the evolutionary trajectories of floral traits, it is critical to disentangle the direct and indirect effects of floral traits on plant fitness in natural populations. We experimentally quantified phenotypic selection on floral traits through female fitness and estimated the casual effects of nectar robbing with different nectar robbing intensities on trait-fitness relationships in both the L- (long-style and short-anther phenotype) and S-morph (short-style and long-anther phenotype) flowers among Primula secundiflora populations. A larger number of flowers and wider corolla tubes had both direct and indirect positive effects on female fitness in the P. secundiflora populations. The indirect effects of these two traits on female fitness were mediated by nectar robbers. The indirect effect of the number of flowers on female fitness increased with increasing nectar robbing intensity. In most populations, the direct and/or indirect effects of floral traits on female fitness were stronger in the S-morph flowers than in the L-morph flowers. In addition, nectar robbers had a direct positive effect on female fitness, but this effect varied between the L- and S-morph flowers. These results show the potential role of nectar robbers in influencing the trait-fitness relationships in this primrose species.

Responses of Primula vulgaris to light quality in the maternal and germination environments.

In the model species Arabidopsis thaliana phytochromes mediate dormancy and germination responses to seasonal cues experienced during seed maturation on the maternal plants. However, the effect of the maternal light environment on seed germination in native wild species has not been well studied. This is particularly important given its practical application in the context of environmental restoration, when there can be marked changes in the canopy. Plants of Primula vulgaris were grown in the field over two vegetative seasons under four shading treatments from low to high ratio of red to far-red light (R:FR). Leaf and seed traits were assessed in response to the light treatments. The germination of seeds from these four maternal environments (pre-dispersal) was investigated at seven light and five temperature treatments (post-dispersal). Thinner leaves, larger leaf area and greater chlorophyll content were found in plants growing in reduced R:FR. Shading in the maternal environment led to increased seed size and yield, although the conditions experienced by the maternal plants had no effect on seed germination. Seeds responded strongly to the cues experienced in their immediate germination environment. Germination was always enhanced under higher R:FR conditions. The observed phenotypic trait variation plays a major role in the ability of P. vulgaris to grow in a wide range of light conditions. However, the increased germination capacity in response to a higher R:FR for all maternal environments suggests potential for seedling establishment under vegetative shade only in the presence of canopy gaps.

Effects of climate warming and prolonged snow cover on phenology of the early life history stages of four alpine herbs on the southeastern Tibetan Plateau.

PREMISE OF THE STUDY: Much research has focused on plant responses to ongoing climate change, but there is relatively little information about how climate change will affect the early plant life history stages. Understanding how global warming and changes in winter snow pattern will affect seed germination and seedling establishment is crucial for predicting future alpine population and vegetation dynamics. METHODS: In a 2-year study, we tested how warming and alteration in the snowmelt regime, both in isolation and combination, influence seedling emergence phenology, first-year growth, biomass allocation, and survival of four native alpine perennial herbs on the southeastern Tibetan Plateau. KEY RESULTS: Warming promoted seedling emergence phenology of all four species and biomass per plant of two species but reduced seedling survival of three species. Prolonged snow cover partly mediated the affects of warming on Primula alpicola (survival and biomass), Pedicularis fletcheri (phenology, biomass, and root:shoot ratio) and Meconopsis integrifolia (survival). For the narrowly distributed species M. racemosa, seedling growth was additively decreased by warming and prolonged snow cover. CONCLUSIONS: Both warming and alteration of the snow cover regime can influence plant recruitment by affecting seedling phenology, growth, and survival, and the effects are largely species-specific. Thus, climate change is likely to affect population dynamics and community structure of the alpine ecosystem. This is the first experimental demonstration of the phenological advancement of seedling emergence in the field by simulated climate warming.

Comparative intra- and interspecific sexual organ reciprocity in four distylous Primula species in the Himalaya-Hengduan Mountains.

Distyly is a mechanism promoting cross-pollination within a balanced polymorphism. Numerous studies show that the degree of inter-morph sexual organ reciprocity (SOR) within species relates to its pollen-mediated gene flow. Similarly, a lower interspecific SOR should promote interspecific isolation when congeners are sympatric, co-blooming and share pollinators. In this comparative study, we address the significance of SOR at both intra- and interspecific levels. Seventeen allopatric and eight sympatric populations representing four Primula species (P. anisodora, P. beesiana, P. bulleyana and P. poissonii) native to the Himalaya-Hengduan Mountains were measured for eight floral traits in both long- and short-styled morphs. GLMM and spatial overlap methods were used to compare intra- and interspecific SOR. While floral morphology differed among four Primula species, SOR within species was generally higher than between species, but in species pairs P. poissonii/P. anisodora and P. beesiana/P. bulleyana, the SOR was high at both intra- and interspecific levels. We did not detect a significant variation in intraspecific SOR or interspecific SOR when comparing allopatric versus sympatric populations for all species studied. As intraspecific SOR increased, disassortative mating may be promoted. As interspecific SOR decreased, interspecific isolation between co-flowering species pairs also may increase. Hybridisation between congeners occurred when interspecific SOR increased in sympatric populations, as confirmed in two species pairs, P. poissonii/P. anisodora and P. beesiana/P. bulleyana.

Pollinator-mediated selection on flowering phenology and floral display in a distylous herb Primula alpicola.

The targets and causes of phenotypic selection are crucial to understanding evolutionary ecology. However, few studies have examined selection quantitatively from multiple sources on the same trait identified the agent of natural selection experimentally. Here we quantified phenotypic selection on traits, including flowering phenology and aspects of floral display via female fitness, in the distylous perennial herb Primula alpicola. To determine the role of pollinators in generating selection effects on floral traits, we compared the phenotypic selection gradients in open-pollinated and hand-pollinated plants. Our results show that pollinator-mediated linear selection on flowering start and correlational selection on the number of flowers and scape height explains most of the net phenotypic selection on these traits suggesting pollinators played an important role in shaping floral diversity. We used path analysis and structural equation modeling (SEM) to examine how herbivores affected the relationship between floral traits and female fitness, but no significant selection was caused by seed predators. These results suggest pollinators, not herbivores maybe the significant agent of selection on flora traits.

Ecological correlates and genetic consequences of evolutionary transitions from distyly to homostyly.

Background and Aims: The outbreeding floral polymorphism heterostyly frequently breaks down, resulting in the evolution of self-fertilization as a result of homostyle formation. Here, the loss of floral polymorphism in distylous Primula oreodoxa, a sub-alpine species restricted to western Sichuan, China, was examined by investigating the ecological correlates and genetic consequences of mating system transitions. Several key questions were addressed. (1) What are the frequencies, geographical distribution and reproductive characteristics of floral morphs in distylous and homostylous populations? (2) Does increased elevation influence pollinator service and the likelihood of inbreeding in populations? (3) How often has homostyly originated and what are the consequences of the breakdown of distyly for the amounts and distribution of genetic diversity in populations? Methods: Fourteen populations throughout the range of P. oreodoxa were sampled, and morph frequencies and floral characteristics were recorded. Polymorphism at microsatellite loci and chloroplast DNA (cpDNA) variation were used to quantify population genetic structure and genetic relationships among populations. Controlled pollinations and studies of pollen tube growth and fertility were conducted to determine the compatibility status of populations and their facility for autonomous self-pollination. Finally, visitation rates of long- and short-tongued pollinators to distylous and homostylous populations at different elevations were compared to determine if increased elevation was associated with deterioration in pollinator service. Key Results: In contrast to most heterostylous species, both distylous and homostylous morphs of P. oreodoxa are highly self-compatible, but only homostyles have the facility for autonomous self-pollination. Homostyles set significantly more fruit and seeds following open pollination than the distylous morphs. Visitation by long-tongued pollinators was significantly lower in homostylous populations, and overall rates of insect visitation decreased with elevation. Genetic diversity was significantly lower in homostylous populations, with evidence of increased inbreeding at higher elevation. Patterns of cpDNA variation were consistent with multiple transitions from distyly to homostyly and limited gene flow among populations. Conclusions: The results of this study support the hypothesis that the multiple loss of floral polymorphism in distylous P. oreodoxa is associated with unsatisfactory pollinator service, with homostyles benefiting from reproductive assurance as a result of autonomous self-pollination.

The measure and mismeasure of reciprocity in heterostylous flowers.

The goal of biological measurement is to capture underlying biological phenomena in numerical form. The reciprocity index applied to heterostylous flowers is meant to measure the degree of correspondence between fertile parts of opposite sex on complementary (inter-compatible) morphs, reflecting the correspondence of locations of pollen placement on, and stigma contact with, pollinators. Pollen of typical heterostylous flowers can achieve unimpeded fertilization only on opposite-morph flowers. Thus, the implicit goal of this measurement is to assess the likelihood of 'legitimate' pollinations between compatible morphs, and hence reproductive fitness. Previous reciprocity metrics fall short of this goal on both empirical and theoretical grounds. We propose a new measure of reciprocity based on theory that relates floral morphology to reproductive fitness. This method establishes a scale based on adaptive inaccuracy, a measure of the fitness cost of the deviation of phenotypes in a population from the optimal phenotype. Inaccuracy allows the estimation of independent contributions of maladaptive bias (mean departure from optimum) and imprecision (within-population variance) to the phenotypic mismatch (inaccuracy) of heterostylous morphs on a common scale. We illustrate this measure using data from three species of Primula (Primulaceae).

How Do Cold-Adapted Plants Respond to Climatic Cycles? Interglacial Expansion Explains Current Distribution and Genomic Diversity in Primula farinosa L.

Understanding the effects of past climatic fluctuations on the distribution and population-size dynamics of cold-adapted species is essential for predicting their responses to ongoing global climate change. In spite of the heterogeneity of cold-adapted species, two main contrasting hypotheses have been proposed to explain their responses to Late Quaternary glacial cycles, namely, the interglacial contraction versus the interglacial expansion hypotheses. Here, we use the cold-adapted plant Primula farinosa to test two demographic models under each of the two alternative hypotheses and a fifth, null model. We first approximate the time and extent of demographic contractions and expansions during the Late Quaternary by projecting species distribution models across the last 72 ka. We also generate genome-wide sequence data using a Reduced Representation Library approach to reconstruct the spatial structure, genetic diversity, and phylogenetic relationships of lineages within P. farinosa. Finally, by integrating the results of climatic and genomic analyses in an Approximate Bayesian Computation framework, we propose the most likely model for the extent and direction of population-size changes in $P$. farinosa through the Late Quaternary. Our results support the interglacial expansion of $P$. farinosa, differing from the prevailing paradigm that the observed distribution of cold-adapted species currently fragmented in high altitude and latitude regions reflects the consequences of postglacial contraction processes.

Do floral and niche shifts favour the establishment and persistence of newly arisen polyploids? A case study in an Alpine primrose.

BACKGROUND AND AIMS: Polyploidization plays a key role in plant evolution. Despite the generally accepted 'minority-cytotype exclusion' theory, the specific mechanisms leading to successful establishment and persistence of new polyploids remain controversial. The majority of newly formed polyploids do not become established, because they are less common, have fewer potential mates or may not be able to compete successfully with co-occurring progenitors at lower ploidy levels. Changes in floral traits and ecological niches have been proposed as important mechanisms to overcome this initial frequency-dependent disadvantage. The aim of this study was to determine whether dodecaploids of the heterostylous P. marginata differ from their hexaploid progenitors in P. marginata and P. allionii for selected floral traits and ecological preferences that might be involved in establishment and persistence, providing a possible explanation for the origin of polyploidized populations. METHODS: Floral morphological traits and ecological niche preferences among dodecaploids and their hexaploid progenitors in P. marginata and P. allionii ,: all restricted to the south-western Alps, were quantified and compared KEY RESULTS: Differences in floral traits were detected between dodecaploids and their closest relatives, but such differences might be too weak to counter the strength of minority cytotype disadvantage and are unlikely to enable the coexistence of different cytotypes. Furthermore, the results suggest the preservation of full distyly and no transition to selfing in dodecaploids. Finally, dodecaploids occur almost exclusively in environments that are predicted to be suitable also for their closest hexaploid relatives. CONCLUSIONS: In light of the results, P. marginata dodecaploids have probably been able to establish and persist by occupying geographical areas not yet filled by their closest relatives without significant evolution in their climatic and pollination niches. Dispersal limitation and minority-cytotype exclusion probably maintain their current range disjunct from those of its close relatives.

Addressing potential local adaptation in species distribution models: implications for conservation under climate change.

Species distribution models (SDMs) have been criticized for involving assumptions that ignore or categorize many ecologically relevant factors such as dispersal ability and biotic interactions. Another potential source of model error is the assumption that species are ecologically uniform in their climatic tolerances across their range. Typically, SDMs treat a species as a single entity, although populations of many species differ due to local adaptation or other genetic differentiation. Not taking local adaptation into account may lead to incorrect range prediction and therefore misplaced conservation efforts. A constraint is that we often do not know the degree to which populations are locally adapted. Lacking experimental evidence, we still can evaluate niche differentiation within a species' range to promote better conservation decisions. We explore possible conservation implications of making type I or type II errors in this context. For each of two species, we construct three separate Max-Ent models, one considering the species as a single population and two of disjunct populations. Principal component analyses and response curves indicate different climate characteristics in the current environments of the populations. Model projections into future climates indicate minimal overlap between areas predicted to be climatically suitable by the whole species vs. population-based models. We present a workflow for addressing uncertainty surrounding local adaptation in SDM application and illustrate the value of conducting population-based models to compare with whole-species models. These comparisons might result in more cautious management actions when alternative range outcomes are considered.

Forest succession and population viability of grassland plants: long repayment of extinction debt in Primula veris.

Time lags in responses of organisms to deteriorating environmental conditions delay population declines and extinctions. We examined how local processes at the population level contribute to extinction debt, and how cycles of habitat deterioration and recovery may delay extinction. We carried out a demographic analysis of the fate of the grassland perennial Primula veris after the cessation of grassland management, where we used either a unidirectional succession model for forest habitat or a rotation model with a period of forest growth followed by a clear-cut and a new successional cycle. The simulations indicated that P. veris populations may have an extinction time of decades to centuries after a detrimental management change. A survey of the current incidence and abundance of P. veris in sites with different histories of afforestation confirmed the simulation results of low extinction rates. P. veris had reduced incidence and abundance only at sites with at least 100 years of forest cover. Time to extinction in simulations was dependent on the duration of the periods with favourable and unfavourable conditions after management cessation, and the population sizes and growth rates in these periods. Our results thus suggest that the ability of a species to survive is a complex function of disturbance regimes, rates of successional change, and the demographic response to environmental changes. Detailed demographic studies over entire successional cycles are therefore essential to identify the environmental conditions that enable long-term persistence and to design management for species experiencing extinction debts.

Asymmetrical disassortative pollination in a distylous primrose: the complementary roles of bumblebee nectar robbers and syrphid flies.

Heterostyly is a floral polymorphism characterized by reciprocal herkogamy maintained through high levels of mating between morphs, serviced by appropriate pollinators. We studied how differential efficiency and abundance of distinct pollinators affect plant female reproduction in self- and intra-morph incompatible distylous Primula secundiflora. Bumblebees and syrphid flies were found to be the most abundant floral visitors. Bumblebees frequently exhibited nectar-robbing behavior. Because the robbing holes were always situated between the high- and low-level organs on both morphs, nectar-robbing bumblebees only pollinated S-styled flowers. L-styled flowers set four times as many seeds as did S-styled flowers after being visited by pollen-collecting syrphid flies. The natural female fecundity and the magnitude of pollen limitation varied between the morphs within populations because of the mosaic distribution of nectar-robbing bumblebees and syrphid flies. L-styled flowers and S-styled flowers set the same number of seeds after supplemental hand pollination, indicating equivalent female reproductive potential. We suggest that bumblebee nectar robbers and syrphid flies play an important role in sustaining the floral dimorphism of heterostyly in P. secundiflora because of their complementary roles in the pollination system.

Reproducing under a warming climate: long winter flowering and extended flower longevity in the only Mediterranean and maritime Primula.

Under the pressure of global warming, general expectations of species migration and evolution of adaptive traits should always be confirmed with species-specific studies. Within this framework, some species can be used as study systems to predict possible consequences of global warming also on other relatives. Unlike its mountain congeneric, Primula palinuri Petagn. has endured all the climatic fluctuations since the Pleistocene, while surviving on Mediterranean coastal cliffs. The aim of this work was to investigate the possible evolution of reproductive biological and ecological traits in P. palinuri adaptation to a warmer environment. Data showed that flowering starts in mid-winter; single flowers remain open for over a month, changing from pendulous to erect. The number of insects visiting flowers of P. palinuri increases during the flowering season, and pollination reduces flower longevity. Overall, the best pollen performances, in terms of viability and germinability, occur at winter temperatures, while pollinator activity prolongs flowering until spring. Moreover, extended longevity of single flowers optimises reproductive success. Both phenotypic plasticity and selective processes might have occurred in P. palinuri. However, we found that reproductive traits of the only Mediterranean Primula remain more associated with cold mountain habitats than warm coastal cliffs. Given the rapid trend of climate warming, migration and new adaptive processes in P. palinuri are unlikely. Response to past climate warming of P. palinuri provides useful indications for future scenarios in other Primula species.

Disruption of the distylous syndrome in Primula veris.

BACKGROUND AND AIMS: Distyly is a floral polymorphism characterized by the presence of two discrete morphs with reciprocal positioning of anthers and stigmas in flowers on different plants within the same population. Although reciprocal herkogamy and associated floral traits are generally thought to be discrete and strict polymorphisms, little is known about variation in floral traits related to the distylous syndrome within and among populations of a single species. In this study, variation in floral morphology and reciprocal positioning of the sexual organs in the distylous Primula veris (cowslip) is quantified. METHODS: Data were collected in ten populations occurring in two contrasting habitat types (grasslands and forests), and for each population the average level of reciprocity was assessed, the strength of the self-incompatibility system was determined, and seed production under natural conditions was quantified. RESULTS: In grassland populations, flowers showed clear distyly with low and symmetric reciprocity indices at both the lower and upper level. In forests, P. veris produced larger flowers that showed strong deviations in stigma-anther separation, especially in the L-morph. This deviation was mainly driven by variation in stigma height, resulting in high and asymmetric reciprocity indices and the occurrence of several short-styled homostylous plants. Self-incompatibility was, however, strict in both habitats, and morph ratios did not differ significantly from isoplethy. The observed shift in reciprocity in forest populations was associated with a significant reduction in seed production in the L-morph. CONCLUSIONS: The results indicate that populations of P. veris show habitat-specific variation in flower morphology. Deviations from perfect reciprocal positioning of stigmas and anthers also translate into reduced seed production, suggesting that small changes in sexual organ reciprocity can have far-reaching ecological and evolutionary implications.

Keeping pace with climate change: stage-structured moving-habitat models.

Life cycles can limit the abilities of species to track changing climatic conditions. We combined age or stage structure and a moving-habitat model to explore the effects of life history on the persistence of populations in the presence of climate change. We studied four dissimilar plant species in moving patches and found that (1) population growth rates, (2) elasticities with respect to the survival (stasis and shrinkage) components of the projection matrix, and (3) the evenness of the elasticities with respect to the components of the projection matrix all decreased as we increased the translational speeds of the patches. In addition, the value of long-distance dispersal increased with patch speed for three of the four species. Our analyses confirm that rapid growth, high fecundity, and long-distance dispersal can benefit species in moving patches. Thus, species with long generation times and limited dispersal ability are especially vulnerable to habitat movement. Stage-structured moving-habitat models can easily incorporate spatial complexity and can help us predict the effects of shifting climatic conditions.

Environmental context influences both the intensity of seed predation and plant demographic sensitivity to attack.

Variation in mutualistic and antagonistic interactions are important sources of variation in population dynamics and natural selection. Environmental heterogeneity can influence the outcome of interactions by affecting the intensity of interactions, but also by affecting the demography of the populations involved. However, little is known about the relative importance of environmental effects on interaction intensities and demographic sensitivity for variation in population growth rates. We investigated how soil depth, soil moisture, soil nutrient composition, and vegetation height influenced the intensity of seed predation as well as host plant demography and sensitivity to seed predation in the perennial herb Primula farinosa. Intensity of seed predation ranged from 0% to 80% of seeds damaged among the 24 study populations and was related to soil moisture in two of four years. The effect of seed predation on plant population growth rate (lambda) ranged from negligible to a reduction in lambda by 0.70. Sensitivity of population growth rate to predation explained as much of the variation in the reductions in population growth rate due to seed predation as did predation intensity. Plant population growth rate in the absence of seed predation and sensitivity to predation were negatively related to soil depth and soil moisture. Both intensity of predation and sensitivity to predation were positively correlated with potential population growth rate and, as a result, there was no significant relationship between predation intensity and realized population growth rate. We conclude that in our study system environmental context influences the effects of seed predation on plant fitness and population dynamics in two important ways: through variation in interaction intensity and through sensitivity to the effects of this interaction. Moreover, our results show that a given abiotic factor can influence population growth rate in different directions through effects on potential growth rate, intensity of biotic interactions, and the sensitivity of population growth rate to interactions.

Unidirectional hybridization and reproductive barriers between two heterostylous primrose species in north-west Yunnan, China.

BACKGROUND AND AIMS: Heteromorphy in flowers has a profound effect on breeding patterns within a species, but little is known about how it affects reproductive barriers between species. The heterostylous genus Primula is very diverse in the Himalaya region, but hybrids there have been little researched. This study examines in detail a natural hybrid zone between P. beesiana and P. bulleyana. METHODS: Chloroplast sequencing, AFLP (amplified fragment length polymorphism) markers and morphological comparisons were employed to characterize putative hybrids in the field, using synthetic F1s from hand pollination as controls. Pollinator visits to parent species and hybrids were observed in the field. Hand pollinations were conducted to compare pollen tube growth, seed production and seed viability for crosses involving different morphs, species and directions of crossing. KEY RESULTS: Molecular data revealed all hybrid derivatives examined to be backcrosses of first or later generations towards P. bulleyana: all had the chloroplast DNA (cpDNA) of this species. Some individuals had morphological traits suggesting they were hybrids, but they were genetically similar to P. bulleyana; they might have been advanced generation backcrosses. Viable F1s could not be produced with P. bulleyana pollen on P. beesiana females, irrespective of the flower morphs used. Within-morph crosses for each species had very low (<10 %) seed viability, whereas crosses between pin P. bulleyana (female) and pin P. beesiana had a higher seed viability of 30 %. Thus genetic incompatibility mechanisms back up mechanical barriers to within-morph crosses in each species, but are not the same between the two species. The two species share their main pollinators, and pollinators were observed to fly between P. bulleyana and hybrids, suggesting that pollinator behaviour may not be an important isolating factor. CONCLUSIONS: Hybridization is strongly asymmetric, with P. bulleyana the only possible mother and all detected hybrids being backcrosses in this direction. Partial ecological isolation and inhibition of heterospecific pollen, and possibly complete barriers to F1 formation on P. beesiana, may be enough to make F1 formation very rare in these species. Therefore, with no F1 detected, this hybrid zone may have a finite life span as successive generations become more similar to P. bulleyana.