[フレーム]
BioOne.org will be down briefly for maintenance on 19 November 2025 between 18:00-21:00 Pacific Time US. We apologize for any inconvenience.
Forgot your password?
Please wait...

We can help you reset your password using the email address linked to your BioOne Complete account.

For additional assistance, contact our Help Desk, helpdesk@bioone.org

Registered users receive a variety of benefits including the ability to customize email alerts, create favorite journals list, and save searches. Please note that a BioOne web account does not automatically grant access to full-text content. An institutional or society member subscription is required to view non-Open Access content. Contact helpdesk@bioone.org with any questions.
All Fields are Required
*
*
*
*
Password Requirements: Minimum 8 characters, must include as least one uppercase, one lowercase letter, and one number or permitted symbol Valid Symbols for password:
~ Tilde
! Exclamation Mark
@ At sign
$ Dollar sign
^ Caret
( Opening Parenthesis
) Closing Parenthesis
_ Underscore
. Period
*
Please wait...
Web Account created successfully
BioOne Complete Homepage
Advanced Search
Home > Journals > Zoological Science > Volume 37 > Issue 4 > Article
How to translate text using browser tools
2 July 2020 Inbreeding between Deep-Sea Snailfishes Careproctus pellucidus and Careproctus rastrinus in the Northwestern Pacific Ocean
Haruka Yashiki, Shogo Takami, Fabienne Spahn, Kay Sakuma, Hajime Itoh, Tomonori Hamatsu, Youji Narimatsu, Takashi Yanagimoto, Shigeaki Kojima
Author Affiliations +
Haruka Yashiki,1,2,* Shogo Takami,1,2,* Fabienne Spahn,3 Kay Sakuma,4 Hajime Itoh,5 Tomonori Hamatsu,6 Youji Narimatsu,7 Takashi Yanagimoto,8 Shigeaki Kojima1,2,**

1Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8563, Japan
2Atmosphere and Ocean Research Institute, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8564, Japan
3Swiss Federal Institute of Technology in Zurich, Rämistrasse 101, 8092 Zurich, Switzerland
4Japan Sea National Research Institute, Japan Fisheries Research and Education Agency, Niigata 951-8121, Japan
5Center for Regional Environmental Research, National Institute for Environmental Studies, 16-2 Onogawa, Tsukuba, Ibaraki 305-8506, Japan
6Kushiro Laboratory, Hokkaido National Fisheries Research Institute, Japan Fisheries Research and Education Agency, 116 Katsurakoi, Kushiro, Hokkaido 085-0802, Japan
7Tohoku National Research Institute, Japan Fisheries Research and Education Agency, 3-27-5 Shinhama, Shiogama, Miyagi 985-0001, Japan
8National Research Institute of Fisheries Science, Japan Fisheries Research and Education Agency, Yokohama, Kanagawa 236-8648, Japan

* These authors contributed equally to the paper.
** Corresponding author. E-mail: kojima@aori.u-tokyo.ac.jp
PERSONAL SIGN IN
Full access may be available with your subscription
Forgot your password?
PURCHASE THIS CONTENT
PURCHASE SINGLE ARTICLE
Price: 10ドル.00 ADD TO CART
Includes PDF & HTML, when available
This article is only available to subscribers.
It is not available for individual sale.
Access to the requested content is limited to institutions that have purchased or subscribe to this BioOne eBook Collection. You are receiving this notice because your organization may not have this eBook access.*
*Shibboleth/Open Athens users-please sign in to access your institution's subscriptions.
Additional information about institution subscriptions can be found here
This will count as one of your downloads.
You will have access to both the presentation and article (if available).
This content is available for download via your institution's subscription. To access this item, please sign in to your personal account.
Forgot your password?
No BioOne Digital Library account? Create an account
My Library
You currently do not have any folders to save your paper to! Create a new folder below.
Abstract

Genetic deviation between two deep-sea snailfishes, Careproctus pellucidus from the Pacific Ocean and Careproctus rastrinus from the Okhotsk Sea, of the Careproctus rastrinus species complex was analyzed, based on the nucleotide sequences of the mitochondrial cytochrome b (Cytb) gene. Our sampling revealed the occurrence of individuals with C. rastrinus-type Cytb genes off northeastern Japan, in the northwestern Pacific. Most of these individuals were collected from an area off Miyagi Prefecture, while few individuals were collected from areas to its north (off Iwate Prefecture) and south (off Fukushima and Ibaraki Prefectures). Phylogeographic analyses based on nucleotide sequences of the first intron region of the nuclear S7 ribosomal protein gene (S7) and 11 microsatellite loci indicated little genetic deviation between individuals with C. pellucidus-type and C. rastrinus-type Cytb genes in the area off Miyagi Prefecture. Significant genetic differences between the Pacific Ocean and the Okhotsk Sea populations may be attributed to unidirectional migration due to a greater current from the latter to the former. In addition, peritoneum and stomach colors of 10 and 17 specimens collected from the Pacific Ocean and the Okhotsk Sea, respectively, were evaluated quantitatively. The colors were significantly different for the two populations; however, they appeared to overlap. Although only one Pacific individual with a C. rastrinus-type Cytb gene was available for color evaluation, it was more similar to Pacific individuals than to Okhotsk Sea individuals. These results suggest inbreeding between C. pellucidus and C. rastrinus in the Pacific Ocean, off Miyagi Prefecture.

