Comparison of different validation methods for single-step genomic evaluations based on a simulated cattle population

The validation of estimated breeding values from single-step genomic BLUP (ssGBLUP) is an important topic, as more and more countries and animal populations are currently changing their genomic prediction to single-step. The objective of this work was to compare different methods to validate single-step genomic breeding values (GEBV). The investigations were carried out using a simulation study based on the German-Austrian-Czech Fleckvieh population. To test the validation methods under different conditions, several biased and unbiased scenarios were simulated. The application of the widely used Interbull GEBV test to the single-step method is only possible to a limited extent, partly because of genomic preselection, which biases conventional estimated breeding values. Alternative validation methods considered in the study are the linear regression method proposed by Legarra and Reverter, the improved genomic validation including additional regressions as suggested by VanRaden and an adaptation of the Interbull GEBV test using daughter yield deviations (DYD) from ssGBLUP instead of pedigree BLUP. The comparison of the different methods for the different scenarios showed that for males the methods based on GEBV estimate the dispersion more accurate and less biased compared with the GEBV test using DYD from ssGBLUP, whereas the standard Interbull GEBV test is highly affected by genomic preselection for males. For females, the GEBV test using yield deviations from ssGBLUP results in better estimations for the true dispersion.

The Combined Use of Automated Milking System and Sensor Data to Improve Detection of Mild Lameness in Dairy Cattle

This study aimed to develop a tool to detect mildly lame cows by combining already existing data from sensors, AMSs, and routinely recorded animal and farm data. For this purpose, ten dairy farms were visited every 30–42 days from January 2020 to May 2021. Locomotion scores (LCS, from one for nonlame to five for severely lame) and body condition scores (BCS) were assessed at each visit, resulting in a total of 594 recorded animals. A questionnaire about farm management and husbandry was completed for the inclusion of potential risk factors. A lameness incidence risk (LCS ≥ 2) was calculated and varied widely between farms with a range from 27.07 to 65.52%. Moreover, the impact of lameness on the derived sensor parameters was inspected and showed no significant impact of lameness on total rumination time. Behavioral patterns for eating, low activity, and medium activity differed significantly in lame cows compared to nonlame cows. Finally, random forest models for lameness detection were fit by including different combinations of influencing variables. The results of these models were compared according to accuracy, sensitivity, and specificity. The best performing model achieved an accuracy of 0.75 with a sensitivity of 0.72 and specificity of 0.78. These approaches with routinely available data and sensor data can deliver promising results for early lameness detection in dairy cattle. While experimental automated lameness detection systems have achieved improved predictive results, the benefit of this presented approach is that it uses results from existing, routinely recorded, and therefore widely available data.

A 1-bp deletion in bovine QRICH2 causes low sperm count and immotile sperm with multiple morphological abnormalities

Background

Semen quality and insemination success are monitored in artificial insemination bulls to ensure high male fertility rates. Only ejaculates that fulfill minimum quality requirements are processed and eventually used for artificial inseminations. We examined 70,990 ejaculates from 1343 Brown Swiss bulls to identify bulls from which all ejaculates were rejected due to low semen quality. This procedure identified a bull that produced 12 ejaculates with an aberrantly small number of sperm (0.2 ± 0.2 × 109 sperm per mL) which were mostly immotile due to multiple morphological abnormalities.

Results

The genome of this bull was sequenced at a 12× coverage to investigate a possible genetic cause. Comparing the sequence variant genotypes of this bull with those from 397 fertile bulls revealed a 1-bp deletion in the coding sequence of the QRICH2 gene which encodes the glutamine rich 2 protein, as a compelling candidate causal variant. This 1-bp deletion causes a frameshift in translation and a premature termination codon (ENSBTAP00000018337.1:p.Cys1644AlafsTer52). The analysis of testis transcriptomes from 76 bulls showed that the transcript with the premature termination codon is subject to nonsense-mediated mRNA decay. The 1-bp deletion resides in a 675-kb haplotype that includes 181 single nucleotide polymorphisms (SNPs) from the Illumina BovineHD Bead chip. This haplotype segregates at a frequency of 5% in the Brown Swiss cattle population. Our analysis also identified another bull that carried the 1-bp deletion in the homozygous state. Semen analyses from the second bull confirmed low sperm concentration and immotile sperm with multiple morphological abnormalities that primarily affect the sperm flagellum and, to a lesser extent, the sperm head.

Conclusions

A recessive loss-of-function allele of the bovine QRICH2 gene likely causes low sperm concentration and immotile sperm with multiple morphological abnormalities. Routine sperm analyses unambiguously identify homozygous bulls for this allele. A direct gene test can be implemented to monitor the frequency of the undesired allele in cattle populations.

Exploring the Potential Genetic Heterogeneity in the Incidence of Hoof Disorders in Austrian Fleckvieh and Braunvieh Cattle

Genetic heterogeneity denotes the situation when different genetic architectures underlying diverse populations result in the same phenotype. In this study, we explore the genetic background underlying differences in the incidence of hoof disorders between Braunvieh and Fleckvieh cattle in the context of genetic heterogeneity between the breeds. Despite potentially higher power of testing due to twice as large sample size, none of the SNPs was significantly associated with the total number of hoof disorders in Fleckvieh, while 15 SNPs were significant in Braunvieh. The most promising candidate genes in Braunvieh were as follows: CBLB on BTA1, which causes arthritis in rats; CAV2 on BTA4, which affects skeletal muscles in mice; PTHLH on BTA5, which causes disease phenotypes related to the skeleton in humans, mice, and zebrafish; and SORCS2 on BTA6, which causes decreased susceptibility to injury in mice. Some of the significant SNPs (BTA1, BTA4, BTA5, BTA13, and BTA16) revealed allelic heterogeneity-i.e., different allele frequencies between Fleckvieh and Braunvieh. Some of the significant regions (BTA1, BTA5, BTA13, and BTA16) correlated to inter-breed differences in linkage disequilibrium (LD) structure and may thus represent false-positive heterogeneity. However, positions on BTA6 (SORCS2), BTA14, and BTA24 mark Braunvieh-specific regions. We hypothesize that the observed genetic heterogeneity of hoof disorders is a by-product of different selection goals defined for the analyzed breeds-toward dairy production in Braunvieh and toward beef production in Fleckvieh. Based on the current dataset, it is not possible to unequivocally confirm or exclude the hypothesis of genetic heterogeneity in the susceptibility to hoof disorders between Fleckvieh and Braunvieh. The main reason for the problem is that the potential heterogeneity was explored through SNP-phenotype associations and not through causal mutations, due to a limited SNP density offered by the SNP-chip. The rationale against genetic heterogeneity comprises a limited power of detection of true associations as well as differences in the length of LD blocks and in linkage phase between breeds. On the other hand, different selection goals defined for the analyzed breeds accompanied by no systematic, genome-wide differences in LD structure between the breeds favor the heterogeneity hypothesis at some smaller genomic regions.

Activation of cryptic splicing in bovine WDR19 is associated with reduced semen quality and male fertility

