Search

Domestication and human selection have formed diverse goat breeds with characteristic phenotypes. This process correlated with the fixation of causative genetic variants controlling breed-specific traits within regions of reduced genetic diversity, so called selection signatures or selective sweeps. Using whole genome sequencing of DNA pools (pool-seq) from 20 genetically diverse modern goat breeds and bezoars, we identified 2,239 putative selection signatures. In two Pakistani goat breeds, Pak Angora and Barbari, we found selection signatures in a region harboring KIT, a gene involved in melanoblast development, migration, and survival. The search for candidate causative variants responsible for these selective sweeps revealed two different copy number variants (CNVs) downstream of KIT that were exclusively present in white Pak Angora and white-spotted Barbari goats. Several Swiss goat breeds selected for specific coat colors showed selection signatures at the ASIP locus encoding the agouti signaling protein. Analysis of these selective sweeps revealed four different CNVs associated with the white or tan (AWt), Swiss markings (Asm), badgerface (Ab), and the newly proposed peacock (Apc) allele. RNA-seq analyses on skin samples from goats with the different CNV alleles suggest that the identified structural variants lead to an altered expression of ASIP between eumelanistic and pheomelanistic body areas. Our study yields novel insights into the genetic control of pigmentation by identifying six functionally relevant CNVs. It illustrates how structural changes of the genome have contributed to phenotypic evolution in domestic goats., Competing Interests: The authors have declared that no competing interests exist.

Inherited bleeding disorders including abnormalities of platelet number and function rarely occur in a variety of dog breeds, but are probably underdiagnosed. Genetically characterized canine forms of platelet disorders provide valuable large animal models for understanding similar platelet disorders in people. Breed-specific disease associated genetic variants in only eight different genes are known to cause intrinsic platelet disorders in dogs. However, the causative genetic variant in many dog breeds has until now remained unknown. Four cases of a mild to severe bleeding disorder in Cocker Spaniel dogs are herein presented. The affected dogs showed a platelet adhesion defect characterized by macrothrombocytopenia with variable platelet counts resembling human Bernard-Soulier syndrome (BSS). Furthermore, the lack of functional GPIb-IX-V was demonstrated by immunocytochemistry. Whole genome sequencing of one affected dog and visual inspection of the candidate genes identified a deletion in the glycoprotein IX platelet (GP9) gene. The GP9 gene encodes a subunit of a platelet surface membrane glycoprotein complex; this functions as a receptor for von Willebrand factor, which initiates the maintenance of hemostasis after injury. Variants in human GP9 are associated with Bernard-Soulier syndrome, type C. The deletion spanned 2460 bp, and included a significant part of the single coding exon of the canine GP9 gene on dog chromosome 20. The variant results in a frameshift and premature stop codon which is predicted to truncate almost two-thirds of the encoded protein. PCR-based genotyping confirmed recessive inheritance. The homozygous variant genotype seen in affected dogs did not occur in 98 control Cocker Spaniels. Thus, it was concluded that the structural variant identified in the GP9 gene was most likely causative for the BSS-phenotype in the dogs examined. These findings provide the first large animal GP9 model for this group of inherited platelet disorders and greatly facilitate the diagnosis and identification of affected and/or normal carriers in Cocker Spaniels., Competing Interests: MET and MD are shareholders and affiliated at Genefast and Genomia, commercial laboratories carrying out genetic testing. This does not alter our adherence to PLOS ONE policies on sharing data and materials. There are no patents products in development or marketed products to declare. All authors are committed to adhere to PLOS ONE policies on sharing data and materials.

Belted cattle have a circular belt of unpigmented hair and skin around their midsection. The belt is inherited as a monogenic autosomal dominant trait. We mapped the causative variant to a 37 kb segment on bovine chromosome 3. Whole genome sequence data of 2 belted and 130 control cattle yielded only one private genetic variant in the critical interval in the two belted animals. The belt-associated variant was a copy number variant (CNV) involving the quadruplication of a 6 kb non-coding sequence located approximately 16 kb upstream of the TWIST2 gene. Increased copy numbers at this CNV were strongly associated with the belt phenotype in a cohort of 333 cases and 1322 controls. We hypothesized that the CNV causes aberrant expression of TWIST2 during neural crest development, which might negatively affect melanoblasts. Functional studies showed that ectopic expression of bovine TWIST2 in neural crest in transgenic zebrafish led to a decrease in melanocyte numbers. Our results thus implicate an unsuspected involvement of TWIST2 in regulating pigmentation and reveal a non-coding CNV underlying a captivating Mendelian character.

One to two percent of all children are born with a developmental disorder requiring pediatric hospital admissions. For many such syndromes, the molecular pathogenesis remains poorly characterized. Parallel developmental disorders in other species could provide complementary models for human rare diseases by uncovering new candidate genes, improving the understanding of the molecular mechanisms and opening possibilities for therapeutic trials. We performed various experiments, e.g. combined genome-wide association and next generation sequencing, to investigate the clinico-pathological features and genetic causes of three developmental syndromes in dogs, including craniomandibular osteopathy (CMO), a previously undescribed skeletal syndrome, and dental hypomineralization, for which we identified pathogenic variants in the canine SLC37A2 (truncating splicing enhancer variant), SCARF2 (truncating 2-bp deletion) and FAM20C (missense variant) genes, respectively. CMO is a clinical equivalent to an infantile cortical hyperostosis (Caffey disease), for which SLC37A2 is a new candidate gene. SLC37A2 is a poorly characterized member of a glucose-phosphate transporter family without previous disease associations. It is expressed in many tissues, including cells of the macrophage lineage, e.g. osteoclasts, and suggests a disease mechanism, in which an impaired glucose homeostasis in osteoclasts compromises their function in the developing bone, leading to hyperostosis. Mutations in SCARF2 and FAM20C have been associated with the human van den Ende-Gupta and Raine syndromes that include numerous features similar to the affected dogs. Given the growing interest in the molecular characterization and treatment of human rare diseases, our study presents three novel physiologically relevant models for further research and therapy approaches, while providing the molecular identity for the canine conditions.

