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

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

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

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.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

A De Novo Mutation in COL1A1 in a Holstein Calf with Osteogenesis Imperfecta Type II

Osteogenesis imperfecta (OI) type II is a genetic connective tissue disorder characterized by bone fragility, severe skeletal deformities and shortened limbs. OI usually causes perinatal death of affected individuals. OI type II diagnosis in humans is established by the identification of heterozygous mutations in genes coding for collagens. The purpose of this study was to characterize the pathological phenotype of an OI type II-affected neonatal Holstein calf and to identify the causative genetic variant by whole-genome sequencing (WGS). The calf had acute as well as intrauterine fractures, abnormally shaped long bones and localized arthrogryposis. Genetic analysis revealed a private heterozygous missense variant in COL1A1 (c.3917T>A) located in the fibrillar collagen NC1 domain (p.Val1306Glu) that most likely occurred de novo. This confirmed the diagnosis of OI type II and represents the first report of a pathogenic variant in the fibrillar collagen NC domain of COL1A1 associated to OI type II in domestic animals. Furthermore, this study highlights the utility of WGS-based precise diagnostics for understanding congenital disorders in cattle and the need for continued surveillance for rare lethal genetic disorders in cattle

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

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

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

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

A de novo start-lost variant in ANKRD28 in a Holstein calf with dwarfism.

N