Ecological Genetics of Mpi and Gpi Polymorphisms in the Acorn Barnacle and the Spatial Scale of Neutral and Non-neutral Variation

Different allozyme genotypes at the mannose phosphate isomerase(Mpi) locus in the northern acorn barnacle (Semibalanus balanoides)show a strong association with distinct intertidal microhabitats.In estuaries along the Maine Coast, the FF homozygote has higherfitness in exposed, high-tide level microhabitats while theSS homozygote has higher fitness under algal cover or at low-tidemicrohabitats. These patterns are consistent with a Levene (1953)model of balancing selection. In these same samples, polymorphismsat the glucose phosphate isomerase locus (Gpi) and mitochondrialDNA (mtDNA) show no fitness differences among microhabitats,providing intra-genomic controls supporting selection at ornear Mpi. Here we report a similar analysis of genotype-by-microhabitatassociations at sites in Narragansett Bay, Rhode Island, closeto the southern range limit of S. balanoides. Genotype zonationat Mpi between high- and low-tide microhabitats is significantlydifferent between Maine and Narragansett Bay due to oppositezonation patterns for the SF and FF genotypes. Enzyme activitydata are consistent with this "reverse" zonation. At Gpi, thereis significant microhabitat zonation in Narragansett Bay, whilethis locus behaves as a neutral marker in Maine. Mt DNA showsno significant microhabitat zonation in either Rhode Islandor Maine. The Mpi data suggest that Levene-type selection foralternative genotypes in alternative habitats may operate atscales of both 10's of meters and 100's of kilometers. The Gpidata show how an apparently neutral locus can exhibit non-neutralvariation under different environmental conditions. We arguethat both Mpi and Gpi provide important genetic variation foradaptation to environmental heterogeneity that is recruitedunder distinct conditions of stress and carbohydrate substrateavailability.     

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.     

Environmental heterogeneity and balancing selection in the acorn barnacle Semibalanus balanoides

The northern acorn barnacle Semibalans banlanoides occupies several intertidal microhabitats which vary greatly in their degree of physical stress. This environmental heterogeneity creates distinct selection regimes which can maintain genetic variation in natural populations. Despite considerable attention placed on the link between spatial variation in fitness and balancing selection at specific loci, experimental manipulations and fitness estimates for molecular polymorphisms have rarely been conducted in the wild. The aim of this transplant experiment was to manipulate the level of physical stress experienced by a cohort of barnacles in the field and then investigate the spatial variation in fitness for genotypes at three loci: two candidate allozymes and the mitochondrial DNA control region. The viability of mannose-6-phosphate isomerase (Mpi) genotypes was dependent on the level of physical stress experienced in the various treatments; alternative homozygotes were favoured in alternative high stress-low stress environments. In contrast, the fitness of genotypes at other loci was equivalent among treatments and unaffected by the manipulation. Evaluated in the light of balancing selection models, these data indicate that the presence of multiple environmental niches is sufficient to promote a stable Mpi polymorphism in barnacle populations and that allelic variation at this locus reflects the process of adaptation to the heterogeneous intertidal landscape.     

Genetic heterogeneity among intertidal habitats in the flat periwinkle, Littorina obtusata

Comparisons among patterns exhibited by functionally distinct genetic markers have been widely used to infer the impacts of demography and selection in structuring genetic variation in natural populations. However, such multilocus comparisons remain an indirect evaluation of selection at particular candidate loci; ideally, the identification of a candidate gene by comparative genetic methodologies should be complemented by functional analyses and experimental manipulations of genotypes in the laboratory or field. We examined genotype frequency variation among replicated intertidal habitats at two spatial scales in the grazing snail Littorina obtusata. Both of the candidate allozyme markers varied predictably with environment, and these patterns were consistent at both spatial scales. Three of four reference loci were spatially homogeneous, but one microsatellite exhibited significant structure at both geographical and mesoscales. To initiate a direct examination of whether the observed genotype frequency variation at one of the candidate markers, mannose-6-phosphate isomerase (MPI), was impacted by differential survivorship of genotypes, we conducted a series of laboratory-based thermal stress assays using snails from two geographically disparate source populations. When snails were exposed to bouts of thermal/desiccation stress, patterns of mortality were nonrandom with respect to MPI genotype. Furthermore, patterns of mortality in the laboratory manipulation coincided with the observed distribution of genotypes in the field. The data suggest the operation of selection at the Mpi or a linked locus, but functional studies and further experimentation are required to establish the relationship between MPI genotype and fitness across heterogeneous intertidal environments.     

