The genetic architecture of reproductive isolation in Louisiana irises: pollination syndromes and pollinator preferences

In animal-pollinated plants, pollinator preferences for divergent floral forms can lead to partial reproductive isolation. We describe regions of plant genomes that affect pollinator preferences for two species of Louisiana Irises, Iris brevicaulis and Iris fulva, and their artificial hybrids. Iris brevicaulis and I. fulva possess bee and bird-pollination syndromes, respectively. Hummingbirds preferred I. fulva and under-visited both I. brevicaulis and backcrosses toward this species. Lepidopterans preferred I. fulva and backcrosses toward I. fulva, but also under-visited I. brevicaulis and I. brevicaulis backcrosses. Bumblebees preferred I. brevicaulis and F1 hybrids and rarely visited I. fulva. Although all three pollen vectors preferred one or the other species, these preferences did not prevent visitation to other hybrid/parental classes. Quantitative trait locus (QTL) mapping, in reciprocal BC1 mapping populations, defined the genetic architecture of loci that affected pollinator behavior. We detected six and nine QTLs that affected pollinator visitation rates in the BCIb and BCIf mapping populations, respectively, with as many as three QTLs detected for each trait. Overall, this study reflects the possible role of quantitative genetic factors in determining (1) reproductive isolation, (2) the pattern of pollinator-mediated genetic exchange, and thus (3) hybrid zone evolution.     

Intimately linked or hardly speaking? The relationship between genotype and environmental gradients in a Louisiana Iris hybrid population

Several models of hybrid zone evolution predict the same spatial patterns of genotypic distribution whether or not structuring is due to environment-dependent or -independent selection. In this study, we tested for evidence of environment-dependent selection in an Iris fulva x Iris brevicaulis hybrid population by examining the distribution of genotypes in relation to environmental gradients. We selected 201 Louisiana Iris plants from within a known hybrid population (80 m x 80 m) and placed them in four different genotypic classes (I. fulva, I. fulva-like hybrid, I. brevicaulis-like hybrid and I. brevicaulis) based on seven species-specific random amplified polymorphic DNA (RAPD) markers and two chloroplast DNA haplotypes. Environmental variables were then measured. These variables included percentage cover by tree canopy, elevation from the high water mark, soil pH and percentage soil organic matter. Each variable was sampled for all 201 plants. Canonical discriminant analysis (CDA) was used to infer the environmental factors most strongly associated with the different genotypic groups. Slight differences in elevation (-0.5 m to +0.4 m) were important for distinguishing habitat distributions described by CDA, even though there were no statistical differences between mean elevations alone. I. brevicaulis occurred in a broad range of habitats, while I. fulva had a narrower distribution. Of all the possible combinations, I. fulva-like hybrids and I. brevicaulis-like hybrids occurred in the most distinct habitat types relative to one another. Each hybrid class was not significantly different from its closest parent with regard to habitat occupied, but was statistically unique from its more distant parental species. Within the hybrid genotypes, most, but not all, RAPD loci were individually correlated with environmental variables. This study suggests that, at a very fine spatial scale, environment-dependent selection contributed to the genetic structuring of this hybrid zone.     

Nectar Production in Louisiana Iris Hybrids

Nectar is an important attractant for pollinators, and a plant's success in sexual reproduction can be influenced by the amount and concentration of nectar produced by its flowers. We studied nectar production over flower lifetime in Iris fulva, Iris brevicaulis, and four classes of hybrids-reciprocal F1's and backcrosses-between these species. Iris fulva produced less concentrated nectar than did I. brevicaulis, whereas I. brevicaulis flowers had a shorter life span. Hybrids were not intermediate, but they had the high nectar concentration of I. brevicaulis combined with the long life span of I. fulva flowers. Nectar production and concentration declined after the first day in all classes, but flowers continued to produce nectar until they were completely wilted. Backcrosses did not show a shift in mean or increased variation for the characters that distinguished the parental species; backcrosses toward I. fulva retained the high nectar concentration of I. brevicaulis, and backcrosses toward I. brevicaulis did not have a reduced flower life span. Overall, F1 hybrid flowers produced the highest amounts of nectar and nectar sugar over their life spans. These results, together with previously obtained data on pollinator choice in mixed arrays of the same flower classes, show that F1 hybrids between these species do not suffer from reduced attractiveness to pollinators. F1 individuals produced more nectar and nectar sugar than did their parents, and thus, they are possibly even more attractive to pollinators that forage for nectar.     

