The extent and distribution of linkage disequilibrium in a multi-hierarchic outbred canine pedigree

Heritability and genome-wide association study of vaccine-induced immune response in Beagles: A pilot study

Both genetic and non-genetic factors contribute to individual variation in the immune response to vaccination. Understanding how genetic background influences variation in both magnitude and persistence of vaccine-induced immunity is vital for improving vaccine development and identifying possible causes of vaccine failure. Dogs provide a relevant biomedical model for investigating mammalian vaccine genetics; canine breed structure and long linkage disequilibrium simplify genetic studies in this species compared to humans. The objective of this study was to estimate the heritability of the antibody response to vaccination against viral and bacterial pathogens, and to identify genes driving variation of the immune response to vaccination in Beagles. Sixty puppies were immunized following a standard vaccination schedule with an attenuated combination vaccine containing antigens for canine adenovirus type 2, canine distemper virus, canine parainfluenza virus, canine parvovirus, and four strains of Leptospira bacteria. Serum antibody measurements for each viral and bacterial component were measured at multiple time points. Heritability estimations and GWAS were conducted using SNP genotypes at 279,902 markers together with serum antibody titer phenotypes. The heritability estimates were: (1) to Leptospira antigens, ranging from 0.178 to 0.628; and (2) to viral antigens, ranging from 0.199 to 0.588. There was not a significant difference between overall heritability of vaccine-induced immune response to Leptospira antigens compared to viral antigens. Genetic architecture indicates that SNPs of low to high effect contribute to immune response to vaccination. GWAS identified two genetic markers associated with vaccine-induced immune response phenotypes. Collectively, these findings indicate that genetic regulation of the immune response to vaccination is antigen-specific and influenced by multiple genes of small effect.

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

Effective population size (N_e) is an important measure to understand population structure and genetic variability in animal species. The objective of this study was to estimate N_e in Sapsaree dogs using the information of rate of inbreeding and genomic data that were obtained from pedigree and the Illumina CanineSNP20 (20K) and CanineHD (170K) beadchips, respectively. Three SNP panels, i.e. Sap134 (20K), Sap60 (170K), and Sap183 (the combined panel from the 20K and 170K), were used to genotype 134, 60, and 183 animal samples, respectively. The N_e estimates based on inbreeding rate ranged from 16 to 51 about five to 13 generations ago. With the use of SNP genotypes, two methods were applied for N_e estimation, i.e. pair-wise r² values using a simple expectation of distance and r² values under a non-linear regression with respective distances assuming a finite population size. The average pair-wise N_e estimates across generations using the pairs of SNPs that were located within 5 Mb in the Sap134, Sap60, and Sap183 panels, were 1,486, 1,025 and 1,293, respectively. Under the non-linear regression method, the average N_e estimates were 1,601, 528, and 1,129 for the respective panels. Also, the point estimates of past N_e at 5, 20, and 50 generations ago ranged between 64 to 75, 245 to 286, and 573 to 646, respectively, indicating a significant N_e reduction in the last several generations. These results suggest a strong necessity for minimizing inbreeding through the application of genomic selection or other breeding strategies to increase N_e, so as to maintain genetic variation and to avoid future bottlenecks in the Sapsaree population.

Differential genetic regulation of canine hip dysplasia and osteoarthritis

BACKGROUND

Canine hip dysplasia (HD) is a common polygenic trait characterized by hip malformation that results in osteoarthritis (OA). The condition in dogs is very similar to developmental dysplasia of the human hip which also leads to OA.

METHODOLOGY/PRINCIPAL FINDINGS

A total of 721 dogs, including both an association and linkage population, were genotyped. The association population included 8 pure breeds (Labrador retriever, Greyhounds, German Shepherd, Newfoundland, Golden retriever, Rottweiler, Border Collie and Bernese Mountain Dog). The linkage population included Labrador retrievers, Greyhounds, and their crosses. Of these, 366 dogs were genotyped at ∼22,000 single nucleotide polymorphism (SNP) loci and a targeted screen across 8 chromosomes with ∼3,300 SNPs was performed on 551 dogs (196 dogs were common to both sets). A mixed linear model approach was used to perform an association study on this combined association and linkage population. The study identified 4 susceptibility SNPs associated with HD and 2 SNPs associated with hip OA.

CONCLUSION/SIGNIFICANCE

The identified SNPs included those near known genes (PTPRD, PARD3B, and COL15A1) reported to be associated with, or expressed in, OA in humans. This suggested that the canine model could provide a unique opportunity to identify genes underlying natural HD and hip OA, which are common and debilitating conditions in both dogs and humans.

Extent of genome-wide linkage disequilibrium in Australian Holstein-Friesian cattle based on a high-density SNP panel

