RESUMEN

BACKGROUND: In this study, we tested whether genotyping both live and dead animals (GSD) realises more genetic gain for post-weaning survival (PWS) in pigs compared to genotyping only live animals (GOS). METHODS: Stochastic simulation was used to estimate the rate of genetic gain realised by GSD and GOS at a 0.01 rate of pedigree-based inbreeding in three breeding schemes, which differed in PWS (95%, 90% and 50%) and litter size (6 and 10). Pedigree-based selection was conducted as a point of reference. Variance components were estimated and then estimated breeding values (EBV) were obtained in each breeding scheme using a linear or a threshold model. Selection was for a single trait, i.e. PWS with a heritability of 0.02 on the observed scale. The trait was simulated on the underlying scale and was recorded as binary (0/1). Selection candidates were genotyped and phenotyped before selection, with only live candidates eligible for selection. Genotyping strategies differed in the proportion of live and dead animals genotyped, but the phenotypes of all animals were used for predicting EBV of the selection candidates. RESULTS: Based on a 0.01 rate of pedigree-based inbreeding, GSD realised 14 to 33% more genetic gain than GOS for all breeding schemes depending on PWS and litter size. GSD increased the prediction accuracy of EBV for PWS by at least 14% compared to GOS. The use of a linear versus a threshold model did not have an impact on genetic gain for PWS regardless of the genotyping strategy and the bias of the EBV did not differ significantly among genotyping strategies. CONCLUSIONS: Genotyping both dead and live animals was more informative than genotyping only live animals to predict the EBV for PWS of selection candidates, but with marginal increases in genetic gain when the proportion of dead animals genotyped was 60% or greater. Therefore, it would be worthwhile to use genomic information on both live and more than 20% dead animals to compute EBV for the genetic improvement of PWS under the assumption that dead animals reflect increased liability on the underlying scale.

Asunto(s)

RESUMEN

BACKGROUND: The objective of this study was to introduce a genome-wide association study (GWAS) in conjunction with segregation analysis on monogenic categorical traits. Genotype probabilities calculated from phenotypes, mode of inheritance and pedigree information, are expressed as the expected allele count (EAC) (range 0 to 2), and are inherited additively, by definition, unlike the original phenotypes, which are non-additive and could be of incomplete penetrance. The EAC are regressed on the single nucleotide polymorphism (SNP) genotypes, similar to an additive GWAS. In this study, horn phenotypes in Merino sheep are used to illustrate the advantages of using the segregation GWAS, a trait believed to be monogenic, affected by dominance, sex-dependent expression and likely affected by incomplete penetrance. We also used simulation to investigate whether incomplete penetrance can cause prediction errors in Merino sheep for horn status. RESULTS: Estimated penetrance values differed between the sexes, where males showed almost complete penetrance, especially for horned and polled phenotypes, while females had low penetrance values for the horned status. This suggests that females homozygous for the 'horned allele' have a horned phenotype in only 22% of the cases while 78% will be knobbed or have scurs. The GWAS using EAC on 4001 animals and 510,174 SNP genotypes from the Illumina Ovine high-density (600k) chip gave a stronger association compared to using actual phenotypes. The correlation between the EAC and the allele count of the SNP with the highest -log10(p-value) was 0.73 in males and 0.67 in females. Simulations using penetrance values found by the segregation analyses resulted in higher correlations between the EAC and the causative mutation (0.95 for males and 0.89 for females, respectively), suggesting that the most predictive SNP is not in full LD with the causative mutation. CONCLUSIONS: Our results show clear differences in penetrance values between males and female Merino sheep for horn status. Segregation analysis for a trait with mutually exclusive phenotypes, non-additive inheritance, and/or incomplete penetrance can lead to considerably more power in a GWAS because the linearized genotype probabilities are additive and can accommodate incomplete penetrance. This method can be extended to any monogenic controlled categorical trait of which the phenotypes are mutually exclusive.

Asunto(s)

RESUMEN

