A genome-wide association study of production traits in a commercial population of Large White pigs: evidence of haplotypes affecting meat quality

Background

Numerous quantitative trait loci (QTL) have been detected in pigs over the past 20 years using microsatellite markers. However, due to the low density of these markers, the accuracy of QTL location has generally been poor. Since 2009, the dense genome coverage provided by the Illumina PorcineSNP60 BeadChip has made it possible to more accurately map QTL using genome-wide association studies (GWAS). Our objective was to perform high-density GWAS in order to identify genomic regions and corresponding haplotypes associated with production traits in a French Large White population of pigs.

Methods

Animals (385 Large White pigs from 106 sires) were genotyped using the PorcineSNP60 BeadChip and evaluated for 19 traits related to feed intake, growth, carcass composition and meat quality. Of the 64 432 SNPs on the chip, 44 412 were used for GWAS with an animal mixed model that included a regression coefficient for the tested SNPs and a genomic kinship matrix. SNP haplotype effects in QTL regions were then tested for association with phenotypes following phase reconstruction based on the Sscrofa10.2 pig genome assembly.

Results

Twenty-three QTL regions were identified on autosomes and their effects ranged from 0.25 to 0.75 phenotypic standard deviation units for feed intake and feed efficiency (four QTL), carcass (12 QTL) and meat quality traits (seven QTL). The 10 most significant QTL regions had effects on carcass (chromosomes 7, 10, 16, 17 and 18) and meat quality traits (two regions on chromosome 1 and one region on chromosomes 8, 9 and 13). Thirteen of the 23 QTL regions had not been previously described. A haplotype block of 183 kb on chromosome 1 (six SNPs) was identified and displayed three distinct haplotypes with significant (0.0001 < P < 0.03) associations with all evaluated meat quality traits.

Conclusions

GWAS analyses with the PorcineSNP60 BeadChip enabled the detection of 23 QTL regions that affect feed consumption, carcass and meat quality traits in a LW population, of which 13 were novel QTL. The proportionally larger number of QTL found for meat quality traits suggests a specific opportunity for improving these traits in the pig by genomic selection.     

A search for quantitative trait loci for milk production traits on chromosome 6 in Finnish Ayrshire cattle

Cattle chromosome 6 was scanned with 11 markers, ten microsatellites and the casein haplotype, to identify quantitative trait loci (QTLs) affecting the following milk production traits: milk yield, fat percentage, fat yield, protein percentage and protein yield. Twelve Finnish Ayrshire half-sib families with a total of 480 sons were genotyped and used in a grand-daughter design. Interval mapping was performed with a multiple-marker regression approach with a one-QTL and a two-QTL model, and the significance threshold values were determined empirically using a permutation test. Across-family analysis with the one-QTL model revealed an effect on protein percentage (P < 0.05) and on milk yield (P < 0.05). The analysis with the two-QTL model identified significant effects (P < 0.05) on protein percentage, milk yield, and fat yield. Comparing these two cases, the results suggest the existence of two QTLs on chromosome 6 with an effect on milk production traits. One of the QTLs was located around the casein genes. As the other QTL was similar in location and effect to a QTL found previously in Holstein-Friesians, an identity-by-descent approach could be applied to fine map this region.     

Detecting parent of origin and dominant QTL in a two-generation commercial poultry pedigree using variance component methodology

Introduction

Variance component QTL methodology was used to analyse three candidate regions on chicken chromosomes 1, 4 and 5 for dominant and parent-of-origin QTL effects. Data were available for bodyweight and conformation score measured at 40 days from a two-generation commercial broiler dam line. One hundred dams were nested in 46 sires with phenotypes and genotypes on 2708 offspring. Linear models were constructed to simultaneously estimate fixed, polygenic and QTL effects. Different genetic models were compared using likelihood ratio test statistics derived from the comparison of full with reduced or null models. Empirical thresholds were derived by permutation analysis.

Results

Dominant QTL were found for bodyweight on chicken chromosome 4 and for bodyweight and conformation score on chicken chromosome 5. Suggestive evidence for a maternally expressed QTL for bodyweight and conformation score was found on chromosome 1 in a region corresponding to orthologous imprinted regions in the human and mouse.

