Skip to main content

QTL variations for growth-related traits in eight distinct families of common carp (Cyprinus carpio)



Comparing QTL analyses of multiple pair-mating families can provide a better understanding of important allelic variations and distributions. However, most QTL mapping studies in common carp have been based on analyses of individual families. In order to improve our understanding of heredity and variation of QTLs in different families and identify important QTLs, we performed QTL analysis of growth-related traits in multiple segregating families.


We completed a genome scan for QTLs that affect body weight (BW), total length (TL), and body thickness (BT) of 522 individuals from eight full-sib families using 250 microsatellites evenly distributed across 50 chromosomes. Sib-pair and half-sib model mapping identified 165 QTLs on 30 linkage groups. Among them, 10 (genome-wide P <0.01 or P < 0.05) and 28 (chromosome-wide P < 0.01) QTLs exhibited significant evidence of linkage, while the remaining 127 exhibited a suggestive effect on the above three traits at a chromosome-wide (P < 0.05) level. Multiple QTLs obtained from different families affect BW, TL, and BT and locate at close or identical positions. It suggests that same genetic factors may control variability in these traits. Furthermore, the results of the comparative QTL analysis of multiple families showed that one QTL was common in four of the eight families, nine QTLs were detected in three of the eight families, and 26 QTLs were found common to two of the eight families. These common QTLs are valuable candidates in marker-assisted selection.


A large number of QTLs were detected in the common carp genome and associated with growth-related traits. Some of the QTLs of different growth-related traits were identified at similar chromosomal regions, suggesting a role for pleiotropy and/or tight linkage and demonstrating a common genetic basis of growth trait variations. The results have set up an example for comparing QTLs in common carp and provided insights into variations in the identified QTLs affecting body growth. Discovery of these common QTLs between families and growth-related traits represents an important step towards understanding of quantitative genetic variation in common carp.


The common carp (Cyprinus carpio) is one of the most widespread freshwater teleost species in the world. It has been domesticated as an important food fish in over 100 countries worldwide with global production exceeding 3.79 million tons in 2012, according to the Food and Agriculture Organization [1]. This globally important aquaculture species is used as a model in many research areas, e.g., ecology, physiology, and evolution. Significant progress has been made in common carp genetics and genome research over the last decade. The current genomic resources available for common carp include polymorphic genetic markers [2, 3], genetic linkage maps [46], QTLs [710], cDNA microarrays [11], bacterial artificial chromosome libraries [12, 13], and physical maps [14]. Furthermore, A draft sequence of the common carp genome has been assembled [15], which provides reference sequences for genomic and comparative genomic studies of all common carp strains and other Cyprinidae species. All of these resources, which are available online (, lay the foundation for future studies on the genetic mechanisms of economic traits in common carp and related species.

Growth is an economically important trait in the fish farming industry because it is directly related to fish production. An improvement in growth rate increases benefits for aquaculture companies because it decreases the raising time at farm facilities, leading to lower costs and higher harvests [16]. Growth-related traits such as body weight, total length, standard length, body depth, and body thickness are quantitative traits influenced by both environmental factors and multiple genes with relatively small effects according to the infinitesimal model [17]. Traditional selective breeding techniques are valuable in achieving improved fish growth [18], but the genetic gain can be increased far faster with MAS [19].

More recently, there has been an increasing number of studies that identify QTLs for growth-related traits in food fish species [2022]. Most have been carried out in salmonids (e.g., Atlantic salmon [23, 24], rainbow trout [25], and coho salmon [26]), tilapia [27], sea bass [28, 29], and turbot [30, 31] as well as in common carp, involving backcross [32] and F1 [9, 33, 34] populations. However, all of these common carp QTL studies were limited to a single segregating family and small population size (~46–190 samples). In our previous study, through comparing the distribution and variation of QTLs in four half-sib families, we found that the major QTLs for growth-related traits were not fixed either between or among families. The study also revealed that both major and minor genes differ in their genetic performance. Therefore, we concluded that the major genes are undergoing change and remain unfixed in these families [35]. Further investigation is needed to identify shared and overlapping QTLs in either populations or families to understand quantitative genetic variation.

In the present work, we conducted a genome scan for QTLs that affect body weight, total length, and body thickness in eight full-sib families containing 522 individuals from a breeding population. The objectives of this current study were: 1) to locate QTLs on linkage groups, 2) to identify which QTLs are either common or specific to all families, 3) to detect common and overlapping QTLs for the three growth-related traits, and 4) to explore the genetic architecture of growth-related traits.


Phenotypic variation

Analysis of the raw phenotypic data in the eight full-sib families revealed that all traits showed substantial levels of phenotypic variation. The average values and associated dispersion measures for BW, TL and BT are summarized in Table 1. The eight families were F234, F275, F4039, F171, F217, F373, F336 and F259, which contained 45–107 progeny. The mean value for BW, TL and BT within each family ranged from 475.9 to 722.9, from 29.8 to 36.1 and from 45.0 to 53.8, respectively. The high dispersal of means across the families resulted from the family effect caused by all of the fish being fed in the same pond. Heritabilities and correlations for the three growth-related traits are shown in Table S1 (Additional file 1). Moderate heritabilities were observed for the three traits, 0.23 for BW, 0.35 for SL and 0.25 for BT. Highly and significantly phenotypic correlations (0.83–0.92) and genetic correlations (0.93–0.98) were observed among the three traits.

Table 1 Phenotypic values of growth-related traits in the eight common carp full-sib families

QTLs for growth-related traits

