Evaluation of updated Feed Saved breeding values developed in Australian Holstein dairy cattle.

Although selection for increased milk production traits has led to a genetic increase in body weight (BW), the genetic gain in milk production has exceeded the gain in BW, so gross feed efficiency has improved. Nonetheless, greater gains may be possible by directly selecting for a measure of feed efficiency. Australia first introduced Feed Saved (FS) estimated breeding value (EBV) in 2015. Feed Saved combines residual feed intake (RFI) genomic EBV and maintenance requirements calculated from mature BW EBV. The FS EBV was designed to enable the selection of cows for reduced energy requirements with similar milk production. In this study, we used a reference population of 3,711 animals in a multivariate analysis including Australian heifers (AUSh), Australian cows (AUSc), and overseas cows (OVEc) to update the Australian EBV for lifetime RFI (i.e., a breeding value that incorporated RFI in growing and lactating cows) and to recalculate the FS EBV in Australian Holstein bulls (AUSb). The estimates of genomic heritabilities using univariate (only AUSc or AUSh) to trivariate (including the OVEc) analyses were similar. Genomic heritabilities for RFI were estimated as 0.18 for AUSc, 0.27 for OVEc, and 0.36 for AUSh. The genomic correlation for RFI between AUSc and AUSh was 0.47 and that between AUSc and OVEc was 0.94, but these estimates were associated with large standard errors (range: 0.18-0.28). The reliability of lifetime RFI (a component of FS) in the trivariate analysis (i.e., including OVEc) increased from 11% to 20% compared with the 2015 model and was greater, by 12%, than in a bivariate analysis in which the reference population included only AUSc and AUSh. By applying the prediction equation of the 2020 model, the average reliability of the FS EBV in 20,816 AUSb that were born between 2010 and 2020 improved from 33% to 43%. Previous selection strategies-that is, using the predecessor of the Balanced Performance Index (Australian Profit Ranking index) that did not include FS-have resulted in an unfavorable genetic trend in FS. However, this unfavorable trend has stabilized since 2015, when FS was included in the Balanced Performance Index, and is expected to move in a favorable direction with selection on Balanced Performance Index or the Health Weighted Index. Doubling the reference population, particularly by incorporating international data for feed efficiency, has improved the reliability of the FS EBV. This could lead to increased genetic gain for feed efficiency in the Australian industry.

Erratum to "Evaluation of updated Feed Saved breeding values developed in Australian Holstein dairy cattle" (JDS Commun. 3:114-119).

[This corrects the article DOI: 10.3168/jdsc.2021-0150.].

The use of milk mid-infrared spectroscopy to improve genomic prediction accuracy of serum biomarkers.

Breeding objectives in the dairy industry have shifted from being solely focused on production to including fertility, animal health, and environmental impact. Increased serum concentrations of candidate biomarkers of health and fertility, such as ß-hydroxybutyric acid (BHB), fatty acids, and urea are difficult and costly to measure, and thus limit the number of records. Accurate genomic prediction requires a large reference population. The inclusion of milk mid-infrared (MIR) spectroscopic predictions of biomarkers may increase genomic prediction accuracy of these traits. Our objectives were to (1) estimate the heritability of, and genetic correlations between, selected serum biomarkers and their respective MIR predictions, and (2) evaluate genomic prediction accuracies of either only measured serum traits, or serum traits plus MIR-predicted traits. The MIR-predicted traits were either fitted in a single trait model, assuming the measured trait and predicted trait were the same trait, or in a multitrait model, where measured and predicted trait were assumed to be correlated traits. We performed all analyses using relationship matrices constructed from pedigree (A matrix), genotypes (G matrix), or both pedigree and genotypes (H matrix). Our data set comprised up to 2,198 and 9,657 Holstein cows with records for serum biomarkers and MIR-predicted traits, respectively. Heritabilities of measured serum traits ranged from 0.04 to 0.07 for BHB, from 0.13 to 0.21 for fatty acids, and from 0.10 to 0.12 for urea. Heritabilities for MIR-predicted traits were not significantly different from those for the measured traits. Genetic correlations between measured traits and MIR-predicted traits were close to 1 for urea. For BHB and fatty acids, genetic correlations were lower and had large standard errors. The inclusion of MIR predicted urea substantially increased prediction accuracy for urea. For BHB, including MIR-predicted BHB reduced the genomic prediction accuracy, whereas for fatty acids, prediction accuracies were similar with either measured fatty acids, MIR-predicted fatty acids, or both. The high genetic correlation between urea and MIR-predicted urea, in combination with the increased prediction accuracy, demonstrated the potential of using MIR-predicted urea for genomic prediction of urea. For BHB and fatty acids, further studies with larger data sets are required to obtain more accurate estimates of genetic correlations.

