Genetic diversity of dog breeds: within-breed diversity comparing genealogical and molecular data

Genomic inbreeding measures applied to a population of mice divergently selected for birth weight environmental variance

This study aimed to compare different inbreeding measures estimated from pedigree and molecular data from two divergent mouse lines selected for environmental birth weight during 26 generations. Furthermore, the performance of different approaches and both molecular and pedigree data sources for estimating Ne were tested in this population. A total of 1,699 individuals were genotyped using a high-density genotyping array. Genomic relationship matrices were used to calculate molecular inbreeding: Nejati-Javaremi (F NEJ), Li and Horvitz (F L&H), Van Raden method 1 (F VR1) and method 2 (F VR2), and Yang (F YAN). Inbreeding based on runs of homozygosity (F ROH) and pedigree inbreeding (F PED) were also computed. F ROH, F NEJ, and F L&H were also adjusted for their average values in the first generation of selection and named F ROH0, F NEJ0, and F L&H0. ∆F was calculated from pedigrees as the individual inbreeding rate between the individual and his parents (∆F PEDt) and individual increases in inbreeding (∆F PEDi). Moreover, individual ∆F was calculated from the different molecular inbreeding coefficients (∆F NEJ0, ∆F L&H, ∆F L&H0, ∆F VR1, ∆F VR2, ∆F YAN, and ∆F ROH0). The Ne was obtained from different ∆F, such as Ne PEDt, Ne PEDi, Ne NEJ0, Ne L&H, Ne L&H0, Ne VR1, Ne VR2, Ne YAN, and Ne ROH0. Comparing with F PED , F ROH , F NEJ and F VR2 overestimated inbreeding while F NEJ0 , F L&H , F L&H0 , F VR1 and F YAN underestimated inbreeding. Correlations between inbreeding coefficients and ∆F were calculated. F ROH had the highest correlation with F PED (0.89); F YAN had correlations >0.95 with all the other molecular inbreeding coefficients. Ne PEDi was more reliable than Ne PEDt and presented similar behaviour to Ne L&H0 and Ne NEJ0. Stable trends in Ne were not observed until the 10th generation. In the 10th generation Ne PEDi was 42.20, Ne L&H0 was 45.04 and Ne NEJ0 was 45.05 and in the last generation these Ne were 35.65, 35.94 and 35.93, respectively F ROH presented the highest correlation with F PED, which addresses the identity by descent probability (IBD). The evolution of Ne L&H0 and Ne NEJ0 was the most similar to that of Ne PEDi. Data from several generations was necessary to reach a stable trend for Ne, both with pedigree and molecular data. This population was useful to test different approaches to computing inbreeding coefficients and Ne using molecular and pedigree data.

A purebred South American breed showing high effective population size and independent breed ancestry: The Chilean Terrier

The Chilean Terrier is a known breed in Chile that has not been genetically assessed despite its distinctive color patterns, agility, and hardiness across the diversity of climates encountered within the Chilean landscape. The population structure and its relatedness with other breeds, as well as the actual origin of the breed, remain unknown. We estimated several population parameters using samples from individuals representing the distribution of the Chilean Terrier across the country. By utilizing the Illumina HD canine genotyping array, we computed the effective population size (Ne ), individual inbreeding, and relatedness to evaluate the genetic diversity of the breed. The results show that linkage disequilibrium was relatively low and decayed rapidly; in fact, Ne was very high when compared to other breeds, and similar to other American indigenous breeds (such as the Chihuahua with values of Ne near 500). These results are in line with the low estimates of genomic inbreeding and relatedness and the relatively large number of effective chromosome segments (Me  = 2467) obtained using the properties of the genomic relationship matrix. Between population analysis (cross-population extended haplotype homozygosity, di ) with other breeds such as the Jack Russell Terrier, the Peruvian-Inca Orchid, and the Chihuahua suggested that candidate regions harboring FGF5, PAX3, and ASIP, probably explained some morphological traits, such as the distinctive color pattern characteristic of the breed. When considering Admixture estimates and phylogenetic analysis, together with other breeds of American and European origin, the Chilean Terrier does not have a recent European ancestry. Overall, the results suggest that the breed has evolved independently in Chile from other terrier breeds, from an unknown European terrier ancestor.

Incorporating Genetic Testing into a Breeding Program

Genetic tests are powerful tools that enable (1) a focus on genetic diversity as mating outcomes can be predicted and thus optimized to minimize or even avoid exclusion and (2) working toward breeding goals by improving a phenotype.

Genetic Diversity and Trends of Ancestral and New Inbreeding in German Sheep Breeds by Pedigree Data