Loci affecting long-term hybrid survivorship in Louisiana irises: implications for reproductive isolation and introgression

Iris fulva and I. brevicaulis are long-lived plant species known to hybridize where they coexist in nature. Year-to-year survival contributes significantly to overall fitness for both species and their hybrid derivatives, and differences in hybrid survivability may have important consequences to interspecific gene flow in nature. We examined the genetic architecture of long-term survivorship of reciprocal backcross I. fulva x I. brevicaulis hybrids in a common-garden, greenhouse environment. Differences in mortality were found between the two backcross (BC1) hybrid classes, with hybrids crossed toward I. fulva (BCIF) revealing twice the mortality of those hybrids backcrossed toward I. brevicaulis (BCIB). Using genomic scans on two separate genetic linkage maps derived from the reciprocal hybrid populations, we found that hybrid survivorship is influenced by several genetic regions. Multiple interval mapping (MIM) revealed four quantitative trait loci (QTLs) in BCIF hybrids that were significantly associated with survivorship. Introgressed I. brevicaulis DNA increased survivorship at three of the four QTLs. For the fourth QTL, introgressed I. brevicaulis DNA was associated with decreased survivorship. No QTLs were detected in BCIB hybrids; however, single-marker analysis revealed five unlinked loci that were significantly associated with survivorship. At all five markers, survivorship was positively associated with introgressed I. fulva DNA. The present findings have important implications for the evolutionary dynamics of naturally occurring hybrid zones. Regions of the genome that increase survivorship when in a heterozygous (i.e., hybrid) state should have an increased likelihood of passing across species boundaries, whereas those that decrease survivorship will be less likely to introgress.     

Detecting adaptive trait introgression between Iris fulva and I. brevicaulis in highly selective field conditions

The idea that natural hybridization has served as an important force in evolutionary and adaptive diversification has gained considerable momentum in recent years. By combining genome analyses with a highly selective field experiment, we provide evidence for adaptive trait introgression between two naturally hybridizing Louisiana Iris species, flood-tolerant Iris fulva and dry-adapted I. brevicaulis. We planted reciprocal backcross (BC1) hybrids along with pure-species plants into natural settings that, due to a flooding event, favored I. fulva. As expected, I. fulva plants survived at much higher rates than I. brevicaulis plants. Backcross hybrids toward I. fulva (BCIF) also survived at significantly higher rates than the reciprocal backcross toward I. brevicaulis (BCIB). Survivorship of BCIB hybrids was strongly influenced by the presence of a number of introgressed I. fulva alleles located throughout the genome, while survivorship in the reciprocal BCIF hybrids was heavily influenced by two epistatically acting QTL of opposite effects. These results demonstrate the potential for adaptive trait introgression between these two species and may help to explain patterns of genetic variation observed in naturally occurring hybrid zones.     

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION IN LOUISIANA IRISES: HYBRID FITNESS IN NATURE

Negative epistasis in hybrid genomes commonly results in postzygotic isolation between divergent lineages. However, some genomic regions may be selectively neutral or adaptive in hybrids and thus may potentially cross species barriers. We examined postzygotic isolation between ecologically similar species of Louisiana Iris: Iris brevicaulis and I. fulva to determine the potential for adaptive introgression in nature. Line-cross analyses allowed us a general overview of the gene action responsible for fitness-related traits. We then used a QTL mapping approach to detect genomic regions responsible for variation in these traits. Although hybrid classes suffered reduced fitness for many traits, hybrid means were equivalent to at least one of the parental species in overall estimates of maternal and paternal fitness during the two years of the field study. The genetic architecture underlying the fitness-related traits varied across field site and year of the study, thus emphasizing the importance of the environment in determining the degree of postzygotic isolation and potential for introgression across natural hybrid zones.     

Genetic mapping of species boundaries in Louisiana irises using IRRE retrotransposon display markers

