Complex disease and phenotype mapping in the domestic dog

Multiple ancestries and shared gene flow among modern livestock guarding dogs

Livestock guarding dogs (LGDs) have been used to protect livestock for millennia. While previous works suggested a single origin of modern LGDs, the degree and source of shared ancestry have not been tested. To address this, we generated genome-wide SNP data from 304 LGDs and combined it with public genomic data from 2,183 modern and 22 ancient dogs. Our findings reveal shared ancestry and extensive gene flow among modern LGD breeds which we attribute to historical livestock migrations. Additionally, admixture between LGDs and free-ranging dogs argues against reproductive isolation as a core mechanism for maintaining the specialized skills of LGDs. Finally, we identify two lineages within modern LGDs and uncover multiple ancestries tracing back to distinct Eurasian ancient dogs, concordant with the absence of a single ancestor. Overall, our work explores the complex evolutionary history of LGDs, offering valuable insights into how human and livestock co-migrations shaped this functional group.

Genome-wide association study for morphological and hunting-behavior traits in Braque Français Type Pyrénées dogs: A preliminary study

High-throughput genotyping offers great potential to increase our understanding of the genomic basis of canid variation. Braque Français Type Pyrénées (BRA) are smart, agile, and friendly dogs originally developed for tracking, hunting, and retrieving feathered game. On a population of 44 unrelated BRA dogs, single nucleotide polymorphism (SNP) genotype data from the CanineHD Whole-Genome Genotyping BeadChip and evaluation scores for 12 traits related to morphology and hunting performance were available. After quality filtering, 95,859 SNPs on the 38 dog autosomes (CFA) were retained. Phenotypic scores were expressed on a scale from 1 (worst) to 6 (best) and were mostly poorly to moderately correlated except for some morphological traits (e.g. r = 0.81 between the conformation of the head and that of the eye). From GWAS, a total of 378 SNP-phenotype associations with posterior odds of association > 1 have been detected. The strongest associations were found for the eye conformation, for the skull/muzzle ratio, and for connection to the hunter. These included both new and previously identified markers and genes potentially involved with type and behavior traits in BRA. Six of the significant markers mapped within SETDB2, a gene known to be related to pointing behavior in dogs. These results advance our understanding of the genetic basis for morphology and hunting behavior in dogs and identify new variants which are potential targets for further research.

Large-scale genomic analysis of the domestic dog informs biological discovery

Recent advances in genomics, coupled with a unique population structure and remarkable levels of variation, have propelled the domestic dog to new levels as a system for understanding fundamental principles in mammalian biology. Central to this advance are more than 350 recognized breeds, each a closed population that has undergone selection for unique features. Genetic variation in the domestic dog is particularly well characterized compared with other domestic mammals, with almost 3000 high-coverage genomes publicly available. Importantly, as the number of sequenced genomes increases, new avenues for analysis are becoming available. Herein, we discuss recent discoveries in canine genomics regarding behavior, morphology, and disease susceptibility. We explore the limitations of current data sets for variant interpretation, tradeoffs between sequencing strategies, and the burgeoning role of long-read genomes for capturing structural variants. In addition, we consider how large-scale collections of whole-genome sequence data drive rare variant discovery and assess the geographic distribution of canine diversity, which identifies Asia as a major source of missing variation. Finally, we review recent comparative genomic analyses that will facilitate annotation of the noncoding genome in dogs.

Imputation of ancient canid genomes reveals inbreeding history over the past 10,000 years

The multi-millenia long history between dogs and humans has placed them at the forefront of archeological and genomic research. Despite ongoing efforts including the analysis of ancient dog and wolf genomes, many questions remain regarding their geographic and temporal origins, and the microevolutionary processes that led to the diversity of breeds today. Although ancient genomes provide valuable information, their use is hindered by low depth of coverage and post-mortem damage, which inhibits confident genotype calling. In the present study, we assess how genotype imputation of ancient dog and wolf genomes, utilising a large reference panel, can improve the resolution provided by ancient datasets. Imputation accuracy was evaluated by down-sampling high coverage dog and wolf genomes to 0.05-2x coverage and comparing concordance between imputed and high coverage genotypes. We measured the impact of imputation on principal component analyses and runs of homozygosity. Our findings show high (R2>0.9) imputation accuracy for dogs with coverage as low as 0.5x and for wolves as low as 1.0x. We then imputed a dataset of 90 ancient dog and wolf genomes, to assess changes in inbreeding during the last 10,000 years of dog evolution. Ancient dog and wolf populations generally exhibited lower inbreeding levels than present-day individuals. Interestingly, regions with low ROH density maintained across ancient and present-day samples were significantly associated with genes related to olfaction and immune response. Our study indicates that imputing ancient canine genomes is a viable strategy that allows for the use of analytical methods previously limited to high-quality genetic data.

Exploring the impact of epidemiological and clinical factors on the progression of canine leishmaniosis by statistical and whole genome analyses: from breed predisposition to comorbidities

Canine leishmaniosis (CanL), caused by Leishmania infantum, is a complex disease of growing importance in Europe. Clinical manifestations result from the down-modulation of the host immune response through multiple host-parasite interactions. Although several factors might influence CanL progression, this is the first known study evaluating risk factors for its different clinical stages in a large referral hospital population (n = 35.669) from an endemic area, over a 20 year period. Genome-wide scans for selection signatures were also conducted to explore the genomic component of clinical susceptibility to L. infantum infection. The prevalence of CanL was 3.2% (16.7% stage I; 43.6% stage II; 32.1% stage III; 7.6% stage IV). Dog breed (crossbreed), bodyweight (<10 kg), living conditions (indoors), regular deworming treatment, and being vaccinated against Leishmania significantly decreased the transmission risk and the risk for developing severe clinical forms. Conversely, the detection of comorbidities was associated with advanced clinical forms, particularly chronic kidney disease, neoplasia, cryptorchidism, infectious tracheobronchitis and urate urolithiasis, although those did not impact the clinical outcome. Significant associations between an increased risk of severe clinical stages and findings in the anamnesis (renal or skin-related manifestations) and physical examination (ocular findings) were also detected, highlighting their diagnostic value in referred cases of CanL. Sixteen breeds were found to be significantly more susceptible to developing severe stages of leishmaniosis (e.g. Great Dane, Rottweiler, English Springer Spaniel, Boxer, American Staffordshire Terrier, Golden Retriever), while 20 breeds displayed a clinical resistantance phenotype and, thus, are more likely to mount an efficient immune response against L. infantum (e.g. Pointer, Samoyed, Spanish Mastiff, Spanish Greyhound, English Setter, Siberian Husky). Genomic analyses of these breeds retrieved 12 regions under selection, 63 candidate genes and pinpointed multiple biological pathways such as the IRE1 branch of the unfolded protein response, which could play a critical role in clinical susceptibility to L. infantum infection.

A variant in the 5'UTR of ERBB4 is associated with lifespan in Golden Retrievers

Genome-wide association studies (GWAS) in long-lived human populations have led to identification of variants associated with Alzheimer's disease and cardiovascular disease, the latter being the most common cause of mortality in people worldwide. In contrast, naturally occurring cancer represents the leading cause of death in pet dogs, and specific breeds like the Golden Retriever (GR) carry up to a 65% cancer-related death rate. We hypothesized that GWAS of long-lived GRs might lead to the identification of genetic variants capable of modifying longevity within this cancer-predisposed breed. A GWAS was performed comparing GR dogs ≥ 14 years to dogs dying prior to age 12 which revealed a significant association to ERBB4, the only member of the epidermal growth factor receptor family capable of serving as both a tumor suppressor gene and an oncogene. No coding variants were identified, however, distinct haplotypes in the 5'UTR were associated with reduced lifespan in two separate populations of GR dogs. When all GR dogs were analyzed together (n = 304), the presence of haplotype 3 was associated with shorter survival (11.8 years vs. 12.8 years, p = 0.024). GRs homozygous for haplotype 3 had the shortest survival, and GRs homozygous for haplotype 1 had the longest survival (11.6 years vs. 13.5 years, p = 0.0008). Sub-analyses revealed that the difference in lifespan for GRs carrying at least 1 copy of haplotype 3 was specific to female dogs (p = 0.009), whereas survival remained significantly different in both male and female GRs homozygous for haplotype 1 or haplotype 3 (p = 0.026 and p = 0.009, respectively). Taken together, these findings implicate a potential role for ERBB4 in GR longevity and provide evidence that within-breed canine lifespan studies could serve as a mechanism to identify favorable or disease-modifying variants important to the axis of aging and cancer.

