Mapping QTL for production traits in segregating Piétrain pig populations using genome‐wide association study results of F2 crosses

In this study, genome‐wide association study (GWAS) results of porcine F2 crosses were used to map QTL in outcross Piétrain populations. For this purpose, two F2 crosses (Piétrain × Meishan, n = 304; Piétrain × Wild Boar, n = 291) were genotyped with the PorcineSNP60v2 BeadChip and phenotyped for the dressing yield, carcass length, daily gain and drip loss traits. GWASs were conducted in the pooled F2 cross applying single marker mixed linear models. For the investigated traits, between two and five (in total 15) QTL core regions, spanning 250 segregating SNPs around a significant trait‐associated peak SNP, were identified. The SNPs within the QTL core regions were subsequently tested for trait association in two outcross Piétrain populations consisting of 771 progeny‐tested boars and 210 sows with their own performance records. In the sow (boar) dataset, five (eight) of the 15 mapped QTL were validated. Hence, many QTL mapped in the F2 crosses (with Piétrain as a common founder breed) are still segregating in the current Piétrain breed. This confirms the usefulness of existing F2 crosses for mapping QTL that are still segregating in the recent founder breed generation. The approach utilizes the high power of an F2 cross to map QTL in a breeding population for which it is not guaranteed that they would be found using a GWAS in this population.     

Linkage disequilibrium pattern and genome‐wide association mapping for meat traits in multiple porcine F2 crosses

In the present study, data from four F2 crosses were analysed and used to study the linkage disequilibrium (LD) structure within and across the crosses. Genome‐wide association analyses (GWASes) for conductivity and dressing out meat traits were conducted using single‐marker and Bayesian multi‐marker models using the pooled data from all F2 crosses. Porcine F2 crosses generated from the distantly related founder breeds Wild Boar, Piétrain and Meishan, as well as from a porcine F2 cross from the closely related founder breed Piétrain and an F1 Large White × Landrace cross were pooled. A total of 2572 F2 animals were genotyped using a 62K SNP chip. The positions of the SNPs were based on genome assembly Sscrofa11.1. After post‐alignment and genotype filtering, approximately 50K SNPs were usable for LD studies and GWASes. The main findings of the present study are that the breakdown of LD was faster in crosses from closely related founder breeds compared to crosses from distantly related founders. The fastest breakdown of LD was observed by pooling the data. Based on the single‐marker results and LD structure, clusters and windows were built for 1‐Mb intervals. For conductivity and dressing out, 183 and 191 nominal significant associations respectively and six and five clusters respectively were found. Dominance was important for conductivity, and considering dominance in GWASes improved the mapping signals. Most clear signals were found for conductivity on SSC6, 8 and 15 and for dressing out on SSC2 and 7. Considering dominance might contribute to the accuracy of genomic selection and serve as a guide for choosing mating pairs with good combining abilities. However, further research is needed to investigate if dominance is also important in crossbreed pig breeding schemes.     

A genome‐wide association study for a proxy of intermuscular fat level in the Italian Large White breed identifies genomic regions affecting an important quality parameter for dry‐cured hams

Intermuscular fat content in protected designations of origin dry‐cured hams is a very important meat quality trait that affects the acceptability of the product by the consumers. An excess in intermuscular fat (defined as the level of fat deposition between leg muscles) is a defect that depreciates the final product. In this study we carried out a genome‐wide association study for visible intermuscular fat (VIF) of hams in the Italian Large White pig breed. This trait was evaluated on the exposed muscles of green legs in 1122 performance‐tested gilts by trained personnel, according to a classification scale useful for routine and cheap evaluation. All animals were genotyped with the Illumina PorcineSNP60 BeadChip. The genome‐wide association study identified three QTL regions on porcine chromosome 1 (SSC1; accounting for ~79% of the SNPs below the 5.0E?04 threshold) and SSC2, two on SSC7 and one each on SSC3, SSC6, SSC9, SSC11, SSC13, SSC15, SSC16 and SSC17. The most significant SNP (ALGA0004143 on SSC1 at 77.3 Mb; P_FDR < 0.05), included in the largest QTL region which spanned about 6.8 Mb on SSC1, is located within the glutamate ionotropic receptor kainate type subunit 2 (GRIK2) gene. Functional annotation of all genes included in QTL regions for VIF suggested that intermuscular fat in the Italian Large White breed is a complex trait apparently influenced by complex biological mechanisms also involving obesity‐related processes. These QTL target mainly chromosome regions different from those affecting subcutaneous and intramuscular fat deposition.     

