Phylogeographic analyses of a circumarctic coastal and a boreal lacustrine mysid crustacean, and evidence of fast postglacial mtDNA rates

Phylogeographic structures of two weakly dispersing Mysis sibling species, one with a circumarctic coastal, the other with a boreal lacustrine-Baltic distribution, were studied from mitochondrial cytochrome c oxidase subunit I (COI) gene sequences. Mysis segerstralei showed high overall diversity and little phylogeographic structure across the Arctic, indicating late-glacial dispersal among coastal and lake populations from Alaska, Siberia and the north of Europe. A strongly divergent refugial lineage was however identified in Beringia. The boreal 'glacial relict'Mysis salemaai in turn displayed clear structuring among postglacially isolated Scandinavian lake populations. The inferred pattern of intralake mitochondrial DNA (mtDNA) monophyly in Scandinavia suggested relatively small population sizes and a remarkably fast postglacial mtDNA divergence rate (0.27% per 10 000 years). Nevertheless, the broader phylogeographic pattern did not support distinct eastern and western glacial refugia in Northern Europe, unlike in some other aquatic taxa. In all, the two species comprised three equidistant mitochondrial lineages (approximately 2% divergence), corresponding to M. salemaai, to the bulk of M. segerstralei, and to the Beringian M. segerstralei lineage. The lack of reciprocal monophyly of the two species in respect to their mitochondrial genealogy could indicate postspeciation mitochondrial introgression, also exemplified by an evidently more recent capture of M. segerstralei mitochondria in a Karelian population of M. salemaai. Overall, the data suggest that the continental boreal M. salemaai has a relatively recent ancestry in arctic coastal waters, whereas two other boreal 'glacial relict'Mysis sibling species in Europe (Mysis relicta) and North America (Mysis diluviana) have colonized inland waters much earlier (approximately 8% COI divergence).     

Nearshore fish (Pholis gunnellus) persists across the North Atlantic through multiple glacial episodes

The intertidal biota of the North Atlantic is characterized by two disjunct communities (North American and European) exposed to different climatic regimes during the Pleistocene and in the Holocene. We collect multilocus DNA sequence data from the nearshore fish Pholis gunnellus to help uncover processes determining biogeographical persistence during periodic coastal glaciations. Coalescent-based estimates from the multilocus DNA sequence data suggest that P. gunnellus persisted on both sides of the North Atlantic throughout the last two glacial maxima (> 202,000 years) with little trans-Atlantic gene flow since divergence, very little structure among populations within Europe (Phi(ST) < 0.05) and some structure within the North American coastline (Phi(ST) = 0.0-0.21). Although the ecological flexibility and high local migration of P. gunnellus could have enhanced this species' survival across the Atlantic, logistic regression did not find a significant determinant of trans-Atlantic persistence when considering 12 other North Atlantic phylogeographical studies from the literature.     

Seasonal variation of mass, condition and gonadosomatic index in fourhorn sculpin (Myoxocephalus quadricornis) with normal and malformed vertebrae in the Gulf of Bothnia

It may be expected that vertehral anomaties would reduer the viability and reproductive success of fish. In the present study. the relationship of vertebral defects to body mass. Iength, age, condition, gonadosomatic index and growth in fourhorn sculpin, Myoxocephalus quadricornis , in the Gulf of Bothnia was studied. There were no differences in growth between normal and defective lish, for either males or females. While normal and defectives males were similar in age and length. defective females were older and longer than normal females. Defective males and females had a significantly greater somatic mass that normal fish. However defective females may be later than in normal fish. Hence, vertebral defects do not necessarly reduce the reproductive potential of fourhorn sculpin.     

Activities of antioxidant enzymes and cytochrome c oxidase in liver of Arctic and temperate teleosts

