Evolutionary adaptation to different thermal environments via transcriptional regulation

Comparison of Gene Expression Between Resistant and Susceptible Families Against VPAHPND and Identification of Biomarkers Used for Resistance Evaluation in Litopenaeus vannamei

Acute hepatopancreatic necrosis disease (AHPND) has caused a heavy loss to shrimp aquaculture since its outbreak. Vibrio parahaemolyticus (VP_AHPND) is regarded as one of the main pathogens that caused AHPND in the Pacific white shrimp Litopenaeus vannamei. In order to learn more about the mechanism of resistance to AHPND, the resistant and susceptible shrimp families were obtained through genetic breeding, and comparative transcriptome approach was used to analyze the gene expression patterns between resistant and susceptible families. A total of 95 families were subjected to VP_AHPND challenge test, and significant variations in the resistance of these families were observed. Three pairs of resistant and susceptible families were selected for transcriptome sequencing. A total of 489 differentially expressed genes (DEGs) that presented in at least two pairwise comparisons were screened, including 196 DEGs highly expressed in the susceptible families and 293 DEGs in the resistant families. Among these DEGs, 16 genes demonstrated significant difference in all three pairwise comparisons. Gene set enrichment analysis (GSEA) of all 27,331 expressed genes indicated that some energy metabolism processes were enriched in the resistant families, while signal transduction and immune system were enriched in the susceptible families. A total of 32 DEGs were further confirmed in the offspring of the detected families, among which 19 genes were successfully verified. The identified genes in this study will be useful for clarifying the genetic mechanism of shrimp resistance against Vibrio and will further provide molecular markers for evaluating the disease resistance of shrimp in the breeding program.

Differential survivorship of congeneric ornamental fishes under forecasted climate changes are related to anaerobic potential

Two Amazonian closely related tetras – cardinal Paracheirodon axelrodi and green neon P. simulans – were artificially acclimatized to environmental chambers mimicking future climate change scenarios (mild, moderate and extreme), using a microcosm facility. P. simulans survived (100%) to all scenarios after 30 days exposure, while P. axelrodi presented decreasing survival percentages according to environmental severity. These differences may be the reflection of distinct natural acclimatization to microhabitats between the species, which differ in thermal conditions. Survival responses might be related to differences in relative gene expression of lactate dehydrogenase (Ldh), suggesting that P. axelrodi anaerobic potential is lower or non-existent compared to P. simulans, not tolerating long-term thermal challenges. Accordingly, increases in temperature and in CO₂ levels caused increases in energy demand and resulted in activation of the anaerobic pathway, as demonstrated by the higher enzyme levels measured in head and tail portions of both species. Sustained anaerobic glycolysis is possible when fish live in challenging environments (low oxygen or high temperature). Our results clearly show that P. simulans has a larger scope for survival to higher energy demands due to its increased anaerobic potential compared to P. axelrodi.

Phylogeography in Nassarius mud snails: Complex patterns in congeneric species

One major goal for phylogeographical studies is to elucidate respective roles of multiple evolutionary and ecological forces that shape the current distribution patterns. In marine and coastal ecosystems, it has been generated a common realization that species with enormous population size and pelagic larval stages can disperse across broad geographical scales, leading to weak or even no phylogeographical structure across large geographical scales. However, the violation of such realization has been frequently reported, and it remains largely unexplored on mechanisms responsible for various phylogeographical patterns observed in different species at varied geographical scales. Here, we used a species-rich genus Nassarius to assess and compare phylogeographical patterns in congeneric species, and discuss causes and consequences underlying varied phylogeographical patterns. Interestingly, we observed complex phylogeographical patterns both within single species and across multiple species, and multiple analyses showed varied levels of genetic heterogeneity among sites within and across species. Available evidence suggests that related species with similar biological characteristics may not be necessary to result in consistent phylogeographical patterns. Multiple factors, including larval ecology, interactions between dispersal and natural selection, and human activity-mediated dispersal, can partially explain the complex patterns observed in this study. Deep investigations should be performed on these factors, particularly their respective roles in determining evolutionary/ecological processes to form phylogeographical patterns in species with high dispersal capacities in marine and coastal ecosystems.

Genetic architecture of gene transcription in two Atlantic salmon (Salmo salar) populations