BACKGROUND

The extent of linkage disequilibrium (LD) within a population determines the number of markers that will be required for successful association mapping and marker-assisted selection. Most studies on LD in cattle reported to date are based on microsatellite markers or small numbers of single nucleotide polymorphisms (SNPs) covering one or only a few chromosomes. This is the first comprehensive study on the extent of LD in cattle by analyzing data on 1,546 Holstein-Friesian bulls genotyped for 15,036 SNP markers covering all regions of all autosomes. Furthermore, most studies in cattle have used relatively small sample sizes and, consequently, may have had biased estimates of measures commonly used to describe LD. We examine minimum sample sizes required to estimate LD without bias and loss in accuracy. Finally, relatively little information is available on comparative LD structures including other mammalian species such as human and mouse, and we compare LD structure in cattle with public-domain data from both human and mouse.

RESULTS

We computed three LD estimates, D', Dvol and r2, for 1,566,890 syntenic SNP pairs and a sample of 365,400 non-syntenic pairs. Mean D' is 0.189 among syntenic SNPs, and 0.105 among non-syntenic SNPs; mean r2 is 0.024 among syntenic SNPs and 0.0032 among non-syntenic SNPs. All three measures of LD for syntenic pairs decline with distance; the decline is much steeper for r2 than for D' and Dvol. The value of D' and Dvol are quite similar. Significant LD in cattle extends to 40 kb (when estimated as r2) and 8.2 Mb (when estimated as D'). The mean values for LD at large physical distances are close to those for non-syntenic SNPs. Minor allelic frequency threshold affects the distribution and extent of LD. For unbiased and accurate estimates of LD across marker intervals spanning < 1 kb to > 50 Mb, minimum sample sizes of 400 (for D') and 75 (for r2) are required. The bias due to small samples sizes increases with inter-marker interval. LD in cattle is much less extensive than in a mouse population created from crossing inbred lines, and more extensive than in humans.

CONCLUSION

For association mapping in Holstein-Friesian cattle, for a given design, at least one SNP is required for each 40 kb, giving a total requirement of at least 75,000 SNPs for a low power whole-genome scan (median r2 > 0.19) and up to 300,000 markers at 10 kb intervals for a high power genome scan (median r2 > 0.62). For estimation of LD by D' and Dvol with sufficient precision, a sample size of at least 400 is required, whereas for r2 a minimum sample of 75 is adequate.

Admixture and patterns of linkage disequilibrium in a free-living vertebrate population

Linkage disequilibrium (LD), a measure of nonrandom association of alleles at different loci, is of great interest to evolutionary geneticists as it can be used to help identify loci that explain phenotypic variation. Surveys of the extent of LD across genomes have been carried out in a number of systems, most notably humans and model organisms. However, studies of natural populations of vertebrates have rarely been performed. Here, we describe an investigation of LD in a free-living island population of red deer Cervus elaphus. Relatively high levels of LD extended several tens of centimorgans, and significant LD was frequently detected between unlinked markers. The magnitude of LD varied depending on how the population was sampled. It also varied across different chromosomes, and was shown to be a function of sample size, intermarker distance and marker heterozygosity. A recent admixture event in the population led to an ephemeral increase in LD. Association mapping may be possible in this population, although a high 'baseline' level of LD could lead to false positive associations between marker loci and a trait of interest.

A random model for mapping imprinted quantitative trait loci in a structured pedigree: an implication for mapping canine hip dysplasia

Genetic imprinting may have played a more notable role in shaping embryonic development of plants, animals, and humans than previously appreciated. Quantitative trait loci that are imprinted (iQTL) exert monoallelic effects, depending on the parent of origin, which is an exception to the laws of Mendelian genetics. In this article, we present a modified random effect-based mapping model to use in a genome-wide scan for the distribution of iQTL that contribute to genetic variance for a complex trait in a structured pedigree. This model, implemented with the maximum likelihood method, capitalizes on a network of relatedness for maternally and paternally derived alleles through identical-by-descent sharing, thus allowing for the discrimination of the genetic variances due to alleles derived from maternal and paternal parents. The model was employed to map iQTL responsible for canine hip dysplasia in a multihierarchical canine pedigree, founded with seven greyhounds and six Labrador retrievers. Of eight significant QTL detected, three, located on CFA1, CFA8, and CF28, were found to trigger significant parent-of-origin effects on the age of femoral capital ossification measured at the left and right hips of a canine. The detected iQTL provide important candidate regions for fine-mapping of imprinted genes and for studying their structure and function in the control of complex traits.

Finding cardiovascular disease genes in the dog