The premise was tested that the additional genetic gain was achieved in the overall breeding objective in a pig breeding program using genomic selection (GS) compared to a conventional breeding program, however, some traits achieved larger gain than other traits. GS scenarios based on different reference population sizes were evaluated. The scenarios were compared using a deterministic simulation model to predict genetic gain in scenarios with and without using genomic information as an additional information source. All scenarios were compared based on selection accuracy and predicted genetic gain per round of selection for objective traits in both sire and dam lines. The results showed that GS scenarios increased overall response in the breeding objectives by 9% to 56% and 3.5% to 27% in the dam and sire lines, respectively. The difference in response resulted from differences in the size of the reference population. Although all traits achieved higher selection accuracy in GS, traits with limited phenotypic information at the time of selection or with low heritability, such as sow longevity, number of piglets born alive, pre- and post-weaning survival, as well as meat and carcass quality traits achieved the largest additional response. This additional response came at the expense of smaller responses for traits that are easy to measure, such as back fat and average daily gain in GS compared to the conventional breeding program. Sow longevity and drip loss percentage did not change in a favourable direction in GS with a reference population of 500 pigs. With a reference population of 1000 pigs or onwards, sow longevity and drip loss percentage began to change in a favourable direction. Despite the smaller responses for average daily gain and back fat thickness in GS, the overall breeding objective achieved additional gain in GS.

RESUMEN

Maintaining genetic diversity and variation in livestock populations is critical for natural and artificial selection promoting genetic improvement while avoiding problems due to inbreeding. In Laos, there are concerns that there has been a decline in genetic diversity and a rise in inbreeding among native goats in their village-based smallholder system. In this study, we investigated the genetic diversity of Lao native goats in Phin, Songkhone and Sepon districts in Central Laos for the first time using Illumina's Goat SNP50 BeadChip. We also explored the genetic relationships between Lao goats with 163 global goat populations from 36 countries. Our results revealled a close genetic relationship between Lao native goats and Chinese, Mongolian and Pakistani goats, sharing ancestries with Guangfen, Jining Grey and Luoping Yellow breeds (China) and Teddi goats (Pakistan). The observed (Ho) and expected (He) heterozygosity were 0.292 and 0.303 (Laos), 0.288 and 0.288 (Sepon), 0.299 and 0.308 (Phin) and 0.289 and 0.305 (Songkhone), respectively. There was low to moderate genetic differentiation (FST: 0.011-0.043) and negligible inbreeding coefficients (FIS: -0.001 to 0.052) between goat districts. The runs of homozygosity (ROH) had an average length of 5.92-6.85 Mb, with short ROH segments (1-5 Mb length) being the most prevalent (66.34%). Longer ROH segments (20-40 and >40 Mb length categories) were less common, comprising only 4.81% and 1.01%, respectively. Lao goats exhibit moderate genetic diversity, low-inbreeding levels and adequate effective population size. Some genetic distinctions between Lao goats may be explained by geographic and cultural features.

Asunto(s)

RESUMEN

Random regression (RR) models are recommended as an alternative to multiple-trait (MT) models for better capturing the variance-covariance structure over a trajectory and hence more accurate genetic evaluation of traits that are repeatedly measured and genetically change gradually over time. However, a limited number of studies have been done to empirically compare RR over a MT model to determine how much extra benefit could be achieved from one method over another. We compared the prediction accuracy of RR and MT models for growth traits of Australian meat sheep measured from 60 to 525 d, using 102,579 weight records from 24,872 animals. Variance components and estimated breeding values (EBVs) estimated at specific ages were compared and validated with forward prediction. The accuracy of EBVs obtained from the MT model was 0.58, 0.51, 0.54, and 0.56 for weaning, postweaning, yearling, and hogget weight stages, respectively. RR model produced accuracy estimates of 0.56, 0.51, 0.54, and 0.54 for equivalent weight stages. Regression of adjusted phenotype on EBVs was very similar between the MT and the RR models (P > 0.05). Although the RR model did not significantly increase the accuracy of predicting future progeny performance, there are other benefits of the model such as no limit to the number of records per animal, estimation of EBVs for early and late growth, no need for age correction. Therefore, RR can be considered a more flexible method for the genetic evaluation of Australian sheep for early and late growth, and no need for age correction.

