Whole genome sequencing reveals a novel deletion variant in the KIT gene in horses with white spotted coat colour phenotypes.

White spotting phenotypes in horses can range in severity from the common white markings up to completely white horses. EDNRB, KIT, MITF, PAX3 and TRPM1 represent known candidate genes for such phenotypes in horses. For the present study, we re-investigated a large horse family segregating a variable white spotting phenotype, for which conventional Sanger sequencing of the candidate genes' individual exons had failed to reveal the causative variant. We obtained whole genome sequence data from an affected horse and specifically searched for structural variants in the known candidate genes. This analysis revealed a heterozygous ~1.9-kb deletion spanning exons 10-13 of the KIT gene (chr3:77,740,239_77,742,136del1898insTATAT). In continuity with previously named equine KIT variants we propose to designate the newly identified deletion variant W22. We had access to 21 horses carrying the W22 allele. Four of them were compound heterozygous W20/W22 and had a completely white phenotype. Our data suggest that W22 represents a true null allele of the KIT gene, whereas the previously identified W20 leads to a partial loss of function. These findings will enable more precise genetic testing for depigmentation phenotypes in horses

Whole Genome Sequencing of Giant Schnauzer Dogs with Progressive Retinal Atrophy Establishes NECAP1 as a Novel Candidate Gene for Retinal Degeneration

Canine progressive retinal atrophies (PRA) are genetically heterogeneous diseases characterized by retinal degeneration and subsequent blindness. PRAs are untreatable and affect multiple dog breeds, significantly impacting welfare. Three out of seven Giant Schnauzer (GS) littermates presented with PRA around four years of age. We sought to identify the causal variant to improve our understanding of the aetiology of this form of PRA and to enable development of a DNA test. Whole genome sequencing of two PRA-affected full-siblings and both unaffected parents was performed. Variants were filtered based on those segregating appropriately for an autosomal recessive disorder and predicted to be deleterious. Successive filtering against 568 canine genomes identified a single nucleotide variant in the gene encoding NECAP endocytosis associated 1 (NECAP1): c.544G>A (p.Gly182Arg). Five thousand one hundred and thirty canids of 175 breeds, 10 cross-breeds and 3 wolves were genotyped for c.544G>A. Only the three PRA-affected GS were homozygous (allele frequency in GS, excluding proband family = 0.015). In addition, we identified heterozygotes belonging to Spitz and Dachshund varieties, demonstrating c.544G>A segregates in other breeds of German origin. This study, in parallel with the known retinal expression and role of NECAP1 in clathrin mediated endocytosis (CME) in synapses, presents NECAP1 as a novel candidate gene for retinal degeneration in dogs and other species

Bayesian model and selection signature analyses reveal risk factors for canine atopic dermatitis

Canine atopic dermatitis is an inflammatory skin disease with clinical similarities to human atopic dermatitis. Several dog breeds are at increased risk for developing this disease but previous genetic associations are poorly defined. To identify additional genetic risk factors for canine atopic dermatitis, we here apply a Bayesian mixture model adapted for mapping complex traits and a cross-population extended haplotype test to search for disease-associated loci and selective sweeps in four dog breeds at risk for atopic dermatitis. We define 15 associated loci and eight candidate regions under selection by comparing cases with controls. One associated locus is syntenic to the major genetic risk locus (Filaggrin locus) in human atopic dermatitis. One selection signal in common type Labrador retriever cases positions across the TBC1D1 gene (body weight) and one signal of selection in working type German shepherd controls overlaps the LRP1B gene (brain), near the KYNU gene (psoriasis). In conclusion, we identify candidate genes, including genes belonging to the same biological pathways across multiple loci, with potential relevance to the pathogenesis of canine atopic dermatitis. The results show genetic similarities between dog and human atopic dermatitis, and future across-species genetic comparisons are hereby further motivated

Autosomal recessive hyposegmentation of granulocytes in Australian Shepherd Dogs indicates a role for LMBR1L in myeloid leukocytes