Recent advances in canine genomics are changing the landscape of veterinary biology, and by default, veterinary medicine. No longer are clinicians locked into traditional methods of diagnoses and therapy. Rather, major advances in canine genetics and genomics from the past five years are now changing the way the veterinarian of the 21st century practices medicine. First, the availability of a dense genome map gives canine genetics a much-needed foothold in comparative medicine, allowing advances made in human and mouse genetics to be applied to companion animals. Second, the recently released 7.5x whole genome sequence of the dog is facilitating the identification of hereditary disease genes. Finally, development of genetic tools for rapid screening of families and populations at risk for inherited disease means that the cost of identifying and testing for disease loci will significantly decrease in coming years. Out of these advances will come major changes in companion animal diagnostics and therapy. Clinicians will be able to offer their clients genetic testing and counseling for a myriad of disorders. In this review we summarize recent findings in canine genomics and discuss their application to the study of canine cardiac health.

Canine genomics and genetics: running with the pack

The domestication of the dog from its wolf ancestors is perhaps the most complex genetic experiment in history, and certainly the most extensive. Beginning with the wolf, man has created dog breeds that are hunters or herders, big or small, lean or squat, and independent or loyal. Most breeds were established in the 1800s by dog fanciers, using a small number of founders that featured traits of particular interest. Popular sire effects, population bottlenecks, and strict breeding programs designed to expand populations with desirable traits led to the development of what are now closed breeding populations, with limited phenotypic and genetic heterogeneity, but which are ideal for genetic dissection of complex traits. In this review, we first discuss the advances in mapping and sequencing that accelerated the field in recent years. We then highlight findings of interest related to disease gene mapping and population structure. Finally, we summarize novel results on the genetics of morphologic variation.

Modeling extent and distribution of zygotic disequilibrium: implications for a multigenerational canine pedigree

Unlike gametic linkage disequilibrium defined for a random-mating population, zygotic disequilibrium describes the nonrandom association between different loci in a nonequilibrium population that deviates from Hardy-Weinberg equilibrium. Zygotic disequilibrium specifies five different types of disequilibria simultaneously that are (1) Hardy-Weinberg disequilibria at each locus, (2) gametic disequilibrium (including two alleles in the same gamete, each from a different locus), (3) nongametic disequilibrium (including two alleles in different gametes, each from a different locus), (4) trigenic disequilibrium (including a zygote at one locus and an allele at the other), and (5) quadrigenic disequilibrium (including two zygotes each from a different locus). However, because of the uncertainty on the phase of the double heterozygote, gametic and nongametic disequilibria need to be combined into a composite digenic disequilibrium and further define a composite quadrigenic disequilibrium together with the quadrigenic disequilibrium. To investigate the extent and distribution of zygotic disequilibrium across the canine genome, a total of 148 dogs were genotyped at 247 microsatellite markers located on 39 pairs of chromosomes for an outbred multigenerational pedigree, initiated with a limited number of unrelated founders. A major portion of zygotic disequilibrium was contributed by the composite digenic and quadrigenic disequilibrium whose values and numbers of significant marker pairs are both greater than those of trigenic disequilibrium. All types of disequilibrium are extensive in the canine genome, although their values tend to decrease with extended map distances, but with a greater slope for trigenic disequilibrium than for the other types of disequilibrium. Considerable variation in the pattern of disequilibrium reduction was observed among different chromosomes. The results from this study provide scientific guidance about the determination of the number of markers used for whole-genome association studies.

Genome-wide linkage disequilibrium in two Japanese beef cattle breeds

There is little knowledge about the degree of linkage disequilibrium (LD) in beef cattle. This study aims to perform a genome-wide search for LD in Japanese Black and Japanese Brown beef cattle and to compare the level of LD between these two breeds. Parameter D' (the LD coefficient) was used as a measure of LD, and LD was tested for significance of allelic associations between syntenic and between non-syntenic marker pairs. Effects of breed, chromosome, genetic map distance and their interactions with D' were tested based on least squares analyses. Both breeds showed high levels of LD, which ranged over several tens of cM and declined as the marker distance increased for syntenic marker pairs. A rapid decline of the D' value was observed between markers that were spaced 5 and 20 cM apart. LD was significant in most cases for marker pairs <40 cM apart but was not significant between non-syntenic loci. The pattern of LD found in these two breeds was similar to that previously published for dairy cattle. The D' value between breeds was not significantly different (P > 0.05), but the interaction between breed and chromosome was highly significant (P < 0.001). Genetic selection seems to have caused the heterogeneity of the D' values among chromosomes within breed. These results indicate that LD mapping is a useful tool for fine-mapping quantitative trait loci of economically important traits in Japanese beef cattle.

Genome sequence, comparative analysis and haplotype structure of the domestic dog