In Germany, many autochthonous sheep breeds have developed, adapted to mountain, heath, moorland, or other marginal sites, but breeds imported from other countries have also contributed to the domestic breeds, particularly improving wool and meat quality. Selective breeding and the intense use of rams may risk losing genetic diversity and increasing rates of inbreeding. On the other hand, breeds with a low number of founder animals and only regional popularity may not leave their endangered status, as the number of breeders interested in the breed is limited. The objective of the present study was to determine demographic measures of genetic diversity and recent as well as ancestral trends of inbreeding in all autochthonous German sheep breeds and sheep of all breeding directions, including wool, meat, and milk. We used pedigree data from 1,435,562 sheep of 35 different breeds and a reference population of 981,093 sheep, born from 2010 to 2020. The mean number of equivalent generations, founders, effective founders, effective ancestors, and effective founder genomes were 5.77, 1669, 123.2, 63.5, and 33.0, respectively. Genetic drift accounted for 69% of the loss of genetic diversity, while loss due to unequal founder contributions was 31%. The mean inbreeding coefficient, individual rate of inbreeding (∆F_i), and realized effective population size across breeds were 0.031, 0.0074, and 91.4, respectively, with a significantly decreasing trend in ∆F_i in 11/35 breeds. New inbreeding, according to Kalinowski, contributed to 71.8% of individual inbreeding, but ancestral inbreeding coefficients showed an increasing trend in all breeds. In conclusion, in our study, all but one of the mountain-stone sheep breeds and the country sheep breed Wald were the most vulnerable populations, with N_e < 50. The next most endangered breeds are exotic, country, and heath breeds, with average N_e of 66, 83, and 89, respectively. The wool, meat, and milk breeds showed the highest genetic diversity, with average N_e of 158, 120, and 111, respectively. The results of our study should help strengthen conservation program efforts for the most endangered sheep breeds and maintain a high genetic diversity in all sheep breeds.

Phylogeographic and population genetic structure of hound-like native dogs of the Mediterranean Basin

The dog was probably the first domesticated animal. Despite extensive archaeological and genetic investigations, the origin and the evolution of the extant dogs are still being debated. Dog breeds that have over time been selected for hunting share common ancestral traits. This study represents the first comprehensive attempt to survey at the genomic and mitochondrial level eight hound-like dogs breeds indigenous to the Mediterranean Basin to determine if they share common ancient origins. Results from the microsatellite analysis indicate that all the dog populations have a low inbreeding value.The Kelb tal-Fenek has a high divergence from the current Egyptian street population, however there is not enough evidence from this study to exclude completely the potential of an ancient common relationship. Overall, the mitochondrial results indicate high frequencies of haplogroups A and B and a low representation of haplogroup C, while only one Egyptian dog could be assigned to haplogroup D. Results reveal identities and shared clades, suggesting the conservation of ancient European mitotypes in the Mediterranean hound-like breeds, especially in the Egyptian population. Although none of the dog populations/breeds participating in this study indicate to be direct descendants of the Egyptian dogs, they still have a very close morphologically resemblance to those iconic Egyptian dogs often depicted in ancient art forms and share some genetic links with the current Egyptian population. Further research is required with other markers such us complete mitogenomes and SNP panels to confirm the complex history of the Mediterranean dogs involved in this study.

Genetic Diversity and Trends of Ancestral and New Inbreeding in Deutsch Drahthaar Assessed by Pedigree Data

Simple Summary

Deutsch Drahthaar (DD) is the most popular hunting dog in Germany, fulfilling all aspects of hunting including searching for trails. This breed was newly created at the beginning of the 20th century from a large number existing versatile hunting dog breeds. The aim of the breed was, and still is, to achieve the best performance in all aspects of hunting. We analyzed pedigrees of DD using demographic measures to quantify genetic diversity such as probabilities of gene origin and degrees of ancestral and individual inbreeding. A large number of genetically diverse founder dogs should open up the opportunity of creating a breed with a high genetic diversity and a low increase of inbreeding per generation. On the other hand, intense use of top sires and dams from a limited number of breeding lines may accelerate breeding progress in hunting abilities but reduce genetic diversity. Monitoring genetic diversity should help to maintain a high diversity of breeding populations. Our analysis of pedigree data from 101,887 DD dogs revealed inbreeding measures (coefficient of inbreeding F = 0.042, individual rate of inbreeding ΔF_i = 0.00551) and effective population size (Ne = 92) in the mean range compared to a wide range of other dog breeds. Ancestral inbreeding had a strong increasing trend, whereas trends in individual inbreeding and rate of individual inbreeding were slightly negative.

Abstract

Loss of genetic diversity and high inbreeding rates confer an increased risk of congenital anomalies and diseases and thus impacting dog breeding. In this study, we analyzed recent and ancestral inbreeding as well as other measures of genetic variability in the Deutsch Drahthaar (DD) dog population. Analyses included pedigree data from 101,887 animals and a reference population with 65,927 dogs born between 2000 and 2020. The mean equivalent complete generations was 8.6 with 69% known ancestors in generation 8. The mean realized effective population size was 92 with an increasing trend from 83 to 108 over birth years. The numbers of founders, effective founders and effective ancestors, as well as founder genomes, were 814, 66, 38 and 16.15, respectively. Thirteen ancestors explained 50% of the genetic diversity. The mean coefficient of inbreeding and individual rate of inbreeding (ΔF_i) were 0.042 and 0.00551, respectively, with a slightly decreasing trend in ΔF_i. Exposure of ancestors to identical-by-descent alleles explored through ancestral coefficients of inbreeding showed a strong increasing trend. Comparisons between new and ancestral inbreeding coefficients according to Kalinowski et al. showed an average relative contribution of 62% of new inbreeding to individual inbreeding. Comparisons among average coancestry within the parental population and average inbreeding in the reference population were not indicative of genetic substructures. In conclusion, the creation of the DD dog breed about 120 years ago resulted in a popular breed with considerable genetic diversity without substructuring into lines or subpopulations. The trend of new inbreeding was declining, while ancestral inbreeding through ancestors who were autozygous at least once in previous generations was increasing.

Skeletal-dental features in 33 bull terrier dogs

Background

The Bull terrier breed has been reported in the veterinary literature to suffer frequent dental and skeletal malocclusions. In this retrospective case series, we report skeletal-dental anomalies in a group of 33 Bull terriers presented for a dental consultation.

Results