Gene expression regulation has an important role in short-term acclimation and long-term adaptation to changing environments. However, the genetic architecture of gene expression has received much less attention than that of traditional phenotypic traits. In this study, we used a 5 × 5 full-factorial breeding design within each of two Atlantic salmon (Salmo salar) populations to characterize the genetic architecture of gene transcription. The two populations (LaHave and Sebago) are being used for reintroduction efforts into Lake Ontario, Canada. We used high-throughput quantitative real-time PCR to measure gene transcription levels for 22 genes in muscle tissue of Atlantic salmon fry. We tested for population differences in gene transcription and partitioned the transcription variance into additive genetic, non-additive genetic and maternal effects within each population. Interestingly, average additive genetic effects for gene transcription were smaller than those reported for traditional phenotypic traits in salmonids, suggesting that the evolutionary potential of gene transcription is lower than that of traditional traits. Contrary to expectations for early life stage traits, maternal effects were small. In general, the LaHave population had higher additive genetic effects for gene transcription than the Sebago population had, indicating that the LaHave fish have a higher adaptive potential to respond to the novel selection pressures associated with reintroduction into a novel environment. This study highlights not only the profound variation in gene transcription possible among salmonid populations but also the among-population variation in the underlying genetic architecture of such traits.

A Cost-Effective Approach to Sequence Hundreds of Complete Mitochondrial Genomes

We present a cost-effective approach to sequence whole mitochondrial genomes for hundreds of individuals. Our approach uses small reaction volumes and unmodified (non-phosphorylated) barcoded adaptors to minimize reagent costs. We demonstrate our approach by sequencing 383 Fundulus sp. mitochondrial genomes (192 F. heteroclitus and 191 F. majalis). Prior to sequencing, we amplified the mitochondrial genomes using 4–5 custom-made, overlapping primer pairs, and sequencing was performed on an Illumina HiSeq 2500 platform. After removing low quality and short sequences, 2.9 million and 2.8 million reads were generated for F. heteroclitus and F. majalis respectively. Individual genomes were assembled for each species by mapping barcoded reads to a reference genome. For F. majalis, the reference genome was built de novo. On average, individual consensus sequences had high coverage: 61-fold for F. heteroclitus and 57-fold for F. majalis. The approach discussed in this paper is optimized for sequencing mitochondrial genomes on an Illumina platform. However, with the proper modifications, this approach could be easily applied to other small genomes and sequencing platforms.

Constraints, independence, and evolution of thermal plasticity: probing genetic architecture of long- and short-term thermal acclimation

Seasonal and daily thermal variation can limit species distributions because of physiological tolerances. Low temperatures are particularly challenging for ectotherms, which use both basal thermotolerance and acclimation, an adaptive plastic response, to mitigate thermal stress. Both basal thermotolerance and acclimation are thought to be important for local adaptation and persistence in the face of climate change. However, the evolutionary independence of basal and plastic tolerances remains unclear. Acclimation can occur over longer (seasonal) or shorter (hours to days) time scales, and the degree of mechanistic overlap is unresolved. Using a midlatitude population of Drosophila melanogaster, we show substantial heritable variation in both short- and long-term acclimation. Rapid cold hardening (short-term plasticity) and developmental acclimation (long-term plasticity) are positively correlated, suggesting shared mechanisms. However, there are independent components of these traits, because developmentally acclimated flies respond positively to short-term acclimation. A strong negative correlation between basal cold tolerance and developmental acclimation suggests that basal cold tolerance may constrain developmental acclimation, whereas a weaker negative correlation between basal cold tolerance and short-term acclimation suggests less constraint. Using genome-wide association mapping, we show the genetic architecture of rapid cold hardening and developmental acclimation responses are nonoverlapping at the SNP and corresponding gene level. However, genes associated with each trait share functional similarities, including genes involved in apoptosis and autophagy, cytoskeletal and membrane structural components, and ion binding and transport. These results indicate substantial opportunity for short-term and long-term acclimation responses to evolve separately from each other and for short-term acclimation to evolve separately from basal thermotolerance.

Gene expression variation in duplicate lactate dehydrogenase genes: do ecological species show distinct responses?