Conclusion

Initial results suggest that variance component analysis can be applied within commercial populations for the direct detection of segregating dominant and parent of origin effects.     

Number and mode of inheritance of QTL influencing backfat thickness on SSC2p in Sino-European pig pedigrees

Background

In the pig, multiple QTL associated with growth and fatness traits have been mapped to chromosome 2 (SSC2) and among these, at least one shows paternal expression due to the IGF2-intron3-G3072A substitution. Previously published results on the position and imprinting status of this QTL disagree between analyses from French and Dutch F2 crossbred pig populations obtained with the same breeds (Meishan crossed with Large White or Landrace).

Methods

To study the role of paternal and maternal alleles at the IGF2 locus and to test the hypothesis of a second QTL affecting backfat thickness on the short arm of SSC2 (SSC2p), a QTL mapping analysis was carried out on a combined pedigree including both the French and Dutch F2 populations, on the progeny of F1 males that were heterozygous (A/G) and homozygous (G/G) at the IGF2 locus. Simulations were performed to clarify the relations between the two QTL and to understand to what extent they can explain the discrepancies previously reported.

Results

The QTL analyses showed the segregation of at least two QTL on chromosome 2 in both pedigrees, i.e. the IGF2 locus and a second QTL segregating at least in the G/G F1 males and located between positions 30 and 51 cM. Statistical analyses highlighted that the maternally inherited allele at the IGF2 locus had a significant effect but simulation studies showed that this is probably a spurious effect due to the segregation of the second QTL.

Conclusions

Our results show that two QTL on SSC2p affect backfat thickness. Differences in the pedigree structures and in the number of heterozygous females at the IGF2 locus result in different imprinting statuses in the two pedigrees studied. The spurious effect observed when a maternally allele is present at the IGF2 locus, is in fact due to the presence of a second closely located QTL. This work confirms that pig chromosome 2 is a major region associated with fattening traits.     

Genome-wide QTL mapping of nine body composition and bone mineral density traits in pigs

Background

Since the pig is one of the most important livestock animals worldwide, mapping loci that are associated with economically important traits and/or traits that influence animal welfare is extremely relevant for efficient future pig breeding. Therefore, the purpose of this study was a genome-wide mapping of quantitative trait loci (QTL) associated with nine body composition and bone mineral traits: absolute (Fat, Lean) and percentage (FatPC, LeanPC) fat and lean mass, live weight (Weight), soft tissue X-ray attenuation coefficient (R), absolute (BMC) and percentage (BMCPC) bone mineral content and bone mineral density (BMD).

Methods

Data on the nine traits investigated were obtained by Dual-energy X-ray absorptiometry for 551 pigs that were between 160 and 200 days old. In addition, all pigs were genotyped using Illumina’s PorcineSNP60 Genotyping BeadChip. Based on these data, a genome-wide combined linkage and linkage disequilibrium analysis was conducted. Thus, we used 44 611 sliding windows that each consisted of 20 adjacent single nucleotide polymorphisms (SNPs). For the middle of each sliding window a variance component analysis was carried out using ASReml. The underlying mixed linear model included random QTL and polygenic effects, with fixed effects of sex, housing, season and age.

Results

Using a Bonferroni-corrected genome-wide significance threshold of P < 0.001, significant peaks were identified for all traits except BMCPC. Overall, we identified 72 QTL on 16 chromosomes, of which 24 were significantly associated with one trait only and the remaining with more than one trait. For example, a QTL on chromosome 2 included the highest peak across the genome for four traits (Fat, FatPC, LeanPC and R). The nearby gene, ZNF608, is known to be associated with body mass index in humans and involved in starvation in Drosophila, which makes it an extremely good candidate gene for this QTL.

Conclusions

Our QTL mapping approach identified 72 QTL, some of which confirmed results of previous studies in pigs. However, we also detected significant associations that have not been published before and were able to identify a number of new and promising candidate genes, such as ZNF608.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0068-2) contains supplementary material, which is available to authorized users.     

Estimation of (co)variances for genomic regions of flexible sizes: application to complex infectious udder diseases in dairy cattle