Currently, multiple-trait (MT) models are used in large-scale genetic evaluation of growth traits, where body weight traits are defined as separate traits at a finite number of fixed ages. Random regression (RR) models are expected to be superior since they can handle repeated measurements of weight and model these as a function of the actual age of measurement. These two models were compared in predicting breeding values for the body weight of Australian meat sheep. Phenotypic variation and estimated breeding values (EBVs) estimated at specific ages between 60 and 525 d with RR and MT models were compared and EBVs were validated in progeny data. The accuracy of EBVs in forecasting the performance of progeny was not statistically different between the two models. Other benefits of the RR model include the use of multiple records per animal, estimation of EBVs for early and late growth, with no need for age correction. Hence, RR models can be useful for the genetic evaluation of growth traits of sheep in Australia, but they do not necessarily predict breeding values at different ages more accurately than MT models.

Asunto(s)

RESUMEN

BACKGROUND: There can be variation between animals in how stable their genetic merit is across different environments due to genotype-by-environment (G×E) interactions. This variation could be used in breeding programs to select robust genotypes that combine high overall performance with stable genetic ranking across environments. There have been few attempts to validate breeding values for robustness in livestock, although this is a necessary step towards their implementation in selection decisions. The objective of this study was to validate breeding values for the robustness of body weight across different growth environments that were estimated using reaction norm models in sheep data. RESULTS: Using threefold cross-validation for the progeny of 337 sires, the average correlation between single-step breeding values for the reaction norm slope and the realised robustness of progeny across different growth environments was 0.21. The correlation between breeding values for the reaction slope estimated independently in two different datasets linked by common sires was close to the expected correlation based on theory. CONCLUSIONS: Slope estimated breeding values (EBV) obtained using reaction norm models were predictive of the phenotypic robustness of progeny across different environments and were consistent for sires with progeny in two different datasets. Selection based on reaction norm EBV could be used to increase the robustness of a population to environmental variation.

Asunto(s)

RESUMEN

BACKGROUND: Commercial poultry production systems follow a pyramidal structure with a nucleus of purebred animals under controlled conditions at the top and crossbred animals under commercial production conditions at the bottom. Genetic correlations between the same phenotypes on nucleus and production animals can therefore be influenced by differences both in purebred-crossbred genotypes and in genotype-by-environment interactions across the two environments, known as macro-genetic environmental sensitivity (GES). Within each environment, genotype-by-environment interactions can also occur due to so-called micro-GES. Micro-GES causes heritable variation in phenotypes and decreases uniformity. In this study, genetic variances of body weight (BW) and of micro-GES of BW and the impacts of purebred-crossbred differences and macro-environmental differences on micro-GES of BW were estimated. The dataset contained three subpopulations of slow-growing broiler chickens: purebred chickens (PB) reared in France, and crossbred chickens reared in France (FR) under the same conditions as PB or reared in Burkina Faso (BF) under local conditions. The crossbred chickens were offspring of the same dam line and had PB as their sire line. RESULTS: Estimates of heritability of BW and micro-GES of BW were 0.54 (SE of 0.02) and 0.06 (0.01), 0.67 (0.03) and 0.03 (0.01), and 0.68 (0.04) and 0.02 (0.01) for the BF, FR, and PB subpopulations, respectively. Estimates of the genetic correlations for BW between the three subpopulations were moderately positive (0.37 to 0.53) and those for micro-GES were weakly to moderately positive (0.01 to 0.44). CONCLUSIONS: The results show that the heritability of the micro-GES of BW varies with macro-environment, which indicates that responses to selection are expected to differ between macro-environments. The weak to moderate positive genetic correlations between subpopulations indicate that both macro-environmental differences and purebred-crossbred differences can cause re-ranking of sires based on their estimated breeding values for micro-GES of BW. Thus, the sire that produces the most variable progeny in one macro-environment may not be the one that produces the most variable offspring in another. Similarly, the sire that produces the most variable purebred progeny may not produce the most variable crossbred progeny. The results highlight the need for investigating micro-GES for all subpopulations included in the selection scheme, to ensure optimal genetic gain in all subpopulations.

Asunto(s)

RESUMEN

