An association study using imputed whole‐genome sequence data identifies novel significant loci for growth‐related traits in a Duroc × Erhualian F2 population

The average daily gain (ADG) and body weight (BW) are very important traits for breeding programs and for the meat production industry, which have attracted many researchers to delineate the genetic architecture behind these traits. In the present study, single‐ and multi‐trait genome‐wide association studies (GWAS) were performed between imputed whole‐genome sequence data and the traits of the ADG and BW at different stages in a large‐scale White Duroc × Erhualian F₂ population. A bioinformatics annotation analysis was used to assist in the identification of candidate genes that are associated with these traits. Five and seven genome‐wide significant quantitative trait loci (QTLs) were identified by single‐ and multi‐trait GWAS, respectively. Furthermore, more than 40 genome‐wide suggestive loci were detected. On the basis of the whole‐genome sequence association study and the bioinformatics analysis, NDUFAF6, TNS1 and HMGA1 stood out as the strongest candidate genes. The presented single‐ and multi‐trait GWAS analysis using imputed whole‐genome sequence data identified several novel QTLs for pig growth‐related traits. Integrating the GWAS with bioinformatics analysis can facilitate the more accurate identification of candidate genes. Higher imputation accuracy, time‐saving algorithms, improved models and comprehensive databases will accelerate the identification of causal genes or mutations, which will contribute to genomic selection and pig breeding in the future.     

Genome‐wide association studies reveal additional related loci for fatty acid composition in a Duroc pig multigenerational population

The aim of the present study was to detect quantitative trait loci affecting fatty acid composition in back fat and intramuscular fat in a Duroc pig population comprising seventh‐generation pedigrees using genome‐wide association studies (GWAS). In total, 305 animals were genotyped using single nucleotide polymorphisms (SNPs) array and five selected SNPs from regions containing known candidate genes related to fatty acid synthesis or metabolism. In total, 24 genome‐wide significant SNP regions were detected in 12 traits, and 76 genome‐wide suggestive SNP regions were detected in 33 traits. The Sus scrofa chromosome (SSC) 7 at 10.3 Mb was significantly associated with C17:0 in intramuscular fat, while the SSC9 at 13.6 Mb was significantly associated with C14:0 in intramuscular fat. The SSC12 at 1.0 Mb was significantly associated with C14:0 in back fat and the SSC14 at 121.0 Mb was significantly associated with C18:0 in intramuscular fat. These regions not only replicated previously reported loci containing some candidate genes involved in fatty acid composition (fatty acid synthase and stearoyl‐CoA desaturase) but also included several additional related loci.     

A genome‐wide SNP scan in a porcine Large White × Minzhu intercross population reveals a locus influencing muscle mass on chromosome 2

A high‐density single nucleotide polymorphism (SNP) array containing 62 163 markers was employed for a genome‐wide association study (GWAS) to identify variants associated with lean meat in ham (LMH, %) and lean meat percentage (LMP, %) within a porcine Large White × Minzhu intercross population. For each individual, LMH and LMP were measured after slaughter at the age of 240 ± 7 days. A total of 557 F2 animals were genotyped. The GWAS revealed that 21 SNPs showed significant genome‐wide or chromosome‐wide associations with LMH and LMP by the Genome‐wide Rapid Association using Mixed Model and Regression‐Genomic Control approach. Nineteen significant genome‐wide SNPs were mapped to the distal end of Sus Scrofa Chromosome (SSC) 2, where a major known gene responsible for muscle mass, IGF2 is located. A conditioned analysis, in which the genotype of the strongest associated SNP is included as a fixed effect in the model, showed that those significant SNPs on SSC2 were derived from a single quantitative trait locus. The two chromosome‐wide association SNPs on SSC1 disappeared after conditioned analysis suggested the association signal is a false association derived from using a F2 population. The present result is expected to lead to novel insights into muscle mass in different pig breeds and lays a preliminary foundation for follow‐up studies for identification of causal mutations for subsequent application in marker‐assisted selection programs for improving muscle mass in pigs.     

Genome‐wide association studies for two exterior traits in Chinese Dongxiang spotted pigs

