Canine reference genome accuracy impacts variant calling: Lessons learned from investigating embryonic lethal variants.

Deficient homozygosity of a variant maintained in a population suggests that the variant may be embryonic lethal. We examined whole genome sequence data from 675 canids to investigate for variants with missing homozygosity and high predicted impact. Our analysis identified 45 variants, in 32 genes. However, further scrutiny of the sequence reads revealed that all but one of these variants were artifacts of the variant calling process when using CanFam3.1, a widely utilized canine reference genome. We demonstrate that the use of multiple, newer reference genomes could reduce artifacts and lead to more accurate variant identification.

Amino Acids in Dog Nutrition and Health.

The dog has assumed a prominent role in human society. Associated with that status, diet choices for companion dogs have begun to reflect the personal preferences of the owners, with greater emphasis on specialty diets such as organic, vegan/vegetarian, and omission or inclusion of specific ingredients. Despite consumer preferences and many marketing strategies employed, the diets must ensure nutritional adequacy for the dog; if not, health becomes compromised, sometimes severely. The most frequent consideration of consumers and dog food manufacturers is protein source and concentration with a growing emphasis on amino acid composition and bioavailability. Amino acids in general play diverse and critical roles in the dog, with specific amino acids being essential. This review covers what is known regarding amino acids in dog nutrition.

Genetic characterization of Addison's disease in Bearded Collies.

