Search

In modern cattle production, we have seen a negative trend for decades in reproduction while productivity and performance have improved. Although considered genetically complex, part of these fecundity, fertility, and rearing success issues are caused by Mendelian monogenic disorders. Traditionally, such disorders are investigated opportunistically based on their sporadic occurrence and through subsequent targeted analysis of affected individuals. This approach is called the forward genetic approach (FGA). Modern genomic technologies, such as single nucleotide polymorphism (SNP) array genotyping and whole-genome sequencing (WGS), allow for straightforward locus mapping and the identification of candidate causal variants in affected individuals or families. Nevertheless, a major drawback is the arbitrary sampling and availability of well-phenotyped individuals for research, especially for mostly invisible defects affecting fecundity, early embryonic death, and abortions. Therefore, the reverse genetic approach (RGA) is applied to screen for underlying recessive lethal or sub-lethal variants. This approach requires the availability of massive population-wide genomic data. By applying a haplotype screen for a significant deviation of the Hardy-Weinberg equilibrium, genomic regions potentially harboring candidate causal variants are identified. The subsequent generation of WGS data of haplotype carriers allows for the mining for pathogenic variants potentially causing a reduction in homozygosity. In the first part of my thesis, I present 18 successful, 1 inconclusive example, and 1 example addressing co-dominant effects of a known disorder. These FGA analyzes include heritable skin (n=7), bone (n=7), neuromuscular (n=1), eye (n=2), as well as syndromic disorders (n=3) in various European cattle breeds. Missense and frameshift variants in the IL17RA, DSP, and FA2H genes were described in three recessive genodermatoses: immunodeficiency with psoriasis-like skin alterations, syndromic ichthyosis, and ichthyosis congenita, respectively. Hypohidrotic ectodermal dysplasia was described as X-linked disorder that is associated with a gross deletion in the EDA gene. In dominant genodermatoses, a missense variant in COL5A2 was shown to lead to classical Ehlers-Danlos syndrome, an in-frame deletion in KRT5 was shown to cause epidermolysis bullosa simplex, and results of a study using an individual case of juvenile angiomatosis remained inconclusive. A recessive disorder described as hemifacial macrosomia was associated with a missense variant in LAMB1. Chondrodysplasia in a single family was shown to be caused by a de novo mutation in the bull leading to a stop-loss of the gene FGFR3. De novo mutations (missense and large deletions) in the COL2A1 and COL1A1 genes were associated with achondrogenesis type II (bulldog calf syndrome), and osteogenesis imperfecta type II, respectively. Another mutation that we found to affect bone morphology was a trisomy in chromosome 29 leading to proportional dwarfism with facial dysplasia. Congenital neuromuscular channelopathy was for the first time associated with a missense variant in KCNG1. Furthermore, a de novo missense variant in ADAMTSL4 and a recessive missense variant in CNGB3 were shown to cause congenital cataract and achromatopsia, respectively. Additionally, cases of pulmonary hypoplasia and anasarca syndrome were analyzed and shown to be caused by trisomy 20 in two unrelated calves and a recessively inherited missense variant in ADAMTS3. Moreover, the fatal syndromic disorder skeletal-cardo-enteric dysplasia was described to be caused by a de novo missense variant in MAP2K2. Finally, I investigated the effects on blood cholesterol and triglyceride levels of heterozygous carriers of the previously described APOB-related cholesterol deficiency. In the second part of my thesis, I present the outcome of the RGA in four main Swiss populations, that was validated with the SWISScow custom array. In the Brown Swiss dairy population, 72 haplotype regions showed significant depletions in homozygosity. Four of these haplotypes (BH6, BH14, BH24, and BH34) were associated with missense and nonsense variants in different genes (MARS2, MRPL55, CPT1C, and ACSL5, respectively). In the Original Braunvieh population, eight haplotype regions were identified. Candidate causal variants included a missense variant in TUBGCP5 gene associated with haplotype OH2, and a splice site frameshift variant in LIG3 gene associated with haplotype OH4. In the Holstein population, 24 haplotype regions were identified with a significant reduction of homozygosity. Subsequently, four novel candidate variants were proposed: a nonsense variant in KIR2DS1 for haplotype HH13, in-frame deletion in the genes NOTCH3 for HH21 haplotype, and RIOX1 for HH25 haplotype, and finally, a missense variant in PCDH15 for HH35 haplotype. In the Simmental population, eleven haplotype regions were detected. The haplotype SH5 was associated with a frameshift variant in DIS3 gene and the haplotypes SH8 and SH9 with missense variants in the CYP2B6 and NUBPL genes, respectively. For the breeds Brown Swiss, Original Braunvieh, and Holstein, association studies were carried out including traits describing fertility, birth, growth, and survival. Thereby most of the described mentioned haplotypes show additive effects. Regardless of the approach, all the described candidate causal variants can be used as a tool of precision diagnostics and represent a step forward towards personalized medicine in cattle. Furthermore, these variants can be easily genotyped and allow for targeted breeding to reduce the number of risk matings, which would lead to a reduction of affected animals and significant improvement in animal health and welfare.

