Genome-wide diversity and runs of homozygosity in the "Braque Français, type Pyrénées" dog breed

Genetic diversity and origin of Kazakh Tobet Dogs

The Kazakh Tobet is an indigenous Kazakh dog breed that has been used to guard livestock since ancient times. To understand the genetic structure and phylogenetic relationship of the Kazakh Tobet breed with other herding and livestock guarding dog breeds, we analysed short tandem repeat data of 107 Kazakh Tobet dogs from different regions of Kazakhstan and Mongolia, as well as whole genome sequencing data from two Kazakh Tobet dogs and 43 dogs from 24 working breeds. Our results indicate a high genetic diversity of the Kazakh Tobet, with the average number of alleles per locus ranging from 6.00 to 10.22 and observed heterozygosity ranging from 76 to 78%. The breed has a complex genetic structure characterised by seven different clusters. The neighbour-joining tree constructed based on 14,668,406 autosomal and the maximum likelihood tree based on mitochondrial D-loop sequences indicate a common genetic heritage between the Kazakh Tobet, the Central Asian Shepherd Dog and the Turkish Akbash. The presence of haplotype A18 in the Kazakh Tobets supports the hypothesis of the ancient origin of the breed, which was previously suggested by archaeological finds and written sources. These results provide an important genetic basis for the ongoing efforts to improve the Kazakh Tobet breed, to ensure its preservation as an independent genetic lineage and to recognise a breed on an international level.

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.

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.

Traces of Human-Mediated Selection in the Gene Pool of Red Deer Populations

In this study, we analysed the effect of human-mediated selection on the gene pool of wild and farmed red deer populations based on genotyping-by-sequencing data. The farmed red deer sample covered populations spread across seven countries and two continents (France, Germany, Hungary, Latvia, New Zealand, Poland, and Slovakia). The Slovak and Spain wild red deer populations (the latter one in a large game estate) were used as control outgroups. The gene flow intensity, relationship and admixture among populations were tested by the Bayesian approach and discriminant analysis of principal components (DAPC). The highest gene diversity (He = 0.19) and the lowest genomic inbreeding (FHOM = 0.04) found in Slovak wild population confirmed our hypothesis that artificial selection accompanied by bottlenecks has led to the increase in overall genomic homozygosity. The Bayesian approach and DAPC consistently identified three separate genetic groups. As expected, the farmed populations were clustered together, while the Slovak and Spanish populations formed two separate clusters. Identified traces of genetic admixture in the gene pool of farmed populations reflected a strong contemporary migration rate between them. This study suggests that even if the history of deer farming has been shorter than traditional livestock species, it may leave significant traces in the genome structure.

Homozygosity mapping in the Kazakh national dog breed Tazy

The Tazy is a breed of sighthound common in Kazakhstan. The identification of runs of homozygosity (ROH) is an informative approach to assessing the history and possible patterns of directional selection pressure. To our knowledge, the present study is the first to provide an overview of the ROH pattern in the Tazy dogs from a genome-wide perspective. The ROH of the Tazy was found to be mainly composed of shorter segments (1-2 Mb), accounting for approximately 67% of the total ROH. The estimated ROH-based inbreeding coefficients (F_ROH) ranged from 0.028 to 0.058 with a mean of 0.057. Five genomic regions under positive selection were identified on chromosomes 18, 22, and 25. The regions on chromosomes 18 and 22 may be breed specific, while the region on chromosome 22 overlaps with regions of hunting traits in other hunting dog breeds. Among the 12 candidate genes located in these regions, the gene CAB39L may be a candidate that affects running speed and endurance of the Tazy dog. Eight genes could belong to an evolutionarily conserved complex as they were clustered in a large protein network with strong linkages. The results may enable effective interventions when incorporated into conservation planning and selection of the Tazy breed.

New Insight into the Genome-Wide Diversity and Admixture of Six Colombian Sheep Populations