Out of 33 dogs examined, 24 cases had full mouth radiography or Cone-beam computed tomography performed. Eruption and development abnormalities observed were as follows: hypodontia in 54.1% (13/24), eruption changes in 29.2% (7/24), and tooth shape abnormalities in 33.3% (8/24). All dogs presented with some type of dental or skeletal malocclusion: neutroclusion was the most common (66.7% of the animals), followed by mandibular mesioclusion (18.8%), maxillo-mandibular asymmetry (9.4%), and mandibular distoclusion (6.3%). Dental abnormalities noted included rotation of mandibular and maxillary premolar teeth, distal displacement of the incisor teeth, lingual displacement of the mandibular canine teeth, and absence of mandibular premolar and molar teeth. Lingual displacement of mandibular canine teeth was associated with malocclusion causing trauma (odds ratio 7.1, 95% confidence interval [1.4 to 36.1], p = 0.024).

Conclusions

Malocclusions and tooth shape abnormalities were found to be the most frequent finding in this group of Bull terriers. Although these findings cannot be generalized to the global population further studies are needed to observe the true expression of these anomalies in the general breed population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-022-03164-0.

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.

Genetic Variability in Polish Lowland Sheepdogs Assessed by Pedigree and Genomic Data

Simple Summary

Dogs are an important part of society. The Polish Lowland Sheepdog (PON) is one of 353 of the world’s largest cynological organization listed dog breeds. Breeds with small population sizes, like the PONs, are often characterized by high inbreeding rates and thus an increased risk of congenital diseases. To examine the endangerment of the PONs, measures for genetic diversity and inbreeding were calculated for the German PON population. The study showed that the PONs had to be classified as a minimally endangered population according to threshold values specified by the European Association for Animal Production. However, the very recent trend showed a slight improvement.

Abstract

Genetic variability of Polish Lowland Sheepdog (PON) population was evaluated using both pedigree and genomic data. The analyzed pedigree encompassed 8628 PONs, including 153 individuals genotyped on the Illumina CanineHD BeadChip. Runs of homozygosity (ROH) were defined for homozygous stretches extending over 60 to 4300 kb. The inbreeding coefficients F_Ped based on pedigree data and F_ROH50 based on ROHs were at 0.18 and 0.31. The correlation between both was 0.41 but 0.52 when excluding animals with less than seven complete generations. The realized effective population size (Ne¯) was 22.2 with an increasing trend over years. Five PONs explained 79% of the genetic diversity of the reference population. The effective population size derived from linkage disequilibrium measured by r^² was 36. PANTHER analysis of genes in ROHs shared by ≥50% of the PONs revealed four highly over- or underrepresented biological processes. One among those is the 7.35 fold enriched “forelimb morphogenesis”. Candidate loci for hip dysplasia and patent ductus arteriosus were discovered in frequently shared ROHs. In conclusion, the inbreeding measures of the PONs were high and the genetic variability small compared to various dog breeds. Regarding Ne¯, PON population was minimally endangered according to the European Association for Animal Production.

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.

A Preliminary Study to Investigate the Genetic Background of Longevity Based on Whole-Genome Sequence Data of Two Methuselah Dogs

Aging is the largest risk factor in many diseases and mortality alike. As the elderly population is expected to increase at an accelerating rate in the future, these phenomena will pose a growing socio-economic burden on societies. To successfully cope with this challenge, a deeper understanding of aging is crucial. In many aspects, the companion dog is an increasingly popular model organism to study aging, with the promise of producing results that are more applicable to humans than the findings that come from the studies of classical model organisms. In this preliminary study we used the whole-genome sequence of two extremely old dogs - age: 22 and 27 years (or 90-135% more, than the average lifespan of dogs) - in order to make the first steps to understand the genetic background of extreme longevity in dogs. We identified more than ∼80 1000 novel SNPs in the two dogs (7500 of which overlapped between them) when compared to three publicly available canine SNP databases, which included SNP information from850 dogs. Most novel mutations (∼52000 SNPs) were identified at non-coding regions, while 4.6% of the remaining SNPs (n∼1600) were at exons, including 670 missense variants - 76 of which overlapped between the two animals - across 472 genes. Based on their gene ontologies, these genes were related - among others - to gene transcription/translation and its regulation, to immune response and the nervous system in general. We also detected 12 loss-of-function mutations, although their actual effect is unclear. Several genetic pathways were also identified, which pathways may be tempting candidates to be investigated in large sample sizes in order to confirm their relevance in extreme longevity in dogs (and possibly, in humans). We hypothesize a possible link between extreme longevity and the regulation of gene transcription/translation, which hypothesis should be further investigated in the future. This phenomenon could define an interesting direction for future research aiming to better understand longevity. The presented preliminary results highlight the utility of the companion dog in the study of the genetic background of longevity and aging.

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

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

Demographic assessment of the Dalmatian dog - effective population size, linkage disequilibrium and inbreeding coefficients

Background

The calculation of demographic measures is a useful tool for evaluating the genomic architecture of dog breeds and enables ranking dog breeds in terms of genetic diversity. To achieve this for the German Dalmatian dog population, 307 purebred animals of this breed were genotyped on the Illumina Canine high density BeadChip. The analysis of pedigree-based inbreeding was performed based on a pedigree with 25,761 dogs including the genotyped dogs.

Results