Clinical, pathological and genetic examination revealed an as yet uncharacterized juvenile-onset neuroaxonal dystrophy (NAD) in Spanish water dogs. Affected dogs presented with various neurological deficits including gait abnormalities and behavioral deficits. Histopathology demonstrated spheroid formation accentuated in the grey matter of the cerebral hemispheres, the cerebellum, the brain stem and in the sensory pathways of the spinal cord. Iron accumulation was absent. Ultrastructurally spheroids contained predominantly closely packed vesicles with a double-layered membrane, which were characterized as autophagosomes using immunohistochemistry. The family history of the four affected dogs suggested an autosomal recessive inheritance. SNP genotyping showed a single genomic region of extended homozygosity of 4.5 Mb in the four cases on CFA 8. Linkage analysis revealed a maximal parametric LOD score of 2.5 at this region. By whole genome re-sequencing of one affected dog, a perfectly associated, single, non-synonymous coding variant in the canine tectonin beta-propeller repeat-containing protein 2 (TECPR2) gene affecting a highly conserved region was detected (c.4009C>T or p.R1337W). This canine NAD form displays etiologic parallels to an inherited TECPR2 associated type of human hereditary spastic paraparesis (HSP). In contrast to the canine NAD, the spinal cord lesions in most types of human HSP involve the sensory and the motor pathways. Furthermore, the canine NAD form reveals similarities to cases of human NAD defined by widespread spheroid formation without iron accumulation in the basal ganglia. Thus TECPR2 should also be considered as candidate gene for human NAD. Immunohistochemistry and the ultrastructural findings further support the assumption, that TECPR2 regulates autophagosome accumulation in the autophagic pathways. Consequently, this report provides the first genetic characterization of juvenile canine NAD, describes the histopathological features associated with the TECPR2 mutation and provides evidence to emphasize the association between failure of autophagy and neurodegeneration.

The identification of quantitative trait loci (QTL) such as height and their underlying causative variants is still challenging and often requires large sample sizes. In humans hundreds of loci with small effects control the heritable portion of height variability. In domestic animals, typically only a few loci with comparatively large effects explain a major fraction of the heritability. We investigated height at withers in Shetland ponies and mapped a QTL to ECA 6 by genome-wide association (GWAS) using a small cohort of only 48 animals and the Illumina equine SNP70 BeadChip. Fine-mapping revealed a shared haplotype block of 793 kb in small Shetland ponies. The HMGA2 gene, known to be associated with height in horses and many other species, was located in the associated haplotype. After closing a gap in the equine reference genome we identified a non-synonymous variant in the first exon of HMGA2 in small Shetland ponies. The variant was predicted to affect the functionally important first AT-hook DNA binding domain of the HMGA2 protein (c.83G>A; p.G28E). We assessed the functional impact and found impaired DNA binding of a peptide with the mutant sequence in an electrophoretic mobility shift assay. This suggests that the HMGA2 variant also affects DNA binding in vivo and thus leads to reduced growth and a smaller stature in Shetland ponies. The identified HMGA2 variant also segregates in several other pony breeds but was not found in regular-sized horse breeds. We therefore conclude that we identified a quantitative trait nucleotide for height in horses.

Four related cows showed hairless streaks on various parts of the body with no correlation to the pigmentation pattern. The stripes occurred in a consistent pattern resembling the lines of Blaschko. The non-syndromic hairlessness phenotype observed occurred across three generations of a single family and was compatible with an X-linked mode of inheritance. Linkage analysis and subsequent whole genome sequencing of one affected female identified two perfectly associated non-synonymous sequence variants in the critical interval on bovine chromosome X. Both variants occurred in complete linkage disequilibrium and were absent in more than 3900 controls. An ERCC6L missense mutation was predicted to cause an amino acid substitution of a non-conserved residue. Analysis in mice showed no specific Ercc6l expression pattern related to hair follicle development and therefore ERCC6L was not considered as causative gene. A point mutation at the 5'-splice junction of exon 5 of the TSR2, 20S rRNA accumulation, homolog (S. cerevisiae), gene led to the production of two mutant transcripts, both of which contain a frameshift and generate a premature stop codon predicted to truncate approximately 25% of the protein. Interestingly, in addition to the presence of both physiological TSR2 transcripts, the two mutant transcripts were predominantly detected in the hairless skin of the affected cows. Immunohistochemistry, using an antibody against the N-terminal part of the bovine protein demonstrated the specific expression of the TSR2 protein in the skin and the hair of the affected and the control cows as well as in bovine fetal skin and hair. The RNA hybridization in situ showed that Tsr2 was expressed in pre- and post-natal phases of hair follicle development in mice. Mammalian TSR2 proteins are highly conserved and are known to be broadly expressed, but their precise in vivo functions are poorly understood. Thus, by dissecting a naturally occurring mutation in a domestic animal species, we identified TSR2 as a regulator of hair follicle development.

