Genetics of immune response to Epstein-Barr virus: prospects for multiple sclerosis pathogenesis

Epstein-Barr virus (EBV) infection has been advocated as a prerequisite for developing multiple sclerosis (MS) and possibly the propagation of the disease. However, the precise mechanisms for such influences are still unclear. A large-scale study investigating the host genetics of EBV serology and related clinical manifestations, such as infectious mononucleosis (IM), may help us better understand the role of EBV in MS pathogenesis. This study evaluates the host genetic factors that influence serological response against EBV and history of IM and cross-evaluates them with MS risk and genetic susceptibility in the Swedish population. Plasma IgG antibody levels against EBV nuclear antigen-1 (EBNA-1, truncated=aa[325-641], peptide=aa[385-420]) and viral capsid antigen p18 (VCAp18) were measured using bead-based multiplex serology for 8744 MS cases and 7229 population-matched controls. The MS risk association for high/low EBV antibody levels and history of IM was compared to relevant clinical measures along with sex, age at sampling, and associated HLA allele variants. Genome-wide and HLA allele association analyses were also performed to identify genetic risk factors for EBV antibody response and IM history. Higher antibody levels against VCAp18 (OR=1.74, 95% CI=1.60-1.88) and EBNA-1, particularly the peptide (OR=3.13, 95% CI=2.93-3.35), were associated with an increased risk for MS. The risk increased with higher anti-EBNA-1 IgG levels up to twelve times the reference risk. We also identified several independent HLA haplotypes associated with EBV serology overlapping with known MS risk alleles (e.g., DRB1*15:01). Although there were several candidates, no variants outside the HLA region reached genome-wide significance. Cumulative HLA risk for anti-EBNA-1 IgG levels, particularly the peptide fragment, was strongly associated with MS. In contrast, the genetic risk for high anti-VCAp18 IgG levels was not as strongly associated with MS risk. IM history was not associated with class II HLA genes but negatively associated with A*02:01, which is protective against MS. Our findings emphasize that the risk association between anti-EBNA-1 IgG levels and MS may be partly due to overlapping HLA associations. Additionally, the increasing MS risk with increasing anti-EBNA-1 levels would be consistent with a pathogenic role of the EBNA-1 immune response, perhaps through molecular mimicry. Given that high anti-EBNA-1 antibodies may reflect a poorly controlled T-cell defense against the virus, our findings would be consistent with DRB1*15:01 being a poor class II antigen in the immune defense against EBV. Lastly, the difference in genetic control of IM supports the independent roles of EBNA-1 and IM in MS susceptibility.

Cross-reactive EBNA1 immunity targets alpha-crystallin B and is associated with multiple sclerosis

Multiple sclerosis (MS) is an inflammatory disease of the central nervous system, for which and Epstein-Barr virus (EBV) infection is a likely prerequisite. Due to the homology between Epstein-Barr nuclear antigen 1 (EBNA1) and alpha-crystallin B (CRYAB), we examined antibody reactivity to EBNA1 and CRYAB peptide libraries in 713 persons with MS (pwMS) and 722 matched controls (Con). Antibody response to CRYAB amino acids 7 to 16 was associated with MS (OR = 2.0), and combination of high EBNA1 responses with CRYAB positivity markedly increased disease risk (OR = 9.0). Blocking experiments revealed antibody cross-reactivity between the homologous EBNA1 and CRYAB epitopes. Evidence for T cell cross-reactivity was obtained in mice between EBNA1 and CRYAB, and increased CRYAB and EBNA1 CD4+ T cell responses were detected in natalizumab-treated pwMS. This study provides evidence for antibody cross-reactivity between EBNA1 and CRYAB and points to a similar cross-reactivity in T cells, further demonstrating the role of EBV adaptive immune responses in MS development.

Bayesian mixed model analysis uncovered 21 risk loci for chronic kidney disease in boxer dogs