Pelger-Huët anomaly (PHA) in humans is an autosomal dominant hematological phenotype without major clinical consequences. PHA involves a characteristic hyposegmentation of granulocytes (HG). Human PHA is caused by heterozygous loss of function variants in the LBR gene encoding lamin receptor B. Bi-allelic variants and complete deficiency of LBR cause the much more severe Greenberg skeletal dysplasia which is lethal in utero and characterized by massive skeletal malformation and gross fetal hydrops. HG phenotypes have also been described in domestic animals and homology to human PHA has been claimed in the literature. We studied a litter of Australian Shepherd Dogs with four stillborn puppies in which both parents had an HG phenotype. Linkage analysis excluded LBR as responsible gene for the stillborn puppies. We then investigated the HG phenotype in Australian Shepherd Dogs independently of the prenatal lethality. Genome-wide association mapped the HG locus to chromosome 27 and established an autosomal recessive mode of inheritance. Whole genome sequencing identified a splice site variant in LMBR1L, c.191+1G>A, as most likely causal variant for the HG phenotype. The mutant allele abrogates the expression of the longer X2 isoform but does not affect transcripts encoding the shorter X1 isoform of the LMBR1L protein. The homozygous mutant LMBR1L genotype associated with HG is common in Australian Shepherd Dogs and was found in 39 of 300 genotyped dogs (13%). Our results point to a previously unsuspected function of LMBR1L in the myeloid lineage of leukocytes

Two MC1R loss-of-function alleles in cream-coloured Australian Cattle Dogs and white Huskies.

Loss-of-function variants in the MC1R gene cause recessive red or yellow coat-colour phenotypes in many species. The canine MC1R:c.916C>T (p.Arg306Ter) variant is widespread and found in a homozygous state in many uniformly yellow- or red-coloured dogs. We investigated cream-coloured Australian Cattle Dogs whose coat colour could not be explained by this variant. A genome-wide association study with 10 cream and 123 red Australian Cattle Dogs confirmed that the cream locus indeed maps to MC1R. Whole-genome sequencing of cream dogs revealed a single nucleotide variant within the MITF binding site of the canine MC1R promoter. We propose to designate the mutant alleles at MC1R:c.916C>T as e1 and at the new promoter variant as e2. Both alleles segregate in the Australian Cattle Dog breed. When we considered both alleles in combination, we observed perfect association between the MC1R genotypes and the cream coat colour phenotype in a cohort of 10 cases and 324 control dogs. Analysis of the MC1R transcript levels in an e /e compound heterozygous dog confirmed that the transcript levels of the e2 allele were markedly reduced with respect to the e1 allele. We further report another MC1R loss-of-function allele in Alaskan and Siberian Huskies caused by a 2-bp deletion in the coding sequence, MC1R:c.816_817delCT. We propose to term this allele e3. Huskies that carry two copies of MC1R loss-of-function alleles have a white coat colour

Variants at the ASIP locus contribute to coat color darkening in Nellore cattle

Background: Nellore cattle (Bos indicus) are well-known for their adaptation to warm and humid environments. Hair length and coat color may impact heat tolerance. The Nellore breed has been strongly selected for white coat, but bulls generally exhibit darker hair ranging from light grey to black on the head, neck, hump, and knees. Given the potential contribution of coat color variation to the adaptation of cattle populations to tropical and sub-tropical environments, our aim was to map positional and functional candidate genetic variants associated with darkness of hair coat (DHC) in Nellore bulls. Results: We performed a genome-wide association study (GWAS) for DHC using data from 432 Nellore bulls that were genotyped for more than 777 k single nucleotide polymorphism (SNP) markers. A single major association signal was detected in the vicinity of the agouti signaling protein gene (ASIP). The analysis of whole-genome sequence (WGS) data from 21 bulls revealed functional variants that are associated with DHC, including a structural rearrangement involving ASIP (ASIP-SV1). We further characterized this structural variant using Oxford Nanopore sequencing data from 13 Australian Brahman heifers, which share ancestry with Nellore cattle; we found that this variant originates from a 1155-bp deletion followed by an insertion of a transposable element of more than 150 bp that may impact the recruitment of ASIP non-coding exons. Conclusions: Our results indicate that the variant ASIP sequence causes darker coat pigmentation on specific parts of the body, most likely through a decreased expression of ASIP and consequently an increased production of eumelanin

