# Heritability and complex segregation analysis of deafness in Jack Russell Terriers

- Thomas R Famula
^{1}Email author, - Edward J Cargill
^{2}and - George M Strain
^{3}

**3**:31

https://doi.org/10.1186/1746-6148-3-31

© Famula et al; licensee BioMed Central Ltd. 2007

**Received: **26 January 2007

**Accepted: **13 November 2007

**Published: **13 November 2007

## Abstract

### Background

The association between patterns of pigmentation and deafness in the dog has a long-documented history, with reports dating back over one hundred years. Long suspected of having a genetic basis, the search for loci with a pronounced influence in the expression of hearing loss in the dog has yet to be successful. No studies in the dog to date have found a possible influence of a specific colour locus associated with deafness. The present study is intended to evaluate the heritability of deafness in the Jack Russell Terrier (JRT), characterize the mode of inheritance, and evaluate the existence of a sex, coat colour, or coat texture influence on the expression of sensorineural deafness.

### Results

The estimation of heritability of deafness in the JRT was 0.22 when deafness was considered a binary (normal/deaf) trait and 0.31 when deafness was considered a three-category (normal/unilateral/bilateral deafness). The influence of coat colour in the incidence of JRT deafness was statistically significant, indicating that dogs with more white are more likely to be deaf. The influence of sex or coat texture was not statistically significant in the incidence of JRT deafness. Complex segregation analysis revealed a model of a single locus with a large effect on the binary measure of hearing loss is not supported.

### Conclusion

This is the first attempt, to our knowledge, to characterize a genetic component responsible for deafness in the JRT. The heritability of deafness in the JRT was found to be 0.22 and 0.31 considering deafness to be a two-category or three-category trait, respectively. There appears to be an influence of coat colour on the expression of deafness. In an attempt to characterize the mode of inheritance of deafness in the JRT, a model of a single locus with a large effect on hearing loss is not supported with this data. Further study is needed to determine if a single locus may be influencing deafness in the JRT. While the absence of a clear mode of inheritance complicates genetic dissection of deafness in the JRT, the assembling of this pedigree provides a tool for eventually defining the genetic bases of this disorder.

## Background

The association between patterns of pigmentation and deafness in the dog has a long-documented history, with reports dating back over one hundred years. Long suspected of having a genetic basis, the search for loci with a pronounced influence in the expression of hearing loss in the dog has yet to be successful. One study [1] has discussed various candidate loci based upon research on deafness in the dog, the human, and the mouse, but to date no specific loci have been shown to influence deafness in the dog. The interested reader is directed to a more comprehensive review of a molecular genetic approach to deafness in dogs, one outlining the application of comparative genomics [2].

Perhaps the most dramatic association between patterns of pigmentation and deafness in the dog can be found in the Dalmatian. In the Dalmatian, iris colour is positively correlated with deafness (dogs with at least one blue iris are more likely to be deaf) and the presence of a colour patch is negatively correlated with deafness (dogs with a colour patch are less likely to be deaf) [3–6]. No studies in the dog to date have found a possible influence of a specific colour locus (*i.e*., black or liver in the Dalmatian) associated with deafness [6].

There is a review of several breeds known to have a high risk for pigment-associated sensorineural deafness [6]. In addition to the Dalmatian, the review reported evidence showing white Bull Terriers are more likely to be deaf than Bull Terriers with colour in their coats. The Jack Russell Terrier (JRT) was also examined, however with records on only 56 dogs it could only be noted that 47 of the 56 dogs had normal hearing, with 4 dogs being unilaterally deaf and 5 bilaterally deaf. More observations would be needed in order to attempt to substantiate a genetic component to hearing loss in the JRT.

*i.e*., 90% or more white) as shown in Figure 1. The JRT standard also describes three different types of coat texture: smooth, rough, and broken (a combination of smooth and rough) [8]. Examples of smooth texture and broken texture are shown in Figure 1. The JRT can also have blue eyes, but a JRT with blue eyes is rare compared to the incidence in other breeds (

*e.g*., the Dalmatian). To the author's knowledge, no previous study has examined the influence of coat colour or coat texture on deafness in the JRT.

The present study is intended to evaluate the heritability of deafness in the JRT, characterize the mode of inheritance, and evaluate the existence of a sex, coat colour, or coat texture influence on the expression of sensorineural deafness.

## Results

Recorded phenotypes^{a} of Jack Russell Terriers.

Phenotype | Female | Male |
---|---|---|

Hearing/Smooth/White | 2 | 4 |

Hearing/Rough/White | 7 | 5 |

Hearing/Broken/White | 2 | 4 |

Hearing/Smooth/Tricolour | 7 | 5 |

Hearing/Rough/Tricolour | 6 | 8 |

Hearing/Broken/Tricolour | 5 | 4 |

Hearing/Smooth/Tan | 5 | 5 |

Hearing/Rough/Tan | 18 | 5 |

Hearing/Broken/Tan | 8 | 12 |

Hearing/Smooth/Black | 2 | 4 |

Hearing/Rough/Black | 0 | 1 |

Hearing/Broken/Black | 3 | 2 |