Lactate dehydrogenase (LDH) has been shown to play an important role in adaptation of several aquatic species to different habitats. The genomes of Daphnia pulex, a pond species, and Daphnia pulicaria, a lake inhabitant, encode two L-LDH enzymes, LDHA and LDHB. We estimated relative levels of Ldh gene expression in these two closely related species and their hybrids in four environmental settings, each characterized by one of two temperatures (10°C or 20°C), and one of two concentrations of dissolved oxygen (DO; 6.5–7 mg/l or 2–3 mg/l). We found that levels of LdhA expression were 4 to 48 times higher than LdhB expression (p<0.005) in all three groups (the two parental species and hybrids). Moreover, levels of LdhB expression differed significantly (p<0.05) between D. pulex and D. pulicaria, but neither species differed from the hybrid. Consistently higher expression of LdhA relative to LdhB in both species and the hybrid suggests that the two isozymes could be performing different functions. No significant differences in levels of gene expression were observed among the four combinations of temperature and dissolved oxygen (p>0.1). Given that Daphnia dwell in environments characterized by fluctuating conditions with long periods of low dissolved oxygen concentration, we suggest that these species could employ regulated metabolic depression to survive in such environments.

Body temperatures in dinosaurs: what can growth curves tell us?

To estimate the body temperature (BT) of seven dinosaurs Gillooly, Alleen, and Charnov (2006) used an equation that predicts BT from the body mass and maximum growth rate (MGR) with the latter preserved in ontogenetic growth trajectories (BT-equation). The results of these authors evidence inertial homeothermy in Dinosauria and suggest that, due to overheating, the maximum body size in Dinosauria was ultimately limited by BT. In this paper, I revisit this hypothesis of Gillooly, Alleen, and Charnov (2006). I first studied whether BTs derived from the BT-equation of today's crocodiles, birds and mammals are consistent with core temperatures of animals. Second, I applied the BT-equation to a larger number of dinosaurs than Gillooly, Alleen, and Charnov (2006) did. In particular, I estimated BT of Archaeopteryx (from two MGRs), ornithischians (two), theropods (three), prosauropods (three), and sauropods (nine). For extant species, the BT value estimated from the BT-equation was a poor estimate of an animal's core temperature. For birds, BT was always strongly overestimated and for crocodiles underestimated; for mammals the accuracy of BT was moderate. I argue that taxon-specific differences in the scaling of MGR (intercept and exponent of the regression line, log-log-transformed) and in the parameterization of the Arrhenius model both used in the BT-equation as well as ecological and evolutionary adaptations of species cause these inaccuracies. Irrespective of the found inaccuracy of BTs estimated from the BT-equation and contrary to the results of Gillooly, Alleen, and Charnov (2006) I found no increase in BT with increasing body mass across all dinosaurs (Sauropodomorpha, Sauropoda) studied. This observation questions that, due to overheating, the maximum size in Dinosauria was ultimately limited by BT. However, the general high inaccuracy of dinosaurian BTs derived from the BT-equation makes a reliable test of whether body size in dinosaurs was ultimately limited by overheating impossible.

Experimental approaches to evaluate the contributions of candidate protein-coding mutations to phenotypic evolution

Identifying mechanisms of molecular adaptation can provide important insights into the process of phenotypic evolution, but it can be exceedingly difficult to quantify the phenotypic effects of specific mutational changes. To verify the adaptive significance of genetically based changes in protein function, it is necessary to document functional differences between the products of derived and wild-type alleles and to demonstrate that such differences impinge on higher-level physiological processes (and ultimately, fitness). In the case of metabolic enzymes, this requires documenting in vivo differences in reaction rate that give rise to differences in flux through the pathway in which the enzymes function. These measured differences in pathway flux should then give rise to differences in cellular or systemic physiology that affect fitness-related variation in whole-organism performance. Efforts to establish these causal connections between genotype, phenotype, and fitness require experiments that carefully control for environmental variation and background genetic variation. Here, we discuss experimental approaches to evaluate the contributions of amino-acid mutations to adaptive phenotypic change. We discuss conceptual and methodological issues associated with in vitro and in vivo studies of protein function, and the evolutionary insights that can be gleaned from such studies. We also discuss the importance of isolating the effects of individual mutations to distinguish between positively selected substitutions that directly contribute to improvements in protein function versus positively selected, compensatory substitutions that mitigate negative pleiotropic effects of antecedent changes.

The effect of short-term hypoxic exposure on metabolic gene expression