Identification of the Bovine Arachnomelia Mutation by Massively Parallel Sequencing Implicates Sulfite Oxidase (SUOX) in Bone Development

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

An ADAMTS3 Missense Variant is Associated with Norwich Terrier Upper Airway Syndrome

In flat-faced dog breeds, air resistance caused by skull conformation is believed to be a major determinant of Brachycephalic Obstructive Airway Syndrome (BOAS). The clinical presentation of BOAS is heterogeneous, suggesting determinants independent of skull conformation contribute to airway disease. Norwich Terriers, a mesocephalic breed, are predisposed to Upper Airway Syndrome (UAS), a disease whose pathological features overlap with BOAS. Our health screening clinic examined and scored the airways of 401 Norwich terriers by laryngoscopy. Genome-wide association analyses of UAS-related pathologies revealed a genetic association on canine chromosome 13 (rs9043975, p = 7.79x10-16). Whole genome resequencing was used to identify causal variant(s) within a 414 kb critical interval. This approach highlighted an error in the CanFam3.1 dog assembly, which when resolved, led to the discovery of a c.2786G>A missense variant in exon 20 of the positional candidate gene, ADAM metallopeptidase with thrombospondin type 1 motif 3 (ADAMTS3). In addition to segregating with UAS amongst Norwich Terriers, the ADAMTS3 c.2786G>A risk allele frequency was enriched among the BOAS-susceptible French and (English) Bulldogs. Previous studies indicate that ADAMTS3 loss of function results in lymphoedema. Our results suggest a new paradigm in the understanding of canine upper airway disease aetiology: airway oedema caused by disruption of ADAMTS3 predisposes dogs to respiratory obstruction. These findings will enhance breeding practices and could refine the prognostics of surgical interventions that are often used to treat airway obstruction

Comparison against 186 canid whole genome sequences reveals survival strategies of an ancient clonally transmissible canine tumor.

Canine transmissible venereal tumor (CTVT) is a parasitic cancer clone that has propagated for thousands of years via sexual transfer of malignant cells. Little is understood about the mechanisms that converted an ancient tumor into the world's oldest known continuously propagating somatic cell lineage. We created the largest existing catalog of canine genome-wide variation and compared it against two CTVT genome sequences, thereby separating alleles derived from the founder's genome from somatic drivers of clonal transmissibility. We show that CTVT has undergone continuous adaptation to its transmissible allograft niche, with overlapping mutations at every step of immunosurveillance, particularly self-antigen presentation and apoptosis. We also identified chronologically early somatic mutations in oncogenesis- and immune-related genes that may represent key initiators of clonal transmissibility. Thus, we provide the first insights into the specific genomic aberrations that underlie CTVT's dogged perseverance in canids around the world

Genetic diversity in the modern horse illustrated from genome-wide SNP data

Horses were domesticated from the Eurasian steppes 5,000-6,000 years ago. Since then, the use of horses for transportation, warfare, and agriculture, as well as selection for desired traits and fitness, has resulted in diverse populations distributed across the world, many of which have become or are in the process of becoming formally organized into closed, breeding populations (breeds). This report describes the use of a genome-wide set of autosomal SNPs and 814 horses from 36 breeds to provide the first detailed description of equine breed diversity. F(ST) calculations, parsimony, and distance analysis demonstrated relationships among the breeds that largely reflect geographic origins and known breed histories. Low levels of population divergence were observed between breeds that are relatively early on in the process of breed development, and between those with high levels of within-breed diversity, whether due to large population size, ongoing outcrossing, or large within-breed phenotypic diversity. Populations with low within-breed diversity included those which have experienced population bottlenecks, have been under intense selective pressure, or are closed populations with long breed histories. These results provide new insights into the relationships among and the diversity within breeds of horses. In addition these results will facilitate future genome-wide association studies and investigations into genomic targets of selection