Cattle are ideally suited to investigate the genetics of male reproduction, because semen quality and fertility are recorded for all ejaculates of artificial insemination bulls. We analysed 26,090 ejaculates of 794 Brown Swiss bulls to assess ejaculate volume, sperm concentration, sperm motility, sperm head and tail anomalies and insemination success. The heritability of the six semen traits was between 0 and 0.26. Genome-wide association testing on 607,511 SNPs revealed a QTL on bovine chromosome 6 that was associated with sperm motility (P = 2.5 x 10⁻²⁷), head (P = 2.0 x 10⁻⁴⁴) and tail anomalies (P = 7.2 x 10⁻⁴⁹) and insemination success (P = 9.9 x 10⁻¹³). The QTL harbors a recessive allele that compromises semen quality and male fertility. We replicated the effect of the QTL on fertility (P = 7.1 x 10⁻³²) in an independent cohort of 2481 Brown Swiss bulls. The analysis of whole-genome sequencing data revealed that a synonymous variant (BTA6:58373887C>T, rs474302732) in WDR19 encoding WD repeat-containing protein 19 was in linkage disequilibrium with the fertility-associated haplotype. WD repeat-containing protein 19 is a constituent of the intraflagellar transport complex that is essential for the physiological function of motile cilia and flagella. Bioinformatic and transcription analyses revealed that the BTA6:58373887 T-allele activates a cryptic exonic splice site that eliminates three evolutionarily conserved amino acids from WDR19. Western blot analysis demonstrated that the BTA6:58373887 T-allele decreases protein expression. We make the remarkable observation that, in spite of negative effects on semen quality and bull fertility, the BTA6:58373887 T-allele has a frequency of 24% in the Brown Swiss population. Our findings are the first to uncover a variant that is associated with quantitative variation in semen quality and male fertility in cattle.

In cattle farming, artificial insemination is the most common method of breeding. To ensure high fertilization rates, ejaculate quality and insemination success are closely monitored in artificial insemination bulls. We analyse semen quality, insemination success and microarray-called genotypes at more than 600,000 genome-wide SNP markers of 794 bulls to identify a recessive allele that compromises semen quality. We take advantage of whole-genome sequencing to pinpoint a variant in the coding sequence of WDR19 encoding WD repeat-containing protein 19 that activates a novel exonic splice site. Our results indicate that cryptic splicing in WDR19 is associated with reduced male reproductive performance. This is the first report of a variant that contributes to quantitative variation in bovine semen quality.

Two-stage genome-wide association study for the identification of causal variants underlying hoof disorders in cattle

Feet and legs disorders influence dairy cattle breeding by their effect on animal welfare, economic losses due to lower production and fertility, costs of treatment, and problems with herd management. In our study, we estimated heritabilities and performed a 2-step GWAS for 3 traits describing hoof health: hoof health status defined by a veterinarian (HSV), hoof health status defined by a claw trimmer (HSC), and the total number of hoof disorders (NHD), scored in 1,998 Fleckvieh and 979 Braunvieh cows. The individuals were genotyped with a high-density (HD) panel consisting of 76,934 SNP. For significant genomic regions, the SNP information was enhanced by SNP imputed from the whole-genome sequence of Fleckvieh and Braunvieh bulls from the 1000 Bulls Genome project. The heritabilities were estimated to be 0.035 for HSV, 0.249 for HSC, and 0.279 for NHD. Based on the first-stage GWAS with SNP from the HD panel, 7 significant genomic regions on 6 chromosomes were defined: (1) 120 SNP spanning 15,522 bp on BTA1, including the TOPBP1 gene; (2) 4,139 SNP spanning 1,426,046 bp on BTA7, including the RIOK2 and RGMB genes; (3) 167 SNP spanning 167,352 bp on BTA13, including the C13H20orf194 gene; (4) 2 regions on BTA14, one harboring 1,071 SNP spanning 380,024 bp, including RRM2B and NCALD, and the other comprising 632 SNP spanning 385,111 bp, including STK3; (5) 328 SNP on BTA15, spanning 235,567 bp between FAM168A and PLEKHB1; and (6) 1,549 SNP on BTA22, spanning 596,101 bp in the neighborhood of PTPRG. Then, we conducted a second-stage GWAS based on SNP from whole-genome sequences within the significant regions obtained in the first stage of the analysis. For HSV, the highest additive effect was estimated for 23 SNP located within a region on BTA15, close to FAM168A, corresponding to a predicted gene sequence. For HSC, the highest additive effect was attributed to 44 SNP located within a region of BTA22 corresponding to 4 predicted gene sequences, with rs135082893 within a sequence encoding a microRNA. Another potential causal mutation for HSC was rs134142607 on BTA13, within the exon of C13H20orf194. For NHD, 33 SNP with the highest estimated effect were located on BTA7 within a region of a predicted gene positioned between RIOK2 and RGMB. On BTA14, all significant SNP were located in introns of STK3, which is responsible for the "abnormal gait" phenotype in mice.

[Congenital flexural deformity in 93 calves - appearance, treatment techniques and results of pedigree analysis]

INTRODUCTION

Contracture of the flexor tendons (CFT) is very common in calves and it is usually diagnosed within the first few days after birth (congenital flexural deformity). However, CFT can appear even in older calves caused by chronic pain. The aetiology of CFT is still unknown. In this study, the distribution of sex, age, breed, the severity of flexural deformity, concurrent presence of other diseases, applied treatment methods for flexural deformity, and the outcome of calves with CFT, which were examined at the University Clinic for Ruminants in Vienna from 2001 to 2016, were evaluated retrospectively. 93 calves were admitted with CFT in the observation period. 70 (75.3%) calves were male and 78 (83.9%) of the affected animals were Simmental calves. The age of calves with CFT varied from one day to 41 days. Twenty-six calves suffered exclusively from CFT, and CFT was diagnosed as an additional finding in 67 calves. 91 animals (97.8%) showed CFT on the front limbs, 79 of them (84.9%) on both front limbs. The distribution of the severity scores was as follows: 69 calves (74.2%) had score 1, 17 calves (18.3%) had score 2, three calves (3.2%) had score 3. Three additional calves (3.2%) had a score 1 CFT on one front limb and a score 2 CFT on the other front limb, and one additional calf showed all three scores on both front limbs and one hind limb. 69 patients (74.2%) could be discharged with a significant improvement in CFT after treatment and 24 calves (25.8%) had to be euthanized due to other severe diseases. The results of the applied pedigree analysis do not show that a single gene mutation is the cause for the development of CFT, but rather a complex hereditary pattern has to be assumed. Depending on the severity of CFT and the presence of other concurrent diseases, an early and consistent therapy has to be carried out to achieve the highest possible success. Since animals with CFT are usually restricted in their movement, sufficient colostrum intake must be ensured within the first hours of life.

On the genomic regions associated with milk lactose in Fleckvieh cattle

