Genome-wide association study reveals candidate markers related to field fertility and semen quality traits in Holstein-Friesian bulls

In vitro assessment of bull semen quality is routinely used in bull semen processing centres in order to ensure that semen destined to be used in the field has passed minimum standards. Despite these stringent quality control checks, individual bulls that pass the quality control checks can still vary in field fertility by up to 25%. A genome-wide association study was undertaken to determine genetic markers associated with prefreeze and post-thaw bull sperm quality traits as well as field fertility. Genome-wide association analysis was performed using a single nucleotide polymorphism (SNP) regression mixed linear model in WOMBAT. Genes within a 250 Kb span of a suggestive (P ≤ 1 × 10^-5) SNP were considered as candidate genes. One SNP was associated with adjusted pregnancy rate, and 21 SNPs were associated across the seven semen quality traits (P ≤ 1 × 10^-5). Functional candidate genes include SIPA1L2 which was associated with adjusted pregnancy rate. This encodes a Rap GTPase-activating protein involved in Rap1 signalling pathway and was previously found to play a role in the process of sperm differentiation. Gene ontology (GO) analysis also identified significantly enriched biological processes involved protein tyrosine kinase activity including genes such as DYRK1A, TEC and TXK that were associated with sperm motility prior to freezing. Another candidate gene associated with post-thaw sperm motility was FHDC1 which coordinates actin filament and microtubule dynamics. The induced 11 GO terms in the ejaculates rejected after freezing trait were related to ATPase, phosphatase and hydrolase activity. These results reveal novel specific genomic regions and candidate genes associated with economically important phenotypes such as field fertility and semen quality traits.

Characterisation of eight cattle with Swyer syndrome by whole-genome sequencing

Swyer syndrome is where an individual has the karyotype of a typical male yet is phenotypically a female. The lack of a (functional) SRY gene located on the Y-chromosome is implicated in some cases of the Swyer syndrome, although many Swyer individuals with an apparently fully functional SRY gene have also been documented. The present study undertook whole genome sequence analyses of eight cattle with suspected Swyer syndrome and compared their genome to that of both a control male and female. Sequence analyses coupled with female phenotypes confirmed that all eight individuals had the 60,XY sex reversal Swyer syndrome. Seven of the eight Swyer syndrome individuals had a deletion on the Y chromosome encompassing the SRY gene (i.e., SRY-). The eighth individual had no obvious mutation in the SRY gene (SRY+) or indeed in any reported gene associated with sex reversal in mammals; a necropsy was performed on this individual. No testicles were detected during the necropsy. Histological examination of the reproductive tract revealed an immature uterine body and horns with inactive glandular tissue of normal histological appearance; both gonads were elongated, a characteristic of most reported cases of Swyer in mammals. The flanking sequence of 11 single nucleotide polymorphisms within 10 kb of the SRY gene are provided to help diagnose some cases of Swyer syndrome. These single nucleotide polymorphisms will not, however, detect all cases of Swyer syndrome since, as evidenced from the present study (and other studies), some individuals with the Swyer condition still contain the SRY gene (i.e., SRY+).

Risk factors for, and genetic association with, intestinal atresia in dairy calves

Intestinal atresia is an under-diagnosed congenital defect in cattle. It results in complete occlusion of the intestinal lumen and, unless surgically corrected, results in death or euthanasia of the affected calf. There is limited information on the incidence of this condition or on risk factors, including predisposing alleles, associated with the defect. In this study, active surveillance of 39 dairy farms over 8 years identified 197 cases of intestinal atresia among 56 454 calves born, an incidence of 0.35%. The majority of cases (83%) had occlusion of the jejunum, although cases with blockage of the colon (14%) or anus (4%) were also identified. The defect was twice as common in male as in female calves (p < 0.0001), and was more common in progeny of older cows than in progeny of first or second lactation cows (p < 0.001). Year and farm of birth were also significantly associated with incidence (p < 0.05). The incidence of intestinal atresia was highest among the progeny of three related Jersey sires, suggesting that a gene for intestinal atresia was segregating within this family. Linkage analysis utilising 28 affected progeny of two half-sib putative carrier sires identified two putative quantitative trait loci associated with the defect, on chromosomes 14 and 26, although no clear candidate genes were identified. There was no evidence of a sire-effect among the progeny of Holstein-Friesian sires. However, a case-control genome-wide association study involving 91 cases and 375 healthy controls, identified 31 SNP in 18 loci as associated with the defect in this breed. These data suggest that intestinal atresia in dairy calves is not a simple Mendelian trait as previously reported but a complex multigenic disorder.

