Estimation of Effective Population Size in the Sapsaree: A Korean Native Dog (Canis familiaris)

Origin and population structure of native dog breeds in the Korean peninsula and East Asia

To study the ancestry and phylogenetic relationships of native Korean dog breeds to other Asian dog populations, we analyzed nucleotide variations in whole-genome sequences of 205 canid individuals. Sapsaree, Northern Chinese indigenous dog, and Tibetan Mastiff were largely related to West Eurasian ancestry. Jindo, Donggyeongi, Shiba, Southern Chinese indigenous (SCHI), Vietnamese indigenous dogs (VIET), and Indonesian indigenous dogs were related to Southeast and East Asian ancestry. Among East Asian dog breeds, Sapsaree presented the highest haplotype sharing with German Shepherds, indicating ancient admixture of European ancestry to modern East Asian dog breeds. SCHI showed greater haplotype sharing with New Guinea singing dogs, VIET, and Jindo than with other Asian breeds. The predicted divergence time of East Asian populations from their common ancestor was approximately 2,000 to 11,000 years ago. Our results expand understanding of the genetic history of dogs in the Korean peninsula to the Asian continent and Oceanic region.

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.

Whole genome resequencing of the Iranian native dogs and wolves to unravel variome during dog domestication

BACKGROUND

Advances in genome technology have simplified a new comprehension of the genetic and historical processes crucial to rapid phenotypic evolution under domestication. To get new insight into the genetic basis of the dog domestication process, we conducted whole-genome sequence analysis of three wolves and three dogs from Iran which covers the eastern part of the Fertile Crescent located in Southwest Asia where the independent domestication of most of the plants and animals has been documented and also high haplotype sharing between wolves and dog breeds has been reported.

RESULTS

Higher diversity was found within the wolf genome compared with the dog genome. A total number of 12.45 million SNPs were detected in all individuals (10.45 and 7.82 million SNPs were identified for all the studied wolves and dogs, respectively) and a total number of 3.49 million small Indels were detected in all individuals (3.11 and 2.24 million small Indels were identified for all the studied wolves and dogs, respectively). A total of 10,571 copy number variation regions (CNVRs) were detected across the 6 individual genomes, covering 154.65 Mb, or 6.41%, of the reference genome (canFam3.1). Further analysis showed that the distribution of deleterious variants in the dog genome is higher than the wolf genome. Also, genomic annotation results from intron and intergenic regions showed that the proportion of variations in the wolf genome is higher than that in the dog genome, while the proportion of the coding sequences and 3'-UTR in the dog genome is higher than that in the wolf genome. The genes related to the olfactory and immune systems were enriched in the set of the structural variants (SVs) identified in this work.

CONCLUSIONS

Our results showed more deleterious mutations and coding sequence variants in the domestic dog genome than those in wolf genome. By providing the first Iranian dog and wolf variome map, our findings contribute to understanding the genetic architecture of the dog domestication.

Genetic diversity and population structure of the Sapsaree, a native Korean dog breed

Background

The Sapsaree is a breed of dog (Canis familiaris) native to Korea, which became perilously close to extinction in the mid-1980s. However, with systematic genetic conservation and restoration efforts, this breed was rescued from extinction and population sizes have been gradually increasing over the past few decades. The aim of this study was to ascertain novel information about the genetic diversity, population structure, and demographic history of the Sapsaree breed using genome-wide single nucleotide polymorphism data. We characterized the genetic profile of the Sapsaree breed by comparison with seven foreign dog breeds with similar morphologies to estimate genetic differentiation within and among these breeds.

Results

The results suggest that Sapsarees have higher genetic variance compared with the other breeds analyzed. The majority of the Sapsarees in this study share a discrete genetic pattern, although some individuals were slightly different, possibly as a consequence of the recent restoration process. Concordant results from analyses of linkage disequilibrium, effective population size, genetic diversity, and population structural analyses illustrate a relationship among the Sapsaree and the Tibetan breeds Tibetan terrier and Lhasa Apso, and a small genetic introgression from European breeds. The effective population size of the Sapsaree has contracted dramatically over the past generations, and is currently insufficient to maintain long-term viability of the breed’s genetic diversity.

Conclusions

This study provides novel insights regarding the genetic diversity and population structure of the native Korean dog breed Sapsaree. Our results suggest the importance of a strategic and systematic approach to ensure the genetic diversity and the authenticity of the Sapsaree breed.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0757-5) contains supplementary material, which is available to authorized users.

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.

Genome-wide analysis of the diversity and ancestry of Korean dogs

