Comprehensive Draft Genome Analyses of Three Rockfishes (Scorpaeniformes, Sebastiscus) via Genome Survey Sequencing.

Sebastiscus species, marine rockfishes, are of essential economic value. However, the genomic data of this genus is lacking and incomplete. Here, whole genome sequencing of all species of Sebastiscus was conducted to provide fundamental genomic information. The genome sizes were estimated to be 802.49 Mb (S. albofasciatus), 786.79 Mb (S. tertius), and 776.00 Mb (S. marmoratus) by using k-mer analyses. The draft genome sequences were initially assembled, and genome-wide microsatellite motifs were identified. The heterozygosity, repeat ratios, and numbers of microsatellite motifs all suggested possibly that S. tertius is more closely related to S. albofasciatus than S. marmoratus at the genetic level. Moreover, the complete mitochondrial genome sequences were assembled from the whole genome data and the phylogenetic analyses genetically supported the validation of Sebastiscus species. This study provides an important genome resource for further studies of Sebastiscus species.

Inbreeding between Deep-Sea Snailfishes Careproctus pellucidus and Careproctus rastrinus in the Northwestern Pacific Ocean.

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.

Phylogeographic analysis of the genus Platycephalus along the coastline of the northwestern Pacific inferred by mitochondrial DNA.

BACKGROUND: Flathead fishes of the genus Platycephalus are economically important demersal fishes that widely inhabit the continental shelves of tropical and temperate sea waters. This genus has a long history of taxonomic revision, and recently four Platycephalus species (Platycephalus sp. 1, Platycephalus sp. 2, P. indicus, and P. cultellatus) in the northwestern Pacific Ocean (NWP) have been recognized and redescribed. However, many aspects of their systematics and evolutionary history are unclear. RESULTS: A total of 411 individuals were sampled from 22 different sites across their distributions in the NWP. Three mitochondrial loci were sequenced to clarify the phylogeny and phylogeographic history of the fishes. The results showed significant differentiation of four Platycephalus species in the NWP with well-supported clades in which Platycephalus sp. 1 and Platycephalus sp. 2 were the closest, clustered with P. cultellatus, while their genetic relationship with P. indicus was the furthest. There were significant genealogical branches corresponding to P. indicus but not to other Platycephalus lineages. We further examined the phylogeographic patterns of 16 Platycephalus sp. 1 populations along the coastlines of China and Japan. A total of 69 haplotypes were obtained, with 23 shared among populations. One dominant haplotypic group, with a modest lineage structure and low levels of haplotype diversity and nucleotide diversity, was observed among Platycephalus sp. 1 populations. The demographic history reconstruction suggested a Platycephalus sp. 1 population expansion event dating back to the late Pleistocene. CONCLUSIONS: Distributional rang variations may be the crucial factors shaping the genetic relationships of the genus Platycephalus. Reproductive schooling and potential egg/larval dispersal ability, coupled with the effects of ocean currents, are responsible for the present phylogeographic pattern of Platycephalus sp. 1.

Restricted Gene Flow for Gadus macrocephalus from Yellow Sea Based on Microsatellite Markers: Geographic Block of Tsushima Current.

The Pacific cod Gadus macrocephalus is a demersal, economically important fish in the family Gadidae. Population genetic differentiation of Pacific cod was examined across its northwestern Pacific range by screening variation of eight microsatellite loci in the present study. All four populations exhibited high genetic diversity. Pairwise fixation index (Fst) suggested a moderate to high level of genetic differentiation among populations. Population of the Yellow Sea (YS) showed higher genetic difference compared to the other three populations based on the results of pairwise Fst, three-dimensional factorial correspondence analysis (3D-FCA) and STRUCTURE, which implied restricted gene flow among them. Wilcoxon signed rank tests suggested no significant heterozygosity excess and no recent genetic bottleneck events were detected. Microsatellite DNA is an effective molecular marker for detecting the phylogeographic pattern of Pacific cod, and these Pacific cod populations should be three management units.

Population genetic structure of chub mackerel Scomber japonicus in the Northwestern Pacific inferred from microsatellite analysis.

Marine pelagic fishes are usually characterized by subtle but complex patterns of genetic differentiation, which are influenced by both historical process and contemporary gene flow. Genetic population differentiation of chub mackerel, Scomber japonicus, was examined across most of its range in the Northwestern Pacific by screening variation of eight microsatellite loci. Our genetic analysis detected a weak but significant genetic structure of chub mackerel, which was characterized by areas of gene flow and isolation by distance. Consistent with previous estimates of stock structure, we found genetic discontinuity between Japan and China samples. Local-scale pattern of genetic differentiation was observed between samples from the Bohai Sea and North Yellow Sea and those from the East China Sea, which we ascribed to differences in spawning time and migratory behavior. Furthermore, the observed homogeneity among collections of chub mackerel from the East and South China Seas could be the result of an interaction between biological characteristics and marine currents. The present study underlies the importance of understanding the biological significance of genetic differentiation to establish management strategies for exploited fish populations.