Creole sheep represent a strategic genetic resource for populations living in marginal areas under financial restrictions on the American continent. Six Colombian sheep breeds (two wool (BCL-Boyacá and NCL-Nariño, 12 and 14 samples) and four hair (OPCE-Ethiopian, 54 samples; OPCS-Sudan, 74 samples; OPCP-Pelibeuy, 59 samples; OPCW-Wayúu, 24 samples) were genotyped using the Illumina Ovine SNP50 BeadChip. Data was also included from international 44 breeds from International Sheep Genomics Consortium (ISGC) and from data published in previous a previous work on the Caribbean and African breeds. Although geographically separated, wool (NCL, BCL) and hair types (OPCE, OPCS, OPCW) presented little genetic differentiation (FST 0.05) at a global level but several groups of animals separated suggesting local clustering due to geographical isolation. The OPCP underwent a recent crossing with Mexican Pelibuey, explaining its differentiation. Findings in this work such as the proximity to West African Djallonké (WAD) and Barbados Black Belly (BBB), suggest different introductions of African type animals from the Caribbean region on a pre-existing genetic basis formed by animals deriving from the first importations coming from Europe in colonial times. As expected, Colombian wool breeds showed, in particular in Admixture software results, a greater genomic component in common with European breeds and in particular with Iberian ones (Churra). This study provides a basis for future research into the genetic diversity within and between the Colombian sheep breeds analysed, and scientific data for policy decisions on Farm Animal Genetic Resources (FAnGR).

Evaluation of genetic diversity and management of disease in Border Collie dogs

Maintaining genetic diversity in dog breeds is an important consideration for the management of inherited diseases. We evaluated genetic diversity in Border Collies using molecular and genealogical methods, and examined changes to genetic diversity when carriers for Trapped Neutrophil Syndrome (TNS) and Neuronal Ceroid Lipofuscinosis (NCL) are removed from the genotyped population. Genotype data for 255 Border Collies and a pedigree database of 83,996 Border Collies were used for analysis. Molecular estimates revealed a mean multi-locus heterozygosity (MLH) of 0.311 (SD 0.027), 20.79% of the genome consisted of runs of homozygosity (ROH ) > 1 Mb, effective population size (N_e) was 84.7, and mean inbreeding (F) was 0.052 (SD 0.083). For 227 genotyped Border Collies that had available pedigree information (GenoPed), molecular and pedigree estimates of diversity were compared. A reference population (dogs born between 2005 and 2015, inclusive; N = 13,523; RefPop) and their ancestors (N = 12,478) were used to evaluate the diversity of the population that are contributing to the current generation. The reference population had a N_e of 123.5, a mean F of 0.095 (SD 0.082), 2276 founders (f), 205.5 effective founders (f_e), 28 effective ancestors (f_a) and 10.65 (SD 2.82) founder genomes (N_g). Removing TNS and NCL carriers from the genotyped population had a small impact on diversity measures (ROH > 1 Mb, MLH, heterozygosity), however, there was a loss of > 10% minor allele frequency for 89 SNPs around the TNS mutation (maximum loss of 12.7%), and a loss of > 5% for 5 SNPs around the NCL mutation (maximum 5.18%). A common ancestor was identified for 38 TNS-affected dogs and 64 TNS carriers, and a different common ancestor was identified for 33 NCL-affected dogs and 28 carriers, with some overlap of prominent individuals between both pedigrees. Overall, Border Collies have a high level of genetic diversity compared to other breeds.

Inbreeding levels in an open-registry pedigreed dog breed: The Australian working kelpie

The depletion in genetic diversity of closed-pedigree dog breeds can be a contentious topic and can lead to calls for open-registry. However, strong evidence in support of proposed open-registry solutions is lacking, with the reproductive isolation of these breeds unlikely to be the sole cause of elevated inbreeding levels. Human-induced limitations, such as popular sire effects, are unlikely to be confined to closed-registry breeds and conceivably play an important role in maintaining genetic diversity within all breeds. Consequently, the aim of the current study was to explore inbreeding levels in an open-registry breed and determine the impact open-registry has on genetic diversity. Complete pedigree records on all Australian working kelpies (AWKs) were provided by the Working Kelpie Council with the cleaned pedigree consisting of 86,671 individuals with a median pedigree depth of 6.6 generations. The average inbreeding coefficient in the AWK population was 0.049 with an increase in inbreeding coefficient of 0.0016/year. This demonstrates that opening a breed registry can have a beneficial impact on the level of inbreeding within a population over the longer-term. However, allowing for a generation length of 5.1 years yielded an effective population size of 61 for AWKs and demonstrated a pattern consistent with closed-registry dog populations of comparable size.

Does Functionality Condition the Population Structure and Genetic Diversity of Endangered Dog Breeds under Island Territorial Isolation?

Simple Summary

Early references to Ca de Rater and Ca de Bestiar endangered autochthonous breeds were witnesses to their exceptional ratting/pet and shepherding/guard skills for centuries. Studbooks or associations promote an increase in the number of effectives and their genealogies. Genetic diversity parameters were evaluated along the history of definition of both dog breeds. Guard or hunting purposes condition an increased registration of genealogical information. Hunting animals have more complete genealogies and need more time to select breeding animals. Male guard dogs are preferred over females due to their suitability for guarding. Selection for performance acts as a diversity promoter and breeding policy driver. The uses or purposes for which certain breeds were selected condition the genetic diversity evolution of endangered breeds, even if these share the same geographic isolation conditions.

Abstract

Despite the undefinition of the origins of Ca de Rater (CR) and Ca de Bestiar (CB) dogs, references to these endangered autochthonous breeds highlighted their ratting/pet and shepherding/guard skills for centuries. Genealogical historical records were traced back to founders. Founder number in the reference population (146 and 53 for CR and CB, respectively), historical and reference maximum generations traced (eight and seven for CR and CB, respectively), and historical average number of complete generations (1.04 for both breeds) were determined. Structure assessment revealed the existence of subpopulations regarding criteria such as breeders (75 and 17), breeder location (32 and eight), owners (368 and 198), and owner location (73 and 51) for CR and CB, respectively. Average inbreeding (F) within breed subpopulations ranged from 0.27–1.20% for CB breeders and the rest of subpopulation criteria for both breeds, respectively, except for CB owners and owner location. F ranged from 0.27–1.41% for CB historical population and CR current population, respectively. The study of genetic diversity revealed a relatively similar genetic background between subpopulations. Average coancestry between and within breeds suggested a similar evolutionary process. However, Mann–Whitney U test determined significant differences for diversity parameters (F, ΔR, coancestry, nonrandom mating degree, maximum, complete, and equivalent generations, ΔF, and genetic conservation index) between breeds and their functionalities. Conclusively, functionality in dog breeds may determine the genetic diversity evolution of endangered breeds, even when these share the same geographic isolation conditions.

Analysis of Genetic Diversity in the Czech Spotted Dog

Loss off genetic diversity negatively affects most of the modern dog breeds. However, no breed created strictly for laboratory purposes has been analyzed so far. In this paper, we sought to explore by pedigree analysis exactly such a breed-the Czech Spotted Dog (CSD). The pedigree contained a total of 2010 individuals registered since the second half of the 20th century. Parameters such as the mean average relatedness, coefficient of inbreeding, effective population size, effective number of founders, ancestors and founder genomes and loss of genetic diversity-which was calculated based on the reference population and pedigree completeness-were used to assess genetic variability. Compared to the founding population, the reference population lost 38.2% of its genetic diversity, of which 26% is due to random genetic drift and 12.2% is due to the uneven contribution of the founders. The reference population is highly inbred and related. The average inbreeding coefficient is 36.45%, and the mean average relatedness is 74.83%. The effective population size calculated based on the increase of inbreeding coefficient is 10.28. Thus, the Czech Spotted Dog suffered significant losses of genetic diversity that threaten its future existence.

The German Shorthair Pointer Dog Breed (Canis lupus familiaris): Genomic Inbreeding and Variability

The German Shorthaired Pointer (GSHP) is a breed worldwide known for its hunting versatility. Dogs of this breed are appreciated as valuable companions, effective trackers, field trailers and obedience athletes. The aim of the present work is to describe the genomic architecture of the GSHP breed and to analyze inbreeding levels under a genomic and a genealogic perspective. A total of 34 samples were collected (24 Italian, 10 USA), and the genomic and pedigree coefficients of inbreeding have been calculated. A total of 3183 runs of homozygosity (ROH) across all 34 dogs have been identified. The minimum and maximum number of Single Nucleotide Polymorphisms (SNPs) defining all ROH are 40 and 3060. The mean number of ROH for the sample was 93.6. ROH were found on all chromosomes. A total of 854 SNPs (TOP_SNPs) defined 11 ROH island regions (TOP_ROH), in which some gene already associated with behavioral and morphological canine traits was annotated. The proportion of averaged observed homozygotes estimated on total number of SNPs was 0.70. The genomic inbreeding coefficient based on ROH was 0.17. The mean inbreeding based on genealogical information resulted 0.023. The results describe a low inbred population with quite a good level of genetic variability.

True Colors: Commercially-acquired morphological genotypes reveal hidden allele variation among dog breeds, informing both trait ancestry and breed potential

Direct-to-consumer canine genetic testing is becoming increasingly popular among dog owners. The data collected therein provides intriguing insight into the current status of morphological variation present within purebred populations. Mars WISDOM PANELTM data from 11,790 anonymized dogs, representing 212 breeds and 4 wild canine species, were evaluated at genes associated with 7 coat color traits and 5 physical characteristics. Frequencies for all tested alleles at these 12 genes were determined by breed and by phylogenetic grouping. A sub-set of the data, consisting of 30 breeds, was divided into separate same-breed populations based on country of collection, body size, coat variation, or lineages selected for working or conformation traits. Significantly different (p ≤ 0.00167) allele frequencies were observed between populations for at least one of the tested genes in 26 of the 30 breeds. Next, standard breed descriptions from major American and international registries were used to determine colors and tail lengths (e.g. genetic bobtail) accepted within each breed. Alleles capable of producing traits incongruous with breed descriptions were observed in 143 breeds, such that random mating within breeds has probabilities of between 4.9e-7 and 0.25 of creating undesirable phenotypes. Finally, the presence of rare alleles within breeds, such as those for the recessive black coloration and natural bobtail, was combined with previously published identity-by-decent haplotype sharing levels to propose pathways by which the alleles may have spread throughout dog breeds. Taken together, this work demonstrates that: 1) the occurrence of low frequency alleles within breeds can reveal the influence of regional or functional selection practices; 2) it is possible to visualize the potential historic connections between breeds that share rare alleles; and 3) the necessity of addressing conflicting ideals in breed descriptions relative to actual genetic potential is crucial.

