Identification of candidate genes associated with body weight in swimming crab (Portunus trituberculatus) according to conjunctive analyses of BSA-seq and RNA-seq

The swimming crab (Portunus trituberculatus) is a valuable crustacean species extensively cultivated in China, owing to its fast development, palatability, and rich nutrition. However, the genetic factors regulating its growth remain elusive, impeding the progress of breeding initiatives for this economically important species. In this study, we combined bulk segregation analysis sequencing (BSA-seq) and transcriptome sequencing (RNA-seq) to identify potential candidate genes related to the body weight of P. trituberculatus using a full-sib family. BSA-seq successfully identified five candidate intervals on chromosomes 31, 38, and 47, spanning a total length of 9.76 Mb, which encompassed 137 annotated genes. By conducting RNA-seq analysis on swimming crab individuals in two extreme weight groups, 1581 genes were mined as differentially expressed genes (DEGs), involving 357 upregulated and 1224 down-regulated genes. GO enrichment analysis demonstrated that DEGs were significantly enriched in proteolysis, obsolete oxidation-reduction process, nucleus, cytoplasm, molecular function, and metal ion binding. KEGG enrichment analysis showed that the proteasome signaling pathway could play a prominent role in the growth difference of P. trituberculatus. The conjunctive analyses of BSA-seq and RNA-seq led to the identification of eight critical DEGs, including peripheral plasma membrane protein CASK (CASK), cytoplasmic dynein 2 heavy chain 1 (Dyh1b), RNA-binding protein 5 (Rbm5), structural maintenance of chromosomes protein 3 (Smc3), hemicentin-2 (Hmcn1), A-kinase anchor protein 13 (Akap13), zinc finger MYM-type protein 2 (Zmym3), and extracellular matrix protein FRAS1 (Fras1), which exert influence over the growth of P. trituberculatus. This study not only enhances our understanding of the molecular mechanism governing the growth of P. trituberculatus but also provides valuable insights for future aquaculture programs aimed at selective breeding of this species.

Automated summary

This study combined BSA-seq and RNA-seq to identify potential candidate genes related to the body weight of the swimming crab. BSA-seq successfully identified five candidate intervals, while RNA-seq analysis yielded 1581 differentially expressed genes. Enrichment analysis showed that these genes were involved in proteolysis, nucleus, and molecular function. Eight critical differentially expressed genes were also identified.