Genomic Characterization of External Morphology Traits in Kelpies Does Not Support Common Ancestry with the Australian Dingo

Genomic Diversity and Runs of Homozygosity in Bernese Mountain Dogs

Bernese mountain dogs are a large dog breed formed in the early 1900s in Switzerland. While originally farm dogs that were used for pulling carts, guarding, and driving cattle, today they are considered multi-purpose companion and family dogs. The breed is predisposed to several complex diseases, such as histiocytic sarcoma, degenerative myelopathy, or hip dysplasia. Using whole-genome sequencing (WGS) data, we assessed the genomic architecture of 33 unrelated dogs from four countries: France, Sweden, Switzerland, and the United States. Analysis of runs of homozygosity (ROH) identified 12,643 ROH with an average length of 2.29 Mb and an average inbreeding coefficient of 0.395. Multidimensional scaling analysis of the genetic relatedness revealed limited clustering of European versus USA dogs, suggesting exchanges of breeding stock between continents. Furthermore, only two mtDNA haplotypes were detected in the 33 studied dogs, both of which are widespread throughout multiple dog breeds. WGS-based ROH analyses revealed several fixed or nearly fixed regions harboring discreet morphological trait-associated as well as disease-associated genetic variants. Several genes involved in the regulation of immune cells were found in the ROH shared by all dogs, which is notable in the context of the breed’s strong predisposition to hematopoietic cancers. High levels of inbreeding and relatedness, strongly exaggerated in the last 30 years, have likely led to the high prevalence of specific genetic disorders in this breed.

Cerebellar Abiotrophy in Australian Working Kelpies Is Associated with Two Major Risk Loci

An autosomal recessive form of inherited cerebellar abiotrophy (CA) that is characterized by a degeneration of Purkinje and granule cells in the cerebellar cortex occurs in the Australian working kelpie dog breed. The clinical signs of CA include ataxia, head tremor, motor in-coordination, wide-based stance, and high-stepping gait. Investigation of clinical and pathological features indicated two closely related diseases with differences in age of onset. A genome-wide association study on 45 CA affected and 290 normal healthy Kelpies identified two significantly associated loci, one on CFA9 and a second on CFA20. Dogs homozygous for the risk haplotype on CFA20 (23 dogs) show clinical signs before ten weeks of age. Missense variants in the sixth exon of disruptor of telomeric silencing 1-like (DOT1Lp.R200Q) and in the only exon of Leucine Rich Repeat And Ig Domain Containing 3 (LINGO3p.R359C), both on CFA20, segregate with the associated risk marker which has incomplete penetrance (42%). Affected dogs homozygous for the risk haplotype on CFA9 have later onset ataxia. A missense variant in exon 5 of Vacuole Membrane Protein 1 (VMP1 p.P160Q) on CFA9 segregates as a fully penetrant Mendelian recessive with later-onset CA. Across mammals, the variety of causative loci so far identified as influencing cerebellar disorders reinforces the complexity of the pathways that contribute to cerebellar development and function, and to the pathophysiological mechanisms that may lead to cerebellar ataxia.

De-novo and genome-wide meta-analyses identify a risk haplotype for congenital sensorineural deafness in Dalmatian dogs

Congenital sensorineural deafness (CSD) has been reported to affect up to 30% of Dalmatian dogs world-wide and while unilaterally deaf dogs can live a close to normal life, dogs suffering bilateral deafness are frequently euthanized. Extreme-white coat patterning as encoded by the gene Melanocyte Inducing Transcription Factor (MITF) has long been postulated as the major risk factor for CSD in the Dalmatian breed. While attempts to identify causative risk variants associated with CSD have been numerous, no genome-wide association study has positively identified MITF as a risk locus for either bilateral or unilateral deafness in the Dalmatian breed to date. In this study, we identified an association with CSD on CFA20 in the vicinity of MITF within Australian Dalmatian dogs. Although not genome-wide significant, the association signal was validated by reanalysing publicly available data and merging the wider data resource with the local data to improve statistical power. The merged data, representing three major global populations of Dalmatian dogs, enabled us to identify a single, well-defined genome-wide significant risk haplotype for CSD. The haplotype was formed by three genome-wide significant associated markers (BICF2G630233852T>C, BICF2G630233861T>C, BICF2G630233888G>A) on CFA20 with 62% of bilaterally deaf dogs homozygous for the risk haplotype (CCA), while 30% of bilaterally deaf and 45% of hearing dogs carried one copy of the risk haplotype. Animals homozygous or heterozygous for the low-risk haplotype were less likely to be unilaterally deaf. While the association between the risk haplotype and deafness is incomplete, animals homozygous for the risk haplotype were 10-times more likely to be bilaterally deaf. Although the underlying causative variants are yet to be discovered, results from this study can now assist with reducing deafness in Dalmatian dogs.