Whereas the genetic background of horn growth in cattle has been studied extensively, little is known about the morphological changes in the developing fetal horn bud. In this study we histologically analyzed the development of horn buds of bovine fetuses between ~70 and ~268 days of pregnancy and compared them with biopsies taken from the frontal skin of the same fetuses. In addition we compared the samples from the wild type (horned) fetuses with samples taken from the horn bud region of age-matched genetically hornless (polled) fetuses. In summary, the horn bud with multiple layers of vacuolated keratinocytes is histologically visible early in fetal life already at around day 70 of gestation and can be easily differentiated from the much thinner epidermis of the frontal skin. However, at the gestation day (gd) 212 the epidermis above the horn bud shows a similar morphology to the epidermis of the frontal skin and the outstanding layers of vacuolated keratinocytes have disappeared. Immature hair follicles are seen in the frontal skin at gd 115 whereas hair follicles below the horn bud are not present until gd 155. Interestingly, thick nerve bundles appear in the dermis below the horn bud at gd 115. These nerve fibers grow in size over time and are prominent shortly before birth. Prominent nerve bundles are not present in the frontal skin of wild type or in polled fetuses at any time, indicating that the horn bud is a very sensitive area. The samples from the horn bud region from polled fetuses are histologically equivalent to samples taken from the frontal skin in horned species. This is the first study that presents unique histological data on bovine prenatal horn bud differentiation at different developmental stages which creates knowledge for a better understanding of recent molecular findings.

Inherited neurodegenerative disorders are debilitating diseases that occur across different species. We have performed clinical, pathological and genetic studies to characterize a novel canine neurodegenerative disease present in the Lagotto Romagnolo dog breed. Affected dogs suffer from progressive cerebellar ataxia, sometimes accompanied by episodic nystagmus and behavioral changes. Histological examination revealed unique pathological changes, including profound neuronal cytoplasmic vacuolization in the nervous system, as well as spheroid formation and cytoplasmic aggregation of vacuoles in secretory epithelial tissues and mesenchymal cells. Genetic analyses uncovered a missense change, c.1288G>A; p.A430T, in the autophagy-related ATG4D gene on canine chromosome 20 with a highly significant disease association (p = 3.8 x 10-136) in a cohort of more than 2300 Lagotto Romagnolo dogs. ATG4D encodes a poorly characterized cysteine protease belonging to the macroautophagy pathway. Accordingly, our histological analyses indicated altered autophagic flux in affected tissues. The knockdown of the zebrafish homologue atg4da resulted in a widespread developmental disturbance and neurodegeneration in the central nervous system. Our study describes a previously unknown canine neurological disease with particular pathological features and implicates the ATG4D protein as an important autophagy mediator in neuronal homeostasis. The canine phenotype serves as a model to delineate the disease-causing pathological mechanism(s) and ATG4D function, and can also be used to explore treatment options. Furthermore, our results reveal a novel candidate gene for human neurodegeneration and enable the development of a genetic test for veterinary diagnostic and breeding purposes.

Dandy-Walker-like malformation (DWLM) is the result of aberrant brain development and mainly characterized by cerebellar hypoplasia. DWLM affected dogs display a non-progressive cerebellar ataxia. Several DWLM cases were recently observed in the Eurasier dog breed, which strongly suggested a monogenic autosomal recessive inheritance in this breed. We performed a genome-wide association study (GWAS) with 9 cases and 11 controls and found the best association of DWLM with markers on chromosome 1. Subsequent homozygosity mapping confirmed that all 9 cases were homozygous for a shared haplotype in this region, which delineated a critical interval of 3.35 Mb. We sequenced the genome of an affected Eurasier and compared it with the Boxer reference genome and 47 control genomes of dogs from other breeds. This analysis revealed 4 private non-synonymous variants in the critical interval of the affected Eurasier. We genotyped these variants in additional dogs and found perfect association for only one of these variants, a single base deletion in the VLDLR gene encoding the very low density lipoprotein receptor. This variant, VLDLR:c.1713delC is predicted to cause a frameshift and premature stop codon (p.W572Gfs*10). Variants in the VLDLR gene have been shown to cause congenital cerebellar ataxia and mental retardation in human patients and Vldlr knockout mice also display an ataxia phenotype. Our combined genetic data together with the functional knowledge on the VLDLR gene from other species thus strongly suggest that VLDLR:c.1713delC is indeed causing DWLM in Eurasier dogs.

Cataract is a known condition leading to opacification of the eye lens causing partial or total blindness. Mutations are known to cause autosomal dominant or recessive inherited forms of cataracts in humans, mice, rats, guinea pigs and dogs. The use of large-sized animal models instead of those using mice for the study of this condition has been discussed due to the small size of rodent lenses. Four juvenile-onset cases of bilateral incomplete immature nuclear cataract were recently observed in Romagnola cattle. Pedigree analysis suggested a monogenic autosomal recessive inheritance. In addition to the cataract, one of the cases displayed abnormal head movements. Genome-wide association and homozygosity mapping and subsequent whole genome sequencing of a single case identified two perfectly associated sequence variants in a critical interval of 7.2 Mb on cattle chromosome 28: a missense point mutation located in an uncharacterized locus and an 855 bp deletion across the exon 19/intron 19 border of the bovine nidogen 1 (NID1) gene (c.3579_3604+829del). RT-PCR showed that NID1 is expressed in bovine lenses while the transcript of the second locus was absent. The NID1 deletion leads to the skipping of exon 19 during transcription and is therefore predicted to cause a frameshift and premature stop codon (p.1164fs27X). The truncated protein lacks a C-terminal domain essential for binding with matrix assembly complexes. Nidogen 1 deficient mice show neurological abnormalities and highly irregular crystal lens alterations. This study adds NID1 to the list of candidate genes for inherited cataract in humans and is the first report of a naturally occurring mutation leading to non-syndromic catarct in cattle provides a potential large animal model for human cataract.

