Evaluating the use of statistical and machine learning methods for estimating breed composition of purebred and crossbred animals in thirteen cattle breeds using genomic information.

Introduction: The ability to accurately predict breed composition using genomic information has many potential uses including increasing the accuracy of genetic evaluations, optimising mating plans and as a parameter for genotype quality control. The objective of the present study was to use a database of genotyped purebred and crossbred cattle to compare breed composition predictions using a freely available software, Admixture, with those from a single nucleotide polymorphism Best Linear Unbiased Prediction (SNP-BLUP) approach; a supplementary objective was to determine the accuracy and general robustness of low-density genotype panels for predicting breed composition. Methods: All animals had genotype information on 49,213 autosomal single nucleotide polymorphism (SNPs). Thirteen breeds were included in the analysis and 500 purebred animals per breed were used to establish the breed training populations. Accuracy of breed composition prediction was determined using a separate validation population of 3,146 verified purebred and 4,330 two and three-way crossbred cattle. Results: When all 49,213 autosomal SNPs were used for breed prediction, a minimal absolute mean difference of 0.04 between Admixture vs. SNP-BLUP breed predictions was evident. For crossbreds, the average absolute difference in breed prediction estimates generated using SNP-BLUP and Admixture was 0.068 with a root mean square error of 0.08. Breed predictions from low-density SNP panels were generated using both SNP-BLUP and Admixture and compared to breed prediction estimates using all 49,213 SNPs (representing the gold standard). Breed composition estimates of crossbreds required more SNPs than predicting the breed composition of purebreds. SNP-BLUP required ≥3,000 SNPs to predict crossbred breed composition, but only 2,000 SNPs were required to predict purebred breed status. The absolute mean (standard deviation) difference across all panels <2,000 SNPs was 0.091 (0.054) and 0.315 (0.316) when predicting the breed composition of all animals using Admixture and SNP-BLUP, respectively compared to the gold standard prediction. Discussion: Nevertheless, a negligible absolute mean (standard deviation) difference of 0.009 (0.123) in breed prediction existed between SNP-BLUP and Admixture once ≥3,000 SNPs were considered, indicating that the prediction of breed composition could be readily integrated into SNP-BLUP pipelines used for genomic evaluations thereby avoiding the necessity for a stand-alone software.

Population structure and breed composition prediction in a multi-breed sheep population using genome-wide single nucleotide polymorphism genotypes.

Knowledge of population structure and breed composition of a population can be advantageous for a number of reasons; these include designing optimal (cross)breeding strategies in order to maximise non-additive genetic effects, maintaining flockbook integrity by authenticating animals being registered and as a quality control measure in the genotyping process. The objectives of the present study were to 1) describe the population structure of 24 sheep breeds, 2) quantify the breed composition of both flockbook-recorded and crossbred animals using single nucleotide polymorphism BLUP (SNP-BLUP), and 3) quantify the accuracy of breed composition prediction from low-density genotype panels containing between 2000 and 6000 SNPs. In total, 9334 autosomal SNPs on 11 144 flockbook-recorded animals and 1172 crossbred animals were used. The population structure of all breeds was characterised by principal component analysis (PCA) as well as the pairwise breed fixation index (Fst). The total number of animals, all of which were purebred, included in the calibration population for SNP-BLUP was 2579 with the number of animals per breed ranging from 9 to 500. The remaining 9559 flockbook-recorded animals, composite breeds and crossbred animals represented the test population; three breeds were excluded from breed composition prediction. The breed composition predicted using SNP-BLUP with 9334 SNPs was considered the gold standard prediction. The pairwise breed Fst ranged from 0.040 (between the Irish Blackface and Scottish Blackface) to 0.282 (between the Border Leicester and Suffolk). Principal component analysis revealed that the Suffolk from Ireland and the Suffolk from New Zealand formed distinct, non-overlapping clusters. In contrast, the Texel from Ireland and that from New Zealand formed integrated, overlapping clusters. Composite animals such as the Belclare clustered close to its founder breeds (i.e., Finn, Galway, Lleyn and Texel). When all 9334 SNPs were used to predict breed composition, an animal that had a majority breed proportion predicted to be ≥0.90 was defined as purebred for the present study. As the panel density decreased, the predicted breed proportion threshold, used to identify animals as purebred, also decreased (≥0.85 with 6000 SNPs to ≥0.60 with 2000 SNPs). In all, results from the study suggest that breed composition for purebred and crossbred animals can be determined with SNP-BLUP using ≥5000 SNPs.