Genome-Wide Association Analysis for Chronic Superficial Keratitis in the Australian Racing Greyhound

Chronic superficial keratitis (CSK) is a progressive inflammatory condition of the eye (cornea) that can cause discomfort and blindness. Differential disease risk across dog breeds strongly suggests that CSK has a genetic basis. In addition to genetic risk, the occurrence of CSK is exacerbated by exposure to ultraviolet light. Genome-wide association analysis considered 109 greyhounds, 70 with CSK and the remainder with normal phenotype at an age over four years. Three co-located variants on CFA18 near the 5' region of the Epidermal Growth Factor Receptor (EGFR) gene were associated with genome-wide significance after multiple-test correction (BICF2P579527, CFA18: 6,068,508, p_raw = 1.77 × 10^-7, p_genome = 0.017; BICF2P1310662, CFA18: 6,077,388, p_raw = 4.09 × 10^-7, p_genome = 0.040; BICF2P160719, CFA18: 6,087,347, p_raw = 4.09 × 10^-7, p_genome = 0.040) (canFam4)). Of the top 10 associated markers, eight were co-located with the significantly associated markers on CFA18. The associated haplotype on CFA18 is protective for the CSK condition. EGFR is known to play a role in corneal healing, where it initiates differentiation and proliferation of epithelial cells that in turn signal the involvement of stromal keratocytes to commence apoptosis. Further validation of the putative functional variants is required prior to their use in genetic testing for breeding programs.

Before Azaria: A Historical Perspective on Dingo Attacks

This paper investigates the origin of the once popular belief in Australian society that wild dingoes do not attack humans. To address this problem, a digital repository of archived newspaper articles and other published texts written between 1788 and 1979 were searched for references to dingoes attacking non-Indigenous people. A total of 52 accounts spanning the period between 1804 and 1928 was identified. A comparison of these historical accounts with the details of modern dingo attacks suggests that at least some of the former are credible. The paper also examined commonly held attitudes towards dingoes in past Australian society based on historical print media articles and other records. Early chroniclers of Australian rural life and culture maintained that dingoes occasionally killed and ate humans out of a predatory motivation. By the early decades of the 20th century, however, an opposing view of this species had emerged: namely, that dingoes were timid animals that continued to pose a danger to livestock, but never to people. This change in the cultural image of dingoes can possibly be linked to more than a century of lethal dingo control efforts greatly reducing the frequency of human–dingo interactions in the most populous parts of the country. This intensive culling may also have expunged the wild genetic pool of dingoes that exhibited bold behaviour around people and/or created a dingo population that was largely wary of humans.

Forensic DNA phenotyping: Canis familiaris breed classification and skeletal phenotype prediction using functionally significant skeletal SNPs and indels

This review highlights a novel application of breed identification and prediction of skeletal traits in forensic investigations using canine DNA evidence. Currently, genotyping methods used for canine breed classification involve the application of highly polymorphic short tandem repeats in addition to larger commercially available SNP arrays. Both applications face technical challenges. An additional approach to breed identification could be through genotyping SNPs and indels that characterise the array of skeletal differences displayed across domestic dog populations. Research has shown that a small number of genetic variants of large effect drive differences in skeletal phenotypes among domestic dog breeds. This feature makes functionally significant canine skeletal variants a cost-effective target for forensic investigators to classify individuals according to their breed. Further analysis of these skeletal variants would enable the prediction of external appearance. To date, functionally significant genes with genetic variants associated with differences in size, bulk, skull shape, ear shape, limb length, digit type, and tail morphology have been uncovered. Recommendations of a cost-effective genotyping method that can be readily designed and applied by forensic investigators have been given. Further advances to improve the field of canine skeletal forensic DNA phenotyping include the refinement of phenotyping methods, further biological validation of the skeletal genetic variants and establishing a publicly available database for storage of allele frequencies of the skeletal genetic variants in the wider domestic dog population.