Congenital hepatic fibrosis has been described as a lethal disease with monogenic autosomal recessive inheritance in the Swiss Franches-Montagnes horse breed. We performed a genome-wide association study with 5 cases and 12 controls and detected an association on chromosome 20. Subsequent homozygosity mapping defined a critical interval of 952 kb harboring 10 annotated genes and loci including the polycystic kidney and hepatic disease 1 (autosomal recessive) gene (PKHD1). PKHD1 represents an excellent functional candidate as variants in this gene were identified in human patients with autosomal recessive polycystic kidney and hepatic disease (ARPKD) as well as several mouse and rat mutants. Whereas most pathogenic PKHD1 variants lead to polycystic defects in kidney and liver, a small subset of the human ARPKD patients have only liver symptoms, similar to our horses with congenital hepatic fibrosis. The PKHD1 gene is one of the largest genes in the genome with multiple alternative transcripts that have not yet been fully characterized. We sequenced the genomes of an affected foal and 46 control horses to establish a comprehensive list of variants in the critical interval. We identified two missense variants in the PKHD1 gene which were strongly, but not perfectly associated with congenital hepatic fibrosis. We speculate that reduced penetrance and/or potential epistatic interactions with hypothetical modifier genes may explain the imperfect association of the detected PKHD1 variants. Our data thus indicate that horses with congenital hepatic fibrosis represent an interesting large animal model for the liver-restricted subtype of human ARPKD.

An inherited polyneuropathy (PN) observed in Leonberger dogs has clinical similarities to a genetically heterogeneous group of peripheral neuropathies termed Charcot-Marie-Tooth (CMT) disease in humans. The Leonberger disorder is a severe, juvenile-onset, chronic, progressive, and mixed PN, characterized by exercise intolerance, gait abnormalities and muscle atrophy of the pelvic limbs, as well as inspiratory stridor and dyspnea. We mapped a PN locus in Leonbergers to a 250 kb region on canine chromosome 16 (Praw = 1.16×10-10, Pgenome, corrected = 0.006) utilizing a high-density SNP array. Within this interval is the ARHGEF10 gene, a member of the rho family of GTPases known to be involved in neuronal growth and axonal migration, and implicated in human hypomyelination. ARHGEF10 sequencing identified a 10 bp deletion in affected dogs that removes four nucleotides from the 3'-end of exon 17 and six nucleotides from the 5'-end of intron 17 (c.1955_1958+6delCACGGTGAGC). This eliminates the 3'-splice junction of exon 17, creates an alternate splice site immediately downstream in which the processed mRNA contains a frame shift, and generates a premature stop codon predicted to truncate approximately 50% of the protein. Homozygosity for the deletion was highly associated with the severe juvenile-onset PN phenotype in both Leonberger and Saint Bernard dogs. The overall clinical picture of PN in these breeds, and the effects of sex and heterozygosity of the ARHGEF10 deletion, are less clear due to the likely presence of other forms of PN with variable ages of onset and severity of clinical signs. This is the first documented severe polyneuropathy associated with a mutation in ARHGEF10 in any species.

Hereditary footpad hyperkeratosis (HFH) represents a palmoplantar hyperkeratosis, which is inherited as a monogenic autosomal recessive trait in several dog breeds, such as e.g. Kromfohrländer and Irish Terriers. We performed genome-wide association studies (GWAS) in both breeds. In Kromfohrländer we obtained a single strong association signal on chromosome 5 (p(raw) = 1.0×10(-13)) using 13 HFH cases and 29 controls. The association signal replicated in an independent cohort of Irish Terriers with 10 cases and 21 controls (p(raw) = 6.9×10(-10)). The analysis of shared haplotypes among the combined Kromfohrländer and Irish Terrier cases defined a critical interval of 611 kb with 13 predicted genes. We re-sequenced the genome of one affected Kromfohrländer at 23.5× coverage. The comparison of the sequence data with 46 genomes of non-affected dogs from other breeds revealed a single private non-synonymous variant in the critical interval with respect to the reference genome assembly. The variant is a missense variant (c.155G>C) in the FAM83G gene encoding a protein with largely unknown function. It is predicted to change an evolutionary conserved arginine into a proline residue (p.R52P). We genotyped this variant in a larger cohort of dogs and found perfect association with the HFH phenotype. We further studied the clinical and histopathological alterations in the epidermis in vivo. Affected dogs show a moderate to severe orthokeratotic hyperplasia of the palmoplantar epidermis. Thus, our data provide the first evidence that FAM83G has an essential role for maintaining the integrity of the palmoplantar epidermis.

During the summer of 2013 seven Italian Tyrolean Grey calves were born with abnormally short limbs. Detailed clinical and pathological examination revealed similarities to chondrodysplastic dwarfism. Pedigree analysis showed a common founder, assuming autosomal monogenic recessive transmission of the defective allele. A positional cloning approach combining genome wide association and homozygosity mapping identified a single 1.6 Mb genomic region on BTA 6 that was associated with the disease. Whole genome re-sequencing of an affected calf revealed a single candidate causal mutation in the Ellis van Creveld syndrome 2 (EVC2) gene. This gene is known to be associated with chondrodysplastic dwarfism in Japanese Brown cattle, and dwarfism, abnormal nails and teeth, and dysostosis in humans with Ellis-van Creveld syndrome. Sanger sequencing confirmed the presence of a 2 bp deletion in exon 19 (c.2993_2994ACdel) that led to a premature stop codon in the coding sequence of bovine EVC2, and was concordant with the recessive pattern of inheritance in affected and carrier animals. This loss of function mutation confirms the important role of EVC2 in bone development. Genetic testing can now be used to eliminate this form of chondrodysplastic dwarfism from Tyrolean Grey cattle.