SNP genotypes reveal breed substructure,selection signatures and highly inbred regions in Piétrain pigs

The Piétrain pig originates from the Belgian village Piétrain some time between 1920 and 1950. Owing to its superior conformation, the Piétrain has spread worldwide since the 1960s. As initial population sizes were limited and close inbreeding was commonplace, the breed’s genetic diversity has been questioned. Therefore, this study examines Piétrain breed substructure, diversity and selection signatures using SNP data in comparison with Duroc, Landrace and Large White populations. Principal component analysis indicated three subpopulations, and F_ST analysis showed that US Piétrains differ most from European Piétrains. Average inbreeding based on runs of homozygosity (ROH) segments larger than 4 Mb ranged between 16.7 and 20.9%. The highest chromosomal inbreeding levels were found on SSC8 (42.7%). ROH islands were found on SSC8, SSC15 and SSC18 in all Piétrain populations, but numerous population-specific ROH islands were also detected. Moreover, a large ROH island on SSC8 (34–126 Mb) appears nearly fixed in all Piétrain populations, with a unique genotype. Chromosomal ROH patterns were similar between Piétrain populations. This study shows that Piétrain populations are genetically diverging, with at least three genetically distinct populations worldwide. Increasing genetic diversity in local Piétrain populations by introgression from other Piétrain populations seems to be only limited. Moreover, a unique 90 Mb region on SSC8 appeared largely fixed in the Piétrain breed, indicating that fixation was already present before the 1960s. We believe that strong selection and inbreeding during breed formation fixed these genomic regions in Piétrains. Finally, we hypothesize that independent coat color selection may have led to large ROH pattern similarities on SSC8 between unrelated pig breeds.     

Molecular Characterization of Five Porcine Candidate Genes for Drip Loss in Pork

Drip loss is the loss of fluid from a piece of meat without mechanical force and represents an important meat quality trait. Previous work revealed a quantitative trait locus (QTL) for drip loss in pork in an experimental Duroc x Piétrain (DUPI) F2 family on SSC 5. Based on functional data indicating their possible involvement in water holding capacity and their expression in skeletal muscle, we selected five positional candidates (ACO2, ADSL, CBY1, KCNJ4, PLA2AG6) out of 130 predicted genes in the QTL interval for further analysis. We performed a mutation analysis of all coding exons and discovered 204 polymorphisms. We genotyped 39 single nucleotide polymorphisms (SNPs) in 192 Piétrain pigs with extreme drip loss phenotypes and detected a possible association with drip loss for one non-coding SNP in the ADSL gene (ss107793818, p_raw = 0.021). Correspondingly, ADSL diplotypes were associated with drip loss and pH1 of M. longissimus dorsi. However, after correction for multiple testing, none of the tested SNPs were significantly associated with drip loss. One possible explanation for these results is that one of the QTL-alleles from the experimental DUPI family may be fixed or nearly fixed in the tested Piétrain population.     

Genome‐wide association study identifies quantitative trait loci affecting hematological traits in an F2 intercross between Landrace and Korean native pigs

Changes affecting the status of health and robustness can bring about physiological alterations including hematological parameters in swine. To identify quantitative trait loci (QTL) associated with eight hematological traits (one leukocyte trait, six erythrocyte traits and one platelet trait), we conducted a genome‐wide association study using the PorcineSNP60K BeadChip in a resource population derived from an intercross between Landrace and Korean native pigs. A total of 36 740 SNPs from 816 F₂ progeny were analyzed for each blood‐related trait after filtering for quality control. Data were analyzed by the genome‐wide rapid association using mixed model and regression (GRAMMAR) approach. A total of 257 significant SNPs (P < 1.36 × 10^?6) on SSC3, 6, 8, 13 and 17 were identified for blood‐related traits in this study. Interestingly, the genomic region between 17.9 and 130 Mb on SSC8 was found to be significantly associated with red blood cell, mean corpuscular volume and mean corpuscular hemoglobin. Our results include the identification of five significant SNPs within five candidate genes (KIT, IL15, TXK, ARAP2 and ERG) for hematopoiesis. Further validation of these identified SNPs could give valuable information for understanding the variation of hematological traits in pigs.     

Genome‐wide association study for ham weight loss at first salting in Italian Large White pigs: towards the genetic dissection of a key trait for dry‐cured ham production