Microbial communities inhabiting the gut have the ability to influence physiological processes contributing to livestock production and performance. Livestock enterprises rely on animal production traits such as growth performance for profit. Previous studies have shown that gut microbiota are correlated to growth performance and could even influence it. The aim of this study was to characterise the faecal microbial profiles of Angus steers with high and low ADG at both weaning and yearling stages by profiling 16S rRNA gene sequences from rectal faecal samples. When microbial profiles were compared in terms of relative abundances, LEfSe analysis, alpha diversity metrics, and beta diversity, at the weaning stage, few significant differences were found between the high and low ADG groups. However, at yearling stage, microbial profiles significantly differed between the high and low ADG groups. The relative abundances of eight phyla and six genera significantly differed between the two groups. Alpha diversity metrics showed a significant decrease (p = 0.001) in species richness in the high ADG group. Similarly, beta diversity analysis showed that samples clustered clearly according to high and low ADG groups at yearling stage, indicating that phylogenetic similarity between the two ADG groups was significantly reduced (p = 0.005).

RESUMEN

The H-matrix best linear unbiased prediction (HBLUP) method has been widely used in livestock breeding programs. It can integrate all information, including pedigree, genotypes, and phenotypes on both genotyped and non-genotyped individuals into one single evaluation that can provide reliable predictions of breeding values. The existing HBLUP method requires hyper-parameters that should be adequately optimised as otherwise the genomic prediction accuracy may decrease. In this study, we assess the performance of HBLUP using various hyper-parameters such as blending, tuning, and scale factor in simulated and real data on Hanwoo cattle. In both simulated and cattle data, we show that blending is not necessary, indicating that the prediction accuracy decreases when using a blending hyper-parameter <1. The tuning process (adjusting genomic relationships accounting for base allele frequencies) improves prediction accuracy in the simulated data, confirming previous studies, although the improvement is not statistically significant in the Hanwoo cattle data. We also demonstrate that a scale factor, α, which determines the relationship between allele frequency and per-allele effect size, can improve the HBLUP accuracy in both simulated and real data. Our findings suggest that an optimal scale factor should be considered to increase prediction accuracy, in addition to blending and tuning processes, when using HBLUP.

RESUMEN

KEY MESSAGE: The reaction norm analysis of stability can be enhanced by partitioning the contribution of different types of G × E to the variation in slope. The slope of regression in a reaction norm model, where the performance of a genotype is regressed over an environmental covariable, is often used as a measure of stability of genotype performance. This method could be developed further by partitioning variation in the slope of regression into the two sources of genotype-by-environment interaction (G × E) which cause it: scale-type G × E (heterogeneity of variance) and rank-type G × E (heterogeneity of correlation). Because the two types of G × E have very different properties, separating their effect would enable a clearer understanding of stability. The aim of this paper was to demonstrate two methods which seek to achieve this in reaction norm models. Reaction norm models were fit to yield data from a multi-environment trial in Barley (Hordeum vulgare), with the adjusted mean yield from each environment used as the environmental covariable. Stability estimated from factor-analytic models, which can disentangle the two types of G × E and estimate stability based on rank-type G × E, was used for comparison. Adjusting the reaction norm slope to account for scale-type G × E using a genetic regression more than tripled the correlation with factor-analytic estimates of stability (0.24-0.26 to 0.80-0.85), indicating that it removed variation in the reaction norm slope that originated from scale-type G × E. A standardisation procedure had a more modest increase (055-0.59) but could be useful when curvilinear reaction norms are required. Analyses which use reaction norms to explore the stability of genotypes could gain additional insight into the mechanisms of stability by applying the methods outlined in this study.

Asunto(s)

RESUMEN