The molecular regulation of horn growth in ruminants is still poorly understood. To investigate this process, we collected 1019 hornless (polled) animals from different cattle breeds. High-density SNP genotyping confirmed the presence of two different polled associated haplotypes in Simmental and Holstein cattle co-localized on BTA 1. We refined the critical region of the Simmental polled mutation to 212 kb and identified an overlapping region of 932 kb containing the Holstein polled mutation. Subsequently, whole genome sequencing of polled Simmental and Holstein cows was used to determine polled associated genomic variants. By genotyping larger cohorts of animals with known horn status we found a single perfectly associated insertion/deletion variant in Simmental and other beef cattle confirming the recently published possible Celtic polled mutation. We identified a total of 182 sequence variants as candidate mutations for polledness in Holstein cattle, including an 80 kb genomic duplication and three SNPs reported before. For the first time we showed that hornless cattle with scurs are obligate heterozygous for one of the polled mutations. This is in contrast to published complex inheritance models for the bovine scurs phenotype. Studying differential expression of the annotated genes and loci within the mapped region on BTA 1 revealed a locus (LOC100848215), known in cow and buffalo only, which is higher expressed in fetal tissue of wildtype horn buds compared to tissue of polled fetuses. This implicates that the presence of this long noncoding RNA is a prerequisite for horn bud formation. In addition, both transcripts associated with polledness in goat and sheep (FOXL2 and RXFP2), show an overexpression in horn buds confirming their importance during horn development in cattle.

Ectodermal dysplasias (EDs) are a large and heterogeneous group of hereditary disorders characterized by abnormalities in structures of ectodermal origin. Incontinentia pigmenti (IP) is an ED characterized by skin lesions evolving over time, as well as dental, nail, and ocular abnormalities. Due to X-linked dominant inheritance IP symptoms can only be seen in female individuals while affected males die during development in utero. We observed a family of horses, in which several mares developed signs of a skin disorder reminiscent of human IP. Cutaneous manifestations in affected horses included the development of pruritic, exudative lesions soon after birth. These developed into wart-like lesions and areas of alopecia with occasional wooly hair re-growth. Affected horses also had streaks of darker and lighter coat coloration from birth. The observation that only females were affected together with a high number of spontaneous abortions suggested an X-linked dominant mechanism of transmission. Using next generation sequencing we sequenced the whole genome of one affected mare. We analyzed the sequence data for non-synonymous variants in candidate genes and found a heterozygous nonsense variant in the X-chromosomal IKBKG gene (c.184C>T; p.Arg62*). Mutations in IKBKG were previously reported to cause IP in humans and the homologous p.Arg62* variant has already been observed in a human IP patient. The comparative data thus strongly suggest that this is also the causative variant for the observed IP in horses. To our knowledge this is the first large animal model for IP.

Highland cattle with congenital crop ears have notches of variable size on the tips of both ears. In some cases, cartilage deformation can be seen and occasionally the external ears are shortened. We collected 40 cases and 80 controls across Switzerland. Pedigree data analysis confirmed a monogenic autosomal dominant mode of inheritance with variable expressivity. All affected animals could be traced back to a single common ancestor. A genome-wide association study was performed and the causative mutation was mapped to a 4 Mb interval on bovine chromosome 6. The H6 family homeobox 1 (HMX1) gene was selected as a positional and functional candidate gene. By whole genome re-sequencing of an affected Highland cattle, we detected 6 non-synonymous coding sequence variants and two variants in an ultra-conserved element at the HMX1 locus with respect to the reference genome. Of these 8 variants, only a non-coding 76 bp genomic duplication (g.106720058_106720133dup) located in the conserved region was perfectly associated with crop ears. The identified copy number variation probably results in HMX1 misregulation and possible gain-of-function. Our findings confirm the role of HMX1 during the development of the external ear. As it is sometimes difficult to phenotypically diagnose Highland cattle with slight ear notches, genetic testing can now be used to improve selection against this undesired trait.

Imerslund-Gräsbeck syndrome (IGS) or selective cobalamin malabsorption has been described in humans and dogs. IGS occurs in Border Collies and is inherited as a monogenic autosomal recessive trait in this breed. Using 7 IGS cases and 7 non-affected controls we mapped the causative mutation by genome-wide association and homozygosity mapping to a 3.53 Mb interval on chromosome 2. We re-sequenced the genome of one affected dog at ∼10× coverage and detected 17 non-synonymous variants in the critical interval. Two of these non-synonymous variants were in the cubilin gene (CUBN), which is known to play an essential role in cobalamin uptake from the ileum. We tested these two CUBN variants for association with IGS in larger cohorts of dogs and found that only one of them was perfectly associated with the phenotype. This variant, a single base pair deletion (c.8392delC), is predicted to cause a frameshift and premature stop codon in the CUBN gene. The resulting mutant open reading frame is 821 codons shorter than the wildtype open reading frame (p.Q2798Rfs*3). Interestingly, we observed an additional nonsense mutation in the MRC1 gene encoding the mannose receptor, C type 1, which was in perfect linkage disequilibrium with the CUBN frameshift mutation. Based on our genetic data and the known role of CUBN for cobalamin uptake we conclude that the identified CUBN frameshift mutation is most likely causative for IGS in Border Collies.