Background

Multi-trait genomic models in a Bayesian context can be used to estimate genomic (co)variances, either for a complete genome or for genomic regions (e.g. per chromosome) for the purpose of multi-trait genomic selection or to gain further insight into the genomic architecture of related traits such as mammary disease traits in dairy cattle.

Methods

Data on progeny means of six traits related to mastitis resistance in dairy cattle (general mastitis resistance and five pathogen-specific mastitis resistance traits) were analyzed using a bivariate Bayesian SNP-based genomic model with a common prior distribution for the marker allele substitution effects and estimation of the hyperparameters in this prior distribution from the progeny means data. From the Markov chain Monte Carlo samples of the allele substitution effects, genomic (co)variances were calculated on a whole-genome level, per chromosome, and in regions of 100 SNP on a chromosome.

Results

Genomic proportions of the total variance differed between traits. Genomic correlations were lower than pedigree-based genetic correlations and they were highest between general mastitis and pathogen-specific traits because of the part-whole relationship between these traits. The chromosome-wise genomic proportions of the total variance differed between traits, with some chromosomes explaining higher or lower values than expected in relation to chromosome size. Few chromosomes showed pleiotropic effects and only chromosome 19 had a clear effect on all traits, indicating the presence of QTL with a general effect on mastitis resistance. The region-wise patterns of genomic variances differed between traits. Peaks indicating QTL were identified but were not very distinctive because a common prior for the marker effects was used. There was a clear difference in the region-wise patterns of genomic correlation among combinations of traits, with distinctive peaks indicating the presence of pleiotropic QTL.

Conclusions

The results show that it is possible to estimate, genome-wide and region-wise genomic (co)variances of mastitis resistance traits in dairy cattle using multivariate genomic models.     

Overlapping mouse subcongenic strains successfully separate two linked body fat QTL on distal MMU 2

Background

Mouse chromosome 2 is linked to growth and body fat phenotypes in many mouse crosses. With the goal to identify the underlying genes regulating growth and body fat on mouse chromosome 2, we developed five overlapping subcongenic strains that contained CAST/EiJ donor regions in a C57BL/6J^hg/hg background (hg is a spontaneous deletion of 500 Kb on mouse chromosome 10). To fine map QTL on distal mouse chromosome 2 a total of 1,712 F2 mice from the five subcongenic strains, plus 278 F2 mice from the HG2D founder congenic strain were phenotyped and analyzed. Interval mapping (IM) and composite IM (CIM) were performed on body weight and body fat traits on a combination of SNP and microsatellite markers, which generated a high-density genotyping panel.

Results

Phenotypic analysis and interval mapping of total fat mass identified two QTL on distal mouse chromosome 2. One QTL between 150 and 161 Mb, Fatq2a, and the second between 173.3 and 175.6 Mb, Fatq2b. The two QTL reside in different congenic strains with significant total fat differences between homozygous cast/cast and b6/b6 littermates. Both of these QTL were previously identified only as a single QTL affecting body fat, Fatq2. Furthermore, through a novel approach referred here as replicated CIM, Fatq2b was mapped to the Gnas imprinted locus.

Conclusions

The integration of subcongenic strains, high-density genotyping, and CIM succesfully partitioned two previously linked QTL 20 Mb apart, and the strongest QTL, Fatq2b, was fine mapped to a ~2.3 Mb region interval encompassing the Gnas imprinted locus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-014-1191-8) contains supplementary material, which is available to authorized users.     

Genome-wide QTL analysis of meat quality-related traits in a large F2 intercross between Landrace and Korean native pigs

Background

We conducted a genome-wide linkage analysis to identify quantitative trait loci (QTL) that influence meat quality-related traits in a large F₂ intercross between Landrace and Korean native pigs. Thirteen meat quality-related traits of the m. longissimus lumborum et thoracis were measured in more than 830 F₂ progeny. All these animals were genotyped with 173 microsatellite markers located throughout the pig genome, and the GridQTL program based on the least squares regression model was used to perform the QTL analysis.

Results