Additional file 1. British blue cross calf with hypohidrotic ectodermal dysplasia. A. Note the areas of erosion and the variation in degree of alopecia on the thorax (A), over the carpi (B), and over the stifle joint (C).

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.

Entropion is a known congenital disorder in sheep presumed to be heritable but no causative genetic variant has been reported. Affected lambs show a variable inward rolling of the lower eyelids leading to blindness in severe cases. In Switzerland, the Swiss White Alpine (SWA) breed showed a significantly higher prevalence for entropion than other breeds. A GWAS using 150 SWA sheep (90 affected lambs and 60 controls), based on 600k SNP data, revealed a genome-wide significant signal on chromosome 15. The 0.2 Mb associated region contains functional candidate genes, SMTNL1 and CTNND1. Pathogenic variants in human CTNND1 cause blepharocheilodontic syndrome 2, a rare disorder including eyelid anomalies, and SMTNL1 regulates contraction and relaxation of skeletal and smooth muscle. WGS of a single entropion-affected lamb revealed two private missense variants in SMTNL1 and CTNND1. Subsequent genotyping of both variants in 231 phenotyped SWA sheep was performed. The SMTNL1 variant p.(Asp452Asn) affects an evolutionary conserved residue within an important domain and represents a rare allele, which occurred also in controls. The p.(Glu943Lys) variant in CTNND1 represents a common variant unlikely to cause entropion as the mutant allele occurred more frequently in non-affected sheep. Therefore, we propose that these protein-changing variants are unlikely to explain the phenotype. Additionally, WGS of three further disconcordant pairs of full siblings was carried out but revealed no obvious causative variant. Finally, we conclude that entropion represents a more complex disease caused by different non-coding regulatory variants.

In 2015, cholesterol deficiency (CD) was reported for the first time as a new recessive defect in Holstein cattle. After GWAS mapping and identification of a disease-associated haplotype, a causative loss-of-function variant in APOB was identified. CD-clinically affected APOB homozygotes showed poor development, intermittent diarrhea and hypocholesterolemia and, consequently, a limited life expectation. Herein, we present a collection of 18 cases clinically diagnosed as CD-affected APOB heterozygotes. CD-clinically affected heterozygotes show reduced cholesterol and triglyceride blood concentrations. The differences in total blood cholesterol and triglycerides between nine CD-clinically affected and 36 non-affected heterozygotes were significant. As only some APOB heterozygotes show the clinical CD phenotype, we assume that the penetrance is reduced in heterozygotes compared to the fully penetrant effect observed in homozygotes. We conclude that APOB-associated CD represents most likely an incomplete dominant inherited metabolic disease with incomplete penetrance in heterozygotes.

Two clinical forms of ichthyosis in cattle have been reported, ichthyosis fetalis and congenital ichthyosis. Ichthyosis poses animal welfare and economic issues and the more severe form, ichthyosis fetalis, is lethal. A Shorthorn calf with ichthyosis fetalis was investigated and a likely causal missense variant on chromosome 2 in the ABCA12 gene (NM_001191294.2:c.6776T>C) was identified by whole genome sequencing. Mutations in the ABCA12 gene are known to cause ichthyosis fetalis in cattle and Harlequin ichthyosis in humans. Sanger sequencing of the affected calf and the dam confirmed the variant was homozygous in the affected calf and heterozygous in the dam. Further genotyping of 130 Shorthorn animals from the same property revealed an estimated allele frequency of 3.8%. The presented findings enable genetic testing for breeding and diagnostics.

A specific white spotting phenotype, termed finching or line-backed spotting, is known for all Pinzgauer cattle and occurs occasionally in Tux-Zillertaler cattle, two Austrian breeds. The so-called Pinzgauer spotting is inherited as an autosomal incompletely dominant trait. A genome-wide association study using 27 white spotted and 16 solid-coloured Tux-Zillertaler cattle, based on 777k SNP data, revealed a strong signal on chromosome 6 at the KIT locus. Haplotype analyses defined a critical interval of 122 kb downstream of the KIT coding region. Whole-genome sequencing of a Pinzgauer cattle and comparison to 338 control genomes revealed a complex structural variant consisting of a 9.4-kb deletion and an inversely inserted duplication of 1.5 kb fused to a 310-kb duplicated segment from chromosome 4. A diagnostic PCR was developed for straightforward genotyping of carriers for this structural variant (KIT ) and confirmed that the variant allele was present in all Pinzgauer and most of the white spotted Tux-Zillertaler cattle. In addition, we detected the variant in all Slovenian Cika, British Gloucester and Spanish Berrenda en negro cattle with similar spotting patterns. Interestingly, the KIT variant occurs in some white spotted animals of the Swiss breeds Evolèner and Eringer. The introgression of the KIT variant confirms admixture and the reported historical relationship of these short-headed breeds with Austrian Tux-Zillertaler and suggests a mutation event, occurring before breed formation.