INTRODUCTION

Careproctus is one of the most species-rich genera in family Liparidae (Orr et al., 2015, 2019), and includes more than 100 species, more than 50 of which are found in the North Pacific (Chernova et al., 2004). Of these, species that are morphologically similar to the salmon snailfish Careproctus rastrinus Gilbert & Burke, 1912 collectively constitute the Careproctus rastrinus species complex. Kai et al. (2011) conducted the first genetic study of this species complex, and molecular phylogenetic analyses based on the 16 S ribosomal RNA and cytochrome b (Cytb) mitochondrial genes revealed that this species complex consisted of nine monophyletic lineages, including one from the Pacific Ocean, one from the Arctic Ocean and Gulf of Alaska, three from the Bering Sea, one from the Okhotsk Sea, and three from the Japan Sea.

Orr et al. (2015) revised the species in the Careproctus rastrinus complex based on their morphological characteristics, and concluded that the complex contained eight species. Of these, the pellucid snailfish Careproctus pellucidus Gilbert & Burke, 1912, from the Cytb lineage PAC1 described by Kai et al. (2011), and C. rastrinus, from the lineage OKH1, were found to be the sole and endemic species of the Pacific Ocean and the Okhotsk Sea, respectively. The diagnostic morphological characteristics distinguishing these two species, as described by Orr et al. (2015), are summarized in Table 1. However, no objective standard is available for color characteristics, and the quantitative characteristics overlap between the two species.

Table 1.

Diagnostic morphological characteristics for discriminating C. pellucidus and C. rastrinus.

img-ApY1_01.gif

Table 2.

List of sampling sites.

img-z2-32_01.gif

Following Kai et al. (2011), Hashiguchi (unpublished data) further identified individuals with the C. rastrinus-type mitochondrial DNA (mtDNA), from the Pacific Ocean, off the Sanriku Coast, in northeastern Japan, based on the nucleotide sequences of the Cytb gene. Therefore, the phylogeographic relationships between C. pellucidus and C. rastrinus need to be confirmed. In this study, the genetic relationships between these two species were analyzed, based on the nucleotide sequences of the mitochondrial Cytb gene, the first intron region of the nuclear S7 ribosomal protein gene (S7), and 11 microsatellite loci, from a total of 156 specimens. In addition, their color characteristics were evaluated quantitatively. The results of this study provide important information for revising the Careproctus rastrinus species complex.

MATERIALS AND METHODS

A total of 100 Careproctus specimens were collected from the Pacific Ocean and the Okhotsk Sea (Table 2). The specimens were preserved either by freezing or in ethanol, at the Atmosphere and Ocean Research Institute, The University of Tokyo (AORI), and the Hokkaido University Museum (HUMZ). Total DNA was extracted from the muscle of each preserved specimen, using a DNeasy Blood and Tissue Kit (Qiagen, Valencia, CA, USA), according to the manufacturer's protocol.

A part of the mitochondrial Cytb gene and the first intron region of the nuclear S7 gene were amplified by polymerase chain reaction (PCR), using the primer sets CareCBfwd (5′-GCTTCCTGCACCCTCTAATATCTCA-3′) and CareCBrev (5′-ATGAAGCCAGCTAGGGGGAATATT-3′) for Cytb (Kai et al., 2011), and S7RPEX1F (5′-TGGCCTCTTCCTTGGC-CGTC-3′) and S7RPEX2R (5′-AACTCGTCTGGCTTTTCGCC-3′) for S7 (Chow and Hazama, 1998). The PCR amplification was performed in a final reaction volume of 20 μl containing < 10 ng of the template DNA, 2 μl of 10 ×ばつ Ex Taq buffer, 2.0 μl of dNTPs (2.5 μM), 0.125 μl of Takara Ex Taq Hot Start Version DNA polymerase (5 unit/μl; Takara Bio Inc., Shiga, Japan), 2.5 μl of 10 μM forward and reverse primers, using a Thermal Cycler (Applied Biosystems, Foster City, CA, USA). The PCR conditions were as follows: an initial denaturation at 94°C for 3 minutes (min), followed by 35 cycles of denaturation at 94°C for 1 min, annealing at 58°C for 1 min (Cytb) or 45 s (S7), and extension at 72°C for 2 min (Cytb) or 1 min (S7), with a final extension at 72°C for 10 min. The PCR products were separated electrophoretically in a 2% agarose gel at 100 V for 25 min. Subsequently the gel was stained with 0.5 μg/ml ethidium bromide for about 20 minutes, and the DNA fragments were visualized under ultraviolet light. To degrade the remaining primers and nucleotides, 1.5 μL of each PCR product was mixed with 0.3 μL of ExoSAP-IT (United States Biochemical, Cleveland, OH, USA) and incubated at 37°C for 40 min and then at 80°C for 15 min. Purified PCR products were used for the cycle sequencing reactions, carried out using the BigDye Terminator Cycle Sequencing Kit, version 3.0 (Applied Biosystems) with the same PCR primers. The nucleotide sequences were determined bidirectionally, using an ABI 3130 automated DNA sequencer (Applied Biosystems).

