Germline mutation within COL2A1 associated with lethal chondrodysplasia in a polled Holstein family

Genomic diversity and signatures of selection in meat and fancy rabbit breeds based on high-density marker data

Background

Domestication of the rabbit (Oryctolagus cuniculus) has led to a multi-purpose species that includes many breeds and lines with a broad phenotypic diversity, mainly for external traits (e.g. coat colours and patterns, fur structure, and morphometric traits) that are valued by fancy rabbit breeders. As a consequence of this human-driven selection, distinct signatures are expected to be present in the rabbit genome, defined as signatures of selection or selective sweeps. Here, we investigated the genome of three Italian commercial meat rabbit breeds (Italian Silver, Italian Spotted and Italian White) and 12 fancy rabbit breeds (Belgian Hare, Burgundy Fawn, Champagne d’Argent, Checkered Giant, Coloured Dwarf, Dwarf Lop, Ermine, Giant Grey, Giant White, Rex, Rhinelander and Thuringian) by using high-density single nucleotide polymorphism data. Signatures of selection were identified based on the fixation index (F_ST) statistic with different approaches, including single-breed and group-based methods, the latter comparing breeds that are grouped based on external traits (different coat colours and body sizes) and types (i.e. meat vs. fancy breeds).

Results

We identified 309 genomic regions that contained signatures of selection and that included genes that are known to affect coat colour (ASIP, MC1R and TYR), coat structure (LIPH), and body size (LCORL/NCAPG, COL11A1 and HOXD) in rabbits and that characterize the investigated breeds. Their identification proves the suitability of the applied methodologies for capturing recent selection events. Other regions included novel candidate genes that might contribute to the phenotypic variation among the analyzed breeds, including genes for pigmentation-related traits (EDNRA, EDNRB, MITF and OCA2) and body size, with a strong candidate for dwarfism in rabbit (COL2A1).

Conclusions

We report a genome-wide view of genetic loci that underlie the main phenotypic differences in the analyzed rabbit breeds, which can be useful to understand the shift from the domestication process to the development of breeds in O. cuniculus. These results enhance our knowledge about the major genetic loci involved in rabbit external traits and add novel information to understand the complexity of the genetic architecture underlying body size in mammals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00696-9.

Bayesian Inference in Y-Linked Two-Sex Branching Processes with Mutations: ABC Approach

A Y-linked two-sex branching process with mutations and blind choice of males is a suitable model for analyzing the evolution of the number of carriers of a Y-linked allele and its mutations. Such a model considers a two-sex monogamous population in which each female chooses her partner from among the male population without caring about his type (i.e., the allele he carries). In this work, we deal with the problem of estimating the main parameters of these models by developing Bayesian inference in a parametric framework. First, we consider as a sample scheme the observation of the total number of females and males up to some generation as well as the number of males of each genotype in the last generation. Subsequently, we introduce the information on the mutated males in only the last generation, obtaining in this way a second sample scheme. For both samples, we apply the Approximate Bayesian Computation (ABC) method to approximate the posterior distributions of the main parameters of the model. The accuracy of the procedure based on these samples is illustrated and discussed by way of simulated examples.

A large deletion in the COL2A1 gene expands the spectrum of pathogenic variants causing bulldog calf syndrome in cattle

Background

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.

Case presentation

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.

Conclusions

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.

A Heterozygous Mutation in the Triple Helical Region of the Alpha 1 (II) Chain of the COL2A1 Protein Causes Non-Lethal Spondyloepiphyseal Dysplasia Congenita

Objective:Heterozygous pathogenic variants in the COL2A1 gene result in several clinical features including impaired skeletal growth, ocular and otolaryngological abnormalities. Missense mutations in the triple helical region of the COL2A1 protein have been associated with lethal spondyloepiphyseal dysplasia (SED). In this study, we aimed to identify the underlying cause of a case of SED congenita (SEDC) in a 27-month-old child. Materials and Methods: A patient who was diagnosed initially with osteochondrodysplasia underwent a detailed clinical and radiological examination to obtain a conclusive diagnosis. The patient did not show any clinical features of hypochondrogenesis. Whole exome sequencing of the COL2A1 gene was carried out to identify the underlying genetic cause of the disorder. Results: Variant annotation and filtration detected a heterozygous missense mutation c.1357G>A (p.G453S) in the exon 21 of the COL2A1 gene of the proband which was confirmed by Sanger sequencing. Neither parent carried the mvariant suggesting this was a new mutation. Conclusion: The COL2A1 mutation (c.1357G>A), identified in this case, results in more mild phenotype than other missense mutations in exon 21 which are known to cause lethal hypochondrogenesis. We showed, for the first time, that a missense mutation (p.G453S) in the triple helical region of the alpha 1 (II) chain of the COL2A1 protein underlies SEDC and is not always lethal.

A recessive lethal chondrodysplasia in a miniature zebu family results from an insertion affecting the chondroitin sulfat domain of aggrecan

Background

Congenital skeletal malformations represent a heterogeneous group of disorders affecting bone and cartilage development. In cattle, particular chondrodysplastic forms have been identified in several miniature breeds. In this study, a phenotypic characterization was performed of an affected Miniature Zebu calf using computed tomography, necropsy and histopathological examinations, whole genome sequencing of the case and its parents on an Illumina NextSeq 500 in 2 × 150 bp paired-end mode and validation using Sanger sequencing and a Kompetitive Allele Specific PCR assay. Samples from the family of an affected Miniature Zebu with bulldog syndrome including parents and siblings, 42 healthy Miniature Zebu not related with members of the herd and 88 individuals from eight different taurine cattle breeds were available for validation.

Results

A bulldog-like Miniature Zebu calf showing a large bulging head, a short and compressed body and extremely short and stocky limbs was delivered after a fetotomy. Computed tomography and necropsy revealed severe craniofacial abnormalities including a shortening of the ventral nasal conchae, a cleft hard palate, rotated limbs as well as malformed and fused vertebrae and ribs. Histopathologic examination showed a disorganization of the physeal cartilage with disorderly arranged chondrocytes in columns and a multifocal closed epiphyseal plate. Whole-genome sequencing of this malformed Miniature Zebu calf, its dam and sire and subsequent comparative sequence analysis revealed a one base pair insertion (ACAN:c.5686insC) located within the cartilage development gene aggrecan (ACAN) exclusively homozygous in the affected calf and heterozygous in its parents. This variant was predicted to cause a frameshift (p.Val1898fsTer9) and thus a truncation of the chondroitin sulfate domain as well as a loss of the C-terminal globular domain of ACAN. It perfectly co-segregated with the lethal bulldog syndrome in Miniature Zebus.

Conclusions

We found a novel mutation in ACAN causing a recessive lethal chondrodysplasia in Miniature Zebu cattle. A diagnostic test for this mutation is now available for Miniature Zebu breeders preventing further cases of bulldog syndrome by targeted matings. To the authors’ best knowledge, this is the first case of a Miniature Zebu associated with an ACAN mutation.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0678-8) contains supplementary material, which is available to authorized users.