The Association Between Genomic Heterozygosity and Carcass Merit in Cattle

The objective of the present study was to quantify the association between both pedigree and genome-based measures of global heterozygosity and carcass traits, and to identify single nucleotide polymorphisms (SNPs) exhibiting non-additive associations with these traits. The carcass traits of interest were carcass weight (CW), carcass conformation (CC) and carcass fat (CF). To define the genome-based measures of heterozygosity, and to quantify the non-additive associations between SNPs and the carcass traits, imputed, high-density genotype data, comprising of 619,158 SNPs, from 27,213 cattle were used. The correlations between the pedigree-based heterosis coefficient and the three defined genomic measures of heterozygosity ranged from 0.18 to 0.76. The associations between the different measures of heterozygosity and the carcass traits were biologically small, with positive associations for CW and CC, and negative associations for CF. Furthermore, even after accounting for the pedigree-based heterosis coefficient of an animal, part of the remaining variability in some of the carcass traits could be captured by a genomic heterozygosity measure. This signifies that the inclusion of both a heterosis coefficient based on pedigree information and a genome-based measure of heterozygosity could be beneficial to limiting bias in predicting additive genetic merit. Finally, one SNP located on Bos taurus (BTA) chromosome number 5 demonstrated a non-additive association with CW. Furthermore, 182 SNPs (180 SNPs on BTA 2 and two SNPs on BTA 21) demonstrated a non-additive association with CC, while 231 SNPs located on BTA 2, 5, 11, 13, 14, 18, 19 and 21 demonstrated a non-additive association with CF. Results demonstrate that heterozygosity both at a global level and at the level of individual loci contribute little to the variability in carcass merit.

Identification of genomic regions that exhibit sexual dimorphism for size and muscularity in cattle

Sexual dimorphism, the phenomenon whereby males and females of the same species are distinctive in some aspect of appearance or size, has previously been documented in cattle for traits such as growth rate and carcass merit using a quantitative genetics approach. No previous study in cattle has attempted to document sexual dimorphism at a genome level; therefore, the objective of the present study was to determine whether genomic regions associated with size and muscularity in cattle exhibited signs of sexual dimorphism. Analyses were undertaken on 10 linear-type traits that describe the muscular and skeletal characteristics of both males and females of five beef cattle breeds: 1,444 Angus (AA), 6,433 Charolais (CH), 1,129 Hereford, 8,745 Limousin (LM), and 1,698 Simmental. Genome-wide association analyses were undertaken using imputed whole-genome sequence data for each sex separately by breed. For each single-nucleotide polymorphism (SNP) that was segregating in both sexes, the difference between the allele substitution effect sizes for each sex, in each breed separately, was calculated. Suggestively (P ≤ 1 × 10-5) sexually dimorphic SNPs that were segregating in both males and females were detected for all traits in all breeds, although the location of these SNPs differed by both trait and breed. Significantly (P ≤ 1 × 10-8) dimorphic SNPs were detected in just three traits in the AA, seven traits in the CH, and three traits in the LM. The vast majority of all segregating autosomal SNPs (86% in AA to 94% in LM) had the same minor allele in both males and females. Differences (P ≤ 0.05) in allele frequencies between the sexes were observed for between 36% (LM) and 66% (AA) of the total autosomal SNPs that were segregating in both sexes. Dimorphic SNPs were located within a number of genes related to muscularity and/or size including the NAB1, COL5A2, and IWS1 genes on BTA2 that are located close to, and thought to be co-inherited with, the MSTN gene. Overall, sexual dimorphism exists in cattle at the genome level, but it is not consistent by either trait or breed.