We identified 23 genome-wide significant QTL in eight chromosome regions (SSC1, 2, 6, 7, 9, 12, 13, and 16) (SSC for Sus Scrofa) and detected 51 suggestive QTL in the 17 chromosome regions. QTL that affect 10 meat quality traits were detected on SSC12 and were highly significant at the genome-wide level. In particular, the QTL with the largest effect affected crude fat percentage and explained 22.5% of the phenotypic variance (F-ratio = 278.0 under the additive model, nominal P = 5.5 × 10⁻⁵⁵). Interestingly, the QTL on SSC12 that influenced meat quality traits showed an obvious trend for co-localization.

Conclusions

Our results confirm several previously reported QTL. In addition, we identified novel QTL for meat quality traits, which together with the associated positional candidate genes improve the knowledge on the genetic structure that underlies genetic variation for meat quality traits in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0080-6) contains supplementary material, which is available to authorized users.     

Genetic variation in the pleiotropic association between physical activity and body weight in mice

Background

A sedentary lifestyle is often assumed to lead to increases in body weight and potentially obesity and related diseases but in fact little is known about the genetic association between physical activity and body weight. We tested for such an association between body weight and the distance, duration, and speed voluntarily run by 310 mice from the F₂ generation produced from an intercross of two inbred lines that differed dramatically in their physical activity levels.

Methods

We used a conventional interval mapping approach with SNP markers to search for QTLs that affected both body weight and activity traits. We also conducted a genome scan to search for relationship QTLs (relQTLs), or chromosomal regions that affected an activity trait variably depending on the phenotypic value of body weight.

Results

We uncovered seven quantitative trait loci (QTLs) affecting body weight, but only one co-localized with another QTL previously found for activity traits. We discovered 19 relQTLs that provided evidence for a genetic (pleiotropic) association of physical activity and body weight. The three genotypes at each of these loci typically exhibited a combination of negative, zero, and positive regressions of the activity traits on body weight, the net effect of which was to produce overall independence of body weight from physical activity. We also demonstrated that the relQTLs produced these varying associations through differential epistatic interactions with a number of other epistatic QTLs throughout the genome.

Conclusion

It was concluded that individuals with specific combinations of genotypes at the relQTLs and epiQTLs might account for some of the variation typically seen in plots of the association of physical activity with body weight.     

Selection in action: dissecting the molecular underpinnings of the increasing muscle mass of Belgian Blue Cattle

Background

Belgian Blue cattle are famous for their exceptional muscular development or “double-muscling”. This defining feature emerged following the fixation of a loss-of-function variant in the myostatin gene in the eighties. Since then, sustained selection has further increased muscle mass of Belgian Blue animals to a comparable extent. In the present paper, we study the genetic determinants of this second wave of muscle growth.

Results

A scan for selective sweeps did not reveal the recent fixation of another allele with major effect on muscularity. However, a genome-wide association study identified two genome-wide significant and three suggestive quantitative trait loci (QTL) affecting specific muscle groups and jointly explaining 8-21% of the heritability. The top two QTL are caused by presumably recent mutations on unique haplotypes that have rapidly risen in frequency in the population. While one appears on its way to fixation, the ascent of the other is compromised as the likely underlying MRC2 mutation causes crooked tail syndrome in homozygotes. Genomic prediction models indicate that the residual additive variance is largely polygenic.

Conclusions

Contrary to complex traits in humans which have a near-exclusive polygenic architecture, muscle mass in beef cattle (as other production traits under directional selection), appears to be controlled by (i) a handful of recent mutations with large effect that rapidly sweep through the population, and (ii) a large number of presumably older variants with very small effects that rise slowly in the population (polygenic adaptation).

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-796) contains supplementary material, which is available to authorized users.     

Persistence of accuracy of genome-wide breeding values over generations when including a polygenic effect

Background

When estimating marker effects in genomic selection, estimates of marker effects may simply act as a proxy for pedigree, i.e. their effect may partially be attributed to their association with superior parents and not be linked to any causative QTL. Hence, these markers mainly explain polygenic effects rather than QTL effects. However, if a polygenic effect is included in a Bayesian model, it is expected that the estimated effect of these markers will be more persistent over generations without having to re-estimate the marker effects every generation and will result in increased accuracy and reduced bias.

