Effective population size, extended linkage disequilibrium and signatures of selection in the rare dog breed lundehund

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.

A Variant-Centric Analysis of Allele Sharing in Dogs and Wolves

Canines are an important model system for genetics and evolution. Recent advances in sequencing technologies have enabled the creation of large databases of genetic variation in canines, but analyses of allele sharing among canine groups have been limited. We applied GeoVar, an approach originally developed to study the sharing of single nucleotide polymorphisms across human populations, to assess the sharing of genetic variation among groups of wolves, village dogs, and breed dogs. Our analysis shows that wolves differ from each other at an average of approximately 2.3 million sites while dogs from the same breed differ at nearly 1 million sites. We found that 22% of the variants are common across wolves, village dogs, and breed dogs, that ~16% of variable sites are common across breed dogs, and that nearly half of the differences between two dogs of different breeds are due to sites that are common in all clades. These analyses represent a succinct summary of allele sharing across canines and illustrate the effects of canine history on the apportionment of genetic variation.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

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.

Genetic Rescue of the Highly Inbred Norwegian Lundehund

Augmenting the genetic diversity of small, inbred populations by the introduction of new individuals is often termed “genetic rescue”. An example is the Norwegian Lundehund, a small spitz dog with inbreeding-related health problems that is being crossed with three Nordic breeds, including the Norwegian Buhund. Conservation breeding decisions for the (typically) small number of outcrossed individuals are vital for managing the rescue process, and we genotyped the Lundehund (n = 12), the Buhund (n = 12), their crosses (F1, n = 7) and first-generation backcrosses to the Lundehund (F2, n = 12) with >170,000 single nucleotide polymorphism loci to compare their levels of genetic diversity. We predicted that genome-wide diversity in F2 dogs would be higher than in the Lundehund but lower than in the F1 and the Buhund, and the heterozygosity values showed the expected patterns. We also found that runs of homozygosity, extended chromosomal regions of homozygous genotypes inherited from a common ancestor, were reduced in F2 individuals compared with Lundehund individuals. Our analyses demonstrate the benefits of outcrossing but indicate that some of the acquired genetic diversity is lost following immediate backcrossing. Additional breeding among F2 crosses could therefore merit from further consideration in genetic rescue management.

Linkage Disequilibrium and Effective Population Size of Buffalo Populations of Iran, Turkey, Pakistan, and Egypt Using a Medium Density SNP Array

Linkage disequilibrium (LD) across the genome provides information to identify the genes and variations related to quantitative traits in genome-wide association studies (GWAS) and for the implementation of genomic selection (GS). LD can also be used to evaluate genetic diversity and population structure and reveal genomic regions affected by selection. LD structure and Ne were assessed in a set of 83 water buffaloes, comprising Azeri (AZI), Khuzestani (KHU), and Mazandarani (MAZ) breeds from Iran, Kundi (KUN) and Nili-Ravi (NIL) from Pakistan, Anatolian (ANA) buffalo from Turkey, and buffalo from Egypt (EGY). The values of corrected r² (defined as the correlation between two loci) of adjacent SNPs for three pooled Iranian breeds (IRI), ANA, EGY, and two pooled Pakistani breeds (PAK) populations were 0.24, 0.28, 0.27, and 0.22, respectively. The corrected r² between SNPs decreased with increasing physical distance from 100 Kb to 1 Mb. The LD values for IRI, ANA, EGY, and PAK populations were 0.16, 0.23, 0.24, and 0.21 for less than 100Kb, respectively, which reduced rapidly to 0.018, 0.042, 0.059, and 0.024, for a distance of 1 Mb. In all the populations, the decay rate was low for distances greater than 2Mb, up to the longest studied distance (15 Mb). The r² values for adjacent SNPs in unrelated samples indicated that the Affymetrix Axiom 90 K SNP genomic array was suitable for GWAS and GS in these populations. The persistency of LD phase (PLDP) between populations was assessed, and results showed that PLPD values between the populations were more than 0.9 for distances of less than 100 Kb. The Ne in the recent generations has declined to the extent that breeding plans are urgently required to ensure that these buffalo populations are not at risk of being lost. We found that results are affected by sample size, which could be partially corrected for; however, additional data should be obtained to be confident of the results.