The effective population size derived from squared correlation coefficients between SNP alleles (r²) was 69. The maximum value of r² was 0.56, resulting in a 50% decay value of 0.28 at a marker distance of 37.5 kb. The effective population size calculated from pedigree data using individual increase in inbreeding over equivalent generations was 116. The pedigree inbreeding coefficient was 0.026. The genomic inbreeding coefficient based on the length of runs of homozygosity (ROH) was calculated for seven length categories of ROHs, and ranged from 0.08 to 0.28. The fixation coefficients F_{IS_PED} and F_{IS_GENO} were at 0.017 and 0.004. PANTHER statistical overrepresentation analysis of genes located in consensus ROHs revealed highly underrepresented biological processes in 50% of the investigated dogs. One of those is the 0.28 fold enriched “immune response”, which might be associated to the high prevalence of allergic dermatitis in the breed. Candidate genes for congenital sensorineural deafness (CCSD, a highly prevalent disease in the Dalmatian) were discovered in consensus ROHs.

Conclusions

The fast decay of r² and the moderate inbreeding coefficients indicate that the German Dalmatian dog population is rather diverse. Pedigree- and genomic-based inbreeding measures were highly correlated and therefore prove good reliability for the given population. Analyses of consensus ROHs with genes coding for deafness and other breed-defining traits, such as hyperuricosuria, indicate that those ROH became fixed in the Dalmatian population about 500 years ago. In case of the Dalmatian dog, a ROH of 40 SNPs length is enough to investigate signatures of selection (e.g. the ROH with the fixed hyperuricosuria mutation) as far back as the breed formation point approximately 500 years ago.

The dog 2.0: Lessons learned from the past

In recent years, concerns have been raised on the diversity, health and welfare of our (pedigree) dog population. Somewhat justified, the popular sire effect, population bottlenecks, the founder effect and inbreeding have left their marks on the dog as we know it. In order to improve the health and welfare of the canine population in general, individual breeding programs should adhere to the concept of ethical breeding (i.e. "the use of healthy animals true to their species in behaviour and looks, and when applicable, showing a sustainable performance") when population-specific breeding goals are defined. Even though every population has its own problems, the approach to get to possible solution(s) is similar. The starting point will always be the identification of which (un)desirable pheno- and genotypes are segregating and what their prevalence is, followed by an evaluation of the genetic diversity. Based on that information and, when applicable, additional criteria like breed standards, breeding goals can be defined. It is of critical importance that these goals are put forward with a long term vision in mind and with consensus from the stakeholders to ensure collaboration. Upon prioritization of the most important goals, when necessary with the help of specifically developed tools, the final step is choosing the most optimal combination of breeding strategies. This paper aims to provide a stepwise approach to identify and tackle population-specific problems encountered in breeding programs.

Population structure and genetic history of Tibetan Terriers

Background

Tibetan Terrier is a popular medium-sized companion dog breed. According to the history of the breed, the western population of Tibetan Terriers includes two lineages, Lamleh and Luneville. These two lineages derive from a small number of founder animals from the native Tibetan Terrier population, which were brought to Europe in the 1920s. For almost a century, the western population of Tibetan Terriers and the native population in Tibet were reproductively isolated. In this study, we analysed the structure of the western population of Tibetan Terriers, the original native population from Tibet and of different crosses between these two populations. We also examined the genetic relationships of Tibetan Terriers with other dog breeds, especially terriers and some Asian breeds, and the within-breed structure of both Tibetan Terrier populations.

Results

Our analyses were based on high-density single nucleotide polymorphism (SNP) array (Illumina HD Canine 170 K) and microsatellite (18 loci) genotypes of 64 Tibetan Terriers belonging to different populations and lineages. For the comparative analysis, we used 348 publicly available SNP array genotypes of dogs from other breeds. We found that the western population of Tibetan Terriers and the native Tibetan Terriers clustered together with other Asian dog breeds, whereas all other terrier breeds were grouped into a separate group. We were also able to differentiate the western Tibetan Terrier lineages (Lamleh and Luneville) from the native Tibetan Terrier population.

Conclusions

Our results reveal the relationships between the western and native populations of Tibetan Terriers and support the hypothesis that Tibetan Terrier belongs to the group of ancient dog breeds of Asian origin, which are close to the ancestors of the modern dog that were involved in the early domestication process. Thus, we were able to reject the initial hypothesis that Tibetan Terriers belong to the group of terrier breeds. The existence of this native population of Tibetan Terriers at its original location represents an exceptional and valuable genetic resource.

Tracing genetic resurrection of pointing dog breeds: Cesky Fousek as both survivor and rescuer

Cesky Fousek is considered to be one of the oldest pointing dog breeds in Europe and has been appreciated for its versatile working skills. Because it faced extinction in the past, the Cesky Fousek was restored from German Wirehaired and Shorthaired Pointers. Additionally, the breed was recently used in the USA with the initial intent of improvement of the Wirehaired Pointing Griffon (synonymous with Korthals Griffon) by the Bohemian Wirehaired Pointing Griffon Club of America. This study evaluates genetic diversity parameters of Cesky Fousek and compares them to the other continental pointing dogs that played a role in the formation of its gene pool. DNA from buccal swab and blood samples (n = 405) were analyzed using 18 microsatellite markers. Parameters of genetic polymorphism show that the Cesky Fousek breed has a comparable rate of variation as other hunting breeds despite the low population size and severe historical bottlenecks. Clustering analyses reveal a unique genetic status as a distinct pointing dog breed and the relatedness of the breeds is in good concordance with historical data. The present study demonstrates that despite historical admixture among lineages, separate pointing breeds constitute genetically differentiated units, mirroring unique breeding stocks and pedigree isolation among specific breed clubs, reflecting differences in breeding programs under each association.

Genetic heterogeneity and diversity of North American golden retrievers using a low density STR marker panel