Limits of Bayesian skyline plot analysis of mtDNA sequences to infer historical demographies in Pacific herring (and other species).

A previous analysis of Pacific herring mitochondrial (mt) DNA with Bayesian skyline plots (BSPs) was interpreted to reflect population growth in the late Pleistocene that was preceded by population stability over several hundred thousand years. Here we use an independent set of mtDNA control region (CR) sequences and simulations to test these hypotheses. The CR haplotype genealogy shows three deeply divided lineages, A, B and C, with divergences ranging from d=1.6% to 1.9% and with similar genetic diversities (h=0.95, 0.96, 0.94; Θ(π)=0.011, 0.012, 0.014, respectively). Lineage A occurs almost exclusively in the NW Pacific and Bering Sea, but lineages B and C are co-distributed in the Northeastern Pacific. This distribution points to a historical allopatric separation between A and B-C across the North Pacific during Pleistocene glaciations. The origins of B and C are uncertain. One hypothesis invokes long-term isolation of lineage C in the Sea of Cortez, but the present-day lack of geographical segregation from lineage B argues for lineage sorting to explain the deep divergence between B and C. BSPs depict rapid population growth in each lineage, but the timing of this growth is uncertain, because of questions about an appropriate molecular clock calibration. We simulated historical demographies under a Pleistocene climate model using observed genetic parameters. BSPs for these sequences showed rapid population growth after the Last Glacial Maximum (LGM) 18-20 k years ago and a flat population history during previous climate fluctuations. Population declines during the LGM appear to have erased signals of previous population fluctuations.

Population genetic structure of the point-head flounder, Cleisthenes herzensteini, in the Northwestern Pacific.

Intraspecific phylogenies can provide useful insights into how populations have been shaped by historical and contemporary processes. To determine the population genetic structure and the demographic and colonization history of Cleisthenes herzensteini in the Northwestern Pacific, one hundred and twenty-one individuals were sampled from six localities along the coastal regions of Japan and the Yellow Sea of China. Mitochondrial DNA variation was analyzed using DNA sequence data from the 5' end of control region. High levels of haplotype diversity (>0.96) were found for all populations, indicating a high level of genetic diversity. No pattern of isolation by distance was detected among the population differentiation throughout the examined range. Analyses of molecular variance (AMOVA) and the conventional population statistic Fst revealed no significant population genetic structure among populations. According to the exact test of differentiation among populations, the null hypothesis that C. herzensteini within the examined range constituted a non-differential mtDNA gene pool was accepted. The demographic history of C. herzensteini was examined using neutrality test and mismatch distribution analyses and results indicated Pleistocene population expansion (about 94-376 kya) in the species, which was consistent with the inference result of nested clade phylogeographical analysis (NCPA) showing contiguous range expansion for C. herzensteini. The lack of phylogeographical structure for the species may reflect a recent range expansion after the glacial maximum and insufficient time to attain migration-drift equilibrium.

Molecular phylogeny of a red-snow-crab species complex using mitochondrial and nuclear DNA markers.

Nucleotide sequence variation of mitochondrial DNA COI and nuclear rRNA gene regions was used to reconstruct phylogenetic relationships for the red-snow-crab species complex, including the red snow crab, Chionoecetes japonicus, its nominal subspecies, C. japonicus pacificus, and the triangle tanner crab, C. angulatus. The topologies of the Bayesian and neighbor-joining (NJ) trees of the COI and of NJ trees of rRNA sequences placed C. japonicus and C. angulatus in a single clade. The net sequence divergence between these taxa was d(net) = 0.000 in COI, and strongly suggests that these taxa represent a single species. In contrast, haplotypes in C. j. pacificus clustered separately from the C. japonicus - C. angulatus clade. Net sequence divergence from C. japonicus - C. angulatus to C. j. pacificus was d(net) = 0.026 in COI, indicating that C. j. pacificus should be elevated to a separate species, C. pacificus. A 165 bp insert appeared in the rRNA gene of C. j. pacificus, but was absent in the remaining species of Chionoecetes. This autapomorphic condition in C. j. pacificus adds support for an independent evolution of this taxon. Evolutionary divergences between these taxa may reflect contrasting evolutionary process influenced by ocean bathymetry.