Genome-wide association study of economically important traits in Charolais and Limousin beef cows

Genomic selection has proven effective for advancing genetic gain for key profit traits in dairy cattle production systems. However, its impact to-date on genetic improvement programs for beef cattle has been less effective. Despite this, the technology is thought to be particularly useful for low heritability traits such as those associated with reproductive efficiency. The objective of this study was to identify genetic variants associated with key determinants of reproductive and overall productive efficiency in beef cows. The analysis employed a large dataset derived from the national genetic evaluation program in Ireland for two of the most predominant beef breeds, viz. Charolais (n = 5 244 cows) and Limousin (n = 7 304 cows). Single nucleotide polymorphisms (SNPs) were identified as being statistically significantly associated (adj. P < 0.05) with both reproductive and productive traits for both breed types. However, there was little across breed commonality, with only two SNPs (rs110240246 and rs110344317; adj. P < 0.05) located within the genomic regions of the LCORL and MSTN genes respectively, identified in both Charolais and Limousin populations, associated with traits including carcass weight, cull-cow weight and live-weight. Significant SNPs within the MSTN gene were also associated with both reproduction and production related traits within each breed. Finally, traits including calving difficulty, calf mortality and calving interval were associated with SNPs within genomic regions comprising genes involved in cellular growth and lipid metabolism. Genetic variants identified as associated with both important reproductive efficiency and production related traits from this study warrant further analyses for their potential incorporation into breeding programmes to support the sustainability of beef cattle production.

A large interactive visual database of copy number variants discovered in taurine cattle

BACKGROUND

Copy number variants (CNVs) contribute to genetic diversity and phenotypic variation. We aimed to discover CNVs in taurine cattle using a large collection of whole-genome sequences and to provide an interactive database of the identified CNV regions (CNVRs) that includes visualizations of sequence read alignments, CNV boundaries, and genome annotations.

RESULTS

CNVs were identified in each of 4 whole-genome sequencing datasets, which together represent >500 bulls from 17 breeds, using a popular multi-sample read-depth-based algorithm, cn.MOPS. Quality control and CNVR construction, performed dataset-wise to avoid batch effects, resulted in 26,223 CNVRs covering 107.75 unique Mb (4.05%) of the bovine genome. Hierarchical clustering of samples by CNVR genotypes indicated clear separation by breeds. An interactive HTML database was created that allows data filtering options, provides graphical and tabular data summaries including Hardy-Weinberg equilibrium tests on genotype proportions, and displays genes and quantitative trait loci at each CNVR. Notably, the database provides sequence read alignments at each CNVR genotype and the boundaries of constituent CNVs in individual samples. Besides numerous novel discoveries, we corroborated the genotypes reported for a CNVR at the KIT locus known to be associated with the piebald coat colour phenotype in Hereford and some Simmental cattle.

CONCLUSIONS

We present a large comprehensive collection of taurine cattle CNVs in a novel interactive visual database that displays CNV boundaries, read depths, and genome features for individual CNVRs, thus providing users with a powerful means to explore and scrutinize CNVRs of interest more thoroughly.

Genomic Regions Associated With Skeletal Type Traits in Beef and Dairy Cattle Are Common to Regions Associated With Carcass Traits, Feed Intake and Calving Difficulty