A Genome-Wide Detection of Copy Number Variations Using SNP Genotyping Arrays in Braque Français Type Pyrénées Dogs

Copy number variants (CNVs) are an important source of genetic variation complementary to single nucleotide polymorphisms (SNPs). Only few studies have been conducted in dogs on CNVs derived from high-density SNP array data, and many canine breeds still remain uncharacterized, e.g., the Braque Français, type Pyrénées breed (BRA). Therefore, in an effort to more comprehensively investigate the canine genome for CNVs, we used a high-density SNP array (170 K) to discover CNVs in BRA. The CNV regions (CNVRs) were identified through the merging of two different CNVRs datasets, obtained separately from SNP data using the PennCNV and SVS software. A total of 45 stringent CNVRs, ranging from 3.5 kb to 458,716 kb in length were detected in 26 dog samples. Results overlapped moderately in comparison with previous studies on CNVs in dogs, leading to the identification of 16 novel CNVRs. A total of 159 genes were annotated in the CNVRs detected with stringent quality criteria in particular high classification stringency and false discovery rate correction. The gene ontology enrichment analysis provided information on biological processes and cellular components related to muscle structure development and muscle cell differentiation. Considering that BRA is a breed used for speed in hunting and retrieval, for the ability to find feathered game, and for pointing, we can hypothesize that selection for such hunting behavior could have driven, at least in part, the presence of these genes into the CNVRs.