Enhanced antioxidant status in polar fishes may occur due to high dissolved oxygen levels and membranes rich in peroxidation-sensitive polyunsaturated fatty acids. To evaluate the importance of antioxidant enzymes in polar fishes, activities of catalase (CAT), superoxide dismutase (SOD), and glutathione reductase (GR), as well as the aerobic enzyme cytochrome c oxidase (CCO), were measured at 6 degrees C and 1 degrees C in livers of confamilial Arctic and temperate teleosts: the Arctic fourhorn sculpin Triglopsis quadricornis (Cottidae) and saddled eelpout Lycodes mucosus (Zoarcidae) vs. the temperate longhorn sculpin (Myoxocephalus octodecimspinosus) (Cottidae) and ocean pout (Zoarces americanus) (Zoarcidae), respectively. At both assay temperatures, CAT activities were substantially lower in both Arctic species, SOD was similar in the cottids but lower in the Arctic zoarcid, and GR was similar in temperate and Arctic fishes. Activities at respective habitat temperatures were always significantly lower in the Arctic fishes. The lower antioxidant enzyme activities in the Arctic fishes cannot be attributed to lower aerobic status because CCO activity was similar or higher in the Arctic fishes; significant negative relationships were found between CCO and CAT and GR (but not SOD) when all species were combined, indicating that a higher apparent aerobic status does not necessarily coincide with higher antioxidant enzyme activities. Antioxidant enzyme activities may not be enhanced as part of cold adaptation in Arctic fishes, at least in the liver.     

Molecular phylogeography and cryptic speciation in the mosses, Mielichhoferia elongata and M. mielichhoferiana (Bryaceae)

Nucleotide sequence variation in the ITS1-5.8S-ITS2 region of nuclear ribosomal DNA (nrDNA) from 70 populations of Mielichhoferia elongata and M. mielichhoferiana, plus two outgroup species, was analysed using maximum parsimony and maximum likelihood methods. High levels of nucleotide substitution and numerous insertion-deletion events were detected within and between the two species. M. elongata is monophyletic with regard to nrDNA variation, but M. mielichhoferiana is paraphyletic. (M. elongata is nested within it.) A clade within M. mielichhoferiana provides evidence of vicariance, with North American and Scandinavian sister groups of populations. Two major clades are resolved in M. elongata by sequence data that are completely congruent with previous isozyme work. One clade includes populations from both North America and Europe whereas the other is strictly North American. These two clades, resolved by multiple independent loci, clearly represent cryptic species within the morphologically uniform M. elongata. Certain geographical areas, most notably southwestern Colorado in Ouray and San Juan Counties, harbour diverse populations of M. elongata with distinct phylogenetic and phylogeographical histories. Morphologically indistinguishable but phylogenetically distant populations were detected a few metres apart at one site. In contrast, all populations collected over hundreds of kilometres in California belong to a single clade. Arctic North American populations belong to a clade that includes disjunct populations in Alaska, northern Ellesmere Island, and the northeastern USA, but not subarctic Swedish populations, which are more closely related to plants from the Rocky Mountains. Morphological uniformity belies complex infraspecific phylogenetic patterns within M. elongata and M. mielichhoferiana.     

Phylogeography and postglacial dispersal of arctic charr Salvelinus alpinus in North America

Restriction fragment length polymorphism analysis of mitochondrial DNA (mtDNA) was used to reconstruct postglacial dispersal routes of arctic charr Salvelinus alpinus in North America. Twelve of 35 restriction enzymes detected polymorphisms among representative populations, revealing two distinct lineages with an estimated nucleotide divergence of 1.32%. Subsequent screening of 869 fish from 54 populations with four diagnostic restriction enzymes showed that these lineages have largely allopatric distributions, suggesting their dispersal from separate northern and eastern glacial refugia. In addition, geographical and genetic structure among eastern populations suggested the existence of a second eastern refuge. Among the three lineages, the most divergent (Arctic) lineage occurred from Alaska east to northern Labrador. Quebec, New Brunswick, and New England were colonized by a second (Laurentian) lineage, and Labrador by a third group. Contact between refugial groups was only detected in two Labrador populations. The Arctic lineage was highly differentiated from eastern North American and European haplotypes, and probably diverged during the early Pleistocene. By contrast, the Laurentian and Labrador groups were similar to Old World charr, suggesting a shared ancestry during the mid-Pleistocene. In addition, the close relationship between Labrador and Laurentian charr indicates their probable divergence during the Wisconsinan glaciation.     

Comparative phylogeography of two North American 'glacial relict' crustaceans