A large deletion on CFA28 omitting ACSL5 gene is associated with intestinal lipid malabsorption in the Australian Kelpie dog breed

Inborn errors of metabolism are genetic conditions that can disrupt intermediary metabolic pathways and cause defective absorption and metabolism of dietary nutrients. In an Australian Kelpie breeding population, 17 puppies presented with intestinal lipid malabsorption. Juvenile dogs exhibited stunted postnatal growth, steatorrhea, abdominal distension and a wiry coat. Using genome-wide association analysis, an associated locus on CFA28 (P_raw = 2.87E^-06) was discovered and validated in a closely related population (P_raw = 1.75E^-45). A 103.3 kb deletion NC_006610.3CFA28:g.23380074_23483377del, containing genes Acyl-CoA Synthetase Long Chain Family Member 5 (ACSL5) and Zinc Finger DHHC-Type Containing 6 (ZDHHC6), was characterised using whole transcriptomic data. Whole transcriptomic sequencing revealed no expression of ACSL5 and disrupted splicing of ZDHHC6 in jejunal tissue of affected Kelpies. The ACSL5 gene plays a key role in long chain fatty acid absorption, a phenotype similar to that of our affected Kelpies has been observed in a knockout mouse model. A PCR-based diagnostic test was developed and confirmed fully penetrant autosomal recessive mode of inheritance. We conclude the structural variant causing a deletion of the ACSL5 gene is the most likely cause for intestinal lipid malabsorption in the Australian Kelpie.

Ear type in sheep is associated with the MSRB3 locus

Ear morphology is an important determinant of sheep breeds. It includes different variable traits such as ear size and erectness, suggesting a polygenic architecture. Here, we performed a comprehensive genome-wide analysis to identify regions under selection for ear morphology in 515 sheep from 17 breeds fixed for diverse ear phenotypes using 34k SNP genotyping data. GWASs for two ear type traits, size and erectness, revealed a single genome-wide significant association on ovine chromosome 3. The derived marker alleles were enriched in sheep with large and/or floppy ears. The GWAS signal harboured the MSRB3 gene encoding methionine sulphoxide reductase B3, which has already been found to be associated with different ear types in other species. We attempted whole-genome resequencing to identify causal variant(s) within a 1 Mb interval around MSRB3. This experiment excluded major copy number variants in the interval, but failed to identify a compelling candidate causal variant. Fine-mapping suggested that the causal variant for large floppy ears most likely resides in a 175 kb interval downstream of the MSRB3 coding region.

Impact of breeding for coat and spotting patterns on the population structure and genetic diversity of an islander endangered dog breed

First references of the endangered autochthonous Majorcan Ca Mè dog date from the 13th century and enhance its skills and adaptability to the orography. Genealogical historical records were traced back to founders. Founder number in the reference population (397), maximum generations traced, and average number of complete generations were 32, 5, and 2.75, respectively. Structure assessment revealed the existence of subpopulations regarding criteria such as owners (402), breeders (55), coat colours (liver, lemon, black and orange) and spotting patterns (piebald, roan, solid colour, tie or star presence). Average inbreeding (F) within colour groups ranged from 6.3-10.4%, for orange and black populations, respectively. F ranged from 9.43-12.22% for roan patterns and star presence, respectively. Tan point markings showed an F coefficient of 5.85%. The study of genetic diversity revealed a slightly different genetic background between subpopulations. Average coancestry between and within coat colours suggested orange and roan traits could be ascribed to the original nuclei, without omitting the high relationships among other subpopulations. Breeding strategy should select breeding pairs holding a relatedness coefficient below 15%. Hence, coat patterns in dog breeds can help preserving the genetic diversity in endangered dogs, even when these are geographically isolated.