Evidence of lactate dehydrogenase-B allozyme effects in the teleost, Fundulus heteroclitus

Mitochondrial Ecophysiology: Assessing the Evolutionary Forces That Shape Mitochondrial Variation

The mitonuclear species concept hypothesizes that incompatibilities between interacting gene products of the nuclear and mitochondrial genomes are a major factor establishing and maintaining species boundaries. However, most of the data available to test this concept come from studies of genetic variation in mitochondrial DNA, and clines in the mitochondrial genome across contact zones can be produced by a variety of forces. Here, we show that using a combination of population genomic analyses of the nuclear and mitochondrial genomes and studies of mitochondrial function can provide insight into the relative roles of neutral processes, adaptive evolution, and mitonuclear incompatibility in establishing and maintaining mitochondrial clines, using Atlantic killifish (Fundulus heteroclitus) as a case study. There is strong evidence for a role of secondary contact following the last glaciation in shaping a steep mitochondrial cline across a contact zone between northern and southern subspecies of killifish, but there is also evidence for a role of adaptive evolution in driving differentiation between the subspecies in a variety of traits from the level of the whole organism to the level of mitochondrial function. In addition, studies are beginning to address the potential for mitonuclear incompatibilities in admixed populations. However, population genomic studies have failed to detect evidence for a strong and pervasive influence of mitonuclear incompatibilities, and we suggest that polygenic selection may be responsible for the complex patterns observed. This case study demonstrates that multiple forces can act together in shaping mitochondrial clines, and illustrates the challenge of disentangling their relative roles.

Intraspecific variation and plasticity in mitochondrial oxygen binding affinity as a response to environmental temperature

Mitochondrial function has been suggested to underlie constraints on whole-organism aerobic performance and associated hypoxia and thermal tolerance limits, but most studies have focused on measures of maximum mitochondrial capacity. Here we investigated whether variation in mitochondrial oxygen kinetics could contribute to local adaptation and plasticity in response to temperature using two subspecies of the Atlantic killifish (Fundulus heteroclitus) acclimated to a range of temperatures (5, 15, and 33 °C). The southern subspecies of F. heteroclitus, which has superior thermal and hypoxia tolerances compared to the northern subspecies, exhibited lower mitochondrial O₂ P50 (higher O₂ affinity). Acclimation to thermal extremes (5 or 33 °C) altered mitochondrial O₂ P50 in both subspecies consistent with the effects of thermal acclimation on whole-organism thermal tolerance limits. We also examined differences between subspecies and thermal acclimation effects on whole-blood Hb O₂-P50 to assess whether variation in oxygen delivery is involved in these responses. In contrast to the clear differences between subspecies in mitochondrial O₂-P50 there were no differences in whole-blood Hb-O₂ P50 between subspecies. Taken together these findings support a general role for mitochondrial oxygen kinetics in differentiating whole-organism aerobic performance and thus in influencing species responses to environmental change.

The Adaptive Cline at LDH (Lactate Dehydrogenase) in Killifish Fundulus heteroclitus Remains Stationary After 40 Years of Warming Estuaries

Since the 1970s, water temperatures along the Atlantic seaboard of the United States have risen by an average of 0.5 °C in summer months and 2.2 °C in winter months. In response, the distribution and abundance of several nearshore species have changed dramatically, but no study has attempted to document whether estuarine populations have evolved greater thermal tolerance. Here, we re-examine the classic latitudinal cline at lactate dehydrogenase (LDH) in the killifish Fundulus heteroclitus that was originally described by Dennis Powers and associates from samples collected between 1970 and 1972. Laboratory and field evidences indicated that northern and southern isozymes at muscle LDH are locally adapted to cold and warm temperatures, respectively. Despite the potential for evolutionary response at this adaptive locus, we detected no significant shift of the LDH cline from 20 to 30 F. heteroclitus collected at each of 13 locations between the early 1970s and 2010. We conclude that the microevolution of LDH-mediated thermal tolerance has not occurred, that shifts in alleles are too incremental to be distinguished from random processes, or that F. heteroclitus uses phenotypic and genetic mechanisms besides LDH to respond to warmer waters.

Phenotypic variation of the housefly, Musca domestica: amounts and patterns of wing shape asymmetry in wild populations and laboratory colonies

