Multiple independent de novo mutations are associated with the development of schistosoma reflexum, a lethal syndrome in cattle

Schistosoma reflexum (SR) is a lethal congenital syndrome characterized by U-shaped dorsal retroflexion of the spine and exposure of abdominal viscera. SR is usually associated with severe dystocia. The syndrome is thought to be inherited as a Mendelian trait. We collected a series of 23 SR-affected calves from four breeds (20 Holstein, one Red Danish, one Limousin, one Romagnola) and performed whole-genome sequencing (WGS). WGS was performed on 51 cattle, including 14 cases with parents (trio-based; Group 1) and nine single cases (solo-based; Group 2). Sequencing-based genome-wide association studies with 20 Holstein cases and 154 controls showed no association (above Bonferroni threshold; P-value<3 ×10-09). Assuming a monogenic recessive inheritance, no region of shared homozygosity was observed, suggesting heterogeneity. Alternatively, the presence of possible dominant acting de novo mutations were assessed. In Group 1, heterozygous private variants, absent in both parents, were found in seven cases. These involved the ACTL6A, FLNA, GLG1, IQSEC2, MAST3, MBTPS2, and MLLT1 genes. In addition, heterozygous private variants affecting the genes DYNC1LI1, PPP2R2B, SCAF8, SUGP1, and UBP1 were identified in five cases from Group 2. The detected frameshift and missense variants are predicted to cause haploinsufficiency. Each of these 12 affected genes belong to the class of haploinsufficient loss-of-function genes or are involved in embryonic and pre-weaning lethality or are known to be associated with severe malformation syndromes in humans and/or mice. This study presents for the first time a detailed genomic evaluation of bovine SR, suggesting that independent de novo mutations may explain the sporadic occurrence of SR in cattle.

Identification of genomic regions associated with differences in fleece type in Huacaya and Suri alpacas (Vicugna pacos)

The difference in fleece type is the distinguishing trait between the two types of alpacas (Vicugna pacos), Huacaya and Suri. The Suri fleece type has been found to be inherited dominantly over the Huacaya type, resulting in offspring with the Suri phenotype. The aim of our study was to map genomic regions associated with the two different fleece types. In this study, 91 alpacas (54 Huacayas and 37 Suris) from Germany and Switzerland were genotyped using the 76k alpaca SNP array. Only 59k chromosome-localised markers map to the alpaca reference assembly VicPac3.1, and after quality control 49 866 SNPs, were retained for population structure assessment and to conduct a genome-wide association study. Both principal component and neighbour-joining tree analysis showed that the two fleece-type cohorts overlapped rather than forming two distinct clusters. Genome-wide significantly associated markers were observed in the scaffold region of chromosome 16 (NW_021964192.1), which contains a cluster of keratin genes. A haplotype predominantly found in Suri alpacas has been identified which supports dominant inheritance. Variant filtering of nine whole-genome sequenced alpacas from both fleece types in the critical interval of 0.4 Mb did not reveal perfect segregation of either fleece type for specific variants. To our knowledge, this is the first study to use the recently developed species-specific SNP array to identify genomic regions associated with differences in fleece type in alpacas. There are still some limitations, such as the preliminary status of the reference assembly and the incomplete annotation of the alpaca genome.

MFSD2A frameshift variant in Kerry Hill sheep with microcephaly

Microcephaly is a rare neurodevelopmental disorder characterized by reduced skull circumference and brain volume that occurs sporadically in farm animals. We investigated an early-onset neurodegenerative disorder observed in seven lambs of purebred Kerry Hill sheep. Clinical signs included inability to stand or severe ataxia, convulsions, and early death. Diagnostic imaging and brain necropsy confirmed microcephaly. The pedigree of the lambs suggested monogenic autosomal recessive inheritance. We sequenced the genome of one affected lamb, and comparison with 115 control genomes revealed a single private protein-changing variant. This frameshift variant, MFSD2A: c.285dupA, p.(Asp96fs*9), represents a 1-bp duplication predicted to truncate 80% of the open reading frame. MFSD2A is a transmembrane protein that is essential for maintaining blood-brain barrier homeostasis and plays a key role in regulating brain lipogenesis. Human MFSD2A pathogenic variants are associated with a neurodevelopmental disorder with progressive microcephaly, spasticity, and brain imaging abnormalities (NEDMISBA, OMIM 616486). Here we present evidence for the occurrence of a recessively inherited form of microcephaly in sheep due to a loss-of-function variant in MFSD2A (OMIA 002371-9940). To the best of our knowledge, this is the first report of a spontaneous MFSD2A variant in domestic animals.

WNT10B: A locus increasing risk of brachygnathia inferior in Brown Swiss cattle

Shortening of the mandible (brachygnathia inferior) is a congenital, often inherited and variably expressed craniofacial anomaly in domestic animals including cattle. Brachygnathia inferior can lead to poorer animal health and welfare and reduced growth, which ultimately affects productivity. Within the course of the systematic conformation scoring, cases with a frequency of about 0.1% were observed in the Brown Swiss cattle population of Switzerland. In contrast, this anomaly is almost unknown in the Original Braunvieh population, representing the breed of origin. Because none of the individually examined 46 living offspring of our study cohort of 145 affected cows showed the trait, we can most likely exclude a monogenic-dominant mode of inheritance. We hypothesized that either a monogenic recessive or a complex mode of inheritance was underlying. Through a genome-wide association study of 145 cases and 509 controls with imputed 624k SNP data, we identified a 4.5 Mb genomic region on bovine chromosome 5 significantly associated with this anomaly. This locus was fine-mapped using whole-genome sequencing data. A run of homozygosity analysis revealed a critical interval of 430 kb. A breed specific frameshift duplication in WNT10B (rs525007739; c.910dupC; p.Arg304ProfsTer14) located in this genomic region was found to be associated with a 21.5-fold increased risk of brachygnathia inferior in homozygous carriers. Consequently, we present for the first time a genetic locus associated with this well-known anomaly in cattle, which allows DNA-based selection of Brown Swiss animals at decreased risk for mandibular shortening. In addition, this study represents the first large animal model of a WNT10B-related inherited developmental disorder in a mammalian species.

MOCOS-associated renal syndrome in a Brown Swiss cattle

BACKGROUND

A recessive form of MOCOS-associated xanthinuria type II is described in Tyrolean grey cattle. A similar case was identified in a 5-month-old Brown Swiss calf with hoof overgrowth, rough coat, urine sediment, and pneumonia.

HYPOTHESIS/OBJECTIVES

To characterize the disease phenotype, to evaluate its genetic etiology, and to determine the prevalence of the deleterious allele in the Brown Swiss population.

ANIMALS

An affected calf, its parents, and 65 441 Swiss dairy cattle.

METHODS

The affected animal was clinically examined and necropsied. Microarray genotyping was used to determine the genotypes and to assess the frequency of the MOCOS allele in a Brown Swiss control cohort.

RESULTS

Ultrasonography revealed hyperechoic renal pyramids with multifocal distal shadowing and echogenic sediment in the urinary bladder. Necropsy revealed suppurative bronchopneumonia and urolithiasis. Histology revealed numerous nephroliths with multifocal chronic lymphohistiocytic interstitial infiltrates, fibrosis, tubular degeneration, chronic multifocal glomerulonephritis with sclerosis, and bilateral hydronephrosis. Dysplastic changes were observed in the corium of the claw and the cornea. Genetic testing identified the homozygous presence of a known MOCOS frameshift variant in the case. Both parents were heterozygous and the prevalence of carriers in genotyped Brown Swiss cattle was 1.4% (342/24337).