DNA methylation and chromatin accessibility predict age in the domestic dog

Across mammals, the epigenome is highly predictive of chronological age. These "epigenetic clocks," most of which have been built using DNA methylation (DNAm) profiles, have gained traction as biomarkers of aging and organismal health. While the ability of DNAm to predict chronological age has been repeatedly demonstrated, the ability of other epigenetic features to predict age remains unclear. Here, we use two types of epigenetic information-DNAm, and chromatin accessibility as measured by ATAC-seq-to develop age predictors in peripheral blood mononuclear cells sampled from a population of domesticated dogs. We measured DNAm and ATAC-seq profiles for 71 dogs, building separate predictive clocks from each, as well as the combined dataset. We also use fluorescence-assisted cell sorting to quantify major lymphoid populations for each sample. We found that chromatin accessibility can accurately predict chronological age (R2 ATAC = 26%), though less accurately than the DNAm clock (R2 DNAm = 33%), and the clock built from the combined datasets was comparable to both (R2 combined = 29%). We also observed various populations of CD62L+ T cells significantly correlated with dog age. Finally, we found that all three clocks selected features that were in or near at least two protein-coding genes: BAIAP2 and SCARF2, both previously implicated in processes related to cognitive or neurological impairment. Taken together, these results highlight the potential of chromatin accessibility as a complementary epigenetic resource for modeling and investigating biologic age.

Mass spectrometry-based proteomics profiling of dogs with and without oral diseases: a systematic review

BACKGROUND

Understanding the distinct proteomics profiles in dogs' oral biofluids enhances diagnostic and therapeutic insights for canine oral diseases, fostering cross-species translational research in dentistry and medicine. This study aimed to conduct a systematic review to investigate the similarities and differences between the oral biofluids' proteomics profile of dogs with and without oral diseases.

METHODS

PubMed, Web of Science, and Scopus were searched with no restrictions on publication language or year to address the following focused question: "What is the proteome signature of healthy versus diseased (oral) dogs' biofluids?" Gene Ontology enrichment and the Kyoto Encyclopedia of Genes and Genomes pathway analyses of the most abundant proteins were performed. Moreover, protein-protein interaction analysis was conducted. The risk of bias (RoB) among the included studies was assessed using the Joanna Briggs Institute (JBI) Critical Appraisal Checklist for Studies Reporting Prevalence Data.

RESULTS

In healthy dogs, the proteomic analysis identified 5,451 proteins, with 137 being the most abundant, predominantly associated with 'innate immune response'. Dogs with oral diseases displayed 6,470 proteins, with distinct associations: 'defense response to bacterium' (periodontal diseases), 'negative regulation of transcription' (dental calculus), and 'positive regulation of transcription' (oral tumors). Clustering revealed significant protein clusters in each case, emphasizing the diverse molecular profiles in health and oral diseases. Only six studies were provided to the JBI tool, as they encompassed case-control evaluations that compared healthy dogs to dogs with oral disease(s). All included studies were found to have low RoB (high quality).

CONCLUSION

Significant differences in the proteomics profiles of oral biofluids between dogs with and without oral diseases were found. The synergy of animal proteomics and bioinformatics offers a promising avenue for cross-species research, despite persistent challenges in result validation.

Genome wide association study in Swedish Labrador retrievers identifies genetic loci associated with hip dysplasia and body weight

Genome wide association studies (GWAS) have been utilized to identify genetic risk loci associated with both simple and complex inherited disorders. Here, we performed a GWAS in Labrador retrievers to identify genetic loci associated with hip dysplasia and body weight. Hip dysplasia scores were available for 209 genotyped dogs. We identified a significantly associated locus for hip dysplasia on chromosome 24, with three equally associated SNPs (p = 4.3 × 10-7) in complete linkage disequilibrium located within NDRG3, a gene which in humans has been shown to be differentially expressed in osteoarthritic joint cartilage. Body weight, available for 85 female dogs, was used as phenotype for a second analysis. We identified two significantly associated loci on chromosome 10 (p = 4.5 × 10-7) and chromosome 31 (p = 2.5 × 10-6). The most associated SNPs within these loci were located within the introns of the PRKCE and CADM2 genes, respectively. PRKCE has been shown to play a role in regulation of adipogenesis whilst CADM2 has been associated with body weight in multiple human GWAS. In summary, we identified credible candidate loci explaining part of the genetic inheritance for hip dysplasia and body weight in Labrador retrievers with strong candidate genes in each locus previously implicated in the phenotypes investigated.

Genomic Analysis of Lymphoma Risk in Bullmastiff Dogs

Lymphoma is the most common haematological malignancy affecting dogs and has a high incidence in the Bullmastiff breed. The aim of this study was to identify risk loci predisposing this breed to the disease. The average age of lymphoma diagnosis in 55 cases was less than 6 years, similar to the median age of 64 cases from our clinical and pathology databases. When fine-scale population structure was explored using NETVIEW, cases were distributed throughout an extended pedigree. When genotyped cases (n = 49) and dogs from the control group (n = 281) were compared in a genome-wide association analysis of lymphoma risk, the most prominent associated regions were detected on CFA13 and CFA33. The top SNPs in a 5.4 Mb region on CFA13 were significant at a chromosome-wide level, and the region was fine-mapped to ~1.2 Mb (CFA13: 25.2-26.4 Mb; CanFam3.1) with four potential functional candidates, including the MYC proto-oncogene bHLH transcription factor (MYC) and a region syntenic with the human and mouse lncRNA Pvt1 oncogene (PVT1). A 380 Kb associated region at CFA33: 7.7-8.1 Mb contained the coding sequence for SUMO specific peptidase7 (SENP7) and NFK inhibitor zeta (NFKBIZ) genes. These genes have annotations related to cancer, amongst others, and both have functional links to MYC regulation. Genomic signatures identified in lymphoma cases suggest that increased risk contributed by the regions identified by GWAS may complement a complex predisposing genetic background.

The TiHoCL panel for canine lymphoma: a feasibility study integrating functional genomics and network biology approaches for comparative oncology targeted NGS panel design

Targeted next-generation sequencing (NGS) enables the identification of genomic variants in cancer patients with high sensitivity at relatively low costs, and has thus opened the era to personalized human oncology. Veterinary medicine tends to adopt new technologies at a slower pace compared to human medicine due to lower funding, nonetheless it embraces technological advancements over time. Hence, it is reasonable to assume that targeted NGS will be incorporated into routine veterinary practice in the foreseeable future. Many animal diseases have well-researched human counterparts and hence, insights gained from the latter might, in principle, be harnessed to elucidate the former. Here, we present the TiHoCL targeted NGS panel as a proof of concept, exemplifying how functional genomics and network approaches can be effectively used to leverage the wealth of information available for human diseases in the development of targeted sequencing panels for veterinary medicine. Specifically, the TiHoCL targeted NGS panel is a molecular tool for characterizing and stratifying canine lymphoma (CL) patients designed based on human non-Hodgkin lymphoma (NHL) research outputs. While various single nucleotide polymorphisms (SNPs) have been associated with high risk of developing NHL, poor prognosis and resistance to treatment in NHL patients, little is known about the genetics of CL. Thus, the ~100 SNPs featured in the TiHoCL targeted NGS panel were selected using functional genomics and network approaches following a literature and database search that shielded ~500 SNPs associated with, in nearly all cases, human hematologic malignancies. The TiHoCL targeted NGS panel underwent technical validation and preliminary functional assessment by sequencing DNA samples isolated from blood of 29 lymphoma dogs using an Ion Torrent™ PGM System achieving good sequencing run metrics. Our design framework holds new possibilities for the design of similar molecular tools applied to other diseases for which limited knowledge is available and will improve drug target discovery and patient care.

Gene purging and the evolution of Neoave metabolism and longevity