The Pleistocene glaciations represent the most recent and dramatic series of habitat changes since the Cretaceous. The impact of these events was particularly acute for aquatic taxa with poor powers of dispersal, but few organisms have evolutionary histories more intimately entwined with the advance and retreat of ice than the 'glacial relicts'. In this study, we used a mitochondrial gene, cytochrome c oxidase subunit I (COI), to examine and compare the phylogeographical structure of two glacial relict crustaceans (Limnocalanus macrurus and members of the Mysis relicta species group) across North America. In both cases, we found a sharp phylogenetic division between populations from inland lakes formed during glacial retreat, and arctic lakes isolated from polar seas via isostatic rebound. However, the depth of this phylogenetic partition varied between taxa. In L. macrurus, nucleotide sequence divergence of 2.2% between these zones is consistent with its current status as a single morphologically variable species, but in Mysis the split occurred among recently described, morphologically conserved species, at a divergence of 8.2%. The disparity in the depth of divergence indicates a history of recurrent freshwater invasions from the arctic seas, in concordance with previous studies of Eurasian glacial relicts. However, we suggest further consideration of a largely overlooked explanation that could account for some of the discrepancies between molecular divergences and glaciation events. Many cladogenetic events could have occurred in arctic seas prior to the transition to inland waters, a possibility supported both by the complex physical and ionic history of arctic seas and by high marine and estuarine lineage diversity in the north.     

Opsin gene sequence variation across phylogenetic and population histories in Mysis (Crustacea: Mysida) does not match current light environments or visual-pigment absorbance spectra

The hypothesis that selection on the opsin gene is efficient in tuning vision to the ambient light environment of an organism was assessed in 49 populations of 12 Mysis crustacean species, inhabiting arctic marine waters, coastal littoral habitats, freshwater lakes ('glacial relicts') and the deep Caspian Sea. Extensive sequence variation was found within and among taxa, but its patterns did not match expectations based on light environments, spectral sensitivity of the visual pigment measured by microspectrophotometry or the history of species and populations. The main split in the opsin gene tree was between lineages I and II, differing in six amino acids. Lineage I was present in marine and Caspian Sea species and in the North American freshwater Mysis diluviana, whereas lineage II was found in the European and circumarctic fresh- and brackish-water Mysis relicta, Mysis salemaai and Mysis segerstralei. Both lineages were present in some populations of M. salemaai and M. segerstralei. Absorbance spectra of the visual pigment in nine populations of the latter three species showed a dichotomy between lake (λ(max) =554-562 nm) and brackish-water (Baltic Sea) populations (λ(max) = 521-535 nm). Judged by the shape of spectra, this difference was not because of different chromophores (A2 vs. A1), but neither did it coincide with the split in the opsin tree (lineages I/II), species identity or current light environments. In all, adaptive evolution of the opsin gene in Mysis could not be demonstrated, but its sequence variation did not conform to a neutral expectation either, suggesting evolutionary constraints and/or unidentified mechanisms of spectral tuning.     

Microsatellite DNA analysis of northern pike (Esox lucius L.) populations: insights into the genetic structure and demographic history of a genetically depauperate species

The northern pike Esox lucius L. is a freshwater fish exhibiting pronounced population subdivision and low genetic variability. However, there is limited knowledge on phylogeographical patterns within the species, and it is not known whether the low genetic variability reflects primarily current low effective population sizes or historical bottlenecks. We analysed six microsatellite loci in ten populations from Europe and North America. Genetic variation was low, with the average number of alleles within populations ranging from 2.3 to 4.0 per locus. Genetic differentiation among populations was high (overall θ_ST = 0.51; overall ρ_ST = 0.50). Multidimensional scaling analysis of genetic distances between populations and spatial analysis of molecular variance suggested a single phylogeographical race within the sampled populations from northern Europe, whereas North American and southern European populations were highly distinct. A population from Ireland was monomorphic at all loci, presumably reflecting founder events associated with introduction of the species to the island in the sixteenth century. Bayesian analysis of demographic parameters showed differences in θ (a product of effective population size and mutation rate) among populations from large and small water bodies, but the relative differences in θ were smaller than expected, which could reflect population subdivision within the larger water bodies. Finally, the analyses showed drastic population declines on a time scale of several thousand years within European populations, which we ascribe to either glacial bottlenecks or postglacial founder events.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 91–101.     

Around or across the Carpathians: colonization model of the Danube basin inferred from genetic diversification of stone loach (Barbatula barbatula) populations