CONCLUSIONS AND CLINICAL IMPORTANCE

The findings were consistent with the diagnosis of a recessive renal syndrome similar to xanthinuria type II described in Tyrolean grey cattle. The prevalence of the deleterious MOCOS allele is low in the Brown Swiss breed. However, mating of carriers should be avoided to prevent further losses.

[DDB2-associated incidence of squamous cell carcinoma in Haflingers: risk minimization by genotyping]

INTRODUCTION

SCC (squamous cell carinomas) are among the most common eye neoplasms in horses. In recent studies Haflinger horses with a homozygous genotype for a missense variant in the DDB2 gene (damage specific DNA binding protein 2) had a significant increased risk of developing ocular SCC. The aims of this study were to determine the frequency of the SCC-associated risk allele in the DDB2 gene in Swiss and Austrian Haflinger populations and to validate the previously described phenotypic correlation. For this purpose, Haflingers presented at various horse clinics in Switzerland (n = 21, including 11 SCC cases), privately kept Haflingers (n = 52, including 1 SCC case), and Haflingers from a stud farm in the Austrian Tyrol (n = 53) were recruited. The individual DDB2 genotype of the animals was determined using a polymerase chain ceaction (PCR) test using hair follicle or whole blood samples. Of the 12 horses suffering from SCC, nine had ocular SCC and three had non-ocular SCC. Six of the nine Haflingers with ocular SCC and one of the three Haflingers with non-ocular SCC were homozygous for the DDB2 variant. Of the 113 clinically normal animals, 7/113 were homozygous (6 %) and 42/113 were heterozygous (37 %), which corresponds to an allele frequency of 24,8 % in the control cohort. The risk of ocular SCC occurring in Haflingers is significantly increased with the homozygous DDB2 genotype. However, not all animals with SCC carry this gene variant and not all DDB2 homozygous animals develop SCC, which can be explained by the multifactorial genesis of the disease. Due to the high frequency of the undesirable allele, we recommend taking the individual DDB2 genotype of breeding animals into account in order to avoid homozygous offspring with a greatly increased SCC risk by excluding high-risk matings.

A missense variant in DGKG as a recessive functional variant for hepatic fibrinogen storage disease in Wagyu cattle

Hepatic fibrinogen storage disease (HFSD) was diagnosed in a 5-month-old Wagyu calf with a history of recurrent respiratory disease. It was characterized by lethargy, dehydration, acidemia, and increased liver enzyme activities. Histologically, disseminated hepatocytes were swollen and showed a single, sharply demarcated, faintly eosinophilic cytoplasmic inclusion with a ground-glass appearance, with the nucleus in an eccentric position. Cytoplasmic inclusions did not stain with the periodic acid-Schiff (PAS) reaction. Using a rabbit polyclonal antibody against fibrinogen, the cytoplasmic vacuoles in the hepatocytes stained intensely. Electron microscopy disclosed hepatocytes with membrane-bound cytoplasmic inclusions filled with fine granular material interspersed with a few coarse-grained electron-dense granules. A trio whole-genome sequencing approach identified a deleterious homozygous missense variant in DGKG (p.Thr721Ile). The allele frequency in 209 genotyped Wagyu was 7.2%. This is a report of a DGKG-related recessive inherited disorder in cattle and adds DGKG to the list of candidate genes for HFSD in other species.

Increased heterozygosity in low-pass sequencing data allows identification of blood chimeras in cattle

In about 90% of multiple pregnancies in cattle, shared blood circulation between fetuses leads to genetic chimerism in peripheral blood and can reduce reproductive performance in heterosexual co-twins. However, the early detection of heterosexual chimeras requires specialized tests. Here, we used low-pass sequencing data with a median coverage of 0.64× generated from blood samples of 322 F1 crosses between beef and dairy cattle and identified 20 putative blood chimeras through increased levels of genome-wide heterozygosity. In contrast, for 77 samples with routine SNP microarray data generated from hair bulbs of the same F1s, we found no evidence of chimerism, simultaneously observing high levels of genotype discordance with sequencing data. Fifteen out of 18 reported twins showed signs of blood chimerism, in line with previous reports, whereas the presence of five alleged singletons with strong signs of chimerism suggests that the in-utero death rate of co-twins is at the upper limit of former estimates. Together, our results show that low-pass sequencing data allow reliable screening for blood chimeras. They further affirm that blood is not recommended as a source of DNA for the detection of germline variants.

[Ocular and non-ocular squamous cell carcinomas in the Haflinger: Eight cases at the Institut Suisse de Médecine Equine (2015-2022)]

INTRODUCTION

This case series describes the clinical course of ocular and non-ocular squamous cell carinoma (SCC) in the Haflinger horse and is intended to raise awareness of the high recurrence rate and tendency to metastasize. Eight Haflingers with histologically confirmed SCC were included, five ocular and three non-ocular, who were presented at the Institut Suisse de Médecine Équine (ISME) Bern between July 2015 and January 2022. The ocular SCC cases were all presented because of an apparent mass, which in most cases was post-treatment recurrence. The occurrence of recurrences was observed between 3 weeks and 16 years after initial therapy. Four of five Haflingers with ocular SCC had an enucleation, three of which were clinically normal at the time of the completion of this study, one case was euthanized due to confirmed metastases and one due to lameness. The result of enucleations for therapy of ocular SCC was good if no metastases occurred. Of the three non-ocular SCC cases, only one case, a penile SCC, had an apparent mass. Therapy was initiated in this case, while the other two cases were euthanized shortly after diagnosis due to the poor prognosis of SCC in the appropriate locations (maxillary sinus, mandible). Metastases occurred three and two years after removal of the primary tumor in ocular SCC in the scapula, liver and lungs and in non-ocular SCC from the penis to the nostrils. Since a postmortem pathological examination was not carried out on all Haflingers, further metastases cannot be ruled out. Haflingers with SCC should be monitored by a veterinarian over the long term, as recurrences and/or metastases can still occur years later.

DYRK1B haploinsufficiency in a Holstein cattle with epilepsy

In this study, epilepsy with focal seizures progressing to generalized seizures was diagnosed in a 6-month-old Holstein heifer. The seizures were characterized by a brief pre-ictal phase with depression and vocalization. During the ictal phase eyelid spasms, tongue contractions, nodding and abundant salivation were observed, rapidly followed by a convulsive phase with bilateral tonic, clonic or tonic-clonic activity and loss of consciousness. Finally, during the postictal phase the heifer was obtunded and disorientated, unable to perceive obstacles and hypermetric, and pressed its head against objects. In the inter-seizure phase, the heifer was clinically normal. Neuropathology revealed axonal degeneration in the brainstem and diffuse astrocytic hypertrophic gliosis. Whole genome sequencing of the affected heifer identified a private heterozygous splice-site variant in DYRK1B (NM_001081515.1: c.-101-1G>A), most likely resulting in haploinsufficiency owing to loss-of-function. This represents a report of a DYRK1B-associated disease in cattle and adds DYRK1B to the candidate genes for epilepsy.

Loss-of-function variant in the ovine TMCO6 gene in North Country Cheviot sheep with motor neuron disease

