Red and blond Mangalitza pigs display a signature of divergent directional selection in the SLC45A2 gene.

The Mangalitza lard-type pig breed is well known for its fat appearance and curly hair, and it is mainly distributed in Eastern Europe. Four main lines were created in the nineteenth century by artificial selection: Blond Mangalitza, Black Mangalitza, Swallow-Belly Mangalitza and Red Mangalitza. The Swallow-Belly line has a black coat combined with yellow-blond throat and underbelly. In the current work, we aimed to investigate if the colourations of Mangalitza pigs are genetically determined by one or a few loci whose frequencies have been modified by artificial selection. The results of selection scans, with HapFLK and BayeScan, and of a GWAS for coat colour highlighted the existence of one region on SSC16 (18-20 Mb) with potential effects on hair pigmentation (Red vs. Blond contrast). The analysis of the gene content of this region allowed us to detect the solute carrier family 45 member 2 (SLC45A2) locus as a candidate gene for this trait. The polymorphism of the SLC45A2 locus has been associated with reduced levels or the absence of melanin in several mammalian species. The genotyping of four missense polymorphisms evidenced that rs341599992:G > A and rs693695020:G > A SNPs are strongly but not fully associated with the red and blond coat colours of Mangalitza pigs, a result that was confirmed by performing a haplotype association test. The near fixation of alternative SLC45A2 genotypes in Red and Blond Mangalitza pigs provides a compelling example of the consequences of a divergent directional selection for coat colour in a domestic species.

Detecting the footprint of selection on the genomes of Murciano-Granadina goats.

Artificial selection is one of the major forces modifying the genetic composition of livestock populations. Identifying genes under selection could be useful to elucidate their impact on phenotypic variation. We aimed to identify genomic regions targeted by selection for dairy and pigmentation traits in Murciano-Granadina goats. Performance of a selection scan based on the integrated haplotype score test in a population of 1183 Murciano-Granadina goats resulted in the identification of 77 candidate genomic regions/SNPs. The most significant selective sweeps mapped to chromosomes 1 (69.86 Mb), 4 (41.80-49.95 Mb), 11 (65.74 Mb), 12 (31.24 and 52.51 Mb), 17 (34.76-37.67 Mb), 22 (31.75 Mb), and 26 (26.69-31.05 Mb). By using previously generated RNA-Seq data, we built a catalogue of 6414 genes that are differentially expressed across goat lactation (i.e. 78 days post-partum, early lactation; 216 days post-partum, late lactation; 285 days post-partum, dry period). Interestingly, 183 of these genes mapped to selective sweeps and several of them display functions related with lipid, protein, and carbohydrate metabolism, insulin signaling, cell proliferation, as well as mammary development and involution. Of particular interest are the CSN3 and CSN1S2 genes, which encode two major milk proteins. Additionally, we found three pigmentation genes (GLI3, MC1R, and MITF) co-localizing with selective sweeps. Performance of a genome-wide association study and Sanger sequencing and TaqMan genotyping experiments revealed that the c.801C>G (p.Cys267Trp) polymorphism in the melanocortin 1 receptor (MC1R) gene is the main determinant of the black (GG or GC genotypes) and brown (CC genotypes) colorations of Murciano-Granadina goats.

Comparing the diversity of the casein genes in the Asian mouflon and domestic sheep.

We aimed to determine whether casein variants that are currently segregating in ovine populations existed before the domestication of sheep or, to the contrary, if their emergence is much more recent. To this end, we have retrieved whole-genome sequences from Iranian and domestic sheep from Africa, Europe, South and East Asia and West Asia. Population structure analysis based on 55,352,935 SNPs revealed a clear separation between Iranian mouflons and domestic sheep. Moreover, we also observed a strong genetic differentiation between Iranian mouflons sampled in geographic areas close to Tehran and Tabriz. Based on sequence data, hundreds of SNPs mapping to the casein αS1 (CSN1S1, 248 SNPs), casein αS2 (CSN1S2, 268 SNPs), casein ß (CSN2, 146 SNPs) and casein κ (CSN3, 112 SNPs) genes were identified. Approximately 25-63.02% of the casein variation was shared between Iranian mouflons and domestic sheep, and the four domestic sheep populations also shared 44.2-57.4% of the casein polymorphic sites. These findings suggest that an important fraction of the casein variation present in domestic sheep was already segregating in the mouflon prior to its domestication. Genomic studies performed in horses and dogs are consistent with this view, suggesting that much of the diversity that we currently detect in domestic animals comes from standing variation already segregating in their wild ancestors.