Thirty-three autosomal short tandem repeat (STR) markers were used to evaluate genetic heterogeneity and diversity in 525 golden retrievers (GRs). This breed was selected because of its popularity and artificial selection for conformation vs. performance phenotypes. Seven additional STRs were used to evaluate the highly polymorphic dog leukocyte antigen (DLA) class I and class II regions. From 3 to 13 alleles were found at each of the 33 loci (mean 7) and the average effective alleles (Ne) was 3.34. The observed heterozygosity was 0.65 and the expected heterozygosity was 0.68. The resulting fixation index was 0.035 indicating that the population was randomly breeding. We found that modern GRs retain 46% of genomic diversity present in all canids and 21/175 (12%) and 20/90 (22%) of the known DLA class I and class II haplotypes, respectively. Selection for performance or conformation led to a narrowing of genomic and DLA diversity with conformation having a greater effect than performance. A comparison was made between coefficient of inbreeding (COI) determined from 10 or 12 generation pedigrees and DNA based internal relatedness values. A weak but significant correlation was observed between IR score and 10 or 12 generation COI (r = 0.38, p<0.0001 and r = 0.40, p<0.0001, respectively). IR values were higher in conformation than performance lines but only significant at p = 0.17. This was supported by 10 and 12 generation COI values that were significantly (p<0.0001) higher in conformation than performance lines. We demonstrate herein that a low density of STR markers can be utilized to study the genetic makeup of GRs.

Elucidating biogeographical patterns in Australian native canids using genome wide SNPs

Dingoes play a strong role in Australia's ecological framework as the apex predator but are under threat from hybridization and agricultural control programs. Government legislation lists the conservation of the dingo as an important aim, yet little is known about the biogeography of this enigmatic canine, making conservation difficult. Mitochondrial and Y chromosome DNA studies show evidence of population structure within the dingo. Here, we present the data from Illumina HD canine chip genotyping for 23 dingoes from five regional populations, and five New Guinea Singing Dogs to further explore patterns of biogeography using genome-wide data. Whole genome single nucleotide polymorphism (SNP) data supported the presence of three distinct dingo populations (or ESUs) subject to geographical subdivision: southeastern (SE), Fraser Island (FI) and northwestern (NW). These ESUs should be managed discretely. The FI dingoes are a known reservoir of pure, genetically distinct dingoes. Elevated inbreeding coefficients identified here suggest this population may be genetically compromised and in need of rescue; current lethal management strategies that do not consider genetic information should be suspended until further data can be gathered. D statistics identify evidence of historical admixture or ancestry sharing between southeastern dingoes and South East Asian village dogs. Conservation efforts on mainland Australia should focus on the SE dingo population that is under pressure from domestic dog hybridization and high levels of lethal control. Further data concerning the genetic health, demographics and prevalence of hybridization in the SE and FI dingo populations is urgently needed to develop evidence based conservation and management strategies.

The companion dog as a model for human aging and mortality

Around the world, human populations have experienced large increases in average lifespan over the last 150 years, and while individuals are living longer, they are spending more years of life with multiple chronic morbidities. Researchers have used numerous laboratory animal models to understand the biological and environmental factors that influence aging, morbidity, and longevity. However, the most commonly studied animal species, laboratory mice and rats, do not experience environmental conditions similar to those to which humans are exposed, nor do we often diagnose them with many of the naturally occurring pathologies seen in humans. Recently, the companion dog has been proposed as a powerful model to better understand the genetic and environmental determinants of morbidity and mortality in humans. However, it is not known to what extent the age‐related dynamics of morbidity, comorbidity, and mortality are shared between humans and dogs. Here, we present the first large‐scale comparison of human and canine patterns of age‐specific morbidity and mortality. We find that many chronic conditions that commonly occur in human populations (obesity, arthritis, hypothyroidism, and diabetes), and which are associated with comorbidities, are also associated with similarly high levels of comorbidity in companion dogs. We also find significant similarities in the effect of age on disease risk in humans and dogs, with neoplastic, congenital, and metabolic causes of death showing similar age trajectories between the two species. Overall, our study suggests that the companion dog may be an ideal translational model to study the many complex facets of human morbidity and mortality.

Assessing the relative importance of health and conformation traits in the cavalier king Charles spaniel

Background

The selection of a future breeding dog is a complicated task, in which disease characteristics and different traits have to be combined and weighed against one another. Truncation selection, that is the exclusion of affected animals, may be very inefficient when selecting on a large number of traits, and may result in a reduction of the genetic diversity in a population or breed. Selection could be facilitated by the use of a selection index that combines multiple traits or breeding values into one score. This however requires a consideration of their relative value according to their economic weight, which is difficult to express in monetary units for health traits. The use of a choice experiment to derive non-market values might be a solution to this problem. This is a pilot study to assess the potential use of choice experiments to ascertain the public preference and relative importance attached to health- and conformation traits in the selection of a Cavalier King Charles spaniel. The focus was on two prevalent disorders, mitral valve disease and syringomyelia, and on several important conformation traits such as muzzle length and eye shape. Based on available prior information, a Bayesian D-optimal design approach was used to develop a choice experiment and the resulting choice sets.

Results

Every participant (breeder or owner) in the choice experiment was presented with a total of 17 choice sets, in which at most four traits could vary to reduce the cognitive burden. A total of 114 respondents participated in the choice experiment and results showed that respondents (breeders/owners) current attitudes were directed towards health (syringomyelia and mitral valve disease), followed by eye shape and level of inbreeding.