Value of sharing cow reference population between countries on reliability of genomic prediction for milk yield traits.

Increasing the reliability of genomic prediction (GP) of economic traits in the pasture-based dairy production systems of New Zealand (NZ) and Australia (AU) is important to both countries. This study assessed if sharing cow phenotype and genotype data of NZ and AU improves the reliability of GP for NZ bulls. Data from approximately 32,000 NZ genotyped cows and their contemporaries were included in the May 2018 routine genetic evaluation of the Australian Dairy cattle in an attempt to provide consistent phenotypes for both countries. After the genetic evaluation, deregressed proofs of cows were calculated for milk yield traits. The April 2018 multiple across-country evaluation of Interbull was also used to calculate deregressed proofs for bulls on the NZ scale. Approximately 1,178 Jersey (Jer) and 6,422 Holstein (Hol) bulls had genotype and phenotype data. In addition to NZ cows, phenotype data of close to 60,000 genotyped Australian (AU) cows from the same genetic evaluation run as NZ cows were used. All AU and NZ females were genotyped using low-density SNP chips (<10K SNP) and were imputed first to 50K and then to ∼600K (referred to as high density; HD). We used up to 98,000 animals in the reference populations, both by expanding the NZ reference set (cow, bull, single breed to multi-breed set) and by adding AU cows. Reliabilities of GP were calculated for 508 Jer and 1,251 Hol bulls whose sires are not included in the reference set (RS) to ensure that real differences are not masked by close relationships. The GP was tested using 50K or high-density SNP chip using genomic BLUP in bivariate (considering country as a trait) or single trait models. The RS that gave the highest reliability for each breed were also tested using a hybrid GP method that combines expectation maximization with Bayes R. The addition of the AU cows to an NZ RS that included either NZ cows only, or cows and bulls, improved the reliability of GP for both NZ Hol and Jer validation bulls for all traits. Using single breed reference populations also increased reliability when NZ crossbred cows were added to reference populations that included only purebred NZ bulls and cows and AU cows. The full multi-breed RS (all NZ cows and bulls and AU cows) provided similar reliabilities in NZ Hol bulls, when compared with the single breed reference with crossbred NZ cows. For Jer validation bulls, the RS that included Jer cows and bulls and crossbred cows from NZ and Jer cows from AU was marginally better than the all-breed, all-country RS. In terms of reliability, the advantage of the HD SNP chip was small but captured more of the genomic variance than the 50K, particularly for Hol. The expectation maximization Bayes R GP method was slightly (up to 3 percentage points) better than genomic BLUP. We conclude that GP of milk production traits in NZ bulls improves by up to 7 percentage points in reliability by expanding the NZ reference population to include AU cows.

Vitamin D deficiency is associated with tuberculosis infection among household contacts in Ulaanbaatar, Mongolia.