The long-term effect of hypoxia is to decrease both the production and use of ATP and thus decrease the reliance on mitochondrial oxidative energy production. Yet, recent studies include more immediate affects of hypoxia on gene expression and these data suggest the maintenance of mitochondrial function. To better understand the short-term physiological response to hypoxia, we quantified metabolic mRNA expression in the heart ventricles and livers of the teleost fish Fundulus grandis exposed to partial oxygen pressure of 2.8 kPa (-13.5% air saturation).Twenty-eight individuals from a single population were exposed to hypoxia for 0, 4, 8, 12, 24, 48, and 96  hr. Liver and cardiac tissues were sampled from the same individuals at 0-48  hr. At 96  hr, only cardiac tissue was assayed. Gene expression was significantly different (ANOVA, P < 0.05) for 17 of 226 metabolic genes (7.5%) in cardiac tissue and for 20 of 256 (7.8%) metabolic genes in hepatic tissue. For the two tissues examined in this study, the maximum response occurred at different times. For cardiac tissue, using Dunnett's post hoc test, most of these significant differences occurred at 96  hr of exposure. For liver, all but one significant difference occurred at 4  hr. Surprisingly, too many (relative to random expectations) of the genes with significant increase in mRNA are involved in the oxidative phosphorylation pathway: 44% of the significant genes at 96  hr in the heart and 33% of the significant genes at 4  hr in the liver are involved in the oxidative phosphorylation pathway. These data indicate that there are tissue-specific differences in the timing of the response to hypoxia, yet both cardiac and hepatic tissues have increases in mRNA that code for enzyme in the oxidative phosphorylation pathway. If these changes in mRNA produce a similar change in protein, then these data suggest that the initial response to hypoxia involves an increase in the oxidative pathway potentially as a mechanism to maintain ATP production.

Altitudinal variation at duplicated β-globin genes in deer mice: effects of selection, recombination, and gene conversion

Spatially varying selection on a given polymorphism is expected to produce a localized peak in the between-population component of nucleotide diversity, and theory suggests that the chromosomal extent of elevated differentiation may be enhanced in cases where tandemly linked genes contribute to fitness variation. An intriguing example is provided by the tandemly duplicated β-globin genes of deer mice (Peromyscus maniculatus), which contribute to adaptive differentiation in blood-oxygen affinity between high- and low-altitude populations. Remarkably, the two β-globin genes segregate the same pair of functionally distinct alleles due to a history of interparalog gene conversion and alleles of the same functional type are in perfect coupling-phase linkage disequilibrium (LD). Here we report a multilocus analysis of nucleotide polymorphism and LD in highland and lowland mice with different genetic backgrounds at the β-globin genes. The analysis of haplotype structure revealed a paradoxical pattern whereby perfect LD between the two β-globin paralogs (which are separated by 16.2 kb) is maintained in spite of the fact that LD within both paralogs decays to background levels over physical distances of less than 1 kb. The survey of nucleotide polymorphism revealed that elevated levels of altitudinal differentiation at each of the β-globin genes drop away quite rapidly in the external flanking regions (upstream of the 5' paralog and downstream of the 3' paralog), but the level of differentiation remains unexpectedly high across the intergenic region. Observed patterns of diversity and haplotype structure are difficult to reconcile with expectations of a two-locus selection model with multiplicative fitness.

The effects of mitochondrial genotype on hypoxic survival and gene expression in a hybrid population of the killifish, Fundulus heteroclitus

The physiological link between oxygen availability and mitochondrial function is well established. However, whether or not fitness variation is associated with mitochondrial genotypes in the field remains a contested topic in evolutionary biology. In this study, we draw on a population of the teleost fish, Fundulus heteroclitus, where functionally distinct subspecies hybridize, likely as a result of past glacial events. We had two specific aims: (i) to determine the effect of mtDNA genotype on survivorship of male and female fish under hypoxic stress and (ii) to determine the effect of hypoxic stress, sex and mtDNA genotype on gene expression. We found an unexpected and highly significant effect of sex on survivorship under hypoxic conditions, but no significant effect of mtDNA genotype. Gene expression analyses revealed hundreds of transcripts differentially regulated by sex and hypoxia. Mitochondrial transcripts and other predicted pathways were among those influenced by hypoxic stress, and a transcript corresponding to the mtDNA control region was the most highly suppressed transcript under the conditions of hypoxia. An RT-PCR experiment on the control region was consistent with microarray results. Effects of mtDNA sequence variation on genome expression were limited; however, a potentially important epistasis between mtDNA sequence and expression of a nuclear-encoded mitochondrial translation protein was discovered. Overall, these results confirm that mitochondrial regulation is a major component of hypoxia tolerance and further suggest that purifying selection has been the predominant selective force on mitochondrial genomes in these two subspecies.

