Lethal chondrodysplasia in a family of Holstein cattle is associated with a de novo splice site variant of COL2A1
© The Author(s). 2016
Received: 4 February 2016
Accepted: 7 June 2016
Published: 13 June 2016
Lethal chondrodysplasia (bulldog syndrome) is a well-known congenital syndrome in cattle and occurs sporadically in many breeds. In 2015, it was noticed that about 12 % of the offspring of the phenotypically normal Danish Holstein sire VH Cadiz Captivo showed chondrodysplasia resembling previously reported bulldog calves. Pedigree analysis of affected calves did not display obvious inbreeding to a common ancestor, suggesting the causative allele was not a rare recessive. The normal phenotype of the sire suggested a dominant inheritance with incomplete penetrance or a mosaic mutation.
Three malformed calves were examined by necropsy, histopathology, radiology, and computed tomography scanning. These calves were morphologically similar and displayed severe disproportionate dwarfism and reduced body weight. The syndrome was characterized by shortening and compression of the body due to reduced length of the spine and the long bones of the limbs. The vicerocranium had severe dysplasia and palatoschisis. The bones had small irregular diaphyses and enlarged epiphyses consisting only of chondroid tissue.
The sire and a total of four affected half-sib offspring and their dams were genotyped with the BovineHD SNP array to map the defect in the genome. Significant genetic linkage was obtained for several regions of the bovine genome including chromosome 5 where whole genome sequencing of an affected calf revealed a COL2A1 point mutation (g.32473300 G > A). This private sequence variant was predicted to affect splicing as it altered the conserved splice donor sequence GT at the 5’-end of COL2A1 intron 36, which was changed to AT. All five available cases carried the mutant allele in heterozygous state and all five dams were homozygous wild type. The sire VH Cadiz Captivo was shown to be a gonadal and somatic mosaic as assessed by the presence of the mutant allele at levels of about 5 % in peripheral blood and 15 % in semen.
The phenotypic and genetic findings are comparable to a previously reported COL2A1 missense mutation underlying lethal chondrodysplasia in the offspring of a mosaic French Holstein sire (Igale Masc). The identified independent spontaneous splice site variant in COL2A1 most likely caused chondrodysplasia and must have occurred during the early foetal development of the sire. This study provides a first example of a dominant COL2A1 splice site variant as candidate causal mutation of a severe lethal chondrodysplasia phenotype. Germline mosaicism is a relatively frequent mechanism in the origin of genetic disorders and explains the prevalence of a certain fraction of affected offspring. Paternal dominant de novo mutations are a risk in cattle breeding, especially because the ratio of defective offspring may be very high and be associated with significant animal welfare problems.
Chondrodysplasia is a developmental bone defect occurring due to disturbed endochondral osteogenesis. This leads to a reduced longitudinal growth of bones such as those of the limbs, spine and face. The most severe forms are found in lethal congenital generalized chondrodysplasia, which was originally reported in Dexter cattle in 1904  and usually referred to as “bulldog calves.” Congenital chondrodysplasia has been described in a number of cattle breeds with varying phenotype and different modes of inheritance (OMIA 000004-9913, OMIA 000187-9913, OMIA 000311-9913, OMIA 000189-9913).
Development in molecular genetics has revolutionized the research in bovine teratology as the methods to identify genetic causes have improved significantly . The molecular basis of a few bovine chondrodysplasias has been determined. In 2007, two independent mutations affecting the coding region of the aggrecan (ACAN) gene causing lethal chondrodysplasia in Dexter cattle (OMIA 001271-9913) were described. These ACAN mutations showed incomplete dominant inheritance, leading to a mild form of dwarfism in heterozygotes, while homozygous animals displayed extreme chondrodysplasia and usually died during gestation . For another lethal chondrodysplasia phenotype (OMIA 001926-9913) reported in Holstein cattle in 2004  it was later demonstrated that the sire (Igale Masc) of affected calves was mosaic for a dominant acting collagen type II (COL2A1) missense mutation (pG600D) disrupting the Gly-X-Y structural motif essential for the assembly of the collagen triple-helix . In 2015, cases of stillborn bulldog-like calves occurring among the offspring of the Danish Holstein sire VH Cadiz Captivo (DK256588) were reported to the Danish bovine genetic disease programme . This study reports the pathological investigation of the condition and the genetic analyses that led to the identification of a candidate causal mutation.
Three Holstein calves were submitted for examination: Case 1: a male delivered at gestation day (GD) 271; Case 2: a female delivered at GD 273; Case 3: a male delivered at GD 269. The cases originated from different herds. Ear tissue from two additional cases, where the herd veterinarian diagnosed the condition, was also submitted. The Holstein bull VH Cadiz Captivo was registered as the sire in all cases. Ethylenediaminetraacetic acid (EDTA) stabilized blood from all five dams was also available as was semen and EDTA stabilized blood of the sire.
Post mortem examinations
Two calves (cases 2 and 3) underwent initial full body computed tomography (CT) scans using a single slice helical CT machine (Emotion, Siemens, Erlangen, Germany) to obtain a full view of the bone malformations. Slice thickness was 3 mm and surface rendered reconstructions, which help illustrate bone malformations, were made using OsiriX software . The calves were then necropsied, which included longitudinal sectioning of the right limbs and the spine. Specimens of heart, lung, liver, kidneys, adrenal glands, thymus, vertebrae, femur and humerus were fixed in 10 % neutral buffered formalin for histology. The left limbs and the head of all calves were frozen at -20 °C and later examined by radiology. The head was sectioned longitudinally though the midline before radiology to obtain better images. The limbs were sectioned longitudinally after radiology to check for bilateral similarities in gross bone morphology. Tissues for histology were processed by routine methods, embedded in paraffin, sectioned at 2–3 μm and stained with haematoxylin and eosin. Bone specimens were decalcified in a 3.3 % formaldehyde/17 % formic acid solution before processing.
Data on the outcomes of inseminations with semen of VH Cadiz Captivo in Denmark were obtained from the Danish Cattle Database. These included offspring recorded as “defective”, “stillborn”, “dead within 24 hours” or “still alive.” A questionnaire was constructed and mailed to owners of offspring recorded as “defective”, “stillborn” or “dead within 24 hours.” The questionnaire included pictures of two chondrodysplastic calves sired by VH Cadiz Captivo, the ear tack number of the dam of the calf and the date of delivery. The owners were then contacted by phone to assess if the offspring was chondrodysplastic or not.
Data on the gestation period for all calves sired by VH Cadiz Captivo were obtained from the database and the length of the gestation period for chondrodysplastic vs. normal calves was compared using the Welch Two-sample t-test.
Genomic DNA was extracted for a total of 11 family members using standard protocols. For one of the affected calves the obtained quality of isolated DNA was very low due to tissue degradation. Therefore only a total of nine samples (four cases with their respective dams and VH Cadiz Captivo) were used for genotyping with the BovineHD BeadChip (Illumina), including 777,962 evenly distributed SNPs, at Geneseek.
MERLIN v 1.1.2 software  was used to analyze the dataset and carry out the linkage analysis. The Whittemore and Halpern non-parametric linkage pair statistics Z-mean and Kong and Cox logarithm of the odds (LOD) score were used to test for linkage using allele sharing among affected pedigree members without any assumptions on the mode of inheritance or on the genetic parameters of a specified model.
A polymerase chain reaction (PCR) free fragment library with a 390 base pair (bp) insert was prepared from one affected calf which was sequenced on half of a lane of an Illumina HiSeq3000 instrument using 2 × 150 bp paired-end reads. The genome sequencing data were deposited in the European Nucleotide Archive (ENA) (http://www.ebi.ac.uk/ena) under accession PRJEB12095. The mapping to the UMD3.1 bovine reference genome assembly and variant calling were undertaken as previously described . During the mutation analysis, 118 sequenced genomes of normal cattle from 20 genetically diverse Bos taurus breeds were used as a local control cohort. The recent sequence variant database containing 1119 already sequenced genomes of the ongoing 1000 bull genomes project  was used as a global control cohort during filtering for private variants of the sequenced affected calf.
Longitudinal sectioning of the spine and long bones of the limbs showed small irregular diaphyses of increased hardness and significantly enlarged epiphyses consisting entirely of a homogenous chondroid rubber-like tissue. Distinct epiphyseal lines could not be identified as the irregular metaphyses opposed the chondroid epiphyses that lacked a center of ossification. The dorsally enlarged epiphyses of the vertebrae caused compression of the spinal cord.
The thorax and abdomen were of reduced volume and the enclosed organs were consequently closely opposed. The heart had biventricular myocardial hypertrophia and dilation of the right ventricle. Slight hydrothorax was present. The lungs had diffuse congenital atelectasis. The liver surface was irregular and the texture of the parenchyma increased. Males had bilateral abdominal cryptorchidism.
Analysis of breeding data showed that HV Cadiz Captivo had 521 recorded offspring of which 412 were born alive and had received an ear tack number, i.e. regarded as normal. The owners of the remaining 109 calves were contacted by questionnaire and phone. Sixty-four owners reported that the offspring had been a “bulldog calf” while 39 reported that the calf was of normal proportions. Six calves were excluded as the owner did not respond or could not recall the morphology. HV Cadiz Captivo therefore produced “bulldog calves” at a ratio of 12.4 % (64/515). Pedigree analysis of affected calves did not display obvious inbreeding to a common ancestor.
The length of the gestation period between cows having normal offspring vs.“bulldog calves” was compared. The mean among the cows with a normal calving was 280.2 days, which was significantly different (P < 0.0001) from the mean among cows giving birth to “bulldog calves” (276.2 days).
The gross, microscopic and radiographic lesions in the examined offspring of the sire VH Cadiz Captivo are consistent with generalized congenital chondrodysplasia, a condition usually referred to as “bulldog calf syndrome.” Pedigree analysis of affected calves did not display obvious inbreeding to a common ancestor, suggesting the causative allele is not a rare recessive. The normal phenotypes of the sire and the dams suggest a dominant inheritance with incomplete penetrance or a mosaic mutation. Due to the lethal effect of the mutation we concentrated on variants that were located within the coding sequences or within the splice sites of the candidate genes in the identified linked region of the bovine genome. Finally, mutation analysis for chondrodysplasia detected a perfect association between the splice site mutation of COL2A1 and the disease phenotype. The COL2A1 variant (c.2463 + 1G > A) is predicted to affect splicing because it alters the conserved splice donor sequence GT at the 5’-end of intron 36, which was changed to AT. The predicted consequences are either skipping of exon 36 or retention of intron 36. Both possible scenarios will lead to aberrant splice variants encoding truncated proteins, but nonsense-mediated decay might be the most likely major consequence of the splice site mutation. Unfortunately, the RNA was too degraded by the time the calves were available for examination so reverse transcription PCR failed to verify the consequence of the sequence variant experimentally.
Heterozygous mutations of COL2A1 (OMIM 120140) have been identified in human patients with various rare autosomal dominant conditions characterized by skeletal dysplasia, short stature, and sensorial defects collectively termed as type II collagenopathies [10, 11]. These disorders not only impair skeletal growth but also cause ocular and otolaryngological abnormalities. The classical phenotypes include the spondyloepiphyseal dysplasia (SED) spectrum with variable severity from lethal SED, including achondrogenesis type II and hypochondrogenesis (OMIM 200610), through congenital SED (OMIM 200160) to late-onset SED (OMIM 271700). One third of human COL2A1 defects are dominant-negative mutations in the triple-helical region of alpha 1 (II) chains which disrupt the collagen triple helix . These predominantly glycine to nonserine residue substitutions exclusively create more severe phenotypes like achondrogenesis type II and hypochondrogenesis and correspond to the bovine COL2A1 allele causing chondrodysplasia in Holstein cattle . In humans, splice-site mutations are assumed to cause haploinsufficiency through nonsense-mediated mRNA decay . The COL2A1 mutations resulting in a premature stop codon are found in less severe phenotypes such as Stickler dysplasia type I (OMIM 108300) characterized by ocular, auditory, skeletal, and orofacial abnormalities or Kniest dysplasia (OMIM 156550) characterized with short stature, restricted joint mobility, and blindness [10, 11]. Only in two cases, splice site mutations, similar to the mutation identified in this study, are shown to produce phenotypes at the lethal end of the type II collagenopathy spectrum [12, 13]. Therefore, it seems to be highly likely that the detected bovine COL2A1 splice site variant explains the phenotype in the affected calves. Probably a functional haploinsufficiency of collagen II, resulting from degradation of mutant transcripts due to nonsense-mediated mRNA decay, leads to a quantitative deficit of structurally normal collagen II causing the major malformations observed in the affected calves.
This lethal syndrome was present in 12.4 % of VH Cadiz Captivo offspring and was associated with the COL2A1 splice site variant being present in his germ cells, presumable in around 15 % corresponding to the prevalence of malformed offspring. Therefore this sire is confirmed as a founder mosaic as shown before for “Solid Gold” for ovine callipyge  and “Campus” for porcine myopathy . On the other hand, the genetic analysis of blood from the sire showed that the mutant COL2A1 allele was present in peripheral blood cells at a low level of about 5 %. These results reveal that VH Cadiz Captivo is a gonadal and somatic mosaic. Germline mosaicism is a relatively frequent mechanism in the origin of genetic disorders . Depending on various factors, such as the gene involved and/or the degree of mosaicism, the carrier of somatic and germline mosaicism may be asymptomatic or may present with various symptoms. There are two possibilities for the existence of such a mosaicism: one is that the mutation occurs in a germ cell that continues to divide. The other possibility is that the mutation occurs very early in a somatic cell before the separation to germinal cells and is therefore present both in somatic and germinal cells . This seems to be the most likely explanation for VH Cadiz Captivo as he carried the mutant allele in both somatic and germinal cells, although he appeared normal.
A similar chondrodysplastic phenotype was observed in the offspring of the Holstein sire Igale Masc but affected only 1–2 % of his offspring [4, 5], which indicates that the mutation probably occurred later in the development of the testicular germ cells than in VH Cadiz Captivo. These findings show that the proportion of malformed calves may vary considerably between offspring of sires carrying germ cell mutations.
Paternal dominant germ cell mutations is a risk in cattle breeding, especially because the ratio of defective offspring may be very high. Presence of a dominant mutation in a certain proportion of spermatozoa has been hypothesized in cases of osteogenesis imperfecta. Denholm and Cole  found this lethal bone disorder in 44 % of a clinically normal Friesian bull mated to unrelated females, while Agerholm et al.  diagnosed osteogenesis imperfecta in 9 % of the offspring of a normal Holstein sire. Paternal dominant germ cell mutations are a challenge to cattle breeding as diseased progeny may occur without inbreeding, which in other ways is the most significant challenge regarding bovine genetic diseases. Reporting systems and disease programs that allow proper diagnostic of defective calves and targeted genetic investigations are highly valuable to limit the occurrence of severe disorders as lethal chondrodysplasia and osteogenesis imperfecta. The sire VH Cadiz Captivo was omitted for breeding by the breeding association as soon as the condition was diagnosed in his offspring.
It is important to recognize that even within a single breed, cases of lethal generalized chondrodysplasia show allelic heterogeneity. The cases reported here had a morphology indistinguishable from the Igale Masc bulldog calves  even by detailed post mortem examinations. A mutation in the same gene was therefore possible to be the cause. However, lethal congenital chondrodysplasia in Holsteins has been reported previously. Although such cases may share the basic lesion, i.e. chondrodysplasia, during examination in the field, significant variation in the phenotype may be found during post mortem examination, e.g. the syndrome reported by Berger and Innes  had severe hydrocephalus, which distinguishes it from the COL2A1 associated “Igale Masc” and “VH Cadiz Captivo” types. Investigations and reports detailing characteristics of chondrodysplasia phenotypes are needed.
Genomic analyses identified a causative spontaneous dominant acting candidate causal mutation in COL2A1. This study provides a first example of a dominant COL2A1 splice site variant associated with a severe lethal chondrodysplasia phenotype. This mutation must have occurred during the early development of the asymptomatic sire VH Cadiz Captivo. Germline mosaicism as a relatively frequent mechanism in the origin of genetic disorders explains the prevalence of a certain fraction of affected offspring. Paternal dominant de novo mutations are a risk in cattle breeding, especially because the ratio of defective offspring may be very high and be associated with significant animal welfare problems. It is therefore important to have surveillance programmes for congenital syndromes in cattle.
ACAN, aggrecan gene; bp, base pair; CT, computed tomography; EDTA, ethylenediaminetraacetic acid; ENA, European Nucleotide Archive; GD, gestation day; LOD, logarithm of the odds; PCR, polymerase chain reaction; SED, spondyloepiphyseal dysplasia
The breeders and their veterinarians are acknowledged for referring the cases. Ms L Pedersen and Ms AB Brandt are thanked for performing the telephone interviews. Mrs H Holm is thanked for preparation of tissues sections and Mrs CS Due and Mrs MK Carlsen for technical assistance in radiology. We acknowledge the Next Generation Sequencing Platform of the University of Bern for performing the sequencing analyses. The study was supported by Seges Cattle, Aarhus, Denmark through funding of the Danish Bovine Genetic Disease Programme.
Availability of data and materials
The genome sequencing data were deposited in the European Nucleotide Archive (ENA) (http://www.ebi.ac.uk/ena) under accession PRJEB12095.
JSA was study leader, examined the defective calves, designed the questionnaire and analysed the breeding data. FJM did the CT scanning and radiology. FM performed linkage analysis, did the COL2A1 genotyping and prepared the genetic figures. VJ did the bioinformatics of the whole genome sequencing experiment. CD supervised the genetic investigations. All authors participated in writing the manuscript and have read and approved the final version.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study did not require official or institutional ethical approval. The animals were handled according to good ethical standards and Danish legislation. The owners transferred the calves through their local veterinarian to the University of Copenhagen for research purposes aiming at identifying the cause of the congenital disorder.
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.
- Seligmann MB. Cretinism in calves. J Pathol Bacteriol. 1904;9:311–22.View ArticleGoogle Scholar
- Nicholas FW, Hobbs M. Mutation discovery for Mendelian traits in non-laboratory animals: a review of achievements up to 2012. Anim Genet. 2014;45:157–70.View ArticlePubMedGoogle Scholar
- Cavanagh JA, Tammen I, Windsor PA, Bateman JF, Savarirayan R, Nicholas FW, Raadsma HW. Bulldog dwarfism in Dexter cattle is caused by mutations in ACAN. Mamm Genome. 2007;18:808–14.View ArticlePubMedGoogle Scholar
- Agerholm JS, Arnbjerg J, Andersen O. Familial chondrodysplasia in Holstein calves. J Vet Diagn Invest. 2004;16:293–8.View ArticlePubMedGoogle Scholar
- 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. doi:10.1038/ng.3034.View ArticlePubMedGoogle Scholar
- Agerholm JS, Basse A, Christensen K. Investigations on the occurrence of hereditary diseases in the Danish cattle population 1989-1991. Acta Vet Scand. 1993;34:245–53.PubMedGoogle Scholar
- Rosset A, Spadola L, Ratib O. OsiriX: An open-source software for navigating in multidimensional DICOM images. J Digit Imaging. 2004;17:205–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Abecasis GR, Cherny SS, Cookson WO, Cardon LR. Merlin-rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2002;30:97–101.View ArticlePubMedGoogle Scholar
- Murgiano L, Shirokova V, Welle MM, Jagannathan V, Plattet P, Oevermann A, et al. Hairless streaks in cattle implicate TSR2 in early hair follicle formation. PLoS Genet. 2015;11:e1005427.View ArticlePubMedPubMed CentralGoogle Scholar
- Nishimura G, Haga N, Kitoh H, Tanaka Y, Sonoda T, Kitamura M, et al. The phenotypic spectrum of COL2A1 mutations. Hum Mutat. 2005;26:36–43.View ArticlePubMedGoogle Scholar
- Barat-Houari M, Sarrabay G, Gatinois V, Fabre A, Dumont B, Genevieve D, et al. Mutation update for COL2A1 gene variants associated with type II collagenopathies. Hum Mutat. 2016;37:7–15.View ArticlePubMedGoogle Scholar
- Körkkö J, Cohn DH, Ala-Kokko L, Krakow D, Prockop DJ. Widely distributed mutations in the COL2A1 gene produce achondrogenesis type II/hypochondrogenesis. Am J Med Genet. 2000;92:95–100.View ArticlePubMedGoogle Scholar
- Mortier GR, Weis M, Nuytinck L, King LM, Wilkin DJ, De Paepe A, et al. Report of five novel and one recurrent COL2A1 mutations with analysis of genotype-phenotype correlation in patients with a lethal type II collagen disorder. J Med Genet. 2000;37:263–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Smit M, Segers K, Carrascosa LG, Shay T, Baraldi F, Gyapay G, et al. Mosaicism of Solid Gold supports the causality of a noncoding A-to-G transition in the determinism of the callipyge phenotype. Genetics. 2003;163:453–6.PubMedPubMed CentralGoogle Scholar
- Murgiano L, Tammen I, Harlizius B, Drögemüller C. A de novo germline mutation in MYH7 causes a progressive dominant myopathy in pigs. BMC Genet. 2012;13:99.View ArticlePubMedPubMed CentralGoogle Scholar
- Zlotogora J. Germ line mosaicism. Hum Genet. 1998;102:381–6.View ArticlePubMedGoogle Scholar
- Denholm LJ, Cole WG. Heritable bone fragility, joint laxity and dysplastic dentin in Friesian calves: a bovine syndrome of osteogenesis imperfecta. Aust Vet J. 1983;60:9–17.View ArticlePubMedGoogle Scholar
- Agerholm JS, Lund AM, Bloch B, Reibel J, Basse A, Arnbjerg J. Osteogenesis imperfecta in Holstein-Friesian calves. Zentralbl Veterinarmed A. 1994;41:128–38.View ArticlePubMedGoogle Scholar
- Berger J, Innes JR. Bull-dog calves (chondrodystrophy, achondroplasia) in a Friesian herd. Vet Rec. 1948;60:57.PubMedGoogle Scholar