A missense mutation in TUBD1 is associated with high juvenile mortality in Braunvieh and Fleckvieh cattle
- Hermann Schwarzenbacher1,
- Johann Burgstaller2,
- Franz R. Seefried3,
- Christine Wurmser4,
- Monika Hilbe5,
- Simone Jung4,
- Christian Fuerst1,
- Nora Dinhopl6,
- Herbert Weissenböck6,
- Birgit Fuerst-Waltl7,
- Marlies Dolezal8,
- Reinhard Winkler9,
- Oskar Grueter10,
- Ulrich Bleul11,
- Thomas Wittek2,
- Ruedi Fries4 and
- Hubert Pausch4Email author
© The Author(s). 2016
Received: 1 March 2016
Accepted: 14 May 2016
Published: 25 May 2016
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.
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−12). 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−15). 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.
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.
Reducing calf mortality is an important objective in cattle breeding populations both for animal welfare and economic reasons. Bovine calf losses including late abortions, stillbirths and diseases during rearing may be in excess of 10 % [1–4]. Most rearing losses are attributable to infectious diarrheal and respiratory diseases . While a number of genetic variants predisposing to high perinatal mortality have been identified in bovine populations (e.g., [6, 7]), the mapping of loci affecting disease susceptibility and rearing success is difficult because of their low heritability [8, 9].
Case-control association testing using genome-wide marker data facilitated the discovery of causal variants for monogenic disorders that result in high juvenile mortality [10–16]. Although the pathophysiology of such conditions is heterogeneous, pertinently affected calves may be apathetic, without vigor, retarded in growth or highly susceptible to infectious disease. Affected calves may be born without typical signs of disease but develop clinical features shortly after birth [10, 11, 14].
The availability of large-scale genotype data enables the identification of haplotypes with homozygosity depletion [14, 17, 18]. Haplotypes with reduced or missing homozygosity in adult animals are likely to harbor recessive alleles that are associated with an increased pre-, peri- or postnatal mortality. Genome-wide scans for homozygous haplotype deficiency in North American  and European Braunvieh  populations uncovered a short segment on bovine chromosome 19 that was associated with high postnatal calf losses (OMIA 001939-9913). The segment associated with calf mortality was denominated 424.49 and BTA19-1 in the North American  and European Braunvieh  populations, respectively. Since 2013, the associated genomic segment is consistently referred to as BH2 with BH being an abbreviation for Braunvieh haplotype .
In the present study, we exploit comprehensive genotyping and sequencing data to detect a missense mutation in TUBD1 that is most likely causal for the high juvenile mortality of BH2 homozygous calves. We validate the missense mutation in an independent cattle population and provide evidence that homozygous calves suffer from chronic respiratory disease.
BH2 compromises juvenile survival in Braunvieh cattle
Harmful effects of two haplotypes with reduced homozygosity in Braunvieh (BH2) and Fleckvieh (BH2FV) cattle
Haplotype frequency (%)
Number of calvings
Difference in survival rate
Two male Braunvieh calves homozygous for BH2 were detected in the Swiss bovine genome database. According to the possessing farmers, both calves were underweight at birth. One calf died for an unknown cause at 72 days of age and was not available for phenotyping. The second calf (BH2hom) was referred to the animal clinic at 66 days of age. During a hospitalization period of 95 days, BH2hom suffered repeatedly from bronchopneumonia. At 161 days of age, BH2hom was euthanized because of suddenly occurring severe dyspnea.
A missense mutation in TUBD1 is compatible with recessive inheritance of BH2
Genotypes for a missense mutation in TUBD1 in 661 adult Braunvieh cattle
BH2 haplotype status
Genotype at rs383232842
The BH2 haplotype affects rearing success in another cattle breed
The analysis of whole-genome sequencing data from 236 animals from six bovine breeds other than Braunvieh revealed that rs383232842 was heterozygous in seven out of 149 Fleckvieh animals (Additional file 1 and 3). To obtain genotypes for rs383232842 in a representative sample of the Fleckvieh population, we genotyped 3807 randomly selected adult animals using a KASP genotyping assay. Among those, 142 were heterozygous and none was homozygous for the C-allele corresponding to a C-allele frequency of 1.86 %. Of 3807 animals that had been genotyped at rs383232842, 1,974 were also genotyped with the Illumina BovineHD Bead chip comprising 777,962 SNPs. The availability of dense marker data enabled us to detect an 1.05 Mb haplotype (BH2FV; from 10,406,017 bp to 11,453,904 bp) that was specific for the rs383232842 C-allele: ninety-two animals that carried the C-allele were heterozygous carriers of BH2FV, whereas 1882 animals that were homozygous for the reference allele did not carry BH2FV. In 37,314 Fleckvieh animals of the national bovine genome databases that had been genotyped with the Illumina BovineSNP50 Bead chip, the frequency of BH2FV was 1.68 %. Only four animals were homozygous for BH2FV although ten were expected (P = 0.09). The analysis of 46,566 calving records revealed that the first-year mortality of descendants from BH2FV risk matings was increased by 6.1 % compared to progeny from non-risk matings (P = 2.2x10−15) (Table 1).
A common haplotype with the rs383232842 C-allele segregates in Braunvieh and Fleckvieh cattle
BH2 encompasses a 1.14 Mb segment from 10,694,269 bp to 11,833,182 bp. It is defined by eleven SNPs of the Illumina BovineSNP50 Bead chip. BH2FV, defined by 249 SNPs of the Illumina BovineHD Bead chip, is located between 10,406,017 bp and 11,453,904 bp. The BH2 alleles are a subset of the BH2FV alleles within a 741 kb segment between 10,694,269 bp and 11,435,269 bp (Additional file 6). Since both BH2 and BH2FV contain the rs383232842 C-allele, a common origin of the two haplotypes is most likely. The analysis of pedigree records enabled us to track BH2 and BH2FV back to the Braunvieh bull Rancho Rustic My Design (birth year 1963) and the Fleckvieh bull Polzer (birth year 1959), respectively. However, we were not able to identify common ancestors in their pedigrees possibly because of missing pedigree information from very distantly related relatives. Cluster analyses using genome-wide marker data of 442 Fleckvieh and 280 Braunvieh animals born between 1970 and 2010 revealed no evidence for an exchange of genetic material between both breeds (Additional file 7).
BH2 segregates in the Holstein breed, but does not contain the rs383232842 C-allele
To test if the common BH2/BH2FV haplotype segment segregated in another dairy breed, we analyzed haplotypes of 8840 adult Holstein animals that had been genotyped with the Illumina BovineSNP50 Bead chip. We identified 917 heterozygotes and 30 homozygotes for BH2 among the genotyped Holstein animals. The observed haplotype distribution did not deviate from the Hardy-Weinberg equilibrium (P = 0.54). We additionally genotyped the rs383232842 polymorphism in 503 Holstein bulls using a KASP genotyping assay. All of them including 74 heterozygotes and two homozygotes for BH2 were homozygous for the rs383232842 reference allele (Additional file 6). Another 311 Holstein animals that had been sequenced within the 1000 bull genomes project  were homozygous for the reference allele indicating that the rs383232842 C-allele does not segregate in this breed. The non-reference alleles at two non-coding variants located within the BH2 interval were only observed at very low frequency in Holstein cattle (Additional file 1 & Additional file 5).
Identification of BH2 homozygous calves
Sampling of calves from BH2 risk matings
Genotype at rs383232842
Number of calves
Died within 10 days
Died within 30 days
During the sampling of BH2 homozygous calves, we encountered a male Fleckvieh calf (FV1) with low birth weight and postnatal growth restriction that descended from a BH2FV carrier bull. At 40 days of age, its weight was 28.5 kg, which is only a third compared to healthy Fleckvieh calves of the same age. FV1 was tested negatively for two recessive mutations that are known to cause retarded growth [14, 22]. Another 16-months old Braunvieh bull (BV5) with poor growth performance and a history of recurrent respiratory disease was reported by a farmer. Inspection of the bull’s pedigree revealed BH2 carriers among its paternal and maternal ancestors. Sanger sequencing of DNA samples from FV1 and BV5 confirmed that both animals were homozygous for the rs383232842 C-allele. Both animals were referred to the animal clinic.
Homozygous animals suffer from chronic airway disease
Characteristics of six hospitalized animals
Age at admission (days)
Weight at admission (kg)
Pulse (heartbeats per minute)
Body temperature (°C)
Respiratory rate (breaths per minute)
Hospitalization period (days)
Age at euthanasia (days)
Weight at euthanasia (kg)
Homozygosity for BH2 may be incompletely penetrant
Breeding consultants reported another three Braunvieh animals that were homozygous for BH2. These homozygous animals were detected during routine genomic evaluation using genotypes of the Illumina BovineSNP50 Bead chip. Sanger sequencing of DNA samples confirmed homozygosity for the rs383232842 C-allele in all animals. While one young bull (BV6) was healthy and normally developed at the age of 370 days, two animals (BV7, BV8) were retarded in growth and had a history of recurrent respiratory disease (Additional file 13). The three homozygous animals were also homozygous for the alternate allele at two non-coding variants in LD with rs383232842 (Additional file 5).
We show that the proximal region of bovine chromosome 19 is associated with high peri- and postnatal calf losses in Braunvieh cattle, corroborating previous reports of a harmful haplotype (Braunvieh haplotype 2, BH2) in European and North-American Braunvieh populations [17, 19]. The joint analysis of homozygosity depletion and insemination and rearing data revealed a 1.14 Mb segment between 10,694,269 bp and 11,833,182 bp on BTA19 as the most likely interval harboring the causal mutation, which agrees with previous findings [17, 19]. The first-year mortality of calves descending from BH2 risk matings is less than expected for a lethal recessive disorder, indicating that some BH2 homozygous calves may reach adulthood. Despite its detrimental effect on juvenile survival, the frequency of BH2 is high in the European Braunvieh populations (6.5 %). Based on the retrospective analysis of inadvertent BH2 riskmatings that happened in the Austrian Braunvieh population (N = 50,000 calvings per year) between 1990 and 2013, we estimated that between 11 and 393 BH2 homozygous animals were born annually in the past 25 years. Deleterious alleles may reach high frequency in cattle populations because of the widespread use of unnoticed carrier bulls in artificial insemination and pleiotropic effects on desirable traits [23–25]. So far, there is no evidence that BH2 has desirable effects on important breeding objectives. Several harmful conditions other than BH2 segregate in the Braunvieh populations [17, 26–29]. Pinpointing the causal mutation allows for the implementation of efficient genome-based mating programs to avoid the inadvertent mating of carrier animals while maintaining genetic diversity and high rates of genetic gain. A missense mutation (rs383232842, p.H210R) in TUBD1 was found to be compatible with recessive inheritance of BH2. Affected, carrier and non-carrier animals were homozygous for the alternate allele, heterozygous and homozygous, respectively, for the reference allele at rs383232842. Although BH2 is in high LD with rs383232842, the identification of carrier animals is less reliable with haplotype information. Five out of 661 animals were misclassified based on the haplotype test indicating both lower specificity and sensitivity of haplotype-based classification of animals into carriers and non-carriers which agrees with previous findings in cattle [10, 18, 26, 30].
Pinpointing causal mutations is often difficult as many sequence variants may be compatible with the presumed pattern of inheritance (e.g., [18, 24]). In 54 sequenced Braunvieh animals of our study, two intergenic variants were in complete LD with rs383232842. Yet we consider rs383232842 as the most likely to be causal because it results in a substitution of a conserved amino acid residue. Moreover, all compatible non-coding variants are far away from annotated transcripts and are thus less likely to be functionally important. However, it cannot be ruled out a priori that such variants affect gene regulation, e.g., if located in a distant enhancer region. A criterion for causality is often that the segregation of the candidate variant is restricted to the population affected by the condition (e.g., [13, 14, 29, 31]). Strictly applying this criterion to our data would have led to the exclusion of rs383232842 as a plausible causal mutation for BH2, because it segregates in Fleckvieh cattle where a previous scan for homozygous haplotype deficiency did not detect homozygosity depletion at the BH2 region . In the present study, we specifically tested the BH2FV haplotype with the rs383232842 C-allele for association with juvenile mortality in more than 37,000 Fleckvieh animals. Homozygosity depletion at BH2FV was not significant (P = 0.09) possibly because most Fleckvieh animals of the present study were genotyped shortly after birth, before typical signs of disease may become manifest. However, the first-year mortality of descendants is significantly higher in BH2FV risk than non-risk matings (P = 2.2x10−15) corroborating a detrimental effect of BH2FV.
The exchange of genetic material between breeds may result in the manifestation of the same recessively inherited conditions in different populations . Kumar Kadri et al. identified a sequence variant that compromises fertility in three admixed Nordic dairy breeds . Drögemüller et al. mapped a sequence variant for an autosomal recessive skin disease in two geographically separated albeit genetically related dog populations . A diagnostic haplotype associated with progressive degenerative myeloencephalopathy (Weaver syndrome) in Braunvieh cattle also segregates in the Holstein cattle breed . However, the most likely causal variant for Weaver syndrome was never found in haplotype carriers from the Holstein cattle breed indicating that an ancestral version of the haplotype without the Weaver mutation persists in Holstein cattle [32, 33]. The discovery of a recessive mutation that compromises juvenile survival in Fleckvieh and Braunvieh cattle demonstrates for the first time, that deleterious alleles may segregate across closed cattle populations, e.g., populations without recent admixture. Fleckvieh and Braunvieh animals with the rs383232842 C-allele share a common haplotype, indicating they inherited the corresponding chromosome segment from a common ancestor. However, our data neither revealed common ancestors nor provided evidence of an exchange of genetic material between both breeds. A version of the ancestral haplotype without the rs383232842 C-allele also segregates in Holstein cattle at relatively high frequency (5.5 %). These findings indicate that the rs383232842 mutation occurred in an ancient ancestor of Braunvieh and Fleckvieh cattle, however, after the divergence from Holstein cattle, which is compatible with the history of those breeds . While BH2 and BH2FV affect juvenile survival in Braunvieh and Fleckvieh cattle, respectively, the haplotype version without the rs383232842 C-allele does not exhibit homozygous haplotype deficiency in Holstein cattle supporting causality of the rs383232842 polymorphism.
To investigate the pathophysiology that may result from the rs383232842 mutation in TUBD1, we examined six homozygous animals. Although the age of the examined animals differed, the disease manifestation was homogeneous: recurrent airway disease resulted in a steadily declining general condition and poor growth performance despite excellent husbandry conditions and medical treatment at the animal clinic. Thus, a more severe disease progression in homozygous animals kept under normal husbandry conditions in conventional farms is likely. The respiratory manifestations resemble disease patterns arising from defective cilia of the respiratory tract [35–37]. We observed microtubular assembly defects in 20-30 % of the cilia of the respiratory epithelium in affected animals which is 4-6 times higher than in healthy individuals [38–40]. Ultrastructural defects of the cilia might also result from infectious or inflammatory diseases . Persistent coughing, tachypnea and airway mucus obstruction in homozygous animals might indicate an impaired mucociliary clearance due to an aberrant ciliary beat pattern [42–45]. We used a cytology brush to collect viable epithelial cells from the nasal mucosa of affected calves to investigate ciliary beat pattern and beat frequency (data not shown). However, all collected samples were excessively covered with mucous exudate and contained numerous inflammatory cells, precluding the analysis of ciliary beating using standard protocols.
Several sequence variants associated with ultrastructural aberrations of respiratory cilia cause chronic airway disease in humans, mice and dogs (e.g., [37, 46, 47]). Respiratory cilia are characterized by a “9x2 + 2” architecture consisting of a central pair of microtubule surrounded by nine outer microtubule doublets. The basal body is located at the base of the cilium where it attaches to the cell body. Tubulin delta 1 is required for proper microtubule polymerization in basal bodies [48–51]. A lack of TUBD1 causes ultrastructural microtubular defects in basal bodies in Chlamydomonas and Paramecium [48, 50, 51]. We also observed disorganized microtubules in airway cilia of animals homozygous for the rs383232842 mutation corroborating a crucial role of TUBD1 for proper assembly of cilia. To our knowledge, this is the first report of a phenotypic effect associated with genomic variation in TUBD1 in a mammalian species. A conserved histidine residue at amino acid position 210 in TUBD1 is involved in microtubule polymerization . Since the rs383232842 C-allele introduces a putatively damaging substitution of the histidine by an arginine, proper microtubule polymerization might be compromised in homozygous animals, resulting in ultrastructural defects of the airway cilia. Besides being involved in cilia assembly, microtubules constitute the bipolar spindle that binds and moves the chromosomes during the different meiotic phases. Disturbed microtubule formation may lead to erroneous segregation of chromosomes with negative effects on growth and differentiation (e.g., ). Such disturbances could be responsible for frequently observed prenatal growth retardation and stillbirth of fetuses from BH2 risk matings.
Most homozygous calves (75 %) were stillborn or died shortly after birth. Six hospitalized calves suffered from airway disease early in life. The general condition of all hospitalized animals declined steadily because the disease was not responsive to medical treatment. However, homozygosity for the rs383232842 mutation does not necessarily have fatal consequences as evidenced by a bull homozygous for the putative disease allele without any apparent signs of disease. Although minimizing exposure to environmental pollutants and respiratory pathogens may contribute to a less severe disease progression , our results suggest that the mutation is highly penetrant because 17 out of 18 homozygous animals were stillborn or diseased shortly after birth. In any case, our findings now enable to implement genome-based mating strategies to avoid the mating of carrier animals thereby preventing the birth of homozygous calves that are likely to suffer from chronic disease.
A missense mutation in TUBD1 is associated with high juvenile mortality in Braunvieh and, unexpectedly, Fleckvieh cattle. Since genomic analyses revealed no evidence for an exchange of genetic material between both breeds, we demonstrate for the first time that the same harmful recessive mutation may segregate in different closed populations, indicating that such mutations may have occurred many generations ago. Homozygous animals suffer from chronic airway disease possibly resulting from defective cilia in the respiratory tract. Our findings now enable the implementation of genome-based mating strategies in order to prevent animal suffering and economic losses.
Animals, genotypes, quality control and phasing
Genotypes of 8501 Braunvieh and 37,314 Fleckvieh animals were extracted from the Austrian and German bovine genome databases. Genotypes of 1765 Braunvieh bulls were provided by breeding organisations from Switzerland (n = 993), Italy (n = 465) and USA (n = 316) based on bilateral genotype exchange agreements. All animals had been genotyped with the Illumina BovineSNP50 bead chip (version 1 and version 2) comprising approximately 54,000 SNPs. The chromosomal position of the SNPs was based on the UMD 3.1 assembly of the bovine genome . Mitochondrial, X-chromosomal, Y-chromosomal SNPs and SNPs that had no known chromosomal position were not considered for further analyses. Animals with a call-rate less than 95 % were excluded as well as animals with mendelian conflicts with their genotyped sires at more than 200 SNPs. We excluded SNPs with a minor allele frequency below 0.5 % as well as SNPs with call rate less than 90 % or more than 0.5 % mendelian conflicts in sire-offspring pairs. The final Braunvieh data consisted of 8,446 animals and 35,382 autosomal SNPs with an average per-individual call-rate of 99.61 %. The final Fleckvieh data consisted of 37,314 animals and 41,027 autosomal SNPs with an average per-individual call-rate of 99.82 %. Haplotypes for BTA19 were inferred and sporadically missing genotypes were imputed using the shapeit software (v2.778) .
Mapping of BH2 in the Braunvieh cattle population
A sliding window with variable size (ranging from 0.75 to 8 Mb) was shifted along chromosome 19 (in steps of 0.5 x window size). Within each sliding window, the observed number of homozygotes was compared to the expected number for all haplotypes with a frequency above 3 % using Fisher exact tests. The expected number of homozygous animals was calculated using haplotype information from sire, maternal grandsire and haplotype frequency. Phenotypic effects were analyzed for all haplotypes that showed significant depletion of homozygous animals (P < 1x10−4).
Mapping of BH2FV in the Fleckvieh cattle population
To identify a haplotype specific for the rs383232842 C-allele, we exploited haplotypes of 10,363 Fleckvieh animals that had been genotyped and (partially) imputed at 652,853 SNPs . We extracted haplotypes of 1947 animals that had also been genotyped at rs383232842 (see below). A haplotype window consisting of five adjacent SNPs was centered at rs383232842 and increased by one SNP on either side as long as a common haplotype (BH2FV) was heterozygous in all animals that were heterozygous at rs383232842. To assess the frequency of BH2FV in a representative sample of the Fleckvieh cattle population, we analysed haplotypes of 37,314 Fleckvieh animals from the bovine genome databases (see above).
Mapping of BH2 in the Holstein cattle population
To test if the common BH2/BH2FV haplotype also segregates in the Holstein cattle breed, we used haplotypes of 8840 Holstein artificial insemination bulls that had been genotyped previously with the Illumina BovineSNP50 Bead chip . BH2 carriers were identified based on eight SNPs (from 10,694,269 bp to 11,435,269 bp) that were common in BH2/BH2FV carriers, respectively.
Phenotypic effects associated with BH2 and BH2FV
Insemination success was analyzed based on 55,837 artificial inseminations that were carried out with Braunvieh bulls that had a known haplotype status for BH2. The insemination success in risk matings (carrier bulls mated to daughters of carrier sires) was compared with non-risk matings (non-carrier bulls mated to daughters of carrier sires) using a logistic regression model. Perinatal and first-year mortality was analyzed in 32,100 and 46,566 Braunvieh and Fleckvieh calvings, respectively, with known haplotype states from sires and maternal grandsires. Calves that left the recording system within 365 days after birth for other reasons than death (e.g., export to other countries) were not considered. A Kaplan-Meier estimator was obtained by comparing the first-year survival rate of calves from risk matings with non-risk matings .
Generation of sequence data
Genomic DNA of 290 animals (149 Fleckvieh, 54 Braunvieh, 51 Holstein, 15 Original Simmental, 12 Gelbvieh, 7 Northern Finncattle, 2 Ayrshire) was prepared from semen and blood samples following standard DNA extraction protocols as previously described . Briefly, paired-end libraries were prepared using the paired-end TruSeq DNA sample prep kit (Illumina inc., San Diego, CA, USA) and sequenced using the Illumina HiSeq 2500 instrument (Illumina inc., San Diego, CA, USA). The resulting reads were processed with the Illumina BaseCaller during the sequencing step. The alignment of the reads to the University of Maryland reference sequence (UMD3.1)  was performed using the BWA mem algorithm [58, 59]. The resulting per individual SAM files were converted into BAM files with SAMtools . Duplicate reads were identified and marked with the MarkDuplicates command of Picard . Part of the sequencing data were previously generated [14, 62] and were contributed to the 1000 bull genomes project .
Variant calling and imputation
SNPs, short insertions and deletions were genotyped in 290 sequenced animals simultaneously using the multi-sample approach implemented in mpileup of SAMtools along with BCFtools . Beagle  phasing and imputation was used to improve the primary genotype calling by SAMtools. The detection of structural variants was performed in 227 sequenced animals with an average genome fold coverage higher than 10x using the Pindel software package . The functional effects of polymorphic sites were predicted based on the annotation of the UMD3.1 assembly of the bovine genome using the Variant Effect Predictor tool from Ensembl [65, 66]. The consequences of non-synonymous amino acid substitutions on protein function were predicted using PolyPhen-2  and SIFT .
Identification of BH2-associated variants
We screened 110,294 polymorphic sites within a 9.71 Mb segment on BTA19 encompassing BH2 (from 3.67 Mb to 13.38 Mb) for variants that were homozygous for the alternate allele in BH2hom, heterozygous in five BH2 carriers and homozygous for the reference allele in 48 control animals of the Braunvieh population. The genotype distributions of 52 variants that fulfilled these criteria were assessed in 236 sequenced animals representing six cattle breeds other than Braunvieh.
Validation of the rs383232842 polymorphism
Genotypes for rs383232842 were obtained in 661 Braunvieh, 3807 Fleckvieh and 503 Holstein animals using a customized KASP™ genotyping assay (LGC Genomics) (FAM/HEX/reverse primer sequence: TGTTCATGAGAACGATGCTGTTCG/CTTGTTCATGAGAACGATGCTGTTCA/TGCTTGATATTCATCAGCTTCACACAGAT).
Seven animals were hospitalized for a period of 19 to 164 days at the animal clinic. Initial examination (including weighing) was performed upon admission. Weight records and blood samples from three animals were collected once or twice a week. Red and white blood cell count was determined using an automatic hematology analyzer (ADVIA 2120i, Siemens Healthcare, Vienna, Austria). Glucose, total protein, albumin, cholesterol, NEFA, total bilirubin, creatinine, potassium, sodium, AST, GLDH, LDH were measured in blood plasma using an automatic analyzer (Cobas 6000/c501, Roche Diagnostics GmbH, Rotkreuz, Switzerland).
All hospitalized animals were euthanized at 32 to 502 days of age because they suffered from chronic disease with no prospect of improvement. Necropsy and collection of tissue samples were done immediately after euthanasia. For histology, samples from nasal mucosa, trachea, lung, heart, liver, kidney, pancreas, small and large intestine, mesenteric lymph node, intestine were fixed in 10 % buffered formalin and embedded in paraffin wax. Sections (3 μm) were stained with hematoxylin and eosin (HE) and alcian blue. For transmission electron microscopy, samples from nasal mucosa, trachea, bronchi and bronchioles were collected from four animals and fixed in 5 % glutaraldehyde (Merck, Darmstadt, Germany) in 0.1 M phosphate buffer (Sigma-Aldrich, Vienna, Austria), pH 7.2, at 4 °C for 3 h and then postfixed in 1 % osmium tetroxide (Merck, Darmstadt, Germany) in the same buffer at 4 °C for 2 h. After dehydration in an alcohol gradient series and propylene oxide (Merck, Darmstadt, Germany), the tissue samples were embedded in glycidyl ether 100 (Serva, Heidelberg, Germany). The ultrathin sections were cut on a Leica Ultramicrotome (Leica Ultracut S, Vienna, Austria) and stained with uranyl acetate (Sigma-Aldrich, Vienna, Austria) and lead citrate (Merck, Darmstadt, Germany). Ultrathin sections were examined with a Zeiss TEM 900 electron microscope (Carl Zeiss, Oberkochen, Germany) operated at 50 kV.
We acknowledge the InterGenomics partners (Arbeitsgemeinschaft der Österreichischen Braunviehzüchter, Austria; Arbeitsgemeinschaft Deutsches Braunvieh, Germany; Associazione Nazionle Allevatori Bovini della Razza Bruna Italiana, Italy; Brown Swiss Association of the US, USA; Brune Genetique Services, France; Braunvieh Schweiz, Switzerland; and Zveza rejcev govedi rjave pasme Slovenije, Slovenia), the Arbeitsgemeinschaft Süddeutscher Rinderzüchter und Besamungsorganisationen e.V. (ASR), the Arbeitsgemeinschaft österreichischer Fleckviehzüchter (AGÖF) and the Förderverein Bioökonomieforschung e.V. (FBF) for providing genotype data. We thank Associazione Nazionale Allevatori Bovini della Razza Bruna Italiana (ANARBI) and Arbeitsgemeinschaft Deutsches Braunvieh for providing pictures of affected animals and support in sample collection.
The calf monitoring project was financially supported by Braunvieh Austria and Braunvieh Schweiz. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of supporting data
Whole-genome sequencing data of BH2hom were deposited in the European Nucleotide Archive under accession number PRJEB12807 (http://www.ebi.ac.uk/ena/data/view/PRJEB12807). The sequencing data of 241 animals from the 1000 bull genomes project are publically available in the Sequence Read Archive of NCBI (http://www.ncbi.nlm.nih.gov/sra) under accession number SRP039339.
HP, HS and RF conceived the study, designed the experiments and analyzed the data; HS and CF analyzed genotype and phenotype data; HP analyzed sequence and genotype data; HS, FRS, RW and OG conceived the calf monitoring project; JB, TW, ND, MH, HW and UB examined homozygous animals at the animal clinic; CW carried out the sequencing experiments, CW and SJ carried out the molecular genetic experiments, BFW and MD contributed to the design of the experiment, HP wrote the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All homozygous animals examined were descendants from matings between carrier animals that happened inadvertently in the Braunvieh and Fleckvieh populations prior to the declaration of BH2 carrier bulls. Ear tissue samples of calves from at-risk matings were collected from trained breeding consultants. Seven diseased homozygous animals were hospitalized and examined at the Clinic for Ruminants of the Veterinary University of Vienna, Austria and at the Clinic of Reproductive Medicine, Vetsuisse-Faculty, University Zurich, Switzerland as part of their regular practice. All hospitalized animals were euthanized because of steadily declining health condition resulting from recurrent airway disease with no prospect of improvement. The decision to euthanize the animals was at the discretion of the veterinarians at the animal clinic. DNA of healthy animals for genotyping and sequencing was prepared from semen samples of artificial insemination (AI) bulls. Semen samples were collected by approved commercial AI stations as part of their regular breeding and reproduction measures in cattle industry. No ethical approval was required for this study.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Gulliksen SM, Lie KI, Løken T, Osterås O. Calf mortality in Norwegian dairy herds. J Dairy Sci. 2009;92:2782–95.View ArticlePubMedGoogle Scholar
- Silva del Río N, Stewart S, Rapnicki P, Chang YM, Fricke PM. An observational analysis of twin births, calf sex ratio, and calf mortality in Holstein dairy cattle. J Dairy Sci. 2007;90:1255–64.View ArticlePubMedGoogle Scholar
- Brickell JS, McGowan MM, Pfeiffer DU, Wathes DC. Mortality in Holstein-Friesian calves and replacement heifers, in relation to body weight and IGF-I concentration, on 19 farms in England. Animal. 2009;3:1175–82.View ArticlePubMedGoogle Scholar
- Fuerst-Waltl B, Fuerst C. Effect of inbreeding depression on survival of Austrian Brown Swiss calves and heifers. J Dairy Sci. 2012;95:6086–92.View ArticlePubMedGoogle Scholar
- Windeyer MC, Leslie KE, Godden SM, Hodgins DC, Lissemore KD, LeBlanc SJ. Factors associated with morbidity, mortality, and growth of dairy heifer calves up to 3 months of age. Prev Vet Med. 2014;113:231–40.View ArticlePubMedGoogle Scholar
- Olsen HG, Hayes BJ, Kent MP, Nome T, Svendsen M, Lien S. A genome wide association study for QTL affecting direct and maternal effects of stillbirth and dystocia in cattle. Anim Genet. 2010;41:273–80.View ArticlePubMedGoogle Scholar
- Pausch H, Flisikowski K, Jung S, Emmerling R, Edel C, Götz K-U, et al. Genome-Wide Association Study Identifies Two Major Loci Affecting Calving Ease and Growth-Related Traits in Cattle. Genetics. 2011;187:289–97.View ArticlePubMedPubMed CentralGoogle Scholar
- Hansen M, Madsen P, Jensen J, Pedersen J, Christensen LG. Genetic parameters of postnatal mortality in Danish Holstein calves. J Dairy Sci. 2003;86:1807–17.View ArticlePubMedGoogle Scholar
- Fuerst-Waltl B, Sørensen MK. Genetic analysis of calf and heifer losses in Danish Holstein. J Dairy Sci. 2010;93:5436–42.View ArticlePubMedGoogle Scholar
- Jung S, Pausch H, Langenmayer MC, Schwarzenbacher H, Majzoub-Altweck M, Gollnick NS, et al. A nonsense mutation in PLD4 is associated with a zinc deficiency-like syndrome in Fleckvieh cattle. BMC Genomics. 2014;15:623.View ArticlePubMedPubMed CentralGoogle Scholar
- Sartelet A, Druet T, Michaux C, Fasquelle C, Géron S, Tamma N, et al. A Splice Site Variant in the Bovine RNF11 Gene Compromises Growth and Regulation of the Inflammatory Response. PLoS Genet. 2012;8, e1002581.View ArticlePubMedPubMed CentralGoogle Scholar
- Hirano T, Kobayashi N, Matsuhashi T, Watanabe D, Watanabe T, Takasuga A, et al. Mapping and Exome Sequencing Identifies a Mutation in the IARS Gene as the Cause of Hereditary Perinatal Weak Calf Syndrome. PLoS One. 2013;8, e64036.View ArticlePubMedPubMed CentralGoogle Scholar
- Venhoranta H, Pausch H, Flisikowski K, Wurmser C, Taponen J, Rautala H, et al. In frame exon skipping in UBE3B is associated with developmental disorders and increased mortality in cattle. BMC Genomics. 2014;15:890.View ArticlePubMedPubMed CentralGoogle Scholar
- Pausch H, Schwarzenbacher H, Burgstaller J, Flisikowski K, Wurmser C, Jansen S, et al. Homozygous haplotype deficiency reveals deleterious mutations compromising reproductive and rearing success in cattle. BMC Genomics. 2015;16:312.View ArticlePubMedPubMed CentralGoogle Scholar
- Kipp S, Segelke D, Reinhardt F, Reents R, Schierenbeck S, Wurmser C, et al. A new Holstein Haplotype affecting calf survival. Interbull Bulletin. 2015;49:49–53.Google Scholar
- Menzi F, Besuchet-Schmutz N, Fragnière M, Hofstetter S, Jagannathan V, Mock T, et al. A transposable element insertion in APOB causes cholesterol deficiency in Holstein cattle. Anim Genet. 2016;47:253–7.View ArticlePubMedPubMed CentralGoogle Scholar
- VanRaden PM, Olson KM, Null DJ, Hutchison JL. Harmful recessive effects on fertility detected by absence of homozygous haplotypes. J Dairy Sci. 2011;94:6153–61.View ArticlePubMedGoogle Scholar
- Fritz S, Capitan A, Djari A, Rodriguez SC, Barbat A, Baur A, et al. Detection of Haplotypes Associated with Prenatal Death in Dairy Cattle and Identification of Deleterious Mutations in GART, SHBG and SLC37A2. PLoS One. 2013;8, e65550.View ArticlePubMedPubMed CentralGoogle Scholar
- Schwarzenbacher H, Fuerst C, Fuerst-Waltl B, Dolezal M. A genome-wide search for harmful recessive haplotypes in Brown Swiss and Fleckvieh cattle. In: Book of Abstracts of the 63rd European Association for Animal Production Annual Meeting Bratislava. Slovakia: EAAP; 2012.Google Scholar
- AIPL Changes to Evaluation System (August 2013). https://www.cdcb.us/reference/changes/eval1308.htm. Accessed 6 April 2016
- Daetwyler HD, Capitan A, Pausch H, Stothard P, van Binsbergen R, Brøndum RF, et al. Whole-genome sequencing of 234 bulls facilitates mapping of monogenic and complex traits in cattle. Nat Genet. 2014;46:858–65.View ArticlePubMedGoogle Scholar
- Schwarzenbacher H, Wurmser C, Flisikowski K, Misurova L, Jung S, Langenmayer MC, et al. A frameshift mutation in GON4L is associated with proportionate dwarfism in Fleckvieh cattle. Genet Sel Evol. 2016;48:25.View ArticlePubMedPubMed CentralGoogle Scholar
- Drögemüller C, Reichart U, Seuberlich T, Oevermann A, Baumgartner M, Kühni Boghenbor K, et al. An Unusual Splice Defect in the Mitofusin 2 Gene (MFN2) Is Associated with Degenerative Axonopathy in Tyrolean Grey Cattle. PLoS One. 2011;6, e18931.View ArticlePubMedPubMed CentralGoogle Scholar
- Pausch H, Kölle S, Wurmser C, Schwarzenbacher H, Emmerling R, Jansen S, et al. A Nonsense Mutation in TMEM95 Encoding a Nondescript Transmembrane Protein Causes Idiopathic Male Subfertility in Cattle. PLoS Genet. 2014;10, e1004044.View ArticlePubMedPubMed CentralGoogle Scholar
- Kadri NK, Sahana G, Charlier C, Iso-Touru T, Guldbrandtsen B, Karim L, et al. A 660-Kb Deletion with Antagonistic Effects on Fertility and Milk Production Segregates at High Frequency in Nordic Red Cattle: Additional Evidence for the Common Occurrence of Balancing Selection in Livestock. PLoS Genet. 2014;10, e1004049.View ArticlePubMedPubMed CentralGoogle Scholar
- Drögemüller C, Tetens J, Sigurdsson S, Gentile A, Testoni S, Lindblad-Toh K, et al. Identification of the Bovine Arachnomelia Mutation by Massively Parallel Sequencing Implicates Sulfite Oxidase (SUOX) in Bone Development. PLoS Genet. 2010;6, e1001079.View ArticlePubMedPubMed CentralGoogle Scholar
- Thomsen B, Nissen PH, Agerholm JS, Bendixen C. Congenital bovine spinal dysmyelination is caused by a missense mutation in the SPAST gene. Neurogenetics. 2010;11:175–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Krebs S, Medugorac I, Röther S, Strässer K, Förster M. A missense mutation in the 3-ketodihydrosphingosine reductase FVT1 as candidate causal mutation for bovine spinal muscular atrophy. Proc Natl Acad Sci U S A. 2007;104:6746–51.View ArticlePubMedPubMed CentralGoogle Scholar
- McClure MC, Bickhart D, Null D, VanRaden P, Xu L, Wiggans G, et al. Bovine Exome Sequence Analysis and Targeted SNP Genotyping of Recessive Fertility Defects BH1, HH2, and HH3 Reveal a Putative Causative Mutation in SMC2 for HH3. PLoS One. 2014;9, e92769.View ArticlePubMedPubMed CentralGoogle Scholar
- Biffani S, Dimauro C, Macciotta N, Rossoni A, Stella A, Biscarini F. Predicting haplotype carriers from SNP genotypes in Bos taurus through linear discriminant analysis. Genet Sel Evol. 2015;47:4.View ArticlePubMedPubMed CentralGoogle Scholar
- Drögemüller M, Jagannathan V, Becker D, Drögemüller C, Schelling C, Plassais J, et al. A Mutation in the FAM83G Gene in Dogs with Hereditary Footpad Hyperkeratosis (HFH). PLoS Genet. 2014;10, e1004370.View ArticlePubMedPubMed CentralGoogle Scholar
- McClure M, Kim E, Bickhart D, Null D, Cooper T, Cole J, et al. Fine Mapping for Weaver Syndrome in Brown Swiss Cattle and the Identification of 41 Concordant Mutations across NRCAM, PNPLA8 and CTTNBP2. PLoS One. 2013;8, e59251.View ArticlePubMedPubMed CentralGoogle Scholar
- Kunz E, Rothammer S, Pausch H, Schwarzenbacher H, Seefried FR, Matiasek K, et al. Confirmation of a non-synonymous SNP in PNPLA8 as a candidate causal mutation for Weaver syndrome in Brown Swiss cattle. Genet Sel Evol. 2016;48:21.View ArticlePubMedPubMed CentralGoogle Scholar
- Decker JE, McKay SD, Rolf MM, Kim J, Molina Alcalá A, Sonstegard TS, et al. Worldwide Patterns of Ancestry, Divergence, and Admixture in Domesticated Cattle. PLoS Genet. 2014;10, e1004254.View ArticlePubMedPubMed CentralGoogle Scholar
- Noone PG, Leigh MW, Sannuti A, Minnix SL, Carson JL, Hazucha M, et al. Primary Ciliary Dyskinesia. Am J Respir Crit Care Med. 2004;169:459–67.View ArticlePubMedGoogle Scholar
- Coren M, Meeks M, Morrison I, Buchdahl R, Bush A. Primary ciliary dyskinesia: age at diagnosis and symptom history. Acta Paediatr. 2002;91:667–9.View ArticlePubMedGoogle Scholar
- Leigh MW, Pittman JE, Carson JL, Ferkol TW, Dell SD, Davis SD, et al. Clinical and genetic aspects of primary ciliary dyskinesia/Kartagener syndrome. Genet Med. 2009;11:473–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Wisseman CL, Simel DL, Spock A, Shelburne JD. The prevalence of abnormal cilia in normal pediatric lungs. Arch Pathol Lab Med. 1981;105:552–5.PubMedGoogle Scholar
- Chilvers M, Rutman A, O’Callaghan C. Functional analysis of cilia and ciliated epithelial ultrastructure in healthy children and young adults. Thorax. 2003;58:333–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Ziętkiewicz E, Bukowy-Bieryłło Z, Voelkel K, Klimek B, Dmeńska H, Pogorzelski A, et al. Mutations in Radial Spoke Head Genes and Ultrastructural Cilia Defects in East-European Cohort of Primary Ciliary Dyskinesia Patients. PLoS One. 2012;7, e33667.View ArticlePubMedPubMed CentralGoogle Scholar
- Carson JL, Collier AM, Hu SS. Acquired ciliary defects in nasal epithelium of children with acute viral upper respiratory infections. N Engl J Med. 1985;312:463–8.View ArticlePubMedGoogle Scholar
- Zariwala MA, Knowles MR, Omran H. Genetic Defects in Ciliary Structure and Function. Annu Rev Physiol. 2007;69:423–50.View ArticlePubMedGoogle Scholar
- Chilvers MA, Rutman A, O’Callaghan C. Ciliary beat pattern is associated with specific ultrastructural defects in primary ciliary dyskinesia. J Allergy Clin Immunol. 2003;112:518–24.View ArticlePubMedGoogle Scholar
- Leigh MW, O’Callaghan C, Knowles MR. The Challenges of Diagnosing Primary Ciliary Dyskinesia. Proc Am Thorac Soc. 2011;8:434–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Olm MAK, Kögler JE, Macchione M, Shoemark A, Saldiva PHN, Rodrigues JC. Primary ciliary dyskinesia: evaluation using cilia beat frequency assessment via spectral analysis of digital microscopy images. J Appl Physiol. 2011;111:295–302.View ArticlePubMedPubMed CentralGoogle Scholar
- Merveille A-C, Davis EE, Becker-Heck A, Legendre M, Amirav I, Bataille G, et al. CCDC39 is required for assembly of inner dynein arms and the dynein regulatory complex and for normal ciliary motility in humans and dogs. Nat Genet. 2011;43:72–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee L, Campagna DR, Pinkus JL, Mulhern H, Wyatt TA, Sisson JH, et al. Primary ciliary dyskinesia in mice lacking the novel ciliary protein Pcdp1. Mol Cell Biol. 2008;28:949–57.View ArticlePubMedPubMed CentralGoogle Scholar
- Dutcher SK, Trabuco EC. The UNI3 Gene Is Required for Assembly of Basal Bodies of Chlamydomonas and Encodes δ-Tubulin, a New Member of the Tubulin Superfamily. Mol Biol Cell. 1998;9:1293–308.View ArticlePubMedPubMed CentralGoogle Scholar
- Inclán YF, Nogales E. Structural models for the self-assembly and microtubule interactions of gamma-, delta- and epsilon-tubulin. J Cell Sci. 2001;114:413–22.PubMedGoogle Scholar
- O’Toole ET, Giddings TH, McIntosh JR, Dutcher SK. Three-dimensional Organization of Basal Bodies from Wild-Type and δ-Tubulin Deletion Strains of Chlamydomonas reinhardtii. Mol Biol Cell. 2003;14:2999–3012.View ArticlePubMedPubMed CentralGoogle Scholar
- Garreau de Loubresse N, Ruiz F, Beisson J, Klotz C. Role of delta-tubulin and the C-tubule in assembly of Paramecium basal bodies. BMC Cell Biol. 2001;2:4.View ArticlePubMedPubMed CentralGoogle Scholar
- Heald R, Khodjakov A. Thirty years of search and capture: The complex simplicity of mitotic spindle assembly. J Cell Biol. 2015;211:1103–11.View ArticlePubMedGoogle Scholar
- Barbato A, Frischer T, Kuehni CE, Snijders D, Azevedo I, Baktai G, et al. Primary ciliary dyskinesia: a consensus statement on diagnostic and treatment approaches in children. Eur Respir J. 2009;34:1264–76.View ArticlePubMedGoogle Scholar
- Zimin AV, Delcher AL, Florea L, Kelley DR, Schatz MC, Puiu D, et al. A whole-genome assembly of the domestic cow. Bos taurus Genome Biol. 2009;10:R42.View ArticlePubMedGoogle Scholar
- Delaneau O, Zagury J-F, Marchini J. Improved whole-chromosome phasing for disease and population genetic studies. Nat Meth. 2013;10:5–6.View ArticleGoogle Scholar
- Pausch H, Wurmser C, Reinhardt F, Emmerling R, Fries R. Short communication: Validation of 4 candidate causative trait variants in 2 cattle breeds using targeted sequence imputation. J Dairy Sci. 2015;98:4162–7.PubMedGoogle Scholar
- Therneau TM, Grambsch PM. Modeling survival data: extending the Cox model. New York: Springer; 2000.View ArticleGoogle Scholar
- Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. Preprint at http://arxiv.org/pdf/1303.3997.pdf. 2013. Accessed 24 May 2016
- Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25:1754–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Picard Tools - By Broad Institute: http://broadinstitute.github.io/picard/. Accessed 21 October 2015
- Jansen S, Aigner B, Pausch H, Wysocki M, Eck S, Benet-Pagès A, et al. Assessment of the genomic variation in a cattle population by re-sequencing of key animals at low to medium coverage. BMC Genomics. 2013;14:446.View ArticlePubMedPubMed CentralGoogle Scholar
- Browning BL, Browning SR. A Unified Approach to Genotype Imputation and Haplotype-Phase Inference for Large Data Sets of Trios and Unrelated Individuals. Am J Hum Genet. 2009;84:210–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Ye K, Schulz MH, Long Q, Apweiler R, Ning Z. Pindel: a pattern growth approach to detect break points of large deletions and medium sized insertions from paired-end short reads. Bioinformatics. 2009;25:2865–71.View ArticlePubMedPubMed CentralGoogle Scholar
- McLaren W, Pritchard B, Rios D, Chen Y, Flicek P, Cunningham F. Deriving the consequences of genomic variants with the Ensembl API and SNP Effect Predictor. Bioinformatics. 2010;26:2069–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Variant Effect Predictor. http://www.ensembl.org/Tools/VEP. Accessed 13 January 2016
- Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7:248–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Ng PC, Henikoff S. SIFT: predicting amino acid changes that affect protein function. Nucl Acids Res. 2003;31:3812–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Kraft W, Dürr UM. Klinische Labordiagnostik in der Tiermedizin. 6th ed. Stuttgart: Schattauer Verlag; 2005.Google Scholar