Methods

Genomic selection using the Bayesian method, ''BayesB'' was evaluated for different marker densities when a polygenic effect is included (GWpEBV) and not included (GWEBV) in the model. Linkage disequilibrium and a mutation drift balance were obtained by simulating a population with a Ne of 100 over 1,000 generations.

Results

Accuracy of selection was slightly higher for the model including a polygenic effect than for the model not including a polygenic effect whatever the marker density. The accuracy decreased in later generations, and this reduction was stronger for lower marker densities. However, no significant difference in accuracy was observed between the two models. The linear regression of TBV on GWEBV and GWpEBV was used as a measure of bias. The regression coefficient was more stable over generations when a polygenic effect was included in the model, and was always between 0.98 and 1.00 for the highest marker density. The regression coefficient decreased more quickly with decreasing marker density.

Conclusions

Including a polygenic effect had no impact on the selection accuracy, but showed reduced bias, which is especially important when estimates of genome-wide markers are used to estimate breeding values over more than one generation.     

Mapping quantitative trait loci (QTL) in sheep. IV. Analysis of lactation persistency and extended lactation traits in sheep

Background

In sheep dairy production, total lactation performance, and length of lactation of lactation are of economic significance. A more persistent lactation has been associated with improved udder health. An extended lactation is defined by a longer period of milkability. This study is the first investigation to examine the presence of quantitative trait loci (QTL) for extended lactation and lactation persistency in sheep.

Methods

An (Awassi × Merino) × Merino single-sire backcross family with 172 ewes was used to map QTL for lactation persistency and extended lactation traits on a framework map of 189 loci across all autosomes. The Wood model was fitted to data from multiple lactations to estimate parameters of ovine lactation curves, and these estimates were used to derive measures of lactation persistency and extended lactation traits of milk, protein, fat, lactose, useful yield, and somatic cell score. These derived traits were subjected to QTL analyses using maximum likelihood estimation and regression analysis.

Results

Overall, one highly significant (LOD > 3.0), four significant (2.0 < LOD < 3.0) and five suggestive (1.7 < LOD < 2.0) QTL were detected across all traits in common by both mapping methods. One additional suggestive QTL was identified using maximum likelihood estimation, and four suggestive (0.01 < P < 0.05) and two significant (P < 0.01) QTL using the regression approach only. All detected QTL had effect sizes in the range of 0.48 to 0.64 SD, corresponding to QTL heritabilities of 3.1 to 8.9%. The comparison of the detected QTL with results in cattle showed conserved linkage regions. Most of the QTL identified for lactation persistency and extended lactation did not coincide. This suggests that persistency and extended lactation for the same as well as different milk yield and component traits are not controlled by the same genes.

Conclusion

This study identified ten novel QTL for lactation persistency and extended lactation in sheep, but results suggest that lactation persistency and extended lactation do not have a major gene in common. These results provide a basis for further validation in extended families and other breeds as well as targeting regions for genome-wide association mapping using high-density SNP arrays.     

Genome-wide mapping of Quantitative Trait Loci for fatness,fat cell characteristics and fat metabolism in three porcine F2 crosses

Background

QTL affecting fat deposition related performance traits have been considered in several studies and mapped on numerous porcine chromosomes. However, activity of specific enzymes, protein content and cell structure in fat tissue probably depend on a smaller number of genes than traits related to fat content in carcass. Thus, in this work traits related to metabolic and cytological features of back fat tissue and fat related performance traits were investigated in a genome-wide QTL analysis. QTL similarities and differences were examined between three F₂ crosses, and between male and female animals.

Methods

A total of 966 F₂ animals originating from crosses between Meishan (M), Pietrain (P) and European wild boar (W) were analysed for traits related to fat performance (11), enzymatic activity (9) and number and volume of fat cells (20). Per cross, 216 (M × P), 169 (W × P) and 195 (W × M) genome-wide distributed marker loci were genotyped. QTL mapping was performed separately for each cross in steps of 1 cM and steps were reduced when the distance between loci was shorter. The additive and dominant components of QTL positions were detected stepwise by using a multiple position model.

