Microsatellite analysis of the phylogeography, Pleistocene history and secondary contact hypotheses for the killifish, Fundulus heteroclitus

The mummichog, Fundulus heteroclitus, exhibits extensive latitudinal clinal variation in a number of physiological and biochemical traits, coupled with phylogeographical patterns at mitochondrial and nuclear DNA loci that suggest a complicated history of spatially variable selection and secondary intergradation. This species continues to serve as a model for understanding local and regional adaptation to variable environments. Resolving the influences of historical processes on the distribution of genetic variation within and among extant populations of F. heteroclitus is crucial to a better understanding of how populations evolve in the context of contemporary environments. In this study, we analysed geographical patterns of genetic variation at eight microsatellite loci among 15 populations of F. heteroclitus distributed throughout the North American range of the species from Nova Scotia to Georgia. Genetic variation in Northern populations was lower than in Southern populations and was strongly correlated with latitude throughout the species range. The most common Northern alleles at all eight loci exhibited concordant latitudinal clinal patterns, and the existence of an abrupt transition zone in allele frequencies between Northern and Southern populations was similar to that observed for mitochondrial DNA and allozyme loci. A significant pattern of isolation by distance was observed both within and between northern and southern regions. This pattern was unexpected, particularly for northern populations, given the recent colonization history of post-Pleistocene habitats, and was inconsistent with either a recent northward population expansion or a geographically restricted northern Pleistocene refugium. The data provided no evidence for recent population bottlenecks, and estimates of historical effective population sizes suggest that post-Pleistocene populations have been large throughout the species distribution. These results suggest that F. heteroclitus was broadly distributed throughout most of its current range during the last glacial event and that the abrupt transition in allele frequencies that separate Northern and Southern populations may reflect regional disequilibrium conditions associated with the post-Pleistocene colonization history of habitats in that region.     

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.     

Biochemical genetics of Fundulus heteroclitus (L.). IV. Spatial variation in gene frequencies of Idh-A,Idh-B, 6-Pgdh-A,and Est-S

The spatial variation in gene frequencies of four unlinked polymorphic loci was studied in the teleost Fundulus heteroclitus. Three loci (Idh-A, Idh-B, and Est-S) exhibit significant north-south clinal variation in allelic frequencies along the Atlantic Coast of North America, while a fourth locus (6-Pgdh-A) shows a modest clinal variation. These data, together with our previous data for Ldh-B, Mdh-A, Gpi-B, and Pgm-A, reveal a pattern of low gene diversity in the colder northern extremes of the species range and high gene diversity in warmer southern latitudes.This work was supported by Grants DEB76-19877 and DEB79-12216 from the National Science Foundation and by Grant P60-80-04 from the State of Maryland. REC and RVB were supported by NIH Training Grant GM07231 to the Department of Biology.Contribution No. 1104 from the Department of Biology, The Johns Hopkins University.     

Adaptive significance of amylase polymorphism in Drosophila I. The geographical pattern of allozyme polymorphism at the amylase locus in Drosophila subobscura

Allelic variation at the Amy locus was studied in eight natural populations from the central and northern range of D. subobscura, and the geographical pattern of Amy polymorphism over the range of this species was described. Even though regional and local differences in gene frequencies were found, in general the same alleles occur at high, intermediate and low frequencies, in nearly all populations. There are no significant differences in allele frequencies, but there is significant difference in the degree of heterozygosity among groups of populations from the northern, central and southern range. An analysis of population subdivision indicates that heterogeneity within populations is higher than between populations. Genetic distance values indicate that there is a variable degree of geographical differentiation between local populations. Variability within and between continental and insular populations is also discussed.     

Allozyme Variation and the Genetic Structure of Populations of the Rare Sedge Carex misera (Cyperaceae)

Abstract Carex misera is a rare sedge, endemic to rocky outcrops and mountain summits within the southern Appalachian Mountains from northern Georgia to northern North Carolina. We assessed allozyme diversity for 406 individuals from nine populations over most of the geographic range. Twenty-seven putative loci were assayed and eight (30%) were found to be polymorphic. Nei's gene diversity statistics ( H _T =0.043, H _S =0.019, G _ST =0.551) indicated low levels of variation but relatively highly differentiated populations, suggesting little gene flow. Significant deviations from genotypic expectations under Hardy-Weinberg equilibrium, high positive fixation indices, and the existence of small genetic neighborhoods within populations suggest that at least some inbreeding occurs. Cluster analysis of Nei's genetic identity statistics and principal component analysis of allele frequency data showed high similarity among the six southern populations with the two northern populations more differentiated from them and from each other. These results suggest that preservation of the northern populations is necessary to conserve the already low levels of genetic diversity within the species.     

