Advances in next-generation sequencing techniques have allowed for the generation of genome-wide sequence data, to gain insight into the dynamics influencing genetic structure and the local adaptation of marine fish. Here, using genotyping-by-sequencing (GBS) technique, we identified 31,119 single nucleotide polymorphisms (SNPs) for Sebastiscus marmoratus in 59 individuals from three populations in Chinese coastal waters. Based on all SNPs, there was little evidence of genetic differentiation among populations. However, outlier tests revealed 329 SNPs putatively under divergent selection across populations. Structural and phylogenetic topology analyses based on the outliers showed clear genetic differentiation among populations. Gene Ontology (GO) annotation results revealed that most of these outliers are known or hypothesized to be involved in metabolic process. Together with previous work using mitochondrial cytochrome b sequences, the present results further suggest that the population structure is strongly influenced by locally adaptive pressure. Overall, adaptive evolution in a heterogeneous environment plays an important role in inducing genetic differentiation among local populations. This study increases understanding of the factors (including gene flow and local adaptation) promoting and constraining population genetic differentiation in marine organisms.Inferring the degree of genetic differentiation among populations of marine fish species is key to successfully managing fishery resources, allowing the identification of management units, assignment of individuals to geographic regions, and detection of product mislabeling and fraud [1][2][3] . Many marine fish species are distributed widely across heterogeneous landscapes, and across these ranges, natural selection can promote genetic differentiation and local adaptation. However, the evolution of ecological divergence may be impeded if high rates of migration homogenize the gene pool among populations 4, 5 . Owing to its distinct geographical features with a series of marginal seas 6 , its wide latitudinal range and its complex geological history 7 , the Northwestern Pacific (NWP) is a good region to investigate the patterns of genetic population structuring in marine fish species. Numerous studies on the phylogeographic patterns of marine fish based on neutral markers have been unable to detect genetic structures over such large geographic distances 8 (e.g., Scomberomorus niphonius 9 , Larimichthys polyactis 10 , Trachurus japonicas 11 ), even for species with weak migration abilities (e.g., Hexagrammos otakii 12 , Sebastes schlegelii 13 ). This failure is understandable for fish species with strong migration ability, because the gene flow is unimpeded by ecological homogeneity and a lack of clear dispersal barriers in marine habitats 14,15 . However, owing to local adaptation and a lack of genetic exchange 16,17 , this observation makes less intuitive sense for settled or short-distance-migration fish species. Local adaptation plays an important r...