In North Country Cheviot lambs with early-onset progressive ataxia and motor neuron degeneration, whole-genome sequencing identified a homozygous loss-of-function variant in the ovine transmembrane and coiled-coil domains (TMCO6) gene. The familial recessive form of motor neuron disease in sheep is due to a pathogenic 4 bp deletion leading to a 50% protein truncation that is assumed to result in the absence of a functional TMCO6. This uncharacterised protein is proposed to interact with ubiquilin 1 which is associated with Alzheimer's disease, whereas sporadic forms of amyotrophic lateral sclerosis are caused by variants in UBQLN2. Our findings provide a first spontaneous animal model for TMCO6, which could have implications in the studies of other comparative neurodegenerative diseases. In addition, these results will allow the design of a genetic test to prevent the occurrence of this fatal disease in the affected sheep population.

SGCD Missense Variant in a Lagotto Romagnolo Dog with Autosomal Recessively Inherited Limb-Girdle Muscular Dystrophy

An 8-month-old female Lagotto Romagnolo dog was presented for a 1-month history of an initial severe reluctance to move, rapidly progressing to a marked stiff gait and progressive muscular weakness and evolving to tetraparesis, which persuaded the owner to request euthanasia. A primary muscle pathology was supported by necropsy and histopathological findings. Macroscopically, the muscles were moderately atrophic, except for the diaphragm and the neck muscles, which were markedly thickened. Histologically, all the skeletal muscles examined showed atrophy, hypertrophy, necrosis with calcification of the fibers, and mild fibrosis and inflammation. On immunohistochemistry, all three dystrophin domains and sarcoglycan proteins were absent. On Western blot analysis, no band was present for delta sarcoglycan. We sequenced the genome of the affected dog and compared the data to more than 900 control genomes of different dog breeds. Genetic analysis revealed a homozygous private protein-changing variant in the SGCD gene encoding delta- sarcoglycan in the affected dog. The variant was predicted to induce a SGCD:p.(Leu242Pro) change in the protein. In silico tools predicted the change to be deleterious. Other 770 Lagotto Romagnolo dogs were genotyped for the variant and all found to be homozygous wild type. Based on current knowledge of gene function in other mammalian species, including humans, hamsters, and dogs, we propose the SGCD missense variant as the causative variant of the observed form of muscular dystrophy in the index case. The absence of the variant allele in the Lagotto Romagnolo breeding population indicates a rare allele that has appeared recently.

Organoids of patient-derived medullary thyroid carcinoma: The first milestone towards a new in vitro model in dogs

Organoid cultures could constitute a valuable in vitro model to explore new treatments for canine (c) medullary thyroid carcinoma (MTC). The study's objectives were to establish and characterize 3D organoid cultures of cMTC using histology and immunohistochemistry (IHC) and to evaluate the effect of antitumor drugs on organoids' viability. Five cMTC tissue samples were used to develop organoid cultures of which one organoid line, named cMTC N°2, could be passaged for an extended period. This cMTC N°2 organoid line was further compared to the primary tumour regarding morphology and IHC expression of thyroid transcription factor-1 (TTF-1), thyroglobulin, calcitonin, synaptophysin, vimentin, Ki-67, cyclooxygenase-2 (COX-2), P-glycoprotein and vascular endothelial growth factor (VEGF). Quality control of the cMTC N°2 organoid line was achieved by a single nucleotide polymorphism (SNP) array of the organoids, primary tumour and healthy blood cells of the same dog. The effect of carboplatin, meloxicam and toceranib phosphate (TOC) on cMTC N°2 organoids' viability was evaluated. The cMTC N°2 organoid line was cultured for 94 days and showed similar histological features with the primary tumour. Immunolabelling for TTF-1, thyroglobulin, calcitonin and VEGF was similar between the primary tumour and cMTC N°2 organoids. Compared to the primary tumour, organoids showed higher immunolabelling for vimentin and Ki-67, and lower immunolabelling for synaptophysin, COX-2 and P-glycoprotein. The SNP genotype was similar for each chromosome between healthy blood cells, primary tumour and cMTC N°2 organoids. Carboplatin, meloxicam and TOC had no effect on cMTC N°2 organoid cell viability within achievable in vivo concentration range. In conclusion, the cMTC N°2 organoid line is a promising first milestone towards an established in vitro organoid model to explore pathophysiology and new treatment modalities in cMTC.

Genomic regions underlying positive selection in local, Alpine cattle breeds

We used genome-wide SNP data from 18 local cattle breeds from six countries of the Alpine region to characterize population structure and identify genomic regions underlying positive selection. The geographically close breeds Evolèner, Eringer, Valdostana Pezzata Nera, and Valdostana Castana were found to differ from all other Alpine breeds. In addition, three breeds, Simmental, and Original Braunvieh from Switzerland and Pinzgauer from Austria built three separate clusters. Of the 18 breeds studied, the intra-alpine Swiss breed Evolèner had the highest average inbreeding based on runs of homozygosity (F_ROH ) and the highest average genomic relationship within the breed. In contrast, Slovenian Cika cattle had the lowest average genomic inbreeding and the lowest average genomic relationship within the breed. We found selection signatures on chromosome 6 near known genes such as KIT and LCORL explaining variation in coat color and body size in cattle. The most prominent selection signatures were similar regardless of marker density and the breeds in the data set. In addition, using available high-density SNP data from 14 of the breeds we identified 47 genome regions as ROH islands. The proportion of homozygous animals was higher in all studied animals of local breeds than in Holstein and Brown Swiss cattle, the two most important commercial breeds in the Alpine region. We report ROH islands near genes related to thermoregulation, coat color, production, and stature. The results of this study serve as a basis for the search for causal variants underlying adaptation to the alpine environment and other specific characteristics selected during the evolution of local Alpine cattle breeds.

The role of trust and risk perception in current German nuclear waste management

One of the lessons learned in various countries that have to deal with spent nuclear fuel is that finding a proper place and siting a repository for high-level nuclear waste (HLW) cannot be achieved without public consent. After decades of obstruction, Germany recently launched a new, participatory, site-selection process for the disposal of HLW in deep geological formations. Nonetheless, significant opposition is assumed. Therefore, citizens' trust in the procedure and the agents involved may be paramount. We conducted an online survey (N ≈ 5000) in March/April 2020 to test a theoretical model on trust, perceived risks and benefits, and acceptance. We differentiated acceptance as a dependent variable according to distinct phases: the procedure, a possible decision on a disposal location, and the repository facility itself. The results show that trust is mainly important for explaining acceptance of the ongoing procedure and less so for the acceptance of the decision or the repository facility itself. Moreover, our investigation of the sample using a cluster analysis reveals characteristic patterns of trust, risk perception, and acceptance by three clusters: a cluster focusing on risk perception, an ambivalent cluster, and an indifferent cluster. Trust is lowest in the risk-focused cluster and highest in the ambivalent cluster.

Differential Analysis of Gly211Val and Gly286Val Mutations Affecting Sarco(endo)plasmic Reticulum Ca2+-ATPase (SERCA1) in Congenital Pseudomyotonia Romagnola Cattle