QTLs that exceeded the suggestive or significant linkage threshold are revealed in Tables 2, 3, 4 and 5 and Additional file 1: Tables S2–S4. We found evidence for QTLs affecting BW on 29 LGs, TL on 28 LGs, and BT on 25 LGs. Thus, QTLs affecting all traits were found on 30 of the 50 common carp chromosomes covered by the linkage maps (Fig. 1). A total of 165 QTLs were detected across all traits by the whole genome scan. Out of the 165 QTLs, ten (six at genome-wide P <0.01 and four at genome-wide P < 0.05) and 28 (chromosome-wide P < 0.01) exhibited significant evidence of linkage, while the remaining 127 exhibited a suggestive effect on the three growth-related traits at the chromosome-wide (P < 0.05) level.

Table 2 Suggestive (P < 0.05) and significant (P < 0.01) QTLs from the sib-pair genome scan
Table 3 Significant (P < 0.01) QTLs for body weight in eight common carp families from half-sib genome scans
Table 4 Significant (P < 0.01) QTLs for total length in eight common carp families from half-sib genome scans
Table 5 Significant (P < 0.01) QTLs for body thickness in eight common carp families from half-sib genome scans
Fig. 1
figure 1

Comparison of QTLs detected among the 30 chromosomes in the eight common carp families. S and D represent the QTL based on sire- and dam half-sib analysis, respectively; combined indicates the QTL based on sib-pair analysis

For BW, the sib-pair analysis using shared markers among all eight families, identified two significant QTLs on LG24 and LG45, and six LGs containing suggestive QTLs for this trait (Table 2). A single QTL, on LG24 at 53 cM reached the genome-wide (P < 0.01) significance threshold. Additionally, the half-sib analysis revealed eight significant QTLs on six LGs (4, 11, 18, 24, 26, and 30), whereas 42 suggestive QTLs were detected on 29 LGs (Table 3, Additional file 1: Table S2). Three out of the eight significant QTLs surpassed the genome-wide significance threshold; of these, two QTLs at the genome-wide (P < 0.01) level located on LG4 in F234 segregated from the dam and LG24 in F259 segregated from the sire, and one QTL located on LG30 in F4039 that segregated from the dam reached the genome-wide (P < 0.05) significance threshold, accounting for 27.2–48.7 % of PVE. Seven QTLs were found to segregate from both sire- and dam-based analysis, these were located on LG1, LG4, and LG22 in F234, LG24, and LG30 in F4039, LG18 in F373, and LG6 in F259. The QTL found on LG30 was common to F336, F259, and F4039, and although they were located close together on the chromosome, the QTL in F336 and F259 segregated from the dam, whereas it segregated from both parents in F4039. The QTL on LG45 was common to F275, F217, and F373, and although they were positioned close together in the LG, the QTL in F217 and F373 segregated from the sire, whereas it segregated from the dam in F275.

Sib-pair QTL analysis for TL based on shared markers for all families identified significant QTLs (genome-wide P < 0.05) on two LGs (24 and 45) and suggestive QTLs (chromosome-wide P < 0.05) on five LGs (9, 13, 30, 31, and 38) (Table 2). The half-sib analysis detected 13 LGs (1, 3, 4, 8, 10, 13, 24, 30, 37, 38, 41, 43, and 45) containing 14 significant and 37 suggestive QTLs in 23 LGs (Table 4, Additional file 1: Table S3). Most of the significant QTLs were detected in only one of the families. Two significant QTLs, on LG4 in F234 that segregated from the dam and LG24 in F259 segregated from the sire reached the genome-wide (P < 0.01) level of significance, accounting for 27 % and 43.2 % of PVE, respectively. The QTL on LG45 was detected in four families (F275, F336, F373, and F4039), but only one of them that segregated from the dam in F275 was significant (chromosome-wide P < 0.01) and contributed 28.4 % of PVE. The QTL on LG1 was common to F234, F259, and F275, and although they were close to 37 cM on the LG, the QTL segregated from the dam in F259 and from the sire in F275, whereas it segregated from both parents in F234. Three families (F259, F336, and F4039) had QTLs segregating from the dam on LG30, and although they were close together on the chromosome, only one significant QTL (chromosome-wide P < 0.01) was found in F4039.

For BT, the sib-pair analysis identified significant QTLs on two LGs (13 and 24) and suggestive QTLs on four LGs (9, 31, 42, and 45) (Table 2). Only the QTL found on LG24 reached the genome-wide significance threshold (P < 0.05). This QTL was not only located on the same position at 54 cM on LG24, but also reached a genome-wide level of significance (P < 0.05) in all three traits (BW, TL, and BT). In the half-sib analysis, a total of 43 QTLs dispersed over 25 LGs were associated with BT, and from these, ten QTLs in 7 LGs (4, 9, 13, 22, 24, 30, and 44) were significant (Table 5) and 33 suggestive QTLs were found on 23 LGs (Additional file 1: Table S4). Only one QTL found on LG24 reached the genome-wide significance threshold (P < 0.01), and contributed 38.6 % of PVE. At two locations, a significant QTL was found at the same position: 17 cM on LG22 in F234 and 0 cM on LG13 in families F373 and F234. The PVE by the individual QTL varied from 13.8 to 33 %. QTLs segregating from the dams in F259, F336, and F4039 were detected at similar positions on LG30, but only one of them reached the chromosome-wide level of significance (P < 0.01) and accounted for 34.6 % of PVE in F4039. The QTL on LG45 was common to F217, F259, and F275, and all of them were suggestive, the QTL in F217 and F259 segregated from the sire whereas it segregated from both parents in F275. QTLs segregating from the sires in F217 and F259 were detected at similar positions on LG25.