Protected designation of origin dry‐cured hams are the most important productions of the Italian heavy pig industry. Hams capable of minimal seasoning losses produce better quality dry‐cured hams. Ham weight loss during the first 7 days in brine (first salting) is highly correlated with the total loss of weight up to the end of seasoning, and it has quite high heritability (0.30–0.61). For these reasons, ham weight loss at first salting has been included as a meat quality trait in the Italian heavy pig selection program. In this work, we carried out a genome‐wide association study for this parameter in the Italian Large White pig breed by genotyping 1365 animals with the Illumina BeadChip PorcineSNP60 chip. A total of 44 single nucleotide polymorphisms (SNPs) had a P_{nominal value} below 5.0E‐04, five of which were below 5.0E‐05 and one of them (ALGA0057985 on chromosome 10) was associated with this trait at a P_Bonferroni threshold of 0.10. These SNPs identified a total of at least 29 putative QTLs that were located on most porcine autosomal chromosomes. This study provides genomic information that could be useful in dissecting this complex trait by identifying potential candidate genes whose function could contribute to understanding the biological mechanisms affecting meat quality for seasoning aptitude.     

The effect of Piétrain sire on the performance of the progeny of two commercial dam breeds: a pig intervention study

Genetic evaluation of Piétrain sires in Flanders occurs under standardized conditions, on test stations with fixed dam breeds, standardized diets and uniform management practices. As environmental conditions vary on commercial farms and differ from the test stations, this study aimed at understanding to what extent the sire, the dam breed and the interaction between both affects the translation of breeding values to practice. Dams of two commercial breeds were inseminated with semen from one of five different sires selected for contrasting breeding values (daily gain, feed conversion ratio and carcass quality). For each sire by dam breed combination, six pen replicates (with three gilts and three barrows per pen) were evaluated for growth performance from 9 weeks of age (20 kg) to slaughter (110 kg), and for carcass and meat quality. In our experimental setup, both sire and dam breed affected growth, carcass and meat quality traits. No significant sire×dam breed interactions on performance could be detected. Though a tendency for interaction on average daily feed intake between 20 and 110 kg (P=0.087), and on pork colour (lightness) (P=0.093) was present. In general, offspring of all tested sires behaved similarly in both dam breeds, indicating that estimated breeding values for Piétrain sires determined in one dam breed are representative in other dam breeds as well.     

A genome‐wide association analysis for carcass traits in a commercial Duroc pig population

We performed a genome‐wide association study to map the genetic determinants of carcass traits in 350 Duroc pigs typed with the Porcine SNP60 BeadChip. Association analyses were carried out using the gemma software. The proportion of phenotypic variance explained by the SNPs ranged between negligible to moderate (= 0.01–0.30) depending on the trait under consideration. At the genome‐wide level, we detected one significant association between backfat thickness between the 3^rd and 4^th ribs and six SNPs mapping to SSC12 (37–40 Mb). We also identified several chromosome‐wide significant associations for ham weight (SSC11: 51–53 Mb, three SNPs; 67–68 Mb, two SNPs), carcass weight (SSC11: 66–68 Mb, two SNPs), backfat thickness between the 3^rd and 4^th ribs (SSC12: 21 Mb, one SNP; 33–40 Mb, 17 SNPs; 51–58 Mb, two SNPs), backfat thickness in the last rib (SSC12: 37 Mb, one SNP; 40–41 Mb, nine SNPs) and lean meat content (SSC13: 34 Mb, three SNPs and SSC16: 45.1 Mb, one SNP; 62–63 Mb, 10 SNPs; 71–75 Mb, nine SNPs). The ham weight trait‐associated region on SSC11 contains two genes (UCHL3 and LMO7) related to muscle development. In addition, the ACACA gene, which encodes an enzyme for the catalysis of fatty acid synthesis, maps to the SSC12 (37–41 Mb) region harbouring trait‐associated regions for backfat thickness traits. Sequencing of these candidate genes may help to uncover the causal mutations responsible for the associations found in the present study.     

Quantitative trait loci associated with AutoFOM grading characteristics, carcass cuts and chemical body composition during growth of Sus scrofa

