The impact of identity by descent on fitness and disease in dogs

Genetic diversity, demographic history, and selective signatures of Silkie chicken

BACKGROUND

Silkie is a traditional Chinese chicken breed characterized by its unique combination of specialized morphological traits. While previous studies have focused on the genetic basis of these traits, the overall genomic characteristics of the Silkie breed remain largely unexplored. In this study, we employed whole genome resequencing data to examine the genetic diversity, selective signals and demographic history of the Silkie breed through comparative analyses with seven other Chinese indigenous breeds (IDGBs), a commercial breed, and the wild ancestor Red Jungle Fowl.

RESULTS

In total, 20.8 million high-quality single nucleotide polymorphisms and 86 large structural variations were obtained. We discovered that Silkie exhibits a relatively high level of inbreeding and is genetically distinct from other IDGBs. Furthermore, our analysis indicated that Silkie has experienced a stronger historical population bottleneck and has a smaller effective population size compared with other IDGBs. We identified 45 putatively selected genes that are enriched in the melanogenesis pathway, which probably is related to the feather color. Among these genes, LMBR1 and PDSS2 have been previously associated with the extra toe and the hookless feathers, respectively. Six of the selected genes (KITLG, GSK3B, SOBP, CTBP1, ELMO2, SNRPN) are known to be associated with neurodevelopment and mental diseases in human, and are possibly related to the distinct behavior of Silkie. We further identified structural variants in Silkie and found previously reported variants linked to hyperpigmentation (END3), muff and beard (HOXB8), and Rose-comb phenotype (MNR2). Additionally, we found a 0.61 Mb inversion overlapping with the GMDS gene, which was previously linked to neurodevelopmental defects in zebrafish and humans. This may also be related to the behavior distinctiveness of Silkie.

CONCLUSIONS

Our study revealed that Silkie is genetically distinct and relatively highly inbred compared to other IDGB chicken populations, possibly attributed to more prolong population bottlenecks and selective breeding practice. These results enhance our understanding of how domestication and selective breeding have shaped the genome of Silkie. These findings contribute to the broader field of domestication and avian genomics, and have implications for the future conservation and breeding efforts.

The Association between the Abundance of Homozygous Deleterious Variants and the Morbidity of Dog Breeds

It is well known that highly inbred dogs are more prone to diseases than less inbred or outbred dogs. This is because inbreeding increases the load of recessive deleterious variants. Using the genomes of 392 dogs belonging to 83 breeds, we investigated the association between the abundance of homozygous deleterious variants and dog health. We used the number of non-routine veterinary care events for each breed to assess the level of morbidity. Our results revealed a highly significant positive relationship between the number of homozygous deleterious variants located within the runs of homozygosity (RoH) tracts of the breeds and the level of morbidity. The dog breeds with low morbidity had a mean of 87 deleterious SNVs within the RoH, but those with very high morbidity had 187 SNVs. A highly significant correlation was also observed for the loss-of-function (LoF) SNVs within RoH tracts. The dog breeds that required more veterinary care had 2.3 times more homozygous LoF SNVs than those that required less veterinary care (112 vs. 50). The results of this study could be useful for understanding the disease burden on breed dogs and as a guide for dog breeding programs.

The Maintenance of Deleterious Variation in Wild Chinese Rhesus Macaques

Understanding how deleterious variation is shaped and maintained in natural populations is important in conservation and evolutionary biology, as decreased fitness caused by these deleterious mutations can potentially lead to an increase in extinction risk. It is known that demographic processes can influence these patterns. For example, population bottlenecks and inbreeding increase the probability of inheriting identical-by-descent haplotypes from a recent common ancestor, creating long tracts of homozygous genotypes called runs of homozygosity (ROH), which have been associated with an accumulation of mildly deleterious homozygotes. Counterintuitively, positive selection can also maintain deleterious variants in a population through genetic hitchhiking. Here, we analyze the whole genomes of 79 wild Chinese rhesus macaques across five subspecies and characterize patterns of deleterious variation with respect to ROH and signals of recent positive selection. We show that the fraction of homozygotes occurring in long ROH is significantly higher for deleterious homozygotes than tolerated ones, whereas this trend is not observed for short and medium ROH. This confirms that inbreeding, by generating these long tracts of homozygosity, is the main driver of the high burden of homozygous deleterious alleles in wild macaque populations. Furthermore, we show evidence that homozygous LOF variants are being purged. Next, we identify seven deleterious variants at high frequency in regions putatively under selection near genes involved with olfaction and other processes. Our results shed light on how evolutionary processes can shape the distribution of deleterious variation in wild nonhuman primates.

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.