Comparison of growth-related trait QTLs

Although most QTLs were either detected in just one family or at different positions on the map, 36 (24.5 %) QTLs common to all families were observed in the half-sib analysis (Fig. 1). For a single trait, one QTL common to four families, nine common to three families and 26 common to two families were identified. However, no QTL was common to more than five families. The QTL for TL located on LG45 was common to four families (F275, F336, F373, and F4039); it segregated from the dam in F275 and from the sire in F336, F373 and F4039. Similarly, the QTL for BW on LG45 was detected in three families (F275, F217, and F373) and segregated from the dam in F275 and from the sire in F217 and F373. The QTL for TL on LG45 was detected in two families (F275 and F217); it also segregated from the dam in F275 and from the sire in F217. However, five QTLs were located close to 50 cM and four QTLs harbored near120 cM on LG45 (Fig. 2). Meanwhile, we found many QTLs (16 QTLs for BW, 14 QTLs for TL and 13 QTLs for BT) that were detected only in one family. There might be many reasons for the presence of family-specific QTLs, including family-specific marker regression coefficients [36], fixation of both the QTL and markers in some families, and the linkage phase between marker and QTL alleles might vary between families.

Fig. 2
figure 2

F-distributions for QTLs on common carp LG45 affecting body weight, total length and body thickness. The F-distributions for the half-sib analyses are plotted against the primary y-axis, and the F-distribution for the sib-pair analysis is plotted against the secondary y-axis. The x-axis indicates the relative positions of the markers on the linkage map (cM). F217, F373, F275, F336, and F4039 are the family names, S and D represent sire- and dam-based half-sib analysis. sp indicates sib-pair analysis

As shown in Fig. 1, the different traits have sharing QTLs in multiple families. Approximately two LGs (3 and 45) with overlapping QTL regions associated with all traits across five families were identified; four LGs (1, 24, 25, and 30) with common QTL regions related to all traits were identified in three families; five LGs (4, 16, 18, 39, and 42) with common QTLs affecting all traits in two families; and 11 LGs with common QTLs affecting all traits in a single family. For example, QTLs influencing all traits in F373 were detected at similar positions on LG18, which segregated from the sire for BT and from both parents for BW and TL. Additionally, another QTL for BW also was found to have common regions on LG18 that segregated from the dam in F234. The QTL associated with all traits had common regions on LG1, which segregated from both parents in F234, from the sire in F275, and from the dam in F259.

Candidate gene markers for growth-related traits

The sequences of 57 loci that appeared next to QTL peaks were BLAST searched against the common carp genome sequences. The results presented that all of the sequences had significant hits on the common carp genome. Further analysis indicated that 20 sequences had very high similarity to common carp annotation genes (Additional file 1: Tables S2–4). One gene (farsb) was associated with growth, three with lipid metabolism (apobl, adiporla, and lpcat1), three with muscle development (mybcp2a, mbn12, and foxk2), and five with ionic transport and physiology (pcdh2ac, clcn3, vps13b, stim1a, and casq2).


Methodological comparison

In our study, the datasets were analyzed by two different methods. First, analyses were performed using a sib-pair model to make use of the full-sib pedigree structure and the genotypes of the parents. Second, we carried out half-sib analyses to identify segregating alleles in sires and dams with a QTL. This two-stage analytical method was used by Gutierrez et al. in their QTL investigation in Atlantic salmon [23]. Together with the data from the eight families into one analysis and each family analysis separately, we can determine which family has a specific QTL. In the sib-pair analysis, the QTL affecting all three traits on LG24 was highly significant, exceeding the 0.01 and 0.05 significance thresholds at the genome-wide level. In comparison, the half-sib analysis detected the same QTL on LG24, which influences all traits shared by three families (F234, F259, and F4039), evidence for this QTL also reached the 0.01 genome-wide significance level in F259 based on sire analysis, and surpassed either 0.01 or 0.05 chromosome-wide significance level in F234, F259, and F4039 based on sire and dam analysis. Similarly, the QTL affecting standard length (SL) identified on LG45 was detected in both the sib-pair and half-sib analyses, exceeding the 0.01 genome-wide threshold in the sib-pair analysis, 0.01 chromosome-wide threshold in F275 based on dam half-sib analysis, and 0.05 chromosome-wide in F336, F373, and F4039 based on sire half-sib analysis. Therefore, there was a better consistency between the sib-pair and half-sib analyses in their ability to identify and position this QTL.

Relationships between the QTLs identified in this study and those in the literature

Growth-related traits (e.g., body weight and length) constitute the main purpose of genetic breeding projects in aquaculture. Earlier studies on common carp have confirmed a heritable basis in these traits [37], and some QTLs affecting growth-related traits have been reported [3235, 38, 39]. For example, Laghari et al. [32] performed analysis looking for BW, TL, and condition factor (K) based on a backcross common carp family, they detected 12 QTLs associated with BW and TL on ten LGs. Laghari et al. [33] based their BW QTL analysis on phenotypic data obtained from 3-month points ranging from 10 to 12 months of culture, and they detected seven QTLs on three LGs. Zhang et al. [34] described the identification of SL QTLs using an F1 common carp family, they located four QTLs affecting SL on four LGs. However, previous QTL studies related to growth traits in common carp have generally been conducted on a single population. The ability to detect QTLs is often limited, because only the QTLs that segregate in either one or both parents can be detected [40]. The QTLs identified may not be representative of the genetic architecture of common carp growth. Therefore, the QTLs detected should be validated in other common carp families before they are used in genetic breeding. Comparing our QTL results with previously detected QTLs in populations with different genetic backgrounds, the QTL on LG22 reached 0.01 chromosome-wide significance for BT, which segregated from both parents in F234, this is similar to the QTL affecting SL found in a full-sib family [38]. Another significant QTL (chromosome-wide P < 0.01) for TL on LG13 in F234, also had an overlapping QTL region with that associated with SL reported by Zheng et al. [38]. Additionally, one QTL for BW on LG30 in F259 was identified as having a shared region with the QTL for body height (qBH30) on this chromosome in an F1 population [39]. Other than these examples, there are no other reports of overlap in the literature on QTL mapping for common carp growth traits.