BACKGROUND: Vitamin D deficiency (VDD) is a known risk factor for tuberculous infection. We investigated if VDD is a risk factor for tuberculous infection among the household contacts (HHCs) of patients with tuberculosis (TB) in Mongolia.MATERIALS and METHOD: All HHCs of TB patients diagnosed in Khan-Uul District, Mongolia, were enrolled. The serum level of 25-hydroxyvitamin D [25(OH)D] was detected and TB infection determined using QuantiFERON-TB Gold Plus (QFT-Plus). A tuberculin skin test (TST) reading >10 mm was considered to be positive. Epidemiological and bacteriological data were collected from routine surveillance of the National Tuberculosis Programme.RESULTS: Among study participants, 48.2% (135/285) were QFT-Plus-positive. Of QFT-positive HHCs, 77.0% (104/135) were TST-positive and the overall concordance of tests was low (κ 0.374, P < 0.001). A low serum level of 25(OH)D was an independent predictor for QFT-Plus positivity (P < 0.001). CD8+ T-cell stimulation measured by QFT-Plus had borderline association with the serum level of 25(OH)D (P = 0.089).CONCLUSION: We showed a high rate of TB infection among HHCs in Mongolia. QFT-Plus could decrease the number of people requiring TB preventive treatment, in addition to aiding detection of new TB infection. A low serum level of vitamin D was an independent predictor of TB infection, but not a predictor of stimulation of CD8+ T cells.

The use of mid-infrared spectra to map genes affecting milk composition.

The aim of this study was to investigate the feasibility of using mid-infrared (MIR) spectroscopy analysis of milk samples to increase the power and precision of genome-wide association studies (GWAS) for milk composition and to better distinguish linked quantitative trait loci (QTL). To achieve this goal, we analyzed phenotypic data of milk composition traits, related MIR spectra, and genotypic data comprising 626,777 SNP on 5,202 Holstein, Jersey, and crossbred cows. We performed a conventional GWAS on protein, lactose, fat, and fatty acid concentrations in milk, a GWAS on individual MIR wavenumbers, and a partial least squares regression (PLS), which is equivalent to a multi-trait GWAS, exploiting MIR data simultaneously to predict SNP genotypes. The PLS detected most of the QTL identified using single-trait GWAS, usually with a higher significance value, as well as previously undetected QTL for milk composition. Each QTL tends to have a different pattern of effects across the MIR spectrum and this explains the increased power. Because SNP tracking different QTL tend to have different patterns of effect, it was possible to distinguish closely linked QTL. Overall, the results of this study suggest that using MIR data through either GWAS or PLS analysis applied to genomic data can provide a powerful tool to distinguish milk composition QTL.

Genomic prediction of reproduction traits for Merino sheep.

Economically important reproduction traits in sheep, such as number of lambs weaned and litter size, are expressed only in females and later in life after most selection decisions are made, which makes them ideal candidates for genomic selection. Accurate genomic predictions would lead to greater genetic gain for these traits by enabling accurate selection of young rams with high genetic merit. The aim of this study was to design and evaluate the accuracy of a genomic prediction method for female reproduction in sheep using daughter trait deviations (DTD) for sires and ewe phenotypes (when individual ewes were genotyped) for three reproduction traits: number of lambs born (NLB), litter size (LSIZE) and number of lambs weaned. Genomic best linear unbiased prediction (GBLUP), BayesR and pedigree BLUP analyses of the three reproduction traits measured on 5340 sheep (4503 ewes and 837 sires) with real and imputed genotypes for 510 174 SNPs were performed. The prediction of breeding values using both sire and ewe trait records was validated in Merino sheep. Prediction accuracy was evaluated by across sire family and random cross-validations. Accuracies of genomic estimated breeding values (GEBVs) were assessed as the mean Pearson correlation adjusted by the accuracy of the input phenotypes. The addition of sire DTD into the prediction analysis resulted in higher accuracies compared with using only ewe records in genomic predictions or pedigree BLUP. Using GBLUP, the average accuracy based on the combined records (ewes and sire DTD) was 0.43 across traits, but the accuracies varied by trait and type of cross-validations. The accuracies of GEBVs from random cross-validations (range 0.17-0.61) were higher than were those from sire family cross-validations (range 0.00-0.51). The GEBV accuracies of 0.41-0.54 for NLB and LSIZE based on the combined records were amongst the highest in the study. Although BayesR was not significantly different from GBLUP in prediction accuracy, it identified several candidate genes which are known to be associated with NLB and LSIZE. The approach provides a way to make use of all data available in genomic prediction for traits that have limited recording.