Despite increasing information about postglacial recolonization of European freshwater systems, very little is known about pre-Pleistocene history. We used data on the recent distribution and phylogenetic relationships of stone loach mitochondrial lineages to reconstruct the initial colonization pattern of the Danube river system, one of the most important refuges for European freshwater ichthyofauna. Fine-scale phylogeography of the Danubian populations revealed five highly divergent lineages of pre-Pleistocene age and suggested the multiple origin of the Danubian stone loach. The mean sequence divergence among lineages extended from 7.0% to 13.4%, which is the highest intraspecific divergence observed so far within this river system. Based on the phylogeographical patterns, we propose the following hypothesis to relate the evolution and dispersal of the studied species with the evolution of the Danube river system and the Carpathian Mountains: (i) during the warmer period in the Miocene, the areas surrounding the uplifting Alps and Carpathians served as mountainous refuges for cold-water adapted fish and promoted the diversification of its populations, and (ii) from these refuges, colonization of the emerging Danube river system may have taken place following the retreat of the Central Paratethys. Co-existence of highly divergent mtDNA lineages in a single river system shows that range shifts in response to climatic changes during the Quaternary did not cause extensive genetic homogenization in the stone loach populations. However, the wide distribution of some mtDNA lineages indicates that the Pleistocene glaciations promoted the dispersal and mixing of populations through the lowlands.     

Extreme intraspecific mitochondrial DNA sequence divergence in Galaxias maculatus (Osteichthys: Galaxiidae), one of the world's most widespread freshwater fish

Biogeographic controversies surrounding the widespread freshwater fish, Galaxias maculatus, were addressed with DNA sequence data. Mitochondrial cytochrome b and 16S rRNA sequences were obtained from representatives of six populations of this species. Substantial levels of cytochrome b (maximum 14.6%) and 16S rRNA sequence divergence (maximum 6.0%) were detected between western Pacific (Tasmania-New Zealand) and South American (Chile-Falkland Islands) haplotypes. A considerable level of divergence was also detected between Tasmanian and New Zealand haplotypes (maximum 5.1%) and within and among Chilean and Falkland Island G. maculatus (maximum 3. 8%). The phylogenetic structure of haplotypes conflicts with the accepted pattern of continental fragmentation. Molecular clock calibrations suggest that haplotype divergences postdate the fragmentation of Gondwana. These findings point to marine dispersal rather than ancient vicariance as an explanation for the wide distribution. The phylogenetic structure of South American haplotypes was not consistent with their geographic distribution. We consider factors such as population divergence, population size, dispersal, secondary contact, and philopatry as potential causes of the high level of mtDNA nucleotide diversity in this species.     

Corynosoma acanthocephalans in their paratenic fish hosts in the northern Baltic Sea

Cystacanth stages of three Corynosoma (Acanthocephala) species, C. strumosum and C. semerme, and a new species in the Bothnian Bay, C. magdaleni, were studied in the fourhorn sculpin (Myoxocephalus quadricornis) in 1996-1997. The length of the trunk and proboscis differentiated the three species in the fish. The stability of the Corynosoma infection was studied by comparing the present results with those collected from the same areas (central and coastal) in 1977-1982 (Valtonen, 1983a). As C. magdaleni and C. strumosum were not separated at that time, the joint infections of these two species (called "C. strumosum") were compared. The stability of infection with C. semerme in the central Bothnian Bay was noticeable (prevalences were 82.9 in earlier period and 81.9% in later period), while the prevalences of "C. strumosum" in the same areas had decreased (21 and 13.5%, respectively). This is suggested to be due to the disappearance of the marine bull-rout, Myoxocephalus scorpius, from the Bothnian Bay during the 1990s due to an overall decline in salinity throughout the Baltic Sea. In the coastal area Corynosoma infection was clearly lower than in central area in both periods.     

Colonization, extinction, and phylogeographic patterning in a freshwater crustacean

Phylogeographic analyses have revealed the importance of Pleistocene vicariance events in shaping the distribution of genetic diversity in freshwater fishes. However, few studies have examined the patterning of variation in freshwater organisms with differing dispersal syndromes and life histories. The present investigation addresses this gap, examining the phylogeography of Sida crystallina, a species whose production of diapausing eggs capable of passive dispersal was thought to constrain its regional genetic differentiation. By contrast, the present analysis has revealed deep allozyme and cytochrome oxidase I mitochondrial DNA divergence between populations from North America and Europe. Moreover, North American populations are separated into four allopatric assemblages, whose distribution suggests their derivation from different Pleistocene refugia. These lineages show higher haplotype diversity and deeper sequence divergence than those of any fish from temperate North America. Its distinctive life history traits have evidently sheltered lineages of Sida from extinction, contributing to a remarkably comprehensive and high resolution phylogeographic record.     

American marten (Martes americana) in the Pacific Northwest: population differentiation across a landscape fragmented in time and space