We describe a mild form of disproportionate dwarfism in Labrador Retrievers, which is not associated with any obvious health problems such as secondary arthrosis. We designate this phenotype as skeletal dysplasia 2 (SD2). It is inherited as a monogenic autosomal recessive trait with incomplete penetrance primarily in working lines of the Labrador Retriever breed. Using 23 cases and 37 controls we mapped the causative mutation by genome-wide association and homozygosity mapping to a 4.44 Mb interval on chromosome 12. We re-sequenced the genome of one affected dog at 30x coverage and detected 92 non-synonymous variants in the critical interval. Only two of these variants, located in the lymphotoxin A (LTA) and collagen alpha-2(XI) chain gene (COL11A2), respectively, were perfectly associated with the trait. Previously described COL11A2 variants in humans or mice lead to skeletal dysplasias and/or deafness. The dog variant associated with disproportionate dwarfism, COL11A2:c.143G>C or p.R48P, probably has only a minor effect on collagen XI function, which might explain the comparatively mild phenotype seen in our study. The identification of this candidate causative mutation thus widens the known phenotypic spectrum of COL11A2 mutations. We speculate that non-pathogenic COL11A2 variants might even contribute to the heritable variation in height.

The first cases of early-onset progressive polyneuropathy appeared in the Alaskan Malamute population in Norway in the late 1970s. Affected dogs were of both sexes and were ambulatory paraparetic, progressing to non-ambulatory tetraparesis. On neurologic examination, affected dogs displayed predominantly laryngeal paresis, decreased postural reactions, decreased spinal reflexes and muscle atrophy. The disease was considered eradicated through breeding programmes but recently new cases have occurred in the Nordic countries and the USA. The N-myc downstream-regulated gene (NDRG1) is implicated in neuropathies with comparable symptoms or clinical signs both in humans and in Greyhound dogs. This gene was therefore considered a candidate gene for the polyneuropathy in Alaskan Malamutes. The coding sequence of the NDRG1 gene derived from one healthy and one affected Alaskan Malamute revealed a non-synonymous G>T mutation in exon 4 in the affected dog that causes a Gly98Val amino acid substitution. This substitution was categorized to be "probably damaging" to the protein function by PolyPhen2 (score: 1.000). Subsequently, 102 Alaskan Malamutes from the Nordic countries and the USA known to be either affected (n = 22), obligate carriers (n = 7) or healthy (n = 73) were genotyped for the SNP using TaqMan. All affected dogs had the T/T genotype, the obligate carriers had the G/T genotype and the healthy dogs had the G/G genotype except for 13 who had the G/T genotype. A protein alignment showed that residue 98 is conserved in mammals and also that the entire NDRG1 protein is highly conserved (94.7%) in mammals. We conclude that the G>T substitution is most likely the mutation that causes polyneuropathy in Alaskan Malamutes. Our characterization of a novel candidate causative mutation for polyneuropathy offers a new canine model that can provide further insight into pathobiology and therapy of human polyneuropathy. Furthermore, selection against this mutation can now be used to eliminate the disease in Alaskan Malamutes.

Hereditary nasal parakeratosis (HNPK), an inherited monogenic autosomal recessive skin disorder, leads to crusts and fissures on the nasal planum of Labrador Retrievers. We performed a genome-wide association study (GWAS) using 13 HNPK cases and 23 controls. We obtained a single strong association signal on chromosome 2 (p(raw) = 4.4×10⁻¹⁴). The analysis of shared haplotypes among the 13 cases defined a critical interval of 1.6 Mb with 25 predicted genes. We re-sequenced the genome of one case at 38× coverage and detected 3 non-synonymous variants in the critical interval with respect to the reference genome assembly. We genotyped these variants in larger cohorts of dogs and only one was perfectly associated with the HNPK phenotype in a cohort of more than 500 dogs. This candidate causative variant is a missense variant in the SUV39H2 gene encoding a histone 3 lysine 9 (H3K9) methyltransferase, which mediates chromatin silencing. The variant c.972T>G is predicted to change an evolutionary conserved asparagine into a lysine in the catalytically active domain of the enzyme (p.N324K). We further studied the histopathological alterations in the epidermis in vivo. Our data suggest that the HNPK phenotype is not caused by hyperproliferation, but rather delayed terminal differentiation of keratinocytes. Thus, our data provide evidence that SUV39H2 is involved in the epigenetic regulation of keratinocyte differentiation ensuring proper stratification and tight sealing of the mammalian epidermis., Competing Interests: I have read the journal's policy and have the following conflicts: A patent concerning the commercial genetic testing of HNPK in dogs has been filed. TL is the inventor on this patent.

We identified a congenital mechanobullous skin disorder in six calves on a single farm of an endangered German cattle breed in 2010. The condition presented as a large loss of skin distal to the fetlocks and at the mucosa of the muzzle. All affected calves were euthanized on humane grounds due to the severity, extent and progression of the skin and oral lesions. Examination of skin samples under light microscopy revealed detachment of the epidermis from the dermis at the level of the dermo epidermal junction, leading to the diagnosis of a subepidermal bullous dermatosis such as epidermolysis bullosa. The pedigree was consistent with monogenic autosomal recessive inheritance. We localized the causative mutation to an 18 Mb interval on chromosome 22 by homozygosity mapping. The COL7A1 gene encoding collagen type VII alpha 1 is located within this interval and COL7A1 mutations have been shown to cause inherited dystrophic epidermolysis bullosa (DEB) in humans. A SNP in the bovine COL7A1 exon 49 (c.4756C>T) was perfectly associated with the observed disease. The homozygous mutant T/T genotype was exclusively present in affected calves and their parents were heterozygous C/T confirming the assumed recessive mode of inheritance. All known cases and genotyped carriers were related to a single cow, which is supposed to be the founder animal. The mutant T allele was absent in 63 animals from 24 cattle breeds. The identified mutation causes a premature stop codon which leads to a truncated protein representing a complete loss of COL7A1 function (p.R1586*). We thus have identified a candidate causative mutation for this genetic disease using only three cases to unravel its molecular basis. Selection against this mutation can now be used to eliminate the mutant allele from the Rotes Höhenvieh breed.

