Genetic structure of Patagonian toothfish (Dissostichus eleginoides) populations on the Patagonian Shelf and Atlantic and western Indian Ocean Sectors of the Southern Ocean

The genetic structure of Patagonian toothfish populations in the Atlantic and western Indian Ocean Sectors of the Southern Ocean (SO) were analysed using partial sequences of the mitochondrial 12S rRNA gene and seven microsatellite loci. Both haplotype frequency data (F _ST>0.906, P<0.01) and microsatellite genotype frequency data (F _ST=0.0141–0.0338, P<0.05) indicated that populations of toothfish from around the Falkland Islands were genetically distinct from those at South Georgia (eastern Atlantic Sector SO), around Bouvet Island (western Atlantic Sector SO) and the Ob Seamount (western Indian Ocean Sector of the SO). Genetic differentiation between these populations is thought to result from hydrographic isolation, as the sites are separated by two, full-depth, ocean-fronts and topographic isolation, as samples are separated by deep water. The South Georgia, Bouvet and Ob Seamount samples were characterised by an identical haplotype. However, microsatellite genotype frequencies showed genetic differentiation between South Georgia samples and those obtained from around Bouvet Island and nearby seamounts (F _ST=0.0037, P<0.05). These areas are separated by large geographic distance and water in excess of 3,000 m deep, below the distributional range of toothfish (<2,200 m). No significant genetic differentiation was detected between samples around Bouvet Island and the Ob Seamount although comparisons may have been influenced by low sample size. These localities are linked by topographic features, including both ridges and seamounts, that may act as oceanic “stepping stones” for migration between these populations. As for other species of deep-sea fish, Patagonian toothfish populations are genetically structured at the regional and sub-regional scales.     

Marine biogeographic disjunction in central New Zealand

We present a phylogeographic analysis of an abundant New Zealand endemic sea-star, Patiriella regularis, to help pinpoint the location of an important biogeographic disjunction in central New Zealand. The analysis incorporates 284 mtDNA control region sequences (approximately 800 bp) of P. regularis from 22 coastal locations around New Zealand. We detected 132 haplotypes, with a mean divergence of 0.96%. AMOVA analysis of New Zealand samples is consistent with a north-south biogeographic disjunction across central New Zealand (among-group genetic variance=6.10%; P=0.0005). Cook Strait, the shallow marine strait separating the main islands, is not correlated with the disjunction: samples from northern South Island are genetically indistinguishable from North Island samples (variance=1.69%; P=0.073). These results are consistent with the hypothesis that upwelling zones south of Cook Strait constitute a significant barrier to larval dispersal.Communicated by M.S. Johnson, Crawley     

Population connectivity in the temperate damselfish <Emphasis Type="Italic">Parma microlepis</Emphasis>: analyses of genetic structure across multiple spatial scales

This study investigated the utility of microsatellite markers for providing information on levels of population connectivity for a low dispersing reef fish in New South Wales (NSW), Australia, at scales ≤400 km. It was hypothesized that the temperate damselfish Parma microlepis, which produces benthic eggs and has limited post-settlement dispersal, would exhibit spatial genetic structure and a significant pattern of isolation-by-distance (IBD). A fully nested hierarchical sampling design incorporating three spatial scales (sites, location and regions, separated by 1–2, 10–50 and 70–80 km respectively) was used to determine genetic variability at seven microsatellite loci. Broad-scale genetic homogeneity and lack of IBD was well supported by single and multi-locus analyses. The proportion of the total genetic variation attributable to differences among regions, locations or sites was effectively zero (Φ/R-statistics ≤0.007). The geographic distribution of genetic diversity and levels of polymorphism (H _E 0.21–0.95) indicate high mutation rates, large effective population sizes, and high rates of gene flow. Significant gene flow may be driven by factors influencing pre-settlement dispersal, including the East Australian Current (EAC) and habitat continuity. Genetic connectivity may not reflect demographically important connectivity, but does imply that P. microlepis populations are well connected from an evolutionary perspective. Total observed genetic diversity was accounted for within 1–2 km of reef and could be represented within small Marine Protected Areas. Reef fishes in NSW which have life histories similar to P. microlepis (e.g. pre-settlement durations ≥2 weeks) are also likely to exhibit genetic homogeneity. Genetic markers are, therefore, most likely to provide information on demographically relevant connectivity for species with lower dispersal capabilities, small population sizes, short life spans, and whose habitats are rare, or patchily distributed along-shore. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Population variation in the mitochondrial cytochrome b gene of the orange roughy Hoplostethus atlanticus and the hoki Macruronus novaezelandiae