American marten (Martes americana) have a close association with mature temperate forests, a habitat that expanded throughout the Pacific Northwest as glaciers receded at the end of the Pleistocene. Similar to several other forest-associated mammals in North America (e.g. black bear), genetic analysis of the marten shows a deep phylogeographical subdivision that reflects populations with distinctive evolutionary histories. Using a suite of 14 microsatellite markers, we explored the genetic structure of marten populations in two reciprocally monophyletic clades in the Pacific Northwest identified previously as M. caurina and M. americana by mitochondrial haplotypes and morphology. Microsatellite phylogeographical patterns were congruent with mitochondrial analyses. These independent data sets shed light upon hybridization patterns, population structure and evolutionary histories. Hybridization between M. caurina and M. americana individuals was documented in two regions of sympatry (Kuiu Island in southeastern Alaska and southern Montana). Northern insular populations of M. caurina exhibited higher differentiation and lower variability relative to northern populations of M. americana. Greater divergence among M. caurina populations may reflect longer isolation and persistence in coastal forest habitat that was fragmented by rising sea level in the early Holocene. Lower differentiation among northern M. americana populations and close relationship to other continental M. americana populations may reflect more recent expansion into the Pacific Northwest and/or continued gene flow among populations. Differentiation among M. caurina populations was attributed to habitat fragmentation (i.e. rising sea level), as opposed to isolation-by-distance; oceanic straits pose significant barriers to gene flow among M. caurina populations and between populations of M. caurina and M. americana.     

Feeding ecology of dominant groundfish in the northern Bering Sea

We investigated current diets of the six most abundant benthic fish in the northern Bering Sea. Our objective was to explore feeding strategies and potential competition with other top predators as ecosystem changes occur in the northern Bering Sea ecosystem. Our approach used stomach content data collected from field sampling during spring 2006 and 2007. Calanoid copepods and ampeliscid amphipods were important prey of Arctic cod (Boreogadus saida) but in different proportions depending upon fish size, feeding location, and local environmental conditions. Snailfish (Liparidae) occupied a broad niche and fed on a variety of benthic amphipods. Arctic alligatorfish (Ulcina olrikii) and Arctic staghorn sculpin (Gymnocanthus tricuspis) consumed ampeliscid amphipods predominantly. Shorthorn sculpin (Myoxocephalus scorpius) had a less-diverse diet, with snow crab (Chionoecetes opilio) most important by weight. Finally, all Bering flounder (Hippoglossoides robustus) sampled had empty stomachs. Our results indicate that ampeliscid amphipods, which have high biomass in the central region of the northern Bering Sea, are the most important prey for the dominant groundfish in the Chirikov Basin. Generally, all dominant benthic fish in the northern Bering Sea had narrow feeding niches, except snailfish. High diet overlap was found among many of the fish species, including Arctic cod and snailfish, snailfish and Arctic alligatorfish, and Arctic alligatorfish and Arctic staghorn sculpin. These findings are consistent with a relatively short food chain for benthic fish that are for the most part specialized feeders with narrow preferences for food and may be affected by changes in benthic prey distributions.     

Holarctic phylogeography of the root vole (Microtus oeconomus): implications for late Quaternary biogeography of high latitudes

A species-wide phylogeographical study of the root vole (Microtus oeconomus) was performed using the whole 1140 base pair mitochondrial (mt) cytochrome b gene. We examined 83 specimens from 52 localities resulting in 65 unique haplotypes. Our results demonstrate that the root vole is divided into four main mtDNA phylogenetic lineages that seem to have largely allopatric distributions. Net divergence estimates (2.0-3.5%) between phylogroups, as well as relatively high nucleotide diversity estimates within phylogroups, indicate that the distinct phylogeographical structure was initiated by historical events that predated the latest glaciation. European root voles are divided into a Northern and a Central mtDNA phylogroup. The mtDNA data in concert with fossil records imply that root voles remained north of the classical refugial areas in southern Europe during the last glacial period. The currently fragmented populations in central Europe belong to a single mtDNA phylogroup. The Central Asian and the North European lineages are separated by the Ural Mountains, a phylogeographical split also found in collared lemmings (Dicrostonyx) and the common vole (M. arvalis). The Beringian lineage occurs from eastern Russia through Alaska to northwestern Canada. This distribution is congruent with the traditional boundaries of the Beringian refugium and with phylogeographical work on other organisms. In conclusion, similarities between the phylogeographical patterns in the root vole and other rodents, such as Arctic and subarctic lemmings, as well as more temperate vole species, indicate that late Quaternary geological and climatic events played a strong role in structuring northern biotic communities.     

