Genetic diversity in the mtDNA control region and population structure in the small yellow croaker <Emphasis Type="Italic">Larimichthys polyactis</Emphasis>

The genetic diversity and population genetic structure of the small yellow croaker (Larimichthys polyactis) were investigated. One hundred and fourteen individuals were sampled from 8 localities of the Yellow Sea and the northern East China Sea. Genetic variation in DNA sequences were examined from the first hypervariable region (HVR-1) of the mitochondrial DNA control region. High levels of haplotype diversity (h = 0.98 ± 0.87%) in the HVR-1 region were detected, indicating a high level of genetic diverstiy. A total of 84 polymorphic sites were found, and 87 haplotypes were defined. The pairwise nucleotide differences between samples ranged from 3.83 ± 2.19 to 6.56 ± 3.25. The demographic history of L. polyactis was examined by using neutrality tests and mismatch distribution analysis, which indicated a Pleistocene population expansion at about 49,300–197,000 years. The star burst structure of the minimum spanning tree also suggestted a very recent origin for most haplotypes. Hierarchical molecular variance analysis (AMOVA) and conventional population Fst comparisons revealed no significant genetic structure throughout the examined range, which is inconsistent with previous findings based on the morphological and ecological studies. Long-term dispersal and high gene flow likely have contributed to the genetically homogeneous population structure of the species. The knowledge on genetic diversity and genetic structure will be crucial to establish appropriate fishery management stocks for the species.     

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD‐seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high‐throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction‐site‐associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long‐term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top F_ST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an F_ST‐based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.     

Lack of genetic structure in endangered large yellow croaker Larimichthys crocea from China inferred from mitochondrial control region sequence data

The large yellow croaker Larimichthys crocea is an important commercial marine fish species in China. However, information about the population structure of this species is limited. In the present study, mitochondrial DNA (mtDNA) control region was sequenced from four populations of the yellow croaker in the southern Yellow Sea, East China Sea and South China Sea to investigate the genetic structure of this species. A total of 54 haplotypes were detected among 62 individuals of large yellow croaker. High levels of population genetic diversity were observed. Among the four populations, the haplotype diversity was between 0.9895 ± 0.0193 (Xiamen) and 1.0000 ± 0.0524 (Lvsi, Zhoushan). The nucleotide diversity ranged from 0.0208 ± 0.0108 (Xuwen) to 0.0246 ± 0.0138 (Lvsi). The results of AMOVA detected no significant differences among populations. The conventional F_ST statistics were negative and insignificant values. These indicated lack of population genetic structure throughout the Yellow Sea, East China Sea and South China Sea, and random mixing of individuals among the samples. Biological characters of large yellow croaker and lack of physical barrier in the studied area could be the reasons for lack of genetic structure in this species.     

A morphometric investigation of the small yellow croaker (Larimichthys polyactis Bleeker, 1877): evidence for subpopulations on the Chinese coast

In this study the morphometric characteristics based on the truss network system were constructed from the body shapes and used to investigate the population structure of the small yellow croaker (Larimichthys polyactis Bleeker, 1877). From April to June 2013, a total of 431 fish samples were collected from commercial catches acquired with bottom trawlers (20 × 20 mm) in seven main spawning grounds along the coast of China. The results of principal component and discriminant function analyses indicated a clear distinction among the regional stocks and within the Haizhou Bay (HZB) stock. Observed differences primarily involved head measurements, body length and tail size, which were associated with environmental factors. Latitudinal gradient of the stock distribution and the relationship among the stocks provided further evidence for the existence of subpopulations and the re‐division of the East China Sea group. Moreover, the temporal inter‐structure of the HZB stock is a potential reason for the contrast in results from previous studies. Additional thorough studies of migration trajectories are required for a conclusive and improved understanding of the population structure.     

Genetic analysis of the populations of Japanese anchovy Engraulis japonicus from the Yellow Sea and East China Sea based on mitochondrial cytochrome b sequence