Comparison of growth-related trait QTLs

Overlapping QTL regions affecting all traits were detected on 22 LGs, for example, one region on LG3, a chromosomal region was associated with three traits among five families. Overlapping QTL regions suggest the possibility that common genetic factors govern all growth traits, which is similar to the results of other researches, e.g., turbot Scophthalmus maximus (growth-related QTL) [30]. Overlapping QTL regions among multiple traits were found on nine LGs (1, 2, 5, 6, 12, 13, and 15–17). All traits were strongly correlated (Additional file 1: Table S1), further suggesting that they may be controlled by some of the common genes. While this might be caused by the wide confidence intervals for the QTL regions detected, it is also possible that this could indicate either tight QTL linkage or pleiotropic effects on growth trait variation. For example, the significant QTL on LG24 close to marker HLJ3988, affecting all three traits, is controlled by the common gene casq2 (calsequestrin 2). In some chromosome segments, QTLs identified for single traits did not overlap with QTLs for other traits, e.g., QTL for TL on LG15 and QTL for BW on LG27, suggesting that these chromosome segments contain genes that may be specifically responsible for a single trait.

The advantages of mapping QTLs in multiple families

A major limitation of QTL analysis in biparental populations is that the estimated effects are specific to that population and QTL results are often not suitable for other populations, thus restricting their use for MAS projects [41]. Therefore, mapping QTLs based on multiple segregating families allows us to identify more QTL locations, which should also allow more powerful ability to estimate the QTL effects across populations. In this study, the analysis of eight different mapping families identified more QTLs than in a single population, showed more alleles, and locate a more precise position of the QTLs that were shared by several families. Comparison of the QTL mapping results across all families showed that one QTL was common to four families, nine were detected in three families, and 26 were found in two families. The common QTLs that exist among multiple families are not only more valuable in MAS, but also probable represent relatively stable gene sites over long terms of species evolution [42]. Except for the seven common QTLs, most common QTL intervals represent existing genes, which is similar to the 60 % common QTL regions involving genes that were found in four common carp families [35]. Therefore, the authors believe that the genes that regulate traits in QTL regions are probably the basis of shared QTLs in different families. Common QTLs also showed variable phenotypic effect in terms of PVE in different families. As is shown in Fig. 3, the PVE explained by common QTLs varied. For examples, the QTL for TL on LG45 accounted for approximately 30, 20, 8.2–9.1 % of PVE in F336, F373 and F4039, and F217, F234, and F259, respectively, but less than 5 % in F171 and F275. However, only the QTL in F336, F373, and F4039 reached the significance level. The QTL for BW on LG24 contributed a high percentage, 48.7 % of the PVE in F259 and 25.2 % of PVE in F4039; however, it accounted for less than 5 % in the other six families. This might be because of the relevant loci in QTL intervals have different heterozygous states in different families, which can also be interpreted as polymorphisms at different markers in the same region contributing to genetic variation in different families. Therefore, to enhance the efficiency of MAS, it is necessary to detect QTLs that are effective across different genetic backgrounds.

Fig. 3
figure 3

Hereditary effects on the variation of QTLs based on sire-based analysis in eight families. Abscissa indicates eight families’ name, the ordinate represents phenotypic variance explained (PVE) of QTLs. BW, body weight; TL, total length; BT, body thickness

Positional candidate genes

By BLAST searching QTL loci sequences, we identified 20 markers with related genes and relevant function annotations. The functions of these genes are mainly involved in energy storage, protein synthesis, protein activity regulation, cellular components, cell differentiation, ion transport, and signal transduction, there were also some unknown functional genes. Interestingly, a particular case is that of the QTL that affects both BW and SL located on LG9 close to the gene-related marker HLJ3833, this marker is located near the farsb gene (phenylalanyl-tRNA synthetase, beta subunit), which is a candidate gene for weaning weight in Canchim beef cattle [43] as well as affecting feed efficiency in chickens [44]. The QTL associated with SL located on LG3 close to HLJ2461, which is located near clcn3 (chloride channel 3), is also a candidate gene for body size in the ninespine stickleback (Pungitius pungitius L.) [45]. We expect that focusing on these markers that are related to traits of interest will be fruitful in investigating genes. This is because the identification of QTLs influencing different traits could increase the efficiency of MAS and enhance genetic progress.


In summary, we have identified 38 significant (P < 0.01) and 127 suggestive (P < 0.05) QTLs associated with three growth traits by analyzing eight full-sib families from a breeding population. The abundant putative QTLs provide a broad view of the genetic structure of growth traits in common carp. Some of the QTLs were identified on similar chromosome segments, suggesting a role for pleiotropy and/or tight linkage in explaining the genetic basis of growth trait variation. In general, our studies lay the foundations for future comparative QTL mapping for growth-related traits and provide deeper insights into understanding and using quantitative genetic variation in common carp.