Genetic mapping studies provide insight into the pattern and extent of genetic incompatibilities affecting hybridization between closely related species. Genetic maps of two species of Louisiana Irises, Iris fulva and I. brevicaulis, were constructed from transposon-based molecular markers segregating in reciprocal backcross (BC1) interspecific hybrids and used to investigate genomic patterns of species barriers inhibiting introgression. Linkage mapping analyses indicated very little genetic incompatibility between I. fulva and I. brevicaulis in the form of map regions exhibiting transmission ratio distortion, and this was confirmed using a Bayesian multipoint mapping analysis. These results demonstrate the utility of transposon-based marker systems for genetic mapping studies of wild plant species and indicate that the genomes of I. fulva and I. brevicaulis are highly permeable to gene flow and introgression from one another via backcrossing.     

Fitness and genetic architecture of parent and hybrid willows in common gardens

Models of hybrid zone dynamics incorporate different patterns of hybrid fitness relative to parental species fitness. An important but understudied source of variation underlying these fitness differences is the environment. We investigated the performance of two willow species and their F1, F2, and backcross hybrids using a common-garden experiment with six replicated gardens that differed in soil moisture. Aboveground biomass, catkin production, seed production per catkin, and seed germination rate were significantly different among genetic classes. For aboveground biomass and catkin production, hybrids generally had intermediate or inferior performance compared to parent species. Salix eriocephala had the highest performance for all performance measures, but in two gardens F, plants had superior or equal performance for aboveground biomass and female catkin production. Salix eriocephala and backcrosses to S. eriocephala had the highest numbers of filled seeds per catkin and the highest estimates of total fitness in all gardens. Measures of filled seeds per catkin and germination rate tend to support the model of endogenous hybrid unfitness, and these two measures had major effects on estimates of total seed production per catkin. We also estimated how the two willow species differ genetically in these fitness measures using line cross analysis. We found a complex genetic architecture underlying the fitness differences between species that involved additive, dominance, and epistatic genetic effects for all fitness measures. The environment was important in the expression of these genetic differences, because the type of epistasis differed among the gardens for above-ground biomass and for female catkin production. These findings suggest that fine-scale environmental variation can have a significant impact on hybrid fitness in hybrid zones where parents and hybrids are widely interspersed.     

The genetic architecture of reproductive isolation in Louisiana irises: flowering phenology

Despite the potential importance of divergent reproductive phenologies as a barrier to gene flow, we know less about the genetics of this factor than we do about any other isolating barrier. Here, we report on the genetic architecture of divergent flowering phenologies that result in substantial reproductive isolation between the naturally hybridizing plant species Iris fulva and I. brevicaulis. I. fulva initiates and terminates flowering significantly earlier than I. brevicaulis. We examined line crosses of reciprocal F1 and backcross (BC1) hybrids and determined that flowering time was polygenic in nature. We further defined quantitative trait loci (QTL) that affect the initiation of flowering in each of these species. QTL analyses were performed separately for two different growing seasons in the greenhouse, as well as in two field plots where experimental plants were placed into nature. For BCIF hybrids (BC1 toward I. fulva), 14 of 17 detected QTL caused flowering to occur later in the season when I. brevicaulis alleles were present, while the remaining 3 caused flowering to occur earlier. In BCIB hybrids (BC1 toward I. brevicaulis), 11 of 15 detected QTL caused flowering to occur earlier in the season when introgressed I. fulva alleles were present, while the remaining 4 caused flowering to occur later. These ratios are consistent with expectations of selection (as opposed to drift) promoting flowering divergence in the evolutionary history of these species. Furthermore, epistatic interactions among the QTL also reflected the same trends, with the majority of epistatic effects causing later flowering than expected in BCIF hybrids and earlier flowering in BCIB hybrids. Overlapping QTL that influenced flowering time across all four habitat/treatment types were not detected, indicating that increasing the sample size of genotyped plants would likely increase the number of significant QTL found in this study.     

Reproductive interactions between hybridizing irises: analyses of pollen-tube growth and fertilization success