Linear type traits describing the skeletal characteristics of an animal are moderately to strongly genetically correlated with a range of other performance traits in cattle including feed intake, reproduction traits and carcass merit; thus, type traits could also provide useful insights into the morphological differences among animals underpinning phenotypic differences in these complex traits. The objective of the present study was to identify genomic regions associated with five subjectively scored skeletal linear traits, to determine if these associated regions are common in multiple beef and dairy breeds, and also to determine if these regions overlap with those proposed elsewhere to be associated with correlated performance traits. Analyses were carried out using linear mixed models on imputed whole genome sequence data separately in 1,444 Angus, 1,129 Hereford, 6,433 Charolais, 8,745 Limousin, 1,698 Simmental, and 4,494 Holstein-Friesian cattle, all scored for the linear type traits. There was, on average, 18 months difference in age at assessment of the beef versus the dairy animals. While the majority of the identified quantitative trait loci (QTL), and thus genes, were both trait-specific and breed-specific, a large-effect pleiotropic QTL on BTA6 containing the NCAPG and LCORL genes was associated with all skeletal traits in the Limousin population and with wither height in the Angus. Other than that, little overlap existed in detected QTLs for the skeletal type traits in the other breeds. Only two QTLs overlapped the beef and dairy breeds; both QTLs were located on BTA5 and were associated with height in both the Angus and the Holstein-Friesian, despite the difference in age at assessment. Several detected QTLs in the present study overlapped with QTLs documented elsewhere that are associated with carcass traits, feed intake, and calving difficulty. While most breeding programs select for the macro-traits like carcass weight, carcass conformation, and feed intake, the higher degree of granularity with selection on the individual linear type traits in a multi-trait index underpinning the macro-level goal traits, presents an opportunity to help resolve genetic antagonisms among morphological traits in the pursuit of the animal with optimum performance metrics.

Genomic regions associated with muscularity in beef cattle differ in five contrasting cattle breeds

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.

Genomic Regions Associated With Gestation Length Detected Using Whole-Genome Sequence Data Differ Between Dairy and Beef Cattle

While many association studies exist that have attempted to relate genomic markers to phenotypic performance in cattle, very few have considered gestation length as a phenotype, and of those that did, none used whole genome sequence data from multiple breeds. The objective of the present study was therefore to relate imputed whole genome sequence data to estimated breeding values for gestation length using 22,566 sires (representing 2,262,706 progeny) of multiple breeds [Angus (AA), Charolais (CH), Holstein-Friesian (HF), and Limousin (LM)]. The associations were undertaken within breed using linear mixed models that accounted for genomic relatedness among sires; a separate association analysis was undertaken with all breeds analysed together but with breed included as a fixed effect in the model. Furthermore, the genome was divided into 500 kb segments and whether or not segments harboured a single nucleotide polymorphism (SNP) with a P ≤ 1 × 10^-4 common to different combinations of breeds was determined. Putative quantitative trait loci (QTL) regions associated with gestation length were detected in all breeds; significant associations with gestation length were only detected in the HF population and in the across-breed analysis of all 22,566 sires. Twenty-five SNPs were significantly associated (P ≤ 5 × 10^-8) with gestation length in the HF population. Of the 25 significant SNPs, 18 were located within three QTLs on Bos taurus autosome number (BTA) 18, six were in two QTL on BTA19, and one was located within a QTL on BTA7. The strongest association was rs381577268, a downstream variant of ZNF613 located within a QTL spanning from 58.06 to 58.19 Mb on BTA18; it accounted for 1.37% of the genetic variance in gestation length. Overall there were 11 HF animals within the edited dataset that were homozygous for the T allele at rs381577268 and these had a 3.3 day longer (P < 0.0001) estimated breeding value (EBV) for gestation length than the heterozygous animals and a 4.7 day longer (P < 0.0001) EBV for gestation length than the homozygous CC animals. The majority of the 500 kb windows harboring a SNP with a P ≤ 1 × 10^-4 were unique to a single breed and no window was shared among all four breeds for gestation length, suggesting any QTLs identified are breed-specific associations.

A mating advice system in dairy cattle incorporating genomic information