Maintenance of the proteasome requires oxidative phosphorylation (ATP) and mitigation of oxidative damage, in an increasingly dysfunctional relationship with aging. SLC3A2 plays a role on both sides of this dichotomy as an adaptor to SLC7A5, a transporter of branched-chain amino acids (BCAA: Leu, Ile, Val), and to SLC7A11, a cystine importer supplying cysteine to the synthesis of the antioxidant glutathione. Endurance in mammalian muscle depends in part on oxidation of BCAA; however, elevated serum levels are associated with insulin resistance and shortened lifespans. Intriguingly, the evolution of modern birds (Neoaves) has entailed the purging of genes including SLC3A2, SLC7A5, -7, -8, -10, and SLC1A4, -5, largely removing BCAA exchangers and their interacting Na+/Gln symporters in pursuit of improved energetics. Additional gene purging included mitochondrial BCAA aminotransferase (BCAT2), pointing to reduced oxidation of BCAA and increased hepatic conversion to triglycerides and glucose. Fat deposits are anhydrous and highly reduced, maximizing the fuel/weight ratio for prolonged flight, but fat accumulation in muscle cells of aging humans contributes to inflammation and senescence. Duplications of the bidirectional α-ketoacid transporters SLC16A3, SLC16A7, the cystine transporters SLC7A9, SLC7A11, and N-glycan branching enzymes MGAT4B, MGAT4C in Neoaves suggests a shift to the transport of deaminated essential amino acid, and stronger mitigation of oxidative stress supported by the galectin lattice. We suggest that Alfred Lotka's theory of natural selection as a maximum power organizer (PNAS 8:151,1922) made an unusually large contribution to Neoave evolution. Further molecular analysis of Neoaves may reveal novel rewiring with applications for human health and longevity.

BRASS: Permutation methods for binary traits in genetic association studies with structured samples

In genetic association analysis of complex traits, permutation testing can be a valuable tool for assessing significance when the distribution of the test statistic is unknown or not well-approximated. This commonly arises, e.g, in tests of gene-set, pathway or genome-wide significance, or when the statistic is formed by machine learning or data adaptive methods. Existing applications include eQTL mapping, association testing with rare variants, inclusion of admixed individuals in genetic association analysis, and epistasis detection among many others. For genetic association testing in samples with population structure and/or relatedness, use of naive permutation can lead to inflated type 1 error. To address this in quantitative traits, the MVNpermute method was developed. However, for association mapping of a binary trait, the relationship between the mean and variance makes both naive permutation and the MVNpermute method invalid. We propose BRASS, a permutation method for binary traits, for use in association mapping in structured samples. In addition to modeling structure in the sample, BRASS allows for covariates, ascertainment and simultaneous testing of multiple markers, and it accommodates a wide range of test statistics. In simulation studies, we compare BRASS to other permutation and resampling-based methods in a range of scenarios that include population structure, familial relatedness, ascertainment and phenotype model misspecification. In these settings, we demonstrate the superior control of type 1 error by BRASS compared to the other 6 methods considered. We apply BRASS to assess genome-wide significance for association analyses in domestic dog for elbow dysplasia (ED) and idiopathic epilepsy (IE). For both traits we detect previously identified associations, and in addition, for ED, we detect significant association with a SNP on chromosome 35 that was not detected by previous analyses, demonstrating the potential of the method.

Identification of genetic variants associated with anterior cruciate ligament rupture and AKC standard coat color in the Labrador Retriever

Canine anterior cruciate ligament (ACL) rupture is a common complex disease. Prevalence of ACL rupture is breed dependent. In an epidemiological study, yellow coat color was associated with increased risk of ACL rupture in the Labrador Retriever. ACL rupture risk variants may be linked to coat color through genetic selection or through linkage with coat color genes. To investigate these associations, Labrador Retrievers were phenotyped as ACL rupture case or controls and for coat color and were single nucleotide polymorphism (SNP) genotyped. After filtering, ~ 697 K SNPs were analyzed using GEMMA and mvBIMBAM for multivariate association. Functional annotation clustering analysis with DAVID was performed on candidate genes. A large 8 Mb region on chromosome 5 that included ACSF3, as well as 32 additional SNPs, met genome-wide significance at P < 6.07E-7 or Log10(BF) = 3.0 for GEMMA and mvBIMBAM, respectively. On chromosome 23, SNPs were located within or near PCCB and MSL2. On chromosome 30, a SNP was located within IGDCC3. SNPs associated with coat color were also located within ADAM9, FAM109B, SULT1C4, RTDR1, BCR, and RGS7. DZIP1L was associated with ACL rupture. Several significant SNPs on chromosomes 2, 3, 7, 24, and 26 were located within uncharacterized regions or long non-coding RNA sequences. This study validates associations with the previous ACL rupture candidate genes ACSF3 and DZIP1L and identifies novel candidate genes. These variants could act as targets for treatment or as factors in disease prediction modeling. The study highlighted the importance of regulatory SNPs in the disease, as several significant SNPs were located within non-coding regions.

Gut Bacterial Community Determines the Therapeutic Effect of Ginsenoside on Canine Inflammatory Bowel Disease by Modulating the Colonic Mucosal Barrier

Inflammatory bowel disease (IBD) comprises systemic inflammatory conditions primarily affecting the gastrointestinal tract, including Crohn's disease and ulcerative colitis. This research aims to analyze the clinical symptoms and pathogenesis of a Dextran sodium sulfate (DSS)-induced canine IBD model and evaluate the restorative effect of ginsenoside from a pathogenesis perspective. We established the DSS-induced canine IBD model and studied the pathological mechanisms. Additionally, we examined the therapeutic effect of ginsenosides by assessing the Canine Inflammatory Bowel Disease Activity Index (CIBDAI), C-reactive protein (CRP) levels, colonic tissue morphology, protein expression, and mucosal bacterial community analysis. Our findings revealed a total ginsenoside content of 22.7% in the ginsenoside extract. Animal experiments demonstrated that dogs with IBD exhibited decreased mental state, significantly increased CIBDAI and CRP levels, disrupted colonic epithelial tissue structure, decreased expression of mucin, tight junctions, and adherens junctions, as well as reduced diversity of the colonic mucosal bacterial community. Furthermore, correlation analysis highlighted a total of 38 bacterial strains correlated with physiological indices. Significantly, ginsenoside treatment could improve these symptoms and reverse the relative abundance of some bacterial communities. In conclusion, alterations in the properties of the colonic mucus layer or the reduction in MUC2, its core component, in dogs with IBD can lead to bacterial penetration of the mucus layer and subsequent contact with intestinal epithelial cells, resulting in inflammation. Remarkably, ginsenoside intervention showcased the capacity to positively influence the relative abundance of bacteria and impact the colonic mucus layer properties, thereby offering promising prospects for IBD management and recovery.

Mechanical osteoarthritis of the hip in a one medicine concept: a narrative review

Human and veterinary medicine have historically presented many medical areas of potential synergy and convergence. Mechanical osteoarthritis (MOA) is characterized by a gradual complex imbalance between cartilage production, loss, and derangement. Any joint instability that results in an abnormal overload of the joint surface can trigger MOA. As MOA has a prevailing mechanical aetiology, treatment effectiveness can only be accomplished if altered joint mechanics and mechanosensitive pathways are normalized and restored. Otherwise, the inflammatory cascade of osteoarthritis will be initiated, and the changes may become irreversible. The management of the disease using non-steroidal anti-inflammatory drugs, analgesics, physical therapy, diet changes, or nutraceuticals is conservative and less effective. MOA is a determinant factor for the development of hip dysplasia in both humans and dogs. Hip dysplasia is a hereditary disease with a high incidence and, therefore, of great clinical importance due to the associated discomfort and significant functional limitations. Furthermore, on account of analogous human and canine hip dysplasia disease and under the One Medicine concept, unifying veterinary and human research could improve the well-being and health of both species, increasing the acknowledgement of shared diseases. Great success has been accomplished in humans regarding preventive conservative management of hip dysplasia and following One Medicine concept, similar measures would benefit dogs. Moreover, animal models have long been used to better understand the different diseases' mechanisms. Current research in animal models was addressed and the role of rabbit models in pathophysiologic studies and of the dog as a spontaneous animal model were highlighted, denoting the inexistence of rabbit functional models to investigate therapeutic approaches in hip MOA.

A comprehensive analysis of gut and skin microbiota in canine atopic dermatitis in Shiba Inu dogs

BACKGROUND