A three-generation full-sib resource family was constructed by crossing two commercial pig lines. Genotypes for 37 molecular markers covering chromosomes SSC1, SSC6, SSC7 and SSC13 were obtained for 315 F2 animals of 49 families and their parents and grandparents. Phenotypic records of traits including carcass characteristics measured by the AutoFOM grading system, dissected carcass cuts and meat quality characteristics were recorded at 140 kg slaughter weight. Furthermore, phenotypic records on live animals were obtained for chemical composition of the empty body, protein and lipid accretion (determined by the deuterium dilution technique), daily gain and feed intake during the course of growth from 30 to 140 kg body weight. Quantitative trait loci (QTL) detection was conducted using least-squares regression interval mapping. Highest significance at the 0.1% chromosome-wise level was obtained for five QTL: AutoFOM belly weight on SSC1; ham lean-meat weight, percentage of fat of primal cuts and daily feed intake between 60 and 90 kg live weight on SSC6; and loin lean-meat weight on SSC13. QTL affecting daily gain and protein accretion were found on SSC1 in the same region. QTL for protein and lipid content of empty body at 60 kg liveweight were located close to the ryanodine receptor 1 (RYR1) locus on SSC6. On SSC13, significant QTL for protein accretion and feed conversion ratio were detected during growth from 60 to 90 kg. In general, additive genetic effects of alleles originating from the Piétrain line were associated with lower fatness and larger muscularity as well as lower daily gain and lower protein accretion rates. Most of the QTL for carcass characteristics were found on SSC6 and were estimated after adjustment for the RYR1 gene. QTL for carcass traits, fatness and growth on SSC7 reported in the literature, mainly detected in crosses of commercial lines x obese breeds, were not obtained in the present study using crosses of only commercial lines, suggesting that these QTL are not segregating in the analysed commercial lines.     

Genomic regions affecting backfat thickness and cannon bone circumference identified by genome‐wide association study in a Duroc pig population

We performed a genome‐wide association study using the porcine 60K SNP array to detect QTL regions for nine traits in a three‐generational Duroc samples (n = 651), viz. generations 1, 2 and 3 from a population selected over five generations using a closed nucleus breeding scheme. We applied a linear mixed model for association mapping to detect SNP effects, adjusting for fixed effects (sex and season) and random polygenic effects (reflecting genetic relatedness), and derived a likelihood ratio statistic for each SNP using the efficient mixed‐model association method. We detected a region on SSC6 for backfat thickness (BFT) and on SSC7 for cannon bone circumference (CANNON), with a genome‐wide significance of P < 0.01 after Bonferroni correction. These regions had been detected previously in other pig populations. Six genes are located in the BFT‐associated region, while the CANNON‐associated region includes 66 genes. In the future, significantly associated SNPs, derived by sequencing the coding regions of the six genes in the BFT region, can be used in marker‐assisted selection of BFT, whereas haplotypes constructed from the SSC7 region with strong LD can be used to select for the CANNON trait in our resource family.     

A genome‐wide association study to detect genetic variation for postpartum dysgalactia syndrome in five commercial pig breeding lines

Postpartum dysgalactia syndrome (PDS) in sows is an important disease after parturition with a relevant economic impact, affecting the health and welfare of both sows and piglets. The genetic background of this disease has been discussed and its heritability estimated, but further genetic analyses are lacking in detail. The aim of the current study was to detect loci affecting the susceptibility to PDS through a genome‐wide association approach. The study was designed as a family‐based association study with matched sampling of affected sows and healthy half‐ or full‐sib control sows on six farms. For the study, 597 sows (322 affected vs. 275 healthy control sows) were genotyped on 62 163 single nucleotide polymorphisms (SNPs) using the Illumina PorcineSNP60 BeadChip. After quality control, 585 sows (314 affected vs. 271 healthy control sows) and 49 740 SNPs remained for further analysis. Statistics were performed mainly with the r package genabel and included a principal component analysis. A statistically significant genome‐wide associated SNP was identified on porcine chromosome (SSC) 17. Further promising results with moderate significance were detected on SSC 13 and on an unplaced scaffold with an older annotation on SSC 15. The PRICKLE2 and NRP2 genes were identified as candidate genes near associated SNPs. Several quantitative trait loci (QTL) have been previously described in these genomic regions, including QTL for mammary gland condition, as teat number and non‐functional nipples QTL, as well as QTL for body temperature and gestation length.     

A Genome-Wide Association Study to Detect QTL for Commercially Important Traits in Swiss Large White Boars

