Microsatellite markers reveal genetic diversity and population structure of Portunus trituberculatus in the Bohai Sea, China

The swimming crab, Portunus trituberculatus, is one of the main aquaculture species in Chinese coastal regions due to its palatability and high economic value. To obtain a better understanding of the genetic diversity of P. trituberculatus in the Bohai Sea, the present study used 40 SSR loci to investigate the genetic diversity and population structure of 420 P. trituberculatus individuals collected from seven populations in the Bohai Sea. Genetic parameters revealed a low level of genetic diversity in the cultured population (SI = 1.374, He = 0.687, and PIC = 0.643) in comparison with wild populations (SI ≥ 1.399, He ≥ 0.692, and PIC ≥ 0.651). The genetic differentiation index (Fst) and gene flow (Nm) ranged from 0.001 to 0.060 (mean: 0.022) and 3.917 to 249.750 (mean: 31.289) respectively, showing a low differentiation among the seven populations of P. trituberculatus. Population structure analysis, phylogenetic tree, and principal component analysis (PCA) demonstrated that the seven groups of P. trituberculatus were divided into four subpopulations (K = 4), but the correlation between genetic structure and geographical distribution was not obvious. These results are expected to provide useful information for the fishery management of wild swimming crabs.

Development of Microsatellite Markers Based on Transcriptome Sequencing and Evaluation of Genetic Diversity in Swimming Crab (Portunus trituberculatus)

P. trituberculatus is an economically important mariculture species in China. Evaluating its genetic diversity and population structure can contribute to the exploration of germplasm resources and promote sustainable aquaculture production. In this study, a total of 246,243 SSRs were generated by transcriptome sequencing of P. trituberculatus. Among the examined 254,746 unigenes, 66,331 had more than one SSR. Among the different SSR motif types, dinucleotide repeats (110,758, 44.98%) were the most abundant. In 173 different base repeats, A/T (96.86%), AC/GT (51.46%), and ACC/GGT (26.20%) were dominant in mono-, di-, and trinucleotide, respectively. GO annotations showed 87,079 unigenes in 57 GO terms. Cellular process, cell, and binding were the most abundant terms in biological process, cellular component, and molecular function categories separately. A total of 34,406 annotated unigenes were classified into 26 functional categories according to the functional annotation analysis of KOG, of which “general function prediction only” was the biggest category (6,028 unigenes, 17.52%). KEGG pathway annotations revealed the clustering of 34,715 unigenes into 32 different pathways. Nineteen SSRs were identified as polymorphic and, thus, used to assess the genetic diversity and structure of 240 P. trituberculatus individuals from four populations in the Bohai Sea. Genetic parameter analysis showed a similar level of genetic diversity within wild populations, and the cultured population indicated a reduction in genetic diversity compared with wild populations. The pairwise F_ST values were between 0.001 and 0.04 with an average of 0.0205 (p < 0.05), suggesting a low but significant level of genetic differentiation among the four populations. Structure analysis demonstrated that the four populations were classified into two groups including the cultured group and other populations. The phylogenetic tree and PCA revealed that a vast number of samples were clustered together and that cultivated individuals were distributed more centrally than wild individuals. The findings contribute to the further assessment of germplasm resources and assist to provide valuable SSRs for marker-assisted breeding of P. trituberculatus in the future.