To describe the genetic relationship among regional populations of two commercially valuable species of marine fish, the orange roughy Hoplostethus atlanticus and the hoki Macruronus novaezelandiae, the mitochondrial (mt) DNA haplotypes of 59 individuals were defined by direct sequencing of the cytochrome b gene. Samples of orange roughy were collected on four fishing grounds around New Zealand, one off the west coast of Tasmania and one near South Africa from 1990–1991. Samples from hoki were collected on three fishing grounds around New Zealand and one off Tasmania during the same period. An average of 252 nucleotides were sequenced from each of 32 orange roughy and an average of 372 nucleotides from each of 27 hoki. Sequence variation allowed the definition of 9 unique orange roughy haplotypes and 5 hoki haplotypes. Genetic variation, as measured by both average sequence divergence and haplotype diversity, was high in the orange roughy (nucleotide diversity=0.590%, haplotype diversity=0.782) and low in the hoki (nucleotide diversity=0.076%, haplotype diversity=0.279) relative to a similar survey of the Atlantic cod, Gadus morhua. Differences in haplotype frequencies of orange roughy from New Zealand, Tasmania and South Africa were not significant, and the most common haplotype was found in similar frequencies in these three geographically distant regions. Differences in haplotype frequencies between the New Zealand and Tasmanian samples of hoki were significant, suggesting restricted gene flow between these two regions. The contrasting patterns of low but regionally subdivided genetic variation in the hoki versus high but geographically undifferentiated genetic variation in the orange roughy may be attributed to the low fecundity, slow maturation and long lifespan of the orange roughy relative to the hoki.     

Complex genetic patterns and a phylogeographic disjunction among New Zealand mud snails Zeacumantus subcarinatus and Z. lutulentus

The comparative analysis of genetic structure and phylogeography of marine species can reveal the relative importance of common biogeographic barriers and species-specific evolutionary histories. In the present study, mitochondrial COI sequences from a total of 724 individuals collected from 38 localities throughout New Zealand during 2006 and 2007 were used to examine the genetic structure and demographic histories of the intertidal gastropods Zeacumantus subcarinatus and Zeacumantus lutulentus. For both species, results revealed isolation of populations located along the northern South Island and southern and western North Island, and extensive genetic structure throughout the remainder of their ranges. Despite this, long-distance dispersal and secondary admixture of divergent haplotypes was evident, especially for the widespread Z. subcarinatus. These findings reveal the importance of common barriers to gene flow and highlight how the interaction of inherent dispersal limitations of direct-developing marine organisms and periodic long-distance movements can produce complex genetic patterns.     

Population structure of the wreckfish Polyprion americanus determined with microsatellite genetic markers

 We examined population structure in the wreckfish, Polyprion americanus, by assaying six microsatellite loci in 422 individuals collected throughout the geographic range. Eighteen hapuku, P. oxygeneios, were assayed at the same loci for use as an outgroup. Expected heterozygosities ranged from 0.49 to 0.88 and averaged 0.66. Allele-frequency distributions at those loci indicated that samples from the eastern North Atlantic, western North Atlantic and the Mediterranean were genetically similar, confirming the pattern seen in a previous analysis of mitochondrial DNA (mtDNA) restriction fragment length polymorphisms (RFLPs). Both mtDNA and microsatellite studies differentiated northern and southern wreckfish stocks. However, in contrast to the mtDNA studies, allelic variation at microsatellite loci clearly differentiated wreckfish from two Southern Hemisphere locations, Brazil and the South Pacific. Far more genetic variation was observed at microsatellite loci than with mtDNA RFLPs (haplotype diversity averaged 0.01), and we saw more evidence of population structure with the microsatellite loci. The differentiation between southern and northern wreckfish supports the distribution records, which indicate that wreckfish do not occur in the tropics. Temperature profiles and current patterns throughout the southern oceans apparently also prevent significant gene flow between the South Pacific and Brazilian samples. Received: 29 January 2000 / Accepted: 27 June 2000     

Global population divergence of the sea star Hippasteria phrygiana corresponds to the onset of the last glacial period of the Pleistocene