Like its human counterpart, canine atopic dermatitis (cAD) is a chronic relapsing condition; thus, most cAD-affected dogs will require lifelong treatment to maintain an acceptable quality of life. A potential intervention is modulation of the composition of gut microbiota, and in fact, probiotic treatment has been proposed and tried in human atopic dermatitis (AD) patients. Since dogs are currently receiving intensive medical care, this will be the same option for dogs, while evidence of gut dysbiosis in cAD is still missing, although skin microbial profiling in cAD has been conducted in several studies. Therefore, we conducted a comprehensive analysis of both gut and skin microbiota in cAD in one specific cAD-predisposed breed, Shiba Inu. Additionally, we evaluated the impact of commonly used medical management on cAD (Janus kinase; JAK inhibitor, oclacitinib) on the gut and skin microbiota. Furthermore, we genotyped the Shiba Inu dogs according to the mitochondrial DNA haplogroup and assessed its association with the composition of the gut microbiota.

RESULTS

Staphylococcus was the most predominant bacterial genus observed in the skin; Escherichia/Shigella and Clostridium sensu stricto were highly abundant in the gut of cAD-affected dogs. In the gut microbiota, Fusobacteria and Megamonas were highly abundant in healthy dogs but significantly reduced in cAD-affected dogs. The abundance of these bacterial taxa was positively correlated with the effect of the treatment and state of the disease. Oclacitinib treatment on cAD-affected dogs shifted the composition of microbiota towards that in healthy dogs, and the latter brought it much closer to healthy microbiota, particularly in the gut. Additionally, even within the same dog breed, the mtDNA haplogroup varied, and there was an association between the mtDNA haplogroup and microbial composition in the gut and skin.

CONCLUSIONS

Dysbiosis of both the skin and the gut was observed in cAD in Shiba Inu dogs. Our findings provide a basis for the potential treatment of cAD by manipulating the gut microbiota as well as the skin microbiota. Video Abstract.

Lymphoma in Border Collies: Genome-Wide Association and Pedigree Analysis

There has been considerable interest in studying cancer in dogs and its potential as a model system for humans. One area of research has been the search for genetic risk variants in canine lymphoma, which is amongst the most common canine cancers. Previous studies have focused on a limited number of breeds, but none have included Border Collies. The aims of this study were to identify relationships between Border Collie lymphoma cases through an extensive pedigree investigation and to utilise relationship information to conduct genome-wide association study (GWAS) analyses to identify risk regions associated with lymphoma. The expanded pedigree analysis included 83,000 Border Collies, with 71 identified lymphoma cases. The analysis identified affected close relatives, and a common ancestor was identified for 54 cases. For the genomic study, a GWAS was designed to incorporate lymphoma cases, putative "carriers", and controls. A case-control GWAS was also conducted as a comparison. Both analyses showed significant SNPs in regions on chromosomes 18 and 27. Putative top candidate genes from these regions included DLA-79, WNT10B, LMBR1L, KMT2D, and CCNT1.

Red Ginseng Dietary Fiber Shows Prebiotic Potential by Modulating Gut Microbiota in Dogs

Red ginseng, widely used in traditional medicine for various conditions, imparts health benefits mainly by modulating the gut microbiota in humans. Given the similarities in gut microbiota between humans and dogs, red ginseng-derived dietary fiber may have prebiotic potential in dogs; however, its effects on the gut microbiota in dogs remain elusive. This double-blinded, longitudinal study investigated the impact of red ginseng dietary fiber on the gut microbiota and host response in dogs. A total of 40 healthy household dogs were randomly assigned to low-dose (n = 12), high-dose (n = 16), or control (n = 12) groups and fed a normal diet supplemented with red ginseng dietary fiber (3 g/5 kg body weight per day, 8 g/5 kg per day, or no supplement, respectively) for 8 weeks. The gut microbiota of the dogs was analyzed at 4 weeks and 8 weeks using 16S rRNA gene sequencing of fecal samples. Alpha diversity was significantly increased at 8 and 4 weeks in the low-dose and high-dose groups, respectively. Moreover, biomarker analysis showed that short-chain fatty acid producers such as Sarcina and Proteiniclasticum were significantly enriched, while potential pathogens such as Helicobacter were significantly decreased, indicating the increased gut health and pathogen resistance by red ginseng dietary fiber. Microbial network analysis showed that the complexity of microbial interactions was increased by both doses, indicating the increased stability of the gut microbiota. These findings suggest that red ginseng-derived dietary fiber could be used as a prebiotic to modulate gut microbiota and improve gut health in dogs. IMPORTANCE The canine gut microbiota is an attractive model for translational studies, as it responds to dietary interventions similarly to those in humans. Investigating the gut microbiota of household dogs that share the environment with humans can produce highly generalizable and reproducible results owing to their representativeness of the general canine population. This double-blind and longitudinal study investigated the impact of dietary fiber derived from red ginseng on the gut microbiota of household dogs. Red ginseng dietary fiber altered the canine gut microbiota by increasing diversity, enriching short-chain fatty acid-producing microbes, decreasing potential pathogens, and increasing the complexity of microbial interactions. These findings indicate that red ginseng-derived dietary fiber may promote canine gut health by modulating gut microbiota, suggesting the possibility of its use as a potential prebiotic.

A case of canine intestinal malakoplakia

Malakoplakia is a rare chronic granulomatous disease usually affecting the urinary bladder and other locations. In humans, the gastrointestinal tract is the second most common location but there are no reports of intestinal malakoplakia in animals. A 10-month-old female French Bulldog was presented with chronic haemorrhagic diarrhoea and anorexia with normochromic-normocytic anaemia and hypoalbuminaemia. Grossly, there was mucosal thickening and ulceration of the caecum, colon and rectum. Microscopically, transmural sheets of foamy macrophages were seen in these tissues. Macrophages were periodic acid-Schiff, vimentin and ionized calcium-binding adaptor molecule 1 positive and contained von Kossa- and Prussian blue-positive Michaelis-Gutmann bodies. Giemsa staining revealed rod-shaped bacterial colonies and fluorescence in-situ hybridization demonstrated Escherichia coli within macrophages. This is the first reported case of intestinal malakoplakia in domestic animals. Pathological features of intestinal malakoplakia share many similarities with ulcerative histiocytic colitis in dogs but it is unclear if they are different forms of the same pathological process or distinct entities.

Genomic analysis and prediction of genomic values for distichiasis in Staffordshire bull terriers

BACKGROUND

Distichiasis is a condition characterized by aberrant hairs along the eyelid margins. The symptoms are usually mild but can lead to ulcerations and lesions of the cornea in severe cases. It is the most frequently noted ocular disorder in Norwegian Staffordshire bull terriers (SBT), with a prevalence above 18% in the adult population. A complex inheritance is assumed, but there is sparse knowledge about the genetic background of distichiasis in dogs. We have performed a genome-wide association study of distichiasis in SBT and used genomic data in an attempt to predict genomic values for the disorder.

RESULTS

We identified four genetic regions on CFA1, CFA18, CFA32 and CFA34 using a mixed linear model association analysis and a Bayesian mixed model analysis. Genomic values were predicted using GBLUP and a Bayesian approach, BayesR. The genomic prediction showed that the 1/4 of dogs with predicted values most likely to acquire distichiasis had a 3.9 -4.0 times higher risk of developing distichiasis compared to the quarter (1/4) of dogs least likely to acquire the disease. There was no significant difference between the two methods used.

CONCLUSION

Four genomic regions associated with distichiasis were discovered in the association analysis, suggesting that distichiasis in SBT is a complex trait involving numerous loci. The four associated regions need to be confirmed in an independent sample. We also used all 95 K SNPs for genomic prediction and showed that genomic prediction can be a helpful tool in selective breeding schemes at breed level aiming at reducing the prevalence of distichiasis in SBTs in the future, even if the predictive value of single dogs may be low.

Integrative mapping of the dog epigenome: Reference annotation for comparative intertissue and cross-species studies