Congenital pseudomyotonia in cattle (PMT) is a rare skeletal muscle disorder, clinically characterized by stiffness and by delayed muscle relaxation after exercise. Muscle relaxation impairment is due to defective content of the Sarco(endo)plasmic Reticulum Ca²⁺ ATPase isoform 1 (SERCA1) protein, caused by missense mutations in the ATP2A1 gene. PMT represents the only mammalian model of human Brody myopathy. In the Romagnola breed, two missense variants occurring in the same allele were described, leading to Gly211Val and Gly286Val (G211V/G286V) substitutions. In this study, we analyzed the consequences of G211V and G286V mutations. Results support that the reduced amount of SERCA1 is a consequence of the G211V mutation, the G286V mutation almost being benign and the ubiquitin–proteasome system (UPS) being involved. After blocking the proteasome using a proteasome inhibitor, we found that the G211V mutant accumulates in cells at levels comparable to those of WT SERCA1. Our conclusion is that G211/286V mutations presumably originate in a folding-defective SERCA1 protein, recognized and diverted to degradation by UPS, although still catalytically functional, and that the main role is played by G211V mutation. Rescue of mutated SERCA1 to the sarcoplasmic reticulum membrane can re-establish resting cytosolic Ca²⁺ concentration and prevent the appearance of pathological signs, paving the way for a possible therapeutic approach against Brody disease.

A Missense Variant in PLP2 in Holstein Cattle with X-Linked Congenital Mast Cell Tumor

Congenital tumors occur infrequently in cattle. The aim of this study was to detail the clinicopathological phenotype of a Holstein calf with a congenital mast cell tumor and to identify the genetic cause by a whole-genome sequencing (WGS) trio-approach. An 18-day-old male Holstein calf was clinically examed and revealed multifocal, alopecic, thick and wrinkled skin lesions over the entire body. At 6 months of age, the general condition of the calf was characterized by retarded growth, poor nutritional status, and ulceration of the skin lesions. Histopathological examination revealed a primary cutaneous, poorly differentiated embryonal mast cell tumor with metastases in the lymph nodes and liver. Genetic analysis revealed a private X-linked variant in the PLP2 gene (chrX:87216480C > T; c.50C > T), which was present only in the genomes of the case (hemizygous) and his mother (heterozygous). It was absent in the sire as well as in 5365 control genomes. The identified missense variant exchanges the encoded amino acid of PLP2 at position 17 (p.Thr17Ile), which is classified as deleterious and affects a protein that plays a role in tumor growth and metastasis. Therefore, we suggested that the detected PLPL2 variant could be a plausible cause for this congenital condition in the affected calf.

A large deletion encompassing exon 2 of the ectodysplasin A (EDA) gene in a British blue crossbred calf with hypohidrotic ectodermal dysplasia

Background

Hypohidrotic ectodermal dysplasia (HED) is a congenital syndrome of mammals affecting organs and tissues of ectodermal origin characterized by absence or hypoplasia of hair, teeth, and eccrine glands. The disorder has been reported in several species, including humans, mice, dogs and cattle, associated with variants in genes affecting the ectodysplasin pathway, including the X-linked ectodysplasin A (EDA) gene. Until now, nine pathogenic variants have been found in the bovine EDA gene. Here we report a novel variant in EDA in a crossbreed male Belgian Blue calf with HED, and provide an overview of the phenotypic and allelic heterogeneity of EDA-related forms of HED in cattle.

Case presentation

A 45-day-old male crossbreed British Blue calf was referred with congenital hypotrichosis, oligodontia and omphalitis. On histopathological examination of the nasal planum, nasolabial glands and ducts were not observed. The density of hair follicles was low, and they were small, with a predominance of telogen-phase hairs, and some serocellular crusts. The phenotype of the calf resembled that of HED. Whole-genome sequencing (WGS) was performed and revealed a 21,899 base-pair deletion encompassing the coding exon 2 of EDA, predicted to result in an altered transcript and aberrant protein.

Conclusions

The clinicopathological and genetic findings were consistent with a case of X-linked HED. A very similar EDA deletion has been previously reported in a family of Holstein cattle with HED. The newly identified hemizygous EDA loss-of-function variant is certainly pathogenic and therefore is the genetic cause for the observed phenotype. This case report provides an additional example of the potential of WGS-based precise diagnostics in livestock species such as cattle to increase the diagnostic yield in rare diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13028-022-00641-2.

Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY, DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes Y1A, Y1B, Y2A and Y2B with a marked geographic partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and 7 wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during the migrations into northern Europe, eastern and southern Asia and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19^th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographic range.

Associated regions for multiple birth in Brown Swiss and Original Braunvieh cattle on chromosomes 15 and 11

Twin and multiple births have negative effects on the performance and health of cows and calves. To decipher the genetic architecture of this trait in the two Swiss Brown Swiss cattle populations, we performed various association analyses based on de-regressed breeding values. Genome-wide association analyses were executed using ~600 K imputed SNPs for the maternal multiple birth trait in ~3500 Original Braunvieh and ~7800 Brown Swiss animals. Significantly associated QTL were observed on different chromosomes for both breeds. We have identified on chromosome 11 a QTL that explains ~6% of the total genetic variance of the maternal multiple birth trait in Original Braunvieh. For the Brown Swiss breed, we have discovered a QTL on chromosome 15 that accounts for ~4% of the total genetic variance. For Original Braunvieh, subsequent haplotype analysis revealed a 90-kb window on chromosome 11 at 88 Mb, where a likely regulatory region is located close to the ID2 gene. In Brown Swiss, a 130-kb window at 75 Mb on chromosome 15 was identified. Analysis of whole-genome sequence data using linkage-disequilibrium estimation revealed possible causal variants for the identified QTL. A presumably regulatory variant in the non-coding 5' region of the ID2 gene was strongly associated with the haplotype for Original Braunvieh. In Brown Swiss, an intron variant in PRDM11, one 3' UTR variant in SYT13 and three intergenic variants 5' upstream of SYT13 were identified as candidate variants for the trait multiple birth maternal. In this study, we report for the first time QTL for the trait of multiple births in Original Braunvieh and Brown Swiss cattle. Moreover, our findings are another step towards a better understanding of the complex genetic architecture of this polygenic trait.

Short communication: Prevalence of deleterious variants causing recessive disorders in Italian Chianina, Marchigiana and Romagnola cattle

In the last two decades, the molecular cause of six monogenic autosomal recessive disorders has been identified in native Italian beef cattle: two different ATP2A1 variants for the pseudomyotonia congenita, the first in Chianina and Romagnola (PMT1) and the second in Romagnola (PMT2); a KDM2B variant for the paunch calf syndrome (PCS) in Marchigiana and Romagnola; a NID1 variant for the congenital cataract (CC) in Romagnola; a LAMB1 variant for the hemifacial microsomia (HFM) in Romagnola; an ABCA12 variant for the ichthyosis fetalis (IF) in Chianina and a FA2H variant for the ichthyosis congenita (IC) in Chianina. The aim of this study was to evaluate the potential impact of these disorders in the affected Italian populations. For this purpose, 3331 Chianina, 2812 Marchigiana and 1680 Romagnola bulls born in the last 40 years were considered. The allelic frequency (AF) of the variant for PMT1 was 1.0% in Romagnola, 4.6% in Marchigiana and 5.9% in Chianina. The AF of the variant for PMT2 was 3.3% in Romagnola and 0% in the other two breeds. The AF of the variant for PCS was 11.7% in Romagnola, 2.0% in Marchigiana and 0% in Chianina. The AF of the variants for CC, HFM, IF and IC resulted below 3%, being the variants detected only in the breed populations in which they were previously reported. Considering a selected male population in the single breed, Chianina showed carrier prevalence of 11.9% for PMT1, 7.7% for IC and 6.4% for IF. Romagnola showed carrier prevalence of 23.4% for PCS, 6.7% for PMT2, 4.1% for HFM, 3.2% for CC and 2.0% for PMT1. Marchigiana showed carrier prevalence of 9.1% for PMT1 and 4.0% for PCS. With respect to the Romagnola cattle, the concerning presence of a total of five defect alleles in the population hampers a general approach based on the prevention of carriers from artificial insemination. However, identification of carriers may allow conscious mating to prevent the risk of homozygous descendants as well as the spread of heterozygous offspring. Therefore, systematic genotyping for all seven known harmful alleles is recommended to prevent risk mating between carriers, in particular to avoid the occurrence of affected offspring.