Pigs show extensive variation in exterior appearance. This variation has been explored as one of the selection targets to form breeding features in the pig industry. In this study, a customized Affymetrix Axiom Pig1.4M array plate was used to conduct genome‐wide association studies (GWAS) for two exterior traits—spotted coat color (SC) and facial type (FT)—in Chinese Dongxiang spotted pigs. Two and 12 single nucleotide polymorphisms (SNPs) were identified to be associated with SC and FT at the 5% genome‐wide significant level, respectively. Of these SNPs, two associated with SC were located around the DACH1 gene on Sus Scrofa chromosome (SSC)11. Eleven SNPs associated with FT were located within a 7.68 Mb region (29.03–36.71 Mb) on SSC7. Based on the GWAS results and the biological functions of genes, we highlight EDNRB as a candidate gene for SC and HMGA1 and RPS10 as the potential genes affecting facial variation. The findings contribute to the final characterization of causative genes and mutations underlying the effects of the SSC7 and SSC11 loci and improve our understanding of the genetic basis of phenotypic variation in Chinese indigenous pigs.     

A genome scan for quantitative trait loci affecting the length of small intestine in a White Duroc × Chinese Erhualian intercross resource population

The small intestine is a vital organ in animal gastrointestinal system, in which a large variety of nutrients are absorbed. To identify quantitative trait loci (QTL) for the length of porcine small intestine, phenotypic values were measured in 1034 individuals at 240 d from a White Duroc × Chinese Erhualian intercross F₂ population. The length of small intestine showed strong correlation with growth traits and carcass length in the F2 population. A whole‐genome scan was performed based on 183 microsatellites covering the pig genome in the F₂ population. A total of 10 QTL for this trait were identified on 8 pig chromosomes (SSC), including four 1% genome‐wide significant QTL on SSC2, 4, 7 and 8, one 5% genome‐wide significant QTL on SSC12, and five 5% chromosome‐wide significant QTL on SSC5, 7, 13 and 14. The Erhualian alleles were generally associated with shorter length of the small intestine except the alleles on SSC7 and 13. The QTL on SSC4 overlapped with the previously reported QTL for the length of small intestine. Several significant QTL on SSC2, 8, and 12 were consistent with previous reports. The significant QTL detected on SSC7 was reported for the first time. All QTL identified in this study corresponded to the known region significantly associated with growth traits, supporting the important role of the length of small intestine in pig growth.     

Unravelling the genetic loci for growth and carcass traits in Chinese Bamaxiang pigs based on a 1.4 million SNP array

Bamaxiang pig is from Guangxi province in China, characterized by its small body size and two‐end black coat colour. It is an important indigenous breed for local pork market and excellent animal model for biomedical research. In this study, we performed genomewide association studies (GWAS) on 43 growth and carcass traits in 315 purebred Bamaxiang pigs based on a 1.4 million SNP array. We observed considerable phenotypic variability in the growth and carcass traits in the Bamaxiang pigs. The corresponding SNP based heritability varied greatly across the 43 traits and ranged from 9.0% to 88%. Through a conditional GWAS, we identified 53 significant associations for 35 traits at p value threshold of 10^?6. Among which, 26 associations on chromosome 3, 7, 14 and X passed a genomewide significance threshold of 5 × 10^?8. The most remarkable loci were at around 30.6 Mb on chromosome 7, which had growth stage‐dependent effects on body lengths and cannon circumferences and showed large effects on multiple carcass traits. We discussed HMGA1 NUDT3, EIF2AK1, TMEM132C and AFF2 that near the lead SNP of significant loci as plausible candidate genes for corresponding traits. We also showed that including phenotypic covariate in GWAS can help to reveal additional significant loci for the target traits. The results provide insight into the genetic architecture of growth and carcass traits in Bamaxiang pigs.     

利用GWAS筛选影响野猪被毛表型变异的候选功能基因

识别野生动物群体内潜在影响动物表型变异的相关基因是进化遗传学研究的主要目的,而动物毛色是研究动物被毛表型形成遗传机制的最佳模型之一。应用Illumina公司提供的猪60 k SNP基因芯片对选取的62只不同被毛表型的野猪个体进行基因分型,利用SNP分析结果,通过对照全基因组关联分析（GWAS）识别影响野猪被毛表型差异的相关变异。结果表明,识别了6个与野猪被毛表型相关的基因组变异区域,分别位于SSC1（ALGA0001794,ASGA0006416）、SSC2（ASGA0011559）、SSC6（H3GA0018683）、SSC7（ASGA0035535）和SSC14（ASGA0060641）;最显著相关的SNP（ALGA0001794）位于猪1号染色体上（SSC1）的27 899 596-27 899 696 bp区间（P=2.96×10-（-5））。该研究初步鉴定了6个与野猪毛色性状相关的易感位点,为进一步研究野猪不同毛色性状的形成机制提供了基础。     