Dogs have become a valuable model in exploring multifaceted diseases and biology relevant to human health. Despite large-scale dog genome projects producing high-quality draft references, a comprehensive annotation of functional elements is still lacking. We addressed this through integrative next-generation sequencing of transcriptomes paired with five histone marks and DNA methylome profiling across 11 tissue types, deciphering the dog's epigenetic code by defining distinct chromatin states, super-enhancer, and methylome landscapes, and thus showed that these regions are associated with a wide range of biological functions and cell/tissue identity. In addition, we confirmed that the phenotype-associated variants are enriched in tissue-specific regulatory regions and, therefore, the tissue of origin of the variants can be traced. Ultimately, we delineated conserved and dynamic epigenomic changes at the tissue- and species-specific resolutions. Our study provides an epigenomic blueprint of the dog that can be used for comparative biology and medical research.

Systematically identifying genetic signatures including novel SNP-clusters, nonsense variants, frame-shift INDELs, and long STR expansions that potentially link to unknown phenotypes existing in dog breeds

BACKGROUND

In light of previous studies that profiled breed-specific traits or used genome-wide association studies to refine loci associated with characteristic morphological features in dogs, the field has gained tremendous genetic insights for known dog traits observed among breeds. Here we aim to address the question from a reserve perspective: whether there are breed-specific genotypes that may underlie currently unknown phenotypes. This study provides a complete set of breed-specific genetic signatures (BSGS). Several novel BSGS with significant protein-altering effects were highlighted and validated.

RESULTS

Using the next generation whole-genome sequencing technology coupled with unsupervised machine learning for pattern recognitions, we constructed and analyzed a high-resolution sequence map for 76 breeds of 412 dogs. Genomic structures including novel single nucleotide polymorphisms (SNPs), SNP clusters, insertions, deletions (INDELs) and short tandem repeats (STRs) were uncovered mutually exclusively among breeds. We also partially validated some novel nonsense variants by Sanger sequencing with additional dogs. Four novel nonsense BSGS were found in the Bernese Mountain Dog, Samoyed, Bull Terrier, and Basset Hound, respectively. Four INDELs resulting in either frame-shift or codon disruptions were found in the Norwich Terrier, Airedale Terrier, Chow Chow and Bernese Mountain Dog, respectively. A total of 15 genomic regions containing three types of BSGS (SNP-clusters, INDELs and STRs) were identified in the Akita, Alaskan Malamute, Chow Chow, Field Spaniel, Keeshond, Shetland Sheepdog and Sussex Spaniel, in which Keeshond and Sussex Spaniel each carried one amino-acid changing BSGS in such regions.

CONCLUSION

Given the strong relationship between human and dog breed-specific traits, this study might be of considerable interest to researchers and all. Novel genetic signatures that can differentiate dog breeds were uncovered. Several functional genetic signatures might indicate potentially breed-specific unknown phenotypic traits or disease predispositions. These results open the door for further investigations. Importantly, the computational tools we developed can be applied to any dog breeds as well as other species. This study will stimulate new thinking, as the results of breed-specific genetic signatures may offer an overarching relevance of the animal models to human health and disease.

Psychological Stress Is Associated with Increased Cancer Risk in Dogs

Although there is evidence that psychological stress may be associated with increased cancer risk, the effect of stress on cancer risk is difficult to study, both in humans, due to socioeconomic factors, and in animal models, due to questionable biological relevance. Here, we test whether heritable canine temperament that increases psychological stress is associated with cancer risk. The study data are breed-specific averages of incidences of multiple cancer types and of temperament classes. The latter are derived from a latent class analysis of behavioral questionnaires completed by owners (C-BARQ). We thus classified the dogs according to whether they are calm vs. reactive within and across breeds. Using meta-analysis approaches, we modeled the risk of multiple cancer types in calm vs. reactive dogs. We adjusted for breed averages of body mass and lifespan, which are common confounders that impact cancer. Our study confirms that body size has a significant effect of on risk of multiple types of cancers in dogs and shows for the first time that temperament also has a moderate effect. These findings suggest dog models of heritable psychological stress are suitable for molecular epidemiological and translational studies on its effects on cancer risk.

Common Orthopedic Traits and Screening for Breeding Programs

Orthopedic diseases are complex traits, meaning genetics and environmental factors affect risk, making identification of genetic associations difficult. In the United States, hip and elbow scores, patellar luxation scores, Legg-Calvé-Perthes disease, and shoulder osteochondrosis affectedness are available in the Orthopedic Foundation for Animals registry. Distraction indices and extended, ventrodorsal hip conformation scores are recorded by PennHIP. Application of estimated breeding values for hip and elbow dysplasia in breeder selection reduces the severity and prevalence of these traits. Genomic prediction and whole-genome sequence technologies and methods should improve knowledge of genetics underlying orthopedic diseases, leading to improved canine orthopedic genetic quality.

​Comparative genomics of Balto, a famous historic dog, captures lost diversity of 1920s sled dogs

We reconstruct the phenotype of Balto, the heroic sled dog renowned for transporting diphtheria antitoxin to Nome, Alaska, in 1925, using evolutionary constraint estimates from the Zoonomia alignment of 240 mammals and 682 genomes from dogs and wolves of the 21st century. Balto shares just part of his diverse ancestry with the eponymous Siberian husky breed. Balto's genotype predicts a combination of coat features atypical for modern sled dog breeds, and a slightly smaller stature. He had enhanced starch digestion compared with Greenland sled dogs and a compendium of derived homozygous coding variants at constrained positions in genes connected to bone and skin development. We propose that Balto's population of origin, which was less inbred and genetically healthier than that of modern breeds, was adapted to the extreme environment of 1920s Alaska.

Dynamic chromatin architectures provide insights into the genetics of cattle myogenesis

BACKGROUND

Sharply increased beef consumption is propelling the genetic improvement projects of beef cattle in China. Three-dimensional genome structure is confirmed to be an important layer of transcription regulation. Although genome-wide interaction data of several livestock species have already been produced, the genome structure states and its regulatory rules in cattle muscle are still limited.

RESULTS

Here we present the first 3D genome data in Longissimus dorsi muscle of fetal and adult cattle (Bos taurus). We showed that compartments, topologically associating domains (TADs), and loop undergo re-organization and the structure dynamics were consistent with transcriptomic divergence during muscle development. Furthermore, we annotated cis-regulatory elements in cattle genome during myogenesis and demonstrated the enrichments of promoter and enhancer in selection sweeps. We further validated the regulatory function of one HMGA2 intronic enhancer near a strong sweep region on primary bovine myoblast proliferation.

CONCLUSIONS

Our data provide key insights of the regulatory function of high order chromatin structure and cattle myogenic biology, which will benefit the progress of genetic improvement of beef cattle.

Identification of Genomic Signatures in Bullmastiff Dogs Using Composite Selection Signals Analysis of 23 Purebred Clades

Dog breeds represent canine sub-populations with distinctive phenotypic features and limited genetic diversity. We have established a resource to study breed-specific genetic diversity. Utilising genetic resources within our laboratory biobank, public domain genotype data and the phylogenetic framework of 23 breed clades, the primary objective for this study was to identify genomic regions that differentiate the Bullmastiff breed. Through application of a composite index analysis (CSS), genomic signatures were identified in Bullmastiffs when compared to the formative breeds, Mastiffs and Bulldogs, and to 22 other breed groups. Significant regions were identified on 15 chromosomes, with the most differentiated regions found on CFA1, CFA9, and CFA18. These regions may reflect genetic drift following establishment of the breed or the effects of selective breeding during development of the modern Bullmastiff. This was supported by analysis of genes from the identified genomic regions, including 458 genes from the multi-clade analysis, which revealed enriched pathways that may be related to characteristic traits and distinct morphology of the breed. The study demonstrates the utility of the CSS method in breed-specific genome analysis and advances our understanding of genetic diversity in Bullmastiff dogs.

Banking on a new understanding: translational opportunities from veterinary biobanks

Current advances in geroscience are due in part to the discovery of biomarkers with high predictive ability in short-lived laboratory animals such as flies and mice. These model species, however, do not always adequately reflect human physiology and disease, highlighting the need for a more comprehensive and relevant model of human aging. Domestic dogs offer a solution to this obstacle, as they share many aspects not only of the physiological and pathological trajectories of their human counterpart, but also of their environment. Furthermore, they age at a considerably faster rate. Studying aging in the companion dog provides an opportunity to better understand the biological and environmental determinants of healthy lifespan in our pets, and to translate those findings to human aging. Biobanking, the systematic collection, processing, storage, and distribution of biological material and associated data has contributed to basic, clinical, and translational research by streamlining the management of high-quality biospecimens for biomarker discovery and validation. In this review, we discuss how veterinary biobanks can support research on aging, particularly when integrated into large-scale longitudinal studies. As an example of this concept, we introduce the Dog Aging Project Biobank.

