The Nao-zhou stock large yellow croaker (Larimichthys crocea) is a unique economic seawater fish species in China and exhibits significant dimorphism in both male and female phenotypes. Cultivating all-female seedlings can significantly improve breeding efficiency. To accelerate the cultivation process of all female seedlings of this species, it is necessary to deeply understand the regulatory mechanisms of sexual differentiation and gonadal development. This study used Illumina high-throughput sequencing to sequence the transcriptome of the testes and ovaries of Nao-zhou stock large yellow croaker to identify genes and molecular functions related to sex determination. A total of 10,536 differentially expressed genes were identified between males and females, including 5682 upregulated and 4854 downregulated genes. Functional annotation screened out 70 important candidate genes related to sex, including 34 genes highly expressed in the testis (including dmrt1, foxm1, and amh) and 36 genes highly expressed in the ovary (including gdf9, hsd3b1, and sox19b). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis found that differentially expressed genes were significantly enriched in nine signaling pathways related to sex determination and gonadal development, including steroid hormone biosynthesis, MAPK signaling pathway, and the TGF-beta signaling pathway. By screening sex-related differentially expressed genes and mapping protein–protein interaction networks, hub genes such as dmrt1, amh, and cyp19a1a were found to be highly connected. The expression levels of 15 sex-related genes, including amh, dmrt1, dmrt2a, foxl1, and zp3b, were determined by qRT–PCR and RNA sequencing. This study screened for differentially expressed genes related to sex determination and differentiation of Nao-zhou stock large yellow croaker and revealed the signaling pathways involved in gonad development of male and female individuals. The results provide important data for future research on sex determination and differentiation mechanisms, thereby providing a scientific basis for the cultivation of all-female seedlings.

The Naozhou stock of large yellow croakers (Larimichthys crocea) exhibits unique phenotypic traits and high genetic diversity, making it a valuable resource for selective breeding and genetic conservation in aquaculture. Despite its importance, simple sequence repeat (SSR) markers have not been developed for this stock, which limits efforts in genetic evaluation, breeding optimization, and sustainable utilization of this commercially important species. In this study, 195,263 SSRs were identified from the genome of the Naozhou stock of large yellow croaker, covering a total length of 16,578,990 bp with a density of 288 bp/Mb. Dinucleotide repeats were the most common, with the AC motif being the most prevalent. The frequency of SSR markers ranged from 245.63 to 346.60 per Mb. A total of 30 primer pairs were synthesized, of which 28 pairs (93.3%) successfully amplified clear and reproducible bands in PCR assays. Among these, 28 SSR markers exhibited distinct and reproducible bands following gel electrophoresis. For eight SSR loci, the number of alleles (Na) ranged from 4 to 22 (mean = 11.375), while the effective number of alleles (Ne) ranged from 1.5401 to 10.4727 (mean = 5.6475). The assembled mitochondrial genome (mtDNA) was 16,467 bp in length and comprised 37 genes, including 13 protein-coding genes (PCGs), 22 tRNA genes, and 2 rRNA genes. The total sequence length of the PCGs was 11,431 bp, accounting for 69.4% of the mtDNA. A large portion of the PCGs (5) used incomplete stop codons (e.g., nad2, nad3, cox2), while others used TAA stop codons (e.g., nad6, nad5, TrnT). The mtDNA encoded a total of 3808 codons, with UAA showing the highest relative synonymous codon usage value. The SSR markers and mtDNA data generated in this study provide valuable tools for future genetic breeding and genomic research on the Naozhou stock of large yellow croakers.

This study examined the digestive and metabolic responses of Naozhou (NZ) stock large yellow croaker (Larimichthys crocea) larvae and juveniles under five developmental feeding stages (DAH3, DAH7, DAH12, DAH19, DAH49) to clarify mechanisms of early nutritional adaptation. Digestive enzyme assays, transcriptome sequencing, and metabolomics were integrated to compare physiological changes across diets. Protease activity increased sharply from DAH7–19 with the introduction of rotifers, Artemia, and copepods, while amylase and lipase activities rose at DAH19–49, reflecting enhanced carbohydrate and lipid utilization during transition to formulated feeds. Transcriptomic analysis showed that differentially expressed genes were enriched in pathways involving protein digestion, lipid and energy metabolism, and cell cycle regulation. The metabolomic analysis further highlighted dynamic changes in amino acid, lipid, carbohydrate, and vitamin metabolism, consistent with transcriptomic findings. The integrated analysis suggests that the coordinated modulation of digestive enzyme activities, gene expression, and metabolite profiles enabled a smooth transition from yolk dependency to live prey feeding and a subsequent use of artificial diets. These findings provide new insights into the early nutritional development of NZ large yellow croaker and provide a scientific basis for the improvement of artificial aquaculture seed production.