Design of a low-density SNP chip for the main Australian sheep breeds and its effect on imputation and genomic prediction accuracy.

Genotyping sheep for genome-wide SNPs at lower density and imputing to a higher density would enable cost-effective implementation of genomic selection, provided imputation was accurate enough. Here, we describe the design of a low-density (12k) SNP chip and evaluate the accuracy of imputation from the 12k SNP genotypes to 50k SNP genotypes in the major Australian sheep breeds. In addition, the impact of imperfect imputation on genomic predictions was evaluated by comparing the accuracy of genomic predictions for 15 novel meat traits including carcass and meat quality and omega fatty acid traits in sheep, from 12k SNP genotypes, imputed 50k SNP genotypes and real 50k SNP genotypes. The 12k chip design included 12 223 SNPs with a high minor allele frequency that were selected with intermarker spacing of 50-475 kb. SNPs for parentage and horned or polled tests also were represented. Chromosome ends were enriched with SNPs to reduce edge effects on imputation. The imputation performance of the 12k SNP chip was evaluated using 50k SNP genotypes of 4642 animals from six breeds in three different scenarios: (1) within breed, (2) single breed from multibreed reference and (3) multibreed from a single-breed reference. The highest imputation accuracies were found with scenario 2, whereas scenario 3 was the worst, as expected. Using scenario 2, the average imputation accuracy in Border Leicester, Polled Dorset, Merino, White Suffolk and crosses was 0.95, 0.95, 0.92, 0.91 and 0.93 respectively. Imputation scenario 2 was used to impute 50k genotypes for 10 396 animals with novel meat trait phenotypes to compare genomic prediction accuracy using genomic best linear unbiased prediction (GBLUP) with real and imputed 50k genotypes. The weighted mean imputation accuracy achieved was 0.92. The average accuracy of genomic estimated breeding values (GEBVs) based on only 12k data was 0.08 across traits and breeds, but accuracies varied widely. The mean GBLUP accuracies with imputed 50k data more than doubled to 0.21. Accuracies of genomic prediction were very similar for imputed and real 50k genotypes. There was no apparent impact on accuracy of GEBVs as a result of using imputed rather than real 50k genotypes, provided imputation accuracy was >90%.

Estimation of genomic breeding values for residual feed intake in a multibreed cattle population.

Residual feed intake (RFI) is a measure of the efficiency of animals in feed utilization. The accuracies of GEBV for RFI could be improved by increasing the size of the reference population. Combining RFI records of different breeds is a way to do that. The aims of this study were to 1) develop a method for calculating GEBV in a multibreed population and 2) improve the accuracies of GEBV by using SNP associated with RFI. An alternative method for calculating accuracies of GEBV using genomic BLUP (GBLUP) equations is also described and compared to cross-validation tests. The dataset included RFI records and 606,096 SNP genotypes for 5,614 Bos taurus animals including 842 Holstein heifers and 2,009 Australian and 2,763 Canadian beef cattle. A range of models were tested for combining genotype and phenotype information from different breeds and the best model included an overall effect of each SNP, an effect of each SNP specific to a breed, and a small residual polygenic effect defined by the pedigree. In this model, the Holsteins and some Angus cattle were combined into 1 "breed class" because they were the only cattle measured for RFI at an early age (6-9 mo of age) and were fed a similar diet. The average empirical accuracy (0.31), estimated by calculating the correlation between GEBV and actual phenotypes divided by the square root of estimated heritability in 5-fold cross-validation tests, was near to that expected using the GBLUP equations (0.34). The average empirical and expected accuracies were 0.30 and 0.31, respectively, when the GEBV were estimated for each breed separately. Therefore, the across-breed reference population increased the accuracy of GEBV slightly, although the gain was greater for breeds with smaller number of individuals in the reference population (0.08 in Murray Grey and 0.11 in Hereford for empirical accuracy). In a second approach, SNP that were significantly (P < 0.001) associated with RFI in the beef cattle genomewide association studies were used to create an auxiliary genomic relationship matrix for estimating GEBV in Holstein heifers. The empirical (and expected) accuracy of GEBV within Holsteins increased from 0.33 (0.35) to 0.39 (0.36) and improved even more to 0.43 (0.50) when using a multibreed reference population. Therefore, a multibreed reference population is a useful resource to find SNP with a greater than average association with RFI in 1 breed and use them to estimate GEBV in another breed.

