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.

Multiple Genetic Loci Associated with Pug Dog Thoracolumbar Myelopathy.

Pug dogs with thoracolumbar myelopathy (PDM) present with a specific clinical phenotype that includes progressive pelvic limb ataxia and paresis, commonly accompanied by incontinence. Vertebral column malformations and lesions, excessive scar tissue of the meninges, and central nervous system inflammation have been described. PDM has a late onset and affects more male than female dogs. The breed-specific presentation of the disorder suggests that genetic risk factors are involved in the disease development. To perform a genome-wide search for PDM-associated loci, we applied a Bayesian model adapted for mapping complex traits (BayesR) and a cross-population extended haplotype homozygosity test (XP-EHH) in 51 affected and 38 control pugs. Nineteen associated loci (harboring 67 genes in total, including 34 potential candidate genes) and three candidate regions under selection (with four genes within or next to the signal) were identified. The multiple candidate genes identified have implicated functions in bone homeostasis, fibrotic scar tissue, inflammatory responses, or the formation, regulation, and differentiation of cartilage, suggesting the potential relevance of these processes to the pathogenesis of PDM.

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.

Bayesian model and selection signature analyses reveal risk factors for canine atopic dermatitis.

Canine atopic dermatitis is an inflammatory skin disease with clinical similarities to human atopic dermatitis. Several dog breeds are at increased risk for developing this disease but previous genetic associations are poorly defined. To identify additional genetic risk factors for canine atopic dermatitis, we here apply a Bayesian mixture model adapted for mapping complex traits and a cross-population extended haplotype test to search for disease-associated loci and selective sweeps in four dog breeds at risk for atopic dermatitis. We define 15 associated loci and eight candidate regions under selection by comparing cases with controls. One associated locus is syntenic to the major genetic risk locus (Filaggrin locus) in human atopic dermatitis. One selection signal in common type Labrador retriever cases positions across the TBC1D1 gene (body weight) and one signal of selection in working type German shepherd controls overlaps the LRP1B gene (brain), near the KYNU gene (psoriasis). In conclusion, we identify candidate genes, including genes belonging to the same biological pathways across multiple loci, with potential relevance to the pathogenesis of canine atopic dermatitis. The results show genetic similarities between dog and human atopic dermatitis, and future across-species genetic comparisons are hereby further motivated.

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.

Molecular mimicry between Anoctamin 2 and Epstein-Barr virus nuclear antigen 1 associates with multiple sclerosis risk.

Multiple sclerosis (MS) is a chronic inflammatory, likely autoimmune disease of the central nervous system with a combination of genetic and environmental risk factors, among which Epstein-Barr virus (EBV) infection is a strong suspect. We have previously identified increased autoantibody levels toward the chloride-channel protein Anoctamin 2 (ANO2) in MS. Here, IgG antibody reactivity toward ANO2 and EBV nuclear antigen 1 (EBNA1) was measured using bead-based multiplex serology in plasma samples from 8,746 MS cases and 7,228 controls. We detected increased anti-ANO2 antibody levels in MS (P = 3.5 × 10-36) with 14.6% of cases and 7.8% of controls being ANO2 seropositive (odds ratio [OR] = 1.6; 95% confidence intervals [95%CI]: 1.5 to 1.8). The MS risk increase in ANO2-seropositive individuals was dramatic when also exposed to 3 known risk factors for MS: HLA-DRB1*15:01 carriage, absence of HLA-A*02:01, and high anti-EBNA1 antibody levels (OR = 24.9; 95%CI: 17.9 to 34.8). Reciprocal blocking experiments with ANO2 and EBNA1 peptides demonstrated antibody cross-reactivity, mapping to ANO2 [aa 140 to 149] and EBNA1 [aa 431 to 440]. HLA gene region was associated with anti-ANO2 antibody levels and HLA-DRB1*04:01 haplotype was negatively associated with ANO2 seropositivity (OR = 0.6; 95%CI: 0.5 to 0.7). Anti-ANO2 antibody levels were not increased in patients from 3 other inflammatory disease cohorts. The HLA influence and the fact that specific IgG production usually needs T cell help provides indirect evidence for a T cell ANO2 autoreactivity in MS. We propose a hypothesis where immune reactivity toward EBNA1 through molecular mimicry with ANO2 contributes to the etiopathogenesis of MS.

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(-7)) 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(-5)), 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 (pdifference = 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.

cgmisc: enhanced genome-wide association analyses and visualization.

UNLABELLED: 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.

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⁻¹4). 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.

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(-8)) 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.

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.

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.

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.