An association analysis for 14 candidate genes mapping to meat quality quantitative trait loci in a Duroc pig population reveals that the ATP1A2 genotype is highly associated with muscle electric conductivity.

In previous GWAS carried out in a Duroc commercial line (Lipgen population), we detected on pig chromosomes 3, 4 and 14 several QTL for gluteus medius muscle redness (GM a*), electric conductivity in the longissimus dorsi muscle (LD CE) and vaccenic acid content in the LD muscle (LD C18:1 n - 7), respectively. We have genotyped, in the Lipgen population, 19 SNPs mapping to 14 genes located within these QTL. Subsequently, association analyses have been performed. After correction for multiple testing, two SNPs in the TGFBRAP1 (rs321173745) and SELENOI (rs330820437) genes were associated with GM a*, whereas ACADSB (rs81449951) and GPR26 (rs343087568) genotypes displayed significant associations with LD vaccenic content. Moreover, the polymorphisms located at the ATP1A2 (rs344748241), ATP8B2 (rs81382410) and CREB3L4 (rs321278469 and rs330133789) genes showed significant associations with LD CE. We made a second round of association analyses including the SNPs mentioned above as well as other SNPs located in the chromosomes to which they map. After performing a correction for multiple testing, the only association that remained significant at the chromosome-wide level was that between the ATP1A2 genotype and LD CE. From a functional point of view, this association is meaningful because this locus encodes a subunit of the Na+ /K+ -ATPase responsible for maintaining an electrochemical gradient across the plasma membrane.

Detection of homozygous genotypes for a putatively lethal recessive mutation in the porcine argininosuccinate synthase 1 (ASS1) gene.

The sequencing of the pig genome revealed the existence of homozygous individuals for a nonsense mutation in the argininosuccinate synthase 1 (ASS1) gene (rs81212146, c.944T>A, L315X). Paradoxically, an AA homozygous genotype for this polymorphism is expected to abolish the function of the ASS1 enzyme that participates in the urea cycle, leading to citrullinemia, hyperammonemia, coma and death. Sequencing of five Duroc boars that sired a population of 350 Duroc barrows revealed the segregation of the c.944T>A polymorphism, so we aimed to investigate its phenotypic consequences. Genotyping of this mutation in the 350 Duroc barrows revealed the existence of seven individuals homozygous (AA) for the nonsense mutation. These AA pigs had a normal weight despite the fact that mild citrullinemia often involves impaired growth. Sequencing of the region surrounding the mutation in TT, TA and AA individuals revealed that the A substitution in the second position of the codon (c.944T>A) is in complete linkage disequilibrium with a C replacement (c.943T>C) in the first position of the codon. This second mutation would compensate for the potentially damaging effect of the c.944T>A replacement. In fact, this is the most probable reason why pigs with homozygous AA genotypes at the 944 site of the ASS1 coding region are alive. Our results illustrate the complexities of predicting the consequences of nonsense mutations on gene function and phenotypes, not only because of annotation issues but also owing to the existence of genetic mechanisms that sometimes limit the penetrance of highly harmful mutations.

Low genome-wide homozygosity in 11 Spanish ovine breeds.