A marker-derived gene network reveals the regulatory role of PPARGC1A, HNF4G, and FOXP3 in intramuscular fat deposition of beef cattle.

High intramuscular fat (IMF) awards price premiums to beef producers and is associated with meat quality and flavor. Studying gene interactions and pathways that affect IMF might unveil causative physiological mechanisms and inform genomic selection, leading to increased accuracy of predictions of breeding value. To study gene interactions and pathways, a gene network was derived from genetic markers associated with direct measures of IMF, other fat phenotypes, feedlot performance, and a number of meat quality traits relating to body conformation, development, and metabolism that might be plausibly expected to interact with IMF biology. Marker associations were inferred from genomewide association studies (GWAS) based on high density genotypes and 29 traits measured on 10,181 beef cattle animals from 3 breed types. For the network inference, SNP pairs were assessed according to the strength of the correlation between their additive association effects across the 29 traits. The co-association inferred network was formed by 2,434 genes connected by 28,283 edges. Topological network parameters suggested a highly cohesive network, in which the genes are strongly functionally interconnected. Pathway and network analyses pointed towards a trio of transcription factors (TF) as key regulators of carcass IMF: PPARGC1A, HNF4G, and FOXP3. Importantly, none of these genes would have been deemed as significantly associated with IMF from the GWAS. Instead, a total of 313 network genes show significant co-association with the 3 TF. These genes belong to a wide variety of biological functions, canonical pathways, and genetic networks linked to IMF-related phenotypes. In summary, our GWAS and network predictions are supported by the current literature and suggest a cooperative role for the 3 TF and other interacting genes including CAPN6, STC2, MAP2K4, EYA1, COPS5, XKR4, NR2E1, TOX, ATF1, ASPH, TGS1, and TTPA as modulators of carcass and meat quality traits in beef cattle.

Evidence for pleiotropism and recent selection in the PLAG1 region in Australian Beef cattle.

A putative functional mutation (rs109231213) near PLAG1 (BTA14) associated with stature was studied in beef cattle. Data from 8199 Bos taurus, Bos indicus and Tropical Composite cattle were used to test the associations between rs109231213 and various phenotypes. Further, 23 496 SNPs located on BTA14 were tested for association with these phenotypes, both independently and fitted together with rs109231213. The C allele of rs109231213 significantly increased hip height, weight, net food intake, age at puberty in males and females and decreased IGF-I concentration in blood and fat depth. When rs109231213 was fitted as a fixed effect in the model, there was an overall reduction in associations between other SNPs and these traits but some SNPs remained associated (P < 10(-4) ). Frequency of the mutant C allele of rs109231213 differed among B. indicus (0.52), B. taurus (0.96) and Tropical Composite (0.68). Most chromosomes carrying the C allele had the same surrounding 10 SNP haplotype, probably because the C allele was introgressed into Brahman from B. taurus cattle. A region of reduced heterozygosity surrounds the C allele; this is small in B. taurus but 20 Mb long in Brahmans, indicating recent and strong selection for the mutant allele. Thus, the C allele appears to mark a mutation that has been selected almost to fixation in the B. taurus breeds studied here and introduced into Brahman cattle during grading up and selected to a frequency of 0.52 despite its negative effects on fertility.