Most carcass and meat quality traits are moderate to highly heritable, indicating that they can be improved through selection. Genetic evaluation for these types of traits is performed using performance data obtained from commercial and progeny testing evaluation. The performance data from commercial farms are available in large volume, however, some drawbacks have been observed. The drawback of the commercial data is mainly due to sorting of animals based on live weight prior to slaughter, and this could lead to bias in the genetic evaluation of later measured traits such as carcass traits. The current study has two components to address the drawback of the commercial data. The first component of the study aimed to estimate genetic parameters for carcass and meat quality traits in Korean Hanwoo cattle using a large sample size of industry-based carcass performance records (n = 469,002). The second component of the study aimed to describe the impact of sorting animals into different contemporary groups based on an early measured trait and then examine the effect on the genetic evaluation of subsequently measured traits. To demonstrate our objectives, we used real performance data to estimate genetic parameters and simulated data was used to assess the bias in genetic evaluation. The results of our first study showed that commercial data obtained from slaughterhouses is a potential source of carcass performance data and useful for genetic evaluation of carcass traits to improve beef cattle performance. However, we observed some harvesting effect which leads to bias in genetic evaluation of carcass traits. This is mainly due to the selection of animal based on their body weight before arrival to slaughterhouse. Overall, the non-random allocation of animals into a contemporary group leads to a biased estimated breeding value in genetic evaluation, the severity of which increases when the evaluation traits are highly correlated.

RESUMEN

The aim of this study was to find significant genomic regions associated with carcass traits in Hanwoo cattle and to compare the benefit of using additional information from non-genotyped animals. Imputed whole-genome sequence data were used along with phenotypic data on 13 715 genotyped animals as well as phenotypes of 440 284 non-genotyped animals that were offspring of 454 genotyped sires. For carcass weight, 15 083 SNPs in 33 QTL regions and 313 candidate genes were identified. We found 410 SNPs in 17 QTL regions containing 122 candidate genes for back fat thickness. In total, 656 SNPs in 19 QTLs with 137 candidate genes for eye muscle area and 79 SNPs in 12 QTL regions with 77 candidate genes were identified for marbling score. The most important candidate genes included ZFAT, TG, PLAG1, CHCHD7, and TOX for carcass weight and eye muscle area, NOG for back fat thickness, and EVOVL5 for marbling score. This study showed that the use of phenotypic records on non-genotyped progeny along with imputed whole-genome sequence data increased the power of detecting new significant genomic regions.

Asunto(s)

RESUMEN

BACKGROUND: Selection of livestock based on their robustness or sensitivity to environmental variation could help improve the efficiency of production systems, particularly in the light of climate change. Genetic variation in robustness arises from genotype-by-environment (G × E) interactions, with genotypes performing differently when animals are raised in contrasted environments. Understanding the nature of this genetic variation is essential to implement strategies to improve robustness. In this study, our aim was to explore the genetics of robustness in Australian sheep to different growth environments using linear reaction norm models (RNM), with post-weaning weight records of 22,513 lambs and 60 k single nucleotide polymorphisms (SNPs). The use of scale-corrected genomic estimated breeding values (GEBV) for the slope to account for scale-type G × E interactions was also investigated. RESULTS: Additive genetic variance was observed for the slope of the RNM, with genetic correlations between low- and high-growth environments indicating substantial re-ranking of genotypes (0.44-0.49). The genetic variance increased from low- to high-growth environments. The heritability of post-weaning body weight ranged from 0.28 to 0.39. The genetic correlation between intercept and slope of the reaction norm for post-weaning body weight was low to moderate when based on the estimated (co)variance components but was much higher when based on back-solved SNP effects. An initial analysis suggested that a region on chromosome 11 affected both the intercept and the slope, but when the GEBV for the slope were conditioned on the GEBV for the intercept to remove the effect of scale-type G × E interactions on SNP effects for robustness, a single genomic region on chromosome 7 was found to be associated with robustness. This region included genes previously associated with growth traits and disease susceptibility in livestock. CONCLUSIONS: This study shows a significant genetic variation in the slope of RNM that could be used for selecting for increased robustness of sheep. Both scale-type and rank-type G × E interactions contributed to variation in the slope. The correction for scale effects of GEBV for the slope should be considered when analysing robustness using RNM. Overall, robustness appears to be a highly polygenic trait.

Asunto(s)

RESUMEN