Genetic variation in the acorn barnacle from allozymes to population genomics

Understanding the patterns of genetic variation within and among populations is a central problem in population and evolutionary genetics. We examine this question in the acorn barnacle, Semibalanus balanoides, in which the allozyme loci Mpi and Gpi have been implicated in balancing selection due to varying selective pressures at different spatial scales. We review the patterns of genetic variation at the Mpi locus, compare this to levels of population differentiation at mtDNA and microsatellites, and place these data in the context of genome-wide variation from high-throughput sequencing of population samples spanning the North Atlantic. Despite considerable geographic variation in the patterns of selection at the Mpi allozyme, this locus shows rather low levels of population differentiation at ecological and trans-oceanic scales (F(ST)?～?5%). Pooled population sequencing was performed on samples from Rhode Island (RI), Maine (ME), and Southwold, England (UK). Analysis of more than 650 million reads identified approximately 335,000 high-quality SNPs in 19 million base pairs of the S. balanoides genome. Much variation is shared across the Atlantic, but there are significant examples of strong population differentiation among samples from RI, ME, and UK. An F(ST) outlier screen of more than 22,000 contigs provided a genome-wide context for interpretation of earlier studies on allozymes, mtDNA, and microsatellites. F(ST) values for allozymes, mtDNA and microsatellites are close to the genome-wide average for random SNPs, with the exception of the trans-Atlantic F(ST) for mtDNA. The majority of F(ST) outliers were unique between individual pairs of populations, but some genes show shared patterns of excess differentiation. These data indicate that gene flow is high, that selection is strong on a subset of genes, and that a variety of genes are experiencing diversifying selection at large spatial scales. This survey of polymorphism in S. balanoides provides a number of genomic tools that promise to make this a powerful model for ecological genomics of the rocky intertidal.     

INTERTIDAL MICROHABITAT AND SELECTION AT MPI: INTERLOCUS CONTRASTS IN THE NORTHERN ACORN BARNACLE,SEMIBALANUS BALANOIDES

Barnacles were sampled from various microhabitats in the rocky intertidal at multiple sites in two years. At sites in which there were large differences among microhabitats in temperature profiles, Mpi genotype frequencies were consistently and significantly different. Genotype frequencies for another allozyme locus (Gpi) as well as a DNA marker shown to be neutral (the mtDNA control region) were statistically homogeneous among thermal microhabitats at all sites in both years. The data indicate that temperature and/or desiccation mediated selection is operating at Mpi or a linked locus and that Mpi genotypes experience differential mortality in the various habitat types. If the relative fitness of genotypes is dependent on habitat type, the Mpi polymorphism may be actively maintained by a Levene model of balancing selection (Levene 1953). Because barnacle larvae are produced in abundance each year and spend five to eight weeks dispersing in the water column, there is little opportunity for the accumulation of adaptive divergence over the environmental grain size relevant in intertidal habitats. The Mpi polymorphism may be an important component of a suite of traits involved in the adaptation of barnacles to heterogeneous environments. Due to the relatively high concentration of mannose in a variety of algal groups, the metabolism of mannose may substantially affect individual performance and fitness in this and other species that feed on algae and phytoplankton. Because the Mpi locus is one of the most strongly polymorphic in marine organisms, these findings may be relevant for a diversity of other such species.     

A fundamental relationship between genotype frequencies and fitnesses

The set of possible postselection genotype frequencies in an infinite, randomly mating population is found. Geometric mean heterozygote frequency divided by geometric mean homozygote frequency equals two times the geometric mean heterozygote fitness divided by geometric mean homozygote fitness. The ratio of genotype frequencies provides a measure of genetic variation that is independent of allele frequencies. When this ratio does not equal two, either selection or population structure is present. Within-population HapMap data show population-specific patterns, while pooled data show an excess of homozygotes.     

Using environmental effects on fecundity to compare the adaptive characteristics of the morphs in a hybrid zone of Littorina saxatilis

We studied the role of ecological processes in the maintenance of a hybrid zone for the intertidal snail Littorina saxatilis by comparing the fecundity responses of the females from different morphs to local environmental variation. This procedure makes it possible to discern if the environments maximizing fecundity are the same in these morphs, and therefore to study adaptive differences between them. We found that the pure morphs were adapted to their respective shore levels, thus confirming that this hybrid zone is associated with the occurrence of an ecotone between two habitats. However, fecundity differences did not explain the precise spatial location and width of the hybrid zone. We also found that the adaptive characteristics of hybrids were very different from those of the pure morphs, in sharp contrast with any expectations based on their intermediate phenotype. The comparison of fitness responses to environmental variation might be an useful procedure in the field study of adaptive differences between groups of individuals.     