Lactose is a sugar uniquely found in mammals' milk and it is the major milk solid in bovines. Lactose yield (LY, kg/d) is responsible for milk volume, whereas lactose percentage (LP) is thought to be more related to epithelial integrity and thus to udder health. There is a paucity of studies that have investigated lactose at the genomic level in dairy cows. This paper aimed to improve our knowledge on LP and LY, providing new insights into the significant genomic regions affecting these traits. A genome-wide association study for LP and LY was carried out in Fleckvieh cattle by using bulls' deregressed estimated breeding values of first lactation as pseudo-phenotypes. Heritabilities of first-lactation test-day LP and LY estimated using linear animal models were 0.38 and 0.25, respectively. A total of 2,854 bulls genotyped with a 54K SNP chip were available for the genome-wide association study; a linear mixed model approach was adopted for the analysis. The significant SNP of LP were scattered across the whole genome, with signals on chromosomes 1, 2, 3, 7, 12, 16, 18, 19, 20, 28, and 29; the top 4 significant SNP explained 4.90% of the LP genetic variance. The signals were mostly in regions or genes with involvement in molecular intra- or extracellular transport; for example, CDH5, RASGEF1C, ABCA6, and SLC35F3. A significant region within chromosome 20 was previously shown to affect mastitis or somatic cell score in cattle. As regards LY, the significant SNP were concentrated in fewer regions (chromosomes 6 and 14), related to mastitis/somatic cell score, immune response, and transport mechanisms. The 5 most significant SNP for LY explained 8.45% of genetic variance and more than one-quarter of this value has to be attributed to the variant within ADGRB1. Significant peaks in target regions remained even after adjustment for the 2 most significant variants previously detected on BTA6 and BTA14. The present study is a prelude for deeper investigations into the biological role of lactose for milk secretion and volume determination, stressing the connection with genes regulating intra- or extracellular trafficking and immune and inflammatory responses in dairy cows. Also, these results improve the knowledge on the relationship between lactose and udder health; they support the idea that LP and its derived traits are potential candidates as indicators of udder health in breeding programs aimed to enhance cows' resistance to mastitis.

The effect of mislabeled phenotypic status on the identification of mutation-carriers from SNP genotypes in dairy cattle

Background

Statistical and machine learning applications are increasingly popular in animal breeding and genetics, especially to compute genomic predictions for phenotypes of interest. Noise (errors) in the data may have a negative impact on the accuracy of predictions. The effects of noisy data have been investigated in genome-wide association studies for case–control experiments, and in genomic predictions for binary traits in plants. No studies have been published yet on the impact of noisy data in animal genomics. In this work, the susceptibility to noise of five classification models (Lasso-penalised logistic regression—Lasso, K-nearest neighbours—KNN, random forest—RF, support vector machines with linear—SVML—or radial—SVMR—kernel) was tested. As illustration, the identification of carriers of a recessive mutation in cattle (Bos taurus) was used. A population of 3116 Fleckvieh animals with SNP genotypes on the same chromosome as the mutation locus (BTA 19) was available. The carrier status (0/1 phenotype) was randomly sampled to generate noise. Increasing proportions of noise—up to 20%— were introduced in the data.

Results

SVMR and Lasso were relatively more robust to noise in the data, with total accuracy still above 0.975 and TPR (true positive rate; accuracy in the minority class) in the range 0.5–0.80 also with 17.5–20% mislabeled observations. The performance of SVML and RF decreased monotonically with increasing noise in the data, while KNN constantly failed to identify mutation carriers (observations in the minority class). The computation time increased with noise in the data, especially for the two support vector machines classifiers.

Conclusions

This work was the first to assess the impact of phenotyping errors on the accuracy of genomic predictions in animal genetics. The choice of the classification method can influence results in terms of higher or lower susceptibility to noise. In the presented problem, SVM with radial kernel performed relatively well even when the proportion of errors in the data reached 12.5%. Lasso was the second best method, while SVML, RF and KNN were very sensitive to noise. Taking into account both accuracy and computation time, Lasso provided the best combination.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-017-2540-x) contains supplementary material, which is available to authorized users.

Use of SNP genotypes to identify carriers of harmful recessive mutations in cattle populations

BACKGROUND

SNP (single nucleotide polymorphisms) genotype data are increasingly available in cattle populations and, among other things, can be used to predict carriers of specific mutations. It is therefore convenient to have a practical statistical method for the accurate classification of individuals into carriers and non-carriers. In this paper, we compared - through cross-validation- five classification models (Lasso-penalized logistic regression -Lasso, Support Vector Machines with either linear or radial kernel -SVML and SVMR, k-nearest neighbors -KNN, and multi-allelic gene prediction -MAG), for the identification of carriers of the TUBD1 recessive mutation on BTA19 (Bos taurus autosome 19), known to be associated with high calf mortality. A population of 3116 Fleckvieh and 392 Brown Swiss animals genotyped with the 54K SNP-chip was available for the analysis.

RESULTS

In general, the use of SNP genotypes proved to be very effective for the identification of mutation carriers. The best predictive models were Lasso, SVML and MAG, with an average error rate, respectively, of 0.2 %, 0.4 % and 0.6 % in Fleckvieh, and 1.2 %, 0.9 % and 1.7 % in Brown Swiss. For the three models, the false positive rate was, respectively, 0.1 %, 0.1 % and 0.2 % in Fleckvieh, and 3.0 %, 2.4 % and 1.6 % in Brown Swiss; the false negative rate was 4.4 %, 7.6 %1.0 % in Fleckvieh, and 0.0 %, 0.1% and 0.8 % in Brown Swiss. MAG appeared to be more robust to sample size reduction: with 25 % of the data, the average error rate was 0.7 % and 2.2 % in Fleckvieh and Brown Swiss, compared to 2.1 % and 5.5 % with Lasso, and 2.6 % and 12.0 % with SVML.

CONCLUSIONS

The use of SNP genotypes is a very effective and efficient technique for the identification of mutation carriers in cattle populations. Very few misclassifications were observed, overall and both in the carriers and non-carriers classes. This indicates that this is a very reliable approach for potential applications in cattle breeding.

A practical approach to detect ancestral haplotypes in livestock populations

Background

The effects of different evolutionary forces are expected to lead to the conservation, over many generations, of particular genomic regions (haplotypes) due to the development of linkage disequilibrium (LD). The detection and identification of early (ancestral) haplotypes can be used to clarify the evolutionary dynamics of different populations as well as identify selection signatures and genomic regions of interest to be used both in conservation and breeding programs. The aims of this study were to develop a simple procedure to identify ancestral haplotypes segregating across several generations both within and between populations with genetic links based on whole-genome scanning. This procedure was tested with simulated and then applied to real data from different genotyped populations of Spanish, Fleckvieh, Simmental and Brown-Swiss cattle.

Results

The identification of ancestral haplotypes has shown coincident patterns of selection across different breeds, allowing the detection of common regions of interest on different bovine chromosomes and mirroring the evolutionary dynamics of the studied populations. These regions, mainly located on chromosomes BTA5, BTA6, BTA7 and BTA21 are related with certain animal traits such as coat colour and milk protein and fat content.

Conclusion