SNP identification, verification, and utility for population genetics in a non-model genus

Background

By targeting SNPs contained in both coding and non-coding areas of the genome, we are able to identify genetic differences and characterize genome-wide patterns of variation among individuals, populations and species. We investigated the utility of 454 sequencing and MassARRAY genotyping for population genetics in natural populations of the teleost, Fundulus heteroclitus as well as closely related Fundulus species (F. grandis, F. majalis and F. similis).

Results

We used 454 pyrosequencing and MassARRAY genotyping technology to identify and type 458 genome-wide SNPs and determine genetic differentiation within and between populations and species of Fundulus. Specifically, pyrosequencing identified 96 putative SNPs across coding and non-coding regions of the F. heteroclitus genome: 88.8% were verified as true SNPs with MassARRAY. Additionally, putative SNPs identified in F. heteroclitus EST sequences were verified in most (86.5%) F. heteroclitus individuals; fewer were genotyped in F. grandis (74.4%), F. majalis (72.9%), and F. similis (60.7%) individuals. SNPs were polymorphic and showed latitudinal clinal variation separating northern and southern populations and established isolation by distance in F. heteroclitus populations. In F. grandis, SNPs were less polymorphic but still established isolation by distance. Markers differentiated species and populations.

Conclusions

In total, these approaches were used to quickly determine differences within the Fundulus genome and provide markers for population genetic studies.

Genetic structure among 50 species of the northeastern Pacific rocky intertidal community

Comparing many species' population genetic patterns across the same seascape can identify species with different levels of structure, and suggest hypotheses about the processes that cause such variation for species in the same ecosystem. This comparative approach helps focus on geographic barriers and selective or demographic processes that define genetic connectivity on an ecosystem scale, the understanding of which is particularly important for large-scale management efforts. Moreover, a multispecies dataset has great statistical advantages over single-species studies, lending explanatory power in an effort to uncover the mechanisms driving population structure. Here, we analyze a 50-species dataset of Pacific nearshore invertebrates with the aim of discovering the most influential structuring factors along the Pacific coast of North America. We collected cytochrome c oxidase I (COI) mtDNA data from populations of 34 species of marine invertebrates sampled coarsely at four coastal locations in California, Oregon, and Alaska, and added published data from 16 additional species. All nine species with non-pelagic development have strong genetic structure. For the 41 species with pelagic development, 13 show significant genetic differentiation, nine of which show striking FST levels of 0.1-0.6. Finer scale geographic investigations show unexpected regional patterns of genetic change near Cape Mendocino in northern California for five of the six species tested. The region between Oregon and Alaska is a second focus of intraspecific genetic change, showing differentiation in half the species tested. Across regions, strong genetic subdivision occurs more often than expected in mid-to-high intertidal species, a result that may reflect reduced gene flow due to natural selection along coastal environmental gradients. Finally, the results highlight the importance of making primary research accessible to policymakers, as unexpected barriers to marine dispersal break the coast into separate demographic zones that may require their own management plans.

Comparative characterization of a temperature responsive gene (lactate dehydrogenase-B, ldh-b) in two congeneric tropical fish, Lates calcarifer and Lates niloticus