Japanese anchovy Engraulis japonicus was not only a commercial species but also a key link in the marine food chain in the Yellow Sea and East China Sea. However, in recent years, E. japonicus experienced a deep recession caused by over fishing. Therefore, determination of population genetic structure is essential to underpin resource recovery and to aid delineating and monitoring populations for fishery management. In this study, sequence variation of a 402-bp fragment of mitochondrial cytochrome b (cyt b) gene from ten populations of E. japonicus in the Yellow Sea and East China Sea was determined, detecting 73 haplotypes. Haplotype diversities and nucleotide diversities for the ten populations ranged from 0.4706 ± 0.1177 (D) to 0.9935 ± 0.210 (WH), and from 0.0014 ± 0.0013 (D) to 0.0090 ± 0.0053 (WH), respectively. Analysis of molecular variance (AMOVA) showed no significant differences at all hierarchical levels, and almost all F_ST values were non-significant, indicating no significant population genetic structure in the Yellow Sea and East China Sea. Two indistinct clades with low support values were detected and the analysis of haplotype frequency distribution showed no significant geographic differences among populations. These results supported the null hypothesis that E. japonicus within the Yellow Sea and East China Sea constituted a panmictic mtDNA gene pool. Based on the present study, the existence of separate genetic stocks in the Yellow Sea and East China Sea were not detected. Strong dispersal capacity of planktonic larvae, and ocean current could be the reasons for genetic homogeneity in this species in the Yellow Sea and East China Sea.     

Factors affecting morphological development of the sagittal otolith in juvenile and adult small yellow croaker (Larimichthys polyactis Bleeker, 1877)

The morphology and morphometrics of the sagittal otoliths of small yellow croaker (Larimichthys polyactis) from the southern Yellow Sea were investigated. Study objectives were to evaluate the shape variability and morphometric variables of sagittae of juveniles and adults as related to developmental changes and habit shift. A total of 152 fish were sampled from April to June of 2012 and 2013, along the coastal waters of the Lüsi spawning ground in the southern Yellow Sea. Changes in otolith shapes from the juvenile to the adult were presented with the rim development through the entire‐lobed‐entire transition and with the curvature of the cauda toward the ventral margin. The otolith elongation in the juvenile stage occurred at 10–20 mm standard length (SL) and was likely associated with the formation of otolith accessory growth centers from larvae to juvenile. The L. polyactis sagittal otoliths acquired their definitive shape at 130 mm SL maturity. Ontogeny on otolith shape might be related, for example, to the factors of metamorphic development, feeding habitat, and ambient water salinity, which varied throughout the growth of L. polyactis.     

Identification of a large SNP dataset in Larimichthys crocea using specific‐locus amplified fragment sequencing

The large yellow croaker, Larimichthys crocea, is a commercially important drum fish (Family: Sciaenidae) native to the East and South China Sea. Habitat deterioration and overfishing have led to significant population decline and the collapse of its fishery over the past decades. Today, the market supply of L. crocea depends solely on stocks produced in hatcheries and farms. Common issues that occur in the culture of L. crocea include germplasm degradation, precocious puberty, elevated disease susceptibility and growth retardation. In this study, we employed SLAF‐seq (specific‐locus amplified fragment sequencing) technology to identify single nucleotide polymorphism (SNP) loci across the L. crocea genome. Sixty samples were selected for SLAF analysis out of 1000 progeny in the same cohort of a cultured stock. Our analysis obtained a total of 151 253 SLAFs, of which 65.88% (99 652) were identified to be polymorphic, scoring a total of 710 567 putative SNPs. Further filtration resulted in a final panel of 1782 SNP loci. The data derived from this work could be beneficial for understanding the genetics of complex phenotypic traits as well as for developing marker‐selection‐assisted breeding programs in L. crocea.     

Significant genetic differentiation between the Yellow Sea and East China Sea populations of cocktail shrimp Trachypenaeus curvirostris revealed by the mitochondrial DNA COI gene

The cocktail shrimp Trachypenaeus curvirostris is an ecologically and economically important shrimp species in the Yellow Sea and East China Sea. However, there is no information about its population genetic structure. The population genetic structure and level of gene flow of T. curvirostris from the Yellow Sea and the East China Sea were studied with a 658-bp segment of mtDNA COI gene. In total, 85 individuals were collected from five locations and 13 haplotypes were obtained. The genetic variation of COI gene in five populations was moderate, giving an overall haplotype diversity of 0.6888 ± 0.0432 and nucleotide diversity of 0.0069 ± 0.0038. Conflicting to our expectation, significant genetic differentiation was detected in this species. The result revealed two genetically divergent lineages, displaying clear different geographical distributions in the studied area. The significant genetic differentiation between the Yellow Sea and East China Sea populations might be caused by the Yangtze River outflow. Mismatch distribution revealed that T. curvirostris had undergone population range expansion, possibly before 103,400–109,700 years ago in the last interglaciation, rejecting the sudden demographic expansion.     