Results

A total of 147 genome-wide significant QTL (76 at P < 0.05 and 71 at P < 0.01) were detected for the three crosses. Most of the QTL were identified on SSC1 (between 76-78 and 87-90 cM), SSC7 (predominantly in the MHC region) and SSCX (in the vicinity of the gene CAPN6). Additional genome-wide significant QTL were found on SSC8, 12, 13, 14, 16, and 18. In many cases, the QTL are mainly additive and differ between F₂ crosses. Many of the QTL profiles possess multiple peaks especially in regions with a high marker density. Sex specific analyses, performed for example on SSC6, SSC7 and SSCX, show that for some traits the positions differ between male and female animals. For the selected traits, the additive and dominant components that were analysed for QTL positions on different chromosomes, explain in combination up to 23% of the total trait variance.

Conclusions

Our results reveal specific and partly new QTL positions across genetically diverse pig crosses. For some of the traits associated with specific enzymes, protein content and cell structure in fat tissue, it is the first time that they are included in a QTL analysis. They provide large-scale information to analyse causative genes and useful data for the pig industry.     

Combined detection and introgression of QTL in outbred populations

Background

Detecting a QTL is only the first step in genetic improvement programs. When a QTL with desirable characteristics is found, e.g. in a wild or unimproved population, it may be interesting to introgress the detected QTL into the commercial population. One approach to shorten the time needed for introgression is to combine both QTL identification and introgression, into a single step. This combines the strengths of fine mapping and backcrossing and paves the way for introgression of desirable but unknown QTL into recipient animal and plant lines.

Methods

The method consisting in combining QTL mapping and gene introgression has been extended from inbred to outbred populations in which QTL allele frequencies vary both in recipient and donor lines in different scenarios and for which polygenic effects are included in order to model background genes. The effectiveness of the combined QTL detection and introgression procedure was evaluated by simulation through four backcross generations.

Results

The allele substitution effect is underestimated when the favourable QTL allele is not fixed in the donor line. This underestimation is proportional to the frequency differences of the favourable QTL allele between the lines. In most scenarios, the estimates of the QTL location are unbiased and accurate. The retained donor chromosome segment and linkage drag are similar to expected values from other published studies.

Conclusions

In general, our results show that it is possible to combine QTL detection and introgression even in outbred species. Separating QTL mapping and introgression processes is often thought to be longer and more costly. However, using a combined process saves at least one generation. With respect to the linkage drag and obligatory drag, the results of the combined detection and introgression scheme are very similar to those of traditional introgression schemes.     

Modelling and estimation of genotype by environment interactions for production traits in French dairy cattle

Background

Genotype by environment interactions are currently ignored in national genetic evaluations of dairy cattle. However, this is often questioned, especially when environment or herd management is wide-ranging. The aim of this study was to assess genotype by environment interactions for production traits (milk, protein, fat yields and fat and protein contents) in French dairy cattle using an original approach to characterize the environments.

Methods

Genetic parameters of production traits were estimated for three breeds (Holstein, Normande and Montbéliarde) using multiple-trait and reaction norm models. Variables derived from Herd Test Day profiles obtained after a test day model evaluation were used to define herd environment.

Results

Multiple-trait and reaction norm models gave similar results. Genetic correlations were very close to unity for all traits, except between some extreme environments. However, a relatively wide range of heritabilities by trait and breed was found across environments. This was more the case for milk, protein and fat yields than for protein and fat contents.

Conclusions

No real reranking of animals was observed across environments. However, a significant scale effect exists: the more intensive the herd management for milk yield, the larger the heritability.     

Epistatic QTL pairs associated with meat quality and carcass composition traits in a porcine Duroc × Pietrain population

Background

Quantitative trait loci (QTL) analyses in pig have revealed numerous individual QTL affecting growth, carcass composition, reproduction and meat quality, indicating a complex genetic architecture. In general, statistical QTL models consider only additive and dominance effects and identification of epistatic effects in livestock is not yet widespread. The aim of this study was to identify and characterize epistatic effects between common and novel QTL regions for carcass composition and meat quality traits in pig.

Methods