Translational veterinary epilepsy: A win-win situation for human and veterinary neurology

Epilepsy is a challenging multifactorial disorder with a complex genetic background. Our current understanding of the pathophysiology and treatment of epilepsy has substantially increased due to animal model studies, including canine studies, but additional basic and clinical research is required. Drug-resistant epilepsy is an important problem in both dogs and humans, since seizure freedom is not achieved with the available antiseizure medications. The evaluation and exploration of pharmacological and particularly non-pharmacological therapeutic options need to remain a priority in epilepsy research. Combined efforts and sharing knowledge and expertise between human medical and veterinary neurologists are important for improving the treatment outcomes or even curing epilepsy in dogs. Such interactions could offer an exciting approach to translate the knowledge gained from people and rodents to dogs and vice versa. In this article, a panel of experts discusses the similarities and knowledge gaps in human and animal epileptology, with the aim of establishing a common framework and the basis for future translational epilepsy research.

Identification of Genetic Risk Factors for Monogenic and Complex Canine Diseases

Advances in DNA sequencing and other technologies have greatly facilitated the identification of genetic risk factors for inherited diseases in dogs. We review recent technological developments based on selected examples from canine disease genetics. The identification of disease-causing variants in dogs with monogenic diseases may become a widely employed diagnostic approach in clinical veterinary medicine in the not-too-distant future. Diseases with complex modes of inheritance continue to pose challenges to researchers but have also become much more tangible than in the past. In addition to strategies for identifying genetic risk factors, we provide some thoughts on the interpretation of sequence variants that are largely inspired by developments in human clinical genetics.

Long-term Small Population Size, Deleterious Variation, and Altitude Adaptation in the Ethiopian Wolf, a Severely Endangered Canid

Ethiopian wolves, a canid species endemic to the Ethiopian Highlands, have been steadily declining in numbers for decades. Currently, out of 35 extant species, it is now one of the world's most endangered canids. Most conservation efforts have focused on preventing disease, monitoring movements and behavior, and assessing the geographic ranges of sub-populations. Here, we add an essential layer by determining the Ethiopian wolf's demographic and evolutionary history using high-coverage (∼40×) whole-genome sequencing from 10 Ethiopian wolves from the Bale Mountains. We observe exceptionally low diversity and enrichment of weakly deleterious variants in the Ethiopian wolves in comparison with two North American gray wolf populations and four dog breeds. These patterns are consequences of long-term small population size, rather than recent inbreeding. We infer the demographic history of the Ethiopian wolf and find it to be concordant with historic records and previous genetic analyses, suggesting Ethiopian wolves experienced a series of both ancient and recent bottlenecks, resulting in a census population size of fewer than 500 individuals and an estimated effective population size of approximately 100 individuals. Additionally, long-term small population size may have limited the accumulation of strongly deleterious recessive mutations. Finally, as the Ethiopian wolves have inhabited high-altitude areas for thousands of years, we searched for evidence of high-altitude adaptation, finding evidence of positive selection at a transcription factor in a hypoxia-response pathway [CREB-binding protein (CREBBP)]. Our findings are pertinent to continuing conservation efforts and understanding how demography influences the persistence of deleterious variation in small populations.

Domestic dog lineages reveal genetic drivers of behavioral diversification

Selective breeding of domestic dogs has generated diverse breeds often optimized for performing specialized tasks. Despite the heritability of breed-typical behavioral traits, identification of causal loci has proven challenging due to the complexity of canine population structure. We overcome longstanding difficulties in identifying genetic drivers of canine behavior by developing a framework for understanding relationships between breeds and the behaviors that define them, utilizing genetic data for over 4,000 domestic, semi-feral, and wild canids and behavioral survey data for over 46,000 dogs. We identify ten major canine genetic lineages and their behavioral correlates and show that breed diversification is predominantly driven by non-coding regulatory variation. We determine that lineage-associated genes converge in neurodevelopmental co-expression networks, identifying a sheepdog-associated enrichment for interrelated axon guidance functions. This work presents a scaffold for canine diversification that positions the domestic dog as an unparalleled system for revealing the genetic origins of behavioral diversity.

Across-breed genetic investigation of canine hip dysplasia, elbow dysplasia, and anterior cruciate ligament rupture using whole-genome sequencing

Here, we report the use of genome-wide association study (GWAS) for the analysis of canine whole-genome sequencing (WGS) repository data using breed phenotypes. Single-nucleotide polymorphisms (SNPs) were called from WGS data from 648 dogs that included 119 breeds from the Dog10K Genomes Project. Next, we assigned breed phenotypes for hip dysplasia (Orthopedic Foundation for Animals (OFA) HD, n = 230 dogs from 27 breeds; hospital HD, n = 279 dogs from 38 breeds), elbow dysplasia (ED, n = 230 dogs from 27 breeds), and anterior cruciate ligament rupture (ACL rupture, n = 279 dogs from 38 breeds), the three most important canine spontaneous complex orthopedic diseases. Substantial morbidity is common with these diseases. Previous within- and between-breed GWAS for HD, ED, and ACL rupture using array SNPs have identified disease-associated loci. Individual disease phenotypes are lacking in repository data. There is a critical knowledge gap regarding the optimal approach to undertake categorical GWAS without individual phenotypes. We considered four GWAS approaches: a classical linear mixed model, a haplotype-based model, a binary case-control model, and a weighted least squares model using SNP average allelic frequency. We found that categorical GWAS was able to validate HD candidate loci. Additionally, we discovered novel candidate loci and genes for all three diseases, including FBX025, IL1A, IL1B, COL27A1, SPRED2 (HD), UGDH, FAF1 (ED), TGIF2 (ED & ACL rupture), and IL22, IL26, CSMD1, LDHA, and TNS1 (ACL rupture). Therefore, categorical GWAS of ancestral dog populations may contribute to the understanding of any disease for which breed epidemiological risk data are available, including diseases for which GWAS has not been performed and candidate loci remain elusive.

Races canines et comportement

It is commonly believed that a dog’s breed predicts its behaviour. A 2019 paper seemed to confirm this, but had some serious limitations. An extensive new study shows that, in fact, breed is not a good predictor of behaviour at the individual level, and has only a weak influence on average characteristics.

Detection of invasive Escherichia coli in dogs with granulomatous colitis using immunohistochemistry

Granulomatous colitis in dogs can be associated with infection of the colonic mucosa by invasive strains of Escherichia coli. To date, fluorescence in situ hybridization (FISH) is the gold-standard method to assess intramucosal and intracellular bacterial invasion. However, FISH requires expensive fluorescence microscopy equipment and is therefore not widely available. We investigated the use of immunohistochemistry (IHC) as an alternative method to detect invasive E. coli in dogs with granulomatous colitis. Archived paraffin-embedded blocks were selected from 26 dogs with colitis, in which FISH had been performed by an outside laboratory. Using a polyclonal antibody, IHC for E. coli was performed on sections cut from the same blocks, and the presence of invasive E. coli was recorded. All 11 specimens in which FISH had detected E. coli were also positive on IHC, with strong immunolabeling in the cytoplasm of macrophages and extracellularly in the lamina propria; all 15 specimens that were negative for invasive bacteria on FISH were also negative on IHC. We found that IHC is a sensitive technique for the detection of invasive E. coli in dogs with granulomatous colitis.

A genome-wide association study to investigate genetic loci associated with primary glaucoma in American Cocker Spaniels

OBJECTIVE

To identify genetic associations with primary glaucoma (PG) in American Cocker Spaniels using a genome-wide association study (GWAS).

ANIMALS

A nationwide ambidirectional case-control cohort study was performed in American Cocker Spaniels that had an ophthalmic examination performed by a veterinarian. Ninety-four dogs with PG (cases) and 111 dogs without glaucoma (controls) met phenotypic criteria and had a blood sample collected after receiving informed owner consent.

PROCEDURES

Genomic DNA was extracted from whole blood samples and genotyped (CanineHD BeadChip, Illumina Inc). A case-control GWAS using a linear mixed model was performed, and 3 significance thresholds were calculated (1) using a Bonferroni correction on all single nucleotide polymorphisms (SNPs) included in the GWAS, (2) using a Bonferroni correction on only the unlinked SNPs from a pruned data set, and (3) using 10,000 random phenotype permutations.