Ecological and genetic associations across a Heliconius hybrid zone

Differences in habitat use can bridge early and late stages of speciation by initiating assortative mating. Heliconius colour pattern races might select habitats over which each pattern confers a relative fitness advantage because signal efficacy of wing patterns can vary by environment. Thus habitat preferences could serve to promote the evolution of mimetic colour patterns for mate choice. Here I compare colour pattern genotype and phenotype frequencies to environmental variation across the H. erato hydara x H. erato erato hybrid zone in French Guiana to determine whether races exhibit habitat preferences. I found that genotype and phenotype frequencies correspond to differences in land cover moreso than to other environmental factors. Temporal shifts in colour pattern genotypes, phenotypes and land cover also were associated at individual sample sites, which further suggests that H. erato races differ in habitat use and that habitat preferences may promote speciation among Heliconius butterflies.     

GxG epistasis in growth and condition and the maintenance of genetic polymorphism in Gambusia holbrooki

Theory on indirect genetic effects (IGEs) indicates that variation in the genetic composition of social groups can generate GxG epistasis that may promote the evolution of stable polymorphisms. Using a livebearing fish with a genetic polymorphism in coloration and associated behavioral differences, we tested whether genotypes of social partners interacted with focal individual genotypes to influence growth and condition over 16 weeks of development. We found that IGEs had a significant influence on patterns of feeding, regardless of focal fish genotype. There was no influence of social environment on juvenile length, but there was significant GxG epistasis for body condition. Each focal juvenile was in better condition when its own genotype was not present in adult social partners. These data are consistent with negative frequency‐dependent selection in which each morph performs better when it is rare. Neither variation in feeding nor activity‐related behaviors explained variation in body condition, suggesting that GxG epistasis for condition was caused by physiological differences between the two genotypes. These findings indicate that GxG epistasis in a given polymorphism can generate fitness landscapes that contribute to the maintenance of that polymorphism and to maintenance of genetic variation for additional fitness‐related traits.     

Regional variation in the spatial scale of selection at MPI* and GPI* in the acorn barnacle Semibalanus balanoides (Crustacea)

Elucidating the ecological processes by which adaptive genetic polymorphism is maintained in heterogeneous environments requires knowledge on the spatial scale at which alternate habitats affect genotype-specific fitness. The general objective of this study was to document patterns of temporal and spatial variation of genetic polymorphism in the acorn barnacle (Semibalanus balanoides) at MPI* and GPI* allozyme loci. A total of 7261 barnacles were sampled in the intertidal at various locations north and south of the Miramichi estuary, New Brunswick, Canada. The results of this study supported the hypothesis that both MPI* and GPI* are under the effect of strong directional selection south of the Miramichi, whereas neutrality cannot be ruled out at sampling sites located north of the estuary. Comparisons between this study and previous ones also question the generality of current hypotheses regarding ecological processes that are responsible for maintaining polymorphism at MPI* and GPI* in the acorn barnacle.     

The maintenance of mitochondrial genetic variation by negative frequency‐dependent selection

Mitochondrial genes generally show high levels of standing genetic variation, which is puzzling given the accumulating evidence for phenotypic effects of mitochondrial genetic variation. Negative frequency‐dependent selection, where the relative fitness of a genotype is inversely related to its frequency in a population, provides a potent and potentially general process that can maintain mitochondrial polymorphism. We assessed the change in mitochondrial haplotype frequencies over 10 generations of experimental evolution in 180 seed beetle populations in the laboratory, where haplotypes competed for propagation to subsequent generations. We found that haplotypes consistently increased in frequency when they were initially rare and decreased in frequency when initially common. Our results have important implications for the use of mtDNA haplotype frequency data to infer population level processes and they revive the general hypothesis that negative frequency‐dependent selection, presumably caused by habitat heterogeneity, may commonly promote polymorphism in ecologically relevant life history genes.     

Genetic variation associated with chronic water and nutrient stress in pinyon pine