Conclusions

This approach identifies the value breeders and owners attach to certain traits in the breeding objective. The resulting relative weights, represented as the logworths obtained from the choice experiment, could be an alternative to economic weights. They could be implemented as a weight when breeding values are available, but more study on this topic will be necessary. A challenge in this approach is to scale up the experiment with additional traits. Moreover, for other traits, the genetic parameters and correlations should be known first, in order to include them in the health selection index as well.

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

Objective

Braque Français, type Pyrénées is a French hunting-dog breed whose origin is traced back to old pointing gun-dogs used to assist hunters in finding and retrieving game. This breed is popular in France, but seldom seen elsewhere. Despite the ancient background, the literature on its genetic characterization is surprisingly scarce. A recent study looked into the demography and inbreeding using pedigree records, but there is yet no report on the use of molecular markers in this breed. The aim of this work was to genotype a population of Braque Français, type Pyrénées dogs with the high-density SNP array to study the genomic diversity of the breed.

Results

The average observed (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_O$$\end{document}HO) and expected (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_E$$\end{document}HE) heterozygosity were 0.371 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm \,0.142$$\end{document}±0.142) and 0.359 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm \,0.124$$\end{document}±0.124). Effective population size (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_e$$\end{document}Ne) was 27.5635 runs of homozygosity (ROH) were identified with average length of 2.16 MB. A ROH shared by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$75\%$$\end{document}75% of the dogs was detected at the beginning of chromosome 22. Inbreeding coefficients from marker genotypes were in the range \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{IS}=[-\,0.127,0.172]$$\end{document}FIS=[-0.127,0.172]. Inbreeding estimated from ROH (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{ROH}$$\end{document}FROH) had mean \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.112\,(\pm \,0.023$$\end{document}0.112(±0.023), with range [0.0526, 0.225]. These results show that the Braque Français, type Pyrénées breed is a relatively inbred population, but with still sufficient genetic variability for conservation and genetic improvement.

Electronic supplementary material

The online version of this article (10.1186/s13104-017-3112-9) contains supplementary material, which is available to authorized users.

Cross-Breeding Is Inevitable to Conserve the Highly Inbred Population of Puffin Hunter: The Norwegian Lundehund

The Norwegian Lundehund is a highly endangered native dog breed. Low fertility and high frequency predisposition to intestinal disorder imply inbreeding depression. We assessed the genetic diversity of the Lundehund population from pedigree data and evaluated the potential of optimal contribution selection and cross-breeding in the long-term management of the Lundehund population. The current Norwegian Lundehund population is highly inbred and has lost 38.8% of the genetic diversity in the base population. Effective population size estimates varied between 13 and 82 depending on the method used. Optimal contribution selection alone facilitates no improvement in the current situation in the Lundehund due to the extremely high relatedness of the whole population. Addition of (replacement with) 10 breeding candidates of foreign breed to 30 Lundehund breeders reduced the parental additive genetic relationship by 40-42% (48-53%). Immediate actions are needed to increase the genetic diversity in the current Lundehund population. The only option to secure the conservation of this rare breed is to introduce individuals from foreign breeds as breeding candidates.

Merging pedigree databases to describe and compare mating practices and gene flow between pedigree dogs in France, Sweden and the UK

Merging pedigree databases across countries may improve the ability of kennel organizations to monitor genetic variability and health-related issues of pedigree dogs. We used data provided by the Société Centrale Canine (France), Svenska Kennelklubben (Sweden) and the Kennel Club (UK) to study the feasibility of merging pedigree databases across countries and describe breeding practices and international gene flow within the following four breeds: Bullmastiff (BMA), English setter (ESE), Bernese mountain dog (BMD) and Labrador retriever (LBR). After merging the databases, genealogical parameters and founder contributions were calculated according to the birth period, breed and registration country of the dogs. Throughout the investigated period, mating between close relatives, measured as the proportion of inbred individuals (considering only two generations of pedigree), decreased or remained stable, with the exception of LBR in France. Gene flow between countries became more frequent, and the origins of populations within countries became more diverse over time. In conclusion, the potential to reduce inbreeding within purebred dog populations through exchanging breeding animals across countries was confirmed by an improved effective population size when merging populations from different countries.

Comparative Analysis of Genome Diversity in Bullmastiff Dogs

Management and preservation of genomic diversity in dog breeds is a major objective for maintaining health. The present study was undertaken to characterise genomic diversity in Bullmastiff dogs using both genealogical and molecular analysis. Genealogical analysis of diversity was conducted using a database consisting of 16,378 Bullmastiff pedigrees from year 1980 to 2013. Additionally, a total of 188 Bullmastiff dogs were genotyped using the 170,000 SNP Illumina CanineHD Beadchip. Genealogical parameters revealed a mean inbreeding coefficient of 0.047; 142 total founders (f); an effective number of founders (fe) of 79; an effective number of ancestors (fa) of 62; and an effective population size of the reference population of 41. Genetic diversity and the degree of genome-wide homogeneity within the breed were also investigated using molecular data. Multiple-locus heterozygosity (MLH) was equal to 0.206; runs of homozygosity (ROH) as proportion of the genome, averaged 16.44%; effective population size was 29.1, with an average inbreeding coefficient of 0.035, all estimated using SNP Data. Fine-scale population structure was analysed using NETVIEW, a population analysis pipeline. Visualisation of the high definition network captured relationships among individuals within and between subpopulations. Effects of unequal founder use, and ancestral inbreeding and selection, were evident. While current levels of Bullmastiff heterozygosity, inbreeding and homozygosity are not unusual, a relatively small effective population size indicates that a breeding strategy to reduce the inbreeding rate may be beneficial.

Why is the dog an ideal model for aging research?

With many caveats to the traditional vertebrate species pertaining to biogerontology investigations, it has been suggested that a most informative model is the one which: 1) examines closely related species, or various members of the same species with naturally occurring lifespan variation, 2) already has adequate medical procedures developed, 3) has a well annotated genome, 4) does not require artificial housing, and can live in its natural environment while being investigated, and 5) allows considerable information to be gathered within a relatively short period of time. The domestic dog unsurprisingly fits each criterion mentioned. The dog has already become a key model system in which to evaluate surgical techniques and novel medications because of the remarkable similarity between human and canine conditions, treatments, and response to therapy. The dog naturally serves as a disease model for study, obviating the need to construct artificial genetically modified examples of disease. Just as the dog offers a natural model for human conditions and diseases, simple observation leads to the conclusion that the canine aging phenotype also mimics that of the human. Genotype information, biochemical information pertaining to the GH/IGF-1 pathway, and some limited longitudinal investigations have begun the establishment of the domestic dog as a model of aging. Although we find that dogs indeed are a model to study aging and there are many independent pieces of canine aging data, there are many more "open" areas, ripe for investigation.