RESULTS

Following genotype data quality control, 89 cases and 93 controls were included in the GWAS. We identified an association on canine chromosome (CFA10); however, it did not reach statistical significance. Potential candidate genes within the surrounding linkage disequilibrium interval include coiled-coil domain containing 85A (CCDC85A) and extracellular growth factor containing fibulin extracellular matrix protein 1 (EFEMP1).

CLINICAL RELEVANCE

Primary glaucoma in the American Cocker Spaniel is a complex heterogeneous disease that may be influenced by a locus on CFA10. The candidate genes CCDC85A and EFEMP1 within the identified linkage disequilibrium interval have been shown to be involved in human open-angle glaucoma.

The genomic landscape of canine osteosarcoma cell lines reveals conserved structural complexity and pathway alterations

The characterization of immortalized canine osteosarcoma (OS) cell lines used for research has historically been based on phenotypic features such as cellular morphology and expression of bone specific markers. With the increasing use of these cell lines to investigate novel therapeutic approaches prior to in vivo translation, a much more detailed understanding regarding the genomic landscape of these lines is required to ensure accurate interpretation of findings. Here we report the first whole genome characterization of eight canine OS cell lines, including single nucleotide variants, copy number variants and other structural variants. Many alterations previously characterized in primary canine OS tissue were observed in these cell lines, including TP53 mutations, MYC copy number gains, loss of CDKN2A, PTEN, DLG2, MAGI2, and RB1 and structural variants involving SETD2, DLG2 and DMD. These data provide a new framework for understanding how best to incorporate in vitro findings generated using these cell lines into the design of future clinical studies involving dogs with spontaneous OS.

Predicting Dog Phenotypes from Genotypes

We analyze dog genotypes (i.e., positions of dog DNA sequences that often vary between different dogs) in order to predict the corresponding phenotypes (i.e., unique observed characteristics). More specifically, given chromosome data from a dog, we aim to predict the breed, height, and weight. We explore a variety of linear and non-linear classification and regression techniques to accomplish these three tasks. We also investigate the use of a neural network (both in linear and non-linear modes) for breed classification and compare the performance to traditional statistical methods. We show that linear methods generally outperform or match the performance of non-linear methods for breed classification. However, we show that the reverse is true for height and weight regression. Finally, we evaluate the results of all of these methods based on the number of input features used in the analysis. We conduct experiments using different fractions of the full genomic sequences, resulting in input sequences ranging from 20 SNPs to ∼200k SNPs. In doing so, we explore the impact of using a very limited number of SNPs for prediction. Our experiments demonstrate that these phenotypes in dogs can be predicted with as few as 0.5% of randomly selected SNPs (i.e., 992 SNPs) and that dog breeds can be classified with 50% balanced accuracy with as few as 0.02% SNPs (i.e., 40 SNPs).

Reconstructing the history of founder events using genome-wide patterns of allele sharing across individuals

Founder events play a critical role in shaping genetic diversity, fitness and disease risk in a population. Yet our understanding of the prevalence and distribution of founder events in humans and other species remains incomplete, as most existing methods require large sample sizes or phased genomes. Thus, we developed ASCEND that measures the correlation in allele sharing between pairs of individuals across the genome to infer the age and strength of founder events. We show that ASCEND can reliably estimate the parameters of founder events under a range of demographic scenarios. We then apply ASCEND to two species with contrasting evolutionary histories: ~460 worldwide human populations and ~40 modern dog breeds. In humans, we find that over half of the analyzed populations have evidence for recent founder events, associated with geographic isolation, modes of sustenance, or cultural practices such as endogamy. Notably, island populations have lower population sizes than continental groups and most hunter-gatherer, nomadic and indigenous groups have evidence of recent founder events. Many present-day groups––including Native Americans, Oceanians and South Asians––have experienced more extreme founder events than Ashkenazi Jews who have high rates of recessive diseases due their known history of founder events. Using ancient genomes, we show that the strength of founder events differs markedly across geographic regions and time––with three major founder events related to the peopling of Americas and a trend in decreasing strength of founder events in Europe following the Neolithic transition and steppe migrations. In dogs, we estimate extreme founder events in most breeds that occurred in the last 25 generations, concordant with the establishment of many dog breeds during the Victorian times. Our analysis highlights a widespread history of founder events in humans and dogs and elucidates some of the demographic and cultural practices related to these events.

A founder event occurs when small numbers of ancestral individuals give rise to a large fraction of the population. Founder events reduce genetic variation and increase the risk of recessive diseases. Despite their importance in evolutionary and disease studies, we still only have a limited comprehension of their prevalence and properties in humans and other species, as most existing methods require large sample sizes or phased genomes. Here, we present a flexible method, ASCEND, to infer the timing and the strength of founder events that is suitable for sparse datasets with few samples or limited coverage. ASCEND provides reliable estimates across a wide range of demographic scenarios. By applying it to data from two species (humans and dogs), we document a widespread history of recent founder events in both species and provide insights about the demographic processes related to these events. Our analysis helps to identify groups with strong founder events that should be prioritized for future studies as they offer a unique opportunity for biological discovery and reducing disease burden through mapping of recessive disease-causing genes and pathways, as previously shown in studies of Ashkenazi Jews and Finns.

Dogs as a Natural Animal Model of Epilepsy

Epilepsy is a common neurological disease in both humans and domestic dogs, making dogs an ideal translational model of epilepsy. In both species, epilepsy is a complex brain disease characterized by an enduring predisposition to generate spontaneous recurrent epileptic seizures. Furthermore, as in humans, status epilepticus is one of the more common neurological emergencies in dogs with epilepsy. In both species, epilepsy is not a single disease but a group of disorders characterized by a broad array of clinical signs, age of onset, and underlying causes. Brain imaging suggests that the limbic system, including the hippocampus and cingulate gyrus, is often affected in canine epilepsy, which could explain the high incidence of comorbid behavioral problems such as anxiety and cognitive alterations. Resistance to antiseizure medications is a significant problem in both canine and human epilepsy, so dogs can be used to study mechanisms of drug resistance and develop novel therapeutic strategies to benefit both species. Importantly, dogs are large enough to accommodate intracranial EEG and responsive neurostimulation devices designed for humans. Studies in epileptic dogs with such devices have reported ictal and interictal events that are remarkably similar to those occurring in human epilepsy. Continuous (24/7) EEG recordings in a select group of epileptic dogs for >1 year have provided a rich dataset of unprecedented length for studying seizure periodicities and developing new methods for seizure forecasting. The data presented in this review substantiate that canine epilepsy is an excellent translational model for several facets of epilepsy research. Furthermore, several techniques of inducing seizures in laboratory dogs are discussed as related to therapeutic advances. Importantly, the development of vagus nerve stimulation as a novel therapy for drug-resistant epilepsy in people was based on a series of studies in dogs with induced seizures. Dogs with naturally occurring or induced seizures provide excellent large-animal models to bridge the translational gap between rodents and humans in the development of novel therapies. Furthermore, because the dog is not only a preclinical species for human medicine but also a potential patient and pet, research on this species serves both veterinary and human medicine.

Complex Feline Disease Mapping Using a Dense Genotyping Array

The current feline genotyping array of 63 k single nucleotide polymorphisms has proven its utility for mapping within breeds, and its use has led to the identification of variants associated with Mendelian traits in purebred cats. However, compared to single gene disorders, association studies of complex diseases, especially with the inclusion of random bred cats with relatively low linkage disequilibrium, require a denser genotyping array and an increased sample size to provide statistically significant associations. Here, we undertook a multi-breed study of 1,122 cats, most of which were admitted and phenotyped for nine common complex feline diseases at the Cornell University Hospital for Animals. Using a proprietary 340 k single nucleotide polymorphism mapping array, we identified significant genome-wide associations with hyperthyroidism, diabetes mellitus, and eosinophilic keratoconjunctivitis. These results provide genomic locations for variant discovery and candidate gene screening for these important complex feline diseases, which are relevant not only to feline health, but also to the development of disease models for comparative studies.

DNA methylation clocks for dogs and humans

