Integrative QTL Mapping and Transcriptomic Profiling Identify Growth-Associated QTL and Candidate Genes in Hong Kong Catfish (Clarias fuscus)

1. Introduction

Growth traits are among the most critical economic features in aquaculture species. Enterprises focusing on fish with superior growth performance can benefit from reduced rearing cycles, lower production costs, and significant increases in yield and economic returns. Consequently, selecting and breeding fish species with exceptional growth potential has become a central goal for the sustainable development of the aquaculture industry [1]. Growth traits are regulated by numerous small-effect polygenes, known as quantitative trait loci (QTL). While most of these loci have minor effects, a few key genes exert a substantial influence on growth [2]. Several genes involved in growth regulation have been identified in fish. Genes within the growth hormone/insulin-like growth factor (GH/IGF) axis , such as growth hormone (gh), growth hormone receptor (ghr), insulin-like growth factor 1 (igf-I), insulin-like growth factor 2 (igf-II), and insulin-like growth factor binding proteins (igfbps), play critical roles in regulating metabolic processes and cell growth [3]. Additionally, members of the myogenic regulatory factors (MRFs) gene family, including myogenin (myog), myogenic differentiation 1 (myod), and myogenic factor 5 (myf-5), are primarily expressed in vertebrate skeletal muscle and are crucial for muscle development through the regulation of cell proliferation and activation of muscle-specific genes [4]. Other genes, such as melanocortin receptor-4 (mc4r), influence growth by regulating energy metabolism, while gastric hunger hormone (ghrelin) plays a role in growth regulation by enhancing feeding and digestive efficiency [5,6]. These genes interact through a complex signaling network that collectively influences growth-related phenotypic variation in fish. Given the species-specific nature of growth regulation, screening growth-related genes across different species is essential for gaining a comprehensive understanding of the genetic mechanisms behind growth regulation in fish.

QTL mapping is a crucial tool for investigating the genetic basis of economic traits in aquaculture species, playing a key role in guiding the identification of candidate genes [7]. Traditional breeding methods rely on the careful selection of individuals with desirable characteristics followed by successive generations of breeding to produce seedlings with superior economic traits. However, this process is time-consuming, costly, and highly susceptible to external factors such as environmental influences, limiting its practical application in production [8]. In contrast, QTL mapping enables the identification of complex traits by pinpointing specific genomic regions, facilitating the early selection of superior individuals. This approach not only significantly shortens the breeding cycle but also improves breeding efficiency and enhances economic returns. Currently, QTL mapping is widely applied in breeding studies of aquaculture species, addressing key economic traits such as sex determination[9,10], hypoxia tolerance [11,12] and body color[13,14]. For growth traits, QTL mapping studies have been conducted in various aquaculture species, including common carp (Cyprinus carpio)[15], mirror carp (Cyprinus carpiovar. specularis)[16]. These studies have enriched genetic resources for fish breeding and provided valuable insights into the genetic basis of growth traits in fish.

While QTL mapping is valuable for identifying chromosomal regions linked to growth traits by analyzing the genetic structure of economic traits, it primarily locates broader chromosomal regions where candidate genes are mapped. Without further fine-mapping, its ability to provide detailed candidate gene information is limited [17]. On the other hand, RNA-seq offers insights into gene expression dynamics, revealing how gene expression varies across individuals of different sizes [18,19]. In recent years, transcriptomics analysis has been widely used to identify differentially expressed genes (DEGs) associated with growth traits in aquaculture species, for example rainbow trout (Oncorhynchus mykiss)[20], mandarin fish (Siniperca chuatsi)[21], and fine flounder (Paralichthys adspersus)[22]. Combining QTL mapping with RNA-seq allows for the narrowing of the candidate gene range through co-mapping strategies like eQTL analysis, leading to more accurate gene mapping. For instance, a co-analysis in giant grouper (Epinephelus Lanceolatus) identified 34 candidate genes related to growth and six growth-related QTL [23]. Similarly, in Pacific whiteleg shrimp (Litopenaeus vannamei), four candidate genes were identified from nine growth-related QTLs using conjoint analysis [24]. This integrated approach improves the identification of growth-related candidate genes and enhances our understanding of the genetic mechanisms underlying growth traits in aquaculture species.

