Abstract The genetic diversity of one of the most abundant species in the Arctic and subarctic oceans, the pelagic snail Limacina helicina, has not yet been characterized in the north Pacific. This species has different ‘forma’ (L. helicina forma helicina, acuta, pacifica and ochotensis), but whether or not the morphological differences between these forma are caused by phenotypic plasticity or genetic differentiation remains unclear. Here, we analysed partial nucleotide sequences of the mitochondrial cytochrome c oxidase subunit I gene in L. helicina from the subarctic western North Pacific Ocean (SWNP; L. helicina f. acuta) and compared them with those from Svalbard (L. helicina f. helicina) and other localities (Beaufort Sea, eastern Pacific, northern Sea of Japan and western Atlantic). The results show clear genetic differentiation between populations in the SWNP and Svalbard (ΦCT = 0.59282, P < 0.001). These genetic differences are consistent with the previous description of the two forma L. h. f. acuta (SWNP) and L. h. f. helicina (Svalbard) based on shell morphology. INTRODUCTION Planktonic organisms (including swimming larvae of benthic organisms) can disperse across relatively wide areas. Their population size, population structure and genetic diversity are influenced by oceanographic events and structural changes in the water mass and currents. For example, many high-latitude marine organisms underwent evolutionary bottleneck events when expansion of polar ice sheets reduced their habitats and relegated them to refugia during the Quaternary ice ages (Hewitt, 2000; Wares & Cunningham, 2001; Marko, 2004; Marko et al., 2010). Subsequent secondary contact between the isolated populations may have increased genetic polymorphisms during interglacial periods (Wilson et al., 2007), allowing some polar organisms to recover genetic diversity. However, it remains unclear whether the low genetic diversity in high-latitude marine species is caused by the small size of former refugia at lower latitudes or by the rapid recolonization of high latitudes and consequent founder effects. Studies of intraspecific genetic variation offer a basic understanding of the biogeographic and/or demographic history of the organism (Avise, 2000) and may be useful for predicting population dynamics in future environments. Global climate change will likely impact the biodiversity of species or populations to varying degrees (Pauls et al., 2013). Responses to changing environmental conditions have been observed as phenotypic plasticity or evolutionary adaptations (Riddle et al., 2008; Scoble & Lowe, 2010; Hoffmann and Sgrò, 2011) and as changes in habitat (Parmesan et al., 1999; Chen et al., 2011). Organisms with low genetic diversity are generally more susceptible to rapid environmental changes than those with high genetic diversity, which also have higher stress resistance (Nowak et al., 2007; Markert, Schelly & Stiassny, 2010) and adaptive potential (Hughes & Stachowicz, 2004). Accordingly, previous studies have demonstrated the importance of genetic diversity for understanding ecological dynamics (Bolnick et al., 2011) and investigations of genetic diversity are considered to be central to understanding the ecosystem impacts of current rapid environmental changes. Ocean acidification caused by increasing atmospheric CO2 dating to the industrial revolution is one of the main environmental stressors in marine ecosystems (Rockström et al., 2009). Because these changes lead to decreased calcium carbonate saturation states (Ω), ongoing ocean chemical changes are expected to have severe effects on marine organisms that have calcifying mechanisms and calcareous structures, such as coccolithophores, foraminiferans, corals, molluscs and sea urchins (Doney et al., 2009). Although the effects of ocean acidification on marine invertebrates have been documented, response patterns differ not only by species but also among populations within species (Fabry et al., 2008). These observations suggest that consideration of genetic background and population structure is needed in order to understand the impact of ocean acidification. The pelagic snail Limacina helicina is widely distributed in the Arctic and subarctic oceans where it can comprise more than 50% of the biomass of polar zooplankton and is one of the key food resources for large zooplankton, fish and seabirds (Lalli & Gilmer, 1989; Hunt et al., 2008, 2010; Comeau et al., 2009). Limacina helicina is a notable species in ocean acidification research (Comeau et al., 2009; Lischka et al., 2011; Lischka & Riebesell, 2012; Bednaršk et al. 2012a, b, 2014; Koh et al., 2015; Johnson & Hofmann, 2016; Peck et al., 2016). It is classified as four ‘forma’ in the northern hemisphere (North Atlantic: L. h. f. helicina Phipps, 1774; western North Pacific: L. h. f. acutaSpoel, 1967; eastern North Pacific: L. h. f. pacifica Dall, 1871; Okhotsk Sea: L. h. f. ochotensis Shkoldina, 1999; Spoel, 1967: 257, 349), defined by shell morphology and habitat (Spoel, 1967: 36–43, 257, 349; Shkoldina, 1999), but their genetic diversity and genetic population structure have not been clarified. Hunt et al. (2010) found that large genetic differences exist between Arctic and Antarctic populations of L. ‘helicina’ and the two morphotypes helicina and antarctica are now considered to be distinct species, L. helicina (Phipps, 1774) and L. antarctica Woodward, 1854 (for which the valid name is L. rangii (d’Orbigny, 1834); World Register of Marine Species, http://www.marinespecies.org accessed 6 October 2017). Sromek, Lasota & Wolowicz (2015) performed population genetic analyses of the mitochondrial cytochrome c oxidase subunit I (COI) gene of L. helicina s. l. and showed lower genetic diversity in the Arctic species (haplotype diversity, H = 0.778; nucleotide diversity, π = 0.0031) than in the Antarctic species (H = 0.967, π = 0.0061). However, L. helicina individuals from only five localities around Svalbard were used in the comparison, so that the resulting data may not be sufficient to understand the genetic structure of this species. Here, we examined more samples from the subarctic western North Pacific (SWNP) in order to clarify further the population structure and genetic diversity of L. helicina. MATERIAL AND METHODS Samples A total of 101 individuals of Limacina helicina f. acuta were collected at five stations (Tsugaru, KNOT, K2, ESM, NWP) in the subarctic western North Pacific Ocean (SWNP; Fig. 1, Table 1) during cruises of the R/V Hakuho Maru (KH15-J01, K2), T/S Ushio Maru (US348, Tsugaru), T/S Osyoro Maru (OS026-Leg2, KNOT and ESM) and R/V Yokosuka (YK16-16, NWP). Samples were obtained from vertical hauls from 150 m to the surface with a ring net (45 cm in diameter, 100-μm mesh). Specimens were fixed in 100% ethanol and stored at −20 °C. Table 1. Station and sampling information for the subarctic western North Pacific population of Limacina helicina. Name Latitude Longitude Date N GenBank Accession Nos Tsugaru 41°83′N 141°22′E 17 Nov. 2015 16 LC185064-LC185073 LC229751-LC229757 KNOT 43°11′N 153°19′E 4–8 July 2016 30 LC185015-LC185032 LC229727-LC229734 LC229766-LC229769 K2 47°00′N 160°00′E 4 July 2015 19 LC185045-LC185063 ESM 44°11′N 170°24′E 7 July 2016 20 LC185033-LC185044 LC229758-LC229765 NWP 41°83′N 141°22′E 13 Nov. 2016 16 LC229735-LC229750 Name Latitude Longitude Date N GenBank Accession Nos Tsugaru 41°83′N 141°22′E 17 Nov. 2015 16 LC185064-LC185073 LC229751-LC229757 KNOT 43°11′N 153°19′E 4–8 July 2016 30 LC185015-LC185032 LC229727-LC229734 LC229766-LC229769 K2 47°00′N 160°00′E 4 July 2015 19 LC185045-LC185063 ESM 44°11′N 170°24′E 7 July 2016 20 LC185033-LC185044 LC229758-LC229765 NWP 41°83′N 141°22′E 13 Nov. 2016 16 LC229735-LC229750 Figure 1. View largeDownload slide Sampling locations for Limacina helicina in the subarctic western North Pacific Ocean. Figure 1. View largeDownload slide Sampling locations for Limacina helicina in the subarctic western North Pacific Ocean. Micro X-ray CT The shell morphology of individuals L. h. f. acuta collected from K2 was visualized using a Microfocus X-ray CT scanner (ScanXmateD160TSS105, Comscantechno Co., Shin-Yokohama, Japan) at Japan Agency for Marine-Earth Science and Technology (JAMSTEC). High-resolution settings (X-ray focus diameter 0.8 μm, voltage 80 KeV, detector array size 1024 × 1024 pixels, 1500 projections in 360° rotations) were applied for 3D quantitative densitometry of individual shells. Geometric resolution of isotropic voxel size was varied from 2 to 4 μm depending on each shell size. We used ConeCTexpress (White Rabbit Corp., Tokyo, Japan) for corrections and reconstruction of tomography data, and the general principle of Feldkamp cone beam reconstruction was followed to reconstruct the image cross-sections based on filtered back-projections. We obtained 2D images from 3D data and measured shell height (H) and diameter (D) using ImageJ v. 1.50i. Morphometric analysis We attempted to evaluate two morphological aspects of coiling patterns and aperture shapes of shells, based on theoretical morphological and geometric-morphometric approaches. The coiling patterns were approximated using the growing tube model (Okamoto, 1988a, b, c) and patterns were characterized by three parameters of the model (see Okamoto, 1988a for details): expansion rate (EG), standardized curvature (CG) and standardized torsion (TG). We estimated the parameters of the growing tube model using Raup’s model, following a method proposed by Noshita (2014). Raup’s model is a theoretical morphological model for describing coiling patterns based on three parameters (Raup 1962, 1966, 1967; Raup & Michelson, 1965): whorl expansion rate (WR), whorl translation rate (TR) and relative width of the umbilicus (DR). These parameters were estimated from 2D images. We analysed the shape of the helicocone bisected by a plane passing through the coiling axis instead of using the changes in aperture shapes through growth stages. These ‘aperture shapes’ were quantified by elliptic Fourier descriptors (EFDs) (Kuhl & Giardina, 1982). To estimate the parameters of the growing tube model and EFDs, we developed the following simplified procedure. First, we resampled a slice image passing through the coiling axis from a 3D CT model with Amira® (FEI). Second, centroids, areas and outlines of ‘apertures’ were digitized from the slice image using the image-processing software Fiji in ImageJ (Schneider, Rasband & Eliceiri, 2012). Third, Raup’s W, T and D (WR, TR, DR) were estimated from areas and centroid coordinate values in a specimen by fitting Raup’s model describing the growth trajectory. Fourth, Raup’s parameters were converted into the parameters of the growing tube model based on an equation by Noshita (2014). Finally, the outlines of ‘apertures’ were quantified with EFDs using Momocs, which is an R package for outline-based geometric morphometrics (Bonhomme et al., 2014). In this way, our simple procedure could be used to evaluate the coiling patterns and aperture shapes of L. h. f. acuta shells and to provide quantitative information characterizing the changes in growth. Coiling patterns and aperture shape analysis For coiling patterns, a pairwise multivariate analysis of variance (MANOVA) of the parameters of the growing tube model was performed to evaluate the effect of growth stage (first half and second half). For analysis of aperture shapes, principal component analysis (PCA) was used to contrast morphological variation in terms of the EFDs of aperture shapes. A pairwise MANOVA of the scores of PC1 to PC18, for which the cumulative contribution ratio reached 99%, was performed to evaluate the effect of growth stage, first half (D = 687–814 μm) vs. second half. In both analyses, Pillai’s trace test was used for this purpose (Hand & Taylor, 1987). To avoid multiple testing problems, P-values were adjusted by the Holm method (Holm, 1979). DNA extraction, amplification and sequencing DNA extraction was performed using GeneReleaser (BioVentures, Murfreesboro, TN) in accordance with the manufacturer’s protocol. Gene sequences for the COI gene were amplified by PCR using the universal primers (LCO1490 and HCO2198) as described by Folmer et al. (1994). PCR products were purified using a commercial kit (ExoSAP-IT for PCR product clean-up; Thermo Fisher Scientific, Waltham, MA). Sequencing reactions were performed using ABI BigDye Terminators v. 3.1 (Applied Biosystems, Foster City, CA) and purified by ethanol precipitation, and sequencing products were sequenced using an ABI 3130 genetic analyser with the same primers as for the PCR. Sequences were deposited in GenBank (for accession numbers see Table 1). Molecular phylogenetic analysis We obtained 101 COI sequences of L. helicina from SWNP and retrieved an additional 44 sequences of L. helicina f. helicina, L. antarctica, L. helicoidea, L. inflata, L. retroversa and Hyalocylis striata obtained from GenBank (Supplementary Material Table S1) for phylogenetic analysis. Sequence alignments were conducted using the online version of MAFFT (v. 7.310; http://mafft.cbrc.jp/alignment/server/index.html; Katoh et al., 2002) based on 503 bp of COI. The best-fit nucleotide substitution model with the lowest Bayesian information criterion (BIC) score was selected using MEGA v. 5.1 (Tamura et al., 2011). Maximum-likelihood (ML) trees were inferred with MEGA v. 5.1 using Tamura’s 3-parameter model + G + I (Tamura, 1992) with 1,000 bootstrap replications. Sequences of H. striata were used as the outgroup. Population genetic analysis The 101 sequences gathered by this study were combined with 77 from Svalbard and other regions derived from previous studies (Supplementary Material Fig. S1, Table S2; Hunt et al., 2010; Jennings et al., 2010; Corse et al., 2013; Layton, Martel, Hebert, 2014; Sromek et al., 2015; Chichvarkhin, 2016). Sequence alignment was conducted using the online version of MAFFT (Katoh et al., 2002) and sequences were trimmed to 503 bp using MEGA v. 5.1. Population genetic analyses (diversity indices and neutrality test) were performed using Arlequin v. 220.127.116.11 (Excoffier & Lischer, 2010). Pairwise ΦST (among localities) and ΦCT (between SWNP and Svalbard) with AMOVA were performed by Arlequin and significance tested by 10,000 permutations. Parsimony networks were constructed using TCS v. 1.21 software with the connection probability set at 95% (Clement, Posada & Crandall, 2000). RESULTS Shell morphology in SWNP D and H of Limacina helicina f. acuta shells from K2 were 3.168 and 3.153 mm, respectively (Supplementary Material Table S3). The average shell height to diameter ratio (H/D) was 0.841 ± 0.080 (SD) (Fig. 2A). The shape of small individuals with D < 0.8 mm differed from that of the larger ones (small: H/D = 0.806 ± 0.036; large: H/D = 0.987 ± 0.196). We then analysed the aperture shapes and growth patterns using 3D data and compared small and large individuals (i.e. early and late stages measured on the same individuals). We found that aperture shapes and growth patterns differed significantly between early and late stages of large individuals (Fig. 2B, C and Supplementary Material Table S4), but not between small individuals and the early stage of large individuals (Fig. 2B, C and Supplementary Material Table S4). Figure 2. View largeDownload slide Shell morphology of Limacina helicina f. acuta at locality K2. A. Shell height to shell diameter ratio (H/D). B. Cube plot of the three estimated parameters of the growing tube model (EG, expansion rate; CG, standardized curvature, TG, standardized torsion). ‘Small’ indicates D < 1 mm. C. Principal component analysis of aperture shape by size and stage. Figure 2. View largeDownload slide Shell morphology of Limacina helicina f. acuta at locality K2. A. Shell height to shell diameter ratio (H/D). B. Cube plot of the three estimated parameters of the growing tube model (EG, expansion rate; CG, standardized curvature, TG, standardized torsion). ‘Small’ indicates D < 1 mm. C. Principal component analysis of aperture shape by size and stage. Molecular phylogeny, phylogeography and population genetic analysis Phylogenetic analysis using the ML method placed all COI sequences of L. helicina in a single clade (Fig. 3). There were 24 variable sites (4.77%) among the 503 bp of the COI sequences. H and π were markedly lower for the SWNP population than for the Svalbard population (Table 2). Based on analyses of pooled samples from SWNP and other regions, we constructed a haplotype network that showed two major haplotypes (H1 and H2) and multiple singletons (Fig. 4). Haplotype H1 was shared by samples from the northern Sea of Japan, eastern North Pacific, Hudson Bay, the Canadian Arctic Ocean and Svalbard, while Haplotype H2 (two mutational steps from H1) was restricted to the Svalbard population (Fig. 4). Most samples from SWNP were haplotype H1 (75/101, 74.2%), while most individuals from Svalbard were haplotype H2 (32/67, 47.76%). Significant population differences were found between SWNP (Tsugaru, KNOT, K2, ESM and NWP) and Svalbard based on pairwise ΦST and ΦCT comparisons (ΦCT = 0.59282, P < 0.001); (Table 3). The haplotype network shows a typical star-like structure with multiple singletons and Tajima’s D (Tajima, 1989) and Fu’s Fs neutrality tests (Fu, 1997) were significantly negative (Table 2). Table 2. Population genetic analyses (diversity indices and neutrality tests) of populations of Limacina helicina from subarctic western North Pacific and Svalbard, based on partial COI sequence data. Subarctic western North Pacific Svalbard Number of individuals 101 67 Number of sites 503 503 Number of haplotypes (k) 25 24 Polymorphic sites (S) 24 25 Haplotype diversity (H) 0.450 (±0.064) 0.764 (±0.054) Nucleotide diversity (π) 0.00121 (±0.00021) 0.00324 (±0.00049) Average number of nucleotide differences (Π) 0.610 1.629 Tajima’s D −2.60025*** −2.268** Fu’s Fs −40.848*** −23.916*** Subarctic western North Pacific Svalbard Number of individuals 101 67 Number of sites 503 503 Number of haplotypes (k) 25 24 Polymorphic sites (S) 24 25 Haplotype diversity (H) 0.450 (±0.064) 0.764 (±0.054) Nucleotide diversity (π) 0.00121 (±0.00021) 0.00324 (±0.00049) Average number of nucleotide differences (Π) 0.610 1.629 Tajima’s D −2.60025*** −2.268** Fu’s Fs −40.848*** −23.916*** *P < 0.05; **P < 0.01; ***P < 0.001. Table 3. Population differentiation in Limacina helicina based on pairwise ΦST values calculated from COI data. Tsugaru KNOT K2 ESM NWP Svalbard Tsugaru − KNOT 0.01186 − K2 0.00120 −0.00479 − ESM −0.00483 −0.00526 −0.01430 − NWP 0.00000 −0.02961 −0.00117 −0.02796 − Svalbard 0.47470*** 0.52570*** 0.51465*** 0.49440*** 0.50745*** − Tsugaru KNOT K2 ESM NWP Svalbard Tsugaru − KNOT 0.01186 − K2 0.00120 −0.00479 − ESM −0.00483 −0.00526 −0.01430 − NWP 0.00000 −0.02961 −0.00117 −0.02796 − Svalbard 0.47470*** 0.52570*** 0.51465*** 0.49440*** 0.50745*** − *P < 0.05; **P < 0.01; ***P < 0.001. For sampling localities see Figure 1 and Supplementary Material Figure S1. Figure 3. View largeDownload slide Molecular phylogenetic analysis of Limacina species based on 503 bp of COI gene. A. Maximum-likelihood trees inferred using Tamura’s 3-parameter model + G + I (Tamura, 1992) with 1,000 bootstrap replications. B. Details of L. helicina clade. Open circles indicate Svalbard samples (from Sromek et al., 2015) and filled circles the subarctic western North Pacific (SWNP) samples (this study). H1 and H2 indicate haplotypes H1 and H2 (Fig. 4). Figure 3. View largeDownload slide Molecular phylogenetic analysis of Limacina species based on 503 bp of COI gene. A. Maximum-likelihood trees inferred using Tamura’s 3-parameter model + G + I (Tamura, 1992) with 1,000 bootstrap replications. B. Details of L. helicina clade. Open circles indicate Svalbard samples (from Sromek et al., 2015) and filled circles the subarctic western North Pacific (SWNP) samples (this study). H1 and H2 indicate haplotypes H1 and H2 (Fig. 4). Figure 4. View large Download slide Parsimony haplotype network of COI sequences from Limacina helicina. Haplotype circle sizes indicate frequency and colours indicate sampling location (Supplementary Material Fig. S1, Table S2). SWNP, subarctic western North Pacific; NSJ, northern Sea of Japan; ENP, eastern North Pacific; HB, Hudson Bay; CAO, Canadian Arctic Ocean. Figure 4. View large Download slide Parsimony haplotype network of COI sequences from Limacina helicina. Haplotype circle sizes indicate frequency and colours indicate sampling location (Supplementary Material Fig. S1, Table S2). SWNP, subarctic western North Pacific; NSJ, northern Sea of Japan; ENP, eastern North Pacific; HB, Hudson Bay; CAO, Canadian Arctic Ocean. DISCUSSION Based on morphological descriptions, four different forma of Limacina helicina have been identified in subarctic and arctic regions (North Atlantic: L. h. f. helicina; western North Pacific: L. h. f. acuta; eastern North Pacific: L. h. f. pacifica; Okhotsk Sea: L. h. f. ochotensis; Spoel, 1967). Whether the morphological differences between these forma are the result of phenotypic plasticity or genetic differentiation is unknown. In this study we characterized the shell ratio H/D of L. helicina collected from SWNP (K2). Two of the forma of L. helicina are recorded from the North Pacific, L. h. f. acuta and L. h. f. pacifica, and are described as having a high-spired and low-spired shell, respectively. The former is known to occur in the SWNP (McGowan, 1963). Our samples could be categorized by shell shape into two groups: large individuals with high-spired shells and small individuals with low-spired shells (Fig. 2A). The large individuals are consistent with previous descriptions of L. h. f. acuta (H/D = 0.85–1.2; D = 0.7–1.8 mm; McGowan, 1963). Our analysis of aperture shapes and growth patterns revealed a change in aperture shape and growth pattern between early and late stages of the same shells, such that the early stage of large individuals showed aperture shapes and growth patterns similar to those of small individuals (Fig. 2B, C). These results suggest that the two different forms found at K2 were simply early and late stages of a common ontogenetic pattern. Most samples from SWNP were haplotype H1 and most from Svalbard were haplotype H2, indicating significant genetic differentiation between the two populations (Table 3). In addition, our SWNP samples had a shell shape different from that of L. h. f. helicina in the North Atlantic Ocean (H/D = 0.75; H up to 6 mm, D up to 8 mm; Spoel, 1967). The genetic results are thus consistent with the previous establishment of two distinct forma, L. h. f. acuta and L. h. f. helicina. Although few COI sequences of L. helicina have been reported from the northern Sea of Japan, eastern North Pacific, Hudson Bay and the Canadian Arctic Ocean, most are of haplotype H1 (Fig. 4). This suggests the possibility of dispersal from the SWNP to the Arctic Ocean via the Bering Strait (McLaughlin et al., 1996). Similarly, the copepod Pseudocalanus newmani shows no regional differentiation between the Gulf of Alaska and the Chukchi Sea, and northward directional gene flow from the Gulf of Alaska to the Beaufort Sea (Questel et al. 2016). Despite the significant differentiation of the SWNP and Svalbard populations of L. helicina, a few of the Svalbard individuals have haplotype H1 (6/67). This could be explained by dispersal among the Chukchi Sea, Beaufort Sea and around Svalbard by the Arctic intermediate layer circulation (Jones, 2001). Genetic data for individuals from the Chukchi and Beaufort Seas is required to test this hypothesis. Our haplotype network analysis showed a typical star-like structure (Fig. 4) and Tajima’s D (Tajima, 1989) and Fu’s Fs neutrality tests (Fu, 1997) were both significantly negative for samples from both the SWNP and Svalbard (Table 2). This suggests a recent dramatic population expansion in both regions, as has been reported for L. h. f. helicina (Sromek et al., 2015) and other high-latitude organisms (Bernatchez & Wilson, 1998; Marko, 2004; Maggs et al., 2008; Grant et al., 2011). These patterns are also consistent with population bottlenecks caused by restriction of Arctic marine organisms to small southern refugia during the ice ages and their subsequent population expansions (Hewitt, 2004; Hardy et al., 2010). MtDNA diversity in L. h. f. acuta is very low, lower than that of L. h. f. helicina (Table 2; Sromek et al., 2015), which may reflect differences in life history and physical dispersal abilities. Sromek et al. (2015) estimated that L. h. f. helicina in Svalbard showed increased population size and genetic diversity at 131 kyr BP (confidence interval 44–225 kyr BP) and this was related to the Eemian interglacial period, based on mismatch distribution analyses. This agrees with simulations of historical changes in the continent ice sheets of Greenland and North American by Abe-Ouchi et al. (2013), which showed shrinkage of ice sheets, ending the isolation between the SWNP, the Arctic Ocean and Svalbard, at 131 ky BP. This and similar geological events likely triggered secondary contact between the SWNP and Svalbard populations (via the Chukchi and Beaufort Seas), leading to an increase in genetic polymorphism during interglacial periods. Importantly, the current through the Bering Strait is in a northward direction, from the SWNP to the Arctic Ocean (McLaughlin et al., 1996), and this directionality permitted genetic divergence between the SWNP population and the Svalbard population. This pattern is not universal. For example, the high-latitude planktonic foraminiferan Neogloboquadrina pachyderma shows evidence of a geographical barrier between the Bering Sea and the Chukchi Sea (Darling, Kucera, Wade, 2007), while some rocky-shore species show significant genetic differentiation even within the northeastern Pacific itself (Marko et al., 2010). These contrasting patterns reflect the dispersal characteristics and habitats of the organisms inhabiting the same sea region. Invertebrates have been reported to show different responses to ocean acidification, not only between species but also between populations within a species (Fabry et al., 2008). Thus, these forma L. h. helicina and L. h. acuta might respond differently to ocean acidification and we must take genetic differences into careful consideration in future studies of responses to environmental change. In the SWNP (K2 and KNOT), wintertime surface waters have a lower calcium carbonate saturation state (Ω about 1.4) and a shallower CaCO3 saturation horizon (about 120 m; Wakita et al., 2013) than at other monitoring stations (Bates et al., 2014). In addition, acidification in winter in this region occurs at a slower rate than the annual mean rate because of a reduced rate of increase of dissolved inorganic carbon and an increase of total alkalinity (Wakita et al., 2017). This increase of total alkalinity is possibly caused by the weakening of calcification by organisms in winter. Thus, the SWNP is expected to be among the regions in which organisms will be most seriously impacted (Orr et al., 2005). It is recommended that future studies of L. helicina and other organisms at high risk in the SWNP should investigate genetic structure in more detail using variable markers such as microsatellites or SNPs regions (i.e. RADseq, Miller et al., 2007 and MIGseq, Suyama & Matsuki, 2015). SUPPLEMENTARY MATERIAL Supplementary material is available at Journal of Molluscan Studies online. ACKNOWLEDGEMENTS We thank Yoshiyuki Abe, Souta Komeda and Marie Maekakuchi (Hokkaido University) for the collection of Limacina helicina. We wish to thank Yuriko Nakamura (JAMSTEC) for help in micro X-ray CT analysis. We thank Minoru Kitamura, Naomi Harada and Masashi Tsuchiya (JAMSTEC) for invaluable discussion and suggestion. We would like to thank the captains, crews and cruise members (R/V Hakuho Maru, KH-15-J01; R/V Yokosuka, YK16-16; T/S Ushio Maru, US348; T/S Oshoro Maru, OS026-Leg2). This study was supported by JSPS Grants-in-Aid for Scientific Research 15H06908 and 16H04961. REFERENCES Abe-Ouchi, A., Saito, F., Kawamura, K., Raymo, M., Okuno, J., Takahashi, K. & Blatter, H. 2013. Insolation driven 100,000-year glacial cycles and hysteresis of ice sheet volume. Nature , 500: 190– 193. Google Scholar CrossRef Search ADS PubMed Avise, J.C. 2000. Phylogeography: the history and formation of species . Harvard University Press, Cambridge, MS. Bates, N.R., Astor, Y.M., Church, M.J., Currie, K., Dore, J.E., González-Dávila, M., Laura, L., Muller-Karger, F., Olaffson, J. & Santana-Casiano, J.M. 2014. A time-series view of changing ocean chemistry due to ocean uptake of anthropogenic CO2 and ocean acidification. Oceanography , 27: 126– 141. Google Scholar CrossRef Search ADS Bednaršk, N., Tarling, G.A., Bakker, D.C., Fielding, S., Cohen, A., Kuzirian, A., Mccorkle, D., Leze, B. & Montagna, R. 2012a. Description and quantification of pteropod shell dissolution: a sensitive bioindicator of ocean acidification. Global Change Biology , 18: 2378– 2388. Google Scholar CrossRef Search ADS Bednaršek, N., Tarling, G.A., Bakker, D.C.E., Fielding, S., Jones, E.M., Venables, H.J., Ward, P., Kuzirian, A., Leze, B., Feely, R.A. & Murphy, E.J. 2012b. Extensive dissolution of live pteropods in the Southern Ocean. Nature Geoscience , 5: 881– 885. Google Scholar CrossRef Search ADS Bednaršek, N., Feely, R.A., Reum, J.C.P., Peterson, W., Menkel, J., Alin, S.R. & Hales, B. 2014. Limacina helicina shell dissolution as an indicator of declining habitat suitability due to ocean acidification in the California Current Ecosystem. Proceedings of the Royal Society of London, Series B , 281: 20140123. Google Scholar CrossRef Search ADS Bernatchez, L. & Wilson, C.C. 1998. Comparative phylogeography of Nearctic and Palearctic fishes. Molecular Ecology , 7: 431– 452. Google Scholar CrossRef Search ADS Bolnick, D.I., Amarasekare, P., Araujo, M.S., Burger, R., Levine, J.M., Novak, M., Rudolf, V.H.W., Schreiber, S.J., Urban, M.C. & Vasseur, D.A. 2011. Why intraspecific trait variation matters in community ecology. Trends in Ecology and Evolution , 26: 183– 192. Google Scholar CrossRef Search ADS PubMed Bonhomme, V., Picq, S., Gaucherel, C. & Claude, J. 2014. Momocs: outline analysis using R. Journal of Statistical Software , 56: 1– 24. Google Scholar CrossRef Search ADS Chen, I.C., Hill, J.K., Ohlemüller, R., Roy, D.B. & Thomas, C.D. 2011. Rapid range shifts of species associated with high levels of climate warming. Science , 333: 1024– 1026. Google Scholar CrossRef Search ADS PubMed Chichvarkhin, A. 2016. Shallow water sea slugs (Gastropoda: Heterobranchia) from the northwestern coast of the Sea of Japan, north of Peter the Great Bay, Russia. PeerJ , 4: e2774. Google Scholar CrossRef Search ADS PubMed Clement, M., Posada, D. & Crandall, K. 2000. TCS: a computer program to estimate gene genealogies. Molecular Ecology , 9: 1657– 1660. Google Scholar CrossRef Search ADS PubMed Comeau, S., Gorsky, G., Jeffree, R., Teyssié, J.L. & Gattuso, J.P. 2009. Impact of ocean acidification on a key arcticpelagic mollusc (Limacina helicina). Biogeosciences , 6: 1877– 1882. Google Scholar CrossRef Search ADS Corse, E., Rampal, J., Cuoc, C., Pech, N., Perez, Y. & Gilles, A. 2013. Phylogenetic analysis of Thecosomata Blainville, 1824 (Holoplanktonic Opisthobranchia) using morphological and molecular data. PLoS One , 8: e59439. Google Scholar CrossRef Search ADS PubMed Darling, K.F., Kucera, M. & Wade, C.M. 2007. Global molecular phylogeography reveals persistent Arctic circumpolar isolation in a marine planktonic protist. Proceedings of the National Academy of Sciences of the USA , 104: 5002– 5007. Google Scholar CrossRef Search ADS PubMed Doney, S.C., Fabry, V.J., Feely, R.A. & Kleypas, J.A. 2009. Ocean acidification: the other CO2 problem. Annual Review of Marine Science , 1: 169– 192. Google Scholar CrossRef Search ADS PubMed Excoffier, L. & Lischer, H.E.L. 2010. Arlequin suite ver. 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources , 10: 564– 567. Google Scholar CrossRef Search ADS PubMed Fabry, V.J., Seibel, B.A., Feely, R.A. & Orr, J.C. 2008. Impacts of ocean acidification on marine fauna and ecosystem processes. ICES Journal of Marine Science , 65: 414– 432. Google Scholar CrossRef Search ADS Folmer, O., Black, M., Hoeh, W., Lutz, R. & Vrijenhoek, R. 1994. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Molecular Marine Ecology and. Biotechnology , 3: 294– 299. Fu, Y.X. 1997. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics , 147: 915– 925. Google Scholar PubMed Grant, W.S., Merkouris, S.E., Kruse, G.H. & Seeb, L.W. 2011. Low allozyme heterozygosity in North Pacific and Bering Sea populations of red king crab (Paralithodes camtschaticus): adaptive specialization, population bottleneck, or metapopulation structure? ICES Journal of Marine Science , 68: 499– 506. Google Scholar CrossRef Search ADS Hand, D.J. & Taylor, C.C. 1987. Multivariate analysis of variance and repeated measures: a practical approach for behavioural scientists . Chapman & Hall, London. Google Scholar CrossRef Search ADS Hardy, S. M., Carr, C. M., Hardman, M., Steinke, D., Corstorphine, E. & Mah, C. 2011. Biodiversity and phylogeography of Arctic marine fauna: insights from molecular tools. Marine Biodiversity , 41: 195– 210. Google Scholar CrossRef Search ADS Hewitt, G. 2000. The genetic legacy of the Quaternary ice ages. Nature , 405: 907– 913. Google Scholar CrossRef Search ADS PubMed Hewitt, G.M. 2004. Genetic consequences of climatic oscillations in the Quaternary. Philosophical Transactions of the Royal Society of London, Series B , 359: 183– 195. Google Scholar CrossRef Search ADS Hoffmann, A.A. & Sgrò, C.M. 2011. Climate change and evolutionary adaptation. Nature , 70: 479– 485. Google Scholar CrossRef Search ADS Holm, S. 1979. A simple sequentially rejective multiple test procedure. Scandinavian Journal of Statistics , 6: 65– 70. Hughes, A.R. & Stachowicz, J.J. 2004. Genetic diversity enhances the resistance of a seagrass ecosystem to disturbance. Proceedings of the National Academy of Sciences of the USA , 101: 8998– 9002. Google Scholar CrossRef Search ADS PubMed Hunt, B., Strugnell, J., Bednarsek, N., Linse, K., Nelson, R.J., Pakhomov, E., Seibel, B., Steinke, D. & Würzberg, L. 2010. Poles apart: the ‘bipolar’ pteropod species Limacina helicina is genetically distinct between the Arctic and Antarctic oceans. PLoS One , 5: e9835. Google Scholar CrossRef Search ADS PubMed Hunt, B.P.V., Pakhomov, E.A., Hosie, G.W., Siegel, V., Ward, P. & Bernard, K. 2008. Pteropods in southern ocean ecosystems. Progress in Oceanography , 78: 193– 221. Google Scholar CrossRef Search ADS Jennings, R.M., Bucklin, A., Ossenbrügger, H. & Hopcroft, R.R. 2010. Species diversity of planktonic gastropods (Pteropoda and Heteropoda) from six ocean regions based on DNA barcode analysis. Deep-Sea Research II , 57: 2199– 2210. Google Scholar CrossRef Search ADS Johnson, K. M. & Hofmann, G. E. 2016. A transcriptome resource for the Antarctic pteropod Limacina helicina antarctica. Marine Genomics , 28: 25– 28. Google Scholar CrossRef Search ADS PubMed Jones, E.P. 2001. Circulation in the arctic ocean. Polar Research , 20: 139– 146. Google Scholar CrossRef Search ADS Katoh, K., Misawa, K., Kuma, K. & Miyata, T. 2002. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Research , 30: 3059– 3066. Google Scholar CrossRef Search ADS PubMed Koh, H. Y., Lee, J. H., Han, S. J., Park, H., Shin, S. C. & Lee, S. G. 2015. A transcriptomic analysis of the response of the arctic pteropod Limacina helicina to carbon dioxide-driven seawater acidification. Polar Biology , 38: 1727– 1740. Google Scholar CrossRef Search ADS Kuhl, F.P. & Giardina, C.R. 1982. Elliptic Fourier features of a closed contour. Computer Graphics and Image Processing , 18: 236– 258. Google Scholar CrossRef Search ADS Lalli, C.M. & Gilmer, R.W. 1989. Pelagic snails. The biology of holoplanktonic gastropod mollusks . Stanford University Press, Stanford, CA. Lischka, S., Büdenbender, J., Boxhammer, T. & Riebesell, U. 2011. Impact of ocean acidification and elevated temperatures on early juveniles of the polar shelled pteropod Limacina helicina: mortality, shell degradation, and shell growth. Biogeosciences , 8: 919– 932. Google Scholar CrossRef Search ADS Lischka, S & Riebesell, U. 2012. Synergistic effects of ocean acidification and warming on overwintering pteropods in the Arctic. Global Change Biology , 18: 3517– 3528. Google Scholar CrossRef Search ADS Layton, K. K., Martel, A. L. & Hebert, P. D. 2014. Patterns of DNA barcode variation in Canadian marine molluscs. PLoS One , 9: e95003. Google Scholar CrossRef Search ADS PubMed Maggs, C.A., Castilho, R., Foltz, D., Henzler, C., Jolly, M.T., Kelly, J., Olsen, J., Perez, K.E., Stam, W., Vainola, R., Viard, F. & Wares, J. 2008. Evaluating signatures of glacial refugia for North Atlantic benthic marine taxa. Ecology , 89: sp11. Google Scholar CrossRef Search ADS Markert, J.A., Schelly, R.C. & Stiassny, M. 2010. Genetic isolation and morphological divergence mediated by high-energy rapids in two cichlid genera from the lower Congo rapids. BMC Evolutionary Biology , 10: 149. Google Scholar CrossRef Search ADS PubMed Marko, P.B. 2004. ‘What’s larvae got to do with it?’ Disparate patterns of post-glacial population structure in two benthic marine gastropods with identical dispersal potential. Molecular Ecology , 13: 597– 611. Google Scholar CrossRef Search ADS PubMed Marko, P. B., Hoffman, J. M., Emme, S. A., McGovern, T. M., Keever, C. C. & Nicole Cox, L. 2010. The ‘Expansion–Contraction’ model of Pleistocene biogeography: rocky shores suffer a sea change? Molecular Ecology , 19: 146– 169. Google Scholar CrossRef Search ADS PubMed Mcgowan, J.A. 1963. Geographical variation in Limacina helicina in the North Pacific. In: Speciation in the sea , ( J.P. Harding & N. Tebble, eds), pp. 109– 128. Publications of the Systematics Association 5. McLaughlin, F.A., Carmack, E.C., Macdonald, R.W. & Bishop, J. 1996. Physical and geochemical properties across the Atlantic/Pacific water mass front in the southern Canadian Basin. Journal of Geophysics Research. , 101: 1183– 1197. Google Scholar CrossRef Search ADS Miller, M. R., Dunham, J. P., Amores, A., Cresko, W. A. & Johnson, E. A. 2007. Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers. Genome Research , 17: 240– 248. Google Scholar CrossRef Search ADS PubMed Noshita, K. 2014. Quantification and geometric analysis of coiling patterns in gastropod shells based on 3D and 2D image data. Journal of Theoretical Biolopgy , 363: 93– 104. Google Scholar CrossRef Search ADS Nowak, C., Jost, D., Vogt, C., Oetken, M., Schwenk, K. & Oehlmann, J. 2007. Consequences of inbreeding and reduced genetic variation on tolerance to cadmium stress in the midge Chironomus riparius. Aquatic Toxicology , 85: 278– 284. Google Scholar CrossRef Search ADS PubMed Okamoto, T. 1988a. Analysis of heteromorph ammonoids by differential geometry. Palaeontology , 31: 35– 52. Okamoto, T. 1988b. Changes in life orientation during the ontogeny of some heteromorph ammonoids. Palaeontology , 31: 281– 294. Okamoto, T. 1988c. Developmental regulation and morphological saltation in the heteromorph ammonite Nipponites. Paleobiology , 14: 272– 286. Google Scholar CrossRef Search ADS Orr, J.C., Fabry, V.J., Aumont, O., Bopp, L., Doney, S.C., Feely, R.A., Gnanadesikan, A, Gruber, N, Ishida, A, Joos, F, Key, RM, Lindsay, K, Maier-Reimer, E, Matear, R, Monfray, P, Mouchet, A, Najjar, RG, Plattner, GK, Rodgers, KB, Sabine, CL, Sarmiento, JL, Schlitzer, R, Slater, RD, Totterdell, IJ, Weirig, MF, Yamanaka, Y & Yool, A 2005. Anthropogenic ocean acidification over the twenty-first century and its impact on calcifying organisms. Nature , 437: 681– 686. Google Scholar CrossRef Search ADS PubMed Parmesan, C., Ryrholm, N., Stefanescu, C. & Hill, J.K. 1999. Poleward shifts in geographical ranges of butterfly species associated with regional warming. Nature , 399: 579– 583. Google Scholar CrossRef Search ADS Pauls, S.U., Nowak, C., Bálint, M. & Pfenninger, M. 2013. The impact of global climate change on genetic diversity within populations and species. Molecular Ecology , 22: 925– 946. Google Scholar CrossRef Search ADS PubMed Peck, V. L., Tarling, G. A., Manno, C., Harper, E. M. & Tynan, E. 2016. Outer organic layer and internal repair mechanism protects pteropod Limacina helicina from ocean acidification. Deep Sea Research II , 127: 41– 52. Google Scholar CrossRef Search ADS Questel, J., Blanco-Bercial, L., Hopcroft, R.R. & Bucklin, A. 2016. Phylogeography and connectivity of the Pseudocalanus (Copepoda: Calanoida) species complex in the eastern North Pacific and the Pacific Arctic Region. Journal of Plankton Research , 38: 610– 623. Google Scholar CrossRef Search ADS PubMed Raup, D.M. 1962. Computer as aid in describing form in gastropod shells. Science , 138: 150– 152. Google Scholar CrossRef Search ADS PubMed Raup, D.M. 1966. Geometric analysis of shell coiling: general problems. Journal of Paleontology , 40: 1178– 1190. Raup, D.M. 1967. Geometric analysis of shell coiling: coiling in ammonoids. Journal of Paleontology , 41: 43– 65. Raup, D.M. & Michelson, A. 1965. Theoretical morphology of the coiled shell. Science , 147: 1294– 1295. Google Scholar CrossRef Search ADS PubMed Riddle, B.R., Dawson, M.N., Hadly, E.A., Hafner, D.J., Hickerson, M.J., Mantooth, S.J. & Yoder, A.D. 2008. The role of molecular genetics in sculpting the future of integrative biogeography. Progress in Physical Geography , 32: 173– 202. Google Scholar CrossRef Search ADS Rockström, J., Steffen, W., Noone, K., Persson, Å., Chapin, F. S., Lambin, E. F., Lenton, TM, Scheffer, M, Folke, C, Schellnhuber, HJ, Nykvist, B, de Wit, CA, Hughes, T, van der Leeuw, S, Rodhe, H, Sörlin, S, Snyder, PK, Costanza, R, Svedin, U, Falkenmark, M, Karlberg, L, Corell, RW, Fabry, VJ, Hansen, J, Walker, B, Liverman, D, Richardson, K, Crutzen, P & Foley, JA. 2009. A safe operating space for humanity. Nature , 461: 472– 475. Google Scholar CrossRef Search ADS PubMed Schneider, C.A., Rasband, W.S. & Eliceiri, K.W. 2012. NIH Image to ImageJ: 25 years of image analysis. Nature Methods , 9: 671– 675. Google Scholar CrossRef Search ADS PubMed Scoble, J. & Lowe, A.J. 2010. A case for incorporating phylogeography and landscape genetics into species distribution modelling approaches to improve climate adaptation and conservation planning. Diversity and Distributions , 16: 343– 353. Google Scholar CrossRef Search ADS Sromek, L., Lasota, R. & Wolowicz, M. 2015. Impact of glaciations on genetic diversity of pelagic mollusks: Antarctic Limacina antarctica and Arctic Limacina helicina. Marine Ecology Progress Series , 525: 143– 152. Google Scholar CrossRef Search ADS Suyama, Y. & Matsuki, Y. 2015. MIG-seq: an effective PCR-based method for genome-wide single-nucleotide polymorphism genotyping using the next-generation sequencing platform. Scientific Reports , 5: 16963. Google Scholar CrossRef Search ADS PubMed Tajima, F. 1989. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics , 123: 585– 595. Google Scholar PubMed Tamura, K. 1992. Estimation of the number of nucleotide substitutions when there are strong transition-transversion and G+C content biases. Molecular Biology and Evolution , 9: 678– 687. Google Scholar PubMed Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M. & Kumar, S. 2011. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Molecular Biology and Evolution , 28: 2731– 2739. Google Scholar CrossRef Search ADS PubMed Van der Spoel, S 1967. Euthecosomata: a group with remarkable developmental stages (Gastropoda, Pteropoda) . J. Noorduijn en Zoon N. V., Gorinchem, The Netherlands. Wakita, M., Watanabe, S., Honda, M., Nagano, A., Kimoto, K., Matsumoto, K., Kitamura, M., Sasaki, K., Kawakami, H., Fujiki, T., Sasaoka, K., Nakano, Y. & Murata, A. 2013. Ocean acidification from 1997 to 2011 in the subarctic western North Pacific Ocean. Biogeosciences , 10: 7817– 7827. Google Scholar CrossRef Search ADS Wakita, M., Nagano, A., Fujiki, T. & Watanabe, S. 2017. Slow acidification of the winter mixed layer in the subarctic western North Pacific. Journal of Geophysical Research , 122 DOI:10.1002/2017JC013002. Wares, J. P. & Cunningham, C. W. 2001. Phylogeography and historical ecology of the North Atlantic intertidal. Evolution , 55: 2455– 2469. Google Scholar CrossRef Search ADS PubMed Wilson, N.G., Hunter, R.L., Lockhart, S.J. & Halanych, K.M. 2007. Multiple lineages and absence of panmixia in the ‘circumpolar’ crinoid Promachocrinus kerguelensis from the Atlantic sector of Antarctica. Marine Biology , 152: 895– 904. Google Scholar CrossRef Search ADS © The Author 2017. Published by Oxford University Press on behalf of The Malacological Society of London, all rights reserved. For Permissions, please email: email@example.com
Journal of Molluscan Studies – Oxford University Press
Published: Feb 1, 2018
It’s your single place to instantly
discover and read the research
that matters to you.
Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.
All for just $49/month
Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly
Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.
Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.
Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.
All the latest content is available, no embargo periods.
“Hi guys, I cannot tell you how much I love this resource. Incredible. I really believe you've hit the nail on the head with this site in regards to solving the research-purchase issue.”Daniel C.
“Whoa! It’s like Spotify but for academic articles.”@Phil_Robichaud
“I must say, @deepdyve is a fabulous solution to the independent researcher's problem of #access to #information.”@deepthiw
“My last article couldn't be possible without the platform @deepdyve that makes journal papers cheaper.”@JoseServera