The Complex and Diverse Genetic Architecture of the Absence of Horns (Polledness) in Domestic Ruminants, including Goats and Sheep

Horns are the most obvious common feature of Bovidae. The naturally occurring absence of horns in these species, also known as polledness, is of surprisingly heterogeneous nature, although they are Mendelian traits. This review compares in detail the molecular differences among the causes of inherited polledness in the domestic ruminant species of cattle, yak, sheep, and goat based on the causal gene variants that have been discovered in recent years. The genetic causes for the lack of horns in small ruminants seem not only to be more complex, e.g., in sheep, breed-specific characteristics are still unexplained, but in goats, there is also the associated disorder of intersexuality—polled intersex syndrome (PIS). In connection with animal welfare and the associated discussion about a legal ban on the dehorning of all farm animals, naturally hornless animals and the causal genetic variants are of increasing research interest in the age of genome editing. However, the low acceptance of genetic engineering in livestock, especially in European societies, limits its use in food-producing animals. Therefore, genotype-based targeted selection of naturally occurring variants is still a widely used method for spreading this desired trait within and across populations, at least in cattle and sheep.

Identification of two new recessive MC1R alleles in red-coloured Evolèner cattle and other breeds

Sequence variations in the melanocortin-1 receptor (MC1R) gene are associated with melanism in different animal species. Six functionally relevant alleles have been described in cattle to date. In a hypothesis-free approach we performed a genome-wide allelic association study with black, red and wild-coloured cattle of three Alpine cattle breeds (Eringer, Evolèner and Valdostana), revealing a single significant association signal close to the MC1R gene. We searched for candidate causative variants by sequencing the entire coding sequence and identified two novel protein-changing variants. We propose designating the mutant alleles at MC1R:c.424C>T as e^v1 and at MC1R:c.263G>A as e^v2 . Both affect conserved amino acid residues in functionally important transmembrane domains (p.Arg142Cys and p.Ser88Asn). Both alleles segregate predominantly in the Swiss Evolèner breed. They occur in other European cattle breeds such as Abondance and Rotes Höhenvieh as well. We observed almost perfect association between the MC1R genotypes and the coat colour phenotype in a cohort of 513 black, red and wild-coloured cattle. Animals carrying two copies of MC1R loss-of-function alleles or that were compound heterozygous for e, e^v1 , or e^v2 have a red to dark red (chestnut-like red) coat colour. These findings expand the spectrum of causal MC1R variants causing recessive red in cattle.

Four novel candidate causal variants for deficient homozygous haplotypes in Holstein cattle

Mendelian variants can determine both insemination success and neonatal survival and thus influence fertility and rearing success of cattle. We present 24 deficient homozygous haplotype regions in the Holstein population of Switzerland and provide an overview of the previously identified haplotypes in the global Holstein breed. This study encompasses massive genotyping, whole-genome sequencing (WGS) and phenotype association analyses. We performed haplotype screenings on almost 53 thousand genotyped animals including 114 k SNP data with two different approaches. We revealed significant haplotype associations to several survival, birth and fertility traits. Within haplotype regions, we mined WGS data of hundreds of bovine genomes for candidate causal variants, which were subsequently evaluated by using a custom genotyping array in several thousand breeding animals. With this approach, we confirmed the known deleterious SMC2:p.Phe1135Ser missense variant associated with Holstein haplotype (HH) 3. For two previously reported deficient homozygous haplotypes that show negative associations to female fertility traits, we propose candidate causative loss-of-function variants: the HH13-related KIR2DS1:p.Gln159* nonsense variant and the HH21-related NOTCH3:p.Cys44del deletion. In addition, we propose the RIOX1:p.Ala133_Glu142del deletion as well as the PCDH15:p.Leu867Val missense variant to explain the unexpected low number of homozygous haplotype carriers for HH25 and HH35, respectively. In conclusion, we demonstrate that with mining massive SNP data in combination with WGS data, we can map several haplotype regions and unravel novel recessive protein-changing variants segregating at frequencies of 1 to 5%. Our findings both confirm previously identified loci and expand the spectrum of undesired alleles impairing reproduction success in Holstein cattle, the world's most important dairy breed.

A de novo variant in the bovine ADAMTSL4 gene in an Original Braunvieh calf with congenital cataract

Inherited forms of cataract are a heterogeneous group of eye disorders known in livestock species. Clinicopathological analysis of a single case of impaired vision in a newborn Original Braunvieh calf revealed nuclear cataract. Whole‐genome sequencing of the parent‐offspring trio revealed a de novo mutation of ADAMTSL4 in this case. The heterozygous p.Arg776His missense variant affects a conserved residue of the ADAMTSL4 gene that encodes a secreted glycoprotein expressed in the lens throughout embryonic development. In humans, ADAMTSL4 genetic variants cause recessively inherited forms of subluxation of the lens. Given that ADAMTSL4 is a functional candidate gene for inherited disorders of the lens, we suggest that heterozygosity for the identified missense variant may have caused the congenital cataract in the affected calf. Cattle populations should be monitored for unexplained cataract cases, with subsequent DNA sequencing a hypothesized pathogenic effect of heterozygous ADAMTSL4 variants could be confirmed.

Runs of homozygosity in Swiss goats reveal genetic changes associated with domestication and modern selection

Background

The domestication of goat (Capra hircus) started 11,000 years ago in the fertile crescent. Breed formation in the nineteenth century, establishment of herd books, and selection for specific traits resulted in 10 modern goat breeds in Switzerland. We analyzed whole-genome sequencing (WGS) data from 217 modern goats and nine wild Bezoar goats (Capra aegagrus). After quality control, 27,728,288 biallelic single nucleotide variants (SNVs) were used for the identification of runs of homozygosity (ROH) and the detection of ROH islands.

Results

Across the 226 caprine genomes from 11 populations, we detected 344 ROH islands that harbor 1220 annotated genes. We compared the ROH islands between the modern breeds and the Bezoar goats. As a proof of principle, we confirmed a signature of selection, which contains the ASIP gene that controls several breed-specific coat color patterns. In two other ROH islands, we identified two missense variants, STC1:p.Lys139Arg and TSHR:p.Ala239Thr, which might represent causative functional variants for domestication signatures.

Conclusions

We have shown that the information from ROH islands using WGS data is suitable for the analysis of signatures of selection and allowed the detection of protein coding variants that may have conferred beneficial phenotypes during goat domestication. We hypothesize that the TSHR:p.Ala239Thr variant may have played a role in changing the seasonality of reproduction in modern domesticated goats. The exact functional significance of the STC1:p.Lys139Arg variant remains unclear and requires further investigation. Nonetheless, STC1 might represent a new domestication gene affecting relevant traits such as body size and/or milk yield in goats.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00695-w.

DSP missense variant in a Scottish Highland calf with congenital ichthyosis, alopecia, acantholysis of the tongue and corneal defects

BACKGROUND