Upwelling and eddies affect connectivity among local populations of the goldeye rockfish,Sebastes thompsoni (Pisces,Scorpaenoidei)

The goldeye rockfish, Sebastes thompsoni, commercial rockfish catch in the Northwest Pacific Ocean, may influence its population structure. To clarify the population genetic structure of Korean S. thompsoni and its degree of hybridization with the most close species, Sebastes joyneri, we analyzed a mitochondrial (mt) DNA control region and eleven polymorphic microsatellite (ms) loci. S. joyneri individuals were clearly distinguished from S. thompsoni by the mtDNA control region and ms loci results, with single interspecific hybridization between two species suggesting no impact on genetic structure of S. thompsoni. Analysis of mtDNA revealed no population structure within S. thompsoni, suggesting the survival of a single population in southern refugia during the glacial period. The ms loci results, in contrast, showed two genetically distinct clusters within S. thompsoni: One was predominant throughout Korean coasts (from the Yellow Sea, via the Korea Strait to the East Sea); the other was predominant at Dokdo Island in the East Sea; and both occurred in similar ratios at Wangdolcho Reef in the East Sea. A possible factor that restricts gene flow between Korean coastal and offshore populations in the East Sea may be related to the complex oceanic current patterns such as eddies and upwelling, which represent impermeable barriers to population connectivity for this species. Our findings highlight that these two populations might be representative of two separate stock within Korean waters and maintain their geographically related genetic structure.     

Single nucleotide polymorphism discovery in albacore and Atlantic bluefin tuna provides insights into worldwide population structure

The optimal management of the commercially important, but mostly over‐exploited, pelagic tunas, albacore (Thunnus alalunga Bonn., 1788) and Atlantic bluefin tuna (BFT; Thunnus thynnus L., 1758), requires a better understanding of population structure than has been provided by previous molecular methods. Despite numerous studies of both species, their population structures remain controversial. This study reports the development of single nucleotide polymorphisms (SNPs) in albacore and BFT and the application of these SNPs to survey genetic variability across the geographic ranges of these tunas. A total of 616 SNPs were discovered in 35 albacore tuna by comparing sequences of 54 nuclear DNA fragments. A panel of 53 SNPs yielded F_ST values ranging from 0.0 to 0.050 between samples after genotyping 460 albacore collected throughout the distribution of this species. No significant heterogeneity was detected within oceans, but between‐ocean comparisons (Atlantic, Pacific and Indian oceans along with Mediterranean Sea) were significant. Additionally, a 17‐SNP panel was developed in Atlantic BFT by cross‐species amplification in 107 fish. This limited number of SNPs discriminated between samples from the two major spawning areas of Atlantic BFT (F_ST = 0.116). The SNP markers developed in this study can be used to genotype large numbers of fish without the need for standardizing alleles among laboratories.     

Population structure and genetic basis of the agronomic traits of upland cotton in China revealed by a genome‐wide association study using high‐density SNPs

Gossypium hirsutum L. represents the largest source of textile fibre, and China is one of the largest cotton‐producing and cotton‐consuming countries in the world. To investigate the genetic architecture of the agronomic traits of upland cotton in China, a diverse and nationwide population containing 503 G. hirsutum accessions was collected for a genome‐wide association study (GWAS) on 16 agronomic traits. The accessions were planted in four places from 2012 to 2013 for phenotyping. The CottonSNP63K array and a published high‐density map based on this array were used for genotyping. The 503 G. hirsutum accessions were divided into three subpopulations based on 11 975 quantified polymorphic single‐nucleotide polymorphisms (SNPs). By comparing the genetic structure and phenotypic variation among three genetic subpopulations, seven geographic distributions and four breeding periods, we found that geographic distribution and breeding period were not the determinants of genetic structure. In addition, no obvious phenotypic differentiations were found among the three subpopulations, even though they had different genetic backgrounds. A total of 324 SNPs and 160 candidate quantitative trait loci (QTL) regions were identified as significantly associated with the 16 agronomic traits. A network was established for multieffects in QTLs and interassociations among traits. Thirty‐eight associated regions had pleiotropic effects controlling more than one trait. One candidate gene, Gh_D08G2376, was speculated to control the lint percentage (LP). This GWAS is the first report using high‐resolution SNPs in upland cotton in China to comprehensively investigate agronomic traits, and it provides a fundamental resource for cotton genetic research and breeding.     