DNA methylation profiles have been used to develop biomarkers of aging known as epigenetic clocks, which predict chronological age with remarkable accuracy and show promise for inferring health status as an indicator of biological age. Epigenetic clocks were first built to monitor human aging, but their underlying principles appear to be evolutionarily conserved, as they have now been successfully developed for many mammalian species. Here, we describe reliable and highly accurate epigenetic clocks shown to apply to 93 domestic dog breeds. The methylation profiles were generated using the mammalian methylation array, which utilizes DNA sequences that are conserved across all mammalian species. Canine epigenetic clocks were constructed to estimate age and also average time to death. We also present two highly accurate human–dog dual species epigenetic clocks (R = 0.97), which may facilitate the ready translation from canine to human use (or vice versa) of antiaging treatments being developed for longevity and preventive medicine. Finally, epigenome-wide association studies here reveal individual methylation sites that may underlie the inverse relationship between breed weight and lifespan. Overall, we describe robust biomarkers to measure aging and, potentially, health status in canines.

Genetic Origins of the Two Canis lupus familiaris (Dog) Freight Dog Populations

Despite periodic drops in popularity, Arctic sled dogs continue to play a vital role in northern societies, providing both freight transit and recreational race activities. In this study, we selected the Mackenzie River Husky, a freight dog of complex history, and the Chinook, an American Kennel Club recognized freight dog breed whose heritage reportedly overlaps that of the MKRH, for detailed population analysis. We tested each to determine their component breeds and used admixture analysis to ascertain their population structure. We utilized haplotype analysis to identify genomic regions shared between each population and their founding breeds. Our data show that the Alaskan Malamutes and modern Greenland sled dog contributed to both populations, but there are also unexpected contributions from the German Shepherd dog and Collie. We used haplotype analysis to identify genomic regions nearing fixation in population type and identify provocative genes in each region. Finally, in response to recent reports regarding the importance of dietary lipid genes in Arctic dogs, we analyzed 8 such genes in a targeted analysis observing signatures of selection in both populations at the MLXIPL gene loci. These data highlight the genetic routes that breeds of similar function have taken toward their occupation as successful sled dogs.

Genetic architectures and selection signatures of body height in Chinese indigenous donkeys revealed by next-generation sequencing

Donkeys are widely distributed labour animals in the world. During the process of the domestication and artificial selection of domestic donkeys, body sizes show significant differences among different breeds of donkeys. Based on the genome resequencing data of 103 Chinese indigenous donkeys from 11 breeds (Biyang, Dezhou, Guangling, Hetian, Jiami, Kulun, Qingyang, Turfan, Tibetan, Xinjiang, and Yunnan), seven Spanish donkeys from two breeds (Zamorano~Leonés and Andalusian), and three wild donkeys, we investigated the population structures of Chinese domestic donkeys with different body sizes. We used F_ST and XP-EHH analyses to explore the selected regions related to body sizes. The results showed that Chinese indigenous donkeys have a closer relationship with African wild donkeys than with Asian wild donkeys. LCORL/NCAPG, FAM184B, TBX3, and IHH were identified as genes with strong signals in analysis of selection signature (F_ST and XP-EHH) in large and small donkeys. The seven identified variants can be served as candidate loci affecting the body size of Chinese donkeys. Five of seven loci were located in intron 9 of FAM184B and were in a haplotype block, and one of the identified variants (Chr03:112664848) located in the CDS region of the LCORL gene was found to cause stop-loss. These candidate genes and variants shed new light on the molecular basis of donkey body size and will facilitate the breeding activities of donkeys.

Ancestry-inclusive dog genomics challenges popular breed stereotypes

Behavioral genetics in dogs has focused on modern breeds, which are isolated subgroups with distinctive physical and, purportedly, behavioral characteristics. We interrogated breed stereotypes by surveying owners of 18,385 purebred and mixed-breed dogs and genotyping 2155 dogs. Most behavioral traits are heritable [heritability (h²) > 25%], and admixture patterns in mixed-breed dogs reveal breed propensities. Breed explains just 9% of behavioral variation in individuals. Genome-wide association analyses identify 11 loci that are significantly associated with behavior, and characteristic breed behaviors exhibit genetic complexity. Behavioral loci are not unusually differentiated in breeds, but breed propensities align, albeit weakly, with ancestral function. We propose that behaviors perceived as characteristic of modern breeds derive from thousands of years of polygenic adaptation that predates breed formation, with modern breeds distinguished primarily by aesthetic traits.

Genomic and Transcriptomic Characterization of Atypical Recurrent Flank Alopecia in the Cesky Fousek

Non-inflammatory alopecia is a frequent skin problem in dogs, causing damaged coat integrity and compromised appearance of affected individuals. In this study, we examined the Cesky Fousek breed, which displays atypical recurrent flank alopecia (aRFA) at a high frequency. This type of alopecia can be quite severe and is characterized by seasonal episodes of well demarcated alopecic areas without hyperpigmentation. The genetic component responsible for aRFA remains unknown. Thus, here we aimed to identify variants involved in aRFA using a combination of histological, genomic, and transcriptomic data. We showed that aRFA is histologically similar to recurrent flank alopecia, characterized by a lack of anagen hair follicles and the presence of severely shortened telogen or kenogen hair follicles. We performed a genome-wide association study (GWAS) using 216 dogs phenotyped for aRFA and identified associations on chromosomes 19, 8, 30, 36, and 21, highlighting 144 candidate genes, which suggests a polygenic basis for aRFA. By comparing the skin cell transcription pattern of six aRFA and five control dogs, we identified 236 strongly differentially expressed genes (DEGs). We showed that the GWAS genes associated with aRFA are often predicted to interact with DEGs, suggesting their joint contribution to the development of the disease. Together, these genes affect four major metabolic pathways connected to aRFA: collagen formation, muscle structure/contraction, lipid metabolism, and the immune system.

A bitesize introduction to canine hematologic malignancies

Hematologic malignancies are frequently diagnosed in dogs and result in a spectrum of clinical signs associated with specific disease types. The most frequently encountered hematologic tumors in dogs include lymphoma, lymphoid and myeloid leukemias, and mast cell, plasma cell, and histiocytic neoplasias. Coupled with the heterogeneous presentations of the different categories and subtypes of canine hematologic malignancies, outcomes for these tumors are also variable. Considering this, appropriate treatment options range from active surveillance to curative intent approaches harnessing surgical, chemotherapeutic, and radiation-based modalities. The underlying pathology of many of these diseases bears remarkable resemblance to that of the corresponding diagnosis made in human patients. We introduce some of the pathogenic drivers of canine hematologic cancers alongside their clinical presentations. An overview of standard-of-care therapies for each of these diseases is also provided. As comparative oncology gains recognition as a valuable setting in which to investigate the pathogenesis of neoplasia and provide powerful, clinically relevant, immunocompetent models for the evaluation of novel therapies, the number of clinicians and scientists participating in cancer research involving dogs is expected to increase. This review aims at providing an introductory overview of canine hematologic malignancies.

Poly(A) RNA sequencing reveals age-related differences in the prefrontal cortex of dogs

Dogs may possess a unique translational potential to investigate neural aging and dementia because they are prone to age-related cognitive decline, including an Alzheimer’s disease–like pathological condition. Yet very little is known about the molecular mechanisms underlying canine cognitive decline. The goal of the current study was to explore the transcriptomic differences between young and old dogs’ frontal cortex, which is a brain region often affected by various forms of age-related dementia in humans. RNA isolates from the frontal cortical brain area of 13 pet dogs, which represented 7 different breeds and crossbreds, were analyzed. The dogs were euthanized for medical reasons, and their bodies had been donated by their owners for scientific purposes. The poly(A) tail RNA subfraction of the total transcriptome was targeted in the sequencing analysis. Cluster analyses, differential gene expression analyses, and gene ontology analyses were carried out to assess which genes and genetic regulatory mechanisms were mostly affected by aging. Age was the most prominent factor in the clustering of the animals, indicating the presence of distinct gene expression patterns related to aging in a genetically variable population. A total of 3436 genes were found to be differentially expressed between the age groups, many of which were linked to neural function, immune system, and protein synthesis. These findings are in accordance with previous human brain aging RNA sequencing studies. Some genes were found to behave more similarly to humans than to rodents, further supporting the applicability of dogs in translational aging research.