We examined how genotypic variation at the glycerate dehydrogenase (GLY) locus in pinyon pine might be affected by environmental stress and herbivory. We compared GLY allelic and genotypic frequencies among mature and juvenile trees growing in stressful cinder soil and adjacent sandy-loam soil. While no association was found with herbivory, three lines of evidence support the hypothesis that GLY slow homozygotes (SS) are selected for under conditions of environmental stress. 1) We found no differences at the GLY locus between juvenile and mature trees growing in sandy-loam soil. 2) However, in the stressful cinder soil we found a lower frequency of SS homozygotes in juvenile trees than in mature trees, suggesting that selection has occurred in one generation. 3) Using biomass as a measure of fitness, SS homozygotes were superior to other GLY genotypes in cinder soil. Two lines of evidence were inconsistent with the hypothesis of selection for SS homozygotes in stressful environments. 1) SS homozygote frequency between habitats was not significantly different, suggesting that too few generations have passed for differentiation to occur, and/or insufficient selection pressure. 2) The low frequency of SS homozygote juvenile trees growing in cinder soil suggests that counterbalancing selection may prevent SS homozygotes from increasing in the population.     

Genetic Variation in a Heterogeneous Environment. I. Temporal Heterogeneity and the Absolute Dominance Model

The conditions for a stable polymorphism and the equilibrium gene frequency in an infinite population are compared when there is spatial or temporal environmental heterogeneity for the absolute dominance model. For temporal variation the conditions for stability are more restrictive and the equilibrium gene frequency is often at a low gene frequency. In a finite population, temporal environmental heterogeneity for the absolute dominance model was found to be quite ineffective in maintaining genetic variation and is often less effective than no selection at all. For comparison, the maximum maintenance for temporal variation is related to the overdominant model. In general, cyclic environmental variation was found to be more effective at maintaining genetic variation than where the environment varies stochastically. The importance of temporal environmental variation and the maintenance of genetic variation is discussed.     

GEOGRAPHIC VARIATION, FREQUENCY-DEPENDENT SELECTION, AND THE MAINTENANCE OF A FEMALE-LIMITED POLYMORPHISM

A central problem in evolutionary biology is to understand how spatial and temporal variation in selection maintain genetic variation within and among populations. Brown anole lizards ( Anolis sagrei ) exhibit a dorsal pattern polymorphism that is expressed only in females, which occur in "diamond,""bar," and intermediate "diamond-bar" morphs. To understand the inheritance of this polymorphism, we conducted a captive breeding study that refuted several single-locus models and supported a two-locus mode of inheritance. To describe geographic variation in morph frequencies, we surveyed 13 populations from two major islands in The Bahamas. Morph frequencies differed substantially between major islands but were highly congruent within each island. Finally, we measured viability selection on each island to test two hypotheses regarding the maintenance of the polymorphism: (1) that spatial variation in selection maintains variation in morph frequencies between islands, and (2) that temporal variation in selection across years maintains variation within islands. Although bar females had relatively lower survival where they were rare, our data do not otherwise suggest that selection varies spatially between islands. However, diamond-bar females were subject to positive frequency-dependent selection across years, and the relative fitness of bar and diamond females alternated across years. We propose that this polymorphism is maintained by temporal variation in selection coupled with the sheltering of alleles via a two-locus inheritance pattern and sex-limited expression.     

Stochastic selection in large and small populations

We describe two models of stochastic variation in selection intensity. In both models the arithmetic mean fitness of all genotypes is equal; in both models the geometric mean fitness of the heterozygous genotype is greater than that of both homozygous genotypes. In one model the correlation between the fitnesses of the homozygous genotypes is +1; in the other it is −1. We show that the expected time to absorption of an allele in a finite population is significantly retarded for all initial gene frequencies in the former model. The expected time to absorption of an allele in the latter model is retarded only at extreme initial gene frequencies; at intermediate initial gene frequencies the expected time to absorption is accelerated. We conclude that the criterion for polymorphism based on the geometric mean of the heterozygote being greater than that of both homozygotes provides only limited information about the fate of gene frequency.     

FINE-GRAINED SPATIAL AND TEMPORAL VARIATION IN SELECTION DOES NOT MAINTAIN GENETIC VARIATION IN ERIGERON ANNUUS

