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.     

Genome‐wide association analysis for growth,muscularity and meat quality in Piétrain pigs

Improvement in growth and meat quality is one of the main objectives in sire line pig breeding programmes. Mapping quantitative trait loci for these traits using experimental crosses and a linkage‐based approach has been performed frequently in the past. The Piétrain breed often was involved as a founder breed to establish the experimental crosses. This breed was selected for muscularity and leanness but shows relatively poor meat quality. It is frequently used as a sire line breed. With the advent of genome‐wide and dense SNP chips in pig genomic research, it is possible to also conduct genome‐wide association studies within the Piétrain breed. In this study, around 500 progeny‐tested sires were genotyped with 60k SNPs. Data filtering showed that around 48k SNPs were useable in this sample. These SNPs were used to conduct a genome‐wide association study for growth, muscularity and meat quality traits. Because it is known that a mutation in the RYR1 gene located on chromosome 6 shows a major effect on meat quality, this mutation was included in the models. Single‐marker and multimarker association analyses were performed. The results revealed between zero and eight significant associations per trait with P < 5 × 10^?5. Of special interest are SNPs located on SSC6, SSC10 and SSC15.     

A two‐step approach to map quantitative trait loci for meat quality in connected porcine F2 crosses considering main and epistatic effects

The aim of this study was to map QTL for meat quality traits in three connected porcine F₂ crosses comprising around 1000 individuals. The three crosses were derived from the founder breeds Chinese Meishan, European Wild Boar and Pietrain. The animals were genotyped genomewide for approximately 250 genetic markers, mostly microsatellites. They were phenotyped for seven meat quality traits (pH at 45 min and 24 h after slaughter, conductivity at 45 min and 24 h after slaughter, meat colour, drip loss and rigour). QTL mapping was conducted using a two‐step procedure. In the first step, the QTL were mapped using a multi‐QTL multi‐allele model that was tailored to analyse multiple connected F₂ crosses. It considered additive, dominance and imprinting effects. The major gene RYR1:g.1843C>T affecting the meat quality on SSC6 was included as a cofactor in the model. The mapped QTL were tested for pairwise epistatic effects in the second step. All possible epistatic effects between additive, dominant and imprinting effects were considered, leading to nine orthogonal forms of epistasis. Numerous QTL were found. The most interesting chromosome was SSC6. Not all genetic variance of meat quality was explained by RYR1:g.1843C>T. A small confidence interval was obtained, which facilitated the identification of candidate genes underlying the QTL. Epistasis was significant for the pairwise QTL on SSC12 and SSC14 for pH24 and for the QTL on SSC2 and SSC5 for rigour. Some evidence for additional pairwise epistatic effects was found, although not significant. Imprinting was involved in epistasis.     

Detection of quantitative trait loci for carcass composition traits in pigs

A quantitative trait locus (QTL) analysis of carcass composition data from a three-generation experimental cross between Meishan (MS) and Large White (LW) pig breeds is presented. A total of 488 F2 males issued from six F1 boars and 23 F1 sows, the progeny of six LW boars and six MS sows, were slaughtered at approximately 80 kg live weight and were submitted to a standardised cutting of the carcass. Fifteen traits, i.e. dressing percentage, loin, ham, shoulder, belly, backfat, leaf fat, feet and head weights, two backfat thickness and one muscle depth measurements, ham + loin and back + leaf fat percentages and estimated carcass lean content were analysed. Animals were typed for a total of 137 markers covering the entire porcine genome. Analyses were performed using a line-cross (LC) regression method where founder lines were assumed to be fixed for different QTL alleles and a half/full sib (HFS) maximum likelihood method where allele substitution effects were estimated within each half-/full-sib family. Additional analyses were performed to search for multiple linked QTL and imprinting effects. Significant gene effects were evidenced for both leanness and fatness traits in the telomeric regions of SSC 1q and SSC 2p, on SSC 4, SSC 7 and SSC X. Additional significant QTL were identified for ham weight on SSC 5, for head weight on SSC 1 and SSC 7, for feet weight on SSC 7 and for dressing percentage on SSC X. LW alleles were associated with a higher lean content and a lower fat content of the carcass, except for the fatness trait on SSC 7. Suggestive evidence of linked QTL on SSC 7 and of imprinting effects on SSC 6, SSC 7, SSC 9 and SSC 17 were also obtained.     

Linkage disequilibrium analysis in the genetically isolated Norfolk Island population