There are various hypotheses on dog domestication based on archeological and genetic studies. Although many studies have been conducted on the origin of dogs, the existing literature about the ancestry, diversity, and population structure of Korean dogs is sparse. Therefore, this study is focused on the origin, diversity and population structure of Korean dogs. The study sample comprised four major categories, including non-dogs (coyotes and wolves), ancient, modern and Korean dogs. Selected samples were genotyped using an Illumina CanineHD array containing 173,662 single nucleotide polymorphisms. The genome-wide data were filtered using quality control parameters in PLINK 1.9. Only autosomal chromosomes were used for further analysis. The negative off-diagonal variance of the genetic relationship matrix analysis depicted, the variability of samples in each population. FIS (inbreeding rate within a population) values indicated, a low level of inbreeding within populations, and the patterns were in concordance with the results of Nei's genetic distance analysis. The lowest FST (inbreeding rate between populations) values among Korean and Chinese breeds, using a phylogenetic tree, multi-dimensional scaling, and a TreeMix likelihood tree showed Korean breeds are highly related to Chinese breeds. The Korean breeds possessed a unique and large diversity of admixtures compared with other breeds. The highest and lowest effective population sizes were observed in Korean Jindo Black (485) and Korean Donggyeong White (109), respectively. The historical effective population size of all Korean dogs showed declining trend from the past to present. It is important to take immediate action to protect the Korean dog population while conserving their diversity. Furthermore, this study suggests that Korean dogs have unique diversity and are one of the basal lineages of East Asian dogs, originating from China.

A method for detecting recent changes in contemporary effective population size from linkage disequilibrium at linked and unlinked loci

Estimation of contemporary effective population size (Ne) from linkage disequilibrium (LD) between unlinked pairs of genetic markers has become an important tool in the field of population and conservation genetics. If data pertaining to physical linkage or genomic position are available for genetic markers, estimates of recombination rate between loci can be combined with LD data to estimate contemporary Ne at various times in the past. We extend the well-known, LD-based method of estimating contemporary Ne to include linkage information and show via simulation that even relatively small, recent changes in Ne can be detected reliably with a modest number of single-nucleotide polymorphism (SNP) loci. We explore several issues important to interpretation of the results and quantify the bias in estimates of contemporary Ne associated with the assumption that all loci in a large SNP data set are unlinked. The approach is applied to an empirical data set of SNP genotypes from a population of a marine fish where a recent, temporary decline in Ne is known to have occurred.

Estimation of effective population size using single-nucleotide polymorphism (SNP) data in Jeju horse

This study was conducted to estimate the effective population size using SNPs data of 240 Jeju horses that had raced at the Jeju racing park. Of the total 61,746 genotyped autosomal SNPs, 17,320 (28.1%) SNPs (missing genotype rate of >10%, minor allele frequency of <0.05 and Hardy–Weinberg equilibrium test P-value of <10^–6) were excluded after quality control processes. SNPs on the X and Y chromosomes and genotyped individuals with missing genotype rate over 10% were also excluded, and finally, 44,426 (71.9%) SNPs were selected and used for the analysis. The measures of the LD, square of correlation coefficient (r²) between SNP pairs, were calculated for each allele and the effective population size was determined based on r² measures. The polymorphism information contents (PIC) and expected heterozygosity (HE) were 0.27 and 0.34, respectively. In LD, the most rapid decline was observed over the first 1 Mb. But r² decreased more slowly with increasing distance and was constant after 2 Mb of distance and the decline was almost linear with log-transformed distance. The average r² between adjacent SNP pairs ranged from 0.20 to 0.31 in each chromosome and whole average was 0.26, while the whole average r² between all SNP pairs was 0.02. We observed an initial pattern of decreasing N_e and estimated values were closer to 41 at 1 ~ 5 generations ago. The effective population size (41 heads) estimated in this study seems to be large considering Jeju horse’s population size (about 2,000 heads), but it should be interpreted with caution because of the technical limitations of the methods and sample size.

Development of Genetic Markers for Triploid Verification of the Pacific Oyster, Crassostrea gigas

The triploid Pacific oyster, which is produced by mating tetraploid and diploid oysters, is favored by the aquaculture industry because of its better flavor and firmer texture, particularly during the summer. However, tetraploid oyster production is not feasible in all oysters; the development of tetraploid oysters is ongoing in some oyster species. Thus, a method for ploidy verification is necessary for this endeavor, in addition to ploidy verification in aquaculture farms and in the natural environment. In this study, a method for ploidy verification of triploid and diploid oysters was developed using multiplex polymerase chain reaction (PCR) panels containing primers for molecular microsatellite markers. Two microsatellite multiplex PCR panels consisting of three markers each were developed using previously developed microsatellite markers that were optimized for performance. Both panels were able to verify the ploidy levels of 30 triploid oysters with 100% accuracy, illustrating the utility of microsatellite markers as a tool for verifying the ploidy of individual oysters.