This study investigated the effects of alternative mating programs that incorporate genomic information on expected progeny herd performance and inbreeding, as well as methods to include un-genotyped animals in such mating programs. A total of 54,535 Holstein-Friesian cattle with imputed high-density genotypes (547,650 SNP after edits) were available. First, to quantify the accuracy of imputing un-genotyped animals (often an issue in populations), a sub-population of 729 genotyped animals had their genotypes masked, and their allele dosages were imputed, using linear regression exploiting information on genotyped relatives. The reference population for imputation included all genotyped animals, excluding the 729 selected animals and their sires, dams, and grandsires, and had either (1) their sires' genotypes, (2) their dams' genotypes (3) both their sires' and their dams' genotypes, or (4) both their sires' and maternal grandsires' genotypes introduced into the reference population. The correlations between true genotypes and the imputed allele dosages ranged from 0.58 (sire only) to 0.68 (both sire and dam). A herd of 100 cows was then simulated (1,000 replicates) from the sub-population of 729 imputed animals. The top 10 bulls from the genotyped population, based on their total genetic merit index (TMI) were selected to be used as sires. Three mating allotment methods were investigated: (1) random mating, (2) sequential mating based on maximizing only the expected TMI of the progeny, and (3) linear programming to maximize a generated index constructed to maximize genetic merit and minimize expected progeny inbreeding as well as intra- and inter-progeny variability in genetic merit. Relationships among candidate parents were calculated using either the pedigree relationship matrix or the genomic relationship matrix; the latter was constructed using either the true genotypes of both parents or the true genotypes of the sire plus the imputed allele dosages of the dam. Using the genomic co-ancestry estimates resulted in lower average herd expected genomic inbreeding levels compared with using the pedigree-based co-ancestry estimates. Additionally, if the dams were not genotyped, using their imputed allele dosages also resulted in lower average herd expected inbreeding levels compared with using the pedigree co-ancestry estimates. The inter-progeny coefficient of variation for selected traits, milk and fertility, estimated breeding values were reduced by 12 to 65% using the linear programing method compared with sequential mating.

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.

Characterization of copy number variants in a large multibreed population of beef and dairy cattle using high-density single nucleotide polymorphism genotype data

Copy number variants (CNVs) are a form of genomic variation that changes the structure of the genome through deletion or duplication of stretches of DNA. The objective of the present study was to characterize CNVs in a large multibreed population of beef and dairy bulls. The CNVs were called on the autosomes of 5,551 cattle from 22 different beef and dairy breeds, using 2 freely available software suites, QuantiSNP and PennCNV. All CNVs were classified into either deletions or duplications. The median concordance between PennCNV and QuantiSNP, per animal, was 18.5% for deletions and 0% for duplications. The low concordance rate between PennCNV and QuantiSNP indicated that neither algorithm, by itself, could identify all CNVs in the population. In total, PennCNV and QuantiSNP collectively identified 747,129 deletions and 432,523 duplications; 80.2% of all duplications and 69.1% of all deletions were present only once in the population. Only 0.154% of all CNVs identified were present in more than 50 animals in the population. The distribution of the percentage of the autosomes that were composed of deletions, per animal, was positively skewed, as was the distribution for the percentage of the autosomes that were composed of duplications, per animal. The first quartile, median, and third quartile of the distribution of the percentage of the autosomes that were composed of deletions were 0.019%, 0.037%, and 0.201%, respectively. The first quartile, median, and third quartile of the distribution of the percentage of the autosomes that were composed of duplications were 0.013%, 0.028%, and 0.076%, respectively. The distributions of the number of deletions and duplications per animal were both positively skewed. The interquartile range for the number of deletions per animal in the population was between 16 and 117, whereas for duplications it was between 8 and 23. Per animal, there tended to be twice as many deletions as duplications. The distribution of the length of deletions was positively skewed, as was the distribution of the length of duplications. The interquartile range for the length of deletions in the population was between 25 and 101 kb, and for duplications the interquartile range was between 46 and 235 kb. Per animal, duplications tended to be twice as long as deletions. This study provides a description of the characteristics and distribution of CNVs in a large multibreed population of beef and dairy cattle.

Genetic covariance components within and among linear type traits differ among contrasting beef cattle breeds