Reaffirmation of known major genes and the identification of novel candidate genes associated with carcass-related metrics based on whole genome sequence within a large multi-breed cattle population.

BACKGROUND: The high narrow sense heritability of carcass traits suggests that the underlying additive genetic potential of an individual should be strongly correlated with both animal carcass quality and quantity, and therefore, by extension, carcass value. Therefore, the objective of the present study was to detect genomic regions associated with three carcass traits, namely carcass weight, conformation and fat cover, using imputed whole genome sequence in 28,470 dairy and beef sires from six breeds with a total of 2,199,926 phenotyped progeny. RESULTS: Major genes previously associated with carcass performance were identified, as well as several putative novel candidate genes that likely operate both within and across breeds. The role of MSTN in carcass performance was re-affirmed with the segregating Q204X mutation explaining 1.21, 1.11 and 5.95% of the genetic variance in carcass weight, fat and conformation, respectively in the Charolais population. In addition, a genomic region on BTA6 encompassing the NCAPG/LCORL locus, which is a known candidate locus associated with body size, was associated with carcass weight in Angus, Charolais and Limousin. Novel candidate genes identified included ZFAT in Angus, and SLC40A1 and the olfactory gene cluster on BTA15 in Charolais. Although the majority of associations were breed specific, associations that operated across breeds included SORCS1 on BTA26, MCTP2 on BTA21 and ARL15 on BTA20; these are of particular interest due to their potential informativeness in across-breed genomic evaluations. Genomic regions affecting all three carcass traits were identified in each of the breeds, although these were mainly concentrated on BTA2 and BTA6, surrounding MSTN and NCAPG/LCORL, respectively. This suggests that although major genes may be associated with all three carcass traits, the majority of genes containing significant variants (unadjusted p-value < 10- 4) may be trait specific associations of small effect. CONCLUSIONS: Although plausible novel candidate genes were identified, the proportion of variance explained by these candidates was minimal thus reaffirming that while carcass performance may be affected by major genes in the form of MSTN and NCAPG/LCORL, the majority of variance is attributed to the additive (and possibly multiplicative) effect of many polymorphisms of small effect.

Anti-Müllerian hormone in grazing dairy cows: Identification of factors affecting plasma concentration, relationship with phenotypic fertility, and genome-wide associations.

The objectives of this study were to (1) characterize the distribution and variability of plasma anti-Müllerian hormone (AMH) concentration; (2) evaluate factors associated with phenotypic variation in plasma AMH; (3) examine the associations between categories of plasma AMH and reproductive outcomes [pregnancy to first artificial insemination (P/AI), and pregnancy rates within 21, 42, and 84 d after the mating start date (MSD)]; (4) estimate pedigree and genomic heritability for plasma AMH; and (5) identify and validate SNP associated with phenotypic variation in plasma AMH. Plasma AMH concentration (pg/mL) was determined from a blood sample collected (mean ± standard deviation) 10 ± 2 d after first insemination at detected estrus (IDE) in 2,628 first- and second-parity Irish dairy cows. Overall, plasma AMH had a positively skewed distribution with mean (± standard deviation), median, minimum, and maximum concentrations of 326 ± 231, 268, 15, and 2,863 pg/mL, respectively. Plasma AMH was greatest for Jersey, followed by Holstein × Jersey, Holstein × Norwegian Red, and Holstein cows (410, 332, 284, and 257 pg/mL, respectively). Second-parity cows had greater plasma AMH than first-parity cows (333 vs. 301 pg/mL, respectively). Samples collected at 7 and 8 d after first IDE had lesser plasma AMH than those collected on d 9, 10, 11, 12, and 13 after first IDE (291 and 297 vs. 317, 319, 331, 337, and 320 pg/mL). Plasma AMH was not associated with either body condition score at first IDE or the interval from calving to MSD. Cows were categorized into low (≤150 pg/mL; n = 526; lowest 20%), intermediate (>150 to ≤461 pg/mL; n = 1,576; intermediate 60%), and high AMH (>461 pg/mL; n = 526; highest 20%) groups based on plasma AMH, and associations with reproductive outcomes were tested. Cows with high and intermediate plasma AMH had 1.42- and 1.51-times-greater odds of becoming pregnant within 84 d after the MSD than those with low plasma AMH (90.3 and 90.8 vs. 86.8%, respectively); however, P/AI and pregnancy rate within 21 and 42 d after the MSD did not differ among AMH categories. Plasma AMH was moderately heritable (pedigree heritability of 0.40 ± 0.06 and genomic heritability of 0.45 ± 0.05), and 68 SNP across Bos taurus autosomes 7 and 11 were associated with phenotypic variation in plasma AMH. Out of 68 SNP, 42 were located in a single quantitative trait locus on Bos taurus autosome 11 that harbored 6 previously identified candidate genes (NR5A1, HSPA5, CRB2, DENND1A, NDUFA8, and PTGS) linked to fertility-related phenotypes in dairy cows.