Five hundred and eighty five F₂ pigs from a Duroc × Pietrain resource population were genotyped using 131 genetic markers (microsatellites and SNP) spread over the 18 pig autosomes. Phenotypic information for 26 carcass composition and meat quality traits was available for all F₂ animals. Linkage analysis was performed in a two-step procedure using a maximum likelihood approach implemented in the QxPak program.

Results

A number of interacting QTL was observed for different traits, leading to the identification of a variety of networks among chromosomal regions throughout the porcine genome. We distinguished 17 epistatic QTL pairs for carcass composition and 39 for meat quality traits. These interacting QTL pairs explained up to 8% of the phenotypic variance.

Conclusions

Our findings demonstrate the significance of epistasis in pigs. We have revealed evidence for epistatic relationships between different chromosomal regions, confirmed known QTL loci and connected regions reported in other studies. Considering interactions between loci allowed us to identify several novel QTL and trait-specific relationships of loci within and across chromosomes.     

QTL fine mapping with Bayes C(π): a simulation study

Background

Accurate QTL mapping is a prerequisite in the search for causative mutations. Bayesian genomic selection models that analyse many markers simultaneously should provide more accurate QTL detection results than single-marker models. Our objectives were to (a) evaluate by simulation the influence of heritability, number of QTL and number of records on the accuracy of QTL mapping with Bayes Cπ and Bayes C; (b) estimate the QTL status (homozygous vs. heterozygous) of the individuals analysed. This study focussed on the ten largest detected QTL, assuming they are candidates for further characterization.

Methods

Our simulations were based on a true dairy cattle population genotyped for 38 277 phased markers. Some of these markers were considered biallelic QTL and used to generate corresponding phenotypes. Different numbers of records (4387 and 1500), heritability values (0.1, 0.4 and 0.7) and numbers of QTL (10, 100 and 1000) were studied. QTL detection was based on the posterior inclusion probability for individual markers, or on the sum of the posterior inclusion probabilities for consecutive markers, estimated using Bayes C or Bayes Cπ. The QTL status of the individuals was derived from the contrast between the sums of the SNP allelic effects of their chromosomal segments.

Results

The proportion of markers with null effect (π) frequently did not reach convergence, leading to poor results for Bayes Cπ in QTL detection. Fixing π led to better results. Detection of the largest QTL was most accurate for medium to high heritability, for low to moderate numbers of QTL, and with a large number of records. The QTL status was accurately inferred when the distribution of the contrast between chromosomal segment effects was bimodal.

Conclusions

QTL detection is feasible with Bayes C. For QTL detection, it is recommended to use a large dataset and to focus on highly heritable traits and on the largest QTL. QTL statuses were inferred based on the distribution of the contrast between chromosomal segment effects.     

Joint QTL analysis of three connected F2-crosses in pigs

Background

Numerous QTL mapping resource populations are available in livestock species. Usually they are analysed separately, although the same founder breeds are often used. The aim of the present study was to show the strength of analysing F₂-crosses jointly in pig breeding when the founder breeds of several F₂-crosses are the same.

Methods

Three porcine F₂-crosses were generated from three founder breeds (i.e. Meishan, Pietrain and wild boar). The crosses were analysed jointly, using a flexible genetic model that estimated an additive QTL effect for each founder breed allele and a dominant QTL effect for each combination of alleles derived from different founder breeds. The following traits were analysed: daily gain, back fat and carcass weight. Substantial phenotypic variation was observed within and between crosses. Multiple QTL, multiple QTL alleles and imprinting effects were considered. The results were compared to those obtained when each cross was analysed separately.

Results

For daily gain, back fat and carcass weight, 13, 15 and 16 QTL were found, respectively. For back fat, daily gain and carcass weight, respectively three, four, and five loci showed significant imprinting effects. The number of QTL mapped was much higher than when each design was analysed individually. Additionally, the test statistic plot along the chromosomes was much sharper leading to smaller QTL confidence intervals. In many cases, three QTL alleles were observed.

Conclusions

The present study showed the strength of analysing three connected F₂-crosses jointly. In this experiment, statistical power was high because of the reduced number of estimated parameters and the large number of individuals. The applied model was flexible and was computationally fast.     