BACKGROUND: South Africa and Australia shares multiple important sheep breeds. For some of these breeds, genomic breeding values are provided to breeders in Australia, but not yet in South Africa. Combining genomic resources could facilitate development for across country selection, but the influence of population structures could be important to the compatability of genomic data from varying origins. The genetic structure within and across breeds, countries and strains was evaluated in this study by population genomic parameters derived from SNP-marker data. Populations were first analysed by breed and country of origin and then by subpopulations of South African and Australian Merinos. RESULTS: Mean estimated relatedness according to the genomic relationship matrix varied by breed (-0.11 to 0.16) and bloodline (-0.08 to 0.06) groups and depended on co-ancestry as well as recent genetic links. Measures of divergence across bloodlines (FST: 0.04-0.12) were sometimes more distant than across some breeds (FST: 0.05-0.24), but the divergence of common breeds from their across-country equivalents was weak (FST: 0.01-0.04). According to mean relatedness, FST, PCA and Admixture, the Australian Ultrafine line was better connected to the SA Cradock Fine Wool flock than with other AUS bloodlines. Levels of linkage disequilibrium (LD) between adjacent markers was generally low, but also varied across breeds (r2: 0.14-0.22) as well as bloodlines (r2: 0.15-0.19). Patterns of LD decay was also unique to breeds, but bloodlines differed only at the absolute level. Estimates of effective population size (Ne) showed genetic diversity to be high for the majority of breeds (Ne: 128-418) but also for bloodlines (Ne: 137-369). CONCLUSIONS: This study reinforced the genetic complexity and diversity of important sheep breeds, especially the Merino breed. The results also showed that implications of isolation can be highly variable and extended beyond breed structures. However, knowledge of useful links across these population substructures allows for a fine-tuned approach in the combination of genomic resources. Isolation across country rarely proved restricting compared to other structures considered. Consequently, research into the accuracy of across-country genomic prediction is recommended.

Asunto(s)

RESUMEN

Economic values for annual milk yield (MY, kg), annual fat yield (FY, kg), annual protein yield (PY, kg), age at first calving (AFC, days), number of services per conception (NSC), calving interval (CI, days) and mastitis episodes (MS) were derived for temperate dairy cattle breeds in tropical Sri Lanka using a bio-economic model. Economic values were calculated on a per cow per year basis. Derived economic values in rupees (LKR) for MY, FY and PY were 107, -162 and -15, while for AFC, NSC, CI and MS, economic values were -59, -270, -84 and -8,303. Economic values for FY and PY further decreased with higher feed prices, and a less negative economic value for FY was obtained with increased price for fat. Negative economic values for FY and PY show that genetic improvement for these traits is not economical due to the high feed costs and/or the insufficient payment for fat and protein. Therefore, revision of milk fat and protein payments is recommended. Furthermore, the breeding objective developed in this study was dominated by milk production and fertility traits. Adaptability and functional traits that are important in a temperate dairy cattle breeding programme in tropical Sri Lanka, such as longevity, feed efficiency, disease resistance and heat tolerance should be recorded to incorporate them in the breeding objective. Continued trait recording of all traits is recommended to ensure dairy cows can be selected more effectively in a tropical environment based on a breeding objective that also includes adaptability and functional traits.

Asunto(s)

RESUMEN

The objective of this study was to investigate the accuracy of genomic prediction of body weight and eating quality traits in a numerically small sheep population (Dorper sheep). Prediction was based on a large multi-breed/admixed reference population and using (a) 50k or 500k single nucleotide polymorphism (SNP) genotypes, (b) imputed whole-genome sequencing data (~31 million), (c) selected SNPs from whole genome sequence data and (d) 50k SNP genotypes plus selected SNPs from whole-genome sequence data. Furthermore, the impact of using a breed-adjusted genomic relationship matrix on accuracy of genomic breeding value was assessed. The selection of genetic variants was based on an association study performed on imputed whole-genome sequence data in an independent population, which was chosen either randomly from the base population or according to higher genetic proximity to the target population. Genomic prediction was based on genomic best linear unbiased prediction (GBLUP), and the accuracy of genomic prediction was assessed according to the correlation between genomic breeding value and corrected phenotypes divided by the square root of trait heritability. The accuracy of genomic prediction was between 0.20 and 0.30 across different traits based on common 50k SNP genotypes, which improved on average by 0.06 (absolute value) on average based on using prioritized genetic markers from whole-genome sequence data. Using prioritized genetic markers from a genetically more related GWAS population resulted in slightly higher prediction accuracy (0.02 absolute value) compared to genetic markers derived from a random GWAS population. Using high-density SNP genotypes or imputed whole-genome sequence data in GBLUP showed almost no improvement in genomic prediction accuracy however, accounting for different marker allele frequencies in reference population according to a breed-adjusted GRM resulted to on average 0.024 (absolute value) increase in accuracy of genomic prediction.