Chronic kidney disease (CKD) affects 10% of the human population, with only a small fraction genetically defined. CKD is also common in dogs and has been diagnosed in nearly all breeds, but its genetic basis remains unclear. Here, we performed a Bayesian mixed model genome-wide association analysis for canine CKD in a boxer population of 117 canine cases and 137 controls, and identified 21 genetic regions associated with the disease. At the top markers from each CKD region, the cases carried an average of 20.2 risk alleles, significantly higher than controls (15.6 risk alleles). An ANOVA test showed that the 21 CKD regions together explained 57% of CKD phenotypic variation in the population. Based on whole genome sequencing data of 20 boxers, we identified 5,206 variants in LD with the top 50 BayesR markers. Following comparative analysis with human regulatory data, 17 putative regulatory variants were identified and tested with electrophoretic mobility shift assays. In total four variants, three intronic variants from the MAGI2 and GALNT18 genes, and one variant in an intergenic region on chr28, showed alternative binding ability for the risk and protective alleles in kidney cell lines. Many genes from the 21 CKD regions, RELN, MAGI2, FGFR2 and others, have been implicated in human kidney development or disease. The results from this study provide new information that may enlighten the etiology of CKD in both dogs and humans.

Identification of four novel T cell autoantigens and personal autoreactive profiles in multiple sclerosis

Multiple sclerosis (MS) is an inflammatory disease of the central nervous system (CNS), in which pathological T cells, likely autoimmune, play a key role. Despite its central importance, the autoantigen repertoire remains largely uncharacterized. Using a novel in vitro antigen delivery method combined with the Human Protein Atlas library, we screened for T cell autoreactivity against 63 CNS-expressed proteins. We identified four previously unreported autoantigens in MS: fatty acid-binding protein 7, prokineticin-2, reticulon-3, and synaptosomal-associated protein 91, which were verified to induce interferon-γ responses in MS in two cohorts. Autoreactive profiles were heterogeneous, and reactivity to several autoantigens was MS-selective. Autoreactive T cells were predominantly CD4⁺ and human leukocyte antigen-DR restricted. Mouse immunization induced antigen-specific responses and CNS leukocyte infiltration. This represents one of the largest systematic efforts to date in the search for MS autoantigens, demonstrates the heterogeneity of autoreactive profiles, and highlights promising targets for future diagnostic tools and immunomodulatory therapies in MS.

The ABCC4 gene is associated with pyometra in golden retriever dogs

Pyometra is one of the most common diseases in female dogs, presenting as purulent inflammation and bacterial infection of the uterus. On average 20% of intact female dogs are affected before 10 years of age, a proportion that varies greatly between breeds (3–66%). The clear breed predisposition suggests that genetic risk factors are involved in disease development. To identify genetic risk factors associated with the disease, we performed a genome-wide association study (GWAS) in golden retrievers, a breed with increased risk of developing pyometra (risk ratio: 3.3). We applied a mixed model approach comparing 98 cases, and 96 healthy controls and identified an associated locus on chromosome 22 (p = 1.2 × 10^–6, passing Bonferroni corrected significance). This locus contained five significantly associated SNPs positioned within introns of the ATP-binding cassette transporter 4 (ABCC4) gene. This gene encodes a transmembrane transporter that is important for prostaglandin transport. Next generation sequencing and genotyping of cases and controls subsequently identified four missense SNPs within the ABCC4 gene. One missense SNP at chr22:45,893,198 (p.Met787Val) showed complete linkage disequilibrium with the associated GWAS SNPs suggesting a potential role in disease development. Another locus on chromosome 18 overlapping the TESMIN gene, is also potentially implicated in the development of the disease.

Transcriptomes from German shepherd dogs reveal differences in immune activity between atopic dermatitis affected and control skin