The improvement of meat quality and production traits has high priority in the pork industry. Many of these traits show a low to moderate heritability and are difficult and expensive to measure. Their improvement by targeted breeding programs is challenging and requires knowledge of the genetic and molecular background. For this study we genotyped 192 artificial insemination boars of a commercial line derived from the Swiss Large White breed using the PorcineSNP60 BeadChip with 62,163 evenly spaced SNPs across the pig genome. We obtained 26 estimated breeding values (EBVs) for various traits including exterior, meat quality, reproduction, and production. The subsequent genome-wide association analysis allowed us to identify four QTL with suggestive significance for three of these traits (p-values ranging from 4.99×10⁻⁶ to 2.73×10⁻⁵). Single QTL for the EBVs pH one hour post mortem (pH1) and carcass length were on pig chromosome (SSC) 14 and SSC 2, respectively. Two QTL for the EBV rear view hind legs were on SSC 10 and SSC 16.     

QTL analysis of body weight and carcass body length traits in an F2 intercross between Landrace and Korean native pigs

Growth traits, such as body weight and carcass body length, directly affect productivity and economic efficiency in the livestock industry. We performed a genome‐wide linkage analysis to detect the quantitative trait loci (QTL) that affect body weight, growth curve parameters and carcass body length in an F₂ intercross between Landrace and Korean native pigs. Eight phenotypes related to growth were measured in approximately 1000 F₂ progeny. All experimental animals were subjected to genotypic analysis using 173 microsatellite markers located throughout the pig genome. The least squares regression approach was used to conduct the QTL analysis. For body weight traits, we mapped 16 genome‐wide significant QTL on SSC1, 3, 5, 6, 8, 9 and 12 as well as 22 suggestive QTL on SSC2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 16 and 17. On SSC12, we identified a major QTL affecting body weight at 140 days of age that accounted for 4.3% of the phenotypic variance, which was the highest test statistic (F‐ratio = 45.6 under the additive model, nominal P = 2.4 × 10^?11) observed in this study. We also showed that there were significant QTL on SSC2, 5, 7, 8, 9 and 12 affecting carcass body length and growth curve parameters. Interestingly, the QTL on SSC2, 3, 5, 6, 8, 9, 10, 12 and 17 influencing the growth‐related traits showed an obvious trend for co‐localization. In conclusion, the identified QTL may play an important role in investigating the genetic structure underlying the phenotypic variation of growth in pigs.     

Four loci differentially expressed in muscle tissue depending on water-holding capacity are associated with meat quality in commercial pig herds

Four genes, VTN, KERA, LYZ, and a non-annotated EST (Affymetrix probe set ID: Ssc.25503.1.S1_at), whose candidacy for traits related to water-holding capacity of meat arises from their trait-dependent differential expression, were selected for candidate gene analysis. Based on in silico analysis SNPs were detected, confirmed by sequencing and used to genotype animals of 4 pig populations including 3 commercial herds of Pietrain (PI), Pietrain × (German Large White × German Landrace) (PIF1), German Landrace (DL) and 1 experimental F₂ population Duroc × Pietrain (DUPI). Comparative and genetic mapping established the location of VTN on SSC12, of LYZ and KERA on SSC5 and of UN on SSC7, coinciding with QTL regions for meat quality traits. VTN showed association with pH1, pH24 and drip loss. LYZ revealed association with conductivity 24, pH1 and drip loss. KERA was associated with pH. UN showed association with pH24 and drip loss, respectively. However, none of the candidate genes showed significant associations for a particular trait across all populations. This may be due to breed specific effects that are related to the differences in meat quality of theses pig breeds. The studies revealed statistic evidence for a link of genetic variation at these loci or close to them and promoted those four candidate genes as functional and/or positional candidate genes for meat quality traits.     

Genome‐wide association analysis identifies quantitative trait loci for growth in a Landrace purebred population

Growth‐related traits are complex and economically important in the livestock industry. The aim of this study was to identify quantitative trait loci (QTL) and the associated positional candidate genes affecting growth in pigs. A genome‐wide association study (GWAS) was performed using the porcine single‐nucleotide polymorphism (SNP) 60K bead chip. A mixed‐effects model and linear regression approach were used for the GWAS. The data used in the study included 490 purebred Landrace pigs. All experimental animals were genotyped with 39 438 SNPs located throughout the pig autosomes. We identified a strong association between a SNP marker on chromosome 16 and body weight at 71 days of age (ALGA0092396, P = 5.35 × 10^?9, Bonferroni adjusted P < 0.05). The SNP marker was located near the genomic region containing IRX4, which encodes iroquois homeobox 4. This SNP marker could be useful in the selective breeding program after validating its effect on other populations.     