Accuracy of prediction of genomic breeding values for residual feed intake and carcass and meat quality traits in Bos taurus, Bos indicus, and composite beef cattle.

The aim of this study was to assess the accuracy of genomic predictions for 19 traits including feed efficiency, growth, and carcass and meat quality traits in beef cattle. The 10,181 cattle in our study had real or imputed genotypes for 729,068 SNP although not all cattle were measured for all traits. Animals included Bos taurus, Brahman, composite, and crossbred animals. Genomic EBV (GEBV) were calculated using 2 methods of genomic prediction [BayesR and genomic BLUP (GBLUP)] either using a common training dataset for all breeds or using a training dataset comprising only animals of the same breed. Accuracies of GEBV were assessed using 5-fold cross-validation. The accuracy of genomic prediction varied by trait and by method. Traits with a large number of recorded and genotyped animals and with high heritability gave the greatest accuracy of GEBV. Using GBLUP, the average accuracy was 0.27 across traits and breeds, but the accuracies between breeds and between traits varied widely. When the training population was restricted to animals from the same breed as the validation population, GBLUP accuracies declined by an average of 0.04. The greatest decline in accuracy was found for the 4 composite breeds. The BayesR accuracies were greater by an average of 0.03 than GBLUP accuracies, particularly for traits with known genes of moderate to large effect mutations segregating. The accuracies of 0.43 to 0.48 for IGF-I traits were among the greatest in the study. Although accuracies are low compared with those observed in dairy cattle, genomic selection would still be beneficial for traits that are hard to improve by conventional selection, such as tenderness and residual feed intake. BayesR identified many of the same quantitative trait loci as a genomewide association study but appeared to map them more precisely. All traits appear to be highly polygenic with thousands of SNP independently associated with each trait.

Validation of single nucleotide polymorphisms associated with milk production traits in dairy cattle.