Canine atopic dermatitis (CAD) is an inflammatory and pruritic allergic skin disease with both genetic and environmental risk factors described. We performed mRNA sequencing of non-lesional axillary skin biopsies from nine German shepherd dogs. Obtained RNA sequences were mapped to the dog genome (CanFam3.1) and a high-quality skin transcriptome was generated with 23,510 expressed gene transcripts. Differentially expressed genes (DEGs) were defined by comparing three controls to five treated CAD cases. Using a leave-one-out analysis, we identified seven DEGs: five known to encode proteins with functions related to an activated immune system (CD209, CLEC4G, LOC102156842 (lipopolysaccharide-binding protein-like), LOC480601 (regakine-1-like), LOC479668 (haptoglobin-like)), one (OBP) encoding an odorant-binding protein potentially connected to rhinitis, and the last (LOC607095) encoding a novel long non-coding RNA. Furthermore, high mRNA expression of inflammatory genes was found in axillary skin from an untreated mild CAD case compared with healthy skin. In conclusion, we define genes with different expression patterns in CAD case skin helping us understand post-treatment atopic skin. Further studies in larger sample sets are warranted to confirm and to transfer these results into clinical practice.

Comparison of cellular location and expression of Plakophilin-2 in epidermal cells from nonlesional atopic skin and healthy skin in German shepherd dogs

Background

Canine atopic dermatitis (CAD) is an inflammatory and pruritic allergic skin disease caused by interactions between genetic and environmental factors. Previously, a genome‐wide significant risk locus on canine chromosome 27 for CAD was identified in German shepherd dogs (GSDs) and Plakophilin‐2 (PKP2) was defined as the top candidate gene. PKP2 constitutes a crucial component of desmosomes and also is important in signalling, metabolic and transcriptional activities.

Objectives

The main objective was to evaluate the role of PKP2 in CAD by investigating PKP2 expression and desmosome structure in nonlesional skin from CAD‐affected (carrying the top GWAS SNP risk allele) and healthy GSDs. We also aimed at defining the cell types in the skin that express PKP2 and its intracellular location.

Animals/Methods

Skin biopsies were collected from nine CAD‐affected and five control GSDs. The biopsies were frozen for immunofluorescence and fixed for electron microscopy immunolabelling and morphology.

Results

We observed the novel finding of PKP2 expression in dendritic cells and T cells in dog skin. Moreover, we detected that PKP2 was more evenly expressed within keratinocytes compared to its desmosomal binding‐partner plakoglobin. PKP2 protein was located in the nucleus and on keratin filaments attached to desmosomes. No difference in PKP2 abundance between CAD cases and controls was observed.

Conclusion

Plakophilin‐2 protein in dog skin is expressed in both epithelial and immune cells; based on its subcellular location its functional role is implicated in both nuclear and structural processes.

Multiple regulatory variants located in cell type-specific enhancers within the PKP2 locus form major risk and protective haplotypes for canine atopic dermatitis in German shepherd dogs

Background

Canine atopic dermatitis (CAD) is a chronic inflammatory skin disease triggered by allergic reactions involving IgE antibodies directed towards environmental allergens. We previously identified a ~1.5 Mb locus on canine chromosome 27 associated with CAD in German shepherd dogs (GSDs). Fine-mapping indicated association closest to the PKP2 gene encoding plakophilin 2.

Results

Additional genotyping and association analyses in GSDs combined with control dogs from five breeds with low-risk for CAD revealed the top SNP 27:19,086,778 (p = 1.4 × 10⁻⁷) and a rare ~48 kb risk haplotype overlapping the PKP2 gene and shared only with other high-risk CAD breeds. We selected altogether nine SNPs (four top-associated in GSDs and five within the ~48 kb risk haplotype) that spanned ~280 kb forming one risk haplotype carried by 35 % of the GSD cases and 10 % of the GSD controls (OR = 5.1, p = 5.9 × 10⁻⁵), and another haplotype present in 85 % of the GSD cases and 98 % of the GSD controls and conferring a protective effect against CAD in GSDs (OR = 0.14, p = 0.0032). Eight of these SNPs were analyzed for transcriptional regulation using reporter assays where all tested regions exerted regulatory effects on transcription in epithelial and/or immune cell lines, and seven SNPs showed allelic differences. The DNA fragment with the top-associated SNP 27:19,086,778 displayed the highest activity in keratinocytes with 11-fold induction of transcription by the risk allele versus 8-fold by the control allele (p_difference = 0.003), and also mapped close (~3 kb) to an ENCODE skin-specific enhancer region.