Phylogeography and the origins of range disjunctions in a north temperate fish,the pygmy whitefish (Prosopium coulterii), inferred from mitochondrial and nuclear DNA sequence analysis

Aim To investigate the degree of phylogeographical divergence within pygmy whitefish (Prosopium coulterii) and to test hypotheses concerning the origin of disjunct populations within North America. Location North America from western Alaska to Lake Superior. Methods Mitochondrial (ATPase subunit VI) and nuclear (ITS‐1, ITS‐2) DNA sequence variation was assessed across the species’ North American range to test for the existence of distinct phylogeographical groupings of pygmy whitefish associated with known glacial refugia. Coalescent simulations of the mitochondrial DNA (mtDNA) data were used to test hypotheses of population structure. Results This species is composed of two monophyletic mitochondrial clades across its North American range. The two mtDNA clades differed by an average 3.3% nucleotide sequence divergence. These clades were also distinguished by ITS‐2, but the relationships among lineages were not resolved by the ITS‐1 analysis. Coalescent analyses rejected the null hypothesis of the current disjunct distributions being a result of fragmentation of a single widespread ancestral lineage across a variety of effective population sizes and divergence times. Main conclusions The current range disjunctions of pygmy whitefish in North America probably resulted from isolation, genetic divergence, and selective dispersal from at least two major Pleistocene glacial refugia: Beringia and Cascadia. More recent isolation and dispersal from an upper Mississippi refugium is suggested by relationships within one of the clades and by distributional evidence from co‐distributed species. The Beringian and Cascadian refugia have played major roles in the zoogeography of Nearctic temperate aquatics, but the roles of smaller refugia appear more variable among other species.     

Atlantic salmon colonization of the Russian Arctic coast: pioneers from North America

Previous studies of the ESTD * isozyme locus in the Atlantic salmon show the * 80 allele to be absent across the species' European range, with the exception of northern Russia, whereas the allele is nearly fixed in North American populations. The allele was found in samples from 15 out of 18 rivers on the Kola Peninsula and White sea coast and had frequencies that ranged from 0·017 to 0·363. Typing of fish in nine of these rivers for mtDNA variation in the ND1 gene region found variation characteristic of North American salmon in three Kola Peninsula populations and is the only part of Europe where such variants have been detected. The study area was completely glaciated during the late Pleistocene period and the restriction of the European distribution of these ESTD * and mtDNA variants to this area suggests that salmon that colonized the Russian Arctic coast rivers included fish of North American origin after the Pleistocene glacier had retreated.     

Genetic divergence at the NAD +-dependent malic enzyme locus in Atlantic salmon from Europe and North America

The distribution of genetic variation for NAD+ malic enzyme ( ME* ) polymorphism in the Atlantic salmon was assessed in both anadromous and resident populations. The analysis revealed a major allelic divergence between North America and Europe. The *80 variant occurred in 39 of 40 North American samples, ranging in frequency from 0 to 0.515. In contrast, it was detected in only three fish from two of the 35 European locations analysed. Eleven fish from two rivers in north-west France had the * 110 variant. While heterogeneity among North American populations was significant, no regional differentiation was apparent and anadromous and resident salmon were not found to differ.     

Population genetics of grotto sculpin (Cottus specus), a new cave-adapted fish species

The grotto sculpin (Cottus specus) is a troglomorphic fish endemic to cave systems in Perry County, Missouri. These fish are state-threatened and were recently listed as federally endangered under the Endangered Species Act (ESA). Due to the unstable nature of the cave environment, grotto sculpins are highly susceptible to pollution via agricultural and waste runoff. Here we provide a summary of population genetic diversity within the range of the grotto sculpin. Samples from twenty-one cave, surface, and resurgence sites within and surrounding the Bois Brule drainage in Perry County, southeast Missouri, USA were collected and sequence data from the mitochondrial control were analyzed using analysis of molecular variance, maximum likelihood, parsimony, and Bayesian techniques. We found a substantial degree of divergence from surface populations of Midlands and Black River banded sculpin races (7.7–8.0 %). Within the Bois Brule drainage, two distinct lineages of grotto sculpin were identified which correspond with northern and southern cave system boundaries and may represent independent cave invasions.