Search

Background: Temperament traits are of high importance across species. In humans, temperament or personality traits correlate with psychological traits and psychiatric disorders. In cattle, they impact animal welfare, product quality and human safety, and are therefore of direct commercial importance. We hypothesized that genetic factors that contribute to variation in temperament among individuals within a species will be shared between humans and cattle. Using imputed whole-genome sequence data from 9223 beef cattle from three cohorts, a series of genome-wide association studies was undertaken on cattle flight time, a temperament phenotype measured as the time taken for an animal to cover a short-fixed distance after release from an enclosure. We also investigated the association of cattle temperament with polymorphisms in bovine orthologs of risk genes for neuroticism, schizophrenia, autism spectrum disorders (ASD), and developmental delay disorders in humans., Results: Variants with the strongest associations were located in the bovine orthologous region that is involved in several behavioural and cognitive disorders in humans. These variants were also partially validated in independent cattle cohorts. Genes in these regions (BARHL2, NDN, SNRPN, MAGEL2, ABCA12, KIFAP3, TOPAZ1, FZD3, UBE3A, and GABRA5) were enriched for the GO term neuron migration and were differentially expressed in brain and pituitary tissues in humans. Moreover, variants within 100 kb of ASD susceptibility genes were associated with cattle temperament and explained 6.5% of the total additive genetic variance in the largest cattle cohort. The ASD genes with the most significant associations were GABRB3 and CUL3. Using the same 100 kb window, a weak association was found with polymorphisms in schizophrenia risk genes and no association with polymorphisms in neuroticism and developmental delay disorders risk genes., Conclusions: Our analysis showed that genes identified in a meta-analysis of cattle temperament contribute to neuron development functions and are differentially expressed in human brain tissues. Furthermore, some ASD susceptibility genes are associated with cattle temperament. These findings provide evidence that genetic control of temperament might be shared between humans and cattle and highlight the potential for future analyses to leverage results between species.

Reducing the incidence of both the degree of assistance required at calving, as well as the extent of perinatal mortality (PM) has both economic and societal benefits. The existence of heritable genetic variability in both traits signifies the presence of underlying genomic variability. The objective of the present study was to locate regions of the genome, and by extension putative genes and mutations, that are likely to be underpinning the genetic variability in direct calving difficulty (DCD), maternal calving difficulty (MCD), and PM. Imputed whole-genome single-nucleotide polymorphism (SNP) data on up to 8,304 Angus (AA), 17,175 Charolais (CH), 16,794 Limousin (LM), and 18,474 Holstein-Friesian (HF) sires representing 5,866,712 calving events from descendants were used. Several putative quantitative trait loci (QTL) regions associated with calving performance both within and across dairy and beef breeds were identified, although the majority were both breed- and trait-specific. QTL surrounding and encompassing the myostatin (MSTN) gene were associated (P < 5 × 10-8) with DCD and PM in both the CH and LM populations. The well-known Q204X mutation was the fifth strongest association with DCD in the CH population and accounted for 5.09% of the genetic variance in DCD. In contrast, none of the 259 segregating variants in MSTN were associated (P > × 10-6) with DCD in the LM population but a genomic region 617 kb downstream of MSTN was associated (P < 5 × 10-8). The genetic architecture for DCD differed in the HF population relative to the CH and LM, where two QTL encompassing ZNF613 on Bos taurus autosome (BTA)18 and PLAG1 on BTA14 were identified in the former. Pleiotropic SNP associated with all three calving performance traits were also identified in the three beef breeds; 5 SNP were pleiotropic in AA, 116 in LM, and 882 in CH but no SNP was associated with more than one trait within the HF population. The majority of these pleiotropic SNP were on BTA2 surrounding MSTN and were associated with both DCD and PM. Multiple previously reported, but also novel QTL, associated with calving performance were detected in this large study. These also included QTL regions harboring SNP with the same direction of allele substitution effect for both DCD and MCD thus contributing to a more effective simultaneous selection for both traits., (© The Author(s) 2020. Published by Oxford University Press on behalf of the American Society of Animal Science.)