The Maintenance of Deleterious Variation in Wild Chinese Rhesus Macaques

Understanding how deleterious variation is shaped and maintained in natural populations is important in conservation and evolutionary biology, as decreased fitness caused by these deleterious mutations can potentially lead to an increase in extinction risk. It is known that demographic processes can influence these patterns. For example, population bottlenecks and inbreeding increase the probability of inheriting identical-by-descent haplotypes from a recent common ancestor, creating long tracts of homozygous genotypes called runs of homozygosity (ROH), which have been associated with an accumulation of mildly deleterious homozygotes. Counter intuitively, positive selection can also maintain deleterious variants in a population through genetic hitchhiking. Here we analyze the whole genomes of 79 wild Chinese rhesus macaques across five subspecies and characterize patterns of deleterious variation with respect to ROH and signals of recent positive selection. We show that the fraction of homozygotes occurring in long ROH is significantly higher for deleterious homozygotes than tolerated ones, whereas this trend is not observed for short and medium ROH. This confirms that inbreeding, by generating these long tracts of homozygosity, is the main driver of the high burden of homozygous deleterious alleles in wild macaque populations. Furthermore, we show evidence that homozygous LOF variants are being purged. Next, we identify 7 deleterious variants at high frequency in regions putatively under selection near genes involved with olfaction and other processes. Our results shed light on how evolutionary processes can shape the distribution of deleterious variation in wild non-human primates.

Differential patterns in runs of homozygosity in two mice lines under divergent selection for environmental variability for birth weight

Runs of homozygosity (ROH) are defined as long continuous homozygous stretches in the genome which are assumed to originate from a common ancestor. It has been demonstrated that divergent selection for variability in mice is possible and that low variability in birth weight is associated with robustness. To analyse ROH patterns and ROH-based genomic inbreeding, two mouse lines that were divergently selected for birth weight variability for 26 generations were used, with: 752 individuals for the high variability line (H-Line), 766 individuals for the low variability line (L-Line) and 74 individuals as a reference population. Individuals were genotyped using the high density Affymetrix Mouse Diversity Genotyping Array. ROH were identified using both the sliding windows (SW) and the consecutive runs (CR) methods. Inbreeding coefficients were calculated based on pedigree (FPED ) information, on ROH identified using the SW method (FROHSW ) and on ROH identified using the CR method (FROHCR ). Differences in genomic inbreeding were not consistent across generations and these parameters did not show clear differences between lines. Correlations between FPED and FROH were high, particularly for FROHSW . Moreover, correlations between FROHSW and FPED were even higher when ROH were identified with no restrictions in the number of heterozygotes per ROH. The comparison of FROH estimates between either of the selected lines were based on significant differences at the chromosome level, mainly in chromosomes 3, 4, 6, 8, 11, 15 and 19. ROH-based inbreeding estimates that were computed using longer homozygous segments had a higher relationship with FPED . Differences in robustness between lines were not attributable to a higher homozygosis in the L-Line, but maybe to the different distribution of ROH at the chromosome level between lines. The analysis identified a set of genomic regions for future research to establish the genomic basis of robustness.

Modeling the effects of consanguinity on autosomal and X-chromosomal runs of homozygosity and identity-by-descent sharing

Runs of homozygosity (ROH) and identity-by-descent (IBD) sharing can be studied in diploid coalescent models by noting that ROH and IBD-sharing at a genomic site are predicted to be inversely related to coalescence times-which in turn can be mathematically obtained in terms of parameters describing consanguinity rates. Comparing autosomal and X-chromosomal coalescent models, we consider ROH and IBD-sharing in relation to consanguinity that proceeds via multiple forms of first-cousin mating. We predict that across populations with different levels of consanguinity, (1) in a manner that is qualitatively parallel to the increase of autosomal IBD-sharing with autosomal ROH, X-chromosomal IBD-sharing increases with X-chromosomal ROH, owing to the dependence of both quantities on consanguinity levels; (2) even in the absence of consanguinity, X-chromosomal ROH and IBD-sharing levels exceed corresponding values for the autosomes, owing to the smaller population size and lower coalescence time for the X chromosome than for autosomes; (3) with matrilateral consanguinity, the relative increase in ROH and IBD-sharing on the X chromosome compared to the autosomes is greater than in the absence of consanguinity. Examining genome-wide SNPs in human populations for which consanguinity levels have been estimated, we find that autosomal and X-chromosomal ROH and IBD-sharing levels generally accord with the predictions. We find that each 1% increase in autosomal ROH is associated with an increase of 2.1% in X-chromosomal ROH, and each 1% increase in autosomal IBD-sharing is associated with an increase of 1.6% in X-chromosomal IBD-sharing. For each calculation, particularly for ROH, the estimate is reasonably close to the increase of 2% predicted by the population-size difference between autosomes and X chromosomes. The results support the utility of coalescent models for understanding patterns of genomic sharing and their dependence on sex-biased processes.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Investigating the relationship between inbreeding and life expectancy in dogs: mongrels live longer than pure breeds