Genetic structure in the marbled rockfish (Sebastiscus marmoratus) across most of the distribution in the northwestern Pacific

The marbled rockfish (Sebastiscus marmoratus) is an important species that is widely distributed across the marginal seas of the northwestern Pacific. Several kinds of DNA markers, such as amplified fragment length polymorphisms (AFLP), mitochondrial DNA (mtDNA) and single nucleotide polymorphisms (SNPs), have been used to assess the population genetic characteristics of this species in previous studies. However, there have been no genetic profiling studies to cover the entire distribution. Here, six highly polymorphic microsatellite markers were applied to 18 populations across a spatially large area (from the Sea of Japan in the north to the South China Sea in the south). The results showed that the whole population exhibits high, stable genetic diversity, low relatedness and large effective population size. There was very weak genetic differentiation overall and no isolation by distance occurred among the populations, which may be attributed to significant contemporary gene flow, each generation. The above results may lead us to infer that the recent S. marmoratus population structure could be considered as one large integrated population. The present study will be beneficial to population conservation and fisheries management of S. marmoratus and will lead to better insights into the genetic characteristics of other species that belong to the Sebastidae family.     

Analysis of whole-genome changes and characteristics of Linepithema humile (Hymenoptera: Formicidae) in the Republic of Korea

The Argentine ant, Linepithema humile (Mayr, 1868), native to the Paraguay River in South America, was first discovered in Korea in 2019. With increasing reports of L. humile, its genetic variation according to domestic growth colonies and its effects on the domestic ecosystem should be studied. Here, the genomes of L. humile specimens found at three locations in Busan were analyzed for genetic changes. First, morphological observation of L. humile samples collected from the three sites showed no phenotypic differences among them. Next, single-nucleotide polymorphism (SNP)/insertion and deletion (INDEL) analyses on the genomic DNA from the three groups showed that in the 1-NIE sample, the most frequent mutations were G → A and C → T. The mutations A → C and T → G were confirmed in the 14-NIE and 19-NIE samples. Although the number of SNPs in the N section was small, sequences of 4681 bp (1-NIE), 4217 bp (14-NIE) and 4631 bp (19-NIE) in length were identified. From the INDEL length distribution of the three samples, most changes were associated with insertions and deletions of 1–2 bp. However, no heterogeneity was found in the population samples analyzed based on SNP data. Comparative analysis of the SNPs investigated revealed that 760 819 (11.72%) of the total 6 492 517 SNPs were found in common, demonstrating that the three groups analyzed had different genetic backgrounds. Overall, we have developed a method for analyzing the genetic diversity of L. humile invading the Republic of Korea, precisely classified its genetic characteristics and obtained genomic data on interspecies mutations, according to the local environment.     

Genetic differentiation between populations of swimming crab <Emphasis Type="Italic">Portunus trituberculatus</Emphasis> along the coastal waters of the East China Sea

Portunus trituberculatus is a commercially important species widely spread in the East China Sea. Intraspecific variation of the mitochondrial DNA cytochrome oxidase subunit I (mtDNA COI) gene was investigated in 213 individuals from six localities (Changjiang Estuary, Shengsi Islands, Zhoushan Islands, Dongtou Islands, Dinghai Bay, and Quanzhou Bay) ranging from north (31°21′N) to south (24°55′N) coastal waters of the East China Sea. Overall, a total of 27 mtDNA haplotypes and 21 variable sites were detected in the 787 bp segment of COI gene. Analysis of mtDNA COI sequence data revealed that crabs from the six localities were characterized by moderately high haplotypic diversity (h = 0.787 ± 0.026), while sequence divergence values between haplotypes were relatively low (π = 0.00241 ± 0.00098). Each population was characterized by a single most frequent haplotype, shared among all six localities, and a small number of rare ones, typically present in only one or two individuals and representative of a specific population. However, neither the neighbor-joining tree nor the minimum spanning network (MSN) based on the haplotype data exhibited geographical patterns of the six populations. Mismatch distribution analysis of P. trituberculatus individuals sampled from the six localities suggested that sudden population expansion might have occurred in CJ and SS population that might be consistent with over-exploitation of the swimming crab. Analysis of molecular variance (AMOVA) and F _ST statistics showed that significant genetic differentiation existed among the SS, ZS, DT, DH, and QZ populations, suggesting that gene flow might be reduced, even between the geographically close sites, despite the high potential of dispersal. The possible causes of the observed genetic heterogeneity among the P. trituberculatus populations and the potential applications of the mtDNA COI marker in the artificial breeding and fisheries management are discussed. Handling editor: C. Sturmbauer     