Unilateral/Smooth/White | 0 | 0 |

Unilateral/Rough/White | 2 | 1 |

Unilateral/Broken/White | 2 | 0 |

Unilateral/Smooth/Tricolour | 2 | 0 |

Unilateral/Rough/Tricolour | 1 | 1 |

Unilateral/Broken/Tricolour | 0 | 1 |

Unilateral/Smooth/Tan | 1 | 0 |

Unilateral/Rough/Tan | 1 | 2 |

Unilateral/Broken/Tan | 0 | 0 |

Unilateral/Smooth/Black | 0 | 0 |

Unilateral/Rough/Black | 0 | 0 |

Unilateral/Broken/Black | 0 | 0 |

Unilateral/Smooth/White | 0 | 0 |

Deaf/Smooth/White | 0 | 1 |

Deaf/Rough/White | 1 | 1 |

Deaf/Broken/White | 0 | 1 |

Deaf/Smooth/Tricolour | 0 | 0 |

Deaf/Rough/Tricolour | 0 | 2 |

Deaf/Broken/Tricolour | 0 | 0 |

Deaf/Smooth/Tan | 1 | 0 |

Deaf/Rough/Tan | 1 | 0 |

Deaf/Broken/Tan | 0 | 1 |

Deaf/Smooth/Black | 0 | 0 |

Deaf/Rough/Black | 0 | 0 |

Deaf/Broken/Black | 0 | 0 |

Estimation of heritability^{a}

Mean | Median | SD | Effective Sample Size | Convergence Score (p-value) | 95% HDR | Relative Risk | |
---|---|---|---|---|---|---|---|

| |||||||

Genetic Variance | 0.44 | 0.28 | 0.55 | 1664.7 | 0.90 (0.37) | 0.06,1.83 | |

Heritability | 0.26 | 0.22 | 0.16 | 1746.4 | 0.94 (0.35) | 0.06,0.65 | |

White – Tricolour | 0.69 | 0.68 | 0.42 | 8930.9 | 1.39 (0.16) | -0.12,1.51 | 1.51 |

White – Tan | 0.83 | 0.82 | 0.40 | 8384.9 | -0.74(0.46) | 0.06,1.62 | 1.58 |

Female-Male | -0.23 | -0.22 | 0.30 | 8619.0 | 0.13 (0.90) | -0.82,0.36 | 0.95 |

| |||||||

Genetic Variance | 0.80 | 0.45 | 1.29 | 1681.2 | 1.07 (0.28) | 0.12,3.67 | |

Heritability | 0.35 | 0.31 | 0.18 | 1786.5 | 1.10 (0.27) | 0.11,0.79 | |

White – Tricolour | 0.74 | 0.74 | 0.40 | 9868.5 | 1.19 (0.23) | -0.05,1.54 | 1.64 |

White – Tan | 0.84 | 0.84 | 0.39 | 8535.9 | 0.62 (0.53) | 0.06,1.61 | 1.63 |

Female-Male | -0.31 | -0.31 | 0.29 | 9496.6 | -0.18 (0.85) | -0.89,0.25 | 0.93 |

Also presented in Table 2 is the influence of coat colour in the incidence of JRT deafness. Consistent for each measure of hearing loss, the difference between mostly white dogs and those with a tan colour pattern indicates that dogs with more white are more likely to be deaf. This conclusion is justified by examination of the empirical 95% Highest Density Region (HDR) generated by the 10,000 Gibbs samples. Specifically, the comparison of white and tan coloured dogs found that 95% of the Gibbs sample estimates of this difference are between 0.06 and 1.61 (on the underlying unobservable continuous scale), an interval that does not overlap with 0. Given that we scored the trait 1 for deaf and 0 for unaffected, larger values of the parameter θ (from equation [e1]) imply a greater risk for deafness. Only the contrast for white and tan coated dogs had an HDR that did not overlap 0. Accordingly, only selected contrasts are presented in Table 2. As mentioned previously, there was only one black and white JRT with deafness (unilateral) in the pedigree, hence there was insufficient data to analyze this separate coat colour category and its possible influence on deafness. However, the contrasting results observed between the mostly white dogs and those with a tan colour pattern are sufficient to conclude a colour influence on the expression of deafness in the JRT. The contrast for coat texture had an HDR overlapping 0 (data not shown). Note, the significant contribution of the tan coat colour stands in contrast to the result found when the genetic contribution to deafness was not considered.

Estimation of heritability and genetic correlation^{a}

Mean | Median | SD | Effective Sample Size | Convergence Score (p-value) | 95% HDR | |
---|---|---|---|---|---|---|

| ||||||

Genetic Variance | 0.54 | 0.46 | 0.31 | 2291.9 | 0.98 (0.33) | 0.13,1.11 |

Heritability | 0.33 | 0.32 | 0.11 | 2209.9 | 1.51 (0.13) | 0.14,0.54 |

Female-Male | -0.23 | -0.22 | 0.29 | 9198.2 | 0.75 (0.45) | -0.79,0.33 |

| ||||||

Genetic Variance | 1.32 | 1.06 | 0.94 | 1206.3 | 1.68 (0.09) | 0.20,3.14 |

Heritability | 0.52 | 0.51 | 0.14 | 1380.9 | 1.52 (0.13) | 0.26,0.79 |

Female-Male | -0.26 | -0.25 | 0.35 | 8497.4 | -0.29 (0.77) | -0.95,0.40 |

| ||||||

r | 0.34 | 0.36 | 0.25 | 3258.1 | 1.55 (0.12) | -0.15,0.79 |

r | 0.22 | 0.23 | 0.12 | 4313.9 | 1.40 (0.16) | -0.10,0.47 |

| ||||||

Genetic Variance | 0.93 | 0.78 | 0.56 | 1408.8 | 0.28 (0.78) | 0.23,1.94 |

Heritability | 0.45 | 0.44 | 0.12 | 1749.8 | 0.23 (0.81) | 0.23,0.68 |

Female-Male | -0.33 | -0.33 | 0.31 | 9379.2 | 0.65 (0.52) | -0.93,0.29 |

| ||||||

Genetic Variance | 1.39 | 1.09 | 1.12 | 922.0 | -0.74 (0.46) | 0.18,3.37 |

Heritability | 0.53 | 0.52 | 0.14 | 1052.4 | -0.45 (0.65) | 0.26,0.80 |

Female-Male | -0.26 | -0.25 | 0.36 | 8427.6 | 1.17 (0.24) | -0.98,0.44 |

| ||||||

r | 0.34 | 0.36 | 0.24 | 2907.3 | 0.10 (0.92) | -0.12,0.76 |

r | 0.25 | 0.25 | 0.13 | 3036.1 | -0.01 (0.99) | -0.01,0.50 |

Mixed-inheritance model parameters^{a} for dichotomous deafness in Jack Russell Terriers

Polygenic Variance | Major Locus Variance | Additive Effect (a) | Dominance Deviation (d) | τ | τ | τ | Frequency (q) | |
---|---|---|---|---|---|---|---|---|

General Major Locus, Mendelian Transmission | ||||||||

Mean | 1.48 | 8.46 | 3.64 | -2.65 | 1.0 | 0.5 | 0.0 | 0.84 |

Median | 1.53 | 9.72 | 3.89 | -2.57 | - | - | - | 0.83 |

SD | 0.84 | 10.21 | 1.18 | 1.35 | - | - | - | 0.11 |

Eff Sample Size | 1417.3 | 1007.3 | 980.2 | 1013.6 | - | - | - | 1198.6 |

Conv Score (p) | 1.32 (0.19) | 1.55 (0.12) | 1.68 (0.09) | 1.71 (0.09) | - | - | - | 1.49 (0.14) |

HDR | 0.00 | 0.00 | 1.48 | -5.21 | - | - | - | 0.38 |

HDR 95% High | 3.22 | 33.14 | 7.77 | 0.06 | - | - | - | 1.00 |

Recessive Major Model, Mendelian Transmission | ||||||||

Mean | 2.21 | 13.19 | 3.68 | -3.68 | 1.0 | 0.5 | 0.0 | 0.74 |

Median | 2.33 | 12.32 | 3.65 | -3.65 | - | - | - | 0.75 |

SD | 0.60 | 6.03 | 0.84 | 0.84 | - | - | - | 0.41 |

Eff Sample Size | 1330.9 | 1121.0 | 1032.6 | 1032.6 | - | - | - | 1064.1 |

Conv Score (p) | -1.27 (0.20) | 1.51 (0.13) | 1.74 (0.08) | 1.74 (0.08) | - | - | - | 1.42 (0.16) |

HDR 95% Low | 0.39 | 0.00 | 2.15 | -6.23 | - | - | - | 0.44 |

HDR 95% High | 3.16 | 45.89 | 7.42 | -1.03 | - | - | - | 1.00 |

General Major Locus, Non-Mendelian Transmission | ||||||||

Mean | 1.46 | 6.53 | 4.13 | -3.01 | 0.63 | 0.52 | 0.18 | 0.94 |

Median | 1.57 | 6.34 | 4.36 | -3.29 | 0.70 | 0.51 | 0.01 | 0.99 |

SD | 0.88 | 8.18 | 1.62 | 1.57 | 0.05 | 0.05 | 0.13 | 0.03 |

Eff Sample Size | 1387.6 | 1109.2 | 1290.3 | 1088.1 | 1172 | 1043 | 941 | 1186.4 |

Conv Score (p) | 1.18 (0.24) | 1.40 (0.16) | 1.37 (0.17) | 1.19 (0.23) | 1.08 (.28) | 1.19 (.23) | 1.75 (.08) | 0.97 (0.33) |

HDR 95% Low | 0.00 | 0.00 | 0.53 | -7.11 | 0.52 | 0.26 | 0.00 | 0.75 |

HDR 95% High | 3.23 | 33.69 | 8.04 | -0.48 | 0.86 | 0.66 | 0.63 | 1.00 |

Not presented is an analysis of the trichotomous definition of deafness because such a phenotype cannot be evaluated legitimately with the iBay software; the software is limited to the evaluation of binary and normally distributed phenotypes. However, fitting the values of normal hearing, unilateral deafness and bilateral deafness as scores of 0, 1, 2, respectively, did support the conclusion of a segregating major locus (through examination of the 95% HDR).

## Discussion

This is the first attempt, to our knowledge, to characterize a genetic component responsible for deafness in the JRT. It is clear from the results presented in Table 2 that deafness in the JRT is hereditary and is influenced by genetic information passed from parent to offspring. As such, the heritability of deafness is of sufficient magnitude that attempts to select against it are potentially successful.

However, the heritability of deafness reported here for the JRT is lower than the estimates of the heritability of deafness in the Dalmatian [3–5], the only other breed with such estimates to the authors' knowledge. While the prevalence of deafness in these respective breeds is also not the same [6], as the JRT is less affected compared to Dalmatians, differing estimates of the heritability of deafness between the breeds could be suggestive of dissimilar mechanisms resulting in a similar phenotype. There are obvious issues (*e.g*., pedigree structure, pedigree size, *etc*) in comparisons between studies analyzing estimations of heritability within one breed. Those issues also apply in comparison of this study to any conducted utilizing another breed. However, since this is the only study to date analyzing deafness in the JRT, comparisons to studies in the Dalmatian are all that are available.

Unlike the Dalmatian, where the heritability of deafness has been estimated to be as high as 0.75 for a trichotomous trait [3], the highest heritability estimate in the JRT in this study of 0.31 for a trichotomous trait is not indicative by itself of a single major locus exerting a large effect. It has been demonstrated [11] that major loci tend to increase the heritability of a trait in a given population and a value greater than 0.70 is comparatively large for many polygenic traits. Because the estimate in the JRT is markedly below 0.70, other loci may be exerting an effect on any major locus responsible for deafness in the JRT. To further support this, Table 4 presents the results of the complex segregation analysis whereby the model of a single locus with a large effect on hearing loss is not supported.

Though we can conclude that deafness is heritable from the results in Table 2, the exact genetic mechanism that leads to expression of this disease cannot be stated with certainty based on the results in Table 4. A manual review of the pattern of inheritance did not support a model of a simple autosomal Mendelian locus. For example, the majority of the affected progeny were the result of matings of two unaffected parents, eliminating models of a single dominant deafness allele. Discarding a model of a single recessive autosomal allele is not possible with the pedigree, because there were no matings of two bilaterally deaf dogs; nor was there any mating of two unilaterally deaf dogs. Once again, the reader is reminded that our evaluation of a major locus may be influenced by possible ascertainment bias, our sample of dogs being built upon affected animals. Regrettably, there is no simple means by which this potential effect can be eliminated from the present analysis. What remains is the need for a cautious interpretation of the results of the several analyses.

Also presented in Table 2 is the influence of coat colour and coat texture in the incidence of JRT deafness. Consistent for each measure of hearing loss, there is a difference between mostly white dogs and those with a tan colour pattern. This indicates that mostly white dogs are more likely to be deaf, supporting that deafness in the JRT (as in other breeds such as the Dalmatian); there is a pigmentation association with deafness. The breed standard for the JRT mandates that the body be at least 51% white [7]. Unfortunately, no JRTs with blue eyes were available with this pedigree to evaluate a possible association between eye colour and deafness as has been observed in the Dalmatian. However, as the multiple trait analysis presented in Table 3 reveals, the association between coat colour and deafness appears to be less pronounced in the JRT than in the Dalmatian. Though the mean estimate of the genetic correlation (taken as the mean of the Gibbs sample) is strong and positive (e.g., a value of 0.34 for dichotomous deafness with the binary coat colour score), the 95% HDR suggests that estimate to be quite imprecise. Perhaps a larger sample of dogs, or a better means of quantifying coat colour would reveal a more precise relationship between hearing loss and pigmentation.

Although deafness in the JRT is clearly inherited, the evidence for the presence of a single major gene affecting the disorder is not persuasive with the data from this pedigree. In a review of complex segregation analysis [12], it was suggested to exercise caution in the interpretation of complex segregation analysis until several sets of data had confirmed or rejected the presence of a Mendelian locus. Further studies will be valuable in this context.

## Conclusion

The objectives of this study were to evaluate the heritability of deafness and the existence of a sex, coat colour, and/or coat texture influence on the expression of deafness in the JRT as well as characterize the mode of inheritance. The heritability of deafness in the JRT was found to be 0.22 and 0.31 considering deafness to be dichotomous and trichotomous, respectively. There appears to be an influence of coat colour on the expression of deafness as a difference was observed between mostly white JRTs and those with a tan colour pattern, indicating that white JRTs are more likely to be deaf. In an attempt to characterize the mode of inheritance of deafness in the JRT, a model of a single locus with a large effect on hearing loss is not supported with this data. Further study is needed to determine if a single locus may be influencing deafness in the JRT.

While the absence of a clear mode of inheritance complicates genetic dissection of deafness in the JRT, the assembling of this pedigree provides a tool for eventually defining the genetic bases of this disorder. Initially, one study [13] reported a heritability estimate of 0.32 in Californian Dalmatians, a value comparable to the estimate presented here in JRTs. However, a subsequent study [4] included more records of Californian Dalmatians in a larger data set and reported a higher heritability estimate of 0.76, a value comparable to the estimates presented by other studies involving Dalmatians [3, 5]. Inclusion of more individuals related to those in this pedigree may provide the opportunity for a more in-depth analysis of the heritability of deafness in JRTs. Also, collection of additional unilaterally or bilaterally deaf JRTs with black and white coat, as well as collection of normal and affected JRTs with blue eyes, will further elucidate the influence of pigmentation on the expression of deafness in this breed.

## Methods

### Phenotypic data

Phenotypes for hearing loss were measured using the brainstem auditory evoked response (BAER), permitting the discrimination between normal hearing dogs and bilaterally or unilaterally deaf dogs. Data were collected on 236 dogs from one large family, starting with a deaf proband and her immediate relatives, then extending broadly to available relatives. The majority of animals originated largely from the eastern United States (US), but subjects are included from across the country. BAER testing was typically done in puppies at five weeks of age, but several were tested as adults. Average test age is not known. Out of the 236 dogs, 201 have a known auditory status. The additional 35 dogs were included to help build appropriate pedigrees despite having an unknown auditory status. BAER measurements were previously performed at various test sites in the US at the owners' initiative. Copies of BAER test results were collected as confirmation of each animal's hearing status. BAER testing followed standard methods [6]. Figure 2 illustrates a subset of 51 dogs with known auditory status from this pedigree. Phenotypic information on sex, coat colour, coat texture, auditory status, and eye colour were collected by use of a standardized form distributed and returned by mail. In addition, pedigree information was recorded, along with the BAER results and phenotypic information of relatives where known. Colour classes were defined as white, tri-coloured, tan, or black, and coat texture was described as smooth, rough or broken [8]. There were no JRTs in the assembled pedigree with blue eyes. Also of note, only one dog with black and white coat colour was collected; a dog also affected by deafness (unilateral). The remaining JRTs with black and white coat colour in the pedigree have normal hearing or unknown auditory status.

An initial evaluation of the data was conducted with the R-program [10] without consideration of the pedigree or relationships among dogs in the data set. The purpose was to evaluate the potential impact of fixed effects such as sex or coat colour on deafness without regard to the genetic contributions to this disease. Analysis of deafness as a binary trait included models with effects for sex, coat colour and coat texture (along with all possible interactions) using a logit link function under the "glm" command of the R-program[10].

### Estimation of heritability

The BAER, used to determine the auditory function of each ear, provides for two possible deafness phenotypes in this pedigree. One phenotype is dichotomous, in which unilaterally deaf and bilaterally deaf dogs are classified as deaf (*i.e*., affected vs. unaffected). An alternative phenotype is trichotomous, with classes for normal hearing, unilateral deafness and bilateral deafness, representing ordered categories of increasing disease.

Most data sets utilized in the study of hereditary diseases are constructed around probands, making correction for ascertainment bias necessary; this set of data is no exception. In estimation of heritability, mixed linear models are capable of accommodating non-randomly sampled data [14]. Accordingly, the estimation of the heritability of deafness should not be biased by family selection, provided that the animals at the top of the pedigree (those animals with no parents identified) can be considered a random sample of JRTs. This is more assumption than assertion because it is not feasible to create or discount a process of selection against deafness or for sampling such animals disproportionately among those animals at the top that have no known auditory status.

We choose to estimate heritability through the use of threshold models [15], an approach typical for the analysis of binary and ordered categorical traits. For example, we consider the phenotype of deafness as a binary trait y_{ijkl} (where y_{ijkl} = 0 when unaffected; 1 when affected) for the *l*-th dog (*l* = 1, 2,...201) of the *i*-th sex (*i* = 1 for males; 2 for females) in the *j*-th coat texture class (*j* = 1 for smooth; 2 for rough; 3 for broken) and the *k*-th coat colour class (*k* = 1 for white; 2 for tricolour; 3 for tan; 4 for black). In threshold models, this binary phenotype is assumed to be related to an underlying, unobservable, normally distributed continuous variable, θ, through a set of three fixed thresholds, [γ_{0} = -∞; γ_{1} = 0; γ_{2} = ∞]; γ_{1} is set to zero for computational convenience, with no loss in generality or impact on subsequent analysis. Accordingly, the combination of continuous genetic and environmental terms thought to control the unobservable θ are translated into a categorical observation through comparison to the fixed thresholds (*i.e*., observe an unaffected dog when γ_{0} ≤ θ < γ_{1} or an affected dog when γ_{1} ≤ θ < γ_{2}).

We also consider deafness as a trichotomous trait, in which normal hearing dogs are scored as a zero, unilaterally deaf dogs scored as a one, and bilaterally deaf dogs are scored as a two. Such a characterization of the auditory phenotype requires only minor modification of the threshold model. Specifically we add a fourth fixed threshold [γ_{0} = -∞; γ_{1} = 0; γ_{2}; γ_{3} = ∞], where γ_{2} must be estimated from the data. Furthermore, normal hearing dogs would be observed when γ_{0} ≤ θ < γ_{1}, unilaterally deaf dogs would be observed when γ_{1} ≤ θ < γ_{2}, and bilaterally deaf dogs would be observed when γ_{2} ≤ θ < γ_{3}.

_{ijkl}= μ + sex

_{i}+ texture

_{j}+ colour

_{k}+ a

_{l}+ e

_{ijkl}

where θ_{ijkl} is an unobservable continuous variate for the *l*-th (*l* = 1, 2,..., 201) dog of the *i*-th sex in the *j*-th class of coat texture and the *k*-th coat colour class. The component μ is an unknown constant while sex_{i} is the contribution of the *i*-th sex to the expression of deafness. Coat texture_{j} and coat colour_{k} are similar contributions of these physical characteristics to the liability for deafness; a_{l} is the additive genetic contribution of the *l*-th animal and e_{ijkl} is an unknown residual. Both a_{l} and e_{ijkl} are assumed to be random effects with zero means and variances of σ_{a}^{2} (the additive genetic variance) and σ_{e}^{2} (the residual variance), respectively. The additive genetic effect for each animal accounts for the covariance in phenotypes of relatives and is assumed to be multivariately-normally distributed, with a covariance structure based upon the additive relationships among all 236 animals. Because the underlying scale is unobservable, the total variance is assumed to be σ_{P}^{2} = σ_{a}^{2} + σ_{e}^{2} where σ_{e}^{2} = 1.0, with no loss of generality [16–18]. The heritability of deafness, on the unobservable continuous scale, can be estimated as h^{2} = σ_{a}^{2}/(σ_{a}^{2} + σ_{e}^{2}).

To estimate the unknown fixed effects and unknown σ_{a}^{2} we used a mixed model Bayesian strategy [18]. An advantage of Bayesian methods is the ability to arrive at not only a point estimate of the unknown parameters (*e.g*., heritability), but also a distributional estimate. Though a more complete description of the statistical aspects of this analysis is available [17], briefly, the assumed prior densities for the fixed effects (sex, coat texture and colour effects) are uninformed, what Bayesian modellers refer to as a "flat" prior density. That is, we assume no prior knowledge of the magnitude of the fixed effects, allowing for the possibility that any value along the real line is a possible value. For the analysis of deafness as a binary observation there is no need to estimate the fixed thresholds. However, for the case of the trichotomous deafness, γ_{2} must be estimated. The assumed prior distribution for this parameter is the uniform with bounds established by γ_{1} and γ_{3}. As for the random contributions to θ, the additive genetic effects are assumed to be multivariately-normally distributed with a null mean and variance-covariance structure consisting of the numerator relationship matrix times the unknown additive genetic variance, σ_{a}^{2}. Similarly the random residuals are assumed to be independently normally distributed with null mean with variance σ_{e}^{2} = 1.0 (with no loss of generality since θ is an unobservable variate). Finally, given our Bayesian approach to this problem, we also must establish a prior density for the unknown variance σ_{a}^{2}. Specifically, we look to the inverted Wishart distribution where the expected prior mean for the additive genetic variance was started at 1.0, with a shape parameter of 3. The shape parameter reflects the degree of certainty we have in the choice of the prior mean for the additive genetic variance (the larger the value the more certain). A value of 3, speaking relatively, would be considered small, reflecting weak prior knowledge of the actual value for the additive genetic variance.

Estimation of the distribution of the unknown parameters employs a technique of numerical integration referred to as Gibbs sampling [19]. The algorithm is based on the iterative generation of a sequence of random variables from the known conditional distributions of the parameters, given the likelihood function of the data. Subsequent estimates of the parameters are found in the analysis of this sequence of random numbers, called the Gibbs sample. A more complete description of the Gibbs sampling process and its theoretical justification is available [18], as well as in the manual of the public domain software, MTGSAM [9], with which this analysis was performed.

In this study, the total length for the Gibbs sampling process was set to 300,000, with the first 50,000 samples discarded from any subsequent analysis (called the "burn-in"). The post-Gibbs analysis was implemented with the packages boa [20] and coda [21], both part of the R-program [10]. Convergence of the Gibbs sampling process was evaluated visually by sample plots and by a diagnostic test contrasting sample means from the first 10% of the sample with the last 50% of the sample [22]. Autocorrelations were calculated within the complete Gibbs sample to arrive at a suggested thinning rate. Gibbs sample statistics, including effective sample size, were calculated with a thinning rate of 25 (chosen based on computation of a maximum autocorrelation at lag 25 of 0.02 for all parameters), creating a final Gibbs sample of 10,000 sample observations (*i.e*., [300,000-50,000]/25 = 10,000). Highest density regions (HDR) were computed as described [23] with public domain software hdrcde, [24] a package within the R-program [10].

To further assess the relationship between deafness and coat colour, we also considered a multiple trait analysis of these two categorical phenotypes. Such an analysis permits estimation of the genetic correlation between hearing loss and coat colour on the underlying unobservable scale. The challenge of such an analysis is in finding the appropriate means to score coat colour. Because there is no way to rank coat colour differences on a presumed scale of colour, we chose a binary definition of coat colour, where dogs classified as white were given a score of 1.0 and all other colour combinations (e.g., tri colour, tan) were given a score of zero (0.0).

_{ijkl}= μ + sex

_{il}+ texture

_{jl}+ a

_{kl}+ e

_{ijkl}

where θ_{ijkl} is an unobservable continuous variate for the *k*-th (*k* = 1, 2,..., 201) dog of the *i*-th sex in the *j*-th class of coat texture for the *l*-th trait (*l* = 1 for deafness and *l* = 2 for binary coat colour). The terms for sex, coat texture and animal effect are as described in model [e1]. This extension of model [e1] requires estimation of two genetic variances, one for deafness and the second for coat colour, a genetic covariance between these two categorical phenotypes and a residual covariance between the two phenotypes. As described above for the univariate deafness model, unknown parameters were estimated from a Gibbs sample of 10,000 values.

### Complex segregation analysis

The possibility that deafness in JRTs is influenced by the action of a segregating locus of large effect can also be examined. This technique, called complex segregation analysis [25], is intended to integrate Mendelian transmission genetics and models of penetrance with the patterns of covariance expected in polygenic inheritance. A more complete description of complex segregation analysis is available [26].

An outline of the criteria that must be satisfied before acceptance of the single major locus model has been provided [27]. Adherence to these criteria reduces the number of false positives. Evaluation of the models necessary for complex segregation analysis was conducted with the package iBay (version 1.0) [28], an extension of the program MaGGIC [29] written to accommodate binary traits (but not trichotomous traits) in pedigrees with inbreeding. The iBay package [28] was recently used to evaluate the contribution of a major locus to osteochondral diseases in pigs [30], where a more complete outline of the Monte Carlo Markov chains approach is detailed.

The goal of this strategy was to simultaneously estimate the posterior density for a polygenic contribution to binary deafness disease along with the contributions of a putative Mendelian locus. Specifically, for this mixed-inheritance model, the strategy allowed the evaluation of a polygenic variance component, the additive and dominance contributions of a single locus (the parameters -a, d, and a for the putative major locus genotypes AA, AB, and BB, respectively) and the frequency of allele A of the putative major locus (defined as *"q"*). Given our scoring of binary phenotypes, where deaf (both unilateral and bilateral) is 1 and normal is scored as zero, the "B" allele represents the putative disease-enhancing allele. Note also that the iBay software models the unobservable scale of this threshold trait such that the residual variance is fixed at 1.0 (i.e., σ_{e}^{2} = 1).

Creation of the Gibbs sample requires several key assumptions about the behaviour of these unknown parameters. Though a variety of models can be considered, all are some variant of the following: sex, coat texture and coat colour as a fixed effects with a flat (i.e., uniform) prior densities, the polygenic variance component with a flat prior density, as well as flat prior densities for the additive, dominance, and allele frequency parameters. A Gibbs sample of 9,000 was generated, beginning with the creation of 350,000 total samples, a "burn-in" of 50,000 and a sampling rate of every 100-th Gibbs value. This process was repeated two additional times, to create three replicate chains. As outlined above, the post-Gibbs analysis was implemented with the packages boa [20] and coda [21], both part of the R-program [10]. Convergence of the Gibbs sampling process was evaluated visually by sample plots and by contrasting sample means from the first 10% of the sample with the last 50% of the sample [22]. Also as outlined above, from the 9,000 Gibbs samples, the mean, standard deviation, median and the upper and lower limits of a 95% HDR was computed for each of the unknown parameters with hdrcde [24].

## Declarations

### Acknowledgements

This work was supported by a grant from the Jack Russell Terrier Research Foundation to GMS. The authors wish to thank Jack Russell Terrier owners for their enthusiastic participation. The authors are also grateful to Mr. And Mrs. Don Wilson for providing photographs of the Jack Russell Terriers shown in Figure 1. Some of the results reported here were obtained by using the program package iBay which was made available by Dr. L. Janss and we would like to thank him for this generous contribution.

## Authors’ Affiliations

## References

- Rak SG, Drogemuller C, Leeb T, Quignon P, Andre C, Scott A, Breen M, Distl O: Chromosomal assignment of 20 candidate genes for canine congenital sensorineural deafness by FISH and RH mapping. Cytogenet Genome Res. 2003, 101: 130-135. 10.1159/000074168.PubMedView ArticleGoogle Scholar
- Rak SG, Distl O: Congenital sensorineural deafness in dogs: A molecular genetic approach toward unravelling the responsible genes. Vet J. 2005, 169: 188-196. 10.1016/j.tvjl.2004.01.015.PubMedView ArticleGoogle Scholar
- Cargill EJ, Famula TR, Strain GM, Murphy KE: Heritability and segregation analysis of deafness in U.S. Dalmatians. Genetics. 2004, 166: 1385-1393. 10.1534/genetics.166.3.1385.PubMedPubMed CentralView ArticleGoogle Scholar
- Famula TR, Oberbauer AM, Williams DC: Gender effects in hearing loss in Dalmatians. Prev Vet Med. 2001, 48: 15-24. 10.1016/S0167-5877(00)00183-5.PubMedView ArticleGoogle Scholar
- Juraschko K, Meyers-Lindenberg A, Nolte I, Distl O: Analysis of systemic effects on congenital sensorineural deafness in German Dalmatian dogs. Vet J. 2003, 166: 164-169. 10.1016/S1090-0233(02)00256-3.PubMedView ArticleGoogle Scholar
- Strain GM: Deafness prevalence and pigmentation and gender associations in dog breeds at risk. Vet J. 2004, 167: 23-32. 10.1016/S1090-0233(03)00104-7.PubMedView ArticleGoogle Scholar
- Jack Russell Terrier Club of America. [http://www.terrier.com/]
- Plummer DB: The Complete Jack Russell Terrier. 1998, McCarthy & Bassett, UKGoogle Scholar
- Van Tassell CP, Van Vleck LD: A Manual for use of MTGSAM: A set of Fortran programs to apply Gibbs sampling to animal models for variance component estimation. 1995, U. S. Department of Agriculture, Agricultural Research ServiceGoogle Scholar
- Team RDC: R: A language and environment for statistical computing. Vienna, Austria:2005.Google Scholar
- Morton NE, MacLean CJ: Analysis of family resemblance. III Complex segregation of quantitative traits. Am J Hum Genet. 1974, 26: 489-503.PubMedPubMed CentralGoogle Scholar
- Jarvik GP: Complex segregation analyses: uses and limitations. Am J Hum Genet. 1998, 63: 942-946. 10.1086/302075.PubMedPubMed CentralView ArticleGoogle Scholar
- Famula TR, Oberbauer AM, Sousa CA: Complex segregation analysis of deafness in Dalmatians. Am J Vet Res. 2000, 61: 550-553. 10.2460/ajvr.2000.61.550.PubMedView ArticleGoogle Scholar
- Henderson CR: Applications of linear models in animal breeding. Guelph, Ontario: University of Guelph: 1984.Google Scholar
- Falconer DS, Mackay TFC: Introduction to quantitative genetics Fourth Edition. Essex, England: Longman Group LTD:1996.Google Scholar
- Gianola D, Foulley JL: Sire evaluation for ordered categorical data with a threshold model. Genet Selec Evol. 1983, 15: 201-224. 10.1051/gse:19830203.View ArticleGoogle Scholar
- Harville DA, Mee RW: A mixed model procedure for analyzing ordered categorical data. Biometrics. 1984, 40: 393-408. 10.2307/2531393.View ArticleGoogle Scholar
- Sorensen DA, Andersen S, Gianola D, Korsgaard I: Bayesian inference in threshold models using Gibbs sampling. Genet Sel Evol. 1995, 27: 229-249. 10.1051/gse:19950303.PubMed CentralView ArticleGoogle Scholar
- Geman S, Geman D: Stochastic relaxation, Gibbs distributions and Bayesian restoration of images. IEEE Trans Patt Anal Mach Intell. 1984, 6: 721-741.View ArticleGoogle Scholar
- Smith BJ: boa: Bayesian Output Analysis Program (BOA) for MCMC. R package version 1.6.6. 2007Google Scholar
- Plummer M, Best N, Cowles K, Vines K: coda: Output analysis and diagnostics for MCMC. R package version 11-2. 2007Google Scholar
- Geweke J: Evaluating the accuracy of sampling-based approaches to calculating posterior moments. Bayesian Statistics 4. Edited by: Bernardo JM, Berger JO, Dawid AP, Smith AFM. Oxford University Press, Oxford:1992.Google Scholar
- Hyndman R: Computing and graphing highest density regions. The Amer Stat. 1996, 50: 120-126. 10.2307/2684423.Google Scholar
- Hyndman R: hdrcde: Highest density regions and conditional density estimation. R package version 2.06. 2007Google Scholar
- Bonney GE: Regressive logistic models for familial disease and other binary traits. Biometrics. 1986, 42: 611-625. 10.2307/2531211.PubMedView ArticleGoogle Scholar
- Lynch M, Walsh B: Genetics and Analysis of Quantitative Traits. Sinauer Associates, Sunderland, MA: 1998.Google Scholar
- Elston RC, Nasmboodrii KK, Glueck CJ, Fallat R, Tsang R, Leuba V: Studies of the genetic transmission of hypercholesterolemia and hypertriglyceridemia in a 195 member kindred. Ann Hum Genet. 1975, 39: 67-87. 10.1111/j.1469-1809.1975.tb00109.x.PubMedView ArticleGoogle Scholar
- Janss (2006) iBay Reference. [http://www.lucjanss.com]
- Janss LLG: Maggic: a package of subroutines for genetic analyses with Gibbs sampling. Proceedings of the Sixth World Congress on Genetics Applied to Livestock Production. Armidale, Australia; 1998:459-460.Google Scholar
- Kadarmideen HN, Janss LLG: Evidence of a major gene from Bayesian segregation analyses of liability to osteochondral diseases in pigs. Genetics. 2005, 171: 1195-1206. 10.1534/genetics.105.040956.PubMedPubMed CentralView ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.