Genetic characterization of four native Italian shepherd dog breeds and analysis of their relationship to cosmopolitan dog breeds using microsatellite markers

Very little research into genetic diversity of Italian native dog breeds has been carried out so far. In this study we aimed to estimate and compare the genetic diversity of four native Italian shepherd dog breeds: the Maremma, Bergamasco, Lupino del Gigante and Oropa shepherds. Therefore, some cosmopolitan dog breeds, which have been widely raised in Italy for a long time past, have also been considered to check possible influence of these dog populations on the Italian autochthonous breeds considered here. A total of 212 individuals, belonging to 10 different dog breeds, were sampled and genotyped using 18 autosomal microsatellite loci. We analyzed the genetic diversity of these breeds, within breed diversity, breed relationship and population structure. The 10 breeds considered in this study were clearly genetically differentiated from each other, regardless of current population sizes and the onset of separate breeding history. The level of genetic diversity explained 20% of the total genetic variation. The level of H E found here is in agreement with that found by other studies. The native Italian breeds showed generally higher genetic diversity compared with the long established, well-defined cosmopolitan dog breeds. As the Border Collie seems closer to the Italian breeds than the other cosmopolitan shepherd dogs considered here, a possible utilization of this breed to improve working performance in Italian traditional working shepherd dogs cannot be ignored. The data and information found here can be utilized in the organization of conservation programs planned to reduce inbreeding and to minimize loss of genetic variability.

An insight into population structure and gene flow within pure-bred cats

Investigation of genetic structure on the basis of pedigree information requires indicators adapted to the specific context of the populations studied. On the basis of pedigree-based estimates of diversity, we analysed genetic diversity, mating practices and gene flow among eight cat populations raised in France, five of them being single breeds and three consisting of breed groups with varieties that may interbreed. When computed on the basis of coancestry rate, effective population sizes ranged from 127 to 1406, while the contribution of founders from other breeds ranged from 0.7 to 16.4%. In the five breeds, FIS ranged between 0.96 and 1.83%, with this result being related to mating practices such as close inbreeding (on average 5% of individuals being inbred within two generations). Within the three groups of varieties studied, FIT ranged from 1.59 to 3%, while FST¯ values were estimated between 0.04 and 0.91%, which was linked to various amounts of gene exchanges between subpopulations at the parental level. The results indicate that cat breeds constitute populations submitted to low selection intensity, contrasting with relatively high individual inbreeding level caused by close inbreeding practices.

Methods to estimate effective population size using pedigree data: Examples in dog, sheep, cattle and horse

Background

Effective population sizes of 140 populations (including 60 dog breeds, 40 sheep breeds, 20 cattle breeds and 20 horse breeds) were computed using pedigree information and six different computation methods. Simple demographical information (number of breeding males and females), variance of progeny size, or evolution of identity by descent probabilities based on coancestry or inbreeding were used as well as identity by descent rate between two successive generations or individual identity by descent rate.

Results

Depending on breed and method, effective population sizes ranged from 15 to 133 056, computation method and interaction between computation method and species showing a significant effect on effective population size (P < 0.0001). On average, methods based on number of breeding males and females and variance of progeny size produced larger values (4425 and 356, respectively), than those based on identity by descent probabilities (average values between 93 and 203). Since breeding practices and genetic substructure within dog breeds increased inbreeding, methods taking into account the evolution of inbreeding produced lower effective population sizes than those taking into account evolution of coancestry. The correlation level between the simplest method (number of breeding males and females, requiring no genealogical information) and the most sophisticated one ranged from 0.44 to 0.60 according to species.

Conclusions

When choosing a method to compute effective population size, particular attention should be paid to the species and the specific genetic structure of the population studied.

Marker panels for genealogy-based mapping, breed demographics, and inference-of-ancestry in the dog