Notes on the esterase variation in Swedish Dactylorhiza incarnata s. l. (Orchidaceae)

Swedish material of Dactylorhiza incarnara s. l. shows little variation at commonly investigated allozyme loci. However, interpretable variation was found at one esterase locus. All plants investigated of D. incarnata var. cruenta from southern Sweden (with spotted leaves) were homozygous for allele a at this locus, whereas all plants investigated of D. incarnata var. ochroleuca were fixed for allele b. In D. incarnata var. incarnata , both alleles were found, although the b allele dominated. In northern Swedish material of D. incarnata s. l., only allele b was found, regardless whether the material had spotted leaves (sometimes referred to as var. cruenta ), or had unspotted leaves.
These results indicate that there is restricted gene flow between var. cruenta and the other varieties in southern Sweden, although they often grow in mixed populations. The northern Swedish material with spotted leaves appears not to be related to the southern var. cruenta .     

Intraspecific variation in gene expression after seawater transfer in gills of the euryhaline killifish Fundulus heteroclitus

Previous research has suggested that northern populations of the euryhaline killifish (Fundulus heteroclitus) are better adapted to freshwater environments than their southern counterparts. In this study, we examined whether this adaptation has come at an ionoregulatory cost in seawater, by comparing published data for northern killifish to newly acquired data on the molecular responses of southern killifish to seawater transfer. After abrupt transfer from brackish water (10 per thousand) to seawater, Na,K-ATPase activity, Na,K-ATPase alpha(1a) mRNA expression, and NKCC1 mRNA expression increased 1 and 4 days after transfer in the gills of southern fish (by 2-3-fold), but increased at 1 day and not 4 days after transfer in northern fish. Small increases in mRNA expression were observed in both populations at 14 days. CFTR expression also increased in southern and northern fish at 1 and 4 days into seawater, and was also elevated at 14 days in northern fish. Because fish from both southern and northern populations maintained plasma Na(+) and Cl(-) balance after seawater transfer, the differences in activity and expression could not have been caused by differences in plasma ion levels. Instead, some other regulatory factor may account for the differences in expression between populations. This study shows that freshwater adaptation in northern populations of killifish has not necessarily come at a significant ionoregulatory cost in seawater, but has altered the molecular responses of their gills to seawater transfer compared to southern killifish.     

Evolutionary History of the Sex-Peptide (Acp70a) Gene Region in Drosophila Melanogaster

In Drosophila the products of the seminal fluid stimulate oviposition and suppress remating in the female. Of all the accessory gland peptides (Acp's) involved in these two responses, the sex-peptide (coded by the Acp70A gene) is among the best characterized at the functional level. A 1.2-kb fragment encompassing the Acp70A gene of nine lines from a natural population of D. melanogaster and one allele of D. sechellia was sequenced to study the forces shaping nucleotide variation within and between species. The coding region of D. simulans and D. mauritiana was also sequenced. A Ser to Ala replacement polymorphism at the last position of the signal peptide was detected in D. melanogaster. The Ser and Ala alleles are at intermediate frequencies. The level of nucleotide variation is lower for the derived Ala allele, which is compatible with a recent origin and an increase in frequency due to positive selection. Variation at the 5' flanking region is structured in two major highly differentiated haplotypes, whose distribution does not conform to neutral expectations. Selective and/or historical factors could contribute to the observed overall patterning of nucleotide variation at the Acp70A region.     

Intraspecific variation in aerobic metabolism and glycolytic enzyme expression in heart ventricles

A previous phylogenetic analysis among 15 taxa of the teleost fish Fundulus suggested that there should be thermal-adaptive differences in heart metabolism among populations. To test this hypothesis, the rate of oxygen consumption and the activities of all 11 glycolytic enzymes were measured in isolated heart ventricle from two populations of Fundulus heteroclitus. Heart ventricular metabolism is greater in a northern population versus a southern population of these fish. Analysis of the amount of glycolytic enzymes indicates that 87% of the variation in cardiac metabolism within and between populations is explained by the variation in three enzymes (pyruvate kinase, glyceraldehyde-3-phosphate dehydrogenase, and lactate dehydrogenase). These enzymes are the same three enzymes that were predicted to be important based on previously determined phylogenetic patterns of expression. Our data indicate that near-equilibrium enzymes, as well as classically defined rate-limiting enzymes, can also influence metabolism.     

Genetic profile of cosmopolitan populations: effects of hidden subdivision

Natural populations of many organisms exhibit excess of rare alleles in comparison with the predictions of the neutral mutation hypothesis. It has been shown before that either a population bottleneck or the presence of slightly deleterious mutations can explain this phenomenon. A third explanation is presented in this work, showing that hidden subdivision within a population can also lead to an excess of rare alleles in the total population when the expectations of the neutral model are based on the allele frequency profile of the entire population data. With two examples (mitochondrial DNA-morph distribution and isozyme allele frequency distributions), it is shown that most cosmopolitan human populations exhibit excess of rare as well as total allele counts, when these are compared with the expectations of the neutral mutation hypothesis. The mitochondrial data demonstrate that such excesses can be detected from genetic variation at a single locus as well, and this is not due to stochastic error of allele frequency distributions. Contrast of the present observations with the allele frequency profiles in agglomerated tribal populations from South and Central America shows that even when the neutral expectations hold for individual subpopulations, if all subpopulations are grouped into a single population, the pooled data exhibit an excess of total number of alleles that is mainly due to the excess of rare alleles. Therefore, a primary cause of the excess number of rare alleles could be the hidden subdivision, and the magnitude of the excess indicates the extent of substructuring. The two components of hidden subdivision are: 1) Number of subpopulations, and 2) the average genetic distance among them. The implications of this observation in estimating mutation rate are discussed indicating the difficulties of comparing mutation rates from different population surveys.     