Ichthyosis describes a localized or generalized hereditary cornification disorder caused by an impaired terminal keratinocyte differentiation resulting in excessive stratum corneum with the formation of more or less adherent scales. Ichthyosis affects humans and animals. Two rare bovine forms are reported, the severe harlequin ichthyosis and the less severe congenital ichthyosis, both characterized by a severe orthokeratotic lamellar hyperkeratosis.

RESULTS

A 2-weeks-old purebred Scottish Highland calf was referred because of a syndrome resembling congenital ichthyosis. The clinical phenotype included diffuse alopecia and a markedly lichenified skin covered with large and excessive scales. Additionally, conjunctivitis and ulceration of the cornea were noted. Post-mortem examination revealed deep fissures in the diffusely thickened tongue and histopathological findings in the skin confirmed the clinical diagnosis. Whole-genome sequencing of the affected calf and comparison of the data with control genomes was performed. A search for private variants in known candidate genes for skin phenotypes including genes related with erosive and hyperkeratotic lesions revealed a single homozygous protein-changing variant, DSP: c.6893 C>A, or p.Ala2298Asp. The variant is predicted to change a highly conserved residue in the C-terminal plakin domain of the desmoplakin protein, which represents a main intracellular component of desmosomes, important intercellular adhesion molecules in various tissues including epidermis. Sanger sequencing confirmed the variant was homozygous in the affected calf and heterozygous in both parents. Further genotyping of 257 Scottish Highland animals from Switzerland revealed an estimated allele frequency of 1.2%. The mutant allele was absent in more than 4800 controls from various other cattle breeds.

CONCLUSIONS

This study represents the first report of combined lesions compatible with congenital ichthyosis, alopecia, acantholysis of the tongue and corneal defects associated with a DSP missense variant as the most likely underlying cause. To the best of our knowledge, this study is also the first report of a DSP-related syndromic form of congenital ichthyosis in domestic animals. The results of our study enable genetic testing to avoid the unintentional occurrence of further affected cattle. The findings were added to the Online Mendelian Inheritance in Animals (OMIA) database (OMIA 002243-9913).

A Nonsense Variant in CCDC65 Gene Causes Respiratory Failure Associated with Increased Lamb Mortality in French Lacaune Dairy Sheep

We recently demonstrated that the Lacaune deficient homozygous haplotype 6 (LDHH6) potentially hosts a recessive perinatal lethal mutation in Lacaune dairy sheep mapped on OAR3. In the present study, we have analyzed the whole-genome sequences of two Lacaune ram heterozygous carriers of LDHH6. After variant calling and filtering against the variants of 86 non-carrier rams, we have identified a single nucleotide variant (SNV) in the two LDHH6 carriers whose variant allele induced a premature stop codon (p.Glu111*) in the Coiled-Coil Domain Containing 65 (CCDC65) gene. CCDC65 is involved in the assembly of the nexin-dynein regulatory complex for the formation of microtubules in ciliated cells. In order to identify the phenotype in homozygous sheep, we generated at-risk matings (n = 17) between rams and ewes heterozygous for the candidate variant in CCDC65. A total of 16 lambs were born alive with five genotyped as homozygous carriers. The homozygous lambs suffered from respiratory problems, and four of them died within the first month of life. At necropsy, we observed a broad hepatization of lung lobes possibly induced by infectious pneumonia. The management of this lethal recessive allele (frequency of 0.06) through reasoned mating in the Lacaune sheep selection schemes could reduce lamb mortality by 2%.

Mining massive genomic data of two Swiss Braunvieh cattle populations reveals six novel candidate variants that impair reproductive success

Background

This study was carried out on the two Braunvieh populations reared in Switzerland, the dairy Brown Swiss (BS) and the dual-purpose Original Braunvieh (OB). We performed a genome-wide analysis of array data of trios (sire, dam, and offspring) from the routine genomic selection to identify candidate regions showing missing homozygosity and phenotypic associations with five fertility, ten birth, and nine growth-related traits. In addition, genome-wide single SNP regression studies based on 114,890 single nucleotide polymorphisms (SNPs) for each of the two populations were performed. Furthermore, whole-genome sequencing data of 430 cattle including 70 putative haplotype carriers were mined to identify potential candidate variants that were validated by genotyping the current population using a custom array.

Results

Using a trio-based approach, we identified 38 haplotype regions for BS and five for OB that segregated at low to moderate frequencies. For the BS population, we confirmed two known haplotypes, BH1 and BH2. Twenty-four variants that potentially explained the missing homozygosity and associated traits were detected, in addition to the previously reported TUBD1:p.His210Arg variant associated with BH2. For example, for BS we identified a stop-gain variant (p.Arg57*) in the MRPL55 gene in the haplotype region on chromosome 7. This region is associated with the ‘interval between first and last insemination’ trait in our data, and the MRPL55 gene is known to be associated with early pregnancy loss in mice. In addition, we discuss candidate missense variants in the CPT1C, MARS2, and ACSL5 genes for haplotypes mapped in BS. In OB, we highlight a haplotype region on chromosome 19, which is potentially caused by a frameshift variant (p.Lys828fs) in the LIG3 gene, which is reported to be associated with early embryonic lethality in mice. Furthermore, we propose another potential causal missense variant in the TUBGCP5 gene for a haplotype mapped in OB.

Conclusions

We describe, for the first time, several haplotype regions that segregate at low to moderate frequencies and provide evidence of causality by trait associations in the two populations of Swiss Braunvieh. We propose a list of six protein-changing variants as potentially causing missing homozygosity. These variants need to be functionally validated and incorporated in the breeding program.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00686-3.

Reverse Genetic Screen for Deleterious Recessive Variants in the Local Simmental Cattle Population of Switzerland

We herein report the result of a large-scale reverse genetic screen in the Swiss Simmental population, a local dual-purpose cattle breed. We aimed to detect possible recessively inherited variants affecting protein-coding genes, as such deleterious variants can impair fertility and rearing success significantly. We used 115,000 phased SNP data of almost 10 thousand cattle with pedigree data. This revealed evidence for 11 genomic regions of 1.17 Mb on average, with haplotypes (SH1 to SH11) showing a significant depletion in homozygosity and an allele frequency between 3.2 and 10.6%. For the proposed haplotypes, it was unfortunately not possible to evaluate associations with fertility traits as no corresponding data were available. For each haplotype region, possible candidate genes were listed based on their known function in development and disease. Subsequent mining of single-nucleotide variants and short indels in the genomes of 23 sequenced haplotype carriers allowed us to identify three perfectly linked candidate causative protein-changing variants: a SH5-related DIS3:p.Ile678fs loss-of-function variant, a SH8-related CYP2B6:p.Ile313Asn missense variant, and a SH9-related NUBPL:p.Ser143Tyr missense variant. None of these variants occurred in homozygous state in any of more than 5200 sequenced cattle of various breeds. Selection against these alleles in order to reduce reproductive failure and animal loss is recommended.

A homozygous missense variant in laminin subunit beta 1 as candidate causal mutation of hemifacial microsomia in Romagnola cattle