Norfolk Island is a human genetic isolate, possessing unique population characteristics that could be utilized for complex disease gene localization. Our intention was to evaluate the extent and strength of linkage disequilibrium (LD) in the Norfolk isolate by investigating markers within Xq13.3 and the NOS2A gene encoding the inducible nitric oxide synthase. A total of six microsatellite markers spanning approximately 11 Mb were assessed on chromosome Xq13.3 in a group of 56 men from Norfolk Island. Additionally, three single nucleotide polymorphisms (SNPs) localizing to the NOS2A gene were analyzed in a subset of the complex Norfolk pedigree. With the exception of two of the marker pairs, one of which is the most distantly spaced marker, all the Xq13.3 marker pairs were found to be in significant LD indicating that LD extends up to 9.5-11.5 Mb in the Norfolk Island population. Also, all SNPs studied showed significant LD in both Norfolk Islanders and Australian Caucasians, with two of the marker pairs in complete LD in the Norfolk population only. The Norfolk Island study population possesses a unique set of characteristics including founder effect, geographical isolation, exhaustive genealogical information and phenotypic data of use to cardiovascular disease risk traits. With LD extending up to 9.5-11 Mb, the Norfolk isolate should be a powerful resource for the localization of complex disease genes.     

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.     

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.     

Development of a single nucleotide polymorphism map of porcine chromosome 2

Single nucleotide polymorphism markers are developed on SSC2, predominantly on the p-arm. Several studies reported a quantitative trait loci (QTL) for backfat thickness in this region. Single nucleotide polymorphisms were identified by comparative re-sequencing of polymerase chain reaction (PCR) products from a panel of eight individuals. The panel consisted of five Large Whites (each from a different Dutch breeding company), a Meishan, a Pietrain and a Wild Boar. In total, 67 different PCR products were sequenced and 301 SNPs were identified in 32,429 bp of consensus sequence, an average of one SNP in every 108 bp. After correction for sample size, this polymorphism rate corresponds to a heterozygosity value of one SNP in every 357 bp. For 63% of the SNPs, there was variation among the five Large Whites, and these SNPs are relevant for linkage and association studies in commercial populations. Comparing the Whites with other breeds revealed higher variation rates with: (i) Meishan, 89%; (ii) Pietrain, 69%; (iii) Wild Boar, 70%. Because many of the experimental populations to identify QTL are based on crosses between these breeds, these SNPs are relevant for the fine mapping of the QTL identified within these crosses.     

Association of the melanocortin 4 receptor with feed intake and daily gain in F2 Mangalitsa x Piétrain pigs

The melanocortin 4 receptor (MC4R) is a key factor in the regulation of energy balance and body weight. Hence it is a candidate for feed intake and energy homeostasis-related traits. Studies in humans and swine have revealed several sequence variants in the gene that are associated with some of these traits. In pigs the coding non-synonymous missense variant Asp298Asn in MC4R has been associated with feed intake, fatness and growth. Here we confirm the association of this Piétrain-derived polymorphism with feed intake and daily gain in the F2 generation of a Mangalitsa x Piétrain cross. In one Piétrain founder animal, we detected an additional non-synonymous missense variant Arg236His. Thus, the MC4R gene could be a useful marker for increased growth in the relatively slow-growing Piétrain breed.     

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.     

Linkage disequilibrium and persistence of phase in Holstein-Friesian, Jersey and Angus cattle

When a genetic marker and a quantitative trait locus (QTL) are in linkage disequilibrium (LD) in one population, they may not be in LD in another population or their LD phase may be reversed. The objectives of this study were to compare the extent of LD and the persistence of LD phase across multiple cattle populations. LD measures r and r(2) were calculated for syntenic marker pairs using genomewide single-nucleotide polymorphisms (SNP) that were genotyped in Dutch and Australian Holstein-Friesian (HF) bulls, Australian Angus cattle, and New Zealand Friesian and Jersey cows. Average r(2) was approximately 0.35, 0.25, 0.22, 0.14, and 0.06 at marker distances 10, 20, 40, 100, and 1000 kb, respectively, which indicates that genomic selection within cattle breeds with r(2) >or= 0.20 between adjacent markers would require approximately 50,000 SNPs. The correlation of r values between populations for the same marker pairs was close to 1 for pairs of very close markers (<10 kb) and decreased with increasing marker distance and the extent of divergence between the populations. To find markers that are in LD with QTL across diverged breeds, such as HF, Jersey, and Angus, would require approximately 300,000 markers.     

Genetic relationships among Vietnamese local pigs investigated using genome‐wide SNP markers