As the nucleotide sequences of some S7 PCR products could not be determined by direct sequencing, TA cloning was performed. The PCR products were subcloned into the pMD20-T vector of the Mighty TA cloning kit (Takara Bio Inc.), following the manufacturer's protocol. Eight clones per sample were used for direct PCR with forward and reverse primers (M13M4 and M13RV, Takara Bio Inc.). The same protocols were used for direct sequencing. No sequence variation was observed among clones obtained from individual specimens.

The obtained sequences were assembled using the MEGA X version 10.0.5 software package (Kumar et al., 2018). The nucleotide sequences obtained in the present study were deposited in the DNA Data Bank of Japan (DDBJ), European Molecular Biology Laboratory (EMBL), and National Center for Biotechnology Information (NCBI) GenBank databases under the accession numbers LC534911–534925 (Cytb) and LC534926–534940 (S7).

To clarify the phylogenetic relationships among the haplotypes, a haplotype network was constructed using the TCS method in PopART software (Leigh and Bryant, 2014). The Cytb sequences of C. pellucidus (PAC1) and C. rastrinus (OKH1), determined by Kai et al. (2011), were retrieved from the DDBJ/Genbank/EMBL databases.

The unbiased fixation index, FST (Weir and Cockerham, 1984), was estimated, and its significance was tested using a nonparametric permutation approach (10,000 permutations) performed in Arlequin version 3.5.1.2. (Excoffier and Lischer, 2010). Differences in the haplotype frequencies between the populations were examined using the exact test of population differentiation (Raymond and Rousset, 1995a) in Arlequin.

Eleven microsatellite loci developed for Acanthopterygii (Orla2-91, Orla3-65, Orla5-131, Orla4-222, Orla8-113, Orla9-204, Orla18-49, Orla20-134, Orla21-231, Orla22-135, and Orla23-61; Gotoh et al., 2013) were used in the microsatellite analyses. Three multiplex PCR panels were designed using Multiplex Manager 1.0 (Holleley and Geerts, 2009); universal probe target sequence (Blacket et al., 2012) was added to the 5′-end of the forward primer and the pigtail sequence (5′-GTTTCTT; Brownstein et al., 1996) was added to the reverse primer, to reduce noise peaks during the fragment analysis. The PCR amplification was performed using a Thermal Cycler (Applied Biosystems). Each 5 μL of the multiplex PCR mixture contained 0.005 μL of each universal probe (100 μM), 0.0004 μL of forward primer (50 μM), 0.0008 μL of reverse primer (50 μM), 2 μL of 2 ×ばつ QIAGEN Multiplex PCR Master Mix (Multiplex PCR Kit, QIAGEN), and 1 μL of DNA. The following conditions were used for PCR: initial denaturation at 95°C for 15 min, followed by 35 cycles of 94°C for 30 s, 53°C for 30 s, and 65°C for 30 s, with a final extension at 60°C for 30 min. Fragments were analyzed by electrophoresis, using an ABI3130 automated DNA sequencer (Applied Biosystems), with GeneScan 500 Liz (Applied Biosystems) as a size marker. Allele size was determined using Peak Scanner Software version 2.0 (Applied Biosystems). Deviation from the Hardy-Weinberg equilibrium (HWE) for each locus and linkage disequilibrium were tested using Genepop version 4.7 (Raymond and Rousset, 1995b; Rousset, 2008), using an exact probability test (Markov chain parameters: 10,000 dememorisations, 100 batches, and 1000 iterations per batch), with the sequential Bonferroni's correction (Rice, 1989).

A Bayesian clustering approach was adopted in Structure version 2.3.4 (Pritchard et al., 2000). The burn-in period was set to 100,000 iterations, following which a Markov chain Monte Carlo (MCMC) program was run 1,000,000 times. The number of clusters (K) was set to 1–7, and 10 independent runs were conducted for each K value. The K value showing the highest ΔK was determined to be the most suitable (Evanno et al., 2005), using the online program Structure Harvester version 0.6.94 (Earl and vonHoldt, 2012). In the run with the most suitable K, the result showing the highest posterior probability (Ln P(D)) was used. The results of the independent runs for each K value were merged using a greedy algorithm in Clumpp 1.1.2 (Jakobsson and Rosenberg, 2007), and a bar plot for each K value was then generated using Distruct (Rosenberg, 2004).