In agreement with previous studies, the detection of ancestral haplotypes provides useful information for the development and comparison of breeding and conservation programs both through the identification of selection signatures and other regions of interest, and as indicator of the general genetic status of the populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0405-2) contains supplementary material, which is available to authorized users.

A missense mutation in TUBD1 is associated with high juvenile mortality in Braunvieh and Fleckvieh cattle

Background

Haplotypes with reduced or missing homozygosity may harbor deleterious alleles that compromise juvenile survival. A scan for homozygous haplotype deficiency revealed a short segment on bovine chromosome 19 (Braunvieh haplotype 2, BH2) that was associated with high juvenile mortality in Braunvieh cattle. However, the molecular genetic underpinnings and the pathophysiology of BH2 remain to be elucidated.

Results

The frequency of BH2 was 6.5 % in 8,446 Braunvieh animals from the national bovine genome databases. Both perinatal and juvenile mortality of BH2 homozygous calves were higher than the average in Braunvieh cattle resulting in a depletion of BH2 homozygous adult animals (P = 9.3x10⁻¹²). The analysis of whole-genome sequence data from 54 Braunvieh animals uncovered a missense mutation in TUBD1 (rs383232842, p.H210R) that was compatible with recessive inheritance of BH2. The availability of sequence data of 236 animals from diverse bovine populations revealed that the missense mutation also segregated at a low frequency (1.7 %) in the Fleckvieh breed. A validation study in 37,314 Fleckvieh animals confirmed high juvenile mortality of homozygous calves (P = 2.2x10⁻¹⁵). Our findings show that the putative disease allele is located on an ancestral haplotype that segregates in Braunvieh and Fleckvieh cattle. To unravel the pathophysiology of BH2, six homozygous animals were examined at the animal clinic. Clinical and pathological findings revealed that homozygous calves suffered from chronic airway disease possibly resulting from defective cilia in the respiratory tract.

Conclusions

A missense mutation in TUBD1 is associated with high perinatal and juvenile mortality in Braunvieh and Fleckvieh cattle. The mutation is located on a common haplotype likely originating from an ancient ancestor of Braunvieh and Fleckvieh cattle. Our findings demonstrate for the first time that deleterious alleles may segregate across closed cattle breeds without recent admixture. Homozygous calves suffer from chronic airway disease resulting in poor growth performance and high juvenile mortality. The respiratory manifestations resemble key features of diseases resulting from impaired function of airway cilia.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2742-y) contains supplementary material, which is available to authorized users.

Confirmation of a non-synonymous SNP in PNPLA8 as a candidate causal mutation for Weaver syndrome in Brown Swiss cattle

Background

Bovine progressive degenerative myeloencephalopathy (Weaver syndrome) is a neurodegenerative disorder in Brown Swiss cattle that is characterized by progressive hind leg weakness and ataxia, while sensorium and spinal reflexes remain unaffected. Although the causal mutation has not been identified yet, an indirect genetic test based on six microsatellite markers and consequent exclusion of Weaver carriers from breeding have led to the complete absence of new cases for over two decades. Evaluation of disease status by imputation of 41 diagnostic single nucleotide polymorphisms (SNPs) and a common haplotype published in 2013 identified several suspected carriers in the current breeding population, which suggests a higher frequency of the Weaver allele than anticipated. In order to prevent the reemergence of the disease, this study aimed at mapping the gene that underlies Weaver syndrome and thus at providing the basis for direct genetic testing and monitoring of today’s Braunvieh/Brown Swiss herds.

Results

Combined linkage/linkage disequilibrium mapping on Bos taurus chromosome (BTA) 4 based on Illumina Bovine SNP50 genotypes of 43 Weaver-affected, 31 Weaver carrier and 86 Weaver-free animals resulted in a maximum likelihood ratio test statistic value at position 49,812,384 bp. The confidence interval (0.853 Mb) determined by the 2-LOD drop-off method was contained within a 1.72-Mb segment of extended homozygosity. Exploitation of whole-genome sequence data from two official Weaver carriers and 1145 other bulls that were sequenced in Run4 of the 1000 bull genomes project showed that only a non-synonymous SNP (rs800397662) within the PNPLA8 gene at position 49,878,773 bp was concordant with the Weaver carrier status. Targeted SNP genotyping confirmed this SNP as a candidate causal mutation for Weaver syndrome. Genotyping for the candidate causal mutation in a random sample of 2334 current Braunvieh animals suggested a frequency of the Weaver allele of 0.26 %.

Conclusions

Through combined use of exhaustive sequencing data and SNP genotyping results, we were able to provide evidence that supports the non-synonymous mutation at position 49,878,773 bp as the most likely causal mutation for Weaver syndrome. Further studies are needed to uncover the exact mechanisms that underlie this syndrome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0201-5) contains supplementary material, which is available to authorized users.

Clinical and biochemical signs in Fleckvieh cattle with genetically confirmed Fanconi-Bickel syndrome (cattle homozygous for Fleckvieh haplotype 2)

Fanconi-Bickel Syndrome (FBS) is an autosomal recessive disorder of the carbohydrate metabolism, which has been reported in human and some animals (OMIA 000366-9913). In Fleckvieh cattle it is caused by mutations in SLC2A2, a gene encoding for glucose transporter protein 2 (GLUT2), which is primarily expressed in liver, kidney, pancreas and intestines. The causal mutation resides in a previously reported Fleckvieh Haplotype 2 (FH-2). FH-2 homozygous individuals are rare, but due to widespread use of heterozygous bulls in artificial insemination, heterozygous animals are likely to be present in a larger number in the cattle population. Two clinical cases of Fleckvieh cattle with a syndrome resembling the phenotypic appearance of FBS are presented in the present study describing the association between the clinical manifestations of FBS and the postulated frameshift mutation in bovine SLC2A2. Clinical examination showed poor growth, retarded development, polyuria, and polydipsia. Laboratory analyses showed an increased plasma glucose but normal insulin concentration and increased renal glucose excretion. Histopathological examination of kidney and liver samples revealed massively increased liver glycogen storage and nephrosis. Sires of both cases were tested positive for being heterozygous carriers for the same frameshift mutation in SLC2A2 as was originally reported in Fleckvieh cattle. DNA of both cases described was analyzed and Sanger sequencing confirmed homozygosity for the frameshift mutation in SLC2A2.

A multi-trait meta-analysis with imputed sequence variants reveals twelve QTL for mammary gland morphology in Fleckvieh cattle

Background

The availability of whole-genome sequence data from key ancestors in bovine populations provides an exhaustive catalogue of polymorphic sites that segregate within and across cattle breeds. Sequence variants identified from the sequenced genome of key ancestors can be imputed into animals that have been genotyped using medium- and high-density genotyping arrays. Association analysis with imputed sequences, particularly when applied to multiple traits simultaneously, is a very powerful approach to detect candidate causal variants that underlie complex phenotypes.

Results