An Indel Polymorphism in the MtnA 3' Untranslated Region Is Associated with Gene Expression Variation and Local Adaptation in Drosophila melanogaster

Insertions and deletions (indels) are a major source of genetic variation within species and may result in functional changes to coding or regulatory sequences. In this study we report that an indel polymorphism in the 3’ untranslated region (UTR) of the metallothionein gene MtnA is associated with gene expression variation in natural populations of Drosophila melanogaster. A derived allele of MtnA with a 49-bp deletion in the 3'' UTR segregates at high frequency in populations outside of sub-Saharan Africa. The frequency of the deletion increases with latitude across multiple continents and approaches 100% in northern Europe. Flies with the deletion have more than 4-fold higher MtnA expression than flies with the ancestral sequence. Using reporter gene constructs in transgenic flies, we show that the 3'' UTR deletion significantly contributes to the observed expression difference. Population genetic analyses uncovered signatures of a selective sweep in the MtnA region within populations from northern Europe. We also find that the 3’ UTR deletion is associated with increased oxidative stress tolerance. These results suggest that the 3'' UTR deletion has been a target of selection for its ability to confer increased levels of MtnA expression in northern European populations, likely due to a local adaptive advantage of increased oxidative stress tolerance.     

Postglacial migration and the loss of allozyme variation in northern populations of Asclepias exaltata (Asclepiadaceae)