Candidate genes associated with the heritable humoral response to Mycobacterium avium ssp. paratuberculosis in dairy cows have factors in common with gastrointestinal diseases in humans.

Infection of cattle with bovine paratuberculosis (i.e., Johne's disease) is caused by Mycobacterium avium ssp. paratuberculosis (MAP) and results in a chronic incurable gastroenteritis. This disease, which has economic ramifications for the cattle industry, is increasing in detected prevalence globally; subclinically infected animals can silently shed the bacterium into the environment for years, exposing contemporaries and hampering disease-control programs. The objective of the present study was to first quantify the genetic parameters for humoral response to MAP in dairy cattle. This was followed by a genome-based association analysis and subsequent downstream bioinformatic analyses from imputed whole genome sequence SNP data. After edits, ELISA test records were available on 136,767 cows; analyses were also undertaken on a subset of 33,818 of these animals from herds with at least 5 MAP ELISA-positive cows, with at least 1 of those positive cows being homebred. Variance components were estimated using univariate animal and sire linear mixed models. The heritability calculated from the animal model for humoral response to MAP using alternative phenotype definitions varied from 0.02 (standard error = 0.003) to 0.05 (standard error = 0.008). The genome-based associations were undertaken within a mixed model framework using weighted deregressed estimated breeding values as a dependent variable on 1,883 phenotyped animals that were ≥87.5% Holstein-Friesian. Putative susceptibility quantitative trait loci (QTL) were identified on Bos taurus autosome 1, 3, 5, 6, 8, 9, 10, 11, 13, 14, 18, 21, 23, 25, 26, 27, and 29; mapping the most significant SNP to genes within and overlapping these QTL revealed that the most significant associations were with the 10 functional candidate genes KALRN, ZBTB20, LPP, SLA2, FI3A1, LRCH3, DNAJC6, ZDHHC14, SNX1, and HAS2. Pathway analysis failed to reveal significantly enriched biological pathways, when both bovine-specific pathway data and human ortholog data were taken into account. The existence of genetic variation for MAP susceptibility in a large data set of dairy cows signifies the potential of breeding programs for reducing MAP susceptibility. Furthermore, the identification of susceptible QTL facilitates greater biological understanding of bovine paratuberculosis and potential therapeutic targets for future investigation. The novel molecular similarities identified between bovine paratuberculosis and human inflammatory bowel disease suggest potential for human therapeutic interventions to be translated to veterinary medicine and vice versa.

Variance components for bovine tuberculosis infection and multi-breed genome-wide association analysis using imputed whole genome sequence data.

Bovine tuberculosis (bTB) is an infectious disease of cattle generally caused by Mycobacterium bovis, a bacterium that can elicit disease humans. Since the 1950s, the objective of the national bTB eradication program in Republic of Ireland was the biological extinction of bTB; that purpose has yet to be achieved. Objectives of the present study were to develop the statistical methodology and variance components to undertake routine genetic evaluations for resistance to bTB; also of interest was the detection of regions of the bovine genome putatively associated with bTB infection in dairy and beef breeds. The novelty of the present study, in terms of research on bTB infection, was the use of beef breeds in the genome-wide association and the utilization of imputed whole genome sequence data. Phenotypic bTB data on 781,270 animals together with imputed whole genome sequence data on 7,346 of these animals' sires were available. Linear mixed models were used to quantify variance components for bTB and EBVs were validated. Within-breed and multi-breed genome-wide associations were undertaken using a single-SNP regression approach. The estimated genetic standard deviation (0.09), heritability (0.12), and repeatability (0.30) substantiate that genetic selection help to eradicate bTB. The multi-breed genome-wide association analysis identified 38 SNPs and 64 QTL regions associated with bTB infection; two QTL regions (both on BTA23) identified in the multi-breed analysis overlapped with the within-breed analyses of Charolais, Limousin, and Holstein-Friesian. Results from the association analysis, coupled with previous studies, suggest bTB is controlled by an infinitely large number of loci, each having a small effect. The methodology and results from the present study will be used to develop national genetic evaluations for bTB in the Republic of Ireland. In addition, results can also be used to help uncover the biological architecture underlying resistance to bTB infection in cattle.

The relationship between anogenital distance and fertility, and genome-wide associations for anogenital distance in Irish Holstein-Friesian cows.

The evaluation of anogenital distance (AGD), the distance from the center of the anus to base of the clitoris, as a potential fertility trait for genetic selection in dairy cows has generated recent interest. The objectives of this cross-sectional observational study were to (1) characterize the distribution and variability of AGD, (2) determine factors associated with AGD, (3) estimate heritability for AGD, (4) identify single nucleotide polymorphisms (SNP) associated with phenotypic variation of AGD, and (5) validate the relationship between categories of AGD and fertility in Irish Holstein-Friesian cows. Anogenital distance was measured using digital calipers in 1,180 Holstein cows (mean ± standard deviation: 225 ± 79 d in milk) from 10 dairy herds located in Munster, Ireland. In addition, age (yr), weight (kg), height at hip (cm), and body condition score (BCS) at the time of AGD measurement were determined in a subset of 281 cows. Genotype information available from 908 cows was subsequently imputed to the Illumina Bovine High Density BeadChip (Illumina Inc., San Diego, CA) for genome-wide association analysis of phenotypic variation in AGD. Overall, AGD had a normal distribution and high variability (mean ± standard deviation; 119.2 ± 11.6 mm). Anogenital distance was weakly but positively associated with cow age, hip height, and body weight, and negatively associated with BCS; the phenotypic variation in AGD that was explainable by these variables was small (coefficient of determination; R2 = 0.09, 0.06, 0.10, and 0.02, respectively). The estimated heritability for AGD was 0.37 (standard error of mean ± 0.08). Six SNP of suggestive significance were identified on Bos taurus autosomes 6, 15, 20, and 26; however, none of these SNP was related to previously identified candidate genes for fertility. Cows were categorized into quartiles (Q1; 86 to 111 mm; n = 311, Q2; 112 to 120 mm; n = 330; Q3; 121 to 127 mm; n = 265, and Q4; 128 to 160 mm; n = 274) based on AGD and the association with reproductive outcomes examined (21-d submission rate, pregnancy to first AI, pregnancy rate within 21, 42 and 84-d after the farm mating start date, and number of times bred). None of the reproductive variables differed significantly between AGD categories. In summary, despite identification of high variability and moderate heritability for AGD in Irish Holstein-Friesian cows, reproductive outcomes did not differ between categories of AGD. This latter result differs from our previous finding of an inverse relationship between AGD and pregnancy outcomes in first- and second-parity Canadian Holstein cows, emphasizing the need to test and validate this new phenotype in diverse cow populations.

The relationship between serum insulin-like growth factor-1 (IGF-1) concentration and reproductive performance, and genome-wide associations for serum IGF-1 in Holstein cows.

The objectives of this study were to determine (1) factors associated with serum concentration of insulin-like growth factor-1 (IGF-1); (2) the relationship between serum IGF-1 concentration during the first week postpartum and ovarian cyclicity status by 35 d postpartum (DPP); (3) an optimum serum IGF-1 concentration threshold predictive of pregnancy to first artificial insemination (P/AI), including its diagnostic values; (4) the associations among categories of serum IGF-1 concentration and reproductive outcomes (P/AI and pregnancy risk up to 150 and 250 DPP); and (5) single nucleotide polymorphisms (SNP) associated with phenotypic variation in serum IGF-1 concentration in dairy cows. Serum IGF-1 concentration was determined at 7 (±2.4; ±standard error of the mean) DPP in 647 lactating Holstein cows (213 primiparous, 434 multiparous) from 7 herds in Alberta, Canada. The overall mean, median, minimum, and maximum serum IGF-1 concentrations during the first week postpartum were 37.8 (±1.23), 31.0, 20.0, and 225.0 ng/mL, respectively. Herd, age, parity, precalving body condition score, and season of blood sampling were all identified as factors associated with serum IGF-1 concentrations. Although serum IGF-1 concentration during the first week postpartum had no association with ovarian cyclicity status by 35 DPP in primiparous cows, it was greater in cyclic than in acyclic multiparous cows (32.2 vs. 27.4 ng/mL, respectively). The optimum serum IGF-1 thresholds predictive of P/AI were 85.0 ng/mL (sensitivity = 31.9%; specificity = 89.1%) and 31.0 ng/mL (sensitivity = 45.5%; specificity = 66.9%) for primiparous and multiparous cows, respectively. When cows were grouped into either high or low IGF-1 categories (greater or less than or equal to 85.0 ng/mL for primiparous cows and greater or less than or equal to 31.0 ng/mL for multiparous cows, respectively), primiparous cows with high IGF-1 had 4.43 times greater odds of P/AI and a tendency for higher pregnancy risk up to 150 DPP than those with low IGF-1, but not up to 250 DPP. Likewise, multiparous cows with high IGF-1 had 1.61 times greater odds of P/AI than those with low IGF-1. Pregnancy risk up to 150 and 250 DPP, however, did not differ between IGF-1 categories in multiparous cows. Moreover, 37 SNP across 10 Bos taurus autosomes were associated with variation in serum IGF-1 concentration, and 4 previously identified candidate genes related to fertility that were in linkage disequilibrium with some of these SNP were also identified.

The relationship between serum anti-Müllerian hormone concentrations and fertility, and genome-wide associations for anti-Müllerian hormone in Holstein cows.

The objectives of this study were to (1) evaluate factors associated with variation in circulating anti-Müllerian hormone (AMH) concentrations, (2) establish an optimum AMH threshold predictive of pregnancy to first artificial insemination (P/AI), (3) examine the relationship between AMH and fertility (P/AI, pregnancy loss between 30 and 60 d after artificial insemination, and pregnancy risk up to 250 d postpartum), and (4) identify quantitative trait loci associated with phenotypic variation of AMH concentrations in dairy cows. Serum AMH concentrations (pg/mL) were determined at 7 ± 2.4 d postpartum in 647 lactating Holstein cows (213 primiparous, 434 multiparous) from 1 research and 6 commercial dairy herds in Alberta, Canada. Of these, 589 cows were genotyped on the 26K Bovine BeadChip (Neogen Inc., Lincoln, NE) and subsequently imputed to the Illumina Bovine High Density BeadChip (Illumina, San Diego, CA) for genome-wide association analysis for variation in serum AMH concentrations. Factors associated with variation in serum AMH concentrations and the relationship between categories of AMH and aforementioned fertility outcomes were evaluated only in a subset of 460 cows that had a complete data set available. The overall mean (±standard error of the mean), median, minimum, and maximum AMH concentrations were 191.1 ± 6.3, 151.7, 13.9, and 1,879.0 pg/mL, respectively. The AMH concentrations were not associated with herd, precalving body condition score, postpartum week, and season of sampling; the lactation number, however, had a quadratic relationship with serum AMH concentrations (116.2, 204.9 204.5, and 157.9 pg/mL for first, second, third, and ≥fourth lactation, respectively). The optimum AMH threshold predictive of P/AI could not be established because the receiver operating characteristic curve analysis model was nonsignificant. Categories of AMH [low (<83.0 pg/mL; n = 92), intermediate (≥83.0 to ≤285.0 pg/mL; n = 276), and high (>285.0 pg/mL; n = 92) based on lowest 20%, intermediate 60%, and highest 20% serum AMH) had no associations with P/AI (34, 43, and 40%), pregnancy loss between 30 and 60 d after artificial insemination (20, 12, and 8%), or pregnancy risk up to 250 d postpartum. One candidate gene associated with AMH production [AMH gene on Bos taurus autosome (BTA) 7] and 4 candidate genes related to embryo development (SCAI and PPP6C genes on BTA11 and FGF18 and EEF2K genes on BTA20 and BTA25, respectively) were in linkage disequilibrium with single nucleotide polymorphisms associated with phenotypic variation in serum AMH in dairy cows.

Aneuploidy in dizygotic twin sheep detected using genome-wide single nucleotide polymorphism data from two commonly used commercial vendors.

Early detection of karyotype abnormalities, including aneuploidy, could aid producers in identifying animals which, for example, would not be suitable candidate parents. Genome-wide genetic marker data in the form of single nucleotide polymorphisms (SNPs) are now being routinely generated on animals. The objective of the present study was to describe the statistics that could be generated from the allele intensity values from such SNP data to diagnose karyotype abnormalities; of particular interest was whether detection of aneuploidy was possible with both commonly used genotyping platforms in agricultural species, namely the Applied BiosystemsTM AxiomTM and the Illumina platform. The hypothesis was tested using a case study of a set of dizygotic X-chromosome monosomy 53,X sheep twins. Genome-wide SNP data were available from the Illumina platform (11 082 autosomal and 191 X-chromosome SNPs) on 1848 male and 8954 female sheep and available from the AxiomTM platform (11 128 autosomal and 68 X-chromosome SNPs) on 383 female sheep. Genotype allele intensity values, either as their original raw values or transformed to logarithm intensity ratio (LRR), were used to accurately diagnose two dizygotic (i.e. fraternal) twin 53,X sheep, both of which received their single X chromosome from their sire. This is the first reported case of 53,X dizygotic twins in any species. Relative to the X-chromosome SNP genotype mean allele intensity values of normal females, the mean allele intensity value of SNP genotypes on the X chromosome of the two females monosomic for the X chromosome was 7.45 to 12.4 standard deviations less, and were easily detectable using either the AxiomTM or Illumina genotype platform; the next lowest mean allele intensity value of a female was 4.71 or 3.3 standard deviations less than the population mean depending on the platform used. Both 53,X females could also be detected based on the genotype LRR although this was more easily detectable when comparing the mean LRR of the X chromosome of each female to the mean LRR of their respective autosomes. On autopsy, the ovaries of the two sheep were small for their age and evidence of prior ovulation was not appreciated. In both sheep, the density of primordial follicles in the ovarian cortex was lower than normally found in ovine ovaries and primary follicle development was not observed. Mammary gland development was very limited. Results substantiate previous studies in other species that aneuploidy can be readily detected using SNP genotype allele intensity values generally already available, and the approach proposed in the present study was agnostic to genotype platform.

The impact of multi-generational genotype imputation strategies on imputation accuracy and subsequent genomic predictions.

The objective of the present study was to quantify, using simulations, the impact of successive generations of genotype imputation on genomic predictions. The impact of using a small reference population of true genotypes versus a larger reference population of imputed genotypes on the accuracy of genomic predictions was also investigated. After construction of a founder population, high-density (HD) genotypes ( = 43,500 single nucleotide polymorphisms, SNP) were simulated across 25 generations ( = 46,800 per generation); a low-density genotype panel ( = 3,000 SNP) was developed from these HD genotypes, which was then used to impute genotypes using 7 alternative imputation strategies. Both low (0.03) and moderately (0.35) heritable phenotypes were simulated. Direct genomic values (DGV) were estimated using imputed genotypes from the investigated scenarios and the accuracy of predicting the simulated true breeding values (TBV) were expressed relative to the accuracy when the true genotypes were used. Mean allele concordance rate and the rate of change in mean allele concordance per generation differed between the imputation strategies investigated. Imputation was most accurate when the true HD genotypes of sires and 50% of the dams of the generation being imputed were included in the reference population; the average allele concordance rate for this scenario across generations was 0.9707. The strongest correlation between the TBV and DGV of the last generation was when the reference population included sequentially imputed HD genotypes of all previous generations, plus the true HD genotypes of all sires of the previous generations (0.987 as efficient as when the true genotypes were used in the reference population). With a moderate heritability, the correlation between the TBV and the DGV using a small reference population of accurate genotypes were, on average, 0.07 units stronger compared to DGV generated using a larger population of imputed genotypes. When the heritability was low, the accuracy of genomic predictions benefited from a larger reference population, even if SNP were imputed. The impact on the accuracy of genomic predictions from the accumulation of imputation errors across generations indicates the need to routinely generate HD genotypes on influential animals to reduce the accumulation of imputation errors over generations.

Characterization of an X-chromosomal non-mosaic monosomy (59, X0) dairy heifer detected using routinely available single nucleotide polymorphism genotype data.

Evidence exists from a range of species on the impact of karyotype abnormalities on reproductive performance. Despite this, cytogenetic analyses of cattle, especially females, are not routinely undertaken. Genome-wide single nucleotide polymorphism (SNP) genotype data are now, however, routinely being generated in many species globally at a relatively low cost. The objective of the present study was to evaluate the potential of routinely available SNP genotype data to identify sex-chromosome aberrations using X chromosome monosomy 59,X0 as a case study for illustration. A single 2.5-yr old Holstein-Friesian heifer was detected with a mean allelic intensity of SNP on the X chromosome almost 17 standard deviations less than the mean of other genotyped females ( = 103,326). Following cytogenetic analysis (10 replicates by karyotyping and a further 140 by FISH), the female was deduced to be a non-mosaic 59,X0. The female had never produced a calf and, although gross examination revealed no physical abnormalities, she was smaller in size than expected based on her breed and age. Given the age of the animal at slaughter, the uterus and uterine tubes appeared immature and inactive. The oviduct appeared normal while the single ovary present contained a markedly reduced number of follicles. There was, however, some evidence of prior ovulation and formation of corpora lutea. The approach proposed in the present study to identify allosome aneuploidy from routinely available genotype data can be used to screen for such abnormalities at no additional cost to the breeder or producer.

Inference of population structure of purebred dairy and beef cattle using high-density genotype data.

Information on the genetic diversity and population structure of cattle breeds is useful when deciding the most optimal, for example, crossbreeding strategies to improve phenotypic performance by exploiting heterosis. The present study investigated the genetic diversity and population structure of the most prominent dairy and beef breeds used in Ireland. Illumina high-density genotypes (777 962 single nucleotide polymorphisms; SNPs) were available on 4623 purebred bulls from nine breeds; Angus (n=430), Belgian Blue (n=298), Charolais (n=893), Hereford (n=327), Holstein-Friesian (n=1261), Jersey (n=75), Limousin (n=943), Montbéliarde (n=33) and Simmental (n=363). Principal component analysis revealed that Angus, Hereford, and Jersey formed non-overlapping clusters, representing distinct populations. In contrast, overlapping clusters suggested geographical proximity of origin and genetic similarity between Limousin, Simmental and Montbéliarde and to a lesser extent between Holstein, Friesian and Belgian Blue. The observed SNP heterozygosity averaged across all loci was 0.379. The Belgian Blue had the greatest mean observed heterozygosity (HO=0.389) among individuals within breed while the Holstein-Friesian and Jersey populations had the lowest mean heterozygosity (HO=0.370 and 0.376, respectively). The correlation between the genomic-based and pedigree-based inbreeding coefficients was weak (r=0.171; P<0.001). Mean genomic inbreeding estimates were greatest for Jersey (0.173) and least for Hereford (0.051). The pair-wise breed fixation index (F st) ranged from 0.049 (Limousin and Charolais) to 0.165 (Hereford and Jersey). In conclusion, substantial genetic variation exists among breeds commercially used in Ireland. Thus custom-mating strategies would be successful in maximising the exploitation of heterosis in crossbreeding strategies.

Justification for setting the individual animal genotype call rate threshold at eighty-five percent.

Data quality of SNP arrays impacts the accuracy and precision of downstream data analyses. One such quality control measure often imposed is a threshold on individual animal call rate. Different call rate thresholds have been applied across studies; little is known, however, about the impact of these thresholds on the quality of the genotype data. The objective of the present study was to investigate the effect of different call rate thresholds on the integrity of the genotypes but also to quantify the contribution of different factors to the variability in animal call rate. Data included 142,342 samples genotyped on a custom Illumina genotype panel from 141,591 dairy and beef cattle; the number of Illumina SNP on the panel was 14,371. The mean animal call rate across all samples was 99.09%; 487 animals had both a low call rate (<99%) and a subsequent high call rate (≥99%) after resampling and regenotyping. Several factors were associated ( < 0.001) with individual call rate including animal sex, the sampling herd, the date of genotyping, the genotyping plate, and the plate well. The genotype and allele concordance between the genotypes of the 487 low- and high-call rate individuals improved at a diminishing rate as mean animal call rate increased. Mean genotype and allele concordance rates of 0.987 and 0.997, respectively, existed when animal call rate was between 85 and 90%, increasing to 0.998 and 0.999, respectively, when animal call rate was between 95 and <99%. The mean within-animal allele concordance rate of rare variants (i.e., minor allele frequency < 0.05) between low and high genotype call rate animals increased when animal call rate improved; an allele concordance rate of 1.00 was achieved when animal call rate was between 85 and <99%. The accuracy of imputation of the nonobserved genotypes in the low-call rate animals improved as animal call rate increased; the mean genotype concordance rate of the imputed nonobserved SNP was 0.41 when animal call rate was <40% but increased to 0.95 when animal call rate was between 95 and <99%. Parentage validation, determined by the count of opposing homozygotes in a parent-progeny pair, was unreliable when animal call rate was <85%. Therefore, to ensure the provision of high-quality genotypes while also considering the cost and inconvenience of resampling and regenotyping, we suggest a minimum animal call rate threshold of 85%.

Genetic characterization of a herd of the endangered Danish Jutland cattle.

In this paper we present results from a genetic characterization of a herd of the Danish Jutland cattle breed named the Kortegaard herd (n = 135; 57 males and 78 females). The herd is genotyped on the Bovine HD BeadChip microarray with 697,548 evenly spaced SNP across the bovine genome. The aim of the study was to characterize the genetic profile of the Kortegaard herd, which has been closed for several generations, by quantifying the degree of genetic homogeneity within the herd and to compare its genetic profile to that of other cattle breeds. A total of 868 animals from the Angus, Belgian Blue, Charolais, Friesian, Hereford, Holstein, Holstein-Friesian crosses, Limousin, and Simmental breeds was used for genetic profile comparisons. The level of genetic variation within the breeds were quantified by the expected heterozygosity (H(E)), observed heterozygosity (H(O)), average minor allele frequency (MAF), the degree of polymorphism, and runs of homozygosity (ROH), which are contiguous lengths of homozygous genotypes of varying length. Interestingly, the Kortegaard herd had the lowest within-breed genetic variation (lowest H(E), H(O), and MAF), showed moderate levels of short ROH (<5 Mb), and had the highest mean long ROH (>5 Mb) compared to all the other breeds. This is possibly due to recent consanguineous matings, a strong founder effect, and a lack of gene flow from other herds and breeds. We further examined whether the observed genetic patterns in the Kortegaard herd can be used to design breeding strategies for the preservation of the genetic pool by focusing on a subset of SNP outside homozygote regions. By calculating the pairwise identical-by-state between all possible matings, we designed a breeding plan that maximized heterozygosity in the short term. The benefits and limitations of such a breeding strategy are discussed.

Genome-wide association study for calving traits in Holstein-Friesian dairy cattle.

Dystocia and perinatal mortality are quantitative traits that significantly impact animal productivity and welfare. Their economic importance is reflected by their inclusion in the national breeding goals of many cattle populations. The genetic architecture that influences these traits, however, has still yet to be thoroughly defined. Regions of the bovine genome associated with calving difficulty (direct and maternal) and perinatal mortality were detected in this study using a Bayesian approach with 43 204 single nucleotide polymorphisms (SNPs) on up to 1970 Holstein-Friesian bulls. Several SNPs on chromosomes 5, 6, 11, 12, 17,18 and 28 were detected to be strongly associated with these calving performance traits. Novel genomic regions with previously reported associations with growth, stature, birth weight and bone morphology were identified in the present study as being associated with the three calving performance traits. Morphological abnormalities are a known contributor to perinatal mortality and the most significantly associated SNP for perinatal mortality in the present study was located in a region in linkage disequilibrium with the gene SLC26A7. This gene, SLC26A7, has similarities and colocalises with SLC4A2, which has previously been associated with osteoporosis and mortality in cattle populations. The HHIP gene that is known to be associated with stature in humans was strongly associated with direct calving difficulty in the present study; large calves are known to, on average, have a greater likelihood of dystocia. A stemloop microRNA, bta-mir-1256, on chromosome 12, involved in post-transcriptional regulation of gene expression was associated with maternal calving difficulty. Previously reported quantitative trait loci associated with calving performance traits in other populations were again identified in this study; with one genomic region on chromosome 18 supporting very strong evidence of an underlying causative mutation and accounting for 2.1% of the genetic variation in direct calving difficulty. Overlapping genomic regions associated with one or more of the calving traits were also detected substantiating the known genetic covariances existing between these traits. Moreover, some genomic regions were only associated with one of the calving traits implying the selective genomic breeding programs exploiting these regions could help resolve genetic antagonisms.