Conclusions

Our experiments indicate that multiple CAD-associated genetic variants located in cell type-specific enhancers are involved in gene regulation in different cells and tissues. No single causative variant alone, but rather multiple variants combined in a risk haplotype likely contribute to an altered expression of the PKP2 gene, and possibly nearby genes, in immune and epithelial cells, and predispose GSDs to CAD.

Electronic supplementary material

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

Naturally Occurring Adrenocortical Insufficiency--An Epidemiological Study Based on a Swedish-Insured Dog Population of 525,028 Dogs

Background

Naturally occurring adrenocortical insufficiency (NOAI) in dogs is considered an uncommon disease with good prognosis with hormonal replacement treatment. However, there are no epidemiological studies with estimates for the general dog population.

Objectives

To investigate the epidemiological characteristics of NOAI in a large population of insured dogs.

Animals

Data were derived from 525,028 client‐owned dogs insured by a Swedish insurance company representing 2,364,652 dog‐years at risk (DYAR) during the period between 1995–2006.

Methods

Retrospective cohort study. Incidence rates, prevalences, and relative risks for dogs with NOAI (AI with no previous claim for hypercortisolism), were calculated for the whole dog population, and for subgroups divided by breed and sex. Mortality rates were calculated and compared in dogs with NOAI and the remaining dogs overall.

Results

In total 534 dogs were identified with NOAI. The overall incidence was 2.3 cases per 10,000 DYAR. The relative risk of disease was significantly higher in the Portuguese Water Dog, Standard Poodle, Bearded Collie, Cairn Terrier, and Cocker Spaniel compared with other breeds combined. Female dogs overall were at higher risk of developing AI than male dogs (RR 1.85; 95% CI, 1.55–2.22; P < .001). The relative risk of death was 1.9 times higher in dogs with NOAI than in dogs overall.

Conclusion and Clinical Importance

The data supports the existence of breed‐specific differences in incidence rates of NOAI in dogs.

cgmisc: enhanced genome-wide association analyses and visualization

High-throughput genotyping and sequencing technologies facilitate studies of complex genetic traits and provide new research opportunities. The increasing popularity of genome-wide association studies (GWAS) leads to the discovery of new associated loci and a better understanding of the genetic architecture underlying not only diseases, but also other monogenic and complex phenotypes. Several softwares are available for performing GWAS analyses, R environment being one of them.

RESULTS

We present cgmisc, an R package that enables enhanced data analysis and visualization of results from GWAS. The package contains several utilities and modules that complement and enhance the functionality of the existing software. It also provides several tools for advanced visualization of genomic data and utilizes the power of the R language to aid in preparation of publication-quality figures. Some of the package functions are specific for the domestic dog (Canis familiaris) data.

AVAILABILITY AND IMPLEMENTATION

The package is operating system-independent and is available from: https://github.com/cgmisc-team/cgmisc

CONTACT

marcin.kierczak@imbim.uu.se.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-Wide Analyses Suggest Mechanisms Involving Early B-Cell Development in Canine IgA Deficiency

Immunoglobulin A deficiency (IgAD) is the most common primary immune deficiency disorder in both humans and dogs, characterized by recurrent mucosal tract infections and a predisposition for allergic and other immune mediated diseases. In several dog breeds, low IgA levels have been observed at a high frequency and with a clinical resemblance to human IgAD. In this study, we used genome-wide association studies (GWAS) to identify genomic regions associated with low IgA levels in dogs as a comparative model for human IgAD. We used a novel percentile groups-approach to establish breed-specific cut-offs and to perform analyses in a close to continuous manner. GWAS performed in four breeds prone to low IgA levels (German shepherd, Golden retriever, Labrador retriever and Shar-Pei) identified 35 genomic loci suggestively associated (p <0.0005) to IgA levels. In German shepherd, three genomic regions (candidate genes include KIRREL3 and SERPINA9) were genome-wide significantly associated (p <0.0002) with IgA levels. A ~20kb long haplotype on CFA28, significantly associated (p = 0.0005) to IgA levels in Shar-Pei, was positioned within the first intron of the gene SLIT1. Both KIRREL3 and SLIT1 are highly expressed in the central nervous system and in bone marrow and are potentially important during B-cell development. SERPINA9 expression is restricted to B-cells and peaks at the time-point when B-cells proliferate into antibody-producing plasma cells. The suggestively associated regions were enriched for genes in Gene Ontology gene sets involving inflammation and early immune cell development.