The population of Spanish sheep has decreased from 24 to 15 million heads in the last 75 years due to multiple social and economic factors. Such a demographic reduction might have caused an increase in homozygosity and inbreeding, thus limiting the viability of local breeds with excellent adaptations to harsh ecosystems. The main goal of our study was to investigate the homozygosity patterns of 11 Spanish ovine breeds and to elucidate the relationship of these Spanish breeds with reference populations from Europe, Africa and the Near East. By using Ovine SNP50 BeadChip data retrieved from previous publications, we have found that the majority of studied Spanish ovine breeds have close genetic relatedness with other European populations; the one exception is the Canaria de Pelo breed, which is similar to North African breeds. Our analysis has also demonstrated that, with few exceptions, the genomes of Spanish sheep harbor fewer than 50 runs of homozygosity (ROH) with a total length of less than 350 Mb. Moreover, the frequencies of very long ROH (>30 Mb) are very low, and the inbreeding coefficients (FROH ) are generally small (FROH  < 0.10), ranging from 0.008 (Rasa Aragonesa) to 0.086 (Canaria de Pelo). The low levels of homozygosity observed in the 11 Spanish sheep under analysis might be due to their extensive management and the high number of small to medium farms.

Genomic analysis of the origins of extant casein variation in goats.

The variation in the casein genes has a major impact on the milk composition of goats. Even though many casein polymorphisms have been identified so far, we do not know yet whether they are evolutionarily ancient (i.e., they existed before domestication) or young (i.e., they emerged after domestication). Herewith, we identified casein polymorphisms in a data set of 106 caprine whole-genome sequences corresponding to bezoars (Capra aegagrus, the ancestor of domestic goats) and 4 domestic goat (Capra hircus) populations from Europe, Africa, the Far East, and the Near East. Domestic and wild goat populations shared a substantial number of casein SNP, from 36.1% (CSN2) to 55.1% (CSN1S2). The comparison of casein variation among bezoars and the 4 domestic goat populations demonstrated that more than 50% of the casein SNP are shared by 2 or more populations, and 18 to 44% are shared by all populations. Moreover, the majority of casein alleles reported in domestic goats also segregate in the bezoar, including several alleles displaying significant associations with milk composition (e.g., the A/B alleles of the CSN1S1 and CSN3 genes, the A allele of the CSN2 gene). We conclude that much of the current diversity of the caprine casein genes comes from ancient standing variation segregating in the ancestor of modern domestic goats.

Goat domestication and breeding: a jigsaw of historical, biological and molecular data with missing pieces.

Domestic goats (Capra hircus) are spread across the five continents with a census of 1 billion individuals. The worldwide population of goats descends from a limited number of bezoars (Capra aegagrus) domesticated 10 000 YBP (years before the present) in the Fertile Crescent. The extraordinary adaptability and hardiness of goats favoured their rapid spread over the Old World, reaching the Iberian Peninsula and Southern Africa 7000 YBP and 2000 YBP respectively. Molecular studies have revealed one major mitochondrial haplogroup A and five less frequent haplogroups B, C, D, F and G. Moreover, the analysis of autosomal and Y-chromosome markers has evidenced an appreciable geographic differentiation. The implementation of new molecular technologies, such as whole-genome sequencing and genome-wide genotyping, allows for the exploration of caprine diversity at an unprecedented scale, thus providing new insights into the evolutionary history of goats. In spite of a number of pitfalls, the characterization of the functional elements of the goat genome is expected to play a key role in understanding the genetic determination of economically relevant traits. Genomic selection and genome editing also hold great potential, particularly for improving traits that cannot be modified easily by traditional selection.

Investigating the genetic regulation of the expression of 63 lipid metabolism genes in the pig skeletal muscle.