Linear type traits describing the skeletal, muscular, and functional characteristics of an animal are routinely scored on live animals in both the dairy and beef cattle industries. Previous studies have demonstrated that genetic parameters for certain performance traits may differ between breeds; no study, however, has attempted to determine if differences exist in genetic parameters of linear type traits among breeds or sexes. Therefore, the objective of the present study was to determine if genetic covariance components for linear type traits differed among five contrasting cattle breeds, and to also investigate if these components differed by sex. A total of 18 linear type traits scored on 3,356 Angus (AA), 31,049 Charolais (CH), 3,004 Hereford (HE), 35,159 Limousin (LM), and 8,632 Simmental (SI) were used in the analysis. Data were analyzed using animal linear mixed models which included the fixed effects of sex of the animal (except in the investigation into the presence of sexual dimorphism), age at scoring, parity of the dam, and contemporary group of herd-date of scoring. Differences (P < 0.05) in heritability estimates, between at least two breeds, existed for 13 out of 18 linear type traits. Differences (P < 0.05) also existed between the pairwise within-breed genetic correlations among the linear type traits. Overall, the linear type traits in the continental breeds (i.e., CH, LM, SI) tended to have similar heritability estimates to each other as well as similar genetic correlations among the same pairwise traits, as did the traits in the British breeds (i.e., AA, HE). The correlation between a linear function of breeding values computed conditional on covariance parameters estimated from the CH breed with a linear function of breeding values computed conditional on covariance parameters estimated from the other breeds was estimated. Replacing the genetic covariance components estimated in the CH breed with those of the LM had least effect but the impact was considerable when the genetic covariance components of the AA were used. Genetic correlations between the same linear type traits in the two sexes were all close to unity (≥0.90) suggesting little advantage in considering these as separate traits for males and females. Results for the present study indicate the potential increase in accuracy of estimated breeding value prediction from considering, at least, the British breed traits separate to continental breed traits.

Genetic parameters of ovarian and uterine reproductive traits in dairy cows

The objective of the study was to estimate genetic parameters of detailed reproductive traits derived from ultrasound examination of the reproductive tract as well as their genetic correlations with traditional reproductive traits. A total of 226,141 calving and insemination records as well as 74,134 ultrasound records from Irish dairy cows were used. Traditional reproductive traits included postpartum interval to first service, conception, and next calving, as well as the interval from first to last service; number of inseminations, pregnancy rate to first service, pregnant within 42 d of the herd breeding season, and submission in the first 21 d of the herd breeding season were also available. Detailed reproductive traits included resumed cyclicity at the time of ultrasound examination, incidence of multiple ovulations, incidence of early postpartum ovulation, heat detection, ovarian cystic structures, embryo loss, and uterine score; the latter was a subjectively assessed on a scale of 1 (little fluid with normal uterine tone) to 4 (large quantity of fluid with a flaccid uterine tone). Variance (and covariance) components were estimated using repeatability animal linear mixed models. Heritability for all reproductive traits were generally low (0.001-0.05), with the exception of traits related to cyclicity postpartum, regardless if defined traditionally (0.07; calving to first service) or from ultrasound examination [resumed cyclicity at the time of examination (0.07) or early postpartum ovulation (0.10)]. The genetic correlations among the detailed reproductive traits were generally favorable. The exception was the genetic correlation (0.29) between resumed cyclicity and uterine score; superior genetic merit for cyclicity postpartum was associated with inferior uterine score. Superior genetic merit for most traditional reproductive traits was associated with superior genetic merit for resumed cyclicity (genetic correlations ranged from -0.59 to -0.36 and from 0.56 to 0.70) and uterine score (genetic correlations ranged from -0.47 to 0.32 and from 0.25 to 0.52). Genetic predisposition to an increased incidence of embryo loss was associated with both an inferior uterine score (0.24) and inferior genetic merit for traditional reproductive traits (genetic correlations ranged from -0.52 to -0.42 and from 0.33 to 0.80). The results from the present study indicate that selection based on traditional reproductive traits, such as calving interval or days open, resulted in improved genetic merit of all the detailed reproductive traits evaluated in this study. Additionally, greater accuracy of selection for calving interval is expected for a relatively small progeny group size when detailed reproductive traits are included in a multitrait genetic evaluation.