The dog as a genetic model for immunoglobulin A (IgA) deficiency: identification of several breeds with low serum IgA concentrations

Immunoglobulin A (IgA) serves as the basis of the secretory immune system by protecting the lining of mucosal sites from pathogens. In both humans and dogs, IgA deficiency (IgAD) is associated with recurrent infections of mucosal sites and immune-mediated diseases. Low concentrations of serum IgA have previously been reported to occur in a number of dog breeds but no generally accepted cut-off value has been established for canine IgAD. The current study represents the largest screening to date of IgA in dogs in terms of both number of dogs (n=1267) and number of breeds studied (n=22). Serum IgA concentrations were quantified by using capture ELISA and were found to vary widely between breeds. We also found IgA to be positively correlated with age (p<0.0001). Apart from the two breeds previously reported as predisposed to low IgA (Shar-Pei and German shepherd), we identified six additional breeds in which ≥ 10% of all tested dogs had very low (<0.07 g/l) IgA concentrations (Hovawart, Norwegian elkhound, Nova Scotia duck tolling retriever, Bullterrier, Golden retriever and Labrador retriever). In addition, we discovered low IgA concentrations to be significantly associated with canine atopic dermatitis (CAD, p<0.0001) and pancreatic acinar atrophy (PAA, p=0.04) in German shepherds.

Thorough investigation of a canine autoinflammatory disease (AID) confirms one main risk locus and suggests a modifier locus for amyloidosis

Autoinflammatory disease (AID) manifests from the dysregulation of the innate immune system and is characterised by systemic and persistent inflammation. Clinical heterogeneity leads to patients presenting with one or a spectrum of phenotypic signs, leading to difficult diagnoses in the absence of a clear genetic cause. We used separate genome-wide SNP analyses to investigate five signs of AID (recurrent fever, arthritis, breed specific secondary dermatitis, otitis and systemic reactive amyloidosis) in a canine comparative model, the pure bred Chinese Shar-Pei. Analysis of 255 DNA samples revealed a shared locus on chromosome 13 spanning two peaks of association. A three-marker haplotype based on the most significant SNP (p<2.6×10⁻⁸) from each analysis showed that one haplotypic pair (H13-11) was present in the majority of AID individuals, implicating this as a shared risk factor for all phenotypes. We also noted that a genetic signature (F_ST) distinguishing the phenotypic extremes of the breed specific Chinese Shar-Pei thick and wrinkled skin, flanked the chromosome 13 AID locus; suggesting that breed development and differentiation has played a parallel role in the genetics of breed fitness. Intriguingly, a potential modifier locus for amyloidosis was revealed on chromosome 14, and an investigation of candidate genes from both this and the chromosome 13 regions revealed significant (p<0.05) renal differential expression in four genes previously implicated in kidney or immune health (AOAH, ELMO1, HAS2 and IL6). These results illustrate that phenotypic heterogeneity need not be a reflection of genetic heterogeneity, and that genetic modifiers of disease could be masked if syndromes were not first considered as individual clinical signs and then as a sum of their component parts.

A mutation in the SUV39H2 gene in Labrador Retrievers with hereditary nasal parakeratosis (HNPK) provides insights into the epigenetics of keratinocyte differentiation