BACKGROUND: Primary hypoadrenocorticism (or Addison's disease, AD) is an autoimmune disease that results in destruction of the adrenal cortex and consequent adrenal insufficiency. The disease has been described in purebred and mixed breed dogs, although some breeds, including the Bearded Collie, are at increased risk for AD. Candidate gene approaches have yielded few associations that appear to be breed-specific. A single other genome-wide association study reported no significant regions of association for AD in Standard Poodles. The present study aimed to identify genomic regions of association for canine AD in Bearded Collies. RESULTS: Our study consists of the first genome-wide association analysis to identify a genome-wide significant region of association with canine AD (CFA18). Peaks of suggestive association were also noted on chromosomes 11, 16 and 29. Logistic regression analysis supported an additive effect of risk genotypes at these smaller effect loci on the probability of disease associated with carrying a risk genotype on CFA18. Potential candidate genes involved in adrenal steroidogenesis, regulation of immune responses and/or inflammation were identified within the associated regions of chromosomes 11 and 16. The gene-poor regions of chromosomes 18 and 29 may, however, harbor regulatory sequences that can modulate gene expression and contribute to disease susceptibility. CONCLUSION: Our findings support the polygenic and complex nature of canine AD and identified a strongly associated locus on CFA18 that, when combined with three other smaller effect loci, was predictive of disease. The results offer progress in the identification of susceptibility loci for canine AD in the Bearded Collie. Further studies are needed to confirm association with the suggested candidate genes and identify actual causative mutations involved with AD susceptibility in this breed.

Gonadectomy effects on the risk of immune disorders in the dog: a retrospective study.

BACKGROUND: Gonadectomy is one of the most common procedures performed on dogs in the United States. Neutering has been shown to reduce the risk for some diseases although recent reports suggest increased prevalence for structural disorders and some neoplasias. The relation between neuter status and autoimmune diseases has not been explored. This study evaluated the prevalence and risk of atopic dermatitis (ATOP), autoimmune hemolytic anemia (AIHA), canine myasthenia gravis (CMG), colitis (COL), hypoadrenocorticism (ADD), hypothyroidism (HYPO), immune-mediated polyarthritis (IMPA), immune-mediated thrombocytopenia (ITP), inflammatory bowel disease (IBD), lupus erythematosus (LUP), and pemphigus complex (PEMC), for intact females, intact males, neutered females, and neutered males. Pyometra (PYO) was evaluated as a control condition. RESULTS: Patient records (90,090) from the William R. Pritchard Veterinary Medical Teaching Hospital at the University of California, Davis from 1995 to 2010 were analyzed in order to determine the risk of immune-mediated disease relative to neuter status in dogs. Neutered dogs had a significantly greater risk of ATOP, AIHA, ADD, HYPO, ITP, and IBD than intact dogs with neutered females being at greater risk than neutered males for all but AIHA and ADD. Neutered females, but not males, had a significantly greater risk of LUP than intact females. Pyometra was a greater risk for intact females. CONCLUSIONS: The data underscore the importance of sex steroids on immune function emphasizing a role of these hormones on tissue self-recognition. Neutering is critically important for population control, reduction of reproductive disorders, and offers convenience for owners. Despite these advantages, the analyses of the present study suggest that neutering is associated with increased risk for certain autoimmune disorders and underscore the need for owners to consult with their veterinary practitioner prior to neutering to evaluate possible benefits and risks associated with such a procedure.

Identification of a common risk haplotype for canine idiopathic epilepsy in the ADAM23 gene.

BACKGROUND: Idiopathic epilepsy is a common neurological disease in human and domestic dogs but relatively few risk genes have been identified to date. The seizure characteristics, including focal and generalised seizures, are similar between the two species, with gene discovery facilitated by the reduced genetic heterogeneity of purebred dogs. We have recently identified a risk locus for idiopathic epilepsy in the Belgian Shepherd breed on a 4.4 megabase region on CFA37. RESULTS: We have expanded a previous study replicating the association with a combined analysis of 157 cases and 179 controls in three additional breeds: Schipperke, Finnish Spitz and Beagle (p(c) = 2.9e-07, p(GWAS) = 1.74E-02). A targeted resequencing of the 4.4 megabase region in twelve Belgian Shepherd cases and twelve controls with opposite haplotypes identified 37 case-specific variants within the ADAM23 gene. Twenty-seven variants were validated in 285 cases and 355 controls from four breeds, resulting in a strong replication of the ADAM23 locus (p(raw) = 2.76e-15) and the identification of a common 28 kb-risk haplotype in all four breeds. Risk haplotype was present in frequencies of 0.49-0.7 in the breeds, suggesting that ADAM23 is a low penetrance risk gene for canine epilepsy. CONCLUSIONS: These results implicate ADAM23 in common canine idiopathic epilepsy, although the causative variant remains yet to be identified. ADAM23 plays a role in synaptic transmission and interacts with known epilepsy genes, LGI1 and LGI2, and should be considered as a candidate gene for human epilepsies.

Single nucleotide polymorphism (SNP) variation of wolves (Canis lupus) in Southeast Alaska and comparison with wolves, dogs, and coyotes in North America.

There is considerable interest in the genetics of wolves (Canis lupus) because of their close relationship to domestic dogs (C. familiaris) and the need for informed conservation and management. This includes wolf populations in Southeast Alaska for which we determined genotypes of 305 wolves at 173662 single nucleotide polymorphism (SNP) loci. After removal of invariant and linked SNP, 123801 SNP were used to quantify genetic differentiation of wolves in Southeast Alaska and wolves, coyotes (C. latrans), and dogs from other areas in North America. There is differentiation of SNP allele frequencies between the species (wolves, coyotes, and dogs), although differentiation is relatively low between some wolf and coyote populations. There are varying levels of differentiation among populations of wolves, including low differentiation of wolves in interior Alaska, British Columbia, and the northern US Rocky Mountains. There is considerable differentiation of SNP allele frequencies of wolves in Southeast Alaska from wolves in other areas. However, wolves in Southeast Alaska are not a genetically homogeneous group and there are comparable levels of genetic differentiation among areas within Southeast Alaska and between Southeast Alaska and other geographic areas. SNP variation and other genetic data are discussed regarding taxonomy and management.

The influence of growth hormone on bone and adipose programming.

In utero growth hormone exposure is associated with distinct immediate growth responses and long term impacts on adult physiological parameters that include obesity, insulin resistance, and bone function. Growth hormone accelerates cellular proliferation in many tissues but is exemplified by increases in the number of cells within the cartilaginous growth plate of bone. In some cases growth hormone also potentiates differentiation as seen in the differentiation of adipocytes that rapidly fill upon withdrawal of growth hormone. Growth hormone provokes these changes either by direct action or through intermediaries such as insulin-like growth factor-I and other downstream effector molecules. The specific mechanism used by growth hormone in programming tissues is not yet fully characterized and likely represents a multipronged approach involving DNA modification, altered adult hormonal milieu, and the development of an augmented stem cell pool capable of future engagement as is seen in adipose accrual. This review summarizes findings of growth hormone's influence on in utero and neonatal cellular and metabolic profiles related to bone and adipose tissue.

Idiopathic Epilepsy Risk Allele Trends in Belgian Tervuren: A Longitudinal Genetic Analysis.

Idiopathic epilepsy (IE) has been known to be inherited in the Belgian Tervuren for many decades. Risk genotypes for IE in this breed have recently been identified on Canis familiaris chromosomes (CFA) 14 and 37. In the current study, the allele frequencies of these loci were analyzed to determine whether dog breeders had employed a purposeful selection against IE, leading to a reduction in risk-associated allele frequency within the breed over time. The allele frequencies of two generational groupings of Belgian Tervuren with and without IE were compared. Allele frequencies for risk-associated alleles on CFA14 were unchanged between 1985 and 2015, whereas those on CFA37 increased during that time in the control population (p < 0.05). In contrast, dogs with IE showed a decrease (p < 0.05) in the IE risk-associated allele frequency at the CFA37 locus. Seizure prevalence in the Belgian Tervuren appears to be increasing. These results suggest that, despite awareness that IE is inherited, selection against IE has not been successful.

Characterization of a dopamine transporter polymorphism and behavior in Belgian Malinois.

BACKGROUND: The Belgian Malinois dog breed (MAL) is frequently used in law enforcement and military environments. Owners have reported seizures and unpredictable behavioral changes including dogs' eyes "glazing over," dogs' lack of response to environmental stimuli, and loss of behavioral inhibition including owner-directed biting behavior. Dogs with severe behavioral changes may be euthanized as they can represent a danger to humans and other dogs. In the dog, the dopamine transporter gene (DAT) contains a 38-base pair variable number tandem repeat (DAT-VNTR); alleles have either one or two copies of the 38-base pair sequence. The objective of this study was to assess frequency of DAT-VNTR alleles, and characterize the association between DAT-VNTR alleles and behavior in MAL and other breeds. RESULTS: In an American sample of 280 dogs comprising 26 breeds, most breeds are predominantly homozygous for the DAT-VNTR two-tandem-repeat allele (2/2). The one-tandem-repeat allele is over-represented in American MAL (AM-MAL) (n = 144), both as heterozygotes (1/2) and homozygotes (1/1). All AM-MAL with reported seizures (n = 5) were 1/1 genotype. For AM-MAL with at least one "1" allele (1/1 or 1/2 genotype, n = 121), owners reported higher levels of attention, increased frequency of episodic aggression, and increased frequency of loss of responsiveness to environmental stimuli. In behavior observations, Belgian Military Working Dogs (MWD) with 1/1 or 1/2 genotypes displayed fewer distracted behaviors and more stress-related behaviors such as lower posture and increased yawning. Handlers' treatment of MWD varied with DAT-VNTR genotype as did dogs' responses to handlers' behavior. For 1/1 or 1/2 genotype MWD, 1) lower posture after the first aversive stimulus given by handlers was associated with poorer obedience performance; 2) increased aversive stimuli during protection exercises were associated with decreased performance; 3) more aversive stimuli during obedience were associated with more aversive stimuli during protection; and 4) handlers used more aversive stimuli in protection compared with obedience exercises. CONCLUSIONS: The single copy allele of DAT-VNTR is associated with owner-reported seizures, loss of responsiveness to environmental stimuli, episodic aggression, and hyper-vigilance in MAL. Behavioral changes are associated with differential treatment by handlers. Findings should be considered preliminary until replicated in a larger sample.

The Histopathological Characteristic of Gastric Carcinoma in the Belgian Tervueren and Groenendael Dog: A Comparison of Two Classification Methods.

Gastric carcinoma is generally considered to be a rare disease in dogs, carrying a grave prognosis. However, in the Tervueren and Groenendael varieties of the Belgian Shepherd dog breed, the disease is highly prevalent. While histopathology is the gold standard for diagnosing gastric carcinoma, there is no general consensus on the methods for histological classification in these cases. Biopsies of a group of 61 dogs with confirmed gastric carcinoma (45 Tervueren and 16 Groenendael) were examined and classified according to World Health Organization (WHO) and Laurén classifications. Kaplan-Meier curves were used to compare survival between the different subtypes and simple and multiple linear regression were used to analyse the association between age of onset and breed variant, sex, neuter status, location of the tumour, inflammation score, and Laurén and WHO classifications. Mean age at diagnosis was significantly different in Groenendael (10.1 ± 2.01) and Tervueren dogs (8.5 ± 1.90). The Laurén classification resulted in 29 (48%) diffuse- and 32 (52%) intestinal-type tumours. Applying the WHO classification resulted in 30 (49%) tubular carcinoma growth patterns and 31 (51%) others. Median survival time was significantly reduced for the diffuse type as compared to the intestinal type according to the Laurén classification, with the same median survival time results for tubular compared to non-tubular subtypes according to the WHO classification (median survival time of 61 vs. 182 days, respectively). Using the WHO and Lauren classification on tumour biopsies may help the practising clinician in the prognostication of gastric carcinoma in Tervueren and Groenendael dogs.

Validation of a Chromosome 14 Risk Haplotype for Idiopathic Epilepsy in the Belgian Shepherd Dog Found to Be Associated with an Insertion in the RAPGEF5 Gene.

An idiopathic epilepsy (IE) risk haplotype on canine chromosome (CFA) 14 has been reported to interact with the CFA37 common risk haplotype in the Belgian shepherd (BS). Additional IE cases and control dogs were genotyped for the risk haplotypes to validate these previous findings. In the new cohort, the interaction between the two regions significantly elevated IE risk. When the haplotypes were analyzed individually, particular haplotypes on both CFA14 (ACTG) and 37 (GG) were associated with elevated IE risk, though only the CFA37 AA was significantly associated (p < 0.003) with reduced risk in the new cohort. However, the CFA14 ACTG risk was statistically significant when the new and previous cohort data were combined. The frequency of the ACTG haplotype was four-fold higher in BS dogs than in other breeds. Whole genome sequence analysis revealed that a 3-base pair predicted disruptive insertion in the RAPGEF5 gene, which is adjacent to the CFA14 risk haplotype. RAPGEF5 is involved in the Wnt-ß-catenin signaling pathway that is crucial for normal brain function. Although this risk variant does not fully predict the likelihood of a BS developing IE, the association with a variant in a candidate gene may provide insight into the genetic control of canine IE.

Handler beliefs affect scent detection dog outcomes.

Our aim was to evaluate how human beliefs affect working dog outcomes in an applied environment. We asked whether beliefs of scent detection dog handlers affect team performance and evaluated relative importance of human versus dog influences on handlers' beliefs. Eighteen drug and/or explosive detection dog/handler teams each completed two sets of four brief search scenarios (conditions). Handlers were falsely told that two conditions contained a paper marking scent location (human influence). Two conditions contained decoy scents (food/toy) to encourage dog interest in a false location (dog influence). Conditions were (1) control; (2) paper marker; (3) decoy scent; and (4) paper marker at decoy scent. No conditions contained drug or explosive scent; any alerting response was incorrect. A repeated measures analysis of variance was used with search condition as the independent variable and number of alerts as the dependent variable. Additional nonparametric tests compared human and dog influence. There were 225 incorrect responses, with no differences in mean responses across conditions. Response patterns differed by condition. There were more correct (no alert responses) searches in conditions without markers. Within marked conditions, handlers reported that dogs alerted more at marked locations than other locations. Handlers' beliefs that scent was present potentiated handler identification of detection dog alerts. Human more than dog influences affected alert locations. This confirms that handler beliefs affect outcomes of scent detection dog deployments.

Genome-Wide Association Studies Reveal Susceptibility Loci for Noninfectious Claw Lesions in Holstein Dairy Cattle.

Sole ulcers (SUs) and white line disease (WLD) are two common noninfectious claw lesions (NICL) that arise due to a compromised horn production and are frequent causes of lameness in dairy cattle, imposing welfare and profitability concerns. Low to moderate heritability estimates of SU and WLD susceptibility indicate that genetic selection could reduce their prevalence. To identify the susceptibility loci for SU, WLD, SU and/or WLD, and any type of noninfectious claw lesion, genome-wide association studies (GWAS) were performed using generalized linear mixed model (GLMM) regression, chunk-based association testing (CBAT), and a random forest (RF) approach. Cows from five commercial dairies in California were classified as controls having no lameness records and ≥6 years old (n = 102) or cases having SU (n = 152), WLD (n = 117), SU and/or WLD (SU + WLD, n = 198), or any type of noninfectious claw lesion (n = 217). The top single nucleotide polymorphisms (SNPs) were defined as those passing the Bonferroni-corrected suggestive and significance thresholds in the GLMM analysis or those that a validated RF model considered important. Effects of the top SNPs were quantified using Bayesian estimation. Linkage disequilibrium (LD) blocks defined by the top SNPs were explored for candidate genes and previously identified, functionally relevant quantitative trait loci. The GLMM and CBAT approaches revealed the same regions of association on BTA8 for SU and BTA13 common to WLD, SU + WLD, and NICL. These SNPs had effects significantly different from zero, and the LD blocks they defined explained a significant amount of phenotypic variance for each dataset (6.1-8.1%, p < 0.05), indicating the small but notable contribution of these regions to susceptibility. These regions contained candidate genes involved in wound healing, skin lesions, bone growth and mineralization, adipose tissue, and keratinization. The LD block defined by the most significant SNP on BTA8 for SU included a SNP previously associated with SU. The RF models were overfitted, indicating that the SNP effects were very small, thereby preventing meaningful interpretation of SNPs and any downstream analyses. These findings suggested that variants associated with various physiological systems may contribute to susceptibility for NICL, demonstrating the complexity of genetic predisposition.

Pleiotropic Loci Associated With Foot Disorders and Common Periparturient Diseases in Holstein Cattle.

Lameness is an animal welfare issue that incurs substantial financial and environmental costs. This condition is commonly caused by digital dermatitis (DD), sole ulcers (SU), and white line disease (WLD). Susceptibility to these three foot disorders is due in part to genetics, indicating that genomic selection against these foot lesions can be used to reduce lameness prevalence. It is unclear whether selection against foot lesions will lead to increased susceptibility to other common diseases such as mastitis and metritis. Thus, the aim of this study was to determine the genetic correlation between causes of lameness and other common health disorders to identify loci contributing to the correlation. Genetic correlation estimates between SU and DD and between SU and WLD were significantly different from zero (p < 0.05), whereas estimates between DD and mastitis, DD and milk fever, and SU and metritis were suggestive (p < 0.1). All five of these genetic correlation estimates were positive. Two-trait genome-wide association studies (GWAS) for each of these five pairs of traits revealed common regions of association on BTA1 and BTA8 for pairs that included DD or SU as one of the traits, respectively. Other regions of association were unique to the pair of traits and not observed in GWAS for other pairs of traits. The positive genetic correlation estimates between foot disorders and other health disorders imply that selection against foot disorders may also decrease susceptibility to other health disorders. Linkage disequilibrium blocks defined around significant and suggestive SNPs from the two-trait GWAS included genes and QTL that were functionally relevant, supporting that these regions included pleiotropic loci.

Enhancing student scientific literacy through participation in citizen science focused on companion animal behavior.

The Covid-19 pandemic served as the impetus to implement activities designed to engage students in the remote instructional environment while simultaneously developing scientific literacy skills. In a high enrollment general education animal science course, numerous activities were designed to improve scientific literacy. These included specifically developed videos covering strategies for reading published science literature, the utilization of topically relevant scientific articles that captured student interest, and engaging students in a citizen science exercise on whether dogs align themselves to the Earth's magnetic field during excretion behavior. Employing pre- and post-self-perception surveys coupled with tasking students to apply their scientific literacy skills in an assessment scenario demonstrated that students' self-perception of their scientific literacy improved 30% (P < 0.05) with approximately 80% of students accurately applying their literacy skills. The citizen science study on excretory behavior was modeled on previously published findings thereby providing an opportunity to validate the published work which had indicated that dogs align their bodies in a North-South axis during excretion. The present study did not demonstrate preferential alignment to any geomagnetic orientation which emphasized to the students the need for scientific replication. Inclusion of simple activities that were relevant to students' daily lives, and providing interpretive context for those activities, resulted in improved self-perceived scientific literacy.

Whole Genome Sequencing Reveals Multiple Linked Genetic Variants on Canine Chromosome 12 Associated with Risk for Symmetrical Lupoid Onychodystrophy (SLO) in the Bearded Collie.

In dogs, symmetrical lupoid onychodystrophy (SLO) results in nail loss and an abnormal regrowth of the claws. In Bearded Collies, an autoimmune nature has been suggested because certain dog leukocyte antigen (DLA) class II haplotypes are associated with the condition. A genome-wide association study of the Bearded Collie revealed two regions of association that conferred risk for disease: one on canine chromosome (CFA) 12 that encompasses the DLA genes, and one on CFA17. Case-control association was employed on whole genome sequencing data to uncover putative causative variants in SLO within the CFA12 and CFA17 associated regions. Genotype imputation was then employed to refine variants of interest. Although no SLO-associated protein-coding variants were identified on CFA17, multiple variants, many with predicted damaging effects, were identified within potential candidate genes on CFA12. Furthermore, many potentially damaging alleles were fully correlated with the presence of DLA class II risk haplotypes for SLO, suggesting that the variants may reflect DLA class II haplotype association with disease or vice versa. Strong linkage disequilibrium in the region precluded the ability to isolate and assess the individual or combined effect of variants on disease development. Nonetheless, all were predictive of risk for SLO and, with judicious assessment, their application in selective breeding may prove useful to reduce the incidence of SLO in the breed.

Genome-wide linkage scan for loci associated with epilepsy in Belgian shepherd dogs.

BACKGROUND: Idiopathic epilepsy in the Belgian shepherd dog is known to have a substantial genetic component. The objective of this study was to identify genomic regions associated with the expression of generalized seizures in the Belgian Tervuren and Sheepdog. RESULTS: DNA from 366 dogs, of which 74 were classified as epileptic, representing two extended families were subjected to a genome-wide linkage scan using 410 microsatellite markers yielding informative coverage averaging 5.95 +/- 0.21 Mb. Though previous studies based on pedigree analyses proposed a major gene of influence, the present study demonstrated the trait to be highly polygenic. Studies of complex disorders in humans indicate that a liberal composite evaluation of genetic linkage is needed to identify underlying quantitative trait loci (QTLs). Four chromosomes yielded tentative linkage based upon LOD scores in excess of 1.0. Possible QTLs within these regions were supported also by analyses of multipoint linkage, allele frequency, TDT, and transmission of haplotype blocks. CONCLUSIONS: Taken together the data tentatively indicate six QTLs, three on CFA 2, and one on each of CFA 6, 12, and 37, that support fine mapping for mutations associated with epilepsy in the Belgian shepherd. The study also underscores the complexity of genomic linkage studies for polygenic disorders.

Owner reports of attention, activity, and impulsivity in dogs: a replication study.

BACKGROUND: When developing behaviour measurement tools that use third party assessments, such as parent report, it is important to demonstrate reliability of resulting scales through replication using novel cohorts. The domestic dog has been suggested as a model to investigate normal variation in attention, hyperactivity, and impulsive behaviours impaired in Attention Deficit Hyperactive Disorder (ADHD). The human ADHD Rating Scale, modified for dogs and using owner-directed surveys, was applied in a European sample. We asked whether findings would be replicated utilizing an Internet survey in a novel sample, where unassisted survey completion, participant attitudes and breeds might affect previous findings. METHODS: Using a slightly modified version of the prior survey, we collected responses (n = 1030, 118 breeds representing 7 breed groups) primarily in the United States and Canada. This study was conducted using an Internet survey mechanism. RESULTS: Reliability analyses confirmed two scales previously identified for dogs (inattention [IA], hyperactivity-impulsivity [HA-IM]). Models including age, training status, and breed group accounted for very little variance in subscales, with no effect of gender. CONCLUSIONS: The factor invariance demonstrated in these findings confirms that owner report, using this modified human questionnaire, provides dog scores according to "inattention" and "hyperactivity-impulsivity" axes. Further characterization of naturally occurring variability of attention, activity, and impulsivity in domestic dogs may provide insight into genetic backgrounds underlying behaviours impaired in attention and associated disorders.