Vietnam is one of the most important countries for pig domestication, and a total of 26 local breeds have been reported. In the present study, genetic relationships among the various pig breeds were investigated using 90 samples collected from local pigs (15 breeds) in 15 distantly separated, distinct areas of the country and six samples from Landrace pigs in Hanoi as an out‐group of a common Western breed. All samples were genotyped using the Illumina Porcine SNP60 v2 Genotyping BeadChip. We used 15 160–15 217 SNPs that showed a high degree of polymorphism in the Vietnamese breeds for identifying genetic relationships among the Vietnamese breeds. Principal components analysis showed that most pigs indigenous to Vietnam formed clusters correlated with their original geographic locations. Some Vietnamese breeds formed a cluster that was genetically related to the Western breed Landrace, suggesting the possibility of crossbreeding. These findings will be useful for the conservation and management of Vietnamese local pig breeds.     

Linkage disequilibrium over short physical distances measured in sheep using a high‐density SNP chip

The extent of linkage disequilibrium (LD) between genetic loci has implications for both association studies and the accuracy of genomic prediction. To characterise the persistence of LD in diverse sheep breeds, two SNP genotyping platforms were used. First, existing SNP genotypes from 63 breeds obtained using the ovine SNP50 BeadChip (49 034 loci) were used to estimate LD decay in populations with contrasting levels of genetic diversity. Given the paucity of marker pairs separated by short physical distances on the SNP50 BeadChip, genotyping was subsequently performed for four breeds using the recently developed ovine HD BeadChip that assays approximately 600 000 SNPs with an average genomic spacing of 5 kb. This facilitated a highly accurate estimate of LD over short genomic distances (<30 kb) and revealed LD varies considerably between sheep breeds. Further, sheep appear to contain generally lower levels of LD than do other domestic species, likely a reflection of aspects of their past population history.     

Effective genetic markers for identifying the Escherichia coli F4ac receptor status of pigs

The F4ac receptor locus (F4acR), which encodes susceptibility or resistance to Escherichia coli diarrhoea, is inherited as an autosomal recessive monogenetic trait. F4acR is localized on pig chromosome 13 (SSC13q41–q44) near the MUC13 gene. Two flanking markers (CHCF1 and ALGA0106330) with a high linkage disequilibrium (LD) with F4acR were found to be effective for the genetic identification of F4ac‐resistant pigs in the Swiss Large White breed (one recombinant out of 2034 genotyped pigs). Three recombinant boars, one each from the Duroc, Swiss Landrace and Piétrain breeds, were genotyped with seven different markers and phenotyped by means of a microscopic adhesion test. Only ALGA0072075, CHCF1 and CHCF3 indicated the correct phenotype. To test the effect of the resistance allele on production traits, 530 Large White pigs from the national test station were investigated. A significant difference existed among the F4acR locus genotypes in the intramuscular fat content of the longissimus dorsi muscle, whereas no other production traits were influenced by the resistance allele. The frequency of the CHCF1‐C and ALGA0106330‐A alleles associated with resistance in the Swiss Large White population was 60%, which is advantageous for implementing this trait in a breeding programme to select for E. coli F4ac‐resistant animals. The selection of resistant pigs should start on the male side due to the inability of resistant sows to produce sufficient amounts of protecting antibodies in the colostrum. Selection of genetically F4ac‐resistant pigs is a sustainable and suitable alternative to decreasing animal loss and antibiotic use due to diarrhoea.     

A genome scan reveals QTL for growth, fatness, leanness and meat quality in a Duroc-Pietrain resource population

We performed a genome-wide QTL scan for production traits in a line cross between Duroc and Pietrain breeds of pigs, which included 585 F(2) progeny produced from 31 full-sib families genotyped with 106 informative microsatellites. A linkage map covering all 18 autosomes and spanning 1987 Kosambi cM was constructed. Thirty-five phenotypic traits including body weight, growth, carcass composition and meat quality traits were analysed using least square regression interval mapping. Twenty-four QTL exceeded the genome-wide significance threshold, while 47 QTL reached the suggestive threshold. These QTL were located at 28 genomic regions on 16 autosomal chromosomes and QTL in 11 regions were significant at the genome-wide level. A QTL affecting pH value in loin was detected on SSC1 between marker-interval S0312-S0113 with strong statistical support (P < 3.0 x 10(-14)); this QTL was also associated with meat colour and conductivity. QTL for carcass composition and average daily gain was also found on SSC1, suggesting multiple QTL. Seventeen genomic segments had only a single QTL that reached at least suggestive significance. Forty QTL exhibited additive inheritance whereas 31 QTL showed (over-) dominance effects. Two QTL for trait backfat thickness were detected on SSC2; a significant paternal effect was found for a QTL in the IGF2 region while another QTL in the middle of SSC2 showed Mendelian expression.     