Here we report a high-quality draft genome sequence of the domestic dog (Canis familiaris), together with a dense map of single nucleotide polymorphisms (SNPs) across breeds. The dog is of particular interest because it provides important evolutionary information and because existing breeds show great phenotypic diversity for morphological, physiological and behavioural traits. We use sequence comparison with the primate and rodent lineages to shed light on the structure and evolution of genomes and genes. Notably, the majority of the most highly conserved non-coding sequences in mammalian genomes are clustered near a small subset of genes with important roles in development. Analysis of SNPs reveals long-range haplotypes across the entire dog genome, and defines the nature of genetic diversity within and across breeds. The current SNP map now makes it possible for genome-wide association studies to identify genes responsible for diseases and traits, with important consequences for human and companion animal health.

Uveodermatologic (VKH-like) syndrome in American Akita dogs is associated with an increased frequency of DQA1*00201

The Akita breed of dog is affected by a number of distinct immune-mediated diseases, including thyroiditis, sebaceous adenitis, pemphigus foliaceus, uveitis, polyarthritis, myasthenia gravis, and uveodermatologic (UV) syndrome. UV syndrome is manifested by progressive uveitis and depigmenting dermatitis that closely resembles the human Vogt - Koyanagi - Harada syndrome. This study examined the allelic diversity of the three DLA class II loci (DRB1, DQA1, and DQB1) in the American Akita dog, and the relationship of specific DLA class II alleles to the UV. Low allelic variation was demonstrated within genes of DLA class II. American Akita dogs possessed six of the reported 16 DQA1 alleles, but only eight of 61 reported alleles in DRB1 and nine of 47 reported alleles in DQB1. Almost one-half of American Akita dogs were homozygous for a single allele at DQA1 and up to a quarter at DRB1 and DQB1. DLA-DQA1*00201 was associated with a significantly higher relative risk (RR = 15.3) or odds ratio (OR = 15.99) for UV syndrome than other DLA class II alleles. No significant association was noted with haplotypes of DRB1, DQB1, and DQA1 alleles; DRB1*03201-DQA1*00201 trended toward significance. This study confirmed loss of DLA genetic diversity in the American Akita dog in common with other pure breeds of dog and suggested a role for certain DLA class II gene alleles in the pathogenesis of UV.

Quantitative trait loci for hip dysplasia in a crossbreed canine pedigree

Canine hip dysplasia is a common developmental inherited trait characterized by hip laxity, subluxation or incongruity of the femoral head and acetabulum in affected hips. The inheritance pattern is complex and the mutations contributing to trait expression are unknown. In the study reported here, 240 microsatellite markers distributed in 38 autosomes and the X chromosome were genotyped on 152 dogs from three generations of a crossbred pedigree based on trait-free Greyhound and dysplastic Labrador Retriever founders. Interval mapping was undertaken to map the QTL underlying the quantitative dysplastic traits of maximum passive hip laxity (the distraction index), the dorsolateral subluxation score, and the Norberg angle. Permutation testing was used to derive the chromosome-wide level of significance at p<0.05 for each QTL. Chromosomes 4, 9, 10, 11 (p<0.01), 16, 20, 22, 25, 29 (p<0.01), 30, 35, and 37 harbor putative QTL for one or more traits. Successful detection of QTL was due to the cross-breed pedigree, multiple-trait measurements, control of environmental background, and marked advancement in canine mapping tools.

Extensive and breed-specific linkage disequilibrium in Canis familiaris

The 156 breeds of registered dogs in the United States offer a unique opportunity to map genes important in disease susceptibility, morphology, and behavior. Linkage disequilibrium (LD) is of current interest for its application in whole genome association mapping, since the extent of LD determines the feasibility of such studies. We have measured LD at five genomic intervals, each 5 Mb in length and composed of five clusters of sequence variants spaced 800 kb-1.6 Mb apart. These intervals are located on canine chromosomes 1, 2, 3, 34, and 37, and none is under obvious selective pressure. Approximately 20 unrelated dogs were assayed from each of five breeds: Akita, Bernese Mountain Dog, Golden Retriever, Labrador Retriever, and Pekingese. At each genomic interval, SNPs and indels were discovered and typed by resequencing. Strikingly, LD in canines is much more extensive than in humans: D' falls to 0.5 at 400-700 kb in Golden Retriever and Labrador Retriever, 2.4 Mb in Akita, and 3-3.2 Mb in Bernese Mountain Dog and Pekingese. LD in dog breeds is up to 100x more extensive than in humans, suggesting that a correspondingly smaller number of markers will be required for association mapping studies in dogs compared to humans. We also report low haplotype diversity within regions of high LD, with 80% of chromosomes in a breed carrying two to four haplotypes, as well as a high degree of haplotype sharing among breeds.