In the last decade breeders of Romagnola cattle observed an outbreak of a new congenital anomaly. This lethal multi-organ developmental dysplasia is mainly characterized by facial deformities, ascites and hepatic fibrosis. Affected stillborn calves were inbred to a single founder sire suggesting autosomal monogenic recessive inheritance. We localized the causative mutation to a 1.2 Mb interval on BTA 17 by genome-wide association and identical by descent mapping. A solution-based method for targeted DNA capture combined with massively parallel sequencing was used to analyze the entire critical region containing 24 genes. Homozygosity for two non-synonymous coding sequence variants affecting the RNF34 and KDM2B genes was detected by evaluating one affected calf. Here we show that the disease phenotype is associated with a KDM2B missense mutation (c.2503G>A) leading to an amino acid exchange (p.D835N) in an evolutionary strongly conserved domain. In addition, the genetic makeup of three inbred cattle strongly supports the causality of the KDM2B mutation. This report of a naturally-occurring spontaneous mutation of a JmjC domain containing histone demethylase gene provides evidence for their important role in the endo- and mesodermal organ development.

During fetal development neural-crest-derived melanoblasts migrate across the entire body surface and differentiate into melanocytes, the pigment-producing cells. Alterations in this precisely regulated process can lead to white spotting patterns. White spotting patterns in horses are a complex trait with a large phenotypic variance ranging from minimal white markings up to completely white horses. The "splashed white" pattern is primarily characterized by an extremely large blaze, often accompanied by extended white markings at the distal limbs and blue eyes. Some, but not all, splashed white horses are deaf. We analyzed a Quarter Horse family segregating for the splashed white coat color. Genome-wide linkage analysis in 31 horses gave a positive LOD score of 1.6 in a region on chromosome 6 containing the PAX3 gene. However, the linkage data were not in agreement with a monogenic inheritance of a single fully penetrant mutation. We sequenced the PAX3 gene and identified a missense mutation in some, but not all, splashed white Quarter Horses. Genome-wide association analysis indicated a potential second signal near MITF. We therefore sequenced the MITF gene and found a 10 bp insertion in the melanocyte-specific promoter. The MITF promoter variant was present in some splashed white Quarter Horses from the studied family, but also in splashed white horses from other horse breeds. Finally, we identified two additional non-synonymous mutations in the MITF gene in unrelated horses with white spotting phenotypes. Thus, several independent mutations in MITF and PAX3 together with known variants in the EDNRB and KIT genes explain a large proportion of horses with the more extreme white spotting phenotypes., Competing Interests: The authors have declared that no competing interests exist.

Tyrolean Grey cattle represent a local breed with a population size of ∼5000 registered cows. In 2003, a previously unknown neurological disorder was recognized in Tyrolean Grey cattle. The clinical signs of the disorder are similar to those of bovine progressive degenerative myeloencephalopathy (weaver syndrome) in Brown Swiss cattle but occur much earlier in life. The neuropathological investigation of an affected calf showed axonal degeneration in the central nervous system (CNS) and femoral nerve. The pedigrees of the affected calves suggested a monogenic autosomal recessive inheritance. We localized the responsible mutation to a 1.9 Mb interval on chromosome 16 by genome-wide association and haplotype mapping. The MFN2 gene located in this interval encodes mitofusin 2, a mitochondrial membrane protein. A heritable human axonal neuropathy, Charcot-Marie-Tooth disease-2A2 (CMT2A2), is caused by MFN2 mutations. Therefore, we considered MFN2 a positional and functional candidate gene and performed mutation analysis in affected and control Tyrolean Grey cattle. We did not find any non-synonymous variants. However, we identified a perfectly associated silent SNP in the coding region of exon 20 of the MFN2 gene. This SNP is located within a putative exonic splice enhancer (ESE) and the variant allele leads to partial retention of the entire intron 19 and a premature stop codon in the aberrant MFN2 transcript. Thus we have identified a highly unusual splicing defect, where an exonic single base exchange leads to the retention of the preceding intron. This splicing defect represents a potential explanation for the observed degenerative axonopathy. Marker assisted selection can now be used to eliminate degenerative axonopathy from Tyrolean Grey cattle.

Dilated cardiomyopathy (DCM) is a heterogeneous group of heart diseases with a strong genetic background. Currently, many human DCM cases exist where no causative mutation can be identified. DCM also occurs with high prevalence in several large dog breeds. In the Doberman Pinscher a specific DCM form characterized by arrhythmias and/or echocardiographic changes has been intensively studied by veterinary cardiologists. We performed a genome-wide association study in Doberman Pinschers. Using 71 cases and 70 controls collected in Germany we identified a genome-wide significant association to DCM on chromosome 5. We validated the association in an independent cohort collected in the United Kingdom. There is no known DCM candidate gene under the association signal. Therefore, DCM in Doberman Pinschers offers the chance of identifying a novel DCM gene that might also be relevant for human health.

Arachnomelia is a monogenic recessive defect of skeletal development in cattle. The causative mutation was previously mapped to a ∼7 Mb interval on chromosome 5. Here we show that array-based sequence capture and massively parallel sequencing technology, combined with the typical family structure in livestock populations, facilitates the identification of the causative mutation. We re-sequenced the entire critical interval in a healthy partially inbred cow carrying one copy of the critical chromosome segment in its ancestral state and one copy of the same segment with the arachnomelia mutation, and we detected a single heterozygous position. The genetic makeup of several partially inbred cattle provides extremely strong support for the causality of this mutation. The mutation represents a single base insertion leading to a premature stop codon in the coding sequence of the SUOX gene and is perfectly associated with the arachnomelia phenotype. Our findings suggest an important role for sulfite oxidase in bone development., Competing Interests: The authors have declared that no competing interests exist.