A principal coordinate analysis (PCoA) was conducted based on genetic distances between individuals (DA) (Nei et al., 1983), using Populations version 1.2.32 (Langella, 1999) and GenAlEx version 6.5 (Peakall and Smouse, 2012). The genotypes were clustered by discriminant analysis of principal components (DAPC; Jombart et al., 2010) executed using the program adegenet version 2.1.0 (Jombart, 2008) for R version 3.4.4 software (R Core Team, 2016).

To evaluate the genetic differentiation between the mtDNA lineages, pairwise FST values were calculated based on microsatellites. The significance of the values was tested using a nonparametric permutation approach (10,000 permutations) in Arlequin.

Peritoneum and stomach colors of specimens of sufficiently sized specimens collected in 2019 were evaluated quantitatively based on three indices in the CIE L*a*b* color space, i.e., lightness (L*), red/green value (a*), and yellow/blue value (b*), which were measured using a CR-300 color difference meter (Konica Minolta Inc., Tokyo, Japan). Differences in the average values were tested, using the t-test if the F-test showed no significant difference in variance or using the Mann-Whitney U test if the F-test showed a significant difference in variance. A principal component analysis (PCA) was conducted based on peritoneum and stomach colors. All of these analyses were conducted using BellCurve for Excel version 2.02 (Social Survey Research Information Co., Tokyo, Japan).

RESULTS

From the 100 Careproctus specimens, 692 bp of the Cytb sequences were determined, resulting in 15 haplotypes with 23 nucleotide substitutions and no indels (insertion/deletion) ( Supplementary Table S1 (zs200011.s1.pdf)). Based on the haplotype network (Fig. 1), the specimens could be classified into the two lineages described by Kai et al. (2011), namely, PAC1 (n = 50) and OKH1 (n = 50). Specimens of the OKH1 lineage were obtained from both the Pacific Ocean (n = 22) and the Okhotsk Sea (n = 28), and are hereafter referred to as OKH1pac and OKH1okh, respectively. Only a single haplotype was obtained from the OKH1pac individuals, and it was the most dominant haplotype in the OKH1okh individuals. The test of FST value (FST = 0.01990) and the exact test based on Cytb did not indicate any significant genetic difference between OKH1pac and OKH1okh individuals (P > 0.05). Half (50.0%) of the individuals collected at stations off Miyagi Prefecture (Stations P1–5) were OKH1pac individuals. However, only a single OKH1pac individual (9.1%) was collected at stations off Iwate Prefecture (Stations 8–9), which were to the north of the Miyagi Prefecture stations, and no OKH1pac individuals were collected off Fukushima Prefecture (Stations 6–7) or Ibaraki Prefecture (Station 10), which were to the south of Miyagi Prefecture ( Supplementary Table S1 (zs200011.s1.pdf)).

Fig. 1.

A statistical parsimony haplotype network of Careproctus pellucidus and Careproctus rastrinus based on nucleotide sequences of the cytochrome b (Cytb) mitochondrial genes. Haplotype names are the same as in Supplementary Table S1 (zs200011.s1.pdf). * and ** denote the same sequences as PAC1 and OKH1 individuals in Kai et al. (2011), respectively. The areas of the circles are proportional to the frequency of the occurrence of the haplotypes. Black and white sections denote the relative occurrence frequency of individuals from the Pacific Ocean and the Okhotsk Sea, respectively. Small closed circles are intermediate haplotypes that were not observed in the present study.

img-z4-1_01.jpg

Fig. 2.

A statistical parsimony haplotype network of Careproctus pellucidus and Careproctus rastrinus based on nucleotide sequences of the first intron region of the nuclear S7 ribosomal protein gene. Haplotype names are the same as in Supplementary Table S2 (zs200011.s1.pdf). The areas of the open circles are proportional to the frequency of the occurrence of the haplotypes. Black, gray, and white sections denote the relative occurrence frequency of PAC1, OKH1pac, and OKH1okh individuals, respectively.

img-z4-3_01.jpg

Table 3.

Pairwise FST values based on the nucleotide sequences of the first intron region of the nuclear S7 ribosomal protein gene. The value indicated in parentheses is based on the individuals collected off Miyagi Prefecture only. Significant differences between groups are indicated by * (P < 0.05) and ** (P < 0.001).

img-z4-6_01.gif

Fig. 3.

Influence of the number of clusters (K) that best fit the data. ΔK is the standardized second order rate of change of the posterior probability of each K.

img-z4-8_01.jpg

Fig. 4.

Bayesian clustering of 84 Careproctus individuals from the Pacific Ocean and the Okhotsk Sea into two genetically distinguishable clusters.

img-z4-10_01.jpg