Detection of heteroplasmy and nuclear mitochondrial pseudogenes in the Japanese spiny lobster Panulirus japonicus.

Partial mtDNA cytochrome oxidase subunit I (COI) fragments and near entire stretch of 12S rDNA (12S) and control region (Dloop) of the Japanese spiny lobster (Panulirus japonicus) (n = 3) were amplified by PCR and used for direct nucleotide sequencing and for clone library-based nucleotide sequence analysis. Nucleotide sequences of a total of 75 clones in COI, 77 in 12S and 92 in Dloop were determined. Haplotypes of the clones matched with those obtained by direct sequencing were determined to be genuine mtDNA sequence of the individual. Phylogenetic analysis revealed several distinct groups of haplotypes in all three regions. Genuine mtDNA sequences were observed to form a group with their closely related variables, and most of these variables may be due to amplification error but a few to be heteroplasmy. Haplotypes determined as nuclear mitochondrial pseudogenes (NUMTs) formed distinct groups. Nucleotide sequence divergence (K2P distance) between genuine haplotypes and NUMTs were substantial (7.169-23.880% for COI, 1.336-23.434% for 12S, and 7.897-71.862% for Dloop). These values were comparable to or smaller than those between species of the genus Panulirus, indicating that integration of mtDNA into the nuclear genome is a continuous and dynamic process throughout pre- and post-speciation events. Double peaks in electropherograms obtained by direct nucleotide sequencing were attributed to common nucleotides shared by multiple NUMTs. Information on the heteroplasmy and NUMTs would be very important for addressing their impact on direct nucleotide sequencing and for quality control of nucleotide sequences obtained.

Molecular Identification of Mid to Final Stage Slipper Lobster Phyllosoma Larvae of the Genus Chelarctus (Crustacea: Decapoda: Scyllaridae) Collected in the Pacific with Descriptions of Their Larval Morphology.

Morphological descriptions of phyllosoma larvae are essential for correct species identification and investigating the spatiotemporal distribution and recruitment process of spiny and slipper lobsters. Species identification of the phyllosoma larvae in the Scyllarinae subfamily is particularly difficult because of the morphological similarities among species and the scarcity of morphological information describing correct species identity. We extracted mid-to final-stage (V to VIII) phyllosoma larvae (n = 12) belonging to the subfamily Scyllarinae from several plankton samples collected in the Pacific and then performed molecular species identification using mitochondrial DNA COI and 16S rDNA sequence analyses. Three larvae collected around the Ryukyu Archipelago were identified as Chelarctus aureus (stage VI to VIII), and four collected around the Ryukyu Archipelago and Ogasawara Islands were identified as C. virgosus (V to VIII). One larva (V) collected in the central South Pacific was determined to be a subspecies of C. crosnieri. DNA barcodes could not be made for the remaining four larvae (V to VIII) collected around the Ryukyu Archipelago (designated by ?Chelarctus sp-1). Based on the morphological characteristics of the C. virgosus phyllosoma described in this study and the adult distributions reported to date, C. cultrifer phyllosomas previously reported in Japanese and Taiwanese waters are likely to be C. virgosus. This paper also presents a set of diagnostic morphological characteristics that can be used to discriminate among these four species of Chelarctus and from other genera in the subfamily Scyllarinae.

Demographic history and population structure of blackfin flounder (Glyptocephalus stelleri) in Japan revealed by mitochondrial control region sequences.

The demographic history and population genetic structure of the blackfin flounder (Glyptocephalus stelleri) along coastal regions of Japan were investigated. Genetic variation in DNA sequences was examined from the first hypervariable region of the mitochondrial DNA control region. A high level of haplotypic diversity (h = 0.99 +/- 0.004) was detected, indicating a high level of intrapopulation genetic diversity. The starburst structure of the minimum spanning tree suggested a very recent origin for most haplotypes. The demographic history of G. stelleri was examined using neutrality tests and mismatch distribution analysis, which also indicated a Pleistocene population expansion at about 124,100-413,400 years ago. Hierarchical molecular variance analysis and conventional population Fst comparisons revealed no significant genetic differentiation throughout the range examined.

Comparative analysis of the complete mitochondrial genomes of three rockfishes (Scorpaeniformes, Sebastiscus) and insights into the phylogenetic relationships of Sebastidae.