Unravelling selection signatures in a single dog breed suggests recent selection for morphological and behavioral traits

Strong selection has resulted in substantial morphological and behavioral diversity across modern dog breeds, which makes dogs interesting model animals to study the underlying genetic architecture of these traits. However, results from between-breed analyses may confound selection signatures for behavior and morphological features that were coselected during breed development. In this study, we assess population genetic differences in a unique resource of dogs of the same breed but with systematic behavioral selection in only one population. We exploit these different breeding backgrounds to identify signatures of recent selection. Selection signatures within populations were found on chromosomes 4 and 19, with the strongest signals in behavior-related genes. Regions showing strong signals of divergent selection were located on chromosomes 1, 24, and 32, and include candidate genes for both physical features and behavior. Some of the selection signatures appear to be driven by loci associated with coat color (Chr 24; ASIP) and length (Chr 32; FGF5), while others showed evidence of association with behavior. Our findings suggest that signatures of selection within dog breeds have been driven by selection for morphology and behavior. Furthermore, we demonstrate that combining selection scans with association analyses is effective for dissecting the traits under selection.

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.

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 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.

Wolf outside, dog inside? The genomic make-up of the Czechoslovakian Wolfdog

BACKGROUND

Genomic methods can provide extraordinary tools to explore the genetic background of wild species and domestic breeds, optimize breeding practices, monitor and limit the spread of recessive diseases, and discourage illegal crossings. In this study we analysed a panel of 170k Single Nucleotide Polymorphisms with a combination of multivariate, Bayesian and outlier gene approaches to examine the genome-wide diversity and inbreeding levels in a recent wolf x dog cross-breed, the Czechoslovakian Wolfdog, which is becoming increasingly popular across Europe.

RESULTS

Pairwise FST values, multivariate and assignment procedures indicated that the Czechoslovakian Wolfdog was significantly differentiated from all the other analysed breeds and also well-distinguished from both parental populations (Carpathian wolves and German Shepherds). Coherently with the low number of founders involved in the breed selection, the individual inbreeding levels calculated from homozygosity regions were relatively high and comparable with those derived from the pedigree data. In contrast, the coefficient of relatedness between individuals estimated from the pedigrees often underestimated the identity-by-descent scores determined using genetic profiles. The timing of the admixture and the effective population size trends estimated from the LD patterns reflected the documented history of the breed. Ancestry reconstruction methods identified more than 300 genes with excess of wolf ancestry compared to random expectations, mainly related to key morphological features, and more than 2000 genes with excess of dog ancestry, playing important roles in lipid metabolism, in the regulation of circadian rhythms, in learning and memory processes, and in sociability, such as the COMT gene, which has been described as a candidate gene for the latter trait in dogs.

CONCLUSIONS

In this study we successfully applied genome-wide procedures to reconstruct the history of the Czechoslovakian Wolfdog, assess individual wolf ancestry proportions and, thanks to the availability of a well-annotated reference genome, identify possible candidate genes for wolf-like and dog-like phenotypic traits typical of this breed, including commonly inherited disorders. Moreover, through the identification of ancestry-informative markers, these genomic approaches could provide tools for forensic applications to unmask illegal crossings with wolves and uncontrolled trades of recent and undeclared wolfdog hybrids.

Studies of modern Italian dog populations reveal multiple patterns for domestic breed evolution