Genetic structure and connectivity of populations of the globally distributed and eurybathic sea star Hippasteria phrygiana (Parelius 1768) were studied in 165 individuals sampled from three oceanic regions: the North Pacific Ocean, the South Pacific Ocean (considered to include the adjacent regions of the Southern Ocean and the southern Indian Ocean) and the North Atlantic Ocean. A nuclear gene region (ATP synthase subunit α intron #5, ATPSα) and a mitochondrial gene region (cytochrome oxidase subunit I, COI) were amplified and sequenced. Significant heterogeneity was detected in an AMOVA analysis among the three sampled oceanic regions for COI, but not for ATPSα. Neither gene showed significant genetic heterogeneity within the North Atlantic, as assessed by ΦST values. Significant heterogeneity was detected for COI (but not ATPSα) in the North Pacific, but the converse was true in the South Pacific. Coalescent simulations suggested that the three regions have been diverging with little or no gene flow for the past 50–75,000 years, a time frame that corresponds to the onset of the last glacial period of the Pleistocene. A possible genetic signature of recent population expansion (or non-neutrality) was detected for each gene in the North Pacific, but not in the other two oceanic regions.     

Genetic structure of natural populations of California red abalone (<Emphasis Type="Italic">Haliotis rufescens</Emphasis>) using multiple genetic markers

Mitochondrial cytochrome oxidase subunit one (COI) sequence, nuclear microsatellites, and amplified fragment length polymorphisms (AFLPs) were used to evaluate connectivity among nine red abalone (Haliotis rufescens) populations sampled between August 1998 and November 2003 along approximately 1,300 km of California coastline from Crescent City (41°46′N, 124°12′W) to San Miguel Island (34°02′N, 120°22′W). COI sequences and microsatellite genotypes did not show significant genetic divergence among nine sampled populations. A subset of five populations spanning the geographic range of the study was scored for 163 polymorphic AFLP markers. Of these, 41 loci showed significant divergence (P < 0.001) among populations. Still, no AFLP markers were diagnostic for any of the study populations, and assignment tests did not consistently assign individuals to the correct population. Although the AFLP data are the first to suggest there is significant genetic differentiation among California red abalone populations, the discordance between the different genetic markers needs further study before unambiguous conclusions can be drawn with respect to connectivity among the populations. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Evidence of stock separation in southern hemisphere organge roughy (Hoplostethus atlanticus,Trachichthyidae) from restriction-enzyme analysis of mitochondrial DNA

Restriction enzyme analysis of mitochondrial DNA (mtDNA) was used to test for genetic homogeneity of orange roughy (Hoplostethus atlanticus) in the southern hemisphere. Two hundred and eighty-six orange roughy specimens were collected from seven general localities: the Great Australian Bight; South Australia (off southeastern Kangaroo Island); the west coast of Tasmania; the east coast of Tasmania; New South Wales; New Zealand and South Africa. Mitochondrial DNA was extracted from developing ovary tissue and analysed with 10 six-base enzymes and 3 four-base enzymes. Both forms of analysis revealed a low level of genetic diversity in this species. The six-base enzyme study found no evidence of reproductively isolated populations of orange roughy in southeastern Australian waters. However, an analysis of 107 fish with 3 four-base enzymes identified at least partial genetic separation of the New South Wales (NSW) sample of orange roughy from South Australian (off southeastern Kangaroo Island) and Tasmanian samples. This finding supports biological evidence for the presence of a distinct subpopulation of orange roughy in NSW waters. The four-base study also provided evidence of the presence of genetically distinct samples of orange roughy occurring in the same localities off southeastern Kangaroo Island from consecutive years. Additional sampling and the use of a greater number of four-base enzymes may be needed to determine if any genetic structuring exists among orange roughy south of New South Wales.     

Microsatellite variation, effective population size, and population genetic structure of vermilion snapper, Rhomboplites aurorubens, off the southeastern USA

Vermilion snapper (Rhomboplites aurorubens) were collected from four sites off the Atlantic coast of the USA and one site in the Gulf of Mexico to evaluate effective population size and genetic stock structure. Previous studies had suggested geographic variation in the ratio of males to females, so this population characteristic was explored in conjunction with the genetic analysis. Sex ratio varied greatly among the five sample sites; males comprised 57% of samples in the Gulf of Mexico, while within the South Atlantic Bight they comprised between 36% (Morehead City, North Carolina) and 53% (Carolina Beach, North Carolina) of samples. No clear geographic trends in the sex ratio emerged; instead, it was found to vary with fish length, the percentage of males decreasing with increasing size. Allelic variation assessed at seven dinucleotide microsatellite loci was large; gene diversities ranged from 0.43 to 0.95 and allelic counts from 7 to 39. Estimates of the effective population size ranged from 24 500 (based on the infinite-alleles model) to 150 500 (based on the stepwise-mutation model). There was evidence for excess homozygosity within samples: estimates of F _IS (the correlation of alleles within individuals) ranged from 0.01 to 0.03 among the seven loci, and three estimates were significantly greater than zero. Differentiation among localities was very weak, as estimates of F _ST (the correlation of alleles within populations) were on the order of 0.001 to 0.002 and genetic distance estimates between localities were not related to geographic distances. This suggested that vermilion snapper in the South Atlantic Bight (Cape Hatteras, North Carolina to Cape Canaveral, Florida) and Gulf of Mexico are likely to consist of one genetic stock. Despite the overall homogeneity, there were indications of a temporally dynamic local structure that would bear further examination. Received: 6 July 1998 / Accepted: 9 February 1999     

Global population structure of albacore (Thunnus alalunga) inferred by RFLP analysis of the mitochondrial ATPase gene

The genetic population structure of the highly migratory albacore (Thunnus alalunga) was investigated using restriction fragment length polymorphism (RFLP) analysis of the mitochondrial ATPase gene amplified by the polymerase chain reaction (PCR). 620 individuals comprising 13 geographically distant samples (ten Pacific, two Atlantic and the Cape of Good Hope) were surveyed between 1991 and 1994 with two restriction endonucleases (Mse I and Rsa I), resulting in seven haplotypes. No heterogeneity was observed in the distribution of haplotypes among the ten samples from the North and South Pacific, nor among the samples from North and South Atlantic and Cape of Good Hope. However, highly significant heterogeneity was evident among Atlantic and Pacific samples. Higher haplotypic diversity (h) was observed in the Pacific samples (0.59 to 0.69) than in the Atlantic and Cape samples (0.22 to 0.43). These results suggest greater gene flow between albacore of the northern and southern hemispheres (within oceans) than between the Atlantic and Pacific Oceans.     

Genetic evidence for limited trans-Atlantic movements of the harbor porpoise Phocoena phocoena

In the North Atlantic, 14 regional subpopulations have been proposed for the harbor porpoise Phocoena phocoena (Linnaeus, 1758). To what degree these populations experience genetic interchange is relatively unknown, particularly on the larger scale of the North Atlantic as a whole. With the recent completion of several regional genetic studies on population structure of the harbor porpoise in the North Atlantic, came the opportunity to combine datasets in an effort to broaden the geographic scope of focus. Three datasets comprised of mitochondrial DNA control region sequences, representing the Northeast and Northwest Atlantic regions were pooled and reanalyzed to examine the degree of trans-Atlantic exchange among harbor porpoise populations, and to examine the evolutionary history of the species in the North Atlantic. The movement of harbor porpoises across the Atlantic appears to occur at a low level. Genetic variability in the Northeast Atlantic is significantly lower than in the Northwest Atlantic, and may indicate a more recent recolonization for the Northeast Atlantic. The star phylogeny of northeastern haplotypes, with a number of rare haplotypes closely related to the most abundant type, indicates a recent population expansion. A disjunction in haplotypic frequencies between the Northeast and Northwest Atlantic probably occurs east of Greenland, but the exact location and source of the disjunction has yet to be determined. Received: 23 March 1998 / Accepted: 3 December 1998     

Role of wax esters in determining buoyancy inHoplostethus atlanticus (Beryciformes: Trachichthyidae)

Specimens of the orange roughyHoplostethus atlanticus were caught at the edge of the Chatham Rise, east of the South Island of New Zealand, in June 1988, during a commercial fishing voyage. Several tissues of the fish were found to contain large quantities of lipid almost entirely as wax esters. These tissues include a fat-invested swimbladder, the skin, skull and spine, and a fat-infiltrated tissue found in the neurocranial cavity. The physical properties of the wax esters have been examined relative to the temperatures that exist where the fish are found. Differential scanning calorimetry and NMR spectroscopy both suggest that at 6°C, the temperature at which the fish are located in the water column, the lipid is partially solid. However, at 14°C, the temperature at the surface where the fish are positively buoyant, the lipid is completely liquid. It is proposed that the increased density of the lipid at 6°C would make the fish neutrally buoyant in its normal habitat. The location of the lipidrich tissues has been examined relative to their possible buoyancy contribution, and a histological examination of the swimbladder has been performed.     

Evidence of limited gene flow in three species of coral reef fishes in the lagoon of New Caledonia

Allozyme variation was used to investigate the genetic structure of the coral reef fishes Stegastes nigricans, Epinephelus merra and Acanthurus triostegus around New Caledonia. Each species was sampled from each of three sites in the ≃1000 km circumference of the lagoon of New Caledonia. Allelic variation was recorded for each species at 14, 13 and 17 loci, respectively, and heterozygosity diversity (H s) was 0.082, 0.065 and 0.116, respectively. Analysis of genetic differentiation between sites produced inconsistent results between species, with spatial heterogeneity in two species (S. nigricans, F st  = 0.038; A. triostegus, F st = 0.049) and homogeneity in one species (E. merra, F st = 0.000). Hydrological and climatic data from the lagoon suggest that the eastern and western sides of the lagoon are isolated, since they lie in water masses of different origin. This may explain the genetic differentiation and restricted gene flow found at a local scale for S. nigricans and A. triostegus. Homogeneity in populations of E. merra is discussed in relation to its low genetic diversity and its reproductive behaviour. Received: 23 April 1997 / Accepted: 25 September 1997     

Mitochondrial DNA phylogeography of the harbour porpoise Phocoena phocoena in the North Pacific

Genetic structure and phylogeography of the harbour porpoise Phocoena phocoena in the North Pacific were examined using 358 bps sequences from the 5′ end of the mitochondrial DNA control region including those reported previously and newly obtained from the west Pacific. AMOVA and pairwise population φ _st estimates clearly revealed genetic differentiation between an east/south and a north/northwest group with the break along the Pacific Rim at British Columbia. In addition, nested clade phylogeographical analysis, neutrality tests, mismatch distribution analysis, genetic diversities and Mantel test, suggested that the observed genetic structure might have been influenced by contiguous range expansion with restricted gene flow in the direction from south to north along the North American coasts and east to west along the Pacific Rim in the middle to late Pleistocene.     

The structure of subtidal algal and invertebrate assemblages at the Chatham Islands,New Zealand

We examined the distribution and abundance of organisms on subtidal rocky reefs at nine sites around the Chatham Islands, a remote group 780 km east of southern New Zealand. We sampled five depth strata ranging from 1 to<16 m to identify spatial patterns in the abundance of algae and invertebrates and to assess their variation within and among sites. This information is used to discuss hypotheses concerning community structure at this remote locality. Several patterns were apparent. The immediate subtidal was occupied by the southern bull kelp Durvillaea spp. A suite of 11 fucalean species were dominant to a depth of 10 m with an average abundance of 28 m^-2, while one species, Carpophyllum flexuosum, occurred mostly in deeper water. Only two laminarian species of algae were present at the islands. The indigenous Lessonia tholiformis was abundant at 2.5 to 15 m and was not found in deeper water, while the giant kelp Macrocystis pyrifera was abundant at two sites in 12 to 18 m. The commercially valuable abalone Haliotis iris was extremely abundant in shallow water, with an overall mean of 6 m^-2 at 5 m. The sea urchin Evechinus chloroticus was common, but reached high densities only in small (<25 m²) patches. The characteristic urchin-dominated zones reported in kelp beds world-wide were not seen. There was considerable site-to-site variation in the occurrence and abundance of individual species. Some differences between sites were associated with shelter from swell (e.g. M. pyrifera was found only in sheltered sites) and physical habitat (e.g. juvenile H. iris were found only beneath boulders inshore), but much of the variation could not be explained by physical or depth-related factors alone. We hypothesize that the differences in these kelp bed assemblages compared to mainland New Zealand are partially due to the high degree of endemism at the Chatham Islands. Local variation cannot be explained by herbivory, and is most likely the result of the various life-history characteristics of the major habitat-forming species, the large brown algae.     

A genetic analysis of the pea crabs (Decapoda: Pinnotheridae) of New Zealand

Seven populations of the pea crabPinnotheres atrinicola Page were sampled from around the North Island of New Zealand from February to October 1987, and individuals were scored at 23 presumptive enzymatic loci. For a brachyuran crab,P. atrinicola revealed high levels of polymorphism and heterozygosity. Of the loci scored, phosphoglucose isomerase (Gpi) and phosphoglucomutase (Pgm) were distinguished by high variability ( =0.602 and 0.526, respectively). A clinal variation in electromorph frequency was evident at several loci, and atGpi in particular. Statistical analyses revealed that, despite relatively small genetic distance separation, a high degree of structuring was present between the geographic populations. The degree of population subdivision observed in this study is atypical of brachyuran crabs. It is suggested that the genetic differentiation observed between pea crab populations is maintained by life-history attributes and current movements which restrict gene flow between populations and, to some extent, by random genetic drift.Publication No. 41 from the Evolutionary Genetics Laboratory, University of Auckland     

MtDNA population structure and gene flow in lingcod (<Emphasis Type="Italic">Ophiodon elongatus</Emphasis>): limited connectivity despite long-lived pelagic larvae

Lingcod, Ophiodon elongatus Girard, have a 3-month pelagic larval stage and are an important recreational and commercial species on the west coast of North America. Cytochrome-c oxidase I sequences from tissue samples were used to characterize population structure and infer patterns of gene flow from California to Alaska. No significant genetic structure was found when estimates of Wright’s F _ST (i.e., Φ_ST) were generated among all populations sampled. Nesting populations within regions, however, indicated that the inner coast of Washington State is distinct, a result corroborating previous allozyme work. Coalescent-based estimates of gene flow indicate that although migration can be high from an evolutionary perspective, nearly half of all comparisons among populations showed no gene flow in at least one direction. From an ecological perspective, moderate migration rates (Nm < 10) among most populations provide surprisingly limited connectivity at large (∼ 1,000 km) and small (∼100 km) spatial scales. Coalescent-based estimates also show that gene flow between the inner and the outer coasts is asymmetric, a result consistent with prevailing surface currents. Because the expected inter-locus variances for coalescent-based estimates of gene flow are likely large, future work will benefit from analyses of nuclear DNA markers. However, limited demographic connectivity on large spatial scales may help explain why stock recovery has been uneven, with greater recovery in the northern (87% rebuilt) than in the southern (24% rebuilt) fishery region, supporting a regional management strategy. These results suggest that despite a 3-month pelagic larval stage, some areas may be effectively closed with respect to both population dynamics and fishery management issues.     

Phylogeographic structure of <Emphasis Type="Italic">Octopus vulgaris</Emphasis> in South Africa revisited: identification of a second lineage near Durban harbour

In a previous study that investigated genetic structure of Octopus vulgaris along the South African coast by sequencing the mitochondrial cytochrome oxidase III gene (COIII), all sequences generated were identical. Such a finding is unusual, because mitochondrial DNA mutates quickly, and several marine invertebrates present in southern Africa show considerable genetic variation and structure. We reanalysed the samples using two different mitochondrial markers, namely cytochrome oxidase I (COI) and the large ribosomal subunit (16S rRNA). Sequences of both these markers showed variation. The conclusion of the previous study, that South Africa’s O. vulgaris population is characterised by a lack of genetic structure along the coast, is rejected. Some specimens from Durban (southeast Africa) were genetically more different from those found in the remainder of the country than were specimens from other regions (Tristan da Cunha and Senegal). We suggest that the lineage in Durban may have been recently introduced.     

Population structure of the planktonic copepod Calanus pacificus in the North Pacific Ocean

The pelagic copepod Calanus pacificus ranges nearly continuously across temperate-boreal regions of the North Pacific Ocean and is currently divided into three subspecies—C. pacificus oceanicus, C. pacificus californicus, C. pacificus pacificus—based on subtle morphological differences and geographic location. The relation between geography and genetic differentiation was examined for 398 C. pacificus individuals sampled from six widely distributed locations across the North Pacific, including an open ocean site and coastal sites on both sides of the North Pacific basin. For each individual copepod, the DNA sequence was determined for a 421-bp region of the mitochondrial coxI gene (mtCOI). A total of sixty-three different mtCOI sequences, or haplotypes, were detected, with a sequence divergence between haplotypes of 0.2–3.1%. The number and distribution of haplotypes varied with sampling location; 12 haplotypes were distributed across multiple sampling locations, and 51 occurred at only one location. Five genetically distinct populations were detected based on F _ST values. Haplotype minimum spanning networks, nucleotide divergence and F _ST values indicated that individuals from coastal sites in the North Pacific Ocean were more closely related to each other than to individuals from the open ocean site at Station P. These results provide genetic support for the designation of two subspecies—a coastal subspecies that consists of what is currently referred to as C. p. pacificus and C. p. californicus and an open ocean subspecies C. p. oceanicus. This work also indicates that planktonic copepods with potentially high dispersal capacity can develop genetically structured populations in the absence of obvious geographic barriers between proximate locales within an ocean basin.