The polyneuropathy of juvenile Greyhound show dogs shows clinical similarities to the genetically heterogeneous Charcot-Marie-Tooth (CMT) disease in humans. The pedigrees containing affected dogs suggest monogenic autosomal recessive inheritance and all affected dogs trace back to a single male. Here, we studied the neuropathology of this disease and identified a candidate causative mutation. Peripheral nerve biopsies from affected dogs were examined using semi-thin histology, nerve fibre teasing and electron microscopy. A severe chronic progressive mixed polyneuropathy was observed. Seven affected and 17 related control dogs were genotyped on the 50k canine SNP chip. This allowed us to localize the causative mutation to a 19.5 Mb interval on chromosome 13 by homozygosity mapping. The NDRG1 gene is located within this interval and NDRG1 mutations have been shown to cause hereditary motor and sensory neuropathy-Lom in humans (CMT4D). Therefore, we considered NDRG1 a positional and functional candidate gene and performed mutation analysis in affected and control Greyhounds. A 10 bp deletion in canine NDRG1 exon 15 (c.1080_1089delTCGCCTGGAC) was perfectly associated with the polyneuropathy phenotype of Greyhound show dogs. The deletion causes a frame shift (p.Arg361SerfsX60) which alters several amino acids before a stop codon is encountered. A reduced level of NDRG1 transcript could be detected by RT-PCR. Western blot analysis demonstrated an absence of NDRG1 protein in peripheral nerve biopsy of an affected Greyhound. We thus have identified a candidate causative mutation for polyneuropathy in Greyhounds and identified the first genetically characterized canine CMT model which offers an opportunity to gain further insights into the pathobiology and therapy of human NDRG1 associated CMT disease. Selection against this mutation can now be used to eliminate polyneuropathy from Greyhound show dogs.

Microphthalmia in sheep is an autosomal recessive inherited congenital anomaly found within the Texel breed. It is characterized by extremely small or absent eyes and affected lambs are absolutely blind. For the first time, we use a genome-wide ovine SNP array for positional cloning of a Mendelian trait in sheep. Genotyping 23 cases and 23 controls using Illumina's OvineSNP50 BeadChip allowed us to localize the causative mutation for microphthalmia to a 2.4 Mb interval on sheep chromosome 22 by association and homozygosity mapping. The PITX3 gene is located within this interval and encodes a homeodomain-containing transcription factor involved in vertebrate lens formation. An abnormal development of the lens vesicle was shown to be the primary event in ovine microphthalmia. Therefore, we considered PITX3 a positional and functional candidate gene. An ovine BAC clone was sequenced, and after full-length cDNA cloning the PITX3 gene was annotated. Here we show that the ovine microphthalmia phenotype is perfectly associated with a missense mutation (c.338G>C, p.R113P) in the evolutionary conserved homeodomain of PITX3. Selection against this candidate causative mutation can now be used to eliminate microphthalmia from Texel sheep in production systems. Furthermore, the identification of a naturally occurring PITX3 mutation offers the opportunity to use the Texel as a genetically characterized large animal model for human microphthalmia.

Osteogenesis imperfecta (OI) is a hereditary disease occurring in humans and dogs. It is characterized by extremely fragile bones and teeth. Most human and some canine OI cases are caused by mutations in the COL1A1 and COL1A2 genes encoding the subunits of collagen I. Recently, mutations in the CRTAP and LEPRE1 genes were found to cause some rare forms of human OI. Many OI cases exist where the causative mutation has not yet been found. We investigated Dachshunds with an autosomal recessive form of OI. Genotyping only five affected dogs on the 50 k canine SNP chip allowed us to localize the causative mutation to a 5.82 Mb interval on chromosome 21 by homozygosity mapping. Haplotype analysis of five additional carriers narrowed the interval further down to 4.74 Mb. The SERPINH1 gene is located within this interval and encodes an essential chaperone involved in the correct folding of the collagen triple helix. Therefore, we considered SERPINH1 a positional and functional candidate gene and performed mutation analysis in affected and control Dachshunds. A missense mutation (c.977C>T, p.L326P) located in an evolutionary conserved domain was perfectly associated with the OI phenotype. We thus have identified a candidate causative mutation for OI in Dachshunds and identified a fifth OI gene., Competing Interests: The authors have declared that no competing interests exist.

Scrapie is a transmissible spongiform encephalopathy (TSE) in sheep and goats. In recent years, atypical scrapie cases were identified that differed from classical scrapie in the molecular characteristics of the disease-associated pathological prion protein (PrP(sc)). In this study, we analyze the molecular and neuropathological phenotype of nine Swiss TSE cases in sheep and goats. One sheep was identified as classical scrapie, whereas six sheep, as well as two goats, were classified as atypical scrapie. The latter revealed a uniform electrophoretic mobility pattern of the proteinase K-resistant core fragment of PrP(sc) distinct from classical scrapie regardless of the genotype, the species, and the neuroanatomical structure. Remarkably different types of neuroanatomical PrP(sc) distribution were observed in atypical scrapie cases by both western immunoblotting and immunohistochemistry. Our findings indicate that the biodiversity in atypical scrapie is larger than expected and thus impacts on current sampling and testing strategies in small ruminant TSE surveillance.