Genomic regions associated with ventro-cranial chronic pleuritis in pig

Ventro‐cranial chronic pleuritis can be a result of pleuropneumonia and enzootic pneumonia. These diseases cause severe losses in intensive pig production worldwide, but host resistance is difficult to breed for. It could be beneficial to use marker‐assisted selection, and a step towards this is to identify genomic regions associated with the trait. For this purpose, 7304 pigs from 11 boar families were analysed for associations between single nucleotide polymorphisms and ventro‐cranial chronic pleuritis. The pigs were genotyped by the use of the iSelect Custom 7 K porcine SNP Chip. Quantitative trait loci (QTL), significant at the chromosome‐wide level, were identified on Sus scrofa chromosomes (SSC) 2, 4, 11, 12 and 13 in four different boar families. The QTL on SSC 4 in family G was also significant at the genome‐wide threshold according to Bonferroni correction. We have identified a number of candidate genes, but the causative mutations still need to be identified. Markers closely associated with the resistance traits have a strong potential for use in breeding towards animals with improved characteristics concerning ventro‐cranial chronic pleuritis.     

Mapping QTL for porcine muscle fibre traits in a White Duroc × Erhualian F2 resource population

Muscle fibre traits are related with meat quality in meat animals. In this study, a whole‐genome scan with 183 microsatellite markers covering the pig genome was performed to identify quantitative trait loci (QTL) for cross‐sectional area, numerical percentage and relative area of type I, IIA and IIB myofibres, fibre number per square centimetre and total fibre number in the longissimus muscle by using 120 F₂ animals in a White Duroc × Erhualian intercross. In total, 20 QTL were mapped on pig chromosomes (SSC) 1, 2, 7, 8, 9, 11, 15, 16 and X, of which eight reached genome‐wide significance levels and explained large proportions (6.53–34.63%) of phenotypic variance. Five QTL detected in this study confirmed the previous QTL reports and the others were detected for the first time. Chinese Erhualian alleles are generally associated with muscle fibre traits favourable for meat quality.     

Combined approaches to identify genomic regions involved in phenotypic differentiation between low divergent breeds: Application in Sardinian sheep populations

Selective breeding has led to modifications in the genome of many livestock breeds. In this study, we identified the genomic regions that may explain some of the phenotypic differences between two closely related breeds from Sardinia. A total of 44 animals, 20 Sardinian Ancestral Black (SAB) and 24 Sardinian White (SW), were genotyped using the Illumina Ovine 50K array. A total of 68, 38 and 15 significant markers were identified using the case–control genome‐wide association study (GWAS), the Bayesian population differentiation analysis (F_ST) and the Rsb metric, respectively. Comparisons among the approaches revealed a total of 22 overlapping markers between GWAS and F_ST and one marker between GWAS and Rsb. Three markers detected by Rsb were also located near (<2 Mb) to highly significant regions identified by GWAS and F_ST analyses. Moreover, one candidate marker identified by GWAS and F_ST approaches was located in a run of homozygosity island that was shared by both breeds. We identified several genes involved in many phenotypic differences (such as stature and growth, reproduction, ear size, coat colour, behaviour) between the two analysed breeds. This study shows that combining several genome‐wide approaches could improve discovery of regions involved in the variability of breeding traits and responsible for the phenotypic diversity even between closely related breeds. Overall, the combination of such genome‐wide methods can be extended to other livestock breeds that share between them a similar genetic background, to understand the process that shapes the patterns of genetic variability between closely related populations.     

Genome‐wide association study for egg production and quality in layer chickens