Short tandem repeat polymorphisms (STRPs) are robust and informative markers for a range of genetic applications. STRPs are advantageous in experimental designs that derive power from sampling many individuals rather than many loci (e.g., pedigree-based studies, fine-scale mapping, and conservation genetics). STRPs have proven useful for vetting samples prior to costly high-density SNP analysis. Here we present validated STRPs (n = 1,012) spanning the canine genome (2.1 +/-1.4 Mb; 2.1 +/-2.1 cM). Standardized design, pre-multiplexing, M13-based dye-labeling, and selection for loci amenable to semi-automated allele-scoring minimize cost and facilitate efficient genotyping. The markers are leveraged from the canine linkage map, and thus are backed by genetic data useful for parametric multipoint analysis and assessment of empiric coverage. We demonstrate several applications with different marker subsets. The complete set provides a genome scan for linkage at ∼5 cM resolution. A subset of the markers measures molecular diversity between domestic and wild canid populations. Another subset reflects ancestry within breeds, uncovering hidden stratification and flagging genetic outliers prior to SNP genotyping. Thus, the markers described here add flexibility and cost effectiveness to several genetic applications in the dog that complement genome-wide SNP genotyping studies. Supplemental material is available for this article. Go to the publisher's online edition of Animal Biotechnology.

Low neutral genetic variability in a specialist puffin hunter: the Norwegian Lundehund

The genetic variability of 125 Norwegian Lundehund and 27 Nova Scotia Duck Tolling Retriever was analysed using a set of 26 microsatellite markers. In Lundehund, the average number of alleles per locus was 1.73, and average observed (H(O)) and expected (H(E)) heterozygosity were 0.07. In Toller, all measures of genetic diversity were much higher than in Lundehund and similar to studies on other dog breeds. The cluster analysis correctly assigned individuals to their respective breed. The low genetic variability in Lundehund was not surprising, given the two strong bottlenecks in the 1940s and the 1960s. The relatedness of Lundehund to other Nordic small spitzes should be investigated in the view of possible outcrossing.

Gene diversity, agroecological structure and introgression patterns among village chicken populations across North, West and Central Africa

Background

Chickens represent an important animal genetic resource for improving farmers’ income in Africa. The present study provides a comparative analysis of the genetic diversity of village chickens across a subset of African countries. Four hundred seventy-two chickens were sampled in 23 administrative provinces across Cameroon, Benin, Ghana, Côte d’Ivoire, and Morocco. Geographical coordinates were recorded to analyze the relationships between geographic distribution and genetic diversity. Molecular characterization was performed with a set of 22 microsatellite markers. Five commercial lines, broilers and layers, were also genotyped to investigate potential gene flow. A genetic diversity analysis was conducted both within and between populations.

Results

High heterozygosity levels, ranging from 0.51 to 0.67, were reported for all local populations, corresponding to the values usually found in scavenging populations worldwide. Allelic richness varied from 2.04 for a commercial line to 4.84 for one population from Côte d’Ivoire. Evidence of gene flow between commercial and local populations was observed in Morocco and in Cameroon, which could be related to long-term improvement programs with the distribution of crossbred chicks. The impact of such introgressions seemed rather limited, probably because of poor adaptation of exotic birds to village conditions, and because of the consumers’ preference for local chickens. No such gene flow was observed in Benin, Ghana, and Côte d’Ivoire, where improvement programs are also less developed. The clustering approach revealed an interesting similarity between local populations found in regions sharing high levels of precipitation, from Cameroon to Côte d’Ivoire. Restricting the study to Benin, Ghana, and Côte d’Ivoire, did not result in a typical breed structure but a south-west to north-east gradient was observed. Three genetically differentiated areas (P < 0.01) were identified, matching with Major Farming Systems (namely Tree Crop, Cereal-Root Crop, and Root Crop) described by the FAO.

Conclusions

Local chickens form a highly variable gene pool constituting a valuable resource for human populations. Climatic conditions, farming systems, and cultural practices may influence the genetic diversity of village chickens in Africa. A higher density of markers would be needed to identify more precisely the relative importance of these factors.

Genetic diversity, structure, and breed relationships in Iberian cattle

In Iberia there are 51 officially recognized cattle breeds of which 15 are found in Portugal and 38 in Spain. We present here a comprehensive analysis of the genetic diversity and structure of Iberian cattle. Forty of these breeds were genotyped with 19 highly polymorphic microsatellite markers. Asturiana de los Valles displayed the greatest allelic diversity and Mallorquina the least. Unbiased heterozygosity values ranged from 0.596 to 0.787. The network based on Reynolds distances was star-shaped with few pairs of interrelated breeds and a clear cluster of 4 breeds (Alistana/Arouquesa/Marinhoa/Mirandesa). The analysis of the genetic structure of Iberian cattle indicated that the most probable number of population clusters included in the study would be 36. Distance results were supported by the STRUCTURE software indicating a relatively recent origin or possible crossbreeding or both between pairs or small groups of breeds. Five clusters included 2 different breeds (Betizu/Pirenaica, Morucha/Avileña, Parda de Montaña/Bruna de los Pirineos, Barrosã/Cachena, and Toro de Lidia/Brava de Lide), 3 breeds (Berrenda en Negro, Negra Andaluza, and Mertolenga) were divided in 2 independent clusters each, and 2 breeds were considered admixed (Asturiana de los Valles and Berrenda en Colorado). Individual assignation to breeds was not possible in the 2 admixed breeds and the pair Parda de Montaña/Bruna de los Pirineos. The relationship between Iberian cattle reflects their geographical origin rather than their morphotypes. Exceptions to this geographic clustering are most probably a consequence of crossbreeding with foreign breeds. The relative genetic isolation within their geographical origin, the consequent genetic drift, the adaptation to specific environment and production systems, and the influence of African and European cattle have contributed to the current genetic status of Iberian cattle, which are grouped according to their geographical origin. The greater degree of admixture observed in some breeds should be taken into account before using molecular markers for genetic assignment of individuals to breeds.