Because interactions among plants are spatially local, the scale of environmental heterogeneity can have large effects on evolutionary dynamics. However, very little is known about the spatial patterns of variation in fitness and the relative magnitude of spatial and temporal variation in selection. Replicates of 12 genotypes of Erigeron annuus (Asteraceae) were planted in 288 locations within a field, separated by distances of 0.1 to 30.0 m, and replicated in two years. In a given year, most spatial variation in relative fitness (genotype-environment [G × E] interactions for fitness) occurred over distances of only 50 cm. Year effects were as large or larger than the spatial variation in fitness; in particular there was a large, three-way, genotype-year-environment interaction at the smallest spatial scale. The genetic correlation of fitness across years at a given location was near zero, 0.03. Thus, the relative fitness of genotypes is spatially unpredictable and a map of the selective environment has constantly shifting locations of peaks and valleys. Including measurements of soil nutrients as covariates in the analysis removed most of the spatial G × E interaction. Vegetation and microtopography had no effect on the G × E terms, suggesting that differential response to soil nutrients is the cause of spatial variation in fitness. However, the slope of response to NH₄ and P0₄ was negative; therefore the soil nutrients are probably just indicators of other, unknown, environmental factors. We explored via simulation the evolutionary consequences of spatial and temporal variation in fitness and showed that, for this system, the spatial scale of variation was too fine grained (by a factor of 3 to 5) to be a powerful force maintaining genetic variation in the population. The inclusion of both spatial and temporal variation in fitness actually reduced the coexistence of genotypes compared to pure spatial models. Thus the presence of spatial or temporal variation in selection does not guarantee that it is an effective evolutionary force maintaining diversity. Instead the pattern of selection favors generalist genotypes.     

Molecular genetic markers provide no evidence for reproductive isolation among retreat building phenotypes of the net-spinning caddisfly Macrostemum carolina

Larvae of the stream-dwelling, filter-feeding caddisfly Macrostemum carolina construct silken catchnets within protective retreats. In the Savannah River, M. carolina individuals make three different retreats, each with a distinct water entrance hole: (i) at the end of a silken tube; (ii) with a approximately 180 degrees silken backstop; and (iii) flush with the top of the retreat. To resolve whether these different retreats represent alternative behavioural phenotypes within a single panmictic population or fixed phenotypes within three genetically distinct populations or species, we compared the allele frequencies at three polymorphic nuclear loci (allozyme electrophoresis for Gpi, Mpi and Pgm) and the mitochondrial DNA (mtDNA) haplotype frequencies among individuals displaying the three retreat morphs. We also calculated pairwise exact tests of population differentiation using the allozyme and mtDNA allele frequencies. No significant genetic differentiation was detected among caddisflies exhibiting the different retreat morphs. Therefore, these morphs apparently represent a single panmictic population in the Savannah River. Consequently, additional study is required to assess whether this retreat polymorphism is a phenotypically plastic trait under conditional control, or is mediated by alternative alleles at a Mendelian gene or genes (or a combination of the two).     

Ecological genetics of adaptive color polymorphism in pocket mice: geographic variation in selected and neutral genes

Patterns of geographic variation in phenotype or genotype may provide evidence for natural selection. Here, we compare phenotypic variation in color, allele frequencies of a pigmentation gene (the melanocortin-1 receptor, Mc1r), and patterns of neutral mitochondrial DNA (mtDNA) variation in rock pocket mice (Chaetodipus intermedius) across a habitat gradient in southern Arizona. Pocket mice inhabiting volcanic lava have dark coats with unbanded, uniformly melanic hairs, whereas mice from nearby light-colored granitic rocks have light coats with banded hairs. This color polymorphism is a presumed adaptation to avoid predation. Previous work has demonstrated that two Mc1r alleles, D and d, differ by four amino acids, and are responsible for the color polymorphism: DD and Dd genotypes are melanic whereas dd genotypes are light colored. To determine the frequency of the two Mc1r allelic classes across the dark-colored lava and neighboring light-colored granite, we sequenced the Mc1r gene in 175 individuals from a 35-km transect in the Pinacate lava region. We also sequenced two neutral mtDNA genes, COIII and ND3, in the same individuals. We found a strong correlation between Mc1r allele frequency and habitat color and no correlation between mtDNA markers and habitat color. Using estimates of migration from mtDNA haplotypes between dark- and light-colored sampling sites and Mc1r allele frequencies at each site, we estimated selection coefficients against mismatched Mc1r alleles, assuming a simple model of migration-selection balance. Habitat-dependent selection appears strong but asymmetric: selection is stronger against light mice on dark rock than against melanic mice on light rock. Together these results suggest that natural selection acts to match pocket mouse coat color to substrate color, despite high levels of gene flow between light and melanic populations.