The S7 nucleotide sequences (422 bp) were determined for 49 specimens collected in 2006, 2011, and 2012. They contained 15 haplotypes with nine nucleotide substitutions and no indels ( Supplementary Table S2 (zs200011.s1.pdf)). The phylogenetic relationships between the S7 haplotypes are shown as a haplotype network in Fig. 2. Haplotypes of the two mitochondrial lineages, i.e., PAC1 and OKH1, did not form exclusive clusters in the network, and five haplotypes were shared between the lineages. The tests of pairwise FST values based on the S7 sequences indicated significant genetic differences between OKH1pac and OKHokh individuals (P < 0.01; Table 3). Genetic differences were also detected between Pacific Ocean individuals with PAC1 and OKH1pac mtDNA (P < 0.05), although the FST value was considerably lower (Table 3). The exact test revealed significant genetic differentiation between OKH1pac and OKHokh individuals (P < 0.01), but not between Pacific Ocean individuals with PAC1 and OKH1pac mtDNA (P > 0.05). In the area off Miyagi Prefecture (Stations P1–5), no significant genetic difference was observed between PAC1 and OKH1pac individuals, as evident from the results of FST value (P > 0.05; Table 3) and exact (P > 0.05) tests.

Fig. 5.

Principal coordinate analysis (PCoA) plot based on 11 microsatellites from 84 Careproctus individuals from the Pacific Ocean and the Okhotsk Sea. Open circles, closed circles, and open triangles denote PAC1, OKH1pac, and OKH1okh individuals, respectively.

img-z5-1_01.jpg

Fig. 6.

Results of the discriminant analysis of principal components (DAPC) showing genetic similarities between PAC1, OKH1pac, and OKH1okh individuals.

img-z5-8_01.jpg

A total of 84 specimens were used for microsatellite analyses. The numbers of alleles per locus and the observed and expected heterozygosities are summarized in Supplementary Table S3 (zs200011.s1.pdf). Significant deviations from HWE after the sequential Bonferroni's correction (P < 0.05) were observed in Orla9-204 for PAC1, Orla18-49 and Orla9-204 for OKH1pac, and Orla3-65 and Orla23-61 for OKH1okh. Significant linkage disequilibrium was not found for any pairs. A Bayesian clustering analysis based on 11 microsatellites was implemented to estimate the genetic population structure of Careproctus individuals from the Pacific Ocean and the Okhotsk Sea. The most probable number of clusters (K) was estimated to be two (Fig. 3). In the Bayesian clustering of individuals, OKH1pac individuals exhibited genetic characteristics intermediate to PAC1 and OKH1okh individuals (Fig. 4). PCoA also revealed a trend similar to that observed in the Bayesian clustering (Fig. 5). The first and second principal coordinates explained 9.14% and 6.20% of the total variation, respectively. In the scatter plot of the first and second scores, the plotted areas of PAC1 and OKH1okh individuals did not overlap with each other; however, plotted areas of OKH1pac individuals overlapped with those of both the aforementioned individuals. The DAPC analysis showed that OKH1pac individuals were more closely related to PAC1 individuals than to OKHokh individuals (Fig. 6). Pairwise FST values based on microsatellites indicated significant genetic differentiation between all pairs among the three groups (P < 0.01); however, the FST value between PAC1 and OKH1pac individuals was slightly lower (Table 4). In the area off Miyagi Prefecture (Sts. P1–5), no significant genetic difference could be detected between the PAC1 and OKH1pac individuals by the tests of FST value (P > 0.05; Table 4) or the exact test (P > 0.05).

Fig. 7.

Peritoneum (A) and stomach (B) colors of Careproctus individuals from the Pacific Ocean and the Okhotsk Sea, expressed as lightness (L*), red/green value (a*), and yellow/blue value (b*). Open circles, open triangles, and closed circles denote PAC1, OKH1pac, and OKH1okh individuals, respectively.

img-z6-1_01.jpg

Table 4.

Pairwise FST values based on 11 microsatellites. The value indicated in parentheses is based on individuals collected off Miyagi Prefecture only. Significant differences between groups are indicated by ** (P < 0.01).

img-A0r6_01.gif

Peritoneum and stomach colors of 9 PAC1, 1 OKH1pac, and 17 OKHokh individuals were quantitatively evaluated (Fig. 7). Significant differences were observed between the individuals from the Pacific Ocean and the Okhotsk Sea in terms of the stomach L* values (Mann-Whitney U test, P < 0.001), and the peritoneum a* (t-test, P < 0.05) and b* (t-test, P < 0.001) values. As described by Orr et al. (2015), lightness (L*) of the stomach of individuals from the Okhotsk Sea was significantly higher than that of individuals from the Pacific Ocean, although the two ranges overlapped. The OKH1pac individual exhibited greater color similarity to most PAC1 individuals than to most OKHokh individuals (Fig. 7).

DISCUSSION

Our results showed that the OKH1 mitochondrial lineage is distributed not only in the Okhotsk Sea, but also in the Pacific Ocean (Fig. 1), which is different from the findings of Kai et al. (2011). It is worth noting that 22 of the 72 Pacific individuals sampled in the present study were classified as OKH1 individuals. These individuals were collected almost exclusively from the area off Miyagi Prefecture, while few of these individuals were collected from areas north (off Iwate Prefecture) and south (off Fukushima and Ibaraki Prefectures) of Miyagi Prefecture. Kai et al. (2011) only used C. pellucidus specimens collected from Iwate and Fukushima Prefectures, which might explain why they did not detect individuals belonging to the OKH1 mitochondrial lineage from the Pacific Ocean.

In contrast to the OKH1 mitochondrial lineage, individuals with the PAC1-type mtDNA have not yet been obtained from the Okhotsk Sea. According to Orr et al. (2015), C. pellucidus and C. rastrinus have been collected from depths of 145–300 and 114–217 m, respectively. As the Okhotsk Sea and the Pacific Ocean are connected by deep straits, namely, the Kruzenshterna Strait (1900 m deep) and the Bussol Strait (2300 m deep), Careproctus habitats in these areas would not have been isolated from each other, even during the Last Glacial Maximum (LGM). In the Bussol Strait, water flow from the Okhotsk Sea to the Pacific Ocean is in excess of the opposing flow, and most of the water mass flowing into the Okhotsk Sea is estimated to flow back to the Pacific Ocean (Katsumata et al., 2004). Because of the unidirectional migration occurring due to the stronger current flowing from the Okhotsk Sea to the Pacific Ocean, individuals with the OKH1-type Cyt b genes could have migrated to the Pacific Ocean, while those with the PAC1-type Cyt b genes would have found it difficult to migrate to the Okhotsk Sea. Unidirectional migration from the Okhotsk Sea to the Pacific Ocean has also been suggested for the Japan Sea eelpout Bothrocara hollandi (Jordan and Hubbs, 1925), based on the contrasting genetic diversity between the populations of the two areas (Kojima et al., 2011). However, Tohkairin et al. (2016) stated that the population of the barred snailfish Crystallichthys matsushimae (Jordan and Snyder, 1902) in the Okhotsk Sea was most likely formed by immigrants from both the Pacific Ocean and the Japan Sea. As suggested by Tohkairin et al. (2016), this difference in the migratory history may be attributed to a shallower distributional depth of C. matsushimae compared to that of the Careproctus rastrinus complex and B. hollandi.

Although significant genetic differentiation between PAC1 and OKH1pac individuals was detected for the nuclear S7 gene and microsatellites, this differentiation was not significant when the analyses were limited to the individuals off Miyagi Prefecture, where OKH1pac individuals were almost exclusively collected (Tables 3 and 4). The DAPC analysis showed that OKH1pac individuals were more closely related to PAC1 individuals than to OKHokh individuals (Fig. 6). Unfortunately, only a single OKH1pac individual was available for color evaluation, and its peritoneum and stomach were more similar to most PAC1 individuals than to most OKHokh individuals (Fig. 7). If complete reproductive isolation had been established between C. pellucidus and C. rastrinus, OKH1pac individuals would exhibit characteristics that are more similar to OKH1okh than to PAC1 individuals. Thus, the present results suggest that C. pellucidus and C. rastrinus are inbreeding in the Pacific Ocean off Miyagi Prefecture.

Another possibility is that the OKH1-type mtDNAs within the Pacific population might have originated from past hybridization between C. pellucidus and C. rastrinus, resulting in the introgression of mtDNA of the latter into the former. The limited geographical distribution of OKH1pac individuals and the assumed dispersal direction previously discussed suggest that hybridization occurred between migrating C. rastrinus and indigenous C. pellucidus within the Pacific Ocean. As organelle genes are often preferentially introgressed during interspecific hybridization, and as introgression almost exclusively occurs from the local to the invading species, irrespective of the relative densities of the two species (Currat et al., 2008), the possibility of an introgression event from C. rastrinus to C. pellucidus seems to be very low.

In the present study, OKH1pac individuals were collected almost exclusively from the area off Miyagi Prefecture. Such a geographic trend might reflect an interspecific difference in the preference of physical environments. However, it can also be attributed to annual variation, as sampling off Iwate and Ibaraki Prefectures was carried out 8 years after that off Miyagi and Fukushima Prefectures. Systematic sampling of additional specimens, especially from areas around deep straits connecting the Pacific Ocean and the Okhotsk Sea, as well as more comprehensive molecular and morphological analyses are necessary in order to reveal the evolutionary relationships between C. pellucidus and C. rastrinus.

ACKNOWLEDGMENTS

We are grateful to Mr. H. Hashiguchi, Graduate School of Frontier Sciences, The University of Tokyo, for providing valuable information on the Cytb sequences of the Pacific specimens, and to Dr. T. Takano, Meguro Parasitological Museum, Tokyo, for useful advice. Dr. T. Yahagi kindly helped with statistical analyses. We thank the captains, officers, and crew members of the research vessels (R/V) Wakataka-Maru of Tohoku National Fishery Institute, Kaiyo-Maru No. 5 of Kaiyo Engineering Co., Ltd., and Shinsei-Maru of the Atmosphere and Ocean Research Institute, The University of Tokyo, for their help with sampling. This study was supported in part by the Fisheries Agency of Japan under the ‘Assessment of Fisheries Stocks in the Waters around Japan’ project. Dr. F. Tashiro, Hokkaido University Museum, Hokkaido University, kindly provided tissue samples preserved in the museum. We are also thankful to Mr. K. Kofuji, Ibaraki Prefectural Oarai Aquarium, and Mr. K. Matsuzaki, Aquamarine Fukushima, for providing specimens collected in 2019 from the Pacific Ocean and the Okhotsk Sea, respectively. Dr. J. Ashi, Graduate School of Frontier Sciences, The University of Tokyo, kindly helped with color evaluations. This research was partly supported by the Tohoku Ecosystem-Associated Marine Sciences (TEAMS) of the Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT).

COMPETING INTERESTS

The authors have no competing interests to declare.

AUTHOR CONTRIBUTIONS

SK designed the study. KS, TH, YN, and SY collected the specimens. HY, ST, FS, and HI analyzed the nucleotide sequences and microsatellites. ST measured the color differences. SK and HY prepared the manuscript. All authors have read and approved the final manuscript.

SUPPLEMENTARY MATERIALS