This study aimed to investigate the establishment of relationship between inbreeding and life expectancy in dogs. A dataset of N = 30,563 dogs sourced from the VetCompass™ Program, UK was made available by the Royal Veterinary College, University of London, containing information about breed and longevity and was subject to survival analysis. A Cox regression proportional hazards model was used to differentiate survivability in three groups of dogs (mongrel, cross-bred and pure breed). The model was found highly significant (p < 0.001) and we found that mongrel dog had the highest life expectancy, followed by cross-bred dogs with only one purebred ancestor and purebred dogs had the lowest life expectancy. A second Cox regression was also found highly significant (p < 0.001) differentiating the lifespan of different dog breed and correlating positively the hazard ratio and the Genetic Illness Severity Index for Dogs (GISID). The results show that survivability is higher in mongrel dogs followed by cross-bred with one of the ancestor only as a pure breed, and pure breed dog have the highest morbidity level. Higher morbidity is associated with higher GISID scores, and therefore, higher inbreeding coefficients. These findings have important implications for dog breeders, owners, and animal welfare organizations seeking to promote healthier, longer-lived dogs.

Genomes of the extinct Bachman's warbler show high divergence and no evidence of admixture with other extant Vermivora warblers

Bachman's warbler¹ (Vermivora bachmanii)-last sighted in 1988-is one of the only North American passerines to recently go extinct.^2,3,4 Given extensive ongoing hybridization of its two extant congeners-the blue-winged warbler (V. cyanoptera) and golden-winged warbler (V. chrysoptera)^5,6,7,8-and shared patterns of plumage variation between Bachman's warbler and hybrids between those extant species, it has been suggested that Bachman's warbler might have also had a component of hybrid ancestry. Here, we use historic DNA (hDNA) and whole genomes of Bachman's warblers collected at the turn of the 20^th century to address this. We combine these data with the two extant Vermivora species to examine patterns of population differentiation, inbreeding, and gene flow. In contrast to the admixture hypothesis, the genomic evidence is consistent with V. bachmanii having been a highly divergent, reproductively isolated species, with no evidence of introgression. We show that these three species have similar levels of runs of homozygosity (ROH), consistent with effects of a small long-term effective population size or population bottlenecks, with one V. bachmanii outlier showing numerous long ROH and a F_ROH greater than 5%. We also found-using population branch statistic estimates-previously undocumented evidence of lineage-specific evolution in V. chrysoptera near a pigmentation gene candidate, CORIN, which is a known modifier of ASIP, which is in turn involved in melanic throat and mask coloration in this family of birds. Together, these genomic results also highlight how natural history collections are such invaluable repositories of information about extant and extinct species.

Comprehensive analysis of geographic and breed-purpose influences on genetic diversity and inherited disease risk in the Doberman dog breed

BACKGROUND

Publicly available phenotype data and genotyping array data from two citizen science projects: "Doberman Health Surveys" and "The Doberman Diversity Project" were analyzed to explore relative homozygosity, diversity, and disorder risk according to geographical locale and breeding purpose in the Doberman.

RESULTS