Asunto(s)

RESUMEN

The objective of this study was to compare the accuracies of genomic prediction for milk yield, fat yield, and protein yield from Philippine dairy buffaloes using genomic best linear unbiased prediction (GBLUP) and single-step GBLUP (ssGBLUP) with the accuracies based on pedigree BLUP (pBLUP). To also assess the bias of the prediction, the regression coefficient (slope) of the adjusted phenotypes on the predicted breeding values (BVs) was also calculated. Two data sets were analyzed. The GENO data consisting of all female buffaloes that have both phenotypes and genotypes (n = 904 with 1,773,305-days lactation records) were analyzed using pBLUP and GBLUP. The ALL data, consisting of the GENO data plus females with phenotypes but not genotyped (n = 1,975 with 3,821,305-days lactation records), were analyzed using pBLUP and ssGBLUP. Animals were genotyped with the Affymetrix 90k buffalo genotyping array. After quality control, 60,827 single-nucleotide polymorphisms were used for downward analysis. A pedigree file containing 2,642 animals was used for pBLUP and ssGBLUP. Accuracy of prediction was calculated as the correlation between the predicted BVs of the test set and adjusted phenotypes, which were corrected for fixed effects, divided by the square root of the heritability of the trait, corrected for the number of lactations used in the test set. To assess the bias of the prediction, the regression coefficient (slope) of the adjusted phenotypes on the predicted BVs was also calculated. Results showed that genomic methods (GBLUP and ssGBLUP) provide more accurate predictions compared to pBLUP. Average GBLUP and ssGBLUP accuracies were 0.24 and 0.29, respectively, whereas average pBLUP accuracies (for GENO and ALL data) were 0.21 and 0.22, respectively. Slopes of the two genomic methods were also closer to one, indicating lesser bias, compared to pBLUP. Average GBLUP and ssGBLUP slopes were 0.89 and 0.84, respectively, whereas the average pBLUP (for GENO and ALL data) slopes were 0.80 and 0.54, respectively.

RESUMEN

BACKGROUND: Imputation to whole-genome sequence is now possible in large sheep populations. It is therefore of interest to use this data in genome-wide association studies (GWAS) to investigate putative causal variants and genes that underpin economically important traits. Merino wool is globally sought after for luxury fabrics, but some key wool quality attributes are unfavourably correlated with the characteristic skin wrinkle of Merinos. In turn, skin wrinkle is strongly linked to susceptibility to "fly strike" (Cutaneous myiasis), which is a major welfare issue. Here, we use whole-genome sequence data in a multi-trait GWAS to identify pleiotropic putative causal variants and genes associated with changes in key wool traits and skin wrinkle. RESULTS: A stepwise conditional multi-trait GWAS (CM-GWAS) identified putative causal variants and related genes from 178 independent quantitative trait loci (QTL) of 16 wool and skin wrinkle traits, measured on up to 7218 Merino sheep with 31 million imputed whole-genome sequence (WGS) genotypes. Novel candidate gene findings included the MAT1A gene that encodes an enzyme involved in the sulphur metabolism pathway critical to production of wool proteins, and the ESRP1 gene. We also discovered a significant wrinkle variant upstream of the HAS2 gene, which in dogs is associated with the exaggerated skin folds in the Shar-Pei breed. CONCLUSIONS: The wool and skin wrinkle traits studied here appear to be highly polygenic with many putative candidate variants showing considerable pleiotropy. Our CM-GWAS identified many highly plausible candidate genes for wool traits as well as breech wrinkle and breech area wool cover.

Asunto(s)