We used whole-genome sequence data from 157 key ancestors of the German Fleckvieh cattle population to impute 20,561,798 sequence variants into 10,363 animals that had (partly imputed) genotypes based on 634,109 single nucleotide polymorphisms (SNPs). Rare variants were more frequent among the sequence-derived than the array-derived genotypes. Association studies with imputed sequence variants were performed using seven correlated udder conformation traits as response variables. The calculation of an approximate multi-trait test statistic enabled us to detect 12 quantitative trait loci (QTL) (P < 2.97 × 10⁻⁹) that affect different morphological features of the mammary gland. Among the tested variants, the most significant associations were found for imputed sequence variants at 11 QTL, whereas the top association signal was observed for an array-derived variant at a QTL on bovine chromosome 14. Seven QTL were associated with multiple phenotypes. Most QTL were located in non-coding regions of the genome but in close proximity of candidate genes that could be involved in mammary gland morphology (SP5, GC, NPFFR2, CRIM1, RXFP2, TBX5, RBM19 and ADAM12).

Conclusions

Using imputed sequence variants in association analyses allows the detection of QTL at maximum resolution. Multi-trait approaches can reveal QTL that are not detected in single-trait association studies. Most QTL for udder conformation traits were located in non-coding regions of the genome, which suggests that mutations in regulatory sequences are the major determinants of variation in mammary gland morphology in cattle.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0190-4) contains supplementary material, which is available to authorized users.

A comparison of methods to calculate a total merit index using stochastic simulation

BACKGROUND

Modern dairy cattle breeding goals include several production and more and more functional traits. Estimated breeding values (EBV) that are combined in the total merit index usually come from single-trait models or from multivariate models for groups of traits. In most cases, a multivariate animal model based on phenotypic data for all traits is not feasible and approximate methods based on selection index theory are applied to derive the total merit index. Therefore, the objective of this study was to compare a full multitrait animal model with two approximate multitrait models and a selection index approach based on simulated data.

METHODS

Three production and two functional traits were simulated to mimic the national Austrian Brown Swiss population. The reference method for derivation of the total merit index was a multitrait evaluation based on all phenotypic data. Two of the approximate methods were variations of an approximate multitrait model that used either yield deviations or de-regressed breeding values. The final method was an adaptation of the selection index method that is used in routine evaluations in Austria and Germany. Three scenarios with respect to residual covariances were set up: residual covariances were equal to zero, or half of or equal to the genetic covariances.

RESULTS

Results of both approximate multitrait models were very close to those of the reference method, with rank correlations of 1. Both methods were nearly unbiased. Rank correlations for the selection index method showed good results when residual covariances were zero but correlations with the reference method decreased when residual covariances were large. Furthermore, EBV were biased when residual covariances were high.

CONCLUSIONS

We applied an approximate multitrait two-step procedure to yield deviations and de-regressed breeding values, which led to nearly unbiased results. De-regressed breeding values gave even slightly better results. Our results confirmed that ignoring residual covariances when a selection index approach is applied leads to remarkable bias. This could be relevant in terms of selection accuracy. Our findings suggest that the approximate multitrait approach applied to de-regressed breeding values can be used in routine genetic evaluation.

Homozygous haplotype deficiency reveals deleterious mutations compromising reproductive and rearing success in cattle

Background

Cattle breeding populations are susceptible to the propagation of recessive diseases. Individual sires generate tens of thousands of progeny via artificial insemination. The frequency of deleterious alleles carried by such sires may increase considerably within few generations. Deleterious alleles manifest themselves often by missing homozygosity resulting from embryonic/fetal, perinatal or juvenile lethality of homozygotes.

Results

A scan for homozygous haplotype deficiency in 25,544 Fleckvieh cattle uncovered four haplotypes affecting reproductive and rearing success. Exploiting whole-genome resequencing data from 263 animals facilitated to pinpoint putatively causal mutations in two of these haplotypes. A mutation causing an evolutionarily unlikely substitution in SUGT1 was perfectly associated with a haplotype compromising insemination success. The mutation was not found in homozygous state in 10,363 animals (P = 1.79 × 10⁻⁵) and is thus likely to cause lethality of homozygous embryos. A frameshift mutation in SLC2A2 encoding glucose transporter 2 (GLUT2) compromises calf survival. The mutation leads to premature termination of translation and activates cryptic splice sites resulting in multiple exon variants also with premature translation termination. The affected calves exhibit stunted growth, resembling the phenotypic appearance of Fanconi-Bickel syndrome in humans (OMIM 227810), which is also caused by mutations in SLC2A2.

Conclusions

Exploiting comprehensive genotype and sequence data enabled us to reveal two deleterious alleles in SLC2A2 and SUGT1 that compromise pre- and postnatal survival in homozygous state. Our results provide the basis for genome-assisted approaches to avoiding inadvertent carrier matings and to improving reproductive and rearing success in Fleckvieh cattle.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1483-7) contains supplementary material, which is available to authorized users.

Genomic regions influencing coat color saturation and facial markings in Fleckvieh cattle

Genomic regions associated with coat color and pigmented areas of the head were identified for Fleckvieh (dual-purpose Simmental), a red-spotted and white-headed cattle breed. Coat color was measured with a chromameter, implementing the CIELAB color space and resulting in numerical representation of lightness, color intensity, red/green and blue/yellow color components, rather than subjective classification. Single marker regression analyses with fixed effects of the sex and barn were applied, and significant regions were determined with the local false discovery rate methodology. The PMEL and ERBB3 genes on chromosome 5 were in the most significant region for the color measurements. In addition to the blue/yellow color component and color intensity, the AP3B2 gene on chromosome 21 was identified. Its function was confirmed for similar traits in a range of model species. The KIT gene on chromosome 6 was found to be strongly associated with the inhibition of circum-ocular pigmentation and pigmented spots on the cheek.

A nonsense mutation in PLD4 is associated with a zinc deficiency-like syndrome in Fleckvieh cattle

BACKGROUND

Bovine hereditary zinc deficiency (BHZD) is an autosomal recessive disorder of cattle, first described in Holstein-Friesian animals. Affected calves suffer from severe skin lesions and show a poor general health status. Recently, eight calves with the phenotypic appearance of BHZD have been reported in the Fleckvieh cattle population.

RESULTS

In spite of the similar disease phenotypes, SLC39A4, the gene responsible for BHZD in Holstein-Friesian was excluded as underlying gene for the disorder in the affected Fleckvieh calves. In order to identify the disease-associated region, genotypes of eight affected calves obtained with the Illumina BovineHD BeadChip comprising 777,962 SNPs were contrasted with the genotypes of 1,339 unaffected animals. A strong association signal was observed on chromosome 21 (P = 5.87 × 10(-89)). Autozygosity mapping in the eight affected animals revealed a common segment of extended homozygosity encompassing 1,023 kb (BTA 21: 70,550,045 - 71,573,501). This region contains 17 genes/transcripts, among them two genes encoding gastro-intestinal zinc transporters (CRIP1, CRIP2). However, no mutation that was compatible with recessive inheritance could be detected in these candidate genes. One of the affected calves was re-sequenced together with 42 unaffected Fleckvieh animals. Analysis of the sequencing data revealed a nonsense mutation (p.W215X) in a phospholipase encoding gene (PLD4) as candidate causal polymorphism. To confirm the causality, genotypes of the p.W215X-mutation were obtained from 3,650 animals representing three different breeds. None of the unaffected animals was homozygous for the defect allele, while all eight affected calves were homozygous. The deleterious effect of the mutation is manifested in a significantly lower survival rate of descendants from risk matings when compared with the survival rate of descendants from non-risk matings. The deleterious allele has an estimated frequency of 1.1% in the Fleckvieh population.