Discovery of genes with large effects on economically important traits has for many years been of interest to breeders. The development of SNP panels which cover the whole genome with high density and, more importantly, that can be genotyped on large numbers of individuals at relatively low cost, has opened new opportunities for genome‐wide association studies (GWAS). The objective of this study was to find genomic regions associated with egg production and quality traits in layers using analysis methods developed for the purpose of whole genome prediction. Genotypes on over 4500 birds and phenotypes on over 13 000 hens from eight generations of a brown egg layer line were used. Birds were genotyped with a custom 42K Illumina SNP chip. Recorded traits included two egg production and 11 egg quality traits (puncture score, albumen height, yolk weight and shell colour) at early and late stages of production, as well as body weight and age at first egg. Egg weight was previously analysed by Wolc et al. ( 2012 ). The Bayesian whole genome prediction model – BayesB (Meuwissen et al. 2001 ) was used to locate 1 Mb regions that were most strongly associated with each trait. The posterior probability of a 1 Mb window contributing to genetic variation was used as the criterion for suggesting the presence of a quantitative trait locus (QTL) in that window. Depending upon the trait, from 1 to 7 significant (posterior probability >0.9) 1 Mb regions were found. The largest QTL, a region explaining 32% of genetic variance, was found on chr4 at 78 Mb for body weight but had pleiotropic effects on other traits. For the other traits, the largest effects were much smaller, explaining <7% of genetic variance, with regions on chromosomes 2, 12 and 17 explaining above 5% of genetic variance for albumen height, shell colour and egg production, respectively. In total, 45 of 1043 1 Mb windows were estimated to have a non‐zero effect with posterior probability > 0.9 for one or more traits.     

Genome-wide linkage and QTL mapping in porcine F2 families generated from Pietrain, Meishan and Wild Boar crosses

Three informative pig F₂ families based on European Wild Boar (W), Meishan (M) and Pietrain (P) crosses have been used for genome‐wide linkage and quantitative trait loci (QTL) analysis. Altogether 129 microsatellites, 56 type I loci and 46 trait definitions (specific to growth, fattening, fat deposition, muscling, meat quality, stress resistance and body conformation) were included in the study. In the linkage maps of M × P, W × P and W × M families, average spacing of markers were 18.4, 19.7 and 18.8 cM, the numbers of informative meioses were 582, 534 and 625, and the total lengths of autosomes measured were 27.3, 26.0 and 26.2 Morgan units, respectively. Maternal maps were on average 1.3 times longer than paternal maps. QTLs contributing more than 3% of F₂ phenotypic variance could be identified at p < 0.05 chromosome‐wide level. Differences in the numbers and positions of QTLs were observed between families. Genome‐wide significant QTL effects were mapped for growth and fattening traits on eight chromosomes (1, 2, 4, 13, 14, 17, 18 and X), for fat deposition traits on seven chromosomes (1, 2, 3, 4, 6, 7 and X), for muscling traits on 11 chromosomes (1, 2, 3, 4, 6, 7, 8, 12, 14, 15 and X), for meat quality and stress resistance traits on seven chromosomes (2, 3, 6, 13, 16, 18 and X), and QTLs for body‐conformation traits were detected on 14 chromosomes. Closely correlated traits showed similar QTL profiles within families. Major QTL effects for meat quality and stress resistance traits were found on SSC6 in the interval RYR1‐A1BG in the W × P and M × P families, and could be attributed to segregation of the RYR1 allele T derived from Pietrain, whereas no effect in the corresponding SSC6 interval was found in family W × M, where Wild Boar and Meishan both contributed the RYR1 allele C. QTL positions were mostly similar in two of the three families for body conformation traits and for growth, fattening, fat deposition and muscling traits, especially on SSC4 (interval SW1073‐NGFB). QTLs with large effects were also mapped on SSC7 in the major histocompatibility complex (MHC) (interval CYP21A2‐S0102) and affected body length, weight of head and many other traits. The identification of DNA variants in genes causative for the QTLs requires further fine mapping of QTL intervals and a positional cloning. However, for these subsequent steps, the genome‐wide QTL mapping in F₂ families represents an essential starting point and is therefore significant for animal breeding.     

Quantitative trait loci analysis for growth and carcass traits in a Meishan x Duroc F2 resource population

We constructed a pig F2 resource population by crossing a Meishan sow and a Duroc boar to locate economically important trait loci. The F2 generation was composed of 865 animals (450 males and 415 females) from four F1 males and 24 F1 females and was genotyped for 180 informative microsatellite markers spanning 2,263.6 cM of the whole pig genome. Results of the genome scan showed evidence for significant quantitative trait loci (<1% genomewise error rate) affecting weight at 30 d and average daily gain on Sus scrofa chromosome (SSC) 6, carcass yield on SSC 7, backfat thickness on SSC 7 and SSC X, vertebra number on SSC 1 and SSC 7, loin muscle area on SSC 1 and SSC 7, moisture on SSC 13, intramuscular fat content on SSC 7, and testicular weight on SSC 3 and SSC X. Moreover, 5% genomewise significant QTL were found for birth weight on SSC 7, average daily gain on SSC 4, carcass length on SSC 6, SSC 7, and SSC X and lightness (L value) on SSC 3. We identified 38 QTL for 28 traits at the 5% genomewise level. Of the 38 QTL, 24 QTL for 17 traits were significant at the 1% genomewise level. Analysis of marker genotypes supported the breed of origin results and provided further evidence that a suggestive QTL for circumference of cannon bone also was segregating within the Meishan parent. We identified genomic regions related with growth and meat quality traits. Fine mapping will be required for their application in introgression programs and gene cloning.     