Through thousands of years of breeding and strong human selection, the dog (Canis lupus familiaris) exists today within hundreds of closed populations throughout the world, each with defined phenotypes. A singular geographic region with broad diversity in dog breeds presents an interesting opportunity to observe potential mechanisms of breed formation. Italy claims 14 internationally recognized dog breeds, with numerous additional local varieties. To determine the relationship among Italian dog populations, we integrated genetic data from 263 dogs representing 23 closed dog populations from Italy, seven Apennine gray wolves, and an established dataset of 161 globally recognized dog breeds, applying multiple genetic methods to characterize the modes by which breeds are formed within a single geographic region. Our consideration of each of five genetic analyses reveals a series of development events that mirror historical modes of breed formation, but with variations unique to the codevelopment of early dog and human populations. Using 142,840 genome-wide SNPs and a dataset of 1,609 canines, representing 182 breeds and 16 wild canids, we identified breed development routes for the Italian breeds that included divergence from common populations for a specific purpose, admixture of regional stock with that from other regions, and isolated selection of local stock with specific attributes.

Genomic data illuminates demography, genetic structure and selection of a popular dog breed

BACKGROUND

Genomic methods have proved to be important tools in the analysis of genetic diversity across the range of species and can be used to reveal processes underlying both short- and long-term evolutionary change. This study applied genomic methods to investigate population structure and inbreeding in a common UK dog breed, the Labrador Retriever.

RESULTS

We found substantial within-breed genetic differentiation, which was associated with the role of the dog (i.e. working, pet, show) and also with coat colour (i.e. black, yellow, brown). There was little evidence of geographical differentiation. Highly differentiated genomic regions contained genes and markers associated with skull shape, suggesting that at least some of the differentiation is related to human-imposed selection on this trait. We also found that the total length of homozygous segments (runs of homozygosity, ROHs) was highly correlated with inbreeding coefficient.

CONCLUSIONS

This study demonstrates that high-density genomic data can be used to quantify genetic diversity and to decipher demographic and selection processes. Analysis of genetically differentiated regions in the UK Labrador Retriever population suggests the possibility of human-imposed selection on craniofacial characteristics. The high correlation between estimates of inbreeding from genomic and pedigree data for this breed demonstrates that genomic approaches can be used to quantify inbreeding levels in dogs, which will be particularly useful where pedigree information is missing.

Genetic rescue of an endangered domestic animal through outcrossing with closely related breeds: A case study of the Norwegian Lundehund

Genetic rescue, outcrossing with individuals from a related population, is used to augment genetic diversity in populations threatened by severe inbreeding and extinction. The endangered Norwegian Lundehund dog underwent at least two severe bottlenecks in the 1940s and 1960s that each left only five inbred dogs, and the approximately 1500 dogs remaining world-wide today appear to descend from only two individuals. The Lundehund has a high prevalence of a gastrointestinal disease, to which all remaining dogs may be predisposed. Outcrossing is currently performed with three Nordic Spitz breeds: Norwegian Buhund, Icelandic Sheepdog, and Norrbottenspets. Examination of single nucleotide polymorphism (SNP) genotypes based on 165K loci in 48 dogs from the four breeds revealed substantially lower genetic diversity for the Lundehund (H_E 0.035) than for other breeds (H_E 0.209–0.284). Analyses of genetic structure with > 15K linkage disequilibrium-pruned SNPs showed four distinct genetic clusters. Pairwise F_ST values between Lundehund and the candidate breeds were highest for Icelandic Sheepdog, followed by Buhund and Norrbottenspets. We assessed the presence of outlier loci among candidate breeds and examined flanking genome regions (1 megabase) for genes under possible selection to identify potential adaptive differences among breeds; outliers were observed in flanking regions of genes associated with key functions including the immune system, metabolism, cognition and physical development. We suggest crossbreeding with multiple breeds as the best strategy to increase genetic diversity for the Lundehund and to reduce the incidence of health problems. For this project, the three candidate breeds were first selected based on phenotypes and then subjected to genetic investigation. Because phenotypes are often paramount for domestic breed owners, such a strategy could provide a helpful approach for genetic rescue and restoration of other domestic populations at risk, by ensuring the involvement of owners, breeders and managers at the start of the project.

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.