From the phenotypic data cohort, life expectancy of a Doberman at birth is 9.1 years. The leading causes of death were heart disease (accounting for 28% of deaths) and cancers (collectively accounting for 14% of deaths). By genotyping, the world Doberman population exists as four major cohorts (European exhibition-bred, Americas exhibition-bred, European work, Americas pet/informal). Considering the entire Doberman population, four genomic regions longer than 500 Kb are fixed in 90% or more of 3,226 dogs included in this study. The four fixed regions reside on two autosomal chromosomes: CFA3:0.8-2.3 Mb (1.55 Mb); CFA3: 57.9-59.8 Mb (1.8 Mb); CFA31:0-1.2 Mb (1.2 Mb); and CFA31:4.80-6.47 Mb (1.67 Mb). Using public variant call files including variants for eight Doberman pinschers, we observed 30 potentially functional alternate variants that were evolutionarily diverged relative to the wider sequenced dog population within the four strongly homozygous chromosomal regions. Effective population size (Ne) is a statistical measure of breed diversity at the time of sampling that approximates the number of unique individuals. The major identified sub-populations of Dobermans demonstrated Ne in the range 70-236. The mean level of inbreeding in the Doberman breed is 40% as calculated by the number of array variants in runs of homozygosity divided by the assayed genome size (excluding the X chromosome). The lowest observed level of inbreeding in the Dobermans assayed was 15% in animals that were first generation mixes of European and USA bred Dobermans. Array variant analysis shows that inter-crossing between European and USA-bred Dobermans has capacity to re-introduce variation at many loci that are strongly homozygous.

CONCLUSIONS

We conclude that efforts to improve breed diversity first should focus on regions with the highest fixation levels, but managers must ensure that mutation loads are not worsened by increasing the frequencies of rarer haplotypes in the identified regions. The analysis of global data identified regions of strong fixation that might impact known disorder risks in the breed. Plausible gene candidates for future analysis of the genetic basis of cardiac disease and cancer were identified in the analysis.

Selection, recombination and population history effects on runs of homozygosity (ROH) in wild red deer (Cervus elaphus)

The distribution of runs of homozygosity (ROH) may be shaped by a number of interacting processes such as selection, recombination and population history, but little is known about the importance of these mechanisms in shaping ROH in wild populations. We combined an empirical dataset of >3000 red deer genotyped at >35,000 genome-wide autosomal SNPs and evolutionary simulations to investigate the influence of each of these factors on ROH. We assessed ROH in a focal and comparison population to investigate the effect of population history. We investigated the role of recombination using both a physical map and a genetic linkage map to search for ROH. We found differences in ROH distribution between both populations and map types indicating that population history and local recombination rate have an effect on ROH. Finally, we ran forward genetic simulations with varying population histories, recombination rates and levels of selection, allowing us to further interpret our empirical data. These simulations showed that population history has a greater effect on ROH distribution than either recombination or selection. We further show that selection can cause genomic regions where ROH is common, only when the effective population size (Ne) is large or selection is particularly strong. In populations having undergone a population bottleneck, genetic drift can outweigh the effect of selection. Overall, we conclude that in this population, genetic drift resulting from a historical population bottleneck is most likely to have resulted in the observed ROH distribution, with selection possibly playing a minor role.

How size and genetic diversity shape lifespan across breeds of purebred dogs

While the lifespan advantage of small body size and mixed breed status has been documented repeatedly, evidence for an effect of genetic diversity across dog breeds is equivocal. We hypothesized that this might be due to a strong right-censoring bias in available breed-specific lifespan estimates where early-dying dogs from birth cohorts that have not died off completely at the time of data collection are sampled disproportionately, especially in breeds with rapidly growing populations. We took advantage of data on owner reported lifespan and cause of death from a large public database to quantify the effect of size and genetic diversity (heterozygosity) on mortality patterns across 118 breeds based on more than 40,000 dogs. After documenting and removing the right-censoring bias from the breed-specific lifespan estimates by including only completed birth cohorts in our analyses, we show that small size and genetic diversity are both linked to a significant increase in mean lifespan across breeds. To better understand the proximate mechanisms underlying these patterns, we then investigated two major mortality causes in dogs - the cumulative pathophysiologies of old age and cancer. Old age lifespan, as well as the percentage of old age mortality, decreased with increasing body size and increased with increasing genetic diversity. The lifespan of dogs dying of cancer followed the same patterns, but while large size significantly increased proportional cancer mortality, we could not detect a significant signal for lowered cancer mortality with increasing diversity. Our findings suggest that outcross programs will be beneficial for breed health and longevity. They also emphasize the need for high-quality mortality data for veterinary epidemiology as well as for developing the dog as a translational model for human geroscience.

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.

Genome scanning of behavioral selection in a canine olfactory detection breeding cohort