A comprehensive and systematic view of the genetic regulation of lipid metabolism genes is still lacking in pigs. Herewith, we have investigated the genetic regulation of 63 porcine genes with crucial roles in the uptake, transport, synthesis and catabolism of lipids. With this aim, we have performed an expression QTL (eQTL) scan in 104 pigs with available genotypes for the Illumina Porcine SNP60 chip and microarray measurements of gene expression in the gluteus medius muscle. Analysis of the data with gemma software revealed 13 cis- and 18 trans-eQTL modulating the expression of 19 loci. Genes regulated by eQTL participated in a wide array of lipid metabolism pathways such as the ß-oxidation of fatty acids, lipid biosynthesis and lipolysis, fatty acid activation and desaturation, lipoprotein uptake, apolipoprotein assembly and cholesterol trafficking. These data provide a first picture of the genetic regulation of loci involved in porcine lipid metabolism.

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 (hSNP2= 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 3rd and 4th 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 3rd and 4th 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.

Variations at regulatory regions of the milk protein genes are associated with milk traits and coagulation properties in the Sarda sheep.

Regulatory variation at the ovine casein genes could have important effects on the composition and coagulation properties of milk. Herewith, we have partially resequenced the promoters and the 3'-UTR of the four casein genes in 25 Sarda sheep. Alignment of these sequences allowed us to identify a total of 29 SNPs. This level of polymorphism (one SNP every 250 bp) is remarkably high if compared with SNP densities estimated in human genic regions (approximately one SNP per bp). The 29 SNPs identified in our resequencing experiment, plus three previously reported SNPs mapping to the lactalbumin, alpha (LALBA) and ß-lactoglobulin (BLG, also known as progestagen-associated endometrial protein, PAEP) genes, were genotyped with a multiplex TaqMan Open Array Real-Time PCR assay in 760 Sarda sheep with records for milk composition and coagulation properties. Association analysis revealed the existence of significant associations of CSN1S2 and CSN3 genotypes with milk protein and casein contents. Moreover, genotypes at CSN1S1 were significantly associated with rennet coagulation time, curd firming time and curd firmness, whereas CSN2 was associated with curd firming time. These results suggest that SNPs mapping to the promoters and 3'-UTRs of ovine casein genes may exert regulatory effects on gene expression and that they could be used for improving sheep milk quality and technological traits at the population level through marker assisted selection.

East African pigs have a complex Indian, Far Eastern and Western ancestry.

In this study, we have characterized the mitochondrial diversity of 81 swine from Uganda. Median-joining network analysis of D-loop sequences from these individuals and others characterized in previous studies allowed us to determine that Ugandan pigs cluster with populations from the West (Europe/North Africa), Far East and India. In addition, partial sequencing of the Y-chromosome UTY locus in 18 Ugandan domestic pigs revealed the segregation of a single HY1 lineage that has a cosmopolitan distribution. A Western and Far Eastern ancestry for East African pigs had been already reported, but this is the first study demonstrating an additional contribution from the Indian porcine gene pool. This result is consistent with the high frequency of zebuine alleles in cattle from East Africa. The geographic coordinates of East Africa, at the crossroads of many trading routes that, through the ages, linked Europe, Africa and Asia, might explain the rich and complex genetic heritage of livestock native to this area.

A mitochondrial analysis reveals distinct founder effect signatures in Canarian and Balearic goats.

In the course of human migrations, domestic animals often have been translocated to islands with the aim of assuring food availability. These founder events are expected to leave a genetic footprint that may be recognised nowadays. Herewith, we have examined the mitochondrial diversity of goat populations living in the Canarian and Balearic archipelagos. Median-joining network analysis produced very distinct network topologies for these two populations. Indeed, a majority of Canarian goats shared a single ancestral haplotype that segregated in all sampled islands, suggesting a single founder effect followed by a stepping-stone pattern of diffusion. This haplotype also was present in samples collected from archaeological assemblies at Gran Canaria and Lanzarote, making evident its widespread distribution in ancient times. In stark contrast, goats from Majorca and Ibiza did not share any mitochondrial haplotypes, indicating the occurrence of two independent founder events. Furthermore, in Majorcan goats, we detected the segregation of the mitochondrial G haplogroup that has only been identified in goats from Egypt, Iran and Turkey. This finding suggests the translocation of Asian and/or African goats to Majorca, possibly as a consequence of the Phoenician and Carthaginian colonisations of this island.

Mining the pig genome to investigate the domestication process.

Pig domestication began around 9000 YBP in the Fertile Crescent and Far East, involving marked morphological and genetic changes that occurred in a relatively short window of time. Identifying the alleles that drove the behavioural and physiological transformation of wild boars into pigs through artificial selection constitutes a formidable challenge that can only be faced from an interdisciplinary perspective. Indeed, although basic facts regarding the demography of pig domestication and dispersal have been uncovered, the biological substrate of these processes remains enigmatic. Considerable hope has been placed on new approaches, based on next-generation sequencing, which allow whole-genome variation to be analyzed at the population level. In this review, we provide an outline of the current knowledge on pig domestication by considering both archaeological and genetic data. Moreover, we discuss several potential scenarios of genome evolution under the complex mixture of demography and selection forces at play during domestication. Finally, we highlight several technical and methodological approaches that may represent significant advances in resolving the conundrum of livestock domestication.

Application of the microarray technology to the transcriptional analysis of muscle phenotypes in pigs.

The transcriptome refers to the collection of all transcripts present in a cell. Gene expression has a very dynamic nature; it acts as a bridge between epigenetic marks, DNA sequence and proteins and changes to accommodate the requirements of the cell at each given time. Recent technological advances have created new opportunities to study complex phenotypes from a global point of view. From an animal production perspective, muscle transcriptomics has been investigated in relation to muscle growth, carcass fattening and meat quality traits. In this review, we discuss the impact of nutritional, anatomical and genetic factors on muscle gene expression and meat quality of pigs assessed by microarray technologies. Altogether, several common themes have been revealed by the in-depth analysis of the current body of knowledge, for instance, the involvement of genes related to energy balance and substrate turnover in the oxidative/glycolytic phenotype of red/white muscle fibre types and in the storage of intramuscular fat. The review also covers recent advances in the discovery of expression QTL and regulatory RNAs in porcine breeds as well as technical developments in the field of deep-sequencing technologies that are expected to substantially increase our knowledge about the genetic architecture of meat quality and production traits.

Technical note: Advantages and limitations of authenticating Palmera goat dairy products by pyrosequencing the melanocortin 1 receptor (MC1R) gene.

Inferring the breed of origin of dairy products can be achieved through molecular analysis of genetic markers with a population-specific pattern of segregation. The goal of the current work was to generate such markers in goats by resequencing several pigmentation genes [melanocortin 1 receptor (MC1R), v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog (KIT), tyrosinase (TYR), and tyrosinase-related protein 2 (TYRP2)]. This experiment revealed 10 single nucleotide polymorphisms (SNP), including 5 missense mutations and 1 nonsense mutation. These markers were genotyped in 560 goats from 18 breeds originally from Italy, the Iberian Peninsula, the Canary Islands, and North Africa. Although the majority of SNP segregated at moderate frequencies in all populations (including 2 additional markers that were used as a source of information), we identified a c.764G>A SNP in MC1R that displayed highly divergent allelic frequencies in the Palmera breed compared with the Majorera and Tinerfeña breeds from the Canary Islands. Thus, we optimized a pyrosequencing-based technique that allowed us to estimate, very accurately, the allele frequencies of this marker in complex DNA mixtures from different individuals. Once validated, we applied this method to generating breed-specific DNA profiles that made it possible to detect fraudulent cheeses in which Palmero cheese was manufactured with milk from Majorera goats. One limitation of this approach, however, is that it cannot be used to detect illegal manufacturing where Palmero dairy products are produced by mixing milk from Palmera and Majorera goats, because the c.764G>A SNP segregates in both breeds.

DECR1 and ME1 genotypes are associated with lipid composition traits in Duroc pigs.

Variation at the porcine DECR1 and ME1 genes has been associated with meat quality traits and backfat thickness in Landrace pigs, respectively. However, it has not been investigated yet whether DECR1 and ME1 genotypes influence lipid composition. With this aim, we have genotyped two missense DECR1 substitutions (c.160G>C and c.437G>C) and one silent ME1 (c.576C>T) polymorphism in 361 Duroc barrows distributed in five half-sib families and phenotyped for serum lipid concentrations and intramuscular fat content and composition traits. At the whole-population level, relevant associations, that is, with a posterior probability of the allele substitution effect to be over or below zero (PPN0) > 0.90, were observed between DECR1 genotype and serum cholesterol (CHOL) (PPN0 = 0.932) and LDL concentrations (PPN0 = 0.945) at 190 days, as well as between ME1 genotype and longissimus dorsi saturated fatty acid content (PPN0 = 0.924). At the within-family level, we found relevant associations between DECR1 and ME1 genotypes and diverse lipid composition traits, but most of them were family-specific. Discrepancies in allele substitution effects estimated in half-sib families might be produced by many factors such as number of individuals, marker allele frequencies and informativeness in each family, unaccounted random genetic and environmental effects, epistasis and family-specific differences in the linkage phase or amount of linkage disequilibrium between causal and marker mutations. This lack of consistency across families, combined with the fact that the ME1 mutation is synonymous and that the two DECR1 polymorphisms are conservative, suggests that the associations found are not causative.

Genetics of serum and muscle lipids in pigs.

Pork meat is one of the most important sources of animal protein in the human diet. Its nutritional properties are partly determined by intramuscular fat content and composition, with existing general consensus about the detrimental effects of cholesterol and saturated fat on cardiovascular health in humans. Because of their physiological resemblance, pigs can be also used as a valuable animal model to study the genetics of human diseases such as atherosclerosis, obesity and dyslipidaemias. Heritability estimates and QTL maps of porcine muscle and serum lipid traits evidence that a considerable amount of genetic variance determining these phenotypes exists, but its molecular basis remains mostly unknown. The recent advent of high-throughput genotyping and sequencing technologies has revolutionised the field of animal genomics. With these powerful tools, finding needles in the genomic haystack has become increasingly feasible. However, these methodological advances should not be deemed as magic bullets. The goal of identifying the many polymorphisms that shape the variability of lipid phenotypes is so challenging that success can be achieved only under the scope of large international consortia.

Copy number variation in the genomes of domestic animals.

Copy number variation (CNV) might be one of the main contributors to phenotypic diversity and evolutionary adaptation in animals and plants, employing a wide variety of mechanisms, such as gene dosage and transcript structure alterations, to modulate organismal plasticity. In the past 4 years, considerable advances have been made in the characterization of the genomic architecture of CNV in domestic species. First, low-resolution CNV maps were produced for cattle, goat, sheep, pig, dog, chicken, duck and turkey, showing that these structural polymorphisms comprise a significant part of these genomes. Furthermore, CNVs have been associated with several pigmentation (white coat in horse, pig and sheep) and morphological (late feathering and pea comb in chicken) traits, as well as with susceptibility to a wide array of diseases and developmental disorders, for example osteopetrosis, anhidrotic ectodermal dysplasia, copper toxicosis, intersexuality, cone degeneration, periodic fever and dermoid sinus, among others. In the future, development of high-resolution tools for CNV detection and typing combined with the implementation of databases integrating CNV, QTL and gene expression data will be essential to identify and measure the impact of this source of structural variation on the many phenotypes that are relevant to animal breeders and veterinary practitioners.

Quantitative trait loci analysis of a Duroc commercial population highlights differences in the genetic determination of meat quality traits at two different muscles.

We performed a whole-genome scan with 110 informative microsatellites in a commercial Duroc population for which growth, fatness, carcass and meat quality phenotypes were available. Importantly, meat quality traits were recorded in two different muscles, that is, gluteus medius (GM) and longissimus thoracis et lumborum (LTL), to find out whether these traits are determined by muscle-specific genetic factors. At the whole-population level, three genome-wide QTL were identified for carcass weight (SSC7, 60 cM), meat redness (SSC13, 84 cM) and yellowness (SSC15, 108 cM). Within-family analyses allowed us to detect genome-wide significant QTL for muscle loin depth between the 3rd and 4th ribs (SSC15, 54 cM), backfat thickness (BFT) in vivo (SSC10, 58 cM), ham weight (SSC9, 69 cM), carcass weight (SSC7, 60 cM; SSC9, 68 cM), BFT on the last rib (SSC11, 48 cM) and GM redness (SSC8, 85 cM; SSC13, 84 cM). Interestingly, there was low positional concordance between meat quality QTL maps obtained for GM and LTL. As a matter of fact, the three genome-wide significant QTL for colour traits (SSC8, SSC13 and SSC15) that we detected in our study were all GM specific. This result suggests that QTL effects might be modulated to a certain extent by genetic and environmental factors linked to muscle function and anatomical location.