The Hong Kong catfish (Clarias fuscus) is the only native species of the Clariidae family in China [25]. This species is highly valued for its desirable traits, including strong adaptability, tender meat, and a high content of essential amino acids with low fat content [26]. It has become one of the important economic aquaculture species in southern China [27]. However, the development of the aquaculture industry for this species is hindered by the lack of targeted selective breeding programs. Therefore, there is an urgent need to enhance key economic traits, particularly growth, and to select seedlings with superior growth potential. Currently, the whole genome sequences and high-density genetic linkage map of C. fuscus have been constructed [25,26,28]. The acquisition of these whole genome sequences and high-density genetic maps has laid the groundwork for identifying economic traits in this species. However, the genetic mechanisms underlying growth traits remain poorly understood. In this study, we performed QTL mapping for eight growth traits of C. fuscus using the established genetic map [28] and integrated muscle RNA-seq data to identify growth-related candidate genes. The findings will provide a scientific foundation for understanding the genetic basis of growth traits in C. fuscus and contribute to the development of C. fuscus lines with improved growth performance.

2. Materials and Methods

2.1. Ethics Statement

All experimental protocols in this study were conducted in accordance with the standards of the Guangdong Ocean University Animal Research and Ethics Committee (NIH Pub. No.85-23, revised in 1996). The subjects of this study are not endangered or protected species.

2.2. Sample Preparation and Phenotyping

The family group used for QTL mapping of growth traits in C. fuscus in this study originates from the same family group previously employed for genetic linkage map construction [28]. Specifically, a mapping family of Hong Kong catfish was provided by the Guangxi Hongtai Aquatic Product Farm in Guangxi Province, China. This family was established by randomly selecting 200 F₁ individuals for the construction of a genetic linkage map [28], which was subsequently used for QTL mapping of growth traits. In this study, eight growth-related traits, including body weight (BW), body height (BH), body length (BL), body width (BWI), orbital diameter (OD), caudal peduncle length (CPL), caudal peduncle height (CPH), and pre-dorsal length (PDL), were measured in 200 individuals of a full-sibling F₁ family group of C. fuscus for which the genetic map was constructed [28]. The raw data for these traits were initially preprocessed using Microsoft Excel. Subsequent analyses, including descriptive statistics, correlation analysis, and analysis of variance (ANOVA), were conducted using SPSS software.

For transcriptome sequencing, a mixed family group of C. fuscus was provided by the Guangxi Hongtai Aquatic Product Farm, Guangxi Province, China. Selected same-size 1-month-old C. fuscus were transferred from the farm to the freshwater aquaculture base of Guangdong Ocean University. They were cultured in the same pond and fed pellet feed three times a day. At 6 months of age, 254 individuals were randomly selected for body weight measurement. From these, the 9 largest and 9 smallest individuals were chosen to form 3 big-sized groups (B1, B2, B3) and 3 small-sized groups (S1, S2, S3). The average body weight of the large-sized groups was 100.13±12.96g, while that of the small-sized groups was 32.29±7.48g. The fish were fasted for 24 h before sample collection. Before dissection, the fish were anesthetized with eugenol (Shandong Keyuan Biochemical Co., Ltd., China). Muscle tissues were harvested, rapidly frozen in liquid nitrogen and stored at -80°C. For each group, muscle tissues from three individuals were mixed to a single sample for subsequent sequencing and analysis.

2.3. QTL Mapping for Growth Traits

QTL mapping for for growth traits was conducted using the genetic linkage map of C. fuscus [28]. Eight growth traits were analyzed using the Permutation Test (PT) in MapQTL 6.0. The PT was performed with 1000 permutations to determine the LOD threshold for each trait at a 95% confidence level. The results from the PT were combined with the actual QTL significance levels obtained from the analysis to establish uniform thresholds for the eight growth traits. Interval mapping (IM) in MapQTL 6.0 was then used to analyze these traits, and the corresponding QTL segments were identified based on the established thresholds [29].

2.4. Screening of Candidate Genes for Growth Traits QTL and Functional Enrichment Analysis