The characterization of candidate loci is a critical step in obtaining insight into adaptation and acclimation of organisms. In this study of two non-model tropical (to sub-tropical) congeneric perciformes (Lates calcarifer and Lates niloticus) we characterized both coding and non-coding regions of lactate dehydrogenase-B (ldh-b), a locus which exhibits temperature-adaptive differences among temperate and sub-tropical populations of the North American killifish Fundulus heteroclitus. Ldh-b was 5,004 and 3,527 bp in length in L. calcarifer and L. niloticus, respectively, with coding regions comprising 1,005 bp in both species. A high level of sequence homology existed between species for both coding and non-coding regions of ldh-b (> 97% homology), corresponding to a 98.5% amino acid sequence homology. All six known functional sites within the encoded protein sequence (LDH-B) were conserved between the two Lates species. Ten simple sequence repeat (SSR) motifs (mono-, di-, tri- and tetranucleotide) and thirty putative microRNA elements (miRNAs) were identified within introns 1, 2, 5 and 6 of both Lates species. Five single nucleotide polymorphisms (SNPs) were also identified within miRNA containing intron regions. Such SNPs are implicated in several complex human conditions and/or diseases (as demonstrated by extensive genome-wide association studies). This novel characterization serves as a platform to further examine how non-model species may respond to changes in their native temperatures, which are expected to increase by up to 6 degrees C over the next century.

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.

Neutral and adaptive variation in gene expression

Variation among populations in gene expression should be related to the accumulation of random-neutral changes and evolution by natural selection. The following evolutionary analysis has general applicability to biological and medical science because it accounts for genetic relatedness and identifies patterns of expression variation that are affected by natural selection. To identify genes evolving by natural selection, we allocate the maximum among-population variation to genetic distance and then examine the remaining variation relative to a hypothesized important ecological parameter (temperature). These analyses measure the expression of metabolic genes in common-gardened populations of the fish Fundulus heteroclitus whose habitat is distributed along a steep thermal gradient. Although much of the variation in gene expression fits a null model of neutral drift, the variation in expression for 22% of the genes that regress with habitat temperature was far greater than could be accounted for by genetic distance alone. The most parsimonious explanation for among-population variation for these genes is evolution by natural selection. In addition, many metabolic genes have patterns of variation incongruent with neutral evolution: They have too much or too little variation. These patterns of biological variation in expression may reflect important physiological or ecological functions.

Adaptation and acclimation of photosynthetic microorganisms to permanently cold environments

Persistently cold environments constitute one of our world's largest ecosystems, and microorganisms dominate the biomass and metabolic activity in these extreme environments. The stress of low temperatures on life is exacerbated in organisms that rely on photoautrophic production of organic carbon and energy sources. Phototrophic organisms must coordinate temperature-independent reactions of light absorption and photochemistry with temperature-dependent processes of electron transport and utilization of energy sources through growth and metabolism. Despite this conundrum, phototrophic microorganisms thrive in all cold ecosystems described and (together with chemoautrophs) provide the base of autotrophic production in low-temperature food webs. Psychrophilic (organisms with a requirement for low growth temperatures) and psychrotolerant (organisms tolerant of low growth temperatures) photoautotrophs rely on low-temperature acclimative and adaptive strategies that have been described for other low-temperature-adapted heterotrophic organisms, such as cold-active proteins and maintenance of membrane fluidity. In addition, photoautrophic organisms possess other strategies to balance the absorption of light and the transduction of light energy to stored chemical energy products (NADPH and ATP) with downstream consumption of photosynthetically derived energy products at low temperatures. Lastly, differential adaptive and acclimative mechanisms exist in phototrophic microorganisms residing in low-temperature environments that are exposed to constant low-light environments versus high-light- and high-UV-exposed phototrophic assemblages.

The effects of diet and physiological stress on the evolutionary dynamics of an enzyme polymorphism

In the northern acorn barnacle Semibalanus balanoides, polymorphism at the mannose-6-phosphate isomerase (Mpi) locus appears to be maintained by distinct selection regimes that vary between intertidal microhabitats. The goal of the present experiment was to elucidate the mechanism of selection at the Mpi locus by examining the relationship between genotype and fitness-related life-history traits in laboratory manipulations. When barnacles were cultured on a mannose-supplemented diet and exposed to thermal stress, different Mpi genotypes exhibited differences in the rate of growth that predicted survivorship. In contrast, no such relationship was observed in control or fructose-supplemented dietary treatments either in the presence or in the absence of stress. Similarly, the phenotype and survivorship of genotypes at another allozyme locus and a presumably neutral mitochondrial DNA marker were homogeneous across all treatments and unaffected by experimental manipulations. These results suggest that the differential survivorship of Mpi genotypes in the field and laboratory results from a differential ability to process mannose-6-phosphate through glycolysis. The widespread polymorphism at Mpi observed in marine taxa may reflect the interaction between dietary composition and environmental heterogeneity in intertidal habitats.

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.

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.