A large deletion in the COL2A1 gene expands the spectrum of pathogenic variants causing bulldog calf syndrome in cattle
Acta Veterinaria Scandinavica volume 62, Article number: 49 (2020)
Congenital bovine chondrodysplasia, also known as bulldog calf syndrome, is characterized by disproportionate growth of bones resulting in a shortened and compressed body, mainly due to reduced length of the spine and the long bones of the limbs. In addition, severe facial dysmorphisms including palatoschisis and shortening of the viscerocranium are present. Abnormalities in the gene collagen type II alpha 1 chain (COL2A1) have been associated with some cases of the bulldog calf syndrome. Until now, six pathogenic single-nucleotide variants have been found in COL2A1. Here we present a novel variant in COL2A1 of a Holstein calf and provide an overview of the phenotypic and allelic heterogeneity of the COL2A1-related bulldog calf syndrome in cattle.
The calf was aborted at gestation day 264 and showed generalized disproportionate dwarfism, with a shortened compressed body and limbs, and dysplasia of the viscerocranium; a phenotype resembling bulldog calf syndrome due to an abnormality in COL2A1. Whole-genome sequence (WGS) data was obtained and revealed a heterozygous 3513 base pair deletion encompassing 10 of the 54 coding exons of COL2A1. Polymerase chain reaction analysis and Sanger sequencing confirmed the breakpoints of the deletion and its absence in the genomes of both parents.
The pathological and genetic findings were consistent with a case of “bulldog calf syndrome”. The identified variant causing the syndrome was the result of a de novo mutation event that either occurred post-zygotically in the developing embryo or was inherited because of low-level mosaicism in one of the parents. The identified loss-of-function variant is pathogenic due to COL2A1 haploinsufficiency and represents the first structural variant causing bulldog calf syndrome in cattle. Furthermore, this case report highlights the utility of WGS-based precise diagnostics for understanding congenital disorders in cattle and the need for continued surveillance for genetic disorders in cattle.
The bulldog calf syndrome (BDS) is a congenital form of bovine chondrodysplasia affecting bones with endochondral osteogenesis. In its most severe form, this syndrome is lethal . The BDS is often exemplified by the Dexter BDS type , which is linked to abnormalities in the aggrecan (ACAN) gene . However other BDS types, which share gross morphology features with Dexter type, are associated with abnormalities in other genes and occur in different cattle breeds. Abnormalities in the collagen type II alpha 1 chain (COL2A1) gene causing BDS have been reported several times during the last 15 years (achondrogenesis/hypochondrogenesis type II in Bos taurus; OMIA 001926-9913; https://omia.org/OMIA001926/9913/). The purpose of this study was to report a variant in the COL2A1 gene leading to BDS and provide an overview of the phenotypic and allelic heterogeneity of COL2A1-related BDS.
A stillborn Holstein male calf with a body weight of 18.1 kg was aborted at gestation day 264 (normal gestation 281 days (mean)). The pregnancy was the result of insemination with semen of a purebred Holstein sire on a Holstein dam. The parents were not related within at least four generations. The calf had moderate autolysis and was frozen at − 20 °C before submission for necropsy and was examined after thawing.
The calf had generalized disproportionate dwarfism resembling a case of BDS (Fig. 1). The body appeared shortened and compact. The limbs showed bilateral symmetric shortening, which especially affected the bones proximal to the phalanges, giving the limbs a compact appearance. The phalanges were slightly rotated medially. The limbs were sawed longitudinally, which confirmed the irregular development of diaphysis and the presence of enlarged chondroid epiphyses without ossification centers (Fig. 2a). Radiological examination prior to sawing revealed normally structured phalangeal bones, but otherwise bones were only seen as irregular diaphyseal segments that could only be identified based on their location (Fig. 2b). Vertebrae had a similar appearance with enlarged chondroid epiphyses and irregular diaphyses. The head had dysplasia of the viscerocranium with shortening of the maxillary bones, palatoschisis, protrusion of the tongue and doming of the calvarium (Fig. 3). Longitudinal sawing of the head through the midline revealed that the direction of the brain axis was elevated due to the abnormally shaped neurocranium (Fig. 3a). Radiological examination highlighted the abnormally shaped bones (Fig. 3b). The thorax was narrow and of reduced volume and was mostly occupied by an enlarged malformed heart. The heart malformation consisted of bilateral ventricular dilation, muscular hypertrophy of the right ventricular wall and dilation of the pulmonary trunk. The lung was hypoplastic and nonaerated. The abdomen was dilated with eventration of intestinal segments and the liver appeared indurated. Due to the level of autolysis and freezing artefacts, histopathology was not performed.
Whole-genome sequencing using the NovaSeq 6000 (illumina) was performed at a read depth of ~ 26× using DNA extracted from skin and cartilage from the ear of the calf. The generated sequences were mapped to the ARS‐UCD1.2 reference genome, and single-nucleotide variants (SNVs) and small indel variants were called. The applied software and steps to process fastq files into binary alignment map (BAM) and genomic variant call format (GVCF) files were in accordance with the latest 1000 Bull Genomes Project processing guidelines (www.1000bullgenomes.com) . Furthermore, CombineGVCFs and CatVariants of GATK v3.8  were used to combine the GVCF files and the VariantFiltration tool of GATK was used to give the variants quality labels based on the standard GATK best practices. Lastly, functional impacts were annotated using SNPEFF v4.3  by integrating the information from the NCBI Annotation Release 106 (https://www.ncbi.nlm.nih.gov/genome/annotation_euk/Bos_taurus/106/). With the resulting GVCF, including all individual variants and their functional predictions, filtering for private variants was performed. We compared the genotypes of the calf with 494 cattle genomes of various breeds that had been sequenced in the course of other ongoing studies. The WGS data of the case can be found on ENA under the sample accession number SAMEA6528902, while a comprehensive list with all ENA accession numbers is shown in Additional file 1. A total of 20 private protein-changing single-nucleotide or short indel variants with a moderate or high predicted impact, located within 19 different genes or loci, were identified (Additional file 2). This list included no variants in COL2A1, the most likely candidate gene for the observed BDS phenotype. Therefore, Integrative Genomics Viewer (IGV)  software was used for visual inspection of the genome region containing COL2A1 on chromosome 5. A heterozygous 3513 base pair (bp) sized deletion from position 32,303,127 to 32,306,640 spanning 10 coding exons of the COL2A1 gene leading to haploinsufficiency of the encoded collagen type II alpha 1 chain protein was observed (Fig. 4). The heterozygous gross deletion variant in COL2A1 is predicted to lead to a loss-of-function of the encoded collagen type II alpha 1 chain protein and was not observed in any of the 494 cattle genomes used for comparison. Therefore, this variant was further investigated as a potentially causative variant for the observed phenotype.
To evaluate whether the deletion in COL2A1 occurred de novo, the affected genomic region was amplified by polymerase chain reaction (PCR) and Sanger sequenced using DNA of the calf and both parents. Genomic DNA was extracted from EDTA blood and semen of the dam and sire respectively, and compared to the calf’s DNA. PCR products were obtained by using primers flanking the detected COL2A1 deletion (forward: 5′-CAGGGGATGGGTCTTCCT-3′ and reverse: 5′-GCGTTAGAGAGGGAGACAGG-3′) and subsequently sequenced on an ABI3730 capillary sequencer (Thermofisher, Darmstadt, Germany). Only with the DNA of the affected calf, a PCR product of 128 bp could be amplified, whereas for both parents the amplification failed. Sanger sequencing of the obtained amplicon confirmed the previously identified breakpoints in combination with the insertion of a 10-bp segment fused in-between (chr5:g.32303127_32306640delinsTCTGGGGAGC).
Discussion and conclusions
Based on the morphology of the presented BDS case, a causative genetic variation in the COL2A1 gene was suspected. As for humans, the morphology of BDS in cattle vary widely both in the overall gross morphology and in bone morphology as exemplified in Fig. 5. It appears that cases of BDS due to abnormalities in the COL2A1 gene share a common morphology that separates them from at least some other types of bovine BDS, although few types of bovine BDS have been characterised to the molecular level. BDS cases due to abnormalities in the COL2A1 gene are delivered at term or during the last 3 weeks of gestation. The affected calves have a significantly reduced body weight with a mean of 22.3 kg (variation 16.3–27.5 kg) for 11 Holstein cases [1, 8, 9]. The body and limbs are short and compressed with the digits being almost half of normal size, but normally shaped. The long bones of the limbs and the vertebrae have small irregular diaphyses and enlarged chondroid epiphyses. The viscerocranium is dysplastic with palatoschisis, the neurocranium doomed causing dorso-caudal rotation of the brain, the heart is malformed due to the narrow-spaced thorax, the lungs compressed and the liver with signs of chronic stasis. Cases of BDS that share this morphology, should be suspected of having a defect in the COL2A1 gene; a suspicion that is helpful when analysing WGS data. However, in this case, filtering for private variants in COL2A1 did not lead to the detection of a private single-nucleotide or short indel variant. Consequently, the genome data was visually inspected for the presence of structural variants in the gene that allowed the detection of a heterozygous gross deletion. It was assumed that it had occurred either post-zygotically in the developing embryo or was inherited from a parent having low-level mosaicism. The former seems to be more likely as amplification of the mutant allele failed in the examined tissues of both parents, especially because the germline of the sire was analysed by extracting DNA from semen. This means that the COL2A1 deletion observed in heterozygous state in the affected offspring was most likely absent in the genome of both parents. Therefore, we can assume that the identified mutation arose indeed de novo in the developing embryo explaining this isolated case.
This pathogenic variant is predicted to affect a large portion of the COL2A1 gene leading to haploinsufficiency. Recent large data from human genome sequencing studies presented in the Genome Aggregation Database (gnomAD)  showed that the probability of loss-of-function intolerance score for COL2A1 was 1 meaning that COL2A1 falls into the class of loss-of-function haploinsufficient genes. Collagens are normally extracellular structural proteins involved in formation of connective tissue structure. The highly conserved sequence predominantly consists of repeated three amino acids with glycine (Gly) followed by two other amino acids (Gly-x-y, where x and y can be any amino acid) but glycine being mandatory for the tight packing of the polyproline II type helices within the triple helix  (Fig. 6). We assume that the pathogenic variant reported in this study disrupted the triple-helical region of alpha 1 (II) chain and caused a dominant-negative effect similar to most of the alterations responsible for achondrogenesis/hypochondrogenesis type II (OMIM 200610; https://www.omim.org/entry/200610) in human patients. In man, variants in COL2A1 are associated with 15 different phenotypes exclusively following dominant inheritance (OMIM 120140).
Interestingly, the OMIA 001926-9913 BDS type occurs either as de novo or inherited from a mosaic parent . Mosaic sires have been found to transmit the dominant genetic abnormality to their offspring at rates ranging from 1 to 21% [12, 13] reflecting at what fetal developmental stage the gene change occurred. As it cannot be predicted if the abnormality is occurring de novo or if it is transmitted from a parent, cases must be analyzed in detail to prevent the birth of large numbers of defective offspring, in particular if the abnormality is transmitted from breeding sire with high generic merit used for artificial breeding.
A total of six independent pathogenic dominant variants in COL2A1, considered to be responsible for BDS, have been previously identified [8, 9, 12,13,14] (Table 1). All these variants involve a single nucleotide; five out of the six reported variants represent missense variants that cause a change in a glycine residue disrupting the Gly-x–y structural motif essential for the assembly of the collagen triple-helix.
This is the first report of a large deletion in the COL2A1 gene associated with BDS. The previous reported single nucleotide variants were missense and splicing. The relevance of this case report is to show that also larger-sized genomic deletions cause a similar congenital phenotype and thereby expanding the knowledge on this condition by emphasizing that different mutations in COL2A1 cause a uniform phenotype. For many genes it is known that the kind of genetic alteration influence the phenotypic outcome, e.g. the severity of a congenital defect varies or differs totally depending on the individual variant. Interestingly for COL2A1 in cattle this seems not to be the case as different kinds of variants always cause an identical phenotype which is of importance for diagnostic pathologists. Furthermore, this report provides an overview of the phenotypic and allelic heterogeneity of the COL2A1-related BDS in cattle. This example highlights the utility of WGS-based precise diagnostics for understanding disorders linked to de novo mutations in animals with an available reference genome sequence and the need for continued surveillance for genetic disorders in cattle breeding. Genome sequencing might improve the precision of the clinicopathological diagnosis as sometimes unexpected variants in genes that were not known to be associated with a certain disorder could be detected.
Availability of data and materials
Whole-genome sequence data generated from the affected calf is available under study accession PRJEB18113 and sample accession SAMEA6528902 from the European Nucleotide Archive (ENA). In addition, further control genomes are listed in Additional file 1 and can also be accessed on ENA.
Binary alignment map
Bulldog calf syndrome
- COL2A1 :
Collagen type II alpha 1 chain
Genome Aggregation Database
Genomic variant call format
Integrative Genomics Viewer
Online Mendelian Inheritance in Animals, https://omia.org/home/
Online Mendelian Inheritance in Man, https://omim.org/
Polymerase chain reaction
Agerholm JS, Arnbjerg J, Andersen O. Familial chondrodysplasia in Holstein calves. J Vet Diagn Invest. 2004;16:293–8.
Seligmann CG. Cretinism in calves. J Pathol. 1904;9:311–22.
Cavanagh JA, Tammen I, Windsor PA, Bateman JF, Savarirayan R, Nicholas FW, et al. Bulldog dwarfism in Dexter cattle is caused by mutations in ACAN. Mamm Genome. 2007;18:808–14.
Hayes BJ, Daetwyler HD. 1000 Bull Genomes Project to map simple and complex genetic traits in cattle: applications and outcomes. Annu Rev Anim Biosci. 2019;7:89–102.
DePristo M, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43:491–8.
Cingolani P, Platts A, le Wang L, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: sNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly. 2012;6:80–92.
Robinson JT, Thorvaldsdóttir H, Wenger AM, Zehir A, Mesirov JP. Variant review with the Integrative Genomics Viewer. Cancer Res. 2017;77:31–4.
Agerholm JS, Menzi F, McEvoy FJ, Jagannathan V, Drögemüller C. Lethal chondrodysplasia in a family of Holstein cattle is associated with a de novo splice site variant of COL2A1. BMC Vet Res. 2016;12:100.
Reinartz S, Mohwinkel H, Sürie C, Hellige M, Feige K, Eikelberg D, et al. Germline mutation within COL2A1 associated with lethal chondrodysplasia in a polled Holstein family. BMC Genomics. 2017;18:762.
Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. 2020;581:434–43.
Shoulders MD, Raines RT. Collagen structure and stability. Annu Rev Biochem. 2009;78:929–58.
Bourneuf E, Otz P, Pausch H, Jagannathan V, Michot P, Grohs C, et al. Rapid discovery of de novo deleterious mutations in cattle enhances the value of livestock as model species. Sci Rep. 2017;7:11466.
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.
Häfliger IM, Behn H, Freick M, Jagannathan V, Drögemüller C. A COL2A1 de novo variant in a Holstein bulldog calf. Anim Genet. 2019;50:113–4.
Andresen E, Christensen K, Jensen PT, Venge O, Rasmussen PG, Basse A, et al. Dværgvækst hos kvæg af Rød Dansk Malkerace (Achondroplasia bovis typus RDM). Nord VetMed. 1974;26:681–91.
The veterinary practice Kvægdyrlægerne Midt Aps, Bording, Denmark is thanked for submitting the calf and the owner for donating it. The authors would like to acknowledge the Next Generation Sequencing Platform of the University of Bern for performing the whole-genome sequencing experiments and the Interfaculty Bioinformatics Unit of the University of Bern for providing high-performance computational infrastructure. The authors are grateful to Nathalie Besuchet-Schmutz for expert technical assistance.
Figure 5a has previously been published in Agerholm, JS. Inherited disorders in Danish cattle. APMIS. 2007;115 (suppl 122): 1–76.
The study was funded by the Danish Bovine Genetic Disease Programme, internal funds for the University of Copenhagen, and the Swiss National Science Foundation.
Ethics approval and consent to participate
This study did not require official or institutional ethical approval as it was not experimental.
Consent for publication
JSA is editor-in-chief of Acta Veterinaria Scandinavia, but has not in any way been involved in or interacted with the journal’s review process or editorial decision-making. The editor was blinded to the review process. The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
EBI Accession numbers of all publicly available genome sequences. We compared the genotypes of the calf with 494 cattle genomes of various breeds that had been sequenced in the course of other ongoing studies and that were publicly available.
List of the remaining private protein-coding variants after comparison of the genotypes of the calf with 494 cattle genome. A total of 20 private protein-changing single-nucleotide or short indel variants with a moderate or high predicted impact, located within 19 different genes or loci, were identified.
About this article
Cite this article
Jacinto, J.G.P., Häfliger, I.M., Letko, A. et al. A large deletion in the COL2A1 gene expands the spectrum of pathogenic variants causing bulldog calf syndrome in cattle. Acta Vet Scand 62, 49 (2020). https://doi.org/10.1186/s13028-020-00548-w