Evaluations of Boar Gonad Development,Spermatogenesis with regard to Semen Characteristics,Libido and Serum Testosterone Levels based on large White Duroc × Chinese Erhualian Crossbred Boars

Chinese Erhualian pigs are known for prolificacy with distinct reproductive traits compared with Western commercial breeds. In this study, a four‐generation intercross resource population was constructed using White Duroc boars and Chinese Erhualian sows as founder animals, and a total of 14 male reproductive traits were recorded in 411 F₂/F₃ boars including the testis and epididymis weights, the seminiferous tubular diameter and spermatogenesis at 60, 90 and 300 days of age, semen characteristics, serum testosterone concentration and libido level at 300 days of age. The White Duroc–Erhualian boars showed remarkable segregations in the traits measured except for the seminiferous tubular diameter and had high ratio (13.9%) of the abnormality of spermatogenesis, providing a good experimental population for detecting quantitative trait loci affecting these male reproductive traits. Furthermore, the correlations among nine male reproductive traits at 300 days of age indicated that the testis weight and the body weight were strongly correlated with the sperm production, supporting the two traits as important parameters for boar selection to increase sperm production and ultimately improve boar fertility. The libido level in the White Duroc–Erhualian boars that was evaluated by a new and easily recorded scoring system showed a significant correlation with serum testosterone concentration. Yet, both libido and serum testosterone concentration were not correlated with the sperm production. Results of this study provided new information on the male reproductive physiology and genetics in Chinese Erhualian and White Duroc boars.     

荣昌猪初产繁殖性状的全基因组关联研究

旨在鉴定荣昌猪初产繁殖性状的重要变异位点和基因,为荣昌猪繁殖性状的遗传改良提供重要的分子标记和基因资源。本研究选取429头荣昌母猪进行猪50K芯片基因分型,经过质量控制和基因型填充后,保留35 046个SNPs用于分析。采用主成分分析法研究群体结构,利用混合线性模型(mixed-linear model, MLM)将出生年、出生月作为固定效应,将主成分值作为协变量对总产仔数、活产仔数、死胎数和初生窝重性状进行全基因组关联分析(GWAS)。结果显示,在全基因组显著水平上鉴定出2个影响荣昌猪初生窝重的SNPs和1个影响荣昌猪死胎数的SNP;在潜在显著水平上鉴定到5个影响荣昌猪总产仔数的SNPs, 3个影响荣昌猪活产仔数的SNPs和10个影响荣昌猪死胎数的SNPs。通过全基因组关联分析筛选到1个显著的SNP(SSC17:57 315 180 bp)同时影响荣昌猪总产仔数、活产仔数和初生窝重,1个显著的SNP(SSC1:279 214 647 bp)同时影响荣昌猪活产仔数和总产仔数,暗示基因在不同性状间具有一因多效性。本研究根据候选基因的相关分子生物学功能,确定BMP7基因为影响荣昌猪总产仔数...     

Genetic variation at the alpha‐1‐fucosyltransferase (FUT1) gene in Asian wild boar and Chinese and Western commercial pig breeds

Escherichia coli F18 bacteria producing enterotoxins and/or shigatoxin (ETEC/STEC) are main pathogens that cause oedema disease and postweaning diarrhoea in piglets, and alpha‐1‐fucosyltransferase (FUT1) gene has been identified as a candidate gene for controlling the expression of ETEC F18 receptor. The genetic variations at nucleotide position 307 in open reading frame of FUT1 gene in one wild boar breed and 20 western commercial and Chinese native pig breeds were investigated by polymerase chain reaction–restriction fragment length polymorphism. The results showed that the genetic polymorphisms of the FUT1 locus were only detected in western pig breeds and the Chinese Taihu (including Meishan pig, Fengjing pig and Erhualian pig), Huai and Lingao pig breeds; only Duroc and Pietrain possessed the resistant AA genotype, while the wild boar and other Chinese pig breeds only presented the susceptible genotype GG. The results indicated that Chinese native pig breeds lack genetic factors providing resistance to ETEC F18 bacteria. The resistant allele to ETEC F18 might originate from European wild boar. It was inferred that oedema and postweaning diarrhoea caused by ETEC F18 have close relationship with the growth rate, which can explain why on the contrary Chinese native pig breeds have stronger resistance to oedema and postweaning diarrhoea in piglets compared with western pig breeds.     