CONCLUSION

Our results provide strong evidence that a newly identified recessive disorder in the Fleckvieh population is caused by a nonsense mutation in PLD4, most likely resulting in an impaired function of the encoded protein. Although the phenotype of affected calves strongly resembles BHZD, a zinc deficiency resulting from malabsorption is unlikely to be responsible for the diseased Fleckvieh calves.

Short communication: Genotyping of cows to speed up availability of genomic estimated breeding values for direct health traits in Austrian Fleckvieh (Simmental) cattle--genetic and economic aspects

The aim of this study was to quantify the impact of genotyping cows with reliable phenotypes for direct health traits on annual monetary genetic gain (AMGG) and discounted profit. The calculations were based on a deterministic approach using ZPLAN software (University of Hohenheim, Stuttgart, Germany). It was assumed that increases in reliability of the total merit index (TMI) of 5, 15, and 25 percentage points were achieved through genotyping 5,000, 25,000, and 50,000 cows, respectively. Costs for phenotyping, genotyping, and genomic estimated breeding values vary between €150 and €20 per cow. The gain in genotyping cows for traits with medium to high heritability is more than for direct health traits with low heritability. The AMGG is increased by 1.5% if the reliability of TMI is 5 percentage points higher (i.e., 5,000 cows genotyped) and 6.53% higher AMGG can be expected when the reliability of TMI is increased by 25 percentage points (i.e., 50,000 cows genotyped). The discounted profit depends not only on the costs of genotyping but also on the population size. This study indicates that genotyping cows with reliable phenotypes is feasible to speed up the availability of genomic estimated breeding values for direct health traits. But, because of the huge amount of valid phenotypes and genotypes needed to establish an efficient genomic evaluation, it is likely that financial constraints will be the main limiting factor for implementation into breeding program such as Fleckvieh Austria.

A nonsense mutation in TMEM95 encoding a nondescript transmembrane protein causes idiopathic male subfertility in cattle

Genetic variants underlying reduced male reproductive performance have been identified in humans and model organisms, most of them compromising semen quality. Occasionally, male fertility is severely compromised although semen analysis remains without any apparent pathological findings (i.e., idiopathic subfertility). Artificial insemination (AI) in most cattle populations requires close examination of all ejaculates before insemination. Although anomalous ejaculates are rejected, insemination success varies considerably among AI bulls. In an attempt to identify genetic causes of such variation, we undertook a genome-wide association study (GWAS). Imputed genotypes of 652,856 SNPs were available for 7962 AI bulls of the Fleckvieh (FV) population. Male reproductive ability (MRA) was assessed based on 15.3 million artificial inseminations. The GWAS uncovered a strong association signal on bovine chromosome 19 (P = 4.08×10⁻⁵⁹). Subsequent autozygosity mapping revealed a common 1386 kb segment of extended homozygosity in 40 bulls with exceptionally poor reproductive performance. Only 1.7% of 35,671 inseminations with semen samples of those bulls were successful. None of the bulls with normal reproductive performance was homozygous, indicating recessive inheritance. Exploiting whole-genome re-sequencing data of 43 animals revealed a candidate causal nonsense mutation (rs378652941, c.483C>A, p.Cys161X) in the transmembrane protein 95 encoding gene TMEM95 which was subsequently validated in 1990 AI bulls. Immunohistochemical investigations evidenced that TMEM95 is located at the surface of spermatozoa of fertile animals whereas it is absent in spermatozoa of subfertile animals. These findings imply that integrity of TMEM95 is required for an undisturbed fertilisation. Our results demonstrate that deficiency of TMEM95 severely compromises male reproductive performance in cattle and reveal for the first time a phenotypic effect associated with genomic variation in TMEM95.

Impaired male fertility is a prevalent condition in many species and is often explained by aberrant semen quality. In some cases, male fertility is severely compromised although semen quality is without any apparent pathological findings (i.e., idiopathic male subfertility). The genetic mechanisms underlying idiopathic male subfertility often remain unexplained. In the present paper, we report a recessively inherited variant of idiopathic male subfertility in a cattle population. We use 650,000 genome-wide SNP markers genotyped in >7900 artificial insemination bulls to pinpoint the underlying genomic region. We take advantage of whole-genome re-sequencing data of 43 animals to identify a causal loss-of-function mutation in TMEM95 encoding a nondescript transmembrane protein. We demonstrate that transmembrane protein 95 is located at the plasma membrane of spermatozoa of fertile animals whereas it is absent in spermatozoa of subfertile animals. Our results indicate that integrity of transmembrane protein 95 is required for an undisturbed fertilisation. This is the first report to reveal a phenotypic effect associated with genomic variation in TMEM95 in any organism.

Hot topic: Effect of breeding strategies using genomic information on fitness and health

A complex deterministic approach was used to model the breeding goal and breeding structure for the Austrian Fleckvieh (dual-purpose Simmental) breed. The reference breeding goal corresponded to the current total merit index (TMI-R), where dairy traits have a relative weight of 37.9% and fitness traits of 43.7% (beef traits 16.5%; milkability 2%). The breeding program was characterized by 280,000 cows under performance recording, 3,200 bull dams, 100 test bulls with a test capacity of 25%, and 15 proven bulls and 8 bull sires per year. The annual monetary genetic gain (AMGG) was generated mainly by increases in milk fat and milk protein yield (80.6%) and only to a small extent by fitness traits (6.6%). The inclusion of direct health traits (early reproductive disorders, cystic ovaries, and mastitis) with their economic weights increased the relative AMGG for fitness traits from 6.6 to 11.2%. The presently slightly negative AMGG for fertility index and udder health changed in a positive direction. Increasing the weight on the direct health traits by 50% resulted in a further shift toward fitness and health. The effect of strategies using genomic information in a total merit index (TMI) with varying weights on fitness and health traits was also analyzed. The conventional progeny-testing scheme was defined as the reference breeding program. A breeding program was considered to be genomically enhanced (GS50) when 50% of inseminations of herdbook cows and of bull dams were from young bulls with a genomic TMI, and a second program (GS100) did not rely on progeny-tested bulls at all. For GS50, a clear shift of the relative gain in AMGG toward fitness and health traits was observed for all 3 TMI scenarios, as a result of larger progeny groups and a shorter generation interval. For GS100, where no gene flow from progeny-tested bulls was assumed, the genetic gain per generation was lower for the fertility and udder health index but higher per year. The results based on natural genetic gain per year showed that no positive genetic response for fertility and udder health index were achieved for TMI-R (without the inclusion of direct health traits) in GS50 and GS100. The direction of the genetic trend was determined by the weights given to fertility and udder health indices within the TMI. When appropriate weights generated a clear positive trend, GS50 and GS100 reinforced this trend.