Hemifacial microsomia (HFM) was diagnosed in a 9‐day‐old Romagnola calf. The condition was characterized by microtia of the left ear, anotia of the right ear, asymmetry of the face, and deafness. Magnetic resonance imaging revealed agenesis of the right pinna and both tympanic bullae, asymmetry of the temporal bones and temporomandibular joints, and right pontine meningocele. Brainstem auditory evoked responses confirmed the impaired auditory capacity. At gross post mortem examination, there was agenesis and hypoplasia of the right and the left external ear, respectively. No histological abnormalities were detected in the inner ears. A trio whole‐genome sequencing approach was carried out and identified a private homozygous missense variant in LAMB1 affecting a conserved residue (p.Arg668Cys). Genotyping of 221 Romagnola bulls revealed a carrier prevalence <2%. This represents a report of a LAMB1‐related autosomal recessive inherited disorder in domestic animals and adds LAMB1 to the candidate genes for HFM.

A frameshift insertion in FA2H causes a recessively inherited form of ichthyosis congenita in Chianina cattle

The aim of this study was to characterize the phenotype and to identify the genetic etiology of a syndromic form of ichthyosis congenita (IC) observed in Italian Chianina cattle and to estimate the prevalence of the deleterious allele in the population. Sporadic occurrence of different forms of ichthyosis including IC have been previously reported in cattle. However, so far, no causative genetic variant has been found for bovine IC. Nine affected cattle presenting congenital xerosis, hyperkeratosis and scaling of the skin as well as urolithiasis and cystitis associated with retarded growth were examined. Skin histopathology revealed a severe, diffuse orthokeratotic hyperkeratosis with mild to moderate epidermal hyperplasia. The pedigree records indicated a monogenic recessive trait. Homozygosity mapping and whole-genome sequencing allowed the identification of a homozygous frameshift 1 bp insertion in the FA2H gene (c.9dupC; p.Ala4ArgfsTer142) located in a 1.92 Mb shared identical-by-descent region on chromosome 18 present in all cases, while the parents were heterozygous as expected for obligate carriers. These findings enable the selection against this sub-lethal allele showing an estimated frequency of ~ 7.5% in Chianina top sires. A sporadic incidence of mild clinical signs in the skin of heterozygous carriers was observed. So far, pathogenic variants affecting the encoded fatty acid 2-hydroxylase catalyzing the synthesis of 2-hydroxysphingolipids have been associated with myelin disorders. In conclusion, this study represents the first report of an FA2H-related autosomal recessive inherited skin disorder in a mammalian species and adds FA2H to the list of candidate genes for ichthyosis in humans and animals. Furthermore, this study provides a DNA-based diagnostic test that enables selection against the identified pathogenic variant in the Chianina cattle population. However, functional studies are needed to better understand the expression of FA2H in IC-affected Chianina cattle.

Introgression of ASIP and TYRP1 Alleles Explains Coat Color Variation in Valais Goats

The Valais Blackneck goat is a Swiss goat breed with a characteristic coat color phenotype. Before the revision of the breed standard in 1938, 4 different color varieties of Valais goats were known. Besides Blackneck animals resembling the modern breed standard, the brown and white Copperneck goat, the white Capra Sempione, and the greyish Grüenochte comprised the historic Valais goats. The brown pigmentation of Copperneck goats had previously been traced back to an introgression of a mutant TYRP1 allele from Toggenburg goats. In the present study, we identified additional introgression events of distinct ASIP alleles causing the remaining 2 rare coat color patterns within the Valais Blackneck goat breed. We identified the introgression of the AWt allele from Appenzell or Saanen goats in white Capra Sempione goats. Similarly, introgression of the Apc allele from Peacock goats resulted in the greyish Grüenochte phenotype. These results demonstrate past hybridization events between breeds that are separated today. A perfect genotype-phenotype association in 393 Valais goats supported the causality of the genotyped variants for the different coat color phenotypes. Our study gives insights into the introgression of functionally relevant copy number variant (CNV) alleles controlling pigmentation between goat breeds with strikingly different coat color patterns.

Analysis of Polycerate Mutants Reveals the Evolutionary Co-option of HOXD1 for Horn Patterning in Bovidae

In the course of evolution, pecorans (i.e., higher ruminants) developed a remarkable diversity of osseous cranial appendages, collectively referred to as “headgear,” which likely share the same origin and genetic basis. However, the nature and function of the genetic determinants underlying their number and position remain elusive. Jacob and other rare populations of sheep and goats are characterized by polyceraty, the presence of more than two horns. Here, we characterize distinct POLYCERATE alleles in each species, both associated with defective HOXD1 function. We show that haploinsufficiency at this locus results in the splitting of horn bud primordia, likely following the abnormal extension of an initial morphogenetic field. These results highlight the key role played by this gene in headgear patterning and illustrate the evolutionary co-option of a gene involved in the early development of bilateria to properly fix the position and number of these distinctive organs of Bovidae.

GWAS Hits for Bilateral Convergent Strabismus with Exophthalmos in Holstein Cattle Using Imputed Sequence Level Genotypes

Bilateral convergent strabismus with exophthalmos (BCSE) is a malformation of the eyes and is recognized as a mild but progressive disorder that affects cattle in the first two years of life. This most likely inherited disorder is rarely described in cattle resembling autosomal dominantly inherited forms of human progressive external ophthalmoplegia (PEO). In German Braunvieh cattle, two linked genome regions were found that could be responsible for the development and/or progression of BCSE. The goal of this study was to phenotypically characterize BCSE in Holstein cattle from Germany and Switzerland as well as to identify associated genome regions by GWAS. The clinicopathological phenotype of 52 BCSE-affected Holstein cattle was in accordance with the phenotype described in German Braunvieh cattle, but in addition, signs of degeneration and cellular infiltration in the eye muscles were found. By using imputed sequence level genotype data, three genome-wide significant GWAS hits were revealed on different chromosomes that were not detected by initial GWAS based on high density SNP array data highlighting the usefulness of this approach for mapping studies. The associated genome regions include the ABCC4 gene as well as markers adjacent to the NCOR2 and DNAJC3 genes all illustrating possible functional candidate genes. Our results challenge a monogenic mode of inheritance and indicate a more complex inheritance of BCSE in Holstein cattle. Furthermore, in comparison to previous results from German Braunvieh cattle, it illustrates an obvious genetic heterogeneity causing BSCE in cattle. Subsequent whole genome sequencing (WGS)-based analyses might elucidate pathogenic variants in the future.

A major QTL at the LHCGR/FSHR locus for multiple birth in Holstein cattle

BACKGROUND

Twin and multiple births are rare in cattle and have a negative impact on the performance and health of cows and calves. Therefore, selection against multiple birth would be desirable in dairy cattle breeds such as Holstein. We applied different methods to decipher the genetic architecture of this trait using de-regressed breeding values for maternal multiple birth of ~ 2500 Holstein individuals to perform genome-wide association analyses using ~ 600 K imputed single nucleotide polymorphisms (SNPs).

RESULTS

In the population studied, we found no significant genetic trend over time of the estimated breeding values for multiple birth, which indicates that this trait has not been selected for in the past. In addition to several suggestive non-significant quantitative trait loci (QTL) on different chromosomes, we identified a major QTL on chromosome 11 for maternal multiple birth that explains ~ 16% of the total genetic variance. Using a haplotype-based approach, this QTL was fine-mapped to a 70-kb window on chromosome 11 between 31.00 and 31.07 Mb that harbors two functional candidate genes (LHCGR and FSHR). Analysis of whole-genome sequence data by linkage-disequilibrium estimation revealed a regulatory variant in the 5'-region of LHCGR as a possible candidate causal variant for the identified major QTL. Furthermore, the identified haplotype showed significant effects on stillbirth and days to first service.

CONCLUSIONS