A panel of polymorphic EST-derived microsatellite loci for the small yellow croaker (Larimichthys polyactis)

Small yellow croaker (Larimichthys polyactis) is an important economic species of marine fishery. We developed and evaluated simple sequence repeat (SSR) markers from expressed sequence tags (ESTs) of Pseudosciaena crocea, Paralichthys olivaceus and Psetta maxima. Characteristics of nine EST–SSR loci were investigated using 46 L. polyactis individuals. The number of alleles per locus ranged from 2 to 6. The observed heterozygosity ranged from 0.0652 to 0.7391, while the expected heterozygosity ranged from 0.0638 to 0.7754. Seven loci departed from Hardy–Weinberg equilibrium (P < 0.01) significantly. These loci and markers will be useful for population genetics and systemic evolution among species of small yellow croaker.     

Outlier SNP markers reveal fine‐scale genetic structuring across European hake populations (Merluccius merluccius)

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome‐wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large‐ and fine‐scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F_CT = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F_CT range 0.275–0.705) and fine‐scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.     

Genetic diversity of Chinese cattle revealed by Y‐SNP and Y‐STR markers

With its vast territory and complex natural environment, China boasts rich cattle genetic resources. To gain the further insight into the genetic diversity and paternal origins of Chinese cattle, we analyzed the polymorphism of Y‐SNPs (UTY19 and ZFY10) and Y‐STRs (INRA189 and BM861) in 34 Chinese cattle breeds/populations, including 606 males representative of 24 cattle breeds/populations collected in this study as well as previously published data for 302 bulls. Combined genotypic data identified 14 Y‐chromosome haplotypes that represented three haplogroups. Y2‐104‐158 and Y2‐102‐158 were the most common taurine haplotypes detected mainly in northern and central China, whereas the indicine haplotype Y3‐88‐156 predominates in southern China. Haplotypes Y2‐108‐158, Y2‐110‐158, Y2‐112‐158 and Y3‐92‐156 were private to Chinese cattle. The population structure revealed by multidimensional scaling analysis differentiated Tibetan cattle from the other three groups of cattle. Analysis of molecular variance showed that the majority of the genetic variation was explained by the genetic differences among groups. Overall, our study indicates that Chinese cattle retain high paternal diversity (H = 0.607 ± 0.016) and probably much of the original lineages that derived from the domestication center in the Near East without strong admixture from commercial cattle carrying Y1 haplotypes.     

The curious case of Hermodice carunculata (Annelida: Amphinomidae): evidence for genetic homogeneity throughout the Atlantic Ocean and adjacent basins

Over the last few decades, advances in molecular techniques have led to the detection of strong geographic population structure and cryptic speciation in many benthic marine taxa, even those with long‐lived pelagic larval stages. Polychaete annelids, in particular, generally show a high degree of population divergence, especially in mitochondrial genes. Rarely have molecular studies confirmed the presence of ‘cosmopolitan’ species. The amphinomid polychaete Hermodice carunculata was long considered the sole species within its genus, with a reported distribution throughout the Atlantic and adjacent basins. However, recent studies have indicated morphological differences, primarily in the number of branchial filaments, between the East and West Atlantic populations; these differences were invoked to re‐instate Hermodice nigrolineata, formerly considered a junior synonym of H. carunculata. We utilized sequence data from two mitochondrial (cytochrome c oxidase subunit I, 16S rDNA) markers and one nuclear (internal transcribed spacer) marker to examine the genetic diversity of Hermodice throughout its distribution range in the Atlantic Ocean, including the Mediterranean Sea, the Caribbean Sea, the Gulf of Mexico and the Gulf of Guinea. Our analyses revealed generally low genetic divergences among collecting localities and between the East and West Atlantic, although phylogenetic trees based on mitochondrial data indicate the presence of a private lineage in the Mediterranean Sea. A re‐evaluation of the number of branchial filaments confirmed differences between East and West Atlantic populations; however, the differences were not diagnostic and did not reflect the observed genetic population structure. Rather, we suspect that the number of branchial filaments is a function of oxygen saturation in the environment. Our results do not support the distinction between H. carunculata in the West Atlantic and H. nigrolineata in the East Atlantic. Instead, they re‐affirm the older notion that H. carunculata is a cohesive species with a broad distribution across the Atlantic Ocean.