The recent Wisconsin glaciation has provided opportunities for examining the effects of postglacial recolonization on the population genetics of plant and animal communities. In this study allozyme Variation was examined in 19 populations of the herbaceous perennial Asclepias exaltata occurring in previously glaciated regions of North America. These northern populations of A. exaltata possess significantly fewer polymorphic loci (46.31 ± 2.7; mean ± 1 SD), alleles per polymorphic locus (1.84 ± 0.24), and expected heterozygosity (0.133 ± 0.031) than populations found in the Pleistocene refugium in the southern Appalachians. Population-level allozyme diversity decreased linearly from south to north and from east to west. Nineteen uncommon alleles previously observed in southern Appalachian populations were undetected in the northern region. Seven common alleles exhibited a clinal change in allele frequency. Of these, only Pgd-1a and Mnr-1c were at low-frequency in the southern Appalachians and increased significantly with increasing latitude and longitude, respectively. Despite this loss of allozyme diversity following postglacial migration, northern populations of A. exaltata have higher allozyme diversity and lower population differentiation (G" = 0.1 17) than mean values for other long-lived herbaceous perennials. Increased habitat fragmentation in northern regions and potential habitat loss in the southern Appalachians are likely to reduce the historically rich gene pool that has provided the genetic stock for postglacial recoveries.     

Major histocompatibility complex variation in insular populations of the Egyptian vulture: inferences about the roles of genetic drift and selection

Insular populations have attracted the attention of evolutionary biologists because of their morphological and ecological peculiarities with respect to their mainland counterparts. Founder effects and genetic drift are known to distribute neutral genetic variability in these demes. However, elucidating whether these evolutionary forces have also shaped adaptive variation is crucial to evaluate the real impact of reduced genetic variation in small populations. Genes of the major histocompatibility complex (MHC) are classical examples of evolutionarily relevant loci because of their well-known role in pathogen confrontation and clearance. In this study, we aim to disentangle the partial roles of genetic drift and natural selection in the spatial distribution of MHC variation in insular populations. To this end, we integrate the study of neutral (22 microsatellites and one mtDNA locus) and MHC class II variation in one mainland (Iberia) and two insular populations (Fuerteventura and Menorca) of the endangered Egyptian vulture (Neophron percnopterus). Overall, the distribution of the frequencies of individual MHC alleles (n=17 alleles from two class II B loci) does not significantly depart from neutral expectations, which indicates a prominent role for genetic drift over selection. However, our results point towards an interesting co-evolution of gene duplicates that maintains different pairs of divergent alleles in strong linkage disequilibrium on islands. We hypothesize that the co-evolution of genes may counteract the loss of genetic diversity in insular demes, maximize antigen recognition capabilities when gene diversity is reduced, and promote the co-segregation of the most efficient allele combinations to cope with local pathogen communities.     

BIOCHEMICAL GENETICS OF FUNDULUS HETEROCLITUS (L.). VI. GEOGRAPHICAL VARIATION IN THE GENE FREQUENCIES OF 15 LOCI

Geographic variation in the gene frequencies corresponding to 15 polymorphic enzymes were studied in the common killifish Fundulus heteroclitus. Aat-A, Est-B, Fum-A, H6pdh-A, Mpi-A and Pgm-B showed clinal variation in allelic frequencies along the Atlantic coast of North America, while Aat-B, Ap-A, and the EST-C phenotypes did not. The clinal allelic variation of six previously examined loci (Ldh-B, Mdh-A, Gpi-B, Idh-A, Pgm-A, and 6-Pgdh-A) was extended to locations farther north. Gene diversity was lowest in the cold waters of northern latitudes and highest in the warmer waters of southern latitudes. The variety of clinal shapes and widths suggests that selection has affected the allelic distributions for at least some of these loci. This hypothesis is discussed with regard to the range contractions and extensions caused by the glacial advances and retreats during the Pleistocene.     

Selection of the sex‐linked inhibitor of apoptosis in mountain pine beetle (Dendroctonus ponderosae) driven by enhanced expression during early overwintering

The mountain pine beetle (Dendroctonus ponderosae) is an insect native to western North America; however, its geographical range has recently expanded north in BC and east into Alberta. To understand the population structure in the areas of expansion, 16 gene‐linked microsatellites were screened and compared to neutral microsatellites using outlier analyses of F_st and F_ct values. One sex‐linked gene, inhibitor of apoptosis (IAP), showed a strong signature of positive selection for neo‐X alleles and was analyzed for evidence of adaptive variation. Alleles of IAP were sequenced, and differences between the neo‐X and neo‐Y alleles were consistent with neutral evolution suggesting that the neo‐Y allele may not be under functional constraints. Neo‐Y alleles were amplified from gDNA, but not effectively from cDNA, suggesting that there was little IAP expression from neo‐Y alleles. There were no differences in overall IAP expression between males and females with the common northern neo‐X allele suggesting that the neo‐X allele in males compensates for the reduced expression of neo‐Y alleles. However, males lacking the most common northern neo‐X allele thought to be selected for in northern populations had reduced overall IAP expression in early October—at a time when beetles are preparing for overwintering. This suggests that the most common allele may have more rapid upregulation. The reduced function of neo‐Y alleles of IAP suggested by both sequence differences and lower levels of expression may foster a highly selective environment for neo‐X alleles such as the common northern allele with more efficient upregulation.     

Phosphoglucomutase polymorphism in the oriental fruit fly, Dacus dorsalis (Diptera, Tephritidae) from Peninsular Malaysia

Seven natural populations of Dacus dorsalis were analysed for phosphoglucomutase by means of horizontal starch-gel electrophoresis. The electrophoretic phenotypes were governed by four codominant Pgm alleles. The commonest allele in all the seven population samples was PgmB which encoded an electrophoretic band with intermediate mobility. The distributions of PGM phenotype were in accordance with Hardy-Weinberg expectations. There was geographic variation in the distribution of Pgm alleles.