Fish and phenotypic measurements

The mapping population comprising eight full-sib families was used for QTL analysis of flesh fat content trait in common carp [7]. Briefly, the progeny came from 60 parents (40 female and 20 male) spawned in 2009 as part of the common carp breeding population, in Songpu Aquaculture Experimental Station, Heilongjiang River Fisheries Research Institute. Thirty full-sib families were constructed by crossing mature brooders with good performance, indicating high allelic diversity and genetic differences. After 50 days post hatching in a hatchery tank, ~2000 fish fry were randomly selected to feed in one pond for 2 years, we identified the genetic diversity and paternity relationships among 30 families containing 991 individuals using 25 microsatellite markers. The population consisting of 522 individuals, comprising eight full-sib families ranging from 45 to 107 progeny, was euthanized to determine the muscle fat content and record the growth-related traits. BW (g), TL (cm), and BT (mm) were measured in 522 progeny according to Part 3: measurement of characters of inspection of germplasm for cultured fishes reference standard (GB/T 18654.3-2008).

Trait correlations and heritabilities

Preliminary statistical analysis of data were completed in SPSS 13.0, and variance-covariance components were calculated using ASReml 3.0 [46]. The genetic parameters of the three traits were estimated using a multiple-trait animal model:

$$ {\mathrm{Y}}_{\mathrm{i}\mathrm{j}} = {\upmu}_{\mathrm{i}} + {\mathrm{a}}_{\mathrm{i}\mathrm{j}} + {\mathrm{e}}_{\mathrm{i}\mathrm{j}} $$

where i represents the traits, Yij is the observation of trait i for animal j, μi is the mean value for trait i, аij is the random effect of trait i for animal j and eij is the random residual error for trait i for animal j.

Heritabilities (h2 = VA/VP) were calculated using ASReml single-trait analyses, where h2 is the estimated heritability, VA is the additive genetic variance, and VP the total phenotypic variance of the trait. Genetic correlations were calculated between the traits using ASReml bivariate analysis.

Genotype data and linkage map

All individuals were genotyped using 250 microsatellite loci across regions of varying length on 50 LGs according to the common carp consensus linkage map [5]. The marker data and linkage map used in the study were described in Kuang et al. [7]. The genotype data of 233 markers, which fell into 50 common carp autosomes, were available after removal of monomorphic and unlink from a total of 250 markers. The average linkage map was used, which spanned a length of 3131.5 cM and the marker spacing per LG ranged from 6.1 to 25.8 cM with an average of 16.8 cM. Most of the LGs were moderately informative with an average information content range of 0.30–0.85. The average information content at the genome level was 0.63 [7].

QTL mapping