Single nucleotide polymorphism (SNP) associations with milk production traits found to be significant in different screening experiments, including SNP in genes hypothesized to be in gene pathways affecting milk production, were tested in a validation population to confirm their association. In total, 423 SNP were genotyped across 411 Holstein bulls, and their association with 6 milk production traits--Australian Selection Index (indicating the profitability of an animal's milk production), protein, fat, and milk yields, and protein and fat composition--were tested using single SNP regressions. Seventy-two SNP were significantly associated with one or more of the traits; their effects were in the same direction as in the screening experiment and therefore their association was considered validated. An over-representation of SNP (43 of the 423) on chromosome 20 was observed, including a SNP in the growth hormone receptor gene previously published as having an association with protein composition and protein and milk yields. The association with protein composition was confirmed in this experiment, but not the association with protein and milk yields. A multiple SNP regression analysis for all SNP on chromosome 20 was performed for all 6 traits, which revealed that this mutation was not significantly associated with any of the milk production traits and that at least 2 other quantitative trait loci were present on chromosome 20.

Short communication: Genomic selection using a multi-breed, across-country reference population.

Three breeds (Fleckvieh, Holstein, and Jersey) were included in a reference population, separately and together, to assess the accuracy of prediction of genomic breeding values in single-breed validation populations. The accuracy of genomic selection was defined as the correlation between estimated breeding values, calculated using phenotypic data, and genomic breeding values. The Holstein and Jersey populations were from Australia, whereas the Fleckvieh population (dual-purpose Simmental) was from Austria and Germany. Both a BLUP with a multi-breed genomic relationship matrix (GBLUP) and a Bayesian method (BayesA) were used to derive the prediction equations. The hypothesis tested was that having a multi-breed reference population increased the accuracy of genomic selection. Minimal advantage existed of either GBLUP or BayesA multi-breed genomic evaluations over single-breed evaluations. However, when the goal was to predict genomic breeding values for a breed with no individuals in the reference population, using 2 other breeds in the reference was generally better than only 1 breed.

A genome-wide association study of meat and carcass traits in Australian cattle.

Chromosomal regions containing DNA variation affecting the traits intramuscular fat percentage (IMF), meat tenderness measured as peak force to shear the LM (LLPF), and rump fat measured at the sacro-iliac crest in the chiller (CHILLP8) were identified using a set of 53,798 SNP genotyped on 940 taurine and indicine cattle sampled from a large progeny test experiment. Of these SNP, 87, 64, and 63 were significantly (P < 0.001) associated with the traits IMF, LLPF, and CHILLP8, respectively. A second, nonoverlapping sample of 1,338 taurine and indicine cattle from the same large progeny test experiment genotyped for 335 SNP, including as a positive control the calpastatin (CAST) c.2832A > G SNP, was used to confirm these locations. In total, 37 SNP were significantly (P < 0.05) associated with the same trait and with the same favorable homozygote in both data sets, representing 27 chromosomal regions. For the trait IMF, the effect of SNP in the confirmation data set was predicted from the discovery set by multiplying the estimated allele effect of each SNP in the discovery set by the number of copies of the reference allele of each SNP in the confirmation set. These weighted effects were then summed over all SNP to generate a molecular breeding value (MBV) for each animal in the confirmation data set. Using a bivariate analysis of MBV and IMF phenotypes of animals in the confirmation set, a panel of 14 SNP explained 5.6 and 15.6% of the phenotypic and genetic variance of IMF, respectively, in the confirmation data set. The amount of variation did not increase as more SNP were added to the MBV and instead decreased to 1.2 and 3.8% of the phenotypic and genetic variance of IMF, respectively, when 329 SNP were included in the analysis.

Detection of chromosome segments of zebu and taurine origin and their effect on beef production and growth.

Cattle in breeds formed by recent crossing of Bos taurus (Bt) and Bos indicus (Bi) subspecies should contain chromosomes that are a composite of Bt and Bt segments. Using data from a 50K SNP chip, we were able to identify whether a chromosome segment of 11 SNP in a composite animal descended from a Bt or a Bi ancestor. When the method was tested in purebred Bt or Brahman cattle, about 94% of segments were assigned correctly. About 10% of the genome in Australian Brahman cattle appears to be of Bt origin, as might be expected from their history. We then examined the effect of the origin of each chromosome segment on BW in a population of 515 Bt × Bi composite cattle and found 67 chromosome segments with a significant (P<0.01) effect. We confirmed these effects by examining these 67 segments in a population of Brahman cattle and in a population of mixed breeds including composite breeds such as Santa Gertrudis and Brahman cattle. About 66% of the 67 segments had an effect in the same direction in the confirmation analyses as in the discovery population. However, the effect on BW and other traits of chromosome segment origin is small, indicating that we had low power to detect these effects with the number of animals available. Consequently, when chromosome segment origin was used in genomic selection to predict BW, the accuracy was low (0.08). Chromosome segments that had a positive effect on BW tend to be at greater frequency in composite breeds than chromosome segments with a negative effect on BW.

Genome-wide association studies for feedlot and growth traits in cattle.

A genome wide-association study for production traits in cattle was carried out using genotype data from the 10K Affymetrix (Santa Clara, CA) and the 50K Illumina (San Diego, CA) SNP chips. The results for residual feed intake (RFI), BW, and hip height in 3 beef breed types (Bos indicus, Bos taurus, and B. indicus × B. taurus), and for stature in dairy cattle, are presented. The aims were to discover SNP associated with all traits studied, but especially RFI, and further to test the consistency of SNP effects across different cattle populations and breed types. The data were analyzed within data sets and within breed types by using a mixed model and fitting 1 SNP at a time. In each case, the number of significant SNP was more than expected by chance alone. A total of 75 SNP from the reference population with 50K chip data were significant (P < 0.001) for RFI, with a false discovery rate of 68%. These 75 SNP were mapped on 24 different BTA. Of the 75 SNP, the 9 most significant SNP were detected on BTA 3, 5, 7, and 8, with P ≤ 6.0 × 10(-5). In a population of Angus cattle divergently selected for high and low RFI and 10K chip data, 111 SNP were significantly (P < 0.001) associated with RFI, with a false discovery rate of 7%. Approximately 103 of these SNP were therefore likely to represent true positives. Because of the small number of SNP common to both the 10K and 50K SNP chips, only 27 SNP were significantly (P < 0.05) associated with RFI in the 2 populations. However, other chromosome regions were found that contained SNP significantly associated with RFI in both data sets, although no SNP within the region showed a consistent effect on RFI. The SNP effects were consistent between data sets only when estimated within the same breed type.

Multivariate analysis of a genome-wide association study in dairy cattle.

Multiple-trait genome-wide association study (GWAS) analyses were compared with single-trait GWAS for power to discover and subsequently validate genetic markers (single nucleotide polymorphisms; SNP) associated with dairy traits. The SNP associations were discovered in 1 Holstein population and validated in both a Holstein population consisting of bulls younger than those in the discovery population and a Jersey population. The multivariate methods used were a principal component analysis and a series of bivariate analyses. The statistical power of detecting associations using multiple-trait GWAS was as good as or better than that of the best single-trait GWAS. Additional SNP associations were found with the multivariate methods that had not been discovered in the single-trait analyses; this was achieved without an increase in the false discovery rate. From the multivariate analysis, 4 common pleiotropic patterns were identified among the putative quantitative trait loci (QTL) affecting the Australian selection index. These patterns could be interpreted as a primary effect of the putative QTL on 1 or more milk components and secondary effects on other components. The multivariate analysis did not appear to increase the precision with which putative QTL were mapped.

A validated genome-wide association study in 2 dairy cattle breeds for milk production and fertility traits using variable length haplotypes.

Genome-wide association studies (GWAS) were used to discover genomic regions explaining variation in dairy production and fertility traits. Associations were detected with either single nucleotide polymorphism (SNP) markers or haplotypes of SNP alleles. An across-breed validation strategy was used to narrow the genomic interval containing causative mutations. There were 39,048 SNP tested in a discovery population of 780 Holstein sires and validated in 386 Holsteins and 364 Jersey sires. Previously identified mutations affecting milk production traits were confirmed. In addition, several novel regions were identified, including a putative quantitative trait loci for fertility on chromosome 18 that was detected only using haplotypes greater than 3 SNP long. It was found that the precision of quantitative trait loci mapping increased with haplotype length as did the number of validated haplotypes discovered, especially across breed. Promising candidate genes have been identified in several of the validated regions.

A quantitative trait locus for faecal worm egg and blood eosinophil counts on chromosome 23 in Australian goats.

Three microsatellite markers on goat chromosome 23 adjacent to the MHC were used to test for quantitative trait loci (QTL) affecting faecal worm egg count (WEC) and leukocyte traits in ten Australian Angora and twelve Australian Cashmere half-sib families (n = 16-57 per family). Data were collected from 280 Angora and 347 Cashmere kids over a 3- and 4-year period. A putative QTL affecting trichostrongyle WEC was found in two small families at the 5% chromosome-wise threshold level. The biggest QTL effect for WEC of 1.65 standard deviations (sigma(p)) was found within the region of OarCP73-BM1258. A significant QTL affecting blood eosinophil counts at the 1% chromosome-wise threshold level was detected at marker BM1258 (at 26 cM) in two Angora and Cashmere families. The magnitude of the putative QTL was 0.69 and 0.85 sigma(p) in Angora and Cashmere families, respectively. Due to the comparatively low power of the study these findings should be viewed as indicative rather than definitive.