QTL detection and subsequent identification of causal variants in livestock species remain challenging in spite of the availability of large-scale genotype and phenotype data. Here, we report for the first time a major QTL for multiple birth in Holstein cattle and provide evidence for a linked variant in the non-coding region of a functional candidate gene. This discovery, which is a first step towards the understanding of the genetic architecture of this polygenic trait, opens the path for future selection against this undesirable trait, and thus contributes to increased animal health and welfare.

A pathogenic HEXA missense variant in wild boars with Tay-Sachs disease

Gangliosidoses are inherited lysosomal storage disorders caused by reduced or absent activity of either a lysosomal enzyme involved in ganglioside catabolism, or an activator protein required for the proper activity of a ganglioside hydrolase, which results in the intra-lysosomal accumulation of undegraded metabolites. We hereby describe morphological, ultrastructural, biochemical and genetic features of GM2 gangliosidosis in three captive bred wild boar littermates. The piglets were kept in a partially-free range farm and presented progressive neurological signs, starting at 6 months of age. Animals were euthanized at approximately one year of age due to their poor conditions. Neuropathogens were excluded as a possible cause of the signs. Gross examination showed a reduction of cerebral and cerebellar consistency. Central (CNS) and peripheral (PNS) nervous system neurons were enlarged and foamy, with severe and diffuse cytoplasmic vacuolization. Transmission electron microscopy (TEM) of CNS neurons demonstrated numerous lysosomes, filled by parallel or concentric layers of membranous electron-dense material, defined as membranous cytoplasmic bodies (MCB). Biochemical composition of gangliosides analysis from CNS revealed accumulation of GM2 ganglioside; furthermore, Hex A enzyme activity was less than 1% compared to control animals. These data confirmed the diagnosis of GM2 gangliosidosis. Genetic analysis identified, at a homozygous level, the presence of a missense nucleotide variant c.1495C > T (p Arg499Cys) in the hexosaminidase subunit alpha gene (HEXA), located within the GH20 hexosaminidase superfamily domain of the encoded protein. This specific HEXA variant is known to be pathogenic and associated with Tay-Sachs disease in humans, but has never been identified in other animal species. This is the first report of a HEXA gene associated Tay-Sachs disease in wild boars and provides a comprehensive description of a novel spontaneous animal model for this lysosomal storage disease.

A Heterozygous Missense Variant in MAP2K2 in a Stillborn Romagnola Calf with Skeletal-Cardio-Enteric Dysplasia

Simple Summary

Skeletal dysplasias encompass a clinical-, pathological- and genetically heterogeneous group of disorders characterized by abnormal cartilage and/or bone formation, growth, and remodeling. They may belong to the so-called RASopathies, congenital conditions caused by heterozygous variants in genes that encode components of the Ras/mitogen-activated protein kinase (MAPK) cell signaling pathway. Herein, an affected calf of the Italian Romagnola breed was reported showing a skeletal-cardio-enteric dysplasia. We identified a most likely disease-causing mutation in the MAP2K2 gene by whole-genome sequencing (WGS). The MAP2K2 gene is known to be related with dominant inherited cardio-facio-cutaneous syndrome in humans, but it was so far unknown to cause a similar disease in domestic animals. We assume that the identified missense variant that was predicted to impair the function of the protein, occurred either within the germline of the dam or post-zygotically in the embryo. Rare lethal diseases such as the skeletal-cardio-enteric dysplasia in livestock are usually not characterized to the molecular level, mainly because of the lack of funds and diagnostic opportunities. Precise WGS-based diagnostics enables the understanding of rare diseases and supports the value of monitoring cattle breeding populations for fatal genetic defects.

Abstract

RASopathies are a group of developmental disorders caused by dominant mutations in genes that encode components of the Ras/mitogen-activated protein kinase (MAPK) cell signaling pathway. The goal of this study was to characterize the pathological phenotype of a Romagnola stillborn calf with skeletal-cardio-enteric dysplasia and to identify a genetic cause by whole-genome sequencing (WGS). The calf showed reduced fetal growth, a short-spine, a long and narrow face, cardiac defects and heterotopy of the spiral colon. Genetic analysis revealed a private heterozygous missense variant in MAP2K2:p.Arg179Trp, located in the protein kinase domain in the calf, and not found in more than 4500 control genomes including its sire. The identified variant affecting a conserved residue was predicted to be deleterious and most likely occurred de novo. This represents the first example of a dominant acting, and most likely pathogenic, variant in MAP2K2 in domestic animals, thereby providing the first MAP2K2-related large animal model, especially in respect to the enteric malformation. In addition, this study demonstrates the utility of WGS-based precise diagnostics for understanding sporadic congenital syndromic anomalies in cattle and the general utility of continuous surveillance for rare hereditary defects in cattle.

Genome-wide association study between copy number variants and hoof health traits in Holstein dairy cattle

Genome-wide association studies based on SNP have been completed for multiple traits in dairy cattle; however, copy number variants (CNV) could add genomic information that has yet to be harnessed. The objectives of this study were to identify CNV in genotyped Holstein animals and assess their association with hoof health traits using deregressed estimated breeding values as pseudophenotypes. A total of 23,256 CNV comprising 1,645 genomic regions were identified in 5,845 animals. Fourteen genomic regions harboring structural variations, including 9 deletions and 5 duplications, were associated with at least 1 of the studied hoof health traits. This group of traits included digital dermatitis, interdigital dermatitis, heel horn erosion, sole ulcer, white line lesion, sole hemorrhage, and interdigital hyperplasia; no regions were associated with toe ulcer. Twenty candidate genes overlapped with the regions associated with these traits including SCART1, NRXN2, KIF26A, GPHN, and OR7A17. In this study, an effect on infectious hoof lesions could be attributed to the PRAME (Preferentially Expressed Antigen in Melanoma) gene. Almost all genes detected in association with noninfectious hoof lesions could be linked to known metabolic disorders. The knowledge obtained considering information of associated CNV to the traits of interest in this study could improve the accuracy of estimated breeding values. This may further increase the genetic gain for these traits in the Canadian Holstein population, thus reducing the involuntary animal losses due to lameness.

Bovine viral diarrhoea virus 1b infection associated with congenital tremor and hypomyelination in Holstein calves

Hypomyelination is a rare consequence of in utero bovine viral diarrhoea virus (BVDV) infection. We describe a BVDV outbreak in a naïve Holstein dairy herd in northern Italy, with an unusually high prevalence of calves with neurological signs, generalised tremors and ataxia. Histological analysis showed that hypomyelination was the predominant lesion and that the most typical BVDV neuropathological findings (e.g. cerebellar hypoplasia) were absent. Virological and molecular analyses showed that non-cytopathic BVDV genotype 1b was associated with the calves' neurological signs and excluded other viruses responsible for congenital infection or neurological disorders. Whole-genome sequencing of BVDVs from the brain of a calf with neurological signs and the whole blood of a persistently infected herd-mate with no such sign showed >99.7 % sequence identity. Analysis of the quasispecies distribution revealed the greatest variation rates in regions coding for the structural proteins E1 and E2. Variation was slightly greater in the brain- than in the blood-derived sequence and occurred at different sites, suggesting the occurrence of distinct evolutionary processes in the two persistently infected calves. Molecular characterisation of BVDV genomes from five other calves with neurological signs from the same farm confirmed that the E1 and E2 regions were the most variable. Several factors, including genetic variability and host factors, appear to have contributed to the observed unique BVDV disease phenotype, characterised by hypomyelination and neurological signs.