Supplementary materials for this article are available online. (URL: https://doi.org/10.2108/zs200011)

Supplementary Figure S1 (zs200011.s1.pdf). Bayesian clustering of 84 Careproctus individuals from the Pacific Ocean and the Okhotsk Sea into two, three, four, five, and six genetically distinguishable clusters.

Supplementary Figure S2 (zs200011.s1.pdf). Results of a principal component analysis (PCA) based on peritoneum colors, stomach colors, and both.

Supplementary Table S1 (zs200011.s1.pdf). Haplotype composition of C. pellucidus and C. rastrinus for the cytochrome b (Cytb) mitochondrial genes.

Supplementary Table S2 (zs200011.s1.pdf). Haplotype composition of C. pellucidus and C.rastrinus for the first intron region of the nuclear S7 ribosomal protein gene.

Supplementary Table S3 (zs200011.s1.pdf). Numbers of alleles per locus and expected and observed heterozygosities for Careproctus individuals from the Pacific Ocean and the Okhotsk Sea.

REFERENCES

1.

Blacket MJ, Robin C, Good RT, Lee SF, Miller AD (2012) Universal primers for fluorescent labelling of PCR fragments—an efficient and cost-effective approach to genotyping by fluorescence. Mol Ecol Resour 12: 456–463 Google Scholar

2.

Brownstein MJ, Carpten JD, Smith JR (1996) Modulation of non-templated nucleotide addition by Taq DNA polymerase: primer modifications that facilitate genotyping. Biotechniques 20: 1004–1010 Google Scholar

3.

Chernova NV (2004) Family Liparidae Scopoli 1777, snailfishes. Calif Acad Sci Annot Checklists Fishes 31: 1–72 Google Scholar

4.

Chow S, Hazama K (1998) Universal PCR primers for S7 ribosomal protein gene introns in fish. Mol Ecol 7: 1255–1256 Google Scholar

5.

Currat M, Ruedi M, Petit RJ, Excoffier L (2008) The hidden side of invasions: massive introgression by local genes. Evolution 62: 1908–1920 Google Scholar

6.

Earl DA, vonHoldt BM (2012) Structure Harvester: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour 4: 359–361 Google Scholar

7.

Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 14: 2611–2620 Google Scholar

8.

Excoffier L, Lischer HE (2010) Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10: 564–567 Google Scholar

9.

Gotoh RO, Tamate S, Yokoyama J, Tamate HB, Hanzawa N (2013) Characterization of comparative genome-derived simple sequence repeats for acanthopterygian fishes. Mol Ecol Resour 13: 461–472 Google Scholar

10.

Holleley CE, Geerts PG (2009) Multiplex Manager 1.0: a cross-platform computer program that plans and optimizes multiplex PCR. BioTechniques 46: 511–517 Google Scholar

11.

Jakobsson M, Rosenberg NA (2007) CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics 23: 1801–1806 Google Scholar

12.

Jombart T (2008) adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24: 1403–1405 Google Scholar

13.

Jombart T, Devillard S, Balloux F (2010) Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet 11: 94 Google Scholar

14.

Kai Y, Orr JW, Sakai K, Nakabo T (2011) Genetic and morphological evidence for cryptic diversity in the Careproctus rastrinus species complex (Liparidae) of the North Pacific. Ichthyol Res 58: 143–154 Google Scholar

15.

Katsumata K, Ohshima KI, Kono T, Itoh M, Yasuda I, Volkov YN, et al. (2004) Water exchange and tidal currents through the Bussol' Strait revealed by direct current measurements. Jour Geophys Res 109: C09S06 Google Scholar

16.

Kojima S, Maeda R, Sakuma K, Kokubu Y, Hagihara S, Itoh M (2011) Genetic characterization of the northwestern Pacific population of a deep-sea demersal fish, Bothrocara hollandi. Plankton Benthos Res 6: 108–114 Google Scholar

17.

Kumar S, Stecher G, Li M, Knyaz C, Tamura K (2018) MEGA X: Molecular Evolutionary Genetics Analysis across computing platforms. Mol Biol Evol 35: 1547–1549 Google Scholar

18.

Langella O (1999) Populations, Ver. 1.2.31: a population genetic software. CNRS UPR9034. http://bioinformatics.org/populations/. Accessed 19 December 2019 Google Scholar

19.

Leigh JW, Bryant D (2014) PopArt: Full feature software for population genetics. http://popart.otago.ac.nz. Accessed 19 December 2019 Google Scholar

20.

Nei M, Tajima F, Tateno Y (1983) Accuracy of estimated phylogenetic trees from molecular data. Jour Mol Evol 19: 153–170 Google Scholar

21.

Orr JW, Kai Y, Nakabo T (2015) Snailfishes of the Careproctus rastrinus complex (Liparidae): redescriptions of seven species in the North Pacific Ocean region, with the description of a new species from the Beaufort Sea. Zootaxa 4018: 301–348 Google Scholar

22.

Orr JW, Spies I, Stevenson DE, Longo GFC, Kai Y, Ghods S, et al. (2019) Molecular phylogenetics of snailfishes (Cottoidei: Liparidae) based on Mt DNA and RADseq genomic analyses, with comments on selected morphological characters. Zootaxa 4642: 1–79 Google Scholar

23.

Peakall R, Smouse PE (2012) GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics 28: 2537–2539 Google Scholar

24.

Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959 Google Scholar

25.

Raymond M, Rousset F (1995a) An exact test for population differentiation. Evolution 49: 1280–1283 Google Scholar

26.

Raymond M, Rousset F (1995b) GENEPOP (ver. 3.3); a population genetics software for exact test and ecumenicism. J Hered 86: 248–249 Google Scholar

27.

R Core Team (2016) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/. Accessed 15 December 2019 Google Scholar

28.

Rice WR (1989) Analyzing tables of statistical tests. Evolution 43: 223–225 Google Scholar

29.

Rosenberg NA (2004) DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes 4: 137–138 Google Scholar

30.

Rousset F (2008) GENEPOP'007: a complete re-implementation of the GENEPOP software for Windows and Linux. Mol Ecol Resour 8: 103–106 Google Scholar

31.

Tohkairin A, Kai Y, Ueda Y, Hamatsu T, Ito M, Nakabo T (2016) Genetic population structure of Crystallichthys matsushimae (Cottoidei: Liparidae) with comments on color variation. Ichthyol Res 63: 370–381 Google Scholar

32.

Weir BS, Cockerham CC (1984) Estimating F-statistics for the analysis of population structure. Evolution 38: 1358–1370 Google Scholar
© 2020 Zoological Society of Japan
Haruka Yashiki, Shogo Takami, Fabienne Spahn, Kay Sakuma, Hajime Itoh, Tomonori Hamatsu, Youji Narimatsu, Takashi Yanagimoto, and Shigeaki Kojima "Inbreeding between Deep-Sea Snailfishes Careproctus pellucidus and Careproctus rastrinus in the Northwestern Pacific Ocean," Zoological Science 37(4), 323-330, (2 July 2020). https://doi.org/10.2108/zs200011
Received: 24 January 2020; Accepted: 12 April 2020; Published: 2 July 2020
My Library
You currently do not have any folders to save your paper to! Create a new folder below.
KEYWORDS
Careproctus rastrinus species complex
inbreeding
intron
microsatellite
Okhotsk Sea
Pacific Ocean
RIGHTS & PERMISSIONS
Get copyright permission
Haruka Yashiki, Shogo Takami, Fabienne Spahn, Kay Sakuma, Hajime Itoh, Tomonori Hamatsu, Youji Narimatsu, Takashi Yanagimoto, Shigeaki Kojima "Inbreeding between Deep-Sea Snailfishes Careproctus pellucidus and Careproctus rastrinus in the Northwestern Pacific Ocean," Zoological Science, 37(4), 323-330, (2 July 2020)
Include:
Format:
Access provided by
Back to Top

KEYWORDS/PHRASES

Keywords
in
Remove
in
Remove
in
Remove
+ Add another field

PUBLICATION TITLE:


COLLECTION:


PUBLICATION YEARS

Range
Single Year

Clear Form

AltStyle によって変換されたページ (->オリジナル) /