Candidate genes and functional annotations within the QTL intervals for growth traits were identified using the reference genome of C. fuscus (https://www.ncbi.nlm.nih.gov/datasets/genome/GCA_030347435.1/). Functional genes within the QTL intervals identified through linkage analysis were selected based on this reference genome. The corresponding sequences within each QTL region were extracted and annotated by BLAST analysis. Subsequently, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed on the annotated genes using TBtools software. The resulting visualizations were generated using the R package clusterProfiler [30].

2.5. Transcriptome Library Creation and Raw Data Processing

Total RNA was extracted from six samples in the big-sized and small-sized individual groups. mRNA was randomly fragmented using NEB Fragmentation Buffer. The first strand of cDNA was synthesized using random oligonucleotides as primers in the M-MuLV reverse transcriptase system. The RNA strand was then degraded with RNaseH, and the second strand of cDNA was synthesized with dNTPs in the DNA polymerase I system. The purified double-stranded cDNA underwent end repair, A-tail addition, and ligation of sequencing junctions. cDNAs of approximately 250-300 bp were selected using AMPure XP beads, followed by PCR amplification, and the PCR products were purified again. The libraries were then constructed [30]. Raw data were filtered to remove adapters and low-quality reads (those with an average read quality score below 5), resulting in clean reads [31]. These high-quality reads were aligned to the reference genome using HISAT2 software, ensuring fast and accurate alignment to obtain the localization information of the reads on the reference genome(https://www.ncbi.nlm.nih.gov/datasets/genome/GCA_046453815.1/) [32]. Finally, new transcripts were assembled using StringTie software [33].

2.6. Screening of DEGs and Functional Enrichment Analysis

The big-sized individual group was used as the experimental group, while the small-sized individual group served as the control group. Differential expression analysis between the two groups was conducted using the DESeq2 R package (version 1.6.3), identifying differentially expressed genes (DEGs) with a p-value < 0.05 and |fold change| ≥ 1. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of the DEGs were performed using TBtools software, and the results were visualized using the R clusterProfiler package.

2.7. QTL and Transcriptome Association Analysis

By integrating the results from the QTL mapping analysis of significant regions and transcriptome differential gene analysis, a set of common genes between the two datasets was identified for further quantitative analysis and expression profiling. For gene expression analysis, expression levels were log2-transformed, and heatmaps were generated using the R package pheatmap. The heatmaps were constructed using a fully interlocked clustering method and the Euclidean distance measurement metric.

2.8. Real-Time Fluorescence Quantitative PCR (qRT-PCR) Validation

The expression of candidate genes in the muscle of C. fuscus was quantified using qRT-PCR. The RNA used for transcriptome sequencing was reverse transcribed into complementary DNA (cDNA) using HiScript III All-in-One RT SuperMix (Vazyme, China), optimized for quantitative PCR (qPCR). qRT-PCR was performed using PerfectStart Green qPCR SuperMix (Code#TG-AQ601-02) on a LightCycler® 480 detection system (Roche Diagnostics Ltd., USA), with each sample analyzed in triplicate. Primers for qPCR were designed using the NCBI tool, with the β-actin gene serving as the internal reference gene. The reaction mixture was prepared in a total volume of 10 μL, consisting of 5.0 μL of 2× PerfectStart Green qPCR SuperMix, 1.0 μL of cDNA, 3.4 μL ddH2O, 0.3 μL of the forward primer, and 0.3 μL of the reverse primer. The reaction procedure was as follows: 1) Preincubation: 94°C for 300s; 2) Amplification: (95°C for 20s, 60°C for 20s, 72°C for 20s) × 40 cycles; 3) Melting: 95°C for 10s, 65°C for 60s, 97°C for 1s; 4) Cooling: 37°C for 300s. The comparative CT (2−(ΔΔCT)) method was used to assess relative gene expression levels. The primers for candidate genes and internal control gene (β-actin) are shown in Table S1.

3. Results

3.1. Morphological Statistics

In this study, eight morphological traits (BW, BL, BH, BWI, OD, CPL, CPH, and PDL) were measured in 200 F1 individuals of C. fuscus. The morphological analyses showed that the frequency distribution of these traits conformed to a normal distribution (Figure 1). The coefficients of variation (CV) for the eight traits ranged from 8.00% to 11.90% (Table S2), indicating substantial variability, which makes these traits suitable candidates for QTL analysis.

3.2. QTL Analysis and Candidate Gene Identification

Using QTL composite interval mapping with the MapQTL 6.0 software, we analyzed eight growth traits of C. fuscus. The analysis identified 17 QTLs (LOD > 3.5) distributed across eight linkage groups, including LG06, LG07, LG09, LG13, LG14, LG17, LG18, and LG20 (Figure 2). The peak positions of these QTLs ranged from 5.62 cM (qCPH-8) to 186.52 cM (qCPH-5). The phenotypic variance explained (PVE) by these QTLs ranged from 8.00% (qBL-4) to 11.90% (qCPH-6), with an average PVE of 9.04% (Table 1). Multiple growth trait QTLs were clustered on LG06, LG09, and LG17. Specifically, three QTLs for CPH, CPL, and BH were identified on LG06, On LG09, four QTLs were identified for BW, BL, PDL, and BWI, with a notable overlap between the significant regions for BW (41.931–45.835 cM) and BWI. Similarly, QTLs for BW, BWI, and PDL were on LG17, with the significant region for BW (13.520–17.098 cM) overlapping with that for BWI. A total of 162 candidate genes were identified within these 17 QTL intervals (Table S3).

3.3. Functional Annotation of QTL Candidate Genes

To explore the functions and associated pathways of the QTL candidate genes, GO and KEGG enrichment analyses were performed. The GO analysis revealed five significantly enriched terms within the biological process (BP) category, including cell adhesion (GO:0007155) and RNA metabolic process (GO:0016070). In the molecular function (MF) category, three terms were significantly enriched, such as calcium ion binding (GO:0005509) and structural molecule activity (GO:0005198). No significant enrichment was observed in the cellular component (CC) category (Figure 3A). KEGG pathway analysis identified nine significantly enriched pathways, including Brite Hierarchies, Propanoate metabolism, and Transcription machinery (Figure 3B).

3.4. Quality Control of Transcriptome Sequencing Data

Following quality control procedures, low-quality reads from the raw sequencing data were removed, resulting in an average of 50,513,416 (98.25%) clean reads for the big-sized group and 44,737,262 (95.63%) clean reads for the small-sized group. The average Q20 and Q30 values for the big-sized individuals were 98.69% and 96.20%, respectively, while for the small-sized individuals, these values were 98.10% and 95.03%, respectively. After quality control, the clean reads from both groups were mapped to the reference genome of the C. fuscus. On average, 90.47% of the clean reads from the big-sized group and 83.57% from the small-sized group were successfully mapped to the reference genome (Table 2). These results indicate that the sequencing results are of high quality and suitable for subsequent analyses.

3.5. Differential Gene Expression Analysis Related to Body Size

Gene expression changes in muscle tissue samples from big-sized and small-sized individuals were compared, identifying a total of 3,824 differentially expressed genes (DEGs) associated with individual size in C. fuscus. Of these, 2,252 genes were up-regulated, and 1,572 genes were down-regulated (p < 0.05 and |log₂(fold change) | ≥ 1, Figure 4A). Table 3 presents the top 20 genes ranked by fold change between the big-sized and small-sized individual groups.

3.6. Functional Classification of DEGs

To further explore the biological significance of the DEGs, GO and KEGG enrichment analyses were conducted. The GO enrichment analysis revealed significant enrichment of DEGs in two functional categories, with 94 subcategories identified within the Biological Process (BP) category and 16 subcategories within the Cellular Component (CC) category. No significant enrichment was observed in the MF category. The DEGs were significantly enriched in terms such as Envelope (GO:0031975), Mitochondrion (GO:0005739), ATP metabolic process (GO:0046034), and Translation (GO:0006412) (Figure 4B). KEGG pathway analysis indicated that DEGs were significantly enriched in nine pathways, including Carbon metabolism, 2-Oxocarboxylic acid metabolism, and Glycolysis/Gluconeogenesis. Among these pathways, the Biosynthesis of Amino Acids pathway was the most significantly enriched, while the Insulin Signaling pathway had the highest number of enriched DEGs (Figure 4C).

3.7. Carbohydrate Metabolism Pathway Analysis

To investigate the molecular mechanisms underlying size variation in C. fuscus, we examined the expression patterns of key DEGs in major carbohydrate metabolism pathways across two groups (Figure 5). The analysis revealed that genes involved in glycolysis and the tricarboxylic acid (TCA) cycle, such as glucose-6-phosphate isomerase (gpi), ATP-dependent 6-phosphofructokinase (pfkm), and fructose-1,6-bisphosphatase isozyme 2 (fbp2), were upregulated in the large-sized individuals. Conversely, genes associated with odd-chain fatty acid catabolism, isoleucine catabolism, and ethanol metabolism, including alcohol dehydrogenase class-3 (adh3) and peroxisomal acyl-coenzyme A oxidase 1 (acox1), exhibited higher expression levels in the small-sized individuals.

3.8. Integration of QTL Mapping and Transcriptome Data

By integrating the results of QTL mapping and transcriptome differential expression analysis, we identified 27 candidate genes associated with growth, including epidermal growth factor receptor (egfr), ataxin-7protein 3 (atxn7), glutamine-rich protein 1 (qrich1), and nucleophosmin (npm1) (Table 4). The expression analysis of the 27 growth candidate genes is shown in Figure 6A. To validate the expression patterns of these growth-related candidate genes, we performed qRT-PCR analysis on five randomly selected genes in muscle tissues from both big-sized and small-sized C. fuscus individuals. The expression patterns observed through qRT-PCR were consistent with those from the transcriptome data, confirming the accuracy and reliability of the transcriptome results (Figure 6B).

4. Discussion

Growth traits are critical economic traits in aquaculture species and have garnered considerable attention [34]. While existing studies have offered valuable insights into the regulatory mechanisms of growth genes in fish, the complexity of growth regulation has limited the identification of related genes. In this study, we identified QTLs associated with growth in C. fuscus, by integrating transcriptome data related to body size, pinpointed candidate genes and relevant pathways. This work provides a foundation for further investigation into the growth regulatory mechanisms in C. fuscus and offers key data for developing fast-growing strains.

4.1. QTL Mapping of Growth Traits in C. fuscus

Numerous studies have demonstrated that QTL mapping is a primary method for identifying quantitative trait loci and understanding the underlying mechanisms of quantitative traits [35]. One critical metric in QTL mapping is the proportion of variance explained (PVE), which quantifies the extent to which a QTL contributes to phenotypic trait variation. A higher PVE value indicates a greater influence of the QTL on the trait, while a lower PVE suggests a lesser effect, implying that other genetic or environmental factors may also significantly shape the trait [36]. In fish, growth traits are classic quantitative traits, with associated QTLs typically exhibiting relatively low PVE values. For example, mirror carp (9.6-24.2%) [16], crucian carp (Carassius auratus) (PVE, 10.1-13.2%) [37], ussuri catfish (Pseudobagrus ussuriensis) (11.9-20.5%) [38], and dusky kob (Argyrosomus japonicus) (9.3-29.5%) [39] all demonstrate this pattern.

Based on the available genetic map of C. fuscus, this study analyzed eight growth traits using QTL mapping. Except for OD, CPL, and BWI, all other traits were mapped to two or more significant QTL intervals, resulting in the identification of 17 growth-related QTLs. The PVE of these 17 QTLs ranged from 8.00 to 11.90%. These results suggest that growth traits in C. fuscus are complex and influenced by multiple genes. Additionally, the QTL intervals for single growth trait in C. fuscus are more concentrated compared to other fish species. For example, two BW-related QTLs in this study were mapped to LG09 and LG17, while in red-tail catfish (Hemibagrus wyckioides), 12 BW-related QTLs are distributed across eight LGs [40], and in crucian carp, eight BW-related QTLs span five LGs [37]. Among these 17 growth QTLs, 11 growth QTLs were distributed on LG06, LG09, and LG17. Notably, two overlapping intervals were identified: BW and BWI on LG09 (44.543 cM–45.835 cM), and BW and BL on LG17 (13.520 cM-14.832 cM). These findings indicate that, although growth traits in C. fuscus are regulated by multiple genes, the regulatory genes are mapped to a few major chromosomes. Similar observations have been made in black carp [1] and pikeperch (Sander lucioperca) [40].

4.2. Transcriptome Enrichment Pathways Associated with body Size in C. fuscus

To further explore the genes associated with the growth of C. fuscus, a comprehensive transcriptome analysis was conducted, identifying 3,824 DEGs. Among these, 2,252 genes were highly expressed, while 1,572 were lowly expressed. The functional roles and biological pathways of these DEGs were examined using GO and KEGG enrichment analyses.

In the GO enrichment analysis, the majority of DEGs were found to be enriched in biological processes (94 out of 110 pathways, 85%), a trend observed in other species as well, such as grass carp (Ctenopharyngodon idella) (42%), black porgy (Acanthopagrus schlegelii) (52%), and Zig-Zag eel (Mastacembelus armatus) (60%) [41,42,43]. These results suggest that DEGs related to biological processes are likely central to the growth rate differences observed among individuals. The most significantly enriched pathways were related to organic metabolism, particularly "pyruvate metabolism" and "ATP metabolism," consistent with the known biological roles of these genes. Pathways related to ATP precursor synthesis, including "ribonucleoside triphosphate biosynthetic process" and "ATP biosynthesis," were predominantly enriched in small-sized individuals. In contrast, big-sized individuals exhibited greater enrichment in metabolic pathways for the synthesis of various organic compounds, such as "pyruvate metabolic process" and "regulation of cellular component biosynthesis." The upregulation of the tricarboxylic acid cycle (TCA cycle) in bigger C. fuscus suggests an enhanced capacity to utilize TCA intermediates for biosynthetic pathways, thus supporting cellular function and homeostasis [44]. Conversely, smaller individuals rely on oxidative phosphorylation to rapidly generate large amounts of ATP, which is crucial for cell proliferation and tissue growth.

KEGG analysis of carbohydrate metabolism pathways revealed that genes related to glycolysis, the TCA cycle, and propionate metabolism were upregulated in bigger individuals, supporting the maintenance of cellular functions and homeostasis [45]. In contrast, smaller individuals showed a greater tendency to convert odd-chain fatty acids and amino acids into TCA cycle intermediates. This observation, consistent with previous studies [46], was supported by the increased expression of key enzymes involved in the catabolism of isoleucine and the β-oxidation of odd-chain fatty acids. These enzymes include acyl-coA oxidase 1 (acox1) [47], lipoamide acyltransferase component of branched-chain alpha-keto acid dehydrogenase (dbt) [48], and ethylmalonyl-CoA decarboxylase (echdc1) [49], which facilitate the conversion of succinyl-CoA to succinate, releasing GTP [50,51]. This observation highlights the need for rapid energy generation in smaller individuals to support growth and proliferation.

4.3. Identification and Characterization of Growth-Associated Genes in in C. fuscus

In this study, we conducted a comprehensive analysis integrating QTL mapping and transcriptome data to identify growth-associated candidate genes in C. fuscus. A total of 27 candidate genes were successfully identified across 17 QTL regions linked to growth traits. Notably, compared to small-sized individuals, several genes including eya4, serca1, f11r, npm1, and tecta were significantly upregulated in big-sized individuals. The eya4 gene, part of the EYA gene family, has been implicated in various biological processes such as organogenesis and cell differentiation [52]. Previous studies have shown that eya4 can influence the development of the zebrafish sensory system by modulating Na+/K ATPase (atp1b2b) activity [53]. Additionally, eya4 has been found to dephosphorylate polo-like kinase 1 (plk1), thereby promoting mitotic progression and enhancing cell differentiation [54]. Similarly, the gene encoding the protein serca1, a key regulator in growth promotion, has been associated with asymmetric embryonic development [55]. Mutations in serca1 disrupt left-right asymmetry during embryogenesis, leading to altered lateralization [56].Other genes such as f11r [57], npm1 [58], and tecta [59] have also been linked to cellular differentiation and organ development. The upregulation of these genes in big-sized individuals suggests their pivotal role in the growth and development of C. fuscus.

In contrast, small-sized individuals exhibited higher expression levels of genes such as mstn, egfr, mto1, ikzf2, and mks1. Mstn, a member of the transforming growth factor-β (TGF-β) superfamily, has been shown to inhibit muscle production, thereby affecting vertebrate growth [60]. Gene knockout studies have demonstrated that disrupting mstn expression can stimulate myoblast proliferation, leading to increased fish weight [61,62]. Additionally, egfr plays a critical role during early development in species such as Japanese medaka (Oryzias latipes), where it is involved in cleavage, protocole formation, and organogenesis [63]. Other studies have shown that egfr is also implicated in the maturation of zebrafish ovaries [64]. The genes mto1[65], ikzf2 [66], and mks1[67] are associated with embryonic development and organogenesis.

Our findings suggest that these genes generally promote embryonic development, cell differentiation, and proliferation, consistent with observations from enrichment analyses, which showed that many genes are enriched in biological process-related pathways. However, a key distinction emerged: genes upregulated in big-sized individuals appeared to primarily support the development of sensory systems, while those upregulated in small-sized individuals favored the development of various tissues and organs. This observation implies that energy allocation differences, stemming from varying cellular differentiation processes, may account for the observed growth and developmental discrepancies between C. fuscus individuals [68].

5. Conclusions

In this study, we conducted QTL mapping for eight growth traits in C. fuscus, identifying 17 growth-related QTLs across eight intervals, which encompass 162 functional genes. RNA-seq analysis of individuals with different sizes revealed 3,824 differentially expressed genes, with the big-sized group showing 2,252 highly expressed and 1,572 lowly expressed genes compared to the small-sized group. By integrating QTL mapping and RNA-seq data, we identified 27 candidate genes, including eya4, serca1, f11r, and npm1. This combined approach offers a more precise list of candidate genes for growth traits in C. fuscus, providing valuable genetic insights to enhance growth performance in this species.