Short-term intermittent hypoxia reduces the severity of acute mountain sickness

Intermittent hypoxia (IH) is a promising approach to induce acclimatization and hence lower the risk of developing acute mountain sickness (AMS). We hypothesized that a short-term IH protocol in normobaric hypoxia (7 × 1 h to 4500 m) effectively increases the hypoxic ventilatory response (HVR) and reduces the incidence and severity of AMS. Therefore, 26 men (25.5 ± 4.4 years), assigned in a double-blinded fashion to the hypoxia group (HG) or placebo group (PG), spent 8 h at 5300 m before (PRE) and 2 days after cessation of the IH protocol (POST). Measurements included the evaluation of the Lake Louise Score (LLS) and the HVR. The severity of AMS decreased from PRE to POST in the HG (from 6.0 ± 2.7 at PRE to 4.1 ± 2.1 at POST), whereas the LLS in the PG stayed high (from 5.7 ± 2.9 to 5.5 ± 2.8, respectively). The HVR in the HG increased from 0.73 ± 0.4 L/min/% at PRE to 1.10 ± 0.5 L/min/% at POST and did not increase in the PG. The reduction of the LLS was inversely related to the changes in the HVR (r = -0.434), but the AMS incidence was not different between the HG and the PG at POST. In conclusion, short-term IH reduced the severity of AMS development during a subsequent 8-h exposure to normobaric hypoxia.

Combining evidence of selection with association analysis increases power to detect regions influencing complex traits in dairy cattle

Background

Hitchhiking mapping and association studies are two popular approaches to map genotypes to phenotypes. In this study we combine both approaches to complement their specific strengths and weaknesses, resulting in a method with higher statistical power and fewer false positive signals. We applied our approach to dairy cattle as they underwent extremely successful selection for milk production traits and since an excellent phenotypic record is available. We performed whole genome association tests with a new mixed model approach to account for stratification, which we validated via Monte Carlo simulations. Selection signatures were inferred with the integrated haplotype score and a locus specific permutation based integrated haplotype score that works with a folded frequency spectrum and provides a formal test of signifance to identify selection signatures.

Results

About 1,600 out of 34,851 SNPs showed signatures of selection and the locus specific permutation based integrated haplotype score showed overall good accordance with the whole genome association study. Each approach provides distinct information about the genomic regions that influence complex traits. Combining whole genome association with hitchhiking mapping yielded two significant loci for the trait protein yield. These regions agree well with previous results from other selection signature scans and whole genome association studies in cattle.

Conclusion

We show that the combination of whole genome association and selection signature mapping based on the same SNPs increases the power to detect loci influencing complex traits. The locus specific permutation based integrated haplotype score provides a formal test of significance in selection signature mapping. Importantly it does not rely on knowledge of ancestral and derived allele states.

A novel mutation in the maternally imprinted PEG3 domain results in a loss of MIMT1 expression and causes abortions and stillbirths in cattle (Bos taurus)

Congenital malformations resulting in late abortions and stillbirths affect the economic wellbeing of producers and the welfare of cattle in breeding programs. An extremely high incidence of stillbirths of “half-sized” calves of normal karyotype and uninflated lungs was diagnosed in the progeny of the Finnish Ayrshire (Bos taurus) bull - YN51. No other visible anatomical abnormalities were apparent in the stillborn calves. We herein describe the positional identification of a 110 kb microdeletion in the maternally imprinted PEG3 domain that results in a loss of paternal MIMT1 expression and causes late term abortion and stillbirth in cattle. Using the BovineSNP50 BeadChip we performed a genome-wide half-sib linkage analysis that identified a 13.3 Mb associated region on BTA18 containing the maternally imprinted PEG3 domain. Within this cluster we found a 110 kb microdeletion that removes a part of the non-protein coding MER1 repeat containing imprinted transcript 1 gene (MIMT1). To confirm the elimination of gene expression in calves inheriting this deletion, we examined the mRNA levels of the three maternally imprinted genes within the PEG3 domain, in brain and cotyledon tissue collected from eight fetuses sired by the proband. None of the fetuses that inherited the microdeletion expressed MIMT1 in either tissue. The mutation, when inherited from the sire, is semi-lethal for his progeny with an observed mortality rate of 85%. The survival of 15% is presumably due to the incomplete silencing of maternally inherited MIMT1 alleles. We designed a PCR-based assay to confirm the existence of the microdeletion in the MIMT1 region that can be used to assist cattle breeders in preventing the stillbirths.

Haplotypes of the porcine peroxisome proliferator-activated receptor delta gene are associated with backfat thickness

BACKGROUND

Peroxisome proliferator-activated receptor delta belongs to the nuclear receptor superfamily of ligand-inducible transcription factors. It is a key regulator of lipid metabolism. The peroxisome proliferator-activated receptor delta gene (PPARD) has been assigned to a region on porcine chromosome 7, which harbours a quantitative trait locus for backfat. Thus, PPARD is considered a functional and positional candidate gene for backfat thickness. The purpose of this study was to test this candidate gene hypothesis in a cross of breeds that were highly divergent in lipid deposition characteristics.

RESULTS

Screening for genetic variation in porcine PPARD revealed only silent mutations. Nevertheless, significant associations between PPARD haplotypes and backfat thickness were observed in the F2 generation of the Mangalitsa x Piétrain cross as well as a commercial German Landrace population. Haplotype 5 is associated with increased backfat in F2 Mangalitsa x Piétrain pigs, whereas haplotype 4 is associated with lower backfat thickness in the German Landrace population. Haplotype 4 and 5 carry the same alleles at all but one SNP. Interestingly, the opposite effects of PPARD haplotypes 4 and 5 on backfat thickness are reflected by opposite effects of these two haplotypes on PPAR-delta mRNA levels. Haplotype 4 significantly increases PPAR-delta mRNA levels, whereas haplotype 5 decreases mRNA levels of PPAR-delta.

CONCLUSION

This study provides evidence for an association between PPARD and backfat thickness. The association is substantiated by mRNA quantification. Further studies are required to clarify, whether the observed associations are caused by PPARD or are the result of linkage disequilibrium with a causal variant in a neighbouring gene.

Variation in neighbouring genes of the dopaminergic and serotonergic systems affects feather pecking behaviour of laying hens