Hereditary nasal parakeratosis (HNPK), an inherited monogenic autosomal recessive skin disorder, leads to crusts and fissures on the nasal planum of Labrador Retrievers. We performed a genome-wide association study (GWAS) using 13 HNPK cases and 23 controls. We obtained a single strong association signal on chromosome 2 (p(raw) = 4.4×10⁻¹⁴). The analysis of shared haplotypes among the 13 cases defined a critical interval of 1.6 Mb with 25 predicted genes. We re-sequenced the genome of one case at 38× coverage and detected 3 non-synonymous variants in the critical interval with respect to the reference genome assembly. We genotyped these variants in larger cohorts of dogs and only one was perfectly associated with the HNPK phenotype in a cohort of more than 500 dogs. This candidate causative variant is a missense variant in the SUV39H2 gene encoding a histone 3 lysine 9 (H3K9) methyltransferase, which mediates chromatin silencing. The variant c.972T>G is predicted to change an evolutionary conserved asparagine into a lysine in the catalytically active domain of the enzyme (p.N324K). We further studied the histopathological alterations in the epidermis in vivo. Our data suggest that the HNPK phenotype is not caused by hyperproliferation, but rather delayed terminal differentiation of keratinocytes. Thus, our data provide evidence that SUV39H2 is involved in the epigenetic regulation of keratinocyte differentiation ensuring proper stratification and tight sealing of the mammalian epidermis.

IgA deficiency in wolves

Low mean concentrations of serum immunoglobulin A (IgA) and an increased frequency of overt IgA deficiency (IgAD) in certain dog breeds raises the question whether it is a breeding-enriched phenomenon or a legacy from the dog's ancestor, the gray wolf (Canis lupus). The IgA concentration in 99 serum samples from 58 free-ranging and 13 captive Scandinavian wolves, was therefore measured by capture ELISA. The concentrations were markedly lower in the wolf serum samples than in the dog controls. Potential differences in the IgA molecule between dogs and wolves were addressed by sequencing the wolf IgA heavy chain constant region encoding gene (IGHA). Complete amino acid sequence homology was found. Detection of wolf and dog IgA was ascertained by showing identity using double immunodiffusion. We suggest that the vast majority of wolves, the ancestor of the dog, are IgA deficient.

Genome-wide analysis in German shepherd dogs reveals association of a locus on CFA 27 with atopic dermatitis

Humans and dogs are both affected by the allergic skin disease atopic dermatitis (AD), caused by an interaction between genetic and environmental factors. The German shepherd dog (GSD) is a high-risk breed for canine AD (CAD). In this study, we used a Swedish cohort of GSDs as a model for human AD. Serum IgA levels are known to be lower in GSDs compared to other breeds. We detected significantly lower IgA levels in the CAD cases compared to controls (p = 1.1 × 10(-5)) in our study population. We also detected a separation within the GSD cohort, where dogs could be grouped into two different subpopulations. Disease prevalence differed significantly between the subpopulations contributing to population stratification (λ = 1.3), which was successfully corrected for using a mixed model approach. A genome-wide association analysis of CAD was performed (n cases = 91, n controls = 88). IgA levels were included in the model, due to the high correlation between CAD and low IgA levels. In addition, we detected a correlation between IgA levels and the age at the time of sampling (corr = 0.42, p = 3.0 × 10(-9)), thus age was included in the model. A genome-wide significant association was detected on chromosome 27 (praw = 3.1 × 10(-7), pgenome = 0.03). The total associated region was defined as a ~1.5-Mb-long haplotype including eight genes. Through targeted re-sequencing and additional genotyping of a subset of identified SNPs, we defined 11 smaller haplotype blocks within the associated region. Two blocks showed the strongest association to CAD. The ~209-kb region, defined by the two blocks, harbors only the PKP2 gene, encoding Plakophilin 2 expressed in the desmosomes and important for skin structure. Our results may yield further insight into the genetics behind both canine and human AD.

Two loci on chromosome 5 are associated with serum IgE levels in Labrador retrievers

