Candidate genes for idiopathic epilepsy in four dog breeds

Background Idiopathic epilepsy (IE) is a naturally occurring and significant seizure disorder affecting all dog breeds. Because dog breeds are genetically isolated populations, it is possible that IE is attributable to common founders and is genetically homogenous within breeds. In humans, a number of mutations, the majority of which are genes encoding ion channels, neurotransmitters, or their regulatory subunits, have been discovered to cause rare, specific types of IE. It was hypothesized that there are simple genetic bases for IE in some purebred dog breeds, specifically in Vizslas, English Springer Spaniels (ESS), Greater Swiss Mountain Dogs (GSMD), and Beagles, and that the gene(s) responsible may, in some cases, be the same as those already discovered in humans. Results Candidate genes known to be involved in human epilepsy, along with selected additional genes in the same gene families that are involved in murine epilepsy or are expressed in neural tissue, were examined in populations of affected and unaffected dogs. Microsatellite markers in close proximity to each candidate gene were genotyped and subjected to two-point linkage in Vizslas, and association analysis in ESS, GSMD and Beagles. Conclusions Most of these candidate genes were not significantly associated with IE in these four dog breeds, while a few genes remained inconclusive. Other genes not included in this study may still be causing monogenic IE in these breeds or, like many cases of human IE, the disease in dogs may be likewise polygenic.


Background
The prevalence of epilepsy in humans is reported to be in the range of 4 -10/1000 in most study settings [1], with idiopathic epilepsy (IE) representing 15 -20% of these cases [2]. It is now generally accepted that IE in humans is due to an underlying genetic origin [2], although causative mutations have been discovered in only a small subset of IEs, mostly in isolated populations. These identified epilepsy mutations are, for the most part, Mendelian or monogenic IEs [3][4][5][6], and are often termed "channelopathies" due to their occurrence in ion channel genes. In a recent review, 16 of 21 susceptibility genes for human epilepsy were ion channels or neurotransmitter receptors [7].
Canine epilepsy is a naturally occurring, spontaneous condition. Canine seizures exhibit a remarkable resemblance to human seizures [8] and the usefulness of naturally occurring canine epilepsy as a translational model to explore potential treatments for human epilepsy was recently proposed by Leppik, et al. [9]. The prevalence of canine epilepsy is estimated to be between 0.5% and 5.7% and it is the most common chronic neurological disorder in dogs [10]. A diagnosis of IE in the canine indicates recurrent seizures for which no cause can be identified and implies a genetic predisposition [11]. While the molecular basis of IE in canines is entirely unknown, a hereditary basis for IE has been suggested in a number of dog breeds including Beagles [12], British Alsatians [13], Keeshonds [14], Labrador retrievers [15], Golden retrievers [16,17], Bernese mountain dogs [18], Belgian Tervurens [19,20], Boxers [21], Shetland sheepdogs [22], Vizslas [23], English Springer Spaniels [24], Irish wolfhounds [25], and Standard Poodles [26]. One inherited progressive seizure disorder that has been identified in dogs is the autosomal recessive Lafora disease in Miniature Wirehaired Dachshunds [27]. This progressive myoclonic epilepsy disorder results from a mutation in the EPM2B (also known as NHLRC1-NHL repeat containing 1) gene.
Modes of inheritance predicted from examination of pedigrees of dogs with IE show evidence for a gene of major effect or autosomal recessive inheritance, both of which could occur due to a founder effect. Examples include the English Springer Spaniel (partially penetrant autosomal recessive or polygenic inheritance) [24], the Vizsla (autosomal recessive or polygenic inheritance) [23], and the Golden Retriever [16] and the Bernese Mountain Dog [18] (both autosomal recessive but controlled by more than one gene). There is enough variation in mode of inheritance and clinical characteristics in dogs to presume that the underlying genetic basis is unlikely to be identical between breeds.
Linebreeding and inbreeding have become commonplace to meet the very specific dog breed standards of the hundreds of dog breeds in existence today. These practices, while creating more uniform offspring, can also result in deleterious diseases due to the concentration of recessive mutations passed from the founders. Almost half of the nearly four hundred genetic disorders that have been described in dogs are recognized exclusively in one or a few breeds [28], presumably due to a strong founder effect.
The aim of the present study was to take advantage of the likely founder effect occurring in purebred dogs with IE in order to investigate an underlying genetic basis for IE in four dog breeds: the Vizsla, the English Springer Spaniel (ESS), the Greater Swiss Mountain Dog (GSMD) and the Beagle. A candidate gene approach was utilized; genes known to be involved in inherited human epilepsy, along with two genes involved in mouse epilepsy models and additional related gene family members, were selected. Most candidate genes tested were ion channels. Microsatellite markers within or very close to these genes were used to generate genotypes which were then subjected to either two-point linkage analysis or association chi-square studies. Using four different breeds, which may each have different forms of IE, increased the chances of potentially finding an associated gene in at least one of the breeds.

Results
Call rates in each sample set for the microsatellite markers that were ultimately used in each analysis were approximately 90% or higher. Heterozygosity was recorded for each marker as a measure of that marker's extent of utility for this study. Average heterozygosity ± standard error of all markers analyzed for Beagles was 0.574 ± 0.024, for Vizslas was 0.569 ± 0.026, for ESS was 0.425 ± 0.030, and for GSMD was 0.341 ± 0.030. As expected, the average heterozygosity varied by breed and markers that were informative in one breed were not necessarily so in another breed. Out of all tested microsatellites, there were greater numbers of low heterozygosity (cut-off of 0.3) markers in the GSMD (27) and ESS (20) compared to the Vizslas (9) and Beagles (6), resulting in the lower average heterozygosities for the GSMD and ESS. For some genes, despite testing several microsatellites, results remain uninformative. For example, GABRD remains inconclusive in the ESS, despite testing three different microsatellites.
From the list of 52 candidate genes, 20 are known to be associated with inherited human epilepsy (indicated by a bolded entry with * in the tables) and two are associated with inherited mouse epilepsy models (indicated by a bolded entry with # in the tables). All 20 of the inherited human epilepsy-associated genes and both mouse model genes were analyzed at least once in all four of the test breeds, except for CHRNA2 in the Vizslas. However, not all of the other 30 genes were examined in every breed.
In the two-point linkage study, the two large Vizsla families were divided into ten smaller families to decrease complexity and inbreeding loops. The power of this ten family pedigree to detect significant linkage was demonstrated with simulated linkage analysis. The highest average LOD score for simulated linkage was 3.63, obtained by using a simulated marker with four alleles, no heterogeneity, and a recombination frequency of 0.001. The maximum attainable LOD score in these simulations was 8.46, obtained by using a simulated marker with three alleles, 5% heterogeneity, and a recombination frequency of 0.001. All of the simulation combinations obtained average LOD scores above 2.0 (suggestive linkage), demonstrating that the pedigrees are sufficient to detect linkage for a monogenic mutation. Results for the Vizslas are shown in Table 1. The majority of the candidate gene loci have very negative LOD scores and 43 genes were successfully excluded from being linked with IE in this breed. All of the human epilepsy-associated genes, with the exception of CHRNA2, and both of the mouse model genes were excluded. Four genes (CACNA1I, CACNB3, CHRNA1, and SCN11A) had low heterozygosity or monomorphism for the markers tested, therefore, these genes cannot be excluded. In addition, two genes (CACNB1 and CACNB2) had LOD values falling between -2.0 and 3.2, a grey zone of being neither linked nor unlinked; thus, conclusions cannot be drawn concerning possible linkage of these genes to IE in Vizslas. The most interesting of these was CACNB1, which had a LOD of 0.820 and a heterozygosity of 0.548.
Results for the English Springer Spaniels are shown in Table 2. Pearson chi-square p-values are reported for this nonparametric association study. The majority of the loci (n = 36) had highly insignificant p-values, decreasing the likelihood that these genes are associated with IE in this breed; eighteen of the inherited human epilepsy genes were in this group, as were both of the mouse epilepsy genes. Nine genes, two of which are human epilepsyassociated genes (GABRD and CHRNA2), had markers with heterozygosities falling below the 0.3 cut-off, prohibiting the drawing of conclusions about these genes. One  Table 3. As with the ESS, Pearson chi-square p-values are reported for this association study. Again, the majority of the loci (n = 35) had highly insignificant p-values and are unlikely to be associated with IE in the GSMD; this included 17 of the 20 inherited human epilepsy-associated genes and both of the mouse model genes. This breed had the lowest overall average heterozygosity and results for eleven genes, three of which are human epilepsy-associated (GABRD, KCNA1, and ME2), could not be determined due to low heterozygosity. The first marker analyzed for KCND2 resulted in a p-value of 0.036, prompting the examination of two additional markers. One of these was uninformative, but the other had a heterozygosity of 0.729 and a p-value of 0.155, rendering this gene not significant and unlikely to be associated to IE in GSMD. Additionally, CHRNB2, a gene that is associated with human epilepsy, resulted in a p-value of 0.055. While this p-value is potentially suggestive of association, when correcting for multiple testing with the Bonferroni correction, it does not approach significance.
Results for the Beagles are shown in Table 4. In this breed association study, highly insignificant p-values for 43 candidate genes indicate that these genes are likely not associated with IE in Beagles. All twenty of the inherited human epilepsy-associated genes and both of the mouse model genes had insignificant p-values in this breed. Only three genes (CACNA1D, CACNG3, and KCNQ5) were inconclusive due to low heterozygosity in this breed. CHRNA1 gave a p-value of 0.05 on its initial marker, leading to the development of primers for additional microsatellites. Two of these, CHRNA1's third set and fourth set of primers, gave p-values of 0.072 and 0.056, respectively, but CHRNA1's second set of primers resulted in a highly insignificant p-value of 0.600, making it unlikely that this locus is associated with IE in this breed. Two other loci, both human epilepsy-associated genes, though insignificant, had lower p-values: KCNQ3 at 0.077 and LGI1 at 0.070. Again, applying correction for multiple testing indicates that these p-values are not significant.

Discussion
This study was designed to use LD to exploit a possible founder effect suspected to be responsible for IE in these breeds by testing as many candidate genes as possible using microsatellite markers. Microsatellites have been shown to be a practical and useful tool in candidate gene studies [29]. Marker locations were selected within average blocks of canine LD so that one marker's result could, in theory, accurately represent the entire gene. The use of multiple breeds allowed for the possibility of assessing whether a discovered mutation was unique to one breed, suggesting the suspected founder effect specific to that breed, or was observed in several breeds, indicating a much older mutation. While no mutations were identified to draw definitive conclusions, it is not unreasonable to speculate that in dogs the underlying genetic basis of IE may vary by breed, as it is known that the rare human monogenic IEs are typically isolated by family [30][31][32][33][34][35]. These experiments incorporated both linkage analysis and association analysis, each having different strengths. While linkage analysis is the more powerful of these two statistical analysis methods for detecting Mendelian disease mutations, it requires samples from appropriate family structures. Conversely, association analysis does not require family-based samples, and it is generally considered more powerful than linkage analysis in detecting polygenic disease associations. Sequencing of the subset of human IE-associated candidate genes was not pursued for this study, as brain Bolded genes marked with a * indicate those associated with human IE, and bolded genes marked with a # indicate those associated with mouse models of epilepsy. The breed cohort was assembled and statistics were performed as described in Materials and Methods. "Designator" indicates each set of primers, first, second, etc., designed for that gene and used in this breed (See Additional files 1 and 2). H = heterozygosity of the marker.  Bolded genes marked with a * indicate those associated with human IE, and bolded genes marked with a # indicate those associated with mouse models of epilepsy. The breed cohort was assembled and statistics were performed as described in Materials and Methods. "Designator" indicates each set of primers, first, second, etc., designed for that gene and used in this breed (See Additional files 1 and 2). H = heterozygosity of the marker.    tissues for obtaining cDNAs for most cost efficient sequencing were not available.
Overall marker informativeness was highly influenced by sample population. The Vizslas and Beagles had excellent average heterozygosities (0.569 and 0.574, respectively), while the ESS and GSMD were much lower (0.425 and 0.341, respectively). This may be due, in part, to the sample cohorts: the ESS and GSMD were set up as discordant full-or half-sibling pairs, whereas the Vizslas were constructed as large pedigrees, and the Beagles were case/control matched pairs that did not share common ancestors to at least the grandparent level. Reduced marker heterozygosity across breeds could also be influenced by a higher degree of inbreeding within the ESS and GSMD sample cohorts in this study, and the lower average heterozygosities of these two breeds may suggest a founder effect that could eventually help uncover associations. However, marker informativeness is not entirely dependent on sample population, as evidenced by CACNA1A, which was polymorphic enough to be confidently insignificant after analyzing one or two microsatellites in ESS, GSMD, and Beagles, but remained inconclusive through three markers in the Vizslas and ultimately required a fourth microsatellite to achieve a suitable heterozygosity in this breed. The accepted marker heterozygosity of 0.3 is perhaps low, especially when analyzing only a single microsatellite for each candidate gene, however, for many of the insignificant microsatellite results in all four breeds the heterozygosity was > 0.5.
Ultimately, 16 of the 20 human epilepsy-associated candidate genes (ARX, CACNA1A, CACNA1H, CACNB4, CHRNA4, CHRNB2, CLCN2, GABRA1, GABRG2, KCNQ2, KCNQ3, LGI1, NHLRC1, SCN1A, SCN1B, and SCN2A) and both of the mouse model genes (CACNA2D2 and CACNG2) were either excluded (with linkage) or demonstrated insignificant association to IE in all four of the breed cohorts. Another two of them (KCNA1 and ME2) were insignificant in three out of four breed cohorts, remaining inconclusive only in the GSMD. One gene, GABRD, was inconclusive in two breeds (ESS and GSMD), and CHRNA2 was inconclusive in ESS and was not tested in Vizslas. Insignificant association to IE was demonstrated in most of the additional candidate genes tested. For every breed, however, there were a handful of markers that were inconclusive, due to low heterozygosity.
The few markers with potentially interesting results (uncorrected p-values of < 0.05) that were not excluded by follow-up markers, including CACNB1 in the Vizsla, CHRNB2 in the GSMD, and KCNQ3 and LGI1 in the Beagle, were of minimal interest when correcting for multiple testing. With the Bonferroni correction, a p-value of 0.05 considered as statistically significant would be lowered to 0.0025 for twenty tests, and further lowered to 0.00125 for forty tests; p-values of this magnitude were not obtained for any marker, and all breeds were tested on more than forty genes. The potential of false positives due to population stratification must also be considered for these few genes with suggestive results in the association studies, despite attempts to control for the degree of relatedness. Haplotype association analysis was not possible with this data since there were not Bolded genes marked with a * indicate those associated with human IE, and bolded genes marked with a # indicate those associated with mouse models of epilepsy. The breed cohort was assembled and statistics were performed as described in Materials and Methods. "Designator" indicates each set of primers, first, second, etc., designed for that gene and used in this breed (See Additional files 1 and 2). H = heterozygosity of the marker.
enough closely-spaced markers to generate haplotypes. Conclusive confirmation or exclusion of these loci can be performed in the future with newer technologies such as whole-genome SNP arrays, which can generate vast amounts of data in less time. It is possible that the populations tested were underpowered to detect association if it existed. Utilizing a siblingpair case-control design could decrease power because the control dog may also carry the risk allele. However, this situation is difficult to avoid in highly inbred dog populations, and the sibling-pair case-control design should aid in avoiding population stratification, which can create false positive results. Although sample size requirements for canine association studies have not been precisely defined, based on average linkage disequilibrium estimates in dog breeds, a starting point of approximately 25 cases and 25 controls appears to be adequate to find statistical significance for a completely penetrant recessive trait [36]. The number increases to approximately 50 cases and 50 controls when a trait is dominant, and for more complex traits as many as 100 dogs in each group may be necessary. Therefore, it is possible that there were not adequate numbers of dogs in the association studies, particularly if IE is not monogenic. A meta-analysis might improve power by pooling data from the three association study breeds (ESS, GSMD, and Beagle). However, this proved to be impractical because 1) many microsatellites had varying informativeness in each breed, so that pooling data would only work for a single microsatellite that was informative for a gene in all three breeds and 2) microsatellites by their very nature are highly polymorphic and there was often little overlap in alleles between breeds. Lastly, if IE has a different genetic basis in each of these breeds, pooling data across breeds seems unlikely to yield interesting results.
In the linkage analysis of Vizsla pedigrees, it is possible to conclude that these candidate gene loci are truly excluded, due to convincingly negative LOD scores at zero centiMorgans. However, if the inheritance and age dependent penetrance assumptions for the Vizsla model are incorrect, then the present linkage study is flawed. These input assumptions were based on previous study of IE in this breed [23] and on the best information available. Conversely, in the association studies of the other three breeds, it is not possible to specifically state that the loci are "excluded", as there is some degree of doubt that one microsatellite marker is adequate to provide convincing evidence for exclusion of a locus with this type of study. The addition of more markers would help to verify these insignificant results and truly exclude a locus. That work is beyond the scope of this study and is better addressed with whole-genome SNP analysis. The latest commercial canine SNP array, with over 170,000 SNPs, became available well after these studies were initiated. Further, the excluded and insignificant genes reported in the present study may only reflect the specific breeds examined, or even lines within these breeds, and results should not necessarily be extrapolated to other breeds.
It is possible that canine IE, like the majority of human IE, is a genetically complex disease in most breeds and that multiple loci contribute to susceptibility in any given breed. The present study would very likely have detected a major contributing locus, and it is unlikely that a truly causative locus is being excluded as a false negative. In humans, the vast majority of IE remains genetically unexplained and is considered to be polygenic [37][38][39]. A recent study by Oberbauer et al. [40] utilized microsatellites in a genome-wide linkage scan for epilepsy loci in the Belgian shepherd dog and concluded that the disease was highly polygenic, reporting a tentative six QTLs. These results further support the conclusion that canine IE is a more complex disease than originally hypothesized. Whole-genome association analyses with SNP arrays are likely the platform of choice for further studies of IE, as they can query tens of thousands of markers simultaneously across the genome and are better able to identify multiple susceptibility loci. Additionally, copy number variants (CNVs) have been increasingly shown to be involved in neurologic disorders such as autism [41,42] and schizophrenia [43,44], as well as epilepsy [45]. CNV studies of canine epilepsy may reveal this as a similar mechanism for disease in both species.
Ultimately, canine IE found to be significantly associated with DNA markers or with a mutation in a specific gene would allow genetic tests to be developed to assist dog breeders with decreasing the incidence of this disease. This would be most effective for monogenic, highly penetrant forms of IE, but if the disease proves to be genetically complex, with multiple genes contributing and less than 100% penetrance, then such a test could still possibly be used to provide a relative risk for IE development within a breed. In any event, discovery of IE associated gene loci in dogs may not only improve the understanding of canine health, but could advance the study of neurobiology and human health as well.

Conclusions
The present study did not identify any obvious association or linkage between IE and the microsatellite markers tested for the various candidate genes in any of the four dog breeds. It is possible that IE may in fact be a polygenic disorder even in highly inbred dog breeds, which underscores the challenging nature of investigating inherited epilepsy. Eventually identifying the genetic causes for IE in dogs would have a significant impact on canine health, as well as providing a useful model for IE in humans.

Case Definition
Since seizure activity can result from numerous discernable abnormalities, and only dogs very likely to be affected with IE were to be included in the study, the case definition was designed to eliminate, as much as possible, all other likely causes of seizures. Therefore, in order to qualify for "case" status for IE, each affected dog must have had two or more seizures with clinical normalcy between episodes, a veterinarian-administered normal physical examination and neurological examination, as well as normal bloodwork, including a complete blood count and blood chemistry profile, and, when possible, a bile acids assay. Since IE has been shown to be statistically much more probable when the dog is aged 1 to 5 years at the time of its first seizure [46], only dogs whose first seizure occurred during this time frame were included. If seizures started prior to 1 year of age or after 5 years of age, extra criteria were then required to enroll a dog, including a normal MRI or CT, and/or normal CSF analysis. EEG was not used as a case criterion due to its limited availability in clinical veterinary practice. To be classified as unaffected, the dog must have reached at least five years of age and be completely free of known seizure activity. In the situation where a case was identified at a young age, and its full-or half-sibling was used as a control, the dogs were followed for years until the unaffected dog had reached at least five years of age.
Dogs in these studies were privately owned and voluntarily enrolled by their owners with informed consent. Cases were solicited through breed publications, presentations at local and national breed club meetings, and word of mouth. For each enrolled dog, a general health survey, a signed owner consent form, and up to ten mLs of whole blood in an EDTA tube were required. Blood samples were obtained via venipuncture either at the dog's regular veterinarian or at specially-organized blood-draws at breed shows across the USA. This study was approved by the University of Minnesota IACUC.
Owners of dogs affected with seizures were asked to fill out an additional questionnaire which asked for information about the dog's seizure characteristics, preictal and postictal signs, previous medical history, and any anti-epileptic drug therapy; these were submitted by phone, fax, postal mail or e-mail. Any additional followup questions were conducted via telephone interview. Finally, medical records for each affected dog were obtained from the dog's regular veterinarian and any specialists (if applicable); these were inspected for bloodwork results and the veterinarian's physical and neurological examination findings. Dogs were assigned case, control, or unknown status after careful evaluation of this information. "Unknown" status was reserved for those dogs displaying suspicious seizure or seizure-like activity, but couldn't be classified as a case for various reasons; for example, if the dog had only one recorded seizure or had suffered head trauma prior to onset of the first seizure.

DNA Isolation
DNA was isolated from whole blood following standard protocols, utilizing the Puregene DNA Isolation Kit from Qiagen Inc. (Valencia, CA).

Breed Cohorts
Cohorts of affected and unaffected dogs for genotyping were assembled for each of the four test breeds -the Vizsla, the GSMD, the ESS, and the Beagle. Ultimately, the available samples determined what type of genetic analysis was conducted for each breed.

Vizslas
For Vizslas, the depth of the available DNA samples from multiple siblings, parents, and many grandparents allowed the development of multi-generation pedigrees ( Figure 1 and Figure 2) suitable for genetic linkage analysis. 96 dogs were included: 31 affected dogs, 60 unaffected dogs, and 5 dogs with unknown phenotype status. These dogs were separated into two large family pedigrees shown in Figures 1 and 2, and were subsequently broken into ten smaller families to decrease inbreeding loops for linkage analysis and calculation of LOD scores.

Greater Swiss Mountain Dogs
For the GSMD, discordant full-sibling or half-sibling pairs of affected (n = 24) and unaffected (n = 24) dogs were assembled and used in association studies. Siblings were used to help avoid spurious association from population stratification. In a few instances, two different pairs of siblings were related to each other, such as sharing one parent.

English Springer Spaniels
Similar to the GSMD, the ESS were assembled into discordant full-or half-sibling pairs of affected (n = 24) and unaffected (n = 24) dogs and were used in association studies. As in the GSMD, occasionally the pairs of siblings were also related to one another; this was impossible to avoid due to the inbred nature of the breed.

Beagles
Association studies were also conducted for the Beagles. Due to lack of available sibling pairs, pedigree analyses were used to ensure that only dogs with no common ? ? ? ? Figure 1 Vizsla linkage family 1 pedigree. Pedigree of Vizsla family 1. Squares and circles represent males and females, respectively. Filled shapes represent cases; those with question marks represent unknown phenotype status. Dogs with arrows were genotyped in this study. The dog marked with an arrow and a + is included on both families 1 & 2, effectively making this one very large family. The Vizsla pedigrees were broken into ten smaller families to decrease inbreeding loops before being analyzed in linkage analysis. Three dogs representing one of the ten sub-families are not shown on either pedigree.
? Figure 2 Vizsla linkage family 2 pedigree. Pedigree of Vizsla family 2. Squares and circles represent males and females, respectively. Filled shapes represent cases; those with question marks represent unknown phenotype status. Dogs with arrows were genotyped in this study. The dog marked with an arrow and a + is included on both families 1 & 2, effectively making this one very large family. The Vizsla pedigrees were broken into ten smaller families to decrease inbreeding loops before being analyzed in linkage analysis. Three dogs representing one of the ten sub-families are not shown on either pedigree. grandparents were utilized. This strategy has previously shown to be adequate for canine genetic mutation detection [47]. Twenty-four affected dogs and 24 unaffected dogs were selected for this dataset.
In addition to this core set of genes, further genes were selected from the same gene families based on their presence in neural tissue, including: calcium channels CACNA1B, CACNA1D, CACNA1E, CACNA1F, CACNA1G, CACNA1I, CACNB1, CACNB2, CACNB3, CACNG3, CACNG4, and CACNG6; sodium channels SCN3A, SCN3B, SCN8A, and SCN11A; GABA receptors GABRA2, and GABRA6; acetylcholine receptors CHRNA1, CHRNA3, CHRNA5, CHRNA9, CHRNA10, CHRNB1, and CHRND as well as CHRNA7, which has had suggested linkage with human IE [57]; and potassium channels KCND2, KCNQ1, and KCNQ5. Finally, the dynamin-1 gene (DNM1) was also examined. This gene was recently found to be mutated in the canine disease termed exercise induced collapse [58] and a mutation in dynamin-1 has recently been characterized in the fitful mouse model, wherein heterozygotes for the mutation experience recurrent seizures and homozygotes for the mutation often have lethal seizures and other neurosensory deficits [59]. Therefore, a total of 52 candidate genes were chosen for this study.

Primer design, PCR, and Genotyping
Human mRNA sequence was obtained for each gene from the National Center for Biotechnology Information [60]. This sequence was BLATed on the University of California -Santa Cruz Genome Browser website [61] against the May 2005 dog genome assembly to identify the resultant genomic position of each gene in the dog. Microsatellites were selected either within the gene or closely flanking the gene utilizing the Variation and Repeat -Microsatellite feature of this browser. After capturing the microsatellite's flanking sequence, primer pairs to amplify each microsatellite were designed using Primer3 [62]. An additional 5' tail sequence was added to each right (reverse) primer [63] to generate fluorescent labels in PCR. All primer sequences are shown in Additional file 1, and Additional file 2 lists the Mb positions of each gene and microsatellite. Every attempt was made to select microsatellites less than 1 Mb distant from the gene, in order to remain within estimated blocks of canine linkage disequilibrium (LD) [36] for each gene. LD is the non-random association of alleles at two loci and, within dog breeds, can extend over distances of several megabases [36]. However, for many of the candidate genes, several microsatellites (and therefore several sets of primers) needed to be tested due to the microsatellite's lack of informativeness (monomorphism or low heterozygosity), and eventually the only useable microsatellite was located further away. Additional file 2 indicates which microsatellites (n = 3) were >1 Mb distant. CHRNA4 and KCNQ2 were close enough to one another in the canine genome that one marker was used for both genes.
Each microsatellite was PCR amplified in 15-μL reactions including 12.5 ng canine DNA, 10X PCR Buffer, 2.5 pmol of the forward primer and 0.75 pmol of the reverse-tailed primer, 0.78 pmol of a fluorescent dyelabeled primer complementary to the tail sequence, 20 μM dNTP mix, and 0.5 units of Taq DNA polymerase (Qiagen Hot StarTaq, Qiagen Inc, Valencia, CA), with dH 2 O to final volume. PCR cycling conditions were 95°C for 20 minutes, then 40 cycles of the following: 94°C for 30 seconds, 58°C for 30 seconds, and 72°C for 30 to 60 seconds (depending on product size), and a final extension of 72°C for 15 minutes. PCR products were prepared for capillary gel electrophoresis according to manufacturer's recommendations and size separated using the Beckman CEQ 8000 automated DNA-fragment analyzer with fluorescence detection (Beckman Coulter Inc, Fullerton CA). Genotype data were analyzed with instrument software, and alleles were manually verified.

Statistical Analysis
Marker heterozygosity was calculated using the following formula: H = 1 -Σp i 2 , where p i is the population frequency of the ith allele, and H is the probability that a random individual is heterozygous for any two alleles at a gene locus with allele frequencies, p i [64]. A cut-off of 0.3 was selected for heterozygosity; markers with heterozygosity falling below this cut-off were not considered informative enough to confidently claim a non-significant result for a gene, since less than one third of all individuals would be expected to be heterozygous at that locus [64].
In order to determine the power of the ten Vizsla family pedigrees to detect linkage, simulated linkage was performed with the FASTSLINK program [65][66][67]. Multiple parameters were used, including combinations of 3 and 4 equally frequent alleles, 0 and 5% heterogeneity (proportion of unlinked families) and recombination frequencies of 0.01 and 0.001. Recombination frequencies were kept very small since the microsatellites tested were within or very near to the candidate genes. Two hundred simulations were performed for each set of parameters, assuming autosomal recessive inheritance.
Parametric two-point linkage analysis was utilized for Vizsla marker data, with significant evidence for linkage a LOD (logarithm of odds) score of 3.2 or greater [68]. A LOD score of less than -2.0 was considered to exclude a locus. Input assumptions were an autosomal recessive mode of inheritance with 95% penetrance for the disease genotype, and a disease allele frequency of 0.20; these assumptions were based on previously published data concerning IE in this breed [23]. Two-point LOD scores were calculated using the FASTLINK version of MLINK [69][70][71][72].
For the GSMD, ESS, and Beagles, the allele frequencies were subjected to a nonparametric chi-square test for independence (association analysis) [73]. Alleles with frequencies of less than 10% (minor alleles) were removed before conducting the chi-square analysis. A Pearson chi-square p-value of 0.05 was first used as the level of significance necessary to detect significantly different allele frequencies between affected and unaffected dogs. Subsequently, a Bonferroni corrected p-value of 0.001, correcting for~52 tests, was considered evidence for strong association. Power analyses were not specifically conducted for these three breeds, as sample size requirements for canine association studies have not been precisely defined.

Additional material
Additional file 1: Primer sequences for candidate gene microsatellites. Sequence information and product sizes are presented for primers designed for microsatellites examined for each candidate gene. Bolded genes marked with a * indicate those associated with human epilepsy, and bolded genes marked with a # indicate those associated with mouse models of epilepsy. Two of the DNM1 marker primer pairs, DNM1 Third Set and DNM1 Fifth Set, are taken from the University of California-Davis Canine Genetic Linkage Map [74], and are named therein 0945 and 0946, respectively. CFA = canis familiaris chromosome.
Additional file 2: Candidate gene locations and primer locations. All pertinent canine chromosomal location information is presented for each gene and all microsatellites examined for that gene. Bolded genes marked with a * indicate those associated with human epilepsy, and bolded genes marked with a # indicate those associated with mouse models of epilepsy. Two of the DNM1 marker primer pairs, DNM1 Third Set and DNM1 Fifth Set, are taken from the University of California-Davis Canine Genetic Linkage Map [74], and are named therein 0945 and 0946, respectively.^Bolded entries in this column indicate microsatellites located outside the gene itself and further than 1 Mb away (n = 3). Mb positions given for microsatellites are the first nucleotide of each microsatellite. Mb positions were taken from Ensembl [75] and the UCSC Genome Browser [61]. CFA = canis familiaris chromosome.