Feather pecking is a behavioural disorder of laying hens and has serious animal welfare and economic implications. One of the several aetiological hypotheses proposes that the disorder results from redirected exploratory behaviour. Variation in the gene encoding the dopamine D4 receptor (DRD4) has been shown to be associated with exploratory behaviour in several species, including in a passerine bird species. We therefore considered DRD4 as a candidate gene for feather pecking. We have annotated DRD4 in the chicken genome and have re-sequenced it in 140 animals belonging to: experimental layer lines divergently selected for high and low propensity to feather pecking; the unselected founder population; and two commercial lines with low and high propensity to feather pecking. We have identified two sub-haplotypes of DRD4 that are highly significantly associated with feather pecking behaviour in the experimental (P = 7.30 x 10(-7)) as well as in the commercial lines (P = 2.78 x 10(-6)). Linkage disequilibrium (LD) extends into a neighbouring gene encoding deformed epidermal autoregulatory factor 1 (DEAF1). The product of DEAF1 regulates the transcription of the gene encoding the serotonin (5-hydroxytryptamine) 1A receptor. Thus, DEAF1 represents another candidate gene for feather pecking. Re-sequencing of five animals homozygous for the 'low-pecking' sub-haplotype and of six animals homozygous for the 'high-pecking' sub-haplotype delineated an LD block of 14 833 bases spanning the two genes. None of the variants in the LD block is obviously functional. However, the haplotype information will be useful to select against the propensity to feather pecking in chicken and to elucidate the functional implications of the variants.

Characterization of a 320-kb region containing the HEXA gene on bovine chromosome 10 and analysis of its association with BSE susceptibility

Bovine spongiform encephalopathy (BSE) belongs to a group of neurodegenerative diseases known as transmissible prion diseases. Recently, variants in the promoter region of the prion protein (PRNP) gene have been shown to have a considerable effect on the susceptibility to BSE. However, a previous genome scan revealed other putative BSE-susceptibility loci. Here, we analysed such a region on BTA10, which contains the functional candidate gene HEXA. Three hundred and twenty kilobases that, besides HEXA, also contain ARIH1, BRUNOL6 and PARP6 were characterized and screened for polymorphisms. Genotyping of 38 SNPs in Holstein-Friesian animals from the UK (350 diseased and 270 controls) revealed two intronic SNPs that were associated with BSE incidence, with experiment-wise P-values of 3.5 x 10(-3) and 7.7 x 10(-3) respectively. Both SNPs were in strong linkage disequilibrium and the rare alleles had a protective effect. These alleles were contained in a haplotype dubbed 'UK-protective' that was significantly overrepresented in the controls with a permuted P-value of 2 x 10(-3). An association study in German Holstein animals (73 diseased and 627 controls) revealed an opposite effect of the 'UK-protective' haplotype in this population, i.e. it was overrepresented in the diseased animals, although not significant after correction for multiple testing. These findings indicate a causal variant for BSE susceptibility on BTA10 in linkage disequilibrium with the markers studied. Candidate gene analyses of the surrounding region and additional association studies will help to clarify the origin of the protective effects and to identify causal variants for BSE susceptibility on BTA10.

Genetic parameters for semen production traits in Austrian dual-purpose Simmental bulls

Genetic parameters were estimated for semen production traits collected in an Austrian AI centre in the years 2000-2004. In total, 12,746 ejaculates from 301 Austrian dual-purpose Simmental (Fleckvieh) AI bulls were examined considering different effects on ejaculate volume, sperm concentration, percentage of viable spermatozoa in the ejaculate, total spermatozoa per ejaculate and motility. The model for genetic parameter estimation included the fixed effects age of bull, collection interval, number of collections on collection day, bull handler, semen collector, year and month of collection, a random additive genetic component and a permanent environmental effect. Correlations between estimated breeding values for semen traits and male fertility from the routine evaluation were calculated. The fertility trait considered in the routine evaluation is non-return rate 90 for the first insemination. All semen production traits were moderately heritable. Heritabilities for volume, concentration, percentage of viable spermatozoa, total number of spermatozoa and motility were 0.18, 0.14, 0.10, 0.22 and 0.04, respectively. Correlations between breeding values for semen quality traits and routinely estimated breeding values for male fertility were low and ranged from 0.08 to 0.17 indicating that semen production traits are rather poor predictors of male fertility.

A major genetic component of BSE susceptibility

Background

Coding variants of the prion protein gene (PRNP) have been shown to be major determinants for the susceptibility to transmitted prion diseases in humans, mice and sheep. However, to date, the effects of polymorphisms in the coding and regulatory regions of bovine PRNP on bovine spongiform encephalopathy (BSE) susceptibility have been considered marginal or non-existent. Here we analysed two insertion/deletion (indel) polymorphisms in the regulatory region of bovine PRNP in BSE affected animals and controls of four independent cattle populations from UK and Germany.

Results

In the present report, we show that two previously reported 23- and 12-bp insertion/deletion (indel) polymorphisms in the regulatory region of bovine PRNP are strongly associated with BSE incidence in cattle. Genotyping of BSE-affected and control animals of UK Holstein, German Holstein, German Brown and German Fleckvieh breeds revealed a significant overrepresentation of the deletion alleles at both polymorphic sites in diseased animals (P = 2.01 × 10^-3 and P = 8.66 × 10^-5, respectively). The main effect on susceptibility is associated with the 12-bp indel polymorphism. Compared with non-carriers, heterozygous and homozygous carriers of the 12-bp deletion allele possess relatively higher risks of having BSE, ranging from 1.32 to 4.01 and 1.74 to 3.65 in the different breeds. These values correspond to population attributable risks ranging from 35% to 53%.

Conclusion

Our results demonstrate a substantial genetic PRNP associated component for BSE susceptibility in cattle. Although the BSE risk conferred by the deletion allele of the 12-bp indel in the regulatory region of PRNP is substantial, the main risk factor for BSE in cattle is environmental, i.e. exposure to feedstuffs contaminated with the infectious agent.

Effects of age and environmental factors on semen production and semen quality of Austrian Simmental bulls

More than 90% of the breeding stock of Austrian dual purpose Simmental cows is artificially inseminated. Knowledge of factors affecting sperm production and semen quality is of importance with regard to reproductive efficiency and thus genetic improvement as well as for the productivity and profitability of AI centers. Hence, semen data from two Austrian AI centres collected in the years 2000 and 2001 were evaluated. In total, 3625 and 3654 ejaculates from 147 and 127 AI bulls, respectively, were analysed regarding ejaculate volume, sperm concentration, percentage of viable spermatozoa in the ejaculate, total spermatozoa per ejaculate and motility. Effects accounted for were the bull (random), age of bull, collection interval, number of collection on collection day, bull handler, semen collector, temperature on day of semen collection, in the course of epididymal maturation (average temperature of days 1-11 before collection) and during spermatogenesis (average temperature of days 12-65 before collection). Age of bull significantly affected all traits (P<0.01 to P<0.001) except motility score in center 2. Ejaculate volume and total number of spermatozoa increased with age of bull while sperm concentration was lower in higher age classes (center 1). The collection team was also found to significantly influence semen quality traits. With increasing collection interval ejaculate volume and total number of spermatozoa increased significantly (P<0.05 to P<0.001) while collection intervals between 4-9 days and 1-6 days were superior with regard to sperm concentration and percentage of viable spermatozoa, respectively (P<0.10 to P<0.001). First ejaculates were superior with respect to ejaculate volumes, sperm concentrations and total number of spermatozoa per ejaculate (P<0.001). Temperature, either on day of semen collection or during epididymal maturation or spermatogenesis, had important but inconsistent effects on semen production and sperm quality. Overall, however, ambient temperatures in the range of 5-15 degrees C were found to be optimal for semen production.