Testing for beneficial reversal of dominance during salinity shifts in the invasive copepod Eurytemora affinis,and implications for the maintenance of genetic variation

Maintenance of genetic variation at loci under selection has profound implications for adaptation under environmental change. In temporally and spatially varying habitats, non‐neutral polymorphism could be maintained by heterozygote advantage across environments (marginal overdominance), which could be greatly increased by beneficial reversal of dominance across conditions. We tested for reversal of dominance and marginal overdominance in salinity tolerance in the saltwater‐to‐freshwater invading copepod Eurytemora affinis. We compared survival of F1 offspring generated by crossing saline and freshwater inbred lines (between‐salinity F1 crosses) relative to within‐salinity F1 crosses, across three salinities. We found evidence for both beneficial reversal of dominance and marginal overdominance in salinity tolerance. In support of reversal of dominance, survival of between‐salinity F1 crosses was not different from that of freshwater F1 crosses under freshwater conditions and saltwater F1 crosses under saltwater conditions. In support of marginal overdominance, between‐salinity F1 crosses exhibited significantly higher survival across salinities relative to both freshwater and saltwater F1 crosses. Our study provides a rare empirical example of complete beneficial reversal of dominance associated with environmental change. This mechanism might be crucial for maintaining genetic variation in salinity tolerance in E. affinis populations, allowing rapid adaptation to salinity changes during habitat invasions.     

A Genome-Wide Association Study Reveals Dominance Effects on Number of Teats in Pigs

Dominance has been suggested as one of the genetic mechanisms explaining heterosis. However, using traditional quantitative genetic methods it is difficult to obtain accurate estimates of dominance effects. With the availability of dense SNP (Single Nucleotide Polymorphism) panels, we now have new opportunities for the detection and use of dominance at individual loci. Thus, the aim of this study was to detect additive and dominance effects on number of teats (NT), specifically to investigate the importance of dominance in a Landrace-based population of pigs. In total, 1,550 animals, genotyped for 32,911 SNPs, were used in single SNP analysis. SNPs with a significant genetic effect were tested for their mode of gene action being additive, dominant or a combination. In total, 21 SNPs were associated with NT, located in three regions with additive (SSC6, 7 and 12) and one region with dominant effects (SSC4). Estimates of additive effects ranged from 0.24 to 0.29 teats. The dominance effect of the QTL located on SSC4 was negative (−0.26 teats). The additive variance of the four QTLs together explained 7.37% of the total phenotypic variance. The dominance variance of the four QTLs together explained 1.82% of the total phenotypic variance, which corresponds to one-fourth of the variance explained by additive effects. The results suggest that dominance effects play a relevant role in the genetic architecture of NT. The QTL region on SSC7 contains the most promising candidate gene: VRTN. This gene has been suggested to be related to the number of vertebrae, a trait correlated with NT.     

Mapping and fine mapping of quantitative trait loci for the number of vertebrae in a White Duroc × Chinese Erhualian intercross resource population

The number of vertebrae is associated with body size and meat production in pigs. To identify quantitative trait loci (QTL) for the number of vertebrae, phenotypic values were measured in 1029 individuals from a White Duroc × Chinese Erhualian intercross F₂ population. A whole genome scan was performed with 194 microsatellite markers in the F₂ population. Four genome‐wide significant QTL and eight chromosome‐wide significant QTL for the number of vertebrae were identified on pig chromosomes (SSC) 1, 2, 6, 7, 10 and 12. The most significant QTL was detected on SSC7 with a confidence interval of 1 cM, explaining 42.32% of the phenotypic variance in the thoracic vertebral number. The significant QTL on SSC1, 2 and 7 confirmed previous reports. A panel of 276 animals representing seven Western and Chinese breeds was genotyped with 34 microsatellite markers in the SSC7 QTL region. No obvious selective sweep effect was observed in the tested breeds, indicating that intensive selection for enlarged body size in Western commercial breeds did not wipe out the genetic variability in the QTL region. The Q alleles for increased vertebral number originated from both Chinese Erhualian and White Duroc founder animals. A haplotype block of approximately 900 kb was found to be shared by all Q‐bearing chromosomes of F₁ sires except for one distinct Q chromosome. The critical region harbours the newly reported VRTN gene associated with vertebral number. Further investigations are required to confirm whether VRTN or two other positional candidate genes, PROX2 and FOS, cause the QTL effect.     

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.