Pollen-tube growth and seed siring ability were measured in crosses between the Louisiana iris species Iris fulva and Iris hexagona and their F₁ and F₂ hybrids. Flowers of the parental species were pollinated with self, outcross intraspecific, and interspecific pollen. Pollen-tube lengths were similar for all three pollen types in I. fulva, but in I. hexagona interspecific pollen tubes were longer than intraspecific pollen tubes. Pollen-tube lengths also differed for F₁ and F₂ flowers pollinated with I. fulva, I. hexagona, and hybrid pollen. For both hybrid classes I. fulva pollen tubes were the shortest while pollen tubes from I. hexagona and hybrids grew the furthest. Mixtures of genetically marked pollen were used to assess the seed siring ability of intra- and interspecific pollen in the parental species by varying the proportion of each pollen type in a replacement series design. For both species, the observed proportions of hybrid seeds were lower than the expected based on the frequency of each pollen type in the mixtures across all treatments. Flowers of I. fulva produced less than 10% hybrid progeny even when 75% of the pollen applied to stigmas was derived from interspecific flowers. The frequency of hybrid seed formation was somewhat greater in I. hexagona, but was still significantly lower than expected across all mixture treatments. Seed set per fruit remained constant across the mixture treatments for both species, but in I. fulva fruit set decreased with an increase in the proportion of interspecific pollen. The data indicate that both pre- and postfertilization processes contribute to discrimination against hybrid seed formation.     

SHADE TOLERANCE AND ITS EFFECT ON THE SEGREGATION OF TWO SPECIES OF LOUISIANA IRIS AND THEIR HYBRIDS

Iris fulva Ker. Gawler and Iris hexagona Walter have overlapping geographic ranges in Louisiana. In areas of overlap hybrids are fairly common. Iris hexagona occupies the borders of freshwater marshes of southern Louisiana while I. fulva can be found farther north along edges of natural levees, canals and swamps. Where the natural levee penetrates the marsh, natural hybridization can occur between I. hexagona and I. fulva. It has been suggested that one principal explanation for the segregation of the two species is that I. fulva grows best in semishade and I. hexagona grows best in full sun. A greenhouse study was conducted using rhizomes collected from the field to test this hypothesis and determine the relative shade tolerance of two natural hybrid types. Iris fulva, I. hexagona, and the two hybrid taxa were grown under 0% (control), 50% (medium shade), and 80% (high shade) reduction of sunlight for 6 months and then harvested. Iris fulva was found to be more tolerant of shading than I. hexagona and the two hybrids. Further, I. fulva was found to grow as well in control as in medium shade. Both hybrid taxa were more shade tolerant than I. hexagona. Iris hexagona was greatly affected by all levels of shade. In general, the results suggest that these hybrids are intermediate to the parental taxa in terms of shade tolerance.     

Leaf physiology reflects environmental differences and cytoplasmic background in Ipomopsis (Polemoniaceae) hybrids

Natural hybridization can produce individuals that vary widely in fitness, depending upon the performance of particular genotypes in a given environment. In a hybrid zone with habitat heterogeneity, differences in physiological responses to abiotic conditions could influence the fitness and spatial distribution of hybrids and parental species. This study compared gas exchange physiology of Ipomopsis aggregata, I. tenuituba, and their natural hybrids in situ and assessed whether physiological differences were consistent with their native environmental conditions. We also produced reciprocal F2s in a greenhouse study to test for cytonuclear effects on water-use efficiency (WUE). The relative performance of natural hybrids and parentals was consistent with their native habitats: I. aggregata at the coolest, wettest locations had the lowest WUE, while hybrids from the most xeric sites had the highest WUE. In hybrids, the mechanism by which both natural and experimental hybrids achieved this high WUE depended on cytotype: those with I. tenuituba cytoplasm had reduced transpiration, while those with I. aggregata cytoplasm had an increased photosynthetic rate, consistent with patterns in the cytoplasmic parent. The high WUE in hybrids may contribute to their high survival in the dry center of the natural hybrid zone, consistent with environment-dependent models of hybrid zone dynamics.     

IRIS NELSONII (IRIDACEAE): ORIGIN AND GENETIC COMPOSITION OF A HOMOPLOID HYBRID SPECIES

Nuclear and chloroplast DNA variation was assayed for two populations of Louisiana irises (Bayou Teche and Young's Coulee) that demonstrated extreme morphological variation and for a sample of the putative hybrid species, Iris nelsonii. The genetic markers examined in this analysis were diagnostic for either Iris fulva. Iris hexagona, or Iris brevicaulis. These data demonstrated that the two morphologically variable populations were hybrid associations involving all three of these species and that all three of these species were involved in the origin of I. nelsonii. The distribution of genetic variation in I. nelsonii was significantly different from that present in either of the two hybrid populations. I. nelsonii demonstrated significantly fewer foreign markers than the two hybrid populations. This finding is in accord with the prediction that I. nelsonii is a hybrid species that has undergone stabilization with regard to genetic recombination and segregation. Although the genotypic makeup of I. nelsonii was significantly different from other parental and hybrid populations, individual plants from this species cannot be unequivocally differentiated from either I. fulva or certain hybrid genotypes. This reflects the paradoxical nature of genotypic variation in hybrid species. Thus, a hybrid species may include genotypes that overlap with both parental and contemporary hybrid populations. In the case of I. nelsonii it is necessary to utilize additional information (morphological, chromosomal, ecological) to identify plants belonging to this taxon. One hybrid population (Young's Coulee) is suggested as a paradigm for the progenitor population that gave rise to I. nelsonii.     

Ecophysiology of first and second generation hybrids in a natural plant hybrid zone

Hybrids between related species vary widely in relative fitness, and that fitness can depend upon the environment. We investigated aspects of physiology that might influence fitness patterns in a plant hybrid zone. Seeds of Ipomopsis aggregata, I. tenuituba, F1 hybrids, F2 hybrids, and offspring of crosses between natural hybrids were planted into the relatively mesic site of origin for I. aggregata and the drier site for natural hybrids. We measured rates of photosynthesis (A _max), transpiration (E), instantaneous (A/E) and long-term (δ¹³C) indices of water use efficiency (WUE), and leaf nitrogen and carbon. We also examined correlations of these traits with plant size. Photosynthetic rate and A/E were higher in vegetative than flowering plants. WUE varied between sites and years, but differences among genotypic classes were spatially and temporally consistent. Instantaneous WUE was higher for F1 hybrids than for the average of the parental species, thereby showing heterosis. There was no evidence of hybrid breakdown, as WUE was no different in the F2 than the average across the F1 and parental species. Nor did WUE depend on cross direction in producing F1 progeny. Carbon isotope discrimination revealed higher long-term water use efficiency in I. tenuituba than I. aggregata. Leaf nitrogen was higher in I. tenuituba than I. aggregata, and higher in offspring of natural hybrids than in the F2. Results indicate heterosis for water use efficiency, with no hybrid breakdown. Heterosis in WUE may help to explain the relatively high survival of both reciprocal F1 hybrids in dry sites within the natural hybrid zone.     

QTL mapping reveals the genetic architecture of loci affecting pre- and post-zygotic isolating barriers in Louisiana Iris

ABSTRACT: BACKGROUND: Hybridization among Louisiana Irises has been well established and the genetic architecture of reproductive isolation is known to affect the potential for and the directionality of introgression between taxa. Here we use co-dominant markers to identify regions where QTL are located both within and between backcross maps to compare the genetic architecture of reproductive isolation and fitness traits across treatments and years. RESULTS: QTL mapping was used to elucidate the genetic architecture of reproductive isolation between Iris fulva and Iris brevicaulis. Homologous co-dominant EST-SSR markers scored in two backcross populations between I. fulva and I. brevicaulis were used to generate genetic linkage maps. These were used as the framework for mapping QTL associated with variation in 11 phenotypic traits likely responsible for reproductive isolation and fitness. QTL were dispersed throughout the genome, with the exception of one region of a single linkage group (LG) where QTL for flowering time, sterility, and fruit production clustered. In most cases, homologous QTL were not identified in both backcross populations, however, homologous QTL for flowering time, number of growth points per rhizome, number of nodes per inflorescence, and number of flowers per node were identified on several linkage groups. CONCLUSION: Two different traits affecting reproductive isolation, flowering time and sterility, exhibit different genetic architectures, with numerous QTL across the Iris genome controlling flowering time and fewer, less distributed QTL affecting sterility. QTL for traits affecting fitness are largely distributed across the genome with occasional overlap, especially on LG 4, where several QTL increasing fitness and decreasing sterility cluster. Given the distribution and effect direction of QTL affecting reproductive isolation and fitness, we have predicted genomic regions where introgression may be more likely to occur (those regions associated with an increase in fitness and unlinked to loci controlling reproductive isolation) and those that are less likely to exhibit introgression (those regions linked to traits decreasing fitness and reproductive isolation).     

Are natural hybrids fit or unfit relative to their parents?

The process of natural hybridization may produce genotypes that establish new evolutionary lineages. However, many authors have concluded that natural hybridization is of little evolutionary importance because hybrids, in general, are unfit relative to their progenitors. Deciding between these alternative conclusions requires that fitness be measured for hybrid classes and parental species. Recent analyses have found that hybrids are not uniformly unfit, but rather are genotypic classes that possess lower, equivalent or higher levels of fitness relative to their parental taxa.     

Interpopulation hybrid breakdown maps to the mitochondrial genome

Hybrid breakdown, or outbreeding depression, is the loss of fitness observed in crosses between genetically divergent populations. The role of maternally inherited mitochondrial genomes in hybrid breakdown has not been widely examined. Using laboratory crosses of the marine copepod Tigriopus californicus, we report that the low fitness of F(3) hybrids is completely restored in the offspring of maternal backcrosses, where parental mitochondrial and nuclear genomic combinations are reassembled. Paternal backcrosses, which result in mismatched mitochondrial and nuclear genomes, fail to restore hybrid fitness. These results suggest that fitness loss in T. californicus hybrids is completely attributable to nuclear-mitochondrial genomic interactions. Analyses of ATP synthetic capacity in isolated mitochondria from hybrid and backcross animals found that reduced ATP synthesis in hybrids was also largely restored in backcrosses, again with maternal backcrosses outperforming paternal backcrosses. The strong fitness consequences of nuclear-mitochondrial interactions have important, and often overlooked, implications for evolutionary and conservation biology.     

POLLEN TRANSFER BY NATURAL HYBRIDS AND PARENTAL SPECIES IN AN IPOMOPSIS HYBRID ZONE

Models of hybrid zones differ in their assumptions about the relative fitnesses of hybrids and the parental species. These fitness relationships determine the form of selection across the hybrid zone and, along with gene flow, the evolutionary dynamics and eventual outcome of natural hybridization. We measured a component of fitness, export and receipt of pollen in single pollinator visits, for hybrids between the herbaceous plants Ipomopsis aggregata and I. tenuituba and for both parental species. In aviary experiments with captive hummingbirds, hybrid flowers outperformed flowers of both parental species by receiving more pollen on the stigma. Although hummingbirds were more effective at removing pollen from anthers of I. aggregata, hybrid flowers matched both parental species in the amount of pollen exported to stigmas of other flowers. These patterns of pollen transfer led to phenotypic stabilizing selection, during that stage of the life cycle, for a stigma position intermediate between that of the two species and to directional selection for exserted anthers. Pollen transfer between the species was high, with flowers of I. aggregata exporting pollen equally successfully to conspecific and I. tenuituba flowers. Although this study showed that natural hybrids enjoy the highest quality of pollinator visits, a previous study found that I. aggregata receives the highest quantity of pollinator visits. Thus, the relative fitness of hybrids changes over the life cycle. By combining the results of both studies, pollinator-mediated selection in this hybrid zone is predicted to be strong and directional, with hybrid fitness intermediate between that of the parental species.     

Repeat swimming performance and its implications for inferring the relative fitness of asexual hybrid dace (Pisces: Phoxinus) and their sexually reproducing parental species

While theories explaining the evolution and maintenance of sex are abundant, empirical data on the costs and benefits of asexual relative to sexual reproduction are less common. Asexually reproducing vertebrates, while few, provide a rare opportunity to measure differences in fitness between asexual and sexual species. All known asexually reproducing vertebrates are of hybrid origin, and hybrid disadvantage (i.e., reduced hybrid fitness) is thought to facilitate long-term coexistence between asexual and sexual species. We used repeat swimming performance as a proxy for fitness to compare the fitness of asexual hybrid dace (Pisces: Phoxinus) and their sexually reproducing parental species, finescale dace (Phoxinus neogaeus) and northern redbelly dace (Phoxinus eos). We tested the prediction that, given the widespread coexistence of these hybrid and parental dace, the parental species should show equivalent and perhaps superior repeat performance relative to hybrids. A repeat constant acceleration test (U(max)) was conducted at both acclimation temperature (16 °C) and at an elevated temperature (25 °C) to simulate the combined influence of a repeat swim and acute temperature change that fish might experience in the wild. The asexual hybrids performed more poorly than at least one of the parental species. There was a negative effect of temperature on repeat swimming performance in all fish, and the repeat performance of hybrids was more severely affected by temperature than that of finescale dace. No difference in the effect of temperature on repeat performance was detected between hybrids and northern redbelly dace. These results suggest that hybrids suffer physiological costs relative to the parentals or at least that the hybrids do not gain advantage from hybrid vigor, which probably contributes to the coexistence of asexual and sexual species in this system.