Heterozygosity testing and multiplex DNA panel screening as a potential tool to monitor health and inbreeding in a small, closed dog population

Background

Selective breeding in populations with a limited effective population size may result in a loss of genetic diversity, which can cause an increased concentration of specific disease liability genes. The Dutch Shepherd Dog (DSD) in the Netherlands is an example of such a breed with a small effective population.

Objective

To evaluate the measurement of genetic diversity and multiplex DNA panel screening for implementation in a breeding strategy for the Dutch Shepherd Dog (DSD) and to investigate the clinical relevance of potentially identified mutations in the multiplex DNA panel screening.

Results

Genome-wide SNP testing showed genetic isolation and reduced genetic diversity within coat variety subgroups of the DSD. Panel screening identified a Von Willebrand's Disease type I mutation. Although decreased Von Willebrand's Factor proteins were significantly lower in DSDs carrying the VWD-I allele compared to the wildtype, clinical follow-up did not show a significant association between the clinical phenotype and VWD-I genotype.

Conclusions

Genetic relationship measurement within a breed population may be a useful tool to enable breeding strategies to conserve genetic diversity. Results from a disease panel screening need to be evaluated for clinical relevance before breed selection restrictions can be considered.

Genomic characterization of the Braque Français type Pyrénées dog and relationship with other breeds