Background: Linear type traits, which reflect the muscular characteristics of an animal, could provide insight into how, in some cases, morphologically very different animals can yield the same carcass weight. Such variability may contribute to differences in the overall value of the carcass since primal cuts vary greatly in price; such variability may also hinder successful genome-based association studies. Therefore, the objective of our study was to identify genomic regions that are associated with five muscularity linear type traits and to determine if these significant regions are common across five different breeds. Analyses were carried out using linear mixed models on imputed whole-genome sequence data in each of the five breeds, separately. Then, the results of the within-breed analyses were used to conduct an across-breed meta-analysis per trait., Results: We identified many quantitative trait loci (QTL) that are located across the whole genome and associated with each trait in each breed. The only commonality among the breeds and traits was a large-effect pleiotropic QTL on BTA2 that contained the MSTN gene, which was associated with all traits in the Charolais and Limousin breeds. Other plausible candidate genes were identified for muscularity traits including PDE1A, PPP1R1C and multiple collagen and HOXD genes. In addition, associated (gene ontology) GO terms and KEGG pathways tended to differ between breeds and between traits especially in the numerically smaller populations of Angus, Hereford, and Simmental breeds. Most of the SNPs that were associated with any of the traits were intergenic or intronic SNPs located within regulatory regions of the genome., Conclusions: The commonality between the Charolais and Limousin breeds indicates that the genetic architecture of the muscularity traits may be similar in these breeds due to their similar origins. Conversely, there were vast differences in the QTL associated with muscularity in Angus, Hereford, and Simmental. Knowledge of these differences in genetic architecture between breeds is useful to develop accurate genomic prediction equations that can operate effectively across breeds. Overall, the associated QTL differed according to trait, which suggests that breeding for a morphologically different (e.g. longer and wider versus shorter and smaller) more efficient animal may become possible in the future.

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., (Copyright © 2019 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

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.

Background: Quantitative genetic studies suggest the existence of variation at the genome level that affects the ability of cattle to resist to parasitic diseases. The objective of the current study was to identify regions of the bovine genome that are associated with resistance to endo-parasites., Methods: Individual cattle records were available for Fasciola hepatica-damaged liver from 18 abattoirs. Deregressed estimated breeding values (EBV) for F. hepatica-damaged liver were generated for genotyped animals with a record for F. hepatica-damaged liver and for genotyped sires with a least one progeny record for F. hepatica-damaged liver; 3702 animals were available. In addition, individual cow records for antibody response to F. hepatica on 6388 genotyped dairy cows, antibody response to Ostertagia ostertagi on 8334 genotyped dairy cows and antibody response to Neospora caninum on 4597 genotyped dairy cows were adjusted for non-genetic effects. Genotypes were imputed to whole-sequence; after edits, 14,190,141 single nucleotide polymorphisms (SNPs) and 16,603,644 SNPs were available for cattle with deregressed EBV for F. hepatica-damaged liver and cows with an antibody response to a parasitic disease, respectively. Association analyses were undertaken using linear regression on one SNP at a time, in which a genomic relationship matrix accounted for the relationships between animals., Results: Genomic regions for F. hepatica-damaged liver were located on Bos taurus autosomes (BTA) 1, 8, 11, 16, 17 and 18; each region included at least one SNP with a p value lower than 10 -6 . Five SNPs were identified as significant (q value < 0.05) for antibody response to N. caninum and were located on BTA21 or 25. For antibody response to F. hepatica and O. ostertagi, six and nine quantitative trait loci (QTL) regions that included at least one SNP with a p value lower than 10 -6 were identified, respectively. Gene set enrichment analysis revealed a significant association between functional annotations related to the olfactory system and QTL that were suggestively associated with endo-parasite phenotypes., Conclusions: A number of novel genomic regions were suggestively associated with endo-parasite phenotypes across the bovine genome and two genomic regions on BTA21 and 25 were associated with antibody response to N. caninum.

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; R 2 = 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., (Copyright © 2019 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2018 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2018 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

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.

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%.

Background: Calving difficulty and perinatal mortality are prevalent in modern-day cattle production systems. It is well-established that there is a genetic component to both traits, yet little is known about their underlying genomic architecture, particularly in beef breeds. Therefore, we performed a genome-wide association study using high-density genotypes to elucidate the genomic architecture of these traits and to identify regions of the bovine genome associated with them., Results: Genomic regions associated with calving difficulty (direct and maternal) and perinatal mortality were detected using two statistical approaches: (1) single-SNP (single nucleotide polymorphism) regression and (2) a Bayesian approach. Data included high-density genotypes on 770 Holstein-Friesian, 927 Charolais and 963 Limousin bulls. Several novel or previously identified genomic regions were detected but associations differed by breed. For example, two genomic associations, one each on chromosomes 18 and 2 explained 2.49 % and 3.13 % of the genetic variance in direct calving difficulty in the Holstein-Friesian and Charolais populations, respectively. Imputed Holstein-Friesian sequence data was used to refine the genomic regions responsible for significant associations. Several candidate genes on chromosome 18 were identified and four highly significant missense variants were detected within three of these genes (SIGLEC12, CTU1, and ZNF615). Nevertheless, only CTU1 contained a missense variant with a putative impact on direct calving difficulty based on SIFT (0.06) and Polyphen (0.95) scores. Using imputed sequence data, we refined a genomic region on chromosome 4 associated with maternal calving difficulty in the Holstein-Friesian population and found the strongest association with an intronic variant in the PCLO gene. A meta-analysis was performed across the three breeds for each calving performance trait to identify common variants associated with these traits in the three breeds. Our results suggest that a portion of the genetic variation in calving performance is common to all three breeds., Conclusion: The genomic architecture of calving performance is complex and mainly influenced by many polymorphisms of small effect. We identified several associations of moderate effect size but the majority were breed-specific, indicating that breed-specific alleles exist for calving performance or that the linkage phase between genotyped allele and causal mutation varies between breeds.

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.

Background: Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes that are present in an individual due to parents transmitting identical haplotypes to their offspring. The extent and frequency of ROHs may inform on the ancestry of an individual and its population. Here we use high density (n = 777,962) bi-allelic SNPs in a range of cattle breed samples to correlate ROH with the pedigree-based inbreeding coefficients and to validate subsequent analyses using 54,001 SNP genotypes. This study provides a first testing of the inference drawn from ROH through comparison with estimates of inbreeding from calculations based on the detailed pedigree data available for several breeds., Results: All animals genotyped on the HD panel displayed at least one ROH that was between 1-5 Mb in length with certain regions of the genome more likely to be involved in a ROH than others. Strong correlations (r = 0.75, p < 0.0001) existed between the pedigree-based inbreeding coefficient and a statistic based on sum of ROH of length > 0.5 KB and suggests that in the absence of an animal's pedigree data, the extent of a genome under ROH may be used to infer aspects of recent population history even from relatively few samples., Conclusions: Our findings suggest that ROH are frequent across all breeds but differing patterns of ROH length and burden illustrate variations in breed origins and recent management.