Investigation of potential genetic factors for growth traits in yellow-feather broilers using weighted single-step genome-wide association study.

Yellow-feather broilers take a large portion of poultry industry in China due to its meat characteristics. Improving the growth traits of yellow-feathered broilers will have great significance for the Chinese poultry market. The current study was designed to investigate the potential genetic factors using the weighted single-step genome-wide association study (wssGWAS) method, which takes consideration of more factors including pedigree, sex, environment and has more accuracy than traditional GWAS. The yellow-feather dwarf chickens from Wens Nanfang Poultry Breeding Co. Ltd. were revolved to recode 9 growth traits: Average daily gain (ADG), body weight (BW) at 45 d, 49 d, 56 d, 63 d, 70 d, 77 d, 84 d, 91 d for analysis. For the results, the region 4.63 to 5.03 Mb on chromosome 15, which was the QTL overlapped in BW45, BW49, BW56, BW63, BW84, might be the crucial genetic region for growth traits. Seven GO terms and 3 KEGG pathways, GO:0005200, GO:0005882, GO:0045111, GO:0099513, GO:0099081, GO:0099512, GO:0099080, KEGG:gga04020, KEGG:gga04540, KEGG:gga04210, were detected to relevant with growth traits. The genes enriched in these biological processes (NRAS, TUBB1, ADORA2B, NTRK3, NGF, TNNC2, F-KER, LOC429492, LOC431325, LOC431324, LOC396480) might have the function in growth of yellow-feather broilers.

Haplotype analysis of genomic prediction by incorporating genomic pathway information based on high-density SNP marker in Chinese yellow-feathered chicken.

Genomic selection using single nucleotide polymorphism (SNP) markers is now intensively investigated in breeding and has been widely utilized for genetic improvement. Currently, several studies have used haplotype (consisting of multiallelic SNPs) for genomic prediction and revealed its performance advantage. In this study, we comprehensively evaluated the performance of haplotype models for genomic prediction in 15 traits, including 6 growth, 5 carcass, and 4 feeding traits in a Chinese yellow-feathered chicken population. We adopted 3 methods to define haplotypes from high-density SNP panels, and our strategy included combining Kyoto Encyclopedia of Genes and Genomes pathway information and considering linkage disequilibrium (LD) information. Our results showed an increase in prediction accuracy due to haplotypes ranging from -0.04∼27.16% in all traits, where the significant improvements were found in 12 traits. The estimates of haplotype epistasis heritability were strongly correlated with the accuracy increase by haplotype models. In addition, incorporating genomic annotation information could further increase the accuracy of the haplotype model, where the further increase in accuracy is significantly relative to the increase of relative haplotype epistasis heritability. The genomic prediction using LD information for constructing haplotypes has the best prediction performance among the 4 traits. These results uncovered that haplotype methods were beneficial for genomic prediction, and the accuracy could be further increased by incorporating genomic annotation information. Moreover, using LD information would potentially improve the performance of genomic prediction.

Genomic Prediction Using LD-Based Haplotypes in Combined Pig Populations.

The size of reference population is an important factor affecting genomic prediction. Thus, combining different populations in genomic prediction is an attractive way to improve prediction ability. However, combining multireference population roughly cannot increase the prediction accuracy as well as expected in pig. This may be due to different linkage disequilibrium (LD) pattern differences between population. In this study, we used the imputed whole-genome sequencing (WGS) data to construct LD-based haplotypes for genomic prediction in combined population to explore the impact of different single-nucleotide polymorphism (SNP) densities, variant representation (SNPs or haplotype alleles), and reference population size on the prediction accuracy for reproduction traits. Our results showed that genomic best linear unbiased prediction (GBLUP) using the WGS data can improve prediction accuracy in multi-population but not within-population. Not only the genomic prediction accuracy of the haplotype method using 80 K chip data in multi-population but also GBLUP for the multi-population (3.4-5.9%) was higher than that within-population (1.2-4.3%). More importantly, we have found that using the haplotype method based on the WGS data in multi-population has better genomic prediction performance, and our results showed that building haploblock in this scenario based on low LD threshold (r 2 = 0.2-0.3) produced an optimal set of variables for reproduction traits in Yorkshire pig population. Our results suggested that whether the use of the haplotype method based on the chip data or GBLUP (individual SNP method) based on the WGS data were beneficial for genomic prediction in multi-population, while simultaneously combining the haplotype method and WGS data was a better strategy for multi-population genomic evaluation.

Impact of linkage disequilibrium heterogeneity along the genome on genomic prediction and heritability estimation.

BACKGROUND: Compared to medium-density single nucleotide polymorphism (SNP) data, high-density SNP data contain abundant genetic variants and provide more information for the genetic evaluation of livestock, but it has been shown that they do not confer any advantage for genomic prediction and heritability estimation. One possible reason is the uneven distribution of the linkage disequilibrium (LD) along the genome, i.e., LD heterogeneity among regions. The aim of this study was to effectively use genome-wide SNP data for genomic prediction and heritability estimation by using models that control LD heterogeneity among regions. METHODS: The LD-adjusted kinship (LDAK) and LD-stratified multicomponent (LDS) models were used to control LD heterogeneity among regions and were compared with the classical model that has no such control. Simulated and real traits of 2000 dairy cattle individuals with imputed high-density (770K) SNP data were used. Five types of phenotypes were simulated, which were controlled by very strongly, strongly, moderately, weakly and very weakly tagged causal variants, respectively. The performances of the models with high- and medium-density (50K) panels were compared to verify that the models that controlled LD heterogeneity among regions were more effective with high-density data. RESULTS: Compared to the medium-density panel, the use of the high-density panel did not improve and even decreased prediction accuracies and heritability estimates from the classical model for both simulated and real traits. Compared to the classical model, LDS effectively improved the accuracy of genomic predictions and unbiasedness of heritability estimates, regardless of the genetic architecture of the trait. LDAK applies only to traits that are mainly controlled by weakly tagged causal variants, but is still less effective than LDS for this type of trait. Compared with the classical model, LDS improved prediction accuracy by about 13% for simulated phenotypes and by 0.3 to ~ 10.7% for real traits with the high-density panel, and by ~ 1% for simulated phenotypes and by - 0.1 to ~ 6.9% for real traits with the medium-density panel. CONCLUSIONS: Grouping SNPs based on regional LD to construct the LD-stratified multicomponent model can effectively eliminate the adverse effects of LD heterogeneity among regions, and greatly improve the efficiency of high-density SNP data for genomic prediction and heritability estimation.