Research on working dogs is growing rapidly due to increasing global demand. Here we report genome scanning of the risk of puppies being eliminated for behavioral reasons prior to entering the training phase of the US Transportation Security Administration's (TSA) canine olfactory detection breeding and training program through 2013. Elimination of dogs for behavioral rather than medical reasons was based on evaluations at three, six, nine and twelve months after birth. Throughout that period, the fostered dogs underwent standardized behavioral tests at TSA facilities, and, for a subset of tests, dogs were tested in four different environments. Using methods developed for family studies, we performed a case-control genome wide association study (GWAS) of elimination due to behavioral observation and testing results in a cohort of 528 Labrador Retrievers (2002-2013). We accounted for relatedness by including the pedigree as a covariate and maximized power by including individuals with phenotype, but not genotype, data (approximately half of this cohort). We determined genome wide significance based on Bonferroni adjustment of two quasi-likelihood score tests optimized for either small or nearly-fully penetrant effect sizes. Six loci were significant and five suggestive, with approximately equal numbers of loci for the two tests and frequencies of loci with single versus multiple mapped markers. Several loci implicate a single gene, including CHD2, NRG3 and PDE1A which have strong relevance to behavior in humans and other species. We briefly discuss how expanded studies of canine breeding programs could advance understanding of learning and performance in the mammalian life course. Although human interactions and other environmental conditions will remain critical, our findings suggest genomic breeding selection could help improve working dog populations.

Inheritance of Monogenic Hereditary Skin Disease and Related Canine Breeds

The plasticity of the genome is an evolutionary factor in all animal species, including canines, but it can also be the origin of diseases caused by hereditary genetic mutation. Genetic changes, or mutations, that give rise to a pathology in most cases result from recessive alleles that are normally found with minority allelic frequency. The use of genetic improvement increases the consanguinity within canine breeds and, on many occasions, also increases the frequency of these recessive alleles, increasing the prevalence of these pathologies. This prevalence has been known for a long time, but mutations differ according to the canine breed. These genetic diseases, including skin diseases, or genodermatosis, which is narrowly defined as monogenic hereditary dermatosis. In this review, we focus on genodermatosis sensu estricto, i.e., monogenic, and hereditary dermatosis, in addition to the clinical features, diagnosis, pathogeny, and treatment. Specifically, this review analyzes epidermolytic and non-epidermolytic ichthyosis, junctional epidermolysis bullosa, nasal parakeratosis, mucinosis, dermoid sinus, among others, in canine breeds, such as Golden Retriever, German Pointer, Australian Shepherd, American Bulldog, Great Dane, Jack Russell Terrier, Labrador Retriever, Shar-Pei, and Rhodesian Ridgeback.

The impact of identity by descent on fitness and disease in dogs

Dogs and humans have coexisted together for thousands of years, but it was not until the Victorian Era that humans practiced selective breeding to produce the modern standards we see today. Strong artificial selection during the breed formation period has simplified the genetic architecture of complex traits and caused an enrichment of identity-by-descent (IBD) segments in the dog genome. This study demonstrates the value of IBD segments and utilizes them to infer the recent demography of canids, predict case-control status for complex traits, locate regions of the genome potentially linked to inbreeding depression, and to identify understudied breeds where there is potential to discover new disease-associated variants.

Domestic dogs have experienced population bottlenecks, recent inbreeding, and strong artificial selection. These processes have simplified the genetic architecture of complex traits, allowed deleterious variation to persist, and increased both identity-by-descent (IBD) segments and runs of homozygosity (ROH). As such, dogs provide an excellent model for examining how these evolutionary processes influence disease. We assembled a dataset containing 4,414 breed dogs, 327 village dogs, and 380 wolves genotyped at 117,288 markers and data for clinical and morphological phenotypes. Breed dogs have an enrichment of IBD and ROH, relative to both village dogs and wolves, and we use these patterns to show that breed dogs have experienced differing severities of bottlenecks in their recent past. We then found that ROH burden is associated with phenotypes in breed dogs, such as lymphoma. We next test the prediction that breeds with greater ROH have more disease alleles reported in the Online Mendelian Inheritance in Animals (OMIA). Surprisingly, the number of causal variants identified correlates with the popularity of that breed rather than the ROH or IBD burden, suggesting an ascertainment bias in OMIA. Lastly, we use the distribution of ROH across the genome to identify genes with depletions of ROH as potential hotspots for inbreeding depression and find multiple exons where ROH are never observed. Our results suggest that inbreeding has played a large role in shaping genetic and phenotypic variation in dogs and that future work on understudied breeds may reveal new disease-causing variation.