Twenty years of artificial directional selection have shaped the genome of the Italian Large White pig breed

In this study, we investigated at the genome‐wide level if 20 years of artificial directional selection based on boar genetic evaluation obtained with a classical BLUP animal model shaped the genome of the Italian Large White pig breed. The most influential boars of this breed (n = 192), born from 1992 (the beginning of the selection program of this breed) to 2012, with an estimated breeding value reliability of >0.85, were genotyped with the Illumina Porcine SNP60 BeadChip. After grouping the boars in eight classes according to their year of birth, filtered single nucleotide polymorphisms (SNPs) were used to evaluate the effects of time on genotype frequency changes using multinomial logistic regression models. Of these markers, 493 had a P_Bonferroni < 0.10. However, there was an increasing number of SNPs with a decreasing level of allele frequency changes over time, representing a continuous profile across the genome. The largest proportion of the 493 SNPs was on porcine chromosome (SSC) 7, SSC2, SSC8 and SSC18 for a total of 204 haploblocks. Functional annotations of genomic regions, including the 493 shifted SNPs, reported a few Gene Ontology terms that might underly the biological processes that contributed to increase performances of the pigs over the 20 years of the selection program. The obtained results indicated that the genome of the Italian Large White pigs was shaped by a directional selection program derived by the application of methodologies assuming the infinitesimal model that captured a continuous trend of allele frequency changes in the boar population.     

Genome‐wide association study for the level of serum electrolytes in Italian Large White pigs

Calcium, magnesium and phosphorus are essential electrolytes involved in a large number of biological processes. Imbalance of these minerals in blood may indicate clinically relevant conditions and are important in inferring acute or chronic pathologies in humans and animals. In this work, we carried out a genome‐wide association study (GWAS) for the level of these three electrolytes in the serum of 843 performance‐tested Italian Large White pigs. All pigs were genotyped with the Illumina PorcineSNP60 BeadChip, and GWAS was carried out using genome‐wide efficient mixed‐model association. For the level of Ca²⁺, eight single nucleotide polymorphisms (SNPs) were significant, considering a false discovery rate (FDR) < 0.05, and another eight were above the moderate association threshold (P_{nominal value} < 5.00E‐05). These SNPs are distributed in four porcine chromosomes (SSC): SSC8, SSC11, SSC12 and SSC13. In particular, a few putative different signals of association detected on SSC13 and one on SSC12 were in genes or close to genes involved in calcium metabolism (P2RY1, RAP2B, SLC9A9, C3orf58, TSC22D2, PLCH1 and CACNB1). Only one SNP (on SSC7) and six SNPs (on SSC2 and SSC7) showed moderate association with the level of magnesium and phosphorus respectively. The association signals for these two latter minerals might identify genes not known thus far for playing a role in their biological functions and regulations. In conclusion, our GWAS contributed to increased knowledge on the role that calcium, magnesium and phosphorus may play in the genetically determined physiological mechanisms affecting the natural variability of mineral levels in mammalian blood.     

Genetic mapping of quantitative trait loci for meat quality and muscle metabolic traits in cattle

A whole‐genome scan was carried out in New Zealand and Australia to detect quantitative trait loci (QTL) for live animal and carcass composition traits and meat quality attributes in cattle. Backcross calves (385 heifers and 398 steers) were generated, with Jersey and Limousin backgrounds. The New Zealand cattle were reared and finished on pasture, whilst Australian cattle were reared on grass and finished on grain for at least 180 days. This paper reports on meat quality traits (tenderness measured as shear force at 4–5 ages on two muscles as well as associated traits of meat colour, pH and cooking loss) and a number of metabolic traits. For meat quality traits, 18 significant QTL (P < 0.05), located in nine linkage groups, were detected on a genome‐wise basis, in combined‐sire (seven QTL) or within‐sire analyses (11 QTL). For metabolic traits, 11 significant QTL (P < 0.05), located in eight linkage groups, were detected on a genome‐wise basis, in combined‐sire (five QTL) or within‐sire analyses (six QTL). BTA2 and BTA3 had QTL for both metabolic traits and meat quality traits. Six significant QTL for meat quality and metabolic traits were found at the proximal end of chromosome 2. BTA2 and BTA29 were the most common chromosomes harbouring QTL for meat quality traits; QTL for improved tenderness were associated with Limousin‐derived and Jersey‐derived alleles on these two chromosomes, respectively.