Mitochondrial genome is a powerful molecule marker to provide information for phylogenetic relationships and revealing molecular evolution in ichthyological studies. Sebastiscus species, a marine rockfish, are of essential economic value. However, the taxonomic status and phylogenetic relationships of Sebastidae have been controversial so far. Here, the mitochondrial genomes (mitogenomes) of three species, S. tertius, S. albofasciatus, and S. marmoratus, were systemically investigated. The lengths of the mitogenomes' sequences of S. tertius, S. albofasciatus, and S. marmoratus were 16910, 17056, and 17580 bp, respectively. It contained 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNA (tRNA) genes, and one identical control region (D-loop) among the three species. The genetic distance and Ka/Ks ratio analyses indicated 13 PCGs were suffering purifying selection and the selection pressures were different from certain deep-sea fishes, which were most likely due to the difference in their living environment. The phylogenetic tree was constructed by Bayesian Inference (BI) and Maximum Likelihood (ML). Most interestingly, the results indicated that Sebastidae and Scorpaenidae were grouped into a separate branch, so the taxonomic status of Sebastidae should be classified into subfamily Sebastinae. Our results may lead to a taxonomic revision of Scorpaenoidei.

Population genetic structure of Marbled Rockfish, Sebastiscusmarmoratus (Cuvier, 1829), in the northwestern Pacific Ocean.

Sebastiscusmarmoratus is an ovoviviparous fish widely distributed in the northwestern Pacific. To examine the gene flow and test larval dispersal strategy of S.marmoratus in Chinese and Japanese coastal waters, 421 specimens were collected from 22 localities across its natural distribution. A 458 base-pair fragment of the mitochondrial DNA (mtDNA) control region was sequenced to examine genetic diversity and population structure. One-hundred-six variable sites defined 166 haplotypes. The populations of S.marmoratus showed high haplotype diversity with a range from 0.8587 to 0.9996, indicating a high level of intrapopulation genetic diversity. Low non-significant genetic differentiation was estimated among populations except those of Hyogo, Behai, and Niiigata, which showed significant genetic differences from the other populations. The demographic history examined by neutrality tests, mismatch distribution analysis, and Bayesian skyline analysis suggested a sudden population expansion dating to the late Pleistocene. Recent population expansion in the last glacial period, wide dispersal of larvae by coastal currents, and the homogeneity of the environment may have important influences on the population genetic pattern. Knowledge of genetic diversity and genetic structure will be crucial to establish appropriate fishery management of S.marmoratus.

Population genomics reveals possible genetic evidence for parallel evolution of Sebastiscus marmoratus in the northwestern Pacific Ocean.

Understanding patterns of population diversity and structuring among marine populations is of great importance for evolutionary biology, and can also directly inform fisheries management and conservation. In this study, genotyping-by-sequencing was used to assess population genetic diversity and connectivity of Sebastiscus marmoratus. Based on 130 individuals sampled from 10 locations in the northwestern Pacific Ocean, we identified and genotyped 17 653 single-nucleotide polymorphisms. The patterns of genetic diversity and population differentiation suggested that the Okinawa Trough might be the ancestral centre of S. marmoratus after the Last Glacial Maximum. A shallow genetic structure was observed among sampled populations based on the implemented structuring approaches. Surprisingly, we detected genetic homogeneity in two population pairs (i.e. Xiamen-Niigata and Zhuhai-Iki Island), in which populations have large geographical and latitudinal intervals. Population structure and allele frequency distribution based on outlier loci also mirrored the observed genetic homogeneity in the above-mentioned population pairs. Integrated with biological, environmental and genomic data, our results provide possible genetic evidence for parallel evolution. Our study also provides new perspectives on the population structure of S. marmoratus, which could have important implications for sound management and conservation of this fishery species.

Levels and patterns of genetic variation in Japanese whiting (Sillago japonica) based on mitochondrial DNA control region.

Japanese whiting (Sillago japonica) is a commercially important demersal fish distributed along the coasts of East Asia and becomes recently an aquaculture species. Despite its commercial importance, the levels and patterns of population genetic variation remain poorly understood. In this study, 346 specimens were collected from 14 localities along the coastal waters of China, Korea and Japan and their genetic variation was analyzed with mtDNA control region (D-loop) sequences. A total of 131 polymorphic sites were detected which determined 294 haplotypes. A pattern with high levels of haplotype diversity (h = 0.999 ± 0.001) and nucleotide diversity (л = 0.030 ± 0.015) was detected in the examined range. Analyses of analysis of molecular variance (AMOVA) and Fst showed that no significant genetic differentiation existed among China, Korea and Japan populations, excepting for the populations between Ise Bay (IBP) sample and the other ones. Minimum spanning tree (MST), neutrality tests, mismatch distribution and Bayesian skyline analyses indicated that the species along coastline of China, Korea and Japan have experienced population expansions originated in its most recent history at about 106-423 kya during the late Pleistocene glaciations and deglaciations periods.