The QTL analysis were performed using regression techniques [47] implemented by the online software package GridQTL [48]. GridQTL is a portlet environment (available at that permits the analysis of computationally intensive datasets.

The QTL analysis was carried out using two mapping strategies. First, the analysis was conducted using a sib-pair (SP) model, which used a variance method to analyze linkage, based on alleles that are identical-by-descent (IBD). Considering the full-sib pedigree information and the absence of genotypes on the parents, the default regression method for QTL linkage analysis was applied [49]. Across-family analyses were carried out for each trait and each LG to test for evidence of QTL segregation in all families for which information was available. A trait showing evidence for a single QTL was tested for the presence of two or more QTLs by fitting a two QTL model, which was followed by within-family analyses to identify specific families segregating for putative QTLs.

Second, the half-sib (HS) model for QTL-mapping analysis was carried out, which takes advantage of the disparity in female and male recombination rates (female:male ratio was 4.2:1) in common carp [5]. By using one parent at a time, we identified the parent that was segregating for the alleles defining the QTL. To test whether or not the QTLs were segregating in the eight full-sib families, a two-stage linear regression approach was applied separately for male and female parents [50]. The QTL mapping method was according to Houston et al. [51]; in their QTL detection of Atlantic salmon, there was a similar disparity between the male and female recombination ratio. Briefly, the detection of QTLs in the initial genome scan was founded on the sire-based analysis because the low recombination rate gives greater power to detect QTLs using low resolution with few markers per LG [52]. Dam-based analyses were subsequently conducted for the LGs that showed evidence of QTL segregation in the male map, because the larger female map had a greater resolving capacity to locate putative QTLs.

F-statistics were generated at 1-cM intervals along each LG to identify the most likely QTL position. Significance thresholds were determined by a permutation test [53], with 10,000 iterations performed to obtain a P < 0.05 and P < 0.01 chromosome-wide significance levels. The genome-wide thresholds were calculated by first obtaining a Bonferroni corrected P-value at the 0.05 and 0.01 significance level (given the 50 independent linkage groups in common carp, adjusted P-value = 0.05/50 and 0.01/50) and then obtaining the genome-wide F-ratio threshold at this adjusted P-value, using 10,000 permutations [50]. F-ratio values of P < 0.01 were considered significant QTLs, whereas those in which 0.01 < P < 0.05 were considered suggestive QTLs. The chromosomes that contained significant QTLs were tested for their significance at a genome-wide level of P < 0.05. The 95 % confidence intervals were estimated using bootstrap analyses with 10,000 iterations [54]. On each LG, regions were defined based on the occurrence of QTLs.

The percentage of the phenotypic variance explained (PVE) was calculated according to Knott et al.’s method [47]. In the sire- or dam-based analysis, the formula is:

$$ \mathrm{P}\mathrm{V}\mathrm{E} = 4\left[1\hbox{-} \left({\mathrm{MSE}}_{\mathrm{full}}/{\mathrm{MSE}}_{\mathrm{reduced}}\right)\right], $$

where MSEfull and MSEreduced are the mean squared error of the full model and mean squared error of the reduced model (parameters fixed), respectively. Accordingly, the PVE was calculated from the combined sire- and dam-based analysis according to the formula:

$$ \mathrm{P}\mathrm{V}\mathrm{E} = 2\ \left(\left[1\hbox{-} {\left({\mathrm{MSE}}_{\mathrm{full}}/{\mathrm{MSE}}_{\mathrm{reduced}}\right)}^{\mathrm{Sire}}\right] + \left[1\hbox{-} {\left({\mathrm{MSE}}_{\mathrm{full}}/{\mathrm{MSE}}_{\mathrm{reduced}}\right)}^{\mathrm{Dam}}\right]\right). $$

Absolute t-value >2 for the test contrasting the effect of alternative alleles within each parent were used to identify putative heterozygous parents for the QTLs.

Candidate genes

To identify candidate genes associated with growth-related traits, we took advantage of the currently available information from the common carp genome sequencing project [15], which is publicly available at the Common Carp Genome Base (

Availability of data and materials

All marker sequences have been deposited in the NCBI database, the GenBank accession number of the sequences were shown in Zhang et al. [5]. The linkage maps used in the study were described in Kuang et al. [7]. QTL data are available within the manuscript and its additional files.

Ethics (and consent to participate)

This study was approved by the Animal Care and Use Committee of the Heilongjiang River Fisheries Research Institute at the Chinese Academy of Fishery Sciences.



Body thickness


Body weight






Linkage group


Marker-assisted selection


National Center for Biotechnology Information


Percentage of the phenotypic variance explained


Quantitative trait locus




Total length


  1. Food and Agriculture Organization. Fao Yearbook Fishery and Aquaculture Statistics 2012. Rome: Food and Agriculture Organization of the United Nations; 2014.

    Google Scholar 

  2. Ji P, Zhang Y, Li C, Zhao Z, Wang J, Li J, et al. High throughput mining and characterization of microsatellites from common carp genome. Int J Mol Sci. 2012;13:9798–807.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Xu J, Ji P, Zhao Z, Zhang Y, Feng J, Wang J, et al. Genome-wide SNP discovery from transcriptome of four common carp strains. PLoS ONE. 2012;7:e48140.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Sun X, Liu D, Zhang X, Li W, Liu H, Hong W, et al. SLAF-seq: an efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. PLoS ONE. 2013;8:58700.

    Article  Google Scholar 

  5. Zhang X, Zhang Y, Zheng X, Kuang Y, Zhao Z, Zhao L, et al. A consensus linkage map provides insights on genome character and evolution in common carp (Cyprinus carpio L.). Mar Biotechnol. 2013;15:275–312.

    Article  CAS  PubMed  Google Scholar 

  6. Zheng X, Kuang Y, Zhang X, Lu C, Cao D, Li C, et al. A genetic linkage map and comparative genome analysis of common carp (Cyprinus carpio L.) using microsatellites and SNPs. Mol Genet Genomics. 2011;286:261–77.

    Article  CAS  PubMed  Google Scholar 

  7. Kuang Y, Zheng X, Lv W, Cao D, Sun X. Mapping quantitative trait loci for flesh fat content in common carp (Cyprinus carpio). Aquaculture. 2015;435:100–5.

    Article  CAS  Google Scholar 

  8. Lu C, Zhang Y, Zheng X, Zhang X, Li C, Kuang Y, et al. Mapping QTLs of caudal fin length in common carp (Cyprinus carpio L.). New Zeal J Mar Fresh. 2015;49:96–105.

    Article  CAS  Google Scholar 

  9. Jin S, Zhang X, Jia Z, Fu H, Zheng X, Sun X. Genetic linkage mapping and genetic analysis of QTL related to eye cross and eye diameter in common carp (Cyprinus carpio L.) using microsatellites and SNPs. Aquaculture. 2012;358–359:176–82.

    Article  Google Scholar 

  10. Zhang Y, Xu P, Lu C, Kuang Y, Zhang X, Cao D, et al. Genetic linkage mapping and analysis of muscle fiber-related QTLs in common carp (Cyprinus carpio L.). Mar Biotechnol. 2011;13:376–92.

    Article  CAS  PubMed  Google Scholar 

  11. Ji P, Liu G, Xu J, Wang X, Li J, Zhao Z, et al. Characterization of common carp transcriptome: sequencing, de novo assembly, annotation and comparative genomics. PLoS ONE. 2012;7:e35152.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. LI Y, Xu P, Zhao Z, Wang J, Zhang Y, Sun XW. Construction and characterization of the BAC library for common carp Cyprinus carpio L. and establishment of microsynteny with zebrafish Danio rerio. Mar Biotechnol. 2011;13:706–12.

    Article  CAS  PubMed  Google Scholar 

  13. Xu P, Li J, LI Y, Cui R, Wang J, Zhang Y, et al. Genomic insight into the common carp (Cyprinus carpio) genome by sequencing analysis of BAC-end sequences. BMC Genomics. 2011;12:188.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Zhao L, Zhang Y, Ji P, Zhang X, Zhao Z, Hou G, et al. A dense genetic linkage map for common carp and its integration with a BAC-based physical map. PLoS ONE. 2013;8:e63928.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Xu P, Zhang X, Wang X, Li J, Liu G, Kuang Y, et al. Genome sequence and genetic diversity of the common carp, Cyprinus carpio. Nat Genet. 2014;46:1212–9.

    Article  CAS  PubMed  Google Scholar 

  16. Gjerde B. Growth and reproduction in fish and shellfish. Aquaculture. 1986;57:37–55.

    Article  Google Scholar 

  17. Mackay TF. The genetic architecture of quantitative traits. Annu Rev Genet. 2001;35:303–39.

    Article  CAS  PubMed  Google Scholar 

  18. Hulata G. Genetic manipulations in aquaculture: a review of stock improvement by classical and modern technologies. Genetica. 2001;111:155–73.

    Article  CAS  PubMed  Google Scholar 

  19. Liu ZJ, Cordes JF. DNA marker technologies and their applications in aquaculture genetics. Aquaculture. 2004;238:1–37.

    Article  CAS  Google Scholar 

  20. Gui J, Zhu Z. Molecular basis and genetic improvement of economically important traits in aquaculture animals. Chinese Sci Bull. 2012;57:1751–60.

    Article  CAS  Google Scholar 

  21. Yue GH. Recent advances of genome mapping and marker-assisted selection in aquaculture. Fish Fish. 2013;15:376–96.

    Article  Google Scholar 

  22. Tong J, Sun X. Genetic and genomic analyses for economically important traits and their applications in molecular breeding of cultured fish. Sci China Life Sci. 2015;58:178–86.

    Article  PubMed  Google Scholar 

  23. Gutierrez AP, Lubieniecki KP, Davidson EA, Lien S, Kent MP, Fukui S, et al. Genetic mapping of quantitative trait loci (QTL) for body-weight in Atlantic salmon (Salmo salar) using a 6.5K SNP array. Aquaculture. 2012;358–359:61–70.

    Article  Google Scholar 

  24. Tsai HY, Hamilton A, Guy DR, Tinch AE, Bishop SC, Houston RD. The genetic architecture of growth and fillet traits in farmed Atlantic salmon (Salmo salar). BMC Genet. 2015;16:117.

    Article  Google Scholar 

  25. Wringe BF, Devlin RH, Ferguson MM, Moghadam HK, Sakhrani D, Danzmann RG. Growth-related quantitative trait loci in domestic and wild rainbow trout (Oncorhynchus mykiss). BMC Genet. 2010;11:63.

    Article  PubMed  PubMed Central  Google Scholar 

  26. McClelland EK, Naish KA. Quantitative trait locus analysis of hatch timing, weight, length and growth rate in coho salmon, Oncorhynchus kisutch. Heredity (Edinb). 2010;105:562–73.

    Article  CAS  Google Scholar 

  27. Liu F, Sun F, Xia JH, Li J, Fu GH, Lin G, et al. A genome scan revealed significant associations of growth traits with a major QTL and GHR2 in tilapia. Sci Rep. 2014;4:7256.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Xia JH, Lin G, He X, Liu P, Liu F, Sun F, et al. Whole genome scanning and association mapping identified a significant association between growth and a SNP in the IFABP-a gene of the Asian seabass. BMC Genomics. 2013;14:295.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Wang CM, Bai ZY, He XP, Lin G, Xia JH, Sun F, et al. A high-resolution linkage map for comparative genome analysis and QTL fine mapping in Asian seabass, Lates calcarifer. BMC Genomics. 2011;12:174.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Sánchez-Molano E, Cerna A, Toro MA, Bouza C, Hermida M, Pardo BG, et al. Detection of growth-related QTL in turbot (Scophthalmus maximus). BMC Genomics. 2011;12:473.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Wang W, Hu Y, Ma Y, Xu L, Guan J, Kong J. High-density genetic linkage mapping in turbot (Scophthalmus maximus L.) based on SNP markers and major sex- and growth-related regions detection. PLoS ONE. 2015;10:e0120410.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Laghari MY, Lashari P, Zhang X, Xu P, Xin B, Zhang Y, et al. Mapping quantitative trait loci (QTL) for body weight, length and condition factor traits in backcross (BC1) family of Common carp (Cyprinus carpio L.). Mol Biol Rep. 2014;41:721–31.

    Article  CAS  PubMed  Google Scholar 

  33. Laghari MY, Zhang Y, Lashari P, Zhang X, Xu P, Xin B, et al. Quantitative trait loci (QTL) associated with growth rate trait in common carp (Cyprinus carpio). Aquacult Int. 2013;21:1373–79.

    Article  Google Scholar 

  34. Zhang Y, Wang S, Li J, Zhang X, Jiang L, Xu P, et al. Primary genome scan for complex body shape-related traits in the common carp Cyprinus carpio. J Fish Biol. 2012;82:125–40.

    Article  PubMed  Google Scholar 

  35. Zheng X, Kuang Y, Lv W, Cao D, Zhang X, Li C, et al. A consensus linkage map of common carp (Cyprinus carpio L.) to compare the distribution and variation of QTLs associated with growth traits. Sci China Life Sci. 2013;56:351–9.

    Article  PubMed  Google Scholar 

  36. Kadarmideen HN, Dekkers JC. Regression on markers with uncertain allele transmission for QTL mapping in half-sib designs. Genet Sel Evol. 1999;31:1.

    Article  Google Scholar 

  37. Vandeputte M, Kocour M, Mauger S, Dupont-Nivet M, De Guerry D, Rodina M, et al. Heritability estimates for growth-related traits using microsatellite parentage assignment in juvenile common carp (Cyprinus carpio L.). Aquaculture. 2004;235:223–36.

    Article  CAS  Google Scholar 

  38. Zheng XH, Kuang YY, Lu CY, Wang XP, Li WS, Lü WH, et al. Quantitative trait locus analysis of standard length, body depth and body thickness in mirror carp (Cyprinus carpio L.). Yichuan. 2011;33:1366–73.

    CAS  Google Scholar 

  39. Laghari MY, Lashari P, Zhang X, Xu P, Narejo NT, Liu Y, et al. Mapping QTLs for swimming ability related traits in Cyprinus carpio L. Mar Biotechnol. 2014;16:629–37.

    Article  CAS  PubMed  Google Scholar 

  40. Wang CM, Lo LC, Zhu ZY, Yue GH. A genome scan for quantitative trait loci affecting growth-related traits in an F 1 family of Asian seabass (Lates calcarifer). BMC Genomics. 2006;7:274.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Würschum T. Mapping QTL, for agronomic traits in breeding populations. Theor Appl Genet. 2012;125:201–10.

    Article  PubMed  Google Scholar 

  42. Barton NH, Keightley PD. Understanding quantitative genetic variation. Nat Rev Genet. 2002;3:11–21.

    Article  CAS  PubMed  Google Scholar 

  43. Buzanskas ME, Grossi DA, Ventura RV, Schenkel FS, Sargolzaei M, Meirelles SLC, et al. Genome-wide association for growth traits in canchim beef cattle. PLoS ONE. 2014;9:e94802.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Lee J, Karnuah AB, Rekaya R, Anthony NB, Aggrey SE. Transcriptomic analysis to elucidate the molecular mechanisms that underlie feed efficiency in meat-type chickens. Mol Genet Genomics. 2015;290:1673–82.

    Article  CAS  PubMed  Google Scholar 

  45. Laine VN, Shikano T, Herczeg G, Vilkki J, Merilä J. Quantitative trait loci for growth and body size in the nine-spined stickleback Pungitius pungitius L. Mol Ecol. 2013;22:5861–76.

    Article  PubMed  Google Scholar 

  46. Gilmour AR, Gogel BJ, Cullis BR, Thompson R. ASReml user guide release 3.0. VSN International Ltd, Hemel Hempstead, HP1 1ES, UK. 2009.

  47. Knott SA, Elsen JM, Haley CS. Methods for multiple-marker mapping of quantitative trait loci in half-sib populations. Theor Appl Genet. 1996;93:71–80.

    Article  CAS  PubMed  Google Scholar 

  48. Seaton G, Hernandez J, Grunchec JA, White I, Allen J, De Koning DJ, et al. GridQTL: a grid portal for QTL mapping of compute intensive datasets. 2006. p. 13–8.

    Google Scholar 

  49. Haseman JK, Elston RC. The investigation of linkage between a quantitative trait and a marker locus. Behav Genet. 1972;2:3–19.

    Article  CAS  PubMed  Google Scholar 

  50. Knott SA, Marklund L, Haley CS, Andersson K, Davies W, Ellegren H, et al. Multiple marker mapping of quantitative trait loci in a cross between outbred wild boar and large white pigs. Genetics. 1998;149:1069–80.

    CAS  PubMed  PubMed Central  Google Scholar 

  51. Houston RD, Haley CS, Hamilton A, Guy DR, Tinch AE, Taggart JB, et al. Major quantitative trait loci affect resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar). Genetics. 2008;178:1109–15.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Hayes BJ, Gjuvsland A, Omholt S. Power of QTL mapping experiments in commercial Atlantic salmon populations, exploiting linkage and linkage disequilibrium and effect of limited recombination in males. Heredity (Edinb). 2006;97:19–26.

    Article  CAS  Google Scholar 

  53. Churchill GA, Doerge RW. Empirical threshold values for quantitative trait mapping. Genetics. 1994;138:963–71.

    CAS  PubMed  PubMed Central  Google Scholar 

  54. Visscher PM, Thompson R, Haley CS. Confidence intervals in QTL mapping by bootstrapping. Genetics. 1996;143:1013–20.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references


Authors thank Dr. Yi Zhou (Boston Children’s Hospital of Harvard Medical School, USA) for English editing. We also appreciate the reviewers and editors for their patience to the work.


This work was funded by the National High Technology Research and Development Program of China (No. 2011AA100402), the Central-Level Non-profit Scientific Research Institutes Special Funds (2015B03XK01), and the National Infrastructure of Fishery Germplasm Resources.

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Yunqin Yan or Xiaowen Sun.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

XS and YY conceived and designed the experiments; WL performed the experiments; XZ and YK analyzed the data; WL and DC performed the measurements; WL and XZ wrote the manuscript; YK, YY, and XS revised the manuscript. All authors have read and approved the final version of the manuscript.

Additional file

Additional file 1: Table S1.

Phenotypic correlations between growth-related traits in eight common. Table S2. List of suggestive (P < 0.05) and significant (P < 0.01) QTLs for body weight (BW) in eight common carp families based on half-sib analysis. Table S3. List of suggestive (P < 0.05) and significant (P < 0.01) QTLs for total length (TL) in eight common carp families based on half-sib analysis. Table S4. List of suggestive (P < 0.05) and significant (P < 0.01) QTLs for body thickness (BT) in eight common carp families based on half-sib analysis. (DOCX 123 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lv, W., Zheng, X., Kuang, Y. et al. QTL variations for growth-related traits in eight distinct families of common carp (Cyprinus carpio). BMC Genet 17, 65 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: