Estimation of variance and genomic prediction using genotypes imputed from low‐density marker subsets for carcass traits in Japanese black cattle

The influence of genotype imputation using low‐density single nucleotide polymorphism (SNP) marker subsets on the genomic relationship matrix (G matrix), genetic variance explained, and genomic prediction (GP) was investigated for carcass weight and marbling score in Japanese Black fattened steers, using genotype data of approximately 40,000 SNPs. Genotypes were imputed using equally spaced SNP subsets of different densities. Two different linear models were used. The first (model 1) incorporated one G matrix, while the second (model 2) used two different G matrices constructed using the selected and remaining SNPs. When using model 1, the estimated additive genetic variance was always larger when using all SNPs obtained via genotype imputation than when using only equally spaced SNP subsets. The correlations between the genomic estimated breeding values obtained using genotype imputation with at least 3,000 SNPs and those using all available SNPs without imputation were higher than 0.99 for both traits. While additive genetic variance was likely to be partitioned with model 2, it did not enhance the accuracy of GP compared with model 1. These results indicate that genotype imputation using an equally spaced low‐density panel of an appropriate size can be used to produce a cost‐effective, valid GP.     

Improving accuracy of genomic prediction in Brangus cattle by adding animals with imputed low‐density SNP genotypes

Reliable genomic prediction of breeding values for quantitative traits requires the availability of sufficient number of animals with genotypes and phenotypes in the training set. As of 31 October 2016, there were 3,797 Brangus animals with genotypes and phenotypes. These Brangus animals were genotyped using different commercial SNP chips. Of them, the largest group consisted of 1,535 animals genotyped by the GGP‐LDV4 SNP chip. The remaining 2,262 genotypes were imputed to the SNP content of the GGP‐LDV4 chip, so that the number of animals available for training the genomic prediction models was more than doubled. The present study showed that the pooling of animals with both original or imputed 40K SNP genotypes substantially increased genomic prediction accuracies on the ten traits. By supplementing imputed genotypes, the relative gains in genomic prediction accuracies on estimated breeding values (EBV) were from 12.60% to 31.27%, and the relative gain in genomic prediction accuracies on de‐regressed EBV was slightly small (i.e. 0.87%–18.75%). The present study also compared the performance of five genomic prediction models and two cross‐validation methods. The five genomic models predicted EBV and de‐regressed EBV of the ten traits similarly well. Of the two cross‐validation methods, leave‐one‐out cross‐validation maximized the number of animals at the stage of training for genomic prediction. Genomic prediction accuracy (GPA) on the ten quantitative traits was validated in 1,106 newly genotyped Brangus animals based on the SNP effects estimated in the previous set of 3,797 Brangus animals, and they were slightly lower than GPA in the original data. The present study was the first to leverage currently available genotype and phenotype resources in order to harness genomic prediction in Brangus beef cattle.     

Accuracy of imputation of single nucleotide polymorphism marker genotypes from low‐density panels in Japanese Black cattle

Using target and reference fattened steer populations, the performance of genotype imputation using lower‐density marker panels in Japanese Black cattle was evaluated. Population imputation was performed using BEAGLE software. Genotype information for approximately 40 000 single nucleotide polymorphism (SNP) markers by Illumina BovineSNP50 BeadChip was available, and imputation accuracy was assessed based on the average concordance rates of the genotypes, varying equally spaced SNP densities, and the number of individuals in the reference population. Two additional statistics were also calculated as indicators of imputation performance. The concordance rates tended to be lower for SNPs with greater minor allele frequencies, or those located near the ends of the chromosomes. Longer autosomes yielded greater imputation accuracies than shorter ones. When SNPs were selected based on linkage disequilibrium information, relative imputation accuracy was slightly improved. When 3000 and 10 000 equally spaced SNPs were used, the imputation accuracies were greater than 90% and approximately 97%, respectively. These results indicate that combining genotyping using a lower‐density SNP chip with genotype imputation based on a population of individuals genotyped using a higher‐density SNP chip is a cost‐effective and valid approach for genomic prediction.     

Identification of genome‐wide selection signatures in the Limousin beef cattle breed

The study is aimed at identifying selection footprints within the genome of Limousin cattle. With the use of Extended Haplotype Homozygosity test, supplemented with correction for variation in recombination rates across the genome, we created map of selection footprints and detected 173 significant (p < 0.01) core haplotypes being potentially under positive selection. Within these regions, a number of candidate genes associated inter alia with skeletal muscle growth (GDF15, BMP7, BMP4 and TGFB3) or postmortem proteolysis and meat maturation (CAPN1 and CAPN5) were annotated. Noticeable clusters of selection footprints were detected on chromosomes 1, 4, 8 and 14, which are known to carry several quantitative trait loci for growth traits and meat quality. The study provides information about the genes and metabolic pathways potentially modified under the influence of directional selection, aimed at improving beef production characteristics in Limousin cattle.     

Across breed multi‐trait random regression genomic predictions in the Nordic Red dairy cattle

The current study evaluates reliability of genomic predictions in selection candidates using multi‐trait random regression model, which accounts for interactions between marker effects and breed of origin in the Nordic Red dairy cattle (RDC). The population structure of the RDC is admixed. Data consisted of individual animal breed proportions calculated from the full pedigree, deregressed proofs (DRP) of published estimated breeding values (EBV) for yield traits and genotypic data for 37 595 single nucleotide polymorphic markers. The analysed data included 3330 bulls in the reference population and 812 bulls that were used for validation. Direct genomic breeding values (DGV) were estimated using the model under study, which accounts for breed effects and also with GBLUP, which assume uniform population. Validation reliability was calculated as a coefficient of determination from weighted regression of DRP on DGV (), scaled by the mean reliability of DRP. Using the breed‐specific model increased the reliability of DGV by 2 and 3% for milk and protein, respectively, when compared to homogeneous population GBLUP. The exception was for fat, where there was no gain in reliability. Estimated validation reliabilities were low for milk (0.32) and protein (0.32) and slightly higher (0.42) for fat.     

Genomic evaluation of regional dairy cattle breeds in single‐breed and multibreed contexts

An important prerequisite for high prediction accuracy in genomic prediction is the availability of a large training population, which allows accurate marker effect estimation. This requirement is not fulfilled in case of regional breeds with a limited number of breeding animals. We assessed the efficiency of the current French routine genomic evaluation procedure in four regional breeds (Abondance, Tarentaise, French Simmental and Vosgienne) as well as the potential benefits when the training populations consisting of males and females of these breeds are merged to form a multibreed training population. Genomic evaluation was 5–11% more accurate than a pedigree‐based BLUP in three of the four breeds, while the numerically smallest breed showed a < 1% increase in accuracy. Multibreed genomic evaluation was beneficial for two breeds (Abondance and French Simmental) with maximum gains of 5 and 8% in correlation coefficients between yield deviations and genomic estimated breeding values, when compared to the single‐breed genomic evaluation results. Inflation of genomic evaluation of young candidates was also reduced. Our results indicate that genomic selection can be effective in regional breeds as well. Here, we provide empirical evidence proving that genetic distance between breeds is only one of the factors affecting the efficiency of multibreed genomic evaluation.     

Effects of a breeding scheme combined by genomic pre‐selection and progeny testing on annual genetic gain in a dairy cattle population

The effectiveness of the incorporation of genomic pre‐selection into dairy cattle progeny testing (GS‐PT) was compared with that of progeny testing (PT) where the fraction of dam to breed bull (DB) selected was 0.01. When the fraction of sires to breed bulls (SB) selected without being progeny tested to produce young bulls (YB) in the next generation was 0.2, the annual genetic gain from GS‐PT was 13% to 43% greater when h² = 0.3 and 16% to 53% greater when h² = 0.1 compared with that from PT. Given h² = 0.3, a selection accuracy of 0.8 for both YB and DB, and selected fractions of 0.117 for YB and 0.04 for DB, GS‐PT produced 40% to 43% greater annual genetic gain than PT. Given h² = 0.1, a selection accuracy of 0.6 for both YB and DB, and selected fractions of 0.117 for YB and 0.04 for DB, annual genetic gain from GS‐PT was 48% to 53% greater than that from PT. When h² = 0.3, progeny testing capacity had little effect on annual genetic gain from GS‐PT. However, when h² = 0.1, annual genetic gain from GS‐PT increased with increasing progeny testing capacity.     

Use of SNP genotyping to determine pedigree and breed composition of dairy cattle in Kenya

High levels of inbreeding in East African dairy cattle are a potential concern because of use of a limited range of imported germplasm coupled with strong selection, especially by disease, and sparse performance recording. To address this, genetic relationships and breed composition in an admixed population of Kenyan dairy cattle were estimated by means of a 50K SNP scan. Genomic DNA from 3 worldwide Holstein and 20 Kenyan bulls, 71 putative cow‐calf pairs, 25 cows from a large ranch and 5 other Kenyan animals were genotyped for 37 238 informative SNPs. Sires were predicted and 89% of putative dam‐calf relationships were supported by genotype data. Animals were clustered with the HapMap population using Structure software to assess breed composition. Cows from a large ranch primarily clustered with Holsteins, while animals from smaller farms were generally crosses between Holstein and Guernsey. Coefficients of relatedness were estimated and showed evidence of heavy use of one AI bull. We conclude that little native germplasm exists within the genotyped populations and mostly European ancestry remains.     

Application of imputation methods to genomic selection in Chinese Holstein cattle

Missing genotypes are a common feature of high density SNP datasets obtained using SNP chip technology and this is likely to decrease the accuracy of genomic selection. This problem can be circumvented by imputing the missing genotypes with estimated genotypes. When implementing imputation, the criteria used for SNP data quality control and whether to perform imputation before or after data quality control need to consider. In this paper, we compared six strategies of imputation and quality control using different imputation methods, different quality control criteria and by changing the order of imputation and quality control, against a real dataset of milk production traits in Chinese Holstein cattle. The results demonstrated that, no matter what imputation method and quality control criteria were used, strategies with imputation before quality control performed better than strategies with imputation after quality control in terms of accuracy of genomic selection. The different imputation methods and quality control criteria did not significantly influence the accuracy of genomic selection. We concluded that performing imputation before quality control could increase the accuracy of genomic selection, especially when the rate of missing genotypes is high and the reference population is small.     

Performance of genomic selection under a single-step approach in autochthonous Spanish beef cattle populations

This study evaluated different strategies for implementing a single-step genomic selection programme in two autochthonous Spanish beef cattle populations (Pirenaica—Pi and Rubia Gallega—RG). The strategies were compared in terms of accuracy attained under different scenarios by simulating genomic data over the known genealogy. Several genotyping approaches were tested, as well as, other factors like marker density, effective population size, mutation rate and heritability of the trait. The results obtained showed gains in accuracy with respect to pedigree BLUP evaluation in all cases. The greatest benefit was obtained when the candidates to selection had their genotypes included in the evaluation. Moreover, genotyping the individuals with the most accurate predictions maximized the gains but other suboptimal strategies also yielded satisfactory results. Furthermore, the gains in accuracy increased with the marker density reaching a plateau at around 50,000 markers. Likewise, the effective population size and the mutation rate have also shown an effect, both increasing the accuracy with decreasing values of these population parameters. Finally, the results obtained for the RG population showed greater gains compared to the Pi population, probably attributed to the wider implantation of artificial insemination.     

Development of genomic predictions for Angus cattle in Brazil incorporating genotypes from related American sires

Genomic prediction has become the new standard for genetic improvement programs, and currently, there is a desire to implement this technology for the evaluation of Angus cattle in Brazil. Thus, the main objective of this study was to assess the feasibility of evaluating young Brazilian Angus (BA) bulls and heifers for 12 routinely recorded traits using single-step genomic BLUP (ssGBLUP) with and without genotypes from American Angus (AA) sires. The second objective was to obtain estimates of effective population size (N_e) and linkage disequilibrium (LD) in the Brazilian Angus population. The dataset contained phenotypic information for up to 277,661 animals belonging to the Promebo breeding program, pedigree for 362,900, of which 1,386 were genotyped for 50k, 77k, and 150k single nucleotide polymorphism (SNP) panels. After imputation and quality control, 61,666 SNPs were available for the analyses. In addition, genotypes from 332 American Angus (AA) sires widely used in Brazil were retrieved from the AA Association database to be used for genomic predictions. Bivariate animal models were used to estimate variance components, traditional EBV, and genomic EBV (GEBV). Validation was carried out with the linear regression method (LR) using young-genotyped animals born between 2013 and 2015 without phenotypes in the reduced dataset and with records in the complete dataset. Validation animals were further split into progeny of BA and AA sires to evaluate if their progenies would benefit by including genotypes from AA sires. The N_e was 254 based on pedigree and 197 based on LD, and the average LD (±SD) and distance between adjacent single nucleotide polymorphisms (SNPs) across all chromosomes were 0.27 (±0.27) and 40743.68 bp, respectively. Prediction accuracies with ssGBLUP outperformed BLUP for all traits, improving accuracies by, on average, 16% for BA young bulls and heifers. The GEBV prediction accuracies ranged from 0.37 (total maternal for weaning weight and tick count) to 0.54 (yearling precocity) across all traits, and dispersion (LR coefficients) fluctuated between 0.92 and 1.06. Inclusion of genotyped sires from the AA improved GEBV accuracies by 2%, on average, compared to using only the BA reference population. Our study indicated that genomic information could help us to improve GEBV accuracies and hence genetic progress in the Brazilian Angus population. The inclusion of genotypes from American Angus sires heavily used in Brazil just marginally increased the GEBV accuracies for selection candidates.     

Accurate genomic predictions for BCWD resistance in rainbow trout are achieved using low‐density SNP panels: Evidence that long‐range LD is a major contributing factor

Previously accurate genomic predictions for Bacterial cold water disease (BCWD) resistance in rainbow trout were obtained using a medium‐density single nucleotide polymorphism (SNP) array. Here, the impact of lower‐density SNP panels on the accuracy of genomic predictions was investigated in a commercial rainbow trout breeding population. Using progeny performance data, the accuracy of genomic breeding values (GEBV) using 35K, 10K, 3K, 1K, 500, 300 and 200 SNP panels as well as a panel with 70 quantitative trait loci (QTL)‐flanking SNP was compared. The GEBVs were estimated using the Bayesian method BayesB, single‐step GBLUP (ssGBLUP) and weighted ssGBLUP (wssGBLUP). The accuracy of GEBVs remained high despite the sharp reductions in SNP density, and even with 500 SNP accuracy was higher than the pedigree‐based prediction (0.50–0.56 versus 0.36). Furthermore, the prediction accuracy with the 70 QTL‐flanking SNP (0.65–0.72) was similar to the panel with 35K SNP (0.65–0.71). Genomewide linkage disequilibrium (LD) analysis revealed strong LD (r² ≥ 0.25) spanning on average over 1 Mb across the rainbow trout genome. This long‐range LD likely contributed to the accurate genomic predictions with the low‐density SNP panels. Population structure analysis supported the hypothesis that long‐range LD in this population may be caused by admixture. Results suggest that lower‐cost, low‐density SNP panels can be used for implementing genomic selection for BCWD resistance in rainbow trout breeding programs.     

Differences in rumen fermentation characteristics between low‐yield and high‐yield dairy cows in early lactation

Relationship between rumen fermentation parameters, blood biochemical profiles and milk production traits in different yielding dairy cows during early lactation was investigated. Twelve dairy cows were divided into two groups based on their milk yield, that is low‐yield (LY) and high‐yield (HY) groups. Rumen fluid and blood were collected at 3 weeks prepartum and 4, 8 and 12 weeks postpartum. Results showed that proportions of acetate, propionate to total short chain fatty acids and acetate : propionate ratio were changed (P < 0.05) in both groups during the peripartum period, whereas butyrate and acetate : butyrate ratio were only altered in the HY group. Blood cholesterol, beta‐hydroxybutyric acid (BHBA) and glutamic oxaloacetic transaminase in the HY group were higher (P < 0.01) than those in the LY group. Principal component analysis revealed that milk yield and milk compositions were differently clustered between groups. These parameters showed similar direction with dry matter intake in the HY group and adverse direction in the LY group. Linear regression analysis indicated that butyrate was positively correlated with BHBA (P < 0.05) in the HY group. This study suggests that cows in the HY group seem to accommodate appropriately to negative energy balance in early lactation through rumen fermentation.     

Comparative digestibility of low‐quality grass hay by two breeds of cattle differing in mature live weight

The digestive capacity of the Portuguese native breed of Barrosão cattle and the Holstein–Friesian breed was measured when fed meadow hay (72–74 g CP and 641–671 g NDF/kg dry matter), offered either alone or supplemented with soya bean meal (150 g/kg dry matter), at maintenance level. Four mature cows of each breed were used. Average initial live weight (LW) was 457 and 635 kg for the Barrosão and the Holstein–Friesian cows respectively. The organic matter digestibility (OMD) of the meadow hay was higher in Holstein–Friesian than in Barrosão cows (p < 0.08) whether supplemented or not. The neutral detergent fibre digestibility of the hay was also higher in Holstein–Frisian when the hay was fed alone (p < 0.08). The soya bean meal supplementation increased the OMD and the NDFD digestibility of the total diet (p < 0.05), but not the OMD digestibility of the hay (p > 0.05). The results obtained in this study suggest a higher ability to digest fibre in the large dairy breed than in the small native breed. This suggestion is also supported by previous findings with sheep breeds largely differing in mature live weight.     

Comparing deregression methods for genomic prediction of test‐day traits in dairy cattle

We aimed to investigate the performance of three deregression methods (VanRaden, VR; Wiggans, WG; and Garrick, GR) of cows’ and bulls’ breeding values to be used as pseudophenotypes in the genomic evaluation of test‐day dairy production traits. Three scenarios were considered within each deregression method: (i) including only animals with reliability of estimated breeding value (REL_EBV ) higher than the average of parent reliability (REL_PA ) in the training and validation populations; (ii) including only animals with REL_EBV higher than 0.50 in the training and REL_EBV higher than REL_PA in the validation population; and (iii) including only animals with REL_EBV higher than 0.50 in both training and validation populations. Individual random regression coefficients of lactation curves were predicted using the genomic best linear unbiased prediction (GBLUP), considering either unweighted or weighted residual variances based on effective records contributions. In summary, VR and WG deregression methods seemed more appropriate for genomic prediction of test‐day traits without need for weighting in the genomic analysis, unless large differences in REL_EBV between training population animals exist.     

Effect of feeding a high‐carbohydrate or a high‐fat diet on subsequent food intake and blood concentration of satiety‐related hormones in dogs

Although studies in rodents and humans have evidenced a weaker effect of fat in comparison to carbohydrates on the suppression of food intake, very few studies have been carried out in this field in dogs. This study investigates the effects of a high‐carbohydrate (HC ) and a high‐fat (HF ) diets on subsequent food intake and blood satiety‐related hormones in dogs. Diets differed mainly in their starch (442 vs. 271 g/kg dry matter) and fat (99.3 vs. 214 g/kg dry matter) contents. Twelve Beagle dogs received the experimental diets at maintenance energy requirements in two experimental periods, following a cross‐over arrangement. In week 7 of each period, blood concentrations of active ghrelin, glucagon‐like peptide (GLP ‐1), peptide YY , insulin, and glucose were determined before and at 30, 60, 120, 180, and 360 min post‐feeding. The following week, intake of a challenge food offered 180 min after the HC and HF diets was recorded over two days. In comparison to the dogs on the HC diet, those on the HF diet had a higher basal concentration of GLP ‐1 (p  = .010) and a higher total area under the curve over 180 min post‐prandial (tAUC _0–180) (p  = .031). Dogs on the HC diet showed a higher elevation of ghrelin at 180 min (p  = .033) and of insulin at 360 min (p  = .041), although ghrelin and insulin tAUC _0–180 did not differ between the two diets (p  ? .10). Diet had no effect on challenge food intake (p  ? .10), which correlated with the tAUC _0–180 of ghrelin (r = .514, p  = .010), insulin (r = ?.595, p  = .002), and glucose (r = ?.516, p  = .010). Feeding a diet high in carbohydrate or fat at these inclusion levels does not affect the feeding response at 180 min post‐prandial, suggesting a similar short‐term satiating capacity.     

The importance of haplotype length and heritability using genomic selection in dairy cattle

Reliabilities for genomic estimated breeding values (GEBV) were investigated by simulation for a typical dairy cattle breeding setting. Scenarios were simulated with different heritabilites ( h ²) and for different haplotype sizes, and seven generations with only genotypes were generated to investigate reliability of GEBV over time. A genome with 5000 single nucleotide polymorphisms (SNP) at distances of 0.1 cM and 50 quantitative trait loci (QTL) was simulated, and a Bayesian variable selection model was implemented to predict GEBV. Highest reliabilities were obtained for 10 SNP haplotypes. At optimal haplotype size, reliabilities in generation 1 without phenotypes ranged from 0.80 for h ² = 0.02 to 0.93 for h ² = 0.30, and in the seventh generation without phenotypes ranged from 0.69 for h ² = 0.02 to 0.86 for h ² = 0.30. Reliabilities of GEBV were found sufficiently high to implement dairy selection schemes without progeny testing in which case a data time-lag of two to three generations may be present. Reliabilities were also relatively high for low heritable traits, implying that genomic selection could be especially beneficial to improve the selection on, e.g. health and fertility.     

Metagenomic profiles of the rumen microbiota during the transition period in low‐yield and high‐yield dairy cows

We investigated potential relationships between rumen microbiota and milk production in dairy cows during the transition period. Twelve dairy cows were divided into a low‐yield (LY) or high‐yield (HY) group based on their milk yield. Rumen samples were taken from dairy cows at 3 weeks before parturition, and at 4, 8, and 12 weeks after parturition. 16S rDNA‐based metagenomic analysis showed that diversities of rumen microbiota in both groups were similar and the number of operational taxonomic units (OTUs) was lower in the postpartum than prepartum period in both groups. The abundance of Bacteroidetes and ratio of Bacteroidetes:Firmicutes was higher in the HY than the LY group. OTUs assigned to Prevotella bryantii, Fibrobacter succinogenes, Ruminococcus albus, Butyrivibrio fibrisolvens, and Succinivibrio sp. were abundant in the HY group. These OTUs were significantly related to the propionate molar proportion of rumen fluids in the HY group. OTUs assigned to Lachnospiraceae, Bifidobacterium sp. and Saccharofermentans were dominant in the LY group. Predictive functional profiling revealed that abundance of gene families involved in amino acid and vitamin metabolism was higher in the HY than the LY group. These results suggest that the community structure and fermentation products of rumen microbiota could be associated with milk production of dairy cows.     

Significance testing and genomic inflation factor using high‐density genotypes or whole‐genome sequence data

Significance testing for genome‐wide association study (GWAS) with increasing SNP density up to whole‐genome sequence data (WGS) is not straightforward, because of strong LD between SNP and population stratification. Therefore, the objective of this study was to investigate genomic control and different significance testing procedures using data from a commercial pig breeding scheme. A GWAS was performed in GCTA with data of 4,964 Large White pigs using medium density, high density or imputed whole‐genome sequence data, fitting a genomic relationship matrix based on a leave‐one–chromosome‐out approach to account for population structure. Subsequently, genomic inflation factors were assessed on whole‐genome level and the chromosome level. To establish a significance threshold, permutation testing, Bonferroni corrections using either the total number of SNPs or the number of independent chromosome fragments, and false discovery rates (FDR) using either the Benjamini–Hochberg procedure or the Benjamini and Yekutieli procedure were evaluated. We found that genomic inflation factors did not differ between different density genotypes but do differ between chromosomes. Also, the leave‐one‐chromosome‐out approach for GWAS or using the pedigree relationships did not account appropriately for population stratification and gave strong genomic inflation. Regarding different procedures for significance testing, when the aim is to find QTL regions that are associated with a trait of interest, we recommend applying the FDR following the Benjamini and Yekutieli approach to establish a significance threshold that is adjusted for multiple testing. When the aim is to pinpoint a specific mutation, the more conservative Bonferroni correction based on the total number of SNPs is more appropriate, till an appropriate method is established to adjust for the number of independent tests.     

Signatures of selection in Charolais beef cattle identified by genome‐wide analysis

Charolais cattle are one of the most important breeds for meat production worldwide; in México, its selection is mainly made by live weight traits. One strategy for mapping important genomic regions that might influence productive traits is the identification of signatures of selection. This type of genomic features contains loci with extended linkage disequilibrium (LD) and homozygosity patterns that are commonly associated with sites of quantitative trait locus (QTL). Therefore, the objective of this study was to identify the signatures of selection in Charolais cattle genotyped with the GeneSeek Genomic Profiler Bovine HD panel consisting of 77 K single nucleotide polymorphisms (SNPs). A total 61,311 SNPs and 819 samples were used for the analysis. Identification of signatures of selection was carried out using the integrated haplotype score (iHS) methodology implemented in the rehh R package. The top ten SNPs with the highest piHS values were located on BTA 4, 5, 6 and 14. By identifying markers in LD with top ten SNPs, the candidate regions defined were mapped to 52.8–59.3 Mb on BTA 4; 67.5–69.3 on BTA 5; 39.5–41.0 Mb on BTA 6; and 26.4–29.6 Mb on BTA 14. The comparison of these candidate regions with the bovine QTLdb effectively confirmed the association (p < 0.05) with QTL related to growth traits and other important productive traits. The genomic regions identified in this study indicated selection for growth traits on the Charolais population via the conservation of haplotypes on various chromosomes. These genomic regions and their associated genes could serve as the basis for haplotype association studies and for the identification of causal genes related to growth traits.