Variant detection and runs of homozygosity in next generation sequencing data elucidate the genetic background of Lundehund syndrome

BACKGROUND

The Lundehund is a highly specialized breed characterized by a unique flexibility of the joints and polydactyly in all four limbs. The extremely small population size and high inbreeding has promoted a high frequency of diseased dogs affected by the Lundehund syndrome (LS), a severe gastro-enteropathic disease.

RESULTS

Comprehensive analysis of bead chip and whole-genome sequencing data for LS in the Lundehund resulted in a genome-wide association signal on CFA 34 and LS-specific runs of homozygosity (ROH) in this region. Filtering analysis for variants with predicted high or moderate effects revealed a missense mutation in LEPREL1 1.2 Mb proximal to the region of the genome-wide association, which was shown to be significantly associated with LS. LS-affected Lundehund harbored the mutant LEPREL1:g.139212C>G genotype A/A whereas all controls of other breeds showed the C/C wild type. In addition, ROH analysis for the Lundehund indicated a high enrichment of genes in potential signatures of selection affecting protein activation and immunoregulatory processes like NOD1 potentially involved in LS breed disposition.

CONCLUSIONS

Sequencing results for Lundehund specific traits reveal a potential causative mutation for LS in the neuropeptide operating gene LEPREL1 and suggests it as a precursor of the inflammatory process. Analyses of ROH regions give an insight into the genetic background of characteristic traits in the Lundehund that remain to be elucidated in the future.

Runs of homozygosity reveal signatures of positive selection for reproduction traits in breed and non-breed horses

Background

Modern horses represent heterogeneous populations specifically selected for appearance and performance. Genomic regions under high selective pressure show characteristic runs of homozygosity (ROH) which represent a low genetic diversity. This study aims at detecting the number and functional distribution of ROHs in different horse populations using next generation sequencing data.

Methods

Next generation sequencing was performed for two Sorraia, one Dülmen Horse, one Arabian, one Saxon-Thuringian Heavy Warmblood, one Thoroughbred and four Hanoverian. After quality control reads were mapped to the reference genome EquCab2.70. ROH detection was performed using PLINK, version 1.07 for a trimmed dataset with 11,325,777 SNPs and a mean read depth of 12. Stretches with homozygous genotypes of >40 kb as well as >400 kb were defined as ROHs. SNPs within consensus ROHs were tested for neutrality. Functional classification was done for genes annotated within ROHs using PANTHER gene list analysis and functional variants were tested for their distribution among breed or non-breed groups.

Results

ROH detection was performed using whole genome sequences of ten horses of six populations representing various breed types and non-breed horses. In total, an average number of 3492 ROHs were detected in windows of a minimum of 50 consecutive homozygous SNPs and an average number of 292 ROHs in windows of 500 consecutive homozygous SNPs. Functional analyses of private ROHs in each horse revealed a high frequency of genes affecting cellular, metabolic, developmental, immune system and reproduction processes. In non-breed horses, 198 ROHs in 50-SNP windows and seven ROHs in 500-SNP windows showed an enrichment of genes involved in reproduction, embryonic development, energy metabolism, muscle and cardiac development whereas all seven breed horses revealed only three common ROHs in 50-SNP windows harboring the fertility-related gene YES1. In the Hanoverian, a total of 18 private ROHs could be shown to be located in the region of genes potentially involved in neurologic control, signaling, glycogen balance and reproduction. Comparative analysis of homozygous stretches common in all ten horses displayed three ROHs which were all located in the region of KITLG, the ligand of KIT known to be involved in melanogenesis, haematopoiesis and gametogenesis.

Conclusions

The results of this study give a comprehensive insight into the frequency and number of ROHs in various horses and their potential influence on population diversity and selection pressures. Comparisons of breed and non-breed horses suggest a significant artificial as well as natural selection pressure on reproduction performance in all types of horse populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1977-3) contains supplementary material, which is available to authorized users.