The evaluation of genetic variability is a useful research tool for the correct management of selection and conservation strategies in dog breeds. In addition to pedigree genealogies, genomic data allow a deeper knowledge of the variability and genetic structure of populations. To date, many dog breeds, such as small regional breeds, still remain uncharacterized. Braque Français type Pyrénées (BRA) is a dog breed originating from a very old type of gun-dog used for pointing the location of game birds to hunters. Despite the ancient background, the knowledge about levels of genetic diversity, degree of inbreeding and population structure is scarce. This may raise concerns on the possibility that few inbred bloodlines may dominate the breed, and on its future health. The aim of this work was therefore to provide a high-resolution representation of the genome-wide diversity and population structure of BRA dogs, using the 170K genome-wide SNP array. Genome-wide polymorphisms in BRA were compared with those of other worldwide dog breeds. Between-dog relationships estimated from genomic data were very similar to pedigree relationships (Pearson correlation r_g,a = 0.92). Results showed that BRA generally presents moderate levels of genetic diversity when compared with the major canine breeds. The estimated effective population size (recent N_e = 51) shows a similar declining pattern over generations as all other dog breeds, pointing at a common demographic history of modern canine breeds, clearly different from the demography of feral wolves. Multidimensional scaling (MDS), Bayesian clustering and Neighbor Joining tree were used to visualize and explore the genetic relationships among breeds, and revealed that BRA was highly differentiated and presented only low levels of admixture with other breeds. Brittany Spaniel, English Setter, Gordon Setter and Weimaraner dogs are the closest breeds to BRA. The exact reason for BRA being so divergent from other dog breeds, based on these results, is not yet clear. Further studies including additional ≪braccoid≫ breeds will be needed to refine the results presented here and to investigate the origin of the BRA breed. Nonetheless, the genome-wide characterization reported here provides a comprehensive insight into the genome diversity and population structure of the Braque Français, type Pyrénées breed.

Genetic diversity and population structure of African village dogs based on microsatellite and immunity-related molecular markers

The village and street dogs represent a unique model of canine populations. In the absence of selective breeding and veterinary care, they are subject mostly to natural selection. Their analyses contribute to understanding general mechanisms governing the genetic diversity, evolution and adaptation. In this study, we analyzed the genetic diversity and population structure of African village dogs living in villages in three different geographical areas in Northern Kenya. Data obtained for neutral microsatellite molecular markers were compared with those computed for potentially non-neutral markers of candidate immunity-related genes. The neutral genetic diversity was similar to other comparable village dog populations studied so far. The overall genetic diversity in microsatellites was higher than the diversity of European pure breeds, but it was similar to the range of diversity observed in a group composed of many European breeds, indicating that the African population has maintained a large proportion of the genetic diversity of the canine species as a whole. Microsatellite marker diversity indicated that the entire population is subdivided into three genetically distinct, although closely related subpopulations. This genetical partitioning corresponded to their geographical separation and the observed gene flow well correlated with the communication patterns among the three localities. In contrast to neutral microsatellites, the genetic diversity in immunity-related candidate SNP markers was similar across all three subpopulations and to the European group. It seems that the genetic structure of this particular population of Kenyan village dogs is mostly determined by geographical and anthropogenic factors influencing the gene flow between various subpopulations rather than by biological factors, such as genetic contribution of original migrating populations and/or the pathogen-mediated selection. On the other hand, the study of oldest surviving dogs suggested a biological mechanism, i.e. a possible advantage of the overal heterozygosity marked by the the microsatellite loci analyzed.