Crosslinking of immunoglobulin E antibodies (IgE) bound at the surface of mast cells and subsequent mediator release is considered the most important trigger for allergic reactions. Therefore, the genetic control of IgE levels is studied in the context of allergic diseases, such as asthma, atopic rhinitis, or atopic dermatitis (AD). We performed genome-wide association studies in 161 Labrador Retrievers with regard to total and allergen-specific immunoglobulin E (IgE) levels. We identified a genome-wide significant association on CFA 5 with the antigen-specific IgE responsiveness to Acarus siro. We detected a second genome-wide significant association with respect to the antigen-specific IgE responsiveness to Tyrophagus putrescentiae at a different locus on chromosome 5. A. siro and T. putrescentiae both belong to the family Acaridae and represent so-called storage or forage mites. These forage mites are discussed as major allergen sources in canine AD. No obvious candidate gene for the regulation of IgE levels is located under the two association signals. Therefore our studies offer a chance of identifying a novel mechanism controlling the host's IgE response.

Total IgE and allergen-specific IgE and IgG antibody levels in sera of atopic dermatitis affected and non-affected Labrador- and Golden retrievers

Canine atopic dermatitis (CAD) is an allergic skin disease associated with IgE and IgG antibodies (Ab) to environmental allergens. The aim of this study was to determine which other factors influence serum Ab levels in CAD-affected and non-affected dogs as this has only been poorly investigated in dogs so far. Total and allergen-specific IgE levels and Dermatophagoides farinae (DF)-specific IgG1 and IgG4 were measured by ELISA in sera of 145 CAD-affected and 271 non-affected Labrador- and Golden retrievers. A multivariable logistic regression analysis including the factors age, breed, gender, castration, clinical CAD status and allergen-specific immunotherapy (ASIT) was performed. Golden retrievers had more frequently total (OR=1.87, 95% CI=1.26-2.87, p<0.01) and specific IgE levels above the threshold value than Labrador retrievers, suggesting that genetic factors influence IgE levels in dogs. Castration was generally associated with low Ab levels (OR=0.43-0.65, p<0.05). Surprisingly, dogs with CAD did not have increased odds for high IgE against any of the allergens tested. ASIT with DF was associated with high DF-specific IgG1 (OR=4.32, 95% CI 1.46-12.8, p<0.01) but was not associated with DF-specific IgG4 or decreased IgE levels. Further studies are needed to understand the role of allergen-specific IgE in CAD and of IgG1 in ASIT.

Comparison of buccal and blood-derived canine DNA, either native or whole genome amplified, for array-based genome-wide association studies

Background

The availability of array-based genotyping platforms for single nucleotide polymorphisms (SNPs) for the canine genome has expanded the opportunities to undertake genome-wide association (GWA) studies to identify the genetic basis for Mendelian and complex traits. Whole blood as the source of high quality DNA is undisputed but often proves impractical for collection of the large numbers of samples necessary to discover the loci underlying complex traits. Further, many countries prohibit the collection of blood from dogs unless medically necessary thereby restricting access to critical control samples from healthy dogs. Alternate sources of DNA, typically from buccal cytobrush extractions, while convenient, have been suggested to have low yield and perform poorly in GWA. Yet buccal cytobrushes provide a cost-effective means of collecting DNA, are readily accepted by dog owners, and represent a large resource base in many canine genetics laboratories. To increase the DNA quantities, whole genome amplification (WGA) can be performed. Thus, the present study assessed the utility of buccal-derived DNA as well as whole genome amplification in comparison to blood samples for use on the most recent iteration of the canine HD SNP array (Illumina).

Findings

In both buccal and blood samples, whether whole genome amplified or not, 97% of the samples had SNP call rates in excess of 80% indicating that the vast majority of the SNPs would be suitable to perform association studies regardless of the DNA source. Similarly, there were no significant differences in marker intensity measurements between buccal and blood samples for copy number variations (CNV) analysis.

Conclusions

All DNA samples assayed, buccal or blood, native or whole genome amplified, are appropriate for use in array-based genome-wide association studies. The concordance between subsets of dogs for which both buccal and blood samples, or those samples whole genome amplified, was shown to average >99%. Thus, the two DNA sources were comparable in the generation of SNP genotypes and intensity values to estimate structural variation indicating the utility for the use of buccal cytobrush samples and the reliability of whole genome amplification for genome-wide association and CNV studies.