A comprehensive linkage map and QTL map for carcass traits in a cross between Giant Grey and New Zealand White rabbits
BMC Genetics volume 16, Article number: 16 (2015)
Genomic resources for the rabbit are still limited compared to many other livestock species. The genomic sequence as well as linkage maps have gaps that hamper their use in rabbit genome research. Therefore, the aims of this study were the improvement of existing linkage maps and the mapping of quantitative trait loci (QTL) for carcass and meat quality traits. The study was performed in a F2 population of an initial cross between Giant Grey (GG) and New Zealand White (NZW) rabbits. The population consisted of 363 F2 animals derived from 9 F1 bucks and 33 F1 does. 186 microsatellite and three SNP markers were informative for mapping.
Out of 189 markers, which could be assigned to linkage groups, 110 markers were genetically mapped for the first time. The average marker distance was 7.8 cM. The map across all autosomes reached a total length of 1419 cM. The maternal linkage map was 1.4 times longer than the paternal. All linkage groups could be anchored to chromosomes. On the basis of the generated genetic map, we identified a highly significant QTL (genome-wide significance p < 0.01) for different carcass weights on chromosome 7 with a peak position at 91 cM (157 Mb), a significant QTL (p < 0.05) for bone mass on chromosome 9 at 61 cM (65 Mb), and another one for drip loss on chromosome 12 at 94 cM (128 Mb). Additional suggestive QTL were found on almost all chromosomes. Several genomic loci affecting the fore, intermediate and hind parts of the carcass were identified. The identified QTL explain between 2.5 to 14.6% of the phenotypic variance in the F2 population.
The results present the most comprehensive genetic map and the first genome-wide QTL mapping study for carcass and meat quality traits in rabbits. The identified QTL, in particular the major QTL on chromosome 7, provide starting points for fine mapping and candidate gene search. The data contribute to linking physical and genetic information in the rabbit genome.
Rabbit meat is healthy and in many countries a delicious protein source for human nutrition. As such, carcass composition and meat quality are of economic importance for rabbit breeders. The effective improvement of breeding requires the understanding of the genomic architecture and genomic information of such complex traits. Compared to other farm animal species, genomic resources for the rabbit are still limited. Although the rabbit genome has been sequenced (http://www.ensembl.org/Oryctolagus_cuniculus/, Ensembl 73, OryCun 2.0), currently only about 82% of the 2.74 Gigabase of the rabbit genome have been anchored to chromosomes. The existing microsatellite based linkage maps for the rabbit were built in two reference populations, one at INRA (France) using three rabbit INRA strains (INRA2066, Castor Orylag and Laghmere, ) and the other at the Utrecht University (Netherlands) using an F2 intercross of the rabbit strains AX/JU and IIIVO/JU [2,3]. These maps do not cover all rabbit chromosomes.
A comprehensive linkage map can help to improve the annotation and sequence assembly of the rabbit genome since it can link existing sequences of non-anchored rabbit bacterial artificial chromosomes (BACs) to the genome assembly. Such a genetic map is also an essential condition for the mapping of QTL in structured pedigrees which is an important step in quantitative trait gene identification. Even in the era of genome-wide association studies (GWAS), which are performed in unstructured populations, linkage mapping in families provides reliable genomic positions which support accurate mapping in GWAS.
This study aimed at building a comprehensive linkage map, which is anchored to the existing physical map, and using this information for mapping QTL for carcass and meat quality traits. The study was performed in an F2 pedigree of a cross between GG and NZW rabbits.
Results and discussion
Pedigree specific linkage map
Out of 387 known microsatellites [2-10], which were initially tested (Additional file 1: Table S1), 186 microsatellite and additionally three SNP markers in the myostatin (MSTN) gene , were informative for the cross between GG and NZW. The myostatin gene is located on OCU7 at 130,429,151 bp (http://www.ensembl.org/Oryctolagus_cuniculus/, Ensembl 73, OryCun 2.0). The physical position corresponds to 75.9 cM in the generated genetic map. The gene resides between the markers D7Utr3 und D7Utr4 that we used in this study. The number of observed alleles for each of the 186 microsatellites varied from two to eight with an average of 3.2 ± 1.0 in the founder animals (Additional file 2: Table S2). The heterozygosity index for all informative markers in the F1 population ranged from 0.13 to 1.00 with an average of 0.65 ± 0.22. For the number of informative meioses a variation between 54 and 775 (468.6 ± 177.6) was observed, among these 0 to 558 meioses (158.3 ± 129.9) had known phases. The mean polymorphism information content per marker in the F2 population was 0.43 ± 0.14 and varied from 0.09 (INRACCDDV0074, INRACCDDV0017) to 0.77 (INRACCDDV0293, data not shown).
The 186 microsatellite and three SNP markers could be assigned to 21 linkage groups. Twenty linkage groups are located on Oryctolagus cuniculus (OCU) chromosomes 1 to 19 and X. Two linkage groups were assigned to OCU4, but could not be linked to each other (Additional file 3: Figure S1). Out of 189 markers, 110 markers were genetically mapped for the first time. For 53 markers, we could confirm their cytogenetic positions. Compared to existing linkage groups, which were generated in other crosses, the cytogenetic positions could not be confirmed for five markers (INRACCDDV0084, INRACCDDV0218, INRACCDDV0219, INRACCDDV0230, INRACCDDV0256) in our cross. Twelve markers, which had sequence information, but had not been previously mapped, neither cytogenetically, genetically nor physically, were mapped in our population for the first time. These markers are INRACCDDV0103, INRACCDDV0165, INRACCDDV0194, INRACCDDV0302, D5L1C3, D6L2B5, D6L2H3, D6L3H10, D12L1C2, D12L1E11, D12L4A1 and OCRLADF4. These markers are particularly valuable for the mapping of BAC clones to the genome assembly to improve the rabbit genomic sequence (Additional file 2: Table S2) [1-4,6-8,10,12-17].
Among 161 markers with known physical position (http://www.ensembl.org/Oryctolagus_cuniculus/, Ensembl 73, OryCun 2.0), 155 markers had consistent positions in the genetic and physical maps. Six markers could not be mapped to their expected genomic positions. Although they were assigned to linkage groups on the expected chromosome (OCU6, 7, 8, 15 and 18) they had different positions on the chromosome. The affected markers are D6Utr4 (OCU6: 25.038137), D7L2F2 (OCU7: 60.680119), INRACCDD0323 (OCU7: 59.306464), INRACCDDV0080 (OCU8: 37.478525), INRACCDDV0143 (OCU15: 108.340756) and INRACCDDV0218 (OCU18: 68.515123). The position of the marker INRACCDDV0143 is inconsistent with respect to different published cytogenetic positions [1,4]. The map position identified in our population corresponds with the cytogenetic position 15q12 .
Seven other microsatellite markers, for which a cytogenetic position was reported [1,4,12,14], could not be linked to neighbouring markers in the expected region, instead they showed a close linkage to loci on other chromosomes in our cross. This refers to markers INRACCDDV0230 (14p11)  and INRACCDDV0219 (6p12-p13) , which were assigned to OCU13, markers INRACCDDV0218 (OCU3q14)  and INRACCDDV0256 (OCUXq12prox) , which were relocated to OCU18, and markers INRACCDDV0213 (OCU6p12prox) [1,12], INRACCDDV0127 (OCU6)  and INRACCDDV0084 (OCU20q12) , which showed an X-linked inheritance in our population. The new marker positions were confirmed by sequence alignments to the rabbit genome (http://www.ensembl.org/Oryctolagus_cuniculus/, Ensembl 73, OryCun 2.0).
The calculated order of the other markers in the corresponding linkage groups on chromosomes 1, 3, 5, 7, 8, 11, 14, 16 and 19 is consistent with previous maps [1-3], although, slight differences with regard to the distances between markers occur. This is consistent with our knowledge about pedigree specific linkage maps [18,19]. Differences in the marker order between our map and previously published maps were identified for linkage groups on chromosomes 4, 6, 9, 13, 15 and 18. Because different markers were used in different mapping populations for chromosomes 2, 10, 12, 17, 20, 21 and X, maps cannot be compared.
Since different resource populations were used to construct linkage maps, deviations in marker distances, marker orders and even positions on different chromosomes can be expected. While mapping errors cannot be excluded entirely, different mapping positions could mainly result from genome reorganizations between the different breeds that were used in the generation of resource populations for mapping the markers. The information about deviant marker positions in different populations is valuable for the genomic assembly of the rabbit genome sequence as well as for genetic and maybe even phenotypic diversity.
Maps calculated from maternal meioses across all autosomes were on average 1.4 times longer than paternal maps. Higher recombination rates in females are also consistent with findings in other species, for example in pigs , cattle , humans , and mice . However, in distinct regions on OCU4 (LG4b), OCU9 and OCU16 maternal maps were shorter than paternal maps, a result that was also observed in pigs, e. g. [18,19]. The ratios of genetic lengths between the female and male maps varied from 0.7 for OCU4 (LG 4b) and OCU16 to 5.4 for OCU11.
With a total genetic length for all autosomes of 1419 cM and an average marker distance of 7.8 cM our genetic map provides the linkage map with the highest marker coverage. Nevertheless, exceptions to good coverage still exist on chromosomes 20, 21 and Y, because no marker mapped to OCU20 and OCUY, and only one marker could be assigned to OCU21.
Phenotypic characteristics and correlation between traits
GG rabbits have about 500 g higher liveweights than NZW rabbits at the age of 84 days (Table 1). Since the total body, carcass and meat weights as well as the portions of head, fore, intermediate and hind parts are important for rabbit breeders, we have analysed all these traits. For all carcass traits, GG rabbits have higher weights compared to NZW rabbits. The F1 and F2 means of hot and reference carcass weights shifted towards the mean value of the GG breed suggesting dominance components in the mode of inheritance. The weight of the intermediate part and the meat weights of fore and intermediate parts in F1 rabbits exceeded the average performances of their parents. As expected, the liveweight and carcass weights as well as the total weight and the meat weight of the different carcass parts were highly correlated (r ≥ 0.89, p < 0.0001, data not shown); bone weights and head weight showed high correlations to the other carcass traits (r ≥ 0.62, p < 0.0001, Table 2). Similar results for highly correlation between carcass traits were observed in other studies [24-26]. Drip loss of the whole carcass showed low correlation with the pH value 24 hours p.m. at M. biceps femoris (r = 0.19, p < 0.0001, Table 2). This is in line with the correlation (r = 0.20) found between the pH value of M. biceps femoris and water holding capacity (WHC) of M. longissimus dorsi in a three-way cross . Most phenotypic correlations between carcass composition and meat quality parameters are low and partially negative (−0.28 ≤ r ≤ 0.31, p < 0.05, Table 2). Studies using principal component analysis indicated that all colour measurements, pH values and fat had low correlations [28,29]
QTL effects on carcass composition traits
The QTL analysis for carcass composition traits identified 13 genome-wide (p < 0.05 corresponding to F = 8.1) significant QTL in five genomic regions (Table 3). Additionally 55 chromosome-wise significant QTL (p < 0.05 corresponding to F > 3.6), which are considered as suggestive at the genome-wide level, were also identified (Additional file 4: Table S3).
The most significant genomic region at the genome-wide highly significance level was mapped for carcass (F-value ≥ 11.02) and meat weights (F-value = 11.49) on OCU7 with a peak position between 91 and 92 cM (157 Mb, Figure 1). The QTL peak positions are located in the distal part of the q-arm of OCU7 within the flanking interval D7L1B10 (90.8 cM, 157.32 Mb) and INRACCDDV0092 (92.4 cM, 157.49 Mb). Consistent with the high phenotypic correlation between traits, this QTL was also highly significant for hot and reference carcass weights as well as for the total weight of the intermediate part and the meat weight of the intermediate part (F-value ≥ 11.02). These QTL accounted between 6.45 to 7.45% of the respective total F2 variance. In addition, the effect was significant for the carcass and meat weights of the fore and hind parts (8.67 ≤ F-value ≤ 9.95) and kidney weight (F-value = 8.36). This region has also a suggestive effect on liveweight.
Although the F-value curve for the different traits suggests the presence of a second QTL in the linkage group (Additional file 5: Figure S2), a two QTL model did not provide statistical evidence for the presence of a second QTL on OCU7. As expected from the differences between parental rabbit breeds the GG alleles increased all carcass and meat traits. The QTL effects were additive (Figure 2, Table 3, Additional file 4: Table S3). The identified major QTL on OCU7 is probably responsible for linear growth. Since the weights of all carcass parts and meat weights are affected by this QTL in the same direction and the correlation between these traits adjusted for the QTL on OCU7 genotypes is high (r > 0.9), pleiotropic effects on the development of carcass and skeletal muscles can be assumed. This assumption is further supported by the finding that the OCU7 QTL for the weights of the carcass parts and meat weight were lost when the reference carcass weight was included as a covariate into the model.
A genome-wide significant QTL for bone weight in the fore part (F-value = 8.94) was identified on OCU9 at 61 cM (65.57 Mb, Figure 1, Table 3). The nearest markers to the peak position of the QTL for bone weight in the fore part on OCU9 were INRACCDV0010/016 (60.2 cM, 64.72 Mb) and INRACCDDV0146 (61.5 cM, 66.06 Mb). QTL alleles of Giant Grey had additive effects on the bone weight in the fore part (Figure 2, Table 3). To the same region additional suggestive effects were mapped for bone weight in the hind part (F-value = 5.58), for the fore (F-value = 6.45) and hind (F-value = 6.96) part weights of the carcass, for liveweight (F-value = 7.34), hot carcass weight (F-value = 5.81), reference carcass weight (F-value = 5.72), and the head weight (F-value = 7.22, Additional file 4: Table S3). The QTL explained between 3.29 to 7.59% of the phenotypic F2 variance of the corresponding traits (Table 3, Additional file 4: Table S3). The QTL on OCU9 affected not only bone weights, but also carcass weights and fat content in M. longissimus dorsi. This QTL particularly influenced the fore and hind parts of the carcass including total mass, bone and meat weights. When reference carcass weight was included as a covariate in the one QTL model for bone weight in the fore part, the position of the highest peak of the bone weight in the fore part shifted from 61 cM (65.57 Mb) to 102 cM (113.67 Mb, Figure 1). The direction and magnitude of the additive effects of the two QTL were consistent (Table 3, Additional file 4: Table S3).
Using the reference carcass as covariate in the model (model 2), genome-wide QTL for hind part weight were observed on OCU2 at 0 cM (29.01 Mb) and OCU19 at 45 cM (41.96 Mb) (Table 3, Figure 1). The OCU2 QTL alleles of GG had additive effects and the OCU19 GG alleles were dominant (Figure 2, Table 3). These QTL explained 5.96% and 5.07% of the phenotypic F2 variance. With the model 2, an additional genome-wide significant QTL was identified for bone weight in the fore part on OCU3 at 90 cM (132.70 Mb, Table 3, Figure 1). QTL alleles of GG had overdominance effects (Figure 2, Table 3). The QTL accounted for 6.03% of the phenotypic F2 variance.
QTL effects on meat quality traits
The QTL analysis for meat quality traits identified one genome-wide (p < 0.05) significant QTL on OCU12 (Table 3). Additionally 13 suggestive QTL at the chromosome-wise significance threshold of p < 0.05 were identified on chromosomes 1, 2, 5, 8, 9, 11, 16, 17 and 18 (Additional file 4: Table S3). The genome-wide scan for meat quality traits identified a significant QTL on OCU12 affecting drip loss of the whole carcass (F-value = 8.16, Figure 1). The peak QTL position is located at the end of the q-arm at 94 cM (127.58 Mb) near the marker INRACCDDV0176. The QTL accounted for 4.78% of the phenotypic F2 variance. The GG QTL allele had negative dominance effects (Figure 2, Table 3). Another QTL for drip loss which was suggestive was mapped on OCU18 (F-value = 5.00, Figure 1, Additional file 4: Table S3).
Since 68 QTL for carcass composition and meat quality traits are suggestive further studies are needed to confirm their effects. Therefore, these QTL are listed in Additional file 4: Table S3, but are not further discussed here.
Candidate gene identification
Previously, a single marker association analyses of the MSTN gene, as a key candidate gene affecting muscle development in different species [30-32], identified association of myostatin variants with several carcass composition traits in rabbits . In our candidate gene study, out of three SNPs in the MSTN gene, only SNP c.373 + 234G > A [GenBank: NM_001109821] showed a significant association, while the SNPs c.-125T > C and c.747 + 34C > T were not significant. The F-value curve pertaining to the linkage analysis across the whole chromosome 7, suggests the presence of a major peak for different carcass traits at the end of the chromosome in addition to the QTL effects on the same traits 6 cM away from the MSTN gene (Additional file 5: Figure S2). However, the two QTL analysis in the examined F2 population did not reach the significance level to provide evidence for the existence of a second QTL different than that identified by the single-QTL analysis at position 90–93 cM.
The estimated confidence intervals for all identified QTL effects were very large and cover almost the whole chromosome. Therefore, the selection of putative candidate genes is difficult and requires further studies to reduce the confidence intervals. Considering the confidence interval (62.0 to 98.0 cM) of the main QTL peak on OCU7, the search in the rabbit genome database (http://www.ensembl.org/Oryctolagus_cuniculus/, Ensembl release 73, OryCun 2.0) provides a list of about 300 genes. For example, the insulin like growth factor binding protein 2 (IGFBP2, 7:158.054321Mb) and the insulin like growth factor binding protein 5 (IGFBP5, 7:158.093549Mb) are located near directly under the peak position of the OCU7 at 158 Mb. These genes are positional and functional candidate genes for effect on carcass weights. IGFBP2 gene effects associated with growth and carcass composition were reported for chicken  and pigs . In addition, an overexpression of IGFBP2 reduces the postnatal body weight gain in transgenic mice . A common QTL region between sheep on OAR2 and cattle on BTA2, which is orthologous to the rabbit OCU7 QTL region, have been previously reported for carcass weight, eye muscle area and retail product yield .
This study provides a comprehensive genetic map of 189 markers for the rabbit genome. The marker linkage map as well as the link to the physical map provides valuable information for the further improvement of the rabbit genomic sequence assembly and a tool for mapping functional effects. In addition, this study was the first QTL analysis in rabbits for carcass composition and meat quality traits. The major QTL on OCU7 for carcass and meat weights, the QTL for bone and carcass weights on OCU9, as well as the QTL for drip loss on OCU12 have not been described before. This genetic information provides an important step in the identification of functional quantitative trait genes. Fine mapping in an advanced intercross population and in particular association mapping in breeding populations using dense SNP markers will facilitate candidate genes identification in the future.
For linkage analyses an F2 intercross population with 363 offspring (183 males and 180 females from 9 F1 bucks and 33 F1 does) was generated from an initial cross between six purebred GG bucks and six purebred NZW does (Additional file 6: Table S4). GG and NZW rabbits were obtained from local breeders. Rabbits were housed under standardized conditions in the experimental station of the Department of Genetics and Animal Breeding of the Agricultural University of Krakow. Adult rabbits were housed in two-storey wooden cages which were placed in a heated hall with lighting and exhaust ventilation. Cages were equipped with a water supply system (nipple drinkers). Offspring were weaned at the age of five weeks and subsequently housed in metal cages arranged in batteries with two rabbits per cage. Rabbits had ad libitum access to feed and water. The feed consisted of 16.5% protein, 14% crude fibre, and 10.2 MJ metabolisable energy. The experiment was approved by the Agricultural University of Krakow.
Phenotypes - carcass composition
The animals were slaughtered at the age of 12 weeks. After removing the skin, the head and the giblets, the weights of liver, kidney, lung, heart, head and hot carcass (without head and giblets) were recorded. Afterwards, the carcass was first kept at room temperature in a ventilated area for 45 min and then at 4°C until 24 h post mortem. After cooling, the carcass was weighted to determine the reference carcass weight. Then the carcass was cut to the fore (cut after the last rib), intermediate (cut after the last lumbar vertebra) and hind part and further dissected to meat, bone and dissectible fat. All carcass parts were weighted. The scapular, perirenal and inguinal fat percentages were calculated as percentage of the appropriated carcass part.
Phenotypes - meat quality
The pH values in the M. biceps femoris were measured at 45 min and 24 h post mortem using a pH meter with an accuracy of 0.01 (HI-9024). Meat colour was measured on the surface of M. biceps femoris according to the CIELab standards (CIE 1976: light source D65 and 8 mm diameter) at room temperature (20°C) at 45 min and 24 h post mortem with a CR-400 Minolta chromometer (Minolta Co., Ltd., Osaka, Japan). The values of lightness (L*), redness (a*) and yellowness (b*) were recorded. The shear force by Warner-Bratzler was determined in a fresh M. longissimus dorsi sample (14mm diameter, 15mm high) with the Texture Analyser TA-XT2 (Stable Micro System, Goldaming, UK) using a triangular knife incision. Drip loss was calculated as percentage of the weight difference between hot and reference carcass weight to the hot carcass weight. Protein and lipid content in M. longissimus dorsi were determined according to ISO standards. The protein content was determined by the method of Kjeldahl (PNA-04018:1975). The lipid content was determined using the method of Soxhlet (PN-ISO-1444:2000). Some phenotypes could not be measured in purebred animals of GG and NZW. These are the traits for meat colour in NZW and for shear force, protein and fat content in GG rabbits.
For genotyping, the DNA was extracted from 200 μL whole EDTA blood using NucleoSpin® Blood kit (Macherey & Nagel, Düren, Germany). Microsatellites were amplified by locus-specific PCR and fragment size was determined by the LI-COR DNA Analyzer 4200 (LI-COR Biosciences, Lincoln, USA) as described in detail, previously . Initially, we tested 387 available rabbit microsatellites [2-10] with all parental animals to identity informative markers for the cross between GG and NZW. Nine markers were fully and 180 partially informative. These markers were genotyped in all F2 animals. 122 markers were uniform and 76 markers could not be successfully amplified or did not give specific fragments (Additional file 1: Table S1). The polymorphism information content in the F2 population was calculated according to Botstein .
Since numerous mutations in the MSTN gene had been associated with growth, muscle mass, and other carcass composition traits in different species [30-32,39-43], the MSTN gene was chosen as a functional candidate gene. We additionally genotyped three SNP (c.-125T > C, c.373 + 234G > A, c.747 + 34C > T) in the rabbit MSTN gene  by allele specific PCR . To detect genotyping errors, the observed F2 genotype frequencies were compared with the expected frequencies using a chi-square test. This analysis revealed three microsatellite markers (INRACCDDV0035, INRACCDDV0157 and INRACCDDV0204) with null alleles, which were excluded from further analyses. Furthermore, we checked recombination frequencies and double recombinations between adjacent markers to detect potential genotyping errors. For checking recombination events and counting the number of informative meioses per locus we used the CHROMPIC and PREPARE options, respectively, from the software package CRI-MAP, version 2.4 .
Construction of a pedigree specific linkage map
The pedigree specific linkage map for the studied population was built on the basis of 186 microsatellite markers and three SNP markers using Kosambi mapping function in CRI-MAP software, version 2.4 . In the first step, a two-point linkage analysis was performed in which all markers were analyzed against each other. In the second step, the marker locus order was calculated with the BUILD option allowing different recombination rates in the intervals between the sexes. The BUILD option was started with the highest informative loci. Subsequently, the other loci were consecutively included for the construction of linkage groups. The FLIPS option was used to confirm the correct order of the marker loci. Finally, we generated sex-specific and sex-averaged genetic maps. The genetic distances are given in centiMorgan (cM) between markers, with the first marker of every linkage group at 0 cM. The physical positions of markers in megabase (Mb) were given according to the respective sequence position in the rabbit genome assembly at ENSEMBL (http://www.ensembl.org/Oryctolagus_cuniculus/, Ensembl 73, OryCun 2.0). Peak QTL positions were translated into physical positions as a linear genetic distance between adjacent markers with physical positions.
Basic statistical analysis
Basic statistics were performed using the PASW software package version 18.0 (SPSS, Inc., Somers, NY, USA). The phenotypic data were checked for normal distribution using the procedure EXAMINE (Kolmogorov-Smirnov test). Pearson’s correlation coefficients between traits were estimated using the CORRELATE procedure. Family (full-sibs), sex and season were detected as factors affecting the phenotypes using the GLM procedure and thus were considered as fixed effects in the QTL model. Genotype effect plots were drawn with least square means (LSM). A t-test with Bonferroni correction for multiple testing was performed to test phenotypic differences between genotype classes of the nearest marker to a QTL peak.
Single marker analyses were performed for the three MSTN SNPs. The model included common litter effects, season, sex, SNP genotype, interaction between season and family as fixed effects, and birth weight as covariate (PASW, Version 18.0). For pairwise comparisons, p-values were adjusted for multiple testing using the Bonferroni procedure .
QTL mapping was performed on the basis of the sex-averaged map using Grid-QTL . The BC-F2 module was used which assumes that founder lines are fixed for alternative alleles at QTL loci. Data were analysed with least squares regression interval mapping method using family (full-sibs, 36 levels), sex (2 levels) and season (4 levels) as fixed effects, and birth weight as an interactive covariate (model 1). Reference carcass weight was highly correlated with the weights of the carcass parts as well as the individual meat, bone and fat weights of the carcass parts (0.23 < r < 0.99, p < 0.001). Therefore, it was included as a covariate in additional analyses (model 2). Genome-wide and chromosome-wise significance thresholds were determined by permutation tests . One thousand permutations were performed for all traits. Threshold values for a given level of significance were calculated as an average of thresholds over all traits. The F-value of 10.0 corresponds to genome-wide highly significance (α = 0.01) and the F-value of 8.1 to genome-wide significance (α = 0.05). Genome-wide suggestive QTL are detected at chromosome-wise significance of α < 0.05, corresponding to F-value thresholds between 3.6 and 6.0 for the different chromosomes. The 95% confidence interval of a QTL was estimated using parametric bootstrap analysis with 1000 iterations . OCUX was analysed as a pseudo-autosome in all analyses as all markers were located in that region. The direction of the genetic effects was given as GG allele effect compared with NZW. QTL positions are given as cM distance of the highest F-value from the first marker on a chromosome. The phenotypic variance explained by a QTL was calculated as reduction of residual sum of squares in the full model (with QTL) compared with the reduced model (without QTL).
A 1 cM grid search was performed in Grid QTL by fitting model to estimate the effects of two QTL at separate positions within the same linkage group simultaneously, examining all possible pairs of markers, to test whether the two-QTL model explained significantly more variation than the best QTL from the one-QTL analysis. Two F-statistics were computed. The two-QTL model was accepted if there was a significant improvement over the best possible one-QTL model at p < 0.05 using a variance ratio (F) test.
New Zealand White
Bacterial artificial chromosomes
Genome-wide association study
Oryctolagus cuniculus chromosome
Single nucleotide polymorphism
Mega base pairs
Quantitative trait locus
Chantry-Darmon C, Urien C, de Rochambeau H, Allain D, Pena B, Hayes H, et al. A first-generation microsatellite-based integrated genetic and cytogenetic map for the European rabbit (Oryctolagus cuniculus) and localization of angora and albino. Anim Genet. 2006;37(4):335–41.
Korstanje R, Gillissen GF, den Bieman MG, Versteeg SA, van Oost B, Fox RR, et al. Mapping of rabbit chromosome 1 markers generated from a microsatellite-enriched chromosome-specific library. Anim Genet. 2001;32(5):308–12.
Korstanje R, Gillissen GF, Versteeg SA, van Oost BA, Bosma AA, Rogel-Gaillard C, et al. Mapping of rabbit microsatellite markers using chromosome-specific libraries. J Hered. 2003;94(2):161–9.
Chantry-Darmon C, Urien C, Hayes H, Bertaud M, Chadi-Taourit S, Chardon P, et al. Construction of a cytogenetically anchored microsatellite map in rabbit. Mamm Genome. 2005;16(6):442–59.
Ferrand N, Azevedo M, Mougel F. A diallelic short tandem repeat (CCCCG)4 or 5, located in intron 1 of rabbit alpha-globin gene. Anim Genet. 2000;31(1):74–5.
van Haeringen WA, den Bieman M, van Zutphen LF, van Lith HA. Polymorphic microsatellite DNA markers in the rabbit (Oryctolagus cuniculus). J Exp Anim Sci. 1996;38(2):49–57.
Korstanje R, Gillissen GF, Kodde LP, Den Bieman M, Lankhorst A, Van Zutphen LF, et al. Mapping of microsatellite loci and association of aorta atherosclerosis with LG VI markers in the rabbit. Physiol Genomics. 2001;6(1):11–8.
Mougel F, Mounolou JC, Monnerot M. Nine polymorphic microsatellite loci in the rabbit. Oryctolagus cuniculus Anim Genet. 1997;28(1):58–9.
Rico C, Rico I, Webb N, Smith S, Bell D, Hewitt G. Four polymorphic microsatellite loci for the European wild rabbit. Oryctolagus cuniculus Anim Genet. 1994;25(5):367.
Surridge AK, Bell DJ, Rico C, Hewitt GM. Polymorphic microsatellite loci in the European rabbit (Oryctolagus cuniculus) are also amplified in other lagomorph species. Anim Genet. 1997;28(4):302–5.
Sternstein I, Reissmann M, Maj D, Bieniek J, Brockmann GA. A new single nucleotide polymorphism in the rabbit (Oryctolagus cuniculus) myostatin (MSTN) gene is associated with carcass composition traits. Anim Genet. 2014;45(4):596–9.
Chantry-Darmon C, Rogel-Gaillard C, Bertaud M, Urien C, Perrocheau M, Chardon P, et al. 133 new gene localizations on the rabbit cytogenetic map. Cytogenet Genome Res. 2003;103(1–2):192–201.
Zijlstra C, de Haan NA, Korstanje R, Rogel-Gaillard C, Piumi F, van Lith HA, et al. Fourteen chromosomal localizations and an update of the cytogenetic map of the rabbit. Cytogenet Genome Res. 2002;97(3–4):191–9.
Hayes H, Rogel-Gaillard C, Zijlstra C, De Haan NA, Urien C, Bourgeaux N, et al. Establishment of an R-banded rabbit karyotype nomenclature by FISH localization of 23 chromosome-specific genes on both G- and R-banded chromosomes. Cytogenet Genome Res. 2002;98(2–3):199–205.
Chantry-Darmon C, Bertaud M, Urien C, Chadi-Taourit S, Perrocheau M, Rogel-Gaillard C, et al. Expanded comparative mapping between man and rabbit and detection of a new conserved segment between HSA22 and OCU4. Cytogenet Genome Res. 2005;111(2):134–9.
Rogel-Gaillard C, Zijlstra C, Bosma AA, Thepot D, Fontaine ML, Devinoy E, et al. Assignment of the rabbit whey acidic protein gene (WAP) to rabbit chromosome 10 by in situ hybridization and description of a large region surrounding this gene. Cytogenet Cell Genet. 2000;89(1–2):107–9.
Rogel-Gaillard C, Piumi F, Billault A, Bourgeaux N, Save JC, Urien C, et al. Construction of a rabbit bacterial artificial chromosome (BAC) library: application to the mapping of the major histocompatibility complex to position 12q.1.1. Mamm Genome. 2001;12(3):253–5.
Geldermann H, Müller E, Moser G, Reiner G, Bartenschlager H, Cepica S, et al. Genome-wide linkage and QTL mapping in porcine F2 families generated from Pietrain, Meishan and Wild Boar crosses. Genomweite Kopplungs- und QTL-Kartierung in porcinen F2-Familien, erzeugt aus Kreuzungen von Pietrain, Meishan und Wildschwein. J Anim Breed Genet. 2003;120(6):363–93.
Tortereau F, Servin B, Frantz L, Megens H-J, Milan D, Rohrer G, et al. A high density recombination map of the pig reveals a correlation between sex-specific recombination and GC content. BMC Genomics. 2012;13(1):586.
Archibald AL, Haley CS, Brown JF, Couperwhite S, McQueen HA, Nicholson D, et al. The PiGMaP consortium linkage map of the pig (Sus scrofa). Mamm Genome. 1995;6(3):157–75.
Kappes SM, Keele JW, Stone RT, McGraw RA, Sonstegard TS, Smith TP, et al. A second-generation linkage map of the bovine genome. Genome Res. 1997;7(3):235–49.
Dracopoli NC, O'Connell P, Elsner TI, Lalouel J-M, White RL, Buetow KH, et al. The CEPH consortium linkage map of human chromosome 1. Genomics. 1991;9(4):686–700.
Cox A, Ackert-Bicknell CL, Dumont BL, Ding Y, Bell JT, Brockmann GA, et al. A new standard genetic map for the laboratory mouse. Genetics. 2009;182(4):1335–44.
Elamin K. Sex effects on carcass and non carcass traits of sudanese mature belladi rabbits. Wayamba J Anim Sci. 2013;5:598–604.
Hernández P, Pla M, Blasco A. Prediction of carcass composition in the rabbit. Meat Sci. 1996;44(1–2):75–83.
Michalik D, Lewczuk A, Wilkiewicz-Wawro E, Brzozowski W. Prediction of the meat content of the carcass and valuable carcass parts in French lop rabbits using some traits measured in vivo and post mortem. Czech J Anim Sci. 2006;51(9):406–15.
Hernandez P, Pla M, Oliver MA, Blasco A. Relationships between meat quality measurements in rabbits fed with three diets of different fat type and content. Meat Sci. 2000;55(4):379–84.
Hernandez P, Pla M, Blasco A. Relationships of meat characteristics of two lines of rabbits selected for litter size and growth rate. J Anim Sci. 1997;75(11):2936–41.
Hernández P, Pla M, Blasco A. Carcass characteristics and meat quality of rabbit lines selected for different objectives: II. Relationships between meat characteristics. Livest Prod Sci. 1998;54(2):125–31.
Boman IA, Klemetsdal G, Blichfeldt T, Nafstad O, Vage DI. A frameshift mutation in the coding region of the myostatin gene (MSTN) affects carcass conformation and fatness in Norwegian White Sheep (Ovis aries). Anim Genet. 2009;40(4):418–22.
Grobet L, Martin LJ, Poncelet D, Pirottin D, Brouwers B, Riquet J, et al. A deletion in the bovine myostatin gene causes the double-muscled phenotype in cattle. Nat Genet. 1997;17(1):71–4.
McPherron AC, Lee SJ. Double muscling in cattle due to mutations in the myostatin gene. Proc Natl Acad Sci U S A. 1997;94(23):12457–61.
Lei MM, Nie QH, Peng X, Zhang DX, Zhang XQ. Single nucleotide polymorphisms of the chicken insulin-like factor binding protein 2 gene associated with chicken growth and carcass traits. Poult Sci. 2005;84(8):1191–8.
Wang W, Meng Q, Hu X, Fei J, Feng J, Liu W, et al. Chromosome location and association of haplotypes of insulin-like growth factor binding protein-2 with production performance in swine. Biochem Genet. 2008;46(7–8):381–91.
Hoeflich A, Wu M, Mohan S, Foll J, Wanke R, Froehlich T, et al. Overexpression of insulin-like growth factor-binding protein-2 in transgenic mice reduces postnatal body weight gain. Endocrinology. 1999;140(12):5488–96.
Cavanagh CR, Jonas E, Hobbs M, Thomson PC, Tammen I, Raadsma HW. Mapping Quantitative Trait Loci (QTL) in sheep. III. QTL for carcass composition traits derived from CT scans and aligned with a meta-assembly for sheep and cattle carcass QTL. Genet Sel Evol. 2010;42:36.
Bevova MR, Aulchenko YS, Aksu S, Renne U, Brockmann GA. Chromosome-wise dissection of the genome of the extremely big mouse line DU6i. Genetics. 2006;172(1):401–10.
Botstein D, White RL, Skolnick M, Davis RW. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am J Hum Genet. 1980;32(3):314–31.
Boman IA, Klemetsdal G, Nafstad O, Blichfeldt T, Vage DI. Impact of two myostatin (MSTN) mutations on weight gain and lamb carcass classification in Norwegian White Sheep (Ovis aries). Genet Sel Evol. 2010;42:4.
Clop A, Marcq F, Takeda H, Pirottin D, Tordoir X, Bibe B, et al. A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep. Nat Genet. 2006;38(7):813–8.
McPherron AC, Lawler AM, Lee SJ. Regulation of skeletal muscle mass in mice by a new TGF-beta superfamily member. Nature. 1997;387(6628):83–90.
Schuelke M, Wagner KR, Stolz LE, Hubner C, Riebel T, Komen W, et al. Myostatin mutation associated with gross muscle hypertrophy in a child. N Engl J Med. 2004;350(26):2682–8.
Stinckens A, Luyten T, Bijttebier J, Van den Maagdenberg K, Dieltiens D, Janssens S, et al. Characterization of the complete porcine MSTN gene and expression levels in pig breeds differing in muscularity. Anim Genet. 2008;39(6):586–96.
Kreuzer S, Reissmann M, Brockmann GA. Gene test to elucidate the ETEC F4ab/F4ac receptor status in pigs. Vet Microbiol. 2013;162(1):293–5.
Green P, Falls K, Crooks S. Documentation for CRIMAP, version 2.4. St. Louis, Missouri USA: Washington University School of Medicine; 1990.
Seaton G, Hernandez J, Grunchec JA, White I, Allen J, De Koning DJ, Wei W, Berry D, Haley C, Knott S: GridQTL: A Grid Portal for QTL Mapping of Compute Intensive Datasets. Belo Horizonte, Brazil: Proceedings of the 8th World Congress on Genetics Applied to Livestock Production,: 13–18 August 2006; 2006
Churchill GA, Doerge RW. Empirical threshold values for quantitative trait mapping. Genetics. 1994;138(3):963–71.
Visscher PM, Thompson R, Haley CS. Confidence intervals in QTL mapping by bootstrapping. Genetics. 1996;143(2):1013–20.
The project was supported by the German Research Foundation (DFG, project STE 1461/2-2). We acknowledge technical support from Nora Thaben for DNA isolation and microsatellite genotyping.
The authors declare that they have no competing interests.
IS designed the experiment, supervised the study, tested markers, performed genotyping and data quality control, carried out the linkage and QTL mapping analyses and drafted the manuscript, GAB contributed to write the manuscript, MR genotyped the SNP, JB generated the pedigree and supervised the phenotyping, DM phenotyped animals and collected blood. All authors read and approved the final manuscript.
Information about tested markers.
Cytogenetic associated linkage map for the F2 cross GG x NZW. The cytogenetic map (left) of every chromosome is connected to the sex averaged pedigree-specific linkage map (right). The numbers on the right hand side of the linkage maps give the estimated distances between loci in cM (Kosambi), the statistical support for the pair-wise order of markers is given on the left hand side. The total genetic length of every linkage group is given at the bottom of each bar. Connecting lines indicate the cytogenetic positions of microsatellites previously mapped by fluorescence in situ hybridization [1-4,7,12-15,17].
Positions and effects of suggestive QTL for carcass traits of the cross between GG and NZW rabbits.
F value curves across OCU7 pertaining to QTL scans for carcass traits (Model 1).
Structure of the F2 pedigree (GG x NZW).
About this article
Cite this article
Sternstein, I., Reissmann, M., Maj, D. et al. A comprehensive linkage map and QTL map for carcass traits in a cross between Giant Grey and New Zealand White rabbits. BMC Genet 16, 16 (2015). https://doi.org/10.1186/s12863-015-0168-1