Musca domestica L. (Diptera: Muscidae) is a vector of a range variety of pathogens infecting humans and animals. During a year, housefly experiences serial population bottlenecks resulted in reduction of genetic diversity. Population structure has also been subjected to different selection regimes created by insect control programs and pest management. Both environmental and genetic disturbances can affect developmental stability, which is often reflected in morphological traits as asymmetry. Since developmental stability is of great adaptive importance, the aim of this study was to examine fluctuating asymmetry (FA), as a measure of developmental instability, in both wild populations and laboratory colonies of M. domestica. The amount and pattern of wing shape FA was compared among samples within each of two groups (laboratory and wild) and between groups. Firstly, the amount of FA does not differ significantly among samples within the group and neither does it differ between groups. Regarding the mean shape of FA, contrary to non-significant difference within the wild population group and among some colonies, the significant difference between groups was found. These results suggest that the laboratory colonies and wild samples differ in buffering mechanisms to perturbations during development. Hence, inbreeding and stochastic processes, mechanisms dominating in the laboratory-bred samples contributed to significant changes in FA of wing shape. Secondly, general patterns of left-right displacements of landmarks across both studied sample groups are consistent. Observed consistent direction of FA implies high degrees of wing integration. Thus, our findings shed light on developmental buffering processes important for population persistence in the environmental change and genetic stress influence on M. domestica.

Geographic variation in Drosophila: From molecules to morphology and back

The examination of spatial variation can act as a substitute for temporal variation in population studies. Of particular use in the effort to understand the selective forces that govern spatial variation are the sibling species Drosophila melanogaster and D. simulans. Recent work at the level of both molecules and morphology has uncovered a great deal of spatial variation within and between these two species. Here we summarize these data, with reference to what they tell us about the history and nature of selection operating in these species.

Genomic approaches with natural fish populations

Natural populations v. inbred stocks provide a much richer resource for identifying the effects of nucleotide substitutions because natural populations have greater polymorphism. Additionally, natural populations offer an advantage over most common research organisms because they are subject to natural selection, and analyses of these adaptations can be used to identify biologically important changes. Among fishes, these analyses are enhanced by having a wide diversity of species (>28 000 species, more than any other group of vertebrates) living in a huge range of environments (from below freezing to > 46 degrees C, in fresh water to salinities >40 ppt.). Moreover, fishes exhibit many different life-history and reproductive strategies and have many different phenotypes and social structures. Although fishes provide numerous advantages over other vertebrate models, there is still a dearth of available genomic tools for fishes. Fishes make up approximately half of all known vertebrate species, yet <0.2% of fish species have significant genomic resources. Nonetheless, genomic approaches with fishes have provided some of the first measures of individual variation in gene expression and insights into environmental and ecological adaptations. Thus, genomic approaches with natural fish populations have the potential to revolutionize fundamental studies of diverse fish species that offer myriad ecological and evolutionary questions.

Classic approach revitalizes genomics: Complete characterization of a candidate gene for thermal adaptation in two coral reef fishes

Lactate dehydrogenase-B (ldh-b) encodes a metabolic enzyme (LDH-B) which plays an important role in maintaining aerobic performance and in thermal acclimation and/or adaptation of fish. As the first step in understanding the effect this enzyme has on the ability of tropical coral reef fishes to cope with thermal stress, we characterized both coding and non-coding regions of ldh-b in two congeneric perciformes, Plectropomus leopardus and Plectropomus laevis. Ldh-b was 4666 and 4539bp in length in P. leopardus and P. laevis, respectively, with coding regions comprising 1005bp in both species. We report a high level of sequence homology between the coding regions of ldh-b in these two species, with 98.1% identity of nucleotides corresponding to 100% amino acid identity between the deduced protein sequences. Comparison between non-coding (intron) regions of both species revealed the presence of several indels, despite the high level of homology observed (95.9% identity of intron nucleotides). Potential regulatory motifs and elements, including twenty-six simple sequence repeat motifs (mono-, di-, tri- and tetranucleotide) and twenty-three putative microRNA elements are identified within the introns of both species, further supporting recent demonstrations that such short motifs and elements exhibit widespread positioning throughout non-coding regions of the genome. This novel characterization of ldh-b in these two coral reef fishes allows for a wide range of future studies (e.g. analytical comparisons of ldh-b and LDH-B among different fish genera from different thermal environments and habitats).

Linking genotypes to phenotypes and fitness: how mechanistic biology can inform molecular ecology

The accessibility of new genomic resources, high-throughput molecular technologies and analytical approaches such as genome scans have made finding genes contributing to fitness variation in natural populations an increasingly feasible task. Once candidate genes are identified, we argue that it is necessary to take a mechanistic approach and work up through the levels of biological organization to fully understand the impacts of genetic variation at these candidate genes. We demonstrate how this approach provides testable hypotheses about the causal links among levels of biological organization, and assists in designing relevant experiments to test the effects of genetic variation on phenotype, whole-organism performance capabilities and fitness. We review some of the research programs that have incorporated mechanistic approaches when examining naturally occurring genetic and phenotypic variation and use these examples to highlight the value of developing a comprehensive understanding of the relationship between genotype and fitness. We give suggestions to guide future research aimed at uncovering and understanding the genetic basis of adaptation and argue that further integration of mechanistic approaches will help molecular ecologists better understand the evolution of natural populations.

Fundulus as the premier teleost model in environmental biology: opportunities for new insights using genomics

A strong foundation of basic and applied research documents that the estuarine fish Fundulus heteroclitus and related species are unique laboratory and field models for understanding how individuals and populations interact with their environment. In this paper we summarize an extensive body of work examining the adaptive responses of Fundulus species to environmental conditions, and describe how this research has contributed importantly to our understanding of physiology, gene regulation, toxicology, and ecological and evolutionary genetics of teleosts and other vertebrates. These explorations have reached a critical juncture at which advancement is hindered by the lack of genomic resources for these species. We suggest that a more complete genomics toolbox for F. heteroclitus and related species will permit researchers to exploit the power of this model organism to rapidly advance our understanding of fundamental biological and pathological mechanisms among vertebrates, as well as ecological strategies and evolutionary processes common to all living organisms.

Evidence of genetic selection for growth in new recruits of a marine fish

Abstract A cohort of Diplodus sargus, a coastal marine fish abundant in the Mediterranean Sea, has been surveyed from its settlement following the pelagic larval stage up to 4 months of age, when the juveniles are moving to adult habitats in order to assess selective processes. We followed the mortality by looking at the decrease in population abundance and, simultaneously, the genetic structure using allozymes and the growth associated with each genotype to test for a relationship between genotype and phenotype. The recruitment survey demonstrated that 80% of individuals arrived within a single night and that they show very similar age providing a discrete pulse of new recruits that we followed for changes in survival and allele frequencies. After 4 months, there was a total mortality of 80.8%, with the disappearance of 181 of 224 fish that initially colonized the rocky barrier. The decrease in number followed a logarithmic model with a maximum decrease in the early period (first 30 days). The model derived from the 4 months of data demonstrates that most of the mortality in the cohort occurs over the first 120 days following settlement and the model predicted a final abundance of 10 individuals after 1 year. Within the same period of 4 months, we observed significant decrease in multilocus heterozygosity. Such a decrease in heterozygosity partly resulted from a purge of the Pgm-80* allele. Together with this major change in a natural population, an aquarium experiment demonstrated that individuals with Pgm-80* alleles show significantly lower growth than other new recruits. We propose that the decrease in frequency of Pgm-80* in the natural environment is the result of targeted predation that eliminates smaller individuals and therefore individuals bearing Pgm-80*. The potential metabolic effect as well as a scenario that could lead to the maintenance of polymorphism is discussed.

Small-scale clinal variation, genetic diversity and environmental heterogeneity in the marine gobies Pomatoschistus minutus and P. lozanoi (Gobiidae, Teleostei)

Genetic variation was assayed at 14 allozyme loci in estuarine, coastal and offshore samples of lozano's goby, Pomatoschistus lozanoi and the sand goby, P. minutus. Samples were taken from locations on the Belgian Continental Shelf and in the Schelde estuary with a range of environmental heterogeneity. We evaluate whether any differences in (1) the degree of genetic variation and (2) allele frequencies at the various loci exist within samples occurring in various habitats on the BCS and in the Schelde estuary. No significant differences in levels of genetic diversity were recorded between estuarine, coastal and offshore samples in either species. A temporally stable clinal gradient in allele frequencies at the two-allele locus GPI-A(*) was observed in P. lozanoi, differentiating the samples in an estuarine, coastal and offshore group. We suggest that these differences might be maintained by balancing selection at locus GPI-A(*).

Variation in gene expression in response to stress in two populations of Fundulus heteroclitus

We used differential display PCR to identify hepatic genes responsive to handling stress and genes that differ in expression between populations of a fish, Fundulus heteroclitus, from different thermal environments. Despite substantial inter-individual variation, we cloned 20 putatively stress-regulated bands from Northern fish, 10 of which had high similarity to genes of known function. We selected five of these genes for further analysis based on their known roles in the stress response. Three of these genes (glucokinase, serine-threonine kinase 10 and cRAF) were confirmed as stress-responsive using real-time PCR. These genes increased in expression in response to a 7-day chronic stress protocol in fish from the Southern population of F. heteroclitus, but did not change significantly in fish from the Northern population. These three genes also differed in expression between populations in control fish, suggesting a link between the response to chronic stress and inter-population differences in gene expression in unstressed laboratory-acclimated fish. Two genes that did not respond to stress (glycogen synthase kinase and warm acclimation-related protein (WAP)) also differed between populations. Expression of WAP was eight-fold higher in Southern than in Northern fish, consistent with a previously suggested role for this gene in thermal acclimation or adaptation in fish.

Environmental adaptations as windows on molecular evolution

Changes in gene regulation may play an important role in adaptive evolution, particularly during adaptation to a changing environment. However, little is known about the molecular mechanisms underlying adaptively significant variation in gene regulation. To address this question, we are using environmental adaptations in populations of a fish, Fundulus heteroclitus as a window into the molecular evolution of gene regulation. F. heteroclitus are found along the East Coast of North America, with populations distributed along a steep thermal gradient. At the extremes of the species range, populations have undergone local adaptation to their habitat temperatures. A variety of genes differ in their regulation between these populations. We have determined the mechanism responsible for changes in lactate dehydrogenase-B (Ldh-B) gene regulation. A limited number of mutations in the regulatory sequence of this gene result in changes in its expression. Both the phenotypic (increased LDH activity) and genotypic (changes in Ldh-B regulatory sequences) differences between populations have been shown to be affected by natural selection, rather than genetic drift. Therefore, even a small number of mutations within important regulatory sequences can provide evolutionarily significant variation and have an impact on environmental adaptation.

Differential performance among LDH-B genotypes in Rana lessonae tadpoles

The European pool frog, Rana lessonae, is widely polymorphic for two common alleles (b,e) at the lactate dehydrogenase-B (LDH-B) locus. We compared fitness-related larval life-history traits among LDH-B genotypes, which originated from segregation in heterozygous parents, in an artificial pond experiment where tadpoles of R. lessonae from a Swiss population were raised together with tadpoles of the hemiclonal hybrid R. esculenta at two densities. In R. lessonae, LDH-B e/e homozygotes at each density had a higher proportion of metamorphs among survivors, reached metamorphosis earlier, and were heavier at metamorphosis than b/b homozygotes; b/e heterozygotes had intermediate values. That e/e individuals were superior to b/b in both time to and mass at metamorphosis is surprising because these two life-history traits are thought to reflect a performance trade-off; e/e genotypes apparently compensated for shorter time to metamorphosis by a higher growth rate. The two alleles showed the same performance ranking when combined in hybrids with a R. ridibunda allele: When R. esculenta from Swiss populations reared in the same ponds had received the e allele rather than the b allele from their R. lessonae parent, they reached metamorphosis earlier, but did not differ in mass at metamorphosis. The degree of linkage disequilibrium in the source population of the eight R. lessonae used as parents of the R. lessonae tadpoles is unknown, so we cannot exclude the possibility that the performance differences are caused by some anonymous tightly linked gene, rather than the LDH-B locus, that constitutes the genomically localized target of natural selection. A causal involvement of LDH-B is plausible, nevertheless, because this enzyme takes part in the central energy-metabolizing processes and has been reported to underlie fitness differences in other animals; also, differential performance of LDH-B genotypes has been observed in R. lessonae larvae from another population. The present results suggest strong directional selection for allele e; the sum of available data, including an independent laboratory experiment, suggests that partial environment-dependent overdominance combined with balancing selection favoring e/e homozygotes under some and b/b homozygotes under other conditions may be partially responsible for the broad maintenance of the LDH-B polymorphism in R. lessonae.

Adaptive variation in lactate dehydrogenase-B gene expression: role of a stress-responsive regulatory element

Although changes in gene regulation may play an important role in adaptive evolution, there have been few attempts to investigate the molecular mechanisms responsible for adaptively significant variation in gene expression. Here we describe the mechanism underlying an adaptive difference in the expression of the lactate dehydrogenase-B gene (Ldh-B) between northern and southern populations of the fish Fundulus heteroclitus. Ldh-B regulatory sequences from northern and southern individuals, coupled to a luciferase reporter gene, were introduced into the livers of live fish. Deletion studies indicated that sequence changes between 400 and 500 bp upstream of the transcription start site resulted in a 2-fold difference in reporter gene transcription. These sequence changes can account for the previously observed 2-fold difference in Ldh-B transcription between populations. Variation in transcription factors did not play an important role. Sequences within the functionally important region resemble a mammary tumor virus glucocorticoid responsive element (MTV-GRE) in southern alleles, whereas northern alleles differ from the consensus by 1 bp. To test the hypothesis that this element is involved in the variation between populations of F. heteroclitus, we exposed transiently transgenic fish containing Ldh-B regulatory sequence/reporter gene constructs to handling stress or injected cortisol. Both treatments increased reporter gene transcription driven by southern alleles but not northern alleles, as expected if an MTV-GRE sequence were involved. This finding suggests that sequence variation in a GRE is the cause of the adaptive differences in Ldh-B gene expression between populations and demonstrates that small changes in gene regulatory sequences can have important evolutionary consequences.

Do discrepancies between microsatellite and allozyme variation reveal differential selection between sea and lagoon in the sea bass (Dicentrarchus labrax)?

In the present study the genetic structure of Dicentrarchus labrax (14 samples from the Mediterranean) was analysed at six microsatellite loci, in order to test the hypothesis that some enzymatic loci undergo selection between marine and lagoon habitat. Eight of the 14 samples were analysed at both microsatellite and allozyme markers. The analysis of the genetic variation among the Mediterranean samples showed that (i) &Fcirc;ST values obtained with the six microsatellite loci were much smaller than those obtained with the 28 allozymes and (ii) microsatellite loci seemed to reflect more the geographical proximity than an ecological one. Thirteen enzymatic loci exhibited moderate to high values compared with microsatellites. This was interpreted as evidence that these allozymes are non-neutral. However, only six loci seemed to be implicated in differentiation between marine and lagoon samples, the causes of selection being unknown for the others. A possible scenario of population dynamics of the sea bass between marine and lagoon habitat is suggested.

Effects of urea on M4-lactate dehydrogenase from elasmobranchs and urea-accumulating Australian desert frogs

We measured the effect of urea on M4-lactate dehydrogenase (M4-LDH) from elasmobranchs and Australian desert frogs (urea accumulators) and from two animals that do not accumulate urea, the axolotl and the rabbit. An analysis of the effect of urea on the Kd(NADH), V, V/K(m(prr)) and V/K(m(NADH)) shows that in all cases the major effect of urea was on the binding of pyruvate, which fits with data in the literature that show that urea acts as a competitive inhibitor of LDH. The characteristics of the elasmobranch enzymes are consistent with a proposed adaptation model, but the situation for the enzymes from the aestivating frogs is equivocal. Urea (400 mM) had less effect on the K(m(prr)) of M4-LDH from the urea accumulators than it did on the non-accumulators, suggesting a general adaptation and that the enzyme produced by the aestivating frogs (urea accumulators) is kinetically different from that of non-aestivating frogs (non-accumulators). A new approach is used to characterize the overall pattern of adaptation to urea. The pattern is similar in an enzyme from an elasmobranch and an aestivating frog despite the temporary presence of urea in the latter and the phylogenetic difference between these animals.

Qualitative and quantitative strategies of thermal adaptation of grass carp (Ctenopharyngodon idella) cytoplasmic malate dehydrogenases

Malate dehydrogenase isozymes of grass carp,Ctenopharyngodon idella, were identified by mitochondrial preparation and thermal denaturation. The structural and kinetic characteristics of chromatographically separated thermostable and thermolabile cMDHs were different in (1) half-life at 42°C, 10 min and 24 h, respectively, (2) optimal substrate, oxaloacetate and malate, concentrations, and (3) the apparent Michaelis-Menten constants of NADH and oxaloacetate.Total MDH activity in white muscle of 11°C-acclimated fish was about twice that of the 30°C-acclimated group. In addition, the ratio of the thermostable to thermolabile cMDH activity in white muscle of 30°C-acclimated fish was significantly higher than that of 11°C-acclimated fish. These results suggest that temperature acclimation can induce temperature compensation in MDH activity and differential expression of thermostable and thermolabile cMDH isozymes in freshwater fish.

Rapid enzyme assays investigating the variation in the glycolytic pathway in field-caught populations of Fundulus heteroclitus

Variation in enzyme expression may be important in evolutionary adaptation, yet is seldom studied. Furthermore, no studies have examined the expression of all enzymes in a defined metabolic pathway. Enzyme concentration is a measure of enzyme expression and was ascertained by assaying maximal activity. Presented here is an analysis of variation of maximal enzyme activity for all the enzymes in a single metabolic pathway, glycolysis, from three clinically distributed populations of the fish, Fundulus heteroclitus. Techniques for rapidly analyzing maximal enzyme activity for all the enzymes of an entire metabolic pathway from many individuals are described. The high degree of repeatability (mean coefficient of variation for replicates, 4.4%) and sensitivity (less than 3 mg of tissue is required to measure all 10 enzymes) of these assays demonstrate the utility of such an approach for analyzing variation among populations for a large numbers of enzymes. Results from these studies indicate that (1) the average coefficient of variation for all enzyme determinations within a population is 45.3% and (2) between populations, the activity of 5 of the 10 glycolytic enzymes are significantly different. This considerable variation occurs even in populations where there is little allelic variation. These data demonstrating substantial variation in enzyme expression support the idea that changes in gene regulation may be as important as, or even more important than, changes in biochemical kinetic parameters in evolutionary processes.

HIERARCHICAL GENETIC STRUCTURE AND GENE FLOW IN MACROGEOGRAPHIC POPULATIONS OF THE EASTERN TENT CATERPILLAR (MALACOSOMA AMERICANUM)

Genetic structure and inferred rates of gene flow in macrogeographic populations of the eastern tent caterpillar Malacosoma americanum were analyzed at two hierarchical scales: local demes and regional subpopulations. Wright's F-statistics were used to estimate population genetic structure using multilocus genotypic data generated electrophoretically. Estimated values of F_ST and the distribution of private alleles were then used to obtain indirect estimates of gene flow. We found modest, though significant, genetic structure at both spatial scales, a pattern consistent with high rates of gene flow over the large distances involved. Modest values obtained for Nei's genetic distance also suggested high levels of gene flow across the range of this species, although some gene-flow restriction resulting from isolation by distance was suggested by a positive regression of genetic distance on geographic distance. The observed homogeneity at enzyme loci across the range of M. americanum parallels the reported uniformity in morphology, suggesting a general absence of local genetic differentiation in this widely distributed species. The genetic homogeneity observed in this wide-ranging insect is discussed in terms of organism-specific environmental experience at different spatial scales. Some organisms occupying apparently heterogeneous environments may ameliorate unsuitable local conditions through microhabitat selection or behavioral modification of their microenvironment. This may be accomplished in M. americanum through group shelter construction and behavioral thermoregulation, closely tying thermoregulation to social biology in this species. If in this way the tent helps produce an effectively homogeneous environment for this species across its extensive range, this system may provide a unique example of how social behavior can influence the distribution of genetic variation in a population.

Regulation of fitness in yeast overexpressing glycolytic enzymes: responses to heat shock and nitrogen starvation

Current models based on the analysis of linear metabolic pathways at steady-state predict that large increases over wild type in the activity of one enzyme will not alter an organism's fitness. This prediction is tested at steps in a highly branched pathway under two conditions known to alter steady-state: heat shock and nitrogen starvation. Saccharomyces cerevisiae transformants overproducing 1 of 4 enzymes in glycolysis (hexokinase B, phosphoglucose isomerase, phosphofructokinase, or pyruvate kinase) were subjected to heat shock in both exponential and stationary phases of growth. In neither phase does enzyme overexpression alter heat shock sensitivity. When starved for nitrogen in acetate medium, transformants overproducing hexokinase, phosphoglucose isomerase, and phosphofructokinase sporulate at the same rate and with the same frequency as cells harbouring only the plasmid vector. Current models therefore correctly predict the relationship between activity and components of fitness for 3 of 4 enzymes. By contrast, cells overexpressing pyruvate kinase sporulate poorly. This defect is not observed among cells transformed with a plasmid containing a Tn5 disrupted copy of the PYK gene. These findings are consistent with reports that implicate the PYK locus in yeast cell cycle control and suggest that it may be challenging to model relations between fitness and activity for multifunctional proteins.