Pool‐based genome‐wide association study identified novel candidate regions on BTA9 and 14 for oleic acid percentage in Japanese Black cattle

Fatty acid composition is an important indicator of beef quality. The objective of this study was to search the potential candidate region for fatty acid composition. We performed pool‐based genome‐wide association studies (GWAS) for oleic acid percentage (C18:1) in a Japanese Black cattle population from the Hyogo prefecture. GWAS analysis revealed two novel candidate regions on BTA9 and BTA14. The most significant single nucleotide polymorphisms (SNPs) in each region were genotyped in a population (n = 899) to verify their effect on C18:1. Statistical analysis revealed that both SNPs were significantly associated with C18:1 (p = .0080 and .0003), validating the quantitative trait loci (QTLs) detected in GWAS. We subsequently selected VNN1 and LYPLA1 genes as candidate genes from each region on BTA9 and BTA14, respectively. We sequenced full‐length coding sequence (CDS) of these genes in eight individuals and identified a nonsynonymous SNP T66M on VNN1 gene as a putative candidate polymorphism. The polymorphism was also significantly associated with C18:1, but the p value (p = .0162) was higher than the most significant SNP on BTA9, suggesting that it would not be responsible for the QTL. Although further investigation will be needed to determine the responsible gene and polymorphism, our findings would contribute to development of selective markers for fatty acid composition in the Japanese Black cattle of Hyogo.     

Detection and analysis of genome-wide copy number variation in the pig genome using an 80 K SNP Beadchip

Copy number variation (CNV) is an important source of genetic variability in human or animal genomes and play key roles in phenotypic diversity and disease susceptibility. In the present study, we performed a genome-wide analysis for CNV detection using SNP genotyping data of 857 Large White pigs. A total of 312 CNV regions (CNVRs) were detected with the PennCNV algorithm, which covered 57.76 Mb of the pig genome and correspond to 2.36% of the genome sequence. The length of the CNVRs on autosomes ranged from 1.77 Kb to 1.76 Mb with an average of 185.11 Kb. Of these, 220 completely or partially overlapped with 1,092 annotated genes, which enriched a wide variety of biological processes. Comparisons with previously reported pig CNVR revealed 92 (29.49%) novel CNVRs. Experimentally, 80% of CNVRs selected randomly were validated by quantitative PCR (qPCR). We also performed an association analysis between some of the CNVRs and reproductive traits, with results demonstrating the potential importance of CNVR61 and CNVR283 associated with litter sizes. Notably, the GPER1 gene located in CNVR61 plays a key role in reproduction. Our study is an important complement to the CNV map in the pig genome and provides valuable information for investigating the association between genomic variation and economic traits.     

Genome‐wide association study for fatty acid composition in Japanese Black cattle

Fatty acid composition is one of the important traits in beef. The aim of this study was to identify candidate genomic regions for fatty acid composition by genome‐wide association study with 50 K single nucleotide polymorphism (SNP) array in Japanese Black cattle. A total of 461 individuals and 40 657 SNPs were used in this study. We applied genome‐wide rapid association using mixed model and regression (GRAMMAR) and genomic control approaches to estimate the associations between genotypes and fatty acid composition. In addition, two SNPs in fatty acid synthase (FASN) (T1952A) and stearoyl‐CoA desaturase (SCD) (V293A) genes were also genotyped. Association analysis revealed that 30 significant SNPs for several fatty acids (C14:0, C14:1, C16:1 and C18:1) were located in the BTA19 FASN gene located within this region but the FASN mutation had no significant effect on any traits. We also detected one significant SNP for C18:1 on BTA23 and two SNPs for C16:0 on BTA25. The region around 17 Mb on BTA26 harbored two significant SNPs for C14:1 and SNP in SCD in this region showed the strongest association with C14:1. This study demonstrated novel candidate regions in BTA19, 23 and 25 for fatty acid composition.