Close-up observations on the spawning behavior of a captive Japanese flying squid (Todarodes pacificus).

The spawning behavior of a Japanese flying squid (Todarodes pacificus) is described based on up-close observation of a captive female. The squid was first transferred from a 10-ton tank to a polystyrene plastic box containing 45 liters of seawater. About one hour later, the mantle-contraction rate increased rapidly, followed by a brief convulsion of the mantle and arms and a whitening of the body. The mantle contractions become shallow and rapid, and several seconds later, semitransparent jelly presumably from the nidamental glands emerged from the funnel and passed between the ventral pair of arms. Approximately 90 seconds after the egg mass first emerged, the female began ejecting oocytes through the funnel into the egg mass using rapid, powerful mantle contractions. Soon after the oocytes were ejected, translucent strands (presumably sperm) emanated from the buccal membrane. The female continued to eject oocytes for approximately two minutes, after which the mantle convulsed, and the mantle-contraction rate decreased slowly for about one minute until the contractions stopped. The squid died soon afterwards.

Molecular diet analysis of Anguilliformes leptocephalus larvae collected in the western North Pacific.

Natural diets of leptocephalus larvae have been enigmatic. In this study, we collected DNA samples from the gut contents and body surface of leptocephali belonging to the five Anguilliform families (Anguillidae, Chlopsidae, Congridae, Muraenidae, and Serrivomeridae) from the northwest Pacific and performed next-generation 18S rDNA sequencing. Wide variety of eukaryotes was detected in both samples, from which eight eukaryotic groups (jellyfish, conoid parasite, tunicate, copepod, krill, segmented worm, fungi, and dinoflagellate) were selected on the basis of abundance. All groups except conoid parasites were common in both the samples. Cnidarian 18S rDNA reads were the most abundant in both the samples; however, the number of samples having cnidarian reads and the read counts were significantly higher in the body surface scraping samples than in the gut content samples, regardless of careful rinsing of the body surface. These results indicate that the cnidarian DNAs are most likely found because of cross contamination from the body surface and/or environment. 18S rDNA read counts of copepod and tunicate in the gut contents were greater than or comparable with those in the body surface scraping samples, which may correspond to the previous observations of fecal pellets and larvacean houses in the leptocephali gut. Thus, the present study supports previous implications that leptocephali utilize detritus materials, so called marine snow.

Assessing the Speciation of a Cold Water Species, Japanese Sand Lance Ammodytes personatus, in the Northwestern Pacific by AFLP Markers.

The use of molecular techniques in biodiversity research increasingly results in the recognition of multiple divergent mitochondrial DNA (mtDNA) lineages below the morphospecies level. However, the overlapping distribution of multiple divergent lineages raises the question of whether some of these lineages are in fact cryptic species. Assessing the status of these divergent lineages, delimiting evolutionarily significant units (ESUs), and identifying the dominant evolutionary and ecological drivers are critical components of successful wildlife conservation and management strategies. Amplified fragment length polymorphism (AFLP) markers were applied to characterize the phylogeography pattern of a cold water species, the Japanese sand lance Ammodytes personatus, in warm and cold ocean currents. A total of 211 individuals sampled from 12 populations through the species' range, including samples from Kuroshio Current, Oyashio Current, Tsushima Current, and Yellow Sea, were analyzed. The Bayesian assignment probability test and Neighbor joining (NJ) analysis divided these populations into two genetically and geographically distinct clades (northern and southern clades) characterized by different sea surface temperatures. The incongruence between nuclear clades and previous mitochondrial lineages suggested that A. personatus is indeed composed of at least two genetically divergent cryptic species. Pleistocene glaciation isolation after secondary contact, local thermal adaptation, and isolation by distance may explain the observed geographic pattern of two cryptic species and genetic structure within clades.

Lack of population genetic differentiation of a marine ovoviviparous fish Sebastes schlegelii in Northwestern Pacific.

Sebastes schlegelii is one of the fishes that aggregate around drifting seaweed during early development. To examine the population genetic structure of S. schlegelii, a 452-bp fragment of the mtDNA control region was sequenced and used to interpret life history characteristics and larval dispersal strategy. Two-hundred and twenty-one individuals from 13 sites across the entire range of S. schlegelii in China, Japan and Korea were analyzed. A neighbor-joining tree and network showed that there were no significant genealogical structures corresponding to sampling locations. AMOVA, pair-wise FST and exact test revealed no significant genetic differentiation among locations. The migration rate among locations was high based on the result of LAMARC. We conclude that larval dispersal with drifting seaweed and the current environmental factors may play an important role in shaping the contemporary phylogeographic pattern and genetic homogeneity of S. schlegelii.