Genome-wide association analyses for meat quality traits in Chinese Erhualian pigs and a Western Duroc?×?(Landrace × Yorkshire) commercial population

Background

Understanding the genetic mechanisms that underlie meat quality traits is essential to improve pork quality. To date, most quantitative trait loci (QTL) analyses have been performed on F₂ crosses between outbred pig strains and have led to the identification of numerous QTL. However, because linkage disequilibrium is high in such crosses, QTL mapping precision is unsatisfactory and only a few QTL have been found to segregate within outbred strains, which limits their use to improve animal performance. To detect QTL in outbred pig populations of Chinese and Western origins, we performed genome-wide association studies (GWAS) for meat quality traits in Chinese purebred Erhualian pigs and a Western Duroc × (Landrace × Yorkshire) (DLY) commercial population.

Methods

Three hundred and thirty six Chinese Erhualian and 610 DLY pigs were genotyped using the Illumina PorcineSNP60K Beadchip and evaluated for 20 meat quality traits. After quality control, 35 985 and 56 216 single nucleotide polymorphisms (SNPs) were available for the Chinese Erhualian and DLY datasets, respectively, and were used to perform two separate GWAS. We also performed a meta-analysis that combined P-values and effects of 29 516 SNPs that were common to Erhualian, DLY, F₂ and Sutai pig populations.

Results

We detected 28 and nine suggestive SNPs that surpassed the significance level for meat quality in Erhualian and DLY pigs, respectively. Among these SNPs, ss131261254 on pig chromosome 4 (SSC4) was the most significant (P = 7.97E-09) and was associated with drip loss in Erhualian pigs. Our results suggested that at least two QTL on SSC12 and on SSC15 may have pleiotropic effects on several related traits. All the QTL that were detected by GWAS were population-specific, including 12 novel regions. However, the meta-analysis revealed seven novel QTL for meat characteristics, which suggests the existence of common underlying variants that may differ in frequency across populations. These QTL regions contain several relevant candidate genes.

Conclusions

These findings provide valuable insights into the molecular basis of convergent evolution of meat quality traits in Chinese and Western breeds that show divergent phenotypes. They may contribute to genetic improvement of purebreds for crossbred performance.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0120-x) contains supplementary material, which is available to authorized users.     

Accounting for trait architecture in genomic predictions of US Holstein cattle using a weighted realized relationship matrix

Background

Genomic BLUP (GBLUP) can predict breeding values for non-phenotyped individuals based on the identity-by-state genomic relationship matrix (G). The G matrix can be constructed from thousands of markers spread across the genome. The strongest assumption of G and consequently of GBLUP is that all markers contribute equally to the genetic variance of a trait. This assumption is violated for traits that are controlled by a small number of quantitative trait loci (QTL) or individual QTL with large effects. In this paper, we investigate the performance of using a weighted genomic relationship matrix (wG) that takes into consideration the genetic architecture of the trait in order to improve predictive ability for a wide range of traits. Multiple methods were used to calculate weights for several economically relevant traits in US Holstein dairy cattle. Predictive performance was tested by k-means cross-validation.

Results

Relaxing the GBLUP assumption of equal marker contribution by increasing the weight that is given to a specific marker in the construction of the trait-specific G resulted in increased predictive performance. The increase was strongest for traits that are controlled by a small number of QTL (e.g. fat and protein percentage). Furthermore, bias in prediction estimates was reduced compared to that resulting from the use of regular G. Even for traits with low heritability and lower general predictive performance (e.g. calving ease traits), weighted G still yielded a gain in accuracy.

Conclusions

Genomic relationship matrices weighted by marker realized variance yielded more accurate and less biased predictions for traits regulated by few QTL. Genome-wide association analyses were used to derive marker weights for creating weighted genomic relationship matrices. However, this can be cumbersome and prone to low stability over generations because of erosion of linkage disequilibrium between markers and QTL. Future studies may include other sources of information, such as functional annotation and gene networks, to better exploit the genetic architecture of traits and produce more stable predictions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0100-1) contains supplementary material, which is available to authorized users.